Search

During the domestication and improvement processes, interspecific introgression from Gossypium hirsutum has reorganized the genomic architecture of Gossypium barbadense; however, the introgression details and the trait-related genetic loci remain largely unknown. Here, we perform a genome-wide population analysis and genetically recategorize 365 G. barbadense accessions into four subgroups which is different from previous categorizations. A total of 315 introgression events from G. hirsutum to G. barbadense, which primarily contributed to population divergence and agronomic trait variation in G. barbadense, are identified. We find that 70% introgression from G. hirsutum have greatly increased the genetic diversity and divergence of G. barbadense. Some loci are identified with divergent haplotype selection for adaptation to the environment at high latitudes. Through genome-wide association study and genome linkage disequilibrium interval haplotyping analyses, two fiber micronaire-related haplotype blocks are detected, one of which (FM2) is introgressed from G. hirsutum. Seven distinguished traits related to growth period, plant architecture, and stronger vegetative growth habit are found to have pleiotropic effects controlled by a single gene in G. barbadense. Thus, this study provides new insights into the breeding history of G. barbadense and highlights introgression is a driver for improving cultivars in G. barbadense.

Background Colored cotton is also called eco-cotton because of its natural color fiber. However, the yield and quality of natural colored cotton are generally lower than white cotton. Now, little is known about the regulatory mechanisms involved in biochemical and metabolic pathways and critical genes between fiber quality and pigment synthesis. The present study used a pair of brown cotton near-isogenic lines (brown cotton ‘Z1-61’ and white cotton ‘RT’) to determine the transcriptomic and proteomic data of fiber (ovule) at 0DPA, 10DPA, and 20DPA. Results Integration analysis of the two-omics data demonstrated that the differential proteins with the same expression trend in the adjacent developmental stages were significantly and positively correlated with fold change. Enrichment analysis of metabolic pathway revealed that Z1-61, enriched in fiber elongation genes related to a flavonoid biosynthesis, phenylalanine metabolism, glutathione metabolism and many more genes (proteins) are up-regulated. Short Time-series Expression Miner analysis inferred that up-regulated genes of Z1-61 were specifically involved in the metabolism of carbohydrates, energy, and lipid, suggesting that these pathways play an important role in fiber elongation and pigmentation. We further analyzed the 164 glycosyltransferases genes were predominantly expressed in the early stage of Z1-61 fiber elongation while they were predominantly expressed during the later period of RT fibers elongation. And 15 MYB-bHLH-WD40 triple complex genes and other transcription factor genes such as C2H2 (12), ERF (11), and NAC (7) were also preferentially expressed in Z1-61. Weighted correlation network analysis exhibited that fatty acid synthesis and energy metabolism were the active metabolic pathways in both cotton genotypes during fiber development. The 10 hub genes obtained from the final screening can be considered to play a critical role in the process of fiber pigment deposition while ensuring that fiber quality does not degraded. Conclusion This is the first complementary analyses of the transcriptome and proteome related with the quality of brown cotton fiber. Those studies can be understood as better playing a major role in fiber quality enhancement and pigmentation during fiber elongation.

Background Heterosis has been extensively utilized in different crops and made a significant contribution to global food security. Genetic distance (GD) is one of the valuable criteria for selecting parents in hybrid breeding. The objectives of this study were to estimate the GD between parents using both simple sequence repeat (SSR) markers and single nucleotide polymorphism (SNP) markers and to investigate the efficiency of the prediction of hybrid performance based on GD. The experiment comprised of four male parents, 282 female parents and 1128 F1, derived from NCII mating scheme. The hybrids, their parents and two check cultivars were evaluated for two years. Performance of F1, mid-parent heterosis (MPH), and best parent heterosis (BPH) were evaluated for ten agronomic and fiber quality traits, including plant height, boll weight, boll number, lint percentage, fiber length, fiber strength, fiber uniformity, fiber elongation ratio, micronaire, and spinning consistent index. Results Heterosis was observed in all hybrids and, the traits like plant height, boll number, boll weight and lint percentage exhibited higher heterosis than the fiber quality traits. Correlations were significant between parental and F1 performances. The F1 performances between three hybrid sets (Elite×Elite, Exotic×Elite, and Historic×Elite) showed significant differences in eight traits, including boll number, lint percentage, fiber length, fiber strength, fiber uniformity, fiber elongation ratio, micronaire, and spinning consistent index. The correlation of the GD assessed by both SSR and SNP markers was significantly positive. The cluster analysis based on GD results estimated using SNP showed that all the female parents divided into five groups and the F1 performance between these five groups showed significant differences in four traits, including lint percentage, micronaire, fiber strength, and fiber elongation ratio. The correlation between GD and F1 performance, MPH and BPH were significant for lint percentage and micronaire. Conclusions Our results suggested that GD between parents could be helpful in heterosis prediction for certain traits. This study reveals that molecular marker analysis can serve as a basis for assigning germplasm into heterotic groups and to provide guidelines for parental selection in hybrid cotton breeding.

SUMMARY An evaluation of combining ability can facilitate the selection of suitable parents and superior F1 hybrids for hybrid cotton breeding, although the molecular genetic basis of combining ability has not been fully characterized. In the present study, 282 female parents were crossed with four male parents in accordance with the North Carolina II mating scheme to generate 1128 hybrids. The parental lines were genotyped based on restriction site‐associated DNA sequencing and 306 814 filtered single nucleotide polymorphisms were used for genome‐wide association analysis involving the phenotypes, general combining ability (GCA) values, and specific combining ability values of eight fiber quality‐ and yield‐related traits. The main results were: (i) all parents could be clustered into five subgroups based on population structure analyses and the GCA performance of the female parents had significant differences between subgroups; (ii) 20 accessions with a top 5% GCA value for more than one trait were identified as elite parents for hybrid cotton breeding; (iii) 120 significant single nucleotide polymorphisms, clustered into 66 quantitative trait loci, such as the previously reported Gh_A07G1769 and GhHOX3 genes, were found to be significantly associated with GCA; and (iv) identified quantitative trait loci for GCA had a cumulative effect on GCA of the accessions. Overall, our results suggest that pyramiding the favorable loci for GCA may improve the efficiency of hybrid cotton breeding., Significance Statement General and specific combining ability values are very important factors for the selection of appropriate parents and superior combinations in hybrid cotton breeding, although the molecular genetic basis has not been fully characterized. In the present study, we identified significant single nucleotide polymorphisms for general and specific combining ability in a large‐scale North Carolina II population through a genome‐wide association study. The results obtained provide new insights with respect to our understanding of the genetic factors related to combining ability.

DNA methylation has a growing potential for use as a biomarker because of its involvement in disease. DNA methylation data have also substantially grown in volume during the past 5 years. To facilitate access to these fragmented data, we proposed DiseaseMeth version 3.0 based on DiseaseMeth version 2.0, in which the number of diseases including increased from 88 to 162 and High-throughput profiles samples increased from 32 701 to 49 949. Experimentally confirmed associations added 448 pairs obtained by manual literature mining from 1472 papers in PubMed. The search, analyze and tools sections were updated to increase performance. In particular, the FunctionSearch now provides for the functional enrichment of genes from localized GO and KEGG annotation. We have also developed a unified analysis pipeline for identifying differentially DNA methylated genes (DMGs) from the original data stored in the database. 22 718 DMGs were found in 99 diseases. These DMGs offer application in disease evaluation using two self-developed online tools, Methylation Disease Correlation and Cancer Prognosis & Co-Methylation. All query results can be downloaded and can also be displayed through a box plot, heatmap or network module according to whichever search section is used. DiseaseMeth version 3.0 is freely available at http://diseasemeth.edbc.org/.

Upland cotton (Gossypium hirsutum L.) is an important economic crop for renewable textile fibers. However, the simultaneous improvement of yield and fiber quality in cotton is difficult as the linkage drag. Compared with breaking the linkage drag, identification of the favorable pleiotropic loci on the genome level by genome-wide association study (GWAS) provides a new way to improve the yield and fiber quality simultaneously. In our study restriction-site-associated DNA sequencing (RAD-seq) was used to genotype 316 cotton accessions. Eight major traits in three categories including yield, fiber quality and maturation were investigated in nine environments (3 sites × 3 years). 231 SNPs associated with these eight traits (− log10(P) > 5.27) were identified, located in 27 genomic regions respectively by linkage disequilibrium analysis. Further analysis showed that four genomic regions (the region 1, 6, 8 and 23) held favorable pleiotropic loci and 6 candidate genes were identified. Through genotyping, 14 elite accessions carrying the favorable loci on four pleiotropic regions were identified. These favorable pleiotropic loci and elite genotypes identified in this study will be utilized to improve the yield and fiber quality simultaneously in future cotton breeding.

For about a century, plant breeding has widely exploited the heterosis phenomenon–often considered as hybrid vigor–to increase agricultural productivity. The ensuing F1 hybrids can substantially outperform their progenitors due to heterozygous combinations that mitigate deleterious mutations occurring in each genome. However, only fragmented knowledge is available concerning the underlying genes and processes that foster heterosis. Although cotton is among the highly valued crops, its improvement programs that involve the exploitation of heterosis are still limited in terms of significant accomplishments to make it broadly applicable in different agro-ecological zones. Here, F1 hybrids were derived from mating a diverse Upland Cotton germplasm with commercially valuable cultivars in the Line × Tester fashion and evaluated across multiple environments for 10 measurable traits. These traits were dissected into five different heterosis types and specific combining ability (SCA). Subsequent genome-wide predictions along-with association analyses uncovered a set of 298 highly significant key single nucleotide polymorphisms (SNPs)/Quantitative Trait Nucleotides (QTNs) and 271 heterotic Quantitative Trait Nucleotides (hQTNs) related to agronomic and fiber quality traits. The integration of a genome wide association study with RNA-sequence analysis yielded 275 candidate genes in the vicinity of key SNPs/QTNs. Fiber micronaire (MIC) and lint percentage (LP) had the maximum number of associated genes, i.e., each with 45 related to QTNs/hQTNs. A total of 54 putative candidate genes were identified in association with HETEROSIS of quoted traits. The novel players in the heterosis mechanism highlighted in this study may prove to be scientifically and biologically important for cotton biologists, and for those breeders engaged in cotton fiber and yield improvement programs.

Large-scale genomic surveys of crop germplasm are important for understanding the genetic architecture of favorable traits. The genomic basis of geographic differentiation and fiber improvement in cultivated cotton is poorly understood. Here, we analyzed 3,248 tetraploid cotton genomes and confirmed that the extensive chromosome inversions on chromosomes A06 and A08 underlies the geographic differentiation in cultivated Gossypium hirsutum. We further revealed that the haplotypic diversity originated from landraces, which might be essential for understanding adaptative evolution in cultivated cotton. Introgression and association analyses identified new fiber quality-related loci and demonstrated that the introgressed alleles from two diploid cottons had a large effect on fiber quality improvement. These loci provided the potential power to overcome the bottleneck in fiber quality improvement. Our study uncovered several critical genomic signatures generated by historical breeding effects in cotton and a wealth of data that enrich genomic resources for the research community.

Improving yield is a primary mission for cotton (Gossypium hirsutum) breeders; development of cultivars with suitable architecture for high planting density (HPDA) can increase yield per unit area. We characterized a natural cotton mutant, AiSheng98 (AS98), which exhibits shorter height, shorter branch length, and more acute branch angle than wild-type. A copy number variant at the HPDA locus on Chromosome D12 (HPDA-D12), encoding a dehydration-responsive element-binding (DREB) transcription factor, GhDREB1B, strongly affects plant architecture in the AS98 mutant. We found an association between a tandem duplication of a c. 13.5 kb segment in HPDA-D12 and elevated GhDREB1B expression resulting in the AS98 mutant phenotype. GhDREB1B overexpression confers a significant decrease in plant height and branch length, and reduced branch angle. Our results suggest that fine-tuning GhDREB1B expression may be a viable engineering strategy for modification of plant architecture favorable to high planting density in cotton.

Upland cotton is the most important natural-fiber crop. The genomic variation of diverse germplasms and alleles underpinning fiber quality and yield should be extensively explored. Here, we resequenced a core collection comprising 419 accessions with 6.55-fold coverage depth and identified approximately 3.66 million SNPs for evaluating the genomic variation. We performed phenotyping across 12 environments and conducted genome-wide association study of 13 fiber-related traits. 7,383 unique SNPs were significantly associated with these traits and were located within or near 4,820 genes; more associated loci were detected for fiber quality than fiber yield, and more fiber genes were detected in the D than the A subgenome. Several previously undescribed causal genes for days to flowering, fiber length, and fiber strength were identified. Phenotypic selection for these traits increased the frequency of elite alleles during domestication and breeding. These results provide targets for molecular selection and genetic manipulation in cotton improvement.

Fiber yield and quality are the most important traits for Upland cotton (Gossypium hirsutum L.). Identifying high yield and good fiber quality genes are the prime concern of researchers in cotton breeding. Association mapping offers an alternative and powerful method for detecting those complex agronomic traits. In this study, 198 simple sequence repeats (SSRs) were used to screen markers associated with fiber yield and quality traits with 302 elite Upland cotton accessions that were evaluated in 12 locations representing the Yellow River and Yangtze River cotton growing regions of China. Three subpopulations were found after the estimation of population structure. The pair-wise kinship values varied from 0 to 0.867. Only 1.59% of the total marker locus pairs showed significant linkage disequilibrium (LD, p

Background Cottonseed is one of the most important raw materials for plant protein, oil and alternative biofuel for diesel engines. Understanding the complex genetic basis of cottonseed traits is requisite for achieving efficient genetic improvement of the traits. However, it is not yet clear about their genetic architecture in genomic level. GWAS has been an effective way to explore genetic basis of quantitative traits in human and many crops. This study aims to dissect genetic mechanism seven cottonseed traits by a GWAS for genetic improvement. Results A genome-wide association study (GWAS) based on a full gene model with gene effects as fixed and gene-environment interaction as random, was conducted for protein, oil and 5 fatty acids using 316 accessions and ~ 390 K SNPs. Totally, 124 significant quantitative trait SNPs (QTSs), consisting of 16, 21, 87 for protein, oil and fatty acids (palmitic, linoleic, oleic, myristic, stearic), respectively, were identified and the broad-sense heritability was estimated from 71.62 to 93.43%; no QTS-environment interaction was detected for the protein, the palmitic and the oleic contents; the protein content was predominantly controlled by epistatic effects accounting for 65.18% of the total variation, but the oil content and the fatty acids except the palmitic were mainly determined by gene main effects and no epistasis was detected for the myristic and the stearic. Prediction of superior pure line and hybrid revealed the potential of the QTSs in the improvement of cottonseed traits, and the hybrid could achieve higher or lower genetic values compared with pure lines. Conclusions This study revealed complex genetic architecture of seven cottonseed traits at whole genome-wide by mixed linear model approach; the identified genetic variants and estimated genetic component effects of gene, gene-gene and gene-environment interaction provide cotton geneticist or breeders new knowledge on the genetic mechanism of the traits and the potential molecular breeding design strategy. Electronic supplementary material The online version of this article (10.1186/s12864-018-4837-0) contains supplementary material, which is available to authorized users.

The ancestors of Gossypium arboreum and Gossypium herbaceum provided the A subgenome for the modern cultivated allotetraploid cotton. Here, we upgraded the G. arboreum genome assembly by integrating different technologies. We resequenced 243 G. arboreum and G. herbaceum accessions to generate a map of genome variations and found that they are equally diverged from Gossypium raimondii. Independent analysis suggested that Chinese G. arboreum originated in South China and was subsequently introduced to the Yangtze and Yellow River regions. Most accessions with domestication-related traits experienced geographic isolation. Genome-wide association study (GWAS) identified 98 significant peak associations for 11 agronomically important traits in G. arboreum. A nonsynonymous substitution (cysteine-to-arginine substitution) of GaKASIII seems to confer substantial fatty acid composition (C16:0 and C16:1) changes in cotton seeds. Resistance to fusarium wilt disease is associated with activation of GaGSTF9 expression. Our work represents a major step toward understanding the evolution of the A genome of cotton.

A second regulatory gene (jadR(1)) is located immediately upstream of the putative repressor gene (jadR(2)) in the jad cluster for biosynthesis of the antibiotic jadomycin B in Streptomyces venezuelae ISP5230. It encodes a 234-amino acid polypeptide with a sequence resembling those of response regulator proteins in two-component control systems. Features in the conserved C-terminal domain of JadR(1) place the protein in the OmpR-PhoB subfamily of response regulators. In mutants where jadR(1) was deleted or disrupted, jadomycin B was not produced, implying that the gene has an essential role in biosynthesis of the antibiotic. Cloning jadR(1) from S. venezuelae in pJV73A, and introducing additional copies of the gene into the wild-type parent by plasmid transformation gave unstable strains with pJV73A integrated into the chromosome. The transformants initially showed increased production of jadomycin B but gave lower titers as excess copies of jadR(1) were lost; mature cultures stabilized with a wild-type level of antibiotic production. The mutant from which jadR(1) had been deleted could not be transformed with pJV73A. Altering the composition of jadR genes in the chromosome by integration of vectors carrying intact and disrupted copies of jadR(1) and jadR(2) provided evidence that the two genes form a regulatory pair different in function from previously reported two-component systems controlling antibiotic biosynthesis in streptomycetes.

Three new genes (jadW 1, jadW 2 and jadW 3) were isolated from a region of the Streptomyces venezuelae ISP5230 chromosome at the left-hand end of the jad cluster for jadomycin B (JdB) biosynthesis. The deduced amino acid sequence of jadW 1 showed strong similarity to gene products associated in several streptomycetes with γ-butyrolactone autoregulators controlling morphological differentiation and secondary metabolism. Examination of JadW1 for conserved domains detected a repeat sequence characteristic of proteins in the AfsA regulatory family. Insertional inactivation of jadW 1 reduced the growth rate of S. venezuelae cultures in aerated liquid media containing complex nitrogen sources and altered growth morphology in minimal medium. It also affected sporulation on agar media. Cultures of jadW 1-disrupted mutants grown under conditions supporting biosynthesis of JdB or chloramphenicol by the wild-type strain failed to produce either of the antibiotics. Complementing the disrupted strain by transformation with pJV435, containing a cloned copy of the gene, improved sporulation and restored antibiotic biosynthesis in transformants to titres close to those of the wild-type similarly transformed with pJV435 as a control. The results are consistent with a role for jadW 1 in regulating morphological and metabolic differentiation. Further sequence analysis of jadR 2, which functions with jadR 1 in stress-induced activation of JdB biosynthesis, indicated that this gene encodes a γ-butyrolactone receptor homologue. The growth-rate-sensitive phenotype of the jadW 1-disrupted mutant, and the proximity of jadW 1 to jadR 2 indicate that this region of the jad gene cluster contains a regulatory mechanism incorporating γ-butyrolactone signalling and sensitivity to environmental stress.

Improving yield is a major objective for cotton breeding schemes, and lint yield and its three component traits (boll number, boll weight and lint percentage) are complex traits controlled by multiple genes and various environments. Association mapping was performed to detect markers associated with these four traits using 651 simple sequence repeats (SSRs). A mixed linear model including epistasis and environmental interaction was used to screen the loci associated with these four yield traits by 323 accessions of Gossypium hirsutum L. evaluated in nine different environments. 251 significant loci were detected to be associated with lint yield and its three components, including 69 loci with individual effects and all involved in epistasis interactions. These significant loci explain ∼ 62.05% of the phenotypic variance (ranging from 49.06% ∼ 72.29% for these four traits). It was indicated by high contribution of environmental interaction to the phenotypic variance for lint yield and boll numbers, that genetic effects of SSR loci were susceptible to environment factors. Shared loci were also observed among these four traits, which may be used for simultaneous improvement in cotton breeding for yield traits. Furthermore, consistent and elite loci were screened with −Log10 (P-value) >8.0 based on predicted effects of loci detected in different environments. There was one locus and 6 pairs of epistasis for lint yield, 4 loci and 10 epistasis for boll number, 15 loci and 2 epistasis for boll weight, and 2 loci and 5 epistasis for lint percentage, respectively. These results provided insights into the genetic basis of lint yield and its components and may be useful for marker-assisted breeding to improve cotton production.