A novel integrative approach elucidates fine-scale dispersal patchiness in marine populations.

Dispersal is one of the main determining factors of population structure. In the marine habitat, well-connected populations with large numbers of reproducing individuals are common but even so population structure can exist on a small-scale. Variation in dispersal patterns between populations or over time is often associated to geographic distance or changing oceanographic barriers. Consequently, detecting structure and variation in dispersal on a fine-scale within marine populations still remains a challenge. Here we propose and use a novel approach of combining a clustering model, early-life history trait information from fish otoliths, spatial coordinates and genetic markers to detect very fine-scale dispersal patterns. We collected 1573 individuals (946 adults and 627 juveniles) of the black-faced blenny across a small-scale (2 km) coastline as well as at a larger-scale area (<50 kms). A total of 178 single nucleotide polymorphism markers were used to evaluate relatedness patterns within this well-connected population. In our clustering models we categorized SHORT-range dispersers to be potential local recruits based on their high relatedness within and low relatedness towards other spatial clusters. Local retention and/or dispersal of this potential local recruitment varied across the 2 km coastline with higher frequency of SHORT-range dispersers towards the southwest of the area for adults. An inverse pattern was found for juveniles, showing an increase of SHORT-range dispersers towards the northeast. As we rule out selective movement and mortality from one year to the next, this pattern reveals a complex but not full genetic mixing, and variability in coastal circulation is most likely the main driver of this fine-scale chaotic genetic patchiness within this otherwise homogeneous population. When focusing on the patterns within one recruitment season, we found large differences in temperatures (from approx. 17 °C to 25 °C) as well as pelagic larval duration (PLD) for juveniles from the beginning of the season and the end of the season. We were able to detect fine-scale differences in LONG-range juvenile dispersers, representing distant migrants, depending on whether they were born at the beginning of the season with a longer PLD, or at the end of the reproductive season. The ability to detect such fine-scale dispersal patchiness will aid in our understanding of the underlying mechanisms of population structuring and chaotic patchiness in a wide range of species even with high potential dispersal abilities.

Adiponectin regulates contextual fear extinction and intrinsic excitability of dentate gyrus granule neurons through AdipoR2 receptors.

Post-traumatic stress disorder (PTSD) is characterized by exaggerated fear expression and impaired fear extinction. The underlying molecular and cellular mechanisms of PTSD are largely unknown. The current pharmacological and non-pharmacological treatments for PTSD are either ineffective or temporary with high relapse rates. Here we report that adiponectin-deficient mice exhibited normal contextual fear conditioning but displayed slower extinction learning. Infusions of adiponectin into the dentate gyrus (DG) of the hippocampus in fear-conditioned mice facilitated extinction of contextual fear. Whole-cell patch-clamp recordings in brain slices revealed that intrinsic excitability of DG granule neurons was enhanced by adiponectin deficiency and suppressed after treatment with the adiponectin mimetic AdipoRon, which were associated with increased input resistance and hyperpolarized resting membrane potential, respectively. Moreover, deletion of AdipoR2, but not AdipoR1 in the DG, resulted in augmented fear expression and reduced extinction, accompanied by intrinsic hyperexcitability of DG granule neurons. Adiponectin and AdipoRon failed to induce facilitation of fear extinction and elicit inhibition of intrinsic excitability of DG neurons in AdipoR2 knockout mice. These results indicated that adiponectin action via AdipoR2 was both necessary and sufficient for extinction of contextual fear and intrinsic excitability of DG granule neurons, implying that enhancing or dampening DG neuronal excitability may cause resistance to or facilitation of extinction. Therefore, our findings provide a functional link between adiponectin/AdipoR2 activation, DG neuronal excitability and contextual fear extinction, and suggest that targeting adiponectin/AdipoR2 may be used to strengthen extinction-based exposure therapies for PTSD.

Kinship analyses identify fish dispersal events on a temperate coastline.

Connectivity is crucial for the persistence and resilience of marine species, the establishment of networks of marine protected areas and the delineation of fishery management units. In the marine environment, understanding connectivity is still a major challenge, due to the technical difficulties of tracking larvae. Recently, parentage analysis has provided a means to address this question effectively. To be effective, this method requires limited adult movement and extensive sampling of parents, which is often not possible for marine species. An alternative approach that is less sensitive to constraints in parental movement and sampling could be the reconstruction of sibships. Here, we directly measure connectivity and larval dispersal in a temperate marine ecosystem through both analytical approaches. We use data from 178 single nucleotide polymorphism markers to perform parentage and sibship reconstruction of the black-faced blenny (Tripterygion delaisi) from an open coastline in the Mediterranean Sea. Parentage analysis revealed a decrease in dispersal success in the focal area over 1 km distance and approximately 6.5% of the juveniles were identified as self-recruits. Sibship reconstruction analysis found that, in general, full siblings did not recruit together to the same location, and that the largest distance between recruitment locations was much higher (11.5 km) than found for parent-offspring pairs (1.2 km). Direct measurements of dispersal are essential to understanding connectivity patterns in different marine habitats, and show the degree of self-replenishment and sustainability of populations of marine organisms. We demonstrate that sibship reconstruction allows direct measurements of dispersal and family structure in marine species while being more easily applied in those species for which the collection of the parental population is difficult or unfeasible.

SNP development from RNA-seq data in a nonmodel fish: how many individuals are needed for accurate allele frequency prediction?

Single nucleotide polymorphisms (SNPs) are rapidly becoming the marker of choice in population genetics due to a variety of advantages relative to other markers, including higher genomic density, data quality, reproducibility and genotyping efficiency, as well as ease of portability between laboratories. Advances in sequencing technology and methodologies to reduce genomic representation have made the isolation of SNPs feasible for nonmodel organisms. RNA-seq is one such technique for the discovery of SNPs and development of markers for large-scale genotyping. Here, we report the development of 192 validated SNP markers for parentage analysis in Tripterygion delaisi (the black-faced blenny), a small rocky-shore fish from the Mediterranean Sea. RNA-seq data for 15 individual samples were used for SNP discovery by applying a series of selection criteria. Genotypes were then collected from 1599 individuals from the same population with the resulting loci. Differences in heterozygosity and allele frequencies were found between the two data sets. Heterozygosity was lower, on average, in the population sample, and the mean difference between the frequencies of particular alleles in the two data sets was 0.135 ± 0.100. We used bootstrap resampling of the sequence data to predict appropriate sample sizes for SNP discovery. As cDNA library production is time-consuming and expensive, we suggest that using seven individuals for RNA sequencing reduces the probability of discarding highly informative SNP loci, due to lack of observed polymorphism, whereas use of more than 12 samples does not considerably improve prediction of true allele frequencies.

Development and evaluation of 200 novel SNP assays for population genetic studies of westslope cutthroat trout and genetic identification of related taxa.

DNA sequence data were collected and screened for single nucleotide polymorphisms (SNPs) in westslope cutthroat trout (Oncorhynchus clarki lewisi) and also for substitutions that could be used to genetically discriminate rainbow trout (O. mykiss) and cutthroat trout, as well as several cutthroat trout subspecies. In total, 260 expressed sequence tag-derived loci were sequenced and allelic discrimination genotyping assays developed from 217 of the variable sites. Another 50 putative SNPs in westslope cutthroat trout were identified by restriction-site-associated DNA sequencing, and seven of these were developed into assays. Twelve O. mykiss SNP assays that were variable within westslope cutthroat trout and 12 previously published SNP assays were also included in downstream testing. A total of 241 assays were tested on six westslope cutthroat trout populations (N = 32 per population), as well as collections of four other cutthroat trout subspecies and a population of rainbow trout. All assays were evaluated for reliability and deviation from Hardy-Weinberg and linkage equilibria. Poorly performing and duplicate assays were removed from the data set, and the remaining 200 assays were used in tests of population differentiation. The remaining markers easily distinguished the various subspecies tested, as evidenced by mean G(ST) of 0.74. A smaller subset of the markers (N = 86; average G(ST) = 0.40) was useful for distinguishing the six populations of westslope cutthroat trout. This study increases by an order of magnitude the number of genetic markers available for the study of westslope cutthroat trout and closely related taxa and includes many markers in genes (developed from ESTs).

Discovery and characterization of a large number of diagnostic markers to discriminate Oncorhynchus mykiss and O. clarkii.

Hybridization of cutthroat trout and steelhead/rainbow trout is ubiquitous where they are sympatric, either naturally or owing to introductions. The ability to detect hybridization and introgression between the two species would be greatly improved by the development of more diagnostic markers validated across the two species' many phylogenetic lineages. Here, we describe 81 novel genetic markers and associated assays for discriminating the genomes of these sister species. These diagnostic nucleotide polymorphisms were discovered by sequencing of rainbow trout expressed sequence tags (ESTs) in a diverse panel of both cutthroat trout and steelhead/rainbow trout. The resulting markers were validated in a large number of lineages of both species, including all extant subspecies of cutthroat trout and most of the lineages of rainbow trout that are found in natural sympatry with cutthroat trout or used in stocking practices. Most of these markers (79%) distinguish genomic regions for all lineages of the two species, but a small number do not reliably diagnose coastal, westslope and/or other subspecies of cutthroat trout. Surveys of natural populations and hatchery strains of trout and steelhead found rare occurrences of the alternative allele, which may be due to either previous introgression or shared polymorphism. The availability of a large number of genetic markers for distinguishing genomic regions originating in these sister species will allow the detection of both recent and more distant hybridization events, facilitate the study of the evolutionary dynamics of hybridization and provide a powerful set of tools for the conservation and management of both species.

Leptin restores adult hippocampal neurogenesis in a chronic unpredictable stress model of depression and reverses glucocorticoid-induced inhibition of GSK-3ß/ß-catenin signaling.

Stress and glucocorticoid stress hormones inhibit neurogenesis, whereas antidepressants increase neurogenesis and block stress-induced decrease in neurogenesis. Our previous studies have shown that leptin, an adipocyte-derived hormone with antidepressant-like properties, promotes baseline neurogenesis in the adult hippocampus. This study aimed to determine whether leptin is able to restore suppression of neurogenesis in a rat chronic unpredictable stress (CUS) model of depression. Chronic treatment with leptin reversed the CUS-induced reduction of hippocampal neurogenesis and depression-like behaviors. Leptin treatment elicited a delayed long-lasting antidepressant-like effect in the forced swim behavioral despair test, and this effect was blocked by ablation of neurogenesis with X-irradiation. The functional isoform of the leptin receptor, LepRb, and the glucocorticoid receptor (GR) were colocalized in hippocampal neural stem/progenitor cells in vivo and in vitro. Leptin treatment reversed the GR agonist dexamethasone (DEX)-induced reduction of proliferation of cultured neural stem/progenitor cells from adult hippocampus. Further mechanistic analysis revealed that leptin and DEX converged on glycogen synthase kinase-3ß (GSK-3ß) and ß-catenin. While DEX decreased Ser9 phosphorylation and increased Tyr216 phosphorylation of GSK-3ß, leptin increased Ser9 phosphorylation and attenuated the effects of DEX at both Ser9 and Tyr216 phosphorylation sites of GSK-3ß. Moreover, leptin increased total level and nuclear translocation of ß-catenin, a primary substrate of GSK-3ß and a key regulator in controlling hippocampal neural progenitor cell proliferation, and reversed the inhibitory effects of DEX on ß-catenin. Taken together, our results suggest that adult neurogenesis is involved in the delayed long-lasting antidepressant-like behavioral effects of leptin, and leptin treatment counteracts chronic stress and glucocorticoid-induced suppression of hippocampal neurogenesis via activating the GSK-3ß/ß-catenin signaling pathway.

Discovery and characterization of single nucleotide polymorphisms in Chinook salmon, Oncorhynchus tshawytscha.

Molecular population genetics of non-model organisms has been dominated by the use of microsatellite loci over the last two decades. The availability of extensive genomic resources for many species is contributing to a transition to the use of single nucleotide polymorphisms (SNPs) for the study of many natural populations. Here we describe the discovery of a large number of SNPs in Chinook salmon, one of the world's most important fishery species, through large-scale Sanger sequencing of expressed sequence tag (EST) regions. More than 3 Mb of sequence was collected in a survey of variation in almost 132 kb of unique genic regions, from 225 separate ESTs, in a diverse ascertainment panel of 24 salmon. This survey yielded 117 TaqMan (5' nuclease) assays, almost all from separate ESTs, which were validated in population samples from five major stocks of salmon from the three largest basins on the Pacific coast of the contiguous United States: the Sacramento, Klamath and Columbia Rivers. The proportion of these loci that was variable in each of these stocks ranged from 86.3% to 90.6% and the mean minor allele frequency ranged from 0.194 to 0.236. There was substantial differentiation between populations with these markers, with a mean F(ST) estimate of 0.107, and values for individual loci ranging from 0 to 0.592. This substantial polymorphism and population-specific differentiation indicates that these markers will be broadly useful, including for both pedigree reconstruction and genetic stock identification applications.

Loss of polyubiquitin gene Ubb leads to metabolic and sleep abnormalities in mice.

AIMS: Ubiquitin performs essential roles in a myriad of signalling pathways required for cellular function and survival. Recently, we reported that disruption of the stress-inducible ubiquitin-encoding gene Ubb reduces ubiquitin content in the hypothalamus and leads to adult-onset obesity coupled with a loss of arcuate nucleus neurones and disrupted energy homeostasis in mice. Neuropeptides expressed in the hypothalamus control both metabolic and sleep behaviours. In order to demonstrate that the loss of Ubb results in broad hypothalamic abnormalities, we attempted to determine whether metabolic and sleep behaviours were altered in Ubb knockout mice. METHODS: Metabolic rate and energy expenditure were measured in a metabolic chamber, and sleep stage was monitored via electroencephalographic/electromyographic recording. The presence of neurodegeneration and increased reactive gliosis in the hypothalamus were also evaluated. RESULTS: We found that Ubb disruption leads to early-onset reduced activity and metabolic rate. Additionally, we have demonstrated that sleep behaviour is altered and sleep homeostasis is disrupted in Ubb knockout mice. These early metabolic and sleep abnormalities are accompanied by persistent reactive gliosis and the loss of arcuate nucleus neurones, but are independent of neurodegeneration in the lateral hypothalamus. CONCLUSIONS: Ubb knockout mice exhibit phenotypes consistent with hypothalamic dysfunction. Our data also indicate that Ubb is essential for the maintenance of the ubiquitin levels required for proper regulation of metabolic and sleep behaviours in mice.

Differentiating salmon populations at broad and fine geographical scales with microsatellites and single nucleotide polymorphisms.

Single nucleotide polymorphisms (SNPs) are appealing genetic markers due to several beneficial attributes, but uncertainty remains about how many of these bi-allelic markers are necessary to have sufficient power to differentiate populations, a task now generally accomplished with highly polymorphic microsatellite markers. In this study, we tested the utility of 37 SNPs and 13 microsatellites for differentiating 29 broadly distributed populations of Chinook salmon (n = 2783). Information content of all loci was determined by In and G'(ST), and the top 12 markers ranked by In were microsatellites, but the 6 highest, and 7 of the top 10 G'(ST) ranked markers, were SNPs. The mean ratio of random SNPs to random microsatellites ranged from 3.9 to 4.1, but this ratio was consistently reduced when only the most informative loci were included. Individual assignment test accuracy was higher for microsatellites (73.1%) than SNPs (66.6%), and pooling all 50 markers provided the highest accuracy (83.2%). When marker types were combined, as few as 15 of the top ranked loci provided higher assignment accuracy than either microsatellites or SNPs alone. Neighbour-joining dendrograms revealed similar clustering patterns and pairwise tests of population differentiation had nearly identical results with each suite of markers. Statistical tests and simulations indicated that closely related populations were better differentiated by microsatellites than SNPs. Our results indicate that both types of markers are likely to be useful in population genetics studies and that, in some cases, a combination of SNPs and microsatellites may be the most effective suite of loci.

Low genetic variability in the highly endangered mediterranean monk seal.

Genetic variability is an important component in the ability of populations to adapt in the face of environmental change. Here we report the first description of nuclear genetic variability in the only remaining sizable colony of the Mediterranean monk seal (Monachus monachus), located at Cap Blanc (Western Sahara, Mauritania), whose estimated size during the study period (1994-May 1997) was about 320 individuals. We tested 42 microsatellite loci isolated from five pinniped species in a sample of 52 pups. Three loci failed to give any product, and of the remaining 39, only 15 were polymorphic, with a maximum of 3 alleles detected. Three loci appeared to be X-linked. No departures from Hardy-Weinberg equilibrium were detected and no genetic structure was found between the two nursing caves currently occupied by the seals. Several analytical methods show that, as a consequence of a severe bottleneck, the population has suffered a decrease in genetic variability over the last few centuries.

Detection of reduction in population size using data from microsatellite loci.

We demonstrate that the mean ratio of the number of alleles to the range in allele size, which we term M, calculated from a population sample of microsatellite loci, can be used to detect reductions in population size. Using simulations, we show that, for a general class of mutation models, the value of M decreases when a population is reduced in size. The magnitude of the decrease is positively correlated with the severity and duration of the reduction in size. We also find that the rate of recovery of M following a reduction in size is positively correlated with post-reduction population size, but that recovery occurs in both small and large populations. This indicates that M can distinguish between populations that have been recently reduced in size and those which have been small for a long time. We employ M to develop a statistical test for recent reductions in population size that can detect such changes for more than 100 generations with the post-reduction demographic scenarios we examine. We also compute M for a variety of populations and species using microsatellite data collected from the literature. We find that the value of M consistently predicts the reported demographic history for these populations. This method, and others like it, promises to be an important tool for the conservation and management of populations that are in need of intervention or recovery.

Derivation of a simple microsatellite locus from a compound ancestor in the genus Mus.

Many microsatellite loci contain more than a single repeat motif. These compound microsatellites are perhaps most easily explained as arising from simple ancestors. We have uncovered a contrary example in mice of the genus Mus. Sequence analysis of the locus D1MIT29 in most of the members of the genus Mus reveals that this locus is compound in all species except M. musculus, in which it is simple. Moreover, phylogenetic analyses of base substitutions in the non-repetitive flanking region gives trees which are consistent with the previously accepted phylogenetic hypothesis that M. musculus is nested within the subgenus Mus. This confirms that the history of this locus is similar to that of molecular markers previously used in phylogenetic studies of this group and, therefore, demonstrates that the simple state in this lineage is derived from a compound ancestor. We also demonstrate the utility of this type of nuclear sequence variation for phylogeographic studies in the genus Mus. Finally, our sequences reveal homoplasy for size, reemphasizing the danger of using microsatellite size variation alone when the individuals under study are not closely related.

An empirical genetic assessment of the severity of the northern elephant seal population bottleneck.

A bottleneck in population size of a species is often correlated with a sharp reduction in genetic variation. The northern elephant seal (Mirounga angustirostris) has undergone at least one extreme bottleneck, having rebounded from 20-100 individuals a century ago to over 175,000 individuals today. The relative lack of molecular-genetic variation in contemporary populations has been documented, but the extent of variation before the late 19th century remains unknown. We have determined the nucleotide sequence of a 179 base-pair segment of the mitochondrial DNA (mtDNA) control region from seals that lived before, during and after a bottleneck low in 1892. A 'primerless' PCR was used to improve the recovery of information from older samples. Only two mtDNA genotypes were present in all 150+ seals from the 1892 bottleneck on, but we discovered four genotypes in five pre-bottleneck seals. This suggests a much greater amount of mtDNA genotypic variation before this bottleneck, and that the persistence of two genotypes today is a consequence of random lineage sampling. We cannot correlate the loss of mtDNA genotypes with a lowered mean fitness of individuals in the species today. However, we show that the species historically possessed additional genotypes to those present now, and that sampling of ancient DNA could elucidate the genetic consequences of severe reductions in population size.

Social structure of the mound-building mouse Mus spicilegus revealed by genetic analysis with microsatellites.

The Mound-building mouse Mus spicilegus possesses a unique behaviour amongst mice. It constructs large earthen mounds and associated nesting chambers which serve to store food for immature individuals during the winter nesting period. We have used genetic analysis of four autosomal and four X-linked microsatellite loci to determine relationships between individuals inhabiting 40 mounds in Bulgaria. We show that, in almost all cases, individuals in a mound are the product of multiple parentage. We estimate the minimum number of males and female parents contributing offspring to each mound and demonstrate that at least two male and two female parents contribute offspring to a minimum of seven mounds. Analyses of relatedness coefficients and allele sharing values demonstrate that parents of different sibships within mounds are more related than if they had been chosen at random from the population and suggest that it is the female parents that contribute this excess relatedness. These results suggest that the mechanism by which individuals congregate to build mounds is kin-based and that the evolution of mound building and communal nesting in M. spicilegus is due in part to kin selection. This study represents a novel approach to the study of mammalian behavioural ecology. We have used a genetic dataset to construct an outline of social structure in the absence of behavioural data. These inferences can now be used to direct further work on this species.

Homoplasy for size at microsatellite loci in humans and chimpanzees.

Homoplasy (convergence in the size of different alleles) at microsatellite loci was examined by sequencing multiple alleles of two compound microsatellites and single copies of alleles of the same size at two other compound loci in both chimpanzees and humans. At one of the two loci for which multiple alleles were sequenced, extensive homoplasy for size was uncovered both within and between species. At the three loci for which alleles of the same size were examined in the two species, sequencing demonstrated different internal structures. These results confirm theoretical predictions that a certain fraction of mutations at microsatellite loci should produce alleles that are identical in size but differ by a number of mutations. The sequence data reveal a previously unrecognized class of variation at microsatellites and open up the possibility that DNA sequencing may allow the extraction of more information from these loci, thus increasing their power as variable markers for genetic mapping studies. Conversely, the data also indicate that the assumption that alleles of the same size are identical in sequence, which is implicit in several methods of analysis, is violated in some cases. Therefore, caution should be used when employing microsatellites in phylogenetic and other studies in which the individuals being examined are separated by a great number of generations from a common ancestor.

Microsatellite allele frequencies in humans and chimpanzees, with implications for constraints on allele size.

The distributions of allele sizes at eight simple-sequence repeat (SSR) or microsatellite loci in chimpanzees are found and compared with the distributions previously obtained from several human populations. At several loci, the differences in average allele size between chimpanzees and humans are sufficiently small that there might be a constraint on the evolution of average allele size. Furthermore, a model that allows for a bias in the mutation process shows that for some loci a weak bias can account for the observations. Several alleles at one of the loci (Mfd 59) were sequenced. Differences between alleles of different lengths were found to be more complex than previously assumed. An 8-base-pair deletion was present in the nonvariable region of the chimpanzee locus. This locus contains a previously unrecognized repeated region, which is imperfect in humans and perfect in chimpanzees. The apparently greater opportunity for mutation conferred by the two perfect repeat regions in chimpanzees is reflected in the higher variance in repeat number at Mfd 59 in chimpanzees than in humans. These data indicate that interspecific differences in allele length are not always attributable to simple changes in the number of repeats.

Mutational processes of simple-sequence repeat loci in human populations.

Mutational processes of simple sequence repeats (SSRs) in complex genomes are poorly understood. We examined these processes by introducing a two-phase mutation model. In this model, most mutations are single-step changes, but infrequent large jumps in repeat number also occur. We used computer simulations to determine expected values of statistics that reflect frequency distributions of allele size for the two-phase model and two alternatives, the one-step and geometric models. The theoretical expectations for each model were tested by comparison with observed values for 10 SSR loci genotyped in the Sardinian population, whose genetic and demographic histories have been previously reconstructed. The two-phase model provided the best fit to the data for most of these loci in this population. In the analysis we assumed that the loci were neutral and that this population had undergone rapid population growth. Recent observations made for unstable trinucleotide repeats support our suggestion that frequent small changes and rare large changes in repeat number represent two distinct classes of mutation at SSR loci. We genotyped the same 10 loci in Egyptian and sub-Saharan African samples to assess the utility of SSRs for studying the divergence of populations and found that estimates of interpopulation distances from SSRs were similar to those derived from analysis of mitochondrial DNA.

A phylogenetic study of the gibbons (Hylobates) using DNA obtained noninvasively from hair.

Variation in a 252-nucleotide segment of the cytochrome b gene from 26 gibbons is described. DNA was extracted from hair, amplified, and directly sequenced. These sequences represent seven of the nine nominal species and three of the four hylobatid subgenera. Variation was observed at 55 sites, 42 of which are phylogenetically informative. Levels of transitional and transversional divergence between the taxa are similar to those reported for homologous mtDNA sequences in other mammals. Parsimony, maximum likelihood, and bootstrap analyses (1) support some traditional phylogenetic hypotheses (monophyly of the concolor gibbons), (2) suggest previously unrecognized affinities between the lar species group and Hylobates klossi and between H. lar and H. agilis unko, and (3) show that this segment does not contain information sufficient for completely resolving gibbon relationships at the subgeneric level. The study demonstrates the great potential of noninvasive DNA sampling for phylogenetic analyses of mammals.