Genomic and GWAS analyses demonstrate phylogenomic relationships of Gossypium barbadense in China and selection for fibre length, lint percentage and Fusarium wilt resistance

Deepening genomic sequences of 1081 Gossypium hirsutum accessions reveals novel SNPs and haplotypes relevant for practical breeding utility

Fiber quality is a major breeding goal in cotton, but phenotypically direct selection is often hindered. In this study, we identified fiber quality and yield related loci using GWAS based on 2.97 million SNPs obtained from 10.65× resequencing data of 1081 accessions. The results showed that 585 novel fiber loci, including two novel stable SNP peaks associated with fiber length on chromosomes At12 and Dt05 and one novel genome regions linked with fiber strength on chromosome Dt12 were identified. Furthermore, by means of gene expression analysis, GhM_A12G0090, GhM_D05G1692, GhM_D12G3135 were identified and GhM_D11G2208 function was identified in Arabidopsis. Additionally, 14 consistent and stable superior haplotypes were identified, and 25 accessions were detected as possessing these 14 superior haplotype in breeding. This study providing fundamental insight relevant to identification of genes associated with fiber quality and yield will enhance future efforts toward improvement of upland cotton.

GbPP2C80 Interacts with GbWAKL14 to Negatively Co-Regulate Resistance to Fusarium and Verticillium wilt via MPK3 and ROS Signaling in Sea Island Cotton

Fusarium wilt (FW) is widespread in global cotton production, but the mechanism underlying FW resistance in superior-fiber-quality Sea Island cotton is unclear. This study reveals that FW resistance has been the target of genetic improvement of Sea Island cotton in China since the 2010s. The key nonsynonymous single nucleotide polymorphism (SNP, T/C) of gene Gbar_D03G001670 encoding protein phosphatase 2C 80 (PP2C80) results in an amino acid shift (L/S), which is significantly associated with FW resistance of Sea Island cotton. Silencing GbPP2C80 increases FW resistance in Sea Island cotton, whereas overexpressing GbPP2C80 reduces FW resistance in Arabidopsis. GbPP2C80 and GbWAKL14 exist synergistically in Sea Island cotton accessions with haplotype forms "susceptible-susceptible" (TA) and "resistant-resistant" (CC), and interact with each other. CRISPR/Cas9-mediated knockout of GbWAKL14 enhances FW and Verticillium wilt (VW) resistance in upland cotton and overexpression of GbWAKL14 and GbPP2C80 weakens FW and VW resistance in Arabidopsis. GbPP2C80 and GbWAKL14 respond to FW and VW by modulating reactive oxygen species (ROS) content via affecting MPK3 expression. In summary, two tandem genes on chromosome D03, GbPP2C80, and GbWAKL14, functions as cooperative negative regulators in cotton wilt disease defense, providing novel genetic resources and molecular markers for the development of resistant cotton cultivars.

Genome and GWAS analysis identified genes significantly related to phenotypic state of Rhododendron bark

As an important horticultural plant, Rhododendron is often used in urban greening and landscape design. However, factors such as the high rate of genetic recombination, frequent outcrossing in the wild, weak linkage disequilibrium, and the susceptibility of gene expression to environmental factors limit further exploration of functional genes related to important horticultural traits, and make the breeding of new varieties require a longer time. Therefore, we choose bark as the target trait which is not easily affected by environmental factors, but also has ornamental properties. Genome-wide association study (GWAS) of Rhododendron delavayi (30 samples), R. irroratum (30 samples) and their F1 generation R. agastum (200 samples) was conducted on the roughness of bark phenotypes. Finally, we obtained 2416.31 Gbp of clean data and identified 5 328 800 high-quality SNPs. According to the P-value and the degree of linkage disequilibrium of SNPs, we further identified 4 out of 11 candidate genes that affect bark roughness. The results of gene differential expression analysis further indicated that the expression levels of Rhdel02G0243600 and Rhdel08G0220700 in different bark phenotypes were significantly different. Our study identified functional genes that influence important horticultural traits of Rhododendron, and illustrated the powerful utility and great potential of GWAS in understanding and exploiting wild germplasm genetic resources of Rhododendron.

Dissecting the major genetic components underlying cotton lint development

Numerous genetic loci and several functionally characterized genes have been linked to determination of lint percentage (lint%), one of the most important cotton yield components, but we still know little about the major genetic components underlying lint%. Here, we first linked the genetic loci containing MYB25-like_At and HD1_At to the fiberless seed trait of 'SL1-7-1' and found that MYB25-like_At and HD1_At were very lowly expressed in 'SL1-7-1' ovules during fiber initiation. We then dissected the genetic components involved in determination of lint% using segregating populations derived from crosses of fuzzless mutants and intermediate segregants with different lint%, which not only confirmed the HD1_At locus but identified the HD1_Dt locus as being the major genetic components contributing to fiber initiation and lint%. The segregating populations also allowed us to evaluate the relative contributions of MYB25-like_At, MYB25-like_Dt, HD1_At, and HD1_Dt to lint%. Haplotype analysis of an Upland cotton (Gossypium hirsutum) population with 723 accessions (including 81 fuzzless seed accessions) showed that lint% of the accessions with the LP allele (higher lint%) at MYB25-like_At, MYB25-like_Dt, or HD1_At was significantly higher than that with the lp allele (lower lint%). The lint% of the Upland cotton accessions with 3 or 4 LP alleles at MYB25-like and HD1 was significantly higher than that with 2 LP alleles. The results prompted us to propose a strategy for breeding high-yielding cotton varieties, i.e. pyramiding the LP alleles of MYB25-like and HD1 with new lint% LP alleles without negative impact on seed size and fiber quality.

Genome and haplotype provide insights into the population differentiation and breeding improvement of Gossypium barbadense

INTRODUCTION

Sea-island cotton (Gossypium barbadense, Gb) is one of the major sources of high-grade natural fiber. Besides the common annual Gb cotton, perennial Gb cotton is also cultivated, but studies on perennial Gb cotton are rare.

OBJECTIVES

We aimed to make a systematic analysis of perennial sea-island cotton and lay a foundation for its utilization in breeding, and try to identify the representative structural variations (SVs) in sea-island cotton, and to reveal the population differentiation and adaptive improvement of sea-island cotton.

METHODS

Through genome assembly of one perennial Gb cotton accession (named Gb_M210936) and comparative genome analysis, variations during Gb cotton domestication were identified by comparing Gb_M210936 with annual Gb accession 3-79 and with wild allotetraploid cotton G. darwinii. Six perennial Gb accessions combining with the resequenced 1,129 cotton accessions were used to conduct population and genetic analysis. Large haplotype blocks (haploblocks), generated from interspecific introgressions and intraspecific inversions, were identified and were used to analyze their effects on population differentiation and agronomic traits of sea-island cotton.

RESULTS

One reference genome of perennial sea-island cotton was assembled. Representative SVs in sea-island cotton were identified, and 31 SVs were found to be associated with agronomic traits. Perennial Gb cotton had a closer kinship with the wild-to-landrace continuum Gb cotton from south America where Gb cotton is originally domesticated. Haploblocks were associated with agronomic traits improvement of sea-island cotton, promoted sea-island cotton differentiation into three subgroups, were suffered from breeding selection, and may drive Gb cotton to be adapted to central Asian.

CONCLUSION

Our study made up the lack of perennial Gb cotton genome, and clarified that exotic introgressions improved the traits of sea-island cotton, promoted the population differentiation, and drove sea-island cotton adaptive to central Asia, which will provide new insights for the genetic breeding improvement of sea-island cottons.

Development and validation of allele-specific PCR-based SNP typing in a gene on chromosome D03 conferring resistance to Fusarium wilt race 4 in Upland cotton (Gossypium hirsutum)

Upland cotton (Gossypium hirsutum) is the most important fiber crop for the global textile industry. Fusarium oxysporum f. sp. vasinfectum (FOV) is one of the most destructive soil-borne fungal pathogens in cotton. Among eight pathogenic races and other strains, FOV race 4 (FOV4) is the most virulent race in US cotton production. A single nucleotide polymorphism (SNP) in a glutamate receptor-like gene (GhGLR4.8) on chromosome D03 was previously identified and validated to confer resistance to FOV race 7, and targeted genome sequencing demonstrated that it was also associated with resistance to FOV4. The objective of this study was to develop an easy and convenient PCR-based marker assay. To target the resistance SNP, a forward primer for the SNP with a mismatch in the 3rd position was designed for both the resistance (R) and susceptibility (S) alleles, respectively, with addition of 20-mer T7 promoter primer to the 5' end of the forward primer for the R allele. The two forward primers, in combination with each of five common reverse primers, were targeted to amplify amplicons of 50-260 bp in size with R and S alleles differing in 20 bp. Results showed that each of three common reverse primers in combination with the two forward primers produced polymorphic markers between R and S plants that were consistent with the targeted genome sequencing results. The polymorphism was distinctly resolved using both polyacrylamide and agarose gel electrophoreses. In addition, a sequence comparative analysis between the resistance gene and homologous sequences in sequenced tetraploid and diploid A and D genome species showed that none of the species possessed the resistance gene allele, suggesting its recent origin from a natural point mutation. The allele-specific PCR-based SNP typing method based on a three-primer combination provides a fast and convenient marker-assisted selection method to search and select for FOV4-resistant Upland cotton.

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.

UAV-based time-series phenotyping reveals the genetic basis of plant height in upland cotton

Plant height (PH) is an important agronomic trait affecting crop architecture, biomass, resistance to lodging and mechanical harvesting. Elucidating the genetic governance of plant height is crucial because of the global demand for high crop yields. However, during the rapid growth period of plants the PH changes a lot on a daily basis, which makes it difficult to accurately phenotype the trait by hand on a large scale. In this study, an unmanned aerial vehicle (UAV)-based remote-sensing phenotyping platform was applied to obtain time-series PHs of 320 upland cotton accessions in three different field trials. The results showed that the PHs obtained from UAV images were significantly correlated with ground-based manual measurements, for three trials (R2  = 0.96, 0.95 and 0.96). Two genetic loci on chromosomes A01 and A11 associated with PH were identified by genome-wide association studies (GWAS). GhUBP15 and GhCUL1 were identified to influence PH in further analysis. We obtained a time series of PH values for three field conditions based on remote sensing with UAV. The key genes identified in this study are of great value for the breeding of ideal plant architecture in cotton.

Dissection of the Genetic Basis of Resistance to Stem Rot in Cultivated Peanuts (Arachis hypogaea L.) through Genome-Wide Association Study

Peanut (Arachis hypogaea) is an important oilseed and cash crop worldwide, contributing an important source of edible oil and protein for human nutrition. However, the incidence of stem rot disease caused by Athelia rolfsii poses a major challenge to peanut cultivation, resulting in significant yield losses. In this study, a panel of 202 peanut accessions was evaluated for their resistance to stem rot by inoculating plants in the field with A. rolfsii-infested oat grains in three environments. The mean disease index value of each environment for accessions in subsp. fasitigiate and subsp. hypogaea showed no significant difference. Accessions from southern China displayed the lowest disease index value compared to those from other ecological regions. We used whole-genome resequencing to analyze the genotypes of the accessions and to identify significant SNPs associated with stem rot resistance through genome-wide association study (GWAS). A total of 121 significant SNPs associated with stem rot resistance in peanut were identified, with phenotypic variation explained (PVE) ranging from 12.23% to 15.51%. A total of 27 candidate genes within 100 kb upstream and downstream of 23 significant SNPs were annotated, which have functions related to recognition, signal transduction, and defense response. These significant SNPs and candidate genes provide valuable information for further validation and molecular breeding to improve stem rot resistance in peanut.

Natural variation in Beauty Mark is associated with UV-based geographical adaptation in Gossypium species

BACKGROUND

Anthocyanins, a class of specialized metabolites that are ubiquitous among plant species, have attracted a great deal of attention from plant biologists due to their chemical diversity. They confer purple, pink, and blue colors that attract pollinators, protect plants from ultraviolet (UV) radiation, and scavenge reactive oxygen species (ROS) to facilitate plant survival during abiotic stress. In a previous study, we identified Beauty Mark (BM) in Gossypium barbadense as an activator of the anthocyanin biosynthesis pathway; this gene also directly led to the formation of a pollinator-attracting purple spot.

RESULTS

Here, we found that a single nucleotide polymorphism (SNP) (C/T) within the BM coding sequence was responsible for variations in this trait. Transient expression assays of BM from G. barbadense and G. hirsutum in Nicotiana benthamiana using luciferase reporter gene also suggested that SNPs in the coding sequence could be responsible for the absent beauty mark phenotype observed in G. hirsutum. We next demonstrated that the beauty mark and UV floral patterns are associated phenotypes and that UV exposure resulted in increased ROS generation in floral tissues; BM thus contributed to ROS scavenging in G. barbadense and wild cotton plants with flowers containing the beauty mark. Furthermore, a nucleotide diversity analysis and Tajima's D Test suggested that there have been strong selective sweeps in the GhBM locus during G. hirsutum domestication.

CONCLUSIONS

Taken together, these results suggest that cotton species differ in their approaches to absorbing or reflecting UV light and thus exhibit variations in floral anthocyanin biosynthesis to scavenge reactive ROS; furthermore, these traits are related to the geographic distribution of cotton species.

Resequencing of Rosa rugosa accessions revealed the history of population dynamics, breed origin, and domestication pathways

BACKGROUND

Rosa rugosa is a shrub that originated in China and has economic and ecological value. However, during the development of R. rugosa, the genetic background was chaotic, and the genetic structure among different wild populations was unclear, as well as wild and cultivated accessions. Here, we report whole-genome resequencing of wild and cultivated R. rugosa accessions.

RESULTS

A total of 19,041,284 SNPs were identified in 188 R. rugosa accessions and 3 R. chinensis accessions by resequencing. Population genetic analysis revealed that cultivated and wild groups were separated very early. All R. rugosa accessions were divided into 8 categories based on genetic structure: (1) Weihai, Yantai, and Liaoning category, (2) Jilin category, and (3) Hammonasset category (above three are wild); (4) traditional varieties, (5) hybrids between R. rugosa and R. chinensis, (6) Zizhi Rose, (7) Kushui Rose, (8) hybrids between R. rugosa and R. multiflora. We found that the heterozygosity and genetic diversity of wild accessions were generally lower than those of cultivated individuals. The genes that were selected during cultivation were identified, and it was found that these genes were mainly related to environmental adaptation and growth.

CONCLUSIONS

The Jilin population was the oldest population and later migrated to Liaoning and then migrated to Yantai and Weihai by sea regression in the Bohai Basin. The Hammonasset naturalized population probably originated from the Jilin population and then experienced separate differentiation. The long-term asexual reproduction pattern of R. rugosa decreased genetic diversity in the wild population. During R. rugosa cultivation, the ancestors of the Jilin population were involved in breeding traditional varieties, after which almost no wild individuals were engaged in breeding. However, in recent decades, cross breeding of R. rugosa started the utilization of wild germplasms. In comparison, some other species play important roles in variety formation. Few genes related to economic traits were selected, suggesting no directional domestication in the R. rugosa cultivation process.

Targeted development of diagnostic SNP markers for resistance to Fusarium wilt race 4 in Upland cotton (Gossypium hirsutum)

Fusarium wilt caused by the soil-borne fungus Fusarium oxysporum f. sp. vasinfectum (FOV) race 4 (FOV4) has become one of the most important emerging diseases in US cotton production. Numerous QTLs have been reported for resistance to FOV; however, no major FOV4-resistance QTL or gene has been identified and used in breeding Upland cotton (Gossypium hirsutum) for FOV4 resistance. In this study, a panel of 223 Chinese Upland cotton accessions was evaluated for FOV4 resistance based on seedling mortality rate (MR) and stem and root vascular discoloration (SVD and RVD). SNP markers were developed based on targeted genome sequencing using AgriPlex Genomics. The chromosome region at 2.130-2.292 Mb on D03 was significantly correlated with both SVD and RVD but not with MR. Based on the two most significant SNP markers, accessions homozygous for AA or TT SNP genotype averaged significantly lower SVD (0.88 vs. 2.54) and RVD (1.46 vs. 3.02) than those homozygous for CC or GG SNP genotype. The results suggested that a gene or genes within the region conferred resistance to vascular discoloration caused by FOV4. The Chinese Upland accessions had 37.22% homozygous AA or TT SNP genotype and 11.66% heterozygous AC or TG SNP genotype, while 32 US elite public breeding lines all had the CC or GG SNP genotype. Among 463 obsolete US Upland accessions, only 0.86% possessed the AA or TT SNP genotype. This study, for the first time, has developed diagnostic SNPs for marker-assisted selection and identified FOV4-resistant Upland germplasms with the SNPs.

Structural variation (SV)-based pan-genome and GWAS reveal the impacts of SVs on the speciation and diversification of allotetraploid cottons

Structural variations (SVs) have long been described as being involved in the origin, adaption, and domestication of species. However, the underlying genetic and genomic mechanisms are poorly understood. Here, we report a high-quality genome assembly of Gossypium barbadense acc. Tanguis, a landrace that is closely related to formation of extra-long-staple (ELS) cultivated cotton. An SV-based pan-genome (Pan-SV) was then constructed using a total of 182 593 non-redundant SVs, including 2236 inversions, 97 398 insertions, and 82 959 deletions from 11 assembled genomes of allopolyploid cotton. The utility of this Pan-SV was then demonstrated through population structure analysis and genome-wide association studies (GWASs). Using segregation mapping populations produced through crossing ELS cotton and the landrace along with an SV-based GWAS, certain SVs responsible for speciation, domestication, and improvement in tetraploid cottons were identified. Importantly, some of the SVs presently identified as associated with the yield and fiber quality improvement had not been identified in previous SNP-based GWAS. In particular, a 9-bp insertion or deletion was found to associate with elimination of the interspecific reproductive isolation between Gossypium hirsutum and G. barbadense. Collectively, this study provides new insights into genome-wide, gene-scale SVs linked to important agronomic traits in a major crop species and highlights the importance of SVs during the speciation, domestication, and improvement of cultivated crop species.

Fine mapping and candidate gene analysis of qFL-A12-5: a fiber length-related QTL introgressed from Gossypium barbadense into Gossypium hirsutum

The fiber length-related qFL-A12-5 identified in CSSLs introgressed from Gossypium barbadense into Gossypium hirsutum was fine-mapped to an 18.8 kb region on chromosome A12, leading to the identification of the GhTPR gene as a potential regulator of cotton fiber length. Fiber length is a key determinant of fiber quality in cotton, and it is a key target of artificial selection for breeding and domestication. Although many fiber length-related quantitative trait loci have been identified, there are few reports on their fine mapping or candidate gene validation, thus hampering efforts to understand the mechanistic basis of cotton fiber development. Our previous study identified the qFL-A12-5 associated with superior fiber quality on chromosome A12 in the chromosome segment substitution line (CSSL) MBI7747 (BC₄F_3:5). A single segment substitution line (CSSL-106) screened from BC₆F₂ was backcrossed to construct a larger segregation population with its recurrent parent CCRI45, thus enabling the fine mapping of 2852 BC₇F₂ individuals using denser simple sequence repeat markers to narrow the qFL-A12-5 to an 18.8 kb region of the genome, in which six annotated genes were identified in Gossypium hirsutum. Quantitative real-time PCR and comparative analyses led to the identification of GH_A12G2192 (GhTPR) encoding a tetratricopeptide repeat-like superfamily protein as a promising candidate gene for qFL-A12-5. A comparative analysis of the protein-coding regions of GhTPR among Hai1, MBI7747, and CCRI45 revealed two non-synonymous mutations. The overexpression of GhTPR resulted in longer roots in Arabidopsis, suggesting that GhTPR may regulate cotton fiber development. These results provide a foundation for future efforts to improve cotton fiber length.

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.

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.

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 analysis of serine carboxypeptidase-like protein (SCPL) family and functional validation of Gh_SCPL42 unchromosome conferring cotton Verticillium der Verticillium wilt stress in Gossypium hirsutum

BACKGROUND

Serine carboxypeptidase-like protein (SCPL) plays an important role in response to stress in plant. However, our knowledge of the function of the SCPL gene family is limited.

RESULTS

In this study, a comprehensive and systematic analysis of SCPL gene family was conducted to explore the phylogeny and evolution of the SCPL gene in Gossypium hirsutum. The phenotype and molecular mechanism of silencing of the Gh_SCPL42 under Verticillium wilt stress was also studied. Our results showed that 96 SCPL genes were observed in genome of G. hirsutum, which distributed on 25 chromosomes and most of them were located in the nucleus. The phylogenetic tree analysis showed that members of SCPL gene family can be divided into three subgroups in G. hirsutum, which are relatively conservative in evolution. SCPL gene has a wide range of tissue expression types in G. hirsutum. Promoter analysis showed that the most cis-acting elements related to MeJA and ABA were contained. Through RNA-seq combined with genotyping, it was found that 11 GhSCPL genes not only had significant expression changes during Verticillium wilt stress but also had differential SNPs in the upstream, downstream, exonic or intronic regions. The expression of these 11 genes in the resistant (Zhongzhimian 2) and susceptible (Junmian 1) materials was further analyzed by qRT-PCR, it was found that 6 genes showed significant expression differences in the two materials. Among them, Gh_SCPL42 has the most obvious expression change. Furthermore, virus-induced gene silencing (VIGS) showed necrosis and yellowing of leaves and significantly higher disease severity index (DSI) and disease severity rate (DSR) values in VIGS plants than in control silenced Gh_SCPL42 plants. Moreover, the expression levels of genes related to the SA and JA pathways were significantly downregulated. These results show that Gh_SCPL42 might improve resistance to Verticillium wilt through the SA and JA pathways in G. hirsutum.

CONCLUSION

In conclusion, our findings indicated that Gh_SCPL42 gene plays an important role in resistance to Verticillium wilt in cotton. It was provided an important theoretical basis for further research on the function of SCPL gene family and the molecular mechanism of resistance to Verticillium wilt in cotton.

Gb_ANR-47 Enhances the Resistance of Gossypium barbadense to Fusarium oxysporum f. sp. vasinfectum (FOV) by Regulating the Content of Proanthocyanidins

Anthocyanidin reductase (ANR) is an important regulator of flavonoid metabolism, and proanthocyanidins, the secondary metabolites of flavonoids, play an important role in the response of plants to pathogenic stress. Therefore, in this study, the expression analysis of the ANR gene family of Gossypium barbadense after inoculation with Fusarium oxysporum f. sp. vasinfectum (FOV) was performed at different time points. It was found that Gb_ANR-47 showed significant differences in the disease-resistant cultivar 06-146 and the susceptible cultivar Xinhai 14, as well as in the highest root expression. It was found that the expression of Gb_ANR-47 in the resistant cultivar was significantly higher than that in the susceptible cultivar by MeJA and SA, and different amounts of methyl jasmonate (MeJA) and salicylic acid (SA) response elements were found in the promoter region of Gb_ANR-47. After silencing GbANR-47 in 06-146 material by VIGS technology, its resistance to FOV decreased significantly. The disease severity index (DSI) was significantly increased, and the anthocyanin content was significantly decreased in silenced plants, compared to controls. Our findings suggest that GbANR-47 is a positive regulator of FOV resistance in Gossypium barbadense. The research results provide an important theoretical basis for in-depth analysis of the molecular mechanism of GbANR-47 and improving the anti-FOV of Gossypium barbadense.

A GWAS identified a major QTL for resistance to Fusarium wilt (Fusarium oxysporum f. sp. vasinfectum) race 4 in a MAGIC population of Upland cotton and a meta-analysis of QTLs for Fusarium wilt resistance

A major QTL conferring resistance to Fusarium wilt race 4 in a narrow region of chromosome D02 was identified in a MAGIC population of 550 RILs of Upland cotton. Numerous studies have been conducted to investigate the genetic basis of Fusarium wilt (FW, caused by Fusarium oxysporum f. sp. vasinfectum, FOV) resistance using bi-parental and association mapping populations in cotton. In this study, a multi-parent advanced generation inter-cross (MAGIC) population of 550 recombinant inbred lines (RILs), together with their 11 Upland cotton (Gossypium hirsutum) parents, was used to identify QTLs for FOV race 4 (FOV4) resistance. Among the parents, Acala Ultima, M-240 RNR, and Stoneville 474 were the most resistant, while Deltapine Acala 90, Coker 315, and Stoneville 825 were the most susceptible. Twenty-two MAGIC lines were consistently resistant to FOV4. Through a genome-wide association study (GWAS) based on 473,516 polymorphic SNPs, a major FOV4 resistance QTL within a narrow region on chromosomes D02 was detected, allowing identification of 14 candidate genes. Additionally, a meta-analysis of 133 published FW resistance QTLs showed a D subgenome and individual chromosome bias and no correlation between homeologous chromosome pairs. This study represents the first GWAS study using a largest genetic population and the most comprehensive meta-analysis for FW resistance in cotton. The results illustrated that 550 lines were not enough for high resolution mapping to pinpoint a candidate gene, and experimental errors in phenotyping cotton for FW resistance further compromised the accuracy and precision in QTL localization and identification of candidate genes. This study identified FOV4-resistant parents and MAGIC lines, and the first major QTL for FOV4 resistance in Upland cotton, providing useful information for developing FOV4-resistant cultivars and further genomic studies towards identification of causal genes for FOV4 resistance in cotton.

Identification of Candidate Genes for Lint Percentage and Fiber Quality Through QTL Mapping and Transcriptome Analysis in an Allotetraploid Interspecific Cotton CSSLs Population

Upland cotton (Gossypium hirsutum) has long been an important fiber crop, but the narrow genetic diversity of modern G. hirsutum limits the potential for simultaneous improvement of yield and fiber quality. It is an effective approach to broaden the genetic base of G. hirsutum through introgression of novel alleles from G. barbadense with excellent fiber quality. In the present study, an interspecific chromosome segment substitution lines (CSSLs) population was established using G. barbadense cultivar Pima S-7 as the donor parent and G. hirsutum cultivar CCRI35 as the recipient parent. A total of 105 quantitative trait loci (QTL), including 85 QTL for fiber quality and 20 QTL for lint percentage (LP), were identified based on phenotypic data collected from four environments. Among these QTL, 25 stable QTL were detected in two or more environments, including four for LP, eleven for fiber length (FL), three for fiber strength (FS), six for fiber micronaire (FM), and one for fiber elongation (FE). Eleven QTL clusters were observed on nine chromosomes, of which seven QTL clusters harbored stable QTL. Moreover, eleven major QTL for fiber quality were verified through analysis of introgressed segments of the eight superior lines with the best comprehensive phenotypes. A total of 586 putative candidate genes were identified for 25 stable QTL associated with lint percentage and fiber quality through transcriptome analysis. Furthermore, three candidate genes for FL, GH_A08G1681 (GhSCPL40), GH_A12G2328 (GhPBL19), and GH_D02G0370 (GhHSP22.7), and one candidate gene for FM, GH_D05G1346 (GhAPG), were identified through RNA-Seq and qRT-PCR analysis. These results lay the foundation for understanding the molecular regulatory mechanism of fiber development and provide valuable information for marker-assisted selection (MAS) in cotton breeding.

A Calmodulin-Like Gene (GbCML7) for Fiber Strength and Yield Improvement Identified by Resequencing Core Accessions of a Pedigree in Gossypium barbadense

Sea Island cotton (Gossypium barbadense) is world-renowned for its superior natural fiber. Although fiber strength is one of the most important fiber quality traits, genes contributing to fiber strength are poorly understood. Production of sea island cotton also is inextricably linked to improving its relatively low yield, thus enhancing the importance of joint improvement of both fiber quality and yield. We used genomic variation to uncover the genetic evidence of trait improvement resulting from pedigree breeding of Sea Island cotton. This pedigree was aimed at improving fiber strength and yielded an elite cultivar, XH35. Using a combination of genome-wide association study (GWAS) and selection screens, we detected 82 putative fiber-strength-related genes. Expression analysis confirmed a calmodulin-like gene, GbCML7, which enhanced fiber strength in a specific haplotype. This gene is a major-effect gene, which interacts with a minor-effect gene, GbTUA3, facilitating the enhancement of fiber strength in a synergistic fashion. Moreover, GbCML7 participates in the cooperative improvement of fiber strength, fiber length, and fiber uniformity, though a slight compromise exists between the first two of these traits and the latter. Importantly, GbCML7 is shown to boost yield in some backgrounds by increasing multiple yield components to varying degrees, especially boll number. Our work provides valuable genomic evidence and a key genetic factor for the joint improvement of fiber quality and yield in Sea Island cotton.

Genome-editing in millets: current knowledge and future perspectives

Millets are small seeded cereal crops predominantly cultivated and consumed by resource-poor farmers in the semi-arid tropics of Asia and Africa. Millets possess rich nutrients and a climate resilience property when compared to the other cereals such as rice and wheat. Millet improvement using modern genetic and genomic tools is falling behind other cereal crops due to their cultivation being restricted to less developed countries. Genome editing tools have been successfully applied to major cereal crops and, as a result, many key traits have been introduced into rice, wheat and maize. However, genome editing tools have not yet been used for most millets although they possess rich nutrients. The foxtail millet is the only millet utilised up to now for genome editing works. Limited genomic resources and lack of efficient transformation systems may slow down genome editing in millets. As millets possess many important traits of agricultural importance, high resolution studies with genome editing tools will help to understand the specific mechanism and transfer such traits to major cereals in the future. This review covers the current status of genome editing studies in millets and discusses the future prospects of genome editing in millets to understand key traits of nutrient fortification and develop climate resilient crops in the future.