MPK homolog GhNTF6 was involved in cotton against Verticillium wilt by interacted with VdEPG1

Mitogen-activated protein kinases (MPKs) are important in regulating plant development and stress response. Rapid activation of MPKs in plants usually depends on its phosphorylated. In view of this situation, a phosphorylated GhNTF6 belonged to MPKs family was screened in cotton roots under Verticillium dahliae challenge by phosphoproteomics analysis. Expression of GhNTF6 in cotton plants was did not induce by V. dahliae infection, while, silencing GhNTF6 results to enhance cotton plants susceptibility to V. dahliae, overexpression - GhNTF6 enhance Arabidopsis plants survivability to V. dahliae. Moreover, the mutation of GhNTF6 at site Thr195 and Thy197 with the phosphorylation decreased the plant resistance to V. dahliae. Therefore, GhNTF6 phosphorylation is important in plants against V. dahliae. Further analysis demonstrated that GhNTF6 interacted with a V. dahliae endopolygalacturonase (VdEPG1) on the cell nucleus. We propose that GhNTF6 is a potential molecular target for improving resistance to Verticillium wilt in cotton.

Identification of Raf-Like Kinases B Subfamily Genes in Gossypium Species Revealed GhRAF42 Enhanced Salt Tolerance in Cotton

Salinity is a critical abiotic factor that significantly reduces agricultural production. Cotton is an important fiber crop and a pioneer on saline soil, hence genetic architecture that underpins salt tolerance should be thoroughly investigated. The Raf-like kinase B-subfamily (RAF) genes were discovered to regulate the salt stress response in cotton plants. However, understanding the RAFs in cotton, such as Enhanced Disease Resistance 1 and Constitutive Triple Response 1 kinase, remains a mystery. This study obtained 29, 28, 56, and 54 RAF genes from G. arboreum, G. raimondii, G. hirsutum, and G. barbadense, respectively. The RAF gene family described allopolyploidy and hybridization events in allotetraploid cotton evolutionary connections. Ka/Ks analysis advocates that cotton evolution was subjected to an intense purifying selection of the RAF gene family. Interestingly, integrated analysis of synteny and gene collinearity suggested dispersed and segmental duplication events involved in the extension of RAFs in cotton. Transcriptome studies, functional validation, and virus-induced gene silencing on salt treatments revealed that GhRAF42 is engaged in salt tolerance in upland cotton. This research might lead to a better understanding of the role of RAFs in plants and the identification of suitable candidate salt-tolerant genes for cotton breeding.

Ethnobotanical and antimicrobial activities of the Gossypium (Cotton) genus: A review

ETHNOPHARMACOLOGICAL RELEVANCE

The Malvaceae family, an important group of plants that have the Gossypium (cotton) genus has been used in folk medicine to treat microbial diseases and symptoms.

AIMS OF THE STUDY

This article aims to understand its ethnobotany expression in communities and scientific elucidation of antimicrobial activities of this genus through literature review.

MATERIALS AND METHODS

The bibliographic survey was carried out from 1999 to 2019 with keywords combinations such as "Gossypium + ethnobotanical", " Gossypium + medicinal ", "Gossypium + the biological activity" in scientific databases as Pubmed, Scopus, Web of Science, DOAJ, Scielo, Bireme.

RESULTS

After data analysis, we found that the Gossypium genus, specifically Gossypium hirsutum, G. barbadense, G. herbaceum, G. arboreum are the species most cited in the treatment of microbial diseases and symptoms in communities all over the world. In light of scientific elucidation of biological activities, the Gossypium genus has been used to treat protozoal, bacterial, fungal, and viral diseases.

CONCLUSIONS

The review demonstrated that the Gossypium genus is a promising source of biological activities against microbial diseases, especially in the treatment of protozoal diseases like malaria.

Genome wide analysis of IQD gene family in diploid and tetraploid species of cotton (Gossypium spp.)

Calmodulin (CaM) is considered as the most significant Ca²⁺ signaling messenger that mediate various biochemical and physiological reactions. IQ domain (IQD) proteins are plant specific CML/CaM calcium binding which are characterized by domains of 67 amino acids. 50, 50, 94, and 99 IQD genes were detected from G. arboreum (A2), G. raimondii (D5), G. barbadense (AD2) and G. hirsutum (AD1) respectively. Existence of more orthologous genes in cotton species than Arabidopsis, advocated that polyploidization produced new cotton specific orthologous gene clusters. Duplication of gene events depicts that IQD gene family of cotton evolution was under strong purifying selection. G. hirsutum exhibited high level synteny. GarIQD25 exhibited high expression in stem, root, flower, ovule and fiber in G. arboreum. In G. raimondii, GraIQD03 demonstrated upregulation across stem, ovule, fiber and seed. GbaIQD11 and GbaIQD62 exhibited upregulation in fiber development in G. barbadense. GhiIQD69 recognized as main candidate genes for plant parts, floral tissues, fiber and ovule development. Promotor analysis identified cis-regulatory elements were involved in plant growth and development. Overwhelmingly, present study paves the way to better understand the evolution of cotton IQD genes and lays a foundation for future investigation of IQD in cotton.

Extensive haplotypes are associated with population differentiation and environmental adaptability in Upland cotton (Gossypium hirsutum)

Three extensive eco-haplotypes associated with population differentiation and environmental adaptability in Upland cotton were identified, with A06_85658585, A08_43734499 and A06_113104285 considered the eco-loci for environmental adaptability. Population divergence is suggested to be the primary force driving the evolution of environmental adaptability in various species. Chromosome inversion increases reproductive isolation between subspecies and accelerates population divergence to adapt to new environments. Although modern cultivated Upland cotton (Gossypium hirsutum L.) has spread worldwide, the noticeable phenotypic differences still existed among cultivars grown in different areas. In recent years, the long-distance migration of cotton cultivation areas throughout China has demanded that breeders better understand the genetic basis of environmental adaptability in Upland cotton. Here, we integrated the genotypes of 419 diverse accessions, long-term environment-associated variables (EAVs) and environment-associated traits (EATs) to evaluate subgroup differentiation and identify adaptive loci in Upland cotton. Two highly divergent genomic regions were found on chromosomes A06 and A08, which likely caused by extensive chromosome inversions. The subgroups could be geographically classified based on distinct haplotypes in the divergent regions. A genome-wide association study (GWAS) also confirmed that loci located in these regions were significantly associated with environmental adaptability in Upland cotton. Our study first revealed the cause of population divergence in Upland cotton, as well as the consequences of variation in its environmental adaptability. These findings provide new insights into the genetic basis of environmental adaptability in Upland cotton, which could accelerate the development of molecular markers for adaptation to climate change in future cotton breeding.

A cotton α1,3-/4-fucosyltransferase-encoding gene, FucT4, plays an important role in cell elongation and is significantly associated with fiber quality

Fucosylation, one of the key posttranslational modifications, plays an important role in plants. It is involved in the development, signal transduction, reproduction, and disease resistance. α1,3-/4-Fucosyltransferase is responsible for transferring L-fucose from GDP-L-fucose to the N-glycan to exert fucosylational functions. However, the roles of the fucosyltransferase gene in cotton remain unknown. This study provided a comprehensive investigation of its possible functions. A genome-wide analysis identified four, four, eight, and eight FucT genes presented in the four sequenced cotton species, diploid Gossypium raimondii, G. arboreum, tetraploid G. hirsutum acc. TM-1, and G. barbadense cv. H7124, respectively. These FucTs were classified into two groups, with FucT4 homologs alone as a group. We isolated FucT4 in TM-1 and H7124, and named it GhFucT4 and GbFucT4, respectively. Quantitative RT-PCR and transcriptome data demonstrated that GhFucT4 had the highest expression levels in fibers among all GhFucT genes. Association studies and QTL co-localization supported the possible involvement of GhFucT4 in cotton fiber development. GhFucT4 and GbFucT4 shared high sequence identities, and FucT4 had higher expression in H7124 fiber tissues compared with TM-1. Furthermore, ectopic expression of FucT4 in transgenic Arabidopsis promoted root cell elongation, upregulated expression of genes related to cell wall loosening, and led to longer primary root. These results collectively indicate that FucT4 plays an important role in promoting cell elongation and modulating fiber development, which could be utilized to improve fiber quality traits in cotton breeding.

Comprehensive characterization and gene expression patterns of LBD gene family in Gossypium

A comprehensive account of the LBD gene family of Gossypium was provided in this work. Expression analysis and functional characterization revealed that LBD genes might play different roles in G. hirsutum and G. barbadense. The Lateral Organ Boundaries Domain (LBD) proteins comprise a plant-specific transcription factor family, which plays crucial roles in physiological processes of plant growth, development, and stress tolerance. In the present work, a systematical analysis of LBD gene family from two allotetraploid cotton species, G. hirsutum and G. barbadense, together with their genomic donor species, G. arboreum and G. raimondii, was conducted. There were 131, 128, 62, and 68 LBDs identified in G. hirsutum, G. barbadense, G. arboreum and G. raimondii, respectively. The LBD proteins could be classified into two main classes, class I and class II, based on the structure of their lateral organ boundaries domain and traits of phylogenetic tree, and class I was further divided into five subgroups. The gene structure and motif composition analyses conducted in both G. hirsutum and G. barbadense revealed that LBD genes kept relatively conserved within the subfamilies. Synteny analysis suggested that segmental duplication acted as an important mechanism in expansion of the cotton LBD gene family. Cis-element analysis predicated the possible functions of LBD genes. Public RNA-seq data were investigated to analyze the expression patterns of cotton LBD genes in various tissues as well as gene expression under abiotic stress treatments. Furthermore, RT-qPCR results found that GhLBDs had various expression regulation under MeJA treatments. Expression analysis indicated the differential functions of cotton LBD genes in response to abiotic stress and hormones.

Drought-induced disturbance of carbohydrate metabolism in anthers and male abortion of two Gossypium hirsutum cultivars differing in drought tolerance

Cotton pollen abortion, under drought stress, was closely associated with changes in anther carbohydrate metabolism, and pollen abortion rate due to drought was higher in drought-sensitive cultivars than drought-tolerant cultivars. Cotton reproductive failure under drought stress is intrinsically connected with altered male fertility, however, studies investigating the effect of drought stress on cotton male fertility are nonexistent. Thus, a drought stress experiment was conducted with two cotton cultivars, differing in drought tolerance, to study pollen fertility and anthers' physiology. Results indicated that drought stress reduced pollen fertility of both cultivars due to decreases in anther starch and adenosine triphosphate (ATP) synthesis. Lower assimilate supply capacity in conjunction with impaired activities of ADP-glucose pyrophosphorylase and soluble starch synthase were the main reasons for the decreased starch levels in drought-stressed anthers. The decreased activities of sucrose synthetase and acid invertase were responsible for the higher sucrose level in drought-stressed anthers than well-watered anthers and the changing trend of sucrose was intensified by the decreased expressions of sucrose synthase genes (GhSusA, GhSusB, GhSusD) and acid invertase genes (GhINV1, GhINV2). However, despite sucrose degradation being limited in drought-stressed anthers, glucose level was higher in droughted anthers than well-watered ones, and that might be attributed to the down-regulated respiration since decreased anther ATP levels were detected in drought-stressed plants. Furthermore, compared to the drought-tolerant cultivar, pollen fertility was more suppressed by drought stress for the drought-sensitive cultivar, and that was attributed to the larger decrease in starch and ATP contents.

QTL analysis and candidate gene identification for plant height in cotton based on an interspecific backcross inbred line population of Gossypium hirsutum × Gossypium barbadense

We constructed the first high-quality and high-density genetic linkage map for an interspecific BIL population in cotton by specific-locus amplified fragment sequencing for QTL mapping. A novel gene GhPIN3 for plant height was identified in cotton. Ideal plant height (PH) is important for improving lint yield and mechanized harvesting in cotton. Most published genetic studies on cotton have focused on fibre yield and quality traits rather than PH. To facilitate the understanding of the genetic basis in PH, an interspecific backcross inbred line (BIL) population of 250 lines derived from upland cotton (Gossypium hirsutum L.) CRI36 and Egyptian cotton (G. barbadense L.) Hai7124 was used to construct a high-density genetic linkage map for quantitative trait locus (QTL) mapping. The high-density genetic map harboured 7,709 genotyping-by-sequencing (GBS)-based single nucleotide polymorphism (SNP) markers that covered 3,433.24 cM with a mean marker interval of 0.67 cM. In total, ten PH QTLs were identified and each explained 4.27-14.92% of the phenotypic variation, four of which were stable as they were mapped in at least two tests or based on best linear unbiased prediction in seven field tests. Based on functional annotation of orthologues in Arabidopsis and transcriptome data for the genes within the stable QTL regions, GhPIN3 encoding for the hormone auxin efflux carrier protein was identified as a candidate gene located in the stable QTL qPH-Dt1-1 region. A qRT-PCR analysis showed that the expression level of GhPIN3 in apical tissues was significantly higher in four short-statured cotton genotypes than that in four tall-statured cotton genotypes. Virus-induced gene silencing cotton has significantly increased PH when the expression of the GhPIN3 gene was suppressed.

Population genomics reveals a fine-scale recombination landscape for genetic improvement of cotton

Recombination breaks up ancestral linkage disequilibrium, creates combinations of alleles, affects the efficiency of natural selection, and plays a major role in crop domestication and improvement. However, there is little knowledge regarding the variation in the population-scaled recombination rate in cotton. We constructed recombination maps and characterized the difference in the genomic landscape of the population-scaled recombination rate between Gossypium hirsutum and G. arboreum and sub-genomes based on the 381 sequenced G. hirsutum and 215 G. arboreum accessions. Comparative genomics identified large structural variations and syntenic genes in the recombination regions, suggesting that recombination was related to structural variation and occurred preferentially in the distal chromosomal regions. Correlation analysis indicated that recombination was only slightly affected by geographical distribution and breeding period. A genome-wide association study (GWAS) was performed with 15 agronomic traits using 267 cotton accessions and identified 163 quantitative trait loci (QTL) and an important candidate gene (Ghir_COL2) for early maturity traits. Comparative analysis of recombination and a GWAS revealed that the QTL of fibre quality traits tended to be more common in high-recombination regions than were those of yield and early maturity traits. These results provide insights into the population-scaled recombination landscape, suggesting that recombination contributed to the domestication and improvement of cotton, which provides a useful reference for studying recombination in other species.

Genome-wide comparative analysis of RNA-binding Glycine-rich protein family genes between Gossypium arboreum and Gossypium raimondii

RB-GRP (RNA-binding Glycine-rich protein gene) family belongs to the fourth subfamily of the GRP (Glycine-rich protein gene) superfamily, which plays a great role in plant growth and development, as well as in abiotic stresses response, while it has not been identified in cotton. Here, we identified 37 and 32 RB-GRPs from two cotton species (Gossypium arboreum and Gossypium raimondii, respectively), which were divided into four distinct subfamilies based on the presence of additional motifs and the arrangement of the glycine repeats. The distribution of RB-GRPs was nonrandom and uneven among the chromosomes both in two cotton species. The expansion of RB-GRP gene family between two cultivars was mainly attributed to segmental and tandem duplication events indicated by synteny analysis, and the tandem duplicated genes were mapped into homologous collinear blocks, indicated that they shared a common ancestral gene in both species. Furthermore, most RB-GRPs in two cotton species undergone stronger negative selective pressure by evolutionary analysis of RB-GRP orthologous genes. Meanwhile, RB-GRPs participated in different abiotic stresses (Abscisic acid, salt and Polyethylene glycol) responses and tissues at different developmental stages between two cotton species were showed by gene expression analysis. This research would provide insight into the evolution and function of the RB-GRPs in Gossypium species.

Different divergence events for three pairs of PEBPs in Gossypium as implied by evolutionary analysis

INTRODUCTION

The phosphatidylethanolamine-binding protein (PEBP) gene family plays a crucial role in seed germination, reproductive transformation, and other important developmental processes in plants, but its distribution in Gossypium genomes or species, evolutionary properties, and the fates of multiple duplicated genes remain unclear.

OBJECTIVES

The primary objectives of this study were to elucidate the distribution and characteristics of PEBP genes in Gossypium, as well as the evolutionary pattern of duplication and deletion, and functional differentiation of PEBPs in plants.

METHODS

Using the PEBP protein sequences in Arabidopsis thaliana as queries, blast alignment was carried out for the identification of PEBP genes in four sequenced cotton species. Using the primers designed according to the PEBP genome sequences, PEBP genes were cloned from 15 representative genomes of Gossypium genus, and the gene structure, CDS sequence, protein sequence and properties were predicted and phylogenetic analysis was performed. Taking PEBP proteins of grape as reference, grouping of orthologous gene, analysis of phylogeny and divergence of PEBPs in nine species were conducted to reconstruct the evolutionary pattern of PEBP genes in plants.

RESULTS

We identified and cloned 160 PEBPs from 15 cotton species, and the phylogenetic analysis showed that the genes could be classified into the following three subfamilies: MFT-like, FT-like and TFL1-like. There were eight single orthologous group (OG) members in each diploid and 16 double OG members in each tetraploid. An analysis of the expression and selective pressure indicated that expression divergence and strong purification selection within the same OG presented in the PEBP gene family.

CONCLUSION

An evolutionary pattern of duplication and deletion of the PEBP family in the evolutionary history of Gossypium was suggested, and three pairs of genes resulted from different whole-genome duplication events.

The response of dominant and rare taxa for fungal diversity within different root environments to the cultivation of Bt and conventional cotton varieties

BACKGROUND

Bacillus thuringiensis (Bt) crops have been cultivated at a large scale over the past several decades, which have raised concern about unintended effects on natural environments. Microbial communities typically contain numerous rare taxa that make up the majority of community populations. However, the response of dominant and rare taxa for fungal diversity to the different root environments of Bt plants remains unclear.

RESULTS

We quantified fungal population sizes and community composition via quantitative PCR of ITS genes and 18S rRNA gene sequencing of, respectively, that were associated with Bt and conventional cotton variety rhizosphere soils from different plant growth stages. qPCR analyses indicated that fungal abundances reached their peak at the seedling stage and that the taproots and lateral root rhizospheres of the Bt cotton SGK321 were significantly different. However, no significant differences in population sizes were detected between the same root zones from Bt and the conventional cotton varieties. The overall patterns of fungal genera abundances followed that of the dominant genera, whereas overall patterns of fungal genera richness followed those of the rare genera. These results suggest that the dominant and rare taxa play different roles in the maintenance of rhizosphere microhabitat ecosystems. Cluster analyses indicated a separation of fungal communities based on the lateral roots or taproots from the three cotton varieties at the seedling stage, suggesting that root microhabitats had marked effects on fungal community composition. Redundancy analyses indicated that pH was more correlated to soil fungal community composition than Bt protein content.

CONCLUSIONS

In conclusion, these results indicate that dominant and rare fungal taxa differentially contribute to community dynamics in different root microhabitats of both Bt and conventional cotton varieties. Moreover, these results showed that the rhizosphere fungal community of Bt cotton did not respond significantly to the presence of Bt protein when compared to the two conventional cotton varieties.

Genome-wide identification and analysis of the evolution and expression patterns of the GATA transcription factors in three species of Gossypium genus

GATA transcription factors (TFs), which bind to DNA in regulatory regions, are involved in cell differentiation and possess a type-IV zinc finger and a DNA-binding domain. GATA genes have been characterized in plant species such as Arabidopsis thaliana, Oryza sativa, and Glycine max, and their functions have been elucidated in A. thaliana. Although many Gossypium quantitative trait loci for fiber quality harbor GATA TFs, GATA genes have not yet been characterized in cotton. In this study, we identified 179 GATA genes from the genomes of three Gossypium species. We analyzed the phylogenetic relationships, chromosomal distribution, gene structure, expression pattern, and predicted promoters of all 179 Gossypium GATA genes (46 in G. raimondii, 46 in G. arboreum, and 87 in G. hirsutum). Phylogenetic analysis grouped the 179 GATA genes into four subfamilies. Domain analysis revealed that GATA domains in subfamilies I, II, and III were located near the C-terminal, whereas those in subfamily IV were adjacent to the N-terminal. RNA-seq and (Real-time PCR) qRT-PCR revealed that 39.1% (34/87) of GATA genes were expressed in growing plant tissues in G. hirsutum, but only 12.6% (11/87) were expressed during fiber development. In addition, 45.7% (21/46) and 26.1% (12/46) of GATA genes were expressed in G. arboreum and G. raimondii, respectively. Our results may be useful for elucidating the evolution, expression patterns, and functional divergence of GATA genes in Gossypium.

Genome-wide identification and characterization of LRR-RLKs reveal functional conservation of the SIF subfamily in cotton (Gossypium hirsutum)

BACKGROUND

As one of the largest subfamilies of the receptor-like protein kinases (RLKs) in plants, Leucine Rich Repeats-RLKs (LRR-RLKs) are involved in many critical biological processes including growth, development and stress responses in addition to various physiological roles. Arabidopsis contains 234 LRR-RLKs, and four members of Stress Induced Factor (SIF) subfamily (AtSIF1-AtSIF4) which are involved in abiotic and biotic stress responses. Herein, we aimed at identification and functional characterization of SIF subfamily in cultivated tetraploid cotton Gossypium hirsutum.

RESULTS

Genome-wide analysis of cotton LRR-RLK gene family identified 543 members and phylogenetic analysis led to the identification of 6 cotton LRR-RLKs with high homology to Arabidopsis SIFs. Of the six SIF homologs, GhSIF1 is highly conserved exhibiting 46-47% of homology with AtSIF subfamily in amino acid sequence. The GhSIF1 was transiently silenced using Virus-Induced Gene Silencing system specifically targeting the 3' Untranslated Region. The transiently silenced cotton seedlings showed enhanced salt tolerance compared to the control plants. Further, the transiently silenced plants showed better growth, lower electrolyte leakage, and higher chlorophyll and biomass contents.

CONCLUSIONS

Overall, 543 LRR-RLK genes were identified using genome-wide analysis in cultivated tetraploid cotton G. hirsutum. The present investigation also demonstrated the conserved salt tolerance function of SIF family member in cotton. The GhSIF1 gene can be knocked out using genome editing technologies to improve salt tolerance in cotton.

Molecular and Physio-Biochemical Characterization of Cotton Species for Assessing Drought Stress Tolerance

Drought stress significantly limits cotton growth and production due to the necessity of water at every stage of crop growth. Hence, it is essential to identify tolerant genetic resources and understand the mechanisms of drought tolerance in economically and socially important plants such as cotton. In this study, molecular and physio-biochemical investigations were conducted by analyzing different parameters by following standard protocols in three different cotton species, namely TM-1 (Gossypium hirsutum), Zhongmian-16 (Gossypium arboreum), and Pima4-S (Gossypium barbadense). Drought stress significantly decreased plant growth, chlorophyll content, net photosynthetic rate (Pn), stomatal conductance (Gs), maximum photochemical efficiency of PSII (Fv/Fm), and relative water content. TM-1 resulted in more tolerance than the other two species. The accumulation of proline, soluble proteins, soluble sugars, hydrogen peroxide (H₂O₂), and superoxide radicals (O₂•-) increased significantly in TM-1. In addition, TM-1 maintained the integrity of the chloroplast structure under drought conditions. The relative expression level of drought-responsive genes including coding for transcription factors and other regulatory proteins or enzymes controlling genes (ERF, ERFB, DREB, WRKY6, ZFP1, FeSOD, CuZnSOD, MAPKKK17, P5CR, and PRP5) were higher in TM-1 under drought, conferring a more tolerant status than in Zhongmian-16 and Pima4-S. The findings of this research could be utilized for predicting a tolerant cotton genotype as well as evaluating prospective cotton species in the variety development program.

Analysis of genetic diversity and population structure in upland cotton (Gossypium hirsutum L.) germplasm using simple sequence repeats

Improvement of cotton fibre yield and quality is challenging due to the narrow genetic base of modern cotton cultivars, which emphasizes the great need to effectively explore the existing germplasm resources. With major objective to assess the genetic diversity and population structure at DNA level, 302 elite upland cotton germplasm accessions (253 Chinese and 49 different exotic origins), were genotyped using 198 simple sequence repeats (SSRs) markers. Each of the 198 markers differed greatly in its ability to detect variations in the panel of cotton germplasm. The SSRs amplified 897 alleles, of which 77.7% were polymorphic. The number of alleles varied from 2 to 12 (mean 4.53). Gene diversity ranged from 0.020 to 0.492 with a mean of 0.279. The polymorphic information content (PIC) values ranged from 0.371 to 0.019 (mean 0.225). Genetic distances in the whole cotton germplasm ranged from 0.451 to 0.052 (mean 0.270), demonstrating relatively wider genetic diversity range. Chinese-origin cotton germplasm showed the highest level of SSR polymorphisms (gene diversity=0.268, PIC=0.218), whereas American-origin revealed the highest mean genetic distance (0.274). Model-based Bayesian analysis clustered the whole cotton germplasm into three subpopulations, and the highest molecular variation ws revealed between subpopulations (4%, P<0.001). The SSRs revealed moderate level of genetic diversity at DNA level, identified three structured subpopulations, suggesting a potential use of these markers for genomewide association mapping studies and for identifying and conserving useful alleles in upland cotton germplasm.

Genomic Comparison of the P-ATPase Gene Family in Four Cotton Species and Their Expression Patterns in Gossypium hirsutum

Plant P-type H⁺-ATPase (P-ATPase) is a membrane protein existing in the plasma membrane that plays an important role in the transmembrane transport of plant cells. To understand the variety and quantity of P-ATPase proteins in different cotton species, we combined four databases from two diploid cotton species (Gossypium raimondii and G. arboreum) and two tetraploid cotton species (G. hirsutum and G. barbadense) to screen the P-ATPase gene family and resolved the evolutionary relationships between the former cotton species. We identified 53, 51, 99 and 98 P-ATPase genes from G. arboretum, G. raimondii, G. barbadense and G. hirsutum, respectively. The structural and phylogenetic analyses revealed that the gene structure was consistent between P-ATPase genes, with a close evolutionary relationship. The expression analysis of P-ATPase genes showed that many P-ATPase genes were highly expressed in various tissues and at different fiber developmental stages in G. hirsutum, suggesting that they have potential functions during growth and fiber development in cotton.

Genome-wide analysis of codon usage bias in four sequenced cotton species

Codon usage bias (CUB) is an important evolutionary feature in a genome which provides important information for studying organism evolution, gene function and exogenous gene expression. The CUB and its shaping factors in the nuclear genomes of four sequenced cotton species, G. arboreum (A₂), G. raimondii (D₅), G. hirsutum (AD₁) and G. barbadense (AD₂) were analyzed in the present study. The effective number of codons (ENC) analysis showed the CUB was weak in these four species and the four subgenomes of the two tetraploids. Codon composition analysis revealed these four species preferred to use pyrimidine-rich codons more frequently than purine-rich codons. Correlation analysis indicated that the base content at the third position of codons affect the degree of codon preference. PR2-bias plot and ENC-plot analyses revealed that the CUB patterns in these genomes and subgenomes were influenced by combined effects of translational selection, directional mutation and other factors. The translational selection (P2) analysis results, together with the non-significant correlation between GC12 and GC3, further revealed that translational selection played the dominant role over mutation pressure in the codon usage bias. Through relative synonymous codon usage (RSCU) analysis, we detected 25 high frequency codons preferred to end with T or A, and 31 low frequency codons inclined to end with C or G in these four species and four subgenomes. Finally, 19 to 26 optimal codons with 19 common ones were determined for each species and subgenomes, which preferred to end with A or T. We concluded that the codon usage bias was weak and the translation selection was the main shaping factor in nuclear genes of these four cotton genomes and four subgenomes.

Evolution and Stress Responses of Gossypium hirsutum SWEET Genes

The SWEET (sugars will eventually be exported transporters) proteins are sugar efflux transporters containing the MtN3_saliva domain, which affects plant development as well as responses to biotic and abiotic stresses. These proteins have not been functionally characterized in the tetraploid cotton, Gossypium hirsutum, which is a widely cultivated cotton species. In this study, we comprehensively analyzed the cotton SWEET gene family. A total of 55 putative G. hirsutumSWEET genes were identified. The GhSWEET genes were classified into four clades based on a phylogenetic analysis and on the examination of gene structural features. Moreover, chromosomal localization and an analysis of homologous genes in Gossypium arboreum, Gossypium raimondii, and G. hirsutum suggested that a whole-genome duplication, several tandem duplications, and a polyploidy event contributed to the expansion of the cotton SWEET gene family, especially in Clade III and IV. Analyses of cis-acting regulatory elements in the promoter regions, expression profiles, and artificial selection revealed that the GhSWEET genes were likely involved in cotton developmental processes and responses to diverse stresses. These findings may clarify the evolution of G. hirsutum SWEET gene family and may provide a foundation for future functional studies of SWEET proteins regarding cotton development and responses to abiotic stresses.

Identification of the group IIa WRKY subfamily and the functional analysis of GhWRKY17 in upland cotton (Gossypium hirsutum L.)

WRKY transcription factors play important roles in plant defense, stress response, leaf senescence, and plant growth and development. Previous studies have revealed the important roles of the group IIa GhWRKY genes in cotton. To comprehensively analyze the group IIa GhWRKY genes in upland cotton, we identified 15 candidate group IIa GhWRKY genes in the Gossypium hirsutum genome. The phylogenetic tree, intron-exon structure, motif prediction and Ka/Ks analyses indicated that most group IIa GhWRKY genes shared high similarity and conservation and underwent purifying selection during evolution. In addition, we detected the expression patterns of several group IIa GhWRKY genes in individual tissues as well as during leaf senescence using public RNA sequencing data and real-time quantitative PCR. To better understand the functions of group IIa GhWRKYs in cotton, GhWRKY17 (KF669857) was isolated from upland cotton, and its sequence alignment, promoter cis-acting elements and subcellular localization were characterized. Moreover, the over-expression of GhWRKY17 in Arabidopsis up-regulated the senescence-associated genes AtWRKY53, AtSAG12 and AtSAG13, enhancing the plant’s susceptibility to leaf senescence. These findings lay the foundation for further analysis and study of the functions of WRKY genes in cotton.

Characterization of the late embryogenesis abundant (LEA) proteins family and their role in drought stress tolerance in upland cotton

BACKGROUND

Late embryogenesis abundant (LEA) proteins are large groups of hydrophilic proteins with major role in drought and other abiotic stresses tolerance in plants. In-depth study and characterization of LEA protein families have been carried out in other plants, but not in upland cotton. The main aim of this research work was to characterize the late embryogenesis abundant (LEA) protein families and to carry out gene expression analysis to determine their potential role in drought stress tolerance in upland cotton. Increased cotton production in the face of declining precipitation and availability of fresh water for agriculture use is the focus for breeders, cotton being the backbone of textile industries and a cash crop for many countries globally.

RESULTS

In this work, a total of 242, 136 and 142 LEA genes were identified in G. hirsutum, G. arboreum and G. raimondii respectively. The identified genes were classified into eight groups based on their conserved domain and phylogenetic tree analysis. LEA 2 were the most abundant, this could be attributed to their hydrophobic character. Upland cotton LEA genes have fewer introns and are distributed in all chromosomes. Majority of the duplicated LEA genes were segmental. Syntenic analysis showed that greater percentages of LEA genes are conserved. Segmental gene duplication played a key role in the expansion of LEA genes. Sixty three miRNAs were found to target 89 genes, such as miR164, ghr-miR394 among others. Gene ontology analysis revealed that LEA genes are involved in desiccation and defense responses. Almost all the LEA genes in their promoters contained ABRE, MBS, W-Box and TAC-elements, functionally known to be involved in drought stress and other stress responses. Majority of the LEA genes were involved in secretory pathways. Expression profile analysis indicated that most of the LEA genes were highly expressed in drought tolerant cultivars Gossypium tomentosum as opposed to drought susceptible, G. hirsutum. The tolerant genotypes have a greater ability to modulate genes under drought stress than the more susceptible upland cotton cultivars.

CONCLUSION

The finding provides comprehensive information on LEA genes in upland cotton, G. hirsutum and possible function in plants under drought stress.

Comparative Genomics of an Unusual Biogeographic Disjunction in the Cotton Tribe (Gossypieae) Yields Insights into Genome Downsizing

Long-distance insular dispersal is associated with divergence and speciation because of founder effects and strong genetic drift. The cotton tribe (Gossypieae) has experienced multiple transoceanic dispersals, generating an aggregate geographic range that encompasses much of the tropics and subtropics worldwide. Two genera in the Gossypieae, Kokia and Gossypioides, exhibit a remarkable geographic disjunction, being restricted to the Hawaiian Islands and Madagascar/East Africa, respectively. We assembled and use de novo genome sequences to address questions regarding the divergence of these two genera from each other and from their sister-group, Gossypium. In addition, we explore processes underlying the genome downsizing that characterizes Kokia and Gossypioides relative to other genera in the tribe. Using 13,000 gene orthologs and synonymous substitution rates, we show that the two disjuncts last shared a common ancestor ∼5 Ma, or half as long ago as their divergence from Gossypium. We report relative stasis in the transposable element fraction. In comparison to Gossypium, there is loss of ∼30% of the gene content in the two disjunct genera and a history of genome-wide accumulation of deletions. In both genera, there is a genome-wide bias toward deletions over insertions, and the number of gene losses exceeds the number of gains by ∼2- to 4-fold. The genomic analyses presented here elucidate genomic consequences of the demographic and biogeographic history of these closest relatives of Gossypium, and enhance their value as phylogenetic outgroups.

Segregation distortion and genome-wide digenic interactions affect transmission of introgressed chromatin from wild cotton species

This study reports transmission genetics of chromosomal segments into Gossypium hirsutum from its most distant euploid relative, Gossypium mustelinum . Mutilocus interactions and structural rearrangements affect introgression and segregation of donor chromatin. Wild allotetraploid relatives of cotton are a rich source of genetic diversity that can be used in genetic improvement, but linkage drag and non-Mendelian transmission genetics are prevalent in interspecific crosses. These problems necessitate knowledge of transmission patterns of chromatin from wild donor species in cultivated recipient species. From an interspecific cross, Gossypium hirsutum × Gossypium mustelinum, we studied G. mustelinum (the most distant tetraploid relative of Upland cotton) allele retention in 35 BC₃F₁ plants and segregation patterns in BC₃F₂ populations totaling 3202 individuals, using 216 DNA marker loci. The average retention of donor alleles across BC₃F₁ plants was higher than expected and the average frequency of G. mustelinum alleles in BC₃F₂ segregating families was less than expected. Despite surprisingly high retention of G. mustelinum alleles in BC₃F₁, 46 genomic regions showed no introgression. Regions on chromosomes 3 and 15 lacking introgression were closely associated with possible small inversions previously reported. Nonlinear two-locus interactions are abundant among loci with single-locus segregation distortion, and among loci originating from one of the two subgenomes. Comparison of the present results with those of prior studies indicates different permeability of Upland cotton for donor chromatin from different allotetraploid relatives. Different contributions of subgenomes to two-locus interactions suggest different fates of subgenomes in the evolution of allotetraploid cottons. Transmission genetics of G. hirsutum × G. mustelinum crosses reveals allelic interactions, constraints on fixation and selection of donor alleles, and challenges with retention of introgressed chromatin for crop improvement.

Genetic and phenotypic effects of chromosome segments introgressed from Gossypium barbadense into Gossypium hirsutum

MBI9915 is an introgression cotton line with excellent fiber quality. It was obtained by advanced backcrossing and continuous inbreeding from an interspecific cross between the upland cotton (Gossypium hirsutum) cultivar CCRI36 as the recurrent parent and the sea island cotton (G. barbadense) cultivar Hai1, as the donor parent. To study the genetic effects of the introgressed chromosome segments in G. hirsutum, an F₂ secondary segregating population of 1537 individuals was created by crossing MBI9915 and CCRI36, and an F_2:3 population was created by randomly selecting 347 individuals from the F₂ generation. Quantitative trait locus (QTL) mapping and interaction for fiber length and strength were identified using IciMapping software. The genotype analysis showed that the recovery rate for MBI9915 was 97.9%, with a total 6 heterozygous segments and 13 homozygous segments. A total of 18 QTLs for fiber quality and 6 QTLs for yield related traits were detected using the two segregating generations. These QTLs were distributed across 7 chromosomes and collectively explained 0.81%–9.51% of the observed phenotypic variations. Six QTLs were consistently detected in two generations and 6 QTLs were identified in previous studies. A total of 13 pairs of interaction for fiber length and 13 pairs of interaction for fiber strength were identified in two generations. Among them, 3 pairs of interaction for fiber length and 3 pairs of interaction for fiber strength could be identified in all generations; 4 pairs of interactions affected fiber length and fiber strength simultaneously. The results clearly showed that 5 chromosome segments (Seg-5-1, Seg-7-1, Seg-8-1, Seg-20-2 and Seg-20-3) have important effects on fiber yield and quality. This study provides the useful information for gene cloning and marker-assisted breeding for excellent fiber related quality.

Characterization of Developmental Immature Fiber ( im) Mutant and Texas Marker-1 (TM-1) Cotton Fibers Using Attenuated Total Reflection Fourier Transform Infrared (ATR FT-IR) Spectroscopy

The immature fiber ( im) mutant is one type of cotton fiber mutant with unique characteristics of non-fluffy cotton bolls. Compared to its near-isogenic wild type Texas Marker-1 (TM-1), im fiber has a thin secondary cell wall and is less mature. In this work, we applied the previously proposed principal component analysis (PCA) and simple algorithms to analyze the attenuated total reflection Fourier transform infrared (ATR FT-IR) spectra of developmental im and TM-1 fibers. The results from these approaches could not effectively and consistently indicate the inherent difference between TM-1 and im fibers at the same developmental stage. The difference between TM-1 and corresponding im fibers was detected when comparing the normalized intensity variations of the 730 cm^-1 bands. The 730 cm^-1 band intensities in developmental im fibers are temporally lower than those in developmental TM-1 fibers although they became similar when the TM-1 and im fibers are fully mature. The observation might imply the likelihood of temporal reduction of amorphous regions in developmental im fibers rather than in developmental TM-1 fibers.

Diversity analysis of cotton (Gossypium hirsutum L.) germplasm using the CottonSNP63K Array

BACKGROUND

Cotton germplasm resources contain beneficial alleles that can be exploited to develop germplasm adapted to emerging environmental and climate conditions. Accessions and lines have traditionally been characterized based on phenotypes, but phenotypic profiles are limited by the cost, time, and space required to make visual observations and measurements. With advances in molecular genetic methods, genotypic profiles are increasingly able to identify differences among accessions due to the larger number of genetic markers that can be measured. A combination of both methods would greatly enhance our ability to characterize germplasm resources. Recent efforts have culminated in the identification of sufficient SNP markers to establish high-throughput genotyping systems, such as the CottonSNP63K array, which enables a researcher to efficiently analyze large numbers of SNP markers and obtain highly repeatable results. In the current investigation, we have utilized the SNP array for analyzing genetic diversity primarily among cotton cultivars, making comparisons to SSR-based phylogenetic analyses, and identifying loci associated with seed nutritional traits.

RESULTS

The SNP markers distinctly separated G. hirsutum from other Gossypium species and distinguished the wild from cultivated types of G. hirsutum. The markers also efficiently discerned differences among cultivars, which was the primary goal when designing the CottonSNP63K array. Population structure within the genus compared favorably with previous results obtained using SSR markers, and an association study identified loci linked to factors that affect cottonseed protein content.

CONCLUSIONS

Our results provide a large genome-wide variation data set for primarily cultivated cotton. Thousands of SNPs in representative cotton genotypes provide an opportunity to finely discriminate among cultivated cotton from around the world. The SNPs will be relevant as dense markers of genome variation for association mapping approaches aimed at correlating molecular polymorphisms with variation in phenotypic traits, as well as for molecular breeding approaches in cotton.

Discovery and identification of candidate genes from the chitinase gene family for Verticillium dahliae resistance in cotton

Verticillium dahliae, a destructive and soil-borne fungal pathogen, causes massive losses in cotton yields. However, the resistance mechanism to V. dahilae in cotton is still poorly understood. Accumulating evidence indicates that chitinases are crucial hydrolytic enzymes, which attack fungal pathogens by catalyzing the fungal cell wall degradation. As a large gene family, to date, the chitinase genes (Chis) have not been systematically analyzed and effectively utilized in cotton. Here, we identified 47, 49, 92, and 116 Chis from four sequenced cotton species, diploid Gossypium raimondii (D5), G. arboreum (A2), tetraploid G. hirsutum acc. TM-1 (AD1), and G. barbadense acc. 3-79 (AD2), respectively. The orthologous genes were not one-to-one correspondence in the diploid and tetraploid cotton species, implying changes in the number of Chis in different cotton species during the evolution of Gossypium. Phylogenetic classification indicated that these Chis could be classified into six groups, with distinguishable structural characteristics. The expression patterns of Chis indicated their various expressions in different organs and tissues, and in the V. dahliae response. Silencing of Chi23, Chi32, or Chi47 in cotton significantly impaired the resistance to V. dahliae, suggesting these genes might act as positive regulators in disease resistance to V. dahliae.

Soil Potassium Deficiency Reduces Cotton Fiber Strength by Accelerating and Shortening Fiber Development

Low potassium (K)-induced premature senescence in cotton has been observed worldwide, but how it affects cotton fiber properties remain unclear. We hypothesized that K deficiency affects cotton fiber properties by causing disordered fiber development, which may in turn be caused by the induction of a carbohydrate acquisition difficulty. To investigate this issue, we employed a low-K-sensitive cotton cultivar Siza 3 and a low-K-tolerant cultivar Simian 3 and planted them in three regions of different K supply. Data concerning lint yield, Pn and main fiber properties were collected from three years of testing. Soil K deficiency significantly accelerated fiber cellulose accumulation and dehydration processes, which, together with previous findings, suggests that the low-K induced carbohydrate acquisition difficulty could cause disordered fiber development by stimulating the expression of functional proteins such as CDKA (cyclin-dependent kinase). As a result, fiber strength and lint weight were reduced by up to 7.8% and 2.1%, respectively. Additional quantitative analysis revealed that the degree of accelerated fiber development negatively correlated with fiber strength. According to the results of this study, it is feasible to address the effects of soil K deficiency on fiber properties using existing cultivation strategies to prevent premature senescence of cotton plants.

Comparative assessment of genetic diversity in cytoplasmic and nuclear genome of upland cotton

The importance of the cytoplasmic genome for many economically important traits is well documented in several crop species, including cotton. There is no report on application of cotton chloroplast specific SSR markers as a diagnostic tool to study genetic diversity among improved Upland cotton lines. The complete plastome sequence information in GenBank provided us an opportunity to report on 17 chloroplast specific SSR markers using a cost-effective data mining strategy. Here we report the comparative analysis of genetic diversity among a set of 42 improved Upland cotton lines using SSR markers specific to chloroplast and nuclear genome, respectively. Our results revealed that low to moderate level of genetic diversity existed in both nuclear and cytoplasm genome among this set of cotton lines. However, the specific estimation suggested that genetic diversity is lower in cytoplasmic genome compared to the nuclear genome among this set of Upland cotton lines. In summary, this research is important from several perspectives. We detected a set of cytoplasm genome specific SSR primer pairs by using a cost-effective data mining strategy. We reported for the first time the genetic diversity in the cytoplasmic genome within a set of improved Upland cotton accessions. Results revealed that the genetic diversity in cytoplasmic genome is narrow, compared to the nuclear genome within this set of Upland cotton accessions. Our results suggested that most of these polymorphic chloroplast SSRs would be a valuable complementary tool in addition to the nuclear SSR in the study of evolution, gene flow and genetic diversity in Upland cotton.

Construction of a high-density genetic map by specific locus amplified fragment sequencing (SLAF-seq) and its application to Quantitative Trait Loci (QTL) analysis for boll weight in upland cotton (Gossypium hirsutum.)

BACKGROUND

Upland Cotton (Gossypium hirsutum) is one of the most important worldwide crops it provides natural high-quality fiber for the industrial production and everyday use. Next-generation sequencing is a powerful method to identify single nucleotide polymorphism markers on a large scale for the construction of a high-density genetic map for quantitative trait loci mapping.

RESULTS

In this research, a recombinant inbred lines population developed from two upland cotton cultivars 0-153 and sGK9708 was used to construct a high-density genetic map through the specific locus amplified fragment sequencing method. The high-density genetic map harbored 5521 single nucleotide polymorphism markers which covered a total distance of 3259.37 cM with an average marker interval of 0.78 cM without gaps larger than 10 cM. In total 18 quantitative trait loci of boll weight were identified as stable quantitative trait loci and were detected in at least three out of 11 environments and explained 4.15-16.70 % of the observed phenotypic variation. In total, 344 candidate genes were identified within the confidence intervals of these stable quantitative trait loci based on the cotton genome sequence. These genes were categorized based on their function through gene ontology analysis, Kyoto Encyclopedia of Genes and Genomes analysis and eukaryotic orthologous groups analysis.

CONCLUSIONS

This research reported the first high-density genetic map for Upland Cotton (Gossypium hirsutum) with a recombinant inbred line population using single nucleotide polymorphism markers developed by specific locus amplified fragment sequencing. We also identified quantitative trait loci of boll weight across 11 environments and identified candidate genes within the quantitative trait loci confidence intervals. The results of this research would provide useful information for the next-step work including fine mapping, gene functional analysis, pyramiding breeding of functional genes as well as marker-assisted selection.

Breeding Potential of Introgression Lines Developed from Interspecific Crossing between Upland Cotton (Gossypium hirsutum) and Gossypium barbadense: Heterosis, Combining Ability and Genetic Effects

Upland cotton (Gossypium hirstum L.), which produces more than 95% of the world natural cotton fibers, has a narrow genetic base which hinders progress in cotton breeding. Introducing germplasm from exotic sources especially from another cultivated tetraploid G. barbadense L. can broaden the genetic base of Upland cotton. However, the breeding potential of introgression lines (ILs) in Upland cotton with G. barbadense germplasm integration has not been well addressed. This study involved six ILs developed from an interspecific crossing and backcrossing between Upland cotton and G. barbadense and represented one of the first studies to investigate breeding potentials of a set of ILs using a full diallel analysis. High mid-parent heterosis was detected in several hybrids between ILs and a commercial cultivar, which also out-yielded the high-yielding cultivar parent in F₁, F₂ and F₃ generations. A further analysis indicated that general ability (GCA) variance was predominant for all the traits, while specific combining ability (SCA) variance was either non-existent or much lower than GCA. The estimated GCA effects and predicted additive effects for parents in each trait were positively correlated (at P<0.01). Furthermore, GCA and additive effects for each trait were also positively correlated among generations (at P<0.05), suggesting that F₂ and F₃ generations can be used as a proxy to F₁ in analyzing combining abilities and estimating genetic parameters. In addition, differences between reciprocal crosses in F₁ and F₂ were not significant for yield, yield components and fiber quality traits. But maternal effects appeared to be present for seed oil and protein contents in F₃. This study identified introgression lines as good general combiners for yield and fiber quality improvement and hybrids with high heterotic vigor in yield, and therefore provided useful information for further utilization of introgression lines in cotton breeding.

Pollen-mediated gene flow from transgenic cotton under greenhouse conditions is dependent on different pollinators

With the large-scale release of genetically modified (GM) crops, there are ecological concerns on transgene movement from GM crops to non-GM counterparts and wild relatives. In this research, we conducted greenhouse experiments to measure pollen-mediated gene flow (PGF) in the absence and presence of pollinators (Bombus ignitus, Apis mellifera and Pieris rapae) in one GM cotton (resistant to the insect Helicoverpa armigera and the herbicide glyphosate) and two non-GM lines (Shiyuan321 and Hai7124) during 2012 and 2013. Our results revealed that: (1) PGF varied depending on the pollinator species, and was highest with B. ignitus (10.83%) and lowest with P. rapae (2.71%); (2) PGF with B. ignitus depended on the distance between GM and non-GM cottons; (3) total PGF to Shiyuan321 (8.61%) was higher than to Hai7124 (4.10%). To confirm gene flow, we tested hybrids carrying transgenes for their resistance to glyphosate and H. armigera, and most hybrids showed strong resistance to the herbicide and insect. Our research confirmed that PGF depended on pollinator species, distance between plants and the receptor plant.

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.

Mapping-by-sequencing of Ligon-lintless-1 (Li 1 ) reveals a cluster of neighboring genes with correlated expression in developing fibers of Upland cotton (Gossypium hirsutum L.)

Mapping-by-sequencing and SNP marker analysis were used to fine map the Ligon-lintless-1 ( Li 1 ) short fiber mutation in tetraploid cotton to a 255-kb region that contains 16 annotated proteins. The Ligon-lintless-1 (Li 1 ) mutant of cotton (Gossypium hirsutum L.) has been studied as a model for cotton fiber development since its identification in 1929; however, the causative mutation has not been identified yet. Here we report the fine genetic mapping of the mutation to a 255-kb region that contains only 16 annotated genes in the reference Gossypium raimondii genome. We took advantage of the incompletely dominant dwarf vegetative phenotype to identify 100 mutants (Li 1 /Li 1 ) and 100 wild-type (li 1 /li 1 ) homozygotes from a mapping population of 2567 F2 plants, which we bulked and deep sequenced. Since only homozygotes were sequenced, we were able to use a high stringency in SNP calling to rapidly narrow down the region harboring the Li 1 locus, and designed subgenome-specific SNP markers to test the population. We characterized the expression of all sixteen genes in the region by RNA sequencing of elongating fibers and by RT-qPCR at seven time points spanning fiber development. One of the most highly expressed genes found in this interval in wild-type fiber cells is 40-fold under-expressed at the day of anthesis (DOA) in the mutant fiber cells.  This gene is a major facilitator superfamily protein, part of the large family of proteins that includes auxin and sugar transporters. Interestingly, nearly all genes in this region were most highly expressed at DOA and showed a high degree of co-expression. Further characterization is required to determine if transport of hormones or carbohydrates is involved in both the dwarf and lintless phenotypes of Li 1 plants.

Development of Gossypium anomalum-derived microsatellite markers and their use for genome-wide identification of recombination between the G. anomalum and G. hirsutum genomes

We reported the first development of Gossypium anomalum -derived microsatellite markers and identification of recombination between sexually incompatible species by a synthesized hexaploid on genome level. To continue to develop improved cotton varieties, it is essential to transfer desired characters from diploid wild cotton species such as Gossypium anomalum to cultivated allotetraploid cotton species. However, interspecific reproductive barriers limit gene transfer between species. In a previous study, we used colchicine treatment to produce a synthesized hexaploid derived from an interspecific hybrid between Gossypium hirsutum and G. anomalum and demonstrated its hybridity and doubled status using morphological, cytological and molecular marker methods. In the current study, to effectively monitor G. anomalum genome components in the G. hirsutum background, we developed 5974 non-redundant G. anomalum-derived SSR primer pairs using RNA-Seq technology, which were combined with a publicly available physical map. Based on this combined map and segregation data from the BC2F1 population, we identified a set of 230 informative G. anomalum-specific SSR markers distributed on the chromosomes, which cover 95.72% of the cotton genome. After analyzing BC2F1 segregation data, 50 recombination types from 357 recombination events were identified, which cover 81.48% of the corresponding G. anomalum genome. A total of 203 recombination events occurred on chromosome 11, accounting for 56.86% of the recombination events on all chromosomes. Recombination hotspots were observed at marker intervals JAAS1148-NAU5100 on chromosome 1 and JAAS0426-NAU998 on chromosome 2. Therefore, all G. anomalum chromosomes are capable of recombining with At chromosomes in G. hirsutum. This study represents an important step towards introgressing desirable traits into cultivated cotton from the wild cotton species G. anomalum.

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.

Molecular characterization of the Gossypium Diversity Reference Set of the US National Cotton Germplasm Collection

A core marker set containing markers developed to be informative within a single commercial cotton species can elucidate diversity structure within a multi-species subset of the Gossypium germplasm collection. An understanding of the genetic diversity of cotton (Gossypium spp.) as represented in the US National Cotton Germplasm Collection is essential to develop strategies for collecting, conserving, and utilizing these germplasm resources. The US collection is one of the largest world collections and includes not only accessions with improved yield and fiber quality within cultivated species, but also accessions possessing sources of abiotic and biotic stress resistance often found in wild species. We evaluated the genetic diversity of a subset of 272 diploid and 1,984 tetraploid accessions in the collection (designated the Gossypium Diversity Reference Set) using a core set of 105 microsatellite markers. Utility of the core set of markers in differentiating intra-genome variation was much greater in commercial tetraploid genomes (99.7 % polymorphic bands) than in wild diploid genomes (72.7 % polymorphic bands), and may have been influenced by pre-selection of markers for effectiveness in the commercial species. Principal coordinate analyses revealed that the marker set differentiated interspecific variation among tetraploid species, but was only capable of partially differentiating among species and genomes of the wild diploids. Putative species-specific marker bands in G. hirsutum (73) and G. barbadense (81) were identified that could be used for qualitative identification of misclassifications, redundancies, and introgression within commercial tetraploid species. The results of this broad-scale molecular characterization are essential to the management and conservation of the collection and provide insight and guidance in the use of the collection by the cotton research community in their cotton improvement efforts.

Ancient gene duplicates in Gossypium (cotton) exhibit near-complete expression divergence

Whole genome duplication (WGD) is widespread in flowering plants and is a driving force in angiosperm diversification. The redundancy introduced by WGD allows the evolution of novel gene interactions and functions, although the patterns and processes of diversification are poorly understood. We identified ∼ 2,000 pairs of paralogous genes in Gossypium raimondii (cotton) resulting from an approximately 60 My old 5- to 6-fold ploidy increase. Gene expression analyses revealed that, in G. raimondii, 99.4% of the gene pairs exhibit differential expression in at least one of the three tissues (petal, leaf, and seed), with 93% to 94% exhibiting differential expression on a per-tissue basis. For 1,666 (85%) pairs, differential expression was observed in all tissues. These observations were mirrored in a time series of G. raimondii seed, and separately in leaf, petal, and seed of G. arboreum, indicating expression level diversification before species divergence. A generalized linear model revealed 92.4% of the paralog pairs exhibited expression divergence, with most exhibiting significant gene and tissue interactions indicating complementary expression patterns in different tissues. These data indicate massive, near-complete expression level neo- and/or subfunctionalization among ancient gene duplicates, suggesting these processes are essential in their maintenance over ∼ 60 Ma.

Genetic diversity and population structure in the US Upland cotton (Gossypium hirsutum L.)

Genetic diversity and population structure in the US Upland cotton was established and core sets of allelic richness were identified for developing association mapping populations in cotton. Elite plant breeding programs could likely benefit from the unexploited standing genetic variation of obsolete cultivars without the yield drag typically associated with wild accessions. A set of 381 accessions comprising 378 Upland (Gossypium hirsutum L.) and 3 G. barbadense L. accessions of the United States cotton belt were genotyped using 120 genome-wide SSR markers to establish the genetic diversity and population structure in tetraploid cotton. These accessions represent more than 100 years of Upland cotton breeding in the United States. Genetic diversity analysis identified a total of 546 alleles across 141 marker loci. Twenty-two percent of the alleles in Upland accessions were unique, specific to a single accession. Population structure analysis revealed extensive admixture and identified five subgroups corresponding to Southeastern, Midsouth, Southwest, and Western zones of cotton growing areas in the United States, with the three accessions of G. barbadense forming a separate cluster. Phylogenetic analysis supported the subgroups identified by STRUCTURE. Average genetic distance between G. hirsutum accessions was 0.195 indicating low levels of genetic diversity in Upland cotton germplasm pool. The results from both population structure and phylogenetic analysis were in agreement with pedigree information, although there were a few exceptions. Further, core sets of different sizes representing different levels of allelic richness in Upland cotton were identified. Establishment of genetic diversity, population structure, and identification of core sets from this study could be useful for genetic and genomic analysis and systematic utilization of the standing genetic variation in Upland cotton.

A cotton LIM domain-containing protein (GhWLIM5) is involved in bundling actin filaments

LIM-domain proteins play important roles in cellular processes in eukaryotes. In this study, a LIM protein gene, GhWLIM5, was identified in cotton. Quantitative RT-PCR analysis showed that GhWLIM5 was expressed widely in different cotton tissues and had a peak in expression during fiber elongation. GFP fluorescence assay revealed that cotton cells expressing GhWLIM5:eGFP fusion gene displayed a network distribution of eGFP fluorescence, suggesting that GhWLIM5 protein is mainly localized to the cell cytoskeleton. When GhWLIM5:eGFP transformed cells were stained with rhodamine-phalloidin there was consistent overlap in eGFP and rhodamine-palloidin signals, demonstrating that GhWLIM5 protein is colocalized with the F-actin cytoskeleton. In addition, high-speed cosedimentation assay verified that GhWLIM5 directly bound actin filaments, while low cosedimentation assay and microscopic observation indicated that GhWLIM5 bundled F-actin in vitro. Increasing amounts of GhWLIM5 protein were able to protect F-actin from depolymerization in vitro in the presence of Lat B (an F-actin depolymerizer). Our results contribute to a better understanding of the biochemical role of GhWLIM5 in modulating the dynamic F-actin network in cotton.

Surface plasmon resonance detection of transgenic Cry1Ac cotton ( Gossypium spp.)

The detection and identification of genetically modified (GM) plants are challenging issues that have arisen from the potential negative impacts of extensive cultivation of transgenic plants. The screening process is a long-term focus and needs specific detection strategies. Surface plasmon resonance (SPR) has been used to detect a variety of biomolecules including proteins and nucleic acids due to its ability to monitor specific intermolecular interactions. In the present study, two high-throughput, label-free, and specific methods based on SPR technology were developed to detect transgenic Cry1Ac cotton ( Gossypium spp.) by separately targeting protein and DNA. In the protein-based detection system, monoclonal anti-Cry1Ac antibodies were immobilized on the surface of a CM5 sensor chip. Conventional cotton samples were used to define the detection threshold. Transgenic cotton was easily identified within 5 min per sample. For the DNA-based model, a 25-mer biotinylated oligonucleotide probe was immobilized on an SA sensor chip. PCR products of Cry1Ac (230 bp) were used to investigate the reaction conditions. The sensitivity of the constructed sensor chip was identified at concentrations as low as 0.1 nM based on its complementary base pairing.

[Ecological regionalization of cotton varieties based on GGE biplot]

By using the heritability-adjusted GGE biplot analysis method, and taking the trial sites Anqing, Nanyang, Huanggang, Jingzhou, Wuhan, Xiangyang, Changde, Yueyang, Nanjing, Nantong, Yancheng, Jiujiang, Jianyang, Shehong, and Cixi as the representative cotton-planting areas in the Yangtze River basin, the ecological regionalization of cotton varieties in the basin was made based on the lint cotton yield, and the regionalization results were adjusted by the information ratio (IR) method, aimed to provide scientific basis for the selection of cotton varieties in the cotton-planting areas of the basin. The cotton-planting areas in the Yangtze River basin could be divided into three ecological regions, i.e., the "Sichuan basin cotton region" with Jianyang and Shehong as the representative, the "Nan-Xiang basin cotton region" with Xiangyang and Nanyang as the representative, and the "majority complex cotton region in the Yangtze River basin" including all the other sites in the basin.

Analysis of complete nucleotide sequences of 12 Gossypium chloroplast genomes: origin and evolution of allotetraploids

BACKGROUND

Cotton (Gossypium spp.) is a model system for the analysis of polyploidization. Although ascertaining the donor species of allotetraploid cotton has been intensively studied, sequence comparison of Gossypium chloroplast genomes is still of interest to understand the mechanisms underlining the evolution of Gossypium allotetraploids, while it is generally accepted that the parents were A- and D-genome containing species. Here we performed a comparative analysis of 13 Gossypium chloroplast genomes, twelve of which are presented here for the first time.

METHODOLOGY/PRINCIPAL FINDINGS

The size of 12 chloroplast genomes under study varied from 159,959 bp to 160,433 bp. The chromosomes were highly similar having >98% sequence identity. They encoded the same set of 112 unique genes which occurred in a uniform order with only slightly different boundary junctions. Divergence due to indels as well as substitutions was examined separately for genome, coding and noncoding sequences. The genome divergence was estimated as 0.374% to 0.583% between allotetraploid species and A-genome, and 0.159% to 0.454% within allotetraploids. Forty protein-coding genes were completely identical at the protein level, and 20 intergenic sequences were completely conserved. The 9 allotetraploids shared 5 insertions and 9 deletions in whole genome, and 7-bp substitutions in protein-coding genes. The phylogenetic tree confirmed a close relationship between allotetraploids and the ancestor of A-genome, and the allotetraploids were divided into four separate groups. Progenitor allotetraploid cotton originated 0.43-0.68 million years ago (MYA).

CONCLUSION

Despite high degree of conservation between the Gossypium chloroplast genomes, sequence variations among species could still be detected. Gossypium chloroplast genomes preferred for 5-bp indels and 1-3-bp indels are mainly attributed to the SSR polymorphisms. This study supports that the common ancestor of diploid A-genome species in Gossypium is the maternal source of extant allotetraploid species and allotetraploids have a monophyletic origin. G. hirsutum AD1 lineages have experienced more sequence variations than other allotetraploids in intergenic regions. The available complete nucleotide sequences of 12 Gossypium chloroplast genomes should facilitate studies to uncover the molecular mechanisms of compartmental co-evolution and speciation of Gossypium allotetraploids.

[Diurnal variation of gas exchange and chlorophyll fluorescence parameters of cotton functional leaves under effects of soil salinity]

A two-year (2007-2008) pot experiment with cotton varieties Sumian 12 (salinity-sensitive) and Zhongmiansuo 44 (salinity-tolerance) was conducted at the Pailou experimental station of Nanjing Agricultural University to study the diurnal variation of the gas exchange and chlorophyll fluorescence parameters of cotton functional leaves under five levels (0, 0.35%, 0.60% , 0.85%, and 1.00%) of soil salinity. With the increase of soil salinity, the concentrations of Na+, Cl-, and Mg2+ in functional leaves increased, whereas the concentrations of K+ and Ca2+ decreased. The salinity level <0. 35% had little effects on the gas exchange and chlorophyll fluorescence parameters, but that >0.35% depressed the net photosynthetic rate (Pn) dramatically. At the salinity level >0.35%, the sensitivity of functional leaves to daytime photon flux density (PFD) and air temperature (Ta) enhanced, which in turn resulted in more severe photo- and temperature inhibition, and changed the diurnal variation patterns of Pn and stomatal conductance (Gs) from a one-peak curve to a constantly decreasing one. Along with the variations of daytime PED and Ta, the diurnal variation patterns of the maximum photochemical efficiency (F(v)/F(m)), quantum yield of electron transport (phi(PS II), and photochemical quenching coefficient (q(P)) of functional leaves presented a V-shaped curve, with the minimum value appeared at 12:00-13:00, while the non-photochemical quenching coefficient (q(N)) showed a single-peak curve. Soil salinity decreased the F(v)/F(m), phi(PS II), and q(P) significantly, but increased the q(N) and enlarged its change trend. The comparatively low concentrations of Na+ and Cl- and the relatively high concentrations of K+ and Ca2+ in salt-tolerant Zhongmiansuo 44 functional leaves benefited the relative stability of PS II, and the maintenance of a relatively high thermal dissipation capacity could be one of the reasons for a high level of Pn at high salinity level.

An L1 box binding protein, GbML1, interacts with GbMYB25 to control cotton fibre development

Transcription factors play key roles in plant development through their interaction with cis-elements and/or other transcription factors. A HD-Zip IV family transcription factor, Gossypium barbadense Meristem Layer 1 (GbML1) has been identified and characterized here. GbML1 specifically bound to the L1 box and the promoters of GbML1 and GbRDL1. GbML1 physically interacted with a key regulator of cotton fibre development, GbMYB25. Truncated and point mutation assays indicated the START-SAD domain was required for the binding to the C terminal domain (CTD) of GbMYB25. GbML1 overexpression in Arabidopsis increased the number of trichomes on stems and leaves and increased the accumulation of anthocyanin in leaves. Taken together, the L1 box binding protein, GbML1 was identified as the first partner for GbMYB25 and the role of START domain was discovered to be a protein binding domain in plants. Our findings will help the improvement of cotton fibre production and the understanding of the key role of HD-Zip family and MYB family in plants.

Identification of the family of aquaporin genes and their expression in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Cotton (Gossypium spp.) is produced in over 30 countries and represents the most important natural fiber in the world. One of the primary factors affecting both the quantity and quality of cotton production is water. A major facilitator of water movement through cell membranes of cotton and other plants are the aquaporin proteins. Aquaporin proteins are present as diverse forms in plants, where they function as transport systems for water and other small molecules. The plant aquaporins belong to the large major intrinsic protein (MIP) family. In higher plants, they consist of five subfamilies including plasma membrane intrinsic proteins (PIP), tonoplast intrinsic proteins (TIP), NOD26-like intrinsic proteins (NIP), small basic intrinsic proteins (SIP), and the recently discovered X intrinsic proteins (XIP). Although a great deal is known about aquaporins in plants, very little is known in cotton.

RESULTS

From a molecular cloning effort, together with a bioinformatic homology search, 71 upland cotton (G. hirsutum) aquaporin genes were identified. The cotton aquaporins consist of 28 PIP and 23 TIP members with high sequence similarity. We also identified 12 NIP and 7 SIP members that showed more divergence. In addition, one XIP member was identified that formed a distinct 5th subfamily. To explore the physiological roles of these aquaporin genes in cotton, expression analyses were performed for a select set of aquaporin genes from each subfamily using semi-quantitative reverse transcription (RT)-PCR. Our results suggest that many cotton aquaporin genes have high sequence similarity and diverse roles as evidenced by analysis of sequences and their expression.

CONCLUSION

This study presents a comprehensive identification of 71 cotton aquaporin genes. Phylogenetic analysis of amino acid sequences divided the large and highly similar multi-gene family into the known 5 aquaporin subfamilies. Together with expression and bioinformatic analyses, our results support the idea that the genes identified in this study represent an important genetic resource providing potential targets to modify the water use properties of cotton.

Rapid DNA loss as a counterbalance to genome expansion through retrotransposon proliferation in plants

Transposable elements, particularly LTR-retrotransposons, comprise the primary vehicle for genome size expansion in plants, while DNA removal through illegitimate recombination and intrastrand homologous recombination serve as the most important counteracting forces to plant genomic obesity. Despite extensive research, the relative impact of these opposing forces and hence the directionality of genome size change remains unknown. In Gossypium (cotton), the 3-fold genome size variation among diploids is due largely to copy number variation of the gypsy-like retrotransposon Gorge3. Here we combine comparative sequence analysis with a modeling approach to study the directionality of genome size change in Gossypium. We demonstrate that the rate of DNA removal in the smaller genomes is sufficient to reverse genome expansion through Gorge3 proliferation. These data indicate that rates of DNA loss can be highly variable even within a single plant genus, and that the known mechanisms of DNA loss can indeed reverse the march toward genomic obesity.

[Low-temperature responses of enzyme activities related to fiber development of two cotton (Gossvpium hirsutum L.) cultivars with different temperature-sensitivity]

Taking two cotton cultivars with difterent temperature-sensitivity during tneir liner strength formation as test materials, a field experiment of different sowing dates was conducted in Nanjing of Jiangsu Province in 2006 and 2007 to study the effects of low temperature on the activities and gene expression of the enzymes related to fiber development. The low temperature induced by late sowing (with the mean daily minimum temperature being 21.1, 20.5, and 18.1 degrees C during fiber development period) had definite effects on the enzyme activities, and accordingly, the fiber strength formation. Low temperature increased the invertase and beta-1, 3-glucanase activities, decreased the sucrose synthase and sucrose phosphate synthase activities, prolonged the time with higher gene expression level of Expansin and sucrose synthase, and delayed the expression peak and decreased the gene expression quantity of beta-1,3-glucanase. There existed significant differences in the low-temperature responses of related enzymes activities between the two cultivars, with the change ranges of the enzyme activities being larger for temperature-sensitive cultivar Sumian 15 than for temperature-insensitive cultivar Kemian 1, which could be the main reasons leading to the different temnerature-sensitivitv of the two cotton cultivars during their fiber strength formation.

[Cotton identification and extraction using near infrared sensor and object-oriented spectral segmentation technique]

The real-time, effective and reliable method of identifying crop is the foundation of scientific management for crop in the precision agriculture. It is also one of the key techniques for the precision agriculture. However, this expectation cannot be fulfilled by the traditional pixel-based information extraction method with respect to complicated image processing and accurate objective identification. In the present study, visible-near infrared image of cotton was acquired using high-resolution sensor. Object-oriented segmentation technique was performed on the image to produce image objects and spatial/spectral features of cotton. Afterwards, nearest neighbor classifier integrated the spectral, shape and topologic information of image objects to precisely identify cotton according to various features. Finally, 300 random samples and an error matrix were applied to undertake the accuracy assessment of identification. Although errors and confusion exist, this method shows satisfying results with an overall accuracy of 96.33% and a KAPPA coefficient of 0.926 7, which can meet the demand of automatic management and decision-making in precision agriculture.