The etiology of autism spectrum disorders (ASD) is very heterogeneous. Mitochondrial dysfunction has been described in ASD; however, primary mitochondrial disease has been genetically proven in a small subset of patients. The main goal of the present study was to investigate correlations between mitochondrial DNA (mtDNA) changes and alterations of genes associated with mtDNA maintenance or ASD. Sixty patients with ASD and sixty healthy individuals were screened for common mtDNA mutations. Next generation sequencing was performed on patients with major mtDNA deletions (mtdel-ASD) using two gene panels to investigate nuclear genes that are associated with ASD or are responsible for mtDNA maintenance. Cohorts of healthy controls, ASD patients without mtDNA alterations, and patients with mitochondrial disorders (non-ASD) harbouring mtDNA deletions served as comparison groups. MtDNA deletions were confirmed in 16.6% (10/60) of patients with ASD (mtdel-ASD). In 90% of this mtdel-ASD children we found rare SNVs in ASD-associated genes (one of those was pathogenic). In the intergenomic panel of this cohort one likely pathogenic variant was present. In patients with mitochondrial disease in genes responsible for mtDNA maintenance pathogenic mutations and variants of uncertain significance (VUS) were detected more frequently than those found in patients from the mtdel-ASD or other comparison groups. In healthy controls and in patients without a mtDNA deletion, only VUS were detected in both panel. MtDNA alterations are more common in patients with ASD than in control individuals. MtDNA deletions are not isolated genetic alterations found in ASD; they coexist either with other ASD-associated genetic risk factors or with alterations in genes responsible for intergenomic communication. These findings indicate that mitochondrial dysfunction is not rare in ASD. The occurring mtDNA deletions in ASD may be mostly a consequence of the alterations of the causative culprit genes for autism or genes responsible for mtDNA maintenance, or because of the harmful effect of environmental factors.

The treatment of acute lymphoblastic leukemia (ALL) and osteosarcoma (OSC) is very effective: the vast majority of patients recover and survive for decades. However, they still need to face serious adverse effects of chemotherapy. One of these is cardiotoxicity which may lead to progressive heart failure in the long term. Cardiotoxicity is contributed mainly to the use of anthracyclines and might have genetic risk factors. Our goal was to test the association between left ventricular function and genetic variations of candidate genes. Echocardiography data from medical records of 622 pediatric ALL and 39 OSC patients were collected from the period 1989–2015. Fractional shortening (FS) and ejection fraction (EF) were determined, 70 single nucleotide polymorphisms (SNPs) in 26 genes were genotyped. Multivariate logistic regression and multi-adjusted general linear model were performed to investigate the influence of genetic polymorphisms on the left ventricular parameters. Bayesian network based Bayesian multilevel analysis of relevance (BN-BMLA) method was applied to test for the potential interaction of the studied cofactors and SNPs. Our results indicate that variations in ABCC2, CYP3A5, NQO1, SLC22A6 and SLC28A3 genes might influence the left ventricular parameters. CYP3A5 rs4646450 TT was 17% among ALL cases with FS lower than 28, and 3% in ALL patients without pathological FS (p = 5.60E-03; OR = 6.94 (1.76–27.39)). SLC28A3 rs7853758 AA was 12% in ALL cases population, while only 1% among controls (p = 6.50E-03; OR = 11.56 (1.98–67.45)). Patients with ABCC2 rs3740066 GG genotype had lower FS during the acute phase of therapy and 5–10 years after treatment (p = 7.38E-03, p = 7.11E-04, respectively). NQO1 rs1043470 rare T allele was associated with lower left ventricular function in the acute phase and 5–10 years after the diagnosis (p = 4.28E-03 and 5.82E-03, respectively), and SLC22A6 gene rs6591722 AA genotype was associated with lower mean FS (p = 1.71E-03), 5–10 years after the diagnosis. Genetic variants in transporters and metabolic enzymes might modulate the individual risk to cardiac toxicity after chemotherapy.

Background: We carried out a candidate gene association study in pediatric acute lymphoblastic leukemia (ALL) to identify possible genetic risk factors in a Hungarian population.Methods: The results were evaluated with traditional statistical methods and with our newly developed Bayesian network based Bayesian multilevel analysis of relevance (BN-BMLA) method. We collected genomic DNA and clinical data from 543 children, who underwent chemotherapy due to ALL, and 529 healthy controls. Altogether 66 single nucleotide polymorphisms (SNPs) in 19 candidate genes were genotyped.Results: With logistic regression, we identified 6 SNPs in theARID5BandIKZF1genes associated with increased risk to B-cell ALL, and two SNPs in theSTAT3gene, which decreased the risk to hyperdiploid ALL. Because the associated SNPs were in linkage in each gene, these associations corresponded to one signal per gene. The odds ratio (OR) associated with the tag SNPs were: OR = 1.69, P = 2.22x10^-7for rs4132601 (IKZF1), OR = 1.53, P = 1.95x10^-5for rs10821936 (ARID5B) and OR = 0.64, P = 2.32x10^-4for rs12949918 (STAT3). With the BN-BMLA we confirmed the findings of the frequentist-based method and received additional information about the nature of the relations between the SNPs and the disease. E.g. the rs10821936 inARID5Band rs17405722 inSTAT3showed a weak interaction, and in case of T-cell lineage sample group, the gender showed a weak interaction with three SNPs in three genes. In the hyperdiploid patient group the BN-BMLA detected a strong interaction among SNPs in theNOTCH1,STAT1,STAT3andBCL2genes. Evaluating the survival rate of the patients with ALL, the BN-BMLA showed that besides risk groups and subtypes, genetic variations in theBAXandCEBPAgenes might also influence the probability of survival of the patients.Conclusions: In the present study we confirmed the roles of genetic variations inARID5BandIKZF1in the susceptibility to B-cell ALL. With the newly developed BN-BMLA method several gene-gene, gene-phenotype and phenotype-phenotype connections were revealed. We showed several advantageous features of the new method, and suggested that in gene association studies the BN-BMLA might be a useful supplementary to the traditional frequentist-based statistical method.

The low concordance between different variant calling methods still poses a challenge for the wide-spread application of next-generation sequencing in research and clinical practice. A wide range of variant annotations can be used for filtering call sets in order to improve the precision of the variant calls, but the choice of the appropriate filtering thresholds is not straightforward. Variant quality score recalibration provides an alternative solution to hard filtering, but it requires large-scale, genomic data. We evaluated germline variant calling pipelines based on BWA and Bowtie 2 aligners in combination with GATK UnifiedGenotyper, GATK HaplotypeCaller, FreeBayes and SAMtools variant callers, using simulated and real benchmark sequencing data (NA12878 with Illumina Platinum Genomes). We argue that these pipelines are not merely discordant, but they extract complementary useful information. We introduce VariantMetaCaller to test the hypothesis that the automated fusion of measurement related information allows better performance than the recommended hard-filtering settings or recalibration and the fusion of the individual call sets without using annotations. VariantMetaCaller uses Support Vector Machines to combine multiple information sources generated by variant calling pipelines and estimates probabilities of variants. This novel method had significantly higher sensitivity and precision than the individual variant callers in all target region sizes, ranging from a few hundred kilobases to whole exomes. We also demonstrated that VariantMetaCaller supports a quantitative, precision based filtering of variants under wider conditions. Specifically, the computed probabilities of the variants can be used to order the variants, and for a given threshold, probabilities can be used to estimate precision. Precision then can be directly translated to the number of true called variants, or equivalently, to the number of false calls, which allows finding problem-specific balance between sensitivity and precision. VariantMetaCaller can be applied to small target regions and whole exomes as well, and it can be used in cases of organisms for which highly accurate variant call sets are not yet available, therefore it can be a viable alternative to hard filtering in cases where variant quality score recalibration cannot be used. VariantMetaCaller is freely available at http://bioinformatics.mit.bme.hu/VariantMetaCaller .

The heterogeneity and complexity of symptom presentation, comorbidities and genetic factors pose challenges to the identification of biological mechanisms underlying complex diseases. Current approaches used to identify biological subtypes of major depressive disorder (MDD) mainly focus on clinical characteristics that cannot be linked to specific biological models. Here, we examined multimorbidities to identify MDD subtypes with distinct genetic and non-genetic factors. We leveraged dynamic Bayesian network approaches to determine a minimal set of multimorbidities relevant to MDD and identified seven clusters of disease-burden trajectories throughout the lifespan among 1.2 million participants from cohorts in the UK, Finland, and Spain. The clusters had clear protective- and risk-factor profiles as well as age-specific clinical courses mainly driven by inflammatory processes, and a comprehensive map of heritability and genetic correlations among these clusters was revealed. Our results can guide the development of personalized treatments for MDD based on the unique genetic, clinical and non-genetic risk-factor profiles of patients.

Background: NKX2-1 related disorders (also known as brain-lung-thyroid syndrome or benign hereditary chorea 1) are associated with a wide spectrum of symptoms. The core features are various movement disorders, characteristically chorea, less frequently myoclonus, dystonia, ataxia; thyroid disease; and lung involvement. The full triad is present in 50% of affected individuals. Numerous additional symptoms may be associated, although many of these were reported only in single cases. Pituitary dysfunction was ambiguously linked to NKX2-1 haploinsufficiency previously. Case Presentation: We examined two members of a family with motor developmental delay, mixed movement disorder (myoclonus, dystonia and chorea) and endocrinological abnormalities (peripheric thyroid disease, and pituitary hormone deficiencies). Dystonia predominated at the father, and myoclonus at the daughter. The father had hypogonadotropic hypogonadism, while the daughter was treated with growth hormone deficiency. Both patients had empty sella on MRI. Candidate gene analyses were negative. Exome sequencing detected a pathogenic stop variation (NM_003317:c.338G>A, p.Trp113*) in the NKX2-1 gene. Conclusions: This case study has two highlights. (1) It draws attention to possible pituitary dysfunction in brain-lung-thyroid syndrome, and provide further evidences that this might be linked to loss of function of the NKX2-1 gene. (2) It underscores the importance of considering NKX2-1 related disorders in the differential diagnosis of myoclonus dystonia.