Genome-wide discovery of DNA polymorphism in Brassica rapa

Single nucleotide polymorphisms (SNPs) and/or insertion/deletions (InDels) are frequent sequence variations in the plant genome, which can be developed as molecular markers for genetic studies on crop improvement. The ongoing Brassica rapa genome sequencing project has generated vast amounts of sequence data useful in genetic research. Here, we report a genome-wide survey of DNA polymorphisms in the B. rapa genome based on the 557 bacterial artificial clone sequences of B. rapa ssp. pekinensis cv. Chiifu. We identified and characterized 21,311 SNPs and 6,753 InDels in the gene space of the B. rapa genome by re-sequencing 1,398 sequence-tagged sites (STSs) in eight genotypes. Comparison of our findings with a B. rapa genetic linkage map confirmed that STS loci were distributed randomly over the B. rapa whole genome. In the 1.4 Mb of aligned sequences, mean nucleotide polymorphism and diversity were θ = 0.00890 and π = 0.00917, respectively. Additionally, the nucleotide diversity in introns was almost three times greater than that in exons, and the frequency of observed InDel was almost 17 times higher in introns than in exons. Information regarding SNPs/InDels obtained here will provide an important resource for genetic studies and breeding programs of B. rapa.     

An Improved Brassica rapa Genetic Linkage Map and Locus-specific Variations in a Doubled Haploid Population

In this study, we describe the construction of an improved Chinese cabbage genetic linkage map by integrating simple sequence repeats (SSRs) and insertion/deletion polymorphisms (InDels) into a previously published map of a doubled haploid (DH) population. The population was derived from a cross between the Chinese cabbage line BY (Brassica rapa ssp. pekinensis) and a European turnip line MM (Brassica rapa L. ssp. rapifera). A total of 629 markers were aligned to ten linkage groups, with a total map length of 1,173.8 cM, and an average distance between markers of 1.87 cm. Of the 126 SSRs and 133 InDels mapped, 46 and 34 were novel, respectively. A comparison of the linkage map with the B. rapa genome showed that more than 93 % of the markers, including 112 SSRs and 129 InDels, could be anchored unambiguously to a specific location on one of the ten chromosomes. In most cases, the order of markers on the linkage map and physical map was similar; however, the majority of linkage groups contained a number of markers whose positions were either transposed or had moved slightly forwards or backwards. During microspore culture, it was observed that 11 SSRs and one InDel showed either variation in size, or the appearance of new marker bands in the DH lines. As a first step to addressing this SSR/InDel marker instability, six SSR and one InDel loci were sequenced, which revealed that the size variation was due mainly to changes in repeat-motif number or to the insertion/deletion of new fragments of DNA.     

Development of Species-Specific InDel Markers in Citrus

Citrus taxonomy is complex owing to the existence of a wide range of species: Poncirus is used mainly for rootstock; Fortunella produces small fruit and edible pericarp; and Citrus comprises the most widespread fruit crop species worldwide. Rapidly increasing genome resources from different citrus species facilitate the development of convenient and genome-wide molecular markers that can be applied to both inter- and intra-species analyses. In this study, by comparing the genome sequences of four citrus species, a set of 1958 InDels were identified and 453 candidate InDels were converted into PCR-based markers. Among these candidate InDels, 268 (65%) exhibited length polymorphisms from 30 bp to 200 bp when applied to seven species from the genera Poncirus, Fortunella and Citrus. Seven InDel markers exhibited high intraspecific polymorphisms in a natural pummelo population. The results showed that the InDel markers are effective for both inter- and intra-specific variation and identification analyses. These InDel markers are expected to be applied to germplasm identification, phylogenetic analysis, genetic diversity evaluation and marker-assisted breeding in citrus.     

Fine genetic and physical mapping of the CRb gene conferring resistance to clubroot disease in Brassica rapa

The use of clubroot resistance (CR) genes is an effective and economical approach for controlling Plasmodiophora brassicae, the causal agent of clubroot disease in Chinese cabbage (Brassica rapa) and other Brassica crops. In a previous study, we identified and mapped the CRb locus on chromosome A03 of B. rapa in the doubled-haploid (DH) line ‘CR Shinki DH line’ of Chinese cabbage. In this study, CRb, a dominant gene conferring resistance to pathotype 4 of P. brassicae, was finely mapped in combination with bulked segregant analysis and bioinformatics analysis (BIA). Using 1,486 highly susceptible individuals and 2,896 individuals from two separate F₂ populations of ‘702-5’ (B. rapa ssp. chinensis) ×  ‘CR Shinki DH line,’ the CRb locus was narrowed to a region of approximately 0.14 cM between two flanking markers, TCR79 and TCR108. The sequences of seven newly developed markers linked to CRb were landed on bacterial artificial chromosome (BAC) of the reference B. rapa ‘Chiifu-401-42’ by BIA, and a physical map consisting of three BAC clones was constructed. The CRb locus was defined as an interval of approximately 83.5 kb on a BAC clone (KBrB085J21). The target interval contained one Toll-interleukin-1 receptor/nucleotide-binding site/leucine-rich repeat (TIR–NBS–LRR) gene, one NBS–LRR gene, and several putative regulatory genes in the B. rapa genome. The CRb gene was tightly linked to two other CR genes, CRa and CRb ^Kato . These results provide useful information for isolation of the CRb gene and tightly linked molecular markers for breeding CR in B. rapa.     

Deep re-sequencing of a widely used maintainer line of hybrid rice for discovery of DNA polymorphisms and evaluation of genetic diversity

Genetic diversity within parental lines of hybrid rice is the foundation of heterosis utilization and yield improvement. Previous studies have suggested that genetic diversity was narrow in cytoplasmic male sterile (CMS/A line) and restorer lines (R line) for Three-line hybrid rice. However, the genetic diversity within maintainer lines (B line), especially at a genome-wide scale, remains largely unknown. In the present study, we performed deep re-sequencing of the elite maintainer line V20B (Oryza sativa L. ssp. indica). We then compared the V20B sequence with the 93-11 (Oryza sativa L. ssp. indica) genome sequence. 112.1 × 106 paired-end reads (PE reads) were generated with approximately 30-fold sequencing depth. The V20B PE reads uniquely covered 87.6 % of the 93-11 genome sequence. Overall, a total of 660,778 single-nucleotide polymorphism (SNPs) and 266,301 insertions and deletions (InDels) were identified, yielding an average of 2.1 SNPs/kb and 0.8 InDels/kb. Genome-wide distribution of the SNPs and InDels was non-random, and variation-rich and variation-poor regions were identified in all chromosomes. A total of 20,562 non-synonymous SNPs spanning 8,854 genes were annotated. Our results identified DNA polymorphisms at the genome-wide scale and uncovered the high level of genetic diversity between V20B and 93-11. Our results proved that next-generation sequencing technologies can be powerful tools to study genome-wide DNA polymorphisms, to query genetic diversity, and to enable molecular improvement efforts with Three-line hybrid rice. Further, our results also indicated that 93-11 could be used as core germplasm for the improvement of wild-abortive CMS lines and the maintainer lines.     

In Silico Identification and Experimental Validation of Insertion–Deletion Polymorphisms in Tomato Genome

Comparative analysis of the genome sequences of Solanum lycopersicum variety Heinz 1706 and S. pimpinellifolium accession LA 1589 using MUGSY software identified 145 695 insertion–deletion (InDel) polymorphisms. A selected set of 3029 candidate InDels (≥2 bp) across the entire tomato genome were subjected to PCR validation, and 82.4% could be verified. Of 2272 polymorphic InDels between LA 1589 and Heinz 1706, 61.6, 45.2, and 31.6% were polymorphic in 8 accessions of S. pimpinellifolium, 4 accessions of S. lycopersicum var. cerasiforme, and 10 varieties of S. lycopersicum, respectively. Genetic distance was 0.216 in S. pimpinellifolium, 0.202 in S. lycopersicum var. cerasiforme, and 0.108 in S. lycopersicum. The data suggested a reduction of genetic variation from S. pimpinellifolium to S. lycopersicum var. cerasiforme and S. lycopersicum. Cluster analysis showed that the 8 accessions of S. pimpinellifolium were in one group, whereas 4 accessions of S. lycopersicum var. cerasiforme and 10 varieties of S. lycopersicum were in the same group.     

The effects of artificial selection on sugar metabolism and transporter genes in grape

Sugar content is a key feature of grape quality. The sugar content of grapes has been significantly improved after nearly a thousand years of artificial selection. However, the mechanism underlying the changes in the grape sugar content during the process of artificial selection remains largely unknown although several genes involved in sugar metabolism and transportation in grape have been identified. In this study, the genomes of 13 wild Vitis species and 14 cultivated Vitis vinifera accessions were resequenced to 2–5 X depth using the Illumina Hiseq2000 platform. Genetic variation of 138 genes involved in sugar biosynthesis and transport was investigated, and 7,690 and 12,717 single nucleotide polymorphisms/insertions and deletions (SNPs/InDel) were identified within the cultivated V. vinifera and wild Vitis species, respectively. The percentages of SNPs/InDels were 0.93 and 1.54 % in cultivated and wild species, respectively, and the wild Vitis species had 1.65-fold more SNPs/InDels than the cultivated V. vinifera. Moreover, the distribution of SNPs/InDels in gene regions was also investigated. Eight genes (HT4, PPFTK4, PPFTK6, PMT3, SPS1, HT8, HT15, SUSy1) showed low level of allelic diversity in cultivated species, suggesting they might have undergone purifying selection during the domestication process of grapes. Our genome DNA resequencing data provided a valuable resource for analyzing the effects of artificial selection on trait-related pathways in grape. The result that eight genes showed lower level of DNA variation in cultivated species than in wild species will be very helpful in understanding sugar accumulation in grapes.     

Development of genome-wide DNA polymorphism database for map-based cloning of rice genes

DNA polymorphism is the basis to develop molecular markers that are widely used in genetic mapping today. A genome-wide rice (Oryza sativa) DNA polymorphism database has been constructed in this work using the genomes of Nipponbare, a cultivar of japonica, and 93-11, a cultivar of indica. This database contains 1,703,176 single nucleotide polymorphisms (SNPs) and 479,406 Insertion/Deletions (InDels), approximately one SNP every 268 bp and one InDel every 953 bp in rice genome. Both SNPs and InDels in the database were experimentally validated. Of 109 randomly selected SNPs, 107 SNPs (98.2%) are accurate. PCR analysis indicated that 90% (97 of 108) of InDels in the database could be used as molecular markers, and 68% to 89% of the 97 InDel markers have polymorphisms between other indica cultivars (Guang-lu-ai 4 and Long-te-pu B) and japonica cultivars (Zhong-hua 11 and 9522). This suggests that this database can be used not only for Nipponbare and 93-11, but also for other japonica and indica cultivars. While validating InDel polymorphisms in the database, a set of InDel markers with each chromosome 3 to 5 marker was developed. These markers are inexpensive and easy to use, and can be used for any combination of japonica and indica cultivars used in this work. This rice DNA polymorphism database will be a valuable resource and important tool for map-based cloning of rice gene, as well as in other various research on rice (http://shenghuan.shnu.edu.cn/ricemarker).     

Cross-transferability and Polymorphic Potential of Genomic STMS Markers of Brassica Species

The present study was carried out with the objective of evaluating genomic STMS markers developed earlier in Brassica napus, B. oleracea, B. rapa and B. nigra for their use in Brassica juncea and B. carinata. Ninety-six of the 100 STMS markers used under standardized annealing temperatures and gel concentrations produced clear reproducible amplification pattern. For majority of the markers 60 °C annealing temperature and 3.5% metaphor agarose gel were found suitable. High cross-transferability of STMS markers to related Brassica species including B. carinata (91.6%) and B. juncea (87.5%) suggested the possibility of utilizing these markers for genome analysis in the species where no such markers are available. The ‘B’ genome derived markers showed lower level of transferability to the ‘A’ and ‘C’ genome Brassica species. The potential of STMS markers to detect polymorphism among Brassica species and genera was 98.9%. The level of inter-specific polymorphism was much higher than the intea-specific polymorphism. The markers capable of revealing polymorphism among Brassica species and genera would be useful in Brassica introgression breeding programme. The polymorphic markers were found efficient in establishing the expected evolutionary relationships among the six different Brassica species and two related genera. Low level of intra-specific polymorphism revealed by these markers suggested use of a large set of such markers for various applications in Brassica genetics, genomics and breeding.     

Sequence variation in 3′UTR region of crtRB1 gene and its effect on β-carotene accumulation in maize kernel

β-carotene fortification of maize has emerged as a potential, long-term and sustainable approach to alleviate vitamin A deficiency in humans. Among the several genes involved in the carotenoid biosynthetic pathway, the 543 bp allele at crtRB1 3′TE (Transposable Element) gene (allele 1, without insertion) is associated with higher β-carotene accumulation. Estimation of β-carotene through high performance liquid chromatography showed that the CIMMYT genotypes with allele 1 had high kernel β-carotene content whereas the Indian inbreds with the same allele had low β-carotene content. To know the reason for this variation, allele 1 of crtRB1 3′TE gene was sequenced from a set of 11 diverse maize inbreds collected from CIMMYT and Indian germplasm. The sequence data of the allele 1 revealed the presence of 13 single nucleotide polymorphisms (SNPs) and 7 insertions and deletions (InDels). Exonic region had two SNPs, intronic region had one SNP and one InDel, whereas 3′-untranslated region (UTR) region of the gene showed 10 SNPs and 6 InDels. Among the several SNPs and InDels, SNP4, SNP13, InDel6 and InDel7 identified in the 3′-UTR region clearly differentiated the high and the low β-carotene genotypes. These 3′-UTR polymorphisms in allele 1 of the crtRB1 3′TE gene could be associated with the variation in kernel β-carotene accumulation by regulating the translation and stability of the mRNA. The SNPs and the InDels associated with higher level of β-carotene will be used as a gene-based marker(s) in selection of genotypes and to develop biofortified maize hybrids to alleviate vitamin A deficiency in humans.     

Construction of a genetic map based on high-throughput SNP genotyping and genetic mapping of a TuMV resistance locus in Brassica rapa

Brassica rapa is a member of the Brassicaceae family and includes vegetables and oil crops that are cultivated worldwide. The introduction of durable resistance against turnip mosaic virus (TuMV) into agronomically important cultivars has been a significant challenge for genetic and horticultural breeding studies of B. rapa. Based on our previous genome-wide analysis of DNA polymorphisms between the TuMV-resistant doubled haploid (DH) line VC40 and the TuMV-susceptible DH line SR5, we constructed a core genetic map of the VCS-13M DH population, which is composed of 83 individuals derived from microspore cultures of a F₁ cross between VC40 and SR5, by analyzing the segregation of 314 sequence-characterized genetic markers. The genetic markers correspond to 221 SNPs and 31 InDels of genes as well as 62 SSRs, covering 1,115.9 cM with an average distance of 3.6 cM between the adjacent marker loci. The alignment and orientation of the constructed map showed good agreement with the draft genome sequence of Chiifu, thus providing an efficient strategy to map genic sequences. Using the genetic map, a novel dominant TuMV resistance locus (TuMV-R) in the VCS-13M DH population was identified as a 0.34 Mb region in the short arm of chromosome A6 in which four CC–NBS–LRR resistance genes and two pathogenesis-related-1 genes reside. The genetic map developed in this study can play an important role in the genetic study of TuMV resistance and the molecular breeding of B. rapa.     

Evaluation of Microsatellite Markers for Genome Mapping in Indian Mustard (Brassica juncea L)

Microsatellites are robust markers for genome mapping, gene tagging and marker assisted selection. The genus Brassica, having a large and complex genome, requires such type of markers for various applications in genetics and breeding. A set of 202 microsatellite markers were used to screen two parental genotypes of Indian mustard (Brassica juncea) namely, ‘Varuna’, an indigenous cultivar and BEC144, an exotic collection from Poland, of which 36 (17.8%) were informative and usable for segregation analysis. The polymorphic markers detected heterozygosity in advanced generation recombinant inbred lines (RILs) developed earlier from the cross Varuna × BEC144 with a varying frequency that ranged from 0% to 23.5%. Normal Mendelian segregation for majority of microsatellite markers was observed. Eleven markers showed significant deviation from the expected 1:1 segregation ratio. Twelve markers were assigned to six different linkage groups of Indian mustard genome map. The level of polymorphism between the parents and the percentage of useful informative markers as observed in this study, suggested that many more markers are needed to achieve a reasonable coverage of mustard genome. This is the first report on the evaluation of microsatellite markers for genome mapping in B. juncea.     

Development of cost-effective single nucleotide polymorphism marker assays for genetic diversity analysis in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Competitive allele-specific PCR (KASPar) assay is a user-friendly system that provides flexibility in the numbers of single nucleotide polymorphisms (SNPs) and genotypes. Based on Illumina-GA-IIx genomic data from 10 genotypes with a broad genetic background, 3183 SNPs were selected for KASPar assays development, and 568 were finally converted and selected for Brassica rapa germplasm characterization (17.8%) on the basis of reproducibility, missing data rate, and uniform genetic distribution. High levels of polymorphism of these markers across 231 B. rapa genotypes were verified, illustrating by high polymorphic information content (averaged 0.34), minor allele frequency (0.37), genetic diversity (0.45), and the low observed heterozygosity (0.10). Based on the SNP dataset, structure and principal coordinates analysis, and neighbor-joining phylogenetic methods were used to examine the population structure and gave highly consistent results. The 231 accessions were divided into the four primary subspecies, representing 99 accessions from B. rapa ssp. pekinensis, 85 from B. rapa ssp. chinensis, 30 from B. rapa ssp. rapifera, and 17 from B. rapa ssp. oleifera and were further subdivided into 12 lower-order clusters according to different morphotypes. The genetic variability and pairwise fixation index analysis revealed that the ssp. pekinensis accessions possess the most extensive genetic variation among the four subspecies. The KASPar system is highly useful for validating SNPs and will be valuable for genetics research and breeding applications in B. rapa.     

Genome-wide single nucleotide polymorphisms and insertion–deletions of Oryza sativa L. subsp. japonica cultivars grown near the northern limit of rice cultivation

Single nucleotide polymorphisms (SNPs) and insertions–deletions (InDels) are valuable molecular markers for molecular breeding among genetically closely related cultivars. Rice (Oryza sativa L. subsp. japonica) cultivars grown in Hokkaido (45–42°N), the northernmost region of rice paddy cultivation in Japan, have been bred for over 100 years for adaptation to low summer temperatures together with high yield and good eating quality. In this study, for 10 closely related rice cultivars released in Hokkaido and cultivar Koshihikari, we identified genome-wide SNPs and InDels by next-generation sequencing. More than 29 million reads from the Hokkaido cultivars, each 101 nucleotides long, were uniquely mapped to the Nipponbare reference genome. The average of the total nucleotide length of all uniquely mapped reads corresponded to 10.9 times (3,978 Mb with genome coverage of 90.7 %) the Nipponbare reference genome. An average of 99,955 putative SNPs (1.8 times the number in Koshihikari) and 14,617 putative InDels (also 1.8 times the number in Koshihikari) were detected in Hokkaido cultivars relative to the Nipponbare genome, which enabled analyses of the inheritance of pedigree haplotypes of four cultivars, SNPs and InDels among closely related Hokkaido cultivars, and haplotype blocks unique to Hokkaido cultivars. The comprehensive SNP and InDel data provide DNA marker resources and will facilitate quantitative trait locus analysis of biparental mapping of very closely related Hokkaido cultivars. Furthermore, the haplotype blocks unique to Hokkaido cultivars represent ideal genetic regions for improvement of cultivars to be grown near the northern and southern limits of rice cultivation.     

Genome-wide identification, classification and expression analysis of genes encoding putative fasciclin-like arabinogalactan proteins in Chinese cabbage (Brassica rapa L.)

Fasciclin-like arabinogalactan proteins (FLAs), a subclass of arabinogalactan proteins (AGPs), have both predicted AGP-like glycosylated regions and putative fasciclin (FAS) domains, which may function in cell adhesion and communication. Previous studies have identified 21, 27, and 34 FLAs in Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), and wheat (Triticum aestivum), respectively. In this study, we identified 33 FLAs in the annotated genome of Chinese cabbage (Brassica rapa ssp. pekinensis line Chiifu-401-42). Sequence analysis indicated that FAS domains each contain two highly conserved regions, named H1 and H2, and that 17 FLAs from B. rapa (BrFLAs) possess both of these regions. Prediction of glycosylphosphatidylinositol (GPI) modification sites suggested that 15 BrFLAs were GPI-anchored to the plasma membrane. Additionally, 25 BrFLAs may have been duplicated during the processes that shaped the triplicated genome of the mesopolyploid B. rapa. Expression analyses indicated that BrFLA1, BrFLA11, BrFLA13, BrFLA28 and BrFLA32 were specifically expressed in inflorescence. Meanwhile, BrFLA9 (homologous to AtFLA12) is specifically expressed in stem, and BrFLA6/22 (homologous to AtFLA11) is also highly expressed in stem, suggesting BrFLA6/9/22 may have the same functions as AtFLA11/12 in A. thaliana. Taken together, the identification and bioinformatic analysis of FLAs in B. rapa will open the way for studying their biological functions in plant growth and development as well as evolutionary history of this gene family from A. thaliana to B. rapa.     

Genome-wide insertion–deletion (InDel) marker discovery and genotyping for genomics-assisted breeding applications in chickpea

We developed 21,499 genome-wide insertion–deletion (InDel) markers (2- to 54-bp in silico fragment length polymorphism) by comparing the genomic sequences of four (desi, kabuli and wild C. reticulatum) chickpea [Cicer arietinum (L.)] accessions. InDel markers showing 2- to 6-bp fragment length polymorphism among accessions were abundant (76.8%) in the chickpea genome. The physically mapped 7,643 and 13,856 markers on eight chromosomes and unanchored scaffolds, respectively, were structurally and functionally annotated. The 4,506 coding (23% large-effect frameshift mutations) and regulatory InDel markers were identified from 3,228 genes (representing 11.7% of total 27,571 desi genes), suggesting their functional relevance for trait association/genetic mapping. High amplification (97%) and intra-specific polymorphic (60–83%) potential and wider genetic diversity (15–89%) were detected by genome-wide 6,254 InDel markers among desi, kabuli and wild accessions using even a simpler cost-effective agarose gel-based assay. This signifies added advantages of this user-friendly genetic marker system for manifold large-scale genotyping applications in laboratories with limited infrastructure and resources. Utilizing 6,254 InDel markers-based high-density (inter-marker distance: 0.212 cM) inter-specific genetic linkage map (ICC 4958 × ICC 17160) of chickpea as a reference, three major genomic regions harboring six flowering and maturity time robust QTLs (16.4–27.5% phenotypic variation explained, 8.1–11.5 logarithm of odds) were identified. Integration of genetic and physical maps at these target QTL intervals mapped on three chromosomes delineated five InDel markers-containing candidate genes tightly linked to the QTLs governing flowering and maturity time in chickpea. Taken together, our study demonstrated the practical utility of developing and high-throughput genotyping of such beneficial InDel markers at a genome-wide scale to expedite genomics-assisted breeding applications in chickpea.     

High‐throughput multiplex cpDNA resequencing clarifies the genetic diversity and genetic relationships among Brassica napus,Brassica rapa and Brassica oleracea

Brassica napus (rapeseed) is a recent allotetraploid plant and the second most important oilseed crop worldwide. The origin of B. napus and the genetic relationships with its diploid ancestor species remain largely unresolved. Here, chloroplast DNA (cpDNA) from 488 B. napus accessions of global origin, 139 B. rapa accessions and 49 B. oleracea accessions were populationally resequenced using Illumina Solexa sequencing technologies. The intraspecific cpDNA variants and their allelic frequencies were called genomewide and further validated via EcoTILLING analyses of the rpo region. The cpDNA of the current global B. napus population comprises more than 400 variants (SNPs and short InDels) and maintains one predominant haplotype (Bncp1). Whole‐genome resequencing of the cpDNA of Bncp1 haplotype eliminated its direct inheritance from any accession of the B. rapa or B. oleracea species. The distribution of the polymorphism information content (PIC) values for each variant demonstrated that B. napus has much lower cpDNA diversity than B. rapa; however, a vast majority of the wild and cultivated B. oleracea specimens appeared to share one same distinct cpDNA haplotype, in contrast to its wild C‐genome relatives. This finding suggests that the cpDNA of the three Brassica species is well differentiated. The predominant B. napus cpDNA haplotype may have originated from uninvestigated relatives or from interactions between cpDNA mutations and natural/artificial selection during speciation and evolution. These exhaustive data on variation in cpDNA would provide fundamental data for research on cpDNA and chloroplasts.     

Genome-Wide Discovery of DNA Polymorphisms in Mei (<Emphasis Type="Italic">Prunus mume</Emphasis> Sieb. et Zucc.), an Ornamental Woody Plant,with Contrasting Tree Architecture and their Functional Relevance for Weeping Trait

Next-generation sequencing technologies provide opportunities to ascertain the genetic basis of phenotypic differences, even in the closely related cultivars via detection of large amount of DNA polymorphisms. In this study, we performed whole-genome re-sequencing of two mei cultivars with contrasting tree architecture. 75.87 million 100 bp pair-end reads were generated, with 92 % coverage of the genome. Re-sequencing data of two former upright mei cultivars were applied for detecting DNA polymorphisms, since we were more interested in variations conferring weeping trait. Applying stringent parameters, 157,317 mutual single nucleotide polymorphisms (SNPs) and 15,064 mutual insertions-deletions (InDels) were detected and found unevenly distributed within and among the mei chromosomes, which lead to the discovery of 220 high-density, 463 low-density SNP regions together with 80 high-density InDel regions. Additionally, 322 large-effect SNPs and 433 large-effect InDels were detected, and 10.09 % of the SNPs were observed in coding regions. 5.25 % SNPs in coding regions resulted in non-synonymous changes. Ninety SNPs were chosen randomly for validation using high-resolution melt analysis. 93.3 % of the candidate SNPs contained the predicted SNPs. Pfam analysis was further conducted to better understand SNP effects on gene functions. DNA polymorphisms of two known QTL loci conferring weeping trait and their functional effect were also analyzed thoroughly. This study highlights promising functional markers for molecular breeding and a whole-genome genetic basis of weeping trait in mei.