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

Transcriptomic and Physiological Studies Unveil that Brassinolide Maintains the Balance of Maize's Multiple Metabolisms under Low-Temperature Stress

Low-temperature (LT) is one of the major abiotic stresses that restrict the growth and development of maize seedlings. Brassinolides (BRs) have been shown to enhance LT tolerance in several plant species; the physiological and molecular mechanisms by which BRs enhance maize tolerance are still unclear. Here, we characterized changes in the physiology and transcriptome of N192 and Ji853 seedlings at the three-leaf stage with or without 2 μM 2,4-epibrassinolide (EBR) application at 25 and 15 °C environments via high-performance liquid chromatography and RNA-Sequencing. Physiological analyses revealed that EBR increased the antioxidant enzyme activities, enhanced the cell membrane stability, decreased the malondialdehyde formation, and inhibited the reactive oxygen species (ROS) accumulation in maize seedlings under 15 °C stress; meanwhile, EBR also maintained hormone balance by increasing indole-3-acetic acid and gibberellin 3 contents and decreasing the abscisic acid level under stress. Transcriptome analysis revealed 332 differentially expressed genes (DEGs) enriched in ROS homeostasis, plant hormone signal transduction, and the mitogen-activated protein kinase (MAPK) cascade. These DEGs exhibited synergistic and antagonistic interactions, forming a complex LT tolerance network in maize. Additionally, weighted gene co-expression network analysis (WGCNA) revealed that 109 hub genes involved in LT stress regulation pathways were discovered from the four modules with the highest correlation with target traits. In conclusion, our findings provide new insights into the molecular mechanisms of exogenous BRs in enhancing LT tolerance of maize at the seedling stage, thus opening up possibilities for a breeding program of maize tolerance to LT stress.

Joint GWAS and WGCNA Identify Genes Regulating the Isoflavone Content in Soybean Seeds

Isoflavone is a secondary metabolite of the soybean phenylpropyl biosynthesis pathway with physiological activity and is beneficial to human health. In this study, the isoflavone content of 205 soybean germplasm resources from 3 locations in 2020 showed wide phenotypic variation. A joint genome-wide association study (GWAS) and weighted gene coexpression network analysis (WGCNA) identified 33 single-nucleotide polymorphisms and 11 key genes associated with soybean isoflavone content. Gene ontology enrichment analysis, gene coexpression, and haplotype analysis revealed natural variations in the Glyma.12G109800 (GmOMT7) gene and promoter region, with Hap1 being the elite haplotype. Transient overexpression and knockout of GmOMT7 increased and decreased the isoflavone content, respectively, in hairy roots. The combination of GWAS and WGCNA effectively revealed the genetic basis of soybean isoflavone and identified potential genes affecting isoflavone synthesis and accumulation in soybean, providing a valuable basis for the functional study of soybean isoflavone.

Comprehensive Transcriptome and Proteome Analyses Reveal the Drought Responsive Gene Network in Potato Roots

The root system plays a decisive role in the growth and development of plants. The water requirement of a root system depends strongly on the plant species. Potatoes are an important food and vegetable crop grown worldwide, especially under irrigation in arid and semi-arid regions. However, the expected impact of global warming on potato yields calls for an investigation of genes related to root development and drought resistance signaling pathways in potatoes. In this study, we investigated the molecular mechanisms of different drought-tolerant potato root systems in response to drought stress under controlled water conditions, using potato as a model. We analyzed the transcriptome and proteome of the drought-sensitive potato cultivar Atlantic (Atl) and the drought-tolerant cultivar Qingshu 9 (Q9) under normal irrigation (CK) and weekly drought stress (D). The results showed that a total of 14,113 differentially expressed genes (DEGs) and 5596 differentially expressed proteins (DEPs) were identified in the cultivars. A heat map analysis of DEGs and DEPs showed that the same genes and proteins in Atl and Q9 exhibited different expression patterns under drought stress. Weighted gene correlation network analysis (WGCNA) showed that in Atl, Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG)-enriched pathways were related to pyruvate metabolism and glycolysis, as well as cellular signaling and ion transmembrane transporter protein activity. However, GO terms and KEGG-enriched pathways related to phytohormone signaling and the tricarboxylic acid cycle were predominantly enriched in Q9. The present study provides a unique genetic resource to effectively explore the functional genes and uncover the molecular regulatory mechanism of the potato root system in response to drought stress.

GWAS and WGCNA Analysis Uncover Candidate Genes Associated with Oil Content in Soybean

Soybean vegetable oil is an important source of the human diet. However, the analysis of the genetic mechanism leading to changes in soybean oil content is still incomplete. In this study, a total of 227 soybean materials were applied and analyzed by a genome-wide association study (GWAS). There are 44 quantitative trait nucleotides (QTNs) that were identified as associated with oil content. A total of six, four, and 34 significant QTN loci were identified in Xiangyang, Hulan, and Acheng, respectively. Of those, 26 QTNs overlapped with or were near the known oil content quantitative trait locus (QTL), and 18 new QTNs related to oil content were identified. A total of 594 genes were located near the peak single nucleotide polymorphism (SNP) from three tested environments. These candidate genes exhibited significant enrichment in tropane, piperidine, and pyridine alkaloid biosynthesiss (ko00960), ABC transporters (ko02010), photosynthesis-antenna proteins (ko00196), and betalain biosynthesis (ko00965). Combined with the GWAS and weighted gene co-expression network analysis (WGCNA), four candidate genes (Glyma.18G300100, Glyma.11G221100, Glyma.13G343300, and Glyma.02G166100) that may regulate oil content were identified. In addition, Glyma.18G300100 was divided into two main haplotypes in the studied accessions. The oil content of haplotype 1 is significantly lower than that of haplotype 2. Our research findings provide a theoretical basis for improving the regulatory mechanism of soybean oil content.

Impacts of salinity stress on crop plants: improving salt tolerance through genetic and molecular dissection

Improper use of water resources in irrigation that contain a significant amount of salts, faulty agronomic practices such as improper fertilization, climate change etc. are gradually increasing soil salinity of arable lands across the globe. It is one of the major abiotic factors that inhibits overall plant growth through ionic imbalance, osmotic stress, oxidative stress, and reduced nutrient uptake. Plants have evolved with several adaptation strategies at morphological and molecular levels to withstand salinity stress. Among various approaches, harnessing the crop genetic variability across different genepools and developing salinity tolerant crop plants offer the most sustainable way of salt stress mitigation. Some important major genetic determinants controlling salinity tolerance have been uncovered using classical genetic approaches. However, its complex inheritance pattern makes breeding for salinity tolerance challenging. Subsequently, advances in sequence based breeding approaches and functional genomics have greatly assisted in underpinning novel genetic variants controlling salinity tolerance in plants at the whole genome level. This current review aims to shed light on physiological, biochemical, and molecular responses under salt stress, defense mechanisms of plants, underlying genetics of salt tolerance through bi-parental QTL mapping and Genome Wide Association Studies, and implication of Genomic Selection to breed salt tolerant lines.

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.

Genome-Wide Meta-Analysis of QTLs Associated with Root Traits and Implications for Maize Breeding

Root system architecture (RSA), also known as root morphology, is critical in plant acquisition of soil resources, plant growth, and yield formation. Many QTLs associated with RSA or root traits in maize have been identified using several bi-parental populations, particularly in response to various environmental factors. In the present study, a meta-analysis of QTLs associated with root traits was performed in maize using 917 QTLs retrieved from 43 mapping studies published from 1998 to 2020. A total of 631 QTLs were projected onto a consensus map involving 19,714 markers, which led to the prediction of 68 meta-QTLs (MQTLs). Among these 68 MQTLs, 36 MQTLs were validated with the marker-trait associations available from previous genome-wide association studies for root traits. The use of comparative genomics approaches revealed several gene models conserved among the maize, sorghum, and rice genomes. Among the conserved genomic regions, the ortho-MQTL analysis uncovered 20 maize MQTLs syntenic to 27 rice MQTLs for root traits. Functional analysis of some high-confidence MQTL regions revealed 442 gene models, which were then subjected to in silico expression analysis, yielding 235 gene models with significant expression in various tissues. Furthermore, 16 known genes viz., DXS2, PHT, RTP1, TUA4, YUC3, YUC6, RTCS1, NSA1, EIN2, NHX1, CPPS4, BIGE1, RCP1, SKUS13, YUC5, and AW330564 associated with various root traits were present within or near the MQTL regions. These results could aid in QTL cloning and pyramiding in developing new maize varieties with specific root architecture for proper plant growth and development under optimum and abiotic stress conditions.

Upregulated ENC1 predicts unfavorable prognosis and correlates with immune infiltration in endometrial cancer

A better knowledge of the molecular process behind uterine corpus endometrial carcinoma (UCEC) is important for prognosis prediction and the development of innovative targeted gene therapies. The purpose of this research is to discover critical genes associated with UCEC. We analyzed the gene expression profiles of TCGA-UCEC and GSE17025, respectively, using Weighted Gene Co-expression Network Analysis (WGCNA) and differential gene expression analysis. From four sets of findings, a total of 95 overlapping genes were retrieved. On the 95 overlapping genes, KEGG pathway and GO enrichment analysis were conducted. Then, we mapped the PPI network of 95 overlapping genes using the STRING database. Twenty hub genes were evaluated using the Cytohubba plugin, including NR3C1, ATF3, KLF15, THRA, NR4A1, FOSB, PER3, HLF, NTRK3, EGR3, MAPK13, ARNTL2, PKM2, SCD, EIF5A, ADHFE1, RERGL, TUB, and ENC1. The expression levels of NR3C1, PKM2, and ENC1 were shown to be adversely linked with the survival time of UCEC patients using univariate Cox regression analysis and Kaplan-Meier survival calculation. ENC1 were also overexpressed in UCEC tumor tissues or cell lines, as shown by quantitative real-time PCR and Western blotting. Then we looked into it further and discovered that ENC1 expression was linked to tumor microenvironment and predicted various immunological checkpoints. In conclusion, our data indicate that ENC1 may be required for the development of UCEC and may serve as a future biomarker for diagnosis and therapy.

A Combination of a Genome-Wide Association Study and a Transcriptome Analysis Reveals circRNAs as New Regulators Involved in the Response to Salt Stress in Maize

Salinization seriously threatens the normal growth of maize, especially at the seedling stage. Recent studies have demonstrated that circular RNAs (circRNAs) play vital roles in the regulation of plant stress resistance. Here, we performed a genome-wide association study (GWAS) on the survival rate of 300 maize accessions under a salt stress treatment. A total of 5 trait-associated SNPs and 86 candidate genes were obtained by the GWAS. We performed RNA sequencing for 28 transcriptome libraries derived from 2 maize lines with contrasting salt tolerance under normal and salt treatment conditions. A total of 1217 highly expressed circRNAs were identified, of which 371 were responsive to a salt treatment. Using PCR and Sanger sequencing, we verified the reliability of these differentially expressed circRNAs. An integration of the GWAS and RNA-Seq analyses uncovered two differentially expressed hub genes (Zm00001eb013650 and Zm00001eb198930), which were regulated by four circRNAs. Based on these results, we constructed a regulation model of circRNA/miRNA/mRNA that mediated salt stress tolerance in maize. By conducting hub gene-based association analyses, we detected a favorable haplotype in Zm00001eb198930, which was responsible for high salt tolerance. These results help to clarify the regulatory relationship between circRNAs and their target genes as well as to develop salt-tolerant lines for maize breeding.

Genome-Wide Association Study Reveals the Genetic Basis of Kernel and Cob Moisture Changes in Maize at Physiological Maturity Stage

Low moisture content (MC) and high dehydration rate (DR) at physiological maturity affect grain mechanical harvest, transport, and storage. In this study, we used an association panel composed of 241 maize inbred lines to analyze ear moisture changes at physiological maturity stage. A genome-wide association study revealed nine significant SNPs and 91 candidate genes. One SNP (SYN38588) was repeatedly detected for two traits, and 15 candidate genes were scanned in the linkage disequilibrium regions of this SNP. Of these, genes Zm00001d020615 and Zm00001d020623 were individually annotated as a polygalacturonase (PG) and a copper transporter 5.1 (COPT5.1), respectively. Candidate gene association analysis showed that three SNPs located in the exons of Zm00001d020615 were significantly associated with the dehydration rate, and AATTAA was determined as the superior haplotype. All these findings suggested that Zm00001d020615 was a key gene affecting moisture changes of maize at the physiological maturity stage. These results have demonstrated the genetic basis of ear moisture changes in maize and indicated a superior haplotype for cultivating maize varieties with low moisture content and high dehydration rates.

Identification of major candidate genes for multiple abiotic stress tolerance at seedling stage by network analysis and their validation by expression profiling in rice (Oryza sativa L.)

A wealth of microarray and RNA-seq data for studying abiotic stress tolerance in rice exists but only limited studies have been carried out on multiple stress-tolerance responses and mechanisms. In this study, we identified 6657 abiotic stress-responsive genes pertaining to drought, salinity and heat stresses from the seedling stage microarray data of 83 samples and used them to perform unweighted network analysis and to identify key hub genes or master regulators for multiple abiotic stress tolerance. Of the total 55 modules identified from the analysis, the top 10 modules with 8-61 nodes comprised 239 genes. From these 10 modules, 10 genes common to all the three stresses were selected. Further, based on the centrality properties and highly dense interactions, we identified 7 intra-modular hub genes leading to a total of 17 potential candidate genes. Out of these 17 genes, 15 were validated by expression analysis using a panel of 4 test genotypes and a pair of standard check genotypes for each abiotic stress response. Interestingly, all the 15 genes showed upregulation under all stresses and in all the genotypes, suggesting that they could be representing some of the core abiotic stress-responsive genes. More pertinently, eight of the genes were found to be co-localized with the stress-tolerance QTL regions. Thus, in conclusion, our study not only provided an effective approach for studying abiotic stress tolerance in rice, but also identified major candidate genes which could be further validated by functional genomics for abiotic stress tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03182-7.