Postcranial skeleton of Goodwin's brush-tailed mouse (Calomyscus elburzensis Goodwin, 1939) (Rodentia: Calomyscidae): Shape, size, function, and locomotor adaptation.

Goodwin's brush-tailed mouse (Calomyscus elburzensis Goodwin, 1939) is a poorly known small rodent that occupies rocky habitats in Iran, Turkmenistan, Afghanistan, Pakistan, Azerbaijan, and Syria. Herein, a detailed description of the shape, size, and function of the postcranial skeleton of this species is presented for the first time. Trapping was carried out in eastern Iran between the years 2013 and 2015. Skeletal parts of 24 adult male specimens were removed using the papain digestion protocol, and several postcranial morphological characteristics and measurements were examined. We attempted to achieve a morpho-functional characterization of Goodwin's brush-tailed mouse and to match morphological specializations with previous information on the ecology, behavior, and phylogenetic inferences of this rodent. Goodwin's brush-tailed mouse has extended transverse processes and long zygapophyses in the first five caudal vertebrae along with a good innervation of the caudal vertebrae, which has resulted in a well-developed basal musculature of the tail. It has extended forelimb, long ilium, and short post-acetabular part of the innominate bone, loose hip joint with high degree of lateral movement of the hindlimb, and long distal elements of the hindlimb. These features have resulted in fast terrestrial movements in open microhabitats, including climbing and jumping. Although superficial scratching of the ground is observed, the species is incapable of digging burrows. Evaluation of postcranial morphological characteristics and character states further indicated the basal radiation of the genus Calomyscus among other Muroidea. Findings constitute a source of information for morpho-functional and phylogenetic comparisons between Calomyscidae and other mouse-like muroids.

Two new species of Peromyscus (Cricetidae: Neotominae) from the Transverse Volcanic Belt of Mexico.

Specimens of the Peromyscus boylii species group distributed in the western and northeastern montane regions of Michoacán, México, historically have been assigned to P. levipes. Previous studies indicated that these specimens possessed mitochondrial DNA haplotypes that were distinct from both P. levipes and P. kilpatricki, a recently named species in the P. boylii species group from northeastern Michoacán and western Morelos. Herein karyotypic, DNA sequence, and morphological data were analyzed from those populations to evaluate their taxonomic affinity. Karyotypic data indicated that individuals from western Michoacán (Dos Aguas and Aguililla) and from a newly discovered population in northeastern Michoacán (Zinapécuaro) were chromosomally similar to P. carletoni (FN = 68) but distinct from other taxa assigned to the P. boylii species group. Analyses of cranial characteristics indicated that, relative to other species in the P. boylii species group, two morphologically distinct groups were present that corresponded to the Dos Aguas/Aguililla and Zinapécuaro populations, respectively. The latter population, although represented by a small sample size (n = 5 specimens), appeared to exhibit some trenchant morphological distinctions compared with other cryptic species in the P. boylii group. Phylogenetic analyses (parsimony, Bayesian, and likelihood) of DNA sequences obtained from the mitochondrial cytochrome-b gene indicated that although the individuals from Dos Aguas/Aguililla and Zinapécuaro formed a sister group relationship, they formed monophyletic clades that differed genetically (2.54%)-a level approaching that seen between other sister species of Peromyscus. Further, the Dos Aguas/Aguililla and Zinapécuaro clade was more closely aligned with a clade containing representatives of P. carletoni and P. levipes instead of with those from closer geographic proximities (P. kilpatricki) located in eastern Michoacán. Together, these results indicated that these two populations seemingly represent two undescribed species in the P. boylii species group for which we propose the names Peromyscus greenbaumi for populations in western Michoacán (circa Dos Aguas and Aguililla) and Peromyscus ensinki for populations in northeastern Michoacán (circa Zinapécuaro).

Los especímenes del grupo de especies de Peromyscus boylii distribuidos en las regiones montañosas occidentales y el noreste de Michoacán, México, históricamente fueron asignados a P. levipes. Sin embargo, estudios previos han indicado que estos especímenes poseen haplotipos de ADN mitocondrial que son distintos de P. levipes y P. kilpatricki, una especie recientemente nombrada en el grupo de especies P. boylii del noreste de Michoacán y el oeste de Morelos. Los datos cariotípicos indicaron que los individuos del oeste de Michoacán (Dos Aguas y Aguililla) y de una población recién descubierta en el noreste de Michoacán (Zinapécuaro) eran cromosómicamente similares a P. carletoni (FN = 68), pero distintos de otros taxones asignados al grupo de especies P. boylii. Los análisis de las características craneales indicaron que en relación con otras especies del grupo de especies P. boylii, dos grupos morfológicamente distintos estaban presentes, y que correspondían a las poblaciones de Dos Aguas/Aguililla y Zinapécuaro, respectivamente. Aunque representada por un tamaño de muestra pequeño (n = 5 especímenes), esta última población pareció exhibir algunas diferencias morfológicas en comparación con las otras especies crípticas en el grupo P. boylii. Los análisis filogenéticos (parsimonia, inferencia bayesiana y verosimilitud) de secuencias de ADN obtenidas del gen mitocondrial citocromo-b, indicaron que, aunque los individuos de Dos Aguas/Aguililla y Zinapécuaro tienen una relación de grupo hermano, estos forman clados monofiléticos que difieren genéticamente (2.54%), nivel que se acerca al observado entre otras especies hermanas de Peromyscus. Además, el clado de Dos Aguas/Aguililla y Zinapécuaro está más estrechamente alineado con un clado que contenía representantes de P. carletoni y P. levipes en lugar de aquellos de proximidades geográficas más cercanas (P. kilpatricki) ubicados en el este de Michoacán. Estos resultados indicaron que estas dos poblaciones aparentemente representan dos especies no descritas en el grupo de especies P. boylii para las cuales proponemos los nombres Peromyscus greenbaumi para poblaciones en el oeste de Michoacán (hacia Dos Aguas y Aguililla) y Peromyscus ensinki para poblaciones en el noreste de Michoacán (circa Zinapécuaro).

Preserve a Voucher Specimen! The Critical Need for Integrating Natural History Collections in Infectious Disease Studies.

Despite being nearly 10 months into the COVID-19 (coronavirus disease 2019) pandemic, the definitive animal host for SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2), the causal agent of COVID-19, remains unknown. Unfortunately, similar problems exist for other betacoronaviruses, and no vouchered specimens exist to corroborate host species identification for most of these pathogens. This most basic information is critical to the full understanding and mitigation of emerging zoonotic diseases. To overcome this hurdle, we recommend that host-pathogen researchers adopt vouchering practices and collaborate with natural history collections to permanently archive microbiological samples and host specimens. Vouchered specimens and associated samples provide both repeatability and extension to host-pathogen studies, and using them mobilizes a large workforce (i.e., biodiversity scientists) to assist in pandemic preparedness. We review several well-known examples that successfully integrate host-pathogen research with natural history collections (e.g., yellow fever, hantaviruses, helminths). However, vouchering remains an underutilized practice in such studies. Using an online survey, we assessed vouchering practices used by microbiologists (e.g., bacteriologists, parasitologists, virologists) in host-pathogen research. A much greater number of respondents permanently archive microbiological samples than archive host specimens, and less than half of respondents voucher host specimens from which microbiological samples were lethally collected. To foster collaborations between microbiologists and natural history collections, we provide recommendations for integrating vouchering techniques and archiving of microbiological samples into host-pathogen studies. This integrative approach exemplifies the premise underlying One Health initiatives, providing critical infrastructure for addressing related issues ranging from public health to global climate change and the biodiversity crisis.

Structural Variation of the Turtle Mitochondrial Control Region.

The present study describes the most comprehensive comparison of turtle mtD-loop regions to date. The primary structure was compared from DNA sequences accessed from GenBank from 48 species in 13 families of extant turtles, and secondary structures of the mtD-loop region were inferred from thermal stabilities, using the program Mfold, for each superfamiliy of turtles. Both primary and secondary structures were found to be highly variable across the order. The Cryptodira showed conservation in the primary structure at conserved sequence blocks (CSBs), but the Pleurodira displayed limited conservation of primary structural characters, other than the coreTAS, a binding site for the helicase TWINKLE, which was highly conserved in the Central and Right Domains across the order. No secondary structure was associated with a TAS, but an AT-rich fold (secondary structure) near the 3' terminus of the mtD-loop region was detected in all turtle superfamilies. Mapping of character states of structural features of the mtD-loop region revealed that most character states were autapomorphies and inferred a number of homoplasies. The Left Domain of turtles, containing no highly conserved structural elements, likely does not serve a functional role; therefore, the Central Domain in turtles is likely equivalent to the Left Domain of mammals. The AT-rich secondary structural element near the 3' terminus of the mtD-loop region may be conserved across turtles because of a functional role, perhaps containing the Light Strand Promotor, or perhaps interacting with the TWINKLE-coreTAS complex in the Central and Right Domains to regulate mtDNA replication and transcription.

Multiscale landscape genetics of American marten at their southern range periphery.

American marten (Martes americana) are a conservation priority in many forested regions of North America. Populations are fragmented at the southern edge of their distribution due to suboptimal habitat conditions. Facilitating gene flow may improve population resilience through genetic and demographic rescue. We used a multiscale approach to estimate the relationship between genetic connectivity and landscape characteristics among individuals at three scales in the northeastern United States: regional, subregional, and local. We integrated multiple modeling techniques and identified top models based on consensus. Top models were used to parameterize resistance surfaces at each scale, and circuit theory was used to identify potential movement corridors. Regional gene flow was affected by forest cover, elevation, developed land cover, and slope. At subregional and local scales, the effects were site specific and included subsets of temperature, elevation, developed land cover, and slope. Developed land cover significantly affected gene flow at each scale. At finer scales, lack of variance in forest cover may have limited the ability to detect a relationship with gene flow. The effect of slope on gene flow was positive or negative, depending on the site examined. Occupancy probability was a relatively poor predictor, and we caution its use as a proxy for landscape resistance. Our results underscore the importance of replication and multiscale approaches in landscape genetics. Climate warming and landscape conversion may reduce the genetic connectivity of marten populations in the northeastern United States, and represent the primary challenges to marten conservation at the southern periphery of their range.

Further assessment of the Genus Neodon and the description of a new species from Nepal.

Recent molecular systematic studies of arvicoline voles of the genera Neodon, Lasiopodomys, Phaiomys, and Microtus from Central Asia suggest the inclusion of Phaiomys leucurus, Microtus clarkei, and Lasiopodomys fuscus into Neodon and moving Neodon juldaschi into Microtus (Blanfordimys). In addition, three new species of Neodon (N. linzhiensis, N. medogensis, and N. nyalamensis) have recently been described from Tibet. Analyses of concatenated mitochondrial (Cytb, COI) and nuclear (Ghr, Rbp3) genes recovered Neodon as a well-supported monophyletic clade including all the recently described and relocated species. Kimura-2-parameter distance between Neodon from western Nepal compared to N. sikimensis (K2P = 13.1) and N. irene (K2P = 13.4) was equivalent to genetic distances observed between recognized species of this genus. The specimens sampled from western Nepal were recovered sister to N. sikimensis in the concatenated analysis. However, analyses conducted exclusively with mitochondrial loci did not support this relationship. The occlusal patterns of the first lower (m1) and third upper (M3) molars were simpler in specimens from western Nepal in comparison to N. sikimensis from eastern Nepal and India. Twelve craniodental characters and four external field measurements were examined from specimens of N. sikimensis from eastern Nepal and India, N. irene, and Neodon from western Nepal. Neodon from western Nepal were significantly different from N. sikimensis from eastern Nepal and India in ten out of 16 characters measured and from N. irene for all characters except ear height. Specimens from western Nepal were smaller in size than N. sikimensis from Eastern Nepal and India and larger than N. irene. Together the results of the molecular and morphological analyses indicate that Neodon from western Nepal are distinct under the phylogenetic, genetic and morpho species concepts.

Landscape connectivity for bobcat (Lynx rufus) and lynx (Lynx canadensis) in the Northeastern United States.

Landscape connectivity is integral to the persistence of metapopulations of wide ranging carnivores and other terrestrial species. The objectives of this research were to investigate the landscape characteristics essential to use of areas by lynx and bobcats in northern New England, map a habitat availability model for each species, and explore connectivity across areas of the region likely to experience future development pressure. A Mahalanobis distance analysis was conducted on location data collected between 2005 and 2010 from 16 bobcats in western Vermont and 31 lynx in northern Maine to determine which variables were most consistent across all locations for each species using three scales based on average 1) local (15 minute) movement, 2) linear distance between daily locations, and 3) female home range size. The bobcat model providing the widest separation between used locations and random study area locations suggests that they cue into landscape features such as edge, availability of cover, and development density at different scales. The lynx model with the widest separation between random and used locations contained five variables including natural habitat, cover, and elevation-all at different scales. Shrub scrub habitat-where lynx's preferred prey is most abundant-was represented at the daily distance moved scale. Cross validation indicated that outliers had little effect on models for either species. A habitat suitability value was calculated for each 30 m2 pixel across Vermont, New Hampshire, and Maine for each species and used to map connectivity between conserved lands within selected areas across the region. Projections of future landscape change illustrated potential impacts of anthropogenic development on areas lynx and bobcat may use, and indicated where connectivity for bobcats and lynx may be lost. These projections provided a guide for conservation of landscape permeability for lynx, bobcat, and species relying on similar habitats in the region.

Tooth Morphogenesis and FGF4 Expression During Development of Molar Tooth in Three Muroid Rodents: Calomyscus elburzensis (Calomyscidae), Mesocricetus auratus (Cricetidae) and Mus musculus (Muridae).

To date, no studies have examined the tooth formation during developmental stages of brush-tailed mice (Calomyscidae) and true hamsters (Cricetidae). Herein, we compared the timing of tooth morphogenesis and FGF4 expression pattern during development of the first lower molar in Goodwin's brush-tailed mouse, Calomyscus elburzensis with two other muroid rodents; the house mouse, Mus musculus (Muridae), model organism for tooth morphogenesis, and the golden hamster, Mesocricetus auratus which shares great similarities in cusp pattern with brush-tailed mice. All three species were bred in captivity and developing embryos were isolated at different embryonic days (E). Histological evaluation of lower molars was performed and spatiotemporal pattern of FGF4 expression was determined by immunohistochemistry. Results indicated that morphogenesis of the tooth cusps starts at the beginning of the cap stage of the first lower molar (E14 in house mouse, about E11.5 in golden hamster and E22 in Goodwin's brush-tailed mouse). During the cap to bell stage (E15 in house mouse, E12 in golden hamster and at about E24 in Goodwin's brush-tailed mouse), a decrease in the expression of FGF4 was observed in the mesenchyme, except for the cusp tips. According to our observations, the developmental process of the first lower molar formation in Goodwin's brush-tailed mouse began much later as compared with the other two species. Despite the differences in the temporal pattern of molar development between these three members of the same superfamily (Muroidea), the correlation in the expression of FGF4 with specific stages of tooth morphogenesis supported its regulatory function. Anat Rec, 300:2138-2149, 2017. © 2017 Wiley Periodicals, Inc.

Hypothesis testing clarifies the systematics of the main Central American Chagas disease vector, Triatoma dimidiata (Latreille, 1811), across its geographic range.

The widespread and diverse Triatoma dimidiata is the kissing bug species most important for Chagas disease transmission in Central America and a secondary vector in Mexico and northern South America. Its diversity may contribute to different Chagas disease prevalence in different localities and has led to conflicting systematic hypotheses describing various populations as subspecies or cryptic species. To resolve these conflicting hypotheses, we sequenced a nuclear (internal transcribed spacer 2, ITS-2) and mitochondrial gene (cytochrome b) from an extensive sampling of T. dimidiata across its geographic range. We evaluated the congruence of ITS-2 and cyt b phylogenies and tested the support for the previously proposed subspecies (inferred from ITS-2) by: (1) overlaying the ITS-2 subspecies assignments on a cyt b tree and, (2) assessing the statistical support for a cyt b topology constrained by the subspecies hypothesis. Unconstrained phylogenies inferred from ITS-2 and cyt b are congruent and reveal three clades including two putative cryptic species in addition to T. dimidiata sensu stricto. Neither the cyt b phylogeny nor hypothesis testing support the proposed subspecies inferred from ITS-2. Additionally, the two cryptic species are supported by phylogenies inferred from mitochondrially-encoded genes cytochrome c oxidase I and NADH dehydrogenase 4. In summary, our results reveal two cryptic species. Phylogenetic relationships indicate T. dimidiata sensu stricto is not subdivided into monophyletic clades consistent with subspecies. Based on increased support by hypothesis testing, we propose an updated systematic hypothesis for T. dimidiata based on extensive taxon sampling and analysis of both mitochondrial and nuclear genes.

Mammals from 'down under': a multi-gene species-level phylogeny of marsupial mammals (Mammalia, Metatheria).

Marsupials or metatherians are a group of mammals that are distinct in giving birth to young at early stages of development and in having a prolonged investment in lactation. The group consists of nearly 350 extant species, including kangaroos, koala, possums, and their relatives. Marsupials are an old lineage thought to have diverged from early therian mammals some 160 million years ago in the Jurassic, and have a remarkable evolutionary and biogeographical history, with extant species restricted to the Americas, mostly South America, and to Australasia. Although the group has been the subject of decades of phylogenetic research, the marsupial tree of life remains controversial, with most studies focusing on only a fraction of the species diversity within the infraclass. Here we present the first Methaterian species-level phylogeny to include 80% of the extant marsupial species and five nuclear and five mitochondrial markers obtained from Genbank and a recently published retroposon matrix. Our primary goal is to provide a summary phylogeny that will serve as a tool for comparative research. We evaluate the extent to which the phylogeny recovers current phylogenetic knowledge based on the recovery of "benchmark clades" from prior studies-unambiguously supported key clades and undisputed traditional taxonomic groups. The Bayesian phylogenetic analyses recovered nearly all benchmark clades but failed to find support for the suborder Phalagiformes. The most significant difference with previous published topologies is the support for Australidelphia as a group containing Microbiotheriidae, nested within American marsupials. However, a likelihood ratio test shows that alternative topologies with monophyletic Australidelphia and Ameridelphia are not significantly different than the preferred tree. Although further data are needed to solidify understanding of Methateria phylogeny, the new phylogenetic hypothesis provided here offers a well resolved and detailed tool for comparative analyses, covering the majority of the known species richness of the group.

Towards a phylogenetic approach to the composition of species complexes in the North and Central American Triatoma, vectors of Chagas disease.

Phylogenetic relationships of insect vectors of parasitic diseases are important for understanding the evolution of epidemiologically relevant traits, and may be useful in vector control. The sub-family Triatominae (Hemiptera:Reduviidae) includes ∼140 extant species arranged in five tribes comprised of 15 genera. The genus Triatoma is the most species-rich and contains important vectors of Trypanosoma cruzi, the causative agent of Chagas disease. Triatoma species were grouped into complexes originally by morphology and more recently with the addition of information from molecular phylogenetics (the four-complex hypothesis); however, without a strict adherence to monophyly. To date, the validity of proposed species complexes has not been tested by statistical tests of topology. The goal of this study was to clarify the systematics of 19 Triatoma species from North and Central America. We inferred their evolutionary relatedness using two independent data sets: the complete nuclear internal transcribed spacer-2 ribosomal DNA (ITS-2 rDNA) and head morphometrics. In addition, we used the Shimodaira-Hasegawa statistical test of topology to assess the fit of the data to a set of competing systematic hypotheses (topologies). An unconstrained topology inferred from the ITS-2 data was compared to topologies constrained based on the four-complex hypothesis or one inferred from our morphometry results. The unconstrained topology represents a statistically significant better fit of the molecular data than either the four-complex or the morphometric topology. We propose an update to the composition of species complexes in the North and Central American Triatoma, based on a phylogeny inferred from ITS-2 as a first step towards updating the phylogeny of the complexes based on monophyly and statistical tests of topologies.

Multilocus systematics and non-punctuated evolution of Holarctic Myodini (Rodentia: Arvicolinae).

The tribe Myodini consists of five genera of forest and alpine voles (Alticola, Caryomys, Eothenomys, Hyperacrius and Myodes) distributed throughout the Holarctic. Because mitochondrial evidence has revealed paraphyly and polyphyly among genera, we apply the first multilocus tests to clarify taxonomy and phylogenetic relationships. Our analyses of 28 of 36 species within Myodini, including three not previously sequenced (A. montosa, A. albicaudus, and H. fertilis), identify four distinct clades and provide the first molecular evidence that Hyperacrius may not belong in Myodini. Myodes is paraphyletic, while polyphyly of Alticola reflects apparent ancient mitochondrial introgression. Diversification in this tribe was hypothesized to be tightly linked to Late Cenozoic climatic events, however, lineage through time analysis indicates diversification over the last 4 My was gradual and not strongly punctuated.

Phylogeographic study of Apodemus ilex (Rodentia: Muridae) in Southwest China.

BACKGROUND: The Mountains of southwest China have complex river systems and a profoundly complex topography and are among the most important biodiversity hotspots in the world. However, only a few studies have shed light on how the mountains and river valleys promote genetic diversity. Apodemus ilex is a fine model for investigating this subject. METHODOLOGY/PRINCIPAL FINDINGS: To assess the genetic diversity and biogeographic patterns of Apodemus ilex, the complete cytochrome b gene sequences (1,140 bp) were determined from 203 samples of A. draco/ilex that were collected from southwest China. The results obtained suggested that A. ilex and A. draco are sistergroups and diverged from each other approximately 2.25 million years ago. A. ilex could be divided into Eastern and Western phylogroups, each containing two sub-groups and being widespread in different geographical regions of the southern Hengduan Mountains and the western Yunnan - Guizhou Plateau. The population expansions of A. ilex were roughly from 0.089 Mya to 0.023 Mya. CONCLUSIONS: Our result suggested that A. ilex is a valid species rather than synonym of A. draco. As a middle-high elevation inhabitant, the phylogenetic pattern of A. ilex was strongly related to the complex geographical structures in southwest China, particularly the existence of deep river valley systems, such as the Mekong and Salween rivers. Also, it appears that the evolutionary history of A. ilex, such as lineage divergences and population expansions were strongly affected by climate fluctuation in the Late Pleistocene.

Landscape genetics of fishers (Martes pennanti) in the Northeast: dispersal barriers and historical influences.

Habitat fragmentation and overtrapping are thought to have resulted in severe population declines for fisher (Martes pennanti) across the northeastern United States, and by the end of the 1930s only 3 remnant populations remained. Subsequent trapping cessation, extensive reintroduction programs, and natural recolonization have helped fishers to reclaim much of their historical range. The degree to which these processes have impacted genetic structure in this species, however, remains unknown. We used 11 microsatellites from tissue samples (n = 432) of fishers to characterize contemporary population structure in light of historical population structure and thus to determine the relative influence of anthropogenic disturbances and natural landscape features in shaping genetic structure of the contemporary population. Our results indicated that 3 well-differentiated contemporary populations are present that correspond well with what would be expected based on their reported history. A course barrier to dispersal appears in the western portion of the study area associated with several lakes including Lake George and Great Sacandaga Lake. Large-scale reintroduction efforts and natural recolonizations have largely had predictable impacts on population structure. An important exception is the substantial impact of the reintroduction of fishers to Vermont.

Agricultural management affects evolutionary processes in a migratory songbird.

Hay harvests have detrimental ecological effects on breeding songbirds, as harvesting results in nest failure. Importantly, whether harvesting also affects evolutionary processes is not known. We explored how hay harvest affected social and genetic mating patterns, and thus, the overall opportunity for sexual selection and evolutionary processes for a ground-nesting songbird, the Savannah sparrow (Passerculus sandwichensis). On an unharvested field, 55% of females were in polygynous associations, and social polygyny was associated with greater rates of extra-pair paternity (EPP). In this treatment, synchrony explained variation in EPP rates, as broods by more synchronous females had more EPP than broods by asynchronous females. In contrast, on a harvested field, simultaneous nest failure caused by haying dramatically decreased the overall incidence of EPP by increasing the occurrence of social monogamy and, apparently, the ability of polygynous males to maintain paternity in their own nests. Despite increased social and genetic monogamy, these haying-mediated changes in mating systems resulted in greater than twofold increase in the opportunity for sexual selection. This effect arose, in part, from a 30% increase in the variance associated with within-pair fertilization success, relative to the unharvested field. This effect was caused by a notable increase (+110%) in variance associated with the quality of social mates following simultaneous nest failure. Because up to 40% of regional habitat is harvested by early June, these data may demonstrate a strong population-level effect on mating systems, sexual selection, and consequently, evolutionary processes.

Tempo and mode in the endocannaboinoid system.

The best-known endocannabinoid ligands, anandamide and 2-AG, signal at least seven receptors and involve ten metabolic enzymes. Genes for the receptors and enzymes were examined for heterogeneities in tempo (relative rate of evolution, RRE) and mode (selection pressure, Ka/Ks) in six organisms with sequenced genomes. BLAST identified orthologs as reciprocal best hits, and nucleotide alignments were performed with ClustalX and MacClade. Two bioinformatics platforms, LiKaKs (a distance-based LWL85 model) and SNAP (a parsimony-based NG86 model) made pairwise comparisons of orthologs in murids (rat and mouse) and primates (human and macaque). Mean RRE of the 18 endocannabinoid genes was significantly greater in murids than primates, whereas mean Ka/Ks did not differ significantly. Next we used FUGE (tree-based maximum-likelihood model) to compute human lineage-specific Ka/Ks calculations for 18 genes, which ranged from 1.11 to 0.00, in rank order from highest to lowest: PTPN22, NAAA, TRPV1, TRPA1, NAPE-PLD, MAGL, PPARgamma, FAAH1, COX2, FAAH2, ABDH4, CB2, GPR55, DAGLbeta, PPARalpha, TRPV4, CB1, DAGLalpha; differences were significant (p < 0.0001). Rat and mouse presented different rank orders (e.g., GPR55 generated the greatest Ka/Ks ratio). The 18 genes were then tested for recent positive selection (within 10,000 yr) using an extended haplotype homozygosity analysis of SNP data from the HapMap database. Significant evidence (p < 0.05) of a positive "selective sweep" was exhibited by PTPN22, TRPV1, NAPE-PLD, and DAGLalpha. In conclusion, the endocannabinoid system is collectively under strong purifying selection, although some genes show evidence of adaptive evolution.

Coevolution between cannabinoid receptors and endocannabinoid ligands.

Genes for receptors and ligands must coevolve to maintain coordinated gene expression and binding affinities. Researchers have debated whether anandamide or 2-arachidonyl glycerol (2-AG) is a more "intrinsic" ligand of cannabinoid receptors. We addressed this debate with a coevolutionary analysis, by examining genes for CB1, CB2, and ten genes that encode ligand metabolic enzymes: abhydrolase domain containing 4 protein, cyclooxygenase 2, diacylglycerol lipase paralogs (DAGLalpha, DAGLbeta), fatty acid amide hydrolase paralogs (FAAH1, FAAH2), monoglyceride lipase, N-acylethanolamine acid amidase, NAPE-selective phospholipase D, and protein tyrosine phosphatase non-receptor type 22. Gene trees (cladograms) of CB1, CB2, and ligand enzymes were obtained by searching for orthologs (tBLASTn) in the genomes of nine phylogenetically diverse species, aligning ortholog sequences with ClustalX, and applying Bayesian analysis (MrBayes). Mirrored cladograms provided evidence of coevolution (i.e., parallel cladogenesis). Next we constructed phylograms of CB1, CB2, and the ten enzymes. Phylogram branch lengths were proportional to three sets of maximum likelihood metrics: all-nucleotide-substitutions and NS/SS ratios (using PAUP()), and Ka/Ks ratios (using FUGE). Spurious correlations in all-nucleotide-substitutions trees (due to phylogenetic bias) and in Ka/Ks ratio trees (due to simplistic modeling) were parsed. Branch lengths from equivalent branches in paired trees were correlated by linear regression. Regression analyses, mirrored cladograms, and phylogenetic profiles produced the same results: close associations between cannabinoid receptors and DAGL enzymes. Therefore we propose that cannabinoid receptors initially coevolved with a fatty acid ester ligand (akin to 2-AG) in ancestral metazoans, and affinity for fatty acid ethanolamide ligands (e.g., AEA) evolved thereafter.

Multiple molecular evidences for a living mammalian fossil.

Laonastes aenigmamus is an enigmatic rodent first described in 2005. Molecular and morphological data suggested that it is the sole representative of a new mammalian family, the Laonastidae, and a member of the Hystricognathi. However, the validity of this family is controversial because fossil-based phylogenetic analyses suggest that Laonastes is a surviving member of the Diatomyidae, a family considered to have been extinct for 11 million years. According to these data, Laonastes and Diatomyidae are the sister clade of extant Ctenodactylidae (i.e., gundies) and do not belong to the Hystricognathi. To solve the phylogenetic position of Laonastes, we conducted a large-scale molecular phylogeny of rodents. The analysis includes representatives of all major rodent taxonomic groups and was based on 5.5 kb of sequence data from four nuclear and two mitochondrial genes. To further validate the obtained results, a short interspersed element insertion analysis including 11 informative loci was also performed. Our molecular data based on sequence and short interspersed element analyses unambiguously placed Laonastes as a sister clade of gundies. All alternative hypotheses were significantly rejected based on Shimodaira-Hasegawa tests, supporting the idea that Laonastes does not belong to the Hystricognathi. Molecular dating analysis also supports an ancient divergence, approximately 44 Mya ago, between Ctenodactylidae and Laonastes. These combined analyses support the hypothesis that Laonastes is indeed a living fossil. Protection of this surviving species would conserve an ancient mammalian family.

A shifted repertoire of endocannabinoid genes in the zebrafish (Danio rerio).

The zebrafish has served as a model organism for developmental biology. Sequencing its genome has expanded zebrafish research into physiology and drug-development testing. Several cannabinoid pharmaceuticals are in development, but expression of endocannabinoid receptors and enzymes remains unknown in this species. We conducted a bioinformatics analysis of the zebrafish genome using 17 human endocannabinoid genes as a reference set. Putative zebrafish orthologs were identified in filtered BLAST searches as reciprocal best hits. Orthology was confirmed by three in silico methods: phylogenetic testing, synteny analysis, and functional mapping. Zebrafish expressed orthologs of cannabinoid receptor 1, transient receptor potential channel vanilloid receptor 4, GPR55 receptor, fatty acid amide hydrolase 1, monoacylglycerol lipase, NAPE-selective phospholipase D, abhydrolase domain-containing protein 4, and diacylglycerol lipase alpha and beta; and paired paralogs of cannabinoid receptor 2, fatty acid amide hydrolase 2, peroxisome proliferator-activated receptor alpha, prostaglandin-endoperoxide synthase 2, and transient receptor potential cation channel subtype A1. Functional mapping suggested the orthologs of transient receptor potential vanilloid receptor 1 and peroxisome proliferator-activated receptor gamma lack specific amino acids critical for cannabinoid ligand binding. No orthologs of N-acylethanolamine acid amidase or protein tyrosine phosphatase, non-receptor type 22 were identified. In conclusion, the zebrafish genome expresses a shifted repertoire of endocannabinoid genes. In vitro analyses are warranted before using zebrafish for cannabinoid development testing.

TOWARD A MOLECULAR PHYLOGENY FOR PEROMYSCUS: EVIDENCE FROM MITOCHONDRIAL CYTOCHROME-b SEQUENCES.

One hundred DNA sequences from the mitochondrial cytochrome-b gene of 44 species of deer mice (Peromyscus (sensu stricto), 1 of Habromys, 1 of Isthmomys, 2 of Megadontomys, and the monotypic genera Neotomodon, Osgoodomys, and Podomys were used to develop a molecular phylogeny for Peromyscus. Phylogenetic analyses (maximum parsimony, maximum likelihood, and Bayesian inference) were conducted to evaluate alternative hypotheses concerning taxonomic arrangements (sensu stricto versus sensu lato) of the genus. In all analyses, monophyletic clades were obtained that corresponded to species groups proposed by previous authors; however, relationships among species groups generally were poorly resolved. The concept of the genus Peromyscus based on molecular data differed significantly from the most current taxonomic arrangement. Maximum-likelihood and Bayesian trees depicted strong support for a clade placing Habromys, Megadontomys, Neotomodon, Osgoodomys, and Podomys within Peromyscus. If Habromys, Megadontomys, Neotomodon, Osgoodomys, and Podomys are regarded as genera, then several species groups within Peromyscus (sensu stricto) should be elevated to generic rank. Isthmomys was associated with the genus Reithrodontomys; in turn this clade was sister to Baiomys, indicating a distant relationship of Isthmomys to Peromyscus. A formal taxonomic revision awaits synthesis of additional sequence data from nuclear markers together with inclusion of available allozymic and karyotypic data.