Joint genetic and network analyses identify loci associated with root growth under NaCl stress in Arabidopsis thaliana

Comprehensive Evaluation and Transcriptome Analysis Reveal the Salt Tolerance Mechanism in Semi-Wild Cotton (Gossypium purpurascens)

Elevated salinity significantly threatens cotton growth, particularly during the germination and seedling stages. The utilization of primitive species of Gossypium hirsutum, specifically Gossypium purpurascens, has the potential to facilitate the restoration of genetic diversity that has been depleted due to selective breeding in modern cultivars. This investigation evaluated 45 G. purpurascens varieties and a salt-tolerant cotton variety based on 34 morphological, physiological, and biochemical indicators and comprehensive salt tolerance index values. This study effectively identified a total of 19 salt-tolerant and two salt-resistant varieties. Furthermore, transcriptome sequencing of a salt-tolerant genotype (Nayanmian-2; NY2) and a salt-sensitive genotype (Sanshagaopao-2; GP2) revealed 2776, 6680, 4660, and 4174 differentially expressed genes (DEGs) under 0.5, 3, 12, and 24 h of salt stress. Gene ontology enrichment analysis indicated that the DEGs exhibited significant enrichment in biological processes like metabolic (GO:0008152) and cellular (GO:0009987) processes. MAPK signaling, plant-pathogen interaction, starch and sucrose metabolism, plant hormone signaling, photosynthesis, and fatty acid metabolism were identified as key KEGG pathways involved in salinity stress. Among the DEGs, including NAC, MYB, WRKY, ERF, bHLH, and bZIP, transcription factors, receptor-like kinases, and carbohydrate-active enzymes were crucial in salinity tolerance. Weighted gene co-expression network analysis (WGCNA) unveiled associations of salt-tolerant genotypes with flavonoid metabolism, carbon metabolism, and MAPK signaling pathways. Identifying nine hub genes (MYB4, MYB105, MYB36, bZIP19, bZIP43, FRS2 SMARCAL1, BBX21, F-box) across various intervals offered insights into the transcriptional regulation mechanism of salt tolerance in G. purpurascens. This study lays the groundwork for understanding the important pathways and gene networks in response to salt stress, thereby providing a foundation for enhancing salt tolerance in upland cotton.

Combining GS-assisted GWAS and transcriptome analysis to mine candidate genes for nitrogen utilization efficiency in Populus cathayana

BACKGROUND

Forest trees such as poplar, shrub willow, et al. are essential natural resources for sustainable and renewable energy production, and their wood can reduce dependence on fossil fuels and reduce environmental pollution. However, the productivity of forest trees is often limited by the availability of nitrogen (N), improving nitrogen use efficiency (NUE) is an important way to address it. Currently, NUE genetic resources are scarce in forest tree research, and more genetic resources are urgently needed.

RESULTS

Here, we performed genome-wide association studies (GWAS) using the mixed linear model (MLM) to identify genetic loci regulating growth traits in Populus cathayana at two N levels, and attempted to enhance the signal strength of single nucleotide polymorphism (SNP) detection by performing genome selection (GS) assistance GWAS. The results of the two GWAS analyses identified 55 and 40 SNPs that were respectively associated with plant height (PH) and ground diameter (GD), and 92 and 69 candidate genes, including 30 overlapping genes. The prediction accuracy of the GS model (rrBLUP) for phenotype exceeds 0.9. Transcriptome analysis of 13 genotypes under two N levels showed that genes related to carbon and N metabolism, amino acid metabolism, energy metabolism, and signal transduction were differentially expressed in the xylem of P. cathayana under N treatment. Furthermore, we observed strong regional patterns in gene expression levels of P. cathayana, with significant differences between different regions. Among them, P. cathayana in Longquan region exhibited the highest response to N. Finally, through weighted gene co-expression network analysis (WGCNA), we identified a module closely related to the N metabolic process and eight hub genes.

CONCLUSIONS

Integrating the GWAS, RNA-seq and WGCNA data, we ultimately identified four key regulatory genes (PtrNAC123, PtrNAC025, Potri.002G233100, and Potri.006G236200) involved in the wood formation process, and they may affect P. cathayana growth and wood formation by regulating nitrogen metabolism. This study will provide strong evidence for N regulation mechanisms, and reliable genetic resources for growth and NUE genetic improvement in poplar.

Weighted gene co-expression network analysis revealed the key pathways and hub genes of potassium regulating cotton root adaptation to salt stress

Soil salinization is one of the main abiotic stresses affecting cotton yield and planting area. Potassium application has been proven to be an important strategy to reduce salt damage in agricultural production. However, the mechanism of potassium regulating the salt adaptability of cotton has not been fully elucidated. In the present research, the appropriate potassium application rate for alleviating salt damage of cotton based on different K⁺/Na⁺ ratios we screened, and a gene co-expression network based on weighted gene co-expression network analysis (WGCNA) using the transcriptome data sets treated with CK (0 mM NaCl), S (150 mM NaCl), and SK8 (150 mM NaCl + 9.38 mM K₂SO₄) was constructed. In this study, four key modules that are highly related to potassium regulation of cotton salt tolerance were identified, and the mitogen-activated protein kinase (MAPK) signaling pathway, tricarboxylic acid (TCA) cycle and glutathione metabolism pathway were identified as the key biological processes and metabolic pathways for potassium to improve cotton root salt adaptability. In addition, 21 hub genes and 120 key candidate genes were identified in this study, suggesting that they may play an important role in the enhancement of salt adaptability of cotton by potassium. The key modules, key biological pathways and hub genes discovered in this study will provide a new understanding of the molecular mechanism of potassium enhancing salinity adaptability in cotton, and lay a theoretical foundation for the improvement and innovation of high-quality cotton germplasm.

WGCNA Analysis Revealed the Hub Genes Related to Soil Cadmium Stress in Maize Kernel (Zea mays L.)

Soil contamination by heavy metals has become a prevalent topic due to their widespread release from industry, agriculture, and other human activities. Great progress has been made in elucidating the uptake and translocation of cadmium (Cd) accumulation in rice. However, there is still little known about corresponding progress in maize. In the current study, we performed a comparative RNA-Seq-based approach to identify differentially expressed genes (DEGs) of maize immature kernel related to Cd stress. In total, 55, 92, 22, and 542 DEGs responsive to high cadmium concentration soil were identified between XNY22-CHS-8 vs. XNY22-YA-8, XNY22-CHS-24 vs. XNY22-YA-24, XNY27-CHS-8 vs. XNY27-YA-8, and XNY27-CHS-24 vs. XNY27-YA-24, respectively. The weighted gene co-expression network analysis (WGCNA) categorized the 9599 Cd stress-responsive hub genes into 37 different gene network modules. Combining the hub genes and DEGs, we obtained 71 candidate genes. Gene Ontology (GO) enrichment analysis of genes in the greenyellow module in XNY27-YA-24 and connectivity genes of these 71 candidate hub genes showed that the responses to metal ion, inorganic substance, abiotic stimulus, hydrogen peroxide, oxidative stress, stimulus, and other processes were enrichment. Moreover, five candidate genes that were responsive to Cd stress in maize kernel were detected. These results provided the putative key genes and pathways to response to Cd stress in maize kernel, and a useful dataset for unraveling the underlying mechanism of Cd accumulation in maize kernel.

How do plants remember drought?

Plants develop both short-term and transgenerational memory of drought stress through epigenetic regulation of transcription for a better response to subsequent exposure. Recurrent spells of droughts are more common than a single drought, with intermittent moist recovery intervals. While the detrimental effects of the first drought on plant structure and physiology are unavoidable, if survived, plants can memorize the first drought to present a more robust response to the following droughts. This includes a partial stomatal opening in the watered recovery interval, higher levels of osmoprotectants and ABA, and attenuation of photosynthesis in the subsequent exposure. Short-term drought memory is regulated by ABA and other phytohormone signaling with transcriptional memory behavior in various genes. High levels of methylated histones are deposited at the drought-tolerance genes. During the recovery interval, the RNA polymerase is stalled to be activated by a pause-breaking factor in the subsequent drought. Drought leads to DNA demethylation near drought-response genes, with genetic control of the process. Progenies of the drought-exposed plants can better adapt to drought owing to the inheritance of particular methylation patterns. However, a prolonged watered recovery interval leads to loss of drought memory, mediated by certain demethylases and chromatin accessibility factors. Small RNAs act as critical regulators of drought memory by altering transcript levels of drought-responsive target genes. Further studies in the future will throw more light on the genetic control of drought memory and the interplay of genetic and epigenetic factors in its inheritance. Plants from extreme environments can give queues to understanding robust memory responses at the ecosystem level.

Genome-Wide Association Study of Salt Tolerance at the Seed Germination Stage in Flax (Linum usitatissimum L.)

Soil salinization seriously affects the growth and distribution of flax. However, there is little information about the salt tolerance of flax. In this study, the salt tolerance of 200 diverse flax accessions during the germination stage was evaluated, and then the Genome-wide Association Study (GWAS) was carried out based on the relative germination rate (RGR), relative shoot length (RSL) and relative root length (RRL), whereby quantitative trait loci (QTLs) related to salt tolerance were identified. The results showed that oil flax had a better salt tolerance than fiber flax. A total of 902 single nucleotide polymorphisms (SNPs) were identified on 15 chromosomes. These SNPs were integrated into 64 QTLs, explaining 14.48 to 29.38% (R²) of the phenotypic variation. In addition, 268 candidate genes were screened by combining previous transcriptome data and homologous gene annotation. Among them, Lus10033213 is a single-point SNP repeat mapping gene, which encodes a Glutathione S-transferase (GST). This study is the first to use GWAS to excavate genes related to salt tolerance during the germination stage of flax. The results of this study provide important information for studying the genetic mechanism of salt tolerance of flax, and also provide the possibility to improve the salt tolerance of flax.

Homologous Cloning of Potassium Channel Genes From the Superior Apple Rootstock Line 12-2, Which is Tolerant to Apple Replant Disease

Potassium channels are important ion channels that are responsible for the absorption of potassium in the plant nutrient uptake system. In this study, we used homologous molecular cloning to obtain 8 K⁺ channel genes from the superior apple rootstock line 12-2 (self-named): MsAKT1-1, MsKAT3-2, MsKAT1-3, MsK2P3-4, MsK2P3-5, MsK2P5-6, MsK2P3-7, and MsK2P3-8. Their lengths varied from 942 bp (MsK2P5-6) to 2625 bp (MsAKT1-1), and the number of encoded amino acids varied from 314 (MsK2P5-6) to 874 (MsAKT1-1). Subcellular localization predictions showed that MsAKT1-1, MsKAT3-2, and MsKAT1-3 were localized on the plasma membrane, and MsK2P3-4, MsK2P3-5, MsK2P5-6, MsK2P3-7, and MsK2P3-8 were localized on the vacuole and plasma membrane. The 8 K⁺ channel proteins contained α helices, extended strands, β turns, and random coils. MsKAT1-3 had four transmembrane structures, MsKAT3-2 had six, and the other six K⁺ channel genes had five. Protein structure domain analysis showed that MsAKT1-1 contained nine protein domains, followed by MsKAT3-2 with four, MsKAT1-3 with three, and the other five two-pore domain K⁺ channel proteins with two. Semi-quantitative RT-PCR detection of the K⁺ channel genes showed that their expression levels were high in roots. qRT-PCR analysis showed that the relative expression levels of the 8 genes changed after exposure to ARD stress. The above results provide a theoretical basis for further research on the functions of potassium channel genes in 12-2 and a scientific basis for the breeding of ARD-resistant rootstock.

Expression GWAS of PGIP1 Identifies STOP1-Dependent and STOP1-Independent Regulation of PGIP1 in Aluminum Stress Signaling in Arabidopsis

To elucidate the unknown regulatory mechanisms involved in aluminum (Al)-induced expression of POLYGALACTURONASE-INHIBITING PROTEIN 1 (PGIP1), which is one of the downstream genes of SENSITIVE TO PROTON RHIZOTOXICITY 1 (STOP1) regulating Al-tolerance genes, we conducted a genome-wide association analysis of gene expression levels (eGWAS) of PGIP1 in the shoots under Al stress using 83 Arabidopsis thaliana accessions. The eGWAS, conducted through a mixed linear model, revealed 17 suggestive SNPs across the genome having the association with the expression level variation in PGIP1. The GWAS-detected SNPs were directly located inside transcription factors and other genes involved in stress signaling, which were expressed in response to Al. These candidate genes carried different expression level and amino acid polymorphisms. Among them, three genes encoding NAC domain-containing protein 27 (NAC027), TRX superfamily protein, and R-R-type MYB protein were associated with the suppression of PGIP1 expression in their mutants, and accordingly, the system affected Al tolerance. We also found the involvement of Al-induced endogenous nitric oxide (NO) signaling, which induces NAC027 and R-R-type MYB genes to regulate PGIP1 expression. In this study, we provide genetic evidence that STOP1-independent NO signaling pathway and STOP1-dependent regulation in phosphoinositide (PI) signaling pathway are involved in the regulation of PGIP1 expression under Al stress.

Identification of Key Transcription Factors Related to Bacterial Spot Resistance in Pepper through Regulatory Network Analyses

Bacterial spot (BS), caused by Xanthomonas campestris pv. Vesicatoria (Xcv), severely affects the quality and yield of pepper. Thus, breeding new pepper cultivars with enhanced resistance to BS can improve economic benefits for pepper production. Identification of BS resistance genes is an essential step to achieve this goal. However, very few BS resistance genes have been well characterized in pepper so far. In this study, we reanalyzed public multiple time points related to RNA-seq data sets from two pepper cultivars, the Xcv-susceptible cultivar ECW and the Xcv-resistant cultivar VI037601, post Xcv infection. We identified a total of 3568 differentially expressed genes (DEGs) between two cultivars post Xcv infection, which were mainly involved in some biological processes, such as Gene Ontology (GO) terms related to defense response to bacterium, immune system process, and regulation of defense response, etc. Through weighted gene co-expression network analysis (WGCNA), we identified 15 hub (Hub) transcription factor (TF) candidates in response to Xcv infection. We further selected 20 TFs from the gene regulatory network (GRN) potentially involved in Xcv resistance response. Finally, we predicted 4 TFs, C3H (p-coumarate 3-hydroxylase), ERF (ethylene-responsive element binding factor), TALE (three-amino-acid-loop-extension), and HSF (heat shock transcription factor), as key factors responsible for BS disease resistance in pepper. In conclusion, our study provides valuable resources for dissecting the underlying molecular mechanism responsible for Xcv resistance in pepper. Additionally, it also provides valuable references for mining transcriptomic data to identify key candidates for disease resistance in horticulture crops.

Genetic mapping of the early responses to salt stress in Arabidopsis thaliana

Salt stress decreases plant growth prior to significant ion accumulation in the shoot. However, the processes underlying this rapid reduction in growth are still unknown. To understand the changes in salt stress responses through time and at multiple physiological levels, examining different plant processes within a single set-up is required. Recent advances in phenotyping has allowed the image-based estimation of plant growth, morphology, colour and photosynthetic activity. In this study, we examined the salt stress-induced responses of 191 Arabidopsis accessions from 1 h to 7 days after treatment using high-throughput phenotyping. Multivariate analyses and machine learning algorithms identified that quantum yield measured in the light-adapted state (F_v' /F_m' ) greatly affected growth maintenance in the early phase of salt stress, whereas the maximum quantum yield (QY_max ) was crucial at a later stage. In addition, our genome-wide association study (GWAS) identified 770 loci that were specific to salt stress, in which two loci associated with QY_max and F_v' /F_m' were selected for validation using T-DNA insertion lines. We characterized an unknown protein kinase found in the QY_max locus that reduced photosynthetic efficiency and growth maintenance under salt stress. Understanding the molecular context of the candidate genes identified will provide valuable insights into the early plant responses to salt stress. Furthermore, our work incorporates high-throughput phenotyping, multivariate analyses and GWAS, uncovering details of temporal stress responses and identifying associations across different traits and time points, which are likely to constitute the genetic components of salinity tolerance.

A novel motif in the 5'-UTR of an orphan gene 'Big Root Biomass' modulates root biomass in sesame

Developing crops with improved root system is crucial in current global warming scenario. Underexploited crops are valuable reservoirs of unique genes that can be harnessed for the improvement of major crops. In this study, we performed genome-wide association studies on seven root traits in sesame (Sesamum indicum L.) and uncovered 409 significant signals, 19 quantitative trait loci containing 32 candidate genes. A peak SNP significantly associated with root number and root dry weight traits was located in the promoter of the gene named 'Big Root Biomass' (BRB), which was subsequently validated in a bi-parental population. BRB has no functional annotation and is restricted to the Lamiales order. We detected the presence of a novel motif 'AACACACAC' located in the 5'-UTR of BRB in single and duplicated copy in accessions with high and small root biomass, respectively. A strong expression level of BRB was negatively correlated with high root biomass, and this was attributed to the gene SiMYB181 which represses the activity of BRB by binding specifically to the single motif but not to the duplicated one. Curiously, the allele that enhanced BRB expression has been intensively selected by modern breeding. Overexpression of BRB in Arabidopsis modulates auxin pathway leading to reduced root biomass, improved yield parameters under normal growth conditions and increased drought stress sensitivity. Overall, BRB represents a solid gene model for improving the performance of sesame and other crops.

Functions and structure of roots and their contributions to salinity tolerance in plants

Soil salinity is an increasing threat to the productivity of glycophytic crops worldwide. The root plays vital roles under various stress conditions, including salinity, as well as has diverse functions in non-stress soil environments. In this review, we focus on the essential functions of roots such as in ion homeostasis mediated by several different membrane transporters and signaling molecules under salinity stress and describe recent advances in the impacts of quantitative trait loci (QTLs) or genetic loci (and their causal genes, if applicable) on salinity tolerance. Furthermore, we introduce important literature for the development of barriers against the apoplastic flow of ions, including Na⁺, as well as for understanding the functions and components of the barrier structure under salinity stress.

Expression genome-wide association study identifies that phosphatidylinositol-derived signalling regulates ALUMINIUM SENSITIVE3 expression under aluminium stress in the shoots of Arabidopsis thaliana

To identify unknown regulatory mechanisms leading to aluminium (Al)-induction of the Al tolerance gene ALS3, we conducted an expression genome-wide association study (eGWAS) for ALS3 in the shoots of 95 Arabidopsis thaliana accessions in the presence of Al. The eGWAS was conducted using a mixed linear model with 145,940 genome-wide single nucleotide polymorphisms (SNPs) and the association results were validated using reverse genetics. We found that many SNPs from the eGWAS were associated with genes related to phosphatidylinositol metabolism as well as stress signal transduction, including Ca²⁺signals, inter-connected in a co-expression network. Of these, PLC9, CDPK32, ANAC071, DIR1, and a hypothetical protein (AT4G10470) possessed amino acid sequence/ gene expression level polymorphisms that were significantly associated with ALS3 expression level variation. Furthermore, T-DNA insertion mutants of PLC9, CDPK32, and ANAC071 suppressed shoot ALS3 expression in the presence of Al. This study clarified the regulatory mechanisms of ALS3 expression in the shoot and provided genetic evidence of the involvement of phosphatidylinositol-derived signal transduction under Al stress.

Identification of Genes Conferring Plant Salt Tolerance using GWAS: Current Success and Perspectives

An understanding of the molecular mechanisms that underlie plant salt tolerance is important for both economic and scientific interests. Genome-wide association study (GWAS) is a promising approach to pinpoint genes that confer plant salt tolerance. With the advancement of supporting technology and methodology, GWAS has enabled the discovery of genes that play central roles in regulating plant salt tolerance in the past decade. Here, I highlight recent successful GWAS work in unveiling the molecular factors underlying plant salt tolerance and discuss the concerns and opportunities in conducting such experiments. It is anticipated that GWAS will be increasingly successful in the identification of key genes that are useful for crop improvement.

Gene Regulatory Network Inference: Connecting Plant Biology and Mathematical Modeling

Plant responses to environmental and intrinsic signals are tightly controlled by multiple transcription factors (TFs). These TFs and their regulatory connections form gene regulatory networks (GRNs), which provide a blueprint of the transcriptional regulations underlying plant development and environmental responses. This review provides examples of experimental methodologies commonly used to identify regulatory interactions and generate GRNs. Additionally, this review describes network inference techniques that leverage gene expression data to predict regulatory interactions. These computational and experimental methodologies yield complex networks that can identify new regulatory interactions, driving novel hypotheses. Biological properties that contribute to the complexity of GRNs are also described in this review. These include network topology, network size, transient binding of TFs to DNA, and competition between multiple upstream regulators. Finally, this review highlights the potential of machine learning approaches to leverage gene expression data to predict phenotypic outputs.

Exploiting natural variation in root system architecture via genome-wide association studies

Root growth and development has become an important research topic for breeders and researchers based on a growing need to adapt plants to changing and more demanding environmental conditions worldwide. Over the last few years, genome-wide association studies (GWASs) became an important tool to identify the link between traits in the field and their genetic background. Here we give an overview of the current literature concerning GWASs performed on root system architecture (RSA) in plants. We summarize which root traits and approaches have been used for GWAS, mentioning their respective success rate towards a successful gene discovery. Furthermore, we zoom in on the current technical hurdles in root phenotyping and GWAS, and discuss future possibilities in this field of research.

Genome-Wide Association Study and Genomic Prediction Elucidate the Distinct Genetic Architecture of Aluminum and Proton Tolerance in Arabidopsis thaliana

Under acid soil conditions, Al stress and proton stress can occur, reducing root growth and function. However, these stressors are distinct, and tolerance to each is governed by multiple physiological processes. To better understand the genes that underlie these coincidental but experimentally separable stresses, a genome-wide association study (GWAS) and genomic prediction (GP) models were created for approximately 200 diverse Arabidopsis thaliana accessions. GWAS and genomic prediction identified 140/160 SNPs associated with Al and proton tolerance, respectively, which explained approximately 70% of the variance observed. Reverse genetics of the genes in loci identified novel Al and proton tolerance genes, including TON1-RECRUITING MOTIF 28 (AtTRM28) and THIOREDOXIN H-TYPE 1 (AtTRX1), as well as genes known to be associated with tolerance, such as the Al-activated malate transporter, AtALMT1. Additionally, variation in Al tolerance was partially explained by expression level polymorphisms of AtALMT1 and AtTRX1 caused by cis-regulatory allelic variation. These results suggest that we successfully identified the loci that regulate Al and proton tolerance. Furthermore, very small numbers of loci were shared by Al and proton tolerance as determined by the GWAS. There were substantial differences between the phenotype predicted by genomic prediction and the observed phenotype for Al tolerance. This suggested that the GWAS-undetectable genetic factors (e.g., rare-allele mutations) contributing to the variation of tolerance were more important for Al tolerance than for proton tolerance. This study provides important new insights into the genetic architecture that produces variation in the tolerance of acid soil.

Weighted Gene Co-expression Network Analysis (WGCNA) Reveals the Hub Role of Protein Ubiquitination in the Acquisition of Desiccation Tolerance in Boea hygrometrica

Boea hygrometrica can survive extreme drought conditions and has been used as a model to study desiccation tolerance. A genome-wide transcriptome analysis of B. hygrometrica showed that the plant can survive rapid air-drying after experiencing a slow soil-drying acclimation phase. In addition, a weighted gene co-expression network analysis was used to study the transcriptomic datasets. A network comprising 22 modules was constructed, and seven modules were found to be significantly related to desiccation response using an enrichment analysis. Protein ubiquitination was observed to be a common process linked to hub genes in all the seven modules. Ubiquitin-modified proteins with diversified functions were identified using immunoprecipitation coupled with mass spectrometry. The lowest level of ubiquitination was noted at the full soil drying priming stage, which coincided the accumulation of dehydration-responsive gene BhLEA2. The highly conserved RY motif (CATGCA) was identified from the promoters of ubiquitin-related genes that were downregulated in the desiccated samples. An in silico gene expression analysis showed that the negative regulation of ubiquitin-related genes is potentially mediated via a B3 domain-containing transcription repressor VAL1. This study suggests that priming may involve the transcriptional regulation of several major processes, and the transcriptional regulation of genes in protein ubiquitination may play a hub role to deliver acclimation signals to posttranslational level in the acquisition of desiccation tolerance in B. hygrometrica.

WGCNA Analysis of Salt-Responsive Core Transcriptome Identifies Novel Hub Genes in Rice

Rice, being a major staple food crop and sensitive to salinity conditions, bears heavy yield losses due to saline soil. Although some salt responsive genes have been identified in rice, their applications in developing salt tolerant cultivars have resulted in limited achievements. Herein, we used bioinformatic approaches to perform a meta-analysis of three transcriptome datasets from salinity and control conditions in order to reveal novel genes and the molecular pathways underlying rice response to salt. From a total of 28,432 expressed genes, we identify 457 core differentially expressed genes (DEGs) constitutively responding to salt, regardless of the stress duration, genotype, or the tissue. Gene co-expression analysis divided the core DEGs into three different modules, each of them contributing to salt response in a unique metabolic pathway. Gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses highlighted key biological processes and metabolic pathways involved in the salt response. We identified important novel hub genes encoding proteins of different families including CAM, DUF630/632, DUF581, CHL27, PP2-13, LEA4-5, and transcription factors, which could be functionally characterized using reverse genetic experiments. This novel repertoire of candidate genes related to salt response in rice will be useful for engineering salt tolerant varieties.

Identification of Loci and Candidate Genes Responsible for Pod Dehiscence in Soybean via Genome-Wide Association Analysis Across Multiple Environments

Pod dehiscence (shattering) is the main cause of serious yield loss during the soybean mechanical harvesting process. A better understanding of the genetic architecture and molecular mechanisms of pod dehiscence is of great significance for soybean breeding. In this study, genome-wide association analysis (GWAS) with NJAU 355K SoySNP array was performed to detect single nucleotide polymorphisms (SNPs) associated with pod dehiscence in an association panel containing 211 accessions across five environments. A total of 163 SNPs were identified as significantly associated with pod dehiscence. Among these markers, 136 SNPs identified on chromosome 16 were located in the known QTL qPDH1. One, one, three, eleven, three, one, three, three and one SNPs were distributed on chromosome 1, 4, 6, 8, 9, 11, 17, 18, and 20, respectively. Favorable SNPs and six haplotypes were identified based on ten functional SNPs; among those Hap2 and Hap3 were considered as optimal haplotypes. In addition, based on GWAS results, the candidate gene Glyma09g06290 was identified. Quantitative real-time PCR (qRT-PCR) results and polymorphism analysis suggested that Glyma09g06290 might be involved in pod dehiscence. Furthermore, a derived cleaved amplified polymorphic sequences (dCAPS) marker for Glyma09g06290 was developed. Overall, the loci and genes identified in this study will be helpful in breeding soybean accessions resistant to pod dehiscence.

Genome-wide association analysis of salinity responsive traits in Medicago truncatula

Salinity stress is an important cause of crop yield loss in many parts of the world. Here, we performed genome-wide association studies of salinity-stress responsive traits in 132 HapMap genotypes of the model legume Medicago truncatula. Plants grown in soil were subjected to a step-wise increase in NaCl concentration, from 0 through 0.5% and 1.0% to 1.5%, and the following traits were measured: vigor, shoot biomass, shoot water content, leaf chlorophyll content, leaf size, and leaf and root concentrations of proline and major ions (Na+ , Cl- , K+ , Ca2+ , etc.). Genome-wide association studies were carried out using 2.5 million single nucleotide polymorphisms, and 12 genomic regions associated with at least four traits each were identified. Transcript-level analysis of the top eight candidate genes in five extreme genotypes revealed association between salinity tolerance and transcript-level changes for seven of the genes, encoding a vacuolar H+ -ATPase, two transcription factors, two proteins involved in vesicle trafficking, one peroxidase, and a protein of unknown function. Earlier functional studies on putative orthologues of two of the top eight genes (a vacuolar H+ -ATPase and a peroxidase) demonstrated their involvement in plant salinity tolerance.

Large Scale Proteomic Data and Network-Based Systems Biology Approaches to Explore the Plant World

The investigation of plant organisms by means of data-derived systems biology approaches based on network modeling is mainly characterized by genomic data, while the potential of proteomics is largely unexplored. This delay is mainly caused by the paucity of plant genomic/proteomic sequences and annotations which are fundamental to perform mass-spectrometry (MS) data interpretation. However, Next Generation Sequencing (NGS) techniques are contributing to filling this gap and an increasing number of studies are focusing on plant proteome profiling and protein-protein interactions (PPIs) identification. Interesting results were obtained by evaluating the topology of PPI networks in the context of organ-associated biological processes as well as plant-pathogen relationships. These examples foreshadow well the benefits that these approaches may provide to plant research. Thus, in addition to providing an overview of the main-omic technologies recently used on plant organisms, we will focus on studies that rely on concepts of module, hub and shortest path, and how they can contribute to the plant discovery processes. In this scenario, we will also consider gene co-expression networks, and some examples of integration with metabolomic data and genome-wide association studies (GWAS) to select candidate genes will be mentioned.

Natural allelic variation of the AZI1 gene controls root growth under zinc-limiting condition

Zinc is an essential micronutrient for all living organisms and is involved in a plethora of processes including growth and development, and immunity. However, it is unknown if there is a common genetic and molecular basis underlying multiple facets of zinc function. Here we used natural variation in Arabidopsis thaliana to study the role of zinc in regulating growth. We identify allelic variation of the systemic immunity gene AZI1 as a key for determining root growth responses to low zinc conditions. We further demonstrate that this gene is important for modulating primary root length depending on the zinc and defence status. Finally, we show that the interaction of the immunity signal azelaic acid and zinc level to regulate root growth is conserved in rice. This work demonstrates that there is a common genetic and molecular basis for multiple zinc dependent processes and that nutrient cues can determine the balance of growth and immune responses in plants.

An Integration of Genome-Wide Association Study and Gene Co-expression Network Analysis Identifies Candidate Genes of Stem Lodging-Related Traits in Brassica napus

Lodging is a persistent problem which severely reduce yield and impair seed quality in rapeseed (Brassica napus L.). Enhancing stem strength (SS) has proven to be an effective approach to decrease lodging risk. In the present study, four interrelated stem lodging-related traits, including stem breaking resistance (SBR), stem diameter (SD), SS, and lodging coefficient (LC), were investigated among 472 rapeseed accessions. A genome-wide association study (GWAS) using Brassica 60K SNP array for stem lodging-related traits identified 67 significantly associated quantitative trait loci (QTLs) and 71 candidate genes. In parallel, a gene co-expression network based on transcriptome sequencing was constructed. The module associated with cellulose biosynthesis was highlighted. By integrating GWAS and gene co-expression network analysis, some promising candidate genes, such as ESKIMO1 (ESK1, BnaC08g26920D), CELLULOSE SYNTHASE 6 (CESA6, BnaA09g06990D), and FRAGILE FIBER 8 (FRA8, BnaC04g39510D), were prioritized for further research. These findings revealed the genetic basis underlying stem lodging and provided worthwhile QTLs and genes information for genetic improvement of stem lodging resistance in B. napus.

Co-expression network analysis of the transcriptomes of rice roots exposed to various cadmium stresses reveals universal cadmium-responsive genes

BACKGROUND

The migration of cadmium (Cd) from contaminated soil to rice is a cause for concern. However, the molecular mechanism underlying the response of rice roots to various Cd stresses remains to be clarified from the viewpoint of the co-expression network at a system-wide scale.

RESULTS

We employed a comparative RNAseq-based approach to identify early Cd-responsive differentially expressed genes (DEGs) in rice 'Nipponbare' seedling roots after 1 h of high-Cd treatment. A multiplicity of the identified 1772 DEGs were implicated in hormone signaling and transcriptional regulation, particularly NACs and WRKYs were all upregulated under Cd stress. All of the 6 Cd-upregulated ABC transporters were pleiotropic drug resistance proteins (PDRs), whereas all of the 6 ZRT/IRT-like proteins (ZIPs) were consistently downregulated by Cd treatment. To further confirm our results of this early transcriptomic response to Cd exposure, we then conducted weighted gene co-expression network analysis (WGCNA) to re-analyze our RNAseq data in combination with other 11 previously published RNAseq datasets for rice roots exposed to diverse concentrations of Cd for extended treatment periods. This integrative approach identified 271 transcripts as universal Cd-regulated DEGs that are key components of the Cd treatment coupled co-expression module. A global view of the 164 transcripts with annotated functions in pathway networks revealed several Cd-upregulated key functional genes, including transporter ABCG36/OsPDR9, 12-oxo-phytodienoic acid reductases (OPRs) for JA synthesis, and ZIM domain proteins JAZs in JA signaling, as well as OsWRKY10, NAC, and ZFP transcription factors. More importantly, 104 of these, including ABCG36/OsPDR9, OsNAC3, as well as several orthologs in group metalloendoproteinase, plastocyanin-like domain containing proteins and pectin methylesterase inhibitor, may respond specifically to various Cd pressures, after subtracting the 60 general stress-responsive genes reported to be commonly upregulated following multiple stresses.

CONCLUSION

An integrative approach was implemented to identify DEGs and co-expression network modules in response to various Cd pressures, and 104 of the 164 annotatable universal Cd-responsive DEGs may specifically respond to various Cd pressures. These results provide insight into the universal molecular mechanisms beneath the Cd response in rice roots, and suggest many promising targets for improving the rice acclimation process against Cd toxicity.

Genome-wide Association Study Reveals that the Aquaporin NIP1;1 Contributes to Variation in Hydrogen Peroxide Sensitivity in Arabidopsis thaliana

Hydrogen peroxide (H₂O₂) is a reactive oxygen species that affects cell signaling in various plant defense responses and induces programmed cell death. To identify the new components associated with H₂O₂ signaling and tolerance, we conducted a genome-wide association study (GWAS) on the root growth of 133 Arabidopsis thaliana accessions grown in the presence of toxic H₂O₂ levels. The most significant SNPs were associated with a cluster of chromosome 4 genes encoding an aquaporin NODULIN 26-LIKE INTRINSIC PROTEIN 1; 1 (NIP1;1), an NB-ARC domain-containing disease resistance protein (AT4G19050), and a putative membrane lipoprotein (AT4G19070). The expression level of NIP1;1 was relatively high in A. thaliana accessions sensitive to H₂O₂. Additionally, overexpression of NIP1;1 in a tolerant accession (e.g., Col-0) increased the sensitivity of transgenic plants to H₂O₂. An in planta β-glucuronidase reporter assay revealed that variations in the NIP1;1 promoter were responsible for the differences of its expression level in H₂O₂-tolerant and -sensitive accessions. Cell death was extensive and H₂O₂ levels were high in the roots of H₂O₂-sensitive and NIP1;1-overexpressing accessions. Together, our results indicate that the aquaporin NIP1;1 is a key determinant of the sensitivity of A. thaliana to H₂O₂, and contributes to the phenotypic variations detected by our GWAS.

The common transcriptional subnetworks of the grape berry skin in the late stages of ripening

BACKGROUND

Wine grapes are important economically in many countries around the world. Defining the optimum time for grape harvest is a major challenge to the grower and winemaker. Berry skins are an important source of flavor, color and other quality traits in the ripening stage. Senescent-like processes such as chloroplast disorganization and cell death characterize the late ripening stage.

RESULTS

To better understand the molecular and physiological processes involved in the late stages of berry ripening, RNA-seq analysis of the skins of seven wine grape cultivars (Cabernet Franc, Cabernet Sauvignon, Merlot, Pinot Noir, Chardonnay, Sauvignon Blanc and Semillon) was performed. RNA-seq analysis identified approximately 2000 common differentially expressed genes for all seven cultivars across four different berry sugar levels (20 to 26 °Brix). Network analyses, both a posteriori (standard) and a priori (gene co-expression network analysis), were used to elucidate transcriptional subnetworks and hub genes associated with traits in the berry skins of the late stages of berry ripening. These independent approaches revealed genes involved in photosynthesis, catabolism, and nucleotide metabolism. The transcript abundance of most photosynthetic genes declined with increasing sugar levels in the berries. The transcript abundance of other processes increased such as nucleic acid metabolism, chromosome organization and lipid catabolism. Weighted gene co-expression network analysis (WGCNA) identified 64 gene modules that were organized into 12 subnetworks of three modules or more and six higher order gene subnetworks. Some gene subnetworks were highly correlated with sugar levels and some subnetworks were highly enriched in the chloroplast and nucleus. The petal R package was utilized independently to construct a true small-world and scale-free complex gene co-expression network model. A subnetwork of 216 genes with the highest connectivity was elucidated, consistent with the module results from WGCNA. Hub genes in these subnetworks were identified including numerous members of the core circadian clock, RNA splicing, proteolysis and chromosome organization. An integrated model was constructed linking light sensing with alternative splicing, chromosome remodeling and the circadian clock.

CONCLUSIONS

A common set of differentially expressed genes and gene subnetworks from seven different cultivars were examined in the skin of the late stages of grapevine berry ripening. A densely connected gene subnetwork was elucidated involving a complex interaction of berry senescent processes (autophagy), catabolism, the circadian clock, RNA splicing, proteolysis and epigenetic regulation. Hypotheses were induced from these data sets involving sugar accumulation, light, autophagy, epigenetic regulation, and fruit development. This work provides a better understanding of berry development and the transcriptional processes involved in the late stages of ripening.

Transcriptomic variation among six Arabidopsis thaliana accessions identified several novel genes controlling aluminium tolerance

Differences in the expression levels of aluminium (Al) tolerance genes are a known determinant of Al tolerance among plant varieties. We combined transcriptomic analysis of six Arabidopsis thaliana accessions with contrasting Al tolerance and a reverse genetic approach to identify Al-tolerance genes responsible for differences in Al tolerance between accession groups. Gene expression variation increased in the signal transduction process under Al stress and in growth-related processes in the absence of stress. Co-expression analysis and promoter single nucleotide polymorphism searching suggested that both trans-acting polymorphisms of Al signal transduction pathway and cis-acting polymorphisms in the promoter sequences caused the variations in gene expression associated with Al tolerance. Compared with the wild type, Al sensitivity increased in T-DNA knockout (KO) lines for five genes, including TARGET OF AVRB OPERATION1 (TAO1) and an unannotated gene (At5g22530). These were identified from 53 Al-inducible genes showing significantly higher expression in tolerant accessions than in sensitive accessions. These results indicate that the difference in transcriptional signalling is partly associated with the natural variation in Al tolerance in Arabidopsis. Our study also demonstrates the feasibility of comparative transcriptome analysis by using natural genetic variation for the identification of genes responsible for Al stress tolerance.