A novel root-specific gene, MIC-3, with increased expression in nematode-resistant cotton (Gossypium hirsutum L.) after root-knot nematode infection

Transcriptome Analysis of Resistant Cotton Germplasm Responding to Reniform Nematodes

Reniform nematode (Rotylenchulus reniformis) is an important microparasite for Upland cotton (Gossypium hirsutum L.) production. Growing resistant cultivars is the most economical management method, but only a few G. barbadense genotypes and some diploid Gossypium species confer high levels of resistance. This study conducted a transcriptome analysis of resistant genotypes to identify genes involved in host plant defense. Seedlings of G. arboreum accessions PI 529728 (A2-100) and PI 615699 (A2-190), and G. barbadense genotypes PI 608139 (GB 713) and PI 163608 (TX 110), were inoculated with the reniform nematode population MSRR04 and root samples were collected on the fifth (D5) and ninth (D9) day after inoculation. Differentially expressed genes (DEGs) were identified by comparing root transcriptomes from inoculated plants with those from non-inoculated plants. Accessions A2-100 and A2-190 showed 52 and 29 DEGs on D5, respectively, with 14 DEGs in common, and 18 DEGs for A2-100 and 11 DEGs for A2-190 on chromosome 5. On D9, four DEGs were found in A2-100 and two DEGs in A2-190. For GB 713, 52 and 43 DEGs were found, and for TX 110, 29 and 117 DEGs were observed on D5 and D9, respectively. Six DEGs were common at the two sampling times for these genotypes. Some DEGs were identified as Meloidogyne-induced cotton (MIC) 3 and 4, resistance gene analogs, or receptor-like proteins. Other DEGs have potential roles in plant defense, such as peroxidases, programmed cell death, pathogenesis related proteins, and systemic acquired resistance. Further research on these DEGs will aid in understanding the mechanisms of resistance to explore new applications for the development of resistant cultivars.

MicroRNA-target gene responses to root knot nematode (Meloidogyne incognita) infection in cotton (Gossypium hirsutum L.)

MicroRNAs (miRNAs) are a large class of small regulatory RNA molecules, however no study has been performed to elucidate the role of miRNAs in cotton (Gossypium hirsutum) response to the root knot nematode (RKN, Meloidogyne incognita) infection. We selected 28 miRNAs and 8 miRNA target genes to investigate the miRNA-target gene response to M. incognita infection. Our results show that RKN infection significantly affected the expression of several miRNAs and their targeted genes. After 10 days of RKN infection, expression fold changes on miRNA expressions ranged from down-regulated by 33% to upregulated by 406%; meanwhile the expression levels of miRNA target genes were 45.8% to 231%. Three miRNA-target pairs, miR159-MYB, miR319-TCP4 and miR167-ARF8, showed inverse expression patterns between gene targets and their corresponded miRNAs, suggesting miRNA-mediated gene regulation in cotton roots in response to RKN infection.

Fine mapping and identification of candidate genes for a QTL affecting Meloidogyne incognita reproduction in Upland cotton

Background

The southern root-knot nematode (Meloidogyne incognita; RKN) is one of the most important economic pests of Upland cotton (Gossypium hirsutum L.). Host plant resistance, the ability of a plant to suppress nematode reproduction, is the most economical, practical, and environmentally sound method to provide protection against this subterranean pest. The resistant line Auburn 623RNR and a number of elite breeding lines derived from it remain the most important source of root-knot nematode (RKN) resistance. Prior genetic analysis has identified two epistatically interacting RKN resistance QTLs, qMi-C11 and qMi-C14, affecting gall formation and RKN reproduction, respectively.

Results

We developed a genetic population segregating only for the qMi-C14 locus and evaluated the genetic effects of this QTL on RKN resistance in the absence of the qMi-C11 locus. The qMi-C14 locus had a LOD score of 12 and accounted for 24.5 % of total phenotypic variation for egg production. In addition to not being significantly associated with gall formation, this locus had a lower main effect on RKN reproduction than found in our previous study, which lends further support to evidence of epistasis with qMi-C11 in imparting RKN resistance in the Auburn 623RNR source. The locus qMi-C14 was fine-mapped with the addition of 16 newly developed markers. By using the reference genome sequence of G. raimondii, we identified 20 candidate genes encoding disease resistance protein homologs in the newly defined 2.3 Mb region flanked by two SSR markers. Resequencing of an RKN resistant and susceptible G. hirsutum germplasm revealed non-synonymous mutations in only four of the coding regions of candidate genes, and these four genes are consequently of high interest.

Conclusions

Our mapping results validated the effects of the qMi-C14 resistance locus, delimiting the QTL to a smaller region, and identified tightly linked SSR markers to improve the efficiency of marker-assisted selection. The candidate genes identified warrant functional studies that will help in identifying and characterizing the actual qMi-C14 defense gene(s) against root-knot nematodes.

Electronic supplementary material

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

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.

Clustering, haplotype diversity and locations of MIC-3: a unique root-specific defense-related gene family in Upland cotton (Gossypium hirsutum L.)

MIC-3 is a recently identified gene family shown to exhibit increased root-specific expression following nematode infection of cotton plants that are resistant to root-knot nematode. Here, we cloned and sequenced MIC-3 genes from selected diploid and tetraploid cotton species to reveal sequence differences at the molecular level and identify chromosomal locations of MIC-3 genes in Gossypium species. Detailed sequence analysis and phylogenetic clustering of MIC-3 genes indicated the presence of multiple MIC-3 gene members in Gossypium species. Haplotypes of a MIC-3 gene family member were discovered by comparative analysis among consensus sequences across genotypes within an individual clade in the phylogram to overcome the problem of duplicated loci in the tetraploid cotton. Deficiency tests of the SNPs delimited six A(t)-genome members of the MIC-3 family clustered to chromosome arm 4sh, and one D(t)-genome member to chromosome 19. Clustering was confirmed by long-PCR amplification of the intergenic regions using A(t)-genome-specific MIC-3 primer pairs. The clustered distribution may have been favored by selection for responsiveness to evolving disease and/or pest pressures, because large variants of the MIC-3 gene family may have been recovered from small physical areas by recombination. This could give a buffer against selection pressure from a broad range of pest and pathogens in the future. To our knowledge, these are the first results on the evolution of clustering and genome-specific haplotype members of a unique cotton gene family associated with resistant response against a major pathogen.

Phenotypic and molecular evaluation of cotton hairy roots as a model system for studying nematode resistance

Agrobacterium rhizogenes-induced cotton (Gossypium hirsutum L.) hairy roots were evaluated as a model system for studying molecular cotton-nematode interactions. Hairy root cultures were developed from the root-knot nematode (RKN) (Meloidogyne incognita [Kofoid and White] Chitwood, race 3)-resistant breeding line M315 and from the reniform nematode (RN) (Rotylenchulus reniformis Linford & Oliveira)-resistant accession GB713 (G. barbadense L.) and compared to a nematode-susceptible culture derived from the obsolete cultivar DPL90. M315, GB713, and DPL90 hairy roots differed significantly in their appearance and growth potential; however, these differences were not correlated with transcript levels of the A. rhizogenes T-DNA genes rolB and aux2 which help regulate hairy root initiation and proliferation. DPL90 hairy roots were found to support both RKN and RN reproduction in tissue culture, whereas M315 and GB713 hairy roots were resistant to RKN and RN, respectively. M315 hairy roots showed constitutive up-regulation of the defense gene MIC3 (Meloidogyne Induced Cotton3) compared to M315 whole-plant roots and DPL90 hairy roots. Our data show the potential use of cotton hairy roots in maintaining monoxenic RKN and RN cultures and suggest hairy roots may be useful in evaluating the effect of manipulated host gene expression on nematode resistance in cotton.

Molecular characterization and temporal expression analyses indicate that the MIC (Meloidogyne Induced Cotton) gene family represents a novel group of root-specific defense-related genes in upland cotton (Gossypium hirsutum L.)

The molecular events underlying the resistance of Upland cotton (Gossypium hirsutum L.) to the root-knot nematode (RKN) are largely unknown. In this report, we further characterize the previously identified MIC3 gene including the identification of 14 related MIC cDNAs in nematode-infected roots of allotetraploid cotton that show >85% identity with MIC3. A time-course analysis of RKN infection in resistant and susceptible cotton lines showed that maximum MIC transcript accumulation occurred immediately prior to the phenotypic manifestation of resistance. MIC expression was not induced by mechanical wounding or by virulent reniform nematode infection. MIC expression was undetectable in cotton leaves undergoing a hypersensitive response to Xanthomonas campestris. A time-course analysis of defense gene expression (PR10, ERF5, CDNS, LOX1, POD4, POD8) in resistant and susceptible cotton roots showed that RKN infection specifically elicits the induction of MIC in resistant roots and not other common defense-signaling pathways. These results suggest that cotton resistance to RKN involves novel defense-signaling pathways and further supports the idea that the MIC genes are intimately involved in this resistance response and represent a group of root-specific defense-related genes in cotton.