ABSTRACT
Inborn errors of human IFN-γ-dependent macrophagic immunity underlie mycobacterial diseases, whereas inborn errors of IFN-α/ß-dependent intrinsic immunity underlie viral diseases. Both types of IFNs induce the transcription factor IRF1. We describe unrelated children with inherited complete IRF1 deficiency and early-onset, multiple, life-threatening diseases caused by weakly virulent mycobacteria and related intramacrophagic pathogens. These children have no history of severe viral disease, despite exposure to many viruses, including SARS-CoV-2, which is life-threatening in individuals with impaired IFN-α/ß immunity. In leukocytes or fibroblasts stimulated in vitro, IRF1-dependent responses to IFN-γ are, both quantitatively and qualitatively, much stronger than those to IFN-α/ß. Moreover, IRF1-deficient mononuclear phagocytes do not control mycobacteria and related pathogens normally when stimulated with IFN-γ. By contrast, IFN-α/ß-dependent intrinsic immunity to nine viruses, including SARS-CoV-2, is almost normal in IRF1-deficient fibroblasts. Human IRF1 is essential for IFN-γ-dependent macrophagic immunity to mycobacteria, but largely redundant for IFN-α/ß-dependent antiviral immunity.
Subject(s)
COVID-19 , Mycobacterium , Child , Humans , Interferon-gamma , SARS-CoV-2 , Interferon-alpha , Interferon Regulatory Factor-1ABSTRACT
We study a patient with the human papilloma virus (HPV)-2-driven "tree-man" phenotype and two relatives with unusually severe HPV4-driven warts. The giant horns form an HPV-2-driven multifocal benign epithelial tumor overexpressing viral oncogenes in the epidermis basal layer. The patients are unexpectedly homozygous for a private CD28 variant. They have no detectable CD28 on their T cells, with the exception of a small contingent of revertant memory CD4+ T cells. T cell development is barely affected, and T cells respond to CD3 and CD2, but not CD28, costimulation. Although the patients do not display HPV-2- and HPV-4-reactive CD4+ T cells in vitro, they make antibodies specific for both viruses in vivo. CD28-deficient mice are susceptible to cutaneous infections with the mouse papillomavirus MmuPV1. The control of HPV-2 and HPV-4 in keratinocytes is dependent on the T cell CD28 co-activation pathway. Surprisingly, human CD28-dependent T cell responses are largely redundant for protective immunity.
Subject(s)
CD28 Antigens/deficiency , Inheritance Patterns/genetics , Papillomaviridae/physiology , Skin/virology , T-Lymphocytes/immunology , Adult , Amino Acid Sequence , Animals , Base Sequence , CD28 Antigens/genetics , CD28 Antigens/metabolism , CD4-Positive T-Lymphocytes/immunology , Child , Endopeptidases/metabolism , Female , Genes, Recessive , HEK293 Cells , Homozygote , Humans , Immunity, Humoral , Immunologic Memory , Jurkat Cells , Keratinocytes/pathology , Male , Mice, Inbred C57BL , Oncogenes , Papilloma/pathology , Papilloma/virology , Pedigree , Protein Sorting Signals , RNA, Messenger/genetics , RNA, Messenger/metabolismABSTRACT
Inborn errors of human interferon gamma (IFN-γ) immunity underlie mycobacterial disease. We report a patient with mycobacterial disease due to inherited deficiency of the transcription factor T-bet. The patient has extremely low counts of circulating Mycobacterium-reactive natural killer (NK), invariant NKT (iNKT), mucosal-associated invariant T (MAIT), and Vδ2+ γδ T lymphocytes, and of Mycobacterium-non reactive classic TH1 lymphocytes, with the residual populations of these cells also producing abnormally small amounts of IFN-γ. Other lymphocyte subsets develop normally but produce low levels of IFN-γ, with the exception of CD8+ αß T and non-classic CD4+ αß TH1∗ lymphocytes, which produce IFN-γ normally in response to mycobacterial antigens. Human T-bet deficiency thus underlies mycobacterial disease by preventing the development of innate (NK) and innate-like adaptive lymphocytes (iNKT, MAIT, and Vδ2+ γδ T cells) and IFN-γ production by them, with mycobacterium-specific, IFN-γ-producing, purely adaptive CD8+ αß T, and CD4+ αß TH1∗ cells unable to compensate for this deficit.
Subject(s)
Adaptive Immunity , Immunity, Innate , Interferon-gamma/immunology , Mycobacterium/immunology , T-Box Domain Proteins/metabolism , Amino Acid Sequence , Base Sequence , Cell Lineage , Child, Preschool , Chromatin/metabolism , CpG Islands/genetics , DNA Methylation/genetics , Dendritic Cells/metabolism , Epigenesis, Genetic , Female , Homozygote , Humans , INDEL Mutation/genetics , Infant , Interferon-gamma/metabolism , Killer Cells, Natural/cytology , Killer Cells, Natural/metabolism , Loss of Function Mutation/genetics , Male , Mycobacterium Infections/genetics , Mycobacterium Infections/immunology , Mycobacterium Infections/microbiology , Pedigree , T-Box Domain Proteins/chemistry , T-Box Domain Proteins/deficiency , T-Box Domain Proteins/genetics , T-Lymphocytes, Helper-Inducer/immunology , Transcriptome/geneticsABSTRACT
Cerebral white matter hyperintensities on MRI are markers of cerebral small vessel disease, a major risk factor for dementia and stroke. Despite the successful identification of multiple genetic variants associated with this highly heritable condition, its genetic architecture remains incompletely understood. More specifically, the role of DNA methylation has received little attention. We investigated the association between white matter hyperintensity burden and DNA methylation in blood at â¼450 000 cytosine-phosphate-guanine (CpG) sites in 9732 middle-aged to older adults from 14 community-based studies. Single CpG and region-based association analyses were carried out. Functional annotation and integrative cross-omics analyses were performed to identify novel genes underlying the relationship between DNA methylation and white matter hyperintensities. We identified 12 single CpG and 46 region-based DNA methylation associations with white matter hyperintensity burden. Our top discovery single CpG, cg24202936 (P = 7.6 × 10-8), was associated with F2 expression in blood (P = 6.4 × 10-5) and co-localized with FOLH1 expression in brain (posterior probability = 0.75). Our top differentially methylated regions were in PRMT1 and in CCDC144NL-AS1, which were also represented in single CpG associations (cg17417856 and cg06809326, respectively). Through Mendelian randomization analyses cg06809326 was putatively associated with white matter hyperintensity burden (P = 0.03) and expression of CCDC144NL-AS1 possibly mediated this association. Differentially methylated region analysis, joint epigenetic association analysis and multi-omics co-localization analysis consistently identified a role of DNA methylation near SH3PXD2A, a locus previously identified in genome-wide association studies of white matter hyperintensities. Gene set enrichment analyses revealed functions of the identified DNA methylation loci in the blood-brain barrier and in the immune response. Integrative cross-omics analysis identified 19 key regulatory genes in two networks related to extracellular matrix organization, and lipid and lipoprotein metabolism. A drug-repositioning analysis indicated antihyperlipidaemic agents, more specifically peroxisome proliferator-activated receptor-alpha, as possible target drugs for white matter hyperintensities. Our epigenome-wide association study and integrative cross-omics analyses implicate novel genes influencing white matter hyperintensity burden, which converged on pathways related to the immune response and to a compromised blood-brain barrier possibly due to disrupted cell-cell and cell-extracellular matrix interactions. The results also suggest that antihyperlipidaemic therapy may contribute to lowering risk for white matter hyperintensities possibly through protection against blood-brain barrier disruption.
Subject(s)
White Matter , Middle Aged , Humans , Aged , White Matter/diagnostic imaging , Genome-Wide Association Study/methods , Brain/diagnostic imaging , DNA Methylation/genetics , Magnetic Resonance Imaging , Epigenesis, Genetic , Protein-Arginine N-Methyltransferases , Repressor ProteinsABSTRACT
Recent reports have identified differences in the mutational spectra across human populations. Although some of these reports have been replicated in other cohorts, most have been reported only in the 1000 Genomes Project (1kGP) data. While investigating an intriguing putative population stratification within the Japanese population, we identified a previously unreported batch effect leading to spurious mutation calls in the 1kGP data and to the apparent population stratification. Because the 1kGP data are used extensively, we find that the batch effects also lead to incorrect imputation by leading imputation servers and a small number of suspicious GWAS associations. Lower quality data from the early phases of the 1kGP thus continue to contaminate modern studies in hidden ways. It may be time to retire or upgrade such legacy sequencing data.
Subject(s)
Human Genome Project , Artifacts , Genome-Wide Association Study , Humans , Japan , MutationABSTRACT
RATIONALE: Familial recurrence studies provide strong evidence for a genetic component to the predisposition to sporadic, nonsyndromic Tetralogy of Fallot (TOF), the most common cyanotic congenital heart disease phenotype. Rare genetic variants have been identified as important contributors to the risk of congenital heart disease, but relatively small numbers of TOF cases have been studied to date. OBJECTIVE: We used whole exome sequencing to assess the prevalence of unique, deleterious variants in the largest cohort of nonsyndromic TOF patients reported to date. METHODS AND RESULTS: Eight hundred twenty-nine TOF patients underwent whole exome sequencing. The presence of unique, deleterious variants was determined; defined by their absence in the Genome Aggregation Database and a scaled combined annotation-dependent depletion score of ≥20. The clustering of variants in 2 genes, NOTCH1 and FLT4, surpassed thresholds for genome-wide significance (assigned as P<5×10-8) after correction for multiple comparisons. NOTCH1 was most frequently found to harbor unique, deleterious variants. Thirty-one changes were observed in 37 probands (4.5%; 95% CI, 3.2%-6.1%) and included 7 loss-of-function variants 22 missense variants and 2 in-frame indels. Sanger sequencing of the unaffected parents of 7 cases identified 5 de novo variants. Three NOTCH1 variants (p.G200R, p.C607Y, and p.N1875S) were subjected to functional evaluation, and 2 showed a reduction in Jagged1-induced NOTCH signaling. FLT4 variants were found in 2.4% (95% CI, 1.6%-3.8%) of TOF patients, with 21 patients harboring 22 unique, deleterious variants. The variants identified were distinct to those that cause the congenital lymphoedema syndrome Milroy disease. In addition to NOTCH1, FLT4 and the well-established TOF gene, TBX1, we identified potential association with variants in several other candidates, including RYR1, ZFPM1, CAMTA2, DLX6, and PCM1. CONCLUSIONS: The NOTCH1 locus is the most frequent site of genetic variants predisposing to nonsyndromic TOF, followed by FLT4. Together, variants in these genes are found in almost 7% of TOF patients.
Subject(s)
Exome , Mutation Rate , Tetralogy of Fallot/genetics , Autoantigens/genetics , Calcium-Binding Proteins/genetics , Cell Cycle Proteins/genetics , Homeodomain Proteins/genetics , Humans , Loss of Function Mutation , Mutation, Missense , Nuclear Proteins/genetics , Receptor, Notch1/genetics , Trans-Activators/genetics , Transcription Factors/genetics , Vascular Endothelial Growth Factor Receptor-3/geneticsABSTRACT
Whole genome sequencing (WGS) at high-depth (30X) allows the accurate discovery of variants in the coding and non-coding DNA regions and helps elucidate the genetic underpinnings of human health and diseases. Yet, due to the prohibitive cost of high-depth WGS, most large-scale genetic association studies use genotyping arrays or high-depth whole exome sequencing (WES). Here we propose a cost-effective method which we call "Whole Exome Genome Sequencing" (WEGS), that combines low-depth WGS and high-depth WES with up to 8 samples pooled and sequenced simultaneously (multiplexed). We experimentally assess the performance of WEGS with four different depth of coverage and sample multiplexing configurations. We show that the optimal WEGS configurations are 1.7-2.0 times cheaper than standard WES (no-plexing), 1.8-2.1 times cheaper than high-depth WGS, reach similar recall and precision rates in detecting coding variants as WES, and capture more population-specific variants in the rest of the genome that are difficult to recover when using genotype imputation methods. We apply WEGS to 862 patients with peripheral artery disease and show that it directly assesses more known disease-associated variants than a typical genotyping array and thousands of non-imputable variants per disease-associated locus.
ABSTRACT
The basal breast cancer subtype is enriched for triple-negative breast cancer (TNBC) and displays consistent large chromosomal deletions. Here, we characterize evolution and maintenance of chromosome 4p (chr4p) loss in basal breast cancer. Analysis of The Cancer Genome Atlas data shows recurrent deletion of chr4p in basal breast cancer. Phylogenetic analysis of a panel of 23 primary tumor/patient-derived xenograft basal breast cancers reveals early evolution of chr4p deletion. Mechanistically we show that chr4p loss is associated with enhanced proliferation. Gene function studies identify an unknown gene, C4orf19, within chr4p, which suppresses proliferation when overexpressed-a member of the PDCD10-GCKIII kinase module we name PGCKA1. Genome-wide pooled overexpression screens using a barcoded library of human open reading frames identify chromosomal regions, including chr4p, that suppress proliferation when overexpressed in a context-dependent manner, implicating network interactions. Together, these results shed light on the early emergence of complex aneuploid karyotypes involving chr4p and adaptive landscapes shaping breast cancer genomes.
Subject(s)
Breast Neoplasms , Gene Regulatory Networks , Humans , Female , Breast Neoplasms/genetics , Breast Neoplasms/pathology , Animals , Mice , Chromosomes, Human, Pair 4/genetics , Cell Proliferation/genetics , Chromosome Aberrations , Cell Line, Tumor , Triple Negative Breast Neoplasms/genetics , Triple Negative Breast Neoplasms/pathologyABSTRACT
Patients with inherited CARMIL2 or CD28 deficiency have defective T cell CD28 signaling, but their immunological and clinical phenotypes remain largely unknown. We show that only one of three CARMIL2 isoforms is produced and functional across leukocyte subsets. Tested mutant CARMIL2 alleles from 89 patients and 52 families impair canonical NF-κB but not AP-1 and NFAT activation in T cells stimulated via CD28. Like CD28-deficient patients, CARMIL2-deficient patients display recalcitrant warts and low blood counts of CD4+ and CD8+ memory T cells and CD4+ TREGs. Unlike CD28-deficient patients, they have low counts of NK cells and memory B cells, and their antibody responses are weak. CARMIL2 deficiency is fully penetrant by the age of 10 yr and is characterized by numerous infections, EBV+ smooth muscle tumors, and mucocutaneous inflammation, including inflammatory bowel disease. Patients with somatic reversions of a mutant allele in CD4+ T cells have milder phenotypes. Our study suggests that CARMIL2 governs immunological pathways beyond CD28.
Subject(s)
CD28 Antigens , Microfilament Proteins , Humans , CD28 Antigens/metabolism , Microfilament Proteins/genetics , Mutation/genetics , Phenotype , CD4-Positive T-LymphocytesABSTRACT
Perivascular space (PVS) burden is an emerging, poorly understood, magnetic resonance imaging marker of cerebral small vessel disease, a leading cause of stroke and dementia. Genome-wide association studies in up to 40,095 participants (18 population-based cohorts, 66.3 ± 8.6 yr, 96.9% European ancestry) revealed 24 genome-wide significant PVS risk loci, mainly in the white matter. These were associated with white matter PVS already in young adults (N = 1,748; 22.1 ± 2.3 yr) and were enriched in early-onset leukodystrophy genes and genes expressed in fetal brain endothelial cells, suggesting early-life mechanisms. In total, 53% of white matter PVS risk loci showed nominally significant associations (27% after multiple-testing correction) in a Japanese population-based cohort (N = 2,862; 68.3 ± 5.3 yr). Mendelian randomization supported causal associations of high blood pressure with basal ganglia and hippocampal PVS, and of basal ganglia PVS and hippocampal PVS with stroke, accounting for blood pressure. Our findings provide insight into the biology of PVS and cerebral small vessel disease, pointing to pathways involving extracellular matrix, membrane transport and developmental processes, and the potential for genetically informed prioritization of drug targets.
Subject(s)
Cerebral Small Vessel Diseases , Stroke , Humans , Endothelial Cells/pathology , Genome-Wide Association Study , Cerebral Small Vessel Diseases/diagnostic imaging , Cerebral Small Vessel Diseases/genetics , Cerebral Small Vessel Diseases/complications , Magnetic Resonance Imaging/methods , GenomicsABSTRACT
The study of the genetic component of early-onset diseases requires investigation into parental genetic effects, particularly those mediated by the mother who can influence the offspring's risk of disease through the effects of her genes acting directly on the intrauterine milieu or indirectly through maternal-gene child-gene interaction effects. An important source of bias that can arise in feto-maternal association testing is the possibility of confounding due to mating asymmetry (MA). However, there is little information on the levels of MA present in human populations and the impact on maternal association testing. In this study, we developed a novel approach to measuring MA and, using HapMap mate-pairs of European and African descent, carried out a genome-wide investigation and characterization of MA. We further investigated the impact of observed levels of MA on maternal association tests through simulation experiments. For the first time, we showed that non-negligible levels of MA are detected in human populations, such that subtle genotype frequency differences between individuals mating in the population are sufficient to induce spurious maternal genotype associations. Though the underlying mechanisms driving the asymmetry within these populations remain elusive, our findings provide consequential evidence for the occurrence of MA in humans and highlight the importance of controlling for MA in maternal association testing.
Subject(s)
Genome-Wide Association Study , Molecular Epidemiology/methods , Africa , Europe , Female , Genetic Predisposition to Disease , Genetics, Population , Genome , Genotype , HapMap Project , Humans , Male , Models, Genetic , Models, Statistical , Sample SizeABSTRACT
Pediatric glioblastomas (GBM) including diffuse intrinsic pontine gliomas (DIPG) are devastating brain tumors with no effective therapy. Here, we investigated clinical and biological impacts of histone H3.3 mutations. Forty-two DIPGs were tested for H3.3 mutations. Wild-type versus mutated (K27M-H3.3) subgroups were compared for HIST1H3B, IDH, ATRX and TP53 mutations, copy number alterations and clinical outcome. K27M-H3.3 occurred in 71 %, TP53 mutations in 77 % and ATRX mutations in 9 % of DIPGs. ATRX mutations were more frequent in older children (p < 0.0001). No G34V/R-H3.3, IDH1/2 or H3.1 mutations were identified. K27M-H3.3 DIPGs showed specific copy number changes, including all gains/amplifications of PDGFRA and MYC/PVT1 loci. Notably, all long-term survivors were H3.3 wild type and this group of patients had better overall survival. K27M-H3.3 mutation defines clinically and biologically distinct subgroups and is prevalent in DIPG, which will impact future therapeutic trial design. K27M- and G34V-H3.3 have location-based incidence (brainstem/cortex) and potentially play distinct roles in pediatric GBM pathogenesis. K27M-H3.3 is universally associated with short survival in DIPG, while patients wild-type for H3.3 show improved survival. Based on prognostic and therapeutic implications, our findings argue for H3.3-mutation testing at diagnosis, which should be rapidly integrated into the clinical decision-making algorithm, particularly in atypical DIPG.
Subject(s)
Brain Stem Neoplasms/genetics , Glioma/genetics , Histones/genetics , Pons/pathology , Adolescent , Brain Stem Neoplasms/mortality , Brain Stem Neoplasms/pathology , Child , Child, Preschool , Female , Gene Expression Profiling , Glioma/mortality , Glioma/pathology , Humans , Infant , Male , Mutation , Prognosis , Survival RateABSTRACT
Despite the availability of highly efficacious vaccines, coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lacks effective drug treatment, which results in a high rate of mortality. To address this therapeutic shortcoming, we applied a systems biology approach to the study of patients hospitalized with severe COVID. We show that, at the time of hospital admission, patients who were equivalent on the clinical ordinal scale displayed significant differential monocyte epigenetic and transcriptomic attributes between those who would survive and those who would succumb to COVID-19. We identified messenger RNA metabolism, RNA splicing, and interferon signaling pathways as key host responses overactivated by patients who would not survive. Those pathways are prime drug targets to reduce mortality of critically ill patients with COVID-19, leading us to identify tacrolimus, zotatifin, and nintedanib as three strong candidates for treatment of severely ill patients at the time of hospital admission.
Subject(s)
COVID-19 Drug Treatment , Antiviral Agents/pharmacology , Antiviral Agents/therapeutic use , Humans , SARS-CoV-2 , Systems BiologyABSTRACT
The precisionFDA Truth Challenge V2 aimed to assess the state of the art of variant calling in challenging genomic regions. Starting with FASTQs, 20 challenge participants applied their variant-calling pipelines and submitted 64 variant call sets for one or more sequencing technologies (Illumina, PacBio HiFi, and Oxford Nanopore Technologies). Submissions were evaluated following best practices for benchmarking small variants with updated Genome in a Bottle benchmark sets and genome stratifications. Challenge submissions included numerous innovative methods, with graph-based and machine learning methods scoring best for short-read and long-read datasets, respectively. With machine learning approaches, combining multiple sequencing technologies performed particularly well. Recent developments in sequencing and variant calling have enabled benchmarking variants in challenging genomic regions, paving the way for the identification of previously unknown clinically relevant variants.
ABSTRACT
The kidney and upper urinary tract develop through reciprocal interactions between the ureteric bud and the surrounding mesenchyme. Ureteric bud branching forms the arborized collecting duct system of the kidney, while ureteric tips promote nephron formation from dedicated progenitor cells. While nephron progenitor cells are relatively well characterized, the origin of ureteric bud progenitors has received little attention so far. It is well established that the ureteric bud is induced from the nephric duct, an epithelial duct derived from the intermediate mesoderm of the embryo. However, the cell state transitions underlying the progression from intermediate mesoderm to nephric duct and ureteric bud remain unknown. Here we show that nephric duct morphogenesis results from the coordinated organization of four major progenitor cell populations. Using single cell RNA-seq and Cluster RNA-seq, we show that these progenitors emerge in time and space according to a stereotypical pattern. We identify the transcription factors Tfap2a/b and Gata3 as critical coordinators of this progenitor cell progression. This study provides a better understanding of the cellular origin of the renal collecting duct system and associated urinary tract developmental diseases, which may inform guided differentiation of functional kidney tissue.
Subject(s)
Nephrons/embryology , Organogenesis/genetics , Stem Cells/physiology , Animals , Cell Differentiation/genetics , Embryo, Mammalian , Female , GATA3 Transcription Factor/metabolism , Gene Expression Regulation, Developmental , Male , Mice , Mice, Transgenic , Models, Animal , RNA-Seq , Single-Cell Analysis , Transcription Factor AP-2/metabolismABSTRACT
BACKGROUND: Coronavirus disease 2019 (COVID-19) primarily affects the lungs, but evidence of systemic disease with multi-organ involvement is emerging. Here, we developed a blood test to broadly quantify cell-, tissue-, and organ-specific injury due to COVID-19. METHODS: Our test leverages genome-wide methylation profiling of circulating cell-free DNA in plasma. We assessed the utility of this test to identify subjects with severe disease in two independent, longitudinal cohorts of hospitalized patients. Cell-free DNA profiling was performed on 104 plasma samples from 33 COVID-19 patients and compared to samples from patients with other viral infections and healthy controls. FINDINGS: We found evidence of injury to the lung and liver and involvement of red blood cell progenitors associated with severe COVID-19. The concentration of cell-free DNA correlated with the World Health Organization (WHO) ordinal scale for disease progression and was significantly increased in patients requiring intubation. CONCLUSIONS: This study points to the utility of cell-free DNA as an analyte to monitor and study COVID-19. FUNDING: This work was supported by NIH grants 1DP2AI138242 (to I.D.V.), R01AI146165 (to I.D.V., M.P.C., F.M.M., and J.R.), 1R01AI151059 (to I.D.V.), K08-CA230156 (to W.G.), and R33-AI129455 to C.Y.C., a Synergy award from the Rainin Foundation (to I.D.V.), a SARS-CoV-2 seed grant at Cornell (to I.D.V.), a National Sciences and Engineering Research Council of Canada fellowship PGS-D3 (to A.P.C.), and a Burroughs-Wellcome CAMS Award (to W.G.). D.C.V. is supported by a Fonds de la Recherche en Sante du Quebec Clinical Research Scholar Junior 2 award. C.Y.C. is supported by the California Initiative to Advance Precision Medicine, and the Charles and Helen Schwab Foundation.
Subject(s)
COVID-19 , Cell-Free Nucleic Acids , Virus Diseases , Humans , Methylation , SARS-CoV-2/geneticsABSTRACT
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.
Subject(s)
COVID-19/transmission , COVID-19/epidemiology , COVID-19/virology , Canada/epidemiology , Europe/epidemiology , Genome, Viral , Humans , Molecular Epidemiology , Pandemics , Phylogeny , Public Health , Quebec/epidemiology , SARS-CoV-2/genetics , SARS-CoV-2/isolation & purification , TravelABSTRACT
We studied a child with severe viral, bacterial, fungal, and parasitic diseases, who was homozygous for a loss-of-function mutation of REL, encoding c-Rel, which is selectively expressed in lymphoid and myeloid cells. The patient had low frequencies of NK, effector memory cells reexpressing CD45RA (Temra) CD8+ T cells, memory CD4+ T cells, including Th1 and Th1*, Tregs, and memory B cells, whereas the counts and proportions of other leukocyte subsets were normal. Functional deficits of myeloid cells included the abolition of IL-12 and IL-23 production by conventional DC1s (cDC1s) and monocytes, but not cDC2s. c-Rel was also required for induction of CD86 expression on, and thus antigen-presenting cell function of, cDCs. Functional deficits of lymphoid cells included reduced IL-2 production by naive T cells, correlating with low proliferation and survival rates and poor production of Th1, Th2, and Th17 cytokines by memory CD4+ T cells. In naive CD4+ T cells, c-Rel is dispensable for early IL2 induction but contributes to later phases of IL2 expression. The patient's naive B cells displayed impaired MYC and BCL2L1 induction, compromising B cell survival and proliferation and preventing their differentiation into Ig-secreting plasmablasts. Inherited c-Rel deficiency disrupts the development and function of multiple myeloid and lymphoid cells, compromising innate and adaptive immunity to multiple infectious agents.
Subject(s)
Genes, rel , Primary Immunodeficiency Diseases/genetics , Primary Immunodeficiency Diseases/immunology , Proto-Oncogene Proteins c-rel/deficiency , Proto-Oncogene Proteins c-rel/genetics , Adaptive Immunity/genetics , Adaptive Immunity/immunology , Child , Consanguinity , Female , Hematopoietic Stem Cell Transplantation , Homozygote , Host Microbial Interactions/genetics , Host Microbial Interactions/immunology , Humans , Immunity, Innate/genetics , Immunity, Innate/immunology , Lymphocyte Activation , Lymphocytes/classification , Lymphocytes/immunology , Mutation , Myeloid Cells/immunology , Primary Immunodeficiency Diseases/therapy , Protein IsoformsABSTRACT
The pathophysiology of adverse events following programmed cell death protein 1 (PD-1) blockade, including tuberculosis (TB) and autoimmunity, remains poorly characterized. We studied a patient with inherited PD-1 deficiency and TB who died of pulmonary autoimmunity. The patient's leukocytes did not express PD-1 or respond to PD-1-mediated suppression. The patient's lymphocytes produced only small amounts of interferon (IFN)-γ upon mycobacterial stimuli, similarly to patients with inborn errors of IFN-γ production who are vulnerable to TB. This phenotype resulted from a combined depletion of Vδ2+ γδ T, mucosal-associated invariant T and CD56bright natural killer lymphocytes and dysfunction of other T lymphocyte subsets. Moreover, the patient displayed hepatosplenomegaly and an expansion of total, activated and RORγT+ CD4-CD8- double-negative αß T cells, similar to patients with STAT3 gain-of-function mutations who display lymphoproliferative autoimmunity. This phenotype resulted from excessive amounts of STAT3-activating cytokines interleukin (IL)-6 and IL-23 produced by activated T lymphocytes and monocytes, and the STAT3-dependent expression of RORγT by activated T lymphocytes. Our work highlights the indispensable role of human PD-1 in governing both antimycobacterial immunity and self-tolerance, while identifying potentially actionable molecular targets for the diagnostic and therapeutic management of TB and autoimmunity in patients on PD-1 blockade.
Subject(s)
Autoimmunity/genetics , Nuclear Receptor Subfamily 1, Group F, Member 3/genetics , Programmed Cell Death 1 Receptor/genetics , STAT3 Transcription Factor/genetics , Tuberculosis/immunology , Autoimmunity/immunology , CD4-Positive T-Lymphocytes/immunology , CD4-Positive T-Lymphocytes/pathology , CD56 Antigen/genetics , CD8-Positive T-Lymphocytes/immunology , CD8-Positive T-Lymphocytes/pathology , Child , Humans , Immune Checkpoint Inhibitors/adverse effects , Immune Checkpoint Inhibitors/therapeutic use , Immunotherapy/adverse effects , Interleukin-23/genetics , Interleukin-6/genetics , Intraepithelial Lymphocytes/immunology , Intraepithelial Lymphocytes/pathology , Male , Mycobacterium tuberculosis/pathogenicity , Neoplasms/complications , Neoplasms/drug therapy , Neoplasms/mortality , Programmed Cell Death 1 Receptor/deficiency , Tuberculosis/genetics , Tuberculosis/mortalityABSTRACT
A recent report found that rare predicted loss-of-function (pLOF) variants across 13 candidate genes in TLR3- and IRF7-dependent type I IFN pathways explain up to 3.5% of severe COVID-19 cases. We performed whole-exome or whole-genome sequencing of 1,864 COVID-19 cases (713 with severe and 1,151 with mild disease) and 15,033 ancestry-matched population controls across 4 independent COVID-19 biobanks. We tested whether rare pLOF variants in these 13 genes were associated with severe COVID-19. We identified only 1 rare pLOF mutation across these genes among 713 cases with severe COVID-19 and observed no enrichment of pLOFs in severe cases compared to population controls or mild COVID-19 cases. We found no evidence of association of rare LOF variants in the 13 candidate genes with severe COVID-19 outcomes.