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.     

Comparative Genome-Wide Analysis of the Malate Dehydrogenase Gene Families in Cotton

Malate dehydrogenases (MDHs) play crucial roles in the physiological processes of plant growth and development. In this study, 13 and 25 MDH genes were identified from Gossypium raimondii and Gossypium hirsutum, respectively. Using these and 13 previously reported Gossypium arboretum MDH genes, a comparative molecular analysis between identified MDH genes from G. raimondii, G. hirsutum, and G. arboretum was performed. Based on multiple sequence alignments, cotton MDHs were divided into five subgroups: mitochondrial MDH, peroxisomal MDH, plastidial MDH, chloroplastic MDH and cytoplasmic MDH. Almost all of the MDHs within the same subgroup shared similar gene structure, amino acid sequence, and conserved motifs in their functional domains. An analysis of chromosomal localization suggested that segmental duplication played a major role in the expansion of cotton MDH gene families. Additionally, a selective pressure analysis indicated that purifying selection acted as a vital force in the evolution of MDH gene families in cotton. Meanwhile, an expression analysis showed the distinct expression profiles of GhMDHs in different vegetative tissues and at different fiber developmental stages, suggesting the functional diversification of these genes in cotton growth and fiber development. Finally, a promoter analysis indicated redundant but typical cis-regulatory elements for the potential functions and stress activity of many MDH genes. This study provides fundamental information for a better understanding of cotton MDH gene families and aids in functional analyses of the MDH genes in cotton fiber development.     

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.     

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.     

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.     

Molecular markers associated with the immature fiber (im) gene affecting the degree of fiber cell wall thickening in cotton (Gossypium hirsutum L.)

Cotton fiber fineness and maturity measured indirectly as micronaire (MIC) are important properties of determining fiber grades in the textile market. To understand the genetic control and molecular mechanisms of fiber fineness and maturity, we studied two near isogenic lines, Gossypium hirsutum, Texas Marker-1 wild type (TM-1) and immature fiber (im) mutant showing a significant difference in MIC values. The fibers from im mutant plants were finer and less mature with lower MIC values than those from the recurrent parent, TM-1. A comprehensive fiber property analysis of TM-1 and im mutant showed that the lower MIC of fibers in im mutant was due to the lower degree of fiber cell wall thickening as compared to the TM-1 fibers. Using an F₂ population comprising 366 progenies derived from a cross between TM-1 and im mutant, we confirmed that the immature fiber phenotype present in a mutant plant was controlled by one single recessive gene im. Furthermore, we identified 13 simple sequence repeat markers that were closely linked to the im gene located on chromosome 3. Molecular markers associated with the im gene will lay the foundation to further investigate genetic information required for improving cotton fiber fineness and maturity.     

DNA Sequence Evolution and Rare Homoeologous Conversion in Tetraploid Cotton

Allotetraploid cotton species are a vital source of spinnable fiber for textiles. The polyploid nature of the cotton genome raises many evolutionary questions as to the relationships between duplicated genomes. We describe the evolution of the cotton genome (SNPs and structural variants) with the greatly improved resolution of 34 deeply re-sequenced genomes. We also explore the evolution of homoeologous regions in the A_T- and D_T-genomes and especially the phenomenon of conversion between genomes. We did not find any compelling evidence for homoeologous conversion between genomes. These findings are very different from other recent reports of frequent conversion events between genomes. We also identified several distinct regions of the genome that have been introgressed between G. hirsutum and G. barbadense, which presumably resulted from breeding efforts targeting associated beneficial alleles. Finally, the genotypic data resulting from this study provides access to a wealth of diversity sorely needed in the narrow germplasm of cotton cultivars.     

Sampling nucleotide diversity in cotton

Background  

Cultivated cotton is an annual fiber crop derived mainly from two perennial species, Gossypium hirsutum L. or upland cotton, and G. barbadense L., extra long-staple fiber Pima or Egyptian cotton. These two cultivated species are among five allotetraploid species presumably derived monophyletically between G. arboreum and G. raimondii. Genomic-based approaches have been hindered by the limited variation within species. Yet, population-based methods are being used for genome-wide introgression of novel alleles from G. mustelinum and G. tomentosum into G. hirsutum using combinations of backcrossing, selfing, and inter-mating. Recombinant inbred line populations between genetics standards TM-1, (G. hirsutum) × 3-79 (G. barbadense) have been developed to allow high-density genetic mapping of traits.     

QTL Mapping for Fiber and Yield Traits in Upland Cotton under Multiple Environments

A population of 178 recombinant inbred lines (RILs) was developed using a single seed descendant from a cross between G. hirsutum. acc DH962 and G. hirsutum. cv Jimian5, was used to construct a genetic map and to map QTL for fiber and yield traits. A total of 644 polymorphic loci were used to construct a final genetic map, containing 616 loci and spanning 2016.44 cM, with an average of 3.27 cM between adjacent markers. Statistical analysis revealed that segregation distortion in the intraspecific population was more serious than that in the interspecific population. The RIL population and the two parents were phenotyped under 8 environments (two locations, six years), revealing a total of 134 QTL, including 64 for fiber qualities and 70 for yield components, independently detected in seven environments, explaining 4.40–15.28% of phenotypic variation (PV). Among the 134 QTL, 9 common QTL were detected in more than one environment, and 22 QTL and 19 new QTL were detected in combined analysis (E9). A total of 26 QTL hotspot regions were observed on 13 chromosomes and 2 larger linkage groups, and some QTL clusters related to fiber qualities or yield components were also observed. The results obtained in the present study suggested that to map accurate QTL in crops with larger plant types, such as cotton, phenotyping under multiple environments is necessary to effectively apply the obtained results in molecular marker-assisted selection breeding and QTL cloning.     

Chromosomal Locations of 5S and 45S rDNA in Gossypium Genus and Its Phylogenetic Implications Revealed by FISH

We investigated the locations of 5S and 45S rDNA in Gossypium diploid A, B, D, E, F, G genomes and tetraploid genome (AD) using multi-probe fluorescent in situ hybridization (FISH) for evolution analysis in Gossypium genus. The rDNA numbers and sizes, and synteny relationships between 5S and 45S were revealed using 5S and 45S as double-probe for all species, and the rDNA-bearing chromosomes were identified for A, D and AD genomes with one more probe that is single-chromosome-specific BAC clone from G. hirsutum (A₁D₁). Two to four 45S and one 5S loci were found in diploid-species except two 5S loci in G . incanum (E₄), the same as that in tetraploid species. The 45S on the 7th and 9th chromosomes and the 5S on the 9th chromosomes seemed to be conserved in A, D and AD genomes. In the species of B, E, F and G genomes, the rDNA numbers, sizes, and synteny relationships were first reported in this paper. The rDNA pattern agrees with previously reported phylogenetic history with some disagreements. Combined with the whole-genome sequencing data from G . raimondii (D₅) and the conserved cotton karyotype, it is suggested that the expansion, decrease and transposition of rDNA other than chromosome rearrangements might occur during the Gossypium evolution.     

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.     

Down-regulating annexin gene GhAnn2 inhibits cotton fiber elongation and decreases Ca2+ influx at the cell apex

Cotton fiber is a single cell that differentiates from the ovule epidermis and undergoes synchronous elongation with high secretion and growth rate. Apart from economic importance, cotton fiber provides an excellent single-celled model for studying mechanisms of cell-growth. Annexins are Ca²⁺- and phospholipid-binding proteins that have been reported to be localized in multiple cellular compartments and involved in control of vesicle secretions. Although several annexins have been found to be highly expressed in elongating cotton fibers, their functional roles in fiber development remain unknown. Here, 14 annexin family members were identified from the fully sequenced diploid G. raimondii (D₅ genome), half of which were expressed in fibers of the cultivated tetraploid species G. hirsutum (cv. YZ1). Among them, GhAnn2 from the D genome of the tetraploid species displayed high expression level in elongating fiber. The expression of GhAnn2 could be induced by some phytohormones that play important roles in fiber elongation, such as IAA and GA₃. RNAi-mediated down-regulation of GhAnn2 inhibited fiber elongation and secondary cell wall synthesis, resulting in shorter and thinner mature fibers in the transgenic plants. Measurement with non-invasive scanning ion-selective electrode revealed that the rate of Ca²⁺ influx from extracellular to intracellular was decreased at the fiber cell apex of GhAnn2 silencing lines, in comparison to that in the wild type. These results indicate that GhAnn2 may regulate fiber development through modulating Ca²⁺ fluxes and signaling.     

Analysis of [Gossypium capitis-viridis × (G.hirsutum × G.australe)2] Trispecific Hybrid and Selected Characteristics

Speciation is always a contentious and challenging issue following with the presence of gene flow. In Gossypium, there are many valuable resources and wild diploid cotton especially C and B genome species possess some excellent traits which cultivated cotton always lacks. In order to explore character transferring rule from wild cotton to upland tetraploid cotton, the [G. capitis-viridis × (G. hirsutum × G. australe)²] triple hybrid was synthesized by interspecies hybridization and chromosome doubling. Morphology comparisons were measured among this hybrid and its parents. It showed that trispecific hybrid F₁ had some intermediate morphological characters like leaf style between its parents and some different characters from its parents, like crawl growth characteristics and two kind flower color. It is highly resistant to insects comparing with other cotton species by four year field investigation. By cytogenetic analysis, triple hybrid was further confirmed by meiosis behavior of pollen mother cells. Comparing with regular meiosis of its three parents, it was distinguished by the occurrence of polyads with various numbers of unbalanced microspores and finally generating various abnormal pollen grains. All this phenomenon results in the sterility of this hybrid. This hybrid was further identified by SSR marker from DNA molecular level. It showed that 98 selected polymorphism primers amplified effective bands in this hybrids and its parents. The genetic proportion of three parents in this hybrid is 47.8% from G. hirsutum, 14.3% from G. australe, 7.0% from G. capitis-viridis, and 30.9% recombination bands respectively. It was testified that wild genetic material has been transferred into cultivated cotton and this new germplasm can be incorporated into cotton breeding program.     

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.