Sequencing of chloroplast genomes from wheat, barley, rye and their relatives provides a detailed insight into the evolution of the Triticeae tribe

Selection of dysfunctional alleles of bHLH1 and MYB1 has produced white grain in the tribe Triticeae

Grain color is a key agronomic trait that greatly determines food quality. The molecular and evolutionary mechanisms that underlie grain-color regulation are also important questions in evolutionary biology and crop breeding. Here, we confirm that both bHLH and MYB genes have played a critical role in the evolution of grain color in Triticeae. Blue grain is the ancestral trait in Triticeae, whereas white grain caused by bHLH or MYB dysfunctions is the derived trait. HvbHLH1 and HvMYB1 have been the targets of selection in barley, and dysfunctions caused by deletion(s), insertion(s), and/or point mutation(s) in the vast majority of Triticeae species are accompanied by a change from blue grain to white grain. Wheat with white grains exhibits high seed vigor under stress. Artificial co-expression of ThbHLH1 and ThMYB1 in the wheat endosperm or aleurone layer can generate purple grains with health benefits and blue grains for use in a new hybrid breeding technology, respectively. Our study thus reveals that white grain may be a favorable derived trait retained through natural or artificial selection in Triticeae and that the ancient blue-grain trait could be regained and reused in molecular breeding of modern wheat.

A barley MLA immune receptor is activated by a fungal nonribosomal peptide effector for disease susceptibility

The barley Mla locus contains functionally diversified genes that encode intracellular nucleotide-binding leucine-rich repeat receptors (NLRs) and confer strain-specific immunity to biotrophic and hemibiotrophic fungal pathogens. In this study, we isolated a barley gene Scs6, which is an allelic variant of Mla genes but confers susceptibility to the isolate ND90Pr (BsND90Pr) of the necrotrophic fungus Bipolaris sorokiniana. We generated Scs6 transgenic barley lines and showed that Scs6 is sufficient to confer susceptibility to BsND90Pr in barley genotypes naturally lacking the receptor. The Scs6-encoded NLR (SCS6) is activated by a nonribosomal peptide (NRP) effector produced by BsND90Pr to induce cell death in barley and Nicotiana benthamiana. Domain swaps between MLAs and SCS6 reveal that the SCS6 leucine-rich repeat domain is a specificity determinant for receptor activation by the NRP effector. Scs6 is maintained in both wild and domesticated barley populations. Our phylogenetic analysis suggests that Scs6 is a Hordeum-specific innovation. We infer that SCS6 is a bona fide immune receptor that is likely directly activated by the nonribosomal peptide effector of BsND90Pr for disease susceptibility in barley. Our study provides a stepping stone for the future development of synthetic NLR receptors in crops that are less vulnerable to modification by necrotrophic pathogens.

The First Complete Chloroplast Genome Sequence of Secale strictum subsp. africanum Stapf (Poaceae), the Putative Ancestor of the Genus Secale

Secale strictum ssp. africanum (synonym Secale africanum), a putative ancestor of the genus Secale, has been classified within Secale strictum, although recent phylogenetic studies suggest that it represents a distinct species. This study reports the first complete chloroplast genome of S. africanum, highlighting its structure, genetic composition, and phylogenetic relationships within Secale and related Triticiceae species. Phylogeny reconstruction based on the maximum-likelihood method reveals notable genetic similarity between S. strictum and S. africanum, supporting their genetic and phylogenetic distinction. Here, we assembled the complete, annotated chloroplast genome sequence of Secale strictum ssp. africanum. The genome is 137,068 base pair (bp) long. It is the first complete chloroplast genome that can be used as a reference genome for further analysis. The genome can be accessed on GenBank with the accession number OQ700974. This work sheds light on the evolutionary history of Secale and contributes to our understanding of chloroplast genomics in cereal ancestors, with potential applications in improving cereal crop resilience, advancing breeding strategies, and informing conservation efforts for genetic diversity.

Comparative analysis of Aegilops speltoides and wheat repetitive elements and development of S genome-specific FISH painting

Aegilops speltoides (2n = 2x = 14, genome SS) is a wild relative of wheat and a donor of useful traits for wheat improvement. Several whole-genome studies compared genic regions of Aegilops from the Sitopsis section and wheat and found that Ae. speltoides is most closely related to the wheat B subgenome but is not its direct progenitor. The results showed that a B subgenome ancestor diverged from Ae. speltoides more than 4 MYA and either has not yet been discovered, or is extinct. To further explore the evolutionary relationship between wheat and Ae. speltoides and develop Ae. speltoides chromosome paints, we performed comparative analysis of repetitive fractions of the S genome and three subgenomes of hexaploid wheat. The low-coverage sequence data were analyzed with RepeatExplorer pipeline to annotate repeats and estimate their content. The LTR-retrotransposons comprised about 80% of repeats in Ae. speltoides and wheat datasets and about two-third of them were LTR/Ty3-Gypsy. Ae. speltoides had 1.5 times more LTR/Ty-Copia repeats and 1.5 times less DNA transposons than wheat subgenomes. Several S genome-specific dispersed repeats were found and annotated. Their sequences were used to develop S genome-specific paints for detecting Ae. speltoides chromatin in the wheat background using fluorescent in situ hybridization.

Tandemly duplicated Rubisco activase genes of cereals show differential evolution and response to heat stress

Heat stress affects various components of photosynthetic machinery of which Rubisco activation inhibition due to heat sensitive Rubisco activase (RCA) is the most prominent. Detailed comparison of RCA coding genes identified a tandem duplication event in the grass family lineage where the duplicated genes showed very different evolutionary pattern. One of the two genes showed high level of sequence conservation whereas the second copy, although present only 1.5 kb away, was highly variable among various plant species because of loss of introns, alternative splicing and loss of the last exon coding redox regulated C-terminal extension domain. Gene specific expression analysis, both at the transcription as well as the protein level, showed very different expression pattern of the two RCA copies. Expression of the highly conserved copy was higher under normal plant growing conditions that decreased many folds under heat stress with substantial genotypic variation, but the variable copy showed much higher expression under heat stress conditions across all grass species. The cultivated rice has only one functional gene as the second copy became nonfunctional due to multiple deletions but Oryza brachyantha and Oryza australiensis still have two functional Rca genes. Detailed analysis of the promoter region of the two copies among various plant species showed insertion of several transposable elements harboring heat responsive elements in the heat inducible copy of the gene. The conserved RCA copy of wheat didn't have any transposable insertions whereas in that of maize has one heat shock element and sorghum had two. It would be interesting to study if the higher level of heat stress tolerance observed in sorghum and maize is associated with the differences observed for RCA. Key message This manuscript is reporting a grass family-specific tandem duplication event in RCA genes of cereals. The duplicated copies underwent neo-functionalization to evolve novel function to deal with heat stress. One copy of the tandem duplication maintained a high level of conservation whereas the second copy showed tremendous divergence to evolve species specific function of the gene. Specific function to respond to heat stress likely evolved via the insertion of various heat responsive elements carried by transposable elements.

Kazakhstan tulips: comparative analysis of complete chloroplast genomes of four local and endangered species of the genus Tulipa L

Species of Tulipa are important ornamental plants used for horticultural purposes in various countries, across Asia, Europe, and North Africa. The present study is the first report on typical features of the complete chloroplast genome sequence of four local and endangered species including T. alberti, T. kaufmanniana, T. greigii, and T. dubia from Kazakhstan using Illumina sequencing technology. The comparative analyses revealed that the complete genomes of four species were highly conserved in terms of total genome size (152. 006 bp - 152. 382 bp), including a pair of inverted repeat regions (26. 330 bp - 26. 371 bp), separated by a large single copy region (82.169 bp - 82,378 bp) and a small copy region (17.172 bp -17.260 bp). Total GC content (36.58-36.62 %), gene number (131), and intron length (540 bp - 2620 bp) of 28 genes. The complete genomes of four species showed nucleotide diversity (π =0,003257). The total number of SSR loci was 159 in T. alberti, 158 in T. kaufmanniana, 174 in T. greigii, and 163 in T. dubia. The result indicated that ten CDS genes, namely rpoC2, cemA, rbcL, rpl36, psbH, rps3, rpl22, ndhF, ycf1, and matK, with effective polymorphic simple sequence repeats (SSRs), high sequence variability (SV) ranging from 2.581 to 6.102, and high nucleotide diversity (Pi) of these loci ranging from 0,004 to 0,010. For all intergenic regions longer than 150 bp, twenty one most variable regions were found with high sequence variability (SV) ranging from 4,848 to 11,862 and high nucleotide diversity (Pi) ranging from 0,01599 to 0,01839. Relative synonymous codon usage (RSCU) analysis was used to identify overrepresented and underrepresented codons for each amino acid. Based on the phylogenic analysis, the sequences clustered into four major groups, reflecting distinct evolutionary lineages corresponding to the subgenera Eriostemons, Tulipa, and Orithyia. Notably, T. greigii was distinctively grouped with species from Orithyia and Eriostemons rather than with other Tulipa species, suggesting a unique evolutionary history potentially shaped by geographical isolation or specific ecological pressures. The complete chloroplast genome of the four Tulipa species provides fundamental information for future research studies, even for designing the high number of available molecular markers.

The complete chloroplast genome sequence and phylogenetic analysis of Hordeum pusillum Nutt., 1818 (Poaceae)

Hordeum pusillum (Hordeum pusillum Nutt., 1818) is an annual barley that is native to the west of the North America and widespread in southern United States and tropical America. In this study, we have provided the first complete characterization of the chloroplast genome of H. pusillum. Our research revealed that the circular chloroplast genome of H. pusillum consists of a large single-copy region (LSC: 80,853 bp), a small single-copy region (SSC: 12,745 bp), and a pair of inverted repeat regions (IRs: 21,587 bp), totaling 136,772 bp in length. Within the chloroplast genome of H. pusillum, 91 protein-coding genes, 38 tRNA genes, and ten rRNA genes were identified. To determine the evolutionary relationship of Hordeum species with reported chloroplast genome sequences, we constructed a phylogenetic tree using the entire chloroplast genome sequences. The evolutionary position of H. pusillum corresponds to its geographical location. The chloroplast genome of H. pusillum provided in this study may have significant implications for the phylogenetic study of Poaceae species.

Chlorophyll Fluorescence in Wheat Breeding for Heat and Drought Tolerance

The constantly growing need to increase the production of agricultural products in changing climatic conditions makes it necessary to accelerate the development of new cultivars that meet the modern demands of agronomists. Currently, the breeding process includes the stages of genotyping and phenotyping to optimize the selection of promising genotypes. One of the most popular phenotypic methods is the pulse-amplitude modulated (PAM) fluorometry, due to its non-invasiveness and high information content. In this review, we focused on the opportunities of using chlorophyll fluorescence (ChlF) parameters recorded using PAM fluorometry to assess the state of plants in drought and heat stress conditions and predict the economically significant traits of wheat, as one of the most important agricultural crops, and also analyzed the relationship between the ChlF parameters and genetic markers.

Evolution of Einkorn wheat centromeres is driven by the mutualistic interplay of two LTR retrotransposons

BACKGROUND

Centromere function is highly conserved across eukaryotes, but the underlying centromeric DNA sequences vary dramatically between species. Centromeres often contain a high proportion of repetitive DNA, such as tandem repeats and/or transposable elements (TEs). Einkorn wheat centromeres lack tandem repeat arrays and are instead composed mostly of the two long terminal repeat (LTR) retrotransposon families RLG_Cereba and RLG_Quinta which specifically insert in centromeres. However, it is poorly understood how these two TE families relate to each other and if and how they contribute to centromere function and evolution.

RESULTS

Based on conservation of diagnostic motifs (LTRs, integrase and primer binding site and polypurine-tract), we propose that RLG_Cereba and RLG_Quinta are a pair of autonomous and non-autonomous partners, in which the autonomous RLG_Cereba contributes all the proteins required for transposition, while the non-autonomous RLG_Quinta contributes GAG protein. Phylogenetic analysis of predicted GAG proteins showed that the RLG_Cereba lineage was present for at least 100 million years in monocotyledon plants. In contrast, RLG_Quinta evolved from RLG_Cereba between 28 and 35 million years ago in the common ancestor of oat and wheat. Interestingly, the integrase of RLG_Cereba is fused to a so-called CR-domain, which is hypothesized to guide the integrase to the functional centromere. Indeed, ChIP-seq data and TE population analysis show only the youngest subfamilies of RLG_Cereba and RLG_Quinta are found in the active centromeres. Importantly, the LTRs of RLG_Quinta and RLG_Cereba are strongly associated with the presence of the centromere-specific CENH3 histone variant. We hypothesize that the LTRs of RLG_Cereba and RLG_Quinta contribute to wheat centromere integrity by phasing and/or placing CENH3 nucleosomes, thus favoring their persistence in the competitive centromere-niche.

CONCLUSION

Our data show that RLG_Cereba cross-mobilizes the non-autonomous RLG_Quinta retrotransposons. New copies of both families are specifically integrated into functional centromeres presumably through direct binding of the integrase CR domain to CENH3 histone variants. The LTRs of newly inserted RLG_Cereba and RLG_Quinta elements, in turn, recruit and/or phase new CENH3 deposition. This mutualistic interplay between the two TE families and the plant host dynamically maintains wheat centromeres.

Genetic erosion in domesticated barley and a hypothesis of a North African centre of diversity

Barley is one of the founder crops of the Neolithic transition in West Asia. While recent advances in genomics have provided a rather detailed picture of barley domestication, there are contradictory views on how the domestication process affected genetic diversity. We set out to revisit this question by integrating public DNA sequencing data from ancient barley and wide collections of extant wild and domesticated accessions. Using two previously overlooked approaches - analyses of chloroplast genomes and genome-wide proportions of private variants - we found that the barley cultivated six millennia ago was genetically unique and more diverse when compared to extant landraces and cultivars. Moreover, the chloroplast genomes revealed a link between the ancient barley, an obscure wild genotype from north-eastern Libya, and a distinct population of barley cultivated in Ethiopia/Eritrea. Based on these results, we hypothesize past existence of a wider North African population that included both wild and cultivated types and suffered from genetic erosion in the past six millennia, likely due to a rapid desertification that ended the Holocene African humid period. Besides providing clues about the origin of Ethiopian landraces, the hypothesis explains the post-domestication loss of diversity observed in barley. Analyses of additional samples will be necessary to resolve the history of African barley and its contribution to the extant cultivated gene pool.

Genome constitution and evolution of Elymus atratus (Poaceae: Triticeae) inferred from cytogenetic and phylogenetic analysis

BACKGROUND

Elymus atratus (Nevski) Hand.-Mazz. is perennial hexaploid wheatgrass. It was assigned to the genus Elymus L. sensu stricto based on morphological characters. Its genome constitution has not been disentangled yet.

OBJECTIVE

To identify the genome constitution and origin of E. atratus.

METHODS

In this study, genomic in situ hybridization and fluorescence in situ hybridization, and phylogenetic analysis based on the Acc1, DMC1 and matK sequences were performed.

RESULTS

Genomic in situ hybridization and fluorescence in situ hybridization results reveal that E. atratus 2n = 6x = 42 is composed of 14 St genome chromosomes, 14 H genome chromosomes, and 14 Y genome chromosomes including two H-Y type translocation chromosomes, suggesting that the genome formula of E. atratus is StStYYHH. The phylogenetic analysis based on Acc1 and DMC1 sequences not only shows that the Y genome originated in a separate diploid, but also suggests that Pseudoroegneria (St), Hordeum (H), and a diploid species with Y genome were the potential donors of E. atratus. Data from chloroplast DNA showed that the maternal donor of E. atratus contains the St genome.

CONCLUSION

Elymus atratus is an allohexaploid species with StYH genome, which may have originated through the hybridization between an allotetraploid Roegneria (StY) species as the maternal donor and a diploid Hordeum (H) species as the paternal donor.

Bread wheat satellitome: a complex scenario in a huge genome

In bread wheat (Triticum aestivum L.), chromosome associations during meiosis are extremely regulated and initiate at the telomeres and subtelomeres, which are enriched in satellite DNA (satDNA). We present the study and characterization of the bread wheat satellitome to shed light on the molecular organization of wheat subtelomeres. Our results revealed that the 2.53% of bread wheat genome is composed by satDNA and subtelomeres are particularly enriched in such DNA sequences. Thirty-four satellite DNA (21 for the first time in this work) have been identified, analyzed and cytogenetically validated. Many of the satDNAs were specifically found at particular subtelomeric chromosome regions revealing the asymmetry in subtelomere organisation among the wheat subgenomes, which might play a role in proper homologous recognition and pairing during meiosis. An integrated physical map of the wheat satellitome was also constructed. To the best of our knowledge, our results show that the combination of both cytogenetics and genome research allowed the first comprehensive analysis of the wheat satellitome, shedding light on the complex wheat genome organization, especially on the polymorphic nature of subtelomeres and their putative implication in chromosome recognition and pairing during meiosis.

The complete chloroplast genome of Hordeum marinum ssp. marinum

Hordeum marinum ssp. marinum (Huds.) R. J. Soreng (2003) is a halophyte wild relative of barley and wheat, which exhibits remarkable salt tolerance characteristics. In this study, we presented the first characterization of the complete chloroplast genome of H. marinum ssp. marinum. Our findings reveal that the chloroplast genome of H. marinum ssp. marinum consists of a small single-copy region (SSC: 12,715 bp), a large single-copy region (LSC: 81,130 bp), and a pair of inverted repeat regions (IRs: 21,517 bp), amounting to a total length of 136,989 bp. The chloroplast genome encodes 139 genes, including 91 protein-coding genes, 38 tRNA genes, and ten rRNA genes. By utilizing phylogenetic analysis, we determine the evolutionary position of H. marinum in Triticeae. Our study provides valuable insights into the chloroplast genome of H. marinum ssp. marinum, which may have important implications for the improvement of cereal crops.

Crops under past diversification and ongoing climate change: more than just selection of nuclear genes for flowering

Diversification and breeding following domestication and under current climate change across the globe are the two most significant evolutionary events experienced by major crops. Diversification of crops from their wild ancestors has favored dramatic changes in the sensitivity of the plants to the environment, particularly significantly in transducing light inputs to the circadian clock, which has allowed the growth of major crops in the relatively short growing season experienced in the Northern Hemisphere. Historically, mutants and the mapping of quantitative trait loci (QTL) have facilitated the identification and the cloning of genes that underlie major changes of the clock and the regulation of flowering. Recent studies have suggested that the thermal plasticity of the circadian clock output, and not just the core genes that follow temperature compensation, has also been under selection during diversification and breeding. Wild alleles that accelerate output rhythmicity could be beneficial for crop resilience. Furthermore, wild alleles with beneficial and flowering-independent effects under stress indicate their possible role in maintaining a balanced source-sink relationship, thereby allowing productivity under climatic change. Because the chloroplast genome also regulates the plasticity of the clock output, mapping populations including cytonuclear interactions should be utilized within an integrated field and clock phenomics framework. In this review, we highlight the need to integrate physiological and developmental approaches (physio-devo) to gain a better understanding when re-domesticating wild gene alleles into modern cultivars to increase their robustness under abiotic heat and drought stresses.

Advancing the Conservation and Utilization of Barley Genetic Resources: Insights into Germplasm Management and Breeding for Sustainable Agriculture

Barley is a very important crop particularly in marginal dry areas, where it often serves as the most viable option for farmers. Additionally, barley carries great significance in the Western world, serving not only as a fundamental crop for animal feed and malting but also as a nutritious food source. The broad adaptability of barley and its ability to withstand various biotic and abiotic stresses often make this species the sole cereal that can be cultivated in arid regions. The collection and utilization of barley genetic resources are crucial for identifying valuable traits to enhance productivity and mitigate the adverse effects of climate change. This review aims to provide an overview of the management and exploitation of barley genetic resources. Furthermore, the review explores the relationship between gene banks and participatory breeding, offering insights into the diversity and utilization of barley genetic resources through some examples such as the initiatives undertaken by ICARDA. Finally, this contribution highlights the importance of these resources for boosting barley productivity, addressing climate change impacts, and meeting the growing food demands in a rapidly changing agriculture. The understanding and utilizing the rich genetic diversity of barley can contribute to sustainable agriculture and ensure the success of this vital crop for future generations globally.

All families of transposable elements were active in the recent wheat genome evolution and polyploidy had no impact on their activity

Bread wheat (Triticum aestivum L.) is a major crop and its genome is one of the largest ever assembled at reference-quality level. It is 15 Gb, hexaploid, with 85% of transposable elements (TEs). Wheat genetic diversity was mainly focused on genes and little is known about the extent of genomic variability affecting TEs, transposition rate, and the impact of polyploidy. Multiple chromosome-scale assemblies are now available for bread wheat and for its tetraploid and diploid wild relatives. In this study, we computed base pair-resolved, gene-anchored, whole genome alignments of A, B, and D lineages at different ploidy levels in order to estimate the variability that affects the TE space. We used assembled genomes of 13 T. aestivum cultivars (6x = AABBDD) and a single genome for Triticum durum (4x = AABB), Triticum dicoccoides (4x = AABB), Triticum urartu (2x = AA), and Aegilops tauschii (2x = DD). We show that 5%-34% of the TE fraction is variable, depending on the species divergence. Between 400 and 13,000 novel TE insertions per subgenome were detected. We found lineage-specific insertions for nearly all TE families in di-, tetra-, and hexaploids. No burst of transposition was observed and polyploidization did not trigger any boost of transposition. This study challenges the prevailing idea of wheat TE dynamics and is more in agreement with an equilibrium model of evolution.

Complete chloroplast genomes of three wild perennial Hordeum species from Central Asia: genome structure, mutation hotspot, phylogenetic relationships, and comparative analysis

Hordeum L. is widely distributed in mountain or plateau of subtropical and warm temperate regions around the world. Three wild perennial Hordeum species, including H. bogdanii, H. brevisubulatum, and H. violaceum, have been used as forage and for grassland ecological restoration in high-altitude areas in recent years. To date, the degree of interspecies sequence variation in the three Hordeum species within existing gene pools is still not well-defined. Herein, we sequenced and assembled chloroplast (cp) genomes of the three species. The results revealed that the cp genome of H. bogdanii showed certain sequence variations compared with the cp genomes of the other two species (H. brevisubulatum and H. violaceum), and the latter two were characterized by a higher relative affinity. Parity rule 2 plot (PR2) analysis illuminated that most genes of all ten Hordeum species were concentrated in nucleotide T and G. Numerous single nucleotide polymorphism (SNP) and insertion/deletion (In/Del) events were detected in the three Hordeum species. A series of hotspots regions (tRNA-GGU ~ tRNA-GCA, tRNA-UGU ~ ndhJ, psbE ~ rps18, ndhF ~ tRNA-UAG, etc.) were identified by mVISTA procedures, and the five highly polymorphic genes (tRNA-UGC, tRNA-UAA, tRNA-UUU, tRNA-UAC, and ndhA) were proved by the nucleotide diversity (Pi). Although the distribution and existence of cp simple sequence repeats (cpSSRs) were predicted in the three Hordeum cp genomes, no rearrangement was found between them. A similar phenomenon has been found in the cp genome of the other seven Hordeum species, which has been published so far. In addition, evolutionary relationships were reappraised based on the currently reported cp genome of Hordeum L. This study offers a framework for gaining a better understanding of the evolutionary history of Hordeum species through the re-examination of their cp genomes, and by identifying highly polymorphic genes and hotspot regions that could provide important insights into the genetic diversity and differentiation of these species.

Comparative genomic analysis of 5Mg chromosome of Aegilops geniculata and 5Uu chromosome of Aegilops umbellulata reveal genic diversity in the tertiary gene pool

Wheat is one of the most important cereal crops for the global food security. Due to its narrow genetic base, modern bread wheat cultivars face challenges from increasing abiotic and biotic stresses. Since genetic improvement is the most sustainable approach, finding novel genes and alleles is critical for enhancing the genetic diversity of wheat. The tertiary gene pool of wheat is considered a gold mine for genetic diversity as novel genes and alleles can be identified and transferred to wheat cultivars. Aegilops geniculata and Ae. umbellulata are the key members of the tertiary gene pool of wheat and harbor important genes against abiotic and biotic stresses. Homoeologous-group five chromosomes (5U^u and 5M^g) have been extensively studied from Ae. geniculata and Ae. umbellulata as they harbor several important genes including Lr57, Lr76, Yr40, Yr70, Sr53 and chromosomal pairing loci. In the present study, using chromosome DNA sequencing and RNAseq datasets, we performed comparative analysis to study homoeologous gene evolution in 5M^g, 5U^u, and group 5 wheat chromosomes. Our findings highlight the diversity of transcription factors and resistance genes, resulting from the differential expansion of the gene families. Both the chromosomes were found to be enriched with the "response to stimulus" category of genes providing resistance against biotic and abiotic stress. Phylogenetic study positioned the M genome closer to the D genome, with higher proximity to the A genome than the B genome. Over 4000 genes were impacted by SNPs on 5D, with 4-5% of those genes displaying non-disruptive variations that affect gene function.

Chromosome-level genome and high nitrogen stress response of the widespread and ecologically important wetland plant Typha angustifolia

Typha angustifolia L., known as narrowleaf cattail, is widely distributed in Eurasia but has been introduced to North America. Typha angustifolia is a semi-aquatic, wetland obligate plant that is widely distributed in Eurasia and North America. It is ecologically important for nutrient cycling in wetlands where it occurs and is used in phytoremediation and traditional medicine. In order to construct a high-quality genome for Typha angustifolia and investigate genes in response to high nitrogen stress, we carried out complete genome sequencing and high-nitrogen-stress experiments. We generated a chromosomal-level genome of T. angustifolia, which had 15 pseudochromosomes, a size of 207 Mb, and a contig N50 length of 13.57 Mb. Genome duplication analyses detected no recent whole-genome duplication (WGD) event for T. angustifolia. An analysis of gene family expansion and contraction showed that T. angustifolia gained 1,310 genes and lost 1,426 genes. High-nitrogen-stress experiments showed that a high nitrogen level had a significant inhibitory effect on root growth and differential gene expression analyses using 24 samples found 128 differentially expressed genes (DEGs) between the nitrogen-treated and control groups. DEGs in the roots and leaves were enriched in alanines, aspartate, and glutamate metabolism, nitrogen metabolism, photosynthesis, phenylpropanoid biosynthesis, plant-pathogen interaction, and mitogen-activated protein kinase pathways, among others. This study provides genomic data for a medicinal and ecologically important herb and lays a theoretical foundation for plant-assisted water pollution remediation.

Molecular structure, comparative and phylogenetic analysis of the complete chloroplast genome sequences of weedy rye Secale cereale ssp. segetale

The complete chloroplast genome of Secale cereale ssp. segetale (Zhuk.) Roshev. (Poaceae: Triticeae) was sequenced and analyzed to better use its genetic resources to enrich rye and wheat breeding. The study was carried out using the following methods: DNA extraction, sequencing, assembly and annotation, comparison with other complete chloroplast genomes of the five Secale species, and multigene phylogeny. As a result of the study, it was determined that the chloroplast genome is 137,042 base pair (bp) long and contains 137 genes, including 113 unique genes and 24 genes which are duplicated in the IRs. Moreover, a total of 29 SSRs were detected in the Secale cereale ssp. segetale chloroplast genome. The phylogenetic analysis showed that Secale cereale ssp. segetale appeared to share the highest degree of similarity with S. cereale and S. strictum. Intraspecific diversity has been observed between the published chloroplast genome sequences of S. cereale ssp. segetale. The genome can be accessed on GenBank with the accession number (OL688773).

Grain Disarticulation in Wild Wheat and Barley

Our industrial-scale crop monocultures, which are necessary to provide grain for large-scale food and feed production, are highly vulnerable to biotic and abiotic stresses. Crop wild relatives have adapted to harsh environmental conditions over millennia; thus, they are an important source of genetic variation and crop diversification. Despite several examples where significant yield increases have been achieved through the introgression of genomic regions from wild relatives, more detailed understanding of the differences between wild and cultivated species for favorable and unfavorable traits is still required to harness these valuable resources. Recently, as an alternative to the introgression of beneficial alleles from the wild into domesticated species, a radical suggestion is to domesticate wild relatives to generate new crops. A first and critical step for the domestication of cereal wild relatives would be to prevent grain disarticulation from the inflorescence at maturity. Discovering the molecular mechanisms and understanding the network of interactions behind grain retention/disarticulation would enable the implementation of approaches to select for this character in targeted species. Brittle rachis 1 and Brittle rachis 2 are major genes responsible for grain disarticulation in the wild progenitors of wheat and barley that were the target of mutations during domestication. These two genes are only found in the Triticeae tribe and are hypothesized to have evolved by a duplication followed by neo-functionalization. Current knowledge gaps include the molecular mechanisms controlling grain retention in cereals and the genomic consequences of strong selection for this essential character.

Complete Chloroplast Genome Sequence of the Endemic and Endangered Plant Dendropanax oligodontus: Genome Structure, Comparative and Phylogenetic Analysis

Dendropanax oligodontus, which belongs to the family Araliaceae, is an endemic and endangered species of Hainan Island, China. It has potential economic and medicinal value owing to the presence of phenylpropanoids, flavonoids, triterpenoids, etc. The analysis of the structure and characteristics of the D. oligodontus chloroplast genome (cpDNA) is crucial for understanding the genetic and phylogenetic evolution of this species. In this study, the cpDNA of D. oligodontus was sequenced for the first time using next-generation sequencing methods, assembled, and annotated. We observed a circular quadripartite structure comprising a large single-copy region (86,440 bp), a small single-copy region (18,075 bp), and a pair of inverted repeat regions (25,944 bp). The total length of the cpDNA was 156,403 bp, and the GC% was 37.99%. We found that the D. oligodontus chloroplast genome comprised 131 genes, with 86 protein-coding genes, 8 rRNA genes, and 37 tRNAs. Furthermore, we identified 26,514 codons, 13 repetitive sequences, and 43 simple sequence repeat sites in the D. oligodontus cpDNA. The most common amino acid encoded was leucine, with a strong A/T preference at the third position of the codon. The prediction of RNA editing sites in the protein-coding genes indicated that RNA editing was observed in 19 genes with a total of 54 editing sites, all of which involved C-to-T transitions. Finally, the cpDNA of 11 species of the family Araliaceae were selected for comparative analysis. The sequences of the untranslated regions and coding regions among 11 species were highly conserved, and minor differences were observed in the length of the inverted repeat regions; therefore, the cpDNAs were relatively stable and consistent among these 11 species. The variable hotspots in the genome included clpP, ycf1, rnK-rps16, rps16-trnQ, atpH-atpI, trnE-trnT, psbM-trnD, ycf3-trnS, and rpl32-trnL, providing valuable molecular markers for species authentication and regions for inferring phylogenetic relationships among them, as well as for evolutionary studies. Evolutionary selection pressure analysis indicated that the atpF gene was strongly subjected to positive environmental selection. Phylogenetic analysis indicated that D. oligodontus and Dendropanax dentiger were the most closely related species within the genus, and D. oligodontus was closely related to the genera Kalopanax and Metapanax in the Araliaceae family. Overall, the cp genomes reported in this study will provide resources for studying the genetic diversity and conservation of the endangered plant D. oligodontus, as well as resolving phylogenetic relationships within the family.

Improving Yield and Yield Stability in Winter Rye by Hybrid Breeding

Rye is the only cross-pollinating small-grain cereal. The unique reproduction biology results in an exceptional complexity concerning genetic improvement of rye by breeding. Rye is a close relative of wheat and has a strong adaptation potential that refers to its mating system, making this overlooked cereal readily adjustable to a changing environment. Rye breeding addresses the emerging challenges of food security associated with climate change. The systematic identification, management, and use of its valuable natural diversity became a feasible option in outbreeding rye only following the establishment of hybrid breeding late in the 20th century. In this article, we review the most recent technological advances to improve yield and yield stability in winter rye. Based on recently released reference genome sequences, SMART breeding approaches are described to counterbalance undesired linkage drag effects of major restorer genes on grain yield. We present the development of gibberellin-sensitive semidwarf hybrids as a novel plant breeding innovation based on an approach that is different from current methods of increasing productivity in rye and wheat. Breeding of new rye cultivars with improved performance and resilience is indispensable for a renaissance of this healthy minor cereal as a homogeneous commodity with cultural relevance in Europe that allows for comparatively smooth but substantial complementation of wheat with rye-based diets, supporting the necessary restoration of the balance between human action and nature.

Grain dispersal mechanism in cereals arose from a genome duplication followed by changes in spatial expression of genes involved in pollen development

KEY MESSAGE

Grain disarticulation in wild progenitor of wheat and barley evolved through a local duplication event followed by neo-functionalization resulting from changes in location of gene expression. One of the most critical events in the process of cereal domestication was the loss of the natural mode of grain dispersal. Grain dispersal in barley is controlled by two major genes, Btr1 and Btr2, which affect the thickness of cell walls around the disarticulation zone. The barley genome also encodes Btr1-like and Btr2-like genes, which have been shown to be the ancestral copies. While Btr and Btr-like genes are non-redundant, the biological function of Btr-like genes is unknown. We explored the potential biological role of the Btr-like genes by surveying their expression profile across 212 publicly available transcriptome datasets representing diverse organs, developmental stages and stress conditions. We found that Btr1-like and Btr2-like are expressed exclusively in immature anther samples throughout Prophase I of meiosis within the meiocyte. The similar and restricted expression profile of these two genes suggests they are involved in a common biological function. Further analysis revealed 141 genes co-expressed with Btr1-like and 122 genes co-expressed with Btr2-like, with 105 genes in common, supporting Btr-like genes involvement in a shared molecular pathway. We hypothesize that the Btr-like genes play a crucial role in pollen development by facilitating the formation of the callose wall around the meiocyte or in the secretion of callase by the tapetum. Our data suggest that Btr genes retained an ancestral function in cell wall modification and gained a new role in grain dispersal due to changes in their spatial expression becoming spike specific after gene duplication.

Genetic characterization and curation of diploid A-genome wheat species

A-genome diploid wheats represent the earliest domesticated and cultivated wheat species in the Fertile Crescent and include the donor of the wheat A sub-genome. The A-genome species encompass the cultivated einkorn (Triticum monococcum L. subsp. monococcum), wild einkorn (T. monococcum L. subsp. aegilopoides (Link) Thell.), and Triticum urartu. We evaluated the collection of 930 accessions in the Wheat Genetics Resource Center (WGRC) using genotyping by sequencing and identified 13,860 curated single-nucleotide polymorphisms. Genomic analysis detected misclassified and genetically identical (>99%) accessions, with most of the identical accessions originating from the same or nearby locations. About 56% (n = 520) of the WGRC A-genome species collections were genetically identical, supporting the need for genomic characterization for effective curation and maintenance of these collections. Population structure analysis confirmed the morphology-based classifications of the accessions and reflected the species geographic distributions. We also showed that T. urartu is the closest A-genome diploid to the A-subgenome in common wheat (Triticum aestivum L.) through phylogenetic analysis. Population analysis within the wild einkorn group showed three genetically distinct clusters, which corresponded with wild einkorn races α, β, and γ described previously. The T. monococcum genome-wide FST scan identified candidate genomic regions harboring a domestication selection signature at the Non-brittle rachis 1 (Btr1) locus on the short arm of chromosome 3Am at ∼70 Mb. We established an A-genome core set (79 accessions) based on allelic diversity, geographical distribution, and available phenotypic data. The individual species core set maintained at least 79% of allelic variants in the A-genome collection and constituted a valuable genetic resource to improve wheat and domesticated einkorn in breeding programs.

The Plastome Sequences of Triticum sphaerococcum (ABD) and Triticum turgidum subsp. durum (AB) Exhibit Evolutionary Changes, Structural Characterization, Comparative Analysis, Phylogenomics and Time Divergence

The mechanism and course of Triticum plastome evolution is currently unknown; thus, it remains unclear how Triticum plastomes evolved during recent polyploidization. Here, we report the complete plastomes of two polyploid wheat species, Triticum sphaerococcum (AABBDD) and Triticum turgidum subsp. durum (AABB), and compare them with 19 available and complete Triticum plastomes to create the first map of genomic structural variation. Both T. sphaerococcum and T. turgidum subsp. durum plastomes were found to have a quadripartite structure, with plastome lengths of 134,531 bp and 134,015 bp, respectively. Furthermore, diploid (AA), tetraploid (AB, AG) and hexaploid (ABD, AGA^m) Triticum species plastomes displayed a conserved gene content and commonly harbored an identical set of annotated unique genes. Overall, there was a positive correlation between the number of repeats and plastome size. In all plastomes, the number of tandem repeats was higher than the number of palindromic and forward repeats. We constructed a Triticum phylogeny based on the complete plastomes and 42 shared genes from 71 plastomes. We estimated the divergence of Hordeum vulgare from wheat around 11.04-11.9 million years ago (mya) using a well-resolved plastome tree. Similarly, Sitopsis species diverged 2.8-2.9 mya before Triticum urartu (AA) and Triticum monococcum (AA). Aegilops speltoides was shown to be the maternal donor of polyploid wheat genomes and diverged ~0.2-0.9 mya. The phylogeny and divergence time estimates presented here can act as a reference framework for future studies of Triticum evolution.

Transposable Element Populations Shed Light on the Evolutionary History of Wheat and the Complex Co-Evolution of Autonomous and Non-Autonomous Retrotransposons

Wheat has one of the largest and most repetitive genomes among major crop plants, containing over 85% transposable elements (TEs). TEs populate genomes much in the way that individuals populate ecosystems, diversifying into different lineages, sub-families and sub-populations. The recent availability of high-quality, chromosome-scale genome sequences from ten wheat lines enables a detailed analysis how TEs evolved in allohexaploid wheat, its diploids progenitors, and in various chromosomal haplotype segments. LTR retrotransposon families evolved into distinct sub-populations and sub-families that were active in waves lasting several hundred thousand years. Furthermore, It is shown that different retrotransposon sub-families were active in the three wheat sub-genomes, making them useful markers to study and date polyploidization events and chromosomal rearrangements. Additionally, haplotype-specific TE sub-families are used to characterize chromosomal introgressions in different wheat lines. Additionally, populations of non-autonomous TEs co-evolved over millions of years with their autonomous partners, leading to complex systems with multiple types of autonomous, semi-autonomous and non-autonomous elements. Phylogenetic and TE population analyses revealed the relationships between non-autonomous elements and their mobilizing autonomous partners. TE population analysis provided insights into genome evolution of allohexaploid wheat and genetic diversity of species, and may have implication for future crop breeding.

Fine structure and transcription dynamics of bread wheat ribosomal DNA loci deciphered by a multi-omics approach

Three out of four RNA components of ribosomes are encoded by 45S ribosomal DNA (rDNA) loci, which are organized as long head-to-tail tandem arrays of nearly identical units, spanning several megabases of sequence. Due to this structure, the rDNA loci are the major sources of gaps in genome assemblies, and gene copy number, sequence composition, and expression status of particular arrays remain elusive, especially in complex genomes harboring multiple loci. Here we conducted a multi-omics study to decipher the 45S rDNA loci in hexaploid bread wheat. Coupling chromosomal genomics with optical mapping, we reconstructed individual rDNA arrays, enabling locus-specific analyses of transcription activity and methylation status from RNA- and bisulfite-sequencing data. We estimated a total of 6,650 rDNA units in the bread wheat genome, with approximately 2,321, 3,910, 253, and 50 gene copies located in short arms of chromosomes 1B, 6B, 5D, and 1A, respectively. Only 1B and 6B loci contributed substantially to rRNA transcription at a roughly 2:1 ratio. The ratio varied among five tissues analyzed (embryo, coleoptile, root tip, primary leaf, mature leaf), being the highest (2.64:1) in mature leaf and lowest (1.72:1) in coleoptile. Cytosine methylation was considerably higher in CHG context in the silenced 5D locus as compared with the active 1B and 6B loci. In conclusion, a fine genomic organization and tissue-specific expression of rDNA loci were deciphered, for the first time, in a complex polyploid species. The results are discussed in the context of wheat evolution and transcription regulation.

Chloroplast Genes Are Involved in The Male-Sterility of K-Type CMS in Wheat

The utilization of crop heterosis can greatly improve crop yield. The sterile line is vital for the heterosis utilization of wheat (Triticum aestivum L.). The chloroplast genomes of two sterile lines and one maintainer were sequenced using second-generation high-throughput technology and assembled. The nonsynonymous mutated genes among the three varieties were identified, the expressed difference was further analyzed by qPCR, and finally, the function of the differentially expressed genes was analyzed by the barley stripe mosaic virus-induced gene silencing (BSMV-VIGS) method. A total of 16 genes containing 31 nonsynonymous mutations between K519A and 519B were identified. There were no base mutations in the protein-encoding genes between K519A and YS3038. The chloroplast genomes of 519B and K519A were closely related to the Triticum genus and Aegilops genus, respectively. The gene expression levels of the six selected genes with nonsynonymous mutation sites for K519A compared to 519B were mostly downregulated at the binucleate and trinucleate stages of pollen development. The seed setting rates of atpB-silenced or ndhH-silenced 519B plants by BSMV-VIGS method were significantly reduced. It can be concluded that atpB and the ndhH are likely to be involved in the reproductive transformation of 519B.

Analysis of Chloroplast Genomes Provides Insights Into the Evolution of Agropyron

Plants of the Agropyron genus are important pasture resources, and they also play important roles in the ecological restoration. Chloroplast genomes are inherited from maternal parents, and they are important for studying species taxonomy and evolution. In this study, we sequenced the complete chloroplast genomes of five typical species of the Agropyron genus (eg., A. cristatum × A. desertorum Fisch. Schult, A. desertorum, A. desertorum Fisch. Schult. cv. Nordan, A. michnoi Roshev, and A. mongolicum Keng) using the Illumina NovaSeq platform. We found that these five chloroplast genomes exhibit a typical quadripartite structure with a conserved genome arrangement and structure. Their chloroplast genomes contain the large single-copy regions (LSC, 79,613 bp-79,634 bp), the small single-copy regions (SSC, 12,760 bp-12,768 bp), and the inverted repeat regions (IR, 43,060 bp-43,090 bp). Each of the five chloroplast genomes contains 129 genes, including 38 tRNA genes, eight rRNA genes, and 83 protein-coding genes. Among them, the genes trnG-GCC, matK, petL, ccsA, and rpl32 showed significant nucleotide diversity in these five species, and they may be used as molecular markers in taxonomic studies. Phylogenetic analysis showed that A. mongolicum is closely related to A. michnoi, while others have a closer genetic relationship with the Triticum genus.

Chloroplast phylogenomics in Camelina (Brassicaceae) reveals multiple origins of polyploid species and the maternal lineage of C. sativa

The genus Camelina (Brassicaceae) comprises 7-8 diploid, tetraploid, and hexaploid species. Of particular agricultural interest is the biofuel crop, C. sativa (gold-of-pleasure or false flax), an allohexaploid domesticated from the widespread weed, C. microcarpa. Recent cytogenetics and genomics work has uncovered the identity of the parental diploid species involved in ancient polyploidization events in Camelina. However, little is known about the maternal subgenome ancestry of contemporary polyploid species. To determine the diploid maternal contributors of polyploid Camelina lineages, we sequenced and assembled 84 Camelina chloroplast genomes for phylogenetic analysis. Divergence time estimation was used to infer the timing of polyploidization events. Chromosome counts were also determined for 82 individuals to assess ploidy and cytotypic variation. Chloroplast genomes showed minimal divergence across the genus, with no observed gene-loss or structural variation. Phylogenetic analyses revealed C. hispida as a maternal diploid parent to the allotetraploid Camelina rumelica, and C. neglecta as the closest extant diploid contributor to the allohexaploids C. microcarpa and C. sativa. The tetraploid C. rumelica appears to have evolved through multiple independent hybridization events. Divergence times for polyploid lineages closely related to C. sativa were all inferred to be very recent, at only ~65 thousand years ago. Chromosome counts confirm that there are two distinct cytotypes within C. microcarpa (2n = 38 and 2n = 40). Based on these findings and other recent research, we propose a model of Camelina subgenome relationships representing our current understanding of the hybridization and polyploidization history of this recently-diverged genus.

Comparative Chloroplast Genome Analyses of the Winter-Blooming Eastern Asian Endemic Genus Chimonanthus (Calycanthaceae) With Implications For Its Phylogeny and Diversification

Chimonanthus of Calycanthaceae is a small endemic genus in China, with unusual winter-blooming sweet flowers widely cultivated for ornamentals and medicinal uses. The evolution of Chimonanthus plastomes and its phylogenetic relationships remain unresolved due to limited availability of genetic resources. Here, we report fully assembled and annotated chloroplast genomes of five Chimonanthus species. The chloroplast genomes of the genus (size range 153,010 – 153,299 bp) reveal high similarities in gene content, gene order, GC content, codon usage, amino acid frequency, simple sequence repeats, oligonucleotide repeats, synonymous and non-synonymous substitutions, and transition and transversion substitutions. Signatures of positive selection are detected in atpF and rpoB genes in C. campanulatus. The correlations among substitutions, InDels, and oligonucleotide repeats reveal weak to strong correlations in distantly related species at the intergeneric levels, and very weak to weak correlations among closely related Chimonanthus species. Chloroplast genomes are used to reconstruct a well-resolved phylogenetic tree, which supports the monophyly of Chimonanthus. Within Chimonanthus, C. praecox and C. campanulatus form one clade, while C. grammatus, C. salicifolius, C. zhejiangensis, and C. nitens constitute another clade. Chimonanthus nitens appears paraphyletic and is closely related to C. salicifolius and C. zhejiangensis, suggesting the need to reevaluate the species delimitation of C. nitens. Chimonanthus and Calycanthus diverged in mid-Oligocene; the radiation of extant Chimonanthus species was dated to the mid-Miocene, while C. grammatus diverged from other Chimonanthus species in the late Miocene. C. salicifolius, C. nitens(a), and C. zhejiangensis are inferred to have diverged in the Pleistocene of the Quaternary period, suggesting recent speciation of a relict lineage in the subtropical forest regions in eastern China. This study provides important insights into the chloroplast genome features and evolutionary history of Chimonanthus and family Calycanthaceae.

Chloroplast Phylogenomic Analyses Reveal a Maternal Hybridization Event Leading to the Formation of Cultivated Peanuts

Peanuts (Arachis hypogaea L.) offer numerous healthy benefits, and the production of peanuts has a prominent role in global food security. As a result, it is in the interest of society to improve the productivity and quality of peanuts with transgenic means. However, the lack of a robust phylogeny of cultivated and wild peanut species has limited the utilization of genetic resources in peanut molecular breeding. In this study, a total of 33 complete peanut plastomes were sequenced, analyzed and used for phylogenetic analyses. Our results suggest that sect. Arachis can be subdivided into two lineages. All the cultivated species are contained in Lineage I with AABB and AA are the two predominant genome types present, while species in Lineage II possess diverse genome types, including BB, KK, GG, etc. Phylogenetic studies also indicate that all allotetraploid cultivated peanut species have been derived from a possible maternal hybridization event with one of the diploid Arachis duranensis accessions being a potential AA sub-genome ancestor. In addition, Arachis monticola, a tetraploid wild species, is placed in the same group with all the cultivated peanuts, and it may represent a transitional species, which has been through the recent hybridization event. This research could facilitate a better understanding of the taxonomic status of various Arachis species/accessions and the evolutionary relationship among them, and assists in the correct and efficient use of germplasm resources in breeding efforts to improve peanuts for the benefit of human beings.

The Dynamism of Transposon Methylation for Plant Development and Stress Adaptation

Plant development processes are regulated by epigenetic alterations that shape nuclear structure, gene expression, and phenotypic plasticity; these alterations can provide the plant with protection from environmental stresses. During plant growth and development, these processes play a significant role in regulating gene expression to remodel chromatin structure. These epigenetic alterations are mainly regulated by transposable elements (TEs) whose abundance in plant genomes results in their interaction with genomes. Thus, TEs are the main source of epigenetic changes and form a substantial part of the plant genome. Furthermore, TEs can be activated under stress conditions, and activated elements cause mutagenic effects and substantial genetic variability. This introduces novel gene functions and structural variation in the insertion sites and primarily contributes to epigenetic modifications. Altogether, these modifications indirectly or directly provide the ability to withstand environmental stresses. In recent years, many studies have shown that TE methylation plays a major role in the evolution of the plant genome through epigenetic process that regulate gene imprinting, thereby upholding genome stability. The induced genetic rearrangements and insertions of mobile genetic elements in regions of active euchromatin contribute to genome alteration, leading to genomic stress. These TE-mediated epigenetic modifications lead to phenotypic diversity, genetic variation, and environmental stress tolerance. Thus, TE methylation is essential for plant evolution and stress adaptation, and TEs hold a relevant military position in the plant genome. High-throughput techniques have greatly advanced the understanding of TE-mediated gene expression and its associations with genome methylation and suggest that controlled mobilization of TEs could be used for crop breeding. However, development application in this area has been limited, and an integrated view of TE function and subsequent processes is lacking. In this review, we explore the enormous diversity and likely functions of the TE repertoire in adaptive evolution and discuss some recent examples of how TEs impact gene expression in plant development and stress adaptation.

Introducing Beneficial Alleles from Plant Genetic Resources into the Wheat Germplasm

Wheat (Triticum sp.) is one of the world's most important crops, and constantly increasing its productivity is crucial to the livelihoods of millions of people. However, more than a century of intensive breeding and selection processes have eroded genetic diversity in the elite genepool, making new genetic gains difficult. Therefore, the need to introduce novel genetic diversity into modern wheat has become increasingly important. This review provides an overview of the plant genetic resources (PGR) available for wheat. We describe the most important taxonomic and phylogenetic relationships of these PGR to guide their use in wheat breeding. In addition, we present the status of the use of some of these resources in wheat breeding programs. We propose several introgression schemes that allow the transfer of qualitative and quantitative alleles from PGR into elite germplasm. With this in mind, we propose the use of a stage-gate approach to align the pre-breeding with main breeding programs to meet the needs of breeders, farmers, and end-users. Overall, this review provides a clear starting point to guide the introgression of useful alleles over the next decade.

RNA-Seq-based DNA marker analysis of the genetics and molecular evolution of Triticeae species

The release of high-quality chromosome-level genome sequences of members of the Triticeae tribe has greatly facilitated genetic and genomic analyses of important crops such as wheat (Triticum aestivum) and barley (Hordeum vulgare). Due to the large diploid genome size of Triticeae plants (ca. 5 Gbp), transcript analysis is an important method for identifying genetic and genomic differences among Triticeae species. In this review, we summarize our results of RNA-Seq analyses of diploid wheat accessions belonging to the genera Aegilops and Triticum. We also describe studies of the molecular relationships among these accessions and provide insight into the evolution of common hexaploid wheat. DNA markers based on polymorphisms within species can be used to map loci of interest. Even though the genome sequence of diploid Aegilops tauschii, the D-genome donor of common wheat, has been released, the diploid barley genome continues to provide key information about the physical structures of diploid wheat genomes. We describe how a series of RNA-Seq analyses of wheat relatives has helped uncover the structural and evolutionary features of genomic and genetic systems in wild and cultivated Triticeae species.

Characterizing chloroplast genomes and inferring maternal divergence of the Triticum-Aegilops complex

The Triticum (wheat)-Aegilops (goatgrass) complex has been extensively studied, but the evolutionary history of polyploid wheats has not been fully elucidated. The chloroplast (cp) with maternal inheritance and homoplasy can simplify the sequence-based evolutionary inferences, but informative inferences would require a complete and accurate cp genome sequence. In this study, 16 cp genomes representing five Aegilops and 11 Triticum species and subspecies were sequenced, assembled and annotated, yielding five novel circular cp genome sequences. Analyzing the assembled cp genomes revealed no marked differences in genome structure and gene arrangement across the assayed species. A polymorphism analysis of 72 published cp genome sequences representing 10 Aegilops and 15 Triticum species and subspecies detected 1183 SNPs and 1881 SSRs. More than 80% SNPs detected resided on the downstream and upstream gene regions and only 2.78% or less SNPs were predicted to be deleterious. The largest nucleotide diversity was observed in the short single-copy genomic region. Relatively weak selection pressure on cp coding genes was detected. Different phylogenetic analyses confirmed that the maternal divergence of the Triticum-Aegilops complex had three deep lineages each representing a diploid species with nuclear A, B, or D genome. Dating the maternal divergence yielded age estimates of divergence that matched well with those reported previously. The divergence between emmer and bread wheats occurred at 8200-11,200 years ago. These findings are useful for further genomic studies, provide insight into cp genome evolvability and allow for better understanding of the maternal divergence of the Triticum-Aegilops complex.

Cellular and transcriptomic analyses reveal two-staged chloroplast biogenesis underpinning photosynthesis build-up in the wheat leaf

BACKGROUND

The developmental gradient in monocot leaves has been exploited to uncover leaf developmental gene expression programs and chloroplast biogenesis processes. However, the relationship between the two is barely understood, which limits the value of transcriptome data to understand the process of chloroplast development.

RESULTS

Taking advantage of the developmental gradient in the bread wheat leaf, we provide a simultaneous quantitative analysis for the development of mesophyll cells and of chloroplasts as a cellular compartment. This allows us to generate the first biologically-informed gene expression map of this leaf, with the entire developmental gradient from meristematic to fully differentiated cells captured. We show that the first phase of plastid development begins with organelle proliferation, which extends well beyond cell proliferation, and continues with the establishment and then the build-up of the plastid genetic machinery. The second phase is marked by the development of photosynthetic chloroplasts which occupy the available cellular space. Using a network reconstruction algorithm, we predict that known chloroplast gene expression regulators are differentially involved across those developmental stages.

CONCLUSIONS

Our analysis generates both the first wheat leaf transcriptional map and one of the most comprehensive descriptions to date of the developmental history of chloroplasts in higher plants. It reveals functionally distinct plastid and chloroplast development stages, identifies processes occurring in each of them, and highlights our very limited knowledge of the earliest drivers of plastid biogenesis, while providing a basis for their future identification.

Green giant-a tiny chloroplast genome with mighty power to produce high-value proteins: history and phylogeny

Free-living cyanobacteria were entrapped by eukaryotic cells ~2 billion years ago, ultimately giving rise to chloroplasts. After a century of debate, the presence of chloroplast DNA was demonstrated in the 1960s. The first chloroplast genomes were sequenced in the 1980s, followed by ~100 vegetable, fruit, cereal, beverage, oil and starch/sugar crop chloroplast genomes in the past three decades. Foreign genes were expressed in isolated chloroplasts or intact plant cells in the late 1980s and stably integrated into chloroplast genomes, with typically maternal inheritance shown in the 1990s. Since then, chloroplast genomes conferred the highest reported levels of tolerance or resistance to biotic or abiotic stress. Although launching products with agronomic traits in important crops using this concept has been elusive, commercial products developed include enzymes used in everyday life from processing fruit juice, to enhancing water absorption of cotton fibre or removal of stains as laundry detergents and in dye removal in the textile industry. Plastid genome sequences have revealed the framework of green plant phylogeny as well as the intricate history of plastid genome transfer events to other eukaryotes. Discordant historical signals among plastid genes suggest possible variable constraints across the plastome and further understanding and mitigation of these constraints may yield new opportunities for bioengineering. In this review, we trace the evolutionary history of chloroplasts, status of autonomy and recent advances in products developed for everyday use or those advanced to the clinic, including treatment of COVID-19 patients and SARS-CoV-2 vaccine.

Mangrove tree (Avicennia marina): insight into chloroplast genome evolutionary divergence and its comparison with related species from family Acanthaceae

Avicennia marina (family Acanthaceae) is a halotolerant woody shrub that grows wildly and cultivated in the coastal regions. Despite its importance, the species suffers from lack of genomic datasets to improve its taxonomy and phylogenetic placement across the related species. Here, we have aimed to sequence the plastid genome of A. marina and its comparison with related species in family Acanthaceae. Detailed next-generation sequencing and analysis showed a complete chloroplast genome of 150,279 bp, comprising 38.6% GC. Genome architecture is quadripartite revealing large single copy (82,522 bp), small single copy (17,523 bp), and pair of inverted repeats (25,117 bp). Furthermore, the genome contains 132 different genes, including 87 protein-coding genes, 8 rRNA, 37 tRNA genes, and 126 simple sequence repeats (122 mononucleotide, 2 dinucleotides, and 2 trinucleotides). Interestingly, about 25 forward, 15 reversed and 14 palindromic repeats were also found in the A. marina. High degree synteny was observed in the pairwise alignment with related genomes. The chloroplast genome comparative assessment showed a high degree of sequence similarity in coding regions and varying divergence in the intergenic spacers among ten Acanthaceae species. The pairwise distance showed that A. marina exhibited the highest divergence (0.084) with Justicia flava and showed lowest divergence with Aphelandra knappiae (0.059). Current genomic datasets are a valuable resource for investigating the population and evolutionary genetics of family Acanthaceae members' specifically A. marina and related species.

Evolution and Adaptation of Forest and Crop Pathogens in the Anthropocene

Anthropocene marks the era when human activity is making a significant impact on earth, its ecological and biogeographical systems. The domestication and intensification of agricultural and forest production systems have had a large impact on plant and tree health. Some pathogens benefitted from these human activities and have evolved and adapted in response to the expansion of crop and forest systems, resulting in global outbreaks. Global pathogen genomics data including population genomics and high-quality reference assemblies are crucial for understanding the evolution and adaptation of pathogens. Crops and forest trees have remarkably different characteristics, such as reproductive time and the level of domestication. They also have different production systems for disease management with more intensive management in crops than forest trees. By comparing and contrasting results from pathogen population genomic studies done on widely different agricultural and forest production systems, we can improve our understanding of pathogen evolution and adaptation to different selection pressures. We find that in spite of these differences, similar processes such as hybridization, host jumps, selection, specialization, and clonal expansion are shaping the pathogen populations in both crops and forest trees. We propose some solutions to reduce these impacts and lower the probability of global pathogen outbreaks so that we can envision better management strategies to sustain global food production as well as ecosystem services.

Variation in Plastome Sizes Accompanied by Evolutionary History in Monogenomic Triticeae (Poaceae: Triticeae)

To investigate the pattern of chloroplast genome variation in Triticeae, we comprehensively analyzed the indels in protein-coding genes and intergenic sequence, gene loss/pseudonization, intron variation, expansion/contraction in inverted repeat regions, and the relationship between sequence characteristics and chloroplast genome size in 34 monogenomic Triticeae plants. Ancestral genome reconstruction suggests that major length variations occurred in four-stem branches of monogenomic Triticeae followed by independent changes in each genus. It was shown that the chloroplast genome sizes of monogenomic Triticeae were highly variable. The chloroplast genome of Pseudoroegneria, Dasypyrum, Lophopyrum, Thinopyrum, Eremopyrum, Agropyron, Australopyrum, and Henradia in Triticeae had evolved toward size reduction largely because of pseudogenes elimination events and length deletion fragments in intergenic. The Aegilops/Triticum complex, Taeniatherum, Secale, Crithopsis, Herteranthelium, and Hordeum in Triticeae had a larger chloroplast genome size. The large size variation in major lineages and their subclades are most likely consequences of adaptive processes since these variations were significantly correlated with divergence time and historical climatic changes. We also found that several intergenic regions, such as petN-trnC and psbE-petL containing unique genetic information, which can be used as important tools to identify the maternal relationship among Triticeae species. Our results contribute to the novel knowledge of plastid genome evolution in Triticeae.

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.

SPInDel Analysis of the Non-Coding Regions of cpDNA as a More Useful Tool for the Identification of Rye (Poaceae: Secale) Species

Secale is a small but very diverse genus from the tribe Triticeae (family Poaceae), which includes annual, perennial, self-pollinating and open-pollinating, cultivated, weedy and wild species of various phenotypes. Despite its high economic importance, classification of this genus, comprising 3–8 species, is inconsistent. This has resulted in significantly reduced progress in the breeding of rye which could be enriched with functional traits derived from wild rye species. Our previous research has suggested the utility of non-coding sequences of chloroplast and mitochondrial DNA in studies on closely related species of the genus Secale. Here we applied the SPInDel (Species Identification by Insertions/Deletions) approach, which targets hypervariable genomic regions containing multiple insertions/deletions (indels) and exhibiting extensive length variability. We analysed a total of 140 and 210 non-coding sequences from cpDNA and mtDNA, respectively. The resulting data highlight regions which may represent useful molecular markers with respect to closely related species of the genus Secale, however, we found the chloroplast genome to be more informative. These molecular markers include non-coding regions of chloroplast DNA: atpB-rbcL and trnT-trnL and non-coding regions of mitochondrial DNA: nad1B-nad1C and rrn5/rrn18. Our results demonstrate the utility of the SPInDel concept for the characterisation of Secale species.

Genome-Edited Triple-Recessive Mutation Alters Seed Dormancy in Wheat

Common wheat has three sets of sub-genomes, making mutations difficult to observe, especially for traits controlled by recessive genes. Here, we produced hexaploid wheat lines with loss of function of homeoalleles of Qsd1, which controls seed dormancy in barley, by Agrobacterium-mediated CRISPR/Cas9. Of the eight transformed wheat events produced, three independent events carrying multiple mutations in wheat Qsd1 homeoalleles were obtained. Notably, one line had mutations in every homeoallele. We crossed this plant with wild-type cultivar Fielder to generate a transgene-free triple-recessive mutant, as revealed by Mendelian segregation. The mutant showed a significantly longer seed dormancy period than wild-type, which may result in reduced pre-harvest sprouting of grains on spikes. PCR, southern blotting, and whole-genome shotgun sequencing revealed that this segregant lacked transgenes in its genomic sequence. This technique serves as a model for trait improvement in wheat, particularly for genetically recessive traits, based on locus information from diploid barley.

Wheat Disease Resistance Genes and Their Diversification Through Integrated Domain Fusions

Plants are in a constant evolutionary arms race with their pathogens. At the molecular level, the plant nucleotide-binding leucine-rich repeat receptors (NLRs) family has coevolved with rapidly evolving pathogen effectors. While many NLRs utilize variable leucine-rich repeats (LRRs) to detect effectors, some have gained integrated domains (IDs) that may be involved in receptor activation or downstream signaling. The major objectives of this project were to identify NLR genes in wheat (Triticum aestivum L.) and assess IDs associated with immune signaling (e.g., kinase and transcription factor domains). We identified 2,151 NLR-like genes in wheat, of which 1,298 formed 547 gene clusters. Among the non-toll/interleukin-1 receptor NLR (non-TNL)-like genes, 1,552 encode LRRs, 802 are coiled-coil (CC) domain-encoding (CC-NBS-LRR or CNL) genes, and three encode resistance to powdery mildew 8 (RPW8) domains (RPW8-NBS-LRR or RNL). The expansion of the NLR gene family in wheat is attributable to its origin by recent polyploidy events. Gene clusters were likely formed by tandem duplications, and wheat NLR phylogenetic relationships were similar to those in barley and Aegilops. We also identified wheat NLR-ID fusion proteins as candidates for NLR functional diversification, often as kinase and transcription factor domains. Comparative analyses of the IDs revealed evolutionary conservation of more than 80% amino acid sequence similarity. Homology assessment indicates that these domains originated as functional non-NLR-encoding genes that were incorporated into NLR-encoding genes through duplication events. We also found that many of the NLR-ID genes encode alternative transcripts that include or exclude IDs, a phenomenon that seems to be conserved among species. To verify this, we have analyzed the alternative transcripts that include or exclude an ID of an NLR-ID from another monocotyledon species, rice (Oryza sativa). This indicates that plants employ alternative splicing to regulate IDs, possibly using them as baits, decoys, and functional signaling components. Genomic and expression data support the hypothesis that wheat uses alternative splicing to include and exclude IDs from NLR proteins.

History and future perspectives of barley genomics

Barley (Hordeum vulgare), one of the most widely cultivated cereal crops, possesses a large genome of 5.1 Gbp. Through various international collaborations, the genome has recently been sequenced and assembled at the chromosome-scale by exploiting available genetic and genomic resources. Many wild and cultivated barley accessions have been collected and preserved around the world. These accessions are crucial to obtain diverse natural and induced barley variants. The barley bioresource project aims to investigate the diversity of this crop based on purified seed and DNA samples of a large number of collected accessions. The long-term goal of this project is to analyse the genome sequences of major barley accessions worldwide. In view of technical limitations, a strategy has been employed to establish the exome structure of a selected number of accessions and to perform high-quality chromosome-scale assembly of the genomes of several major representative accessions. For the future project, an efficient annotation pipeline is essential for establishing the function of genomes and genes as well as for using this information for sequence-based digital barley breeding. In this article, the author reviews the existing barley resources along with their applications and discuss possible future directions of research in barley genomics.

Illegitimate Recombination Between Homeologous Genes in Wheat Genome

Polyploidies produce a large number of duplicated regions and genes in genomes, which have a long-term impact and stimulate genetic innovation. The high similarity between homeologous chromosomes, forming different subgenomes, or homologous regions after genome repatterning, may permit illegitimate DNA recombination. Here, based on gene colinearity, we aligned the (sub)genomes of common wheat (Triticum aestivum, AABBDD genotype) and its relatives, including Triticum urartu (AA), Aegilops tauschii (DD), and T. turgidum ssp. dicoccoides (AABB) to detect the homeologous (paralogous or orthologous) colinear genes within and between (sub)genomes. Besides, we inferred more ancient paralogous regions produced by a much ancient grass-common tetraploidization. By comparing the sequence similarity between paralogous and orthologous genes, we assumed abnormality in the topology of constructed gene trees, which could be explained by gene conversion as a result of illegitimate recombination. We found large numbers of inferred converted genes (>2,000 gene pairs) suggested long-lasting genome instability of the hexaploid plant, and preferential donor roles by DD genes. Though illegitimate recombination was much restricted, duplicated genes produced by an ancient whole-genome duplication, which occurred millions of years ago, also showed evidence of likely gene conversion. As to biological function, we found that ~40% catalytic genes in colinearity, including those involved in starch biosynthesis, were likely affected by gene conversion. The present study will contribute to understanding the functional and structural innovation of the common wheat genome.

Global Role of Crop Genomics in the Face of Climate Change

The development of climate change resilient crops is necessary if we are to meet the challenge of feeding the growing world's population. We must be able to increase food production despite the projected decrease in arable land and unpredictable environmental conditions. This review summarizes the technological and conceptual advances that have the potential to transform plant breeding, help overcome the challenges of climate change, and initiate the next plant breeding revolution. Recent developments in genomics in combination with high-throughput and precision phenotyping facilitate the identification of genes controlling critical agronomic traits. The discovery of these genes can now be paired with genome editing techniques to rapidly develop climate change resilient crops, including plants with better biotic and abiotic stress tolerance and enhanced nutritional value. Utilizing the genetic potential of crop wild relatives (CWRs) enables the domestication of new species and the generation of synthetic polyploids. The high-quality crop plant genome assemblies and annotations provide new, exciting research targets, including long non-coding RNAs (lncRNAs) and cis-regulatory regions. Metagenomic studies give insights into plant-microbiome interactions and guide selection of optimal soils for plant cultivation. Together, all these advances will allow breeders to produce improved, resilient crops in relatively short timeframes meeting the demands of the growing population and changing climate.