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.

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.

Discovery of mitochondrial DNA variants associated with genome-wide blood cell gene expression: a population-based mtDNA sequencing study.

The effect of mitochondrial DNA (mtDNA) variation on peripheral blood transcriptomics in health and disease is not fully known. Sex-specific mitochondrially controlled gene expression patterns have been shown in Drosophila melanogaster but in humans, evidence is lacking. Functional variation in mtDNA may also have a role in the development of type 2 diabetes and its precursor state, i.e. prediabetes. We examined the associations between mitochondrial single-nucleotide polymorphisms (mtSNPs) and peripheral blood transcriptomics with a focus on sex- and prediabetes-specific effects. The genome-wide blood cell expression data of 19 637 probes, 199 deep-sequenced mtSNPs and nine haplogroups of 955 individuals from a population-based Young Finns Study cohort were used. Significant associations were identified with linear regression and analysis of covariance. The effects of sex and prediabetes on the associations between gene expression and mtSNPs were studied using random-effect meta-analysis. Our analysis identified 53 significant expression probe-mtSNP associations after Bonferroni correction, involving 7 genes and 31 mtSNPs. Eight probe-mtSNP signals remained independent after conditional analysis. In addition, five genes showed differential expression between haplogroups. The meta-analysis did not show any significant differences in linear model effect sizes between males and females but identified the association between the OASL gene and mtSNP C16294T to show prediabetes-specific effects. This study pinpoints new independent mtSNPs associated with peripheral blood transcriptomics and replicates six previously reported associations, providing further evidence of the mitochondrial genetic control of blood cell gene expression. In addition, we present evidence that prediabetes might lead to perturbations in mitochondrial control.

Somatic STAT3 mutations in large granular lymphocytic leukemia.

BACKGROUND: T-cell large granular lymphocytic leukemia is a rare lymphoproliferative disorder characterized by the expansion of clonal CD3+CD8+ cytotoxic T lymphocytes (CTLs) and often associated with autoimmune disorders and immune-mediated cytopenias. METHODS: We used next-generation exome sequencing to identify somatic mutations in CTLs from an index patient with large granular lymphocytic leukemia. Targeted resequencing was performed in a well-characterized cohort of 76 patients with this disorder, characterized by clonal T-cell-receptor rearrangements and increased numbers of large granular lymphocytes. RESULTS: Mutations in the signal transducer and activator of transcription 3 gene (STAT3) were found in 31 of 77 patients (40%) with large granular lymphocytic leukemia. Among these 31 patients, recurrent mutational hot spots included Y640F in 13 (17%), D661V in 7 (9%), D661Y in 7 (9%), and N647I in 3 (4%). All mutations were located in exon 21, encoding the Src homology 2 (SH2) domain, which mediates the dimerization and activation of STAT protein. The amino acid changes resulted in a more hydrophobic protein surface and were associated with phosphorylation of STAT3 and its localization in the nucleus. In vitro functional studies showed that the Y640F and D661V mutations increased the transcriptional activity of STAT3. In the affected patients, downstream target genes of the STAT3 pathway (IFNGR2, BCL2L1, and JAK2) were up-regulated. Patients with STAT3 mutations presented more often with neutropenia and rheumatoid arthritis than did patients without these mutations. CONCLUSIONS: The SH2 dimerization and activation domain of STAT3 is frequently mutated in patients with large granular lymphocytic leukemia; these findings suggest that aberrant STAT3 signaling underlies the pathogenesis of this disease. (Funded by the Academy of Finland and others.).

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.

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.].

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.

TaME-seq: An efficient sequencing approach for characterisation of HPV genomic variability and chromosomal integration.

HPV genomic variability and chromosomal integration are important in the HPV-induced carcinogenic process. To uncover these genomic events in an HPV infection, we have developed an innovative and cost-effective sequencing approach named TaME-seq (tagmentation-assisted multiplex PCR enrichment sequencing). TaME-seq combines tagmentation and multiplex PCR enrichment for simultaneous analysis of HPV variation and chromosomal integration, and it can also be adapted to other viruses. For method validation, cell lines (n = 4), plasmids (n = 3), and HPV16, 18, 31, 33 and 45 positive clinical samples (n = 21) were analysed. Our results showed deep HPV genome-wide sequencing coverage. Chromosomal integration breakpoints and large deletions were identified in HPV positive cell lines and in one clinical sample. HPV genomic variability was observed in all samples allowing identification of low frequency variants. In contrast to other approaches, TaME-seq proved to be highly efficient in HPV target enrichment, leading to reduced sequencing costs. Comprehensive studies on HPV intra-host variability generated during a persistent infection will improve our understanding of viral carcinogenesis. Efficient identification of both HPV variability and integration sites will be important for the study of HPV evolution and adaptability and may be an important tool for use in cervical cancer diagnostics.

Genome-wide DNA methylation in saliva and body size of adolescent girls.

AIM: We performed an epigenome-wide association study within the Finnish Health in Teens cohort to identify differential DNA methylation and its association with BMI in adolescents. MATERIALS & METHODS: Differential DNA methylation analyses of 3.1 million CpG sites were performed in saliva samples from 50 lean and 50 heavy adolescent girls by genome-wide targeted bisulfite-sequencing. RESULTS: We identified 100 CpG sites with p-values < 0.000524, seven regions by 'bumphunting' and five CpG islands that differed significantly between the two groups. The ten CpG sites and regions most strongly associated with BMI substantially overlapped with obesity- and insulin-related genes, including MC2R, IGFBPL1, IP6K1 and IGF2BP1. CONCLUSION: Our findings suggest an association between the saliva methylome and BMI in adolescence.

Enrichment of rare variants in population isolates: single AICDA mutation responsible for hyper-IgM syndrome type 2 in Finland.

Antibody class-switch recombination and somatic hypermutation critically depend on the function of activation-induced cytidine deaminase (AID). Rare variants in its gene AICDA have been reported to cause autosomal recessive AID deficiency (autosomal recessive hyper-IgM syndrome type 2 (HIGM2)). Exome sequencing of a multicase Finnish family with an HIGM2 phenotype identified a rare, homozygous, variant (c.416T>C, p.(Met139Thr)) in the AICDA gene, found to be significantly enriched in the Finnish population compared with other populations of European origin (38.56-fold, P<0.001). The population history of Finland, characterized by a restricted number of founders, isolation and several population bottlenecks, has caused enrichment of certain rare disease-causing variants and losses of others, as part of a phenomenon called the Finnish Disease Heritage. Accordingly, rare founder mutations cause the majority of observed Finnish cases in these mostly autosomal recessive disorders that consequently are more frequent in Finland than elsewhere. Screening of all currently known Finnish patients with an HIGM2 phenotype showed them to be homozygous for p.(Met139Thr). All the Finnish p.(Met139Thr) carriers with available data on their geographic descent originated from the eastern and northeastern parts of Finland. They were observed to share more of their genome identity by descent (IBD) than Finns in general (P<0.001), and they all carried a 207.5-kb ancestral haplotype containing the variant. In conclusion, the identified p.(Met139Thr) variant is significantly enriched in Finns and explains all thus far found AID deficiencies in Finland.

mtDNA Mutagenesis Disrupts Pluripotent Stem Cell Function by Altering Redox Signaling.

mtDNA mutagenesis in somatic stem cells leads to their dysfunction and to progeria in mouse. The mechanism was proposed to involve modification of reactive oxygen species (ROS)/redox signaling. We studied the effect of mtDNA mutagenesis on reprogramming and stemness of pluripotent stem cells (PSCs) and show that PSCs select against specific mtDNA mutations, mimicking germline and promoting mtDNA integrity despite their glycolytic metabolism. Furthermore, mtDNA mutagenesis is associated with an increase in mitochondrial H2O2, reduced PSC reprogramming efficiency, and self-renewal. Mitochondria-targeted ubiquinone, MitoQ, and N-acetyl-L-cysteine efficiently rescued these defects, indicating that both reprogramming efficiency and stemness are modified by mitochondrial ROS. The redox sensitivity, however, rendered PSCs and especially neural stem cells sensitive to MitoQ toxicity. Our results imply that stem cell compartment warrants special attention when the safety of new antioxidants is assessed and point to an essential role for mitochondrial redox signaling in maintaining normal stem cell function.

Individualized systems medicine strategy to tailor treatments for patients with chemorefractory acute myeloid leukemia.

UNLABELLED: We present an individualized systems medicine (ISM) approach to optimize cancer drug therapies one patient at a time. ISM is based on (i) molecular profiling and ex vivo drug sensitivity and resistance testing (DSRT) of patients' cancer cells to 187 oncology drugs, (ii) clinical implementation of therapies predicted to be effective, and (iii) studying consecutive samples from the treated patients to understand the basis of resistance. Here, application of ISM to 28 samples from patients with acute myeloid leukemia (AML) uncovered five major taxonomic drug-response subtypes based on DSRT profiles, some with distinct genomic features (e.g., MLL gene fusions in subgroup IV and FLT3-ITD mutations in subgroup V). Therapy based on DSRT resulted in several clinical responses. After progression under DSRT-guided therapies, AML cells displayed significant clonal evolution and novel genomic changes potentially explaining resistance, whereas ex vivo DSRT data showed resistance to the clinically applied drugs and new vulnerabilities to previously ineffective drugs. SIGNIFICANCE: Here, we demonstrate an ISM strategy to optimize safe and effective personalized cancer therapies for individual patients as well as to understand and predict disease evolution and the next line of therapy. This approach could facilitate systematic drug repositioning of approved targeted drugs as well as help to prioritize and de-risk emerging drugs for clinical testing.

Comparison of solution-based exome capture methods for next generation sequencing.

BACKGROUND: Techniques enabling targeted re-sequencing of the protein coding sequences of the human genome on next generation sequencing instruments are of great interest. We conducted a systematic comparison of the solution-based exome capture kits provided by Agilent and Roche NimbleGen. A control DNA sample was captured with all four capture methods and prepared for Illumina GAII sequencing. Sequence data from additional samples prepared with the same protocols were also used in the comparison. RESULTS: We developed a bioinformatics pipeline for quality control, short read alignment, variant identification and annotation of the sequence data. In our analysis, a larger percentage of the high quality reads from the NimbleGen captures than from the Agilent captures aligned to the capture target regions. High GC content of the target sequence was associated with poor capture success in all exome enrichment methods. Comparison of mean allele balances for heterozygous variants indicated a tendency to have more reference bases than variant bases in the heterozygous variant positions within the target regions in all methods. There was virtually no difference in the genotype concordance compared to genotypes derived from SNP arrays. A minimum of 11× coverage was required to make a heterozygote genotype call with 99% accuracy when compared to common SNPs on genome-wide association arrays. CONCLUSIONS: Libraries captured with NimbleGen kits aligned more accurately to the target regions. The updated NimbleGen kit most efficiently covered the exome with a minimum coverage of 20×, yet none of the kits captured all the Consensus Coding Sequence annotated exons.