TIR-NBS-LRR genes are rare in monocots: evidence from diverse monocot orders

High-Copy Transposons from a Pathogen Give Rise to a Conserved sRNA Family with a Novel Host Immunity Target

Small RNAs (sRNAs) are involved in gene silencing in multiple ways, including through cross-kingdom transfers from parasites to their hosts. Little is known about the evolutionary mechanisms enabling eukaryotic microbes to evolve functional mimics of host small regulatory RNAs. Here, we describe the identification and functional characterization of SINE_sRNA1, an sRNA family derived from highly abundant short interspersed nuclear element (SINE) retrotransposons in the genome of the wheat powdery mildew pathogen. SINE_sRNA1 is encoded by a sequence motif that is conserved in multiple SINE families and corresponds to a functional plant microRNA (miRNA) mimic targeting Tae_AP1, a wheat gene encoding an aspartic protease only found in monocots. Tae_AP1 has a novel function enhancing both pattern-triggered immunity (PTI) and effector-triggered immunity (ETI), thereby contributing to the cross activation of plant defenses. We conclude that SINE_sRNA1 and Tae_AP1 are functional innovations, suggesting the contribution of transposons to the evolutionary arms race between a parasite and its host. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Transcriptome analysis of genes involved in the pathogenesis mechanism of potato virus Y in potato cultivar YouJin

Introduction

Potatoes (Solanum tuberosum L.) can be infected by various viruses, but out of all of viruses, the potato virus Y (PVY) is the most detrimental. Research shows that the potato cultivar YouJin is especially vulnerable to PVY and displays severe symptoms, including leaf vein chlorosis, curled leaf margins, large necrotic spots on the leaf blades, and the growth of small new leaves.

Methods

PVY infection in potato cultivar YouJin was confirmed through symptom observation, RT-PCR, and Western blot analysis. Transcriptome sequencing was used to analyze the genes associated with PVY pathogenesis in this cultivar.

Result

Transcriptome analysis of differential genes was conducted in this study to examine the pathogenesis of PVY on YouJin. The results showed that 1,949 genes were differentially regulated, including 853 upregulated genes and 1,096 downregulated genes. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that carbohydrate synthesis and metabolism pathways were suppressed, and electron transferase and hydrolase activities were reduced. Moreover, there were increased expression levels of protein kinase genes. By focusing on plant-pathogen interaction pathways, six core genes all upregulating the WARK family of transcription factors were obtained. Additionally, a constructed PPI network revealed the identification of key modular differential genes, such as downregulated photosynthesis-related protein genes and upregulated AP2/ERF-ERF transcription factors. Functional network enrichment analysis revealed that PVY infection limited RNA metabolism, glutathionylation, and peroxiredoxin activity while triggering the expression of associated defense genes in YouJin. After analyzing the above, 26 DEGs were screened and 12 DEGs were confirmed via RT-qPCR.

Conclusion

These results establish a hypothetical framework for clarifying the pathogenesis of PVY in the YouJin variety of potatoes, which will help design the disease resistance of YouJin.

The American Cherimoya Genome Reveals Insights into the Intra-Specific Divergence, the Evolution of Magnoliales, and a Putative Gene Cluster for Acetogenin Biosynthesis

Annona cherimola (cherimoya) is a species renowned for its delectable fruit and medicinal properties. In this study, we developed a chromosome-level genome assembly for the cherimoya 'Booth' cultivar from the United States. The genome assembly has a size of 794 Mb with a N50 = 97.59 Mb. The seven longest scaffolds account for 87.6% of the total genome length, which corresponds to the seven pseudo-chromosomes. A total of 45,272 protein-coding genes (≥30 aa) were predicted with 92.9% gene content completeness. No recent whole genome duplications were identified by an intra-genome collinearity analysis. Phylogenetic analysis supports that eudicots and magnoliids are more closely related to each other than to monocots. Moreover, the Magnoliales was found to be more closely related to the Laurales than the Piperales. Genome comparison revealed that the 'Booth' cultivar has 200 Mb less repeats than the Spanish cultivar 'Fino de Jete', despite their highly similar (>99%) genome sequence identity and collinearity. These two cultivars were diverged during the early Pleistocene (1.93 Mya), which suggests a different origin and domestication of the cherimoya. Terpene/terpenoid metabolism functions were found to be enriched in Magnoliales, while TNL (Toll/Interleukin-1-NBS-LRR) disease resistance gene has been lost in Magnoliales during evolution. We have also identified a gene cluster that is potentially responsible for the biosynthesis of acetogenins, a class of natural products found exclusively in Annonaceae. The cherimoya genome provides an invaluable resource for supporting characterization, conservation, and utilization of Annona genetic resources.

Identification of QTLs for resistance to 10 pathotypes of Plasmodiophora brassicae in Brassica oleracea cultivar ECD11 through genotyping-by-sequencing

KEY MESSAGE

Two major quantitative trait loci (QTLs) and five minor QTLs for 10 pathotypes were identified on chromosomes C01, C03, C04 and C08 through genotyping-by-sequencing from Brassica oleracea. Clubroot caused by Plasmodiophora brassicae is an important disease in brassica crops. Managing clubroot disease of canola on the Canadian prairie is challenging due to the continuous emergence of new pathotypes. Brassica oleracea is considered a major source of quantitative resistance to clubroot. Genotyping-by-sequencing (GBS) was performed in the parental lines; T010000DH3 (susceptible), ECD11 (resistant) and 124 BC1 plants. A total of 4769 high-quality polymorphic SNP loci were obtained and distributed on 9 chromosomes of B. oleracea. Evaluation of 124 BC1S1 lines for resistance to 10 pathotypes: 3A, 2B, 5C, 3D, 5G, 3H, 8J, 5K, 5L and 3O of P. brassicae, was carried out. Seven QTLs, 5 originating from ECD11 and 2 from T010000DH3, were detected. One major QTL designated as Rcr_C03-1 on C03 contributed 16.0-65.6% of phenotypic variation explained (PVE) for 8 pathotypes: 2B, 5C, 5G, 3H, 8J, 5K, 5L and 3O. Another major QTL designated as Rcr_C08-1 on C08 contributed 8.3 and 23.5% PVE for resistance to 8J and 5K, respectively. Five minor QTLs designated as Rcr_C01-1, Rcr_C03-2, Rcr_C03-3, Rcr_C04-1 and Rcr_C08-2 were detected on chromosomes C01, C03, C04 and C08 that contributed 8.3-23.5% PVE for 5 pathotypes each of 3A, 2B, 3D, 8J and 5K. There were 1, 10 and 4 genes encoding TIR-NBS-LRR/CC-NBS-LRR class disease resistance proteins in the Rcr_C01-1, Rcr_C03-1 and Rcr_C08-1 flanking regions. The syntenic regions of the two major QTLs Rcr_C03-1 and Rcr_C08-1 in the B. rapa genome 'Chiifu' were searched.

Comparative genomics profiling revealed multi-stress responsive roles of the CC-NBS-LRR genes in three mango cultivars

The nucleotide-binding site-leucine-rich repeat (NBS-LRR) gene family is the largest group of disease resistance (R) genes in plants and is active in response to viruses, bacteria, and fungi usually involved in effector-triggered immunity (ETI). Pangenome-wide studies allow researchers to analyze the genetic diversity of multiple species or their members simultaneously, providing a comprehensive understanding of the evolutionary relationships and diversity present among them. The draft pan-genome of three Mangifera indica cultivars (Alphonso, Hong Xiang Ya, and Tommy atkins) was constructed and Presence/absence variants (PAVs) were filtered through the ppsPCP pipeline. As a result, 2823 genes and 5907 PAVs from H. Xiang Ya, and 1266 genes and 2098 PAVs from T. atkins were added to the reference genome. For the identification of CC-NBS-LRR (CNL) genes in these mango cultivars, this draft pan-genome study has successfully identified 47, 27, and 36 members in Alphonso, H. Xiang Ya, and T. atkins respectively. The phylogenetic analysis divided MiCNL proteins into four distinct subgroups. All MiCNL genes are unevenly distributed on chromosomes. Both tandem and segmental duplication events played a significant role in the expansion of the CNL gene family. These genes contain cis-elements related to light, stress, hormone, and development. The analysis of protein-protein interactions (PPI) revealed that MiCNL proteins interacted with other defense-responsive proteins. Gene Ontology (GO) analysis indicated that MiCNL genes play a role in defense mechanisms within the organism. The expression level of the identified genes in fruit peel was observed under disease and cold stress which showed that Mi_A_CNL13 and 14 were up-regulated while Mi_A_CNL15, 25, 30, 31, and 40 were down-regulated in disease stress. On the other hand, Mi_A_CNL2, 14, 41, and 45 were up-regulated and Mi_A_CNL47 is down-regulated in cold stress. Subsequently, the Random Forest (RF) classifier was used to assess the multi-stress response of MiCNLs. It was found that Mi_A_CNL14 is a gene that responds to multiple stress conditions. The CNLs have similar protein structures which show that they are involved in the same function. The above findings provide a foundation for a deeper understanding of the functional characteristics of the mango CNL gene family.

The genome of Acorus deciphers insights into early monocot evolution

Acorales is the sister lineage to all the other extant monocot plants. Genomic resource enhancement of this genus can help to reveal early monocot genomic architecture and evolution. Here, we assemble the genome of Acorus gramineus and reveal that it has ~45% fewer genes than the majority of monocots, although they have similar genome size. Phylogenetic analyses based on both chloroplast and nuclear genes consistently support that A. gramineus is the sister to the remaining monocots. In addition, we assemble a 2.2 Mb mitochondrial genome and observe many genes exhibit higher mutation rates than that of most angiosperms, which could be the reason leading to the controversies of nuclear genes- and mitochondrial genes-based phylogenetic trees existing in the literature. Further, Acorales did not experience tau (τ) whole-genome duplication, unlike majority of monocot clades, and no large-scale gene expansion is observed. Moreover, we identify gene contractions and expansions likely linking to plant architecture, stress resistance, light harvesting, and essential oil metabolism. These findings shed light on the evolution of early monocots and genomic footprints of wetland plant adaptations.

The Role of Plant Transcription Factors in the Fight against Plant Viruses

Plants are vulnerable to the challenges of unstable environments and pathogen infections due to their immobility. Among various stress conditions, viral infection is a major threat that causes significant crop loss. In response to viral infection, plants undergo complex molecular and physiological changes, which trigger defense and morphogenic pathways. Transcription factors (TFs), and their interactions with cofactors and cis-regulatory genomic elements, are essential for plant defense mechanisms. The transcriptional regulation by TFs is crucial in establishing plant defense and associated activities during viral infections. Therefore, identifying and characterizing the critical genes involved in the responses of plants against virus stress is essential for the development of transgenic plants that exhibit enhanced tolerance or resistance. This article reviews the current understanding of the transcriptional control of plant defenses, with a special focus on NAC, MYB, WRKY, bZIP, and AP2/ERF TFs. The review provides an update on the latest advances in understanding how plant TFs regulate defense genes expression during viral infection.

Genome-wide analysis of NBS-LRR genes revealed contribution of disease resistance from Saccharum spontaneum to modern sugarcane cultivar

INTRODUCTION

During plant evolution, nucleotide-binding sites (NBS) and leucine-rich repeat (LRR) genes have made significant contributions to plant disease resistance. With many high-quality plant genomes sequenced, identification and comprehensive analyses of NBS-LRR genes at whole genome level are of great importance to understand and utilize them.

METHODS

In this study, we identified the NBS-LRR genes of 23 representative species at whole genome level, and researches on NBS-LRR genes of four monocotyledonous grass species, Saccharum spontaneum, Saccharum officinarum, Sorghum bicolor and Miscanthus sinensis, were focused.

RESULTS AND DISCUSSION

We found that whole genome duplication, gene expansion, and allele loss could be factors affecting the number of NBS-LRR genes in the species, and whole genome duplication is likely to be the main cause of the number of NBS-LRR genes in sugarcane. Meanwhile, we also found a progressive trend of positive selection on NBS-LRR genes. These studies further elucidated the evolutionary pattern of NBS-LRR genes in plants. Transcriptome data from multiple sugarcane diseases revealed that more differentially expressed NBS-LRR genes were derived from S. spontaneum than from S. officinarum in modern sugarcane cultivars, and the proportion was significantly higher than the expected. This finding reveals that S. spontaneum has a greater contribution to disease resistance for modern sugarcane cultivars. In addition, we observed allelespecific expression of seven NBS-LRR genes under leaf scald, and 125 NBS-LRR genes responding to multiple diseases were identified. Finally, we built a plant NBS-LRR gene database to facilitate subsequent analysis and use of NBSLRR genes obtained here. In conclusion, this study complemented and completed the research of plant NBS-LRR genes, and discussed how NBS-LRR genes responding to sugarcane diseases, which provided a guide and genetic resources for further research and utilization of NBS-LRR genes.

Legumes Regulate Symbiosis with Rhizobia via Their Innate Immune System

Plant roots are constantly exposed to a diverse microbiota of pathogens and mutualistic partners. The host's immune system is an essential component for its survival, enabling it to monitor nearby microbes for potential threats and respond with a defence response when required. Current research suggests that the plant immune system has also been employed in the legume-rhizobia symbiosis as a means of monitoring different rhizobia strains and that successful rhizobia have evolved to overcome this system to infect the roots and initiate nodulation. With clear implications for host-specificity, the immune system has the potential to be an important target for engineering versatile crops for effective nodulation in the field. However, current knowledge of the interacting components governing this pathway is limited, and further research is required to build on what is currently known to improve our understanding. This review provides a general overview of the plant immune system's role in nodulation. With a focus on the cycles of microbe-associated molecular pattern-triggered immunity (MTI) and effector-triggered immunity (ETI), we highlight key molecular players and recent findings while addressing the current knowledge gaps in this area.

Advances in Biological Control and Resistance Genes of Brassicaceae Clubroot Disease-The Study Case of China

Clubroot disease is a soil-borne disease caused by Plasmodiophora brassicae. It occurs in cruciferous crops exclusively, and causes serious damage to the economic value of cruciferous crops worldwide. Although different measures have been taken to prevent the spread of clubroot disease, the most fundamental and effective way is to explore and use disease-resistance genes to breed resistant varieties. However, the resistance level of plant hosts is influenced both by environment and pathogen race. In this work, we described clubroot disease in terms of discovery and current distribution, life cycle, and race identification systems; in particular, we summarized recent progress on clubroot control methods and breeding practices for resistant cultivars. With the knowledge of these identified resistance loci and R genes, we discussed feasible strategies for disease-resistance breeding in the future.

Bioinformatics analysis and function prediction of NBS-LRR gene family in Broussonetia papyrifera

Most of the currently available disease resistance (R) genes have NBS (nucleotide-binding site) and LRR (leucine-rich-repeat) domain which belongs to the NBS-LRR gene family. The whole genome sequencing of Broussonetia papyrifera provides an important bioinformatics database for the study of the NBS-LRR gene family. In this study, 328 NBS-LRR family genes were identified and classified in B. papyrifera according to different classification schemes, where there are 92 N types, 47 CN type, 54 CNL type, 29 NL types, 55 TN type, and 51 TNL type. Subsequently, we conducted bioinformatics analysis of the NBS-LRR gene family. Classification, motif analysis of protein sequences, and phylogenetic tree studies of the NBS-LRR genes in B. papyrifera provide important basis for the functional study of NBS-LRR family genes. Additionally, we performed structural analysis of the chromosomal location, physicochemical properties, and sequences identified by genetic characterization. In addition, through the analysis of GO enrichment, it was found that NBS-LRR genes were involved in defense responses and were significantly enriched in biological stimulation, immune response, and abiotic stress. In addition, we found that Bp06g0955 was the most sensitive to low temperature and encoded the RPM1 protein by analyzing the low temperature transcriptome data of B. papyrifera. Quantitative results of gene expression after 48 h of Fusarium infection showed that Bp01g3293 increased 14 times after infection, which encodes RPM1 protein. The potential of NBS-LRR gene responsive to biotic and abiotic stresses can be exploited to improve the resistance of B. papyrifera.

Thirty years of resistance: Zig-zag through the plant immune system

Understanding the plant immune system is crucial for using genetics to protect crops from diseases. Plants resist pathogens via a two-tiered innate immune detection-and-response system. The first plant Resistance (R) gene was cloned in 1992 . Since then, many cell-surface pattern recognition receptors (PRRs) have been identified, and R genes that encode intracellular nucleotide-binding leucine-rich repeat receptors (NLRs) have been cloned. Here, we provide a list of characterized PRRs and NLRs. In addition to immune receptors, many components of immune signaling networks were discovered over the last 30 years. We review the signaling pathways, physiological responses, and molecular regulation of both PRR- and NLR-mediated immunity. Recent studies have reinforced the importance of interactions between the two immune systems. We provide an overview of interactions between PRR- and NLR-mediated immunity, highlighting challenges and perspectives for future research.

Identification and characterization of genome-wide resistance gene analogs (RGAs) of durian (Durio zibethinus L.)

BACKGROUND

Durian (Durio zibethinus L.) is a tropical fruit crop which is popular in Southeast Asia but recently gaining popularity in other parts of the world. In this study, we analyzed the resistance gene analogs (RGAs) of durian through mining of the currently available reference genome of its 'Musang King' cultivar (PRJNA400310).

RESULTS

A total of 2586 RGAs were identified in the durian genome consisting of 47 nucleotide binding site proteins (NBS), 158 NBS-leucine rich repeat proteins (NL), 400 coiled-coil NBS-LRR (CNL), 72 toll/interleukin-1 receptor NBS-LRR (TNL), 54 coiled-coil NBS (CN), 10 toll/interleukin-1 receptor NBS (TN), 19 toll/interleukin-1 receptor with unknown domain (TX), 246 receptor-like proteins (RLP), 1,377 receptor-like kinases (RLK), 185 TM-CC, and 18 other NBS-containing proteins with other domains. These RGAs were functionally annotated and characterized via gene ontology (GO) analysis. Among the RGAs with the highest copies in durian genome include the putative disease resistance RPP13-like protein 1, disease resistance protein At4g27190, disease resistance protein RPS6, Probable disease resistance protein At4g27220, and putative disease resistance protein RGA3, while 35 RGAs were found to be novel. Phylogenetic analyses revealed that the genome-wide RGAs were broadly clustered into four major clades based on their domain classification.

CONCLUSION

To our knowledge, this is the most comprehensive analysis of durian RGAs which provides a valuable resource for genetic, agronomic, and other biological research of this important tropical fruit crop.

NB-LRR genes: characteristics in three Solanum species and transcriptional response to Ralstonia solanacearum in tomato

NB-LRR genes in the three Solanum species showed specific constitution characteristics and evolved multiple clusters and duplicates. Some genes could respond to biotic stresses such as tomato bacterial wilt. Nucleotide-binding and leucine-rich repeat (NB-LRR, NLR) is a largest resistance gene family in plants, which plays a key role in response to biotic stresses. In this study, NB-LRR genes in cultivated tomato Solanum lycopersicum (Sl) and its wild relatives S. pennellii (Spe) and S. pimpinellifolium (Spi) were analyzed using bioinformatics approaches. In total, 238, 202 and 217 NB-LRR genes of 8 different types were found in Sl, Spe and Spi, respectively. The three species showed similar genomic characteristics. The NB-LRR genes were mainly distributed on chromosomes 4, 5 and 11 and located at the distal zones, forming multiple clusters and tandem duplicates. A large number of homologs appeared through gene expansion, with most Ka/Ks values being less than 1, indicating that purifying selection had occurred in evolution. These genes were mainly expressed in root and could respond to different biotic stresses. RT-qPCR analysis revealed that SlNLR genes could respond to tomato bacterial wilt, with SlNLR1 probably involved in the resistance response, whereas others being the opposite. The transcription factors (TFs) and interaction proteins that regulate target genes were mainly Dof, NAC and MYB families and kinases. The results provide a basis for the isolation and application of related genes in plant disease resistance breeding.

Plant NLR diversity: the known unknowns of pan-NLRomes

Plants and pathogens constantly adapt to each other. As a consequence, many members of the plant immune system, and especially the intracellular nucleotide-binding site leucine-rich repeat receptors, also known as NOD-like receptors (NLRs), are highly diversified, both among family members in the same genome, and between individuals in the same species. While this diversity has long been appreciated, its true extent has remained unknown. With pan-genome and pan-NLRome studies becoming more and more comprehensive, our knowledge of NLR sequence diversity is growing rapidly, and pan-NLRomes provide powerful platforms for assigning function to NLRs. These efforts are an important step toward the goal of comprehensively predicting from sequence alone whether an NLR provides disease resistance, and if so, to which pathogens.

Genome-Wide Identification and Analysis of CC-NBS-LRR Family in Response to Downy Mildew and Black Rot in Chinese Cabbage

The nucleotide-binding site-leucine-rich repeat (NBS-LRR) gene family is the largest group of plant disease resistance (R) genes widespread in response to viruses, bacteria, and fungi usually involved in effector triggered immunity (ETI). Forty members of the Chinese cabbage CC type NBS-LRR family were investigated in this study. Gene and protein characteristics, such as distributed locations on chromosomes and gene structures, were explored through comprehensive analysis. CC-NBS-LRR proteins were classified according to their conserved domains, and the phylogenetic relationships of CC-NBS-LRR proteins in Brassica rapa, Arabidopsis thaliana, and Oryza sativa were compared. Moreover, the roles of BrCC-NBS-LRR genes involved in pathogenesis-related defense were studied and analyzed. First, the expression profiles of BrCC-NBS-LRR genes were detected by inoculating with downy mildew and black rot pathogens. Second, sensitive and resistant Chinese cabbage inbred lines were screened by downy mildew and black rot. Finally, the differential expression levels of BrCC-NBS-LRR genes were monitored at 0, 1, 3, 6, 12 and 24 h for short and 0, 3, 5, 7, 10 and 14 days for long inoculation time. Our study provides information on BrCC-NBS-LRR genes for the investigation of the functions and mechanisms of CC-NBS-LRR genes in Chinese cabbage.

Integrity of the Post-LRR Domain Is Required for TIR-NB-LRR Function

Plants trigger appropriate defense responses, notably, through intracellular nucleotide-binding (NB) and leucine-rich repeat (LRR)-containing receptors (NLRs) that detect secreted pathogen effector proteins. In NLR resistance genes, the toll/interleukin-1 receptor (TIR)-NB-LRR proteins (TNLs) are an important subfamily, out of which approximately half the members carry a post-LRR (PL) domain of unknown role. We first investigated the requirement of the PL domain for TNL-mediated immune response by mutating the most conserved amino acids across PL domains of Arabidopsis thaliana TNLs. We identified several amino acids in the PL domain of RPS4, required for its ability to trigger a hypersensitive response to AvrRps4 in a Nicotiana tabacum transient assay. Mutating the corresponding amino acids within the PL domain of the tobacco TNL gene N also affected its function. Consequently, our results indicate that the integrity of the PL domain at conserved positions is crucial for at least two unrelated TNLs. We then tested the PL domain specificity for function by swapping PL domains between the paralogs RPS4 and RPS4B. Our results suggest that the PL domain is involved in their TNL pair specificity, 'off state' stability, and NLR complex activation. Considering genetically paired Arabidopsis TNLs, we finally compared the PL and TIR domains of their sensor and executor sequences, respectively. While TIR and PL domains from executors present complete motifs, sensors showed a lack of conservation with degenerated motifs. We here provide a contribution to the functional analysis of the PL domain in order to decipher its role for TNL function.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Evolution of NLR Resistance Genes in Magnoliids: Dramatic Expansions of CNLs and Multiple Losses of TNLs

Magnoliids are the third-largest group of angiosperms and occupy a critical position in angiosperm evolution. In the past years, due to the lack of sequenced genomes, the disease resistance gene (R gene) profile of magnoliids remains poorly understood. By the genome-wide identification of 1,832 NLR genes from seven magnoliid genomes, we built a framework for the evolution of magnoliid R genes. TNL genes were completely absent from five magnoliids, presumably due to immune pathway deficiencies. A total of 74 ancestral R genes (70 CNLs, 3 TNLs, and 1 RNL) were recovered in a common ancestor of magnoliids, from which all current NLR gene repertoires were derived. Tandem duplication served as the major drive for NLR genes expansion in seven magnoliid genomes, as most surveyed angiosperms. Due to recent rapid expansions, most magnoliids exhibited "a first expansion followed by a slight contraction and a further stronger expansion" evolutionary pattern, while both Litsea cubeba and Persea americana showed a two-times-repeated pattern of "expansion followed by contraction." The transcriptome analysis of seven different tissues of Saururus chinensis revealed a low expression of most NLR genes, with some R genes displaying a relatively higher expression in roots and fruits. Overall, our study sheds light on the evolution of NLR genes in magnoliids, compensates for insufficiency in major angiosperm lineages, and provides an important reference for a better understanding of angiosperm NLR genes.

Assembly and Annotation of the Nuclear and Organellar Genomes of a Dwarf Coconut (Chowghat Green Dwarf) Possessing Enhanced Disease Resistance

Coconut (Cocos nucifera L.), an important source of vegetable oil, nutraceuticals, functional foods, and housing materials, provides raw materials for a repertoire of industries engaged in the manufacture of cosmetics, soaps, detergents, paints, varnishes, and emulsifiers, among other products. The palm plays a vital role in maintaining and promoting the sustainability of farming systems of the fragile ecosystems of islands and coastal regions of the tropics. In this study, we present the genome of a dwarf coconut variety "Chowghat Green Dwarf" (CGD) from India, possessing enhanced resistance to root (wilt) disease. Utilizing short reads from the Illumina HiSeq 4000 platform and long reads from the Pacific Biosciences RSII platform, we have assembled the draft genome assembly of 1.93 Gb. The genome is distributed over 26,855 scaffolds, with ∼81.56% of the assembled genome present in scaffolds of lengths longer than 50 kb. About 77.29% of the genome was composed of transposable elements and repeats. Gene prediction yielded 51,953 genes, which upon stringent filtering, based on Annotation Edit Distance, resulted in 13,707 genes, which coded for 11,181 proteins. Among these, we gathered transcript level evidence for a total of 6828 predicted genes based on the RNA-Seq data from different coconut tissues, since they presented assembled transcripts within the genome annotation coordinates. A total of 112 nucleotide-binding and leucine-rich repeat loci, belonging to six classes, were detected. We have also undertaken the assembly and annotation of the CGD chloroplast and mitochondrial genomes. The availability of the dwarf coconut genome shall prove invaluable for deducing the origin of dwarf coconut cultivars, dissection of genes controlling plant habit and fruit color, and accelerated breeding for improved agronomic traits.

A CNL protein in wild emmer wheat confers powdery mildew resistance

Powdery mildew, a fungal disease caused by Blumeria graminis f. sp. tritici (Bgt), has a serious impact on wheat production. Loss of resistance in cultivars prompts a continuing search for new sources of resistance. Wild emmer wheat (Triticum turgidum ssp. dicoccoides, WEW), the progenitor of both modern tetraploid and hexaploid wheats, harbors many powdery mildew resistance genes. We report here the positional cloning and functional characterization of Pm41, a powdery mildew resistance gene derived from WEW, which encodes a coiled-coil, nucleotide-binding site and leucine-rich repeat protein (CNL). Mutagenesis and stable genetic transformation confirmed the function of Pm41 against Bgt infection in wheat. We demonstrated that Pm41 was present at a very low frequency (1.81%) only in southern WEW populations. It was absent in other WEW populations, domesticated emmer, durum, and common wheat, suggesting that the ancestral Pm41 was restricted to its place of origin and was not incorporated into domesticated wheat. Our findings emphasize the importance of conservation and exploitation of the primary WEW gene pool, as a valuable resource for discovery of resistance genes for improvement of modern wheat cultivars.

Taro Genome Assembly and Linkage Map Reveal QTLs for Resistance to Taro Leaf Blight

Taro (Colocasia esculenta) is a food staple widely cultivated in the humid tropics of Asia, Africa, Pacific and the Caribbean. One of the greatest threats to taro production is Taro Leaf Blight caused by the oomycete pathogen Phytophthora colocasiae Here we describe a de novo taro genome assembly and use it to analyze sequence data from a Taro Leaf Blight resistant mapping population. The genome was assembled from linked-read sequences (10x Genomics; ∼60x coverage) and gap-filled and scaffolded with contigs assembled from Oxford Nanopore Technology long-reads and linkage map results. The haploid assembly was 2.45 Gb total, with a maximum contig length of 38 Mb and scaffold N50 of 317,420 bp. A comparison of family-level (Araceae) genome features reveals the repeat content of taro to be 82%, >3.5x greater than in great duckweed (Spirodela polyrhiza), 23%. Both genomes recovered a similar percent of Benchmarking Universal Single-copy Orthologs, 80% and 84%, based on a 3,236 gene database for monocot plants. A greater number of nucleotide-binding leucine-rich repeat disease resistance genes were present in genomes of taro than the duckweed, ∼391 vs. ∼70 (∼182 and ∼46 complete). The mapping population data revealed 16 major linkage groups with 520 markers, and 10 quantitative trait loci (QTL) significantly associated with Taro Leaf Blight disease resistance. The genome sequence of taro enhances our understanding of resistance to TLB, and provides markers that may accelerate breeding programs. This genome project may provide a template for developing genomic resources in other understudied plant species.

Convergent Loss of an EDS1/PAD4 Signaling Pathway in Several Plant Lineages Reveals Coevolved Components of Plant Immunity and Drought Response

Plant innate immunity relies on nucleotide binding leucine-rich repeat receptors (NLRs) that recognize pathogen-derived molecules and activate downstream signaling pathways. We analyzed the variation in NLR gene copy number and identified plants with a low number of NLR genes relative to sister species. We specifically focused on four plants from two distinct lineages, one monocot lineage (Alismatales) and one eudicot lineage (Lentibulariaceae). In these lineages, the loss of NLR genes coincides with loss of the well-known downstream immune signaling complex ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1)/PHYTOALEXIN DEFICIENT 4 (PAD4). We expanded our analysis across whole proteomes and found that other characterized immune genes were absent only in Lentibulariaceae and Alismatales. Additionally, we identified genes of unknown function that were convergently lost together with EDS1/PAD4 in five plant species. Gene expression analyses in Arabidopsis (Arabidopsis thaliana) and Oryza sativa revealed that several homologs of the candidates are differentially expressed during pathogen infection, drought, and abscisic acid treatment. Our analysis provides evolutionary evidence for the rewiring of plant immunity in some plant lineages, as well as the coevolution of the EDS1/PAD4 pathway and drought responses.

QTL mapping of web blotch resistance in peanut by high-throughput genome-wide sequencing

BACKGROUND

Web blotch is one of the most important foliar diseases worldwide in peanut (Arachis hypogaea L.). The identification of quantitative trait loci (QTLs) for peanut web blotch resistance represents the basis for gene mining and the application of molecular breeding technologies.

RESULTS

In this study, a peanut recombinant inbred line (RIL) population was used to map QTLs for web blotch resistance based on high-throughput genome-wide sequencing. Frequency distributions of disease grade and disease index in five environments indicated wide phenotypic variations in response to web blotch among RILs. A high-density genetic map was constructed, containing 3634 bin markers distributed on 20 peanut linkage groups (LGs) with an average genetic distance of 0.5 cM. In total, eight QTLs were detected for peanut web blotch resistance in at least two environments, explaining from 2.8 to 15.1% of phenotypic variance. Two major QTLs qWBRA04 and qWBRA14 were detected in all five environments and were linked to 40 candidate genes encoding nucleotide-binding site leucine-rich repeat (NBS-LRR) or other proteins related to disease resistances.

CONCLUSIONS

The results of this study provide a basis for breeding peanut cultivars with web blotch resistance.

Genome- Wide Analysis of the Nucleotide Binding Site Leucine-Rich Repeat Genes of Four Orchids Revealed Extremely Low Numbers of Disease Resistance Genes

Orchids are one of the most diverse flowering plant families, yet possibly maintain the smallest number of the nucleotide-binding site-leucine-rich repeat (NBS-LRR) type plant resistance (R) genes among the angiosperms. In this study, a genome-wide search in four orchid taxa identified 186 NBS-LRR genes. Furthermore, 214 NBS-LRR genes were identified from seven orchid transcriptomes. A phylogenetic analysis recovered 30 ancestral lineages (29 CNL and one RNL), far fewer than other angiosperm families. From the genetics aspect, the relatively low number of ancestral R genes is unlikely to explain the low number of R genes in orchids alone, as historical gene loss and scarce gene duplication has continuously occurred, which also contributes to the low number of R genes. Due to recent sharp expansions, Phalaenopsis equestris and Dendrobium catenatum having 52 and 115 genes, respectively, and exhibited an “early shrinking to recent expanding” evolutionary pattern, while Gastrodia elata and Apostasia shenzhenica both exhibit a “consistently shrinking” evolutionary pattern and have retained only five and 14 NBS-LRR genes, respectively. RNL genes remain in extremely low numbers with only one or two copies per genome. Notably, all of the orchid RNL genes belong to the ADR1 lineage. A separate lineage, NRG1, was entirely absent and was likely lost in the common ancestor of all monocots. All of the TNL genes were absent as well, coincident with the RNL NRG1 lineage, which supports the previously proposed notion that a potential functional association between the TNL and RNL NRG1 genes.

Genome survey of resistance gene analogs in sugarcane: genomic features and differential expression of the innate immune system from a smut-resistant genotype

BACKGROUND

Resistance genes composing the two-layer immune system of plants are thought as important markers for breeding pathogen-resistant crops. Many have been the attempts to establish relationships between the genomic content of Resistance Gene Analogs (RGAs) of modern sugarcane cultivars to its degrees of resistance to diseases such as smut. However, due to the highly polyploid and heterozygous nature of sugarcane genome, large scale RGA predictions is challenging.

RESULTS

We predicted, searched for orthologs, and investigated the genomic features of RGAs within a recently released sugarcane elite cultivar genome, alongside the genomes of sorghum, one sugarcane ancestor (Saccharum spontaneum), and a collection of de novo transcripts generated for six modern cultivars. In addition, transcriptomes from two sugarcane genotypes were obtained to investigate the roles of RGAs differentially expressed (RGADE) in their distinct degrees of resistance to smut. Sugarcane references lack RGAs from the TNL class (Toll-Interleukin receptor (TIR) domain associated to nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains) and harbor elevated content of membrane-associated RGAs. Up to 39% of RGAs were organized in clusters, and 40% of those clusters shared synteny. Basically, 79% of predicted NBS-encoding genes are located in a few chromosomes. S. spontaneum chromosome 5 harbors most RGADE orthologs responsive to smut in modern sugarcane. Resistant sugarcane had an increased number of RGAs differentially expressed from both classes of RLK (receptor-like kinase) and RLP (receptor-like protein) as compared to the smut-susceptible. Tandem duplications have largely contributed to the expansion of both RGA clusters and the predicted clades of RGADEs.

CONCLUSIONS

Most of smut-responsive RGAs in modern sugarcane were potentially originated in chromosome 5 of the ancestral S. spontaneum genotype. Smut resistant and susceptible genotypes of sugarcane have a distinct pattern of RGADE. TM-LRR (transmembrane domains followed by LRR) family was the most responsive to the early moment of pathogen infection in the resistant genotype, suggesting the relevance of an innate immune system. This work can help to outline strategies for further understanding of allele and paralog expression of RGAs in sugarcane, and the results should help to develop a more applied procedure for the selection of resistant plants in sugarcane.

The utility of NBS-profiling for characterization of yellow rust resistance in an F6 durum wheat population

Seedling and adult plant (field) resistance to yellow rust in the durum wheat (Triticum turgidum ssp. durum) cross Kunduru-1149 x Cham-1 was characterized using a functionally-targeted DNA marker system, NBS-profiling. Chi-squared analysis indicated a four gene model conferring seedling yellow rust resistance against Puccinia striiformis f. sp. tritici isolate WYR85/22 (virulent on Yr2, Yr6, Yr7 and Yr9). Interval mapping located two QTL for yellow rust resistance on the long arm of chromosome 1B, while Kruskal-Wallis single marker regression identified a number of additional marker loci associated with seedling and/or adult plant, field resistance to yellow rust. These results suggested that much of the yellow rust resistance seen in the field may be due to seedling expressed resistance (R) genes. Characterization of the DNA sequence of three NBS marker loci indicated that all showed significant homology to functionally-characterized R-genes and resistance gene analogues (RGAs), with the greatest homology being NBS-LRR-type R-genes and RGAs from cereal species.

De Novo Genome Sequence Assembly of Dwarf Coconut (Cocos nucifera L. 'Catigan Green Dwarf') Provides Insights into Genomic Variation Between Coconut Types and Related Palm Species

We report the first whole genome sequence (WGS) assembly and annotation of a dwarf coconut variety, 'Catigan Green Dwarf' (CATD). The genome sequence was generated using the PacBio SMRT sequencing platform at 15X coverage of the expected genome size of 2.15 Gbp, which was corrected with assembled 50X Illumina paired-end MiSeq reads of the same genome. The draft genome was improved through Chicago sequencing to generate a scaffold assembly that results in a total genome size of 2.1 Gbp consisting of 7,998 scaffolds with N50 of 570,487 bp. The final assembly covers around 97.6% of the estimated genome size of coconut 'CATD' based on homozygous k-mer peak analysis. A total of 34,958 high-confidence gene models were predicted and functionally associated to various economically important traits, such as pest/disease resistance, drought tolerance, coconut oil biosynthesis, and putative transcription factors. The assembled genome was used to infer the evolutionary relationship within the palm family based on genomic variations and synteny of coding gene sequences. Data show that at least three (3) rounds of whole genome duplication occurred and are commonly shared by these members of the Arecaceae family. A total of 7,139 unique SSR markers were designed to be used as a resource in marker-based breeding. In addition, we discovered 58,503 variants in coconut by aligning the Hainan Tall (HAT) WGS reads to the non-repetitive regions of the assembled CATD genome. The gene markers and genome-wide SSR markers established here will facilitate the development of varieties with resilience to climate change, resistance to pests and diseases, and improved oil yield and quality.

Identification and expression profiling analysis of NBS-LRR genes involved in Fusarium oxysporum f.sp. conglutinans resistance in cabbage

As one of the most important resistance (R) gene families in plants, the NBS-LRR genes, encoding proteins with nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains, play significant roles in resisting pathogens. The published genomic data for cabbage (Brassica oleracea L.) provide valuable data to identify and characterize the genomic organization of cabbage NBS-LRR genes. Ultimately, we identified 105 TIR (N-terminal Toll/interleukin-1 receptor)-NBS-LRR (TNL) genes and 33 CC (coiled-coil)-NBS-LRR (CNL) genes. Further research indicated that 50.7% of the 138 NBS-LRR genes exist in 27 clusters and there are large differences among the gene structures and protein characteristics. Conserved motif and phylogenetic analysis showed that the structures of TNLs and CNLs were similar, with some differences. These NBS-LRRs are evolved under negative selection and mostly arose from whole-genome duplication events during evolution. Tissue-expression profiling of NBS-LRR genes revealed that 37.1% of the TNL genes are highly or specifically expressed in roots, especially the genes on chromosome 7 (76.5%). Digital gene expression and reverse transcription PCR analyses revealed the expression patterns of the NBS-LRR genes upon challenge by Fusarium oxysporum f.sp. conglutinans: nine genes were upregulated, and five were downregulated. The major resistance gene Foc1 probably works together with the other four genes in the same cluster to resist F. oxysporum infection.

Out of Water: The Origin and Early Diversification of Plant R-Genes

During plant-pathogen interactions, plants use intracellular proteins with nucleotide-binding site and Leu-rich repeat (NBS-LRR) domains to detect pathogens. NBS-LRR proteins represent a major class of plant disease resistance genes (R-genes). Whereas R-genes have been well characterized in angiosperms, little is known about their origin and early diversification. Here, we perform comprehensive evolutionary analyses of R-genes in plants and report the identification of R-genes in basal-branching streptophytes, including charophytes, liverworts, and mosses. Phylogenetic analyses suggest that plant R-genes originated in charophytes and R-proteins diversified into TIR-NBS-LRR proteins and non-TIR-NBS-LRR proteins in charophytes. Moreover, we show that plant R-proteins evolved in a modular fashion through frequent gain or loss of protein domains. Most of the R-genes in basal-branching streptophytes underwent adaptive evolution, indicating an ancient involvement of R-genes in plant-pathogen interactions. Our findings provide novel insights into the origin and evolution of R-genes and the mechanisms underlying colonization of terrestrial environments by plants.

The Brassicaceae Family Displays Divergent, Shoot-Skewed NLR Resistance Gene Expression

Nucleotide-binding site leucine-rich repeat resistance genes (NLRs) allow plants to detect microbial effectors. We hypothesized that NLR expression patterns could reflect organ-specific differences in effector challenge and tested this by carrying out a meta-analysis of expression data for 1,235 NLRs from nine plant species. We found stable NLR root/shoot expression ratios within species, suggesting organ-specific hardwiring of NLR expression patterns in anticipation of distinct challenges. Most monocot and dicot plant species preferentially expressed NLRs in roots. In contrast, Brassicaceae species, including oilseed rape (Brassica napus) and the model plant Arabidopsis (Arabidopsis thaliana), were unique in showing NLR expression skewed toward the shoot across multiple phylogenetically distinct groups of NLRs. The Brassicaceae are also outliers in the sense that they have lost the common symbiosis signaling pathway, which enables intracellular infection by root symbionts. While it is unclear if these two events are related, the NLR expression shift identified here suggests that the Brassicaceae may have evolved unique pattern-recognition receptors and antimicrobial root metabolites to substitute for NLR protection. Such innovations in root protection could potentially be exploited in crop rotation schemes or for enhancing root defense systems of non-Brassicaceae crops.

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.

Resistance Genes in Piper colubrinum: In Silico Survey From Leaf Transcriptome and Expression Studies Upon Challenge Inoculation with Phytophthora capsici

The oomycetes, Phytophthora capsici, cause foot rot disease in black pepper. Piper colubrinum Link, a distant relative of cultivated black pepper, is highly resistant to this destructive pathogen. Identification of resistance (R) genes in P. colubrinum and the study of its expression profile during interaction with the pathogen can help in understanding the resistance mechanism involved. In the present study, 1289 R gene-related transcripts were mined from P. colubrinum transcriptome, clustered, and classified according to the conserved motifs and domains. Transcripts belonging to four major R gene classes were identified in P. colubrinum, but TIR-NBS-LRR-type R genes were absent. The relative expression of 12 selected R genes was studied using two virulent isolates of P. capsici, and these were found to be upregulated in the initial hours of plant pathogen interaction. The R genes studied were expressed even in aseptically maintained tissue-cultured plants and uninoculated greenhouse-grown plants at basal level suggesting that the plants are geared up with the R gene all the time and are under continuous surveillance for the pathogen and basal level of R gene expression do not require a pathogen trigger. ACT, ATUB, and EIF3E were identified as the most stable reference genes that can be used for real-time PCR study. The present study identified promising R genes in P. colubrinum which can be used in developing Phytophthora-resistant black pepper.

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.

Prediction and analysis of three gene families related to leaf rust (Puccinia triticina) resistance in wheat (Triticum aestivum L.)

BACKGROUND

The resistance to leaf rust (Lr) caused by Puccinia triticina in wheat (Triticum aestivum L.) has been well studied over the past decades with over 70 Lr genes being mapped on different chromosomes and numerous QTLs (quantitative trait loci) being detected or mapped using DNA markers. Such resistance is often divided into race-specific and race-nonspecific resistance. The race-nonspecific resistance can be further divided into resistance to most or all races of the same pathogen and resistance to multiple pathogens. At the molecular level, these three types of resistance may cover across the whole spectrum of pathogen specificities that are controlled by genes encoding different protein families in wheat. The objective of this study is to predict and analyze genes in three such families: NBS-LRR (nucleotide-binding sites and leucine-rich repeats or NLR), START (Steroidogenic Acute Regulatory protein [STaR] related lipid-transfer) and ABC (ATP-Binding Cassette) transporter. The focus of the analysis is on the patterns of relationships between these protein-coding genes within the gene families and QTLs detected for leaf rust resistance.

RESULTS

We predicted 526 ABC, 1117 NLR and 144 START genes in the hexaploid wheat genome through a domain analysis of wheat proteome. Of the 1809 SNPs from leaf rust resistance QTLs in seedling and adult stages of wheat, 126 SNPs were found within coding regions of these genes or their neighborhood (5 Kb upstream from transcription start site [TSS] or downstream from transcription termination site [TTS] of the genes). Forty-three of these SNPs for adult resistance and 18 SNPs for seedling resistance reside within coding or neighboring regions of the ABC genes whereas 14 SNPs for adult resistance and 29 SNPs for seedling resistance reside within coding or neighboring regions of the NLR gene. Moreover, we found 17 nonsynonymous SNPs for adult resistance and five SNPs for seedling resistance in the ABC genes, and five nonsynonymous SNPs for adult resistance and six SNPs for seedling resistance in the NLR genes. Most of these coding SNPs were predicted to alter encoded amino acids and such information may serve as a starting point towards more thorough molecular and functional characterization of the designated Lr genes. Using the primer sequences of 99 known non-SNP markers from leaf rust resistance QTLs, we found candidate genes closely linked to these markers, including Lr34 with distances to its two gene-specific markers being 1212 bases (to cssfr1) and 2189 bases (to cssfr2).

CONCLUSION

This study represents a comprehensive analysis of ABC, NLR and START genes in the hexaploid wheat genome and their physical relationships with QTLs for leaf rust resistance at seedling and adult stages. Our analysis suggests that the ABC (and START) genes are more likely to be co-located with QTLs for race-nonspecific, adult resistance whereas the NLR genes are more likely to be co-located with QTLs for race-specific resistance that would be often expressed at the seedling stage. Though our analysis was hampered by inaccurate or unknown physical positions of numerous QTLs due to the incomplete assembly of the complex hexaploid wheat genome that is currently available, the observed associations between (i) QTLs for race-specific resistance and NLR genes and (ii) QTLs for nonspecific resistance and ABC genes will help discover SNP variants for leaf rust resistance at seedling and adult stages. The genes containing nonsynonymous SNPs are promising candidates that can be investigated in future studies as potential new sources of leaf rust resistance in wheat breeding.

Genome analysis of Hibiscus syriacus provides insights of polyploidization and indeterminate flowering in woody plants

Hibiscus syriacus (L.) (rose of Sharon) is one of the most widespread garden shrubs in the world. We report a draft of the H. syriacus genome comprised of a 1.75 Gb assembly that covers 92% of the genome with only 1.7% (33 Mb) gap sequences. Predicted gene modeling detected 87,603 genes, mostly supported by deep RNA sequencing data. To define gene family distribution among relatives of H. syriacus, orthologous gene sets containing 164,660 genes in 21,472 clusters were identified by OrthoMCL analysis of five plant species, including H. syriacus, Arabidopsis thaliana, Gossypium raimondii, Theobroma cacao and Amborella trichopoda. We inferred their evolutionary relationships based on divergence times among Malvaceae plant genes and found that gene families involved in flowering regulation and disease resistance were more highly divergent and expanded in H. syriacus than in its close relatives, G. raimondii (DD) and T. cacao. Clustered gene families and gene collinearity analysis revealed that two recent rounds of whole-genome duplication were followed by diploidization of the H. syriacus genome after speciation. Copy number variation and phylogenetic divergence indicates that WGDs and subsequent diploidization led to unequal duplication and deletion of flowering-related genes in H. syriacus and may affect its unique floral morphology.

RGAugury: a pipeline for genome-wide prediction of resistance gene analogs (RGAs) in plants

Background

Resistance gene analogs (RGAs), such as NBS-encoding proteins, receptor-like protein kinases (RLKs) and receptor-like proteins (RLPs), are potential R-genes that contain specific conserved domains and motifs. Thus, RGAs can be predicted based on their conserved structural features using bioinformatics tools. Computer programs have been developed for the identification of individual domains and motifs from the protein sequences of RGAs but none offer a systematic assessment of the different types of RGAs. A user-friendly and efficient pipeline is needed for large-scale genome-wide RGA predictions of the growing number of sequenced plant genomes.

Results

An integrative pipeline, named RGAugury, was developed to automate RGA prediction. The pipeline first identifies RGA-related protein domains and motifs, namely nucleotide binding site (NB-ARC), leucine rich repeat (LRR), transmembrane (TM), serine/threonine and tyrosine kinase (STTK), lysin motif (LysM), coiled-coil (CC) and Toll/Interleukin-1 receptor (TIR). RGA candidates are identified and classified into four major families based on the presence of combinations of these RGA domains and motifs: NBS-encoding, TM-CC, and membrane associated RLP and RLK. All time-consuming analyses of the pipeline are paralleled to improve performance. The pipeline was evaluated using the well-annotated Arabidopsis genome. A total of 98.5, 85.2, and 100 % of the reported NBS-encoding genes, membrane associated RLPs and RLKs were validated, respectively. The pipeline was also successfully applied to predict RGAs for 50 sequenced plant genomes. A user-friendly web interface was implemented to ease command line operations, facilitate visualization and simplify result management for multiple datasets.

Conclusions

RGAugury is an efficiently integrative bioinformatics tool for large scale genome-wide identification of RGAs. It is freely available at Bitbucket: https://bitbucket.org/yaanlpc/rgaugury.

Electronic supplementary material

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

A resistance locus in the American heirloom rice variety Carolina Gold Select is triggered by TAL effectors with diverse predicted targets and is effective against African strains of Xanthomonas oryzae pv. oryzicola

The rice pathogens Xanthomonas oryzae pathovar (pv.) oryzae and pv. oryzicola produce numerous transcription activator-like (TAL) effectors that increase bacterial virulence by activating expression of host susceptibility genes. Rice resistance mechanisms against TAL effectors include polymorphisms that prevent effector binding to susceptibility gene promoters, or that allow effector activation of resistance genes. This study identifies, in the heirloom variety Carolina Gold Select, a third mechanism of rice resistance involving TAL effectors. This resistance manifests through strong suppression of disease development in response to diverse TAL effectors from both X. oryzae pathovars. The resistance can be triggered by an effector with only 3.5 central repeats, is independent of the composition of the repeat variable di-residues that determine TAL effector binding specificity, and is independent of the transcriptional activation domain. We determined that the resistance is conferred by a single dominant locus, designated Xo1, that maps to a 1.09 Mbp fragment on chromosome 4. The Xo1 interval also confers complete resistance to the strains in the African clade of X. oryzae pv. oryzicola, representing the first dominant resistance locus against bacterial leaf streak in rice. The strong phenotypic similarity between the TAL effector-triggered resistance conferred by Xo1 and that conferred by the tomato resistance gene Bs4 suggests that monocots and dicots share an ancient or convergently evolved mechanism to recognize analogous TAL effector epitopes.

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

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

TNL genes in peach: insights into the post-LRR domain

Background

Plants develop sustainable defence responses to pathogen attacks through resistance (R) genes contributing to effector-triggered immunity (ETI). TIR-NB-LRR genes (TNL genes) constitute a major family of ETI R genes in dicots. The putative functions or roles of the TIR, NB and LRR domains of the proteins they encode (TNLs) are well documented, but TNLs also have a poorly characterised C-terminal region, the function of which is unknown in most cases. We characterised this prevalent stress-response protein family in a perennial plant, using the genome of peach (Prunus persica), the model Prunus species. The first TNL gene from this genus to be cloned, the Ma gene, confers complete-spectrum resistance to root-knot nematodes (RKNs) and encodes a protein with a huge C-terminal region with five duplicated post-LRR (PL) domains. This gene was the cornerstone of this study.

Results

We investigated the role of this C-terminal region, by first describing the frequency, distribution and structural characteristics of i) TNL genes and ii) their PL domains in the peach genome, using the v1.0 Sanger sequence together with the v2.0 sequence, which has better genome annotation due to the incorporation of transcriptomic data. We detected 195 predicted TNL genes from the eight peach chromosomes: 85 % of these genes mapped to chromosomes 1, 2, 7 and 8. We reconstructed the putative structure of the predicted exons of all the TNL genes identified, and it was possible to retrieve the PL domains among two thirds of the TNL genes. We used our predicted TNL gene sequences to develop an annotation file for use with the Gbrowse tool in the v2.0 genome. The use of these annotation data made it possible to detect transcribed PL sequences in two Prunus species. We then used consensus sequences defined on the basis of 124 PL domains to design specific motifs, and we found that the use of these motifs significantly increased the numbers of PL domains and correlative TNL genes detected in diverse dicot genomes. Based on PL signatures, we showed that TNL genes with multiple PL domains were rare in peach and the other plants screened. The five-PL domain pattern is probably unique to Ma and its orthologues within Prunus and closely related genera from the Rosaceae and was probably inherited from the common ancestor of these plants in the subfamily Spiraeoideae.

Conclusions

The first physical TNL gene map for Prunus species can be used for the further investigation of R genes in this genus. The PL signature motifs are a complementary tool for the detection of TNL R genes in dicots. The low degree of similarity between PL domains and the neighbouring LRR exons and the specificity of PL signature motifs suggest that PL and LRR domains have different origins, with PL domains being specific to TNL genes, and possibly essential to the functioning of these genes in some cases. Investigations of the role of the oversized Ma PL region, in ligand binding or intramolecular interactions for example, may help to enrich our understanding of NB-LRR-mediated plant immunity to RKNs.

Electronic supplementary material

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

Large-Scale Analyses of Angiosperm Nucleotide-Binding Site-Leucine-Rich Repeat Genes Reveal Three Anciently Diverged Classes with Distinct Evolutionary Patterns

Nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes make up the largest plant disease resistance gene family (R genes), with hundreds of copies occurring in individual angiosperm genomes. However, the expansion history of NBS-LRR genes during angiosperm evolution is largely unknown. By identifying more than 6,000 NBS-LRR genes in 22 representative angiosperms and reconstructing their phylogenies, we present a potential framework of NBS-LRR gene evolution in the angiosperm. Three anciently diverged NBS-LRR classes (TNLs, CNLs, and RNLs) were distinguished with unique exon-intron structures and DNA motif sequences. A total of seven ancient TNL, 14 CNL, and two RNL lineages were discovered in the ancestral angiosperm, from which all current NBS-LRR gene repertoires were evolved. A pattern of gradual expansion during the first 100 million years of evolution of the angiosperm clade was observed for CNLs. TNL numbers remained stable during this period but were eventually deleted in three divergent angiosperm lineages. We inferred that an intense expansion of both TNL and CNL genes started from the Cretaceous-Paleogene boundary. Because dramatic environmental changes and an explosion in fungal diversity occurred during this period, the observed expansions of R genes probably reflect convergent adaptive responses of various angiosperm families. An ancient whole-genome duplication event that occurred in an angiosperm ancestor resulted in two RNL lineages, which were conservatively evolved and acted as scaffold proteins for defense signal transduction. Overall, the reconstructed framework of angiosperm NBS-LRR gene evolution in this study may serve as a fundamental reference for better understanding angiosperm NBS-LRR genes.

Uncovering the dynamic evolution of nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes in Brassicaceae

Plant genomes harbor dozens to hundreds of nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes; however, the long-term evolutionary history of these resistance genes has not been fully understood. This study focuses on five Brassicaceae genomes and the Carica papaya genome to explore changes in NBS-LRR genes that have taken place in this Rosid II lineage during the past 72 million years. Various numbers of NBS-LRR genes were identified from Arabidopsis lyrata (198), A. thaliana (165), Brassica rapa (204), Capsella rubella (127), Thellungiella salsuginea (88), and C. papaya (51). In each genome, the identified NBS-LRR genes were found to be unevenly distributed among chromosomes and most of them were clustered together. Phylogenetic analysis revealed that, before and after Brassicaceae speciation events, both toll/interleukin-1 receptor-NBS-LRR (TNL) genes and non-toll/interleukin-1 receptor-NBS-LRR (nTNL) genes exhibited a pattern of first expansion and then contraction, suggesting that both subclasses of NBS-LRR genes were responding to pathogen pressures synchronically. Further, by examining the gain/loss of TNL and nTNL genes at different evolutionary nodes, this study revealed that both events often occurred more drastically in TNL genes. Finally, the phylogeny of nTNL genes suggested that this NBS-LRR subclass is composed of two separate ancient gene types: RPW8-NBS-LRR and Coiled-coil-NBS-LRR.

The Eucalyptus grandis NBS-LRR Gene Family: Physical Clustering and Expression Hotspots

Eucalyptus grandis is a commercially important hardwood species and is known to be susceptible to a number of pests and pathogens. Determining mechanisms of defense is therefore a research priority. The published genome for E. grandis has aided the identification of one important class of resistance (R) genes that incorporate nucleotide binding sites and leucine-rich repeat domains (NBS-LRR). Using an iterative search process we identified NBS-LRR gene models within the E. grandis genome. We characterized the gene models and identified their genomic arrangement. The gene expression patterns were examined in E. grandis clones, challenged with a fungal pathogen (Chrysoporthe austroafricana) and insect pest (Leptocybe invasa). One thousand two hundred and fifteen putative NBS-LRR coding sequences were located which aligned into two large classes, Toll or interleukin-1 receptor (TIR) and coiled-coil (CC) based on NB-ARC domains. NBS-LRR gene-rich regions were identified with 76% organized in clusters of three or more genes. A further 272 putative incomplete resistance genes were also identified. We determined that E. grandis has a higher ratio of TIR to CC classed genes compared to other woody plant species as well as a smaller percentage of single NBS-LRR genes. Transcriptome profiles indicated expression hotspots, within physical clusters, including expression of many incomplete genes. The clustering of putative NBS-LRR genes correlates with differential expression responses in resistant and susceptible plants indicating functional relevance for the physical arrangement of this gene family. This analysis of the repertoire and expression of E. grandis putative NBS-LRR genes provides an important resource for the identification of novel and functional R-genes; a key objective for strategies to enhance resilience.

Identification of Immune Related LRR-Containing Genes in Maize (Zea mays L.) by Genome-Wide Sequence Analysis

A large number of immune receptors consist of nucleotide binding site-leucine rich repeat (NBS-LRR) proteins and leucine rich repeat-receptor-like kinases (LRR-RLK) that play a crucial role in plant disease resistance. Although many NBS-LRR genes have been previously identified in Zea mays, there are no reports on identifying NBS-LRR genes encoded in the N-terminal Toll/interleukin-1 receptor (TIR) motif and identifying genome-wide LRR-RLK genes. In the present study, 151 NBS-LRR genes and 226 LRR-RLK genes were identified after performing bioinformatics analysis of the entire maize genome. Of these identified genes, 64 NBS-LRR genes and four TIR-NBS-LRR genes were identified for the first time. The NBS-LRR genes are unevenly distributed on each chromosome with gene clusters located at the distal end of each chromosome, while LRR-RLK genes have a random chromosomal distribution with more paired genes. Additionally, six LRR-RLK/RLPs including FLS2, PSY1R, PSKR1, BIR1, SERK3, and Cf5 were characterized in Zea mays for the first time. Their predicted amino acid sequences have similar protein structures with their respective homologues in other plants, indicating that these maize LRR-RLK/RLPs have the same functions as their homologues act as immune receptors. The identified gene sequences would assist in the study of their functions in maize.

Arabidopsis BNT1, an atypical TIR-NBS-LRR gene, acting as a regulator of the hormonal response to stress

During their life cycle, plants have to cope with fluctuating environmental conditions. The perception of the stressful environmental conditions induces a specific stress hormone signature specifying a proper response with an efficient fitness. By reverse genetics, we isolated and characterized a novel mutation in Arabidopsis, associated with environmental stress responses, that affects the At5g11250/BURNOUT1 (BNT1) gene which encode a Toll/Interleukin1 receptor-nucleotide binding site leucine-rich repeat (TIR-NBS-LRR) protein. The knock-out bnt1 mutants displayed, in the absence of stress conditions, a multitude of growth and development defects, suchas severe dwarfism, early senescence and flower sterility, similar to those observed in vitro in wild type plants upon different biotic and/or abiotic stresses. The disruption of BNT1 causes also a drastic increase of the jasmonic, salicylic and abscisic acids as well as ethylene levels. Which was consistent with the expression pattern observed in bnt1 showing an over representation of genes involved in the hormonal response to stress? Therefore, a defect in BNT1 forced the plant to engage in an exhausting general stress response, which produced frail, weakened and poorly adapted plants expressing "burnout" syndromes. Furthermore, by in vitro phenocopying experiments, physiological, chemical and molecular analyses, we propose that BNT1 could represent a molecular link between stress perception and specific hormonal signature.

Disease Resistance Gene Analogs (RGAs) in Plants

Plants have developed effective mechanisms to recognize and respond to infections caused by pathogens. Plant resistance gene analogs (RGAs), as resistance (R) gene candidates, have conserved domains and motifs that play specific roles in pathogens' resistance. Well-known RGAs are nucleotide binding site leucine rich repeats, receptor like kinases, and receptor like proteins. Others include pentatricopeptide repeats and apoplastic peroxidases. RGAs can be detected using bioinformatics tools based on their conserved structural features. Thousands of RGAs have been identified from sequenced plant genomes. High-density genome-wide RGA genetic maps are useful for designing diagnostic markers and identifying quantitative trait loci (QTL) or markers associated with plant disease resistance. This review focuses on recent advances in structures and mechanisms of RGAs, and their identification from sequenced genomes using bioinformatics tools. Applications in enhancing fine mapping and cloning of plant disease resistance genes are also discussed.

Identification and distribution of the NBS-LRR gene family in the Cassava genome

BACKGROUND

Plant resistance genes (R genes) exist in large families and usually contain both a nucleotide-binding site domain and a leucine-rich repeat domain, denoted NBS-LRR. The genome sequence of cassava (Manihot esculenta) is a valuable resource for analysing the genomic organization of resistance genes in this crop.

RESULTS

With searches for Pfam domains and manual curation of the cassava gene annotations, we identified 228 NBS-LRR type genes and 99 partial NBS genes. These represent almost 1% of the total predicted genes and show high sequence similarity to proteins from other plant species. Furthermore, 34 contained an N-terminal toll/interleukin (TIR)-like domain, and 128 contained an N-terminal coiled-coil (CC) domain. 63% of the 327 R genes occurred in 39 clusters on the chromosomes. These clusters are mostly homogeneous, containing NBS-LRRs derived from a recent common ancestor.

CONCLUSIONS

This study provides insight into the evolution of NBS-LRR genes in the cassava genome; the phylogenetic and mapping information may aid efforts to further characterize the function of these predicted R genes.