Characterization and Promoter Analysis of a Cotton RING-Type Ubiquitin Ligase (E3) Gene

Genomic and GWAS analyses demonstrate phylogenomic relationships of Gossypium barbadense in China and selection for fibre length, lint percentage and Fusarium wilt resistance

Sea Island cotton (Gossypium barbadense) is the source of the world's finest fibre quality cotton, yet relatively little is understood about genetic variations among diverse germplasms, genes underlying important traits and the effects of pedigree selection. Here, we resequenced 336 G. barbadense accessions and identified 16 million SNPs. Phylogenetic and population structure analyses revealed two major gene pools and a third admixed subgroup derived from geographical dissemination and interbreeding. We conducted a genome-wide association study (GWAS) of 15 traits including fibre quality, yield, disease resistance, maturity and plant architecture. The highest number of associated loci was for fibre quality, followed by disease resistance and yield. Using gene expression analyses and VIGS transgenic experiments, we confirmed the roles of five candidate genes regulating four key traits, that is disease resistance, fibre length, fibre strength and lint percentage. Geographical and temporal considerations demonstrated selection for the superior fibre quality (fibre length and fibre strength), and high lint percentage in improving G. barbadense in China. Pedigree selection breeding increased Fusarium wilt disease resistance and separately improved fibre quality and yield. Our work provides a foundation for understanding genomic variation and selective breeding of Sea Island cotton.

Map-based cloning of qBWT-c12 discovered brassinosteroid-mediated control of organ size in cotton

Assuring fiber yield stability is the primary objective for cotton breeders since the world population is on the rise, and the demand for cotton fiber is increasing every year. Thus, enhancing average cotton boll weight (BWT) could improve seed cotton production, and ultimately to increase cotton fiber yield. This study accomplished the map-based cloning of a novel boll weight regulating locus, qBWT-c12, in cotton. Bulk segregation analysis detected linked markers, aided in the detection of a stable BWT regulating locus, qBWT-c12, on Chr12 in a novel boll size mutant, BS41. Progeny evaluation confined the qBWT-c12 to a 0.89 cM interval between the AD-A12_07 and AD-FM_44 markers in recombinant derived F₃ and F₄ populations. Homology mapping detected a 40 bp insertion-deletion (InDel) site in the AD-FM_44 clone sequence situated +341 downstream of GhBRH1_A12, which showed complete linkage to the BWT phenotype. The suppressed expression of GhBRH1_A12 suggested its putative involvement during early boll development events in BS41. Although brassinosteroid (BR) biosynthesis and signaling pathway genes were up regulated in different tissues, but the organ growth was suppressed leading to dwarf plants, smaller leaves, and de-morphed smaller bolls in BS41. Thus, a disruption in the BR signal cascade is anticipated and could be related to lower GhBRH1_A12 expression in BS41.This study firstly reported the genetic dissection of boll size regulation of G. barbadense in G. hirsutum background using map-based cloning of a BWT regulating locus, qBWT-c12. Moreover, it also emphasized the putative role GhBRH1_A12 in regulating BR homeostasis and its potential to modulate plant growth and boll development in cotton.

Long non-coding RNAs and their potential functions in Ligon-lintless-1 mutant cotton during fiber development

Background

Long non-coding RNAs (LncRNAs) are part of genes, which are not translated into proteins and play a vital role in plant growth and development. Nevertheless, the presence of LncRNAs and how they functions in Ligon-lintless-1 mutant during the early cessation of cotton fiber development are still not well understood. In order to investigate the function of LncRNAs in cotton fiber development, it is necessary and important to identify LncRNAs and their potential roles in fiber cell development.

Results

In this work, we identified 18,333 LncRNAs, with the proportion of long intergenic noncoding RNAs (LincRNAs) (91.5%) and anti-sense LncRNAs (8.5%), all transcribed from Ligon-lintless-1 (Li1) and wild-type (WT). Expression differences were detected between Ligon-lintless-1 and wild-type at 0 and 8 DPA (day post anthesis). Pathway analysis and Gene Ontology based on differentially expressed LncRNAs on target genes, indicated fatty acid biosynthesis and fatty acid elongation being integral to lack of fiber in mutant cotton. The result of RNA-seq and RT-qPCR clearly singles out two potential LncRNAs, LNC_001237 and LNC_017085, to be highly down-regulated in the mutant cotton. The two LncRNAs were found to be destabilized or repressed by ghr-miR2950. Both RNA-seq analysis and RT-qPCR results in Ligon-lintless-1 mutant and wild-type may provide strong evidence of LNC_001237, LNC_017085 and ghr-miR2950 being integral molecular elements participating in various pathways of cotton fiber development.

Conclusion

The results of this study provide fundamental evidence for the better understanding of LncRNAs regulatory role in the molecular pathways governing cotton fiber development. Further research on designing and transforming LncRNAs will help not only in the understanding of their functions but will also in the improvement of fiber quality.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5978-5) contains supplementary material, which is available to authorized users.

GhHUB2, a ubiquitin ligase, is involved in cotton fiber development via the ubiquitin-26S proteasome pathway

Cotton fibers, which are extremely elongated single cells of epidermal seed trichomes and have highly thickened cell walls, constitute the most important natural textile material worldwide. However, the regulation of fiber development is not well understood. Here, we report that GhHUB2, a functional homolog of AtHUB2, controls fiber elongation and secondary cell wall (SCW) deposition. GhHUB2 is ubiquitously expressed, including within fibers. Overexpression of GhHUB2 in cotton increased fiber length and SCW thickness, while RNAi knockdown of GhHUB2 resulted in shortened fibers and thinner cell walls. We found that GhHUB2 interacted with GhKNL1, a transcriptional repressor predominantly expressed in developing fibers, and that GhHUB2 ubiquitinated and degraded GhKNL1 via the ubiquitin-26S proteasome pathway. GhHUB2 negatively regulated GhKNL1 protein levels and lead to the disinhibition of genes such as GhXTH1, Gh1,3-β-G, GhCesA4, GhAGP4, GhCTL1, and GhCOBL4, thus promoting fiber elongation and enhancing SCW biosynthesis. We found that GhREV-08, a transcription factor that participates in SCW deposition and auxin signaling pathway, was a direct target of GhKNL1. In conclusion, our study uncovers a novel function of HUB2 in plants in addition to its monoubiquitination of H2B. Moreover, we provide evidence for control of the fiber development by the ubiquitin-26S proteasome pathway.

GhMYB1 regulates SCW stage-specific expression of the GhGDSL promoter in the fibres of Gossypium hirsutum L

Secondary cell wall (SCW) biosynthesis is an important stage of the cotton fibre development, and its transcriptional regulation is poorly understood. We selected the Gossypium hirsutum GDSL (GhGDSL) lipase/hydrolase gene (CotAD_74480), which is expressed during SCW biosynthesis (19 through to 25 days postanthesis; DPA), for study. T₁ -transgenic cotton lines expressing the β-glucuronidase (gus) reporter under the control of a 1026-bp promoter fragment of GhGDSL (P_GhGDSL ) showed 19 DPA stage-specific increase in GUS expression. 5' deletion indicated that the 194-bp fragment between -788 and -594 relative to the transcription start site was essential for this stage-specific expression. Site-directed mutagenesis of eight transcription factor binding sites within P_GhGDSL demonstrated that the MYB1AT motif (AAACCA) at -603/-598 was critical for the 19 DPA-specific reporter gene expressions. Yeast one-hybrid (Y1H) analysis identified nine proteins, including GhMYB1 (CotAD_64719) that bound to the P_GhGDSL promoter. Further, Y1H experiments using the 5' promoter deletions and individually mutated promoter motifs indicated that GhMYB1 interacted with P_GhGDSL at MYB1AT sequence. GhMYB1 was expressed specifically in fibre from 19 DPA, overlapping with the sharp rise in GhGDSL expression, indicating that it could regulate GhGDSL during fibre development. Analysis of genes co-expressed with GhMYB1 showed that it potentially regulates a number of other 19-25 DPA-specific genes in networks including those functioning in the cell wall and precursor synthesis, but not the major polysaccharide and protein components of the fibre SCW. GhGDSL and its promoter are therefore potential tools for the improvement of cotton fibre quality traits.

Genome-Wide Survey Reveals Transcriptional Differences Underlying the Contrasting Trichome Phenotypes of Two Sister Desert Poplars

Trichomes, which are widely used as an important diagnostic characteristic in plant species delimitation, play important roles in plant defense and adaptation to adverse environments. In this study, we used two sister poplar species, Populus pruinosa and Populus euphratica—which have, respectively, dense and sparse trichomes—to examine the genetic differences associated with these contrasting phenotypes. The results showed that 42 and 45 genes could be identified as candidate genes related to trichomes in P. pruinosa and P. euphratica, respectively; most of these genes possessed high degrees of diversification in their coding sequences, but they were similar in intron/exon structure in the two species. We also found that most of the candidate trichome genes were expressed at higher levels in P. pruinosa, which has dense trichomes, than in P. euphratica, where there are few trichomes. Based on analyses of transcriptional profiles, a total of 195 genes, including many transcription factors, were found to show distinct differences in expression. The results of gene function annotation suggested that the genes identified as having contrasting levels of expression level are mainly associated with trichome elongation, ATPase activity, and hormone transduction. Changes in the expression of these and other related genes with high sequence diversification may have contributed to the contrast in the pattern of trichome phenotypes between the two species.

Overexpression of GhSARP1 encoding a E3 ligase from cotton reduce the tolerance to salt in transgenic Arabidopsis

Ubiquitination plays a very important role in the response to abiotic stresses of plant. To identify key regulators of salt stress, a gene GhSARP1(Salt-Associated Ring finger Protein)encoding C3H2C3-type E3 ligase, was cloned from cotton. Transcription level of GhSARP1 was high in leaf, flower and fiber of 24,27 and 27DPA (Days Post-Anthesis), but low in root and stem. Except PEG6000 treatment, the expression of GhSARP1 was down-regulated by NaCl, cold and ABA after being treated for 1 h. GhSARP1-GFP fusion protein located on the plasma membrane, which was dependent on trans-membrane motif. In vitro ubiquitination assay showed that GhSARP1 had E3 ligase activity. Heterogeneous overexpression of GhSARP1reduced salt tolerance of transgenic Arabidopsis in germination and post-germination stage. Our results suggested that the GhSARP1 might negatively regulate the response to salt stress mediated by the ubiquitination in cotton.

Detection and validation of one stable fiber strength QTL on c9 in tetraploid cotton

Fiber strength is an essential trait of fiber property in cotton, and it is quantitatively inherited. Identification of stable quantitative trait loci (QTL) contributing to fiber strength would provide the key basis for marker-assisted selection (MAS) in cotton breeding. In this study, four interspecific hybridization populations were established with a common G. barbadense parent Pima 90-53 and two G. hirsutum parents (CCRI 8 and Handan 208), each of which had fiber strength characteristic. Based on the phenotypic data of fiber strength from seven environments, a stable QTL, qFS-c9-1, was detected and validated on c9 in a marker interval between SSR markers NAU2395 and NAU1092. The QTL explaining 14.4-17.9 % of the phenotypic variation was firstly detected in two populations (CCRI 8 × Pima 90-53, BC1F1 and BC1F2) and its derived lines in four environments. And it accounting for 12.1-14.8 % of the phenotypic variation was further confirmed in two populations (Handan 208 × Pima 90-53, BC1F1, and F2) under one environment. In silico mapping using three sequenced cotton genomes indicated that homologous genes, anchored by NAU2395 and NAU1092, were aligned to the G. arboreum genome within a physical distance between 81.10 Mbps and 87.07 Mbps. In that interval, several genes were confirmed in literatures to associate with fiber development. Among these genes, seven genes were further selected for an expression analysis through fiber development transcriptome database, revealing unique expression patterns across different stages of fiber development between CCRI 8 and Pima 90-53. The genes underlying qFS-c9-1 were favorable to fine mapping and cloning. The current study results provided valuable evidence for mapping stable QTL of fiber strength utilizing multiple populations and environments, as well as map-based cloning the candidate gene underlying the QTL for future prospective research directions.

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.