The maintenance of extreme amino acid diversity at the disease resistance gene, RPP13, in Arabidopsis thaliana

Identification of candidate single-nucleotide polymorphisms (SNPs) and genes associated with sugarcane leaf scald disease

Leaf scald, caused by Xanthomonas albilineans, is a severe disease affecting sugarcane worldwide. One of the most practical ways to control it is by developing resistant sugarcane cultivars. It is essential to identify genes associated with the response to leaf scald. A panel of 170 sugarcane genotypes was evaluated for resistance to leaf scald in field conditions for 2 years, followed by a 1-year greenhouse experiment. The phenotypic evaluation data showed a wide continuous distribution, with heritability values ranging from 0.58 to 0.84. Thirteen single nucleotide polymorphisms (SNPs) were identified, significantly associated with leaf scald resistance. Among these, eight were stable across multiple environments and association models. The candidate genes identified and validated based on RNA-seq and qRT-PCR included two genes that encode NB-ARC leucine-rich repeat (LRR)-containing domain disease-resistance protein. These findings provide a basis for developing marker-assisted selection strategies in sugarcane breeding programs.

Single amino acid change alters specificity of the multi-allelic wheat stem rust resistance locus SR9

Most rust resistance genes thus far isolated from wheat have a very limited number of functional alleles. Here, we report the isolation of most of the alleles at wheat stem rust resistance gene locus SR9. The seven previously reported resistance alleles (Sr9a, Sr9b, Sr9d, Sr9e, Sr9f, Sr9g, and Sr9h) are characterised using a synergistic strategy. Loss-of-function mutants and/or transgenic complementation are used to confirm Sr9b, two haplotypes of Sr9e (Sr9e_h1 and Sr9e_h2), Sr9g, and Sr9h. Each allele encodes a highly related nucleotide-binding site leucine-rich repeat (NB-LRR) type immune receptor, containing an unusual long LRR domain, that confers resistance to a unique spectrum of isolates of the wheat stem rust pathogen. The only SR9 protein effective against stem rust pathogen race TTKSK (Ug99), SR9H, differs from SR9B by a single amino acid. SR9B and SR9G resistance proteins are also distinguished by only a single amino acid. The SR9 allelic series found in the B subgenome are orthologs of wheat stem rust resistance gene Sr21 located in the A subgenome with around 85% identity in protein sequences. Together, our results show that functional diversification of allelic variants at the SR9 locus involves single and multiple amino acid changes that recognize isolates of wheat stem rust.

Integrating omics reveals that miRNA-guided genetic regulation on plant hormone level and defense response pathways shape resistance to Cladosporium fulvum in the tomato Cf-10-gene-carrying line

Invasion of C. fulvum causes the most serious diseases affecting the reproduction of tomatoes. Cf-10-gene-carrying line showed remarkable resistance to Cladosporium fulvum. To exploit its defense response mechanism, we performed a multiple-omics profiling of Cf-10-gene-carrying line and a susceptible line without carrying any resistance genes at non-inoculation and 3 days post-inoculation (dpi) of C. fulvum. We detected 54 differentially expressed miRNAs (DE-miRNAs) between the non-inoculation and 3 dpi in the Cf-10-gene-carrying line, which potentially regulated plant-pathogen interaction pathways and hormone signaling pathways. We also revealed 3,016 differentially expressed genes (DEGs) between the non-inoculated and 3 dpi in the Cf-10-gene-carrying line whose functions enriched in pathways that were potentially regulated by the DE-miRNAs. Integrating DE-miRNAs, gene expression and plant-hormone metabolites indicated a regulation network where the downregulation of miRNAs at 3 dpi activated crucial resistance genes to trigger host hypersensitive cell death, improved hormone levels and upregulated the receptors/critical responsive transcription factors (TFs) of plant hormones, to shape immunity to the pathogen. Notably, our transcriptome, miRNA and hormone metabolites profiling and qPCR analysis suggested that that the downregulation of miR9472 potentially upregulated the expression of SAR Deficient 1 (SARD1), a key regulator for ICS1 (Isochorismate Synthase 1) induction and salicylic acid (SA) synthesis, to improve the level of SA in the Cf-10-gene-carrying line. Our results exploited potential regulatory network and new pathways underlying the resistance to C. fulvum in Cf-10-gene-carrying line, providing a more comprehensive genetic circuit and valuable gene targets for modulating resistance to the virus.

The evolutionary patterns, expression profiles, and genetic diversity of expanded genes in barley

Gene duplication resulting from whole-genome duplication (WGD), small-scale duplication (SSD), or unequal hybridization plays an important role in the expansion of gene families. Gene family expansion can also mediate species formation and adaptive evolution. Barley (Hordeum vulgare) is the world's fourth largest cereal crop, and it contains valuable genetic resources due to its ability to tolerate various types of environmental stress. In this study, 27,438 orthogroups in the genomes of seven Poaceae were identified, and 214 of them were significantly expanded in barley. The evolutionary rates, gene properties, expression profiles, and nucleotide diversity between expanded and non-expanded genes were compared. Expanded genes evolved more rapidly and experienced lower negative selection. Expanded genes, including their exons and introns, were shorter, they had fewer exons, their GC content was lower, and their first exons were longer compared with non-expanded genes. Codon usage bias was also lower for expanded genes than for non-expanded genes; the expression levels of expanded genes were lower than those of non-expanded genes, and the expression of expanded genes showed higher tissue specificity than that of non-expanded genes. Several stress-response-related genes/gene families were identified, and these genes could be used to breed barley plants with greater resistance to environmental stress. Overall, our analysis revealed evolutionary, structural, and functional differences between expanded and non-expanded genes in barley. Additional studies are needed to clarify the functions of the candidate genes identified in our study and evaluate their utility for breeding barley plants with greater stress resistance.

ANTHRACNOSE RESISTANCE GENE2 confers fungal resistance in sorghum

Sorghum is an important food and feed crop globally; its production is hampered by anthracnose disease caused by the fungal pathogen Colletotrichum sublineola (Cs). Here, we report identification and characterization of ANTHRACNOSE RESISTANCE GENE 2 (ARG2) encoding a nucleotide-binding leucine-rich repeat (NLR) protein that confers race-specific resistance to Cs strains. ARG2 is one of a cluster of several NLR genes initially identified in the sorghum differential line SC328C that is resistant to some Cs strains. This cluster shows structural and copy number variations in different sorghum genotypes. Different sorghum lines carrying independent ARG2 alleles provided the genetic validation for the identity of the ARG2 gene. ARG2 expression is induced by Cs, and chitin induces ARG2 expression in resistant but not in susceptible lines. ARG2-mediated resistance is accompanied by higher expression of defense and secondary metabolite genes at early stages of infection, and anthocyanin and zeatin metabolisms are upregulated in resistant plants. Interestingly, ARG2 localizes to the plasma membrane when transiently expressed in Nicotiana benthamiana. Importantly, ARG2 plants produced higher shoot dry matter than near-isogenic lines carrying the susceptible allele suggesting an absence of an ARG2 associated growth trade-off. Furthermore, ARG2-mediated resistance is stable at a wide range of temperatures. Our observations open avenues for resistance breeding and for dissecting mechanisms of resistance.

PM2b, a CC-NBS-LRR protein, interacts with TaWRKY76-D to regulate powdery mildew resistance in common wheat

Powdery mildew caused by Blumeria graminis f. sp. tritici (Bgt) is a destructive disease of wheat throughout the world. Host resistance is considered the most sustainable way to control this disease. Powdery mildew resistance gene Pm2b was mapped to the same genetic interval with Pm2a and PmCH1357 cloned previously, but showed different resistance spectra from them, indicating that they might be caused by different resistance genes or alleles. In this study, Pm2b was delimited to a 1.64 Mb physical interval using a large segregating population containing 4,354 F_2:3 families of resistant parent KM2939 and susceptible cultivar Shimai 15. In this interval, TraesCS5D03G0111700 encoding the coiled-coil nucleotide-binding site leucine-rich repeat protein (CC-NBS-LRR) was determined as the candidate gene of Pm2b. Silencing by barley stripe mosaic virus-induced gene silencing (BSMV-VIGS) technology and two independent mutants analysis in KM2939 confirmed the candidate gene TraesCS5D03G0111700 was Pm2b. The sequence of Pm2b was consistent with Pm2a/PmCH1357. Subcellular localization showed Pm2b was located on the cell nucleus and plasma membrane. Pm2b had the highest expression level in leaves and was rapidly up-regulated after inoculating with Bgt isolate E09. The yeast two-hybrid (Y2H) and luciferase complementation imaging assays (LCI) showed that PM2b could self-associate through the NB domain. Notably, we identified PM2b interacting with the transcription factor TaWRKY76-D, which depended on the NB domain of PM2b and WRKY domain of TaWRKY76-D. TaWRKY76-D negatively regulated the resistance to powdery mildew in wheat. The specific KASP marker K529 could take the advantage of high-throughput and high-efficiency for detecting Pm2b and be useful in molecular marker assisted-selection breeding. In conclusion, cloning and disease resistance mechanism analysis of Pm2b provided an example to emphasize a need of the molecular isolation of resistance genes, which has implications in marker assisted wheat breeding.

Mining of Cloned Disease Resistance Gene Homologs (CDRHs) in Brassica Species and Arabidopsis thaliana

Simple Summary

Developing cultivars with resistance genes (R genes) is an effective strategy to support high yield and quality in Brassica crops. The availability of clone R gene and genomic sequences in Brassica species and Arabidopsis thaliana provide the opportunity to compare genomic regions and survey R genes across genomic databases. In this paper, we aim to identify genes related to cloned genes through sequence identity, providing a repertoire of species-wide related R genes in Brassica crops. The comprehensive list of candidate R genes can be used as a reference for functional analysis.

Abstract

Various diseases severely affect Brassica crops, leading to significant global yield losses and a reduction in crop quality. In this study, we used the complete protein sequences of 49 cloned resistance genes (R genes) that confer resistance to fungal and bacterial diseases known to impact species in the Brassicaceae family. Homology searches were carried out across Brassica napus, B. rapa, B. oleracea, B. nigra, B. juncea, B. carinata and Arabidopsis thaliana genomes. In total, 660 cloned disease R gene homologs (CDRHs) were identified across the seven species, including 431 resistance gene analogs (RGAs) (248 nucleotide binding site-leucine rich repeats (NLRs), 150 receptor-like protein kinases (RLKs) and 33 receptor-like proteins (RLPs)) and 229 non-RGAs. Based on the position and distribution of specific homologs in each of the species, we observed a total of 87 CDRH clusters composed of 36 NLR, 16 RLK and 3 RLP homogeneous clusters and 32 heterogeneous clusters. The CDRHs detected consistently across the seven species are candidates that can be investigated for broad-spectrum resistance, potentially providing resistance to multiple pathogens. The R genes identified in this study provide a novel resource for the future functional analysis and gene cloning of Brassicaceae R genes towards crop improvement.

Evolution of resistance (R) gene specificity

Plant resistance (R) genes are members of large gene families with significant within and between species variation. It has been hypothesised that a variety of processes have shaped R gene evolution and the evolution of R gene specificity. In this review, we illustrate the main mechanisms that generate R gene diversity and provide examples of how they can change R gene specificity. Next, we explain which evolutionary mechanisms are at play and how they determine the fate of new R gene alleles and R genes. Finally, we place this in a larger context by comparing the diversity and evolution of R gene specificity within and between species scales.

Genetic diversity and disease: The past, present, and future of an old idea

Why do infectious diseases erupt in some host populations and not others? This question has spawned independent fields of research in evolution, ecology, public health, agriculture, and conservation. In the search for environmental and genetic factors that predict variation in parasitism, one hypothesis stands out for its generality and longevity: genetically homogeneous host populations are more likely to experience severe parasitism than genetically diverse populations. In this perspective piece, I draw on overlapping ideas from evolutionary biology, agriculture, and conservation to capture the far-reaching implications of the link between genetic diversity and disease. I first summarize the development of this hypothesis and the results of experimental tests. Given the convincing support for the protective effect of genetic diversity, I then address the following questions: (1) Where has this idea been put to use, in a basic and applied sense, and how can we better use genetic diversity to limit disease spread? (2) What new hypotheses does the established disease-diversity relationship compel us to test? I conclude that monitoring, preserving, and augmenting genetic diversity is one of our most promising evolutionarily informed strategies for buffering wild, domesticated, and human populations against future outbreaks.

Multiple Sources of Introduction of North American Arabidopsis thaliana from across Eurasia

Large-scale movement of organisms across their habitable range, or migration, is an important evolutionary process that can shape genetic diversity and influence the adaptive spread of alleles. Although human migrations have been studied in great detail with modern and ancient genomes, recent anthropogenic influence on reducing the biogeographical constraints on the migration of nonnative species has presented opportunities in several study systems to ask the questions about how repeated introductions shape genetic diversity in the introduced range. We present an extensive overview of population structure of North American Arabidopsis thaliana by studying a set of 500 whole-genome sequenced and over 2,800 RAD-seq genotyped individuals in the context of global diversity represented by Afro-Eurasian genomes. We use methods based on haplotype and rare-allele sharing as well as phylogenetic modeling to identify likely sources of introductions of extant N. American A. thaliana from the native range in Africa and Eurasia. We find evidence of admixture among the introduced lineages having increased haplotype diversity and reduced mutational load. We also detect signals of selection in immune-system-related genes that may impart qualitative disease resistance to pathogens of bacterial and oomycete origin. We conclude that multiple introductions to a nonnative range can rapidly enhance the adaptive potential of a colonizing species by increasing haplotypic diversity through admixture. Our results lay the foundation for further investigations into the functional significance of admixture.

Characterization of PmHHXM, a New Broad-Spectrum Powdery Mildew Resistance Gene in Chinese Wheat Landrace Honghuaxiaomai

Powdery mildew, caused by fungal pathogen Blumeria graminis f. sp. tritici, is an agronomically important and widespread wheat disease causing severe yield losses. Deployment of broad-spectrum disease resistance genes is the preferred strategy to prevent this pathogen. Chinese wheat landrace Honghuaxiaomai (HHXM) was resistant to all 23 tested B. graminis f. sp. tritici isolates at the seedling stage. The F₁, F₂, and F_2:3 progenies derived from the cross HHXM × Yangmai 158 were used in this study, and genetic analysis revealed that a single dominant gene, designated PmHHXM, conferred resistance to B. graminis f. sp. tritici isolate E09. Bulked segregant analysis and molecular mapping initially located PmHHXM to the distal region of chromosome 4AL. To fine map PmHHXM, we identified two critical recombinants from 592 F₂ plants and delimited PmHHXM to a 0.18-cM Xkasp475200 to Xhnu552 interval covering 1.77 Mb, in which a number of disease resistance-related gene clusters were annotated. Comparative mapping of this interval revealed a perturbed synteny among Triticeae species. This study reports the new powdery mildew resistance gene PmHHXM, which seems different from three known quantitative trait loci/genes identified on chromosome 4AL and has significant values for further genetic improvement. Analysis of the polymorphisms of 13 cosegregating markers between HHXM and 170 modern wheat cultivars indicates that Xhnu227 and Xsts478700 developed here are ideal for marker-assisted introgression of this locus in wheat breeding.

Genotype-Specific Expression and NLR Repertoire Contribute to Phenotypic Resistance Diversity in Plantago lanceolata

High levels of phenotypic variation in resistance appears to be nearly ubiquitous across natural host populations. Molecular processes contributing to this variation in nature are still poorly known, although theory predicts resistance to evolve at specific loci driven by pathogen-imposed selection. Nucleotide-binding leucine-rich repeat (NLR) genes play an important role in pathogen recognition, downstream defense responses and defense signaling. Identifying the natural variation in NLRs has the potential to increase our understanding of how NLR diversity is generated and maintained, and how to manage disease resistance. Here, we sequenced the transcriptomes of five different Plantago lanceolata genotypes when inoculated by the same strain of obligate fungal pathogen Podosphaera plantaginis. A de novo transcriptome assembly of RNA-sequencing data yielded 24,332 gene models with N50 value of 1,329 base pairs and gene space completeness of 66.5%. The gene expression data showed highly varying responses where each plant genotype demonstrated a unique expression profile in response to the pathogen, regardless of the resistance phenotype. Analysis on the conserved NB-ARC domain demonstrated a diverse NLR repertoire in P. lanceolata consistent with the high phenotypic resistance diversity in this species. We find evidence of selection generating diversity at some of the NLR loci. Jointly, our results demonstrate that phenotypic resistance diversity results from a crosstalk between different defense mechanisms. In conclusion, characterizing the architecture of resistance in natural host populations may shed unprecedented light on the potential of evolution to generate variation.

Analysis of intraspecies diversity reveals a subset of highly variable plant immune receptors and predicts their binding sites

The evolution of recognition specificities by the immune system depends on the generation of receptor diversity and on connecting the binding of new antigens with the initiation of downstream signaling. In plant immunity, the innate Nucleotide-Binding Leucine-Rich Repeat (NLR) receptor family enables antigen binding and immune signaling. In this study, we surveyed the NLR complements of 62 ecotypes of Arabidopsis thaliana and 54 lines of Brachypodium distachyon and identified a limited number of NLR subfamilies that show high allelic diversity. We show that the predicted specificity-determining residues cluster on the surfaces of Leucine-Rich Repeat domains, but the locations of the clusters vary among NLR subfamilies. By comparing NLR phylogeny, allelic diversity, and known functions of the Arabidopsis NLRs, we formulate a hypothesis for the emergence of direct and indirect pathogen-sensing receptors and of the autoimmune NLRs. These findings reveal the recurring patterns of evolution of innate immunity and can inform NLR engineering efforts.

Whole-genome assembly of Corylus avellana cv'Tonda Gentile delle Langhe' using linked-reads (10X Genomics)

The European hazelnut (Corylus avellana L.; 2n = 2x = 22) is a worldwide economically important tree nut that is cross-pollinated due to sporophytic incompatibility. Therefore, any individual plant is highly heterozygous. Cultivars are clonally propagated using mound layering, rooted suckers, and micropropagation. In recent years, the interest in this crop has increased, due to a growing demand related to the recognized health benefits of nut consumption. C. avellana cv “Tonda Gentile delle Langhe” (“TGdL”) is well-known for its high kernel quality, and the premium price paid for this cultivar is an economic benefit for producers in northern Italy. Assembly of a high-quality genome is a difficult task in many plant species because of the high level of heterozygosity. We assembled a chromosome-level genome sequence of “TGdL” with a two-step approach. First, 10X Genomics Chromium Technology was used to create a high-quality sequence, which was then assembled into scaffolds with cv “Tombul” genome as the reference. Eleven pseudomolecules were obtained, corresponding to 11 chromosomes. A total of 11,046 scaffolds remained unplaced, representing 11% of the genome (46,504,161 bp). Gene prediction, performed with Maker-P software, identified 27,791 genes (AED ≤0.4 and 92% of BUSCO completeness), whose function was analyzed with BlastP and InterProScan software. To characterize “TGdL” specific genetic mechanisms, Orthofinder was used to detect orthologs between hazelnut and closely related species. The “TGdL” genome sequence is expected to be a powerful tool to understand hazelnut genetics and allow detection of markers/genes for important traits to be used in targeted breeding programs.

Comparative Genomics and Functional Studies of Wheat BED-NLR Loci

Nucleotide-binding leucine-rich-repeat (LRR) receptors (NLRs) with non-canonical integrated domains (NLR-IDs) are widespread in plant genomes. Zinc-finger BED (named after the Drosophila proteins Boundary Element-Associated Factor and DNA Replication-related Element binding Factor, named BED hereafter) are among the most frequently found IDs. Five BED-NLRs conferring resistance against bacterial and fungal pathogens have been characterized. However, it is unknown whether BED-NLRs function in a manner similar to other NLR-IDs. Here, we used chromosome-level assemblies of wheat to explore the Yr7 and Yr5a genomic regions and show that, unlike known NLR-ID loci, there is no evidence for a NLR-partner in their vicinity. Using neighbor-network analyses, we observed that BED domains from BED-NLRs share more similarities with BED domains from single-BED proteins and from BED-containing proteins harboring domains that are conserved in transposases. We identified a nuclear localization signal (NLS) in Yr7, Yr5, and the other characterized BED-NLRs. We thus propose that this is a feature of BED-NLRs that confer resistance to plant pathogens. We show that the NLS was functional in truncated versions of the Yr7 protein when expressed in N. benthamiana. We did not observe cell-death upon the overexpression of Yr7 full-length, truncated, and 'MHD' variants in N. benthamiana. This suggests that either this system is not suitable to study BED-NLR signaling or that BED-NLRs require additional components to trigger cell death. These results define novel future directions to further understand the role of BED domains in BED-NLR mediated resistance.

What natural variation can teach us about resistance durability

Breeding a crop variety to be resistant to a pathogen usually takes years. This is problematic because pathogens, with short generation times and fluid genomes, adapt quickly to overcome resistance. The triumph of the pathogen is not inevitable, however, as there are numerous examples of durable resistance, particularly in wild plants. Which factors then contribute to such resistance stability over millennia? We review current knowledge of wild and agricultural pathosystems, detailing the importance of genetic, species and spatial heterogeneity in the prevention of pathogen outbreaks. We also highlight challenges associated with increasing resistance diversity in crops, both in light of pathogen (co-)evolution and breeding practices. Historically it has been difficult to incorporate heterogeneity into agriculture due to reduced efficiency in harvesting. Recent advances implementing computer vision and automation in agricultural production may improve our ability to harvest mixed genotype and mixed species plantings, thereby increasing resistance durability.

Identification of vital candidate microRNA/mRNA pairs regulating ovule development using high-throughput sequencing in hazel

BACKGROUND

Hazels (Corylus spp.) are economically important nut-producing species in which ovule development determines seed plumpness, one of the key parameters reflecting nut quality. microRNAs (miRNAs) play important roles in RNA silencing and the post-transcriptional regulation of gene expression. However, very little is currently known regarding the miRNAs involved in regulating ovule growth and development.

RESULTS

In this study, we accordingly sought to determine the important miRNAs involved in ovule development and growth in hazel. We examined ovules at four developmental stages, namely, ovule formation (Ov1), early ovule growth (Ov2), rapid ovule growth (Ov3), and ovule maturity (Ov4). On the basis of small RNA and mRNA sequencing using the Illumina sequencing platform, we identified 970 miRNAs in hazel, of which 766 and 204 were known and novel miRNAs, respectively. In Ov1-vs-Ov2, Ov1-vs-Ov3, Ov1-vs-Ov4, Ov2-vs-Ov3, Ov2-vs-Ov4, and Ov3-vs-Ov4 paired comparisons, 471 differentially expressed microRNAs (DEmiRNAs) and their 3117 target differentially expressed messenger RNAs (DEmRNAs) formed 11,199 DEmiRNA/DEmRNA pairs, with each DEmiRNA changing the expression of an average of 6.62 target mRNAs. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of all DEmRNAs revealed 29 significantly enriched KEGG pathways in the six paired comparisons, including protein export (ko03060), fatty acid elongation (ko00062), starch and sucrose metabolism (ko00500), fatty acid biosynthesis (ko00061), and amino sugar and nucleotide sugar metabolism (ko00520). Our results indicate that DEmiRNA/DEmRNA pairs showing opposite change trends were related to stress tolerance, embryo and seed development, cell proliferation, auxin transduction, and the biosynthesis of proteins, starch, and fats may participate in ovule growth and development.

CONCLUSIONS

These findings contribute to a better understanding of ovule development at the level of post-transcriptional regulation, and lay the foundation for further functional analyses of hazelnut ovule growth and development.

Transcriptome Analysis by RNA-Seq Reveals Genes Related to Plant Height in Two Sets of Parent-hybrid Combinations in Easter lily (Lilium longiflorum)

In this study, two different hybrids of Easter lily (Lilium longiflorum), obtained from two cross combinations, along with their four parents were sequenced by high–throughput RNA–sequencing (RNA–Seq) to find out differentially expressed gene in parent-hybrid combinations. The leaf mRNA profiles of two hybrids and their four parents were RNA–sequenced with a view to identify the potential candidate genes related to plant height heterosis. In both cross combinations, based to morphological traits mid–parent heterosis (MPH) was higher than high–parent heterosis (HPH) for plant height, leaf length, and number of flowers whereas HPH was higher than MPH for flowering time. A total of 4,327 differentially expressed genes (DEGs) were identified through RNA–Seq between the hybrids and their parents based on fold changes (FC) ≥ 2 for up– and ≤ –2 for down–regulation. Venn diagram analysis revealed that there were 703 common DEGs in two hybrid combinations, those were either up– or down–regulated. Most of the commonly expressed DEGs exhibited higher non–additive effects especially overdominance (75.9%) rather than additive (19.4%) and dominance (4.76%) effects. Among the 384 functionally annotated DEGs identified through Blast2GO tool, 12 DEGs were up–regulated and 16 of them were down–regulated in a similar fashion in both hybrids as revealed by heat map analysis. These 28 universally expressed DEGs were found to encode different types of proteins and enzymes those might regulate heterosis by modulating growth, development and stress–related functions in lily. In addition, gene ontology (GO) analysis of 260 annotated DEGs revealed that biological process might play dominant role in heterotic expression. In this first report of transcriptome sequencing in Easter lily, the notable universally up-regulated DEGs annotated ABC transporter A family member–like, B3 domain–containing, disease resistance RPP13/1, auxin–responsive SAUR68–like, and vicilin–like antimicrobial peptides 2–2 proteins those were perhaps associated with plant height heterosis. The genes expressed universally due to their overdominace function perhaps influenced MPH for greater plant height― largely by modulating biological processes involved therein. The genes identified in this study might be exploited in heterosis breeding for plant height of L. longiflorum.

Three previously characterized resistances to yellow rust are encoded by a single locus Wtk1

The wild emmer wheat (Triticum turgidum ssp. dicoccoides; WEW) yellow (stripe) rust resistance genes Yr15, YrG303, and YrH52 were discovered in natural populations from different geographic locations. They all localize to chromosome 1B but were thought to be non-allelic based on differences in resistance response. We recently cloned Yr15 as a Wheat Tandem Kinase 1 (WTK1) and show here that these three resistance loci co-segregate in fine-mapping populations and share an identical full-length genomic sequence of functional Wtk1. Independent ethyl methanesulfonate (EMS)-mutagenized susceptible yrG303 and yrH52 lines carried single nucleotide mutations in Wtk1 that disrupted function. A comparison of the mutations for yr15, yrG303, and yrH52 mutants showed that while key conserved residues were intact, other conserved regions in critical kinase subdomains were frequently affected. Thus, we concluded that Yr15-, YrG303-, and YrH52-mediated resistances to yellow rust are encoded by a single locus, Wtk1. Introgression of Wtk1 into multiple genetic backgrounds resulted in variable phenotypic responses, confirming that Wtk1-mediated resistance is part of a complex immune response network. WEW natural populations subjected to natural selection and adaptation have potential to serve as a good source for evolutionary studies of different traits and multifaceted gene networks.

Evolutionary dynamics and impacts of chromosome regions carrying R-gene clusters in rice

To elucidate R-gene evolution, we compared the genomic compositions and structures of chromosome regions carrying R-gene clusters among cultivated and wild rice species. Map-based sequencing and gene annotation of orthologous genomic regions (1.2 to 1.9 Mb) close to the terminal end of the long arm of rice chromosome 11 revealed R-gene clusters within six cultivated and ancestral wild rice accessions. NBS-LRR R-genes were much more abundant in Asian cultivated rice (O. sativa L.) than in its ancestors, indicating that homologs of functional genes involved in the same pathway likely increase in number because of tandem duplication of chromosomal segments and were selected during cultivation. Phylogenetic analysis using amino acid sequences indicated that homologs of paired Pikm1–Pikm2 (NBS-LRR) genes conferring rice-blast resistance were likely conserved among all cultivated and wild rice species we examined, and the homolog of Xa3/Xa26 (LRR-RLK) conferring bacterial blight resistance was lacking only in Kasalath.

Dominant, Heritable Resistance to Stewart's Wilt in Maize Is Associated with an Enhanced Vascular Defense Response to Infection with Pantoea stewartii

Vascular wilt bacteria such as Pantoea stewartii, the causal agent of Stewart's bacterial wilt of maize (SW), are destructive pathogens that are difficult to control. These bacteria colonize the xylem, where they form biofilms that block sap flow leading to characteristic wilting symptoms. Heritable forms of SW resistance exist and are used in maize breeding programs but the underlying genes and mechanisms are mostly unknown. Here, we show that seedlings of maize inbred lines with pan1 mutations are highly resistant to SW. However, current evidence suggests that other genes introgressed along with pan1 are responsible for resistance. Genomic analyses of pan1 lines were used to identify candidate resistance genes. In-depth comparison of P. stewartii interaction with susceptible and resistant maize lines revealed an enhanced vascular defense response in pan1 lines characterized by accumulation of electron-dense materials in xylem conduits visible by electron microscopy. We propose that this vascular defense response restricts P. stewartii spread through the vasculature, reducing both systemic bacterial colonization of the xylem network and consequent wilting. Though apparently unrelated to the resistance phenotype of pan1 lines, we also demonstrate that the effector WtsE is essential for P. stewartii xylem dissemination, show evidence for a nutritional immunity response to P. stewartii that alters xylem sap composition, and present the first analysis of maize transcriptional responses to P. stewartii infection.

Population Genetics of the Highly Polymorphic RPP8 Gene Family

Plant nucleotide-binding domain and leucine-rich repeat containing (NLR) genes provide some of the most extreme examples of polymorphism in eukaryotic genomes, rivalling even the vertebrate major histocompatibility complex. Surprisingly, this is also true in Arabidopsis thaliana, a predominantly selfing species with low heterozygosity. Here, we investigate how gene duplication and intergenic exchange contribute to this extraordinary variation. RPP8 is a three-locus system that is configured chromosomally as either a direct-repeat tandem duplication or as a single copy locus, plus a locus 2 Mb distant. We sequenced 48 RPP8 alleles from 37 accessions of A. thaliana and 12 RPP8 alleles from Arabidopsis lyrata to investigate the patterns of interlocus shared variation. The tandem duplicates display fixed differences and share less variation with each other than either shares with the distant paralog. A high level of shared polymorphism among alleles at one of the tandem duplicates, the single-copy locus and the distal locus, must involve both classical crossing over and intergenic gene conversion. Despite these polymorphism-enhancing mechanisms, the observed nucleotide diversity could not be replicated under neutral forward-in-time simulations. Only by adding balancing selection to the simulations do they approach the level of polymorphism observed at RPP8. In this NLR gene triad, genetic architecture, gene function and selection all combine to generate diversity.

Oomycetes Used in Arabidopsis Research

Arabidopsis plants in their natural environment are susceptible to infection by oomycete pathogens, in particular to downy mildew and white rust diseases. These naturally occurring infectious agents have imposed evolutionary pressures on Arabidopsis populations and are therefore highly relevant for the study of host-pathogen co-evolution. In addition, the study of oomycete diseases, including infections caused by several Phytophthora species, has led to many scientific discoveries on Arabidopsis immunity and disease. Herein, we describe the major oomycete species used for experiments on Arabidopsis, and how these pathosystems have been used to provide significant insights into mechanistic and evolutionary aspects of plant-oomycete interactions. We also highlight understudied aspects of plant-oomycete interactions, as well as translational approaches, that can be productively addressed using the reference pathosystems described in this article.

A Species-Wide Inventory of NLR Genes and Alleles in Arabidopsis thaliana

Infectious disease is both a major force of selection in nature and a prime cause of yield loss in agriculture. In plants, disease resistance is often conferred by nucleotide-binding leucine-rich repeat (NLR) proteins, intracellular immune receptors that recognize pathogen proteins and their effects on the host. Consistent with extensive balancing and positive selection, NLRs are encoded by one of the most variable gene families in plants, but the true extent of intraspecific NLR diversity has been unclear. Here, we define a nearly complete species-wide pan-NLRome in Arabidopsis thaliana based on sequence enrichment and long-read sequencing. The pan-NLRome largely saturates with approximately 40 well-chosen wild strains, with half of the pan-NLRome being present in most accessions. We chart NLR architectural diversity, identify new architectures, and quantify selective forces that act on specific NLRs and NLR domains. Our study provides a blueprint for defining pan-NLRomes.

A Species-Wide Inventory of NLR Genes and Alleles in Arabidopsis thaliana

Infectious disease is both a major force of selection in nature and a prime cause of yield loss in agriculture. In plants, disease resistance is often conferred by nucleotide-binding leucine-rich repeat (NLR) proteins, intracellular immune receptors that recognize pathogen proteins and their effects on the host. Consistent with extensive balancing and positive selection, NLRs are encoded by one of the most variable gene families in plants, but the true extent of intraspecific NLR diversity has been unclear. Here, we define a nearly complete species-wide pan-NLRome in Arabidopsis thaliana based on sequence enrichment and long-read sequencing. The pan-NLRome largely saturates with approximately 40 well-chosen wild strains, with half of the pan-NLRome being present in most accessions. We chart NLR architectural diversity, identify new architectures, and quantify selective forces that act on specific NLRs and NLR domains. Our study provides a blueprint for defining pan-NLRomes.

Identifying candidate genes for Phytophthora capsici resistance in pepper (Capsicum annuum) via genotyping-by-sequencing-based QTL mapping and genome-wide association study

Phytophthora capsici (Leon.) is a globally prevalent, devastating oomycete pathogen that causes root rot in pepper (Capsicum annuum). Several studies have identified quantitative trait loci (QTL) underlying resistance to P. capsici root rot (PcRR). However, breeding for pepper cultivars resistant to PcRR remains challenging due to the complexity of PcRR resistance. Here, we combined traditional QTL mapping with GWAS to broaden our understanding of PcRR resistance in pepper. Three major-effect loci (5.1, 5.2, and 5.3) conferring broad-spectrum resistance to three isolates of P. capsici were mapped to pepper chromosome P5. In addition, QTLs with epistatic interactions and minor effects specific to isolate and environment were detected on other chromosomes. GWAS detected 117 significant SNPs across the genome associated with PcRR resistance, including SNPs on chromosomes P5, P7, and P11 that colocalized with the QTLs identified here and in previous studies. Clusters of candidate nucleotide-binding site-leucine-rich repeat (NBS-LRR) and receptor-like kinase (RLK) genes were predicted within the QTL and GWAS regions; such genes often function in disease resistance. These candidate genes lay the foundation for the molecular dissection of PcRR resistance. SNP markers associated with QTLs for PcRR resistance will be useful for marker-assisted breeding and genomic selection in pepper breeding.

The AvrPm3-Pm3 effector-NLR interactions control both race-specific resistance and host-specificity of cereal mildews on wheat

The wheat Pm3 resistance gene against the powdery mildew pathogen occurs as an allelic series encoding functionally different immune receptors which induce resistance upon recognition of isolate-specific avirulence (AVR) effectors from the pathogen. Here, we describe the identification of five effector proteins from the mildew pathogens of wheat, rye, and the wild grass Dactylis glomerata, specifically recognized by the PM3B, PM3C and PM3D receptors. Together with the earlier identified AVRPM3A2/F2, the recognized AVRs of PM3B/C, (AVRPM3B2/C2), and PM3D (AVRPM3D3) belong to a large group of proteins with low sequence homology but predicted structural similarities. AvrPm3b2/c2 and AvrPm3d3 are conserved in all tested isolates of wheat and rye mildew, and non-host infection assays demonstrate that Pm3b, Pm3c, and Pm3d are also restricting the growth of rye mildew on wheat. Furthermore, divergent AVR homologues from non-adapted rye and Dactylis mildews are recognized by PM3B, PM3C, or PM3D, demonstrating their involvement in host specificity.

The Elaboration of miRNA Regulation and Gene Regulatory Networks in Plant⁻Microbe Interactions

Plants are exposed to diverse abiotic and biotic stimuli. These require fast and specific integrated responses. Such responses are coordinated at the protein and transcript levels and are incorporated into larger regulatory networks. Here, we focus on the evolution of transcriptional regulatory networks involved in plant–pathogen interactions. We discuss the evolution of regulatory networks and their role in fine-tuning plant defense responses. Based on the observation that many of the cornerstones of immune signaling in angiosperms are also present in streptophyte algae, it is likely that some regulatory components also predate the origin of land plants. The degree of functional conservation of many of these ancient components has not been elucidated. However, ongoing functional analyses in bryophytes show that some components are conserved. Hence, some of these regulatory components and how they are wired may also trace back to the last common ancestor of land plants or earlier. Of course, an understanding of the similarities and differences during the evolution of plant defense networks cannot ignore the lineage-specific coevolution between plants and their pathogens. In this review, we specifically focus on the small RNA regulatory networks involved in fine-tuning of the strength and timing of defense responses and highlight examples of pathogen exploitation of the host RNA silencing system. These examples illustrate well how pathogens frequently target gene regulation and thereby alter immune responses on a larger scale. That this is effective is demonstrated by the diversity of pathogens from distinct kingdoms capable of manipulating the same gene regulatory networks, such as the RNA silencing machinery.

Origin and evolution of the plant immune system

Contents Summary 70 I. Introduction 70 II. Ancient associations between plants and microbes 72 III. Evolutionary dynamics of plant-pathogen interactions 74 IV. Evolutionary signature of plant-pathogen interactions 74 V. Origin and evolution of RLK proteins 75 VI. Origin and evolution of NLR proteins 77 VII. Origin and evolution of SA signaling 78 VIII. Origin and evolution of RNA-based defense 79 IX. Perspectives 79 Acknowledgements 80 References 80 SUMMARY: Microbes have engaged in antagonistic associations with plants for hundreds of millions of years. Plants, in turn, have evolved diverse immune strategies to combat microbial pathogens. The conflicts between plants and pathogens result in everchanging coevolutionary cycles known as 'Red Queen' dynamics. These ancient and ongoing plant-pathogen interactions have shaped the evolution of both plant and pathogen genomes. With the recent explosion of plant genome-scale data, comparative analyses provide novel insights into the coevolutionary dynamics of plants and pathogens. Here, we discuss the ancient associations between plants and microbes as well as the evolutionary principles underlying plant-pathogen interactions. We synthesize and review the current knowledge on the origin and evolution of key components of the plant immune system. We also highlight the importance of studying algae and nonflowering land plants in understanding the evolution of the plant immune system.

Long-term balancing selection drives evolution of immunity genes in Capsella

Genetic drift is expected to remove polymorphism from populations over long periods of time, with the rate of polymorphism loss being accelerated when species experience strong reductions in population size. Adaptive forces that maintain genetic variation in populations, or balancing selection, might counteract this process. To understand the extent to which natural selection can drive the retention of genetic diversity, we document genomic variability after two parallel species-wide bottlenecks in the genus Capsella. We find that ancestral variation preferentially persists at immunity related loci, and that the same collection of alleles has been maintained in different lineages that have been separated for several million years. By reconstructing the evolution of the disease-related locus MLO2b, we find that divergence between ancient haplotypes can be obscured by referenced based re-sequencing methods, and that trans-specific alleles can encode substantially diverged protein sequences. Our data point to long-term balancing selection as an important factor shaping the genetics of immune systems in plants and as the predominant driver of genomic variability after a population bottleneck.

Small RNA molecules and their role in plant disease

All plant species are subject to disease. Plant diseases are caused by parasites, e.g. viruses, bacteria, oomycetes, parasitic plants, fungi, or nematodes. In all organisms, gene expression is tightly regulated and underpins essential functions and physiology. The coordination and regulation of both host and pathogen gene expression is essential for pathogens to infect and cause disease. One mode of gene regulation is RNA silencing. This biological process is widespread in the natural world, present in plants, animals and several pathogens. In RNA silencing, small (20-40 nucleotides) non-coding RNAs (small-RNAs, sRNAs) accumulate and regulate gene expression transcriptionally or post-transcriptionally in a sequence-specific manner. Regulation of sRNA molecules provides a fast mode to alter gene activity of multiple gene transcripts. RNA silencing is an ancient mechanism that protects the most sensitive part of an organism: its genetic code. sRNA molecules emerged as regulators of plant development, growth and plant immunity. sRNA based RNA silencing functions both within and between organisms. Here we present the described sRNAs from plants and pathogens and discuss how they regulate host immunity and pathogen virulence. We speculate on how sRNA molecules can be exploited to develop disease resistant plants. Finally, the activity of sRNA molecules can be prevented by proteins that suppress RNA silencing. This counter silencing response completes the dialog between plants and pathogens controlling plant disease or resistance outcome on the RNA (controlling gene expression) and protein level.

Interhomolog polymorphism shapes meiotic crossover within the Arabidopsis RAC1 and RPP13 disease resistance genes

During meiosis, chromosomes undergo DNA double-strand breaks (DSBs), which can be repaired using a homologous chromosome to produce crossovers. Meiotic recombination frequency is variable along chromosomes and tends to concentrate in narrow hotspots. We mapped crossover hotspots located in the Arabidopsis thaliana RAC1 and RPP13 disease resistance genes, using varying haplotypic combinations. We observed a negative non-linear relationship between interhomolog divergence and crossover frequency within the hotspots, consistent with polymorphism locally suppressing crossover repair of DSBs. The fancm, recq4a recq4b, figl1 and msh2 mutants, or lines with increased HEI10 dosage, are known to show increased crossovers throughout the genome. Surprisingly, RAC1 crossovers were either unchanged or decreased in these genetic backgrounds, showing that chromosome location and local chromatin environment are important for regulation of crossover activity. We employed deep sequencing of crossovers to examine recombination topology within RAC1, in wild type, fancm, recq4a recq4b and fancm recq4a recq4b backgrounds. The RAC1 recombination landscape was broadly conserved in the anti-crossover mutants and showed a negative relationship with interhomolog divergence. However, crossovers at the RAC1 5'-end were relatively suppressed in recq4a recq4b backgrounds, further indicating that local context may influence recombination outcomes. Our results demonstrate the importance of interhomolog divergence in shaping recombination within plant disease resistance genes and crossover hotspots.

Comparative Analysis of Impatiens Leaf Transcriptomes Reveal Candidate Genes for Resistance to Downy Mildew Caused by Plasmopara obducens

Impatiens downy mildew (IDM) is a devastating disease to garden impatiens. A good understanding of IDM resistance in New Guinea impatiens is essential for improving garden impatiens resistance to this disease. The present study was conducted to sequence, assemble, annotate and compare the leaf transcriptomes of two impatiens cultivars differing in resistance to IDM, reveal sequence polymorphisms and identify candidate genes for IDM resistance. RNA-Seq was performed on cultivars Super Elfin® XP Pink (SEP) and SunPatiens® Compact Royal Magenta (SPR). De novo assembly of obtained sequence reads resulted in 121,497 unigenes with an average length of 1156 nucleotides and N50 length of 1778 nucleotides. Searching the non-redundant protein and non-redundant nucleotide, Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes and Clusters of Orthologous Groups and Gene Ontology databases, resulted in annotation of 57.7% to 73.6% of the unigenes. Fifteen unigenes were highly similar to disease resistance genes and more abundant in the IDM-resistant cultivar than in the susceptible cultivar. A total of 22,484 simple sequence repeats (SSRs) and 245,936 and 120,073 single nucleotide polymorphisms (SNPs) were identified from SPR and SEP respectively. The assembled transcripts and unigenes, identified disease resistance genes and SSRs and SNPs sites will be a valuable resource for improving impatiens and its IDM resistance.

Restriction associated DNA-genotyping at multiple spatial scales in Arabidopsis lyrata reveals signatures of pathogen-mediated selection

Background

Genome scans based on outlier analyses have revolutionized detection of genes involved in adaptive processes, but reports of some forms of selection, such as balancing selection, are still limited. It is unclear whether high throughput genotyping approaches for identification of single nucleotide polymorphisms have sufficient power to detect modes of selection expected to result in reduced genetic differentiation among populations. In this study, we used Arabidopsis lyrata to investigate whether signatures of balancing selection can be detected based on genomic smoothing of Restriction Associated DNA sequencing (RAD-seq) data. We compared how different sampling approaches (both within and between subspecies) and different background levels of polymorphism (inbreeding or outcrossing populations) affected the ability to detect genomic regions showing key signatures of balancing selection, specifically elevated polymorphism, reduced differentiation and shifts towards intermediate allele frequencies. We then tested whether candidate genes associated with disease resistance (R-gene analogs) were detected more frequently in these regions compared to other regions of the genome.

Results

We found that genomic regions showing elevated polymorphism contained a significantly higher density of R-gene analogs predicted to be under pathogen-mediated selection than regions of non-elevated polymorphism, and that many of these also showed evidence for an intermediate site-frequency spectrum based on Tajima’s D. However, we found few genomic regions that showed both elevated polymorphism and reduced F_ST among populations, despite strong background levels of genetic differentiation among populations. This suggests either insufficient power to detect the reduced population structure predicted for genes under balancing selection using sparsely distributed RAD markers, or that other forms of diversifying selection are more common for the R-gene analogs tested.

Conclusions

Genome scans based on a small number of individuals sampled from a wide range of populations were sufficient to confirm the relative scarcity of signatures of balancing selection across the genome, but also identified new potential disease resistance candidates within genomic regions showing signatures of balancing selection that would be strong candidates for further sequencing efforts.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4806-7) contains supplementary material, which is available to authorized users.

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.

Genome-Wide Association Mapping Uncovers Fw1, a Dominant Gene Conferring Resistance to Fusarium Wilt in Strawberry

Fusarium wilt, a soil-borne disease caused by the fungal pathogen Fusarium oxysporum f. sp. fragariae, threatens strawberry (Fragaria × ananassa) production worldwide. The spread of the pathogen, coupled with disruptive changes in soil fumigation practices, have greatly increased disease pressure and the importance of developing resistant cultivars. While resistant and susceptible cultivars have been reported, a limited number of germplasm accessions have been analyzed, and contradictory conclusions have been reached in earlier studies to elucidate the underlying genetic basis of resistance. Here, we report the discovery of Fw1, a dominant gene conferring resistance to Fusarium wilt in strawberry. The Fw1 locus was uncovered in a genome-wide association study of 565 historically and commercially important strawberry accessions genotyped with 14,408 SNP markers. Fourteen SNPs in linkage disequilibrium with Fw1 physically mapped to a 2.3 Mb segment on chromosome 2 in a diploid F. vesca reference genome. Fw1 and 11 tightly linked GWAS-significant SNPs mapped to linkage group 2C in octoploid segregating populations. The most significant SNP explained 85% of the phenotypic variability and predicted resistance in 97% of the accessions tested-broad-sense heritability was 0.96. Several disease resistance and defense-related gene homologs, including a small cluster of genes encoding nucleotide-binding leucine-rich-repeat proteins, were identified in the 0.7 Mb genomic segment predicted to harbor Fw1 DNA variants and candidate genes identified in the present study should facilitate the development of high-throughput genotyping assays for accurately predicting Fusarium wilt phenotypes and applying marker-assisted selection.

Distinct domains of the AVRPM3A2/F2 avirulence protein from wheat powdery mildew are involved in immune receptor recognition and putative effector function

Recognition of the AVRPM3A2/F2 avirulence protein from powdery mildew by the wheat PM3A/F immune receptor induces a hypersensitive response after co-expression in Nicotiana benthamiana. The molecular determinants of this interaction and how they shape natural AvrPm3a2/f2 allelic diversity are unknown. We sequenced the AvrPm3a2/f2 gene in a worldwide collection of 272 mildew isolates. Using the natural polymorphisms of AvrPm3a2/f2 as well as sequence information from related gene family members, we tested 85 single-residue-altered AVRPM3A2/F2 variants with PM3A, PM3F and PM3FL456P/Y458H (modified for improved signaling) in Nicotiana benthamiana for effects on recognition. An intact AvrPm3a2/f2 gene was found in all analyzed isolates and the protein variant recognized by PM3A/F occurred globally at high frequencies. Single-residue alterations in AVRPM3A2/F2 mostly disrupted, but occasionally enhanced, the recognition response by PM3A, PM3F and PM3FL456P/Y458H . Residues enhancing hypersensitive responses constituted a protein domain separate from both naturally occurring polymorphisms and positively selected residues of the gene family. These results demonstrate the utility of using gene family sequence diversity to screen residues for their role in recognition. This approach identified a putative interaction surface in AVRPM3A2/F2 not polymorphic in natural alleles. We conclude that molecular mechanisms besides recognition drive AvrPm3a2/f2 diversification.

Proteomic Analysis Reveals Coordinated Regulation of Anthocyanin Biosynthesis through Signal Transduction and Sugar Metabolism in Black Rice Leaf

Black rice (Oryza sativa L.) is considered to be a healthy food due to its high content of anthocyanins in the pericarp. The synthetic pathway of anthocyanins in black rice grains has been identified, however, the proteomic profile of leaves during grain development is still unclear. Here, isobaric Tags Relative and Absolute Quantification (iTRAQ) MS/MS was carried out to identify statistically significant changes of leaf proteome in the black rice during grain development. Throughout three sequential developmental stages, a total of 3562 proteins were detected and 24 functional proteins were differentially expressed 3–10 days after flowering (DAF). The detected proteins are known to be involved in various biological processes and most of these proteins were related to gene expression regulatory (33.3%), signal transduction (16.7%) and developmental regulation and hormone-like proteins (12.5%). The coordinated changes were consistent with changes in regulatory proteins playing a leading role in leaves during black rice grain development. This indicated that signal transduction between leaves and grains may have an important role in anthocyanin biosynthesis and accumulation during grain development of black rice. In addition, four identified up-regulated proteins associated with starch metabolism suggested that the remobilization of nutrients for starch synthesis plays a potential role in anthocyanin biosynthesis of grain. The mRNA transcription for eight selected proteins was validated with quantitative real-time PCR. Our results explored the proteomics of the coordination between leaf and grain in anthocyanins biosynthesis of grain, which might be regulated by signal transduction and sugar metabolism in black rice leaf.

Promoter variants of Xa23 alleles affect bacterial blight resistance and evolutionary pattern

Bacterial blight, caused by Xanthomonas oryzae pv. oryzae (Xoo), is the most important bacterial disease in rice (Oryza sativa L.). Our previous studies have revealed that the bacterial blight resistance gene Xa23 from wild rice O. rufipogon Griff. confers the broadest-spectrum resistance against all the naturally occurring Xoo races. As a novel executor R gene, Xa23 is transcriptionally activated by the bacterial avirulence (Avr) protein AvrXa23 via binding to a 28-bp DNA element (EBE_AvrXa23) in the promoter region. So far, the evolutionary mechanism of Xa23 remains to be illustrated. Here, a rice germplasm collection of 97 accessions, including 29 rice cultivars (indica and japonica) and 68 wild relatives, was used to analyze the evolution, phylogeographic relationship and association of Xa23 alleles with bacterial blight resistance. All the ~ 473 bp DNA fragments consisting of promoter and coding regions of Xa23 alleles in the germplasm accessions were PCR-amplified and sequenced, and nine single nucleotide polymorphisms (SNPs) were detected in the promoter regions (~131 bp sequence upstream from the start codon ATG) of Xa23/xa23 alleles while only two SNPs were found in the coding regions. The SNPs in the promoter regions formed 5 haplotypes (Pro-A, B, C, D, E) which showed no significant difference in geographic distribution among these 97 rice accessions. However, haplotype association analysis indicated that Pro-A is the most favored haplotype for bacterial blight resistance. Moreover, SNP changes among the 5 haplotypes mostly located in the EBE/ebe regions (EBE_AvrXa23 and corresponding ebes located in promoters of xa23 alleles), confirming that the EBE region is the key factor to confer bacterial blight resistance by altering gene expression. Polymorphism analysis and neutral test implied that Xa23 had undergone a bottleneck effect, and selection process of Xa23 was not detected in cultivated rice. In addition, the Xa23 coding region was found highly conserved in the Oryza genus but absent in other plant species by searching the plant database, suggesting that Xa23 originated along with the diversification of the Oryza genus from the grass family during evolution. This research offers a potential for flexible use of novel Xa23 alleles in rice breeding programs and provide a model for evolution analysis of other executor R genes.

NLR diversity, helpers and integrated domains: making sense of the NLR IDentity

Plant innate immunity relies on genetically predetermined repertoires of immune receptors to detect pathogens and trigger an effective immune response. A large proportion of these receptors are from the Nucletoide Binding Leucine Rich Repeat (NLR) gene family. As plants live longer than most pathogens, maintaining diversity of NLRs and deploying efficient 'pathogen traps' is necessary to withstand the evolutionary battle. In this review, we summarize the sources of diversity in NLR plant immune receptors giving an overview of genomic, regulatory as well as functional studies, including the latest concepts of NLR helpers and NLRs with integrated domains.

Ecological genomics of tropical trees: how local population size and allelic diversity of resistance genes relate to immune responses, cosusceptibility to pathogens, and negative density dependence

In tropical forests, rarer species show increased sensitivity to species-specific soil pathogens and more negative effects of conspecific density on seedling survival (NDD). These patterns suggest a connection between ecology and immunity, perhaps because small population size disproportionately reduces genetic diversity of hyperdiverse loci such as immunity genes. In an experiment examining seedling roots from six species in one tropical tree community, we found that smaller populations have reduced amino acid diversity in pathogen resistance (R) genes but not the transcriptome in general. Normalized R gene amino acid diversity varied with local abundance and prior measures of differences in sensitivity to conspecific soil and NDD. After exposure to live soil, species with lower R gene diversity had reduced defence gene induction, more cosusceptibility of maternal cohorts to colonization by potentially pathogenic fungi, reduced root growth arrest (an R gene-mediated response) and their root-associated fungi showed lower induction of self-defence (antioxidants). Local abundance was not related to the ability to induce immune responses when pathogen recognition was bypassed by application of salicylic acid, a phytohormone that activates defence responses downstream of R gene signalling. These initial results support the hypothesis that smaller local tree populations have reduced R gene diversity and recognition-dependent immune responses, along with greater cosusceptibility to species-specific pathogens that may facilitate disease transmission and NDD. Locally rare species may be less able to increase their equilibrium abundance without genetic boosts to defence via immigration of novel R gene alleles from a larger and more diverse regional population.

Recombination Rate Heterogeneity within Arabidopsis Disease Resistance Genes

Meiotic crossover frequency varies extensively along chromosomes and is typically concentrated in hotspots. As recombination increases genetic diversity, hotspots are predicted to occur at immunity genes, where variation may be beneficial. A major component of plant immunity is recognition of pathogen Avirulence (Avr) effectors by resistance (R) genes that encode NBS-LRR domain proteins. Therefore, we sought to test whether NBS-LRR genes would overlap with meiotic crossover hotspots using experimental genetics in Arabidopsis thaliana. NBS-LRR genes tend to physically cluster in plant genomes; for example, in Arabidopsis most are located in large clusters on the south arms of chromosomes 1 and 5. We experimentally mapped 1,439 crossovers within these clusters and observed NBS-LRR gene associated hotspots, which were also detected as historical hotspots via analysis of linkage disequilibrium. However, we also observed NBS-LRR gene coldspots, which in some cases correlate with structural heterozygosity. To study recombination at the fine-scale we used high-throughput sequencing to analyze ~1,000 crossovers within the RESISTANCE TO ALBUGO CANDIDA1 (RAC1) R gene hotspot. This revealed elevated intragenic crossovers, overlapping nucleosome-occupied exons that encode the TIR, NBS and LRR domains. The highest RAC1 recombination frequency was promoter-proximal and overlapped CTT-repeat DNA sequence motifs, which have previously been associated with plant crossover hotspots. Additionally, we show a significant influence of natural genetic variation on NBS-LRR cluster recombination rates, using crosses between Arabidopsis ecotypes. In conclusion, we show that a subset of NBS-LRR genes are strong hotspots, whereas others are coldspots. This reveals a complex recombination landscape in Arabidopsis NBS-LRR genes, which we propose results from varying coevolutionary pressures exerted by host-pathogen relationships, and is influenced by structural heterozygosity.

Pooled Enrichment Sequencing Identifies Diversity and Evolutionary Pressures at NLR Resistance Genes within a Wild Tomato Population

Nod-like receptors (NLRs) are nucleotide-binding domain and leucine-rich repeats containing proteins that are important in plant resistance signaling. Many of the known pathogen resistance (R) genes in plants are NLRs and they can recognize pathogen molecules directly or indirectly. As such, divergence and copy number variants at these genes are found to be high between species. Within populations, positive and balancing selection are to be expected if plants coevolve with their pathogens. In order to understand the complexity of R-gene coevolution in wild nonmodel species, it is necessary to identify the full range of NLRs and infer their evolutionary history. Here we investigate and reveal polymorphism occurring at 220 NLR genes within one population of the partially selfing wild tomato species Solanum pennellii. We use a combination of enrichment sequencing and pooling ten individuals, to specifically sequence NLR genes in a resource and cost-effective manner. We focus on the effects which different mapping and single nucleotide polymorphism calling software and settings have on calling polymorphisms in customized pooled samples. Our results are accurately verified using Sanger sequencing of polymorphic gene fragments. Our results indicate that some NLRs, namely 13 out of 220, have maintained polymorphism within our S. pennellii population. These genes show a wide range of πN/πS ratios and differing site frequency spectra. We compare our observed rate of heterozygosity with expectations for this selfing and bottlenecked population. We conclude that our method enables us to pinpoint NLR genes which have experienced natural selection in their habitat.

R-gene variation across Arabidopsis lyrata subspecies: effects of population structure, selection and mating system

Background

Examining allelic variation of R-genes in closely related perennial species of Arabidopsis thaliana is critical to understanding how population structure and ecology interact with selection to shape the evolution of innate immunity in plants. We finely sampled natural populations of Arabidopsis lyrata from the Great Lakes region of North America (A. l. lyrata) and broadly sampled six European countries (A. l. petraea) to investigate allelic variation of two R-genes (RPM1 and WRR4) and neutral genetic markers (Restriction Associated DNA sequences and microsatellites) in relation to mating system, phylogeographic structure and subspecies divergence.

Results

Fine-scale sampling of populations revealed strong effects of mating system and population structure on patterns of polymorphism for both neutral loci and R-genes, with no strong evidence for selection. Broad geographic sampling revealed evidence of balancing selection maintaining polymorphism in R-genes, with elevated heterozygosity and diversity compared to neutral expectations and sharing of alleles among diverged subspecies. Codon-based tests detected both positive and purifying selection for both R-genes, as commonly found for animal immune genes.

Conclusions

Our results highlight that combining fine and broad-scale sampling strategies can reveal the multiple factors influencing polymorphism and divergence at potentially adaptive genes such as R-genes.

Electronic supplementary material

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

Diversity and Evolution of Disease Resistance Genes in Barley (Hordeum vulgare L.)

Plant disease resistance genes (R-genes) play a critical role in the defense response to pathogens. Barley is one of the most important cereal crops, having a genome recently made available, for which the diversity and evolution of R-genes are not well understood. The main objectives of this research were to conduct a genome-wide identification of barley Coiled-coil, Nucleotide-binding site, Leucine-rich repeat (CNL) genes and elucidate their evolutionary history. We employed a Hidden Markov Model using 52 Arabidopsis thaliana CNL reference sequences and analyzed for phylogenetic relationships, structural variation, and gene clustering. We identified 175 barley CNL genes nested into three clades, showing (a) evidence of an expansion of the CNL-C clade, primarily due to tandem duplications; (b) very few members of clade CNL-A and CNL-B; and (c) a complete absence of clade CNL-D. Our results also showed that several of the previously identified mildew locus A (MLA) genes may be allelic variants of two barley CNL genes, MLOC_66581 and MLOC_10425, which respond to powdery mildew. Approximately 23% of the barley CNL genes formed 15 gene clusters located in the extra-pericentromeric regions on six of the seven chromosomes; more than half of the clustered genes were located on chromosomes 1H and 7H. Higher average numbers of exons and multiple splice variants in barley relative to those in Arabidopsis and rice may have contributed to a diversification of the CNL-C members. These results will help us understand the evolution of R-genes with potential implications for developing durable resistance in barley cultivars.

The Genetics Underlying Natural Variation in the Biotic Interactions of Arabidopsis thaliana: The Challenges of Linking Evolutionary Genetics and Community Ecology

In the context of global change, predicting the responses of plant communities in an ever-changing biotic environment calls for a multipronged approach at the interface of evolutionary genetics and community ecology. However, our understanding of the genetic basis of natural variation involved in mediating biotic interactions, and associated adaptive dynamics of focal plants in their natural communities, is still in its infancy. Here, we review the genetic and molecular bases of natural variation in the response to biotic interactions (viruses, bacteria, fungi, oomycetes, herbivores, and plants) in the model plant Arabidopsis thaliana as well as the adaptive value of these bases. Among the 60 identified genes are a number that encode nucleotide-binding site leucine-rich repeat (NBS-LRR)-type proteins, consistent with early examples of plant defense genes. However, recent studies have revealed an extensive diversity in the molecular mechanisms of defense. Many types of genetic variants associate with phenotypic variation in biotic interactions, even among the genes of large effect that tend to be identified. In general, we found that (i) balancing selection rather than directional selection explains the observed patterns of genetic diversity within A. thaliana and (ii) the cost/benefit tradeoffs of adaptive alleles can be strongly dependent on both genomic and environmental contexts. Finally, because A. thaliana rarely interacts with only one biotic partner in nature, we highlight the benefit of exploring diffuse biotic interactions rather than tightly associated host-enemy pairs. This challenge would help to improve our understanding of coevolutionary quantitative genetics within the context of realistic community complexity.