Integrating Coexpression Networks with GWAS to Prioritize Causal Genes in Maize

Prioritization of causal genes from genome-wide association studies by Bayesian data integration across loci

MOTIVATION

Genome-wide association studies (GWAS) have identified genetic variants, usually single-nucleotide polymorphisms (SNPs), associated with human traits, including disease and disease risk. These variants (or causal variants in linkage disequilibrium with them) usually affect the regulation or function of a nearby gene. A GWAS locus can span many genes, however, and prioritizing which gene or genes in a locus are most likely to be causal remains a challenge. Better prioritization and prediction of causal genes could reveal disease mechanisms and suggest interventions.

RESULTS

We describe a new Bayesian method, termed SigNet for significance networks, that combines information both within and across loci to identify the most likely causal gene at each locus. The SigNet method builds on existing methods that focus on individual loci with evidence from gene distance and expression quantitative trait loci (eQTL) by sharing information across loci using protein-protein and gene regulatory interaction network data. In an application to cardiac electrophysiology with 226 GWAS loci, only 46 (20%) have within-locus evidence from Mendelian genes, protein-coding changes, or colocalization with eQTL signals. At the remaining 180 loci lacking functional information, SigNet selects 56 genes other than the minimum distance gene, equal to 31% of the information-poor loci and 25% of the GWAS loci overall. Assessment by pathway enrichment demonstrates improved performance by SigNet. Review of individual loci shows literature evidence for genes selected by SigNet, including PMP22 as a novel causal gene candidate.

Association of the benzoxazinoid pathway with boron homeostasis in maize

Both deficiency and toxicity of the micronutrient boron lead to severe reductions in crop yield. Despite this agricultural importance, the molecular basis underlying boron homeostasis in plants remains unclear. To identify molecular players involved in boron homeostasis in maize (Zea mays L.), we measured boron levels in the Goodman-Buckler association panel and performed genome-wide association studies. These analyses identified a benzoxazinless (bx) gene, bx3, involved in the biosynthesis of benzoxazinoids, such as 2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one (DIMBOA), which are major defense compounds in maize. Genes involved in DIMBOA biosynthesis are all located in close proximity in the genome, and benzoxazinoid biosynthesis mutants, including bx3, are all DIMBOA deficient. We determined that leaves of the bx3 mutant have a greater boron concentration than those of B73 control plants, which corresponded with enhanced leaf tip necrosis, a phenotype associated with boron toxicity. By contrast, other DIMBOA-deficient maize mutants did not show altered boron levels or the leaf tip necrosis phenotype, suggesting that boron is not associated with DIMBOA. Instead, our analyses suggest that the accumulation of boron is linked to the benzoxazinoid intermediates indolin-2-one (ION) and 3-hydroxy-ION. Therefore, our results connect boron homeostasis to the benzoxazinoid plant defense pathway through bx3 and specific intermediates, rendering the benzoxazinoid biosynthesis pathway a potential target for crop improvement under inadequate boron conditions.

Prioritization of candidate genes for major QTLs governing yield traits employing integrated multi-omics approach in rice (Oryza sativa L.)

Rapidly identifying candidate genes underlying major QTLs is crucial for improving rice (Oryza sativa L.). In this study, we developed a workflow to rapidly prioritize candidate genes underpinning 99 major QTLs governing yield component traits. This workflow integrates multiomics databases, including sequence variation, gene expression, gene ontology, co-expression analysis, and protein-protein interaction. We predicted 206 candidate genes for 99 reported QTLs governing ten economically important yield-contributing traits using this approach. Among these, transcription factors belonging to families of MADS-box, WRKY, helix-loop-helix, TCP, MYB, GRAS, auxin response factor, and nuclear transcription factor Y subunit were promising. Validation of key prioritized candidate genes in contrasting rice genotypes for sequence variation and differential expression identified Leucine-Rich Repeat family protein (LOC_Os03g28270) and cytochrome P450 (LOC_Os02g57290) as candidate genes for the major QTLs GL1 and pl2.1, which govern grain length and panicle length, respectively. In conclusion, this study demonstrates that our workflow can significantly narrow down a large number of annotated genes in a QTL to a very small number of the most probable candidates, achieving approximately a 21-fold reduction. These candidate genes have potential implications for enhancing rice yield.

Identifying hub genes and key functional modules in leaf tissue of Populus species based on WGCNA

As one of the most important parts of plants, the genetic mechanisms of photosynthesis or the response of leaf to a single abiotic and biotic stress have been well studied. However, few researches have involved in the integration of data analysis from system level in leaf tissue under multiple abiotic stresses by utilizing biological networks. In this study, the weighted gene co-expression network analysis (WGCNA) strategy was used to integrate multiple data in leaf tissue of Populus species under different sample treatments. The gene co-expression networks were constructed and functional modules were identified by selecting the suitable soft threshold power β in the procedure of WGCNA. The identified hub genes and gene modules were annotated by agriGO, NetAffx Analysis Center, The Plant Genome Integrative Explorer (PlantGenIE) and other annotation tools. The annotation results have displayed that the highly correlated modules and hub genes are involved in the important biological processes or pathways related to module traits. The efficiency of the WGCNA strategy can generate comprehensive understanding of gene module-traits associations in leaf tissue, which will provide novel insight into the genetic mechanism of Populus species.

Integrating GWAS with a gene co-expression network better prioritizes candidate genes associated with root metaxylem phenes in maize

Root metaxylems are phenotypically diverse structures whose function is particularly important under drought stress. Significant research has dissected the genetic machinery underlying metaxylem phenotypes in dicots, but that of monocots are relatively underexplored. In maize (Zea mays), a robust pipeline integrated a genome-wide association study (GWAS) of root metaxylem phenes under well-watered and water-stress conditions with a gene co-expression network to prioritize the strongest gene candidates. We identified 244 candidate genes by GWAS, of which 103 reside in gene co-expression modules most relevant to xylem development. Several candidate genes may be involved in biosynthetic processes related to the cell wall, hormone signaling, oxidative stress responses, and drought responses. Of those, six gene candidates were detected in multiple root metaxylem phenes in both well-watered and water-stress conditions. We posit that candidate genes that are more essential to network function based on gene co-expression (i.e., hubs or bottlenecks) should be prioritized and classify 33 essential genes for further investigation. Our study demonstrates a new strategy for identifying promising gene candidates and presents several gene candidates that may enhance our understanding of vascular development and responses to drought in cereals.

Integrated analysis of transcriptomics and defense-related phytohormones to discover hub genes conferring maize Gibberella ear rot caused by Fusarium Graminearum

BACKGROUND

Gibberella ear rot (GER) is one of the most devastating diseases in maize growing areas, which directly reduces grain yield and quality. However, the underlying defense response of maize to pathogens infection is largely unknown.

RESULTS

To gain a comprehensive understanding of the defense response in GER resistance, two contrasting inbred lines 'Nov-82' and 'H10' were used to explore transcriptomic profiles and defense-related phytohormonal alterations during Fusarium graminearum infection. Transcriptomic analysis revealed 4,417 and 4,313 differentially expressed genes (DEGs) from the Nov-82 and H10, respectively, and 647 common DEGs between the two lines. More DEGs were obviously enriched in phenylpropanoid biosynthesis, secondary metabolites biosynthesis, metabolic process and defense-related pathways. In addition, the concentration of the defense-related phytohormones, jasmonates (JAs) and salicylates (SAs), was greatly induced after the pathogen infection. The level of JAs in H10 was more higher than in Nov-82, whereas an opposite pattern for the SA between the both lines. Integrated analysis of the DEGs and the phytohormones revealed five vital modules based on co-expression network analysis according to their correlation. A total of 12 hub genes encoding fatty acid desaturase, subtilisin-like protease, ethylene-responsive transcription factor, 1-aminocyclopropane-1-carboxylate oxidase, and sugar transport protein were captured from the key modules, indicating that these genes might play unique roles in response to pathogen infection, CONCLUSIONS: Overall, our results indicate that large number DEGs related to plant disease resistance and different alteration of defensive phytohormones were activated during F. graminearum infection, providing new insight into the defense response against pathogen invasion, in addition to the identified hub genes that can be further investigated for enhancing maize GER resistance.

Prioritizing Maize Metabolic Gene Regulators through Multi-Omic Network Integration

Elucidating gene regulatory networks is a major area of study within plant systems biology. Phenotypic traits are intricately linked to specific gene expression profiles. These expression patterns arise primarily from regulatory connections between sets of transcription factors (TFs) and their target genes. Here, we integrated 46 co-expression networks, 283 protein-DNA interaction (PDI) assays, and 16 million SNPs used to identify expression quantitative trait loci (eQTL) to construct TF-target networks. In total, we analyzed ∼4.6M interactions to generate four distinct types of TF-target networks: co-expression, PDI, trans -eQTL, and cis -eQTL combined with PDIs. To functionally annotate TFs based on their target genes, we implemented three different network integration strategies. We evaluated the effectiveness of each strategy through TF loss-of function mutant inspection and random network analyses. The multi-network integration allowed us to identify transcriptional regulators of several biological processes. Using the topological properties of the fully integrated network, we identified potential functionally redundant TF paralogs. Our findings retrieved functions previously documented for numerous TFs and revealed novel functions that are crucial for informing the design of future experiments. The approach here-described lays the foundation for the integration of multi-omic datasets in maize and other plant systems.

GRAPHICAL ABSTRACT

Genome-wide screening of meta-QTL and candidate genes controlling yield and yield-related traits in barley (Hordeum vulgare L.)

Increasing yield is an important goal of barley breeding. In this study, 54 papers published from 2001-2022 on QTL mapping for yield and yield-related traits in barley were collected, which contained 1080 QTLs mapped to the barley high-density consensus map for QTL meta-analysis. These initial QTLs were integrated into 85 meta-QTLs (MQTL) with a mean confidence interval (CI) of 2.76 cM, which was 7.86-fold narrower than the CI of the initial QTL. Among these 85 MQTLs, 68 MQTLs were validated in GWAS studies, and 25 breeder's MQTLs were screened from them. Seventeen barley orthologs of yield-related genes in rice and maize were identified within the hcMQTL region based on comparative genomics strategy and were presumed to be reliable candidates for controlling yield-related traits. The results of this study provide useful information for molecular marker-assisted breeding and candidate gene mining of yield-related traits in barley.

Deep learning the cis-regulatory code for gene expression in selected model plants

Elucidating the relationship between non-coding regulatory element sequences and gene expression is crucial for understanding gene regulation and genetic variation. We explored this link with the training of interpretable deep learning models predicting gene expression profiles from gene flanking regions of the plant species Arabidopsis thaliana, Solanum lycopersicum, Sorghum bicolor, and Zea mays. With over 80% accuracy, our models enabled predictive feature selection, highlighting e.g. the significant role of UTR regions in determining gene expression levels. The models demonstrated remarkable cross-species performance, effectively identifying both conserved and species-specific regulatory sequence features and their predictive power for gene expression. We illustrated the application of our approach by revealing causal links between genetic variation and gene expression changes across fourteen tomato genomes. Lastly, our models efficiently predicted genotype-specific expression of key functional gene groups, exemplified by underscoring known phenotypic and metabolic differences between Solanum lycopersicum and its wild, drought-resistant relative, Solanum pennellii.

Integrated multi-omics analyses and genome-wide association studies reveal prime candidate genes of metabolic and vegetative growth variation in canola

Genome-wide association studies (GWAS) identified thousands of genetic loci associated with complex plant traits, including many traits of agronomical importance. However, functional interpretation of GWAS results remains challenging because of large candidate regions due to linkage disequilibrium. High-throughput omics technologies, such as genomics, transcriptomics, proteomics and metabolomics open new avenues for integrative systems biological analyses and help to nominate systems information supported (prime) candidate genes. In the present study, we capitalise on a diverse canola population with 477 spring-type lines which was previously analysed by high-throughput phenotyping of growth-related traits and by RNA sequencing and metabolite profiling for multi-omics-based hybrid performance prediction. We deepened the phenotypic data analysis, now providing 123 time-resolved image-based traits, to gain insight into the complex relations during early vegetative growth and reanalysed the transcriptome data based on the latest Darmor-bzh v10 genome assembly. Genome-wide association testing revealed 61 298 robust quantitative trait loci (QTL) including 187 metabolite QTL, 56814 expression QTL and 4297 phenotypic QTL, many clustered in pronounced hotspots. Combining information about QTL colocalisation across omics layers and correlations between omics features allowed us to discover prime candidate genes for metabolic and vegetative growth variation. Prioritised candidate genes for early biomass accumulation include A06p05760.1_BnaDAR (PIAL1), A10p16280.1_BnaDAR, C07p48260.1_BnaDAR (PRL1) and C07p48510.1_BnaDAR (CLPR4). Moreover, we observed unequal effects of the Brassica A and C subgenomes on early biomass production.

Superior haplotypes of key drought-responsive genes reveal opportunities for the development of climate-resilient rice varieties

Haplotype-based breeding is an emerging and innovative concept that enables the development of designer crop varieties by exploiting and exploring superior alleles/haplotypes among target genes to create new traits in breeding programs. In this regard, whole-genome re-sequencing of 399 genotypes (landraces and breeding lines) from the 3000 rice genomes panel (3K-RG) is mined to identify the superior haplotypes for 95 drought-responsive candidate genes. Candidate gene-based association analysis reveals 69 marker-trait associations (MTAs) in 16 genes for single plant yield (SPY) under drought stress. Haplo-pheno analysis of these 16 genes identifies superior haplotypes for seven genes associated with the higher SPY under drought stress. Our study reveals that the performance of lines possessing superior haplotypes is significantly higher (p ≤ 0.05) as measured by single plant yield (SPY), for the OsGSK1-H4, OsDSR2-H3, OsDIL1-H22, OsDREB1C-H3, ASR3-H88, DSM3-H4 and ZFP182-H4 genes as compared to lines without the superior haplotypes. The validation results indicate that a superior haplotype for the DREB transcription factor (OsDREB1C) is present in all the drought-tolerant rice varieties, while it was notably absent in all susceptible varieties. These lines carrying the superior haplotypes can be used as potential donors in haplotype-based breeding to develop high-yielding drought-tolerant rice varieties.

RicePilaf: a post-GWAS/QTL dashboard to integrate pangenomic, coexpression, regulatory, epigenomic, ontology, pathway, and text-mining information to provide functional insights into rice QTLs and GWAS loci

BACKGROUND

As the number of genome-wide association study (GWAS) and quantitative trait locus (QTL) mappings in rice continues to grow, so does the already long list of genomic loci associated with important agronomic traits. Typically, loci implicated by GWAS/QTL analysis contain tens to hundreds to thousands of single-nucleotide polmorphisms (SNPs)/genes, not all of which are causal and many of which are in noncoding regions. Unraveling the biological mechanisms that tie the GWAS regions and QTLs to the trait of interest is challenging, especially since it requires collating functional genomics information about the loci from multiple, disparate data sources.

RESULTS

We present RicePilaf, a web app for post-GWAS/QTL analysis, that performs a slew of novel bioinformatics analyses to cross-reference GWAS results and QTL mappings with a host of publicly available rice databases. In particular, it integrates (i) pangenomic information from high-quality genome builds of multiple rice varieties, (ii) coexpression information from genome-scale coexpression networks, (iii) ontology and pathway information, (iv) regulatory information from rice transcription factor databases, (v) epigenomic information from multiple high-throughput epigenetic experiments, and (vi) text-mining information extracted from scientific abstracts linking genes and traits. We demonstrate the utility of RicePilaf by applying it to analyze GWAS peaks of preharvest sprouting and genes underlying yield-under-drought QTLs.

CONCLUSIONS

RicePilaf enables rice scientists and breeders to shed functional light on their GWAS regions and QTLs, and it provides them with a means to prioritize SNPs/genes for further experiments. The source code, a Docker image, and a demo version of RicePilaf are publicly available at https://github.com/bioinfodlsu/rice-pilaf.

Genome-wide family prediction unveils molecular mechanisms underlying the regulation of agronomic traits in Urochloa ruziziensis

Tropical forage grasses, particularly those belonging to the Urochloa genus, play a crucial role in cattle production and serve as the main food source for animals in tropical and subtropical regions. The majority of these species are apomictic and tetraploid, highlighting the significance of U. ruziziensis, a sexual diploid species that can be tetraploidized for use in interspecific crosses with apomictic species. As a means to support breeding programs, our study investigates the feasibility of genome-wide family prediction in U. ruziziensis families to predict agronomic traits. Fifty half-sibling families were assessed for green matter yield, dry matter yield, regrowth capacity, leaf dry matter, and stem dry matter across different clippings established in contrasting seasons with varying available water capacity. Genotyping was performed using a genotyping-by-sequencing approach based on DNA samples from family pools. In addition to conventional genomic prediction methods, machine learning and feature selection algorithms were employed to reduce the necessary number of markers for prediction and enhance predictive accuracy across phenotypes. To explore the regulation of agronomic traits, our study evaluated the significance of selected markers for prediction using a tree-based approach, potentially linking these regions to quantitative trait loci (QTLs). In a multiomic approach, genes from the species transcriptome were mapped and correlated to those markers. A gene coexpression network was modeled with gene expression estimates from a diverse set of U. ruziziensis genotypes, enabling a comprehensive investigation of molecular mechanisms associated with these regions. The heritabilities of the evaluated traits ranged from 0.44 to 0.92. A total of 28,106 filtered SNPs were used to predict phenotypic measurements, achieving a mean predictive ability of 0.762. By employing feature selection techniques, we could reduce the dimensionality of SNP datasets, revealing potential genotype-phenotype associations. The functional annotation of genes near these markers revealed associations with auxin transport and biosynthesis of lignin, flavonol, and folic acid. Further exploration with the gene coexpression network uncovered associations with DNA metabolism, stress response, and circadian rhythm. These genes and regions represent important targets for expanding our understanding of the metabolic regulation of agronomic traits and offer valuable insights applicable to species breeding. Our work represents an innovative contribution to molecular breeding techniques for tropical forages, presenting a viable marker-assisted breeding approach and identifying target regions for future molecular studies on these agronomic traits.

Integration of genome-wide association study, linkage analysis, and population transcriptome analysis to reveal the TaFMO1-5B modulating seminal root growth in bread wheat

Bread wheat, one of the keystone crops for global food security, is challenged by climate change and resource shortage. The root system plays a vital role in water and nutrient absorption, making it essential for meeting the growing global demand. Here, using an association-mapping population composed of 406 accessions, we identified QTrl.Rs-5B modulating seminal root development with a genome-wide association study and validated its genetic effects with two F5 segregation populations. Transcriptome-wide association study prioritized TaFMO1-5B, a gene encoding the flavin-containing monooxygenases, as the causal gene for QTrl.Rs-5B, whose expression levels correlate negatively with the phenotyping variations among our population. The lines silenced for TaFMO1-5B consistently showed significantly larger seminal roots in different genetic backgrounds. Additionally, the agriculture traits measured in multiple environments showed that QTrl.Rs-5B also affects yield component traits and plant architecture-related traits, and its favorable haplotype modulates these traits toward that of modern cultivars, suggesting the application potential of QTrl.Rs-5B for wheat breeding. Consistently, the frequency of the favorable haplotype of QTrl.Rs-5B increased with habitat expansion and breeding improvement of bread wheat. In conclusion, our findings identified and demonstrated the effects of QTrl.Rs-5B on seminal root development and illustrated that it is a valuable genetic locus for wheat root improvement.

Leveraging prior biological knowledge improves prediction of tocochromanols in maize grain

With an essential role in human health, tocochromanols are mostly obtained by consuming seed oils; however, the vitamin E content of the most abundant tocochromanols in maize (Zea mays L.) grain is low. Several large-effect genes with cis-acting variants affecting messenger RNA (mRNA) expression are mostly responsible for tocochromanol variation in maize grain, with other relevant associated quantitative trait loci (QTL) yet to be fully resolved. Leveraging existing genomic and transcriptomic information for maize inbreds could improve prediction when selecting for higher vitamin E content. Here, we first evaluated a multikernel genomic best linear unbiased prediction (MK-GBLUP) approach for modeling known QTL in the prediction of nine tocochromanol grain phenotypes (12-21 QTL per trait) within and between two panels of 1,462 and 242 maize inbred lines. On average, MK-GBLUP models improved predictive abilities by 7.0-13.6% when compared with GBLUP. In a second approach with a subset of 545 lines from the larger panel, the highest average improvement in predictive ability relative to GBLUP was achieved with a multi-trait GBLUP model (15.4%) that had a tocochromanol phenotype and transcript abundances in developing grain for a few large-effect candidate causal genes (1-3 genes per trait) as multiple response variables. Taken together, our study illustrates the enhancement of prediction models when informed by existing biological knowledge pertaining to QTL and candidate causal genes.

Lipid Metabolism and Improvement in Oilseed Crops: Recent Advances in Multi-Omics Studies

Oilseed crops are rich in plant lipids that not only provide essential fatty acids for the human diet but also play important roles as major sources of biofuels and indispensable raw materials for the chemical industry. The regulation of lipid metabolism genes is a major factor affecting oil production. In this review, we systematically summarize the metabolic pathways related to lipid production and storage in plants and highlight key research advances in characterizing the genes and regulatory factors influencing lipid anabolic metabolism. In addition, we integrate the latest results from multi-omics studies on lipid metabolism to provide a reference to better understand the molecular mechanisms underlying oil anabolism in oilseed crops.

Genetic architecture of ear traits based on association mapping and co-expression networks in maize inbred lines and hybrids

Ear traits are key contributors to grain yield in maize; therefore, exploring their genetic basis facilitates the improvement of grain yield. However, the underlying molecular mechanisms of ear traits remain obscure in both inbred lines and hybrids. Here, two association panels, respectively, comprising 203 inbred lines (IP) and 246 F1 hybrids (HP) were employed to identify candidate genes for six ear traits. The IP showed higher phenotypic variation and lower phenotypic mean than the HP for all traits, except ear tip-barrenness length. By conducting a genome-wide association study (GWAS) across multiple environments, 101 and 228 significant single-nucleotide polymorphisms (SNPs) associated with six ear traits were identified in the IP and HP, respectively. Of these significant SNPs identified in the HP, most showed complete-incomplete dominance and over-dominance effects for each ear trait. Combining a gene co-expression network with GWAS results, 186 and 440 candidate genes were predicted in the IP and HP, respectively, including known ear development genes ids1 and sid1. Of these, nine candidate genes were detected in both populations and expressed in maize ear and spikelet tissues. Furthermore, two key shared genes (GRMZM2G143330 and GRMZM2G171139) in both populations were found to be significantly associated with ear traits in the maize Goodman diversity panel with high-density variations. These findings advance our knowledge of the genetic architecture of ear traits between inbred lines and hybrids and provide a valuable resource for the genetic improvement of ear traits in maize.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01426-9.

Hybrid performance evaluation and genome-wide association analysis of root system architecture in a maize association population

KEY MESSAGE

The genetic architecture of RSA traits was dissected by GWAS and coexpression networks analysis in a maize association population. Root system architecture (RSA) is a crucial determinant of water and nutrient uptake efficiency in crops. However, the maize genetic architecture of RSA is still poorly understood due to the challenges in quantifying root traits and the lack of dense molecular markers. Here, an association mapping panel including 356 inbred lines were crossed with a common tester, Zheng58, and the test crosses were phenotyped for 12 RSA traits in three locations. We observed a 1.3 ~ sixfold phenotypic variation for measured RSA in the association panel. The association panel consisted of four subpopulations, non-stiff stalk (NSS) lines, stiff stalk (SS), tropical/subtropical (TST), and mixed. Zheng58 × TST has a 2.1% higher crown root number (CRN) and 8.6% less brace root number (BRN) than Zheng58 × NSS and Zheng58 × SS, respectively. Using a genome-wide association study (GWAS) with 1.25 million SNPs and correction for population structure, 191 significant SNPs were identified for root traits. Ninety (47%) of the significant SNPs showed positive allelic effects, and 101 (53%) showed negative effects. Each locus could explain 0.39% to 11.8% of phenotypic variation. By integrating GWAS results and comparing coexpression networks, 26 high-priority candidate genes were identified. Gene GRMZM2G377215, which belongs to the COBRA-like gene family, affected root growth and development. Gene GRMZM2G468657 encodes the aspartic proteinase nepenthesin-1, related to root development and N-deficient response. Collectively, our research provides progress in the genetic dissection of root system architecture. These findings present the further possibility for the genetic improvement of root traits in maize.

Genome-wide association study of soybean (Glycine max [L.] Merr.) germplasm for dissecting the quantitative trait nucleotides and candidate genes underlying yield-related traits

Soybean (Glycine max [L.] Merr.) is one of the most significant crops in the world in terms of oil and protein. Owing to the rising demand for soybean products, there is an increasing need for improved varieties for more productive farming. However, complex correlation patterns among quantitative traits along with genetic interactions pose a challenge for soybean breeding. Association studies play an important role in the identification of accession with useful alleles by locating genomic sites associated with the phenotype in germplasm collections. In the present study, a genome-wide association study was carried out for seven agronomic and yield-related traits. A field experiment was conducted in 2015/2016 at two locations that include 155 diverse soybean germplasm. These germplasms were genotyped using SoySNP50K Illumina Infinium Bead-Chip. A total of 51 markers were identified for node number, plant height, pods per plant, seeds per plant, seed weight per plant, hundred-grain weight, and total yield using a multi-locus linear mixed model (MLMM) in FarmCPU. Among these significant SNPs, 18 were putative novel QTNs, while 33 co-localized with previously reported QTLs. A total of 2,356 genes were found in 250 kb upstream and downstream of significant SNPs, of which 17 genes were functional and the rest were hypothetical proteins. These 17 candidate genes were located in the region of 14 QTNs, of which ss715580365, ss715608427, ss715632502, and ss715620131 are novel QTNs for PH, PPP, SDPP, and TY respectively. Four candidate genes, Glyma.01g199200, Glyma.10g065700, Glyma.18g297900, and Glyma.14g009900, were identified in the vicinity of these novel QTNs, which encode lsd one like 1, Ergosterol biosynthesis ERG4/ERG24 family, HEAT repeat-containing protein, and RbcX2, respectively. Although further experimental validation of these candidate genes is required, several appear to be involved in growth and developmental processes related to the respective agronomic traits when compared with their homologs in Arabidopsis thaliana. This study supports the usefulness of association studies and provides valuable data for functional markers and investigating candidate genes within a diverse germplasm collection in future breeding programs.

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.

Genome wide association and haplotype analyses for the crease depth trait in bread wheat (Triticum aestivum L.)

Wheat grain has a complex structure that includes a crease on one side, and tissues within the crease region play an important role in nutrient transportation during wheat grain development. However, the genetic architecture of the crease region is still unclear. In this study, 413 global wheat accessions were resequenced and a method was developed for evaluating the phenotypic data of crease depth (CD). The CD values exhibited continuous and considerable large variation in the population, and the broad-sense heritability was 84.09%. CD was found to be positively correlated with grain-related traits and negatively with quality-related traits. Analysis of differentiation of traits between landraces and cultivars revealed that grain-related traits and CD were simultaneously improved during breeding improvement. Moreover, 2,150.8-Mb genetic segments were identified to fall within the selective sweeps between the landraces and cultivars; they contained some known functional genes for quality- and grain-related traits. Genome-wide association study (GWAS) was performed using around 10 million SNPs generated by genome resequencing and 551 significant SNPs and 18 QTLs were detected significantly associated with CD. Combined with cluster analysis of gene expression, haplotype analysis, and annotated information of candidate genes, two promising genes TraesCS3D02G197700 and TraesCS5A02G292900 were identified to potentially regulate CD. To the best of our knowledge, this is the first study to provide the genetic basis of CD, and the genetic loci identified in this study may ultimately assist in wheat breeding programs.

Design, execution, and interpretation of plant RNA-seq analyses

Genomics has transformed our understanding of the genetic architecture of traits and the genetic variation present in plants. Here, we present a review of how RNA-seq can be performed to tackle research challenges addressed by plant sciences. We discuss the importance of experimental design in RNA-seq, including considerations for sampling and replication, to avoid pitfalls and wasted resources. Approaches for processing RNA-seq data include quality control and counting features, and we describe common approaches and variations. Though differential gene expression analysis is the most common analysis of RNA-seq data, we review multiple methods for assessing gene expression, including detecting allele-specific gene expression and building co-expression networks. With the production of more RNA-seq data, strategies for integrating these data into genetic mapping pipelines is of increased interest. Finally, special considerations for RNA-seq analysis and interpretation in plants are needed, due to the high genome complexity common across plants. By incorporating informed decisions throughout an RNA-seq experiment, we can increase the knowledge gained.

Mining genes regulating root system architecture in maize based on data integration analysis

KEY MESSAGE

A new strategy that integrated multiple public data resources was established to construct root gene co-expression network and mine genes regulating root system architecture in maize. A root gene co-expression network, containing 13,874 genes, was constructed. A total of 53 root hub genes and 16 priority root candidate genes were identified. One priority root candidate was further functionally verified using overexpression transgenic maize lines. Root system architecture (RSA) is crucial for crops productivity and stress tolerance. In maize, few RSA genes are functionally cloned, and effective discovery of RSA genes remains a great of challenge. In this work, we established a strategy to mine maize RSA genes by integrating functionally characterized root genes, root transcriptome, weighted gene co-expression network analysis (WGCNA) and genome-wide association analysis (GWAS) of RSA traits based on public data resources. A total of 589 maize root genes were collected by searching well-characterized root genes in maize or homologous genes of other species. We performed WGCNA to construct a maize root gene co-expression network containing 13874 genes based on public available root transcriptome data, and further discovered the 53 hub genes related to root traits. In addition, by the prediction function of obtained root gene co-expression network, a total of 1082 new root candidate genes were explored. By further overlapping the obtained new root candidate gene with the root-related GWAS of RSA candidate genes, 16 priority root candidate genes were identified. Finally, a priority root candidate gene, Zm00001d023379 (encodes pyruvate kinase 2), was validated to modulate root open angle and shoot-borne roots number using its overexpression transgenic lines. Our results develop an integration analysis method for effectively exploring regulatory genes of RSA in maize and open a new avenue to mine the candidate genes underlying complex traits.

Charting plant gene functions in the multi-omics and single-cell era

Despite the increased access to high-quality plant genome sequences, the set of genes with a known function remains far from complete. With the advent of novel bulk and single-cell omics profiling methods, we are entering a new era where advanced and highly integrative functional annotation strategies are being developed to elucidate the functions of all plant genes. Here, we review different multi-omics approaches to improve functional and regulatory gene characterization and highlight the power of machine learning and network biology to fully exploit the complementary information embedded in different omics layers. Finally, we discuss the potential of emerging single-cell methods and algorithms to further increase the resolution, allowing generation of functional insights about plant biology.

Genetic architecture of protein and oil content in soybean seed and meal

Soybean is grown primarily for the protein and oil extracted from its seed and its value is influenced by these components. The objective of this study was to map marker-trait associations (MTAs) for the concentration of seed protein, oil, and meal protein using the soybean nested association mapping (SoyNAM) population. The composition traits were evaluated on seed harvested from over 5000 inbred lines of the SoyNAM population grown in 10 field locations across 3 years. Estimated heritabilities were at least 0.85 for all three traits. The genotyping of lines with single nucleotide polymorphism markers resulted in the identification of 107 MTAs for the three traits. When MTAs for the three traits that mapped within 5 cM intervals were binned together, the MTAs were mapped to 64 intervals on 19 of the 20 soybean chromosomes. The majority of the MTA effects were small and of the 107 MTAs, 37 were for protein content, 39 for meal protein, and 31 for oil content. For cases where a protein and oil MTAs mapped to the same interval, most (94%) significant effects were opposite for the two traits, consistent with the negative correlation between these traits. A coexpression analysis identified candidate genes linked to MTAs and 18 candidate genes were identified. The large number of small effect MTAs for the composition traits suggest that genomic prediction would be more effective in improving these traits than marker-assisted selection.

Identification of hub genes regulating isoflavone accumulation in soybean seeds via GWAS and WGCNA approaches

INTRODUCTION

Isoflavones are the secondary metabolites synthesized by the phenylpropanoid biosynthesis pathway in soybean that benefits human and plant health.

METHODS

In this study, we have profiled seed isoflavone content by HPLC in 1551 soybean accessions grown in Beijing and Hainan for two consecutive years (2017 and 2018) and in Anhui for one year (2017).

RESULTS

A broad range of phenotypic variations was observed for individual and total isoflavone (TIF) content. The TIF content ranged from 677.25 to 5823.29 µg g^-1 in the soybean natural population. Using a genome-wide association study (GWAS) based on 6,149,599 single nucleotide polymorphisms (SNPs), we identified 11,704 SNPs significantly associated with isoflavone contents; 75% of them were located within previously reported QTL regions for isoflavone. Two significant regions on chromosomes 5 and 11 were associated with TIF and malonylglycitin across more than 3 environments. Furthermore, the WGCNA identified eight key modules: black, blue, brown, green, magenta, pink, purple, and turquoise. Of the eight co-expressed modules, brown (r = 0.68***), magenta (r = 0.64***), and green (r = 0.51**) showed a significant positive association with TIF, as well as with individual isoflavone contents. By combining the gene significance, functional annotation, and enrichment analysis information, four hub genes Glyma.11G108100, Glyma.11G107100, Glyma.11G106900, and Glyma.11G109100 encoding, basic-leucine zipper (bZIP) transcription factor, MYB4 transcription factor, early responsive to dehydration, and PLATZ transcription factor respectively were identified in brown and green modules. The allelic variation in Glyma.11G108100 significantly influenced individual and TIF accumulation.

DISCUSSION

The present study demonstrated that the GWAS approach, combined with WGCNA, could efficiently identify isoflavone candidate genes in the natural soybean population.

Genome wide association analysis of root hair traits in rice reveals novel genomic regions controlling epidermal cell differentiation

BACKGROUND

Genome wide association (GWA) studies demonstrate linkages between genetic variants and traits of interest. Here, we tested associations between single nucleotide polymorphisms (SNPs) in rice (Oryza sativa) and two root hair traits, root hair length (RHL) and root hair density (RHD). Root hairs are outgrowths of single cells on the root epidermis that aid in nutrient and water acquisition and have also served as a model system to study cell differentiation and tip growth. Using lines from the Rice Diversity Panel-1, we explored the diversity of root hair length and density across four subpopulations of rice (aus, indica, temperate japonica, and tropical japonica). GWA analysis was completed using the high-density rice array (HDRA) and the rice reference panel (RICE-RP) SNP sets.

RESULTS

We identified 18 genomic regions related to root hair traits, 14 of which related to RHD and four to RHL. No genomic regions were significantly associated with both traits. Two regions overlapped with previously identified quantitative trait loci (QTL) associated with root hair density in rice. We identified candidate genes in these regions and present those with previously published expression data relevant to root hair development. We re-phenotyped a subset of lines with extreme RHD phenotypes and found that the variation in RHD was due to differences in cell differentiation, not cell size, indicating genes in an associated genomic region may influence root hair cell fate. The candidate genes that we identified showed little overlap with previously characterized genes in rice and Arabidopsis.

CONCLUSIONS

Root hair length and density are quantitative traits with complex and independent genetic control in rice. The genomic regions described here could be used as the basis for QTL development and further analysis of the genetic control of root hair length and density. We present a list of candidate genes involved in root hair formation and growth in rice, many of which have not been previously identified as having a relation to root hair growth. Since little is known about root hair growth in grasses, these provide a guide for further research and crop improvement.

Metabolite-Based Genome-Wide Association Studies of Large-Scale Metabolome Analysis to Illustrate Alterations in the Metabolite Landscape of Plants upon Responses to Stresses

Given that anthropogenic activities are evoking a profound effect on the climate resulting in more extreme events such as severe drought and heat waves while global demand for food is ever-increasing, understanding plant responses to stresses is critical. As metabolites are fundamental for plant growth regulation and plant lifespan and an important component of yield, illustrating how the metabolite landscape of plant changes following stress will supply important clues as to how to improve the plant resistance to stress. Recently, billions of single-nucleotide polymorphisms (SNPs) have been obtained and used to identify the associations between genetic variants of genomes and relevant crop agronomic traits through different genetic methods such as genome-wide association studies (GWAS). Therefore, in this chapter, we provide comprehensive guidelines concerning the experimental design, metabolite profiling, and metabolite-based genome-wide association studies (mGWAS) of large-scale metabolome analysis to accelerate the future identification of the valuable stress-resistant genes and metabolites.

Leaf physiology variations are modulated by natural variations that underlie stomatal morphology in Populus

Stomata are essential for photosynthesis and abiotic stress tolerance. Here, we used multiomics approaches to dissect the genetic architecture and adaptive mechanisms that underlie stomatal morphology in Populus tomentosa juvenile natural population (303 accessions). We detected 46 candidate genes and 15 epistatic gene-pairs, associated with 5 stomatal morphologies and 18 leaf development and photosynthesis traits, through genome-wide association studies. Expression quantitative trait locus mapping revealed that stomata-associated gene loci were significantly associated with the expression of leaf-related genes; selective sweep analysis uncovered significant differentiation in the allele frequencies of genes that underlie stomatal variations. An allelic regulatory network operating under drought stress and adequate precipitation conditions, with three key regulators (DUF538, TRA2 and AbFH2) and eight interacting genes, was identified that might regulate leaf physiology via modulation of stomatal shape and density. Validation of candidate gene variations in drought-tolerant and F₁ hybrid populations of P. tomentosa showed that the DUF538, TRA2 and AbFH2 loci cause functional stabilisation of spatiotemporal regulatory, whose favourable alleles can be faithfully transmitted to offspring. This study provides insights concerning leaf physiology and stress tolerance via the regulation of stomatal determination in perennial plants.

A multi-omics integrative network map of maize

Networks are powerful tools to uncover functional roles of genes in phenotypic variation at a system-wide scale. Here, we constructed a maize network map that contains the genomic, transcriptomic, translatomic and proteomic networks across maize development. This map comprises over 2.8 million edges in more than 1,400 functional subnetworks, demonstrating an extensive network divergence of duplicated genes. We applied this map to identify factors regulating flowering time and identified 2,651 genes enriched in eight subnetworks. We validated the functions of 20 genes, including 18 with previously unknown connections to flowering time in maize. Furthermore, we uncovered a flowering pathway involving histone modification. The multi-omics integrative network map illustrates the principles of how molecular networks connect different types of genes and potential pathways to map a genome-wide functional landscape in maize, which should be applicable in a wide range of species.

Identification of genes for seed isoflavones based on bulk segregant analysis sequencing in soybean natural population

We identified four hub genes for isoflavone biosynthesis based on BSA-seq and WGCNA methods and validated that GmIE3-1 positively contribute to isoflavone accumulation in soybean. Soybean isoflavones are secondary metabolites of great interest owing to their beneficial impact on human health. Herein, we profiled the seed isoflavone content by HPLC in 1551 soybean accessions grown in two locations for two years and constructed two extreme pools with high (4065.1 µg g^-1) and low (1427.23 µg g^-1) isoflavone contents to identify candidate genes involved in isoflavone biosynthesis pathways using bulk segregant analysis sequencing (BSA-seq) approach. The results showed that the average sequencing depths were 50.3× and 65.7× in high and low pools, respectively. A total of 23,626 polymorphic SNPs and 5299 InDels were detected between both pools and 1492 genes with different variations were identified. Based on differential genes in BSA-seq and weighted gene co-expression network analysis (WGCNA), four hub genes, Glyma.06G290400 (designated as GmIE3-1), Glyma.01G239200, Glyma.01G241500, Glyma.13G256100 were identified, encoding E3 ubiquitin-protein ligase, arm repeat protein interacting with ABF2, zinc metallopeptidase EGY3, and dynamin-related protein 3A, respectively. The allelic variation in GmIE3-1 showed a significant influence on isoflavone accumulation. The virus-induced gene silencing (VIGS) and RNAi hairy root transformation of GmIE3-1 revealed partial suppression of this gene could cause a significant decrease (P < 0.0001) of total isoflavone content, suggesting GmIE3-1 is a positive regulator for isoflavones. The present study demonstrated that the BSA-seq approach combined with WGCNA, VIGS and hairy root transformation can efficiently identify isoflavone candidate genes in soybean natural population.

Co-expression pan-network reveals genes involved in complex traits within maize pan-genome

BACKGROUND

With the advances in the high throughput next generation sequencing technologies, genome-wide association studies (GWAS) have identified a large set of variants associated with complex phenotypic traits at a very fine scale. Despite the progress in GWAS, identification of genotype-phenotype relationship remains challenging in maize due to its nature with dozens of variants controlling the same trait. As the causal variations results in the change in expression, gene expression analyses carry a pivotal role in unraveling the transcriptional regulatory mechanisms behind the phenotypes.

RESULTS

To address these challenges, we incorporated the gene expression and GWAS-driven traits to extend the knowledge of genotype-phenotype relationships and transcriptional regulatory mechanisms behind the phenotypes. We constructed a large collection of gene co-expression networks and identified more than 2 million co-expressing gene pairs in the GWAS-driven pan-network which contains all the gene-pairs in individual genomes of the nested association mapping (NAM) population. We defined four sub-categories for the pan-network: (1) core-network contains the highest represented ~ 1% of the gene-pairs, (2) near-core network contains the next highest represented 1-5% of the gene-pairs, (3) private-network contains ~ 50% of the gene pairs that are unique to individual genomes, and (4) the dispensable-network contains the remaining 50-95% of the gene-pairs in the maize pan-genome. Strikingly, the private-network contained almost all the genes in the pan-network but lacked half of the interactions. We performed gene ontology (GO) enrichment analysis for the pan-, core-, and private- networks and compared the contributions of variants overlapping with genes and promoters to the GWAS-driven pan-network.

CONCLUSIONS

Gene co-expression networks revealed meaningful information about groups of co-regulated genes that play a central role in regulatory processes. Pan-network approach enabled us to visualize the global view of the gene regulatory network for the studied system that could not be well inferred by the core-network alone.

Genome-wide dissection of changes in maize root system architecture during modern breeding

Appropriate root system architecture (RSA) can improve maize yields in densely planted fields, but little is known about its genetic basis in maize. Here we performed root phenotyping of 14,301 field-grown plants from an association mapping panel to study the genetic architecture of maize RSA. A genome-wide association study identified 81 high-confidence RSA-associated candidate genes and revealed that 28 (24.3%) of known root-related genes were selected during maize domestication and improvement. We found that modern maize breeding has selected for a steeply angled root system. Favourable alleles related to steep root system angle have continuously accumulated over the course of modern breeding, and our data pinpoint the root-related genes that have been selected in different breeding eras. We confirm that two auxin-related genes, ZmRSA3.1 and ZmRSA3.2, contribute to the regulation of root angle and depth in maize. Our genome-wide identification of RSA-associated genes provides new strategies and genetic resources for breeding maize suitable for high-density planting.

eQTLs play critical roles in regulating gene expression and identifying key regulators in rice

The regulation of gene expression plays an essential role in both the phenotype and adaptation of plants. Transcriptome sequencing enables simultaneous identification of exonic variants and quantification of gene expression. Here, we sequenced the leaf transcriptomes of 287 rice accessions from around the world and obtained a total of 177 853 high-quality single nucleotide polymorphisms after filtering. Genome-wide association study identified 44 354 expression quantitative trait loci (eQTLs), which regulate the expression of 13 201 genes, as well as 17 local eQTL hotspots and 96 distant eQTL hotspots. Furthermore, a transcriptome-wide association study screened 21 candidate genes for starch content in the flag leaves at the heading stage. HS002 was identified as a significant distant eQTL hotspot with five downstream genes enriched for diterpene antitoxin synthesis. Co-expression analysis, eQTL analysis, and linkage mapping together demonstrated that bHLH026 acts as a key regulator to activate the expression of downstream genes. The transgenic assay revealed that bHLH026 is an important regulator of diterpenoid antitoxin synthesis and enhances the disease resistance of rice. These findings improve our knowledge of the regulatory mechanisms of gene expression variation and complex regulatory networks of the rice genome and will facilitate genetic improvement of cultivated rice varieties.

Mining synergistic genes for nutrient utilization and disease resistance in maize based on co-expression network and consensus QTLs

Nutrient restrictions and large-scale emergence of diseases are threatening the maize production. Recent findings demonstrated that there is a certain synergistic interaction between nutrition and diseases pathways in model plants, however there are few studies on the synergistic genes of nutrients and diseases in maize. Thus, the transcriptome data of nitrogen (N) and phosphorus (P) nutrients and diseases treatments in maize, rice, wheat and Arabidopsis thaliana were collected in this study, and four and 22 weighted co-expression modules were obtained by using Weighted Gene Co-expression Network Analysis (WGCNA) in leaf and root tissues, respectively. With a total of 5252 genes, MFUZZ cluster analysis screened 26 clusters with the same expression trend under nutrition and disease treatments. In the meantime, 1427 genes and 22 specific consensus quantitative trait loci (scQTLs) loci were identified by meta-QTL analysis of nitrogen and phosphorus nutrition and disease stress in maize. Combined with the results of cluster analysis and scQTLs, a total of 195 consistent genes were screened, of which six genes were shown to synergistically respond to nutrition and disease both in roots and leaves. Moreover, the six candidate genes were found in scQTLs associated with gray leaf spot (GLS) and corn leaf blight (CLB). In addition, subcellular localization and bioinformatics analysis of the six candidate genes revealed that they were primarily expressed in endoplasmic reticulum, mitochondria, nucleus and plasma membrane, and were involved in defense and stress, MeJA and abscisic acid response pathways. The fluorescence quantitative PCR confirmed their responsiveness to nitrogen and phosphorus nutrition as well as GLS treatments. Taken together, findings of this study indicated that the nutrition and disease have a significant synergistic response in maize.

Transcriptional responses to gibberellin in the maize tassel and control by DELLA domain proteins

The plant hormone gibberellin (GA) impacts plant growth and development differently depending on the developmental context. In the maize (Zea mays) tassel, application of GA alters floral development, resulting in the persistence of pistils. GA signaling is achieved by the GA-dependent turnover of DELLA domain transcription factors, encoded by dwarf8 (d8) and dwarf9 (d9) in maize. The D8-Mpl and D9-1 alleles disrupt GA signaling, resulting in short plants and normal tassel floret development in the presence of excess GA. However, D9-1 mutants are unable to block GA-induced pistil development. Gene expression in developing tassels of D8-Mpl and D9-1 mutants and their wild-type siblings was determined upon excess GA₃ and mock treatments. Using GA-sensitive transcripts as reporters of GA signaling, we identified a weak loss of repression under mock conditions in both mutants, with the effect in D9-1 being greater. D9-1 was also less able to repress GA signaling in the presence of excess GA₃ . We treated a diverse set of maize inbred lines with excess GA₃ and measured the phenotypic consequences on multiple aspects of development (e.g., height and pistil persistence in tassel florets). Genotype affected all GA-regulated phenotypes but there was no correlation between any of the GA-affected phenotypes, indicating that the complexity of the relationship between GA and development extends beyond the two-gene epistasis previously demonstrated for GA and brassinosteroid biosynthetic mutants.

GWAS and Transcriptome Analysis Reveal Key Genes Affecting Root Growth under Low Nitrogen Supply in Maize

Nitrogen (N) is one of the most important factors affecting crop production. Root morphology exhibits a high degree of plasticity to nitrogen deficiency. However, the mechanisms underlying the root foraging response under low-N conditions remain poorly understood. In this study, we analyzed 213 maize inbred lines using hydroponic systems and regarding their natural variations in 22 root traits and 6 shoot traits under normal (2 mM nitrate) and low-N (0 mM nitrate) conditions. Substantial phenotypic variations were detected for all traits. N deficiency increased the root length and decreased the root diameter and shoot related traits. A total of 297 significant marker-trait associations were identified by a genome-wide association study involving different N levels and the N response value. A total of 51 candidate genes with amino acid variations in coding regions or differentially expressed under low nitrogen conditions were identified. Furthermore, a candidate gene ZmNAC36 was resequenced in all tested lines. A total of 38 single nucleotide polymorphisms and 12 insertions and deletions were significantly associated with lateral root length of primary root, primary root length between 0 and 0.5 mm in diameter, primary root surface area, and total length of primary root under a low-N condition. These findings help us to improve our understanding of the genetic mechanism of root plasticity to N deficiency, and the identified loci and candidate genes will be useful for the genetic improvement of maize tolerance cultivars to N deficiency.

NetREx: Network-based Rice Expression Analysis Server for abiotic stress conditions

Recent focus on transcriptomic studies in food crops like rice, wheat and maize provide new opportunities to address issues related to agriculture and climate change. Re-analysis of such data available in public domain supplemented with annotations across molecular hierarchy can be of immense help to the plant research community, particularly co-expression networks representing transcriptionally coordinated genes that are often part of the same biological process. With this objective, we have developed NetREx, a Network-based Rice Expression Analysis Server, that hosts ranked co-expression networks of Oryza sativa using publicly available messenger RNA sequencing data across uniform experimental conditions. It provides a range of interactable data viewers and modules for analysing user-queried genes across different stress conditions (drought, flood, cold and osmosis) and hormonal treatments (abscisic and jasmonic acid) and tissues (root and shoot). Subnetworks of user-defined genes can be queried in pre-constructed tissue-specific networks, allowing users to view the fold change, module memberships, gene annotations and analysis of their neighbourhood genes and associated pathways. The web server also allows querying of orthologous genes from Arabidopsis, wheat, maize, barley and sorghum. Here, we demonstrate that NetREx can be used to identify novel candidate genes and tissue-specific interactions under stress conditions and can aid in the analysis and understanding of complex phenotypes linked to stress response in rice. Database URL: https://bioinf.iiit.ac.in/netrex/index.html.

Protecting plant nutrition from the effects of climate change

Climate change adversely affects plant nutrition, which serves as a major hurdle in the production of enough nutritious food to meet the needs of the growing global population. Here, we discuss how various climatic stressors impact nutrient homeostasis and how natural variation studies can yield resilient crop production systems to ensure future food security.

Genome-wide mediation analysis: an empirical study to connect phenotype with genotype via intermediate transcriptomic data in maize

Mapping genotype to phenotype is an essential topic in genetics and genomics research. As the Omics data become increasingly available, 2-variable methods have been widely applied to associate genotype with the phenotype (genome-wide association study), gene expression with the phenotype (transcriptome-wide association study), and genotype with gene expression. However, signals detected by these 2-variable association methods suffer from low mapping resolution or inexplicit causality between genotype and phenotype, making it challenging to interpret and validate the molecular mechanisms of the underlying genomic variations and the candidate genes. Under the context of genetics research, we hypothesized a causal chain from genotype to phenotype partially mediated by intermediate molecular processes, i.e. gene expression. To test this hypothesis, we applied the high-dimensional mediation analysis, a class of causal inference method with an assumed causal chain from the exposure to the mediator to the outcome, and implemented it with a maize association panel (N = 280 lines). Using 40 publicly available agronomy traits, 66 newly generated metabolite traits, and published RNA-seq data from 7 different tissues, our empirical study detected 736 unique mediating genes. Noticeably, 83/736 (11%) genes were identified in mediating more than 1 trait, suggesting the prevalence of pleiotropic mediating effects. We demonstrated that several identified mediating genes are consistent with their known functions. In addition, our results provided explicit hypotheses for functional validation and suggested that the mediation analysis is a powerful tool to integrate Omics data to connect genotype to phenotype.

Multi-omics analysis dissects the genetic architecture of seed coat content in Brassica napus

Background

Brassica napus is an important vegetable oil source worldwide. Seed coat content is a complex quantitative trait that negatively correlates with the seed oil content in B. napus.

Results

Here we provide insights into the genetic basis of natural variation of seed coat content by transcriptome-wide association studies (TWAS) and genome-wide association studies (GWAS) using 382 B. napus accessions. By population transcriptomic analysis, we identify more than 700 genes and four gene modules that are significantly associated with seed coat content. We also characterize three reliable quantitative trait loci (QTLs) controlling seed coat content by GWAS. Combining TWAS and correlation networks of seed coat content-related gene modules, we find that BnaC07.CCR-LIKE (CCRL) and BnaTT8s play key roles in the determination of the trait by modulating lignin biosynthesis. By expression GWAS analysis, we identify a regulatory hotspot on chromosome A09, which is involved in controlling seed coat content through BnaC07.CCRL and BnaTT8s. We then predict the downstream genes regulated by BnaTT8s using multi-omics datasets. We further experimentally validate that BnaCCRL and BnaTT8 positively regulate seed coat content and lignin content. BnaCCRL represents a novel identified gene involved in seed coat development. Furthermore, we also predict the key genes regulating carbon allocation between phenylpropane compounds and oil during seed development in B. napus.

Conclusions

This study helps us to better understand the complex machinery of seed coat development and provides a genetic resource for genetic improvement of seed coat content in B. napus breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02647-5.

Multiomics approach reveals a role of translational machinery in shaping maize kernel amino acid composition

Maize (Zea mays) seeds are a good source of protein, despite being deficient in several essential amino acids. However, eliminating the highly abundant but poorly balanced seed storage proteins has revealed that the regulation of seed amino acids is complex and does not rely on only a handful of proteins. In this study, we used two complementary omics-based approaches to shed light on the genes and biological processes that underlie the regulation of seed amino acid composition. We first conducted a genome-wide association study to identify candidate genes involved in the natural variation of seed protein-bound amino acids. We then used weighted gene correlation network analysis to associate protein expression with seed amino acid composition dynamics during kernel development and maturation. We found that almost half of the proteome was significantly reduced during kernel development and maturation, including several translational machinery components such as ribosomal proteins, which strongly suggests translational reprogramming. The reduction was significantly associated with a decrease in several amino acids, including lysine and methionine, pointing to their role in shaping the seed amino acid composition. When we compared the candidate gene lists generated from both approaches, we found a nonrandom overlap of 80 genes. A functional analysis of these genes showed a tight interconnected cluster dominated by translational machinery genes, especially ribosomal proteins, further supporting the role of translation dynamics in shaping seed amino acid composition. These findings strongly suggest that seed biofortification strategies that target the translation machinery dynamics should be considered and explored further.

Genome-Wide Association Study of Root System Architecture in Maize

Roots are important plant organs for the absorption of water and nutrients. To date, there have been few genome-wide association studies of maize root system architecture (RSA) in the field. The genetic basis of maize RSA is poorly understood, and the maize RSA-related genes that have been cloned are very limited. Here, 421 maize inbred lines of an association panel were planted to measure the root systems at the maturity stage, and a genome-wide association study was performed. There was a strong correlation among eight RSA traits, and the RSA traits were highly correlated with the aboveground plant architecture traits (e.g., plant height and ear leaf length, r = 0.13–0.25, p < 0.05). The RSA traits of the stiff stalk subgroup (SS) showed lower values than those of the non-stiff stalk subgroup (NSS) and tropical/subtropical subgroup (TST). Using the RSA traits, the genome-wide association study identified 63 SNPs and 189 candidate genes. Among them, nine candidate genes co-localized between RSA and aboveground architecture traits. A further co-expression analysis identified 88 candidate genes having high confidence levels. Furthermore, we identified four highly reliable RSA candidate genes, GRMZM2G099797, GRMZM2G354338, GRMZM2G085042, and GRMZM5G812926. This research provides theoretical support for the genetic improvement of maize root systems, and it identified candidate genes that may act as genetic resources for breeding.

Gene Co-expression Network Analysis

Gene co-expression analysis is a data analysis technique that helps identify groups of genes with similar expression patterns across several different conditions. By means of these techniques, different groups have been able to assign putative metabolic pathways and functions to understudied genes and to identify novel metabolic regulation networks for different metabolites. Some groups have even used network comparative studies to understand the evolution of these networks from green algae to land plants. In this chapter, we will go over the basic definitions required to understand network topology and gene module identification. Additionally, we offer the reader a walk-through a standard analysis pipeline as implemented in the package WGCNA that takes as input raw fastq files and obtains co-expressed gene clusters and representative gene expression patterns from each module for downstream applications.

Plant Microbiome-Based Genome-Wide Association Studies

Plants form intimate associations with microorganisms, and these associations are directly impacted by the host genotype. However, identifying specific host genetic pathways that influence these microbial interactions has proved challenging. Genome-wide association-based approaches that use features of microbiome composition as a quantitative trait represent a novel and underutilized strategy to identify such pathways. Several recent studies have demonstrated the potential utility of plant microbiome-based genome-wide association studies (GWAS). In this chapter, we describe the process of implementing GWAS using the plant microbiome as the primary quantitative trait, considering experimental design, sample harvest, and processing, but with an emphasis on data filtering, data normalization, and statistical analyses.

Joint GWAS and WGCNA uncover the genetic control of calcium accumulation under salt treatment in maize seedlings

Soil salinization is an important factor threatening the yield and quality of maize. Ca²⁺ plays a considerable role in regulating plant growth under salt stress. Herein, we examined the shoot Ca²⁺ concentrations, root Ca²⁺ concentrations, and transport coefficients of seedlings in an association panel composed of 305 maize inbred lines under normal and salt conditions. A genome-wide association study was conducted by using the investigated phenotypes and 46,408 single-nucleotide polymorphisms of the panel. As a result, 53 significant SNPs were specifically detected under salt treatment, and 544 genes were identified in the linkage disequilibrium regions of these SNPs. According to the expression data of the 544 genes, we carried out a weighted coexpression network analysis. Combining the enrichment analyses and functional annotations, four hub genes (GRMZM2G051032, GRMZM2G004314, GRMZM2G421669, and GRMZM2G123314) were finally determined, which were then used to evaluate the genetic variation effects by gene-based association analysis. Only GRMZM2G123314, which encodes a pentatricopeptide repeat protein, was significantly associated with Ca²⁺ transport and the haplotype G-CT was identified as the superior haplotype. Our study brings novel insights into the genetic and molecular mechanisms of salt stress response and contributes to the development of salt-tolerant varieties in maize.

Integration of genome-wide association studies and gene coexpression networks unveils promising soybean resistance genes against five common fungal pathogens

Soybean is one of the most important legume crops worldwide. However, soybean yield is dramatically affected by fungal diseases, leading to economic losses of billions of dollars yearly. Here, we integrated publicly available genome-wide association studies and transcriptomic data to prioritize candidate genes associated with resistance to Cadophora gregata, Fusarium graminearum, Fusarium virguliforme, Macrophomina phaseolina, and Phakopsora pachyrhizi. We identified 188, 56, 11, 8, and 3 high-confidence candidates for resistance to F. virguliforme, F. graminearum, C. gregata, M. phaseolina and P. pachyrhizi, respectively. The prioritized candidate genes are highly conserved in the pangenome of cultivated soybeans and are heavily biased towards fungal species-specific defense responses. The vast majority of the prioritized candidate resistance genes are related to plant immunity processes, such as recognition, signaling, oxidative stress, systemic acquired resistance, and physical defense. Based on the number of resistance alleles, we selected the five most resistant accessions against each fungal species in the soybean USDA germplasm. Interestingly, the most resistant accessions do not reach the maximum theoretical resistance potential. Hence, they can be further improved to increase resistance in breeding programs or through genetic engineering. Finally, the coexpression network generated here is available in a user-friendly web application ( https://soyfungigcn.venanciogroup.uenf.br/ ) and an R/Shiny package ( https://github.com/almeidasilvaf/SoyFungiGCN ) that serve as a public resource to explore soybean-pathogenic fungi interactions at the transcriptional level.

Unravelling Rubber Tree Growth by Integrating GWAS and Biological Network-Based Approaches

Hevea brasiliensis (rubber tree) is a large tree species of the Euphorbiaceae family with inestimable economic importance. Rubber tree breeding programs currently aim to improve growth and production, and the use of early genotype selection technologies can accelerate such processes, mainly with the incorporation of genomic tools, such as marker-assisted selection (MAS). However, few quantitative trait loci (QTLs) have been used successfully in MAS for complex characteristics. Recent research shows the efficiency of genome-wide association studies (GWAS) for locating QTL regions in different populations. In this way, the integration of GWAS, RNA-sequencing (RNA-Seq) methodologies, coexpression networks and enzyme networks can provide a better understanding of the molecular relationships involved in the definition of the phenotypes of interest, supplying research support for the development of appropriate genomic based strategies for breeding. In this context, this work presents the potential of using combined multiomics to decipher the mechanisms of genotype and phenotype associations involved in the growth of rubber trees. Using GWAS from a genotyping-by-sequencing (GBS) Hevea population, we were able to identify molecular markers in QTL regions with a main effect on rubber tree plant growth under constant water stress. The underlying genes were evaluated and incorporated into a gene coexpression network modelled with an assembled RNA-Seq-based transcriptome of the species, where novel gene relationships were estimated and evaluated through in silico methodologies, including an estimated enzymatic network. From all these analyses, we were able to estimate not only the main genes involved in defining the phenotype but also the interactions between a core of genes related to rubber tree growth at the transcriptional and translational levels. This work was the first to integrate multiomics analysis into the in-depth investigation of rubber tree plant growth, producing useful data for future genetic studies in the species and enhancing the efficiency of the species improvement programs.

Addressing noise in co-expression network construction

Gene co-expression networks (GCNs) provide multiple benefits to molecular research including hypothesis generation and biomarker discovery. Transcriptome profiles serve as input for GCN construction and are derived from increasingly larger studies with samples across multiple experimental conditions, treatments, time points, genotypes, etc. Such experiments with larger numbers of variables confound discovery of true network edges, exclude edges and inhibit discovery of context (or condition) specific network edges. To demonstrate this problem, a 475-sample dataset is used to show that up to 97% of GCN edges can be misleading because correlations are false or incorrect. False and incorrect correlations can occur when tests are applied without ensuring assumptions are met, and pairwise gene expression may not meet test assumptions if the expression of at least one gene in the pairwise comparison is a function of multiple confounding variables. The ‘one-size-fits-all’ approach to GCN construction is therefore problematic for large, multivariable datasets. Recently, the Knowledge Independent Network Construction toolkit has been used in multiple studies to provide a dynamic approach to GCN construction that ensures statistical tests meet assumptions and confounding variables are addressed. Additionally, it can associate experimental context for each edge of the network resulting in context-specific GCNs (csGCNs). To help researchers recognize such challenges in GCN construction, and the creation of csGCNs, we provide a review of the workflow.

Genetic architecture of the metabolic pathway of salicylic acid biosynthesis in Populus

Salicylic acid (SA) is a vital hormone for adaptive responses to biotic and abiotic stresses, which facilitates growth-immunity trade-offs in plants. However, the genetic regulatory networks underlying the metabolic pathway of SA biosynthesis in perennial species remain unclear. Here, we integrated genome-wide association study (GWAS) with metabolite and expression profiling methodologies to dissect the genetic architecture of SA biosynthesis in Populus. First, we quantified nine intermediate metabolites of SA biosynthesis in 300 unrelated Populus tomentosa Carr. individuals. Then, we used a systematic genetic strategy to identify candidate genes for constructing the genetic regulatory network of SA biosynthesis. We focused on WRKY70, an efficient transcription factor, as the key causal gene in the regulatory network, and combined the novel genes coordinating the accumulation of SA. Finally, we identified eight GWAS signals and eight expression quantitative trait loci situated in a selective sweep, and showed the presence of large allele frequency differences among the three geographic populations, revealing that candidate genes subject to selection were involved in SA biosynthesis. This study provides an integrated strategy for dissecting the genetic architecture of the metabolic pathway of SA biosynthesis in Populus, thereby enhancing our understanding of genetic regulation of SA biosynthesis in trees, and accelerating marker-assisted breeding efforts toward high-resistance elite varieties of Populus.