Structure, Evolution, and Mitochondrial Genome Analysis of Mussel Species (Bivalvia, Mytilidae).

Based on the nucleotide sequences of the mitochondrial genome (mitogenome) of specimens taken from two mussel species (Arcuatula senhousia and Mytilus coruscus), an investigation was performed by means of the complex approaches of the genomics, molecular phylogenetics, and evolutionary genetics. The mitogenome structure of studied mussels, like in many other invertebrates, appears to be much more variable than in vertebrates and includes changing gene order, duplications, and deletions, which were most frequent for tRNA genes; the mussel species' mitogenomes also have variable sizes. The results demonstrate some of the very important properties of protein polypeptides, such as hydrophobicity and its determination by the purine and pyrimidine nucleotide ratio. This fact might indirectly indicate the necessity of purifying natural selection for the support of polypeptide functionality. However, in accordance with the widely accepted and logical concept of natural cutoff selection for organisms living in nature, which explains its action against deleterious nucleotide substitutions in the nonsynonymous codons (mutations) and its holding of the active (effective) macromolecules of the polypeptides in a population, we were unable to get unambiguous evidence in favor of this concept in the current paper. Here, the phylogeny and systematics of mussel species from one of the largest taxons of bivalve mollusks are studied, the family known as Mytilidae. The phylogeny for Mytilidae (order Mytilida), which currently has no consensus in terms of systematics, is reconstructed using a data matrix of 26-27 mitogenomes. Initially, a set of 100 sequences from GenBank were downloaded and checked for their gender: whether they were female (F) or male (M) in origin. Our analysis of the new data confirms the known drastic differences between the F/M mitogenome lines in mussels. Phylogenetic reconstructions of the F-lines were performed using the combined set of genetic markers, reconstructing only protein-coding genes (PCGs), only rRNA + tRNA genes, and all genes. Additionally, the analysis includes the usage of nucleotide sequences composed of other data matrices, such as 20-68 mitogenome sequences. The time of divergence from MRCA, estimated via BEAST2, for Mytilidae is close to 293 Mya, suggesting that they originate in the Silurian Period. From all these data, a consensus for the phylogeny of the subfamily of Mytilinae and its systematics is suggested. In particular, the long-debated argument on mussel systematics was resolved as to whether Mytilidae, and the subfamily of Mytilinae, are monophyletic. The topology signal, which was strongly resolved in this paper and in the literature, has refuted the theory regarding the monophyly of Mytilinae.

Temporal evolution and global spread of hepatitis B virus genotype G.

Hepatitis B virus (HBV) infection is considered a major health problem in the world. HBV is classified into genotypes A to J disseminated worldwide. Genotypes A, D and F are the most frequent in the Western World, B and C are predominant in the East, and E, F, H and J are infrequent and restricted to specific regions. HBV-G is a rare genotype, but it has been detected in different continents. This study aimed to report the temporal evolution and global spread of HBV-G comparing whole-genome sequences of this genotype from different regions in the world. Bayesian coalescent analysis was performed to estimate the time to the most recent common ancestor (tMRCA) and the population dynamics in the last decades. The results demonstrated that tMRCA of all HBV-Gs dated back to 1855 (95% highest posterior density interval [HPD 95%]: 1778 - 1931). This genotype has a possible origin in North America and it was disseminated to other continents (South and Central America, Europe, Asia and Africa) more than one century later (around the 1970s). The viral population demonstrated constant spreading from 1855 to the 1980s, followed by an increase in the 1990s and reached a plateau after the 2000s. Wide spreading at the beginning of the 1990s was probably associated with the dissemination by highly sexual active groups and injecting drug users. In conclusion, the present study demonstrated that HBV-G was originated in the 19th century with main events of spread at the end of the 20th century.

Evolutionary history of hepatitis B virus genotype H.

Hepatitis B virus genotype H (HBV-H) molecular evolution was studied by comparing all published whole-genome sequences. Bayesian coalescent analysis was performed to estimate phylogenetic relationships, time to the most recent common ancestor (tMRCA), and viral population dynamics along the time. Phylogenetic tree demonstrated two main clades or lineages: HBV-H I (with sequences from Central and North America) and HBV-H II (with sequences from North and South America, and Asia). HBV-H II had more genome sequences (n = 26; 83.9%), including one specific subclade with all sequences outside of the Americas. Overall HBV-H tMRCA dated back to 1933 (95% highest posterior density interval [HPD 95%]: 1875-1957) with a very probable origin in Mexico and posterior dissemination to other American and Asian countries. The temporal analysis demonstrated that HBV-H I spread only in Mexico and the neighbor country of Nicaragua probably in the 1960s to the 1970s (1968; HPD 95%: 1908-1981), while HBV-II disseminated to other American and Asian countries around one decade later (1977; HPD 95%: 1925-1985). The phylogeographic analysis reinforced the Mexican origin of this genotype. The whole HBV-H population increased from the 1980s to the 2000s. In conclusion, HBV-H has two main lineages with a common origin in Mexico approximately nine decades ago.

Spatial and Temporal Scales of Range Expansion in Wild Phaseolus vulgaris.

The wild progenitor of common-bean has an exceptionally large distribution from northern Mexico to northwestern Argentina, unusual among crop wild progenitors. This research sought to document major events of range expansion that led to this distribution and associated environmental changes. Through the use of genotyping-by-sequencing (∼20,000 SNPs) and geographic information systems applied to a sample of 246 accessions of wild Phaseolus vulgaris, including 157 genotypes of the Mesoamerican, 77 of the southern Andean, and 12 of the Northern Peru-Ecuador gene pools, we identified five geographically distinct subpopulations. Three of these subpopulations belong to the Mesoamerican gene pool (Northern and Central Mexico, Oaxaca, and Southern Mexico, Central America and northern South America) and one each to the Northern Peru-Ecuador (PhI) and the southern Andean gene pools. The five subpopulations were distributed in different floristic provinces of the Neotropical seasonally dry forest and showed distinct distributions for temperature and rainfall resulting in decreased local potential evapotranspiration (PhI and southern Andes groups) compared with the two Mexican groups. Three of these subpopulations represent long-distance dispersal events from Mesoamerica into Northern Peru-Ecuador, southern Andes, and Central America and Colombia, in chronological order. Of particular note is that the dispersal to Northern Peru-Ecuador markedly predates the dispersal to the southern Andes (∼400 vs. ∼100 ky), consistent with the ancestral nature of the phaseolin seed protein and chloroplast sequences observed in the PhI group. Seed dispersal in common bean can be, therefore, described at different spatial and temporal scales, from localized, annual seed shattering to long-distance, evolutionarily rare migration.

Population genomic analysis reveals that homoploid hybrid speciation can be a lengthy process.

An increasing number of species are thought to have originated by homoploid hybrid speciation (HHS), but in only a handful of cases are details of the process known. A previous study indicated that Picea purpurea, a conifer in the Qinghai-Tibet Plateau (QTP), originated through HHS from P. likiangensis and P. wilsonii. To investigate this origin in more detail, we analysed transcriptome data for 114 individuals collected from 34 populations of the three Picea species from their core distributions in the QTP. Phylogenetic, principal component and admixture analyses of nuclear SNPs showed the species to be delimited genetically and that P. purpurea was admixed with approximately 60% of its ancestry derived from P. wilsonii and 40% from P. likiangensis. Coalescent simulations revealed the best-fitting model of origin involved formation of an intermediate hybrid lineage between P. likiangensis and P. wilsonii approximately 6 million years ago (mya), which backcrossed to P. wilsonii to form P. purpurea approximately one mya. The intermediate hybrid lineage no longer exists and is referred to as a "ghost" lineage. Our study emphasizes the power of population genomic analysis combined with coalescent analysis for reconstructing the stages involved in the origin of a homoploid hybrid species over an extended period. In contrast to other studies, we show that these stages can in some instances span a relatively long period of evolutionary time.

Peripatric speciation associated with genome expansion and female-biased sex ratios in the moss genus Ceratodon.

PREMISE OF THE STUDY: A period of allopatry is widely believed to be essential for the evolution of reproductive isolation. However, strict allopatry may be difficult to achieve in some cosmopolitan, spore-dispersed groups, like mosses. We examined the genetic and genome size diversity in Mediterranean populations of the moss Ceratodon purpureus s.l. to evaluate the role of allopatry and ploidy change in population divergence. METHODS: We sampled populations of the genus Ceratodon from mountainous areas and lowlands of the Mediterranean region, and from Western and Central Europe. We performed phylogenetic and coalescent analyses on sequences from five nuclear introns and a chloroplast locus to reconstruct their evolutionary history. We also estimated genome size using flow cytometry (employing propidium iodide) and determined the sex of samples using a sex-linked PCR marker. KEY RESULTS: Two well-differentiated clades were resolved, discriminating two homogeneous groups: the widespread C. purpureus and a local group mostly restricted to the mountains in Southern Spain. The latter also possessed a genome size 25% larger than the widespread C. purpureus, and the samples of this group consist entirely of females. We also found hybrids, and some of them had a genome size equivalent to the sum of the C. purpureus and Spanish genome, suggesting that they arose by allopolyploidy. CONCLUSIONS: These data suggest that a new species of Ceratodon arose via peripatric speciation, potentially involving a genome size change and a strong female-biased sex ratio. The new species has hybridized in the past with C. purpureus.

Assessing the potential of RAD-sequencing to resolve phylogenetic relationships within species radiations: The fly genus Chiastocheta (Diptera: Anthomyiidae) as a case study.

Determining phylogenetic relationships among recently diverged species has long been a challenge in evolutionary biology. Cytoplasmic DNA markers, which have been widely used, notably in the context of molecular barcoding, have not always proved successful in resolving such phylogenies. However, with the advent of next-generation-sequencing technologies and associated techniques of reduced genome representation, phylogenies of closely related species have been resolved at a much higher detail in the last couple of years. Here we examine the potential and limitations of one of such techniques-Restriction-site Associated DNA (RAD) sequencing, a method that produces thousands of (mostly) anonymous nuclear markers, in disentangling the phylogeny of the fly genus Chiastocheta (Diptera: Anthomyiidae). In Europe, this genus encompasses seven species of seed predators, which have been widely studied in the context of their ecological and evolutionary interactions with the plant Trollius europaeus (Ranunculaceae). So far, phylogenetic analyses using mitochondrial markers failed to resolve monophyly of most of the species from this recently diversified genus, suggesting that their taxonomy may need a revision. However, relying on a single, non-recombining marker and ignoring potential incongruences between mitochondrial and nuclear loci may provide an incomplete account of the lineage history. In this study, we applied both classical Sanger sequencing of three mtDNA regions and RAD-sequencing, for reconstructing the phylogeny of the genus. Contrasting with results based on mitochondrial markers, RAD-sequencing analyses retrieved the monophyly of all seven species, in agreement with the morphological species assignment. We found robust nuclear-based species assignment of individual samples, and low levels of estimated contemporary gene flow among them. However, despite recovering species' monophyly, interspecific relationships varied depending on the set of RAD loci considered, producing contradictory topologies. Moreover, coalescence-based phylogenetic analyses revealed low supports for most of the interspecific relationships. Our results indicate that despite the higher performance of RAD-sequencing in terms of species trees resolution compared to cytoplasmic markers, reconstructing inter-specific relationships among recently-diverged lineages may lie beyond the possibilities offered by large sets of RAD-sequencing markers in cases of strong gene tree incongruence.

A large historical refugium explains spatial patterns of genetic diversity in a Neotropical savanna tree species.

BACKGROUND AND AIMS: The relative role of Pleistocene climate changes in driving the geographic distribution and genetic diversity of South American species is not well known, especially from open biomes such as the Cerrado, the most diverse tropical savanna, encompassing high levels of endemism. Here the effects of Quaternary climatic changes on demographic history, distribution dynamics and genetic diversity of Dimorphandra mollis, an endemic tree species widely distributed in the Cerrado, were investigated. METHODS: A total of 38 populations covering most of the distribution of D. mollis were analysed using internal transcribed spacer (ITS) sequences and nuclear microsatellite variation [simple sequence repeats (SSRs)]. The framework incorporated statistical phylogeography, coalescent analyses and ecological niche modelling (ENM). KEY RESULTS: Different signatures of Quaternary climatic changes were found for ITS sequences and SSRs corresponding to different time slices. Coalescent analyses revealed large and constant effective population sizes, with high historical connectivity among the populations for ITS sequences and low effective population sizes and gene flow with recent population retraction for SSRs. ENMs indicated a slight geographical range retraction during the Last Glacial Maximum. A large historical refugium across central Brazil was predicted. Spatially explicit analyses showed a spatial cline pattern in genetic diversity related to the paleodistribution of D. mollis and to the centre of its historical refugium. CONCLUSIONS: The complex genetic patterns found in D. mollis are the result of a slight geographical range retraction during the Last Glacial Maximum followed by population expansion to the east and south from a large refugium in the central part of the Cerrado. This historical refugium is coincident with an area predicted to be climatically stable under future climate scenarios. The identified refugium should be given high priority in conservation polices to safeguard the evolutionary potential of the species under predicted future climatic changes.

The Irish potato famine pathogen Phytophthora infestans originated in central Mexico rather than the Andes.

Phytophthora infestans is a destructive plant pathogen best known for causing the disease that triggered the Irish potato famine and remains the most costly potato pathogen to manage worldwide. Identification of P. infestan's elusive center of origin is critical to understanding the mechanisms of repeated global emergence of this pathogen. There are two competing theories, placing the origin in either South America or in central Mexico, both of which are centers of diversity of Solanum host plants. To test these competing hypotheses, we conducted detailed phylogeographic and approximate Bayesian computation analyses, which are suitable approaches to unraveling complex demographic histories. Our analyses used microsatellite markers and sequences of four nuclear genes sampled from populations in the Andes, Mexico, and elsewhere. To infer the ancestral state, we included the closest known relatives Phytophthora phaseoli, Phytophthora mirabilis, and Phytophthora ipomoeae, as well as the interspecific hybrid Phytophthora andina. We did not find support for an Andean origin of P. infestans; rather, the sequence data suggest a Mexican origin. Our findings support the hypothesis that populations found in the Andes are descendants of the Mexican populations and reconcile previous findings of ancestral variation in the Andes. Although centers of origin are well documented as centers of evolution and diversity for numerous crop plants, the number of plant pathogens with a known geographic origin are limited. This work has important implications for our understanding of the coevolution of hosts and pathogens, as well as the harnessing of plant disease resistance to manage late blight.

Babesia microti from humans and ticks hold a genomic signature of strong population structure in the United States.

BACKGROUND: Babesia microti is an emerging tick-borne apicomplexan parasite with increasing geographic range and incidence in the United States. The rapid expansion of B. microti into its current distribution in the northeastern USA has been due to the range expansion of the tick vector, Ixodes scapularis, upon which the causative agent is dependent for transmission to humans. RESULTS: To reconstruct the history of B. microti in the continental USA and clarify the evolutionary origin of human strains, we used multiplexed hybrid capture of 25 B. microti isolates obtained from I. scapularis and human blood. Despite low genomic variation compared with other Apicomplexa, B. microti was strongly structured into three highly differentiated genetic clusters in the northeastern USA. Bayesian analyses of the apicoplast genomes suggest that the origin of the current diversity of B. microti in northeastern USA dates back 46 thousand years with a signature of recent population expansion in the last 1000 years. Human-derived samples belonged to two rarely intermixing clusters, raising the possibility of highly divergent infectious phenotypes in humans. CONCLUSIONS: Our results validate the multiplexed hybrid capture strategy for characterizing genome-wide diversity and relatedness of B. microti from ticks and humans. We find strong population structure in B. microti samples from the Northeast indicating potential barriers to gene flow.

Region-wide and ecotype-specific differences in demographic histories of threespine stickleback populations, estimated from whole genome sequences.

We analysed 81 whole genome sequences of threespine sticklebacks from Pacific North America, Greenland and Northern Europe, representing 16 populations. Principal component analysis of nuclear SNPs grouped populations according to geographical location, with Pacific populations being more divergent from each other relative to European and Greenlandic populations. Analysis of mitogenome sequences showed Northern European populations to represent a single phylogeographical lineage, whereas Greenlandic and particularly Pacific populations showed admixture between lineages. We estimated demographic history using a genomewide coalescence with recombination approach. The Pacific populations showed gradual population expansion starting >100 Kya, possibly reflecting persistence in cryptic refuges near the present distributional range, although we do not rule out possible influence of ancient admixture. Sharp population declines ca. 14-15 Kya were suggested to reflect founding of freshwater populations by marine ancestors. In Greenland and Northern Europe, demographic expansion started ca. 20-25 Kya coinciding with the end of the Last Glacial Maximum. In both regions, marine and freshwater populations started to show different demographic trajectories ca. 8-9 Kya, suggesting that this was the time of recolonization. In Northern Europe, this estimate was surprisingly late, but found support in subfossil evidence for presence of several freshwater fish species but not sticklebacks 12 Kya. The results demonstrate distinctly different demographic histories across geographical regions with potential consequences for adaptive processes. They also provide empirical support for previous assumptions about freshwater populations being founded independently from large, coherent marine populations, a key element in the Transporter Hypothesis invoked to explain the widespread occurrence of parallel evolution across freshwater stickleback populations.

Historical Relationships among Wild Boar Populations of the Ryukyu Archipelago and Other Eurasian regions, as Inferred from Mitochondrial Cytochrome b Gene Sequences.

The Ryukyu wild boar (Sus scrofa riukiuanus) is an endemic, morphologically defined subspecies of the Eurasian wild boar (S. scrofa) found on five islands of the Ryukyu Archipelago (a group of small islands stretching from mainland Japan to Taiwan). Two hypothetical scenarios have been proposed regarding the origin of the current Ryukyu wild boar populations: 1) natural dispersal and 2) transportation and subsequent release by prehistoric humans. To test these two hypotheses, we compared the mitochondrial cytochrome b gene sequence (1140 base pairs) in 352 individual wild boar samples that included representatives of all five insular populations of the Ryukyu wild boar and populations of other conspecific subspecies in insular East and Southeast Asia and the Eurasian Continent. A total of 68 haplotypes were recognized, of which 12 were unique to the Ryukyu wild boar populations. The results of Bayesian phylogenetic analyses supported monophyly of the five Ryukyu populations (posterior probability value of 92), confirming the validity of the subspecies as a natural group. Coalescent analysis estimated the divergence times between the Ryukyu wild boar and the other conspecific subspecies as 144-465 thousand years ago (Kya), with a 95% HPD (highest posterior density) range of 51-837 Kya, and with no significant migration. Taking the broadly accepted date of initial human migration to the Ryukyus (no earlier than 50 Kya) into consideration, our results strongly suggest that the ancestral form of the Ryukyu wild boar first entered the Ryukyu Archipelago by natural dispersal prior to the arrival of prehistoric humans.

Evolutionary and Ecological Characterization of Mayaro Virus Strains Isolated during an Outbreak, Venezuela, 2010.

In 2010, an outbreak of febrile illness with arthralgic manifestations was detected at La Estación village, Portuguesa State, Venezuela. The etiologic agent was determined to be Mayaro virus (MAYV), a reemerging South American alphavirus. A total of 77 cases was reported and 19 were confirmed as seropositive. MAYV was isolated from acute-phase serum samples from 6 symptomatic patients. We sequenced 27 complete genomes representing the full spectrum of MAYV genetic diversity, which facilitated detection of a new genotype, designated N. Phylogenetic analysis of genomic sequences indicated that etiologic strains from Venezuela belong to genotype D. Results indicate that MAYV is highly conserved genetically, showing ≈17% nucleotide divergence across all 3 genotypes and 4% among genotype D strains in the most variable genes. Coalescent analyses suggested genotypes D and L diverged ≈150 years ago and genotype diverged N ≈250 years ago. This virus commonly infects persons residing near enzootic transmission foci because of anthropogenic incursions.

Target capture and genome skimming for plant diversity studies.

Recent technological advances in long-read high-throughput sequencing and assembly methods have facilitated the generation of annotated chromosome-scale whole-genome sequence data for evolutionary studies; however, generating such data can still be difficult for many plant species. For example, obtaining high-molecular-weight DNA is typically impossible for samples in historical herbarium collections, which often have degraded DNA. The need to fast-freeze newly collected living samples to conserve high-quality DNA can be complicated when plants are only found in remote areas. Therefore, short-read reduced-genome representations, such as target capture and genome skimming, remain important for evolutionary studies. Here, we review the pros and cons of each technique for non-model plant taxa. We provide guidance related to logistics, budget, the genomic resources previously available for the target clade, and the nature of the study. Furthermore, we assess the available bioinformatic analyses, detailing best practices and pitfalls, and suggest pathways to combine newly generated data with legacy data. Finally, we explore the possible downstream analyses allowed by the type of data generated using each technique. We provide a practical guide to help researchers make the best-informed choice regarding reduced genome representation for evolutionary studies of non-model plants in cases where whole-genome sequencing remains impractical.

Variation in intraspecific demography drives localised concordance but species-wide discordance in response to past climatic change.

BACKGROUND: Understanding how species biology may facilitate resilience to climate change remains a critical factor in detecting and protecting species at risk of extinction. Many studies have focused on the role of particular ecological traits in driving species responses, but less so on demographic history and levels of standing genetic variation. Additionally, spatial variation in the interaction of demographic and adaptive factors may further complicate prediction of species responses to environmental change. We used environmental and genomic datasets to reconstruct the phylogeographic histories of two ecologically similar and largely co-distributed freshwater fishes, the southern (Nannoperca australis) and Yarra (N. obscura) pygmy perches, to assess the degree of concordance in their responses to Plio-Pleistocene climatic changes. We described contemporary genetic diversity, phylogenetic histories, demographic histories, and historical species distributions across both species, and statistically evaluated the degree of concordance in co-occurring populations. RESULTS: Marked differences in contemporary genetic diversity, historical distribution changes and historical migration were observed across the species, with a distinct lack of genetic diversity and historical range expansion suggested for N. obscura. Although several co-occurring populations within a shared climatic refugium demonstrated concordant demographic histories, idiosyncratic population size changes were found at the range edges of the more spatially restricted species. Discordant responses between species were associated with low standing genetic variation in peripheral populations. This might have hindered adaptive potential, as documented in recent demographic declines and population extinctions for the two species. CONCLUSION: Our results highlight both the role of spatial scale in the degree of concordance in species responses to climate change, and the importance of standing genetic variation in facilitating range shifts. Even when ecological traits are similar between species, long-term genetic diversity and historical population demography may lead to discordant responses to ongoing and future climate change.

A high-resolution genome of an euryhaline and eurythermal rhinogoby (Rhinogobius similis Gill 1895).

Rhinogobius similis is distributed in East and Southeast Asia. It is an amphidromous species found mostly in freshwater and sometimes brackish waters. We have obtained a high-resolution assembly of the R. similis genome using nanopore sequencing, high-throughput chromosome conformation capture (Hi-C), and transcriptomic data. The assembled genome was 890.10 Mb in size and 40.15% in GC content. Including 1373 contigs with contig N50 is 1.54 Mb, and scaffold N50 is 41.51 Mb. All of the 1373 contigs were anchored on 22 pairs of chromosomes. The BUSCO evaluation score was 93.02% indicating high quality of genome assembly. The repeat sequences accounted for 34.92% of the whole genome, with retroelements (30.13%), DNA transposons (1.64%), simple repeats (2.34%), and so forth. A total of 31,089 protein-coding genes were predicted in the genome and functionally annotated using Maker, of those genes, 26,893 (86.50%) were found in InterProScan5. There were 1910 gene families expanded in R. similis, 1171 gene families contracted and 170 gene families rapidly evolving. We have compared one rapidly change gene family (PF05970) commonly found in four species (Boleophthalmus pectinirostris, Neogobius melanostomus, Periophthalmus magnuspinnatus, and R. similis), which was found probably related to the lifespan of those species. During 400-10 Ka, the period of the Guxiang Ice Age, the population of R. similis decreased drastically, and then increased gradually following the last interglacial period. A high-resolution genome of R. similis should be useful to study taxonomy, biogeography, comparative genomics, and adaptive evolution of the most speciose freshwater goby genus, Rhinogobius.

Isolation-driven divergence: speciation in a widespread North American songbird (Aves: Certhiidae).

Lineage, or true 'species', trees may differ from gene trees because of stochastic processes in molecular evolution leading to gene-tree heterogeneity. Problems with inferring species trees because of excessive incomplete lineage sorting may be exacerbated in lineages with rapid diversification or recent divergences necessitating the use of multiple loci and individuals. Many recent multilocus studies that investigate divergence times identify lineage splitting to be more recent than single-locus studies, forcing the revision of biogeographic scenarios driving divergence. Here, we use 21 nuclear loci from regional populations to re-evaluate hypotheses identified in an mtDNA phylogeographic study of the Brown Creeper (Certhia americana), as well as identify processes driving divergence. Nuclear phylogeographic analyses identified hierarchical genetic structure, supporting a basal split at approximately 32°N latitude, splitting northern and southern populations, with mixed patterns of genealogical concordance and discordance between data sets within the major lineages. Coalescent-based analyses identify isolation, with little to no gene flow, as the primary driver of divergence between lineages. Recent isolation appears to have caused genetic bottlenecks in populations in the Sierra Madre Oriental and coastal mountain ranges of California, which may be targets for conservation concerns.

Tracing back hepatitis B virus genotype D introduction and dissemination in South Brazil.

Hepatitis B virus (HBV) infection is widespread and it is considered a major health problem in the world. HBV is classified into genotypes and subgenotypes. HBV genotype D (HBV-D) has been detected worldwide with high prevalence in some specific regions from Europe and South America. In Brazil, this genotype is very frequent in the South region and its introduction and dissemination have been associated with European immigration (mainly Italian). The present study aimed to trace back the introduction and dissemination of HBV-D in South Brazil. Fifty-two chronic hepatitis B patients from two cities with an early history of Italian immigration in South Brazil were selected for the present study. HBV-DNA was detected, quantified and a partial genomic region (S/P overlapped genes) was amplified by PCR and sequenced for the determination of HBV genotypes/subgenotypes. HBV complete genome sequences of some selected samples were further obtained. Bayesian coalescent analyses were performed to estimate the HBV-D evolutionary dynamics. Phylogenetic analysis demonstrated the occurrence of three genotypes according the tree topology: HBV-D (n = 49; 94.2%), HBV-A (n = 2; 3.9%) and HBV-G (n = 1; 1.9%). In addition, HBV-D presented three subgenotypes: HBV-D3 (n = 39; 79.6%), HBV-D2 (n = 8; 16.3%), and HBV-D1 (n = 2; 4.1%). The Bayesian coalescent analysis demonstrated that the HBV-D was introduced in the 20th century. HBV-D3 was the first to be introduced in South Brazil, probably between 1904 and 1942. HBV-D2 and HBV-D1 were introduced later; HBV-D2 between 1946 and 1953 and HBV-D1 between 1954 and 1969. HBV-D3 spread at a high rate from the 1920s to the 1980s, while HBV-D2 showed a slower growth from the 1960s to the 1990s and HBV-D1 infections demonstrated low and constant population size across time. After the 2000s, a stationary growth was detected for all these three-D subgenotypes. HBV-D showed a high prevalence in South Brazil and this is possibly associated with the first introduction and dissemination of HBV-D3 at the beginning of the 20th century.

Phylodynamics of sunflower chlorotic mottle virus, an emerging pathosystem.

Distribution and epidemiological patterns of sunflower chlorotic mottle virus (SCMoV) in sunflower (Helianthus annuus L.) growing areas in Argentina were studied from 2006 to 2017. The virus was detected exclusively in the Pampas region (Entre Ríos, Santa Fe, Córdoba, La Pampa and Buenos Aires provinces). Phylodynamic analyses performed using the coat protein gene of SCMoV isolates from sunflower and weeds dated the most recent common ancestor (MRCA) back to 1887 (HPD95% = 1572-1971), which coincides with the dates of sunflower introduction in Argentina. The MRCA was located in the south of Buenos Aires province and was associated with sunflower host (posterior probability for the ancestral host, ppah = 0.98). The Bayesian phylodynamic analyses revealed the dispersal patterns of SCMoV, suggesting a link between natural host diversity, crop displacement by human activities and virus spread.

Hepatitis E virus subtype 3f strains isolated from Japanese hepatitis patients with no history of travel to endemic areas - The origin analyzed by molecular evolution.

Hepatitis E virus subtype 3f (HEV-3f) strains are usually isolated in Europe and Thailand. Recently, HEV-3f strains were detected from six acute hepatitis E patients in Japan, none of whom had a history of travel to endemic areas. We inferred the origin and transmission route of the six HEV-3f strains. A time-scaled phylogenetic tree of the six strains with reference strains was constructed using a Bayesian statistical inference framework. The time-scaled tree indicated that the six strains independently derived from similar European strains between 2008 and 2014. The pattern suggested recent inflow of multiple HEV-3f strains from Europe to Japan. Japan imports a substantial amount of pork from European countries every year. The emergence of acute hepatitis cases caused by HEV-3f strains in Japan, in patients with no history of travel abroad, might be influenced by the increased opportunities to consume pork products imported from European countries.