Insights into transcriptional characteristics and homoeolog expression bias of embryo and de-embryonated kernels in developing grain through RNA-Seq and Iso-Seq

Evolutionary trajectory of organelle-derived nuclear DNAs in the Triticum/Aegilops complex species

Organelle-derived nuclear DNAs, nuclear plastid DNAs (NUPTs), and nuclear mitochondrial DNAs (NUMTs) have been identified in plants. Most, if not all, genes residing in NUPTs/NUMTs (NUPGs/NUMGs) are known to be inactivated and pseudogenized. However, the role of epigenetic control in silencing NUPGs/NUMGs and the dynamic evolution of NUPTs/NUMTs with respect to organismal phylogeny remain barely explored. Based on the available nuclear and organellar genomic resources of wheat (genus Triticum) and goat grass (genus Aegilops) within Triticum/Aegilops complex species, we investigated the evolutionary fates of NUPTs/NUMTs in terms of their epigenetic silencing and their dynamic occurrence rates in the nuclear diploid genomes and allopolyploid subgenomes. NUPTs and NUMTs possessed similar genomic atlas, including (i) predominantly located in intergenic regions and preferential integration to gene regulation regions and (ii) generating sequence variations in the nuclear genome. Unlike nuclear indigenous genes, the alien NUPGs/NUMGs were associated with repressive epigenetic signals, namely high levels of DNA methylation and low levels of active histone modifications. Phylogenomic analyses suggested that the species-specific and gradual accumulation of NUPTs/NUMTs accompanied the speciation processes. Moreover, based on further pan-genomic analyses, we found significant subgenomic asymmetry in the NUPT/NUMT occurrence, which accumulated during allopolyploid wheat evolution. Our findings provide insight into the dynamic evolutionary fates of organelle-derived nuclear DNA in plants.

Kernel Transcriptome Profiles of Susceptible Wheat Genotypes in Response to Wheat Dwarf Bunt

Wheat dwarf bunt is caused by Tilletia controversa J. G. Kühn (TCK), which is a serious fungal diseases affecting kernels of wheat. In order to identify candidate genes involved in the abnormal development of kernels in wheat, we used RNA sequencing technology to analyze the transcriptome of the abnormal and healthy kernels of a susceptible variety (Yili053) at the mid-filling stage, late-filling stage, and maturity stage, respectively. The differentially expressed genes (DEGs) were analyzed, and there were 3930 DEGs, 28,422 DEGs, and 20,874 DEGs found at the mid-filling stage, late-filling stage, and maturity stage in Yili053, respectively. A total of 1592 DEGs (506 DEGs up-regulated) showed continuously differential expression in the three stages. Gene ontology analysis showed that these DEGs were related to biological regulation, metabolic processes, and the response to stimulus. Kyoto Encyclopedia of Genes and Genomes enrichment analysis showed that these DEGs play major roles in pathways including photosynthesis, carbon metabolism, carbon fixation in photosynthetic organisms, and glyoxylate and dicarboxylate metabolism. Moreover, we predicted that 13 MADS-MIKC transcription factors, which were continuously up-regulated, were crucial for regulating the maturation and senescence of eukaryotes. Some 21 genes related to the plant hormone signaling transduction pathway and 61 genes related to the response to stimulus were analyzed. A total of 26 of them were successful validated with a qPCR analysis. These genes were thought to be involved in the abnormal development of kernels infected by TCK. A transcriptomics analysis of wheat kernels in response to TCK will contribute to understanding the interaction of TCK and wheat, and may provide a basis for knowledge of molecular events in the abnormal development of kernels, which will be helpful for more efficient TCK management.

Cold stress induces differential gene expression of retained homeologs in Camelina sativa cv Suneson

Camelina sativa (L.) Crantz, a member of the Brassicaceae, has potential as a biofuel feedstock which is attributable to the production of fatty acids in its seeds, its fast growth cycle, and low input requirements. While a genome assembly is available for camelina, it was generated from short sequence reads and is thus highly fragmented in nature. Using long read sequences, we generated a chromosome-scale, highly contiguous genome assembly (644,491,969 bp) for the spring biotype cultivar 'Suneson' with an N50 contig length of 12,031,512 bp and a scaffold N50 length of 32,184,682 bp. Annotation of protein-coding genes revealed 91,877 genes that encode 133,355 gene models. We identified a total of 4,467 genes that were significantly up-regulated under cold stress which were enriched in gene ontology terms associated with "response to cold" and "response to abiotic stress". Coexpression analyses revealed multiple coexpression modules that were enriched in genes differentially expressed following cold stress that had putative functions involved in stress adaptation, specifically within the plastid. With access to a highly contiguous genome assembly, comparative analyses with Arabidopsis thaliana revealed 23,625 A. thaliana genes syntenic with 45,453 Suneson genes. Of these, 24,960 Suneson genes were syntenic to 8,320 A. thaliana genes reflecting a 3 camelina homeolog to 1 Arabidopsis gene relationship and retention of all three homeologs. Some of the retained triplicated homeologs showed conserved gene expression patterns under control and cold-stressed conditions whereas other triplicated homeologs displayed diverged expression patterns revealing sub- and neo-functionalization of the homeologs at the transcription level. Access to the chromosome-scale assembly of Suneson will enable both basic and applied research efforts in the improvement of camelina as a sustainable biofuel feedstock.

Alternative Splicing Variation: Accessing and Exploiting in Crop Improvement Programs

Alternative splicing (AS) is a gene regulatory mechanism modulating gene expression in multiple ways. AS is prevalent in all eukaryotes including plants. AS generates two or more mRNAs from the precursor mRNA (pre-mRNA) to regulate transcriptome complexity and proteome diversity. Advances in next-generation sequencing, omics technology, bioinformatics tools, and computational methods provide new opportunities to quantify and visualize AS-based quantitative trait variation associated with plant growth, development, reproduction, and stress tolerance. Domestication, polyploidization, and environmental perturbation may evolve novel splicing variants associated with agronomically beneficial traits. To date, pre-mRNAs from many genes are spliced into multiple transcripts that cause phenotypic variation for complex traits, both in model plant Arabidopsis and field crops. Cataloguing and exploiting such variation may provide new paths to enhance climate resilience, resource-use efficiency, productivity, and nutritional quality of staple food crops. This review provides insights into AS variation alongside a gene expression analysis to select for novel phenotypic diversity for use in breeding programs. AS contributes to heterosis, enhances plant symbiosis (mycorrhiza and rhizobium), and provides a mechanistic link between the core clock genes and diverse environmental clues.

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.

Molecular prospective on the wheat grain development

Wheat grain development is an important biological process to determine grain yield and quality, which is controlled by the interplay of genetic, epigenetic, and environmental factors. Wheat grain development has been extensively characterized at the phenotypic and genetic levels. The advent of innovative molecular technologies allows us to characterize genes, proteins, and regulatory factors involved in wheat grain development, which have enhanced our understanding of the wheat seed development process. However, wheat is an allohexaploid with a large genome size, the molecular mechanisms underlying the wheat grain development have not been well understood as those in diploids. Understanding grain development, and how it is regulated, is of fundamental importance for improving grain yield and quality through conventional breeding or genetic engineering. Herein, we review the current discoveries on the molecular mechanisms underlying wheat grain development. Notably, only a handful of genes that control wheat grain development have, thus far, been well characterized, their interplay underlying the grain development remains elusive. The synergistic network-integrated genomics and epigenetics underlying wheat grain development and how the subgenome divergence dynamically and precisely regulates wheat grain development are unknown.

Genome-wide association study reveals a NAC transcription factor TaNAC074 linked to pre-harvest sprouting tolerance in wheat

Twelve QTL associated with pre-harvest sprouting tolerance were identified using association analysis in wheat. Two markers were validated and a candidate gene TaNAC074 for Qgpf.cas-3B.2 was verified using Agrobacterium-mediated transformation. Pre-harvest sprouting (PHS) is a considerable global threat to wheat yield and quality. Due to this threat, breeders must identify quantitative trait loci (QTL) and genes conferring PHS-tolerance (PHST) to reduce the negative effects of PHS caused by low seed dormancy. In this study, we evaluated a panel of 302 diverse wheat genotypes for PHST in four environments and genotyped the panel with a high-density wheat 660 K SNP array. By using a genome-wide association study (GWAS), we identified 12 stable loci significantly associated with PHST (P < 0.0001), explaining 3.34 - 9.88% of the phenotypic variances. Seven of these loci co-located with QTL and genes reported previously. Five loci (Qgpf.cas-3B.2, Qgpf.cas-3B.3, Qgpf.cas-3B.4, Qgpf.cas-7B.2, and Qgpf.cas-7B.3), located in genomic regions with no known PHST QTL or genes, are likely to be new QTL conferring PHST. Additionally, two molecular markers were developed for Qgpf.cas-3A and Qgpf.cas-7B.3, and validated using a different set of 233 wheat accessions. Finally, the PHST-related function of candidate gene TaNAC074 for Qgpf.cas-3B.2 was confirmed by CAPS (cleaved amplified polymorphic sequences) marker association analysis in 233 wheat accessions and by expression and phenotypic analysis of transgenic wheat. Overexpression of TaNAC074 significantly reduced seed dormancy in wheat. This study contributes to broaden the genetic basis and molecular marker-assisted breeding of PHST.