Identification of Chromosome Segment Substitution Lines of Gossypium barbadense Introgressed in G. hirsutum and Quantitative Trait Locus Mapping for Fiber Quality and Yield Traits

Unravelling the genetic basis and regulation networks related to fibre quality improvement using chromosome segment substitution lines in cotton

The elucidation of genetic architecture and molecular regulatory networks underlying complex traits remains a significant challenge in life science, largely due to the substantial background effects that arise from epistasis and gene-environment interactions. The chromosome segment substitution line (CSSL) is an ideal material for genetic and molecular dissection of complex traits due to its near-isogenic properties; yet a comprehensive analysis, from the basic identification of substitution segments to advanced regulatory network, is still insufficient. Here, we developed two cotton CSSL populations on the Gossypium hirsutum background, representing wide adaptation and high lint yield, with introgression from G. barbadense, representing superior fibre quality. We sequenced 99 CSSLs that demonstrated significant differences from G. hirsutum in fibre, and characterized 836 dynamic fibre transcriptomes in three crucial developmental stages. We developed a workflow for precise resolution of chromosomal substitution segments; the genome sequencing revealed substitutions collectively representing 87.25% of the G. barbadense genome. Together, the genomic and transcriptomic survey identified 18 novel fibre-quality-related quantitative trait loci with high genetic contributions and the comprehensive landscape of fibre development regulation. Furthermore, analysis determined unique cis-expression patterns in CSSLs to be the driving force for fibre quality alteration; building upon this, the co-expression regulatory network revealed biological relationships among the noted pathways and accurately described the molecular interactions of GhHOX3, GhRDL1 and GhEXPA1 during fibre elongation, along with reliable predictions for their interactions with GhTBA8A5. Our study will enhance more strategic employment of CSSL in crop molecular biology and breeding programmes.

A wheat chromosome segment substitution line series supports characterization and use of progenitor genetic variation

Genome-wide introgression and substitution lines have been developed in many plant species, enhancing mapping precision, gene discovery, and the identification and exploitation of variation from wild relatives. Created over multiple generations of crossing and/or backcrossing accompanied by marker-assisted selection, the resulting introgression lines are a fixed genetic resource. In this study we report the development of spring wheat (Triticum aestivum L.) chromosome segment substitution lines (CSSLs) generated to systematically capture genetic variation from tetraploid (T. turgidum ssp. dicoccoides) and diploid (Aegilops tauschii) progenitor species. Generated in a common genetic background over four generations of backcrossing, this is a base resource for the mapping and characterization of wheat progenitor variation. To facilitate further exploitation the final population was genetically characterized using a high-density genotyping array and a range of agronomic and grain traits assessed to demonstrate the potential use of the populations for trait localization in wheat.

Lint percentage and boll weight QTLs in three excellent upland cotton (Gossypium hirsutum): ZR014121, CCRI60, and EZ60

BACKGROUND

Upland cotton (Gossypium hirsutum L.) is the most economically important species in the cotton genus (Gossypium spp.). Enhancing the cotton yield is a major goal in cotton breeding programs. Lint percentage (LP) and boll weight (BW) are the two most important components of cotton lint yield. The identification of stable and effective quantitative trait loci (QTLs) will aid the molecular breeding of cotton cultivars with high yield.

RESULTS

Genotyping by target sequencing (GBTS) and genome-wide association study (GWAS) with 3VmrMLM were used to identify LP and BW related QTLs from two recombinant inbred line (RIL) populations derived from high lint yield and fiber quality lines (ZR014121, CCRI60 and EZ60). The average call rate of a single locus was 94.35%, and the average call rate of an individual was 92.10% in GBTS. A total of 100 QTLs were identified; 22 of them were overlapping with the reported QTLs, and 78 were novel QTLs. Of the 100 QTLs, 51 QTLs were for LP, and they explained 0.29-9.96% of the phenotypic variation; 49 QTLs were for BW, and they explained 0.41-6.31% of the phenotypic variation. One QTL (qBW-E-A10-1, qBW-C-A10-1) was identified in both populations. Six key QTLs were identified in multiple-environments; three were for LP, and three were for BW. A total of 108 candidate genes were identified in the regions of the six key QTLs. Several candidate genes were positively related to the developments of LP and BW, such as genes involved in gene transcription, protein synthesis, calcium signaling, carbon metabolism, and biosynthesis of secondary metabolites. Seven major candidate genes were predicted to form a co-expression network. Six significantly highly expressed candidate genes of the six QTLs after anthesis were the key genes regulating LP and BW and affecting cotton yield formation.

CONCLUSIONS

A total of 100 stable QTLs for LP and BW in upland cotton were identified in this study; these QTLs could be used in cotton molecular breeding programs. Putative candidate genes of the six key QTLs were identified; this result provided clues for future studies on the mechanisms of LP and BW developments.

Inheritance, QTLs, and Candidate Genes of Lint Percentage in Upland Cotton

Cotton (Gossypium spp.) is an important natural fiber plant. Lint percentage (LP) is one of the most important determinants of cotton yield and is a typical quantitative trait with high variation and heritability. Many cotton LP genetic linkages and association maps have been reported. This work summarizes the inheritance, quantitative trait loci (QTLs), and candidate genes of LP to facilitate LP genetic study and molecular breeding. More than 1439 QTLs controlling LP have been reported. Excluding replicate QTLs, 417 unique QTLs have been identified on 26 chromosomes, including 243 QTLs identified at LOD >3. More than 60 are stable, major effective QTLs that can be used in marker-assisted selection (MAS). More than 90 candidate genes for LP have been reported. These genes encode MYB, HOX, NET, and other proteins, and most are preferentially expressed during fiber initiation and elongation. A putative molecular regulatory model of LP was constructed and provides the foundation for the genetic study and molecular breeding of LP.

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 F₂ and F_2: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.

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.

High-Temperature and Drought-Resilience Traits among Interspecific Chromosome Substitution Lines for Genetic Improvement of Upland Cotton

Upland cotton (Gossypium hirsutum L.) growth and development during the pre-and post-flowering stages are susceptible to high temperature and drought. We report the field-based characterization of multiple morpho-physiological and reproductive stress resilience traits in 11 interspecific chromosome substitution (CS) lines isogenic to each other and the inbred G. hirsutum line TM-1. Significant genetic variability was detected (p < 0.001) in multiple traits in CS lines carrying chromosomes and chromosome segments from CS-B (G. barbadense) and CS-T (G. tomentosum). Line CS-T15sh had a positive effect on photosynthesis (13%), stomatal conductance (33%), and transpiration (24%), and a canopy 6.8 °C cooler than TM-1. The average pollen germination was approximately 8% greater among the CS-B than CS-T lines. Based on the stress response index, three CS lines are identified as heat- and drought-tolerant (CS-T07, CS-B15sh, and CS-B18). The three lines demonstrated enhanced photosynthesis (14%), stomatal conductance (29%), transpiration (13%), and pollen germination (23.6%) compared to TM-1 under field conditions, i.e., traits that would expectedly enhance performance in stressful environments. The generated phenotypic data and stress-tolerance indices on novel CS lines, along with phenotypic methods, would help in developing new cultivars with improved resilience to the effects of global warming.

Localization of quantitative trait loci for cucumber fruit shape by a population of chromosome segment substitution lines

Cucumber fruit shape, a significant agronomic trait, is controlled by quantitative trait loci (QTLs). Feasibility of chromosome segment substitution lines (CSSLs) is well demonstrated to map QTLs, especially the minor-effect ones. To detect and identify QTLs with CSSLs can provide new insights into the underlying mechanisms regarding cucumber fruit shape. In the present study, 71 CSSLs were built from a population of backcross progeny (BC4F2) by using RNS7 (a round-fruit cucumber) as the recurrent parent and CNS21 (a long-stick-fruit cucumber) as the donor parent in order to globally detect QTLs for cucumber fruit shape. With the aid of 114 InDel markers covering the whole cucumber genome, 21 QTLs were detected for fruit shape-related traits including ovary length, ovary diameter, ovary shape index, immature fruit length, immature fruit diameter, immature fruit shape index, mature fruit length, mature fruit diameter and mature fruit shape index, and 4 QTLs for other traits including fruit ground and flesh color, and seed size were detected as well. Together our results provide important resources for the subsequent theoretical and applied researches on cucumber fruit shape and other traits.

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

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.

Map-based cloning of qBWT-c12 discovered brassinosteroid-mediated control of organ size in cotton

Assuring fiber yield stability is the primary objective for cotton breeders since the world population is on the rise, and the demand for cotton fiber is increasing every year. Thus, enhancing average cotton boll weight (BWT) could improve seed cotton production, and ultimately to increase cotton fiber yield. This study accomplished the map-based cloning of a novel boll weight regulating locus, qBWT-c12, in cotton. Bulk segregation analysis detected linked markers, aided in the detection of a stable BWT regulating locus, qBWT-c12, on Chr12 in a novel boll size mutant, BS41. Progeny evaluation confined the qBWT-c12 to a 0.89 cM interval between the AD-A12_07 and AD-FM_44 markers in recombinant derived F₃ and F₄ populations. Homology mapping detected a 40 bp insertion-deletion (InDel) site in the AD-FM_44 clone sequence situated +341 downstream of GhBRH1_A12, which showed complete linkage to the BWT phenotype. The suppressed expression of GhBRH1_A12 suggested its putative involvement during early boll development events in BS41. Although brassinosteroid (BR) biosynthesis and signaling pathway genes were up regulated in different tissues, but the organ growth was suppressed leading to dwarf plants, smaller leaves, and de-morphed smaller bolls in BS41. Thus, a disruption in the BR signal cascade is anticipated and could be related to lower GhBRH1_A12 expression in BS41.This study firstly reported the genetic dissection of boll size regulation of G. barbadense in G. hirsutum background using map-based cloning of a BWT regulating locus, qBWT-c12. Moreover, it also emphasized the putative role GhBRH1_A12 in regulating BR homeostasis and its potential to modulate plant growth and boll development in cotton.

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.

Combined Linkage Mapping and BSA to Identify QTL and Candidate Genes for Plant Height and the Number of Nodes on the Main Stem in Soybean

Soybean is one of the most important food and oil crops in the world. Plant height (PH) and the number of nodes on the main stem (NNMS) are quantitative traits closely related to soybean yield. In this study, we used 208 chromosome segment substitution lines (CSSL) populations constructed using "SN14" and "ZYD00006" for quantitative trait locus (QTL) mapping of PH and NNMS. Combined with bulked segregant analysis (BSA) by extreme materials, 8 consistent QTLs were identified. According to the gene annotation of the QTL interval, a total of 335 genes were obtained. Five of which were associated with PH and NNMS, potentially representing candidate genes. RT-qPCR of these 5 candidate genes revealed two genes with differential relative expression levels in the stems of different materials. Haplotype analysis showed that different single nucleotide polymorphisms (SNPs) between the excellent haplotypes in Glyma.04G251900 and Glyma.16G156700 may be the cause of changes in these traits. These results provide the basis for research on candidate genes and marker-assisted selection (MAS) in soybean breeding.

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 BC₄F_3:5 of CCRI45 (Gossypium hirsutum) × Hai1 (Gossypium barbadense) was used to construct 3 secondary segregating populations with 2 generations (BC₅F₂ and BC₅F_2: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 (BC₁F₁, BC₁S₁, BC₂F₁, and BC₃F₀) 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 BC₁F₁, BC₁S₁, BC₂F₁ and BC₃F₀ populations in the field and from the BC₂F₁ population in an artificial disease nursery using a high-density genetic linkage map with 2292 marker loci covering 5115.16 centimorgans (cM) from the BC₁F₁ 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.

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 BC₄F_3:5 of CCRI45 (Gossypium hirsutum) × Hai1 (Gossypium barbadense) was used to construct 3 secondary segregating populations with 2 generations (BC₅F₂ and BC₅F_2: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.

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.

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.

Development and use of chromosome segment substitution lines as a genetic resource for crop improvement

CSSLs are a complete library of introgression lines with chromosomal segments of usually a distant genotype in an adapted background and are valuable genetic resources for basic and applied research on improvement of complex traits. Chromosome segment substitution lines (CSSLs) are genetic stocks representing the complete genome of any genotype in the background of a cultivar as overlapping segments. Ideally, each CSSL has a single chromosome segment from the donor with a maximum recurrent parent genome recovered in the background. CSSL development program requires population-wide backcross breeding and genome-wide marker-assisted selection followed by selfing. Each line in a CSSL library has a specific marker-defined large donor segment. CSSLs are evaluated for any target phenotype to identify lines significantly different from the parental line. These CSSLs are then used to map quantitative trait loci (QTLs) or causal genes. CSSLs are valuable prebreeding tools for broadening the genetic base of existing cultivars and harnessing the genetic diversity from the wild- and distant-related species. These are resources for genetic map construction, mapping QTLs, genes or gene interactions and their functional analysis for crop improvement. In the last two decades, the utility of CSSLs in identification of novel genomic regions and QTL hot spots influencing a wide range of traits has been well demonstrated in food and commercial crops. This review presents an overview of how CSSLs are developed, their status in major crops and their use in genomic studies and gene discovery.

Identification of circularRNAs and their targets in Gossypium under Verticillium wilt stress based on RNA-seq

Circular RNAs (circRNAs), a class of recently discovered non-coding RNAs, play a role in biological and developmental processes. A recent study showed that circRNAs exist in plants and play a role in their environmental stress responses. However, cotton circRNAs and their role in Verticillium wilt response have not been identified up to now. In this study, two CSSLs (chromosome segment substitution lines) of G.barbadense introgressed into G. hirsutum, CSSL-1 and CSSL-4 (a resistant line and a susceptible line to Verticillium wilt, respectively), were inoculated with V. dahliae for RNA-seq library construction and circRNA analysis. A total of 686 novel circRNAs were identified. CSSL-1 and CSSL-4 had similar numbers of circRNAs and shared many circRNAs in common. However, CSSL-4 differentially expressed approximately twice as many circRNAs as CSSL-1, and the differential expression levels of the common circRNAs were generally higher in CSSL-1 than in CSSL-4. Moreover, two C-RRI comparisons, C-RRI-vs-C-RRM and C-RRI-vs-C-RSI, possessed a large proportion (approximately 50%) of the commonly and differentially expressed circRNAs. These results indicate that the differentially expressed circRNAs may play roles in the Verticillium wilt response in cotton. A total of 280 differentially expressed circRNAs were identified. A Gene Ontology analysis showed that most of the ‘stimulus response’ term source genes were NBS family genes, of which most were the source genes from the differentially expressed circRNAs, indicating that NBS genes may play a role in Verticillium wilt resistance and might be regulated by circRNAs in the disease-resistance process in cotton.

Quantitative Trait Locus Mapping and Candidate Gene Analysis for Verticillium Wilt Resistance Using Gossypium barbadense Chromosomal Segment Introgressed Line

Verticillium wilt (VW) is a soil-borne fungal disease that is caused by Verticillium dahliae Kleb and seriously damages cotton production annually in China. To date, many efforts have been made to improve the resistance of upland cotton against VW, but little progress has been achieved because of a lack of resistant upland cotton to VW. G. barbadense is known to carry high resistance to VW; however, it is difficult to transfer the resistance trait from G. barbadense to upland cotton because of linkage drag and distortion in the interspecific hybrid. In this study, a chromosomal segment introgression line (CSIL), SuVR043, containing a single and homozygous chromosome segment of G. barbadense cv. H7124 D04 (Chr 22), was created and used to construct an F₂ population for mapping of VW resistance quantitative trait loci (QTLs) in the greenhouse. Two major resistance QTLs against nondefoliating V. dahliae isolate Bp2, called qVW-Bp2-1 and qVW-Bp2-2, which were flanked by the markers cgr6409-ZHX37 and ZHX57-ZHX70 and explained an average of 16.38 and 22.36% of the observed phenotypic variation, respectively, were detected in three independent replicate experiments. The genetic distances from cgr6409 to ZHX37 and from ZHX57 to ZHX70 were 2.4 and 0.8 cM, respectively. By analyzing the genome sequence of the qVW-Bp2-1 and qVW-Bp2-2 regions, we determined that the accurate physical distances from cgr6409 to ZHX37 and from ZHX57 to ZHX70 in the G. barbadense genome are 254 and 140 kb, and that those spans 36 and 20 putative genes, respectively. The results of the expression analysis showed significant differences in the expression profiles of GbCYP450, GbTMEM214, and GbRLK among G. barbadense cv. H7124, CSIL SuVR043 and G. hirsutum acc. Sumian 8 at different times after inoculation with V. dahliae isolate Bp2. Virus-induced gene silencing (VIGS) analysis showed that silencing of GbCYP450 and GbTMEM214 decreased H7124 and CSIL SuVR043 resistance to VW. These results form a solid foundation for fine mapping and cloning of resistance genes in the substituted segment and will provide valuable assistance in future efforts to breed for VW resistance.

Genetic and phenotypic effects of chromosome segments introgressed from Gossypium barbadense into Gossypium hirsutum

MBI9915 is an introgression cotton line with excellent fiber quality. It was obtained by advanced backcrossing and continuous inbreeding from an interspecific cross between the upland cotton (Gossypium hirsutum) cultivar CCRI36 as the recurrent parent and the sea island cotton (G. barbadense) cultivar Hai1, as the donor parent. To study the genetic effects of the introgressed chromosome segments in G. hirsutum, an F₂ secondary segregating population of 1537 individuals was created by crossing MBI9915 and CCRI36, and an F_2:3 population was created by randomly selecting 347 individuals from the F₂ generation. Quantitative trait locus (QTL) mapping and interaction for fiber length and strength were identified using IciMapping software. The genotype analysis showed that the recovery rate for MBI9915 was 97.9%, with a total 6 heterozygous segments and 13 homozygous segments. A total of 18 QTLs for fiber quality and 6 QTLs for yield related traits were detected using the two segregating generations. These QTLs were distributed across 7 chromosomes and collectively explained 0.81%–9.51% of the observed phenotypic variations. Six QTLs were consistently detected in two generations and 6 QTLs were identified in previous studies. A total of 13 pairs of interaction for fiber length and 13 pairs of interaction for fiber strength were identified in two generations. Among them, 3 pairs of interaction for fiber length and 3 pairs of interaction for fiber strength could be identified in all generations; 4 pairs of interactions affected fiber length and fiber strength simultaneously. The results clearly showed that 5 chromosome segments (Seg-5-1, Seg-7-1, Seg-8-1, Seg-20-2 and Seg-20-3) have important effects on fiber yield and quality. This study provides the useful information for gene cloning and marker-assisted breeding for excellent fiber related quality.

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.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4077-8) contains supplementary material, which is available to authorized users.

Detection of Favorable QTL Alleles and Candidate Genes for Lint Percentage by GWAS in Chinese Upland Cotton

Improving cotton yield is a major breeding goal for Chinese upland cotton. Lint percentage is an important yield component and a critical economic index for cotton cultivars, and raising the lint percentage has a close relationship to improving cotton lint yield. To investigate the genetic architecture of lint percentage, a diversity panel consisting of 355 upland cotton accessions was grown, and the lint percentage was measured in four different environments. Genotyping was performed with specific-locus amplified fragment sequencing (SLAF-seq). Twelve single-nucleotide polymorphisms (SNPs) associated with lint percentage were detected via a genome-wide association study (GWAS), in which five SNP loci distributed on chromosomes A_t3 (A02) and A_t4 (A08) and contained two major-effect QTLs, which were detected in the best linear unbiased predictions (BLUPs) and in more than three environments simultaneously. Furthermore, favorable haplotypes (FHs) of two major-effect QTLs and 47 putative candidate genes in the two linkage disequilibrium (LD) blocks of these associated loci were identified. The expression levels of these putative candidate genes were estimated using RNA-seq data from ten upland cotton tissues. We found that Gh_A02G1268 was very highly expressed during the early fiber development stage, whereas the gene was poorly expressed in the seed. These results implied that Gh_A02G1268 may determine the lint percentage by regulating seed and fiber development. The favorable QTL alleles and candidate genes for lint percentage identified in this study will have high potential for improving lint yield in future Chinese cotton breeding programs.