Cotton Genome Research Articles

Response of cotton genotypes to water and heat stress: from field to genes

Cotton is a crop of tropical and subtropical regions but the seed cotton yield is highly influenced by abiotic stresses like drought and heat. Response of cotton genome to abiotic stresses is highly complex and involve many genes. A comprehensive study, involving cotton genotypes developed through conventional and synthetic tetraploid method, was designed to (i) study the introgression of heat and water stress tolerance by using wild relatives (ii) evaluate genetic markers for marker assisted selection against water and heat stress. Two separate experiments for water and heat stress tolerance with a common control were established. Treatments in each experiment include a control and a stress treatment. Heat stress was applied by sowing crop two month earlier than the control treatment, whereas water stress was imposed by withholding alternate irrigation. Analyses of variance depicted highly significant (P ≤ 0.01) effect of genotypes and both stresses on boll retention, boll weight and seed cotton yield. Interaction of genotypes with stress in both experiments was also highly significant (P ≤ 0.01). Genotypes derived from interspecific crosses performed consistently in stress conditions compared to control which prove it a reliable method to introgress stress related genes from wild parents. Four genes reported for water stress tolerance and five genes reported for heat stress tolerance were evaluated by field results for efficient marker assisted selection (MAS). Results verified drought stress genes but heat stress genes could not explain genetic variability caused by heat stress. It is concluded from the results that separate genes may be responsible for heat stress tolerance for vegetative and reproductive stages, therefore, selection criteria should include both the traits.

A comparative genomics approach revealed evolutionary dynamics of microsatellite imperfection and conservation in genus Gossypium

BackgroundOngoing molecular processes in a cell could target microsatellites, a kind of repetitive DNA, owing to length variations and motif imperfection. Mutational mechanisms underlying such kind of genetic variations have been extensively investigated in diverse organisms. However, obscure impact of ploidization, an evolutionary process of genome content duplication prevails mostly in plants, on non-coding DNA is poorly understood.ResultsGenome sequences of diversely originated plant species were examined for genome-wide motif imperfection pattern, and various analytical tools were employed to canvass characteristic relationships among repeat density, imperfection and length of microsatellites. Moreover, comparative genomics approach aided in exploration of microsatellites conservation footprints in Gossypium evolution. Based on our results, motif imperfection in repeat length was found intricately related to genomic abundance of imperfect microsatellites among 13 genomes. Microsatellite decay estimation depicted slower decay of long motif repeats which led to predominant abundance of 5-nt repeat motif in Gossypium species. Short motif repeats exhibited rapid decay through the evolution of Gossypium lineage ensuing drastic decrease of 2-nt repeats, of which, “AT” motif type dilapidated in cultivated tetraploids of cotton.ConclusionThe outcome could be a directive to explore comparative evolutionary footprints of simple non-coding genetic elements i.e., repeat elements, through the evolution of genus-specific characteristics in cotton genomes.

Microdissection of the Ah01 chromosome in upland cotton and microcloning of resistance gene anologs from the single chromosome

BackgroundChromosome microdissection is one of the most important techniques in molecular cytogenetic research. Cotton (Gossypium Linnaeus, 1753) is the main natural fiber crop in the world. The resistance gene analog (RGA) cloning after its single chromosome microdissection can greatly promote cotton genome research and breeding.ResultsUsing the linker adaptor PCR (LA-PCR) with the primers of rice disease-resistance homologues, three nucleotide sequences PS016 (KU051681), PS054 (KU051682), and PS157 (KU051680) were obtained from the chromosome Ah01 of upland cotton (cv. TM-1). The Blast results showed that the three sequences are the nucleotide binding site-leucine rich repeat (NBS-LRR) type RGAs. Clustering results indicated that they are homologous to these published RGAs. Thus, the three RGAs can definitely be confirmed as NBS-LRR class of RGAs in upland cotton.ConclusionsUsing single chromosome microdissection technique, DNA libraries containing cotton RGAs were obtained. This technique can promote cotton gene cloning, marker development and even the improvement of cotton genome research and breeding.

Nucleotide diversity in the two co-resident genomes of allopolyploid cotton

Genetic diversity within and among populations lies at the heart of evolution. Unraveling the extent to which each intrinsic or extrinsic factor determines levels of diversity among genes, populations, and species is challenging, given the difficulty of isolating any single potentially important variable from all others. Allopolyploid species provide an opportunity to disentangle external and intrinsic factors, as the two (or more) homoeologous genomes co-occur in the same nucleus, often exhibiting high collinearity along homoeologous chromosomes. Here we evaluate the pace of molecular evolution and intraspecific, intragenomic diversity in two species of allopolyploid Gossypium, G. hirsutum and G. barbadense, using several hundred genes sequenced from multiple accessions of each species. Genic diversity in both species is low, having been influenced both by the polyploid bottleneck and a domestication bottleneck (for cultivated accessions), but with a directional bias in homoeolog diversity favoring the same genome in both allopolyploids. Total diversity is remarkably similar for the two homoeologous genomes overall, but the two copies of many gene pairs have accumulated statistically different diversity levels, and in a biased fashion with respect to genome. Domesticated accessions show reduced diversity in both genomes, as expected, but with a much greater reduction in one of the two homoeologous genomes. Furthermore, this biased reduction affects opposite homoeologous genomes in the two species. Interspecific introgression has played a role in shaping diversity within each species. Introgression was only detected for certain accessions, and only from G. barbadense into G. hirsutum in one of the two co-resident genomes.

Quantitative Trait Locus Mapping for Verticillium wilt Resistance in an Upland Cotton Recombinant Inbred Line Using SNP-Based High Density Genetic Map.

Verticillium wilt (VW) caused by Verticillium dahlia Kleb is one of the most destructive diseases of cotton. Numerous efforts have been made to improve the resistance of upland cotton against VW, with little progress achieved due to the paucity of upland cotton breeding germplasms with high level of resistance to VW. Gossypium barbadense was regarded as more resistant compared to upland cotton; however, it is difficult to apply the resistance from G. barbadense to upland cotton improvement because of the hybrid breakdown and the difficulty to fix resistant phenotype in their interspecific filial. Here we reported QTLs related to VW resistance identified in upland cotton based on 1 year experiment in greenhouse with six replications and 4 years investigations in field with two replications each year. In total, 119 QTLs of disease index (DI) and of disease incidence (DInc) were identified on 25 chromosome of cotton genome except chromosome 13 (c13). For DI, 62 QTLs explaining 3.7–12.2% of the observed phenotypic variations were detected on 24 chromosomes except c11 and c13. For DInc, 59 QTLs explaining 2.3–21.30% of the observed PV were identified on 19 chromosomes except c5, c8, c12-c13, c18-c19, and c26. Seven DI QTLs were detected to be stable in at least environments, among which six have sGK9708 alleles, while 28 DInc QTLs were detected to be stable in at least environments. Eighteen QTL clusters containing 40 QTLs were identified on 13 chromosomes (c1-c4, c6-c7, c10, c14, c17 c20-c22, and c24-c25). Most of the stable QTLs aggregated into these clusters. These QTLs and clusters identification can be an important step toward Verticillium wilt resistant gene cloning in upland cotton and provide useful information to understand the complex genetic bases of Verticillium wilt resistance.

The crosstalk between Target of Rapamycin (TOR) and Jasmonic Acid (JA) signaling existing in Arabidopsis and cotton

Target of rapamycin (TOR) acts as an important regulator of cell growth, development and stress responses in most examined diploid eukaryotes. However, little is known about TOR in tetraploid species such as cotton. Here, we show that TORC1-S6K-RPS6, the major signaling components, are conserved and further expanded in cotton genome. Though the cotton seedlings are insensitive to rapamycin, AZD8055, the second-generation inhibitor of TOR, can significantly suppress the growth in cotton. Global transcriptome analysis revealed that genes associated with jasmonic acid (JA) biosynthesis and transduction were significantly altered in AZD8055 treated cotton seedlings, suggesting the potential crosstalk between TOR and JA signaling. Pharmacological and genetic approaches have been employed to get further insights into the molecular mechanism of the crosstalk between TOR and JA. Combination of AZD8055 with methyl jasmonate can synergistically inhibit cotton growth, and additionally JA levels were significantly increased when cotton seedlings were subjected to AZD8055. JA biosynthetic and signaling mutants including jar1, coi1-2 and myc2-2 displayed TOR inhibitor-resistant phenotypes, whereas COI1 overexpression transgenic lines and jaz10 exhibited sensitivity to AZD8055. Consistently, cotton JAZ can partially rescue TOR-suppressed phenotypes in Arabidopsis. These evidences revealed that the crosstalk between TOR and JA pathway operates in cotton and Arabidopsis.

Uncovering SNP and indel variations of tetraploid cottons by SLAF-seq

BackgroundCotton (Gossypium spp.), as the world’s most utilized textile fibre source, is an important, economically valuable crop worldwide. Understanding the genomic variation of tetraploid cotton species is important for exploitation of the excellent characteristics of wild cotton and for improving the diversity of cotton in breeding. However, the discovery of DNA polymorphisms in tetraploid cotton genomes has lagged behind other important crops.ResultsA total of 111,795,823 reads, 467,735 specific length amplified fragment (SLAF) tags and 139,176 high-quality DNA polymorphisms were identified using specific length amplified fragment sequencing (SLAF-seq), including 132,880 SNPs and 6,296 InDels between the reference genome (TM-1) and the five tetraploid cotton species. Intriguingly, gene ontology (GO) enrichment analysis revealed that a number of significant terms were related to reproduction in G. barbadense acc. 3–79. Based on the new data sets, we reconstructed phylogenetic trees that showed a high concordance to the phylogeny of diploid and polyploid cottons. A large amount of interspecific genetic variations were identified, and some of them were validated by the single-strand conformation polymorphism (SSCP) method, which will be applied in introgression genetics and breeding with G. hirsutum cv. Emian22 as the receptor and the other species as donors.ConclusionsUsing SLAF-seq, a large number of DNA polymorphisms were identified. The comprehensive analysis of DNA polymorphisms provided invaluable insights into the different tetraploid cotton species. More importantly, the identification of numerous interspecific genetic variations provides the basis and is very practical for future introgression breeding. The results presented herein provide a valuable genomic resource for new insights into the genetics and breeding of cotton.

Targeted mutagenesis in cotton (Gossypium hirsutum L.) using the CRISPR/Cas9 system

The CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)/Cas9 system has been widely used for genome editing in various plants because of its simplicity, high efficiency and design flexibility. However, to our knowledge, there is no report on the application of CRISPR/Cas9-mediated targeted mutagenesis in cotton. Here, we report the genome editing and targeted mutagenesis in upland cotton (Gossypium hirsutum L., hereafter cotton) using the CRISPR/Cas9 system. We designed two guide RNAs to target distinct sites of the cotton Cloroplastos alterados 1 (GhCLA1) and vacuolar H+-pyrophosphatase (GhVP) genes. Mutations in these two genes were detected in cotton protoplasts. Most of the mutations were nucleotide substitutions, with one nucleotide insertion and one substitution found in GhCLA1 and one deletion found in GhVP in cotton protoplasts. Subsequently, the two vectors were transformed into cotton shoot apexes through Agrobacterium-mediated transformation, resulting in efficient target gene editing. Most of the mutations were nucleotide deletions, and the mutation efficiencies were 47.6–81.8% in transgenic cotton plants. Evaluation using restriction-enzyme-PCR assay and sequence analysis detected no off-target mutations. Our results indicated that the CRISPR/Cas9 system was an efficient and specific tool for targeted mutagenesis of the cotton genome.

CRISPR/Cas9-mediated targeted mutagenesis in upland cotton (Gossypium hirsutum L.).

The clustered, regularly interspaced, short palindromic repeats (CRISPR)/CRISPR associated (Cas)9 protein system has emerged as a simple and efficient tool for genome editing in eukaryotic cells. It has been shown to be functional in several crop species, yet there are no reports on the application of this or any other genome editing technologies in the cotton plant. Cotton is an important crop that is grown mainly for its fiber, but its seed also serves as a useful source of edible oil and feed protein. Most of the commercially-grown cotton is tetraploid, thus making it much more difficult to target both sets of homeologous alleles. Therefore, in order to understand the efficacy of the CRISPR/Cas9 system to target a gene within the genome of cotton, we made use of a transgenic cotton line previously generated in our laboratory that had a single copy of the green fluorescent protein (GFP) gene integrated into its genome. We demonstrate, for the first time, the use of this powerful new tool in targeted knockout of a gene residing in the cotton genome. By following the loss of GFP fluorescence, we were able to observe the cells that had undergone targeted mutations as a result of CRISPR/Cas9 activity. In addition, we provide examples of the different types of indels obtained by Cas9-mediated cleavage of the GFP gene, guided by three independent sgRNAs. The results provide useful information that will help us target important native genes in the cotton plant in future.

A high-efficiency CRISPR/Cas9 system for targeted mutagenesis in Cotton (Gossypium hirsutum L.)

The complex allotetraploid genome is one of major challenges in cotton for repressing gene expression. Developing site-specific DNA mutation is the long-term dream for cotton breeding scientists. The clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) system is emerging as a robust biotechnology for targeted-DNA mutation. In this study, two sgRNAs, GhMYB25-like-sgRNA1 and GhMYB25-like-sgRNA2, were designed in the identical genomic regions of GhMYB25-like A and GhMYB25-like D, which were encoded by cotton A subgenome and the D subgenome, respectively, was assembled to direct Cas9-mediated allotetraploid cotton genome editing. High proportion (14.2–21.4%) CRISPR/Cas9-induced specific truncation events, either from GhMYB25-like A DNA site or from GhMYB25-like D DNA site, were detected in 50% examined transgenic cotton through PCR amplification assay and sequencing analyses. Sequencing results also demonstrated that 100% and 98.8% mutation frequency were occurred on GhMYB25-like-sgRNA1 and GhMYB25-like-sgRNA2 target site respectively. The off-target effect was evaluated by sequencing two putative off-target sites, which have 3 and 1 mismatched nucleotides with GhMYB25-like-sgRNA1 and GhMYB25-like-sgRNA2, respectively; all the examined samples were not detected any off-target-caused mutation events. Thus, these results demonstrated that CRISPR/Cas9 is qualified for generating DNA level mutations on allotetraploid cotton genome with high-efficiency and high-specificity.

Targeted identification of association between cotton fiber quality traits and microsatellite markers

Primitive and exotic accessions of cotton are potential sources of favorable alleles for genetic improvement, enriching diversity in the genetically constricted gene pool of elite cultivars. Three exotic accessions of cotton (MDN101, MDN063 and MDN257), collected from different parts of Central America and converted to day-neutral flowering; and four elite cultivars (PD94042, DES56, PMHS200 and Acala Maxxa) representing the US cotton gene pool were used as parents to create experimental populations. The corresponding F2 and F2:3 progenies of these populations were grown in two successive years (i.e., some in 2011–2012, some in 2012–2013) and phenotypes were scored in both F2 and F2:3 progenies in all 3 years (2011–2012–2013). These populations were screened with 113 polymorphic microsatellite markers selected from “hotspots” for fiber quality quantitative trait loci in the cotton genome and single marker analyses were performed to identify significant associations of the markers with six fiber quality traits. A total of 134 nominal marker-trait associations were identified, among which 15 were significant after Bonferroni correction for multiple comparisons. In 67 of 134 nominal associations and 4 of 15 significant associations, the exotic parents contributed favorable alleles to multiple backgrounds and for multiple traits, in addition to the traits for which they were selected. These results indicate that utilization of exotic and wild accessions of cotton is useful in introducing favorable alleles into the cultivated cotton gene pool for genetic improvement.

Functional characterization of AGAMOUS-subfamily members from cotton during reproductive development and in response to plant hormones

Expression analysis of the AG -subfamily members from G. hirsutum during flower and fruit development. Reproductive development in cotton, including the fruit and fiber formation, is a complex process; it involves the coordinated action of gene expression regulators, and it is highly influenced by plant hormones. Several studies have reported the identification and expression of the transcription factor family MADS-box members in cotton ovules and fibers; however, their roles are still elusive during the reproductive development in cotton. In this study, we evaluated the expression profiles of five MADS-box genes (GhMADS3, GhMADS4, GhMADS5, GhMADS6 and GhMADS7) belonging to the AGAMOUS-subfamily in Gossypium hirsutum. Phylogenetic and protein sequence analyses were performed using diploid (G. arboreum, G. raimondii) and tetraploid (G. barbadense, G. hirsutum) cotton genomes, as well as the AG-subfamily members from Arabidopsis thaliana, Petunia hybrida and Antirrhinum majus. qPCR analysis showed that the AG-subfamily genes had high expression during flower and fruit development in G. hirsutum. In situ hybridization analysis also substantiates the involvement of AG-subfamily members on reproductive tissues of G. hirsutum, including ovule and ovary. The effect of plant hormones on AG-subfamily genes expression was verified in cotton fruits treated with gibberellin, auxin and brassinosteroid. All the genes were significantly regulated in response to auxin, whereas only GhMADS3, GhMADS4 and GhMADS7 genes were also regulated by brassinosteroid treatment. In addition, we have investigated the GhMADS3 and GhMADS4 overexpression effects in Arabidopsis plants. Interestingly, the transgenic plants from both cotton AG-like genes in Arabidopsis significantly altered the fruit size compared to the control plants. This alteration suggests that cotton AG-like genes might act regulating fruit formation. Our results demonstrate that members of the AG-subfamily in G. hirsutum present a conserved expression profile during flower development, but also demonstrate their expression during fruit development and in response to phytohormones.

Genome-wide identification and expression analysis of polyamine oxidase genes in upland cotton (Gossypium hirsutum L.)

Many reports showed that polyamine oxidase (PAO) genes play crucial roles in various abiotic and biotic stresses response. In this study, 12 putative PAO genes (GhPAO1–GhPAO12) were identified in upland cotton (Gossypium hirsutum L.) genome using the released genome database. A total of 349 cis-acting factors involved in hormone and stress response were predicted in 1.5-kb promoter regions at the upstream region of GhPAOs. These GhPAOs can be divided into four major groups when compared with their counterparts in Arabidopsis thaliana. Transcriptional level of GhPAOs were detected by qRT-PCR showing the highest levels in roots. Although GhPAOs exhibited various responses to NaCl treatment, the highest expression level was observed from 3 or 6 h after treatment. To further elucidate GhPAOs functions, GhPAO3 were selected for genetic transformation in Arabidopsis. The content of spermine was significantly decreased in transgenic Arabidopsis lines, whereas the content of putrescine was enhanced, which indicates a potential role of GhPAO3 in the conversion of spermidine and spermine. Discrepancy of salt resistance were observed in germination and seedling stage, which might be controlled by H2O2 concentration threshold. This study provides a valuable insight of the PAO genes in upland cotton response to abiotic stress.

Genomic, evolutionary and expression profile analysis of Hsp70 superfamily in A and D genome of cotton (Gossypium spp.) under the challenge of Verticillium dahliae

In this study, we comparatively analyzed the 115 Hsp70 genes identified in Gossypium raimondii, Gossypium hirsutum and Gossypium arboreum genomes. Those Hsp70 genes unequally distributed among chromosomes in A and D genome of cotton (Gossypium spp.), and were classified into 29 groups according to the homology of them. Based on the localization information of the orthologs in Arabidopsis, the Hsp70 proteins were predicted to locate in cytosol, endoplasmic reticulum, mitochondrion or chloroplast. Homologous analysis indicated the evolutionary conservation of Hsp70 in cotton. In addition, those Hsp70 genes were differently expressed in Suyuan-045, Hai-7124 and TM-1, which were highly resistant, resistant, and sensitive to Verticillium dahliae respectively. The expressions of 26 Hsp70 genes were induced by Verticillium dahliae except for Hsp70-07/16/25/26, and the result suggested the potential involvement of them in responding to Verticillium wilt. Hsp70-08/30/31 was highly expressed in both Suyuan-045 and Hai-7124, and it was hypothesized that they might be involved in the resistance to the invasion of Verticillium dahliae. 144h after inoculation with Verticillium dahliae, the expression of Hsp70-13/14/15 was only up-regulated in Suyuan-045, and it was assumed that they might be involved in resistance to the extension of Verticillium dahliae. Further study on those Hsp70 genes would be valuable to reveal the role of them in Verticillium wilt resistance.

Expansion and stress responses of the AP2/EREBP superfamily in cotton

BackgroundThe allotetraploid cotton originated from one hybridization event between an extant progenitor of Gosssypium herbaceum (A1) or G. arboreum (A2) and another progenitor, G. raimondii Ulbrich (D5) 1–1.5 million years ago (Mya). The APETALA2/ethylene-responsive element binding protein (AP2/EREBP) transcription factors constitute one of the largest and most conserved gene families in plants. They are characterized by their AP2 domain, which comprises 60–70 amino acids, and are classified into four main subfamilies: the APETALA2 (AP2), Related to ABI3/VP1 (RAV), Dehydration-Responsive Element Binding protein (DREB) and Ethylene-Responsive Factor (ERF) subfamilies. The AP2/EREBP genes play crucial roles in plant growth, development and biotic and abiotic stress responses. Hence, understanding the molecular characteristics of cotton stress tolerance and gene family expansion would undoubtedly facilitate cotton resistance breeding and evolution research.ResultsA total of 269 AP2/EREBP genes were identified in the G. raimondii (D5) cotton genome. The protein domain architecture and intron/exon structure are simple and relatively conserved within each subfamily. They are distributed throughout all chromosomes but are clustered on various chromosomes due to genomic tandem duplication. We identified 73 tandem duplicated genes and 221 segmental duplicated gene pairs which contributed to the expansion of AP2/EREBP superfamily. Of them, tandem duplication was the most important force of the expansion of the B3 group. Transcriptome analysis showed that 504 AP2/EREBP genes were expressed in at least one tested G. hirsutum TM-1 tissues. In G. hirsutum, 151 non-repeated genes of the DREB and ERF subfamily genes were responsive to different stresses: 132 genes were induced by cold, 63 genes by drought and 94 genes by heat. qRT-PCR confirmed that 13 GhDREB and 15 GhERF genes were induced by cold and/or drought. No transcripts detected for 53 of the 111 tandem duplicated genes in TM-1. In addition, some homoeologous genes showed biased expression toward either A-or D-subgenome.ConclusionsThe AP2/EREBP genes were obviously expanded in Gossypium. The GhDREB and GhERF genes play crucial roles in cotton stress responses. Our genome-wide analysis of AP2/EREBP genes in cotton provides valuable information for characterizing the molecular functions of AP2/EREBP genes and reveals insights into their evolution in polyploid plants.

Genome-wide analysis of salinity-stress induced DNA methylation alterations in cotton (Gossypium hirsutum L.) using the Me-DIP sequencing technology.

Cytosine DNA methylation is a significant form of DNA modification closely associated with gene expression in eukaryotes, fungi, animals, and plants. Although the reference genomes of cotton (Gossypium hirsutum L.) have been publically available, the salinity-stress-induced DNA methylome alterations in cotton are not well understood. Here, we constructed a map of genome-wide DNA methylation characteristics of cotton leaves under salt stress using the methylated DNA immunoprecipitation sequencing method. The results showed that the methylation reads on chromosome 9 were most comparable with those on the other chromosomes, but the greatest changes occurred on chromosome 8 under salt stress. The DNA methylation pattern analysis indicated that a relatively higher methylation density was found in the upstream2k and downstream2k elements of the CDS region and CG-islands. Almost 94% of the reads belonged to LTR-gspsy and LTR-copia, and the number of methylation reads in LTR-gypsy was four times greater than that in LTR-copia in both control and stressed samples. The analysis of differentially methylated regions (DMRs) showed that the gene elements upstream2k, intron, and downstream2k were hypomethylated, but the CDS regions were hypermethylated. The GO (Gene Ontology) analysis suggested that the methylated genes were most enriched in cellular processes, metabolic processes, cell parts and catalytic activities, which might be closely correlated with response to NaCl stress. In this study, we completed a genomic DNA methylation profile and conducted a DMR analysis under salt stress, which provided valuable information for the better understanding of epigenetics in response to salt stress in cotton.

<i>AtDREB2A-CA</i> Influences Root Architecture and Increases Drought Tolerance in Transgenic Cotton

Drought is a major environmental factor limiting cotton (Gossypium hirsutum L.) productivity worldwide and projected climate changes could increase their negative effects in the future. Thus, targeting the molecular mechanisms correlated with drought tolerance without reducing productivity is a challenge for plant breeding. In this way, we evaluated the effects of water deficit progress on AtDREB2A-CA transgenic cotton plant responses, driven by the stress-inducible rd29 promoter. Besides shoot and root morphometric traits, gas exchange and osmotic adjustment analyses were also included. Here, we present how altered root traits shown by transgenic plants impacted on physiological acclimation responses when submitted to severe water stress. The integration of AtDREB2A-CA into the cotton genome increased total root volume, surface area and total root length, without negatively affecting shoot morphometric growth parameters and nor phenotypic evaluated traits. Additionally, when compared to wild-type plants, transgenic plants (17-T0 plants and its progeny) highlighted a gradual pattern of phenotypic plasticity tosome photosynthetic parameters such as photosynthetic rate and stomatal conductance with water deficit progress. Transgene also promoted greater shoot development and root robustness (greater and deeper root mass) allowing roots to grow into deeper soil layers. The same morpho-physiological trend was observed in the subsequent generation (17.6-T2). Our results suggest that the altered root traits shown by transgenic plants are the major contributors to higher tolerance response, allowing the AtDRE2A-CA-cotton plants to maintain elevated stomatal conductance and assimilate rates and, consequently, reducing their metabolic costs involved in the antioxidant responses activation. These results also suggest that these morpho-physiological changes increased the number of reproductive structures retained per plant (26% higher) when compared with its non-transgenic counterpart. This is the first report of cotton plants overexpressing the AtDRE2A-CA transcription factor, demonstrating a morpho-physiological and yield advantages under drought stress, without displaying any yield penalty under irrigated conditions. The mechanisms by which the root traits influenced the acclimation of the transgenic plants to severe water deficit conditions are also discussed. These data present an opportunity to use this strategy in cotton breeding programs in order to improve drought adaptation toward better rooting features.

GrTEdb: the first web-based database of transposable elements in cotton (Gossypium raimondii).

Although several diploid and tetroploid Gossypium species genomes have been sequenced, the well annotated web-based transposable elements (TEs) database is lacking. To better understand the roles of TEs in structural, functional and evolutionary dynamics of the cotton genome, a comprehensive, specific, and user-friendly web-based database, Gossypium raimondii transposable elements database (GrTEdb), was constructed. A total of 14 332 TEs were structurally annotated and clearly categorized in G. raimondii genome, and these elements have been classified into seven distinct superfamilies based on the order of protein-coding domains, structures and/or sequence similarity, including 2929 Copia-like elements, 10 368 Gypsy-like elements, 299 L1, 12 Mutators, 435 PIF-Harbingers, 275 CACTAs and 14 Helitrons. Meanwhile, the web-based sequence browsing, searching, downloading and blast tool were implemented to help users easily and effectively to annotate the TEs or TE fragments in genomic sequences from G. raimondii and other closely related Gossypium species. GrTEdb provides resources and information related with TEs in G. raimondii, and will facilitate gene and genome analyses within or across Gossypium species, evaluating the impact of TEs on their host genomes, and investigating the potential interaction between TEs and protein-coding genes in Gossypium species. Database URL: http://www.grtedb.org/

Genetic dissection of lint yield and fiber quality traits of G. hirsutum in G. barbadense background

Gossypium hirsutum L. is a widely cultivated species characterized by its high yield and wide environmental adaptability, while Gossypium barbadense is well known for its superior fiber quality. In the present report, we, for the first time, developed G. hirsutum chromosome segment introgression lines (ILs) in a G. barbadense background (GhILs_Gb) and genetically dissected the inheritance of lint yield and fiber quality of G. hirsutum in G. barbadense background. The GhILs_Gb contains introgressed segments spanning 4121.20 cM, which represents 82.20% of the tetraploid cotton genome, with an average length of 18.65 cM. A total of 39 quantitative trait loci (QTLs) for six traits are identified in this IL population planted in Xinjiang. Four QTL clusters are detected. Of them, however, three clusters have deleterious effects on fiber length and strength and boll weight, and only one cluster on Chr. D9 can be used in marker-assisted selection (MAS) to increase lint percentage and decrease micronaire value in G. barbadense. QTL mapping showed that most of yield-related QTLs detected have positive effects and increase lint yield in G. barbadense, while most of fiber quality-related QTLs have deleterious effects except for micronaire. It suggested that G. hirsutum evolved to have a high lint yield. Several lines improved in lint percentage and boll size in G. barbadense by introgressed one fragment of G. hirsutum have been developed from the GhILs_Gb. The ILs developed, and the analyses presented here will enhance the understanding of the genetics of lint yield and fiber quality in G. hirsutum and facilitate further molecular breeding to improve lint yield in G. barbadense.

Genome-wide identification, phylogeny, and expression analysis of pectin methylesterases reveal their major role in cotton fiber development.

BackgroundPectin methylesterase (PME, EC 3.1.1.11) is a hydrolytic enzyme that utilizes pectin as substrates, and plays a significant role in regulating pectin reconstruction thereby regulating plant growth. Pectin is one of the important components of the plant cell wall, which forms the main structural material of cotton fiber. In this research, cotton genome information was used to identify PMEs.ResultsWe identified 80 (GaPME01-GaPME80) PME genes from diploid G. arboreum (A genome), 78 (GrPME01-GrPME78) PME genes from G. raimondii (D genome), and 135 (GhPME001-GhPME135) PME genes from tetraploid cotton G. hirsutum (AD genome). We further analyzed their gene structure, conserved domain, gene expression, and systematic evolution to lay the foundation for deeper research on the function of PMEs. Phylogenetic data indicated that members from the same species demonstrated relatively high sequence identities and genetic similarities. Analysis of gene structures showed that most of the PMEs genes had 2–3 exons, with a few having a variable number of exons from 4 to 6. There are nearly no differences in the gene structure of PMEs among the three (two diploid and one tetraploid) cotton species. Selective pressure analysis showed that the Ka/Ks value for each of the three cotton species PME families was less than one.ConclusionConserved domain analysis showed that PMEs members had a relatively conserved C-terminal pectinesterase domain (PME) while the N-terminus was less conserved. Moreover, some of the family members contained a pectin methylesterase inhibitor (PMEI) domain. The Ka/Ks ratios suggested that the duplicated PMEs underwent purifying selection after the duplication events. This study provided an important basis for further research on the functions of cotton PMEs. Results from qRT-PCR indicated that the expression level of different PMEs at various fiber developmental stages was different. Moreover, some of the PMEs showed fiber predominant expression in secondary wall thickening indicating tissue-specific expression patterns.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3365-z) contains supplementary material, which is available to authorized users.