Genomic variants affecting homoeologous gene expression dosage contribute to agronomic trait variation in allopolyploid wheat

MRI-Seed-Wizard: combining deep learning algorithms with magnetic resonance imaging enables advanced seed phenotyping

Evaluation of relevant seed traits is an essential part of most plant breeding and biotechnology programmes. There is a need for non-destructive, three-dimensional assessment of the morphometry, composition, and internal features of seeds. Here, we introduce a novel tool, MRI-Seed-Wizard, which integrates deep learning algorithms with non-invasive magnetic resonance imaging (MRI) for use in a new domain-plant MRI. The tool enabled in vivo quantification of 23 grain traits, including volumetric parameters of inner seed structure. Several of these features cannot be assessed using conventional techniques, including X-ray computed tomography. MRI-Seed-Wizard was designed to automate the manual processes of identifying, labeling, and analysing digital MRI data. We further provide advanced MRI protocols that allow the evaluation of multiple seeds simultaneously to increase throughput. The versatility of MRI-Seed-Wizard in seed phenotyping is demonstrated for wheat (Triticum aestivum) and barley (Hordeum vulgare) grains, and it is applicable to a wide range of crop seeds. Thus, artificial intelligence, combined with the most versatile imaging modality, MRI, opens up new perspectives in seed phenotyping and crop improvement.

Regulatory networks of coresident subgenomes during rapid fiber cell elongation in upland cotton

Cotton, an intriguing plant species shaped by polyploidization, evolution, and domestication, holds particular interest due to the complex mechanisms governing fiber traits across its two subgenomes. However, the regulatory elements or transcriptional networks between subgenomes during fiber elongation remain to be fully clarified. Here, we analyzed 1462 cotton fiber samples to reconstruct the gene-expression regulatory networks that influence fiber cell elongation. Inter-subgenome expression quantitative trait loci (eQTLs) largely dictate gene transcription, with a notable tendency for the D subgenome to regulate A-subgenome eGenes. This regulation reveals synchronized homoeologous gene expression driven by co-localized eQTLs and divergent patterns that diminish genetic correlations, thus leading to preferential expression in the A and D subgenomes. Hotspot456 emerged as a key regulator of fiber initiation and elongation, and artificial selection of trans-eQTLs in hotspot456 that positively regulate KCS1 has facilitated cell elongation. Experiments designed to clarify the roles of trans-eQTLs in improved fiber breeding confirmed the inhibition of GhTOL9 by a specific trans-eQTL via GhWRKY28, which negatively affects fiber elongation. We propose a model in which the GhWRKY28-GhTOL9 module regulates this process through the ESCRT (endosomal sorting complex required for transport) pathway. This research significantly advances our understanding of cotton's evolutionary and domestication processes and the intricate regulatory mechanisms that underlie significant plant traits.

Population-wide DNA methylation polymorphisms at single-nucleotide resolution in 207 cotton accessions reveal epigenomic contributions to complex traits

DNA methylation plays multiple regulatory roles in crop development. However, the relationships of methylation polymorphisms with genetic polymorphisms, gene expression, and phenotypic variation in natural crop populations remain largely unknown. Here, we surveyed high-quality methylomes, transcriptomes, and genomes obtained from the 20-days-post-anthesis (DPA) cotton fibers of 207 accessions and extended the classical framework of population genetics to epigenetics. Over 287 million single methylation polymorphisms (SMPs) were identified, 100 times more than the number of single nucleotide polymorphisms (SNPs). These SMPs were significantly enriched in intragenic regions while depleted in transposable elements. Association analysis further identified a total of 5,426,782 cis-methylation quantitative trait loci (cis-meQTLs), 5078 cis-expression quantitative trait methylation (cis-eQTMs), and 9157 expression quantitative trait loci (eQTLs). Notably, 36.39% of cis-eQTM genes were not associated with genetic variation, indicating that a large number of SMPs associated with gene expression variation are independent of SNPs. In addition, out of the 1715 epigenetic loci associated with yield and fiber quality traits, only 36 (2.10%) were shared with genome-wide association study (GWAS) loci. The construction of multi-omics regulatory networks revealed 43 cis-eQTM genes potentially involved in fiber development, which cannot be identified by GWAS alone. Among these genes, the role of one encoding CBL-interacting protein kinase 10 in fiber length regulation was successfully validated through gene editing. Taken together, our findings prove that DNA methylation data can serve as an additional resource for breeding purposes and can offer opportunities to enhance and expedite the crop improvement process.

A Genome-Wide Alternative Splicing Analysis of Gossypium arboreum and Gossypium raimondii During Fiber Development

Alternative splicing (AS) is a crucial post-transcriptional regulatory mechanism that contributes to proteome complexity and versatility in different plant species. However, detailed AS exploration in diploid cotton during fiber development has not been reported. In this study, we comparatively analyzed G. arboreum and G. raimondii AS events during fiber development using transcriptome data and identified 9690 and 7617 AS events that were distributed in 6483 and 4859 genes, respectively. G. arboreum had more AS genes and AS events than G. raimondii, and most AS genes were distributed at both ends of all 13 chromosomes in both diploid cotton species. Four major AS types, including IR, SE, A3SS, and A5SS, were all experimentally validated through RT-PCR assays. G. arboreum and G. raimondii had only 1888 AS genes in common, accounting for one-third and one-half of the total number of AS genes, respectively. Furthermore, we found a lysine-specific demethylase coding gene with a different AS mechanism in G. arboreum and G. raimondii, in which AS isoforms lacked part of a key conserved domain. Our findings may provide new directions for the discovery of functional genes involved in cotton species differentiation.

Genetic Regulatory Perturbation of Gene Expression Impacted by Genomic Introgression in Fiber Development of Allotetraploid Cotton

Interspecific genomic introgression is an important evolutionary process with respect to the generation of novel phenotypic diversity and adaptation. A key question is how gene flow perturbs gene expression networks and regulatory interactions. Here, an introgression population of two species of allopolyploid cotton (Gossypium) to delineate the regulatory perturbations of gene expression regarding fiber development accompanying fiber quality change is utilized. De novo assembly of the recipient parent (G. hirsutum Emian22) genome allowed the identification of genomic variation and introgression segments (ISs) in 323 introgression lines (ILs) from the donor parent (G. barbadense 3-79). It documented gene expression dynamics by sequencing 1,284 transcriptomes of developing fibers and characterized genetic regulatory perturbations mediated by genomic introgression using a multi-locus model. Introgression of individual homoeologous genes exhibiting extreme low or high expression bias can lead to a parallel expression bias in their non-introgressed duplicates, implying a shared yet divergent regulatory fate of duplicated genes following allopolyploidy. Additionally, the IL N182 with improved fiber quality is characterized, and the candidate gene GhFLAP1 related to fiber length is validated. This study outlines a framework for understanding introgression-mediated regulatory perturbations in polyploids, and provides insights for targeted breeding of superior upland cotton fiber.

Genome-wide characterization and expression analysis of the CINNAMYL ALCOHOL DEHYDROGENASE gene family in Triticum aestivum

BACKGROUND

CINNAMYL ALCOHOL DEHYDROGENASE (CAD) catalyzes the NADPH-dependent reduction of cinnamaldehydes into cinnamyl alcohols and is a key enzyme found at the final step of the monolignol pathway. Cinnamyl alcohols and their conjugates are subsequently polymerized in the secondary cell wall to form lignin. CAD genes are typically encoded by multi-gene families and thus traditionally organized into general classifications of functional relevance.

RESULTS

In silico analysis of the hexaploid Triticum aestivum genome revealed 47 high confidence TaCAD copies, of which three were determined to be the most significant isoforms (class I) considered bone fide CADs. Class I CADs were expressed throughout development both in RNAseq data sets as well as via qRT-PCR analysis. Of the 37 class II TaCADs identified, two groups were observed to be significantly co-expressed with class I TaCADs in developing tissue and under chitin elicitation in RNAseq data sets. These co-expressed class II TaCADs were also found to be phylogenetically unrelated to a separate clade of class II TaCADs previously reported to be an influential resistance factor to pathogenic fungal infection. Lastly, two groups were phylogenetically identified as class III TaCADs, which possess distinct conserved gene structures. However, the lack of data supporting their catalytic activity for cinnamaldehydes and their bereft transcriptional presence in lignifying tissues challenges their designation and function as CADs.

CONCLUSIONS

Taken together, our comprehensive transcriptomic analyses suggest that TaCAD genes contribute to overlapping but nonredundant functions during T. aestivum growth and development across a wide variety of agroecosystems and provide tolerance to various stressors.

Epigenomic identification of vernalization cis-regulatory elements in winter wheat

BACKGROUND

Winter wheat undergoes vernalization, a process activated by prolonged exposure to low temperatures. During this phase, flowering signals are generated and transported to the apical meristems, stimulating the transition to the inflorescence meristem while inhibiting tiller bud elongation. Although some vernalization genes have been identified, the key cis-regulatory elements and precise mechanisms governing this process in wheat remain largely unknown.

RESULTS

In this study, we construct extensive epigenomic and transcriptomic profiling across multiple tissues-leaf, axillary bud, and shoot apex-during the vernalization of winter wheat. Epigenetic modifications play a crucial role in eliciting tissue-specific responses and sub-genome-divergent expressions during vernalization. Notably, we observe that H3K27me3 primarily regulates vernalization-induced genes and has limited influence on vernalization-repressed genes. The integration of these datasets enables the identification of 10,600 putative vernalization-related regulatory elements including distal accessible chromatin regions (ACRs) situated 30Kb upstream of VRN3, contributing to the construction of a comprehensive regulatory network. Furthermore, we discover that TaSPL7/15, integral components of the aging-related flowering pathway, interact with the VRN1 promoter and VRN3 distal regulatory elements. These interactions finely regulate their expressions, consequently impacting the vernalization process and flowering.

CONCLUSIONS

Our study offers critical insights into wheat vernalization's epigenomic dynamics and identifies the putative regulatory elements crucial for developing wheat germplasm with varied vernalization characteristics. It also establishes a vernalization-related transcriptional network, and uncovers that TaSPL7/15 from the aging pathway participates in vernalization by directly binding to the VRN1 promoter and VRN3 distal regulatory elements.

Combined GWAS and eGWAS reveals the genetic basis underlying drought tolerance in emmer wheat (Triticum turgidum L.)

Drought is one of the major environmental constraints for wheat production world-wide. As the progenitor and genetic reservoir of common wheat, emmer wheat is considered as an invaluable gene pool for breeding drought-tolerant wheat. Combining GWAS and eGWAS analysis of 107 accessions, we identified 86 QTLs, 105 462 eQTLs as well as 68 eQTL hotspots associating with drought tolerance (DT) in emmer wheat. A complex regulatory network composed of 185 upstream regulator and 2432 downstream drought-responsive candidates was developed, of which TtOTS1 was found to play a negative effect in determining DT through affecting root development. This study sheds light on revealing the genetic basis underlying DT, which will provide the indispensable genes and germplasm resources for elite drought tolerance wheat improvement and breeding.

Drought response revealed by chromatin organization variation and transcriptional regulation in cotton

BACKGROUND

Cotton is a major world cash crop and an important source of natural fiber, oil, and protein. Drought stress is becoming a restrictive factor affecting cotton production. To facilitate the development of drought-tolerant cotton varieties, it is necessary to study the molecular mechanism of drought stress response by exploring key drought-resistant genes and related regulatory factors.

RESULTS

In this study, two cotton varieties, ZY007 (drought-sensitive) and ZY168 (drought-tolerant), showing obvious phenotypic differences under drought stress, were selected. A total of 25,898 drought-induced genes were identified, exhibiting significant enrichment in pathways related to plant stress responses. Under drought induction, At subgenome expression bias was observed at the whole-genome level, which may be due to stronger inhibition of Dt subgenome expression. A gene co-expression module that was significantly associated with drought resistance was identified. About 90% of topologically associating domain (TAD) boundaries were stable, and 6613 TAD variation events were identified between the two varieties under drought. We identified 92 genes in ZY007 and 98 in ZY168 related to chromatin 3D structural variation and induced by drought stress. These genes are closely linked to the cotton response to drought stress through canonical hormone-responsive pathways, modulation of kinase and phosphatase activities, facilitation of calcium ion transport, and other related molecular mechanisms.

CONCLUSIONS

These results lay a foundation for elucidating the molecular mechanism of the cotton drought response and provide important regulatory locus and gene resources for the future molecular breeding of drought-resistant cotton varieties.

Dissecting the molecular basis of spike traits by integrating gene regulatory networks and genetic variation in wheat

Spike architecture influences both grain weight and grain number per spike, which are the two major components of grain yield in bread wheat (Triticum aestivum L.). However, the complex wheat genome and the influence of various environmental factors pose challenges in mapping the causal genes that affect spike traits. Here, we systematically identified genes involved in spike trait formation by integrating information on genomic variation and gene regulatory networks controlling young spike development in wheat. We identified 170 loci that are responsible for variations in spike length, spikelet number per spike, and grain number per spike through genome-wide association study and meta-QTL analyses. We constructed gene regulatory networks for young inflorescences at the double ridge stage and the floret primordium stage, in which the spikelet meristem and the floret meristem are predominant, respectively, by integrating transcriptome, histone modification, chromatin accessibility, eQTL, and protein-protein interactome data. From these networks, we identified 169 hub genes located in 76 of the 170 QTL regions whose polymorphisms are significantly associated with variation in spike traits. The functions of TaZF-B1, VRT-B2, and TaSPL15-A/D in establishment of wheat spike architecture were verified. This study provides valuable molecular resources for understanding spike traits and demonstrates that combining genetic analysis and developmental regulatory networks is a robust approach for dissection of complex traits.

Introgressions lead to reference bias in wheat RNA-seq analysis

BACKGROUND

RNA-seq is a fundamental technique in genomics, yet reference bias, where transcripts derived from non-reference alleles are quantified less accurately, can undermine the accuracy of RNA-seq quantification and thus the conclusions made downstream. Reference bias in RNA-seq analysis has yet to be explored in complex polyploid genomes despite evidence that they are often a complex mosaic of wild relative introgressions, which introduce blocks of highly divergent genes.

RESULTS

Here we use hexaploid wheat as a model complex polyploid, using both simulated and experimental data to show that RNA-seq alignment in wheat suffers from widespread reference bias which is largely driven by divergent introgressed genes. This leads to underestimation of gene expression and incorrect assessment of homoeologue expression balance. By incorporating gene models from ten wheat genome assemblies into a pantranscriptome reference, we present a novel method to reduce reference bias, which can be readily scaled to capture more variation as new genome and transcriptome data becomes available.

CONCLUSIONS

This study shows that the presence of introgressions can lead to reference bias in wheat RNA-seq analysis. Caution should be exercised by researchers using non-sample reference genomes for RNA-seq alignment and novel methods, such as the one presented here, should be considered.

Asymmetric and parallel subgenome selection co-shape common carp domestication

BACKGROUND

The common carp (Cyprinus carpio) might best represent the domesticated allopolyploid animals. Although subgenome divergence which is well-known to be a key to allopolyploid domestication has been comprehensively characterized in common carps, the link between genetic architecture underlying agronomic traits and subgenome divergence is unknown in the selective breeding of common carps globally.

RESULTS

We utilized a comprehensive SNP dataset in 13 representative common carp strains worldwide to detect genome-wide genetic variations associated with scale reduction, vibrant skin color, and high growth rate in common carp domestication. We identified numerous novel candidate genes underlie the three agronomically most desirable traits in domesticated common carps, providing potential molecular targets for future genetic improvement in the selective breeding of common carps. We found that independently selective breeding of the same agronomic trait (e.g., fast growing) in common carp domestication could result from completely different genetic variations, indicating the potential advantage of allopolyploid in domestication. We observed that candidate genes associated with scale reduction, vibrant skin color, and/or high growth rate are repeatedly enriched in the immune system, suggesting that domestication of common carps was often accompanied by the disease resistance improvement.

CONCLUSIONS

In common carp domestication, asymmetric subgenome selection is prevalent, while parallel subgenome selection occurs in selective breeding of common carps. This observation is not due to asymmetric gene retention/loss between subgenomes but might be better explained by reduced pleiotropy through transposable element-mediated expression divergence between ohnologs. Our results demonstrate that domestication benefits from polyploidy not only in plants but also in animals.

The Dissection of Nitrogen Response Traits Using Drone Phenotyping and Dynamic Phenotypic Analysis to Explore N Responsiveness and Associated Genetic Loci in Wheat

Inefficient nitrogen (N) utilization in agricultural production has led to many negative impacts such as excessive use of N fertilizers, redundant plant growth, greenhouse gases, long-lasting toxicity in ecosystem, and even effect on human health, indicating the importance to optimize N applications in cropping systems. Here, we present a multiseasonal study that focused on measuring phenotypic changes in wheat plants when they were responding to different N treatments under field conditions. Powered by drone-based aerial phenotyping and the AirMeasurer platform, we first quantified 6 N response-related traits as targets using plot-based morphological, spectral, and textural signals collected from 54 winter wheat varieties. Then, we developed dynamic phenotypic analysis using curve fitting to establish profile curves of the traits during the season, which enabled us to compute static phenotypes at key growth stages and dynamic phenotypes (i.e., phenotypic changes) during N response. After that, we combine 12 yield production and N-utilization indices manually measured to produce N efficiency comprehensive scores (NECS), based on which we classified the varieties into 4 N responsiveness (i.e., N-dependent yield increase) groups. The NECS ranking facilitated us to establish a tailored machine learning model for N responsiveness-related varietal classification just using N-response phenotypes with high accuracies. Finally, we employed the Wheat55K SNP Array to map single-nucleotide polymorphisms using N response-related static and dynamic phenotypes, helping us explore genetic components underlying N responsiveness in wheat. In summary, we believe that our work demonstrates valuable advances in N response-related plant research, which could have major implications for improving N sustainability in wheat breeding and production.

Deciphering the evolution and complexity of wheat germplasm from a genomic perspective

Bread wheat provides an essential fraction of the daily calorific intake for humanity. Due to its huge and complex genome, progress in studying on the wheat genome is substantially trailed behind those of the other two major crops, rice and maize, for at least a decade. With rapid advances in genome assembling and reduced cost of high-throughput sequencing, emerging de novo genome assemblies of wheat and whole-genome sequencing data are leading to a paradigm shift in wheat research. Here, we review recent progress in dissecting the complex genome and germplasm evolution of wheat since the release of the first high-quality wheat genome. New insights have been gained in the evolution of wheat germplasm during domestication and modern breeding progress, genomic variations at multiple scales contributing to the diversity of wheat germplasm, and complex transcriptional and epigenetic regulations of functional genes in polyploid wheat. Genomics databases and bioinformatics tools meeting the urgent needs of wheat genomics research are also summarized. The ever-increasing omics data, along with advanced tools and well-structured databases, are expected to accelerate deciphering the germplasm and gene resources in wheat for future breeding advances.

Dynamic gene regulatory networks improving spike fertility through regulation of floret primordia fate in wheat

The developmental process of spike is critical for spike fertility through affecting floret primordia fate in wheat; however, the genetic regulation of this dynamic and complex developmental process remains unclear. Here, we conducted a high temporal-resolution analysis of spike transcriptomes and monitored the number and morphology of floret primordia within spike. The development of all floret primordia in a spike was clearly separated into three distinct phases: differentiation, pre-dimorphism and dimorphism. Notably, we identified that floret primordia with meiosis ability at the pre-dimorphism phase usually develop into fertile floret primordia in the next dimorphism phase. Compared to control, increasing plant space treatment achieved the maximum increasement range (i.e., 50%) in number of fertile florets by accelerating spike development. The process of spike fertility improvement was directed by a continuous and dynamic regulatory network involved in transcription factor and genes interaction. This was based on the coordination of genes related to heat shock protein and jasmonic acid biosynthesis during differentiation phase, and genes related to lignin, anthocyanin and chlorophyll biosynthesis during dimorphism phase. The multi-dimensional association with high temporal-resolution approach reported here allows rapid identification of genetic resource for future breeding studies to realise the maximum spike fertility potential in more cereal crops.

Integration of eQTL and machine learning to dissect causal genes with pleiotropic effects in genetic regulation networks of seed cotton yield

The dissection of a gene regulatory network (GRN) that complements the genome-wide association study (GWAS) locus and the crosstalk underlying multiple agronomical traits remains a major challenge. In this study, we generate 558 transcriptional profiles of lint-bearing ovules at one day post-anthesis from a selective core cotton germplasm, from which 12,207 expression quantitative trait loci (eQTLs) are identified. Sixty-six known phenotypic GWAS loci are colocalized with 1,090 eQTLs, forming 38 functional GRNs associated predominantly with seed yield. Of the eGenes, 34 exhibit pleiotropic effects. Combining the eQTLs within the seed yield GRNs significantly increases the portion of narrow-sense heritability. The extreme gradient boosting (XGBoost) machine learning approach is applied to predict seed cotton yield phenotypes on the basis of gene expression. Top-ranking eGenes (NF-YB3, FLA2, and GRDP1) derived with pleiotropic effects on yield traits are validated, along with their potential roles by correlation analysis, domestication selection analysis, and transgenic plants.

Sequencing 4.3 million mutations in wheat promoters to understand and modify gene expression

Wheat is an important contributor to global food security, and further improvements are required to feed a growing human population. Functional genetics and genomics tools can help us to understand the function of different genes and to engineer beneficial changes. In this study, we used a promoter capture assay to sequence 2-kb regions upstream of all high-confidence annotated genes from 1,513 mutagenized plants from the tetraploid wheat variety Kronos. We identified 4.3 million induced mutations with an accuracy of 99.8%, resulting in a mutation density of 41.9 mutations per kb. We also remapped Kronos exome capture reads to Chinese Spring RefSeq v1.1, identified 4.7 million mutations, and predicted their effects on annotated genes. Using these predictions, we identified 59% more nonsynonymous substitutions and 49% more truncation mutations than in the original study. To show the biological value of the promoter dataset, we selected two mutations within the promoter of the VRN-A1 vernalization gene. Both mutations, located within transcription factor binding sites, significantly altered VRN-A1 expression, and one reduced the number of spikelets per spike. These publicly available sequenced mutant datasets provide rapid and inexpensive access to induced variation in the promoters and coding regions of most wheat genes. These mutations can be used to understand and modulate gene expression and phenotypes for both basic and commercial applications, where limited governmental regulations can facilitate deployment. These mutant collections, together with gene editing, provide valuable tools to accelerate functional genetic studies in this economically important crop.

Alterations in Growth Habit to Channel End-of-Season Perennial Reserves towards Increased Yield and Reduced Regrowth after Defoliation in Upland Cotton (Gossypium hirsutum L.)

Cotton (Gossypium spp.) is the primary source of natural textile fiber in the U.S. and a major crop in the Southeastern U.S. Despite constant efforts to increase the cotton fiber yield, the yield gain has stagnated. Therefore, we undertook a novel approach to improve the cotton fiber yield by altering its growth habit from perennial to annual. In this effort, we identified genotypes with high-expression alleles of five floral induction and meristem identity genes (FT, SOC1, FUL, LFY, and AP1) from an Upland cotton mini-core collection and crossed them in various combinations to develop cotton lines with annual growth habit, optimal flowering time, and enhanced productivity. To facilitate the characterization of genotypes with the desired combinations of stacked alleles, we identified molecular markers associated with the gene expression traits via genome-wide association analysis using a 63 K SNP Array. Over 14,500 SNPs showed polymorphism and were used for association analysis. A total of 396 markers showed associations with expression traits. Of these 396 markers, 159 were mapped to genes, 50 to untranslated regions, and 187 to random genomic regions. Biased genomic distribution of associated markers was observed where more trait-associated markers mapped to the cotton D sub-genome. Many quantitative trait loci coincided at specific genomic regions. This observation has implications as these traits could be bred together. The analysis also allowed the identification of candidate regulators of the expression patterns of these floral induction and meristem identity genes whose functions will be validated.

Duplicate Genes Contribute to Variability in Abiotic Stress Resistance in Allopolyploid Wheat

Gene duplication is a universal biological phenomenon that drives genomic variation and diversity, plays a crucial role in plant evolution, and contributes to innovations in genetic engineering and crop development. Duplicated genes participate in the emergence of novel functionality, such as adaptability to new or more severe abiotic stress resistance. Future crop research will benefit from advanced, mechanistic understanding of the effects of gene duplication, especially in the development and deployment of high-performance, stress-resistant, elite wheat lines. In this review, we summarize the current knowledge of gene duplication in wheat, including the principle of gene duplication and its effects on gene function, the diversity of duplicated genes, and how they have functionally diverged. Then, we discuss how duplicated genes contribute to abiotic stress response and the mechanisms of duplication. Finally, we have a future prospects section that discusses the direction of future efforts in the short term regarding the elucidation of replication and retention mechanisms of repetitive genes related to abiotic stress response in wheat, excellent gene function research, and practical applications.

The translational landscape of bread wheat during grain development

The dynamics of gene expression in crop grains has typically been investigated at the transcriptional level. However, this approach neglects translational regulation, a widespread mechanism that rapidly modulates gene expression to increase the plasticity of organisms. Here, we performed ribosome profiling and polysome profiling to obtain a comprehensive translatome data set of developing bread wheat (Triticum aestivum) grains. We further investigated the genome-wide translational dynamics during grain development, revealing that the translation of many functional genes is modulated in a stage-specific manner. The unbalanced translation between subgenomes is pervasive, which increases the expression flexibility of allohexaploid wheat. In addition, we uncovered widespread previously unannotated translation events, including upstream open reading frames (uORFs), downstream open reading frames (dORFs), and open reading frames (ORFs) in long noncoding RNAs, and characterized the temporal expression dynamics of small ORFs. We demonstrated that uORFs act as cis-regulatory elements that can repress or even enhance the translation of mRNAs. Gene translation may be combinatorially modulated by uORFs, dORFs, and microRNAs. In summary, our study presents a translatomic resource that provides a comprehensive and detailed overview of the translational regulation in developing bread wheat grains. This resource will facilitate future crop improvements for optimal yield and quality.

A Survey of the Transcriptomic Resources in Durum Wheat: Stress Responses, Data Integration and Exploitation

Durum wheat (Triticum turgidum subsp. durum (Desf.) Husn.) is an allotetraploid cereal crop of worldwide importance, given its use for making pasta, couscous, and bulgur. Under climate change scenarios, abiotic (e.g., high and low temperatures, salinity, drought) and biotic (mainly exemplified by fungal pathogens) stresses represent a significant limit for durum cultivation because they can severely affect yield and grain quality. The advent of next-generation sequencing technologies has brought a huge development in transcriptomic resources with many relevant datasets now available for durum wheat, at various anatomical levels, also focusing on phenological phases and environmental conditions. In this review, we cover all the transcriptomic resources generated on durum wheat to date and focus on the corresponding scientific insights gained into abiotic and biotic stress responses. We describe relevant databases, tools and approaches, including connections with other "omics" that could assist data integration for candidate gene discovery for bio-agronomical traits. The biological knowledge summarized here will ultimately help in accelerating durum wheat breeding.

A wheat integrative regulatory network from large-scale complementary functional datasets enables trait-associated gene discovery for crop improvement

Gene regulation is central to all aspects of organism growth, and understanding it using large-scale functional datasets can provide a whole view of biological processes controlling complex phenotypic traits in crops. However, the connection between massive functional datasets and trait-associated gene discovery for crop improvement is still lacking. In this study, we constructed a wheat integrative gene regulatory network (wGRN) by combining an updated genome annotation and diverse complementary functional datasets, including gene expression, sequence motif, transcription factor (TF) binding, chromatin accessibility, and evolutionarily conserved regulation. wGRN contains 7.2 million genome-wide interactions covering 5947 TFs and 127 439 target genes, which were further verified using known regulatory relationships, condition-specific expression, gene functional information, and experiments. We used wGRN to assign genome-wide genes to 3891 specific biological pathways and accurately prioritize candidate genes associated with complex phenotypic traits in genome-wide association studies. In addition, wGRN was used to enhance the interpretation of a spike temporal transcriptome dataset to construct high-resolution networks. We further unveiled novel regulators that enhance the power of spike phenotypic trait prediction using machine learning and contribute to the spike phenotypic differences among modern wheat accessions. Finally, we developed an interactive webserver, wGRN (http://wheat.cau.edu.cn/wGRN), for the community to explore gene regulation and discover trait-associated genes. Collectively, this community resource establishes the foundation for using large-scale functional datasets to guide trait-associated gene discovery for crop improvement.

High-quality haplotype-resolved genome assembly of cultivated octoploid strawberry

Cultivated strawberry (Fragaria × ananassa), a perennial herb belonging to the family Rosaceae, is a complex octoploid with high heterozygosity at most loci. However, there is no research on the haplotype of the octoploid strawberry genome. Here we aimed to obtain a high-quality genome of the cultivated strawberry cultivar, "Yanli", using single molecule real-time sequencing and high-throughput chromosome conformation capture technology. The "Yanli" genome was 823 Mb in size, with a long terminal repeat assembly index of 14.99. The genome was phased into two haplotypes, Hap1 (825 Mb with contig N50 of 26.70 Mb) and Hap2 (808 Mb with contig N50 of 27.51 Mb). Using the combination of Hap1 and Hap2, we obtained for the first time a haplotype-resolved genome with 56 chromosomes for the cultivated octoploid strawberry. We identified a ~ 10 Mb inversion and translocation on chromosome 2-1. 104 957 and 102 356 protein-coding genes were annotated in Hap1 and Hap2, respectively. Analysis of the genes related to the anthocyanin biosynthesis pathway revealed the structural diversity and complexity in the expression of the alleles in the octoploid F. × ananassa genome. In summary, we obtained a high-quality haplotype-resolved genome assembly of F. × ananassa, which will provide the foundation for investigating gene function and evolution of the genome of cultivated octoploid strawberry.

Multiplexed promoter and gene editing in wheat using a virus-based guide RNA delivery system

The low efficiency of genetic transformation and gene editing across diverse cultivars hinder the broad application of CRISPR technology for crop improvement. The development of virus-based methods of CRISPR-Cas system delivery into the plant cells holds great promise to overcome these limitations. Here, we perform direct inoculation of wheat leaves with the barley stripe mosaic virus (BSMV) transcripts to deliver guide RNAs (sgRNA) into the Cas9-expressing wheat. We demonstrate that wheat inoculation with the pool of BSMV-sgRNAs could be used to generate heritable precise deletions in the promoter region of a transcription factor and to perform multiplexed editing of agronomic genes. We transfer the high-expressing locus of Cas9 into adapted spring and winter cultivars by marker-assisted introgression and use of the BSMV-sgRNAs to edit two agronomic genes. A strategy presented in our study could be applied to any adapted cultivar for creating new cis-regulatory diversity or large-scale editing of multiple genes in biological pathways or QTL regions, opening possibilities for the effective engineering of crop genomes, and accelerating gene discovery and trait improvement efforts.

Comprehensive transcriptional variability analysis reveals gene networks regulating seed oil content of Brassica napus

BACKGROUND

Regulation of gene expression plays an essential role in controlling the phenotypes of plants. Brassica napus (B. napus) is an important source for the vegetable oil in the world, and the seed oil content is an important trait of B. napus.

RESULTS

We perform a comprehensive analysis of the transcriptional variability in the seeds of B. napus at two developmental stages, 20 and 40 days after flowering (DAF). We detect 53,759 and 53,550 independent expression quantitative trait loci (eQTLs) for 79,605 and 76,713 expressed genes at 20 and 40 DAF, respectively. Among them, the local eQTLs are mapped to the adjacent genes more frequently. The adjacent gene pairs are regulated by local eQTLs with the same open chromatin state and show a stronger mode of expression piggybacking. Inter-subgenomic analysis indicates that there is a feedback regulation for the homoeologous gene pairs to maintain partial expression dosage. We also identify 141 eQTL hotspots and find that hotspot87-88 co-localizes with a QTL for the seed oil content. To further resolve the regulatory network of this eQTL hotspot, we construct the XGBoost model using 856 RNA-seq datasets and the Basenji model using 59 ATAC-seq datasets. Using these two models, we predict the mechanisms affecting the seed oil content regulated by hotspot87-88 and experimentally validate that the transcription factors, NAC13 and SCL31, positively regulate the seed oil content.

CONCLUSIONS

We comprehensively characterize the gene regulatory features in the seeds of B. napus and reveal the gene networks regulating the seed oil content of B. napus.

A multi-omics framework reveals strawberry flavor genes and their regulatory elements

Flavor is essential to consumer preference of foods and is an increasing focus of plant breeding programs. In fruit crops, identifying genes underlying volatile organic compounds has great promise to accelerate flavor improvement, but polyploidy and heterozygosity in many species have slowed progress. Here we use octoploid cultivated strawberry to demonstrate how genomic heterozygosity, transcriptomic intricacy and fruit metabolomic diversity can be treated as strengths and leveraged to uncover fruit flavor genes and their regulatory elements. Multi-omics datasets were generated including an expression quantitative trait loci map with 196 diverse breeding lines, haplotype-phased genomes of a highly-flavored breeding selection, a genome-wide structural variant map using five haplotypes, and volatile genome-wide association study (GWAS) with > 300 individuals. Overlaying regulatory elements, structural variants and GWAS-linked allele-specific expression of numerous genes to variation in volatile compounds important to flavor. In one example, the functional role of anthranilate synthase alpha subunit 1 in methyl anthranilate biosynthesis was supported via fruit transient gene expression assays. These results demonstrate a framework for flavor gene discovery in fruit crops and a pathway to molecular breeding of cultivars with complex and desirable flavor.

Robustness and the generalist niche of polyploid species: Genome shock or gradual evolution?

The prevalence of polyploidy in wild and crop species has stimulated debate over its evolutionary advantages and disadvantages. Previous studies have focused on changes occurring at the polyploidization events, including genome-wide changes termed "genome shock," as well as ancient polyploidy. Recent bioinformatics advances and empirical studies of Arabidopsis and wheat relatives are filling a research gap: the functional evolutionary study of polyploid species using RNA-seq, DNA polymorphism, and epigenomics. Polyploid species can become generalists in natura through environmental robustness by inheriting and merging parental stress responses. Their evolvability is enhanced by mutational robustness working on inherited standing variation. The identification of key genes responsible for gradual adaptive evolution will encourage synthetic biological approaches to transfer polyploid advantages to other species.

Transcriptional regulation of photomorphogenesis in seedlings of Brassica napus under different light qualities

This study displayed the transcriptional regulation network of key regulators and downstream pathway in seedling morphogenesis of Brassica napus under different light quality. Plants undergo photomorphogenesis upon the presence of light, mediated by different light (e.g., blue, red, and far-red) signaling pathways. Although the light signaling pathway has been well documented in Arabidopsis, the underlying mechanisms were studied to a less extent in other plant species including Brassica napus. In this study, we investigated the effect of different light qualities (white, blue, red, and far-red light) on the hypocotyl elongation in B. napus, and performed the transcriptomic analysis of seedlings in response to different light qualities. The results showed that hypocotyl elongation was slightly inhibited by red light, while it was strongly inhibited by blue/far-red light. Transcriptome analysis identified 9748 differentially expressed genes (DEGs) among treatments. Gene ontology (GO) enrichment analysis of DEGs showed that light-responsive and photosynthesis-related genes were highly expressed in response to blue/far-red light rather than in red light. Furthermore, the key genes in light signaling (i.e., PHYB, HY5, HYH, HFR1, and PIF3) exhibited distinct expression patterns between blue/far-red and red light treatments. In addition, subgenome dominant expression of homoeologous genes were observed for some genes, such as PHYA, PHYB, HFR1, and BBXs. The current study displayed a comprehensive dissection of light-mediated transcriptional regulation network, including light signaling, phytohormone, and cell elongation/modification, which improved the understanding on the underlying mechanism of light-regulated hypocotyl growth in B. napus.

TCP Transcription Factors Involved in Shoot Development of Ma Bamboo (Dendrocalamus latiflorus Munro)

Ma bamboo (Dendrocalamus latiflorus Munro) is the most widely cultivated clumping bamboo in Southern China and is valuable for both consumption and wood production. The development of bamboo shoots involving the occurrence of lateral buds is unique, and it affects both shoot yield and the resulting timber. Plant-specific TCP transcription factors are involved in plant growth and development, particularly in lateral bud outgrowth and morphogenesis. However, the comprehensive information of the TCP genes in Ma bamboo remains poorly understood. In this study, 66 TCP transcription factors were identified in Ma bamboo at the genome-wide level. Members of the same subfamily had conservative gene structures and conserved motifs. The collinear analysis demonstrated that segmental duplication occurred widely in the TCP transcription factors of Ma bamboo, which mainly led to the expansion of a gene family. Cis-acting elements related to growth and development and stress response were found in the promoter regions of DlTCPs. Expression patterns revealed that DlTCPs have tissue expression specificity, which is usually highly expressed in shoots and leaves. Subcellular localization and transcriptional self-activation experiments demonstrated that the five candidate TCP proteins were typical self-activating nuclear-localized transcription factors. Additionally, the transcriptome analysis of the bamboo shoot buds at different developmental stages helped to clarify the underlying functions of the TCP members during the growth of bamboo shoots. DlTCP12-C, significantly downregulated as the bamboo shoots developed, was selected to further verify its molecular function in Arabidopsis. The DlTCP12-C overexpressing lines exhibited a marked reduction in the number of rosettes and branches compared with the wild type in Arabidopsis, suggesting that DlTCP12-C conservatively inhibits lateral bud outgrowth and branching in plants. This study provides useful insights into the evolutionary patterns and molecular functions of the TCP transcription factors in Ma bamboo and provides a valuable reference for further research on the regulatory mechanism of bamboo shoot development and lateral bud growth.