A recurrent de novo splice site variant involving DNM1 exon 10a causes developmental and epileptic encephalopathy through a dominant-negative mechanism.

Heterozygous pathogenic variants in DNM1 cause developmental and epileptic encephalopathy (DEE) as a result of a dominant-negative mechanism impeding vesicular fission. Thus far, pathogenic variants in DNM1 have been studied with a canonical transcript that includes the alternatively spliced exon 10b. However, after performing RNA sequencing in 39 pediatric brain samples, we find the primary transcript expressed in the brain includes the downstream exon 10a instead. Using this information, we evaluated genotype-phenotype correlations of variants affecting exon 10a and identified a cohort of eleven previously unreported individuals. Eight individuals harbor a recurrent de novo splice site variant, c.1197-8G>A (GenBank: NM_001288739.1), which affects exon 10a and leads to DEE consistent with the classical DNM1 phenotype. We find this splice site variant leads to disease through an unexpected dominant-negative mechanism. Functional testing reveals an in-frame upstream splice acceptor causing insertion of two amino acids predicted to impair oligomerization-dependent activity. This is supported by neuropathological samples showing accumulation of enlarged synaptic vesicles adherent to the plasma membrane consistent with impaired vesicular fission. Two additional individuals with missense variants affecting exon 10a, p.Arg399Trp and p.Gly401Asp, had a similar DEE phenotype. In contrast, one individual with a missense variant affecting exon 10b, p.Pro405Leu, which is less expressed in the brain, had a correspondingly less severe presentation. Thus, we implicate variants affecting exon 10a as causing the severe DEE typically associated with DNM1-related disorders. We highlight the importance of considering relevant isoforms for disease-causing variants as well as the possibility of splice site variants acting through a dominant-negative mechanism.

Clinical genomics and precision medicine.

Precision Medicine emerges from the genomic paradigm of health and disease. For precise molecular diagnoses of genetic diseases, we must analyze the Whole Exome (WES) or the Whole Genome (WGS). By not needing exon capture, WGS is more powerful to detect single nucleotide variants and copy number variants. In healthy individuals, we can observe monogenic highly penetrant variants, which may be causally responsible for diseases, and also susceptibility variants, associated with common polygenic diseases. But there is the major problem of penetrance. Thus, there is the question of whether it is worthwhile to perform WGS in all healthy individuals as a step towards Precision Medicine. The genetic architecture of disease is consistent with the fact that they are all polygenic. Moreover, ancestry adds another layer of complexity. We are now capable of obtaining Polygenic Risk Scores for all complex diseases using only data from new generation sequencing. Yet, review of available evidence does not at present favor the idea that WGS analyses are sufficiently developed to allow reliable predictions of the risk components for monogenic and polygenic hereditary diseases in healthy individuals. Probably, it is still better for WGS to remain reserved for the diagnosis of pathogenic variants of Mendelian diseases.

Assessing the utility of long-read nanopore sequencing for rapid and efficient characterization of mobile element insertions.

Short-read next generation sequencing (NGS) has become the predominant first-line technique used to diagnose patients with rare genetic conditions. Inherent limitations of short-read technology, notably for the detection and characterization of complex insertion-containing variants, are offset by the ability to concurrently screen many disease genes. "Third-generation" long-read sequencers are increasingly being deployed as an orthogonal adjunct technology, but their full potential for molecular genetic diagnosis has yet to be exploited. Here, we describe three diagnostic cases in which pathogenic mobile element insertions were refractory to characterization by short-read sequencing. To validate the accuracy of the long-read technology, we first used Sanger sequencing to confirm the integration sites and derive curated benchmark sequences of the variant-containing alleles. Long-read nanopore sequencing was then performed on locus-specific amplicons. Pairwise comparison between these data and the previously determined benchmark alleles revealed 100% identity of the variant-containing sequences. We demonstrate a number of technical advantages over existing wet-laboratory approaches, including in silico size selection of a mixed pool of amplification products, and the relative ease with which an automated informatics workflow can be established. Our findings add to a growing body of literature describing the diagnostic utility of long-read sequencing.

De novo loss-of-function variants in X-linked MED12 are associated with Hardikar syndrome in females.

PURPOSE: Hardikar syndrome (MIM 612726) is a rare multiple congenital anomaly syndrome characterized by facial clefting, pigmentary retinopathy, biliary anomalies, and intestinal malrotation, but with preserved cognition. Only four patients have been reported previously, and none had a molecular diagnosis. Our objective was to identify the genetic basis of Hardikar syndrome (HS) and expand the phenotypic spectrum of this disorder. METHODS: We performed exome sequencing on two previously reported and five unpublished female patients with a clinical diagnosis of HS. X-chromosome inactivation (XCI) studies were also performed. RESULTS: We report clinical features of HS with previously undescribed phenotypes, including a fatal unprovoked intracranial hemorrhage at age 21. We additionally report the discovery of de novo pathogenic nonsense and frameshift variants in MED12 in these seven individuals and evidence of extremely skewed XCI in all patients with informative testing. CONCLUSION: Pathogenic missense variants in the X-chromosome gene MED12 have previously been associated with Opitz-Kaveggia syndrome, Lujan syndrome, Ohdo syndrome, and nonsyndromic intellectual disability, primarily in males. We propose a fifth, female-specific phenotype for MED12, and suggest that nonsense and frameshift loss-of-function MED12 variants in females cause HS. This expands the MED12-associated phenotype in females beyond intellectual disability.

Structural analysis of new compound heterozygous variants in PEPD gene identified in a patient with Prolidase Deficiency diagnosed by exome sequencing.

Prolidase Deficiency (PD) is an autosomal recessive rare disorder caused by loss or reduction of prolidase enzymatic activity due to variants in the PEPD gene. PD clinical features vary among affected individuals: skin ulcerations, recurrent infections, and developmental delay are common. In this study, we describe a 16-year-old boy with a mild PD phenotype comprising chronic eczema, recurrent infections and elevated IgE. Whole exome sequencing analysis revealed three PEPD variants: c.575T>C p.(Leu192Pro) inherited from the mother, and c.692_694del p.(Tyr231del) and c.1409G>A p.(Arg470His), both inherited from the father. The variant p.(Tyr231del) has been previously characterized by high-resolution X-ray structure analysis as altering protein dynamics/flexibility. In order to study the effects of the other two prolidase variants, we performed site directed mutagenesis purification and crystallization studies. A high-resolution X-ray structure could only be obtained for the p.(Arg470His) variant, which showed no significant structural differences in comparison to WT prolidase. On the other hand, the p.(Leu192Pro) variant led to significant protein destabilization. Hence, we conclude that the maternal p.(Leu192Pro) variant was likely causally associated with the proband´s disease, together with the known pathogenic paternal variant p.(Tyr231del). Our results demonstrated the utility of exome sequencing to perform diagnosis in PD cases with mild phenotype.

Genetic admixture in Brazil.

We review studies from our laboratories using different molecular tools to characterize the Amerindian, European and African ancestry of Brazilians. Initially we used uniparental DNA markers to investigate the contribution of distinct Y chromosome and mitochondrial DNA lineages to present-day populations. High levels of genetic admixture and strong directional mating between European males and Amerindian and African females were unraveled. We next analyzed different types of biparental autosomal polymorphisms. Especially useful was a set of 40 insertion-deletion polymorphisms (indels) that when studied worldwide proved exquisitely sensitive in discriminating between Amerindians, Europeans and Sub-Saharan Africans. When applied to the study of Brazilians these markers confirmed extensive genomic admixture. We then studied ancestry differences in different regions by statistically controlling them to eliminate color considerations. The European ancestry was predominant in all regions studied, with proportions ranging from 60.6% in the Northeast to 77.7% in the South. We propose that the immigration of 6 million Europeans to Brazil in the 19th and 20th centuries is in large part responsible for dissipating previous ancestry dissimilarities that reflected region-specific population histories. Brazilians should be assessed individually, as 210 million human beings, and not as members of specific regions or color groups.

Mendel,MD: A user-friendly open-source web tool for analyzing WES and WGS in the diagnosis of patients with Mendelian disorders.

Whole exome and whole genome sequencing have both become widely adopted methods for investigating and diagnosing human Mendelian disorders. As pangenomic agnostic tests, they are capable of more accurate and agile diagnosis compared to traditional sequencing methods. This article describes new software called Mendel,MD, which combines multiple types of filter options and makes use of regularly updated databases to facilitate exome and genome annotation, the filtering process and the selection of candidate genes and variants for experimental validation and possible diagnosis. This tool offers a user-friendly interface, and leads clinicians through simple steps by limiting the number of candidates to achieve a final diagnosis of a medical genetics case. A useful innovation is the "1-click" method, which enables listing all the relevant variants in genes present at OMIM for perusal by clinicians. Mendel,MD was experimentally validated using clinical cases from the literature and was tested by students at the Universidade Federal de Minas Gerais, at GENE-Núcleo de Genética Médica in Brazil and at the Children's University Hospital in Dublin, Ireland. We show in this article how it can simplify and increase the speed of identifying the culprit mutation in each of the clinical cases that were received for further investigation. Mendel,MD proved to be a reliable web-based tool, being open-source and time efficient for identifying the culprit mutation in different clinical cases of patients with Mendelian Disorders. It is also freely accessible for academic users on the following URL: https://mendelmd.org.

Clinical spectrum and genotype-phenotype associations of KCNA2-related encephalopathies.

Recently, de novo mutations in the gene KCNA2, causing either a dominant-negative loss-of-function or a gain-of-function of the voltage-gated K+ channel Kv1.2, were described to cause a new molecular entity within the epileptic encephalopathies. Here, we report a cohort of 23 patients (eight previously described) with epileptic encephalopathy carrying either novel or known KCNA2 mutations, with the aim to detail the clinical phenotype associated with each of them, to characterize the functional effects of the newly identified mutations, and to assess genotype-phenotype associations. We identified five novel and confirmed six known mutations, three of which recurred in three, five and seven patients, respectively. Ten mutations were missense and one was a truncation mutation; de novo occurrence could be shown in 20 patients. Functional studies using a Xenopus oocyte two-microelectrode voltage clamp system revealed mutations with only loss-of-function effects (mostly dominant-negative current amplitude reduction) in eight patients or only gain-of-function effects (hyperpolarizing shift of voltage-dependent activation, increased amplitude) in nine patients. In six patients, the gain-of-function was diminished by an additional loss-of-function (gain-and loss-of-function) due to a hyperpolarizing shift of voltage-dependent activation combined with either decreased amplitudes or an additional hyperpolarizing shift of the inactivation curve. These electrophysiological findings correlated with distinct phenotypic features. The main differences were (i) predominant focal (loss-of-function) versus generalized (gain-of-function) seizures and corresponding epileptic discharges with prominent sleep activation in most cases with loss-of-function mutations; (ii) more severe epilepsy, developmental problems and ataxia, and atrophy of the cerebellum or even the whole brain in about half of the patients with gain-of-function mutations; and (iii) most severe early-onset phenotypes, occasionally with neonatal onset epilepsy and developmental impairment, as well as generalized and focal seizures and EEG abnormalities for patients with gain- and loss-of-function mutations. Our study thus indicates well represented genotype-phenotype associations between three subgroups of patients with KCNA2 encephalopathy according to the electrophysiological features of the mutations.

Unexplained early onset epileptic encephalopathy: Exome screening and phenotype expansion.

Early onset epileptic encephalopathies (EOEEs) represent a significant diagnostic challenge. Newer genomic approaches have begun to elucidate an increasing number of responsible single genes as well as emerging diagnostic strategies. In this single-center study, we aimed to investigate a cohort of children with unexplained EOEE. We performed whole-exome sequencing (WES), targeting a list of 137 epilepsy-associated genes on 50 children with unexplained EOEE. We characterized all phenotypes in detail and classified children according to known electroclinical syndromes where possible. Infants with previous genetic diagnoses, causative brain malformations, or inborn errors of metabolism were excluded. We identified disease-causing variants in 11 children (22%) in the following genes: STXBP1 (n = 3), KCNB1 (n = 2), KCNT1, SCN1A, SCN2A, GRIN2A, DNM1, and KCNA2. We also identified two further variants (in GRIA3 and CPA6) in two children requiring further investigation. Eleven variants were de novo, and in one paternal testing was not possible. Phenotypes were broadened for some variants identified. This study demonstrates that WES is a clinically useful screening tool for previously investigated unexplained EOEE and allows for reanalysis of data as new genes are being discovered. Detailed phenotyping allows for expansion of specific gene disorders leading to epileptic encephalopathy and emerging sub-phenotypes.

1p13.2 deletion displays clinical features overlapping Noonan syndrome, likely related to NRAS gene haploinsufficiency.

Deletion-induced hemizygosity may unmask deleterious autosomal recessive variants and be a cause of the phenotypic variability observed in microdeletion syndromes. We performed complete exome sequencing (WES) analysis to examine this possibility in a patient with 1p13.2 microdeletion. Since the patient displayed clinical features suggestive of Noonan Syndrome (NS), we also used WES to rule out the presence of pathogenic variants in any of the genes associated with the different types of NS. We concluded that the clinical findings could be attributed solely to the 1p13.2 haploinsufficiency. Retrospective analysis of other nine reported patients with 1p13.2 microdeletions showed that six of them also presented some characteristics of NS. In all these cases, the deleted segment included the NRAS gene. Gain-of-function mutations of NRAS gene are causally related to NS type 6. Thus, it is conceivable that NRAS haploinsufficiency and gain-of-function mutations may have similar clinical consequences. The same phenomenon has been described for two other genes belonging to the Ras/MAPK pathway: MAP2K2 and SHOC2. In conclusion, we here report genotype-phenotype correlations in patients with chromosome 1p13.2 microdeletions and we propose that NRAS may be a critical gene for the NS characteristics in the patients.

Rapid Detection of Rare Deleterious Variants by Next Generation Sequencing with Optional Microarray SNP Genotype Data.

Autozygosity mapping is a powerful technique for the identification of rare, autosomal recessive, disease-causing genes. The ease with which this category of disease gene can be identified has greatly increased through the availability of genome-wide SNP genotyping microarrays and subsequently of exome sequencing. Although these methods have simplified the generation of experimental data, its analysis, particularly when disparate data types must be integrated, remains time consuming. Moreover, the huge volume of sequence variant data generated from next generation sequencing experiments opens up the possibility of using these data instead of microarray genotype data to identify disease loci. To allow these two types of data to be used in an integrated fashion, we have developed AgileVCFMapper, a program that performs both the mapping of disease loci by SNP genotyping and the analysis of potentially deleterious variants using exome sequence variant data, in a single step. This method does not require microarray SNP genotype data, although analysis with a combination of microarray and exome genotype data enables more precise delineation of disease loci, due to superior marker density and distribution.

Simple and efficient identification of rare recessive pathologically important sequence variants from next generation exome sequence data.

Massively parallel ("next generation") DNA sequencing (NGS) has quickly become the method of choice for seeking pathogenic mutations in rare uncharacterized monogenic diseases. Typically, before DNA sequencing, protein-coding regions are enriched from patient genomic DNA, representing either the entire genome ("exome sequencing") or selected mapped candidate loci. Sequence variants, identified as differences between the patient's and the human genome reference sequences, are then filtered according to various quality parameters. Changes are screened against datasets of known polymorphisms, such as dbSNP and the 1000 Genomes Project, in the effort to narrow the list of candidate causative variants. An increasing number of commercial services now offer to both generate and align NGS data to a reference genome. This potentially allows small groups with limited computing infrastructure and informatics skills to utilize this technology. However, the capability to effectively filter and assess sequence variants is still an important bottleneck in the identification of deleterious sequence variants in both research and diagnostic settings. We have developed an approach to this problem comprising a user-friendly suite of programs that can interactively analyze, filter and screen data from enrichment-capture NGS data. These programs ("Agile Suite") are particularly suitable for small-scale gene discovery or for diagnostic analysis.

Tracing the genomic ancestry of Peruvians reveals a major legacy of pre-Columbian ancestors.

In order to investigate the underlying genetic structure and genomic ancestry proportions of Peruvian subpopulations, we analyzed 551 human samples of 25 localities from the Andean, Amazonian, and Coastal regions of Peru with a set of 40 ancestry informative insertion-deletion polymorphisms. Using genotypes of reference populations from different continents for comparison, our analysis indicated that populations from all 25 Peruvian locations had predominantly Amerindian genetic ancestry. Among populations from the Titicaca Lake islands of Taquile, Amantani, Anapia, and Uros, and the Yanque locality from the southern Peruvian Andes, there was no significant proportion of non-autochthonous genomes, indicating that their genetic background is effectively derived from the first settlers of South America. However, the Andean populations from San Marcos, Cajamarca, Characato and Chogo, and coastal populations from Lambayeque and Lima displayed a low but significant European ancestry proportion. Furthermore, Amazonian localities of Pucallpa, Lamas, Chachapoyas, and Andean localities of Ayacucho and Huancayo displayed intermediate levels of non-autochthonous ancestry, mostly from Europe. These results are in close agreement with the documented history of post-Columbian immigrations in Peru and with several reports suggesting a larger effective size of indigenous inhabitants during the formation of the current country's population.

Automated Genotyping of a Highly Informative Panel of 40 Short Insertion-Deletion Polymorphisms Resolved in Polyacrylamide Gels for Forensic Identification and Kinship Analysis.

OBJECTIVE: Short insertion-deletion polymorphisms (indels) are the second most abundant form of genetic variations in humans after SNPs. Since indel alleles differ in size, they can be typed using the same methodological approaches and equipment currently utilized for microsatellite genotyping, which is already operational in forensic laboratories. We have previously shown that a panel of 40 carefully chosen indels has excellent potential for forensic identification, with combined probability of identity (match probability) of 7.09 × 10(-17) for Europeans. METHODS: We describe the successful development of a multiplex system for genotyping the 40-indel panel in long thin denaturing polyacrylamide gels with silver staining. We also demonstrate that the system can be easily fully automated with a simple large scanner and commercial software. RESULTS AND CONCLUSION: The great advantage of the new system of typing is its very low cost. The total price for laboratory equipment is less than EUR 10,000.-, and genotyping of an individual patient will cost less than EUR 10.- in supplies. Thus, the 40-indel panel described here and the newly developed 'low-tech' analysis platform represent useful new tools for forensic identification and kinship analysis in laboratories with limited budgets, especially in developing countries.

Mild Phenotype of Arthrogryposis, Renal Dysfunction, and Cholestasis Syndrome 1 Caused by a Novel VPS33B Variant.

The arthrogryposis, renal dysfunction, and cholestasis syndrome (ARCS) is an autosomal recessive multisystem disease caused by variants in VPS33B or VIPAS39. The classical presentation includes congenital joint contractures, renal tubular dysfunction, cholestasis, and early death. Additional features include ichthyosis, central nervous system malformations, platelet dysfunction, and severe failure to thrive. We studied three patients with cholestasis, increased aminotransferases, normal gamma-glutamyl transferase, and developmental and language delay. Whole exome sequencing analysis identified VPS33B variants in all patients: patients 1 and 2 presented a novel homozygous variant at position c.1148T>A. p.(Ile383Asn), and patient 3 was compound heterozygous for the same c.1148T>A. variant, in addition to the c.940-2A>G. variant. ARCS is compatible with the symptomatology presented by the studied patients. However, most patients that have been described in the literature with ARCS had severe failure to thrive and died in the first 6 months of life. The three patients studied here have a mild ARCS phenotype with prolonged survival. Consequently, we believe that the molecular analysis of the VPS33B and VIPAS39 should be considered in patients with normal gamma-glutamyl transferase cholestasis.

Estimating the degree of identity by descent in consanguineous couples.

In some clinical and research settings, it is often necessary to identify the true level of "identity by descent" (IBD) between two individuals. However, as the individuals become more distantly related, it is increasingly difficult to accurately calculate this value. Consequently, we have developed a computer program that uses genome-wide SNP genotype data from related individuals to estimate the size and extent of IBD in their genomes. In addition, the software can compare a couple's IBD regions with either the autozygous regions of a relative affected by an autosomal recessive disease of unknown cause, or the IBD regions in the parents of the affected relative. It is then possible to calculate the probability of one of the couple's children suffering from the same disease. The software works by finding SNPs that exclude any possible IBD and then identifies regions that lack these SNPs, while exceeding a minimum size and number of SNPs. The accuracy of the algorithm was established by estimating the pairwise IBD between different members of a large pedigree with varying known coefficients of genetic relationship (CGR).

Differential tissue tropism of Trypanosoma cruzi strains: an in vitro study.

We have previously demonstrated selection favoring the JG strain of Trypanosoma cruzi in hearts of BALB/c mice that were chronically infected with an equal mixture of the monoclonal JG strain and a clone of the Colombian strain, Col1.7G2. To evaluate whether cell invasion efficiency drives this selection, we infected primary cultures of BALB/c cardiomyocytes using these same T. cruzi populations. Contrary to expectation, Col1.7G2 parasites invaded heart cell cultures in higher numbers than JG parasites; however, intracellular multiplication of JG parasites was more efficient than that of Col1.7G2 parasites. This phenomenon was only observed for cardiomyocytes and not for cultured Vero cells. Double infections (Col1.7G2 + JG) showed similar results. Even though invasion might influence tissue selection, our data strongly suggest that intracellular development is important to determine parasite tissue tropism.

Trypanosoma cruzi: ancestral genomes and population structure.

Although the genome of Trypanosoma cruzi has been completely sequenced, little is known about its population structure and evolution. Since 1999, two major evolutionary lineages presenting distinct epidemiological characteristics have been recognised: T. cruzi I and T. cruzi II. We describe new and important aspects of the population structure of the parasite, and unequivocally characterise a third ancestral lineage that we propose to name T. cruzi III. Through a careful analysis of haplotypes (blocks of genes that are stably transmitted from generation to generation of the parasite), we inferred at least two hybridisation events between the parental lineages T. cruzi II and T. cruzi III. The strain CL Brener, whose genome was sequenced, is one such hybrid. Based on these results, we propose a simple evolutionary model based on three ancestral genomes, T. cruzi I, T. cruzi II and T. cruzi III. At least two hybridisation events produced evolutionarily viable progeny, and T. cruzi III was the cytoplasmic donor for the resulting offspring (as identified by the mitochondrial clade of the hybrid strains) in both events. This model should be useful to inform evolutionary and pathogenetic hypotheses regarding T. cruzi.

The phylogeography of African Brazilians.

BACKGROUND/AIMS: Approximately four million Africans were taken as slaves to Brazil, where they interbred extensively with Amerindians and Europeans. We have previously shown that while most White Brazilians carry Y chromosomes of European origin, they display high proportions of African and Amerindian mtDNA lineages, because of sex-biased genetic admixture. METHODS: We studied the Y chromosome and mtDNA haplogroup structure of 120 Black males from Sao Paulo, Brazil. RESULTS: Only 48% of the Y chromosomes, but 85% of the mtDNA haplogroups were characteristic of sub-Saharan Africa, confirming our previous observation of sexually biased mating. We mined literature data for mtDNA and Y chromosome haplogroup frequencies for African native populations from regions involved in Atlantic Slave Trade. Principal Components Analysis and Bayesian analysis of population structure revealed no genetic differentiation of Y chromosome marker frequencies between the African regions. However, mtDNA examination unraveled considerable genetic structure, with three clusters at Central-West Africa, West Africa and Southeast Africa. A hypothesis is proposed to explain this structure. CONCLUSION: Using these mtDNA data we could obtain for the first time an estimate of the relative ancestral contribution of Central-West (0.445), West (0.431) and Southeast Africa (0.123) to African Brazilians from Sao Paulo. These estimates are consistent with historical information.