Clusters of resistance genes in plants evolve by divergent selection and a birth-and-death process

Comparative Genomics Analysis in Grass Species Reveals Two Distinct Evolutionary Strategies Adopted by R Genes

Resistance genes play an important role in the defense of plants against the invasion of pathogens. In Setaria italica and closely related grass species, R genes have been identified through genetic mapping and genome-wide homologous/domain searching. However, there has been to date no systematic analysis of the evolutionary features of R genes across all sequenced grass genomes. Here, we determined and comprehensively compared R genes in all 12 assembled grass genomes and an outgroup species (Arabidopsis thaliana) through synteny and selection analyses of multiple genomes. We found that the two groups of nucleotide binding site (NBS) domains containing R genes—R tandem duplications (TD) and R singletons—adopted different strategies and showed different features in their evolution. Based on K_a/K_s analysis between syntenic R loci pairs of TDs or singletons, we conclude that R singletons are under stronger purifying selection to be conserved among different grass species than R TDs, while R genes located at TD arrays have evolved much faster through diversifying selection. Furthermore, using the variome datasets of S. italica populations, we scanned for selection signals on genes and observed that a part of R singleton genes have been under purifying selection in populations of S. italica, which is consistent with the pattern observed in syntenic R singletons among different grass species. Additionally, we checked the synteny relationships of reported R genes in grass species and found that the functionally mapped R genes for novel resistance traits are prone to appear in TDs and are heavily divergent from their syntenic orthologs in other grass species, such the black streak R gene Rxo1 in Z. mays and the blast R gene Pi37 in O. sativa. These findings indicate that the R genes from TDs adopted tandem duplications to evolve faster and accumulate more mutations to facilitate functional innovation to cope with variable threats from a fluctuating environment, while R singletons provide a way for R genes to maintain sequence stability and retain conservation of function.

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.

Detached Leaf Assays for Resistance to Crown Rust Reveal Diversity Within Populations of Avena sterilis

Crown rust is the most widespread and damaging disease of oat (Avena species). Genetic resistance to the pathogen is the preferred method for crop protection but widespread deployment of limited numbers of major effect genes has promoted the rapid emergence and spread of pathogen races that are able to overcome these genes. Combining genes with even partial resistance may help develop durable cultivars that are less vulnerable to changes in pathogen virulence. Partial resistance is expected to be relatively common in populations of wild species where constant pathogen pressure encourages diversity in host resistance mechanisms, but it may be discarded in conventional screens for major gene resistance. Here, we used a detached leaf assay to detect resistance to the crown rust pathogen, Puccinia coronata Cda. f. sp. avenae, in previously uncharacterized collections of the hexaploid wild oat relative A. sterilis made by the Polish National Centre for Plant Genetic Resources. Many of the accessions were collected in Morocco, the center of diversity for the Avena genus. The detached leaf assessment allowed individual plants to be challenged with multiple pathotypes and their responses compared with 34 known differentials. Broad-spectrum resistance was identified within accession PL 51855, which behaved as a single major locus on crossing to three cultivars. The locus provided resistance to over 50 rust pathotypes, a greater range than seen for any of the known host resistance (Pc) genes. Strong resistance was identified in other accessions, and heterogeneity in response within accessions was common. Several accessions show multiple partial resistance responses that may be of value for developing durable resistance in cultivars. Because the sources of resistance in all but two differential lines were collected outside of Morocco, resistance in all accessions tested here are potentially novel. This study demonstrates that diversity within A. sterilis accessions collected in Morocco could be a very valuable source of resistance to crown rust, and it provides new germplasm for use in resistance breeding programs. Detached leaf assessment provides a valuable first step in the identification of promising candidates in complex gene bank accessions.

Variation in abundance of predicted resistance genes in the Brassica oleracea pangenome

Brassica oleracea is an important agricultural species encompassing many vegetable crops including cabbage, cauliflower, broccoli and kale; however, it can be susceptible to a variety of fungal diseases such as clubroot, blackleg, leaf spot and downy mildew. Resistance to these diseases is meditated by specific disease resistance genes analogs (RGAs) which are differently distributed across B. oleracea lines. The sequenced reference cultivar does not contain all B. oleracea genes due to gene presence/absence variation between individuals, which makes it necessary to search for RGA candidates in the B. oleracea pangenome. Here we present a comparative analysis of RGA candidates in the pangenome of B. oleracea. We show that the presence of RGA candidates differs between lines and suggests that in B. oleracea, SNPs and presence/absence variation drive RGA diversity using separate mechanisms. We identified 59 RGA candidates linked to Sclerotinia, clubroot, and Fusarium wilt resistance QTL, and these findings have implications for crop breeding in B. oleracea, which may also be applicable in other crops species.

Divergence between bread wheat and Triticum militinae in the powdery mildew resistance QPm.tut-4A locus and its implications for cloning of the resistance gene

A segment of Triticum militinae chromosome 7G harbors a gene(s) conferring powdery mildew resistance which is effective at both the seedling and the adult plant stages when transferred into bread wheat (T. aestivum). The introgressed segment replaces a piece of wheat chromosome arm 4AL. An analysis of segregating materials generated to positionally clone the gene highlighted that in a plant heterozygous for the introgression segment, only limited recombination occurs between the introgressed region and bread wheat 4A. Nevertheless, 75 genetic markers were successfully placed within the region, thereby confining the gene to a 0.012 cM window along the 4AL arm. In a background lacking the Ph1 locus, the localized rate of recombination was raised 33-fold, enabling the reduction in the length of the region containing the resistance gene to a 480 kbp stretch harboring 12 predicted genes. The substituted segment in the reference sequence of bread wheat cv. Chinese Spring is longer (640 kbp) and harbors 16 genes. A comparison of the segments' sequences revealed a high degree of divergence with respect to both their gene content and nucleotide sequence. Of the 12 T. militinae genes, only four have a homolog in cv. Chinese Spring. Possible candidate genes for the resistance have been identified based on function predicted from their sequence.

Reconstruction of evolutionary trajectories of chromosomes unraveled independent genomic repatterning between Triticeae and Brachypodium

BACKGROUND

After polyploidization, a genome may experience large-scale genome-repatterning, featuring wide-spread DNA rearrangement and loss, and often chromosome number reduction. Grasses share a common tetraploidization, after which the originally doubled chromosome numbers reduced to different chromosome numbers among them. A telomere-centric reduction model was proposed previously to explain chromosome number reduction. With Brachpodium as an intermediate linking different major lineages of grasses and a model plant of the Pooideae plants, we wonder whether it mediated the evolution from ancestral grass karyotype to Triticeae karyotype.

RESULTS

By inferring the homology among Triticeae, rice, and Brachpodium chromosomes, we reconstructed the evolutionary trajectories of the Triticeae chromosomes. By performing comparative genomics analysis with rice as a reference, we reconstructed the evolutionary trajectories of Pooideae plants, including Ae. Tauschii (2n = 14, DD), barley (2n = 14), Triticum turgidum (2n = 4x = 28, AABB), and Brachypodium (2n = 10). Their extant Pooidea and Brachypodium chromosomes were independently produced after sequential nested chromosome fusions in the last tens of millions of years, respectively, after their split from rice. More frequently than would be expected by chance, in Brachypodium, the 'invading' and 'invaded' chromosomes are homoeologs, originating from duplication of a common ancestral chromosome, that is, with more extensive DNA-level correspondence to one another than random chromosomes, nested chromosome fusion events between homoeologs account for three of seven cases in Brachypodium (P-value≈0.00078). However, this phenomenon was not observed during the formation of other Pooideae chromosomes.

CONCLUSIONS

Notably, we found that the Brachypodium chromosomes formed through exclusively distinctive trajectories from those of Pooideae plants, and were well explained by the telomere-centric model. Our work will contribute to understanding the structural and functional innovation of chromosomes in different Pooideae lineages and beyond.

Comparative genomics and transcriptomics of host-pathogen interactions in insects: evolutionary insights and future directions

Classical evolutionary studies of protein-coding genes have established that genes in the canonical immune system are often among the most rapidly evolving within and between species. As more genomes and transcriptomes across insects are sequenced, it is becoming clear that duplications and losses of immune genes are also a likely consequence of host-pathogen interactions. Furthermore, particular species respond to diverse pathogenic challenges with a wide range of challenge-specific responses that are still poorly understood. Transcriptional studies, using RNA-seq to characterize the infection-regulated transcriptome of diverse insects, are crucial for additional progress in understanding the ecology and evolution of the full complexity of the host response.

Quantitative trait loci controlling Phytophthora cactorum resistance in the cultivated octoploid strawberry (Fragaria × ananassa)

The cultivated strawberry, Fragaria × ananassa  (Fragaria spp.) is the most economically important global soft fruit. Phytophthora cactorum, a water-borne oomycete causes economic losses in strawberry production globally. A bi-parental cross of octoploid cultivated strawberry segregating for resistance to P. cactorum, the causative agent of crown rot disease, was screened using artificial inoculation. Multiple putative resistance quantitative trait loci (QTL) were identified and mapped. Three major effect QTL (FaRPc6C, FaRPc6D and FaRPc7D) explained 37% of the variation observed. There were no epistatic interactions detected between the three major QTLs. Testing a subset of the mapping population progeny against a range of P. cactorum isolates revealed no significant interaction (p = 0.0593). However, some lines showed higher susceptibility than predicted, indicating that additional undetected factors may affect the expression of some quantitative resistance loci. Using historic crown rot disease score data from strawberry accessions, a preliminary genome-wide association study (GWAS) of 114 individuals revealed an additional locus associated with resistance to P. cactorum. Mining of the Fragaria vesca Hawaii 4 v1.1 genome revealed candidate resistance genes in the QTL regions.

Physical Localization of the Root-Knot Nematode (Meloidogyne incognita) Resistance Locus Me7 in Pepper (Capsicum annuum)

The root-knot nematode (RKN) Meloidogyne incognita severely reduces yields of pepper (Capsicum annuum) worldwide. A single dominant locus, Me7, conferring RKN resistance was previously mapped on the long arm of pepper chromosome P9. In the present study, the Me7 locus was fine mapped using an F₂ population of 714 plants derived from a cross between the RKN-susceptible parent C. annuum ECW30R and the RKN-resistant parent C. annuum CM334. CM334 exhibits suppressed RKN juvenile movement, suppressed feeding site enlargement and significant reduction in gall formation compared with ECW30R. RKN resistance screening in the F₂ population identified 558 resistant and 156 susceptible plants, which fit a 3:1 ratio confirming that this RKN resistance was controlled by a single dominant gene. Using the C. annuum CM334 reference genome and BAC library sequencing, fine mapping of Me7 markers was performed. The Me7 locus was delimited between two markers G21U3 and G43U3 covering a physical interval of approximately 394.7 kb on the CM334 chromosome P9. Nine markers co-segregated with the Me7 gene. A cluster of 25 putative nucleotide-binding site and leucine-rich repeat (NBS-LRR)-type disease resistance genes were predicted in the delimited Me7 region. We propose that RKN resistance in CM334 is mediated by one or more of these NBS-LRR class R genes. The Me7-linked markers identified here will facilitate marker-assisted selection (MAS) for RKN resistance in pepper breeding programs, as well as functional analysis of Me7 candidate genes in C. annuum.

Mapping a gene on wheat chromosome 4BL involved in a complementary interaction with adult plant leaf rust resistance gene LrSV2

A complementary gene to LrSV2 for specific adult plant leaf rust resistance in wheat was mapped on chromosome 4BL, tightly linked to Lr12 / 31. LrSV2 is a race-specific adult plant leaf rust (Puccinia triticina) resistance gene on subdistal chromosome 3BS detected in the cross of the traditional Argentinean wheat (Triticum aestivum) variety Sinvalocho MA and the experimental line Gama6. The analysis of the cross of R46 [recombinant inbred line (RIL) derived from Sinvalocho MA carrying LrSV2 gene and the complementary gene Lrc-SV2 identified in the current paper] and the commercial variety Relmo Siriri (not carrying neither of these two genes) allowed the detection of the unlinked complementary gene Lrc-SV2 because the presence of one dominant allele of both is necessary to express the LrSV2-specific adult plant resistance. Lrc-SV2 was mapped within a 1-cM interval on chromosome 4BL using 100 RILs from the cross Sinvalocho MA × Purple Straw. This genetic system resembles the Lr27+31 seedling resistance reported in the Australian varieties Gatcher and Timgalen where interacting genes map at similar chromosomal positions. However, in high-resolution maps, Lr27 and LrSV2 were already mapped to adjacent intervals on 3BS and Lrc-SV2 map position on 4BL is distal to the reported Lr12/31-flanking microsatellites.

Genome-wide scans of selection highlight the impact of biotic and abiotic constraints in natural populations of the model grass Brachypodium distachyon

Grasses are essential plants for ecosystem functioning. Quantifying the selective pressures that act on natural variation in grass species is therefore essential regarding biodiversity maintenance. In this study, we investigate the selection pressures that act on two distinct populations of the grass model Brachypodium distachyon without prior knowledge about the traits under selection. We took advantage of whole-genome sequencing data produced for 44 natural accessions of B. distachyon and used complementary genome-wide selection scans (GWSS) methods to detect genomic regions under balancing and positive selection. We show that selection is shaping genetic diversity at multiple temporal and spatial scales in this species, and affects different genomic regions across the two populations. Gene ontology annotation of candidate genes reveals that pathogens may constitute important factors of positive and balancing selection in B. distachyon. We eventually cross-validated our results with quantitative trait locus data available for leaf-rust resistance in this species and demonstrate that, when paired with classical trait mapping, GWSS can help pinpointing candidate genes for further molecular validation. Thanks to a near base-perfect reference genome and the large collection of freely available natural accessions collected across its natural range, B. distachyon appears as a prime system for studies in ecology, population genomics and evolutionary biology.

Ma Orthologous Genes in Prunus spp. Shed Light on a Noteworthy NBS-LRR Cluster Conferring Differential Resistance to Root-Knot Nematodes

Root-knot nematodes (RKNs) are considerable polyphagous pests that severely challenge plants worldwide and especially perennials. The specific genetic resistance of plants mainly relies on the NBS-LRR genes that are pivotal factors for pathogens control. In Prunus spp., the Ma plum and RMja almond genes possess different spectra for resistance to RKNs. While previous works based on the Ma gene allowed to clone it and to decipher its peculiar TIR-NBS-LRR (TNL) structure, we only knew that the RMja gene mapped on the same chromosome as Ma. We carried out a high-resolution mapping using an almond segregating F2 progeny of 1448 seedlings from resistant (R) and susceptible (S) parental accessions, to locate precisely RMja on the peach genome, the reference sequence for Prunus species. We showed that the RMja gene maps in the Ma resistance cluster and that the Ma ortholog is the best candidate for RMja. This co-localization is a crucial step that opens the way to unravel the molecular determinants involved in the resistance to RKNs. Then we sequenced both almond parental NGS genomes and aligned them onto the RKN susceptible reference peach genome. We produced a BAC library of the R parental accession and, from two overlapping BAC clones, we obtained a 336-kb sequence encompassing the RMja candidate region. Thus, we could benefit from three Ma orthologous regions to investigate their sequence polymorphism, respectively, within plum (complete R spectrum), almond (incomplete R spectrum) and peach (null R spectrum). We showed that the Ma TNL cluster has evolved orthologs with a unique conserved structure comprised of five repeated post-LRR (PL) domains, which contain most polymorphism. In addition to support the Ma and RMja orthologous relationship, our results suggest that the polymorphism contained in the PL sequences might underlie differential resistance interactions with RKNs and an original immune mechanism in woody perennials. Besides, our study illustrates how PL exon duplications and losses shape TNL structure and give rise to atypical PL domain repeats of yet unknown role.

Evolutionary Divergence of TNL Disease-Resistant Proteins in Soybean (Glycine max) and Common Bean (Phaseolus vulgaris)

Disease-resistant genes (R genes) encode proteins that are involved in protecting plants from their pathogens and pests. Availability of complete genome sequences from soybean and common bean allowed us to perform a genome-wide identification and analysis of the Toll interleukin-1 receptor-like nucleotide-binding site leucine-rich repeat (TNL) proteins. Hidden Markov model (HMM) profiling of all protein sequences resulted in the identification of 117 and 77 regular TNL genes in soybean and common bean, respectively. We also identified TNL gene homologs with unique domains, and signal peptides as well as nuclear localization signals. The TNL genes in soybean formed 28 clusters located on 10 of the 20 chromosomes, with the majority found on chromosome 3, 6 and 16. Similarly, the TNL genes in common bean formed 14 clusters located on five of the 11 chromosomes, with the majority found on chromosome 10. Phylogenetic analyses of the TNL genes from Arabidopsis, soybean and common bean revealed less divergence within legumes relative to the divergence between legumes and Arabidopsis. Syntenic blocks were found between chromosomes Pv10 and Gm03, Pv07 and Gm10, as well as Pv01 and Gm14. The gene expression data revealed basal level expression and tissue specificity, while analysis of available microRNA data showed 37 predicted microRNA families involved in targeting the identified TNL genes in soybean and common bean.

The tomato subtilase family includes several cell death-related proteinases with caspase specificity

Phytaspases are Asp-specific subtilisin-like plant proteases that have been likened to animal caspases with respect to their regulatory function in programmed cell death (PCD). We identified twelve putative phytaspase genes in tomato that differed widely in expression level and tissue-specific expression patterns. Most phytaspase genes are tandemly arranged on tomato chromosomes one, four, and eight, and many belong to taxon-specific clades, e.g. the P69 clade in the nightshade family, suggesting that these genes evolved by gene duplication after speciation. Five tomato phytaspases (SlPhyts) were expressed in N. benthamiana and purified to homogeneity. Substrate specificity was analyzed in a proteomics assay and with a panel of fluorogenic peptide substrates. Similar to animal caspases, SlPhyts recognized an extended sequence motif including Asp at the cleavage site. Clear differences in cleavage site preference were observed implying different substrates in vivo and, consequently, different physiological functions. A caspase-like function in PCD was confirmed for five of the seven tested phytaspases. Cell death was triggered by ectopic expression of SlPhyts 2, 3, 4, 5, 6 in tomato leaves by agro-infiltration, as well as in stably transformed transgenic tomato plants. SlPhyts 3, 4, and 5 were found to contribute to cell death under oxidative stress conditions.

Assembly and annotation of a draft genome sequence for Glycine latifolia, a perennial wild relative of soybean

Glycine latifolia (Benth.) Newell & Hymowitz (2n = 40), one of the 27 wild perennial relatives of soybean, possesses genetic diversity and agronomically favorable traits that are lacking in soybean. Here, we report the 939-Mb draft genome assembly of G. latifolia (PI 559298) using exclusively linked-reads sequenced from a single Chromium library. We organized scaffolds into 20 chromosome-scale pseudomolecules utilizing two genetic maps and the Glycine max (L.) Merr. genome sequence. High copy numbers of putative 91-bp centromere-specific tandem repeats were observed in consecutive blocks within predicted pericentromeric regions on several pseudomolecules. No 92-bp putative centromeric repeats, which are abundant in G. max, were detected in G. latifolia or Glycine tomentella. Annotation of the assembled genome and subsequent filtering yielded a high confidence gene set of 54 475 protein-coding loci. In comparative analysis with five legume species, genes related to defense responses were significantly overrepresented in Glycine-specific orthologous gene families. A total of 304 putative nucleotide-binding site (NBS)-leucine-rich-repeat (LRR) genes were identified in this genome assembly. Different from other legume species, we observed a scarcity of TIR-NBS-LRR genes in G. latifolia. The G. latifolia genome was also predicted to contain genes encoding 367 LRR-receptor-like kinases, a family of proteins involved in basal defense responses and responses to abiotic stress. The genome sequence and annotation of G. latifolia provides a valuable source of alternative alleles and novel genes to facilitate soybean improvement. This study also highlights the efficacy and cost-effectiveness of the application of Chromium linked-reads in diploid plant genome de novo assembly.

The island cotton NBS-LRR gene GbaNA1 confers resistance to the non-race 1 Verticillium dahliae isolate Vd991

Wilt caused by Verticillium dahliae significantly reduces cotton yields, as host resistance in commercially cultivated Gossypium species is lacking. Understanding the molecular basis of disease resistance in non-commercial Gossypium species could galvanize the development of Verticillium wilt resistance in cultivated species. Nucleotide-binding site leucine-rich repeat (NBS-LRR) proteins play a central role in plant defence against pathogens. In this study, we focused on the relationship between a locus enriched with eight NBS-LRR genes and Verticillium wilt resistance in G. barbadense. Independent virus-induced gene silencing of each of the eight NBS-LRR genes in G. barbadense cultivar Hai 7124 revealed that silencing of GbaNA1 alone compromised the resistance of G. barbadense to V. dahliae isolate Vd991. In cultivar Hai 7124, GbaNA1 could be induced by V. dahliae isolate Vd991 and by ethylene, jasmonic acid and salicylic acid. Nuclear protein localization of GbaNA1 was demonstrated by transient expression. Sequencing of the GbaNA1 orthologue in nine G. hirsutum accessions revealed that all carried a non-functional allele, caused by a premature peptide truncation. In addition, all 10 G. barbadense and nine G. hirsutum accessions tested carried a full-length (∼1140 amino acids) homologue of the V. dahliae race 1 resistance gene Gbve1, although some sequence polymorphisms were observed. Verticillium dahliae Vd991 is a non-race 1 isolate that lacks the Ave1 gene. Thus, the resistance imparted by GbaNA1 appears to be mediated by a mechanism distinct from recognition of the fungal effector Ave1.

Identification and rapid mapping of a gene conferring broad-spectrum late blight resistance in the diploid potato species Solanum verrucosum through DNA capture technologies

A broad-spectrum late blight disease-resistance gene from Solanum verrucosum has been mapped to potato chromosome 9. The gene is distinct from previously identified-resistance genes. We have identified and characterised a broad-spectrum resistance to Phytophthora infestans from the wild Mexican species Solanum verrucosum. Diagnostic resistance gene enrichment (dRenSeq) revealed that the resistance is not conferred by previously identified nucleotide-binding, leucine-rich repeat genes. Utilising the sequenced potato genome as a reference, two complementary enrichment strategies that target resistance genes (RenSeq) and single/low-copy number genes (Generic-mapping enrichment Sequencing; GenSeq), respectively, were deployed for the rapid, SNP-based mapping of the resistance through bulked-segregant analysis. Both approaches independently positioned the resistance, referred to as Rpi-ver1, to the distal end of potato chromosome 9. Stringent post-enrichment read filtering identified a total of 64 informative SNPs that corresponded to the expected ratio for significant polymorphisms in the parents as well as the bulks. Of these, 61 SNPs are located on potato chromosome 9 and reside within 27 individual genes, which in the sequenced potato clone DM locate to positions 45.9 to 60.9 Mb. RenSeq- and GenSeq-derived SNPs within the target region were converted into allele-specific PCR-based KASP markers and further defined the position of the resistance to a 4.3 Mb interval at the bottom end of chromosome 9 between positions 52.62-56.98 Mb.

Identification and Mapping of the Clubroot Resistance Gene CRd in Chinese Cabbage (Brassica rapa ssp. pekinensis)

The rapid spread of clubroot disease, which is caused by Plasmodiophora brassicae, threatens Brassicaceae crop production worldwide. Breeding plants that have broad-spectrum disease resistance is one of the best ways to prevent clubroot. In the present study, eight Chinese cabbage germplasms were screened using published clubroot-resistant (CR) loci-/gene-linked markers. A CR gene Crr3 potential carrier "85-74" was detected which linked to marker BRSTS61; however, "85-74" shows different responses to local pathogens "LAB-19," "LNND-2," and "LAB-10" from "CR-73" which harbors Crr3. We used a next-generation sequencing-based bulked segregant analysis approach combined with genetic mapping to detect CR genes in an F₂ segregant population generated from a cross between the Chinese cabbage inbred lines "85-74" (CR) and "BJN3-1" (clubroot susceptible). The "85-74" line showed resistance to a local pathogen "LAB-19" which was identified as race 4; a genetic analysis revealed that the resistance was conferred by a single dominant gene. The CR gene which we named CRd was mapped to a 60 kb (1 cM) region between markers yau389 and yau376 on chromosome A03. CRd is located upstream of Crr3 which was confirmed based on the physical positions of Crr3 linked markers. The identification of CRd linked markers can be applied to marker-assisted selection in the breeding of new CR cultivars of Chinese cabbage and other Brassica crops.

Identification and Mapping of the Clubroot Resistance Gene CRd in Chinese Cabbage (Brassica rapa ssp. pekinensis)

The rapid spread of clubroot disease, which is caused by Plasmodiophora brassicae, threatens Brassicaceae crop production worldwide. Breeding plants that have broad-spectrum disease resistance is one of the best ways to prevent clubroot. In the present study, eight Chinese cabbage germplasms were screened using published clubroot-resistant (CR) loci-/gene-linked markers. A CR gene Crr3 potential carrier "85-74" was detected which linked to marker BRSTS61; however, "85-74" shows different responses to local pathogens "LAB-19," "LNND-2," and "LAB-10" from "CR-73" which harbors Crr3. We used a next-generation sequencing-based bulked segregant analysis approach combined with genetic mapping to detect CR genes in an F₂ segregant population generated from a cross between the Chinese cabbage inbred lines "85-74" (CR) and "BJN3-1" (clubroot susceptible). The "85-74" line showed resistance to a local pathogen "LAB-19" which was identified as race 4; a genetic analysis revealed that the resistance was conferred by a single dominant gene. The CR gene which we named CRd was mapped to a 60 kb (1 cM) region between markers yau389 and yau376 on chromosome A03. CRd is located upstream of Crr3 which was confirmed based on the physical positions of Crr3 linked markers. The identification of CRd linked markers can be applied to marker-assisted selection in the breeding of new CR cultivars of Chinese cabbage and other Brassica crops.

Transcript-level expression control of plant NLR genes

Plant NLR genes encode sensitive immune receptors that can mediate the specific recognition of pathogen avirulence effectors and activate a strong defence response, termed effector-triggered immunity. The expression of NLRs requires strict regulation, as their ability to trigger immunity is dependent on their dose, and overexpression of NLRs results in autoimmunity and massive fitness costs. An elaborate interplay of different mechanisms controlling NLR transcript levels allows plants to maximize their defence capacity, whilst limiting negative impact on their fitness. Global suppression of NLR transcripts may be a prerequisite for the fast evolution of new NLR variants and the expansion of this gene family. Here, we summarize recent progress made towards a comprehensive understanding of NLR transcript-level expression control. Multiple mechanistic steps, including transcription as well as co-/post-transcriptional processing and transcript turn-over, contribute to balanced base levels of NLR transcripts and allow for dynamic adjustments to defence situations.

Exploring and exploiting the boundaries of host specificity using the cereal rust and mildew models

Individual plants encounter a vast number of microbes including bacteria, viruses, fungi and oomycetes through their growth cycle, yet few of these pathogens are able to infect them. Plant species have diverged over millions of years, co-evolving with few specific pathogens. The host boundaries of most pathogen species can be clearly defined. In general, the greater the genetic divergence from the preferred host, the less likely that pathogen would be able to infect that plant species. Co-evolution and divergence also occur within pathogen species, leading to highly specialized subspecies with narrow host ranges. For example, cereal rust and mildew pathogens (Puccinia and Blumeria spp.) display high host specificity as a result of ongoing co-evolution with a narrow range of grass species. In rare cases, however, some plant species are in a transition from host to nonhost or are intermediate hosts (near nonhost). Barley was reported as a useful model for genetic and molecular studies of nonhost resistance due to rare susceptibility to numerous heterologous rust and mildew fungi. This review evaluates host specificity in numerous Puccinia/Blumeria-cereal pathosystems and discusses various approaches for transferring nonhost resistance (NHR) genes between crop species to reduce the impact of important diseases in food production.

Refining a major QTL controlling spotted wilt disease resistance in cultivated peanut (Arachis hypogaea L.) and evaluating its contribution to the resistance variations in peanut germplasm

BACKGROUND

Spotted wilt, caused by tomato spotted wilt virus (TSWV), has been one of major diseases in cultivated peanut grown in the southeastern United States (US) since 1990. Previously a major quantitative trait locus (QTL) controlling spotted wilt disease resistance was mapped to an interval of 2.55 cM genetic distance corresponding to a physical distance of 14.4 Mb on chromosome A01 of peanut by using a segregating F₂ population. The current study focuses on refining this major QTL region and evaluating its contributions in the US peanut mini-core germplasm.

RESULTS

Two simple sequence repeat (SSR) markers associated with the major QTL were used to genotype F₅ individuals, and 25 heterozygous individuals were selected and developed into an F₆ segregating population. Based on visual evaluation in the field, a total of 194 susceptible F₆ individuals were selected and planted into F₇ generation for phenotyping. Nine SSR markers were used to genotype the 194 F₆ individuals, and QTL analysis revealed that a confidence interval of 15.2 Mb region had the QTL with 22.8% phenotypic variation explained (PVE). This QTL interval was further genotyped using the Amplicon-seq method. A total of 81 non-redundant single nucleotide polymorphism (SNP) and eight InDel markers were detected. No recombinant was detected among the F₆ individuals. Two InDel markers were integrated into the linkage group and helped to refine the confidence interval of this QTL into a 0.8 Mb region. To test the QTL contributes to the resistance variance in US peanut mini-core germplasm, two flanking SSR markers were used to genotype 107 mini-core germplasm accessions. No statistically significant association was observed between the genotype at the QTL region and spotted wilt resistance in the mini-core germplasm, which indicated that the resistance allelic region at this QTL didn't contribute to the resistance variance in the US peanut mini-core germplasm, thus was a unique resistance source.

CONCLUSION

A major QTL related to spotted wilt disease resistance in peanut was refined to a 0.8 Mb region on A01 chromosome, which didn't relate to spotted wilt disease resistance in the US peanut mini-core germplasm and might be a unique genetic source.

Next-generation sequencing to identify candidate genes and develop diagnostic markers for a novel Phytophthora resistance gene, RpsHC18, in soybean

KEY MESSAGE

A novel Phytophthora sojae resistance gene RpsHC18 was identified and finely mapped on soybean chromosome 3. Two NBS-LRR candidate genes were identified and two diagnostic markers of RpsHC18 were developed. Phytophthora root rot caused by Phytophthora sojae is a destructive disease of soybean. The most effective disease-control strategy is to deploy resistant cultivars carrying Phytophthora-resistant Rps genes. The soybean cultivar Huachun 18 has a broad and distinct resistance spectrum to 12 P. sojae isolates. Quantitative trait loci sequencing (QTL-seq), based on the whole-genome resequencing (WGRS) of two extreme resistant and susceptible phenotype bulks from an F2:3 population, was performed, and one 767-kb genomic region with ΔSNP-index ≥ 0.9 on chromosome 3 was identified as the RpsHC18 candidate region in Huachun 18. The candidate region was reduced to a 146-kb region by fine mapping. Nonsynonymous SNP and haplotype analyses were carried out in the 146-kb region among ten soybean genotypes using WGRS. Four specific nonsynonymous SNPs were identified in two nucleotide-binding sites-leucine-rich repeat (NBS-LRR) genes, RpsHC18-NBL1 and RpsHC18-NBL2, which were considered to be the candidate genes. Finally, one specific SNP marker in each candidate gene was successfully developed using a tetra-primer ARMS-PCR assay, and the two markers were verified to be specific for RpsHC18 and to effectively distinguish other known Rps genes. In this study, we applied an integrated genomic-based strategy combining WGRS with traditional genetic mapping to identify RpsHC18 candidate genes and develop diagnostic markers. These results suggest that next-generation sequencing is a precise, rapid and cost-effective way to identify candidate genes and develop diagnostic markers, and it can accelerate Rps gene cloning and marker-assisted selection for breeding of P. sojae-resistant soybean cultivars.

Molecular characterization of NBS-LRR genes in the soybean Rsv3 locus reveals several divergent alleles that likely confer resistance to the soybean mosaic virus

KEY MESSAGE

The divergence patterns of NBS - LRR genes in soybean Rsv3 locus were deciphered and several divergent alleles ( NBS_C, NBS_D and Columbia NBS_E ) were identified as the likely functional candidates of Rsv3. The soybean Rsv3 locus, which confers resistance to the soybean mosaic virus (SMV), has been previously mapped to a region containing five nucleotide binding site-leucine-rich repeats (NBS-LRR) genes (referred to as nbs_A-E) in Williams 82. In resistant cultivars, however, the number of NBS-LRR genes in this region and their divergence from susceptible alleles remain unclear. In the present study, we constructed and screened a bacterial artificial chromosome (BAC) library for an Rsv3-possessing cultivar, Zaoshu 18. Sequencing two positive BAC inserts on the Rsv3 locus revealed that Zaoshu 18 possesses the same gene content and order as Williams 82, but two of the NBS-LRR genes, NBS_C and NBS_D, exhibit distinct features that were not observed in the Williams 82 alleles. Obtaining these NBS-LRR genes from eight additional cultivars demonstrated that the NBS_A-D genes diverged into two different alleles: the nbs_A-D alleles were associated with the rsv3-type cultivars, whereas the NBS_A-D alleles were associated with the Rsv3-possessing cultivars. For the NBS_E gene, the cultivar Columbia possesses an allele (NBS_E) that differed from that in Zaoshu 18 and rsv3-type cultivars (nbs_E). Exchanged fragments were further detected on alleles of the NBS_C-E genes, suggesting that recombination is a major force responsible for allele divergence. Also, the LRR domains of the NBS_C-E genes exhibited extremely strong signals of positive selection. Overall, the divergence patterns of the NBS-LRR genes in Rsv3 locus elucidated by this study indicate that not only NBS_C but also NBS_D and Columbia NBS_E are likely functional alleles that confer resistance to SMV.

Genotyping by Sequencing Highlights a Polygenic Resistance to Ralstonia pseudosolanacearum in Eggplant (Solanum melongena L.)

Eggplant cultivation is limited by numerous diseases, including the devastating bacterial wilt (BW) caused by the Ralstonia solanacearum species complex (RSSC). Within the RSSC, Ralstonia pseudosolanacearum (including phylotypes I and III) causes severe damage to all solanaceous crops, including eggplant. Therefore, the creation of cultivars resistant to R. pseudosolanacearum strains is a major goal for breeders. An intraspecific eggplant population, segregating for resistance, was created from the cross between the susceptible MM738 and the resistant EG203 lines. The population of 123 doubled haploid lines was challenged with two strains belonging to phylotypes I (PSS4) and III (R3598), which both bypass the published EBWR9 BW-resistance quantitative trait locus (QTL). Ten and three QTLs of resistance to PSS4 and to R3598, respectively, were detected and mapped. All were strongly influenced by environmental conditions. The most stable QTLs were found on chromosomes 3 and 6. Given their estimated physical position, these newly detected QTLs are putatively syntenic with BW-resistance QTLs in tomato. In particular, the QTLs' position on chromosome 6 overlaps with that of the major broad-spectrum tomato resistance QTL Bwr-6. The present study is a first step towards understanding the complex polygenic system, which underlies the high level of BW resistance of the EG203 line.

Diversification of defensins and NLRs in Arabidopsis species by different evolutionary mechanisms

BACKGROUND

Genes encoding proteins underlying host-pathogen co-evolution and which are selected for new resistance specificities frequently are under positive selection, a process that maintains diversity. Here, we tested the contribution of natural selection, recombination and transcriptional divergence to the evolutionary diversification of the plant defensins superfamily in three Arabidopsis species. The intracellular NOD-like receptor (NLR) family was used for comparison because positive selection has been well documented in its members. Similar to defensins, NLRs are encoded by a large and polymorphic gene family and many of their members are involved in the immune response.

RESULTS

Gene trees of Arabidopsis defensins (DEFLs) show a high prevalence of clades containing orthologs. This indicates that their diversity dates back to a common ancestor and species-specific duplications did not significantly contribute to gene family expansion. DEFLs are characterized by a pervasive pattern of neutral evolution with infrequent positive and negative selection as well as recombination. In comparison, most NLR alignment groups are characterized by frequent occurrence of positive selection and recombination in their leucine-rich repeat (LRR) domain as well negative selection in their nucleotide-binding (NB-ARC) domain. While major NLR subgroups are expressed in pistils and leaves both in presence or absence of pathogen infection, the members of DEFL alignment groups are predominantly transcribed in pistils. Furthermore, conserved groups of NLRs and DEFLs are differentially expressed in response to Fusarium graminearum regardless of whether these genes are under positive selection or not.

CONCLUSIONS

The present analyses of NLRs expands previous studies in Arabidopsis thaliana and highlights contrasting patterns of purifying and diversifying selection affecting different gene regions. DEFL genes show a different evolutionary trend, with fewer recombination events and significantly fewer instances of natural selection. Their heterogeneous expression pattern suggests that transcriptional divergence probably made the major contribution to functional diversification. In comparison to smaller families encoding pathogenesis-related (PR) proteins under positive selection, DEFLs are involved in a wide variety of processes that altogether might pose structural and functional trade-offs to their family-wide pattern of evolution.

Modification of meiotic recombination by natural variation in plants

Meiosis is a specialized cell division that produces haploid gametes required for sexual reproduction. During the first meiotic division, homologous chromosomes pair and undergo reciprocal crossing over, which recombines linked sequence variation. Meiotic recombination frequency varies extensively both within and between species. In this review, we will examine the molecular basis of meiotic recombination rate variation, with an emphasis on plant genomes. We first consider cis modification caused by polymorphisms at the site of recombination, or elsewhere on the same chromosome. We review cis effects caused by mismatches within recombining joint molecules, the effect of structural hemizygosity, and the role of specific DNA sequence motifs. In contrast, trans modification of recombination is exerted by polymorphic loci encoding diffusible molecules, which are able to modulate recombination on the same and/or other chromosomes. We consider trans modifiers that act to change total recombination levels, hotspot locations, or interactions between homologous and homeologous chromosomes in polyploid species. Finally, we consider the significance of genetic variation that modifies meiotic recombination for adaptation and evolution of plant species.

Mapping Loci That Control Tuber and Foliar Symptoms Caused by PVY in Autotetraploid Potato (Solanum tuberosum L.)

Potato tuber necrotic ringspot disease (PTNRD) is a tuber deformity associated with infection by the tuber necrotic strain of Potato virus Y (PVYNTN). PTNRD negatively impacts tuber quality and marketability, and poses a serious threat to seed and commercial potato production worldwide. PVYNTN symptoms differ in the cultivars Waneta and Pike: Waneta expresses severe PTNRD and foliar mosaic with vein and leaf necrosis, whereas Pike does not express PTNRD and mosaic is the only foliar symptom. To map loci that influence tuber and foliar symptoms, 236 F1 progeny of a cross between Waneta and Pike were inoculated with PVYNTN isolate NY090029 and genotyped using 12,808 potato SNPs. Foliar symptom type and severity were monitored for 10 wk, while tubers were evaluated for PTNRD expression at harvest and again after 60 d in storage. Pairwise correlation analyses indicate a strong association between PTNRD and vein necrosis (τ = 0.4195). QTL analyses revealed major-effect QTL on chromosomes 4 and 5 for mosaic, 4 for PTNRD, and 5 for foliar necrosis symptoms. Locating QTL associated with PVY-related symptoms provides a foundation for breeders to develop markers that can be used to eliminate potato clones with undesirable phenotypes, e.g., those likely to develop PTNRD or to be symptomless carriers of PVY.

Fine mapping of Pi57(t) conferring broad spectrum resistance against Magnaporthe oryzae in introgression line IL-E1454 derived from Oryza longistaminata

Wild species of the genus Oryza are excellent gene pools for improvement of agronomic traits of Asian cultivated rice. The blast resistance gene Pi57(t) in the introgression line IL-E1454 derived from Oryza longistaminata was previously mapped on rice chromosome 12. Inoculation with 322 Magnaporthe oryzae isolates collected from 6 countries indicated that Pi57(t) conferred broad spectrum resistance against M. oryzae. Two mapping populations consisting of 29070 and 10375 F2 plants derived from the crosses of resistant donor IL-E1454 with susceptible parents RD23 and Lijiangxintuanheigu respectively, were used for fine mapping of Pi57(t) locus. Based on genotyping and phenotyping results of recombinants screened from the two crosses, Pi57(t) was finally mapped to a 51.7-kb region flanked by two molecular markers (STS57-320 and STS57-372) on the short arm and close to the centromere of chromosome 12. Six candidate resistance genes were predicted in the target region according to the reference sequence of Nipponbare. These results could facilitate both marker-assisted selection for disease-resistant breeding and gene cloning.

A New Resistance Gene against Potato Late Blight Originating from Solanum pinnatisectum Located on Potato Chromosome 7

Late blight, caused by the pathogen Phytophthora infestans, is one of the most devastating diseases of potato. Here, we describe a new single dominant resistance gene, Rpi2, from the Mexican diploid wild species Solanum pinnatisectum that confers high level and broad spectrum resistance to late blight. The Rpi2 locus confers full resistance to complex isolates of P. infestans, for which race specificity has not yet been demonstrated. This new gene, flanked by the RFLP-derived marker SpT1756 and AFLP-derived marker SpAFLP2 with a minimal genetic distance of 0.8 cM, was mapped to potato chromosome 7. Using the genomic sequence data of potato, we estimated that the physical distance of the nearest marker to the resistance gene was about 27 kb. The map location and other evidence indicated that this resistance locus was different from the previously reported resistance locus Rpi1 on the same chromosome from S. pinnatisectum. The presence of other reported resistance genes in the target region, such as Gro1-4, I-3, and three NBS-LLR like genes, on a homologous tomato genome segment indicates the Rpi2-related region is a hotspot for resistance genes. Comparative sequence analysis showed that the order of nine markers mapped to the Rpi2 genetic map was highly conserved on tomato chromosome 7; however, some rearrangements were observed in the potato genome sequence. Additional markers and potential resistance genes will promote accurate location of the site of Rpi2 and help breeders to introduce this resistance gene into different cultivars by marker-aided selection.

Comparative Genomics of Non-TNL Disease Resistance Genes from Six Plant Species

Disease resistance genes (R genes), as part of the plant defense system, have coevolved with corresponding pathogen molecules. The main objectives of this project were to identify non-Toll interleukin receptor, nucleotide-binding site, leucine-rich repeat (nTNL) genes and elucidate their evolutionary divergence across six plant genomes. Using reference sequences from Arabidopsis, we investigated nTNL orthologs in the genomes of common bean, Medicago, soybean, poplar, and rice. We used Hidden Markov Models for sequence identification, performed model-based phylogenetic analyses, visualized chromosomal positioning, inferred gene clustering, and assessed gene expression profiles. We analyzed 908 nTNL R genes in the genomes of the six plant species, and classified them into 12 subgroups based on the presence of coiled-coil (CC), nucleotide binding site (NBS), leucine rich repeat (LRR), resistance to Powdery mildew 8 (RPW8), and BED type zinc finger domains. Traditionally classified CC-NBS-LRR (CNL) genes were nested into four clades (CNL A-D) often with abundant, well-supported homogeneous subclades of Type-II R genes. CNL-D members were absent in rice, indicating a unique R gene retention pattern in the rice genome. Genomes from Arabidopsis, common bean, poplar and soybean had one chromosome without any CNL R genes. Medicago and Arabidopsis had the highest and lowest number of gene clusters, respectively. Gene expression analyses suggested unique patterns of expression for each of the CNL clades. Differential gene expression patterns of the nTNL genes were often found to correlate with number of introns and GC content, suggesting structural and functional divergence.

Plant STAND P-loop NTPases: a current perspective of genome distribution, evolution, and function : Plant STAND P-loop NTPases: genomic organization, evolution, and molecular mechanism models contribute broadly to plant pathogen defense

STAND P-loop NTPase is the common weapon used by plant and other organisms from all three kingdoms of life to defend themselves against pathogen invasion. The purpose of this study is to review comprehensively the latest finding of plant STAND P-loop NTPase related to their genomic distribution, evolution, and their mechanism of action. Earlier, the plant STAND P-loop NTPase known to be comprised of only NBS-LRRs/AP-ATPase/NB-ARC ATPase. However, recent finding suggests that genome of early green plants comprised of two types of STAND P-loop NTPases: (1) mammalian NACHT NTPases and (2) NBS-LRRs. Moreover, YchF (unconventional G protein and members of P-loop NTPase) subfamily has been reported to be exceptionally involved in biotic stress (in case of Oryza sativa), thereby a novel member of STAND P-loop NTPase in green plants. The lineage-specific expansion and genome duplication events are responsible for abundance of plant STAND P-loop NTPases; where "moderate tandem and low segmental duplication" trajectory followed in majority of plant species with few exception (equal contribution of tandem and segmental duplication). Since the past decades, systematic research is being investigated into NBS-LRR function supported the direct recognition of pathogen or pathogen effectors by the latest models proposed via 'integrated decoy' or 'sensor domains' model. Here, we integrate the recently published findings together with the previous literature on the genomic distribution, evolution, and distinct models proposed for functional molecular mechanism of plant STAND P-loop NTPases.

Evidence-based gene models for structural and functional annotations of the oil palm genome

Background

Oil palm is an important source of edible oil. The importance of the crop, as well as its long breeding cycle (10-12 years) has led to the sequencing of its genome in 2013 to pave the way for genomics-guided breeding. Nevertheless, the first set of gene predictions, although useful, had many fragmented genes. Classification and characterization of genes associated with traits of interest, such as those for fatty acid biosynthesis and disease resistance, were also limited. Lipid-, especially fatty acid (FA)-related genes are of particular interest for the oil palm as they specify oil yields and quality. This paper presents the characterization of the oil palm genome using different gene prediction methods and comparative genomics analysis, identification of FA biosynthesis and disease resistance genes, and the development of an annotation database and bioinformatics tools.

Results

Using two independent gene-prediction pipelines, Fgenesh++ and Seqping, 26,059 oil palm genes with transcriptome and RefSeq support were identified from the oil palm genome. These coding regions of the genome have a characteristic broad distribution of GC₃ (fraction of cytosine and guanine in the third position of a codon) with over half the GC₃-rich genes (GC₃ ≥ 0.75286) being intronless. In comparison, only one-seventh of the oil palm genes identified are intronless. Using comparative genomics analysis, characterization of conserved domains and active sites, and expression analysis, 42 key genes involved in FA biosynthesis in oil palm were identified. For three of them, namely EgFABF, EgFABH and EgFAD3, segmental duplication events were detected. Our analysis also identified 210 candidate resistance genes in six classes, grouped by their protein domain structures.

Conclusions

We present an accurate and comprehensive annotation of the oil palm genome, focusing on analysis of important categories of genes (GC₃-rich and intronless), as well as those associated with important functions, such as FA biosynthesis and disease resistance. The study demonstrated the advantages of having an integrated approach to gene prediction and developed a computational framework for combining multiple genome annotations. These results, available in the oil palm annotation database (http://palmxplore.mpob.gov.my), will provide important resources for studies on the genomes of oil palm and related crops.

Reviewers

This article was reviewed by Alexander Kel, Igor Rogozin, and Vladimir A. Kuznetsov.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-017-0191-4) contains supplementary material, which is available to authorized users.

Heterogeneous gene duplications can be adaptive because they permanently associate overdominant alleles

Gene duplications are widespread in genomes, but their role in contemporary adaptation is not fully understood. Although mostly deleterious, homogeneous duplications that associate identical repeats of a locus often increase the quantity of protein produced, which can be selected in certain environments. However, another type exists: heterogeneous gene duplications, which permanently associate two (or more) alleles of a single locus on the same chromosome. They are far less studied, as only few examples of contemporary heterogeneous duplications are known. Haldane proposed in 1954 that they could be adaptive in situations of heterozygote advantage, or overdominance, but this hypothesis was never tested. To assess its validity, we took advantage of the well-known model of insecticide resistance in mosquitoes. We used experimental evolution to estimate the fitnesses associated with homozygous and heterozygous genotypes in different selection regimes. It first showed that balanced antagonist selective pressures frequently induce overdominance, generating stable polymorphic equilibriums. The frequency of equilibrium moreover depends on the magnitude of two antagonistic selective pressures, the survival advantage conferred by the resistant allele versus the selective costs it induces. We then showed that heterogeneous duplications are selected over single-copy alleles in such contexts. They allow the fixation of the heterozygote phenotype, providing an alternative and stable intermediate fitness trade-off. By allowing the rapid fixation of divergent alleles, this immediate advantage could contribute to the rarity of overdominance. More importantly, it also creates new material for long-term genetic innovation, making a crucial but underestimated contribution to the evolution of new genes and gene families.

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 History of Subtilases in Land Plants and Their Involvement in Symbiotic Interactions

Subtilases, a family of proteases involved in a variety of developmental processes in land plants, are also involved in both mutualistic symbiosis and host-pathogen interactions in different angiosperm lineages. We examined the evolutionary history of subtilase genes across land plants through a phylogenetic analysis integrating amino acid sequence data from full genomes, transcriptomes, and characterized subtilases of 341 species of diverse green algae and land plants along with subtilases from 12 species of other eukaryotes, archaea, and bacteria. Our analysis reconstructs the subtilase gene phylogeny and identifies 11 new gene lineages, six of which have no previously characterized members. Two large, previously unnamed, subtilase gene lineages that diverged before the origin of angiosperms accounted for the majority of subtilases shown to be associated with symbiotic interactions. These lineages expanded through both whole-genome and tandem duplication, with differential neofunctionalization and subfunctionalization creating paralogs associated with different symbioses, including nodulation with nitrogen-fixing bacteria, arbuscular mycorrhizae, and pathogenesis in different plant clades. This study demonstrates for the first time that a key gene family involved in plant-microbe interactions proliferated in size and functional diversity before the explosive radiation of angiosperms.

Distinct Patterns of Gene Gain and Loss: Diverse Evolutionary Modes of NBS-Encoding Genes in Three Solanaceae Crop Species

Plant resistance conferred by nucleotide binding site (NBS)-encoding resistance genes plays a key role in the defense against various pathogens throughout the entire plant life cycle. However, comparative analyses for the systematic evaluation and determination of the evolutionary modes of NBS-encoding genes among Solanaceae species are rare. In this study, 447, 255, and 306 NBS-encoding genes were identified from the genomes of potato, tomato, and pepper, respectively. These genes usually clustered as tandem arrays on chromosomes; few existed as singletons. Phylogenetic analysis indicated that three subclasses [TNLs (TIR-NBS-LRR), CNLs (CC-NBS-LRR), and RNLs (RPW8-NBS-LRR)] each formed a monophyletic clade and were distinguished by unique exon/intron structures and amino acid motif sequences. By comparing phylogenetic and systematic relationships, we inferred that the NBS-encoding genes in the present genomes of potato, tomato, and pepper were derived from 150 CNL, 22 TNL, and 4 RNL ancestral genes, and underwent independent gene loss and duplication events after speciation. The NBS-encoding genes therefore exhibit diverse and dynamic evolutionary patterns in the three Solanaceae species, giving rise to the discrepant gene numbers observed today. Potato shows a “consistent expansion” pattern, tomato exhibits a pattern of “first expansion and then contraction,” and pepper presents a “shrinking” pattern. The earlier expansion of CNLs in the common ancestor led to the dominance of this subclass in gene numbers. However, RNLs remained at low copy numbers due to their specific functions. Along the evolutionary process of NBS-encoding genes in Solanaceae, species-specific tandem duplications contributed the most to gene expansions.

Resolution of deep eudicot phylogeny and their temporal diversification using nuclear genes from transcriptomic and genomic datasets

Explosive diversification is widespread in eukaryotes, making it difficult to resolve phylogenetic relationships. Eudicots contain c. 75% of extant flowering plants, are important for human livelihood and terrestrial ecosystems, and have probably experienced explosive diversifications. The eudicot phylogenetic relationships, especially among those of the Pentapetalae, remain unresolved. Here, we present a highly supported eudicot phylogeny and diversification rate shifts using 31 newly generated transcriptomes and 88 other datasets covering 70% of eudicot orders. A highly supported eudicot phylogeny divided Pentapetalae into two groups: one with rosids, Saxifragales, Vitales and Santalales; the other containing asterids, Caryophyllales and Dilleniaceae, with uncertainty for Berberidopsidales. Molecular clock analysis estimated that crown eudicots originated c. 146 Ma, considerably earlier than earliest tricolpate pollen fossils and most other molecular clock estimates, and Pentapetalae sequentially diverged into eight major lineages within c. 15 Myr. Two identified increases of diversification rate are located in the stems leading to Pentapetalae and asterids, and lagged behind the gamma hexaploidization. The nuclear genes from newly generated transcriptomes revealed a well-resolved eudicot phylogeny, sequential separation of major core eudicot lineages and temporal mode of diversifications, providing new insights into the evolutionary trend of morphologies and contributions to the diversification of eudicots.

Genome sequence of Plasmopara viticola and insight into the pathogenic mechanism

Plasmopara viticola causes downy mildew disease of grapevine which is one of the most devastating diseases of viticulture worldwide. Here we report a 101.3 Mb whole genome sequence of P. viticola isolate 'JL-7-2' obtained by a combination of Illumina and PacBio sequencing technologies. The P. viticola genome contains 17,014 putative protein-coding genes and has ~26% repetitive sequences. A total of 1,301 putative secreted proteins, including 100 putative RXLR effectors and 90 CRN effectors were identified in this genome. In the secretome, 261 potential pathogenicity genes and 95 carbohydrate-active enzymes were predicted. Transcriptional analysis revealed that most of the RXLR effectors, pathogenicity genes and carbohydrate-active enzymes were significantly up-regulated during infection. Comparative genomic analysis revealed that P. viticola evolved independently from the Arabidopsis downy mildew pathogen Hyaloperonospora arabidopsidis. The availability of the P. viticola genome provides a valuable resource not only for comparative genomic analysis and evolutionary studies among oomycetes, but also enhance our knowledge on the mechanism of interactions between this biotrophic pathogen and its host.

The causes and consequences of DNA methylome variation in plants

Epigenetic variation - polymorphisms at the level of DNA methylation or histone modifications - modulates chromatin accessibility, which can perturb transcriptional activity and spur phenotypic variation. Determining the origin, frequency spectrum, and consequences of epigenetic variants is key to understanding the role of this variation in generating stable phenotypic variation in plants. Here we review recent literature on DNA methylation variation in both model and crop plant species with a focus on the link between genotype, epigenotype, and transcription. We highlight population epigenomics studies that explore the relationship between epigenetic variants and genetic diversity. Moreover, we provide an overview of relevant studies that together advocate a minor, albeit significant role for epigenetic variation in directing specific transcriptional changes.

Exploring structural variation and gene family architecture with De Novo assemblies of 15 Medicago genomes

BACKGROUND

Previous studies exploring sequence variation in the model legume, Medicago truncatula, relied on mapping short reads to a single reference. However, read-mapping approaches are inadequate to examine large, diverse gene families or to probe variation in repeat-rich or highly divergent genome regions. De novo sequencing and assembly of M. truncatula genomes enables near-comprehensive discovery of structural variants (SVs), analysis of rapidly evolving gene families, and ultimately, construction of a pan-genome.

RESULTS

Genome-wide synteny based on 15 de novo M. truncatula assemblies effectively detected different types of SVs indicating that as much as 22% of the genome is involved in large structural changes, altogether affecting 28% of gene models. A total of 63 million base pairs (Mbp) of novel sequence was discovered, expanding the reference genome space for Medicago by 16%. Pan-genome analysis revealed that 42% (180 Mbp) of genomic sequences is missing in one or more accession, while examination of de novo annotated genes identified 67% (50,700) of all ortholog groups as dispensable - estimates comparable to recent studies in rice, maize and soybean. Rapidly evolving gene families typically associated with biotic interactions and stress response were found to be enriched in the accession-specific gene pool. The nucleotide-binding site leucine-rich repeat (NBS-LRR) family, in particular, harbors the highest level of nucleotide diversity, large effect single nucleotide change, protein diversity, and presence/absence variation. However, the leucine-rich repeat (LRR) and heat shock gene families are disproportionately affected by large effect single nucleotide changes and even higher levels of copy number variation.

CONCLUSIONS

Analysis of multiple M. truncatula genomes illustrates the value of de novo assemblies to discover and describe structural variation, something that is often under-estimated when using read-mapping approaches. Comparisons among the de novo assemblies also indicate that different large gene families differ in the architecture of their structural variation.

Draft Genome Sequence of Mentha longifolia and Development of Resources for Mint Cultivar Improvement

The genus Mentha encompasses mint species cultivated for their essential oils, which are formulated into a vast array of consumer products. Desirable oil characteristics and resistance to the fungal disease Verticillium wilt are top priorities for the mint industry. However, cultivated mints have complex polyploid genomes and are sterile. Breeding efforts, therefore, require the development of genomic resources for fertile mint species. Here, we present draft de novo genome and plastome assemblies for a wilt-resistant South African accession of Mentha longifolia (L.) Huds., a diploid species ancestral to cultivated peppermint and spearmint. The 353 Mb genome contains 35 597 predicted protein-coding genes, including 292 disease resistance gene homologs, and nine genes determining essential oil characteristics. A genetic linkage map ordered 1397 genome scaffolds on 12 pseudochromosomes. More than two million simple sequence repeats were identified, which will facilitate molecular marker development. The M. longifolia genome is a valuable resource for both metabolic engineering and molecular breeding. This is exemplified by employing the genome sequence to clone and functionally characterize the promoters in a peppermint cultivar, and demonstrating the utility of a glandular trichome-specific promoter to increase expression of a biosynthetic gene, thereby modulating essential oil composition.

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.

Discovery and introgression of the wild sunflower-derived novel downy mildew resistance gene Pl 19 in confection sunflower (Helianthus annuus L.)

A new downy mildew resistance gene, Pl ₁₉ , was identified from wild Helianthus annuus accession PI 435414, introduced to confection sunflower, and genetically mapped to linkage group 4 of the sunflower genome. Wild Helianthus annuus accession PI 435414 exhibited resistance to downy mildew, which is one of the most destructive diseases to sunflower production globally. Evaluation of the 140 BC₁F_2:3 families derived from the cross of CMS CONFSCLB1 and PI 435414 against Plasmopara halstedii race 734 revealed that a single dominant gene controls downy mildew resistance in the population. Bulked segregant analysis conducted in the BC₁F₂ population with 860 simple sequence repeat (SSR) markers indicated that the resistance derived from wild H. annuus was associated with SSR markers located on linkage group (LG) 4 of the sunflower genome. To map and tag this resistance locus, designated Pl ₁₉ , 140 BC₁F₂ individuals were used to construct a linkage map of the gene region. Two SSR markers, ORS963 and HT298, were linked to Pl ₁₉ within a distance of 4.7 cM. After screening 27 additional single nucleotide polymorphism (SNP) markers previously mapped to this region, two flanking SNP markers, NSA_003564 and NSA_006089, were identified as surrounding the Pl ₁₉ gene at a distance of 0.6 cM from each side. Genetic analysis indicated that Pl ₁₉ is different from Pl ₁₇ , which had previously been mapped to LG4, but is closely linked to Pl ₁₇ . This new gene is highly effective against the most predominant and virulent races of P. halstedii currently identified in North America and is the first downy mildew resistance gene that has been transferred to confection sunflower. The selected resistant germplasm derived from homozygous BC₂F₃ progeny provides a novel gene for use in confection sunflower breeding programs.

Analysis of Plant Pan-Genomes and Transcriptomes with GET_HOMOLOGUES-EST, a Clustering Solution for Sequences of the Same Species

The pan-genome of a species is defined as the union of all the genes and non-coding sequences found in all its individuals. However, constructing a pan-genome for plants with large genomes is daunting both in sequencing cost and the scale of the required computational analysis. A more affordable alternative is to focus on the genic repertoire by using transcriptomic data. Here, the software GET_HOMOLOGUES-EST was benchmarked with genomic and RNA-seq data of 19 Arabidopsis thaliana ecotypes and then applied to the analysis of transcripts from 16 Hordeum vulgare genotypes. The goal was to sample their pan-genomes and classify sequences as core, if detected in all accessions, or accessory, when absent in some of them. The resulting sequence clusters were used to simulate pan-genome growth, and to compile Average Nucleotide Identity matrices that summarize intra-species variation. Although transcripts were found to under-estimate pan-genome size by at least 10%, we concluded that clusters of expressed sequences can recapitulate phylogeny and reproduce two properties observed in A. thaliana gene models: accessory loci show lower expression and higher non-synonymous substitution rates than core genes. Finally, accessory sequences were observed to preferentially encode transposon components in both species, plus disease resistance genes in cultivated barleys, and a variety of protein domains from other families that appear frequently associated with presence/absence variation in the literature. These results demonstrate that pan-genome analyses are useful to explore germplasm diversity.

Complex Evolutionary Dynamics of Massively Expanded Chemosensory Receptor Families in an Extreme Generalist Chelicerate Herbivore

While mechanisms to detoxify plant produced, anti-herbivore compounds have been associated with plant host use by herbivores, less is known about the role of chemosensory perception in their life histories. This is especially true for generalists, including chelicerate herbivores that evolved herbivory independently from the more studied insect lineages. To shed light on chemosensory perception in a generalist herbivore, we characterized the chemosensory receptors (CRs) of the chelicerate two-spotted spider mite, Tetranychus urticae, an extreme generalist. Strikingly, T. urticae has more CRs than reported in any other arthropod to date. Including pseudogenes, 689 gustatory receptors were identified, as were 136 degenerin/Epithelial Na+ Channels (ENaCs) that have also been implicated as CRs in insects. The genomic distribution of T. urticae gustatory receptors indicates recurring bursts of lineage-specific proliferations, with the extent of receptor clusters reminiscent of those observed in the CR-rich genomes of vertebrates or C. elegans Although pseudogenization of many gustatory receptors within clusters suggests relaxed selection, a subset of receptors is expressed. Consistent with functions as CRs, the genomic distribution and expression of ENaCs in lineage-specific T. urticae expansions mirrors that observed for gustatory receptors. The expansion of ENaCs in T. urticae to > 3-fold that reported in other animals was unexpected, raising the possibility that ENaCs in T. urticae have been co-opted to fulfill a major role performed by unrelated CRs in other animals. More broadly, our findings suggest an elaborate role for chemosensory perception in generalist herbivores that are of key ecological and agricultural importance.

Dissimilar evolutionary histories of two resistance gene families in the genus Solanum

Genomic analyses have shown that most genes in eukaryotic lineages belong to families. Gene families vary in terms of number of members, nucleotide similarity, gene integrity, expression, and function. Often, the members of gene families are arranged in clusters, which contribute to maintaining similarity among gene copies and also to generate duplicates through replication errors. Gene families offer us an opportunity to examine the forces involved in the evolution of the genomes and to study recombination events and genomic rearrangements. In this work, we focused on the evolution of two plant resistance gene families, Sw5 and Mi-1, and analyzed the completely sequenced nuclear genomes of potato and tomato. We first noticed that the potato genome carries larger resistance gene families than tomato, but all gene copies are pseudogenes. Second, phylogenetic analyses indicated that Sw5 and Mi-1 gene families had dissimilar evolutionary histories. In contrast to Sw5, Mi-1 homologues suffered repeated gene conversion events among the gene copies, particularly in the tomato genome.