The divergence of two independent lineages of an endemic Chinese gecko,Gekko swinhonis,launched by the Qinling orogenic belt

The genetic structure and demographic history of an endemic Chinese gecko, Gekko swinhonis, were investigated by analysing the mitochondrial cytochrome b gene and 10 microsatellite loci for samples collected from 27 localities. Mitochondrial DNA data provided a detailed distribution of two highly divergent evolutionary lineages, between which the average pairwise distance achieved was 0.14. The geographic division of the two lineages coincided with a plate boundary consisting of the Qinling and Taihang Mts, suggesting a historical vicariant pattern. The orogeny of the Qinling Mts, a dispersal and major climatic barrier of the region, may have launched the independent lineage divergence. Both lineages have experienced recent expansion, and the current sympatric localities comprised the region of contact between the lineages. Individual‐based phylogenetic analyses of nucDNA and Bayesian‐clustering approaches revealed a deep genetic structure analogous to mtDNA. Incongruence between nucDNA and mtDNA at the individual level at localities outside of the contact region can be explained by the different inheritance patterns and male‐biased dispersal in this species. High genetic divergence, long‐term isolation and ecological adaptation, as well as the morphological differences, suggest the presence of a cryptic species.     

Distribution patterns of land snails in Ugandan rain forests support the existence of Pleistocene forest refugia

Aim We investigated whether the biogeographical patterns expected if the East African fauna was affected by cycles of contraction to refugia and expansion of ranges, as has been previously hypothesized, can be found in the land snail fauna of rain forests in Uganda. Location The Albertine Rift region and the Lake Victoria forest belt in Uganda. Methods Snails and slugs were sampled in 60 plots in 13 rain forests, and small species were extracted from 5‐L leaf litter samples. Relative species richness was calculated by rarefaction. The influence of putative determinants of species richness was examined by bivariate correlation and multiple regression. Clustering and nestedness were tested by Monte Carlo simulations with a null model that considers the range size distribution, the species richness distribution of the forests, and the spatial autocorrelation of the occurrences of each species. Biotic elements were determined by model‐based Gaussian clustering. Results A total of 169 land snail species were recorded from 13 Ugandan rain forests. Relative species richness increases with rainfall and altitude, and decreases with evaporation and distance from the putative East Congolian refugia. Mean annual rainfall and distance from the putative East Congolian refugia were included in the best multiple regression model. The distribution areas of the Ugandan land snails are significantly clustered. Two montane, two lowland and a northern biotic element were found. The mean range extension increases with increasing distance from the putative East Congolian refugia. Moreover, the ranges of the Ugandan land snails are significantly nested. The centre of the sets of nested subsets is in the Virunga Mountains, close to the putative East Congolian refugia. Main conclusions The decrease of diversity with increasing distance from the putative East Congolian refugia, the clustering and nestedness of ranges, and the range size increase with increasing distance from the refugia indicate that the East African land snail fauna was affected by cycles of contraction to refugia and expansion of ranges. The significant clustering and nestedness cannot be explained by current environmental conditions. Given the environmental history, it can be supposed that the lowland elements expanded post‐glacially, whereas the ranges of the montane species are probably currently contracting.     

Global phylogeography and seascape genetics of the lemon sharks (genus Negaprion)

Seascapes are complex environments, and populations are often isolated by factors other than distance. Here we investigate the role of coastal habitat preference and philopatry in shaping the distribution and population structure of lemon sharks. The genus Negaprion comprises the amphiatlantic lemon shark (N. brevirostris), with a relict population in the eastern Pacific, and its Indo‐West Pacific sister species, the sicklefin lemon shark (N. acutidens). Analyzing 138 individuals throughout the range of N. brevirostris (N = 80) and N. acutidens (N = 58) at microsatellite loci (nine and six loci, respectively) and the mitochondrial control region, we find evidence of allopatric speciation corresponding to the Tethys Sea closure (10–14 million years ago) and isolation of the eastern Pacific N. brevirostris population via the emergence of the Isthmus of Panama (~3.5 million years ago). There is significant isolation by oceanic distance (R² = 0.89, P = 0.005), defined as the maximum distance travelled at depths greater than 200 m. We find no evidence for contemporary transatlantic gene flow (m, M = 0.00) across an oceanic distance of ~2400 km. Negaprion acutidens populations in Australia and French Polynesia, separated by oceanic distances of at least 750 km, are moderately differentiated (F_ST = 0.070–0.087, P≤ 0.001; Φ_ST = 0.00, P = 0.99), with South Pacific archipelagos probably serving as stepping stones for rare dispersal events. Migration between coastally linked N. brevirostris populations is indicated by nuclear (m = 0.31) but not mitochondrial (m < 0.001) analyses, possibly indicating female natal site fidelity. However, philopatry is equivocal in N. acutidens, which has the lowest control region diversity (h = 0.28) of any shark yet studied. Restricted oceanic dispersal and high coastal connectivity stress the importance of both local and international conservation efforts for these threatened sharks.     

Historical biogeography and phenotype‐phylogeny of Chroodiscus (lichenized Ascomycota: Ostropales: Graphidaceae)

Aim To analyse the historical biogeography of the lichen genus Chroodiscus using a phenotype‐based phylogeny in the context of continental drift and evolution of tropical rain forest vegetation. Location All tropical regions (Central and South America, Africa, India, Southeast Asia, north‐east Australia). Methods We performed a phenotype‐based phylogenetic analysis and ancestral character state reconstruction of 14 species of the lichen genus Chroodiscus, using paup * and mesquite ; dispersal–vicariance analysis (DIVA) and dispersal–extinction–cladogenesis (DEC) modelling to trace the geographical origin of individual clades; and ordination and clustering by means of pc‐ord , based on a novel similarity index, to visualize the biogeographical relationships of floristic regions in which Chroodiscus occurs. Results The 14 species of Chroodiscus show distinctive distribution patterns, with one pantropical and one amphi‐Pacific taxon and 12 species each restricted to a single continent. The genus comprises four clades. DIVA and DEC modelling suggest a South American origin of Chroodiscus in the mid to late Cretaceous (120–100 Ma), with subsequent expansion through a South American–African–Indian–Southeast Asian–Australian dispersal route and late diversification of the argillaceus clade in Southeast Asia. Based on the abundance of extant taxa, the probability of speciation events in Chroodiscus is shown to be extremely low. Slow dispersal of foliicolous rain forest understorey lichens is consistent with estimated phylogenetic ages of individual species and with average lengths of biological species intervals in fungi (10–20 Myr). Main conclusions The present‐day distribution of Chroodiscus can be explained by vicariance and mid‐distance dispersal through the interconnection or proximity of continental shelves, without the need for recent, trans‐oceanic long‐distance dispersal. Phylogenetic reconstruction and age estimation for Chroodiscus are consistent with the ‘biotic ferry’ hypothesis: a South American origin and subsequent eastward expansion through Africa towards Southeast Asia and north‐eastern Australia via the Indian subcontinent. The present‐day pantropical distributions of many clades and species of foliicolous lichens might thus be explained by eastward expansion through continental drift, along with the evolution of modern rain forests starting 120 Ma, rather than by the existence of a hypothetical continuous area of pre‐modern rain forest spanning South America, Africa and Southeast Asia during the mid and late Cretaceous.     

Phylogeographic patterns of trans-Amazonian vicariants and Amazonian biogeography: the Neotropical rattlesnake (Crotalus durissus complex) as an example

Aim To investigate the phylogeography and execute a historical‐demographic analysis of the Neotropical rattlesnake, Crotalus durissus, thereby testing the hypothesis of a Pleistocene central Amazon corridor of dry forest or savanna that partitioned the Amazonian rain forest into western and eastern portions. Location South America. Methods Using sequences of three mitochondrial genes, we estimated the phylogeography, gene and nucleotide diversity across the South American range of C. durissus. Tree topology tests were used to test alternative biogeographical hypotheses, and tests of population genetic structure and statistical parsimony networks and nested clade phylogeographic analysis (NCPA) were used to infer connectivity and historical population processes on both sides of the Amazon basin. Results Tree topology tests rejected the hypothesis of a coastal dispersal in favour of a central corridor scenario. Gene diversity was similar on both sides of the Amazon basin. Nucleotide diversity indicated that the populations from north of the Amazon basin represented ancestral populations. Analysis of molecular variance (amova ) showed that intra‐population molecular variation was greater than between regions. Historical‐demographic statistics showed significant population expansion south of the Amazon, and little differentiation in the north, indicating moderate past gene flow between north and south of the Amazon. The parsimony network connected clades from the Roraima and Guyana populations with Mato Grosso, suggesting an Amazonian central corridor, and NCPA supported allopatric fragmentation between north and south of the Amazon. Main conclusions The distribution of C. durissus on both sides of the Amazon basin is evidence of changes in the distribution of rain forest vegetation during the Pleistocene. Our results suggest a formerly continuous distribution of this rattlesnake along a central Amazonian corridor during the middle Pleistocene. Allopatric fragmentation inferred from NCPA is consistent with vicariance resulting from a subsequent closure of this habitat corridor. This study emphasizes the potential of trans‐Amazonian open formation species to inform the debate on the past distribution of rain forests in the Amazon Basin.     

Molecular dating of the 'Gondwanan' plant family Proteaceae is only partially congruent with the timing of the break-up of Gondwana

Aim  The flowering plant family Proteaceae is putatively of Gondwanan age, with modern and fossil lineages found on all southern continents. Here we test whether the present distribution of Proteaceae can be explained by vicariance caused by the break-up of Gondwana.
Location  Africa, especially southern Africa, Australia, New Zealand, South America, New Caledonia, New Guinea, Southeast Asia, Sulawesi, Tasmania.
Methods  We obtained chloroplast DNA sequence data from the rbc L gene, the rbc L- atp B spacer, and the atp B gene from leaf samples of forty-five genera collected from the field and from living collections. We analysed these data using Bayesian phylogenetic and molecular dating methods, with five carefully selected fossil calibration points to obtain age estimates for the nodes within the family.
Results  Four of eight trans-continental disjunctions of sister groups within our sample of the Proteaceae post-date the break-up of Gondwana. These involve independent lineages, two with an Africa-Australia disjunction, one with an Africa–South America disjunction, and one with a New Zealand–Australasia disjunction. The date of the radiation of the bird-pollinated Embothriinae corresponds approximately to the hypothesized date of origin of nectar-feeding birds in Australia.
Main conclusions  The findings suggest that disjunct distributions in Proteaceae result from both Gondwanan vicariance and transoceanic dispersal. Our results imply that ancestors of some taxa dispersed across oceans rather than rafting with Gondwanan fragments as previously thought. This finding agrees with other studies of Gondwanan plants in dating the divergence of Australian, New Zealand and New Caledonian taxa in the Eocene, consistent with the existence of a shared, ancestral Eocene flora but contrary to a vicariance scenario based on accepted geological knowledge.     

Molecular and fossil evidence place the origin of cichlid fishes long after Gondwanan rifting

Cichlid fishes are a key model system in the study of adaptive radiation, speciation and evolutionary developmental biology. More than 1600 cichlid species inhabit freshwater and marginal marine environments across several southern landmasses. This distributional pattern, combined with parallels between cichlid phylogeny and sequences of Mesozoic continental rifting, has led to the widely accepted hypothesis that cichlids are an ancient group whose major biogeographic patterns arose from Gondwanan vicariance. Although the Early Cretaceous (ca 135 Ma) divergence of living cichlids demanded by the vicariance model now represents a key calibration for teleost molecular clocks, this putative split pre-dates the oldest cichlid fossils by nearly 90 Myr. Here, we provide independent palaeontological and relaxed-molecular-clock estimates for the time of cichlid origin that collectively reject the antiquity of the group required by the Gondwanan vicariance scenario. The distribution of cichlid fossil horizons, the age of stratigraphically consistent outgroup lineages to cichlids and relaxed-clock analysis of a DNA sequence dataset consisting of 10 nuclear genes all deliver overlapping estimates for crown cichlid origin centred on the Palaeocene (ca 65–57 Ma), substantially post-dating the tectonic fragmentation of Gondwana. Our results provide a revised macroevolutionary time scale for cichlids, imply a role for dispersal in generating the observed geographical distribution of this important model clade and add to a growing debate that questions the dominance of the vicariance paradigm of historical biogeography.     

Investigating population stratification and admixture using eigenanalysis of dense genotypes

Principal components analysis of genetic data is used to avoid inflation in type I error rates in association testing due to population stratification by covariate adjustment using the top eigenvectors and to estimate cluster or group membership independent of self-reported or ethnic identities. Eigendecomposition transforms correlated variables into an equal number of uncorrelated variables. Numerous stopping rules have been developed to identify which principal components should be retained. Recent developments in random matrix theory have led to a formal hypothesis test of the top eigenvalue, providing another way to achieve dimension reduction. In this study, I compare Velicer's minimum average partial test to a test on the basis of Tracy-Widom distribution as implemented in EIGENSOFT, the most widely used implementation of principal components analysis in genome-wide association analysis. By computer simulation of vicariance on the basis of coalescent theory, EIGENSOFT systematically overestimates the number of significant principal components. Furthermore, this overestimation is larger for samples of admixed individuals than for samples of unadmixed individuals. Overestimating the number of significant principal components can potentially lead to a loss of power in association testing by adjusting for unnecessary covariates and may lead to incorrect inferences about group differentiation. Velicer's minimum average partial test is shown to have both smaller bias and smaller variance, often with a mean squared error of 0, in estimating the number of principal components to retain. Velicer's minimum average partial test is implemented in R code and is suitable for genome-wide genotype data with or without population labels.