Prediction ability of genome-wide markers in Pinus taeda L. within and between population is affected by relatedness to the training population and trait genetic architecture.

Genomic prediction has the potential to significantly increase the rate of genetic gain in tree breeding programs. In this study, a clonally replicated population (n = 2063) was used to train a genomic prediction model. The model was validated both within the training population and in a separate population (n = 451). The prediction abilities from random (20% vs 80%) cross validation within the training population were 0.56 and 0.78 for height and stem form, respectively. Removal of all full-sib relatives within the training population resulted in ∼50% reduction in their genomic prediction ability for both traits. The average prediction ability for all 451 individual trees was 0.29 for height and 0.57 for stem form. The degree of genetic linkage (full-sib family, half sib family, unrelated) between the training and validation sets had a strong impact on prediction ability for stem form but not for height. A dominant dwarfing allele, the first to be reported in a conifer species, was discovered via genome-wide association studies on linkage Group 5 that conferred a 0.33-m mean height reduction. However, the QTL was family specific. The rapid decay of linkage disequilibrium, large genome size, and inconsistencies in marker-QTL linkage phase suggest that large, diverse training populations are needed for genomic selection in Pinus taeda L.

Toward genomic selection in Pinus taeda: Integrating resources to support array design in a complex conifer genome.

PREMISE: An informatics approach was used for the construction of an Axiom genotyping array from heterogeneous, high-throughput sequence data to assess the complex genome of loblolly pine (Pinus taeda). METHODS: High-throughput sequence data, sourced from exome capture and whole genome reduced-representation approaches from 2698 trees across five sequence populations, were analyzed with the improved genome assembly and annotation for the loblolly pine. A variant detection, filtering, and probe design pipeline was developed to detect true variants across and within populations. From 8.27 million variants, a total of 642,275 were evaluated and 423,695 of those were screened across a range-wide population. RESULTS: The final informatics and screening approach delivered an Axiom array representing 46,439 high-confidence variants to the forest tree breeding and genetics community. Based on the annotated reference genome, 34% were located in or directly upstream or downstream of genic regions. DISCUSSION: The Pita50K array represents a genome-wide resource developed from sequence data for an economically important conifer, loblolly pine. It uniquely integrates independent projects that assessed trees sampled across the native range. The challenges associated with the large and repetitive genome are addressed in the development of this resource.

Major QTL confer race-nonspecific resistance in the co-evolved Cronartium quercuum f. sp. fusiforme-Pinus taeda pathosystem.

Fusiform rust disease, caused by the endemic fungus Cronartium quercuum f. sp. fusiforme, is the most damaging disease affecting economically important pine species in the southeast United States. Unlike the major epidemics of agricultural crops, the co-evolved pine-rust pathosystem is characterized by steady-state dynamics and high levels of genetic diversity within environments. This poses a unique challenge and opportunity for the deployment of large-effect resistance genes. We used trait dissection to study the genetic architecture of disease resistance in two P. taeda parents that showed high resistance across multiple environments. Two mapping populations (full-sib families), each with ~1000 progeny, were challenged with a complex inoculum consisting of 150 pathogen isolates. High-density linkage mapping revealed three major-effect QTL distributed on two linkage groups. All three QTL were validated using a population of 2057 cloned pine genotypes in a 6-year-old multi-environmental field trial. As a complement to the QTL mapping approach, bulked segregant RNAseq analysis revealed a small number of candidate nucleotide binding leucine-rich repeat genes harboring SNP associated with disease resistance. The results of this study show that in P. taeda, a small number of major QTL can provide effective resistance against genetically diverse mixtures of an endemic pathogen. These QTL vary in their impact on disease liability and exhibit additivity in combination.

Characterizing the oligogenic architecture of plant growth phenotypes informs genomic selection approaches in a common wheat population.

BACKGROUND: Genetic variation in growth over the course of the season is a major source of grain yield variation in wheat, and for this reason variants controlling heading date and plant height are among the best-characterized in wheat genetics. While the major variants for these traits have been cloned, the importance of these variants in contributing to genetic variation for plant growth over time is not fully understood. Here we develop a biparental population segregating for major variants for both plant height and flowering time to characterize the genetic architecture of the traits and identify additional novel QTL. RESULTS: We find that additive genetic variation for both traits is almost entirely associated with major and moderate-effect QTL, including four novel heading date QTL and four novel plant height QTL. FT2 and Vrn-A3 are proposed as candidate genes underlying QTL on chromosomes 3A and 7A, while Rht8 is mapped to chromosome 2D. These mapped QTL also underlie genetic variation in a longitudinal analysis of plant growth over time. The oligogenic architecture of these traits is further demonstrated by the superior trait prediction accuracy of QTL-based prediction models compared to polygenic genomic selection models. CONCLUSIONS: In a population constructed from two modern wheat cultivars adapted to the southeast U.S., almost all additive genetic variation in plant growth traits is associated with known major variants or novel moderate-effect QTL. Major transgressive segregation was observed in this population despite the similar plant height and heading date characters of the parental lines. This segregation is being driven primarily by a small number of mapped QTL, instead of by many small-effect, undetected QTL. As most breeding populations in the southeast U.S. segregate for known QTL for these traits, genetic variation in plant height and heading date in these populations likely emerges from similar combinations of major and moderate effect QTL. We can make more accurate and cost-effective prediction models by targeted genotyping of key SNPs.

Sequence-based mapping identifies a candidate transcription repressor underlying awn suppression at the B1 locus in wheat.

Awns are stiff, hair-like structures which grow from the lemmas of wheat (Triticum aestivum) and other grasses that contribute to photosynthesis and play a role in seed dispersal. Variation in awn length in domesticated wheat is controlled primarily by three major genes, most commonly the dominant awn suppressor Tipped1 (B1). This study identifies a transcription repressor responsible for awn inhibition at the B1 locus. Association mapping was combined with analysis in biparental populations to delimit B1 to a distal region of 5AL colocalized with QTL for number of spikelets per spike, kernel weight, kernel length, and test weight. Fine-mapping located B1 to a region containing only two predicted genes, including C2H2 zinc finger transcriptional repressor TraesCS5A02G542800 upregulated in developing spikes of awnless individuals. Deletions encompassing this candidate gene were present in awned mutants of an awnless wheat. Sequence polymorphisms in the B1 coding region were not observed in diverse wheat germplasm whereas a nearby polymorphism was highly predictive of awn suppression. Transcriptional repression by B1 is the major determinant of awn suppression in global wheat germplasm. It is associated with increased number of spikelets per spike and decreased kernel size.