Haplotype diversity of the mitochondrial DNA D-loop region in Calomys musculinus (Rodentia,Muridae) detected by PCR-RFLP

In order to contribute to knowledge of colonization patterns in the rodent Calomys musculinus, a natural reservoir of the virus producing Argentine hemorrhagic fever (AHF), we studied the haplotype diversity of the mitochondrial DNA D-loop region in five natural populations from central Argentina. Digestion with eight restriction enzymes (RsaI, MseI, Tsp509I, AluI, AciI, HaeIII, NlaIII, and AseI) revealed polymorphism in the 1300 bp fragment amplified by PCR. Twenty different composite haplotypes were detected. Hierarchical analyses indicated that almost all variation (94%) is contained within local populations. Haplotypes 1 and 2, shared by all populations, were the most frequent. Nonsignificant genetic differentiation was found among populations of the endemic and nonendemic areas of AHF: All locations sampled presented exclusive haplotypes in spite of their geographic proximity, which would support previous observations indicating restricted gene flow among C. musculinus populations.     

Population structure in the American oyster as inferred by nuclear gene genealogies

Multiple haplotypes from each of three nuclear loci were isolated and sequenced from geographic populations of the American oyster, Crassostrea virginica. In tests of alternative phylogeographic hypotheses for this species, nuclear gene genealogies constructed for these haplotypes were compared to one another, to a mitochondrial gene tree, and to patterns of allele frequency variation in nuclear restriction site polymorphisms (RFLPs) and allozymes. Oyster populations from the Atlantic versus the Gulf of Mexico are not reciprocally monophyletic in any of the nuclear gene trees, despite considerable genetic variation and despite large allele frequency differences previously reported in several other genetic assays. If these populations were separated vicariantly in the past, either insufficient time has elapsed for neutral lineage sorting to have achieved monophyly at most nuclear loci, or balancing selection may have inhibited lineage extinction, or secondary gene flow may have moved haplotypes between regions. These and other possibilities are examined in light of available genetic evidence, and it is concluded that no simple explanation can account for the great variety of population genetic patterns across loci displayed by American oysters. Regardless of the source of this heterogeneity, this study provides an empirical demonstration that different sequences of DNA within the same organismal pedigree can have quite different phylogeographic histories.      

Combining geometric morphometrics and pattern recognition to identify interspecific patterns of skull variation: case study in sympatric Argentinian species of the genus Calomys (Rodentia: Cricetidae: Sigmodontinae)

Sympatric species of vesper mice Calomys laucha and Calomys musculinus are difficult to discriminate, especially in natural history collections where they are identified by external body measurements and cranial characteristics. Accurate identification of these two species can be important because only one of them, C. musculinus , is a Junin virus reservoir, the aetiological agent of the Argentine Hemorragic Fever. Research has focused into the development of molecular techniques to unambiguously identify these species. We apply statistical procedures from the field of pattern recognition to three-dimensional geometric morphometric data based on skull landmarks to identify sympatric species C. laucha , C. musculinus and Calomys venustus . Pattern recognition techniques correctly identified the three species without any prior information on specimen identity. By contrast to expectations, C. venustus differed from the other two species mainly on the basis of shape and not by its centroid size. The main sources of difference between C. laucha and C. musculinus were of shape, specifically localized at the landmarks defined by: (1) the sutures between the premaxillaries, the nasals and the frontals; (2) the sutures between the parietals, the frontals and the squamosals; and (3) the suture between the parietals and the interparietal. Nevertheless, allometries dominate the patterns of interspecific variation between these latter species and may partly explain past identification difficulties. Morphological evolution is discussed. The need for objective methods to define phenotypic clusters is highlighted with respect to the need for fast and precise biodiversity assessments.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 94 , 365–378.     

Assessment of concordance among genealogical reconstructions from various mtDNA segments in three species of Pacific salmon (genus Oncorhynchus)

Seven segments of mitochondrial DNA (mtDNA), comprising 97% of the mitochondrial genome, were amplified by polymerase chain reaction (PCR) and examined for restriction site variation using 13 restriction endonucleases in three species of Pacific salmon: pink (Oncorhynchus gorbuscha), chum (O. keta) and sockeye (O. nerka) salmon. The distribution of variability across the seven mtDNA segments differed substantially among species. Little similarity in the distribution of variable restriction sites was found even between the mitochondrial genomes of the even- and odd-year broodlines of pink salmon. Significantly different levels of nucleotide diversity were detected among three groups of genes: six NADH-dehydrogenase genes had the highest; two rRNA genes had the lowest; and a group that included genes for ATPase and cytochrome oxidase subunits, the cytochrome b gene, and the control region had intermediate levels of nucleotide diversity. Genealogies of mtDNA haplotypes were reconstructed for each species, based on the variation in all mtDNA segments. The contributions of variation within different segments to resolution of the genealogical trees were compared within each species. With the exception of sockeye salmon, restriction site data from different genome segments tended to produce rather different trees (and hence rather different genealogies). In the majority of cases, genealogical information in different segments of mitochondrial genome was additive rather than congruent. This finding has a relevance to phylogeographic studies of other organisms and emphasizes the importance of not relying on a limited segment of the mtDNA genome to derive a phylogeographic structure.     

Recent range expansion and low levels of contemporary gene flow in Calomys musculinus: its relationship with the emergence and spread of Argentine haemorrhagic fever

The geographic distribution of haplotype diversity in the rodent Calomys musculinus, sampled from 16 wild populations of Argentina, was analysed on two geographical scales. The species is the natural reservoir of the Junin virus, the etiological agent of the Argentine haemorrhagic fever (AHF). In all, 24 composite haplotypes were recognised in the mtDNA D-loop region. Haplotypes 1 and 2, internal in the network, were the most frequent and were present in almost all populations. The absence of large genetic gaps between widely distributed haplotypes, the existence of exclusive haplotypes in more than 50% of the sampled populations and the absence of isolation by distance at a macrogeographical scale are in support of the hypothesis of a recent range expansion of the populations of the Humid Pampa, with low to moderate current gene flow. The dispersal of this opportunistic species would have been favoured by the explosive increments in density after agriculture was introduced. When only nearby populations within the endemic area of AHF were considered, a pattern of isolation by distance was detected. At present, genetic drift appears to be the main force acting to randomly differentiate C. musculinus populations, which would also lead to random differentiation of Junin virus strains and a reduction in the virulence of the pathogen in 'historic' AHF areas. The knowledge of migration patterns of the reservoir populations facilitates reliable prediction of the potential spread of the human disease.     

The kinetics of oocyte activation and dynamics of polar body formation in Calomys musculinus and Calomys laucha (Rodentia-Sigmodontinae)

The objective of this study was to determine whether Calomys laucha and Calomys musculinus superovulated oocytes undergo parthenogenetic activation following activation stimuli. Cumulus-intact or denuded oocytes were treated with medium containing ethanol (7%), medium containing strontium chloride, or medium alone. They were then incubated for 6-8 h to allow for activation. A group of oocytes was fixed immediately after maturation to serve as a control. The nuclear status of the oocytes was examined after staining with Hoechst 33342, to determine the timing of pronuclear progression from metaphase II to anaphase II or telophase II or to the pronuclear stage. The proportion of oocytes that underwent activation was higher for oocytes treated with ethanol or strontium chloride than in those incubated in medium alone, for the two species studied (p < 0.001). There was little evidence of spontaneous activation occurring in oocytes during the treatments. Most of the activated oocytes contained a single haploid pronucleus, but it was possible to find immediate cleavage and two pronuclei. The different classes of activated oocytes were cultured for 5 days. The type of activating treatment had a marked effect on the ability of the resulting C. musculinus and C. laucha parthenogenetic embryos to develop to the preimplantation stages. Incubation with ethanol produced only 8-cell embryos while the embryos induced with strontium chloride reached the blastocyst stage. This is the first report of parthenogenesis in C. musculinus and C. laucha. The ability of strontium ions to induce matured secondary oocytes to initiate parthenogenesis and obtain further development of Calomys provides opportunities to use Calomys oocytes in vitro and, therefore, to study the genetics, cell biology and virology of development.     

Successful capacitation and homologous fertilization in vitro in Calomys musculinus and Calomys laucha (Rodentia - sigmodontinae)

Small South American rodents of the genus Calomys have been used extensively for virology and ecological research. Previous studies have demonstrated that Calomys musculinus and Calomys laucha have a relatively short oestrous cycle and that superovulation and parthenogenetic activation can be induced. The purpose of this study was to determine the requirements for in vitro manipulation of the male gamete and in vitro fertilization. Two culture media and different concentrations of spermatozoa were tested for their ability to support sperm motility, hyperactivation and the acrosome reaction. The ability of capacitated Calomys spermatozoa to penetrate zona-free hamster eggs was also evaluated. In vitro fertilization was assessed by examining attachment and binding to the zona pellucida, second polar body extrusion, pronucleus formation and the fertilizing sperm tail. The results of the study showed that: (i) Tyrode's albumin lactate pyruvate (TALP) medium was more effective than T6 medium for maintaining sperm motility in vitro; (ii) hyperactivation was achieved with TALP but not with T6; (iii) the acrosome reaction was easily distinguished by light microscopy and depends on time and sperm concentration; (iv) capacitated spermatozoa are able to penetrate zona-free hamster eggs; and (v) superovulated oocytes can be fertilized in vitro. This is the first report of capacitation and in vitro fertilization for Calomys sp. These results provide opportunities to use C. musculinus and C. laucha as new laboratory animals for research into reproductive biology.     

Glacial-induced altitudinal migrations in Armeria (Plumbaginaceae) inferred from patterns of chloroplast DNA haplotype sharing

In contrast to northern European areas where large-scale migrations occurred to recolonize territories after glacial periods, species in southern regions survived and diverged without large geographical displacements. As a result of the importance of orography in much of the southern areas, such displacements must have involved populations ascending or descending mountains. The present study provides support for glacial-induced altitudinal migrations from chloroplast phylogeographic patterns in Armeria (Plumbaginaceae) in southeast Spain. One hundred and five sequences of the trnL-F spacer were obtained from seven species. Fifteen different haplotypes were recognized, their genealogy was inferred, and associations with geography were explored using nested clade analysis. Seven instances were detected in which the same haplotype is shared by two or three species within a particular massif. In all the cases, at least one of the species involved displayed different haplotypes in other areas; in most, the haplotype shared is predominant either in one of the species involved or in the massif. These patterns of haplotype sharing strongly suggest horizontal transfer between species. In one of the massifs (Sierra Nevada) the three species involved in haplotype sharing (A. splendens, A. filicaulis ssp. nevadensis, A. villosa ssp. bernisii) occur at markedly different altitudinal belts. It is argued that altitudinal migrations within the contraction-expansion model provide the best explanation for the current pattern, and that at least in one case it resulted in the formation of a new hybrid taxon, A. filicaulis ssp. nevadensis.     

Phylogeographic inferences from chloroplast DNA: quantifying the effects of mutations in repetitive and non-repetitive sequences

Phylogeographic inference can be a powerful tool in reconstructing species’ evolutionary histories; however, although inferred phylogeographic patterns should depend in part on the underlying types and rates of mutations, the effects of different types of mutations have seldom been quantified. In this study we identified two chloroplast minisatellites in the common reed Phragmites australis, and showed that these are more variable than chloroplast microsatellites. We then recreated parsimony networks of the global phylogeography of P. australis based on data that either included or excluded repetitive sequences (minisatellites and microsatellites), thereby illustrating the influence that these repetitive sequences can have on large‐scale phylogeographic inference. The resulting networks differed in the numbers of mutational steps, degrees of uncertainty, and total numbers of haplotypes. In addition, the suggested ancestor‐descendant relationships among lineages changed substantially depending on whether repetitive sequences were included. We therefore caution against the inclusion of repetitive sequences in large‐scale networks because of their high potential for homoplasy. Nevertheless, we advocate the inclusion of repetitive sequences in other analyses: specifically, we show that the ratio of mutations in repetitive vs. non‐repetitive regions can provide insight into the relative ages of lineages.     

Genetic diversity of Garra rufa Heckel, 1843 (Teleostei: Cyprinidae) in Anatolia

In the present study, mitochondrial DNA polymerase chain reaction-restriction-fragment length polymorphism (PCR-RFLP) and nuclear DNA inter-simple sequence repeat (ISSR) assays were used to assess the phylogenetic and phylogeographic relationships among Garra rufa samples from Anatolia. The complete mtDNA NADH 3/4 dehydrogenase (ND-3/4) gene amplified by PCR was digested with eight restriction enzymes. These enzymes produced 20 composite haplotypes for G. rufa populations. All the mtDNA haplotypes detected were highly diverged from each other and each lineage had a unique genetic profile. The evaluation of mtDNA PCR-RFLP data coupled with geological history of Anatolia indicated a deep genetic divergence among the mtDNA haplotypes of G. rufa populations from drainages of the Mediterranean and Persian Gulf, suggesting an early isolation of Tigris-Euphrates with Orontes river and other rivers draining into the Mediterranean Sea. In general, data from both mtDNA and nDNA were congruent.     

Patterns of intraspecific diversification in the Piriqueta caroliniana complex in southeastern North America and the Bahamas

Abstract.— We use chloroplast DNA (cpDNA) variation and nested clade phylogeographic analyses to infer the historical processes that have contributed to the high level of morphological and ecological diversification present in a group of herbaceous perennials (the Piriqueta caroliniana complex) in North America and the Bahamas. The presence of morphologically distinct and intercompatible varieties (morphotypes) that can be distinguished based on suites of taxonomic characters (e.g., leaf shape, pubescence type, stature) and contrasting habitat affinities (from marshes to dry pinelands and sand scrub) makes this group particularly appropriate for studies of intraspecific diversification. To examine the distribution of haplotypes among populations, we sampled 467 individuals from 55 locations in Florida, Georgia, and the northern Bahamas (Grand Bahama and Abaco) and screened each individual for cpDNA variation using restriction fragment length polymorphism (RFLP) and heteroduplex analyses. We develop a one-step haplotype phylogeny for this group and use the geographic distributions of haplotypes and clades to test specific phylogeographic hypotheses using the methods developed by Templeton and his colleagues (Templeton 1998). In general, the distribution of haplotypes was strongly influenced by limited dispersal distances, with the more recently derived haplotypes having much lower levels of dispersion and lower frequencies in populations than the ancestral haplotypes. The patterns of clade and haplotype dispersion and displacement and the distribution of morphotypes imply at least three cases of long-distance dispersal and one case of historical fragmentation. The historical patterns inferred for populations of Piriqueta are consistent with known biogeographical events, historical vegetation change, and the concordant patterns of multiple Pleistocene refugia that have been observed for a number of other taxa in southeastern North America.     

Mitochondrial DNA phylogeography in northern grasshopper mice (Onychomys leucogaster) — the influence of Quaternary climatic oscillations on population dispersion and divergence

The details of mitochondrial DNA (mtDNA) phylogenetic structure of the northern grasshopper mouse Onychomys leucogaster were examined using populations from a postulated area of endemism that includes three arid regions (Colorado Plateaus, Interior Plains, and Wyoming Basins) in western North America. Fifteen tetra- and heptanucleotide restriction enzymes were used to assay restriction-site variation in a 2150-bp PCR-amplified fragment of mtDNA representing the ND2 and part of the COI gene regions. A total of 18 mtDNA haplotypes were detected. Although overall genetic divergence among these haplotypes was low (average = 1.1%), phylogeographic structuring was apparent. Notably, a clear phylogenetic split separated one group of haplotypes restricted to the Wyoming Basins from all others. This phylogenetic split was further corroborated by examination of nucleotide sequence variation from a 270-bp stretch of the mtDNA cytochrome b gene. Overall geographic and phylogenetic patterns suggest a complex history of geographic structuring and subsequent mixing of populations of grasshopper mice throughout the late Pleistocene. These patterns of variation are evaluated relative to alternative hypotheses about biotic responses to Quaternary climatic oscillations in western North American arid regions.     

Isolation and characterization of microsatellite markers in Calomys musculinus (Muridae,Sigmodontinae, Phyllotini), the natural reservoir of Junin virus

Calomys musculinus, one of the most abundant rodent species in central Argentina, is the reservoir of Junin virus, the aetiological agent of Argentine hemorrhagic fever. We isolated six polymorphic microsatellite loci for microgeographical studies of population structure in this species. Amplification of these loci in 36 individuals from three natural populations revealed five to 14 alleles per locus and expected heterozygosities from 0.426 to 0.868. Cross‐species amplifications suggest that primers designed for these loci may be useful in other closely related species of the tribe Phyllotini, but not in species of other more distantly related tribes of the subfamily Sigmodontinae.     

Evolutionary relationships among eight species of South American phyllotine rodents (Rodentia: Muridae) based on allozymic data

Rodents of the tribe Phyllotini represents one of the main radiations of South American Sigmodontini. Phylogenetic relationships among species of this highly diversified group are poorly known. In this paper we analyse evolutionary relationships among eight phyllotine species belonging to the genera Calomys , Graomys , Phyllotis and Eligmodontia , on the basis of allozymic polymorphisms. Most of the differences among species were in allele frequencies and not of allele class. Neighbour-joining and maximum likelihood methods place P. xanthopygus in the same group as E. typus and G. griseoflavus , in agreement with results obtained by several authors on the basis of morphological characters. Parsimony analysis of 0–9 coded data suggest that the genus Calomys is paraphyletic, but with a low bootstrap support. In the tree based on genetic distance data, the genus also appears as paraphyletic. The maximum likelihood method yields a tree where Calomys is monophyletic, but this phylogeny is supported by only two out of 78 alleles analysed. Calomys hummelincki and C. venustus occupy a basal position among Calomys species. Calomys musculinus and C. lepidus are the most closely related species of the genus, with C. laucha as sister to them. These relationships are strongly supported by bootstrap percentages.     

Inheritance of random amplified polymorphic DNA (RAPD-PCR) markers and their use in population genetic studies of Calomys musculinus (Rodentia,Muridae), the reservoir of Argentine hemorrhagic fever

In order to assess the reliability of RAPD markers in the estimation of the genetic structure of natural populations of the murid rodent Calomys musculinus (reservoir of Junin virus, ethiological agent of Argentine Hemorrhagic Fever), we have analyzed the heritability of RAPD bands in 10 parents and their offspring (33 individuals). Fourteen out of a total of 119 bands obtained were absent in the parental patterns, but consistently amplified in offspring from some families. These bands can be eliminated from analyses. Overall degree of band sharing between individuals, including non-parental bands, correctly grouped members of a family in the same cluster in a UPGMA tree, with a high bootstrap percentage. Results support the usefulness of RAPDs as hereditable markers. One hundred polymorphic RAPD loci were identified in three natural populations of C. musculinus. Mean expected heterozygosity in three natural populations ranged from 0.206 to 0.220. Allele frequency based and phenotype based measures of genetic differentiation among natural populations of C. musculinus gave similar results (Weir and Cockerham's theta = 0.133; Excoffier et al.'s phi = 0.127). These values were considerably higher than those found previously using allozymes as genetic markers, and are compatible with moderate to low levels of gene flow among populations.     

Cryptic clonal lineages and genetic diversity in the loach Misgurnus anguillicaudatus (Teleostei: Cobitidae) inferred from nuclear and mitochondrial DNA analyses

In the loach Misgurnus anguillicaudatue, the asexual lineage, which produces unreduced clonal diploid eggs, has been identified. Among 833 specimens collected from 54 localities in Japan and two localities in China, 82 candidates of other lineage(s) of cryptic clones were screened by examining RFLP (restriction fragment length polymorphism)-PCR haplotypes in the control region of mtDNA. This analysis was performed because triploid loaches arise from the accidental incorporation of the sperm nucleus into unreduced diploid eggs of a clone. The categorization of members belonging to three newly identified lineages (clones 2–4) and the previously identified clonal lineage (clone 1) was verified by evaluating the genetic identity between two or more individuals from each clonal lineage based on RAPD (random amplified polymorphic DNA)-PCR and multilocus DNA fingerprints. We detected 75 haplotypes by observing the nucleotide status at variable sites from the control region of mtDNA. Phylogenic trees constructed from such sequences showed two highly diversified clades, A and B, that were beyond the level common for interspecific genetic differentiation. That result suggests that M. anguillicaudatus in Japan is not a single species entity. Two clone-specific mtDNA sequences were included in clade A, and the loaches with such sequences may be the maternal origin of the clones. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Accurate identification of three cryptic species of rodents of the genus Calomys using RAPD-PCR and mitDNA RFLP markers

Calomys musculinus, Calomys laucha and Calomys venustus are cryptic species with overlapping distribution ranges. C. musculinus is the natural reservoir of Junin virus (Arenaviridae), the etiological agent of Argentine Hemorrhagic Fever. In epidemiological studies it is very important to unequivocally identify the species of individuals collected in the field in order to test virus infection. The purpose of this work was to describe molecular markers allowing a prompt and clear characterization of individuals of the three species. We studied the D-loop region of mitochondrial DNA by restriction fragment length polymorphism (RFLP). This region was amplified by PCR and the product was digested with nine restriction endonucleases (RE). Eight recognize 4 bp restriction sites (Taqα I, Tsp509 I, Aci I, Rsa I, Alu I, Nla III, Hae III and Mse I) and one recognized a 6 bp sequence (Ase I). Two of them (Aci I and Hae III) did not distinguish any of the three species. Alu I did not discriminate between C. musculinus and C. laucha, but clearly distinguished both from C. venustus. Taq I did not distinguish C. laucha from C. venustus, but differentiated both species from C. musculinus. Mse I distinguished the three species, but some of the polymorphisms of C. musculinus are very similar to C. laucha's restriction pattern. The enzyme Nla III distinguished the three species, but it is highly polymorphic within species. The enzymes Tsp509 I, Rsa I and Ase I clearly discriminated the three species, and patterns obtained with the three of them are recommended for reliable identification of individuals collected in the field. The same DNA samples were used to obtain Random Amplified Polymorphic DNAs (RAPDs) patterns. Several bands produced with primers A02, A06, A08, A09, B09 and OPA02 are species specific and could also be used for identification.     

Parallel evolution of facial stripe patterns in the Neolamprologus brichardi/pulcher species complex endemic to Lake Tanganyika

Colour pattern diversity can be due to random processes or to natural or sexual selection. Consequently, similarities in colour patterns are not always correlated with common ancestry, but may result from convergent evolution under shared selection pressures or drift. Neolamprologus brichardi and Neolamprologus pulcher have been described as two distinct species based on differences in the arrangement of two dark bars on the operculum. Our study uses DNA sequences of the mitochondrial control region to show that relatedness of haplotypes disagrees with species assignment based on head colour pattern. This suggests repeated parallel evolution of particular stripe patterns. The complete lack of shared haplotypes between populations of the same or different phenotypes reflects strong philopatric behaviour, possibly induced by the cooperative breeding mode in which offspring remain in their natal territory and serve as helpers until they disperse to nearby territories or take over a breeding position. Concordant phylogeographic patterns between N. brichardi/N. pulcher populations and other rock-dwelling cichlids suggest that the same colonization routes have been taken by sympatric species and that these routes were affected by lake level fluctuations in the past.     

Analysis of MHC class I genes across horse MHC haplotypes

The genomic sequences of 15 horse major histocompatibility complex (MHC) class I genes and a collection of MHC class I homozygous horses of five different haplotypes were used to investigate the genomic structure and polymorphism of the equine MHC. A combination of conserved and locus-specific primers was used to amplify horse MHC class I genes with classical and nonclassical characteristics. Multiple clones from each haplotype identified three to five classical sequences per homozygous animal and two to three nonclassical sequences. Phylogenetic analysis was applied to these sequences, and groups were identified which appear to be allelic series, but some sequences were left ungrouped. Sequences determined from MHC class I heterozygous horses and previously described MHC class I sequences were then added, representing a total of ten horse MHC haplotypes. These results were consistent with those obtained from the MHC homozygous horses alone, and 30 classical sequences were assigned to four previously confirmed loci and three new provisional loci. The nonclassical genes had few alleles and the classical genes had higher levels of allelic polymorphism. Alleles for two classical loci with the expected pattern of polymorphism were found in the majority of haplotypes tested, but alleles at two other commonly detected loci had more variation outside of the hypervariable region than within. Our data indicate that the equine major histocompatibility complex is characterized by variation in the complement of class I genes expressed in different haplotypes in addition to the expected allelic polymorphism within loci.