Bayesian Coalescent Approach Research Articles

Risk Factors Associated with Hepatitis C Subtypes and the Evolutionary History of Subtype 1a in Mexico.

The Hepatitis C Virus (HCV), with its diverse genotypes and subtypes, has significantly impacted the health of millions of people worldwide. Analyzing the risk factors is essential to understanding the spread of the disease and developing appropriate prevention strategies. This study aimed to identify risk factors associated with HCV subtype transmission and calculate the emergence time of subtype 1a in Mexico. A cross-sectional study was conducted from January 2014 to December 2018, involving 260 HCV-infected adults. HCV infection was confirmed via Enzyme-Linked Immunosorbent Assay, and viral load was measured by real-time PCR. Genotyping/subtyping tools were the Line Probe Assay and Sanger sequencing of the non-structural region 5B (NS5B). The most frequent HCV subtype was 1a (58.5%), followed by subtypes 1b (19.2%), 3a (13.1%), 2b (5.4%), 2a/2c (2.7%), 2a (0.8%), and 4a (0.4%). Intravenous drug use and tattoos were significant risk factors for subtypes 1a and 3a, while hemodialysis and blood transfusion were linked with subtype 1b. For the evolutionary analysis, 73 high-quality DNA sequences of the HCV subtype 1a NS5B region were used, employing a Bayesian coalescent analysis approach. This analysis suggested that subtype 1a was introduced to Mexico in 1976, followed by a diversification event in the mid-1980s. An exponential increase in cases was observed from 1998 to 2006, stabilizing by 2014. In conclusion, this study found that HCV subtypes follow distinct transmission routes, emphasizing the need for targeted prevention strategies. Additionally, the findings provide valuable insights into the origin of HCV subtype 1a. By analyzing the history, risk factors, and dynamics of the HCV epidemic, we have identified these measures: limiting the harm of intravenous drug trafficking, enhancing medical training and infrastructure, and ensuring universal access to antiviral treatments. The successful implementation of these strategies could lead to an HCV-free future in Mexico.

Small, odd and old: The mysterious Tarsius pumilus is the most basal Sulawesi tarsier.

In this study, we present the first genetic evidence of the phylogenetic position of Tarsius pumilus, the mountain tarsier of Sulawesi, Indonesia. This mysterious primate is the only Eastern tarsier species that occurs exclusively in cloud forests above 1800 m.a.s.l. It exhibits striking morphological peculiarities—most prominently its extremely reduced body size, which led to the common name of ‘pygmy tarsier’. However, our results indicate that T. pumilus is not an aberrant form of a lowland tarsier, but in fact, the most basal of all Sulawesi tarsiers. Applying a Bayesian multi-locus coalescent approach, we dated the divergence between the T. pumilus lineage and the ancestor of all other extant Sulawesi tarsiers to 9.88 Mya. This is as deep as the split between the two other tarsier genera Carlito (Philippine tarsiers) and Cephalopachus (Western tarsiers), and predates further tarsier diversification on Sulawesi by around 7 Myr. The date coincides with the deepening of the marine environment between eastern and western Sulawesi, which likely led to allopatric speciation between T. pumilus or its predecessor in the west and the ancestor of all other Sulawesi tarsiers in the east. As the split preceded the emergence of permanent mountains in western Sulawesi, it is unlikely that the shift to montane habitat has driven the formation of the T. pumilus lineage.

Deep genetic structure of a ground-herb along contrasting environments of seasonally dry understories in Amazonia and Cerrado as revealed from targeted genomic sequencing

AbstractCentral South America primarily comprises seasonally dry tropical domains, where a sharp rainy season defines phenological and yearly growth patterns. Here we examine if the genetic structure of the understorey ground-herb Chamaecostus subsessilis species complex, from South American seasonally dry forests, is defined by intrinsic contrasting environments of the Cerrado and southern Amazonia. We analysed a dataset of c. 20 kb of sequence data from targeted capture high-throughput sequencing using PCR-generated probes. We used an admixture model and population pairwise FST comparisons to interpret genetic structure. We also used a Bayesian coalescent approach with a relaxed clock to estimate divergence times between populations and haplotypes. Genetic structure reflects separate populations from Amazonia and from the Cerrado, whereas a haplotype-based phylogenetic tree provides further spatiotemporal evidence of this strong phylogeographical structure. Spatial arrangements among genetic lineages of the C. subsessilis species complex reveal that intraspecific evolution in Central South America seasonally dry understories is closely linked to environmental gradients.

Base-substitution mutation rate across the nuclear genome of Alpheus snapping shrimp and the timing of isolation by the Isthmus of Panama

BackgroundThe formation of the Isthmus of Panama and final closure of the Central American Seaway (CAS) provides an independent calibration point for examining the rate of DNA substitutions. This vicariant event has been widely used to estimate the substitution rate across mitochondrial genomes and to date evolutionary events in other taxonomic groups. Nuclear sequence data is increasingly being used to complement mitochondrial datasets for phylogenetic and evolutionary investigations; these studies would benefit from information regarding the rate and pattern of DNA substitutions derived from the nuclear genome.ResultsTo estimate the genome-wide neutral mutation rate (µ), genotype-by-sequencing (GBS) datasets were generated for three transisthmian species pairs in Alpheus snapping shrimp. A range of bioinformatic filtering parameters were evaluated in order to minimize potential bias in mutation rate estimates that may result from SNP filtering. Using a Bayesian coalescent approach (G-PhoCS) applied to 44,960 GBS loci, we estimated µ to be 2.64E−9 substitutions/site/year, when calibrated with the closure of the CAS at 3 Ma. Post-divergence gene flow was detected in one species pair. Failure to account for this post-split migration inflates our substitution rate estimates, emphasizing the importance of demographic methods that can accommodate gene flow.ConclusionsResults from our study, both parameter estimates and bioinformatic explorations, have broad-ranging implications for phylogeographic studies in other non-model taxa using reduced representation datasets. Our best estimate of µ that accounts for coalescent and demographic processes is remarkably similar to experimentally derived mutation rates in model arthropod systems. These results contradicted recent suggestions that the closure of the Isthmus was completed much earlier (around 10 Ma), as mutation rates based on an early calibration resulted in uncharacteristically low genomic mutation rates. Also, stricter filtering parameters resulted in biased datasets that generated lower mutation rate estimates and influenced demographic parameters, serving as a cautionary tale for the adherence to conservative bioinformatic strategies when generating reduced-representation datasets at the species level. To our knowledge this is the first use of transisthmian species pairs to calibrate the rate of molecular evolution from GBS data.

Biogeographic origins of southern African Silene (Caryophyllaceae)

Silene (Caryophyllaceae) is distributed predominantly in the northern Hemisphere, where it is most diverse around the Mediterranean Basin. The genus is also well represented in North Africa, extending into tropical, sub-Saharan and southern Africa. Eight native species are recognized in southern Africa, taxonomically placed in two sections: Elisanthe and Silene s.l. Although the taxonomy of the southern African taxa has recently been revised, their phylogenetic relationships and biogeographic history remain unclear. This study aims to infer the phylogenetic position and geographic origins of the southern African taxa. We generated DNA sequences of nuclear and plastid loci from several individuals belonging to all eight species of Silene recognized from southern Africa, and combined our DNA sequences with existing data representing species from major clades (i.e. sections) based on the recently revised Silene infrageneric taxonomy. We used a Bayesian coalescent species tree continuous diffusion approach to co-estimate the species tree and the ancestral areas of representative members of the genus. Our results show that the perennial southern African members of section Elisanthe form a strongly-supported clade with the Eurasian annual S. noctiflora and the Central Asian perennial S. turkestanica. The rest of the perennial species form a strongly-supported clade together with the annual S. aethiopica, which is nested in a larger Mediterranean clade comprising mostly annual species classified in section Silene s.l. Estimates of ancestral areas indicate a late Pleistocene dispersal to southern Africa from central and East Africa for the sub-Saharan members of section Silene s.l. The Elisanthe clade is inferred to have colonized southern Africa through long-distance dispersal from Eurasia during the late Pleistocene. Our findings support the hypothesis of a relatively recent colonization into southern Africa resulting from two independent dispersal events during the Pleistocene.

Multi-character taxonomic review, systematics, and biogeography of the Black-capped/Tawny-bellied Screech Owl (emMegascops/em ematricapilla/em-emM. watsonii/em) complex (Aves: Strigidae).

Megascops is the most species-rich owl genus in the New World, with 21 species currently recognized. Phylogenetic relationships within this genus are notoriously difficult to establish due to the considerable plumage similarity among species and polymorphism within species. Previous studies have suggested that the widespread lowland Amazonian M. watsonii might include more than one species, and that the Atlantic Forest endemic M. atricapilla is closely related to the M. watsonii complex, but these relationships are as yet poorly understood. A recently published phylogeny of Megascops demonstrated that M. watsonii is paraphyletic with respect to M. atricapilla and that genetic divergences among some populations of M. watsonii are equal to or surpass the degree of differentiation between some M. watsonii and M. atricapilla. To shed light on the taxonomic status of these species and populations within them, we conducted a multi-character study based on molecular, morphological, and vocal characters. We sequenced three mitochondrial (cytb, CO1 and ND2) and three nuclear genes (BF5, CHD and MUSK) for 49 specimens, covering most of the geographic ranges of M. watsonii and M. atricapilla, and used these sequences to estimate phylogenies under alternative Bayesian, Maximum Likelihood, and multilocus coalescent species tree approaches. We studied 252 specimens and vocal parameters from 83 recordings belonging to 65 individuals, distributed throughout the ranges of M. watsonii and M. atricapilla. We used Discriminant Function Analysis (DFA) to analyze both morphometric and vocal data, and a pairwise diagnostic test to evaluate the significance of vocal differences between distinct genetic lineages. Phylogenetic analyses consistently recovered six statistically well-supported clades whose relationships are not entirely in agreement with currently recognized species limits in M. watsonii and M. atricapilla. Morphometric analyses did not detect significant differences among clades. High plumage variation among individuals within clades was usually associated with the presence of two or more color morphs. By contrast, vocal analyses detected significant differentiation among some clades but considerable overlap among others, with some lineages (particularly the most widespread one) exhibiting significant regional variation. The combined results allow for a redefinition of species limits in both M. watsonii and M. atricapilla, with the recognition of four additional species, two of which we describe here as new. We estimated most cladogenesis in the Megascops atricapilla-M. watsonii complex as having taken place during the Plio-Pleistocene, with the development of the modern Amazonian and São Francisco drainages and the expansion and retraction of forest biomes during interglacial and glacial periods as likely events accounting for this relatively recent burst of diversification.

Phylogenetic relationships, genetic diversity and biogeography of menhadens, genus Brevoortia (Clupeiformes, Clupeidae)

Brevoortia Gill 1861 is a genus of the Clupeidae (Teleostei) that includes six species of fishes commonly known as menhadens in eastern North America and “savelhas” or “saracas” in southeastern South America. Species of Brevoortia are important components of the marine food web of coastal ecosystems in the Atlantic and contribute significantly to fisheries. In this study, the first phylogenetic and biogeographic hypotheses including all species of Brevoortia are presented. A total of 113 specimens were analyzed using three molecular markers (two mitochondrial: COI and 16s; and one nuclear: RAG2). Maximum Parsimony, Maximum Likelihood and Bayesian Inference were employed to estimate phylogenetic relationships. A Bayesian multispecies coalescent approach was used to estimate a dated phylogeny, which supported biogeographic analyses of ancestral geographic ranges. Results corroborate previous hypotheses that the four North Atlantic species are grouped in two clades, one composed of B. tyrannus (Latrobe, 1802) and B. patronus Goode, 1878, and the second including B. smithi Hildebrand, 1941 and B. gunteri Hildebrand, 1948. The South Atlantic B. aurea (Spix and Agassiz, 1829) and B. pectinata (Jenyns, 1842) form a third clade, which is sister to the clade composed of B. smithi and B. gunteri. The monophyly and validity of the six nominal species of Brevoortia were not supported. Results also indicate that Brevoortia originated in the North Atlantic during the middle Miocene (about 15 Mya). A cooling event of the tropical Atlantic at around 10 Mya likely facilitated the range expansion of the genus to the South Atlantic, whereas a significant warming of the tropical Atlantic waters during the late Miocene at 6–7 Mya possibly promoted the isolation between the northern and southern counterparts of that ancestral lineage. The relevance of the Florida Peninsula in association with sea level fluctuations for the diversification within Brevoortia is also discussed.

Hepatitis C Virus Transmission Among Men Who Have Sex With Men in Amsterdam: External Introductions May Complicate Microelimination Efforts.

It is unclear whether unrestricted access and high uptake of direct-acting antivirals (DAAs) is sufficient to eliminate hepatitis C virus (HCV) in high-risk populations such as men who have sex with men (MSM). This study presents historic trends and current dynamics of HCV transmission among MSM in Amsterdam based on sequence data collected between 1994 and 2019. Hypervariable region 1 sequences of 232 primary HCV infections and 56 reinfections were obtained from 244 MSM in care in Amsterdam. Maximum-likelihood phylogenies were constructed for HCV genotypes separately, and time-scaled phylogenies were constructed using a Bayesian coalescent approach. Transmission clusters were determined by Phydelity and trends in the proportion of unclustered sequences over time were evaluated using logistic regression. Seventy-six percent (218/288) of sequences were part of 21 transmission clusters and 13 transmission pairs. Transmission cluster sizes ranged from 3 to 44 sequences. Most clusters were introduced between the late 1990s and early 2010s and no new clusters were introduced after 2012. The proportion of unclustered sequences of subtype 1a, the most prevalent subtype in this population, fluctuated between 0% and 20% in 2009-2012, after which an increase occurred from 0% in 2012 to 50% in 2018. The proportion of external introductions of HCV infections among MSM in Amsterdam has recently increased, coinciding with high DAA uptake. Frequent international transmission events will likely complicate local microelimination efforts. Therefore, international collaboration combined with international scale-up of prevention, testing, and treatment of HCV infections (including reinfections) is warranted, in particular for local microelimination efforts.

A biogeographic–ecological approach to disentangle reticulate evolution in the Triatoma phyllosoma species group (Heteroptera: Triatominae), vectors of Chagas disease

Introgression and incomplete lineage sorting (ILS) are two of the main sources of gene-tree incongruence; both can confound the assessment of phylogenetic relationships among closely related species. The Triatoma phyllosoma species group is a clade of partially co-distributed and cross-fertile Chagas disease vectors. Despite previous efforts, the phylogeny of this group remains unresolved, largely because of substantial gene-tree incongruence. Here, we sequentially address introgression and ILS to provide a robust phylogenetic hypothesis for the T. phyllosoma species group. To identify likely instances of introgression prior to molecular scrutiny, we assessed biogeographic data and information on fertility of inter-specific crosses. We first derived a few explicit hybridization hypotheses by considering the degree of spatial overlap within each species pair. Then, we assessed the plausibility of these hypotheses in the light of each species pair's cross-fertility. Using this contextual information, we evaluated mito-nuclear (cyt b, ITS-2) gene-tree incongruence and found evidence suggesting introgression within two species pairs. Finally, we modeled ILS using a Bayesian multispecies coalescent approach and either (a) a “complete” dataset with all the specimens in our sample, or (b) a “filtered” dataset without putatively introgressed specimens. The “filtered tree” had higher posterior-probability support, as well as more plausible topology and divergence times, than the “complete tree.” Detecting and filtering out introgression and modeling ILS allowed us to derive an improved phylogenetic hypothesis for the T. phyllosoma species group. Our results illustrate how biogeographic and ecological-reproductive contextual information can help clarify the systematics and evolution of recently diverged taxa prone to introgression and ILS.

Evolutionary units delimitation and continental multilocus phylogeny of the hyperdiverse catfish genus Hypostomus

With 149 currently recognized species, Hypostomus is one of the most species-rich catfish genera in the world, widely distributed over most of the Neotropical region. To clarify the evolutionary history of this genus, we reconstructed a comprehensive phylogeny of Hypostomus based on four nuclear and two mitochondrial markers. A total of 206 specimens collected from the main Neotropical rivers were included in the present study. Combining morphology and a Bayesian multispecies coalescent (MSC) approach, we recovered 85 previously recognized species plus 23 putative new species, organized into 118 ‘clusters’. We presented the Cluster Credibility (CC) index that provides numerical support for every hypothesis of cluster delimitation, facilitating delimitation decisions. We then examined the correspondence between the morphologically identified species and their inter-specific COI barcode pairwise divergence. The mean COI barcode divergence between morphological sisters species was 1.3 ± 1.2%, and only in 11% of the comparisons the divergence was ≥2%. This indicates that the COI barcode threshold of 2% classically used to delimit fish species would seriously underestimate the number of species in Hypostomus, advocating for a taxon-specific COI-based inter-specific divergence threshold to be used only when approximations of species richness are needed. The phylogeny of the 108 Hypostomus species, together with 35 additional outgroup species, confirms the monophyly of the genus. Four well-supported main lineages were retrieved, hereinafter called super-groups: Hypostomus cochliodon, H. hemiurus, H. auroguttatus, and H. plecostomus super-groups. We present a compilation of diagnostic characters for each super-group. Our phylogeny lays the foundation for future studies on biogeography and on macroevolution to better understand the successful radiation of this Neotropical fish genus.

Deciphering phylogenetic relationships and delimiting species boundaries using a Bayesian coalescent approach in protists: A case study of the ciliate genus Spirostomum (Ciliophora, Heterotrichea)

The ciliate genus Spirostomum comprises eight morphospecies, inhabiting diverse aquatic environments worldwide, where they can be used as water quality indicators. Although Spirostomum species are relatively easily identified using morphological methods, the previous nuclear rDNA-based phylogenies indicated several conflicts in morphospecies delineation. Moreover, the single locus phylogenies and previous analytical approaches could not unambiguously resolve phylogenetic relationships among Spirostomum morphospecies. Here, we attempt to investigate species boundaries and evolutionary history of Spirostomum taxa, using 166 new sequences from multiple populations employing one mitochondrial locus (CO1 gene) and two nuclear loci (rRNA operon and alpha-tubulin gene). In accordance with previous studies, relationships among the eight Spirostomum morphospecies were poorly supported statistically in individual gene trees. To overcome this problem, we utilised for the first time in ciliates the Bayesian coalescent approach, which accounts for ancestral polymorphisms, incomplete lineage sorting, and recombination. This strategy enabled us to robustly resolve deep relationships between Spirostomum species and to support the hypothesis that taxa with compact macronucleus and taxa with moniliform macronucleus each form a distinct lineage. Bayesian coalescent-based delimitation analyses strongly statistically supported the traditional morphospecies concept but also indicated that there are two S. minus-like cryptic species and S. teres is non-monophyletic. Spirostomum teres was very likely defined by a set of ancestral features of lineages that also gave rise to S. yagiui and S. dharwarensis. However, molecular data from type populations of the morphospecies S. minus and S. teres are required to unambiguously resolve the taxonomic problems.

The First Complete Genome Sequences of Hepatitis C Virus Subtype 2b from Latin America: Molecular Characterization and Phylogeographic Analysis.

The hepatitis C virus (HCV) has remarkable genetic diversity and exists as eight genotypes (1 to 8) with distinct geographic distributions. No complete genome sequence of HCV subtype 2b (HCV-2b) is available from Latin American countries, and the factors underlying its emergence and spread within the continent remain unknown. The present study was conducted to determine the first full-length genomic sequences of HCV-2b isolates from Latin America and reconstruct the spatial and temporal diversification of this subtype in Brazil. Nearly complete HCV-2b genomes isolated from two Brazilian patients were obtained by direct sequencing of long PCR fragments and analyzed together with reference sequences using the Bayesian coalescent and phylogeographic framework approaches. The two HCV-2b genomes were 9318 nucleotides (nt) in length (nt 37–9354). Interestingly, the long RT-PCR technique was able to detect co-circulation of viral variants that contained an in-frame deletion of 2022 nt encompassing E1, E2, and p7 proteins. Spatiotemporal reconstruction analyses suggest that HCV-2b had a single introduction in Brazil during the early 1980s, displaying an epidemic history characterized by a low and virtually constant population size until the present time. These results coincide with epidemiological data in Brazil and may explain the low national prevalence of this subtype.

Tracing the epidemic history of hepatitis C virus genotype 1b in Tunisia and in the world, using a Bayesian coalescent approach

Tracing the epidemic history of hepatitis C virus genotype 1b in Tunisia and in the world, using a Bayesian coalescent approach

Population genetics of Alaska Common Raven show dispersal and isolation in the world's largest songbird

The Common Raven (Corvus corax) is widespread at high latitudes and a noted disperser to remote islands. We hypothesized that, given their dispersal abilities, Alaskan populations would be genetically similar and maintain genetic diversity across thousands of kilometers. We sampled 134 ravens from 9 areas in Alaska, including 6 populations we considered as mainland and 3 as island (Kodiak, Adak, and Attu islands). Using 8 microsatellite loci, we found that in most Alaskan populations gene flow is sufficient to counter population divergence due to genetic drift. Our dataset suggests that conservatively there are just 2 raven populations in Alaska (K = 2), although the second Aleutian population, from Adak Island, was significantly differentiated, too. The Attu Island population was quite different from the other sampled populations, and it was the only location that neither received immigrants nor sent emigrants to other sampled locations. In addition, a comparison of gene flow models using a Bayesian coalescent approach most strongly supported an Attu Island isolation model. Prior work suggested that the uniqueness of the Attu Island population is due to isolation in a glacial refugium during the last glacial maximum. Postglacial dispersal from the Alaskan mainland has not erased this signal, suggesting that there could be limits to the species' dispersal abilities (i.e. remote colonization is rare). The isolated Attu Island population underwent a substantial decline ∼4 Kya, which roughly coincides with the arrival of the ultimate invasive species, humans. Traditional cultural values suggest that any surmised effect of humans would likely be indirect rather than direct.

Origin and spread of Thoroughbred racehorses inferred from complete mitochondrial genome sequences: Phylogenomic and Bayesian coalescent perspectives.

The Thoroughbred horse breed was developed primarily for racing, and has a significant contribution to the qualitative improvement of many other horse breeds. Despite the importance of Thoroughbred racehorses in historical, cultural, and economical viewpoints, there was no temporal and spatial dynamics of them using the mitogenome sequences. To explore this topic, the complete mitochondrial genome sequences of 14 Thoroughbreds and two Przewalski’s horses were determined. These sequences were analyzed together along with 151 previously published horse mitochondrial genomes from a range of breeds across the globe using a Bayesian coalescent approach as well as Bayesian inference and maximum likelihood methods. The racing horses were revealed to have multiple maternal origins and to be closely related to horses from one Asian, two Middle Eastern, and five European breeds. Thoroughbred horse breed was not directly related to the Przewalski’s horse which has been regarded as the closest taxon to the all domestic horses and the only true wild horse species left in the world. Our phylogenomic analyses also supported that there was no apparent correlation between geographic origin or breed and the evolution of global horses. The most recent common ancestor of the Thoroughbreds lived approximately 8,100–111,500 years ago, which was significantly younger than the most recent common ancestor of modern horses (0.7286 My). Bayesian skyline plot revealed that the population expansion of modern horses, including Thoroughbreds, occurred approximately 5,500–11,000 years ago, which coincide with the start of domestication. This is the first phylogenomic study on the Thoroughbred racehorse in association with its spatio-temporal dynamics. The database and genetic history information of Thoroughbred mitogenomes obtained from the present study provide useful information for future horse improvement projects, as well as for the study of horse genomics, conservation, and in association with its geographical distribution.

Disjunct distribution and distinct intraspecific diversification of Eothenomys melanogaster in South China

BackgroundSouth China encompasses complex and diverse landforms, giving rise to high biological diversity and endemism from the Hengduan Mountains to Taiwan Island. Many species are widely distributed across South China with similar disjunct distribution patterns. To explore the causes of these disjunct distribution patterns and their genetic consequences, we investigated the endemic species Père David’s Chinese Vole (Eothenomys melanogaster) by integrating geological and ecological factors. We analysed the genetic structure and divergence time of E. melanogaster based on fast-evolving mitochondrial and nuclear markers using Bayesian trees and coalescent species tree approaches. Historical scenarios of distribution range and demography were reconstructed based on spatial interpolations of genetic diversity and distance, extended Bayesian skyline plots, phylogeographic diffusion analysis, and ecological niche modelling (ENM) during different periods. We also assessed the relationships between geographical distance/ecological vicariance and genetic distance (isolation by distance, IBD; isolation by environment, IBE).ResultsThe genetic analysis revealed three deeply divergent clades—Southeast, Southwest and Central clades, centred on the Wuyi Mountains, the Yunnan-Guizhou Plateau (YGP) and the mountains around the Sichuan Basin, respectively—that have mostly developed since the Pleistocene. IBD played an important role in early divergence, and geological events (sedimentation of plains and linking of palaeo-rivers) and IBE further reinforced genetic differentiation. ENM shows the importance of suitable habitats and elevations.ConclusionsOur results suggest that the primary cause of the disjunct distribution in E. melanogaster is the high dependence on middle-high-altitude habitat in the current period. Mountains in the occurence range have served as “sky islands” for E. melanogaster and hindered gene flow. Pleistocene climatic cycles facilitated genetic admixture in cold periods and genetic diversification in warm periods for inland clades. During cold periods, these cycles led to multiple colonization events between the mainland and Taiwan and erased genetic differentiation.

Tracing the spatio-temporal dynamics of endangered fin whales (Balaenoptera physalus) within baleen whale (Mysticeti) lineages: a mitogenomic perspective

To explore the spatio-temporal dynamics of endangered fin whales (Balaenoptera physalus) within the baleen whale (Mysticeti) lineages, we analyzed 148 published mitochondrial genome sequences of baleen whales. We used a Bayesian coalescent approach as well as Bayesian inferences and maximum likelihood methods. The results showed that the fin whales had a single maternal origin, and that there is a significant correlation between geographic location and evolution of global fin whales. The most recent common female ancestor of this species lived approximately 9.88million years ago (Mya). Here, North Pacific fin whales first appeared about 7.48Mya, followed by a subsequent divergence in Southern Hemisphere approximately 6.63Mya and North Atlantic about 4.42 Mya. Relatively recently, approximately 1.76 and 1.42 Mya, there were two additional occurrences of North Pacific populations; one originated from the Southern Hemisphere and the other from an uncertain location. The evolutionary rate of this species was 1.002 × 10-3 substitutions/site/My. Our Bayesian skyline plot illustrates that the fin whale population has the rapid expansion event since ~ 2.5Mya, during the Quaternary glaciation stage. Additionally, this study indicates that the fin whale has a sister group relationship with humpback whale (Meganoptera novaeangliae) within the baleen whale lineages. Of the 16 genomic regions, NADH5 showed the most powerful signal for baleen whale phylogenetics. Interestingly, fin whales have 16 species-specific amino acid residues in eight mitochondrial genes: NADH2, COX2, COX3, ATPase6, ATPase8, NADH4, NADH5, and Cytb.

Species delimitation in the Caribbean Gesneria viridiflora complex (Gesneriaceae) reveals unsuspected endemism

AbstractAccurate taxonomy and species boundaries are of great importance in biodiversity hotspots with high species endemicity like the Caribbean. Indeed, inaccurate species delimitations can affect biodiversity estimates and influence the decisions taken on conservation issues. The genera Gesneria and Rhytidophyllum constitute the main representatives of the Caribbean Gesneriaceae and contain a few complexes with unclear taxonomic boundaries that are characterized by a confusing history of taxonomic changes and by the presence of several recognized varieties and subspecies. Gesneria viridiflora is a good example of such a complex. Four geographically isolated subspecies are recognized that have similar but variable vegetative and reproductive characters, and numerous taxonomic changes have been suggested in this group over the years. In this study, we used multivariate approaches to delimit distinct clusters of individuals using morphological (qualitative and quantitative characters) and molecular data from four nuclear markers. These groups are then tested using a Bayesian coalescent species delimitation approach and compared using multivariate analyses of bioclimatic variables obtained from occurrence data. The results suggest the presence of four distinct species in this complex according to the unified species concept: G. quisqueyana, G.sintenisii ,G . sylvicola andG .viridiflora . We also maintain G. viridiflora subsp. acrochordonanthe that does not fulfill our criteria for a species but that shows morphological variation associated to a specific geographical area. A distribution map, an identification key to the species, and a taxonomic treatment are provided. The new taxonomy proposed in this study shows an unsuspected species endemism in some regions of the Caribbean and underlines the importance of investigating species delimitation in diversified groups containing taxonomically complex taxa with poorly defined boundaries.

Probing the compartmentalization of HIV-1 in the central nervous system through its neutralization properties.

Compartmentalization of HIV-1 has been observed in the cerebrospinal fluid (CSF) of patients at different clinical stages. Considering the low permeability of the blood-brain barrier, we wondered if a reduced selective pressure by neutralizing antibodies (NAb) in the central nervous system (CNS) could favor the evolution of NAb-sensitive viruses in this compartment. Single genome amplification (SGA) was used to sequence full-length HIV-1 envelope variants (453 sequences) from paired CSF and blood plasma samples in 9 subjects infected by HIV variants of various clades and suffering from diverse neurologic disorders. Dynamics of viral evolution were evaluated with a bayesian coalescent approach for individuals with longitudinal samples. Pseudotyped viruses expressing envelope glycoproteins variants representative of the quasi-species present in each compartment were generated, and their sensitivity to autologous neutralization, broadly neutralizing antibodies (bNAbs) and entry inhibitors was assessed. Significant compartmentalization of HIV populations between blood and CSF were detected in 5 out of 9 subjects. Some of the previously described genetic determinants for compartmentalization in the CNS were observed regardless of the HIV-1 clade. There was no difference of sensitivity to autologous neutralization between blood- and CSF-variants, even for subjects with compartmentalization, suggesting that selective pressure by autologous NAb is not the main driver of HIV evolution in the CNS. However, we observed major differences of sensitivity to sCD4 or to at least one bNAb targeting either the N160-V1V2 site, the N332-V3 site or the CD4bs, between blood- and CSF-variants in all cases. In particular, HIV-1 variants present in the CSF were more resistant to bNAbs than their blood counterpart in some cases. Considering the possible migration from CSF to blood, the CNS could be a reservoir of bNAb resistant viruses, an observation that should be considered for immunotherapeutic approaches.

Evolutionary History of the Morocco lizard-Fingered Geckos of the Saurodactylus brosseti Complex

Studies of biodiversity in the Maghreb have revealed high genetic diversity and divergent genetic lineages among many taxa including squamates. Geographic barriers such as the Atlas Mountains are one of the key factors promoting genetic differentiation and the high levels of endemism. The lizard-fingered gecko Saurodactylus brosseti is endemic to Morocco. Its range includes both sides of the Atlas Mountains, and although high diversity was previously identified within the species, much of the range was unsampled. To understand the evolutionary and biogeographical history of this species, we used mitochondrial and nuclear DNA sequence data from 64 populations sampled across the entire species range. We employed phylogenetic methods based on gene trees and species trees as well as a time calibrated Bayesian genealogy and coalescent species delimitation approaches. We uncovered four highly divergent and allopatric mitochondrial lineages that did not share haplotypes at variable nuclear loci, suggesting the four groups have been evolving independently since the Miocene, according to our molecular dating estimates. These results coupled with the geographic pattern of genetic diversity suggest a possible role of the Atlas Mountains for the divergence observed between the four lineages of S. brosseti, while each lineage probably later underwent several episodes of fragmentation followed by (re-) expansion during Pleistocene climatic oscillations. Bayesian species delimitation analysis indicates that the four lineages may well be distinct species but we suggest that detailed morphological analyses are needed prior to taxonomic changes. The four lineages represent ancient independent evolutionary units, and deserve conservation management as distinct entities.