Genomic analyses reveal the genetic basis of early maturity and identification of loci and candidate genes in upland cotton (Gossypium hirsutum L.)

Genome-Wide Identification and Bioinformatics Analysis of the FK506 Binding Protein Family in Rice

The FK506 Binding Protein (FKBP), ubiquitously present across diverse species, is characterized by its evolutionarily conserved FK506 binding domain (FKBd). In plants, evidence suggests that this gene family plays integral roles in regulating growth, development, and responses to environmental stresses. Notably, research on the identification and functionality of FKBP genes in rice remains limited. Therefore, this study utilized bioinformatic tools to identify 30 FKBP-encoding genes in rice. It provides a detailed analysis of their chromosomal locations, evolutionary relationships with the Arabidopsis thaliana FKBP family, and gene structures. Further analysis of the promoter elements of these rice FKBP genes revealed a high presence of stress-responsive elements. Quantitative PCR assays under drought and heat stress conditions demonstrated that genes OsFKBP15-2, OsFKBP15-3, OsFKBP16-3, OsFKBP18, and OsFKBP42b are inducible by these adverse conditions. These findings suggest a significant role for the rice FKBP gene family in stress adaptation. This research establishes a critical foundation for deeper explorations of the functional roles of the OsFKBP genes in rice.

Genome-wide identification and mining elite allele variation of the Monoacylglycerol lipase (MAGL) gene family in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Monoacylglycerol lipase (MAGL) genes belong to the alpha/beta hydrolase superfamily, catalyze the terminal step of triglyceride (TAG) hydrolysis, converting monoacylglycerol (MAG) into free fatty acids and glycerol.

RESULTS

In this study, 30 MAGL genes in upland cotton have been identified, which have been classified into eight subgroups. The duplication of GhMAGL genes in upland cotton was predominantly influenced by segmental duplication events, as revealed through synteny analysis. Furthermore, all GhMAGL genes were found to contain light-responsive elements. Through comprehensive association and haplotype analyses using resequencing data from 355 cotton accessions, GhMAGL3 and GhMAGL6 were detected as key genes related to lipid hydrolysis processes, suggesting a negative regulatory effect.

CONCLUSIONS

In summary, MAGL has never been studied in upland cotton previously. This study provides the genetic mechanism foundation for the discover of new genes involved in lipid metabolism to improve cottonseed oil content, which will provide a strategic avenue for marker-assisted breeding aimed at incorporating desirable traits into cultivated cotton varieties.

Multistage temporal transcriptomic atlas unveils major contributor to reproductive phase in upland cotton

Flowering is a major developmental transition in plants, but asynchronous flowering hinders the utilization of wild cotton relatives in breeding programs. We performed comparative transcriptomic profiling of early- and late-flowering Gossypium hirsutum genotypes to elucidate genetic factors influencing reproductive timing. Shoot apices were sampled from the photoperiod-sensitive landrace G. hirsutum purpurascens (GhP) and early-maturing variety ZhongMianSuo (ZMS) at five time points following the emergence of sympodial nodes. RNA-sequencing revealed extensive transcriptional differences during floral transition. Numerous flowering-associated genes exhibited genotype-specific expression, including FLOWERING LOCUS T (FT) homologs upregulated in ZMS. FT-interacting factors like SOC1 and CO-like also showed higher expression in ZMS, implicating florigen pathways in early flowering. Additionally, circadian clock and light signalling components were misregulated between varieties, suggesting altered photoperiod responses in GhP. Weighted co-expression network analysis specifically linked a module enriched for circadian-related genes to GhP's late flowering. Through an integrated transcriptome analysis, we defined a regulatory landscape of reproductive phase change in cotton. Differentially expressed genes related to photoperiod, circadian clock, and light signalling likely contribute to delayed flowering in wild cottons. Characterization of upstream flowering regulators will enable modifying photoperiod sensitivity and expand germplasm use for cotton improvement. This study provides candidate targets for elucidating interactive mechanisms that control cotton flowering time across diverse genotypes.

Identification and characterization of ACR gene family in maize for salt stress tolerance

Background

Members of the ACR gene family are commonly involved in various physiological processes, including amino acid metabolism and stress responses. In recent decades, significant progress has been made in the study of ACR genes in plants. However, little is known about their characteristics and function in maize.

Methods

In this study, ACR genes were identified from the maize genome, and their molecular characteristics, gene structure, gene evolution, gene collinearity analysis, cis-acting elements were analyzed. qRT-PCR technology was used to verify the expression patterns of the ZmACR gene family in different tissues under salt stress. In addition, Ectopic expression technique of ZmACR5 in Arabidopsis thaliana was utilized to identify its role in response to salt stress.

Results

A total of 28 ZmACR genes were identified, and their molecular characteristics were extensively described. Two gene pairs arising from segmented replication events were detected in maize, and 18 collinear gene pairs were detected between maize and 3 other species. Through phylogenetic analysis, three subgroups were revealed, demonstrating distinct divergence between monocotyledonous and dicotyledonous plants. Analysis of ZmACR cis-acting elements revealed the optional involvement of ZmACR genes in light response, hormone response and stress resistance. Expression analysis of 8 ZmACR genes under salt treatment clearly revealed their role in the response to salt stress. Ectopic overexpression of ZmACR5 in Arabidopsis notably reduced salt tolerance compared to that of the wild type under salt treatment, suggesting that ZmACR5 has a negative role in the response to salt stress.

Conclusion

Taken together, these findings confirmed the involvement of ZmACR genes in regulating salt stress and contributed significantly to our understanding of the molecular function of ACR genes in maize, facilitating further research in this field.

Genomic loci associated with leaf abscission contribute to machine picking and environmental adaptability in upland cotton (Gossypium hirsutum L.)

INTRODUCTION

Defoliation by applying defoliants before machine picking is an important agricultural practice that enhances harvesting efficiency and leads to increased raw cotton purity. However, the fundamental characteristics of leaf abscission and the underlying genetic basis in cotton are not clearly understood.

OBJECTIVES

In this study, we aimed to (1) reveal the phenotypic variations in cotton leaf abscission, (2) discover the whole-genome differentiation sweeps and genetic loci related to defoliation, (3) identify and verify the functions of key candidate genes associated with defoliation, and (4) explore the relationship between haplotype frequency of loci and environmental adaptability.

METHODS

Four defoliation-related traits of 383 re-sequenced Gossypium hirsutum accessions were investigated in four environments. The genome-wide association study (GWAS), linkage disequilibrium (LD) interval genotyping and functional identification were conducted. Finally, the haplotype variation related to environmental adaptability and defoliation traits was revealed.

RESULTS

Our findings revealed the fundamental phenotypic variations of defoliation traits in cotton. We showed that defoliant significantly increased the defoliation rate without incurring yield and fiber quality penalties. The strong correlations between defoliation traits and growth period traits were observed. A genome-wide association study of defoliation traits identified 174 significant SNPs. Two loci (RDR7 on A02 and RDR13 on A13) that significantly associated with the relative defoliation rate were described, and key candidate genes GhLRR and GhCYCD3;1, encoding a leucine-rich repeat (LRR) family protein and D3-type cell cyclin 1 protein respectively, were functional verified by expression pattern analysis and gene silencing. We found that combining of two favorable haplotypes (HapRDR7 and HapRDR13) improved sensitivity to defoliant. The favorable haplotype frequency generally increased in high latitudes in China, enabling adaptation to the local environment.

CONCLUSION

Our findings lay an important foundation for the potentially broad application of leveraging key genetic loci in breeding machine-pickable cotton.

Genetic interrogation of phenotypic plasticity informs genome-enabled breeding in cotton

Phenotypic plasticity, or the ability to adapt to and thrive in changing climates and variable environments, is essential for developmental programs in plants. Despite its importance, the genetic underpinnings of phenotypic plasticity for key agronomic traits remain poorly understood in many crops. In this study, we aim to fill this gap by using genome-wide association studies to identify genetic variations associated with phenotypic plasticity in upland cotton (Gossypium hirsutum L.). We identified 73 additive quantitative trait loci (QTLs), 32 dominant QTLs, and 6799 epistatic QTLs associated with 20 traits. We also identified 117 additive QTLs, 28 dominant QTLs, and 4691 epistatic QTLs associated with phenotypic plasticity in 19 traits. Our findings reveal new genetic factors, including additive, dominant, and epistatic QTLs, that are linked to phenotypic plasticity and agronomic traits. Meanwhile, we find that the genetic factors controlling the mean phenotype and phenotypic plasticity are largely independent in upland cotton, indicating the potential for simultaneous improvement. Additionally, we envision a genomic design strategy by utilizing the identified QTLs to facilitate cotton breeding. Taken together, our study provides new insights into the genetic basis of phenotypic plasticity in cotton, which should be valuable for future breeding.

Synergistic optimization of crops by combining early maturation with other agronomic traits

Many newly created early maturing varieties exhibit poor stress resistance and low yield, whereas stress-resistant varieties are typically late maturing. For this reason, the polymerization of early maturity and other desired agronomic qualities requires overcoming the negative connection between early maturity, multi-resistance, and yield, which presents a formidable challenge in current breeding techniques. We review the most salient constraints of early maturity breeding in current crop planting practices and the molecular mechanisms of different maturation timeframes in diverse crops from their origin center to production areas. We explore current breeding tactics and the future direction of crop breeding and the issues that must be resolved to accomplish the polymerization of desirable traits in light of the current obstacles and limitations.

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.

Identification of novel candidate loci and genes for seed vigor-related traits in upland cotton (Gossypium hirsutum L.) via GWAS

Seed vigor (SV) is a crucial trait determining the quality of crop seeds. Currently, over 80% of China's cotton-planting area is in Xinjiang Province, where a fully mechanized planting model is adopted, accounting for more than 90% of the total fiber production. Therefore, identifying SV-related loci and genes is crucial for improving cotton yield in Xinjiang. In this study, three seed vigor-related traits, including germination potential, germination rate, and germination index, were investigated across three environments in a panel of 355 diverse accessions based on 2,261,854 high-quality single-nucleotide polymorphisms (SNPs). A total of 26 significant SNPs were detected and divided into six quantitative trait locus regions, including 121 predicted candidate genes. By combining gene expression, gene annotation, and haplotype analysis, two novel candidate genes (Ghir_A09G002730 and Ghir_D03G009280) within qGR-A09-1 and qGI/GP/GR-D03-3 were associated with vigor-related traits, and Ghir_A09G002730 was found to be involved in artificial selection during cotton breeding by population genetic analysis. Thus, understanding the genetic mechanisms underlying seed vigor-related traits in cotton could help increase the efficiency of direct seeding by molecular marker-assisted selection breeding.

Genome-wide identification of R-SNARE gene family in upland cotton and function analysis of GhVAMP72l response to drought stress

Soluble N-ethylmaleimide-sensitive factor attachment protein receptors (R-SNAREs) mainly promoted the assembly of the SNARE complex to drive the final membrane fusion step of membrane transport. Previous research on R-SNAREs has mainly focused on development and growth and has rarely been involved in abiotic stress, especially in cotton. Here, we performed a comprehensive analysis of R-SNARE genes in upland cotton. In total, 51 Gh-R-SNARE genes across six phylogenetic groups were unevenly distributed on 21 chromosomes. Cis elements related to plant growth and response to abiotic stress responses were found in the promoter region of Gh-R-SNAREs. Nine Gh-R-SNARE genes were obviously upregulated under drought stress conditions by RNA-seq and qRT-PCR analysis. Among them, GhVAMP72l might be the key candidate gene contributing to drought stress tolerance in cotton by virus-induced gene silencing (VIGS) assay. These results provide valuable insights for the functional analysis of cotton R-SNAREs in response to drought stress and highlight potential beneficial genes for genetic improvement and breeding in 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.

First Report on the Genetic Diversity of Populations of Gossypium barbadense L. and Gossypium hirsutun L. in the Amazonian Native Communities, Cusco-Peru

The genus Gossypium has important ethnobotanical and economic value for Amazonian Native Communities (A.N.C.). However, little research has been undertaken on the distribution and genetic diversity of cotton populations maintained in the Peruvian rainforest. This work aims to present the first report on the genetic diversity of Gossypium spp. populations in the A.N.C. of the province of La Convención, Cusco-Peru. The methodology was based on exploring, collecting, identifying, and characterizing the Gossypium populations present in the A.N.C. Twenty-six descriptors were evaluated (9 quantitative and 17 qualitative), and with this information, distribution, correlation, and principal component (PC) analyses were carried out. As a result, plants of two species [G. barbadense L. (44 samples) and G. hirsutum L. (19 samples)], one variety [G. barbadense var. brasiliensis (75 samples)], and three previously unidentified variations (9 samples) were identified. Altogether, 147 samples were collected. G. barbadense var. brasiliensis, which was always found in association with other economic crops within an altitude range of 338 to 1086 m, was the most predominant (51%), distributed in eleven A.N.C. and always in small plots (up to 2 ha). G. barbadense L. was cultivated between 397 and 1137 m of altitude in eight A.N.C. in plots of up to 3 ha in marginal lands. G. hirsutum L., with a smaller distribution (13%), was found between 334 and 497 m of altitude in only three communities; this species is cultivated in marginal areas throughout the year. The variability found for the first two PCs when considering the quantitative and qualitative descriptors was high (74.7%) and moderate (48.2%), respectively. When combining all the descriptors, the analysis showed that the first two PCs accounted for 51.8% of the total variability of the data. The PCs of the two types of data and their combination confirmed that the three populations found were grouped. The nine undefined samples were close to or intermediate between the described ones, showing that these samples may be the result of spontaneous crosses; as such, these samples need to be better evaluated with other tools for further definition. The information obtained shows that in the A.N.C. of Cusco-Peru, there is variability conserved by the inhabitants, who have been able to maintain and use these genotypes, even from their Amazonian indigenous ancestry, and the environment has been able to generate variability among the species, as will be highlighted in future works.

Dissection of the Pearl of Csaba pedigree identifies key genomic segments related to early ripening in grape

Pearl of Csaba (PC) is a valuable backbone parent for early-ripening grapevine (Vitis vinifera) breeding, from which many excellent early ripening varieties have been bred. However, the genetic basis of the stable inheritance of its early ripening trait remains largely unknown. Here, the pedigree, consisting of 40 varieties derived from PC, was re-sequenced for an average depth of ∼30×. Combined with the resequencing data of 24 other late-ripening varieties, 5,795,881 high-quality single nucleotide polymorphisms (SNPs) were identified following a strict filtering pipeline. The population genetic analysis showed that these varieties could be distinguished clearly, and the pedigree was characterized by lower nucleotide diversity and stronger linkage disequilibrium than the non-pedigree varieties. The conserved haplotypes (CHs) transmitted in the pedigree were obtained via identity-by-descent analysis. Subsequently, the key genomic segments were identified based on the combination analysis of haplotypes, selective signatures, known ripening-related quantitative trait loci (QTLs), and transcriptomic data. The results demonstrated that varieties with a superior haplotype, H1, significantly (one-way ANOVA, P < 0.001) exhibited early grapevine berry development. Further analyses indicated that H1 encompassed VIT_16s0039g00720 encoding a folate/biopterin transporter protein (VvFBT) with a missense mutation. VvFBT was specifically and highly expressed during grapevine berry development, particularly at veraison. Exogenous folate treatment advanced the veraison of "Kyoho". This work uncovered core haplotypes and genomic segments related to the early ripening trait of PC and provided an important reference for the molecular breeding of early-ripening grapevine varieties.

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.

Recent progression and future perspectives in cotton genomic breeding

Upland cotton is an important global cash crop for its long seed fibers and high edible oil and protein content. Progress in cotton genomics promotes the advancement of cotton genetics, evolutionary studies, functional genetics, and breeding, and has ushered cotton research and breeding into a new era. Here, we summarize high-impact genomics studies for cotton from the last 10 years. The diploid Gossypium arboreum and allotetraploid Gossypium hirsutum are the main focus of most genetic and genomic studies. We next review recent progress in cotton molecular biology and genetics, which builds on cotton genome sequencing efforts, population studies, and functional genomics, to provide insights into the mechanisms shaping abiotic and biotic stress tolerance, plant architecture, seed oil content, and fiber development. We also suggest the application of novel technologies and strategies to facilitate genome-based crop breeding. Explosive growth in the amount of novel genomic data, identified genes, gene modules, and pathways is now enabling researchers to utilize multidisciplinary genomics-enabled breeding strategies to cultivate "super cotton", synergistically improving multiple traits. These strategies must rise to meet urgent demands for a sustainable cotton 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 F₂ plants and derivative F_2:3 and F_2:4 populations by genotyping by sequencing (GBS). BC₁F₂ population was constructed by backcrossing one F_2: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 F_2:3 and F_2: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.

SEP-like genes of Gossypium hirsutum promote flowering via targeting different loci in a concentration-dependent manner

SEP genes are famous for their function in the morphological novelty of bisexual flowers. Although the diverse functions of SEP genes were reported, only the regulatory mechanisms underlying floral organ development have been addressed. In this study, we identified SEP-like genes in Gossypium and found that SEP3 genes were duplicated in diploid cotton varieties. GhSEP4.1 and GhSEP4.2 were abundantly transcribed in the shoot apical meristem (SAM), but only GhSEP4.2 was expressed in the leaf vasculature. The expression pattern of GhSEPs in floral organs was conserved with that of homologs in Arabidopsis, except for GhSEP2 that was preponderantly expressed in ovules and fibers. The overexpression and silencing of each single GhSEP gene suggested their distinct role in promoting flowering via direct binding to GhAP1 and GhLFY genomic regions. The curly leaf and floral defects in overexpression lines with a higher expression of GhSEP genes revealed the concentration-dependent target gene regulation of GhSEP proteins. Moreover, GhSEP proteins were able to dimerize and interact with flowering time regulators. Together, our results suggest the dominant role of GhSEP4.2 in leaves to promote flowering via GhAP1-A04, and differently accumulated GhSEP proteins in the SAM alternately participate in forming the dynamic tetramer complexes to target at the different loci of GhAP1 and GhLFY to maintain reproductive growth. The regulatory roles of cotton SEP genes reveal their conserved and diversified functions.

GhCDPK60 positively regulates drought stress tolerance in both transgenic Arabidopsis and cotton by regulating proline content and ROS level

Calcium-Dependent Protein Kinases (CDPKs) involved in regulating downstream components of calcium signaling pathways play a role in tolerance to abiotic stresses and seed development in plants. However, functions of only a few cotton CDPKs have been clarified at present. In this study, 80 conserved CDPKs in Gossypium hirsutum L. were identified and characterized, which was divided into four subgroups. Among them, the transcript level of GhCDPK60 was significantly upregulated under drought and several hormone treatments. And we found that the expression levels of several stress-inducible genes down-regulated in GhCDPK60-silence cotton and up-regulated in GhCDPK60-overexpressing Arabidopsis. In addition, physiological analyses demonstrated that GhCDPK60 improved drought stress tolerance by improving the osmotic adjustment ability and reducing the accumulation of reactive oxygen species (ROS) in plants. These findings broaden our understanding of the biological roles of GhCDPK60 and mechanisms underlying drought stress tolerance in cotton.

Genome-wide association study reveals novel QTLs and candidate genes for seed vigor in rice

Highly seed vigor (SV) is essential for rice direct seeding (DS). Understanding the genetic mechanism of SV-related traits could contribute to increasing the efficiency of DS. However, only a few genes responsible for SV have been determined in rice, and the regulatory network of SV remains obscure. In this study, the seed germination rate (GR), seedling shoot length (SL), and shoot fresh weight (FW) related to SV traits were measured, and a genome-wide association study (GWAS) was conducted to detect high-quality loci responsible for SV using a panel of 346 diverse accessions. A total of 51 significant SNPs were identified and arranged into six quantitative trait locus (QTL) regions, including one (qGR1-1), two (qSL1-1, qSL1-2), and three (qFW1-1, qFW4-1, and qFW7-1) QTLs associated with GR, SL, and FW respectively, which were further validated using chromosome segment substitution lines (CSSLs). Integrating gene expression, gene annotation, and haplotype analysis, we found 21 strong candidate genes significantly associated with SV. In addition, the SV-related functions of LOC_Os01g11270 and LOC_Os01g55240 were further verified by corresponding CRISPR/Cas9 gene-edited mutants. Thus, these results provide clues for elucidating the genetic basis of SV control. The candidate genes or QTLs would be helpful for improving DS by molecular marker-assisted selection (MAS) breeding in rice.

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.

Identification of miRNAs and their target genes associated with improved maize seed vigor induced by gibberellin

High seed vigor is crucial for agricultural production owing to its potential in high quality and yield of crops and a better understanding of the molecular mechanism associated with maize seed vigor is highly necessary. To better understand the involvement and regulatory mechanism of miRNAs correlated with maize seed vigor, small RNAs and degradome sequencing of two inbred lines Yu537A and Yu82 were performed. A total of 791 mature miRNAs were obtained with different expressions, among of which 505 miRNAs were newly identified and the rest miRNAs have been reported before by comparing the miRNAs with the sequences in miRbase database. Analysis of miRNA families showed maize seeds contain fewer miRNA families and larger miRNA families compared with animals, indicating that functions of miRNAs in maize seeds were more synergistic than animals. Degradome sequencing was used to identify the targets of miRNAs and the results showed a total of 6,196 targets were obtained. Function analysis of differentially expressed miRNAs and targets showed Glycan degradation and galactose metabolism were closely correlated with improved maize seed vigor. These findings provide valuable information to understand the involvement of miRNAs with maize seed vigor and these putative genes will be valuable resources for improving the seed vigor in future maize breeding.

Genetic basis of sorghum leaf width and its potential as a surrogate for transpiration efficiency

Leaf width was correlated with plant-level transpiration efficiency and associated with 19 QTL in sorghum, suggesting it could be a surrogate for transpiration efficiency in large breeding program. Enhancing plant transpiration efficiency (TE) by reducing transpiration without compromising photosynthesis and yield is a desirable selection target in crop improvement programs. While narrow individual leaf width has been correlated with greater intrinsic water use efficiency in C₄ species, the extent to which this translates to greater plant TE has not been investigated. The aims of this study were to evaluate the correlation of leaf width with TE at the whole-plant scale and investigate the genetic control of leaf width in sorghum. Two lysimetry experiments using 16 genotypes varying for stomatal conductance and three field trials using a large sorghum diversity panel (n = 701 lines) were conducted. Negative associations of leaf width with plant TE were found in the lysimetry experiments, suggesting narrow leaves may result in reduced plant transpiration without trade-offs in biomass accumulation. A wide range in width of the largest leaf was found in the sorghum diversity panel with consistent ranking among sorghum races, suggesting that environmental adaptation may have a role in modifying leaf width. Nineteen QTL were identified by genome-wide association studies on leaf width adjusted for flowering time. The QTL identified showed high levels of correspondence with those in maize and rice, suggesting similarities in the genetic control of leaf width across cereals. Three a priori candidate genes for leaf width, previously found to regulate dorsoventrality, were identified based on a 1-cM threshold. This study provides useful physiological and genetic insights for potential manipulation of leaf width to improve plant adaptation to diverse environments.

Insights into the Genomic Regions and Candidate Genes of Senescence-Related Traits in Upland Cotton via GWAS

Senescence is the last stage of plant development and is controlled by both internal and external factors. Premature senescence significantly affects the yield and quality of cotton. However, the genetic architecture underlying cotton senescence remains unclear. In this study, genome-wide association studies (GWAS) were performed based on 3,015,002 high-quality SNP markers from the resequencing data of 355 upland cotton accessions to detect genomic regions for cotton senescence. A total of 977 candidate genes within 55 senescence-related genomic regions (SGRs), SGR1-SGR55, were predicted. Gene ontology (GO) analysis of candidate genes revealed that a set of biological processes was enriched, such as salt stress, ethylene processes, and leaf senescence. Furthermore, in the leaf senescence GO term, one candidate gene was focused on: Gohir.A12G270900 (GhMKK9), located in SGR36, which encodes a protein of the MAP kinase kinase family. Quantitative real-time PCR (qRT-PCR) analysis showed that GhMKK9 was up-regulated in old cotton leaves. Overexpression of GhMKK9 in Arabidopsis accelerated natural leaf senescence. Virus-induced gene silencing (VIGS) of GhMKK9 in cotton increased drought tolerance. These results suggest that GhMKK9 is a positive regulator and might be involved in drought-induced senescence in cotton. The results provide new insights into the genetic basis of cotton senescence and will be useful for improving cotton breeding in the future.

Genomewide Identification and Characterization of the Genes Involved in the Flowering of Cotton

Flowering is a prerequisite for flowering plants to complete reproduction, and flowering time has an important effect on the high and stable yields of crops. However, there are limited reports on flowering-related genes at the genomic level in cotton. In this study, genomewide analysis of the evolutionary relationship of flowering-related genes in different cotton species shows that the numbers of flowering-related genes in the genomes of tetraploid cotton species Gossypium hirsutum and Gossypium barbadense were similar, and that these numbers were approximately twice as much as the number in diploid cotton species Gossypium arboretum. The classification of flowering-related genes shows that most of them belong to the photoperiod and circadian clock flowering pathway. The distribution of flowering-related genes on the chromosomes of the At and Dt subgenomes was similar, with no subgenomic preference detected. In addition, most of the flowering-related core genes in Arabidopsis thaliana had homologs in the cotton genome, but the copy numbers and expression patterns were disparate; moreover, flowering-related genes underwent purifying selection throughout the evolutionary and selection processes. Although the differentiation and reorganization of many key genes of the cotton flowering regulatory network occurred throughout the evolutionary and selection processes, most of them, especially those involved in the important flowering regulatory networks, have been relatively conserved and preferentially selected.

Detection of Stable Elite Haplotypes and Potential Candidate Genes of Boll Weight Across Multiple Environments via GWAS in Upland Cotton

Boll weight (BW) is a key determinant of yield component traits in cotton, and understanding the genetic mechanism of BW could contribute to the progress of cotton fiber yield. Although many yield-related quantitative trait loci (QTLs) responsible for BW have been determined, knowledge of the genes controlling cotton yield remains limited. Here, association mapping based on 25,169 single-nucleotide polymorphisms (SNPs) and 2,315 insertions/deletions (InDels) was conducted to identify high-quality QTLs responsible for BW in a global collection of 290 diverse accessions, and BW was measured in nine different environments. A total of 19 significant markers were detected, and 225 candidate genes within a 400 kb region (± 200 kb surrounding each locus) were predicted. Of them, two major QTLs with highly phenotypic variation explanation on chromosomes A08 and D13 were identified among multiple environments. Furthermore, we found that two novel candidate genes (Ghir_A08G009110 and Ghir_D13G023010) were associated with BW and that Ghir_D13G023010 was involved in artificial selection during cotton breeding by population genetic analysis. The transcription level analyses showed that these two genes were significantly differentially expressed between high-BW accession and low-BW accession during the ovule development stage. Thus, these results reveal valuable information for clarifying the genetic basics of the control of BW, which are useful for increasing yield by molecular marker-assisted selection (MAS) breeding in cotton.

Characterization and Functional Analysis of GhNAC82, A NAM Domain Gene, Coordinates the Leaf Senescence in Upland Cotton (Gossypium hirsutum L.)

In the process of growth and development, cotton exhibits premature senescence under various abiotic stresses, impairing yield and fiber quality. NAC (NAM, ATAF1,2, and CUC2) protein widely distributed in land plants, play the critical role in responding to abiotic stress and regulating leaf senescence. We have identified and functional analyzed a NAM domain gene GhNAC82 in upland cotton, it was located on the A11 chromosome 4,921,702 to 4,922,748 bp, only containing one exon. The spatio-temporal expression pattern analysis revealed that it was highly expressed in root, torus, ovule and fiber development stage. The results of qRT-PCR validated that GhNAC82 negatively regulated by salt stress, drought stress, H₂O₂ stress, IAA treatment, and ethylene treatment, positively regulated by the ABA and MeJA treatment. Moreover, heterologous overexpression of GhNAC82 results in leaf premature senescence and delays root system development in Arabidopsis thaliana. The phenotype of delayed-senescence was performed after silencing GhNAC82 by VIGS in premature cotton. Taken together, GhNAC82 was involved in different abiotic stress pathways and play important roles in negatively regulating leaf premature senescence.

GhSOC1s Evolve to Respond Differently to the Environmental Cues and Promote Flowering in Partially Independent Ways

Gossypium hirsutum is most broadly cultivated in the world due to its broader adaptation to the environment and successful breeding of early maturity varieties. However, how cotton responds to environmental cues to adjust flowering time to achieve reproductive success is largely unknown. SOC1 functions as an essential integrator for the endogenous and exogenous signals to maximize reproduction. Thus we identified six SOC1-like genes in Gossypium that clustered into two groups. GhSOC1-1 contained a large intron and clustered with monocot SOC1s, while GhSOC1-2/3 were close to dicot SOC1s. GhSOC1s expression gradually increased during seedling development suggesting their conserved function in promoting flowering, which was supported by the early flowering phenotype of 35S:GhSOC1-1 Arabidopsis lines and the delayed flowering of cotton silencing lines. Furthermore, GhSOC1-1 responded to short-day and high temperature conditions, while GhSOC1-2 responded to long-day conditions. GhSOC1-3 might function to promote flowering in response to low temperature and cold. Taken together, our results demonstrate that GhSOC1s respond differently to light and temperature and act cooperatively to activate GhLFY expression to promote floral transition and enlighten us in cotton adaptation to environment that is helpful in improvement of cotton maturity.

Transcriptional Landscape of Cotton Fiber Development and Its Alliance With Fiber-Associated Traits

Cotton fiber development is still an intriguing question to understand fiber commitment and development. At different fiber developmental stages, many genes change their expression pattern and have a pivotal role in fiber quality and yield. Recently, numerous studies have been conducted for transcriptional regulation of fiber, and raw data were deposited to the public repository for comprehensive integrative analysis. Here, we remapped > 380 cotton RNAseq data with uniform mapping strategies that span ∼400 fold coverage to the genome. We identified stage-specific features related to fiber cell commitment, initiation, elongation, and Secondary Cell Wall (SCW) synthesis and their putative cis-regulatory elements for the specific regulation in fiber development. We also mined Exclusively Expressed Transcripts (EETs) that were positively selected during cotton fiber evolution and domestication. Furthermore, the expression of EETs was validated in 100 cotton genotypes through the nCounter assay and correlated with different fiber-related traits. Thus, our data mining study reveals several important features related to cotton fiber development and improvement, which were consolidated in the "CottonExpress-omics" database.

Insights Into MicroRNA-Mediated Regulation of Flowering Time in Cotton Through Small RNA Sequencing

The timing of flowering is a key determinant for plant reproductive. It has been demonstrated that microRNAs (miRNAs) play an important role in transition from the vegetative to reproductive stage in cotton; however, knowledge remains limited about the regulatory role of miRNAs involved in flowering time regulation in cotton. To elucidate the molecular basis of miRNAs in response to flowering time in cotton, we performed high-throughput small RNA sequencing at the fifth true leaf stage. We identified 56 and 43 miRNAs that were significantly up- and downregulated in two elite early flowering cultivars (EFC) compared with two late flowering cultivars (LFC), respectively. The miRNA targets by RNA sequencing analysis showed that GhSPL4 in SBP transcription factor family targeted by GhmiR156 was significantly upregulated in EFCs. Co-expression regulatory network analysis (WGCNA) revealed that GhSOC1, GhAP1, GhFD, GhCOL3, and GhAGL16 act as node genes in the auxin- and gibberellin-mediated flowering time regulatory networks in cotton. Therefore, elucidation of miRNA-mediated flowering time regulatory network will contribute to our understanding of molecular mechanisms underlying flowering time in cotton.

High-Altitude Genetic Selection and Genome-Wide Association Analysis of Yield-Related Traits in Elymus sibiricus L. Using SLAF Sequencing

The genetic adaptations to harsh climatic conditions in high altitudes and genetic basis of important agronomic traits are poorly understood in Elymus sibiricus L. In this study, an association population of 210 genotypes was used for population structure, selective sweep analysis, and genome-wide association study (GWAS) based on 88,506 single nucleotide polymorphisms (SNPs). We found 965 alleles under the natural selection of high altitude, which included 7 hub genes involved in the response to UV, and flavonoid and anthocyanin biosynthetic process based on the protein-protein interaction (PPI) analysis. Using a mixed linear model (MLM), the GWAS test identified a total of 1,825 significant loci associated with 12 agronomic traits. Based on the gene expression data of two wheat cultivars and the PPI analysis, we finally identified 12 hub genes. Especially, in plant height traits, the top hub gene (TOPLESS protein) encoding auxins and jasmonic acid signaling pathway, shoot apical meristem specification, and xylem and phloem pattern formation was highly overexpressed. These genes might play essential roles in controlling the growth and development of E. sibiricus. Therefore, this study provides fundamental insights relevant to hub genes and will benefit molecular breeding and improvement in E. sibiricus and other Elymus species.

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.

Non-functional GoFLA19s are responsible for the male sterility caused by hybrid breakdown in cotton (Gossypium spp.)

Hybrid breakdown (HB) functions as a common reproductive barrier and reduces hybrid fitness in many species, including cotton. However, the related genes and the underlying genetic mechanisms of HB in cotton remain unknown. Here, we found that the photosensitive genetic male sterile line CCRI9106 was a hybrid progeny of Gossypium hirsutum and Gossypium barbadense and probably a product of HB. Fine mapping with F₂ s (CCRI9106 × G. hirsutum/G. barbadense lines) identified a pair of male sterility genes GoFLA19s (encoding fasciclin-like arabinogalactan family protein) located on chromosomes A12 and D12. Crucial variations occurring in the fasciclin-like domain and the arabinogalactan protein domain were predicted to cause the non-functionalization of GbFLA19-D and GhFLA19-A. CRISPR/Cas9-mediated knockout assay confirmed the effects of GhFLA19s on male sterility. Sequence alignment analyses showed that variations in GbFLA19-D and GhFLA19-A likely occurred after the formation of allotetraploid cotton species. GoFLA19s are specifically expressed in anthers and contribute to tapetal development, exine assembly, intine formation, and pollen grain maturation. RNA-sequencing and quantitative reverse transcriptase-polymerase chain reaction analyses illustrated that genes related to these biological processes were significantly downregulated in the mutant. Our research on male sterility genes, GoFLA19s, improves the understanding of the molecular characteristics and evolutionary significance of HB in interspecific hybrid breeding.

Parallel and Intertwining Threads of Domestication in Allopolyploid Cotton

The two cultivated allopolyploid cottons, Gossypium hirsutum and Gossypium barbadense, represent a remarkable example of parallel independent domestication, both involving dramatic morphological transformations under selection from wild perennial plants to annualized row crops. Deep resequencing of 643 newly sampled accessions spanning the wild-to-domesticated continuum of both species, and their allopolyploid relatives, are combined with existing data to resolve species relationships and elucidate multiple aspects of their parallel domestication. It is confirmed that wild G. hirsutum and G. barbadense were initially domesticated in the Yucatan Peninsula and NW South America, respectively, and subsequently spread under domestication over 4000-8000 years to encompass most of the American tropics. A robust phylogenomic analysis of infraspecific relationships in each species is presented, quantify genetic diversity in both, and describe genetic bottlenecks associated with domestication and subsequent diffusion. As these species became sympatric over the last several millennia, pervasive genome-wide bidirectional introgression occurred, often with striking asymmetries involving the two co-resident genomes of these allopolyploids. Diversity scans revealed genomic regions and genes unknowingly targeted during domestication and additional subgenomic asymmetries. These analyses provide a comprehensive depiction of the origin, divergence, and adaptation of cotton, and serve as a rich resource for cotton improvement.

Uncovering Novel Genomic Regions and Candidate Genes for Senescence-Related Traits by Genome-Wide Association Studies in Upland Cotton (Gossypium hirsutum L.)

Senescence in plants is a complex trait, which is controlled by both genetic and environmental factors and can affect the yield and quality of cotton. However, the genetic basis of cotton senescence remains relatively unknown. In this study, we reported genome-wide association studies (GWAS) based on 185 accessions of upland cotton and 26,999 high-quality single-nucleotide polymorphisms (SNPs) to reveal the genetic basis of cotton senescence. To determine cotton senescence, we evaluated eight traits/indices. Our results revealed a high positive correlation (r>0.5) among SPAD value 20 days after topping (SPAD20d), relative difference of SPAD (RSPAD), nodes above white flower on topping day (NAWF0d), nodes above white flower 7 days after topping (NAWF7d), and number of open bolls on the upper four branches (NB), and genetic analysis revealed that all traits had medium or high heritability ranging from 0.53 to 0.86. Based on a multi-locus method (FASTmrMLM), a total of 63 stable and significant quantitative trait nucleotides (QTNs) were detected, which represented 50 genomic regions (GWAS risk loci) associated with cotton senescence. We observed three reliable loci located on chromosomes A02 (A02_105891088_107196428), D03 (D03_37952328_38393621) and D13 (D13_59408561_60730103) because of their high repeatability. One candidate gene (Ghir_D03G011060) was found in the locus D03_37952328_38393621, and its Arabidopsis thaliana homologous gene (AT5G23040) encodes a cell growth defect factor-like protein (CDF1), which might be involved in chlorophyll synthesis and cell death. Moreover, qRT-PCR showed that the transcript level of Ghir_D03G011060 was down-regulated in old cotton leaves, and virus-induced gene silencing (VIGS) indicated that silencing of Ghir_D03G011060 resulted in leaf chlorosis and promoted leaf senescence. In addition, two candidate genes (Ghir_A02G017660 and Ghir_D13G021720) were identified in loci A02_105891088_107196428 and D13_59408561_60730103, respectively. These results provide new insights into the genetic basis of cotton senescence and will serve as an important reference for the development and implementation of strategies to prevent premature senescence in cotton breeding programs.

Genetic Diversity, QTL Mapping, and Marker-Assisted Selection Technology in Cotton (Gossypium spp.)

Cotton genetic resources contain diverse economically important traits that can be used widely in breeding approaches to create of high-yielding elite cultivars with superior fiber quality and adapted to biotic and abiotic stresses. Nevertheless, the creation of new cultivars using conventional breeding methods is limited by the cost and proved to be time consuming process, also requires a space to make field observations and measurements. Decoding genomes of cotton species greatly facilitated generating large-scale high-throughput DNA markers and identification of QTLs that allows confirmation of candidate genes, and use them in marker-assisted selection (MAS)-based breeding programs. With the advances of quantitative trait loci (QTL) mapping and genome-wide-association study approaches, DNA markers associated with valuable traits significantly accelerate breeding processes by replacing the selection with a phenotype to the selection at the DNA or gene level. In this review, we discuss the evolution and genetic diversity of cotton Gossypium genus, molecular markers and their types, genetic mapping and QTL analysis, application, and perspectives of MAS-based approaches in cotton breeding.