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 D_t (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.     

Comprehensive Analysis of the COBRA-Like (COBL) Gene Family in Gossypium Identifies Two COBLs Potentially Associated with Fiber Quality

COBRA-Like (COBL) genes, which encode a plant-specific glycosylphosphatidylinositol (GPI) anchored protein, have been proven to be key regulators in the orientation of cell expansion and cellulose crystallinity status. Genome-wide analysis has been performed in A. thaliana, O. sativa, Z. mays and S. lycopersicum, but little in Gossypium. Here we identified 19, 18 and 33 candidate COBL genes from three sequenced cotton species, diploid cotton G. raimondii, G. arboreum and tetraploid cotton G. hirsutum acc. TM-1, respectively. These COBL members were anchored onto 10 chromosomes in G. raimondii and could be divided into two subgroups. Expression patterns of COBL genes showed highly developmental and spatial regulation in G. hirsutum acc. TM-1. Of them, GhCOBL9 and GhCOBL13 were preferentially expressed at the secondary cell wall stage of fiber development and had significantly co-upregulated expression with cellulose synthase genes GhCESA4, GhCESA7 and GhCESA8. Besides, GhCOBL9 D_t and GhCOBL13 D_t were co-localized with previously reported cotton fiber quality quantitative trait loci (QTLs) and the favorable allele types of GhCOBL9 D_t had significantly positive correlations with fiber quality traits, indicating that these two genes might play an important role in fiber development.     

Genome-wide identification and characterization of HSP70 gene family in four species of cotton

Heat shock proteins (HSPs) are important elements of the cellular group of molecular chaperones. Specifically, HSP70 proteins protect cells from being damaged when plants are exposed to environmental stresses. These proteins are catalysts that manage the correct folding of other proteins, and they play a key role in the development of tolerance against biotic and abiotic stresses. In the present study, 113 HSP70 genes were retrieved from the available genome assemblies of four cotton species, including Gossypium hirsutum, G. barbadense, G. arboreum, and G. raimondii. The HSP70 genes were clustered into 11 subfamilies based on phylogeny. One hundred and nine (109) gene duplications were found across these four species. Localization of genes revealed that several HSP70 genes reside in the cytoplasm. Synonymous and non-synonymous substitution rates revealed that functional segregation of HSP70 genes in cotton is due to purifying selection. Furthermore, HSP70 genes in cotton are expressed constitutively during developmental stages. These findings are valuable to understand the complex mechanism of HSP70 gene regulation that occurs in signaling pathways in response to plant stress.     

Identification and expression analysis of phosphatidy ethanolamine-binding protein (PEBP) gene family in cotton

The phosphatidy ethanolamine-binding proteins (PEBP) play an important role in controlling flower development and phase change. Here, a total of 61 PEBP genes were identified, in which 20, 21, 10, and 10 were from tetraploid Gossypium hirsutum (AD1) and G. barbadense (AD2), and diploid G. raimondii (D5) and G. arboreum (A2), respectively. In G. hirsutum, 20 identified PEBP genes were unevenly distributed on 12 chromosomes. The identified PEBP genes were classified into four groups (TFL1, MFT, FT and FT-like). Among those, FT-like group are unique to cotton. The majority of PEBP genes had similar intron/exon distribution, whereas the divergence of PEBP genes suggests the possibility of functional diversification. The expression of PEBP genes varied among different tissues. This study brings new insights into the integrated genome-wide identification of PEBP genes in cotton and provides a foundation for breeding cotton cultivars with early maturation.     

<Emphasis Type="Italic">SnRK2</Emphasis> Homologs in <Emphasis Type="Italic">Gossypium</Emphasis> and <Emphasis Type="Italic">GhSnRK2.6</Emphasis> Improved Salt Tolerance in Transgenic Upland Cotton and <Emphasis Type="Italic">Arabidopsis</Emphasis>

SnRK2s are a large family of plant-specific protein kinases, which play important roles in multiple abiotic stress responses in various plant species. But the family in Gossypium has not been well studied. Here, we identified 13, 10, and 13 members of the SnRK2 family from Gossypium raimondii, Gossypium arboreum, and Gossypium hirsutum, respectively, and analyzed the locations of SnRK2 homologs in chromosomes based on genome data of cotton species. Phylogenetic tree analysis of SnRK2 proteins showed that these families were classified into three groups. All SnRK2 genes were comprised of nine exons and eight introns, and the exon distributions and the intron phase of homolog genes among different cotton species were analogous. Moreover, GhSnRK2.6 was overexpressed in Arabidopsis and upland cotton, respectively. Under salt treatment, overexpressed Arabidopsis could maintain higher biomass accumulation than wild-type plants, and GhSnRK2.6 overexpression in cotton exhibited higher germination rate than the control. So, the gene GhSnRK2.6 could be utilized in cotton breeding for salt tolerance.     

A New Synthetic Amphiploid (AADDAA) between Gossypium hirsutum and G. arboreum Lays the Foundation for Transferring Resistances to Verticillium and Drought

Gossypium arboreum, a cultivated cotton species (2n = 26, AA) native to Asia, possesses invaluable characteristics unavailable in the tetraploid cultivated cotton gene pool, such as resistance to pests and diseases and tolerance to abiotic stresses. However, it is quite difficult to transfer favorable traits into Upland cotton through conventional methods due to the cross-incompatibility of G. hirsutum (2n = 52, AADD) and G. arboreum. Here, we improved an embryo rescue technique to overcome the cross-incompatibility between these two parents for transferring favorable genes from G. arboreum into G. hirsutum. Our results indicate that MSB2K supplemented with 0.5 mgl^-1 kinetin and 250 mg^-1 casein hydrolysate is an efficient initial medium for rescuing early (3 d after pollination) hybrid embryos. Eight putative hybrids were successfully obtained, which were further verified and characterized by cytology, molecular markers and morphological analysis. The putative hybrids were subsequently treated with different concentrations of colchicine solution to double their chromosomes. The results demonstrate that four putative hybrid plants were successfully chromosome-doubled by treatment with 0.1% colchicine for 24 h and become amphiploid, which were confirmed by cytological observation, self-fertilization and backcrossing. Preliminary assessments of resistance at seedling stage indicate that the synthetic amphiploid showed highly resistant to Verticillium and drought. The synthetic amphiploid between G. hirsutum × G. arboreum would lay the foundation for developing G. arboreum-introgressed lines with the uniform genetic background of G. hirsutum acc TM-1, which would greatly enhance and simplify the mining, isolation, characterization, cloning and use of G. arboreum-specific desirable genes in future cotton breeding programs.     

Genome-wide identification and analysis of the aldehyde dehydrogenase (ALDH) gene superfamily of Gossypium raimondii

Background

Aldehyde dehydrogenases (ALDHs) are members of the NAD(P)⁺-dependent protein superfamily which catalyzes aliphatic and aromatic aldehyde oxidation to non-toxic carboxylic acids. ALDH genes may offer promise for improving plant adaptation to environmental stress. Recently, elucidated genome sequences of Gossypium raimondii provide a foundation for systematic identification and analysis of ALDH genes. To date, this has been accomplished for many plant species except G. raimondii.

Results

In this study, thirty unique ALDH sequences that code for 10 ALDH families were identified in the G. raimondii genome. Phylogenetic analysis revealed that ALDHs were split into six clades in G. raimondii, and ALDH proteins from the same families were clustered together. Phylogenetic relationships of ALDHs from 11 plant species suggest that ALDHs in G. raimondii shared the highest protein homology with ALDHs from poplar. Members within ALDH families possessed homologous exon–intron structures. Chromosomal distribution of ALDH did not occur evenly in the G. raimondii genome and many ALDH genes were involved in the syntenic region as documented by identification of physical locations among single chromosomes. In addition, syntenic analysis revealed that homologues of many G. raimondii ALDHs appeared in corresponding Arabidopsis and poplar syntenic blocks, indicating that these genes arose prior to G. raimondii, Arabidopsis and poplar speciation. Finally, based on gene expression analysis of microarray and RNA-seq, we can speculate that some G. raimondii ALDH genes might respond to drought or waterlogging stresses.

Conclusion

Genome-wide identification and analysis of the evolution and expression of ALDH genes in G. raimondii laid a foundation for studying this gene superfamily and offers new insights into the evolution history and speculated roles in Gossypium. These data can be used to inform functional genomic studies and molecular breeding in cotton.     

Expansion of MIR482/2118 by a class‐II transposable element in cotton

Some plant microRNA (miRNA) families contain multiple members generating identical or highly similar mature miRNA variants. Mechanisms underlying the expansion of miRNA families remain elusive, although tandem and/or segmental duplications have been proposed. In this study of two tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, and their extant diploid progenitors, Gossypium arboreum and Gossypium raimondii, we investigated the gain and loss of members of the miR482/2118 superfamily, which modulates the expression of nucleotide‐binding site leucine‐rich repeat (NBS‐LRR) disease resistance genes. We found significant expansion of MIR482/2118d in G. barbadense, G. hirsutum and G. raimondii, but not in G. arboreum. Several newly expanded MIR482/2118d loci have mutated to produce different miR482/2118 variants with altered target‐gene specificity. Based on detailed analysis of sequences flanking these MIR482/2118 loci, we found that this expansion of MIR482/2118d and its derivatives resulted from an initial capture of an MIR482/2118d by a class‐II DNA transposable element (TE) in G. raimondii prior to the tetraploidization event, followed by transposition to new genomic locations in G. barbadense, G. hirsutum and G. raimondii. The ‘GosTE’ involved in the capture and proliferation of MIR482/2118d and its derivatives belongs to the PIF/Harbinger superfamily, generating a 3‐bp target site duplication upon insertion at new locations. All orthologous MIR482/2118 loci in the two diploids were retained in the two tetraploids, but mutation(s) in miR482/2118 were observed across all four species as well as in different cultivars of both G. barbadense and G. hirsutum, suggesting a dynamic co‐evolution of miR482/2118 and its NBS‐LRR targets. Our results provide fresh insights into the mechanisms contributing to MIRNA proliferation and enrich our knowledge on TEs.     

Construction of a complete set of alien chromosome addition lines from Gossypium australe in Gossypium hirsutum: morphological,cytological, and genotypic characterization

Key message

We report the first complete set of alien addition lines of G. hirsutum . The characterized lines can be used to introduce valuable traits from G. australe into cultivated cotton.

Abstract

Gossypium australe is a diploid wild cotton species (2n = 26, GG) native to Australia that possesses valuable characteristics unavailable in the cultivated cotton gene pool, such as delayed pigment gland morphogenesis in the seed and resistances to pests and diseases. However, it is very difficult to directly transfer favorable traits into cultivated cotton through conventional gene recombination due to the absence of pairing and crossover between chromosomes of G. australe and Gossypium hirsutum (2n = 52, AADD). To enhance the transfer of favorable genes from wild species into cultivated cotton, we developed a set of hirsutum–australe monosomic alien chromosome addition lines (MAAL) using a combination of morphological survey, microsatellite marker-assisted selection, and molecular cytogenetic analysis. The amphidiploid (2n = 78, AADDGG) of G. australe and G. hirsutum was consecutively backcrossed with upland cotton to develop alien addition lines of individual G. australe chromosomes in G. hirsutum. From these backcross progeny, we generated the first complete set of chromosome addition lines in cotton; 11 of 13 lines are monosomic additions, and chromosomes 7G^a and 13G^a are multiple additions. MAALs of 1G^a and 11G^a were the first to be isolated. The chromosome addition lines can be employed as bridges for the transfer of desired genes from G. australe into G. hirsutum, as well as for gene assignment, isolation of chromosome-specific probes, flow sorting and microdissection of chromosome, development of chromosome-specific ‘‘paints’’ for fluorochrome-labeled DNA fragments, physical mapping, and selective isolation and mapping of cDNAs for a particular G. australe chromosome.     

Fine-mapping qFS07.1 controlling fiber strength in upland cotton (<Emphasis Type="Italic">Gossypium hirsutum</Emphasis> L.)

Key message

qFS07.1 controlling fiber strength was fine-mapped to a 62.6-kb region containing four annotated genes. RT-qPCR and sequence of candidate genes identified an LRR RLK gene as the most likely candidate.

Abstract

Fiber strength is an important component of cotton fiber quality and is associated with other properties, such as fiber maturity, fineness, and length. Stable QTL qFS07.1, controlling fiber strength, had been identified on chromosome 7 in an upland cotton recombinant inbred line (RIL) population from a cross (CCRI35?×?Yumian1) described in our previous studies. To fine-map qFS07.1, an F₂ population with 2484 individual plants from a cross between recombinant line RIL014 and CCRI35 was established. A total of 1518 SSR primer pairs, including 1062, designed from chromosome 1 of the Gossypium raimondii genome and 456 from chromosome 1 of the G. arboreum genome (corresponding to the QTL region) were used to fine-map qFS07.1, and qFS07.1 was mapped into a 62.6-kb genome region which contained four annotated genes on chromosome A07 of G. hirsutum. RT-qPCR and comparative analysis of candidate genes revealed a leucine-rich repeat protein kinase (LRR RLK) family protein to be a promising candidate gene for qFS07.1. Fine mapping and identification of the candidate gene for qFS07.1 will play a vital role in marker-assisted selection (MAS) and the study of mechanism of cotton fiber development.

     

Distribution and characterization of simple sequence repeats in Gossypium raimondii genome

Simple sequence repeats (SSRs) can be derived from the complete genome sequence. These markers are important for gene mapping as well as marker-assisted selection (MAS). To develop SSRs for cotton gene mapping, we selected the complete genome sequence of Gossypium raimondii, which consisted of 4447 non-redundant scaffolds. Out of 775.2 Mb sequence examined, a total of 136,345 microsatellites were identified with a density of 5.69 kb per SSR in the G. raimondii genome leading to development of 112,177 primer pairs. The distributions of SSRs in the genome were non-random. Among the different motifs ranging from 1 to 6 bp, penta-nucleotide repeats were most abundant (30.5%), followed by tetra-nucleotide repeats (18.2%) and di-nucleotide repeats (16.9%). Among all identified 457 motif types, the most frequently occurring repeat motifs were poly-AT/TA, which accounted for 79.8% of the total di-nt SSRs, followed by AAAT/TTTA with 51.5% of the total tetra-nucleotede. Further, 18,834 microsatellites were detected from the protein-coding genes, and the frequency of gene containing SSRs was 46.0% in 40,976 genes of G. raimondii. These genome-based SSRs developed in the present study will lay the groundwork for developing large numbers of SSR markers for genetic mapping, gene discovery, genetic diversity analysis, and MAS breeding in cotton.     

Characterization of 19 Genes Encoding Membrane-Bound Fatty Acid Desaturases and their Expression Profiles in Gossypium raimondii Under Low Temperature

To produce unsaturated fatty acids, membrane-bound fatty acid desaturases (FADs) can be exploited to introduce double bonds into the acyl chains of fatty acids. In this study, 19 membrane-bound FAD genes were identified in Gossypium raimondii through database searches and were classified into four different subfamilies based on phylogenetic analysis. All 19 membrane-bound FAD proteins shared three highly conserved histidine boxes, except for GrFAD2.1, which lost the third histidine box in the C-terminal region. In the G. raimondii genome, tandem duplication might have led to the increasing size of the FAD2 cluster in the Omega Desaturase subfamily, whereas segmental duplication appeared to be the dominant mechanism for the expansion of the Sphingolipid and Front-end Desaturase subfamilies. Gene expression analysis showed that seven membrane-bound FAD genes were significantly up-regulated and that five genes were greatly suppressed in G. raimondii leaves exposed to low temperature conditions.     

An integrative analysis of four CESA isoforms specific for fiber cellulose production between Gossypium hirsutum and Gossypium barbadense

Cotton fiber is an excellent model system of cellulose biosynthesis; however, it has not been widely studied due to the lack of information about the cellulose synthase (CESA) family of genes in cotton. In this study, we initially identified six full-length CESA genes designated as GhCESA5–GhCESA10. Phylogenetic analysis and gene co-expression profiling revealed that CESA1, CESA2, CESA7, and CESA8 were the major isoforms for secondary cell wall biosynthesis, whereas CESA3, CESA5, CESA6, CESA9, and CESA10 should involve in primary cell wall formation for cotton fiber initiation and elongation. Using integrative analysis of gene expression patterns, CESA protein levels, and cellulose biosynthesis in vivo, we detected that CESA8 could play an enhancing role for rapid and massive cellulose accumulation in Gossypium hirsutum and Gossypium barbadense. We found that CESA2 displayed a major expression in non-fiber tissues and that CESA1, a housekeeping gene like, was predominantly expressed in all tissues. Further, a dynamic alteration was observed in cell wall composition and a significant discrepancy was observed between the cotton species during fiber elongation, suggesting that pectin accumulation and xyloglucan reduction might contribute to cell wall transition. In addition, we discussed that callose synthesis might be regulated in vivo for massive cellulose production during active secondary cell wall biosynthesis in cotton fibers.     

A BIL Population Derived from G. hirsutum and G. barbadense Provides a Resource for Cotton Genetics and Breeding

To provide a resource for cotton genetics and breeding, an interspecific hybridization between Gossypium hirsutum cv. Emian22 and G. barbadense acc. 3–79 was made. A population of 54 BILs (backcross inbred lines, BC₁F₈) was developed with the aim of transferring G. barbadense genes into G. hirsutum in order to genetically analyze these genes’ function in a G. hirsutum background and create new germplasms for breeding. Preliminary investigation of the morphological traits showed that the BILs had diverse variations in plant architecture, seed size, and fuzz color; the related traits of yield and fiber quality evaluated in 4 environments also showed abundant phenotypic variation. In order to explore the molecular diversity of the BIL population, 446 SSR markers selected at an average genetic distance of 10 cM from our interspecific linkage map were used to genotype the BIL population. A total of 393 polymorphic loci accounting for 84.4% MAF (major allele frequency) > 0.05 and 922 allele loci were detected, and the Shannon diversity index (I) was 0.417 per locus. The average introgression segment length was 16.24 cM, and an average of 29.53 segments were introgressed in each BIL line with an average background recovery of 79.8%. QTL mapping revealed 58 QTL associated with fiber quality and yield traits, and 47 favored alleles derived from the donor parent were discovered. This study demonstrated that the interspecific BIL population was enriched with much phenotypic and molecular variation which could be a resource for cotton genetics and breeding.