Members of the MYBMIXTA-like transcription factors may orchestrate the initiation of fiber development in cotton seeds

Deciphering the complex cotton genome for improving fiber traits and abiotic stress resilience in sustainable agriculture

BACKGROUND

Understanding the complex cotton genome is of paramount importance in devising a strategy for sustainable agriculture. Cotton is probably the most economically important cash crop known for its cellulose-rich fiber content. The cotton genome has become an ideal model for deciphering polyploidization due to its polyploidy, setting it apart from other major crops. However, the main challenge in understanding the functional and regulatory functions of many genes in cotton is still the complex cotton polyploidy genome, which is not limited to a single role. Cotton production is vulnerable to the sensitive effects of climate change, which can alter or aggravate soil, pests, and diseases. Thus, conventional plant breeding coupled with advanced technologies has led to substantial progress being made in cotton production.

GENOMICS APPROACHES IN COTTON

In the frontier areas of genomics research, cotton genomics has gained momentum accomplished by robust high-throughput sequencing platforms combined with novel computational tools to make the cotton genome more tractable. Advances in long-read sequencing have allowed for the generation of the complete set of cotton gene transcripts giving incisive scientific knowledge in cotton improvement. In contrast, the integration of the latest sequencing platforms has been used to generate multiple high-quality reference genomes in diploid and tetraploid cotton. While pan-genome and 3D genomic studies are still in the early stages in cotton, it is anticipated that rapid advances in sequencing, assembly algorithms, and analysis pipelines will have a greater impact on advanced cotton research.

CONCLUSIONS

This review article briefly compiles substantial contributions in different areas of the cotton genome, which include genome sequencing, genes, and their molecular regulatory networks in fiber development and stress tolerance mechanism. This will greatly help us in understanding the robust genomic organization which in turn will help unearth candidate genes for functionally important agronomic traits.

Poplar coma morphogenesis and miRNA regulatory networks by combining ovary tissue sectioning and deep sequencing

Poplar coma, commonly referred to as "seed hairs", is a tuft of trichomes attached to the seed coat that helps seed dispersal. However, they can also trigger health impacts for humans, including sneezing, shortness of breath, and skin irritation. Despite efforts to study the regulatory mechanism of herbaceous trichome formation, poplar coma remains poorly understood. In this study, we showed that the epidermal cells of the funiculus and placenta are the origin of poplar coma based on observations of paraffin sections. Small RNA (sRNA) and degradome libraries were also constructed at three stages of poplar coma development, including initiation and elongation stages. Based on 7,904 miRNA-target pairs identified by small RNA and degradome sequencing, we constructed a miRNA-transcript factor and a stage-specific miRNA regulatory network. By combining paraffin section observation and deep sequencing, our research will provide greater insight into the molecular mechanisms of poplar coma development.

Genome-wide analysis of R2R3-MYB genes in cultivated peanut (Arachis hypogaea L.): Gene duplications, functional conservation, and diversification

The cultivated Peanut (Arachis hypogaea L.), an important oilseed and edible legume, are widely grown worldwide. The R2R3-MYB transcription factor, one of the largest gene families in plants, is involved in various plant developmental processes and responds to multiple stresses. In this study we identified 196 typical R2R3-MYB genes in the genome of cultivated peanut. Comparative phylogenetic analysis with Arabidopsis divided them into 48 subgroups. The motif composition and gene structure independently supported the subgroup delineation. Collinearity analysis indicated polyploidization, tandem, and segmental duplication were the main driver of the R2R3-MYB gene amplification in peanut. Homologous gene pairs between the two subgroups showed tissue specific biased expression. In addition, a total of 90 R2R3-MYB genes showed significant differential expression levels in response to waterlogging stress. Furthermore, we identified an SNP located in the third exon region of AdMYB03-18 (AhMYB033) by association analysis, and the three haplotypes of the SNP were significantly correlated with total branch number (TBN), pod length (PL) and root-shoot ratio (RS ratio), respectively, revealing the potential function of AdMYB03-18 (AhMYB033) in improving peanut yield. Together, these studies provide evidence for functional diversity in the R2R3-MYB genes and will contribute to understanding the function of R2R3-MYB genes in peanut.

Single-cell RNA-seq reveals fate determination control of an individual fibre cell initiation in cotton (Gossypium hirsutum)

Cotton fibre is a unicellular seed trichome, and lint fibre initials per seed as a factor determines fibre yield. However, the mechanisms controlling fibre initiation from ovule epidermis are not understood well enough. Here, with single-cell RNA sequencing (scRNA-seq), a total of 14 535 cells were identified from cotton ovule outer integument of Xu142_LF line at four developmental stages (1.5, 1, 0.5 days before anthesis and the day of anthesis). Three major cell types, fibre, non-fibre epidermis and outer pigment layer were identified and then verified by RNA in situ hybridization. A comparative analysis on scRNA-seq data between Xu142 and its fibreless mutant Xu142 fl further confirmed fibre cluster definition. The developmental trajectory of fibre cell was reconstructed, and fibre cell was identified differentiated at 1 day before anthesis. Gene regulatory networks at four stages revealed the spatiotemporal pattern of core transcription factors, and MYB25-like and HOX3 were demonstrated played key roles as commanders in fibre differentiation and tip-biased diffuse growth respectively. A model for early development of a single fibre cell was proposed here, which sheds light on further deciphering mechanism of plant trichome and the improvement of cotton fibre yield.

The R2R3-MYB gene CgMYB4 is involved in the regulation of cell differentiation and fiber development in the stamens of Chelone glabra L

A Plantaginaceae flowering plant, Chelone glabra, is different from Arabidopsis thaliana and cotton (Gossypium hirsutum), as it produces fibers on the anther surface. However, the evolutionary molecular mechanism of how fiber development is controlled in the stamen is unclear. MYB genes are essential transcription factors for trichome and fiber development in plants. In this study, we isolated 29 MYB domain-containing sequences using early-stage anthers and several sets of degenerated primers conserved in the R2R3 domain of the MYB transcription factor. Among them, CgMYB4 is an R2R3-MYB gene encoding 281 amino acids. Phylogenetic analysis showed that CgMYB4 is closely related to GhMYB25L/AmMIXTA, which controls fiber initiation and development in cotton and epidermal cell differentiation in the petals of Antirrhinum. Semiquantitative RT-PCR analysis showed that CgMYB4 is strongly expressed at the stamens and carpels. Overexpression of CgMYB4 significantly enhanced root hair formation in transformed hairy roots, contrary to the root hair numbers, which were reduced in silenced CgMYB4 hairy roots. Moreover, overexpression of CgMYB4 also evidently promoted fiber development at filaments and conical cell-like epidermal cell increases at the anther wall. Our results showed that CgMYB4 is an R2R3-MYB gene and is positively involved in regulating cell division and fiber differentiation in the early stages of stamen development in C. glabra.

Conical petal epidermal cells, regulated by the MYB transcription factor MIXTA, have an ancient origin within the angiosperms

Conical epidermal cells occur on the tepals (perianth organs, typically petals and/or sepals) of the majority of animal-pollinated angiosperms, where they play both visual and tactile roles in pollinator attraction, providing grip to foraging insects, and enhancing colour, temperature, and hydrophobicity. To explore the evolutionary history of conical epidermal cells in angiosperms, we surveyed the tepal epidermis in representative species of the ANA-grade families, the early-diverging successive sister lineages to all other extant angiosperms, and analysed the function of a candidate regulator of cell outgrowth from Cabomba caroliniana (Nymphaeales). We identified conical cells in at least two genera from different families (Austrobaileya and Cabomba). A single SBG9 MYB gene was isolated from C. caroliniana and found to induce strong differentiation of cellular outgrowth, including conical cells, when ectopically expressed in Nicotiana tabacum. Ontogenetic analysis and quantitative reverse transcription-PCR established that CcSBG9A1 is spatially and temporally expressed in a profile which correlates with a role in conical cell development. We conclude that conical or subconical cells on perianth organs are ancient within the angiosperms and most probably develop using a common genetic programme initiated by a SBG9 MYB transcription factor.

A Retrotransposon Insertion in GhMML3_D12 Is Likely Responsible for the Lintless Locus li3 of Tetraploid Cotton

Cotton (Gossypium) seed fibers can be divided into lint (long) or fuzz (very short). Using fiberless (fuzzless-lintless) mutants, the lint initiation gene Li₃ was identified by map-based cloning. The gene is an R2R3-MYB transcription factor located on chromosome D12 (GhMML3_D12). Sequence analysis revealed that li₃ is a loss-of-function allele containing a retrotransposon insertion in the second exon that completely blocks the gene’s expression. The genetic loci n₂ and n₃ underlying the recessive fuzzless phenotype in Gossypium hirsutum were also mapped. The genomic location of n₃ overlapped with that of the dominant fuzzless locus N₁, and n₃ appeared to be a loss-of-function allele caused by a single nucleotide polymorphism (SNP) mutation in the coding region of GhMML3_A12. The n₂ allele was found to be co-located with li₃ and originated from G. babardense. n₂ and li₃ are possibly the multiple alleles of the GhMML3_D12 gene. Genetic analysis showed that Li₃ and N₃ are a pair of homologs with additive effects for the initiation of fibers (fuzz or lint). In addition, the presence of another locus was speculated, and it appeared to show an inhibitory effect on the expression of GhMML3. These findings provide new information about the genetic factors affecting the initiation of fibers in cotton.

Elucidation of sequence polymorphism in fuzzless-seed cotton lines

Most commercially produced cotton cultivars have two types of fibers on the seed coat, short fuzz and long lint. Lint fiber is used in the textile industry, while fuzz is considered an undesirable trait. Both types of fibers are believed to be controlled by the same regulators; however, their mechanisms of actions are still obscure. Cotton fiber mutants provide an excellent system to study the genes that regulate fiber development. Here we described four uncharacterized and three previously reported cotton mutants with fuzzless seed phenotypes. To evaluate whether or not the genes previously associated with fuzzless seed phenotypes have mutations we sequenced whole genomic DNA of seven mutants and wild type varieties. We identified multiple polymorphic changes among the tested genes. Non-synonymous SNPs in the coding region of the MML3-A gene was common in the six mutant lines tested in this study, showing both dominant and recessive fuzzless phenotypes. We have mapped the locus of the causative mutation for one of the uncharacterized fuzzless lines using an F₂ population that originated from a cross between the dominant fuzzless mutant and a wild type. Further, we have clarified the current knowledge about the causative n₂ mutations by analyzing the sequence data and previously reported mapping data. The key genes and possible mechanisms of fiber differentiation are discussed in this study.

Improved lint yield under field conditions in cotton over-expressing transcription factors regulating fibre initiation

Only a few transcription factors (TFs) regulating which cells of the ovule epidermis differentiate into lint fibres have been identified in cotton (Gossypium hirsutum L.). In this study, the effect on lint yield and fibre quality of over-expressing three TFs in cotton, GhHD-1, GhMYB25 and GhMYB25Like, and their double and triple combinations, were evaluated in field experiments over two seasons. The expression of single or stacked TFs were all driven either by an ovule-specific promoter, FBP 7, or a constitutive promoter, Stunt 7, in a Coker 315 background. TF type, either singly or in combination, was found to be the most significant factor affecting lint yield. Among 64 transgenic lines tested, seven were higher yielding than null segregant lines in one or both seasons and were all from the sets with single and double over-expressed TF combinations. A reduced yield was associated with the set of triple combinations. The two most stable high yielding lines across the seasons recorded 12-22% higher yields than the nulls, although were not competitive to locally adapted commercial controls. Over-expression of TFs singly or in combination did not significantly alter fibre length and strength, but sometimes increased fibre micronaire. There were positive relationships between lint yield and lint percentage and lint yield and fibre density amongst the transgenic lines. Our preliminary results suggest that manipulating TF expression, either singly or in pairs, can increase the density of fibres initiated on developing seeds and fibre yields under field conditions while maintaining overall fibre quality.

Transcriptome Analysis Reveals a Gene Expression Pattern Associated with Fuzz Fiber Initiation Induced by High Temperature in Gossypium barbadense

Gossypium barbadense is an important source of natural textile fibers, as is Gossypium hirsutum. Cotton fiber development is often affected by various environmental factors, such as abnormal temperature. However, little is known about the underlying mechanisms of temperature regulating the fuzz fiber initiation. In this study, we reveal that high temperatures (HT) accelerate fiber development, improve fiber quality, and induced fuzz initiation of a thermo-sensitive G. barbadense variety L7009. It was proved that fuzz initiation was inhibited by low temperature (LT), and 4 dpa was the stage most susceptible to temperature stress during the fuzz initiation period. A total of 43,826 differentially expressed genes (DEGs) were identified through comparative transcriptome analysis. Of these, 9667 were involved in fiber development and temperature response with 901 transcription factor genes and 189 genes related to plant hormone signal transduction. Further analysis of gene expression patterns revealed that 240 genes were potentially involved in fuzz initiation induced by high temperature. Functional annotation revealed that the candidate genes related to fuzz initiation were significantly involved in the asparagine biosynthetic process, cell wall biosynthesis, and stress response. The expression trends of sixteen genes randomly selected from the RNA-seq data were almost consistent with the results of qRT-PCR. Our study revealed several potential candidate genes and pathways related to fuzz initiation induced by high temperature. This provides a new view of temperature-induced tissue and organ development in Gossypium barbadense.

Trichomes at the Base of the Petal Are Regulated by the Same Transcription Factors as Cotton Seed Fibers

Many polypetalous plants have a constriction at the base of the petal that leaves a small gap that can provide entry into the young flower bud before the reproductive organs are fully developed. In cotton (Gossypium hirsutum L.), this gap is occluded by tufts of short unicellular trichomes superficially resembling the fibers found on cotton seeds. We are just beginning to understand the developmental regulation of the seed fibers and have previously characterized several MIXTA-like MYB transcription factors (TFs) that are critical for correct seed fiber development but know little about the molecular regulation of other types of cotton trichomes. Here, using RNAi or dominant suppression transgenic cotton lines and natural fiber mutants, we investigated the development and regulation of the petal base trichomes. Petal base trichomes and seed trichomes were also examined across several different species within and outside of the Malvoideae. We found that the petal base trichomes are regulated by the same MYB TFs as cotton seed fibers and, since they are more widely distributed across different taxa than the seed fibers, could have preceded them in the evolution of these important textile fibers produced by some cotton species.

Open-Bud Duplicate Loci Are Identified as MML10s, Orthologs of MIXTA-Like Genes on Homologous Chromosomes of Allotetraploid Cotton

The open-bud (ob) mutants in cotton display abnormal flower buds with the stigma and upper anthers exposed before blooming. This characteristic is potentially useful for the efficient production of hybrid seeds. The recessive inheritance pattern of the ob phenotype in allotetraploid cotton is determined by duplicated recessive loci (ob1ob1ob2ob2). In this study, ob1, which is a MIXTA-like MYB gene on chromosome D13 (MML10_Dt), was identified by map-based cloning. In Gossypium barbadense (Gb) acc. 3-79, a single nucleotide polymorphism (SNP) (G/A) at the splice site of the first intron and an 8-bp deletion in the third exon of MML10_Dt were found, which are the causative mutations at the ob1 loci. A 1783-bp deletion that leads to the loss of the third exon and accounts for the causal variation at the ob2 loci was found in MML10_At of Gossypium hirsutum (Gh) acc. TM-1. The ob phenotype results from the combination of these two loss-of-function loci. Genotyping assays showed that the ob1 and ob2 loci appeared after the formation of allotetraploid cotton and were specific for Gb and Gh, respectively. All Gb lines and most Gh cultivars carry the single corresponding mutant alleles. Genome-wide transcriptome analysis showed that some of the MYB genes and genes related to cell wall biogenesis, trichome differentiation, cytokinin signal transduction, and cell division were repressed in the ob mutants, which may lead to suppression of petal growth. These findings should be of value for breeding superior ob lines in cotton.

A Comparative Genome-Wide Analysis of the R2R3-MYB Gene Family Among Four Gossypium Species and Their Sequence Variation and Association With Fiber Quality Traits in an Interspecific G. hirsutum × G. barbadense Population

Cotton (Gossypium spp.) is the most important natural fiber crop in the world. The R2R3-MYB gene family is a large gene family involved in many plant functions including cotton fiber development. Although previous studies have reported its phylogenetic relationships, gene structures, and expression patterns in tetraploid G. hirsutum and diploid G. raimondii, little is known about the sequence variation of the members between G. hirsutum and G. barbadense and their involvement in the natural quantitative variation in fiber quality and yield. In this study, a comprehensive genome-wide comparative analysis was performed among the four Gossypium species using whole genome sequences, i.e., tetraploid G. hirsutum (AD1) and G. barbadense (AD2) as well as their likely ancestral diploid extants G. raimondii (D5) and G. arboreum (A2), leading to the identification of 406, 393, 216, and 213 R2R3-MYB genes, respectively. To elucidate whether the R2R3-MYB genes are genetically associated with fiber quality traits, 86 R2R3-MYB genes were co-localized with quantitative trait loci (QTL) hotspots for fiber quality and yield, including 42 genes localized within the fiber length QTL hotspots, in interspecific G. hirsutum × G. barbadense populations. There were 20 interspecific nonsynonymous single-nucleotide polymorphism (SNP) sites between the two tetraploid cultivated species, of which 16 developed from 11 R2R3-MYB genes were significantly correlated with fiber quality and yield in a backcross inbred population (BIL) of G. hirsutum × G. barbadense in at least one of the four field tests. Taken together, these results indicate that the sequence variation in these 11 R2R3-MYB genes is associated with the natural variation (i.e., QTL) in fiber quality and yield. Moreover, the functional SNPs of five R2R3-MYB allele pairs from the AD1 and AD2 genomes were significantly correlated with the gene expression related to fiber quality in fiber development. The results will be useful in further elucidating the role of the R2R3-MYB genes during fiber development.

Genetic dissection of the fuzzless seed trait in Gossypium barbadense

Cotton fibres are single-celled trichomes arising from the epidermal cells of the seed coat and may be either long (lint) or very short (fuzz). The dominant fuzzless N1 of Gossypium hirsutum is a defective allele of the At-subgenome homoeolog of MYB25-like, but the genetic components underlying the recessive fuzzless trait from G. barbadense (Gb) are unknown. We have identified five genetic loci, including a major contributing locus containing MYB25-like_Dt, associated with Gb fuzzless seeds based on genotyping of fuzzy and fuzzless near isogenic lines (NILs) from an interspecies cross (G. barbadense × G. hirsutum). At 3 d post-anthesis when fuzz fibres are initiating, expression of MYB25-like_Dt was significantly lower in fuzzless NILs than in fuzzy seeded NILs, while higher MYB25-like_Dt expression was associated with more seed fuzz across different cotton genotypes. Phenotypic and genotypic analysis of MYB25-like homoeoalleles in cottons showing different fibre phenotypes and their crossing progeny indicated that both MYB25-like_At and MYB25-like_Dt are associated with lint development, and that fuzz development is mainly determined by the expression level of MYB25-like_Dt at ~3 d post-anthesis. Expression of Gb fuzzless seeds depends on genetic background and interactions amongst the multiple loci identified. MYB25-like_Dt is one of the best candidates for N2.

Genetics and evolution of MIXTA genes regulating cotton lint fiber development

Cotton, with cellulose-enriched mature fibers, is the largest source of natural textiles. Through a map-based cloning strategy, we isolated an industrially important lint fiber development gene (Li₃ ) that encodes an MYB-MIXTA-like transcription factor (MML) on chromosome D12 (GhMML4_D12). Virus-induced gene silencing or decreasing the expression of the GhMML4_D12 gene in n₂ NSM plants resulted in a significant reduction in epidermal cell prominence and lint fiber production. GhMML4_D12 is arranged in tandem with GhMML3, another MIXTA gene responsible for fuzz fiber development. These two very closely related MIXTA genes direct fiber initiation production in two specialized cell forms: lint and fuzz fibers. They may control the same metabolic pathways in different cell types. The MIXTAs expanded in Malvaceae during their evolution and produced a Malvaceae-specific family that regulates epidermal cell differentiation, different from the gene family that regulates leaf hair trichome development. Cotton has developed a unique transcriptional regulatory network for fiber development. Characterization of target genes regulating fiber production has provided insights into the molecular mechanisms underlying cotton fiber development and has allowed the use of genetic engineering to increase lint yield by inducing more epidermal cells to develop into lint rather than fuzz fibers.

Identification and expression analysis under abiotic stress of the R2R3-MYB genes in Ginkgo biloba L

The R2R3-MYB gene family is the largest MYB subfamily in plants and is involved in the regulation of plant secondary metabolism and specific morphogenesis, as well as the response to biotic and abiotic stress. However, a systematic identification and characterization of this gene family has not been carried out in Ginkgo biloba. In this study, we performed a transcriptome-wide survey from four tissues of G. biloba to determine the genetic variation and expression pattern of the R2R3-MYB genes. We analyzed 45 GbMYBs and identified 42 with a complete coding sequence via conserved motif searches. The MYB domain and other motifs in GbMYBs are highly conserved with Arabidopsis thaliana AtMYBs. Phylogenetic analysis of the GbMYBs and AtMYBs categorized the R2R3-MYBs into 26 subgroups, of which 11 subgroups included proteins from both G. biloba and Arabidopsis, and 1 subgroup was specific to G. biloba. Moreover, the GbMYBs expression patterns were analyzed in different tissues and abiotic stress conditions. The results revealed that GbMYBs were differentially expressed in various tissues and following abiotic stresses and phytohormone treatments, indicating their possible roles in biological processes and abiotic stress tolerance and adaptation. Our study demonstrated the functional diversity of the GbMYBs and will provide a foundation for future research into their biological and molecular functions.

Identification of Norway Spruce MYB-bHLH-WDR Transcription Factor Complex Members Linked to Regulation of the Flavonoid Pathway

Transcription factors (TFs) forming MYB-bHLH-WDR complexes are known to regulate the biosynthesis of specialized metabolites in angiosperms through an intricate network. These specialized metabolites participate in a wide range of biological processes including plant growth, development, reproduction as well as in plant immunity. Studying the regulation of their biosynthesis is thus essential. While MYB (TFs) have been previously shown to control specialized metabolism (SM) in gymnosperms, the identity of their partners, in particular bHLH or WDR members, has not yet been revealed. To gain knowledge about MYB-bHLH-WDR transcription factor complexes in gymnosperms and their regulation of SW, we identified two bHLH homologs of AtTT8, six homologs of the MYB transcription factor AtTT2 and one WDR ortholog of AtTTG1 in Norway spruce. We investigated the expression levels of these genes in diverse tissues and upon treatments with various stimuli including methyl-salicylate, methyl-jasmonate, wounding or fungal inoculation. In addition, we also identified protein-protein interactions among different homologs of MYB, bHLH and WDR. Finally, we generated transgenic spruce cell lines overexpressing four of the Norway spruce AtTT2 homologs and observed differential regulation of genes in the flavonoid pathway and flavonoid contents.

Trichomes control flower bud shape by linking together young petals

Trichomes are widespread in plants and develop from surface cells on different tissues(1). They have many forms and functions, from defensive spines to physical barriers that trap layers of air to insulate against desiccation, but there is growing evidence that trichomes can also have developmental roles in regulating flower structure(2,3). We report here that the trichomes on petals of cotton, Gossypium hirsutum L., are essential for correct flower bud shape through a mechanical entanglement of the trichomes on adjacent petals that anchor the edges to counter the opposing force generated by asymmetric expansion of overlapping petals. Silencing a master regulator of petal trichomes, GhMYB-MIXTA-Like10 (GhMYBML10), by RNA interference (RNAi) suppressed petal trichome growth and resulted in flower buds forming into abnormal corkscrew shapes that exposed developing anthers and stigmas to desiccation damage. Artificially gluing petal edges together could partially restore correct bud shape and fertility. Such petal 'Velcro' is present in other Malvaceae and perhaps more broadly in other plant families, although it is not ubiquitous. This mechanism for physical association between separate organs to regulate flower shape and function is different from the usual organ shape control(4) exerted through cell-to-cell communication and differential cell expansion within floral tissues(5,6).

Small interfering RNAs from bidirectional transcripts of GhMML3_A12 regulate cotton fiber development

Natural antisense transcripts (NATs) are commonly observed in eukaryotic genomes, but only a limited number of such genes have been identified as being involved in gene regulation in plants. In this research, we investigated the function of small RNA derived from a NAT in fiber cell development. Using a map-based cloning strategy for the first time in tetraploid cotton, we cloned a naked seed mutant gene (N1 ) encoding a MYBMIXTA-like transcription factor 3 (MML3)/GhMYB25-like in chromosome A12, GhMML3_A12, that is associated with fuzz fiber development. The extremely low expression of GhMML3_A12 in N1 is associated with NAT production, driven by its 3' antisense promoter, as indicated by the promoter-driven histochemical staining assay. In addition, small RNA deep sequencing analysis suggested that the bidirectional transcriptions of GhMML3_A12 form double-stranded RNAs and generate 21-22 nt small RNAs. Therefore, in a fiber-specific manner, small RNA derived from the GhMML3_A12 locus can mediate GhMML3_A12 mRNA self-cleavage and result in the production of naked seeds followed by lint fiber inhibition in N1 plants. The present research reports the first observation of gene-mediated NATs and siRNA directly controlling fiber development in cotton.

AtRAV1 and AtRAV2 overexpression in cotton increases fiber length differentially under drought stress and delays flowering

There is a longstanding problem of an inverse relationship between cotton fiber qualities versus high yields. To better understand drought stress signaling and adaptation in cotton (Gossypium hirsutum) fiber development, we expressed the Arabidopsis transcription factors RELATED_TO_ABA-INSENSITIVE3/VIVIPAROUS1/(RAV1) and AtRAV2, which encode APETALA2-Basic3 domain proteins shown to repress transcription of FLOWERING_LOCUS_T (FT) and to promote stomatal opening cell-autonomously. In three years of field trials, we show that AtRAV1 and AtRAV2-overexpressing cotton had ∼5% significantly longer fibers with only marginal decreases in yields under well-watered or drought stress conditions that resulted in 40-60% yield penalties and 3-7% fiber length penalties in control plants. The longer transgenic fibers from drought-stressed transgenics could be spun into yarn which was measurably stronger and more uniform than that from well-watered control fibers. The transgenic AtRAV1 and AtRAV2 lines flowered later and retained bolls at higher nodes, which correlated with repression of endogenous GhFT-Like (FTL) transcript accumulation. Elevated expression early in development of ovules was observed for GhRAV2L, GhMYB25-Like (MYB25L) involved in fiber initiation, and GhMYB2 and GhMYB25 involved in fiber elongation. Altered expression of RAVs controlling critical nodes in developmental and environmental signaling hierarchies has the potential for phenotypic modification of crops.

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.