A genome-wide association study of early-maturation traits in upland cotton based on the CottonSNP80K array

Seed treatment: an alternative and sustainable approach to cotton seed delinting

This review article delves into the vital aspects of cotton, emphasizing its global significance as a crucial agricultural commodity. The paper comprehensively explores the composition of cotton and surveys the diverse methods employed for the removal of cotton lint from seeds. Conventional delinting methods, including mechanical and chemical approaches, are scrutinized in terms of their advantages and drawbacks. However, the primary focus of this review is on highlighting the emerging significance of biological delinting methods. By harnessing the power of microbial enzymes and organisms, biological approaches offer a promising alternative for efficient lint removal. The authors discuss the environmental advantages associated with biological delinting, positioning it as a sustainable solution that mitigates the ecological impact of traditional methods. Furthermore, the article contextualizes these delinting methods within the framework of Sustainable Development Goals (SDGs) and underscores the importance of adopting eco-friendly practices in the cotton industry to align with SDG goals. By accentuating the potential of biological delinting in contributing to sustainable agriculture and responsible production, the review advocates for a paradigm shift towards more environmentally conscious approaches in the cotton sector. Overall, the article aims to provide a comprehensive perspective on cotton delinting methods, emphasizing the pivotal role of biological alternatives in fostering a sustainable and goal-oriented future for the cotton industry.

Integrating RTM-GWAS and meta‑QTL data revealed genomic regions and candidate genes associated with the first fruit branch node and its height in upland cotton

KEY MESSAGE

Two genomic regions associated with FFBN and HFFBN and a potential regulatory gene (GhE6) of HFFBN were identified through the integration of RTM-GWAS and meta‑QTL analyses. Abstract The first fruit branch node (FFBN) and the height of the first fruit branch node (HFFBN) are two important traits that are related to plant architecture and early maturation in upland cotton. Several studies have been conducted to elucidate the genetic basis of these traits in cotton using biparental and natural populations. In this study, by using 9,244 SNP linkage disequilibrium block (SNPLDB) loci from 315 upland cotton accessions, we carried out restricted two-stage multilocus and multiallele genome-wide association studies (RTM-GWASs) and identified promising haplotypes/alleles of the four stable and true major SNPLDB loci that were significantly associated with FFBN and HFFBN. Additionally, a meta-quantitative trait locus (MQTL) analysis was conducted on 274 original QTLs that were reported in 27 studies, and 40 MQTLs associated with FFBN and HFFBN were identified. Through the integration of the RTM-GWAS and meta‑QTL analyses, two stable and true major SNPLDBs (LDB_5_15144433 and LDB_16_37952328) that were distributed in the two MQTLs were identified. Ultimately, 142 genes in the two genomic regions were annotated, and three candidate genes associated with FFBN and HFFBN were identified in the genomic region (A05:14.64-15.64 Mb) via RNA-Seq and qRT‒PCR. The results of virus-induced gene silencing (VIGS) experiments indicated that GhE6 was a key gene related to HFFBN and that GhDRM1 and GhGES were important genes associated with early flowering in upland cotton. These findings will aid in the future identification of molecular markers and genetic resources for developing elite early-maturing cultivars with ideal plant characteristics.

Genome-wide association analysis identified molecular markers and candidate genes for flower traits in Chinese orchid (Cymbidium sinense)

The orchid, the champagne of flowers, brings luxury, elegance, and novelty to nature. Cymbidium sinense is a symbol of gigantic floral variability on account of wavering shapes and sizes of floral organs, although marker-trait association (MTA) has not been studied for its floral traits. We evaluated markers associated with 14 floral traits of C. sinense through a genome-wide association study (GWAS) of 195 accessions. A total of 65 318 522 single-nucleotide polymorphisms (SNPs) and 3 906 176 insertion/deletion (InDel) events were identified through genotyping-by-sequencing. Among these, 4694 potential SNPs and 477 InDels were identified as MTAs at -log10 P > 5. The genes related to these SNPs and InDels were largely associated with floral regulators, hormonal pathways, cell division, and metabolism, playing essential roles in tailoring floral morphology. Moreover, 20 candidate SNPs/InDels linked to 11 genes were verified, 8 of which were situated on exons, one was located in the 5'-UTR and two were positioned in introns. Here, the multitepal trait-related gene RABBIT EARS (RBE) was found to be the most crucial gene. We analyzed the role of CsRBE in the regulation of flower-related genes via efficient transient overexpression in C. sinense protoplasts, and found that the floral homeotic genes CsAP3 and CsPI, as well as organ boundary regulators, including CsCUC and CsTCP genes, were regulated by CsRBE. Thus, we obtained key gene loci for important ornamental traits of orchids using genome-wide association analysis of populations with natural variation. The findings of this study can do a great deal to expedite orchid breeding programs for shape variability.

Association Analysis of Tiller-Related Traits with EST-SSR Markers in Psathyrostachys juncea

Psathyrostachys juncea is a long-lived perennial Gramineae grass with dense basal tillers and soft leaves. It is used widely in cold and dry areas of Eurasia and North America to establish grazing pasture and is even used as an ideal plant for revegetation and ecological restoration. Plant architecture, especially tillering traits, is critical for bunch grasses in breeding programs, and these traits in plants are mostly quantitative traits. In this study, the genetic diversity, population structure, and linkage disequilibrium of 480 individual lines were analyzed using 127 pairs of the EST-SSR marker, and a significant association between ten plant-architecture-related traits of P. juncea and molecular markers was found. The results of the genetic diversity analysis showed that the number of observed alleles was 1.957, the number of effective alleles was 1.682, Shannon's information index was 0.554, observed heterozygosity was 0.353, expected heterozygosity was 0.379, and the polymorphism information content was 0.300. A total of 480 individual lines were clustered into five groups based on population genetic structure, principal coordinate analysis, and unweighted pair group method with arithmetic mean analysis (UPGMA). The linkage disequilibrium coefficient (r2) was between 0.00 and 0.68, with an average of 0.04, which indicated a relatively low level of linkage disequilibrium among loci. The results of the association analysis revealed 55 significant marker-trait associations (MTA). Moreover, nine SSR markers were associated with multiple traits. This study provides tools with promising applications in the molecular selection and breeding of P. juncea germplasm.

A comprehensive overview of cotton genomics, biotechnology and molecular biological studies

Cotton is an irreplaceable economic crop currently domesticated in the human world for its extremely elongated fiber cells specialized in seed epidermis, which makes it of high research and application value. To date, numerous research on cotton has navigated various aspects, from multi-genome assembly, genome editing, mechanism of fiber development, metabolite biosynthesis, and analysis to genetic breeding. Genomic and 3D genomic studies reveal the origin of cotton species and the spatiotemporal asymmetric chromatin structure in fibers. Mature multiple genome editing systems, such as CRISPR/Cas9, Cas12 (Cpf1) and cytidine base editing (CBE), have been widely used in the study of candidate genes affecting fiber development. Based on this, the cotton fiber cell development network has been preliminarily drawn. Among them, the MYB-bHLH-WDR (MBW) transcription factor complex and IAA and BR signaling pathway regulate the initiation; various plant hormones, including ethylene, mediated regulatory network and membrane protein overlap fine-regulate elongation. Multistage transcription factors targeting CesA 4, 7, and 8 specifically dominate the whole process of secondary cell wall thickening. And fluorescently labeled cytoskeletal proteins can observe real-time dynamic changes in fiber development. Furthermore, research on the synthesis of cotton secondary metabolite gossypol, resistance to diseases and insect pests, plant architecture regulation, and seed oil utilization are all conducive to finding more high-quality breeding-related genes and subsequently facilitating the cultivation of better cotton varieties. This review summarizes the paramount research achievements in cotton molecular biology over the last few decades from the above aspects, thereby enabling us to conduct a status review on the current studies of cotton and provide strong theoretical support for the future direction.

Quantitative trait loci and candidate genes for yield-related traits of upland cotton revealed by genome-wide association analysis under drought conditions

BACKGROUND

Due to the influence of extreme weather, the environment in China's main cotton-producing areas is prone to drought stress conditions, which affect the growth and development of cotton and lead to a decrease in cotton yield.

RESULTS

In this study, 188 upland cotton germplasm resources were phenotyped for data of 8 traits (including 3 major yield traits) under drought conditions in three environments for two consecutive years. Correlation analysis revealed significant positive correlations between the three yield traits. Genetic analysis showed that the estimated heritability of the seed cotton index (SC) under drought conditions was the highest (80.81%), followed by that of boll weight (BW) (80.64%) and the lint cotton index (LC) (70.49%) With genome-wide association study (GWAS) analysis, a total of 75 quantitative trait loci (QTLs) were identified, including two highly credible new QTL hotspots. Three candidate genes (Gh_D09G064400, Gh_D10G261000 and Gh_D10G254000) located in the two new QTL hotspots, QTL51 and QTL55, were highly expressed in the early stage of fiber development and showed significant correlations with SC, LC and BW. The expression of three candidate genes in two extreme materials after drought stress was analyzed by qRT-PCR, and the expression of these two materials in fibers at 15, 20 and 25 DPA. The expression of these three candidate genes was significantly upregulated after drought stress and was significantly higher in drought-tolerant materials than in drought-sensitive materials. In addition, the expression levels of the three candidate genes were higher in the early stage of fiber development (15 DPA), and the expression levels in drought-tolerant germplasm were higher than those in drought-sensitive germplasm. These three candidate genes may play an important role in determining cotton yield under drought conditions.

CONCLUSIONS

This study is helpful for understanding the regulatory genes affecting cotton yield under drought conditions and provides germplasm and candidate gene resources for breeding high-yield cotton varieties under these conditions.

Genetic dissection of lint percentage in short-season cotton using combined QTL mapping and RNA-seq

KEY MESSAGE

In total, 17 QTLs for lint percentage in short-season cotton, including three stable QTLs, were detected. Twenty-eight differentially expressed genes located within the stable QTLs were identified, and two genes were validated by qRT-PCR. The breeding and use of short-season cotton have significant values in addressing the question of occupying farmlands with either cotton or cereals. However, the fiber yields of short-season cotton varieties are significantly lower than those of middle- and late-maturing varieties. How to effectively improve the fiber yield of short-season cotton has become a focus of cotton research. Here, a high-density genetic map was constructed using genome resequencing and an RIL population generated from the hybridization of two short-season cotton accessions, Dong3 and Dong4. The map contained 4960 bin markers across the 26 cotton chromosomes and spanned 3971.08 cM, with an average distance of 0.80 cM between adjacent markers. Based on the genetic map, quantitative trait locus (QTL) mapping for lint percentage (LP, %), an important yield component trait, was performed. In total, 17 QTLs for LP, including three stable QTLs, qLP-A02, qLP-D04, and qLP-D12, were detected. Three out of 11 non-redundant QTLs overlapped with previously reported QTLs, whereas the other eight were novel QTLs. A total of 28 differentially expressed genes associated with the three stable QTLs were identified using RNA-seq of ovules and fibers at different seed developmental stages from the parental materials. The two genes, Ghir_A02G017640 and Ghir_A02G018500, may be related to LP as determined by further qRT-PCR validation. This study provides useful information for the genetic dissection of LP and promotes the molecular breeding of short-season cotton.

GhAP1-D3 positively regulates flowering time and early maturity with no yield and fiber quality penalties in upland cotton

Flowering time (FTi) is a major factor determining how quickly cotton plants reach maturity. Early maturity greatly affects lint yield and fiber quality and is crucial for mechanical harvesting of cotton in northwestern China. Yet, few quantitative trait loci (QTLs) or genes regulating early maturity have been reported in cotton, and the underlying regulatory mechanisms are largely unknown. In this study, we characterized 152, 68, and 101 loci that were significantly associated with the three key early maturity traits-FTi, flower and boll period (FBP) and whole growth period (WGP), respectively, via four genome-wide association study methods in upland cotton (Gossypium hirsutum). We focused on one major early maturity-related genomic region containing three single nucleotide polymorphisms on chromosome D03, and determined that GhAP1-D3, a gene homologous to Arabidopsis thaliana APETALA1 (AP1), is the causal locus in this region. Transgenic plants overexpressing GhAP1-D3 showed significantly early flowering and early maturity without penalties for yield and fiber quality compared to wild-type (WT) plants. By contrast, the mutant lines of GhAP1-D3 generated by genome editing displayed markedly later flowering than the WT. GhAP1-D3 interacted with GhSOC1 (SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1), a pivotal regulator of FTi, both in vitro and in vivo. Changes in GhAP1-D3 transcript levels clearly affected the expression of multiple key flowering regulatory genes. Additionally, DNA hypomethylation and high levels of H3K9ac affected strong expression of GhAP1-D3 in early-maturing cotton cultivars. We propose that epigenetic modifications modulate GhAP1-D3 expression to positively regulate FTi in cotton through interaction of the encoded GhAP1 with GhSOC1 and affecting the transcription of multiple flowering-related genes. These findings may also lay a foundation for breeding early-maturing cotton varieties in the future.

Genome-wide association study identifies GhSAL1 affects cold tolerance at the seedling emergence stage in upland cotton (Gossypium hirsutum L.)

Genomic analysis of upland cotton revealed that cold tolerance was associated with ecological distribution. GhSAL1 on chromosome D09 negatively regulated cold tolerance of upland cotton. Cotton can undergo low-temperature stress at the seedling emergence stage, which adversely affects growth and yield; however, the regulatory mechanism underlying cold tolerance remains nebulous. Here, we analyze the phenotypic and physiological parameters in 200 accessions from 5 ecological distributions under constant chilling (CC) and diurnal variation of chilling (DVC) stresses at the seedling emergence stage. All accessions were clustered into four groups, of which Group IV, with most germplasms from the northwest inland region (NIR), had better phenotypes than Groups I-III under the two kinds of chilling stresses. A total of 575 significantly associated single-nucleotide polymorphism (SNP) were identified, and 35 stable genetic quantitative trait loci (QTL) were obtained, of which 5 were associated with traits under CC and DVC stress, respectively, while the remaining 25 were co-associated. The accumulation of dry weight (DW) of seedling was associated with the flavonoid biosynthesis process regulated by Gh_A10G0500. The emergence rate (ER), DW, and total length of seedling (TL) under CC stress were associated with the SNPs variation of Gh_D09G0189 (GhSAL1). GhSAL1^HapB was the elite haplotype, which increased ER, DW, and TL by 19.04%, 11.26%, and 7.69%, respectively, compared with that of GhSAL1^HapA. The results of virus-induced gene silencing (VIGS) experiment and determination of metabolic substrate content preliminarily illustrated that GhSAL1 negatively regulated cotton cold tolerance through IP₃-Ca²⁺ signaling pathway. The elite haplotypes and candidate genes identified in this study could be used to improve cold tolerance at the seedling emergence stage in future upland cotton breeding.

Recent advances and future perspectives in early-maturing cotton research

Cotton's fundamental requirements for long periods of growth and specific seasonal temperatures limit the global arable areas that can be utilized to cultivate cotton. This constraint can be alleviated by breeding for early-maturing varieties. By delaying the sowing dates without impacting the boll-opening time, early-maturing varieties not only mitigate the yield losses brought on by unfavorable weathers in early spring and late autumn but also help reducing the competition between cotton and other crops for arable land, thereby optimizing the cropping system. This review presents studies and breeding efforts for early-maturing cotton, which efficiently pyramid early maturity, high-quality, multiresistance traits, and suitable plant architecture by leveraging pleiotropic genes. Attempts are also made to summarize our current understanding of the molecular mechanisms underlying early maturation, which involves many pathways such as epigenetic, circadian clock, and hormone signaling pathways. Moreover, new avenues and effective measures are proposed for fine-scale breeding of early-maturing crops to ensure the healthy development of the agricultural industry.

QTL mapping and BSA-seq map a major QTL for the node of the first fruiting branch in cotton

Understanding the genetic basis of the node of the first fruiting branch (NFFB) improves early-maturity cotton breeding. Here we report QTL mapping on 200 F2 plants and derivative F2:3 and F2:4 populations by genotyping by sequencing (GBS). BC1F2 population was constructed by backcrossing one F2:4 line with the maternal parent JF914 and used for BSA-seq for further QTL mapping. A total of 1,305,642 SNPs were developed between the parents by GBS, and 2,907,790 SNPs were detected by BSA-seq. A high-density genetic map was constructed containing 11,488 SNPs and spanning 4,202.12 cM in length. A total of 13 QTL were mapped in the 3 tested populations. JF914 conferred favorable alleles for 11 QTL, and JF173 conferred favorable alleles for the other 2 QTL. Two stable QTL were repeatedly mapped in F2:3 and F2:4, including qNFFB-D3-1 and qNFFB-D6-1. Only qNFFB-D3-1 contributed more than 10% of the phenotypic variation. This QTL covered about 24.7 Mb (17,130,008-41,839,226 bp) on chromosome D3. Two regions on D3 (41,779,195-41,836,120 bp, 41,836,768-41,872,287 bp) were found by BSA-seq and covered about 92.4 Kb. This 92.4 Kb region overlapped with the stable QTL qNFFB-D3-1 and contained 8 annotated genes. By qRT-PCR, Ghir_D03G012430 showed a lower expression level from the 1- to 2-leaf stage and a higher expression level from the 3- to 6-leaf stage in the buds of JF173 than that of JF914. Ghir_D03G012390 reached the highest level at the 3- and 5-leaf stages in the buds of JF173 and JF914, respectively. As JF173 has lower NFFB and more early maturity than JF914, these two genes might be important in cell division and differentiation during NFFB formation in the seedling stage. The results of this study will facilitate a better understanding of the genetic basis of NFFB and benefit cotton molecular breeding for improving earliness traits.

Status and prospects of genome-wide association studies in cotton

Over the last two decades, the use of high-density SNP arrays and DNA sequencing have allowed scientists to uncover the majority of the genotypic space for various crops, including cotton. Genome-wide association study (GWAS) links the dots between a phenotype and its underlying genetics across the genomes of populations. It was first developed and applied in the field of human disease genetics. Many areas of crop research have incorporated GWAS in plants and considerable literature has been published in the recent decade. Here we will provide a comprehensive review of GWAS studies in cotton crop, which includes case studies on biotic resistance, abiotic tolerance, fiber yield and quality traits, current status, prospects, bottlenecks of GWAS and finally, thought-provoking question. This review will serve as a catalog of GWAS in cotton and suggest new frontiers of the cotton crop to be studied with this important tool.

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

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

Genome-Wide Dissection of the Genetic Basis for Drought Tolerance in Gossypium hirsutum L. Races

Drought seriously threats the growth and development of Gossypium hirsutum L. To dissect the genetic basis for drought tolerance in the G. hirsutum L. germplasm, a population, consisting of 188 accessions of G. hirsutum races and a cultivar (TM-1), was genotyped using the Cotton80KSNP biochip, and 51,268 high-quality single-nucleotide polymorphisms (SNPs) were obtained. Based on the phenotypic data of eight drought relative traits from four environments, we carried out association mapping with five models using GAPIT software. In total, thirty-six SNPs were detected significantly associated at least in two environments or two models. Among these SNPs, 8 and 28 (including 24 SNPs in 5 peak regions) were distributed in the A and D subgenome, respectively; eight SNPs were found to be distributed within separate genes. An SNP, TM73079, located on chromosome D10, was simultaneously associated with leaf fresh weight, leaf wilted weight, and leaf dry weight. Another nine SNPs, TM47696, TM33865, TM40383, TM10267, TM59672, TM59675, TM59677, TM72359, and TM72361, on chromosomes A13, A10, A12, A5, D6, and D9, were localized within or near previously reported quantitative trait loci for drought tolerance. Moreover, 520 genes located 200 kb up- and down-stream of 36 SNPs were obtained and analyzed based on gene annotation and transcriptome sequencing. The results showed that three candidate genes, Gh_D08G2462, Gh_A03G0043, and Gh_A12G0369, may play important roles in drought tolerance. The current GWAS represents the first investigation into mapping QTL for drought tolerance in G. hirsutum races and provides important information for improving cotton cultivars.

The Pivotal Role of Major Chromosomes of Sub-Genomes A and D in Fiber Quality Traits of Cotton

Lack of precise information about the candidate genes involved in a complex quantitative trait is a major obstacle in the cotton fiber quality improvement, and thus, overall genetic gain in conventional phenotypic selection is low. Recent molecular interventions and advancements in genome sequencing have led to the development of high-throughput molecular markers, quantitative trait locus (QTL) fine mapping, and single nucleotide polymorphisms (SNPs). These advanced tools have resolved the existing bottlenecks in trait-specific breeding. This review demonstrates the significance of chromosomes 3, 7, 9, 11, and 12 of sub-genomes A and D carrying candidate genes for fiber quality. However, chromosome 7 carrying SNPs for stable and potent QTLs related to fiber quality provides great insights for fiber quality-targeted research. This information can be validated by marker-assisted selection (MAS) and transgene in Arabidopsis and subsequently in cotton.

Identification of small effect quantitative trait loci of plant architectural, flowering, and early maturity traits in reciprocal interspecific introgression population in cotton

Plant architecture, flowering time and maturity traits are important determinants of yield and fiber quality of cotton. Genetic dissection of loci determining these yield and quality components is complicated by numerous loci with alleles conferring small differences. Therefore, mapping populations segregating for smaller numbers and sizes of introgressed segments is expected to facilitate dissection of these complex quantitative traits. At an advanced stage in the development of reciprocal advanced backcross populations from crosses between elite Gossypium hirsutum cultivar 'Acala Maxxa' (GH) and G. barbadense 'Pima S6' (GB), we undertook mapping of plant architectural traits, flowering time and maturity. A total of 284 BC₄F₁ and BC₄F₂ progeny rows, 120 in GH and 164 in GB background, were evaluated for phenotype, with only 4 and 3 (of 7) traits showing significant differences among progenies. Genotyping by sequencing yielded 3,186 and 3,026 SNPs, respectively, that revealed a total of 27 QTLs in GH background and 22 in GB, for plant height, days to flowering, residual flowering at maturity and maturity. More than of 90% QTLs identified in both backgrounds had small effects (%PV < 10), supporting the merit of this population structure to reduce background noise and small effect QTLs. Germplasm developed in this study may serve as potential pre-breeding material to develop improved cotton cultivars.

QTL and candidate gene identification of the node of the first fruiting branch (NFFB) by QTL-seq in upland cotton (Gossypium hirsutum L.)

Background

The node of the first fruiting branch (NFFB) is an important precocious trait in cotton. Many studies have been conducted on the localization of quantitative trait loci (QTLs) and genes related to fiber quality and yield, but there has been little attention to traits related to early maturity, especially the NFFB, in cotton.

Results

To identify the QTL associated with the NFFB in cotton, a BC₄F₂ population comprising 278 individual plants was constructed. The parents and two DNA bulks for high and low NFFB were whole genome sequenced, and 243.8 Gb of clean nucleotide data were generated. A total of 449,302 polymorphic SNPs and 135,353 Indels between two bulks were identified for QTL-seq. Seventeen QTLs were detected and localized on 11 chromosomes in the cotton genome, among which two QTLs (qNFFB-Dt2–1 and qNFFB-Dt3–3) were located in hotspots. Two candidate genes (GhAPL and GhHDA5) related to the NFFB were identified using quantitative real-time PCR (qRT-PCR) and virus-induced gene silencing (VIGS) experiments in this study. Both genes exhibited higher expression levels in the early-maturing cotton material RIL182 during flower bud differentiation, and the silencing of GhAPL and GhHDA5 delayed the flowering time and increased the NFFB compared to those of VA plants in cotton.

Conclusions

Our study preliminarily found that GhAPL and GhHDA5 are related to the early maturity in cotton. The findings provide a basis for the further functional verification of candidate genes related to the NFFB and contribute to the study of early maturity in cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08164-2.

Identification of hub genes through co-expression network of major QTLs of fiber length and strength traits in multiple RIL populations of cotton

The present study demonstrated a de novo correlation among fiber quality genes in multiple RIL populations including sGK9708 × 0-153, LMY22 × LY343 and Lumianyan28 × Xinluzao24. The current study was conducted to identify the major common QTLs including fiber length and strength, and to identify the co-expression networks of fiber length and strength QTLs harbored genes to target the hub genes. The RNA-seq data of sGK9708 × 0-153 population highlighted 50 and 48 candidate genes of fiber length and fiber strength QTLs. A total of 29 and 21 hub genes were identified in fiber length and strength co-expression network modules. The absolute values of correlation coefficient close to 1 resulted highly positive correlation among hub genes. Results also suggested that the gene correlation significantly influence the gene expression at different fiber development stages. These results might provide useful reference for further experiments in multiple RIL populations and suggest potential candidate genes for functional studies in cotton.

Genome-wide association analysis identified molecular markers associated with important tea flavor-related metabolites

The characteristic secondary metabolites in tea (theanine, caffeine, and catechins) are important factors contributing to unique tea flavors. However, there has been relatively little research on molecular markers related to these metabolites. Thus, we conducted a genome-wide association analysis of the levels of these tea flavor-related metabolites in three seasons. The theanine, caffeine, and catechin levels in Population 1 comprising 191 tea plant germplasms were examined, which revealed that their heritability exceeded 0.5 in the analyzed seasons, with the following rank order (highest to lowest heritabilities): (+)-catechin > (-)-gallocatechin gallate > caffeine = (-)-epicatechin > (-)-epigallocatechin-3-gallate > theanine > (-)-epigallocatechin > (-)-epicatechin-3-gallate > catechin gallate > (+)-gallocatechin. The SNPs detected by amplified-fragment SNP and methylation sequencing divided Population 1 into three groups and seven subgroups. An association analysis yielded 307 SNP markers related to theanine, caffeine, and catechins that were common to all three seasons. Some of the markers were pleiotropic. The functional annotation of 180 key genes at the SNP loci revealed that FLS, UGT, MYB, and WD40 domain-containing proteins, as well as ATP-binding cassette transporters, may be important for catechin synthesis. KEGG and GO analyses indicated that these genes are associated with metabolic pathways and secondary metabolite biosynthesis. Moreover, in Population 2 (98 tea plant germplasm resources), 30 candidate SNPs were verified, including 17 SNPs that were significantly or extremely significantly associated with specific metabolite levels. These results will provide a foundation for future research on important flavor-related metabolites and may help accelerate the breeding of new tea varieties.

Genome-wide association studies of fruit quality traits in jujube germplasm collections using genotyping-by-sequencing

Chinese jujube (Ziziphus jujuba Mill.) is an important fruit crop and harbors many highly diverse traits of potential economic importance. Fruit size, stone size, and fruit cracking have an important influence on the commercial value of jujube. This study is the first to conduct a genome-wide association study (GWAS) on 180 accessions of jujube and focuses on locating single-nucleotide polymorphisms (SNPs) associated with nine important fruit quality traits. Genotyping was performed using genotyping-by-sequencing and 4651 high-quality SNPs were identified. A genetic diversity analysis revealed the presence of three distinct groups, and rapid linkage disequilibrium decay was observed in this jujube population. Using a mixed linear model, a total of 45 significant SNP-trait associations were detected, among which 33 SNPs had associations with fruit size-related traits, nine were associated with stone size-related traits, and three with fruit cracking-related traits. In total, 21 candidate genes involved in cell expansion, abiotic stress responses, hormone signaling, and growth development were identified from the genome sequences of jujube. These results are useful as basic data for GWAS of other jujube traits, and these significant SNP loci and candidate genes should aid marker-assisted breeding and genomic selection of improved jujube cultivars.

Genome-wide association study reveals the genetic basis of fiber quality traits in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Fiber quality is an important economic trait of cotton, and its improvement is a major goal of cotton breeding. To better understand the genetic mechanisms responsible for fiber quality traits, we conducted a genome-wide association study to identify and mine fiber-quality-related quantitative trait loci (QTLs) and genes.

RESULTS

In total, 42 single nucleotide polymorphisms (SNPs) and 31 QTLs were identified as being significantly associated with five fiber quality traits. Twenty-five QTLs were identified in previous studies, and six novel QTLs were firstly identified in this study. In the QTL regions, 822 genes were identified and divided into four clusters based on their expression profiles. We also identified two pleiotropic SNPs. The SNP locus i52359Gb was associated with fiber elongation, strength, length and uniformity, while i11316Gh was associated with fiber strength and length. Moreover, these two SNPs were nonsynonymous and located in genes Gh_D09G2376 and Gh_D06G1908, respectively. RT-qPCR analysis revealed that these two genes were preferentially expressed at one or more stages of cotton fiber development, which was consistent with the RNA-seq data. Thus, Gh_D09G2376 and Gh_D06G1908 may be involved in fiber developmental processes.

CONCLUSIONS

The findings of this study provide insights into the genetic bases of fiber quality traits, and the identified QTLs or genes may be applicable in cotton breeding to improve fiber quality.

Dissecting the genetic architecture of seed-cotton and lint yields in Upland cotton using genome-wide association mapping

Seed-cotton yield (SY) and lint yield (LY) are the most important yield traits of cotton. Thus, it is critical to dissect their genetic architecture. Upland cotton (Gossypium hirsutum) is widely grown worldwide. In this study, a genome-wide association mapping was performed based on the CottonSNP80K array to dissect the genetic architecture of SY and LY in Upland cotton. Twenty-three significant associations were detected within four environments, including 11 associated with SY and 12 associated with LY. Seven single nucleotide polymorphisms (SNPs), TM234, TM237, TM247, TM255, TM256, TM263, and TM264, were co-associated with the two traits, which may indicate pleiotropy or intergenic tight linkages. Five SNPs, TM13332, TM39771, TM57119, TM81653, and TM81660, were coincided with those of previous reports and could be used in marker-assisted selection. Combining functional annotations with expression analyses of the genes identified within 400 kb of the significantly associated SNPs, we hypothesize that the three genes, Gh_D05G1077 and Gh_D13G1571 for SY, and Gh_A11G0775 for LY, may have the potential to increase cotton yield. The results would provide useful information for understanding the genetic basis of yield traits in Upland cotton and for facilitating its high-yield breeding through molecular design.

Dissection of the genetic variation and candidate genes of lint percentage by a genome-wide association study in upland cotton

A genome-wide associated study identified six novel QTLs for lint percentage. Two candidate genes underlying this trait were also detected. Increasing lint percentage (LP) is a core goal of cotton breeding. To better understand the genetic basis of LP, a genome-wide association study (GWAS) was conducted using 276 upland cotton accessions planted in multiple environments and genotyped with a CottonSNP63K array. After filtering, 10,660 high-quality single-nucleotide polymorphisms (SNPs) were retained. Population structure, principal component and neighbor-joining phylogenetic tree analyses divided the accessions into two subpopulations. These results along with linkage disequilibrium decay indicated accessions were not highly structured and exhibited weak relatedness. GWAS uncovered 23 polymorphic SNPs and 15 QTLs significantly associated with LP, with six new QTLs identified. Two candidate genes, Gh_D05G0313 and Gh_D05G1124, both contained one significant SNP, highly expressed during ovule and fiber development stages, implying that the two genes may act as the most promising regulators of LP. Furthermore, the phenotypic value of LP was found to be positively correlated with the number of favorable SNP alleles. These favorable alleles for LP identified in the study may be useful for improving lint yield.