Quantitative Trait Locus Mapping for Plant Height and Branch Number in CCRI70 Recombinant Inbred Line Population of Upland Cotton (Gossypium hirsutum).

Upland cotton accounts for a high percentage (95%) of the world's cotton production. Plant height (PH) and branch number (BN) are two important agronomic traits that have an impact on improving the level of cotton mechanical harvesting and cotton yield. In this research, a recombinant inbred line (RIL) population with 250 lines developed from the variety CCRI70 was used for constructing a high-density genetic map and identification of quantitative trait locus (QTL). The results showed that the map harbored 8298 single nucleotide polymorphism (SNP) markers, spanning a total distance of 4876.70 centimorgans (cMs). A total of 69 QTLs for PH (9 stable) and 63 for BN (11 stable) were identified and only one for PH was reported in previous studies. The QTLs for PH and BN harbored 495 and 446 genes, respectively. Combining the annotation information, expression patterns and previous studies of these genes, six genes could be considered as potential candidate genes for PH and BN. The results could be helpful for cotton researchers to better understand the genetic mechanism of PH and BN development, as well as provide valuable genetic resources for cotton breeders to manipulate cotton plant architecture to meet future demands.

Genome-wide analyses of member identification, expression pattern, and protein-protein interaction of EPF/EPFL gene family in Gossypium.

BACKGROUND: Epidermal patterning factor / -like (EPF/EPFL) gene family encodes a class of cysteine-rich secretory peptides, which are widelyfound in terrestrial plants.Multiple studies has indicated that EPF/EPFLs might play significant roles in coordinating plant development and growth, especially as the morphogenesis processes of stoma, awn, stamen, and fruit skin. However, few research on EPF/EPFL gene family was reported in Gossypium. RESULTS: We separately identified 20 G. raimondii, 24 G. arboreum, 44 G. hirsutum, and 44 G. barbadense EPF/EPFL genes in the 4 representative cotton species, which were divided into four clades together with 11 Arabidopsis thaliana, 13 Oryza sativa, and 17 Selaginella moellendorffii ones based on their evolutionary relationships. The similar gene structure and common motifs indicated the high conservation among the EPF/EPFL members, while the uneven distribution in chromosomes implied the variability during the long-term evolutionary process. Hundreds of collinearity relationships were identified from the pairwise comparisons of intraspecifc and interspecific genomes, which illustrated gene duplication might contribute to the expansion of cotton EPF/EPFL gene family. A total of 15 kinds of cis-regulatory elements were predicted in the promoter regions, and divided into three major categories relevant to the biological processes of development and growth, plant hormone response, and abiotic stress response. Having performing the expression pattern analyses with the basic of the published RNA-seq data, we found most of GhEPF/EPFL and GbEPF/EPFL genes presented the relatively low expression levels among the 9 tissues or organs, while showed more dramatically different responses to high/low temperature and salt or drought stresses. Combined with transcriptome data of developing ovules and fibers and quantitative Real-time PCR results (qRT-PCR) of 15 highly expressed GhEPF/EPFL genes, it could be deduced that the cotton EPF/EPFL genes were closely related with fiber development. Additionally, the networks of protein-protein interacting among EPF/EPFLs concentrated on the cores of GhEPF1 and GhEPF7, and thosefunctional enrichment analyses indicated that most of EPF/EPFLs participate in the GO (Gene Ontology) terms of stomatal development and plant epidermis development, and the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways of DNA or base excision repair. CONCLUSION: Totally, 132 EPF/EPFL genes were identified for the first time in cotton, whose bioinformatic analyses of cis-regulatory elements and expression patterns combined with qRT-PCR experiments to prove the potential functions in the biological processes of plant growth and responding to abiotic stresses, specifically in the fiber development. These results not only provide comprehensive and valuable information for cotton EPF/EPFL gene family, but also lay solid foundation for screening candidate EPF/EPFL genes in further cotton breeding.

Hub Genes in Stable QTLs Orchestrate the Accumulation of Cottonseed Oil in Upland Cotton via Catalyzing Key Steps of Lipid-Related Pathways.

Upland cotton is the fifth-largest oil crop in the world, with an average supply of nearly 20% of vegetable oil production. Cottonseed oil is also an ideal alternative raw material to be efficiently converted into biodiesel. However, the improvement in kernel oil content (KOC) of cottonseed has not received sufficient attention from researchers for a long time, due to the fact that the main product of cotton planting is fiber. Previous studies have tagged QTLs and identified individual candidate genes that regulate KOC of cottonseed. The regulatory mechanism of oil metabolism and accumulation of cottonseed are still elusive. In the current study, two high-density genetic maps (HDGMs), which were constructed based on a recombinant inbred line (RIL) population consisting of 231 individuals, were used to identify KOC QTLs. A total of forty-three stable QTLs were detected via these two HDGM strategies. Bioinformatic analysis of all the genes harbored in the marker intervals of the stable QTLs revealed that a total of fifty-one genes were involved in the pathways related to lipid biosynthesis. Functional analysis via coexpression network and RNA-seq revealed that the hub genes in the co-expression network that also catalyze the key steps of fatty acid synthesis, lipid metabolism and oil body formation pathways (ACX4, LACS4, KCR1, and SQD1) could jointly orchestrate oil accumulation in cottonseed. This study will strengthen our understanding of oil metabolism and accumulation in cottonseed and contribute to KOC improvement in cottonseed in the future, enhancing the security and stability of worldwide food supply.

QTL Verification and Candidate Gene Screening of Fiber Quality and Lint Percentage in the Secondary Segregating Population of Gossypium hirsutum.

Fiber quality traits, especially fiber strength, length, and micronaire (FS, FL, and FM), have been recognized as critical fiber attributes in the textile industry, while the lint percentage (LP) was an important indicator to evaluate the cotton lint yield. So far, the genetic mechanism behind the formation of these traits is still unclear. Quantitative trait loci (QTL) identification and candidate gene validation provide an effective methodology to uncover the genetic and molecular basis of FL, FS, FM, and LP. A previous study identified three important QTL/QTL cluster loci, harboring at least one of the above traits on chromosomes A01, A07, and D12 via a recombinant inbred line (RIL) population derived from a cross of Lumianyan28 (L28) × Xinluzao24 (X24). A secondary segregating population (F2) was developed from a cross between L28 and an RIL, RIL40 (L28 × RIL40). Based on the population, genetic linkage maps of the previous QTL cluster intervals on A01 (6.70-10.15 Mb), A07 (85.48-93.43 Mb), and D12 (0.40-1.43 Mb) were constructed, which span 12.25, 15.90, and 5.56 cM, with 2, 14, and 4 simple sequence repeat (SSR) and insertion/deletion (Indel) markers, respectively. QTLs of FL, FS, FM, and LP on these three intervals were verified by composite interval mapping (CIM) using WinQTL Cartographer 2.5 software via phenotyping of F2 and its derived F2:3 populations. The results validated the previous primary QTL identification of FL, FS, FM, and LP. Analysis of the RNA-seq data of the developing fibers of L28 and RIL40 at 10, 20, and 30 days post anthesis (DPA) identified seven differentially expressed genes (DEGs) as potential candidate genes. qRT-PCR verified that five of them were consistent with the RNA-seq result. These genes may be involved in regulating fiber development, leading to the formation of FL, FS, FM, and LP. This study provides an experimental foundation for further exploration of these functional genes to dissect the genetic mechanism of cotton fiber development.

Comparative transcriptional and co-expression network analysis of two upland cotton accessions with extreme phenotypic differences reveals molecular mechanisms of fiber development.

Introduction: Upland cotton (Gossypium hirsutum) is the main source of natural fiber in the global textile industry, and thus its fiber quality and yield are important parameters. In this study, comparative transcriptomics was used to analyze differentially expressed genes (DEGs) due to its ability to effectively screen candidate genes during the developmental stages of cotton fiber. However, research using this method is limited, particularly on fiber development. The aim of this study was to uncover the molecular mechanisms underlying the whole period of fiber development and the differences in transcriptional levels. Methods: Comparative transcriptomes are used to analyze transcriptome data and to screen for differentially expressed genes. STEM and WGCNA were used to screen for key genes involved in fiber development. qRT-PCR was performed to verify gene expression of selected DEGs and hub genes. Results: Two accessions of upland cotton with extreme phenotypic differences, namely EZ60 and ZR014121, were used to carry out RNA sequencing (RNA-seq) on fiber samples from different fiber development stages. The results identified 704, 376, 141, 269, 761, and 586 genes that were upregulated, and 1,052, 476, 355, 259, 702, and 847 genes that were downregulated at 0, 5, 10, 15, 20, and 25 days post anthesis, respectively. Similar expression patterns of DEGs were monitored using short time-series expression miner (STEM) analysis, and associated pathways of DEGs within profiles were investigated. In addition, weighted gene co-expression network analysis (WGCNA) identified five key modules in fiber development and screened 20 hub genes involved in the development of fibers. Discussion: Through the annotation of the genes, it was found that the excessive expression of resistance-related genes in the early fiber development stages affects the fiber yield, whereas the sustained expression of cell elongation-related genes is critical for long fibers. This study provides new information that can be used to improve fibers in newly developed upland cotton genotypes.

Genome-Wide Analysis and Functional Characterization of LACS Gene Family Associated with Lipid Synthesis in Cotton (Gossypium spp.).

Cotton (Gossypium spp.) is the fifth largest oil crop in the world, and cottonseed provides abundant vegetable oil resources and industrial bioenergy fuels for people; therefore, it is of practical significance to increase the oil content of cotton seeds for improving the oil yield and economic benefits of planting cotton. Long-chain acyl-coenzyme A (CoA) synthetase (LACS) capable of catalyzing the formation of acyl-CoAs from free fatty acids has been proven to significantly participate in lipid metabolism, of which whole-genome identification and functional characterization of the gene family have not yet been comprehensively analyzed in cotton. In this study, a total of sixty-five LACS genes were confirmed in two diploid and two tetraploid Gossypium species, which were divided into six subgroups based on phylogenetic relationships with twenty-one other plants. An analysis of protein motif and genomic organizations displayed structural and functional conservation within the same group but diverged among the different group. Gene duplication relationship analysis illustrates the LACS gene family in large scale expansion through WGDs/segmental duplications. The overall Ka/Ks ratio indicated the intense purifying selection of LACS genes in four cotton species during evolution. The LACS genes promoter elements contain numerous light response cis-elements associated with fatty acids synthesis and catabolism. In addition, the expression of almost all GhLACS genes in high seed oil were higher compared to those in low seed oil. We proposed LACS gene models and shed light on their functional roles in lipid metabolism, demonstrating their engineering potential for modulating TAG synthesis in cotton, and the genetic engineering of cottonseed oil provides a theoretical basis.

Genome-wide characterization of trichome birefringence-like genes provides insights into fiber yield improvement.

Cotton is an important fiber crop. The cotton fiber is an extremely long trichome that develops from the epidermis of an ovule. The trichome is a general and multi-function plant organ, and trichome birefringence-like (TBL) genes are related to trichome development. At the genome-wide scale, we identified TBLs in four cotton species, comprising two cultivated tetraploids (Gossypium hirsutum and G. barbadense) and two ancestral diploids (G. arboreum and G. raimondii). Phylogenetic analysis showed that the TBL genes clustered into six groups. We focused on GH_D02G1759 in group IV because it was located in a lint percentage-related quantitative trait locus. In addition, we used transcriptome profiling to characterize the role of TBLs in group IV in fiber development. The overexpression of GH_D02G1759 in Arabidopsis thaliana resulted in more trichomes on the stems, thereby confirming its function in fiber development. Moreover, the potential interaction network was constructed based on the co-expression network, and it was found that GH_D02G1759 may interact with several genes to regulate fiber development. These findings expand our knowledge of TBL family members and provide new insights for cotton molecular breeding.

Fine mapping and candidate gene analysis of qFL-A12-5: a fiber length-related QTL introgressed from Gossypium barbadense into Gossypium hirsutum.

KEY MESSAGE: The fiber length-related qFL-A12-5 identified in CSSLs introgressed from Gossypium barbadense into Gossypium hirsutum was fine-mapped to an 18.8 kb region on chromosome A12, leading to the identification of the GhTPR gene as a potential regulator of cotton fiber length. Fiber length is a key determinant of fiber quality in cotton, and it is a key target of artificial selection for breeding and domestication. Although many fiber length-related quantitative trait loci have been identified, there are few reports on their fine mapping or candidate gene validation, thus hampering efforts to understand the mechanistic basis of cotton fiber development. Our previous study identified the qFL-A12-5 associated with superior fiber quality on chromosome A12 in the chromosome segment substitution line (CSSL) MBI7747 (BC4F3:5). A single segment substitution line (CSSL-106) screened from BC6F2 was backcrossed to construct a larger segregation population with its recurrent parent CCRI45, thus enabling the fine mapping of 2852 BC7F2 individuals using denser simple sequence repeat markers to narrow the qFL-A12-5 to an 18.8 kb region of the genome, in which six annotated genes were identified in Gossypium hirsutum. Quantitative real-time PCR and comparative analyses led to the identification of GH_A12G2192 (GhTPR) encoding a tetratricopeptide repeat-like superfamily protein as a promising candidate gene for qFL-A12-5. A comparative analysis of the protein-coding regions of GhTPR among Hai1, MBI7747, and CCRI45 revealed two non-synonymous mutations. The overexpression of GhTPR resulted in longer roots in Arabidopsis, suggesting that GhTPR may regulate cotton fiber development. These results provide a foundation for future efforts to improve cotton fiber length.

Genome-Wide Identification and Characterization of the PPO Gene Family in Cotton (Gossypium) and Their Expression Variations Responding to Verticillium Wilt Infection.

Polyphenol oxidases (PPOs) are copper-binding metalloproteinases encoded by nuclear genes, ubiquitously existing in the plastids of microorganisms, plants, and animals. As one of the important defense enzymes, PPOs have been reported to participate in the resistant processes that respond to diseases and insect pests in multiple plant species. However, PPO gene identification and characterization in cotton and their expression patterns under Verticillium wilt (VW) treatment have not been clearly studied. In this study, 7, 8, 14, and 16 PPO genes were separately identified from Gossypium arboreum, G. raimondii, G. hirsutum, and G. barbadense, respectively, which were distributed within 23 chromosomes, though mainly gathered in chromosome 6. The phylogenetic tree manifested that all the PPOs from four cotton species and 14 other plants were divided into seven groups, and the analyses of the conserved motifs and nucleotide sequences showed highly similar characteristics of the gene structure and domains in the cotton PPO genes. The dramatically expressed differences were observed among the different organs at various stages of growth and development or under the diverse stresses referred to in the published RNA-seq data. Quantitative real-time PCR (qRT-PCR) experiments were also performed on the GhPPO genes in the roots, stems, and leaves of VW-resistant MBI8255 and VW-susceptible CCRI36 infected with Verticillium dahliae V991, proving the strong correlation between PPO activity and VW resistance. A comprehensive analysis conducted on cotton PPO genes contributes to the screening of the candidate genes for subsequent biological function studies, which is also of great significance for the in-depth understanding of the molecular genetic basis of cotton resistance to VW.

Genome-wide artificial introgressions of Gossypium barbadense into G. hirsutum reveal superior loci for simultaneous improvement of cotton fiber quality and yield traits.

INTRODUCTION: The simultaneous improvement of fiber quality and yield for cotton is strongly limited by the narrow genetic backgrounds of Gossypium hirsutum (Gh) and the negative genetic correlations among traits. An effective way to overcome the bottlenecks is to introgress the favorable alleles of Gossypium barbadense (Gb) for fiber quality into Gh with high yield. OBJECTIVES: This study was to identify superior loci for the improvement of fiber quality and yield. METHODS: Two sets of chromosome segment substitution lines (CSSLs) were generated by crossing Hai1 (Gb, donor-parent) with cultivar CCRI36 (Gh) and CCRI45 (Gh) as genetic backgrounds, and cultivated in 6 and 8 environments, respectively. The kmer genotyping strategy was improved and applied to the population genetic analysis of 743 genomic sequencing data. A progeny segregating population was constructed to validate genetic effects of the candidate loci. RESULTS: A total of 68,912 and 83,352 genome-wide introgressed kmers were identified in the CCRI36 and CCRI45 populations, respectively. Over 90 % introgressions were homologous exchanges and about 21 % were reverse insertions. In total, 291 major introgressed segments were identified with stable genetic effects, of which 66(22.98 %), 64(21.99 %), 35(12.03 %), 31(10.65 %) and 18(6.19 %) were beneficial for the improvement of fiber length (FL), strength (FS), micronaire, lint-percentage (LP) and boll-weight, respectively. Thirty-nine introgression segments were detected with stable favorable additive effects for simultaneous improvement of 2 or more traits in Gh genetic background, including 6 could increase FL/FS and LP. The pyramiding effects of 3 pleiotropic segments (A07:C45Clu-081, D06:C45Clu-218, D02:C45Clu-193) were further validated in the segregating population. CONCLUSION: The combining of genome-wide introgressions and kmer genotyping strategy showed significant advantages in exploring genetic resources. Through the genome-wide comprehensive mining, a total of 11 clusters (segments) were discovered for the stable simultaneous improvement of FL/FS and LP, which should be paid more attention in the future.

AAQSP increases mapping resolution of stable QTLs through applying NGS-BSA in multiple genetic backgrounds.

KEY MESSAGE: In this study, we present AAQSP as an extension of existing NGS-BSA applications for identifying stable QTLs at high resolution. GhPAP16 and GhIQD14 fine mapped on chromosome D09 of upland cotton are identified as important candidate genes for lint percentage (LP). Bulked segregant analysis combined with next generation sequencing (NGS-BSA) allows rapid identification of genome sequence differences responsible for phenotypic variation. The NGS-BSA approach applied to crops mainly depends on comparing two bulked DNA samples of individuals from an F2 population. Since some F2 individuals still maintain high heterozygosity, heterosis will exert complications in pursuing NGS-BSA in such populations. In addition, the genetic background influences the stability of gene expression in crops, so some QTLs mapped in one segregating population may not be widely applied in crop improvement. The AAQSP (Association Analysis of QTL-seq on Semi-homologous Populations) reported in our study combines the optimized scheme of constructing BSA bulks with NGS-BSA analysis in two (or more) different parental genetic backgrounds for isolating the stable QTLs. With application of AAQSP strategy and construction of a high-density linkage map, we have successfully identified a QTL significantly related to lint percentage (LP) in cultivated upland cotton, followed by map-based cloning to dissect two candidate genes, GhPAP16 and GhIQD14. This study demonstrated that AAQSP can efficiently identify stable QTLs for complex traits of interest, and thus accelerate the genetic improvement of upland cotton and other crop plants.

Genome-wide association study reveals novel quantitative trait loci and candidate genes of lint percentage in upland cotton based on the CottonSNP80K array.

KEY MESSAGE: Thirty-four SNPs corresponding with 22 QTLs for lint percentage, including 13 novel QTLs, was detected via GWAS. Two candidate genes underlying this trait were also identified. Cotton (Gossypium spp.) is an important natural textile fiber and oilseed crop cultivated worldwide. Lint percentage (LP, %) is one of the important yield components, and increasing LP is a core goal of cotton breeding improvement. However, the genetic and molecular mechanisms underlying LP in upland cotton remain unclear. Here, we performed a genome-wide association study (GWAS) for LP based on 254 upland cotton accessions in four environments as well as the best linear unbiased predictors using the high-density CottonSNP80K array. In total, 41,413 high-quality single-nucleotide polymorphisms (SNPs) were screened, and 34 SNPs within 22 quantitative trait loci (QTLs) were significantly associated with LP. In total, 175 candidate genes were identified from two major genomic loci (GR1 and GR2), and 50 hub genes were identified through GO enrichment and weighted gene co-expression network analysis. Two candidate genes (Gh_D01G0162 and Gh_D07G0463), which may participate in early fiber development to affect the number of fiber protrusions and LP, were also identified. Their genetic variation and expression were verified by linkage disequilibrium blocks, haplotypes, and quantitative real-time polymerase chain reaction, respectively. The weighted gene interaction network analysis showed that the expression of Gh_D07G0463 was significantly correlated with that of Gh_D01G0162. These identified SNPs, QTLs and candidate genes provide important insights into the genetic and molecular mechanisms underlying variations in LP and serve as a foundation for LP improvement via marker-assisted breeding.

Linkage and association analyses reveal that hub genes in energy-flow and lipid biosynthesis pathways form a cluster in upland cotton.

Upland cotton is an important allotetraploid crop that provides both natural fiber for the textile industry and edible vegetable oil for the food or feed industry. To better understand the genetic mechanism that regulates the biosynthesis of storage oil in cottonseed, we identified the genes harbored in the major quantitative trait loci/nucleotides (QTLs/QTNs) of kernel oil content (KOC) in cottonseed via both multiple linkage analyses and genome-wide association studies (GWAS). In 'CCRI70' RILs, six stable QTLs were simultaneously identified by linkage analysis of CHIP and SLAF-seq strategies. In '0-153' RILs, eight stable QTLs were detected by consensus linkage analysis integrating multiple strategies. In the natural panel, thirteen and eight loci were associated across multiple environments with two algorithms of GWAS. Within the confidence interval of a major common QTL on chromosome 3, six genes were identified as participating in the interaction network highly correlated with cottonseed KOC. Further observations of gene differential expression showed that four of the genes, LtnD, PGK, LPLAT1, and PAH2, formed hub genes and two of them, FER and RAV1, formed the key genes in the interaction network. Sequence variations in the coding regions of LtnD, FER, PGK, LPLAT1, and PAH2 genes may support their regulatory effects on oil accumulation in mature cottonseed. Taken together, clustering of the hub genes in the lipid biosynthesis interaction network provides new insights to understanding the mechanism of fatty acid biosynthesis and TAG assembly and to further genetic improvement projects for the KOC in cottonseeds.

Multi-environment Evaluations Across Ecological Regions Reveal That the Kernel Oil Content of Cottonseed Is Equally Determined by Genotype and Environment.

Cotton is the fifth-largest oil crop in the world. A high kernel oil content (KOC) and high stability are important cottonseed attributes for food security. In this study, the phenotype of KOC and the genotype-by-environment interaction factors were collectively dissected using 250 recombinant inbred lines, their parental cultivars sGK156 and 901-001, and CCRI70 across multi-environments. ANOVA and correlation analysis showed that both genotype and environment contributed significantly to KOC accumulation. Analyses of additive main effect multiplicative interaction and genotype-by-environment interaction biplot models presented the effects of genotype, environment, and genotype by environment on KOC performance and the stability of the experimental materials. Interaction network analysis revealed that meteorological and geographical factors explained 38% of the total KOC variance, with average daily rainfall contributing the largest positive impact and cumulative rainfall having the largest negative impact on KOC accumulation. This study provides insight into KOC accumulation and could direct selection strategies for improved KOC and field management of cottonseed in the future.

QTL mapping and candidate gene prediction for fiber yield and quality traits in a high-generation cotton chromosome substitution line with Gossypium barbadense segments.

Gossypium provides the foremost natural fiber for supporting the rapid development of the textile industry. Quantitative trait locus (QTL) mapping of fiber yield and quality traits is, thus, of great significance for providing a foundation for the genetic improvement of key target traits in cotton production. In this study, a superior chromosome segment substitution line (CSSL), MBI8255, with high yield and premium fiber quality characteristics was cultivated from the BC5F3:5 lineage derived from G. barbadense Hai1 and G. hirsutum CCRI36, and was chosen to construct a segregation population containing 123 F2 individuals with CCRI36. A total of 71 polymorphic SSR (simple sequence repeat) markers were identified based on a previous high-density linkage map, and 17 QTLs distributed on five chromosomes were detected, of which 10 QTLs for cotton yield explained 0.26-15.41% of phenotypic variations, while 7 QTLs for fiber quality explained 0.84-9.38% of phenotypic variations, separately containing four and one stable QTLs detected from over two environments. Among three identified QTL clusters, only the Chr19 QTL cluster harbored two stable and one unstable QTL for three different traits, and hence this significant region, which included 1546 genes, was subjected to functional enrichment and transcriptome expression analyses, ultimately screening eight candidate genes relevant to fiber development. This study not only provides useful information for the further fine-mapping and functional verification of candidate genes, but also offers a solid foundation for revealing the molecular mechanisms of fiber formation.

Identification and analysis of oil candidate genes reveals the molecular basis of cottonseed oil accumulation in Gossypium hirsutum L.

KEY MESSAGE: Based on the integration of QTL-mapping and regulatory network analyses, five high-confidence stable QTL regions, six candidate genes and two microRNAs that potentially affect the cottonseed oil content were discovered. Cottonseed oil is increasingly becoming a promising target for edible oil with its high content of unsaturated fatty acids. In this study, a recombinant inbred line (RIL) cotton population was constructed to detect quantitative trait loci (QTLs) for the cottonseed oil content. A total of 39 QTLs were detected across eight different environments, of which five QTLs were stable. Forty-three candidate genes potentially involved in carbon metabolism, fatty acid synthesis and triacylglycerol biosynthesis processes were further obtained in the stable QTL regions. Transcriptome analysis showed that nineteen of these candidate genes expressed during  the developing cottonseed ovules and may affect the cottonseed oil content. Besides, transcription factor (TF) and microRNA (miRNA) co-regulatory network analyses based on the nineteen candidate genes suggested that six genes, two core miRNAs (ghr-miR2949b and ghr-miR2949c), and one TF GhHSL1 were considered to be closely associated with the cottonseed oil content. Moreover, four vital genes were validated by quantitative real-time PCR (qRT-PCR). These results provide insights into the oil accumulation mechanism in developing cottonseed ovules through the construction of a detailed oil accumulation model.

Identification of Candidate Cotton Genes Associated With Fiber Length Through Quantitative Trait Loci Mapping and RNA-Sequencing Using a Chromosome Segment Substitution Line.

Fiber length is an important determinant of fiber quality, and it is a quantitative multi-genic trait. Identifying genes associated with fiber length is of great importance for efforts to improve fiber quality in the context of cotton breeding. Integrating transcriptomic information and details regarding candidate gene regions can aid in candidate gene identification. In the present study, the CCRI45 line and a chromosome segment substitution line (CSSL) with a significantly higher fiber length (MBI7747) were utilized to establish F2 and F2:3 populations. Using a high-density genetic map published previously, six quantitative trait loci (QTLs) associated with fiber length and two QTLs associated with fiber strength were identified on four chromosomes. Within these QTLs, qFL-A07-1, qFL-A12-2, qFL-A12-5, and qFL-D02-1 were identified in two or three environments and confirmed by a meta-analysis. By integrating transcriptomic data from the two parental lines and through qPCR analyses, four genes associated with these QTLs including Cellulose synthase-like protein D3 (CSLD3, GH_A12G2259 for qFL-A12-2), expansin-A1 (EXPA1, GH_A12G1972 for qFL-A12-5), plasmodesmata callose-binding protein 3 (PDCB3, GH_A12G2014 for qFL-A12-5), and Polygalacturonase (At1g48100, GH_D02G0616 for qFL-D02-1) were identified as promising candidate genes associated with fiber length. Overall, these results offer a robust foundation for further studies regarding the molecular basis for fiber length and for efforts to improve cotton fiber quality.

Co-expression network and comparative transcriptome analysis for fiber initiation and elongation reveal genetic differences in two lines from upland cotton CCRI70 RIL population.

Upland cotton is the most widely planted for natural fiber around the world, and either lint percentage (LP) or fiber length (FL) is the crucial component tremendously affecting cotton yield and fiber quality, respectively. In this study, two lines MBZ70-053 and MBZ70-236 derived from G. hirsutum CCRI70 recombinant inbred line (RIL) population presenting different phenotypes in LP and FL traits were chosen to conduct RNA sequencing on ovule and fiber samples, aiming at exploring the differences of molecular and genetic mechanisms during cotton fiber initiation and elongation stages. As a result, 249/128, 369/206, 4296/1198 and 3547/2129 up-/down- regulated differentially expressed genes (DGEs) in L2 were obtained at -3, 0, 5 and 10 days post-anthesis (DPA), respectively. Seven gene expression profiles were discriminated using Short Time-series Expression Miner (STEM) analysis; seven modules and hub genes were identified using weighted gene co-expression network analysis. The DEGs were mainly enriched into energetic metabolism and accumulating as well as auxin signaling pathway in initiation and elongation stages, respectively. Meanwhile, 29 hub genes were identified as 14-3-3ω , TBL35, GhACS, PME3, GAMMA-TIP, PUM-7, etc., where the DEGs and hub genes revealed the genetic and molecular mechanisms and differences during cotton fiber development.

Genome-wide quantitative trait loci mapping on Verticillium wilt resistance in 300 chromosome segment substitution lines from Gossypium hirsutum × Gossypium barbadense.

Cotton Verticillium wilt (VW) is a devastating disease seriously affecting fiber yield and quality, and the most effective and economical prevention measure at present is selection and extension of Gossypium varieties harboring high resistance to VW. However, multiple attempts to improve the VW resistance of the most widely cultivated upland cottons have made little significant progress. The introduction of chromosome segment substitution lines (CSSLs) provide the practical solutions for merging the superior genes related with high yield and wide adaptation from Gossypium hirsutum and VW resistance and the excellent fiber quality from Gossypium barbadense. In this study, 300 CSSLs were chosen from the developed BC5F3:5 CSSLs constructed from CCRI36 (G. hirsutum) and Hai1 (G. barbadense) to conduct quantitative trait locus (QTL) mapping of VW resistance, and a total of 40 QTL relevant to VW disease index (DI) were identified. Phenotypic data were obtained from a 2-year investigation in two fields with two replications per year. All the QTL were distributed on 21 chromosomes, with phenotypic variation of 1.05%-10.52%, and 21 stable QTL were consistent in at least two environments. Based on a meta-analysis, 34 novel QTL were identified, while 6 loci were consistent with previously identified QTL. Meanwhile, 70 QTL hotspot regions were detected, including 44 novel regions. This study concentrates on QTL identification and screening for hotspot regions related with VW in the 300 CSSLs, and the results lay a solid foundation not only for revealing the genetic and molecular mechanisms of VW resistance but also for further fine mapping, gene cloning and molecular designing in breeding programs for resistant cotton varieties.

Quantitative Trait Loci and Transcriptome Analysis Reveal Genetic Basis of Fiber Quality Traits in CCRI70 RIL Population of Gossypium hirsutum.

Upland cotton (Gossypium hirsutum) is widely planted around the world for its natural fiber, and producing high-quality fiber is essential for the textile industry. CCRI70 is a hybrid cotton plant harboring superior yield and fiber quality, whose recombinant inbred line (RIL) population was developed from two upland cotton varieties (sGK156 and 901-001) and were used here to investigate the source of high-quality related alleles. Based on the material of the whole population, a high-density genetic map was constructed using specific locus-amplified fragment sequencing (SLAF-seq). It contained 24,425 single nucleotide polymorphism (SNP) markers, spanning a distance of 4,850.47 centimorgans (cM) over 26 chromosomes with an average marker interval of 0.20 cM. In evaluating three fiber quality traits in nine environments to detect multiple environments stable quantitative trait loci (QTLs), we found 289 QTLs, of which 36 of them were stable QTLs and 18 were novel. Based on the transcriptome analysis for two parents and two RILs, 24,941 unique differentially expressed genes (DEGs) were identified, 473 of which were promising genes. For the fiber strength (FS) QTLs, 320 DEGs were identified, suggesting that pectin synthesis, phenylpropanoid biosynthesis, and plant hormone signaling pathways could influence FS, and several transcription factors may regulate fiber development, such as GAE6, C4H, OMT1, AFR18, EIN3, bZIP44, and GAI. Notably, the marker D13_56413025 in qFS-chr18-4 provides a potential basis for enhancing fiber quality of upland cotton via marker-assisted breeding and gene cloning of important fiber quality traits.