Association mapping in forest trees and fruit crops

Ultrahigh-density linkage map for cultivated cucumber (Cucumis sativus L.) using a single-nucleotide polymorphism genotyping array

Genotyping arrays are tools for high-throughput genotyping, which is beneficial in constructing saturated genetic maps and therefore high-resolution mapping of complex traits. Since the report of the first cucumber genome draft, genetic maps have been constructed mainly based on simple-sequence repeats (SSRs) or on combinations of SSRs and sequence-related amplified polymorphism (SRAP). In this study, we developed the first cucumber genotyping array consisting of 32,864 single-nucleotide polymorphisms (SNPs). These markers cover the cucumber genome with a median interval of ~2 Kb and have expected genotype calls in parents/F₁ hybridizations as a training set. The training set was validated with Fluidigm technology and showed 96% concordance with the genotype calls in the parents/F₁ hybridizations. Application of the genotyping array was illustrated by constructing a 598.7 cM genetic map based on a ‘9930’ × ‘Gy14’ recombinant inbred line (RIL) population comprised of 11,156 SNPs. Marker collinearity between the genetic map and reference genomes of the two parents was estimated at R² = 0.97. We also used the array-derived genetic map to investigate chromosomal rearrangements, regional recombination rate, and specific regions with segregation distortions. Finally, 82% of the linkage-map bins were polymorphic in other cucumber variants, suggesting that the array can be applied for genotyping in other lines. The genotyping array presented here, together with the genotype calls of the parents/F₁ hybridizations as a training set, should be a powerful tool in future studies with high-throughput cucumber genotyping. An ultrahigh-density linkage map constructed by this genotyping array on RIL population may be invaluable for assembly improvement, and for mapping important cucumber QTLs.

Developing single nucleotide polymorphism (SNP) markers from transcriptome sequences for identification of longan (Dimocarpus longan) germplasm

Longan (Dimocarpus longan Lour.) is an important tropical fruit tree crop. Accurate varietal identification is essential for germplasm management and breeding. Using longan transcriptome sequences from public databases, we developed single nucleotide polymorphism (SNP) markers; validated 60 SNPs in 50 longan germplasm accessions, including cultivated varieties and wild germplasm; and designated 25 SNP markers that unambiguously identified all tested longan varieties with high statistical rigor (P<0.0001). Multiple trees from the same clone were verified and off-type trees were identified. Diversity analysis revealed genetic relationships among analyzed accessions. Cultivated varieties differed significantly from wild populations (F st=0.300; P<0.001), demonstrating untapped genetic diversity for germplasm conservation and utilization. Within cultivated varieties, apparent differences between varieties from China and those from Thailand and Hawaii indicated geographic patterns of genetic differentiation. These SNP markers provide a powerful tool to manage longan genetic resources and breeding, with accurate and efficient genotype identification.

Multiple QTLs Linked to Agro-Morphological and Physiological Traits Related to Drought Tolerance in Potato

Dissection of the genetic architecture of adaptation and abiotic stress-related traits is highly desirable for developing drought-tolerant potatoes and enhancing the resilience of existing cultivars, particularly as agricultural production in rain-fed areas may be reduced by up to 50 % by 2020. The "DMDD" potato progeny was developed at International Potato Center (CIP) by crossing the sequenced double monoploid line DM and a diploid cultivar of the Solanum tuberosum diploid Andigenum Goniocalyx group. Recently, a high-density integrated genetic map based on single nucleotide polymorphism (SNP), diversity array technology (DArT), simple sequence repeats (SSRs), and amplified fragment length polymorphism (AFLP) markers was also made available for this population. Two trials were conducted, in greenhouse and field, for drought tolerance with two treatments each, well-watered and terminal drought, in which watering was suspended 60 days after planting. The DMDD population was evaluated for agro-morphological and physiological traits before and after initiation of stress, at multiple time points. Two dense parental genetic maps were constructed using published genotypic data, and quantitative trait locus (QTL) analysis identified 45 genomic regions associated with nine traits in well-watered and terminal drought treatments and 26 potentially associated with drought stress. In this study, the strong influence of environmental factors besides water shortage on the expression of traits and QTLs reflects the multigenic control of traits related to drought tolerance. This is the first study to our knowledge in potato identifying QTLs for drought-related traits in field and greenhouse trials, giving new insights into genetic architecture of drought-related traits. Many of the QTLs identified have the potential to be used in potato breeding programs for enhanced drought tolerance.

Association mapping for kernel phytosterol content in almond

Almond kernels are a rich source of phytosterols, which are important compounds for human nutrition. The genetic control of phytosterol content has not yet been documented in almond. Association mapping (AM), also known as linkage disequilibrium (LD), was applied to an almond germplasm collection in order to provide new insight into the genetic control of total and individual sterol contents in kernels. Population structure analysis grouped the accessions into two principal groups, the Mediterranean and the non-Mediterranean. There was a strong subpopulation structure with LD decaying with increasing genetic distance, resulting in lower levels of LD between more distant markers. A significant impact of population structure on LD in the almond cultivar groups was observed. The mean r(2) -value for all intra-chromosomal loci pairs was 0.040, whereas, the r(2) for the inter-chromosomal loci pairs was 0.036. For analysis of association between the markers and phenotypic traits five models were tested. The mixed linear model (MLM) approach using co-ancestry values from population structure and kinship estimates (K model) as covariates identified a maximum of 13 significant associations. Most of the associations found appeared to map within the interval where many candidate genes involved in the sterol biosynthesis pathway are predicted in the peach genome. These findings provide a valuable foundation for quality gene identification and molecular marker assisted breeding in almond.

Western white pine SNP discovery and high-throughput genotyping for breeding and conservation applications

BACKGROUND

Western white pine (WWP, Pinus monticola Douglas ex D. Don) is of high interest in forest breeding and conservation because of its high susceptibility to the invasive disease white pine blister rust (WPBR, caused by the fungus Cronartium ribicola J. C. Fisch). However, WWP lacks genomic resource development and is evolutionarily far away from plants with available draft genome sequences. Here we report a single nucleotide polymorphism (SNP) study by bulked segregation-based RNA-Seq analysis.

RESULTS

A collection of resistance germplasm was used for construction of cDNA libraries and SNP genotyping. Approximately 36-89 million 2 × 100-bp reads were obtained per library and de-novo assembly generated the first shoot-tip reference transcriptome containing a total of 54,661 unique transcripts. Bioinformatic SNP detection identified >100,000 high quality SNPs in three expressed candidate gene groups: Pinus highly conserved genes (HCGs), differential expressed genes (DEGs) in plant defense response, and resistance gene analogs (RGAs). To estimate efficiency of in-silico SNP discovery, genotyping assay was developed by using Sequenom iPlex and it unveiled SNP success rates from 40.1% to 61.1%. SNP clustering analyses consistently revealed distinct populations, each composed of multiple full-sib seed families by parentage assignment in the WWP germplasm collection. Linkage disequilibrium (LD) analysis identified six genes in significant association with major gene (Cr2) resistance, including three RGAs (two NBS-LRR genes and one receptor-like protein kinase -RLK gene), two HCGs, and one DEG. At least one SNP locus provided an excellent marker for Cr2 selection across P. monticola populations.

CONCLUSIONS

The WWP shoot tip transcriptome and those validated SNP markers provide novel genomic resources for genetic, evolutionary and ecological studies. SNP loci of those candidate genes associated with resistant phenotypes can be used as positional and functional variation sites for further characterization of WWP major gene resistance against C. ribicola. Our results demonstrate that integration of RNA-seq-based transcriptome analysis and high-throughput genotyping is an effective approach for discovery of a large number of nucleotide variations and for identification of functional gene variants associated with adaptive traits in a non-model species.

The complex jujube genome provides insights into fruit tree biology

The jujube (Ziziphus jujuba Mill.), a member of family Rhamnaceae, is a major dry fruit and a traditional herbal medicine for more than one billion people. Here we present a high-quality sequence for the complex jujube genome, the first genome sequence of Rhamnaceae, using an integrated strategy. The final assembly spans 437.65 Mb (98.6% of the estimated) with 321.45 Mb anchored to the 12 pseudo-chromosomes and contains 32,808 genes. The jujube genome has undergone frequent inter-chromosome fusions and segmental duplications, but no recent whole-genome duplication. Further analyses of the jujube-specific genes and transcriptome data from 15 tissues reveal the molecular mechanisms underlying some specific properties of the jujube. Its high vitamin C content can be attributed to a unique high level expression of genes involved in both biosynthesis and regeneration. Our study provides insights into jujube-specific biology and valuable genomic resources for the improvement of Rhamnaceae plants and other fruit trees.

The jujube is a major dry fruit crop in China and is commonly used for medicinal purposes. Here the authors sequence the genome and transcriptome of the most widely cultivated jujube cultivar, Dongzao, and highlight the genetic and molecular basis of agronomically important jujube traits, such as vitamin C content.

Development and validation of a 20K single nucleotide polymorphism (SNP) whole genome genotyping array for apple (Malus × domestica Borkh)

High-density SNP arrays for genome-wide assessment of allelic variation have made high resolution genetic characterization of crop germplasm feasible. A medium density array for apple, the IRSC 8K SNP array, has been successfully developed and used for screens of bi-parental populations. However, the number of robust and well-distributed markers contained on this array was not sufficient to perform genome-wide association analyses in wider germplasm sets, or Pedigree-Based Analysis at high precision, because of rapid decay of linkage disequilibrium. We describe the development of an Illumina Infinium array targeting 20K SNPs. The SNPs were predicted from re-sequencing data derived from the genomes of 13 Malus × domestica apple cultivars and one accession belonging to a crab apple species (M. micromalus). A pipeline for SNP selection was devised that avoided the pitfalls associated with the inclusion of paralogous sequence variants, supported the construction of robust multi-allelic SNP haploblocks and selected up to 11 entries within narrow genomic regions of ±5 kb, termed focal points (FPs). Broad genome coverage was attained by placing FPs at 1 cM intervals on a consensus genetic map, complementing them with FPs to enrich the ends of each of the chromosomes, and by bridging physical intervals greater than 400 Kbps. The selection also included ∼3.7K validated SNPs from the IRSC 8K array. The array has already been used in other studies where ∼15.8K SNP markers were mapped with an average of ∼6.8K SNPs per full-sib family. The newly developed array with its high density of polymorphic validated SNPs is expected to be of great utility for Pedigree-Based Analysis and Genomic Selection. It will also be a valuable tool to help dissect the genetic mechanisms controlling important fruit quality traits, and to aid the identification of marker-trait associations suitable for the application of Marker Assisted Selection in apple breeding programs.

A combined linkage and regional association mapping validation and fine mapping of two major pleiotropic QTLs for seed weight and silique length in rapeseed (Brassica napus L.)

BACKGROUND

Seed weight (SW) and silique length (SL) are important determinants of the yield potential in rapeseed (Brassica napus L.). However, the genetic basis of both traits is poorly understood. The main objectives of this study were to dissect the genetic basis of SW and SL in rapeseed through the preliminary mapping of quantitative trait locus (QTL) by linkage analysis and fine mapping of the target major QTL by regional association analysis.

RESULTS

Preliminary linkage mapping identified thirteen and nine consensus QTLs for SW and SL, respectively. These QTLs explained 0.7-67.1% and 2.1-54.4% of the phenotypic variance for SW and SL, respectively. Of these QTLs, three pairs of SW and SL QTLs were co-localized and integrated into three unique QTLs. In addition, the significance level and genetic effect of the three co-localized QTLs for both SW and SL showed great variation before and after the conditional analysis. Moreover, the allelic effects of the three QTLs for SW were highly consistent with those for SL. Two of the three co-localized QTLs, uq.A09-1 (mean R(2) = 20.1% and 19.0% for SW and SL, respectively) and uq.A09-3 (mean R(2) = 13.5% and 13.2% for SW and SL, respectively), were detected in all four environments and showed the opposite additive-effect direction. These QTLs were validated and fine mapped (their confidence intervals were narrowed down from 5.3 cM to 1 cM for uq.A09-1 and 13.2 cM to 2.5 cM for uq.A09-3) by regional association analysis with a panel of 576 inbred lines, which has a relatively rapid linkage disequilibrium decay (0.3 Mb) in the target QTL region.

CONCLUSIONS

A few QTLs with major effects and several QTLs with moderate effects might contribute to the natural variation of SW and SL in rapeseed. The meta-, conditional and allelic effect analyses suggested that pleiotropy, rather than tight linkage, was the genetic basis of the three pairs of co-localized of SW and SL QTLs. Regional association analysis was an effective and highly efficient strategy for the direct fine mapping of target major QTL identified by preliminary linkage mapping.

Genome-wide distribution of genetic diversity and linkage disequilibrium in a mass-selected population of maritime pine

Background

The accessibility of high-throughput genotyping technologies has contributed greatly to the development of genomic resources in non-model organisms. High-density genotyping arrays have only recently been developed for some economically important species such as conifers. The potential for using genomic technologies in association mapping and breeding depends largely on the genome wide patterns of diversity and linkage disequilibrium in current breeding populations. This study aims to deepen our knowledge regarding these issues in maritime pine, the first species used for reforestation in south western Europe.

Results

Using a new map merging algorithm, we first established a 1,712 cM composite linkage map (comprising 1,838 SNP markers in 12 linkage groups) by bringing together three already available genetic maps. Using rigorous statistical testing based on kernel density estimation and resampling we identified cold and hot spots of recombination. In parallel, 186 unrelated trees of a mass-selected population were genotyped using a 12k-SNP array. A total of 2,600 informative SNPs allowed to describe historical recombination, genetic diversity and genetic structure of this recently domesticated breeding pool that forms the basis of much of the current and future breeding of this species. We observe very low levels of population genetic structure and find no evidence that artificial selection has caused a reduction in genetic diversity. By combining these two pieces of information, we provided the map position of 1,671 SNPs corresponding to 1,192 different loci. This made it possible to analyze the spatial pattern of genetic diversity (H_e) and long distance linkage disequilibrium (LD) along the chromosomes. We found no particular pattern in the empirical variogram of H_e across the 12 linkage groups and, as expected for an outcrossing species with large effective population size, we observed an almost complete lack of long distance LD.

Conclusions

These results are a stepping stone for the development of strategies for studies in population genomics, association mapping and genomic prediction in this economical and ecologically important forest tree species.

Identification of genetic loci associated with fire blight resistance in Malus through combined use of QTL and association mapping

Fire blight, incited by the enterobacterium Erwinia amylovora, is a destructive disease of Rosaceae, particularly of apples and pears. There are reports on the molecular mechanisms underlying E. amylovora pathogenesis and how the host activates its resistance mechanism. The host's resistance mechanism is quantitatively controlled, although some major genes might also be involved. Thus far, quantitative trait loci (QTL) mapping and differential expression studies have been used to elucidate those genes and/or genomic regions underlying quantitative resistance present in the apple genome. In this study, an effort is undertaken to dissect the genetic basis of fire blight resistance in apple using both QTL and genome-wide association mapping. On the basis of an F1 pedigree of 'Coop 16' × 'Coop 17' and a genome-wide association study (GWAS) mapping population of Malus accessions (species, old and new cultivars and selections), new QTLs and associations have been identified. A total of three QTLs for resistance to fire blight, with above 95% significant logarithm of odds threshold value of 2.5, have been identified on linkage groups (LGs) 02, 06, and 15 of the apple genome with phenotypic variation explained values of 14.7, 20.1 and 17.4, respectively. Although elevated P-values with signals for marker-trait associations are observed for some LGs, these are not found to be significant. However, a total of 34 significant associations, with P-values ≥0.02, have been detected including 8 for lesion length at 7 days following inoculation (PL1), 14 for lesion length at 14 days following inoculation (PL2), and 12 for shoot length.

A multi-population consensus genetic map reveals inconsistent marker order among maps likely attributed to structural variations in the apple genome

Genetic maps serve as frameworks for determining the genetic architecture of quantitative traits, assessing structure of a genome, as well as aid in pursuing association mapping and comparative genetic studies. In this study, a dense genetic map was constructed using a high-throughput 1,536 EST-derived SNP GoldenGate genotyping platform and a global consensus map established by combining the new genetic map with four existing reliable genetic maps of apple. The consensus map identified markers with both major and minor conflicts in positioning across all five maps. These major inconsistencies among marker positions were attributed either to structural variations within the apple genome, or among mapping populations, or genotyping technical errors. These also highlighted problems in assembly and anchorage of the reference draft apple genome sequence in regions with known segmental duplications. Markers common across all five apple genetic maps resulted in successful positioning of 2875 markers, consisting of 2033 SNPs and 843 SSRs as well as other specific markers, on the global consensus map. These markers were distributed across all 17 linkage groups, with an average of 169±33 marker per linkage group and with an average distance of 0.70±0.14 cM between markers. The total length of the consensus map was 1991.38 cM with an average length of 117.14±24.43 cM per linkage group. A total of 569 SNPs were mapped onto the genetic map, consisting of 140 recombinant individuals, from our recently developed apple Oligonucleotide pool assays (OPA). The new functional SNPs, along with the dense consensus genetic map, will be useful for high resolution QTL mapping of important traits in apple and for pursuing comparative genetic studies in Rosaceae.