H+-slip correlated to rotor free-wheeling as cause of F1FO-ATPase dysfunction in primary mitochondrial disorders.

Inborn errors of metabolism are related to mitochondrial disorders caused by dysfunction of the oxidative phosphorylation (OXPHOS) system. Congenital hypermetabolism in the infant is a rare disease belonging to Luft syndrome, nonthyroidal hypermetabolism, arising from a singular example of a defect in OXPHOS. The mitochondria lose coupling of mitochondrial substrates oxidation from the ADP phosphorylation. Since Luft syndrome is due to uncoupled cell respiration responsible for deficient in ATP production that originates in the respiratory complexes, a de novo heterozygous variant in the catalytic subunit of mitochondrial F1FO-ATPase arises as the main cause of an autosomal dominant syndrome of hypermetabolism associated with dysfunction in ATP production, which does not involve the respiratory complexes. The F1FO-ATPase works as an embedded molecular machine with a rotary action using two different motor engines. The FO, which is an integral domain in the membrane, dissipates the chemical potential difference for H+, a proton motive force (Δp), across the inner membrane to generate a torsion. The F1 domain-the hydrophilic portion responsible for ATP turnover-is powered by the molecular rotary action to synthesize ATP. The structural and functional coupling of F1 and FO domains support the energy transduction for ATP synthesis. The dissipation of Δp by means of an H+ slip correlated to rotor free-wheeling of the F1FO-ATPase has been discovered to cause enzyme dysfunction in primary mitochondrial disorders. In this insight, we try to offer commentary and analysis of the molecular mechanism in these impaired mitochondria.

The Impairment of Cell Metabolism by Cardiovascular Toxicity of Doxorubicin Is Reversed by Bergamot Polyphenolic Fraction Treatment in Endothelial Cells.

The beneficial effects of bergamot polyphenolic fraction (BPF) on the mitochondrial bioenergetics of porcine aortic endothelial cells (pAECs) were verified under the cardiotoxic action of doxorubicin (DOX). The cell viability of pAECs treated for 24 h with different concentrations of DOX was reduced by 50%, but the negative effect of DOX was reversed in the presence of increasing doses of BPF (100 µg/mL and 200 µg/mL BPF). An analysis of the protective effect of BPF on the toxic action of DOX was also carried out on cell respiration. We observed the inhibition of the mitochondrial activity at 10 µM DOX, which was not restored by 200 µg/mL BPF. Conversely, the decrease in basal respiration and ATP production caused by 0.5 or 1.0 µM DOX were improved in the presence of 100 or 200 µg/mL BPF, respectively. After 24 h of cell recovery with 100 µg/mL or 200 µg/mL BPF on pAECs treated with 0.5 µM or 1.0 µM DOX, respectively, the mitochondrial parameters of oxidative metabolism impaired by DOX were re-boosted.

Victor McKusick and his role in the founding of the European School of Genetic Medicine.

Between 1988 and 2007, during the courses of the European School of Genetic Medicine, many of us had the opportunity to appreciate the tolerant and open-minded personality of Victor McKusick. He was gifted with a unique foresight for the innovations introduced into medicine through the development of the Human Genome Project. The aim of our separate contributions in this article is to document how his insights had an important impact on the European medical training system.

Autozygosity-driven genetic diagnosis in consanguineous families from Italy and the Greater Middle East.

Autozygosity-driven exome analysis has been shown effective for identification of genes underlying recessive diseases especially in countries of the so-called Greater Middle East (GME), where high consanguinity unravels the phenotypic effects of recessive alleles and large family sizes facilitate homozygosity mapping. In Italy, as in most European countries, consanguinity is estimated low. Nonetheless, consanguineous Italian families are not uncommon in publications of genetic findings and are often key to new associations of genes with rare diseases. We collected 52 patients from 47 consanguineous families with suspected recessive diseases, 29 originated in GME countries and 18 of Italian descent. We performed autozygosity-driven exome analysis by detecting long runs of homozygosity (ROHs > 1.5 Mb) and by prioritizing candidate clinical variants within. We identified a pathogenic synonymous variant that had been previously missed in NARS2 and we increased an initial high diagnostic rate (47%) to 55% by matchmaking our candidate genes and including in the analysis shorter ROHs that may also happen to be autozygous. GME and Italian families contributed to diagnostic yield comparably. We found no significant difference either in the extension of the autozygous genome, or in the distribution of candidate clinical variants between GME and Italian families, while we showed that the average autozygous genome was larger and the mean number of candidate clinical variants was significantly higher (p = 0.003) in mutation-positive than in mutation-negative individuals, suggesting that these features influence the likelihood that the disease is autozygosity-related. We highlight the utility of autozygosity-driven genomic analysis also in countries and/or communities, where consanguinity is not widespread cultural tradition.

Associations between psychoactive substance use and sensation seeking behavior among drivers in Norway.

BACKGROUND/AIM: Drug use and risky driving is associated with sensation seeking. The aim of this study was to investigate the association between use of psychoactive substances and levels of the sensation seeking personality trait as measured with the Brief Sensation Seeking Scale 4 among drivers in Norway. METHOD: A cross-sectional design was applied to estimate the association between psychoactive substance use and sensation seeking behavior. Drivers in normal traffic were included in two roadside surveys: one in the north (September 2014 - October 2015) and the other in the south-east of Norway (April 2016 - April 2017). Oral fluid was analyzed for alcohol and psychoactive drugs, and data on sex, age and time of participation were recorded. Participants filled in the Brief Sensation Seeking Scale 4 questionnaire. RESULTS: A total of 8053 drivers were included, of which 32% were women and 62% were under 40 years. The prevalence of alcohol was 0.3%, stimulants 0.6%, tetrahydrocannabinol 1.4%, benzodiazepines and/or z-hypnotics 2.0% and polydrug use 0.6%. Associations were found between the use of tetrahydrocannabinol or benzodiazepines and/or z-hypnotics and a low score on the "thrill and adventure seeking" domain of the Brief Sensation Seeking Scale 4 (OR = 1.723, 95% C.I. = 1.001-2.966). Associations were also found between the use of stimulants and the highest scores on the "experience seeking" (OR = 2.085, 95% C.I. = 1.084-4.009) and "disinhibition" (OR = 4.791, 95% C.I. =1.748-13.135) domains of the Brief Sensation Seeking Scale 4. No associations were found between sensation seeking behavior and alcohol or polydrug use. CONCLUSION: A high degree of sensation seeking was found among drivers who had used stimulating drugs, in contrast to drives who had used tetrahydrocannabinol and benzodiazepines and/or z-hypnotics who showed a low degree of sensation seeking. The combination of sensation seeking behavior and the use of stimulants might lead to increased risky behavior and thus traffic crashes.

Mutant MYO1F alters the mitochondrial network and induces tumor proliferation in thyroid cancer.

Familial aggregation is a significant risk factor for the development of thyroid cancer and familial non-medullary thyroid cancer (FNMTC) accounts for 5-7% of all NMTC. Whole exome sequencing analysis in the family affected by FNMTC with oncocytic features where our group previously identified a predisposing locus on chromosome 19p13.2, revealed a novel heterozygous mutation (c.400G > A, NM_012335; p.Gly134Ser) in exon 5 of MYO1F, mapping to the linkage locus. In the thyroid FRTL-5 cell model stably expressing the mutant MYO1F p.Gly134Ser protein, we observed an altered mitochondrial network, with increased mitochondrial mass and a significant increase in both intracellular and extracellular reactive oxygen species, compared to cells expressing the wild-type (wt) protein or carrying the empty vector. The mutation conferred a significant advantage in colony formation, invasion and anchorage-independent growth. These data were corroborated by in vivo studies in zebrafish, since we demonstrated that the mutant MYO1F p.Gly134Ser, when overexpressed, can induce proliferation in whole vertebrate embryos, compared to the wt one. MYO1F screening in additional 192 FNMTC families identified another variant in exon 7, which leads to exon skipping, and is predicted to alter the ATP-binding domain in MYO1F. Our study identified for the first time a role for MYO1F in NMTC.

A mosaic intragenic microduplication of LAMA1 and a constitutional 18p11.32 microduplication in a patient with keratosis pilaris and intellectual disability.

The application of array-based comparative genomic hybridization and next-generation sequencing has identified many chromosomal microdeletions and microduplications in patients with different pathological phenotypes. Different copy number variations are described within the short arm of chromosome 18 in patients with skin diseases. In particular, full or partial monosomy 18p has also been associated with keratosis pilaris. Here, for the first time, we report a young male patient with intellectual disability, diabetes mellitus (type I), and keratosis pilaris, who exhibited a de novo 45-kb microduplication of exons 4-22 of LAMA1, located at 18p11.31, and a 432-kb 18p11.32 microduplication of paternal origin containing the genes METTL4, NDC80, and CBX3P2 and exons 1-15 of the SMCHD1 gene. The microduplication of LAMA1 was identified in skin fibroblasts but not in lymphocytes, whereas the larger microduplication was present in both tissues. We propose LAMA1 as a novel candidate gene for keratosis pilaris. Although inherited from a healthy father, the 18p11.32 microduplication, which included relevant genes, could also contribute to phenotype manifestation.

Model-based bootstrapping when correcting for measurement error with application to logistic regression.

When fitting regression models, measurement error in any of the predictors typically leads to biased coefficients and incorrect inferences. A plethora of methods have been proposed to correct for this. Obtaining standard errors and confidence intervals using the corrected estimators can be challenging and, in addition, there is concern about remaining bias in the corrected estimators. The bootstrap, which is one option to address these problems, has received limited attention in this context. It has usually been employed by simply resampling observations, which, while suitable in some situations, is not always formally justified. In addition, the simple bootstrap does not allow for estimating bias in non-linear models, including logistic regression. Model-based bootstrapping, which can potentially estimate bias in addition to being robust to the original sampling or whether the measurement error variance is constant or not, has received limited attention. However, it faces challenges that are not present in handling regression models with no measurement error. This article develops new methods for model-based bootstrapping when correcting for measurement error in logistic regression with replicate measures. The methodology is illustrated using two examples, and a series of simulations are carried out to assess and compare the simple and model-based bootstrap methods, as well as other standard methods. While not always perfect, the model-based approaches offer some distinct improvements over the other methods.

Outcome of ankle arthrodesis and ankle prosthesis: a review of the current status.

INTRODUCTION: In advanced stages of ankle osteoarthritis (OA), ankle arthrodesis (AA) or total ankle arthroplasty (TAR) may be necessary. Our purpose is to compare AA and total ankle replacement for the surgical management of end stage ankle OA. SOURCES OF DATA: We conducted a literature search of PubMed, Medline, CINAHL, Cochrane, Embase and Google Scholar databases using the terms 'ankle' in combination with 'OA', 'arthrodesis', 'arthroplasty', 'joint fusion', 'joint replacement'. Studies where treatment was exclusively total ankle replacement or AA were excluded. Treatment characteristics and outcome parameters (overall postoperative outcome and complication rate) were reviewed. AREAS OF AGREEMENT: When counseling patients who are considering their options with regard to ankle arthritis treatment, surgeons should determine on an individual basis which procedure is more suitable. AREAS OF CONTROVERSY: TAR has become an accepted treatment for end-stage OA, but revision rates for TAR are significant higher than for AA (odds ratio 2.28 95% confidence interval [CI], 1.63-3.19; P < 0.0001). GROWING POINTS: The results of TAA are gradually improving, but the procedure cannot yet be recommended for the routine management of ankle OA. AREAS TIMELY FOR DEVELOPING RESEARCH: Although there is some evidence to support TAR to conserve ankle motion and offer improved function and decreased pain with high satisfaction rates, revision rates for TAR are significantly higher than revision rates for AA. Proper patient selection should be better addressed in future studies for successful treatment of end-stage ankle OA. LEVEL OF EVIDENCE: Systematic review, level III.

A novel founder MYO15A frameshift duplication is the major cause of genetic hearing loss in Oman.

The increased risk for autosomal recessive disorders is one of the most well-known medical implications of consanguinity. In the Sultanate of Oman, a country characterized by one of the highest rates of consanguineous marriages worldwide, prevalence of genetic hearing loss (GHL) is estimated to be 6/10 000. Families of GHL patients have higher consanguinity rates than the general Omani population, indicating a major role for recessive forms. Mutations in GJB2, the most commonly mutated GHL gene, have been sporadically described. We collected 97 DNA samples of GHL probands, affected/unaffected siblings and parents from 26 Omani consanguineous families. Analyzing a first family by whole-exome sequencing, we identified a novel homozygous frameshift duplication (c.1171_1177dupGCCATCT) in MYO15A, the gene linked to the deafness locus DFNB3. This duplication was then found in a total of 8/26 (28%) families, within a 849 kb founder haplotype. Reconstruction of haplotype structure at MYO15A surrounding genomic regions indicated that the founder haplotype branched out in the past two to three centuries from a haplotype present worldwide. The MYO15A duplication emerges as the major cause of GHL in Oman. These findings have major implications for the design of GHL diagnosis and prevention policies in Oman.

Intra-individual purifying selection on mitochondrial DNA variants during human oogenesis.

STUDY QUESTION: Does selection for mtDNA mutations occur in human oocytes? SUMMARY ANSWER: We provide statistical evidence in favor of the existence of purifying selection for mtDNA mutations in human oocytes acting between the expulsion of the first and second polar bodies (PBs). WHAT IS KNOWN ALREADY: Several lines of evidence in Metazoa, including humans, indicate that variation within the germline of mitochondrial genomes is under purifying selection. The presence of this internal selection filter in the germline has important consequences for the evolutionary trajectory of mtDNA. However, the nature and localization of this internal filter are still unclear while several hypotheses are proposed in the literature. STUDY DESIGN, SIZE, DURATION: In this study, 60 mitochondrial genomes were sequenced from 17 sets of oocytes, first and second PBs, and peripheral blood taken from nine women between 38 and 43 years of age. PARTICIPANTS/MATERIALS, SETTING, METHODS: Whole genome amplification was performed only on the single cell samples and Sanger sequencing was performed on amplicons. The comparison of variant profiles between first and second PB sequences showed no difference in substitution rates but displayed instead a sharp difference in pathogenicity scores of protein-coding sequences using three different metrics (MutPred, Polyphen and SNPs&GO). MAIN RESULTS AND THE ROLE OF CHANCE: Unlike the first, second PBs showed no significant differences in pathogenic scores with blood and oocyte sequences. This suggests that a filtering mechanism for disadvantageous variants operates during oocyte development between the expulsion of the first and second PB. LARGE SCALE DATA: N/A. LIMITATIONS, REASONS FOR CAUTION: The sample size is small and further studies are needed before this approach can be used in clinical practice. Studies on a model organism would allow the sample size to be increased. WIDER IMPLICATIONS OF THE FINDINGS: This work opens the way to the study of the correlation between mtDNA mutations, mitochondrial capacity and viability of oocytes. STUDY FUNDING/COMPETING INTEREST(S): This work was supported by a SISMER grant. Laboratory facilities and skills were freely provided by SISMER, and by the Alma Mater Studiorum, University of Bologna. The authors have no conflict of interest to disclose.

A mutation in PAK3 with a dual molecular effect deregulates the RAS/MAPK pathway and drives an X-linked syndromic phenotype.

Loss-of-function mutations in PAK3 contribute to non-syndromic X-linked intellectual disability (NS-XLID) by affecting dendritic spine density and morphology. Linkage analysis in a three-generation family with affected males showing ID, agenesis of corpus callosum, cerebellar hypoplasia, microcephaly and ichthyosis, revealed a candidate disease locus in Xq21.33q24 encompassing over 280 genes. Subsequent to sequencing all coding exons of the X chromosome, we identified a single novel variant within the linkage region, affecting a conserved codon of PAK3. Biochemical studies showed that, similar to previous NS-XLID-associated lesions, the predicted amino acid substitution (Lys389Asn) abolished the kinase activity of PAK3. In addition, the introduced residue conferred a dominant-negative function to the protein that drives the syndromic phenotype. Using a combination of in vitro and in vivo studies in zebrafish embryos, we show that PAK3(N389) escapes its physiologic degradation and is able to perturb MAPK signaling via an uncontrolled kinase-independent function, which in turn leads to alterations of cerebral and craniofacial structures in vivo. Our data expand the spectrum of phenotypes associated with PAK3 mutations, characterize a novel mechanism resulting in a dual molecular effect of the same mutation with a complex PAK3 functional deregulation and provide evidence for a direct functional impact of aberrant PAK3 function on MAPK signaling.

Mutations in RAD21 disrupt regulation of APOB in patients with chronic intestinal pseudo-obstruction.

BACKGROUND & AIMS: Chronic intestinal pseudo-obstruction (CIPO) is characterized by severe intestinal dysmotility that mimics a mechanical subocclusion with no evidence of gut obstruction. We searched for genetic variants associated with CIPO to increase our understanding of its pathogenesis and to identify potential biomarkers. METHODS: We performed whole-exome sequencing of genomic DNA from patients with familial CIPO syndrome. Blood and lymphoblastoid cells were collected from patients and controls (individuals without CIPO); levels of messenger RNA (mRNA) and proteins were analyzed by quantitative reverse-transcription polymerase chain reaction, immunoblot, and mobility shift assays. Complementary DNAs were transfected into HEK293 cells. Expression of rad21 was suppressed in zebrafish embryos using a splice-blocking morpholino (rad21a). Gut tissues were collected and analyzed. RESULTS: We identified a homozygous mutation (p.622, encodes Ala>Thr) in RAD21 in patients from a consanguineous family with CIPO. Expression of RUNX1, a target of RAD21, was reduced in cells from patients with CIPO compared with controls. In zebrafish, suppression of rad21a reduced expression of runx1; this phenotype was corrected by injection of human RAD21 mRNA, but not with the mRNA from the mutated p.622 allele. rad21a Morpholino zebrafish had delayed intestinal transit and greatly reduced numbers of enteric neurons, similar to patients with CIPO. This defect was greater in zebrafish with suppressed expression of ret and rad21, indicating their interaction in the regulation of gut neurogenesis. The promoter region of APOB bound RAD21 but not RAD21 p.622 Ala>Thr; expression of wild-type RAD21 in HEK293 cells repressed expression of APOB, compared with control vector. The gut-specific isoform of APOB (APOB48) is overexpressed in sera from patients with CIPO who carry the RAD21 mutation. APOB48 also is overexpressed in sporadic CIPO in sera and gut biopsy specimens. CONCLUSIONS: Some patients with CIPO carry mutations in RAD21 that disrupt the ability of its product to regulate genes such as RUNX1 and APOB. Reduced expression of rad21 in zebrafish, and dysregulation of these target genes, disrupts intestinal transit and the development of enteric neurons.

Validation of the Italian version of the Oxford Ankle Foot Questionnaire for children.

PURPOSE: The purpose of this study was to translate the Oxford Ankle Foot Questionnaire (OAFQ) into Italian, to perform a cross-cultural adaptation and to evaluate its psychometric properties. METHODS: The Italian OAFQ was developed according to the recommended forward/backward translation protocol and evaluated in pediatric patients treated for symptomatic flatfoot deformity. Feasibility, reliability, internal consistency, construct validity [comparing OAFQ domains with Child Health Questionnaire (CHQ) domains] and responsiveness to surgical treatment were assessed. RESULTS: A total of 61 children and their parents were enrolled in the study. Results showed satisfactory levels of internal consistency for both children and parent forms. The test-retest reliability was confirmed by high ICC values for both child and parents subscales. Good construct validity was showed by patterns of relationships consistent with theoretically related domains of the CHQ. After surgery, the mean OAFQ scores improved in all the domains after treatment with the subtalar arthroereisis, for both children and parent scales (p < 0.01). Effect size ranged from small to moderate for almost all domains. CONCLUSIONS: The Italian version of the OAFQ might be a reliable and valid instrument in order to evaluate interventions used to treat children's foot or ankle problem, but needs further study on different clinical settings.

Positive selection of lactase persistence among people of Southern Arabia.

OBJECTIVE: Frequency patterns of the lactase persistence (LP)-associated -13,915 G allele and archaeological records pointing to substantial role played by southern regions in the peopling and domestication processes that involved the Arabian Peninsula suggest that Southern Arabia plausibly represented the center of diffusion of such adaptive variant. Nevertheless, a well-defined scenario for evolution of Arabian LP is still to be elucidated and the burgeoning archaeological picture of complex human migrations occurred through the peninsula is not matched by an equivalent high-resolution description of genetic variation underlying this adaptive trait. To fill this gap, we investigated diversity at a wide genomic interval surrounding the LCT gene in different Southern Arabian populations. METHODS: 40 SNPs were genotyped to characterize LCT profiles of 630 Omani and Yemeni individuals to perform population structure, linkage disequilibrium, population differentiation-based and haplotype-based analyses. RESULTS: Typical Arabian LP-related variation was found in Dhofaris and Yemenis, being characterized by private haplotypes carrying the -13,915 G allele, unusual differentiation with respect to northern groups and conserved homozygous haplotype-blocks, suggesting that the adaptive allele was likely introduced in the Arabian gene pool in southern populations and was then subjected to prolonged selective pressure. CONCLUSION: By pointing to Yemen as one of the best candidate centers of diffusion of the Arabian-specific adaptive variant, obtained results indicate the spread of indigenous groups as the main process underlying dispersal of LP along the Arabian Peninsula, supporting a refugia model for Arabian demic movements occurred during the Terminal Pleistocene and Early Holocene.

Characterization and identification of hidden rare variants in the human genome.

BACKGROUND: By examining the genotype calls generated by the 1000 Genomes Project we discovered that the human reference genome GRCh37 contains almost 20,000 loci in which the reference allele has never been observed in healthy individuals and around 70,000 loci in which it has been observed only in the heterozygous state. RESULTS: We show that a large fraction of this rare reference allele (RRA) loci belongs to coding, functional and regulatory elements of the genome and could be linked to rare Mendelian disorders as well as cancer. We also demonstrate that classical germline and somatic variant calling tools are not capable to recognize the rare allele when present in these loci. To overcome such limitations, we developed a novel tool, named RAREVATOR, that is able to identify and call the rare allele in these genomic positions. By using a small cancer dataset we compared our tool with two state-of-the-art callers and we found that RAREVATOR identified more than 1,500 germline and 22 somatic RRA variants missed by the two methods and which belong to significantly mutated pathways. CONCLUSIONS: These results show that, to date, the investigation of around 100,000 loci of the human genome has been missed by re-sequencing experiments based on the GRCh37 assembly and that our tool can fill the gap left by other methods. Moreover, the investigation of the latest version of the human reference genome, GRCh38, showed that although the GRC corrected almost all insertions and a small part of SNVs and deletions, a large number of functionally relevant RRAs still remain unchanged. For this reason, also future resequencing experiments, based on GRCh38, will benefit from RAREVATOR analysis results. RAREVATOR is freely available at http://sourceforge.net/projects/rarevator .

Mitochondrial DNA analysis and numerical chromosome condition in human oocytes and polar bodies.

To investigate the mitochondrial DNA (mtDNA) segregation in human oocytes, the level of heteroplasmy in the three products of meioses, polar bodies (PBs) and corresponding oocytes, was assessed by studying the hypervariable region I (HVRI) of the D-loop region. The DNA from 122 PBs and 51 oocytes from 16 patients was amplified by whole genome amplification (WGA). An aliquot of the WGA product was used to assess aneuploidy, and another aliquot to study mtDNA. The HVRI was amplified and sequenced with an efficiency of 75.4 and 63%, respectively, in PBs, and of 100% in oocytes. The comparison with the mtDNA sequences from blood of the individual donors showed full correspondence of polymorphisms with the matching oocytes, whilst in PBs the degree of concordance dropped to 89.6%. Haplogroups were inferred for all 16 patients. Of the 89 diagnosed PBs from the 13 patients belonging to macrohaplogroup R, 23 were euploid and 66 aneuploid. The incidence of total anomalies was significantly lower in haplogroup H (6.5%) when compared with haplogroups J and T (17.6 and 13.4% respectively; P < 0.001). In haplogroup J, hypoaneuploidy occurred more frequently than hyperaneuploidy. In the three patients belonging to haplogroup N*, 81% of PBs were aneuploid with similar rates of chromosome hypoaneuploidy and hyperaneuploidy. The presence of mtDNA base changes confined to PBs could reflect a selection mechanism against severe mtDNA mutations, while permitting a high evolution rate that could result in bioenergetic diversity. The different susceptibility to aneuploidy by some haplogroups strongly supports this hypothesis.

H3M2: detection of runs of homozygosity from whole-exome sequencing data.

MOTIVATION: Runs of homozygosity (ROH) are sizable chromosomal stretches of homozygous genotypes, ranging in length from tens of kilobases to megabases. ROHs can be relevant for population and medical genetics, playing a role in predisposition to both rare and common disorders. ROHs are commonly detected by single nucleotide polymorphism (SNP) microarrays, but attempts have been made to use whole-exome sequencing (WES) data. Currently available methods developed for the analysis of uniformly spaced SNP-array maps do not fit easily to the analysis of the sparse and non-uniform distribution of the WES target design. RESULTS: To meet the need of an approach specifically tailored to WES data, we developed [Formula: see text], an original algorithm based on heterogeneous hidden Markov model that incorporates inter-marker distances to detect ROH from WES data. We evaluated the performance of [Formula: see text] to correctly identify ROHs on synthetic chromosomes and examined its accuracy in detecting ROHs of different length (short, medium and long) from real 1000 genomes project data. [Formula: see text] turned out to be more accurate than GERMLINE and PLINK, two state-of-the-art algorithms, especially in the detection of short and medium ROHs. AVAILABILITY AND IMPLEMENTATION: [Formula: see text] is a collection of bash, R and Fortran scripts and codes and is freely available at https://sourceforge.net/projects/h3m2/. CONTACT: albertomagi@gmail.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A mutation screening of oncogenes, tumor suppressor gene TP53 and nuclear encoded mitochondrial complex I genes in oncocytic thyroid tumors.

BACKGROUND: Thyroid neoplasias with oncocytic features represent a specific phenotype in non-medullary thyroid cancer, reflecting the unique biological phenomenon of mitochondrial hyperplasia in the cytoplasm. Oncocytic thyroid cells are characterized by a prominent eosinophilia (or oxyphilia) caused by mitochondrial abundance. Although disruptive mutations in the mitochondrial DNA (mtDNA) are the most significant hallmark of such tumors, oncocytomas may be envisioned as heterogeneous neoplasms, characterized by multiple nuclear and mitochondrial gene lesions. We investigated the nuclear mutational profile of oncocytic tumors to pinpoint the mutations that may trigger the early oncogenic hit. METHODS: Total DNA was extracted from paraffin-embedded tissues from 45 biopsies of oncocytic tumors. High-resolution melting was used for mutation screening of mitochondrial complex I subunits genes. Specific nuclear rearrangements were investigated by RT-PCR (RET/PTC) or on isolated nuclei by interphase FISH (PAX8/PPARγ). Recurrent point mutations were analyzed by direct sequencing. RESULTS: In our oncocytic tumor samples, we identified rare TP53 mutations. The series of analyzed cases did not include poorly- or undifferentiated thyroid carcinomas, and none of the TP53 mutated cases had significant mitotic activity or high-grade features. Thus, the presence of disruptive TP53 mutations was completely unexpected. In addition, novel mutations in nuclear-encoded complex I genes were identified. CONCLUSIONS: These findings suggest that nuclear genetic lesions altering the bioenergetics competence of thyroid cells may give rise to an aberrant mitochondria-centered compensatory mechanism and ultimately to the oncocytic phenotype.