Serologic and Genomic Investigation of West Nile Virus in Kosovo

The prevalence of West Nile virus (WNV) is increasing across Europe, with cases emerging in previously unaffected countries. Kosovo is situated in a WNV-endemic region where the seroepidemiological data on WNV in humans remains absent. To address this issue, we have conducted a seroepidemiological investigation of 453 randomly selected sera from a hospital in Kosovo, revealing a 1.55% anti-WNV IgG seroprevalence. Comparative and phylogeographic analyses of the WNV genomes obtained by sequencing archived samples from patients with West Nile fever indicate at least two recent and distinct introductions of WNV lineage 2 into Kosovo from neighboring countries. These findings confirm the eco-epidemiological status of WNV in southeast Europe, where long- and short-range dispersion of lineage 2 strains contributes to a wider circulation via central Europe. Our results suggest an increasing risk for WNV spreading in Kosovo, underscoring the need for an integrated national surveillance program targeting vectors and avian populations for early epidemic detection, as well as the screening of blood donors to gauge the impact of virus circulation on the human population.

The Cauchy Process on Phylogenies: A Tractable Model for Pulsed Evolution

Phylogenetic comparative methods use random processes, such as the Brownian Motion, to model the evolution of continuous traits on phylogenetic trees. Growing evidence for non-gradual evolution motivated the development of complex models, often based on Lévy processes. However, their statistical inference is computationally intensive and currently relies on approximations, high-dimensional sampling, or numerical integration. We consider here the Cauchy Process (CP), a particular pure-jump Lévy process in which the trait increment along each branch follows a centered Cauchy distribution with a dispersion proportional to its length. In this work, we derive an exact algorithm to compute both the joint probability density of the tip trait values of a phylogeny under a CP and the ancestral trait values and branch increments posterior densities in quadratic time. A simulation study shows that the CP generates patterns in comparative data that are distinct from any Gaussian process, and that restricted maximum likelihood parameter estimates and root trait reconstruction are unbiased and accurate for trees with 200 tips or less. The CP has only two parameters but is rich enough to capture complex-pulsed evolution. It can reconstruct posterior ancestral trait distributions that are multimodal, reflecting the uncertainty associated with the inference of the evolutionary history of a trait from extant taxa only. Applied on empirical datasets taken from the Evolutionary Ecology and Virology literature, the CP suggests nuanced scenarios for the body size evolution of Greater Antilles Lizards and for the geographical spread of the West Nile Virus epidemics in North America, both consistent with previous studies using more complex models. The method is efficiently implemented in C with an R interface in package cauphy, which is open source and freely available online.

Evolutionary Analysis of a Parrot Bornavirus 2 Detected in a Sulphur-Crested Cockatoo (Cacatua galerita) Suggests a South American Ancestor

Parrot bornavirus (PaBV) is an RNA virus that causes Proventricular Dilatation Disease (PDD), neurological disorders, and death in Psittaciformes. Its diversity in South America is poorly known. We examined a Cacatua galerita presenting neuropathies, PDD, and oculopathies as the main signs. We detected PaBV through reverse transcription polymerase chain reaction (RT-PCR) and partial sequencing of the nucleoprotein (N) and matrix (M) genes. Maximum likelihood and Bayesian phylogenetic inferences classified it as PaBV-2. The nucleotide identity of the sequenced strain ranged from 88.3% to 90.3% against genotype PaBV-2 and from 80.2% to 84.4% against other genotypes. Selective pressure analysis detected signs of episodic diversifying selection in both the N and M genes. No recombination events were detected. Phylodynamic analysis estimated the time to the most recent common ancestor (TMRCA) as the year 1758 for genotype PaBV-2 and the year 1049 for the Orthobornavirus alphapsittaciforme species. Substitution rates were estimated at 2.73 × 10-4 and 4.08 × 10-4 substitutions per year per site for N and M, respectively. The analysis of population dynamics showed a progressive decline in the effective population size during the last century. Timescale phylogeographic analysis revealed a potential South American ancestor as the origin of genotypes 1, 2, and 8. These results contribute to our knowledge of the evolutionary origin, diversity, and dynamics of PaBVs in South America and the world. Additionally, it highlights the importance of further studies in captive Psittaciformes and the potential impact on endangered wild birds.

Conflicting Evidence between Clinical Perception and Molecular Epidemiology: The Case of Fowl Adenovirus D

Fowl adenoviruses (FAdVs, species FAdV-A/-E) are responsible for several clinical syndromes reported with increasing frequency in poultry farms in the last decades. In the present study, a phylodynamic analysis was performed on a group of FAdV-D Hexon sequences with adequate available metadata. The obtained results demonstrated the long-term circulation of this species, at least several decades before the first identification of the disease. After a period of progressive increase, the viral population showed a high-level circulation from approximately the 1960s to the beginning of the new millennium, mirroring the expansion of intensive poultry production and animal trade. At the same time, strain migration occurred mainly from Europe to other continents, although other among-continent connections were estimated. Thereafter, the viral population declined progressively, likely due to the improved control measures, potentially including the development and application of FAdV vaccines. An increase in the viral evolutionary rate featured this phase. A role of vaccine-induced immunity in shaping viral evolution could thus be hypothesized. Accordingly, several sites of the Hexon, especially those targeted by the host response were proven under a significant pervasive or episodic diversifying selection. The present study results demonstrate the role of intensive poultry production and market globalization in the rise of FAdV. The applied control strategies, on the other hand, were effective in limiting viral circulation and shaping its evolution.

Analysis of Genetic Diversity in Adzuki Beans (Vigna angularis): Insights into Environmental Adaptation and Early Breeding Strategies for Yield Improvement

Adzuki beans are widely cultivated in East Asia and are one of the earliest domesticated crops. In order to gain a deeper understanding of the genetic diversity and domestication history of adzuki beans, we conducted Genotyping by Sequencing (GBS) analysis on 366 landraces originating from Korea, China, and Japan, resulting in 6586 single-nucleotide polymorphisms (SNPs). Population structure analysis divided these 366 landraces into three subpopulations. These three subpopulations exhibited distinctive distributions, suggesting that they underwent extended domestication processes in their respective regions of origin. Phenotypic variance analysis of the three subpopulations indicated that the Korean-domesticated subpopulation exhibited significantly higher 100-seed weights, the Japanese-domesticated subpopulation showed significantly higher numbers of grains per pod, and the Chinese-domesticated subpopulation displayed significantly higher numbers of pods per plant. We speculate that these differences in yield-related traits may be attributed to varying emphases placed by early breeders in these regions on the selection of traits related to yield. A large number of genes related to biotic/abiotic stress resistance and defense were found in most quantitative trait locus (QTL) for yield-related traits using genome-wide association studies (GWAS). Genomic sliding window analysis of Tajima's D and a genetic differentiation coefficient (Fst) revealed distinct domestication selection signatures and genotype variations on these QTLs within each subpopulation. These findings indicate that each subpopulation would have been subjected to varied biotic/abiotic stress events in different origins, of which these stress events have caused balancing selection differences in the QTL of each subpopulation. In these balancing selections, plants tend to select genotypes with strong resistance under biotic/abiotic stress, but reduce the frequency of high-yield genotypes to varying degrees. These biotic/abiotic stressors impact crop yield and may even lead to selection purging, resulting in the loss of several high-yielding genotypes among landraces. However, this also fuels the flow of crop germplasms. Overall, balancing selection appears to have a more significant impact on the three yield-related traits compared to breeder-driven domestication selection. These findings are crucial for understanding the impact of domestication selection history on landraces and yield-related traits, aiding in the improvement of adzuki bean varieties.

Molecular insights on the historical dispersion of Piezodorus guildinii (Hemiptera: Pentatomidae) in Brazil

Piezodorus guildinii (Westwood, 1837) (Hemiptera: Pentatomidae) is an important arthropod pest of soybean (Glycine max (L.) Merr.) throughout American continents. However, the historical events associated with its dispersion are poorly understood. In this study, we employed a phylogeographic approach to investigate the origin and demographic history of P. guildinii in Brazil. We analyzed the cytochrome c oxidase subunit I and Cytb gene sequences of P. guildinii individuals collected in Brazil's 5 soybean production macro-regions and cross-referenced this information with sequences available in public databases. Our findings support an older Caribbean basin establishment for the current genealogical strains of P. guildinii, with subsequent dispersion to Brazil around 0.97 Mya. No secondary dispersion of this species from the Caribbean region to soybean areas in Brazil was identified. The Brazilian populations of P. guildinii are genetically structured across the country's soybean macro-regions and show strong signals of continuous demographic and spatial expansion in Brazil, which may be accelerated by the soybean cropping landscape in the country. The populations from the northern region (MR5) are older than the Central and South populations. The signs of demographic expansion indicate that P. guildinii populations are increasing their effective size in soybean regions, which could reflect its importance as a soybean pest in the coming years.

Molecular Epidemiology of Respiratory Syncytial Virus during 2019-2022 and Surviving Genotypes after the COVID-19 Pandemic in Japan

To evaluate the changes in respiratory syncytial virus (RSV) collected between 2019 and 2022, we analyzed RSV-A and RSV-B strains from various prefectures in Japan before and after the COVID-19 pandemic. RT-PCR-positive samples collected from children with rapid test positivity at outpatient clinics in 11 prefectures in Japan were sequenced for the ectodomain of the G gene to determine the genotype. Time-aware phylogeographic analyses were performed using the second hypervariable region (HVR) of the G gene from 2012 to 2022. Of 967 samples, 739 (76.4%) were found to be RSV-positive using RT-PCR. RSV peaked in September 2019 but was not detected in 2020, except in Okinawa. Nationwide epidemics occurred with peaks in July 2021 and 2022. The genotype remained the same, ON1 for RSV-A and BA9 for RSV-B during 2019-2022. Phylogeographic analysis of HVR revealed that at least seven clusters of RSV-A had circulated previously but decreased to two clusters after the pandemic, whereas RSV-B had a single monophyletic cluster over the 10 years. Both RSV-A and RSV-B were transferred from Okinawa into other prefectures after the pandemic. The RSV epidemic was suppressed due to pandemic restrictions; however, pre-pandemic genotypes spread nationwide after the pandemic.

Genomic landscape of introgression from the ghost lineage in a gobiid fish uncovers the generality of forces shaping hybrid genomes

Extinct lineages can leave legacies in the genomes of extant lineages through ancient introgressive hybridization. The patterns of genomic survival of these extinct lineages provide insight into the role of extinct lineages in current biodiversity. However, our understanding on the genomic landscape of introgression from extinct lineages remains limited due to challenges associated with locating the traces of unsampled 'ghost' extinct lineages without ancient genomes. Herein, we conducted population genomic analyses on the East China Sea (ECS) lineage of Chaenogobius annularis, which was suspected to have originated from ghost introgression, with the aim of elucidating its genomic origins and characterizing its landscape of introgression. By combining phylogeographic analysis and demographic modelling, we demonstrated that the ECS lineage originated from ancient hybridization with an extinct ghost lineage. Forward simulations based on the estimated demography indicated that the statistic γ of the HyDe analysis can be used to distinguish the differences in local introgression rates in our data. Consistent with introgression between extant organisms, we found reduced introgression from extinct lineage in regions with low recombination rates and with functional importance, thereby suggesting a role of linked selection that has eliminated the extinct lineage in shaping the hybrid genome. Moreover, we identified enrichment of repetitive elements in regions associated with ghost introgression, which was hitherto little known but was also observed in the re-analysis of published data on introgression between extant organisms. Overall, our findings underscore the unexpected similarities in the characteristics of introgression landscapes across different taxa, even in cases of ghost introgression.

Phylogeny and spatial distribution of Japanese encephalitis virus vector species in Cambodia

In Southeast Asia, despite the use of Japanese encephalitis vaccines and vaccination coverage, Japanese encephalitis (JE) transmission is still a major public health issue. The main vectors of this virus are mosquitoes from the genus Culex, which diversity and density are important in Southeast Asia. The main vector species of Japanese encephalitis virus (JEV) in Cambodia belong to the Vishnui subgroup. However, their morphological identification solely based on the adult stage remains challenging, making their segregation and detection difficult. In order to identify and describe the distribution of the three main JEV vector species in Cambodia, namely Culex vishnui, Cx. pseudovishnui and Cx. tritaeniorhynchus, mosquito samplings were carried out throughout the country in different environments. Phylogenetic analysis of the cytochrome c oxidase subunit I (coI) gene using maximum-likelihood tree with ultrafast bootstrap and phylogeographic analysis were performed. The three main Culex species are phylogenetically separated, and represent two distinct clades, one with Cx. tritaeniorhynchus and the second with Cx. vishnui and Cx. pseudovishnui, the latter appearing as a subgroup of Cx. vishnui. The phylogeographic analysis shows a distribution of the Vishnui subgroup on the entire Cambodian territory with an overlapped distribution areas leading to a sympatric distribution of these species. The three JEV vector species are geographically well-defined with a strong presence of Cx. pseudovishnui in the forest. Combined with the presence of Cx. tritaeniorhynchus and Cx. vishnui in rural, peri-urban, and urban areas, the presence of JEV-competent vectors is widespread in Cambodia.

Approximate Bayesian computation and ecological niche models elucidate the demographic history and current fragmented population distribution of a Korean endemic shrub

Climatic fluctuations and geological events since the LGM are believed to have significantly impacted the population size, distribution, and mobility of many species that we observe today. In this paper, we determined the processes driving the phylogeographic structure of the Korean endemic white forsythia by combining the use of genome-wide SNPs and predicting paleoclimatic habitats during the LGM (21 kya), Early Holocene (10 kya), Mid-Holocene (6 kya), and Late Holocene (3 kya). Using a maximum of 1897 SNPs retrieved from 124 samples across nine wild populations, five environmental predictors, and the species' natural occurrence records, we aimed to infer the species' demographic history and reconstruct its possible paleodistributions with the use of approximate Bayesian computation and ecological niche models, respectively. Under this integrated framework, we found strong evidence for patterns of range shift and expansion, and population divergence events from the onset of the Holocene, resulting in the formation of its five distinct genetic units. The most highly supported model inferred that after the split of an ancestral population into the southern group and a larger central metapopulation lineage, the latter gave rise to the eastern and northern clusters, before finally dividing into two sub-central groups. While the use of molecular data allowed us to identify and refine the (phylo)genetic relationships of the species' lineages and populations, the use of ecological data helped us infer a past LGM refugium and the directions of post-glacial range dynamics. The time frames of these demographic events were shown to be congruent with climatic and geological events that affected the central Korean Peninsula during these periods. These findings gave us a better understanding of the consequences of past spatiotemporal factors that may have resulted in the current fragmented population distribution of this endangered plant.

Nuclear and chloroplast DNA phylogeography reveals high genetic diversity and postglacial range expansion in Quercus mexicana

PREMISE

Phylogeographical studies are fundamental for understanding factors that influence the spatial distribution of genetic lineages within species. Population expansions and contractions, distribution shifts, and climate changes are among the most important factors shaping the genetic compositions of populations.

METHODS

We investigated the phylogeography of an endemic oak, Quercus mexicana (Fagaceae), which has a restricted distribution in northeastern Mexico along the Sierra Madre Oriental and adjacent areas. Nuclear and chloroplast DNA microsatellite markers were used to describe the genetic diversity and structure of 39 populations of Q. mexicana along its entire distribution area. We tested whether population expansion or contraction events influenced the genetic diversity and structure of the species. We also modeled the historical distributional range of Q. mexicana (for the Mid Holocene, the Last Glacial Maximum, and the Last Interglacial) to estimate the extent to which climate fluctuations have impacted the distribution of this oak species.

RESULTS

Our results revealed high genetic diversity and low genetic structure in Q. mexicana populations. Ecological niche models suggested historical fluctuations in the distributional range of Q. mexicana. Historical range changes, gene flow, and physical barriers seem to have played an important role in shaping the phylogeographic structure of Q. mexicana.

CONCLUSIONS

Our study indicates that the genetic structure of Q. mexicana may have been the result of responses of oak trees not only to heterogeneous environments present in the Sierra Madre Oriental and adjacent areas, but also to elevational and latitudinal shifts in response to climate changes in the past.

Is there a massive glacial-Holocene flora continuity in Central Europe?

The prevailing paradigm about the Quaternary ecological and evolutionary history of Central European ecosystems is that they were repeatedly impoverished by regional extinctions of most species during the glacial periods, followed by massive recolonizations from southern and eastern refugia during interglacial periods. Recent literature partially contradicts this view and provides evidence to re-evaluate this Postglacial Recolonization Hypothesis and develop an alternative one. We examined the long-term history of the flora of the Carpathian (Pannonian) Basin by synthesising recent advances in ecological, phylogeographical, palaeoecological and palaeoclimatological research, and analysing the cold tolerance of the native flora of a test area (Hungary, the central part of the Carpathian Basin). We found that (1) many species have likely occurred there continuously since before the Last Glacial Maximum (LGM); (2) most of the present-day native flora (1404 species, about 80%) can occur in climates as cold as or colder than the LGM (mean annual temperature ≤+3.5°C); and (3) grasslands and forests can be species-rich under an LGM-like cold climate. These arguments support an alternative hypothesis, which we call the Flora Continuity Hypothesis. It states that long-term continuity of much of the flora in the Carpathian Basin is more plausible than regional extinctions during the LGM followed by massive postglacial recolonizations. The long-term continuity of the region's flora may have fundamental implications not only for understanding local biogeography and ecology (e.g. the temporal scale of processes), but also for conservation strategies focusing on protecting ancient species-rich ecosystems and local gene pools.

Genomic epidemiology sheds light on the emergence and spread of Mycobacterium bovis Eu2 Clonal Complex in Portugal

ABSTRACTAnimal tuberculosis (TB) remains a serious concern for animal and human health. Mycobacterium bovis circulates in multi-host systems, dominated by the European 2 clonal complex (Eu2) in Iberia. In this work, we use genomic epidemiology to infer the emergence, spread, and spatiotemporal patterns of Eu2 in the official epidemiological risk area of animal TB in Portugal. Phylogenetic analysis of 144 M. bovis whole-genome sequences from cattle, wild boar, and red deer, representing the 2002-2021 period, distinguished three Eu2 clades that evolved independently. The major Eu2 clade underwent phylodynamic inferences to estimate the time and location of outbreaks, host transitions, and spatial diffusion as well. The origin of this Eu2 clade was attributed to the red deer population in the Castelo Branco district, near the border with Spain. Most host transitions were intraspecific (80%), while interspecific transmissions between wildlife species (wild boar-red deer), and between wild boar and cattle, were highly supported. Phylogeographic reconstruction evidenced that most transitions (82%) occur within municipalities, highlighting local transmission corridors.Our study indicates that M. bovis continues to spread at the cattle-wildlife interface within the animal TB hotspot area, possibly driven by the foraging behaviour of wild boar near agricultural lands. Red deer seems to be an important driver of TB within wildlife hosts, while the wild boar links the multi-host wildlife community and livestock. This work highlights the value of combining genomic epidemiology with phylodynamic inference to resolve host jumps and spatial patterns of M. bovis, providing real-time clues about points of intervention.

Genetic diversity and microevolution in clinical Cryptococcus isolates from Cameroon

Cryptococcal meningitis is the second most common cause of death in people living with HIV/AIDS, yet we have a limited understanding of how cryptococcal isolates change over the course of infection. Cryptococcal infections are environmentally acquired, and the genetic diversity of these infecting isolates can also be geographically linked. Here, we employ whole genome sequences for 372 clinical Cryptococcus isolates from 341 patients with HIV-associated cryptococcal meningitis obtained via a large clinical trial, across both Malawi and Cameroon, to enable population genetic comparisons of isolates between countries. We see that isolates from Cameroon are highly clonal, when compared to those from Malawi, with differential rates of disruptive variants in genes with roles in DNA binding and energy use. For a subset of patients (22) from Cameroon, we leverage longitudinal sampling, with samples taken at days 7 and 14 post-enrollment, to interrogate the genetic changes that arise over the course of infection, and the genetic diversity of isolates within patients. We see disruptive variants arising over the course of infection in several genes, including the phagocytosis-regulating transcription factor GAT204. In addition, in 13% of patients sampled longitudinally, we see evidence for mixed infections. This approach identifies geographically linked genetic variation, signatures of microevolution, and evidence for mixed infections across a clinical cohort of patients affected by cryptococcal meningitis in Central Africa.

Complete genome reconstruction of the global and European regional dispersal history of the lumpy skin disease virus

IMPORTANCE

Lumpy skin disease virus (LSDV) has a complex epidemiology involving multiple strains, recombination, and vaccination. Its DNA genome provides limited genetic variation to trace outbreaks in space and time. Sequencing of LSDV whole genomes has also been patchy at global and regional scales. Here, we provide the first fine-grained whole genome sequence sampling of a constrained LSDV outbreak (southeastern Europe, 2015-2017), which we analyze along with global publicly available genomes. We formally evaluate the past occurrence of recombination events as well as the temporal signal that is required for calibrating molecular clock models and subsequently conduct a time-calibrated spatially explicit phylogeographic reconstruction. Our study further illustrates the importance of accounting for recombination events before reconstructing global and regional dynamics of DNA viruses. More LSDV whole genomes from endemic areas are needed to obtain a comprehensive understanding of global LSDV dispersal dynamics.

Speciation in Coastal Basins Driven by Staggered Headwater Captures: Dispersal of a Species Complex, Leporinus bahiensis, as Revealed by Genome-wide SNP Data

Past sea level changes and geological instability along watershed boundaries have largely influenced fish distribution across coastal basins, either by dispersal via palaeodrainages now submerged or by headwater captures, respectively. Accordingly, the South American Atlantic coast encompasses several small and isolated drainages that share a similar species composition, representing a suitable model to infer historical processes. Leporinus bahiensis is a freshwater fish species widespread along adjacent coastal basins over narrow continental shelf with no evidence of palaeodrainage connections at low sea level periods. Therefore, this study aimed to reconstruct its evolutionary history to infer the role of headwater captures in the dispersal process. To accomplish this, we employed molecular-level phylogenetic and population structure analyses based on Sanger sequences (5 genes) and genome-wide SNP data. Phylogenetic trees based on Sanger data were inconclusive, but SNPs data did support the monophyletic status of L. bahiensis. Both COI and SNP data revealed structured populations according to each hydrographic basin. Species delimitation analyses revealed from 3 (COI) to 5 (multilocus approach) MOTUs, corresponding to the sampled basins. An intricate biogeographic scenario was inferred and supported by Approximate Bayesian Computation (ABC) analysis. Specifically, a staggered pattern was revealed and characterized by sequential headwater captures from basins adjacent to upland drainages into small coastal basins at different periods. These headwater captures resulted in dispersal throughout contiguous coastal basins, followed by deep genetic divergence among lineages. To decipher such recent divergences, as herein represented by L. bahiensis populations, we used genome-wide SNPs data. Indeed, the combined use of genome-wide SNPs data and ABC method allowed us to reconstruct the evolutionary history and speciation of L. bahiensis. This framework might be useful in disentangling the diversification process in other neotropical fishes subject to a reticulate geological history.

Genetic signatures of lineage fusion closely resemble population decline

Accurate interpretation of the genetic signatures of past demographic events is crucial for reconstructing evolutionary history. Lineage fusion (complete merging, resulting in a single panmictic population) is a special case of secondary contact that is seldom considered. Here, the circumstances under which lineage fusion can be distinguished from population size constancy, growth, bottleneck, and decline were investigated. Multi-locus haplotype data were simulated under models of lineage fusion with different divergence versus sampling lag times (D:L ratios). These pseudo-observed datasets also differed in their allocation of a fixed amount of sequencing resources (number of sampled alleles, haplotype length, number of loci). Distinguishability of lineage fusion versus each of 10 untrue non-fusion scenarios was quantified based on six summary statistics (neutrality tests). Some datasets were also analyzed using extended Bayesian skyline plots. Results showed that signatures of lineage fusion very closely resemble those of decline-high distinguishability was generally limited to the most favorable scenario (D:L = 9), using the most sensitive summary statistics (F S and Z nS), coupled with the optimal sequencing resource allocation (maximizing number of loci). Also, extended Bayesian skyline plots often erroneously inferred population decline. Awareness of the potential for lineage fusion to carry the hallmarks of population decline is critical.

Evolution of horn length and lifting strength in the Japanese rhinoceros beetle Trypoxylus dichotomus

What limits the size of nature's most extreme structures? For weapons like beetle horns, one possibility is a tradeoff associated with mechanical levers: as the output arm of the lever system-the beetle horn-gets longer, it also gets weaker. This "paradox of the weakening combatant" could offset reproductive advantages of additional increases in weapon size. However, in contemporary populations of most heavily weaponed species, males with the longest weapons also tend to be the strongest, presumably because selection drove the evolution of compensatory changes to these lever systems that ameliorated the force reductions of increased weapon size. Therefore, we test for biomechanical limits by reconstructing the stages of weapon evolution, exploring whether initial increases in weapon length first led to reductions in weapon force generation that were later ameliorated through the evolution of mechanisms of mechanical compensation. We describe phylogeographic relationships among populations of a rhinoceros beetle and show that the "pitchfork" shaped head horn likely increased in length independently in the northern and southern radiations of beetles. Both increases in horn length were associated with dramatic reductions to horn lifting strength-compelling evidence for the paradox of the weakening combatant-and these initial reductions to horn strength were later ameliorated in some populations through reductions to horn length or through increases in head height (the input arm for the horn lever system). Our results reveal an exciting geographic mosaic of weapon size, weapon force, and mechanical compensation, shedding light on larger questions pertaining to the evolution of extreme structures.

Molecular Epidemiology of Rotavirus A in Calves: Evolutionary Analysis of a Bovine G8P[11] Strain and Spatio-Temporal Dynamics of G6 Lineages in the Americas

Rotavirus A (RVA) causes diarrhea in calves and frequently possesses the G6 and P[5]/P[11] genotypes, whereas G8 is less common. We aimed to compare RVA infections and G/P genotypes in beef and dairy calves from major livestock regions of Argentina, elucidate the evolutionary origin of a G8 strain and analyze the G8 lineages, infer the phylogenetic relationship of RVA field strains, and investigate the evolution and spatio-temporal dynamics of the main G6 lineages in American countries. Fecal samples (n = 422) from diarrheic (beef, 104; dairy, 137) and non-diarrheic (beef, 78; dairy, 103) calves were analyzed by ELISA and semi-nested multiplex RT-PCR. Sequencing, phylogenetic, phylodynamic, and phylogeographic analyses were performed. RVA infections were more frequent in beef (22.0%) than in dairy (14.2%) calves. Prevalent genotypes and G6 lineages were G6(IV)P[5] in beef (90.9%) and G6(III)P[11] (41.2%) or mixed genotypes (23.5%) in dairy calves. The only G8 strain was phylogenetically related to bovine and artiodactyl bovine-like strains. Re-analyses inside the G8 genotype identified G8(I) to G8(VIII) lineages. Of all G6 strains characterized, the G6(IV)P[5](I) strains from "Cuenca del Salado" (Argentina) and Uruguay clustered together. According to farm location, a clustering pattern for G6(IV)P[5] strains of beef farms was observed. Both G6 lineage strains together revealed an evolutionary rate of 1.24 × 10-3 substitutions/site/year, and the time to the most recent common ancestor was dated in 1853. The most probable ancestral locations were Argentina in 1981 for G6(III) strains and the USA in 1940 for G6(IV) strains. The highest migration rates for both G6 lineages together were from Argentina to Brazil and Uruguay. Altogether, the epidemiology, genetic diversity, and phylogeny of RVA in calves can differ according to the production system and farm location. We provide novel knowledge about the evolutionary origin of a bovine G8P[11] strain. Finally, bovine G6 strains from American countries would have originated in the USA nearly a century before its first description.

The Temporal and Geographical Dynamics of Potato Virus Y Diversity in Russia

Potato virus Y, an important viral pathogen of potato, has several genetic variants and geographic distributions which could be affected by environmental factors, aphid vectors, and reservoir plants. PVY is transmitted to virus-free potato plants by aphids and passed on to the next vegetative generations through tubers, but the effects of tuber transmission in PVY is largely unknown. By using high-throughput sequencing, we investigated PVY populations transmitted to potato plants by aphids in different climate zones of Russia, namely the Moscow and Astrakhan regions. We analyzed sprouts from the tubers produced by field-infected plants to investigate the impact of tuber transmission on PVY genetics. We found a significantly higher diversity of PVY isolates in the Astrakhan region, where winters are shorter and milder and summers are warmer compared to the Moscow region. While five PVY types, NTNa, NTNb, N:O, N-Wi, and SYR-I, were present in both regions, SYRI-II, SYRI-III, and 261-4 were only found in the Astrakhan region. All these recombinants were composed of the genome sections derived from PVY types O and N, but no full-length sequences of such types were present. The composition of the PVY variants in the tuber sprouts was not always the same as in their parental plants, suggesting that tuber transmission impacts PVY genetics.

Phylodynamic approaches to studying avian influenza virus

Avian influenza viruses can cause severe disease in domestic and wild birds and are a pandemic threat. Phylodynamics is the study of how epidemiological, evolutionary, and immunological processes can interact to shape viral phylogenies. This review summarizes how phylodynamic methods have and could contribute to the study of avian influenza viruses. Specifically, we assess how phylodynamics can be used to examine viral spread within and between wild or domestic bird populations at various geographical scales, identify factors associated with virus dispersal, and determine the order and timing of virus lineage movement between geographic regions or poultry production systems. We discuss factors that can complicate the interpretation of phylodynamic results and identify how future methodological developments could contribute to improved control of the virus.

Population genetic structure of a Chihuahuan Desert endemic mammal, the desert pocket gopher, Geomys arenarius

The biogeographic history of the Chihuahuan Desert is complex, driven by numerous physiographic events and climatic changes. This dynamic history would have influenced the flora and fauna of the region including the desert pocket gopher, Geomys arenarius, a subterranean rodent endemic to the northern Chihuahuan Desert. G. arenarius is restricted to sandy soils and are considered to have a disjunct distribution. Two subspecies are recognized: G. a. arenarius and G. a. brevirostris. We used multilocus nuclear (amplified fragment length polymorphisms) and mitochondrial DNA (ND2) sequence data to uncover patterns of genetic diversity within and among populations of G. arenarius. We evaluated correspondence of genetic patterns to traditionally accepted subspecies boundaries, mapped the distribution of potentially suitable soils to identify barriers or corridors to dispersal and to guide future survey efforts, provided evidence that could be used to recognize distinct population segments, and quantified genetic diversity within populations. Both datasets were largely concordant and demonstrated hierarchical patterns of genetic divergence. The greatest divergence was consistent with the two recognized subspecies. Mapping of potentially habitable soils revealed likely barriers to dispersal contributing to the allopatric pattern of geographic distribution and areas, which may be occupied by G. arenarius but not yet documented. Because G. arenarius is restricted to soils with high sand content, and these habitable soils are disjunct within the region occupied by this species, historical factors that impacted soil deposition and deflation likely contributed to the observed patterns of genetic divergence. Genetic diversity was higher within populations of the southern subspecies (G. a. arenarius) compared to G. a. brevirostris. This may be due to a greater availability of continuous suitable soils within the range of G. a. arenarius or higher density due to greater food availability (currently or historically)-both of which could allow for a higher effective population size.

Loss of species and genetic diversity during colonization: Insights from acanthocephalan parasites in northern European seals

Studies on host-parasite systems that have experienced distributional shifts, range fragmentation, and population declines in the past can provide information regarding how parasite community richness and genetic diversity will change as a result of anthropogenic environmental changes in the future. Here, we studied how sequential postglacial colonization, shifts in habitat, and reduced host population sizes have influenced species richness and genetic diversity of Corynosoma (Acanthocephala: Polymorphidae) parasites in northern European marine, brackish, and freshwater seal populations. We collected Corynosoma population samples from Arctic, Baltic, Ladoga, and Saimaa ringed seal subspecies and Baltic gray seals, and then applied COI barcoding and triple-enzyme restriction-site associated DNA (3RAD) sequencing to delimit species, clarify their distributions and community structures, and elucidate patterns of intraspecific gene flow and genetic diversity. Our results showed that Corynosoma species diversity reflected host colonization histories and population sizes, with four species being present in the Arctic, three in the Baltic Sea, two in Lake Ladoga, and only one in Lake Saimaa. We found statistically significant population-genetic differentiation within all three Corynosoma species that occur in more than one seal (sub)species. Genetic diversity tended to be high in Corynosoma populations originating from Arctic ringed seals and low in the landlocked populations. Our results indicate that acanthocephalan communities in landlocked seal populations are impoverished with respect to both species and intraspecific genetic diversity. Interestingly, the loss of genetic diversity within Corynosoma species seems to have been less drastic than in their seal hosts, possibly due to their large local effective population sizes resulting from high infection intensities and effective intra-host population mixing. Our study highlights the utility of genomic methods in investigations of community composition and genetic diversity of understudied parasites.

Two Lineages of Oribatid Mites Morphologically Correspond to the Circumpolar Species Ameronothrus nigrofemoratus (Acari, Oribatida) but Differ Genetically as Distinct Species Are Revealed on the Kolguev Island

Data reported from Northern Canada were until recently the only available data on the genetic characteristics of the oribatid mite Ameronothrus nigrofemoratus, which has a circumpolar distribution on the coasts of Arctic seas. A partial cytochrome oxidase I gene (COI) mtDNA sequence was examined in mites morphologically assigned to this species from the Kolguev Island. Two highly divergent phylogenetic lineages of A. nigrofemoratus (7% divergence) were revealed, neither of which was found on the Canadian coast. Four COI amino acid substitutions distinguished one of the lineages from North American A. nigrofemoratus, corresponding to the degree of difference between A. nigrofemoratus and its sister species A. lineatus.

Evolutionary dynamics of Euphorbia carniolica suggest a complex Plio-Pleistocene history of understorey species of deciduous forest in southeastern Europe

Deciduous forests form the dominant natural vegetation of Europe today, but were restricted to small refugia during Pleistocene cold stages, implying an evolutionary past shaped by recurrent range contractions and expansions. Cold-stage forest refugia were probably widespread in southern and central Europe, with the northwestern Balkan Peninsula being of particular importance. However, the actual number and location of deciduous forest refugia, as well as the connections between them, remain disputed. Here, we address the evolutionary dynamics of the deciduous forest understorey species Euphorbia carniolica as a proxy for past forest dynamics. To do so, we obtained genomic and morphometric data from populations representing the species' entire range, investigated phylogenetic position and intraspecific genetic variation, tested explicit demographic scenarios and applied species distribution models. Our data support two disjoint groups linked to separate refugia on the northwestern and central Balkan Peninsula. We find that genetic differentiation between groups started in the early Pleistocene via vicariance, suggesting a larger distribution in the past. Both refugia acted as sources for founder events to the southeastern Alps and the Carpathians; the latter were likely colonised before the last cold stage. In line with traditional views on the pre-Pleistocene origin of many southeastern European deciduous forest species, the origin of E. carniolica was dated to the late Pliocene. The fact that E. carniolica evolved at a time when a period of continuous forestation was ending in much of Eurasia provides an interesting biogeographical perspective on the past links between Eurasian deciduous forests and their biota.