Codon usage bias (CUB) is an important evolutionary feature in a genome which provides important information for studying organism evolution, gene function and exogenous gene expression. The CUB and its shaping factors in the nuclear genomes of four sequenced cotton species, G. arboreum (A2), G. raimondii (D5), G. hirsutum (AD1) and G. barbadense (AD2) were analyzed in the present study. The effective number of codons (ENC) analysis showed the CUB was weak in these four species and the four subgenomes of the two tetraploids. Codon composition analysis revealed these four species preferred to use pyrimidine-rich codons more frequently than purine-rich codons. Correlation analysis indicated that the base content at the third position of codons affect the degree of codon preference. PR2-bias plot and ENC-plot analyses revealed that the CUB patterns in these genomes and subgenomes were influenced by combined effects of translational selection, directional mutation and other factors. The translational selection (P2) analysis results, together with the non-significant correlation between GC12 and GC3, further revealed that translational selection played the dominant role over mutation pressure in the codon usage bias. Through relative synonymous codon usage (RSCU) analysis, we detected 25 high frequency codons preferred to end with T or A, and 31 low frequency codons inclined to end with C or G in these four species and four subgenomes. Finally, 19 to 26 optimal codons with 19 common ones were determined for each species and subgenomes, which preferred to end with A or T. We concluded that the codon usage bias was weak and the translation selection was the main shaping factor in nuclear genes of these four cotton genomes and four subgenomes.

Background Salinity is a major abiotic stress seriously hindering crop yield. Development and utilization of tolerant varieties is the most economical way to address soil salinity. Upland cotton is a major fiber crop and pioneer plant on saline soil and thus its genetic architecture underlying salt tolerance should be extensively explored. Results In this study, genome-wide association analysis and RNA sequencing were employed to detect salt-tolerant qualitative-trait loci (QTLs) and candidate genes in 196 upland cotton genotypes at the germination stage. Using comprehensive evaluation values of salt tolerance in four environments, we identified 33 significant single-nucleotide polymorphisms (SNPs), including 17 and 7 SNPs under at least two and four environments, respectively. The 17 stable SNPs were located within or near 98 candidate genes in 13 QTLs, including 35 genes that were functionally annotated to be involved in salt stress responses. RNA-seq analysis indicated that among the 98 candidate genes, 13 were stably differentially expressed. Furthermore, 12 of the 13 candidate genes were verified by qRT-PCR. RNA-seq analysis detected 6640, 3878, and 6462 differentially expressed genes at three sampling time points, of which 869 were shared. Conclusions These results, including the elite cotton accessions with accurate salt tolerance evaluation, the significant SNP markers, the candidate genes, and the salt-tolerant pathways, could improve our understanding of the molecular regulatory mechanisms under salt stress tolerance and genetic manipulation for cotton improvement.

Cotton (Gossypium spp.) is a leading natural fiber crop and an important source of vegetable protein and oil for humans and livestock. To investigate the genetic architecture of seed nutrients in upland cotton, a genome-wide association study (GWAS) was conducted in a panel of 196 germplasm resources under three environments using a CottonSNP80K chip of 77,774 loci. Relatively high genetic diversity (average gene diversity being 0.331) and phenotypic variation (coefficient of variation, CV, exceeding 3.9%) were detected in this panel. Correlation analysis revealed that the well-documented negative association between seed protein (PR) and oil may be to some extent attributable to the negative correlation between oleic acid (OA) and PR. Linkage disequilibrium (LD) was unevenly distributed among chromosomes and subgenomes. It ranged from 0.10–0.20 Mb (Chr19) to 5.65–5.75 Mb (Chr25) among the chromosomes and the range of Dt-subgenomes LD decay distances was smaller than At-subgenomes. This panel was divided into two subpopulations based on the information of 41,815 polymorphic single-nucleotide polymorphism (SNP) markers. The mixed linear model considering both Q-matrix and K-matrix [MLM(Q+K)] was employed to estimate the association between the SNP markers and the seed nutrients, considering the false positives caused by population structure and the kinship. A total of 47 SNP markers and 28 candidate quantitative trait loci (QTLs) regions were found to be significantly associated with seven cottonseed nutrients, including protein, total fatty acid, and five main fatty acid compositions. In addition, the candidate genes in these regions were analyzed, which included three genes, Gh_D12G1161, Gh_D12G1162, and Gh_D12G1165 that were most likely involved in the control of cottonseed protein concentration. These results improved our understanding of the genetic control of cottonseed nutrients and provided potential molecular tools to develop cultivars with high protein and improved fatty acid compositions in cotton breeding programs through marker-assisted selection.

No abstract

The HVA22 family of genes, induced by abscisic acid and stress, encodes a class of stress response proteins with a conserved TB2/DP1/HVA22 domain that are unique among eukaryotes. Previous studies have shown that HVA22s play an important role in plant responses to abiotic stresses. In the present study, 34, 32, 16, and 17 HVA22s were identified in G. barbadense, G. hirsutum, G. arboreum, and G. raimondii, respectively. These HVA22 genes were classified into nine subgroups, randomly distributed on the chromosomes. Synteny analysis showed that the amplification of the HVA22s were mainly due to segmental duplication or whole genome replication (WGD). Most HVA22s promoter sequences contain a large number of drought response elements (MYB), defense and stress response elements (TC-rich repeats), and hormone response elements (ABRE, ERE, SARE, etc.), suggesting that HVA22s may respond to adversity stresses. Expression profiling demonstrated that most GhHVA22s showed a constitutive expression pattern in G. hirsutum and could respond to abiotic stresses such as salt, drought, and low temperature. Overexpression of GhHVA22E1D (GH_D07G0564) in Arabidopsis thaliana enhances salt and drought tolerance in Arabidopsis. Virus-induced gene silencing of GhHVA22E1D reduced salt and drought tolerance in cotton. This indicates that GhHVA22E1D plays an active role in the plant response to salt stress and drought stress. GhHVA22E1D may act in plant response to adversity by altering the antioxidant capacity of plants. This study provides valuable information for the functional genomic study of the HVA22 gene family in cotton. It also provides a reference for further elucidation of the functional studies of HVA22 in plant resistance to abiotic stress response.