Haplotype information of large neuromuscular disease genes provided by linked-read sequencing has a potential to increase diagnostic yield.

Rare or novel missense variants in large genes such as TTN and NEB are frequent in the general population, which hampers the interpretation of putative disease-causing biallelic variants in patients with sporadic neuromuscular disorders. Often, when the first initial genetic analysis is performed, the reconstructed haplotype, i.e. phasing information of the variants is missing. Segregation analysis increases the diagnostic turnaround time and is not always possible if samples from family members are lacking. To overcome this difficulty, we investigated how well the linked-read technology succeeded to phase variants in these large genes, and whether it improved the identification of structural variants. Linked-read sequencing data of nemaline myopathy, distal myopathy, and proximal myopathy patients were analyzed for phasing, single nucleotide variants, and structural variants. Variant phasing was successful in the large muscle genes studied. The longest continuous phase blocks were gained using high-quality DNA samples with long DNA fragments. Homozygosity increased the number of phase blocks, especially in exome sequencing samples lacking intronic variation. In our cohort, linked-read sequencing added more information about the structural variation but did not lead to a molecular genetic diagnosis. The linked-read technology can support the clinical diagnosis of neuromuscular and other genetic disorders.

Congenital hypogonadotropic hypogonadism in a patient with a de novo POGZ mutation.

OBJECTIVE: Congenital hypogonadotropic hypogonadism (CHH) is a rare, genetically heterogeneous reproductive disorder caused by gonadotropin-releasing hormone (GnRH) deficiency. Approximately half of CHH patients also have decreased or absent sense of smell, that is, Kallmann syndrome (KS). We describe a patient with White-Sutton syndrome (developmental delay and autism spectrum disorder) and KS due to a heterozygous de novo mutation in POGZ (c.2857C>T, p.(Gln953*)), a gene encoding pogo transposable element derived with zinc finger domain, which acts as a transcriptomic regulator of neuronal networks. DESIGN AND METHODS: We modeled the role of POGZ in CHH by generating 2 clonal human pluripotent stem cell lines with CRISPR/Cas9, carrying either the heterozygous patient mutation (H11 line) or a homozygous mutation (c.2803-2906del; p.E935Kfs*7 encoding a truncated POGZ protein; F6del line). RESULTS: During the differentiation to GnRH neurons, neural progenitors derived from F6del line displayed severe proliferation defect, delayed wound-healing capacity, downregulation of intermediate progenitor neuron genes TBR1 and TBR2, and immature neuron markers PAX6 and TUBB3 and gave rise to fewer neurons with shorter neurites and less neurite branch points compared to the WT and H11 lines (P < .005). Both lines, however, could be successfully differentiated to GnRH neurons. CONCLUSIONS: In conclusion, this is the first report on the overlap between White-Sutton syndrome and CHH. POGZ mutations do not hinder GnRH neuron formation but may cause CHH/KS by affecting the size and motility of the anterior neural progenitor pool and neurite outgrowth.

Presentation and diagnosis of childhood-onset combined pituitary hormone deficiency: A single center experience from over 30 years.

Background: Childhood-onset combined pituitary hormone deficiency (CPHD) has a wide spectrum of etiologies and genetic causes for congenital disease. We aimed to describe the clinical spectrum and genetic etiologies of CPHD in a single tertiary center and estimate the population-level incidence of congenital CPHD. Methods: The retrospective clinical cohort comprised 124 CPHD patients (48 with congenital CPHD) treated at the Helsinki University Hospital (HUH) Children's Hospital between 1985 and 2018. Clinical data were collected from the patient charts. Whole exome sequencing was performed in 21 patients with congenital CPHD of unknown etiology. Findings: The majority (61%;76/124) of the patients had acquired CPHD, most frequently due to craniopharyngiomas and gliomas. The estimated incidence of congenital CPHD was 1/16 000 (95%CI, 1/11 000-1/24 000). The clinical presentation of congenital CPHD in infancy included prolonged/severe neonatal hypoglycaemia, prolonged jaundice, and/or micropenis/bilateral cryptorchidism in 23 (66%) patients; despite these clinical cues, only 76% of them were referred to endocrine investigations during the first year of life. The median delay between the first violation of the growth screening rules and the initiation of GH Rx treatment among all congenital CPHD patients was 2·2 years, interquartile range 1·2-3·7 years. Seven patients harbored pathogenic variants in PROP1, SOX3, TBC1D32, OTX2, and SOX2, and one patient carried a likely pathogenic variant in SHH (c.676G>A, p.(Ala226Thr)). Interpretation: Our study suggests that congenital CPHD can occur in 1/16 000 children, and that patients frequently exhibit neonatal cues of hypopituitarism and early height growth deflection. These results need to be corroborated in future studies and might inform clinical practice. Funding: Päivikki and Sakari Sohlberg Foundation, Biomedicum Helsinki Foundation, and Emil Aaltonen Foundation research grants.

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.

Kallmann syndrome in a patient with Weiss-Kruszka syndrome and a de novo deletion in 9q31.2.

Patients with deletions on chromosome 9q31.2 may exhibit delayed puberty, craniofacial phenotype including cleft lip/palate, and olfactory bulb hypoplasia. We report a patient with congenital HH with anosmia (Kallmann syndrome, KS) and a de novo 2.38 Mb heterozygous deletion in 9q31.2. The deletion breakpoints (determined with whole-genome linked-read sequencing) were in the FKTN gene (9:108,331,353) and in a non-coding area (9:110,707,332) (hg19). The deletion encompassed six protein-coding genes (FKTN, ZNF462, TAL2, TMEM38B, RAD23B, and KLF4). ZNF462 haploinsufficiency was consistent with the patient's Weiss-Kruszka syndrome (craniofacial phenotype, developmental delay, and sensorineural hearing loss), but did not explain his KS. In further analyses, he did not carry rare sequence variants in 32 known KS genes in whole-exome sequencing and displayed no aberrant splicing of 15 KS genes that were expressed in peripheral blood leukocyte transcriptome. The deletion was 1.8 Mb upstream of a KS candidate gene locus (PALM2AKAP2) but did not suppress its expression. In conclusion, this is the first report of a patient with Weiss-Kruszka syndrome and KS. We suggest that patients carrying a microdeletion in 9q31.2 should be evaluated for the presence of KS and KS-related features.

Loss of DIAPH1 causes SCBMS, combined immunodeficiency, and mitochondrial dysfunction.

BACKGROUND: Homozygous loss of DIAPH1 results in seizures, cortical blindness, and microcephaly syndrome (SCBMS). We studied 5 Finnish and 2 Omani patients with loss of DIAPH1 presenting with SCBMS, mitochondrial dysfunction, and immunodeficiency. OBJECTIVE: We sought to further characterize phenotypes and disease mechanisms associated with loss of DIAPH1. METHODS: Exome sequencing, genotyping and haplotype analysis, B- and T-cell phenotyping, in vitro lymphocyte stimulation assays, analyses of mitochondrial function, immunofluorescence staining for cytoskeletal proteins and mitochondria, and CRISPR-Cas9 DIAPH1 knockout in heathy donor PBMCs were used. RESULTS: Genetic analyses found all Finnish patients homozygous for a rare DIAPH1 splice-variant (NM_005219:c.684+1G>A) enriched in the Finnish population, and Omani patients homozygous for a previously described pathogenic DIAPH1 frameshift-variant (NM_005219:c.2769delT;p.F923fs). In addition to microcephaly, epilepsy, and cortical blindness characteristic to SCBMS, the patients presented with infection susceptibility due to defective lymphocyte maturation and 3 patients developed B-cell lymphoma. Patients' immunophenotype was characterized by poor lymphocyte activation and proliferation, defective B-cell maturation, and lack of naive T cells. CRISPR-Cas9 knockout of DIAPH1 in PBMCs from healthy donors replicated the T-cell activation defect. Patient-derived peripheral blood T cells exhibited impaired adhesion and inefficient microtubule-organizing center repositioning to the immunologic synapse. The clinical symptoms and laboratory tests also suggested mitochondrial dysfunction. Experiments with immortalized, patient-derived fibroblasts indicated that DIAPH1 affects the amount of complex IV of the mitochondrial respiratory chain. CONCLUSIONS: Our data demonstrate that individuals with SCBMS can have combined immune deficiency and implicate defective cytoskeletal organization and mitochondrial dysfunction in SCBMS pathogenesis.

Somatic mutation profiles as molecular classifiers of ulcerative colitis-associated colorectal cancer.

Ulcerative colitis increases colorectal cancer risk by mechanisms that remain incompletely understood. We approached this question by determining the genetic and epigenetic profiles of colitis-associated colorectal carcinomas (CA-CRC). The findings were compared to Lynch syndrome (LS), a different form of cancer predisposition that shares the importance of immunological factors in tumorigenesis. CA-CRCs (n = 27) were investigated for microsatellite instability, CpG island methylator phenotype and somatic mutations of 999 cancer-relevant genes ("Pan-cancer" panel). A subpanel of "Pan-cancer" design (578 genes) was used for LS colorectal tumors (n = 28). Mutational loads and signatures stratified CA-CRCs into three subgroups: hypermutated microsatellite-unstable (Group 1, n = 1), hypermutated microsatellite-stable (Group 2, n = 9) and nonhypermutated microsatellite-stable (Group 3, n = 17). The Group 1 tumor was the only one with MLH1 promoter hypermethylation and exhibited the mismatch repair deficiency-associated Signatures 21 and 15. Signatures 30 and 32 characterized Group 2, whereas no prominent single signature existed in Group 3. TP53, the most common mutational target in CA-CRC (16/27, 59%), was similarly affected in Groups 2 and 3, but DNA repair genes and Wnt signaling genes were mutated significantly more often in Group 2. In LS tumors, the degree of hypermutability exceeded that of the hypermutated CA-CRC Groups 1 and 2, and somatic mutational profiles and signatures were different. In conclusion, Groups 1 (4%) and 3 (63%) comply with published studies, whereas Group 2 (33%) is novel. The existence of molecularly distinct subgroups within CA-CRC may guide clinical management, such as therapy options.

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.

Recessive MYH3 variants cause "Contractures, pterygia, and variable skeletal fusions syndrome 1B" mimicking Escobar variant multiple pterygium syndrome.

The multiple pterygium syndromes (MPS) are rare disorders with disease severity ranging from lethal to milder forms. The nonlethal Escobar variant MPS (EVMPS) is characterized by multiple pterygia and arthrogryposis, as well as various additional features including congenital anomalies. The genetic etiology of EVMPS is heterogeneous and the diagnosis has been based either on the detection of pathogenic CHRNG variants (~23% of patients), or suggestive clinical features. We describe four patients with a clinical suspicion of EVMPS who manifested with multiple pterygia, mild flexion contractures of several joints, and vertebral anomalies. We revealed recessively inherited MYH3 variants as the underlying cause in all patients: two novel variants, c.1053C>G, p.(Tyr351Ter) and c.3102+5G>C, as compound heterozygous with the hypomorphic MYH3 variant c.-9+1G>A. Recessive MYH3 variants have been previously associated with spondylocarpotarsal synostosis syndrome. Our findings now highlight multiple pterygia as an important feature in patients with recessive MYH3 variants. Based on all patients with recessive MYH3 variants reported up to date, we consider that this disease entity should be designated as "Contractures, pterygia, and variable skeletal fusions syndrome 1B," as recently suggested by OMIM. Our findings underline the importance of analyzing MYH3 in the differential diagnosis of EVMPS, particularly as the hypomorphic MYH3 variant might remain undetected by routine exome sequencing.

Molecular Basis of Mismatch Repair Protein Deficiency in Tumors from Lynch Suspected Cases with Negative Germline Test Results.

Some 10-50% of Lynch-suspected cases with abnormal immunohistochemical (IHC) staining remain without any identifiable germline mutation of DNA mismatch repair (MMR) genes. MMR proteins form heterodimeric complexes, giving rise to distinct IHC patterns when mutant. Potential reasons for not finding a germline mutation include involvement of an MMR gene not predicted by the IHC pattern, epigenetic mechanism of predisposition, primary mutation in another DNA repair or replication-associated gene, and double somatic MMR gene mutations. We addressed these possibilities by germline and tumor studies in 60 Lynch-suspected cases ascertained through diagnostics (n = 55) or research (n = 5). All cases had abnormal MMR protein staining in tumors but no point mutation or large rearrangement of the suspected MMR genes in the germline. In diagnostic practice, MSH2/MSH6 (MutS Homolog 2/MutS Homolog 6) deficiency prompts MSH2 mutation screening; in our study, 3/11 index individuals (27%) with this IHC pattern revealed pathogenic germline mutations in MSH6. Individuals with isolated absence of MSH6 are routinely screened for MSH6 mutations alone; we found a predisposing mutation in MSH2 in 1/7 such cases (14%). Somatic deletion of the MSH2-MSH6 region, joint loss of MSH6 and MSH3 (MutS Homolog 3) proteins, and hindered MSH2/MSH6 dimerization offered explanations to misleading IHC patterns. Constitutional epimutation hypothesis was pursued in the MSH2 and/or MSH6-deficient cases plus 38 cases with MLH1 (MutL Homolog 1)-deficient tumors; a primary MLH1 epimutation was identified in one case with an MLH1-deficient tumor. We conclude that both MSH2 and MSH6 should be screened in MSH2/6- and MSH6-deficient cases. In MLH1-deficient cases, constitutional epimutations of MLH1 warrant consideration.

Does breast carcinoma belong to the Lynch syndrome tumor spectrum? - Somatic mutational profiles vs. ovarian and colorectal carcinomas.

Inherited DNA mismatch repair (MMR) defects cause predisposition to colorectal, endometrial, ovarian, and other cancers occurring in Lynch syndrome (LS). It is unsettled whether breast carcinoma belongs to the LS tumor spectrum. We approached this question through somatic mutational analysis of breast carcinomas from LS families, using established LS-spectrum tumors for comparison. Somatic mutational profiles of 578 cancer-relevant genes were determined for LS-breast cancer (LS-BC, n = 20), non-carrier breast cancer (NC-BC, n = 10), LS-ovarian cancer (LS-OC, n = 16), and LS-colorectal cancer (LS-CRC, n = 18) from the National LS Registry of Finland. Microsatellite and MMR protein analysis stratified LS-BCs into MMR-deficient (dMMR, n = 11) and MMR-proficient (pMMR, n = 9) subgroups. All NC-BCs were pMMR and all LS-OCs and LS-CRCs dMMR. All but one dMMR LS-BCs were hypermutated (> 10 non-synonymous mutations/Mb; average 174/Mb per tumor) and the frequency of MMR-deficiency-associated signatures 6, 20, and 26 was comparable to that in LS-OC and LS-CRC. LS-BCs that were pMMR resembled NC-BCs with respect to somatic mutational loads (4/9, 44%, hypermutated with average mutation count 33/Mb vs. 3/10, 30%, hypermutated with average 88 mutations/Mb), whereas mutational signatures shared features of dMMR LS-BC, LS-OC, and LS-CRC. Epigenetic regulatory genes were significantly enriched as mutational targets in LS-BC, LS-OC, and LS-CRC. Many top mutant genes of our LS-BCs have previously been identified as drivers of unselected breast carcinomas. In conclusion, somatic mutational signatures suggest that conventional MMR status of tumor tissues is likely to underestimate the significance of the predisposing MMR defects as contributors to breast tumorigenesis in LS.

JAK/STAT-Activating Genomic Alterations Are a Hallmark of T-PLL.

T-cell prolymphocytic leukemia (T-PLL) is a rare and poor-prognostic mature T-cell leukemia. Recent studies detected genomic aberrations affecting JAK and STAT genes in T-PLL. Due to the limited number of primary patient samples available, genomic analyses of the JAK/STAT pathway have been performed in rather small cohorts. Therefore, we conducted-via a primary-data based pipeline-a meta-analysis that re-evaluated the genomic landscape of T-PLL. It included all available data sets with sequence information on JAK or STAT gene loci in 275 T-PLL. We eliminated overlapping cases and determined a cumulative rate of 62.1% of cases with mutated JAK or STAT genes. Most frequently, JAK1 (6.3%), JAK3 (36.4%), and STAT5B (18.8%) carried somatic single-nucleotide variants (SNVs), with missense mutations in the SH2 or pseudokinase domains as most prevalent. Importantly, these lesions were predominantly subclonal. We did not detect any strong association between mutations of a JAK or STAT gene with clinical characteristics. Irrespective of the presence of gain-of-function (GOF) SNVs, basal phosphorylation of STAT5B was elevated in all analyzed T-PLL. Fittingly, a significant proportion of genes encoding for potential negative regulators of STAT5B showed genomic losses (in 71.4% of T-PLL in total, in 68.4% of T-PLL without any JAK or STAT mutations). They included DUSP4, CD45, TCPTP, SHP1, SOCS1, SOCS3, and HDAC9. Overall, considering such losses of negative regulators and the GOF mutations in JAK and STAT genes, a total of 89.8% of T-PLL revealed a genomic aberration potentially explaining enhanced STAT5B activity. In essence, we present a comprehensive meta-analysis on the highly prevalent genomic lesions that affect genes encoding JAK/STAT signaling components. This provides an overview of possible modes of activation of this pathway in a large cohort of T-PLL. In light of new advances in JAK/STAT inhibitor development, we also outline translational contexts for harnessing active JAK/STAT signaling, which has emerged as a 'secondary' hallmark of T-PLL.

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.

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.

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.

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

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.

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.

Genetic Instability of Influenza pH1N1 Viruses.

Here, we report full-length genome sequences of influenza pH1N1 viruses obtained prior to and after propagation in MDCK cells. Paired comparisons of the genomes showed that each strain acquired 1.0 to 18.8 mutations per genome per replication cycle, which corresponds to 0.5 to 5.8 mutations per virus proteome per replication cycle. Our analysis indicates that pH1N1 viruses accumulated adaptive mutations among others in response to propagation in cell culture. These results could be important for vaccine and drug-sensitivity surveillance studies, as well as for vaccine and antiviral drug development programs where cell cultures are used for influenza propagation.