Genome structure of cotton revealed by a genome-wide SSR genetic map constructed from a BC1 population between gossypium hirsutum and G. barbadense

Genetic analysis and molecular validation of gene conferring petal spot phenotype in interspecific crosses of cotton

Cotton (Gossypium species) has received considerable interest from the geneticists, cytologists and evolutionary biologists since the last more than a century. Here, we explore the genetics of petal spot in the interspecific derivatives involving tetraploid and diploid cottons; and confirm the location of gene governing petal spot phenotype on chromosome A7 by demonstrating co-segregation of SSR marker NAU 2186 with petal spot phenotype. The presence of petal spot was observed to be dominant over its absence. Petal spot inheritance showed significant deviation from the expected Mendelian ratio in all the segregating populations indicating segregation distortion. The distortion was biased towards the hirsutum parent which has important implications from introgression point of view. We also report a strong association between petal spot and petal margin coloration phenotypes. Extant American cotton varieties generally lack petal spot and margin coloration phenotypes. These petal characteristics can serve as morphological markers during germplasm characterization.

Genomic and Cytogenetic Analysis of Synthetic Polyploids between Diploid and Tetraploid Cotton (Gossypium) Species

Cotton (Gossypium spp.) is the most important natural fiber source in the world. The genetic potential of cotton can be successfully and efficiently exploited by identifying and solving the complex fundamental problems of systematics, evolution, and phylogeny, based on interspecific hybridization of cotton. This study describes the results of interspecific hybridization of G. herbaceum L. (A1-genome) and G. mustelinum Miers ex Watt (AD4-genome) species, obtaining fertile hybrids through synthetic polyploidization of otherwise sterile triploid forms with colchicine (C22H25NO6) treatment. The fertile F1C hybrids were produced from five different cross combinations: (1) G. herbaceum subsp. frutescens × G. mustelinum; (2) G. herbaceum subsp. pseudoarboreum × G. mustelinum; (3) G. herbaceum subsp. pseudoarboreum f. harga × G. mustelinum; (4) G. herbaceum subsp. africanum × G. mustelinum; (5) G. herbaceum subsp. euherbaceum (variety A-833) × G. mustelinum. Cytogenetic analysis discovered normal conjugation of bivalent chromosomes in addition to univalent, open, and closed ring-shaped quadrivalent chromosomes at the stage of metaphase I in the F1C and F2C hybrids. The setting of hybrid bolls obtained as a result of these crosses ranged from 13.8-92.2%, the fertility of seeds in hybrid bolls from 9.7-16.3%, and the pollen viability rates from 36.6-63.8%. Two transgressive plants with long fiber of 35.1-37.0 mm and one plant with extra-long fiber of 39.1-41.0 mm were identified in the F2C progeny of G. herbaceum subsp. frutescens × G. mustelinum cross. Phylogenetic analysis with 72 SSR markers that detect genomic changes showed that tetraploid hybrids derived from the G. herbaceum × G. mustelinum were closer to the species G. mustelinum. The G. herbaceum subsp. frutescens was closer to the cultivated form, and its subsp. africanum was closer to the wild form. New knowledge of the interspecific hybridization and synthetic polyploidization was developed for understanding the genetic mechanisms of the evolution of tetraploid cotton during speciation. The synthetic polyploids of cotton obtained in this study would provide beneficial genes for developing new cotton varieties of the G. hirsutum species, with high-quality cotton fiber and strong tolerance to biotic or abiotic stress. In particular, the introduction of these polyploids to conventional and molecular breeding can serve as a bridge of transferring valuable genes related to high-quality fiber and stress tolerance from different cotton species to the new cultivars.

Validation of candidate gene-based EST-SSR markers for sugar yield in sugarcane

Sugarcane (Saccharum spp.) is a widely cultivated crop that fulfils approximately 75% of the sucrose demand worldwide. Owing to its polyploidy and complex genetic nature, it is difficult to identify and map genes related to complex traits, such as sucrose content. However, association mapping is one of the alternatives for identifying genes or markers for marker-assisted selection. In the present study, EST-SSR primers were obtained from in silico studies. The functionality of each primer was tested using Blast2Go software, and 30 EST-SSR primers related to sugar content were selected. These markers were validated using association analysis. A total of 70 F1 diverse genotypes for sugar content were phenotypes with two check lines. All parameters related to sugar content were recorded. The results showed a significant variation between the genotypes for sugar yield traits such as Brix value, purity, and sucrose content, etc. Correlation studies revealed that the Brix%, sucrose content, and sucrose recovery were significantly correlated. An association analysis was performed using mixed linear model to avoid false positive associations. The association analysis revealed that the SEM 407 marker was significantly associated with Brix% and sucrose content. The SEM 407 primers are putatively related to diphosphate-fructose-6-phosphate 1-phosphotransferase which is associated with Brix% and sucrose content. This functional marker can be used for marker-assisted selection for sugar yield traits in sugarcane that could accelerate the sugarcane breeding program.

Comparative transmission genetics of introgressed chromatin in reciprocal advanced backcross populations in Gossypium (cotton) polyploids

Introgression is a potential source of valuable genetic variation and interspecific introgression lines are important resources for plant breeders to access novel alleles. Experimental advanced-generation backcross populations contain individuals with genomic compositions similar to those resulting from natural interspecific hybridization and provide opportunities to study the nature and transmission pattern of donor chromatin in recipient genomes. Here, we analyze transmission of donor chromatin in reciprocal backcrosses between G. hirsutum and G. barbadense. Across the genome, recurrent backcrossing in both backgrounds yielded donor chromatin at slightly higher frequencies than the Mendelian expectation in BC5F1 plants, while the average frequency of donor alleles in BC5F2 segregating families was less than expected. In the two subgenomes of polyploid cotton, the rate of donor chromatin introgression was similar. Although donor chromatin was tolerated over much of the recipient genomes, 21 regions recalcitrant to donor alleles were identified. Only limited correspondence is observed between the recalcitrant regions in the two backgrounds, suggesting the effect of species background on introgression of donor segments. Genetic breakdown was progressive, with floral abscission and seed inviability ongoing during backcrossing cycles. Regions of either high or low introgression tended to be in terminal chromosomal regions that are generally rich in both genes and crossover events, with long stretches around the centromere having limited crossover activity resulting in relatively constant low introgression frequencies. Constraints on fixation and selection of donor alleles highlights the challenges of utilizing introgression breeding in crop improvement.

Genomic confirmation of Gossypium barbadense introgression into G. hirsutum and a subsequent MAGIC population

Introgression of superior fiber traits from Pima cotton (Gossypium barbadense, GB) into high yield Upland cotton (G. hirsutum) has been a breeding objective for many years in a few breeding programs in the world. However, progress has been very slow due to introgression barriers resulting from whole genome hybridization between the two species. To minimize such barriers, chromosome substitution lines (CS-B) from Pima cotton 3-79 in an Upland cotton cultivar TM-1 were developed. A multiparent advanced generation inter-cross (MAGIC) population consisting of 180 recombinant inbred lines (RILs) was subsequently made using the 18 CS-B lines and three Upland cotton cultivars as parents. In this research, we sequenced the whole genomes of the 21 parents and 180 RILs to examine the G. barbadense introgression. Of the 18 CS-B lines, 11 contained the target GB chromosome or chromosome segment, two contained more than two GB chromosomes, and five did not have the expected introgression. Residual introgression in non-target chromosomes was prevalent in all CS-B lines. A clear structure existed in the MAGIC population and the 180 RILs were distributed into three groups, i.e., high, moderate, and low GB introgression. Large blocks of GB chromosome introgression were still present in some RILs after five cycles of random-mating, an indication of recombination suppression or other unknown reasons present in the population. Identity by descent analysis revealed that the MAGIC RILs contained less introgression than expected. This research presents an insight on understanding the complex problems of introgression between cotton species.

Construction of a high-density genetic map and identification of QTLs related to agronomic and physiological traits in an interspecific (Gossypium hirsutum × Gossypium barbadense) F2 population

Background

Advances in genome sequencing technology, particularly restriction-site associated DNA sequence (RAD-seq) and whole-genome resequencing, have greatly aided the construction of cotton interspecific genetic maps based on single nucleotide polymorphism (SNPs), Indels, and other types of markers. High-density genetic maps can improve accuracy of quantitative trait locus (QTL) mapping, narrow down location intervals, and facilitate identification of the candidate genes.

Result

In this study, 249 individuals from an interspecific F2 population (TM-1 and Hai7124) were re-sequenced, yielding 6303 high-confidence bin markers spanning 5057.13 cM across 26 cotton chromosomes. A total of 3380 recombination hot regions RHRs were identified which unevenly distributed on the 26 chromosomes. Based on this map, 112 QTLs relating to agronomic and physiological traits from seedling to boll opening stage were identified, including 15 loci associated with 14 traits that contained genes harboring nonsynonymous SNPs. We analyzed the sequence and expression of these ten candidate genes and discovered that GhRHD3 (GH_D10G0500) may affect fiber yield while GhGPAT6 (GH_D04G1426) may affect photosynthesis efficiency.

Conclusion

Our research illustrates the efficiency of constructing a genetic map using binmap and QTL mapping on the basis of a certain size of the early-generation population. High-density genetic map features high recombination exchanges in number and distribution. The QTLs and the candidate genes identified based on this high-density genetic map may provide important gene resources for the genetic improvement of cotton.

A method for low-coverage single-gamete sequence analysis demonstrates adherence to Mendel's first law across a large sample of human sperm

Recently published single-cell sequencing data from individual human sperm (n=41,189; 969-3377 cells from each of 25 donors) offer an opportunity to investigate questions of inheritance with improved statistical power, but require new methods tailored to these extremely low-coverage data (∼0.01× per cell). To this end, we developed a method, named rhapsodi, that leverages sparse gamete genotype data to phase the diploid genomes of the donor individuals, impute missing gamete genotypes, and discover meiotic recombination breakpoints, benchmarking its performance across a wide range of study designs. We then applied rhapsodi to the sperm sequencing data to investigate adherence to Mendel's Law of Segregation, which states that the offspring of a diploid, heterozygous parent will inherit either allele with equal probability. While the vast majority of loci adhere to this rule, research in model and non-model organisms has uncovered numerous exceptions whereby 'selfish' alleles are disproportionately transmitted to the next generation. Evidence of such 'transmission distortion' (TD) in humans remains equivocal in part because scans of human pedigrees have been under-powered to detect small effects. After applying rhapsodi to the sperm data and scanning for evidence of TD, our results exhibited close concordance with binomial expectations under balanced transmission. Together, our work demonstrates that rhapsodi can facilitate novel uses of inferred genotype data and meiotic recombination events, while offering a powerful quantitative framework for testing for TD in other cohorts and study systems.

Novel Expressed Sequence Tag-Derived and Other Genomic Simple Sequence Repeat Markers Revealed Genetic Diversity in Ethiopian Finger Millet Landrace Populations and Cultivars

Finger millet (Eleusine coracana (L.) Geartn.) is a self-pollinating amphidiploid crop cultivated with minimal input for food and feed, as well as a source of income for small-scale farmers. To efficiently assess its genetic diversity for conservation and use in breeding programs, polymorphic DNA markers that represent its complex tetraploid genome have to be developed and used. In this study, 13 new expressed sequence tag-derived simple sequence repeat (EST-SSR) markers were developed based on publicly available finger millet ESTs. Using 10 polymorphic SSR markers (3 genomic and 7 novel EST-derived), the genetic diversity of 55 landrace accessions and 5 cultivars of finger millet representing its major growing areas in Ethiopia was assessed. In total, 26 alleles were detected across the 10 loci, and the average observed number of alleles per locus was 5.6. The polymorphic information content (PIC) of the loci ranged from 0.045 (Elco-48) to 0.71 (UGEP-66). The level of genetic diversity did not differ much between the accessions with the mean gene diversity estimates ranging only from 0.44 (accession 216054) to 0.68 (accession 237443). Similarly, a narrow range of variation was recorded at the level of regional states ranging from 0.54 (Oromia) to 0.59 (Amhara and Tigray). Interestingly, the average gene diversity of the landrace accessions (0.57) was similar to that of the cultivars (0.58). The analysis of molecular variance (AMOVA) revealed significant genetic variation both within and among accessions. The variation among the accessions accounted for 18.8% of the total variation (F _ST = 0.19; P < 0.001). Similarly, significant genetic variation was obtained among the geographic regions, accounting for 6.9% of the total variation (P < 0.001). The results of the cluster, principal coordinate, and population structure analyses suggest a poor correlation between the genetic makeups of finger millet landrace populations and their geographic regions of origin, which in turn suggests strong gene flow between populations within and across geographic regions. This study contributed novel EST-SSR markers for their various applications, and those that were monomorphic should be tested in more diverse finger millet genetic resources.

Comparative transcriptomics provides a strategy for phylogenetic analysis and SSR marker development in Chaenomeles

The genus Chaenomeles has long been considered an important ornamental, herbal and cash crop and is widely cultivated in East Asia. Traditional studies of Chaenomeles mainly focus on evolutionary relationships at the phenotypic level. In this study, we conducted RNA-seq on 10 Chaenomeles germplasms supplemented with one outgroup species, Docynia delavayi (D. delavayi), on the Illumina HiSeq2500 platform. After de novo assemblies, we generated from 40,084 to 49,571 unigenes for each germplasm. After pairwise comparison of the orthologous sequences, 9,659 orthologues within the 11 germplasms were obtained, with 6,154 orthologous genes identified as single-copy genes. The phylogenetic tree was visualized to reveal evolutionary relationships for these 11 germplasms. GO and KEGG analyses were performed for these common single-copy genes to compare their functional similarities and differences. Selective pressure analysis based on 6,154 common single-copy genes revealed that 45 genes were under positive selection. Most of these genes are involved in building the plant disease defence system. A total of 292 genes containing simple sequence repeats (SSRs) were used to develop SSR markers and compare their functions in secondary metabolism pathways. Finally, 10 primers were chosen as SSR marker candidates for Chaenomeles germplasms by comprehensive standards. Our research provides a new methodology and reference for future related research in Chaenomeles and is also useful for improvement, breeding and selection projects in other related species.

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 Locus Analysis and Identification of Candidate Genes for Micronaire in an Interspecific Backcross Inbred Line Population of Gossypium hirsutum × Gossypium barbadense

Cotton is the most important fiber crop and provides indispensable natural fibers for the textile industry. Micronaire (MIC) is determined by fiber fineness and maturity and is an important component of fiber quality. Gossypium barbadense L. possesses long, strong and fine fibers, while upland cotton (Gossypium hirsutum L.) is high yielding with high MIC and widely cultivated worldwide. To identify quantitative trait loci (QTLs) and candidate genes for MIC in G. barbadense, a population of 250 backcross inbred lines (BILs), developed from an interspecific cross of upland cotton CRI36 × Egyptian cotton (G. barbadense) Hai7124, was evaluated in 9 replicated field tests. Based on a high-density genetic map with 7709 genotyping-by-sequencing (GBS)-based single-nucleotide polymorphism (SNP) markers, 25 MIC QTLs were identified, including 12 previously described QTLs and 13 new QTLs. Importantly, two stable MIC QTLs (qMIC-D03-2 on D03 and qMIC-D08-1 on D08) were identified. Of a total of 338 genes identified within the two QTL regions, eight candidate genes with differential expression between TM-1 and Hai7124 were identified. Our research provides valuable information for improving MIC in cotton breeding.

Improved reconstruction and comparative analysis of chromosome 12 to rectify Mis-assemblies in Gossypium arboreum

BACKGROUND

Genome sequencing technologies have been improved at an exponential pace but precise chromosome-scale genome assembly still remains a great challenge. The draft genome of cultivated G. arboreum was sequenced and assembled with shotgun sequencing approach, however, it contains several misassemblies. To address this issue, we generated an improved reassembly of G. arboreum chromosome 12 using genetic mapping and reference-assisted approaches and evaluated this reconstruction by comparing with homologous chromosomes of G. raimondii and G. hirsutum.

RESULTS

In this study, we generated a high quality assembly of the 94.64 Mb length of G. arboreum chromosome 12 (A_A12) which comprised of 144 scaffolds and contained 3361 protein coding genes. Evaluation of results using syntenic and collinear analysis of reconstructed G. arboreum chromosome A_A12 with its homologous chromosomes of G. raimondii (D_D08) and G. hirsutum (AD_A12 and AD_D12) confirmed the significant improved quality of current reassembly as compared to previous one. We found major misassemblies in previously assembled chromosome 12 (A_Ca9) of G. arboreum particularly in anchoring and orienting of scaffolds into a pseudo-chromosome. Further, homologous chromosomes 12 of G. raimondii (D_D08) and G. arboreum (A_A12) contained almost equal number of transcription factor (TF) related genes, and showed good collinear relationship with each other. As well, a higher rate of gene loss was found in corresponding homologous chromosomes of tetraploid (AD_A12 and AD_D12) than diploid (A_A12 and D_D08) cotton, signifying that gene loss is likely a continuing process in chromosomal evolution of tetraploid cotton.

CONCLUSION

This study offers a more accurate strategy to correct misassemblies in sequenced draft genomes of cotton which will provide further insights towards its genome organization.

Genetic dissection of an allotetraploid interspecific CSSLs guides interspecific genetics and breeding in cotton

Background

The low genetic diversity of Upland cotton limits the potential for genetic improvement. Making full use of the genetic resources of Sea-island cotton will facilitate genetic improvement of widely cultivated Upland cotton varieties. The chromosome segments substitution lines (CSSLs) provide an ideal strategy for mapping quantitative trait loci (QTL) in interspecific hybridization.

Results

In this study, a CSSL population was developed by PCR-based markers assisted selection (MAS), derived from the crossing and backcrossing of Gossypium hirsutum (Gh) and G. barbadense (Gb), firstly. Then, by whole genome re-sequencing, 11,653,661 high-quality single nucleotide polymorphisms (SNPs) were identified which ultimately constructed 1211 recombination chromosome introgression segments from Gb. The sequencing-based physical map provided more accurate introgressions than the PCR-based markers. By exploiting CSSLs with mutant morphological traits, the genes responding for leaf shape and fuzz-less mutation in the Gb were identified. Based on a high-resolution recombination bin map to uncover genetic loci determining the phenotypic variance between Gh and Gb, 64 QTLs were identified for 14 agronomic traits with an interval length of 158 kb to 27 Mb. Surprisingly, multiple alleles of Gb showed extremely high value in enhancing cottonseed oil content (SOC).

Conclusions

This study provides guidance for studying interspecific inheritance, especially breeding researchers, for future studies using the traditional PCR-based molecular markers and high-throughput re-sequencing technology in the study of CSSLs. Available resources include candidate position for controlling cotton quality and quantitative traits, and excellent breeding materials. Collectively, our results provide insights into the genetic effects of Gb alleles on the Gh, and provide guidance for the utilization of Gb alleles in interspecific breeding.

Identification of candidate genes for key fibre-related QTLs and derivation of favourable alleles in Gossypium hirsutum recombinant inbred lines with G. barbadense introgressions

Fine mapping QTLs and identifying candidate genes for cotton fibre-quality and yield traits would be beneficial to cotton breeding. Here, we constructed a high-density genetic map by specific-locus amplified fragment sequencing (SLAF-seq) to identify QTLs associated with fibre-quality and yield traits using 239 recombinant inbred lines (RILs), which was developed from LMY22 (a high-yield Gossypium hirsutumL. cultivar) × LY343 (a superior fibre-quality germplasm with G. barbadenseL. introgressions). The genetic map spanned 3426.57 cM, including 3556 SLAF-based SNPs and 199 SSR marker loci. A total of 104 QTLs, including 67 QTLs for fibre quality and 37 QTLs for yield traits, were identified with phenotypic data collected from 7 environments. Among these, 66 QTLs were co-located in 19 QTL clusters on 12 chromosomes, and 24 QTLs were detected in three or more environments and determined to be stable. We also investigated the genomic components of LY343 and their contributions to fibre-related traits by deep sequencing the whole genome of LY343, and we found that genomic components from G. hirsutum races (which entered LY343 via its G. barbadense parent) contributed more favourable alleles than those from G. barbadense. We further identified six putative candidate genes for stable QTLs, including Gh_A03G1147 (GhPEL6), Gh_D07G1598 (GhCSLC6) and Gh_D13G1921 (GhTBL5) for fibre-length QTLs and Gh_D03G0919 (GhCOBL4), Gh_D09G1659 (GhMYB4) and Gh_D09G1690 (GhMYB85) for lint-percentage QTLs. Our results provide comprehensive insight into the genetic basis of the formation of fibre-related traits and would be helpful for cloning fibre-development-related genes as well as for marker-assisted genetic improvement in cotton.

Detection and application of genome-wide variations in peach for association and genetic relationship analysis

BACKGROUND

Peach (Prunus persica L.) is a diploid species and model plant of the Rosaceae family. In the past decade, significant progress has been made in peach genetic research via DNA markers, but the number of these markers remains limited.

RESULTS

In this study, we performed a genome-wide DNA markers detection based on sequencing data of six distantly related peach accessions. A total of 650,693~1,053,547 single nucleotide polymorphisms (SNPs), 114,227~178,968 small insertion/deletions (InDels), 8386~12,298 structure variants (SVs), 2111~2581 copy number variants (CNVs) and 229,357~346,940 simple sequence repeats (SSRs) were detected and annotated. To demonstrate the application of DNA markers, 944 SNPs were filtered for association study of fruit ripening time and 15 highly polymorphic SSRs were selected to analyze the genetic relationship among 221 accessions.

CONCLUSIONS

The results showed that the use of high-throughput sequencing to develop DNA markers is fast and effective. Comprehensive identification of DNA markers, including SVs and SSRs, would be of benefit to genetic diversity evaluation, genetic mapping, and molecular breeding of peach.

Enhancing Upland cotton for drought resilience, productivity, and fiber quality: comparative evaluation and genetic dissection

To provision the world sustainably, modern society must increase overall crop production, while conserving and preserving natural resources. Producing more with diminishing water resources is an especially daunting endeavor. Toward the goal of genetically improving drought resilience of cultivated Upland cotton (Gossypium hirsutum L.), this study addresses the genetics of differential yield components referred to as productivity and fiber quality traits under regular-water versus low-water (LW) field conditions. We used ten traits to assess water stress deficit, which included six productivity and four fiber quality traits on two recombinant inbred line (RIL) populations from reciprocally crossed cultivars, Phytogen 72 and Stoneville 474. To facilitate genetic inferences, we genotyped RILs with the CottonSNP63K array, assembled high-density linkage maps of over 7000 SNPs and then analyzed quantitative trait variations. Analysis of variance revealed significant differences for all traits (p < 0.05) in these RIL populations. Although the LW irrigation regime significantly reduced all traits, except lint percent, the RILs exhibited a broad phenotypic spectrum of heritable differences across the water regimes. Transgressive segregation occurred among the RILs, suggesting the possibility of genetic gain through phenotypic selection for drought resilience and perhaps through marker-based selection. Analyses revealed more than 150 quantitative trait loci (QTLs) associated with productivity and fiber quality traits (p < 0.005) on different genomic regions of the cotton genome. The multiple-QTL models analysis with LOD > 3.0 detected 21 QTLs associated with productivity and 22 QTLs associated with fiber quality. For fiber traits, strong clustering and QTL associations occurred in c08 and its homolog c24 as well as c10, c14, and c21. Using contemporary genome sequence assemblies and bioinformatically related information, the identification of genomic regions associated with responses to plant stress/drought elevates the possibility of using marker-assisted and omics-based selection to enhance breeding for drought resilient cultivars and identifying candidate genes and networks. RILs with different responses to drought indicated that it is possible to maintain high fiber quality under LW conditions or reduce the of LW impact on quality. The heritable variation among elite bi-parental RILs for productivity and quality under field drought conditions, and their association of QTLs, and thus specific genomic regions, indicate opportunities for breeding-based gains in water resource conservation, i.e., enhancing cotton's agricultural sustainability.

Reference genome sequences of two cultivated allotetraploid cottons, Gossypium hirsutum and Gossypium barbadense

Allotetraploid cotton species (Gossypium hirsutum and Gossypium barbadense) have long been cultivated worldwide for natural renewable textile fibers. The draft genome sequences of both species are available but they are highly fragmented and incomplete^1-4. Here we report reference-grade genome assemblies and annotations for G. hirsutum accession Texas Marker-1 (TM-1) and G. barbadense accession 3-79 by integrating single-molecule real-time sequencing, BioNano optical mapping and high-throughput chromosome conformation capture techniques. Compared with previous assembled draft genomes^1,3, these genome sequences show considerable improvements in contiguity and completeness for regions with high content of repeats such as centromeres. Comparative genomics analyses identify extensive structural variations that probably occurred after polyploidization, highlighted by large paracentric/pericentric inversions in 14 chromosomes. We constructed an introgression line population to introduce favorable chromosome segments from G. barbadense to G. hirsutum, allowing us to identify 13 quantitative trait loci associated with superior fiber quality. These resources will accelerate evolutionary and functional genomic studies in cotton and inform future breeding programs for fiber improvement.

SSR-Linkage map of interspecific populations derived from Gossypium trilobum and Gossypium thurberi and determination of genes harbored within the segregating distortion regions

Wild cotton species have significant agronomic traits that can be introgressed into elite cultivated varieties. The use of a genetic map is important in exploring, identification and mining genes which carry significant traits. In this study, 188 F₂mapping individuals were developed from Gossypium thurberi (female) and Gossypium trilobum (male), and were genotyped by using simple sequence repeat (SSR) markers. A total of 12,560 simple sequence repeat (SSR) markers, developed by Southwest University, thus coded SWU were screened out of which only 994 were found to be polymorphic, and 849 markers were linked in all the 13 chromosomes. The map had a length of 1,012.458 cM with an average marker distance of 1.193 cM. Segregation distortion regions (SDRs) were observed on Chr01, Chr02, Chr06, Chr07 Chr09, Chr10 and Chr11 with a large proportion of the SDR regions segregating towards the heterozygous allele. There was good syntenic block formation that revealed good collinearity between the genetic and physical map of G. raimondii, compared to the Dt_sub genome of the G. hirsutum and G. barbadense. A total of 2,496 genes were mined within the SSR related regions. The proteins encoding the mined genes within the SDR had varied physiochemical properties; their molecular weights ranged from 6.586 to 252.737 kDa, charge range of -39.5 to 52, grand hydropathy value (GRAVY) of -1.177 to 0.936 and isoelectric (pI) value of 4.087 to 12.206. The low GRAVY values detected showed that the proteins encoding these genes were hydrophilic in nature, a property common among the stress responsive genes. The RNA sequence analysis revealed more of the genes were highly upregulated in various stages of fiber development for instance; Gorai.002G241300 was highly up regulated at 5, 10, 20 and 25 day post anthesis (DPA). Validation through RT-qPCR further revealed that these genes mined within the SDR regions might be playing a significant role under fiber development stages, therefore we infer that Gorai.007G347600 (TFCA), Gorai.012G141600 (FOLB1), Gorai.006G024500 (NMD3), Gorai.002G229900 (LST8) and Gorai.002G235200 (NSA2) are significantly important in fiber development and in turn the quality, and further researches needed to be done to elucidate their exact roles in the fiber development process. The construction of the genetic map between the two wild species paves away for the mapping of quantitative trait loci (QTLs) since the average distance between the markers is small, and mining of genes on the SSR regions will provide an insight in identifying key genes that can be introgressed into the cultivated cotton cultivars.

A high-density genetic map of extra-long staple cotton (Gossypium barbadense) constructed using genotyping-by-sequencing based single nucleotide polymorphic markers and identification of fiber traits-related QTL in a recombinant inbred line population

Background

Gossypium barbadense (Sea Island, Egyptian or Pima cotton) cotton has high fiber quality, however, few studies have investigated the genetic basis of its traits using molecular markers. Genome complexity reduction approaches such as genotyping-by-sequencing have been utilized to develop abundant markers for the construction of high-density genetic maps to locate quantitative trait loci (QTLs).

Results

The Chinese G. barbadense cultivar 5917 and American Pima S-7 were used to develop a recombinant inbred line (RIL) population with 143 lines. The 143 RILs together with their parents were tested in three replicated field tests for lint yield traits (boll weight and lint percentage) and fiber quality traits (fiber length, fiber elongation, fiber strength, fiber uniformity and micronaire) and then genotyped using GBS to develop single-nucleotide polymorphism (SNP) markers. A high-density genetic map with 26 linkage groups (LGs) was constructed using 3557 GBS SNPs spanning a total genetic distance of 3076.23 cM at an average density of 1.09 cM between adjacent markers. A total of 42 QTLs were identified, including 24 QTLs on 12 LGs for fiber quality and 18 QTLs on 7 LGs for lint yield traits, with LG1 (9 QTLs), LG10 (7 QTLs) and LG14 (6 QTLs) carrying more QTLs. Common QTLs for the same traits and overlapping QTLs for different traits were detected. Each individual QTLs explained 0.97 to 20.7% of the phenotypic variation.

Conclusions

This study represents one of the first genetic mapping studies on the fiber quality and lint yield traits in a RIL population of G. barbadense using GBS-SNPs. The results provide important information for the subsequent fine mapping of QTLs and the prediction of candidate genes towards map-based cloning and marker-assisted selection in cotton.

Quantitative Trait Loci Mapping of Multiple Independent Loci for Resistance to Fusarium oxysporum f. sp. vasinfectum Races 1 and 4 in an Interspecific Cotton Population

Fusarium wilt, caused by the soilborne fungal pathogen Fusarium oxysporum f. sp. vasinfectum, is a vascular disease of cotton (Gossypium spp.). F. oxysporum f. sp. vasinfectum race 1 (FOV1) causes major plant injury and yield loss in G. hirsutum cultivars with coinfection with root-knot nematode (Meloidogyne incognita), while F. oxysporum f. sp. vasinfectum race 4 (FOV4) causes plant damage without nematode coinfection in G. hirsutum and in G. barbadense cultivars. Quantitative trait loci (QTL) analysis of the interspecific cross G. barbadense Pima S-7 × G. hirsutum Acala NemX revealed separate multiple loci determining resistance to FOV1 and FOV4, confirming that race specificity occurs in F. oxysporum f. sp. vasinfectum. Based on the area under the disease progress stairs, six major QTLs on chromosomes (Chrs) 1, 2, 12, 15 (2), and 21 contributing 7 to 15% to FOV1 resistance and two major QTLs on Chrs 14 and 17 contributing 12 to 33% to FOV4 resistance were identified. Minor-effect QTLs contributing to resistance to both FOV1 and FOV4 were also identified. These results define and establish a pathosystem of race-specific resistance under polygenic control. This research also validates the importance of previously reported markers and chromosome regions and adds new information for the location of F. oxysporum f. sp. vasinfectum resistance genes. Some F₈ recombinant inbred lines have resistance to both FOV1 and FOV4 and also to root-knot nematode, providing multiple resistance sources for breeding.

Genetic Map Construction and Fiber Quality QTL Mapping Using the CottonSNP80K Array in Upland Cotton

Cotton fiber quality traits are controlled by multiple quantitative trait loci (QTL), and the improvement of these traits requires extensive germplasm. Herein, an Upland cotton cultivar from America, Acala Maxxa, was crossed with a local high fiber quality cultivar, Yumian 1, and 180 recombinant inbred lines (RILs) were obtained. In order to dissect the genetic basis of fiber quality differences between these parents, a genetic map containing 12116 SNP markers was constructed using the CottonSNP80K assay, which covered 3741.81 cM with an average distance of 0.31 cM between markers. Based on the genetic map and growouts in three environments, we detected a total of 104 QTL controlling fiber quality traits. Among these QTL, 25 were detected in all three environments and 35 in two environments. Meanwhile, 19 QTL clusters were also identified, and nine contained at least one stable QTL (detected in three environments for a given trait). These stable QTL or QTL clusters are priorities for fine mapping, identifying candidate genes, elaborating molecular mechanisms of fiber development, and application in cotton breeding programs by marker-assisted selection (MAS).

QTL delineation for five fiber quality traits based on an intra-specific Gossypium hirsutum L. recombinant inbred line population

Gossypium hirsutum L. is the most important fiber crop worldwide and contributes to more than 95% of global cotton production. Marker-assisted selection (MAS) is an effective approach for improving fiber quality, and quantitative trait loci (QTL) mapping of fiber quality traits is important for cotton breeding. In this study, a permanent intra-specific recombinant inbred line (RIL) population containing 137 families was used for fiber quality testing. Based on a previously reported high-density genetic map with an average marker distance of 0.63 cM, 186 additive QTLs were obtained for five fiber quality traits over five consecutive years, including 39 for fiber length (FL), 36 for fiber strength (FS), 50 for fiber uniformity (FU), 33 for micronaire (MC) and 28 for fiber elongation (FE). Three stable QTLs, qMC-A4-1, qMC-D2-3 and qFS-D9-1, were detected in four datasets, and another eight stable QTLs, qMC-A4-2, qMC-D11-2, qFU-A9-1, qFU-A10-4, qFS-D11-1, qFL-D9-2, qFL-D11-1 and qFE-A3-2, were detected in three datasets. The annotated genes in these 11 stable QTLs were collected, and these genes included many transcription factors with functions during fiber development. 33 QTL coincidence regions were found, and these involved nearly half of the total QTLs. Four chromosome regions containing at least 6 QTLs were promising for fine mapping. In addition, 41 pairs of epistatic QTLs (e-QTLs) were screened, including 6 for FL, 30 for FS, 2 for FU and 3 for MC. The identification of stable QTLs adds valuable information for further QTL fine mapping and gene positional cloning for fiber quality genetic detection and provides useful markers for further molecular breeding in enhancing fiber quality.

QTL Mapping of Fiber Quality and Yield-Related Traits in an Intra-Specific Upland Cotton Using Genotype by Sequencing (GBS)

Fiber quality and yield improvement are crucial for cotton domestication and breeding. With the transformation in spinning techniques and multiplicity needs, the development of cotton fiber quality and yield is of great importance. A genetic map of 5178 Single Nucleotide Polymorphism (SNP) markers were generated using 277 F2:3 population, from an intra-specific cross between two upland cotton accessions, CCRI35 a high fiber quality as female and Nan Dan Ba Di Da Hua (NH), with good yield properties as male parent. The map spanned 4768.098 cM with an average distance of 0.92 cM. A total of 110 Quantitative Traits Loci (QTLs) were identified for 11 traits, but only 30 QTLs were consistent in at least two environments. The highest percentage of phenotypic variance explained by a single QTL was 15.45%. Two major cluster regions were found, cluster 1 (chromosome17-D03) and cluster 2 (chromosome26-D12). Five candidate genes were identified in the two QTL cluster regions. Based on GO functional annotation, all the genes were highly correlated with fiber development, with functions such as protein kinase and phosphorylation. The five genes were associated with various fiber traits as follows: Gh_D03G0889 linked to qFM-D03_cb, Gh_D12G0093, Gh_D12G0410, Gh_D12G0435 associated with qFS-D12_cb and Gh_D12G0969 linked to qFY-D12_cb. Further structural annotation and fine mapping is needed to determine the specific role played by the five identified genes in fiber quality and yield related pathway.

Insights Into Upland Cotton (Gossypium hirsutum L.) Genetic Recombination Based on 3 High-Density Single-Nucleotide Polymorphism and a Consensus Map Developed Independently With Common Parents

High-density linkage maps are vital to supporting the correct placement of scaffolds and gene sequences on chromosomes and fundamental to contemporary organismal research and scientific approaches to genetic improvement, especially in paleopolyploids with exceptionally complex genomes, eg, upland cotton (Gossypium hirsutum L., "2n = 52"). Three independently developed intraspecific upland mapping populations were analyzed to generate 3 high-density genetic linkage single-nucleotide polymorphism (SNP) maps and a consensus map using the CottonSNP63K array. The populations consisted of a previously reported F2, a recombinant inbred line (RIL), and reciprocal RIL population, from "Phytogen 72" and "Stoneville 474" cultivars. The cluster file provided 7417 genotyped SNP markers, resulting in 26 linkage groups corresponding to the 26 chromosomes (c) of the allotetraploid upland cotton (AD)1 arisen from the merging of 2 genomes ("A" Old World and "D" New World). Patterns of chromosome-specific recombination were largely consistent across mapping populations. The high-density genetic consensus map included 7244 SNP markers that spanned 3538 cM and comprised 3824 SNP bins, of which 1783 and 2041 were in the At and Dt subgenomes with 1825 and 1713 cM map lengths, respectively. Subgenome average distances were nearly identical, indicating that subgenomic differences in bin number arose due to the high numbers of SNPs on the Dt subgenome. Examination of expected recombination frequency or crossovers (COs) on the chromosomes within each population of the 2 subgenomes revealed that COs were also not affected by the SNPs or SNP bin number in these subgenomes. Comparative alignment analyses identified historical ancestral At-subgenomic translocations of c02 and c03, as well as of c04 and c05. The consensus map SNP sequences aligned with high congruency to the NBI assembly of Gossypium hirsutum. However, the genomic comparisons revealed evidence of additional unconfirmed possible duplications, inversions and translocations, and unbalance SNP sequence homology or SNP sequence/loci genomic dominance, or homeolog loci bias of the upland tetraploid At and Dt subgenomes. The alignments indicated that 364 SNP-associated previously unintegrated scaffolds can be placed in pseudochromosomes of the NBI G hirsutum assembly. This is the first intraspecific SNP genetic linkage consensus map assembled in G hirsutum with a core of reproducible mendelian SNP markers assayed on different populations and it provides further knowledge of chromosome arrangement of genic and nongenic SNPs. Together, the consensus map and RIL populations provide a synergistically useful platform for localizing and identifying agronomically important loci for improvement of the cotton crop.

QTL Analysis of Transgressive Nematode Resistance in Tetraploid Cotton Reveals Complex Interactions in Chromosome 11 Regions

Transgressive segregation in cotton (Gossypium spp.) provides an important approach to enhance resistance to the major pest root-knot nematode (RKN) Meloidogyne incognita. Our previous studies reported transgressive RKN resistance in an intraspecific Gossypium hirsutum resistant NemX × susceptible SJ-2 recombinant inbred line (RIL) population and early generations of interspecific cross Gossypium barbadense (susceptible Pima S-7) × G. hirsutum (NemX). However, the underlying functional mechanisms for this phenomenon are not known. In this study, the region of RKN resistance gene rkn1 on chromosome (Chr) 11 and its homoeologous Chr 21 was fine mapped with G. raimondii D₅ genome reference sequence. Transgressive resistance was found in the later generation of a new RIL population F_2:7 (Pima S-7 × NemX) and one interspecific F₂ (susceptible Pima S-7 × susceptible SJ-2). QTL analysis revealed similar contributions to root-galling and egg-production resistance phenotypes associated with SSR marker CIR316 linked to resistance gene rkn1 in NemX on Chr 11 in all seven populations analyzed. In testcross NemX × F₁ (Pima S-7 × SJ-2) marker allele CIR069-271 from Pima S-7 linked to CIR316 contributed 63% of resistance to galling phenotype in the presence of rkn1. Similarly, in RIL population F_2:8 (NemX × SJ-2), SJ-2 markers closely linked to CIR316 contributed up to 82% of resistance to root-galling. These results were confirmed in BC₁F₁ SJ-2 × F₁ (NemX × SJ-2), F₂ (NemX × SJ-2), and F₂ (Pima S-7 × SJ-2) populations in which up to 44, 36, and 15% contribution in resistance to galling was found, respectively. Transgressive segregation for resistance was universal in all intra- and inter-specific populations, although stronger transgressive resistance occurred in later than in early generations in the intraspecific cross compared with the interspecific cross. Transgressive effects on progeny from susceptible parents are possibly provided in the rkn1 resistance region of chromosome 11 by tandemly arrayed allele (TAA) or gene (TAG) interactions contributing to transgressive resistance. Complex TAA and TAG recombination and interactions in the rkn1 resistance region provide three genes and a model to study disease and transgressive resistance in polyploid plants, and novel genotypes for plant breeding.

Chromosome structural variation of two cultivated tetraploid cottons and their ancestral diploid species based on a new high-density genetic map

A high-resolution genetic map is a useful tool for assaying genomic structural variation and clarifying the evolution of polyploid cotton. A total of 36956 SSRs, including 11289 released in previous studies and 25567 which were newly developed based on the genome sequences of G. arboreum and G. raimondii, were utilized to construct a new genetic map. The new high-density genetic map includes 6009 loci and spanned 3863.97 cM with an average distance of 0.64 cM between consecutive markers. Four inversions (one between Chr08 and Chr24, one between Chr09 and Chr23 and two between Chr10 and Chr20) were identified by homology analysis. Comparative genomic analysis between genetic map and two diploid cottons showed that structural variations between the A genome and At subgenome are more extensive than between D genome and Dt subgenome. A total of 17 inversions, seven simple translocations and two reciprocal translocations were identified between genetic map and G. raimondii. Good colinearity was revealed between the corresponding chromosomes of tetraploid G. hirsutum and G. barbadense genomes, but a total of 16 inversions were detected between them. These results will accelerate the process of evolution analysis of Gossipium genus.

Enriching an intraspecific genetic map and identifying QTL for fiber quality and yield component traits across multiple environments in Upland cotton (Gossypium hirsutum L.)

Cotton is a significant commercial crop that plays an indispensable role in many domains. Constructing high-density genetic maps and identifying stable quantitative trait locus (QTL) controlling agronomic traits are necessary prerequisites for marker-assisted selection (MAS). A total of 14,899 SSR primer pairs designed from the genome sequence of G. raimondii were screened for polymorphic markers between mapping parents CCRI 35 and Yumian 1, and 712 SSR markers showing polymorphism were used to genotype 180 lines from a (CCRI 35 × Yumian 1) recombinant inbred line (RIL) population. Genetic linkage analysis was conducted on 726 loci obtained from the 712 polymorphic SSR markers, along with 1379 SSR loci obtained in our previous study, and a high-density genetic map with 2051 loci was constructed, which spanned 3508.29 cM with an average distance of 1.71 cM between adjacent markers. Marker orders on the linkage map are highly consistent with the corresponding physical orders on a G. hirsutum genome sequence. Based on fiber quality and yield component trait data collected from six environments, 113 QTLs were identified through two analytical methods. Among these 113 QTLs, 50 were considered stable (detected in multiple environments or for which phenotypic variance explained by additive effect was greater than environment effect), and 18 of these 50 were identified with stability by both methods. These 18 QTLs, including eleven for fiber quality and seven for yield component traits, could be priorities for MAS.

Extrapolative microRNA precursor based SSR mining from tea EST database in respect to agronomic traits

Tea (Camellia sinensis, (L.) Kuntze) is considered as most popular drink across the world and it is widely consumed beverage for its several health-benefit characteristics. These positive traits primarily rely on its regulatory networks of different metabolic pathways. Development of microsatellite markers from the conserved genomic regions are being worthwhile for reviewing the genetic diversity of closely related species or self-pollinated species. Although several SSR markers have been reported, in tea, the trait-specific Simple Sequence Repeat (SSR) markers, leading to be useful in marker assisted breeding technique, are yet to be identified. Micro RNAs are short, non-coding RNA molecules, involved in post transcriptional mode of gene regulation and thus effects on related phenotype. Present study deals with identification of the microsatellite motifs within the reported and predicted miRNA precursors that are effectively followed by designing of primers from SSR flanking regions in order to PCR validation. In addition to the earlier reports, two new miRNAs are predicting here from tea expressed tag sequence database. Furthermore, 18 SSR motifs are found to be in 13 of all 33 predicted miRNAs. Trinucleotide motifs are most abundant among all followed by dinucleotides. Since, miRNA based SSR markers are evidenced to have significant role on genetic fingerprinting study, these outcomes would pave the way in developing novel markers for tagging tea specific agronomic traits as well as substantiating non-conventional breeding program.

Fine mapping and candidate gene analysis of qFL-chr1, a fiber length QTL in cotton

A fiber length QTL, qFL-chr1, was fine mapped to a 0.9 cM interval of cotton chromosome 1. Two positional candidate genes showed positive correlation between gene expression level and fiber length. Prior analysis of a backcross-self mapping population derived from a cross between Gossypium hirsutum L. and G. barbadense L. revealed a QTL on chromosome 1 associated with increased fiber length (qFL-chr1), which was confirmed in three independent populations of near-isogenic introgression lines (NIILs). Here, a single NIIL, R01-40-08, was used to develop a large population segregating for the target region. Twenty-two PCR-based polymorphic markers used to genotype 1672 BC₄F₂ plants identified 432 recombinants containing breakpoints in the target region. Substitution mapping using 141 informative recombinants narrowed the position of qFL-chr1 to a 1.0-cM interval between SSR markers MUSS084 and CIR018. To exclude possible effects of non-target introgressions on fiber length, different heterozygous BC₄F₃ plants introgressed between SSR markers NAU3384 and CGR5144 were selected to develop sub-NILs. The qFL-chr1 was further mapped at 0.9-cM interval between MUSS422 and CIR018 by comparisons of sub-NIL phenotype, and increased fiber length by ~1 mm. The 2.38-Mb region between MUSS422 and CIR018 in G. barbadense contained 19 annotated genes. Expression levels of two of these genes, GOBAR07705 (encoding 1-aminocyclopropane-1-carboxylate synthase) and GOBAR25992 (encoding amino acid permease), were positively correlated with fiber length in a small F₂ population, supporting these genes as candidates for qFL-chr1.

Microdissection of the Ah01 chromosome in upland cotton and microcloning of resistance gene anologs from the single chromosome

BACKGROUND

Chromosome microdissection is one of the most important techniques in molecular cytogenetic research. Cotton (Gossypium Linnaeus, 1753) is the main natural fiber crop in the world. The resistance gene analog (RGA) cloning after its single chromosome microdissection can greatly promote cotton genome research and breeding.

RESULTS

Using the linker adaptor PCR (LA-PCR) with the primers of rice disease-resistance homologues, three nucleotide sequences PS016 (KU051681), PS054 (KU051682), and PS157 (KU051680) were obtained from the chromosome A_h01 of upland cotton (cv. TM-1). The Blast results showed that the three sequences are the nucleotide binding site-leucine rich repeat (NBS-LRR) type RGAs. Clustering results indicated that they are homologous to these published RGAs. Thus, the three RGAs can definitely be confirmed as NBS-LRR class of RGAs in upland cotton.

CONCLUSIONS

Using single chromosome microdissection technique, DNA libraries containing cotton RGAs were obtained. This technique can promote cotton gene cloning, marker development and even the improvement of cotton genome research and breeding.

Identification of genes related to salt stress tolerance using intron-length polymorphic markers, association mapping and virus-induced gene silencing in cotton

Intron length polymorphisms (ILPs), a type of gene-based functional marker, could themselves be related to the particular traits. Here, we developed a genome-wide cotton ILPs based on orthologs annotation from two sequenced diploid species, A-genome Gossypium arboreum and D-genome G. raimondii. We identified 10,180 putative ILP markers from 5,021 orthologous genes. Among these, 535 ILP markers from 9 gene families related to stress were selected for experimental verification. Polymorphic rates were 72.71% between G. arboreum and G. raimondii and 36.45% between G. hirsutum acc. TM-1 and G. barbadense cv. Hai7124. Furthermore, 14 polymorphic ILP markers were detected in 264 G. hirsutum accessions. Coupled with previous simple sequence repeats (SSRs) evaluations and salt tolerance assays from the same individuals, we found a total of 25 marker-trait associations involved in nine ILPs. The nine genes, temporally named as C1 to C9, showed the various expressions in different organs and tissues, and five genes (C3, C4, C5, C7 and C9) were significantly upregulated after salt treatment. We verified that the five genes play important roles in salt tolerance. Particularly, silencing of C4 (encodes WRKY DNA-binding protein) and C9 (encodes Mitogen-activated protein kinase) can significantly enhance cotton susceptibility to salt stress.

Uncovering SNP and indel variations of tetraploid cottons by SLAF-seq

Background

Cotton (Gossypium spp.), as the world’s most utilized textile fibre source, is an important, economically valuable crop worldwide. Understanding the genomic variation of tetraploid cotton species is important for exploitation of the excellent characteristics of wild cotton and for improving the diversity of cotton in breeding. However, the discovery of DNA polymorphisms in tetraploid cotton genomes has lagged behind other important crops.

Results

A total of 111,795,823 reads, 467,735 specific length amplified fragment (SLAF) tags and 139,176 high-quality DNA polymorphisms were identified using specific length amplified fragment sequencing (SLAF-seq), including 132,880 SNPs and 6,296 InDels between the reference genome (TM-1) and the five tetraploid cotton species. Intriguingly, gene ontology (GO) enrichment analysis revealed that a number of significant terms were related to reproduction in G. barbadense acc. 3–79. Based on the new data sets, we reconstructed phylogenetic trees that showed a high concordance to the phylogeny of diploid and polyploid cottons. A large amount of interspecific genetic variations were identified, and some of them were validated by the single-strand conformation polymorphism (SSCP) method, which will be applied in introgression genetics and breeding with G. hirsutum cv. Emian22 as the receptor and the other species as donors.

Conclusions

Using SLAF-seq, a large number of DNA polymorphisms were identified. The comprehensive analysis of DNA polymorphisms provided invaluable insights into the different tetraploid cotton species. More importantly, the identification of numerous interspecific genetic variations provides the basis and is very practical for future introgression breeding. The results presented herein provide a valuable genomic resource for new insights into the genetics and breeding of cotton.

Electronic supplementary material

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

High-density linkage map construction and QTL analysis for earliness-related traits in Gossypium hirsutum L

Background

Gossypium hirsutum L., or upland cotton, is an important renewable resource for textile fiber. To enhance understanding of the genetic basis of cotton earliness, we constructed an intra-specific recombinant inbred line population (RIL) containing 137 lines, and performed linkage map construction and quantitative trait locus (QTL) mapping.

Results

Using restriction-site associated DNA sequencing, a genetic map composed of 6,434 loci, including 6,295 single nucleotide polymorphisms and 139 simple sequence repeat loci, was developed from RIL population. This map spanned 4,071.98 cM, with an average distance of 0.63 cM between adjacent markers. A total of 247 QTLs for six earliness-related traits were detected in 6 consecutive years. In addition, 55 QTL coincidence regions representing more than 60 % of total QTLs were found on 22 chromosomes, which indicated that several earliness-related traits might be simultaneously improved. Fine-mapping of a 2-Mb region on chromosome D3 associated with five stable QTLs between Marker25958 and Marker25963 revealed that lines containing alleles derived from CCRI36 in this region exhibited smaller phenotypes and earlier maturity. One candidate gene (EMF2) was predicted and validated by quantitative real-time PCR in early-, medium- and late-maturing cultivars from 3- to 6-leaf stages, with highest expression level in early-maturing cultivar, CCRI74, lowest expression level in late-maturing cultivar, Bomian1.

Conclusions

We developed an SNP-based genetic map, and this map is the first high-density genetic map for short-season cotton and has the potential to provide deeper insights into earliness. Cotton earliness-related QTLs and QTL coincidence regions will provide useful materials for QTL fine mapping, gene positional cloning and MAS. And the gene, EMF2, is promising for further study.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3269-y) contains supplementary material, which is available to authorized users.

Genome-Wide SNP Linkage Mapping and QTL Analysis for Fiber Quality and Yield Traits in the Upland Cotton Recombinant Inbred Lines Population

It is of significance to discover genes related to fiber quality and yield traits and tightly linked markers for marker-assisted selection (MAS) in cotton breeding. In this study, 188 F8 recombinant inbred lines (RILs), derived from a intraspecific cross between HS46 and MARCABUCAG8US-1-88 were genotyped by the cotton 63K single nucleotide polymorphism (SNP) assay. Field trials were conducted in Sanya, Hainan Province, during the 2014-2015 cropping seasons under standard conditions. Results revealed significant differences (P < 0.05) among RILs, environments and replications for fiber quality and yield traits. Broad-sense heritabilities of all traits including fiber length, fiber uniformity, micronaire, fiber elongation, fiber strength, boll weight, and lint percentage ranged from 0.26 to 0.66. A 1784.28 cM (centimorgans) linkage map, harboring 2618 polymorphic SNP markers, was constructed, which had 0.68 cM per marker density. Seventy-one quantitative trait locus (QTLs) for fiber quality and yield traits were detected on 21 chromosomes, explaining 4.70∼32.28% phenotypic variance, in which 16 were identified as stable QTLs across two environments. Meanwhile, 12 certain regions were investigated to be involved in the control of one (hotspot) or more (cluster) traits, mainly focused on Chr05, Chr09, Chr10, Chr14, Chr19, and Chr20. Nineteen pairs of epistatic QTLs (e-QTLs) were identified, of which two pairs involved in two additive QTLs. These additive QTLs, e-QTLs, and QTL clusters were tightly linked to SNP markers, which may serve as target regions for map-based cloning, gene discovery, and MAS in cotton breeding.

A Genetic Map Between Gossypium hirsutum and the Brazilian Endemic G. mustelinum and Its Application to QTL Mapping

Among the seven tetraploid cotton species, little is known about transmission genetics and genome organization in Gossypium mustelinum, the species most distant from the source of most cultivated cotton, G. hirsutum. In this research, an F₂ population was developed from an interspecific cross between G. hirsutum and G. mustelinum (HM). A genetic linkage map was constructed mainly using simple sequence repeat (SSRs) and restriction fragment length polymorphism (RFLP) DNA markers. The arrangements of most genetic loci along the HM chromosomes were identical to those of other tetraploid cotton species. However, both major and minor structural rearrangements were also observed, for which we propose a parsimony-based model for structural divergence of tetraploid cottons from common ancestors. Sequences of mapped markers were used for alignment with the 26 scaffolds of the G. hirsutum draft genome, and showed high consistency. Quantitative trait locus (QTL) mapping of fiber elongation in advanced backcross populations derived from the same parents demonstrated the value of the HM map. The HM map will serve as a valuable resource for QTL mapping and introgression of G. mustelinum alleles into G. hirsutum, and help clarify evolutionary relationships between the tetraploid cotton genomes.

Genome-wide SSR-based association mapping for fiber quality in nation-wide upland cotton inbreed cultivars in China

Background

Since upland cotton was introduced into China during the 1920s–1950s, hundreds of inbreed cultivars have been developed. To explore the molecular diversity, population structure and elite alleles, 503 inbred cultivars developed in China and some foreign cultivars from the United States and the Soviet Union were collected and analyzed by 494 genome-wide SSRs (Simple Sequence Repeats).

Methods

Four hundred and ninety-four pairs of SSRs with high polymorphism and uniform distribution on 26 chromosomes were used to scan polymorphisms in 503 nation-wide upland cottons. The programming language R was used to make boxplots for the phenotypic traits in different environments. Molecular marker data and 6 fiber quality traits were analyzed by the method of MLM (mixed linear model) (P + G + Q + K) in the TASSEL software package on the basis of the population structure and linkage disequilibrium analysis. The loci of elite allelic variation and typical materials carrying elite alleles were identified based on phenotypic effect values.

Results

A total of 179 markers were polymorphic and generated 426 allele loci; the population based on molecular diversity was classified into seven subpopulations corresponding to pedigree origin, ecological and geographical distribution. The attenuation distance of linkage disequilibrium dropped significantly up to 0–5 cM. Association mapping for fiber quality showed that 216 marker loci were associated with fiber quality traits (P < 0.05) explaining 0.58 % ~ 5.12 % of the phenotypic variation, with an average of 2.70 %. Thirteen marker loci were coincident with other studies, and three were detected for the same trait. Seven quantitative trait loci were related to known genes in fiber development. Based on phenotypic effects, 48 typical materials that contained the elite allele loci related to fiber quality traits were identified and are widely used in practical breeding.

Conclusions

The molecular diversity and population structure of 503 nation-wide upland cottons in China were evaluated by 494 genome-wide SSRs, and association mapping for fiber quality revealed known and novel elite alleles. The molecular diversity provides a guide for parental mating in cotton breeding, and the association mapping results will aid in the fine-mapping genes related to fiber quality traits and facilitate further studies on candidate genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2662-x) contains supplementary material, which is available to authorized users.

Detection and validation of one stable fiber strength QTL on c9 in tetraploid cotton

Fiber strength is an essential trait of fiber property in cotton, and it is quantitatively inherited. Identification of stable quantitative trait loci (QTL) contributing to fiber strength would provide the key basis for marker-assisted selection (MAS) in cotton breeding. In this study, four interspecific hybridization populations were established with a common G. barbadense parent Pima 90-53 and two G. hirsutum parents (CCRI 8 and Handan 208), each of which had fiber strength characteristic. Based on the phenotypic data of fiber strength from seven environments, a stable QTL, qFS-c9-1, was detected and validated on c9 in a marker interval between SSR markers NAU2395 and NAU1092. The QTL explaining 14.4-17.9 % of the phenotypic variation was firstly detected in two populations (CCRI 8 × Pima 90-53, BC1F1 and BC1F2) and its derived lines in four environments. And it accounting for 12.1-14.8 % of the phenotypic variation was further confirmed in two populations (Handan 208 × Pima 90-53, BC1F1, and F2) under one environment. In silico mapping using three sequenced cotton genomes indicated that homologous genes, anchored by NAU2395 and NAU1092, were aligned to the G. arboreum genome within a physical distance between 81.10 Mbps and 87.07 Mbps. In that interval, several genes were confirmed in literatures to associate with fiber development. Among these genes, seven genes were further selected for an expression analysis through fiber development transcriptome database, revealing unique expression patterns across different stages of fiber development between CCRI 8 and Pima 90-53. The genes underlying qFS-c9-1 were favorable to fine mapping and cloning. The current study results provided valuable evidence for mapping stable QTL of fiber strength utilizing multiple populations and environments, as well as map-based cloning the candidate gene underlying the QTL for future prospective research directions.

Structure, evolution, and comparative genomics of tetraploid cotton based on a high-density genetic linkage map

A high-density linkage map was constructed using 1,885 newly obtained loci and 3,747 previously published loci, which included 5,152 loci with 4696.03 cM in total length and 0.91 cM in mean distance. Homology analysis in the cotton genome further confirmed the 13 expected homologous chromosome pairs and revealed an obvious inversion on Chr10 or Chr20 and repeated inversions on Chr07 or Chr16. In addition, two reciprocal translocations between Chr02 and Chr03 and between Chr04 and Chr05 were confirmed. Comparative genomics between the tetraploid cotton and the diploid cottons showed that no major structural changes exist between D_T and D chromosomes but rather between A_T and A chromosomes. Blast analysis between the tetraploid cotton genome and the mixed genome of two diploid cottons showed that most AD chromosomes, regardless of whether it is from the A_T or D_T genome, preferentially matched with the corresponding homologous chromosome in the diploid A genome, and then the corresponding homologous chromosome in the diploid D genome, indicating that the diploid D genome underwent converted evolution by the diploid A genome to form the D_T genome during polyploidization. In addition, the results reflected that a series of chromosomal translocations occurred among Chr01/Chr15, Chr02/Chr14, Chr03/Chr17, Chr04/Chr22, and Chr05/Chr19.

Genome Wide SSR High Density Genetic Map Construction from an Interspecific Cross of Gossypium hirsutum × Gossypium tomentosum

A high density genetic map was constructed using F2 population derived from an interspecific cross of G. hirsutum × G. tomentosum. The map consisted of 3093 marker loci distributed across all the 26 chromosomes and covered 4365.3 cM of cotton genome with an average inter-marker distance of 1.48 cM. The maximum length of chromosome was 218.38 cM and the minimum was 122.09 cM with an average length of 167.90 cM. A sub-genome covers more genetic distance (2189.01 cM) with an average inter loci distance of 1.53 cM than D sub-genome which covers a length of 2176.29 cM with an average distance of 1.43 cM. There were 716 distorted loci in the map accounting for 23.14% and most distorted loci were distributed on D sub-genome (25.06%), which were more than on A sub-genome (21.23%). In our map 49 segregation hotspots (SDR) were distributed across the genome with more on D sub-genome as compared to A genome. Two post-polyploidization reciprocal translocations of "A2/A3 and A4/A5" were suggested by seven pairs of duplicate loci. The map constructed through these studies is one of the three densest genetic maps in cotton however; this is the first dense genome wide SSR interspecific genetic map between G. hirsutum and G. tomentosum.

Genomic heterozygosity and hybrid breakdown in cotton (Gossypium): different traits, different effects

Background

Hybrid breakdown has been well documented in various species. Relationships between genomic heterozygosity and traits-fitness have been extensively explored especially in the natural populations. But correlations between genomic heterozygosity and vegetative and reproductive traits in cotton interspecific populations have not been studied. In the current study, two reciprocal F₂ populations were developed using Gossypium hirsutum cv. Emian 22 and G. barbadense acc. 3–79 as parents to study hybrid breakdown in cotton. A total of 125 simple sequence repeat (SSR) markers were used to genotype the two F₂ interspecific populations.

Results

To guarantee mutual independence among the genotyped markers, the 125 SSR markers were checked by the linkage disequilibrium analysis. To our knowledge, this is a novel approach to evaluate the individual genomic heterozygosity. After marker checking, 83 common loci were used to assess the extent of genomic heterozygosity. Hybrid breakdown was found extensively in the two interspecific F₂ populations particularly on the reproductive traits because of the infertility and the bare seeds. And then, the relationships between the genomic heterozygosity and the vegetative reproductive traits were investigated. The only relationships between hybrid breakdown and heterozygosity were observed in the (Emian22 × 3–79) F₂ population for seed index (SI) and boll number per plant (BN). The maternal cytoplasmic environment may have a significant effect on genomic heterozygosity and on correlations between heterozygosity and reproductive traits.

Conclusions

A novel approach was used to evaluate genomic heterozygosity in cotton; and hybrid breakdown was observed in reproductive traits in cotton. These findings may offer new insight into hybrid breakdown in allotetraploid cotton interspecific hybrids, and may be useful for the development of interspecific hybrids for cotton genetic improvement.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0366-5) contains supplementary material, which is available to authorized users.

A Gossypium BAC clone contains key repeat components distinguishing sub-genome of allotetraploidy cottons

Background

Dissecting genome organization is indispensable for further functional and applied studies. As genome sequences data shown, cotton genomes contain more than 60 % repetitive sequences, so study on repetitive sequences composition, structure, and distribution is the key step to dissect cotton genome.

Results

In this study, a bacterial artificial chromosome (BAC) clone enriched in repetitive sequences, was discovered initiatively by fluorescence in situ hybridization (FISH). FISHing with allotetraploidy cotton as target DNA, dispersed signals on most regions of all A sub-genome chromosomes, and only middle regions of all D sub-genome chromosomes were detected. Further FISHing with other cotton species bearing A or D genome as target DNA, specific signals were viewed. After BAC sequencing and bioinformational analysis, 129 repeat elements, size about 57,172 bp were found, accounting for more than 62 % of the BAC sequence (91,238 bp). Among them, a type of long terminal repeat-retrotransposon (LTR-RT), LTR/Gypsy was the key element causing the specific FISH results. Using the fragments of BAC matching with the identified Gypsy-like LTR as probes, the BAC-57I23-like FISH signals were reappeared. Running BLASTN, the fragments had good match with all chromosomes of G. arboreum (A₂) genome and A sub-genome of G. hirsutum (AD₁), and had relatively inferior match with all chromosomes of D sub-genome of AD₁, but had little match with the chromosomes of G. raimondii (D₅) genome, which was consistent with the FISH results.

Conclusion

A repeats-enriched cytogenetic marker to identify A and D sub-genomes of Gossypium was discovered by FISH. Combined sequences analysis with FISH verification, the assembly quality of repetitive sequences in the allotetraploidy cotton draft genome was assessed, and better chromosome belonging was verified. We also found the genomic distribution of the identified Gypsy-LTR-RT was similar to the distribution of heterochromatin. The expansion of this type of Gypsy-LTR-RT in heterochromatic regions may be one of the major reasons for the size gap between A and D genome. The findings showed here will help to understand the composition, structure, and evolution of cotton genome, and contribute to the further perfection of the draft genomes of cotton.

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.

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.

Comparative genomic de-convolution of the cotton genome revealed a decaploid ancestor and widespread chromosomal fractionation

The 'apparently' simple genomes of many angiosperms mask complex evolutionary histories. The reference genome sequence for cotton (Gossypium spp.) revealed a ploidy change of a complexity unprecedented to date, indeed that could not be distinguished as to its exact dosage. Herein, by developing several comparative, computational and statistical approaches, we revealed a 5× multiplication in the cotton lineage of an ancestral genome common to cotton and cacao, and proposed evolutionary models to show how such a decaploid ancestor formed. The c. 70% gene loss necessary to bring the ancestral decaploid to its current gene count appears to fit an approximate geometrical model; that is, although many genes may be lost by single-gene deletion events, some may be lost in groups of consecutive genes. Gene loss following cotton decaploidy has largely just reduced gene copy numbers of some homologous groups. We designed a novel approach to deconvolute layers of chromosome homology, providing definitive information on gene orthology and paralogy across broad evolutionary distances, both of fundamental value and serving as an important platform to support further studies in and beyond cotton and genomics communities.

Identification and Characterization of Segregation Distortion Loci on Cotton Chromosome 18

Segregation distortion is commonly detected via genetic mapping and this phenomenon has been reported in many species. However, the genetic causes of the segregation distortion regions in a majority of species are still unclear. To genetically dissect the SD on chromosome 18 in cotton, eight reciprocal backcross populations and two F₂ populations were developed. Eleven segregation distortion loci (SDL) were detected in these ten populations. Comparative analyses among populations revealed that SDL18.1 and SDL18.9 were consistent with male gametic competition; whereas SDL18.4 and SDL18.11 reflected female gametic selection. Similarly, other SDL could reflect zygotic selection. The surprising finding was that SDL18.8 was detected in all populations, and the direction was skewed towards heterozygotes. Consequently, zygotic selection or heterosis could represent the underlying genetic mechanism for SDL18.8. Among developed introgression lines, SDL18.8 was introgressed as a heterozygote, further substantiating that a heterozygote state was preferred under competition. Six out of 11 SDL on chromosome 18 were dependent on the cytoplasmic environment. These results indicated that different SDL showed varying responses to the cytoplasmic environment. Overall, the results provided a novel strategy to analyze the molecular mechanisms, which could be further exploited in cotton interspecific breeding programs.

Partial Dominance, Overdominance and Epistasis as the Genetic Basis of Heterosis in Upland Cotton (Gossypium hirsutum L.)

Determination of genetic basis of heterosis may promote hybrid production in Upland cotton (Gossypium hirsutum L.). This study was designed to explore the genetic mechanism of heterosis for yield and yield components in F2: 3 and F2: 4 populations derived from a hybrid 'Xinza No. 1'. Replicated yield field trials of the progenies were conducted in 2008 and 2009. Phenotypic data analyses indicated overdominance in F1 for yield and yield components. Additive and dominance effects at single-locus level and digenic epistatic interactions at two-locus level were analyzed by 421 marker loci spanning 3814 cM of the genome. A total of 38 and 49 QTLs controlling yield and yield components were identified in F2: 3 and F2: 4 populations, respectively. Analyses of these QTLs indicated that the effects of partial dominance and overdominance contributed to heterosis in Upland cotton simultaneously. Most of the QTLs showed partial dominance whereas 13 QTLs showing overdominance in F2:3 population, and 19 QTLs showed overdominance in F2:4. Among them, 21 QTLs were common in both F2: 3 and F2: 4 populations. A large number of two-locus interactions for yield and yield components were detected in both generations. AA (additive × additive) epistasis accounted for majority portion of epistatic effects. Thirty three complementary two-locus homozygotes (11/22 and 22/11) were the best genotypes for AA interactions in terms of bolls per plant. Genotypes of double homozygotes, 11/22, 22/11 and 22/22, performed best for AD/DA interactions, while genotype of 11/12 performed best for DD interactions. These results indicated that (1) partial dominance and overdominance effects at single-locus level and (2) epistasis at two-locus level elucidated the genetic basis of heterosis in Upland cotton.

A BIL Population Derived from G. hirsutum and G. barbadense Provides a Resource for Cotton Genetics and Breeding

To provide a resource for cotton genetics and breeding, an interspecific hybridization between Gossypium hirsutum cv. Emian22 and G. barbadense acc. 3–79 was made. A population of 54 BILs (backcross inbred lines, BC₁F₈) was developed with the aim of transferring G. barbadense genes into G. hirsutum in order to genetically analyze these genes’ function in a G. hirsutum background and create new germplasms for breeding. Preliminary investigation of the morphological traits showed that the BILs had diverse variations in plant architecture, seed size, and fuzz color; the related traits of yield and fiber quality evaluated in 4 environments also showed abundant phenotypic variation. In order to explore the molecular diversity of the BIL population, 446 SSR markers selected at an average genetic distance of 10 cM from our interspecific linkage map were used to genotype the BIL population. A total of 393 polymorphic loci accounting for 84.4% MAF (major allele frequency) > 0.05 and 922 allele loci were detected, and the Shannon diversity index (I) was 0.417 per locus. The average introgression segment length was 16.24 cM, and an average of 29.53 segments were introgressed in each BIL line with an average background recovery of 79.8%. QTL mapping revealed 58 QTL associated with fiber quality and yield traits, and 47 favored alleles derived from the donor parent were discovered. This study demonstrated that the interspecific BIL population was enriched with much phenotypic and molecular variation which could be a resource for cotton genetics and breeding.

Construction of an EST-SSR-based interspecific transcriptome linkage map of fibre development in cotton

Quantitative trait locus (QTL) mapping is an important method in marker-assisted selection breeding. Many studies on the QTLs focus on cotton fibre yield and quality; however, most are conducted at the DNA level, which may reveal null QTLs. Hence, QTL mapping based on transcriptome maps at the cDNA level is often more reliable. In this study, an interspecific transcriptome map of allotetraploid cotton was developed based on an F2 population (Emian22 x 3-79) by amplifying cDNA using EST-SSRs. The map was constructed using cDNA obtained from developing fibres at five days post anthesis (DPA). A total of 1270 EST-SSRs were screened for polymorphisms between the mapping parents. The resulting transcriptome linkage map contained 242 markers that were distributed in 32 linkage groups (26 chromosomes). The full length of this map is 1938.72 cM with a mean marker distance of 8.01 cM. The functions of some ESTs have been annotated by exploring homologous sequences. Some markers were related to the differentiation and elongation of cotton fibre, while most were related to the basic metabolism. This study demonstrates that constructing a transcriptome linkage map by amplifying cDNAs using EST-SSRs is a simple and practical method as well as a powerful tool to map eQTLs for fibre quality and other traits in cotton.

A high-density SSR genetic map constructed from a F2 population of Gossypium hirsutum and Gossypium darwinii

The cultivated allotetraploid species Gossypium hirsutum, accounts for 90% of the world cotton production, has narrow genetic basis that's why its yield, quality or stress resistance breeding is stagnant. It is therefore, essential to explore desirable genes from Gossypium darwinii which has enviable traits such as high fiber fineness, drought tolerance, fusarium and verticillium resistance. We used G. darwinii as primary plant materials in this study not only to enrich the genetic diversity of exiting germplasm but also to better understand its genome structure. An interspecific high density linkage map of allotetraploid cotton was constructed using F2 population (G. hirsutum×G. darwinii). The map was based entirely on genome-wide simple sequence repeat (SSR) markers. A total of 2763 markers were mapped in 26 linkage groups (chromosomes) covering a genome length of 4176.7cM with an average inter-locus distance of 1.5cM. The length of the chromosomes ranged from 84.7 to 238.5cM with an average length of 160.6cM. At subgenome length was 2160.7cM with an average distance of 1.6cM, where as Dt genome length was 2016cM with an average distance of 1.4cM. There were 601 distorted SSR loci. Less number of segregation distortion loci were located in At subgenome than in Dt subgenome. Two post-polyploidization reciprocal translocations of "A2/A3 and A4/A5" were suggested by 44 pairs of duplicate loci.