Glaciation effects on the phylogeographic structure of Oligoryzomys longicaudatus (Rodentia: Sigmodontinae) in the southern Andes

Climate influences the genetic structure and niche differentiation among populations of the olive field mouse Abrothrix olivacea (Cricetidae: Abrotrichini)

Even when environmental variation over time and space is commonly considered as an important driver of population divergence, few evaluations of intraspecific genetic variation explicitly assess whether observed structure has been caused by or is correlated with landscape heterogeneity. Several phylogeographic studies have characterized the mitochondrial diversity of Abrothrix olivacea, but none has incorporated landscape genetics analyses and ecological niche modeling, leaving a gap in the understanding of the species evolutionary history. Here, these aspects were addressed based on 186 single nucleotide polymorphisms, extracted from sequences of 801 bp of Cytb gene, gathered from 416 individuals collected at 103 localities in Argentina and Chile. Employing multivariate statistical analyses (gPCA, Mantel and Partial Mantel Tests, Procrustes Analysis, and RDA), associations between genetic differences and geographic and climatic distances were evaluated. Presence data was employed to estimate the potential geographic distribution of this species during historical and contemporary climatic scenarios, and to address differences among the climatic niches of their main mitochondrial lineages. The significant influence of landscape features in structuring mitochondrial variability was evidenced at different spatial scales, as well as the role of past climatic dynamics in driving geographic range shifts, mostly associated to Quaternary glaciations. Overall, these results suggest that throughout geographic range gene flow is unevenly influenced by climatic dissimilarity and the geographic distancing, and that studied lineages do not exhibit distributional signals of climatic niche conservatism. Additionally, genetic differentiation occurred by more complex evolutionary processes than mere disruption of gene flow or drift.

Genetic Diversity and Demographic History of the Shaggy Soft-Haired Mouse Abrothrix hirta (Cricetidae; Abrotrichini)

Genetic information on species can inform decision making regarding conservation of biodiversity since the response of organisms to changing environments depend, in part, on their genetic makeup. Territories of central-southern Chile and Argentina have undergone a varying degree of impact during the Quaternary, where the response of local fauna and flora was rather species-specific. Here, we focus on the sigmodontine rodent Abrothrix hirta, distributed from 35° S in Chile and Argentina to northern Tierra del Fuego. Based on 119,226 transcriptome-derived SNP loci from 46 individuals of A. hirta, we described the geographic distribution of the genetic diversity of this species using a maximum likelihood tree, principal component and admixture analyses. We also addressed the demographic history of the main intraspecific lineages of A. hirta using GADMA. We found that A. hirta exhibited four allopatric intraspecific lineages. Three main genetic groups were identified by a Principal Component Analysis and by Ancestry analysis. The demographic history of A. hirta was characterized by recent population stability for populations at the northernmost part of the range, while southern populations experienced a recent population expansion.

Genetic Diversity, Population Structure, and Migration Scenarios of the Marsupial "Monito del Monte" in South-Central Chile

In this study, we quantified the 3 pivotal genetic processes (i.e., genetic diversity, spatial genetic structuring, and migration) necessary for a better biological understanding and management of the singular "living-fossil" and near-threatened mouse opossum marsupial Dromiciops gliroides, the "Monito del Monte," in south-central Chile. We used 11 microsatellite loci to genotype 47 individuals distributed on the mainland and northern Chiloé Island. Allelic richness, observed and expected heterozygosity, inbreeding coefficient, and levels of genetic differentiation were estimated. The genetic structure was assessed based on Bayesian clustering methods. In addition, potential migration scenarios were evaluated based on a coalescent theory framework and Bayesian approach to parameter estimations. Microsatellites revealed moderate to high levels of genetic diversity across sampled localities. Moreover, such molecular markers suggested that at least 2 consistent genetic clusters could be identified along the D. gliroides distribution ("Northern" and "Southern" cluster). However, general levels of genetic differentiation observed among localities and between the 2 genetic clusters were relatively low. Migration analyses showed that the most likely routes of migration of D. gliroides occurred 1) from the Southern cluster to the Northern cluster and 2) from the Mainland to Chiloé Island. Our results could represent critical information for future conservation programs and for a recent proposal about the taxonomic status of this unique mouse opossum marsupial.

Genetic variation of the Chilean endemic long-haired mouse Abrothrix longipilis (Rodentia, Supramyomorpha, Cricetidae) in a geographical and environmental context

Quaternary climate and associated vegetational changes affected the fauna of the Chilean Mediterranean ecosystem. Here we studied the genetic variation of the long-haired mouse, Abrothrix longipilis, a sigmodontine rodent endemic to this area. Within an environmentally explicit context, we examined the geographic distribution of the genetic diversity and demographic history of the species based on sequences of the mitochondrial Cytochrome-b gene of 50 individuals from 13 localities and a large panel of single nucleotide polymorphisms of 17 individuals from 6 localities. The gene genealogy of A. longipilis revealed three intraspecific lineages that are allopatric and latitudinally segregated (northern, central, and southern lineages) with an estimated crown age for the whole species clade of 552.3 kyr B.P. A principal component analysis based on 336,596 SNP loci is in line with the information given by the the mitochondrial gene genealogy. Along its complete distributional range, A. longipilis showed patterns of isolation by distance and also isolation by environment. The general pattern of historical demography showed stability for most intraspecific lineages of A. longipilis. Northern and central lineages showed signals of historical demographic stability, while the southern lineage showed contrasting signals. In agreement with this, the niche models performed showed that in the northern range of A. longipilis, areas of high suitability for this species increased towards the present time; areas of central range would have remained relatively stable, while southern areas would have experienced more change through time. In summary, our study shows three distinct allopatric lineages of A. longipilis, each showing slightly different demographic history.

A 19 Year Analysis of Small Mammals Associated with Human Hantavirus Cases in Chile

Small mammals present in areas where hantavirus cardiopulmonary syndrome (HCPS) cases had occurred in central and southern Chile were captured and analyzed to evaluate the abundance of rodents and seroprevalence rates of antibodies to Andes orthohantavirus (ANDV). Sampling areas ranged from the Coquimbo to Aysén regions (30–45° S approx.) regions. Ninety-two sites in peridomestic and countryside areas were evaluated in 19 years of sampling. An antibody against ANDV was detected by strip immunoassay in 58 of 1847 specimens captured using Sherman traps. Of the eleven species of rodents sampled, Abrothrix olivacea, Oligoryzomys longicaudatus and Abrothrix hirta were the most frequently trapped. O. longicaudatus had the highest seropositivity rate, and by logistic regression analysis, O. longicaudatus of at least 60 g had 80% or higher probability to be seropositive. Sex, age and wounds were significantly related to seropositivity only for O. longicaudatus. Across administrative regions, the highest seropositivity was found in the El Maule region (34.8–36.2° S), and the highest number of HCPS cases was registered in the Aysén region. Our results highlight the importance of long term and geographically extended studies, particularly for highly fluctuating pathogens and their reservoirs, to understand the implications of the dynamics and transmission of zoonotic diseases in human populations.

Genetic variability and structure of the Olive Field Mouse: a sigmodontine rodent in a biodiversity hotspot of southern Chile

The temperate rainforests of southern Chile, a recognized biodiversity hotspot, were significantly affected by Pleistocene glacial cycles in their southern portion and have been severely disrupted mainly due to recent human activities. Additionally, the landscape is characterized by a series of potential barriers to gene flow, such as the Chacao Channel, Cordillera de Piuche in Chiloé and both the Ancud and the Corcovado gulfs. We used mitochondrial DNA sequences and microsatellite data across several populations to evaluate the genetic variability and structure of the sigmodontine rodent Abrothrix olivacea brachiotis, one of the most common species of small mammals and an inhabitant of these biodiverse forests. Sequencing data showed that along with the recovery of high haplotype variation for this species, there was a low nucleotide diversity between haplotypes, showing no genetic differences between the Chiloé Island and continental populations in southern Chile or through any other geographic barrier in the study area. However, microsatellite data exhibited some level of population structuring. The most evident clusterings were those of the Chiloé Island and that of North Patagonia. These findings are corroborated by a barrier analysis that showed a genetic barrier in the latter areas, whereas the Chacao Channel was not a significant barrier for this rodent. Overall, the genetic variability and structure of A. o. brachiotis was concordant with historical factors, such as the Last Glacial Maximum and the presence of geographic elements that isolate populations.

Local persistence of Mann's soft-haired mouse Abrothrix manni (Rodentia, Sigmodontinae) during Quaternary glaciations in southern Chile

Quaternary climatic oscillations have impacted Patagonian sigmodontine fauna, leaving traceable genetic footprints. In southern Chile, changes in the landscape included transitions to different vegetation formations as well as the extension of ice sheets. In this study, we focus on the Valdivian forest endemic and recently described sigmodontine species Abrothrix manni. We aim to assess the genetic structure of this species, testing for the existence of intraspecific lineages, and inferring the recent demographic history of the species. Analyses were based on the first 801 bp of the mitochondrial gene Cytocrhome-b from 49 individuals of A. manni collected at 10 localities that covers most part of its geographic distribution. Genealogical analyses recovered two main intraspecific lineages that are geographically segregated and present an intermediate site of secondary contact. Historical demography shows signal of recent population decrease. Based on these results, we proposed that current genetic diversity of A. manni differentiated in at least two distinct refugial areas in southern Chile. This scenario, in addition to be unique among those uncovered for the so far studied Valdivian forest rodents, is noteworthy because of the reduced geographic scale inhabited by the species.

Does silvoagropecuary landscape fragmentation affect the genetic diversity of the sigmodontine rodent Oligoryzomys longicaudatus?

BACKGROUND

Fragmentation of native forests is a highly visible result of human land-use throughout the world. In this study, we evaluated the effects of landscape fragmentation and matrix features on the genetic diversity and structure of Oligoryzomys longicaudatus, the natural reservoir of Hantavirus in southern South America. We focused our work in the Valdivian Rainforest where human activities have produced strong change of natural habitats, with an important number of human cases of Hantavirus.

METHODS

We sampled specimens of O. longicaudatus from five native forest patches surrounded by silvoagropecuary matrix from Panguipulli, Los Rios Region, Chile. Using the hypervariable domain I (mtDNA), we characterized the genetic diversity and evaluated the effect of fragmentation and landscape matrix on the genetic structure of O. longicaudatus. For the latter, we used three approaches: (i) Isolation by Distance (IBD) as null model, (ii) Least-cost Path (LCP) where genetic distances between patch pairs increase with cost-weighted distances, and (iii) Isolation by Resistance (IBR) where the resistance distance is the average number of steps that is needed to commute between the patches during a random walk.

RESULTS

We found low values of nucleotide diversity (π) for the five patches surveyed, ranging from 0.012 to 0.015, revealing that the 73 sampled specimens of this study belong to two populations but with low values of genetic distance (γST ) ranging from 0.022 to 0.099. Likewise, we found that there are no significant associations between genetic distance and geographic distance for IBD and IBR. However, we found for the LCP approach, a significant positive relationship (r = 0.737, p = 0.05), with shortest least-cost paths traced through native forest and arborescent shrublands.

DISCUSSION

In this work we found that, at this reduced geographical scale, Oligoryzomys longicaudatus shows genetic signs of fragmentation. In addition, we found that connectivity between full growth native forest remnants is mediated by the presence of dense shrublands and native forest corridors. In this sense, our results are important because they show how native forest patches and associated routes act as source of vector species in silvoagropecuary landscape, increasing the infection risk on human population. This study is the first approach to understand the epidemiological spatial context of silvoagropecuary risk of Hantavirus emergence. Further studies are needed to elucidate the effects of landscape fragmentation in order to generate new predictive models based on vector intrinsic attributes and landscape features.

Genetics, Gene Flow, and Glaciation: The Case of the South American Limpet Nacella mytilina

Glacial episodes of the Quaternary, and particularly the Last Glacial Maximum (LGM) drastically altered the distribution of the Southern-Hemisphere biota, principally at higher latitudes. The irregular coastline of Patagonia expanding for more than 84.000 km constitutes a remarkable area to evaluate the effect of Quaternary landscape and seascape shifts over the demography of near-shore marine benthic organisms. Few studies describing the biogeographic responses of marine species to the LGM have been conducted in Patagonia, but existing data from coastal marine species have demonstrated marked genetic signatures of post-LGM recolonization and expansion. The kelp-dweller limpet Nacella mytilina is broadly distributed along the southern tip of South America and at the Falkland/Malvinas Islands. Considering its distribution, abundance, and narrow bathymetry, N. mytilina represents an appropriate model to infer how historical and contemporary processes affected the distribution of intraspecific genetic diversity and structure along the southern tip of South America. At the same time, it will be possible to determine how life history traits and the ecology of the species are responsible for the current pattern of gene flow and connectivity across the study area. We conducted phylogeographic and demographic inference analyses in N. mytilina from 12 localities along Pacific Patagonia (PP) and one population from the Falkland/Malvinas Islands (FI). Analyses of the mitochondrial gene COI in 300 individuals of N. mytilina revealed low levels of genetic polymorphism and the absence of genetic differentiation along PP. In contrast, FI showed a strong and significant differentiation from Pacific Patagonian populations. Higher levels of genetic diversity were also recorded in the FI population, together with a more expanded genealogy supporting the hypothesis of glacial persistence of the species in these islands. Haplotype genealogy, and mismatch analyses in the FI population recognized an older and more complex demographic history than in PP. Demographic reconstructions along PP suggest a post-LGM expansion process (7.5 ka), also supported by neutrality tests, mismatch distribution and maximum parsimony haplotype genealogies. Migration rate estimations showed evidence of asymmetrical gene flow from PP to FI. The absence of genetic differentiation, the presence of a single dominant haplotype, high estimated migration rates, and marked signal of recent demographic growth, support the hypothesis of rapid post-glacial expansion in N. mytilina along PP. This expansion could have been sustained by larval and rafting-mediated dispersal of adults from northernmost populations following the Cape Horn Current System. Marked genetic differentiation between PP and FI could be explained through differences in their respective glacial histories. During the LGM, Pacific Patagonia (PP) was almost fully covered by the Patagonian Ice Sheet, while sheet coverage in the FI ice was restricted to small cirques and valleys. As previously recorded in the sister-species N. magellanica, the FI rather than represent a classical glacial refugium for N. mytilina, seems to represent a sink area and/or a secondary contact zone. Accordingly, historical and contemporary processes, contrasting glacial histories between the analyzed sectors, as well as life history traits constitute the main factors explaining the current biogeographical patterns of most shallow Patagonian marine benthic organisms.

Phylogeography in Galaxias maculatus (Jenyns, 1848) along Two Biogeographical Provinces in the Chilean Coast

Major geologic and climatic changes during the Quaternary exerted a major role in shaping past and contemporary distribution of genetic diversity and structure of aquatic organisms in southern South America. In fact, the northern glacial limit along the Pacific coast, an area of major environmental changes in terms of topography, currents, and water salinity, represents a major biogeographic transition for marine and freshwater species. We used mitochondrial DNA sequences (D-loop) to investigate the consequences of Quaternary glacial cycles over the pattern of genetic diversity and structure of G. maculatus (Pisces: Galaxiidae) along two biogeographical provinces in the Chilean coast. Extreme levels of genetic diversity and strong phylogeographic structure characterize the species suggesting a low amount of influence of the last glacial cycle over its demography. However, we recognized contrasting patterns of genetic diversity and structure between main biogeographical areas here analyzed. Along the Intermediate Area (38°-41° S) each estuarine population constitutes a different unit. In contrast, Magellanic populations (43°-53° S) exhibited low levels of genetic differentiation. Contrasting patterns of genetic diversity and structure recorded in the species between the analyzed biogeographic areas are consistent with the marked differences in abiotic factors (i.e., different coastal configurations, Quaternary glacial histories, and oceanographic regimes) and to inherent characteristics of the species (i.e., salt-tolerance, physiology, and reproductive behavior).

Comparative phylogeography of co-distributed Phrygilus species (Aves, Thraupidae) from the Central Andes

The Neotropical ecoregion has been an important place of avian diversification where dispersal and allopatric events coupled with periods of active orogeny and climate change (Late Pliocene-Pleistocene) have shaped the biogeography of the region. In the Neotropics, avian population structure has been sculpted not only by geographical barriers, but also by non-allopatric factors such as natural selection and local adaptation. We analyzed the genetic variation of six co-distributed Phrygilus species from the Central Andes, based on mitochondrial and nuclear markers in conjunction with morphological differentiation. We examined if Phrygilus species share patterns of population structure and historical demography, and reviewed the intraspecific taxonomy in part of their geographic range. Our results showed different phylogeographic patterns between species, even among those belonging to the same phylogenetic clade. P. alaudinus, P. atriceps, and P. unicolor showed genetic differentiation mediated by allopatric mechanisms in response to specific geographic barriers; P. gayi showed sympatric lineages in northern Chile, while P. plebejus and P. fruticeti showed a single genetic group. We found no relationship between geographic range size and genetic structure. Additionally, a signature of expansion was found in three species related to the expansion of paleolakes in the Altiplano region and the drying phase of the Atacama Desert. Morphological analysis showed congruence with molecular data and intraspecific taxonomy in most species. While we detected genetic and phenotypic patterns that could be related to natural selection and local adaptation, our results indicate that allopatric events acted as a major factor in the population differentiation of Phrygilus species.

Rapid enzyme-linked immunosorbent assay for the detection of hantavirus-specific antibodies in divergent small mammals

We assessed the utility of an enzyme-linked immunosorbent assay (ELISA) for the detection of hantavirus-specific antibodies from sera of Oligoryzomys longicaudatus, the principal reservoir of Andes virus (ANDV), using an antigen previously developed for detection of antibodies to Sin Nombre virus (SNV) in sera from Peromyscus maniculatus. The assay uses a protein A/G horseradish peroxidase conjugate and can be performed in as little as 1.5 hours. Serum samples from Oligoryzomys longicaudatus collected in central-south Chile were used and the assay identified several that were antibody positive. This assay can be used for the rapid detection of antibodies to divergent hantaviruses from geographically and phylogenetically distant rodent species.

Towards a model of postglacial biogeography in shallow marine species along the Patagonian Province: lessons from the limpet Nacella magellanica (Gmelin, 1791)

BACKGROUND

Patagonia extends for more than 84,000 km of irregular coasts is an area especially apt to evaluate how historic and contemporary processes influence the distribution and connectivity of shallow marine benthic organisms. The true limpet Nacella magellanica has a wide distribution in this province and represents a suitable model to infer the Quaternary glacial legacy on marine benthic organisms. This species inhabits ice-free rocky ecosystems, has a narrow bathymetric range and consequently should have been severely affected by recurrent glacial cycles during the Quaternary. We performed phylogeographic and demographic analyses of N. magellanica from 14 localities along its distribution in Pacific Patagonia, Atlantic Patagonia, and the Falkland/Malvinas Islands.

RESULTS

Mitochondrial (COI) DNA analyses of 357 individuals of N. magellanica revealed an absence of genetic differentiation in the species with a single genetic unit along Pacific Patagonia. However, we detected significant genetic differences among three main groups named Pacific Patagonia, Atlantic Patagonia and Falkland/Malvinas Islands. Migration rate estimations indicated asymmetrical gene flow, primarily from Pacific Patagonia to Atlantic Patagonia (Nem=2.21) and the Falkland/Malvinas Islands (Nem=16.6). Demographic reconstruction in Pacific Patagonia suggests a recent recolonization process (< 10 ka) supported by neutrality tests, mismatch distribution and the median-joining haplotype genealogy.

CONCLUSIONS

Absence of genetic structure, a single dominant haplotype, lack of correlation between geographic and genetic distance, high estimated migration rates and the signal of recent demographic growth represent a large body of evidence supporting the hypothesis of rapid postglacial expansion in this species in Pacific Patagonia. This expansion could have been sustained by larval dispersal following the main current system in this area. Lower levels of genetic diversity in inland sea areas suggest that fjords and channels represent the areas most recently colonized by the species. Hence recolonization seems to follow a west to east direction to areas that were progressively deglaciated. Significant genetic differences among Pacific, Atlantic and Falkland/Malvinas Islands populations may be also explained through disparities in their respective glaciological and geological histories. The Falkland/Malvinas Islands, more than representing a glacial refugium for the species, seems to constitute a sink area considering the strong asymmetric gene flow detected from Pacific to Atlantic sectors. These results suggest that historical and contemporary processes represent the main factors shaping the modern biogeography of most shallow marine benthic invertebrates inhabiting the Patagonian Province.