The GaKAN2, a KANADI transcription factor, modulates stem trichomes in Gossypium arboreum.

KEY MESSAGE: GaKAN2, a member of the KANADI family, was found to be widely expressed in the cotton tissues and regulates trichome development through complex pathways. Cotton trichomes are believed to be the defense barrier against insect pests. Cotton fiber and trichomes are single-cell epidermal extensions with shared regulatory mechanisms. Despite several studies underlying mechanism of trichome development remains elusive. The KANADI is one of the key transcription factors (TFs) family, regulating Arabidopsis trichomes growth. However, the function of KANADI genes in cotton remains unknown. In the current study genome-wide scanning, transcriptomic analysis, gene silencing, subcellular localization, and yeast two-hybrid techniques were employed to decipher the function of KANADI TFs family genes in cotton crop. A total of 7 GaKAN genes were found in the Gossypium arboreum. Transcriptomic data revealed that these genes were significantly expressed in stem and root. Moreover, GaKAN2 was widely expressed in other tissues also. Subsequently, we selected GaKAN2 to validate the function of KANADI genes. Silencing of GaKAN2 resulted in a 24.99% decrease in single-cell trichomes and an 11.33% reduction in internodal distance, indicating its potential role in regulating trichomes and plant growth. RNA-Seq analysis elucidated that GaSuS and GaERS were the downstream genes of GaKAN2. The transcriptional activation and similarity in silencing phenotype between GaKAN2 and GaERS suggested that GaKAN2 regulates trichomes development through GaERS. Moreover, KEGG analysis revealed that a significant number of genes were enriched in the biosynthesis of secondary metabolites and plant hormone signal transduction pathways, thereby suggesting that GaKAN2 regulates the stem trichomes and plant growth. The GFP subcellular localization and yeast transcriptional activation analysis elucidated that GaKAN2 was located in the nucleus and capable of regulating the transcription of downstream genes. This study elucidated the function and characteristics of the KANADI gene family in cotton, providing a fundamental basis for further research on GaKAN2 gene in cotton plant trichomes and plant developmental processes.

Single-cell transcriptomic analysis reveals the developmental trajectory and transcriptional regulatory networks of pigment glands in Gossypium bickii.

Comprehensive utilization of cottonseeds is limited by the presence of pigment glands and its inclusion gossypol. The ideal cotton has glandless seeds but a glanded plant, a trait found in only a few Australian wild cotton species, including Gossypium bickii. Introgression of this trait into cultivated species has proved to be difficult. Understanding the biological processes toward pigment gland morphogenesis and the associated underlying molecular mechanisms will facilitate breeding of cultivated cotton varieties with the trait of glandless seeds and glanded plant. In this study, single-cell RNA sequencing (scRNA-seq) was performed on 12 222 protoplasts isolated from cotyledons of germinating G. bickii seeds 48 h after imbibition. Clustered into 14 distinct clusters unsupervisedly, these cells could be grouped into eight cell populations with the assistance of known cell marker genes. The pigment gland cells were well separated from others and could be separated into pigment gland parenchyma cells, secretory cells, and apoptotic cells. By integrating the pigment gland cell developmental trajectory, transcription factor regulatory networks, and core transcription factor functional validation, we established a model for pigment gland formation. In this model, light and gibberellin were verified to promote the formation of pigment glands. In addition, three novel genes, GbiERF114 (ETHYLENE RESPONSE FACTOR 114), GbiZAT11 (ZINC FINGER OF ARABIDOPSIS THALIANA 11), and GbiNTL9 (NAC TRANSCRIPTION FACTOR-LIKE 9), were found to affect pigment gland formation. Collectively, these findings provide new insights into pigment gland morphogenesis and lay the cornerstone for future cotton scRNA-seq investigations.

A reference-grade genome assembly for Gossypium bickii and insights into its genome evolution and formation of pigment glands and gossypol.

The pigment gland is a morphological characteristic of Gossypium and its related genera. Gossypium bickii (G1) is characterized by delayed pigment gland morphogenesis in the cotyledons. In this study, a reference-grade genome of G1 was generated, and comparative genomics analysis showed that G1 was closest to Gossypium australe (G2), followed by A- and D-genome species. Two large fragment translocations in chromosomes 5 and 13 were detected between the G genome and other Gossypium genomes and were unique to the G1 and G2 genomes. Compared with the G2 genome, two large fragment inversions in chromosomes 12 and 13 were detected in G1. According to the phylogeny, divergence time, and similarity analysis of nuclear and chloroplast genomes, G1 was formed by hybridization between Gossypium sturtianum (C1) and a common ancestor of G2 and Gossypium nelsonii (G3). The coordinated expression patterns of pigment gland formation (GoPGF) and gossypol biosynthesis genes in G1 were verified to be consistent with its phenotype, and nine genes that were related to the process of pigment gland formation were identified. A novel gene, GbiCYP76B6, regulated by GoPGF, was found to affect gossypol biosynthesis. These findings offer insights into the origin and evolution of G1 and its mechanism of pigment gland formation and gossypol biosynthesis.

Deficiency of a peroxisomal NADP-isocitrate dehydrogenase leads to dwarf plant and defect seed in upland cotton.

The NADP-isocitrate dehydrogenase-encoded gene GH_D13G1452 with a C-terminus tripeptide Proline-Lysine-Leucine was localized in the peroxisome. It was highly expressed in stems and ovules of 15 days post-anthesis and responded to multiple external stimuli in upland cotton. An upland cotton mutant (Ghpericdh) was identified by flanking sequence amplification and genome variation detection that exogenous sequence was inserted in the middle of the 12th intron of GH_D13G1452, resulting in the deficiency of gene expression. The Ghpericdh mutant displayed a dwarf plant phenotype when grown under field or greenhouse conditions, and GH_D13G1452 functioned as an incomplete dominance on plant height. The germination rate of mutant seed from greenhouse-grown plants was dramatically lower than that from field-grown plants, which indicated that GhperICDH plays a critical role in seed maturation and germination. Therefore, GH_D13G1452 is indispensable in the development of stems and seeds and functions in the adaptability of cotton to the environment. The Ghpericdh mutant provides insight into the function of peroxisomal ICDH and may contribute to the genetic improvement in cotton.

A Modified Actin (Gly65Val Substitution) Expressed in Cotton Disrupts Polymerization of Actin Filaments Leading to the Phenotype of Ligon Lintless-1 (Li1) Mutant.

Dynamic remodeling of the actin cytoskeleton plays a central role in the elongation of cotton fibers, which are the most important natural fibers in the global textile industry. Here, a high-resolution mapping approach combined with comparative sequencing and a transgenic method revealed that a G65V substitution in the cotton actin Gh_D04G0865 (GhACT17D in the wild-type) is responsible for the Gossypium hirsutum Ligon lintless-1 (Li1) mutant (GhACT17DM). In the mutant GhACT17DM from Li1 plant, Gly65 is substituted with valine on the lip of the nucleotide-binding domain of GhACT17D, which probably affects the polymerization of F-actin. Over-expression of GhACT17DM, but not GhACT17D, driven by either a CaMV35 promoter or a fiber-specific promoter in cotton produced a Li1-like phenotype. Compared with the wild-type control, actin filaments in Li1 fibers showed higher growth and shrinkage rates, decreased filament skewness and parallelness, and increased filament density. Taken together, our results indicate that the incorporation of GhACT17DM during actin polymerization disrupts the establishment and dynamics of the actin cytoskeleton, resulting in defective fiber elongation and the overall dwarf and twisted phenotype of the Li1 mutant.

Differentiation in the genetic basis of stem trichome development between cultivated tetraploid cotton species.

BACKGROUND: Cotton stem trichomes and seed fibers are each single celled structures formed by protrusions of epidermal cells, and were found sharing the overlapping molecular mechanism. Compared with fibers, cotton stem trichomes are more easily observed, but the molecular mechanisms underlying their development are still poorly understood. RESULTS: In this study, Gossypium hirsutum (Gh) and G. barbadense (Gb) were found to differ greatly in percentages of varieties/accessions with glabrous stems and in trichome density, length, and number per trichopore. Gh varieties normally had long singular and clustered trichomes, while Gb varieties had short clustered trichomes. Genetic mapping using five F2 populations from crosses between glabrous varieties and those with different types of stem trichomes revealed that much variation among stem trichome phenotypes could be accounted for by different combinations of genes/alleles on Chr. 06 and Chr. 24. The twenty- six F1 generations from crosses between varieties with different types of trichomes had varied phenotypes, further suggesting that the trichomes of tetraploid cotton were controlled by different genes/alleles. Compared to modern varieties, a greater proportion of Gh wild accessions were glabrous or had shorter and denser trichomes; whereas a smaller proportion of Gb primitive accessions had glabrous stems. A close correlation between fuzz fiber number and stem trichome density was observed in both Gh and Gb primitive accessions and modern varieties. CONCLUSION: Based on these findings, we hypothesize that stem trichomes evolved in parallel with seed fibers during the domestication of cultivated tetraploid cotton. In addition, the current results illustrated that stem trichome can be used as a morphological index of fiber quality in cotton conventional breeding.

GaHD1, a candidate gene for the Gossypium arboreum SMA-4 mutant, promotes trichome and fiber initiation by cellular H2O2 and Ca2+ signals.

Cotton fibers are initiated from the epidermal cells of the ovule before or on the day of anthesis. Gossypium arboreum SMA-4 mutant contains recessive mutation (sma-4(ha)) and has the phenotypes of fibreless seeds and glabrous stems. In this study, fine mapping and alternative splicing analysis indicated a nucleotide substitution (AG → AC) at splicing site in a homeodomain-leucine zipper IV family gene (GaHD1) might cause gene A3S (Alternative 3' splicing) mistake, suggested that GaHD1 was the candidate gene of sma-4(ha). Many genes related to the fiber initiation are identified to be differentially expressed in the mutant which could result in the blocked fiber initiation signals such as H2O2, or Ca in the mutant. Further comparative physiological analysis of H2O2 production and Ca2+ flux in the SMA-4 and wide type cotton confirmed that H2O2 and Ca were important fiber initiation signals and regulated by GaHD1. The in vitro ovule culture of the mutant with hormones recovered the fibered phenotype coupled with the restoration of these signals. Overexpressing of GaHD1 in Arabidopsis increased trichome densities on the sepal, leaf, and stem tissues while transient silencing of the GaHD1 gene in G. arboreum reduced the trichome densities. These phenotypes indicated that GaHD1 is the candidate gene of SMA-4 with a crucial role in acting upstream molecular switch of signal transductions for cotton trichome and fiber initiations.

Preferential insertion of a Ty1 LTR-retrotransposon into the A sub-genome's HD1 gene significantly correlated with the reduction in stem trichomes of tetraploid cotton.

Stem trichomes and seed fibers originate from epidermal cells and partially share a regulatory pathway at the molecular level. In Gossypium barbadense, two insertions of a Ty1 long-terminal repeat-retrotransposon [transposable element TE1 and TE2] in a homeodomain-leucine zipper gene (HD1) result in glabrous stems. The primers used to identify the TE insertions in G. barbadense were applied to screen for the same events in 81 modern G. hirsutum varieties and 31 wild races. Three wild races were found carrying the same TEs as G. barbadense. However, the TE insertions in two of these wild races occurred at different sites (4th exon), therefore, named TE3, while the TE in the other wild race occurred at the same site as TE2. An RNA sequencing and qRT-PCR analysis indicated that the loss of HD1 function was caused by the TE insertion. Genetic mapping revealed a strong association between glabrous stems and TE3 insertions, confirming that HD1 is a critical gene for stem trichome initiation in G. hirsutum, as in G. barbadense. Using the long-terminal repeat sequence as a query to search against the Texas Marker-1 reference genome sequence, we found that the TE occurred after tetraploid cotton formation and evolved at different rates in G. hirsutum and G. barbadense. Interestingly, at least three independent insertion events of the same retrotransposon occurred preferentially in the A sub-genome's HD1 gene, but not in the D sub-genome of G. hirsutum or G. barbadense, suggesting that an unknown TE insertion mechanism and resultant gene function changes may have hastened cotton speciation.

The Hairless Stem Phenotype of Cotton (Gossypium barbadense) Is Linked to a Copia-Like Retrotransposon Insertion in a Homeodomain-Leucine Zipper Gene (HD1).

Cotton (Gossypium) stem trichomes are mostly single cells that arise from stem epidermal cells. In this study, a homeodomain-leucine zipper gene (HD1) was found to cosegregate with the dominant trichome locus previously designated as T1 and mapped to chromosome 6. Characterization of HD1 orthologs revealed that the absence of stem trichomes in modern Gossypium barbadense varieties is linked to a large retrotransposon insertion in the ninth exon, 2565 bp downstream from the initial codon in the At subgenome HD1 gene (At-GbHD1). In both the At and Dt subgenomes, reduced transcription of GbHD1 genes is caused by this insertion. The disruption of At-HD1 further affects the expression of downstream GbMYB25 and GbHOX3 genes. Analyses of primitive cultivated accessions identified another retrotransposon insertion event in the sixth exon of At-GbHD1 that might predate the previously identified retrotransposon in modern varieties. Although both retrotransposon insertions results in similar phenotypic changes, the timing of these two retrotransposon insertion events fits well with our current understanding of the history of cotton speciation and dispersal. Taken together, the results of genetics mapping, gene expression and association analyses suggest that GbHD1 is an important component that controls stem trichome development and is a promising candidate gene for the T1 locus. The interspecific phenotypic difference in stem trichome traits also may be attributable to HD1 inactivation associated with retrotransposon insertion.

Genetic fine mapping and candidate gene analysis of the Gossypium hirsutum Ligon lintless-1 (Li1) mutant on chromosome 22(D).

Ligon lintless-1 (Li1) is a Gossypium hirsutum mutant that is controlled by a dominant gene that arrests the development of cotton fiber after anthesis. Two F2 mapping populations were developed from mutant (Li1 × H7124) F1 plants in 2012 and 2013; each was composed of 142 and 1024 plants, respectively. Using these populations, Li1 was mapped to a 0.3-cM region in which nine single-strand conformation polymorphism markers co-segregated with the Li1 locus. In the published G. raimondii genome, these markers were mapped to a region of about 1.2 Mb (the Li1 region) and were separated by markers that flanked the Li1 locus in the genetic map, dividing the Li1 region into three segments. Thirty-six genes were annotated by the gene prediction software FGENESH (Softberry) in the Li1 region. Twelve genes were candidates of Li1, while the remaining 24 genes were identified as transposable elements, DNA/RNA polymerase superfamily or unknown function genes. Among the 12 candidate genes, those encoding ribosomal protein s10, actin protein, ATP synthase, and beta-tubulin 5 were the most-promising candidates of the Li1 mutant because the function of these genes is closely related to fiber development. High-throughput RNA sequencing and quantitative PCR revealed that these candidate genes had obvious differential gene expression between mutant and wild-type plants at the fiber elongation stage, strengthening the inference that they could be the most likely candidate gene of the Li1 mutant phenotype.

Molecular characterization of a transcriptionally active Ty1/copia-like retrotransposon in Gossypium.

KEY MESSAGE: A transcriptionally active Ty1/copia -like retrotransposon was identified in the genome of Gossypium barbadense. The different heat activation of this element was observed in two tetraploid cotton species. Most retrotransposons from plants are transcriptionally silent, or activated under certain conditions. Only a small portion of elements are transcriptionally active under regular condition. A long terminal repeat (LTR) retrotransposon was isolated from the cultivated Sea Island cotton (H7124) genome during the investigation of the function of a homeodomain leucine zipper gene (HD1) in trichome growth. Insertion of this element in HD1 gene of At sub-genome was related to the trichomeless stem in Gossypium barbadense. The element, named as GBRE-1, had all features of a typical Ty1/copia retrotransposon and possessed high similarity to the members of ONSEN retrotransposon family. It was 4997 bp long, comprising a single 4110 bp open reading frame, which encoded 1369 amino acids including the conserved domains of gag and pol. The expression of GBRE-1 was detected under regular condition in G. barbadense and G. hirsutum, and its expression level was increased under heat-stress condition in G. hirsutum. Besides, its expression pattern was similar to that of the ONSEN retrotransposon. Abundant cis-regulatory motifs related to stress-response and transcriptional regulation were found in the LTR sequence. These results suggested that GBRE-1 was a transcriptionally active retrotransposon in Gossypium. To our knowledge, this is the first report of the isolation of a complete Ty1/copia-type retrotransposon with present-day transcriptional activity in cotton.

Genome-wide investigation and transcriptome analysis of the WRKY gene family in Gossypium.

WRKY transcription factors play important roles in various stress responses in diverse plant species. In cotton, this family has not been well studied, especially in relation to fiber development. Here, the genomes and transcriptomes of Gossypium raimondii and Gossypium arboreum were investigated to identify fiber development related WRKY genes. This represents the first comprehensive comparative study of WRKY transcription factors in both diploid A and D cotton species. In total, 112 G. raimondii and 109 G. arboreum WRKY genes were identified. No significant gene structure or domain alterations were detected between the two species, but many SNPs distributed unequally in exon and intron regions. Physical mapping revealed that the WRKY genes in G. arboreum were not located in the corresponding chromosomes of G. raimondii, suggesting great chromosome rearrangement in the diploid cotton genomes. The cotton WRKY genes, especially subgroups I and II, have expanded through multiple whole genome duplications and tandem duplications compared with other plant species. Sequence comparison showed many functionally divergent sites between WRKY subgroups, while the genes within each group are under strong purifying selection. Transcriptome analysis suggested that many WRKY genes participate in specific fiber development processes such as fiber initiation, elongation and maturation with different expression patterns between species. Complex WRKY gene expression such as differential Dt and At allelic gene expression in G. hirsutum and alternative splicing events were also observed in both diploid and tetraploid cottons during fiber development process. In conclusion, this study provides important information on the evolution and function of WRKY gene family in cotton species.

Gene expression profile analysis of Ligon lintless-1 (Li1) mutant reveals important genes and pathways in cotton leaf and fiber development.

Ligon lintless-1 (Li1) is a monogenic dominant mutant of Gossypium hirsutum (upland cotton) with a phenotype of impaired vegetative growth and short lint fibers. Despite years of research involving genetic mapping and gene expression profile analysis of Li1 mutant ovule tissues, the gene remains uncloned and the underlying pathway of cotton fiber elongation is still unclear. In this study, we report the whole genome-level deep-sequencing analysis of leaf tissues of the Li1 mutant. Differentially expressed genes in leaf tissues of mutant versus wild-type (WT) plants are identified, and the underlying pathways and potential genes that control leaf and fiber development are inferred. The results show that transcription factors AS2, YABBY5, and KANDI-like are significantly differentially expressed in mutant tissues compared with WT ones. Interestingly, several fiber development-related genes are found in the downregulated gene list of the mutant leaf transcriptome. These genes include heat shock protein family, cytoskeleton arrangement, cell wall synthesis, energy, H2O2 metabolism-related genes, and WRKY transcription factors. This finding suggests that the genes are involved in leaf morphology determination and fiber elongation. The expression data are also compared with the previously published microarray data of Li1 ovule tissues. Comparative analysis of the ovule transcriptomes of Li1 and WT reveals that a number of pathways important for fiber elongation are enriched in the downregulated gene list at different fiber development stages (0, 6, 9, 12, 15, 18dpa). Differentially expressed genes identified in both leaf and fiber samples are aligned with cotton whole genome sequences and combined with the genetic fine mapping results to identify a list of candidate genes for Li1.

Genomic characterization of a novel dsRNA virus detected in the phytopathogenic fungus Verticillium dahliae Kleb.

Four novel double-stranded RNA segments were detected in a Verticillium dahliae Kleb. strain (V. dahliae isolate 0-21), a causal fungal agent of Verticillium wilt disease of cotton. Each dsRNA genome segment contains a single large open reading frame (ORF) that encodes a distinctive protein with modest levels of sequence similarities to the corresponding putative proteins in the genus Chrysovirus. These include an RNA-dependent RNA polymerase (RdRp), a coat protein, an undefined replication-related protein and an ovarian tumor domain peptidase. Phylogenetic analysis of the four putative proteins unanimously indicated that they are evolutionarily related to viruses in Chrysovirus. The 5'- and 3'-untranslated regions of the four dsRNAs share highly similar internal sequence and contain conserved sequence stretches of UGAUAAAAAA(/U)UG(/U)AAAAA- (in the 5'-UTR) and -UUUACUACU (in the 3'-UTR), indicating that they have a common virus origin. Indeed, isometric virus-like particles (VLPs) with a diameter of approximately 34nm were extracted from the fungal mycelia, and the four dsRNA segments were also detected in the virus-like particle (VLP) fraction. These results suggest that the mycovirus with four different dsRNA genome segments from the fungal isolate 0-21 is a new member of the genus Chrysovirus. We named the virus Verticillium dahliae chrysovirus 1 (VdCV1).