Genetic Conservation of CBS Domain Containing Protein Family in Oryza Species and Their Association with Abiotic Stress Responses

Tracing the intraspecies expansion of glyoxalase genes and their expanding roles across the genus Oryza

The genus Oryza is of utmost importance to human civilization as two of its species became agronomically productive and widely cultivated, and also because wild rice is a treasure trove of beneficial alleles that can be used for crop improvement. Most of the wild rice genotypes are known for their stress tolerance several times more than the domesticated rice varieties. In this study, we aimed to carry out an exhaustive genomic survey to identify glyoxalase I (GLYI) and glyoxalase II (GLYII) genes across the 11 rice genomes sequenced so far. Notably, we found the putatively functional metal-dependent GLYI and GLYII enzymes to be conserved throughout domestication and a few homologous pairs to have undergone beneficial mutations to drive positive selection, and thus, acquire newer functions. Interestingly, we also report four newly identified GLYII members in O. sativa subsp. japonica in addition to the three previously reported GLYII genes. The presence of different types of cis-elements in the promoter region of the glyoxalase genes gives insights into their role and regulation under various developmental processes besides stress adaptation. Publicly available data suggests the role of glyoxalase genes particularly in salinity stress in both wild and cultivated rice as is also confirmed through qRT-PCR. Interestingly, we found less accumulation of MG and concurrently higher enzymatic activity of GLYI and GLYII proteins in stressed seedlings of selected wild rice genotypes indicating that glyoxalases indeed contribute to the intrinsic stress tolerance of wild rice.

Comparative proteomic discovery of salt stress response in alfalfa roots and overexpression of MsANN2 confers salt tolerance

Soil salinity constrains growth, development and yield of alfalfa (Medicago sativa L.). To illustrate the molecular mechanisms responsible for salt tolerance, a comparative proteome analysis was explored to characterize protein profiles of alfalfa seedling roots exposed to 100 and 200 mM NaCl for three weeks. There were 52 differentially expressed proteins identified, among which the mRNA expressions of 12 were verified by Real-Time-PCR analysis. The results showed increase in abundance of ascorbate peroxidase, POD, CBS protein and PR-10 in salt-stressed alfalfa, suggesting an effectively antioxidant and defense systems. Alfalfa enhanced protein quality control system to refold or degrade abnormal proteins induced by salt stress through upregulation of unfolded protein response (UPR) marker PDIs and molecular chaperones (eg. HSP70, TCP-1, and GroES) as well as the ubiquitin-proteasome system (UPS) including ubiquitin ligase enzyme (E3) and proteasome subunits. Upregulation of proteins responsible for calcium signal transduction including calmodulin and annexin helped alfalfa adapt to salt stress. Specifically, annexin (MsANN2), a key Ca2+-binding protein, was selected for further characterization. The heterologous of the MsANN2 in Arabidopsis conferred salt tolerance. These results provide detailed information for salt-responsive root proteins and highlight the importance of MsANN2 in adapting to salt stress in alfalfa.

A cystathionine beta-synthase domain containing protein, OsCBSCBS4, interacts with OsSnRK1A and OsPKG and functions in abiotic stress tolerance in rice

The Cystathionine-β-Synthase (CBS) domain-containing proteins (CDCPs) constitute a functionally diverse protein superfamily, sharing an evolutionary conserved CBS domain either in pair or quad. Rice genome (Oryza sativa subsp. indica) encodes 42 CDCPs; their functions remain largely unexplored. This study examines OsCBSCBS4, a quadruple CBS domain containing protein towards its role in regulating the abiotic stress tolerance in rice. Gene expression analyses revealed upregulation of OsCBSCBS4 in response to diverse abiotic stresses. Further, the cytoplasm-localised OsCBSCBS4 showed interaction with two different kinases, a cytoplasmic localised cGMP-dependant protein kinase (OsPKG) and the nucleo-cytoplasmic catalytic subunit of sucrose-nonfermentation 1-related protein kinase 1 (OsSnRK1A). The interaction with the latter assisted in trafficking of OsCBSCBS4 to the nucleus as well. Overexpression of OsCBSCBS4 in rice resulted in enhanced tolerance to drought and salinity stress, via maintaining better physiological parameters and antioxidant activity. Additionally, OsCBSCBS4-overexpressing rice plants exhibited reduced yield penalty under stress conditions. The in silico docking and in vitro binding analyses of OsCBSCBS4 with ATP suggest its involvement in cellular energy balance. Overall, this study provides novel insight into the unexplored functions of OsCBSCBS4 and demonstrates it as a new promising target for augmenting crop resilience.

The Potential Role of Genic-SSRs in Driving Ecological Adaptation Diversity in Caragana Plants

Caragana, a xerophytic shrub genus widely distributed in northern China, exhibits distinctive geographical substitution patterns and ecological adaptation diversity. This study employed transcriptome sequencing technology to investigate 12 Caragana species, aiming to explore genic-SSR variations in the Caragana transcriptome and identify their role as a driving force for environmental adaptation within the genus. A total of 3666 polymorphic genic-SSRs were identified across different species. The impact of these variations on the expression of related genes was analyzed, revealing a significant linear correlation (p < 0.05) between the length variation of 264 polymorphic genic-SSRs and the expression of associated genes. Additionally, 2424 polymorphic genic-SSRs were located in differentially expressed genes among Caragana species. Through weighted gene co-expression network analysis, the expressions of these genes were correlated with 19 climatic factors and 16 plant functional traits in various habitats. This approach facilitated the identification of biological processes associated with habitat adaptations in the studied Caragana species. Fifty-five core genes related to functional traits and climatic factors were identified, including various transcription factors such as MYB, TCP, ARF, and structural proteins like HSP90, elongation factor TS, and HECT. The roles of these genes in the ecological adaptation diversity of Caragana were discussed. Our study identified specific genomic components and genes in Caragana plants responsive to heterogeneous habitats. The results contribute to advancements in the molecular understanding of their ecological adaptation, lay a foundation for the conservation and development of Caragana germplasm resources, and provide a scientific basis for plant adaptation to global climate change.

GWAS and Meta-QTL Analysis of Yield-Related Ear Traits in Maize

Maize ear traits are an important component of yield, and the genetic basis of ear traits facilitates further yield improvement. In this study, a panel of 580 maize inbred lines were used as the study material, eight ear-related traits were measured through three years of planting, and whole genome sequencing was performed using the maize 40 K breeding chip based on genotyping by targeted sequencing (GBTS) technology. Five models were used to conduct a genome-wide association study (GWAS) on best linear unbiased estimate (BLUE) of ear traits to find the best model. The FarmCPU (Fixed and random model Circulating Probability Unification) model was the best model for this study; a total of 104 significant single nucleotide polymorphisms (SNPs) were detected, and 10 co-location SNPs were detected simultaneously in more than two environments. Through gene function annotation and prediction, a total of nine genes were identified as potentially associated with ear traits. Moreover, a total of 760 quantitative trait loci (QTL) associated with yield-related traits reported in 37 different articles were collected. Using the collected 760 QTL for meta-QTL analysis, a total of 41 MQTL (meta-QTL) associated with yield-related traits were identified, and 19 MQTL detected yield-related ear trait functional genes and candidate genes that have been reported in maize. Five significant SNPs detected by GWAS were located within these MQTL intervals, and another three significant SNPs were close to MQTL (less than 1 Mb). The results provide a theoretical reference for the analysis of the genetic basis of ear-related traits and the improvement of maize yield.

QTL Verification and Candidate Gene Screening of Fiber Quality and Lint Percentage in the Secondary Segregating Population of Gossypium hirsutum

Fiber quality traits, especially fiber strength, length, and micronaire (FS, FL, and FM), have been recognized as critical fiber attributes in the textile industry, while the lint percentage (LP) was an important indicator to evaluate the cotton lint yield. So far, the genetic mechanism behind the formation of these traits is still unclear. Quantitative trait loci (QTL) identification and candidate gene validation provide an effective methodology to uncover the genetic and molecular basis of FL, FS, FM, and LP. A previous study identified three important QTL/QTL cluster loci, harboring at least one of the above traits on chromosomes A01, A07, and D12 via a recombinant inbred line (RIL) population derived from a cross of Lumianyan28 (L28) × Xinluzao24 (X24). A secondary segregating population (F2) was developed from a cross between L28 and an RIL, RIL40 (L28 × RIL40). Based on the population, genetic linkage maps of the previous QTL cluster intervals on A01 (6.70-10.15 Mb), A07 (85.48-93.43 Mb), and D12 (0.40-1.43 Mb) were constructed, which span 12.25, 15.90, and 5.56 cM, with 2, 14, and 4 simple sequence repeat (SSR) and insertion/deletion (Indel) markers, respectively. QTLs of FL, FS, FM, and LP on these three intervals were verified by composite interval mapping (CIM) using WinQTL Cartographer 2.5 software via phenotyping of F2 and its derived F2:3 populations. The results validated the previous primary QTL identification of FL, FS, FM, and LP. Analysis of the RNA-seq data of the developing fibers of L28 and RIL40 at 10, 20, and 30 days post anthesis (DPA) identified seven differentially expressed genes (DEGs) as potential candidate genes. qRT-PCR verified that five of them were consistent with the RNA-seq result. These genes may be involved in regulating fiber development, leading to the formation of FL, FS, FM, and LP. This study provides an experimental foundation for further exploration of these functional genes to dissect the genetic mechanism of cotton fiber development.

Integration of GWAS and RNA-Seq Analysis to Identify SNPs and Candidate Genes Associated with Alkali Stress Tolerance at the Germination Stage in Mung Bean

Soil salt-alkalization seriously impacts crop growth and productivity worldwide. Breeding and applying tolerant varieties is the most economical and effective way to address soil alkalization. However, genetic resources for breeders to improve alkali tolerance are limited in mung bean. Here, a genome-wide association study (GWAS) was performed to detect alkali-tolerant genetic loci and candidate genes in 277 mung bean accessions during germination. Using the relative values of two germination traits, 19 QTLs containing 32 SNPs significantly associated with alkali tolerance on nine chromosomes were identified, and they explained 3.6 to 14.6% of the phenotypic variance. Moreover, 691 candidate genes were mined within the LD intervals containing significant trait-associated SNPs. Transcriptome sequencing of alkali-tolerant accession 132-346 under alkali and control conditions after 24 h of treatment was conducted, and 2565 DEGs were identified. An integrated analysis of the GWAS and DEGs revealed six hub genes involved in alkali tolerance responses. Moreover, the expression of hub genes was further validated by qRT-PCR. These findings improve our understanding of the molecular mechanism of alkali stress tolerance and provide potential resources (SNPs and genes) for the genetic improvement of alkali tolerance in mung bean.

Genomic Survey of Flavin Monooxygenases in Wild and Cultivated Rice Provides Insight into Evolution and Functional Diversities

The flavin monooxygenase (FMO) enzyme was discovered in mammalian liver cells that convert a carcinogenic compound, N-N'-dimethylaniline, into a non-carcinogenic compound, N-oxide. Since then, many FMOs have been reported in animal systems for their primary role in the detoxification of xenobiotic compounds. In plants, this family has diverged to perform varied functions like pathogen defense, auxin biosynthesis, and S-oxygenation of compounds. Only a few members of this family, primarily those involved in auxin biosynthesis, have been functionally characterized in plant species. Thus, the present study aims to identify all the members of the FMO family in 10 different wild and cultivated Oryza species. Genome-wide analysis of the FMO family in different Oryza species reveals that each species has multiple FMO members in its genome and that this family is conserved throughout evolution. Taking clues from its role in pathogen defense and its possible function in ROS scavenging, we have also assessed the involvement of this family in abiotic stresses. A detailed in silico expression analysis of the FMO family in Oryza sativa subsp. japonica revealed that only a subset of genes responds to different abiotic stresses. This is supported by the experimental validation of a few selected genes using qRT-PCR in stress-sensitive Oryza sativa subsp. indica and stress-sensitive wild rice Oryza nivara. The identification and comprehensive in silico analysis of FMO genes from different Oryza species carried out in this study will serve as the foundation for further structural and functional studies of FMO genes in rice as well as other crop types.

Association mapping unravels the genetics controlling seedling drought stress tolerance in winter wheat

Drought is a major constraint in wheat (Triticum aestivum L.) grain yield. The present work aimed to identify quantitative trait nucleotides (QTNs)/ candidate genes influencing drought tolerance-related traits at the seedling stage in 261 accessions of a diverse winter wheat panel. Seeds from three consecutive years were exposed to polyethylene glycol 12% (PEG-6000) and a control treatment (distilled water). The Farm-CPU method was used for the association analysis with 17,093 polymorphic SNPs. PEG treatment reduced shoot length (SL) (-36.3%) and root length (RL) (-11.3%) compared with control treatments, while the coleoptile length (CL) was increased by 11% under drought conditions, suggesting that it might be considered as an indicator of stress-tolerance. Interestingly, we revealed 70 stable QTN across 17 chromosomes. Eight QTNs related to more than one trait were detected on chromosomes 1B, 2A (2), 2B, 2D, 4B, 7A, and 7B and located nearby or inside candidate genes within the linkage disequilibrium (LD) interval. For instance, the QTN on chromosome 2D is located inside the gene TraesCS2D02G133900 that controls the variation of CL_S and SL_C. The allelic variation at the candidate genes showed significant influence on the associated traits, demonstrating their role in controlling the natural variation of multi-traits of drought stress tolerance. The gene expression of these candidate genes under different stress conditions validates their biological role in stress tolerance. Our findings offer insight into understanding the genetic factors and diverse mechanisms in response to water shortage conditions that are important for wheat improvement and adaptation at early developmental stages.

Genome-wide identification of cystathionine beta synthase genes in wheat and its relationship with anther male sterility under heat stress

Cystathionine beta synthase (CBS) domains containing proteins (CDCPs) plays an important role in plant development through regulation of the thioredoxin system, as well as its ability to respond to biotic and abiotic stress conditions. Despite this, no systematic study has examined the wheat CBS gene family and its relation to high temperature-induced male sterility. In this study, 66 CBS family members were identified in the wheat genome, and their gene or protein sequences were used for subsequent analysis. The TaCBS gene family was found to be unevenly distributed on 21 chromosomes, and they were classified into four subgroups according to their gene structure and phylogeny. The results of collinearity analysis showed that there were 25 shared orthologous genes between wheat, rice and Brachypodium distachyon, and one shared orthologous gene between wheat, millet and barley. The cis-regulatory elements of the TaCBS were related to JA, IAA, MYB, etc. GO and KEGG pathway analysis identified these TaCBS genes to be associated with pollination, reproduction, and signaling and cellular processes, respectively. A heatmap of wheat plants based on transcriptome data showed that TaCBS genes were expressed to a higher extent in spikelets relative to other tissues. In addition, 29 putative tae-miRNAs were identified, targeting 41 TaCBS genes. Moreover, qRT-PCR validation of six TaCBS genes indicated their critical role in anther development, as five of them were expressed at lower levels in heat-stressed male sterile anthers than in Normal anthers. Together with anther phenotypes, paraffin sections, starch potassium iodide staining, and qRT-PCR data, we hypothesized that the TaCBS gene has a very important connection with the heat-stressed sterility process in wheat, and these data provide a basis for further insight into their relationship.

Genome-wide characterization of cys-tathionine-β-synthase domain-containing proteins in sugarcane reveals their role in defense responses under multiple stressors

Cys-tathionine-β-synthase (CBS) domain-containing proteins (CDCPs) are essential for regulating plant responses to various biotic and abiotic stressors. This study describes the systematic identification and characterization of CDCP family genes in Saccharum spontaneum. A total of 95 SsCDCP genes and eight phylogenetic groups were identified that were distributed over 29 chromosomes of the AP85-441 genome. Most (78/95) SsCDCPs underwent fragment duplication events, and 64 gene pairs were located in synteny blocks. Expression profiling of nine ShCDCPs was also carried out in the Saccharum spp. cultivars ROC22 and MT11-611 that are resistant and susceptible to red stripe, respectively, in response to: (i) Infection by the bacterial pathogen Acidovorax avenue subsp. avenae (Aaa); (ii) abiotic stressors (drought and salinity); and (iii) exogenous salicylic acid (SA) treatment. Members of one gene pair (ShCBSD-PB1-5A and ShCBSD-PB1-7A-1) with a fragment duplication event acted as negative regulators in sugarcane under four stresses, as supported by the significantly decreased expression levels of ShCBSD-PB1-5A (23-83%) and ShCBSD-PB1-7A-1 (15-75%) at all-time points, suggesting that they have functional redundancy. Genes in another pair, ShCBS-4C and ShCBS-4D-1, which have a fragment duplication event, play opposing regulatory roles in sugarcane exposed to multiple stresses, particularly Aaa and NaCl treatments. ShCBS-4C expression was significantly decreased by 32-77%, but ShCBS-4D-1 expression was dramatically upregulated by 1.2-6.2-fold in response to Aaa treatment of both cultivars across all-time points. This result suggested that both genes exhibited functional divergence. Meanwhile, the expression of SsCBSDCBS-5A was significantly upregulated in ROC22 by 1.4-4.6-fold in response to the four stressors. These findings provide important clues for further elucidating the function of ShCDCP genes in sugarcane responding to a diverse range of stresses.