Genetic Diversity of Transcription Factor Genes in Triticum and Mining for Promising Haplotypes for Beneficial Agronomic Traits

A system genetics analysis uncovers the regulatory variants controlling drought response in wheat

Plants activate a variable response to drought stress by modulating transcription of key genes. However, our knowledge of genetic variations governing gene expression in response to drought stress remains limited in natural germplasm. Here, we performed a comprehensive analysis of the transcriptional variability of 200 wheat accessions in response to drought stress by using a systems genetics approach integrating pan-transcriptome, co-expression networks, transcriptome-wide association study (TWAS), and expression quantitative trait loci (eQTLs) mapping. We identified 1621 genes and eight co-expression modules significantly correlated with wheat drought tolerance. We also defined 620 664 and 654 798 independent eQTLs associated with the expression of 17 429 and 18 080 eGenes under normal and drought stress conditions. Focusing on dynamic regulatory variants, we further identified 572 eQTL hotspots and constructed transcription factors governed drought-responsive network by the XGBoost model. Subsequently, by combining with genome-wide association study (GWAS), we uncovered a 369-bp insertion variant in the TaKCS3 promoter containing multiple cis-regulatory elements recognized by eQTL hotspot-associated transcription factors that enhance its transcription. Further functional analysis indicated that elevating TaKCS3 expression affects cuticular wax composition to reduce water loss during drought stress, and thereby increase drought tolerance. This study sheds light on the genome-wide genetic variants that influence dynamic transcriptional changes during drought stress and provides a valuable resource for the mining of drought-tolerant genes in the future.

Comprehensive analysis of computational approaches in plant transcription factors binding regions discovery

Transcription factors (TFs) are regulatory proteins which bind to a specific DNA region known as the transcription factor binding regions (TFBRs) to regulate the rate of transcription process. The identification of TFBRs has been made possible by a number of experimental and computational techniques established during the past few years. The process of TFBR identification involves peak identification in the binding data, followed by the identification of motif characteristics. Using the same binding data attempts have been made to raise computational models to identify such binding regions which could save time and resources spent for binding experiments. These computational approaches depend a lot on what way they learn and how. These existing computational approaches are skewed heavily around human TFBRs discovery, while plants have drastically different genomic setup for regulation which these approaches have grossly ignored. Here, we provide a comprehensive study of the current state of the matters in plant specific TF discovery algorithms. While doing so, we encountered several software tools' issues rendering the tools not useable to researches. We fixed them and have also provided the corrected scripts for such tools. We expect this study to serve as a guide for better understanding of software tools' approaches for plant specific TFBRs discovery and the care to be taken while applying them, especially during cross-species applications. The corrected scripts of these software tools are made available at https://github.com/SCBB-LAB/Comparative-analysis-of-plant-TFBS-software.

Genome-wide identification of m6A-related gene family and the involvement of TdFIP37 in salt stress in wild emmer wheat

KEY MESSAGE

The genomic organization, phylogenetic relationship, expression patterns, and genetic variations of m6A-related genes were systematically investigated in wild emmer wheat and the function of TdFIP37 regulating salt tolerance was preliminarily determined. m6A modification is one of the most abundant and crucial RNA modifications in eukaryotics, playing the indispensable role in growth and development as well as stress response in plants. However, its significance in wild emmer wheat remains elusive. Here, a genome-wide search of m6A-related genes was conducted in wild emmer wheat to obtain 64 candidates, including 21 writers, 17 erasers, and 26 readers. Phylogenetic and collinearity analysis demonstrated that segmental duplication and polyploidization contributed mainly to the expansion of m6A-related genes in wild emmer. A number of cis-acting elements involving in stress and hormonal regulation were found in the promoter regions of them, such as MBS, LTR, and ABRE. Genetic variation of them was also investigated using resequencing data and obvious genetic bottleneck was occurred on them during wild emmer wheat domestication process. Furthermore, the salt-responsive candidates were investigated through RNA-seq data and qRT-PCR validation using the salt-tolerant and -sensitive genotypes and the co-expression analysis showed that they played the hub role in regulating salt stress response. Finally, the loss-function mutant of Tdfip37 displayed the significantly higher salt-sensitive compared to WT and then RNA-seq analysis demonstrated that FIP37 mediated the MAPK pathway, hormone signal transduction, as well as transcription factor to regulate salt tolerance. This study provided the potential m6A genes for functional analysis, which will contribute to better understand the regulatory roles of m6A modification and also improve the salt tolerance from the perspective of epigenetic approach in emmer wheat and other crops.

Genetic insights into superior grain number traits: a QTL analysis of wheat-Agropyron cristatum derivative pubing3228

BACKGROUND

Agropyron cristatum (L.) is a valuable genetic resource for expanding the genetic diversity of common wheat. Pubing3228, a novel wheat-A. cristatum hybrid germplasm, exhibits several desirable agricultural traits, including high grain number per spike (GNS). Understanding the genetic architecture of GNS in Pubing3228 is crucial for enhancing wheat yield. This study aims to analyze the specific genetic regions and alleles associated with high GNS in Pubing3228.

METHODS

The study employed a recombination inbred line (RIL) population derived from a cross between Pubing3228 and Jing4839 to investigate the genetic regions and alleles linked to high GNS. Quantitative Trait Loci (QTL) analysis and candidate gene investigation were utilized to explore these traits.

RESULTS

A total of 40 QTLs associated with GNS were identified across 16 chromosomes, accounting for 4.25-17.17% of the total phenotypic variation. Five QTLs (QGns.wa-1D, QGns.wa-5 A, QGns.wa-7Da.1, QGns.wa-7Da.2 and QGns.wa-7Da.3) accounter for over 10% of the phenotypic variation in at least two environments. Furthermore, 94.67% of the GNS QTL with positive effects originated from Pubing3228. Candidate gene analysis of stable QTLs identified 11 candidate genes for GNS, including a senescence-associated protein gene (TraesCS7D01G148000) linked to the most significant SNP (AX-108,748,734) on chromosome 7D, potentially involved in reallocating nutrients from senescing tissues to developing seeds.

CONCLUSION

This study provides new insights into the genetic mechanisms underlying high GNS in Pubing3228, offering valuable resources for marker-assisted selection in wheat breeding to enhance yield.

Genome-Wide Identification, Expression and Evolution Analysis of m6A Writers, Readers and Erasers in Aegilops_tauschii

N6-methyladenosine modifications (m6A) is one of the most abundant and prevalent post-transcriptional RNA modifications in plants, playing the crucial role in plant growth and development and stress adaptation. However, the m6A regulatory machinery in Aegilops_tauschii, the D genome progenitor of common wheat, is not well understood at present. Here, we systematically identified the m6A-related genes in Aegilops with a genome-wide search approach. In total, 25 putative m6A genes composed of 5 writers, 13 readers and 7 erasers were obtained. A phylogenetic analysis clearly grouped them into three subfamilies with the same subfamily showing similar gene structures and conserved domains. These m6A genes were found to contain a large number of cis-acting elements associating with plant hormones, regulation of growth and development as well as stress response, suggesting their widespread regulation function. Furthermore, the expression profiling of them was investigated using RNA-seq data to obtain stress-responsive candidates, of which 5 were further validated with a qPCR analysis. Finally, the genetic variation of m6A-related genes was investigated between Aegilops and D subgenome of wheat based on re-sequencing data, and an obvious genetic bottleneck occurred on them during the wheat domestication process. The promising haplotype association with domestication and agronomic traits was also detected. This study provided some insights on the genomic organization and evolutionary features of m6A-related genes in Aegilops, which will facilitate the further functional study and also contribute to broaden the genetic basis for genetic improvement in wheat and other crops.

The evolutionary patterns, expression profiles, and genetic diversity of expanded genes in barley

Gene duplication resulting from whole-genome duplication (WGD), small-scale duplication (SSD), or unequal hybridization plays an important role in the expansion of gene families. Gene family expansion can also mediate species formation and adaptive evolution. Barley (Hordeum vulgare) is the world's fourth largest cereal crop, and it contains valuable genetic resources due to its ability to tolerate various types of environmental stress. In this study, 27,438 orthogroups in the genomes of seven Poaceae were identified, and 214 of them were significantly expanded in barley. The evolutionary rates, gene properties, expression profiles, and nucleotide diversity between expanded and non-expanded genes were compared. Expanded genes evolved more rapidly and experienced lower negative selection. Expanded genes, including their exons and introns, were shorter, they had fewer exons, their GC content was lower, and their first exons were longer compared with non-expanded genes. Codon usage bias was also lower for expanded genes than for non-expanded genes; the expression levels of expanded genes were lower than those of non-expanded genes, and the expression of expanded genes showed higher tissue specificity than that of non-expanded genes. Several stress-response-related genes/gene families were identified, and these genes could be used to breed barley plants with greater resistance to environmental stress. Overall, our analysis revealed evolutionary, structural, and functional differences between expanded and non-expanded genes in barley. Additional studies are needed to clarify the functions of the candidate genes identified in our study and evaluate their utility for breeding barley plants with greater stress resistance.

The barley DIR gene family: An expanded gene family that is involved in stress responses

Gene family expansion plays a central role in adaptive divergence and, ultimately, speciation is influenced by phenotypic diversity in different environments. Barley (Hordeum vulgare) is the fourth most important cereal crop in the world and is used for brewing purposes, animal feed, and human food. Systematic characterization of expanded gene families is instrumental in the research of the evolutionary history of barley and understanding of the molecular function of their gene products. A total of 31,750 conserved orthologous groups (OGs) were identified using eight genomes/subgenomes, of which 1,113 and 6,739 were rapidly expanded and contracted OGs in barley, respectively. Five expanded OGs containing 20 barley dirigent genes (HvDIRs) were identified. HvDIRs from the same OG were phylogenetically clustered with similar gene structure and domain organization. In particular, 7 and 5 HvDIRs from OG0000960 and OG0001516, respectively, contributed greatly to the expansion of the DIR-c subfamily. Tandem duplication was the driving force for the expansion of the barley DIR gene family. Nucleotide diversity and haplotype network analysis revealed that the expanded HvDIRs experienced severe bottleneck events during barley domestication, and can thus be considered as potential domestication-related candidate genes. The expression profile and co-expression network analysis revealed the critical roles of the expanded HvDIRs in various biological processes, especially in stress responses. HvDIR18, HvDIR19, and HvDIR63 could serve as excellent candidates for further functional genomics studies to improve the production of barley products. Our study revealed that the HvDIR family was significantly expanded in barley and might be involved in different developmental processes and stress responses. Thus, besides providing a framework for future functional genomics and metabolomics studies, this study also identified HvDIRs as candidates for use in improving barley crop resistance to biotic and abiotic stresses.

Genome-Wide Identification, Evolution and Expressional Analysis of OSCA Gene Family in Barley (Hordeum vulgare L.)

The hyperosmolality-gated calcium-permeable channel gene family (OSCA) is one kind of conserved osmosensors, playing a crucial role in maintaining ion and water homeostasis and protecting cellular stability from the damage of hypertonic stress. Although it has been systematically characterized in diverse plants, it is necessary to explore the role of the OSCA family in barley, especially its importance in regulating abiotic stress response. In this study, a total of 13 OSCA genes (HvOSCAs) were identified in barley through an in silico genome search method, which were clustered into 4 clades based on phylogenetic relationships with members in the same clade showing similar protein structures and conserved motif compositions. These HvOSCAs had many cis-regulatory elements related to various abiotic stress, such as MBS and ARE, indicating their potential roles in abiotic stress regulation. Furthermore, their expression patterns were systematically detected under diverse stresses using RNA-seq data and qRT-PCR methods. All of these 13 HvOSCAs were significantly induced by drought, cold, salt and ABA treatment, demonstrating their functions in osmotic regulation. Finally, the genetic variations of the HvOSCAs were investigated using the re-sequencing data, and their nucleotide diversity in wild barley and landrace populations were 0.4966 × 10-3 and 0.391 × 10-3, respectively, indicating that a genetic bottleneck has occurred in the OSCA family during the barley evolution process. This study evaluated the genomic organization, evolutionary relationship and genetic expression of the OSCA family in barley, which not only provides potential candidates for further functional genomic study, but also contributes to genetically improving stress tolerance in barley and other crops.

Population transcriptomic analysis identifies the comprehensive lncRNAs landscape of spike in wheat (Triticum aestivum L.)

BACKGROUND

Long noncoding RNAs (lncRNAs) are emerging as the important regulators involving in growth and development as well as stress response in plants. However, current lncRNA studies were mainly performed at the individual level and the significance of it is not well understood in wheat.

RESULTS

In this study, the lncRNA landscape of wheat spike was characterized through analysing a total of 186 spike RNA-seq datasets from 93 wheat genotypes. A total of 35,913 lncRNAs as well as 1,619 lncRNA-mRNA pairs comprised of 443 lncRNAs and 464 mRNAs were obtained. Compared to coding genes, these lncRNAs displayed rather low conservation among wheat and other gramineous species. Based on re-sequencing data, the genetic variations of these lncRNA were investigated and obvious genetic bottleneck were found on them during wheat domestication process. Furthermore, 122 lncRNAs were found to act as ceRNA to regulate endogenous competition. Finally, association and co-localization analysis of the candidate lncRNA-mRNA pairs identified 170 lncRNAs and 167 target mRNAs significantly associated with spike-related traits, including lncRNA.127690.1/TraesCS2A02G518500.1 (PMEI) and lncRNA.104854.1/TraesCS6A02G050300.1 (ATG5) associated with heading date and spike length, respectively.

CONCLUSIONS

This study reported the lncRNA landscape of wheat spike through the population transcriptome analysis, which not only contribute to better understand the wheat evolution from the perspective of lncRNA, but also lay the foundation for revealing roles of lncRNA playing in spike development.