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.

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.

SARS-CoV-2 rapidly evolves lineage-specific phenotypic differences when passaged repeatedly in immune-naïve mice.

The persistence of SARS-CoV-2 despite the development of vaccines and a degree of herd immunity is partly due to viral evolution reducing vaccine and treatment efficacy. Serial infections of wild-type (WT) SARS-CoV-2 in Balb/c mice yield mouse-adapted strains with greater infectivity and mortality. We investigate if passaging unmodified B.1.351 (Beta) and B.1.617.2 (Delta) 20 times in K18-ACE2 mice, expressing the human ACE2 receptor, in a BSL-3 laboratory without selective pressures, drives human health-relevant evolution and if evolution is lineage-dependent. Late-passage virus causes more severe disease, at organism and lung tissue scales, with late-passage Delta demonstrating antibody resistance and interferon suppression. This resistance co-occurs with a de novo spike S371F mutation, linked with both traits. S371F, an Omicron-characteristic mutation, is co-inherited at times with spike E1182G per Nanopore sequencing, existing in different within-sample viral variants at others. Both S371F and E1182G are linked to mammalian GOLGA7 and ZDHHC5 interactions, which mediate viral-cell entry and antiviral response. This study demonstrates SARS-CoV-2's tendency to evolve with phenotypic consequences, its evolution varying by lineage, and suggests non-dominant quasi-species contribution.

Genetic variation in the Y chromosome and sex-biased DNA methylation in somatic cells in the mouse.

Several lines of evidence suggest that the presence of the Y chromosome influences DNA methylation of autosomal loci. To better understand the impact of the Y chromosome on autosomal DNA methylation patterns and its contribution to sex bias in methylation, we identified Y chromosome dependent differentially methylated regions (yDMRs) using whole-genome bisulfite sequencing methylation data from livers of mice with different combinations of sex-chromosome complement and gonadal sex. Nearly 90% of the autosomal yDMRs mapped to transposable elements (TEs) and most of them had lower methylation in XY compared to XX or XO mice. Follow-up analyses of four reporter autosomal yDMRs showed that Y-dependent methylation levels were consistent across most somatic tissues but varied in strains with different origins of the Y chromosome, suggesting that genetic variation in the Y chromosome influenced methylation levels of autosomal regions. Mice lacking the q-arm of the Y chromosome (B6.NPYq-2) as well as mice with a loss-of-function mutation in Kdm5d showed no differences in methylation levels compared to wild type mice. In conclusion, the Y-linked modifier of TE methylation is likely to reside on the short arm of Y chromosome and further studies are required to identify this gene.

More than a 'Hundred Days War': Persistent SARS-CoV-2 infection in a patient with ANCA-associated vasculitis.

BACKGROUND: Few reports exist on the characteristics and outcomes of persistent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in immunocompromised hosts. METHODS: A 49-year-old patient with granulomatosis with polyangiitis (GPA) and a renal transplant experienced multiple hospitalizations for coronavirus disease 2019 (COVID-19) pneumonia and relapses between October 2020 and February 2021. Careful chart review of medical history, hospitalizations, and microbiological testing including SARS-CoV-2 cycle threshold values, therapies, and imaging was undertaken. SARS-CoV-2 genome sequencing was performed in five viral samples to distinguish persistent infection from re-infection with a different strain. RESULTS: Sequencing confirmed that all samples tested were from the same viral lineage, indicating a long-term, persistent infection rather than re-infection with a new strain. The patient ultimately stabilized after two courses of remdesivir plus dexamethasone, replacement intravenous immunoglobulin, and bamlanivimab. Rituximab maintenance therapy for vasculitis remains on hold. CONCLUSIONS: SARS-CoV-2 may persist for several months in immunocompromised hosts and may go unrecognized as an ongoing active infection. More studies are needed to determine how to optimize COVID-19 treatment in this vulnerable population.

HISTORIQUE: Il existe peu de rapports sur les caractéristiques et les issues de l'infection par le coronavirus 2 du syndrome respiratoire aigu sévère (SRAS-CoV-2) chez les hôtes immunodéprimés. MÉTHODOLOGIE: UNE PATIENTE de 49 ans receveuse d'une transplantation rénale atteinte d'une granulomatose avec polyangéite a été hospitalisée à de multiples reprises à cause d'une pneumonie à maladie à coronavirus 2019 (COVID-19) et de récidives entre octobre 2020 et février 2021. Les chercheurs ont exécuté une analyse attentive du dossier pour connaître l'histoire médicale de la patiente, les hospitalisations et les tests microbiologiques effectués, y compris les valeurs seuils du cycle du SRAS-CoV-2, les traitements et les techniques d'imagerie. Ils ont procédé au séquençage du génome du SRAS-CoV-2 dans cinq prélèvements viraux pour distinguer l'infection persistante de la réinfection par une souche différente. RÉSULTATS : Le séquençage a confirmé que tous les prélèvements effectués provenaient de la même lignée virale, ce qui détermine une infection persistante prolongée plutôt qu'une réinfection par une nouvelle souche. L'état de la patiente a fini par se stabiliser après deux traitements au remdésivir combiné à de la dexaméthasone, une thérapie de substitution par immunoglobuline intraveineuse et du bamlanivimab. Un traitement d'entretien de la vasculite au rituximab demeure en suspens. CONCLUSIONS: Le SRAS-CoV-2 peut persister plusieurs mois chez les hôtes immunodéprimés, et un état d'infection active continue peut passer inaperçu. Plus d'études devront être réalisées pour déterminer le moyen d'optimiser le traitement de la COVID-19 dans cette population vulnérable.

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.

Distinct roles of androgen receptor, estrogen receptor alpha, and BCL6 in the establishment of sex-biased DNA methylation in mouse liver.

Sexual dimorphism in gene regulation, including DNA methylation, is the main driver of sexual dimorphism in phenotypes. However, the questions of how and when sex shapes DNA methylation remain unresolved. Recently, using mice with different combinations of genetic and phenotypic sex, we identified sex-associated differentially methylated regions (sDMRs) that depended on the sex phenotype. Focusing on a panel of validated sex-phenotype dependent male- and female-biased sDMRs, we tested the developmental dynamics of sex bias in liver methylation and the impacts of mutations in the androgen receptor, estrogen receptor alpha, or the transcriptional repressor Bcl6 gene. True hermaphrodites that carry both unilateral ovaries and contralateral testes were also tested. Our data show that sex bias in methylation either coincides with or follows sex bias in the expression of sDMR-proximal genes, suggesting that sex bias in gene expression may be required for demethylation at certain sDMRs. Global ablation of AR, ESR1, or a liver-specific loss of BCL6, all alter sDMR methylation, whereas presence of both an ovary and a testis delays the establishment of male-type methylation levels in hermaphrodites. Moreover, the Bcl6-LKO shows dissociation between expression and methylation, suggesting a distinct role of BCL6 in demethylation of intragenic sDMRs.

Ultrafast functional profiling of RNA-seq data for nonmodel organisms.

Computational time and cost remain a major bottleneck for RNA-seq data analysis of nonmodel organisms without reference genomes. To address this challenge, we have developed Seq2Fun, a novel, all-in-one, ultrafast tool to directly perform functional quantification of RNA-seq reads without transcriptome de novo assembly. The pipeline starts with raw read quality control: sequencing error correction, removing poly(A) tails, and joining overlapped paired-end reads. It then conducts a DNA-to-protein search by translating each read into all possible amino acid fragments and subsequently identifies possible homologous sequences in a well-curated protein database. Finally, the pipeline generates several informative outputs including gene abundance tables, pathway and species hit tables, an HTML report to visualize the results, and an output of clean reads annotated with mapped genes ready for downstream analysis. Seq2Fun does not have any intermediate steps of file writing and loading, making I/O very efficient. Seq2Fun is written in C++ and can run on a personal computer with a limited number of CPUs and memory. It can process >2,000,000 reads/min and is >120 times faster than conventional workflows based on de novo assembly, while maintaining high accuracy in our various test data sets.

Sex Chromosomes and Sex Phenotype Contribute to Biased DNA Methylation in Mouse Liver.

Sex biases in the genome-wide distribution of DNA methylation and gene expression levels are some of the manifestations of sexual dimorphism in mammals. To advance our understanding of the mechanisms that contribute to sex biases in DNA methylation and gene expression, we conducted whole genome bisulfite sequencing (WGBS) as well as RNA-seq on liver samples from mice with different combinations of sex phenotype and sex-chromosome complement. We compared groups of animals with different sex phenotypes, but the same genetic sexes, and vice versa, same sex phenotypes, but different sex-chromosome complements. We also compared sex-biased DNA methylation in mouse and human livers. Our data show that sex phenotype, X-chromosome dosage, and the presence of Y chromosome shape the differences in DNA methylation between males and females. We also demonstrate that sex bias in autosomal methylation is associated with sex bias in gene expression, whereas X-chromosome dosage-dependent methylation differences are not, as expected for a dosage-compensation mechanism. Furthermore, we find partial conservation between the repertoires of mouse and human genes that are associated with sex-biased methylation, an indication that gene function is likely to be an important factor in this phenomenon.

Regulatory interaction between the ZPBP2-ORMDL3/Zpbp2-Ormdl3 region and the circadian clock.

Genome-wide association study (GWAS) loci for several immunity-mediated diseases (early onset asthma, inflammatory bowel disease (IBD), primary biliary cholangitis, and rheumatoid arthritis) map to chromosomal region 17q12-q21. The predominant view is that association between 17q12-q21 alleles and increased risk of developing asthma or IBD is due to regulatory variants. ORM sphingolipid biosynthesis regulator (ORMDL3) residing in this region is the most promising gene candidate for explaining association with disease. However, the relationship between 17q12-q21 alleles and disease is complex suggesting contributions from other factors, such as trans-acting genetic and environmental modifiers or circadian rhythms. Circadian rhythms regulate expression levels of thousands of genes and their dysregulation is implicated in the etiology of several common chronic inflammatory diseases. However, their role in the regulation of the 17q12-q21 genes has not been investigated. Moreover, the core clock gene nuclear receptor subfamily 1, group D, member 1 (NR1D1) resides about 200 kb distal to the GWAS region. We hypothesized that circadian rhythms influenced gene expression levels in 17q12-q21 region and conversely, regulatory elements in this region influenced transcription of the core clock gene NR1D1 in cis. To test these hypotheses, we examined the diurnal expression profiles of zona pellucida binding protein 2 (ZPBP2/Zpbp2), gasdermin B (GSDMB), and ORMDL3/Ormdl3 in human and mouse tissues and analyzed the impact of genetic variation in the ZPBP2/Zpbp2 region on NR1D1/Nr1d1 expression. We found that Ormdl3 and Zpbp2 were controlled by the circadian clock in a tissue-specific fashion. We also report that deletion of the Zpbp2 region altered the expression profile of Nr1d1 in lungs and ileum in a time-dependent manner. In liver, the deletion was associated with enhanced expression of Ormdl3. We provide the first evidence that disease-associated genes Zpbp2 and Ormdl3 are regulated by circadian rhythms and the Zpbp2 region influences expression of the core clock gene Nr1d1.

The nitrogen responsive transcriptome in potato (Solanum tuberosum L.) reveals significant gene regulatory motifs.

Nitrogen (N) is the most important nutrient for the growth of potato (Solanum tuberosum L.). Foliar gene expression in potato plants with and without N supplementation at 180 kg N ha(-1) was compared at mid-season. Genes with consistent differences in foliar expression due to N supplementation over three cultivars and two developmental time points were examined. In total, thirty genes were found to be over-expressed and nine genes were found to be under-expressed with supplemented N. Functional relationships between over-expressed genes were found. The main metabolic pathway represented among differentially expressed genes was amino acid metabolism. The 1000 bp upstream flanking regions of the differentially expressed genes were analysed and nine overrepresented motifs were found using three motif discovery algorithms (Seeder, Weeder and MEME). These results point to coordinated gene regulation at the transcriptional level controlling steady state potato responses to N sufficiency.