Henrikki Almusa

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

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

BackgroundTechniques 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.ResultsWe 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.ConclusionsLibraries 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.

Thymidine kinase 2 mutations in autosomal recessive progressive external ophthalmoplegia with multiple mitochondrial DNA deletions

Autosomal-inherited progressive external ophthalmoplegia (PEO) is an adult-onset disease characterized by the accumulation of multiple mitochondrial DNA (mtDNA) deletions in post-mitotic tissues. Mutations in six different genes have been described to cause the autosomal dominant form of the disease, but only mutations in the DNA polymerase gamma gene are known to cause autosomal recessive PEO (arPEO), leaving the genetic background of arPEO mostly unknown. Here we used whole-exome sequencing and identified compound heterozygous mutations, leading to two amino acid alterations R225W and a novel T230A in thymidine kinase 2 (TK2) in arPEO patients. TK2 is an enzyme of the mitochondrial nucleotide salvage pathway and its loss-of-function mutations have previously been shown to underlie the early-infantile myopathic form of mtDNA depletion syndrome (MDS). Our TK2 activity measurements of patient fibroblasts and mutant recombinant proteins show that the combination of the identified arPEO variants, R225W and T230A, leads to a significant reduction in TK2 activity, consistent with the late-onset phenotype, whereas homozygosity for R225W, previously associated with MDS, leads to near-total loss of activity. Our finding identifies a new genetic cause of arPEO with multiple mtDNA deletions. Furthermore, MDS and multiple mtDNA deletion disorders are manifestations of the same pathogenic pathways affecting mtDNA replication and repair, indicating that MDS-associated genes should be studied when searching for genetic background of PEO disorders.

A Single-Nucleotide Substitution Mutator Phenotype Revealed by Exome Sequencing of Human Colon Adenomas

Oncogene-induced DNA replication stress is thought to drive genomic instability in cancer. In particular, replication stress can explain the high prevalence of focal genomic deletions mapping within very large genes in human tumors. However, the origin of single-nucleotide substitutions (SNS) in nonfamilial cancers is strongly debated. Some argue that cancers have a mutator phenotype, whereas others argue that the normal DNA replication error rates are sufficient to explain the number of observed SNSs. Here, we sequenced the exomes of 24, mostly precancerous, colon polyps. Analysis of the sequences revealed mutations in the APC, CTNNB1, and BRAF genes as the presumptive cancer-initiating events and many passenger SNSs. We used the number of SNSs in the various lesions to calculate mutation rates for normal colon and adenomas and found that colon adenomas exhibit a mutator phenotype. Interestingly, the SNSs in the adenomas mapped more often than expected within very large genes, where focal deletions in response to DNA replication stress also map. We propose that single-stranded DNA generated in response to oncogene-induced replication stress compromises the repair of deaminated cytosines and other damaged bases, leading to the observed SNS mutator phenotype.

The Finnish disease heritage database (FinDis) update-a database for the genes mutated in the Finnish disease heritage brought to the next-generation sequencing era.

The Finnish Disease Heritage Database (FinDis) (http://findis.org) was originally published in 2004 as a centralized information resource for rare monogenic diseases enriched in the Finnish population. The FinDis database originally contained 405 causative variants for 30 diseases. At the time, the FinDis database was a comprehensive collection of data, but since 1994, a large amount of new information has emerged, making the necessity to update the database evident. We collected information and updated the database to contain genes and causative variants for 35 diseases, including six more genes and more than 1,400 additional disease‐causing variants. Information for causative variants for each gene is collected under the LOVD 3.0 platform, enabling easy updating. The FinDis portal provides a centralized resource and user interface to link information on each disease and gene with variant data in the LOVD 3.0 platform. The software written to achieve this has been open‐sourced and made available on GitHub (http://github.com/findis‐db), allowing biomedical institutions in other countries to present their national data in a similar way, and to both contribute to, and benefit from, standardized variation data. The updated FinDis portal provides a unique resource to assist patient diagnosis, research, and the development of new cures.

Full-Genome Sequences of Influenza A(H1N1)pdm09 Viruses Isolated from Finnish Patients from 2009 to 2013

ABSTRACT Here we report full-length sequencing of the first large set of influenza A(H1N1)pdm09 virus genomes isolated in Finland between the years 2009 and 2013 and discuss the advantages and needs of influenza virus sequencing efforts.

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.

Monitoring therapy responses at the leukemic subclone level by ultra-deep amplicon resequencing in acute myeloid leukemia

In our individualized systems medicine program, personalized treatment options are identified and administered to chemorefractory acute myeloid leukemia (AML) patients based on exome sequencing and ex vivo drug sensitivity and resistance testing data. Here, we analyzed how clonal heterogeneity affects the responses of 13 AML patients to chemotherapy or targeted treatments using ultra-deep (average 68 000 × coverage) amplicon resequencing. Using amplicon resequencing, we identified 16 variants from 4 patients (frequency 0.54–2%) that were not detected previously by exome sequencing. A correlation-based method was developed to detect mutation-specific responses in serial samples across multiple time points. Significant subclone-specific responses were observed for both chemotherapy and targeted therapy. We detected subclonal responses in patients where clinical European LeukemiaNet (ELN) criteria showed no response. Subclonal responses also helped to identify putative mechanisms underlying drug sensitivities, such as sensitivity to azacitidine in DNMT3A mutated cell clones and resistance to cytarabine in a subclone with loss of NF1 gene. In summary, ultra-deep amplicon resequencing method enables sensitive quantification of subclonal variants and their responses to therapies. This approach provides new opportunities for designing combinatorial therapies blocking multiple subclones as well as for real-time assessment of such treatments.

Complete Genome Sequences of Influenza A/H1N1 Strains Isolated from Patients during the 2013-2014 Epidemic Season in Finland

ABSTRACT Here, we report 40 complete genome sequences of influenza A/H1N1 strains isolated from 33 nonhospitalized and 7 hospitalized patients during the 2013-2014 epidemic season in Finland. An analysis of the aligned sequences revealed no oseltamivir-resistant genotypes. As a whole, the recent viruses have drifted from the prototype A/California/7/2009 virus by ca. 1.3%.

Influenza pH1N1 Virus Accumulated H275Y Mutation in Neuraminidase during Propagation in MDCK Cells

ABSTRACT Here, we sequenced the genome of the influenza A/Finland/741 M/2014(H1N1) virus and found that the virus accumulated oseltamivir resistance mutation H275Y in its neuraminidase during propagation in cell culture. This indicates that propagation in cell culture modifies virus genomes. The instability of influenza genomes should be taken into consideration during drug-sensitivity studies.