Clinical utility of regions of homozygosity (ROH) identified in exome sequencing: when to pursue confirmatory uniparental disomy testing for imprinting disorders?

OBJECTIVES

Regions of homozygosity (ROH) could implicate uniparental disomy (UPD) on specific chromosomes associated with imprinting disorders. Though the algorithms for ROH detection in exome sequencing (ES) have been developed, optimal reporting thresholds and when to pursue confirmatory UPD testing for imprinting disorders remain in ambiguity. This study used a data-driven approach to assess optimal reporting thresholds of ROH in clinical practice.

METHODS

ROH analysis was performed using Automap in a retrospective cohort of 8,219 patients and a prospective cohort of 1,964 patients with ES data. Cases with ROH on imprinting-disorders related chromosomes were selected for additional methylation-specific confirmatory testing. The diagnostic yield, the ROH pattern of eventually diagnosed cases and optimal thresholds for confirmatory testing were analyzed.

RESULTS

In the retrospective analysis, 15 true UPD cases of imprinting disorders were confirmed among 51 suspected cases by ROH detection. Pattern of ROH differed between confirmed UPD and non-UPD cases. Maximized yield and minimized false discovery rate of confirmatory UPD testing was achieved at the thresholds of >20 Mb or >25 % chromosomal coverage for interstitial ROH, and >5 Mb for terminal ROH. Current recommendation by ACMG was nearly optimal, though refined thresholds as proposed in this study could reduce the workload by 31 % without losing any true UPD diagnosis. Our refined thresholds remained optimal after independent evaluation in a prospective cohort.

CONCLUSIONS

ROH identified in ES could implicate the presence of clinically relevant UPD. This study recommended size and coverage thresholds for confirmatory UPD testing after ROH detection in ES, contributing to the development of evidence-based reporting guidelines.

C19ORF84 connects piRNA and DNA methylation machineries to defend the mammalian germ line

In the male mouse germ line, PIWI-interacting RNAs (piRNAs), bound by the PIWI protein MIWI2 (PIWIL4), guide DNA methylation of young active transposons through SPOCD1. However, the underlying mechanisms of SPOCD1-mediated piRNA-directed transposon methylation and whether this pathway functions to protect the human germ line remain unknown. We identified loss-of-function variants in human SPOCD1 that cause defective transposon silencing and male infertility. Through the analysis of these pathogenic alleles, we discovered that the uncharacterized protein C19ORF84 interacts with SPOCD1. DNMT3C, the DNA methyltransferase responsible for transposon methylation, associates with SPOCD1 and C19ORF84 in fetal gonocytes. Furthermore, C19ORF84 is essential for piRNA-directed DNA methylation and male mouse fertility. Finally, C19ORF84 mediates the in vivo association of SPOCD1 with the de novo methylation machinery. In summary, we have discovered a conserved role for the human piRNA pathway in transposon silencing and C19ORF84, an uncharacterized protein essential for orchestrating piRNA-directed DNA methylation.

New indicators for monitoring genetic diversity applied to alpine brown trout populations using whole genome sequence data

International policy recently adopted commitments to maintain genetic diversity in wild populations to secure their adaptive potential, including metrics to monitor temporal trends in genetic diversity - so-called indicators. A national programme for assessing trends in genetic diversity was recently initiated in Sweden. Relating to this effort, we systematically assess contemporary genome-wide temporal trends (40 years) in wild populations using the newly adopted indicators and whole genome sequencing (WGS). We use pooled and individual WGS data from brown trout (Salmo trutta) in eight alpine lakes in protected areas. Observed temporal trends in diversity metrics (nucleotide diversity, Watterson's ϴ and heterozygosity) lie within proposed acceptable threshold values for six of the lakes, but with consistently low values in lakes above the tree line and declines observed in these northern-most lakes. Local effective population size is low in all lakes, highlighting the importance of continued protection of interconnected systems to allow genetic connectivity for long-term viability of these populations. Inbreeding (FROH ) spans 10%-30% and is mostly represented by ancient (<1 Mb) runs of homozygosity, with observations of little change in mutational load. We also investigate adaptive dynamics over evolutionarily short time frames (a few generations); identifying putative parallel selection across all lakes within a gene pertaining to skin pigmentation as well as candidates of selection unique to specific lakes and lake systems involved in reproduction and immunity. We demonstrate the utility of WGS for systematic monitoring of natural populations, a priority concern if genetic diversity is to be protected.

Recessive MECR pathogenic variants cause an LHON-like optic neuropathy

BACKGROUND

Leber's hereditary optic neuropathy (LHON) is a mitochondrial disorder characterised by complex I defect leading to sudden degeneration of retinal ganglion cells. Although typically associated with pathogenic variants in mitochondrial DNA, LHON was recently described in patients carrying biallelic variants in nuclear genes DNAJC30, NDUFS2 and MCAT. MCAT is part of mitochondrial fatty acid synthesis (mtFAS), as also MECR, the mitochondrial trans-2-enoyl-CoA reductase. MECR mutations lead to a recessive childhood-onset syndromic disorder with dystonia, optic atrophy and basal ganglia abnormalities.

METHODS

We studied through whole exome sequencing two sisters affected by sudden and painless visual loss at young age, with partial recovery and persistent central scotoma. We modelled the candidate variant in yeast and studied mitochondrial dysfunction in yeast and fibroblasts. We tested protein lipoylation and cell response to oxidative stress in yeast.

RESULTS

Both sisters carried a homozygous pathogenic variant in MECR (p.Arg258Trp). In yeast, the MECR-R258W mutant showed an impaired oxidative growth, 30% reduction in oxygen consumption rate and 80% decrease in protein levels, pointing to structure destabilisation. Fibroblasts confirmed the reduced amount of MECR protein, but failed to reproduce the OXPHOS defect. Respiratory complexes assembly was normal. Finally, the yeast mutant lacked lipoylation of key metabolic enzymes and was more sensitive to H2O2 treatment. Lipoic Acid supplementation partially rescued the growth defect.

CONCLUSION

We report the first family with homozygous MECR variant causing an LHON-like optic neuropathy, which pairs the recent MCAT findings, reinforcing the impairment of mtFAS as novel pathogenic mechanism in LHON.

Functional Characterization of a Novel SLC4A4 Variant and Uniparental Isodisomy in Proximal Renal Tubular Acidosis Patient

SLC4A4 variants are the etiologies of inherited proximal renal tubular acidosis (pRTA), which results in metabolic acidosis, hypokalemia, glaucoma, band keratopathy, and cataract. This study aims to characterize SLC4A4 variant and uniparental isodisomy of chromosome 4 in a patient, and analyse the functional characterization of SLC4A4 variants. This study analyzed renal tubular acidosis disease genes by whole exome sequencing (WES). H3M2 algorithm was used to analyze the run of homozygosity region in chromosomal regions in trio-WES data. The pathogenicity analysis of variants was performed using bioinformatics tools. Additionally, protein stability was analyzed by cycloheximide chase assay. Whole-cell patch clamping was used to examine the electrophysiological properties of NBCe1-A. A novel homozygous SLC4A4 variant was identified in the patient: a missense variant c.496C > T, p. Arg166Trp (NM_003759.4). But the father was heterozygous variant carrier, and the mother did not detect the variant. The H3M2 and UPDio algorithm revealed paternal uniparental isodisomy on chromosome 4 in the patient. SIFT, Poly Phen-2, FATHMM and Mutant Taster showed that the variant might be pathogenic. The tertiary structure analysis showed that the variant could cause structural damage to NBCe1 protein. Foldx results showed that the protein stability of the variant was slightly reduced. Cycloheximide chase assay demonstrated that the variant affects protein stability. The result of electrophysiological studies showed that the variant altered Na+/HCO3- cotransport activity of protein. In conclusion, the study is the first to report a pRTA patient with Arg166Trp variant with UPiD (4) pat and analyze the function of Arg166Trp variant.

Detection of Mosaic Absence of Heterozygosity (AOH) Using Low-Pass Whole Genome Sequencing in Prenatal Diagnosis: A Preliminary Report

Objective: Mosaicism is a common biological phenomenon in organisms and has been reported in many types of chromosome abnormalities, including the absence of heterozygosity (AOH). Due to the detection limitations of the sequencing approach, mosaic AOH events are rarely assessed in clinical cases. Herein, we report the performance of mosaic AOH identification using a low-pass (5~8-fold) WGS method (termed 'CMA-seq', an abbreviation for 'Chromosome Analysis by Sequencing') in fetal genetic diagnosis. Methods: Thirty AOH-negative, eleven constitutional AOH, and three mosaic AOH samples were collected as training data sets to develop the algorithm and evaluate the suitable thresholds for distinguishing mosaic AOH. Twenty-four new chromosomal aberrant cases, along with sixteen constitutional AOH samples, which were previously ascertained via the SNP-array-based method, were used as a validation data set to measure the performance in terms of sensitivity and specificity of this algorithm. Results: A new statistic, 'D-value', was implemented to identify and distinguish constitutional and mosaic AOH events. The reporting thresholds for constitutional and mosaic AOH were also established. In the validation set consisting of 24 new cases, seven constitutional AOH cases and 1 mosaic AOH case were successfully identified, indicating that the results were consistent with those of the SNP-array-based method. The results of all sixteen constitutional AOH validation samples also met the threshold requirements. Conclusions: In this study, we developed a new bioinformatic algorithm to accurately distinguish mosaic AOH from constitutional AOH by low-pass WGS. However, due to the small sample size of the training data set, the algorithm proposed in this manuscript still needs further refinements.

Exome sequencing in 116 patients with inherited thrombocytopenia that remained of unknown origin after systematic phenotype-driven diagnostic workup

Inherited thrombocytopenias (IT) are genetic diseases characterized by low platelet count, sometimes associated with congenital defects or a predisposition to develop additional conditions. Next-generation sequencing has substantially improved our knowledge of IT, with more than 40 genes identified so far, but obtaining a molecular diagnosis remains a challenge especially for patients with non-syndromic forms, having no clinical or functional phenotypes that raise suspicion about specific genes. We performed exome sequencing (ES) in a cohort of 116 IT patients (89 families), still undiagnosed after a previously validated phenotype-driven diagnostic algorithm including a targeted analysis of suspected genes. ES achieved a diagnostic yield of 36%, with a gain of 16% over the diagnostic algorithm. This can be explained by genetic heterogeneity and unspecific genotype-phenotype relationships that make the simultaneous analysis of all the genes, enabled by ES, the most reasonable strategy. Furthermore, ES disentangled situations that had been puzzling because of atypical inheritance, sex-related effects or false negative laboratory results. Finally, ES-based copy number variant analysis disclosed an unexpectedly high prevalence of RUNX1 deletions, predisposing to hematologic malignancies. Our findings demonstrate that ES, including copy number variant analysis, can substantially contribute to the diagnosis of IT and can solve diagnostic problems that would otherwise remain a challenge.

Uniparental disomy for chromosome 1 with POMGNT1 splice-site variant causes muscle-eye-brain disease

POMGNT1, encoding protein O-mannose beta-1,2-N-acetylglucosaminyltransferase 1, is one of the genes responsible for dystroglycanopathy (DGP), which includes multiple phenotypes such as muscle-eye-brain disease (MEB), congenital muscular dystrophy with intellectual disability, and limb-girdle muscular dystrophy Here, we report a case of MEB that is the result of a homozygous variant of POMGNT1 that is revealed through uniparental disomy (UPD). An 8-month-old boy was admitted with mental and motor retardation, hypotonia, esotropia, early onset severe myopia, and structural brain abnormalities. A panel testing of genetic myopathy-related genes was used to identify a homozygous c.636C>T (p.Phe212Phe) variant in exon 7 of POMGNT1 in the patient, a heterozygous c.636C>T variant in the father, and the wild type in the mother. Quantitative polymerase chain reaction (q-PCR) revealed no abnormal copy numbers in exon 7. Trio-based whole-exome sequencing (trio-WES) revealed a possible paternal UPD on chromosome 1 of the patient. Chromosomal microarray analysis (CMA) revealed a 120,451 kb loss of heterozygosity (LOH) on 1p36.33-p11.2, encompassing POMGNT1, and a 99,319 kb loss of heterozygosity on 1q21.2-q44, which indicated UPD. Moreover, RNA sequencing (RNA-seq) verified that the c.636C>T variant was a splice-site variant, leading to skipping of exon 7 (p.Asp179Valfs*23). In conclusion, to the best of our knowledge, we present the first case of MEB caused by UPD, providing valuable insights into the genetic mechanisms underlying this condition.

Diagnosis of Prader-Willi syndrome and Angelman syndrome by targeted nanopore long-read sequencing

The CpG island flanking the promoter region of SNRPN on chromosome 15q11.2 contains CpG sites that are completely methylated in the maternally derived allele and unmethylated in the paternally derived allele. Both unmethylated and methylated alleles are observed in normal individuals. Only the methylated allele is observed in patients with Prader-Willi syndrome, whereas only the unmethylated allele is observed in those with Angelman syndrome. Hence, detection of aberrant methylation at the differentially methylated region is fundamental to the molecular diagnosis of Prader-Willi syndrome and Angelman syndromes. Traditionally, bisulfite treatment and methylation-sensitive restriction enzyme treatment or methylation-specific multiplex ligation-dependent probe amplification (MS-MLPA) have been used. We here developed a long-read sequencing assay that can distinguish methylated and unmethylated CpG sites at 15q11.2 by the difference in current intensity generated from nanopore reads. We successfully diagnosed 4 Prader-Willi syndrome patients and 3 Angelman syndrome patients by targeting differentially methylated regions. Concurrent copy number analysis, homozygosity analysis, and structural variant analysis also allowed us to precisely delineate the underlying pathogenic mechanisms, including gross deletion, uniparental heterodisomy, uniparental isodisomy, or imprinting defect. Furthermore, we showed allele-specific methylation in imprinting-related differentially methylated regions on chromosomes 6, 7, 11, 14, and 20 in a normal individual together with 4 Prader-Willi patients and 3 Angelman syndrome patients. Hence, presently reported method is likely to be applicable to the diagnosis of imprinting disorders other than Prader-Willi syndrome and Angelman syndrome as well.

Diverse monogenic subforms of human spermatogenic failure

Non-obstructive azoospermia (NOA) is the most severe form of male infertility and typically incurable. Defining the genetic basis of NOA has proven challenging, and the most advanced classification of NOA subforms is not based on genetics, but simple description of testis histology. In this study, we exome-sequenced over 1000 clinically diagnosed NOA cases and identified a plausible recessive Mendelian cause in 20%. We find further support for 21 genes in a 2-stage burden test with 2072 cases and 11,587 fertile controls. The disrupted genes are primarily on the autosomes, enriched for undescribed human "knockouts", and, for the most part, have yet to be linked to a Mendelian trait. Integration with single-cell RNA sequencing data shows that azoospermia genes can be grouped into molecular subforms with synchronized expression patterns, and analogs of these subforms exist in mice. This analysis framework identifies groups of genes with known roles in spermatogenesis but also reveals unrecognized subforms, such as a set of genes expressed across mitotic divisions of differentiating spermatogonia. Our findings highlight NOA as an understudied Mendelian disorder and provide a conceptual structure for organizing the complex genetics of male infertility, which may provide a rational basis for disease classification.

Clinical and Genetic Characteristics of a Cohort with Distal Vaginal Atresia

Distal vaginal atresia is a rare abnormality of female reproductive tract in which the vagina is closed or absent. The distal vagina may be replaced by fibrous tissue and the condition is often not diagnosed until a girl fails to begin having periods at puberty. Although it is a congenital disorder, potential genetic causes of distal vaginal atresia are still unknown. We recruited a cohort of 39 patients with distal vaginal atresia and analyzed their phenotypic and genetic features. In addition to the complaint of distal vaginal atresia, approximately 17.9% (7/39) of the patients had other Müllerian anomalies, and 17.9% (7/39) of the patients had other structural abnormalities, including renal-tract, skeletal and cardiac anomalies. Using genome sequencing, we identified two fragment duplications on 17q12 encompassing HNF1B and LHX1, two dosage-sensitive genes with candidate pathogenic variants, in two unrelated patients. A large fragment of uniparental disomy was detected in another patient, affecting genes involved in cell morphogenesis and connective tissue development. Additionally, we reported two variants on TBX3 and AXL, leading to distal vaginal atresia in mutated mouse model, in our clinical subjects for the first time. Essential biological functions of these detected genes with pathogenic variants included regulating reproductive development and cell fate and patterning during embryogenesis. We displayed the comprehensive clinical and genetic characteristic of a cohort with distal vaginal atresia and they were highly heterogeneous both phenotypically and genetically. The duplication of 17q12 in our cohort could help to expand its phenotypic spectrum and potential contribution to the distal vaginal atresia. Our findings of pathogenic genetic variants and associated phenotypes in our cohort could provide evidence and new insight for further research attempting to reveal genetic causes of distal vaginal atresia.

Whole-exome analysis of 177 pediatric patients with undiagnosed diseases

Recently, whole-exome sequencing (WES) has been used for genetic diagnoses of patients who remain otherwise undiagnosed. WES was performed in 177 Japanese patients with undiagnosed conditions who were referred to the Tokai regional branch of the Initiative on Rare and Undiagnosed Diseases (IRUD) (TOKAI-IRUD). This study included only patients who had not previously received genome-wide testing. Review meetings with specialists in various medical fields were held to evaluate the genetic diagnosis in each case, which was based on the guidelines of the American College of Medical Genetics and Genomics. WES identified diagnostic single-nucleotide variants in 66 patients and copy number variants (CNVs) in 11 patients. Additionally, a patient was diagnosed with Angelman syndrome with a complex clinical phenotype upon detection of a paternally derived uniparental disomy (UPD) [upd(15)pat] wherein the patient carried a homozygous DUOX2 p.E520D variant in the UPD region. Functional analysis confirmed that this DUOX2 variant was a loss-of-function missense substitution and the primary cause of congenital hypothyroidism. A significantly higher proportion of genetic diagnoses was achieved compared to previous reports (44%, 78/177 vs. 24-35%, respectively), probably due to detailed discussions and the higher rate of CNV detection.

Genome Analysis in Sick Neonates and Infants: High-yield Phenotypes and Contribution of Small Copy Number Variations

OBJECTIVE

To delineate the diagnostic efficacy of medical exome, whole exome, and whole genome sequencing according to primary symptoms, the contribution of small copy number variations, and the impact of molecular diagnosis on clinical management.

STUDY DESIGN

This was a prospective study of 17 tertiary care centers in Japan, conducted between April 2019 and March 2021. Critically ill neonates and infants less than 6 months of age were recruited in neonatal intensive care units and in outpatient clinics. The patients underwent medical exome, whole exome, or whole genome sequencing as the first tier of testing. Patients with negative results after medical exome or whole exome sequencing subsequently underwent whole genome sequencing. The impact of molecular diagnosis on clinical management was evaluated through contacting primary care physicians.

RESULTS

Of the 85 patients, 41 (48%) had positive results. Based on the primary symptoms, patients with metabolic phenotypes had the highest diagnostic yield (67%, 4/6 patients), followed by renal (60%, 3/5 patients), and neurologic phenotypes (58%, 14/24 patients). Among them, 4 patients had pathogenic small copy number variations identified using whole genome sequencing. In the 41 patients with a molecular diagnosis, 20 (49%) had changes in clinical management.

CONCLUSIONS

Genome analysis for critically ill neonates and infants had a high diagnostic yield for metabolic, renal, and neurologic phenotypes. Small copy number variations detected using whole genome sequencing contributed to the overall molecular diagnosis in 5% of all the patients. The resulting molecular diagnoses had a significant impact on clinical management.

Clinical exome sequencing - mistakes and caveats

Massive parallel sequencing technology has become the predominant technique for genetic diagnostics and research. Many genetic laboratories have wrestled with the challenges of setting up genetic testing workflows based on a completely new technology. The learning curve we went through as a laboratory was accompanied by growing pains while we gained new knowledge and expertise. Here we discuss some important mistakes that have been made in our laboratory through 10 years of clinical exome sequencing but that have given us important new insights on how to adapt our working methods. We provide these examples and the lessons that we learned to help other laboratories avoid to make the same mistakes.

Maternal Uniparental Isodisomy of Chromosome 4 and 8 in Patients with Retinal Dystrophy: SRD5A3-Congenital Disorders of Glycosylation and RP1-Related Retinitis Pigmentosa

Purpose: Uniparental disomy (UPD) is a rare chromosomal abnormality. We performed whole-exosome sequencing (WES) in cases of early-onset retinal dystrophy and identified two cases likely caused by UPD. Herein, we report these two cases and attempt to clarify the clinical picture of retinal dystrophies caused by UPD. Methods: WES analysis was performed for two patients and their parents, who were not consanguineous. Functional analysis was performed in cases suspected of congenital disorders of glycosylation (CDG). We obtained clinical case data and reviewed the literature. Results: In case 1, a novel c.57G>C, p.(Trp19Cys) variant in SRD5A3 was detected homozygously. Genetic analysis suggested a maternal UPD on chromosome 4, and functional analysis confirmed CDG. Clinical findings showed early-onset retinal dystrophy, intellectual disability, and epilepsy. In case 2, an Alu insertion (c.4052_4053ins328, p.[Tyr1352Alafs]) in RP1 was detected homozygously. Maternal UPD on chromosome 8 was suspected. The clinical picture was consistent with RP1-related retinitis pigmentosa. Although the clinical features of retinal dystrophy by UPD may vary, most cases present with childhood onset. Conclusions: There have been limited reports of retinal dystrophy caused by UPD, suggesting that it is rare. Genetic counseling may be encouraged in pediatric cases of retinal dystrophy.

Diagnosis of SLC25A46-related pontocerebellar hypoplasia in two siblings with fulminant neonatal course: role of postmortem CT and whole genomic analysis: a case report

Background

Pontocerebellar hypoplasia (PCH) is increasingly known as a degenerative disease rather than simple “hypoplasia”. At least 21 disease-causing genes have been identified for PCH so far. Because PCH is very heterogenous, prognostic prediction based solely on clinical or radiologic findings is not feasible.

Case presentation

Here, we report two siblings who had a fulminant neonatal course. The documentation of pontocerebellar hypoplasia by postmortem brain CT imaging in one of the siblings and a subsequent complex and comprehensive whole genome analysis established that both siblings had bi-allelic compound heterozygous variants (a splicing variant and a deletion) in the SLC25A46 gene which encodes a solute carrier protein essential for mitochondrial function. Long-read whole genome sequencing was required to confirm the presence of the deletion. The fulminant courses suggest that SLC25A46-related PCH is an acutely progressive degenerative condition starting in utero, rather than a simple static hypoplasia.

Conclusion

The genomic analysis was instrumental and essential to solving the enigma of the unexplained neonatal deaths of these two siblings and to provide accurate genetic counseling.

ROHMM-A flexible hidden Markov model framework to detect runs of homozygosity from genotyping data

Runs of long homozygous (ROH) stretches are considered to be the result of consanguinity and usually contain recessive deleterious disease-causing mutations. Several algorithms have been developed to detect ROHs. Here, we developed a simple alternative strategy by examining X chromosome non-pseudoautosomal region to detect the ROHs from next-generation sequencing data utilizing the genotype probabilities and the hidden Markov model algorithm as a tool, namely ROHMM. It is implemented purely in java and contains both a command line and a graphical user interface. We tested ROHMM on simulated data as well as real population data from the 1000G Project and a clinical sample. Our results have shown that ROHMM can perform robustly producing highly accurate homozygosity estimations under all conditions thereby meeting and even exceeding the performance of its natural competitors.

Charting differentially methylated regions in cancer with Rocker-meth

Differentially DNA methylated regions (DMRs) inform on the role of epigenetic changes in cancer. We present Rocker-meth, a new computational method exploiting a heterogeneous hidden Markov model to detect DMRs across multiple experimental platforms. Through an extensive comparative study, we first demonstrate Rocker-meth excellent performance on synthetic data. Its application to more than 6,000 methylation profiles across 14 tumor types provides a comprehensive catalog of tumor type-specific and shared DMRs, and agnostically identifies cancer-related partially methylated domains (PMD). In depth integrative analysis including orthogonal omics shows the enhanced ability of Rocker-meth in recapitulating known associations, further uncovering the pan-cancer relationship between DNA hypermethylation and transcription factor deregulation depending on the baseline chromatin state. Finally, we demonstrate the utility of the catalog for the study of colorectal cancer single-cell DNA-methylation data.

Matteo Benelli et al. present Rocker-meth, a new Hidden Markov Model (HMM)-based method, to robustly identify differentially methylated regions (DMRs). They use Rocker-meth to analyse more than 6000 methylation profiles across 14 cancer types, providing a catalog of tumor-specific and shared DMRs.

Predicting novel candidate human obesity genes and their site of action by systematic functional screening in Drosophila

The discovery of human obesity-associated genes can reveal new mechanisms to target for weight loss therapy. Genetic studies of obese individuals and the analysis of rare genetic variants can identify novel obesity-associated genes. However, establishing a functional relationship between these candidate genes and adiposity remains a significant challenge. We uncovered a large number of rare homozygous gene variants by exome sequencing of severely obese children, including those from consanguineous families. By assessing the function of these genes in vivo in Drosophila, we identified 4 genes, not previously linked to human obesity, that regulate adiposity (itpr, dachsous, calpA, and sdk). Dachsous is a transmembrane protein upstream of the Hippo signalling pathway. We found that 3 further members of the Hippo pathway, fat, four-jointed, and hippo, also regulate adiposity and that they act in neurons, rather than in adipose tissue (fat body). Screening Hippo pathway genes in larger human cohorts revealed rare variants in TAOK2 associated with human obesity. Knockdown of Drosophila tao increased adiposity in vivo demonstrating the strength of our approach in predicting novel human obesity genes and signalling pathways and their site of action.

Whole-genome resequencing confirms reproductive isolation between sympatric demes of brown trout (Salmo trutta) detected with allozymes

The sympatric existence of genetically distinguishable populations of the same species remains a puzzle in ecology. Coexisting salmonid fish populations are known from over 100 freshwater lakes. Most studies of sympatric populations have used limited numbers of genetic markers making it unclear if genetic divergence involves certain parts of the genome. We returned to the first reported case of salmonid sympatry, initially detected through contrasting homozygosity at a single allozyme locus (coding for lactate dehydrogenase A) in brown trout in the small Lakes Bunnersjöarna, Sweden. First, we verified the existence of the two coexisting demes using a 96-SNP fluidigm array. We then applied whole-genome resequencing of pooled DNA to explore genome-wide diversity within and between these demes; nucleotide diversity was higher in deme I than in deme II. Strong genetic divergence is observed with genome-wide F_ST  ≈ 0.2. Compared with data from populations of similar small lakes, this divergence is of similar magnitude as that between reproductively isolated populations. Individual whole-genome resequencing of two individuals per deme suggests higher inbreeding in deme II versus deme I, indicating different degree of isolation. We located two gene-copies for LDH-A and found divergence between demes in a regulatory section of one of these genes. However, we did not find a perfect fit between the sequence data and previous allozyme results, and this will require further research. Our data demonstrates genome-wide divergence governed mostly by genetic drift but also by diversifying selection in coexisting populations. This type of hidden biodiversity needs consideration in conservation management.

Contribution of uniparental disomy in a clinical trio exome cohort of 2675 patients

Background

Uniparental disomy (UPD) is the inheritance of two homologous chromosomes from the same parent. UPD may result in clinical phenotypes when occurring on chromosomes with specific imprinting pattern, when leading to homozygosity of a deleterious recessive allele inherited from one carrier parent, or when associated with a mosaic aneuploidy. Due to the importance of UPD in genetic disease etiology, UPD analysis has started to be implemented in the context of exome sequencing (ES) or genome sequencing.

Methods

We developed an in‐house algorithm TRIPS (Trio Parentage/UPD Studies) to identify UPD events in trio ES cases. This method identifies regions with uniparental inheritance by utilizing the trio genotyping data obtained from the concurrent SNP array to delineate the parental origin of the SNPs in the proband.

Results

We identified 16 UPD events from 2675 ES trios. Among those, four events led to imprinting disorders, seven unmasked a pathogenic/likely pathogenic variant in a recessive disease gene, and two were consistent with a mosaic genome wide paternal UPD pattern. Twelve of these UPD events directly contributed to the molecular diagnosis of the patients.

Conclusion

Our study demonstrated the contribution of UPD to the molecular diagnosis in one clinical ES cohort, thus UPD analysis should be incorporated into routine clinical ES interpretation.

Technological Improvements in the Genetic Diagnosis of Rett Syndrome Spectrum Disorders

Rett syndrome (RTT) is a severe neurodevelopmental disorder that constitutes the second most common cause of intellectual disability in females worldwide. In the past few years, the advancements in genetic diagnosis brought by next generation sequencing (NGS), have made it possible to identify more than 90 causative genes for RTT and significantly overlapping phenotypes (RTT spectrum disorders). Therefore, the clinical entity known as RTT is evolving towards a spectrum of overlapping phenotypes with great genetic heterogeneity. Hence, simultaneous multiple gene testing and thorough phenotypic characterization are mandatory to achieve a fast and accurate genetic diagnosis. In this review, we revise the evolution of the diagnostic process of RTT spectrum disorders in the past decades, and we discuss the effectiveness of state-of-the-art genetic testing options, such as clinical exome sequencing and whole exome sequencing. Moreover, we introduce recent technological advancements that will very soon contribute to the increase in diagnostic yield in patients with RTT spectrum disorders. Techniques such as whole genome sequencing, integration of data from several “omics”, and mosaicism assessment will provide the tools for the detection and interpretation of genomic variants that will not only increase the diagnostic yield but also widen knowledge about the pathophysiology of these disorders.

Variant PNLDC1, Defective piRNA Processing, and Azoospermia

BACKGROUND

P-element-induced wimpy testis (PIWI)-interacting RNAs (piRNAs) are short (21 to 35 nucleotides in length) and noncoding and are found almost exclusively in germ cells, where they regulate aberrant expression of transposable elements and postmeiotic gene expression. Critical to the processing of piRNAs is the protein poly(A)-specific RNase-like domain containing 1 (PNLDC1), which trims their 3' ends and, when disrupted in mice, causes azoospermia and male infertility.

METHODS

We performed exome sequencing on DNA samples from 924 men who had received a diagnosis of nonobstructive azoospermia. Testicular-biopsy samples were analyzed by means of histologic and immunohistochemical tests, in situ hybridization, reverse-transcriptase-quantitative-polymerase-chain-reaction assay, and small-RNA sequencing.

RESULTS

Four unrelated men of Middle Eastern descent who had nonobstructive azoospermia were found to carry mutations in PNLDC1: the first patient had a biallelic stop-gain mutation, p.R452Ter (rs200629089; minor allele frequency, 0.00004); the second, a novel biallelic missense variant, p.P84S; the third, two compound heterozygous mutations consisting of p.M259T (rs141903829; minor allele frequency, 0.0007) and p.L35PfsTer3 (rs754159168; minor allele frequency, 0.00004); and the fourth, a novel biallelic canonical splice acceptor site variant, c.607-2A→T. Testicular histologic findings consistently showed error-prone meiosis and spermatogenic arrest with round spermatids of type Sa as the most advanced population of germ cells. Gene and protein expression of PNLDC1, as well as the piRNA-processing proteins PIWIL1, PIWIL4, MYBL1, and TDRKH, were greatly diminished in cells of the testes. Furthermore, the length distribution of piRNAs and the number of pachytene piRNAs was significantly altered in men carrying PNLDC1 mutations.

CONCLUSIONS

Our results suggest a direct mechanistic effect of faulty piRNA processing on meiosis and spermatogenesis in men, ultimately leading to male infertility. (Funded by Innovation Fund Denmark and others.).

Loss of function mutations in CCDC32 cause a congenital syndrome characterized by craniofacial, cardiac and neurodevelopmental anomalies

Despite the wide use of genomics to investigate the molecular basis of rare congenital malformations, a significant fraction of patients remains bereft of diagnosis. As part of our continuous effort to recruit and perform genomic and functional studies on such cohorts, we investigated the genetic and mechanistic cause of disease in two independent consanguineous families affected by overlapping craniofacial, cardiac, laterality and neurodevelopmental anomalies. Using whole exome sequencing, we identified homozygous frameshift CCDC32 variants in three affected individuals. Functional analysis in a zebrafish model revealed that ccdc32 depletion recapitulates the human phenotypes. Because some of the patient phenotypes overlap defects common to ciliopathies, we asked if loss of CCDC32 might contribute to the dysfunction of this organelle. Consistent with this hypothesis, we show that ccdc32 is required for normal cilia formation in zebrafish embryos and mammalian cell culture, arguing that ciliary defects are at least partially involved in the pathomechanism of this disorder.

Genome-Wide Detection of Runs of Homozygosity in Laiwu Pigs Revealed by Sequencing Data

Laiwu pigs, distinguished by their high intramuscular fat of 7-9%, is an indigenous pig breed of China, and recent studies also found that Laiwu pigs showed high resistance to Porcine circovirus type 2. However, with the introduction of commercial varieties, the population of Laiwu pigs has declined, and some lineages have disappeared, which could result in inbreeding. Runs of homozygosity (ROH) can be used as a good measure of individual inbreeding status and is also normally used to detect selection signatures so as to map the candidate genes associated with economically important traits. In this study, we used data from Genotyping by Genome Reducing and Sequencing to investigate the number, length, coverage, and distribution patterns of ROH in 93 Chinese Laiwu pigs and identified genomic regions with a high ROH frequency. The average inbreeding coefficient calculated by pedigree was 0.021, whereas that estimated by all detected ROH segments was 0.133. Covering 13.4% of the whole genome, a total of 7,508 ROH segments longer than 1 Mb were detected, whose average length was 3.76 Mb, and short segments (1-5 Mb) dominated. For individuals, the coverage was in the range between 0.56 and 36.86%. For chromosomes, SSC6 had the largest number (n = 688), and the number of ROH in SSC12 was the lowest (n = 215). Thirteen ROH islands were detected in our study, and 86 genes were found within those regions. Some of these genes were correlated with economically important traits, such as meat quality (ECI1, LRP12, NDUFA4L2, GIL1, and LYZ), immunity capacity (IL23A, STAT2, STAT6, TBK1, IFNG, and ITH2), production (DCSTAMP, RDH16, and GDF11), and reproduction (ODF1 and CDK2). A total of six significant Gene Ontology terms and nine significant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were identified, most of which were correlated with disease resistance and biosynthesis processes, and one KEGG pathway was related to lipid metabolism. In addition, we aligned all of the ROH islands to the pig quantitative trait loci (QTL) database and finally found eight QTL related to the intramuscular fat trait. These results may help us understand the characteristics of Laiwu pigs and provide insight for future breeding strategies.

Biallelic variants in TSPOAP1, encoding the active-zone protein RIMBP1, cause autosomal recessive dystonia

Dystonia is a debilitating hyperkinetic movement disorder, which can be transmitted as a monogenic trait. Here, we describe homozygous frameshift, nonsense, and missense variants in TSPOAP1, which encodes the active-zone RIM-binding protein 1 (RIMBP1), as a genetic cause of autosomal recessive dystonia in 7 subjects from 3 unrelated families. Subjects carrying loss-of-function variants presented with juvenile-onset progressive generalized dystonia, associated with intellectual disability and cerebellar atrophy. Conversely, subjects carrying a pathogenic missense variant (p.Gly1808Ser) presented with isolated adult-onset focal dystonia. In mice, complete loss of RIMBP1, known to reduce neurotransmission, led to motor abnormalities reminiscent of dystonia, decreased Purkinje cell dendritic arborization, and reduced numbers of cerebellar synapses. In vitro analysis of the p.Gly1808Ser variant showed larger spike-evoked calcium transients and enhanced neurotransmission, suggesting that RIMBP1-linked dystonia can be caused by either reduced or enhanced rates of spike-evoked release in relevant neural networks. Our findings establish a direct link between dysfunction of the presynaptic active zone and dystonia and highlight the critical role played by well-balanced neurotransmission in motor control and disease pathogenesis.

AutoMap is a high performance homozygosity mapping tool using next-generation sequencing data

Homozygosity mapping is a powerful method for identifying mutations in patients with recessive conditions, especially in consanguineous families or isolated populations. Historically, it has been used in conjunction with genotypes from highly polymorphic markers, such as DNA microsatellites or common SNPs. Traditional software performs rather poorly with data from Whole Exome Sequencing (WES) and Whole Genome Sequencing (WGS), which are now extensively used in medical genetics. We develop AutoMap, a tool that is both web-based or downloadable, to allow performing homozygosity mapping directly on VCF (Variant Call Format) calls from WES or WGS projects. Following a training step on WES data from 26 consanguineous families and a validation procedure on a matched cohort, our method shows higher overall performances when compared with eight existing tools. Most importantly, when tested on real cases with negative molecular diagnosis from an internal set, AutoMap detects three gene-disease and multiple variant-disease associations that were previously unrecognized, projecting clear benefits for both molecular diagnosis and research activities in medical genetics.

Homozygosity mapping is a useful tool for identifying candidate mutations in recessive conditions, however application to next generation sequencing data has been sub-optimal. Here, the authors present AutoMap, which efficiently identifies runs of homozygosity in whole exome/genome sequencing data.

Expanding Phenotype of Schimke Immuno-Osseous Dysplasia: Congenital Anomalies of the Kidneys and of the Urinary Tract and Alteration of NK Cells

Schimke immuno-osseous dysplasia (SIOD) is a rare multisystemic disorder with a variable clinical expressivity caused by biallelic variants in SMARCAL1. A phenotype-genotype correlation has been attempted and variable expressivity of biallelic SMARCAL1 variants may be associated with environmental and genetic disturbances of gene expression. We describe two siblings born from consanguineous parents with a diagnosis of SIOD revealed by whole exome sequencing (WES). Results: A homozygous missense variant in the SMARCAL1 gene (c.1682G>A; p.Arg561His) was identified in both patients. Despite carrying the same variant, the two patients showed substantial renal and immunological phenotypic differences. We describe features not previously associated with SIOD-both patients had congenital anomalies of the kidneys and of the urinary tract and one of them succumbed to a classical type congenital mesoblastic nephroma. We performed an extensive characterization of the immunophenotype showing combined immunodeficiency characterized by a profound lymphopenia, lack of thymic output, defective IL-7Rα expression, and disturbed B plasma cells differentiation and immunoglobulin production in addition to an altered NK-cell phenotype and function. Conclusions: Overall, our results contribute to extending the phenotypic spectrum of features associated with SMARCAL1 mutations and to better characterizing the underlying immunologic disorder with critical implications for therapeutic and management strategies.

Diagnostic yield of additional exome sequencing after the detection of long continuous stretches of homozygosity (LCSH) in SNP arrays

Long continuous stretches of homozygosity (LCSH) are associated with risk of recessive disorders. Though LCSH can be detected by SNP microarrays, additional testing is necessary to clarify the clinical significance. This study is to assess the yield of additional exome sequencing (ES) after LCSH detection and inform the likelihood of eventual diagnosis. In 2226 patients referred to SNP microarrays, 35 patients met the criteria of indicative LCSH. These patients were recruited and went through additional ES. The diagnostic yield was analyzed, and the LCSH pattern was compared between eventually diagnosed cases and those undiagnosed. The results showed additional ES attained a diagnostic yield of 31.4% (11/35), but only one-third of the yield (11.4%, 4/35) was relevant to LCSH. In contrast, two-thirds of the diagnostic variants (20%, 7/35) were de novo or dominantly inherited, irrelevant to the original LCSH finding. No particular LCSH pattern, including the chromosomal coverage or LCSH size, was found to associate with the diagnostic outcome. We concluded that additional ES after LCSH detection could reveal diagnostic variants, but it is strongly recommended to consider all possible inheritance mode, as the diagnostic variants may be irrelevant to the original LCSH finding.

Bi-allelic Mutations in M1AP Are a Frequent Cause of Meiotic Arrest and Severely Impaired Spermatogenesis Leading to Male Infertility

Male infertility affects ∼7% of men, but its causes remain poorly understood. The most severe form is non-obstructive azoospermia (NOA), which is, in part, caused by an arrest at meiosis. So far, only a few validated disease-associated genes have been reported. To address this gap, we performed whole-exome sequencing in 58 men with unexplained meiotic arrest and identified the same homozygous frameshift variant c.676dup (p.Trp226LeufsTer4) in M1AP, encoding meiosis 1 associated protein, in three unrelated men. This variant most likely results in a truncated protein as shown in vitro by heterologous expression of mutant M1AP. Next, we screened four large cohorts of infertile men and identified three additional individuals carrying homozygous c.676dup and three carrying combinations of this and other likely causal variants in M1AP. Moreover, a homozygous missense variant, c.1166C>T (p.Pro389Leu), segregated with infertility in five men from a consanguineous Turkish family. The common phenotype between all affected men was NOA, but occasionally spermatids and rarely a few spermatozoa in the semen were observed. A similar phenotype has been described for mice with disruption of M1ap. Collectively, these findings demonstrate that mutations in M1AP are a relatively frequent cause of autosomal recessive severe spermatogenic failure and male infertility with strong clinical validity.

Improved Diagnosis of Rare Disease Patients through Systematic Detection of Runs of Homozygosity

Autozygosity is associated with an increased risk of genetic rare disease, thus being a relevant factor for clinical genetic studies. More than 2400 exome sequencing data sets were analyzed and screened for autozygosity on the basis of detection of >1 Mbp runs of homozygosity (ROHs). A model was built to predict if an individual is likely to be a consanguineous offspring (accuracy, 98%), and probability of consanguinity ranges were established according to the total ROH size. Application of the model resulted in the reclassification of the consanguinity status of 12% of the patients. The analysis of a subset of 79 consanguineous cases with the Rare Disease (RD)–Connect Genome-Phenome Analysis Platform, combining variant filtering and homozygosity mapping, enabled a 50% reduction in the number of candidate variants and the identification of homozygous pathogenic variants in 41 patients, with an overall diagnostic yield of 52%. The newly defined consanguinity ranges provide, for the first time, specific ROH thresholds to estimate inbreeding within a pedigree on disparate exome sequencing data, enabling confirmation or (re)classification of consanguineous status, hence increasing the efficiency of molecular diagnosis and reporting on secondary consanguinity findings, as recommended by American College of Medical Genetics and Genomics guidelines.

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.

Consecutive medical exome analysis at a tertiary center: Diagnostic and health-economic outcomes

The utility of whole exome analysis has been extensively demonstrated in research settings, but its clinical utility as a first-tier genetic test has not been well documented from diagnostic and health economic standpoints in real-life clinical settings. We performed medical exome analyses focusing on a clinically interpretable portion of the genome (4,813 genes) as a first-tier genetic test for 360 consecutive patients visiting a genetics clinic at a tertiary children's hospital in Japan, over a 3-year period. Bioinformatics analyses were conducted using standard software. A molecular diagnosis was made in 171 patients involving a total of 107 causative genes. Among these 107 causative genes, 57 genes were classified as genes with potential organ-specific interventions and management strategies. Clinically relevant results were obtained in 26% of the total cohort and 54% of the patients with a definitive molecular diagnosis. Performing the medical exome analysis at the time of the initial visit to the tertiary center, rather than after visits to pertinent specialists, brain MRI examination, and G-banded chromosome testing, would have reduced the financial cost by 197 euros according to retrospective calculation under multiple assumption. The present study demonstrated a high diagnostic yield (47.5%) for singleton medical exome analysis as a first-tier test in a real-life setting. Medical exome analysis yielded clinically relevant information in a quarter of the total patient cohort. The application of genomic testing during the initial visit to a tertiary medical center could be a rational approach to the diagnosis of patients with suspected genetic disorders.

A framework for high-resolution phenotyping of candidate male infertility mutants: from human to mouse

Male infertility is a heterogeneous condition of largely unknown etiology that affects at least 7% of men worldwide. Classical genetic approaches and emerging next-generation sequencing studies support genetic variants as a frequent cause of male infertility. Meanwhile, the barriers to transmission of this disease mean that most individual genetic cases will be rare, but because of the large percentage of the genome required for spermatogenesis, the number of distinct causal mutations is potentially large. Identifying bona fide causes of male infertility thus requires advanced filtering techniques to select for high-probability candidates, including the ability to test causality in animal models. The mouse remains the gold standard for defining the genotype-phenotype connection in male fertility. Here, we present a best practice guide consisting of (a) major points to consider when interpreting next-generation sequencing data performed on infertile men, and, (b) a systematic strategy to categorize infertility types and how they relate to human male infertility. Phenotyping infertility in mice can involve investigating the function of multiple cell types across the testis and epididymis, as well as sperm function. These findings will feed into the diagnosis and treatment of male infertility as well as male health broadly.

Bioinformatics and Computational Tools for Next-Generation Sequencing Analysis in Clinical Genetics

Clinical genetics has an important role in the healthcare system to provide a definitive diagnosis for many rare syndromes. It also can have an influence over genetics prevention, disease prognosis and assisting the selection of the best options of care/treatment for patients. Next-generation sequencing (NGS) has transformed clinical genetics making possible to analyze hundreds of genes at an unprecedented speed and at a lower price when comparing to conventional Sanger sequencing. Despite the growing literature concerning NGS in a clinical setting, this review aims to fill the gap that exists among (bio)informaticians, molecular geneticists and clinicians, by presenting a general overview of the NGS technology and workflow. First, we will review the current NGS platforms, focusing on the two main platforms Illumina and Ion Torrent, and discussing the major strong points and weaknesses intrinsic to each platform. Next, the NGS analytical bioinformatic pipelines are dissected, giving some emphasis to the algorithms commonly used to generate process data and to analyze sequence variants. Finally, the main challenges around NGS bioinformatics are placed in perspective for future developments. Even with the huge achievements made in NGS technology and bioinformatics, further improvements in bioinformatic algorithms are still required to deal with complex and genetically heterogeneous disorders.

AUDACITY: A comprehensive approach for the detection and classification of Runs of Homozygosity in medical and population genomics

Runs of Homozygosity (RoHs) are popular among geneticists as the footprint of demographic processes, evolutionary forces and inbreeding in shaping our genome, and are known to confer risk of Mendelian and complex diseases. Notwithstanding growing interest in their study, there is unmet need for reliable and rapid methods for genomic analyses in large data sets. AUDACITY is a tool integrating novel RoH detection algorithm and autozygosity prediction score for prioritization of mutation-surrounding regions. It processes data in VCF file format, and outperforms existing methods in identifying RoHs of any size. Simulations and analysis of real exomes/genomes show its potential to foster future RoH studies in medical and population genomics.

Accurate detection of clinically relevant uniparental disomy from exome sequencing data

Purpose

Uniparental disomy (UPD) is the rare occurrence of two homologous chromosomes originating from the same parent and is typically identified by marker analysis or single-nucleotide polymorphism (SNP)-based microarrays. UPDs may lead to disease due to imprinting effects, underlying homozygous pathogenic variants, or low-level mosaic aneuploidies. In this study we detected clinically relevant UPD events in both trio and single exome sequencing (ES) data.

Methods

UPD was detected by applying a method based on Mendelian inheritance errors to a cohort of 4912 ES trios (all UPD types) and by using median absolute deviation–scaled regions of homozygosity to a cohort of 29,723 single ES samples (isodisomy only).

Results

As positive controls, we accurately identified three mixed UPD, three isodisomy, as well as two segmental UPD events that were all previously reported by SNP-based microarrays. In addition, we identified three segmental UPD and 11 isodisomy events. This resulted in a novel diagnosis based on imprinting for one patient, and adjusted genetic counseling for another patient.

Conclusion

UPD can easily be identified using both single and trio ES and may be clinically relevant to patients. UPD analysis should become routine in clinical ES, because it increases the diagnostic yield and could affect genetic counseling.

Iranome: A catalog of genomic variations in the Iranian population

Considering the application of human genome variation databases in precision medicine, population-specific genome projects are continuously being developed. However, the Middle Eastern population is underrepresented in current databases. Accordingly, we established Iranome database (www.iranome.com) by performing whole exome sequencing on 800 individuals from eight major Iranian ethnic groups representing the second largest population of Middle East. We identified 1,575,702 variants of which 308,311 were novel (19.6%). Also, by presenting higher frequency for 37,384 novel or known rare variants, Iranome database can improve the power of molecular diagnosis. Moreover, attainable clinical information makes this database a good resource for classifying pathogenicity of rare variants. Principal components analysis indicated that, apart from Iranian-Baluchs, Iranian-Turkmen, and Iranian-Persian Gulf Islanders, who form their own clusters, rest of the population were genetically linked, forming a super-population. Furthermore, only 0.6% of novel variants showed counterparts in "Greater Middle East Variome Project", emphasizing the value of Iranome at national level by releasing a comprehensive catalog of Iranian genomic variations and also filling another gap in the catalog of human genome variations at international level. We introduce Iranome as a resource which may also be applicable in other countries located in neighboring regions historically called Greater Iran (Persia).

VCF.Filter: interactive prioritization of disease-linked genetic variants from sequencing data

Next generation sequencing is widely used to link genetic variants to diseases, and it has massively accelerated the diagnosis and characterization of rare genetic diseases. After initial bioinformatic data processing, the interactive analysis of genome, exome, and panel sequencing data typically starts from lists of genetic variants in VCF format. Medical geneticists filter and annotate these lists to identify variants that may be relevant for the disease under investigation, or to select variants that are reported in a clinical diagnostics setting. We developed VCF.Filter to facilitate the search for disease-linked variants, providing a standalone Java program with a user-friendly interface for interactive variant filtering and annotation. VCF.Filter allows the user to define a broad range of filtering criteria through a graphical interface. Common workflows such as trio analysis and cohort-based filtering are pre-configured, and more complex analyses can be performed using VCF.Filter's support for custom annotations and filtering criteria. All filtering is documented in the results file, thus providing traceability of the interactive variant prioritization. VCF.Filter is an open source tool that is freely and openly available at http://vcffilter.rarediseases.at.

POLG mutations presenting as Charcot-Marie-Tooth disease

We report on two patients, with different POLG mutations, in whom axonal neuropathy dominated the clinical picture. One patient presented with late onset sensory axonal neuropathy caused by a homozygous c.2243G>C (p.Trp748Ser) mutation that resulted from uniparental disomy of the long arm of chromosome 15. The other patient had a complex phenotype that included early onset axonal Charcot-Marie-Tooth disease (CMT) caused by compound heterozygous c.926G>A (p.Arg309His) and c.2209G>C (p.Gly737Arg) mutations.

Whole-Exome Sequencing in the Isolated Populations of Cilento from South Italy

The present study describes the genetic architecture of the isolated populations of Cilento, through the analysis of exome sequence data of 245 representative individuals of these populations. By annotating the exome variants and cataloguing them according to their frequency and functional effects, we identified 347,684 variants, 67.4% of which are rare and low frequency variants, and 1% of them (corresponding to 319 variants per person) are classified as high functional impact variants; also, 39,946 (11.5% of the total) are novel variants, for which we determined a significant enrichment for deleterious effects. By comparing the allele frequencies in Cilento with those from the Tuscan population from the 1000 Genomes Project Phase 3, we highlighted an increase in allele frequency in Cilento especially for variants which map to genes involved in extracellular matrix formation and organization. Furthermore, among the variants showing increased frequency we identified several known rare disease-causing variants. By different population genetics analyses, we corroborated the status of the Cilento populations as genetic isolates. Finally, we showed that exome data of Cilento represents a useful local reference panel capable of improving the accuracy of genetic imputation, thus adding power to genetic studies of human traits in these populations.

Mutations in ACTL6B, coding for a subunit of the neuron-specific chromatin remodeling complex nBAF, cause early onset severe developmental and epileptic encephalopathy with brain hypomyelination and cerebellar atrophy

Developmental and epileptic encephalopathies (DEEs) are genetically heterogenous conditions, often characterized by early onset, EEG interictal epileptiform abnormalities, polymorphous and drug-resistant seizures, and neurodevelopmental impairments. In this study, we investigated the genetic defects in two siblings who presented with severe DEE, microcephaly, spastic tetraplegia, diffuse brain hypomyelination, cerebellar atrophy, short stature, and kyphoscoliosis. Whole exome next-generation sequencing (WES) identified in both siblings a homozygous non-sense variant in the ACTL6B gene (NM_016188:c.820C>T;p.Gln274*) coding for a subunit of the neuron-specific chromatin remodeling complex nBAF. To further support these findings, a targeted ACTL6B sequencing assay was performed on a cohort of 85 unrelated DEE individuals, leading to the identification of a homozygous missense variant (NM_016188:c.1045G>A;p.Gly349Ser) in a patient. This variant did not segregate in the unaffected siblings in this family and was classified as deleterious by several prediction softwares. Interestingly, in both families, homozygous patients shared a rather homogeneous phenotype. Very few patients with ACTL6B gene variants have been sporadically reported in WES cohort studies of patients with neurodevelopmental disorders and/or congenital brain malformations. However, the limited number of patients with incomplete clinical information yet reported in the literature did not allow to establish a strong gene-disease association. Here, we provide additional genetic and clinical data on three new cases that support the pathogenic role of ACTL6B gene mutation in a syndromic form of DEE.

Perspectives on the Genomics of HSP Beyond Mendelian Inheritance

Hereditary Spastic Paraplegia is an extraordinarily heterogeneous disease caused by over 50 Mendelian genes. Recent applications of next-generation sequencing, large scale data analysis, and data sharing/matchmaking, have discovered a quickly expanding set of additional HSP genes. Since most recently discovered HSP genes are rare causes of the disease, there is a growing concern of a persisting diagnostic gap, estimated at 30-40%, and even higher for sporadic cases. This missing heritability may not be fully closed by classic Mendelian mutations in protein coding genes. Here we show strategies and published examples of broadening areas of attention for Mendelian and non-Mendelian causes of HSP. We suggest a more inclusive perspective on the potential final architecture of HSP genomics. Efforts to narrow the heritability gap will ultimately lead to more precise and comprehensive genetic diagnoses, which is the starting point for emerging, highly specific gene therapies.

Exome sequencing revealed C1Q homozygous mutation in Pediatric Systemic Lupus Erythematosus

INTRODUCTION AND OBJECTIVES

Pediatric Systemic Lupus Erythematosus (pSLE) is an autoimmune disorder of children. Early disease onset raises the probability of genetic etiology and it is more severe than adult SLE.

PATIENTS AND METHODS

Herein an eight-year-old girl with pSLE from consanguineous parents is reported.

RESULTS

Although she was diagnosed as pSLE since the age of two years, Whole Exome Sequencing (WES) revealed a rare stop-gained C>T mutation in C1QA gene. The variant was validated and segregated in patient and the family. Furthermore, serum levels of the C1q protein were measured and found to be much lower than normal ranges.

CONCLUSIONS

This study indicated that C1Q deficiency should be considered as a differential diagnosis of pSLE. Therefore, measurement of C1q should be recommended in all cases with pSLE.

Biallelic sequence and structural variants in RAX2 are a novel cause for autosomal recessive inherited retinal disease

Purpose

RAX2 encodes a homeobox-containing transcription factor, in which four monoallelic pathogenic variants have been described in autosomal dominant cone-dominated retinal disease.

Methods

Exome sequencing in a European cohort with inherited retinal disease (IRD) (n = 2086) was combined with protein structure modeling of RAX2 missense variants, bioinformatics analysis of deletion breakpoints, haplotyping of RAX2 variant c.335dup, and clinical assessment of biallelic RAX2-positive cases and carrier family members.

Results

Biallelic RAX2 sequence and structural variants were found in five unrelated European index cases, displaying nonsyndromic autosomal recessive retinitis pigmentosa (ARRP) with an age of onset ranging from childhood to the mid-40s (average mid-30s). Protein structure modeling points to loss of function of the novel recessive missense variants and to a dominant-negative effect of the reported dominant RAX2 alleles. Structural variants were fine-mapped to disentangle their underlying mechanisms. Haplotyping of c.335dup in two cases suggests a common ancestry.

Conclusion

This study supports a role for RAX2 as a novel disease gene for recessive IRD, broadening the mutation spectrum from sequence to structural variants and revealing a founder effect. The identification of biallelic RAX2 pathogenic variants in five unrelated families shows that RAX2 loss of function may be a nonnegligible cause of IRD in unsolved ARRP cases.