Network-based hierarchical population structure analysis for large genomic data sets.

Analysis of population structure in natural populations using genetic data is a common practice in ecological and evolutionary studies. With large genomic data sets of populations now appearing more frequently across the taxonomic spectrum, it is becoming increasingly possible to reveal many hierarchical levels of structure, including fine-scale genetic clusters. To analyze these data sets, methods need to be appropriately suited to the challenges of extracting multilevel structure from whole-genome data. Here, we present a network-based approach for constructing population structure representations from genetic data. The use of community-detection algorithms from network theory generates a natural hierarchical perspective on the representation that the method produces. The method is computationally efficient, and it requires relatively few assumptions regarding the biological processes that underlie the data. We show the approach by analyzing population structure in the model plant species Arabidopsis thaliana and in human populations. These examples illustrate how network-based approaches for population structure analysis are well-suited to extracting valuable ecological and evolutionary information in the era of large genomic data sets.

Phenotypic plasticity and local adaptations to dissolved oxygen in larvae fire salamander (Salamandra infraimmaculata).

A key environmental factor that varies both spatially and temporally in surface waters is dissolved oxygen (DO). In stagnant ephemeral freshwater ponds, DO can fluctuate diurnally and seasonally, while the constant mixing of water in streams typically maintain DO levels close to saturation with only minor fluctuations. Larvae of the Near Eastern fire salamander (Salamandra infraimmaculata) develop in a range of waterbodies that vary in flow and permanence. To study inter-population variation in larval response to environmental change, we translocated larvae between stream and pond habitats and exposed larvae sampled from different habitat types to hypoxic and normoxic conditions in the laboratory. Larvae transferred from stream to pond retain gill size, while larvae transferred from pond to stream show a reduction in gill size. Larvae that were caged within their native habitat, either stream or pond, display a decrease in gill size similar to larvae transferred from pond to stream. When exposed to experimentally manipulated levels of DO in the laboratory larvae, respectively, increase and decrease gill size under hypoxic and normoxic conditions. Habitat-type origin had a significant effect on the degree of change in gill size with larvae from permanent streams demonstrating the lowest absolute variation in gill size. There was no interaction between DO level (hypoxic/normoxic) and the larvae habitat-type origin. These results suggest that S. infraimmaculata larvae are locally adapted to their aquatic breeding habitat through the plastic ability to respond to the prevailing respiratory conditions by rapidly decreasing or increasing gill size.

Revealing life-history traits by contrasting genetic estimations with predictions of effective population size.

Effective population size, a central concept in conservation biology, is now routinely estimated from genetic surveys and can also be theoretically predicted from demographic, life-history, and mating-system data. By evaluating the consistency of theoretical predictions with empirically estimated effective size, insights can be gained regarding life-history characteristics and the relative impact of different life-history traits on genetic drift. These insights can be used to design and inform management strategies aimed at increasing effective population size. We demonstrated this approach by addressing the conservation of a reintroduced population of Asiatic wild ass (Equus hemionus). We estimated the variance effective size (Nev ) from genetic data (N ev =24.3) and formulated predictions for the impacts on Nev of demography, polygyny, female variance in lifetime reproductive success (RS), and heritability of female RS. By contrasting the genetic estimation with theoretical predictions, we found that polygyny was the strongest factor affecting genetic drift because only when accounting for polygyny were predictions consistent with the genetically measured Nev . The comparison of effective-size estimation and predictions indicated that 10.6% of the males mated per generation when heritability of female RS was unaccounted for (polygyny responsible for 81% decrease in Nev ) and 19.5% mated when female RS was accounted for (polygyny responsible for 67% decrease in Nev ). Heritability of female RS also affected Nev ; hf2=0.91 (heritability responsible for 41% decrease in Nev ). The low effective size is of concern, and we suggest that management actions focus on factors identified as strongly affecting Nev, namely, increasing the availability of artificial water sources to increase number of dominant males contributing to the gene pool. This approach, evaluating life-history hypotheses in light of their impact on effective population size, and contrasting predictions with genetic measurements, is a general, applicable strategy that can be used to inform conservation practice.

Stochastic modelling of shifts in allele frequencies reveals a strongly polygynous mating system in the re-introduced Asiatic wild ass.

Small populations are prone to loss of genetic variation and hence to a reduction in their evolutionary potential. Therefore, studying the mating system of small populations and its potential effects on genetic drift and genetic diversity is of high importance for their viability assessments. The traditional method for studying genetic mating systems is paternity analysis. Yet, as small populations are often rare and elusive, the genetic data required for paternity analysis are frequently unavailable. The endangered Asiatic wild ass (Equus hemionus), like all equids, displays a behaviourally polygynous mating system; however, the level of polygyny has never been measured genetically in wild equids. Combining noninvasive genetic data with stochastic modelling of shifts in allele frequencies, we developed an alternative approach to paternity analysis for studying the genetic mating system of the re-introduced Asiatic wild ass in the Negev Desert, Israel. We compared the shifts in allele frequencies (as a measure of genetic drift) that have occurred in the wild ass population since re-introduction onset to simulated scenarios under different proportions of mating males. We revealed a strongly polygynous mating system in which less than 25% of all males participate in the mating process each generation. This strongly polygynous mating system and its potential effect on the re-introduced population's genetic diversity could have significant consequences for the long-term persistence of the population in the Negev. The stochastic modelling approach and the use of allele-frequency shifts can be further applied to systems that are affected by genetic drift and for which genetic data are limited.

Allele-specific network reveals combinatorial interaction that transcends small effects in psoriasis GWAS.

Hundreds of genetic markers have shown associations with various complex diseases, yet the "missing heritability" remains alarmingly elusive. Combinatorial interactions may account for a substantial portion of this missing heritability, but their discoveries have been impeded by computational complexity and genetic heterogeneity. We present BlocBuster, a novel systems-level approach that efficiently constructs genome-wide, allele-specific networks that accurately segregate homogenous combinations of genetic factors, tests the associations of these combinations with the given phenotype, and rigorously validates the results using a series of unbiased validation methods. BlocBuster employs a correlation measure that is customized for single nucleotide polymorphisms and returns a multi-faceted collection of values that captures genetic heterogeneity. We applied BlocBuster to analyze psoriasis, discovering a combinatorial pattern with an odds ratio of 3.64 and Bonferroni-corrected p-value of 5.01×10(-16). This pattern was replicated in independent data, reflecting robustness of the method. In addition to improving prediction of disease susceptibility and broadening our understanding of the pathogenesis underlying psoriasis, these results demonstrate BlocBuster's potential for discovering combinatorial genetic associations within heterogeneous genome-wide data, thereby transcending the limiting "small effects" produced by individual markers examined in isolation.

Landscape influences on dispersal behaviour: a theoretical model and empirical test using the fire salamander, Salamandra infraimmaculata.

When populations reside within a heterogeneous landscape, isolation by distance may not be a good predictor of genetic divergence if dispersal behaviour and therefore gene flow depend on landscape features. Commonly used approaches linking landscape features to gene flow include the least cost path (LCP), random walk (RW), and isolation by resistance (IBR) models. However, none of these models is likely to be the most appropriate for all species and in all environments. We compared the performance of LCP, RW and IBR models of dispersal with the aid of simulations conducted on artificially generated landscapes. We also applied each model to empirical data on the landscape genetics of the endangered fire salamander, Salamandra infraimmaculata, in northern Israel, where conservation planning requires an understanding of the dispersal corridors. Our simulations demonstrate that wide dispersal corridors of the low-cost environment facilitate dispersal in the IBR model, but inhibit dispersal in the RW model. In our empirical study, IBR explained the genetic divergence better than the LCP and RW models (partial Mantel correlation 0.413 for IBR, compared to 0.212 for LCP, and 0.340 for RW). Overall dispersal cost in salamanders was also well predicted by landscape feature slope steepness (76%), and elevation (24%). We conclude that fire salamander dispersal is well characterised by IBR predictions. Together with our simulation findings, these results indicate that wide dispersal corridors facilitate, rather than hinder, salamander dispersal. Comparison of genetic data to dispersal model outputs can be a useful technique in inferring dispersal behaviour from population genetic data.

Evolution of the human gastrokine locus and confounding factors regarding the pseudogenicity of GKN3.

In a screen for genes expressed specifically in gastric mucous neck cells, we identified GKN3, the recently discovered third member of the gastrokine family. We present confirmatory mouse data and novel porcine data showing that mouse GKN3 expression is confined to mucous cells of the corpus neck and antrum base and is prominently expressed in metaplastic lesions. GKN3 was proposed originally to be expressed in some human populations and a pseudogene in others. To investigate that hypothesis, we studied human GKN3 evolution in the context of its paralogous genomic neighbors, GKN1 and GKN2. Haplotype analysis revealed that GKN3 mimics GKN2 in patterns of exonic SNP allocation, whereas GKN1 appeared to be more stringently selected. GKN3 showed signatures of both directional selection and population based selective sweeps in humans. One such selective sweep includes SNP rs10187256, originally identified as an ancestral tryptophan to premature STOP codon mutation. The derived (nonancestral) allele went to fixation in Asia. We show that another SNP, rs75578132, identified 5 bp downstream of rs10187256, exhibits a second selective sweep in almost all Europeans, some Latinos, and some Africans, possibly resulting from a reintroduction of European genes during African colonization. Finally, we identify a mutation that would destroy the splice donor site in the putative exon3-intron3 boundary, which occurs in all human genomes examined to date. Our results highlight a stomach-specific human genetic locus, which has undergone various selective sweeps across European, Asian, and African populations and thus reflects geographic and ethnic patterns in genome evolution.

Genetic restoration in the eastern collared lizard under prescribed woodland burning.

Eastern collared lizards of the Ozarks live in glades--open, rocky habitats embedded in a woodland matrix. Past fire suppression had made the woodlands a barrier to dispersal, leading to habitat destruction, fragmentation and local extinction. Reintroduced populations of lizards were subjected to 10 years of habitat fragmentation under continued fire suppression followed by twelve years of landscape restoration with prescribed burns. Prior to prescribed burning, genetic diversity decreased within glades and differentiation increased among glades. With woodland burning, genetic diversity within glades first decreased during an expanding colonization phase, but then increased as a dynamically stable metapopulation was established. Population differentiation among glades also stabilized in the metapopulation under weak isolation-by-distance. This study is one of the first to examine the genetic changes in a species of conservation concern throughout all the stages of decline and recovery and shows the importance of landscape-level restoration for maintaining the genetic integrity of populations. This study also demonstrates how mark-recapture and genetic data together can yield detailed insight into metapopulation dynamics that would be impossible from just one type of data alone.

Coherent and incoherent inference in phylogeography and human evolution.

A hypothesis is nested within a more general hypothesis when it is a special case of the more general hypothesis. Composite hypotheses consist of more than one component, and in many cases different composite hypotheses can share some but not all of these components and hence are overlapping. In statistics, coherent measures of fit of nested and overlapping composite hypotheses are technically those measures that are consistent with the constraints of formal logic. For example, the probability of the nested special case must be less than or equal to the probability of the general model within which the special case is nested. Any statistic that assigns greater probability to the special case is said to be incoherent. An example of incoherence is shown in human evolution, for which the approximate Bayesian computation (ABC) method assigned a probability to a model of human evolution that was a thousand-fold larger than a more general model within which the first model was fully nested. Possible causes of this incoherence are identified, and corrections and restrictions are suggested to make ABC and similar methods coherent. Another coalescent-based method, nested clade phylogeographic analysis, is coherent and also allows the testing of individual components of composite hypotheses, another attribute lacking in ABC and other coalescent-simulation approaches. Incoherence is a highly undesirable property because it means that the inference is mathematically incorrect and formally illogical, and the published incoherent inferences on human evolution that favor the out-of-Africa replacement hypothesis have no statistical or logical validity.

Structural variants identified using non-Mendelian inheritance patterns advance the mechanistic understanding of autism spectrum disorder.

The heritability of autism spectrum disorder (ASD), based on 680,000 families and five countries, is estimated to be nearly 80%, yet heritability reported from SNP-based studies are consistently lower, and few significant loci have been identified with genome-wide association studies. This gap in genomic information may reside in rare variants, interaction among variants (epistasis), or cryptic structural variation (SV) and may provide mechanisms that underlie ASD. Here we use a method to identify potential SVs based on non-Mendelian inheritance patterns in pedigrees using parent-child genotypes from ASD families and demonstrate that they are enriched in ASD-risk genes. Most are in non-coding genic space and are over-represented in expression quantitative trait loci, suggesting that they affect gene regulation, which we confirm with their overlap of differentially expressed genes in postmortem brain tissue of ASD individuals. We then identify an SV in the GRIK2 gene that alters RNA splicing and a regulatory region of the ACMSD gene in the kynurenine pathway as significantly associated with a non-verbal ASD phenotype, supporting our hypothesis that these currently excluded loci can provide a clearer mechanistic understanding of ASD. Finally, we use an explainable artificial intelligence approach to define subgroups demonstrating their use in the context of precision medicine.

Plasma proteomics of SARS-CoV-2 infection and severity reveals impact on Alzheimer's and coronary disease pathways.

Identification of proteins dysregulated by COVID-19 infection is critically important for better understanding of its pathophysiology, building prognostic models, and identifying new targets. Plasma proteomic profiling of 4,301 proteins was performed in two independent datasets and tested for the association for three COVID-19 outcomes (infection, ventilation, and death). We identified 1,449 proteins consistently associated in both datasets with any of these three outcomes. We subsequently created highly accurate models that distinctively predict infection, ventilation, and death. These proteins were enriched in specific biological processes including cytokine signaling, Alzheimer's disease, and coronary artery disease. Mendelian randomization and gene network analyses identified eight causal proteins and 141 highly connected hub proteins including 35 with known drug targets. Our findings provide distinctive prognostic biomarkers for two severe COVID-19 outcomes, reveal their relationship to Alzheimer's disease and coronary artery disease, and identify potential therapeutic targets for COVID-19 outcomes.

Plasma proteomics of SARS-CoV-2 infection and severity reveals impact on Alzheimer and coronary disease pathways.

Identification of the plasma proteomic changes of Coronavirus disease 2019 (COVID-19) is essential to understanding the pathophysiology of the disease and developing predictive models and novel therapeutics. We performed plasma deep proteomic profiling from 332 COVID-19 patients and 150 controls and pursued replication in an independent cohort (297 cases and 76 controls) to find potential biomarkers and causal proteins for three COVID-19 outcomes (infection, ventilation, and death). We identified and replicated 1,449 proteins associated with any of the three outcomes (841 for infection, 833 for ventilation, and 253 for death) that can be query on a web portal ( https://covid.proteomics.wustl.edu/ ). Using those proteins and machine learning approached we created and validated specific prediction models for ventilation (AUC>0.91), death (AUC>0.95) and either outcome (AUC>0.80). These proteins were also enriched in specific biological processes, including immune and cytokine signaling (FDR ≤ 3.72×10 -14 ), Alzheimer's disease (FDR ≤ 5.46×10 -10 ) and coronary artery disease (FDR ≤ 4.64×10 -2 ). Mendelian randomization using pQTL as instrumental variants nominated BCAT2 and GOLM1 as a causal proteins for COVID-19. Causal gene network analyses identified 141 highly connected key proteins, of which 35 have known drug targets with FDA-approved compounds. Our findings provide distinctive prognostic biomarkers for two severe COVID-19 outcomes (ventilation and death), reveal their relationship to Alzheimer's disease and coronary artery disease, and identify potential therapeutic targets for COVID-19 outcomes.

The transition from isolated patches to a metapopulation in the eastern collared lizard in response to prescribed fires.

Habitat fragmentation often arises from human-induced alterations to the matrix that reduce or eliminate dispersal between habitat patches. Elimination of dispersal increases local extinction and decreases recolonization. These phenomena were observed in the eastern collared lizard (Crotaphytus collaris collaris), which lives in the mid-continental highland region of the Ozarks (Missouri, USA) on glades: habitats of exposed bedrock that form desert-like habitats imbedded in a woodland matrix. With the onset of woodland fire suppression, glade habitats degenerated and the woodland matrix was altered to create a strong barrier to dispersal. By 1980, lizard populations in the Ozarks were rapidly going extinct. In response to this decline, some glades were restored by clearing and burning. Starting in 1984, collared lizard populations were translocated onto these restored habitats. The translocated populations persisted but did not colonize nearby glades or disperse among one another. In 1994 prescribed woodland fires were initiated, which unleashed much dispersal and colonizing behavior. Dispersal was highly nonrandom by both intrinsic variables (age, gender) and extrinsic variables (overall demography, glade population sizes, glade areas, landscape features), resulting in different classes of lizards being dominant in creating demographic cohesiveness among glades, colonizing new glades on a mountain, and colonizing new mountain systems. A dramatic transition was documented from isolated fragments, to a nonequilibrium colonizing metapopulation, and finally to a stable metapopulation. This transition is characterized by the convergence of rates of extinction and recolonization and a major alteration of dispersal probabilities and pattern in going from the nonequilibrium to stable metapopulation states.

Mitochondrial bioenergetics as a major motive force of speciation.

Mitochondrial bioenergetics plays a key role in multiple basic cellular processes, such as energy production, nucleotide biosynthesis, and iron metabolism. It is an essential system for animals' life and death (apoptosis) and it is required for embryo development. This, in conjunction with its being subjected to adaptive processes in multiple species and its gene products being involved in the formation of reproductive barriers in animals, raises the possibility that mitochondrial bioenergetics could be a candidate genetic mechanism of speciation. Here, we discuss genetic and biochemical evidence for the possible involvement of this unique system, encoded by two genomes (the mitochondrial and nuclear genomes), that differ by an order of magnitude in their mutation rates in processes leading to speciation events.

What Determines Paternity in Wild Lizards? A Spatiotemporal Analysis of Behavior and Morphology.

Mating behavior in animals can be understood as a sequence of events that begins with individuals encountering one another and ends with the production of offspring. Behavioral descriptions of animal interactions characterize early elements of this sequence, and genetic descriptions use offspring parentage to characterize the final outcome, with behavioral and physiological assessments of mates and mechanisms of copulation and fertilization comprising intermediate steps. However, behavioral and genetic descriptions of mating systems are often inconsistent with one another, complicating expectations for crucial aspects of mating biology, such as the presence of multiple mating. Here, we use behavioral and genetic data from a wild population of the lizard Anolis cristatellus to characterize female multiple mating and the potential for sexual selection through female mate choice in this species. We find that 48% of sampled females bore offspring sired by multiple males. Moreover, spatiotemporal proximity between males and females was associated with whether a male sired a female's offspring, and if yes, how many offspring he sired. Additionally, male body size, but not display behavior, was associated with reproductive outcomes for male-female pairs. While much remains to be learned about the mechanisms of mating and targets of sexual selection in A. cristatellus, it is clear that female multiple mating is a substantial component of this species' mating system in nature.

Subspecies hybridization as a potential conservation tool in species reintroductions.

Reintroductions are a powerful tool for the recovery of endangered species. However, their long-term success is strongly influenced by the genetic diversity of the reintroduced population. The chances of population persistence can be improved by enhancing the population's adaptive ability through the mixing of individuals from different sources. However, where source populations are too diverse the reintroduced population could also suffer from outbreeding depression or unsuccessful admixture due to behavioural or genetic barriers. For the reintroduction of Asiatic wild ass Equus hemionus ssp. in Israel, a breeding core was created from individuals of two different subspecies (E. h. onager & E. h. kulan). Today the population comprises approximately 300 individuals and displays no signs of outbreeding depression. The aim of this study was a population genomic evaluation of this conservation reintroduction protocol. We used maximum likelihood methods and genetic clustering analyses to investigate subspecies admixture and test for spatial autocorrelation based on subspecies ancestry. Further, we analysed heterozygosity and effective population sizes in the breeding core prior to release and the current wild population. We discovered high levels of subspecies admixture in the breeding core and wild population, consistent with a significant heterozygote excess in the breeding core. Furthermore, we found no signs of spatial autocorrelation associated with subspecies ancestry in the wild population. Inbreeding and variance effective population size estimates were low. Our results indicate no genetic or behavioural barriers to admixture between the subspecies and suggest that their hybridization has led to greater genetic diversity in the reintroduced population. The study provides rare empirical evidence of the successful application of subspecies hybridization in a reintroduction. It supports use of intraspecific hybridization as a tool to increase genetic diversity in conservation translocations.

How frugal is Mother Nature with haplotypes?

MOTIVATION: Inference of haplotypes from genotype data is crucial and challenging for many vitally important studies. The first, and most critical step, is the ascertainment of a biologically sound model to be optimized. Many models that have been proposed rely partially or entirely on reducing the number of unique haplotypes in the solution. RESULTS: This article examines the parsimony of haplotypes using known haplotypes as well as genotypes from the HapMap project. Our study reveals that there are relatively few unique haplotypes, but not always the least possible, for the datasets with known solutions. Furthermore, we show that there are frequently very large numbers of parsimonious solutions, and the number increases exponentially with increasing cardinality. Moreover, these solutions are quite varied, most of which are not consistent with the true solutions. These results quantify the limitations of the Pure Parsimony model and demonstrate the imperative need to consider additional properties for haplotype inference models. At a higher level, and with broad applicability, this article illustrates the power of combinatorial methods to tease out imperfections in a given biological model.

The reality and importance of founder speciation in evolution.

A founder event occurs when a new population is established from a small number of individuals drawn from a large ancestral population. Mayr proposed that genetic drift in an isolated founder population could alter the selective forces in an epistatic system, an observation supported by recent studies. Carson argued that a period of relaxed selection could occur when a founder population is in an open ecological niche, allowing rapid population growth after the founder event. Selectable genetic variation can actually increase during this founder-flush phase due to recombination, enhanced survival of advantageous mutations, and the conversion of non-additive genetic variance into additive variance in an epistatic system, another empirically confirmed prediction. Templeton combined the theories of Mayr and Carson with population genetic models to predict the conditions under which founder events can contribute to speciation, and these predictions are strongly confirmed by the empirical literature. Much of the criticism of founder speciation is based upon equating founder speciation to an adaptive peak shift opposed by selection. However, Mayr, Carson and Templeton all modeled a positive interaction of selection and drift, and Templeton showed that founder speciation is incompatible with peak-shift conditions. Although rare, founder speciation can have a disproportionate importance in adaptive innovation and radiation, and examples are given to show that "rare" does not mean "unimportant" in evolution. Founder speciation also interacts with other speciation mechanisms such that a speciation event is not a one-dimensional process due to either selection alone or drift alone.

The diverse applications of cladistic analysis of molecular evolution, with special reference to nested clade analysis.

The genetic variation found in small regions of the genomes of many species can be arranged into haplotype trees that reflect the evolutionary genealogy of the DNA lineages found in that region and the accumulation of mutations on those lineages. This review demonstrates some of the many ways in which clades (branches) of haplotype trees have been applied in recent years, including the study of genotype/phenotype associations at candidate loci and in genome-wide association studies, the phylogeographic history of species, human evolution, the conservation of endangered species, and the identification of species.

Phylogeography of the common vampire bat (Desmodus rotundus): marked population structure, Neotropical Pleistocene vicariance and incongruence between nuclear and mtDNA markers.

BACKGROUND: The common vampire bat Desmodus rotundus is an excellent model organism for studying ecological vicariance in the Neotropics due to its broad geographic range and its preference for forested areas as roosting sites. With the objective of testing for Pleistocene ecological vicariance, we sequenced a mitocondrial DNA (mtDNA) marker and two nuclear markers (RAG2 and DRB) to try to understand how Pleistocene glaciations affected the distribution of intraspecific lineages in this bat. RESULTS: Five reciprocally monophyletic clades were evident in the mitochondrial gene tree, and in most cases with high bootstrap support: Central America (CA), Amazon and Cerrado (AMC), Pantanal (PAN), Northern Atlantic Forest (NAF) and Southern Atlantic Forest (SAF). The Atlantic forest clades formed a monophyletic clade with high bootstrap support, creating an east/west division for this species in South America. On the one hand, all coalescent and non-coalescent estimates point to a Pleistocene time of divergence between the clades. On the other hand, the nuclear markers showed extensive sharing of haplotypes between distant localities, a result compatible with male-biased gene flow. In order to test if the disparity between the mitochondrial and nuclear markers was due to the difference in mutation rate and effective size, we performed a coalescent simulation to examine the feasibility that, given the time of separation between the observed lineages, even with a gene flow rate close to zero, there would not be reciprocal monophyly for a neutral nuclear marker. We used the observed values of theta and an estimated mutation rate for the nuclear marker gene to perform 1000 iterations of the simulation. The results of this simulation were inconclusive: the number of iterations with and without reciprocal monophyly of one or more clades are similar. CONCLUSIONS: We therefore conclude that the pattern exhibited by the common vampire bat, with marked geographical structure for a mitochondrial marker and no phylogeographic structure for nuclear markers is compatible with a historical scenario of complete isolation of refuge-like populations during the Pleistocene. The results on demographic history on this species is compatible with the Carnaval-Moritz model of Pleistocene vicariance, with demographic expansions in the southern Atlantic forest.