Validation of plasma protein glycation and oxidation biomarkers for the diagnosis of autism.

Autism Spectrum Disorder (ASD) is a common neurodevelopmental disorder in children. It is currently diagnosed by behaviour-based assessments made by observation and interview. In 2018 we reported a discovery study of a blood biomarker diagnostic test for ASD based on a combination of four plasma protein glycation and oxidation adducts. The test had 88% accuracy in children 5-12 years old. Herein, we present an international multicenter clinical validation study (N = 478) with application of similar biomarkers to a wider age range of 1.5-12 years old children. Three hundred and eleven children with ASD (247 male, 64 female; age 5.2 ± 3.0 years) and 167 children with typical development (94 male, 73 female; 4.9 ± 2.4 years) were recruited for this study at Sidra Medicine and Hamad Medical Corporation hospitals, Qatar, and Hospital Regional Universitario de Málaga, Spain. For subjects 5-12 years old, the diagnostic algorithm with features, advanced glycation endproducts (AGEs)-Nε-carboxymethyl-lysine (CML), Nω-carboxymethylarginine (CMA) and 3-deoxyglucosone-derived hydroimidazolone (3DG-H), and oxidative damage marker, o,o'-dityrosine (DT), age and gender had accuracy 83% (CI 79 - 89%), sensitivity 94% (CI 90-98%), specificity 67% (CI 57-76%) and area-under-the-curve of receiver operating characteristic plot (AUROC) 0.87 (CI 0.84-0.90). Inclusion of additional plasma protein glycation and oxidation adducts increased the specificity to 74%. An algorithm with 12 plasma protein glycation and oxidation adduct features was optimum for children of 1.5-12 years old: accuracy 74% (CI 70-79%), sensitivity 75% (CI 63-87%), specificity 74% (CI 58-90%) and AUROC 0.79 (CI 0.74-0.84). We conclude that ASD diagnosis may be supported using an algorithm with features of plasma protein CML, CMA, 3DG-H and DT in 5-12 years-old children, and an algorithm with additional features applicable for ASD screening in younger children. ASD severity, as assessed by ADOS-2 score, correlated positively with plasma protein glycation adducts derived from methylglyoxal, hydroimidazolone MG-H1 and Nε(1-carboxyethyl)lysine (CEL). The successful validation herein may indicate that the algorithm modifiable features are mechanistic risk markers linking ASD to increased lipid peroxidation, neuronal plasticity and proteotoxic stress.

Analysis of incidental findings in Qatar genome participants reveals novel functional variants in LMNA and DSP.

In order to report clinically actionable incidental findings in genetic testing, the American College of Medical Genetics and Genomics (ACMG) recommended the evaluation of variants in 59 genes associated with highly penetrant mutations. However, there is a lack of epidemiological data on medically actionable rare variants in these genes in Arab populations. We used whole genome sequencing data from 6045 participants from the Qatar Genome Programme and integrated it with phenotypic data collected by the Qatar Biobank. We identified novel putative pathogenic variants in the 59 ACMG genes by filtering previously unrecorded variants based on computational prediction of pathogenicity, variant rarity and segregation evidence. We assessed the phenotypic associations of candidate variants in genes linked to cardiovascular diseases. Finally, we used a zebrafish knockdown and synthetic human mRNA co-injection assay to functionally characterize two of these novel variants. We assessed the zebrafish cardiac function in terms of heart rate, rhythm and hemodynamics, as well as the heart structure. We identified 52 492 novel variants, which have not been reported in global and disease-specific databases. A total of 74 novel variants were selected with potentially pathogenic effect. We prioritized two novel cardiovascular variants, DSP c.1841A > G (p.Asp614Gly) and LMNA c.326 T > G (p.Val109Gly) for functional characterization. Our results showed that both variants resulted in abnormal zebrafish heart rate, rhythm and structure. This study highlights medically actionable variants that are specific to the Middle Eastern Qatari population.

A variant in sperm-specific glycolytic enzyme enolase 4 (ENO4) causes human male infertility.

BACKGROUND: Although defects in sperm morphology and physiology lead to male infertility, in many instances, the exact disruption of molecular pathways in a given patient is often unknown. The glycolytic pathway is an essential process to supply energy in sperm cell motility. Enolase 4 (ENO4) is crucial for the glycolytic process, which provides the energy for sperm cells in motility. ENO4 is located in the sperm principal piece and is essential for the motility and organization of the sperm flagellum. In the present study, we characterized a family with asthenozoospermia and abnormal sperm morphology as a result of a variant in the enolase 4 (ENO4) gene. METHODS: Computer-assisted semen analysis, papanicolaou smear staining and scanning electron microscopy were used to examine sperm motility and morphology for semen analysis in patients. For genetic analysis, whole-exome sequencing followed by Sanger sequencing was performed. RESULTS: Two brothers in a consanguineous family were being clinically investigated for sperm motility and morphology issues. Genetic analysis by whole-exome sequencing revealed a homozygous variant [c.293A>G, p.(Lys98Arg)] in the ENO4 gene that segregated with infertility in the family, shared by affected but not controls. CONCLUSIONS: In view of the association of asthenozoospermia and abnormal sperm morphology in Eno4 knockout mice, we consider this to be the first report describing the involvement of ENO4 gene in human male infertility. We also explore the possible involvement of another variant in explaining other phenotypic features in this family.

Conserved genes regulating human sex differentiation, gametogenesis and fertilization.

The study of the functional genome in mice and humans has been instrumental for describing the conserved molecular mechanisms regulating human reproductive biology, and for defining the etiologies of monogenic fertility disorders. Infertility is a reproductive disorder that includes various conditions affecting a couple's ability to achieve a healthy pregnancy. Recent advances in next-generation sequencing and CRISPR/Cas-mediated genome editing technologies have facilitated the identification and characterization of genes and mechanisms that, if affected, lead to infertility. We report established genes that regulate conserved functions in fundamental reproductive processes (e.g., sex determination, gametogenesis, and fertilization). We only cover genes the deletion of which yields comparable fertility phenotypes in both rodents and humans. In the case of newly-discovered genes, we report the studies demonstrating shared cellular and fertility phenotypes resulting from loss-of-function mutations in both species. Finally, we introduce new model systems for the study of human reproductive biology and highlight the importance of studying human consanguineous populations to discover novel monogenic causes of infertility. The rapid and continuous screening and identification of putative genetic defects coupled with an efficient functional characterization in animal models can reveal novel mechanisms of gene function in human reproductive tissues.

A loss-of-function AGTR1 variant in a critically-ill infant with renal tubular dysgenesis: case presentation and literature review.

BACKGROUND: Renal tubular dysgenesis (RTD) is a severe disorder with poor prognosis significantly impacting the proximal tubules of the kidney while maintaining an anatomically normal gross structure. The genetic origin of RTD, involving variants in the ACE, REN, AGT, and AGTR1 genes, affects various enzymes or receptors within the Renin angiotensin system (RAS). This condition manifests prenatally with oligohydramninos and postnatally with persistent anuria, severe refractory hypotension, and defects in skull ossification. CASE PRESENTATION: In this report, we describe a case of a female patient who, despite receiving multi vasopressor treatment, experienced persistent hypotension, ultimately resulting in early death at five days of age. While there was a history of parental consanguinity, no reported family history of renal disease existed. Blood samples from the parents and the remaining DNA sample of the patient underwent Whole Genome Sequencing (WGS). The genetic analysis revealed a rare homozygous loss of function variant (NM_000685.5; c.415C > T; p.Arg139*) in the Angiotensin II Receptor Type 1 (AGTR1) gene. CONCLUSION: This case highlights the consequence of loss-of-function variants in AGTR1 gene leading to RTD, which is characterized by high mortality rate at birth or during the neonatal period. Furthermore, we provide a comprehensive review of previously reported variants in the AGTR1 gene, which is the least encountered genetic cause of RTD, along with their associated clinical features.

The link between glycemic control measures and eye microvascular complications in a clinical cohort of type 2 diabetes with microRNA-223-3p signature.

BACKGROUND: Type 2 diabetes (T2D) is a critical healthcare challenge and priority in Qatar which is listed amongst the top 10 countries in the world, with its prevalence presently at 17% double the global average. MicroRNAs (miRNAs) are implicated in the pathogenesis of (T2D) and long-term microvascular complications including diabetic retinopathy (DR). METHODS: In this study, a T2D cohort that accurately matches the characteristics of the general population was employed to find microRNA (miRNA) signatures that are correlated with glycemic and ß cell function measurements. Targeted miRNA profiling was performed in (471) T2D individuals with or without DR and (491) (non-diabetic) healthy controls from the Qatar Biobank. Discovery analysis identified 20 differentially expressed miRNAs in T2D compared to controls, of which miR-223-3p was significantly upregulated (fold change:5.16, p = 3.6e-02) and positively correlated with glucose and hemoglobin A1c (HbA1c) levels (p-value = 9.88e-04 and 1.64e-05, respectively), but did not show any significant associations with insulin or C-peptide. Accordingly, we performed functional validation using a miR-223-3p mimic (overexpression) under control and hyperglycemia-induced conditions in a zebrafish model. RESULTS: Over-expression of miR-223-3p alone was associated with significantly higher glucose (42.7 mg/dL, n = 75 vs 38.7 mg/dL, n = 75, p = 0.02) and degenerated retinal vasculature, and altered retinal morphology involving changes in the ganglion cell layer and inner and outer nuclear layers. Assessment of retinal angiogenesis revealed significant upregulation in the expression of vascular endothelial growth factor and its receptors, including kinase insert domain receptor. Further, the pancreatic markers, pancreatic and duodenal homeobox 1, and the insulin gene expressions were upregulated in the miR-223-3p group. CONCLUSION: Our zebrafish model validates a novel correlation between miR-223-3p and DR development. Targeting miR-223-3p in T2D patients may serve as a promising therapeutic strategy to control DR in at-risk individuals.

Biallelic variants in SLC38A3 encoding a glutamine transporter cause epileptic encephalopathy.

The solute carrier (SLC) superfamily encompasses >400 transmembrane transporters involved in the exchange of amino acids, nutrients, ions, metals, neurotransmitters and metabolites across biological membranes. SLCs are highly expressed in the mammalian brain; defects in nearly 100 unique SLC-encoding genes (OMIM: https://www.omim.org) are associated with rare Mendelian disorders including developmental and epileptic encephalopathy and severe neurodevelopmental disorders. Exome sequencing and family-based rare variant analyses on a cohort with neurodevelopmental disorders identified two siblings with developmental and epileptic encephalopathy and a shared deleterious homozygous splicing variant in SLC38A3. The gene encodes SNAT3, a sodium-coupled neutral amino acid transporter and a principal transporter of the amino acids asparagine, histidine, and glutamine, the latter being the precursor for the neurotransmitters GABA and glutamate. Additional subjects with a similar developmental and epileptic encephalopathy phenotype and biallelic predicted-damaging SLC38A3 variants were ascertained through GeneMatcher and collaborations with research and clinical molecular diagnostic laboratories. Untargeted metabolomic analysis was performed to identify novel metabolic biomarkers. Ten individuals from seven unrelated families from six different countries with deleterious biallelic variants in SLC38A3 were identified. Global developmental delay, intellectual disability, hypotonia, and absent speech were common features while microcephaly, epilepsy, and visual impairment were present in the majority. Epilepsy was drug-resistant in half. Metabolomic analysis revealed perturbations of glutamate, histidine, and nitrogen metabolism in plasma, urine, and CSF of selected subjects, potentially representing biomarkers of disease. Our data support the contention that SLC38A3 is a novel disease gene for developmental and epileptic encephalopathy and illuminate the likely pathophysiology of the disease as perturbations in glutamine homeostasis.

Assessing the genetic burden of familial hypercholesterolemia in a large middle eastern biobank.

BACKGROUND: The genetic architecture underlying Familial Hypercholesterolemia (FH) in Middle Eastern Arabs is yet to be fully described, and approaches to assess this from population-wide biobanks are important for public health planning and personalized medicine. METHODS: We evaluate the pilot phase cohort (n = 6,140 adults) of the Qatar Biobank (QBB) for FH using the Dutch Lipid Clinic Network (DLCN) criteria, followed by an in-depth characterization of all genetic alleles in known dominant (LDLR, APOB, and PCSK9) and recessive (LDLRAP1, ABCG5, ABCG8, and LIPA) FH-causing genes derived from whole-genome sequencing (WGS). We also investigate the utility of a globally established 12-SNP polygenic risk score to predict FH individuals in this cohort with Arab ancestry. RESULTS: Using DLCN criteria, we identify eight (0.1%) 'definite', 41 (0.7%) 'probable' and 334 (5.4%) 'possible' FH individuals, estimating a prevalence of 'definite or probable' FH in the Qatari cohort of ~ 1:125. We identify ten previously known pathogenic single-nucleotide variants (SNVs) and 14 putatively novel SNVs, as well as one novel copy number variant in PCSK9. Further, despite the modest sample size, we identify one homozygote for a known pathogenic variant (ABCG8, p. Gly574Arg, global MAF = 4.49E-05) associated with Sitosterolemia 2. Finally, calculation of polygenic risk scores found that individuals with 'definite or probable' FH have a significantly higher LDL-C SNP score than 'unlikely' individuals (p = 0.0003), demonstrating its utility in Arab populations. CONCLUSION: We design and implement a standardized approach to phenotyping a population biobank for FH risk followed by systematically identifying known variants and assessing putative novel variants contributing to FH burden in Qatar. Our results motivate similar studies in population-level biobanks - especially those with globally under-represented ancestries - and highlight the importance of genetic screening programs for early detection and management of individuals with high FH risk in health systems.

Patterns and distribution of de novo mutations in multiplex Middle Eastern families.

While de novo mutations (DNMs) are key to genetic diversity, they are also responsible for a high number of rare disorders. To date, no study has systematically examined the rate and distribution of DNMs in multiplex families in highly consanguineous populations. Leveraging WGS profiles of 645 individuals in 146 families, we implemented a combinatorial approach using 3 complementary tools for DNM discovery in 353 unique trio combinations. We found a total of 27,168 DNMs (median: 70 single-nucleotide and 6 insertion-deletions per individual). Phasing revealed around 80% of DNMs were paternal in origin. Notably, using whole-genome methylation data of spermatogonial stem cells, these DNMs were significantly more likely to occur at highly methylated CpGs (OR: 2.03; p value = 6.62 × 10-11). We then examined the effects of consanguinity and ethnicity on DNMs, and found that consanguinity does not seem to correlate with DNM rate, and special attention has to be considered while measuring such a correlation. Additionally, we found that Middle-Eastern families with Arab ancestry had fewer DNMs than African families, although not significant (p value = 0.16). Finally, for families with diseased probands, we examined the difference in DNM counts and putative impact across affected and unaffected siblings, but did not find significant differences between disease groups, likely owing to the enrichment for recessive disorders in this part of the world, or the small sample size per clinical condition. This study serves as a reference for DNM discovery in multiplex families from the globally under-represented populations of the Middle-East.

An ancestral 10-bp repeat expansion in VWA1 causes recessive hereditary motor neuropathy.

The extracellular matrix comprises a network of macromolecules such as collagens, proteoglycans and glycoproteins. VWA1 (von Willebrand factor A domain containing 1) encodes a component of the extracellular matrix that interacts with perlecan/collagen VI, appears to be involved in stabilizing extracellular matrix structures, and demonstrates high expression levels in tibial nerve. Vwa1-deficient mice manifest with abnormal peripheral nerve structure/function; however, VWA1 variants have not previously been associated with human disease. By interrogating the genome sequences of 74 180 individuals from the 100K Genomes Project in combination with international gene-matching efforts and targeted sequencing, we identified 17 individuals from 15 families with an autosomal-recessive, non-length dependent, hereditary motor neuropathy and rare biallelic variants in VWA1. A single disease-associated allele p.(G25Rfs*74), a 10-bp repeat expansion, was observed in 14/15 families and was homozygous in 10/15. Given an allele frequency in European populations approaching 1/1000, the seven unrelated homozygote individuals ascertained from the 100K Genomes Project represents a substantial enrichment above expected. Haplotype analysis identified a shared 220 kb region suggesting that this founder mutation arose >7000 years ago. A wide age-range of patients (6-83 years) helped delineate the clinical phenotype over time. The commonest disease presentation in the cohort was an early-onset (mean 2.0 ± 1.4 years) non-length-dependent axonal hereditary motor neuropathy, confirmed on electrophysiology, which will have to be differentiated from other predominantly or pure motor neuropathies and neuronopathies. Because of slow disease progression, ambulation was largely preserved. Neurophysiology, muscle histopathology, and muscle MRI findings typically revealed clear neurogenic changes with single isolated cases displaying additional myopathic process. We speculate that a few findings of myopathic changes might be secondary to chronic denervation rather than indicating an additional myopathic disease process. Duplex reverse transcription polymerase chain reaction and immunoblotting using patient fibroblasts revealed that the founder allele results in partial nonsense mediated decay and an absence of detectable protein. CRISPR and morpholino vwa1 modelling in zebrafish demonstrated reductions in motor neuron axonal growth, synaptic formation in the skeletal muscles and locomotive behaviour. In summary, we estimate that biallelic variants in VWA1 may be responsible for up to 1% of unexplained hereditary motor neuropathy cases in Europeans. The detailed clinical characterization provided here will facilitate targeted testing on suitable patient cohorts. This novel disease gene may have previously evaded detection because of high GC content, consequential low coverage and computational difficulties associated with robustly detecting repeat-expansions. Reviewing previously unsolved exomes using lower QC filters may generate further diagnoses.

Transcriptome Profile Identifies Actin as an Essential Regulator of Cardiac Myosin Binding Protein C3 Hypertrophic Cardiomyopathy in a Zebrafish Model.

Variants in cardiac myosin-binding protein C (cMyBP-C) are the leading cause of inherited hypertrophic cardiomyopathy (HCM), demonstrating the key role that cMyBP-C plays in the heart's contractile machinery. To investigate the c-MYBPC3 HCM-related cardiac impairment, we generated a zebrafish mypbc3-knockout model. These knockout zebrafish displayed significant morphological heart alterations related to a significant decrease in ventricular and atrial diameters at systolic and diastolic states at the larval stages. Immunofluorescence staining revealed significant hyperplasia in the mutant's total cardiac and ventricular cardiomyocytes. Although cardiac contractility was similar to the wild-type control, the ejection fraction was significantly increased in the mypbc3 mutants. At later stages of larval development, the mutants demonstrated an early cardiac phenotype of myocardium remodeling, concurrent cardiomyocyte hyperplasia, and increased ejection fraction as critical processes in HCM initiation to counteract the increased ventricular myocardial wall stress. The examination of zebrafish adults showed a thickened ventricular cardiac wall with reduced heart rate, swimming speed, and endurance ability in both the mypbc3 heterozygous and homozygous groups. Furthermore, heart transcriptome profiling showed a significant downregulation of the actin-filament-based process, indicating an impaired actin cytoskeleton organization as the main dysregulating factor associated with the early ventricular cardiac hypertrophy in the zebrafish mypbc3 HCM model.

Point of Care Exome Sequencing Reveals Allelic and Phenotypic Heterogeneity Underlying Mendelian disease in Qatar.

The effectiveness of next generation sequencing at solving genetic disease has motivated the rapid adoption of this technology into clinical practice around the world. In this study, we use whole exome sequencing (WES) to assess 48 patients with Mendelian disease from 30 serial families as part of the "Qatar Mendelian Disease pilot program" - a coordinated multi-center effort to build capacity and clinical expertise in genetic medicine in Qatar. By enrolling whole families (parents plus available siblings), we demonstrate significantly improved discriminatory power for candidate variant identification over trios for both de novo and recessive inheritance patterns. For the same index cases, we further demonstrate that even in the absence of families, variant prioritization is improved up to 8-fold when a modest set of population-matched controls is used vs large public databases, stressing the poor representation of Middle Eastern alleles in presently available databases. Our in-house pipeline identified candidate disease variants in 27 of 30 families (90%), 23 of which (85%) harbor novel pathogenic variants in known disease genes, pointing to significant allelic heterogeneity and founder mutations underlying Mendelian disease in the Middle East. For 6 of these families, the clinical presentation was only partially explained by the candidate gene, suggesting phenotypic expansion of known syndromes. Our pilot study demonstrates the utility of WES for Middle Eastern populations, the dramatic improvement in variant prioritization conferred by enrolling population-matched controls and/or enrolling additional unaffected siblings at the point-of-care, and 25 novel disease-causing alleles, relevant to newborn and premarital screening panels in regional populations.

Diagnosis and treatment of type 1 diabetes at the dawn of the personalized medicine era.

Type 1 diabetes affects millions of people globally and requires careful management to avoid serious long-term complications, including heart and kidney disease, stroke, and loss of sight. The type 1 diabetes patient cohort is highly heterogeneous, with individuals presenting with disease at different stages and severities, arising from distinct etiologies, and overlaying varied genetic backgrounds. At present, the "one-size-fits-all" treatment for type 1 diabetes is exogenic insulin substitution therapy, but this approach fails to achieve optimal blood glucose control in many individuals. With advances in our understanding of early-stage diabetes development, diabetes stratification, and the role of genetics, type 1 diabetes is a promising candidate for a personalized medicine approach, which aims to apply "the right therapy at the right time, to the right patient". In the case of type 1 diabetes, great efforts are now being focused on risk stratification for diabetes development to enable pre-clinical detection, and the application of treatments such as gene therapy, to prevent pancreatic destruction in a sub-set of patients. Alongside this, breakthroughs in stem cell therapies hold great promise for the regeneration of pancreatic tissues in some individuals. Here we review the recent initiatives in the field of personalized medicine for type 1 diabetes, including the latest discoveries in stem cell and gene therapy for the disease, and current obstacles that must be overcome before the dream of personalized medicine for all type 1 diabetes patients can be realized.

Functional characterization of human myosin-binding protein C3 variants associated with hypertrophic cardiomyopathy reveals exon-specific cardiac phenotypes in zebrafish model.

Myosin-binding protein C 3 (MYBPC3) variants are the most common cause of hypertrophic cardiomyopathy (HCM). HCM is a complex cardiac disorder due to its significant genetic and clinical heterogeneity. MYBPC3 variants genotype-phenotype associations remain poorly understood. We investigated the impact of two novel human MYBPC3 splice-site variants: V1: c.654+2_654+4dupTGG targeting exon 5 using morpholino MOe5i5; and V2: c.772+1G>A targeting exon 6 using MOe6i6; located within C1 domain of cMyBP-C protein, known to be critical in regulating sarcomere structure and contractility. Zebrafish MOe5i5 and MOe6i6 morphants recapitulated typical characteristics of human HCM with cardiac phenotypes of varying severity, including reduced cardiomyocyte count, thickened ventricular myocardial wall, a drastic reduction in heart rate, stroke volume, and cardiac output. Analysis of all cardiac morphological and functional parameters demonstrated that V2 cardiac phenotype was more severe than V1. Coinjection with synthetic human MYBPC3 messenger RNA (mRNA) partially rescued disparate cardiac phenotypes in each zebrafish morphant. While human MYBPC3 mRNA partially restored the decreased heart rate in V1 morphants and displayed increased percentages of ejection fraction, fractional shortening, and area change, it failed to revert the V1 ventricular myocardial thickness. These results suggest a possible V1 impact on cardiac contractility. In contrast, attempts to rescue V2 morphants only restored the ventricular myocardial wall hypertrophy phenotype but had no significant effect on impaired heart rate, suggesting a potential V2 impact on the cardiac structure. Our study provides evidence of an association between MYBPC3 exon-specific cardiac phenotypes in the zebrafish model providing important insights into how these genetic variants contribute to HCM disease.

Willingness to participate in genome testing: a survey of public attitudes from Qatar.

Genomics has the potential to revolutionize medical approaches to disease prevention, diagnosis, and treatment, but it does not come without challenges. The success of a national population-based genome program, like the Qatar Genome Program (QGP), depends on the willingness of citizens to donate samples and take up genomic testing services. This study explores public attitudes of the Qatari population toward genetic testing and toward participating in the QGP. A representative sample of 837 adult Qataris was surveyed in May 2016. Approximately 71% of respondents surveyed reported that they were willing to participate in the activities of the QGP. Willingness to participate was significantly associated with basic literacy in genetics, a family history of genetic diseases, and previous experience with genetic testing through premarital screening. Respondents cited the desire to know more about their health status as the principle motivation for participating, while lack of time and information were reported as the most important barriers. With QGP plans to ramp up the scale of its national operation toward more integration into clinical care settings, it is critical to understand public attitudes and their determinants. The results demonstrate public support but also identify the need for more education and individual counseling that not only provide information on the process, challenges, and benefits of genomic testing, but that also address concerns about information security.

De novo mutations in histone-modifying genes in congenital heart disease.

Congenital heart disease (CHD) is the most frequent birth defect, affecting 0.8% of live births. Many cases occur sporadically and impair reproductive fitness, suggesting a role for de novo mutations. Here we compare the incidence of de novo mutations in 362 severe CHD cases and 264 controls by analysing exome sequencing of parent-offspring trios. CHD cases show a significant excess of protein-altering de novo mutations in genes expressed in the developing heart, with an odds ratio of 7.5 for damaging (premature termination, frameshift, splice site) mutations. Similar odds ratios are seen across the main classes of severe CHD. We find a marked excess of de novo mutations in genes involved in the production, removal or reading of histone 3 lysine 4 (H3K4) methylation, or ubiquitination of H2BK120, which is required for H3K4 methylation. There are also two de novo mutations in SMAD2, which regulates H3K27 methylation in the embryonic left-right organizer. The combination of both activating (H3K4 methylation) and inactivating (H3K27 methylation) chromatin marks characterizes 'poised' promoters and enhancers, which regulate expression of key developmental genes. These findings implicate de novo point mutations in several hundreds of genes that collectively contribute to approximately 10% of severe CHD.

Point-of-care whole-exome sequencing of idiopathic male infertility.

PURPOSE: Nonobstructive azoospermia (NOA) affects 1% of the male population; however, despite state-of-the-art clinical assessment, for most patients the cause is unknown. We capitalized on an analysis of multiplex families in the Middle East to identify highly penetrant genetic causes. METHODS: We used whole-exome sequencing (WES) in 8 consanguineous families and combined newly discovered genes with previously reported ones to create a NOA gene panel, which was used to identify additional variants in 75 unrelated idiopathic NOA subjects and 74 fertile controls. RESULTS: In five of eight families, we identified rare deleterious recessive variants in CCDC155, NANOS2, SPO11, TEX14, and WNK3 segregating with disease. These genes, which are novel to human NOA, have remarkable testis-specific expression, and murine functional evidence supports roles for them in spermatogenesis. Among 75 unrelated NOA subjects, we identified 4 (~5.3%) with additional recessive variants in these newly discovered genes and 6 with deleterious variants in previously reported NOA genes, yielding an overall genetic etiology for 13.3% subjects versus 0 fertile controls (p = 0.001). CONCLUSION: NOA affects millions of men, many of whom remain idiopathic despite extensive laboratory evaluation. The genetic etiology for a substantial fraction of these patients (>50% familial and >10% sporadic) may be discovered by WES at the point of care.

Two hits in one: whole genome sequencing unveils LIG4 syndrome and urofacial syndrome in a case report of a child with complex phenotype.

BACKGROUND: Ligase IV syndrome, a hereditary disease associated with compromised DNA damage response mechanisms, and Urofacial syndrome, caused by an impairment of neural cell signaling, are both rare genetic disorders, whose reports in literature are limited. We describe the first case combining both disorders in a specific phenotype. CASE PRESENTATION: We report a case of a 7-year old girl presenting with a complex phenotype characterized by multiple congenital abnormalities and dysmorphic features, microcephaly, short stature, combined immunodeficiency and severe vesicoureteral reflux. Whole Genome Sequencing was performed and a novel ligase IV homozygous missense c.T1312C/p.Y438H mutation was detected, and is believed to be responsible for most of the clinical features of the child, except vesicoureteral reflux which has not been previously described for ligase IV deficiency. However, we observed a second rare damaging (nonsense) homozygous mutation (c.C2125T/p.R709X) in the leucine-rich repeats and immunoglobulin-like domains 2 gene that encodes a protein implicated in neural cell signaling and oncogenesis. Interestingly, this mutation has recently been reported as pathogenic and causing urofacial syndrome, typically displaying vesicoureteral reflux. Thus, this second mutation completes the missing genetic explanation for this intriguing clinical puzzle. We verified that both mutations fit an autosomal recessive inheritance model due to extensive consanguinity. CONCLUSIONS: We successfully identified a novel ligase IV mutation, causing ligase IV syndrome, and an additional rare leucine-rich repeats and immunoglobulin-like domains 2 gene nonsense mutation, in the context of multiple autosomal recessive conditions due to extensive consanguinity. This work demonstrates the utility of Whole Genome Sequencing data in clinical diagnosis in such cases where the combination of multiple rare phenotypes results in very intricate clinical pictures. It also reports a novel causative mutation and a clinical phenotype, which will help in better defining the essential features of both ligase IV and leucine-rich repeats and immunoglobulin-like domains 2 deficiency syndromes.

Copy number variations in the genome of the Qatari population.

BACKGROUND: The populations of the Arabian Peninsula remain the least represented in public genetic databases, both in terms of single nucleotide variants and of larger genomic mutations. We present the first high-resolution copy number variation (CNV) map for a Gulf Arab population, using a hybrid approach that integrates array genotyping intensity data and next-generation sequencing reads to call CNVs in the Qatari population. METHODS: CNVs were detected in 97 unrelated Qatari individuals by running two calling algorithms on each of two primary datasets: high-resolution genotyping (Illumina Omni 2.5M) and high depth whole-genome sequencing (Illumina PE 100bp). The four call-sets were integrated to identify high confidence CNV regions, which were subsequently annotated for putative functional effect and compared to public databases of CNVs in other populations. The availability of genome sequence was leveraged to identify tagging SNPs in high LD with common deletions in this population, enabling their imputation from genotyping experiments in the future. RESULTS: Genotyping intensities and genome sequencing data from 97 Qataris were analyzed with four different algorithms and integrated to discover 16,660 high confidence CNV regions (CNVRs) in the total population, affecting ~28 Mb in the median Qatari genome. Up to 40% of all CNVs affected genes, including novel CNVs affecting Mendelian disease genes, segregating at different frequencies in the 3 major Qatari subpopulations, including those with Bedouin, Persian/South Asian, and African ancestry. Consistent with high consanguinity levels in the Bedouin subpopulation, we found an increased burden for homozygous deletions in this group. In comparison to known CNVs in the comprehensive Database of Genomic Variants, we found that 5% of all CNVRs in Qataris were completely novel, with an enrichment of CNVs affecting several known chromosomal disorder loci and genes known to regulate sugar metabolism and type 2 diabetes in the Qatari cohort. Finally, we leveraged the availability of genome sequence to find suitable tagging SNPs for common deletions in this population. CONCLUSION: We combine four independently generated datasets from 97 individuals to study CNVs for the first time at high-resolution in a Gulf Arab population.

Rare copy number variations in congenital heart disease patients identify unique genes in left-right patterning.

Dominant human genetic diseases that impair reproductive fitness and have high locus heterogeneity constitute a problem for gene discovery because the usual criterion of finding more mutations in specific genes than expected by chance may require extremely large populations. Heterotaxy (Htx), a congenital heart disease resulting from abnormalities in left-right (LR) body patterning, has features suggesting that many cases fall into this category. In this setting, appropriate model systems may provide a means to support implication of specific genes. By high-resolution genotyping of 262 Htx subjects and 991 controls, we identify a twofold excess of subjects with rare genic copy number variations in Htx (14.5% vs. 7.4%, P = 1.5 × 10(-4)). Although 7 of 45 Htx copy number variations were large chromosomal abnormalities, 38 smaller copy number variations altered a total of 61 genes, 22 of which had Xenopus orthologs. In situ hybridization identified 7 of these 22 genes with expression in the ciliated LR organizer (gastrocoel roof plate), a marked enrichment compared with 40 of 845 previously studied genes (sevenfold enrichment, P < 10(-6)). Morpholino knockdown in Xenopus of Htx candidates demonstrated that five (NEK2, ROCK2, TGFBR2, GALNT11, and NUP188) strongly disrupted both morphological LR development and expression of pitx2, a molecular marker of LR patterning. These effects were specific, because 0 of 13 control genes from rare Htx or control copy number variations produced significant LR abnormalities (P = 0.001). These findings identify genes not previously implicated in LR patterning.