Single-cell functional genomics reveals determinants of sensitivity and resistance to natural killer cells in blood cancers.

Cancer cells can evade natural killer (NK) cell activity, thereby limiting anti-tumor immunity. To reveal genetic determinants of susceptibility to NK cell activity, we examined interacting NK cells and blood cancer cells using single-cell and genome-scale functional genomics screens. Interaction of NK and cancer cells induced distinct activation and type I interferon (IFN) states in both cell types depending on the cancer cell lineage and molecular phenotype, ranging from more sensitive myeloid to less sensitive B-lymphoid cancers. CRISPR screens in cancer cells uncovered genes regulating sensitivity and resistance to NK cell-mediated killing, including adhesion-related glycoproteins, protein fucosylation genes, and transcriptional regulators, in addition to confirming the importance of antigen presentation and death receptor signaling pathways. CRISPR screens with a single-cell transcriptomic readout provided insight into underlying mechanisms, including regulation of IFN-γ signaling in cancer cells and NK cell activation states. Our findings highlight the diversity of mechanisms influencing NK cell susceptibility across different cancers and provide a resource for NK cell-based therapies.

The phylodynamics of SARS-CoV-2 during 2020 in Finland.

Background: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused millions of infections and fatalities globally since its emergence in late 2019. The virus was first detected in Finland in January 2020, after which it rapidly spread among the populace in spring. However, compared to other European nations, Finland has had a low incidence of SARS-CoV-2. To gain insight into the origins and turnover of SARS-CoV-2 lineages circulating in Finland in 2020, we investigated the phylogeographic and -dynamic history of the virus. Methods: The origins of SARS-CoV-2 introductions were inferred via Travel-aware Bayesian time-measured phylogeographic analyses. Sequences for the analyses included virus genomes belonging to the B.1 lineage and with the D614G mutation from countries of likely origin, which were determined utilizing Google mobility data. We collected all available sequences from spring and fall peaks to study lineage dynamics. Results: We observed rapid turnover among Finnish lineages during this period. Clade 20C became the most prevalent among sequenced cases and was replaced by other strains in fall 2020. Bayesian phylogeographic reconstructions suggested 42 independent introductions into Finland during spring 2020, mainly from Italy, Austria, and Spain. Conclusions: A single introduction from Spain might have seeded one-third of cases in Finland during spring in 2020. The investigations of the original introductions of SARS-CoV-2 to Finland during the early stages of the pandemic and of the subsequent lineage dynamics could be utilized to assess the role of transboundary movements and the effects of early intervention and public health measures.

NCOR2 is a novel candidate gene for migraine-epilepsy phenotype.

HYPOTHESIS: To identify genetic factors predisposing to migraine-epilepsy phenotype utilizing a multi-generational family with known linkage to chr12q24.2-q24.3. METHODS: We used single nucleotide polymorphism (SNP) genotyping and next-generation sequencing technologies to perform linkage, haplotype, and variant analyses in an extended Finnish migraine-epilepsy family (n = 120). In addition, we used a large genome-wide association study (GWAS) dataset of migraine and two biobank studies, UK Biobank and FinnGen, to test whether variants within the susceptibility region associate with migraine or epilepsy related phenotypes in a population setting. RESULTS: The family showed the highest evidence of linkage (LOD 3.42) between rs7966411 and epilepsy. The haplotype shared among 12 out of 13 epilepsy patients in the family covers almost the entire NCOR2 and co-localizes with one of the risk loci of the recent GWAS on migraine. The haplotype harbors nine low-frequency variants with potential regulatory functions. Three of them, in addition to two common variants, show nominal associations with neurological disorders in either UK Biobank or FinnGen. CONCLUSION: We provide several independent lines of evidence supporting association between migraine-epilepsy phenotype and NCOR2. Our study suggests that NCOR2 may have a role in both migraine and epilepsy and thus would provide evidence for shared pathophysiology underlying these two diseases.

Mitochondrial genome-wide analysis of nuclear DNA methylation quantitative trait loci.

Mitochondria have a complex communication network with the surrounding cell and can alter nuclear DNA methylation (DNAm). Variation in the mitochondrial DNA (mtDNA) has also been linked to differential DNAm. Genome-wide association studies have identified numerous DNAm quantitative trait loci, but these studies have not examined the mitochondrial genome. Herein, we quantified nuclear DNAm from blood and conducted a mitochondrial genome-wide association study of DNAm, with an additional emphasis on sex- and prediabetes-specific heterogeneity. We used the Young Finns Study (n = 926) with sequenced mtDNA genotypes as a discovery sample and sought replication in the Ludwigshafen Risk and Cardiovascular Health study (n = 2317). We identified numerous significant associations in the discovery phase (P < 10-9), but they were not replicated when accounting for multiple testing. In total, 27 associations were nominally replicated with a P < 0.05. The replication analysis presented no evidence of sex- or prediabetes-specific heterogeneity. The 27 associations were included in a joint meta-analysis of the two cohorts, and 19 DNAm sites associated with mtDNA variants, while four other sites showed haplogroup associations. An expression quantitative trait methylation analysis was performed for the identified DNAm sites, pinpointing two statistically significant associations. This study provides evidence of a mitochondrial genetic control of nuclear DNAm with little evidence found for sex- and prediabetes-specific effects. The lack of a comparable mtDNA data set for replication is a limitation in our study and further studies are needed to validate our results.

Molecular features encoded in the ctDNA reveal heterogeneity and predict outcome in high-risk aggressive B-cell lymphoma.

Inadequate molecular and clinical stratification of the patients with high-risk diffuse large B-cell lymphoma (DLBCL) is a clinical challenge hampering the establishment of personalized therapeutic options. We studied the translational significance of liquid biopsy in a uniformly treated trial cohort. Pretreatment circulating tumor DNA (ctDNA) revealed hidden clinical and biological heterogeneity, and high ctDNA burden determined increased risk of relapse and death independently of conventional risk factors. Genomic dissection of pretreatment ctDNA revealed translationally relevant phenotypic, molecular, and prognostic information that extended beyond diagnostic tissue biopsies. During therapy, chemorefractory lymphomas exhibited diverging ctDNA kinetics, whereas end-of-therapy negativity for minimal residual disease (MRD) characterized cured patients and resolved clinical enigmas, including false residual PET positivity. Furthermore, we discovered fragmentation disparities in the cell-free DNA that characterize lymphoma-derived ctDNA and, as a proof-of-concept for their clinical application, used machine learning to show that end-of-therapy fragmentation patterns predict outcome. Altogether, we have discovered novel molecular determinants in the liquid biopsy that can noninvasively guide treatment decisions.

Incidence Trends for SARS-CoV-2 Alpha and Beta Variants, Finland, Spring 2021.

Severe acute respiratory syndrome coronavirus 2 Alpha and Beta variants became dominant in Finland in spring 2021 but had diminished by summer. We used phylogenetic clustering to identify sources of spreading. We found that outbreaks were mostly seeded by a few introductions, highlighting the importance of surveillance and prevention policies.

The CRCbiome study: a large prospective cohort study examining the role of lifestyle and the gut microbiome in colorectal cancer screening participants.

BACKGROUND: Colorectal cancer (CRC) screening reduces CRC incidence and mortality. However, current screening methods are either hampered by invasiveness or suboptimal performance, limiting their effectiveness as primary screening methods. To aid in the development of a non-invasive screening test with improved sensitivity and specificity, we have initiated a prospective biomarker study (CRCbiome), nested within a large randomized CRC screening trial in Norway. We aim to develop a microbiome-based classification algorithm to identify advanced colorectal lesions in screening participants testing positive for an immunochemical fecal occult blood test (FIT). We will also examine interactions with host factors, diet, lifestyle and prescription drugs. The prospective nature of the study also enables the analysis of changes in the gut microbiome following the removal of precancerous lesions. METHODS: The CRCbiome study recruits participants enrolled in the Bowel Cancer Screening in Norway (BCSN) study, a randomized trial initiated in 2012 comparing once-only sigmoidoscopy to repeated biennial FIT, where women and men aged 50-74 years at study entry are invited to participate. Since 2017, participants randomized to FIT screening with a positive test result have been invited to join the CRCbiome study. Self-reported diet, lifestyle and demographic data are collected prior to colonoscopy after the positive FIT-test (baseline). Screening data, including colonoscopy findings are obtained from the BCSN database. Fecal samples for gut microbiome analyses are collected both before and 2 and 12 months after colonoscopy. Samples are analyzed using metagenome sequencing, with taxonomy profiles, and gene and pathway content as primary measures. CRCbiome data will also be linked to national registries to obtain information on prescription histories and cancer relevant outcomes occurring during the 10 year follow-up period. DISCUSSION: The CRCbiome study will increase our understanding of how the gut microbiome, in combination with lifestyle and environmental factors, influences the early stages of colorectal carcinogenesis. This knowledge will be crucial to develop microbiome-based screening tools for CRC. By evaluating biomarker performance in a screening setting, using samples from the target population, the generalizability of the findings to future screening cohorts is likely to be high. TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT01538550 .

Corrigendum: Adult-Onset Anti-Citrullinated Peptide Antibody-Negative Destructive Rheumatoid Arthritis Is Characterized by a Disease-Specific CD8+ T Lymphocyte Signature.

[This corrects the article DOI: 10.3389/fimmu.2020.578848.].

Somatic mutations in lymphocytes in patients with immune-mediated aplastic anemia.

The prevalence and functional impact of somatic mutations in nonleukemic T cells is not well characterized, although clonal T-cell expansions are common. In immune-mediated aplastic anemia (AA), cytotoxic T-cell expansions are shown to participate in disease pathogenesis. We investigated the mutation profiles of T cells in AA by a custom panel of 2533 genes. We sequenced CD4+ and CD8+ T cells of 24 AA patients and compared the results to 20 healthy controls and whole-exome sequencing of 37 patients with AA. Somatic variants were common both in patients and healthy controls but enriched to AA patients' CD8+ T cells, which accumulated most mutations on JAK-STAT and MAPK pathways. Mutation burden was associated with CD8+ T-cell clonality, assessed by T-cell receptor beta sequencing. To understand the effect of mutations, we performed single-cell sequencing of AA patients carrying STAT3 or other mutations in CD8+ T cells. STAT3 mutated clone was cytotoxic, clearly distinguishable from other CD8+ T cells, and attenuated by successful immunosuppressive treatment. Our results suggest that somatic mutations in T cells are common, associate with clonality, and can alter T-cell phenotype, warranting further investigation of their role in the pathogenesis of AA.

Examining the effect of mitochondrial DNA variants on blood pressure in two Finnish cohorts.

High blood pressure (BP) is a major risk factor for many noncommunicable diseases. The effect of mitochondrial DNA single-nucleotide polymorphisms (mtSNPs) on BP is less known than that of nuclear SNPs. We investigated the mitochondrial genetic determinants of systolic, diastolic, and mean arterial BP. MtSNPs were determined from peripheral blood by sequencing or with genome-wide association study SNP arrays in two independent Finnish cohorts, the Young Finns Study and the Finnish Cardiovascular Study, respectively. In total, over 4200 individuals were included. The effects of individual common mtSNPs, with an additional focus on sex-specificity, and aggregates of rare mtSNPs grouped by mitochondrial genes were evaluated by meta-analysis of linear regression and a sequence kernel association test, respectively. We accounted for the predicted pathogenicity of the rare variants within protein-encoding and the tRNA regions. In the meta-analysis of 87 common mtSNPs, we did not observe significant associations with any of the BP traits. Sex-specific and rare-variant analyses did not pinpoint any significant associations either. Our results are in agreement with several previous studies suggesting that mtDNA variation does not have a significant role in the regulation of BP. Future studies might need to reconsider the mechanisms thought to link mtDNA with hypertension.

Identifying haplotypes in recessive inherited retinal dystrophies using whole-genome linked-read sequencing.

Conventional next-generation sequencing methods, used in most gene panels, cannot separate maternally and paternally derived sequence information of distant variants. In recessive diseases, two or more equally plausible causative variants with unsolved phase information prevent accurate molecular diagnosis. In reality, close relatives might be unavailable for segregation analysis. Here, we utilized whole genome linked-read sequencing to assign variants to haplotypes in two patients with inherited retinal dystrophies. Patient 1 with macular dystrophy had variants c.3442T>C, p.(Cys1148Arg), c.4209G>T, p.(Glu1403Asp), and c.1182C>T, p.(Cys394=) in CRB1, and Patient 2 with nonsyndromic retinitis pigmentosa had c.1328T>A, p.(Val443Asp) and c.3032C>G, p.(Ser1011*) in AHI1. The relatives were not available for genotyping. Using whole genome linked-read sequencing we phased the variants to haplotypes providing genetic background for the retinal dystrophies. In future, when the price of sequencing methods that provides long-read data decreases and their read-depth and accuracy increases, they are probably considered the primary or adjunctive sequencing methods in genetic testing, allowing the immediate collection of phase information and thus obviating the need for the carrier testing and segregation analysis.

Adult-Onset Anti-Citrullinated Peptide Antibody-Negative Destructive Rheumatoid Arthritis Is Characterized by a Disease-Specific CD8+ T Lymphocyte Signature.

Rheumatoid arthritis (RA) is a complex autoimmune disease targeting synovial joints. Traditionally, RA is divided into seropositive (SP) and seronegative (SN) disease forms, the latter consisting of an array of unrelated diseases with joint involvement. Recently, we described a severe form of SN-RA that associates with characteristic joint destruction. Here, we sought biological characteristics to differentiate this rare but aggressive anti-citrullinated peptide antibody-negative destructive RA (CND-RA) from early seropositive (SP-RA) and seronegative rheumatoid arthritis (SN-RA). We also aimed to study cytotoxic CD8+ lymphocytes in autoimmune arthritis. CND-RA, SP-RA and SN-RA were compared to healthy controls to reveal differences in T-cell receptor beta (TCRß) repertoire, cytokine levels and autoantibody repertoires. Whole-exome sequencing (WES) followed by single-cell RNA-sequencing (sc-RNA-seq) was performed to study somatic mutations in a clonally expanded CD8+ lymphocyte population in an index patient. A unique TCRß signature was detected in CND-RA patients. In addition, CND-RA patients expressed higher levels of the bone destruction-associated TNFSF14 cytokine. Blood IgG repertoire from CND-RA patients recognized fewer endogenous proteins than SP-RA patients' repertoires. Using WES, we detected a stable mutation profile in the clonally expanded CD8+ T-cell population characterized by cytotoxic gene expression signature discovered by sc-RNA-sequencing. Our results identify CND-RA as an independent RA subset and reveal a CND-RA specific TCR signature in the CD8+ lymphocytes. Improved classification of seronegative RA patients underlines the heterogeneity of RA and also, facilitates development of improved therapeutic options for the treatment resistant patients.

Somatic mutations and T-cell clonality in patients with immunodeficiency.

Common variable immunodeficiency and other late-onset immunodeficiencies often co-manifest with autoimmunity and lymphoproliferation. The pathogenesis of most cases is elusive, as only a minor subset harbors known monogenic germline causes. The involvement of both B and T cells is however implicated. To study whether somatic mutations in CD4+ and CD8+ T cells associate with immunodeficiency, we recruited 17 patients and 21 healthy controls. Eight patients had late-onset common variable immunodeficiency and nine patients other immunodeficiency and/or severe autoimmunity. In total, autoimmunity occurred in 94% and lymphoproliferation in 65%. We performed deep sequencing of 2533 immune-associated genes from CD4+ and CD8+ cells. Deep T-cell receptor beta sequencing was used to characterize CD4+ and CD8+ T-cell receptor repertoires. The prevalence of somatic mutations was 65% in all immunodeficiency patients, 75% in common variable immunodeficiency and 48% in controls. Clonal hematopoiesis-associated variants in both CD4+ and CD8+ cells occurred in 24% of immunodeficiency patients. Results demonstrated mutations in known tumor suppressors, oncogenes, and genes that are critical for immune- and proliferative functions, such as STAT5B (two patients), C5AR1 (two patients), KRAS (one patient), and NOD2 (one patient). Additionally, as a marker of T-cell receptor repertoire perturbation, common variable immunodeficiency patients harbored increased frequencies of clones with identical complementarity determining region 3 sequences despite unique nucleotide sequences when compared to controls. In conclusion, somatic mutations in genes implicated for autoimmunity and lymphoproliferation are common in CD4+ and CD8+ cells of patients with immunodeficiency. They may contribute to immune dysregulation in a subset of immunodeficiency patients.

Somatic mTOR mutation in clonally expanded T lymphocytes associated with chronic graft versus host disease.

Graft versus host disease (GvHD) is the main complication of allogeneic hematopoietic stem cell transplantation (HSCT). Here we report studies of a patient with chronic GvHD (cGvHD) carrying persistent CD4+ T cell clonal expansion harboring somatic mTOR, NFKB2, and TLR2 mutations. In the screening cohort (n = 134), we detect the mTOR P2229R kinase domain mutation in two additional cGvHD patients, but not in healthy or HSCT patients without cGvHD. Functional analyses of the mTOR mutation indicate a gain-of-function alteration and activation of both mTORC1 and mTORC2 signaling pathways, leading to increased cell proliferation and decreased apoptosis. Single-cell RNA sequencing and real-time impedance measurements support increased cytotoxicity of mutated CD4+ T cells. High throughput drug-sensitivity testing suggests that mutations induce resistance to mTOR inhibitors, but increase sensitivity for HSP90 inhibitors. Our findings imply that somatic mutations may contribute to aberrant T cell proliferations and persistent immune activation in cGvHD, thereby paving the way for targeted therapies.

IDH1 Expression via the R132H Mutation-Specific Antibody in Adrenocortical Neoplasias-Prognostic Impact in Carcinomas.

CONTEXT: Mutations to isocitrate dehydrogenase (IDH) appear to play a prognostic or predictive role in several neoplasias. Immunohistochemical staining designed to detect a specific R132H mutation to IDH1 showed expression in the normal adrenal cortex, raising interest to study the potential role of IDH1 in the pathogenesis of adrenocortical tumors. OBJECTIVE: The objective of this work is to study the role of IDH1 and its mutations in adrenocortical tumors. DESIGN AND PATIENTS: IDH1 R132H immunohistological staining was performed on a cohort of 197 adrenocortical tumors. The exon of the IDH1 gene was sequenced in 16 tumors. RESULTS: Positive IDH1 R132H immunohistochemical staining correlated with a better prognosis among patients with a malignant adrenocortical tumor. However, IDH1 R132H immunohistochemistry did not distinguish between local and metastasized tumors. We were unable to identify IDH1 mutations among our adrenocortical tumors using a targeted next-generation sequencing panel or via exon sequencing. CONCLUSIONS: Among adrenocortical carcinomas, IDH1 R132H immunopositivity correlated with a better prognosis. Thus, IDH1 R132H immunohistochemical staining could serve as a prognostic or as a potential predictive marker in adrenocortical carcinomas. Further research is needed to identify the possible alterations in IDH1 that could explain our findings, because we identified no known mutations to the IDH1 gene.

STAT3 Mutation Is Associated with STAT3 Activation in CD30+ ALK- ALCL.

Peripheral T-cell lymphomas (PTCL) are a heterogeneous, and often aggressive group of non-Hodgkin lymphomas. Recent advances in the molecular and genetic characterization of PTCLs have helped to delineate differences and similarities between the various subtypes, and the JAK/STAT pathway has been found to play an important oncogenic role. Here, we aimed to characterize the JAK/STAT pathway in PTCL subtypes and investigate whether the activation of the pathway correlates with the frequency of STAT gene mutations. Patient samples from AITL (n = 30), ALCL (n = 21) and PTCL-NOS (n = 12) cases were sequenced for STAT3, STAT5B, JAK1, JAK3, and RHOA mutations using amplicon sequencing and stained immunohistochemically for pSTAT3, pMAPK, and pAKT. We discovered STAT3 mutations in 13% of AITL, 13% of ALK+ ALCL, 38% of ALK- ALCL and 17% of PTCL-NOS cases. However, no STAT5B mutations were found and JAK mutations were only present in ALK- ALCL (15%). Concurrent mutations were found in all subgroups except ALK+ ALCL where STAT3 mutations were always seen alone. High pY-STAT3 expression was observed especially in AITL and ALCL samples. When studying JAK-STAT pathway mutations, pY-STAT3 expression was highest in PTCLs harboring either JAK1 or STAT3 mutations and CD30+ phenotype representing primarily ALK- ALCLs. Further investigation is needed to elucidate the molecular mechanisms of JAK-STAT pathway activation in PTCL.

Integrated drug profiling and CRISPR screening identify essential pathways for CAR T-cell cytotoxicity.

Chimeric antigen receptor (CAR) T-cell therapy has proven effective in relapsed and refractory B-cell malignancies, but resistance and relapses still occur. Better understanding of mechanisms influencing CAR T-cell cytotoxicity and the potential for modulation using small-molecule drugs could improve current immunotherapies. Here, we systematically investigated druggable mechanisms of CAR T-cell cytotoxicity using >500 small-molecule drugs and genome-scale CRISPR-Cas9 loss-of-function screens. We identified several tyrosine kinase inhibitors that inhibit CAR T-cell cytotoxicity by impairing T-cell signaling transcriptional activity. In contrast, the apoptotic modulator drugs SMAC mimetics sensitized B-cell acute lymphoblastic leukemia and diffuse large B-cell lymphoma cells to anti-CD19 CAR T cells. CRISPR screens identified death receptor signaling through FADD and TNFRSF10B (TRAIL-R2) as a key mediator of CAR T-cell cytotoxicity and elucidated the RIPK1-dependent mechanism of sensitization by SMAC mimetics. Death receptor expression varied across genetic subtypes of B-cell malignancies, suggesting a link between mechanisms of CAR T-cell cytotoxicity and cancer genetics. These results implicate death receptor signaling as an important mediator of cancer cell sensitivity to CAR T-cell cytotoxicity, with potential for pharmacological targeting to enhance cancer immunotherapy. The screening data provide a resource of immunomodulatory properties of cancer drugs and genetic mechanisms influencing CAR T-cell cytotoxicity.

SLC18A3 variants lead to fetal akinesia deformation sequence early in pregnancy.

Fetal akinesia deformation sequence (FADS) and lethal multiple pterygium syndrome (LMPS) are clinically overlapping syndromes manifesting with reduced or absent fetal movement, arthrogryposis, and several anomalies during fetal life. The etiology of these syndromes is heterogeneous, and in many cases it remains unknown. In order to determine the genetic etiology of FADS in two fetuses with fetal akinesia, arthrogryposis, edema, and partial cleft palate, we utilized exome sequencing. Our investigations revealed a homozygous nonsense variant [c.1116C>A, p.(Cys372Ter)] in the SLC18A3 gene, which encodes for the vesicular acetylcholine transporter (VAChT) responsible for active transport of acetylcholine in the neuromuscular junction. This is the first description of a nonsense variant in the SLC18A3 gene, as only missense variants and whole gene deletions have been previously identified in patients. The previously detected SLC18A3 defects have been associated with congenital myasthenic syndromes, and therefore our findings extend the clinical spectrum of SLC18A3 defects to severe prenatal phenotypes. Our findings suggest that nonsense variants in SLC18A3 cause a more severe phenotype than missense variants and are in line with previous studies showing a lethal phenotype in VAChT knockout mice. Our results underline the importance of including SLC18A3 sequencing in the differential diagnostics of fetuses with arthrogryposis, FADS, or LMPS of unknown etiology.