Microsatellite polymorphisms of Sichuan golden monkeys

Previous study using protein electrophoresis shows no polymorphism in 44 nuclear loci of Sichuan golden monkey (Rhinopithecus roxellana), which limits our understandings of its population genetic patterns in the nuclear genome. In order to obtain sufficient information, we scanned 14 microsatellite loci in a sample of 32 individuals from its three major habitats (Minshan, Qinling and Shennongiia). A considerable amount of polymorphisms were detected. The average heterozygosities in the local populations were all above 0.5. The differentiations among local populations were significant. There was evidence of geneflow among subpopullations, but geneflow between Qinling and Shennongjia local populations was the weakest. Minshan and Qinling populations might have gone through recent bottlenecks. The estimation of the ratio of the effective population sizes among local populations was close to that from census sizes. Comparisons to available mitochondria data suggested that R. roxellana's social structures played an important role in shaping its population genetic patterns. Our study showed that the polymorphism level of R. roxellana was no higher than other endangered species; therefore, measures should be taken to preserve genetic diversity of this species.

Live discrimination of Calanus glacialis and C. finmarchicus females: can we trust phenological differences?

Two key players in the Arctic and subarctic marine ecosystem are the calanoid copepods, Calanus finmarchicus and C. glacialis. Although morphologically very similar, these sibling species have different life cycles and roles in the Arctic pelagic marine ecosystem. Considering that the distribution of C. glacialis corresponds to Arctic water masses and C. finmarchicus to Atlantic water masses, the species are frequently used as climate indicators. Consequently, correct identification of the two species is essential if we want to understand climate-impacted changes on Calanus-dominated marine ecosystems such as the Arctic. Here, we present a novel morphological character (redness) to distinguish live females of C. glacialis and C. finmarchicus and compare it to morphological (prosome length) and genetic identification. The characters are tested on 300 live females of C. glacialis and C. finmarchicus from Disko Bay, western Greenland. Our analysis confirms that length cannot be used as a stand-alone criterion for separation. The results based on the new morphological character were verified genetically using a single mitochondrial marker (16S) and nuclear loci (six microsatellites and 12 InDels). The pigmentation criterion was also used on individuals (n = 89) from Young Sound fjord, northeast Greenland to determine whether the technique was viable in different geographical locations. Genetic markers based on mitochondrial and nuclear loci were corroborative in their identification of individuals and revealed no hybrids. Molecular identification confirmed that live females of the two species from Greenlandic waters, both East and West, can easily be separated by the red pigmentation of the antenna and somites of C. glacialis in contrast to the pale opaque antenna and somites of C. finmarchicus, confirming that the pigmentation criterion is valid for separation of the two species

Phylogeographical analysis of two cold-tolerant plants with disjunct Lusitanian distributions does not support in situ survival during the last glaciation

Aim: We used a combination of modelling and genetic approaches to investigate whether Pinguicula grandiflora and Saxifraga spathularis, two species that exhibit disjunct Lusitanian distributions, may have persisted through the Last Glacial Maximum (LGM, c. 21 ka) in separate northern and southern refugia.

Location: Northern and eastern Spain and south-western Ireland.

Methods: Palaeodistribution modelling using maxent was used to identify putative refugial areas for both species at the LGM, as well as to estimate their distributions during the Last Interglacial (LIG, c. 120 ka). Phylogeographical analysis of samples from across both species' ranges was carried out using one chloroplast and three nuclear loci for each species.

Results: The palaeodistribution models identified very limited suitable habitat for either species during the LIG, followed by expansion during the LGM. A single, large refugium across northern Spain and southern France was postulated for P. grandiflora. Two suitable regions were identified for S. spathularis: one in northern Spain, corresponding to the eastern part of the species' present-day distribution in Iberia, and the other on the continental shelf off the west coast of Brittany, south of the limit of the British–Irish ice sheet. Phylogeographical analyses indicated extremely reduced levels of genetic diversity in Irish populations of P. grandiflora relative to those in mainland Europe, but comparable levels of diversity between Irish and mainland European populations of S. spathularis, including the occurrence of private hapotypes in both regions.

Main conclusions: Modelling and phylogeographical analyses indicate that P. grandiflora persisted through the LGM in a southern refugium, and achieved its current Irish distribution via northward dispersal after the retreat of the ice sheets. Although the results for S. spathularis are more equivocal, a similar recolonization scenario also seems the most likely explanation for the species' current distribution.

Population genetic analyses reveal distinct geographical blooms of the jellyfish Rhizostoma octopus (Scyhpozoa)

Understanding the spatial integrity and connectivity of jellyfish blooms is important for ecologists and coastal stakeholders alike. Previous studies have shown that the distribution of jellyfish blooms can display a marked consistency in space and time, suggesting that such patterns cannot be attributed to passive processes alone. In the present study, we used a combination of microsatellite markers and mitochondrial cytochrome oxidase I sequences to investigate genetic structuring of the scyphozoan jellyfish Rhizostoma octopus in the Irish and Celtic Seas. The mitochondrial data indicated far higher levels of population differentiation than the microsatellites: ΦST[MT] = 0.300 vs. ΦST[NUC] = 0.013. Simulation studies indicated that the low levels of nuclear differentiation were not the result of limited power because of low levels of polymorphism. These findings, supported by palaeodistribution modelling and mismatch distribution analysis, are consistent with expansion of R. octopus from a single, limited refugium after the Last Glacial Maximum, followed by subsequent isolation, and that the discrepancy between the mitochondrial and nuclear markers is a result of the nuclear loci taking longer to reach mutation–drift equilibrium following the expansion as a result of their four-fold larger effective population size. The populations studied are probably not well connected via gene flow, and thus genetically as well as geographically distinct, although our findings also highlight the need to use a combination of organellar and nuclear markers to enable a more complete understanding of population demography and structure, particularly for species with large effective population sizes.

Statistical evaluation of alternative models of human evolution

An appropriate model of recent human evolution is not only important to understand our own history, but it is necessary to disentangle the effects of demography and selection on genome diversity. Although most genetic data support the view that our species originated recently in Africa, it is still unclear if it completely replaced former members of the Homo genus, or if some interbreeding occurred during its range expansion. Several scenarios of modern human evolution have been proposed on the basis of molecular and paleontological data, but their likelihood has never been statistically assessed. Using DNA data from 50 nuclear loci sequenced in African, Asian and Native American samples, we show here by extensive simulations that a simple African replacement model with exponential growth has a higher probability (78%) as compared with alternative multiregional evolution or assimilation scenarios. A Bayesian analysis of the data under this best supported model points to an origin of our species approximate to 141 thousand years ago (Kya), an exit out-of-Africa approximate to 51 Kya, and a recent colonization of the Americas approximate to 10.5 Kya. We also find that the African replacement model explains not only the shallow ancestry of mtDNA or Y-chromosomes but also the occurrence of deep lineages at some autosomal loci, which has been formerly interpreted as a sign of interbreeding with Homo erectus.

Cryptic speciation and recombination in the fungus Paracoccidioides brasiliensis as revealed by gene genealogies

Paracoccidioides brasiliensis is the etiologic agent of paracoccidioidomycosis, a disease confined to Latin America and of marked importance in the endemic areas due to its frequency and severity. This species is considered to be clonal according to mycological criteria and has been shown to vary in virulence. To characterize natural genetic variation and reproductive mode in this fungus, we analyzed P. brasiliensis phylogenetically in search of cryptic species and possible recombination using concordance and nondiscordance of gene genealogies with respect to phylogenies of eight regions in five nuclear loci. Our data indicate that this fungus consists of at least three distinct, previously unrecognized species: S1 (species 1 with 38 isolates), PS2 (phylogenetic species 2 with six isolates), and PS3 (phylogenetic species 3 with 21 isolates). Genealogies of four of the regions studied strongly supported the PS2 clade, composed of five Brazilian and one Venezuelan isolate. The second clade, PS3, composed solely of 21 Colombian isolates, was strongly supported by the alpha-tubulin genealogy. The remaining 38 individuals formed S1. Two of the three lineages of P. brasiliensis, S1 and PS2, are sympatric across their range, suggesting barriers to gene flow other than geographic isolation. Our study provides the first evidence for possible sexual reproduction in P. brasiliensis S1, but does not rule it out in the other two species.

Phylogenetic systematics of dart-poison frogs and their relatives (Amphibia: Athesphatanura: Dendrobatidae)

The known diversity of dart-poison frog species has grown from 70 in the 1960s to 247 at present, with no sign that the discovery of new species will wane in the foreseeable future. Although this growth in knowledge of the diversity of this group has been accompanied by detailed investigations of many aspects of the biology of dendrobatids, their phylogenetic relationships remain poorly understood. This study was designed to test hypotheses of dendrobatid diversification by combining new and prior genotypic and phenotypic evidence in a total evidence analysis. DNA sequences were sampled for five mitochondrial and six nuclear loci (approximately 6,100 base pairs [bp]; x=3,740 bp per terminal; total dataset composed of approximately 1.55 million bp), and 174 phenotypic characters were scored from adult and larval morphology, alkaloid profiles, and behavior. These data were combined with relevant published DNA sequences. Ingroup sampling targeted several previously unsampled species, including Aromobates nocturnus, which was hypothesized previously to be the sister of all other dendrobatids. Undescribed and problematic species were sampled from multiple localities when possible. The final dataset consisted of 414 terminals: 367 ingroup terminals of 156 species and 47 outgroup terminals of 46 species. Direct optimization parsimony analysis of the equally weighted evidence resulted in 25,872 optimal trees. Forty nodes collapse in the strict consensus, with all conflict restricted to conspecific terminals. Dendrobatids were recovered as monophyletic, and their sister group consisted of Crossodactylus, Hylodes, and Megaelosia, recognized herein as Hylodidae. Among outgroup taxa, Centrolenidae was found to be the sister group of all athesphatanurans except Hylidae, Leptodactyidae was polyphyletic, Thoropa was nested within Cycloramphidae, and Ceratophryinae was paraphyletic with respect to Telmatobiinae. Among dendrobatids, the monophyly and content of Mannophryne and Phyllobates were corroborated. Aromobates nocturnus and Colostethus saltuensis were found to be nested within Nephelobates, and Minyobates was paraphyletic and nested within Dendrobates. Colostethus was shown to be rampantly nonmonophyletic, with most species falling into two unrelated cis- and trans-Andean clades. A morphologically and behaviorally diverse clade of median lingual process-possessing species was discovered. In light of these findings and the growth in knowledge of the diversity of this large clade over the past 40 years, we propose a new, monophyletic taxonomy for dendrobatids, recognizing the inclusive clade as a superfamily (Dendrobatoidea) composed of two families (one of which is new), six subfamilies (three new), and 16 genera (four new). Although poisonous frogs did not form a monophyletic group, the three poisonous lineages are all confined to the revised family Dendrobatidae, in keeping with the traditional application of this name. We also propose changes to achieve a monophyletic higher-level taxonomy for the athesphatanuran outgroup taxa. Analysis of character evolution revealed multiple origins of phytotelm-breeding, parental provisioning of nutritive oocytes for larval consumption (larval oophagy), and endotrophy. Available evidence indicates that transport of tadpoles on the dorsum of parent nurse frogs-a dendrobatid synapomorphy-is carried out primitively by male nurse frogs, with three independent origins of female transport and five independent origins of biparental transport. Reproductive amplexus is optimally explained as having been lost in the most recent common ancestor of Dendrobatoidea, with cephalic amplexus arising independently three times. © American Museum of Natural History 2006.

Molecular phylogenetic relationships and phenotypic diversity in miniaturized toadlets, genus Brachycephalus (Amphibia: Anura: Brachycephalidae)

Toadlets of the genus Brachycephalus are endemic to the Atlantic rainforests of southeastern and southern Brazil. The 14 species currently described have snout-vent lengths less than 18. mm and are thought to have evolved through miniaturization: an evolutionary process leading to an extremely small adult body size. Here, we present the first comprehensive phylogenetic analysis for Brachycephalus, using a multilocus approach based on two nuclear (Rag-1 and Tyr) and three mitochondrial (Cyt b, 12S, and 16S rRNA) gene regions. Phylogenetic relationships were inferred using a partitioned Bayesian analysis of concatenated sequences and the hierarchical Bayesian method (BEST) that estimates species trees based on the multispecies coalescent model. Individual gene trees showed conflict and also varied in resolution. With the exception of the mitochondrial gene tree, no gene tree was completely resolved. The concatenated gene tree was completely resolved and is identical in topology and degree of statistical support to the individual mtDNA gene tree. On the other hand, the BEST species tree showed reduced significant node support relative to the concatenate tree and recovered a basal trichotomy, although some bipartitions were significantly supported at the tips of the species tree. Comparison of the log likelihoods for the concatenated and BEST trees suggests that the method implemented in BEST explains the multilocus data for Brachycephalus better than the Bayesian analysis of concatenated data. Landmark-based geometric morphometrics revealed marked variation in cranial shape between the species of Brachycephalus. In addition, a statistically significant association was demonstrated between variation in cranial shape and genetic distances estimated from the mtDNA and nuclear loci. Notably, B. ephippium and B. garbeana that are predicted to be sister-species in the individual and concatenated gene trees and the BEST species tree share an evolutionary novelty, the hyperossified dorsal plate. © 2011 Elsevier Inc.

The first molecular phylogeny of Chilodontidae (Teleostei: Ostariophysi: Characiformes) reveals cryptic biodiversity and taxonomic uncertainty

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

Metodi molecolari nello studio della variabilità genetica di Tuber borchii

Tuber borchii (Ascomycota, order Pezizales) is highly valued truffle sold in local markets in Italy. Despite its economic importance, knowledge on its distribution and population variation is scarce. The objective of this work was to investigate the evolutionary forces shaping the genetic structure of this fungus using coalescent and phylogenetic methods to reconstruct the evolutionary history of populations in Italy. To assess population structure, 61 specimens were collected from 11 different Provinces of Italy. Sampling was stratified across hosts and habitats to maximize coverage in native oak and pine stands and both mychorrizae and fruiting bodies were collected. Samples were identified considering anatomo-morphological characters. DNA was extracted and both multilocus (AFLP) and single-locus (18 loci from rDNA, nDNA, and mtDNA) approaches were used to look for polymorphisms. Screening AFLP profiles, both Jaccard and Dice coefficients of similarity were utilized to transform binary matrix into a distance matrix and then to desume Neighbour-Joining trees. Though these are only preliminary examinations, phylogenetic trees were totally concordant with those deriving from single locus analyses. Phylogenetic analyses of the nuclear loci were performed using maximum likelihood with PAUP and a combined phylogenetic inference, using Bayesian estimation with all nuclear gene regions, was carried out. To reconstruct the evolutionary history, we estimated recurrent migration, migration across the history of the sample, and estimated the mutation and approximate age of mutations in each tree using SNAP Workbench. The combined phylogenetic tree using Bayesian estimation suggests that there are two main haplotypes that are difficult to be differentiated on the basis of morphology, of ecological parameters and symbiontic tree. Between these two lineages, that occur in sympatry within T. borchii populations, there is no evidence of recurrent migration. However, migration over the history of the sample was asymmetrical suggesting that isolation was a result of interrupted gene flow followed by range expansion. Low levels of divergence between the haplotypes indicate that there are likely to be two cryptic species within the T. borchii population sampled. Our results suggest that isolation between populations of T. borchii could have led to reproductive isolation between two lineages. This isolation is likely due to sympatric speciation caused by a multiple colonization from different refugia or a recent isolation. In attempting to determinate whether these haplotypes represent separate species or a partition of the same species we applied Biological and Mechanistic species Concepts. Notwithstanding, further analyses are necessary to evaluate if selection favoured premating or post-mating isolation.

Statistical evaluation of alternative models of human evolution

An appropriate model of recent human evolution is not only important to understand our own history, but it is necessary to disentangle the effects of demography and selection on genome diversity. Although most genetic data support the view that our species originated recently in Africa, it is still unclear if it completely replaced former members of the Homo genus, or if some interbreeding occurred during its range expansion. Several scenarios of modern human evolution have been proposed on the basis of molecular and paleontological data, but their likelihood has never been statistically assessed. Using DNA data from 50 nuclear loci sequenced in African, Asian and Native American samples, we show here by extensive simulations that a simple African replacement model with exponential growth has a higher probability (78%) as compared with alternative multiregional evolution or assimilation scenarios. A Bayesian analysis of the data under this best supported model points to an origin of our species approximately 141 thousand years ago (Kya), an exit out-of-Africa approximately 51 Kya, and a recent colonization of the Americas approximately 10.5 Kya. We also find that the African replacement model explains not only the shallow ancestry of mtDNA or Y-chromosomes but also the occurrence of deep lineages at some autosomal loci, which has been formerly interpreted as a sign of interbreeding with Homo erectus.

Persistence of neutral polymorphisms in Lake Victoria cichlid fish

Phylogenetic trees for groups of closely related species often have different topologies, depending on the genes used. One explanation for the discordant topologies is the persistence of polymorphisms through the speciation phase, followed by differential fixation of alleles in the resulting species. The existence of transspecies polymorphisms has been documented for alleles maintained by balancing selection but not for neutral alleles. In the present study, transspecific persistence of neutral polymorphisms was tested in the endemic haplochromine species flock of Lake Victoria cichlid fish. Putative noncoding region polymorphisms were identified at four randomly selected nuclear loci and tested on a collection of 12 Lake Victoria species and their putative riverine ancestors. At all loci, the same polymorphism was found to be present in nearly all the tested species, both lacustrine and riverine. Different polymorphisms at these loci were found in cichlids of other East African lakes (Malawi and Tanganyika). The Lake Victoria polymorphisms must have therefore arisen after the flocks now inhabiting the three great lakes diverged from one another, but before the riverine ancestors of the Lake Victoria flock colonized the Lake. Calculations based on the mtDNA clock suggest that the polymorphisms have persisted for about 1.4 million years. To maintain neutral polymorphisms for such a long time, the population size must have remained large throughout the entire period.

When vicars meet: A narrow contact zone between morphologically cryptic phylogeographic lineages of the rainforest skink, Carlia rubrigularis

Phylogeographic analyses of the fauna of the Australian wet tropics rainforest have provided strong evidence for long-term isolation of populations among allopatric refugia, yet typically there is no corresponding divergence in morphology. This system provides an opportunity to examine the consequences of geographic isolation, independent of morphological divergence, and thus to assess the broader significance of historical subdivisions revealed through mitochondrial DNA phylogeography. We have located and characterized a zone of secondary contact between two long isolated (mtDNA divergence > 15%) lineages of the skink Carlia rubrigularis using one mitochondrial and eight nuclear (two intron, six microsatellite) markers. This revealed a remarkably narrow (width < 3 km) hybrid zone with substantial linkage disequilibrium and strong deficits of heterozygotes at two of three nuclear loci with diagnostic alleles. Cline centers were coincident across loci. Using a novel form of likelihood analysis, we were unable to distinguish between sigmoidal and stepped cline shapes except at one nuclear locus for which the latter was inferred. Given estimated dispersal rates of 90-133 m x gen(-1/2) and assuming equilibrium, the observed cline widths suggest effective selection against heterozygotes of at least 22-49% and possibly as high as 70%. These observations reveal substantial postmating isolation, although the absence of consistent deviations from Hardy-Weinberg equilibrium at diagnostic loci suggests that there is little accompanying premating isolation. The tight geographic correspondence between transitions in mtDNA and those for nuclear genes and corresponding evidence for selection against hybrids indicates that these morphologically cryptic phylogroups could be considered as incipient species. Nonetheless, we caution against the use of mtDNA phylogeography as a sole criterion for defining species boundaries.