Population structure between environmentally transmitted vibrios and bobtail squids using nested clade analysis

Squids from the genus Euprymna (Cephalopoda: Sepiolidae) and their symbiotic bacteria Vibrio fischeri form a mutualism in which vibrios inhabit a complex light organ within the squid host. A host-mediated daily expulsion event seeds surrounding seawater with symbiotically capable V. fischeri that environmentally colonize newly hatched axenic Euprymna juveniles. Competition experiments using native and non-native Vibrio have shown that this expulsion/re-colonization phenomenon has led to cospeciation in this system in the Pacific Ocean; however, the genetic architecture of these symbiotic populations has not been determined. Using genetic diversity and nested clade analyses we have examined the variation and history of three allopatric Euprymna squid species (E. scolopes of Hawaii, E. hyllebergi of Thailand, and E. tasmanica from Australia) and their respective Vibrio symbionts. Euprymna populations appear to be very genetically distinct from each other, exhibiting little or no migration over large geographical distances. In contrast, Vibrio symbiont populations contain more diverse haplotypes, suggesting both host presence and unidentified factors facilitating long-distance migration structure in Pacific Vibrio populations. Findings from this study highlight the importance of how interactions between symbiotic organisms can unexpectedly shape population structure in phylogeographical studies.     

Interspecific phylogenetic analysis enhances intraspecific phylogeographical inference: a case study in Pinus lambertiana

Pinus lambertiana (sugar pine) is an economically and ecologically important conifer with a 1600-km latitudinal range extending from Oregon, USA, to northern Baja California, Mexico. Like all North American white pines (subsect. Strobus), sugar pine is highly susceptible to white pine blister rust, a disease caused by the fungus Cronartium ribicola. We conducted a chloroplast DNA (cpDNA) survey of Pinus subsect. Strobus with comprehensive geographical sampling of P. lambertiana. Sequence analysis of 12 sugar pine individuals revealed strong geographical differentiation for two chloroplast haplotypes. A diagnostic restriction site survey of an additional 72 individuals demarcated a narrow 150-km contact zone in northeastern California. In the contact zone, maternal (megagametophtye) and paternal (embryo) haplotypes were identified in 31 single seeds, demonstrating bidirectional pollen flow extending beyond the range of maternal haplotypes. The frequencies of the Cr1 allele for white pine blister rust major gene resistance, previously determined for 41 seed zones, differ significantly among seed zones that are fixed for the alternate haplotypes, or contain a mixture of both haplotypes. Interspecific phylogenetic analysis reveals that the northern sugar pine haplotype belongs to a clade that includes Pinus albicaulis (whitebark pine) and all of the East Asian white pines. Furthermore, there is little cpDNA divergence between northern sugar pine and whitebark pine (dS = 0.00058). These results are consistent with a Pleistocene migration of whitebark pine into North America and subsequent chloroplast introgression from whitebark pine to sugar pine. This study demonstrates the importance of placing phylogeographical results in a broader phylogenetic context.     

Phylogeographical patterns in chloroplast DNA variation within the Acacia acuminata (Leguminosae: Mimosoideae) complex in Western Australia

Abstract The Acacia acuminata complex includes three taxa, A. acuminata ssp. acuminata, A. acuminata ssp. burkittii and A. oldfieldii, along with several informal variants of A. acuminata. It is widespread throughout southern Australia with the centre of diversity in south‐west Western Australia. Phylogeographical patterns in the complex were investigated using a nested clade analysis of cpDNA RFLPs from 25 populations in Western Australia. Except for A. oldfieldii that was clearly identified as a distinct entity, haplotypes were not restricted to sub‐specific taxa or variants within A. acuminata. There was significant association between phylogenetic position of many haplotypes and their geographical distribution. The fine‐scale phylogeographical patterns were complex but at deeper levels in the phylogeny there was evidence of divergence between two lineages. The pattern of shared haplotypes between lineages suggests retention of ancestral polymorphism as a result of incomplete lineage sorting. The divergence of these lineages is consistent with fragmentation caused by climatic instability during the Pleistocene.     

Genetic connectivity in the Florida reef system: comparative phylogeography of commensal invertebrates with contrasting reproductive strategies

Effective spatial management of coral reefs including design of marine protected areas requires an understanding of interpopulation genetic connectivity. We assessed gene flow along 355 km of the Florida reef system and between Florida and Belize in three commensal invertebrates occupying the same host sponge (Callyspongia vaginalis) but displaying contrasting reproductive dispersal strategies: the broadcast-spawning brittle star Ophiothrix lineata and two brooding amphipods Leucothoe kensleyi and Leucothoe ashleyae. Multiple analytical approaches to sequence variation in the mitochondrial COI gene demonstrated a high degree of overall connectivity for all three species along the Florida reef system. Ophiothrix lineata showed significant genetic structuring between Florida and Belize, and a pattern of isolation by distance but no significant genetic structuring along the Florida coastline. Bayesian estimates of migration detected a strong southerly dispersal bias for O. lineata along the Florida reef system, contrary to the general assumption of northerly gene flow in this region based on the direction of the Florida Current. Both amphipods, despite direct development, also showed high gene flow along the Florida reef system. Multiple inferences of long-distance dispersal from a nested clade analysis support the hypothesis that amphipod transport, possibly in detached sponge fragments, could generate the high levels of overall gene flow observed. However, this transport mechanism appears much less effective across deep water as connectivity between Florida and Belize (1072 km) is highly restricted.     

Nested clade and phylogeographic analyses of the chub, Leuciscus cephalus (Teleostei, cyprinidae), in Greece: implications for Balkan Peninsula biogeography

Phylogenetic relationships among Greek populations of the chub, Leuciscus cephalus, were investigated using 600 bp of the cytochrome b gene. The aim of this study was to test the assumption that the main difference in ichthyological composition between both sides of the Balkan Peninsula is directly linked to differences in the dispersion mechanisms used by fish in order to extend their distribution range. Phylogenetic and nested clade analyses clearly showed that populations in Greece are significantly differentiated. Greek populations were found to descend from three lineages in three geographical provinces: Western, Central, and Eastern Greece. The chub reached Western Greece at the beginning of the Pleistocene and Eastern Greece during the mid-Pleistocene. Chub dispersion occurred mainly by river confluence due to sea level lowering and river capture in Western Greece and sea dispersal with low-salinity conditions within the Aegean Sea in Eastern Greece. However, in Central Greece, the original mtDNA lineage has presumably been lost owing to a genetic introgression following a second invasion from the Danube during the final stage of the last glaciation. This study provides new elements for a better understanding of the composition of the contemporary ichthyofauna in Greece and highlights possible evolutionary mechanisms responsible for the high endemism rate in the Western Greek biogeographic province.     

Bayesian analysis of molecular variance in pyrosequences quantifies population genetic structure across the genome of Lycaeides butterflies

The distribution of genetic variation within and among populations is commonly used to infer their demographic and evolutionary histories. This endeavour has the potential to benefit substantially from high‐throughput next‐generation sequencing technologies through a rapid increase in the amount of data available and a corresponding increase in the precision of parameter estimation. Here we report the results of a phylogeographic study of the North American butterfly genus Lycaeides using 454 sequence data. This study serves the dual purpose of demonstrating novel molecular and analytical methods for population genetic analyses with 454 sequence data and expanding our knowledge of the phylogeographic history of Lycaeides. We obtained 341 045 sequence reads from 12 populations that we were able to assemble into 15 262 contigs (most of which were variable), representing one of the largest population genetic data sets for a non‐model organism to date. We examined patterns of genetic variation using a hierarchical Bayesian analysis of molecular variance model, which provides precise estimates of genome‐level φ_ST while appropriately modelling uncertainty in locus‐specific φ_ST. We found that approximately 36% of sequence variation was partitioned among populations, suggesting historical or current isolation among the sampled populations. Estimates of pairwise genome‐level φ_ST were largely consistent with a previous phylogeographic model for Lycaeides, suggesting fragmentation into two to three refugia during Pleistocene glacial cycles followed by post‐Pleistocene range expansion and secondary contact leading to introgressive hybridization. This study demonstrates the potential of using genome‐level data to better understand the phylogeographic history of populations.     

In silico model‐based inference: A contemporary approach for hypothesis testing in network biology

Inductive inference plays a central role in the study of biological systems where one aims to increase their understanding of the system by reasoning backwards from uncertain observations to identify causal relationships among components of the system. These causal relationships are postulated from prior knowledge as a hypothesis or simply a model. Experiments are designed to test the model. Inferential statistics are used to establish a level of confidence in how well our postulated model explains the acquired data. This iterative process, commonly referred to as the scientific method, either improves our confidence in a model or suggests that we revisit our prior knowledge to develop a new model. Advances in technology impact how we use prior knowledge and data to formulate models of biological networks and how we observe cellular behavior. However, the approach for model‐based inference has remained largely unchanged since Fisher, Neyman and Pearson developed the ideas in the early 1900s that gave rise to what is now known as classical statistical hypothesis (model) testing. Here, I will summarize conventional methods for model‐based inference and suggest a contemporary approach to aid in our quest to discover how cells dynamically interpret and transmit information for therapeutic aims that integrates ideas drawn from high performance computing, Bayesian statistics, and chemical kinetics. © 2014 American Institute of Chemical Engineers Biotechnol. Prog., 30:1247–1261, 2014     

A Second Look at Mitochondrial DNA Variability in European Anchovy (Engraulis encrasicolus): Assessing Models of Population Structure and the Black Sea Isolation Hypothesis

Genetic architectures of marine fishes are generally shallow because of the large potential for gene flow in the sea. European anchovy, however, are unusual among small pelagic fishes in showing large differences among sub-basins and in harbouring two mtDNA phylogroups (‘A’ & ‘B’), representing 1.1–1.85 million years of separation. Here the mtDNA RFLP dataset of Magoulas et al. [1996, Mol. Biol. Evol. 13: 178–190] is re-examined to assess population models accounting for this subdivided population structure and to evaluate the zoogeographical origins of the two major phylogroups. Haplotype and nucleotide diversities are highest in the Ionian Sea and lowest in the Aegean and Black seas. However, this gradient is absent when ‘A’ and ‘B’ haplotypes are examined separately. Neither the self-sustaining nor the basin population models adequately describe anchovy population behaviour. Tests for neutrality, mismatch and nested clade analyses are concordant in depicting recent expansions of both phylogroups. Unimodel mismatch distributions and haplotype coalescences dating to the last (Eemian) interglacial (‘B’) and the Weichselian pleniglacial period (‘A’) indicate separate colonizations of the Mediterranean Basin. Phylogroup ‘A’ is unlikely to have arisen through continuous long-term isolation in the Black Sea because of climate extremes from displaced subpolar weather systems during the ice ages. Ancestors of both groups appear to have colonized the Mediterranean from the Atlantic in the late Pleistocene. Hence, zoogeographic models of anchovy in the Mediterranean must also include the eastern (and possibly southern) Atlantic.     

Response to selection for litter size in Danish Landrace pigs: a Bayesian analysis

A replicated selection experiment aimed at increasing litter size (total number of pigs born per litter) in Danish Landrace pigs was conducted from 1984 to 1991. The experiment included two selection and two control lines. In each generation, 30 and 14 first litters were produced in selection and control lines, respectively, and dams produced two litters. Each replicate, consisting of one selection and one control line, was founded from 60 families chosen randomly from the population at large. Family selection was practiced, and the criterion was the predicted breeding value for litter size computed using a repeatability animal model, and taking into account all available information. The data consisted of 947 records from 523 dams (424 dams had two litters) representing five cycles of selection of increased litter size. Data were analyzed from a Bayesian perspective, based on marginal posterior distributions of genetic parameters of interest. Marginalization was achieved using Gibbs sampling, with a single chain length of 1 205 000. After discarding the first 5 000 iterations, a sample was drawn every ten iterations, so 120 000 samples in total were saved. Densities were estimated and plotted, and summary statistics were computed from the estimated densities. The posterior means (± standard error) of heritability and repeatability were 0.22 ± 0.06 and 0.32 ± 0.05, respectively. These point estimates of genetic parameters were within the range of literature values, although on the high side. The posterior mean (± standard error) of genetic response to selection, defined as the difference between the mean breeding values of the selected lines and that of the base population, was 1.37 ± 0.43 pigs after five cycles of selection. The regression (through the origin) of breeding values in the selected lines on generation was 0.25 ± 0.08 pigs. Several informative priors constructed from information obtained with field data in this population were used to examine their influence on inferences. The priors were influential because of the relatively small scale of the experiment. An analysis excluding data from one of the control lines gave smaller genetic variance and heritability, and a smaller response to selection. However, it appears that selection for litter size is effective, but that the true rate of response is probably smaller than data from this experiment suggest.     

Statistical convenience vs biological insight: consequences of data transformation for the analysis of fitness variation in heterogeneous environments

In plants, more favourable environmental conditions can lead to dramatic increases in both mean fitness and variance in fitness. This results in data that violate the equality-of-variance assumption of anova, a problem that most empiricists would address by log-transforming fitness values. Using heuristic data sets and simple simulations, we show that anova on log-transformed fitness consistently fails to match the outcome of selection in a heterogeneous environment or its sensitivity to environmental frequency. Only anova based on relative fitness within environments accurately predicts the sensitivity of genotype selection to the frequency of alternative environments. Parallel analyses of variance based on absolute fitness and relative fitness can bracket the expected success of alternative genotypes under hard and soft selection, respectively. For example, for Sinapis arvensis growing in full sun and partial shade treatments, families achieving high fitness in the best environment are favoured under hard selection, whereas soft selection favours different families that achieve consistently good performance across environments. Based on these findings, we recommend that log-transformation of fitness should no longer be standard practice in ecological genetics studies. Weighted anova is a preferable method for dealing with unequal variances, and investigators should also make greater use of techniques such as quantile regression or resampling to describe and evaluate fitness variation across heterogeneous environments.