Subgenome-aware analyses reveal the genomic consequences of ancient allopolyploid hybridizations throughout the cotton family.

Malvaceae comprise some 4,225 species in 243 genera and nine subfamilies and include economically important species, such as cacao, cotton, durian, and jute, with cotton an important model system for studying the domestication of polyploids. Here, we use chromosome-level genome assemblies from representatives of five or six subfamilies (depending on the placement of Ochroma) to differentiate coexisting subgenomes and their evolution during the family's deep history. The results reveal that the allohexaploid Helicteroideae partially derive from an allotetraploid Sterculioideae and also form a component of the allodecaploid Bombacoideae and Malvoideae. The ancestral Malvaceae karyotype consists of 11 protochromosomes. Four subfamilies share a unique reciprocal chromosome translocation, and two other subfamilies share a chromosome fusion. DNA alignments of single-copy nuclear genes do not yield the same relationships as inferred from chromosome structural traits, probably because of genes originating from different ancestral subgenomes. These results illustrate how chromosome-structural data can unravel the evolutionary history of groups with ancient hybrid genomes.

Haplotype-resolved assembly of auto-polyploid genomes via combining Hi-C and gametic data.

Haplotype-resolved genome assembly plays a crucial role in understanding allele-specific functions. However, obtaining haplotype-resolved assembly for auto-polyploid genomes remains challenging. Existing methods can be classified into reference-based phasing, assembly-based phasing, and gamete binning. Nevertheless, there is a lack of cost-effective and efficient methods for haplotyping auto-polyploid genomes. In this study, we propose a novel phasing algorithm called PolyGH, which combines Hi-C and gametic data. We conducted experiments on tetraploid potato cultivars and divided the method into three steps. Firstly, gametic data was utilized to bin non-collapsed contigs, followed by merging adjacent fragments of the same type within the same contig. Secondly, accurate Hi-C signals related to differential genomic regions were acquired using unique k-mers. Finally, collapsed fragments were assigned to haplotigs based on combined Hi-C and gametic signals. Comparing PolyGH with Hi-C-based and gametic data-based methods, we found that PolyGH exhibited superior performance in haplotyping auto-polyploid genomes when integrating both data types. This approach has the potential to enhance haplotype-resolved assembly for auto-polyploid genomes.

Genome assemblies of 11 bamboo species highlight diversification induced by dynamic subgenome dominance.

Polyploidy (genome duplication) is a pivotal force in evolution. However, the interactions between parental genomes in a polyploid nucleus, frequently involving subgenome dominance, are poorly understood. Here we showcase analyses of a bamboo system (Poaceae: Bambusoideae) comprising a series of lineages from diploid (herbaceous) to tetraploid and hexaploid (woody), with 11 chromosome-level de novo genome assemblies and 476 transcriptome samples. We find that woody bamboo subgenomes exhibit stunning karyotype stability, with parallel subgenome dominance in the two tetraploid clades and a gradual shift of dominance in the hexaploid clade. Allopolyploidization and subgenome dominance have shaped the evolution of tree-like lignified culms, rapid growth and synchronous flowering characteristic of woody bamboos as large grasses. Our work provides insights into genome dominance in a remarkable polyploid system, including its dependence on genomic context and its ability to switch which subgenomes are dominant over evolutionary time.

The genetics of cardiomyocyte polyploidy.

The regulation of ploidy in cardiomyocytes is a complex and tightly regulated aspect of cardiac development and function. Cardiomyocyte ploidy can range from diploid (2N) to 8N or even 16N, and these states change during key stages of development and disease progression. Polyploidization has been associated with cellular hypertrophy to support normal growth of the heart, increased contractile capacity, and improved stress tolerance in the heart. Conversely, alterations to ploidy also occur during cardiac pathogenesis of diseases, such as ischemic and non-ischemic heart failure and arrhythmia. Therefore, understanding which genes control and modulate cardiomyocyte ploidy may provide mechanistic insight underlying cardiac growth, regeneration, and disease. This chapter summarizes the current knowledge regarding the genes involved in the regulation of cardiomyocyte ploidy. We discuss genes that have been directly tested for their role in cardiomyocyte polyploidization, as well as methodologies used to identify ploidy alterations. These genes encode cell cycle regulators, transcription factors, metabolic proteins, nuclear scaffolding, and components of the sarcomere, among others. The general physiological and pathological phenotypes in the heart associated with the genetic manipulations described, and how they coincide with the respective cardiomyocyte ploidy alterations, are further discussed in this chapter. In addition to being candidates for genetic-based therapies for various cardiac maladies, these genes and their functions provide insightful evidence regarding the purpose of widespread polyploidization in cardiomyocytes.

Polyploidy and mTOR signaling: a possible molecular link.

Polyploidy is typically described as the condition wherein a cell or organism has more than two complete sets of chromosomes. Occurrence of polyploidy is a naturally occurring phenomenon in the body's development and differentiation processes under normal physiological conditions. However, in pathological conditions, the occurrence of polyploidy is documented in numerous disorders, including cancer, aging and diabetes. Due to the frequent association that the polyploidy has with these pathologies and physiological process, understanding the cause and consequences of polyploidy would be beneficial to develop potential therapeutic applications. Many of the genetic and epigenetic alterations leading to cancer, diabetes and aging are linked to signaling pathways. Nonetheless, the specific signaling pathway associated with the cause and consequences of polyploidy still remains largely unknown. Mammalian/mechanistic target of rapamycin (mTOR) plays a key role in the coordination between eukaryotic cell growth and metabolism, thereby simultaneously respond to various environmental inputs including nutrients and growth factors. Extensive research over the past two decades has established a central role for mTOR in the regulation of many fundamental cellular processes that range from protein synthesis to autophagy. Dysregulated mTOR signaling has been found to be implicated in various disease progressions. Importantly, there is a strong correlation between the hallmarks of polyploidy and dysregulated mTOR signaling. In this review, we explore and discuss the molecular connection between mTOR signaling and polyploidy along with its association with cancer, diabetes and aging. Additionally, we address some unanswered questions and provide recommendations to further advance our understanding of the intricate relationship between mTOR signaling and polyploidy.

Autotetraploid Induction of Three A-Genome Wild Peanut Species, Arachis cardenasii, A. correntina, and A. diogoi.

A-genome Arachis species (AA; 2n = 2x = 20) are commonly used as secondary germplasm sources in cultivated peanut breeding, Arachis hypogaea L. (AABB; 2n = 4x = 40), for the introgression of various biotic and abiotic stress resistance genes. Genome doubling is critical to overcoming the hybridization barrier of infertility that arises from ploidy-level differences between wild germplasm and cultivated peanuts. To develop improved genome doubling methods, four trials of various concentrations of the mitotic inhibitor treatments colchicine, oryzalin, and trifluralin were tested on the seedlings and seeds of three A-genome species, A. cardenasii, A. correntina, and A. diogoi. A total of 494 seeds/seedlings were treated in the present four trials, with trials 1 to 3 including different concentrations of the three chemical treatments on seedlings, and trial 4 focusing on the treatment period of 5 mM colchicine solution treatment of seeds. A small number of tetraploids were produced from the colchicine and oryzalin gel treatments of seedlings, but all these tetraploid seedlings reverted to diploid or mixoploid states within six months of treatment. In contrast, the 6-h colchicine solution treatment of seeds showed the highest tetraploid conversion rate (6-13% of total treated seeds or 25-40% of surviving seedlings), and the tetraploid plants were repeatedly tested as stable tetraploids. In addition, visibly and statistically larger leaves and flowers were produced by the tetraploid versions of these three species compared to their diploid versions. As a result, stable tetraploid plants of each A-genome species were produced, and a 5 mM colchicine seed treatment is recommended for A-genome and related wild Arachis species genome doubling.

Physiological and Molecular Modulations to Drought Stress in the Brassica Species.

Climate change, particularly drought stress, significantly impacts plant growth and development, necessitating the development of resilient crops. This study investigated physiological and molecular modulations to drought stress between diploid parent species and their polyploid progeny in the Brassica species. While no significant phenotypic differences were observed among the six species, drought stress reduced growth parameters by 2.4% and increased oxidative stress markers by 1.4-fold. Drought also triggered the expression of genes related to stress responses and led to the accumulation of specific metabolites. We also conducted the first study of perfluorooctane sulfonic acid (PFOS) levels in leaves as a drought indicator. Lower levels of PFOS accumulation were linked to plants taking in less water under drought conditions. Both diploid and polyploid species responded to drought stress similarly, but there was a wide range of variation in their responses. In particular, responses were less variable in polyploid species than in diploid species. This suggests that their additional genomic components acquired through polyploidy may improve their flexibility to modulate stress responses. Despite the hybrid vigor common in polyploid species, Brassica polyploids demonstrated intermediate responses to drought stress. Overall, this study lays the framework for future omics-level research, including transcriptome and proteomic studies, to deepen our understanding of drought tolerance mechanisms in Brassica species.

Copy number variation sequencing for the products of conception: What is the optimal testing strategy.

BACKGROUND: Copy number variation sequencing (CNV-seq) is crucial in prenatal diagnosis, but its limitations in detecting polyploidy, maternal cell contamination (MCC), and uniparental disomy (UPD) restrict its application in the analysis of products of conception (POCs). This study aimed to investigate an optimal genetic testing strategy for POCs in the era of CNV-seq. METHODS: CNV-seq and quantitative fluorescent polymerase chain reaction (QF-PCR) were performed in all 4,211 spontaneous miscarriage cases. Different testing strategies were compared and the optimal testing strategies were proposed. RESULTS: Of the 4,211 cases, 2561 (60.82%) exhibited clinically significant chromosomal abnormalities. CNV-seq alone, without QF-PCR, might misdiagnose 311 (7.39%) cases, including 278 polyploidy, 13 UPD, and 20 MCC. In 20 MCC cases identified by QF-PCR, CNV-seq successfully pinpointed the cause of miscarriage in 13 cases. Furthermore, in cases where QF-PCR suggested polyploidy, CNV-seq improved the diagnostic accuracy in 54 (1.28%) hypo/hypertriploidy cases. After comparing four different strategies, the sequential approach (initiating with CNV-seq followed by QF-PCR if necessary) emerged as advantageous, reducing approximately 70% of the cost associated with QF-PCR while maintaining result accuracy. CONCLUSIONS: We propose an initial CNV-seq followed by QF-PCR if needed-an efficient and cost-effective strategy for the genetic analysis of POCs.

Genome-wide DNA methylation dynamics following recent polyploidy in the allotetraploid Tragopogon miscellus (Asteraceae).

Polyploidy is an important evolutionary force, yet epigenetic mechanisms, such as DNA methylation, that regulate genome-wide expression of duplicated genes remain largely unknown. Here, we use Tragopogon (Asteraceae) as a model system to discover patterns and temporal dynamics of DNA methylation in recently formed polyploids. The naturally occurring allotetraploid Tragopogon miscellus formed in the last 95-100 yr from parental diploids Tragopogon dubius and T. pratensis. We profiled the DNA methylomes of these three species using whole-genome bisulfite sequencing. Genome-wide methylation levels in T. miscellus were intermediate between its diploid parents. However, nonadditive CG and CHG methylation occurred in transposable elements (TEs), with variation among TE types. Most differentially methylated regions (DMRs) showed parental legacy, but some novel DMRs were detected in the polyploid. Differentially methylated genes (DMGs) were also identified and characterized. This study provides the first assessment of both overall and locus-specific patterns of DNA methylation in a recent natural allopolyploid and shows that novel methylation variants can be generated rapidly after polyploid formation. Together, these results demonstrate that mechanisms to regulate duplicate gene expression may arise soon after allopolyploid formation and that these mechanisms vary among genes.

Dynamics of accessible chromatin regions and subgenome dominance in octoploid strawberry.

Subgenome dominance has been reported in diverse allopolyploid species, where genes from one subgenome are preferentially retained and are more highly expressed than those from other subgenome(s). However, the molecular mechanisms responsible for subgenome dominance remain poorly understood. Here, we develop genome-wide map of accessible chromatin regions (ACRs) in cultivated strawberry (2n = 8x = 56, with A, B, C, D subgenomes). Each ACR is identified as an MNase hypersensitive site (MHS). We discover that the dominant subgenome A contains a greater number of total MHSs and MHS per gene than the submissive B/C/D subgenomes. Subgenome A suffers fewer losses of MHS-related DNA sequences and fewer MHS fragmentations caused by insertions of transposable elements. We also discover that genes and MHSs related to stress response have been preferentially retained in subgenome A. We conclude that preservation of genes and their cognate ACRs, especially those related to stress responses, play a major role in the establishment of subgenome dominance in octoploid strawberry.

Polyploid goldback and silverback ferns (Pentagramma) occupy a wider, colder, and wetter bioclimatic niche than diploid counterparts.

PREMISE: The western North American fern genus Pentagramma (Pteridaceae) is characterized by complex patterns of ploidy variation, an understanding of which is critical to comprehending both the evolutionary processes within the genus and its current diversity. METHODS: We undertook a cytogeographic study across the range of the genus, using a combination of chromosome counts and flow cytometry to infer ploidy level. Bioclimatic variables and elevation were used to compare niches. RESULTS: We found that diploids and tetraploids are common and widespread, and triploids are rare and sporadic; in contrast with genome size inferences in earlier studies, no hexaploids were found. Diploids and tetraploids show different geographic ranges: only tetraploids were found in the northernmost portion of the range (Washington, Oregon, and British Columbia) and only diploids were found in the Sierra Nevada of California. Diploid, triploid, and tetraploid cytotypes were found to co-occur in relatively few localities: in the southern (San Diego County, California) and desert Southwest (Arizona) parts of the range, and along the Pacific Coast of California. CONCLUSIONS: Tetraploids occupy a wider bioclimatic niche than diploids both within P. triangularis and at the genus-wide scale. It is unknown whether the wider niche of tetraploids is due to their expansion upon the diploid niche, if diploids have contracted their niche due to competition or changing abiotic conditions, or if this wider niche occupancy is due to multiple origins of tetraploids.

polyGBLUP: a modified genomic best linear unbiased prediction improved the genomic prediction efficiency for autopolyploid species.

Given the universality of autopolyploid species in nature, it is crucial to develop genomic selection methods that consider different allele dosages for autopolyploid breeding. However, no method has been developed to deal with autopolyploid data regardless of the ploidy level. In this study, we developed a modified genomic best linear unbiased prediction (GBLUP) model (polyGBLUP) through constructing additive and dominant genomic relationship matrices based on different allele dosages. polyGBLUP could carry out genomic prediction for autopolyploid species regardless of the ploidy level. Through comprehensive simulations and analysis of real data of autotetraploid blueberry and guinea grass and autohexaploid sweet potato, the results showed that polyGBLUP achieved higher prediction accuracy than GBLUP and its superiority was more obvious when the ploidy level of autopolyploids is high. Furthermore, when the dominant effect was added to polyGBLUP (polyGDBLUP), the greater the dominance degree, the more obvious the advantages of polyGDBLUP over the diploid models in terms of prediction accuracy, bias, mean squared error and mean absolute error. For real data, the superiority of polyGBLUP over GBLUP appeared in blueberry and sweet potato populations and a part of the traits in guinea grass population due to the high correlation coefficients between diploid and polyploidy genomic relationship matrices. In addition, polyGDBLUP did not produce higher prediction accuracy than polyGBLUP for most traits of real data as dominant genetic variance was not captured for these traits. Our study will be a significant promising method for genomic prediction of autopolyploid species.

Origin and diversity of Capsella bursa-pastoris from the genomic point of view.

BACKGROUND: Capsella bursa-pastoris, a cosmopolitan weed of hybrid origin, is an emerging model object for the study of early consequences of polyploidy, being a fast growing annual and a close relative of Arabidopsis thaliana. The development of this model is hampered by the absence of a reference genome sequence. RESULTS: We present here a subgenome-resolved chromosome-scale assembly and a genetic map of the genome of Capsella bursa-pastoris. It shows that the subgenomes are mostly colinear, with no massive deletions, insertions, or rearrangements in any of them. A subgenome-aware annotation reveals the lack of genome dominance-both subgenomes carry similar number of genes. While most chromosomes can be unambiguously recognized as derived from either paternal or maternal parent, we also found homeologous exchange between two chromosomes. It led to an emergence of two hybrid chromosomes; this event is shared between distant populations of C. bursa-pastoris. The whole-genome analysis of 119 samples belonging to C. bursa-pastoris and its parental species C. grandiflora/rubella and C. orientalis reveals introgression from C. orientalis but not from C. grandiflora/rubella. CONCLUSIONS: C. bursa-pastoris does not show genome dominance. In the earliest stages of evolution of this species, a homeologous exchange occurred; its presence in all present-day populations of C. bursa-pastoris indicates on a single origin of this species. The evidence coming from whole-genome analysis challenges the current view that C. grandiflora/rubella was a direct progenitor of C. bursa-pastoris; we hypothesize that it was an extinct (or undiscovered) species sister to C. grandiflora/rubella.

Introgressions lead to reference bias in wheat RNA-seq analysis.

BACKGROUND: RNA-seq is a fundamental technique in genomics, yet reference bias, where transcripts derived from non-reference alleles are quantified less accurately, can undermine the accuracy of RNA-seq quantification and thus the conclusions made downstream. Reference bias in RNA-seq analysis has yet to be explored in complex polyploid genomes despite evidence that they are often a complex mosaic of wild relative introgressions, which introduce blocks of highly divergent genes. RESULTS: Here we use hexaploid wheat as a model complex polyploid, using both simulated and experimental data to show that RNA-seq alignment in wheat suffers from widespread reference bias which is largely driven by divergent introgressed genes. This leads to underestimation of gene expression and incorrect assessment of homoeologue expression balance. By incorporating gene models from ten wheat genome assemblies into a pantranscriptome reference, we present a novel method to reduce reference bias, which can be readily scaled to capture more variation as new genome and transcriptome data becomes available. CONCLUSIONS: This study shows that the presence of introgressions can lead to reference bias in wheat RNA-seq analysis. Caution should be exercised by researchers using non-sample reference genomes for RNA-seq alignment and novel methods, such as the one presented here, should be considered.

MITF regulates the subcellular location of HIF1α through SUMOylation to promote the invasion and metastasis of daughter cells derived from polyploid giant cancer cells.

High concentrations of cobalt chloride (CoCl2) can induce the formation of polyploid giant cancer cells (PGCCs) in various tumors, which can produce daughter cells with strong proliferative, migratory and invasive abilities via asymmetric division. To study the role of hypoxia­inducible factor (HIF) 1α in the formation of PGCCs, colon cancer cell lines Hct116 and LoVo were used as experimental subjects. Western blotting, nuclear and cytoplasmic protein extraction and immunocytochemical experiments were used to compare the changes in the expression and subcellular localization of HIF1α, microphthalmia­associated transcription factor (MITF), protein inhibitor of activated STAT protein 4 (PIAS4) and von Hippel­Lindau disease tumor suppressor (VHL) after treatment with CoCl2. The SUMOylation of HIFα was verified by co­immunoprecipitation assay. After inhibiting HIF1α SUMOylation, the changes in proliferation, migration and invasion abilities of Hct116 and LoVo were compared by plate colony formation, wound healing and Transwell migration and invasion. In addition, lysine sites that led to SUMOylation of HIF1α were identified through site mutation experiments. The results showed that CoCl2 can induce the formation of PGCCs with the expression level of HIF1α higher in treated cells than in control cells. HIF1α was primarily located in the cytoplasm of control cell. Following CoCl2 treatment, the subcellular localization of HIF1α was primarily in the nuclei of PGCCs with daughter cells (PDCs). After treatment with SUMOylation inhibitors, the nuclear HIF1α expression in PDCs decreased. Furthermore, their proliferation, migration and invasion abilities also decreased. After inhibiting the expression of MITF, the expression of HIF1α decreased. MITF can regulate HIF1α SUMOylation. Expression and subcellular localization of VHL and HIF1α did not change following PIAS4 knockdown. SUMOylation of HIF1α occurs at the amino acid sites K391 and K477 in PDCs. After mutation of the two sites, nuclear expression of HIF1α in PDCs was reduced, along with a significant reduction in the proliferation, migration and invasion abilities. In conclusion, the post­translation modification regulated the subcellular location of HIF1α and the nuclear expression of HIF1α promoted the proliferation, migration and invasion abilities of PDCs. MITF could regulate the transcription and protein levels of HIF1α and participate in the regulation of HIF1α SUMOylation.

The emerging importance of cross-ploidy hybridisation and introgression.

Natural hybridisation is now recognised as pervasive in its occurrence across the Tree of Life. Resurgent interest in natural hybridisation fuelled by developments in genomics has led to an improved understanding of the genetic factors that promote or prevent species cross-mating. Despite this body of work overturning many widely held assumptions about the genetic barriers to hybridisation, it is still widely thought that ploidy differences between species will be an absolute barrier to hybridisation and introgression. Here, we revisit this assumption, reviewing findings from surveys of polyploidy and hybridisation in the wild. In a case study in the British flora, 203 hybrids representing 35% of hybrids with suitable data have formed via cross-ploidy matings, while a wider literature search revealed 59 studies (56 in plants and 3 in animals) in which cross-ploidy hybridisation has been confirmed with genetic data. These results show cross-ploidy hybridisation is readily overlooked, and potentially common in some groups. General findings from these studies include strong directionality of hybridisation, with introgression usually towards the higher ploidy parent, and cross-ploidy hybridisation being more likely to involve allopolyploids than autopolyploids. Evidence for adaptive introgression across a ploidy barrier and cases of cross-ploidy hybrid speciation shows the potential for important evolutionary outcomes.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication.

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

Polyploidization: A Biological Force That Enhances Stress Resistance.

Organisms with three or more complete sets of chromosomes are designated as polyploids. Polyploidy serves as a crucial pathway in biological evolution and enriches species diversity, which is demonstrated to have significant advantages in coping with both biotic stressors (such as diseases and pests) and abiotic stressors (like extreme temperatures, drought, and salinity), particularly in the context of ongoing global climate deterioration, increased agrochemical use, and industrialization. Polyploid cultivars have been developed to achieve higher yields and improved product quality. Numerous studies have shown that polyploids exhibit substantial enhancements in cell size and structure, physiological and biochemical traits, gene expression, and epigenetic modifications compared to their diploid counterparts. However, some research also suggested that increased stress tolerance might not always be associated with polyploidy. Therefore, a more comprehensive and detailed investigation is essential to complete the underlying stress tolerance mechanisms of polyploids. Thus, this review summarizes the mechanism of polyploid formation, the polyploid biochemical tolerance mechanism of abiotic and biotic stressors, and molecular regulatory networks that confer polyploidy stress tolerance, which can shed light on the theoretical foundation for future research.

An insight into the gene expression evolution in Gossypium species based on the leaf transcriptomes.

BACKGROUND: Gene expression pattern is associated with biological phenotype and is widely used in exploring gene functions. Its evolution is also crucial in understanding species speciation and divergence. The genus Gossypium is a bona fide model for studying plant evolution and polyploidization. However, the evolution of gene expression during cotton species divergence has yet to be extensively discussed. RESULTS: Based on the seedling leaf transcriptomes, this work analyzed the transcriptomic content and expression patterns across eight cotton species, including six diploids and two natural tetraploids. Our findings indicate that, while the biological function of these cotton transcriptomes remains largely conserved, there has been significant variation in transcriptomic content during species divergence. Furthermore, we conducted a comprehensive analysis of expression distances across cotton species. This analysis lends further support to the use of G. arboreum as a substitute for the A-genome donor of natural cotton polyploids. Moreover, our research highlights the evolution of stress-responsive pathways, including hormone signaling, fatty acid degradation, and flavonoid biosynthesis. These processes appear to have evolved under lower selection pressures, presumably reflecting their critical role in the adaptations of the studied cotton species to diverse environments. CONCLUSIONS: In summary, this study provided insights into the gene expression variation within the genus Gossypium and identified essential genes/pathways whose expression evolution was closely associated with the evolution of cotton species. Furthermore, the method of characterizing genes and pathways under unexpected high or slow selection pressure can also serve as a new strategy for gene function exploration.