High-copy transposons from a pathogen give rise to a conserved sRNA family with a novel host immunity target.

Small RNAs (sRNAs) are involved in gene silencing in multiple ways including through cross-kingdom transfers from parasites to their hosts. Little is known about the evolutionary mechanisms enabling eukaryotic microbes to evolve functional mimics of host small regulatory RNAs. Here, we describe the identification and functional characterization of SINE_sRNA1, a sRNA family derived from highly abundant SINE retrotransposons in the genome of the wheat powdery mildew pathogen. SINE_sRNA1 is encoded by a sequence motif that is conserved in multiple SINE families and corresponds to a functional plant miRNA mimic targeting Tae_AP1, a wheat gene encoding an aspartic protease only found in monocots. Tae_AP1 has a novel function enhancing both pattern-triggered immunity (PTI) and effector-triggered immunity (ETI) thus contributing to the cross activation of plant defenses. We conclude that SINE_sRNA1 and Tae_AP1 are functional innovations suggesting the contribution of transposons to the evolutionary arms race between a parasite and its host.

A survey of lineage-specific genes in Triticeae reveals de novo gene evolution from genomic raw material.

Diploid plant genomes typically contain ~35,000 genes, almost all belonging to highly conserved gene families. Only a small fraction are lineage-specific, which are found in only one or few closely related species. Little is known about how genes arise de novo in plant genomes and how often this occurs; however, they are believed to be important for plants diversification and adaptation. We developed a pipeline to identify lineage-specific genes in Triticeae, using newly available genome assemblies of wheat, barley, and rye. Applying a set of stringent criteria, we identified 5942 candidate Triticeae-specific genes (TSGs), of which 2337 were validated as protein-coding genes in wheat. Differential gene expression analyses revealed that stress-induced wheat TSGs are strongly enriched in putative secreted proteins. Some were previously described to be involved in Triticeae non-host resistance and cold response. Additionally, we show that 1079 TSGs have sequence homology to transposable elements (TEs), ~68% of them deriving from regulatory non-coding regions of Gypsy retrotransposons. Most importantly, we demonstrate that these TSGs are enriched in transmembrane domains and are among the most highly expressed wheat genes overall. To summarize, we conclude that de novo gene formation is relatively rare and that Triticeae probably possess ~779 lineage-specific genes per haploid genome. TSGs, which respond to pathogen and environmental stresses, may be interesting candidates for future targeted resistance breeding in Triticeae. Finally, we propose that non-coding regions of TEs might provide important genetic raw material for the functional innovation of TM domains and the evolution of novel secreted proteins.

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.

Comparative Transcriptome Analysis of Wheat Lines in the Field Reveals Multiple Essential Biochemical Pathways Suppressed by Obligate Pathogens.

Mildew and rust are the most devastating cereal pathogens, and in wheat they can cause up to 50% yield loss every year. Wheat lines containing resistance genes are used to effectively control fungal diseases, but the molecular mechanisms underlying the interaction between wheat and its fungal pathogens are poorly understood. Here, we used RNA sequencing (RNA-Seq) to compare the transcriptomic landscape of susceptible and resistant wheat lines to identify genes and pathways that are targeted by obligate biotrophic fungal pathogens. The five lines differed in the expression of thousands of genes under infection as well as control conditions. Generally, mixed infection with powdery mildew and leaf rust resulted in downregulation of numerous genes in susceptible lines. Interestingly, transcriptomic comparison between the nearly isogenic lines Thatcher and Thatcher-Lr34 identified 753 genes that are uniquely downregulated in the susceptible line upon infection. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis, revealed the suppression of six major biochemical pathways, namely nuclear transport, alternative splicing, DNA damage response, ubiquitin-mediated proteolysis, phosphoinositol signaling, and photosynthesis. We conclude that powdery mildew and leaf rust evade the wheat defense system by suppression of programmed cell death (PCD) and responses to cellular damage. Considering the broad range of the induced changes, we propose that the pathogen targets "master regulators" at critical steps in the respective pathways. Identification of these wheat genes targeted by the pathogen could inspire new directions for future wheat breeding.

A membrane-bound ankyrin repeat protein confers race-specific leaf rust disease resistance in wheat.

Plasma membrane-associated and intracellular proteins and protein complexes play a pivotal role in pathogen recognition and disease resistance signaling in plants and animals. The two predominant protein families perceiving plant pathogens are receptor-like kinases and nucleotide binding-leucine-rich repeat receptors (NLR), which often confer race-specific resistance. Leaf rust is one of the most prevalent and most devastating wheat diseases. Here, we clone the race-specific leaf rust resistance gene Lr14a from hexaploid wheat. The cloning of Lr14a is aided by the recently published genome assembly of ArinaLrFor, an Lr14a-containing wheat line. Lr14a encodes a membrane-localized protein containing twelve ankyrin (ANK) repeats and structural similarities to Ca2+-permeable non-selective cation channels. Transcriptome analyses reveal an induction of genes associated with calcium ion binding in the presence of Lr14a. Haplotype analyses indicate that Lr14a-containing chromosome segments were introgressed multiple times into the bread wheat gene pool, but we find no variation in the Lr14a coding sequence itself. Our work demonstrates the involvement of an ANK-transmembrane (TM)-like type of gene family in race-specific disease resistance in wheat. This forms the basis to explore ANK-TM-like genes in disease resistance breeding.

Domestication of High-Copy Transposons Underlays the Wheat Small RNA Response to an Obligate Pathogen.

Plant genomes have evolved several evolutionary mechanisms to tolerate and make use of transposable elements (TEs). Of these, transposon domestication into cis-regulatory and microRNA (miRNA) sequences is proposed to contribute to abiotic/biotic stress adaptation in plants. The wheat genome is derived at 85% from TEs, and contains thousands of miniature inverted-repeat transposable elements (MITEs), whose sequences are particularly prone for domestication into miRNA precursors. In this study, we investigate the contribution of TEs to the wheat small RNA immune response to the lineage-specific, obligate powdery mildew pathogen. We show that MITEs of the Mariner superfamily contribute the largest diversity of miRNAs to the wheat immune response. In particular, MITE precursors of miRNAs are wide-spread over the wheat genome, and highly conserved copies are found in the Lr34 and QPm.tut-4A mildew resistance loci. Our work suggests that transposon domestication is an important evolutionary force driving miRNA functional innovation in wheat immunity.

A chromosome-scale genome assembly reveals a highly dynamic effector repertoire of wheat powdery mildew.

Blumeria graminis f. sp. tritici (B.g. tritici) is the causal agent of the wheat powdery mildew disease. The highly fragmented B.g. tritici genome available so far has prevented a systematic analysis of effector genes that are known to be involved in host adaptation. To study the diversity and evolution of effector genes we produced a chromosome-scale assembly of the B.g. tritici genome. The genome assembly and annotation was achieved by combining long-read sequencing with high-density genetic mapping, bacterial artificial chromosome fingerprinting and transcriptomics. We found that the 166.6 Mb B.g. tritici genome encodes 844 candidate effector genes, over 40% more than previously reported. Candidate effector genes have characteristic local genomic organization such as gene clustering and enrichment for recombination-active regions and certain transposable element families. A large group of 412 candidate effector genes shows high plasticity in terms of copy number variation in a global set of 36 isolates and of transcription levels. Our data suggest that copy number variation and transcriptional flexibility are the main drivers for adaptation in B.g. tritici. The high repeat content may play a role in providing a genomic environment that allows rapid evolution of effector genes with selection as the driving force.