Identification and characterization of shared duplications between rice and wheat provide new insight into grass genome evolution

Integrative gene duplication and genome-wide analysis as an approach to facilitate wheat reverse genetics: An example in the TaCIPK family

INTRODUCTION

Reverse genetic studies conducted in the plant with a complex or polyploidy genome enriched with large gene families (like wheat) often meet challenges in identifying the key candidate genes related to important traits and prioritizing the genes for functional experiments.

OBJECTIVE

To overcome the above-mentioned challenges of reverse genetics, this work aims to establish an efficient multi-species strategy for genome-wide gene identification and prioritization of the key candidate genes.

METHODS

We established the integrative gene duplication and genome-wide analysis (iGG analysis) as a strategy for pinpointing key candidate genes deserving functional research. The iGG captures the evolution, and the expansion/contraction of large gene families across phylogeny-related species and integrates spatial-temporal expression information for gene function inference. Transgenic approaches were also employed to functional validation.

RESULTS

As a proof-of-concept for the iGG analysis, we took the wheat calcineurin B-like protein-interacting protein kinases (CIPKs) family as an example. We identified CIPKs from seven monocot species, established the orthologous relationship of CIPKs between rice and wheat, and characterized Triticeae-specific CIPK duplicates (e.g., CIPK4 and CIPK17). Integrated with our analysis of CBLs and CBL-CIPK interaction, we revealed that divergent expressions of TaCBLs and TaCIPKs could play an important role in keeping the stoichiometric balance of CBL-CIPK. Furthermore, we validated the function of TaCIPK17-A2 in the regulation of drought tolerance by using transgenic approaches. Overexpression of TaCIPK17 enhanced antioxidant capacity and improved drought tolerance in wheat.

CONCLUSION

The iGG analysis leverages evolutionary and comparative genomics of crops with large genomes to rapidly highlight the duplicated genes potentially associated with speciation, domestication and/or particular traits that deserve reverse-genetic functional studies. Through the identification of Triticeae-specific TaCIPK17 duplicates and functional validation, we demonstrated the effectiveness of the iGG analysis and provided a new target gene for improving drought tolerance in wheat.

Patterns of the Predicted Mutation Burden in 19,778 Domesticated Barley Accessions Conserved Ex Situ

Long-term conservation of more than 7 million plant germplasm accessions in 1750 genebanks worldwide is a challenging mission. The extent of deleterious mutations present in conserved germplasm and the genetic risk associated with accumulative mutations are largely unknown. This study took advantage of published barley genomic data to predict sample-wise mutation burdens for 19,778 domesticated barley (Hordeum vulgare L.) accessions conserved ex situ. It was found that the conserved germplasm harbored 407 deleterious mutations and 337 (or 82%) identified deleterious alleles were present in 20 (or 0.1%) or fewer barley accessions. Analysis of the predicted mutation burdens revealed significant differences in mutation burden for several groups of barley germplasm (landrace > cultivar (or higher burden estimate in landrace than in cultivar); winter barley > spring barley; six-rowed barley > two-rowed barley; and 1000-accession core collection > non-core germplasm). Significant differences in burden estimate were also found among seven major geographical regions. The sample-wise predicted mutation burdens were positively correlated with the estimates of sample average pairwise genetic difference. These findings are significant for barley germplasm management and utilization and for a better understanding of the genetic risk in conserved plant germplasm.

The additive function of YIGE2 and YIGE1 in regulating maize ear length

Ear length (EL) is a key trait that greatly contributes to yield in maize. Although dozens of EL quantitative trait loci have been mapped, very few causal genes have been cloned, and the molecular mechanisms remain largely unknown. Our previous study showed that YIGE1 is involved in sugar and auxin pathways to regulate ear inflorescence meristem (IM) development and thus affects EL in maize. Here, we reveal that YIGE2, the paralog of YIGE1, regulates maize ear development and EL through auxin pathway. Knockout of YIGE2 causes a significant decrease of auxin level, IM length, floret number, EL, and grain yield. yige1 yige2 double mutants had even shorter IM and ears implying that these two genes redundantly regulate IM development and EL. The genes controlling auxin levels are differential expressed in yige1 yige2 double mutants, leading to lower auxin level. These results elucidated the critical role of YIGE2 and the redundancy between YIGE2 and YIGE1 in maize ear development, providing a new genetic resource for maize yield improvement.

Genome-Wide Identification and Expression Analysis of Catalase Gene Families in Triticeae

Aerobic metabolism in plants results in the production of hydrogen peroxide (H2O2), a significant and comparatively stable non-radical reactive oxygen species (ROS). H2O2 is a signaling molecule that regulates particular physiological and biological processes (the cell cycle, photosynthesis, plant growth and development, and plant responses to environmental challenges) at low concentrations. Plants may experience oxidative stress and ultimately die from cell death if excess H2O2 builds up. Triticum dicoccoides, Triticum urartu, and Triticum spelta are different ancient wheat species that present different interesting characteristics, and their importance is becoming more and more clear. In fact, due to their interesting nutritive health, flavor, and nutritional values, as well as their resistance to different parasites, the cultivation of these species is increasingly important. Thus, it is important to understand the mechanisms of plant tolerance to different biotic and abiotic stresses by studying different stress-induced gene families such as catalases (CAT), which are important H2O2-metabolizing enzymes found in plants. Here, we identified seven CAT-encoding genes (TdCATs) in Triticum dicoccoides, four genes in Triticum urartu (TuCATs), and eight genes in Triticum spelta (TsCATs). The accuracy of the newly identified wheat CAT gene members in different wheat genomes is confirmed by the gene structures, phylogenetic relationships, protein domains, and subcellular location analyses discussed in this article. In fact, our analysis showed that the identified genes harbor the following two conserved domains: a catalase domain (pfam00199) and a catalase-related domain (pfam06628). Phylogenetic analyses showed that the identified wheat CAT proteins were present in an analogous form in durum wheat and bread wheat. Moreover, the identified CAT proteins were located essentially in the peroxisome, as revealed by in silico analyses. Interestingly, analyses of CAT promoters in those species revealed the presence of different cis elements related to plant development, maturation, and plant responses to different environmental stresses. According to RT-qPCR, Triticum CAT genes showed distinctive expression designs in the studied organs and in response to different treatments (salt, heat, cold, mannitol, and ABA). This study completed a thorough analysis of the CAT genes in Triticeae, which advances our knowledge of CAT genes and establishes a framework for further functional analyses of the wheat gene family.

Subgenomic Stability of Progenitor Genomes During Repeated Allotetraploid Origins of the Same Grass Brachypodium hybridum

Both homeologous exchanges and homeologous expression bias are generally found in most allopolyploid species. Whether homeologous exchanges and homeologous expression bias differ between repeated allopolyploid speciation events from the same progenitor species remains unknown. Here, we detected a third independent and recent allotetraploid origin for the model grass Brachypodium hybridum. Our homeologous exchange with replacement analyses indicated the absence of significant homeologous exchanges in any of the three types of wild allotetraploids, supporting the integrity of their progenitor subgenomes and the immediate creation of the amphidiploids. Further homeologous expression bias tests did not uncover significant subgenomic dominance in different tissues and conditions of the allotetraploids. This suggests a balanced expression of homeologs under similar or dissimilar ecological conditions in their natural habitats. We observed that the density of transposons around genes was not associated with the initial establishment of subgenome dominance; rather, this feature is inherited from the progenitor genome. We found that drought response genes were highly induced in the two subgenomes, likely contributing to the local adaptation of this species to arid habitats in the third allotetraploid event. These findings provide evidence for the consistency of subgenomic stability of parental genomes across multiple allopolyploidization events that led to the same species at different periods. Our study emphasizes the importance of selecting closely related progenitor species genomes to accurately assess homeologous exchange with replacement in allopolyploids, thereby avoiding the detection of false homeologous exchanges when using less related progenitor species genomes.

Genetic Basis of Grain Size and Weight in Rice, Wheat, and Barley

Grain size is a key component of grain yield in cereals. It is a complex quantitative trait controlled by multiple genes. Grain size is determined via several factors in different plant development stages, beginning with early tillering, spikelet formation, and assimilates accumulation during the pre-anthesis phase, up to grain filling and maturation. Understanding the genetic and molecular mechanisms that control grain size is a prerequisite for improving grain yield potential. The last decade has brought significant progress in genomic studies of grain size control. Several genes underlying grain size and weight were identified and characterized in rice, which is a model plant for cereal crops. A molecular function analysis revealed most genes are involved in different cell signaling pathways, including phytohormone signaling, transcriptional regulation, ubiquitin-proteasome pathway, and other physiological processes. Compared to rice, the genetic background of grain size in other important cereal crops, such as wheat and barley, remains largely unexplored. However, the high level of conservation of genomic structure and sequences between closely related cereal crops should facilitate the identification of functional orthologs in other species. This review provides a comprehensive overview of the genetic and molecular bases of grain size and weight in wheat, barley, and rice, focusing on the latest discoveries in the field. We also present possibly the most updated list of experimentally validated genes that have a strong effect on grain size and discuss their molecular function.

Exploring the genetic landscape of nitrogen uptake in durum wheat: genome-wide characterization and expression profiling of NPF and NRT2 gene families

Nitrate uptake by plants primarily relies on two gene families: Nitrate transporter 1/peptide transporter (NPF) and Nitrate transporter 2 (NRT2). Here, we extensively characterized the NPF and NRT2 families in the durum wheat genome, revealing 211 NPF and 20 NRT2 genes. The two families share many Cis Regulatory Elements (CREs) and Transcription Factor binding sites, highlighting a partially overlapping regulatory system and suggesting a coordinated response for nitrate transport and utilization. Analyzing RNA-seq data from 9 tissues and 20 cultivars, we explored expression profiles and co-expression relationships of both gene families. We observed a strong correlation between nucleotide variation and gene expression within the NRT2 gene family, implicating a shared selection mechanism operating on both coding and regulatory regions. Furthermore, NPF genes showed highly tissue-specific expression profiles, while NRT2s were mainly divided in two co-expression modules, one expressed in roots (NAR2/NRT3 dependent) and the other induced in anthers and/ovaries during maturation. Our evidences confirmed that the majority of these genes were retained after small-scale duplication events, suggesting a neo- or sub-functionalization of many NPFs and NRT2s. Altogether, these findings indicate that the expansion of these gene families in durum wheat could provide valuable genetic variability useful to identify NUE-related and candidate genes for future breeding programs in the context of low-impact and sustainable agriculture.

GhCKX14 responding to drought stress by modulating antioxi-dative enzyme activity in Gossypium hirsutum compared to CKX family genes

BACKGROUND

Cytokinin oxidase/dehydrogenase (CKX) plays a vital role in response to abiotic stress through modulating the antioxidant enzyme activities. Nevertheless, the biological function of the CKX gene family has yet to be reported in cotton.

RESULT

In this study, a total of 27 GhCKXs were identified by the genome-wide investigation and distributed across 18 chromosomes. Phylogenetic tree analysis revealed that CKX genes were clustered into four clades, and most gene expansions originated from segmental duplications. The CKXs gene structure and motif analysis displayed remarkably well conserved among the four groups. Moreover, the cis-acting elements related to the abiotic stress, hormones, and light response were identified within the promoter regions of GhCKXs. Transcriptome data and RT-qPCR showed that GhCKX genes demonstrated higher expression levels in various tissues and were involved in cotton's abiotic stress and phytohormone response. The protein-protein interaction network indicates that the CKX family probably participated in redox regulation, including oxidoreduction or ATP levels, to mediate plant growth and development. Functionally identified via virus-induced gene silencing (VIGS) found that the GhCKX14 gene improved drought resistance by modulating the antioxidant-related activitie.

CONCLUSIONS

In this study, the CKX gene family members were analyzed by bioinformatics, and validates the response of GhCKX gene to various phytohormone treatment and abiotic stresses. Our findings established the foundation of GhCKXs in responding to abiotic stress and GhCKX14 in regulating drought resistance in cotton.

Genome-wide characterization of the VQ genes in Triticeae and their functionalization driven by polyploidization and gene duplication events in wheat

Valine-glutamine motif-containing (VQ) proteins are transcriptional cofactors widely involved in plant growth, development, and response to various stresses. Although the VQ family has been genome-wide identified in some species, but the knowledge regarding duplication-driven functionalization of VQ genes among evolutionarily related species is still lacking. Here, 952 VQ genes have been identified from 16 species, emphasizing seven Triticeae species including the bread wheat. Comprehensive phylogenetic and syntenic analyses allow us to establish the orthologous relationship of VQ genes from rice (Oryza sativa) to bread wheat (Triticum aestivum). The evolutionary analysis revealed that whole-genome duplication (WGD) drives the expansion of OsVQs, while TaVQs expansion is associated with a recent burst of gene duplication (RBGD). We also analyzed the motif composition and molecular properties of TaVQ proteins, enriched biological functions, and expression patterns of TaVQs. We demonstrate that WGD-derived TaVQs have become divergent in both protein motif composition and expression pattern, while RBGD-derived TaVQs tend to adopt specific expression patterns, suggesting their functionalization in certain biological processes or in response to specific stresses. Furthermore, some RBGD-derived TaVQs are found to be associated with salt tolerance. Several of the identified salt-related TaVQ proteins were located in the cytoplasm and nucleus and their salt-responsive expression patterns were validated by qPCR analysis. Yeast-based functional experiments confirmed that TaVQ27 may be a new regulator to salt response and regulation. Overall, this study lays the foundation for further functional validation of VQ family members within the Triticeae species.

A multi-omics integrative network map of maize

Networks are powerful tools to uncover functional roles of genes in phenotypic variation at a system-wide scale. Here, we constructed a maize network map that contains the genomic, transcriptomic, translatomic and proteomic networks across maize development. This map comprises over 2.8 million edges in more than 1,400 functional subnetworks, demonstrating an extensive network divergence of duplicated genes. We applied this map to identify factors regulating flowering time and identified 2,651 genes enriched in eight subnetworks. We validated the functions of 20 genes, including 18 with previously unknown connections to flowering time in maize. Furthermore, we uncovered a flowering pathway involving histone modification. The multi-omics integrative network map illustrates the principles of how molecular networks connect different types of genes and potential pathways to map a genome-wide functional landscape in maize, which should be applicable in a wide range of species.

Breeding for Higher Yields of Wheat and Rice through Modifying Nitrogen Metabolism

Wheat and rice produce nutritious grains that provide 32% of the protein in the human diet globally. Here, we examine how genetic modifications to improve assimilation of the inorganic nitrogen forms ammonium and nitrate into protein influence grain yield of these crops. Successful breeding for modified nitrogen metabolism has focused on genes that coordinate nitrogen and carbon metabolism, including those that regulate tillering, heading date, and ammonium assimilation. Gaps in our current understanding include (1) species differences among candidate genes in nitrogen metabolism pathways, (2) the extent to which relative abundance of these nitrogen forms across natural soil environments shape crop responses, and (3) natural variation and genetic architecture of nitrogen-mediated yield improvement. Despite extensive research on the genetics of nitrogen metabolism since the rise of synthetic fertilizers, only a few projects targeting nitrogen pathways have resulted in development of cultivars with higher yields. To continue improving grain yield and quality, breeding strategies need to focus concurrently on both carbon and nitrogen assimilation and consider manipulating genes with smaller effects or that underlie regulatory networks as well as genes directly associated with nitrogen metabolism.

WGDI: A user-friendly toolkit for evolutionary analyses of whole-genome duplications and ancestral karyotypes

Evidence of whole-genome duplications (WGDs) and subsequent karyotype changes has been detected in most major lineages of living organisms on Earth. To clarify the complex resulting multi-layered patterns of gene collinearity in genome analyses, there is a need for convenient and accurate toolkits. To meet this need, we developed WGDI (Whole-Genome Duplication Integrated analysis), a Python-based command-line tool that facilitates comprehensive analysis of recursive polyploidization events and cross-species genome alignments. WGDI supports three main workflows (polyploid inference, hierarchical inference of genomic homology, and ancestral chromosome karyotyping) that can improve the detection of WGD and characterization of WGD-related events based on high-quality chromosome-level genomes. Significantly, it can extract complete synteny blocks and facilitate reconstruction of detailed karyotype evolution. This toolkit is freely available at GitHub (https://github.com/SunPengChuan/wgdi). As an example of its application, WGDI convincingly clarified karyotype evolution in Aquilegia coerulea and Vitis vinifera following WGDs and rejected the hypothesis that Aquilegia contributed as a parental lineage to the allopolyploid origin of core dicots.

MIKC-type MADS-box transcription factor gene family in peanut: Genome-wide characterization and expression analysis under abiotic stress

Peanut (Arachis hypogaea) is one of the most important economic crops around the world, especially since it provides vegetable oil and high-quality protein for humans. Proteins encoded by MADS-box transcription factors are widely involved in regulating plant growth and development as well as responses to abiotic stresses. However, the MIKC-type MADS-box TFs in peanut remains currently unclear. Hence, in this study, 166 MIKC-type MADS-box genes were identified in both cultivated and wild-type peanut genomes, which were divided into 12 subfamilies. We found a variety of development-, hormone-, and stress-related cis-acting elements in the promoter region of peanut MIKC-type MADS-box genes. The chromosomal distribution of peanut MADS-box genes was not random, and gene duplication contributed to the expansion of the MADS-box gene family. The interaction network of the peanut AhMADS proteins was established. Expression pattern analysis showed that AhMADS genes were specifically expressed in tissues and under abiotic stresses. It was further confirmed via the qRT-PCR technique that five selected AhMADS genes could be induced by abiotic and hormone treatments and presented different expressive profiles under various stresses. Taken together, these findings provide valuable information for the exploration of candidate genes in molecular breeding and further study of AhMADS gene functions.

An alternative angiosperm DGAT1 topology and potential motifs in the N-terminus

The highly variable cytoplasmic N-terminus of the plant diacylglycerol acyltransferase 1 (DGAT1) has been shown to have roles in oligomerization as well as allostery; however, the biological significance of the variation within this region is not understood. Comparing the coding sequences over the variable N-termini revealed the Poaceae DGAT1s contain relatively high GC compositional gradients as well as numerous direct and inverted repeats in this region. Using a variety of reciprocal chimeric DGAT1s from angiosperms we show that related N-termini had similar effects (positive or negative) on the accumulation of the recombinant protein in Saccharomyces cerevisiae. When expressed in Camelina sativa seeds the recombinant proteins of specific chimeras elevated total lipid content of the seeds as well as increased seed size. In addition, we combine N- and C-terminal as well as internal tags with high pH membrane reformation, protease protection and differential permeabilization. This led us to conclude the C-terminus is in the ER lumen; this contradicts earlier reports of the cytoplasmic location of plant DGAT1 C-termini.

Phylogenomics and Systematics of Overlooked Mesoamerican and South American Polyploid Broad-Leaved Festuca Grasses Differentiate F. sects. Glabricarpae and Ruprechtia and F. subgen. Asperifolia, Erosiflorae, Mallopetalon and Coironhuecu (subgen. nov.)

Allopolyploidy is considered a driver of diversity in subtribe Loliinae. We investigate the evolution and systematics of the poorly studied Mesoamerican and South American polyploid broad-leaved Festuca L. species of uncertain origin and unclear taxonomy. A taxonomic study of seven diagnostic morphological traits was conducted on a representation of 22 species. Phylogenomic analyses were performed on a representation of these supraspecific taxa and all other Loliinae lineages using separate data from the entire plastome, nuclear rDNA 45S and 5S genes, and repetitive DNA elements. F. subgen. Mallopetalon falls within the fine-leaved (FL) Loliinae clade, whereas the remaining taxa are nested within the broad-leaved (BL) Loliinae clade forming two separate Mexico–Central–South American (MCSAI, MCSAII) lineages. MCSAI includes representatives of F. sect. Glabricarpae and F. subgen. Asperifolia plus F. superba, and MCSAII of F. subgen. Erosiflorae and F. sect. Ruprechtia plus F. argentina. MCSAII likely had a BL Leucopoa paternal ancestor, MCSAI and MCSAII a BL Meso-South American maternal ancestor, and Mallopetalon FL, American I–II ancestors. Plastome vs. nuclear topological discordances corroborated the hybrid allopolyploid origins of these taxa, some of which probably originated from Northern Hemisphere ancestors. The observed data indicate rapid reticulate radiations in the Central–South American subcontinent. Our systematic study supports the reclassification of some studied taxa in different supraspecific Festuca ranks.

Wheat genomic study for genetic improvement of traits in China

Bread wheat (Triticum aestivum L.) is a major crop that feeds 40% of the world's population. Over the past several decades, advances in genomics have led to tremendous achievements in understanding the origin and domestication of wheat, and the genetic basis of agronomically important traits, which promote the breeding of elite varieties. In this review, we focus on progress that has been made in genomic research and genetic improvement of traits such as grain yield, end-use traits, flowering regulation, nutrient use efficiency, and biotic and abiotic stress responses, and various breeding strategies that contributed mainly by Chinese scientists. Functional genomic research in wheat is entering a new era with the availability of multiple reference wheat genome assemblies and the development of cutting-edge technologies such as precise genome editing tools, high-throughput phenotyping platforms, sequencing-based cloning strategies, high-efficiency genetic transformation systems, and speed-breeding facilities. These insights will further extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture in China and throughout the world.

The Original Form of C4-Photosynthetic Phosphoenolpyruvate Carboxylase Is Retained in Pooids but Lost in Rice

Poaceae is the most prominent monocot family that contains the primary cereal crops wheat, rice, and maize. These cereal species exhibit physiological diversity, such as different photosynthetic systems and environmental stress tolerance. Phosphoenolpyruvate carboxylase (PEPC) in Poaceae is encoded by a small multigene family and plays a central role in C4-photosynthesis and dicarboxylic acid metabolism. Here, to better understand the molecular basis of the cereal species diversity, we analyzed the PEPC gene family in wheat together with other grass species. We could designate seven plant-type and one bacterial-type grass PEPC groups, ppc1a, ppc1b, ppc2a, ppc2b, ppc3, ppc4, ppcC4, and ppc-b, respectively, among which ppc1b is an uncharacterized type of PEPC. Evolutionary inference revealed that these PEPCs were derived from five types of ancient PEPCs (ppc1, ppc2, ppc3, ppc4, and ppc-b) in three chromosomal blocks of the ancestral Poaceae genome. C4-photosynthetic PEPC (ppcC4 ) had evolved from ppc1b, which seemed to be arisen by a chromosomal duplication event. We observed that ppc1b was lost in many Oryza species but preserved in Pooideae after natural selection. In silico analysis of cereal RNA-Seq data highlighted the preferential expression of ppc1b in upper ground organs, selective up-regulation of ppc1b under osmotic stress conditions, and nitrogen response of ppc1b. Characterization of wheat ppc1b showed high levels of gene expression in young leaves, transcriptional responses under nitrogen and abiotic stress, and the presence of a Dof1 binding site, similar to ppcC4 in maize. Our results indicate the evolving status of Poaceae PEPCs and suggest the functional association of ppc1-derivatives with adaptation to environmental changes.

Comparative analysis reveals chromosome number reductions in the evolution of African bermudagrass (Cynodon transvaalensis Burtt-Davy)

African bermudagrass (Cynodon transvaalensis Burtt-Davy) (2n = 2x = 18) belongs to the genus Cynodon, tribe Cynodonteae, subfamily Chloridoideae in the grass family Poaceae. The species is frequently crossed with common bermudagrass (Cynodon dactylon Pers.) in developing high-quality hybrid turf cultivars. Molecular resources for C. transvaalensis are scarce; thus, its genomic evolution is unknown. Recently, a linkage map consisting of 1278 markers provided a powerful tool for African bermudagrass genomic research. The objective of this study was to investigate chromosome number reduction events that resulted in the nine haploid chromosomes in this species. Tag sequences of mapped single nucleotide polymorphism markers in C. transvaalensis were compared against genome sequences of Oropetium thomaeum (L.f.) Trin. (2n = 2x = 20), a genomic model in the Cynodonteae tribe. The comparative genomic analyses revealed broad collinearity between the genomes of these two species. The analyses further revealed that two major interchromosomal rearrangements of the paleochromosome ρ12 (ρ1-ρ12-ρ1 and ρ6-ρ12-ρ6) resulted in nine chromosomes in the genome of C. transvaalensis. The findings provide novel information regarding the formation of the initial diploid species in the Cynodon genus.

Evolution and origin of bread wheat

Bread wheat (Triticum aestivum, genome BBAADD) is a young hexaploid species formed only 8,500-9,000 years ago through hybridization between a domesticated free-threshing tetraploid progenitor, genome BBAA, and Aegilops tauschii, the diploid donor of the D subgenome. Very soon after its formation, it spread globally from its cradle in the fertile crescent into new habitats and climates, to become a staple food of humanity. This extraordinary global expansion was probably enabled by allopolyploidy that accelerated genetic novelty through the acquisition of new traits, new intergenomic interactions, and buffering of mutations, and by the attractiveness of bread wheat's large, tasty, and nutritious grain with high baking quality. New genome sequences suggest that the elusive donor of the B subgenome is a distinct (unknown or extinct) species rather than a mosaic genome. We discuss the origin of the diploid and tetraploid progenitors of bread wheat and the conflicting genetic and archaeological evidence on where it was formed and which species was its free-threshing tetraploid progenitor. Wheat experienced many environmental changes throughout its evolution, therefore, while it might adapt to current climatic changes, efforts are needed to better use and conserve the vast gene pool of wheat biodiversity on which our food security depends.

Evolutionary Dynamics of the Repeatome Explains Contrasting Differences in Genome Sizes and Hybrid and Polyploid Origins of Grass Loliinae Lineages

The repeatome is composed of diverse families of repetitive DNA that keep signatures on the historical events that shaped the evolution of their hosting species. The cold seasonal Loliinae subtribe includes worldwide distributed taxa, some of which are the most important forage and lawn species (fescues and ray-grasses). The Loliinae are prone to hybridization and polyploidization. It has been observed a striking two-fold difference in genome size between the broad-leaved (BL) and fine-leaved (FL) Loliinae diploids and a general trend of genome reduction of some high polyploids. We have used genome skimming data to uncover the composition, abundance, and potential phylogenetic signal of repetitive elements across 47 representatives of the main Loliinae lineages. Independent and comparative analyses of repetitive sequences and of 5S rDNA loci were performed for all taxa under study and for four evolutionary Loliinae groups [Loliinae, Broad-leaved (BL), Fine-leaved (FL), and Schedonorus lineages]. Our data showed that the proportion of the genome covered by the repeatome in the Loliinae species was relatively high (average ∼ 51.8%), ranging from high percentages in some diploids (68.7%) to low percentages in some high-polyploids (30.7%), and that changes in their genome sizes were likely caused by gains or losses in their repeat elements. Ty3-gypsy Retand and Ty1-copia Angela retrotransposons were the most frequent repeat families in the Loliinae although the relatively more conservative Angela repeats presented the highest correlation of repeat content with genome size variation and the highest phylogenetic signal of the whole repeatome. By contrast, Athila retrotransposons presented evidence of recent proliferations almost exclusively in the Lolium clade. The repeatome evolutionary networks showed an overall topological congruence with the nuclear 35S rDNA phylogeny and a geographic-based structure for some lineages. The evolution of the Loliinae repeatome suggests a plausible scenario of recurrent allopolyploidizations followed by diploidizations that generated the large genome sizes of BL diploids as well as large genomic rearrangements in highly hybridogenous lineages that caused massive repeatome and genome contractions in the Schedonorus and Aulaxyper polyploids. Our study has contributed to disentangling the impact of the repeatome dynamics on the genome diversification and evolution of the Loliinae grasses.

Chromosome-Scale Genome Assembly Provides Insights into Speciation of Allotetraploid and Massive Biomass Accumulation of Elephant Grass (Pennisetum purpureum Schum.)

Elephant grass (Pennisetum purpureum Schum) is an important forage, biofuels and industrial plant widely distributed in tropical and subtropical areas globally. It is characterized with robust growth and high biomass. We sequenced its allopolyploid genome and assembled 2.07 Gb into A' and B sub-genomes of 14 chromosomes with scaffold N50 of 8.47 Mb, yielding a total of 77,139 genes. The allotetraploid speciation occurred approximately 15 MYA after the divergence between Setaria italica and Pennisetum glaucum, according to a phylogenetic analysis of Pennisetum species. Double whole-genome duplication (WGD) and polyploidization events resulted in large scale gene expansion, especially in the key steps of growth and biomass accumulation. Integrated transcriptome profiling revealed the functional divergence between sub-genomes A' and B. A' sub-genome mainly contributed to plant growth, development and photosynthesis, whereas the B sub-genome was primarily responsible for effective transportation and resistance to stimulation. Some key gene families related to cellulose biosynthesis were expanded and highly expressed in stems, which could explain the high cellulose content in elephant grass. Our findings provide deep insights into genetic evolution of elephant grass and will aid future biological research and breeding, even for other grasses in the family Poaceae.

Conservation and Divergence of SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) Gene Family between Wheat and Rice

The SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) gene family affects plant architecture, panicle structure, and grain development, representing key genes for crop improvements. The objective of the present study is to utilize the well characterized SPLs’ functions in rice to facilitate the functional genomics of TaSPL genes. To achieve these goals, we combined several approaches, including genome-wide analysis of TaSPLs, comparative genomic analysis, expression profiling, and functional study of TaSPL3 in rice. We established the orthologous relationships of 56 TaSPL genes with the corresponding OsSPLs, laying a foundation for the comparison of known SPL functions between wheat and rice. Some TaSPLs exhibited different spatial–temporal expression patterns when compared to their rice orthologs, thus implicating functional divergence. TaSPL2/6/8/10 were identified to respond to different abiotic stresses through the combination of RNA-seq and qPCR expression analysis. Additionally, ectopic expression of TaSPL3 in rice promotes heading dates, affects leaf and stem development, and leads to smaller panicles and decreased yields per panicle. In conclusion, our work provides useful information toward cataloging of the functions of TaSPLs, emphasized the conservation and divergence between TaSPLs and OsSPLs, and identified the important SPL genes for wheat improvement.

Ancient rapid functional differentiation and fixation of the duplicated members in rice Dof genes after whole genome duplication

Whole genome duplication (WGD) in plants is typically followed by genomic downsizing, where large portions of the new genome are lost. Whether this downsizing is accompanied by increased or decreased evolutionary rates of the remaining genes is poorly known, not least because homeolog pairings are often obscured by chromosomal rearrangement. Here, we use the newly published genome from a sedge, namely Kobresia littledalei, and CRISPR/Cas-9 editing to investigate how the Rho WGD event 70 million years ago (MYA) affected transcription factor evolutionary rates, fates, and function in rice (Oryza sativa) and sorghum (Sorghum bicolor). We focus on the 30-member DNA-binding with one zinc finger (Dof) transcription factor family in both crops due to their agronomic importance. Using the known speciation dates of rice from Kobresia (97 MYA) and sorghum (50 MYA), we find that rates of amino acid substitution in the critical Dof domain region were over twofold higher during the 20-million-year period following the WGD than before or afterward. Through comparison of synteny blocks, we report that at least 11% of Dof genes were purged from 70 to 50 MYA, while only 6% have been lost in the most recent 50-million-year interval. CRISPR/Cas9 editing revealed widespread fitness-related defects in flowering and lack of redundancy of paired members, as well as significant differences in expression between gene pairs. Together these findings demonstrate the strength of Dof genes as a model for deep evolutionary study and offer one of the most detailed portraits yet of the Rho WGD impact on a gene lineage.

The Streptochaeta Genome and the Evolution of the Grasses

In this work, we sequenced and annotated the genome of Streptochaeta angustifolia, one of two genera in the grass subfamily Anomochlooideae, a lineage sister to all other grasses. The final assembly size is over 99% of the estimated genome size. We find good collinearity with the rice genome and have captured most of the gene space. Streptochaeta is similar to other grasses in the structure of its fruit (a caryopsis or grain) but has peculiar flowers and inflorescences that are distinct from those in the outgroups and in other grasses. To provide tools for investigations of floral structure, we analyzed two large families of transcription factors, AP2-like and R2R3 MYBs, that are known to control floral and spikelet development in rice and maize among other grasses. Many of these are also regulated by small RNAs. Structure of the gene trees showed that the well documented whole genome duplication at the origin of the grasses (ρ) occurred before the divergence of the Anomochlooideae lineage from the lineage leading to the rest of the grasses (the spikelet clade) and thus that the common ancestor of all grasses probably had two copies of the developmental genes. However, Streptochaeta (and by inference other members of Anomochlooideae) has lost one copy of many genes. The peculiar floral morphology of Streptochaeta may thus have derived from an ancestral plant that was morphologically similar to the spikelet-bearing grasses. We further identify 114 loci producing microRNAs and 89 loci generating phased, secondary siRNAs, classes of small RNAs known to be influential in transcriptional and post-transcriptional regulation of several plant functions.

A high-quality chromosome-scale assembly of the centipedegrass [Eremochloa ophiuroides (Munro) Hack.] genome provides insights into chromosomal structural evolution and prostrate growth habit

Centipedegrass [Eremochloa ophiuroides (Munro) Hack.], a member of the Panicoideae subfamily, is one of the most important warm-season turfgrasses originating from China. This grass has an extremely developed prostrate growth habit and has been widely used in transitional and warm climatic regions. To better understand the genetic basis of important biological characteristics, such as prostrate growth and seed yield, in warm-season turfgrasses, we present a high-quality reference genome for centipedegrass and use PacBio, BioNano, and Hi-C technologies to anchor the 867.43 Mb genome assembly into nine pseudochromosomes, with a scaffold N50 of 86.05 Mb and 36,572 annotated genes. Centipedegrass was most closely related to sorghum and diverged from their common ancestor ~16.8 Mya. We detected a novel chromosome reshuffling event in centipedegrass, namely, the nest chromosome fusion event in which fusion of chromosomes 8 and 10 of sorghum into chromosome 3 of centipedegrass likely occurred after the divergence of centipedegrass from sorghum. The typical prostrate growth trait in centipedegrass may be linked to the expansion of candidate PROSTRATE GROWTH 1 (PROG1) genes on chromosome 2. Two orthologous genes of OsPROG1, EoPROG1, and EoPROG2, were confirmed to increase the stem number and decrease the stem angle in Arabidopsis. Collectively, our assembled reference genome of centipedegrass offers new knowledge and resources to dissect the genome evolution of Panicoideae and accelerate genome-assisted breeding and improvement of plant architecture in turf plants.

Whole-genome assembly and annotation of northern wild rice, Zizania palustris L., supports a whole-genome duplication in the Zizania genus

Zizania palustris L. (northern wild rice, NWR) is an aquatic grass native to North America that is notable for its nutritious grain. This is an important species with ecological, cultural and agricultural significance, specifically in the Great Lakes region of the USA. Using flow cytometry, we first estimated the NWR genome size to be 1.8 Gb. Using long- and short-range sequencing, Hi-C scaffolding and RNA-seq data from eight tissues, we generated an annotated whole-genome de novo assembly of NWR. The assembly was 1.29 Gb in length, highly repetitive (approx. 76.0%) and contained 46 421 putative protein-coding genes. The expansion of retrotransposons within the genome and a whole-genome duplication (WGD) after the Zizania-Oryza speciation event have both led to an increase in the genome size of NWR in comparison with Oryza sativa L. and Zizania latifolia. Both events depict a genome rapidly undergoing change over a short evolutionary time. Comparative analyses revealed the conservation of large syntenic blocks between NWR and O. sativa, which were used to identify putative seed-shattering genes. Estimates of divergence times revealed that the Zizania genus diverged from Oryza approximately 26-30 million years ago (26-30 MYA), whereas NWR and Z. latifolia diverged from one another approximately 6-8 MYA. Comparative genomics confirmed evidence of a WGD in the Zizania genus and provided support that the event occurred prior to the NWR-Z. latifolia speciation event. This genome assembly and annotation provides a valuable resource for comparative genomics in the Oryzeae tribe and provides an important resource for future conservation and breeding efforts of NWR.

Genome-wide identification of alcohol dehydrogenase (ADH) gene family under waterlogging stress in wheat (Triticum aestivum)

Background

Alcohol dehydrogenase (ADH) plays an important role in plant survival under anaerobic conditions. Although some research about ADH in many plants have been carried out, the bioinformatics analysis of the ADH gene family from Triticum aestivum and their response to abiotic stress is unclear.

Methods

A total of 22 ADH genes were identified from the wheat genome, and these genes could be divided into two subfamilies (subfamily I and subfamily II). All TaADH genes belonged to the Medium-chain ADH subfamily. Sequence alignment analysis showed that all TaADH proteins contained a conservative GroES-like domain and Zinc-binding domain. A total of 64 duplicated gene pairs were found, and the Ka/Ks value of these gene pairs was less than 1, which indicated that these genes were relatively conservative and did not change greatly in the process of duplication.

Results

The organizational analysis showed that nine TaADH genes were highly expressed in all organs, and the rest of TaADH genes had tissue specificity. Cis-acting element analysis showed that almost all of the TaADH genes contained an anaerobic response element. The expression levels of ADH gene in waterlogging tolerant and waterlogging sensitive wheat seeds were analyzed by quantitative real-time PCR (qRT-PCR). This showed that some key ADH genes were significantly responsive to waterlogging stress at the seed germination stage, and the response of waterlogging tolerant and waterlogging sensitive wheat seeds to waterlogging stress was regulated by different ADH genes. The results may be helpful to further study the function of TaADH genes and to determine the candidate gene for wheat stress resistance breeding.

Introduction of glucan synthase into the cytosol in wheat endosperm causes massive maltose accumulation and represses starch synthesis

We expressed a bacterial glucan synthase (Agrobacterium GlgA) in the cytosol of developing endosperm cells in wheat grains, to discover whether it could generate a glucan from cytosolic ADP-glucose. Transgenic lines had high glucan synthase activity during grain filling, but did not accumulate glucan. Instead, grains accumulated very high concentrations of maltose. They had large volumes during development due to high water content, and very shrivelled grains at maturity. Starch synthesis was severely reduced. We propose that cytosolic glucan synthesized by the glucan synthase was immediately hydrolysed to maltose by cytosolic β-amylase(s). Maltose accumulation resulted in a high osmotic potential in developing grain, drawing in excess water that stretched the seed coat and pericarp. Loss of water during grain maturation then led to shrinkage when the grains matured. Maltose accumulation is likely to account for the reduced starch synthesis in transgenic grains, through signalling and toxic effects. Using bioinformatics, we identify an isoform of β-amylase likely to be responsible for maltose accumulation. Removal of this isoform through identification of TILLING mutants or genome editing, combined with co-expression of heterologous glucan synthase and a glucan branching enzyme, may in future enable elevated yields of carbohydrate through simultaneous accumulation of starch and cytosolic glucan.

The Pharus latifolius genome bridges the gap of early grass evolution

The grass family (Poaceae) includes all commercial cereal crops and is a major contributor to biomass in various terrestrial ecosystems. The ancestry of all grass genomes includes a shared whole-genome duplication (WGD), named rho (ρ) WGD, but the evolutionary significance of ρ-WGD remains elusive. We sequenced the genome of Pharus latifolius, a grass species (producing a true spikelet) in the subfamily Pharoideae, a sister lineage to the core Poaceae including the (Panicoideae, Arundinoideae, Chloridoideae, Micrairoideae, Aristidoideae, and Danthonioideae (PACMAD) and Bambusoideae, Oryzoideae, and Pooideae (BOP) clades. Our results indicate that the P. latifolius genome has evolved slowly relative to cereal grass genomes, as reflected by moderate rates of molecular evolution, limited chromosome rearrangements and a low rate of gene loss for duplicated genes. We show that the ρ-WGD event occurred approximately 98.2 million years ago (Ma) in a common ancestor of the Pharoideae and the PACMAD and BOP grasses. This was followed by contrasting patterns of diploidization in the Pharus and core Poaceae lineages. The presence of two FRIZZY PANICLE-like genes in P. latifolius, and duplicated MADS-box genes, support the hypothesis that the ρ-WGD may have played a role in the origin and functional diversification of the spikelet, an adaptation in grasses related directly to cereal yields. The P. latifolius genome sheds light on the origin and early evolution of grasses underpinning the biology and breeding of cereals.

Transcriptome Analysis of Tetraploid and Octoploid Common Reed (Phragmites australis)

Polyploidization in plants is thought to have occurred as coping mechanism with environmental stresses. Polyploidization-driven adaptation is often achieved through interplay of gene networks involved in differentially expressed genes, which triggers the plant to evolve special phenotypic traits for survival. Phragmites australis is a cosmopolitan species with highly variable phenotypic traits and high adaptation capacity to various habitats. The species' ploidy level varies from 3x to 12x, thus it is an ideal organism to investigate the molecular evolution of polyploidy and gene regulation mediated by different numbers of chromosome copies. In this study, we used high-throughput RNAseq data as a tool, to analyze the gene expression profiles in tetraploid and octoploid P. australis. The estimated divergence time between tetraploid and octoploid P. australis was dated to the border between Pliocene and Pleistocene. This study identified 439 up- and 956 down-regulated transcripts in tetraploids compared to octoploids. Gene ontology and pathway analysis revealed that tetraploids tended to express genes responsible for reproduction and seed germination to complete the reproduction cycle early, and expressed genes related to defense against UV-B light and fungi, whereas octoploids expressed mainly genes related to thermotolerance. Most differentially expressed genes were enriched in chaperones, folding catalysts and protein processing in endoplasmic reticulum pathways. Multiple biased isoform usage of the same gene was detected in differentially expressed genes, and the ones upregulated in octoploids were related to reduced DNA methylation. Our study provides new insights into the role of polyploidization on environmental responses and potential stress tolerance in grass species.

The Landscape of the Genomic Distribution and the Expression of the F-Box Genes Unveil Genome Plasticity in Hexaploid Wheat during Grain Development and in Response to Heat and Drought Stress

FBX proteins are subunits of the SCF complex (Skp1-cullin-FBX) belonging to the E3 ligase family, which is involved in the ubiquitin-proteasome 26S (UPS) pathway responsible for the post-translational protein turnover. By targeting, in a selective manner, key regulatory proteins for ubiquitination and 26S proteasome degradation, FBX proteins play a major role in plant responses to diverse developmental and stress conditions. Although studies on the genomic organization of the FBX gene family in various species have been reported, knowledge related to bread wheat (Triticum aestivum) is scarce and needs to be broadened. Using the latest assembly of the wheat genome, we identified 3670 TaFBX genes distributed non-homogeneously within the three subgenomes (A, B and D) and between the 21 chromosomes, establishing it as one of the richest gene families among plant species. Based on the presence of the five different chromosomal regions previously identified, the present study focused on the genomic distribution of the TaFBX family and the identification of differentially expressed genes during the embryogenesis stages and in response to heat and drought stress. Most of the time, when comparing the expected number of genes (taking into account the formal gene distribution on the entire wheat genome), the TaFBX family harbors a different pattern at the various stratum of observation (subgenome, chromosome, chromosomal regions). We report here that the local gene expansion of the TaFBX family must be the consequence of multiple and complex events, including tandem and small-scale duplications. Regarding the differentially expressed TaFBX genes, while the majority of the genes are localized in the distal chromosomal regions (R1 and R3), differentially expressed genes are more present in the interstitial regions (R2a and R2b) than expected, which could be an indication of the preservation of major genes in those specific chromosomal regions.

Genetics and genomics of root system variation in adaptation to drought stress in cereal crops

Cereals are important crops worldwide that help meet food demands and nutritional needs. In recent years, cereal production has been challenged globally by frequent droughts and hot spells. A plant's root is the most relevant organ for the plant adaptation to stress conditions, playing pivotal roles in anchorage and the acquisition of soil-based resources. Thus, dissecting root system variations and trait selection for enhancing yield and sustainability under drought stress conditions should aid in future global food security. This review highlights the variations in root system attributes and their interplay with shoot architecture features to face water scarcity and maintain thus yield of major cereal crops. Further, we compile the root-related drought responsive quantitative trait loci/genes in cereal crops including their interspecies relationships using microsynteny to facilitate comparative genomic analyses. We then discuss the potential of an integrated strategy combining genomics and phenomics at genetic and epigenetic levels to explore natural genetic diversity as a basis for knowledge-based genome editing. Finally, we present an outline to establish innovative breeding leads for the rapid and optimized selection of root traits necessary to develop resilient crop varieties.

A high-quality genome assembly highlights rye genomic characteristics and agronomically important genes

Rye is a valuable food and forage crop, an important genetic resource for wheat and triticale improvement and an indispensable material for efficient comparative genomic studies in grasses. Here, we sequenced the genome of Weining rye, an elite Chinese rye variety. The assembled contigs (7.74 Gb) accounted for 98.47% of the estimated genome size (7.86 Gb), with 93.67% of the contigs (7.25 Gb) assigned to seven chromosomes. Repetitive elements constituted 90.31% of the assembled genome. Compared to previously sequenced Triticeae genomes, Daniela, Sumaya and Sumana retrotransposons showed strong expansion in rye. Further analyses of the Weining assembly shed new light on genome-wide gene duplications and their impact on starch biosynthesis genes, physical organization of complex prolamin loci, gene expression features underlying early heading trait and putative domestication-associated chromosomal regions and loci in rye. This genome sequence promises to accelerate genomic and breeding studies in rye and related cereal crops.

CRISPR-based assessment of genomic structure in the conserved SQUAMOSA promoter-binding-like gene clusters in rice

Although SQUAMOSA promoter-binding-like (SPL) transcription factors are important regulators of development in rice (Oryza sativa), prior assessments of the SPL family have been limited to single genes. A functional comparison across the full gene family in standardized genetic backgrounds has not been reported previously. Here, we demonstrate that the SPL gene family in rice is enriched due to the most recent whole genome duplication (WGD). Notably, 10 of 19 rice SPL genes (52%) cluster in four units that have persisted for at least 50 million years. We show that SPL gene grouping and retention following WGD is widespread in angiosperms, suggesting the conservatism and importance of this gene arrangement. We used Cas9 editing to generate transformation lines for all 19 SPL genes in a common set of backgrounds, and found that knockouts of 14 SPL genes exhibited defects in plant height, 10 exhibited defects in panicle size, and nine had altered grain lengths. We observed subfunctionalization of genes in the paleoduplicated pairs, but little evidence of neofunctionalization. Expression of OsSPL3 was negatively correlated with that of its closest neighbor in its synteny group, OsSPL4, and its sister paired gene, OsSPL12, in the opposing group. Nucleotide diversity was lower in eight of the nine singleton genes in domesticated rice, relative to wild rice, whereas the reverse was true for the paired genes. Together, these results provide functional information on eight previously unexamined OsSPL family members and suggest that paleoduplicate pair redundancy benefits plant survival and innovation.

Evolution of Homeologous Gene Expression in Polyploid Wheat

Polyploidization has played a prominent role in the evolutionary history of plants. Two recent and sequential allopolyploidization events have resulted in the formation of wheat species with different ploidies, and which provide a model to study the effects of polyploidization on the evolution of gene expression. In this study, we identified differentially expressed genes (DEGs) between four BBAA tetraploid wheats of three different ploidy backgrounds. DEGs were found to be unevenly distributed among functional categories and duplication modes. We observed more DEGs in the extracted tetraploid wheat (ETW) than in natural tetraploid wheats (TD and TTR13) as compared to a synthetic tetraploid (AT2). Furthermore, DEGs showed higher Ka/Ks ratios than those that did not show expression changes (non-DEGs) between genotypes, indicating DEGs and non-DEGs experienced different selection pressures. For A-B homeolog pairs with DEGs, most of them had only one differentially expressed copy, however, when both copies of a homeolog pair were DEGs, the A and B copies were more likely to be regulated to the same direction. Our results suggest that both cis- and inter-subgenome trans-regulatory changes are important drivers in the evolution of homeologous gene expression in polyploid wheat, with ploidy playing a significant role in the process.

Molecular characterization of a RING E3 ligase SbHCI1 in sorghum under heat and abscisic acid stress

Molecular function ofRING E3 ligase SbHCI1is involved in ABA-mediated basal heat stress tolerancein sorghum. Global warming generally reduces plant survival, owing to the negative effects of high temperatures on plant development. However, little is known about the role of Really Interesting New Gene (RING) E3 ligase in the heat stress responses of plants. As such, the aim of the present study was to characterize the molecular functions of the Sorghum bicolor ortholog of the Oryza sativa gene for Heat- and Cold-Induced RING finger protein 1 (SbHCI1). Subcellular localization revealed that SbHCI1 was mainly associated with the cytosol and that it moved to the Golgi apparatus under heat stress conditions. The fluorescent signals of SbHCI1 substrate proteins were observed to migrate to the cytoplasm under heat stress conditions. Bimolecular fluorescence complementation (BiFC) and yeast two-hybrid (Y2H) assays revealed that SbHCI1 physically interacted with OsHCI1 ortholog partner proteins in the cytoplasm. Moreover, an in vitro ubiquitination assay revealed that SbHCI1 polyubiquitinated each of the three interacting proteins. The ectopic overexpression of SbHCI1 in Arabidopsis revealed that the protein was capable of inducing abscisic acid (ABA)-hypersensitivity and basal heat stress tolerance. Therefore, SbHCI1 possesses E3 ligase activity and may function as a positive regulator of heat stress responses through the modulation of interacting proteins.

The MYB transcription factor TaPHR3-A1 is involved in phosphate signaling and governs yield-related traits in bread wheat

Improved inorganic phosphate (Pi) use efficiency in crops will be important for sustainable agriculture. Exploring molecular mechanisms that regulate Pi uptake could provide useful information for breeding wheat with improved Pi use efficiency. Here, a TaPHR3-A1 (Gene ID: TraesCS7A02G415800) ortholog of rice OsPHR3 that functions in transcriptional regulation of Pi signaling was cloned from wheat chromosome 7A. Ectopic expression of TaPHR3-A1 in Arabidopsis and rice produced enhanced vegetative growth and more seeds. Overexpression in transgenic rice led to increased biomass, grain number, and primary panicle branching by 61.23, 42.12, and 36.34% compared with the wild type. Transgenic wheat lines with down-regulation of TaPHR3-A1 exhibited retarded growth and root hair development at the seedling stage, and showed yield-related effects at the adult stage when grown in both low- and sufficient Pi conditions, indicating that TaPHR3-A1 positively regulated tolerance to low Pi. Introgression lines further confirmed the effect of TaPHR3-A1 in improving grain number. The Chinese wheat mini core collection and a recombinant inbred line analysis demonstrated that the favorable allele TaPHR3-A1-A associated with higher grain number was positively selected in breeding. A TaPHR3-A1-derived cleaved amplified polymorphic sequence marker effectively identified haplotype TaPHR3-A1-A. Our results suggested that TaPHR3-A1 was a functional regulatory factor for Pi uptake and provided useful information for marker-assisted selection for high yield in wheat.

Segmental duplications are hot spots of copy number variants affecting barley gene content

Copy number variants (CNVs) are pervasive in several animal and plant genomes and contribute to shaping genetic diversity. In barley, there is evidence that changes in gene copy number underlie important agronomic traits. The recently released reference sequence of barley represents a valuable genomic resource for unveiling the incidence of CNVs that affect gene content and for identifying sequence features associated with CNV formation. Using exome sequencing and read count data, we detected 16 605 deletions and duplications that affect barley gene content by surveying a diverse panel of 172 cultivars, 171 landraces, 22 wild relatives and other 32 uncategorized domesticated accessions. The quest for segmental duplications (SDs) in the reference sequence revealed many low-copy repeats, most of which overlap predicted coding sequences. Statistical analyses revealed that the incidence of CNVs increases significantly in SD-rich regions, indicating that these sequence elements act as hot spots for the formation of CNVs. The present study delivers a comprehensive genome-wide study of CNVs affecting barley gene content and implicates SDs in the molecular mechanisms that lead to the formation of this class of CNVs.

Analysis of local genome rearrangement improves resolution of ancestral genomic maps in plants

BACKGROUND

Computationally inferred ancestral genomes play an important role in many areas of genome research. We present an improved workflow for the reconstruction from highly diverged genomes such as those of plants.

RESULTS

Our work relies on an established workflow in the reconstruction of ancestral plants, but improves several steps of this process. Instead of using gene annotations for inferring the genome content of the ancestral sequence, we identify genomic markers through a process called genome segmentation. This enables us to reconstruct the ancestral genome from hundreds of thousands of markers rather than the tens of thousands of annotated genes. We also introduce the concept of local genome rearrangement, through which we refine syntenic blocks before they are used in the reconstruction of contiguous ancestral regions. With the enhanced workflow at hand, we reconstruct the ancestral genome of eudicots, a major sub-clade of flowering plants, using whole genome sequences of five modern plants.

CONCLUSIONS

Our reconstructed genome is highly detailed, yet its layout agrees well with that reported in Badouin et al. (2017). Using local genome rearrangement, not only the marker-based, but also the gene-based reconstruction of the eudicot ancestor exhibited increased genome content, evidencing the power of this novel concept.

Flow cytometric characterisation of the complex polyploid genome of Saccharum officinarum and modern sugarcane cultivars

Sugarcane (Saccharum spp.) is a globally important crop for sugar and bioenergy production. Its highly polyploid, complex genome has hindered progress in understanding its molecular structure. Flow cytometric sorting and analysis has been used in other important crops with large genomes to dissect the genome into component chromosomes. Here we present for the first time a method to prepare suspensions of intact sugarcane chromosomes for flow cytometric analysis and sorting. Flow karyotypes were generated for two S. officinarum and three hybrid cultivars. Five main peaks were identified and each genotype had a distinct flow karyotype profile. The flow karyotypes of S. officinarum were sharper and with more discrete peaks than the hybrids, this difference is probably due to the double genome structure of the hybrids. Simple Sequence Repeat (SSR) markers were used to determine that at least one allelic copy of each of the 10 basic chromosomes could be found in each peak for every genotype, except R570, suggesting that the peaks may represent ancestral Saccharum sub genomes. The ability to flow sort Saccharum chromosomes will allow us to isolate and analyse chromosomes of interest and further examine the structure and evolution of the sugarcane genome.

Genome-wide identification and expression profiling of glutathione transferase gene family under multiple stresses and hormone treatments in wheat (Triticum aestivum L.)

BACKGROUND

Glutathione transferases (GSTs), the ancient, ubiquitous and multi-functional proteins, play significant roles in development, metabolism as well as abiotic and biotic stress responses in plants. Wheat is one of the most important crops, but the functions of GST genes in wheat were less studied.

RESULTS

A total of 330 TaGST genes were identified from the wheat genome and named according to the nomenclature of rice and Arabidopsis GST genes. They were classified into eight classes based on the phylogenetic relationship among wheat, rice, and Arabidopsis, and their gene structure and conserved motif were similar in the same phylogenetic class. The 43 and 171 gene pairs were identified as tandem and segmental duplication genes respectively, and the Ka/Ks ratios of tandem and segmental duplication TaGST genes were less than 1 except segmental duplication gene pair TaGSTU24/TaGSTU154. The 59 TaGST genes were identified to have syntenic relationships with 28 OsGST genes. The expression profiling involved in 15 tissues and biotic and abiotic stresses suggested the different expression and response patterns of the TaGST genes. Furthermore, the qRT-PCR data showed that GST could response to abiotic stresses and hormones extensively in wheat.

CONCLUSIONS

In this study, a large GST family with 330 members was identified from the wheat genome. Duplication events containing tandem and segmental duplication contributed to the expansion of TaGST family, and duplication genes might undergo extensive purifying selection. The expression profiling and cis-elements in promoter region of 330 TaGST genes implied their roles in growth and development as well as adaption to stressful environments. The qRT-PCR data of 14 TaGST genes revealed that they could respond to different abiotic stresses and hormones, especially salt stress and abscisic acid. In conclusion, this study contributed to the further functional analysis of GST genes family in wheat.

Molecular Functions of Rice Cytosol-Localized RING Finger Protein 1 in Response to Salt and Drought and Comparative Analysis of Its Grass Orthologs

In higher plants, the post-translational modification of target proteins via the attachment of molecules such as ubiquitin (Ub) mediates a variety of cellular functions via the Ub/26S proteasome system. Here, a really interesting new gene (RING)-H2 type E3 ligase, which regulates target proteins via the Ub/26S proteasome system, was isolated from a rice plant, and its other grass orthologs were examined to determine the evolution of its molecular function during speciation. The gene encoding Oryza sativa cytoplasmic-localized RING finger protein 1 (OsCLR1) was highly expressed under salt and drought stresses. By contrast, the three grass orthologs, SbCLR1 from Sorghum bicolor, ZmCLR1 from Zea mays and TaCLR1 from Triticum aestivum, showed different responses to these stresses. Despite these differences, all four orthologs exhibited E3 ligase activity with cytosol-targeted localization, demonstrating conserved molecular functions. Although OsCLR1-overexpressing plants showed higher survival rates under both salt and drought stresses than that of the wild type (WT) plants, this pattern was not observed in the other orthologs. In addition, OsCLR1-overexpressing plants exhibited lower germination rates in ABA than that of WT plants, whereas the three ortholog CLR1-overexpressing plants showed rates similar to the WT plants. These results indicate the positive regulation of OsCLR1 in response to salt and drought in an ABA-dependent manner. Despite the molecular functions of the three CLR1 orthologs remaining largely unknown, our results provide an insight into the evolutionary fate of CLR1 grass orthologs during speciation after the divergence from a common ancestor.

Sequence properties of certain GC rich avian genes, their origins and absence from genome assemblies: case studies

Background

More and more eukaryotic genomes are sequenced and assembled, most of them presented as a complete model in which missing chromosomal regions are filled by Ns and where a few chromosomes may be lacking. Avian genomes often contain sequences with high GC content, which has been hypothesized to be at the origin of many missing sequences in these genomes. We investigated features of these missing sequences to discover why some may not have been integrated into genomic libraries and/or sequenced.

Results

The sequences of five red jungle fowl cDNA models with high GC content were used as queries to search publicly available datasets of Illumina and Pacbio sequencing reads. These were used to reconstruct the leptin, TNFα, MRPL52, PCP2 and PET100 genes, all of which are absent from the red jungle fowl genome model. These gene sequences displayed elevated GC contents, had intron sizes that were sometimes larger than non-avian orthologues, and had non-coding regions that contained numerous tandem and inverted repeat sequences with motifs able to assemble into stable G-quadruplexes and intrastrand dyadic structures. Our results suggest that Illumina technology was unable to sequence the non-coding regions of these genes. On the other hand, PacBio technology was able to sequence these regions, but with dramatically lower efficiency than would typically be expected.

Conclusions

High GC content was not the principal reason why numerous GC-rich regions of avian genomes are missing from genome assembly models. Instead, it is the presence of tandem repeats containing motifs capable of assembling into very stable secondary structures that is likely responsible.

Comparative genome analysis of the SPL gene family reveals novel evolutionary features in maize

SPLs are plant-specific transcription factors that play important regulatory roles in plant growth and development. Systematic analysis of the SPL family has been performed in numerous plants, such as Arabidopsis, rice, and Populus. However, no comparative analysis has been performed across different species to examine evolutionary features. In this study, we present a comparative analysis of SPLs in different species. The results showed that 84 SPLs of the four species can be divided into six groups according to phylogeny. We found that most of the SPL-containing regions in maize showed extensive conservation with duplicated regions of rice and sorghum. A gene duplication analysis in maize indicated that ZmSPLs showed a significant excess of segmental duplication. The Ka/Ks analysis indicated that 9 out of 18 duplicated pairs in maize experienced positive selection, while SPL gene pairs of rice and sorghum mainly evolved under purifying selection, suggesting novel evolutionary features for ZmSPLs. The 31 ZmSPLs were further analyzed by describing their gene structure, phylogenetic relationships, chromosomal location, and expression, Among the ZmSPLs, 13 were predicated to be targeted by miR156s and involved in drought stress response. These results provide the foundation for future functional analyses of ZmSPLs.

New insights into the origin and evolution of α-amylase genes in green plants

Gene duplication is a source of genetic materials and evolutionary changes, and has been associated with gene family expansion. Functional divergence of duplicated genes is strongly directed by natural selections such as organism diversification and novel feature acquisition. We show that, plant α-amylase gene family (AMY) is comprised of six subfamilies (AMY1-AMY6) that fell into two ancient phylogenetic lineages (AMY3 and AMY4). Both AMY1 and AMY2 are grass-specific and share a single-copy ancestor, which is derived from grass AMY3 genes that have undergone massive tandem and whole-genome duplications during evolution. Ancestral features of AMY4 and AMY5/AMY6 genes have been retained among four green algal sequences (Chrein_08.g362450, Vocart_0021s0194, Dusali_0430s00012 and Monegl_16464), suggesting a gene duplication event following Chlorophyceae diversification. The observed horizontal gene transfers between plant and bacterial AMYs, and chromosomal locations of AMY3 and AMY4 genes in the most ancestral green body (C. reinhardtii), provide evidences for the monophyletic origin of plant AMYs. Despite subfamily-specific sequence divergence driven by natural selections, the active site and SBS1 are well-conserved across different AMY isoforms. The differentiated electrostatic potentials and hydrogen bands-forming residue polymorphisms, further imply variable digestive abilities for a broad substrates in particular tissues or subcellular localizations.

A study of CCD8 genes/proteins in seven monocots and eight dicots

In plants, the enzyme CCD8 (carotenoid cleavage dioxygenase 8) is involved in the synthesis of an important hormone, strigolactone, and therefore, plays an important role in controlling growth and development. Using cDNA and protein sequence derived from the gene ZmCCD8 from maize, we identified putative orthologs of the gene encoding CCD8 in six other monocots and eight dicots; the sequence similarity ranged from 52–75.9% at the gene level and 60.9–93.7% at the protein level. The average length of the gene was ~3.3 kb (range: 2.08 to 3.98 kb), although the number of introns within the genes differed (4 or 5 in dicots and 3 or 4 in monocots, except in T. urartu with 6 introns). Several cis-acting regulatory elements were identified in the promoters of CCD8 genes, which are known to respond to biotic and abiotic stresses. The N-terminal end (up to ~70 amino acids) of CCD8 proteins was highly variable due to insertions, deletions and mismatches. The variation in genes and proteins were particularly conspicuous in T. urartu and Ae. tauschii among the monocots and A. thaliana and P. persica among the dicots. In CCD8 proteins, 12 motifs were also identified, of which 6 were novel; 4 of these novel motifs occurred in all the 15 species. The 3D structures of proteins had the characteristic features of the related enzyme apocarotenoid oxygenase (ACO) of Synechocystis (a representative of cyanobacteria). The results of qRT-PCR in wheat revealed that under phosphorous (P)-starved condition (relative to expression under optimum P used as control), the expression of TaCCD8 genes increased ~37 fold in root tissue of the cultivar C306 and ~33 fold in shoot tissue of the cultivar HUW468 (the two cultivars differed in their P-use efficiency). This suggested that expression of TaCCD8 genes is genotype-dependent and tissue-specific and is regulated under different levels of P supply.