Modeling of transcriptomic variation among subgenomes in 25 accessions of common wheat reveals cis- and trans-regulation architectures

Yasuyuki Nomura,Moeko Okada,Toshiaki Tameshige,Shotaro Takenaka,Kentaro K Shimizu,Shuhei Nasuda,Atsushi J Nagano
DOI: https://doi.org/10.1101/2024.10.09.617331
2024-10-13
Abstract:Common wheat is an allohexaploid plant, thus making it difficult to obtain homoeolog-distinguished transcriptome data. Lasy-Seq, a type of 3' RNA-seq, is efficient for obtaining homoeolog-distinguished transcriptomes and can thus overcome this measurement difficulty. This study obtained transcriptome data from the seedlings, second leaves, and root tips of 25 lines from mainly eastern transmitted area using Lasy-Seq. Roots and seedlings exhibited similar transcriptome profiles; however, they were different from those of the leaves. We determined the effects of subgenomes, lines and their interactions with leaves, roots, and seedlings on the expression levels of each homoeolog triad. Of the 19,805 homoeolog triads, 50.9-55.4%, 24.2-29.5%, and 7.7-9.0% showed significant effects on their expression levels from subgenome, line, and interaction, respectively. 51-55% and 24-30% have genetic variation in the cis- and trans-regulation. Hierarchical clustering and co-trans regulation network analysis of homoeolog triads revealed that the patterns of expression polymorphisms among the lines were shared in different genes. The triads in which the statistical model detected as line effects imply that expression variation between lines is caused by changes in a smaller number of common trans-factors. We assigned gene ontology (GO) terms of the Arabidopsis orthologs to wheat homoeolog triads via reciprocal BLAST between common wheat and Arabidopsis, thus improving the percentage of gene-assigned GO terms to all analyzed GO terms from 19.1% to 90.6%. GO term enrichment analysis revealed that GO terms related to each tissue type function were enriched in genes expressed in the leaves and roots. Our information provides fundamental knowledge for the future breeding of plants possessing complex gene regulatory networks such as common wheat.
Plant Biology
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to model the transcriptome variation in different sub - genomes of common wheat (Triticum aestivum) in order to reveal the cis - and trans - regulatory structures. Specifically, the research aims to: 1. **Obtain transcriptome data to distinguish homologous triplets**: Use Lasy - Seq technology to obtain transcriptome data from the seedlings, the second leaves and root tips of 25 common wheat lines, which are mainly from the eastern spread area. 2. **Analyze the transcriptome characteristics of different tissue types**: The study found that the transcriptome profiles of roots and seedlings are similar, but different from those of leaves. 3. **Determine the effects of sub - genomes, lines and their interactions on the expression levels of homologous triplets**: Through the analysis of 19,805 homologous triplets, it was found that 50.9% - 55.4% of homologous triplets have significant expression differences among sub - genomes, 24.2% - 29.5% have significant expression differences among lines, and 7.7% - 9.0% have significant expression differences under the interaction of sub - genomes and lines. 4. **Explore the genetic variation of cis - and trans - regulation**: 51% - 55% of homologous triplets have genetic variation in cis - regulation, and 24% - 30% have genetic variation in trans - regulation. 5. **Construct a co - expression network of homologous triplets**: Through hierarchical clustering and co - expression network analysis, the patterns of expression polymorphism among different lines were revealed, and it was found that these patterns are shared among different genes. 6. **Improvement of Gene Ontology (GO) annotation**: By reciprocal BLAST with Arabidopsis thaliana, GO terms were assigned to the homologous triplets of wheat, thus increasing the proportion of GO terms in gene annotation from 19.1% to 90.6%. 7. **GO term enrichment analysis**: GO term enrichment analysis was carried out on differentially expressed genes (DEGs) and differentially expressed homologous triplets (DETs), and it was found that GO terms related to the functions of specific tissue types were enriched in genes expressed in leaves and roots. In general, this study provides basic knowledge about the complex gene regulation network of common wheat and is helpful for future plant breeding work.