The automation and evaluation of nested clade phylogeographic analysis

Nested clade phylogeographic analysis (NCPA) is a popular method for reconstructing the demographic history of spatially distributed populations from genetic data. Although some parts of the analysis are automated, there is no unique and widely followed algorithm for doing this in its entirety, beginning with the data, and ending with the inferences drawn from the data. This article describes a method that automates NCPA, thereby providing a framework for replicating analyses in an objective way. To do so, a number of decisions need to be made so that the automated implementation is representative of previous analyses. We review how the NCPA procedure has evolved since its inception and conclude that there is scope for some variability in the manual application of NCPA. We apply the automated software to three published datasets previously analyzed manually and replicate many details of the manual analyses, suggesting that the current algorithm is representative of how a typical user will perform NCPA. We simulate a large number of replicate datasets for geographically distributed, but entirely random-mating, populations. These are then analyzed using the automated NCPA algorithm. Results indicate that NCPA tends to give a high frequency of false positives. In our simulations we observe that 14% of the clades give a conclusive inference that a demographic event has occurred, and that 75% of the datasets have at least one clade that gives such an inference. This is mainly due to the generation of multiple statistics per clade, of which only one is required to be significant to apply the inference key. We survey the inferences that have been made in recent publications and show that the most commonly inferred processes (restricted gene flow with isolation by distance and contiguous range expansion) are those that are commonly inferred in our simulations. However, published datasets typically yield a richer set of inferences with NCPA than obtained in our random-mating simulations, and further testing of NCPA with models of structured populations is necessary to examine its accuracy.

The automation of nested clade phylogeographic analysis

ANeCA is a fully automated implementation of Nested Clade Phylogeographic Analysis. This was originally developed by Templeton and colleagues, and has been used to infer, from the pattern of gene sequence polymorphisms in a geographically structured population, the historical demographic processes that have shaped its evolution. Until now it has been necessary to perform large parts of the procedure manually. We provide a program that will take data in Nexus sequential format, and directly output a set of inferences. The software also includes TCS v1.18 and GeoDis v2.2 as part of automation.

The phylogeography and connectivity of the latitudinally widespread scleractinian coral Plesiastrea versipora in the Western Pacific

Whereas terrestrial animal populations might show genetic connectivity within a continent, marine species, such as hermatypic corals, may have connectivity stretching to all corners of the planet. We quantified the genetic variability within and among populations of the widespread scleractinian coral, Plesiastrea versipora along the eastern Australian seaboard (4145 km) and the Ryukyu Archipelago (Japan, 681 km) using sequences of internal transcribed spacers (ITS1-2) from ribosomal DNA. Geographic patterns in genetic variability were deduced from a nested clade analysis (NCA) performed on a parsimony network haplotype. This analysis allowed the establishment of geographical associations in the distribution of haplotypes within the network cladogram, therefore allowing us to deduce phylogeographical patterns based under models of restricted gene flow, fragmentation and range expansion. No significant structure was found among Ryukyu Archipelago populations. The lack of an association between the positions of haplotypes in the cladogram with geographical location of these populations may be accounted for by a high level of gene flow of P. versipora within this region, probably due to the strong Kuroshio Current. In contrast, strong geographical associations were apparent among populations of P. versipora along the south-east coast of Australia. This pattern of restricted genetic connectivity among populations of P. versipora on the eastern seaboard of Australia seems to be associated with the present surface ocean current (the East Australian Current) on this side of the south-western Pacific Ocean.

Genetic variability of Aedes aegypti in the Americas using a mitochondrial gene: evidence of multiple introductions

To analyze the genetic relatedness and phylogeographic structure of Aedes aegypti, we collected samples from 36 localities throughout the Americas (Brazil, Peru, Venezuela, Guatemala, US), three from Africa (Guinea, Senegal, Uganda), and three from Asia (Singapore, Cambodia, Tahiti). Amplification and sequencing of a fragment of the mitochondrial NADH dehydrogenase subunit 4 gene identified 20 distinct haplotypes, of which 14 are exclusive to the Americas, four to African/Asian countries, one is common to the Americas and Africa, and one to the Americas and Asia. Nested clade analysis (NCA), pairwise distribution, statistical parsimony, and maximum parsimony analyses were used to infer evolutionary and historic processes, and to estimate phylogenetic relationships among haplotypes. Two clusters were found in all the analyses. Haplotypes clustered in the two clades were separated by eight mutational steps. Phylogeographic structure detected by the NCA was consistent with distant colonization within one clade and fragmentation followed by range expansion via long distance dispersal in the other. Three percent of nucleotide divergence between these two clades is suggestive of a gene pool division that may support the hypothesis of occurrence of two subspecies of Ae. aegypti in the Americas.

Phylogeography of Mabuya maculilabris (Reptilia) from São Tomé Island (Gulf of Guinea) inferred from mtDNA sequences

The pattern of genetic variation of the lizard Mabuya maculilabris from São Tomé Island (Gulf of Guinea) was investigated using a combination of three mitochondrial DNA gene fragments. Forty-eight haplotypes were recovered among 66 individuals covering the whole island. The genealogy inferred from the most parsimonious network of haplotypes allows us to detect two main and long branches departing from the putative group of oldest haplotypes. The tips of these branches exhibit star-like phylogenies, which may indicate of recently expanded populations, most probably from a small number of founders. A nested clade analysis suggests a complex pattern of past events that gave rise to the extant geographical pattern found in the haplotype distribution: past and allopatric fragmentation, range expansion, restricted gene Xow and long-distance dispersal. These results are consistent with the complex geological history of the island where important volcanic activity with extensive lava Xows has occurred during several periods. Mismatch- distribution analysis and AMOVA also support these conclusions. Substantial genetic structuring among these lizards was detected as well as high levels of diVerentiation between the southern edge populations (particularly those from the Rolas Islet) and the remaining ones. However, variation is low relative to the geological age of the island. Our results indicate that patterns of variation observed in reptiles in other oceanic islands are not indicative of those observed in the islands of the Gulf of Guinea.

Phylogeographic analysis of the chloroplast DNA variation in wild common bean (Phaseolus vulgaris L.) in the Americas

The wild common bean (Phaseolus vulgaris) is widely but discontinuously distributed from northern Mexico to northern Argentina on both sides of the Isthmus of Panama. Little is known on how the species has reached its current disjunct distribution. In this research, chloroplast DNA polymorphisms in seven non-coding regions were used to study the history of migration of wild P. vulgaris between Mesoamerica and South America. A penalized likelihood analysis was applied to previously published Leguminosae ITS data to estimate divergence times between P. vulgaris and its sister taxa from Mesoamerica, and divergence times of populations within P. vulgaris. Fourteen chloroplast haplotypes were identified by PCR-RFLP and their geographical associations were studied by means of a Nested Clade Analysis and Mantel Tests. The results suggest that the haplotypes are not randomly distributed but occupy discrete parts of the geographic range of the species. The current distribution of haplotypes may be explained by isolation by distance and by at least two migration events between Mesoamerica and South America: one from Mesoamerica to South America and another one from northern South America to Mesoamerica. Age estimates place the divergence of P. vulgaris from its sister taxa from Mesoamerica at or before 1.3 Ma, and divergence of populations from Ecuador-northern Peru at or before 0.6 Ma. As these ages are taken as minimum divergence times, the influence of past events, such as the closure of the Isthmus of Panama and the final uplift of the Andes, on the migration history and population structure of this species cannot be disregarded.

Statistical inferences in phylogeography

In conventional phylogeographic studies, historical demographic processes are elucidated from the geographical distribution of individuals represented on an inferred gene tree. However, the interpretation of gene trees in this context can be difficult as the same demographic/geographical process can randomly lead to multiple different genealogies. Likewise, the same gene trees can arise under different demographic models. This problem has led to the emergence of many statistical methods for making phylogeographic inferences. A popular phylogeographic approach based on nested clade analysis is challenged by the fact that a certain amount of the interpretation of the data is left to the subjective choices of the user, and it has been argued that the method performs poorly in simulation studies. More rigorous statistical methods based on coalescence theory have been developed. However, these methods may also be challenged by computational problems or poor model choice. In this review, we will describe the development of statistical methods in phylogeographic analysis, and discuss some of the challenges facing these methods.

Genus Paracoccidioides: Species Recognition and Biogeographic Aspects

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Phylogeography of the Solanaceae-infecting Basidiomycota fungus Rhizoctonia solani AG-3 based on sequence analysis of two nuclear DNA loci

Background: the soil fungus Rhizoctonia solani anastomosis group 3 (AG-3) is an important pathogen of cultivated plants in the family Solanaceae. Isolates of R. solani AG-3 are taxonomically related based on the composition of cellular fatty acids, phylogenetic analysis of nuclear ribosomal DNA (rDNA) and beta-tubulin gene sequences, and somatic hyphal interactions. Despite the close genetic relationship among isolates of R. solani AG-3, field populations from potato and tobacco exhibit comparative differences in their disease biology, dispersal ecology, host specialization, genetic diversity and population structure. However, little information is available on how field populations of R. solani AG-3 on potato and tobacco are shaped by population genetic processes. In this study, two field populations of R. solani AG-3 from potato in North Carolina (NC) and the Northern USA; and two field populations from tobacco in NC and Southern Brazil were examined using sequence analysis of two cloned regions of nuclear DNA (pP42F and pP89).Results: Populations of R. solani AG-3 from potato were genetically diverse with a high frequency of heterozygosity, while limited or no genetic diversity was observed within the highly homozygous tobacco populations from NC and Brazil. Except for one isolate (TBR24), all NC and Brazilian isolates from tobacco shared the same alleles. No alleles were shared between potato and tobacco populations of R. solani AG-3, indicating no gene flow between them. To infer historical events that influenced current geographical patterns observed for populations of R. solani AG-3 from potato, we performed an analysis of molecular variance (AMOVA) and a nested clade analysis (NCA). Population differentiation was detected for locus pP89 (Phi(ST) = 0.257, significant at P < 0.05) but not for locus pP42F (Phi(ST) = 0.034, not significant). Results based on NCA of the pP89 locus suggest that historical restricted gene flow is a plausible explanation for the geographical association of clades. Coalescent-based simulations of genealogical relationships between populations of R. solani AG-3 from potato and tobacco were used to estimate the amount and directionality of historical migration patterns in time, and the ages of mutations of populations. Low rates of historical movement of genes were observed between the potato and tobacco populations of R. solani AG-3.Conclusion: the two sisters populations of the basidiomycete fungus R. solani AG-3 from potato and tobacco represent two genetically distinct and historically divergent lineages that have probably evolved within the range of their particular related Solanaceae hosts as sympatric species.

Phylogeography of Hypostomus strigaticeps (Siluriformes: Loricariidae) inferred by mitochondrial DNA reveals its distribution in the upper Paraná River basin

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Genetic characterisation of populations of the critically endangered Goliath grouper (Epinephelus itajara, Serranidae) from the Northern Brazilian coast through analyses of mtDNA

The Goliath grouper (Epinephelus itajara) is one of the most endangered species of fish of the subfamily Epinephelinae. Slow to develop and mature, and dependent on mangrove habitats for breeding, the species also suffers intense harvesting, which has reduced drastically in numbers in many areas. To contribute to the understanding of the characteristics of E. itajara populations, we conducted a molecular genetics study of the species, focusing on populations from the Northern Brazilian coast. The mtDNA control region (D-loop) of 116 individuals from five localities (Bragança, Ajuruteua, Parnaíba, Fortaleza and Natal) was analysed, and a sequence of 499 base pairs identified. Analyses of the sequences indicated that genetic variability was generally lower in E. itajara than in other endangered species of the genus. AMOVA found no significant grouping structure among the populations. Nested Clade Analysis revealed a significant association between genetic variability and geographic distribution among only three populations (Ajuruteua, Parnaíba and Natal). Genetic diversity was higher in populations from the Amazon region, which may be related to the better conservation of mangrove habitats in this area. Therefore, the present study could be used for the implementation of conservation and management measures in order to protect and consolidate these populations.

Filogeografia de Nasutitermes corniger (Isoptera: Termitidae: Nasutitermitinae) na região Neotropical

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Molecular phylogeography of the Woodland Ringlet (Erebia medusa [Denis and Schiffermüller] 1775) in Europe

Phylogeography is a recent field of biological research that links phylogenetics to biogeography through deciphering the imprint that evolutionary history has left on the genetic structure of extant populations. During the cold phases of the successive ice ages, which drastically shaped species’ distributions since the Pliocene, populations of numerous species were isolated in refugia where many of them evolved into different genetic lineages. My dissertation deals with the phylogeography of the Woodland Ringlet (Erebia medusa [Denis and Schiffermüller] 1775) in Central and Eastern Europe. This Palaearctic butterfly species is currently distributed from central France and south eastern Belgium over large parts of Central Europe and southern Siberia to the Pacific. It is absent from those parts of Europe with mediterranean, oceanic and boreal climates. It was supposed to be a Siberian faunal element with a rather homogeneous population structure in Central Europe due to its postglacial expansion out of a single eastern refugium. An already existing evolutionary scenario for the Woodland Ringlet in Central and Eastern Europe is based on nuclear data (allozymes). To know if this is corroborated by organelle evolutionary history, I sequenced two mitochondrial markers (part of the cytochrome oxydase subunit one and the control region) for populations sampled over the same area. Phylogeography largely relies on the construction of networks of uniparentally inherited haplotypes that are compared to geographic haplotype distribution thanks to recent developed methods such as nested clade phylogeographic analysis (NCPA). Several ring-shaped ambiguities (loops) emerged from both haplotype networks in E. medusa. They can be attributed to recombination and homoplasy. Such loops usually avert the straightforward extraction of the phylogeographic signal contained in a gene tree. I developed several new approaches to extract phylogeographic information in the presence of loops, considering either homoplasy or recombination. This allowed me to deduce a consistent evolutionary history for the species from the mitochondrial data and also adds plausibility for the occurrence of recombination in E. medusa mitochondria. Despite the fact that the control region is assumed to have a lack of resolving power in other species, I found a considerable genetic variation of this marker in E. medusa which makes it a useful tool for phylogeographic studies. In combination with the allozyme data, the mitochondrial genome supports the following phylogeographic scenario for E. medusa in Europe: (i) a first vicariance, due to the onset of the Würm glaciation, led to the formation of several major lineages, and is mirrored in the NCPA by restricted gene flow, (ii) later on further vicariances led to the formation of two sub-lineages in the Western lineage and two sub-lineages in the Eastern lineage during the Last Glacial Maximum or Older Dryas; additionally the NCPA supports a restriction of gene flow with isolation by distance, (iii) finally, vicariance resulted in two secondary sub-lineages in the area of Germany and, maybe, to two other secondary sub-lineages in the Czech Republic. The last postglacial warming was accompanied by strong range expansions in most of the genetic lineages. The scenario expected for a presumably Siberian faunal element such as E. medusa is a continuous loss of genetic diversity during postglacial westward expansion. Hence, the pattern found in this thesis contradicts a typical Siberian origin of E. medusa. In contrast, it corroboratess the importance of multiple extra-Mediterranean refugia for European fauna as it was recently assumed for other continental species.

Phylogeography of an east Australian wet-forest bird, the satin bowerbird (Ptilonorhynchus violaceus), derived from mtDNA, and its relationship to morphology

Australian wet forests have undergone a contraction in range since the mid-Tertiary, resulting in a fragmented distribution along the east Australian coast incorporating several biogeographical barriers. Variation in mitochondrial DNA and morphology within the satin bowerbird was used to examine biogeographical structure throughout almost the entire geographical extent of these wet forest fragments. We used several genetic analysis techniques, nested clade and barrier analyses, that use patterns inherent in the data to describe the spatial structuring. We also examined the validity of the two previously described satin bowerbird subspecies that are separated by well-defined biogeographical barriers and tested existing hypotheses that propose divergence occurs within each subspecies across two other barriers, the Black Mountain corridor and the Hunter Valley. Our data showed that the two subspecies were genetically and morphologically divergent. The northern subspecies, found in the Wet Tropics region of Queensland, showed little divergence across the Black Mountain corridor, a barrier found to be significant in other Wet Tropics species. Biogeographical structure was found through southeastern Australia; three geographically isolated populations showed genetic differentiation, although minimal divergence was found across the proposed Hunter Valley barrier. A novel barrier was found separating inland and coastal populations in southern New South Wales. Little morphological divergence was observed within subspecies, bar a trend for birds to be larger in the more southerly parts of the species' range. The results from both novel and well-established genetic analyses were similar, providing greater confidence in the conclusions about spatial divergence and supporting the validity of these new techniques.