Statistical phylogeography: methods of evaluating and minimizing inference errors

Population genetics of the naked mole-rat Heterocephalus glaber: The role of rivers in shaping genetic structure

We used nested clade phylogeographic analysis (NCPA) of mitochondrial DNA sequence data to examine the processes contributing to population structure in naked mole-rats. We examined sequence variation in the (1097 bp) control region D-loop of the mitochondrial genome in 303 individuals from 174 colonies of naked mole-rats (Heterocephalus glaber) located mainly within an 870 km2 area in Meru National Park, Kenya. Four rivers were found to be correlated to a significant fragmentation inference in the NCPA. The largest pairwise divergence between haplotypes from populations separated by rivers was 2.74%, which was well over half of the divergence reported between the extremes of the distribution from southern Ethiopia to southern Kenya (4.6%). However, the size of the river (measured in current discharge) was not a good predictor of the amount of sequence divergence between populations separated by a river. It appears that a large-scale historical fragmentation event may have conflated fragmentation patterns on a smaller scale, when recent colonization and range expansion brought two old lineages together at a location with a relatively small river separating them.

Phylogeny of Drosophila saltans group (Diptera: Drosophilidae) based on morphological and molecular evidence

Drosophila saltans group belongs to the subgenus Sophophora (family Drosophilidae), and it is subdivided into five subgroups, with 23 species. The species in this group are widely distributed in the Americas, primarily in the Neotropics. In the literature, the phylogenetic reconstruction of this group has been performed with various markers, but many inconsistencies remain. Here, we present a phylogenetic reconstruction of the saltans group with a greater number of species, 16 species, which is the most complete to date for the saltans group and includes all subgroups, in a combined analysis with morphological and molecular markers. We incorporated 48 morphological characters of male terminalia, the highest number used to date, and molecular markers based on mitochondrial genes COI and COII. Based on the results, which have recovered the five subgroups as distinct lineages, we propose a new hypothesis regarding the phylogenetic relationships among the subgroups of the saltans group. The relationships of the species within the sturtevanti and elliptica subgroups were well supported. The saltans subgroup showed several polytomies, but the relationship between the sibling species D. austrosaltans and D. saltans and their close relation with D. nigrosaltans were well supported in the molecular and total evidence analyses. The morphological analysis additionally supported the formation of the clade D. nigrosaltans—D. pseudosaltans. The observed polytomies may represent synchronous radiations or have resulted from speciation rates that have been too fast relative to the pace of substitution accumulation.

Phylogeography of Cedros and Tiburón Island Mule Deer in North America's Desert Southwest

Though mule deer (Odocoileus hemionus) persist in robust populations throughout most of their North American distribution, habitat loss, unregulated hunting, and other factors have reduced their historical range in México. Two of the 6 putative subspecies inhabiting México's deserts and Baja California peninsula are of conservation concern, occupying islands in the Pacific Ocean (Odocoileus hemionus cerrosensis on Cedros Island: endangered) and Sea of Cortés (Odocoileus hemionus sheldoni on Tiburón Island: threatened). Focusing on the desert southwest (n = 448), we sampled Tiburón (n = 22) and Cedros (n = 15) Island mule deer using contemporary samples and natural history museum specimens to complete a phylogeographic evaluation of the species complex, and assess the phylogeography of these insular subspecies. Both insular subspecies formed endemic haplotype lineages, consistent with island biogeographic theory. Bayesian skyline plots were consistent with Holocene demographic expansion. Cedros Island deer were genetically most similar to adjacent mainland Baja California deer, but exhibited a suite of unique haplotypes and reduced genetic variation. Tiburón Island deer haplotypes unexpectedly nested within a mainland lineage found in distant New Mexico, rather than the adjacent mainland Sonoran lineage. Such findings suggest the importance of postglacial climate fluctuations and biotic community turnover in the phylogeographic history of mule deer in the desert southwest. Our genetic data corroborates cultural, archaeological, and phenotypic evidence supporting Cedros and Tiburón deer endemicity and subspecies status. Reduced genetic variation, divergence from mainland populations, and demographic trends on both islands indicate that conservation, monitoring, and management are critical to ensure persistence of these endemic insular subspecies.

Biological and biochemical diversity in different biotypes of spotted stem borer, Chilo partellus (Swinhoe) in India

Because of variation in incidence and severity of damage by Chilo partellus (Swinhoe) in different geographical regions, it is difficult to identify stable sources of resistance against this pest. Therefore, the present studies were undertaken on biological attributes (damage in resistant and susceptible genotypes, survival and development) and biochemical profiles (amino acids and lipophilic compound) of C. partellus populations from eight geographical regions to understand it's population structure in India. There was a significant variation in biological attributes and biochemical profiles of C. partellus populations from different geographical regions. Based on virulence and biological attributes, similarity index placed the C. partellus populations in five groups. Likewise, lipophilic and amino acid profiling also placed the C. partellus populations in five groups. However, the different clusters based on biological and biochemical attributes did not include populations from the same regions. Similarity index based on virulence, biological attributes, and amino acids and lipophilic profiles placed the C. partellus populations in six groups. The C. partellus populations from Hisar, Hyderabad, Parbhani and Coimbatore were distinct from each other, indicating that there are four biotypes of C. partellus in India. The results suggested that sorghum and maize genotypes need to be tested against these four populations to identify stable sources of resistance. However, there is a need for further studies to establish the restriction in gene flow through molecular approaches across geographical regions to establish the distinctiveness of different biotypes of C. partellus in India.

Waves Out of the Korean Peninsula and Inter- and Intra-Species Replacements in Freshwater Fishes in Japan

The Japanese archipelago is located at the periphery of the continent of Asia. Rivers in the Japanese archipelago, separated from the continent of Asia by about 17 Ma, have experienced an intermittent exchange of freshwater fish taxa through a narrow land bridge generated by lowered sea level. As the Korean Peninsula and Japanese archipelago were not covered by an ice sheet during glacial periods, phylogeographical analyses in this region can trace the history of biota that were, for a long time, beyond the last glacial maximum. In this study, we analyzed the phylogeography of four freshwater fish taxa, Hemibarbus longirostris, dark chub Nipponocypris temminckii, Tanakia ssp. and Carassius ssp., whose distributions include both the Korean Peninsula and Western Japan. We found for each taxon that a small component of diverse Korean clades of freshwater fishes migrated in waves into the Japanese archipelago to form the current phylogeographic structure of biota. The replacements of indigenous populations by succeeding migrants may have also influenced the phylogeography.

Evolutionary origins of taro (Colocasia esculenta) in Southeast Asia

As an ancient clonal root and leaf crop, taro (Colocasia esculenta, Araceae) is highly polymorphic with uncertain genetic and geographic origins. We explored chloroplast DNA diversity in cultivated and wild taros, and closely related wild taxa, and found cultivated taro to be polyphyletic, with tropical and temperate clades that appear to originate in Southeast Asia sensu lato. A third clade was found exclusively in wild populations from Southeast Asia to Australia and Papua New Guinea. Our findings do not support the hypothesis of taro domestication in Papua New Guinea, despite archaeological evidence for early use or cultivation there, and the presence of apparently natural wild populations in the region (Australia and Papua New Guinea).

Landscape Genetics of Plants: Challenges and Opportunities

Dispersal is one of the most important but least understood processes in plant ecology and evolutionary biology. Dispersal of seeds maintains and establishes populations, and pollen and seed dispersal are responsible for gene flow within and among populations. Traditional views of dispersal and gene flow assume models that are governed solely by geographic distance and do not account for variation in dispersal vector behavior in response to heterogenous landscapes. Landscape genetics integrates population genetics with Geographic Information Systems (GIS) to evaluate the effects of landscape features on gene flow patterns (effective dispersal). Surprisingly, relatively few landscape genetic studies have been conducted on plants. Plants present advantages because their populations are stationary, allowing more reliable estimates of the effects of landscape features on effective dispersal rates. On the other hand, plant dispersal is intrinsically complex because it depends on the habitat preferences of the plant and its pollen and seed dispersal vectors. We discuss strategies to assess the separate contributions of pollen and seed movement to effective dispersal and to delineate the effects of plant habitat quality from those of landscape features that affect vector behavior. Preliminary analyses of seed dispersal for three species indicate that isolation by landscape resistance is a better predictor of the rates and patterns of dispersal than geographic distance. Rates of effective dispersal are lower in areas of high plant habitat quality, which may be due to the effects of the shape of the dispersal kernel or to movement behaviors of biotic vectors. Landscape genetic studies in plants have the potential to provide novel insights into the process of gene flow among populations and to improve our understanding of the behavior of biotic and abiotic dispersal vectors in response to heterogeneous landscapes.

Phylogeography of Daphnia magna Straus (Crustacea: Cladocera) in Northern Eurasia: Evidence for a deep longitudinal split between mitochondrial lineages

Species with a large geographic distributions present a challenge for phylogeographic studies due to logistic difficulties of obtaining adequate sampling. For instance, in most species with a Holarctic distribution, the majority of studies has concentrated on the European or North American part of the distribution, with the Eastern Palearctic region being notably understudied. Here, we study the phylogeography of the freshwater cladoceran Daphnia magna Straus, 1820 (Crustacea: Cladocera), based on partial mitochondrial COI sequences and using specimens from populations spread longitudinally from westernmost Europe to easternmost Asia, with many samples from previously strongly understudied regions in Siberia and Eastern Asia. The results confirm the previously suspected deep split between Eastern and Western mitochondrial haplotype super-clades. We find a narrow contact zone between these two super-clades in the eastern part of Western Siberia, with proven co-occurrence in a single lake in the Novosibirsk region. However, at present there is no evidence suggesting that the two mitochondrial super-clades represent cryptic species. Rather, they may be explained by secondary contact after expansion from different refugia. Interestingly, Central Siberia has previously been found to be an important contact zone also in other cladoceran species, and may thus be a crucial area for understanding the Eurasian phylogeography of freshwater invertebrates. Together, our study provides an unprecedented complete, while still not global, picture of the phylogeography of this important model species.

Phylogeography of Nasutitermes corniger (Isoptera: Termitidae) in the Neotropical Region

BACKGROUND

The Neotropical Region is known for its biodiversity and ranks third in number of known termite species. However, biogeographic and phylogeographic information of termites of this region is limited compared to other world geographic regions. Nasutitermes corniger is widely distributed in the region and is of considerable economic importance. The goal of this study was to describe the phylogeography of N. corniger in the Neotropical Region, to better understand its evolutionary processes.

RESULTS

The sampled populations of N. corniger showed high genetic variation. Results indicated strong geographic structure among N. corniger populations, with most haplotypes not broadly shared among separated locations. Phylogeographic analyses showed a dispersal route for N. corniger from Central America into South America via the Isthmus of Panama, with subsequent dispersal through the highlands east of the Andes and into eastern regions of the continent. The majority of haplotypes were limited in distribution to proximal regions, corresponding to particular biomes (Atlantic Forest, Amazonia, Chaco, Cerrado and Caatinga).

CONCLUSIONS

Nasutitermes corniger is suggested to be a good model for biogeographic and phylogeographic studies in the Neotropical Region. This study clarified the phylogeographic history of N. corniger and can contribute to the understanding of biogeographic dispersion processes in the Neotropical Region.

Molecular epidemiology of infectious diseases

Molecular epidemiology (ME) is a branch of epidemiology developed by merging molecular biology into epidemiological studies. In this paper, the authors try to discuss the ways that molecular epidemiology studies identify infectious diseases' causation and pathogenesis, and unravel infectious agents' sources, reservoirs, circulation pattern, transmission pattern, transmission probability, and transmission order. They bring real-world examples of research works in each area to make each study design more understandable. They also address some research areas and study design aspects that need further attention in future. They close with some thoughts about future directions in this field and emphasize on the need for training competent molecular epidemiology specialists that are capable of dealing with rapid advances in the field.

Phylogeography, genetic diversity and population structure in a Patagonian endemic plant

Phylogeographical patterns provide valuable insight into the historical processes underlying diversification, and may provide a better understanding of biodiversity, dispersal modes, diversification times, extinctions, refuge areas and other species-/population-level processes. Here, we examine the genetic structure of Oxalis laciniata, a representative of Oxalis, which is an important emerging model in evolutionary biology and phylogenetic studies. We analyse genetic correlation, haplotype diversity and genetic structure. In this study, cpDNA reveals the presence of 16 haplotypes, connected in tree networks. Genetic diversity is high and polymorphism is low among populations based on ISSR markers. Both clustering and analysis of the structure of the population indicate two different groups. Distribution modelling predicts two potential distribution areas. Our main conclusions are: (i) The phylogeographical pattern demonstrates non-random organization of genetic variability since two distinct groups can be distinguished; (ii) two refugia are proposed: one is situated in the SE, holding the most ancestral haplotypes; and the second one is situated in the SW; (iii) we propose an in situ diversification hypothesis for the populations located in the steppe; (iv) the centre of diversification coincides with the centre of the distribution; (v) distribution modelling shows a strong correspondence with the distribution of the species but it also suggests the possibility of occurrence in the Central Andes.

Islands within an island: Population genetic structure of the endemic Sardinian newt, Euproctus platycephalus

The identification of historic and contemporary barriers to dispersal is central to the conservation of endangered amphibians, but may be hindered by their complex life history and elusive nature. The complementary information generated by mitochondrial (mtDNA) and microsatellite markers generates a valuable tool in elucidating population structure and the impact of habitat fragmentation. We applied this approach to the study of an endangered montane newt, Euproctus platycephalus. Endemic to the Mediterranean island of Sardinia, it is threatened by anthropogenic activity, disease, and climate change. We have demonstrated a clear hierarchy of structure across genetically divergent and spatially distinct subpopulations. Divergence between three main mountain regions dominated genetic partitioning with both markers. Mitochondrial phylogeography revealed a deep division dating to ca. 1 million years ago (Mya), isolating the northern region, and further differentiation between the central and southern regions ca. 0.5 Mya, suggesting an association with Pleistocene severe glacial oscillations. Our findings are consistent with a model of southward range expansion during glacial periods, with postglacial range retraction to montane habitat and subsequent genetic isolation. Microsatellite markers revealed further strong population structure, demonstrating significant divergence within the central region, and partial differentiation within the south. The northern population showed reduced genetic diversity. Discordance between mitochondrial and microsatellite markers at this scale indicated a further complexity of population structure, in keeping with male-biased dispersal and female philopatry. Our study underscores the need to elucidate cryptic population structure in the ecology and conservation strategies for endangered island-restricted amphibians, especially in the context of disease and climate change.

Demographic Structure and Evolutionary History of Drosophila ornatifrons (Diptera, Drosophilidae) from Atlantic Forest of Southern Brazil

Drosoph1la ornatifrons of the guarani group (Diptera: Drosophilidae) is found mainly in humid areas of the Atlantic Forest biome, especially in the southern region of Brazil. Historical and contemporary fragmentation events influenced species diversity and distribution in this biome, although the role of paleoclimatic and paleogeographic events remain to be verified. The objective of the present study was to evaluate the demographic structure of D. ornatifrons from collection sites that are remnants of Atlantic Forest in southern Brazil, in order to contribute to the understanding of the processes that affected the patterns of genetic variability in this species. To achieve this goal, we sequenced 51 individuals from nine localities and 64 individuals from six localities for the mitochondrial genes Cytochrome Oxidase I and II, respectively. Our results indicate that D. ornatifrons may have experienced a demographic expansion event from the southernmost locations of its distribution, most likely from those located next to the coast and in fragments of Atlantic Forest inserted in the Pampa biome (South 2 group), towards the interior (South 1 group). This expansion probably started after the last glacial maximum, between 20,000 and 18,000 years ago, and was intensified near the Pleistocene-Holocene transition, around 12,000 years ago, when temperature started to rise. In this work we discuss how the haplotypes found barriers to gene flow and dispersal, influenced by the biogeographic pattern of Atlantic Forest.

Phylogeography and palaeodistribution modelling of Nassauvia subgenus Strongyloma (Asteraceae): exploring phylogeographical scenarios in the Patagonian steppe

The Patagonian steppe is an immense, cold, arid region, yet phylogeographically understudied. Nassauvia subgen. Strongyloma is a characteristic element of the steppe, exhibiting a continuum of morphological variation. This taxon provides a relevant phylogeographical model not only to understand how past environmental changes shaped the genetic structure of its populations, but also to explore phylogeographical scenarios at the large geographical scale of the Patagonian steppe. Here, we (1) assess demographic processes and historical events that shaped current geographic patterns of haplotypic diversity; (2) analyze hypotheses of isolation in refugia, fragmentation of populations, and/or colonization of available areas during Pleistocene glaciations; and (3) model extant and palaeoclimatic distributions to support inferred phylogeographical patterns. Chloroplast intergenic spacers, rpl32–trnL and trnQ–5′rps16, were sequenced for 372 individuals from 63 populations. Nested clade analysis, analyses of molecular variance, and neutrality tests were performed to assess genetic structure and range expansion. The present potential distribution was modelled and projected onto a last glacial maximum (LGM) model. Of 41 haplotypes observed, ten were shared among populations associated with different morphological variants. Populations with highest haplotype diversity and private haplotypes were found in central-western and south-eastern Patagonia, consistent with long-term persistence in refugia during Pleistocene. Palaeomodelling suggested a shift toward the palaeoseashore during LGM; new available areas over the exposed Atlantic submarine platform were colonized during glaciations with postglacial retraction of populations. A scenario of fragmentation and posterior range expansion may explain the observed patterns in the center of the steppe, which is supported by palaeomodelling. Northern Patagonian populations were isolated from southern populations by the Chubut and the Deseado river basins during glaciations. Pleistocene glaciations indirectly impacted the distribution, demography, and diversification of subgen. Strongyloma through decreased winter temperatures and water availability in different areas of its range.

Host generalists and specialists emerging side by side: an analysis of evolutionary patterns in the cosmopolitan chewing louse genus Menacanthus

Parasites with wide host spectra provide opportunities to study the ecological parameters of speciation, as well as the process of the evolution of host specificity. The speciose and cosmopolitan louse genus Menacanthus comprises both multi-host and specialised species, allowing exploration of the ecological and historical factors affecting the evolution of parasites using a comparative approach. We used phylogenetic analysis to reconstruct evolutionary relationships in 14 species of Menacanthus based on the sequences of one mitochondrial and one nuclear gene. The results allowed us to validate species identification based on morphology, as well as to explore host distribution by assumed generalist and specialist species. Our analyses confirmed a narrow host use for several species, however in some cases, the supposed host specialists had a wider host spectrum than anticipated. In one case a host generalist (Menacanthus eurysternus) was clustered terminally on a clade almost exclusively containing host specialists. Such a clade topology indicates that the process of host specialisation may not be irreversible in parasite evolution. Finally, we compared patterns of population genetic structure, geographic distribution and host spectra between two selected species, M. eurysternus and Menacanthus camelinus, using haplotype networks. Menacanthus camelinus showed limited geographical distribution in combination with monoxenous host use, whereas M. eurysternus showed a global distribution and lack of host specificity. It is suggested that frequent host switching maintains gene flow between M. eurysternus populations on unrelated hosts in local populations. However, gene flow between geographically distant localities was restricted, suggesting that geography rather than host-specificity is the main factor defining the global genetic diversity of M. eurysternus.

Adaptation to ephemeral habitat may overcome natural barriers and severe habitat fragmentation in a fire-dependent species, the Bachman's Sparrow (Peucaea aestivalis)

Bachman's Sparrow (Peucaea aestivalis) is a fire-dependent species that has undergone range-wide population declines in recent decades. We examined genetic diversity in Bachman's Sparrows to determine whether natural barriers have led to distinct population units and to assess the effect of anthropogenic habitat loss and fragmentation. Genetic diversity was examined across the geographic range by genotyping 226 individuals at 18 microsatellite loci and sequencing 48 individuals at mitochondrial and nuclear genes. Multiple analyses consistently demonstrated little genetic structure and high levels of genetic variation, suggesting that populations are panmictic. Based on these genetic data, separate management units/subspecies designations or translocations to promote gene flow among fragmented populations do not appear to be necessary. Panmixia in Bachman's Sparrow may be a consequence of an historical range expansion and retraction. Alternatively, high vagility in Bachman's Sparrow may be an adaptation to the ephemeral, fire-mediated habitat that this species prefers. In recent times, high vagility also appears to have offset inbreeding and loss of genetic diversity in highly fragmented habitat.

Phylogeography of Cephalotaxus oliveri (Cephalotaxaceae) in relation to habitat heterogeneity, physical barriers and the uplift of the Yungui Plateau

Habitat heterogeneity, physical barriers, and the uplift of the Yungui Plateau were found to deeply affect the phylogeographic pattern and evolutionary history of Cephalotaxus oliveri, a perennial conifer endemic to China. In this study, we explored the phylogeography using three chloroplast sequences (trnL-trnF, trnT-trnD and atpB-rbcL) in 22 natural populations of C. oliveri distributed throughout its range. The Yungui Plateau populations of C. oliveri were revealed to origin ca. 9.15Ma by molecular clock estimation, which is consistent with rapid uplift of the Qinghai-Tibetan Plateau (QTP) ca. 8-10Ma. Additionally, geological effects of the Yungui Plateau were suggested to promote the rapid intra-specific differentiation of C. oliveri in the Pliocene and Early Pleistocene. The relatively low level of genetic diversity (h=0.719, θ=1.17×10(-3)) and high population differentiation (NST=0.771 and GST=0.642) implied restricted gene flow among populations, which was confirmed by the Nested Clade Analysis (NCA). Mismatch distribution and haplotypes network provided evidences of recent demographic population expansion. Furthermore, the statistical dispersal-vicariance analysis indicated that the center of origin was in Central China. The comparison of haplotype distribution patterns indicated that the regions of HNHPS and HBLD were the potential refugia during the Pleistocene ice ages. Our results highlighted that habitat heterogeneity and physical barriers presenting in a species range can predict genetic patterns.

Biological races in humans

Races may exist in humans in a cultural sense, but biological concepts of race are needed to access their reality in a non-species-specific manner and to see if cultural categories correspond to biological categories within humans. Modern biological concepts of race can be implemented objectively with molecular genetic data through hypothesis-testing. Genetic data sets are used to see if biological races exist in humans and in our closest evolutionary relative, the chimpanzee. Using the two most commonly used biological concepts of race, chimpanzees are indeed subdivided into races but humans are not. Adaptive traits, such as skin color, have frequently been used to define races in humans, but such adaptive traits reflect the underlying environmental factor to which they are adaptive and not overall genetic differentiation, and different adaptive traits define discordant groups. There are no objective criteria for choosing one adaptive trait over another to define race. As a consequence, adaptive traits do not define races in humans. Much of the recent scientific literature on human evolution portrays human populations as separate branches on an evolutionary tree. A tree-like structure among humans has been falsified whenever tested, so this practice is scientifically indefensible. It is also socially irresponsible as these pictorial representations of human evolution have more impact on the general public than nuanced phrases in the text of a scientific paper. Humans have much genetic diversity, but the vast majority of this diversity reflects individual uniqueness and not race.

Quaternary range and demographic expansion of Liolaemus darwinii (Squamata: Liolaemidae) in the Monte Desert of Central Argentina using Bayesian phylogeography and ecological niche modelling

Until recently, most phylogeographic approaches have been unable to distinguish between demographic and range expansion processes, making it difficult to test for the possibility of range expansion without population growth and vice versa. In this study, we applied a Bayesian phylogeographic approach to reconstruct both demographic and range expansion in the lizard Liolaemus darwinii of the Monte Desert in Central Argentina, during the Late Quaternary. Based on analysis of 14 anonymous nuclear loci and the cytochrome b mitochondrial DNA gene, we detected signals of demographic expansion starting at ~55 ka based on Bayesian Skyline and Skyride Plots. In contrast, Bayesian relaxed models of spatial diffusion suggested that range expansion occurred only between ~95 and 55 ka, and more recently, diffusion rates were very low during demographic expansion. The possibility of population growth without substantial range expansion could account for the shared patterns of demographic expansion during the Last Glacial Maxima (OIS 2 and 4) in fish, small mammals and other lizards of the Monte Desert. We found substantial variation in diffusion rates over time, and very high rates during the range expansion phase, consistent with a rapidly advancing expansion front towards the southeast shown by palaeo-distribution models. Furthermore, the estimated diffusion rates are congruent with observed dispersal rates of lizards in field conditions and therefore provide additional confidence to the temporal scale of inferred phylogeographic patterns. Our study highlights how the integration of phylogeography with palaeo-distribution models can shed light on both demographic and range expansion processes and their potential causes.

Influence of parameter settings in automated scoring of AFLPs on population genetic analysis

The use of procedures for the automated scoring of amplified fragment length polymorphisms (AFLP) fragments has recently increased. Corresponding software does not only automatically score the presence or absence of AFLP fragments, but also allows an evaluation of how different settings of scoring parameters influence subsequent population genetic analyses. In this study, we used the automated scoring package rawgeno to evaluate how five scoring parameters influence the number of polymorphic bins and estimates of pairwise genetic differentiation between populations (F(st)). Steps were implemented in r to automatically run the scoring process in rawgeno for a set of different parameter combinations. While we found the scoring parameters minimum bin width and minimum number of samples per bin to have only weak influence on pairwise F(st) values, maximum bin width and bin reproducibility had much stronger effects. The minimum average bin fluorescence scoring parameter affected F(st) values in an only moderate way. At a range of scoring parameters around the default settings of rawgeno, the number of polymorphic bins as well as pairwise F(st) values stayed rather constant. This study thus shows the particularities of AFLP scoring, be it either manual or automatical, can have profound effects on subsequent population genetic analysis.

Genetic diversity shaped by historical and recent factors in the live-bearing twoline skiffia Neotoca bilineata

The endangered twoline skiffia Neotoca bilineata, a viviparous fish of the subfamily Goodeinae, endemic to central Mexico (inhabiting two basins, Cuitzeo and Lerma-Santiago) was evaluated using genetic and habitat information. The genetic variation of all remaining populations of the species was analysed using both mitochondrial and microsatellite markers and their habitat conditions were assessed using a water quality index (I(WQ)). An 80% local extinction was found across the distribution of N. bilineata. The species was found in three of the 16 historical localities plus one previously unreported site. Most areas inhabited by the remaining populations had I(WQ) scores unsuitable for the conservation of freshwater biodiversity. Populations showed low but significant genetic differentiation with both markers (mtDNA φ(ST) = 0.076, P < 0.001; microsatellite F(ST) = 0.314, P < 0.001). Borbollon, in the Cuitzeo Basin, showed the highest level of differentiation and was identified as a single genetic unit by Bayesian assignment methods. Rio Grande de Morelia and Salamanca populations showed the highest genetic diversity and also a high migration rate facilitated by an artificial channel that connected the two basins. Overall, high genetic diversity values were observed compared with other freshwater fishes (average N(a) = 16 alleles and loci and mean ±S.D. H(o) = 0.63 ± 0.10 and nucleotide diversity π = 0.006). This suggests that the observed genetic diversity has not diminished as rapidly as the species' habitat destruction. No evidence of correlation between habitat conditions and genetic diversity was found. The current pattern of genetic diversity may be the result of both historical factors and recent modifications of the hydrological system. The main threat to the species may be the rapid habitat deterioration and associated demographic stochasticity rather than genetic factors.

Phylogeny and biogeography of the Anderson's crocodile newt, Echinotriton andersoni (Amphibia: Caudata), as revealed by mitochondrial DNA sequences

The Anderson's crocodile newt, Echinotriton andersoni, is considered a relic and endangered species distributed in the Central Ryukyus. To elucidate phylogenetic relationships and detailed genetic structures among populations, we analyzed variation in the mitochondrial cytochrome b gene. Results strongly support a primary dichotomy between populations from the Amami and Okinawa Island Groups with substantial genetic divergence, favoring a primary divergence between the two island groups. Within the latter, populations from the southern part of Okinawajima Island are shown to be more closely related to those from Tokashikijima Island than to those from the northern and central parts of Okinawajima. The prominent genetic divergence between the two island groups of the Central Ryukyus seems to have initiated in the Miocene, i.e., prior to formation of the strait that has consistently separated these island groups since the Pleistocene. The ancestor of the southern Okinawajima-Tokashikijima is estimated to have migrated from the northern and central parts of Okinawajima into southern Okinawajima at the Pleistocene, and dispersed into Tokashikijima subsequently.

Leve diferenciación genética entre los límites occidental y oriental de distribución de Astroides calycularis (Pallas, 1776) (Anthozoa, Scleractinia, Dendrophylliidae), inferida a partir de secuencias de COI e ITS

El estudio de la estructura de las poblaciones y su diferenciación a nivel genético es de gran utilidad para la elaboración de planes de manejo y conservación de especies amenazadas. En este estudio, utilizamos marcadores nucleares y mitocondriales (espaciadores internos de genes ribosomales -ITS y citocromo oxidasa, subunidad I -COI) y métodos de análisis filogenéticos y de clados anidados (NCA), para realizar la primera valoración de la estructura genética del coral naranja Astroides calycularis (Pallas, 1766), una especie amenazada del Mediterráneo, a partir de muestras de 12 localidades a lo largo de su área de distribución. En las localidades situadas en la región más occidental del Mediterráneo se encontró cierta homogeneidad genética, mientras que al comparar estas localidades con las de las cuencas argelina y del mar Tirreno se observó una ligera diferenciación.

Phylogeography of Camellia taliensis (Theaceae) inferred from chloroplast and nuclear DNA: insights into evolutionary history and conservation

Background

As one of the most important but seriously endangered wild relatives of the cultivated tea, Camellia taliensis harbors valuable gene resources for tea tree improvement in the future. The knowledge of genetic variation and population structure may provide insights into evolutionary history and germplasm conservation of the species.

Results

Here, we sampled 21 natural populations from the species' range in China and performed the phylogeography of C. taliensis by using the nuclear PAL gene fragment and chloroplast rpl32-trnL intergenic spacer. Levels of haplotype diversity and nucleotide diversity detected at rpl32-trnL (h = 0.841; π = 0.00314) were almost as high as at PAL (h = 0.836; π = 0.00417). Significant chloroplast DNA population subdivision was detected (G_ST = 0.988; N_ST = 0.989), suggesting fairly high genetic differentiation and low levels of recurrent gene flow through seeds among populations. Nested clade phylogeographic analysis of chlorotypes suggests that population genetic structure in C. taliensis has been affected by habitat fragmentation in the past. However, the detection of a moderate nrDNA population subdivision (G_ST = 0.222; N_ST = 0.301) provided the evidence of efficient pollen-mediated gene flow among populations and significant phylogeographical structure (N_ST > G_ST; P < 0.01). The analysis of PAL haplotypes indicates that phylogeographical pattern of nrDNA haplotypes might be caused by restricted gene flow with isolation by distance, which was also supported by Mantel’s test of nrDNA haplotypes (r = 0.234, P < 0.001). We found that chlorotype C1 was fixed in seven populations of Lancang River Region, implying that the Lancang River might have provided a corridor for the long-distance dispersal of the species.

Conclusions

We found that C. taliensis showed fairly high genetic differentiation resulting from restricted gene flow and habitat fragmentation. This phylogeographical study gives us deep insights into population structure of the species and conservation strategies for germplasm sampling and developing in situ conservation of natural populations.

Conflicting evolutionary patterns due to mitochondrial introgression and multilocus phylogeography of the Patagonian freshwater crab Aegla neuquensis

BACKGROUND

Multiple loci and population genetic methods were employed to study the phylogeographic history of the Patagonian freshwater crab Aegla neuquensis (Aeglidae: Decopoda). This taxon occurs in two large river systems in the Patagonian Steppe, from the foothills of the Andes Mountains east to the Atlantic Ocean.

METHODOLOGY/PRINCIPAL FINDINGS

A nuclear phylogeny and multilocus nested clade phylogeographic analysis detected a fragmentation event between the Negro and Chico-Chubut river systems. This event occurred approximately 137 thousand years ago. An isolation-with-migration analysis and maximum-likelihood estimates of gene flow showed asymmetrical exchange of genetic material between these two river systems exclusively in their headwaters. We used information theory to determine the best-fit demographic history between these two river systems under an isolation-with-migration model. The best-fit model suggests that the Negro and the ancestral populations have the same effective population sizes; whereas the Chico-Chubut population is smaller and shows that gene flow from the Chico-Chubut into the Negro is four times higher than in the reverse direction. Much of the Chico-Chubut system appears to have only been recently colonized while the Negro populations appear to have been in place for most of the evolutionary history of this taxon.

CONCLUSIONS/SIGNIFICANCE

Due to mitochondrial introgression, three nuclear loci provided different phylogeographic resolution than the three mitochondrial genes for an ancient fragmentation event observed in the nuclear phylogeny. However, the mitochondrial locus provided greater resolution on more recent evolutionary events. Our study, therefore, demonstrates the need to include both nuclear and mitochondrial loci for a more complete understanding of evolutionary histories and associated phylogeographic events. Our results suggest that gene flow between these systems, before and after fragmentation was through periodic paleolakes that formed in the headwaters region. Fragmentation between the Negro and Chico-Chubut systems was driven by the disappearance of these paleolakes during the Patagonian Glaciation.

AFLPMax: a user-friendly application for computing the optimal number of amplified fragment length polymorphism markers needed in phylogenetic reconstruction

Amplified fragment length polymorphisms (AFLPs) are widely used for phylogenetic inference especially in non-model species. Frequently, trees obtained with other nuclear or mitochondrial markers or with morphological information need additional resolution, increased branch support, or independent data sources (i.e. unlinked loci). In such cases, the use of AFLPs is a quick and cheap option. Computer simulation has shown that dominant AFLP markers lead to less accurate tree topologies than bi-allelic codominant markers such as SNPs, but this difference becomes negligible for shallow trees when using AFLP data sets that include a sufficiently large number of characters. Thus, determining how many AFLP characters are required to recover a given phylogeny is a key issue regarding the appropriateness of AFLPs for phylogenetic reconstruction. Here, we present a user-friendly, java-based graphical interface, AFLPMax, which executes an automatic pipeline of different programs providing the user with the optimal number of AFLP characters needed to recover a given phylogeny with high accuracy and support. Executables for Windows, linux and MacOS X operating systems, source code and user manual are available from: http://webs.uvigo.es/acraaj/AFLPMax.htm.

Intraspecific variation in Podostemum ceratophyllum (Podostemaceae): evidence of refugia and colonization since the last glacial maximum

PREMISE OF THE STUDY

Intraspecific variation among 20 populations of Podostemum ceratophyllum Michx. was investigated to test the hypothesis of range expansion from southern refugia since the last glacial maximum.

METHODS

Six noncoding regions of chloroplast DNA were sequenced in 60 individuals. Populations were divided into two groups, north and south of the glacial boundary, in addition to isolated populations in Arkansas and Honduras. Variation in populations north of the boundary was compared with variation in populations to the south and in the isolated populations.

KEY RESULTS

Nucleotide diversity was an order of magnitude lower in populations north of the glacial boundary than in those to the south. The Arkansas and Honduras populations showed no variation. The predominant haplotype in northern populations was also found in a Virginia population.

CONCLUSIONS

Reduced variation north of the glacial boundary suggests a founder event associated with range expansion since the last glacial maximum. Colonization probably occurred from populations in refugia located several hundred kilometers south of the glacial boundary. The results provide insight into the effects of past and current climate change on patterns of geographic distribution and genetic variation in aquatic plants.

Divergence and dispersal of the nematode-trapping fungus Arthrobotrys oligospora from China

Nematode-trapping fungi are of significant agricultural, forestry and ecological importance. However, relatively little is known about the patterns of genetic variation for any nematode-trapping fungus through its broad geographic and ecological contexts. Here, we analysed DNA sequence variation among strains of the nematode-trapping fungus Arthrobotrys oligospora from China. Our analyses revealed that the Chinese A. oligospora is a species complex with at least three divergent lineages (cryptic species). In addition, there was significant geographic structuring with unambiguous evidence for localized recombination within two of the three lineages in nature. However, evidence for clonal reproduction was also found. We discuss the implications of our results to the conservation and biocontrol application of A. oligospora in agriculture and forestry.

The Use of Bioinformatics for Studying HIV Evolutionary and Epidemiological History in South America

The South American human immunodeficiency virus type 1 (HIV-1) epidemic is driven by several subtypes (B, C, and F1) and circulating and unique recombinant forms derived from those subtypes. Those variants are heterogeneously distributed around the continent in a country-specific manner. Despite some inconsistencies mainly derived from sampling biases and analytical constrains, most of studies carried out in the area agreed in pointing out specificities in the evolutionary dynamics of the circulating HIV-1 lineages. In this paper, we covered the theoretical basis, and the application of bioinformatics methods to reconstruct the HIV spatial-temporal dynamics, unveiling relevant information to understand the origin, geographical dissemination and the current molecular scenario of the HIV epidemic in the continent, particularly in the countries of Southern Cone.

Population structure and demographic inferences concerning the endangered onychophoran species Epiperipatus acacioi (Onychophora: Peripatidae)

Epiperipatus acacioi (Onychophora: Peripatidae) is an endemic species of the Atlantic rainforest in southeastern Brazil, with a restricted known distribution, found only in two nearby areas (Tripuí and Itacolomi). Mitochondrial gene COI sequences of 93 specimens collected across the known range of E. acacioi were used to assess the extant genetic diversity and patterns of genetic structure, as well as to infer the demographic history of this species. We found considerable variability within the populations, even though there has been recent environmental disturbance in these habitats. The samples from the two areas where this species is found showed significantly different COI sequences and constitute two distinct populations [exact test of sample differentiation (P = 0.0008) and pairwise F(ST) analyses (F(ST) = 0.214, P < 0.00001)]. However, there was little genetic differentiation among samples from different sampling sites within populations, suggesting that the potential for dispersal of E. acacioi greater than would have been expected, based on their cryptic behavior and reduced vagility. Mismatch analyses and neutrality tests revealed evidence of recent population expansion processes for both populations, possibly related to variations in the past distribution of this species.