Multiplex qPCR discriminates variants of concern to enhance global surveillance of SARS-CoV-2

With the emergence of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants that may increase transmissibility and/or cause escape from immune responses, there is an urgent need for the targeted surveillance of circulating lineages. It was found that the B.1.1.7 (also 501Y.V1) variant, first detected in the United Kingdom, could be serendipitously detected by the Thermo Fisher TaqPath COVID-19 PCR assay because a key deletion in these viruses, spike Δ69–70, would cause a “spike gene target failure” (SGTF) result. However, a SGTF result is not definitive for B.1.1.7, and this assay cannot detect other variants of concern (VOC) that lack spike Δ69–70, such as B.1.351 (also 501Y.V2), detected in South Africa, and P.1 (also 501Y.V3), recently detected in Brazil. We identified a deletion in the ORF1a gene (ORF1a Δ3675–3677) in all 3 variants, which has not yet been widely detected in other SARS-CoV-2 lineages. Using ORF1a Δ3675–3677 as the primary target and spike Δ69–70 to differentiate, we designed and validated an open-source PCR assay to detect SARS-CoV-2 VOC. Our assay can be rapidly deployed in laboratories around the world to enhance surveillance for the local emergence and spread of B.1.1.7, B.1.351, and P.1.

Surveillance for SARS-CoV-2 variants is very important, but sequencing is not always practical or affordable. This study presents a multiplex qPCR that is able to distinguish among different SARS-CoV-2 variants of concern that are currently circulating.     

Type III effector diversification via both pathoadaptation and horizontal transfer in response to a coevolutionary arms race

The concept of the coevolutionary arms race holds a central position in our understanding of pathogen–host interactions. Here we identify the molecular mechanisms and follow the stepwise progression of an arms race in a natural system. We show how the evolution and function of the HopZ family of type III secreted effector proteins carried by the plant pathogen Pseudomonas syringae are influenced by a coevolutionary arms race between pathogen and host. We surveyed 96 isolates of P. syringae and identified three homologs (HopZ1, HopZ2, and HopZ3) distributed among ~45% of the strains. All alleles were sequenced and their expression was confirmed. Evolutionary analyses determined that the diverse HopZ1 homologs are ancestral to P. syringae, and have diverged via pathoadaptive mutational changes into three functional and two degenerate forms, while HopZ2 and HopZ3 have been brought into P. syringae via horizontal transfer from other ecologically similar bacteria. A PAML selection analysis revealed that the C terminus of HopZ1 is under strong positive selection. Despite the extensive genetic variation observed in this family, all three homologs have cysteine–protease activity, although their substrate specificity may vary. The introduction of the ancestral hopZ1 allele into strains harboring alternate alleles results in a resistance protein-mediated defense response in their respective hosts, which is not observed with the endogenous allele. These data indicate that the P. syringae HopZ family has undergone allelic diversification via both pathoadaptive mutational changes and horizontal transfer in response to selection imposed by the host defense system. This genetic diversity permits the pathogen to avoid host defenses while still maintaining a virulence-associated protease, thereby allowing it to thrive on its current host, while simultaneously impacting its host range.     

Natural selection in the evolution of SARS-CoV-2 in bats created a generalist virus and highly capable human pathogen

Virus host shifts are generally associated with novel adaptations to exploit the cells of the new host species optimally. Surprisingly, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has apparently required little to no significant adaptation to humans since the start of the Coronavirus Disease 2019 (COVID-19) pandemic and to October 2020. Here we assess the types of natural selection taking place in Sarbecoviruses in horseshoe bats versus the early SARS-CoV-2 evolution in humans. While there is moderate evidence of diversifying positive selection in SARS-CoV-2 in humans, it is limited to the early phase of the pandemic, and purifying selection is much weaker in SARS-CoV-2 than in related bat Sarbecoviruses. In contrast, our analysis detects evidence for significant positive episodic diversifying selection acting at the base of the bat virus lineage SARS-CoV-2 emerged from, accompanied by an adaptive depletion in CpG composition presumed to be linked to the action of antiviral mechanisms in these ancestral bat hosts. The closest bat virus to SARS-CoV-2, RmYN02 (sharing an ancestor about 1976), is a recombinant with a structure that includes differential CpG content in Spike; clear evidence of coinfection and evolution in bats without involvement of other species. While an undiscovered “facilitating” intermediate species cannot be discounted, collectively, our results support the progenitor of SARS-CoV-2 being capable of efficient human–human transmission as a consequence of its adaptive evolutionary history in bats, not humans, which created a relatively generalist virus.

A study of the natural origins of SARS-CoV-2 reveals very little adaptive evolution occurring since it emerged in humans, but strong evolutionary signals in the bat virus lineage from which SARS-CoV-2 arose. Evolution in bats involved lineage-specific depletion of CpG nucleotides (linked to host anti-viral molecules), and clear evidence of recombination across these lineages, supporting bat host species’ influence on the ancestral viruses.     

Intra‐specific body size variation in Polistes paper wasps as a response to social parasite pressure

1. Like avian brood parasites, obligate insect social parasites exploit the parental care of a host species to rear their brood, causing an evident loss of host reproductive success. This fitness cost imposes selective pressure on the host to reduce the parasite effect. A possible outcome of an evolutionary arms race is the selection of host morphological counter‐adaptations to resist parasite attacks. 2. We studied host–parasite pairs of Polistes wasps in which the fighting equipment of the parasite's body allows it to enter the host colony. 3. We searched for host morphological traits related to fighting ability that could be considered counter‐adaptations. As a host–parasite co‐evolutionary arms race can only occur where the two lineages co‐exist, we compared morphological traits of hosts belonging to populations with or without parasite pressure. We report that host foundresses belonging to populations under strong parasite pressure have a larger body size than those belonging to populations without parasite pressure. 4. Behavioural experiments carried out to test if an increase in host body size is useful to oppose parasite usurpation show that large body size foundresses exhibit a greater ability of nest defence.     

The proximal proteome of 17 SARS-CoV-2 proteins links to disrupted antiviral signaling and host translation

Viral proteins localize within subcellular compartments to subvert host machinery and promote pathogenesis. To study SARS-CoV-2 biology, we generated an atlas of 2422 human proteins vicinal to 17 SARS-CoV-2 viral proteins using proximity proteomics. This identified viral proteins at specific intracellular locations, such as association of accessary proteins with intracellular membranes, and projected SARS-CoV-2 impacts on innate immune signaling, ER-Golgi transport, and protein translation. It identified viral protein adjacency to specific host proteins whose regulatory variants are linked to COVID-19 severity, including the TRIM4 interferon signaling regulator which was found proximal to the SARS-CoV-2 M protein. Viral NSP1 protein adjacency to the EIF3 complex was associated with inhibited host protein translation whereas ORF6 localization with MAVS was associated with inhibited RIG-I 2CARD-mediated IFNB1 promoter activation. Quantitative proteomics identified candidate host targets for the NSP5 protease, with specific functional cleavage sequences in host proteins CWC22 and FANCD2. This data resource identifies host factors proximal to viral proteins in living human cells and nominates pathogenic mechanisms employed by SARS-CoV-2.     

Coexistence of Bacterial Sulfide Oxidizers, Sulfate Reducers, and Spirochetes in a Gutless Worm (Oligochaeta) from the Peru Margin

Olavius crassitunicatus is a small symbiont-bearing worm that occurs at high abundance in oxygen-deficient sediments in the East Pacific Ocean. Using comparative 16S rRNA sequence analysis and fluorescence in situ hybridization, we examined the diversity and phylogeny of bacterial symbionts in two geographically distant O. crassitunicatus populations (separated by 385 km) on the Peru margin (water depth, ~300 m). Five distinct bacterial phylotypes co-occurred in all specimens from both sites: two members of the γ-Proteobacteria (Gamma 1 and 2 symbionts), two members of the δ-Proteobacteria (Delta 1 and 2 symbionts), and one spirochete. A sixth phylotype belonging to the δ-Proteobacteria (Delta 3 symbiont) was found in only one of the two host populations. Three of the O. crassitunicatus bacterial phylotypes are closely related to symbionts of other gutless oligochaete species; the Gamma 1 phylotype is closely related to sulfide-oxidizing symbionts of Olavius algarvensis, Olavius loisae, and Inanidrilus leukodermatus, the Delta 1 phylotype is closely related to sulfate-reducing symbionts of O. algarvensis, and the spirochete is closely related to spirochetal symbionts of O. loisae. In contrast, the Gamma 2 phylotype and the Delta 2 and 3 phylotypes belong to novel lineages that are not related to other bacterial symbionts. Such a phylogenetically diverse yet highly specific and stable association in which multiple bacterial phylotypes coexist within a single host has not been described previously for marine invertebrates.     

Genetics of tan spot resistance in wheat

Tan spot is a devastating foliar disease of wheat caused by the necrotrophic fungal pathogen Pyrenophora tritici-repentis. Much has been learned during the past two decades about the genetics of wheat–P. tritici-repentis interactions. Research has shown that the fungus produces at least three host-selective toxins (HSTs), known as Ptr ToxA, Ptr ToxB, and Ptr ToxC, that interact directly or indirectly with the products of the dominant host genes Tsn1, Tsc2, and Tsc1, respectively. The recent cloning and characterization of Tsn1 provided strong evidence that the pathogen utilizes HSTs to subvert host resistance mechanisms to cause disease. However, in addition to host–HST interactions, broad-spectrum, race non-specific resistance QTLs and recessively inherited qualitative ‘resistance’ genes have been identified. Molecular markers suitable for marker-assisted selection against HST sensitivity genes and for race non-specific resistance QTLs have been developed and used to generate adapted germplasm with good levels of tan spot resistance. Future research is needed to identify novel HSTs and corresponding host sensitivity genes, determine if the recessively inherited resistance genes are HST insensitivities, extend the current race classification system to account for new HSTs, and determine the molecular basis of race non-specific resistance QTLs and their relationships with host–HST interactions at the molecular level. Necrotrophic pathogens such as P. tritici-repentis are likely to become increasingly significant under a changing global climate making it imperative to further characterize the wheat–P. tritici-repentis pathosystem and develop tan spot resistant wheat varieties.     

Clinical and genomic data of sars-cov-2 detected in maternal–fetal interface during the first wave of infection in Brazil

Brazil has the highest SARS-CoV-2 case-fatality rate in pregnant women in the Americas. In this study, clinical and virological findings of five mildly symptomatic pregnant women and their infected fetuses/newborns treated at a referral hospital for COVID19-pregnant women in Midwestern Brazil are reported. Mother and fetal samples were tested by RT-qPCR, ECLIA and Illumina MiSeq sequencing. From the five cases, one resulted in spontaneous abortion, one was stillborn, two were preterm births and one full-term birth. Maternal and fetal placenta, newborn and stillborn secretions were SARS-CoV-2+; one neonate developed ground-glass opacities in his lungs. One neonate's umbilical cord was IgG+ and all were IgM negative upon hospital discharge. Genomes recovered from two placentas belong to the B.1.1.28 and B.1.1.33 lineages and present nonsynonymous mutations associated with virus fitness and infectivity; other not frequently reported mutations (B.1.1.33: NSP3 V2090G, M A2S and ORF3ab S253P and Y264N; B.1.1.28: NSP3 E995D, NSP12 R240K, NSP14H1897Y and in ORF7b V21F) were found in proteins involved in viral replication, viral induction of apoptosis, viral interference on interferon and on NF-Κβ pathways. Phylogeny indicates the south of Brazil as the possible origin of these lineages circulating in Mato Grosso State. These findings contribute to describe SARS-CoV-2 infection and outcomes in pregnant women and their fetuses, at any stage of gestation and even in mild symptomatic cases.     

Diversification of Two Lineages of Symbiotic Photobacterium

Understanding of processes driving bacterial speciation requires examination of closely related, recently diversified lineages. To gain an insight into diversification of bacteria, we conducted comparative genomic analysis of two lineages of bioluminescent symbionts, Photobacterium leiognathi and ‘P. mandapamensis’. The two lineages are evolutionary and ecologically closely related. Based on the methods used in bacterial taxonomy for classification of new species (DNA-DNA hybridization and ANI), genetic relatedness of the two lineages is at a cut-off point for species delineation. In this study, we obtained the whole genome sequence of a representative P. leiognathi strain lrivu.4.1, and compared it to the whole genome sequence of ‘P. mandapamensis’ svers.1.1. Results of the comparative genomic analysis suggest that P. leiognathi has a more plastic genome and acquired genes horizontally more frequently than ‘P. mandapamensis’. We predict that different rates of recombination and gene acquisition contributed to diversification of the two lineages. Analysis of lineage-specific sequences in 25 strains of P. leiognathi and ‘P. mandapamensis’ found no evidence that bioluminescent symbioses with specific host animals have played a role in diversification of the two lineages.     

Distribution of Barley yellow dwarf virus-PAV in the Sub-Antarctic Kerguelen Islands and Characterization of Two New Luteovirus Species

A systematic search for viral infection was performed in the isolated Kerguelen Islands, using a range of polyvalent genus-specific PCR assays. Barley yellow dwarf virus (BYDV) was detected in both introduced and native grasses such as Poa cookii. The geographical distribution of BYDV and its prevalence in P. cookii were analyzed using samples collected from various sites of the archipelago. We estimate the average prevalence of BYDV to be 24.9% in P. cookii, with significant variability between sites. BYDV genetic diversity was assessed using sequence information from two genomic regions: the P3 open reading frame (ORF) (encoding the coat protein) and the hypervariable P6 ORF region. The phylogenetic analysis in the P3 region showed that BYDV sequences segregate into three major lineages, the most frequent of which (Ker-I cluster) showed close homology with BYDV-PAV-I isolates and had very low intra-lineage diversity (0.6%). A similarly low diversity was also recorded in the hypervariable P6 region, suggesting that Ker-I isolates derive from the recent introduction of BYDV-PAV-I. Divergence time estimation suggests that BYDV-PAV-I was likely introduced in the Kerguelen environment at the same time frame as its aphid vector, Rhopalosiphum padi, whose distribution shows good overlap with that of BYDV-Ker-I. The two other lineages show more than 22% amino acid divergence in the P3 region with other known species in the BYDV species complex, indicating that they represent distinct BYDV species. Using species-specific amplification primers, the distribution of these novel species was analyzed. The high prevalence of BYDV on native Poaceae and the presence of the vector R. padi, raises the question of its impact on the vulnerable plant communities of this remote ecosystem.     

Different rates of nucleotide substitutions in Wolbachia endosymbionts of arthropods and nematodes: arms race or host shifts?

The genus Wolbachia encompasses intracellular bacteria found in arthropods and in filarial nematodes. In arthropods, Wolbachia is primarily a reproductive parasite and shows relatively frequent horizontal transfer between host species, while in nematodes it appears to be a mutualist and is strictly vertically transmitted. We can expect that different selective pressures are acting on their genomes. Here we present an analysis of three Wolbachia genes, wsp, ftsZ and dnaA. In wsp of arthropod Wolbachia, an excess of non-synonymous substitutions was observed, providing evidence for positive selection. In nematode Wolbachia, no evidence for positive selection was found. Pressure for amino acid variation in wsp of arthropod Wolbachia could derive either from an arms race with the host or from the occurrence of more frequent hosts shifts due to horizontal transmission. In nematode Wolbachia, the lack of positively selected sites could result from the absence of an arms race, or from the homogeneity of the biochemical environment they exist in (ensured by strict vertical transmission). In ftsZ minor differences in substitution patterns were observed between arthropod and nematode Wolbachia, only in the 3'-portion of the gene. dnaA showed comparable patterns of variation in both lineages, with evidence for strong conservation.     

Host associations and evolutionary relationships of avian blood parasites from West Africa

The host specificity of blood parasites recovered from a survey of 527 birds in Cameroon and Gabon was examined at several levels within an evolutionary framework. Unique mitochondrial lineages of Haemoproteus were recovered from an average of 1.3 host species (maximum = 3) and 1.2 host families (maximum = 3) while lineages of Plasmodium were recovered from an average of 2.5 species (maximum = 27) and 1.6 families (maximum = 9). Averaged within genera, lineages of both Plasmodium and Haemoproteus were constrained in their host distribution relative to random expectations. However, while several individual lineages within both genera exhibited significant host constraint, host breadth varied widely among related lineages, particularly within the genus Plasmodium. Several lineages of Plasmodium exhibited extreme generalist host-parasitism strategies while other lineages appeared to have been constrained to certain host families over recent evolutionary history. Sequence data from two nuclear genes recovered from a limited sample of Plasmodium parasites indicated that, at the resolution of this study, inferences regarding host breadth were unlikely to be grossly affected by the use of parasite mitochondrial lineages as a proxy for biological species. The use of divergent host-parasitism strategies among closely related parasite lineages suggests that host range is a relatively labile character. Since host specificity may also influence parasite virulence, these results argue for considering the impact of haematozoa on avian hosts on a lineage-specific basis.     

Stepwise Evolution and Exceptional Conservation of ORF1a/b Overlap in Coronaviruses

The programmed frameshift element (PFE) rerouting translation from ORF1a to ORF1b is essential for the propagation of coronaviruses. The combination of genomic features that make up PFE—the overlap between the two reading frames, a slippery sequence, as well as an ensemble of complex secondary structure elements—places severe constraints on this region as most possible nucleotide substitution may disrupt one or more of these elements. The vast amount of SARS-CoV-2 sequencing data generated within the past year provides an opportunity to assess the evolutionary dynamics of PFE in great detail. Here, we performed a comparative analysis of all available coronaviral genomic data available to date. We show that the overlap between ORF1a and ORF1b evolved as a set of discrete 7, 16, 22, 25, and 31 nucleotide stretches with a well-defined phylogenetic specificity. We further examined sequencing data from over 1,500,000 complete genomes and 55,000 raw read data sets to demonstrate exceptional conservation and detect signatures of selection within the PFE region.     

Looking for pathways related to COVID-19: confirmation of pathogenic mechanisms by SARS-CoV-2–host interactome

In the last months, many studies have clearly described several mechanisms of SARS-CoV-2 infection at cell and tissue level, but the mechanisms of interaction between host and SARS-CoV-2, determining the grade of COVID-19 severity, are still unknown. We provide a network analysis on protein–protein interactions (PPI) between viral and host proteins to better identify host biological responses, induced by both whole proteome of SARS-CoV-2 and specific viral proteins. A host-virus interactome was inferred, applying an explorative algorithm (Random Walk with Restart, RWR) triggered by 28 proteins of SARS-CoV-2. The analysis of PPI allowed to estimate the distribution of SARS-CoV-2 proteins in the host cell. Interactome built around one single viral protein allowed to define a different response, underlining as ORF8 and ORF3a modulated cardiovascular diseases and pro-inflammatory pathways, respectively. Finally, the network-based approach highlighted a possible direct action of ORF3a and NS7b to enhancing Bradykinin Storm. This network-based representation of SARS-CoV-2 infection could be a framework for pathogenic evaluation of specific clinical outcomes. We identified possible host responses induced by specific proteins of SARS-CoV-2, underlining the important role of specific viral accessory proteins in pathogenic phenotypes of severe COVID-19 patients.Subject terms: Protein-protein interaction networks, Viral infection     

Structural modeling of protein ensembles between E3 RING ligases and SARS-CoV-2: The role of zinc binding domains

BackgroundThe ubiquitin system is a modification process with many different cellular functions including immune signaling and antiviral functions. E3 ubiquitin ligases are enzymes that recruit an E2 ubiquitin-conjugating enzyme bound to ubiquitin in order to catalyze the transfer of ubiquitin from the E2 to a protein substrate. The RING E3s, the most abundant type of ubiquitin ligases, are characterized by a zinc (II)-binding domain called RING (Really Interesting New Gene). Viral replication requires modifying and hijacking key cellular pathways within host cells such as cellular ubiquitination. There are well-established examples where a viral proteins bind to RING E3s, redirecting them to degrade otherwise long-lived host proteins or inhibiting E3’s ubiquitination activity. Recently, three binary interactions between SARS-CoV-2 proteins and innate human immune signaling Ε3 RING ligases: NSP15-RNF41, ORF3a-TRIM59 and NSP9-MIB1 have been experimentally established.MethodsIn this work, we have investigated the mode of the previous experimentally supported NSP15-RNF41, ORF3a,-TRIM59 and NSP9-MIB1 binary interactions by in silico methodologies intending to provide structural insights of E3-virus interplay that can help identify potential inhibitors that could block SARS-CoV-2 infection of immune cells.ConclusionIn silico methodologies have shown that the above human E3 ligases interact with viral partners through their Zn(II) binding domains. This RING mediated formation of stable SARS-CoV-2-E3 complexes indicates a critical structural role of RING domains in immune system disruption by SARS-CoV-2-infection.Data AvailabilityThe data used to support the findings of this research are included within the article and are labeled with references.     

Wheat Domestication Accelerated Evolution and Triggered Positive Selection in the β-Xylosidase Enzyme of Mycosphaerella graminicola

Plant cell wall degrading enzymes (PCWDEs) of plant pathogens are receiving increasing interest for their potential to trigger plant defense reactions. In an antagonistic co-evolutionary arms race between host and pathogen, PCWDEs could be under strong selection. Here, we tested the hypothesis that PCWDEs in the fungal wheat pathogen Mycosphaerella graminicola have been positively selected by analyzing ratios of non-synonymous and synonymous nucleotide changes in the genes encoding these enzymes. Analyses of five PCWDEs demonstrated that one (β-xylosidase) has been under strong positive selection and experienced an accelerated rate of evolution. In contrast, PCWDEs in the closest relatives of M. graminicola collected from wild grasses did not show evidence for selection or deviation from a molecular clock. Since the genealogical divergence of M. graminicola from these latter species coincided with the onset of agriculture, we hypothesize that the recent domestication of the host plant and/or agricultural practices triggered positive selection in β-xylosidase and that this enzyme played a key role in the emergence of a host-specialized pathogen.     

Evaluation of the relative virulence of novel SARS-CoV-2 variants: a retrospective cohort study in Ontario,Canada

Background:Between February and June 2021, the initial wild-type strains of SARS-CoV-2 were supplanted in Ontario, Canada, by new variants of concern (VOCs), first those with the N501Y mutation (i.e., Alpha/B1.1.17, Beta/B.1.351 and Gamma/P.1 variants) and then the Delta/B.1.617 variant. The increased transmissibility of these VOCs has been documented, but knowledge about their virulence is limited. We used Ontario’s COVID-19 case data to evaluate the virulence of these VOCs compared with non-VOC SARS-CoV-2 strains, as measured by risk of hospitalization, intensive care unit (ICU) admission and death.Methods:We created a retrospective cohort of people in Ontario who tested positive for SARS-CoV-2 and were screened for VOCs, with dates of test report between Feb. 7 and June 27, 2021. We constructed mixed-effect logistic regression models with hospitalization, ICU admission and death as outcome variables. We adjusted models for age, sex, time, vaccination status, comorbidities and pregnancy status. We included health units as random intercepts.Results:Our cohort included 212 326 people. Compared with non-VOC SARS-CoV-2 strains, the adjusted elevation in risk associated with N501Y-positive variants was 52% (95% confidence interval [CI] 42%–63%) for hospitalization, 89% (95% CI 67%–117%) for ICU admission and 51% (95% CI 30%–78%) for death. Increased risk with the Delta variant was more pronounced at 108% (95% CI 78%–140%) for hospitalization, 235% (95% CI 160%–331%) for ICU admission and 133% (95% CI 54%–231%) for death.Interpretation:The increasing virulence of SARS-CoV-2 VOCs will lead to a considerably larger, and more deadly, pandemic than would have occurred in the absence of the emergence of VOCs.

Novel SARS-CoV-2 variants of concern (VOCs), including viral lineages carrying the N501Y (Alpha/B.1.1.7) or both the N501Y and E484K mutations (Beta/B.1.351 and Gamma/P.1), were first identified in Ontario, Canada, in December 2020.¹ Although initially uncommon in Ontario, these VOCs outcompeted earlier SARS-CoV-2 lineages and, as of late April 2021, were responsible for almost all new infections in Ontario, with Alpha the most prevalent lineage.¹ In April 2021, the B.1.617.2 variant, now known as Delta under the revised nomenclature from the World Health Organization, emerged in the province, outcompeted earlier VOCs and, by July 2021, represented most infections in the province.^2,3This serial replacement by emerging variants reflects progressively higher effective reproduction numbers that allow novel variants to outcompete previously dominant strains in the face of identical measures to control spread of infection.^4–6 However, VOCs are also concerning because emerging evidence points to increased virulence, with increased risk of hospitalization, intensive care unit (ICU) admission and death, after adjustment for age and other predictive factors among patients with VOC infections.^7–10 Although the increased virulence of strains with the N501Y mutation relative to strains that lack this mutation has been described,^7–9 only limited information is available on the virulence of infection with the Delta variant, relative to earlier N501Y-positive VOCs (i.e., Alpha, Beta and Gamma).^10–12 Our objectives were to evaluate the virulence of N501Y-positive variants relative to earlier SARS-CoV-2 lineages and to evaluate the virulence of the Delta variant of SARS-CoV-2 relative to N501Y-positive VOCs using Ontario’s COVID-19 case data.     

Global Population Genomics of Two Subspecies of Cryptosporidium hominis during 500 Years of Evolution

Cryptosporidiosis is a major global health problem and a primary cause of diarrhea, particularly in young children in low- and middle-income countries (LMICs). The zoonotic Cryptosporidium parvum and anthroponotic Cryptosporidium hominis cause most human infections. Here, we present a comprehensive whole-genome study of C. hominis, comprising 114 isolates from 16 countries within five continents. We detect two lineages with distinct biology and demography, which diverged circa 500 years ago. We consider these lineages two subspecies and propose the names C. hominis hominis and C. hominis aquapotentis (gp60 subtype IbA10G2). In our study, C. h. hominis is almost exclusively represented by isolates from LMICs in Africa and Asia and appears to have undergone recent population contraction. In contrast, C. h. aquapotentis was found in high-income countries, mainly in Europe, North America, and Oceania, and appears to be expanding. Notably, C. h. aquapotentis is associated with high rates of direct human-to-human transmission, which may explain its success in countries with well-developed environmental sanitation infrastructure. Intriguingly, we detected genomic regions of introgression following secondary contact between the subspecies. This resulted in high diversity and divergence in genomic islands of putative virulence genes, including muc5 (CHUDEA2_430) and a hypothetical protein (CHUDEA6_5270). This diversity is maintained by balancing selection, suggesting a co-evolutionary arms race with the host. Finally, we find that recent gene flow from C. h. aquapotentis to C. h. hominis, likely associated with increased human migration, maybe driving the evolution of more virulent C. hominis variants.