Comparative phylogeographic summary statistics for testing simultaneous vicariance

Testing for simultaneous vicariance across comparative phylogeographic data sets is a notoriously difficult problem hindered by mutational variance, the coalescent variance, and variability across pairs of sister taxa in parameters that affect genetic divergence. We simulate vicariance to characterize the behaviour of several commonly used summary statistics across a range of divergence times, and to characterize this behaviour in comparative phylogeographic datasets having multiple taxon-pairs. We found Tajima's D to be relatively uncorrelated with other summary statistics across divergence times, and using simple hypothesis testing of simultaneous vicariance given variable population sizes, we counter-intuitively found that the variance across taxon pairs in Nei and Li's net nucleotide divergence (pi(net)), a common measure of population divergence, is often inferior to using the variance in Tajima's D across taxon pairs as a test statistic to distinguish ancient simultaneous vicariance from variable vicariance histories. The opposite and more intuitive pattern is found for testing more recent simultaneous vicariance, and overall we found that depending on the timing of vicariance, one of these two test statistics can achieve high statistical power for rejecting simultaneous vicariance, given a reasonable number of intron loci (> 5 loci, 400 bp) and a range of conditions. These results suggest that components of these two composite summary statistics should be used in future simulation-based methods which can simultaneously use a pool of summary statistics to test comparative the phylogeographic hypotheses we consider here.

The stability of P in coral reef fishes

The constancy of phenotypic variation and covariation is an assumption that underlies most recent investigations of past selective regimes and attempts to predict future responses to selection. Few studies have tested this assumption of constancy despite good reasons to expect that the pattern of phenotypic variation and covariation may vary in space and time. We compared phenotypic variance-covariance matrices (P) estimated for Populations of six species of distantly related coral reef fishes sampled at two locations on Australia's Great Barrier Reef separated by more than 1000 km. The intraspecific similarity between these matrices was estimated using two methods: matrix correlation and common principal component analysis. Although there was no evidence of equality between pairs of P, both statistical approaches indicated a high degree of similarity in morphology between the two populations for each species. In general, the hierarchical decomposition of the variance-covariance structure of these populations indicated that all principal components of phenotypic variance-covariance were shared but that they differed in the degree of variation associated with each of these components. The consistency of this pattern is remarkable given the diversity of morphologies and life histories encompassed by these species. Although some phenotypic instability was indicated, these results were consistent with a generally conserved pattern of multivariate selection between populations.

Genetic and phenotypic evidence supports evolutionary divergence of the American flamingo (Phoenicopterus ruber) population in the Galapagos islands

The Galapagos archipelago is characterized by a high degree of endemism across many taxa, linked to the archpelago's oceanic origin and distance from other colonizing land masses. A population of ~ 500 American Flamingos (Phoenicopterus ruber) resides in Galapagos, which is thought to share an historical origin with the American Flamingo currently found in the Caribbean region. Genetic and phenotypic parameters in American Flamingos from Galapagos and from the Caribbean were investigated. Microsatellite and microchondrial DNA markers data showed that the American Flamingo population in Galapagos differs genetically from that in the Caribbean. American Flamingos in Galapagos form a clade which differs by a single common nucleotide substitution from American Flamingos in the Caribbean. The genetic differentiation is also evident from nuclear DNA in that microsatellite data reveal a number of private alleles for the American Flamingo in Galapagos. Analysis of skeletal measurements showed that American Flamingos in Galapagos are smaller than those in the Caribbean primarily due to shorter tarsus length, and differences in body shape sexual dimorphism. American Flamingo eggs from Galapagos have smaller linear dimensions and volumes than those from the Caribbean. The findings are consistent with reproductive isolation of American Flamingo population in Galapagos.

Late Pleistocene divergence and subsequent population expansion of two closely related fish species, Japanese anchovy (Engraulis japonicus) and Australian anchovy (Engraulis australis)

Climatic oscillations during the Pleistocene ice ages produced great changes in species' geographical distribution and abundance, which could be expected to have genetic consequences. Living in the temperate upwelling zones of the northwestern Pacific, Ja

Coalescent analysis of mtDNA indicates Pleistocene divergence among three species of howler monkey (Alouatta spp.) and population subdivision within the Atlantic Coastal Forest species, A. guariba

We have used coalescent analysis of mtDNA cytochrome b (cyt b) sequences to estimate times of divergence of three species of Alouatta-A. caraya, A. belzebul, and A. guariba-which are in close geographic proximity. A. caraya is inferred to have diverged from the A. guariba/A. belzebul clade approximately 3.83 million years ago (MYA), with the later pair diverging approximately 1.55 MYA. These dates are much more recent than previous dates based on molecular-clock methods. In addition, analyses of new sequences from the Atlantic Coastal Forest species A. guariba indicate the presence of two distinct haplogroups corresponding to northern and southern populations with both haplogroups occurring in sympatry within Sao Paulo state. The time of divergence of these two haplogroups is estimated to be 1.2 MYA and so follows quite closely after the divergence of A. guariba and A. belzebul. These more recent dates point to the importance of Pleistocene environmental events as important factors in the diversification of A. belzebul and A. guariba. We discuss the diversification of the three Alouatta species in the context of recent models of climatic change and with regard to recent molecular phylogeographic analyses of other animal groups distributed in Brazil.

Divergence Involving Global Regulatory Gene Mutations in an Escherichia coli Population Evolving under Phosphate Limitation

Many of the important changes in evolution are regulatory in nature. Sequenced bacterial genomes point to flexibility in regulatory circuits but we do not know how regulation is remodeled in evolving bacteria. Here, we study the regulatory changes that emerge in populations evolving under controlled conditions during experimental evolution of Escherichia coli in a phosphate-limited chemostat culture. Genomes were sequenced from five clones with different combinations of phenotypic properties that coexisted in a population after 37 days. Each of the distinct isolates contained a different mutation in 1 of 3 highly pleiotropic regulatory genes (hfq, spoT, or rpoS). The mutations resulted in dissimilar proteomic changes, consistent with the documented effects of hfq, spoT, and rpoS mutations. The different mutations do share a common benefit, however, in that the mutations each redirect cellular resources away from stress responses that are redundant in a constant selection environment. The hfq mutation lowers several individual stress responses as well the small RNA-dependent activation of rpoS translation and hence general stress resistance. The spoT mutation reduces ppGpp levels, decreasing the stringent response as well as rpoS expression. The mutations in and upstream of rpoS resulted in partial or complete loss of general stress resistance. Our observations suggest that the degeneracy at the core of bacterial stress regulation provides alternative solutions to a common evolutionary challenge. These results can explain phenotypic divergence in a constant environment and also how evolutionary jumps and adaptive radiations involve altered gene regulation.

Genomics of Divergence along a Continuum of Parapatric Population Differentiation

The patterns of genomic divergence during ecological speciation are shaped by a combination of evolutionary forces. Processes such as genetic drift, local reduction of gene flow around genes causing reproductive isolation, hitchhiking around selected variants, variation in recombination and mutation rates are all factors that can contribute to the heterogeneity of genomic divergence. On the basis of 60 fully sequenced three-spined stickleback genomes, we explore these different mechanisms explaining the heterogeneity of genomic divergence across five parapatric lake and river population pairs varying in their degree of genetic differentiation. We find that divergent regions of the genome are mostly specific for each population pair, while their size and abundance are not correlated with the extent of genome-wide population differentiation. In each pair-wise comparison, an analysis of allele frequency spectra reveals that 25–55% of the divergent regions are consistent with a local restriction of gene flow. Another large proportion of divergent regions (38–75%) appears to be mainly shaped by hitchhiking effects around positively selected variants. We provide empirical evidence that alternative mechanisms determining the evolution of genomic patterns of divergence are not mutually exclusive, but rather act in concert to shape the genome during population differentiation, a first necessary step towards ecological speciation.

Wavelet analysis and hierarchical clustering of ofm phenotypic population's variations originated from southern america and europe and attempts of a first relation to molecular divergence.

In this work we compare Grapholita molesta Busck (Lepidoptera: Tortricidae) populations originated from Brazil, Chile, Spain, Italy and Greece using power spectral density and phylogenetic analysis to detect any similarities between the population macro- and the molecular micro-level. Log-transformed population data were normalized and AR(p) models were developed to generate for each case population time series of equal lengths. The time-frequency/scale properties of the population data were further analyzed using wavelet analysis to detect any population dynamics frequency changes and cluster the populations. Based on the power spectral of each population time series and the hierarchical clustering schemes, populations originated from Southern America (Brazil and Chile) exhibit similar rhythmic properties and are both closer related with populations originated from Greece. Populations from Spain and especially Italy, have higher distance by terms of periodic changes on their population dynamics. Moreover, the members within the same cluster share similar spectral information, therefore they are supposed to participate in the same temporally regulated population process. On the contrary, the phylogenetic approach revealed a less structured pattern that bears indications of panmixia, as the two clusters contain individuals from both Europe and South America. This preliminary outcome will be further assessed by incorporating more individuals and likely employed a second molecular marker.

Time dependency of molecular rate estimates and systematic overestimation of recent divergence times

Studies of molecular evolutionary rates have yielded a wide range of rate estimates for various genes and taxa. Recent studies based on population-level and pedigree data have produced remarkably high estimates of mutation rate, which strongly contrast with substitution rates inferred in phylogenetic (species-level) studies. Using Bayesian analysis with a relaxed-clock model, we estimated rates for three groups of mitochondrial data: avian protein-coding genes, primate protein-coding genes, and primate d-loop sequences. In all three cases, we found a measurable transition between the high, short-term (<1–2 Myr) mutation rate and the low, long-term substitution rate. The relationship between the age of the calibration and the rate of change can be described by a vertically translated exponential decay curve, which may be used for correcting molecular date estimates. The phylogenetic substitution rates in mitochondria are approximately 0.5% per million years for avian protein-coding sequences and 1.5% per million years for primate protein-coding and d-loop sequences. Further analyses showed that purifying selection offers the most convincing explanation for the observed relationship between the estimated rate and the depth of the calibration. We rule out the possibility that it is a spurious result arising from sequence errors, and find it unlikely that the apparent decline in rates over time is caused by mutational saturation. Using a rate curve estimated from the d-loop data, several dates for last common ancestors were calculated: modern humans and Neandertals (354 ka; 222–705 ka), Neandertals (108 ka; 70–156 ka), and modern humans (76 ka; 47–110 ka). If the rate curve for a particular taxonomic group can be accurately estimated, it can be a useful tool for correcting divergence date estimates by taking the rate decay into account. Our results show that it is invalid to extrapolate molecular rates of change across different evolutionary timescales, which has important consequences for studies of populations, domestication, conservation genetics, and human evolution.

Pronounced genetic population structure in a potentially vagile fish species (Pristipomoides multidens, Teleostei; Perciformes; Lutjanidae) from the East Indies triangle

The East Indies triangle, bordered by the Phillipines, Malay Peninsula and New Guinea, has a high level of tropical marine species biodiversity. Pristipomoides multidens is a large, long-lived, fecund snapper species that is distributed throughout the East Indies and Indo-Pacific. Samples were analysed from central and eastern Indonesia and northern Australia to test for genetic discontinuities in population structure. Fish (n = 377) were collected from the Indonesian islands of Bali, Sumbawa, Flores, West Timor, Tanimbar and Tual along with 131 fish from two northern Australian locations (Arafura and Timor Seas) from a previous study. Genetic variation in the control region of the mitochondrial genome was assayed using restriction fragment length polymorphism and direct sequencing. Haplotype diversity was high (0.67-0.82), as was intraspecific sequence divergence (range 0-5.8%). FST between pairs of populations ranged from 0 to 0.2753. Genetic subdivision was apparent on a small spatial scale; FST was 0.16 over 191 km (Bali/Sumbawa) and 0.17 over 491 km (Bali/Flores). Constraints to dispersal that contribute to, and maintain, the observed degree of genetic subdivision are experienced presumably by all life history stages of this tropical marine finfish. The constraints may include (1) little or no movement of eggs or larvae, (2) little or no home range or migratory movement of adults and (3) loss of larval cohorts due to transport of larvae away from suitable habitat by prevailing currents

Genetic population structure of red snappers (Lutjanus malabaricus Bloch & Schneider, 1801 and Lutjanus erythropterus Bloch, 1790) in central and eastern Indonesia and northern Australia

The genetic population structure of red snapper Lutjanus malabaricus and Lutjanus erythropterus in eastern Indonesia and northern Australia was investigated by allozyme electrophoresis and sequence variation in the control region of mtDNA. Samples were collected from eight sites in Indonesia and four sites in northern Australia for both species. A total of 13 allozyme loci were scored. More variable loci were observed in L. malabaricus than in L. erythropterus. Sequence variation in the control region (left domain) of the mitochondrial genome was assessed by RFLP and direct sequencing. MtDNA haplotype diversity was high (L. erythropterus, 0.95 and L. malabaricus, 0.97), as was intraspecific sequence divergence, (L. erythropterus, 0.0-12.5% and L. malabaricus, 0.0-9.5%). The pattern of mtDNA haplotype frequencies grouped both species into two broad fisheries stocks with a genetic boundary either between Kupang and Sape (L. malabaricus) or between Kupang and Australian Timor Sea (L. erythropertus). The allozyme analyses revealed similar boundaries for L. erythropterus. Seven allozymes stocks compared to two mtDNA stocks of L. malabaricus including Ambon, which was not sampled with mtDNA, however, were reported. Possible reasons for differences in discrimination between the methods include: i) increased power of multiple allozyme loci over the single mtDNA locus, ii) insufficient gene sampling in the mtDNA control region and iii) relative evolutionary dynamics of nuclear (allozyme loci) and mitochondrial DNA in these taxa. Allozyme and haplotype data did not distinguish separate stocks among the four Australian locations nor the central Indonesian (Bali and Sape locations) for both L. malabaricus and L. erythropterus.

Population Expansion and Genetic Structure in Carcharhinus brevipinna in the Southern Indo-Pacific

Background:Quantifying genetic diversity and metapopulation structure provides insights into the evolutionary history of a species and helps develop appropriate management strategies. We provide the first assessment of genetic structure in spinner sharks (Carcharhinus brevipinna), a large cosmopolitan carcharhinid, sampled from eastern and northern Australia and South Africa. Methods and Findings:Sequencing of the mitochondrial DNA NADH dehydrogenase subunit 4 gene for 430 individuals revealed 37 haplotypes and moderately high haplotype diversity (h = 0.6770 ±0.025). While two metrics of genetic divergence (ΦST and FST) revealed somewhat different results, subdivision was detected between South Africa and all Australian locations (pairwise ΦST, range 0.02717–0.03508, p values ≤ 0.0013; pairwise FST South Africa vs New South Wales = 0.04056, p = 0.0008). Evidence for fine-scale genetic structuring was also detected along Australia’s east coast (pairwise ΦST = 0.01328, p < 0.015), and between south-eastern and northern locations (pairwise ΦST = 0.00669, p < 0.04).Conclusions: The Indian Ocean represents a robust barrier to contemporary gene flow in C. brevipinna between Australia and South Africa. Gene flow also appears restricted along a continuous continental margin in this species, with data tentatively suggesting the delineation of two management units within Australian waters. Further sampling, however, is required for a more robust evaluation of the latter finding. Evidence indicates that all sampled populations were shaped by a substantial demographic expansion event, with the resultant high genetic diversity being cause for optimism when considering conservation of this commercially-targeted species in the southern Indo-Pacific.

Phylogeography and Adaptive Divergence of Three-spined Stickleback Populations

Genetic studies on phylogeography and adaptive divergence in Northern Hemisphere fish species such as three-spined stickleback (Gasterosteus aculeatus) provide an excellent opportunity to investigate genetic mechanisms underlying population differentiation. According to the theory, the process of population differentiation results from a complex interplay between random and deterministic processes as well historical factors. The main scope in this thesis was to study how historical factors like the Pleistocene ice ages have shaped the patterns molecular diversity in three-spined stickleback populations in Europe and how this information could be utilized in the conservation genetic context. Furthermore, identifying footprints of natural selection at the DNA level might be used in identifying genes involved in evolutionary change. Overall, the results from phylogeographic studies indicate that the three-spined stickleback has colonized the Atlantic basin relatively recently but constitutes three major evolutionary lineages in Europe. In addition, the colonization of freshwater appears to result from multiple and independent invasions by the marine conspecifics. Molecular data together with morphology suggest that the most divergent freshwater populations are located in the Balkan Peninsula and these populations deserve a special conservation genetic status without warranting further taxonomical classification. In order to investigate the adaptive divergence in Fennoscandian three-spined stickleback populations several approaches were used. First, sequence variability in the Eda-gene, coding for the number of lateral plates, was concordant with the previously observed global pattern. Full plated allele is in high frequencies among marine populations whereas low plated allele dominates in the freshwater populations. Second, a microsatellite based genome scan identified both indications of balancing and directional selection in the three-spined stickleback genome, i.e. loci with unusually similar or unusually different allele frequencies over populations. The directionally selected loci were mainly associated with the adaptation to freshwater. A follow up study conducting a more detailed analysis in a chromosome region containing a putatively selected gene locus identified a fairly large genomic region affected by natural selection. However, this region contained several gene predictions, all of which might be the actual target of natural selection. All in all, the phylogeographic and adaptive divergence studies indicate that most of the genetic divergence has occurred in the freshwater populations whereas the marine populations have remained relatively uniform.