Equus roundworms (Parascaris spp.) are undergoing divergence due to natural and anthropogenic factors

The evolution of parasites is often directly affected by the host's environment or behavior. Studies on the evolution of the same parasites in different hosts are extremely attractive and highly relevant to our understanding of divergence and speciation. Here we presented the first molecular evidence of divergence of Equus roundworms in different hosts (horses, zebras and donkeys). At the genetic level, Equus roundworms were mainly separated into two clades (Horse-derived and Zebra & Donkey-derived). This divergence began at 600–1500 years ago, which interestingly coincided with the domestication history of horses. We found that most of the key enzymes related to glycolysis were under strong positive selection in zebra & donkey-derived roundworms, indicating that the evolution of the metabolic level was one of the main reasons for the divergence. In addition, we conducted a selective scan of resistance-related genes and found that the three populations were under different degrees of selection. This prompted us to pay attention to the possible impact of drugs on divergence, not just the drug resistance. This work supports that divergence or speciation is a continuous and dynamic process, and continuous monitoring of environmental factors is conducive to further understanding the adaptive evolution of roundworms.

Large genetic diversity and strong positive selection in F-box and GPCR genes among the wild isolates of Caenorhabditis elegans

The F-box and chemosensory GPCR (csGPCR) gene families are greatly expanded in nematodes, including the model organism Caenorhabditis elegans, compared to insects and vertebrates. However, the intraspecific evolution of these two gene families in nematodes remain unexamined. In this study, we analyzed the genomic sequences of 330 recently sequenced wild isolates of C. elegans using a range of population genetics approaches. We found that F-box and csGPCR genes, especially the Srw family csGPCRs, showed much more diversity than other gene families. Population structure analysis and phylogenetic analysis divided the wild strains into eight non-Hawaiian and three Hawaiian subpopulations. Some Hawaiian strains appeared to be more ancestral than all other strains. F-box and csGPCR genes maintained a great amount of the ancestral variants in the Hawaiian subpopulation and their divergence among the non-Hawaiian subpopulations contributed significantly to population structure. F-box genes are mostly located at the chromosomal arms and high recombination rate correlates with their large polymorphism. Moreover, using both neutrality tests and Extended Haplotype Homozygosity analysis, we identified signatures of strong positive selection in the F-box and csGPCR genes among the wild isolates, especially in the non-Hawaiian population. Accumulation of high-frequency derived alleles in these genes was found in non-Hawaiian population, leading to divergence from the ancestral genotype. In summary, we found that F-box and csGPCR genes harbour a large pool of natural variants, which may be subjected to positive selection. These variants are mostly mapped to the substrate-recognition domains of F-box proteins and the extracellular and intracellular regions of csGPCRs, possibly resulting in advantages during adaptation by affecting protein degradation and the sensing of environmental cues, respectively.

Population Genomics Insights into the First Wave of COVID-19

Full-genome-sequence computational analyses of the SARS-coronavirus (CoV)-2 genomes allow us to understand the evolutionary events and adaptability mechanisms. We used population genetics analyses on human SARS-CoV-2 genomes available on 2 April 2020 to infer the mutation rate and plausible recombination events between the Betacoronavirus genomes in nonhuman hosts that may have contributed to the evolution of SARS-CoV-2. Furthermore, we localized the targets of recent and strong, positive selection during the first pandemic wave. The genomic regions that appear to be under positive selection are largely co-localized with regions in which recombination from nonhuman hosts took place. Our results suggest that the pangolin coronavirus genome may have contributed to the SARS-CoV-2 genome by recombination with the bat coronavirus genome. However, we find evidence for additional recombination events that involve coronavirus genomes from other hosts, i.e., hedgehogs and sparrows. We further infer that recombination may have recently occurred within human hosts. Finally, we estimate the parameters of a demographic scenario involving an exponential growth of the size of the SARS-CoV-2 populations that have infected European, Asian, and Northern American cohorts, and we demonstrate that a rapid exponential growth in population size from the first wave can support the observed polymorphism patterns in SARS-CoV-2 genomes.

Phylogenomics Analysis of SARS-CoV2 Genomes Reveals Distinct Selection Pressure on Different Viral Strains

We are witnessing a tremendous outbreak of a novel coronavirus (SARS-CoV2) across the globe. Upon exposure to different population and changing environment, the viral strain might experience different mutational bias that leads to genetic diversity among the viral population. Also, the diversification can be influenced by distinct selection pressure on different viral genomes. We have carried out a comparative genomic analysis of 82 SARS-CoV2 genomes. We have evaluated their evolutionary divergence, substitution pattern, and rates. Viral genomes under distinct selection pressure have been identified. Sites that experience strong selection pressure also have been identified. Our result shows that the translational preference of a few codons is strongly correlated with the mutational bias imposed by genome compositional constraint and influenced by natural selection. Few genomes are evolving with a higher mutational rate with a distinct signature of nucleotide substitution in comparison to others. Four viral strains are under the effect of purifying selection, while nine SARS-CoV2 genomes are under strong positive selection bias. Site analysis indicates a strong positive selection pressure on two codon positions at 3606th and 8439th positions. Our study elucidates adaptation of few SARS-CoV2 viral strain during the outbreak shaping by natural selection and genomic compositional constraints.

Innovation of heterochromatin functions drives rapid evolution of essential ZAD-ZNF genes in Drosophila

Contrary to dogma, evolutionarily young and dynamic genes can encode essential functions. We find that evolutionarily dynamic ZAD-ZNF genes, which encode the most abundant class of insect transcription factors, are more likely to encode essential functions in Drosophila melanogaster than ancient, conserved ZAD-ZNF genes. We focus on the Nicknack ZAD-ZNF gene, which is evolutionarily young, poorly retained in Drosophila species, and evolves under strong positive selection. Yet we find that it is necessary for larval development in D. melanogaster. We show that Nicknack encodes a heterochromatin-localizing protein like its paralog Oddjob, also an evolutionarily dynamic yet essential ZAD-ZNF gene. We find that the divergent D. simulans Nicknack protein can still localize to D. melanogaster heterochromatin and rescue viability of female but not male Nicknack-null D. melanogaster. Our findings suggest that innovation for rapidly changing heterochromatin functions might generally explain the essentiality of many evolutionarily dynamic ZAD-ZNF genes in insects.

The early evolution of oral poliovirus vaccine is shaped by strong positive selection and tight transmission bottlenecks

The evolution of circulating vaccine-derived polioviruses (cVDPV) from components of the live-attenuated oral poliovirus vaccine (OPV) presents a major challenge to global polio eradication. This process has largely been characterized by consensus sequencing of isolates collected from routine surveillance, and little is known about the early evolution of OPV within vaccinated hosts. These early events are critical steps in the progression of OPV to cVDPV. Here, we use whole genome, high depth of coverage sequencing to define the evolutionary trajectories of monovalent type 2 OPV in a cluster-randomized trial of polio vaccines in Matlab, Bangladesh. By sequencing 416 longitudinal samples from 219 mOPV2 recipients and 81 samples from 52 household contacts, we were able to examine the extent of convergent evolution in vaccine recipients and track the amount of viral diversity transmitted to new hosts. Using time-series data from a synchronized point of vaccine administration, we identify strong positive selection of reversion mutations at three known attenuating sites within two months post-vaccination. Beyond these three recognized “gate-keeper” mutations, we identify 19 mutations that exhibit significant parallelism across vaccine recipients, providing evidence for early positive selection not previously detected by phylogenetic inference. An analysis of shared genetic variants in samples from vaccinated individuals and their household contacts suggests a tight effective bottleneck during transmission. The absence of positively selected variants among household contacts across the cohort suggests that this tight bottleneck limits the transmission of these early adaptive mutations. Together, our results highlight the distinct evolutionary dynamics of live attenuated virus vaccines and have important implications for the success of novel OPV2 and other next generation approaches.SignificanceThe emergence of circulating vaccine-derived polioviruses (cVDPV) through evolution of the oral polio vaccine (OPV) poses a significant obstacle to global eradication. Understanding the genetic changes in OPV that occur as it evolves and transmits in populations is important for preventing future cVDPV outbreaks. Little is known about the early events in VDPV evolution and the selective forces that drive them. We used high depth-of-coverage genome sequencing to assess the within-host evolutionary dynamics of monovalent type 2 OPV in a vaccine trial in Matlab, Bangladesh. We leverage longitudinal sampling from vaccine recipients and household contacts to identify mutations that arise in parallel across individuals and estimate the size of the transmission bottleneck. We find evidence for strong positive selection on key sites in the capsid and the 5’ noncoding region, many of which have not been previously identified. Our results also suggest that narrow transmission bottlenecks can constrain the spread of mutations selected within individuals. These results provide important insights into how OPV variants spread in populations and are highly relevant for ongoing poliovirus surveillance and the design of improved polio vaccines.

Large genetic diversity and strong positive selection in F-box and GPCR genes among the wild isolates of Caenorhabditis elegans

The F-box and chemosensory GPCR (csGPCR) gene families are greatly expanded in nematodes, including the model organism Caenorhabditis elegans, compared to insects and vertebrates. However, the intraspecific evolution of these two gene families in nematodes remain unexamined. In this study, we analyzed the genomic sequences of 330 recently sequenced wild isolates of C. elegans using a range of population genetics approaches. We found that F-box and csGPCR genes, especially the Srw family csGPCRs, showed much more diversity than other gene families. Population structure analysis and phylogenetic analysis divided the wild strains into eight non-Hawaiian and three Hawaiian subpopulations. Some Hawaiian strains appeared to be more ancestral than all other strains. F-box and csGPCR genes maintained a great amount of the ancestral variants in the Hawaiian subpopulation and their divergence among the non-Hawaiian subpopulations contributed significantly to population structure. These genes are mostly located at the chromosomal arms and high recombination rate correlates with their large polymorphism. Gene flow might also contribute to their diversity. Moreover, we identified signatures of strong positive selection in the F-box and csGPCR genes in the non-Hawaiian population using both neutrality tests and Extended Haplotype Homozygosity analysis. Accumulation of high frequency derived alleles in these genes were found in non-Hawaiian population, leading to divergence from the ancestral genotype found in Hawaiian strains. In summary, we found that F-box and csGPCR genes harbour a large pool of natural variants, which may be subjected to positive selection during the recent selective sweep in non-Hawaiian population. These variants are mostly mapped to the substrate-recognition domains of F-box proteins and the extracellular regions of csGPCRs, possibly resulting in advantages during adaptation by affecting protein degradation and the sensing of environmental cues, respectively.Significance statementThe small nematode Caenorhabditis elegans has emerged as an important organism in studying the genetic mechanisms of evolution. F-box and chemosensory GPCR are two of the largest gene families in C. elegans, but their intraspecific evolution within C. elegans was not studied before. In this work, using the nonsynonymous SNV data of 330 C. elegans wild isolates, we found that F-box and chemosensory GPCR genes showed larger polymorphisms and stronger positive selection than other genes. The large diversity is likely the result of rapid gene family expansion, high recombination rate, and gene flow. Analysis of subpopulation suggests that positive selection of these genes occurred most strongly in the non-Hawaiian population, which underwent a selective sweep possibly linked to human activities.

Innovation of heterochromatin functions drives rapid evolution of essential ZAD-ZNF genes in Drosophila

Contrary to prevailing dogma, evolutionarily young and dynamic genes can encode essential functions. Here, we investigate genetic innovation in ZAD-ZNF genes, which encode the most abundant class of insect transcription factors. We find that evolutionarily dynamic ZAD-ZNF genes are more likely to encode essential functions in Drosophila melanogaster than ancient, conserved ZAD-ZNF genes. To understand the basis of this unexpected correlation, we focus on the Nicknack ZAD-ZNF gene. Nicknack is an evolutionarily young, poorly retained in Drosophila species, and evolves under strong positive selection, yet we find that it is necessary for larval development in D. melanogaster. We show that Nicknack encodes a heterochromatin-localizing protein like its closely related paralog Oddjob, also an evolutionarily dynamic, essential ZAD-ZNF gene. We find that the divergent D. simulans Nicknack protein can still localize to D. melanogaster heterochromatin and rescue viability of female but not male Nicknack-null D. melanogaster. Our findings suggest that innovation for rapidly changing heterochromatin functions might provide a general explanation for the essential functions of many evolutionarily dynamic ZAD-ZNF genes in insects.

Sequence diversity and evolution of an iflavirus family associated with ticks

We studied a family of iflaviruses, a group of RNA viruses frequently found in arthropods, focusing on viruses associated with ticks. Our aim was to bring insight on the evolutionary dynamics of this group of viruses, which may interact with the biology of ticks. We explored systematically de novo RNA-Seq assemblies available for species of ticks which allowed to identify nine new genomes of iflaviruses. The phylogeny of virus sequences was not congruent with that of the tick hosts, suggesting recurrent host changes across tick genera along evolution. We identified five different variants with a complete or near-complete genome in Ixodes ricinus. These sequences were closely related, which allowed a fine-scale estimation of patterns of substitutions: we detected a strong excess of synonymous mutations suggesting evolution under strong positive selection. ISIV, a sequence found in the ISE6 cell line of Ixodes scapularis, was unexpectedly nearidentical with I. ricinus variants, suggesting a contamination of this cell line by I. ricinus material. Overall, our work constitutes a step in the understanding of the interactions between this family of viruses and ticks.