Upland cotton (
) is the world's leading fiber crop and one of the most important oilseed crops. Genetic improvement of cotton has primarily focused on fiber yield and quality. However, there is an ...increased interest and demand for enhanced cottonseed traits, including protein, oil, fatty acids, and amino acids for broad food, feed and biofuel applications. As a byproduct of cotton production, cottonseed is an important source of edible oil in many countries and could also be a vital source of protein for human consumption. The focus of cotton breeding on high yield and better fiber quality has substantially reduced the natural genetic variation available for effective cottonseed quality improvement within Upland cotton. However, genetic variation in cottonseed oil and protein content exists within the genus of
and cultivated cotton. A plethora of genes and quantitative trait loci (QTLs) (associated with cottonseed oil, fatty acids, protein and amino acids) have been identified, providing important information for genetic improvement of cottonseed quality. Genetic engineering in cotton through RNA interference and insertions of additional genes of other genetic sources, in addition to the more recent development of genome editing technology has achieved considerable progress in altering the relative levels of protein, oil, fatty acid profile, and amino acids composition in cottonseed for enhanced nutritional value and expanded industrial applications. The objective of this review is to summarize and discuss the cottonseed oil biosynthetic pathway and major genes involved, genetic basis of cottonseed oil and protein content, genetic engineering, genome editing through CRISPR/Cas9, and QTLs associated with quantity and quality enhancement of cottonseed oil and protein.
Identification of stable quantitative trait loci (QTLs) across different environments and mapping populations is a prerequisite for marker-assisted selection (MAS) for cotton yield and fiber quality. ...To construct a genetic linkage map and to identify QTLs for fiber quality and yield traits, a backcross inbred line (BIL) population of 146 lines was developed from a cross between Upland cotton (
Gossypium hirsutum
) and Egyptian cotton (
Gossypium barbadense
) through two generations of backcrossing using Upland cotton as the recurrent parent followed by four generations of self pollination. The BIL population together with its two parents was tested in five environments representing three major cotton production regions in China. The genetic map spanned a total genetic distance of 2,895 cM and contained 392 polymorphic SSR loci with an average genetic distance of 7.4 cM per marker. A total of 67 QTLs including 28 for fiber quality and 39 for yield and its components were detected on 23 chromosomes, each of which explained 6.65–25.27 % of the phenotypic variation. Twenty-nine QTLs were located on the At subgenome originated from a cultivated diploid cotton, while 38 were on the Dt subgenome from an ancestor that does not produce spinnable fibers. Of the eight common QTLs (12 %) detected in more than two environments, two were for fiber quality traits including one for fiber strength and one for uniformity, and six for yield and its components including three for lint yield, one for seedcotton yield, one for lint percentage and one for boll weight. QTL clusters for the same traits or different traits were also identified. This research represents one of the first reports using a permanent advanced backcross inbred population of an interspecific hybrid population to identify QTLs for fiber quality and yield traits in cotton across diverse environments. It provides useful information for transferring desirable genes from
G. barbadense
to
G. hirsutum
using MAS.
Upland cotton (Gossypium hirstum L.), which produces more than 95% of the world natural cotton fibers, has a narrow genetic base which hinders progress in cotton breeding. Introducing germplasm from ...exotic sources especially from another cultivated tetraploid G. barbadense L. can broaden the genetic base of Upland cotton. However, the breeding potential of introgression lines (ILs) in Upland cotton with G. barbadense germplasm integration has not been well addressed. This study involved six ILs developed from an interspecific crossing and backcrossing between Upland cotton and G. barbadense and represented one of the first studies to investigate breeding potentials of a set of ILs using a full diallel analysis. High mid-parent heterosis was detected in several hybrids between ILs and a commercial cultivar, which also out-yielded the high-yielding cultivar parent in F1, F2 and F3 generations. A further analysis indicated that general ability (GCA) variance was predominant for all the traits, while specific combining ability (SCA) variance was either non-existent or much lower than GCA. The estimated GCA effects and predicted additive effects for parents in each trait were positively correlated (at P<0.01). Furthermore, GCA and additive effects for each trait were also positively correlated among generations (at P<0.05), suggesting that F2 and F3 generations can be used as a proxy to F1 in analyzing combining abilities and estimating genetic parameters. In addition, differences between reciprocal crosses in F1 and F2 were not significant for yield, yield components and fiber quality traits. But maternal effects appeared to be present for seed oil and protein contents in F3. This study identified introgression lines as good general combiners for yield and fiber quality improvement and hybrids with high heterotic vigor in yield, and therefore provided useful information for further utilization of introgression lines in cotton breeding.
Cotton (Gossypium spp.) fibers are single-celled elongated trichomes, the molecular aspects of genetic variation in fiber length (FL) among genotypes are currently unknown. In this study, two ...backcross inbred lines (BILs), i.e., NMGA-062 ("Long") and NMGA-105 ("Short") with 32.1 vs. 27.2 mm in FL, respectively, were chosen to perform RNA-Seq on developing fibers at 10 days post anthesis (DPA). The two BILs differed in 4 quantitative trait loci (QTL) for FL and were developed from backcrosses between G. hirsutum as the recurrent parent and G. barbadense.
In total, 51.7 and 54.3 million reads were obtained and assembled to 49,508 and 49,448 transcripts in the two genotypes, respectively. Of 1551 differentially expressed genes (DEGs) between the two BILs, 678 were up-regulated and 873 down-regulated in "Long"; and 703 SNPs were identified in 339 DEGs. Further physical mapping showed that 8 DEGs were co-localized with the 4 FL QTL identified in the BIL population containing the two BILs. Four SNP markers in 3 DEGs that showed significant correlations with FL were developed. Among the three candidate genes encoding for proline-rich protein, D-cysteine desulfhydrase, and thaumatin-like protein, a SNP of thaumatin-like protein gene showed consistent correlations with FL across all testing environments.
This study represents one of the first investigations of positional candidate gene approach of QTL in cotton in integrating transcriptome and SNP identification based on RNA-Seq with linkage and physical mapping of QTL and genes, which will facilitate eventual cloning and identification of genes responsible for FL QTL. The candidate genes may serve as the foundation for further in-depth studies of the molecular mechanism of natural variation in fiber elongation.
Calcineurin B-like protein-interacting protein kinases (
), as key regulators, play an important role in plant growth and development and the response to various stresses. In the present study, we ...identified 80 and 78
genes in the
and
, respectively. The phylogenetic and gene structure analysis divided the cotton
genes into five groups which were classified into an exon-rich clade and an exon-poor clade. A synteny analysis showed that segmental duplication contributed to the expansion of
gene family, and purifying selection played a major role in the evolution of the gene family in cotton. Analyses of expression profiles showed that
genes had temporal and spatial specificity and could be induced by various abiotic stresses. Fourteen
genes were found to contain 17 non-synonymous single nucleotide polymorphisms (SNPs) and co-localized with oil or protein content quantitative trait loci (QTLs). Additionally, five SNPs from four
were found to be significantly associated with oil content in one of the three field tests. Although most
genes were not associated with natural variations in cotton oil content, the overexpression of the
gene reduced the oil content and increased C18:1 and C18:1+C18:1d6 in transgenic cotton as compared to wild-type plants. In addition, we predicted the potential molecular regulatory mechanisms of the
genes. In brief, these results enhance our understanding of the roles of
genes in oil synthesis and stress responses.
Small auxin-up RNA (SAUR) gene family is the largest family of early auxin response genes in higher plants, which have been implicated in the regulation of multiple biological processes. However, no ...comprehensive analysis of SAUR genes has been reported in cotton (Gossypium spp.).
In the study, we identified 145, 97, 214, and 176 SAUR homologous genes in the sequenced genomes of G. raimondii, G. arboreum, G. hirsutum, and G. barbadense, respectively. A phylogenetic analysis revealed that the SAUR genes can be classified into 10 groups. A further analysis of chromosomal locations and gene duplications showed that tandem duplication and segmental duplication events contributed to the expansion of the SAUR gene family in cotton. An exon-intron organization and motif analysis revealed the conservation of SAUR-specific domains, and the auxin responsive elements existed in most of the upstream sequences. The expression levels of 16 GhSAUR genes in response to an exogenous application of IAA were determined by a quantitative RT-PCR analysis. The genome-wide RNA-seq data and qRT-PCR analysis of selected SAUR genes in developing fibers revealed their differential expressions. The physical mapping showed that 20 SAUR genes were co-localized with fiber length quantitative trait locus (QTL) hotspots. Single nucleotide polymorphisms (SNPs) were detected for 12 of these 20 genes between G. hirsutum and G. barbadense, but no SNPs were identified between two backcross inbred lines with differing fiber lengths derived from a cross between the two cultivated tetraploids.
This study provides an important piece of genomic information for the SAUR genes in cotton and lays a solid foundation for elucidating the functions of SAUR genes in auxin signaling pathways to regulate cotton growth.
As the longest known single-celled trichomes, cotton (Gossypium L.) fibers constitute a classic model system to investigate cell initiation and elongation. In this study, we used a high-throughput ...transcriptome sequencing technology to identify fiber-initiation-related single nucleotide polymorphism (SNP) markers and differentially expressed genes (DEGs) between the wild-type (WT) Upland cotton (G. hirsutum) Xuzhou 142 and its natural fuzzless-lintless mutant Xuzhou 142 fl. Approximately 700 million high-quality cDNA reads representing over 58 Gb of sequences were obtained, resulting in the identification of 28,610 SNPs--of which 17,479 were novel--from 13,960 expressed genes. Of these SNPs, 50% of SNPs in fl were identical to those of G. barbadense, which suggests the likely origin of the fl mutant from an interspecific hybridization between Xuzhou 142 and an unknown G. barbadense genotype. Of all detected SNPs, 15,555, 12,750, and 305 were classified as non-synonymous, synonymous, and pre-terminated ones, respectively. Moreover, 1,352 insertion/deletion polymorphisms (InDels) were also detected. A total of 865 DEGs were identified between the WT and fl in ovules at -3 and 0 days post-anthesis, with 302 candidate SNPs selected from these DEGs for validation by a high-resolution melting analysis and Sanger sequencing in seven cotton genotypes. The number of genotypic pairwise polymorphisms varied from 43 to 302, indicating that the identified SNPs are reliable. These SNPs should serve as good resources for breeding and genetic studies in cotton.
In upland cotton (Gossypium hirsutum L.), genotypes with the same mature fiber length (FL) might possess different genes and exhibit differential expression of genes related to fiber elongation at ...different fiber developmental stages. However, there is a lack of information on the genetic variation influencing fiber length and its quantitative trait loci (QTLs) during the fiber elongation stage. In this study, a subset of upland cotton accessions was selected based on a previous GWAS conducted in China and grown in multiple environments to determine the dynamic fiber length at 10, 15, 20, and 25 days post-anthesis (DPA) and maturity. The germplasm lines were genotyped with the Cotton 63 K Illumina single-nucleotide polymorphism (SNP) array for GWAS.
A total of 25, 38, 57, 89 and 88 SNPs showed significant correlations with fiber length at 10, 15, 20 and 25 DPA and maturity, respectively. In addition, 60 more promising SNPs were detected in at least two tests and two FL developmental time points, and 20 SNPs were located within the confidence intervals of QTLs identified in previous studies. The fastest fiber-length growth rates were obtained at 10 to 15 DPA in 69 upland cotton lines and at 15 to 20 DPA in 14 upland cotton accessions, and 10 SNPs showed significant correlations with the fiber-length growth rate. A combined transcriptome and qRT-PCR analysis revealed that two genes (D10G1008 and D13G2037) showed differential expression between two long-fiber genotypes and two short-fiber genotypes.
This study provides important new insights into the genetic basis of the time-dependent fiber-length trait and reveals candidate SNPs and genes for improving fiber length in upland cotton.
Cotton (Gossypium spp.) is the most important natural fiber crop worldwide, and cottonseed oil is its most important byproduct. Phospholipid: diacylglycerol acyltransferase (PDAT) is important in TAG ...biosynthesis, as it catalyzes the transfer of a fatty acyl moiety from the sn-2 position of a phospholipid to the sn-3 position of sn-1, 2-diacylglyerol to form triacylglycerol (TAG) and a lysophospholipid. However, little is known about the genes encoding PDATs involved in cottonseed oil biosynthesis.
A comprehensive genome-wide analysis of G. hirsutum, G. barbadense, G. arboreum, and G. raimondii herein identified 12, 11, 6 and 6 PDATs, respectively. These genes were divided into 3 subfamilies, and a PDAT-like subfamily was identified in comparison with dicotyledonous Arabidopsis. All GhPDATs contained one or two LCAT domains at the C-terminus, while most GhPDATs contained a preserved single transmembrane region at the N-terminus. A chromosomal distribution analysis showed that the 12 GhPDAT genes in G. hirsutum were distributed in 10 chromosomes. However, none of the GhPDATs was co-localized with quantitative trait loci (QTL) for cottonseed oil content, suggesting that their sequence variations are not genetically associated with the natural variation in cottonseed oil content. Most GhPDATs were expressed during the cottonseed oil accumulation stage. Ectopic expression of GhPDAT1d increased Arabidopsis seed oil content.
Our comprehensive genome-wide analysis of the cotton PDAT gene family provides a foundation for further studies into the use of PDAT genes to increase cottonseed oil content through biotechnology.
Extensive studies on floral transition in model species have revealed a network of regulatory interactions between proteins that transduce and integrate developmental and environmental signals to ...promote or inhibit the transition to flowering. Previous studies indicated FLOWERING PROMOTING FACTOR 1 (FPF1) gene was involved in the promotion of flowering, but the molecular mechanism was still unclear. Here, FPF1 homologous sequences were screened from diploid Gossypium raimondii L. (D-genome, n = 13) and Gossypium arboreum L. genome (A-genome, n = 13) databases. Orthologous genes from the two species were compared, suggesting that distinctions at nucleic acid and amino acid levels were not equivalent because of codon degeneracy. Six FPF1 homologous genes were identified from the cultivated allotetraploid Gossypium hirsutum L. (AD-genome, n = 26). Analysis of relative transcripts of the six genes in different tissues revealed that this gene family displayed strong tissue-specific expression. GhFPF1, encoding a 12.0-kDa protein (Accession No: KC832319) exerted more transcripts in floral apices of short-season cotton, hinting that it could be involved in floral regulation. Significantly activated APETALA 1 and suppressed FLOWERING LOCUS C expression were induced by over-expression of GhFPF1 in the Arabidopsis Columbia-0 ecotype. In addition, transgenic Arabidopsis displayed a constitutive shade-avoiding phenotype that is characterized by long hypocotyls and petioles, reduced chlorophyll content, and early flowering. We propose that GhFPF1 may be involved in flowering time control and shade-avoidance responses.