Next Article in Journal
Inducible Expression of Several Drosophila melanogaster Genes Encoding Juvenile Hormone Binding Proteins by a Plant Diterpene Secondary Metabolite, Methyl Lucidone
Next Article in Special Issue
New Species, New Records, and Updated Key to the Paravelia (Hemiptera: Heteroptera: Veliidae) from Brazil
Previous Article in Journal
Effects of Trichoderma harzianum Strain T22 on the Arthropod Community Associated with Tomato Plants and on the Crop Performance in an Experimental Field
Previous Article in Special Issue
A New Species in Notobitiella Hsiao in China Confirmed by Morphological and Molecular Data (Hemiptera: Heteroptera: Coreidae)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Testing Phylogeographic Hypotheses in Mepraia (Hemiptera: Reduviidae) Suggests a Complex Spatio-Temporal Colonization in the Coastal Atacama Desert

by
Ricardo Campos-Soto
1,*,
Evelyn Rodríguez-Valenzuela
2,
Gabriel Díaz-Campusano
2,
Dusan Boric-Bargetto
2,
Álvaro Zúñiga-Reinoso
3,
Franco Cianferoni
2 and
Fernando Torres-Pérez
2
1
Escuela de Ciencias Agrícolas y Veterinarias, Universidad Viña del Mar, Viña del Mar 2572007, Chile
2
Instituto de Biología, Facultad de Ciencias, Pontificia Universidad Católica de Valparaíso, Valparaíso 2373223, Chile
3
Institute for Zoology, University of Cologne, 50674 Cologne, Germany
*
Author to whom correspondence should be addressed.
Insects 2022, 13(5), 419; https://doi.org/10.3390/insects13050419
Submission received: 28 February 2022 / Revised: 21 April 2022 / Accepted: 25 April 2022 / Published: 29 April 2022

Abstract

:

Simple Summary

Mepraia is a blood-sucking bug endemic to Chile and a vector of the parasite that causes Chagas disease. Different colonization routes have been suggested for this bug; therefore, we tested different colonization routes using DNA sequences and bioinformatics approaches to select the most probable route. Our results suggest that, after the split of Triatoma, Mepraia divided into two main groups ~2.1 Mya. The northern group would have speciated between 1.7–1.4 Mya, giving rise to M. parapatrica, M. gajardoi and to a new, still undescribed lineage (Mepraia sp.). The southern group formed M. spinolai ~1.68 Mya. We suggest that Mepraia originated from the north-central Andes due to the last Andes uplift and hyperaridity. The hyperarid cycle would have separated the southern and northern groups. Then, within the northern group, colonization would have occurred from the centre to the north and south through corridors influenced by Pleistocene climatic changes. The habitat colonized by the southern clade led to only one species (M. spinolai). Fluctuations in climatic changes probably influenced speciation strongly in this kissing bug in the Atacama Desert.

Abstract

Mepraia is a genus (Triatominae) endemic to Chile and a vector of Trypanosoma cruzi. Alternative phylogeographic hypotheses have been suggested for Mepraia. We tested different colonization routes hypothesized using mitochondrial sequences and phylogeographic approaches to select the best-supported hypothesis. Our results suggest that, after the split from the sister genus Triatoma at ~4.3 Mya, Mepraia formed two main clades at ~2.1 Mya. The northern clade diverged from Mepraia sp. ~1.7 Mya, giving rise to M. parapatrica and M. gajardoi about ~1.4 Mya. The southern clade originated M. spinolai ~1.68 Mya. We suggest that Mepraia had an origin in the north-central Andes along with orogenic processes, reinforced by hyperaridity during the Pliocene. The hyperarid cycle would have separated the southern and northern clades. Then, in the northern clade, dispersal occurred north and south from the centre through corridors during the Pleistocene Climatic Oscillations. Climate changes may have induced a major speciation process in the Atacama Desert, while the more homogeneous habitat colonized by the southern clade led to only one, but structured, species.

1. Introduction

One of the most important historical events that shaped the current biodiversity in South America is, without a doubt, the Andes uplift, because this huge mountain range is a clear east–west barrier and because both climate and landscape changed as a consequence of the uplift, producing high diversification in this continent [1]. These landscape changes allowed some taxa to disperse to other places because new niches were available [2], whereas other taxa were enclosed in small, suitable patches [3], possibly triggering species pump events. The origin and diversification of the tribe Triatomini seem to be linked with the events explained above, which allowed the colonization of southern South America from a Mesoamerican origin [4]. The genus Mepraia is a Triatomini with a sub-Andean distribution and could be one of the latest groups to diversify. This taxon is endemic to the arid and semiarid regions of Chile; it is known widely because Mepraia species are blood-sucking bugs responsible for the transmission of Trypanosoma cruzi in the wild cycle, causing Chagas disease [5,6]. Mepraia is composed of three species: M. gajardoi, M. parapatrica and M. spinolai. M. gajardoi and M. parapatrica have been found exclusively in the coastal Atacama Desert between 18°30′ S and 26°51′ S [7,8], whereas M. spinolai inhabits from 26° to 34° S, mainly in the interior of valleys and, to a lesser extent, in coastal areas [7,8] (Figure 1). In addition to these species, a new coastal lineage has been reported in Chile (M. sp. from now on), which inhabits between the distributions of M. parapatrica and M. gajardoi (23°25′ to 23°28′ S) [9]. The coastal species inhabit landscapes with low vegetation under stones and in rock crevices, associated with marine bird nests, rodent burrows and lizards [9,10]. Although each of the Mepraia species forms well-supported clades, the phylogenetic relationships between these lineages are still unresolved, leading to controversy regarding the colonization routes of Mepraia species in Chile [7,9]. One of the colonization hypotheses of Mepraia states that, due to the divergence of an ancestral population by the uplift of the Andes Range, Mepraia would have separated from the eastern slope species Triatoma eratyrusiformis and T. breyeri. Later, the process of speciation in Mepraia would have followed latitudinal colonization from north to south [8]. According to this hypothesis, the most northerly species M. gajardoi would be sister to the remaining species (M. parapatrica and M. spinolai) and descended from the most ancient ancestor. However, this hypothesis is not consistent with the time-calibrated phylogeny, which suggests that the most southerly species M. spinolai is sister to the remaining species (M. gajardoi, M. parapatrica and M. sp.) [7,9]. Therefore, new alternative colonization scenarios need to be evaluated. In this study, we used sequences of three mitochondrial genes to assess spatio-temporal scenarios that could explain the historical colonization routes of Mepraia.

2. Materials and Methods

2.1. Samples, Localities and Obtaining Mitochondrial Gene Sequences

We used a set of DNA samples extracted from previous studies [7,9], which were stored at −80 °C. These samples were from localities that covered the distribution ranges of all Mepraia species, and a T. eratyrusiformis sample used as outgroup (Figure 1, Table S1). Segments of 631 bp for cytochrome oxidase subunit I (COI, called COIA) and 563 bp for NADH dehydrogenase subunit 4 (ND4) of each selected sample were amplified by PCR using the polymerase SapphireAmp® fast PCR Master Mix (Takara, Kusatsu shi, Japan). We used the primers COIF (5′-CCTGCAGGAGGAGGAGAYCC) and CO10R (5′-TAAGCGTCTGGGTAGTCTGARTAKCG) used in [11] for the COIA segment, and the primers ND4-F (5′-TCAACATGAGCCCTTGGAAG) and ND4-R (5′-TAATTCGTTGTCATGGTAATG) described in [12] for the ND4 segment. The same cycling conditions were used to amplify the COIA and ND4 segments using a Bioer model TC-96/G/H(b)C LifeEco® thermocycler (Hangzhou, China). The cycling conditions were 30 s at 94 °C, followed by 40 cycles of 94 °C for 30 s, 55 °C for 30 s and 72 °C for 1 min, according to the manufacturer’s instructions. Amplification was confirmed by electrophoresis in 1% agarose gels. The sequencing of both strands was performed by Macrogen, Inc. (Seoul, Korea) using the same PCR primers. Sequences were edited to obtain consensus sequences using BioEdit 7.0.4.1 [13] and aligned using Clustal W [14] incorporated in BioEdit. After alignment, sites that showed nucleotide substitutions were re-examined by the visual inspection of each individual’s raw chromatogram. The COI and ND4 gene sequences were deposited in GenBank with the accession numbers MN117859-MN117888. The samples sequenced for each locality and GenBank accession numbers are described in Table S1. The COIA and ND4 sequences of each sample were manually concatenated with a previous matrix [9] formed by another COI segment designated as COIB and a cyt b segment, obtaining a final matrix of 2237 pb (Table S1) composed of four segments from three genes. New COIA and ND4 and previous sequences [7,9] were obtained from the same kissing bugs for each locality. The nucleotide substitution model was selected with the Bayesian Information Criterion (BIC) using Smart Model Selection [15] in the ATGC bioinformatics platform [16]. The best model was HKY + G for both COI segments and TN93 + G for the cyt b and ND4 segments.

2.2. Phylogenetic Analyses and Estimates of Time of the Most Recent Common Ancestor (TMRCA)

Phylogenetic reconstruction was inferred with the final concatenated matrix (2237 bp), which was partitioned by gene and codon position; trees were linked using BEAUti, implemented in BEAST v.2.6.4 [18]. Phylogenetic relationships and TMRCA were estimated with the Bayesian Markov Chain Monte Carlo (BMCMC) method available in BEAST v.2.6.4 [18]. A Yule model was used as prior, run first using the uncorrelated lognormal relaxed molecular clock, which was calibrated using the mutation rates previously estimated for each gene. A mutation rate of 0.016 substitutions/site/million years (Myr) was used for the COI gene segment (COIA + COIB) [19], 0.023 for cyt b [20] and 0.11 for ND4 [21]. These mutation rates were used as means with a normal distribution in the priors for the molecular clock to estimate the best mutation rates for each segment. The estimated mutation rates (substitutions/site/Myr) were 0.016 for COI, 0.023 for cyt b and 0.047 for ND4. Using these rates and the Yule model as prior, uncorrelated lognormal relaxed and strict molecular clocks were run. The best molecular clock was selected using the Nested Sampling package [22] available in BEAST v.2.6.4. using 100 particles in each model run. This analysis estimated the log Marginal Likelihood (log(ML)) for each model. The log Bayes Factors (log (BF)) were calculated using the log(ML) of each model in the equation: log(BF) = log(ML1) − log(ML2). If log(BF) is greater than 1, the first model is favoured [23]. The best resulting model estimated was the strict molecular clock. With the Yule model and selected molecular clock, two separate MCMC runs of 5 × 107 generations were run using the BEAST program, sampling every 5000 generations. Mixing and convergence were verified based on the effective sample size (ESS) values of each run in Tracer v.1.7 [24]; only ESSs > 500 were accepted in each run. In both runs, the initial 10% of trees were discarded and the rest were combined using LogCombiner. The statistical uncertainty was depicted in values of the 95% Highest Probability Density (HPD) and node support was evaluated with posterior probabilities (pp); only pp > 0.95 were accepted as supported. The consensus tree was calculated with the maximum clade credibility and mean node heights option using TreeAnnotator. The consensus tree was visualized and edited in FigTree v1.4.4. [25].

2.3. Test for Phylogeographic Scenarios Using Coalescent Simulations

In order to better understand the diversification process of Mepraia in Chile, we used the coalescent-based Approximate Bayesian Computation (ABC) method implemented in DIYABC v2.2 [26]. In this approach, different models (scenarios) are compared to judge how well each fits the observed data by generating simulated data using historical and mutational parameters. ABC analyses were conducted on the sequence data of the COIB and cyt b segments of all individuals (Table S1). Based on the phylogenetic results, we considered four groups (M. gajardoi, M. sp., M. parapatrica and M. spinolai) and constructed five phylogeographic scenarios. Scenario 1 represents the hypothesis generated by the BMCMC method [9] that suggests colonization from the Coquimban Biogeographical Province (centre-south) [17] to the north (Figure 1a). In this scenario, the most southerly species M. spinolai is sister to the northern clade formed by M. gajardoi, M. parapatrica and M. sp., and, within this northern clade, M. sp. is sister to the remaining species M. gajardoi and M. parapatrica (Figure 2). Scenario 2 represents the hypothesis of colonization from north to south proposed by Frías-Lasserre (2010) [8] (Figure 2). In scenario 3 M. spinolai is sister to the northern clade; however, within this northern clade, M. gajardoi is sister to the remaining species M. sp. and M. parapatrica (Figure 2). Scenario 4 suggests colonization from south to north following a parapatric speciation model, where M. spinolai is sister to the northern clade, which includes M. parapatrica as sister to the remaining species M. gajardoi and M. sp. (Figure 2). Scenario 5 proposes a southern clade formed by M. spinolai and M. sp. and a northern clade formed by M. gajardoi and M. parapatrica (Figure 2). We used the estimations calculated in this study by the BMCMC method to set the divergence times for the scenarios (Table S2), and the population sizes were set as constant, as suggested by [7]. We used HKY as the mutational model and a mutation rate of 1.00 × 10−8–5.00 × 10−8 per site per generation calculated according to the mutation rate of mitochondrial markers (10−2 subs/site/my) and a generation time of one year for Mepraia species (Table S2). For each scenario, we simulated 1,000,000 datasets and performed pre-evaluations to confirm that our models were suitable for comparison against the observed data. The best scenario was chosen by comparing the posterior probabilities (pp), which were calculated using a logistic regression approach on 1% of the simulated datasets closest to the observed dataset.

2.4. Test of Phylogenetic Hypotheses Using Bayes Factor

The phylogeographic scenarios were also compared using BEAST v.2.6.4 to provide more evidence on the diversification process of Mepraia. With the final concatenated matrix (COI, cyt b and ND4), five files (xml) were formed with different phylogenetic relationships congruent with scenarios 1 to 5 (Figure 2) using the option “Add Prior” in BEAUti and keeping all other parameters unchanged. Each scenario was run separately with 100 particles and log(ML) was calculated using the Nested Sampling package [22] available in BEAST v.2.6.5. Paired comparisons between scenarios were performed with log(BF), which was calculated with log(ML); if log(BF) is larger than 1, the first scenario is favoured [23]. We assumed the topology generated by the BMCMC method as the first scenario, which was compared with the other scenarios.

3. Results

3.1. Phylogenetic Relationships and Divergences Time into Mepraia Species

The phylogenetic BMCMC analyses showed that Mepraia split from the Argentinian taxon about 4.3 Mya (HPD: 3.88–4.76). Mepraia species were clustered in two supported main clades north and south (Figure 3a). The northern clade was divided into three supported subclades that included M. gajardoi, M. parapatrica and M. sp.; all occurred along the long coastal strip of northern Chile. The southern clade contained only M. spinolai (Figure 3a). The estimated TMRCA between the two main clades was about 2 Mya (HPD: 1.85–2.25). Within the northern clade, the estimated TMRCA for Mepraia species was 0.62 Mya (HPD: 0.48–0.78) for M. gajardoi, 0.45 Mya (HPD: 0.36–0.55) for M. parapatrica and 0.16 Mya (HPD: 0.08–0.23) for M. sp. M. gajardoi was recovered as a sister taxon of M. parapatrica and the divergence time of these sister groups was 1.4 Mya (HPD: 1.22–1.59). The TMRCA for the southern M. spinolai clade was 1.68 Mya (HPD: 1.5–1.9).

3.2. Test of Phylogeographic Hypotheses with Coalescent Simulations

Approximate Bayesian Computation analyses showed that the scenario with the highest posterior probability was scenario 1, with pp = 0.8 (0.7891–0.8490), followed by scenario 4, pp = 0.112 (0.0868–0.1374), scenario, 5 pp = 0.033 (0.0194–0.0455), scenario 3, pp = 0.0297 (0.0217–0.0378), and finally scenario 2, pp = 0.0066 (0.0017–0.0116) (Figure 4).

3.3. Test of Phylogenetic Hypotheses

Paired comparisons between the phylogenetic scenarios showed scenario 1 as the best phylogenetic model. The calculated log(BF) for scenario 1 had positive values above 8.8 in all comparisons, which can be interpreted as overwhelming support in favour of scenario 1. The estimated log(ML) and paired comparisons between scenarios using log(BF) are given in Table 1. The calculated standard deviations were 1.2 for all estimated log(ML).

4. Discussion

The effects of the climate changes in South America after the Andes uplift and, in particular, the Plio-Pleistocene climate oscillations (PPCO) were notorious; they caused a strong impact on the landscape and biota in different ecoregions of the continent [28,29]. Although the effects on the biota in the desert areas of Chile have been less studied, the PPCO is considered as one of the drivers of the diversification process, generating high endemism of lineages [3,29,30,31,32]. The hypothesis that was most strongly supported by our results shows that the split between Mepraia and Argentinean Triatoma species (~4.3 Mya, HPD: 3.88–4.76) likely occurred after the last Andes uplift (Figure 3b), as suggested by Frias-Lasserre (2010) [8]; our results partially support this idea. The last uplift of the Andes was estimated at ~5–3 Mya during the late Pliocene [33]; however, during the Miocene, the Andes already had an altitude similar to today [34]. Therefore, this divergence could have been reinforced by the Pliocene climate changes, which made the weather colder and drier [35], even going through a hyperarid period in the Atacama Desert [27] (Figure 3b). These hyperarid or arid periods were also suggested as the cause of divergence for clades of other insects, reinforcing the role of tectonic and climate changes in the separation of main groups [3,30]. This divergent event would have occurred in the north-central Andes between ~26°–30° S, as in Psectrascelis darkling beetles [3] and plant species of the genus Argylia [36].
After this basal split, the history of Mepraia (~2.1 Mya) began during the Plio-Pleistocene (Figure 3d,e). The climate changes were particularly harsh in the Atacama Desert, switching from arid to hyperarid in the last 25 Mya [27]. Today, the Atacama Desert is hyperarid and most insects avoid the core of the Atacama, occupying other environments, such as the Altiplano [3] and/or, like Mepraia, the narrow coastal band. From 28° S, the climate becomes more “friendly”, and roughly from here to the south, M. spinolai can occupy inland areas. Today, M. spinolai has inhabited the most humid area since one of the most important hyperarid periods at the end of the Pliocene. This could have pushed this species to survive in small suitable habitats, definitely separating M. spinolai from the northern clade ancestor. Therefore, our results strongly support the colonization of Mepraia from the central highland area to the north (Figure 3d), as was expected from the results of several earlier studies [7,9]. This colonization pathway also was suggested for other insects, where the ancestor could have occupied the Coquimban Biogeographical Province (where M. spinolai is distributed) and then colonized the northern Atacama during the Pliocene [3,32]. In fact, the Coquimban Biogeographical Province was demonstrated as the central distribution area for M. spinolai, with more diversity than the peripheral populations [37]. The distribution pattern of Mepraia today nicely represents the effect of hyperaridity on the distribution of these triatomines; the climate cycling during the Plio-Pleistocene surely affected the diversification of Mepraia. However, Mepraia are fully dependent on other animals to complete their biological cycle. Mepraia have opportunistic feeding behaviour, i.e., they mainly feed on the most abundant or available hosts, including mammals, lizards and birds [38,39]. Therefore, the phylogeographic patterns of Mepraia are not only affected by abiotic forces [37], but also depend on the status and distribution of the animal populations on which they feed. In other words, the phylogeography of Mepraia is also dependent on the availability and distribution of hosts. High genetic structure is expected under these conditions, and populations of Mepraia are highly structured [7]. This pattern is congruent with the other blood-sucking bug (Triatoma infestans) [40], which is the only species of Triatoma in Chile and is sympatric with M. spinolai [41]. The Pleistocene Climatic Oscillation (PCO) generated the extension (in humid periods) and contraction (in dry periods) of semiarid vegetation [42,43]. It is well known that the semiarid vegetation harboured vertebrate fauna during the humid periods of the PCO [44,45] that could be a source of food for triatomines that followed the structuring process. Climatic oscillations in the distribution range of the southern clade had a minor impact, creating a more homogeneous habitat [32] that produced fewer splits or speciation processes. This habitat could have favoured the connection of populations, but, even under those favourable conditions, a deep genetic structure was found within M. spinolai [7,37].
The PCO would have had a greater impact on the northern clade than on the southern clade. This greater impact was expressed as a major diversification between 1.8 and 1.4 Mya, which drove the origin of three species (Figure 3f), each one also genetically quite structured. All intraspecific structuring processes would have occurred in the last 0.9 Mya, which is congruent with the other colonization process in surrounding areas of the Atacama Desert [3,30].
The dry period of climatic oscillations in the area of the northern clade would have caused major population structuring, separating populations by hyperarid barriers, the contraction of semiarid vegetation, and connection by corridors [3,30,31,36,45,46]. Several vegetated corridors, such as the “Lomas” on the coast of the Atacama Desert, have been suggested [36,47]. These vegetated corridors were established and maintained by the humidity of the coastal fog known as “camanchaca”, and allowed the colonization and development of flora and fauna [36,46,47,48]. These corridors were already established at the time of the divergence of the northern and southern clades. Mepraia of the northern clade may have occupied these corridors to colonize north and south from a central origin (M. sp. Figure 3d–g), feeding on the fauna that occupied these foggy, isolated habitats. This kind of habitat has been present in coastal Atacama from the Middle Neogene [28,49,50]. Today, those three species mainly inhabit these fog oases, such as Paposo-Pan de Azúcar (M. parapatrica), Morro Moreno-Isla Santa María (M. sp.) and Alto Patache-Quebrada Vitor (M. gajardoi) (Figure 1a). This distribution pattern would reinforce the inter–intraspecific genetic structuring by hyperarid barriers even more.
PCO also affected the marine terraces due to rising sea levels, which were frequent in the late Pleistocene [51,52,53]. These changes in sea level and marine terraces had an impact on the structure of these habitats, probably separating populations on coastal islands [9] and causing the extinction of some populations that inhabited the littoral plains during high-sea-level periods. The extinctions caused by these changes left unoccupied habitats that could have been recolonized by ancient populations that were left in refuges. In particular, M. sp. could have originated from an ancestral refuge, but the diversification was very recent (~0.16 Mya), which is congruent with the estimated dates of changes that impacted the coast in this zone [52,54]. The extinction dynamics by sea-level changes left a very recent fingerprint on the genetic diversity of this probable cryptic lineage (M. sp.), even though it is the oldest lineage of the northern clade, and can also explain its small distribution range enclosed between Santa María Island and Morro Moreno, where only females have been collected (Figure 1a).
We acknowledge the limitations of our results that only reflect the maternal history and that it is necessary to integrate nuclear genes to produce a more robust evolutionary history. However, these results provide valuable first spatio-temporal colonization evidence that can be tested with additional nuclear markers in the near future.

5. Conclusions

This study presents the first spatio-temporal analysis to show the colonization history of a blood-sucking insect in the Atacama Desert. After the last Andes uplift, Mepraia was definitively separated from its Argentinean sister group, which was reinforced by colder and more arid periods during the Pliocene. Subsequently, our results support the colonization of Mepraia from the Coquimban Biogeographical Province to the north. The current coastal distribution range of Mepraia nicely illustrates the border of the hyperarid core of the Atacama, which is a real-time picture of the effects of hyperaridity on the distribution of those insects. The colonization of those hyperarid environments should have occurred during the late Pliocene; later, by allopatry, a hyperarid cycle would have separated definitively the southern and northern clades (Figure 3). Today, the southern clade inhabits a more suitable environment, allowing it to occupy inland habitats, while the northern clade is pressured by the hyperarid core of the Atacama to inhabit a narrow coastal strip. In this coastal strip, there are a few fog oases that have a suitable habitat for surviving in this harsh environment. The isolation between these oases triggered population and species differentiation. The genetic structure is also likely affected by biotic factors, because Mepraia is a blood-sucking insect, but their effects would be less than those of the abiotic factors. Hyperarid barriers repeated several times during the PPCO, shaping the current diversity and genetic structure of Mepraia. The effect of these oscillations was greater in the Atacama Desert, where the diversity is higher, suggesting that it caused a major speciation process. This is because colonization should have occurred from the centre to the north and south following vegetated corridors during the humid periods of the Pleistocene, which were interrupted by long hyperarid periods. Finally, these results could provide relevant knowledge to understand the infection history of the protozoan Trypanosoma cruzi under a probable biogeographic scenario and help to understand the eco-epidemiology of Chagas disease in Chile.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/insects13050419/s1, Table S1: Sampling localities and accession numbers. Table S2: Historical and mutational parameters used in the coalescent simulations analysis.

Author Contributions

Conceptualization, investigation, methodology, formal analysis, supervision, writing—original draft, writing—review & editing, funding acquisition, project administration, R.C.-S.; Investigation, methodology, formal analysis, writing—review & editing, E.R.-V.; Investigation, methodology, writing—review & editing, project administration, G.D.-C.; Methodology, writing—review & editing, investigation, D.B.-B., Á.Z.-R. and F.C.; Conceptualization, resources, writing—review & editing, F.T.-P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Chilean National Agency for Research and Development (ANID) program National Fund for Scientific and Technological Development (FONDECYT) grant number 11170643 and Universidad Viña del Mar Project FIIUVM-CTC-2211, both to R. Campos-Soto.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the manuscript and in its Supplementary Material (Tables S1 and S2).

Acknowledgments

We would like to thank Carezza Botto-Mahan for sharing photographs of Mepraia specimens and Oscar Inostroza-Michel for helping with the analyses.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Graham, A. The Andes: A geological overview from a biological perspective. Ann. Mo. Bot. Gard. 2009, 96, 371–385. [Google Scholar] [CrossRef]
  2. Hughes, C.; Eastwood, R. Island radiation on a continental scale: Exceptional rates of plant diversification after uplift of the Andes. Proc. Natl. Acad. Sci. USA 2006, 103, 10334–10339. [Google Scholar] [CrossRef] [Green Version]
  3. Zúñiga-Reinoso, A.; Ritter, B.; Predel, R. The colonization of the Puna and Atacama Biogeographic Province by sister clades of Psectrascelis (Coleoptera: Tenebrionidae): Synchronous expansion without spatial overlap. J. Biogeogr. 2021, 48, 1930–1940. [Google Scholar] [CrossRef]
  4. Monteiro, F.; Weirauch, C.; Felix, M.; Lazoski, C.; Abad-Franch, F. Evolution, Systematics, and Biogeography of the Triatominae, Vectors of Chagas Disease. In Advances in Parasitology; Elsevier: Amsterdam, The Netherlands, 2018; pp. 265–344. [Google Scholar]
  5. Campos-Soto, R.; Ortiz, S.; Ivan, C.; Bruneau, N.; Botto-Mahan, C.; Solari, A. Interactions between Trypanosoma cruzi the Chagas disease parasite and naturally infected wild Mepraia vectors of Chile. Vector-Borne Zoonotic Dis. 2016, 16, 165–171. [Google Scholar] [CrossRef] [Green Version]
  6. Botto-Mahan, C.; Bacigalupo, A.; Correa, J.P.; Fontúrbel, F.E.; Cattan, P.E.; Solari, A. Prevalence, infected density or individual probability of infection? Assessing vector infection risk in the wild transmission of Chagas disease. Proc. R. Soc. B Biol. Sci. 2020, 287, 20193018. [Google Scholar] [CrossRef] [Green Version]
  7. Campos, R.; Torres-Pérez, F.; Botto-Mahan, C.; Coronado, X.; Solari, A. High phylogeographic structure in sylvatic vectors of Chagas disease of the genus Mepraia (Hemiptera: Reduviidae). Infect. Genet. Evol. 2013, 19, 280–286. [Google Scholar] [CrossRef] [PubMed]
  8. Frías-Lasserre, D. A new species and karyotype variation in the bordering distribution of Mepraia spinolai (Porter) and Mepraia gajardoi Frías et al. (Hemiptera: Reduviidae: Triatominae) in Chile and its parapatric model of speciation. Neotrop. Entomol. 2010, 39, 572–583. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Campos-Soto, R.; Díaz-Campusano, G.; Cianferoni, F.; Torres-Pérez, F. Biogeographic origin and phylogenetic relationships of Mepraia (Hemiptera, Reduviidae) on islands of northern Chile. PLoS ONE 2020, 15, e0234056. [Google Scholar] [CrossRef] [PubMed]
  10. Campos-Soto, R.; Díaz-Campusano, G.; Quiroga, N.; Muñoz-San Martín, C.; Rives-Blanchard, N.; Torres-Pérez, F. Trypanosoma cruzi-infected triatomines and rodents co-occur in a coastal island of northern Chile. PeerJ 2020, 8, e9967. [Google Scholar] [CrossRef] [PubMed]
  11. Calleros, L.; Panzera, F.; Bargues, M.D.; Monteiro, F.A.; Klisiowicz, D.R.; Zuriaga, M.A.; Mas-Coma, S.; Pérez, R. Systematics of Mepraia (Hemiptera-Reduviidae): Cytogenetic and molecular variation. Infect. Genet. Evol. 2010, 10, 221–228. [Google Scholar] [CrossRef]
  12. Grisales, N.; Triana, O.; Angulo, V.; Jaramillo, N.; Parra-Henao, G.; Panzera, F.; Gómez-Palacio, A. Diferenciación genética de tres poblaciones colombianas de Triatoma dimidiata (Latreille, 1811) mediante análisis molecular del gen mitocondrial ND4. Biomédica 2010, 30, 207–214. [Google Scholar] [CrossRef] [Green Version]
  13. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. 1999, 41, 95–98. [Google Scholar] [CrossRef] [Green Version]
  14. Thompson, J.D.; Higgins, D.G.; Gibson, T.J. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22, 4673–4680. [Google Scholar] [CrossRef] [Green Version]
  15. Lefort, V.; Longueville, J.E.; Gascuel, O. SMS: Smart Model Selection in PhyML. Mol. Biol. Evol. 2017, 34, 2422–2424. [Google Scholar] [CrossRef] [Green Version]
  16. Guindon, S.; Gascuel, O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol. 2003, 52, 696–704. [Google Scholar] [CrossRef] [Green Version]
  17. Morrone, J.J. Biogeographical regionalisation of the Andean region. Zootaxa 2015, 3936, 207–236. [Google Scholar] [CrossRef] [Green Version]
  18. Bouckaert, R.; Vaughan, T.G.; Barido-Sottani, J.; Duchêne, S.; Fourment, M.; Gavryushkina, A.; Heled, J.; Jones, G.; Kühnert, D.; De Maio, N.; et al. BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis. PLoS Comput. Biol. 2019, 15, e1006650. [Google Scholar] [CrossRef] [Green Version]
  19. Arensburger, P.; Buckley, T. Biogeography and phylogeny of the New Zealand cicada genera (Hemiptera: Cicadidae) based on nuclear and mitochondrial DNA data. J. Biogeogr. 2004, 31, 557–569. [Google Scholar] [CrossRef] [Green Version]
  20. Monteiro, F.; Barrett, T.V.; Fitzpatrick, S.; Cordon-Rosales, C.; Feliciangeli, D.; Beard, C.B. Molecular phylogeography of the Amazonian Chagas disease vectors Rhodnius prolixus and R. robustus. Mol. Ecol. 2003, 12, 997–1006. [Google Scholar] [CrossRef]
  21. Pech-May, A.; Mazariegos-Hidalgo, C.J.; Izeta-Alberdi, A.; López-Cancino, S.A.; Tun-Ku, E.; de la Cruz-Félix, K.; Ibarra-Cerdeña, C.N.; González Ittig, R.E.; Ramsey, J.M. Genetic variation and phylogeography of the Triatoma dimidiata complex evidence a potential center of origin and recent divergence of haplogroups having differential Trypanosoma cruzi and DTU infections. PLoS Negl. Trop. Dis. 2019, 13, e0007044. [Google Scholar] [CrossRef] [Green Version]
  22. Maturana, P.; Brewer, B.J.; Klaere, S.; Bouckaert, R. Model selection and parameter inference in phylogenetics using Nested Sampling. Syst. Biol. 2019, 68, 219–233. [Google Scholar] [CrossRef]
  23. Kass, R.E.; Raftery, A.E. Bayes Factors. J. Am. Stat. Assoc. 1995, 90, 773–795. [Google Scholar] [CrossRef]
  24. Rambaut, A.; Drummond, A.J.; Xie, D.; Baele, G.; Suchard, M.A. Posterior summarization in Bayesian phylogenetics using Tracer 1.7. Syst. Biol. 2018, 67, 901–904. [Google Scholar] [CrossRef] [Green Version]
  25. Rambaut, A. Molecular Evolution, Phylogenetics and Epidemiology. FigTree. Tree Figure Drawing Tool Software v1.4.4. 2018. Available online: http://tree.bio.ed.ac.uk/software/figtree/ (accessed on 18 February 2022).
  26. Cornuet, J.-M.; Pudlo, P.; Veyssier, J.; Dehne-Garcia, A.; Gautier, M.; Leblois, R.; Marin, J.-M.; Estoup, A. DIYABC v2.0: A software to make approximate Bayesian computation inferences about population history using single nucleotide polymorphism, DNA sequence and microsatellite data. Bioinformatics 2014, 30, 1187–1189. [Google Scholar] [CrossRef] [Green Version]
  27. Jordan, T.E.; Kirk-Lawlor, N.E.; Blanco, N.P.; Rech, J.A.; Cosentino, N.J. Landscape modification in response to repeated onset of hyperarid paleoclimate states since 14 Ma, Atacama Desert, Chile. GSA Bull. 2014, 126, 1016–1046. [Google Scholar] [CrossRef]
  28. Hinojosa, L.F.; Villagrán, C. Historia de los bosques del sur de Sudamérica, I: Antecedentes paleobotánicos, geológicos y climáticos del Terciario del cono sur de América. Rev. Chil. Hist. Nat. 1997, 70, 225–239. [Google Scholar]
  29. Latorre, C.; Moreno, P.I.; Vargas, G.; Maldonado, A.; Villa-Martínez, R.; Armesto, J.J.; Villagrán, C.; Pino, M.; Núñez, L.; Grosjean, M. Late Quaternary environments and palaeoclimate. In The Geology of Chile; Moreno, T., Gibbons, W., Eds.; Geological Society of London: London, UK, 2007; pp. 309–328. [Google Scholar]
  30. Zúñiga-Reinoso, Á.; Predel, R. Past climatic changes and their effects on the phylogenetic pattern of the Gondwanan relict Maindronia (Insecta: Zygentoma) in the Chilean Atacama Desert. Glob. Planet. Chang. 2019, 182, 103007. [Google Scholar] [CrossRef]
  31. Boric-Bargetto, D.; Zúñiga-Reinoso, Á.; Inostroza-Michel, O.; Rodríguez-Serrano, E.; González-Acuña, D.; Palma, R.E.; Hernández, C.E. A comprehensive overview of the genetic diversity in Thylamys elegans (Didelphimorphia: Didelphidae): Establishing the phylogeographic determinants. Rev. Chil. Hist. Nat. 2021, 94, 5. [Google Scholar] [CrossRef]
  32. Ceccarelli, F.S.; Ojanguren-Affilastro, A.A.; Ramírez, M.J.; Ochoa, J.A.; Mattoni, C.I.; Prendini, L. Andean uplift drives diversification of the bothriurid scorpion genus Brachistosternus. J. Biogeogr. 2016, 43, 1942–1954. [Google Scholar] [CrossRef]
  33. Ortiz-Jaureguizar, E.; Cladera, G.A. Paleoenvironmental evolution of southern South America during the Cenozoic. J. Arid Environ. 2006, 66, 498–532. [Google Scholar] [CrossRef]
  34. Garzione, C.N.; McQuarrie, N.; Perez, N.D.; Ehlers, T.A.; Beck, S.L.; Kar, N.; Eichelberger, N.; Chapman, A.D.; Ward, K.M.; Ducea, M.N.; et al. Tectonic Evolution of the Central Andean Plateau and Implications for the Growth of Plateaus. Annu. Rev. Earth Planet. Sci. 2017, 45, 529–559. [Google Scholar] [CrossRef]
  35. Zachos, J.; Pagani, M.; Sloan, L.; Thomas, E.; Billups, K. Trends, Rhythms, and Aberrations in Global Climate 65 Ma to Present. Science 2001, 292, 686–693. [Google Scholar] [CrossRef]
  36. Glade-Vargas, N.S.; Rojas, C.; Jara-Arancio, P.; Vidal, P.; Arroyo, M.T.K.; Hinojosa, L.F. Biogeography of Argylia D. Don (Bignoniaceae): Diversification, Andean Uplift and Niche Conservatism. Front. Plant Sci. 2021, 12, 2270. [Google Scholar] [CrossRef]
  37. San Juan, E.; Araya-Donoso, R.; Véliz, D.; Quiroga, N.; Botto-Mahan, C. Genetic diversity in a restricted-dispersal kissing bug: The centre–periphery hypothesis halfway. Mol. Ecol. 2021, 30, 4660–4672. [Google Scholar] [CrossRef]
  38. Canals, M.; Cruzat, L.; Molina, M.C.; Ferreira, A.; Cattan, P.E. Blood host sources of Mepraia spinolai (Heteroptera: Reduviidae), wild vector of Chagas Disease in Chile. J. Med. Entomol. 2001, 38, 303–307. [Google Scholar] [CrossRef] [Green Version]
  39. Quiroga, N.; Correa, J.P.; Campos-Soto, R.; San Juan, E.; Araya-Donoso, R.; Díaz-Campusano, G.; González, C.R.; Botto-Mahan, C. Blood-Meal Sources and Trypanosoma cruzi Infection in Coastal and Insular Triatomine Bugs from the Atacama Desert of Chile. Microorganisms 2022, 10, 785. [Google Scholar] [CrossRef]
  40. Torres-Pérez, F.; Acuna-Retamar, M.; Cook, J.A.; Bacigalupo, A.; García, A.; Cattan, P.E. Statistical phylogeography of Chagas disease vector Triatoma infestans: Testing biogeographic hypotheses of dispersal. Infect. Genet. Evol. 2011, 11, 167–174. [Google Scholar] [CrossRef]
  41. Bacigalupo, A.; Torres-Pérez, F.; Segovia, V.; García, A.; Correa, J.P.; Moreno, L.; Arroyo, P.; Cattan, P.E. Sylvatic foci of the Chagas disease vector Triatoma infestans in Chile: Description of a new focus and challenges for control programs. Memórias do Instituto Oswaldo Cruz 2010, 105, 633–641. [Google Scholar] [CrossRef] [Green Version]
  42. Baranzelli, M.C.; Johnson, L.A.; Cosacov, A.; Sérsic, A.N. Historical and ecological divergence among populations of Monttea chilensis (Plantaginaceae), an endemic endangered shrub bordering the Atacama Desert, Chile. Evol. Ecol. 2014, 28, 751–774. [Google Scholar] [CrossRef]
  43. Mayle, F.E.; Burbridge, R.; Killeen, T.J. Millennial-Scale Dynamics of Southern Amazonian Rain Forests. Science 2000, 290, 2291–2294. [Google Scholar] [CrossRef] [Green Version]
  44. Dézerald, O.; Latorre, C.; Betancourt, J.L.; Brito Vera, G.A.; González, A.L. Ecological fidelity and spatiotemporal resolution of arthropod death assemblages from rodent middens in the central Atacama Desert (northern Chile). Quat. Sci. Rev. 2019, 210, 15–25. [Google Scholar] [CrossRef]
  45. Latorre, C.; Betancourt, J.L.; Rylander, K.A.; Quade, J. Vegetation invasions into absolute desert: A 45;th000 yr rodent midden record from the Calama–Salar de Atacama basins, northern Chile (lat 22°–24°S). GSA Bull. 2002, 114, 349–366. [Google Scholar] [CrossRef]
  46. Zúñiga-Reinoso, Á.; Pinto, P.; Predel, R. A New Species of Gyriosomus Guérinmeneville (Coleoptera: Tenebrionidae) from the Chilean Atacama Desert. Ann. Zool. 2019, 69, 105–112. [Google Scholar] [CrossRef]
  47. Guerrero, P.C.; Rosas, M.; Arroyo, M.T.K.; Wiens, J.J. Evolutionary lag times and recent origin of the biota of an ancient desert (Atacama–Sechura). Proc. Natl. Acad. Sci. USA 2013, 110, 11469–11474. [Google Scholar] [CrossRef] [Green Version]
  48. Rundel, P.; Dillon, M.; Palma, B.; Mooney, H.; Gulmon, S.; Ehleringer, J. The Phytogeography and Ecology of the Coastal Atacama and Peruvian Deserts. Aliso A J. Syst. Florist. Bot. 1991, 13, 1–49. [Google Scholar] [CrossRef] [Green Version]
  49. Garreaud, R.; Barichivich, J.; Christie, D.A.; Maldonado, A. Interannual variability of the coastal fog at Fray Jorge relict forests in semiarid Chile. J. Geophys. Res. 2008, 113, G04011. [Google Scholar] [CrossRef] [Green Version]
  50. Garreaud, R.; Molina, A.; Farias, M. Andean uplift, ocean cooling and Atacama hyperaridity: A climate modeling perspective. Earth Planet. Sci. Lett. 2010, 292, 39–50. [Google Scholar] [CrossRef]
  51. Pedoja, K.; Husson, L.; Johnson, M.E.; Melnick, D.; Witt, C.; Pochat, S.; Nexer, M.; Delcaillau, B.; Pinegina, T.; Poprawski, Y.; et al. Earth-Science Reviews Coastal staircase sequences reflecting sea-level oscillations and tectonic uplift during the Quaternary and Neogene. Earth Sci. Rev. 2014, 132, 13–38. [Google Scholar] [CrossRef]
  52. Ortlieb, L.; Zazo, C.; Goy, J.L.; Hillaire-Marcel, C.; Ghaleb, B.; Cournoyer, L. Coastal deformation and sea-level changes in the northern Chile subduction area (23 degrees S) during the last 330 ky. Quat. Sci. Rev. 1996, 15, 819–831. [Google Scholar] [CrossRef]
  53. Quezada, J.; González, G.; Dunai, T.; Jensen, A.; Juez-Larré, J. Alzamiento litoral Pleistoceno del norte de Chile: Edades 21Ne de la terraza costera más alta del área de Caldera-Bahía Inglesa. Rev. Geológica de Chile 2007, 34, 81–96. [Google Scholar] [CrossRef]
  54. Binnie, A.; Dunai, T.J.; Binnie, S.A.; Victor, P.; González, G.; Bolten, A. Accelerated late quaternary uplift revealed by 10Be exposure dating of marine terraces, Mejillones Peninsula, northern Chile. Quat. Geochronol. 2016, 36, 12–27. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Collection sites of Mepraia species with the biogeographic characteristics of the current habitats. (b) Sampling site of M. gajardoi in the extreme north of the coastal Atacama Desert (Caleta Vitor), photograph by Ricardo Campos-Soto. (c) Sampling site of M. sp. on the coast of the Mejillones Peninsula (Santa María Island) in the Atacama Desert, photograph by Ricardo Campos-Soto. (d) Habitat of M. parapatrica in the Pan de Azúcar National Park in the Atacama Desert, photograph by Ricardo Campos-Soto. (e) Sampling site of M. spinolai in the interior valley of central Chile in the Coquimban Biogeographical Province [17], photograph by Ricardo Campos-Soto. (f) Adult M. gajardoi, male on the left and female on the right, photographs by Vicente Valdés. (g) Female adult M. sp., photograph by Ricardo Campos-Soto (h) Adults of M. parapatrica, male on the left and female on the right, photographs by Vicente Valdés. (i) Adult M. spinolai, male on the left and female on the right, photographs by Vicente Valdés. The vertical black bar represents one centimetre.
Figure 1. (a) Collection sites of Mepraia species with the biogeographic characteristics of the current habitats. (b) Sampling site of M. gajardoi in the extreme north of the coastal Atacama Desert (Caleta Vitor), photograph by Ricardo Campos-Soto. (c) Sampling site of M. sp. on the coast of the Mejillones Peninsula (Santa María Island) in the Atacama Desert, photograph by Ricardo Campos-Soto. (d) Habitat of M. parapatrica in the Pan de Azúcar National Park in the Atacama Desert, photograph by Ricardo Campos-Soto. (e) Sampling site of M. spinolai in the interior valley of central Chile in the Coquimban Biogeographical Province [17], photograph by Ricardo Campos-Soto. (f) Adult M. gajardoi, male on the left and female on the right, photographs by Vicente Valdés. (g) Female adult M. sp., photograph by Ricardo Campos-Soto (h) Adults of M. parapatrica, male on the left and female on the right, photographs by Vicente Valdés. (i) Adult M. spinolai, male on the left and female on the right, photographs by Vicente Valdés. The vertical black bar represents one centimetre.
Insects 13 00419 g001
Figure 2. Phylogeographic hypotheses tested in this study. Scenario 1: hypothesis generated by the BMCMC proposed colonization from the Coquimban Biogeographical Province (Figure 1a) to the north, where the most southerly species M. spinolai is sister to the northern clade formed by M. gajardoi, M. parapatrica and M. sp, and, within this clade, M. sp. is sister to the remaining species M. gajardoi and M. parapatrica [9]. Scenario 2 represents the colonization hypothesis from north to south proposed by Frias-Lasserre (2010) [8]. Scenario 3: M. spinolai is sister to a northern clade that includes M. gajardoi as sister to the remaining species M. sp. and M. parapatrica. Scenario 4 suggests colonization from south to north following a parapatric speciation model, where M. spinolai is sister to a clade that includes M. parapatrica as sister to the remaining species M. gajardoi and M. sp. Scenario 5 shows a southern clade formed by M. spinolai as sister to M. sp. and a northern clade formed by M. gajardoi and M. parapatrica.
Figure 2. Phylogeographic hypotheses tested in this study. Scenario 1: hypothesis generated by the BMCMC proposed colonization from the Coquimban Biogeographical Province (Figure 1a) to the north, where the most southerly species M. spinolai is sister to the northern clade formed by M. gajardoi, M. parapatrica and M. sp, and, within this clade, M. sp. is sister to the remaining species M. gajardoi and M. parapatrica [9]. Scenario 2 represents the colonization hypothesis from north to south proposed by Frias-Lasserre (2010) [8]. Scenario 3: M. spinolai is sister to a northern clade that includes M. gajardoi as sister to the remaining species M. sp. and M. parapatrica. Scenario 4 suggests colonization from south to north following a parapatric speciation model, where M. spinolai is sister to a clade that includes M. parapatrica as sister to the remaining species M. gajardoi and M. sp. Scenario 5 shows a southern clade formed by M. spinolai as sister to M. sp. and a northern clade formed by M. gajardoi and M. parapatrica.
Insects 13 00419 g002
Figure 3. Most probable colonization scenario of Mepraia species inferred in BEAST and DIYABC using mitochondrial gene sequences (COI, cyt b and ND4). (a) Bayesian consensus tree with divergence times. The numbers at the nodes are the divergence times in millions of years; the black bars at the nodes show the 95% HPD for node age. The black circles at the nodes indicate node support above 0.95 pp. The colour coding of the branches and areas in the map corresponds to the species shown in Figure 1. The timing of the hyperarid cycles is shown in the yellowish shaded area on the tree according to [27] and as brown arrows on the map. (bg) illustrate the consecutive biogeographical events. The triangles on the maps illustrate the Andes Range after the last uplift. The red and blue dotted lines on map (b) delimit the Atacama and Coquimban Biogeographical Provinces, respectively, according to [17]. The black arrows on the maps indicate the colonization direction.
Figure 3. Most probable colonization scenario of Mepraia species inferred in BEAST and DIYABC using mitochondrial gene sequences (COI, cyt b and ND4). (a) Bayesian consensus tree with divergence times. The numbers at the nodes are the divergence times in millions of years; the black bars at the nodes show the 95% HPD for node age. The black circles at the nodes indicate node support above 0.95 pp. The colour coding of the branches and areas in the map corresponds to the species shown in Figure 1. The timing of the hyperarid cycles is shown in the yellowish shaded area on the tree according to [27] and as brown arrows on the map. (bg) illustrate the consecutive biogeographical events. The triangles on the maps illustrate the Andes Range after the last uplift. The red and blue dotted lines on map (b) delimit the Atacama and Coquimban Biogeographical Provinces, respectively, according to [17]. The black arrows on the maps indicate the colonization direction.
Insects 13 00419 g003
Figure 4. Graph of the posterior probabilities of each scenario estimated using a logistic regression approach.
Figure 4. Graph of the posterior probabilities of each scenario estimated using a logistic regression approach.
Insects 13 00419 g004
Table 1. Paired comparisons between different phylogenetic hypotheses using log (BF).
Table 1. Paired comparisons between different phylogenetic hypotheses using log (BF).
S1S2S3S4S5
S1−8647.238.5 *23 *18.1 *8.8 *
S2−38.5−8685.7−15.5−20.4−29.7
S3−2315.5−8670.2−4.9−14.2
S4−18.120.44.9−8665.3−9.3
S5−8.829.714.29.3−8656
S1–S5: different phylogenetic scenarios congruent with Figure 2. Italic numbers in the diagonal shows the log(ML) values for each scenario. Numbers with * indicate the best log(BF) values, with overwhelming support for scenario 1.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Campos-Soto, R.; Rodríguez-Valenzuela, E.; Díaz-Campusano, G.; Boric-Bargetto, D.; Zúñiga-Reinoso, Á.; Cianferoni, F.; Torres-Pérez, F. Testing Phylogeographic Hypotheses in Mepraia (Hemiptera: Reduviidae) Suggests a Complex Spatio-Temporal Colonization in the Coastal Atacama Desert. Insects 2022, 13, 419. https://doi.org/10.3390/insects13050419

AMA Style

Campos-Soto R, Rodríguez-Valenzuela E, Díaz-Campusano G, Boric-Bargetto D, Zúñiga-Reinoso Á, Cianferoni F, Torres-Pérez F. Testing Phylogeographic Hypotheses in Mepraia (Hemiptera: Reduviidae) Suggests a Complex Spatio-Temporal Colonization in the Coastal Atacama Desert. Insects. 2022; 13(5):419. https://doi.org/10.3390/insects13050419

Chicago/Turabian Style

Campos-Soto, Ricardo, Evelyn Rodríguez-Valenzuela, Gabriel Díaz-Campusano, Dusan Boric-Bargetto, Álvaro Zúñiga-Reinoso, Franco Cianferoni, and Fernando Torres-Pérez. 2022. "Testing Phylogeographic Hypotheses in Mepraia (Hemiptera: Reduviidae) Suggests a Complex Spatio-Temporal Colonization in the Coastal Atacama Desert" Insects 13, no. 5: 419. https://doi.org/10.3390/insects13050419

APA Style

Campos-Soto, R., Rodríguez-Valenzuela, E., Díaz-Campusano, G., Boric-Bargetto, D., Zúñiga-Reinoso, Á., Cianferoni, F., & Torres-Pérez, F. (2022). Testing Phylogeographic Hypotheses in Mepraia (Hemiptera: Reduviidae) Suggests a Complex Spatio-Temporal Colonization in the Coastal Atacama Desert. Insects, 13(5), 419. https://doi.org/10.3390/insects13050419

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop