Paleo-evolutionary plasticity of plant disease resistance genes

Integrated analysis of Populus deltoides PR1 genes uncovered a PdePR1 as a defense marker against foliar rust

Pathogenesis-related protein 1 (PR1), a hallmark of plant disease resistance, plays pivotal roles in defense signaling. In this study, we identified 16 intronless PR1 genes in Populus deltoides, all classified within the CAP superfamily (cysteine-rich secretory protein, antigen 5, and pathogenesis-related 1) and characterized by conserved N-terminal signal peptides, caveolin-binding motifs, and CAP-derived peptides. Phylogenomic reconstruction of 231 PR1 homologs across 15 plant species traced their origin to Chara braunii, with lineage-specific expansions driven by gene duplication. Evolutionary analyses revealed strong purifying selection acting on ancestral PR1 paralogs to confer a selective advantage for disease resistance. Integrated transcriptomic profiling and quantitative RT-PCR analyses identified PdePR1_10 as a key marker gene for defense activation, exhibiting significant induction at two days post-inoculation in resistant poplars. Co-expression network analysis indicated that PdePR1_10 interacts with several defense-related genes, including NBS-LRR resistance genes, signaling kinases, and hormone biosynthesis enzymes. Specifically, the W-box cis-regulatory element in the PdePR1_10 promoter is hypothesized to interact with WRKY transcription factors, activating PdePR1_10 expression through a salicylic acid (SA)-dependent signaling pathway. Transgenic poplars overexpressing PdePR1_10 exhibited significantly enhanced rust resistance, confirming its critical in defense response. In summary, we thoroughly elucidated the biological functions and regulatory mechanisms of PR1 genes in rust resistance and provided a valuable transgenic poplar line for future studies.

Evolutionary signatures of the erosion of sexual reproduction genes in domesticated cassava (Manihot esculenta)

Centuries of clonal propagation in cassava (Manihot esculenta) have reduced sexual recombination, leading to the accumulation of deleterious mutations. This has resulted in both inbreeding depression affecting yield and a significant decrease in reproductive performance, creating hurdles for contemporary breeding programs. Cassava is a member of the Euphorbiaceae family, including notable species such as rubber tree (Hevea brasiliensis) and poinsettia (Euphorbia pulcherrima). Expanding upon preliminary draft genomes, we annotated 7 long-read genome assemblies and aligned a total of 52 genomes, to analyze selection across the genome and the phylogeny. Through this comparative genomic approach, we identified 48 genes under relaxed selection in cassava. Notably, we discovered an overrepresentation of floral expressed genes, especially focused at 6 pollen-related genes. Our results indicate that domestication and a transition to clonal propagation have reduced selection pressures on sexually reproductive functions in cassava leading to an accumulation of mutations in pollen-related genes. This relaxed selection and the genome-wide deleterious mutations responsible for inbreeding depression are potential targets for improving cassava breeding, where the generation of new varieties relies on recombining favorable alleles through sexual reproduction.

Phylogenomic analysis reveals exceptions to the co-evolution of ZAR1 and ZRK immune gene families in plants

BACKGROUND

HOPZ-ACTIVATED RESISTANCE 1 (ZAR1) is a nucleotide-binding leucine-rich repeat (NLR) protein functioning as a recognition hub to initiate effector-triggered immunity against bacterial pathogens. To initiate defense, ZAR1 associates with different HOPZ-ETI-DEFICIENT 1 (ZED1)-Related Kinases (ZRKs) to form resistosomes to indirectly perceive effector-induced perturbations. Few studies have focused on the phylogenomic characteristics of ZAR1 and ZRK immune gene families and their evolutionary relationships. To trace the origin and divergence of ZAR1 and ZRK immune gene families across the plant kingdom, we performed phylogenomic analyses using an extended set of plant genomes.

RESULTS

Genome-wide identification of ZAR1 and ZRK immune gene families by blast similarity searches combined with phylogenetic analysis showed that these two gene families have experienced frequent gene losses in massive lineages. Gene distribution patterns across the plant kingdom revealed that ZAR1 and ZRK emerged after the divergence of most angiosperms from Amborella and before the split of magnoliids, monocots, and eudicots. Co-occurrence of ZAR1-A and ZRKs was found in various plant species belonging to different angiosperm orders, but both genes were found to be absent in chlorophyta, bryophytes, lycophytes, ferns, and gymnosperms. We also detected a large number of concerted gene losses in angiosperms, especially within the orders Fabales, Cucurbitales, Asterales, and Apiales. All analysed monocot genomes thus far examined, except for the aroid Colocasia esculenta, were previously reported to lack both ZAR1-A and ZRKs. Here we now report other exceptions on the concerted ZAR1-A-ZRKs presence-absence pattern within several early diverging monocot lineages, including the genome of Acorus tatarinowii-a species representing the first branching monocot lineage. We also revealed strong variation in ZAR1-A-ZRKs co-occurrence within the asterid order Ericales, suggesting patterns of de-coevolution in angiosperms. Our research further shows that both gene families experienced significant diversification through various duplication events. Additionally, their evolutionary paths have been shaped by frequent gene losses and lineage-specific transposition.

CONCLUSION

This study provides novel findings on the evolution of ZAR1 and ZRK immune gene families across a wide range of plant species, suggesting that more potential exceptions can be expected when expanding the list of sequenced genomes from distinct orders. Our results provide new hypotheses about the origin and diversification of these critical immune genes for future functional studies.

Integrative Analysis of Transcriptome and Metabolome Reveals the Pivotal Role of the NAM Family Genes in Oncidium hybridum Lodd. Pseudobulb Growth

Oncidium hybridum Lodd. is an important ornamental flower that is used as both a cut flower and a potted plant around the world; additionally, its pseudobulbs serve as essential carriers for floral organs and flower development. The NAM gene family is crucial for managing responses to various stresses as well as regulating growth in plants. However, the mechanisms by which NAM genes regulate the development of pseudobulbs remain unclear. In this study, a total of 144 NAM genes harboring complete structural domains were identified in O. hybridum. The 144 NAM genes were systematically classified into 14 distinct subfamilies via phylogenetic analysis. Delving deeper into the conserved motifs revealed that motifs 1-6 exhibited remarkable conservation, while motifs 7-10 presented in a few NAM genes only. Notably, NAM genes sharing identical specific motifs were classified into the same subfamily, indicating functional relatedness. Furthermore, the examination of occurrences of gene duplication indicated that the NAM genes display 16 pairs of tandem duplications along with five pairs of segmental duplications, suggesting their role in genetic diversity and potential adaptive evolution. By conducting a correlation analysis integrating transcriptomics and metabolomics at four stages of pseudobulb development, we found that OhNAM023, OhNAM030, OhNAM007, OhNAM019, OhNAM083, OhNAM047, OhNAM089, and OhNAM025 exhibited significant relationships with the endogenous plant hormones jasmonates (JAs), hinting at their potential involvement in hormonal signaling. Additionally, OhNAM089, OhNAM025, OhNAM119, OhNAM055, and OhNAM136 showed strong links with abscisic acid (ABA) and abscisic acid glucose ester (ABA-GE), suggesting the possible regulatory function of these NAM genes in plant growth and stress responses. The 144 NAM genes identified in this study provide a basis for subsequent research and contribute to elucidating the intricate molecular mechanisms of NAM genes in Oncidium and potentially in other species.

Universal signatures of transposable element compartmentalization across eukaryotic genomes

The evolutionary mechanisms that drive the emergence of genome architecture remain poorly understood but can now be assessed with unprecedented power due to the massive accumulation of genome assemblies spanning phylogenetic diversity1,2. Transposable elements (TEs) are a rich source of large-effect mutations since they directly and indirectly drive genomic structural variation and changes in gene expression3. Here, we demonstrate universal patterns of TE compartmentalization across eukaryotic genomes spanning ~1.7 billion years of evolution, in which TEs colocalize with gene families under strong predicted selective pressure for dynamic evolution and involved in specific functions. For non-pathogenic species these genes represent families involved in defense, sensory perception and environmental interaction, whereas for pathogenic species, TE-compartmentalized genes are highly enriched for pathogenic functions. Many TE-compartmentalized gene families display signatures of positive selection at the molecular level. Furthermore, TE-compartmentalized genes exhibit an excess of high-frequency alleles for polymorphic TE insertions in fruit fly populations. We postulate that these patterns reflect selection for adaptive TE insertions as well as TE-associated structural variants. This process may drive the emergence of a shared TE-compartmentalized genome architecture across diverse eukaryotic lineages.

The chromosome-scale reference genome and transcriptome analysis of Solanum torvum provides insights into resistance to root-knot nematodes

Solanum torvum (Swartz) (2n = 24) is a wild Solanaceae plant with high economic value that is used as a rootstock in grafting for Solanaceae plants to improve the resistance to a soil-borne disease caused by root-knot nematodes (RKNs). However, the lack of a high-quality reference genome of S. torvum hinders research on the genetic basis for disease resistance and application in horticulture. Herein, we present a chromosome-level assembly of genomic sequences for S. torvum combining PacBio long reads (HiFi reads), Illumina short reads and Hi-C scaffolding technology. The assembled genome size is ~1.25 Gb with a contig N50 and scaffold N50 of 38.65 Mb and 103.02 Mb, respectively as well as a BUSCO estimate of 98%. GO enrichment and KEGG pathway analysis of the unique S. torvum genes, including NLR and ABC transporters, revealed that they were involved in disease resistance processes. RNA-seq data also confirmed that 48 NLR genes were highly expressed in roots and fibrous roots and that three homologous NLR genes (Sto0288260.1, Sto0201960.1 and Sto0265490.1) in S. torvum were significantly upregulated after RKN infection. Two ABC transporters, ABCB9 and ABCB11 were identified as the hub genes in response to RKN infection. The chromosome-scale reference genome of the S. torvum will provide insights into RKN resistance.

Exploring the Seasonal Dynamics and Molecular Mechanism of Wood Formation in Gymnosperm Trees

Forests, comprising 31% of the Earth's surface, play pivotal roles in regulating the carbon, water, and energy cycles. Despite being far less diverse than angiosperms, gymnosperms account for over 50% of the global woody biomass production. To sustain growth and development, gymnosperms have evolved the capacity to sense and respond to cyclical environmental signals, such as changes in photoperiod and seasonal temperature, which initiate growth (spring and summer) and dormancy (fall and winter). Cambium, the lateral meristem responsible for wood formation, is reactivated through a complex interplay among hormonal, genetic, and epigenetic factors. Temperature signals perceived in early spring induce the synthesis of several phytohormones, including auxins, cytokinins, and gibberellins, which in turn reactivate cambium cells. Additionally, microRNA-mediated genetic and epigenetic pathways modulate cambial function. As a result, the cambium becomes active during the summer, resulting in active secondary xylem (i.e., wood) production, and starts to become inactive in autumn. This review summarizes and discusses recent findings regarding the climatic, hormonal, genetic, and epigenetic regulation of wood formation in gymnosperm trees (i.e., conifers) in response to seasonal changes.

Genome-wide identification, characterization, and evolutionary analysis of NBS genes and their association with disease resistance in Musa spp

Banana is an important food crop that is susceptible to a wide range of pests and diseases that can reduce yield and quality. The primary objective of banana breeding programs is to increase disease resistance, which requires the identification of resistance (R) genes. Despite the fact that resistant sources have been identified in bananas, the genes, particularly the nucleotide-binding site (NBS) family, which play an important role in protecting plants against pathogens, have received little attention. As a result, this study included a thorough examination of the NBS disease resistance gene family's classification, phylogenetic analysis, genome organization, evolution, cis-elements, differential expression, regulation by microRNAs, and protein-protein interaction. A total of 116 and 43 putative NBS genes from M. acuminata and M. balbisiana, respectively, were identified and characterized, and were classified into seven sub-families. Structural analysis of NBS genes revealed the presence of signal peptides, their sub-cellular localization, molecular weight and pI. Eight commonly conserved motifs were found, and NBS genes were unevenly distributed across multiple chromosomes, with the majority of NBS genes being located in chr3 and chr1 of the A and B genomes, respectively. Tandem duplication occurrences have helped bananas' NBS genes spread throughout evolution. Transcriptome analysis of NBS genes revealed significant differences in expression between resistant and susceptible cultivars of fusarium wilt, eumusae leaf spot, root lesion nematode, and drought, implying that they can be used as candidate resistant genes. Ninety miRNAs were discovered to have targets in 104 NBS genes from the A genome, providing important insights into NBS gene expression regulation. Overall, this study offers a valuable genomic resource and understanding of the function and evolution of NBS genes in relation to rapidly evolving pathogens, as well as providing breeders with selection targets for fast-tracking breeding of banana varieties with more durable resistance to pathogens.

Genome-Wide Identification, Characterization, and Comparative Analysis of NLR Resistance Genes in Coffea spp

The largest family of disease resistance genes in plants are nucleotide-binding site leucine-rich repeat genes (NLRs). The products of these genes are responsible for recognizing avirulence proteins (Avr) of phytopathogens and triggering specific defense responses. Identifying NLRs in plant genomes with standard gene annotation software is challenging due to their multidomain nature, sequence diversity, and clustered genomic distribution. We present the results of a genome-wide scan and comparative analysis of NLR loci in three coffee species (Coffea canephora, Coffea eugenioides and their interspecific hybrid Coffea arabica). A total of 1311 non-redundant NLR loci were identified in C. arabica, 927 in C. canephora, and 1079 in C. eugenioides, of which 809, 562, and 695 are complete loci, respectively. The NLR-Annotator tool used in this study showed extremely high sensitivities and specificities (over 99%) and increased the detection of putative NLRs in the reference coffee genomes. The NLRs loci in coffee are distributed among all chromosomes and are organized mostly in clusters. The C. arabica genome presented a smaller number of NLR loci when compared to the sum of the parental genomes (C. canephora, and C. eugenioides). There are orthologous NLRs (orthogroups) shared between coffee, tomato, potato, and reference NLRs and those that are shared only among coffee species, which provides clues about the functionality and evolutionary history of these orthogroups. Phylogenetic analysis demonstrated orthologous NLRs shared between C. arabica and the parental genomes and those that were possibly lost. The NLR family members in coffee are subdivided into two main groups: TIR-NLR (TNL) and non-TNL. The non-TNLs seem to represent a repertoire of resistance genes that are important in coffee. These results will support functional studies and contribute to a more precise use of these genes for breeding disease-resistant coffee cultivars.

Epigenetics: a catalyst of plant immunity against pathogens

The plant immune system protects against pests and diseases. The recognition of stress-related molecular patterns triggers localised immune responses, which are often followed by longer-lasting systemic priming and/or up-regulation of defences. In some cases, this induced resistance (IR) can be transmitted to following generations. Such transgenerational IR is gradually reversed in the absence of stress at a rate that is proportional to the severity of disease experienced in previous generations. This review outlines the mechanisms by which epigenetic responses to pathogen infection shape the plant immune system across expanding time scales. We review the cis- and trans-acting mechanisms by which stress-inducible epigenetic changes at transposable elements (TEs) regulate genome-wide defence gene expression and draw particular attention to one regulatory model that is supported by recent evidence about the function of AGO1 and H2A.Z in transcriptional control of defence genes. Additionally, we explore how stress-induced mobilisation of epigenetically controlled TEs acts as a catalyst of Darwinian evolution by generating (epi)genetic diversity at environmentally responsive genes. This raises questions about the long-term evolutionary consequences of stress-induced diversification of the plant immune system in relation to the long-held dichotomy between Darwinian and Lamarckian evolution.

The whole-genome and expression profile analysis of WRKY and RGAs in Dactylis glomerata showed that DG6C02319.1 and DgWRKYs may cooperate in the immunity against rust

Orchardgrass (Dactylis glomerata) is one of the top four perennial forages worldwide and, despite its large economic advantages, often threatened by various environmental stresses. WRKY transcription factors (TFs) can regulate a variety of plant processes, widely participate in plant responses to biotic and abiotic stresses, and are one of the largest gene families in plants. WRKYs can usually bind W-box elements specifically. In this study, we identified a total of 93 DgWRKY genes and 281 RGAs, including 65, 169 and 47 nucleotide-binding site-leucine-rich repeats (NBS-LRRs), leucine-rich repeats receptor-like protein kinases (LRR-RLKs), and leucine-rich repeats receptor-like proteins (LRR-RLPs), respectively. Through analyzing the expression of DgWRKY genes in orchardgrass under different environmental stresses, it was found that many DgWRKY genes were differentially expressed under heat, drought, submergence, and rust stress. In particular, it was found that the greatest number of genes were differentially expressed under rust infection. Consistently, GO and KEGG enrichment analysis of all genes showed that 78 DgWRKY TFs were identified in the plant–pathogen interaction pathway, with 59 of them differentially expressed. Through cis-acting element prediction, 154 RGAs were found to contain W-box elements. Among them, DG6C02319.1 (a member of the LRR-RLK family) was identified as likely to interact with 14 DGWRKYs. Moreover, their expression levels in susceptible plants after rust inoculation were first up-regulated and then down-regulated, while those in the resistant plants were always up-regulated. In general, DgWRKYs responded to both biotic stress and abiotic stress. DgWRKYs and RGAs may synergistically respond to the response of orchardgrass to rust. This study provides meaningful insight into the molecular mechanisms of WRKY proteins in orchardgrass.

Comparative Analysis of Re-Annotated Genes Provides Insight into Evolutionary Divergence and Expressions of Aquaporin Family in Pepper

Aquaporins (AQPs) are known to have a vital role in water transport in all living organisms including agriculturally important crops, but a comprehensive genomic study of AQPs in pepper has not been implemented. Here, we updated previous gene annotations and generated a total of 259 AQP genes from five plants, including pepper. Phylogenetic and motif analyses revealed that a large proportion of pepper AQP genes belong to the specific subgroup of tonoplast intrinsic protein (TIP) subfamily, TIP4. Chromosomal localization and estimated duplication times illustrated that genes in TIP4 formed a tandem array on the short arm of chromosome 1, resulting from pepper-specific expansion after its divergence with Solanaceae species. Transcriptome analyses under various abiotic stress conditions revealed that transport-, photosystem-, and thylakoid-related genes were generally enriched in expression clusters containing AQP genes in pepper. These results provide valuable genomic resources and insight into the evolutionary mechanism that generate genomic diversity of the AQP gene family in pepper.

Building a cluster of NLR genes conferring resistance to pests and pathogens: the story of the Vat gene cluster in cucurbits

Most molecularly characterized plant resistance genes (R genes) belong to the nucleotide-binding-site-leucine-rich-repeat (NLR) receptor family and are prone to duplication and transposition with high sequence diversity. In this family, the Vat gene in melon is one of the few R genes known for conferring resistance to insect, i.e., Aphis gossypii, but it has been misassembled and/or mispredicted in the whole genomes of Cucurbits. We examined 14 genomic regions (about 400 kb) derived from long-read assemblies spanning Vat-related genes in Cucumis melo, Cucumis sativus, Citrullus lanatus, Benincasa hispida, Cucurbita argyrosperma, and Momordica charantia. We built the phylogeny of those genes. Investigating the paleohistory of the Vat gene cluster, we revealed a step by step process beginning from a common ancestry in cucurbits older than 50 my. We highlighted Vat exclusively in the Cucumis genera, which diverged about 20 my ago. We then focused on melon, evaluating a minimum duplication rate of Vat in 80 wild and cultivated melon lines using generalist primers; our results suggested that duplication started before melon domestication. The phylogeny of 44 Vat-CDS obtained from 21 melon lines revealed gain and loss of leucine-rich-repeat domains along diversification. Altogether, we revealed the high putative recognition scale offered in melon based on a combination of SNPs, number of leucine-rich-repeat domains within each homolog and number of homologs within each cluster that might jointly confer resistance to a large pest and pathogen spectrum. Based on our findings, we propose possible avenues for breeding programs.

Comparison and evolutionary analysis of Brassica nucleotide binding site leucine rich repeat (NLR) genes and importance for disease resistance breeding

The Brassica genus contains many agriculturally significant oilseed and vegetable crops, however the crop yield is threatened by a range of fungal and bacterial pathogens. Nucleotide Binding Site Leucine Rich Repeat (NLR) genes play important roles in plant innate immunity. The evolution of NLR genes is influenced by genomic processes and pathogen selection. At the whole genome level, whole genome duplications (WGDs) generate abundant gene copies, most of which are lost during genome fractionation. At sub-genomic levels, some retained copies undergo duplication forming clusters which facilitate rapid evolution through recombination. The number, distribution and genetic variations of the NLR genes vary among Brassica species and within populations suggesting differential selection pressure exerted by pathogen populations throughout the evolutionary history. A study of the evolution of disease resistance genes in agriculturally important plants such as Brassicas helps gain insights into their function and inform the identification of resistance genes for breeding of resistant lines.

Genome sequence analysis provides insights on genomic variation and late blight resistance genes in potato somatic hybrid (parents and progeny)

Wild Solanum species are the important resources for potato improvement. With the availability of potato genome and sequencing progress, knowledge about genomic resources is essential for novel genes discovery. Hence, the aim of this study was to decipher draft genome sequences of unique potato genotypes i.e. somatic hybrid P8 (J1), wild species S. pinnatisectum (J2), progeny MSH/14-112 (P8 × cv. Kufri Jyoti) (J3), and S. tuberosum dihaploid C-13 (J4). Draft genome sequencing using Illumina platform and reference-based assemblies with the potato genome yielded genome assembly size of 725.01 Mb (J1), 724.95 Mb (J2), 725.01 Mb (J3), and 809.59 Mb (J4). Further, 39,260 (J1), 25,711 (J2), 39,730 (J3) and 30,241 (J4) genes were identified and 17,411 genes were found common in the genotypes particularly late blight resistance genes (R3a, RGA2, RGA3, R1B-16, Rpi-blb2, Rpi and Rpi-vnt1). Gene ontology (GO) analysis showed that molecular function was predominant and signal transduction was major KEGG pathways. Further, gene enrichment analysis revealed dominance of metabolic process (GO: 0008152) in all the samples. Phylogeny analysis showed relatedness with potato and other plant species. Heterozygous single nucleotide polymorphism (SNP) was more than homozygous, and SNP in genic region was more than inter-genic region. Copy number variation (CNV) analysis indicated greater number of deletions than duplications. Sequence diversity and conserved motifs analysis revealed variation for late blight resistance genes. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis showed differential expression of late blight resistance genes. Our study provides insights on genome sequence, structural variation and late blight resistance genes in potato somatic hybrid (parents and progeny) for future research.

Evolution of NLR resistance genes with noncanonical N-terminal domains in wild tomato species

Nucleotide-binding and leucine-rich repeat immune receptors (NLRs) provide resistance against diverse pathogens. To create comparative NLR resources, we conducted resistance gene enrichment sequencing (RenSeq) with single-molecule real-time sequencing of PacBio for 18 accessions in Solanaceae, including 15 accessions of five wild tomato species. We investigated the evolution of a class of NLRs, CNLs with extended N-terminal sequences previously named Solanaceae Domain. Through comparative genomic analysis, we revealed that the extended CNLs (exCNLs) anciently emerged in the most recent common ancestor between Asterids and Amaranthaceae, far predating the Solanaceae family. In tomatoes, the exCNLs display exceptional modes of evolution in a clade-specific manner. In the clade G3, exCNLs have substantially elongated their N-termini through tandem duplications of exon segments. In the clade G1, exCNLs have evolved through recent proliferation and sequence diversification. In the clade G6, an ancestral exCNL has lost its N-terminal domains in the course of evolution. Our study provides high-quality NLR gene models for close relatives of domesticated tomatoes that can serve as a useful resource for breeding and molecular engineering for disease resistance. Our findings regarding the exCNLs offer unique backgrounds and insights for future functional studies of the NLRs.

Evolution of PHAS loci in the young spike of Allohexaploid wheat

BACKGROUND

PhasiRNAs (phased secondary siRNAs) play important regulatory roles in the development processes and biotic or abiotic stresses in plants. Some of phasiRNAs involve in the reproductive development in grasses, which include two categories, 21-nt (nucleotide) and 24-nt phasiRNAs. They are triggered by miR2118 and miR2275 respectively, in premeiotic and meiotic anthers of rice, maize and other grass species. Wheat (Triticum aestivum) with three closely related subgenomes (subA, subB and subD), is a model of allopolyploid in plants. Knowledge about the role of phasiRNAs in the inflorescence development of wheat is absent until now, and the evolution of PHAS loci in polyploid plants is also unavailable.

RESULTS

Using 261 small RNA expression datasets from various tissues, a batch of PHAS (phasiRNA precursors) loci were identified in the young spike of wheat, most of which were regulated by miR2118 and miR2275 in their target site regions. Dissection of PHAS and their trigger miRNAs among the diploid (AA and DD), tetraploid (AABB) and hexaploid (AABBDD) genomes of Triticum indicated that distribution of PHAS loci were dominant randomly in local chromosomes, while miR2118 was dominant only in the subB genome. The diversity of PHAS loci in the three subgenomes of wheat and their progenitor genomes (AA, DD and AABB) suggested that they originated or diverged at least before the occurrence of the tetraploid AABB genome. The positive correlation between the PHAS loci or the trigger miRNAs and the ploidy of genome indicated the expansion of genome was the major drive force for the increase of PHAS loci and their trigger miRNAs in Triticum. In addition, the expression profiles of the PHAS transcripts suggested they responded to abiotic stresses such as cold stress in wheat.

CONCLUSIONS

Altogether, non-coding phasiRNAs are conserved transcriptional regulators that display quick plasticity in Triticum genome. They may be involved in reproductive development and abiotic stress in wheat. It could be referred to molecular research on male reproductive development in Triticum.

Comparative Analysis of the YABBY Gene Family of Bienertia sinuspersici, a Single-Cell C4 Plant

The emergence and expression of the YABBY gene family (YGF) coincided with the evolution of leaves in seed plants, and was integral to the early evidence of lamina followed by reproductive development. YGF contains six subclasses, i.e., CRC, INO, FIL, YAB2, YAB3, and YAB5. This study aims to extract the genome sequences of the YGF in Bienertia sinuspersici, an important model plant for single-cell C₄ (SCC₄), non-Kranz photosynthesis. A comparative genomic analysis was undertaken with Vitis vinefera, Arabidopsis thaliana, Brassica rapa, and Chenopodium quinoa. Six copies of YGF were present in B. sinuspersici and A. thaliana with a single copy of each YGF subgroup. V. vinefera possessed seven copies of YGF with duplicates in FIL and YAB2 subgroups, but no YAB3. B. rapa and C. quinoa after whole genome duplication contained additional copies of YGF. The gene structure and conserved motifs were analyzed among the YGF. In addition, the relative quantification of YGF was analyzed in the leaves, reproductive developmental stages such as the bud, and the pre-anthesis and anthesis stages in B. sinuspersici, A. thaliana, and B. rapa. CRC and INO possessed conserved floral-specific expression. Temporal and perpetual changes in the expression of YGF orthologs were observed in the leaves and reproductive developmental stages. The results of this study provide an overview of YGF evolution, copy number, and its differential expression in B. sinuspersici. Further studies are required to shed light on the roles of YABBY genes in the evolution of SCC₄ plants and their distinct physiologies.

Surviving in a Hostile World: Plant Strategies to Resist Pests and Diseases

As primary producers, plants are under constant pressure to defend themselves against potentially deadly pathogens and herbivores. In this review, we describe short- and long-term strategies that enable plants to cope with these stresses. Apart from internal immunological strategies that involve physiological and (epi)genetic modifications at the cellular level, plants also employ external strategies that rely on recruitment of beneficial organisms. We discuss these strategies along a gradient of increasing timescales, ranging from rapid immune responses that are initiated within seconds to (epi)genetic adaptations that occur over multiple plant generations. We cover the latest insights into the mechanistic and evolutionary underpinnings of these strategies and present explanatory models. Finally, we discuss how knowledge from short-lived model species can be translated to economically and ecologically important perennials to exploit adaptive plant strategies and mitigate future impacts of pests and diseases in an increasingly interconnected and changing world.

Evolution of Disease Defense Genes and Their Regulators in Plants

Biotic stresses do damage to the growth and development of plants, and yield losses for some crops. Confronted with microbial infections, plants have evolved multiple defense mechanisms, which play important roles in the never-ending molecular arms race of plant–pathogen interactions. The complicated defense systems include pathogen-associated molecular patterns (PAMP) triggered immunity (PTI), effector triggered immunity (ETI), and the exosome-mediated cross-kingdom RNA interference (CKRI) system. Furthermore, plants have evolved a classical regulation system mediated by miRNAs to regulate these defense genes. Most of the genes/small RNAs or their regulators that involve in the defense pathways can have very rapid evolutionary rates in the longitudinal and horizontal co-evolution with pathogens. According to these internal defense mechanisms, some strategies such as molecular switch for the disease resistance genes, host-induced gene silencing (HIGS), and the new generation of RNA-based fungicides, have been developed to control multiple plant diseases. These broadly applicable new strategies by transgene or spraying ds/sRNA may lead to reduced application of pesticides and improved crop yield.

Genome-wide identification, characterization, and evolutionary analysis of NBS-encoding resistance genes in barley

In this study, a systematic analysis of Nucleotide-Binding Site (NBS) disease resistance (R) gene family in the barley, Hordeum vulgare L. cv. Bowman, genome was performed. Using multiple computational analyses, we could identify 96 regular NBS-encoding genes and characterize them on the bases of structural diversity, conserved protein signatures, genomic distribution, gene duplications, differential expression, selection pressure, codon usage, regulation by microRNAs and phylogenetic relationships. Depending on the presence or absence of CC and LRR domains; the identified NBS genes were assigned to four distinct groups; NBS-LRR (53.1%), CC-NBS-LRR (14.6%), NBS (26%), and CC-NBS (6.3%). NBS-associated domain analysis revealed the presence of signal peptides, zinc fingers, diverse kinases, and other structural features. Eighty-five of the identified NBS-encoding genes were mapped onto the seven barley chromosomes, revealing that 50% of them were located on chromosomes 7H, 2H, and 3H, with a tendency of NBS genes to be clustered in the distal telomeric regions of the barley chromosomes. Nine gene clusters, representing 22.35% of total mapped barley NBS-encoding genes, were found, suggesting that tandem duplication stands for an important mechanism in the expansion of this gene family in barley. Phylogenetic analysis determined 31 HvNBS orthologs from rice and Brachypodium. 87 out of 96 HvNBSs were supported by expression evidence, exhibiting various and quantitatively uneven expression patterns across distinct tissues, organs, and development stages. Fourteen potential miRNA-R gene target pairs were further identified, providing insight into the regulation of NBS genes expression. These findings offer candidate target genes to engineer disease-resistant barley genotypes, and promote our understanding of the evolution of NBS-encoding genes in Poaceae crops.

Nucleotide-binding resistance gene signatures in sugar beet, insights from a new reference genome

Nucleotide-binding (NB-ARC), leucine-rich-repeat genes (NLRs) account for 60.8% of resistance (R) genes molecularly characterized from plants. NLRs exist as large gene families prone to tandem duplication and transposition, with high sequence diversity among crops and their wild relatives. This diversity can be a source of new disease resistance, but difficulty in distinguishing specific sequences from homologous gene family members hinders characterization of resistance for improving crop varieties. Current genome sequencing and assembly technologies, especially those using long-read sequencing, are improving resolution of repeat-rich genomic regions and clarifying locations of duplicated genes, such as NLRs. Using the conserved NB-ARC domain as a model, 231 tentative NB-ARC loci were identified in a highly contiguous genome assembly of sugar beet, revealing diverged and truncated NB-ARC signatures as well as full-length sequences. The NB-ARC-associated proteins contained NLR resistance gene domains, including TIR, CC and LRR, as well as other integrated domains. Phylogenetic relationships of partial and complete domains were determined, and patterns of physical clustering in the genome were evaluated. Comparison of sugar beet NB-ARC domains to validated R-genes from monocots and eudicots suggested extensive Beta vulgaris-specific subfamily expansions. The NLR landscape in the rhizomania resistance conferring Rz region of Chromosome 3 was characterized, identifying 26 NLR-like sequences spanning 20 MB. This work presents the first detailed view of NLR family composition in a member of the Caryophyllales, builds a foundation for additional disease resistance work in B. vulgaris, and demonstrates an additional nucleic-acid-based method for NLR prediction in non-model plant species.

Copy number variation and disease resistance in plants

Plant genome diversity varies from single nucleotide polymorphisms to large-scale deletions, insertions, duplications, or re-arrangements. These re-arrangements of sequences resulting from duplication, gains or losses of DNA segments are termed copy number variations (CNVs). During the last decade, numerous studies have emphasized the importance of CNVs as a factor affecting human phenotype; in particular, CNVs have been associated with risks for several severe diseases. In plants, the exploration of the extent and role of CNVs in resistance against pathogens and pests is just beginning. Since CNVs are likely to be associated with disease resistance in plants, an understanding of the distribution of CNVs could assist in the identification of novel plant disease-resistance genes. In this paper, we review existing information about CNVs; their importance, role and function, as well as their association with disease resistance in plants.

Conservation and Dispersion of Genes Conferring Resistance to Tomato Begomoviruses between Tomato and Pepper Genomes

In the present climate change scenario, controlling plant disease through exploitation of host plant resistance could contribute toward the sustainable crop production and global food security. In this respect, the identification of new sources of resistance and utilization of genetic diversity within the species may help in the generation of cultivars with improved disease resistance. Begomoviruses namely, Tomato yellow leaf curl virus (TYLCV) and Chilli leaf curl virus (ChLCV) are known to cause major yield losses in several economically important crop plants of the family Solanaceae. Though co-occurrence, association and synergistic interactions among these viruses in the host plants is reported, whether orthologous genetic loci in related host plants could be responsible for conferring resistance to these viruses has not been investigated yet. Several loci including Ty1, Ty2, Ty3, Ty4, and ty5 have been reported to confer resistance to leaf curl viruses in tomato. Here, we examined the pepper orthologous markers, corresponding to these QTL regions, for polymorphism between ChLCV susceptible and resistant genotypes of pepper. Further, to examine if the polymorphic markers are segregating with the disease resistance, Bulk Segregant Analysis (BSA) was performed on F₂ population derived from crosses between resistant and susceptible lines. However, none of the markers showed polymorphism in BSA suggesting that the tested markers are not linked to genes/QTLs responsible for conferring resistance to ChLCV in the selected genotypes. In silico analysis was performed to study the synteny and collinearity of genes located within these QTL regions in tomato and pepper genomes, which revealed that more than 60% genes located in Ty2 and Ty4, 13.71% genes in Ty1, 23.07% in Ty3, and 44.77% genes located within ty5 QTL region in tomato are conserved in pepper genome. However, despite such a high conservation in gene content, the linkage relationship in these regions seems to be greatly affected by gross rearrangements in both the species.

Transposable Elements Direct The Coevolution between Plants and Microbes

Transposable elements are powerful drivers of genome evolution in many eukaryotes. Although they are mostly considered as 'selfish' genetic elements, increasing evidence suggests that they contribute to genetic variability; particularly under stress conditions. Over the past few years, the role of transposable elements during host-microbe interactions has been recognised. It has been proposed that many pathogenic microbes have evolved a 'two-speed' genome with regions that show increased variability and that are enriched in transposable elements and pathogenicity-related genes. Plants similarly display structured genomes with transposable-element-rich regions that mediate accelerated evolution. Immune receptor genes typically reside in such regions. Various mechanisms have recently been identified through which transposable elements contribute to the coevolution between plants and their associated microbes.

Investigation of orthologous pathogen recognition gene-rich regions in solanaceous species

Pathogen receptor proteins such as receptor-like protein (RLP), receptor-like kinase (RLK), and nucleotide-binding leucine-rich repeat (NLR) play a leading role in plant immunity activation. The genome architecture of such genes has been extensively investigated in several plant species. However, we still know little about their elaborate reorganization that arose during the plant speciation process. Using recently released pepper and eggplant genome sequences, we were able to identify 1097 pathogen recognition genes (PRGs) in the cultivated pepper Zunla-1 and 775 in the eggplant line Nakate-Shinkuro. The retrieved genes were analysed for their tendency to cluster, using different methods to infer the means of grouping. Orthologous relationships among clustering loci were found, and interesting reshuffling within given loci was observed for each analysed species. The information obtained was integrated into a comparative map to highlight the evolutionary dynamics in which the PRG loci were involved. Diversification of 14 selected PRG-rich regions was also explored using a DNA target-enrichment approach. A large number of gene variants were found as well as rearrangements of sequences encoding single protein domain and changes in chromosome gene order among species. Gene duplication and transposition activity have clearly influenced plant genome R-gene architecture and diversification. Our findings contribute to addressing several biological questions concerning the parallel evolution that occurred between genomes of the family Solanaceae. Moreover, the integration of different methods proved a powerful approach to reconstruct the evolutionary history in plant families and to transfer important biology findings among plant genomes.

NLR network mediates immunity to diverse plant pathogens

Both plants and animals rely on nucleotide-binding domain and leucine-rich repeat-containing (NLR) proteins to respond to invading pathogens and activate immune responses. An emerging concept of NLR function is that "sensor" NLR proteins are paired with "helper" NLRs to mediate immune signaling. However, our fundamental knowledge of sensor/helper NLRs in plants remains limited. In this study, we discovered a complex NLR immune network in which helper NLRs in the NRC (NLR required for cell death) family are functionally redundant but display distinct specificities toward different sensor NLRs that confer immunity to oomycetes, bacteria, viruses, nematodes, and insects. The helper NLR NRC4 is required for the function of several sensor NLRs, including Rpi-blb2, Mi-1.2, and R1, whereas NRC2 and NRC3 are required for the function of the sensor NLR Prf. Interestingly, NRC2, NRC3, and NRC4 redundantly contribute to the immunity mediated by other sensor NLRs, including Rx, Bs2, R8, and Sw5. NRC family and NRC-dependent NLRs are phylogenetically related and cluster into a well-supported superclade. Using extensive phylogenetic analysis, we discovered that the NRC superclade probably emerged over 100 Mya from an NLR pair that diversified to constitute up to one-half of the NLRs of asterids. These findings reveal a complex genetic network of NLRs and point to a link between evolutionary history and the mechanism of immune signaling. We propose that this NLR network increases the robustness of immune signaling to counteract rapidly evolving plant pathogens.

Evolutionary plasticity of restorer-of-fertility-like proteins in rice

Hybrid seed production in rice relies on cytoplasmic male sterility (CMS) induced by specific mitochondrial proteins, whose deleterious effects are suppressed by nuclear Restorer of Fertility (RF) genes. The majority of RF proteins belong to a specific clade of the RNA-binding pentatricopeptide repeat protein family. We have characterised 'restorer-of-fertility-like' (RFL) sequences from 13 Oryza genomes and the Brachypodium distachyon genome. The majority of the RFL sequences are found in genomic clusters located at two or three chromosomal loci with only a minor proportion being present as isolated genes. The RFL genomic cluster located on Oryza chromosome 10, the location of almost all known active rice RF genes, shows extreme variation in structure and gene content between species. We show evidence for homologous recombination events as an efficient mechanism for generating the huge repertoire of RNA sequence recognition motifs within RFL proteins and a major driver of RFL sequence evolution. The RFL sequences identified here will improve our understanding of the molecular basis of CMS and fertility restoration in plants and will accelerate the development of new breeding strategies.

Genetic mapping of Pinus flexilis major gene (Cr4) for resistance to white pine blister rust using transcriptome-based SNP genotyping

Background

Linkage of DNA markers with phenotypic traits provides essential information to dissect clustered genes with potential phenotypic contributions in a target genome region. Pinus flexilis E. James (limber pine) is a keystone five-needle pine species in mountain-top ecosystems of North America. White pine blister rust (WPBR), caused by a non-native fungal pathogen Cronartium ribicola (J.C. Fisch.), has resulted in mortality in this conifer species and is still spreading through the distribution. The objective of this research was to develop P. flexilis transcriptome-wide single nucleotide polymorphism (SNP) markers using RNA-seq analysis for genetic mapping of the major gene (Cr4) that confers complete resistance to C. ribicola.

Results

Needle tissues of one resistant and two susceptible seedling families were subjected to RNA-seq analysis. In silico SNP markers were uncovered by mapping the RNA-seq reads back to the de novo assembled transcriptomes. A total of 110,573 in silico SNPs and 2,870 indels were identified with an average of 3.7 SNPs per Kb. These SNPs were distributed in 17,041 unigenes. Of these polymorphic P. flexilis unigenes, 6,584 were highly conserved as compared to the genome sequence of P. taeda L (loblolly pine). High-throughput genotyping arrays were designed and were used to search for Cr4-linked genic SNPs in megagametophyte populations of four maternal trees by haploid-segregation analysis. A total of 32 SNP markers in 25 genes were localized on the Cr4 linkage group (LG). Syntenic relationships of this Cr4-LG map with the model conifer species P. taeda anchored Cr4 on Pinus consensus LG8, indicating that R genes against C. ribicola have evolved independently in different five-needle pines. Functional genes close to Cr4 were annotated and their potential roles in Cr4-mediated resistance were further discussed.

Conclusions

We demonstrated a very effective, low-cost strategy for developing a SNP genetic map of a phenotypic trait of interest. SNP discovery through transcriptome comparison was integrated with high-throughput genotyping of a small set of in silico SNPs. This strategy may be applied to mapping any trait in non-model plant species that have complex genomes. Whole transcriptome sequencing provides a powerful tool for SNP discovery in conifers and other species with complex genomes, for which sequencing and annotation of complex genomes is still challenging. The genic SNP map for the consensus Cr4-LG may help future molecular breeding efforts by enabling both Cr4 positional characterization and selection of this gene against WPBR.

Electronic supplementary material

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

The Diversification of Plant NBS-LRR Defense Genes Directs the Evolution of MicroRNAs That Target Them

High expression of plant nucleotide binding site leucine-rich repeat (NBS-LRR) defense genes is often lethal to plant cells, a phenotype perhaps associated with fitness costs. Plants implement several mechanisms to control the transcript level of NBS-LRR defense genes. As negative transcriptional regulators, diverse miRNAs target NBS-LRRs in eudicots and gymnosperms. To understand the evolutionary benefits of this miRNA-NBS-LRR regulatory system, we investigated the NBS-LRRs of 70 land plants, coupling this analysis with extensive small RNA data. A tight association between the diversity of NBS-LRRs and miRNAs was found. The miRNAs typically target highly duplicated NBS-LRRs In comparison, families of heterogeneous NBS-LRRs were rarely targeted by miRNAs in Poaceae and Brassicaceae genomes. We observed that duplicated NBS-LRRs from different gene families periodically gave birth to new miRNAs. Most of these newly emerged miRNAs target the same conserved, encoded protein motif of NBS-LRRs, consistent with a model of convergent evolution for these miRNAs. By assessing the interactions between miRNAs and NBS-LRRs, we found nucleotide diversity in the wobble position of the codons in the target site drives the diversification of miRNAs. Taken together, we propose a co-evolutionary model of plant NBS-LRRs and miRNAs hypothesizing how plants balance the benefits and costs of NBS-LRR defense genes.

Evolution of Gene Duplication in Plants

Ancient duplication events and a high rate of retention of extant pairs of duplicate genes have contributed to an abundance of duplicate genes in plant genomes. These duplicates have contributed to the evolution of novel functions, such as the production of floral structures, induction of disease resistance, and adaptation to stress. Additionally, recent whole-genome duplications that have occurred in the lineages of several domesticated crop species, including wheat (Triticum aestivum), cotton (Gossypium hirsutum), and soybean (Glycine max), have contributed to important agronomic traits, such as grain quality, fruit shape, and flowering time. Therefore, understanding the mechanisms and impacts of gene duplication will be important to future studies of plants in general and of agronomically important crops in particular. In this review, we survey the current knowledge about gene duplication, including gene duplication mechanisms, the potential fates of duplicate genes, models explaining duplicate gene retention, the properties that distinguish duplicate from singleton genes, and the evolutionary impact of gene duplication.

Molecular phylogeny and dynamic evolution of disease resistance genes in the legume family

Background

Legumes are the second-most important crop family in agriculture for its economic and nutritional values. Disease resistance (R-) genes play an important role in responding to pathogen infections in plants. To further increase the yield of legume crops, we need a comprehensive understanding of the evolution of R-genes in the legume family.

Results

In this study, we developed a robust pipeline and identified a total of 4,217 R-genes in the genomes of seven sequenced legume species. A dramatic diversity of R-genes with structural variances indicated a rapid birth-and-death rate during the R-gene evolution in legumes. The number of R-genes transiently expanded and then quickly contracted after whole-genome duplications, which meant that R-genes were sensitive to subsequent diploidization. R proteins with the Coiled-coil (CC) domain are more conserved than others in legumes. Meanwhile, other types of legume R proteins with only one or two typical domains were subjected to higher rates of loss during evolution. Although R-genes evolved quickly in legumes, they tended to undergo purifying selection instead of positive selection during evolution. In addition, domestication events in some legume species preferentially selected for the genes directly involved in the plant-pathogen interaction pathway while suppressing those R-genes with low occurrence rates.

Conclusions

Our results provide insights into the dynamic evolution of R-genes in the legume family, which will be valuable for facilitating genetic improvements in the disease resistance of legume cultivars.

Electronic supplementary material

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

Genome-Wide Analysis of NBS-LRR Genes in Sorghum Genome Revealed Several Events Contributing to NBS-LRR Gene Evolution in Grass Species

The nucleotide-binding site (NBS)–leucine-rich repeat (LRR) gene family is crucially important for offering resistance to pathogens. To explore evolutionary conservation and variability of NBS-LRR genes across grass species, we identified 88, 107, 24, and 44 full-length NBS-LRR genes in sorghum, rice, maize, and Brachypodium, respectively. A comprehensive analysis was performed on classification, genome organization, evolution, expression, and regulation of these NBS-LRR genes using sorghum as a representative of grass species. In general, the full-length NBS-LRR genes are highly clustered and duplicated in sorghum genome mainly due to local duplications. NBS-LRR genes have basal expression levels and are highly potentially targeted by miRNA. The number of NBS-LRR genes in the four grass species is positively correlated with the gene clustering rate. The results provided a valuable genomic resource and insights for functional and evolutionary studies of NBS-LRR genes in grass species.

Non-host Resistance Induced by the Xanthomonas Effector XopQ Is Widespread within the Genus Nicotiana and Functionally Depends on EDS1

Most Gram-negative plant pathogenic bacteria translocate effector proteins (T3Es) directly into plant cells via a conserved type III secretion system, which is essential for pathogenicity in susceptible plants. In resistant plants, recognition of some T3Es is mediated by corresponding resistance (R) genes or R proteins and induces effector triggered immunity (ETI) that often results in programmed cell death reactions. The identification of R genes and understanding their evolution/distribution bears great potential for the generation of resistant crop plants. We focus on T3Es from Xanthomonas campestris pv. vesicatoria (Xcv), the causal agent of bacterial spot disease on pepper and tomato plants. Here, 86 Solanaceae lines mainly of the genus Nicotiana were screened for phenotypical reactions after Agrobacterium tumefaciens-mediated transient expression of 21 different Xcv effectors to (i) identify new plant lines for T3E characterization, (ii) analyze conservation/evolution of putative R genes and (iii) identify promising plant lines as repertoire for R gene isolation. The effectors provoked different reactions on closely related plant lines indicative of a high variability and evolution rate of potential R genes. In some cases, putative R genes were conserved within a plant species but not within superordinate phylogenetical units. Interestingly, the effector XopQ was recognized by several Nicotiana spp. lines, and Xcv infection assays revealed that XopQ is a host range determinant in many Nicotiana species. Non-host resistance against Xcv and XopQ recognition in N. benthamiana required EDS1, strongly suggesting the presence of a TIR domain-containing XopQ-specific R protein in these plant lines. XopQ is a conserved effector among most xanthomonads, pointing out the XopQ-recognizing R_xopQ as candidate for targeted crop improvement.

Genome-Wide Architecture of Disease Resistance Genes in Lettuce

Genome-wide motif searches identified 1134 genes in the lettuce reference genome of cv. Salinas that are potentially involved in pathogen recognition, of which 385 were predicted to encode nucleotide binding-leucine rich repeat receptor (NLR) proteins. Using a maximum-likelihood approach, we grouped the NLRs into 25 multigene families and 17 singletons. Forty-one percent of these NLR-encoding genes belong to three families, the largest being RGC16 with 62 genes in cv. Salinas. The majority of NLR-encoding genes are located in five major resistance clusters (MRCs) on chromosomes 1, 2, 3, 4, and 8 and cosegregate with multiple disease resistance phenotypes. Most MRCs contain primarily members of a single NLR gene family but a few are more complex. MRC2 spans 73 Mb and contains 61 NLRs of six different gene families that cosegregate with nine disease resistance phenotypes. MRC3, which is 25 Mb, contains 22 RGC21 genes and colocates with Dm13. A library of 33 transgenic RNA interference tester stocks was generated for functional analysis of NLR-encoding genes that cosegregated with disease resistance phenotypes in each of the MRCs. Members of four NLR-encoding families, RGC1, RGC2, RGC21, and RGC12 were shown to be required for 16 disease resistance phenotypes in lettuce. The general composition of MRCs is conserved across different genotypes; however, the specific repertoire of NLR-encoding genes varied particularly of the rapidly evolving Type I genes. These tester stocks are valuable resources for future analyses of additional resistance phenotypes.

Small RNA and degradome sequencing reveals microRNAs and their targets involved in tomato pedicel abscission

We constructed small RNA and degradome sequencing libraries to identify miRNAs and targets involved in tomato pedicel abscission, and confirmed their roles via quantitative real-time PCR. MicroRNAs (miRNAs) are endogenous small RNAs which play crucial negatively regulatory roles at both the transcriptional and post-transcriptional levels in plants; however, limited knowledge is available on the expression profiles of miRNAs and their target genes during tomato pedicel abscission. Taking advantage of small RNA (sRNA) and degradome sequencing technology, a total of 56 known and 11 novel candidate miRNAs targeting 223 mRNA genes were confirmed during pedicel abscission. Gene ontology annotation and KEGG pathway analysis showed that these target genes were significantly enriched in intracellular, membrane-bounded organelle-related biological processes as well as in metabolic, plant-pathogen interaction and hormone signaling pathways. We screened 17 miRNA/target pairs for further analysis and performed quantitative real-time PCR to identify the roles. Cluster analysis of selected miRNAs revealed that the expression profiles of miRNAs varied in different stages of abscission and could be impacted by ethylene treatment. In the present study, the correlations between miRNAs and targets suggested a complex regulatory network of miRNA-mediated target interaction during pedicel abscission. Additionally, the expression profiles of miRNAs and their targets changed by ethylene might be a considerable reason why ethylene promotes pedicel abscission. Our study provides new insights into the expression and regulatory profiles of miRNAs during tomato pedicel abscission.

Resistance gene analogs involved in tolerant cassava--geminivirus interaction that shows a recovery phenotype

The current literature describes recovery from virus-induced symptoms as a RNA silencing defense, but immunity-related genes, including the structurally specific resistance gene analogs (RGAs) that may play a key role in tolerance and recovery is not yet reported. In this study, the transcriptome data of tolerant cassava TME3 (which exhibits a recovery phenotype) and susceptible cassava T200 infected with South African cassava mosaic virus were explored for RGAs. Putative resistance protein analogs (RPAs) with amide-like indole-3-acetic acid-Ile-Leu-Arg (IAA-ILR) and leucine-rich repeat (LRR)-kinase conserved domains were unique to TME3. Common responsive RPAs in TME3 and T200 were the dirigent-like protein, coil-coil nucleotide-binding site (NBS) and toll-interleukin-resistance, disease resistance zinc finger chromosome condensation-like protein (DZC), and NBS-apoptosis repressor with caspase recruitment (ARC)-LRR domains. Mutations in RPAs in the MHD motif of the NBS-ARC2 subdomain associated with the recovery phase in TME3 were observed. Additionally, a cohort of 25 RGAs mined solely during the recovery process in TME3 was identified. Phylogenetic and expression analyses support that diverse RGAs are differentially expressed during tolerance and recovery. This study reveals that in cassava, a perennial crop, RGAs participate in tolerance and differentially accumulate during recovery as a complementary defense mechanism to natural occurring RNA silencing to impair viral replication.

Stress responsive miRNAs and isomiRs in cereals

Abiotic and biotic stress conditions are vital determinants in the production of cereals, the major caloric source in human nutrition. Small RNAs, miRNAs and isomiRs are central to post-transcriptional regulation of gene expression in a variety of cellular processes including development and stress responses. Several miRNAs have been identified using new technologies and have roles in stress responses in plants, including cereals. The overall knowledge about the cereal miRNA repertoire, as well as an understanding of complex miRNA mediated mechanisms of target regulation in response to stress conditions, is far from complete. Ongoing efforts that add to our understanding of complex miRNA machinery have implications in plant response to stress conditions. Additionally, sequence variants of miRNAs (isomiRNAs or isomiRs), regulation of their expression through dissection of upstream regulatory elements, the role of Processing-bodies (P-bodies) in miRNA exerted gene regulation and yet unveiled organellar plant miRNAs are newly emerging topics, which will contribute to the elucidation of the miRNA machinery and its role in cereal tolerance against abiotic and biotic stresses.

The Solanum commersonii Genome Sequence Provides Insights into Adaptation to Stress Conditions and Genome Evolution of Wild Potato Relatives

Here, we report the draft genome sequence of Solanum commersonii, which consists of ∼830 megabases with an N50 of 44,303 bp anchored to 12 chromosomes, using the potato (Solanum tuberosum) genome sequence as a reference. Compared with potato, S. commersonii shows a striking reduction in heterozygosity (1.5% versus 53 to 59%), and differences in genome sizes were mainly due to variations in intergenic sequence length. Gene annotation by ab initio prediction supported by RNA-seq data produced a catalog of 1703 predicted microRNAs, 18,882 long noncoding RNAs of which 20% are shown to target cold-responsive genes, and 39,290 protein-coding genes with a significant repertoire of nonredundant nucleotide binding site-encoding genes and 126 cold-related genes that are lacking in S. tuberosum. Phylogenetic analyses indicate that domesticated potato and S. commersonii lineages diverged ∼2.3 million years ago. Three duplication periods corresponding to genome enrichment for particular gene families related to response to salt stress, water transport, growth, and defense response were discovered. The draft genome sequence of S. commersonii substantially increases our understanding of the domesticated germplasm, facilitating translation of acquired knowledge into advances in crop stability in light of global climate and environmental changes.

A genetic map of cassava (Manihot esculenta Crantz) with integrated physical mapping of immunity-related genes

BACKGROUND

Cassava, Manihot esculenta Crantz, is one of the most important crops world-wide representing the staple security for more than one billion of people. The development of dense genetic and physical maps, as the basis for implementing genetic and molecular approaches to accelerate the rate of genetic gains in breeding program represents a significant challenge. A reference genome sequence for cassava has been made recently available and community efforts are underway for improving its quality. Cassava is threatened by several pathogens, but the mechanisms of defense are far from being understood. Besides, there has been a lack of information about the number of genes related to immunity as well as their distribution and genomic organization in the cassava genome.

RESULTS

A high dense genetic map of cassava containing 2,141 SNPs has been constructed. Eighteen linkage groups were resolved with an overall size of 2,571 cM and an average distance of 1.26 cM between markers. More than half of mapped SNPs (57.4%) are located in coding sequences. Physical mapping of scaffolds of cassava whole genome sequence draft using the mapped markers as anchors resulted in the orientation of 687 scaffolds covering 45.6% of the genome. One hundred eighty nine new scaffolds are anchored to the genetic cassava map leading to an extension of the present cassava physical map with 30.7 Mb. Comparative analysis using anchor markers showed strong co-linearity to previously reported cassava genetic and physical maps. In silico based searching for conserved domains allowed the annotation of a repertory of 1,061 cassava genes coding for immunity-related proteins (IRPs). Based on physical map of the corresponding sequencing scaffolds, unambiguous genetic localization was possible for 569 IRPs.

CONCLUSIONS

This is the first study reported so far of an integrated high density genetic map using SNPs with integrated genetic and physical localization of newly annotated immunity related genes in cassava. These data build a solid basis for future studies to map and associate markers with single loci or quantitative trait loci for agronomical important traits. The enrichment of the physical map with novel scaffolds is in line with the efforts of the cassava genome sequencing consortium.

The plasticity of NBS resistance genes in sorghum is driven by multiple evolutionary processes

BACKGROUND

Increased disease resistance is a key target of cereal breeding programs, with disease outbreaks continuing to threaten global food production, particularly in Africa. Of the disease resistance gene families, the nucleotide-binding site plus leucine-rich repeat (NBS-LRR) family is the most prevalent and ancient and is also one of the largest gene families known in plants. The sequence diversity in NBS-encoding genes was explored in sorghum, a critical food staple in Africa, with comparisons to rice and maize and with comparisons to fungal pathogen resistance QTL.

RESULTS

In sorghum, NBS-encoding genes had significantly higher diversity in comparison to non NBS-encoding genes and were significantly enriched in regions of the genome under purifying and balancing selection, both through domestication and improvement. Ancestral genes, pre-dating species divergence, were more abundant in regions with signatures of selection than in regions not under selection. Sorghum NBS-encoding genes were also significantly enriched in the regions of the genome containing fungal pathogen disease resistance QTL; with the diversity of the NBS-encoding genes influenced by the type of co-locating biotic stress resistance QTL.

CONCLUSIONS

NBS-encoding genes are under strong selection pressure in sorghum, through the contrasting evolutionary processes of purifying and balancing selection. Such contrasting evolutionary processes have impacted ancestral genes more than species-specific genes. Fungal disease resistance hot-spots in the genome, with resistance against multiple pathogens, provides further insight into the mechanisms that cereals use in the "arms race" with rapidly evolving pathogens in addition to providing plant breeders with selection targets for fast-tracking the development of high performing varieties with more durable pathogen resistance.