Cloning, characterization, and evolution of the NBS-LRR-encoding resistance gene analogue family in polyploid cotton (Gossypium hirsutum L.)

Expression of an Antiviral Gene GmRUN1 from Soybean Is Regulated via Intron-Mediated Enhancement (IME)

Most of R (resistance) genes encode the protein containing NBS-LRR (nucleotide binding site and leucine-rich repeat) domains. Here, N. benthamiana plants were used for transient expression assays at 3-4 weeks of age. We identified a TNL (TIR-NBS-LRR) encoding gene GmRUN1 that was resistant to both soybean mosaic virus (SMV) and tobacco mosaic virus (TMV). Truncation analysis indicated the importance of all three canonical domains for GmRUN1-mediated antiviral activity. Promoter-GUS analysis showed that GmRUN1 expression is inducible by both salicylic acid (SA) and a transcription factor GmDREB3 via the cis-elements as-1 and ERE (ethylene response element), which are present in its promoter region. Interestingly, GmRUN1 gDNA (genomic DNA) shows higher viral resistance than its cDNA (complementary DNA), indicating the existence of intron-mediated enhancement (IME) for GmRUN1 regulation. We provided evidence that intron2 of GmRUN1 increased the mRNA level of native gene GmRUN1, a soybean antiviral gene SRC7 and also a reporter gene Luciferase, indicating the general transcriptional enhancement of intron2 in different genes. In summary, we identified an antiviral TNL type soybean gene GmRUN1, expression of which was regulated at different layers. The investigation of GmRUN1 gene regulatory network would help to explore the mechanism underlying soybean-SMV interactions.

Isolation and characterization of fusarium wilt resistance gene analogs in radish

The resistance gene analog (RGA)-based marker strategy is an effective supplement for current marker reservoir of radish disease-resistance breeding. In this study, we identified RGAs based on the conserved nucleotide-binding site (NBS) and S-receptor-like kinase (SRLK) domains. A total of 68 NBS-RGAs and 46 SRLK-RGAs were isolated from two FW-resistant radish inbred lines, B2 and YR31, and one susceptible line, YR15. A BLASTx search revealed that the NBS-RGAs contained six conserved motifs (i.e., P loop, RNBS-A, Kinase-2, RNBS-B, RNBS-C, and GLPL) and the SRLK-RGAs, contained two conserved motifs (i.e., G-type lectin and PAN-AP). A phylogenetic analysis indicated that the NBS-RGAs could be separated into two classes (i.e., toll/interleukin receptor and coiled-coil types), with six subgroups, and the SRLK-RGAs were divided into three subgroups. Moreover, we designed RGA-specific markers from data-mining approach in radish databases. Based on marker analysis, 24 radish inbred lines were clustered into five main groups with a similarity index of 0.44 and showing genetic diversity with resistance variation in those radish inbred lines. The development of RGA-specific primers would be valuable for marker-assisted selection during the breeding of disease-resistant radish cultivars.

Genome-wide association study discovered candidate genes of Verticillium wilt resistance in upland cotton (Gossypium hirsutum L.)

Verticillium wilt (VW), caused by infection by Verticillium dahliae, is considered one of the most yield-limiting diseases in cotton. To examine the genetic architecture of cotton VW resistance, we performed a genome-wide association study (GWAS) using a panel of 299 accessions and 85 630 single nucleotide polymorphisms (SNPs) detected using the specific-locus amplified fragment sequencing (SLAF-seq) approach. Trait-SNP association analysis detected a total of 17 significant SNPs at P < 1.17 × 10-5 (P = 1/85 630, -log10 P = 4.93); the peaks of SNPs associated with VW resistance on A10 were continuous and common in three environments (RDIG2015, RDIF2015 and RDIF2016). Haplotype block structure analysis predicted 22 candidate genes for VW resistance based on A10_99672586 with a minimum P-value (-log10 P = 6.21). One of these genes (CG02) was near the significant SNP A10_99672586 (0.26 Mb), located in a 372-kb haplotype block, and its Arabidopsis AT3G25510 homologues contain TIR-NBS-LRR domains that may be involved in disease resistance response. Real-time quantitative PCR and virus-induced gene silencing (VIGS) analysis showed that CG02 was specific to up-regulation in the resistant (R) genotype Zhongzhimian2 (ZZM2) and that silenced plants were more susceptible to V. dahliae. These results indicate that CG02 is likely the candidate gene for resistance against V. dahliae in cotton. The identified locus or gene may serve as a promising target for genetic engineering and selection for improving resistance to VW in cotton.

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.

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.

Dramatic Number Variation of R Genes in Solanaceae Species Accounted for by a Few R Gene Subfamilies

Most disease resistance genes encode nucleotide-binding-site (NBS) and leucine-rich-repeat (LRR) domains, and the NBS-LRR encoding genes are often referred to as R genes. Using newly developed approach, 478, 485, 1,194, 1,665, 2,042 and 374 R genes were identified from the genomes of tomato Heinz1706, wild tomato LA716, potato DM1-3, pepper Zunla-1 and wild pepper Chiltepin and tobacco TN90, respectively. The majority of R genes from Solanaceae were grouped into 87 subfamilies, including 16 TIR-NBS-LRR (TNL) and 71 non-TNL subfamilies. Each subfamily was annotated manually, including identification of intron/exon structure and intron phase. Interestingly, TNL subfamilies have similar intron phase patterns, while the non-TNL subfamilies have diverse intron phase due to frequent gain of introns. Prevalent presence/absence polymorphic R gene loci were found among Solanaceae species, and an integrated map with 427 R loci was constructed. The pepper genome (2,042 in Chiltepin) has at least four times of R genes as in tomato (478 in Heinz1706). The high number of R genes in pepper genome is due to the amplification of R genes in a few subfamilies, such as the Rpi-blb2 and BS2 subfamilies. The mechanism underlying the variation of R gene number among different plant genomes is discussed.

Comparative analysis of resistance gene analogues encoding NBS-LRR domains in cotton

BACKGROUND

Plant production is severely affected by biotic and abiotic stresses R-genes exhibit resistance against a range of diseases and pathogens in plants. The nucleotide binding site and leucine rich repeat (NBS-LRR) class of R-genes is the most comprehensively studied in terms of sequence evolution and genome distribution. The differential response for resistance against biotic and abiotic stress has been observed in cultivated and wild relatives of the genus Gossypium.

RESULTS

Efforts have been made to address the recent evolution of NBS-LRR sequences within Gossypium hirsutum and resistance gene analogue (RGA) sequences derived from G. arboreum and G. raimondii. The % identity and phylogenetic analysis of NBS-LRR-encoded RGAs from tetraploid New World cotton and its diploid ancestors G. raimondii and G. arboreum suggest that the evolution of NBS-LRR-encoding sequences in G. hirsutum occurred by gradual accumulation of mutants that led to positive selection and a slow rate of divergence within distinct R-gene families.

CONCLUSION

The allotetraploid genome of cotton, after separating from its diploid parents, experienced polyploidisation, natural and artificial selection, hybrid necrosis, duplication and recombination which became the reason to shed off and evolve new genes for its survival. These driving forces influenced the development of genomic architecture that make it susceptible to diseases and pathogens as compared to donor parents.

Mining, genetic mapping and expression analysis of EST-derived resistance gene homologs (RGHs) in cotton

BACKGROUND

Cotton is the dominant textile crop and also serves as an important oil crop. An estimated 15% economic loss associated with cotton production in China has been caused by diseases, and no resistance genes have been cloned in this crop. Molecular markers developed from resistance gene homologues (RGHs) might be tightly linked with target genes and could be used for marker-assisted selection (MAS) or gene cloning.

RESULTS

To genetically map expressed RGHs, 100 potential pathogenesis-related proteins (PRPs) and 215 resistance gene analogs (RGAs) were identified in the cotton expressed sequence tag database, and 347 specific primers were developed. Meanwhile, 61 cotton genome-derived RGA markers and 24 resistance gene analog polymorphism (RGAP) markers from published papers were included to view their genomic distribution. As a result, 38 EST-derived and 17 genome-derived RGH markers were added to our interspecific genetic map. These 55 markers were distributed on 18 of the 26 cotton chromosomes, with 34 markers on 6 chromosomes (Chr03, Chr04, Chr11, Chr17, Chr19 and Chr26). Homologous RGHs tended to be clustered; RGH clusters appeared on 9 chromosomes, with larger clusters on Chr03, Chr04 and Chr19, which suggests that RGH clusters are widely distributed in the cotton genome. Expression analysis showed that 19 RGHs were significantly altered after inoculation with the V991 stain of Verticillium dahliae. Comparative mapping showed that four RGH markers were linked with mapped loci for Verticillium wilt resistance.

CONCLUSIONS

The genetic mapping of RGHs confirmed their clustering in cotton genome. Expression analysis and comparative mapping suggest that EST-derived RGHs participate in cotton resistance. RGH markers are seemed to be useful tools to detected resistance loci and identify candidate resistance genes in cotton.

In-silico mining, type and frequency analysis of genic microsatellites of finger millet (Eleusine coracana (L.) Gaertn.): a comparative genomic analysis of NBS-LRR regions of finger millet with rice

In recent years, the increased availability of the DNA sequences has given the possibility to develop and explore the expressed sequence tags (ESTs) derived SSR markers. In the present study, a total of 1956 ESTs of finger millet were used to find the microsatellite type, distribution, frequency and developed a total of 545 primer pairs from the ESTs of finger millet. Thirty-two EST sequences had more than two microsatellites and 1357 sequences did not have any SSR repeats. The most frequent type of repeats was trimeric motif, however the second place was occupied by dimeric motif followed by tetra-, hexa- and penta repeat motifs. The most common dimer repeat motif was GA and in case of trimeric SSRs, it was CGG. The EST sequences of NBS-LRR region of finger millet and rice showed higher synteny and were found on nearly same positions on the rice chromosome map. A total of eight, out of 15 EST based SSR primers were polymorphic among the selected resistant and susceptible finger millet genotypes. The primer FMBLEST5 could able to differentiate them into resistant and susceptible genotypes. The alleles specific to the resistant and susceptible genotypes were sequenced using the ABI 3130XL genetic analyzer and found similarity to NBS-LRR regions of rice and finger millet and contained the characteristic kinase-2 and kinase 3a motifs of plant R-genes belonged to NBS-LRR region. The In-silico and comparative analysis showed that the genes responsible for blast resistance can be identified, mapped and further introgressed through molecular breeding approaches for enhancing the blast resistance in finger millet.

miR482 regulation of NBS-LRR defense genes during fungal pathogen infection in cotton

In this study, we characterized the miR482 family in cotton using existing small RNA datasets and the recently released draft genome sequence of Gossypium raimondii, a diploid cotton species whose progenitor is the putative contributor of the Dt (representing the D genome of tetraploid) genome of the cultivated tetraploid cotton species G. hirsutum and G. barbadense. Of the three ghr-miR482 members reported in G. hirsutum, ghr-miR482a has no homolog in G. raimondii, ghr-miR482b and ghr-miR482c each has a single homolog in G. raimondii. Gra-miR482d has five homologous loci (gra-miR482d, f-i) in G. raimondii and also exists in G. hirsutum (ghr-miR482d). A variant, miR482.2 that is a homolog of miR2118 in other species, is produced from several GHR-MIR482 loci in G. hirsutum. Approximately 12% of the G. raimondii NBS-LRR genes were predicted targets of various members of the gra-miR482 family. Based on the rationale that the regulatory relationship between miR482 and NBS-LRR genes will be conserved in G. raimondii and G. hirsutum, we investigated this relationship using G. hirsutum miR482 and G. raimondii NBS-LRR genes, which are not currently available in G. hirsutum. Ghr-miR482/miR482.2-mediated cleavage was confirmed for three of the four NBS-LRR genes analysed. As in tomato, miR482-mediated cleavage of NBS-LRR genes triggered production of phased secondary small RNAs in cotton. In seedlings of the susceptible cultivar Sicot71 (G. hirsutum) infected with the fungal pathogen Verticillium dahliae, the expression levels of ghr-miR482b/miR482b.2, ghr-miR482c and ghr-miR482d.2 were down-regulated, and several NBS-LRR targets of ghr-miR482c and ghr-miR482d were up-regulated. These results imply that, like tomato plants infected with viruses or bacteria, cotton plants are able to induce expression of NBS-LRR defence genes by suppression of the miRNA-mediated gene silencing pathway upon fungal pathogen attack.

Plant nucleotide binding site-leucine-rich repeat (NBS-LRR) genes: active guardians in host defense responses

The most represented group of resistance genes are those of the nucleotide binding site-leucine-rich repeat (NBS-LRR) class. These genes are very numerous in the plant genome, and they often occur in clusters at specific loci following gene duplication and amplification events. To date, hundreds of resistance genes and relatively few quantitative trait loci for plant resistance to pathogens have been mapped in different species, with some also cloned. When these NBS-LRR genes have been physically or genetically mapped, many cases have shown co-localization between resistance loci and NBS-LRR genes. This has allowed the identification of candidate genes for resistance, and the development of molecular markers linked to R genes. This review is focused on recent genomics studies that have described the abundance, distribution and evolution of NBS-LRR genes in plant genomes. Furthermore, in terms of their expression, NBS-LRR genes are under fine regulation by cis- and trans-acting elements. Recent findings have provided insights into the roles of alternative splicing, the ubiquitin/ proteasome system, and miRNAs and secondary siRNAs in the regulation of NBS-LRR gene expression at the post-transcriptional, post-translational and epigenetic levels. The possibility to use this knowledge for genetic improvement of plant resistance to pathogens is discussed.

Construction of a plant-transformation-competent BIBAC library and genome sequence analysis of polyploid Upland cotton (Gossypium hirsutum L.)

Background

Cotton, one of the world’s leading crops, is important to the world’s textile and energy industries, and is a model species for studies of plant polyploidization, cellulose biosynthesis and cell wall biogenesis. Here, we report the construction of a plant-transformation-competent binary bacterial artificial chromosome (BIBAC) library and comparative genome sequence analysis of polyploid Upland cotton (Gossypium hirsutum L.) with one of its diploid putative progenitor species, G. raimondii Ulbr.

Results

We constructed the cotton BIBAC library in a vector competent for high-molecular-weight DNA transformation in different plant species through either Agrobacterium or particle bombardment. The library contains 76,800 clones with an average insert size of 135 kb, providing an approximate 99% probability of obtaining at least one positive clone from the library using a single-copy probe. The quality and utility of the library were verified by identifying BIBACs containing genes important for fiber development, fiber cellulose biosynthesis, seed fatty acid metabolism, cotton-nematode interaction, and bacterial blight resistance. In order to gain an insight into the Upland cotton genome and its relationship with G. raimondii, we sequenced nearly 10,000 BIBAC ends (BESs) randomly selected from the library, generating approximately one BES for every 250 kb along the Upland cotton genome. The retroelement Gypsy/DIRS1 family predominates in the Upland cotton genome, accounting for over 77% of all transposable elements. From the BESs, we identified 1,269 simple sequence repeats (SSRs), of which 1,006 were new, thus providing additional markers for cotton genome research. Surprisingly, comparative sequence analysis showed that Upland cotton is much more diverged from G. raimondii at the genomic sequence level than expected. There seems to be no significant difference between the relationships of the Upland cotton D- and A-subgenomes with the G. raimondii genome, even though G. raimondii contains a D genome (D₅).

Conclusions

The library represents the first BIBAC library in cotton and related species, thus providing tools useful for integrative physical mapping, large-scale genome sequencing and large-scale functional analysis of the Upland cotton genome. Comparative sequence analysis provides insights into the Upland cotton genome, and a possible mechanism underlying the divergence and evolution of polyploid Upland cotton from its diploid putative progenitor species, G. raimondii.

Evolutionary meta-analysis of solanaceous resistance gene and solanum resistance gene analog sequences and a practical framework for cross-species comparisons

Cross-species comparative genomics approaches have been employed to map and clone many important disease resistance (R) genes from Solanum species-especially wild relatives of potato and tomato. These efforts will increase with the recent release of potato genome sequence and the impending release of tomato genome sequence. Most R genes belong to the prominent nucleotide binding site-leucine rich repeat (NBS-LRR) class and conserved NBS-LRR protein motifs enable survey of the R gene space of a plant genome by generation of resistance gene analogs (RGA), polymerase chain reaction fragments derived from R genes. We generated a collection of 97 RGA from the disease-resistant wild potato S. bulbocastanum, complementing smaller collections from other Solanum species. To further comparative genomics approaches, we combined all known Solanum RGA and cloned solanaceous NBS-LRR gene sequences, nearly 800 sequences in total, into a single meta-analysis. We defined R gene diversity bins that reflect both evolutionary relationships and DNA cross-hybridization results. The resulting framework is amendable and expandable, providing the research community with a common vocabulary for present and future study of R gene lineages. Through a series of sequence and hybridization experiments, we demonstrate that all tested R gene lineages are of ancient origin, are shared between Solanum species, and can be successfully accessed via comparative genomics approaches.

Molecular evolution of the clustered MIC-3 multigene family of Gossypium species

The Gossypium MIC-3 (Meloidogyne Induced Cotton-3) gene family is of great interest for molecular evolutionary studies because of its uniqueness to Gossypium species, multi-gene content, clustered localization, and root-knot nematode resistance-associated features. Molecular evolution of the MIC-3 gene family was studied in 15 tetraploid and diploid Gossypium genotypes that collectively represent seven phylogenetically distinct genomes. Synonymous (d(S)) and non-synonymous (d(N)) nucleotide substitution rates suggest that the second of the two exons of the MIC-3 genes has been under strong positive selection pressure, while the first exon has been under strong purifying selection to preserve function. Based on nucleotide substitution rates, we conclude that MIC-3 genes are evolving by a birth-and-death process and that a 'gene amplification' mechanism has helped to retain all duplicate copies, which best fits with the "bait and switch" model of R-gene evolution. The data indicate MIC-3 gene duplication events occurred at various rates, once per 1 million years (MY) in the allotetraploids, once per ~2 MY in the A/F genome clade, and once per ~8 MY in the D-genome clade. Variations in the MIC-3 gene family seem to reflect evolutionary selection for increased functional stability, while also expanding the capacity to develop novel "switch" pockets for responding to diverse pests and pathogens. Such evolutionary roles are congruent with the hypothesis that members of this unique resistance gene family provide fitness advantages in Gossypium.

Molecular marker-assisted genotyping of mungbean yellow mosaic India virus resistant germplasms of mungbean and urdbean

Mungbean Yellow Mosaic India Virus (MYMIV) belonging to the genus begomovirus causes the yellow mosaic disease in a number of economically important edible grain legumes including mungbean (Vigna radiata), urdbean (Vigna mungo) and soybean (Glycine max). The disease is severe, critical, open spread and inflicts heavy yield losses annually. The objective of this study is to develop molecular markers linked to MYMIV-resistance to facilitate genotyping of urdbean and mungbean germplasms for MYMIV-reaction. Resistance-linked molecular markers were successfully developed from consensus motifs of other resistance (R) gene or R gene homologue sequences. Applying linked marker-assisted genotyping, plant breeders can carry out repeated genotyping throughout the growing season in absence of any disease incidence. Two MYMIV-resistance marker loci, YR4 and CYR1, were identified and of these two CYR1 is completely linked with MYMIV-resistant germplasms and co-segregating with MYMIV-resistant F₂, F₃ progenies of urdbean. The present study demonstrated that these two markers could be efficiently employed together in a multiplex-PCR-reaction for genotyping both V. mungo and V. radiata germplasms from field grown plants and also directly from the seed stock. This method of genotyping would save time and labour during the introgression of MYMIV-resistance through molecular breeding, as methods of phenotyping against begomoviruses are tedious, labour and time intensive.

Numbers of genes in the NBS and RLK families vary by more than four-fold within a plant species and are regulated by multiple factors

Many genes exist in the form of families; however, little is known about their size variation, evolution and biology. Here, we present the size variation and evolution of the nucleotide-binding site (NBS)-encoding gene family and receptor-like kinase (RLK) gene family in Oryza, Glycine and Gossypium. The sizes of both families vary by numeral fold, not only among species, surprisingly, also within a species. The size variations of the gene families are shown to correlate with each other, indicating their interactions, and driven by natural selection, artificial selection and genome size variation, but likely not by polyploidization. The numbers of genes in the families in a polyploid species are similar to those of one of its diploid donors, suggesting that polyploidization plays little roles in the expansion of the gene families and that organisms tend not to maintain their ‘surplus’ genes in the course of evolution. Furthermore, it is found that the size variations of both gene families are associated with organisms’ phylogeny, suggesting their roles in speciation and evolution. Since both selection and speciation act on organism’s morphological, physiological and biological variation, our results indicate that the variation of gene family size provides a source of genetic variation and evolution.

Phylogenetic analyses of peanut resistance gene candidates and screening of different genotypes for polymorphic markers

The nucleotide-binding-site-leucine-rich-repeat (NBS-LRR)-encoding gene family has attracted much research interest because approximately 75% of the plant disease resistance genes that have been cloned to date are from this gene family. Here, we describe a collection of peanut NBS-LRR resistance gene candidates (RGCs) isolated from peanut (Arachis) species by mining Gene Bank data base. NBS-LRR sequences assembled into TIR-NBS-LRR (75.4%) and non-TIR-NBS-LRR (24.6%) subfamilies. Total of 20 distinct clades were identified and showed a high level of sequence divergence within TIR-NBS and non-TIR-NBS subfamilies. Thirty-four primer pairs were designed from these RGC sequences and used for screening different genotypes belonging to wild and cultivated peanuts. Therefore, peanut RGC identified in this study will provide useful tools for developing DNA markers and cloning the genes for resistance to different pathogens in peanut.

Genetic diversity and genomic distribution of homologs encoding NBS-LRR disease resistance proteins in sunflower

Three-fourths of the recognition-dependent disease resistance genes (R-genes) identified in plants encode nucleotide binding site (NBS) leucine-rich repeat (LRR) proteins. NBS-LRR homologs have only been isolated on a limited scale from sunflower (Helianthus annuus L.), and most of the previously identified homologs are members of two large NBS-LRR clusters harboring downy mildew R-genes. We mined the sunflower EST database and used comparative genomics approaches to develop a deeper understanding of the diversity and distribution of NBS-LRR homologs in the sunflower genome. Collectively, 630 NBS-LRR homologs were identified, 88 by mining a database of 284,241 sunflower ESTs and 542 by sequencing 1,248 genomic DNA amplicons isolated from common and wild sunflower species. DNA markers were developed from 196 unique NBS-LRR sequences and facilitated genetic mapping of 167 NBS-LRR loci. The latter were distributed throughout the sunflower genome in 44 clusters or singletons. Wild species ESTs were a particularly rich source of novel NBS-LRR homologs, many of which were tightly linked to previously mapped downy mildew, rust, and broomrape R-genes. The DNA sequence and mapping resources described here should facilitate the discovery and isolation of recognition-dependent R-genes guarding sunflower from a broad spectrum of economically important diseases. Sunflower nucleotide and amino acid sequences have been deposited in DDBJ/EMBL/GenBank under accession numbers EF 560168-EF 559378 and ABQ 58077-ABQ 57529.

Characterization of resistance gene analogs with a nucleotide binding site isolated from a triploid white poplar

The majority of cloned plant disease resistance genes (R genes) encode a nucleotide binding site (NBS) and a leucine-rich repeat (LRR) domain. In this study, to better understand the R genes in white poplar, 59 resistance gene analogues (RGAs) were identified from a triploid white poplar [(Populus tomentosa x Populus bolleana) x P. tomentosa], based on conserved NBS regions. The 59 RGAs were phylogenetically classified into 10 subfamilies, and 54 RGAs with open-reading frames (ORFs) were further grouped into two classes, toll and interleukin-1 receptor (TIR) and non-TIR. BLAST searches with reference to the genomic sequence of Populus trichocarpa found 96 highly homologous regions distributed in 37 loci, suggesting the abundance and divergence of NBS-encoding genes in the triploid poplar genome. Within subfamilies 1-3, the average non-synonymous/synonymous substitution (omega) rates were < 1, indicating purifying selection on these RGAs, but some sites were clearly under diversifying selection with omega > 1. Many intergenic exchanges were also detected among these RGAs, indicating a probable role in homogenising NBS domains. Quantitative real-time PCR analysis revealed dramatic variations in the transcript level of 18 RGAs in the mature leaves, bark and roots of the triploid poplar, and identified two RGAs that had significantly higher level of transcripts in bark, four RGAs in mature leaves, and 14 in the above-ground portion of poplars, suggesting their probable roles in resistance against diseases attacking the organs. Our results shed light on genetic resources of poplar resistance and will be useful for further resistance gene isolation and exploitation.

Positive selection and expression divergence following gene duplication in the sunflower CYCLOIDEA gene family

Members of the CYCLOIDEA (CYC)/TEOSINTE-BRANCHED1 (TB1) group of transcription factors have been implicated in the evolution of zygomorphic (i.e., bilaterally symmetric) flowers in Antirrhinum and Lotus and the loss of branching phenotype during the domestication of maize. The composite inflorescences of sunflower (Helianthus annuus L. Asteraceae) contain both zygomorphic and actinomorphic (i.e., radially symmetric) florets (rays and disks, respectively), and the cultivated sunflower has evolved an unbranched phenotype in response to domestication from its highly branched wild progenitor; hence, genes related to CYC/TB1 are of great interest in this study system. We identified 10 members of the CYC/TB1 gene family in sunflower, which is more than found in any other species investigated to date. Phylogenetic analysis indicates that these genes occur in 3 distinct clades, consistent with previous research in other eudicot species. A combination of dating the duplication events and linkage mapping indicates that only some of the duplications were associated with polyploidization. Cosegregation between CYC-like genes and branching-related quantitative trait loci suggest a minor, if any, role for these genes in conferring differences in branching. However, the expression patterns of one gene suggest a possible role in the development of ray versus disk florets. Molecular evolutionary analyses reveal that residues in the conserved domains were the targets of positive selection following gene duplication. Taken together, these results indicate that gene duplication and functional divergence have played a major role in diversification of the sunflower CYC gene family.

The CC-NBS-LRR Subfamily in Pinus monticola: Targeted Identification, Gene Expression, and Genetic Linkage with Resistance to Cronartium ribicola

ABSTRACT To investigate disease resistance gene analogs (RGAs) encoding coiled-coil-nucelotide-binding site-leucine-rich repeats (CC-NBS-LRR) proteins in western white pine, degenerate primers targeting the conserved motifs in the NBS domain were designed to amplify RGAs from genomic DNA and cDNA. Sixty-one distinct RGAs were identified with identities to well-known R proteins of the CC-NBS-LRR subfamily. These RGAs exhibited variation of putative amino acid sequences from 13% to 98%, representing a complex CC-NBS-LRR subfamily. A phylogenetic tree constructed from the amino acid sequence alignment revealed that these 61 RGAs were grouped with other CC-NBS-LRR members from angiosperms, and could be further divided into six classes with an identity threshold of 68%. To map RGAs, RGA polymorphisms and a modified amplified fragment length polymorphism (AFLP) method with incorporated sequences from the NBS domain were used to reveal NBS or NBS-AFLP markers. RGA polymorphism study revealed that three off the identified RGAs were not linked to the Cr2 gene imparting resistance to white pine blister rust. However, the AFLP strategy, using bulk segregant analysis (BSA) and haploid segregation analysis, identified 11 NBS-AFLP markers localized in the Cr2 linkage, the closest two to the gene being 0.41 cM and 1.22 cM away at either side. Eight of these markers showed significant amino acid sequence homologies with RGAs.

Isolation and characterization of resistance and defense gene analogs in cotton (Gossypium barbadense L.)

Plant disease resistance gene (R gene) and defense response gene encode some conserved motifs. In the present work, a PCR strategy was used to clone resistance gene analogs (RGAs) and defense gene analogs (DGAs) from Sea-island cotton variety Hai7124 using oligonucleotide primers based on the nucleotide-binding site (NBS) and serine/threonine kinase (STK) in the R-gene and pathogenesis-related proteins of class 2 (PR2) of defense response gene. 79 NBS sequences, 21 STK sequences and 11 DGAs were cloned from disease-resistance cotton. Phylogenic analysis of 79 NBS-RGAs and NBS-RGAs nucleotide sequences of cotton already deposited in GenBank identified one new sub-cluster. The deduced amino acid sequences of NBS-RGAs and STK-RGAs were divided into two distinct groups respectively: Toll/Interleukin-1 receptor (TIR) group and non-TIR group, A group and B group. The expression of RGAs and DGAs having consecutive open reading frame (ORF) was also investigated and it was found that 6 NBS-RGAs and 1 STK-RGA were induced, and 1 DGA was up-regulated by infection of Verticillium dahliae strain VD8. 4 TIR-NBS-RGAs and 4 non-TIR-NBS-RGAs were arbitrarily used as probes for Southern-blotting. There existed 2-10 blotted bands. In addition, since three non-TIR-NBS-RGAs have the same hybridization pattern, we conjecture that these three RGAs form a cluster distribution in the genome.

Resistance gene analogues of wheat: molecular genetic analysis of ESTs

Using two divergent nucleotide binding site (NBS) regions from wheat sequences of the NBS-LRR (leucine rich repeat) class, we retrieved 211 wheat and barley NBS-containing resistance gene analogue (RGA) expressed sequence tags (ESTs). These ESTs were grouped into 129 gene sequence groups that contained ESTs that were at least 70% identical at the DNA level over at least 200 bp. Probes were obtained for 89 of these RGA families and chromosome locations were determined for 72 of these probes using nullitetrasomic Chinese Spring wheat lines. RFLP analysis of 49 of these RGA probes revealed 65 mappable polymorphic bands in the doubled haploid Cranbrook x Halberd wheat population (C x H). These bands mapped to 49 loci in C x H. RGA loci were detected on all 21 chromosomes using the nullitetrasomic lines and on 18 chromosomes (linkage groups) in the C x H map. This identified a set of potential markers that could be developed further for use in mapping and ultimately cloning NBS-LRR-type disease resistance genes in wheat.

Identification and mapping of microsatellite markers linked to a root-knot nematode resistance gene (rkn1) in Acala NemX cotton (Gossypium hirsutum L.)

Host-plant resistance is the most economic and effective strategy for root-knot nematode (RKN) Meloidogyne incognita control in cotton (Gossypium hirsutum L.). Molecular markers linked to resistance are important for incorporating resistance genes into elite cultivars. To screen for microsatellite markers (SSR) closely linked to RKN resistance in G. hirsutum cv. Acala NemX, F1, F2, BC1F1, and F2:7 recombinant inbred lines (RILs) from intraspecific crosses and an F2 from an interspecific cross with G. barbadense cv. Pima S-7 were used. Screening of 284 SSR markers, which cover all the known identified chromosomes and most linkage groups of cotton, was performed by bulked segregant analysis, revealing informative SSRs. The informative SSRs were then mapped on the above populations. One co-dominant SSR marker CIR316 was identified tightly linked to a major resistance gene (designated as rkn1), producing amplified DNA fragments of approximately 221 bp (CIR316a) and 210 bp (CIR316c) in Acala NemX and susceptible Acala SJ-2, respectively. The linkage between CIR316a marker and resistance gene rkn1 in Acala NemX had an estimated distance of 2.1-3.3 cM depending on the population used. Additional markers, including BNL1231 with loose linkage to rkn1 (map distance 25.1-27.4 cM), BNL1066, and CIR003 allowed the rkn1 gene to be mapped to cotton linkage group A03. This is the first report in cotton with a closely linked major gene locus determining nematode resistance, and informative SSRs may be used for marker-assisted selection.

Isolation of TIR and non-TIR NBS--LRR resistance gene analogues and identification of molecular markers linked to a powdery mildew resistance locus in chestnut rose (Rosa roxburghii Tratt)

Toll and interleukin-1 receptor (TIR) and non-TIR nucleotide binding site-leucine rich repeat (NBS-LRR) resistance gene analogues (RGAs) were obtained from chestnut rose (Rosa roxburghii Tratt) by two PCR-based amplification strategies (direct amplification and overlap extension amplification) with degenerate primers designed to the conserved P-loop, kinase-2, and Gly-Leu-Pro-Leu (GLPL) motifs within the NBS domain of plant resistance gene (R gene) products. Thirty-four of 65 cloned PCR fragments contained a continuous open reading frame (ORF) and their predicted protein products showed homology to the NBS-LRR class R proteins in the GenBank database. These 34 predicted protein sequences exhibited a wide range (19.5--99.4%) of sequence identity among them and were classified into two distinct groups by phylogenetic analysis. The first group consisted of 23 sequences and seemed to belong to the non-TIR NBS-LRR RGAs, since they contained group specific motifs (RNBS-A-non-TIR motif) that are often present in the coiled-coil domain of the non-TIR NBS-LRR class R genes. The second group comprised 11 sequences that contained motifs found in the TIR domain of TIR NBS-LRR class R genes. Restriction fragment length polymorphic (RFLP) markers were developed from some of the RGAs and used for mapping powdery mildew resistance genes in chestnut rose. Three markers, RGA 22 C, RGA 4 A, and RGA 7 B, were identified to be linked to a resistance gene locus, designated CRPM 1 for chestnut rose powdery mildew resistance 1, which accounted for 72% of the variation in powdery mildew resistance phenotype in an F1 segregating population. To our knowledge, this is the first report on isolation, phylogenetic analysis and potential utilization as genetic markers of RGAs in chestnut rose.

Resistance gene analogue markers are mapped to homeologous chromosomes in cultivated tetraploid cotton

Degenerate primers designed from conserved motifs of known plant resistance gene products were used to amplify genomic DNA sequences from the root-knot nematode (Meloidogyne incognita) resistance genetic source, Upland cotton (Gossypium hirsutum) cultivar Auburn 634 RNR. A total of 165 clones were isolated, and sequence analysis revealed 57 of the clones to be novel nucleotide sequences, many containing the resistance (R)-protein nucleotide-binding site motif. A cluster analysis was performed with resistance gene analogue (RGA) nucleotide sequences isolated in this study, in addition to 99 cotton RGA nucleotide sequences already deposited in GenBank, to generate a phylogenetic tree of cotton R genes. The cotton RGA nucleotide sequences were arranged into 11 groups and 56 sub-groups, based on genetic distances. Multiple sequence alignments were performed on the RGA sequences of each sub-group, and either the consensus sequences or individual RGA sequences were used to design 61 RGA-sequence-tagged site primers. A recombinant inbred line (RIL) population of cultivated tetraploid cotton was genotyped using RGA-specific primers that amplified polymorphic fragments between the two RIL parents. Nine RGA markers were mapped to homeologous chromosomes 12 and 26, based on linkage to existing markers that are located on these chromosomes.

Resistance of cotton towards Xanthomonas campestris pv. malvacearum

Interactions between Gossypium spp. and the bacterial pathogen Xanthomonas campestris pv. malvacearum are understood in the context of the gene-for-gene concept. Reviewed here are the genetic basis for cotton resistance, with reference to resistance genes, resistance gene analogs, and bacterial avirulence genes, together with the physiological mechanisms involved in the hypersensitive response to the pathogen, including production of signaling hormones, synthesis of antimicrobial molecules and alteration of host cell structures. This host-pathogen interaction represents the most complex resistance gene/avr gene system yet known and is one of the few in which phytoalexins are known to be specifically localized in HR cells at anti-microbial concentrations.