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

Homoeologous exchanges occur through intragenic recombination generating novel transcripts and proteins in wheat and other polyploids

Recombination between homeologous chromosomes, also known as homeologous exchange (HE), plays a significant role in shaping genome structure and gene expression in interspecific hybrids and allopolyploids of several plant species. However, the molecular mechanisms that govern HEs are not well understood. Here, we studied HE events in the progeny of a nascent allotetraploid (genome AADD) derived from two diploid progenitors of hexaploid bread wheat using cytological and whole-genome sequence analyses. In total, 37 HEs were identified and HE junctions were mapped precisely. HEs exhibit typical patterns of homologous recombination hotspots, being biased toward low-copy, subtelomeric regions of chromosome arms and showing association with known recombination hotspot motifs. But, strikingly, while homologous recombination preferentially takes place upstream and downstream of coding regions, HEs are highly enriched within gene bodies, giving rise to novel recombinant transcripts, which in turn are predicted to generate new protein fusion variants. To test whether this is a widespread phenomenon, a dataset of high-resolution HE junctions was analyzed for allopolyploid Brassica, rice, Arabidopsis suecica, banana, and peanut. Intragenic recombination and formation of chimeric genes was detected in HEs of all species and was prominent in most of them. HE thus provides a mechanism for evolutionary novelty in transcript and protein sequences in nascent allopolyploids.

Drought Stress Tolerance and Photosynthetic Activity of Alloplasmic Lines T. dicoccum x T. aestivum

Tetraploid species T. dicoccum Shuebl is a potential source of drought tolerance for cultivated wheat, including common wheat. This paper describes the genotyping of nine stable allolines isolated in the offspring from crossing of T. dicoccum x T. aestivum L. using 21 microsatellite (simple sequence repeats-SSR) markers and two cytoplasmic mitochondrial markers to orf256, rps19-p genes; evaluation of drought tolerance of allolines at different stages of ontogenesis (growth parameters, relative water content, quantum efficiency of Photosystem II, electron transport rate, energy dissipated in Photosystem II); and the study of drought tolerance regulator gene Dreb-1 with allele-specific PCR (AS-MARKER) and partial sequence analysis. Most allolines differ in genomic composition and T. dicoccum introgressions. Four allolines-D-b-05, D-d-05, D-d-05b, and D-41-05-revealed signs of drought tolerance of varying degrees. The more drought tolerant D-41-05 line was also characterized by Dreb-B1 allele introgression from T. dicoccum. A number of non-specific patterns and significant differences in allolines in regulation of physiological parameters in drought conditions is identified. Changes in photosynthetic activity in stress-drought are shown to reflect the level of drought tolerance of the forms studied. The contribution of different combinations of nuclear/cytoplasmic genome and alleles of Dreb-1 gene in allolines to the formation of stress tolerance and photosynthetic activity is discussed.

Evolutionary patterns of plastome uncover diploid-polyploid maternal relationships in Triticeae

To investigate the diploid-polyploid relationships and the role of maternal progenitors in establishment of polyploid richness in Triticeae, 35 polyploids representing almost all genomic constitutions together with 48 diploid taxa representing 20 basic genomes in the tribe were analyzed. Phylogenomic reconstruction, genetic distance matrix, and nucleotide diversity patterns of plastome sequences indicated that (1) The maternal donor of the annual polyploid species with the U- and D-genome are related to extant Ae. umbellulata and Ae. tauschii, respectively. The maternal donor to the annual polyploid species with the S-, G-, and B-genome originated from the species of Sitopsis section of the genus Aegilops. The annual species with the Xe-containing polyploids were donated by Eremopyrum as the female parent; (2) Pseudoroegneria and Psathyrostachys were the maternal donor of perennial species with the St- and Ns-containing polyploids, respectively; (3) The Lophopyrum, Thinopyrum and Dasypyrum genomes contributed cytoplasm genome to Pseudoroegneria species as a result of incomplete lineage sorting and/or chloroplast captures, and these lineages were genetically transmitted to the St-containing polyploid species via polyploidization; (4) There is a reticulate relationship among the St-containing polyploid species. It can be suggested that genetic heterogeneity might associate with the richness of the polyploids in Triticeae.

Complete Chloroplast Genome Sequence of Hibiscus cannabinus and Comparative Analysis of the Malvaceae Family

Kenaf (Hibiscus cannabinus) is one of the most fast-growing bast in the world and belongs to the family Malvaceae. However, the systematic classification and chloroplast (cp) genome of kenaf has not been reported to date. In this study, we sequenced the cp genome of kenaf and conducted phylogenetic and comparative analyses in the family of Malvaceae. The sizes of H. cannabinus cp genomes were 162,903 bp in length, containing 113 unique genes (79 protein-coding genes, four rRNA genes, and 30 tRNA genes). Phylogenetic analysis indicated that the cp genome sequence of H. cannabinus has closer relationships with Talipariti hamabo and Abelmoschus esculentus than with Hibiscus syriacus, which disagrees with the taxonomical relationship. Further analysis obtained a new version of the cp genome annotation of H. syriacus and found that the orientation variation of small single copy (SSC) region exists widely in the family of Malvaceae. The highly variable ycf1 and the highly conserved gene rrn32 were identified among the family of Malvaceae. In particular, the explanation for two different SSC orientations in the cp genomes associated with phylogenetic analysis is discussed. These results provide insights into the systematic classification of the Hibiscus genus in the Malvaceae family.

The complete chloroplast genomes of seventeen Aegilops tauschii: genome comparative analysis and phylogenetic inference

The D genome progenitor of bread wheat, Aegilops tauschii Cosson (DD, 2n = 2x = 14), which is naturally distributed in Central Eurasia, ranging from northern Syria and Turkey to western China, is considered a potential genetic resource for improving bread wheat. In this study, the chloroplast (cp) genomes of 17 Ae. tauschii accessions were reconstructed. The cp genome sizes ranged from 135,551 bp to 136,009 bp and contained a typical quadripartite structure of angiosperms. Within these genomes, we identified a total of 124 functional genes, including 82 protein-coding genes, 34 transfer RNA genes and eight ribosomal RNA genes, with 17 duplicated genes in the IRs. Although the comparative analysis revealed that the genomic structure (gene order, gene number and IR/SC boundary regions) is conserved, a few variant loci were detected, predominantly in the non-coding regions (intergenic spacer regions). The phylogenetic relationships determined based on the complete genome sequences were consistent with the hypothesis that Ae. tauschii populations in the Yellow River region of China originated in South Asia not Xinjiang province or Iran, which could contribute to more effective utilization of wild germplasm resources. Furthermore, we confirmed that Ae. tauschii was derived from monophyletic speciation rather than hybrid speciation at the cp genome level. We also identified four variable genomic regions, rpl32-trnL-UAG, ccsA-ndhD, rbcL-psaI and rps18-rpl20, showing high levels of nucleotide polymorphisms, which may accordingly prove useful as cpDNA markers in studying the intraspecific genetic structure and diversity of Ae. tauschii.

Nuclear Integrants of Organellar DNA Contribute to Genome Structure and Evolution in Plants

The transfer of genetic material from the mitochondria and plastid to the nucleus gives rise to nuclear integrants of mitochondrial DNA (NUMTs) and nuclear integrants of plastid DNA (NUPTs). This frequently occurring DNA transfer is ongoing and has important evolutionary implications. In this review, based on previous studies and the analysis of NUMT/NUPT insertions of more than 200 sequenced plant genomes, we analyzed and summarized the general features of NUMTs/NUPTs and highlighted the genetic consequence of organellar DNA insertions. The statistics of organellar DNA integrants among various plant genomes revealed that organellar DNA-derived sequence content is positively correlated with the nuclear genome size. After integration, the nuclear organellar DNA could undergo different fates, including elimination, mutation, rearrangement, fragmentation, and proliferation. The integrated organellar DNAs play important roles in increasing genetic diversity, promoting gene and genome evolution, and are involved in sex chromosome evolution in dioecious plants. The integrating mechanisms, involving non-homologous end joining at double-strand breaks were also discussed.

Thermal plasticity of the circadian clock is under nuclear and cytoplasmic control in wild barley

Temperature compensation, expressed as the ability to maintain clock characteristics (mainly period) in face of temperature changes, that is, robustness, is considered a key feature of circadian clock systems. In this study, we explore the genetic basis for lack of robustness, that is, plasticity, of circadian clock as reflected by photosynthesis rhythmicity. The clock rhythmicity of a new wild barley reciprocal doubled haploid population was analysed with a high temporal resolution of pulsed amplitude modulation of chlorophyll fluorescence under optimal (22°C) and high (32°C) temperature. This comparison between two environments pointed to the prevalence of clock acceleration under heat. Genotyping by sequencing of doubled haploid lines indicated a rich recombination landscape with minor fixation (less than 8%) for one of the parental alleles. Quantitative genetic analysis included genotype by environment interactions and binary-threshold models. Variation in the circadian rhythm plasticity phenotypes, expressed as change (delta) of period and amplitude under two temperatures, was associated with maternal organelle genome (the plasmotype), as well as with several nuclear loci. This first reported rhythmicity driven by nuclear loci and plasmotype with few identified variants, paves the way for studying impact of cytonuclear variations on clock robustness and on plant adaptation to changing environments.

Gene Expression Profiling and Fine Mapping Identifies a Gibberellin 2-Oxidase Gene Co-segregating With the Dominant Dwarfing Gene Ddw1 in Rye (Secale cereale L.)

The gibberellin (GA)-sensitive dwarfing gene Ddw1 provides an opportunity to genetically reduce plant height in rye. Genetic analysis in a population of recombinant inbred lines confirmed a monogenetic dominant inheritance of Ddw1. Significant phenotypic differences in PH between homo- and heterozygotic genotypes indicate an incomplete dominance of Ddw1. De novo transcriptome sequencing of Ddw1 mutant as well as tall genotypes resulted in 113,547 contigs with an average length of 318 bp covering 36.18 Mbp rye DNA. A hierarchical cluster analysis based on individual groups of rye homologs of functionally characterized rice genes controlling morphological or physiological traits including plant height, flowering time, and source activity identified the gene expression profile of stems at the begin of heading to most comprehensively mirror effects of Ddw1. Genome-wide expression profiling identified 186 transcripts differentially expressed between semi-dwarf and tall genotypes in stems. In total, 29 novel markers have been established and mapped to a 27.2 cM segment in the distal part of the long arm of chromosome 5R. Ddw1 could be mapped within a 0.4 cM interval co-segregating with a marker representing the C20-GA2-oxidase gene ScGA2ox12, that is up-regulated in stems of Ddw1 genotypes. The increased expression of ScGA2ox12 observed in semi-dwarf rye as well as structural alterations in transcript sequences associated with the ScGA2ox12 gene implicate, that Ddw1 is a dominant gain-of-function mutant. Integration of the target interval in the wheat reference genome sequence indicated perfect micro-colinearity between the Ddw1 locus and a 831 kb segment on chromosome 5A, which resides inside of a 11.21 Mb interval carrying the GA-sensitive dwarfing gene Rht12 in wheat. The potential of Ddw1 as a breeder's option to improve lodging tolerance in rye is discussed.

Lineage-Specific Evolved MicroRNAs Regulating NB-LRR Defense Genes in Triticeae

Disease resistance genes encoding proteins with nucleotide binding sites and Leucine-Rich Repeat (NB-LRR) domains include many members involved in the effector-triggered immunity pathway in plants. The transcript levels of these defense genes are negatively regulated by diverse microRNAs (miRNAs) in angiosperms and gymnosperms. In wheat, using small RNA expression datasets and degradome datasets, we identified five miRNA families targeting NB-LRR defense genes in monocots, some of which arose in the Triticeae species era. These miRNAs regulate different types of NB-LRR genes, most of them with coil-coiled domains, and trigger the generation of secondary small interfering RNAs (siRNA) as a phased pattern in the target site regions. In addition to acting in response to biotic stresses, they are also responsive to abiotic stresses such as heat, drought, salt, and light stress. Their copy number and expression variation in Triticeae suggest a rapid birth and death frequency. Altogether, non-conserved miRNAs as conserved transcriptional regulators in gymnosperms and angiosperms regulating the disease resistance genes displayed quick plasticity including the variations of sequences, gene copy number, functions, and expression level, which accompanied with NB-LRR genes may be tune-regulated to plants in natural environments with various biotic and abiotic stresses.

Cytonuclear Coevolution following Homoploid Hybrid Speciation in Aegilops tauschii

The diploid D-genome lineage of the Triticum/Aegilops complex has an evolutionary history involving genomic contributions from ancient A- and B/S-genome species. We explored here the possible cytonuclear evolutionary responses to this history of hybridization. Phylogenetic analysis of chloroplast DNAs indicates that the D-genome lineage has a maternal origin of the A-genome or some other closely allied lineage. Analyses of the nuclear genome in the D-genome species Aegilops tauschii indicate that accompanying and/or following this ancient hybridization, there has been biased maintenance of maternal A-genome ancestry in nuclear genes encoding cytonuclear enzyme complexes (CECs). Our study provides insights into mechanisms of cytonuclear coevolution accompanying the evolution and eventual stabilization of homoploid hybrid species. We suggest that this coevolutionary process includes likely rapid fixation of A-genome CEC orthologs as well as biased retention of A-genome nucleotides in CEC homologs following population level recombination during the initial generations.

Btr1-A Induces Grain Shattering and Affects Spike Morphology and Yield-Related Traits in Wheat

Spike brittleness represents an important domestication trait in crops. Although the brittle rachis of wild wheat was cloned, however, the molecular mechanism underlying spike brittleness is yet to be elucidated. Here, we identified a single dominant brittle rachis gene Br-Ab on chromosome arm 3AbS using an F2 population of diploid wheat and designated Btr1-Ab. Sequence analysis of the Btr1-A gene in 40 diploid wheat accessions, 80 tetraploid wheat accessions and 38 hexaploid wheat accessions showed that two independent mutations (Ala119Thr for diploid and Gly97* for polyploids) in the Btr1-A coding region resulting in the nonbrittle rachis allele. Overexpression of Btr1-Ab in nonbrittle hexaploid wheat led to brittle rachis in transgenic plants. RNA-Seq analysis revealed that Btr1-A represses the expression of cell wall biosynthesis genes during wheat rachis development. In addition, we found that Btr1-A can modify spike morphology and reduce threshability, grain size and thousand grain weight in transgenic wheat. These results demonstrated that Btr1-A reduces cell wall synthesis in rachis nodes, resulting in natural spikelet shattering, and that the transition from Btr1-A to btr1-A during wheat domestication had profound effects on evolution of spike morphology and yield-related traits.

Broadening the bread wheat D genome

Although Ae. tauschii has been extensively utilised for wheat breeding, the D-genome-containing allopolyploids have largely remained unexploited. In this review, we discuss approaches that can be used to exploit the D genomes of the different Aegilops species for the improvement of bread wheat. The D genome of allohexaploid bread wheat (Triticum aestivum, 2n = AABBDD) is the least diverse of the three wheat genomes and is unarguably less diverse than that of diploid progenitor Aegilops tauschii (2n = DD). Useful genetic variation and phenotypic traits also exist within each of the wheat group species containing a copy of the D genome: allopolyploid Aegilops species Ae. cylindrica (2n = D^cD^cC^cC^c), Ae. crassa 4x (2n = D¹D¹X^crX^cr), Ae. crassa 6x (2n = D¹D¹X^crX^crD^crD^cr), Ae. ventricosa (2n = D^vD^vN^vN^v), Ae. vavilovii (2n = D¹D¹X^crX^crS^vS^v) and Ae. juvenalis (2n = D¹D¹X^crX^crU^jU^j). Although Ae. tauschii has been extensively utilised for wheat breeding, the D-genome-containing allopolyploids have largely remained unexploited. Some of these D genomes appear to be modified relative to the bread wheat and Ae. tauschii D genomes, and others present in the allopolyploids may also contain useful variation as a result of adaptation to an allopolyploid, multi-genome environment. We summarise the genetic relationships, karyotypic variation and phenotypic traits known to be present in each of the D genome species that could be of relevance for bread wheat improvement and discuss approaches that can be used to exploit the D genomes of the different Aegilops species for the improvement of bread wheat. Better understanding of factors controlling chromosome inheritance and recombination in wheat group interspecific hybrids, as well as effective utilisation of new and developing genetics and genomics technologies, have great potential to improve the agronomic potential of the bread wheat D genome.

Genome constitution and evolution of Elytrigia lolioides inferred from Acc1, EF-G, ITS, TrnL-F sequences and GISH

BACKGROUND

Elytrigia lolioides (Kar. et Kir.) Nevski, which is a perennial, cross-pollinating wheatgrass that is distributed in Russia and Kazakhstan, is classified into Elytrigia, Elymus, and Lophopyrum genera by taxonomists on the basis of different taxonomic classification systems. However, the genomic constitution of E. lolioides is still unknown. To identify the genome constitution and evolution of E. lolioides, we used single-copy nuclear genes acetyl-CoA carboxylase (Acc1) and elongation factor G (EF-G), multi-copy nuclear gene internal transcribed space (ITS), chloroplast gene trnL-F together with fluorescence and genomic in situ hybridization.

RESULTS

Despite the widespread homogenization of ITS sequences, two distinct lineages (genera Pseudoroegneria and Hordeum) were identified. Acc1 and EF-G sequences suggested that in addition to Pseudoroegneria and Hordeum, unknown genome was the third potential donor of E. lolioides. Data from chloroplast DNA showed that Pseudoroegneria is the maternal donor of E. lolioides. Data from specific FISH marker for St genome indicated that E. lolioides has two sets of St genomes. Both genomic in situ hybridization (GISH) and fluorescence in situ hybridization (FISH) results confirmed the presence of Hordeum genome in this species. When E genome was used as the probe, no signal was found in 42 chromosomes. The E-like copy of Acc1 sequences was detected in E. lolioides possibly due to the introgression from E genome species. One of the H chromosomes in the accession W6-26586 from Kazakhstan did not hybridize H genome signals but had St genome signals on the pericentromeric regions in the two-color GISH.

CONCLUSIONS

Phylogenetic and in situ hybridization indicated the presence of two sets of Pseudoroegneria and one set of Hordeum genome in E. lolioides. The genome formula of E. lolioides was designed as StStStStHH. E. lolioides may have originated through the hybridization between tetraploid Elymus (StH) and diploid Pseudoroegneria species. E and unknown genomes may participate in the speciation of E. lolioides through introgression. According to the genome classification system, E. lolioides should be transferred into Elymus L. and renamed as Elymus lolioidus (Kar. er Kir.) Meld.

RNAi-suppression of barley caffeic acid O-methyltransferase modifies lignin despite redundancy in the gene family

Caffeic acid O-methyltransferase (COMT), the lignin biosynthesis gene modified in many brown-midrib high-digestibility mutants of maize and sorghum, was targeted for downregulation in the small grain temperate cereal, barley (Hordeum vulgare), to improve straw properties. Phylogenetic and expression analyses identified the barley COMT orthologue(s) expressed in stems, defining a larger gene family than in brachypodium or rice with three COMT genes expressed in lignifying tissues. RNAi significantly reduced stem COMT protein and enzyme activity, and modestly reduced stem lignin content while dramatically changing lignin structure. Lignin syringyl-to-guaiacyl ratio was reduced by ~50%, the 5-hydroxyguaiacyl (5-OH-G) unit incorporated into lignin at 10--15-fold higher levels than normal, and the amount of p-coumaric acid ester-linked to cell walls was reduced by ~50%. No brown-midrib phenotype was observed in any RNAi line despite significant COMT suppression and altered lignin. The novel COMT gene family structure in barley highlights the dynamic nature of grass genomes. Redundancy in barley COMTs may explain the absence of brown-midrib mutants in barley and wheat. The barley COMT RNAi lines nevertheless have the potential to be exploited for bioenergy applications and as animal feed.

Genetic diversity and relationship between cultivated, weedy and wild rye species as revealed by chloroplast and mitochondrial DNA non-coding regions analysis

The genus Secale is small but very diverse. Despite the high economic importance, phylogenetic relationships of rye species have not been fully determined, and they are extremely important for the process of breeding of new cultivars that can be enriched with functional traits derived from wild rye species. The study analyzed the degree of relationship of 35 accessions of the genus Secale, representing 13 most often distinguished species and subspecies, originating from various seed collections in the world, based on the analysis of non-coding regions of the chloroplast (cpDNA) and mitochondrial genome (mtDNA), widely used in phylogenetic and population plant studies, because of a higher rate of evolution than the coding regions. There was no clear genetic structure between different species and subspecies, which may indicated the introgression between these taxa. The obtained data confirmed that S. vavilovii was very similar to S. cereale, which confirmed the assumption that they might share a common ancestor. The results also confirmed the divergence of S. sylvestre from other species and subspecies of rye. Areas that may be useful molecular markers in studies on closely related species of the genus Secale were also indicated.

Complete Chloroplast Genome Sequences of Kaempferia Galanga and Kaempferia Elegans: Molecular Structures and Comparative Analysis

Kaempferia galanga and Kaempferia elegans, which belong to the genus Kaempferia family Zingiberaceae, are used as valuable herbal medicine and ornamental plants, respectively. The chloroplast genomes have been used for molecular markers, species identification and phylogenetic studies. In this study, the complete chloroplast genome sequences of K. galanga and K. elegans are reported. Results show that the complete chloroplast genome of K. galanga is 163,811 bp long, having a quadripartite structure with large single copy (LSC) of 88,405 bp and a small single copy (SSC) of 15,812 bp separated by inverted repeats (IRs) of 29,797 bp. Similarly, the complete chloroplast genome of K. elegans is 163,555 bp long, having a quadripartite structure in which IRs of 29,773 bp length separates 88,020 bp of LSC and 15,989 bp of SSC. A total of 111 genes in K. galanga and 113 genes in K. elegans comprised 79 protein-coding genes and 4 ribosomal RNA (rRNA) genes, as well as 28 and 30 transfer RNA (tRNA) genes in K. galanga and K. elegans, respectively. The gene order, GC content and orientation of the two Kaempferia chloroplast genomes exhibited high similarity. The location and distribution of simple sequence repeats (SSRs) and long repeat sequences were determined. Eight highly variable regions between the two Kaempferia species were identified and 643 mutation events, including 536 single-nucleotide polymorphisms (SNPs) and 107 insertion/deletions (indels), were accurately located. Sequence divergences of the whole chloroplast genomes were calculated among related Zingiberaceae species. The phylogenetic analysis based on SNPs among eleven species strongly supported that K. galanga and K. elegans formed a cluster within Zingiberaceae. This study identified the unique characteristics of the entire K. galanga and K. elegans chloroplast genomes that contribute to our understanding of the chloroplast DNA evolution within Zingiberaceae species. It provides valuable information for phylogenetic analysis and species identification within genus Kaempferia.

Green Biotechnology: A Brief Update on Plastid Genome Engineering

Plant genetic engineering has become an inevitable tool in the molecular breeding of crops. Significant progress has been made in the generation of novel plastid transformation vectors and optimized transformation protocols. There are several advantages of plastid genome engineering over conventional nuclear transformation. Some of the advantages include multigene engineering by expression of biosynthetic pathway genes as operons, extremely high-level expression of protein accumulation, lack of transgene silencing, etc. Transgene containment owing to maternal inheritance is another important advantage of plastid genome engineering. Chloroplast genome modification usually results in alteration of several thousand plastid genome copies in a cell. Several therapeutic proteins, edible vaccines, antimicrobial peptides, and industrially important enzymes have been successfully expressed in chloroplasts so far. Here, we critically recapitulate the latest developments in plastid genome engineering. Latest advancements in plastid genome sequencing are briefed. In addition, advancement of extending the toolbox for plastid engineering for selected applications in the area of molecular farming and production of industrially important enzyme is briefed.

Silencing of a Unique Integrated Domain Nucleotide-Binding Leucine-Rich Repeat Gene in Wheat Abolishes Diuraphis noxia Resistance

Plants respond in a similar manner to aphid feeding as to pathogen attack. Diuraphis noxia is a specialist aphid, feeding only on selected grasses that include wheat, barley, and oats. The wheat-Diuraphis noxia interaction is characterized by responses very similar to those seen in wheat-pathogen interactions with none of the underlying resistance pathways and genes characterized yet. From wheat harboring the Dn1 resistance gene, we have identified a nucleotide-binding leucine-rich repeat (NLR) gene containing two integrated domains (IDs). These are three C-terminus ankyrin repeat domains and an N-terminus WRKY domain. The NLR core of the gene can be traced through speciation events within the grass family, with a recent WRKY domain integration that is Triticum-specific. Virus-induced gene silencing of the gene in a resistant wheat line resulted in the abolishment of the resistance response and induced a highly susceptible phenotype. Silenced plants supported a higher number of aphids, similar to the susceptible near-isogenic line (NIL), and the intrinsic rate of increase of the aphids matched that of aphids feeding on the susceptible NIL. The presence of the gene is necessary for Dn1 resistance and we have named the gene Associated with Dn resistance 1 (Adnr1) to reflect this function.

Impact of transposable elements on genome structure and evolution in bread wheat

BACKGROUND

Transposable elements (TEs) are major components of large plant genomes and main drivers of genome evolution. The most recent assembly of hexaploid bread wheat recovered the highly repetitive TE space in an almost complete chromosomal context and enabled a detailed view into the dynamics of TEs in the A, B, and D subgenomes.

RESULTS

The overall TE content is very similar between the A, B, and D subgenomes, although we find no evidence for bursts of TE amplification after the polyploidization events. Despite the near-complete turnover of TEs since the subgenome lineages diverged from a common ancestor, 76% of TE families are still present in similar proportions in each subgenome. Moreover, spacing between syntenic genes is also conserved, even though syntenic TEs have been replaced by new insertions over time, suggesting that distances between genes, but not sequences, are under evolutionary constraints. The TE composition of the immediate gene vicinity differs from the core intergenic regions. We find the same TE families to be enriched or depleted near genes in all three subgenomes. Evaluations at the subfamily level of timed long terminal repeat-retrotransposon insertions highlight the independent evolution of the diploid A, B, and D lineages before polyploidization and cases of concerted proliferation in the AB tetraploid.

CONCLUSIONS

Even though the intergenic space is changed by the TE turnover, an unexpected preservation is observed between the A, B, and D subgenomes for features like TE family proportions, gene spacing, and TE enrichment near genes.

Structural variation and rates of genome evolution in the grass family seen through comparison of sequences of genomes greatly differing in size

Homology was searched with genes annotated in the Aegilops tauschii pseudomolecules against genes annotated in the pseudomolecules of tetraploid wild emmer wheat, Brachypodium distachyon, sorghum and rice. Similar searches were performed with genes annotated in the rice pseudomolecules. Matrices of collinear genes and rearrangements in their order were constructed. Optical BioNano genome maps were constructed and used to validate rearrangements unique to the wild emmer and Ae. tauschii genomes. Most common rearrangements were short paracentric inversions and short intrachromosomal translocations. Intrachromosomal translocations outnumbered segmental intrachromosomal duplications. The densities of paracentric inversion lengths were approximated by exponential distributions in all six genomes. Densities of collinear genes along the Ae. tauschii chromosomes were highly correlated with meiotic recombination rates but those of rearrangements were not, suggesting different causes of the erosion of gene collinearity and evolution of major chromosome rearrangements. Frequent rearrangements sharing breakpoints suggested that chromosomes have been rearranged recurrently at some sites. The distal 4 Mb of the short arms of rice chromosomes Os11 and Os12 and corresponding regions in the sorghum, B. distachyon and Triticeae genomes contain clusters of interstitial translocations including from 1 to 7 collinear genes. The rates of acquisition of major rearrangements were greater in the large wild emmer wheat and Ae. tauschii genomes than in the lineage preceding their divergence or in the B. distachyon, rice and sorghum lineages. It is suggested that synergy between large quantities of dynamic transposable elements and annual growth habit have been the primary causes of the fast evolution of the Triticeae genomes.

Two Rye Genes Responsible for Abnormal Development of Wheat⁻Rye Hybrids Are Linked in the Vicinity of an Evolutionary Translocation on Chromosome 6R

The post-zygotic reproductive isolation (RI) in plants is frequently based on the negative interaction of the parental genes involved in plant development. Of special interest is the study of such types of interactions in crop plants, because of the importance of distant hybridization in plant breeding. This study is devoted to map rye genes that are incompatible with wheat, determining the development of the shoot apical meristem in wheat⁻rye hybrids. Linkage analysis of microsatellite loci, as well as genes of embryo lethality (Eml-R1) and hybrid dwarfness (Hdw-R1) was carried out in hybrids of Chinese Spring wheat with recombinant inbred lines as well as interline rye hybrids. Eml-R1 and Hdw-R1 could be mapped proximal and distal of two closely linked EST-SSR markers, Xgrm902 and Xgrm959, on rye chromosome 6R. Both rye genes are located on a segment of chromosome 6R that contains a breakpoint of evolutionary translocation between the ancestral chromosomes of homeologous groups 6 and 3. The obtained results are discussed in relation to genes interacting in developmental pathways as a class of causal genes of RI.

Next generation sequencing reveals packaging of host RNAs by brome mosaic virus

Although RNA viruses evolved the mechanisms of specific encapsidation, miss-packaging of cellular RNAs has been reported in such RNA virus systems as flock house virus or cucumber necrosis virus. To find out if brome mosaic virus (BMV), a tripartite RNA virus, can package cellular RNAs, BMV was propagated in barley and in Nicotiana benthamiana hosts, purified by cesium chloride (CsCl) gradient ultracentrifugation followed by nuclease treatment to remove any contaminating cellular (host) RNAs. The extracted virion RNA was then sequenced by using next-generation sequencing (NGS RNA-Seq) with the Illumina protocol. Bioinformatic analysis revealed the content of host RNAs ranging from 0.07% for BMV extracted from barley to 0.10% for the virus extracted from N. benthamiana. The viruses from two sources appeared to co-encapsidate different patterns of host-RNAs, including ribosomal RNAs (rRNAs), messenger RNAs (mRNAs) but also mitochondrial and plastid RNAs and, interestingly, transposable elements, both transposons and retrotransposons. Our data reveal that BMV virions can carry host RNAs, having a potential to mediate horizontal gene transfer (HGT) in plants.

Regulation and Evolution of NLR Genes: A Close Interconnection for Plant Immunity

NLR (NOD-like receptor) genes belong to one of the largest gene families in plants. Their role in plants' resistance to pathogens has been clearly described for many members of this gene family, and dysregulation or overexpression of some of these genes has been shown to induce an autoimmunity state that strongly affects plant growth and yield. For this reason, these genes have to be tightly regulated in their expression and activity, and several regulatory mechanisms are described here that tune their gene expression and protein levels. This gene family is subjected to rapid evolution, and to maintain diversity at NLRs, a plethora of genetic mechanisms have been identified as sources of variation. Interestingly, regulation of gene expression and evolution of this gene family are two strictly interconnected aspects. Indeed, some examples have been reported in which mechanisms of gene expression regulation have roles in promotion of the evolution of this gene family. Moreover, co-evolution of the NLR gene family and other gene families devoted to their control has been recently demonstrated, as in the case of miRNAs.

Comparative plastome genomics and phylogenomics of Brachypodium: flowering time signatures, introgression and recombination in recently diverged ecotypes

Few pan-genomic studies have been conducted in plants, and none of them have focused on the intraspecific diversity and evolution of their plastid genomes. We address this issue in Brachypodium distachyon and its close relatives B. stacei and B. hybridum, for which a large genomic data set has been compiled. We analyze inter- and intraspecific plastid comparative genomics and phylogenomic relationships within a family-wide framework. Major indel differences were detected between Brachypodium plastomes. Within B. distachyon, we detected two main lineages, a mostly Extremely Delayed Flowering (EDF+) clade and a mostly Spanish (S+) - Turkish (T+) clade, plus nine chloroplast capture and two plastid DNA (ptDNA) introgression and micro-recombination events. Early Oligocene (30.9 million yr ago (Ma)) and Late Miocene (10.1 Ma) divergence times were inferred for the respective stem and crown nodes of Brachypodium and a very recent Mid-Pleistocene (0.9 Ma) time for the B. distachyon split. Flowering time variation is a main factor driving rapid intraspecific divergence in B. distachyon, although it is counterbalanced by repeated introgression between previously isolated lineages. Swapping of plastomes between the three different genomic groups, EDF+, T+, S+, probably resulted from random backcrossing followed by stabilization through selection pressure.

Reconstructing the Evolutionary History of Powdery Mildew Lineages (Blumeria graminis) at Different Evolutionary Time Scales with NGS Data

Blumeria graminis (Ascomycota) includes fungal pathogens that infect numerous grasses and cereals. Despite its economic impact on agriculture and its scientific importance in plant–pathogen interaction studies, the evolution of different lineages with different host ranges is poorly understood. Moreover, the taxonomy of grass powdery mildew is rather exceptional: there is only one described species (B. graminis) subdivided in different formae speciales (ff.spp.), which are defined by their host range. In this study we applied phylogenomic and population genomic methods to whole genome sequence data of 31 isolates of B. graminis belonging to different ff.spp. and reconstructed the evolutionary relationships between different lineages. The results of the phylogenomic analysis support a pattern of co-evolution between some of the ff.spp. and their host plant. In addition, we identified exceptions to this pattern, namely host jump events and the recent radiation of a clade less than 280,000 years ago. Furthermore, we found a high level of gene tree incongruence localized in the youngest clade. To distinguish between incomplete lineage sorting and lateral gene flow, we applied a coalescent-based method of demographic inference and found evidence of horizontal gene flow between recently diverged lineages. Overall we found that different processes shaped the diversification of B. graminis, co-evolution with the host species, host jump and fast radiation. Our study is an example of how genomic data can resolve complex evolutionary histories of cryptic lineages at different time scales, dealing with incomplete lineage sorting and lateral gene flow.

Oat evolution revealed in the maternal lineages of 25 Avena species

Cultivated hexaploid oat has three different sets of nuclear genomes (A, C, D), but its evolutionary history remains elusive. A multiplexed shotgun sequencing procedure was explored to acquire maternal phylogenetic signals from chloroplast and mitochondria genomes of 25 Avena species. Phylogenetic analyses of the acquired organelle SNP data revealed a new maternal pathway towards hexaploids of oat genome evolution involving three diploid species (A. ventricosa, A. canariensis and A. longiglumis) and two tetraploid species (A. insularis and A. agadiriana). Cultivated hexaploid A. sativa acquired its maternal genome from an AC genome tetraploid close to A. insularis. Both AC genome A. insularis and AB genome A. agadiriana obtained a maternal genome from an ancient A, not C, genome diploid close to A. longiglumis. Divergence dating showed the major divergences of C genome species 19.9–21.2 million years ago (Mya), of the oldest A genome A. canariensis 13–15 Mya, and of the clade with hexaploids 8.5–9.5 Mya. These findings not only advance our knowledge on oat genome evolution, but also have implications for oat germplasm conservation and utilization in breeding.

A 250 plastome phylogeny of the grass family (Poaceae): topological support under different data partitions

The systematics of grasses has advanced through applications of plastome phylogenomics, although studies have been largely limited to subfamilies or other subgroups of Poaceae. Here we present a plastome phylogenomic analysis of 250 complete plastomes (179 genera) sampled from 44 of the 52 tribes of Poaceae. Plastome sequences were determined from high throughput sequencing libraries and the assemblies represent over 28.7 Mbases of sequence data. Phylogenetic signal was characterized in 14 partitions, including (1) complete plastomes; (2) protein coding regions; (3) noncoding regions; and (4) three loci commonly used in single and multi-gene studies of grasses. Each of the four main partitions was further refined, alternatively including or excluding positively selected codons and also the gaps introduced by the alignment. All 76 protein coding plastome loci were found to be predominantly under purifying selection, but specific codons were found to be under positive selection in 65 loci. The loci that have been widely used in multi-gene phylogenetic studies had among the highest proportions of positively selected codons, suggesting caution in the interpretation of these earlier results. Plastome phylogenomic analyses confirmed the backbone topology for Poaceae with maximum bootstrap support (BP). Among the 14 analyses, 82 clades out of 309 resolved were maximally supported in all trees. Analyses of newly sequenced plastomes were in agreement with current classifications. Five of seven partitions in which alignment gaps were removed retrieved Panicoideae as sister to the remaining PACMAD subfamilies. Alternative topologies were recovered in trees from partitions that included alignment gaps. This suggests that ambiguities in aligning these uncertain regions might introduce a false signal. Resolution of these and other critical branch points in the phylogeny of Poaceae will help to better understand the selective forces that drove the radiation of the BOP and PACMAD clades comprising more than 99.9% of grass diversity.

Complete structure and variation of the chloroplast genome of Agropyron cristatum (L.) Gaertn

Agropyron cristatum (L.) Gaertner, a perennial grass in the tribe Triticeae (Poaceae), is a wild relative of cereal crops that is suitable for genetic improvement. In this study, we first sequenced the complete chloroplast (cp) genome of Ag. cristatum using Hiseq4000 PE150. The Ag. cristatum chloroplast genome is 135,554bp in length, has a typical quadripartite structure and contains 76 protein-coding genes, 29 tRNA genes and four rRNA genes. The cp genome of Ag. cristatum was used for comparison with other seven Triticeae species. One large variable region (800bp), which primarily contained the rpl23 (non-reciprocally translocated from IRs) and accD genes, was detected between rbcL gene and psaI gene within LSC region. The deletion of the accD and translocated rpl23 genes in Ag. cristatum indicated an independent gene-loss events or an additional divergence in Triticeae. Analyses of the dn/ds ratio and K2-P's genetic distance for 76 protein-coding genes showed that genes with evolutionary divergence might suffer from the effect of sequence regional constraints or gene functional constraints in Triticeae species. Our research will generally contribute to the knowledge of plastid genome evolution in Triticeae.

On the Origin of the Non-brittle Rachis Trait of Domesticated Einkorn Wheat

Einkorn and emmer wheat together with barley were among the first cereals domesticated by humans more than 10,000 years ago, long before durum or bread wheat originated. Domesticated einkorn wheat differs from its wild progenitor in basic morphological characters such as the grain dispersal system. This study identified the Non-brittle rachis 1 (btr1) and Non-brittle rachis 2 (btr2) in einkorn as homologous to barley. Re-sequencing of the Btr1 and Btr2 in a collection of 53 lines showed that a single non-synonymous amino acid substitution (alanine to threonine) at position 119 at btr1, is responsible for the non-brittle rachis trait in domesticated einkorn. Tracing this haplotype variation back to wild einkorn samples provides further evidence that the einkorn progenitor came from the Northern Levant. We show that the geographical origin of domesticated haplotype coincides with the non-brittle domesticated barley haplotypes, which suggest the non-brittle rachis phenotypes of einkorn and barley were fixed in same geographic area in today's South-east Turkey.

Gene Duplication and Evolution Dynamics in the Homeologous Regions Harboring Multiple Prolamin and Resistance Gene Families in Hexaploid Wheat

Improving end-use quality and disease resistance are important goals in wheat breeding. The genetic loci controlling these traits are highly complex, consisting of large families of prolamin and resistance genes with members present in all three homeologous A, B, and D genomes in hexaploid bread wheat. Here, orthologous regions harboring both prolamin and resistance gene loci were reconstructed and compared to understand gene duplication and evolution in different wheat genomes. Comparison of the two orthologous D regions from the hexaploid wheat Chinese Spring and the diploid progenitor Aegilops tauschii revealed their considerable difference due to the presence of five large structural variations with sizes ranging from 100 kb to 2 Mb. As a result, 44% of the Ae. tauschii and 71% of the Chinese Spring sequences in the analyzed regions, including 79 genes, are not shared. Gene rearrangement events, including differential gene duplication and deletion in the A, B, and D regions, have resulted in considerable erosion of gene collinearity in the analyzed regions, suggesting rapid evolution of prolamin and resistance gene families after the separation of the three wheat genomes. We hypothesize that this fast evolution is attributed to the co-evolution of the two gene families dispersed within a high recombination region. The identification of a full set of prolamin genes facilitated transcriptome profiling and revealed that the A genome contributes the least to prolamin expression because of its smaller number of expressed intact genes and their low expression levels, while the B and D genomes contribute similarly.

Evolution of the S-Genomes in Triticum-Aegilops Alliance: Evidences From Chromosome Analysis

Five diploid Aegilops species of the Sitopsis section: Ae. speltoides, Ae. longissima, Ae. sharonensis, Ae. searsii, and Ae. bicornis, two tetraploid species Ae. peregrina (= Ae. variabilis) and Ae. kotschyi (Aegilops section) and hexaploid Ae. vavilovii (Vertebrata section) carry the S-genomes. The B- and G-genomes of polyploid wheat are also the derivatives of the S-genome. Evolution of the S-genome species was studied using Giemsa C-banding and fluorescence in situ hybridization (FISH) with DNA probes representing 5S (pTa794) and 18S-5.8S-26S (pTa71) rDNAs as well as nine tandem repeats: pSc119.2, pAesp_SAT86, Spelt-1, Spelt-52, pAs1, pTa-535, and pTa-s53. To correlate the C-banding and FISH patterns we used the microsatellites (CTT)₁₀ and (GTT)₉, which are major components of the C-banding positive heterochromatin in wheat. According to the results obtained, diploid species split into two groups corresponding to Emarginata and Truncata sub-sections, which differ in the C-banding patterns, distribution of rDNA and other repeats. The B- and G-genomes of polyploid wheat are most closely related to the S-genome of Ae. speltoides. The genomes of allopolyploid wheat have been evolved as a result of different species-specific chromosome translocations, sequence amplification, elimination and re-patterning of repetitive DNA sequences. These events occurred independently in different wheat species and in Ae. speltoides _. The 5S rDNA locus of chromosome 1S was probably lost in ancient Ae. speltoides prior to formation of Timopheevii wheat, but after the emergence of ancient emmer. Evolution of Emarginata species was associated with an increase of C-banding and (CTT)₁₀-positive heterochromatin, amplification of Spelt-52, re-pattering of the pAesp_SAT86, and a gradual decrease in the amount of the D-genome-specific repeats pAs1, pTa-535, and pTa-s53. The emergence of Ae. peregrina and Ae. kotschyi did not lead to significant changes of the S^*-genomes. However, partial elimination of 45S rDNA repeats from 5S^* and 6S^* chromosomes and alterations of C-banding and FISH-patterns have been detected. Similarity of the S^v-genome of Ae. vavilovii with the S^s genome of diploid Ae. searsii confirmed the origin of this hexaploid. A model of the S-genome evolution is suggested.

Conserved molecular structure of the centromeric histone CENH3 in Secale and its phylogenetic relationships

It has been repeatedly demonstrated that the centromere-specific histone H3 (CENH3), a key component of the centromere, shows considerable variability between species within taxa. We determined the molecular structure and phylogenetic relationships of CENH3 in 11 Secale species and subspecies that possess distinct pollination systems and are adapted to a wide range of abiotic and biotic stresses. The rye (Secale cereale) genome encodes two paralogous CENH3 genes, which differ in intron-exon structure and are transcribed into two main forms of the protein, αCENH3 and βCENH3. These two forms differ in size and amino acid substitutions. In contrast to the reported differences in CENH3 structure between species within other taxa, the main forms of this protein in Secale species and subspecies have a nearly identical structure except some nonsynonymous substitutions. The CENH3 proteins are strictly controlled by genetic factors responsible for purifying selection. A comparison between Hordeum, Secale and Triticum species demonstrates that the structure of CENH3 in the subtribes Hordeinae and Triticinae evolved at different rates. The assumption that reticulate evolution served as a factor stabilizing the structure and evolutionary rate of CENH3 and that this factor was more powerful within Secale and Triticum than in Hordeum, is discussed.

Chromosomal assignment of centromere-specific histone CENH3 genes in rye (Secale cereale L.) and their phylogeny

Centromeres are essential for correct chromosome segregation during cell division and are determined by the presence of centromere-specific histone 3 (CENH3). Most of the diploid plant species, in which the structure and copy number of CENH3 genes have been determined, have this gene as a singleton; however, some cereal species in the tribe Triticeae have been found to have CENH3 in two variants. In this work, using the set of the wheat-rye addition lines we wanted to establish the chromosomal assignment of the CENH3 genes in the cultivated rye, Secale cereale (Linnaeus, 1753), in order to expand our knowledge about synteny conservation in the most important cereal species and about their chromosome evolution. To this end, we have also analyzed data in available genome sequencing databases. As a result, the αCENH3 and βCENH3 forms have been assigned to rye chromosomes 1R and 6R: specifically, the commonest variants αCENH3v1 and βCENH3v1 to chromosome 1R, and the rare variants, αCENH3v2 and probably βCENH3v2, to chromosome 6R. No other CENH3 variants have been found by analysis of the rye genome sequencing databases. Our chromosomal assignment of CENH3 in rye has been found to be the same as that in barley, suggesting that both main forms of CENH3 appeared in a Triticeae species before the barley and wheatrye lineages split.

Integral Phylogenomic Approach over Ilex L. Species from Southern South America

The use of molecular markers with inadequate variation levels has resulted in poorly resolved phylogenetic relationships within Ilex. Focusing on southern South American and Asian species, we aimed at contributing informative plastid markers. Also, we intended to gain insights into the nature of morphological and physiological characters used to identify species. We obtained the chloroplast genomes of I.paraguariensis and I. dumosa, and combined these with all the congeneric plastomes currently available to accomplish interspecific comparisons and multilocus analyses. We selected seven introns and nine IGSs as variable non-coding markers that were used in phylogenomic analyses. Eight extra IGSs were proposed as candidate markers. Southern South American species formed one lineage, except for I. paraguariensis, I. dumosa and I. argentina, which occupied intermediate positions among sampled taxa; Euroasiatic species formed two lineages. Some concordant relationships were retrieved from nuclear sequence data. We also conducted integral analyses, involving a supernetwork of molecular data, and a simultaneous analysis of quantitative and qualitative morphological and phytochemical characters, together with molecular data. The total evidence tree was used to study the evolution of non-molecular data, evidencing fifteen non-ambiguous synapomorphic character states and consolidating the relationships among southern South American species. More South American representatives should be incorporated to elucidate their origin.

Dated tribe-wide whole chloroplast genome phylogeny indicates recurrent hybridizations within Triticeae

BACKGROUND

Triticeae, the tribe of wheat grasses, harbours the cereals barley, rye and wheat and their wild relatives. Although economically important, relationships within the tribe are still not understood. We analysed the phylogeny of chloroplast lineages among nearly all monogenomic Triticeae taxa and polyploid wheat species aiming at a deeper understanding of the tribe's evolution. We used on- and off-target reads of a target-enrichment experiment followed by Illumina sequencing.

RESULTS

The read data was used to assemble the plastid locus ndhF for 194 individuals and the whole chloroplast genome for 183 individuals, representing 53 Triticeae species and 15 genera. We conducted Bayesian and multispecies coalescent analyses to infer relationships and estimate divergence times of the taxa. We present the most comprehensive dated Triticeae chloroplast phylogeny and review previous hypotheses in the framework of our results. Monophyly of Triticeae chloroplasts could not be confirmed, as either Bromus or Psathyrostachys captured a chloroplast from a lineage closely related to a Bromus-Triticeae ancestor. The most recent common ancestor of Triticeae occurred approximately between ten and 19 million years ago.

CONCLUSIONS

The comparison of the chloroplast phylogeny with available nuclear data in several cases revealed incongruences indicating past hybridizations. Recent events of chloroplast capture were detected as individuals grouped apart from con-specific accessions in otherwise monopyhletic groups.

A systematic review of rye (Secale cereale L.) as a source of resistance to pathogens and pests in wheat (Triticum aestivum L.)

Wheat is globally one of the most important crops. With the current human population growth rate, there is an increasing need to raise wheat productivity by means of plant breeding, along with development of more efficient and sustainable agricultural systems. Damage by pathogens and pests, in combination with adverse climate effects, need to be counteracted by incorporating new germplasm that makes wheat more resistant/tolerant to such stress factors. Rye has been used as a source for improved resistance to pathogens and pests in wheat during more than 50 years. With new devastating stem and yellow rust pathotypes invading wheat at large acreage globally, along with new biotypes of pest insects, there is renewed interest in using rye as a source of resistance. Currently the proportion of wheat cultivars with rye chromatin varies between countries, with examples of up to 34%. There is mainly one rye source, Petkus, that has been widely exploited and that has contributed considerably to raise yields and increase disease resistance in wheat. Successively, the multiple disease resistances conferred by this source has been overcome by new pathotypes of leaf rust, yellow rust, stem rust and powdery mildew. However, there are several other rye sources reported to make wheat more resistant to various biotic constraints when their rye chromatin has been transferred to wheat. There is also development of knowledge on how to produce new rye translocation, substitution and addition lines. Here we compile information that may facilitate decision making for wheat breeders aiming to transfer resistance to biotic constraints from rye to elite wheat germplasm.

Wheat Ms2 encodes for an orphan protein that confers male sterility in grass species

Male sterility is a valuable trait for plant breeding and hybrid seed production. The dominant male-sterile gene Ms2 in common wheat has facilitated the release of hundreds of breeding lines and cultivars in China. Here, we describe the map-based cloning of the Ms2 gene and show that Ms2 confers male sterility in wheat, barley and Brachypodium. MS2 appears as an orphan gene within the Triticinae and expression of Ms2 in anthers is associated with insertion of a retroelement into the promoter. The cloning of Ms2 has substantial potential to assemble practical pipelines for recurrent selection and hybrid seed production in wheat.

Long noncoding miRNA gene represses wheat β-diketone waxes

The cuticle of terrestrial plants functions as a protective barrier against many biotic and abiotic stresses. In wheat and other Triticeae, β-diketone waxes are major components of the epicuticular layer leading to the bluish-white glaucous trait in reproductive-age plants. Glaucousness in durum wheat is controlled by a metabolic gene cluster at the WAX1 (W1) locus and a dominant suppressor INHIBITOR of WAX1 (Iw1) on chromosome 2B. The wheat D subgenome from progenitor Aegilops tauschii contains W2 and Iw2 paralogs on chromosome 2D. Here we identify the Iw1 gene from durum wheat and demonstrate the unique regulatory mechanism by which Iw1 acts to suppress a carboxylesterase-like protein gene, W1-COE, within the W1 multigene locus. Iw1 is a long noncoding RNA (lncRNA) containing an inverted repeat (IR) with >80% identity to W1-COE The Iw1 transcript forms a miRNA precursor-like long hairpin producing a 21-nt predominant miRNA, miRW1, and smaller numbers of related sRNAs associated with the nonglaucous phenotype. When Iw1 was introduced into glaucous bread wheat, miRW1 accumulated, W1-COE and its paralog W2-COE were down-regulated, and the phenotype was nonglaucous and β-diketone-depleted. The IR region of Iw1 has >94% identity to an IR region on chromosome 2 in Ae. tauschii that also produces miRW1 and lies within the marker-based location of Iw2 We propose the Iw loci arose from an inverted duplication of W1-COE and/or W2-COE in ancestral wheat to form evolutionarily young miRNA genes that act to repress the glaucous trait.

AvrPm2 encodes an RNase-like avirulence effector which is conserved in the two different specialized forms of wheat and rye powdery mildew fungus

There is a large diversity of genetically defined resistance genes in bread wheat against the powdery mildew pathogen Blumeria graminis (B. g.) f. sp. tritici. Many confer race-specific resistance to this pathogen, but until now only the mildew avirulence gene AvrPm3a2/f2 that is recognized by Pm3a/f was known molecularly. We performed map-based cloning and genome-wide association studies to isolate a candidate for the mildew avirulence gene AvrPm2. We then used transient expression assays in Nicotiana benthamiana to demonstrate specific and strong recognition of AvrPm2 by Pm2. The virulent AvrPm2 allele arose from a conserved 12 kb deletion, while there is no protein sequence diversity in the gene pool of avirulent B. g. tritici isolates. We found one polymorphic AvrPm2 allele in B. g. triticale and one orthologue in B. g. secalis and both are recognized by Pm2. AvrPm2 belongs to a small gene family encoding structurally conserved RNase-like effectors, including Avra13 from B. g. hordei, the cognate Avr of the barley resistance gene Mla13. These results demonstrate the conservation of functional avirulence genes in two cereal powdery mildews specialized on different hosts, thus providing a possible explanation for successful introgression of resistance genes from rye or other grass relatives to wheat.

Alignment of Common Wheat and Other Grass Genomes Establishes a Comparative Genomics Research Platform

Grass genomes are complicated structures as they share a common tetraploidization, and particular genomes have been further affected by extra polyploidizations. These events and the following genomic re-patternings have resulted in a complex, interweaving gene homology both within a genome, and between genomes. Accurately deciphering the structure of these complicated plant genomes would help us better understand their compositional and functional evolution at multiple scales. Here, we build on our previous research by performing a hierarchical alignment of the common wheat genome vis-à-vis eight other sequenced grass genomes with most up-to-date assemblies, and annotations. With this data, we constructed a list of the homologous genes, and then, in a layer-by-layer process, separated their orthology, and paralogy that were established by speciations and recursive polyploidizations, respectively. Compared with the other grasses, the far fewer collinear outparalogous genes within each of three subgenomes of common wheat suggest that homoeologous recombination, and genomic fractionation should have occurred after its formation. In sum, this work contributes to the establishment of an important and timely comparative genomics platform for researchers in the grass community and possibly beyond. Homologous gene list can be found in Supplemental material.

Genome origin, historical hybridization and genetic differentiation in Anthosachne australasica (Triticeae; Poaceae), inferred from chloroplast rbcL, trnH-psbA and nuclear Acc1 gene sequences

BACKGROUND AND AIMS

Anthosachne Steudel is a group of allopolyploid species that was derived from hexaploidization between the Asian StY genome Roegneria entity and the Australasia W genome Australopyrum species. Polyploidization and apomixis contribute to taxonomic complexity in Anthosachne Here, a study is presented on the phylogeny and evolutionary history of Anthosachne australasica The aims are to demonstrate the process of polyploidization events and to explore the differentiation patterns of the St genome following geographic isolation.

METHODS

Chloroplast rbcL and trnH-psbA and nuclear Acc1 gene sequences of 60 Anthosachne taxa and nine Roegneria species were analysed with those of 33 diploid taxa representing 20 basic genomes in Triticeae. The phylogenetic relationships were reconstructed. A time-calibrated phylogeny was generated to estimate the evolutionary history of A. australasica Nucleotide diversity patterns were used to assess the divergence within A. australasica and between Anthosachne and its putative progenitors.

KEY RESULTS

Three homoeologous copies of the Acc1 sequences from Anthosachne were grouped with the Acc1 sequences from Roegneria, Pseudoroegneria, Australopyrum, Dasypyrum and Peridictyon The chloroplast sequences of Anthosachne were clustered with those from Roegneria and Pseudoroegneria Divergence time for Anthosachne was dated to 4·66 million years ago (MYA). The level of nucleotide diversity in Australasian Anthosachne was higher than that in continental Roegneria A low level of genetic differentiation within the A. australasica complex was found.

CONCLUSIONS

Anthosachne originated from historical hybridization between Australopyrum species and a Roegneria entity colonized from Asia to Australasia via South-east Asia during the late Miocene. The St lineage served as the maternal donor during the speciation of Anthosachne A contrasting pattern of population genetic structure exists in the A. australasica complex. Greater diversity in island Anthosachne compared with continental Roegneria might be associated with mutation, polyploidization, apomixis and expansion. It is reasonable to consider that A. australasica var. scabra and A. australasica var. plurinervisa should be included in the A. australasica complex.

Identification and Analysis of RNA Editing Sites in the Chloroplast Transcripts of Aegilops tauschii L

RNA editing is an important way to convert cytidine (C) to uridine (U) at specific sites within RNA molecules at a post-transcriptional level in the chloroplasts of higher plants. Although it has been systematically studied in many plants, little is known about RNA editing in the wheat D genome donor Aegilops tauschii L. Here, we investigated the chloroplast RNA editing of Ae. tauschii and compared it with other wheat relatives to trace the evolution of wheat. Through bioinformatics prediction, a total of 34 C-to-U editing sites were identified, 17 of which were validated using RT-PCR product sequencing. Furthermore, 60 sites were found by the RNA-Seq read mapping approach, 24 of which agreed with the prediction and six were validated experimentally. The editing sites were biased toward tCn or nCa trinucleotides and 5′-pyrimidines, which were consistent with the flanking bases of editing sites of other seed plants. Furthermore, the editing events could result in the alteration of the secondary structures and topologies of the corresponding proteins, suggesting that RNA editing might impact the function of target genes. Finally, comparative analysis found some evolutionarily conserved editing sites in wheat and two species-specific sites were also obtained. This study is the first to report on RNA editing in Aegilops tauschii L, which not only sheds light on the evolution of wheat from the point of view of RNA editing, but also lays a foundation for further studies to identify the mechanisms of C-to-U alterations.

Mitochondrial genome sequences from wild and cultivated barley (Hordeum vulgare)

Background

Sequencing analysis of mitochondrial genomes is important for understanding the evolution and genome structures of various plant species. Barley is a self-pollinated diploid plant with seven chromosomes comprising a large haploid genome of 5.1 Gbp. Wild barley (Hordeum vulgare ssp. spontaneum) and cultivated barley (H. vulgare ssp. vulgare) have cross compatibility and closely related genomes, although a significant number of nucleotide polymorphisms have been reported between their genomes.

Results

We determined the complete nucleotide sequences of the mitochondrial genomes of wild and cultivated barley. Two independent circular maps of the 525,599 bp barley mitochondrial genome were constructed by de novo assembly of high-throughput sequencing reads of barley lines H602 and Haruna Nijo, with only three SNPs detected between haplotypes. These mitochondrial genomes contained 33 protein-coding genes, three ribosomal RNAs, 16 transfer RNAs, 188 new ORFs, six major repeat sequences and several types of transposable elements. Of the barley mitochondrial genome-encoded proteins, NAD6, NAD9 and RPS4 had unique structures among grass species.

Conclusions

The mitochondrial genome of barley was similar to those of other grass species in terms of gene content, but the configuration of the genes was highly differentiated from that of other grass species. Mitochondrial genome sequencing is essential for annotating the barley nuclear genome; our mitochondrial sequencing identified a significant number of fragmented mitochondrial sequences in the reported nuclear genome sequences. Little polymorphism was detected in the barley mitochondrial genome sequences, which should be explored further to elucidate the evolution of barley.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3159-3) contains supplementary material, which is available to authorized users.

Allelic barley MLA immune receptors recognize sequence-unrelated avirulence effectors of the powdery mildew pathogen

Disease-resistance genes encoding intracellular nucleotide-binding domain and leucine-rich repeat proteins (NLRs) are key components of the plant innate immune system and typically detect the presence of isolate-specific avirulence (AVR) effectors from pathogens. NLR genes define the fastest-evolving gene family of flowering plants and are often arranged in gene clusters containing multiple paralogs, contributing to copy number and allele-specific NLR variation within a host species. Barley mildew resistance locus a (Mla) has been subject to extensive functional diversification, resulting in allelic resistance specificities each recognizing a cognate, but largely unidentified, AVR_a gene of the powdery mildew fungus, Blumeria graminis f. sp. hordei (Bgh). We applied a transcriptome-wide association study among 17 Bgh isolates containing different AVR_a genes and identified AVR_a1 and AVR_a13, encoding candidate-secreted effectors recognized by Mla1 and Mla13 alleles, respectively. Transient expression of the effector genes in barley leaves or protoplasts was sufficient to trigger Mla1 or Mla13 allele-specific cell death, a hallmark of NLR receptor-mediated immunity. AVR_a1 and AVR_a13 are phylogenetically unrelated, demonstrating that certain allelic MLA receptors evolved to recognize sequence-unrelated effectors. They are ancient effectors because corresponding loci are present in wheat powdery mildew. AVR_A1 recognition by barley MLA1 is retained in transgenic Arabidopsis, indicating that AVR_A1 directly binds MLA1 or that its recognition involves an evolutionarily conserved host target of AVR_A1 Furthermore, analysis of transcriptome-wide sequence variation among the Bgh isolates provides evidence for Bgh population structure that is partially linked to geographic isolation.

Molecular genetics and evolution of disease resistance in cereals

Contents 320 I. 320 II. 321 III. 321 IV. 322 V. 324 VI. 328 VII. 329 330 References 330 SUMMARY: Cereal crops produce a large part of the globally consumed food and feed. Because of the constant presence of devastating pathogens, the molecular characterization of disease resistance is a major research area and highly relevant for breeding. There has been recent and accelerating progress in the understanding of three distinct resistance mechanisms in cereals: resistance conferred by plasma membrane-localized receptor proteins; race-specific resistance conferred by intracellular immune receptors; and quantitative disease resistance. Intracellular immune receptors provide a particularly rich source for evolutionary studies, and have, for example, resulted in the recent discovery of a novel detection mechanism based on integrated decoy domains. Evolutionary studies have also revealed the origins of active resistance genes in both wild progenitors of today's cereals as well as in cultivated forms. In addition, independent evolution of orthologous genes in related cereals has resulted in resistance to different pathogen species. Quantitative resistance genes have been best characterized in wheat. The quantitative resistance genes identified so far in wheat encode transporter proteins or unusual kinase proteins. The recent discoveries in these three different resistance mechanisms have contributed to the basic molecular understanding of cereal immunity against pathogens and have suggested novel applications for resistance breeding.

DNA transposon activity is associated with increased mutation rates in genes of rice and other grasses

DNA (class 2) transposons are mobile genetic elements which move within their ‘host' genome through excising and re-inserting elsewhere. Although the rice genome contains tens of thousands of such elements, their actual role in evolution is still unclear. Analysing over 650 transposon polymorphisms in the rice species Oryza sativa and Oryza glaberrima, we find that DNA repair following transposon excisions is associated with an increased number of mutations in the sequences neighbouring the transposon. Indeed, the 3,000 bp flanking the excised transposons can contain over 10 times more mutations than the genome-wide average. Since DNA transposons preferably insert near genes, this is correlated with increases in mutation rates in coding sequences and regulatory regions. Most importantly, we find this phenomenon also in maize, wheat and barley. Thus, these findings suggest that DNA transposon activity is a major evolutionary force in grasses which provide the basis of most food consumed by humankind.

DNA transposons are numerous in plant genomes. Here, Wicker et al. analyse transposon polymorphisms in rice and other grasses and show that sequences flanking excision sites contain up to 10 times more mutations than average, suggesting transposons are a major factor shaping the evolution of grass genomes.

Chloroplast genomes: diversity, evolution, and applications in genetic engineering

Chloroplasts play a crucial role in sustaining life on earth. The availability of over 800 sequenced chloroplast genomes from a variety of land plants has enhanced our understanding of chloroplast biology, intracellular gene transfer, conservation, diversity, and the genetic basis by which chloroplast transgenes can be engineered to enhance plant agronomic traits or to produce high-value agricultural or biomedical products. In this review, we discuss the impact of chloroplast genome sequences on understanding the origins of economically important cultivated species and changes that have taken place during domestication. We also discuss the potential biotechnological applications of chloroplast genomes.

The Cer-cqu gene cluster determines three key players in a β-diketone synthase polyketide pathway synthesizing aliphatics in epicuticular waxes

Aliphatic compounds on plant surfaces, called epicuticular waxes, are the first line of defense against pathogens and pests, contribute to reducing water loss and determine other important phenotypes. Aliphatics can form crystals affecting light refraction, resulting in a color change and allowing identification of mutants in their synthesis or transport. The present study discloses three such Eceriferum (cer) genes in barley - Cer-c, Cer-q and Cer-u - known to be tightly linked and functioning in a biochemical pathway forming dominating amounts of β-diketone and hydroxy-β-diketones plus some esterified alkan-2-ols. These aliphatics are present in many Triticeae as well as dicotyledons such as Eucalyptus and Dianthus. Recently developed genomic resources and mapping populations in barley defined these genes to a small region on chromosome arm 2HS. Exploiting Cer-c and -u potential functions pinpointed five candidates, of which three were missing in apparent cer-cqu triple mutants. Sequencing more than 50 independent mutants for each gene confirmed their identification. Cer-c is a chalcone synthase-like polyketide synthase, designated diketone synthase (DKS), Cer-q is a lipase/carboxyl transferase and Cer-u is a P450 enzyme. All were highly expressed in pertinent leaf sheath tissue of wild type. A physical map revealed the order Cer-c, Cer-u, Cer-q with the flanking genes 101kb apart, confirming they are a gene cluster, Cer-cqu. Homology-based modeling suggests that many of the mutant alleles affect overall protein structure or specific active site residues. The rich diversity of identified mutations will facilitate future studies of three key enzymes involved in synthesis of plant apoplast waxes.