Molecular Traits and Functional Analysis of the CLAVATA3/Endosperm Surrounding Region-Related Small Signaling Peptides in Three Species of Gossypium Genus

The CLAVATA3/endosperm surrounding region-related (CLE) small peptides are a group of C-terminally encoded and post-translationally modified signal molecules involved in regulating the growth and development of various plants. However, the function and evolution of these peptides have so far remained elusive in cotton. In this study, 55, 56, and 86 CLE genes were identified in the Gossypium raimondii, Gossypium arboreum, and Gossypium hirsutum genomes, respectively, and all members were divided into seven groups. These groups were distinctly different in their protein characteristics, gene structures, conserved motifs, and multiple sequence alignment. Whole genome or segmental duplications played a significant role in the expansion of the CLE family in cotton, and experienced purifying selection during the long evolutionary process in cotton. Cis-acting regulatory elements and transcript profiling revealed that the CLE genes of cotton exist in different tissues, developmental stages, and respond to abiotic stresses. Protein properties, structure prediction, protein interaction network prediction of GhCLE2, GhCLE33.2, and GhCLE28.1 peptides were, respectively, analyzed. In addition, the overexpression of GhCLE2, GhCLE33.2, or GhCLE28.1 in Arabidopsis, respectively, resulted in a distinctive shrub-like dwarf plant, slightly purple leaves, large rosettes with large malformed leaves, and lack of reproductive growth. This study provides important insights into the evolution of cotton CLEs and delineates the functional conservatism and divergence of CLE genes in the growth and development of cotton.

GrTCP11, a Cotton TCP Transcription Factor, Inhibits Root Hair Elongation by Down-Regulating Jasmonic Acid Pathway in Arabidopsis thaliana

TCP transcription factors play important roles in diverse aspects of plant development as transcriptional activators or repressors. However, the functional mechanisms of TCPs are not well understood, especially in cotton fibers. Here, we identified a total of 37 non-redundant TCP proteins from the diploid cotton (Gossypium raimondii), which showed great diversity in the expression profile. GrTCP11, an ortholog of AtTCP11, was preferentially expressed in cotton anthers and during fiber initiation and secondary cell wall synthesis stages. Overexpression of GrTCP11 in Arabidopsis thaliana reduced root hair length and delayed flowering. It was found that GrTCP11 negatively regulated genes involved in jasmonic acid (JA) biosynthesis and response, such as AtLOX4, AtAOS, AtAOC1, AtAOC3, AtJAZ1, AtJAZ2, AtMYC2, and AtERF1, which resulted in a decrease in JA concentration in the overexpressed transgenic lines. As with the JA-deficient mutant dde2-2, the transgenic line 4-1 was insensitive to 50 μM methyl jasmonate, compared with the wild-type plants. The results suggest that GrTCP11 may be an important transcription factor for cotton fiber development, by negatively regulating JA biosynthesis and response.

Embryogenic Calli Induction and Salt Stress Response Revealed by RNA-Seq in Diploid Wild Species Gossypium sturtianum and Gossypium raimondii

Wild cotton species can contribute to a valuable gene pool for genetic improvement, such as genes related to salt tolerance. However, reproductive isolation of different species poses an obstacle to produce hybrids through conventional breeding. Protoplast fusion technology for somatic cell hybridization provides an opportunity for genetic manipulation and targeting of agronomic traits. Transcriptome sequencing analysis of callus under salt stress is conducive to study salt tolerance genes. In this study, calli were induced to provide materials for extracting protoplasts and also for screening salt tolerance genes. Calli were successfully induced from leaves of Gossypium sturtianum (C₁ genome) and hypocotyls of G. raimondii (D₅ genome), and embryogenic calli of G. sturtianum and G. raimondii were induced on a differentiation medium with different concentrations of 2, 4-D, KT, and IBA, respectively. In addition, embryogenic calli were also induced successfully from G. raimondii through suspension cultivation. Transcriptome sequencing analysis was performed on the calli of G. raimondii and G. sturtianum, which were treated with 200 mM NaCl at 0, 6, 12, 24, and 48 h, and a total of 12,524 genes were detected with different expression patterns under salt stress. Functional analysis showed that 3,482 genes, which were differentially expressed in calli of G. raimondii and G. sturtianum, were associated with biological processes of nucleic acid binding, plant hormone (such as ABA) biosynthesis, and signal transduction. We demonstrated that DEGs or TFs which related to ABA metabolism were involved in the response to salt stress, including xanthoxin dehydrogenase genes (ABA2), sucrose non-fermenting 1-related protein kinases (SnRK2), NAM, ATAT1/2, and CUC2 transcription factors (NAC), and WRKY class of zinc-finger proteins (WRKY). This research has successfully induced calli from two diploid cotton species and revealed new genes responding to salt stress in callus tissue, which will lay the foundation for protoplast fusion for further understanding of salt stress responses in cotton.

Chromosome Painting Based on Bulked Oligonucleotides in Cotton

Chromosome painting is one of the key technologies in cytogenetic research, which can accurately identify chromosomes or chromosome regions. Oligonucleotide (oligo) probes designed based on genome sequences have both flexibility and specificity, which would be ideal probes for fluorescence in situ hybridization (FISH) analysis of genome structure. In this study, the bulked oligos of the two arms of chromosome seven of cotton were developed based on the genome sequence of Gossypium raimondii (DD, 2n = 2× = 26), and each arm contains 12,544 oligos. Chromosome seven was easily identified in both D genome and AD genome cotton species using the bulked chromosome-specific painting probes. Together with 45S ribosomal DNA (rDNA) probe, the chromosome-specific painting probe was also successfully used to correct the chromosomal localization of 45S rDNA in G. raimondii. The study reveals that bulked oligos specific to a chromosome is a useful tool for chromosome painting in cotton.

De Novo Genome Sequence Assemblies of Gossypium raimondii and Gossypium turneri

Cotton is an agriculturally important crop. Because of its importance, a genome sequence of a diploid cotton species (Gossypium raimondii, D-genome) was first assembled using Sanger sequencing data in 2012. Improvements to DNA sequencing technology have improved accuracy and correctness of assembled genome sequences. Here we report a new de novo genome assembly of G. raimondii and its close relative G. turneri. The two genomes were assembled to a chromosome level using PacBio long-read technology, HiC, and Bionano optical mapping. This report corrects some minor assembly errors found in the Sanger assembly of G. raimondii. We also compare the genome sequences of these two species for gene composition, repetitive element composition, and collinearity. Most of the identified structural rearrangements between these two species are due to intra-chromosomal inversions. More inversions were found in the G. turneri genome sequence than the G. raimondii genome sequence. These findings and updates to the D-genome sequence will improve accuracy and translation of genomics to cotton breeding and genetics.

Genome-Wide Analysis of Multidrug and Toxic Compound Extrusion (MATE) Family in Gossypium raimondii and Gossypium arboreum and Its Expression Analysis Under Salt, Cadmium, and Drought Stress

The extrusion of toxins and substances at a cellular level is a vital life process in plants under abiotic stress. The multidrug and toxic compound extrusion (MATE) gene family plays a large role in the exportation of toxins and other substrates. We carried out a genome-wide analysis of MATE gene families in Gossypium raimondii and Gossypium arboreum and assessed their expression levels under salt, cadmium and drought stresses. We identified 70 and 68 MATE genes in G. raimondii and G. arboreum, respectively. The majority of the genes were predicted to be localized within the plasma membrane, with some distributed in other cell parts. Based on phylogenetic analysis, the genes were subdivided into three subfamilies, designated as M1, M2 and M3. Closely related members shared similar gene structures, and thus were highly conserved in nature and have mainly evolved through purifying selection. The genes were distributed in all chromosomes. Twenty-nine gene duplication events were detected, with segmental being the dominant type. GO annotation revealed a link to salt, drought and cadmium stresses. The genes exhibited differential expression, with GrMATE18, GrMATE34, GaMATE41 and GaMATE51 significantly upregulated under drought, salt and cadmium stress, and these could possibly be the candidate genes. Our results provide the first data on the genome-wide and functional characterization of MATE genes in diploid cotton, and are important for breeders of more stress-tolerant cotton genotypes.

Genome-wide analysis of the MADS-box gene family in polyploid cotton (Gossypium hirsutum) and in its diploid parental species (Gossypium arboreum and Gossypium raimondii)

The MADS-box gene family encodes transcription factors that share a highly conserved domain known to bind to DNA. Members of this family control various processes of development in plants, from root formation to fruit ripening. In this work, a survey of diploid (Gossypium raimondii and Gossypium arboreum) and tetraploid (Gossypium hirsutum) cotton genomes found a total of 147, 133 and 207 MADS-box genes, respectively, distributed in the MIKC, Mα, Mβ, Mγ, and Mδ subclades. A comparative phylogenetic analysis among cotton species, Arabidopsis, poplar and grapevine MADS-box homologous genes allowed us to evaluate the evolution of each MADS-box lineage in cotton plants and identify sequences within well-established subfamilies. Chromosomal localization and phylogenetic analysis revealed that G. raimondii and G. arboreum showed a conserved evolution of the MIKC subclade and a distinct pattern of duplication events in the Mα, Mγ and Mδ subclades. Additionally, G. hirsutum showed a combination of its parental subgenomes followed by a distinct evolutionary history including gene gain and loss in each subclade. qPCR analysis revealed the expression patterns of putative homologs in the AP1, AP3, AGL6, SEP4, AGL15, AG, AGL17, TM8, SVP, SOC and TT16 subfamilies of G. hirsutum. The identification of putative cotton orthologs is discussed in the light of evolution and gene expression data from other plants. This analysis of the MADS-box genes in Gossypium species opens an avenue to understanding the origin and evolution of each gene subfamily within diploid and polyploid species and paves the way for functional studies in cotton species.

Genome-wide Identification and analysis of the stress-resistance function of the TPS (Trehalose-6-Phosphate Synthase) gene family in cotton

BACKGROUND

Trehalose (a-D-glucopyranosyl a-D-glucopyranoside) is a nonreducing disaccharide and is widely distributed in bacteria, fungi, algae, plants and invertebrates. In the study, the identification of trehalose-6-phosphate synthase (TPS) genes stress-related in cotton, and the genetic structure analysis and molecular evolution analysis of TPSs were conducted with bioinformatics methods, which could lay a foundation for further research of TPS functions in cotton.

RESULTS

The genome information of Gossypium raimondii (group D), G. arboreum L. (group A), and G. hirsutum L. (group AD) was used in the study. Fifty-three TPSs were identified comprising 15 genes in group D, 14 in group A, and 24 in group AD. Bioinformatics methods were used to analyze the genetic structure and molecular evolution of TPSs. Real-time PCR analysis was performed to investigate the expression patterns of gene family members. All TPS family members in cotton can be divided into two subfamilies: Class I and Class II. The similarity of the TPS sequence is high within the same species and close within their family relatives. The genetic structures of two TPS subfamily members are different, with more introns and a more complicated gene structure in Class I. There is a TPS domain(Glyco transf_20) at the N-terminal in all TPS family members and a TPP domain(Trehalose_PPase) at the C-terminal in all except GrTPS6, GhTPS4, and GhTPS9. All Class II members contain a UDP-forming domain. The responses to environmental stresses showed that stresses could induce the expression of TPSs but the expression patterns vary with different stresses.

CONCLUSIONS

The distribution of TPSs varies with different species but is relatively uniform on chromosomes. Genetic structure varies with different gene members, and expression levels vary with different stresses and exhibit tissue specificity. The upregulated genes in upland cotton TM-1 is significantly more than that in G. raimondii and G. arboreum L. Shixiya 1.

Expression profiles of miRNAs in Gossypium raimondii

miRNAs are a class of conserved, small, endogenous, and non-protein-coding RNA molecules with 20-24 nucleotides (nt) in length that function as post-transcriptional modulators of gene expression in eukaryotic cells. Functional studies have demonstrated that plant miRNAs are involved in the regulation of a wide range of plant developmental processes. To date, however, no research has been carried out to study the expression profiles of miRNAs in Gossypium raimondii, a model cotton species. We selected 16 miRNAs to profile their tissue-specific expression patterns in G. raimondii four different tissues, and these miRNAs are reported to play important roles in plant growth and development. Our results showed that the expression levels of these miRNAs varied significantly from one to another in a tissue-dependent manner. Eight miRNAs, including miR-159, miR-162, miR-164, miR-172, miR-390, miR-395, miR-397, and miR-398, exhibited exclusively high expression levels in flower buds, suggesting that these miRNAs may play significant roles in floral development. The expression level of miR-164 was relatively high in shoots beside flower buds, implying that the function of miR-164 is not only limited to floral development but it may also play an important role in shoot development. Certain miRNAs such as miR-166 and miR-160 were extremely highly expressed in all of the four tissues tested compared with other miRNAs investigated, suggesting that they may play regulatory roles at multiple development stages. This study will contribute to future studies on the functional characterization of miRNAs in cotton.

Distribution and characterization of simple sequence repeats in Gossypium raimondii genome

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