To date, no miRNAs have been identified in the important diploid cotton species although there are several reports on miRNAs in upland cotton. In this study, we identified 73 miRNAs, belonging to 49 families, from Asiatic cotton using a well-developed comparative genome-based homologue search. Several of the predicted miRNAs were validated using quantitative real time PCR (qRT-PCR). The length of miRNAs varied from 18 to 22 nt with an average of 20 nt. The length of miRNA precursors also varied from 46 to 684 nt with an average of 138 ±120 nt. For a majority of Asiatic cotton miRNAs, there is only one member per family; however, multiple members were identified for miRNA 156, 414, 837, 838, 1044, 1533, 2902, 2868, 5021 and 5142 families. Nucleotides A and U were dominant, accounted for 62.95%, in the Asiatic cotton pre-miRNAs. The Asiatic cotton pre-miRNAs had high negative minimal folding free energy (MFE) and adjusted MFE (AMFE) and high MFE index (MFEI). Many miRNAs identified in Asiatic cotton suggest that miRNAs also play a similar regulatory mechanism in diploid cotton.

Lint percentage is one of the most essential yield components and an important economic index for cotton planting. Improving lint percentage is an effective way to achieve high-yield in cotton breeding worldwide, especially upland cotton (Gossypium hirsutum L.). However, the genetic basis controlling lint percentage has not yet been systematically understood. Here, we performed a genome-wide association mapping for lint percentage using a natural population consisting of 189 G. hirsutum accessions (188 accessions of G. hirsutum races and one cultivar TM-1). The results showed that 274 single-nucleotide polymorphisms (SNPs) significantly associated with lint percentage were detected, and they were distributed on 24 chromosomes. Forty-five SNPs were detected at least by two models or at least in two environments, and their 5 Mb up- and downstream regions included 584 makers related to lint percentage identified in previous studies. In total, 11 out of 45 SNPs were detected at least in two environments, and their 550 Kb up- and downstream region contained 335 genes. Through RNA sequencing, gene annotation, qRT-PCR, protein–protein interaction analysis, the cis-elements of the promotor region, and related miRNA prediction, Gh_D12G0934 and Gh_A08G0526 were selected as key candidate genes for fiber initiation and elongation, respectively. These excavated SNPs and candidate genes could supplement marker and gene information for deciphering the genetic basis of lint percentage and facilitate high-yield breeding programs of G. hirsutum ultimately.

Sedoheptulose-1,7-bisphosphatase (SBPase, EC 3.1.3.37) is a key enzyme in the plant Calvin cycle and one of the main rate-limiting enzymes in the plant photosynthesis pathway. Many studies have demonstrated that the SBPase gene plays an important role in plant photosynthetic efficiency, yield, and stress responses; however, few studies have been conducted on the function and expression of the GhSBPase gene in upland cotton. In this study, our results showed that the coding sequence (CDS) of GhSBPase gene was 1182 bp, encoding a protein with 393 amino acids. The GhSBPase protein had adenosine monophosphate (AMP) binding site and a FIG (FBPase/IMPase/glpX) domain, and had six Cys residues and a CGGT(A/Q)C motif that were involved in redox regulation in plants. Evolutionarily, the GhSBPase protein clustered into the dicotyledon subgroup and was most closely related to the tomato SlSBPase protein. Western-blot analysis further indicated that the GhSBPase gene was indeed the gene encoding the SBPase protein in upland cotton. The GhSBPase protein was localized in chloroplast, which was consistent with its function as a key enzyme in photosynthesis. The GhSBPase gene was specifically highly expressed in leaves, and its expression level was significantly lower in a yellow-green leaf mutant than in the wild type. Moreover, the GhSBPase expression was in response to drought, salt, high- and low-temperature stress, and exhibits different expression patterns. The GhSBPase promoter had the cis-acting elements in response to abiotic stress, phytohormone, and light. In addition, the GhSBPase expression was positively correlated with the chlorophyll fluorescence parameters, suggesting that changes in the expression of the GhSBPase had potential applicability in breeding for enhanced cotton photosynthetic efficiency. These results will help to understand the function of the GhSBPase gene in photosynthesis and the adaptability of plants to external stress and provide important gene information for the high-yield breeding of crops in the future.

CIPK (calcineurin B-like-interacting protein kinase) is a kind of serine/threonine protein kinase widely existing in plants, and it plays an important role in plant growth and development and stress response. To better understand the biological functions of the GhCIPK23 gene in upland cotton, the coding sequence (CDS) of the GhCIPK23 gene was cloned in upland cotton, and its protein sequence, evolutionary relationship, subcellular localization, expression pattern and cis-acting elements in the promoter region were analyzed. Our results showed that the full-length CDS of GhCIPK23 was 1368 bp, encoding a protein with 455 amino acids. The molecular weight and isoelectric point of this protein were 50.83 KDa and 8.94, respectively. The GhCIPK23 protein contained a conserved N-terminal protein kinase domain and C-terminal regulatory domain of the CIPK gene family member. Phylogenetic tree analysis demonstrated that GhCIPK23 had a close relationship with AtCIPK23, followed by OsCIPK23, and belonged to Group A with AtCIPK23 and OsCIPK23. The subcellular localization experiment indicated that GhCIPK23 was located in the plasma membrane. Tissue expression analysis showed that GhCIPK23 had the highest expression in petals, followed by sepals, and the lowest in fibers. Stress expression analysis showed that the expression of the GhCIPK23 gene was in response to drought, salt, low-temperature and exogenous abscisic acid (ABA) treatment, and had different expression patterns under different stress conditions. Further cis-acting elements analysis showed that the GhCIPK23 promoter region had cis-acting elements in response to abiotic stress, phytohormones and light. These results established a foundation for understanding the function of GhCIPK23 and breeding varieties with high-stress tolerance in cotton.

"Pre-summer boll, summer boll and autumn boll" have long been regarded as an important index for prematurity and high-yield in cotton. In this study, the prematurity and high-yield cotton cultivar, Baimian 2, was used as the central parent to cross separately with the middle-late-matury lines TM-1 and CIR12, and then two populations of F2 and F2:3 family lines were obtained, which was used to construct two genetic linkage maps. These maps were comprised of 269 and 127 marker loci with the total length 1837.8 cM and 1244.3 cM, respectively. Results of QTL location showed that a total of 29 QTLs were detected in the two combinations, including 16 suggestive QTLs and 13 significant QTLs, of which 5 significant QTLs had higher LOD values that was not only greater than 3 but also greater than the threshold calculated by permutation test. The contribution rate of 16 QTLs explained 10.9%-44.5% of the phenotypic variations. Four common QTLs, qPSB-17 for pre-summer boll, qSB-17 (qSB-17a/17b) for summer boll, and qAB-17 and qAB-12/26 for autumn boll, were detected close to common markers of the same chromosome in the two combinations, which could be applied in marker-assisted selection. Moreover, the contribution rate of qSB-17 (qSB-17a/17b) for summer boll in the two combinations was greater than 10%, and that of qAB-17, qAB-12/26 for autumn boll in one combination were greater than 10%. These common QTLs with greater contribution rates should take into consideration firstly in marker-assisted selection.