Exploring Lignin Biosynthesis Genes in Rice: Evolution, Function, and Expression.

Lignin is nature's second most abundant vascular plant biopolymer, playing significant roles in mechanical support, water transport, and stress responses. This study identified 90 lignin biosynthesis genes in rice based on phylogeny and motif constitution, and they belong to PAL, C4H, 4CL, HCT, C3H, CCoAOMT, CCR, F5H, COMT, and CAD families. Duplication events contributed largely to the expansion of these gene families, such as PAL, CCoAOMT, CCR, and CAD families, mainly attributed to tandem and segmental duplication. Microarray data of 33 tissue samples covering the entire life cycle of rice suggested fairly high PAL, HCT, C3H, CCoAOMT, CCR, COMT, and CAD gene expressions and rather variable C4H, 4CL, and F5H expressions. Some members of lignin-related genes (OsCCRL11, OsHCT1/2/5, OsCCoAOMT1/3/5, OsCOMT, OsC3H, OsCAD2, and OsPAL1/6) were expressed in all tissues examined. The expression patterns of lignin-related genes can be divided into two major groups with eight subgroups, each showing a distinct co-expression in tissues representing typically primary and secondary cell wall constitutions. Some lignin-related genes were strongly co-expressed in tissues typical of secondary cell walls. Combined HPLC analysis showed increased lignin monomer (H, G, and S) contents from young to old growth stages in five genotypes. Based on 90 genes' microarray data, 27 genes were selected for qRT-PCR gene expression analysis. Four genes (OsPAL9, OsCAD8C, OsCCR8, and OsCOMTL4) were significantly negatively correlated with lignin monomers. Furthermore, eleven genes were co-expressed in certain genotypes during secondary growth stages. Among them, six genes (OsC3H, OsCAD2, OsCCR2, OsCOMT, OsPAL2, and OsPAL8) were overlapped with microarray gene expressions, highlighting their importance in lignin biosynthesis.

The genomes of Australian wild limes.

Australian wild limes occur in highly diverse range of environments and are a unique genetic resource within the genus Citrus. Here we compare the haplotype-resolved genome assemblies of six Australian native limes, including four new assemblies generated using PacBio HiFi and Hi-C sequencing data. The size of the genomes was between 315 and 391 Mb with contig N50s from 29.5 to 35 Mb. Gene completeness of the assemblies was estimated to be from 98.4 to 99.3% and the annotations from 97.7 to 98.9% based upon BUSCO, confirming the high contiguity and completeness of the assembled genomes. High collinearity was observed among the genomes and the two haplotype assemblies for each species. Gene duplication and evolutionary analysis demonstrated that the Australian citrus have undergone only one ancient whole-genome triplication event during evolution. The highest number of species-specific and expanded gene families were found in C. glauca and they were primarily enriched in purine, thiamine metabolism, amino acids and aromatic amino acids metabolism which might help C. glauca to mitigate drought, salinity, and pathogen attacks in the drier environments in which this species is found. Unique genes related to terpene biosynthesis, glutathione metabolism, and toll-like receptors in C. australasica, and starch and sucrose metabolism genes in both C. australis and C. australasica might be important candidate genes for HLB tolerance in these species. Expanded gene families were not lineage specific, however, a greater number of genes related to plant-pathogen interactions, predominantly disease resistant protein, was found in C. australasica and C. australis.

dgfr: an R package to assess sequence diversity of gene families.

BACKGROUND: Gene families are groups of homologous genes that often have similar biological functions. These families are formed by gene duplication events throughout evolution, resulting in multiple copies of an ancestral gene. Over time, these copies can acquire mutations and structural variations, resulting in members that may vary in size, motif ordering and sequence. Multigene families have been described in a broad range of organisms, from single-celled bacteria to complex multicellular organisms, and have been linked to an array of phenomena, such as host-pathogen interactions, immune evasion and embryonic development. Despite the importance of gene families, few approaches have been developed for estimating and graphically visualizing their diversity patterns and expression profiles in genome-wide studies. RESULTS: Here, we introduce an R package named dgfr, which estimates and enables the visualization of sequence divergence within gene families, as well as the visualization of secondary data such as gene expression. The package takes as input a multi-fasta file containing the coding sequences (CDS) or amino acid sequences from a multigene family, performs a pairwise alignment among all sequences, and estimates their distance, which is subjected to dimension reduction, optimal cluster determination, and gene assignment to each cluster. The result is a dataset that allows for the visualization of sequence divergence and expression within the gene family, an approximation of the number of clusters present in the family. CONCLUSIONS: dgfr provides a way to estimate and study the diversity of gene families, as well as visualize the dispersion and secondary profile of the sequences. The dgfr package is available at https://github.com/lailaviana/dgfr under the GPL-3 license.

Genome-wide identification of walnut (Juglans regia) PME gene family members and expression analysis during infection with Cryptosphaeria pullmanensis pathogens.

Walnuts exhibit a higher resistance to diseases, though they are not completely immune. This study focuses on the Pectin methylesterase (PME) gene family to investigate whether it is involved in disease resistance in walnuts. These 21 genes are distributed across 12 chromosomes, with four pairs demonstrating homology. Variations in conserved motifs and gene structures suggest diverse functions within the gene family. Phylogenetic and collinear gene pairs of the PME family indicate that the gene family has evolved in a relatively stable way. The cis-acting elements and gene ontology enrichment of these genes, underscores their potential role in bolstering walnuts' defense mechanisms. Transcriptomic analyses were conducted under conditions of Cryptosphaeria pullmanensis infestation and verified by RT-qPCR. The results showed that certain JrPME family genes were activated in response, leading to the hypothesis that some members may confer resistance to the disease.

Under Stress: Searching for Genes Involved in the Response of Abies pinsapo Boiss to Climate Change.

Currently, Mediterranean forests are experiencing the deleterious effects of global warming, which mainly include increased temperatures and decreased precipitation in the region. Relict Abies pinsapo fir forests, endemic in the southern Iberian Peninsula, are especially sensitive to these recent environmental disturbances, and identifying the genes involved in the response of this endangered tree species to climate-driven stresses is of paramount importance for mitigating their effects. Genomic resources for A. pinsapo allow for the analysis of candidate genes reacting to warming and aridity in their natural habitats. Several members of the complex gene families encoding late embryogenesis abundant proteins (LEAs) and heat shock proteins (HSPs) have been found to exhibit differential expression patterns between wet and dry seasons when samples from distinct geographical locations and dissimilar exposures to the effects of climate change were analyzed. The observed changes were more perceptible in the roots of trees, particularly in declining forests distributed at lower altitudes in the more vulnerable mountains. These findings align with previous studies and lay the groundwork for further research on the molecular level. Molecular and genomic approaches offer valuable insights for mitigating climate stress and safeguarding this endangered conifer.

A two-sequence motif-based method for the inventory of gene families in fragmented and poorly annotated genome sequences.

Databases of genome sequences are growing exponentially, but, in some cases, assembly is incomplete and genes are poorly annotated. For evolutionary studies, it is important to identify all members of a given gene family in a genome. We developed a method for identifying most, if not all, members of a gene family from raw genomes in which assembly is of low quality, using the P-type ATPase superfamily as an example. The method is based on the translation of an entire genome in all six reading frames and the co-occurrence of two family-specific sequence motifs that are in close proximity to each other. To test the method's usability, we first used it to identify P-type ATPase members in the high-quality annotated genome of barley (Hordeum vulgare). Subsequently, after successfully identifying plasma membrane H+-ATPase family members (P3A ATPases) in various plant genomes of varying quality, we tested the hypothesis that the number of P3A ATPases correlates with the ability of the plant to tolerate saline conditions. In 19 genomes of glycophytes and halophytes, the total number of P3A ATPase genes was found to vary from 7 to 22, but no significant difference was found between the two groups. The method successfully identified P-type ATPase family members in raw genomes that are poorly assembled.

Chromosome-level genome of Ambrosia trifida provides insights into adaptation and the evolution of pollen allergens.

Ambrosia trifida (giant ragweed) is an invasive plant that can cause serious damage to natural ecosystems and severe respiratory allergies. However, the genomic basis of invasive adaptation and pollen allergens in Ambrosia species remain largely unknown. Here, we present a 1.66 Gb chromosome-scale reference genome for giant ragweed and identified multiple types of genome duplications, which are responsible for its rapid environmental adaptation and pollen development. The largest copies number and species-specific expansions of resistance-related gene families compared to Heliantheae alliance might contribute to resist stresses, pathogens and rapid adaptation. To extend the knowledge of evolutionary process of allergic pollen proteins, we predicted 26 and 168 potential pollen allergen candidates for giant ragweed and other Asteraceae plant species by combining machine learning and identity screening. Interestingly, we observed a specific tandemly repeated array for potential allergenic pectate lyases among Ambrosia species. Rapid evolutionary rates on putative pectate lyase allergens may imply a crucial role of nonsynonymous mutations on amino acid residues for plant biological function and allergenicity. Altogether, this study provides insight into the molecular ecological adaptation and putative pollen allergens prediction that will be helpful in promoting invasion genomic research and evolution of putative pollen allergy in giant ragweed.

Overexpression of BnMYBL2-1 improves plant drought tolerance via the ABA-dependent pathway.

Drought stress is a major environmental challenge that poses considerable threats to crop survival and growth. Previous research has indicated anthocyanins play a crucial role in alleviating oxidative damage, photoprotection, membrane stabilization, and water retention under drought stress. However, the presence of MYBL2 (MYELOBBLASTOSIS LIKE 2), an R3-MYB transcription factor (TF) which known to suppress anthocyanin biosynthesis. In this study, four BnMYBL2 members were cloned from Brassica napus L, and BnMYBL2-1 was overexpressed in Triticum aestivum L (No BnMYBL2 homologous gene was detected in wheat). Subsequently, the transgenic wheat lines were treated with drought, ABA and anthocyanin. Results showed that transgenic lines exhibited greater drought tolerance compared to the wild-type (WT), characterized by improved leaf water content (LWC), elevated levels of soluble sugars and chlorophyll, and increased antioxidant enzyme activity. Notably, transgenic lines also exhibited significant upregulation in abscisic acid (ABA) content, along with the transcriptional levels of key enzymes involved in ABA signalling under drought. Results also demonstrated that BnMYBL2-1 promoted the accumulation of ABA and anthocyanins in wheat. Overall, the study highlights the positive role of BnMYBL2-1 in enhancing crop drought tolerance through ABA signalling and establishes its close association with anthocyanin biosynthesis. These findings offer valuable insights for the development of drought-resistant crop varieties and enhance the understanding of the molecular mechanisms underlying plant responses to drought stress.

PEO: Plant Expression Omnibus - a comparative transcriptomic database for 103 Archaeplastida.

The Plant Expression Omnibus (PEO) is a web application that provides biologists with access to gene expression insights across over 100 plant species, ~60 000 manually annotated RNA-seq samples, and more than 4 million genes. The tool allows users to explore the expression patterns of genes across different organs, identify organ-specific genes, and discover top co-expressed genes for any gene of interest. PEO also provides functional annotations for each gene, allowing for the identification of genetic modules and pathways. PEO is designed to facilitate comparative kingdom-wide gene expression analysis and provide a valuable resource for plant biology research. We provide two case studies to demonstrate the utility of PEO in identifying candidate genes in pollen coat biosynthesis in Arabidopsis and investigating the biosynthetic pathway components of capsaicin in Capsicum annuum. The database is freely available at https://expression.plant.tools/.

Comparative genomic analysis of chemosensory-related gene families in gastropods.

Chemoreception is critical for the survival and reproduction of animals. Except for a reduced group of insects and chelicerates, the molecular identity of chemosensory proteins is poorly understood in invertebrates. Gastropoda is the extant mollusk class with the greatest species richness, including marine, freshwater, and terrestrial lineages, and likely, highly diverse chemoreception systems. Here, we performed a comprehensive comparative genome analysis taking advantage of the chromosome-level information of two Gastropoda species, one of which belongs to a lineage that underwent a whole genome duplication event. We identified thousands of previously uncharacterized chemosensory-related genes, the majority of them encoding G protein-coupled receptors (GPCR), mostly organized into clusters distributed across all chromosomes. We also detected gene families encoding degenerin epithelial sodium channels (DEG-ENaC), ionotropic receptors (IR), sensory neuron membrane proteins (SNMP), Niemann-Pick type C2 (NPC2) proteins, and lipocalins, although with a lower number of members. Our phylogenetic analysis of the GPCR gene family across protostomes revealed: (i) remarkable gene family expansions in Gastropoda; (ii) clades including members from all protostomes; and (iii) species-specific clades with a substantial number of receptors. For the first time, we provide new and valuable knowledge into the evolution of the chemosensory gene families in invertebrates other than arthropods.

PanDelos-frags: A methodology for discovering pangenomic content of incomplete microbial assemblies.

Pangenomics was originally defined as the problem of comparing the composition of genes into gene families within a set of bacterial isolates belonging to the same species. The problem requires the calculation of sequence homology among such genes. When combined with metagenomics, namely for human microbiome composition analysis, gene-oriented pangenome detection becomes a promising method to decipher ecosystem functions and population-level evolution. Established computational tools are able to investigate the genetic content of isolates for which a complete genomic sequence is available. However, there is a plethora of incomplete genomes that are available on public resources, which only a few tools may analyze. Incomplete means that the process for reconstructing their genomic sequence is not complete, and only fragments of their sequence are currently available. However, the information contained in these fragments may play an essential role in the analyses. Here, we present PanDelos-frags, a computational tool which exploits and extends previous results in analyzing complete genomes. It provides a new methodology for inferring missing genetic information and thus for managing incomplete genomes. PanDelos-frags outperforms state-of-the-art approaches in reconstructing gene families in synthetic benchmarks and in a real use case of metagenomics. PanDelos-frags is publicly available at https://github.com/InfOmics/PanDelos-frags.

Chromosome-level genome assembly of sea scallop Placopecten magellanicus provides insights into the genetic characteristics and adaptive evolution of large scallops.

Placopecten magellanicus (Gmelin, 1791), a deep-sea Atlantic scallop, holds significant commercial value as a benthic marine bivalve along the northwest Atlantic coast. Recognizing its economic importance, the need to reconstruct its genome assembly becomes apparent, fostering insights into natural resources and generic breeding potential. This study reports a high-quality chromosome-level genome of P. magellanicus, achieved through the integration of Illumina short read sequencing, PacBio HiFi sequencing, and Hi-C sequencing techniques. The resulting assembly spans 1778 Mb with a scaffold N50 of 86.71 Mb. An intriguing observation arises - the genome size of P. magellanicus surpasses that of its Pectinidae family peers by 1.80 to 2.46 times. Within this genome, 28,111 protein-coding genes were identified. Comparative genomic analysis involving five scallop species unveils the critical determinant of this expanded genome: the proliferation of repetitive sequences recently inserted, contributing to its enlarged size. The landscape of whole genome collinearity sheds light on the relationships among scallop species, enhancing our broader understanding of their genomic framework. This genome provides genomic resources for future molecular biology research on scallops and serves as a guide for the exploration of longevity-related genes in scallops.

U2 and U4 snDNA Comparative Chromosomal Mapping in the Neotropical Fish Genera Apareiodon and Parodon (Characiformes: Parodontidae).

Small nuclear DNA (snDNA) are valuable cytogenetic markers for comparative studies in chromosome evolution because different distribution patterns were found among species. Parodontidae, a Neotropical fish family, is known to have female heterogametic sex chromosome systems in some species. The U2 and U4 snDNA sites have been found to be involved in Z and W chromosome differentiation in Apareiodon sp., Apareiodon affinis, and Parodon hilarii. However, few studies have evaluated snDNA sites as propulsors of chromosome diversification among closely related fish species. In this study, we investigated the distribution of U2 and U4 snDNA clusters in the chromosomes of 10 populations/species belonging to Apareiodon and Parodon, aiming to identify chromosomal homeologies or diversification. In situ localization data revealed a submetacentric pair carrying the U2 snDNA site among the populations/species analyzed. Furthermore, all studied species demonstrated homeology in the location of U4 snDNA cluster in the proximal region of metacentric pair 1, besides an additional signal showing up with a divergence in Apareiodon. Comparative chromosomal mapping of U4 snDNA also helped to reinforce the proposal of the ZZ/ZW1W2 sex chromosome system origin in an A. affinis population. According to cytogenetic data, the study corroborates the diversification in Parodontidae paired species with uncertain taxonomy.

Mechanisms of Karyotypic Diversification in Ancistrus (Siluriformes, Loricariidae): Inferences from Repetitive Sequence Analysis.

Ancistrus is a highly diverse neotropical fish genus that exhibits extensive chromosomal variability, encompassing karyotypic morphology, diploid chromosome number (2n = 34-54), and the evolution of various types of sex chromosome systems. Robertsonian rearrangements related to unstable chromosomal sites are here described. Here, the karyotypes of two Ancistrus species were comparatively analyzed using classical cytogenetic techniques, in addition to isolation, cloning, sequencing, molecular characterization, and fluorescence in situ hybridization of repetitive sequences (i.e., 18S and 5S rDNA; U1, U2, and U5 snDNA; and telomere sequences). The species analyzed here have different karyotypes: Ancistrus sp. 1 (2n = 38, XX/XY) and Ancistrus cirrhosus (2n = 34, no heteromorphic sex chromosomes). Comparative mapping showed different organizations for the analyzed repetitive sequences: 18S and U1 sequences occurred in a single site in all populations of the analyzed species, while 5S and U2 sequences could occur in single or multiple sites. A sequencing analysis confirmed the identities of the U1, U2, and U5 snDNA sequences. Additionally, a syntenic condition for U2-U5 snDNA was found in Ancistrus. In a comparative analysis, the sequences of rDNA and U snDNA showed inter- and intraspecific chromosomal diversification. The occurrence of Robertsonian rearrangements and other dispersal mechanisms of repetitive sequences are discussed.

Evolution of Three-Finger Toxin Genes in Neotropical Colubrine Snakes (Colubridae).

Snake venom research has historically focused on front-fanged species (Viperidae and Elapidae), limiting our knowledge of venom evolution in rear-fanged snakes across their ecologically diverse phylogeny. Three-finger toxins (3FTxs) are a known neurotoxic component in the venoms of some rear-fanged snakes (Colubridae: Colubrinae), but it is unclear how prevalent 3FTxs are both in expression within venom glands and more broadly among colubrine species. Here, we used a transcriptomic approach to characterize the venom expression profiles of four species of colubrine snakes from the Neotropics that were dominated by 3FTx expression (in the genera Chironius, Oxybelis, Rhinobothryum, and Spilotes). By reconstructing the gene trees of 3FTxs, we found evidence of putative novel heterodimers in the sequences of Chironius multiventris and Oxybelis aeneus, revealing an instance of parallel evolution of this structural change in 3FTxs among rear-fanged colubrine snakes. We also found positive selection at sites within structural loops or "fingers" of 3FTxs, indicating these areas may be key binding sites that interact with prey target molecules. Overall, our results highlight the importance of exploring the venoms of understudied species in reconstructing the full evolutionary history of toxins across the tree of life.

The Fish Pathogen "Candidatus Clavichlamydia salmonicola"-A Missing Link in the Evolution of Chlamydial Pathogens of Humans.

Chlamydiae like Chlamydia trachomatis and Chlamydia psittaci are well-known human and animal pathogens. Yet, the chlamydiae are a much larger group of evolutionary ancient obligate intracellular bacteria that includes predominantly symbionts of protists and diverse animals. This makes them ideal model organisms to study evolutionary transitions from symbionts in microbial eukaryotes to pathogens of humans. To this end, comparative genome analysis has served as an important tool. Genome sequence data for many chlamydial lineages are, however, still lacking, hampering our understanding of their evolutionary history. Here, we determined the first high-quality draft genome sequence of the fish pathogen "Candidatus Clavichlamydia salmonicola", representing a separate genus within the human and animal pathogenic Chlamydiaceae. The "Ca. Clavichlamydia salmonicola" genome harbors genes that so far have been exclusively found in Chlamydia species suggesting that basic mechanisms important for the interaction with chordate hosts have evolved stepwise in the history of chlamydiae. Thus, the genome sequence of "Ca. Clavichlamydia salmonicola" allows to constrain candidate genes to further understand the evolution of chlamydial virulence mechanisms required to infect mammals.

Expansion of photoreception-related gene families may drive ecological adaptation of the dominant diatom species Skeletonema marinoi.

Diatom species of the genus Skeletonema are dominant in global coastal waters with important roles in marine primary production and global biogeochemical cycling. Many Skeletonema species have been extensively studied also because they can cause harmful algae blooms (HABs) with negative impacts on marine ecosystems and aquaculture. In this study, the first chromosome-level assembly of the genome of Skeletonema marinoi was constructed. The genome size was 64.99 Mb with a contig N50 of 1.95 Mb. Up to 97.12 % of contigs were successfully anchored on 24 chromosomes. Analysis of the annotated genes revealed 28 large syntenic blocks with 2397 collinear gene pairs in the genome of S. marinoi, suggesting large-scale segmental duplication events in evolution. Substantial expansion of light-harvesting genes encoding fucoxanthin-chlorophyll a/c binding proteins, as well as expansion of photoreceptor gene families encoding aureochromes and cyptochromes (CRY) in S. marinoi were found, which may have shaped ecological adaptation of S. marinoi. In conclusion, the construction of the first high-quality Skeletonema genome assembly offers valuable clues on the ecological and evolutionary characteristics of this dominant coastal diatom species.

The VGNC: expanding standardized vertebrate gene nomenclature.

The Vertebrate Gene Nomenclature Committee (VGNC) was established in 2016 as a sister project to the HUGO Gene Nomenclature Committee, to approve gene nomenclature in vertebrate species without an existing dedicated nomenclature committee. The VGNC aims to harmonize gene nomenclature across selected vertebrate species in line with human gene nomenclature, with orthologs assigned the same nomenclature where possible. This article presents an overview of the VGNC project and discussion of key findings resulting from this work to date. VGNC-approved nomenclature is accessible at https://vertebrate.genenames.org and is additionally displayed by the NCBI, Ensembl, and UniProt databases.

Genome-Wide Identification of the MAPK and MAPKK Gene Families in Response to Cold Stress in Prunus mume.

Protein kinases of the MAPK cascade family (MAPKKK-MAPKK-MAPK) play an essential role in plant stress response and hormone signal transduction. However, their role in the cold hardiness of Prunus mume (Mei), a class of ornamental woody plant, remains unclear. In this study, we use bioinformatic approaches to assess and analyze two related protein kinase families, namely, MAP kinases (MPKs) and MAPK kinases (MKKs), in wild P. mume and its variety P. mume var. tortuosa. We identify 11 PmMPK and 7 PmMKK genes in the former species and 12 PmvMPK and 7 PmvMKK genes in the latter species, and we investigate whether and how these gene families contribute to cold stress responses. Members of the MPK and MKK gene families located on seven and four chromosomes of both species are free of tandem duplication. Four, three, and one segment duplication events are exhibited in PmMPK, PmvMPK, and PmMKK, respectively, suggesting that segment duplications play an essential role in the expansion and evolution of P. mume and its gene variety. Moreover, synteny analysis suggests that most MPK and MKK genes have similar origins and involved similar evolutionary processes in P. mume and its variety. A cis-acting regulatory element analysis shows that MPK and MKK genes may function in P. mume and its variety's development, modulating processes such as light response, anaerobic induction, and abscisic acid response as well as responses to a variety of stresses, such as low temperature and drought. Most PmMPKs and PmMKKs exhibited tissue-specifific expression patterns, as well as time-specific expression patterns that protect them through cold. In a low-temperature treatment experiment with the cold-tolerant cultivar P. mume 'Songchun' and the cold-sensitive cultivar 'Lve', we find that almost all PmMPK and PmMKK genes, especially PmMPK3/5/6/20 and PmMKK2/3/6, dramatically respond to cold stress as treatment duration increases. This study introduces the possibility that these family members contribute to P. mume's cold stress response. Further investigation is warranted to understand the mechanistic functions of MAPK and MAPKK proteins in P. mume development and response to cold stress.

Comparative genomics revealed drastic gene difference in two small Chinese perches, Siniperca undulata and S. obscura.

Siniperca undulata and S. obscura (Centrarchiformes: Sinipercidae) are small Chinese perches, living in creeks and streams in southern China. While they have sympatric distribution and occupy similar macrohabitat, their body sizes and ecological niches have many differences. Determining the genome sequences of S. undulata and S. obscura would provide us an essential data set for better understanding their genetic makeup and differences that may play important roles in their adaptation to different niches. We determined the genome sequences of both S. undulata and S. obscura using 10× genomics technology and the next-generation sequencing. The assembled genomes of S. undulata and S. obscura were 744 and 733 Mb, respectively. Gene family analysis revealed that there were no overlap between S. undulata and S. obscura in terms of rapid expanding and rapid contracting genes families, which were related to growth, immunity, and mobility. Positive selection analyses also cooperated that the function of selected genes involve growth, athletic ability, and immunity, which may explain the preference of different niches by S. undulata and S. obscura. Pairwise sequentially Markovian coalescent analyses for the two species suggested that populations of both S. undulata and S. obscura showed a rising trend between 90 and 70 Ka probably due to the mild environment during the last interglacial period. A stage of population shrinking occurred from 70 to 20 Ka, which was in with the Tali glacial period in eastern China (57-16 Ka).