Species delimitation and integrative taxonomy of the Reithrodontomys mexicanus (Rodentia: Cricetidae) cryptic complex

Species boundaries are difficult to establish in groups with very similar morphology. As an alternative, it has been suggested to integrate multiple sources of data to clarify taxonomic problems in taxa where cryptic speciation processes have been reported. This is the case of the harvest mouse Reithrodontomys mexicanus, which has a problematic taxonomy history as it is considered a complex species. Here, we evaluate the cryptic diversity of R. mexicanus using an integrative taxonomy approach in order to detect candidate lineages at the species level. The molecular analysis used one mitochondrial (cytb) and two nuclear (Fgb-I7 and IRBP) genes. Species hypotheses were suggested based on three molecular delimitation methods (mPTP, bGMYC, and STACEY) and cytb genetic distance values. Skull and environmental space differences between the delimited species were also tested to complement the discrimination of candidate species. Based on the consensus across the delimitation methods and genetic distance values, four species were proposed, which were mostly supported by morphometric and ecological data: R. mexicanus clade I, R. mexicanus clade IIA, R. mexicanus clade IIIA, and R. mexicanus clade IIIB. In addition, the evolutionary relationships between the species that comprise the R. mexicanus group were discussed from a phylogenetic approach. Our findings present important taxonomic implications for Reithrodontomys, as the number of known species for this genus increases. Furthermore, we highlight the importance of the use of multiple sources of data in systematic studies to establish robust delimitations between species considered taxonomically complex.

Pliocene sigmodontine rodents (Mammalia: Cricetidae) in northernmost South America: test of biogeographic hypotheses and revised evolutionary scenarios

We document the first occurrence of Sigmodontinae (Mammalia, Rodentia, Cricetidae) from the Pliocene of northern South America, from the San Gregorio Formation of northwestern Venezuela. The recovered isolated molars are identified as Oligoryzomys sp. and Zygodontomys sp., two currently widespread sigmodontines in South America. These records constitute the oldest representatives of these genera, potentially new species, and the first Pliocene occurrence for Oryzomyini and the whole subfamily outside Argentina. Hypotheses on the historical biogeography of sigmodontines have been constructed almost exclusively using genetic data and the fossils we report provide a new kind of evidence. The occurrence of Oligoryzomys sp. and Zygodontomys sp. in Venezuela provides novel information for the diversification models suggested for Oligoryzomys, by supporting a potential eastern corridor of open environments from northern to southern South America. The presence of sigmodontines from the locality home of the new reports, Norte Casa Chiguaje, is consistent with the palaeoenvironmental conditions originally proposed for it based on mammals and botanical records, being characterized as mixed open grassland/forest areas surrounding permanent freshwater systems. The new sigmodontine evidence is used to discuss the putative scenarios of the ancient evolution of the subfamily in South America, favouring a model in which open areas (savannahs) to the east of the Andes played crucial role aiding or obstructing Late Miocene-Pliocene sigmodontine dispersion southwards.

Modeling potential risk areas of Orthohantavirus transmission in Northwestern Argentina using an ecological niche approach

BACKGROUND

Hantavirus Pulmonary Syndrome (HPS) is a rodent-borne zoonosis in the Americas, with up to 50% mortality rates. In Argentina, the Northwestern endemic area presents half of the annually notified HPS cases in the country, transmitted by at least three rodent species recognized as reservoirs of Orthohantavirus. The potential distribution of reservoir species based on ecological niche models (ENM) can be a useful tool to establish risk areas for zoonotic diseases. Our main aim was to generate an Orthohantavirus risk transmission map based on ENM of the reservoir species in northwest Argentina (NWA), to compare this map with the distribution of HPS cases; and to explore the possible effect of climatic and environmental variables on the spatial variation of the infection risk.

METHODS

Using the reservoir geographic occurrence data, climatic/environmental variables, and the maximum entropy method, we created models of potential geographic distribution for each reservoir in NWA. We explored the overlap of the HPS cases with the reservoir-based risk map and a deforestation map. Then, we calculated the human population at risk using a census radius layer and a comparison of the environmental variables' latitudinal variation with the distribution of HPS risk.

RESULTS

We obtained a single best model for each reservoir. The temperature, rainfall, and vegetation cover contributed the most to the models. In total, 945 HPS cases were recorded, of which 97,85% were in the highest risk areas. We estimated that 18% of the NWA population was at risk and 78% of the cases occurred less than 10 km from deforestation. The highest niche overlap was between Calomys fecundus and Oligoryzomys chacoensis.

CONCLUSIONS

This study identifies potential risk areas for HPS transmission based on climatic and environmental factors that determine the distribution of the reservoirs and Orthohantavirus transmission in NWA. This can be used by public health authorities as a tool to generate preventive and control measures for HPS in NWA.

Geometric morphometrics and ecological niche modelling for delimitation of Triatoma pallidipennis (Hemiptera: Reduviidae: Triatominae) haplogroups

A recent phylogenetic analysis of Triatoma pallidipennis, an important Chagas disease vector in Mexico, based on molecular markers, revealed five monophyletic haplogroups with validity as cryptic species. Here, we compare T. pallidipennis haplogroups using head and pronotum features, environmental characteristics of their habitats, and ecological niche modeling. To analyze variation in shape, images of the head and pronotum of the specimens were obtained and analyzed using methods based on landmarks and semi-landmarks. Ecological niche models were obtained from occurrence data, as well as a set of bioclimatic variables that characterized the environmental niche of each analyzed haplogroup. Deformation grids for head showed a slight displacement towards posterior region of pre-ocular landmarks. Greatest change in head shape was observed with strong displacement towards anterior region of antenniferous tubercle. Procrustes ANOVA and pairwise comparisons showed differences in mean head shape in almost all haplogroups. However, pairwise comparisons of mean pronotum shape only showed differences among three haplogroups. Correct classification of all haplogroups was not possible using discriminant analysis. Important differences were found among the environmental niches of the analyzed haplogroups. Ecological niche models of each haplogroup did not predict the climatic suitability areas of the other haplogroups, revealing differences in environmental conditions. Significant differences were found between at least two haplogroups, demonstrating distinct environmental preferences among them. Our results show how the analysis of morphometric variation and the characterization of the environmental conditions that define the climatic niche can be used to improve the delimitation of T. pallidipennis haplogroups that constitute cryptic species.

A comparison of gene organisations and phylogenetic relationships of all 22 squamate species listed in South Korea using complete mitochondrial DNA

Studies using complete mitochondrial genome data have the potential to increase our understanding on gene organisations and evolutionary species relationships. In this study, we compared complete mitochondrial genomes between all 22 squamate species listed in South Korea. In addition, we constructed Maximum Parsimony (MP), Maximum Likelihood (ML) and Bayesian Inference (BI) phylogenetic trees using 13 mitochondrial protein-coding genes. The mitochondrial genes for all six species in the suborder Sauria followed the same organisation as the sequenced Testudines (turtle) outgroup. In contrast, 16 snake species in the suborder Serpentes contained some gene organisational variations. For example, all snake species contained a second control region (CR2), while three species in the family Viperidae had a translocated tRNA-Pro gene region. In addition, the snake species, Elapheschrenckii, carried a tRNA-Pro pseudogene. We were also able to identify a translocation of a tRNA-Asn gene within the five tRNA (WANCY gene region) gene clusters for two true sea snake species in the subfamily Hydrophiinae. Our BI phylogenetic tree was also well fitted against currently known Korean squamate phylogenetic trees, where each family and genus unit forms monophyletic clades and the suborder Sauria is paraphyletic to the suborder Serpentes. Our results may form the basis for future northeast Asian squamate phylogenetic studies.

Phylogeography and evolutionary lineage diversity in the small-eared greater galago, Otolemur garnettii (Primates: Galagidae)

Assessing the true lineage diversity in elusive nocturnal organisms is particularly challenging due to their subtle phenotypic variation in diagnostic traits. The cryptic small-eared greater galago (Otolemur garnettii) offers a great opportunity to test if currently recognized subspecies, suggested by discontinuities in coat colour pattern and geographic barriers, represent distinct evolutionary lineages. To answer this question, we conducted the first population-level phylogeographic study of the species, sampling wild specimens from across almost its entire latitudinal range, including the Zanzibar Archipelago. We applied five species-delimitation algorithms to investigate the genetic diversity and distribution pattern of mitochondrial DNA across the geographic range of three out of four subspecies. Our results suggest that far-northern populations of O. g. lasiotis potentially represent an independently evolving lineage, but populations assigned to O. g. garnettii from Zanzibar Island and of O. g panganiensis from mainland Tanzania do not constitute two independent lineages. A dated phylogeny suggests that this northern clade diverged from all remaining samples approximately 4 Mya. Such old divergence age is in line with the split between many galagid species. This northern lineage could potentially represent an incipient species; however, there is not yet enough evidence to support a new taxonomic status for this unique mitochondrial group.

Molecular data confirm Triatoma pallidipennis Stål, 1872 (Hemiptera: Reduviidae: Triatominae) as a novel cryptic species complex

Triatoma pallidipennis constitues one of the most important Chagas disease vector in Mexico. Previous studies based on molecular data suggest T. pallidipennis as a complex of cryptic species. For that reason, we analyzed the phylogenetic relationships of T. pallidipennis using DNA sequences from the mitochondrial ND4 gene and the ITS-2 gene. In addition, the divergence times were estimated, and possible new taxa were delimited with three species delimitation methods. Finally, genetic distances and possible connectivity routes based on shared haplotypes were obtained among the T. pallidipennis populations. Five haplogroups (possible cryptic species) were found, based on delimitation methods and genetic distances. Haplogroup divergence began about 3 Ma, in the Pleistocene. Moreover, none of the haplogroups showed potential connectivity routes between them, evidencing lack of gene flow. Our results suggest the existence of a new cryptic species complex within what is currently recognized as a T. pallidipennis.

Molecular systematics of the Reithrodontomys tenuirostris group (Rodentia: Cricetidae) highlighting the Reithrodontomys microdon species complex

The Reithrodontomys tenuirostris species group is considered "the most specialized" within the genus Reithrodontomys from morphological and ecological perspectives. Previous studies based on molecular data recommended changes in the taxonomy of the group. In particular, R. microdon has been the most taxonomically questioned, with the suggestion that it constitutes a complex of cryptic species. We analyzed the phylogenetic relationships of the R. tenuirostris species group using DNA sequences from the mitochondrial Cytochrome b gene and Intron 7 of the nuclear beta fibrinogen gene. In addition, divergence times were estimated, and possible new taxa delimited with three widely used species delimitation methods. Finally, possible connectivity routes based on shared haplotypes were tested among the R. microdon populations. All species were recovered as monophyletic with the exception of R. microdon, whose individuals were grouped into four different haplogroups, one of which included specimens of R. bakeri. Diversification within the R. tenuirostris species group began about 3 Ma, in the Pleistocene. The bGMYC and STACEY delimitation methods were congruent with each other, delimiting at the species-level each haplogroup within R. microdon, while the mPTP suggested a greater number of species. Moreover, none of the haplogroups showed potential connectivity routes between them, evidencing lack of gene flow. Our results suggest the existence of a higher number of species in the R. tenuirostris group, because we show that there are four species within what is currently recognized as R. microdon.

Expanded distribution of orthohantavirus Pergamino genotype: First record in Entre Rios province, Argentina

Hantavirus Pulmonary Syndrome (HPS) is an emerging infectious disease of the Americas. Eight native rodent species have been identified as HPS virus reservoirs in Argentina. The aim of this work was to detect the orthohantavirus genotypes present in a rodent community that inhabits a zone where a fatal HPS case occurred within an endemic locality of Central Argentina. We captured 27 rodents with a trapping effort of 723 trap nights. We detected 14.3% of infected Akodon azarae with the Pergamino genotype. This result expands the known distribution of this orthohantavirus. Although the Pergamino genotype has not been associated with human cases, the information about its distribution is relevant for risk assessment against potential changes in the virus infectivity.

Identifying Chagas disease vectors using elliptic Fourier descriptors of body contour: a case for the cryptic dimidiata complex

BACKGROUND

Triatoma dimidiata (Reduviidae: Triatominae) is an important vector of Chagas disease in various countries in the Americas. Phylogenetic studies have defined three lineages in Mexico and part of Central America. While there is a marked genetic differentiation, methods for identifying them using morphometric analyses with landmarks have not yet been fully resolutive. Elliptical Fourier descriptors (EFDs), which mathematically describe the shape of any closed two-dimensional contours, could be a potentially useful alternative method. Our objective was to validate the use of EFDs for the identification of three lineages of this species complex.

METHOD

A total of 84 dorsal view images of individuals of the three lineages were used. Body contours were described with EFDs using between 5 and 30 harmonics. The number of obtained coefficients was reduced by a principal components analysis and the first axis scores were used as shape variables. A linear discriminant function analysis and an ordination plot of the discriminant analysis were performed using the shape variables. A confusion matrix of the ordination plot of the discriminant analysis was obtained to estimate the classification errors, the first five PC scores were statistically compared, and a neural network were then performed using the shape variables.

RESULTS

The first principal component explained 50% of the variability, regardless the number of harmonics used. The results of discriminant analysis get improved by increasing the number of harmonics and components considered. With 25 harmonics and 30 components, the identification of haplogroups was achieved with an overall efficiency greater than 97%. The ordering diagram showed the correct discrimination of haplogroups, with only one error of discrimination corroborated by the confusion matrix. When comparing the first five PC scores, significant differences were found among at least two haplogroups. The 30 multilayer perceptron neural networks were also efficient in identification, reaching 91% efficiency with the validation data.

CONCLUSIONS

The use of EFD is a simple and useful method for the identification of the main lineages of Triatoma dimidiata, with high values of correct identification.

New data about home range and movements of Oligoryzomys flavescens (Rodentia: Cricetidae) help to understand the spread and transmission of Andes virus that causes Hantavirus Pulmonary Syndrome

Hantavirus pulmonary syndrome is an emerging infectious disease caused by viruses of the genus Orthohantavirus. The rodent Oligoryzomys flavescens is distributed along four countries of South America. In Argentina, O. flavescens acts as a reservoir of three genotypes of ANDV orthohantavirus. The aims of this work were to estimate home range size and movements-with spool-and-line and radiotelemetry-of infected and non-infected O. flavescens in order to understand the spread and transmission of the virus. O. flavescens use a wide area to satisfice its requirements, reaching a home range of 1.82 ha during spring. Orthohantavirus infection did not change the behaviour of individuals. We observed a great overlapping in the home range of infected and non-infected individuals resulting in a high probability of virus dispersion on rodent population. These results show that human health risks could be high on island environments and knowledge about the movement ecology of O. flavescens provides useful information on prevention.

DNA barcoding of the rodent genus Oligoryzomys (Cricetidae: Sigmodontinae): mitogenomic-anchored database and identification of nuclear mitochondrial translocations (Numts)

DNA barcoding has become a standard method for species identification in taxonomically complex groups. An important step of the barcoding process is the construction of a library of voucher-based material that was properly identified by independent methods, free of inaccurate identification, and paralogs. We provide here a cytochrome oxidase I (mt-Co1) DNA barcode database for species of the genus Oligoryzomys, based on type material and karyotyped specimens, and anchored on the mitochondrial genome of one species of Oligoryzomys, O. stramineus. To evaluate the taxonomic determination of new COI sequences, we assessed species intra/interspecific genetic distances (barcode gap), performed the General Mixed Yule Coalescent method (GMYC) for lineages' delimitation, and identified diagnostic nucleotides for each species of Oligoryzomys. Phylogenetic analyses of Oligoryzomys were performed on 2 datasets including 14 of the 23 recognized species of this genus: a mt-Co1 only matrix, and a concatenated matrix including mt-Co1, cytochrome b (mt-Cytb), and intron 7 of the nuclear fibrinogen beta chain gene (i7Fgb). We recovered nuclear-mitochondrial translocated (Numts) pseudogenes on our samples and identified several published sequences that are cases of Numts. We analyzed the rate of non-synonymous and synonymous substitution, which were higher in Numts in comparison to mtDNA sequences. GMYC delimitations and DNA barcode gap results highlight the need for further work that integrate molecular, karyotypic, and morphological analyses, as well as additional sampling, to tackle persistent problems in the taxonomy of Oligoryzomys.

Orthohantavirus genotype Lechiguanas in Oligoryzomys nigripes (Rodentia: Cricetidae): New evidence of host-switching

To identify and predict situations of increased risk of orthohantavirus infection in humans, it is necessary to study the relationships between the virus and its rodent hosts. The present study investigated orthohantavirus infection in an assemblage of wild Sigmodontinae rodents of the Paraná Delta, Argentina, and providing new evidence of host-switching events. Rodents belonging to the species Oxymycterus rufus (n = 187), Akodon azarae (n = 82), Oligoryzomys flavescens (n = 80), Oligoryzomys nigripes (n = 47), Scapteromys aquaticus (n = 38), Deltamys kempi (n = 7) and Holochilus brasiliensis (n = 2) were captured at 4 sampling sites during 20 trapping sessions. Blood samples were analyzed by IgG ELISA and livers by a nested reverse transcription PCR for the diagnosis of orthohantavirus infection. The amplified products of the S and M orthohantavirus genomes were sequenced and analyzed to determine similarities with species of the Orthohantavirus genus. The species of the Oligoryzomys positive to the virus were confirmed by amplifying and sequencing the complete cyt b gene. Of the 443 serum samples analyzed by IgG ELISA, A. azarae presented the highest host-specific prevalence value (10/82, 12.2%) followed by Ol. nigripes (4/47, 8.5%) and Ox. rufus (1/187, 0.5%). All the sero-positive Ol. nigripes (n = 4) were positive to the amplification of the S and M segments of the Lechiguanas genotype (98% nucleotide identity for both segments). This is surprising given that Ol. nigripes has been previously associated with Juquitiba genotype, not Lechiguanas. The latter is generally associated with Ol. flavescens, which in our study were all sero-negative. In addition, the association Ox. rufus - Pergamino genotype found here is, to our knowledge, novel and another potential evidence of host-switching considering that Pergamino has been originally associated with A. azarae. These findings contribute to the building evidence that contradicts the one-genotype-one-reservoir species premise in the association between rodent reservoirs and orthohantaviruses, and supports the hypothesis that the community structure of sympatric host species may contribute to orthohantavirus dynamics.