Molecular Phylogenetics and Taxonomy of the =Subgenus pika (Ochotona, Lagomorpha)

Sequencing and analysis of the complete mitochondrial genome of Ochotona hyperborea from China and its phylogenetic analysis

The complete mitogenome sequence of Ochotona hyperborea was determined using long PCR. The genome was 17,063 bp in length and contained 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, one origin of L strand replication, and one control region. The overall base composition of the heavy strand is A (31.1%), C (28.7%), T (26.3%), and G (13.9%). The base compositions present clearly the A-T skew, which is most obvious in the control region and protein-coding genes. Mitochondrial genome analyses based on MP, ML, NJ, and Bayesian analyses yielded identical phylogenetic trees. This study verifies the evolutionary status of Ochotona hyperborea in Ochotonidae at the molecular level. The mitochondrial genome would be a significant supplement for the Ochotona hyperborea genetic background. The eight Ochotona species formed a monophyletic group with the high bootstrap value (100%) in all examinations.

Intraspecific variation and taxonomy of northern pika Ochotona hyperborea (Mammalia, Lagomorpha)

Rock-dwelling pikas of the genus Ochotona have fragmented montane geographical distribution presumably associated with geological and climatic events of Late Quaternary. Within this genus, the northern pika, O. hyperborea has the widest distribution. Taxonomic composition and subspecific structure of this species has been controversial during the past century. Sixteen available names were combined in different manners into a variable number of weakly overlapping subspecies by previous authors. A comparison of data sets from the entire distribution range of northern pikas is presented here for the first time, including: craniometric, bioacoustic, and multilocus genetic data. We examined 301 skulls, and the mitochondrial COI gene and two nuclear introns in 79 specimens, as well as the structure of alarm calls from 107 individuals. Our results show that the six subspecies within O. hyperborea correspond to the six genetic parapatric lineages. Five of these are hypothesized to be involved in hybridization at the edges of their distribution. Three acoustic races completely correspond with the six genetic lineages: each acoustic race consists of two lineages. Morphometric data do not display any phylogenetic signal in our study. The subspecies from the Mountains of Khabarovsk Territory, north of Amur River, is described de novo as O. h. fedoseevi ssp. n. The type locality of O. hyperborea is defined more accurately. A neotype is designated for O. h. cinereofusca to stabilize use of names of pikas from the Amur region.

Altai pika (Ochotona alpina) alarm calls: individual acoustic variation and the phenomenon of call-synchronous ear folding behavior

Non-hibernating pikas collect winter food reserves and store them in hay piles. Individualization of alarm calls might allow discrimination between colony members and conspecifics trying to steal food items from a colony pile. We investigated vocal posture, vocal tract length, and individual acoustic variation of alarm calls, emitted by wild-living Altai pikas Ochotona alpina toward a researcher. Recording started when a pika started calling and lasted as long as possible. The alarm call series of 442 individual callers from different colonies consisted of discrete short (0.073-0.157 s), high-frequency (7.31-15.46 kHz), and frequency-modulated calls separated by irregular intervals. Analysis of 442 discrete calls, the second of each series, revealed that 44.34% calls lacked nonlinear phenomena, in 7.02% nonlinear phenomena covered less than half of call duration, and in 48.64% nonlinear phenomena covered more than half of call duration. Peak frequencies varied among individuals but always fitted one of three maxima corresponding to the vocal tract resonance frequencies (formants) calculated for an estimated 45-mm oral vocal tract. Discriminant analysis using variables of 8 calls per series of 36 different callers, each from a different colony, correctly assigned over 90% of the calls to individuals. Consequently, Altai pika alarm calls are individualistic and nonlinear phenomena might further increase this acoustic individualization. Additionally, video analysis revealed a call-synchronous, very fast (0.13-0.23 s) folding, depression, and subsequent re-expansion of the pinna confirming an earlier report of this behavior that apparently contributes to protecting the hearing apparatus from damage by the self-generated high-intensity alarm calls.

Multilocus approaches reveal underestimated species diversity and inter-specific gene flow in pikas (Ochotona) from southwestern China

The phylogeny of living pikas (Ochotonidae, Ochotona) remains obscure, and pika species diversity in southwestern China has never been well explored. In this study, 96 tissue samples from 11 valid species in three classified subgenera (Pika, Ochotona and Conothoa) from 23 locations were characterized using multilocus sequences of 7031bp. Two mitochondrial (CYT B and COI) and five nuclear gene segments (RAG1, RAG2, TTN, OXAIL and IL1RAPL1) were sequenced. We analysed evolutionary histories using maximum likelihood (RAxML) and Bayesian analyses (BEAST), and we also used molecular species delimitation analyses (BPP) to explore species diversity. Our study supported O. syrinx (O. huangensis) as a distinct clade from all named subgenera. Relationships among subgenera were not fully resolved, which may be due to a rapid diversification in the middle Miocene (∼13.90Ma). Conflicting gene trees implied mitochondrial introgression from O. cansus to O. curzoniae. We uncovered three cryptic species from Shaanxi, Sichuan and Yunnan with strong support, suggesting an underestimation of species diversity in the "sky-island" mountains of southwest China.

The phylogeny of pikas (Ochotona) inferred from a multilocus coalescent approach

The clarification of the systematics of pikas (genus Ochotona) has been hindered by largely overlapping morphological characters among species and the lack of a comprehensive molecular phylogeny. Here we estimate the first multilocus phylogeny of the genus to date, by analysing 12 nuclear DNA markers (total of 7.5Kb) in 11 species of pikas from the four classified subgenera (Pika, Ochotona, Lagotona and Conothoa) using a multispecies coalescent-based framework. The species-tree confirmed the subgeneric classification by retrieving as monophyletic the subgenera represented here by more than one species. Contrary to previous phylogenies based on mtDNA alone, Lagotona was found to be sister to Pika. Also, support for the monophyly of the alpina group was not strong, thus caution should be used in future analyses of this group. A relaxed molecular clock calibrated using the Ochotonidae-Leporidae divergence resulted in more recent estimates of divergence times relative to previous studies. Strong concordance with inferences based on fossil records was found, suggesting that the initial diversification of the genus took place by the end of late Miocene. Finally, this work sets up methodologies and gathers molecular markers that can be used to extend the understanding of the evolutionary history of the genus.

An intriguing shift occurs in the novel protein phosphatase 1 binding partner, TCTEX1D4: evidence of positive selection in a pika model

T-complex testis expressed protein 1 domain containing 4 (TCTEX1D4) contains the canonical phosphoprotein phosphatase 1 (PPP1) binding motif, composed by the amino acid sequence RVSF. We identified and validated the binding of TCTEX1D4 to PPP1 and demonstrated that indeed this protein is a novel PPP1 interacting protein. Analyses of twenty-one mammalian species available in public databases and seven Lagomorpha sequences obtained in this work showed that the PPP1 binding motif ₉₀RVSF₉₃ is present in all of them and is flanked by a palindromic sequence, PLGS, except in three species of pikas (Ochotona princeps, O. dauurica and O. pusilla). Furthermore, for the Ochotona species an extra glycosylation site, motif ₉₆NLS₉₈, and the loss of the palindromic sequence were observed. Comparison with other lagomorphs suggests that this event happened before the Ochotona radiation. The d_N/d_S for the sequence region comprising the PPP1 binding motif and the flanking palindrome highly supports the hypothesis that for Ochotona species this region has been evolving under positive selection. In addition, mutational screening shows that the ability of pikas TCTEX1D4 to bind to PPP1 is maintained, although the PPP1 binding motif is disrupted, and the N- and C-terminal surrounding residues are also abrogated. These observations suggest pika as an ideal model to study novel PPP1 complexes regulatory mechanisms.

The chromosomes of Tsing-Ling pika, Ochotona huangensis Matschie, 1908 (Lagomorpha, Ochotonidae)

The karyotype of the Tsing-Ling (Huanghe) pika, Ochotona huangensis Matschie, 1908 from the forest habitats of the Qinling Mountains (Shaanxi Province, China) was described for the first time. The chromosome set contains 42 chromosomes (NFa=80). The autosomes are 15 meta-submetacentric pairs and 5 subtelocentric pairs. The X chromosome is a medium sized submetacentric; the Y chromosome is a small sized acrocentric. C-banding revealed a localization of heterochromatin in the pericentromeric regions of all autosomes.

Return to Beringia: parasites reveal cryptic biogeographic history of North American pikas

Traditional concepts of the Bering Land Bridge as a zone of predominantly eastward expansion from Eurasia and a staging area for subsequent colonization of lower latitudes in North America led to early inferences regarding biogeographic histories of North American faunas, many of which remain untested. Here we apply a host-parasite comparative phylogeographical (HPCP) approach to evaluate one such history, by testing competing biogeographic hypotheses for five lineages of host-specific parasites shared by the collared pika (Ochotona collaris) and American pika (Ochotona princeps) of North America. We determine whether the southern host species (O. princeps) was descended from a northern ancestor or vice versa. Three parasite phylogenies revealed patterns consistent with the hypothesis of a southern origin, which is corroborated by four additional parasite lineages restricted to O. princeps. This finding reverses the traditional narrative for the origins of North American pikas and highlights the role of dispersal from temperate North America into Beringia in structuring northern diversity considerably prior to the Holocene. By evaluating multiple parasite lineages simultaneously, the study demonstrates the power of HPCP for resolving complex biogeographic histories that are not revealed by characteristics of the host alone.

DNA barcoding in mammals

DNA barcoding provides an operational framework for mammalian taxonomic identification and cryptic species discovery. Focused effort to build a reference library of genetic data has resulted in the assembly of over 35 K mammalian cytochrome c oxidase subunit I sequences and outlined the scope of mammal-related barcoding projects. Based on the above experience, this chapter recounts three typical methodological pathways involved in mammalian barcoding: routine methods aimed at assembling the reference sequence library from high quality samples, express approaches used to attain cheap and fast taxonomic identifications for applied purposes, and forensic techniques employed when dealing with degraded material. Most of the methods described are applicable to a range of vertebrate taxa outside Mammalia.

DNA barcodes effectively identify the morphologically similar Common Opossum (Didelphis marsupialis) and Virginia Opossum (Didelphis virginiana) from areas of sympatry in Mexico

Two morphologically similar species of opossum from the genus Didelphis-Didelphis virginiana and Didelphis marsupialis-cooccur sympatrically in Mexico. High intraspecific variation complicates their morphological discrimination, under both field and museum conditions. This study aims to evaluate the utility and reliability of using DNA barcodes (short standardized genome fragments used for DNA-based identification) to distinguish these two species. Sequences of the cytochrome c oxidase subunit I (Cox1) mitochondrial gene were obtained from 12 D. marsupialis and 29 D. virginiana individuals and were compared using the neighbor-joining (NJ) algorithm with Kimura's two-parameter (K2P) model of nucleotide substitution. Average K2P distances were 1.56% within D. virginiana and 1.65% in D. marsupialis. Interspecific distances between D. virginiana and D. marsupialis varied from 7.8 to 9.3% and their barcode sequences formed distinct non-overlapping clusters on NJ trees. All sympatric specimens of both species were effectively discriminated, confirming the utility of Cox1 barcoding as a tool for taxonomic identification of these morphologically similar taxa.