Machine learning approaches to explore digenic inheritance.

Some rare genetic disorders, such as retinitis pigmentosa or Alport syndrome, are caused by the co-inheritance of DNA variants at two different genetic loci (digenic inheritance). To capture the effects of these disease-causing variants and their possible interactive effects, various statistical methods have been developed in human genetics. Analogous developments have taken place in the field of machine learning, particularly for the field that is now called Big Data. In the past, these two areas have grown independently and have started to converge only in recent years. We discuss an overview of each of the two fields, paying special attention to machine learning methods for uncovering the combined effects of pairs of variants on human disease.

Strategic validation of variants of uncertain significance in ECHS1 genetic testing.

BACKGROUND: Enoyl-CoA hydratase short-chain 1 (ECHS1) is an enzyme involved in the metabolism of branched chain amino acids and fatty acids. Mutations in the ECHS1 gene lead to mitochondrial short-chain enoyl-CoA hydratase 1 deficiency, resulting in the accumulation of intermediates of valine. This is one of the most common causative genes in mitochondrial diseases. While genetic analysis studies have diagnosed numerous cases with ECHS1 variants, the increasing number of variants of uncertain significance (VUS) in genetic diagnosis is a major problem. METHODS: Here, we constructed an assay system to verify VUS function for ECHS1 gene. A high-throughput assay using ECHS1 knockout cells was performed to index these phenotypes by expressing cDNAs containing VUS. In parallel with the VUS validation system, a genetic analysis of samples from patients with mitochondrial disease was performed. The effect on gene expression in cases was verified by RNA-seq and proteome analysis. RESULTS: The functional validation of VUS identified novel variants causing loss of ECHS1 function. The VUS validation system also revealed the effect of the VUS in the compound heterozygous state and provided a new methodology for variant interpretation. Moreover, we performed multiomics analysis and identified a synonymous substitution p.P163= that results in splicing abnormality. The multiomics analysis complemented the diagnosis of some cases that could not be diagnosed by the VUS validation system. CONCLUSIONS: In summary, this study uncovered new ECHS1 cases based on VUS validation and omics analysis; these analyses are applicable to the functional evaluation of other genes associated with mitochondrial disease.

Impact of measuring heteroplasmy of a pathogenic mitochondrial DNA variant at the single-cell level in individuals with mitochondrial disease.

Pathogenic mitochondrial DNA heteroplasmy has mainly been assessed with bulk sequencing in individuals with mitochondrial disease. However, the distribution of heteroplasmy at the single-cell level in skin fibroblasts obtained from individuals, together with detailed clinical and biochemical information, remains to be investigated. We used the mitochondrial DNA single-cell assay for the transposase-accessible chromatin sequencing method. Skin fibroblasts were obtained from six individuals with mitochondrial disease and pathogenic m.3243A>G variants of differing severity. Different distributions of heteroplasmy at the single-cell level were identified in skin fibroblasts from all six individuals. Four individuals with different outcomes showed similar averaged heteroplasmy rates with normal mitochondrial respiratory chain enzyme activity, while the distribution of single-cell heteroplasmy patterns differed among the individuals. This study showed different heteroplasmy distribution patterns at the single-cell level in individuals with the m.3243A>G variant, who had a similar averaged heteroplasmy rates with normal mitochondrial respiratory chain enzyme activity. Whether such different heteroplasmy distribution patterns explain the different clinical outcomes should be assessed further in future studies. Measuring heteroplasmy of pathogenic mitochondrial DNA variants at the single-cell level could be important in individuals with mitochondrial disease.

Population genetics: past, present, and future.

We present selected topics of population genetics and molecular phylogeny. As several excellent review articles have been published and generally focus on European and American scientists, here, we emphasize contributions by Japanese researchers. Our review may also be seen as a belated 50-year celebration of Motoo Kimura's early seminal paper on the molecular clock, published in 1968.

Case Report of Left Ventricular Noncompaction Cardiomyopathy Characterized by Undulating Phenotypes in Adult Patients.

Left ventricular noncompaction cardiomyopathy (LVNC) is a heart muscle disorder morphologically characterized by reticulated trabeculations and intertrabecular recesses in the left ventricular (LV) cavity. LVNC is a genetically and phenotypically heterogeneous condition, which has been increasingly recognized with the accumulation of evidence provided by genotype-phenotype correlation analyses. Here, we report 2 sporadic adult cases of LVNC; both developed acute heart failure as an initial clinical manifestation and harbored causal sarcomere gene mutations. One case was a 57-year-old male with digenic heterozygote mutations, p.R1344Q in myosin heavy chain 7 (MYH7) and p.R144W in troponin T2, cardiac type (TNNT2), who showed morphological characteristics of LVNC in the lateral to apical regions of the LV together with a comorbidity of non-transmural myocardial infarction, resulting from a coronary artery stenosis. After the removal of ischemic insult and standard heart failure treatment, LVNC became less clear, and LV function gradually improved. The other case was a 36-year-old male with a heterozygote mutation, p.E334K in myosin binding protein C3 (MYBPC3), who exhibited cardiogenic shock on admission with morphological characteristics of LVNC being most prominent in the apical segment of the LV. The dosage of beta-blocker was deliberately increased in an outpatient clinic over 6 months following hospitalization, which remarkably improved the LV ejection fraction from 21% to 54.3%. Via a combination of imaging and histopathological and genetic tests, we have found that these cases are not compatible with a persistent phenotype of primary cardiomyopathy, but their morphological features are changeable in response to treatment. Thus, we point out phenotypic plasticity or undulation as a noticeable element of LVNC in this case report.

Mortality of Japanese patients with Leigh syndrome: Effects of age at onset and genetic diagnosis.

Leigh syndrome is a major phenotype of mitochondrial diseases in children. With new therapeutic options being proposed, assessing the mortality and clinical condition of Leigh syndrome patients is crucial for evaluating therapeutics. As data are scarce in Japan, we analysed the mortality rate and clinical condition of Japanese Leigh syndrome patients that we diagnosed since 2007. Data from 166 Japanese patients diagnosed with Leigh syndrome from 2007 to 2017 were reviewed. Patients' present status, method of ventilation and feeding, and degree of disability as of April 2018 was analysed. Overall, 124 (74.7%) were living, 40 (24.1%) were deceased, and 2 (1.2%) were lost to follow-up. Median age of living patients was 8 years (1-39 years). Median length of disease course was 91 months for living patients and 23.5 months for deceased patients. Nearly 90% of deaths occurred by age 6. Mortality rate of patients with onset before 6 months of age was significantly higher than that of onset after 6 months. All patients with neonatal onset were either deceased or bedridden. MT-ATP6 deficiency caused by m.8993T>G mutation and MT-ND5 deficiency induced a severe form of Leigh syndrome. Patients with NDUFAF6, ECHS1, and SURF1 deficiency had relatively mild symptoms and better survival. The impact of onset age on prognosis varied across the genetic diagnoses. The clinical condition of many patients was poor; however, few did not require mechanical ventilation or tube-feeding and were not physically dependent. Early disease onset and genetic diagnosis may have prognostic value.

Heterozygosity mapping for human dominant trait variants.

Homozygosity mapping is a well-known technique to identify runs of homozygous variants that are likely to harbor genes responsible for autosomal recessive disease, but a comparable method for autosomal dominant traits has been lacking. We developed an approach to map dominant disease genes based on heterozygosity frequencies of sequence variants in the immediate vicinity of a dominant trait. We demonstrate through theoretical analysis that DNA variants surrounding an inherited dominant disease variant tend to have increased heterozygosity compared with variants elsewhere in the genome. We confirm existence of this phenomenon in sequence data with known dominant pathogenic variants obtained on family members and in unrelated population controls. A computer-based approach to estimating empirical significance levels associated with our test statistics shows genome-wide p-values smaller than 0.05 for many but not all of the individuals carrying a pathogenic variant.

Barth Syndrome: Different Approaches to Diagnosis.

The diagnosis of Barth syndrome is challenging owing to the wide phenotypic spectrum with allelic heterogeneity. Here we report 3 cases of Barth syndrome with phenotypic and allelic heterogeneity that were diagnosed by different approaches, including whole exome sequencing and final confirmation by reverse-transcription polymease chain reaction.

HDR-del: A tool based on Hamming distance for prioritizing pathogenic chromosomal deletions in exome sequencing.

High-density oligonucleotide arrays have widely been used to detect pathogenic chromosomal deletions. In addition to high-density oligonucleotide arrays, programs using whole-exome sequencing have become available for estimating copy-number variations using depth of coverage. Here, we propose a new statistical method, HDR-del, to prioritize pathogenic chromosomal deletions based on Hamming distance in exome sequencing. In vcf (variant call format) files generated from exome sequencing, hemizygous chromosomal deletion regions lack heterozygous variants and lead to apparent long runs of homozygosity (ROH). In our Hamming distance ratio (HDR)-del approach, we calculate the "difference" in heterozygous status between an affected individual and control individuals using the HDR over all candidate chromosomal deletion regions defined as ROH longer than 1Mbp. Using a suitable test statistic, which is expected to be large for a true pathogenic deletion region, we prioritize candidate chromosomal deletion regions based on this statistic. In our approach, we were able to considerably narrow down true pathogenic chromosomal deletion regions, which were confirmed by high-density oligonucleotide arrays in four mitochondrial disease patients. Our HDR-del approach represents an easy method for detecting chromosomal deletions.