Next Article in Journal
Genome-Wide Characterization and Sequence Polymorphism Analyses of Glycine max Fibrillin (FBN) Revealed Its Role in Response to Drought Condition
Next Article in Special Issue
Ancient Mitochondrial Genomes Provide New Clues to the Origin of Domestic Cattle in China
Previous Article in Journal
Genome-Wide Identification of NAC Gene Family and Expression Analysis under Abiotic Stresses in Avena sativa
Previous Article in Special Issue
Deciphering the West Eurasian Genetic Footprints in Ancient South India
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ancient Mitogenomes Reveal Stable Genetic Continuity of the Holocene Serows

1
School of Environmental Studies, China University of Geosciences, Wuhan 430078, China
2
State Key Laboratory of Biogeology and Environmental Geology, China University of Geosciences, Wuhan 430078, China
3
School of Earth Science, China University of Geosciences, Wuhan 430074, China
4
Guizhou Institute of Geological Survey, Guiyang 550081, China
5
Palaeontological Fossil Conservation Center, Qinggang County, Suihua 151600, China
6
Faculty of Materials Science and Chemistry, China University of Geosciences, Wuhan 430078, China
*
Author to whom correspondence should be addressed.
Genes 2023, 14(6), 1187; https://doi.org/10.3390/genes14061187
Submission received: 3 May 2023 / Revised: 24 May 2023 / Accepted: 26 May 2023 / Published: 29 May 2023
(This article belongs to the Special Issue Advances in Ancient Genomes)

Abstract

:
As one of the remaining species of Caprinae only found in Asia, serows (Capricornis) and their classification and conservation have received increasing attention in recent years. However, their evolutionary history and population dynamics are not yet clear. To shed light on these topics, we report the first near-complete ancient mitochondrial genomes from two serow sub-fossils (CADG839 and CADG946) dating to 8860 ± 30 years and 2450 ± 30 years, and incorporate the newly obtained mitogenomes into the dataset of living serows (18 complete mitochondrial genomes drawn from National Center for Biotechnology Information, NCBI) to investigate their relationships and evolution. Phylogenetic results support four clades of serows that can be further divided into five subclades, indicating higher genetic diversity than previously thought. Notably, our two ancient samples do not form a separate branch but belong to Capricornis sumatraensis clade A together with modern individuals, which suggests genetic continuity between ancient and modern serows. Furthermore, our results suggest that the maternal divergences of serows occurred at the beginning of the Pleistocene. Bayesian estimation indicates that the first divergence among all serows happened approximately 2.37 Ma (95% highest posterior density, HPD: 2.74–2.02 Ma) when Japanese serow (Capricornis crispus) appeared, while the last divergence occurred within the Sumatran serow (C. sumatraensis clade A and B) around 0.37–0.25 Ma. Additionally, we found the effective maternal population size of C. sumatraensis increased around 225–160 and 90–50 ka, then remained stable since 50 ka. Overall, our study provides new insights into serow phylogeny and evolutionary history.

1. Introduction

Serows belongs to the genus Capricornis Ogilby, 1837, which is within the Caprinae Gray, 1821, family and the Bovidae Gray, 1821, family. This genus is typically represented by tropical and subtropical species. It possibly originated in Central Asia in the mid to late Pliocene (about 3.3–2.5 Ma). The extant serows mostly live in China, Japan and Southeast Asia [1,2,3,4,5]. As the most widely distributed member of serows in China and Southeast Asia, the Sumatran serow (C. sumatraensis) once reached Northern China (39° N), while the northernmost group of extant mainland serow lives at 30° N in China [1,3]. Unfortunately, population isolation due to habitat fragmentation and human hunting has made their population size reduce quickly, and it is listed as a vulnerable (VU A2; C1) species on the International Union for the Conservation of Nature Red List [1].
Based on morphological characteristics, serows were divided into seven species: C. crispus, Capricornis swinhoei, Capricornis rubidus, C. sumatraensis, Capricornis milnedwardsii, Capricornis thar, and Capricornis maritimus [6,7,8]. However, recent investigation based on molecular phylogeny suggested that only four species exist in Capricornis, i.e., C. crispus, C. swinhoei, C. rubidus, and C. sumatraensis [9]. The classification controversy between morphological and molecular explorations mainly centers on the definition of C. sumatraensis. Specifically, in molecular studies, C. milnedwardsii, C. thar, and C. maritimus, which were previously classified as separate species based on morphology, are now included within the C. sumatraensis group [9]. Moreover, fossil evidence indicates that C. sumatraensis once existed in Northern China (Beijing) during the Pleistocene (about 2.58 Ma–11.7 ka), but the current range in China has been contracted to the southern area. Additionally, morphological investigation of C. sumatraensis indicated that their body size tended to increase from the mid-Pleistocene to the late Pleistocene, and then decreased in the Holocene [3].
At the molecular level, both the phylogeny and the classification of different serows have been investigated using mitochondrial sequences [9,10,11]. However, these molecular investigations are limited in extant specimens. The phylogenetic relationship between Quaternary C. sumatraensis and its extant counterparts remains poorly understood. Additionally, it is unclear whether the genetic diversity of C. sumatraensis was affected during the regional extirpation event that occurred at least in Northern China.
In this study, we aim to clarify the interspecific relationships of ancient and modern serows at the molecular level and therefore to explore their population history. We obtained two near-complete mitochondrial genomes (15,577 bp and 15,597 bp, both are out of 16,524 bp) from two ancient serow individuals and performed phylogenetic analyses of modern and ancient serows and updated the divergence time of Capricornis, as well as investigating the effective maternal population of C. sumatraensis. The new findings in this study provide further understanding into the evolutionary history of Capricornis.

2. Materials and Methods

2.1. Samples

Two ancient serow samples were collected from Beijing and Guizhou, China (Figure 1). One sample (CADG839) is stored at the Institute of Vertebrate Paleontology and Paleoanthropology. The other sample (CADG946) is stored at the China University of Geosciences (Wuhan). Two samples were sent for accelerator mass spectrometry (AMS) radiocarbon dating at the BETA laboratory, Miami, USA. The age of CADG839 is 9884–10,160 cal BP, and CADG946 is 2361–2703 cal BP (Supplementary Table S1).

2.2. Sample Handling and DNA Extraction

Prior to DNA extraction, the samples were cleaned with a brush and soaked in 0.5% bleach for 10 min, and then rinsed with nuclear-free water [13]. Each sample was ground into powder with a mortar and approximately 125 mg of powder was added into a 15 mL centrifuge tube that contained 3.0 mL EDTA (0.5 M, pH = 8) and 40 µL proteinase K (20 mg/mL). The extraction buffer with the sample powder was incubated overnight in a rotating hybridization oven at 37 °C. The incubation mix was centrifuged at 7000 rpm for 10 min, and the supernatant was transferred into an ultrafiltration tube (Millipore, Darmstadt, Germany) and condensed to about 100 µL by centrifuging at 7000 rpm for 39 min. Finally, DNA was purified and eluted in 50 µL EB using a MinElute PCR Purification Kit (Qiagen, Hilden, Germany), according to the manual instructions.

2.3. DNA Library Preparation and Sequencing

Each of 20 µL extractions was used for double-stranded Illumina libraries construction, modified from Meyer and Kircher [14]. In the blunt-end repair process, the following components were used: 5 µL NEB buffer 2 (New England Biolabs, Ipswich, UK), 5 µL ATP (New England Biolabs, Ipswich, UK), 2 µL BSA (New England Biolabs, Ipswich, UK), 2 µL T4 DNA polymerase (New England Biolabs, Ipswich, UK), 0.4 µL T4 polynucleotide kinase (New England Biolabs, Ipswich, UK), and 2 µL dNTPs (Tiangen, Beijing, China). The mixture was then supplemented with 20 µL DNA template and nuclease-free water to a final volume of 50 µL. The reaction was incubated at 15 °C for 15 min, followed by an additional incubation at 25 °C for 15 min using a PCR machine. Subsequently, the mixed reaction was purified using the MinElute PCR Purification Kit (Qiagen, Hilden, Germany), and the DNA was eluted in 21 µL of TET buffer. In the adapter ligation step, the following components were used: 1 µL Adapter (dilution 1:20) (Sangon Biotech, Shanghai, China), 1 µL Quick Ligase (New England Biolabs, Ipswich, UK), 20 µL Quick Ligase buffer (New England Biolabs, Ipswich, UK), and 18 µL of template. The reaction mixture was incubated at 22 °C for 15 min using a PCR machine. Subsequently, the mixed reaction was purified using the MinElute PCR purification kit (Qiagen, Hilden, Germany), and the DNA was eluted in 23 µL TET buffer. In addition, for the adaptor fill-in reaction, we used the following components: 4 µL isothermal buffer (New England Biolabs, Ipswich, UK), 2 µL Bst polymerase (New England Biolabs, Ipswich, UK), 2 µL dNTPs, 20 µL DNA template, and water to achieve a final reaction volume of 40 µL. The reaction mixture was incubated at 37 °C for 20 min using a PCR machine, followed by a temperature increase to 80 °C, and incubated for 20 min. The resulting reaction product was used for PCR amplification. For the indexing PCR amplifications, we prepared the reactions with a Q5 Hot Start High-Fidelity 2× Master Mix (New England Biolabs, Ipswich, UK) using the following cycling protocol: an initial denaturation step at 98 °C for 30 s, followed by 17 cycles of denaturation at 98 °C for 10 s, annealing at 60 °C for 75 s, and extension at 60 °C for 6 min. The library quality was measured using Qubit 4.0 (Invitrogen, Carlsbad, CA, USA) and TapeStation 4150 (Agilent, Santa Clara, CA, USA). To detect contamination, blank controls were set during the DNA extraction, library construction, and PCR amplification procedures. Finally, all libraries were sequenced on an Illumina NovaSeq6000 platform at Annoroad Gene Technology Co., Ltd., Beijing, China.

2.4. Data Analysis

The adapter sequence of raw reads was excised with fastp v 0.21.0 [15]. After we had compared the reads mapping to four Capricornis mitochondrial genomes (Supplementary Table S2), we finally mapped the trimmed reads to the C. sumatraensis mitochondrial genome (GenBank No. NC020629) using the “aln” algorithm in a Burrows–Wheeler aligner (BWA v 0.7.17) with default parameters [16]. Reads with a mapping quality below 30 were discarded using “view”, duplicate reads were removed using “rmdup”, and alignments were sorted using “sort” in SAMtools v 1.9 [17]. All BAM files of libraries were then merged, and duplications were removed to produce a combined BAM file using SAMtools v 1.9. Finally, consensus sequences were constructed using “doFasta 3” in ANGSD v. 0.921 [18]. Read coverage across the reference was calculated using Qualimap v 2.2.1 [19]. The base damage to DNA fragments was analyzed via mapDamage2 with default parameters [20]. Please refer to the supplementary material (Supplementary Table S3) for detailed information on the sequencing quality of all libraries.

2.5. Bioinformatics Analysis

To investigate the phylogenetic relationship within the genus Capricornis, our two newly obtained near-complete mitochondrial sequences and 23 mitochondrial genomes retrieved from GenBank (18 available Capricornis, and 5 outgroup sequences that include 3 Naemorhedus goral, 1 Ovibos moschatus, 1 Ovis aries; see Supplementary Tables S4 and S5) were aligned with MAFFT on the CIPRES portal and manually inspected after alignment [21,22]. Homologous sequences were obtained after we had removed the high variable control region (CR). Phylogenetic trees were constructed using neighbor-joining (NJ) and maximum-likelihood (ML) methods in MEGA 11 [23]. The NJ tree was constructed based on a Jukes–Cantor distance. The ML tree was constructed with the heuristic method of the nearest-neighbor-interchange, and the optimal substitution model ”GTR + I + G” was selected through comparison of Bayesian Information Criterion (BIC) scores in a jModelTest v2.1 [24]. A total of 1000 bootstrap replicates was used to ensure the reliability of nodes [25].
In order to reveal the interspecies relationship of serows, a median-joining network was reconstructed with Popart 1.7 using 2 near-complete mitochondrial DNA sequences in this study and 18 complete mitochondrial DNA sequences of Capricornis from the GenBank (Supplementary Table S5) [26]. Only sites covered by all individuals were retained to obtain a length of 14,034 bp dataset.
Bayesian analysis was conducted with the same dataset of the Capricornis phylogeny that includes 25 sequences, using the Monte Carlo Markov chains (MCMC) method implemented in BEAST 1.8.4 [27]. Due to the lack of direct fossil records for calibration, we used a comprehensive evolutionary framework of bovidae to estimate the divergence time between “Capricornis/Naemorhedus” and the most recent common ancestor (TMRCA) of caprini (Supplementary Figure S2), drawing on the method of Bibi [28] and Pérez et al. [29]. Two node ages were used to correct our BEAST analysis: the divergence time between Capricornis and Naemorhedus with a mean of 4.91 Ma, and the TMRCA of caprini with a mean of 9.8 Ma. Meanwhile, the strict clock was set with reference to the mitochondrial substitution rate of caprini (substitution rate of 3.8–5.4 × 10−8 substitutions per nucleotide per year) [30], and a constant population size in parameter settings. MCMC ran for 75 million iterations, sampling every 1000 steps. Tracer v1.6 was used to test the effective sample size (ESS > 200), discarding the first 10% of samples as burn-in [31]. A maximum clade credibility (MCC) tree was visualized using TreeAnnotator v1.5.4 [27], and modified in FigTree v1.4.4 (http://tree.bio.ed.ac.uk/software/figtree, accessed on 1 May 2023). The Bayesian skyline plots (BSPs) of C. sumatraensis were simulated with identical parameters in BEAST v1.8.4 and visualized in Tracer v1.6.3.

3. Results

After adapter trimming and short read removal, a total number of 2102 and 2934 unique reads for CADG839 and CADG946 was successfully mapped onto a C. sumatraensis mitochondrial sequence (NC020629). Two near-complete mitochondrial genomes were obtained from two serow subfossils (15,577 bp for CADG839, 15,597 bp for CADG946) with mean depths of 8.13 and 7.98 folds, respectively. The length of the DNA fragments of our samples was between 35 and 70 bp, with significant C to T misincorporation at the 5′ end of the reads, which is consistent with the typical ancient DNA damage pattern (Supplementary Figure S3) [32].
Both our Bayesian phylogenetic MCC tree (Figure 2) and the ML phylogenetic tree (Supplementary Figure S4) indicate that there are four clades at the molecular level in Capricornis that correspond to four species, i.e., C. crispus, C. swinhoei, C. rubidus, and C. sumatraensis. Additionally, our analysis revealed that the ancient serows in this study are within the maternal diversity of modern ones, i.e., one of our two samples, CADG839, is settled in C. sumatraensis subclade A2 and close to a modern individual from Sichuan province, China; while the other sample, CADG946, is placed in the basal position of subclade A1 that mainly consists of samples from South China.
Regarding the divergence time of different serow clades, we estimated that the first divergence occurred in the C. crispus clade from the serow common ancestor around 2.37 Ma (95% HPD: 2.02–2.74 Ma), followed by the split of the C. swinhoei clade around 1.7 Ma (95% HPD: 1.5–1.96 Ma). The divergence time between the C. rubidus clade and the C. sumatraensis clade was approximately 0.7 Ma (95% HPD: 0.58–0.85 Ma), whereas differentiation within the C. sumatraensis clade began circa 0.3 Ma (95% HPD: 0.25–0.37 Ma) (Figure 2).
The result of haplogroup network analysis further supports the finding of phylogenetic analysis, indicating that the four major haplogroups correspond to the four species delineated at the molecular level, i.e., C. crispus, C. swinhoei, C. rubidus, and C. sumatraensis. Additionally, the C. sumatraensis haplogroup can be further subdivided into two secondary haplogroups (C. sumatraensis A and C. sumatraensis B) (Figure 3).
The Bayesian skyline plot analysis indicates that the effective maternal population of C. sumatraensis increased slightly from about 225 Kya to 160 Kya, followed by a stable stage until 90 Kya, and then went to another increase period until 50 Kya and remained stable in the most recent 50,000 years (Figure 4).

4. Discussion

4.1. Capricornis Phylogeny and Their Divergences

Capricornis have been morphologically classified into seven species, while only four species have been defined at the molecular level [6,7,8,9]. In this study, with the adding of the ancient Capricornis individuals, we observed four Capricornis mitochondrial clades that agree with previous molecular investigation of four species by Mori et al. [9] (Figure 2 and Figure S4). However, considering that different species or populations could be mixed in the mitochondrial phylogeny in other mammals, such as the African forest elephant and the straight-tusked elephant, the extinct Eurasian cave hyena and the African spotted hyena [33,34], we suspect that the exploration of nuclear genomes of different Capricornis groups may provide better phylogenetic distinctions between serows.
Geist (1987) suggested that C. crispus and C. swinhoei are likely the product of adaptive radiation of mainland serows [35]. It is possible that they were transported to Japan and Taiwan island through a medium such as a land bridge, and subsequently became geographically isolated on the islands [36,37,38]. A previous molecular study based on partial mitochondrial sequence has estimated that the first separation occurred between 0.67 Ma and 1.07 Ma, when C. crispus diverged from the ancestor group [10]. A recent study by Dou et al. utilizing complete mitochondrial genome estimated that Capricornis and Naemorhedus shared a common ancestor approximately 3.79 Ma, and the divergence times of C. crispus and C. swinhoei were estimated around 1.88 Ma and 1.27 Ma, respectively [11]. In this study, we estimated earlier time scales than those noted in Dou et al., with regard to both the Capricornis origin and its divergences. Specifically, we estimated that Capricornis and Naemorhedus had the most recent common ancestor approximately 4.83 Ma, and the divergence times of C. crispus and C. swinhoei were calculated approximately 2.41 Ma and 1.73 Ma, respectively (Figure 2). We suppose that the inconsistent results may be caused by different fossil calibration node selections. We used the updated fossil nodes of bovidae, i.e., 9.8 Ma for the TMRCA of caprini, and 4.91 Ma between Capricornis and Naemorhedus, which were considered to be exhibiting a smaller variance in estimated ages and thus to better reflect the fossil record [28], while Dou et al. chose a fossil calibration node of 4.1 Ma between O. moschatus and C. sumatraensis [11]. Furthermore, we found that the divergence times of these two species that we obtained are close to the times when Japan and Taiwan island became disconnected from the mainland. During the early Pleistocene, land bridges may have existed between the Eurasian continent and the Japanese archipelago due to geological and glacial activities, allowing terrestrial biota to migrate between these two areas [37]. It has been reported that Taiwan island became separated from the mainland in the late Pleistocene due to rising sea levels [36,38,39]. Therefore, we believe it is reasonable that the differentiation of C. crispus and C. swinhoei occurred after serows arrived in Japan and Taiwan island.

4.2. Phylogeography of C. sumatraensis

In both MCC and ML trees, there are two subclades within the C. sumatraensis clade (Figure 2 and Figure S4). Subclade A consists of the Sumatran serow from the Eurasian continent, while most individuals in subclade B are from Southwestern China and Northern Myanmar. Both of our two ancient samples (CADG839 and CADG946, dated 9884–10,160 cal BP, and 2361–2703 cal BP) cluster together with modern serow individuals in clade A, suggesting genetic continuity between Chinese ancient and modern serows. Moreover, fossil evidence indicates that Capricornis arrived in Northern China during the late Pleistocene, while the living range of the extant serows generally does not reach beyond 30° N (Figure 1) [1,3]. The phylogeny of CADG839, which is from Beijing, Northern China, is closely related to individuals from Sichuan, Southwestern China (MH155202) and Cambodia (NC020629), indicating that even the regional extirpation of the serow occurred in Northern China, the genetic diversity was not influenced, at least in the Holocene. Besides, we found that CADG946, the ancient individual from Guizhou, and all modern serow individuals from the same area are in different subclades. This means that there may be a regional replacement in Guizhou in terms of maternal lineage. Subclade A, which the ancient individual from Guizhou belongs to, contains both ancient and modern samples, ranging from Northern China to Cambodia, indicating that this subclade was once widely distributed on the mainland at the beginning of the Holocene.
It has been suggested that porcupine experienced local extirpation in Northern China in the mid-Holocene, while populations in the south survived, resulting in its contemporary distribution pattern [40]. Serow and porcupine shared similar habitats and experienced the same distribution changes according to fossil evidence. We therefore suppose that, similar to the strategy that porcupine adopted when it encountered frequent climate fluctuations, Holocene serows may make themselves adapt to environmental changes by moving southward from Northern China, accompanied by inner regional replacement afterwards.

4.3. Maternal Demographic History of the C. sumatraensis

Our BSP analysis reveals that the effective maternal population size of C. sumatraensis steadily increased from ~225 ka to ~160 ka and from ~90 ka to ~50 ka, and then remained stable for the past 50,000 years (Figure 4), indicating that there had been no reduction in its genetic diversity since MIS7, despite the occurrence of multiple glacial and interglacial cycles [41,42]. Moreover, the population growth in C. sumatraensis occurred in two cold periods (MIS6e–MIS6d and MIS5a–MIS4) [43,44], which may be linked to their adaptive ability. This may be due to their adaptive ability to different altitudes [45], such as migrating to lower altitude valleys during cold periods to alleviate survival pressures and gain stronger inter-species competitive advantages.
Notably, in contrast to the increased or stable mitochondrial genetic diversity indicated in our BSP analysis, recent studies suggest that both environmental changes and human activities have led to a decline in individual number in populations of C. sumatraensis [46,47,48]. We can infer that the decrease in individual numbers in its populations either has not affected its mitochondrial genetic diversity, or that there is a temporal latency to show its effect on genetic components. Thus far, it has been difficult to detect what the truth is due to limited molecular investigations of both ancient and modern C. sumatraensis. However, strengthening the protection of this vulnerable species, such as by reducing damage to its habitat and prohibiting illegal hunting, will help to alleviate the risk of turning it into an endangered species in both situations.

5. Conclusions

In summary, we presented the first molecular information on Holocene serows. We obtained four mitochondrial clades corresponding to four species initiated from a previous molecular investigation. We estimated earlier divergence times of different serows than previously suggested and detected connections between our estimation and the geographic isolation of Japan and Taiwan island. We confirmed that there is genetic continuity between Chinese ancient and modern C. sumatraensis individuals, although this species experienced regional extirpation in Northern China. Additionally, we found that there is no decrease in the maternal genetic diversity of the Sumatran serow even when the number of individuals in its population has declined under the pressure of environmental changes and human activities.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes14061187/s1, Figure S1: Photos of two serow sub-fossil specimens. (a): CADG946, (b): CADG839; Figure S2: Bayesian phylogenetic tree of mitochondrial genomes in the Bovidae family. The calibration points are indicated with an asterisk. Node values are divergence times and node bars are the 95% bounds of the highest posterior density (95% HPD); Figure S3: DNA damage signal of two serow sub-fossil specimens (a) Mitochondrial read length distribution of two serow sub-fossil specimens. (b) Mitochondrial DNA damage plots for four two serow sub-fossil specimens. X axis represents position from 5′ (left) and 3′ (right) read ends. Red line corresponds to C to T substitutions and blue line corresponds to G to A substitutions; Figure S4: The phylogenetic tree constructed using the maximum-likelihood (ML) and neighbor-joining (NJ) methods, with support values for the main nodes indicated below the nodes (ML/NJ); Table S1: Information on two serow sub-fossil specimens; Table S2: Number of reads that match different reference genomes of serows; Table S3: Quality information on the sequencing libraries; Table S4: Information on data set used in this study; Table S5: Mitochondrial genome information on serows used in this study. *** means “No information found”.

Author Contributions

Conceptualization, S.S., X.L. and G.S.; methodology, S.S., G.S. and J.Y.; software, B.X., S.S. and G.S.; validation, S.S., J.H., Z.D. and H.L.; formal analysis, G.S. and J.H.; investigation, S.S. and B.X.; resources, K.X., D.P. and X.H.; data curation, S.S.; writing—original draft preparation, S.S., G.S., B.X., J.H. and X.L.; writing–review and editing, S.S. and G.S.; visualization, S.S. and B.X.; supervision, G.S.; project administration, G.S. and J.Y.; funding acquisition, G.S. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by the National Natural Science Foundation of China (grant No. 42172027).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Two ancient mitochondrial genomes have been submitted to the GenBank database under accession number OQ915127 (CADG839) and OQ915128 (CADG946). The data will be released on 28 April 2024. The address is as follows: GenBank www.ncbi.nlm.nih.gov/genbank, accessed on 28 April 2023.

Acknowledgments

We warmly thank Haowen Tong for his contribution to the sampling collection.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Phan, T.D.; Nijhawan, S.; Li, S. Capricornis sumatraensis. The IUCN Red List of Threatened Species 2020: E.T162916735A162916910. Available online: https://www.iucnredlist.org/species/162916735/162916910 (accessed on 20 April 2023).
  2. Schaller, G.B. Mountain Monarchs: Wild Sheep and Goats of the Himalaya; The University of Chicago Press: Chicago, IL, USA, 1977. [Google Scholar]
  3. Huang, W.B.; Zhu, X.W.; Wang, X.Y. A fossil Ailuropoda and Capricornis from Guilin County, Kwangsi. Vert. PalAsiat 1983, 21, 151–159. [Google Scholar] [CrossRef]
  4. Jass, C.N.; Mead, J.I.  Capricornis crispus. Mamm. Species 2004, 750, 1–10. [Google Scholar] [CrossRef] [PubMed]
  5. Shepherd, C.R. Capricornis rubidus (Amended Version of 2021 Assessment). The IUCN Red List of Threatened Species 2022: E.T3815A214430673. Available online: https://www.iucnredlist.org/species/3815/214430673 (accessed on 20 April 2023).
  6. Grubb, P. Artiodactyla: Bovidae: Caprinae. In Mammal Species of the World: A Taxonomic and Geographic Reference; Wilson, D.E., Reeder, D.M., Eds.; Johns Hopkins University Press: Baltimore, MD, USA, 2005; Volume 1, pp. 703–706. [Google Scholar]
  7. Groves, C.P.; Grubb, P.J. Reclassification of the Serows and Gorals (Nemorhaedus: Bovidae); Croom Helm Editions: London, UK, 1985; pp. 45–50. [Google Scholar]
  8. Valdez, R. Subfamily Caprinae. In Handbook of the Mammals of the World: Vol. 2: Hoofed Mammals; Wilson, D.E., Mittermeier, R.A., Eds.; Lynx Editions: Barcelona, Spain, 2011; pp. 742–759. [Google Scholar]
  9. Mori, E.; Nerva, L.; Lovari, S. Reclassification of the serows and gorals: The end of a neverending story? Mammal Rev. 2019, 49, 256–262. [Google Scholar] [CrossRef]
  10. Liu, W.; Yao, Y.F.; Yu, Q.; Ni, Q.Y.; Zhang, M.W.; Yang, J.D.; Mai, M.M.; Xu, H.L. Genetic variation and phylogenetic relationship between three serow species of the genus Capricornis based on the complete mitochondrial DNA control region sequences. Mol. Biol. Rep. 2013, 40, 6793–6802. [Google Scholar] [CrossRef] [PubMed]
  11. Dou, H.L.; Zhang, Y.S.; Feng, L.M. Complete mitochondrial genome of the Himalayan serow (Capricornis thar) and its phylogenetic status within the genus Capricornis. Biochem. Syst. Ecol. 2016, 65, 115–123. [Google Scholar] [CrossRef]
  12. Tokida, K. Capricornis crispus. The IUCN Red List of Threatened Species 2020: E. T3811A22151909. Available online: https://www.iucnredlist.org/species/3811/22151909 (accessed on 20 April 2023).
  13. Basler, N.; Xenikoudakis, G.; Westbury, M.V.; Song, L.F.; Sheng, G.L.; Barlow, A. Reduction of the contaminant fraction of DNA obtained from an ancient giant panda bone. BMC Res. Notes. 2017, 10, 754. [Google Scholar] [CrossRef]
  14. Meyer, M.; Kircher, M. Illumina sequencing library preparation for highly multiplexed target capture and sequencing. Cold Spring Harb. Protoc. 2010, 2010, pdb.prot5448. [Google Scholar] [CrossRef]
  15. Chen, S.F.; Zhou, Y.Q.; Chen, Y.R.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef]
  16. Li, H.; Durbin, R. Fast and accurate long-read alignment with Burrows–Wheeler transform. Bioinformatics 2010, 26, 589–595. [Google Scholar] [CrossRef]
  17. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef]
  18. Korneliussen, T.S.; Albrechtsen, A.; Nielsen, R. ANGSD: Analysis of Next Generation Sequencing Data. BMC Bioinform. 2014, 15, 356. [Google Scholar] [CrossRef]
  19. Okonechnikov, K.; Conesa, A.; García-Alcalde, F. Qualimap 2: Advanced multi-sample quality control for high-throughput sequencing data. Bioinformatics 2016, 32, 292–294. [Google Scholar] [CrossRef] [PubMed]
  20. Jónsson, H.; Ginolhac, A.; Schubert, M.; Johnson, P.L.; Orlando, L. mapDamage2.0: Fast approximate Bayesian estimates of ancient DNA damage parameters. Bioinformatics 2013, 29, 1682–1684. [Google Scholar] [CrossRef] [PubMed]
  21. Miller, M.A.; Pfeiffer, W.; Schwartz, T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. In Proceedings of the 2010 Gateway Computing Environments Workshop (GCE), New Orleans, LA, USA, 14 November 2010; pp. 1–8. [Google Scholar]
  22. Katoh, K.; Misawa, K.; Kuma, K.; Miyata, T. MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic. Acids. Res. 2002, 30, 3059–3066. [Google Scholar] [CrossRef] [PubMed]
  23. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef] [PubMed]
  24. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods. 2012, 9, 772. [Google Scholar] [CrossRef]
  25. Felsenstein, J. Confidence limits on phylogenies: An approach using the bootstrap. Evolution 1985, 39, 783–791. [Google Scholar] [CrossRef]
  26. Leigh, J.W.; Bryant, D. popart: Full-feature software for haplotype network construction. Methods Ecol. Evol. 2015, 6, 1110–1116. [Google Scholar] [CrossRef]
  27. Drummond, A.J.; Suchard, M.A.; Xie, D.; Rambaut, A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol. Biol. Evol. 2012, 29, 1969–1973. [Google Scholar] [CrossRef]
  28. Bibi, F. A multi-calibrated mitochondrial phylogeny of extant Bovidae (Artiodactyla, Ruminantia) and the importance of the fossil record to systematics. BMC Evol. Biol. 2013, 13, 166. [Google Scholar] [CrossRef]
  29. Pérez, T.; González, I.; Essler, S.E.; Fernández, M.; Domínguez, A. The shared mitochondrial genome of Rupicapra pyrenaica ornata and Rupicapra rupicapra cartusiana: Old remains of a common past. Mol. Phylogenet. Evol. 2014, 79, 375–379. [Google Scholar] [CrossRef] [PubMed]
  30. Luikart, G.; Gielly, L.; Excoffier, L.; Vigne, J.D.; Bouvet, J.; Taberlet, P. Multiple maternal origins and weak phylogeographic structure in domestic goats. Proc. Natl. Acad. Sci. USA 2001, 98, 5927–5932. [Google Scholar] [CrossRef]
  31. Rambaut, A.; Drummond, A.J.; Xie, D.; Baele, G.; Suchard, M.A. Posterior Summarization in Bayesian Phylogenetics Using Tracer 1.7. Syst. Biol. 2018, 67, 901–904. [Google Scholar] [CrossRef]
  32. Sawyer, S.; Krause, J.; Guschanski, K.; Savolainen, V.; Paabo, S. Temporal patterns of nucleotide misincorporations and DNA fragmentation in ancient DNA. PLoS ONE 2012, 7, e34131. [Google Scholar] [CrossRef]
  33. Palkopoulou, E.; Lipson, M.; Mallick, S.; Nielsen, S.; Rohland, N.; Baleka, S.; Karpinski, E.; Ivancevic, A.M.; To, T.; Kortschak, R.D.; et al. A comprehensive genomic history of extinct and living elephants. Proc. Natl. Acad. Sci. USA 2018, 115, E2566–E2574. [Google Scholar] [CrossRef]
  34. Westbury, M.V.; Hartmann, S.; Barlow, A.; Preick, M.; Ridush, B.; Nagel, D.; Rathgeber, T.; Ziegler, R.; Baryshnikov, G.; Sheng, G.; et al. Hyena paleogenomes reveal a complex evolutionary history of cross-continental gene flow between spotted and cave hyena. Sci. Adv. 2020, 6, eaay0456. [Google Scholar] [CrossRef]
  35. Geist, V. On the evolution of the Caprinae. In The Biology and Management of Capricornis and Related Mountain Antelopes; Soma, H., Ed.; Croom Helm Editions: London, UK, 1987; pp. 3–40. [Google Scholar]
  36. Zeng, W.B. The passageway of the flora migration on both sides of the Taiwan strait in Pleistocene epoch. Acta Bot Brasilica 1993, 16, 107–110. [Google Scholar]
  37. Liu, Y.G.; Shi, X.F.; Lv, H.L. Advances in Paleoceanographic Studies on the Japan Sea, Okhotsk Sea and Bering Sea. Adv. Mar. Biol. 2004, 22, 519–530. [Google Scholar]
  38. Zhao, Z.B. A preliminary study on the evolution of Taiwan strait. Taiwan Strait. 1982, 1, 20–24. [Google Scholar]
  39. Jia, C.Z.; He, D.F.; Lu, J.M. Episodes and geodynamic setting of Himalayan movement in China. Oil Gas Geol. 2004, 25, 121–125. [Google Scholar]
  40. Sheng, G.L.; Hu, J.M.; Tong, H.W.; Llamas, B.; Yuan, J.X.; Hou, X.D.; Chen, S.G.; Xiao, B.; Lai, X.L. Ancient DNA of northern China Hystricidae sub-fossils reveals the evolutionary history of old world porcupines in the Late Pleistocene. BMC Evol. Biol. 2020, 20, 88. [Google Scholar] [CrossRef] [PubMed]
  41. Hewitt, G.M. Some genetic consequences of ice ages, and their role in divergence and speciation. Biol. J. Linn. Soc. Lond. 1996, 58, 247–276. [Google Scholar] [CrossRef]
  42. Hewitt, G. The genetic legacy of the Quaternary ice ages. Nature 2000, 405, 907–913. [Google Scholar] [CrossRef]
  43. Weiss, T.L.; Linsley, B.K.; Gordon, A.L.; Rosenthal, Y.; Dannenmann-Di Palma, S. Constraints on Marine Isotope Stage 3 and 5 Sea Level From the Flooding History of the Karimata Strait in Indonesia. Paleoceanogr. Paleoclimatol. 2022, 37, e2021PA004361. [Google Scholar] [CrossRef]
  44. Nowaczyk, N.R.; Liu, J.B.; Plessen, B.; Wegwerth, A.; Arz, H.W. A High-Resolution Paleosecular Variation Record for Marine Isotope Stage 6 From Southeastern Black Sea Sediments. J. Geophys. Res. Solid Earth 2021, 126, e2020JB021350. [Google Scholar] [CrossRef]
  45. Wan, T.; Oaks, J.R.; Jiang, X.L.; Huang, H.; Knowles, L.L. Differences in Quaternary co-divergence reveals community-wide diversification in the mountains of southwest China varied among species. Proc. Biol. Sci. 2021, 288, 20202567. [Google Scholar] [CrossRef]
  46. Shepherd, C.R.; Gomez, L.; Rachmayuningtyas, B.A. The illegal trade of the Sumatran serow Capricornis sumatraensis sumatraensis for traditional medicine in Indonesia. InJAST 2022, 3, 98–104. [Google Scholar] [CrossRef]
  47. Lovari, S.; Mori, E.; Procaccio, E.L. On the Behavioural Biology of the Mainland Serow: A Comparative Study. Animals 2020, 10, 1669. [Google Scholar] [CrossRef]
  48. Isarankura Na Ayudhya, J.; Merceron, G.; Wannaprasert, T.; Jaeger, J.-J.; Chaimanee, Y.; Shoocongdej, R.; Suraprasit, K. Dental mesowear and microwear for the dietary reconstruction of Quaternary Southeast Asian serows and gorals. Front. Ecol. Evol. 2022, 10, 1000168. [Google Scholar] [CrossRef]
Figure 1. The geographical distribution of extant serows (dark grey) and locations of our two ancient samples (red triangles). Serow individuals with geographic information used in subsequent analyses are geographically marked with dots in different colors. Geographical distribution refers to the IUCN Red List [1,5,12]. We used the World Light Gray Base map provided by the ArcGIS platform as the basemap for our study area. The basemap is composed of a series of map tiles. The map tiles were used to provide a geographical reference for our study results, and we extracted a portion of Asia from it.
Figure 1. The geographical distribution of extant serows (dark grey) and locations of our two ancient samples (red triangles). Serow individuals with geographic information used in subsequent analyses are geographically marked with dots in different colors. Geographical distribution refers to the IUCN Red List [1,5,12]. We used the World Light Gray Base map provided by the ArcGIS platform as the basemap for our study area. The basemap is composed of a series of map tiles. The map tiles were used to provide a geographical reference for our study results, and we extracted a portion of Asia from it.
Genes 14 01187 g001
Figure 2. Mitochondrial phylogenetic relationships shown by MCC tree of 20 serows and 3 outgroup species. Two green dots indicate the nodes for calibration used in this study. The labels below the nodes indicate divergence age with blue bars showing 95% highest posterior density. Numbers above the nodes represent the Bayesan posterior probabilities. Our samples’ names are shown in red font. Location information is noted based on previous studies (see Supplementary Table S5).
Figure 2. Mitochondrial phylogenetic relationships shown by MCC tree of 20 serows and 3 outgroup species. Two green dots indicate the nodes for calibration used in this study. The labels below the nodes indicate divergence age with blue bars showing 95% highest posterior density. Numbers above the nodes represent the Bayesan posterior probabilities. Our samples’ names are shown in red font. Location information is noted based on previous studies (see Supplementary Table S5).
Genes 14 01187 g002
Figure 3. The haplogroup network analysis of serows. The number of samples is represented by the size of the circle. Different clades are represented by different colors. The numbers represent the count of mutation sites in the haplotype network, and our samples are shown in red font.
Figure 3. The haplogroup network analysis of serows. The number of samples is represented by the size of the circle. Different clades are represented by different colors. The numbers represent the count of mutation sites in the haplotype network, and our samples are shown in red font.
Genes 14 01187 g003
Figure 4. Bayesian skyline plots reconstructed for C. sumatraensis. The dataset used for BSP analysis consisted of all individuals within the C. sumatraensis (see Supplementary Table S5). Dotted lines show the boundaries of different marine isotope stages (MIS).
Figure 4. Bayesian skyline plots reconstructed for C. sumatraensis. The dataset used for BSP analysis consisted of all individuals within the C. sumatraensis (see Supplementary Table S5). Dotted lines show the boundaries of different marine isotope stages (MIS).
Genes 14 01187 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Song, S.; Xiao, B.; Hu, J.; Lin, H.; Du, Z.; Xiang, K.; Pan, D.; Hou, X.; Yuan, J.; Lai, X.; et al. Ancient Mitogenomes Reveal Stable Genetic Continuity of the Holocene Serows. Genes 2023, 14, 1187. https://doi.org/10.3390/genes14061187

AMA Style

Song S, Xiao B, Hu J, Lin H, Du Z, Xiang K, Pan D, Hou X, Yuan J, Lai X, et al. Ancient Mitogenomes Reveal Stable Genetic Continuity of the Holocene Serows. Genes. 2023; 14(6):1187. https://doi.org/10.3390/genes14061187

Chicago/Turabian Style

Song, Shiwen, Bo Xiao, Jiaming Hu, Haifeng Lin, Zhicheng Du, Kunpeng Xiang, Dong Pan, Xindong Hou, Junxia Yuan, Xulong Lai, and et al. 2023. "Ancient Mitogenomes Reveal Stable Genetic Continuity of the Holocene Serows" Genes 14, no. 6: 1187. https://doi.org/10.3390/genes14061187

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop