Recent expansion of metabolic versatility in Diplonema papillatum, the model species of a highly speciose group of marine eukaryotes.

BACKGROUND: Diplonemid flagellates are among the most abundant and species-rich of known marine microeukaryotes, colonizing all habitats, depths, and geographic regions of the world ocean. However, little is known about their genomes, biology, and ecological role. RESULTS: We present the first nuclear genome sequence from a diplonemid, the type species Diplonema papillatum. The ~ 280-Mb genome assembly contains about 32,000 protein-coding genes, likely co-transcribed in groups of up to 100. Gene clusters are separated by long repetitive regions that include numerous transposable elements, which also reside within introns. Analysis of gene-family evolution reveals that the last common diplonemid ancestor underwent considerable metabolic expansion. D. papillatum-specific gains of carbohydrate-degradation capability were apparently acquired via horizontal gene transfer. The predicted breakdown of polysaccharides including pectin and xylan is at odds with reports of peptides being the predominant carbon source of this organism. Secretome analysis together with feeding experiments suggest that D. papillatum is predatory, able to degrade cell walls of live microeukaryotes, macroalgae, and water plants, not only for protoplast feeding but also for metabolizing cell-wall carbohydrates as an energy source. The analysis of environmental barcode samples shows that D. papillatum is confined to temperate coastal waters, presumably acting in bioremediation of eutrophication. CONCLUSIONS: Nuclear genome information will allow systematic functional and cell-biology studies in D. papillatum. It will also serve as a reference for the highly diverse diplonemids and provide a point of comparison for studying gene complement evolution in the sister group of Kinetoplastida, including human-pathogenic taxa.

Impact of temperature on the affinity of SARS-CoV-2 Spike glycoprotein for host ACE2.

The seasonal nature of outbreaks of respiratory viral infections with increased transmission during low temperatures has been well established. Accordingly, temperature has been suggested to play a role on the viability and transmissibility of SARS-CoV-2, the virus responsible for the COVID-19 pandemic. The receptor-binding domain (RBD) of the Spike glycoprotein is known to bind to its host receptor angiotensin-converting enzyme 2 (ACE2) to initiate viral fusion. Using biochemical, biophysical, and functional assays to dissect the effect of temperature on the receptor-Spike interaction, we observed a significant and stepwise increase in RBD-ACE2 affinity at low temperatures, resulting in slower dissociation kinetics. This translated into enhanced interaction of the full Spike glycoprotein with the ACE2 receptor and higher viral attachment at low temperatures. Interestingly, the RBD N501Y mutation, present in emerging variants of concern (VOCs) that are fueling the pandemic worldwide (including the B.1.1.7 (α) lineage), bypassed this requirement. This data suggests that the acquisition of N501Y reflects an adaptation to warmer climates, a hypothesis that remains to be tested.

Changing Patterns of Disease Severity in Blastomyces dermatitidis Infection, Quebec, Canada.

This retrospective multicenter cohort study assessed temporal changes in the severity and mortality rate of blastomycosis in Quebec, Canada, and identified risk factors for death in patients with blastomycosis in 1988-2016. The primary outcome was 90-day all-cause deaths. Among 185 patients, 122 (66%) needed hospitalization and 30 (16%) died. We noted increases in the proportion of severe cases, in age at diagnosis and in the proportion of diabetic and immunocompromised patients over time. Independent risk factors for death were age (adjusted odds ratio [aOR] 1.04, 95% CI 1.00-1.07), immunosuppression (aOR 4.2, 95% CI 1.5-11.6), and involvement of >2 lung lobes (aOR 5.3, 95% CI 1.9-14.3). There was no association between the Blastomyces genotype group and all-cause mortality. The proportion of severe cases of blastomycosis has increased in Quebec over the past 30 years, partially explained by the higher number of immunosuppressed patients.

Small RNA-mediated regulation of DNA dosage in the ciliate Oxytricha.

Dicer-dependent small noncoding RNAs play important roles in gene regulation in a wide variety of organisms. Endogenous small interfering RNAs (siRNAs) are part of an ancient pathway of transposon control in plants and animals. The ciliate, Oxytricha trifallax, has approximately 16,000 gene-sized chromosomes in its somatic nucleus. Long noncoding RNAs establish high ploidy levels at the onset of sexual development, but the factors that regulate chromosome copy numbers during cell division and growth have been a mystery. We report a novel function of a class of Dicer (Dcl-1)- and RNA-dependent RNA polymerase (RdRP)-dependent endogenous small RNAs in regulating chromosome copy number and gene dosage in O. trifallax Asexually growing populations express an abundant class of 21-nt sRNAs that map to both coding and noncoding regions of most chromosomes. These sRNAs are bound to chromatin and their levels surprisingly do not correlate with mRNA levels. Instead, the levels of these small RNAs correlate with genomic DNA copy number. Reduced sRNA levels in dcl-1 or rdrp mutants lead to concomitant reduction in chromosome copy number. Furthermore, these cells show no signs of transposon activation, but instead display irregular nuclear architecture and signs of replication stress. In conclusion, Oxytricha Dcl-1 and RdRP-dependent small RNAs that derive from the somatic nucleus contribute to the maintenance of gene dosage, possibly via a role in DNA replication, offering a novel role for these small RNAs in eukaryotes.

Prevalence and mechanisms of azole resistance in clinical isolates of Aspergillus section Fumigati species in a Canadian tertiary care centre, 2000 to 2013.

OBJECTIVES: Azole resistance among Aspergillus fumigatus isolates is a growing concern worldwide. Induction of mutations during azole therapy, environment-acquired mutations caused by azole fungicides and intrinsic resistance of cryptic Fumigati species all contribute to the burden of resistance. However, there is a lack of data in Canada on this emerging threat. METHODS: To gain insights into the magnitude and mechanisms of resistance, a 14 year collection of Aspergillus section Fumigati comprising 999 isolates from 807 patients at a Montreal hospital was screened for azole resistance, and resistance mechanisms were investigated with the combined use of genome sequencing, 3D modelling and phenotypic efflux pump assays. RESULTS: Overall azole resistance was low (4/807 patients; 0.5%). A single azole-resistant A. fumigatus sensu stricto strain, isolated from a patient with pulmonary aspergillosis, displayed efflux-pump-mediated resistance. Three patients were colonized or infected with azole-resistant cryptic Fumigati species (one Aspergillus thermomutatus, one Aspergillus lentulus and one Aspergillus turcosus). Evidence is presented that azole resistance is efflux-pump-mediated in the A. turcosus isolate, but not in the A. lentulus and A. thermomutatus isolates. CONCLUSIONS: Azole resistance is rare in our geographic area and currently driven by cryptic Fumigati species. Continued surveillance of emergence of resistance is warranted.

Capture of complete ciliate chromosomes in single sequencing reads reveals widespread chromosome isoforms.

BACKGROUND: Whole-genome shotgun sequencing, which stitches together millions of short sequencing reads into a single genome, ushered in the era of modern genomics and led to a rapid expansion of the number of genome sequences available. Nevertheless, assembly of short reads remains difficult, resulting in fragmented genome sequences. Ultimately, only a sequencing technology capable of capturing complete chromosomes in a single run could resolve all ambiguities. Even "third generation" sequencing technologies produce reads far shorter than most eukaryotic chromosomes. However, the ciliate Oxytricha trifallax has a somatic genome with thousands of chromosomes averaging only 3.2 kbp, making it an ideal candidate for exploring the benefits of sequencing whole chromosomes without assembly. RESULTS: We used single-molecule real-time sequencing to capture thousands of complete chromosomes in single reads and to update the published Oxytricha trifallax JRB310 genome assembly. In this version, over 50% of the completed chromosomes with two telomeres derive from single reads. The improved assembly includes over 12,000 new chromosome isoforms, and demonstrates that somatic chromosomes derive from variable rearrangements between somatic segments encoded up to 191,000 base pairs away. However, while long reads reduce the need for assembly, a hybrid approach that supplements long-read sequencing with short reads for error correction produced the most complete and accurate assembly, overall. CONCLUSIONS: This assembly provides the first example of complete eukaryotic chromosomes captured by single sequencing reads and demonstrates that traditional approaches to genome assembly can mask considerable structural variation.

Genes in Hiding.

Unrecognizable genes are an unsettling problem in genomics. Here, we survey the various types of cryptic genes and the corresponding deciphering strategies employed by cells. Encryption that renders genes substantially different from homologs in other species includes sequence substitution, insertion, deletion, fragmentation plus scrambling, and invasion by mobile genetic elements. Cells decode cryptic genes at the DNA, RNA or protein level. We will focus on a recently discovered case of unparalleled encryption involving massive gene fragmentation and nucleotide deletions and substitutions, occurring in the mitochondrial genome of a poorly understood protist group, the diplonemids. This example illustrates that comprehensive gene detection requires not only auxiliary sequence information - transcriptome and proteome data - but also knowledge about a cell's deciphering arsenal.

Transformation of Diplonema papillatum, the type species of the highly diverse and abundant marine microeukaryotes Diplonemida (Euglenozoa).

Diplonema papillatum is the type species of diplonemids, which are among the most abundant and diverse heterotrophic microeukaryotes in the world's oceans. Diplonemids are also known for a unique form of post-transcriptional processing in mitochondria. However, the lack of reverse genetics methodologies in these protists has hampered elucidation of their cellular and molecular biology. Here we report a protocol for D. papillatum transformation. We have identified several antibiotics to which D. papillatum is sensitive and thus are suitable selectable markers, and focus in particular on puromycin. Constructs were designed encoding antibiotic resistance markers, fluorescent tags, and additional genomic sequences from D. papillatum to facilitate vector integration into chromosomes. We established conditions for effective electroporation, and demonstrate that electroporated constructs can be stably integrated in the D. papillatum nuclear genome. In D. papillatum transformants, the heterologous puromycin resistance gene is transcribed into mRNA and translated into protein, as determined by Southern hybridization, reverse transcription, and Western blot analyses. This is the first documented case of transformation in a euglenozoan protist outside the well-studied kinetoplastids, making D. papillatum a genetically tractable organism and potentially a model system for marine microeukaryotes.

Novel modes of RNA editing in mitochondria.

Gene structure and expression in diplonemid mitochondria are unparalleled. Genes are fragmented in pieces (modules) that are separately transcribed, followed by the joining of module transcripts to contiguous RNAs. Some instances of unique uridine insertion RNA editing at module boundaries were noted, but the extent and potential occurrence of other editing types remained unknown. Comparative analysis of deep transcriptome and genome data from Diplonema papillatum mitochondria reveals ∼220 post-transcriptional insertions of uridines, but no insertions of other nucleotides nor deletions. In addition, we detect in total 114 substitutions of cytosine by uridine and adenosine by inosine, amassed into unusually compact clusters. Inosines in transcripts were confirmed experimentally. This is the first report of adenosine-to-inosine editing of mRNAs and ribosomal RNAs in mitochondria. In mRNAs, editing causes mostly amino-acid additions and non-synonymous substitutions; in ribosomal RNAs, it permits formation of canonical secondary structures. Two extensively edited transcripts were compared across four diplonemids. The pattern of uridine-insertion editing is strictly conserved, whereas substitution editing has diverged dramatically, but still rendering diplonemid proteins more similar to other eukaryotic orthologs. We posit that RNA editing not only compensates but also sustains, or even accelerates, ultra-rapid evolution of genome structure and sequence in diplonemid mitochondria.

Post-transcriptional mending of gene sequences: Looking under the hood of mitochondrial gene expression in diplonemids.

The instructions to make proteins and structural RNAs are laid down in gene sequences. Yet, in certain instances, these primary instructions need to be modified considerably during gene expression, most often at the transcript level. Here we review a case of massive post-transcriptional revisions via trans-splicing and RNA editing, a phenomenon occurring in mitochondria of a recently recognized protist group, the diplonemids. As of now, the various post-transcriptional steps have been cataloged in detail, but how these processes function is still unknown. Since genetic manipulation techniques such as gene replacement and RNA interference have not yet been established for these organisms, alternative strategies have to be deployed. Here, we discuss the experimental and bioinformatics approaches that promise to unravel the molecular machineries of trans-splicing and RNA editing in Diplonema mitochondria.

Trans-splicing and RNA editing of LSU rRNA in Diplonema mitochondria.

Mitochondrial ribosomal RNAs (rRNAs) often display reduced size and deviant secondary structure, and sometimes are fragmented, as are their corresponding genes. Here we report a mitochondrial large subunit rRNA (mt-LSU rRNA) with unprecedented features. In the protist Diplonema, the rnl gene is split into two pieces (modules 1 and 2, 534- and 352-nt long) that are encoded by distinct mitochondrial chromosomes, yet the rRNA is continuous. To reconstruct the post-transcriptional maturation pathway of this rRNA, we have catalogued transcript intermediates by deep RNA sequencing and RT-PCR. Gene modules are transcribed separately. Subsequently, transcripts are end-processed, the module-1 transcript is polyuridylated and the module-2 transcript is polyadenylated. The two modules are joined via trans-splicing that retains at the junction ∼ 26 uridines, resulting in an extent of insertion RNA editing not observed before in any system. The A-tail of trans-spliced molecules is shorter than that of mono-module 2, and completely absent from mitoribosome-associated mt-LSU rRNA. We also characterize putative antisense transcripts. Antisense-mono-modules corroborate bi-directional transcription of chromosomes. Antisense-mt-LSU rRNA, if functional, has the potential of guiding concomitantly trans-splicing and editing of this rRNA. Together, these findings open a window on the investigation of complex regulatory networks that orchestrate multiple and biochemically diverse post-transcriptional events.

Three-dimensional structure model and predicted ATP interaction rewiring of a deviant RNA ligase 2.

BACKGROUND: RNA ligases 2 are scarce and scattered across the tree of life. Two members of this family are well studied: the mitochondrial RNA editing ligase from the parasitic trypanosomes (Kinetoplastea), a promising drug target, and bacteriophage T4 RNA ligase 2, a workhorse in molecular biology. Here we report the identification of a divergent RNA ligase 2 (DpRNL) from Diplonema papillatum (Diplonemea), a member of the kinetoplastids' sister group. METHODS: We identified DpRNL with methods based on sensitive hidden Markov Model. Then, using homology modeling and molecular dynamics simulations, we established a three dimensional structure model of DpRNL complexed with ATP and Mg2+. RESULTS: The 3D model of Diplonema was compared with available crystal structures from Trypanosoma brucei, bacteriophage T4, and two archaeans. Interaction of DpRNL with ATP is predicted to involve double π-stacking, which has not been reported before in RNA ligases. This particular contact would shift the orientation of ATP and have considerable consequences on the interaction network of amino acids in the catalytic pocket. We postulate that certain canonical amino acids assume different functional roles in DpRNL compared to structurally homologous residues in other RNA ligases 2, a reassignment indicative of constructive neutral evolution. Finally, both structure comparison and phylogenetic analysis show that DpRNL is not specifically related to RNA ligases from trypanosomes, suggesting a unique adaptation of the latter for RNA editing, after the split of diplonemids and kinetoplastids. CONCLUSION: Homology modeling and molecular dynamics simulations strongly suggest that DpRNL is an RNA ligase 2. The predicted innovative reshaping of DpRNL's catalytic pocket is worthwhile to be tested experimentally.

Intra-Host Evolution Analyses in an Immunosuppressed Patient Supports SARS-CoV-2 Viral Reservoir Hypothesis.

Throughout the SARS-CoV-2 pandemic, several variants of concern (VOCs) have been identified, many of which share recurrent mutations in the spike glycoprotein's receptor-binding domain (RBD). This region coincides with known epitopes and can therefore have an impact on immune escape. Protracted infections in immunosuppressed patients have been hypothesized to lead to an enrichment of such mutations and therefore drive evolution towards VOCs. Here, we present the case of an immunosuppressed patient that developed distinct populations with immune escape mutations throughout the course of their infection. Notably, by investigating the co-occurrence of substitutions on individual sequencing reads in the RBD, we found quasispecies harboring mutations that confer resistance to known monoclonal antibodies (mAbs) such as S:E484K and S:E484A. These mutations were acquired without the patient being treated with mAbs nor convalescent sera and without them developing a detectable immune response to the virus. We also provide additional evidence for a viral reservoir based on intra-host phylogenetics, which led to a viral substrain that evolved elsewhere in the patient's body, colonizing their upper respiratory tract (URT). The presence of SARS-CoV-2 viral reservoirs can shed light on protracted infections interspersed with periods where the virus is undetectable, and potential explanations for long-COVID cases.

The Canadian VirusSeq Data Portal & Duotang: open resources for SARS-CoV-2 viral sequences and genomic epidemiology.

The COVID-19 pandemic led to a large global effort to sequence SARS-CoV-2 genomes from patient samples to track viral evolution and inform public health response. Millions of SARS-CoV-2 genome sequences have been deposited in global public repositories. The Canadian COVID-19 Genomics Network (CanCOGeN - VirusSeq), a consortium tasked with coordinating expanded sequencing of SARS-CoV-2 genomes across Canada early in the pandemic, created the Canadian VirusSeq Data Portal, with associated data pipelines and procedures, to support these efforts. The goal of VirusSeq was to allow open access to Canadian SARS-CoV-2 genomic sequences and enhanced, standardized contextual data that were unavailable in other repositories and that meet FAIR standards (Findable, Accessible, Interoperable and Reusable). In addition, the Portal data submission pipeline contains data quality checking procedures and appropriate acknowledgement of data generators that encourages collaboration. From inception to execution, the portal was developed with a conscientious focus on strong data governance principles and practices. Extensive efforts ensured a commitment to Canadian privacy laws, data security standards, and organizational processes. This Portal has been coupled with other resources like Viral AI and was further leveraged by the Coronavirus Variants Rapid Response Network (CoVaRR-Net) to produce a suite of continually updated analytical tools and notebooks. Here we highlight this Portal, including its contextual data not available elsewhere, and the 'Duotang', a web platform that presents key genomic epidemiology and modeling analyses on circulating and emerging SARS-CoV-2 variants in Canada. Duotang presents dynamic changes in variant composition of SARS-CoV-2 in Canada and by province, estimates variant growth, and displays complementary interactive visualizations, with a text overview of the current situation. The VirusSeq Data Portal and Duotang resources, alongside additional analyses and resources computed from the Portal (COVID-MVP, CoVizu), are all open-source and freely available. Together, they provide an updated picture of SARS-CoV-2 evolution to spur scientific discussions, inform public discourse, and support communication with and within public health authorities. They also serve as a framework for other jurisdictions interested in open, collaborative sequence data sharing and analyses.

Antigenicity of the Mu (B.1.621) and A.2.5 SARS-CoV-2 Spikes.

The rapid emergence of SARS-CoV-2 variants is fueling the recent waves of the COVID-19 pandemic. Here, we assessed ACE2 binding and antigenicity of Mu (B.1.621) and A.2.5 Spikes. Both these variants carry some mutations shared by other emerging variants. Some of the pivotal mutations such as N501Y and E484K in the receptor-binding domain (RBD) detected in B.1.1.7 (Alpha), B.1.351 (Beta) and P.1 (Gamma) are now present within the Mu variant. Similarly, the L452R mutation of B.1.617.2 (Delta) variant is present in A.2.5. In this study, we observed that these Spike variants bound better to the ACE2 receptor in a temperature-dependent manner. Pseudoviral particles bearing the Spike of Mu were similarly neutralized by plasma from vaccinated individuals than those carrying the Beta (B.1.351) and Delta (B.1.617.2) Spikes. Altogether, our results indicate the importance of measuring critical parameters such as ACE2 interaction, plasma recognition and neutralization ability of each emerging variant.

SARS-CoV-2 Omicron Spike recognition by plasma from individuals receiving BNT162b2 mRNA vaccination with a 16-week interval between doses.

Continuous emergence of SARS-CoV-2 variants of concern (VOCs) is fueling the COVID-19 pandemic. Omicron (B.1.1.529) rapidly spread worldwide. The large number of mutations in its Spike raise concerns about a major antigenic drift that could significantly decrease vaccine efficacy and infection-induced immunity. A long interval between BNT162b2 mRNA doses elicits antibodies that efficiently recognize Spikes from different VOCs. Here, we evaluate the recognition of Omicron Spike by plasma from a cohort of SARS-CoV-2 naive and previously infected individuals who received their BNT162b2 mRNA vaccine 16 weeks apart. Omicron Spike is recognized less efficiently than D614G, Alpha, Beta, Gamma, and Delta Spikes. We compare with plasma activity from participants receiving a short (4 weeks) interval regimen. Plasma from individuals of the long-interval cohort recognize and neutralize better the Omicron Spike compared with those who received a short interval. Whether this difference confers any clinical benefit against Omicron remains unknown.

Addressing Privacy Concerns in Sharing Viral Sequences and Minimum Contextual Data in a Public Repository During the COVID-19 Pandemic.

COVID-19 was declared to be a pandemic in March 2020 by the World Health Organization. Timely sharing of viral genomic sequencing data accompanied by a minimal set of contextual data is essential for informing regional, national, and international public health responses. Such contextual data is also necessary for developing, and improving clinical therapies and vaccines, and enhancing the scientific community's understanding of the SARS-CoV-2 virus. The Canadian COVID-19 Genomics Network (CanCOGeN) was launched in April 2020 to coordinate and upscale existing genomics-based COVID-19 research and surveillance efforts. CanCOGeN is performing large-scale sequencing of both the genomes of SARS-CoV-2 virus samples (VirusSeq) and affected Canadians (HostSeq). This paper addresses the privacy concerns associated with sharing the viral sequence data with a pre-defined set of contextual data describing the sample source and case attribute of the sequence data in the Canadian context. Currently, the viral genome sequences are shared by provincial public health laboratories and their healthcare and academic partners, with the Canadian National Microbiology Laboratory and with publicly accessible databases. However, data sharing delays and the provision of incomplete contextual data often occur because publicly releasing such data triggers privacy and data governance concerns. The CanCOGeN Ethics and Governance Expert Working Group thus has investigated several privacy issues cited by CanCOGeN data providers/stewards. This paper addresses these privacy concerns and offers insights primarily in the Canadian context, although similar privacy considerations also exist in other jurisdictions. We maintain that sharing viral sequencing data and its limited associated contextual data in the public domain generally does not pose insurmountable privacy challenges. However, privacy risks associated with reidentification should be actively monitored due to advancements in reidentification methods and the evolving pandemic landscape. We also argue that during a global health emergency such as COVID-19, privacy should not be used as a blanket measure to prevent such genomic data sharing due to the significant benefits it provides towards public health responses and ongoing research activities.

Outbreak investigation of SARS-CoV-2 transmission in an emergency childcare centre.

OBJECTIVE: To characterize SARS-CoV-2 transmission following a COVID-19 outbreak in an emergency childcare centre (ECCC) in April 2020 in Quebec, Canada. METHODS: The study population consisted of all the children and employees who attended the ECCC as well as household contacts of the confirmed COVID-19 cases. Of the 120 individuals in the study, five cases were confirmed by epidemiological link and 25 were identified as COVID-19 by RT-PCR among which 19 were analyzed by viral whole genome sequencing. Descriptive epidemiology, social network visualization, and phylogenetic analysis were used to characterize viral transmission. RESULTS: Phylogenetic analysis identified two separate introductions of distinct lineages of SARS-CoV-2 and estimated an average effective reproductive number of Re = 1.9 (range 0.9-4.9) with a mean doubling time of 3.2 days (range 2.1-5.2). The first and most prevalent lineage was introduced by two asymptomatic children who were likely infected by their parent, a confirmed COVID-19 case working in a long-term care centre. Among infected household adults, attack rates were significantly higher in mothers than in fathers (risk ratio = 4.5; 95% CI 1.1-18.7). The extent of transmission makes it one of the largest documented outbreaks in a daycare in Canada. CONCLUSION: The analyses carried out showed the probable origin and direction of the transmission of the infection (adult-child, child-adult, and child-child), thus highlighting how asymptomatic children can efficiently transmit SARS-CoV-2.

RéSUMé: OBJECTIF: Caractériser la transmission du SRAS-CoV-2 à la suite d'une éclosion de COVID-19 dans un service de garde d'urgence en milieu scolaire (SGUMS) en avril 2020 au Québec, Canada. MéTHODES: La population à l'étude était composée de tous les enfants et employés ayant fréquenté le SGUMS ainsi que les contacts familiaux des cas confirmés de COVID-19. Sur les 120 personnes à l'étude, cinq cas ont été confirmés par lien épidémiologique et 25 par RT-PCR. Parmi ces derniers, 19 ont été analysés par séquençage viral du génome entier. La caractérisation de la transmission a été réalisée à l'aide d'analyses descriptives et phylogénétiques ainsi que de la visualisation de réseaux sociaux. RéSULTATS: L'analyse phylogénétique a identifié deux introductions de lignées distinctes du SRAS-CoV-2 et un taux de reproduction net Re = 1,9 (étendue 0,9­4,9) avec un temps moyen de doublement de 3,2 jours (étendue 2,1­5,2). La première lignée, et la plus répandue, a été introduite par deux enfants asymptomatiques qui ont probablement été infectés par leur parent, un travailleur de la santé atteint de COVID-19. Dans les noyaux familiaux, les taux d'attaque étaient significativement plus élevés chez les mères que chez les pères (rapport de risque = 4,5 ; IC à 95 % 1,1­18,7). L'ampleur de la transmission en fait de celle-ci la plus importante éclosion documentée dans un service de garde au Canada. CONCLUSION: Cette étude a permis de déterminer l'origine et la direction probables de la transmission de l'infection (adulte-enfant, enfant-adulte et enfant-enfant) et démontrer que les enfants asymptomatiques peuvent transmettre le SRAS-CoV-2.

A small number of early introductions seeded widespread transmission of SARS-CoV-2 in Québec, Canada.

BACKGROUND: Québec was the Canadian province most impacted by COVID-19, with 401,462 cases as of September 24th, 2021, and 11,347 deaths due mostly to a very severe first pandemic wave. In April 2020, we assembled the Coronavirus Sequencing in Québec (CoVSeQ) consortium to sequence SARS-CoV-2 genomes in Québec to track viral introduction events and transmission within the province. METHODS: Using genomic epidemiology, we investigated the arrival of SARS-CoV-2 to Québec. We report 2921 high-quality SARS-CoV-2 genomes in the context of > 12,000 publicly available genomes sampled globally over the first pandemic wave (up to June 1st, 2020). By combining phylogenetic and phylodynamic analyses with epidemiological data, we quantify the number of introduction events into Québec, identify their origins, and characterize the spatiotemporal spread of the virus. RESULTS: Conservatively, we estimated approximately 600 independent introduction events, the majority of which happened from spring break until 2 weeks after the Canadian border closed for non-essential travel. Subsequent mass repatriations did not generate large transmission lineages (> 50 sequenced cases), likely due to mandatory quarantine measures in place at the time. Consistent with common spring break and "snowbird" destinations, most of the introductions were inferred to have originated from Europe via the Americas. Once introduced into Québec, viral lineage sizes were overdispersed, with a few lineages giving rise to most infections. Consistent with founder effects, the earliest lineages to arrive tended to spread most successfully. Fewer than 100 viral introductions arrived during spring break, of which 7-12 led to the largest transmission lineages of the first wave (accounting for 52-75% of all sequenced infections). These successful transmission lineages dispersed widely across the province. Transmission lineage size was greatly reduced after March 11th, when a quarantine order for returning travellers was enacted. While this suggests the effectiveness of early public health measures, the biggest transmission lineages had already been ignited prior to this order. CONCLUSIONS: Combined, our results reinforce how, in the absence of tight travel restrictions or quarantine measures, fewer than 100 viral introductions in a week can ensure the establishment of extended transmission chains.

Impact of temperature on the affinity of SARS-CoV-2 Spike for ACE2.

The seasonal nature in the outbreaks of respiratory viral infections with increased transmission during low temperatures has been well established. The current COVID-19 pandemic makes no exception, and temperature has been suggested to play a role on the viability and transmissibility of SARS-CoV-2. The receptor binding domain (RBD) of the Spike glycoprotein binds to the angiotensin-converting enzyme 2 (ACE2) to initiate viral fusion. Studying the effect of temperature on the receptor-Spike interaction, we observed a significant and stepwise increase in RBD-ACE2 affinity at low temperatures, resulting in slower dissociation kinetics. This translated into enhanced interaction of the full Spike to ACE2 receptor and higher viral attachment at low temperatures. Interestingly, the RBD N501Y mutation, present in emerging variants of concern (VOCs) that are fueling the pandemic worldwide, bypassed this requirement. This data suggests that the acquisition of N501Y reflects an adaptation to warmer climates, a hypothesis that remains to be tested.