Genetic control of juvenile growth and botanical architecture in an ornamental woody plant, Prunus mumeSieb. et Zucc. as revealed by a high-density linkage map
BMC Geneticsvolume 15, Article number: S1 (2014)
Mei, Prunus mume Sieb. et Zucc., is an ornamental plant popular in East Asia and, as an important member of genus Prunus, has played a pivotal role in systematic studies of the Rosaceae. However, the genetic architecture of botanical traits in this species remains elusive. This paper represents the first genome-wide mapping study of quantitative trait loci (QTLs) that affect stem growth and form, leaf morphology and leaf anatomy in an intraspecific cross derived from two different mei cultivars. Genetic mapping based on a high-density linkage map constricted from 120 SSRs and 1,484 SNPs led to the detection of multiple QTLs for each trait, some of which exert pleiotropic effects on correlative traits. Each QTL explains 3-12% of the phenotypic variance. Several leaf size traits were found to share common QTLs, whereas growth-related traits and plant form traits might be controlled by a different set of QTLs. Our findings provide unique insights into the genetic control of tree growth and architecture in mei and help to develop an efficient breeding program for selecting superior mei cultivars.
Mei, Prunus mume Sieb. et Zucc., a species of genus Prunus, is a popular ornamental plant, originated in Southwest China  and widely cultivated in the entire East Asia [1, 2]. Its cold hardiness by blooming in winter or early spring, plus its many prominent ornamental features, such as colorful corollas, varying flower forms, and pleasant fragrance, have made it a symbol of spirit in Chinese culture, favorably praised by litterateurs and ordinary people . Fruits of mei have also been used as raw material to make Chinese herbal medicine beneficial for human health . As an important member of sub-family Prunoideae, mei is a key step towards constructing a phylogenetic tree for family Rosaceae, thought to play a pivotal role in understanding the evolution of woody plants .
Because of its significant value in biological research and practical cultivation, mei has received increasing attention during the past years [3–8]. Fang et al.  developed a set of molecular markers, such as AFLP and SNP, to investigate the genetic relatedness and diversity of 50 cultivars of fruiting mei from China and Japan. Similar work using AFLP markers was conducted by Yang et al.  to analyze the genetic diversity of ornamental mei and compare it with that of other related species. Li et al.  developed more informative multiallelic microsatellite markers, i.e., simple sequence repeats (SSRs), from 20 mei plants, particularly useful to study the genetic structure of natural populations in mei. Using two mei cultivars, Fenban and Kouzi Yudie, and five segregating progeny randomly chosen from the F1 intraspecific hybrid family of these two cultivars, Sun et al.  performed the genome-wide characterization of SSRs in the mei genome and construct a first genetic linkage map of mei using 144 SSR markers. Despite these progresses in mei genetic studies, almost nothing is known about the genetic control of its botanical traits of ornamental and biological value.
More recently, with the advent and widespread application of next-generation technologies, the genetic studies of mei have entered a new era in which multiscale data collected at the molecular, cell and organ levels enables geneticists to characterize the genetic architecture of complex phenotypes and construct a genotype to phenotype map for mei. Right after genomes of apple and strawberry, both belonging to Rosaceae, have been sequenced [9, 10], Zhang et al.  have for the first time sequenced the mei genome, providing an impetus for studying the comparative genomics of Rosaceae and mapping important genes that contribute to mei traits. Based on the mei reference genome, Sun et al.  were able to identify hundreds of thousands of SNPs for cultivars Fenban and Kouzi Yudie. The F1 family of these two cultivars was genotyped for a couple of thousands of SNPs. By adding these segregating SNPs to the SSR linkage map, previously reported in Sun et al. , a high-density genetic map for mei has been generated.
In this article, we report on the detection of quantitative trait loci (QTLs) that affect stem growth, stem form and leaf morphological traits in the juvenile seedlings of mei using a segregating F1 family derived from cultivars Fenban and Kouzi Yudie . Early growth and its morphological components, such as leaf size, in the first year of growth in the field are important traits associated with the ability of mei to build itself to tolerate and resist to adverse conditions, particularly low temperature and drought in early spring. Our high-density genetic map allows the genome-wide mapping and identification of QTLs responsible for the early performance of mei in the field. Results from QTL mapping are not only useful for marker-assisted selection and breeding of rigorous growth traits in mei, but also help to explore the commonality of genetic control for growth traits in Rosaceae through comparing with QTL discoveries in other species such as apple and strawberry.
Unlike an inbred line homozygous for all loci, an outcrossing line is complex in genetic composition, in which some loci are homozygous whereas the others are heterozygous . Thus, the F1 cross of two outcrossing parents, like Fenban and Kouzi Yudie, may generate a segregating progeny as long as one parent is heterozygous for some loci. There are four possible types of segregating markers for an outcrossing family [12, 13]: (1) multi (3 or 4)-allelic intercross markers, (2) bi-allelic intercross markers, (3) testcross markers that are heterozygous for one parent but homozygous for the other, and (4) testcross markers that are opposite to (3). These marker types produce four, three, two and two distinguishable genotypes in the progeny, respectively. Statistical models implemented with different numbers of effect parameters are used to identify QTLs from these markers .
The juvenile phenotypic traits studied in the mapping population of mei are classified into three categories: (1) stem growth and form, described by stem height, stem diameter, and stem slenderness (measured as diameter/height ratio), (2) leaf morphology, including leaf length, length width, single leaf area, petiole length, and leaf shape (measured as length width/length ratio), and leaf structure, including the number of veins per leaf and the number of veins per unit area of a leaf (vein density). Figure 1 illustrates the histograms of each trait in the mapping population, showing an approximate normal distribution and pronounced variability. Each trait was associated with individual markers by a maximum likelihood approach. Plots of log-likelihood ratios for each trait over linkage groups are given in Figures 2, 3, 4, in which the genomic distribution of significant QTLs is shown. We did not detect many significant QTLs for stem growth, only with two for height growth jointly accounting for 7% of the phenotypic variation and three for diameter growth explaining 16% of the phenotypic variation together (Table 1; Figure 2). One diameter QTL on linkage group 8 is an intercross type, acting in an overdominant manner (d/a = 15).
Growth component traits, like leaf length, leaf width and leaf area, were observed to involve larger genetic components explained by QTLs (Table 1; Figure 3). Four QTLs contribute jointly to 20% of the phenotypic variation for leaf length, whereas over a half of the phenotypic variation for leaf width is explained by five QTLs. For leaf area, two QTLs detected account for 15% of its phenotypic variation. One QTL associated with marker PMSNP00307 on linkage group 5 pleiotropically affect both leaf length and width. Two multiallelic intercross QTLs also on linkage group 5 are pleiotropic QTLs for leaf width and area. Although these two traits are controlled by intercross QTLs, the dominant effects are relatively small, compared with their additive effects. A total of three QTLs explain about 9% of the phenotypic variation for leaf petiole.
Relative to growth traits, we found more QTLs involved in form traits; for example, five for stem taper and six for leaf shape (Table 1; Figure 2 and 3). A total of 15% and 23% of the phenotypic variation are due to these QTLs for the two traits, respectively. The number of veins per leaf is controlled by multiple QTLs from different linkage groups 1, 4, 5, and 6, totally explaining 20% of its phenotypic variation (Table 1; Figure 4). Five QTLs, all on linkage group 5, were detected to affect the density of veins, with 17% of the phenotypic variation explained.
Genetic mapping has proven to be a powerful tool in studying the genetic architecture and complex traits and designing marker-assisted selection programs for many species. Genetic linkage maps are generally constructed using a segregating progeny, such as the backcross, F2, or recombinant inbred lines, derived from two homozygous inbred lines. For perennial trees, however, it is difficult or impossible to obtain such inbred lines owing to their long-generation interval, high heterozygosity, and high inbreeding depression . On the other hand, because of their high heterozygosity, the F1 full-sib family produced by crossing two trees may provide an adequate amount of information for linkage analysis . In such a family, there are many types of segregating markers. Earlier pseudo-test backcross designs by Grattapalia and Sederoff [16, 17] can make use of markers that are heterozygous in one parent but homozygous in the other, taking advantage of linkage analysis models developed for the backcross population. Since a more sophisticated model for linkage mapping has been developed [11–13, 18], any type of markers segregating in a full-sib family can be analyzed by simultaneously estimating the linkage and linkage phases. This method was successfully used in poplar tree , sugarcane [19, 20], a yellow passion fruit population , rubber tree  and peach .
This is a first study for mapping QTLs that control botanical traits in mei. By crossing two mei cultivars, a full-sib family was generated as a mapping population. Zhang et al.'s  sequencing result provides sufficient information for genotyping the mei genome. The high-density linkage map constructed by SSR markers, SNP markers and InDels  allows mei QTLs to be identified and estimated. In this mapping population, dramatic differences at phenotypic and genetic levels were observed in growth-related traits and botanical form traits in mei. Although the linkage map used covers a large portion of the mei genome, we did not identify many QTLs for stem growth traits. This may be due to two reasons. First, the mei trees are in their early stage of establishment in the field, when environmental perturbations are a major factor to affect tree growth. As trees are established, genes play an increasing role in growth and growth component traits. Such a transition pattern of genetic control after the establishment was observed in an experimental plantation of poplar hybrids [24, 25]. Second, growth is a complex trait that is likely to involve a complex network of genetic interactions. We expect that epistasis due to different genes which main effects are not significant may contribute to the phenotypic variation of stem growth traits. A more powerful model that can analyze and estimate the genetic effects of all markers at the same time is crucial to confirm this speculation .
Although growth and its components, such as ones related to leaf size, have been mapped in many woody plants, QTL mapping of several important botanical traits, like stem taper, leaf shape and leaf anatomy, has received little attention. To our best knowledge, this is the first study aimed to map QTLs that control the number and density of veins. As physiological pipelines that transport water, nutrients and hormones, leaf veins have been thought to be associated with plant growth and adaptability . The vein QTLs identified from this study will help to understand the genetic variation of leaf venation. Wu et al.  presented one of the first studies that map QTLs for leaf shape in forest trees, and found different patterns of action of QTLs on leaf size and leaf shape. In Wu's  study, QTLs for stem form were found in juvenile poplar trees. Given its ornamental value, botanical form traits in mei are part of breeding objectives. This study has for the first time reported on the detection of QTLs that control stem shape and leaf shape, providing useful information for marker-assisted selection of good-shaped mei cultivars. It is noted that different genomic regions control growth and shape, suggesting different genetic machineries that generate the phenotypic variation of these two types of traits.
We detected the pleiotropic control of the same QTL over two allometrically related traits, leaf length and leaf width. Similar pleiotropic QTLs were also observed for leaf width and leaf area. All these findings are of significant value to unveil the genetic basis of morphological and developmental integration as a mechanism for plants to adapt to environmental changes. Our data was about juvenile treesm, in which there is limitation to make a strong inference about developmental mechanisms. Yet, the current result from young trees, plus those from subsequent years, will enable us to link genes and development into a platform of interplay at which we are in a better position to chart the genotype-phenotype map through developmental trajectories .
Materials and methods
F1hybrids and DNA extraction
Two mei cultivars, P. mume 'Fenban' (BJFU1210120013) and P. mume 'Kouzi Yudie' (BJFU1210120022), were selected from the Qingdao Mei Garden, Qingdao, China (36°04′N, 120°20′E), differing in many growth and morphological traits. The cross of the two cultivars generated a segregating F1 population, of which 190 seedlings (Voucher specimen accession number: BJFU1210120025-0214) were grown in the Xiao Tangshan Horticultural Trial, Beijing, China (40°02′N, 115°50′E). Total DNA was extracted from fresh young leaves of each seedling according to the plant genomic DNA extraction Kit (TIANGEN, Beijing, China) following the manufacturer's instructions.
Marker genotyping and map construction
Sun et al.  described a procedure of identifying and genotyping SSR markers for the F1 hybrids of mei, including SSR primer design and screening and PCR amplification. A total of 144 multiallelic intercross markers were genotyped for the F1 hybrid population. The description of the procedure to identify SNP markers and InDels for mei was shown in Sun et al. . To the end, 105 multiallelic intercross markers, 395 biallelic intercross markers and 1004 testcross markers segregating in the hybrids were generated.
Sun et al.  used JoinMap version 4  to construct a genetic linkage map from SSR markers. The estimated recombination fractions between markers were converted to map distance in centiMorgan using the Kosambi map function. The map is composed of eight linkage groups paralleling to the haploid chromosome number of the mei genome. The total length of the map is 670 cM, with an average marker distance of 5 cM. The positions of SNPs, InDels and SSRs were identified as CDS, intron, 5′UTR, 3′UTR and intergenic regions according to mei genome GFF files. Thus, relative positions of SNPs and InDels on the SSR linkage map can be determined.
During the fast-growing season of mei, usually in July or August, we measured leaf size and morphology for each F1 seedling. Three representative leaves chosen for phenotyping from the same tree are those located at the main stem with leaf plastochron index of 10 to 12. For each chosen leaf, leafblade length and width were measured, from which leaf areas were calculated. The number of stomata was counted for each leaf. At the end of the first growing season in the field, each seedling was evaluated for its main stem height and stem base diameter. Growth and its component traits used for QTL mapping are the height (HT) and diameter (DIA) of the main stem, leaf length (LL), leaf width (LW) and leaf area (LA). The botanical form traits of mei were derived from measured traits, including stem shape (measured by the ratio of DIA to HT) and leaf shape (measured by the ratio of LW to LL). Also, the density of veins on the leaf was calculated as the ratio of veins to leaf area. The averages of trait values over three measured leaves were used for QTL analysis.
Since a high-density map was used for genetic mapping, we directly associated marker genotypes with phenotypic traits to detect significant QTLs using a likelihood approach. In this particular full-sib family, there are multiple marker types, testcross, biallelic intercross and multiallelic intercross. Here, we describe the model to analyze the genetic effects of a multiallelic intercross QTL [31, 32]. Assume two alleles P1 and P2 for parent P and two alleles Q1 and Q2 for parent Q, which generate four progeny genotypes, along with genotypic values (μ11, μμ12, μ21, μ22), expressed as
where μ is the overall mean, a1 is the allelic (additive) effect contributed by parent P, a2 is the allelic effect contributed by parent Q, and d is the dominant effect due to the interaction between alleles from the two different parents. The quantitative genetic analysis of testcross QTLs and biallelic intercross QTLs have been available in previous studies [17, 28, 29]. Biallelic intercross QTLs generate three genotypes in the progeny, allowing one additive effect (a) and one dominant effect (d) to be estimated. For testcross QTLs with two progeny genotypes, only one additive effect (a) can be estimated. In each case, the proportion of the total phenotypic variance explained by each QTL was calculated.
The significance of QTLs was tested by calculating the log-ratio of likelihoods under the null hypothesis (there is no QTL) and alternative hypothesis (there is a QTL) and comparing it with the chi-square distribution. When multiple SNPs were included for QTL mapping, the significance of SNP needs to be adjusted using the false positive rate (FDR) approach. The genome-wide threshold of significance was determined after the FDR adjustment.
Publication of this work is supported by grants from the Ministry of Science and Technology (2011AA100207, 2013AA102607), the State Forestry Administration of China (201004012), the Fundamental Research Funds for the Central Universities (NO.BLX2013011), and "One-thousand Person Plan" Award.
Chen JY: Chinese Mei Flowers (in Chinese). 1996, Hainan Publishing House, Haikou, China, 14-58.
Chu MY: China Fruit Records - Mei(in Chinese). 1999, China Forestry Press, Beijing, 12-60.
Zhang Q, Chen W, Sun L, Zhao F, Huang B, et al: The genome of Prunus mume. Nat Commun. 2012, 3: 1318-
Fang J, Twito T, Zhang Z, Chao CT: Genetic relationships among fruiting-mei (Prunus mume Sieb. et Zucc.) cultivars evaluated with AFLP and SNP markers. Genome. 2006, 49: 1256-1264. 10.1139/g06-097.
Yang C-D, Zhang J-W, Yan X-L, Bao M-Z: Genetic relatedness and genetic diversity of ornamental mei (Prunus mume Sieb. et Zucc.) as analysed by AFLP markers. Tree Genet Genom. 2008, 4: 255-262. 10.1007/s11295-007-0106-0.
Li X, Shangguan L, Song C, Wang C, Gao Z, et al: Analysis of expressed sequence tags from Prunus mume flower and fruit and development of simple sequence repeat markers. BMC Genet. 2010, 11: 66-
Sun L, Yang W, Zhang Q, Cheng T, Pan H, Xu Z, Zhang J, Chen C: Genome-wide characterization and linkage mapping of simple sequence repeats in mei (Prunus mume Sieb. et Zucc.). PloS ONE. 2013, 8 (3): e59562-10.1371/journal.pone.0059562.
Sun L, Zhang Q, Xu Z, Yang W, Guo Y, Lu J, Pan H, Chneg T, Cai M: Genome-wide DNA polymorphisms in two cultivars of mei (Prunus mume sieb. et zucc.). BMC Genet. 2013a, 14: 98-
Velasco R, Zharkikh A, Affourtit J, Dhingra A, Cestaro A, et al: The genome of the domesticated apple (Malus × domestica Borkh.). Nat Genet. 2010, 42: 833-839. 10.1038/ng.654.
Shulaev V, Sargent DJ, Crowhurst RN, Mockler TC, Folkerts O, et al: The genome of woodland strawberry (Fragaria vesca). Nat Genet. 2011, 43: 109-116. 10.1038/ng.740.
Maliepaard C, Jansen J, van Ooijen JW: Linkage analysis in a full-sib family of an outbreeding plant species: Overview and consequences for applications. Genet Res. 1997, 70: 237-250. 10.1017/S0016672397003005.
Wu R, Ma C, Painter I, Zeng Z: Simultaneous maximum likelihood estimation of linkage and linkage phases in outcrossing species. Theor Pop Biol. 2002, 61: 349-363. 10.1006/tpbi.2002.1577.
Wu R, Ma C-X, Wu SS, Zeng Z-B: Linkage mapping of sex-specific differences. Genet Res. 2002, 79: 85-96. 10.1017/S0016672301005389.
Lin M, Lou X, Chang M, Wu R: A general statistical framework for mapping quantitative trait loci in nonmodel systems: Issue for characterizing linkage phases. Genetics. 2003, 165: 901-913.
Wu RL, Zeng Z-B, McKend SE, O'Malley DM: The case for molecular mapping in forest tree breeding. Plant Breed Rev. 2000, 19: 41-68.
Grattapaglia D, Sederoff R: Genetic linkage maps of Eucalyptus grandis and Eucalyptus urophylla using a pseudo-testcross: Mapping strategy and RAPD markers. Genetics. 1994, 137: 1121-1137.
Grattapaglia D, Bertolucci FL, Sederoff RR: Genetic mapping of QTLs controlling vegetative propagation in Eucalyptus grandis and E. urophylla using a pseudo-testcross mapping strategy and RAPD markers. Theor Appl Genet. 1995, 90: 933-947.
Ooijen JV: JoinMap®4, software for the calculation of genetic linkage maps in experimental populations. 2006, Wageningen, The Netherlands, 1-55.
Garcia AAF, Kido EA, Meza AN, Souza HMB, Pinto LR, et al: Development of an integrated genetic map of a sugarcane (Saccharum spp.) commercial cross, based on a maximum-likelihood approach for estimation of linkage and linkage phases. Theor Appl Genet. 2006, 112: 298-314. 10.1007/s00122-005-0129-6.
Oliveira KM, Pinto LR, Marconi TG, Margarido GRA, Pastina MM, et al: Functional integrated genetic linkage map based on EST-markers for a sugarcane (Saccharum spp.) commercial cross. Mol Breed. 2007, 20: 189-208. 10.1007/s11032-007-9082-1.
Oliveira EJ, Vieira MLC, Garcia AAF, Munhoz CF, Margarido GRA, et al: An integrated molecular map of yellow passion fruit based on simultaneous maximum-likehood estimation of linkage and linkage phases. J Am Soc Hort Sci. 2008, 133: 35-41.
Souza LM, Gazaffi R, Mantello CC, Silva CC, Garcia D, et al: QTL mapping of growth-related traits in a full-sib family of rubber tree (Hevea brasiliensis) evaluated in a sub-tropical climate. PLoS ONE. 2013, 8 (4): e61238-10.1371/journal.pone.0061238.
Romeu JF, Monforte AJ, Sánchez G, Granell A, García-Brunton J, Badenes ML, Ríos G: Quantitative trait loci affecting reproductive phenology in peach. BMC Plant Biol. 2014, 14 (1): p1-10.1186/1471-2229-14-1.
Wu RL, Stettler RF: Quantitative genetics of growth, development in Populus. I. A three-generation comparison of tree architecture during the first two years of growth. Theor Appl Genet. 1994, 89: 1046-1054.
Wu R, Stettler RF: The genetic dissection of juvenile canopy structure and function in a three-generation pedigree of Populus. Trees Struc Fun. 1996, 11: 99-108.
Li JH, Das K, Fu G, Li R, Wu RL: The Bayesian lasso for genome-wide association studies. Bioinformatics. 2011, 27 (4): 516-523. 10.1093/bioinformatics/btq688.
Sack L, Scoffoni C: Leaf venation: structure, function, development, evolution, ecology and applications in the past, present and future. New Phytol. 2013, 198: 983-1000. 10.1111/nph.12253.
Wu RL, Bradshaw HD, Stettler RF: Molecular genetics of growth and development in Populus. V. Mapping quantitative trait loci affecting leaf variation. Am J Bot. 1997, 84: 143-153. 10.2307/2446076.
Wu RL: Genetic mapping of QTLs affecting tree growth and architecture in Populus: Implications for ideotype breeding. Theor Appl Genet. 1998, 96: 447-457. 10.1007/s001220050761.
Wu RL, Lin M: Functional mapping - how to map and study the genetic architecture of dynamic complex traits. Nat Rev Genet. 2006, 7: 229-237. 10.1038/nrg1804.
Tong C, Wang Z, Zhang B, Shi J, Wu R: 3FunMap: full-sib family functional mapping of dynamic traits. Bioinformatics. 2011, 27: 2006-2008. 10.1093/bioinformatics/btr306.
Tong C, Zhang B, Wang Z, Xu M, Pang XM, Si JN, Huang MR, Wu RL: Multiallelic epistatic model for an outbred cross and mapping algorithm of interactive quantitative trait loci. BMC Plant Biol. 2011, 11: 148-10.1186/1471-2229-11-148.
This article has been published as part of BMC Genetics Volume 15 Supplement 1, 2014: Selected articles from the International Symposium on Quantitative Genetics and Genomics of Woody Plants. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenet/supplements/15/S1.
The authors declared that they have no competing interests.
Conceived and designed the experiments: QZ. Performed the experiments: LS. Analyzed the data: YW CZ XL. Contributed reagents/materials/analysis tools: LS XY TC KM WY HP JW QZ. Wrote the paper: LS RW.
Lidan Sun, Yaqun Wang contributed equally to this work.