Population and species boundaries in the South American subterranean rodent Ctenomys in a dynamic environment

Underground speciation: Unraveling the systematics and evolution of the highly diverse tuco-tucos (genus Ctenomys) with genomic data

Subterranean rodents of the genus Ctenomys (tuco-tucos) are endemic to South America and have experienced relatively recent radiation. There are about 67 recognized species that originated in approximately 1-2 MY. They stand out for their species richness, extraordinary chromosomal diversity, and wide range of habitat they occupy in the continent. Phylogenetic relationships among species of tuco-tucos have been challenging to resolve. Groups of closely-related species have been suggested, but their relationships must be resolved. This study estimates the phylogeny of the genus using massive sequencing, generating thousands of independent molecular markers obtained by RADseq, with a taxonomic sampling that includes 66% of the recognized species. The sequences obtained were mapped against the C. sociabilis genome, recovering up to 1,215 widely shared RAD loci with more than 19,000 polymorphic sites. Our new phylogenetic hypothesis corroborated the species groups previously proposed with cytochrome b gene sequences and provided a much greater resolution of the relationships among species groups. The frater group is sister to all other tuco-tucos, whereas some of the earlierliest proposals placed the sociabilis group as sister to all other tuco-tucos. Ctenomys leucodon, previously proposed as an independent lineage, is associated with the frater group with moderate statistical support. The magellanicus and mendocinus are sister groups in a major clade formed by the boliviensis, talarum, tucumanus, torquatus, and opimus groups. Ctenomys viperinus, included in the phylogeny for the first time, belongs to the tucumanus group. This multi-locus phylogenetic hypothesis provides insights into the historical biogeography of understanding this highly diverse genus.

Topographic barriers drive the pronounced genetic subdivision of a range-limited fossorial rodent

Due to their limited dispersal ability, fossorial species with predominantly belowground activity usually show increased levels of population subdivision across relatively small spatial scales. This may be exacerbated in harsh mountain ecosystems, where landscape geomorphology limits species' dispersal ability and leads to small effective population sizes, making species relatively vulnerable to environmental change. To better understand the environmental drivers of species' population subdivision in remote mountain ecosystems, particularly in understudied high-elevation systems in Africa, we studied the giant root-rat (Tachyoryctes macrocephalus), a fossorial rodent confined to the afro-alpine ecosystem of the Bale Mountains in Ethiopia. Using mitochondrial and low-coverage nuclear genomes, we investigated 77 giant root-rat individuals sampled from nine localities across its entire ~1000 km2 range. Our data revealed a distinct division into a northern and southern group, with no signs of gene flow, and higher nuclear genetic diversity in the south. Landscape genetic analyses of the mitochondrial and nuclear genomes indicated that population subdivision was driven by slope and elevation differences of up to 500 m across escarpments separating the north and south, potentially reinforced by glaciation of the south during the Late Pleistocene (~42,000-16,000 years ago). Despite this landscape-scale subdivision between the north and south, weak geographic structuring of sampling localities within regions indicated gene flow across distances of at least 16 km at the local scale, suggesting high, aboveground mobility for relatively long distances. Our study highlights that despite the potential for local-scale gene flow in fossorial species, topographic barriers can result in pronounced genetic subdivision. These factors can reduce genetic variability, which should be considered when developing conservation strategies.

The contribution of incomplete lineage sorting and introgression to the evolutionary history of the fast-evolving genus Ctenomys (Rodentia, Ctenomyidae)

Incomplete lineage sorting and introgression have been increasingly recognized as important processes involved in biological differentiation. Both incomplete lineage sorting and introgression result in incongruences between gene trees and species trees, consequently causing difficulties in phylogenetic reconstruction. This is particularly the case for rapid radiations, as short internodal distances and incomplete reproductive isolation increase the likelihood of both ILS and introgression. Estimation of the relative frequency of these processes requires assessments across many genomic regions. We use transcriptomics to test for introgression and estimate the frequency of incomplete lineage sorting in a set of three closely related and geographically adjacent South American tuco-tucos species (Ctenomys), a genus comprising 64 species resulting from recent, rapid radiation. After cleaning and filtering, 5764 orthologous genes strongly support paraphyly of C. pearsoni relative to C. brasiliensis (putatively represented by the population of Villa Serrana). In line with earlier phylogenetic work, the C. pearsoni - C. brasiliensis pair is closely related to C. torquatus, whereas C. rionegrensis is more distantly related to these three nominal species. Classical Patterson's D-statistic shows significant signals of introgression from C. torquatus into C. brasiliensis. However, a 5-taxon test shows no significant results. Incomplete lineage sorting was estimated to have involved about 9% of the loci, suggesting it represents an important process in the incipient diversification of tuco-tucos.

Redefining the Distributional Boundaries and Phylogenetic Relationships for Ctenomids From Central Argentina

With about 68 recognized living species, subterranean rodents of the genus Ctenomys are found in a multiplicity of habitats, from the dunes of the Atlantic coast to the Andes Mountains, including environments ranging from humid steppes of Pampas to the dry deserts of Chaco region. However, this genus needs an exhaustive reevaluation of its systematic and phylogenetic relationships regarding the different groups that compose it. This knowledge is essential to propose biodiversity conservation strategies both at species level and at higher hierarchical levels. In order to clarify the taxonomy and the recent evolutionary history from populations of Ctenomys in the Pampas region, Argentina, phylogenetic relationships among them were evaluated using mitochondrial DNA sequences: gene encoding cytochrome b protein (1,140 bp) and the non-coding D-loop region (434 bp). To infer the divergence times inside the Ctenomys clade, a Bayesian calibrate tree using fossil remains data from different families within Caviomorpha was performed at first. Secondly, that calibration data was used as priors in a new Bayesian phylogenetic inference within the genus Ctenomys. This phylogenetic tree emphasized on species currently distributed on the Pampas region, more precisely considering both the talarum and mendocinus groups. Bayesian inferences (BI) were integrated with the results of a Maximum Likelihood approach (ML). Based on these results, the distributional limits of the mendocinus and talarum groups appear to be related to the physiognomy of the Pampas region soils. On the other hand, the validity of C. pundti complex as a differentiated species of C. talarum is debated. According to previous evidence from morphological and chromosomal studies, these results show a very low divergence between those species that originally were classified within the talarum group. Mitochondrial DNA sequences from populations associated with these putative species have not recovered as reciprocal monophyletic groups in the phylogenetic analyses. In conclusion, C. talarum and C. pundti complex might be considered as the same biological species, or lineages going through a recent or incipient differentiation process. The results obtained in this study have important implications for conservation policies and practices, since both species are currently categorized as Vulnerable and Endangered, respectively.

Genetic diversity and conservation of the endemic tuco-tuco Ctenomys ibicuiensis (Rodentia: Ctenomyidae)

Endemic, small range species are susceptible to environmental changes and landscape modification. Understanding genetic diversity and distributional patterns is important for implementation of effective conservation measures. In this context, genetic diversity was evaluated to update the conservation status of an endemic tuco-tuco, Ctenomys ibicuiensis. Phylogeographic and population genetic analyses of mitochondrial DNA and microsatellite loci were carried out using 46 individuals sampled across the species’ distribution. Ctenomys ibicuiensis presented moderate to high genetic diversity and highly structured populations with low levels of gene flow and isolation by distance. Anthropogenic landscape changes threaten this restricted-range tuco-tuco. Considering its limited geographic distribution and highly structured populations with low gene flow, we consider C. ibicuiensis to be at significant risk of extinction.

Evolutionary and ecological patterns within the South African Bathyergidae: Implications for taxonomy

The family Bathyergidae (comprising six genera) is a group of subterranean rodents endemic to sub-Saharan Africa. Our understanding of the evolution and species richness of the South African bathyergid genera Georychus, Bathyergus and Cryptomys is limited, with the majority of species listed as Least Concern by the IUCN Red List of Threatened Species. Genetic data suggest that several cryptic species may be present in these genera. To explore genetic and ecological distinctiveness, and evaluate taxonomic richness across the ranges of Georychus, Bathyergus and to a lesser degree, Cryptomys, as well as evaluate possible scenarios which have historically influenced evolutionary patterns, we employed four protein coding markers (one mitochondrial and three nuclear) along with distribution wide sampling schemes and large sample sizes. In addition, possible ecological differences among the different intra-generic clades were explored. Genera appear to have originated in the north-eastern interior of South Africa, following novel habitats created through the Post-African I erosion cycle and dramatic changes in climate and phytogeography. In each genus, multiple geographically discrete genetic lineages (clades) are supported by both the mitochondrial and nuclear data. These lineages bear signature of the fragmentation of wider historical distributions through major environmental changes since the middle Miocene (major uplift events, Post-African II erosion cycle, drainage evolution of major river systems, sea-level fluctuations as well as climatic changes and vegetation shifts), thereby leading to long-term isolation. Along with protracted periods of separation, it appears that ecological differences further delimit the lineages in relation to geology, phytogeographic preference, elevation, rainfall and temperature. As such, two lineages in Georychus (Clades 1 and 2) and one lineages in Cryptomys (Clade I) occur at higher elevations above the Great Escarpment (in older deposits harbouring grassland vegetation, with higher rainfall and lower daily temperatures), with the remaining lineages within these genera (Clades 3, 4 and 5 in Georychus and Clades III and IV in Cryptomys) occupying a low-land distribution with contrasting climatic and geological characteristics. Although significant differences in ecological variables were also observed between Bathyergus clades, these were not consistent, given their largely low-land distributions. Our results corroborate and expand previous suggestions that several cryptic species are present within the South African Bathyergidae.

Integrative analysis of chromosome banding, telomere localization and molecular genetics in the highly variable Ctenomys of the Corrientes group (Rodentia; Ctenomyidae)

The genus Ctenomys comprises about 70 species with great chromosome diversity. The Corrientes group is one of the most chromosomally variable lineages in the genus, where the diploid number (2n) varies from 41 to 70. In this group, three nominal species and numerous polymorphic and polytypic populations have been described. In order to get insight into the chromosomal evolution of this species complex, we applied different banding and molecular cytogenetic techniques. The results were interpreted in an evolutionary context, based on mitochondrial cytochrome b analyses. Studied samples are representative of the broad chromosomal variability in the group, including specimens with 2n = 42 to 2n = 70. Heterochromatin was scarce but concentrated in a few chromosomes. Centromeric DAPI-negative heterochromatin was observed in some autosomal pairs, which differed among populations. Location and amount of DAPI-neutral heterochromatin within the Y chromosome varied among populations. The variable distribution of heterochromatin indicates its dynamic behavior. NORs were detected in one pair of autosomes, which also differed among some populations. Telomeric FISH signals were observed in all complements only at the chromosome ends. The Corrientes group belongs to a clade that also includes C. pearsoni, C. lami, C. minutus, C. ibicuiensis and C. torquatus. Almost all of these species are variable at the chromosomal level, suggesting that this is the ancestral condition of the clade. Within the Corrientes group, the observed low genetic divergence, in contrast with its high chromosomal variability, is indicative of decoupling between the rates of chromosomal and mitochondrial evolution.

Spatial genetic diversity in the Cape mole-rat, Georychus capensis: Extreme isolation of populations in a subterranean environment

The subterranean niche harbours animals with extreme adaptations. These adaptations decrease the vagility of taxa and, along with other behavioural adaptations, often result in isolated populations characterized by small effective population sizes, high inbreeding, population bottlenecks, genetic drift and consequently, high spatial genetic structure. Although information is available for some species, estimates of genetic diversity and whether this variation is spatially structured, is lacking for the Cape mole-rat (Georychus capensis). By adopting a range-wide sampling regime and employing two variable mitochondrial markers (cytochrome b and control region), we report on the effects that life-history, population demography and geographic barriers had in shaping genetic variation and population genetic patterns in G. capensis. We also compare our results to information available for the sister taxon of the study species, Bathyergus suillus. Our results show that Georychus capensis exhibits low genetic diversity relative to the concomitantly distributed B. suillus, most likely due to differences in habitat specificity, habitat fragmentation and historical population declines. In addition, the isolated nature of G. capensis populations and low levels of population connectivity has led to small effective population sizes and genetic differentiation, possibly aided by genetic drift. Not surprisingly therefore, G. capensis exhibits pronounced spatial structure across its range in South Africa. Along with geographic distance and demography, other factors shaping the genetic structure of G. capensis include the historical and contemporary impacts of mountains, rivers, sea-level fluctuations and elevation. Given the isolation and differentiation among G. capensis populations, the monotypic genus Georychus may represent a species complex.

Landscape genetics in the subterranean rodent Ctenomys "chasiquensis" associated with highly disturbed habitats from the southeastern Pampas region, Argentina

Studies of genetic differentiation in fragmented environments help us to identify those landscape features that most affect gene flow and dispersal patterns. Particularly, the assessment of the relative significance of intrinsic biological and environmental factors affecting the genetic structure of populations becomes crucial. In this work, we assess the current dispersal patterns and population structure of Ctenomys "chasiquensis", a vulnerable and endemic subterranean rodent distributed on a small area in Central Argentina, using 9 polymorphic microsatellite loci. We use landscape genetics approaches to assess the relationship between genetic connectivity among populations and environmental attributes. Our analyses show that populations of C. "chasiquensis" are moderately to highly structured at a regional level. This pattern is most likely the outcome of substantial gene flow on the more homogeneous sand dune habitat of the Northwest of its distributional range, in conjunction with an important degree of isolation of eastern and southwestern populations, where the optimal habitat is surrounded by a highly fragmented landscape. Landscape genetics analysis suggests that habitat quality and longitude were the environmental factors most strongly associated with genetic differentiation/uniqueness of populations. In conclusion, our results indicate an important genetic structure in this species, even at a small spatial scale, suggesting that contemporary habitat fragmentation increases population differentiation.

Integrative lineage delimitation in rodents of the Ctenomys Corrientes group

The tuco-tucos rodents (genus Ctenomys) of the Corrientes group comprise several populations that inhabit the vast area under the influence of the Iberá wetland. Lineage delimitation within the recently diverged Corrientes group is a challenging task as morphological differentiation is not conspicuous between populations. However, delimitation is crucial for evolutionary studies and conservation issues. In this study, we performed a phylogenetic analysis including cytochrome b (cyt-b) sequences from taxa that had never been studied in a comprehensive context. We integrated previously published chromosomal studies, mitochondrial phylogenies and simple sequence repeat (SSR) variability analyses, and applied a delimitation criterion over the basis of chromosomal incompatibilities and genetic exclusivity. Under this integrative approach seven independently evolving lineages were delimited in the Corrientes group: Ctenomys roigi, which conserves its former definition, Ctenomys dorbignyi and Ctenomys perrensi complex which were redefined, Sarandicito which includes the population of Paraje Sarandicito and probably a group of nearby poorly studied populations, and Iberá i, ii and iii distributed at both sides of the Iberá wetland. We discuss future perspectives to evaluate the proposed lineages and conservation issues concerning these tuco-tucos.

Wet soils affect habitat selection of a solitary subterranean rodent (Ctenomys minutus) in a Neotropical region

Subterranean rodents are characterized by limited individual mobility and patchy distribution of local populations. Habitat patches where the species is absent may either be unoccupied, but suitable patches, or in fact unsuitable as habitat due to specific habitat features. In the coastal plain of southern Brazil, Ctenomys minutus (Ctenomyidae) inhabits sandy grasslands and dunes. The aim of our study was to analyze the relationship between the distribution of this subterranean rodent and its environment. We considered vegetation and soil features. Our results showed that habitat occupancy patterns of C. minutus are determined by soil moisture and vegetation (food supply), but there might be other habitat features that regulate the habitat choice of this mammal. Habitat discontinuities in the coastal plain of southern Brazil, with numerous barriers to dispersal, are responsible for restricting individuals to their natal areas. This might result in increased intraspecific competition within each population and, consequently, in inbreeding.

Os roedores subterrâneos são caracterizados pela baixa mobilidade dos indivíduos e por apresentarem populações distribuídas em manchas de habitat. As manchas em que os indivíduos estão ausentes podem ser locais desocupados, porém que apresentam condições adequadas para a ocupação, ou de fato serem locais inadequados devido a algumas características ambientais específicas. Na planície costeira do sul do Brasil, a espécie Ctenomys minutus (Ctenomyidae) habita regiões de campos arenosos e dunas costeiras. O objetivo do presente estudo foi analisar a relação entre a distribuição deste roedor subterrâneo e o ambiente em que ele habita. Para tal, consideramos características da vegetação e do solo. Nossos resultados demonstraram que o padrão de ocupação de C. minutus é determinado pela umidade do solo e vegetação (recurso alimentar), mas podem existir outras características que regulam a escolha do habitat deste mamífero. As descontinuidades de habitat na planície costeira do sul do Brasil, com inúmeras barreiras para a dispersão da espécie, são responsáveis por restringir os indivíduos as suas áreas natais. Isso pode resultar no aumento da competição intra-específica dentro das populações e, consequentemente, na endogamia.

Local and regional scale genetic variation in the Cape dune mole-rat, Bathyergus suillus

The distribution of genetic variation is determined through the interaction of life history, morphology and habitat specificity of a species in conjunction with landscape structure. While numerous studies have investigated this interplay of factors in species inhabiting aquatic, riverine, terrestrial, arboreal and saxicolous systems, the fossorial system has remained largely unexplored. In this study we attempt to elucidate the impacts of a subterranean lifestyle coupled with a heterogeneous landscape on genetic partitioning by using a subterranean mammal species, the Cape dune mole-rat (Bathyergus suillus), as our model. Bathyergus suillus is one of a few mammal species endemic to the Cape Floristic Region (CFR) of the Western Cape of South Africa. Its distribution is fragmented by rivers and mountains; both geographic phenomena that may act as geographical barriers to gene-flow. Using two mitochondrial fragments (cytochrome b and control region) as well as nine microsatellite loci, we determined the phylogeographic structure and gene-flow patterns at two different spatial scales (local and regional). Furthermore, we investigated genetic differentiation between populations and applied Bayesian clustering and assignment approaches to our data. Nearly every population formed a genetically unique entity with significant genetic structure evident across geographic barriers such as rivers (Berg, Verlorenvlei, Breede and Gourits Rivers), mountains (Piketberg and Hottentots Holland Mountains) and with geographic distance at both spatial scales. Surprisingly, B. suillus was found to be paraphyletic with respect to its sister species, B. janetta-a result largely overlooked by previous studies on these taxa. A systematic revision of the genus Bathyergus is therefore necessary. This study provides a valuable insight into how the biology, life-history and habitat specificity of animals inhabiting a fossorial system may act in concert with the structure of the surrounding landscape to influence genetic distinctiveness and ultimately speciation.

Bioclimatic, ecological, and phenotypic intermediacy and high genetic admixture in a natural hybrid of octoploid strawberries

PREMISE OF THE STUDY

Hybrid zones provide "natural laboratories" for understanding the processes of selection, reinforcement, and speciation. We sought to gain insight into the degree of introgression and the extent of ecological-phenotypic intermediacy in the natural hybrid strawberry, Fragaria × ananassa subsp. cuneifolia. •

METHODS

We used whole-plastome sequencing to identify parental species-specific (Fragaria chiloensis and F. virginiana) chloroplast single-nucleotide polymorphisms and combined the use of these with nuclear microsatellite markers to genetically characterize the hybrid zone. We assessed the potential role of selection in the observed geographic patterns by bioclimatically characterizing the niche of the hybrid populations and phenotypically characterizing hybrid individuals of known genomic constitution. •

KEY RESULTS

Significant admixture and little overall maternal bias in chloroplast or nuclear genomes suggest a high degree of interfertility among the parental and hybrid species and point to a long history of backcrossing and genetic mixing in the hybrid zone. Even though hybrids were phenotypically intermediate to the parental species, there was a discernible fingerprint of the parental genotype within hybrid individuals. Thus, although the pattern of introgression observed suggests geographic limitations to gene flow, it may be reinforced by selection for specific parental traits in the bioclimatically intermediate habitat occupied by the hybrid. •

CONCLUSIONS

This work uncovered the genetic complexity underlying the hybrid zone of the wild relatives of the cultivated strawberry. It lays the foundation for experimental dissection of the causes of genomic introgression and nuclear-cytoplasmic disassociation, and for understanding other parts of Fragaria evolutionary history.

Human impact in naturally patched small populations: genetic structure and conservation of the burrowing rodent, tuco-tuco (Ctenomys lami)

Isolated or semi-isolated small populations are commonly found among species, due to a naturally patchy occupancy of suitable habitats or also as a result of habitat alterations. These populations are subject to an increased risk of local extinction because they are more vulnerable to demographic, genetic, and environmental stochasticity. Considering that natural areas have been becoming progressively more fragmented and smaller, understanding the genetic structure and evolutionary dynamics of small populations is critical. Ctenomys lami has 26 karyotypes distributed in a small area (936 km(2)) continually modified by human actions. We assessed the genetic geographical structure of this species, examining 178 specimens sampled on a fine scale, using information from chromosomal variability, mitochondrial DNA control region and cytochrome c oxidase subunit I sequences, and 14 microsatellite loci. The observed isolation-by-distance pattern and a clinal genetic variation suggest a stepping-stone population model. The results did not indicate genetic structuring associated with distinct karyotypes. However, mitochondrial and nuclear molecular markers demonstrated the existence of 2 demes, which are not completely isolated but are probably reinforced by a geographical barrier. The vulnerability of C. lami is greater than previously supposed, and our data support the designation of one Evolutionary Significant Unit and one Management Unit, and also the inclusion of this species' conservation status as vulnerable.

Evidence of environmental niche differentiation in the striped mouse (Rhabdomys sp.): inference from its current distribution in southern Africa

The aim of this study was to characterize environmental differentiation of lineages within Rhabdomys and provide hypotheses regarding potential areas of contact between them in the Southern African subregion, including the Republic of South Africa, Lesotho, and Namibia. Records of Rhabdomys taxa across the study region were compiled and georeferenced from the literature, museum records, and field expeditions. Presence records were summarized within a 10 × 10 km grid covering the study area. Environmental information regarding climate, topography, land use, and vegetation productivity was gathered at the same resolution. Multivariate statistics were used to characterize the current environmental niche and distribution of the whole genus as well as of three mitochondrial lineages known to occur in southern Africa. Distribution modeling was carried out using MAXENT in order to generate hypotheses regarding current distribution of each taxa and their potential contact zones. Results indicate that the two species within Rhabdomys appear to have differentiated across the precipitation/temperature gradient present in the region from east to west. R. dilectus occupies the wettest areas in eastern southern Africa, while R. pumilio occupies the warmer and drier regions in the west, but also penetrates in the more mesic central part of the region. We provide further evidence of environmental differentiation within two lineages of R. dilectus. Contact zones between lineages appear to occur in areas of strong environmental gradients and topographic complexity, such as the transition zones between major biomes and the escarpment area where a sharp altitudinal gradient separates coastal and plateau areas, but also within more homogeneous areas such as within grassland and savannah biomes. Our results indicate that Rhabdomys may be more specialized than previously thought when considering current knowledge regarding mitochondrial lineages. The genus appears to have differentiated along two major environmental axes in the study region, but results also suggest dispersal limitations and biological interactions having a role in limiting current distribution boundaries. Furthermore, the projection of the potential geographic distribution of the different lineages suggests several contact zones that may be interesting study fields for understanding the interplay between ecological and evolutionary processes during speciation.