QTL and candidate gene identification of the node of the first fruiting branch (NFFB) by QTL-seq in upland cotton (Gossypium hirsutum L.).

BACKGROUND: The node of the first fruiting branch (NFFB) is an important precocious trait in cotton. Many studies have been conducted on the localization of quantitative trait loci (QTLs) and genes related to fiber quality and yield, but there has been little attention to traits related to early maturity, especially the NFFB, in cotton. RESULTS: To identify the QTL associated with the NFFB in cotton, a BC4F2 population comprising 278 individual plants was constructed. The parents and two DNA bulks for high and low NFFB were whole genome sequenced, and 243.8 Gb of clean nucleotide data were generated. A total of 449,302 polymorphic SNPs and 135,353 Indels between two bulks were identified for QTL-seq. Seventeen QTLs were detected and localized on 11 chromosomes in the cotton genome, among which two QTLs (qNFFB-Dt2-1 and qNFFB-Dt3-3) were located in hotspots. Two candidate genes (GhAPL and GhHDA5) related to the NFFB were identified using quantitative real-time PCR (qRT-PCR) and virus-induced gene silencing (VIGS) experiments in this study. Both genes exhibited higher expression levels in the early-maturing cotton material RIL182 during flower bud differentiation, and the silencing of GhAPL and GhHDA5 delayed the flowering time and increased the NFFB compared to those of VA plants in cotton. CONCLUSIONS: Our study preliminarily found that GhAPL and GhHDA5 are related to the early maturity in cotton. The findings provide a basis for the further functional verification of candidate genes related to the NFFB and contribute to the study of early maturity in cotton.

High-resolution temporal dynamic transcriptome landscape reveals a GhCAL-mediated flowering regulatory pathway in cotton (Gossypium hirsutum L.).

The transition from vegetative to reproductive growth is very important for early maturity in cotton. However, the genetic control of this highly dynamic and complex developmental process remains unclear. A high-resolution tissue- and stage-specific transcriptome profile was generated from six developmental stages using 72 samples of two early-maturing and two late-maturing cotton varieties. The results of histological analysis of paraffin sections showed that flower bud differentiation occurred at the third true leaf stage (3TLS) in early-maturing varieties, but at the fifth true leaf stage (5TLS) in late-maturing varieties. Using pairwise comparison and weighted gene co-expression network analysis, 5312 differentially expressed genes were obtained, which were divided into 10 gene co-expression modules. In the MElightcyan module, 46 candidate genes regulating cotton flower bud differentiation were identified and expressed at the flower bud differentiation stage. A novel key regulatory gene related to flower bud differentiation, GhCAL, was identified in the MElightcyan module. Anti-GhCAL transgenic cotton plants exhibited late flower bud differentiation and flowering time. GhCAL formed heterodimers with GhAP1-A04/GhAGL6-D09 and regulated the expression of GhAP1-A04 and GhAGL6-D09. GhAP1-A04- and GhAGL6-D09-silenced plants also showed significant late flowering. Finally, we propose a new flowering regulatory pathway mediated by GhCAL. This study elucidated the molecular mechanism of cotton flowering regulation and provides good genetic resources for cotton early-maturing breeding.

Genome-wide identification and characterization of multiple C2 domains and transmembrane region proteins in Gossypium hirsutum.

BACKGROUND: Multiple C2 domains and transmembrane region proteins (MCTPs) may act as transport mediators of other regulators. Although increased number of MCTPs in higher plants implies their diverse and specific functions in plant growth and development, only a few plant MCTPs have been studied and no study on the MCTPs in cotton has been reported. RESULTS: In this study, we identified 31 MCTPs in G. hirsutum, which were classified into five subfamilies according to the phylogenetic analysis. GhMCTPs from subfamily V exhibited isoelectric points (pIs) less than 7, whereas GhMCTPs from subfamily I, II, III and IV exhibited pIs more than 7.5, implying their distinct biological functions. In addition, GhMCTPs within subfamily III, IV and V exhibited more diverse physicochemical properties, domain architectures and expression patterns than GhMCTPs within subfamily I and II, suggesting that GhMCTPs within subfamily III, IV and V diverged to perform more diverse and specific functions. Analyses of conserved motifs and pIs indicated that the N-terminus was more divergent than the C-terminus and GhMCTPs' functional divergence might be mainly contributed by the N-terminus. Furthermore, yeast two-hybrid assay indicated that the N-terminus was responsible to interact with target proteins. Phylogenetic analysis classified multiple N-terminal C2 domains into four subclades, suggesting that these C2 domains performed different molecular functions in mediating the transport of target proteins. CONCLUSIONS: Our systematic characterization of MCTPs in G. hirsutum will provide helpful information to further research GhMCTPs' molecular roles in mediating other regulators' transport to coordinate growth and development of various cotton tissues.

Genome-wide identification and expression patterns analysis of the RPD3/HDA1 gene family in cotton.

BACKGROUND: Histone deacetylases (HDACs) catalyze histone deacetylation and suppress gene transcription during various cellular processes. Within the superfamily of HDACs, RPD3/HDA1-type HDACs are the most studied, and it is reported that RPD3 genes play crucial roles in plant growth and physiological processes. However, there is a lack of systematic research on the RPD3/HDA1 gene family in cotton. RESULTS: In this study, genome-wide analysis identified 9, 9, 18, and 18 RPD3 genes in Gossypium raimondii, G. arboreum, G. hirsutum, and G. barbadense, respectively. This gene family was divided into 4 subfamilies through phylogenetic analysis. The exon-intron structure and conserved motif analysis revealed high conservation in each branch of the cotton RPD3 genes. Collinearity analysis indicated that segmental duplication was the primary driving force during the expansion of the RPD3 gene family in cotton. There was at least one presumed cis-element related to plant hormones in the promoter regions of all GhRPD3 genes, especially MeJA- and ABA-responsive elements, which have more members than other hormone-relevant elements. The expression patterns showed that most GhRPD3 genes had relatively high expression levels in floral organs and performed higher expression in early-maturity cotton compared with late-maturity cotton during flower bud differentiation. In addition, the expression of GhRPD3 genes could be significantly induced by one or more abiotic stresses as well as exogenous application of MeJA or ABA. CONCLUSIONS: Our findings reveal that GhRPD3 genes may be involved in flower bud differentiation and resistance to abiotic stresses, which provides a basis for further functional verification of GhRPD3 genes in cotton development and a foundation for breeding better early-maturity cotton cultivars in the future.

A comprehensive analysis of cotton VQ gene superfamily reveals their potential and extensive roles in regulating cotton abiotic stress.

BACKGROUND: Valine-glutamine (VQ) motif-containing proteins play important roles in plant growth, development and abiotic stress response. For many plant species, the VQ genes have been identified and their functions have been described. However, little is known about the origin, evolution, and functions (and underlying mechanisms) of the VQ family genes in cotton. RESULTS: In this study, we comprehensively analyzed the characteristics of 268 VQ genes from four Gossypium genomes and found that the VQ proteins evolved into 10 clades, and each clade had a similar structural and conservative motif. The expansion of the VQ gene was mainly through segmental duplication, followed by dispersal. Expression analysis revealed that many GhVQs might play important roles in response to salt and drought stress, and GhVQ18 and GhVQ84 were highly expressed under PEG and salt stress. Further analysis showed that GhVQs were co-expressed with GhWRKY transcription factors (TFs), and microRNAs (miRNAs) could hybridize to their cis-regulatory elements. CONCLUSIONS: The results in this study broaden our understanding of the VQ gene family in plants, and the analysis of the structure, conserved elements, and expression patterns of the VQs provide a solid foundation for exploring their specific functions in cotton responding to abiotic stresses. Our study provides significant insight into the potential functions of VQ genes in cotton.

High-density linkage map construction and QTL analysis for earliness-related traits in Gossypium hirsutum L.

BACKGROUND: Gossypium hirsutum L., or upland cotton, is an important renewable resource for textile fiber. To enhance understanding of the genetic basis of cotton earliness, we constructed an intra-specific recombinant inbred line population (RIL) containing 137 lines, and performed linkage map construction and quantitative trait locus (QTL) mapping. RESULTS: Using restriction-site associated DNA sequencing, a genetic map composed of 6,434 loci, including 6,295 single nucleotide polymorphisms and 139 simple sequence repeat loci, was developed from RIL population. This map spanned 4,071.98 cM, with an average distance of 0.63 cM between adjacent markers. A total of 247 QTLs for six earliness-related traits were detected in 6 consecutive years. In addition, 55 QTL coincidence regions representing more than 60 % of total QTLs were found on 22 chromosomes, which indicated that several earliness-related traits might be simultaneously improved. Fine-mapping of a 2-Mb region on chromosome D3 associated with five stable QTLs between Marker25958 and Marker25963 revealed that lines containing alleles derived from CCRI36 in this region exhibited smaller phenotypes and earlier maturity. One candidate gene (EMF2) was predicted and validated by quantitative real-time PCR in early-, medium- and late-maturing cultivars from 3- to 6-leaf stages, with highest expression level in early-maturing cultivar, CCRI74, lowest expression level in late-maturing cultivar, Bomian1. CONCLUSIONS: We developed an SNP-based genetic map, and this map is the first high-density genetic map for short-season cotton and has the potential to provide deeper insights into earliness. Cotton earliness-related QTLs and QTL coincidence regions will provide useful materials for QTL fine mapping, gene positional cloning and MAS. And the gene, EMF2, is promising for further study.

Gh4CL20/20A involved in flavonoid biosynthesis is essential for male fertility in cotton (Gossypium hirsutum L.).

Flavonoids have been shown to play an essential role in plant growth and fertility. 4-Coumarate CoA ligase (4CL) is one of the indispensable enzymes involved in the biosynthesis of flavonoids. However, the role of 4CL and flavonoids in impact on cotton fertility is still unknown. In this study, on the basis of identification of an additional Gh4CL gene, Gh4CL20A, by using an updated G. hirsutum genome, we found that Gh4CL20A and its homologous Gh4CL20 were preferentially expressed in petals and stamens. The petals of the loss-of-function Gh4CL20/Gh4CL20A mutant generated by CRISPR/Cas9 gene editing remained white until wilting. Notably, the mutant showed indehiscent anthers, reduced number of pollen grains and pollen viability, leading to male sterility. Histological analysis revealed that abnormal degradation of anther tapetum at the tetrad stage and abnormal pollen grain development at the mature stage caused male sterility of the gene editing mutant. Analysis of the anther transcriptome identified a total of 10574 and 11962 genes up- and down-regulated in the mutant, respectively, compared to the wild-type. GO, KEGG, and WGCNA analyses linked the abnormality of the mutant anthers to the defective flavonoid biosynthetic pathway, leading to decreased activity of 4CL and chalcone isomerase (CHI) and reduced accumulation of flavonoids in the mutant. These results imply a role of Gh4CL20/Gh4CL20A in assuring proper development of cotton anthers by regulating flavonoid metabolism. This study elucidates a molecular mechanism underlying cotton anther development and provides candidate genes for creating cotton male sterile germplasm that has the potential to be used in production of hybrid seeds.

The MADS transcription factor GhFYF is involved in abiotic stress responses in upland cotton (Gossypium hirsutum L.).

Cotton is an important textile industry raw material crops, which plays a critical role in the development of society. MADS transcription factors (TFs) play a key role about the flowering time, flower development, and abiotic stress responses in plants, but little is known about their functions on abiotic stress in cotton. In this study, a MIKCC subfamily gene from cotton, GhFYF (FOREVER YOUNG FLOWER), was isolated and characterized. Our data showed that GhFYF localized to the nucleus. A ß-glucuronidase (GUS) activity assay revealed that the promoter of GhFYF was mainly expressed in the flower and seed of ProGhFYF::GUS transgenic A. thaliana plants. The GUS staining of flowers and seeds was deepened after drought, salt treatment, and the expression level of the GUS gene and corresponding stress genes AtERD10, AtAnnexin1 are up-regulated in the inflorescence. Overexpression GhFYF in A. thaliana could promote the seed germination and growth under different salt concentrations, and determin the proline content. Yeast two-hybrid (Y2H) assays showed that GhFYF interacted with the HAD-like protein GhGPP2, which has responds to abiotic stress. Our findings indicate that GhFYF is involved in abiotic stress responses, especially for salt stress. This work establishes a solid foundation for further functional analysis of the GhFYF gene in cotton.

3D biomaterial P scaffolds carrying umbilical cord mesenchymal stem cells improve biointegration of keratoprosthesis.

Biointegration of a keratoprosthesis (KPro) is critical for the device stability and long-term retention. Biointegration of the KPro device and host tissue takes place between the surrounding corneal graft and the central optic (made by poly (methyl methacrylate)). Our previous clinical results showed that auricular cartilage reinforcement is able to enhance the KPro biointegration. However, the auricular cartilage is non-renewable and difficult to acquire. In this study, we developed a novel type of biomaterial using a three-dimensional porous polyethylene glycol acrylate scaffold (3D biological P-scaffold) carrier with chondrocytes differentiated from induced human umbilical cord mesenchymal stem cells (hUC-MSCs) and tested in rabbit corneas. The results showed hUC-MSCs bear stem cell properties and coule be induced into chondrocytes, P-scaffold is beneficial to the growth and differentiation of hUC-MSCs bothin vivoandin vitro. Besides, after implanting the P-scaffold into the corneal stroma, no serious immune rejection response, such as corneal ulcer or perforation were seen, suggested a good biocompatibility of P-scaffold with the corneal tissue. Moreover, after implanting P-scaffold in together with the differentiated chondrocytes into the rabbit corneal stroma, they significantly increased corneal thickness and strengthened the host cornea, and chondrocytes could stably persist inside the cornea. In summary, the 3D biological P-scaffold carrying differentiated hUC-MSCs could be the preferable material for KPro reinforcement.

Silencing of GhKEA4 and GhKEA12 Revealed Their Potential Functions Under Salt and Potassium Stresses in Upland Cotton.

The K+ efflux antiporter (KEA) mediates intracellular K+ and H+ homeostasis to improve salt tolerance in plants. However, the knowledge of KEA gene family in cotton is largely absent. In the present study, 8, 8, 15, and 16 putative KEA genes were identified in Gossypium arboreum, G. raimondii, G. hirsutum, and G. barbadense, respectively. These KEA genes were classified into three subfamilies, and members from the same subfamilies showed similar motif compositions and gene structure characteristics. Some hormone response elements and stress response elements were identified in the upstream 2000 bp sequence of GhKEAs. Transcriptome data showed that most of the GhKEAs were highly expressed in roots and stems. The quantificational real-time polymerase chain reaction (qRT-PCR) results showed that most of the GhKEAs responded to low potassium, salt and drought stresses. Virus-induced gene silencing (VIGS) experiments demonstrated that under salt stress, after silencing genes GhKEA4 and GhKEA12, the chlorophyll content, proline content, soluble sugar content, peroxidase (POD) activity and catalase (CAT) activity were significantly decreased, and the Na+/K+ ratio was extremely significantly increased in leaves, leading to greater salt sensitivity. Under high potassium stress, cotton plants silenced for the GhKEA4 could still maintain a more stable Na+ and K+ balance, and the activity of transporting potassium ions from roots into leaves was reduced silenced for GhKEA12. Under low potassium stress, silencing the GhKEA4 increased the activity of transporting potassium ions to shoots, and silencing the GhKEA12 increased the ability of absorbing potassium ions, but accumulated more Na+ in leaves. These results provided a basis for further studies on the biological roles of KEA genes in cotton development and adaptation to stress conditions.

The MADS transcription factor GhAP1.7 coordinates the flowering regulatory pathway in upland cotton (Gossypium hirsutum L.).

MADS-box gene family plays an important role in the molecular regulatory network of flower development. APETALA1 (AP1), a MADS-box gene, plays an important role in the development of flower organs. Although many studies about MADS-box family genes have been reported, the function of AP1 is still not clear in cotton. In this study, GhAP1.7 (Gh_D03G0922), a candidate gene for cotton flower time and plant height obtained from our previous studies, was cloned from CCRI50 cotton variety and functionally characterized. Subcellular localization demonstrated that GhAP1.7 was located in nucleus. Infection test of Arabidopsis revealed that GhAP1.7 could cause precocious flowering and virus-induced gene silence (VIGS) assay demonstrated that GhAP1.7 could lead to delayed flowering of cotton plants. Yeast one-hybrid assays and transient dual-luciferase assays suggested that floral meristem identity control gene LEAFY (LFY) can bind the promoter of GhAP1.7 and negatively regulate it. Our research indicated that GhAP1.7 might work as a positive regulator in plant flowering. Moreover, GhAP1.7 may negatively regulated by GhLFY in the regulatory pathways. This work laid the foundation for subsequent functional studies of GhAP1.7.

Pectate lyase-like Gene GhPEL76 regulates organ elongation in Arabidopsis and fiber elongation in cotton.

Pectate lyases (PELs) play important roles in plant growth and development, mainly by degrading the pectin in primary cell walls. However, the role of PELs in cotton fiber elongation, which also involves changes in cellular structure and components, is poorly understood. Therefore, we aimed to isolate and characterize GhPEL76, as we suspected it to contribute to the regulation of fiber elongation. Expression analysis (qRT-PCR) revealed that GhPEL76 is predominately expressed in cotton fiber, with significantly different expression levels in long- and short-fiber cultivars, and that GhPEL76 expression is responsive to gibberellic acid and indoleacetic acid treatment. Furthermore, GhPEL76 promoter-driven ß-glucuronidase activity was detected in the roots, hypocotyls, and leaves of transgenic Arabidopsis plants, and the overexpression of GhPEL76 in transgenic Arabidopsis promoted the elongation of several organs, including petioles, hypocotyls, primary roots, and trichomes. Additionally, the virus-induced silencing of GhPEL76 in cotton reduced fiber length, and both yeast one-hybrid and transient dual-luciferase assays suggested that GhbHLH13, a bHLH transcription factor that is up-regulated during fiber elongation, activates GhPEL76 expression by binding to the G-box of the GhPEL76 promoter region. Therefore, these results suggest GhPEL76 positively regulates fiber elongation and provide a basis for future studies of cotton fiber development.

Genome-Wide Identification and Characterization of Glycosyltransferase Family 47 in Cotton.

The glycosyltransferase (GT) 47 family is involved in the biosynthesis of xylose, pectin and xyloglucan and plays a significant role in maintaining the normal morphology of the plant cell wall. However, the functions of GT47s are less well known in cotton. In the present study, a total of 53, 53, 105 and 109 GT47 genes were detected by genome-wide identification in Gossypium arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. All the GT47s were classified into six major groups via phylogenetic analysis. The exon/intron structure and protein motifs indicated that each branch of the GT47 genes was highly conserved. Collinearity analysis showed that GT47 gene family expansion occurred in Gossypium spp. mainly through whole-genome duplication and that segmental duplication mainly promoted GT47 gene expansion within the A and D subgenomes. The Ka/Ks values suggested that the GT47 gene family has undergone purifying selection during the long-term evolutionary process. Transcriptomic data and qRT-PCR showed that GhGT47 genes exhibited different expression patterns in each tissue and during fiber development. Our results suggest that some genes in the GhGT47 family might be associated with fiber development and the abiotic stress response, which could promote further research involving functional analysis of GT47 genes in cotton.

Synthesis and Electrochemical Performance of Molybdenum Disulfide-Reduced Graphene Oxide-Polyaniline Ternary Composites for Supercapacitors.

Molybdenum disulfide/reduced graphene oxide/polyaniline ternary composites (MoS2/rGO/PANI) were designed and synthesized by a facile two-step approach including hydrothermal and in situ polymerization process. The MoS2/rGO/PANI composites presented an interconnected 3D network architecture, in which PANI uniformly coated the outer surface of the MoS2/rGO binary composite. The MoS2/rGO/PANI composites with a weight percent of 80% (MGP-80) exhibits the best specific capacitance (570 F g-1 at 1 A g-1) and cycling stabilities (78.6% retained capacitance after 500 cycles at 1 A g-1). The excellent electrochemical capacitive performance is attributed to its 3D network structure and the synergistic effects among the three components that make the composites obtain both pseudocapacitance and double-layer capacitance.