Evolutionary divergence of duplicated genomes in newly described allotetraploid cottons.

Allotetraploid cotton (Gossypium) species represents a model system for the study of plant polyploidy, molecular evolution, and domestication. Here, chromosome-scale genome sequences were obtained and assembled for two recently described wild species of tetraploid cotton, Gossypium ekmanianum [(AD)6, Ge] and Gossypium stephensii [(AD)7, Gs], and one early form of domesticated Gossypium hirsutum, race punctatum [(AD)1, Ghp]. Based on phylogenomic analysis, we provide a dated whole-genome level perspective for the evolution of the tetraploid Gossypium clade and resolved the evolutionary relationships of Gs, Ge, and domesticated G. hirsutum. We describe genomic structural variation that arose during Gossypium evolution and describe its correlates-including phenotypic differentiation, genetic isolation, and genetic convergence-that contributed to cotton biodiversity and cotton domestication. Presence/absence variation is prominent in causing cotton genomic structural variations. A presence/absence variation-derived gene encoding a phosphopeptide-binding protein is implicated in increasing fiber length during cotton domestication. The relatively unimproved Ghp offers the potential for gene discovery related to adaptation to environmental challenges. Expanded gene families enoyl-CoA δ isomerase 3 and RAP2-7 may have contributed to abiotic stress tolerance, possibly by targeting plant hormone-associated biochemical pathways. Our results generate a genomic context for a better understanding of cotton evolution and for agriculture.

Genome-wide identification and expression-pattern analysis of sulfate transporter (SULTR) gene family in cotton under multiple abiotic stresses and fiber development.

Sulfate transporter (SULTR) proteins are in charge of the transport and absorption on sulfate substances, and have been reported to play vital roles in the biological processes of plant growth and stress response. However, there were few reports of genome-wide identification and expression-pattern analysis of SULTRs in Hibiscus mutabilis. Gossypium genus is a ideal model for studying the allopolyploidy, therefore two diploid species (G. raimondii and G. arboreum) and two tetraploid species (G. hirsutum and G. barbadense) were chosen in this study to perform bioinformatic analyses, identifying 18, 18, 35, and 35 SULTR members, respectively. All the 106 cotton SULTR genes were utilized to construct the phylogenetic tree together with 11 Arabidopsis thaliana, 13 Oryza sativa, and 8 Zea mays ones, which was divided into Group1-Group4. The clustering analyses of gene structures and 10 conserved motifs among the cotton SULTR genes showed the consistent evolutionary relationship with the phylogenetic tree, and the results of gene-duplication identification among the four representative Gossypium species indicated that genome-wide or segment duplication might make main contributions to the expansion of SULTR gene family in cotton. Having conducted the cis-regulatory element analysis in promoter region, we noticed that the existing salicylic acid (SA), jasmonic acid (JA), and abscisic acid (ABA) elements could have influences with expression levels of cotton SULTR genes. The expression patterns of GhSULTR genes were also investigated on the 7 different tissues or organs and the developing ovules and fibers, most of which were highly expressed in root, stem, sepal, receptacel, ovule at 10 DPA, and fiber at 20 and 25 DPA. In addition, more active regulatory were observed in GhSULTR genes responding to multiple abiotic stresses, and 12 highly expressed genes showed the similar expression patterns in the quantitative Real-time PCR experiments under cold, heat, salt, and drought treatments. These findings broaden our insight into the evolutionary relationships and expression patterns of the SULTR gene family in cotton, and provide the valuable information for further screening the vital candidate genes on trait improvement.

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.

Genome-wide chromatin accessibility landscape and dynamics of transcription factor networks during ovule and fiber development in cotton.

BACKGROUND: The development of cotton fiber is regulated by the orchestrated binding of regulatory proteins to cis-regulatory elements associated with developmental genes. The cis-trans regulatory dynamics occurred throughout the course of cotton fiber development are elusive. Here we generated genome-wide high-resolution DNase I hypersensitive sites (DHSs) maps to understand the regulatory mechanisms of cotton ovule and fiber development. RESULTS: We generated DNase I hypersensitive site (DHS) profiles from cotton ovules at 0 and 3 days post anthesis (DPA) and fibers at 8, 12, 15, and 18 DPA. We obtained a total of 1185 million reads and identified a total of 199,351 DHSs through ~ 30% unique mapping reads. It should be noted that more than half of DNase-seq reads mapped multiple genome locations and were not analyzed in order to achieve a high specificity of peak profile and to avoid bias from repetitive genomic regions. Distinct chromatin accessibilities were observed in the ovules (0 and 3 DPA) compared to the fiber elongation stages (8, 12, 15, and 18 DPA). Besides, the chromatin accessibility during ovules was particularly elevated in genomic regions enriched with transposable elements (TEs) and genes in TE-enriched regions were involved in ovule cell division. We analyzed cis-regulatory modules and revealed the influence of hormones on fiber development from the regulatory divergence of transcription factor (TF) motifs. Finally, we constructed a reliable regulatory network of TFs related to ovule and fiber development based on chromatin accessibility and gene co-expression network. From this network, we discovered a novel TF, WRKY46, which may shape fiber development by regulating the lignin content. CONCLUSIONS: Our results not only reveal the contribution of TEs in fiber development, but also predict and validate the TFs related to fiber development, which will benefit the research of cotton fiber 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.

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.

Protoplast Dissociation and Transcriptome Analysis Provides Insights to Salt Stress Response in Cotton.

As one of the pioneer crops widely planted in saline-alkaline areas, Gossypium provides daily necessities, including natural fiber, vegetable proteins, and edible oils. However, cotton fiber yield and quality are highly influenced by salt stress. Therefore, elucidating the molecular mechanisms of cotton in response to salinity stress is importance to breed new cultivars with high tolerance. In this study, we first developed a method for single-cell RNA-seq based on isolating protoplast from cotton root tips; then, we studied the impact of salinity stress on gene expression profiling and their dynamic changes using the developed high-efficiency method for protoplast dissociation suitable for single-cell RNA-seq. A total of 3391 and 2826 differentially expressed genes (DEGs) were identified in salt-treated samples before and after protoplast dissociation, respectively, which were enriched into several molecular components, including response to stimulus, response to stress, and cellular macromolecule metabolic process by gene ontology (GO) analysis. Plant hormone signal transduction, phenylpropanoid biosynthesis, and MAPK signaling pathway were found to be enriched via Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Twenty-two and nine salinity-responsive DEGs participated in plant hormone signaling and MAPK signaling in roots, before and after protoplast dissociation, respectively; six upregulated DEGs were involved in ABA signaling transduction, namely, Ga04G2111, Ga07G0142, Ga09G2061, Ga10G0262, Ga01G0063, and Ga08G1915 which indicates their potential functions on plants adapting to salt stress. Additionally, 384 and 257 transcription factors (TFs) were differentially expressed in salt-stress roots before and after protoplast dissociation, respectively, of which significantly up-regulated TFs mainly belonged to the AP2/ERF-ERF family, which implied their potential roles responding to salt stress. These results not only provide novel insights to reveal the regulatory networks in plant's root response to salt stress, but also lay the solid foundation for further exploration on cellular heterogeneity by single-cell transcriptome sequencing.

Micro-FTIR combined with curve fitting method to study cellulose crystallinity of developing cotton fibers.

This study aimed to use micro-FTIR with transmission mode to investigate cellulose crystallinity of developing cotton fibers. Compared with ATR-FTIR method, we found that micro-FTIR can obtain more information of cellulose inside of the developing cotton fibers, especially in high wavenumber of 2800-3000 cm-1 region. Combined with curve fitting method, a new IR crystallinity index (CI) method named wax crystallinity index (WCI) was introduced to evaluate the cellulose crystallinity in the development of cotton fibers based on the peak and area ratios of 2900 cm-1/2850 cm-1 and 2900 cm-1/2920 cm-1. The obtained WCI values demonstrated an excellent coefficient of determination with X-ray diffraction (XRD) CI method with the value up to 0.99. This study suggested that micro-FTIR was an effective technique to qualitatively analyze the crystallinity in developing cotton fibers combined with curve fitting method.

Disequilibrium evolution of the Fructose-1,6-bisphosphatase gene family leads to their functional biodiversity in Gossypium species.

BACKGROUND: Fructose-1,6-bisphosphatase (FBP) is a key enzyme in the plant sucrose synthesis pathway, in the Calvin cycle, and plays an important role in photosynthesis regulation in green plants. However, no systemic analysis of FBPs has been reported in Gossypium species. RESULTS: A total of 41 FBP genes from four Gossypium species were identified and analyzed. These FBP genes were sorted into two groups and seven subgroups. Results revealed that FBP family genes were under purifying selection pressure that rendered FBP family members as being conserved evolutionarily, and there was no tandem or fragmental DNA duplication in FBP family genes. Collinearity analysis revealed that a FBP gene was located in a translocated DNA fragment and the whole FBP gene family was under disequilibrium evolution that led to a faster evolutionary progress of the members in G. barbadense and in At subgenome than those in other Gossypium species and in the Dt subgenome, respectively, in this study. Through RNA-seq analyses and qRT-PCR verification, different FBP genes had diversified biological functions in cotton fiber development (two genes in 0 DPA and 1DPA ovules and four genes in 20-25 DPA fibers), in plant responses to Verticillium wilt onset (two genes) and to salt stress (eight genes). CONCLUSION: The FBP gene family displayed a disequilibrium evolution pattern in Gossypium species, which led to diversified functions affecting not only fiber development, but also responses to Verticillium wilt and salt stress. All of these findings provide the foundation for further study of the function of FBP genes in cotton fiber development and in environmental adaptability.

Genome-wide quantitative trait loci reveal the genetic basis of cotton fibre quality and yield-related traits in a Gossypium hirsutum recombinant inbred line population.

Cotton is widely cultivated globally because it provides natural fibre for the textile industry and human use. To identify quantitative trait loci (QTLs)/genes associated with fibre quality and yield, a recombinant inbred line (RIL) population was developed in upland cotton. A consensus map covering the whole genome was constructed with three types of markers (8295 markers, 5197.17 centimorgans (cM)). Six fibre yield and quality traits were evaluated in 17 environments, and 983 QTLs were identified, 198 of which were stable and mainly distributed on chromosomes 4, 6, 7, 13, 21 and 25. Thirty-seven QTL clusters were identified, in which 92.8% of paired traits with significant medium or high positive correlations had the same QTL additive effect directions, and all of the paired traits with significant medium or high negative correlations had opposite additive effect directions. In total, 1297 genes were discovered in the QTL clusters, 414 of which were expressed in two RNA-Seq data sets. Many genes were discovered, 23 of which were promising candidates. Six important QTL clusters that included both fibre quality and yield traits were identified with opposite additive effect directions, and those on chromosome 13 (qClu-chr13-2) could increase fibre quality but reduce yield; this result was validated in a natural population using three markers. These data could provide information about the genetic basis of cotton fibre quality and yield and help cotton breeders to improve fibre quality and yield simultaneously.

Examining two sets of introgression lines across multiple environments reveals background-independent and stably expressed quantitative trait loci of fiber quality in cotton.

KEY MESSAGE: Background-independent (BI) and stably expressed (SE) quantitative trait loci (QTLs) were identified using two sets of introgression lines across multiple environments. Genetic background more greatly affected fiber quality traits than environmental factors. Sixty-one SE-QTLs, including two BI-QTLs, were novel and 48 SE-QTLs, including seven BI-QTLs, were previously reported. Cotton fiber quality traits are controlled by QTLs and are susceptible to environmental influence. Fiber quality improvement is an essential goal in cotton breeding but is hindered by limited knowledge of the genetic basis of fiber quality traits. In this study, two sets of introgression lines of Gossypium hirsutum × G. barbadense were used to dissect the QTL stability of three fiber quality traits (fiber length, strength and micronaire) across environments using 551 simple sequence repeat markers selected from our high-density genetic map. A total of 76 and 120 QTLs were detected in the CCRI36 and CCRI45 backgrounds, respectively. Nine BI-QTLs were found, and 78 (41.71%) of the detected QTLs were reported previously. Thirty-nine and 79 QTLs were SE-QTLs in at least two environments in the CCRI36 and CCRI45 backgrounds, respectively. Forty-eight SE-QTLs, including seven BI-QTLs, were confirmed in previous reports, and 61 SE-QTLs, including two BI-QTLs, were considered novel. These results indicate that genetic background more strongly impacts on fiber quality traits than environmental factors. Twenty-three clusters with BI- and/or SE-QTLs were identified, 19 of which harbored favorable alleles from G. barbadense for two or three fiber quality traits. This study is the first report using two sets of introgression lines to identify fiber quality QTLs across environments in cotton, providing insights into the effect of genetic backgrounds and environments on the QTL expression of fiber quality and important information for the genetic basis underlying fiber quality traits toward QTL cloning and molecular breeding.

Transcriptomic and biochemical analysis of upland cotton (Gossypium hirsutum) and a chromosome segment substitution line from G. hirsutum × G. barbadense in response to Verticillium dahliae infection.

BACKGROUND: Verticillium wilt (VW), also known as "cotton cancer," is one of the most destructive diseases in global cotton production that seriously impacts fiber yield and quality. Despite numerous attempts, little significant progress has been made in improving the VW resistance of upland cotton. The development of chromosome segment substitution lines (CSSLs) from Gossypium hirsutum × G. barbadense has emerged as a means of simultaneously developing new cotton varieties with high-yield, superior fiber, and resistance to VW. RESULTS: In this study, VW-resistant investigations were first conducted in an artificial greenhouse, a natural field, and diseased nursery conditions, resulting in the identification of one stably VW-resistant CSSL, MBI8255, and one VW-susceptible G. hirsutum, CCRI36, which were subsequently subjected to biochemical tests and transcriptome sequencing during V991 infection (0, 1, and 2 days after inoculation). Eighteen root samples with three replications were collected to perform multiple comparisons of enzyme activity and biochemical substance contents. The findings indicated that VW resistance was positively correlated with peroxidase and polyphenol oxidase activity, but negatively correlated with malondialdehyde content. Additionally, RNA sequencing was used for the same root samples, resulting in a total of 77,412 genes, of which 23,180 differentially expressed genes were identified from multiple comparisons between samples. After Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis on the expression profiles identified using Short Time-series Expression Miner, we found that the metabolic process in the biological process, as well as the pathways of phenylpropanoid biosynthesis and plant hormone signal transduction, participated significantly in the response to VW. Gene functional annotation and expression quantity analysis indicated the important roles of the phenylpropanoid metabolic pathway and oxidation-reduction process in response to VW, which also provided plenty of candidate genes related to plant resistance. CONCLUSIONS: This study concentrates on the preliminary response to V991 infection by comparing the VW-resistant CSSL and its VW-susceptible recurrent parent. Not only do our findings facilitate the culturing of new resistant varieties with high yield and superior performance, but they also broaden our understanding of the mechanisms of cotton resistance to VW.

QTL mapping and genetic effect of chromosome segment substitution lines with excellent fiber quality from Gossypium hirsutum × Gossypium barbadense.

Chromosome segment substitution lines (CSSLs) are ideal materials for identifying genetic effects. In this study, CSSL MBI7561 with excellent fiber quality that was selected from BC4F3:5 of CCRI45 (Gossypium hirsutum) × Hai1 (Gossypium barbadense) was used to construct 3 secondary segregating populations with 2 generations (BC5F2 and BC5F2:3). Eighty-one polymorphic markers related to 33 chromosome introgressive segments on 18 chromosomes were finally screened using 2292 SSR markers which covered the whole tetraploid cotton genome. A total of 129 quantitative trait loci (QTL) associated with fiber quality (103) and yield-related traits (26) were detected on 17 chromosomes, explaining 0.85-30.35% of the phenotypic variation; 39 were stable (30.2%), 53 were common (41.1%), 76 were new (58.9%), and 86 had favorable effects on the related traits. More QTL were distributed in the Dt subgenome than in the At subgenome. Twenty-five stable QTL clusters (with stable or common QTL) were detected on 22 chromosome introgressed segments. Finally, the 6 important chromosome introgressed segments (Seg-A02-1, Seg-A06-1, Seg-A07-2, Seg-A07-3, Seg-D07-3, and Seg-D06-2) were identified as candidate chromosome regions for fiber quality, which should be given more attention in future QTL fine mapping, gene cloning, and marker-assisted selection (MAS) breeding.

Dissecting the genetic basis of fiber quality and yield traits in interspecific backcross populations of Gossypium hirsutum × Gossypium barbadense.

Fiber quality and yield are important traits of cotton. Quantitative trait locus (QTL) mapping is a prerequisite for marker-assisted selection (MAS) in cotton breeding. To identify QTLs for fiber quality and yield traits, 4 backcross-generation populations (BC1F1, BC1S1, BC2F1, and BC3F0) were developed from an interspecific cross between CCRI36 (Gossypium hirsutum L.) and Hai1 (G. barbadense L.). A total of 153 QTLs for fiber quality and yield traits were identified based on data from the BC1F1, BC1S1, BC2F1 and BC3F0 populations in the field and from the BC2F1 population in an artificial disease nursery using a high-density genetic linkage map with 2292 marker loci covering 5115.16 centimorgans (cM) from the BC1F1 population. These QTLs were located on 24 chromosomes, and each could explain 4.98-19.80% of the observed phenotypic variations. Among the 153 QTLs, 30 were consistent with those identified previously. Specifically, 23 QTLs were stably detected in 2 or 3 environments or generations, 6 of which were consistent with those identified previously and the other 17 of which were stable and novel. Ten QTL clusters for different traits were found and 9 of them were novel, which explained the significant correlations among some phenotypic traits in the populations. The results including these stable or consensus QTLs provide valuable information for marker-assisted selection (MAS) in cotton breeding and will help better understand the genetic basis of fiber quality and yield traits, which can then be used in QTL cloning.

Comparative transcriptome analysis of cotton fiber development of Upland cotton (Gossypium hirsutum) and Chromosome Segment Substitution Lines from G. hirsutum × G. barbadense.

BACKGROUND: How to develop new cotton varieties possessing high yield traits of Upland cotton and superior fiber quality traits of Sea Island cotton remains a key task for cotton breeders and researchers. While multiple attempts bring in little significant progresses, the development of Chromosome Segment Substitution Lines (CSSLs) from Gossypium barbadense in G. hirsutum background provided ideal materials for aforementioned breeding purposes in upland cotton improvement. Based on the excellent fiber performance and relatively clear chromosome substitution segments information identified by Simple Sequence Repeat (SSR) markers, two CSSLs, MBI9915 and MBI9749, together with the recurrent parent CCRI36 were chosen to conduct transcriptome sequencing during the development stages of fiber elongation and Secondary Cell Wall (SCW) synthesis (from 10DPA and 28DPA), aiming at revealing the mechanism of fiber development and the potential contribution of chromosome substitution segments from Sea Island cotton to fiber development of Upland cotton. RESULTS: In total, 15 RNA-seq libraries were constructed and sequenced separately, generating 705.433 million clean reads with mean GC content of 45.13% and average Q30 of 90.26%. Through multiple comparisons between libraries, 1801 differentially expressed genes (DEGs) were identified, of which the 902 up-regulated DEGs were mainly involved in cell wall organization and response to oxidative stress and auxin, while the 898 down-regulated ones participated in translation, regulation of transcription, DNA-templated and cytoplasmic translation based on GO annotation and KEGG enrichment analysis. Subsequently, STEM software was performed to explicate the temporal expression pattern of DEGs. Two peroxidases and four flavonoid pathway-related genes were identified in the "oxidation-reduction process", which could play a role in fiber development and quality formation. Finally, the reliability of RNA-seq data was validated by quantitative real-time PCR of randomly selected 20 genes. CONCLUSIONS: The present report focuses on the similarities and differences of transcriptome profiles between the two CSSLs and the recurrent parent CCRI36 and provides novel insights into the molecular mechanism of fiber development, and into further exploration of the feasible contribution of G. barbadense substitution segments to fiber quality formation, which will lay solid foundation for simultaneously improving fiber yield and quality of upland cotton through CSSLs.

Quantitative trait loci analysis of Verticillium wilt resistance in interspecific backcross populations of Gossypium hirsutum × Gossypium barbadense.

BACKGROUND: Verticillium wilt (VW) caused by Verticillium dahliae (Kleb) is one of the most destructive diseases of cotton. The identification of highly resistant QTLs or genes in the whole cotton genome is quite important for developing a VW-resistant variety and for further molecular design breeding. RESULTS: In the present study, BC1F1, BC1S1, and BC2F1 populations derived from an interspecific backcross between the highly resistant line Hai1 (Gossypium barbadense L.) and the susceptible variety CCRI36 (G. hirsutum L.) as the recurrent parent were constructed. Quantitative trait loci (QTL) related to VW resistance were detected in the whole cotton genome using a high-density simple sequence repeat (SSR) genetic linkage map from the BC1F1 population, with 2292 loci covering 5115.16 centiMorgan (cM) of the cotton (AD) genome, and the data concerning VW resistance that were obtained from four dates of BC2F1 in the artificial disease nursery and one date of BC1S1 and BC2F1 in the field. A total of 48 QTLs for VW resistance were identified, and 37 of these QTLs had positive additive effects, which indicated that the G. barbadense alleles increased resistance to VW and decreased the disease index (DI) by about 2.2-10.7. These QTLs were located on 19 chromosomes, in which 33 in the A subgenome and 15 QTLs in the D subgenome. The 6 QTLs were found to be stable. The 6 QTLs were consistent with those identified previously, and another 42 were new, unreported QTLs, of which 31 QTLs were from G. barbadense. By meta-analysis, 17 QTL hotspot regions were identified and 10 of them were new, unreported hotspot regions. 29 QTLs in this paper were in 12 hotspot regions and were all from G. barbadense. CONCLUSIONS: These stable or consensus QTL regions warrant further investigation to better understand the genetics and molecular mechanisms underlying VW resistance. This study provides useful information for further comparative analysis and marker-assisted selection in the breeding of disease-resistant cotton. It may also lay an important foundation for gene cloning and further molecular design breeding for the entire cotton genome.

Genome-wide identification, phylogeny, and expression analysis of pectin methylesterases reveal their major role in cotton fiber development.

BACKGROUND: Pectin methylesterase (PME, EC 3.1.1.11) is a hydrolytic enzyme that utilizes pectin as substrates, and plays a significant role in regulating pectin reconstruction thereby regulating plant growth. Pectin is one of the important components of the plant cell wall, which forms the main structural material of cotton fiber. In this research, cotton genome information was used to identify PMEs. RESULTS: We identified 80 (GaPME01-GaPME80) PME genes from diploid G. arboreum (A genome), 78 (GrPME01-GrPME78) PME genes from G. raimondii (D genome), and 135 (GhPME001-GhPME135) PME genes from tetraploid cotton G. hirsutum (AD genome). We further analyzed their gene structure, conserved domain, gene expression, and systematic evolution to lay the foundation for deeper research on the function of PMEs. Phylogenetic data indicated that members from the same species demonstrated relatively high sequence identities and genetic similarities. Analysis of gene structures showed that most of the PMEs genes had 2-3 exons, with a few having a variable number of exons from 4 to 6. There are nearly no differences in the gene structure of PMEs among the three (two diploid and one tetraploid) cotton species. Selective pressure analysis showed that the Ka/Ks value for each of the three cotton species PME families was less than one. CONCLUSION: Conserved domain analysis showed that PMEs members had a relatively conserved C-terminal pectinesterase domain (PME) while the N-terminus was less conserved. Moreover, some of the family members contained a pectin methylesterase inhibitor (PMEI) domain. The Ka/Ks ratios suggested that the duplicated PMEs underwent purifying selection after the duplication events. This study provided an important basis for further research on the functions of cotton PMEs. Results from qRT-PCR indicated that the expression level of different PMEs at various fiber developmental stages was different. Moreover, some of the PMEs showed fiber predominant expression in secondary wall thickening indicating tissue-specific expression patterns.

Identification of stable quantitative trait loci (QTLs) for fiber quality traits across multiple environments in Gossypium hirsutum recombinant inbred line population.

BACKGROUND: The identification of quantitative trait loci (QTLs) that are stable and consistent across multiple environments and populations plays an essential role in marker-assisted selection (MAS). In the present study, we used 28,861 simple sequence repeat (SSR) markers, which included 12,560 Gossypium raimondii (D genome) sequence-based SSR markers to identify polymorphism between two upland cotton strains 0-153 and sGK9708. A total of 851 polymorphic primers were finally selected and used to genotype 196 recombinant inbred lines (RIL) derived from a cross between 0 and 153 and sGK9708 and used to construct a linkage map. The RIL population was evaluated for fiber quality traits in six locations in China for five years. Stable QTLs identified in this intraspecific cross could be used in future cotton breeding program and with fewer obstacles. RESULTS: The map covered a distance of 4,110 cM, which represents about 93.2 % of the upland cotton genome, and with an average distance of 5.2 cM between adjacent markers. We identified 165 QTLs for fiber quality traits, of which 47 QTLs were determined to be stable across multiple environments. Most of these QTLs aggregated into clusters with two or more traits. A total of 30 QTL clusters were identified which consisted of 103 QTLs. Sixteen clusters in the At sub-genome comprised 44 QTLs, whereas 14 clusters in the Dt sub-genome that included 59 QTLs for fiber quality were identified. Four chromosomes, including chromosome 4 (c4), c7, c14, and c25 were rich in clusters harboring 5, 4, 5, and 6 clusters respectively. A meta-analysis was performed using Biomercator V4.2 to integrate QTLs from 11 environmental datasets on the RIL populations of the above mentioned parents and previous QTL reports. Among the 165 identified QTLs, 90 were identified as common QTLs, whereas the remaining 75 QTLs were determined to be novel QTLs. The broad sense heritability estimates of fiber quality traits were high for fiber length (0.93), fiber strength (0.92), fiber micronaire (0.85), and fiber uniformity (0.80), but low for fiber elongation (0.27). Meta-clusters on c4, c7, c14 and c25 were identified as stable QTL clusters and were considered more valuable in MAS for the improvement of fiber quality of upland cotton. CONCLUSION: Multiple environmental evaluations of an intraspecific RIL population were conducted to identify stable QTLs. Meta-QTL analyses identified a common chromosomal region that plays an important role in fiber development. Therefore, QTLs identified in the present study are an ideal candidate for MAS in cotton breeding programs to improve fiber quality.

Construction of a high-density genetic map by specific locus amplified fragment sequencing (SLAF-seq) and its application to Quantitative Trait Loci (QTL) analysis for boll weight in upland cotton (Gossypium hirsutum.).

BACKGROUND: Upland Cotton (Gossypium hirsutum) is one of the most important worldwide crops it provides natural high-quality fiber for the industrial production and everyday use. Next-generation sequencing is a powerful method to identify single nucleotide polymorphism markers on a large scale for the construction of a high-density genetic map for quantitative trait loci mapping. RESULTS: In this research, a recombinant inbred lines population developed from two upland cotton cultivars 0-153 and sGK9708 was used to construct a high-density genetic map through the specific locus amplified fragment sequencing method. The high-density genetic map harbored 5521 single nucleotide polymorphism markers which covered a total distance of 3259.37 cM with an average marker interval of 0.78 cM without gaps larger than 10 cM. In total 18 quantitative trait loci of boll weight were identified as stable quantitative trait loci and were detected in at least three out of 11 environments and explained 4.15-16.70 % of the observed phenotypic variation. In total, 344 candidate genes were identified within the confidence intervals of these stable quantitative trait loci based on the cotton genome sequence. These genes were categorized based on their function through gene ontology analysis, Kyoto Encyclopedia of Genes and Genomes analysis and eukaryotic orthologous groups analysis. CONCLUSIONS: This research reported the first high-density genetic map for Upland Cotton (Gossypium hirsutum) with a recombinant inbred line population using single nucleotide polymorphism markers developed by specific locus amplified fragment sequencing. We also identified quantitative trait loci of boll weight across 11 environments and identified candidate genes within the quantitative trait loci confidence intervals. The results of this research would provide useful information for the next-step work including fine mapping, gene functional analysis, pyramiding breeding of functional genes as well as marker-assisted selection.

Exploring the plasmodesmata callose-binding protein gene family in upland cotton: unraveling insights for enhancing fiber length.

Plasmodesmata are transmembrane channels embedded within the cell wall that can facilitate the intercellular communication in plants. Plasmodesmata callose-binding (PDCB) protein that associates with the plasmodesmata contributes to cell wall extension. Given that the elongation of cotton fiber cells correlates with the dynamics of the cell wall, this protein can be related to the cotton fiber elongation. This study sought to identify PDCB family members within the Gossypium. hirsutum genome and to elucidate their expression profiles. A total of 45 distinct family members were observed through the identification and screening processes. The analysis of their physicochemical properties revealed the similarity in the amino acid composition and molecular weight across most members. The phylogenetic analysis facilitated the construction of an evolutionary tree, categorizing these members into five groups mainly distributed on 20 chromosomes. The fine mapping results facilitated a tissue-specific examination of group V, revealing that the expression level of GhPDCB9 peaked five days after flowering. The VIGS experiments resulted in a marked decrease in the gene expression level and a significant reduction in the mature fiber length, averaging a shortening of 1.43-4.77 mm. The results indicated that GhPDCB9 played a pivotal role in the cotton fiber development and served as a candidate for enhancing cotton yield.