Genetic footprints of late Quaternary climate change in the diversity of Patagonian-Fueguian rodents

Species are impacted by climate change at both ecological and evolutionary time scales. Studies in northern continents have provided abundant evidence of dramatic shifts in distributions of species subsequent to the last glacial maximum (LGM), particularly at high latitudes. However, little is known about the history of southern continents, especially at high latitudes. South America is the only continent, other than Antarctica, that extends beyond 40 degrees S. Genetic studies of a few Patagonian species have provided seemingly conflicting results, indicating either postglacial colonization from restricted glacial refugia or persistence through glacial cycles and in situ differentiation. Using mitochondrial DNA sequences of 14 species of sigmodontine rodents, a major faunal ensemble of Patagonia and Tierra del Fuego, we show that at least nine of these species bear genetic footprints of demographic expansion from single restricted sources. However, timing of demographic expansion precedes the LGM in most of these species. Four species are fragmented phylogeographically within the region. Our results indicate that (i) demographic instability in response to historical climate change has been widespread in the Patagonian-Fueguian region, and is generally more pronounced at high latitudes in both southern and northern continents; (ii) colonization from lower latitudes is an important component of current Patagonian-Fueguian diversity; but (iii) in situ differentiation has also contributed to species diversity.

Dispersal and population structure at different spatial scales in the subterranean rodent Ctenomys australis

Background

The population genetic structure of subterranean rodent species is strongly affected by demographic (e.g. rates of dispersal and social structure) and stochastic factors (e.g. random genetic drift among subpopulations and habitat fragmentation). In particular, gene flow estimates at different spatial scales are essential to understand genetic differentiation among populations of a species living in a highly fragmented landscape. Ctenomys australis (the sand dune tuco-tuco) is a territorial subterranean rodent that inhabits a relatively secure, permanently sealed burrow system, occurring in sand dune habitats on the coastal landscape in the south-east of Buenos Aires province, Argentina. Currently, this habitat is threatened by urban development and forestry and, therefore, the survival of this endemic species is at risk. Here, we assess population genetic structure and patterns of dispersal among individuals of this species at different spatial scales using 8 polymorphic microsatellite loci. Furthermore, we evaluate the relative importance of sex and habitat configuration in modulating the dispersal patterns at these geographical scales.

Results

Our results show that dispersal in C. australis is not restricted at regional spatial scales (~ 4 km). Assignment tests revealed significant population substructure within the study area, providing support for the presence of two subpopulations from three original sampling sites. Finally, male-biased dispersal was found in the Western side of our study area, but in the Eastern side no apparent philopatric pattern was found, suggesting that in a more continuous habitat males might move longer distances than females.

Conclusions

Overall, the assignment-based approaches were able to detect population substructure at fine geographical scales. Additionally, the maintenance of a significant genetic structure at regional (~ 4 km) and small (less than 1 km) spatial scales despite apparently moderate to high levels of gene flow between local sampling sites could not be explained simply by the linear distance among them. On the whole, our results support the hypothesis that males disperse more frequently than females; however they do not provide support for strict philopatry within females.

Phylogeographical structure in the subterranean tuco-tuco Ctenomys talarum (Rodentia: Ctenomyidae): contrasting the demographic consequences of regional and habitat-specific histories

In this work we examined the phylogeography of the South American subterranean herbivorous rodent Ctenomys talarum (Talas tuco-tuco) using mitochondrial DNA (mtDNA) control region (D-loop) sequences, and we assessed the geographical genetic structure of this species in comparison with that of subterranean Ctenomys australis, which we have shown previously to be parapatric to C. talarum and to also live in a coastal sand dune habitat. A significant apportionment of the genetic variance among regional groups indicated that putative geographical barriers, such as rivers, substantially affected the pattern of genetic structure in C. talarum. Furthermore, genetic differentiation is consistent with a simple model of isolation by distance, possibly evidencing equilibrium between gene flow and local genetic drift. In contrast, C. australis showed limited hierarchical partitioning of genetic variation and departed from an isolation-by-distance pattern. Mismatch distributions and tests of neutrality suggest contrasting histories of these two species: C. talarum appears to be characterized by demographic stability and no significant departures from neutrality, whereas C. australis has undergone a recent demographic expansion and/or departures from strict neutrality in its mtDNA.

Phylogeography of the deer mouse (Peromyscus maniculatus) provides a predictive framework for research on hantaviruses

Phylogeographical partitioning of Sin Nombre and Monongahela viruses (hantaviruses) may reflect that of their primary rodent host, the deer mouse (Peromyscus maniculatus). Lack of a comprehensive assessment of phylogeographical variation of the host has precluded the possibility of predicting spatial limits of existing strains of these viruses or geographical regions where novel viral strains might emerge. The complete cytochrome b gene was sequenced for 206 deer mice collected from sites throughout North America to provide a foundation for future studies of spatial structure and evolution of this ubiquitous host. Bayesian analyses of these sequences partitioned deer mice into six largely allopatric lineages, some of which may represent unrecognized species. The geographical distributions of these lineages were probably shaped by Quaternary climatic events. Populations of mice were apparently restricted to refugia during glacial advances, where they experienced genetic divergence. Expansion of these populations, following climatic amelioration, brought genetically distinctive forms into contact. Occurrence of parallel changes in virus strains can now be explored in appropriate regions. In New Mexico, for example, near the location where Sin Nombre virus was first discovered, there are three genetically distinctive lineages of deer mice whose geographical ranges need to be delineated precisely. The phylogeography of P. maniculatus provides a framework for interpreting geographical variability, not only in hosts, but also in associated viral variants and disease transmission, and an opportunity to predict the potential geographical distribution of newly emerging viral strains.

Adaptive evolution of the insulin gene in caviomorph rodents

Insulin is a conservative molecule among mammals, maintaining both its structure and function. Rodents that belong to the Suborder Hystricognathi represent an exception, having a very divergent molecule with unusual physiological properties. In this work, we analyzed the evolutionary pattern of the insulin gene in caviomorph rodents (South American hystricomorph rodents). We found that these rodents have higher rates of nonsynonymous:synonymous substitutions (d(N)/d(S)) than nonhystricomorph rodents and that values are heterogeneous inside the group. We estimated codons under positive selection, specifically the second binding site (A13 and B17) and others related with hexamerization (B18, B20, and B22). In the monomer structure, all selected sites formed a single patch around the second binding site. In the hexamer structure, these amino acids were grouped into three major patches. In this structure, contacts between B chains involved all selected sites (except B18), and between faces in the center of the molecule, all contacts were among selected sites. While there is no clear hypothesis regarding the cause of this drastic change, experimental evidence does show that this group of rodents has some peculiarities in growth function, and, whether coincidental or not, these changes appeared together with important changes in life-history traits.

Genetic footprints of demographic expansion in North America, but not Amazonia, during the Late Quaternary

The biotic consequences of climate change have attracted considerable attention. In particular, the "refugial debate" centers on the possible retraction of habitats to limited areas that may have served as refuges for many associated species, especially during glaciations of the Quaternary. One prediction of such scenarios is that populations must have experienced substantial growth accompanying climatic amelioration and the occupation of newly expanded habitats. We used coalescence theory to examine the genetic evidence, or lack thereof, for late Pleistocene refugia of boreal North American and tropical Amazonian mammals. We found substantial and concordant evidence of demographic expansion in North American mammals, particularly at higher latitudes. In contrast, small mammals from western Amazonia appear to have experienced limited or no demographic expansion after the Late Pleistocene. Thus, demographic responses to climate change can be tracked genetically and appear to vary substantially across the latitudinal gradient of biotic diversity.

Genetic and geographic differentiation in the Rio Negro tuco-tuco (Ctenomys rionegrensis): inferring the roles of migration and drift from multiple genetic markers

Among tuco-tucos, Ctenomys rionegrensis is especially amenable to the study of the forces driving population differentiation because of the restricted geographic range it occupies in Uruguay. Within this limited area, the Rio Negro tuco-tuco is limited to sandy soils. It nonetheless exhibits remarkable variation in pelage color, including melanic, agouti, and dark-backed individuals. Two hypotheses have been put forth to explain this pattern: (1) local differentiation and fixation of alternative pelage types by genetic drift under limited gene flow; or (2) fixation by natural selection that may take place even in the presence of gene flow. A previous allozyme study rejected the genetic drift hypothesis on the basis of high inferred levels of migration. New estimates of gene flow from microsatellites and mitochondrial cytochrome b sequences were obtained for C. rionegrensis populations to further test these hypotheses. Much lower levels of gene flow were estimated with these more sensitive markers. Microsatellite-based estimates of gene flow are close to zero and may come closest to estimating current levels of migration. A lack of equilibrium between migration and genetic drift is also strongly suggested by the absence of an isolation-by-distance pattern found in all three genetic datasets. The microsatellite genotype data show that the species is strongly structured geographically, with subpopulations constituting distinct genetic entities. If current levels of gene flow are very low, as indicated by the new data, the local fixation of alternative alleles, including those responsible for pelage color polymorphism, is possible by drift alone. A scenario is thus proposed in which the species expanded in the recent past from a more restricted geographic range and has subsequently differentiated in near isolation, with genetic drift possibly playing a primary role in overall genetic differentiation. The local fixation of pelage color types could also be due to drift, but selection on this trait cannot be ruled out without direct analysis.

Molecular evolution of Aldolase A pseudogenes in mice: multiple origins, subsequent duplications, and heterogeneity of evolutionary rates

The Aldolase multigene family comprises three functional genes (A, B, and C) with tissue-specific expression regulated during ontogeny. DGGE analysis and nucleotide sequencing reveal a family of retropseudogenes of type A in species of MUS: Significant variation in rates of evolution of Aldolase A retropseudogenes is apparent. Our analyses demonstrate that (1) multiple events of retrotransposition are needed to account for the diversity of Aldolase A processed pseudogenes found in mice; (2) some of these sequences have undergone further duplication subsequent to the original retrotransposition event; (3) the patterns of nucleotide substitution are broadly comparable with previous estimates; and (4) estimates of rates of divergence for this array of sequences are up to four times higher than those reported in the literature.

Recurrent amplifications and deletions of satellite DNA accompanied chromosomal diversification in South American tuco-tucos (genus Ctenomys, Rodentia: Octodontidae): a phylogenetic approach

We investigated the relationship between satellite copy number and chromosomal evolution in tuco-tucos (genus Ctenomys), a karyotypically diverse clade of rodents. To explore phylogenetic relationships among 23 species and 5 undescribed forms, we sequenced the complete mitochondrial cytochrome b genes of 27 specimens and incorporated 27 previously published sequences. We then used quantitative dot-blot techniques to assess changes in the copy number of the major Ctenomys satellite DNA (satDNA), named RPCS. Our analysis of the relationship between variation in copy number of RPCS and chromosomal changes employed a maximum-likelihood approach to infer the copy number of the satellite RPCS in the ancestors of each clade. We found that amplifications and deletions of RPCS were associated with extensive chromosomal rearrangements even among closely related species. In contrast, RPCS copy number stability was observed within clades characterized by chromosomal stability. This example reinforces the suspected role of amplification, deletion, and intragenomic movement of satDNA in promoting extensive chromosomal evolution.

The molecular phylogenetics of tuco-tucos (genus Ctenomys, Rodentia: Octodontidae) suggests an early burst of speciation

Variation in the nucleotide sequence of the entire mitochondrial cytochrome b gene (1140 bp) was examined for 27 individuals representing 13 species of South American rodents of the genera Ctenomys, Octodontomys, Tympanoctomys, and Spalacopus. Representatives of the family Echimyidae, Euryzygomatomys and Mesomys, were used as outgroups to test the monophyly of the Octodontinae and Ctenomyinae. Relationships among species of tuco-tucos (genus Ctenomys) were also examined including representatives of the three described subgenera and the two sperm morphs. Reciprocal monophyly of the Octodontinae and Ctenomyinae is strongly supported. Several basal relationships among species of the genus Ctenomys are poorly resolved, suggesting the possibility of a hard polytomy due to a rapid and potentially simultaneous radiation early in the history of the genus. In other cases, clades within the Ctenomyinae previously identified on the basis of allozymes, chromosomes, parasites, or skull morphology were supported. Calibrations based on the fossil record suggest that the mitochondrial cytochrome b of these caviomorphs has evolved at a rapid rate, comparable to those proposed for Mus-Rattus, and three to four times higher than ungulate rates.

Screening techniques for detecting allelic variation in DNA sequences

This article reviews four 'DNA screening techniques', namely heteroduplex analysis, single-strand conformational polymorphism (SSCP), denaturing gradient gel electrophoresis (DGGE) and temperature gradient gel electrophoresis (TGGE) as tools for the study of allelic variation in natural populations. The resolving power, advantages, and limitations of each technique are discussed and compared. We also provide some criteria for choosing among techniques and illustrate some practical issues with examples taken primarily from our own laboratory experience.

Rapid surveying of DNA sequence variation in natural populations

DNA sequencing can be costly and time consuming for population studies because of the relative rarity of variation along exons. These problems can be substantially reduced by the use of the polymerase chain reaction on introns using primers from the exon region. These problems can be further reduced by the use of denaturing gradient gel electrophoresis to identify those alleles in need of sequencing.