Species-specific duplications driving the recent expansion of NBS-LRR genes in five Rosaceae species

Chromosome-level genome of Ambrosia trifida provides insights into adaptation and the evolution of pollen allergens

Ambrosia trifida (giant ragweed) is an invasive plant that can cause serious damage to natural ecosystems and severe respiratory allergies. However, the genomic basis of invasive adaptation and pollen allergens in Ambrosia species remain largely unknown. Here, we present a 1.66 Gb chromosome-scale reference genome for giant ragweed and identified multiple types of genome duplications, which are responsible for its rapid environmental adaptation and pollen development. The largest copies number and species-specific expansions of resistance-related gene families compared to Heliantheae alliance might contribute to resist stresses, pathogens and rapid adaptation. To extend the knowledge of evolutionary process of allergic pollen proteins, we predicted 26 and 168 potential pollen allergen candidates for giant ragweed and other Asteraceae plant species by combining machine learning and identity screening. Interestingly, we observed a specific tandemly repeated array for potential allergenic pectate lyases among Ambrosia species. Rapid evolutionary rates on putative pectate lyase allergens may imply a crucial role of nonsynonymous mutations on amino acid residues for plant biological function and allergenicity. Altogether, this study provides insight into the molecular ecological adaptation and putative pollen allergens prediction that will be helpful in promoting invasion genomic research and evolution of putative pollen allergy in giant ragweed.

The devastating oomycete phytopathogen Phytophthora cactorum: Insights into its biology and molecular features

Phytophthora cactorum is one of the most economically important soilborne oomycete pathogens in the world. It infects more than 200 plant species spanning 54 families, most of which are herbaceous and woody species. Although traditionally considered to be a generalist, marked differences of P. cactorum isolates occur in degree of pathogenicity to different hosts. As the impact of crop loss caused by this species has increased recently, there has been a tremendous increase in the development of new tools, resources, and management strategies to study and combat this devastating pathogen. This review aims to integrate recent molecular biology analyses of P. cactorum with the current knowledge of the cellular and genetic basis of its growth, development, and host infection. The goal is to provide a framework for further studies of P. cactorum by highlighting important biological and molecular features, shedding light on the functions of pathogenicity factors, and developing effective control measures.

TAXONOMY

P. cactorum (Leb. & Cohn) Schröeter: kingdom Chromista; phylum Oomycota; class Oomycetes; order Peronosporales; family Peronosporaceae; genus Phytophthora.

HOST RANGE

Infects about 200 plant species in 154 genera representing 54 families. Economically important host plants include strawberry, apple, pear, Panax spp., and walnut.

DISEASE SYMPTOMS

The soilborne pathogen often causes root, stem, collar, crown, and fruit rots, as well as foliar infection, stem canker, and seedling damping off.

Pan-genome and transcriptome analyses provide insights into genomic variation and differential gene expression profiles related to disease resistance and fatty acid biosynthesis in eastern black walnut (Juglans nigra)

Walnut (Juglans) species are used as nut crops worldwide. Eastern black walnut (EBW, Juglans nigra), a diploid, horticultural important woody species is native to much of eastern North America. Although it is highly valued for its wood and nut, there are few resources for understanding EBW genetics. Here, we present a high-quality genome assembly of J. nigra based on Illumina, Pacbio, and Hi-C technologies. The genome size was 540.8 Mb, with a scaffold N50 size of 35.1 Mb, and 99.0% of the assembly was anchored to 16 chromosomes. Using this genome as a reference, the resequencing of 74 accessions revealed the effective population size of J. nigra declined during the glacial maximum. A single whole-genome duplication event was identified in the J. nigra genome. Large syntenic blocks among J. nigra, Juglans regia, and Juglans microcarpa predominated, but inversions of more than 600 kb were identified. By comparing the EBW genome with those of J. regia and J. microcarpa, we detected InDel sizes of 34.9 Mb in J. regia and 18.3 Mb in J. microcarpa, respectively. Transcriptomic analysis of differentially expressed genes identified five presumed NBS-LRR (NUCLEOTIDE BINDING SITE-LEUCINE-RICH REPEAT) genes were upregulated during the development of walnut husks and shells compared to developing embryos. We also identified candidate genes with essential roles in seed oil synthesis, including FAD (FATTY ACID DESATURASE) and OLE (OLEOSIN). Our work advances the understanding of fatty acid bioaccumulation and disease resistance in nut crops, and also provides an essential resource for conducting genomics-enabled breeding in walnut.

Evolution and functional analysis of the GRAS family genes in six Rosaceae species

BACKGROUND

GRAS genes formed one of the important transcription factor gene families in plants, had been identified in several plant species. The family genes were involved in plant growth, development, and stress resistance. However, the comparative analysis of GRAS genes in Rosaceae species was insufficient.

RESULTS

In this study, a total of 333 GRAS genes were identified in six Rosaceae species, including 51 in strawberry (Fragaria vesca), 78 in apple (Malus domestica), 41 in black raspberry (Rubus occidentalis), 59 in European pear (Pyrus communis), 56 in Chinese rose (Rosa chinensis), and 48 in peach (Prunus persica). Motif analysis showed the VHIID domain, SAW motif, LR I region, and PFYRE motif were considerably conserved in the six Rosaceae species. All GRAS genes were divided into 10 subgroups according to phylogenetic analysis. A total of 15 species-specific duplicated clades and 3 lineage-specific duplicated clades were identified in six Rosaceae species. Chromosomal localization presented the uneven distribution of GRAS genes in six Rosaceae species. Duplication events contributed to the expression of the GRAS genes, and Ka/Ks analysis suggested the purification selection as a major force during the evolution process in six Rosaceae species. Cis-acting elements and GO analysis revealed that most of the GRAS genes were associated with various environmental stress in six Rosaceae species. Coexpression network analysis showed the mutual regulatory relationship between GRAS and bZIP genes, suggesting the ability of the GRAS gene to regulate abiotic stress in woodland strawberry. The expression pattern elucidated the transcriptional levels of FvGRAS genes in various tissues and the drought and salt stress in woodland strawberry, which were verified by RT-qPCR analysis.

CONCLUSIONS

The evolution and functional analysis of GRAS genes provided insights into the further understanding of GRAS genes on the abiotic stress of Rosaceae species.

Genome-wide analysis of NBS-LRR genes in Rosaceae species reveals distinct evolutionary patterns

The nucleotide-binding site and leucine-rich repeat (NBS-LRR) genes, one of the largest gene families in plants, are evolving rapidly and playing a critical role in plant resistance to pathogens. In this study, a genome-wide search in 12 Rosaceae genomes screened out 2188 NBS-LRR genes, with the gene number varied distinctively across different species. The reconciled phylogeny revealed 102 ancestral genes (7 RNLs, 26 TNLs, and 69 CNLs), which underwent independent gene duplication and loss events during the divergence of the Rosaceae. The NBS-LRR genes exhibited dynamic and distinct evolutionary patterns in the 12 Rosaceae species due to independent gene duplication/loss events, which resulted the discrepancy of NBS-LRR gene number among Rosaceae species. Specifically, Rubus occidentalis, Potentilla micrantha, Fragaria iinumae and Gillenia trifoliata, displayed a “first expansion and then contraction” evolutionary pattern; Rosa chinensis exhibited a “continuous expansion” pattern; F. vesca had a “expansion followed by contraction, then a further expansion” pattern, three Prunus species and three Maleae species shared a “early sharp expanding to abrupt shrinking” pattern. Overall, this study elucidated the dynamic and complex evolutionary patterns of NBS-LRR genes in the 12 Rosaceae species, and could assist further investigation of mechanisms driving these evolutionary patterns.

Chromosome-level genome assembly of Mentha longifolia L. reveals gene organization underlying disease resistance and essential oil traits

Mentha longifolia (L.) Huds., a wild, diploid mint species, has been developed as a model for mint genetic and genomic research to aid breeding efforts that target Verticillium wilt disease resistance and essential oil monoterpene composition. Here, we present a near-complete, chromosome-scale mint genome assembly for M. longifolia USDA accession CMEN 585. This new assembly is an update of a previously published genome draft, with dramatic improvements. A total of 42,107 protein-coding genes were annotated and placed on 12 chromosomal scaffolds. One hundred fifty-three genes contained conserved sequence domains consistent with nucleotide binding site-leucine-rich-repeat plant disease resistance genes. Homologs of genes implicated in Verticillium wilt resistance in other plant species were also identified. Multiple paralogs of genes putatively involved in p-menthane monoterpenoid biosynthesis were identified and several cases of gene clustering documented. Heterologous expression of candidate genes, purification of recombinant target proteins, and subsequent enzyme assays allowed us to identify the genes underlying the pathway that leads to the most abundant monoterpenoid volatiles. The bioinformatic and functional analyses presented here are laying the groundwork for using marker-assisted selection in improving disease resistance and essential oil traits in mints.

Different scales of gene duplications occurring at different times have jointly shaped the NBS-LRR genes in Prunus species

In this study, genome-wide identification, phylogenetic relationships, duplication time and selective pressure of the NBS-LRR genes, an important group of plant disease-resistance genes (R genes), were performed to uncover their genetic evolutionary patterns in the six Prunus species. A total of 1946 NBS-LRR genes were identified; specifically, 589, 361, 284, 281, 318, and 113 were identified in Prunus yedoensis, P. domestica, P. avium, P. dulcis, P. persica and P. yedoensis var. nudiflora, respectively. Two NBS-LRR gene subclasses, TIR-NBS-LRR (TNL) and non-TIR-NBS-LRR (non-TNL), were also discovered. In total, 435 TNL and 1511 non-TNL genes were identified and could be classified into 30/55/75 and 103/158/191 multi-gene families, respectively, according to three different criteria. Higher Ks and Ka/Ks values were detected in TNL gene families than in non-TNL gene families. These results indicated that the TNL genes had more members involved in relatively ancient duplications and were affected by stronger selection pressure than the non-TNL genes. In general, the NBS-LRR genes were shaped by species-specific duplications, and lineage-specific duplications occurred at recent and relatively ancient periods among the six Prunus species. Therefore, different duplicated copies of NBS-LRRs can resist specific pathogens and will provide an R-gene library for resistance breeding in Prunus species.

Adaptive evolution driving the young duplications in six Rosaceae species

Background

In plant genomes, high proportions of duplicate copies reveals that gene duplications play an important role in the evolutionary processes of plant species. A series of gene families under positive selection after recent duplication events in plant genomes indicated the evolution of duplicates driven by adaptive evolution. However, the genome-wide evolutionary features of young duplicate genes among closely related species are rarely reported.

Results

In this study, we conducted a systematic survey of young duplicate genes at genome-wide levels among six Rosaceae species, whose whole-genome sequencing data were successively released in recent years. A total of 35,936 gene families were detected among the six species, in which 60.25% were generated by young duplications. The 21,650 young duplicate gene families could be divided into two expansion types based on their duplication patterns, species-specific and lineage-specific expansions. Our results showed the species-specific expansions advantaging over the lineage-specific expansions. In the two types of expansions, high-frequency duplicate domains exhibited functional preference in response to environmental stresses.

Conclusions

The functional preference of the young duplicate genes in both the expansion types showed that they were inclined to respond to abiotic or biotic stimuli. Moreover, young duplicate genes under positive selection in both species-specific and lineage-specific expansions suggested that they were generated to adapt to the environmental factors in Rosaceae species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07422-7.

Contrasting genetic variation and positive selection followed the divergence of NBS-encoding genes in Asian and European pears

BACKGROUND

The NBS disease-related gene family coordinates the inherent immune system in plants in response to pathogen infections. Previous studies have identified NBS-encoding genes in Pyrus bretschneideri ('Dangshansuli', an Asian pear) and Pyrus communis ('Bartlett', a European pear) genomes, but the patterns of genetic variation and selection pressure on these genes during pear domestication have remained unsolved.

RESULTS

In this study, 338 and 412 NBS-encoding genes were identified from Asian and European pear genomes. This difference between the two pear species was the result of proximal duplications. About 15.79% orthologous gene pairs had Ka/Ks ratio more than one, indicating two pear species undergo strong positive selection after the divergence of Asian and European pear. We identified 21 and 15 NBS-encoding genes under fire blight and black spot disease-related QTL, respectively, suggesting their importance in disease resistance. Domestication caused decreased nucleotide diversity across NBS genes in Asian cultivars (cultivated 6.23E-03; wild 6.47E-03), but opposite trend (cultivated 6.48E-03; wild 5.91E-03) appeared in European pears. Many NBS-encoding coding regions showed Ka/Ks ratio of greater than 1, indicating the role of positive selection in shaping diversity of NBS-encoding genes in pear. Furthermore, we detected 295 and 122 significantly different SNPs between wild and domesticated accessions in Asian and European pear populations. Two NBS genes (Pbr025269.1 and Pbr019876.1) with significantly different SNPs showed >5x upregulation between wild and cultivated pear accessions, and > 2x upregulation in Pyrus calleryana after inoculation with Alternaria alternata. We propose that positively selected and significantly different SNPs of an NBS-encoding gene (Pbr025269.1) regulate gene expression differences in the wild and cultivated groups, which may affect resistance in pear against A. alternata.

CONCLUSION

Proximal duplication mainly led to the different number of NBS-encoding genes in P. bretschneideri and P. communis genomes. The patterns of genetic diversity and positive selection pressure differed between Asian and European pear populations, most likely due to their independent domestication events. This analysis helps us understand the evolution, diversity, and selection pressure in the NBS-encoding gene family in Asian and European populations, and provides opportunities to study mechanisms of disease resistance in pear.

Identification and Analysis of NBS-LRR Genes in Actinidia chinensis Genome

Nucleotide-binding site and leucine-rich repeat (NBS-LRR) genes represent the most important disease resistance genes in plants. The genome sequence of kiwifruit (Actinidia chinensis) provides resources for the characterization of NBS-LRR genes and identification of new R-genes in kiwifruit. In the present study, we identified 100 NBS-LRR genes in the kiwifruit genome and they were grouped into six distinct classes based on their domain architecture. Of the 100 genes, 79 are truncated non-regular NBS-LRR genes. Except for 37 NBS-LRR genes with no location information, the remaining 63 genes are distributed unevenly across 18 kiwifruit chromosomes and 38.01% of them are present in clusters. Seventeen families of cis-acting elements were identified in the promoters of the NBS-LRR genes, including AP2, NAC, ERF and MYB. Pseudomonas syringae pv. actinidiae (pathogen of the kiwifruit bacterial canker) infection induced differential expressions of 16 detected NBS-LRR genes and three of them are involved in plant immunity responses. Our study provides insight of the NBS-LRR genes in kiwifruit and a resource for the identification of new R-genes in the fruit.

The R2R3-MYB gene family in banana (Musa acuminata): Genome-wide identification, classification and expression patterns

The R2R3-MYB genes comprise one of the largest transcription factor gene families in plants, playing regulatory roles in plant-specific developmental processes, defense responses and metabolite accumulation. To date MYB family genes have not yet been comprehensively identified in the major staple fruit crop banana. In this study, we present a comprehensive, genome-wide analysis of the MYB genes from Musa acuminata DH-Pahang (A genome). A total of 285 R2R3-MYB genes as well as genes encoding three other classes of MYB proteins containing multiple MYB repeats were identified and characterised with respect to structure and chromosomal organisation. Organ- and development-specific expression patterns were determined from RNA-Seq data. For 280 M. acuminata MYB genes for which expression was found in at least one of the analysed samples, a variety of expression patterns were detected. The M. acuminata R2R3-MYB genes were functionally categorised, leading to the identification of seven clades containing only M. acuminata R2R3-MYBs. The encoded proteins may have specialised functions that were acquired or expanded in Musa during genome evolution. This functional classification and expression analysis of the MYB gene family in banana establishes a solid foundation for future comprehensive functional analysis of MaMYBs and can be utilized in banana improvement programmes.

Resistance Gene Analogs in the Brassicaceae: Identification, Characterization, Distribution, and Evolution

The Brassicaceae consists of a wide range of species, including important Brassica crop species and the model plant Arabidopsis (Arabidopsis thaliana). Brassica spp. crop diseases impose significant yield losses annually. A major way to reduce susceptibility to disease is the selection in breeding for resistance gene analogs (RGAs). Nucleotide binding site-leucine rich repeats (NLRs), receptor-like kinases (RLKs), and receptor-like proteins (RLPs) are the main types of RGAs; they contain conserved domains and motifs and play specific roles in resistance to pathogens. Here, all classes of RGAs have been identified using annotation and assembly-based pipelines in all available genome annotations from the Brassicaceae, including multiple genome assemblies of the same species where available (total of 32 genomes). The number of RGAs, based on genome annotations, varies within and between species. In total 34,065 RGAs were identified, with the majority being RLKs (21,691), then NLRs (8,588) and RLPs (3,786). Analysis of the RGA protein sequences revealed a high level of sequence identity, whereby 99.43% of RGAs fell into several orthogroups. This study establishes a resource for the identification and characterization of RGAs in the Brassicaceae and provides a framework for further studies of RGAs for an ultimate goal of assisting breeders in improving resistance to plant disease.

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

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

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.

Segmental and Tandem Duplications Driving the Recent NBS-LRR Gene Expansion in the Asparagus Genome

Garden asparagus is an important horticultural plant worldwide. It is, however, susceptible to a variety of diseases, which can affect the potential yield, spear quality, and lifespan of production fields. Screening studies have identified resistant germplasm. The genetic resistance is usually complex, and the genes underlying that resistance are still unknown. Most often, disease resistance is determined by resistance genes (R). The most predominant R-genes contain nucleotide binding site and leucine-rich repeat (NBS-LRR) domains. Using bioinformatics and data mining approaches, we identified and characterized 68 NBS predicted proteins encoded by 49 different loci in the asparagus genome. The NBS-encoding genes were grouped into seven distinct classes based on their domain architecture. The NBS genes are unevenly distributed through the genome and nearly 50% of the genes are present in clusters. Chromosome 6 is significantly NBS-enriched and one single cluster hosts 10% of the genes. Phylogenetic analysis points to their diversification into three families during their evolution. Recent duplications are likely to have dominated the NBS expansion with both tandem genes and duplication events across multiple chromosomes. Transcriptome sequencing data provided evidence for their transcription and tissue-specific expression. The total number of cis-regulatory elements as well as their relative positions within the NBS promoters suggests a complex transcriptional network regulating defense responses. Our study provides a strong groundwork for the isolation of candidate R-genes in garden asparagus.

Sorbitol Modulates Resistance to Alternaria alternata by Regulating the Expression of an NLR Resistance Gene in Apple

In plant-microbe interactions, plant sugars produced by photosynthesis are not only a carbon source for pathogens, but may also act as signals that modulate plant defense responses. Here, we report that decreasing sorbitol synthesis in apple (Malus domestica) leaves by antisense suppression of ALDOSE-6-PHOSPHATE REDUCTASE (A6PR) leads to downregulation of 56 NUCLEOTIDE BINDING/LEUCINE-RICH REPEAT (NLR) genes and converts the phenotypic response to Alternaria alternata from resistant to susceptible. We identified a resistance protein encoded by the apple MdNLR16 gene and a small protein encoded by the fungal HRIP1 gene that interact in both a yeast two-hybrid assay and a bimolecular fluorescence complementation assay. Deletion of HRIP1 in A. alternata enables gain of virulence on the wild-type control plant. Overexpression of MdNLR16 in two antisense A6PR lines increases resistance, whereas RNAi suppression of MdNLR16 in the wild-type control decreases resistance against A. alternata MdWRKY79 transcriptionally regulates MdNLR16 by binding to the promoter of MdNLR16 in response to sorbitol, and exogenous sorbitol feeding partially restores resistance of the antisense A6PR lines to A. alternata These findings indicate that sorbitol modulates resistance to A. alternata via the MdNLR16 protein that interacts with the fungal effector in a classic gene-for-gene manner in apple.

Genome-scale examination of NBS-encoding genes in blueberry

Blueberry is an important crop worldwide. It is, however, susceptible to a variety of diseases, which can lead to losses in yield and fruit quality. Although screening studies have identified resistant germplasm for some important diseases, still little is known about the molecular basis underlying that resistance. The most predominant type of resistance (R) genes contains nucleotide binding site and leucine rich repeat (NBS-LRR) domains. The identification and characterization of such a gene family in blueberry would enhance the foundation of knowledge needed for its genetic improvement. In this study, we searched for and found a total of 106 NBS-encoding genes (including 97 NBS-LRR) in the current blueberry genome. The NBS genes were grouped into eleven distinct classes based on their domain architecture. More than 22% of the NBS genes are present in clusters. Ten genes were mapped onto seven linkage groups. Phylogenetic analysis grouped these genes into two major clusters based on their structural variation, the first cluster having toll and interleukin-1 like receptor (TIR) domains and most of the second cluster containing a coiled-coil domain. Our study provides new insight into the NBS gene family in blueberry and is an important resource for the identification of functional R-genes.

Lineage-specific duplications of NBS-LRR genes occurring before the divergence of six Fragaria species

Background

Plant disease resistance (R) genes are evolving rapidly and play a critical role in the innate immune system of plants. The nucleotide binding sites-leucine rich repeat (NBS-LRR) genes are one of the largest classes in plant R genes. Previous studies have focused on the NBS-LRR genes from one or several species of different genera, and the sequenced genomes of the genus Fragaria offer the opportunity to study the evolutionary processes of these R genes among the closely related species.

Results

In this study, 325, 155, 190, 187, and 133 NBS-LRRs were discovered from F. x ananassa, F. iinumae, F. nipponica, F. nubicola, and F. orientalis, respectively. Together with the 144 NBS-LRR genes from F. vesca, a total of 1134 NBS-LRRs containing 866 multi-genes comprised 184 gene families across the six Fragaria genomes. Extremely short branch lengths and shallow nodes were widely present in the phylogenetic tree constructed with all of the NBS-LRR genes of the six strawberry species. The identities of the orthologous genes were highly significantly greater than those of the paralogous genes, while the Ks ratios of the former were very significantly lower than those of the latter in all of the NBS-LRR gene families. In addition, the Ks and Ka/Ks values of the TIR-NBS-LRR genes (TNLs) were significantly greater than those of the non-TIR-NBS-LRR genes (non-TNLs). Furthermore, the expression patterns of the NBS-LRR genes revealed that the same gene expressed differently under different genetic backgrounds in response to pathogens.

Conclusions

These results, combined with the shared hotspot regions of the duplicated NBS-LRRs on the chromosomes, indicated that the lineage-specific duplication of the NBS-LRR genes occurred before the divergence of the six Fragaria species. The Ks and Ka/Ks ratios suggested that the TNLs are more rapidly evolving and driven by stronger diversifying selective pressures than the non-TNLs.

Electronic supplementary material

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

Genome-wide analysis suggests high level of microsynteny and purifying selection affect the evolution of EIN3/EIL family in Rosaceae

The ethylene-insensitive3/ethylene-insensitive3-like (EIN3/EIL) proteins are a type of nuclear-localized protein with DNA-binding activity in plants. Although the EIN3/EIL gene family has been studied in several plant species, little is known about comprehensive study of the EIN3/EIL gene family in Rosaceae. In this study, ten, five, four, and five EIN3/EIL genes were identified in the genomes of pear (Pyrus bretschneideri), mei (Prunus mume), peach (Prunus persica) and strawberry (Fragaria vesca), respectively. Twenty-eight chromosomal segments of EIL/EIN3 gene family were found in four Rosaceae species, and these segments could form seven orthologous or paralogous groups based on interspecies or intraspecies gene colinearity (microsynteny) analysis. Moreover, the highly conserved regions of microsynteny were found in four Rosaceae species. Subsequently it was found that both whole genome duplication and tandem duplication events significantly contributed to the EIL/EIN3 gene family expansion. Gene expression analysis of the EIL/EIN3 genes in the pear revealed subfunctionalization for several PbEIL genes derived from whole genome duplication. It is noteworthy that according to environmental selection pressure analysis, the strong purifying selection should dominate the maintenance of the EIL/EIN3 gene family in four Rosaceae species. These results provided useful information on Rosaceae EIL/EIN3 genes, as well as insights into the evolution of this gene family in four Rosaceae species. Furthermore, high level of microsynteny in the four Rosaceae plants suggested that a large-scale genome duplication event in the EIL/EIN3 gene family was predated to speciation.

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.

A unique RPW8-encoding class of genes that originated in early land plants and evolved through domain fission, fusion, and duplication

Duplication, lateral gene transfer, domain fusion/fission and de novo domain creation play a key role in formation of initial common ancestral protein. Abundant protein diversities are produced by domain rearrangements, including fusions, fissions, duplications, and terminal domain losses. In this report, we explored the origin of the RPW8 domain and examined the domain rearrangements that have driven the evolution of RPW8-encoding genes in land plants. The RPW8 domain first emerged in the early land plant, Physcomitrella patens, and it likely originated de novo from a non-coding sequence or domain divergence after duplication. It was then incorporated into the NBS-LRR protein to create a main sub-class of RPW8-encoding genes, the RPW8-NBS-encoding genes. They evolved by a series of genetic events of domain fissions, fusions, and duplications. Many species-specific duplication events and tandemly duplicated clusters clearly demonstrated that species-specific and tandem duplications played important roles in expansion of RPW8-encoding genes, especially in gymnosperms and species of the Rosaceae. RPW8 domains with greater Ka/Ks values than those of the NBS domains indicated that they evolved faster than the NBS domains in RPW8-NBSs.

Expression of resistance gene analogs in woodland strawberry (Fragaria vesca) during infection with Phytophthora cactorum

Important losses in strawberry production are often caused by the oomycete Phytophthora cactorum, the causal agent of crown rot. However, very limited studies at molecular levels exist of the mechanisms related to strawberry resistance against this pathogen. To begin to rectify this situation, a PCR-based approach (NBS profiling) was used to isolate strawberry resistance gene analogs (RGAs) with altered expression in response to P. cactorum during a time course (2, 4, 6, 24, 48, 96 and 192 h post-infection). Twenty-three distinct RGA fragments of the NB-LRR type were identified from a resistance genotype (Bukammen) of the wild species Fragaria vesca. The gene transcriptional profiles after infection showed that the response of most RGAs was quicker and stronger in the resistance genotype (Bukammen) than in the susceptible one (FDP821) during the early infection stage. The transcriptional patterns of one RGA (RGA109) were further monitored and compared during the P. cactorum infection of two pairs of resistant and susceptible genotype combinations (Bukammen/FDP821 and FDR1218/1603). The 5′ end sequence was cloned, and its putative protein was characteristic of NBS-LRR R protein. Our results yielded a first insight into the strawberry RGAs responding to P. cactorum infection at molecular level.

Resistance to root-knot nematodes Meloidogyne spp. in woody plants

I. 42 II. 43 III. 44 IV. 47 V. 49 VI. 50 VII. 50 VIII. 50 IX. 52 52 References 52 SUMMARY: Root-knot nematodes (RKNs) Meloidogyne spp. cause major damage to cultivated woody plants. Among them, Prunus, grapevine and coffee are the crops most infested by worldwide polyphagous species and species with a more limited distribution and/or narrower host range. The identification and characterization of natural sources of resistance are important steps to develop RKN control strategies. In woody crops, resistant rootstocks genetically different from the scion of agronomical interest may be engineered. We describe herein the interactions between RKNs and different woody crops, and highlight the plant species in which resistance and corresponding resistance (R) genes have been discovered. Even though grapevine and, to a lesser extent, coffee have a history of rootstock selection for RKN resistance, few cases of resistance have been documented. By contrast, in Prunus, R genes with different spectra have been mapped in plums, peach and almond and can be pyramided for durable resistance in interspecific rootstocks. We particularly discuss here the Ma Toll/interleukin-1 receptor-like-nucleotide binding-leucine-rich repeat gene from Myrobalan plum, one of the longest plant R genes cloned to date, due to its unique biological and structural properties. RKN R genes in Prunus will enable us to carry out molecular studies aimed at improving our knowledge of plant immunity in woody plants.

Species-specific duplications of NBS-encoding genes in Chinese chestnut (Castanea mollissima)

The disease resistance (R) genes play an important role in protecting plants from infection by diverse pathogens in the environment. The nucleotide-binding site (NBS)-leucine-rich repeat (LRR) class of genes is one of the largest R gene families. Chinese chestnut (Castanea mollissima) is resistant to Chestnut Blight Disease, but relatively little is known about the resistance mechanism. We identified 519 NBS-encoding genes, including 374 NBS-LRR genes and 145 NBS-only genes. The majority of Ka/Ks were less than 1, suggesting the purifying selection operated during the evolutionary history of NBS-encoding genes. A minority (4/34) of Ka/Ks in non-TIR gene families were greater than 1, showing that some genes were under positive selection pressure. Furthermore, Ks peaked at a range of 0.4 to 0.5, indicating that ancient duplications arose during the evolution. The relationship between Ka/Ks and Ks indicated greater selective pressure on the newer and older genes with the critical value of Ks = 0.4–0.5. Notably, species-specific duplications were detected in NBS-encoding genes. In addition, the group of RPW8-NBS-encoding genes clustered together as an independent clade located at a relatively basal position in the phylogenetic tree. Many cis-acting elements related to plant defense responses were detected in promoters of NBS-encoding genes.