Model averaging and Bayes factor calculation of relaxed molecular clocks in Bayesian phylogenetics

Modeling Substitution Rate Evolution across Lineages and Relaxing the Molecular Clock

Relaxing the molecular clock using models of how substitution rates change across lineages has become essential for addressing evolutionary problems. The diversity of rate evolution models and their implementations are substantial, and studies have demonstrated their impact on divergence time estimates can be as significant as that of calibration information. In this review, we trace the development of rate evolution models from the proposal of the molecular clock concept to the development of sophisticated Bayesian and non-Bayesian methods that handle rate variation in phylogenies. We discuss the various approaches to modeling rate evolution, provide a comprehensive list of available software, and examine the challenges and advancements of the prevalent Bayesian framework, contrasting them to faster non-Bayesian methods. Lastly, we offer insights into potential advancements in the field in the era of big data.

Immune imprinting in early life shapes cross-reactivity to influenza B virus haemagglutinin

Influenza exposures early in life are believed to shape future susceptibility to influenza infections by imprinting immunological biases that affect cross-reactivity to future influenza viruses. However, direct serological evidence linked to susceptibility is limited. Here we analysed haemagglutination-inhibition titres in 1,451 cross-sectional samples collected between 1992 and 2020, from individuals born between 1917 and 2008, against influenza B virus (IBV) isolates from 1940 to 2021. We included testing of 'future' isolates that circulated after sample collection. We show that immunological biases are conferred by early life IBV infection and result in lineage-specific cross-reactivity of a birth cohort towards future IBV isolates. This translates into differential estimates of susceptibility between birth cohorts towards the B/Yamagata and B/Victoria lineages, predicting lineage-specific birth-cohort distributions of observed medically attended IBV infections. Our data suggest that immunological measurements of imprinting could be important in modelling and predicting virus epidemiology.

Climate change will lead to range shifts and genetic diversity losses of dung beetles in the Gobi Desert and Mongolian Steppe

Desertification is known to be a major threat to biodiversity, yet our understanding of the consequent decline in biodiversity remains insufficient. Here, we predicted climate change-induced range shifts and genetic diversity losses in three model dung beetles: Colobopterus erraticus, Cheironitis eumenes, and Gymnopleurus mopsus, distributed across the Gobi Desert and Mongolian Steppe, areas known for desertification. Phylogeographic analyses of mitochondrial COI sequences and species distribution modeling, based on extensive field investigations spanning 14 years, were performed. Species confined to a single biome were predicted to contract and shift their distribution in response to climate change, whereas widespread species was predicted to expand even if affected by range shifts. We indicated that all species are expected to experience significant haplotype losses, yet the presence of high singleton frequencies and low genetic divergence across geographic configurations and lineages mitigate loss of genetic diversity. Notably, Cheironitis eumenes, a desert species with low genetic diversity, appears to be the most vulnerable to climate change due to the extensive degradation in the Gobi Desert. This is the first study to predict the response of insects to desertification in the Gobi Desert. Our findings highlight that dung beetles in the Gobi Desert and Mongolian Steppe might experience high rates of occupancy turnover and genetic loss, which could reshuffle the species composition.

Species boundaries of the whiptail catfish Rineloricaria (Siluriformes: Loricariidae) from the Paraíba do Sul River drainage, southeastern Brazil, with species redescriptions and description of a new species

Species of the catfish genus Rineloricaria are common in the Paraíba do Sul River basin, in southeastern Brazil; here we present a revision of the taxonomic diversity and geographic distribution of the species of the genus inhabiting the basin, based on novel morphologic and molecular data. Five species delimitation methods based on cytochrome C oxidase subunit 1 nucleotide sequences yielded comparable molecular operational taxonomic units. The automatic barcode gap discovery, assemble species by automatic partitioning, barcode index number, and Bayesian implementation of the Poisson tree process methods supported the recognition of five evolutionary lineages. These taxonomic units were assigned to the previously described Rineloricaria nigricauda, Rineloricaria steindachneri, Rineloricaria zawadzkii, and Rineloricaria nudipectoris, and an additional undescribed species. R. zawadzkii was further divided into two intraspecific geographically structured lineages using the generalized mixed Yule coalescent delimitation method. A maximum likelihood phylogenetic analysis revealed that the five lineages from the Paraíba do Sul have closer relationships to different species from southern and southeastern Brazil (Ribeira de Iguape, Lagoa dos Patos, Uruguay, Paraguay, and Parana river basins) than to each other. Based on the analysis of lectotypes, recently collected material, and specimens from ichthyological collections, the poorly described R. nigricauda and R. steindachneri are redescribed following current descriptive standards. The undescribed species from the middle and upper Paraíba do Sul River basin is formally described. The description of a new species, along with the description of species boundaries in R. nigricauda and R. steindachneri, contributes to the knowledge of the ichthyofauna of the Paraíba do Sul River basin and adjacent coastal drainages of southeastern Brazil. An identification key for the species of Rineloricaria occurring in the Paraíba do Sul River basin is provided.

Molecular and morphometric data provide evidence of intraspecific variation in shape and pigmentation pattern in Otocinclus cocama (Siluriformes: Loricariidae) across major river drainages

Otocinclus cocama, a uniquely colored species of the loricariid catfish genus Otocinclus described solely from the type locality in the lower Ucayali River in northern Peru, is reported occurring in the Tigre River, a tributary to the Marañón River that drains a different section of the Andean Mountain range in the western Amazon. Both populations differ in the number of dark bars spanning the flanks of the body, and we investigated whether these morphotypes constitute distinct species. The body shapes of populations from the Tigre and Ucayali rivers were compared using geometric morphometrics. Although principal component analysis detected a broad overlap between populations, multivariate analysis of variance and linear driscriminat analysis revealed a subtle differentiation between the populations of the two hydrographic basins. Average body shape of the Ucayali River population tend to be slightly higher than that of the Tigre River, with the caudal peduncle stretched vertically in the Ucayali population. Multivariate regression of shape and centroid size revealed an allometric effect of 10.7% (p < 0.001), suggesting that the variation between Tigre and Ucayali populations was purely shape variation. Molecular data of coI, cytb, nd2, and 16S mitochondrial genes indicated a nucleotide diversity range from 0.001 to 0.003, and haplotypic diversity range from 0.600 ± 0.11 to 0.79 ± 0.07. The median-joining haplotype network for the concatenated matrix exhibited two divergent haplogroups related to the geographic area and separated by <10 mutational steps. The molecular species delimitation methods based on distance (automatic barcode gap discovery and assemble species by automatic partitioning) recovered two molecular lineages evolving independently, being one of the lineages formed by individuals from both populations. Tree-based methods (generalized mixed Yule coalescent and Bayesian implementation of the Poisson tree process) recovered similar topologies and supported single lineage recognition. Methods of molecular delimitation of species disclosed the high similarity between the two populations of Otocinclus cocama, further supported by the presence of old haplotypes common to both groups which could indicate that the populations still maintain gene flow. Although the morphological data reveal a subtle variation between both river basins, the molecular data suggest a weak population structuration based on hydrographic areas, but not different species lineages, therefore Otocinclus cocama is composed of a single lineage with two distinct morphotypes.

Recent genetic, phenetic and ecological divergence across the Mesoamerican highlands: a study case with Diglossa baritula (Aves: Thraupidae)

The topographical, geological, climatic and biodiversity complexity of Mesoamerica has made it a primary research focus. The Mesoamerican highlands is a region with particularly high species richness and within-species variation. The Cinnamon-bellied Flowerpiercer, Diglossa baritula (Wagler, 1832), is a species endemic to the Mesoamerican highlands, with three allopatric subspecies currently recognized. To characterize divergence within this species, we integrated genomics, morphology, coloration and ecological niche modeling approaches, obtained from sampling individuals across the entire geographic distribution of the species. Our results revealed a clear genomic divergence between the populations to the east versus the west of the Isthmus of Tehuantepec. In contrast to the genomic results, morphology and coloration analyses showed intermediate levels of differentiation, indicating that population groups within D. baritula have probably been under similar selective pressures. Our morphology results indicated that the only sexually dimorphic morphological variable is the wing chord, with males having a longer wing chord than females. Finally, ecological data indicated that there are differences in ecological niche within D. baritula. Our data suggest that D. baritula could contain two or more incipient species at the intermediate phase of the speciation continuum. These results highlight the importance of the geographical barrier of the Isthmus of Tehuantepec and Pleistocene climatic events in driving isolation and population divergence in D. baritula. The present investigation illustrates the speciation potential of the D. baritula complex and the capacity of Mesoamerican highlands to create cryptic biodiversity and endemism.

Biogeographic factors contributing to the diversification of Euphoniinae (Aves, Passeriformes, Fringillidae): a phylogenetic and ancestral areas analysis

Factors such as the Andean uplift, Isthmus of Panama, and climate changes have influenced bird diversity in the Neotropical region. Studying bird species that are widespread in Neotropical highlands and lowlands can help us understand the impact of these factors on taxa diversification. Our main objectives were to determine the biogeographic factors that contributed to the diversification of Euphoniinae and re-evaluate their phylogenetic relationships. The nextRAD and mitochondrial data were utilized to construct phylogenies. The ancestral distribution range was then estimated using a time-calibrated phylogeny, current species ranges, and neotropical regionalization. The phylogenies revealed two main Euphoniinae clades, Chlorophonia and Euphonia, similar to previous findings. Furthermore, each genus has distinctive subclades corresponding to morphology and geography. The biogeographic results suggest that the Andean uplift and the establishment of the western Amazon drove the vicariance of Chlorophonia and Euphonia during the Miocene. The Chlorophonia lineage originated in the Andes mountains and spread to Central America and the Mesoamerican highlands after the formation of the Isthmus of Panama. Meanwhile, the ancestral area of Euphonia was the Amazonas, from which it spread to trans-Andean areas during the Pliocene and Pleistocene due to the separation of the west lowlands from Amazonas due to the Northern Andean uplift. Chlorophonia and Euphonia species migrated to the Atlantic Forest during the Pleistocene through corridors from the East Andean Humid Forest and Amazonas. These two genera had Caribbean invasions with distinct geographic origins and ages. Finally, we suggested taxonomic changes in the genus Euphonia based on the study's phylogenetic, morphological, and biogeographic findings.

Detecting Episodic Evolution through Bayesian Inference of Molecular Clock Models

Molecular evolutionary rate variation is a key aspect of the evolution of many organisms that can be modeled using molecular clock models. For example, fixed local clocks revealed the role of episodic evolution in the emergence of SARS-CoV-2 variants of concern. Like all statistical models, however, the reliability of such inferences is contingent on an assessment of statistical evidence. We present a novel Bayesian phylogenetic approach for detecting episodic evolution. It consists of computing Bayes factors, as the ratio of posterior and prior odds of evolutionary rate increases, effectively quantifying support for the effect size. We conducted an extensive simulation study to illustrate the power of this method and benchmarked it to formal model comparison of a range of molecular clock models using (log) marginal likelihood estimation, and to inference under a random local clock model. Quantifying support for the effect size has higher sensitivity than formal model testing and is straight-forward to compute, because it only needs samples from the posterior and prior distribution. However, formal model testing has the advantage of accommodating a wide range molecular clock models. We also assessed the ability of an automated approach, known as the random local clock, where branches under episodic evolution may be detected without their a priori definition. In an empirical analysis of a data set of SARS-CoV-2 genomes, we find "very strong" evidence for episodic evolution. Our results provide guidelines and practical methods for Bayesian detection of episodic evolution, as well as avenues for further research into this phenomenon.

Propagating uncertainty about molecular evolution models and prior distributions to phylogenetic trees

Phylogenetic trees constructed from molecular sequence data rely on largely arbitrary assumptions about the substitution model, the distribution of substitution rates across sites, the version of the molecular clock, and, in the case of Bayesian inference, the prior distribution. Those assumptions affect results reported in the form of clade probabilities and error bars on divergence times and substitution rates. Overlooking the uncertainty in the assumptions leads to overly confident conclusions in the form of inflated clade probabilities and short confidence intervals or credible intervals. This paper demonstrates how to propagate that uncertainty by combining the models considered along with all of their assumptions, including their prior distributions. The combined models incorporate much more of the uncertainty than Bayesian model averages since the latter tend to settle on a single model due to the higher-level assumption that one of the models is true. Nucleotide sequence data illustrates the proposed model combination method.

The evolution, characterization and phylogeography of avian influenza H9N2 viruses from India

The low pathogenic avian influenza H9N2 virus is a significant zoonotic agent and contributes genes to highly pathogenic avian influenza (HPAI) viruses. H9N2 viruses are prevalent in India with a reported human case. We elucidate the spatio-temporal origins of the H9N2 viruses from India. A total of 30H9N2 viruses were isolated from poultry and environmental specimens (years 2015-2020). Genome sequences of H9N2 viruses (2003-2020) from India were analyzed, revealing several substitutions. We found five reassortant genotypes. The HA, NA and PB2 genes belonged to the Middle-Eastern B sublineage; NP and M to the classical G1 lineage; PB1, PA and NS showed resemblance to genes from either HPAI-H7N3/H5N1 viruses. Molecular clock and phylogeography revealed that the introduction of all the genes to India took place around the year 2000. This is the first report of the genesis and evolution of the H9N2 viruses from India, and highlights the need for surveillance.

Evolutionary features of a prolific subtype of avian influenza A virus in European waterfowl

Avian influenza A virus (AIV) is ubiquitous in waterfowl and is detected annually at high prevalence in waterfowl during the Northern Hemisphere autumn. Some AIV subtypes are globally common in waterfowl, such as H3N8, H4N6, and H6N2, and are detected in the same populations at a high frequency, annually. In order to investigate genetic features associated to the long-term maintenance of common subtypes in migratory ducks, we sequenced 248 H4 viruses isolated across 8 years (2002-9) from mallards (Anas platyrhynchos) sampled in southeast Sweden. Phylogenetic analyses showed that both H4 and N6 sequences fell into three distinct lineages, structured by year of isolation. Specifically, across the 8 years of the study, we observed lineage replacement, whereby a different HA lineage circulated in the population each year. Analysis of deduced amino acid sequences of the HA lineages illustrated key differences in regions of the globular head of hemagglutinin that overlap with established antigenic sites in homologous hemagglutinin H3, suggesting the possibility of antigenic differences among these HA lineages. Beyond HA, lineage replacement was common to all segments, such that novel genome constellations were detected across years. A dominant genome constellation would rapidly amplify in the duck population, followed by unlinking of gene segments as a result of reassortment within 2-3 weeks following introduction. These data help reveal the evolutionary dynamics exhibited by AIV on both annual and decadal scales in an important reservoir host.

Australia as a global sink for the genetic diversity of avian influenza A virus

Most of our understanding of the ecology and evolution of avian influenza A virus (AIV) in wild birds is derived from studies conducted in the northern hemisphere on waterfowl, with a substantial bias towards dabbling ducks. However, relevant environmental conditions and patterns of avian migration and reproduction are substantially different in the southern hemisphere. Through the sequencing and analysis of 333 unique AIV genomes collected from wild birds collected over 15 years we show that Australia is a global sink for AIV diversity and not integrally linked with the Eurasian gene pool. Rather, AIV are infrequently introduced to Australia, followed by decades of isolated circulation and eventual extinction. The number of co-circulating viral lineages varies per subtype. AIV haemagglutinin (HA) subtypes that are rarely identified at duck-centric study sites (H8-12) had more detected introductions and contemporary co-circulating lineages in Australia. Combined with a lack of duck migration beyond the Australian-Papuan region, these findings suggest introductions by long-distance migratory shorebirds. In addition, on the available data we found no evidence of directional or consistent patterns in virus movement across the Australian continent. This feature corresponds to patterns of bird movement, whereby waterfowl have nomadic and erratic rainfall-dependant distributions rather than consistent intra-continental migratory routes. Finally, we detected high levels of virus gene segment reassortment, with a high diversity of AIV genome constellations across years and locations. These data, in addition to those from other studies in Africa and South America, clearly show that patterns of AIV dynamics in the Southern Hemisphere are distinct from those in the temperate north.

Phylogeography and evolutionary dynamics analysis of porcine delta-coronavirus with host expansion to humans

From 2003 onwards, three pandemics have been caused by coronaviruses: severe acute respiratory syndrome coronavirus (SARS-CoV); middle east respiratory syndrome coronavirus (MERS-CoV); and, most recently, SARS-CoV-2. Notably, all three were transmitted from animals to humans. This would suggest that animals are potential sources of epidemics for humans. The emerging porcine delta-coronavirus was reported to infect children. This is a red flag that marks the ability of PDCoV to break barriers of cross-species transmission to humans. Therefore, we conducted molecular genetic analysis of global clade PDCoV to characterize spatio-temporal patterns of viral diffusion and genetic diversity. PDCoV was classified into three major lineages, according to distribution and phylogenetic analysis of PDCoV. It can be inferred based on the analysis results of the currently known PDCoV strains that PDCoV might originate in Asia. We also selected six special spike amino acid sequences to align and analyze to find seven significant mutation sites. The accumulation of these mutations may enhance dynamic movements, accelerating spike protein membrane fusion events and transmission. Altogether, our study offers a novel insight into the diversification, evolution, and interspecies transmission and origin of PDCoV and emphasizes the need to study the zoonotic potential of the PDCoV and comprehensive surveillance and enhanced biosecurity precautions for PDCoV. This article is protected by copyright. All rights reserved.

Chilli leaf curl disease populations in India are highly recombinant, and rapidly segregated

Capsicum annuum, a valuable spice and vegetable crop belonging to the Solanaceae family, is extensively grown across the Indian subcontinent. Chilli production is restricted by a begomoviral infection named as chilli leaf curl disease (ChiLCD) mainly in tropical and subtropical regions which leads to considerable economic losses, thus affecting chilli cultivation. Here, we studied the genetic diversity with structural evaluation of chilli leaf curl disease and satellite molecules infecting Chilli in India. We retrieved 121 reference sequences of ChiLCD including DNA-A, DNA-B, beta-satellite and alpha-satellites from GenBank reported from India. The population diversity and genetic variation were estimated through various parameters which decipher the four major groups of phylogenetic divergence for DNA-A and five groups of beta-satellite showing percentage similarity with isolates within and across India. Further, transitional and transversional bias for ORFs were observed highest in C4 and REn genes, respectively, and for DNA-A and DNA-B, these values were 1.07 and 1.22, respectively. The recombination breakpoints for DNA-A were estimated 49 majorly in V1, C1,C2 and C4 genome region and highest 22 breakpoints were determined for Rep (AC1) of ORFs, similarly 9 events for beta-satellite were found less around βC1ORF. Moreover, the evolution and genetic variability were also contributed through parameters such as nucleotide substitution which were found within the range of RNA viruses for DNA-A, DNA-B, for all 6 ORFs (relaxed clock) and beta-satellite. Additionally, total numbers of mutations (η) for DNA-A, DNA-B, alpha-satellites and beta-satellites were 2505, 419, 807 and 1288 detected, respectively, while it was found 987 highest for Rep gene among all ORFs. Further, neutrality tests determine the dominant nature of population expansion and purifying selection for all the genes of begomovirus associated with ChiLCD and satellite molecules supporting conserved nature of gene. The combined Tajima's D and Fu and Li'S D* negative values in tests indicated that population are under purified selection and an excess of low-frequency polymorphism. Our analysis indicates the potential contribution of genetic mutations and recombination of ChiLCD which leads to rapid adaptation and evolution of begomovirus and its satellite molecules accelerating its host range and diversity within and across the Indian subcontinent.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-022-03139-w.

Evolutionary Dynamics of Type 2 Porcine Reproductive and Respiratory Syndrome Virus by Whole-Genome Analysis

Porcine reproductive and respiratory syndrome virus (PRRSV), an important pathogen in the swine industry, is a genetically highly diverse RNA virus. However, the phylogenetic and genomic recombination properties of this virus are not yet fully understood. In this study, we performed an integrated analysis of all available whole-genome sequences of type 2 PRRSV (n = 901) to reveal its evolutionary dynamics. The results showed that there were three distinct phylogenetic lineages of PRRSV in their distribution patterns. We identified that sublineage 2.7 (L2.7), associated with a NADC30 cluster, had the highest substitution rate and higher viral genetic diversity, and inter-lineage recombination is observed more frequently in L2.7 PRRSV compared to other sublineages. Most inter-lineage recombination events detected are observed between L2.7 PRRSVs (as major parents) and L3.4 (a JXA1-R-related cluster)/L3.7 (a WUH3-related cluster) PRRSVs (as minor parents). Moreover, the recombination hotspots are located in the structural protein gene ORF2 and ORF4, or in the non-structural protein gene nsp7. In addition, a GM2-related cluster, L3.2, shows inconsistent recombination modes compared to those of L2.7, suggesting that it may have undergone extensive and unique recombination in their evolutionary history. We also identified several amino acids under positive selection in GP2, GP4 and GP5, the major glycoproteins of PRRSV, showing the driving force behind adaptive evolution. Taken together, our results provide new insights into the evolutionary dynamics of PPRSV that contribute to our understanding of the critical factors involved in its evolution and guide future efforts to develop effective preventive measures against PRRSV.

Trends and Intensity of Rhinovirus Invasions in Kilifi, Coastal Kenya, Over a 12-Year Period, 2007-2018

Background

Rhinoviruses (RVs) are ubiquitous pathogens and the principal etiological agents of common cold. Despite the high frequency of RV infections, data describing their long-term epidemiological patterns in a defined population remain limited.

Methods

Here, we analyzed 1070 VP4/VP2 genomic region sequences sampled at Kilifi County Hospital on the Kenya coast. The samples were collected between 2007 and 2018 from hospitalized pediatric patients (<60 months of age) with acute respiratory illness.

Results

Of 7231 children enrolled, RV was detected in 1497 (20.7%) and VP4/VP2 sequences were recovered from 1070 samples (71.5%). A total of 144 different RV types were identified (67 Rhinovirus A, 18 Rhinovirus B, and 59 Rhinovirus C) and at any month, several types co-circulated with alternating predominance. Within types, multiple genetically divergent variants were observed. Ongoing RV infections through time appeared to be a combination of (1) persistent types (observed up to 7 consecutive months), (2) reintroduced genetically distinct variants, and (3) new invasions (average of 8 new types annually).

Conclusions

Sustained RV presence in the Kilifi community is mainly due to frequent invasion by new types and variants rather than continuous transmission of locally established types/variants.

Propagating clade and model uncertainty to confidence intervals of divergence times and branch lengths

Confidence intervals of divergence times and branch lengths do not reflect uncertainty about their clades or about the prior distributions and other model assumptions on which they are based. Uncertainty about the clade may be propagated to a confidence interval by multiplying its confidence level by the bootstrap proportion of its clade or by another probability that the clade is correct. (If the confidence level is 95% and the bootstrap proportion is 90%, then the uncertainty-adjusted confidence level is (0.95)(0.90) = 86%.) Uncertainty about the model can be propagated to the confidence interval by reporting the union of the confidence intervals from all the plausible models. Unless there is no overlap between the confidence intervals, that results in an uncertainty-adjusted interval that has as its lower and upper limits the most extreme limits of the models. The proposed methods of uncertainty quantification may be used together.

Complete chloroplast genome of the medicinal plant Evolvulus alsinoides: comparative analysis, identification of mutational hotspots and evolutionary dynamics with species of Solanales

Evolvulus alsinoides, belonging to the family Convolvulaceae, is an important medicinal plant widely used as a nootropic in the Indian traditional medicine system. In the genus Evolvulus, no research on the chloroplast genome has been published. Hence, the present study focuses on annotation, characterization, identification of mutational hotspots, and phylogenetic analysis in the complete chloroplast genome (cp) of E. alsinoides. Genome comparison and evolutionary dynamics were performed with the species of Solanales. The cp genome has 114 genes (80 protein-coding genes, 30 transfer RNA, and 4 ribosomal RNA genes) that were unique with total genome size of 157,015 bp. The cp genome possesses 69 RNA editing sites and 44 simple sequence repeats (SSRs). Predicted SSRs were randomly selected and validated experimentally. Six divergent hotspots such as trnQ-UUG, trnF-GAA, psaI, clpP, ndhF, and ycf1 were discovered from the cp genome. These microsatellites and divergent hot spot sequences of the Taxa 'Evolvulus' could be employed as molecular markers for species identification and genetic divergence investigations. The LSC area was found to be more conserved than the SSC and IR region in genome comparison. The IR contraction and expansion studies show that nine genes rpl2, rpl23, ycf1, ycf2, ycf1, ndhF, ndhA, matK, and psbK were present in the IR-LSC and IR-SSC boundaries of the cp genome. Fifty-four protein-coding genes in the cp genome were under negative selection pressure, indicating that they were well conserved and were undergoing purifying selection. The phylogenetic analysis reveals that E. alsinoides is closely related to the genus Cressa with some divergence from the genus Ipomoea. This is the first time the chloroplast genome of the genus Evolvulus has been published. The findings of the present study and chloroplast genome data could be a valuable resource for future studies in population genetics, genetic diversity, and evolutionary relationship of the family Convolvulaceae.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01051-w.

Divergence time and species delimitation of microbivalves in the Southern Ocean: the case of Kidderia species

The systematics of Subantarctic and Antarctic near-shore marine benthic invertebrates requires major revision and highlights the necessity to incorporate additional sources of information in the specimen identification chart in the Southern Ocean (SO). In this study, we aim to improve our understanding of the biodiversity of Kidderia (Dall 1876) through molecular and morphological comparisons of Antarctic and Subantarctic taxa. The microbivalves of the genus Kidderia are small brooding organisms that inhabit intertidal and shallow subtidal rocky ecosystems. This genus represents an interesting model to test the vicariance and dispersal hypothesis in the biogeography of the SO. However, the description of Kidderia species relies on a few morphological characters and biogeographic records that raise questions about the true diversity in the group. Here we will define the specimens collected with genetic tools, delimiting their respective boundaries across provinces of the SO, validating the presence of two species of Kidderia. Through the revision of taxonomic issues and species delimitation, it was possible to report that the Antarctic species is Kidderia subquadrata and the species recorded in the Subantarctic islands Diego Ramirez, South Georgia and the Kerguelen Archipelago is Kidderia minuta. The divergence time estimation suggests the origin and diversification of Kidderia lineages are related to historical vicariant processes probably associated with the separation of the continental landmasses close to the late Eocene.

Paleogenomics illuminates the evolutionary history of the extinct Holocene "horned" crocodile of Madagascar, Voay robustus

Ancient DNA is transforming our ability to reconstruct historical patterns and mechanisms shaping modern diversity and distributions. In particular, molecular data from extinct Holocene island faunas have revealed surprising biogeographic scenarios. Here, we recovered partial mitochondrial (mt) genomes for 1300-1400 year old specimens (n = 2) of the extinct "horned" crocodile, Voay robustus, collected from Holocene deposits in southwestern Madagascar. Phylogenetic analyses of partial mt genomes and tip-dated timetrees based on molecular, fossil, and stratigraphic data favor a sister group relationship between Voay and Crocodylus (true crocodiles). These well supported trees conflict with recent morphological systematic work that has consistently placed Voay within Osteolaeminae (dwarf crocodiles and kin) and provide evidence for likely homoplasy in crocodylian cranial anatomy and snout shape. The close relationship between Voay and Crocodylus lends additional context for understanding the biogeographic origins of these genera and refines competing hypotheses for the recent extinction of Voay from Madagascar.

Molecular Phylogeny and Evolution of Amazon Parrots in the Greater Antilles

Amazon parrots (Amazona spp.) colonized the islands of the Greater Antilles from the Central American mainland, but there has not been a consensus as to how and when this happened. Today, most of the five remaining island species are listed as endangered, threatened, or vulnerable as a consequence of human activity. We sequenced and annotated full mitochondrial genomes of all the extant Amazon parrot species from the Greater Antillean (A. leucocephala (Cuba), A. agilis, A. collaria (both from Jamaica), A. ventralis (Hispaniola), and A. vittata (Puerto Rico)), A. albifrons from mainland Central America, and A. rhodocorytha from the Atlantic Forest in Brazil. The assembled and annotated mitogenome maps provide information on sequence organization, variation, population diversity, and evolutionary history for the Caribbean species including the critically endangered A. vittata. Despite the larger number of available samples from the Puerto Rican Parrot Recovery Program, the sequence diversity of the A. vittata population in Puerto Rico was the lowest among all parrot species analyzed. Our data support the stepping-stone dispersal and speciation hypothesis that has started approximately 3.47 MYA when the ancestral population arrived from mainland Central America and led to diversification across the Greater Antilles, ultimately reaching the island of Puerto Rico 0.67 MYA. The results are presented and discussed in light of the geological history of the Caribbean and in the context of recent parrot evolution, island biogeography, and conservation. This analysis contributes to understating evolutionary history and empowers subsequent assessments of sequence variation and helps design future conservation efforts in the Caribbean.

Adaptive dating and fast proposals: Revisiting the phylogenetic relaxed clock model

Relaxed clock models enable estimation of molecular substitution rates across lineages and are widely used in phylogenetics for dating evolutionary divergence times. Under the (uncorrelated) relaxed clock model, tree branches are associated with molecular substitution rates which are independently and identically distributed. In this article we delved into the internal complexities of the relaxed clock model in order to develop efficient MCMC operators for Bayesian phylogenetic inference. We compared three substitution rate parameterisations, introduced an adaptive operator which learns the weights of other operators during MCMC, and we explored how relaxed clock model estimation can benefit from two cutting-edge proposal kernels: the AVMVN and Bactrian kernels. This work has produced an operator scheme that is up to 65 times more efficient at exploring continuous relaxed clock parameters compared with previous setups, depending on the dataset. Finally, we explored variants of the standard narrow exchange operator which are specifically designed for the relaxed clock model. In the most extreme case, this new operator traversed tree space 40% more efficiently than narrow exchange. The methodologies introduced are adaptive and highly effective on short as well as long alignments. The results are available via the open source optimised relaxed clock (ORC) package for BEAST 2 under a GNU licence (https://github.com/jordandouglas/ORC).

Factors influencing estimates of HIV-1 infection timing using BEAST

While large datasets of HIV-1 sequences are increasingly being generated, many studies rely on a single gene or fragment of the genome and few comparative studies across genes have been done. We performed genome-based and gene-specific Bayesian phylogenetic analyses to investigate how certain factors impact estimates of the infection dates in an acute HIV-1 infection cohort, RV217. In this cohort, HIV-1 diagnosis corresponded to the first RNA positive test and occurred a median of four days after the last negative test, allowing us to compare timing estimates using BEAST to a narrow window of infection. We analyzed HIV-1 sequences sampled one week, one month and six months after HIV-1 diagnosis in 39 individuals. We found that shared diversity and temporal signal was limited in acute infection, and insufficient to allow timing inferences in the shortest HIV-1 genes, thus dated phylogenies were primarily analyzed for env, gag, pol and near full-length genomes. There was no one best-fitting model across participants and genes, though relaxed molecular clocks (73% of best-fitting models) and the Bayesian skyline (49%) tended to be favored. For infections with single founders, the infection date was estimated to be around one week pre-diagnosis for env (IQR: 3–9 days) and gag (IQR: 5–9 days), whilst the genome placed it at a median of 10 days (IQR: 4–19). Multiply-founded infections proved problematic to date. Our ability to compare timing inferences to precise estimates of HIV-1 infection (within a week) highlights that molecular dating methods can be applied to within-host datasets from early infection. Nonetheless, our results also suggest caution when using uniform clock and population models or short genes with limited information content.

Molecular dating using phylogenetics allows us to estimate the date of an infection from time-stamped within-host sequences alone. There are large datasets of HIV-1 sequences, but genome and gene analyses are not often performed in parallel and rarely with the possibility to compare results against a known narrow window of infection. We showed that all but the longest genes are near-clonal in acute infection, with little information for dating purposes. For infections with single founders, we estimated the eclipse phase—the time between HIV-1 exposure and the first positive diagnostic test—to last between one and two weeks using env, gag, pol and near full-length genomes. This approach could be used to narrow the date of suspected infection in ongoing clinical trials for the prevention of HIV-1 infection.

Detection of Ehrlichia sp. in Amblyomma sculptum parasitizing horses from Brazilian Pantanal wetland

Bacteria of the genus Ehrlichia are transmitted by ticks and also are an important cause of infection in wild and domestic mammals. Infection with Ehrlichia spp. has been reported in horses, especially in the USA, Nicaragua and Brazil. In this study, we report the parasitism by Amblyomma sculptum, Rhipicephalus microplus and Dermacentor nitens ticks in horses from a ranch located in south Pantanal wetland. Molecular and serological analyzes to determine infection by Ehrlichia spp. in horses and their respective ticks were carried out. A total of 12 horses were submitted to blood collection to investigate antibodies by indirect immunofluorescence assay (IFA) using Ehrlichia canis crude antigens and to be tested by polymerase chain reaction (PCR) in order to amplify fragments of the 16S rRNA, dsb, groEL and sodB gene of Ehrlichia spp. A total of 164 tick specimens were removed from horses, stored in isopropanol and later identified as D. nitens, A. sculptum and R. microplus. DNA from ticks were extracted and subjected to the same PCR assays to detect Ehrlichia spp. Anti-Ehrlichia spp. antibodies were detected in five/12 (41.7 %) horses by IFA, with antibody titers ranging from 40 to 160. All horse DNA samples were negative for the 16S rRNA, dsb, groEL and sodB of Ehrlichia spp. One A. sculptum female was positive to all target genes of Ehrlichia. This tick was parasitizing an Ehrlichia-seropositive horse with antibody titer of 80. Nucleotide sequences of 16S rRNA, dsb, groEL and sodB genes showed close relationship with different strains of Ehrlichia detected in wild mammals, Amblyomma ticks and horses from Brazil and Argentina. Detection of anti-Ehrlichia sp. antibodies suggests that horses have been exposed to an ehrlichial agent in the Pantanal. Future studies on Ehrlichia infection should be carried out to better elucidate and to bring new information about equine ehrlichiosis, since these animals are important hosts of ticks in the Brazilian Pantanal wetlands.

Hydrography rather than lip morphology better explains the evolutionary relationship between Gymnogeophagus labiatus and G. lacustris in Southern Brazil (Cichlidae: Geophagini)

ABSTRACT Gymnogeophagus labiatus and G. lacustris have been long recognized as sister species exhibiting different ecological requirements. Gymnogeophagus labiatus occurs in rock bottom rivers in the hydrographic basins of Patos Lagoon (HBP) and Tramandaí River (HBT), while G. lacustris is exclusive from sand bottom coastal lagoons of the HBT. In this study, we used molecular markers, morphological measurements and data from nuptial male coloration to investigate the evolutionary relationship between these species in each hydrographic basin. We found, for all data sets, a closer relationship between G. labiatus and G. lacustris from the HBT than between G. labiatus populations from HBT and HBP. In particular, lip area had a large intraspecific plasticity, being uninformative to diagnose G. lacustris from G. labiatus. Molecular clock-based estimates suggest a recent divergence between species in the HBT (17,000 years ago), but not between G. labiatus from HBP and HBT (3.6 millions of years ago). Finally, we also found a divergent G. labiatus genetic lineage from the Camaquã River, in the HBP. These results show that the current taxonomy of G. labiatus and G. lacustris does not properly represent evolutionary lineages in these species.

Miocene geologic dynamics of the Australian Sahul Shelf determined the biogeographic patterns of freshwater planorbid snails (Miratestinae) in the Indo-Australian Archipelago

The complex geological and climatic processes that have shaped the Indo-Australian Archipelago since the Cenozoic likely also gave rise to its species-rich biota. Strictly freshwater organisms might be particularly suitable for understanding the influence of these abiotic factors on their biogeography in such a insular setting as their distribution may reflect past abiotic events at large and small geographical scales. We here investigate the historical biogeography of the Miratestinae, a subfamily of Planorbidae. These freshwater gastropods are widely distributed in the eastern IAA from Australia, New Guinea, the Moluccas, and Sulawesi to the Philippines. The first comprehensive molecular phylogeny of the Miratestinae was inferred based on two mitochondrial and two nuclear genetic markers using maximum likelihood and Bayesian inference. Four species delimitation methods were applied to identify molecular operational taxonomic units (MOTUs). Divergence times were inferred using an uncorrelated lognormal relaxed-clock model by applying a taxon- and marker-specific substitution rate. Ancestral geographic ranges were estimated based on the dated phylogeny using BioGeoBEARS. The species delimitation revealed a total of 23 MOTUs, 16 of which might represent species new to science. The BioGeoBEARS analyses suggest an Australian origin for the Miratestinae at c. 22 Ma and identified jump dispersal to be the main process of colonization. The first colonization events from Australia to the IAA occurred in the Middle-Late Miocene (12-13 Ma), whereas intra-island diversification took mainly place since the Late Miocene-Pliocene. Colonization and diversification events remarkably coincide with major geologic events that shaped the geography of the region. The increasing availability of landmasses along the Sahul Shelf likely promoted stepping-stone dispersal to New Guinea, Sulawesi and the Philippines as early as the islands emerged. Major geological and climatic events such as the amalgamation of the island Sulawesi, the regional aridification in Australia or the uplift of massive mountain ranges in New Guinea likely played a considerable role for intra-island diversification.

Phylogeographic analysis and species distribution modelling of the wood frog Batrachyla leptopus (Batrachylidae) reveal interglacial diversification in south western Patagonia

Background

The evolutionary history of southern South American organisms has been strongly influenced by Pleistocene climate oscillations. Amphibians are good models to evaluate hypotheses about the influence of these climate cycles on population structure and diversification of the biota, because they are sensitive to environmental changes and have restricted dispersal capabilities. We test hypotheses regarding putative forest refugia and expansion events associated with past climatic changes in the wood frog Batrachyla leptopus distributed along ∼1,000 km of length including glaciated and non-glaciated areas in southwestern Patagonia.

Methods

Using three mitochondrial regions (D-loop, cyt b, and coI) and two nuclear loci (pomc and crybA1), we conducted multilocus phylogeographic analyses and species distribution modelling to gain insights of the evolutionary history of this species. Intraspecific genealogy was explored with maximum likelihood, Bayesian, and phylogenetic network approaches. Diversification time was assessed using molecular clock models in a Bayesian framework, and demographic scenarios were evaluated using approximate Bayesian computation (ABC) and extended Bayesian skyline plot (EBSP). Species distribution models (SDM) were reconstructed using climatic and geographic data.

Results

Population structure and genealogical analyses support the existence of four lineages distributed north to south, with moderate to high phylogenetic support (Bootstrap > 70%; BPP > 0.92). The diversification time of B. leptopus’ populations began at ∼0.107 mya. The divergence between A and B lineages would have occurred by the late Pleistocene, approximately 0.068 mya, and divergence between C and D lineages was approximately 0.065 mya. The ABC simulations indicate that lineages coalesced at two different time periods, suggesting the presence of at least two glacial refugia and a postglacial colonization route that may have generated two southern lineages (p = 0.93, type I error: <0.094, type II error: 0.134). EBSP, mismatch distribution and neutrality indexes suggest sudden population expansion at ∼0.02 mya for all lineages. SDM infers fragmented distributions of B. leptopus associated with Pleistocene glaciations. Although the present populations of B. leptopus are found in zones affected by the last glacial maximum (∼0.023 mya), our analyses recover an older history of interglacial diversification (0.107–0.019 mya). In addition, we hypothesize two glacial refugia and three interglacial colonization routes, one of which gave rise to two expanding lineages in the south.

Survey of vector-borne and nematode parasites involved in the etiology of anemic syndrome in sheep from Southern Brazil

Although anemia has been historically linked to Haemonchus contortus infection, other infectious agents, such as hemotropic mycoplasmas and tick-borne disease pathogens, may also lead to anemic crisis in sheep. This study has aimed to investigate infections related to anemia in a sheep herd from Bandeirantes City, Paraná State, southern Brazil. Seven out of forty-two (16.6%; 95% CI: 8.32-30.6%) sheep were positive for hemoplasmas by a PCR targeting the 16S rRNA gene and all tested negative for A. marginale/A. ovis and Babesia/Theileria spp. by PCR based on msp4 and 18S rRNA genes, respectively. Two (4.7%; 95% CI: 1.32-15.79%) animals were infested with Rhipicephalus microplus ticks. Fecal egg counting was performed in 38 sheep and 24 (63.15%; 95% CI: 47.2-76.6%) presented > 500 eggs per gram. Phylogenetic analysis of partial sequences of the detected hemotropic Mycoplasma sp. 16S and 23S rRNA genes confirmed that the animals were infected with Mycoplasma ovis. Polymorphism analysis of partial 16S rRNA sequences showed three different genotypes of M. ovis infecting sheep assessed in the present study. Mycoplasma ovis and gastrointestinal nematodes occurs in sheep from the northern region of Paraná State.

Dataset on the epidemiology and genetic diversification of dengue virus (DENV) serotypes and genotypes in Mexico

Dengue virus (DENV) evolution has had a significant impact on disease pathogenesis, virulence, and epidemiology in Mexico. Novel genotypic variation in DENV serotypes and genotypes may influence the magnitude and severity of dengue epidemics, as evidenced by 2009 data from Veracruz State. The data presented herein is related to the publication entitled "Epidemiological Implications of the Genetic Diversification of Dengue Virus (DENV) Serotypes and Genotypes in Mexico" [1]. Raw data and trees provide epidemiological data on DENV prevalence and a comprehensive phylogeny of both representative sequences collected from an NCBI repository, and 28 additional isolates from acute-phase plasma samples diagnosed with dengue fever or severe dengue (Raw sequencing data is hosted in the public repository Mendeley Data (http://dx.doi.org/10.17632/bf2kdhhf6x.2). Phylogenetic trees for each DENV serotype (DENV-1, -2, -3 and -4) were constructed using these sequences by a maximum likelihood methodology as well as a Bayesian Markov chain Monte Carlo (MCMC) integration approach. Phylogenetic trees exhibited: (1) DENV-1, genotype V, (2) the DENV-2 Asian/American and Asian II genotypes, (3) DENV-3, genotype III, and (4) DENV-4, genotype I. This data can be beneficial for future analyses on DENV serotype and genotype structure and the introduction of novel DENV genotype sequences in the Americas, for the further elucidation of dengue etiology.

Multilocus, phenotypic, behavioral, and ecological niche analyses provide evidence for two species within Euphonia affinis (Aves, Fringillidae)

The integration of genetic, morphological, behavioral, and ecological information in the analysis of species boundaries has increased, allowing integrative systematics that better reflect the evolutionary history of biological groups. In this context, the goal of this study was to recognize independent evolutionary lineages within Euphonia affinis at the genetic, morphological, and ecological levels. Three subspecies have been described: E. affinis godmani, distributed in the Pacific slope from southern Sonora to Guerrero; E. affinis affinis, from Oaxaca, Chiapas and the Yucatan Peninsula to Costa Rica; and E. affinis olmecorum from Tamaulipas and San Luis Potosi east to northern Chiapas (not recognized by some authors). A multilocus analysis was performed using mitochondrial and nuclear genes. These analyses suggest two genetic lineages: E. godmani and E. affinis, which diverged between 1.34 and 4.3 My, a period in which the ice ages and global cooling fragmented the tropical forests throughout the Neotropics. To analyze morphometric variations, six morphometric measurements were taken, and the Wilcoxon Test was applied to look for sexual dimorphism and differences between the lineages. Behavioral information was included, by performing vocalization analysis which showed significant differences in the temporal characteristics of calls. Finally, Ecological Niche Models were estimated with MaxEnt, and then compared using the method of Broennimann. These analyses showed that the lineage distributed in western Mexico (E. godmani) has a more restricted niche than the eastern lineage (E. affinis) and thus we rejected the hypotheses of niche equivalence and similarity. Based on the combined evidence from genetic, morphological, behavioral, and ecological data, it is concluded that E. affinis (with E. olmecorum as its synonym) and E. godmani represent two independent evolutionary lineages.

Post K-Pg diversification of the mammalian order Eulipotyphla as suggested by phylogenomic analyses of ultra-conserved elements

The origin of the mammalian order Eulipotyphla has been debated intensively with arguments around whether they began diversifying before or after the Cretaceous-Palaeogene (K-Pg) boundary at 66 Ma. Here, we used an in-solution nucleotide capture method and next generation DNA sequencing to determine the sequence of hundreds of ultra-conserved elements (UCEs), and conducted phylogenomic and molecular dating analyses for the four extant eulipotyphlan lineages-Erinaceidae, Solenodontidae, Soricidae, and Talpidae. Concatenated maximum-likelihood analyses with single or partitioned models and a coalescent species-tree analysis showed that divergences among the four major eulipotyphlan lineages occurred within a short period of evolutionary time, but did not resolve the interrelationships among them. Alternative suboptimal phylogenetic hypotheses received consistently the same amount of support from different UCE loci, and were not significantly different from the maximum likelihood tree topology, suggesting the prevalence of stochastic lineage sorting. Molecular dating analyses that incorporated among-lineage evolutionary rate differences supported a scenario where the four eulipotyphlan families diversified between 57.8 and 63.2 Ma. Given short branch lengths with low support values, traces of rampant genome-wide stochastic lineage sorting, and post K-Pg diversification, we concluded that the crown eulipotyphlan lineages arose through a rapid diversification after the K-Pg boundary when novel niches were created by the mass extinction of species.

Origin and evolution of papillomavirus (onco)genes and genomes

Papillomaviruses (PVs) are ancient viruses infecting vertebrates, from fishes to mammals. Although the genomes of PVs are small and show conserved synteny, PVs display large genotypic diversity and ample variation in the phenotypic presentation of the infection. Most PV genomes contain two small early genes E6 and E7. In a bunch of closely related human papillomaviruses (HPVs), the E6 and E7 proteins provide the viruses with oncogenic potential. The recent discoveries of PVs without E6 and E7 in different fish species place a new root on the PV tree, and suggest that ancestral PVs consisted of the minimal PV backbone E1-E2-L2-L1. Bayesian phylogenetic analyses date the most recent common ancestor of the PV backbone to 424 million years ago (Ma). Common ancestry tests on extant E6 and E7 genes indicate that they share a common ancestor dating back to at least 184 Ma. In AlphaPVs infecting Old World monkeys and apes, the appearance of the E5 oncogene 53-58 Ma concurred with (i) a significant increase in substitution rate, (ii) a basal radiation and (iii) key gain of functions in E6 and E7. This series of events was instrumental to construct the extant phenotype of oncogenic HPVs. Our results assemble the current knowledge on PV diversity and present an ancient evolutionary timeline punctuated by evolutionary innovations in the history of this successful viral family. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

Hidden species diversity in Pachyhynobius: A multiple approaches species delimitation with mitogenomes

The lack of distinct morphological features of cryptic species is a hard problem for taxonomy, especially when the taxa are closely related with considerable amounts of ancestral polymorphism. Lately, intensive coalescent-based analyses involving multiple loci have become the preferred method to assess the extent of genetic distinctiveness in otherwise phenotypically similar populations. Previously, phylogenetic studies on Pachyhynobius shangchengensis uncovered five extremely deeply divergent clades, which suggested that this species may be a cryptic species complex. In this study, we used the complete mitochondrial genome data and samples from the entire range of stout salamander (Pachyhynobius), as well as publicly available mitochondrial genomes to assess species boundaries within this genus using a suite of diverse methodologies (e.g. general mixed Yule coalescent model, Automatic Barcode Gap Discovery). The phylogenetic relationships recovered two major groups within P. shangchengensis, with one group formed by four of the six extant populations and corresponding to the central and eastern range of the Dabie mountains, while the other group encompassed two other lineages in the north west of the Dabie mountain range. The species delimitation comparison within Pachyhynobius supported the presence of recognized species within the genus, and consensus was observed across methods for the existence of up to five cryptic species within what has been traditionally considered to be P. shangchengensis. While this implies the existence of four taxa in addition to the described P. shangchengensis species, morphological data and life history information are further required to contribute to the species definition. The observed pattern of genetic variation is likely the outcome of a discontinuous habitat combined with niche conservatism, which produced the sky-island effect observed in Pachyhynobius, and which led to formation of a hidden species diversity in this genus.

Phylodynamic Model Adequacy Using Posterior Predictive Simulations

Rapidly evolving pathogens, such as viruses and bacteria, accumulate genetic change at a similar timescale over which their epidemiological processes occur, such that, it is possible to make inferences about their infectious spread using phylogenetic time-trees. For this purpose it is necessary to choose a phylodynamic model. However, the resulting inferences are contingent on whether the model adequately describes key features of the data. Model adequacy methods allow formal rejection of a model if it cannot generate the main features of the data. We present TreeModelAdequacy, a package for the popular BEAST2 software that allows assessing the adequacy of phylodynamic models. We illustrate its utility by analyzing phylogenetic trees from two viral outbreaks of Ebola and H1N1 influenza. The main features of the Ebola data were adequately described by the coalescent exponential-growth model, whereas the H1N1 influenza data were best described by the birth–death susceptible-infected-recovered model.

Unusual AT-skew of Sinorhodeus microlepis mitogenome provides new insights into mitogenome features and phylogenetic implications of bitterling fishes

Sinorhodeus microlepis (S. microlepis) is recently described as a new species and represents a new genus Sinorhodeu of the subfamily Acheilognathinae. In this study, we first sequenced the complete mitogenome of S. microlepis and compared with the other 29 bitterling mitogenomes. The S. microlepis mitogenome is 16,591 bp in length and contains 37 genes. Gene distribution pattern is identical among 30 bitterling mitogenomes. A significant linear correlation between A+T% and AT-skew were found among 29 bitterling mitogenomes, except S. microlepis shows unusual AT-skew with slightly negative in tRNAs and PCGs. Bitterling mitogenomes exhibit highly conserved usage bias of start codon, relative synonymous codons and amino acids, overlaps and non-coding intergenic spacers. Phylogenetic trees constructed by 13 PCGs strongly support the polyphyly of the genus Acheilognathus and the paraphyly of Rhodeus and Tanakia. Together with the unusual characters of S. microlepis mitogenomes and phylogenetic trees, S. microlepis should be a sister species to the genus Rhodeu that might diverge about 13.69 Ma (95% HPD: 12.96-14.48 Ma).

Population genetic structure and geographical variation in Neotricula aperta (Gastropoda: Pomatiopsidae), the snail intermediate host of Schistosoma mekongi (Digenea: Schistosomatidae)

Background

Neotricula aperta is the snail-intermediate host of the parasitic blood-fluke Schistosoma mekongi which causes Mekong schistosomiasis in Cambodia and the Lao PDR. Despite numerous phylogenetic studies only one DNA-sequence based population-genetic study of N. aperta had been published, and the origin, structure and persistence of N. aperta were poorly understood. Consequently, a phylogenetic and population genetic study was performed, with addition of new data to pre-existing DNA-sequences for N. aperta from remote and inaccessible habitats, including one new taxon from Laos and 505 bp of additional DNA-sequence for all sampled taxa,.

Principal findings

Spatial Principal Component Analysis revealed the presence of significant spatial-genetic clustering. Genetic-distance-based clustering indicated four populations with near perfect match to a priori defined ecogeographical regions. Spring-dwelling taxa were found to form an ecological isolate relative to other N. aperta. The poor dispersal capabilities suggested by spatial-genetic analyses were confirmed by Bayesian inference of migration rates. Population divergence time estimation implied a mid-Miocene colonisation of the present range, with immediate and rapid radiation in each ecogeographical region. Estimated effective population sizes were large (120–310 thousand).

Conclusions

The strong spatial-genetic structure confirmed the poor dispersal capabilities of N. aperta—suggesting human-mediated reintroduction of disease to controlled areas as the primary reason for control failure. The isolation of the spring-dwelling taxa and ecogeographical structure suggests adaptation of sub-populations to different habitats; the epidemiological significance of this needs investigation. The large effective population sizes indicate that the high population densities observed in surveyed habitats are also present in inaccessible areas; affording great potential for recrudescence driven by animal-reservoir transmission in remote streams. Mid-Miocene colonisation implies heterochronous evolution of these snails and associated schistosomes and suggests against coevolution of snail and parasite. Heterochronicity favours ecological factors as shapers of host-parasite specificity and greater potential for escape from schistosomiasis control through host-switching.

The disease Mekong schistosomiasis poses a threat to the health of about 1.5 million people living near the Mekong river and its tributaries in Cambodia and Laos. It is a water-borne parasite transmitted by direct contact with water in which freshwater snails of the species Neotricula aperta live. Control of the snails is an effective approach to control of the parasite; however, because many suitable habitats for N. aperta occur in remote and inaccessible areas, knowledge of N. aperta population sizes and interconnectivity is insufficient for the design of effective snail control interventions. Although much of the region is difficult to survey by conventional means, population genetics can be used to estimate population structure and total size from small samples of accessible populations. The study added to existing data-sets, to give more population samples and longer DNA-sequences, together with improved analytical approaches to provide a better overview of N. aperta. The findings suggest that N. aperta in different kinds of habitats are also genetically different, with very low levels of migration between them; this genetic clustering is greater than expected from spatial distance alone. Further work is needed to determine if these different clusters vary in ability to transmit the parasite. The overall population size estimates were very large; thus suggesting that high snail population densities observed in accessible habitats are also characteristic of inaccessible populations—parasites are therefore more likely to return after disease control by immigration from remote areas. Finally, the timing of evolutionary events for snails and parasites was found to differ; this implies that the parasite may not be as strongly restricted to one species of snail as originally thought, which has implications for avoidance of parasite control by host-switching.

Speciation in a biodiversity hotspot: Phylogenetic relationships, species delimitation, and divergence times of Patagonian ground frogs from the Eupsophus roseus group (Alsodidae)

The alsodid ground frogs of the Eupsophus genus are divided into two groups, the roseus (2n = 30) and vertebralis (2n = 28), which are distributed throughout the temperate Nothofagus forests of South America. Currently, the roseus group is composed by four species, while the vertebralis group consists of two. Phylogenetic relationships and species delimitation within each group are controversial. In fact, previous analyses considered that the roseus group was composed of between four to nine species. In this work, we evaluated phylogenetic relationships, diversification times, and species delimitation within the roseus group using a multi-locus dataset. For this purpose, mitochondrial (D-loop, Cyt b, and COI) and nuclear (POMC and CRYBA1) partial sequences from 164 individuals were amplified, representing all species. Maximum Likelihood (ML) and Bayesian approaches were used to reconstruct phylogenetic relationships. Species tree was estimated using BEAST and singular value decomposition scores for species quartets (SVDquartets). Species limits were evaluated with six coalescent approaches. Diversification times were estimated using mitochondrial and nuclear rates with LogNormal relaxed clock in BEAST. Nine well-supported monophyletic lineages were recovered in Bayesian, ML, and SVDquartets, including eight named species and a lineage composed by specimens from the Villarrica population (Bootstrap:>70, PP:> 0.99). Single-locus species delimitation analyses overestimated the species number in E. migueli, E. calcaratus, and E. roseus lineages, while multi-locus analyses recovered as species the nine lineages observed in phylogenetic analyses (Ctax = 0.69). It is hypothesized that Eupsophus diversification occurred during Mid-Pleistocene (0.42-0.14 Mya), with most species having originated after the Last Southern Patagonian Glaciation (0.18 Mya). Our results revitalize the hypothesis that the E. roseus group is composed of eight species and support the Villarrica lineage as a new putative species.

Morphometrical and molecular evidence suggests cryptic diversity among hookworms (Nematoda: Uncinaria) that parasitize pinnipeds from the south-eastern Pacific coasts

Hookworms of the genus Uncinaria parasitize pinniped pups in various locations worldwide. Four species have been described, two of which parasitize pinniped pups in the southern hemisphere: Uncinaria hamiltoni parasitizes Otaria flavescens and Arctocephalus australis from the South American coast, and Uncinaria sanguinis parasitizes Neophoca cinerea from the Australian coast. However, their geographical ranges and host specificity are unknown. Uncinaria spp. are morphologically similar, but molecular analyses have allowed the recognition of new species in the genus Uncinaria. We used nuclear genetic markers (internal transcribed spacer (ITS) and large subunit (LSU) rDNA) and a mitochondrial genetic marker (cytochrome c oxidase subunit I (COI)) to evaluate the phylogenetic relationships of Uncinaria spp. parasitizing A. australis and O. flavescens from South American coasts (Atlantic and Pacific coasts). We compared our sequences with published Uncinaria sequences. A Generalized Mixed Yule Coalescent (GMYC) analysis was also used to delimit species, and principal component analysis was used to compare morphometry among Uncinaria specimens. Parasites were sampled from A. australis from Peru (12°S), southern Chile (42°S), and the Uruguayan coast, and from O. flavescens from northern Chile (24°S) and the Uruguayan coast. Morphometric differences were observed between Uncinaria specimens from both South American coasts and between Uncinaria specimens from A. australis in Peru and southern Chile. Phylogenetic and GMYC analyses suggest that south-eastern Pacific otariid species harbour U. hamiltoni and an undescribed putative species of Uncinaria. However, more samples from A. australis and O. flavescens are necessary to understand the phylogenetic patterns of Uncinaria spp. across the South Pacific.

Posterior Summarization in Bayesian Phylogenetics Using Tracer 1.7

Bayesian inference of phylogeny using Markov chain Monte Carlo (MCMC) plays a central role in understanding evolutionary history from molecular sequence data. Visualizing and analyzing the MCMC-generated samples from the posterior distribution is a key step in any non-trivial Bayesian inference. We present the software package Tracer (version 1.7) for visualizing and analyzing the MCMC trace files generated through Bayesian phylogenetic inference. Tracer provides kernel density estimation, multivariate visualization, demographic trajectory reconstruction, conditional posterior distribution summary, and more. Tracer is open-source and available at http://beast.community/tracer.

Where do all the subtypes go? Temporal dynamics of H8-H12 influenza A viruses in waterfowl

Influenza A virus (IAV) is ubiquitous in waterfowl. In the northern hemisphere IAV prevalence is highest during the autumn and coincides with a peak in viral subtype diversity. Although haemagglutinin subtypes H1-H12 are associated with waterfowl hosts, subtypes H8-H12 are detected very infrequently. To better understand the role of waterfowl in the maintenance of these rare subtypes, we sequenced H8-H12 viruses isolated from Mallards (Anas platyrhynchos) from 2002 to 2009. These rare viruses exhibited varying ecological and phylodynamic features. The Eurasian clades of H8 and H12 phylogenies were dominated by waterfowl sequences; mostly viruses sequenced in this study. H11, once believed to be a subtype that infected charadriiformes (shorebirds), exhibited patterns more typical of common virus subtypes. Finally, subtypes H9 and H10, which have maintained lineages in poultry, showed markedly different patterns: H10 was associated with all possible NA subtypes and this drove HA lineage diversity within years. Rare viruses belonging to subtypes H8-H12 were highly reassorted, indicating that these rare subtypes are part of the broader IAV pool. Our results suggest that waterfowl play a role in the maintenance of these rare subtypes, but we recommend additional sampling of non-traditional hosts to better understand the reservoirs of these rare viruses.

Reappraisal of the systematics of Microglanis cottoides (Siluriformes, Pseudopimelodidae), a catfish from southern Brazil

The southern region of Brazil is characterized by high species diversity and endemism of freshwater fishes distributed across geographically isolated river basins. Microglanis cottoides has a widespread range across these river basins and occurs in sympatry with other endemic species of the genus (e.g. M. cibelae, M. eurystoma, and M. malabarbai). Herein we tested the monophyly of M. cottoides and presented for the first time information about the molecular phylogeny of species in the genus. The results suggest that M. cottoides currently forms a non-monophyletic group which includes populations endemic to the Uruguay River basin that are more closely related to M. malabarbai, and excludes M. cibelae, found to be nested within M. cottoides. Based on an integrative approach using morphological and molecular data, we propose M. cibelae as a junior synonym of M. cottoides, and the populations of the Uruguay River basin previously assigned to M. cottoides in fact belong to M. malabarbai. Our molecular phylogeny shows that M. cottoides is sister to M. parahybae, which is also a coastal species, and M. malabarbai is sister of M. garavelloi, both endemic to inland river basins. The time-calibrated phylogeny indicates that the separation between inland and the coastal clades occurred in the Tertiary period, and that the species within the coastal basins diverged in the Pliocene, which overlaps with the diversification times estimated for the two inland species as well. This pattern of diversification corroborates some previous studies with other fishes from the same region.

Does adaptation to vertebrate codon usage relate to flavivirus emergence potential?

Codon adaptation index (CAI) is a measure of synonymous codon usage biases given a usage reference. Through mutation, selection, and drift, viruses can optimize their replication efficiency and produce more offspring, which could increase the chance of secondary transmission. To evaluate how higher CAI towards the host has been associated with higher viral titers, we explored temporal trends of several historic and extensively sequenced zoonotic flaviviruses and relationships within the genus itself. To showcase evolutionary and epidemiological relationships associated with silent, adaptive synonymous changes of viruses, we used codon usage tables from human housekeeping and antiviral immune genes, as well as tables from arthropod vectors and vertebrate species involved in the flavivirus maintenance cycle. We argue that temporal trends of CAI changes could lead to a better understanding of zoonotic emergences, evolutionary dynamics, and host adaptation. CAI appears to help illustrate historically relevant trends of well-characterized viruses, in different viral species and genetic diversity within a single species. CAI can be a useful tool together with in vivo and in vitro kinetics, phylodynamics, and additional functional genomics studies to better understand species trafficking and viral emergence in a new host.