Importance of genetic architecture in marker selection decisions for genomic prediction

KEY MESSAGE

We demonstrate potential for improved multi-environment genomic prediction accuracy using structural variant markers. However, the degree of observed improvement is highly dependent on the genetic architecture of the trait. Breeders commonly use genetic markers to predict the performance of untested individuals as a way to improve the efficiency of breeding programs. These genomic prediction models have almost exclusively used single nucleotide polymorphisms (SNPs) as their source of genetic information, even though other types of markers exist, such as structural variants (SVs). Given that SVs are associated with environmental adaptation and not all of them are in linkage disequilibrium to SNPs, SVs have the potential to bring additional information to multi-environment prediction models that are not captured by SNPs alone. Here, we evaluated different marker types (SNPs and/or SVs) on prediction accuracy across a range of genetic architectures for simulated traits across multiple environments. Our results show that SVs can improve prediction accuracy, but it is highly dependent on the genetic architecture of the trait and the relative gain in accuracy is minimal. When SVs are the only causative variant type, 70% of the time SV predictors outperform SNP predictors. However, the improvement in accuracy in these instances is only 1.5% on average. Further simulations with predictors in varying degrees of LD with causative variants of different types (e.g., SNPs, SVs, SNPs and SVs) showed that prediction accuracy increased as linkage disequilibrium between causative variants and predictors increased regardless of the marker type. This study demonstrates that knowing the genetic architecture of a trait in deciding what markers to use in large-scale genomic prediction modeling in a breeding program is more important than what types of markers to use.

Linking genetic and environmental factors through marker effect networks to understand trait plasticity

Understanding how plants adapt to specific environmental changes and identifying genetic markers associated with phenotypic plasticity can help breeders develop plant varieties adapted to a rapidly changing climate. Here, we propose the use of marker effect networks as a novel method to identify markers associated with environmental adaptability. These marker effect networks are built by adapting commonly used software for building gene coexpression networks with marker effects across growth environments as the input data into the networks. To demonstrate the utility of these networks, we built networks from the marker effects of ∼2,000 nonredundant markers from 400 maize hybrids across 9 environments. We demonstrate that networks can be generated using this approach, and that the markers that are covarying are rarely in linkage disequilibrium, thus representing higher biological relevance. Multiple covarying marker modules associated with different weather factors throughout the growing season were identified within the marker effect networks. Finally, a factorial test of analysis parameters demonstrated that marker effect networks are relatively robust to these options, with high overlap in modules associated with the same weather factors across analysis parameters. This novel application of network analysis provides unique insights into phenotypic plasticity and specific environmental factors that modulate the genome.

Evaluation of Corn Germplasm Accessions for Resistance to Clavibacter nebraskensis, Causal Agent of Goss's Bacterial Wilt and Leaf Blight

Goss's bacterial wilt and leaf blight of corn (Zea mays), caused by Clavibacter nebraskensis, is a reemerging disease in the Midwestern United States. From 2011 to 2013, field studies and a greenhouse study were conducted to assess the University of Illinois maize inbred collection for putative sources of resistance to Goss's bacterial wilt and leaf blight. This inbred collection consisted of over 2,000 diverse inbred corn lines that have been collected from all over the world. An initial field screen of over 1,000 inbred lines from the collection was conducted in Urbana, IL in 2011. These lines were inoculated with a C. nebraskensis cell suspension and rated for Goss's bacterial wilt and leaf blight severity using a 1-to-9 scale, with a score of 1 being most resistant. Means for Goss's bacterial wilt and leaf blight ratings ranged from 1 to 8.5. The initial screen identified over 150 lines that had high levels of resistance (severity score of ≤2.5). In total, 177 lines were used in the second stage of field screening. In the second stage, average Goss's bacterial wilt and leaf blight severity ranged from 1.1 to 7.4. Nine lines with high levels of resistance in 2011 and 2012 were advanced to the third stage of field screening. The mean Goss's bacterial wilt and leaf blight severity rating of the resistant lines in the last stage was 1.9, while the susceptible check had a mean score of 6.4. These nine lines were also used in the greenhouse to assess whether resistance varied based on inoculating roots, stems, or leaves. Disease severity was significantly (P ≤ 0.05) less when roots were inoculated compared with both leaf and stem inoculations, which were not significantly different from each other. Lines having high levels of field resistance were also found to be resistant in greenhouse screening regardless of inoculation method. Clustering of pedigree distance of the 34 resistant lines (severity score of ≤2.5) with known pedigree information found that 21 clustered with the Lancaster heterotic family, 4 were related to the Iowa Stiff Stalk Synthetic family, and 9 did not cluster with an identifiable heterotic family. These results show that the Lancaster family is an excellent source of Goss's wilt resistance, and that fewer sources of resistance were found in other families. The most resistant lines identified from this research are potential sources of resistance to Goss's bacterial wilt and leaf blight, and their lineage can be used in corn breeding programs to develop resistant hybrids.

Multiple Maize Reference Genomes Impact the Identification of Variants by Genome-Wide Association Study in a Diverse Inbred Panel

Use of a single reference genome for genome-wide association studies (GWAS) limits the gene space represented to that of a single accession. This limitation can complicate identification and characterization of genes located within presence-absence variations (PAVs). In this study, we present the draft de novo genome assembly of 'PHJ89', an 'Oh43'-type inbred line of maize ( L.). From three separate reference genome assemblies ('B73', 'PH207', and PHJ89) that represent the predominant germplasm groups of maize, we generated three separate whole-seedling gene expression profiles and single nucleotide polymorphism (SNP) matrices from a panel of 942 diverse inbred lines. We identified 34,447 (B73), 39,672 (PH207), and 37,436 (PHJ89) transcripts that are not present in the respective reference genome assemblies. Genome-wide association studies were conducted in the 942 inbred panel with both the SNP and expression data values to map (SCMV) resistance. Highlighting the impact of alternative reference genomes in gene discovery, the GWAS results for SCMV resistance with expression values as a surrogate measure of PAV resulted in robust detection of the physical location of a known resistance gene when the B73 reference that contains the gene was used, but not the PH207 reference. This study provides the valuable resource of the Oh43-type PHJ89 genome assembly as well as SNP and expression data for 942 individuals generated from three different reference genomes.

Multiple Maize Reference Genomes Impact the Identification of Variants by Genome-Wide Association Study in a Diverse Inbred Panel

Use of a single reference genome for genome-wide association studies (GWAS) limits the gene space represented to that of a single accession. This limitation can complicate identification and characterization of genes located within presence-absence variations (PAVs). In this study, we present the draft de novo genome assembly of 'PHJ89', an 'Oh43'-type inbred line of maize ( L.). From three separate reference genome assemblies ('B73', 'PH207', and PHJ89) that represent the predominant germplasm groups of maize, we generated three separate whole-seedling gene expression profiles and single nucleotide polymorphism (SNP) matrices from a panel of 942 diverse inbred lines. We identified 34,447 (B73), 39,672 (PH207), and 37,436 (PHJ89) transcripts that are not present in the respective reference genome assemblies. Genome-wide association studies were conducted in the 942 inbred panel with both the SNP and expression data values to map (SCMV) resistance. Highlighting the impact of alternative reference genomes in gene discovery, the GWAS results for SCMV resistance with expression values as a surrogate measure of PAV resulted in robust detection of the physical location of a known resistance gene when the B73 reference that contains the gene was used, but not the PH207 reference. This study provides the valuable resource of the Oh43-type PHJ89 genome assembly as well as SNP and expression data for 942 individuals generated from three different reference genomes.

Wild emmer genome architecture and diversity elucidate wheat evolution and domestication

Wheat (Triticum spp.) is one of the founder crops that likely drove the Neolithic transition to sedentary agrarian societies in the Fertile Crescent more than 10,000 years ago. Identifying genetic modifications underlying wheat's domestication requires knowledge about the genome of its allo-tetraploid progenitor, wild emmer (T. turgidum ssp. dicoccoides). We report a 10.1-gigabase assembly of the 14 chromosomes of wild tetraploid wheat, as well as analyses of gene content, genome architecture, and genetic diversity. With this fully assembled polyploid wheat genome, we identified the causal mutations in Brittle Rachis 1 (TtBtr1) genes controlling shattering, a key domestication trait. A study of genomic diversity among wild and domesticated accessions revealed genomic regions bearing the signature of selection under domestication. This reference assembly will serve as a resource for accelerating the genome-assisted improvement of modern wheat varieties.

Foliar Glyphosate Treatment Alters Transcript and Hormone Profiles in Crown Buds of Leafy Spurge and Induces Dwarfed and Bushy Phenotypes throughout its Perennial Lifecycle

Leafy spurge ( L.) is an invasive weed of North America and its perennial nature attributed to underground adventitious buds (UABs) that undergo seasonal cycles of para-, endo-, and ecodormancy. Recommended rates of glyphosate (∼1 kg ha) destroy aboveground shoots but plants still regenerate vegetatively; therefore, it is considered glyphosate-tolerant. However, foliar application of glyphosate at higher rates (2.2-6.7 kg ha) causes sublethal effects that induce UABs to produce stunted, bushy phenotypes. We investigated the effects of glyphosate treatment (±2.24 kg ha) on vegetative growth, phytohormone, and transcript profiles in UABs under controlled environments during one simulated seasonal cycle. Because shoots derived from UABs of foliar glyphosate-treated plants produced stunted, bushy phenotypes, we could not directly determine if these UABs transitioned through seasonally induced endo- and ecodormancy. However, transcript abundance for leafy spurge dormancy marker genes and principal component analyses suggested that UABs of foliar glyphosate-treated plants transitioned through endo- and ecodormancy. Glyphosate treatment increased shikimate abundance in UABs 7 d after treatment; however, the abundance of shikimate gradually decreased as UABs transitioned through endo- and ecodormancy. The dissipation of shikimate over time suggests that glyphosate's target site was no longer affected, but these changes did not reverse the altered phenotypes observed from UABs of foliar glyphosate-treated leafy spurge. Transcript profiles further indicated that foliar glyphosate treatment significantly affected phytohormone biosynthesis and signaling, particularly auxin transport; gibberellic acid, abscisic acid and jasmonic acid biosynthesis; ethylene responses; and detoxification and cell cycle processes in UABs. These results correlated well with the available phytohormone profiles and altered phenotypes.

Draft Assembly of Elite Inbred Line PH207 Provides Insights into Genomic and Transcriptome Diversity in Maize

Intense artificial selection over the last 100 years has produced elite maize (Zea mays) inbred lines that combine to produce high-yielding hybrids. To further our understanding of how genome and transcriptome variation contribute to the production of high-yielding hybrids, we generated a draft genome assembly of the inbred line PH207 to complement and compare with the existing B73 reference sequence. B73 is a founder of the Stiff Stalk germplasm pool, while PH207 is a founder of Iodent germplasm, both of which have contributed substantially to the production of temperate commercial maize and are combined to make heterotic hybrids. Comparison of these two assemblies revealed over 2500 genes present in only one of the two genotypes and 136 gene families that have undergone extensive expansion or contraction. Transcriptome profiling revealed extensive expression variation, with as many as 10,564 differentially expressed transcripts and 7128 transcripts expressed in only one of the two genotypes in a single tissue. Genotype-specific genes were more likely to have tissue/condition-specific expression and lower transcript abundance. The availability of a high-quality genome assembly for the elite maize inbred PH207 expands our knowledge of the breadth of natural genome and transcriptome variation in elite maize inbred lines across heterotic pools.

Glyphosate's impact on vegetative growth in leafy spurge identifies molecular processes and hormone cross-talk associated with increased branching

BACKGROUND

Leafy spurge (Euphorbia esula) is a perennial weed that is considered glyphosate tolerant, which is partially attributed to escape through establishment of new vegetative shoots from an abundance of underground adventitious buds. Leafy spurge plants treated with sub-lethal concentrations of foliar-applied glyphosate produce new vegetative shoots with reduced main stem elongation and increased branching. Processes associated with the glyphosate-induced phenotype were determined by RNAseq using aerial shoots derived from crown buds of glyphosate-treated and -untreated plants. Comparison between transcript abundance and accumulation of shikimate or phytohormones (abscisic acid, auxin, cytokinins, and gibberellins) from these same samples was also done to reveal correlations.

RESULTS

Transcriptome assembly and analyses confirmed differential abundance among 12,918 transcripts (FDR ≤ 0.05) and highlighted numerous processes associated with shoot apical meristem maintenance and stem growth, which is consistent with the increased number of actively growing meristems in response to glyphosate. Foliar applied glyphosate increased shikimate abundance in crown buds prior to decapitation of aboveground shoots, which induces growth from these buds, indicating that 5-enolpyruvylshikimate 3-phosphate (EPSPS) the target site of glyphosate was inhibited. However, abundance of shikimate was similar in a subsequent generation of aerial shoots derived from crown buds of treated and untreated plants, suggesting EPSPS is no longer inhibited or abundance of shikimate initially observed in crown buds dissipated over time. Overall, auxins, gibberellins (precursors and catabolites of bioactive gibberellins), and cytokinins (precursors and bioactive cytokinins) were more abundant in the aboveground shoots derived from glyphosate-treated plants.

CONCLUSION

Based on the overall data, we propose that the glyphosate-induced phenotype resulted from complex interactions involving shoot apical meristem maintenance, hormone biosynthesis and signaling (auxin, cytokinins, gibberellins, and strigolactones), cellular transport, and detoxification mechanisms.

Glyphosate's impact on vegetative growth in leafy spurge identifies molecular processes and hormone cross-talk associated with increased branching

BACKGROUND

Leafy spurge (Euphorbia esula) is a perennial weed that is considered glyphosate tolerant, which is partially attributed to escape through establishment of new vegetative shoots from an abundance of underground adventitious buds. Leafy spurge plants treated with sub-lethal concentrations of foliar-applied glyphosate produce new vegetative shoots with reduced main stem elongation and increased branching. Processes associated with the glyphosate-induced phenotype were determined by RNAseq using aerial shoots derived from crown buds of glyphosate-treated and -untreated plants. Comparison between transcript abundance and accumulation of shikimate or phytohormones (abscisic acid, auxin, cytokinins, and gibberellins) from these same samples was also done to reveal correlations.

RESULTS

Transcriptome assembly and analyses confirmed differential abundance among 12,918 transcripts (FDR ≤ 0.05) and highlighted numerous processes associated with shoot apical meristem maintenance and stem growth, which is consistent with the increased number of actively growing meristems in response to glyphosate. Foliar applied glyphosate increased shikimate abundance in crown buds prior to decapitation of aboveground shoots, which induces growth from these buds, indicating that 5-enolpyruvylshikimate 3-phosphate (EPSPS) the target site of glyphosate was inhibited. However, abundance of shikimate was similar in a subsequent generation of aerial shoots derived from crown buds of treated and untreated plants, suggesting EPSPS is no longer inhibited or abundance of shikimate initially observed in crown buds dissipated over time. Overall, auxins, gibberellins (precursors and catabolites of bioactive gibberellins), and cytokinins (precursors and bioactive cytokinins) were more abundant in the aboveground shoots derived from glyphosate-treated plants.

CONCLUSION

Based on the overall data, we propose that the glyphosate-induced phenotype resulted from complex interactions involving shoot apical meristem maintenance, hormone biosynthesis and signaling (auxin, cytokinins, gibberellins, and strigolactones), cellular transport, and detoxification mechanisms.

Halomonas sulfidaeris-dominated microbial community inhabits a 1.8 km-deep subsurface Cambrian Sandstone reservoir

A low-diversity microbial community, dominated by the γ-proteobacterium Halomonas sulfidaeris, was detected in samples of warm saline formation porewater collected from the Cambrian Mt. Simon Sandstone in the Illinois Basin of the North American Midcontinent (1.8 km/5872 ft burial depth, 50°C, pH 8, 181 bars pressure). These highly porous and permeable quartz arenite sandstones are directly analogous to reservoirs around the world targeted for large-scale hydrocarbon extraction, as well as subsurface gas and carbon storage. A new downhole low-contamination subsurface sampling probe was used to collect in situ formation water samples for microbial environmental metagenomic analyses. Multiple lines of evidence suggest that this H. sulfidaeris-dominated subsurface microbial community is indigenous and not derived from drilling mud microbial contamination. Data to support this includes V1-V3 pyrosequencing of formation water and drilling mud, as well as comparison with previously published microbial analyses of drilling muds in other sites. Metabolic pathway reconstruction, constrained by the geology, geochemistry and present-day environmental conditions of the Mt. Simon Sandstone, implies that H. sulfidaeris-dominated subsurface microbial community may utilize iron and nitrogen metabolisms and extensively recycle indigenous nutrients and substrates. The presence of aromatic compound metabolic pathways suggests this microbial community can readily adapt to and survive subsurface hydrocarbon migration.

Development and annotation of perennial Triticeae ESTs and SSR markers

Triticeae contains hundreds of species of both annual and perennial types. Although substantial genomic tools are available for annual Triticeae cereals such as wheat and barley, the perennial Triticeae lack sufficient genomic resources for genetic mapping or diversity research. To increase the amount of sequence information available in the perennial Triticeae, three expressed sequence tag (EST) libraries were developed and annotated for Pseudoroegneria spicata , a mixture of both Elymus wawawaiensis and E. lanceolatus , and a Leymus cinereus  × L. triticoides interspecific hybrid. The ESTs were combined into unigene sets of 8 780 unigenes for P. spicata, 11 281 unigenes for Leymus, and 7 212 unigenes for Elymus. Unigenes were annotated based on putative orthology to genes from rice, wheat, barley, other Poaceae, Arabidopsis, and the non-redundant database of the NCBI. Simple sequence repeat (SSR) markers were developed, tested for amplification and polymorphism, and aligned to the rice genome. Leymus EST markers homologous to rice chromosome 2 genes were syntenous on Leymus homeologous groups 6a and 6b (previously 1b), demonstrating promise for in silico comparative mapping. All ESTs and SSR markers are available on an EST information management and annotation database ( http://titan.biotec.uiuc.edu/triticeae/ ).

Characterization of an 18,166 EST dataset for cassava (Manihot esculenta Crantz) enriched for drought-responsive genes

Cassava (Manihot esculenta Crantz) is a staple food for over 600 million people in the tropics and subtropics and is increasingly used as an industrial crop for starch production. Cassava has a high growth rate under optimal conditions but also performs well in drought-prone areas and on marginal soils. To increase the tools for understanding and manipulating drought tolerance in cassava, we generated expressed sequence tags (ESTs) from normalized cDNA libraries prepared from dehydration-stressed and control well-watered tissues. Analysis of a total of 18,166 ESTs resulted in the identification of 8,577 unique gene clusters (5,383 singletons and 3,194 clusters). Functional categories could be assigned to 63% of the unigenes, while another approximately 11% were homologous to hypothetical genes with unclear functions. The remaining approximately 26% were not significantly homologous to sequences in public databases suggesting that some may be novel and putatively specific to cassava. The dehydration-stressed library uncovered numerous ESTs with recognized roles in drought-responses, including those that encode late-embryogenesis-abundant proteins thought to confer osmoprotective functions during water stress, transcription factors, heat-shock proteins as well as proteins involved in signal transduction and oxidative stress. The unigene clusters were screened for short tandem repeats for further development as microsatellite markers. A total of 592 clusters contained 646 repeats, representing 3.3% of the ESTs queried. The ESTs presented here are the first dehydration stress transcriptome of cassava and can be utilized for the development of microarrays and gene-derived molecular markers to further dissect the molecular basis of drought tolerance in cassava.