Genome-wide identification of NBS-encoding resistance genes in Brassica rapa

Disease resistance gene count increases with rainfall in Silphium integrifolium

Intracellular plant defense against pathogens is mediated by a class of disease resistance genes known as NB-LRRs or NLRs (R genes). Many of the diseases these genes protect against are more prevalent in regions of higher rainfall, which provide better growth conditions for the pathogens. As such, we expect a higher selective pressure for the maintenance and proliferation of R genes in plants adapted to wetter conditions. In this study, we enriched libraries for R genes using RenSeq from baits primarily developed from the common sunflower (Helianthus annuus) reference genome. We sequenced the R gene libraries of Silphium integrifolium Michx, a perennial relative of sunflower, from 12 prairie remnants across a rainfall gradient in the Central Plains of the United States, with both Illumina short-read (n = 99) and PacBio long-read (n = 10) approaches. We found a positive relationship between the mean effective annual precipitation of a plant's source prairie remnant and the number of R genes in its genome, consistent with intensity of plant pathogen coevolution increasing with precipitation. We show that RenSeq can be applied to the study of ecological hypotheses in non-model relatives of model organisms.

In silico structural-functional characterization of three differentially expressed resistance gene analogs identified in Dalbergia sissoo against dieback disease reveals their role in immune response regulation

Plant immunity includes enemy recognition, signal transduction, and defensive response against pathogens. We experimented to identify the genes that contribute resistance against dieback disease to Dalbergia sissoo, an economically important timber tree. In this study, we investigated the role of three differentially expressed genes identified in the dieback-induced transcriptome in Dalbergia sissoo. The transcriptome was probed using DOP-rtPCR analysis. The identified RGAs were characterized in silico as the contributors of disease resistance that switch on under dieback stress. Their predicted fingerprints revealed involvement in stress response. Ds-DbRCaG-02-Rga.a, Ds-DbRCaG-04-Rga.b, and Ds-DbRCaG-06-Rga.c showed structural homology with the Transthyretin-52 domain, EAL associated YkuI_C domain, and Src homology-3 domain respectively, which are the attributes of signaling proteins possessing a role in regulating immune responses in plants. Based on in-silico structural and functional characterization, they were predicted to have a role in immune response regulation in D. sissoo.

Bulk segregant analysis coupled with transcriptomics and metabolomics revealed key regulators of bacterial leaf blight resistance in rice

BACKGROUND

Bacterial leaf blight (BLB) is a highly destructive disease, causing significant yield losses in rice (Oryza sativa). Genetic variation is contemplated as the most effective measure for inducing resistance in plants. The mutant line T1247 derived from R3550 (BLB susceptible) was highly resistant to BLB. Therefore, by utilizing this valuable source, we employed bulk segregant analysis (BSA) and transcriptome profiling to identify the genetic basis of BLB resistance in T1247.

RESULTS

The differential subtraction method in BSA identified a quantitative trait locus (QTL) on chromosome 11 spanning a 27-27.45 Mb region with 33 genes and 4 differentially expressed genes (DEGs). Four DEGs (P < 0.01) with three putative candidate genes, OsR498G1120557200, OsR498G1120555700, and OsR498G1120563600,0.01 in the QTL region were identified with specific regulation as a response to BLB inoculation. Moreover, transcriptome profiling identified 37 resistance analogs genes displaying differential regulation.

CONCLUSIONS

Our study provides a substantial addition to the available information regarding QTLs associated with BLB, and further functional verification of identified candidate genes can broaden the scope of understanding the BLB resistance mechanism in rice.

Genome-wide identification and analysis of the NLR gene family in Medicago ruthenica

Medicago ruthenica, important forage in the legume family, possesses high nutritional value and carries abundant tolerance genes. This study used whole-genome data of M. ruthenica to perform a genome-wide analysis of the nucleotide-binding site-leucine-rich repeat receptor (NLR) gene family, which is the largest family of plant disease resistance genes (R genes). A total of 338 NLR genes were identified in the M. ruthenica genome, including 160 typical genes that contained 80 coiled-coil (CC)-NBS-LRR (CNL) genes, 76 toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) genes, four resistance to powdery mildew 8 (RPW8)-NBS-LRR (RNL) subclass genes, and 178 atypical NLR genes encoding proteins without at least one important domain. Among its eight chromosomes, M. ruthenica chromosomes 3 and 8 contained most of the NLR genes. More than 40% of all NLR genes were located on these two chromosomes, mainly in multigene clusters. The NLR proteins of M. ruthenica had six highly conserved motifs: P-loop, GLPL, RNBS-D, kinase-2, RNBS-C, and MHDV. Phylogenetic analysis revealed that the NLR genes of M. ruthenica formed three deeply separated clades according to the N-terminal domain of the proteins encoded by these genes. Gene duplication and syntenic analysis suggested four gene duplication types in the NLR genes of M. ruthenica, namely, tandem, proximal, dispersed, and segmental duplicates, which involved 189, 49, 59, and 41 genes, respectively. A total of 41 segmental duplication genes formed 23 NLR gene pairs located on syntenic chromosomal blocks mainly between chromosomes 6 and 7. In addition, syntenic analysis between M. truncatula and M. ruthenica revealed 193 gene pairs located on syntenic chromosomal blocks of the two species. The expression analysis of M. ruthenica NLR genes showed that 303 (89.6%) of the NLR genes were expressed in different varieties. Overall, this study described the full NLR profile of the M. ruthenica genome to provide an important resource for mining disease-resistant genes and disease-resistant breeding.

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.

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

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

Editing of a Novel Cd Uptake-Related Gene CUP1 Contributes to Reducing Cd Accumulations in Arabidopsis thaliana and Brassica napus

Brassica napus is a Cd hyperaccumulator, which is a serious threat to food and fodder safety. However, no related studies on developing Cd-safe B. napus have been reported yet. Here, we screened out a novel Cd uptake-related gene, AtCUP1, from the major facilitator superfamily in Arabidopsis thaliana. The mutation of AtCUP1 decreased Cd accumulation, both in roots and shoots of A. thaliana. Furthermore, the disruption of the AtCUP1 gene by the CRISPR/Cas9 system significantly reduced Cd accumulation in A. thaliana. Interestingly, the disruption of the BnCUP1 gene, an orthologous gene of AtCUP1, by the CRISPR/Cas9 system also diminished Cd accumulation in both roots and shoots of B. napus based on the hydroponics assay. Furthermore, for the field experiment, the Cd accumulations of BnCUP1-edited lines were reduced by 52% in roots and 77% in shoots compared to that of wild-type (WT) lines, and the biomass and yield of BnCUP1-edited lines increased by 42% and 47% of that of WT, respectively. Noteworthily, agronomic characteristics of B. napus were not apparently affected by BnCUP1-editing. Thus, BnCUP1-edited lines are excellent non-transgenic germplasm resources for reducing Cd accumulation without a distinct compromise in yield, which could be applied to agricultural production in Cd-contaminated soils.

Genome-wide identification and characterization of NBS-encoding genes in the sweet potato wild ancestor Ipomoea trifida (H.B.K.)

The most predominant type of resistance (R) genes contain nucleotide-binding sites and leucine-rich repeat (NBS-LRR) domains, characterization of which is helpful for plant resistance improvement. However, the NBS genes of Ipomoea trifida (H.B.K.) remain insufficient to date. In this study, a genome-wide analysis of the NBS-encoding gene in I. trifida (H.B.K.) was carried out. A total of 442 NBS encoding genes were identified, amounting to 1.37% of the total genes of I. trifida (H.B.K.). Based on the analysis of the domains, the identified ItfNBS genes were further classified into seven groups: CNL, NL, CN, N, TNL, TN, and RNL. Phylogenetic analysis showed that the I. trifida NBS genes clustered into three independent clades: RNL, TNL, and CNL. Chromosome location analysis revealed that the distribution of ItfNBS genes in chromosomes was uneven, with a number ranging from 3 to 45. Multiple stress-related regulatory elements were detected in the promoters of the NBS-encoding genes, and their expression profiles were obtained. The qRT-PCR analysis revealed that IbNBS10, IbNBS20, IbNBS258, and IbNBS88 responded to stem nematode infection. These results provide critical proof for further characterization and analysis of NBS-encoding genes with important functions.

The Arabidopsis effector-triggered immunity landscape is conserved in oilseed crops

The bacterial phytopathogen Pseudomonas syringae causes disease on a wide array of plants, including the model plant Arabidopsis thaliana and its agronomically important relatives in the Brassicaceae family. To cause disease, P. syringae delivers effector proteins into plant cells through a type III secretion system. In response, plant nucleotide-binding leucine-rich repeat proteins recognize specific effectors and mount effector-triggered immunity (ETI). While ETI is pervasive across A. thaliana, with at least 19 families of P. syringae effectors recognized in this model species, the ETI landscapes of crop species have yet to be systematically studied. Here, we investigated the conservation of the A. thaliana ETI landscape in two closely related oilseed crops, Brassica napus (canola) and Camelina sativa (false flax). We show that the level of immune conservation is inversely related to the degree of evolutionary divergence from A. thaliana, with the more closely related C. sativa losing ETI responses to only one of the 19 P. syringae effectors tested, while the more distantly related B. napus loses ETI responses to four effectors. In contrast to the qualitative conservation of immune response, the quantitative rank order is not as well-maintained across the three species and diverges increasingly with evolutionary distance from A. thaliana. Overall, our results indicate that the A. thaliana ETI profile is qualitatively conserved in oilseed crops, but quantitatively distinct.

Identification and Characterization of NBS Resistance Genes in Akebia trifoliata

Akebia trifoliata is an important multiuse perennial plant that often suffers attacks from various pathogens due to its long growth cycle, seriously affecting its commercial value. The absence of research on the resistance (R) genes of A. trifoliata has greatly limited progress in the breeding of resistant varieties. Genes encoding proteins containing nucleotide binding sites (NBSs) and C-terminal leucine-rich repeats (LRRs), the largest family of plant resistance (R) genes, are vital for plant disease resistance. A comprehensive genome-wide analysis showed that there were only 73 NBS genes in the A. trifoliata genome, including three main subfamilies (50 coiled coil (CC)-NBS-LRR (CNL), 19 Toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) and four resistance to powdery mildew8 (RPW8)-NBS-LRR (RNL) genes). Additionally, 64 mapped NBS candidates were unevenly distributed on 14 chromosomes, most of which were assigned to the chromosome ends; 41 of these genes were located in clusters, and the remaining 23 genes were singletons. Both the CNLs and TNLs were further divided into four subgroups, and the CNLs had fewer exons than the TNLs. Structurally, all eight previously reported conserved motifs were identified in the NBS domains, and both their order and their amino acid sequences exhibited high conservation. Evolutionarily, tandem and dispersed duplications were shown to be the two main forces responsible for NBS expansion, producing 33 and 29 genes, respectively. A transcriptome analysis of three fruit tissues at four developmental stages showed that NBS genes were generally expressed at low levels, while a few of these genes showed relatively high expression during later development in rind tissues. Overall, this research is the first to identify and characterize A. trifoliata NBS genes and is valuable for both the development of new resistant cultivars and the study of molecular mechanisms of resistance.

Characterization of Nucleotide Binding -Site-Encoding Genes in Sweetpotato, Ipomoea batatas(L.) Lam., and Their Response to Biotic and Abiotic Stresses

Sweetpotato, Ipomoea batatas (L.) Lam., is an important and widely grown crop, yet its production is affected severely by biotic and abiotic stresses. The nucleotide binding site (NBS)-encoding genes have been shown to improve stress tolerance in several plant species. However, the characterization of NBS-encoding genes in sweetpotato is not well-documented to date. In this study, a comprehensive analysis of NBS-encoding genes has been conducted on this species by using bioinformatics and molecular biology methods. A total of 315 NBS-encoding genes were identified, and 260 of them contained all essential conserved domains while 55 genes were truncated. Based on domain architectures, the 260 NBS-encoding genes were grouped into 6 distinct categories. Phylogenetic analysis grouped these genes into 3 classes: TIR, CC (I), and CC (II). Chromosome location analysis revealed that the distribution of NBS-encoding genes in chromosomes was uneven, with a number ranging from 1 to 34. Multiple stress-related regulatory elements were detected in the promoters, and the NBS-encoding genes' expression profiles under biotic and abiotic stresses were obtained. According to the bioinformatics analysis, 9 genes were selected for RT-qPCR analysis. The results revealed that IbNBS75, IbNBS219, and IbNBS256 respond to stem nematode infection; Ib-NBS240, IbNBS90, and IbNBS80 respond to cold stress, while IbNBS208, IbNBS71, and IbNBS159 respond to 30% PEG treatment. We hope these results will provide new insights into the evolution of NBS-encoding genes in the sweetpotato genome and contribute to the molecular breeding of sweetpotato in the future.

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.

Identification of resistance loci against new pathotypes of Plasmodiophora brassicae in Brassica napus based on genome-wide association mapping

Genetic resistance is a successful strategy for management of clubroot (Plasmodiophora brassicae) of brassica crops, but resistance can break down quickly. Identification of novel sources of resistance is especially important when new pathotypes arise. In the current study, the reaction of 177 accessions of Brassica napus to four new, virulent pathotypes of P. brassicae was assessed. Each accession was genotyped using genotyping by sequencing to identify and map novel sources of clubroot resistance using mixed linear model (MLM) analysis. The majority of accessions were highly susceptible (70–100 DSI), but a few accessions exhibited strong resistance (0–20 DSI) to pathotypes 5X (21 accessions), 3A (8), 2B (7), and 3D (15), based on the Canadian Clubroot Differential system. In total, 301,753 SNPs were mapped to 19 chromosomes. Population structure analysis indicated that the 177 accessions belong to seven major populations. SNPs were associated with resistance to each pathotype using MLM. In total, 13 important SNP loci were identified, with 9 SNPs mapped to the A-genome and 4 to the C-genome. The SNPs were associated with resistance to pathotypes 5X (2 SNPs), 3A (4), 2B (5) and 3D (6). A Blast search of 1.6 Mb upstream and downstream from each SNP identified 13 disease-resistance genes or domains. The distance between a SNP locus and the nearest resistance gene ranged from 0.04 to 0.74 Mb. The resistant lines and SNP markers identified in this study can be used to breed for resistance to the most prevalent new pathotypes of P. brassicae in Canada.

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

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

Genome-wide identification and characterization of NBS-encoding genes in Raphanus sativus L. and their roles related to Fusarium oxysporum resistance

BACKGROUND

The nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes are important for plant development and disease resistance. Although genome-wide studies of NBS-encoding genes have been performed in several species, the evolution, structure, expression, and function of these genes remain unknown in radish (Raphanus sativus L.). A recently released draft R. sativus L. reference genome has facilitated the genome-wide identification and characterization of NBS-encoding genes in radish.

RESULTS

A total of 225 NBS-encoding genes were identified in the radish genome based on the essential NB-ARC domain through HMM search and Pfam database, with 202 mapped onto nine chromosomes and the remaining 23 localized on different scaffolds. According to a gene structure analysis, we identified 99 NBS-LRR-type genes and 126 partial NBS-encoding genes. Additionally, 80 and 19 genes respectively encoded an N-terminal Toll/interleukin-like domain and a coiled-coil domain. Furthermore, 72% of the 202 NBS-encoding genes were grouped in 48 clusters distributed in 24 crucifer blocks on chromosomes. The U block on chromosomes R02, R04, and R08 had the most NBS-encoding genes (48), followed by the R (24), D (23), E (23), and F (17) blocks. These clusters were mostly homogeneous, containing NBS-encoding genes derived from a recent common ancestor. Tandem (15 events) and segmental (20 events) duplications were revealed in the NBS family. Comparative evolutionary analyses of orthologous genes among Arabidopsis thaliana, Brassica rapa, and Brassica oleracea reflected the importance of the NBS-LRR gene family during evolution. Moreover, examinations of cis-elements identified 70 major elements involved in responses to methyl jasmonate, abscisic acid, auxin, and salicylic acid. According to RNA-seq expression analyses, 75 NBS-encoding genes contributed to the resistance of radish to Fusarium wilt. A quantitative real-time PCR analysis revealed that RsTNL03 (Rs093020) and RsTNL09 (Rs042580) expression positively regulates radish resistance to Fusarium oxysporum, in contrast to the negative regulatory role for RsTNL06 (Rs053740).

CONCLUSIONS

The NBS-encoding gene structures, tandem and segmental duplications, synteny, and expression profiles in radish were elucidated for the first time and compared with those of other Brassicaceae family members (A. thaliana, B. oleracea, and B. rapa) to clarify the evolution of the NBS gene family. These results may be useful for functionally characterizing NBS-encoding genes in radish.

Frontiers in Dissecting and Managing Brassica Diseases: From Reference-Based RGA Candidate Identification to Building Pan-RGAomes

The Brassica genus contains abundant economically important vegetable and oilseed crops, which are under threat of diseases caused by fungal, bacterial and viral pathogens. Resistance gene analogues (RGAs) are associated with quantitative and qualitative disease resistance and the identification of candidate RGAs associated with disease resistance is crucial for understanding the mechanism and management of diseases through breeding. The availability of Brassica genome assemblies has greatly facilitated reference-based quantitative trait loci (QTL) mapping for disease resistance. In addition, pangenomes, which characterise both core and variable genes, have been constructed for B. rapa, B. oleracea and B. napus. Genome-wide characterisation of RGAs using conserved domains and motifs in reference genomes and pangenomes reveals their clustered arrangements and presence of structural variations. Here, we comprehensively review RGA identification in important Brassica genome and pangenome assemblies. Comparison of the RGAs in QTL between resistant and susceptible individuals allows for efficient identification of candidate disease resistance genes. However, the reference-based QTL mapping and RGA candidate identification approach is restricted by the under-represented RGA diversity characterised in the limited number of Brassica assemblies. The species-wide repertoire of RGAs make up the pan-resistance gene analogue genome (pan-RGAome). Building a pan-RGAome, through either whole genome resequencing or resistance gene enrichment sequencing, would effectively capture RGA diversity, greatly expanding breeding resources that can be utilised for crop improvement.

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.

Genome-wide identification of the NPR1-like gene family in Brassica napus and functional characterization of BnaNPR1 in resistance to Sclerotinia sclerotiorum

The BnaNPR1-like gene family was identified in B. napus, and it was revealed that repression of BnaNPR1 significantly reduces resistance toS. sclerotiorum, intensifies ROS accumulation, and changes the expression of genes associated with SA and JA/ET signaling in response to this pathogen. The NONEXPRESSOR OF PATHOGENESIS-RELATED GENES1 (NPR1) and related NPR1-like genes play an important role in regulating plant defense. Oilseed rape (Brassica napus L.) is an important oilseed crop; however, little is known about the B. napus (Bna) NPR1-like gene family. Here, a total of 19 BnaNPR1-like genes were identified in the B. napus genome, and then named according to their respective best match in Arabidopsis thaliana (At), which led to the determination of B. napus homologs of every AtNPR1-like gene. Analysis of important protein domains and functional motifs indicated the conservation and variation among these homologs. Phylogenetic analysis of these BnaNPR1-like proteins and their Arabidopsis homologs revealed six distinct sub-clades, consequently indicating that their name classification totally conformed to their phylogenetic relationships. Further, B. napus transcriptomic data showed that the expression of three BnaNPR1s was significantly down-regulated in response to infection with Sclerotinia sclerotiorum, the most important pathogen of this crop, whereas BnaNPR2/3/4/5/6s did not show the expression differences in general. Further, we generated B. napus BnaNPR1-RNAi lines to interpret the effect of the down-regulated expression of BnaNPR1s on resistance to S. sclerotiorum. The results showed that BnaNPR1-RNAi significantly decreased this resistance. Further experiments revealed that BnaNPR1-RNAi intensified ROS production and changed defense responses in the interaction of plants with this pathogen. These results indicated that S. sclerotiorum might use BnaNPR1 to regulate specific physiological processes of B. napus, such as ROS production and SA defense response, for the infection.

Characterization of disease resistance genes in the Brassica napus pangenome reveals significant structural variation

Methods based on single nucleotide polymorphism (SNP), copy number variation (CNV) and presence/absence variation (PAV) discovery provide a valuable resource to study gene structure and evolution. However, as a result of these structural variations, a single reference genome is unable to cover the entire gene content of a species. Therefore, pangenomics analysis is needed to ensure that the genomic diversity within a species is fully represented. Brassica napus is one of the most important oilseed crops in the world and exhibits variability in its resistance genes across different cultivars. Here, we characterized resistance gene distribution across 50 B. napus lines. We identified a total of 1749 resistance gene analogs (RGAs), of which 996 are core and 753 are variable, 368 of which are not present in the reference genome (cv. Darmor-bzh). In addition, a total of 15 318 SNPs were predicted within 1030 of the RGAs. The results showed that core R-genes harbour more SNPs than variable genes. More nucleotide binding site-leucine-rich repeat (NBS-LRR) genes were located in clusters than as singletons, with variable genes more likely to be found in clusters. We identified 106 RGA candidates linked to blackleg resistance quantitative trait locus (QTL). This study provides a better understanding of resistance genes to target for genomics-based improvement and improved disease resistance.

Genomic Organization and Comparative Phylogenic Analysis of NBS-LRR Resistance Gene Family in Solanum pimpinellifolium and Arabidopsis thaliana

NBS-LRR (nucleotide-binding site and leucine-rich repeat) is one of the largest resistance gene families in plants. The completion of the genome sequencing of wild tomato Solanum pimpinellifolium provided an opportunity to conduct a comprehensive analysis of the NBS-LRR gene superfamily at the genome-wide level. In this study, gene identification, chromosome mapping, and phylogenetic analysis of the NBS-LRR gene family were analyzed using the bioinformatics methods. The results revealed 245 NBS-LRRs in total, similar to that in the cultivated tomato. These genes are unevenly distributed on 12 chromosomes, and ~59.6% of them form gene clusters, most of which are tandem duplications. Phylogenetic analysis divided the NBS-LRRs into 2 subfamilies (CNL-coiled-coil NBS-LRR and TNL-TIR NBS-LRR), and the expansion of the CNL subfamily was more extensive than the TNL subfamily. Novel conserved structures were identified through conserved motif analysis between the CNL and TNL subfamilies. Compared with the NBS-LRR sequences from the model plant Arabidopsis thaliana, wide genetic variation occurred after the divergence of S. pimpinellifolium and A thaliana. Species-specific expansion was also found in the CNL subfamily in S. pimpinellifolium. The results of this study provide the basis for the deeper analysis of NBS-LRR resistance genes and contribute to mapping and isolation of candidate resistance genes in S. pimpinellifolium.

DNA Methylation: Toward Crop Disease Resistance Improvement

Crop diseases, in conjunction with climate change, are a major threat to global crop production. DNA methylation is an epigenetic mark and is involved in plants' biological processes, including development, stress adaptation, and genome evolution. By providing a new source of variation, DNA methylation introduces novel direction to both scientists and breeders with its potential in disease resistance enhancement. Here, we discuss the impact of pathogen-induced DNA methylation modifications on a host's transcriptome reprogramming and genome stability, as part of the plant's defense mechanisms. We also highlight the knowledge gaps that need to be investigated for understanding the entire role of DNA methylation in plant pathogen interactions. This will ultimately assist breeders toward improving resistance and decreasing yield losses.

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.

Genome-Wide Analysis of Dynamin Gene Family in cassava (Manihot esculenta Crantz) and Transcriptional Regulation of Family Members ARC5 in Hormonal Treatments

The Dynamin gene family play a significance role in many physiological processes, especially ARC5 (Accumulation and replication of chloroplasts 5) in the process of plastid division. We performed a genome-wide analysis of the cassava Dynamin family based on the published cassava genome sequence and identified ARC5. 23 cassava Dynamins (MeDynamins) were identified and renamed. 23 MeDynamins were further divided into five major groups based on their structural and phylogenetic characteristics. The segmental duplication events have a significant impact on the expansion of MeDynamins. ARC5 expression analysis showed that there were differences between leaves and roots of cassava at different developmental stages. The tissue-specific expression analysis of the MeDynamins showed that most of MeDynamins were expressed in stem apical meristem and embryogenesis, whereas ARC5 was mainly expressed in leaves. The processing of IAA (Indole-3-acetic Acid) and MeJA (Methyl Jasmonate) verified the prediction results of cis-elements, and ACR5 was closely related to plant growth and positively correlated. It also indicated that high concentrations of MeJA treatment caused the cassava defense mechanism to function in advance. In conclusion, these findings provide basic insights for functional validation of the ARC5 genes in exogenous hormonal treatments.

The large repertoire of conifer NLR resistance genes includes drought responsive and highly diversified RNLs

The NLRs or NBS-LRRs (nucleotide-binding, leucine-rich-repeat) form the largest resistance gene family in plants, with lineage-specific contingents of TNL, CNL and RNL subfamilies and a central role in resilience to stress. The origin, evolution and distribution of NLR sequences has been unclear owing in part to the variable size and diversity of the RNL subfamily and a lack of data in Gymnosperms. We developed, searched and annotated transcriptomes assemblies of seven conifers and identified a resource of 3816 expressed NLR sequences. Our analyses encompassed sequences data spanning the major groups of land plants and determinations of NLR transcripts levels in response to drought in white spruce. We showed that conifers have among the most diverse and numerous RNLs in tested land plants. We report an evolutionary swap in the formation of RNLs, which emerged from the fusion of an RPW8 domain to a NB-ARC domain of CNL. We uncovered a quantitative relationship between RNLs and TNLs across all land plants investigated, with an average ratio of 1:10. The conifer RNL repertoire harbours four distinct groups, with two that differ from Angiosperms, one of which contained several upregulated sequences in response to drought while the majority of responsive NLRs are downregulated.

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

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

Genome-wide identification and characterization of Chitinase gene family in Brassica juncea and Camelina sativa in response to Alternaria brassicae

Chitinases belong to the group of Pathogenesis-related (PR) proteins that provides protection against fungal pathogens. This study presents the, genome-wide identification and characterization of chitinase gene family in two important oilseed crops B. juncea and C. sativa belonging to family Brassicaceae. We have identified 47 and 79 chitinase genes in the genomes of B. juncea and C. sativa, respectively. Phylogenetic analysis of chitinases in both the species revealed four distinct sub-groups, representing different classes of chitinases (I-V). Microscopic and biochemical study reveals the role of reactive oxygen species (ROS) scavenging enzymes in disease resistance of B. juncea and C. sativa. Furthermore, qRT-PCR analysis showed that expression of chitinases in both B. juncea and C. sativa was significantly induced after Alternaria brassicae infection. However, the fold change in chitinase gene expression was considerably higher in C. sativa compared to B. juncea, which further proves their role in C. sativa disease resistance to A. brassicae. This study provides comprehensive analysis on chitinase gene family in B. juncea and C. sativa and in future may serve as a potential candidate for improving disease resistance in B. juncea through transgenic approach.

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

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

Genome-wide characterization revealed role of NBS-LRR genes during powdery mildew infection in Vitis vinifera

NBS-LRR comprises a large class of disease resistance (R) proteins that play a widespread role in plant protection against pathogens. In grapevine, powdery mildew cause significant losses in its productivity and efforts are being directed towards finding of resistance loci or genes imparting resistance/tolerance against such fungal diseases. In the present study, we performed genome-wide analysis of NBS-LRR genes during PM infection in grapevine. We identified 18, 23, 12, 16, 10, 10, 9, 20 and 14 differentially expressed NBS-LRR genes in response to PM infection in seven partially PM-resistant (DVIT3351.27, Husseine, Karadzhandal, Khalchili, Late vavilov, O34-16, Sochal) and 2 PM-susceptible (Carignan and Thompson seedless) V. vinifera accessions. Further, the identified sequences were characterized based on chromosomal locations, physicochemical properties, gene structure and motif analysis, and functional annotation by Gene Ontology (GO) mapping. The NBS-LRR genes responsive to powdery mildew could potentially be exploited to improve resistance in grapes.

NBS-Encoding Genes in Brassica napus Evolved Rapidly After Allopolyploidization and Co-localize With Known Disease Resistance Loci

Genes containing nucleotide-binding sites (NBS) play an important role in pathogen resistance in plants. However, the evolutionary fate of NBS-encoding genes after formation of allotetraploid Brassica napus (A_nA_nC_nC_n, 2n = 38) is still unknown. We performed a genome-wide comparison of putatively functional NBS-encoding genes in B. napus and its progenitor species Brassica rapa (A_rA_r, 2n = 20) and Brassica oleracea (C_oC_o, 2n = 18), identifying 464, 202, and 146 putatively functional NBS-encoding genes respectively, with genes unevenly distributed in several clusters. The A_n-subgenome of B. napus possessed similar numbers of NBS-encoding genes (191 genes) to the A_r genome of B. rapa (202 genes) and similar clustering patterns. However, the C_n genome of B. napus had many more genes (273) than the B. oleracea C_o genome (146), with different clustering trends. Only 97 NBS-encoding genes (66.4%) in B. oleracea were homologous with NBS-encoding genes in B. napus, while 176 NBS-encoding genes (87.1%) were homologous between B. rapa and B. napus. These results suggest a greater diversification of NBS-encoding genes in the C genome may have occurred after formation of B. napus. Although most NBS-encoding genes in B. napus appeared to derive from the progenitors, the birth and death of several NBS-encoding genes was also putatively mediated by non-homologous recombination. The Ka/Ks values of most homologous pairs between B. napus and the progenitor species were less than 1, suggesting purifying selection during B. napus evolution. The majority of NBS-encoding genes (60% in all species) showed higher expression levels in root tissue (out of root, leaf, stem, seed and flower tissue types). Comparative analysis of NBS-encoding genes with mapped resistance QTL against three major diseases of B. napus (blackleg, clubroot and Sclerotinia stem rot) found 204 NBS-encoding genes in B. napus located within 71 resistance QTL intervals. The majority of NBS-encoding genes were co-located with resistance QTLs against a single disease, while 47 genes were co-located with QTLs against two diseases and 3 genes were co-located with QTLs against all three. Our results revealed significant variation as well as interesting evolutionary trajectories of NBS-encoding genes in the different Brassica subgenomes, while co-localization of NBS-encoding genes and resistance QTL may facilitate resistance breeding in oilseed rape.

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

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

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

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

Genome-Wide Identification of NBS-Encoding Resistance Genes in Sunflower (Helianthus annuus L.)

Nucleotide Binding Site—Leucine-Rich Repeat (NBS-LRR) genes encode disease resistance proteins involved in plants’ defense against their pathogens. Although sunflower is affected by many diseases, only a few molecular details have been uncovered regarding pathogenesis and resistance mechanisms. Recent availability of sunflower whole genome sequences in publicly accessible databases allowed us to accomplish a genome-wide identification of Toll-interleukin-1 receptor-like Nucleotide-binding site Leucine-rich repeat (TNL), Coiled Coil (CC)-NBS-LRR (CNL), Resistance to powdery mildew8 (RPW8)-NBS-LRR (RNL) and NBS-LRR (NL) protein encoding genes. Hidden Markov Model (HMM) profiling of 52,243 putative protein sequences from sunflower resulted in 352 NBS-encoding genes, among which 100 genes belong to CNL group including 64 genes with RX_CC like domain, 77 to TNL, 13 to RNL, and 162 belong to NL group. We also identified signal peptides and nuclear localization signals present in the identified genes and their homologs. We found that NBS genes were located on all chromosomes and formed 75 gene clusters, one-third of which were located on chromosome 13. Phylogenetic analyses between sunflower and Arabidopsis NBS genes revealed a clade-specific nesting pattern in CNLs, with RNLs nested in the CNL-A clade, and species-specific nesting pattern for TNLs. Surprisingly, we found a moderate bootstrap support (BS = 50%) for CNL-A clade being nested within TNL clade making both the CNL and TNL clades paraphyletic. Arabidopsis and sunflower showed 87 syntenic blocks with 1049 high synteny hits between chromosome 5 of Arabidopsis and chromosome 6 of sunflower. Expression data revealed functional divergence of the NBS genes with basal level tissue-specific expression. This study represents the first genome-wide identification of NBS genes in sunflower paving avenues for functional characterization and potential crop improvement.

Oak genome reveals facets of long lifespan

Oaks are an important part of our natural and cultural heritage. Not only are they ubiquitous in our most common landscapes¹ but they have also supplied human societies with invaluable services, including food and shelter, since prehistoric times². With 450 species spread throughout Asia, Europe and America³, oaks constitute a critical global renewable resource. The longevity of oaks (several hundred years) probably underlies their emblematic cultural and historical importance. Such long-lived sessile organisms must persist in the face of a wide range of abiotic and biotic threats over their lifespans. We investigated the genomic features associated with such a long lifespan by sequencing, assembling and annotating the oak genome. We then used the growing number of whole-genome sequences for plants (including tree and herbaceous species) to investigate the parallel evolution of genomic characteristics potentially underpinning tree longevity. A further consequence of the long lifespan of trees is their accumulation of somatic mutations during mitotic divisions of stem cells present in the shoot apical meristems. Empirical⁴ and modelling⁵ approaches have shown that intra-organismal genetic heterogeneity can be selected for⁶ and provides direct fitness benefits in the arms race with short-lived pests and pathogens through a patchwork of intra-organismal phenotypes⁷. However, there is no clear proof that large-statured trees consist of a genetic mosaic of clonally distinct cell lineages within and between branches. Through this case study of oak, we demonstrate the accumulation and transmission of somatic mutations and the expansion of disease-resistance gene families in trees.

Identification of NBS-encoding genes linked to black rot resistance in cabbage (Brassica oleracea var. capitata)

Heading cabbage is a nutritionally rich and economically important cruciferous vegetable. Black rot disease, caused by the bacterium Xanthomonas campestris pv. campestris, reduces both the yield and quality of the cabbage head. Nucleotide binding site (NBS)-encoding resistance (R) genes play a vital role in the plant immune response to various pathogens. In this study, we analyzed the expression and DNA sequence variation of 31 NBS-encoding genes in cabbage (Brassica oleracea var. capitata). These genes encoded TIR, NBS, LRR and RPW8 protein domains, all of which are known to be involved in disease resistance. RNA-seq revealed that these 31 genes were differentially expressed in leaf, root, silique, and stem tissues. Furthermore, qPCR analyses revealed that several of these genes were more highly expressed in resistant compared to susceptible cabbage lines, including Bol003711, Bol010135, Bol010559, Bol022784, Bol029866, Bol042121, Bol031422, Bol040045 and Bol042095. Further analysis of these genes promises to yield both practical benefits, such as molecular markers for marker-assisted breeding, and fundamental insights to the mechanisms of resistance to black rot in cabbage.

Genome-wide identification and expression analysis of chitinase gene family in Brassica rapa reveals its role in clubroot resistance

Chitinases, a category of pathogenesis-related proteins, are responsible for catalyzing the hydrolysis of chitin into the N-acetyl-d-glucosamine. Therefore, chitinases are believed to function as a guardian against chitin-containing pathogens. Here, we examined the role of the Brassica rapa chitinase family genes in clubroot disease. A total of 33 chitinase genes were identified and grouped into five classes based on their conserved domain. They were distributed unevenly across eight chromosomes in B. rapa, and 31 of them contained few introns (≤2). In addition, the expression of these genes was organ-specific, and 14 genes were expressed differentially in response to Plasmodiophora brassicae challenge of clubroot-susceptible (CS NIL) and resistant (CR NIL) lines. Furthermore, reduced pathogen DNA content and clubroot symptoms were observed in the CS NILs after their treatment with chitin oligosaccharides 24 h prior to inoculation with P. brassicae. The findings indicate that chitinases play a crucial role in pathogen resistance of the host plants. The results offer an insight into the role of chitinase in B. rapa-P. brassicae interaction.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Abiotic Stress Responsive miRNA-Target Network and Related Markers (SNP, SSR) in Brassica juncea

Abiotic stress is one of the major factors responsible for huge yield loss in crop plants. MicroRNAs play a key role in adaptive responses of plants under abiotic stress conditions through post-transcriptional gene regulations. In present study, 95 potential miRNAs were predicted in Brassica juncea using comparative genomics approach. It was noted that these miRNAs, target several transcription factors (TFs), transporter family proteins, signaling related genes, and protease encoding genes. Nineteen distinct miRNA-target regulatory networks were observed with significant involvement in regulation of transcription, response to stimulus, hormone and auxin mediated signaling pathway related gene ontology (GO) term. The sucrose-starch metabolism and pentose-gluconate interconversion pathways were found significantly enriched for these target genes. Molecular markers such as Simple Sequence Repeats (SSR) and Single Nucleotide Polymorphism (SNPs) were identified on miRNAs (miR-SSRs and miR-SNPs) and their target genes in B. juncea. Notably, one of the miR-SNP (C/T) was found at the 5^th position on mature region of miR2926. This C/T transition led to the distorted and unstable hairpin structure of miR2926, consequently complete loss of target function. Hence, findings from this study will lay a foundation for marker assisted breeding for abiotic stress tolerant varieties of B. juncea.

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.

Genome-wide identification and tissue-specific expression analysis of nucleotide binding site-leucine rich repeat gene family in Cicer arietinum (kabuli chickpea)

The nucleotide binding site-leucine rich repeat (NBS-LRR) proteins play an important role in the defense mechanisms against pathogens. Using bioinformatics approach, we identified and annotated 104 NBS-LRR genes in chickpea. Phylogenetic analysis points to their diversification into two families namely TIR-NBS-LRR and non-TIR-NBS-LRR. Gene architecture revealed intron gain/loss events in this resistance gene family during their independent evolution into two families. Comparative genomics analysis elucidated its evolutionary relationship with other fabaceae species. Around 50% NBS-LRRs reside in macro-syntenic blocks underlining positional conservation along with sequence conservation of NBS-LRR genes in chickpea. Transcriptome sequencing data provided evidence for their transcription and tissue-specific expression. Four cis-regulatory elements namely WBOX, DRE, CBF, and GCC boxes, that commonly occur in resistance genes, were present in the promoter regions of these genes. Further, the findings will provide a strong background to use candidate disease resistance NBS-encoding genes and identify their specific roles in chickpea.

Genome-Wide Association Study Identifies NBS-LRR-Encoding Genes Related with Anthracnose and Common Bacterial Blight in the Common Bean

Nucleotide-binding site and leucine-rich repeat (NBS-LRR) genes represent the largest and most important disease resistance genes in plants. The genome sequence of the common bean (Phaseolus vulgaris L.) provides valuable data for determining the genomic organization of NBS-LRR genes. However, data on the NBS-LRR genes in the common bean are limited. In total, 178 NBS-LRR-type genes and 145 partial genes (with or without a NBS) located on 11 common bean chromosomes were identified from genome sequences database. Furthermore, 30 NBS-LRR genes were classified into Toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) types, and 148 NBS-LRR genes were classified into coiled-coil (CC)-NBS-LRR (CNL) types. Moreover, the phylogenetic tree supported the division of these PvNBS genes into two obvious groups, TNL types and CNL types. We also built expression profiles of NBS genes in response to anthracnose and common bacterial blight using qRT-PCR. Finally, we detected nine disease resistance loci for anthracnose (ANT) and seven for common bacterial blight (CBB) using the developed NBS-SSR markers. Among these loci, NSSR24, NSSR73, and NSSR265 may be located at new regions for ANT resistance, while NSSR65 and NSSR260 may be located at new regions for CBB resistance. Furthermore, we validated NSSR24, NSSR65, NSSR73, NSSR260, and NSSR265 using a new natural population. Our results provide useful information regarding the function of the NBS-LRR proteins and will accelerate the functional genomics and evolutionary studies of NBS-LRR genes in food legumes. NBS-SSR markers represent a wide-reaching resource for molecular breeding in the common bean and other food legumes. Collectively, our results should be of broad interest to bean scientists and breeders.

Functional differentiation and spatial-temporal co-expression networks of the NBS-encoding gene family in Jilin ginseng, Panax ginseng C.A. Meyer

Ginseng, Panax ginseng C.A. Meyer, is one of the most important medicinal plants for human health and medicine. It has been documented that over 80% of genes conferring resistance to bacteria, viruses, fungi and nematodes are contributed by the nucleotide binding site (NBS)-encoding gene family. Therefore, identification and characterization of NBS genes expressed in ginseng are paramount to its genetic improvement and breeding. However, little is known about the NBS-encoding genes in ginseng. Here we report genome-wide identification and systems analysis of the NBS genes actively expressed in ginseng (PgNBS genes). Four hundred twelve PgNBS gene transcripts, derived from 284 gene models, were identified from the transcriptomes of 14 ginseng tissues. These genes were classified into eight types, including TNL, TN, CNL, CN, NL, N, RPW8-NL and RPW8-N. Seven conserved motifs were identified in both the Toll/interleukine-1 receptor (TIR) and coiled-coil (CC) typed genes whereas six were identified in the RPW8 typed genes. Phylogenetic analysis showed that the PgNBS gene family is an ancient family, with a vast majority of its genes originated before ginseng originated. In spite of their belonging to a family, the PgNBS genes have functionally dramatically differentiated and been categorized into numerous functional categories. The expressions of the across tissues, different aged roots and the roots of different genotypes. However, they are coordinating in expression, forming a single co-expression network. These results provide a deeper understanding of the origin, evolution and functional differentiation and expression dynamics of the NBS-encoding gene family in plants in general and in ginseng particularly, and a NBS gene toolkit useful for isolation and characterization of disease resistance genes and for enhanced disease resistance breeding in ginseng and related species.

Genome reconstruction in Cynara cardunculus taxa gains access to chromosome-scale DNA variation

The genome sequence of globe artichoke (Cynara cardunculus L. var. scolymus, 2n = 2x = 34) is now available for use. A survey of C. cardunculus genetic resources is essential for understanding the evolution of the species, carrying out genetic studies and for application of breeding strategies. We report on the resequencing analyses (~35×) of four globe artichoke genotypes, representative of the core varietal types, as well as a genotype of the related taxa cultivated cardoon. The genomes were reconstructed at a chromosomal scale and structurally/functionally annotated. Gene prediction indicated a similar number of genes, while distinctive variations in miRNAs and resistance gene analogues (RGAs) were detected. Overall, 23,5 M SNP/indel were discovered (range 6,34 M –14,50 M). The impact of some missense SNPs on the biological functions of genes involved in the biosynthesis of phenylpropanoid and sesquiterpene lactone secondary metabolites was predicted. The identified variants contribute to infer on globe artichoke domestication of the different varietal types, and represent key tools for dissecting the path from sequence variation to phenotype. The new genomic sequences are fully searchable through independent Jbrowse interfaces (www.artichokegenome.unito.it), which allow the analysis of collinearity and the discovery of genomic variants, thus representing a one-stop resource for C. cardunculus genomics.

Genome-wide identification and characterization of NB-ARC resistant genes in wheat (Triticum aestivum L.) and their expression during leaf rust infection

NB-ARC domain-containing resistance genes from the wheat genome were identified, characterized and localized on chromosome arms that displayed differential yet positive response during incompatible and compatible leaf rust interactions. Wheat (Triticum aestivum L.) is an important cereal crop; however, its production is affected severely by numerous diseases including rusts. An efficient, cost-effective and ecologically viable approach to control pathogens is through host resistance. In wheat, high numbers of resistance loci are present but only few have been identified and cloned. A comprehensive analysis of the NB-ARC-containing genes in complete wheat genome was accomplished in this study. Complete NB-ARC encoding genes were mined from the Ensembl Plants database to predict 604 NB-ARC containing sequences using the HMM approach. Genome-wide analysis of orthologous clusters in the NB-ARC-containing sequences of wheat and other members of the Poaceae family revealed maximum homology with Oryza sativa indica and Brachypodium distachyon. The identification of overlap between orthologous clusters enabled the elucidation of the function and evolution of resistance proteins. The distributions of the NB-ARC domain-containing sequences were found to be balanced among the three wheat sub-genomes. Wheat chromosome arms 4AL and 7BL had the most NB-ARC domain-containing contigs. The spatio-temporal expression profiling studies exemplified the positive role of these genes in resistant and susceptible wheat plants during incompatible and compatible interaction in response to the leaf rust pathogen Puccinia triticina. Two NB-ARC domain-containing sequences were modelled in silico, cloned and sequenced to analyze their fine structures. The data obtained in this study will augment isolation, characterization and application NB-ARC resistance genes in marker-assisted selection based breeding programs for improving rust resistance in wheat.

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

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

Genome-wide comparative analysis of NBS-encoding genes in four Gossypium species

Background

Nucleotide binding site (NBS) genes encode a large family of disease resistance (R) proteins in plants. The availability of genomic data of the two diploid cotton species, Gossypium arboreum and Gossypium raimondii, and the two allotetraploid cotton species, Gossypium hirsutum (TM-1) and Gossypium barbadense allow for a more comprehensive and systematic comparative study of NBS-encoding genes to elucidate the mechanisms of cotton disease resistance.

Results

Based on the genome assembly data, 246, 365, 588 and 682 NBS-encoding genes were identified in G. arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. The distribution of NBS-encoding genes among the chromosomes was nonrandom and uneven, and was tended to form clusters. Gene structure analysis showed that G. arboreum and G. hirsutum possessed a greater proportion of CN, CNL, and N genes and a lower proportion of NL, TN and TNL genes compared to that of G. raimondii and G. barbadense, while the percentages of RN and RNL genes remained relatively unchanged. The percentage changes among them were largest for TNL genes, about 7 times. Exon statistics showed that the average exon numbers per NBS gene in G. raimondii and G. barbadense were all greater than that in G. arboretum and G. hirsutum. Phylogenetic analysis revealed that the TIR-NBS genes of G. barbadense were closely related with that of G. raimondii. Sequence similarity analysis showed that diploid cotton G. arboreum possessed a larger proportion of NBS-encoding genes similar to that of allotetraploid cotton G. hirsutum, while diploid G. raimondii possessed a larger proportion of NBS-encoding genes similar to that of allotetraploid cotton G. barbadense. The synteny analysis showed that more NBS genes in G. raimondii and G. arboreum were syntenic with that in G. barbadense and G. hirsutum, respectively.

Conclusions

The structural architectures, amino acid sequence similarities and synteny of NBS-encoding genes between G. arboreum and G. hirsutum, and between G. raimondii and G. barbadense were the highest among comparisons between the diploid and allotetraploid genomes, indicating that G. hirsutum inherited more NBS-encoding genes from G. arboreum, while G. barbadense inherited more NBS-encoding genes from G. raimondii. This asymmetric evolution of NBS-encoding genes may help to explain why G. raimondii and G. barbadense are more resistant to Verticillium wilt, whereas G. arboreum and G. hirsutum are more susceptible to Verticillium wilt. The disease resistances of the allotetraploid cotton were related to their NBS-encoding genes especially in regard from which diploid progenitor they were derived, and the TNL genes may have a significant role in disease resistance to Verticillium wilt in G. raimondii and G. barbadense.

Electronic supplementary material

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

Comparative analysis of NBS-LRR genes and their response to Aspergillus flavus in Arachis

Studies have demonstrated that nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes respond to pathogen attack in plants. Characterization of NBS-LRR genes in peanut is not well documented. The newly released whole genome sequences of Arachis duranensis and Arachis ipaënsis have allowed a global analysis of this important gene family in peanut to be conducted. In this study, we identified 393 (AdNBS) and 437 (AiNBS) NBS-LRR genes from A. duranensis and A. ipaënsis, respectively, using bioinformatics approaches. Full-length sequences of 278 AdNBS and 303 AiNBS were identified. Fifty-one orthologous, four AdNBS paralogous, and six AiNBS paralogous gene pairs were predicted. All paralogous gene pairs were located in the same chromosomes, indicating that tandem duplication was the most likely mechanism forming these paralogs. The paralogs mainly underwent purifying selection, but most LRR 8 domains underwent positive selection. More gene clusters were found in A. ipaënsis than in A. duranensis, possibly owing to tandem duplication events occurring more frequently in A. ipaënsis. The expression profile of NBS-LRR genes was different between A. duranensis and A. hypogaea after Aspergillus flavus infection. The up-regulated expression of NBS-LRR in A. duranensis was continuous, while these genes responded to the pathogen temporally in A. hypogaea.

Genetic Analysis of NBS-LRR Gene Family in Chickpea and Their Expression Profiles in Response to Ascochyta Blight Infection

Ascochyta blight is one of the major diseases of chickpea worldwide. The genetic resistance to ascochyta blight in chickpea is complex and governed by multiple QTLs. However, the molecular mechanism of quantitative disease resistance to ascochyta blight and the genes underlying these QTLs are still unknown. Most often disease resistance is determined by resistance (R) genes. The most predominant R-genes contain nucleotide binding site and leucine rich repeat (NBS-LRR) domains. A total of 121 NBS-LRR genes were identified in the chickpea genome. Ninety-eight of these genes contained all essential conserved domains while 23 genes were truncated. The NBS-LRR genes were grouped into eight distinct classes based on their domain architecture. Phylogenetic analysis grouped these genes into two major clusters based on their structural variation, the first cluster with toll or interleukin-1 like receptor (TIR) domain and the second cluster either with or without a coiled-coil domain. The NBS-LRR genes are distributed unevenly across the eight chickpea chromosomes and nearly 50% of the genes are present in clusters. Thirty of the NBS-LRR genes were co-localized with nine of the previously reported ascochyta blight QTLs and were tested as potential candidate genes for ascochyta blight resistance. Expression pattern of these genes was studied in two resistant (CDC Corinne and CDC Luna) and one susceptible (ICCV 96029) genotypes at different time points after ascochyta blight infection using real-time quantitative PCR. Twenty-seven NBS-LRR genes showed differential expression in response to ascochyta blight infection in at least one genotype at one time point. Among these 27 genes, the majority of the NBS-LRR genes showed differential expression after inoculation in both resistant and susceptible genotypes which indicates the involvement of these genes in response to ascochyta blight infection. Five NBS-LRR genes showed genotype specific expression. Our study provides a new insight of NBS-LRR gene family in chickpea and the potential involvement of NBS-LRR genes in response to ascochyta blight infection.

Identification, characterization, and gene expression analysis of nucleotide binding site (NB)-type resistance gene homologues in switchgrass

BACKGROUND

Switchgrass (Panicum virgatum L.) is a warm-season perennial grass that can be used as a second generation bioenergy crop. However, foliar fungal pathogens, like switchgrass rust, have the potential to significantly reduce switchgrass biomass yield. Despite its importance as a prominent bioenergy crop, a genome-wide comprehensive analysis of NB-LRR disease resistance genes has yet to be performed in switchgrass.

RESULTS

In this study, we used a homology-based computational approach to identify 1011 potential NB-LRR resistance gene homologs (RGHs) in the switchgrass genome (v 1.1). In addition, we identified 40 RGHs that potentially contain unique domains including major sperm protein domain, jacalin-like binding domain, calmodulin-like binding, and thioredoxin. RNA-sequencing analysis of leaf tissue from 'Alamo', a rust-resistant switchgrass cultivar, and 'Dacotah', a rust-susceptible switchgrass cultivar, identified 2634 high quality variants in the RGHs between the two cultivars. RNA-sequencing data from field-grown cultivar 'Summer' plants indicated that the expression of some of these RGHs was developmentally regulated.

CONCLUSIONS

Our results provide useful insight into the molecular structure, distribution, and expression patterns of members of the NB-LRR gene family in switchgrass. These results also provide a foundation for future work aimed at elucidating the molecular mechanisms underlying disease resistance in this important bioenergy crop.

Oilseed rape: learning about ancient and recent polyploid evolution from a recent crop species

Oilseed rape (Brassica napus) is one of our youngest crop species, arising several times under cultivation in the last few thousand years and completely unknown in the wild. Oilseed rape originated from hybridisation events between progenitor diploid species B. rapa and B. oleracea, both important vegetable species. The diploid progenitors are also ancient polyploids, with remnants of two previous polyploidisation events evident in the triplicated genome structure. This history of polyploid evolution and human agricultural selection makes B. napus an excellent model with which to investigate processes of genomic evolution and selection in polyploid crops. The ease of de novo interspecific hybridisation, responsiveness to tissue culture, and the close relationship of oilseed rape to the model plant Arabidopsis thaliana, coupled with the recent availability of reference genome sequences and suites of molecular cytogenetic and high-throughput genotyping tools, allow detailed dissection of genetic, genomic and phenotypic interactions in this crop. In this review we discuss the past and present uses of B. napus as a model for polyploid speciation and evolution in crop species, along with current and developing analysis tools and resources. We further outline unanswered questions that may now be tractable to investigation.