Search

Many neurodevelopmental disorders (NDDs) share common learning and behavioural impairments, as well as features such as dysregulation of the oxytocin hormone. Here, we examined DNA methylation (DNAm) in the 1st intron of the oxytocin receptor gene, "OXTR," in patients with autism spectrum (ASD), attention deficit and hyperactivity (ADHD) and obsessive compulsive (OCD) disorders. DNAm of "OXTR" was assessed for cohorts of ASD (blood), ADHD (saliva), OCD (saliva), which uncovered sex-specific DNAm differences compared to neurotypical, tissue-matched controls. Individuals with ASD or ADHD exhibiting extreme DNAm values had lower IQ and more social problems, respectively, than those with DNAm within normative ranges. This suggests that "OXTR" DNAm patterns are altered across NDDs and may be correlated with common clinical outcomes.

Anaplastic thyroid carcinoma is arguably the most lethal human malignancy. It often co-occurs with differentiated thyroid cancers, yet the molecular origins of its aggressivity are unknown. We sequenced tumor DNA from 329 regions of thyroid cancer, including 213 from patients with primary anaplastic thyroid carcinomas. We also whole genome sequenced 9 patients using multi-region sequencing of both differentiated and anaplastic thyroid cancer components. Using these data, we demonstrate thatanaplastic thyroid carcinomas have a higher burden of mutations than other thyroid cancers, with distinct mutational signatures and molecular subtypes. Further, different cancer driver genes are mutated in anaplastic and differentiated thyroid carcinomas, even those arising in a single patient. Finally, we unambiguously demonstrate that anaplastic thyroid carcinomas share a genomic origin with co-occurring differentiated carcinomas and emerge from a common malignant field through acquisition of characteristic clonal driver mutations.

OBJECTIVE: Suicidal behavior is heritable and is a major cause of death worldwide. Two large-scale genome-wide association studies (GWASs) recently discovered and cross-validated genome-wide significant (GWS) loci for suicide attempt (SA). The present study leveraged the genetic cohorts from both studies to conduct the largest GWAS meta-analysis of SA to date. Multi-ancestry and admixture-specific meta-analyses were conducted within groups of significant African, East Asian, and European ancestry admixtures. METHODS: This study comprised 22 cohorts, including 43,871 SA cases and 915,025 ancestry-matched controls. Analytical methods across multi-ancestry and individual ancestry admixtures included inverse variance-weighted fixed-effects meta-analyses, followed by gene, gene-set, tissue-set, and drug-target enrichment, as well as summary-data-based Mendelian randomization with brain expression quantitative trait loci data, phenome-wide genetic correlation, and genetic causal proportion analyses. RESULTS: Multi-ancestry and European ancestry admixture GWAS meta-analyses identified 12 risk loci at p values

Cytogenetic abnormalities in the Prader-Willi/Angelman syndrome (PWS/AS) critical region have been described in individuals with autism. Maternal duplications and linkage disequilibrium in families with autism suggest the existence of a susceptibility locus at 15q11-q13. Here, we describe a 6-year-old girl diagnosed with autism, developmental delay, and delayed expressive and receptive language. The karyotype was designated de novo 47, XX, idic(15)(q13). Fluorescence in situ hybridization (FISH) and molecular analysis with 15q11-q13 markers revealed an additional copy of the region being of maternal origin. Duplication of the 15q11-q13 segment represents the most consistent known chromosomal abnormality reported in association with autism. This present case report reinforces the hypothesis that additional copies of this chromosome segment are causally related to autism.

Posttraumatic stress disorder (PTSD) is a heritable (h2 = 24-71%) psychiatric illness. Copy number variation (CNV) is a form of rare genetic variation that has been implicated in the etiology of psychiatric disorders, but no large-scale investigation of CNV in PTSD has been performed. We present an association study of CNV burden and PTSD symptoms in a sample of 114,383 participants (13,036 cases and 101,347 controls) of European ancestry. CNVs were called using two calling algorithms and intersected to a consensus set. Quality control was performed to remove strong outlier samples. CNVs were examined for association with PTSD within each cohort using linear or logistic regression analysis adjusted for population structure and CNV quality metrics, then inverse variance weighted meta-analyzed across cohorts. We examined the genome-wide total span of CNVs, enrichment of CNVs within specified gene-sets, and CNVs overlapping individual genes and implicated neurodevelopmental regions. The total distance covered by deletions crossing over known neurodevelopmental CNV regions was significant (beta = 0.029, SE = 0.005, P = 6.3 × 10-8). The genome-wide neurodevelopmental CNV burden identified explains 0.034% of the variation in PTSD symptoms. The 15q11.2 BP1-BP2 microdeletion region was significantly associated with PTSD (beta = 0.0206, SE = 0.0056, P = 0.0002). No individual significant genes interrupted by CNV were identified. 22 gene pathways related to the function of the nervous system and brain were significant in pathway analysis (FDR q

The additional sex combs-like (ASXL) gene family-encoded by ASXL1, ASXL2, and ASXL3-is crucial for mammalian development. Pathogenic variants in the ASXL gene family are associated with three phenotypically distinct neurodevelopmental syndromes. Our previous work has shown that syndromic conditions caused by pathogenic variants in epigenetic regulatory genes show consistent patterns of genome-wide DNA methylation (DNAm) alterations, i.e., DNAm signatures in peripheral blood. Given the role of ASXL1 in chromatin modification, we hypothesized that pathogenic ASXL1 variants underlying Bohring-Opitz syndrome (BOS) have a unique DNAm signature. We profiled whole-blood DNAm for 17 ASXL1 variants, and 35 sex- and age-matched typically developing individuals, using Illumina's Infinium EPIC array. We identified 763 differentially methylated CpG sites in individuals with BOS. Differentially methylated sites overlapped 323 unique genes, including HOXA5 and HOXB4, supporting the functional relevance of DNAm signatures. We used a machine-learning classification model based on the BOS DNAm signature to classify variants of uncertain significance in ASXL1, as well as pathogenic ASXL2 and ASXL3 variants. The DNAm profile of one individual with the ASXL2 variant was BOS-like, whereas the DNAm profiles of three individuals with ASXL3 variants were control-like. We also used Horvath's epigenetic clock, which showed acceleration in DNAm age in individuals with pathogenic ASXL1 variants, and the individual with the pathogenic ASXL2 variant, but not in individuals with ASXL3 variants. These studies enhance our understanding of the epigenetic dysregulation underpinning ASXL gene family-associated syndromes.

Rare genomic disorders (RGDs) confer elevated risk for neurodevelopmental psychiatric disorders. In this era of intense genomics discoveries, the landscape of RGDs is rapidly evolving. However, there has not been comparable progress to date in scalable, harmonized phenotyping methods. As a result, beyond associations with categorical diagnoses, the effects on dimensional traits remain unclear for many RGDs. The nature and specificity of RGD effects on cognitive and behavioral traits is an area of intense investigation: RGDs are frequently associated with more than one psychiatric condition, and those studied to date affect, to varying degrees, a broad range of developmental and cognitive functions. Although many RGDs have large effects, phenotypic expression is typically influenced by additional genomic and environmental factors. There is emerging evidence that using polygenic risk scores in individuals with RGDs offers opportunities to refine prediction, thus allowing for the identification of those at greatest risk of psychiatric illness. However, translation into the clinic is hindered by roadblocks, which include limited genetic testing in clinical psychiatry, and the lack of guidelines for following individuals with RGDs, who are at high risk of developing psychiatric symptoms. The Genes to Mental Health Network (G2MH) is a newly funded National Institute of Mental Health initiative that will collect, share, and analyze large-scale data sets combining genomics and dimensional measures of psychopathology spanning diverse populations and geography. The authors present here the most recent understanding of the effects of RGDs on dimensional behavioral traits and risk for psychiatric conditions and discuss strategies that will be pursued within the G2MH network, as well as how expected results can be translated into clinical practice to improve patient outcomes.

BackgroundSuicide is a leading cause of death worldwide, and nonfatal suicide attempts, which occur far more frequently, are a major source of disability and social and economic burden. Both have substantial genetic etiology, which is partially shared and partially distinct from that of related psychiatric disorders.MethodsWe conducted a genome-wide association study (GWAS) of 29,782 suicide attempt (SA) cases and 519,961 controls in the International Suicide Genetics Consortium (ISGC). The GWAS of SA was conditioned on psychiatric disorders using GWAS summary statistics via multitrait-based conditional and joint analysis, to remove genetic effects on SA mediated by psychiatric disorders. We investigated the shared and divergent genetic architectures of SA, psychiatric disorders, and other known risk factors.ResultsTwo loci reached genome-wide significance for SA: the major histocompatibility complex and an intergenic locus on chromosome 7, the latter of which remained associated with SA after conditioning on psychiatric disorders and replicated in an independent cohort from the Million Veteran Program. This locus has been implicated in risk-taking behavior, smoking, and insomnia. SA showed strong genetic correlation with psychiatric disorders, particularly major depression, and also with smoking, pain, risk-taking behavior, sleep disturbances, lower educational attainment, reproductive traits, lower socioeconomic status, and poorer general health. After conditioning on psychiatric disorders, the genetic correlations between SA and psychiatric disorders decreased, whereas those with nonpsychiatric traits remained largely unchanged.ConclusionsOur results identify a risk locus that contributes more strongly to SA than other phenotypes and suggest a shared underlying biology between SA and known risk factors that is not mediated by psychiatric disorders.

The particularly interdisciplinary nature of human microbiome research makes the organization and reporting of results spanning epidemiology, biology, bioinformatics, translational medicine and statistics a challenge. Commonly used reporting guidelines for observational or genetic epidemiology studies lack key features specific to microbiome studies. Therefore, a multidisciplinary group of microbiome epidemiology researchers adapted guidelines for observational and genetic studies to culture-independent human microbiome studies, and also developed new reporting elements for laboratory, bioinformatics and statistical analyses tailored to microbiome studies. The resulting tool, called 'Strengthening The Organization and Reporting of Microbiome Studies' (STORMS), is composed of a 17-item checklist organized into six sections that correspond to the typical sections of a scientific publication, presented as an editable table for inclusion in supplementary materials. The STORMS checklist provides guidance for concise and complete reporting of microbiome studies that will facilitate manuscript preparation, peer review, and reader comprehension of publications and comparative analysis of published results.

Truncating variants in exons 33 and 34 of the SNF2-related CREBBP activator protein (SRCAP) gene cause the neurodevelopmental disorder (NDD) Floating-Harbor syndrome (FLHS), characterized by short stature, speech delay, and facial dysmorphism. Here, we present a cohort of 33 individuals with clinical features distinct from FLHS and truncating (mostly de novo) SRCAP variants either proximal (n = 28) or distal (n = 5) to the FLHS locus. Detailed clinical characterization of the proximal SRCAP individuals identified shared characteristics: developmental delay with or without intellectual disability, behavioral and psychiatric problems, non-specific facial features, musculoskeletal issues, and hypotonia. Because FLHS is known to be associated with a unique set of DNA methylation (DNAm) changes in blood, a DNAm signature, we investigated whether there was a distinct signature associated with our affected individuals. A machine-learning model, based on the FLHS DNAm signature, negatively classified all our tested subjects. Comparing proximal variants with typically developing controls, we identified a DNAm signature distinct from the FLHS signature. Based on the DNAm and clinical data, we refer to the condition as "non-FLHS SRCAP-related NDD." All five distal variants classified negatively using the FLHS DNAm model while two classified positively using the proximal model. This suggests divergent pathogenicity of these variants, though clinically the distal group presented with NDD, similar to the proximal SRCAP group. In summary, for SRCAP, there is a clear relationship between variant location, DNAm profile, and clinical phenotype. These results highlight the power of combined epigenetic, molecular, and clinical studies to identify and characterize genotype-epigenotype-phenotype correlations.

The 22q11.2 deletion syndrome (22q11DS) is associated with a 20-25% risk of schizophrenia. In a cohort of 962 individuals with 22q11DS, we examined the shared genetic basis between schizophrenia and schizophrenia-related early trajectory phenotypes: sub-threshold symptoms of psychosis, low baseline intellectual functioning and cognitive decline. We studied the association of these phenotypes with two polygenic scores, derived for schizophrenia and intelligence, and evaluated their use for individual risk prediction in 22q11DS. Polygenic scores were not only associated with schizophrenia and baseline intelligence quotient (IQ), respectively, but schizophrenia polygenic score was also significantly associated with cognitive (verbal IQ) decline and nominally associated with sub-threshold psychosis. Furthermore, in comparing the tail-end deciles of the schizophrenia and IQ polygenic score distributions, 33% versus 9% of individuals with 22q11DS had schizophrenia, and 63% versus 24% of individuals had intellectual disability. Collectively, these data show a shared genetic basis for schizophrenia and schizophrenia-related phenotypes and also highlight the future potential of polygenic scores for risk stratification among individuals with highly, but incompletely, penetrant genetic variants.

Weaver syndrome (WS), an overgrowth/intellectual disability syndrome (OGID), is caused by pathogenic variants in the histone methyltransferase EZH2, which encodes a core component of the Polycomb repressive complex-2 (PRC2). Using genome-wide DNA methylation (DNAm) data for 187 individuals with OGID and 969 control subjects, we show that pathogenic variants in EZH2 generate a highly specific and sensitive DNAm signature reflecting the phenotype of WS. This signature can be used to distinguish loss-of-function from gain-of-function missense variants and to detect somatic mosaicism. We also show that the signature can accurately classify sequence variants in EED and SUZ12, which encode two other core components of PRC2, and predict the presence of pathogenic variants in undiagnosed individuals with OGID. The discovery of a functionally relevant signature with utility for diagnostic classification of sequence variants in EZH2, EED, and SUZ12 supports the emerging paradigm shift for implementation of DNAm signatures into diagnostics and translational research.

The IPDGC (The International Parkinson Disease Genomics Consortium) and EADB (Alzheimer Disease European DNA biobank) are listed correctly as an author to the article, however, they were incorrectly listed more than once.

Hydrodynamic tail vein injection (HTV) is the “gold standard” for delivering naked DNA vectors to mouse liver, thereby transfecting predominately perivenous hepatocytes. While HTV corrects metabolic liver defects such as phenylketonuria or cystathionine β-synthase deficiency, correction of spfash mice with ornithine transcarbamylase (OTC) deficiency was not possible despite overexpression in the liver, as the OTC enzyme is primarily expressed in periportal hepatocytes. To target periportal hepatocytes, we established hydrodynamic retrograde intrabiliary injection (HRII) in mice and optimized minicircle (MC) vector delivery using luciferase as a marker gene. HRII resulted in a transfection efficiency below 1%, 100-fold lower than HTV. While HRII induced minimal liver toxicity compared with HTV, overexpression of luciferase by both methods, but not of a natural liver-specific enzyme, elicited an immune response that led to the elimination of luciferase expression. Further testing of MC vectors delivered via HRII in spfash mice did not result in sufficient therapeutic efficacy and needs further optimization and/or selection of the corrected cells. This study reveals that luciferase expression is toxic for the liver. Furthermore, physical delivery of MC vectors via the bile duct has the potential to treat defects restricted to periportal hepatocytes, which opens new doors for non-viral liver-directed gene therapy.

Abstract Single-cell DNA methylation profiling currently suffers from excessive noise and/or limited cellular throughput. We developed scTAM-seq, a targeted bisulfite-free method for profiling up to 650 CpGs in up to 10,000 cells per experiment, with a dropout rate as low as 7%. We demonstrate that scTAM-seq can resolve DNA methylation dynamics across B-cell differentiation in blood and bone marrow, identifying intermediate differentiation states that were previously masked. scTAM-seq additionally queries surface-protein expression, thus enabling integration of single-cell DNA methylation information with cell atlas data. In summary, scTAM-seq is a high-throughput, high-confidence method for analyzing DNA methylation at single-CpG resolution across thousands of single cells.

The underpinnings of mild to moderate neurodevelopmental delay remain elusive, often leading to late diagnosis and interventions. Here, we present data on exome and genome sequencing as well as array analysis of 13 individuals that point to pathogenic, heterozygous, mostly de novo variants in WDFY3 (significant de novo enrichment P = 0.003) as a monogenic cause of mild and non-specific neurodevelopmental delay. Nine variants were protein-truncating and four missense. Overlapping symptoms included neurodevelopmental delay, intellectual disability, macrocephaly, and psychiatric disorders (autism spectrum disorders/attention deficit hyperactivity disorder). One proband presented with an opposing phenotype of microcephaly and the only missense-variant located in the PH-domain of WDFY3. Findings of this case are supported by previously published data, demonstrating that pathogenic PH-domain variants can lead to microcephaly via canonical Wnt-pathway upregulation. In a separate study, we reported that the autophagy scaffolding protein WDFY3 is required for cerebral cortical size regulation in mice, by controlling proper division of neural progenitors. Here, we show that proliferating cortical neural progenitors of human embryonic brains highly express WDFY3, further supporting a role for this molecule in the regulation of prenatal neurogenesis. We present data on Wnt-pathway dysregulation in Wdfy3-haploinsufficient mice, which display macrocephaly and deficits in motor coordination and associative learning, recapitulating the human phenotype. Consequently, we propose that in humans WDFY3 loss-of-function variants lead to macrocephaly via downregulation of the Wnt pathway. In summary, we present WDFY3 as a novel gene linked to mild to moderate neurodevelopmental delay and intellectual disability and conclude that variants putatively causing haploinsufficiency lead to macrocephaly, while an opposing pathomechanism due to variants in the PH-domain of WDFY3 leads to microcephaly.

The genetic variant rs72824905-G (minor allele) in the PLCG2 gene was previously associated with a reduced Alzheimer's disease risk (AD). The role of PLCG2 in immune system signaling suggests it may also protect against other neurodegenerative diseases and possibly associates with longevity. We studied the effect of the rs72824905-G on seven neurodegenerative diseases and longevity, using 53,627 patients, 3,516 long-lived individuals and 149,290 study-matched controls. We replicated the association of rs72824905-G with reduced AD risk and we found an association with reduced risk of dementia with Lewy bodies (DLB) and frontotemporal dementia (FTD). We did not find evidence for an effect on Parkinson's disease (PD), amyotrophic lateral sclerosis (ALS) and multiple sclerosis (MS) risks, despite adequate sample sizes. Conversely, the rs72824905-G allele was associated with increased likelihood of longevity. By-proxy analyses in the UK Biobank supported the associations with both dementia and longevity. Concluding, rs72824905-G has a protective effect against multiple neurodegenerative diseases indicating shared aspects of disease etiology. Our findings merit studying the PLCγ2 pathway as drug-target.

High red meat intake is associated with the risk of colorectal cancer (CRC), whereas dietary fibers, such as resistant starch (RS) seemed to protect against CRC. The aim of this study was to determine whether high-amylose potato starch (HAPS), high-amylose maize starch (HAMS), and butyrylated high-amylose maize starch (HAMSB)-produced by an organocatalytic route-could oppose the negative effects of a high-protein meat diet (HPM), in terms of fermentation pattern, cecal microbial composition, and colonic biomarkers of CRC. Rats were fed a HPM diet or an HPM diet where 10% of the maize starch was substituted with either HAPS, HAMS, or HAMSB, for 4 weeks. Feces, cecum digesta, and colonic tissue were obtained for biochemical, microbial, gene expression (oncogenic microRNA), and immuno-histochemical (O6-methyl-2-deoxyguanosine (O6MeG) adduct) analysis. The HAMS and HAMSB diets shifted the fecal fermentation pattern from protein towards carbohydrate metabolism. The HAMSB diet also substantially increased fecal butyrate concentration and the pool, compared with the other diets. All three RS treatments altered the cecal microbial composition in a diet specific manner. HAPS and HAMSB showed CRC preventive effects, based on the reduced colonic oncogenic miR17-92 cluster miRNA expression, but there was no significant diet-induced differences in the colonic O6MeG adduct levels. Overall, HAMSB consumption showed the most potential for limiting the negative effects of a high-meat diet.

Identification of genetic biomarkers associated with autism spectrum disorders (ASDs) could improve recurrence prediction for families with a child with ASD. Here, we describe clinical microarray findings for 253 longitudinally phenotyped ASD families from the Baby Siblings Research Consortium (BSRC), encompassing 288 infant siblings. By age 3, 103 siblings (35.8%) were diagnosed with ASD and 54 (18.8%) were developing atypically. Thirteen siblings have copy number variants (CNVs) involving ASD-relevant genes: 6 with ASD, 5 atypically developing, and 2 typically developing. Within these families, an ASD-related CNV in a sibling has a positive predictive value (PPV) for ASD or atypical development of 0.83; the Simons Simplex Collection of ASD families shows similar PPVs. Polygenic risk analyses suggest that common genetic variants may also contribute to ASD. CNV findings would have been pre-symptomatically predictive of ASD or atypical development in 11 (7%) of the 157 BSRC siblings who were eventually diagnosed clinically.

Despite promising findings, quantitative PCR (qPCR)-based tests for RNA quantification have experienced serious limitations in their clinical application. The noticeable lack of technical standardization remains a huge obstacle in the translation of qPCR-based tests. The incorporation of qPCR-based tests into the clinic will benefit from guidelines for clinical research assay validation. This will ultimately impact the clinical management of the patient, including diagnosis, prognosis, prediction, monitoring of the therapeutic response, and evaluation of toxicity. However, clear assay validation protocols for biomarker investigation in clinical trials using molecular assays are currently lacking. Here, we will focus on the necessary steps, including sample acquisition, processing and storage, RNA purification, target selection, assay design, and experimental design, that need to be taken toward the appropriate validation of qRT-PCR assays in clinical research. These recommendations can fill the gap between research use only (RUO) and in vitro diagnostics (IVD). Our contribution provides a tool for basic and clinical research for the development of validated assays in the intermediate steps of biomarker research. These guidelines are based on the current understanding and consensus within the EU-CardioRNA COST Action consortium (www.cardiorna.eu). Their applicability encompasses all clinical areas.

Hepatic gene therapy by delivering non-integrating therapeutic vectors in newborns remains challenging due to the risk of dilution and loss of efficacy in the growing liver. Previously we reported on hepatocyte transfection in piglets by intraportal injection of naked DNA vectors. Here, we established delivery of naked DNA vectors to target periportal hepatocytes in weaned pigs by hydrodynamic retrograde intrabiliary injection (HRII). The surgical procedure involved laparotomy and transient isolation of the liver. For vector delivery, a catheter was placed within the common bile duct by enterotomy. Under optimal conditions, no histological abnormalities were observed in liver tissue upon pressurized injections. The transfection of hepatocytes in all tested liver samples was observed with vectors expressing luciferase from a liver-specific promoter. However, vector copy number and luciferase expression were low compared to hydrodynamic intraportal injection. A 10-fold higher number of vector genomes and luciferase expression was observed in pigs using a non-integrating naked DNA vector with the potential for replication. In summary, the HRII application was less efficient (i.e., lower luciferase activity and vector copy numbers) than the intraportal delivery method but was significantly less distressful for the piglets and has the potential for injection (or re-injection) of vector DNA by endoscopic retrograde cholangiopancreatography.

Abstract Cardiomyopathy (CMP) is a heritable disorder. Over 50% of cases are gene-elusive on clinical gene panel testing. The contribution of variants in non-coding DNA elements that result in cryptic splicing and regulate gene expression has not been explored. We analyzed whole-genome sequencing (WGS) data in a discovery cohort of 209 pediatric CMP patients and 1953 independent replication genomes and exomes. We searched for protein-coding variants, and non-coding variants predicted to affect the function or expression of genes. Thirty-nine percent of cases harbored pathogenic coding variants in known CMP genes, and 5% harbored high-risk loss-of-function (LoF) variants in additional candidate CMP genes. Fifteen percent harbored high-risk regulatory variants in promoters and enhancers of CMP genes (odds ratio 2.25, p = 6.70 × 10−7 versus controls). Genes involved in α-dystroglycan glycosylation (FKTN, DTNA) and desmosomal signaling (DSC2, DSG2) were most highly enriched for regulatory variants (odds ratio 6.7–58.1). Functional effects were confirmed in patient myocardium and reporter assays in human cardiomyocytes, and in zebrafish CRISPR knockouts. We provide strong evidence for the genomic contribution of functionally active variants in new genes and in regulatory elements of known CMP genes to early onset CMP.

Abstract Background Variant interpretation is the main bottleneck in medical genomic sequencing efforts. This usually involves genome analysts manually searching through a multitude of independent databases, often with the aid of several, mostly independent, computational tools. To streamline variant interpretation, we developed the GeneTerpret platform which collates data from current interpretation tools and databases, and applies a phenotype-driven query to categorize the variants identified in the genome(s). The platform assigns quantitative validity scores to genes by query and assembly of the genotype–phenotype data, sequence homology, molecular interactions, expression data, and animal models. It also uses the American College of Medical Genetics and Genomics (ACMG) criteria to categorize variants into five tiers of pathogenicity. The final output is a prioritized list of potentially causal variants/genes. Results We tested GeneTerpret by comparing its performance to expert-curated genes (ClinGen’s gene-validity database) and variant pathogenicity reports (DECIPHER database). Output from GeneTerpret was 97.2% and 83.5% concordant with the expert-curated sources, respectively. Additionally, similar concordance was observed when GeneTerpret’s performance was compared with our internal expert-interpreted clinical datasets. Conclusions GeneTerpret is a flexible platform designed to streamline the genome interpretation process, through a unique interface, with improved ease, speed and accuracy. This modular and customizable system allows the user to tailor the component-programs in the analysis process to their preference. GeneTerpret is available online at https://geneterpret.com .

Abstract Autism spectrum disorder (ASD) is a collection of neurodevelopmental disorders characterized by deficits in social communication and restricted, repetitive patterns of behavior or interests. ASD is highly heritable, but genetically and phenotypically heterogeneous, reducing the power to identify causative genes. We performed whole genome sequencing (WGS) in an ASD cohort of 68 individuals from 22 families enriched for recent shared ancestry. We identified an average of 3.07 million variants per genome, of which an average of 112,512 were rare. We mapped runs of homozygosity (ROHs) in affected individuals and found an average genomic homozygosity of 9.65%, consistent with expectations for multiple generations of consanguineous unions. We identified potentially pathogenic rare exonic or splice site variants in 12 known (including KMT2C, SCN1A, SPTBN1, SYNE1, ZNF292) and 12 candidate (including CHD5, GRB10, PPP1R13B) ASD genes. Furthermore, we annotated noncoding variants in ROHs with brain-specific regulatory elements and identified putative disease-causing variants within brain-specific promoters and enhancers for 5 known ASD and neurodevelopmental disease genes (ACTG1, AUTS2, CTNND2, CNTNAP4, SPTBN4). We also identified copy number variants in two known ASD and neurodevelopmental disease loci in two affected individuals. In total we identified potentially etiological variants in known ASD or neurodevelopmental disease genes for ~61% (14/23) of affected individuals. We combined WGS with homozygosity mapping and regulatory element annotations to identify candidate ASD variants. Our analyses add to the growing number of ASD genes and variants and emphasize the importance of leveraging recent shared ancestry to map disease variants in complex neurodevelopmental disorders.

Abstract Background Reproducible detection of inherited variants with whole genome sequencing (WGS) is vital for the implementation of precision medicine and is a complicated process in which each step affects variant call quality. Systematically assessing reproducibility of inherited variants with WGS and impact of each step in the process is needed for understanding and improving quality of inherited variants from WGS. Results To dissect the impact of factors involved in detection of inherited variants with WGS, we sequence triplicates of eight DNA samples representing two populations on three short-read sequencing platforms using three library kits in six labs and call variants with 56 combinations of aligners and callers. We find that bioinformatics pipelines (callers and aligners) have a larger impact on variant reproducibility than WGS platform or library preparation. Single-nucleotide variants (SNVs), particularly outside difficult-to-map regions, are more reproducible than small insertions and deletions (indels), which are least reproducible when > 5 bp. Increasing sequencing coverage improves indel reproducibility but has limited impact on SNVs above 30×. Conclusions Our findings highlight sources of variability in variant detection and the need for improvement of bioinformatics pipelines in the era of precision medicine with WGS.

Abstract Copy number variants (CNVs) are recognized as a crucial genetic cause of neurodevelopmental disorders (NDDs). Chromosomal microarray analysis (CMA), the first-tier diagnostic test for individuals with NDDs, has been utilized to detect CNVs in clinical practice, but most reports are still from populations of European ancestry. To contribute more worldwide clinical genomics data, we investigated the genetic etiology of 410 Han Chinese patients with NDDs (151 with autism and 259 with unexplained intellectual disability (ID) and developmental delay (DD)) using CMA (Affymetrix) after G-banding karyotyping. Among all the NDD patients, 109 (26.6%) carried clinically relevant CNVs or uniparental disomies (UPDs), and 8 (2.0%) had aneuploidies (6 with trisomy 21 syndrome, 1 with 47,XXY, 1 with 47,XYY). In total, we found 129 clinically relevant CNVs and UPDs, including 32 CNVs in 30 ASD patients, and 92 CNVs and 5 UPDs in 79 ID/DD cases. When excluding the eight patients with aneuploidies, the diagnostic yield of pathogenic and likely pathogenic CNVs and UPDs was 20.9% for all NDDs (84/402), 3.3% in ASD (5/151), and 31.5% in ID/DD (79/251). When aneuploidies were included, the diagnostic yield increased to 22.4% for all NDDs (92/410), and 33.6% for ID/DD (87/259). We identified a de novo CNV in 14.9% (60/402) of subjects with NDDs. Interestingly, a higher diagnostic yield was observed in females (31.3%, 40/128) compared to males (16.1%, 44/274) for all NDDs (P = 4.8 × 10−4), suggesting that a female protective mechanism exists for deleterious CNVs and UPDs.

The genetic basis of autism spectrum disorder (ASD) is known to consist of contributions from de novo mutations in variant-intolerant genes. We hypothesize that rare inherited structural variants in cis-regulatory elements (CRE-SVs) of these genes also contribute to ASD. We investigated this by assessing the evidence for natural selection and transmission distortion of CRE-SVs in whole genomes of 9274 subjects from 2600 families affected by ASD. In a discovery cohort of 829 families, structural variants were depleted within promoters and untranslated regions, and paternally inherited CRE-SVs were preferentially transmitted to affected offspring and not to their unaffected siblings. The association of paternal CRE-SVs was replicated in an independent sample of 1771 families. Our results suggest that rare inherited noncoding variants predispose children to ASD, with differing contributions from each parent.

Gene therapy by engineering patient's own blood cells to confer HIV resistance can potentially lead to a functional cure for AIDS. Toward this goal, we have previously developed an anti-HIV lentivirus vector that deploys a combination of shRNA, ribozyme and RNA decoy. To further improve this therapeutic vector against viral escape, we sought an additional reagent to target HIV integrase. Here, we report the development of a new strategy for selection and expression of aptamer for gene therapy. We developed a SELEX protocol (multi-tag SELEX) for selecting RNA aptamers against proteins with low solubility or stability, such as integrase. More importantly, we expressed these aptamers in vivo by incorporating them in the terminal loop of shRNAs. This novel strategy allowed efficient expression of the shRNA-aptamer fusions that targeted RNAs and proteins simultaneously. Expressed shRNA-aptamer fusions targeting HIV integrase or reverse transcriptase inhibited HIV replication in cell cultures. Viral inhibition was further enhanced by combining an anti-integrase aptamer with an anti-HIV Tat-Rev shRNA. This construct exhibited efficacy comparable to that of integrase inhibitor Raltegravir. Our strategy for the selection and expression of RNA aptamers can potentially extend to other gene therapy applications.

Abstract Background Overlapping genes (OLGs) with long protein-coding overlapping sequences are disallowed by standard genome annotation programs, outside of viruses. Recently however they have been discovered in Archaea, diverse Bacteria, and Mammals. The biological factors underlying life’s ability to create overlapping genes require more study, and may have important applications in understanding evolution and in biotechnology. A previous study claimed that protein domains from viruses were much better suited to forming overlaps than those from other cellular organisms - in this study we assessed this claim, in order to discover what might underlie taxonomic differences in the creation of gene overlaps. Results After overlapping arbitrary Pfam domain pairs and evaluating them with Hidden Markov Models we find OLG construction to be much less constrained than expected. For instance, close to 10% of the constructed sequences cannot be distinguished from typical sequences in their protein family. Most are also indistinguishable from natural protein sequences regarding identity and secondary structure. Surprisingly, contrary to a previous study, virus domains were much less suitable for designing OLGs than bacterial or eukaryotic domains were. In general, the amount of amino acid change required to force a domain to overlap is approximately equal to the variation observed within a typical domain family. The resulting high similarity between natural sequences and those altered so as to overlap is mostly due to the combination of high redundancy in the genetic code and the evolutionary exchangeability of many amino acids. Conclusions Synthetic overlapping genes which closely resemble natural gene sequences, as measured by HMM profiles, are remarkably easy to construct, and most arbitrary domain pairs can be altered so as to overlap while retaining high similarity to the original sequences. Future work however will need to assess important factors not considered such as intragenic interactions which affect protein folding. While the analysis here is not sufficient to guarantee functional folding proteins, further analysis of constructed OLGs will improve our understanding of the origin of these remarkable genetic elements across life and opens up exciting possibilities for synthetic biology.

Abstract Background Cytosine modifications in DNA such as 5-methylcytosine (5mC) underlie a broad range of developmental processes, maintain cellular lineage specification, and can define or stratify types of cancer and other diseases. However, the wide variety of approaches available to interrogate these modifications has created a need for harmonized materials, methods, and rigorous benchmarking to improve genome-wide methylome sequencing applications in clinical and basic research. Here, we present a multi-platform assessment and cross-validated resource for epigenetics research from the FDA’s Epigenomics Quality Control Group. Results Each sample is processed in multiple replicates by three whole-genome bisulfite sequencing (WGBS) protocols (TruSeq DNA methylation, Accel-NGS MethylSeq, and SPLAT), oxidative bisulfite sequencing (TrueMethyl), enzymatic deamination method (EMSeq), targeted methylation sequencing (Illumina Methyl Capture EPIC), single-molecule long-read nanopore sequencing from Oxford Nanopore Technologies, and 850k Illumina methylation arrays. After rigorous quality assessment and comparison to Illumina EPIC methylation microarrays and testing on a range of algorithms (Bismark, BitmapperBS, bwa-meth, and BitMapperBS), we find overall high concordance between assays, but also differences in efficiency of read mapping, CpG capture, coverage, and platform performance, and variable performance across 26 microarray normalization algorithms. Conclusions The data provided herein can guide the use of these DNA reference materials in epigenomics research, as well as provide best practices for experimental design in future studies. By leveraging seven human cell lines that are designated as publicly available reference materials, these data can be used as a baseline to advance epigenomics research.

Abstract Autism Spectrum Disorder (ASD) is genetically complex with ~100 copy number variants and genes involved. To try to establish more definitive genotype and phenotype correlations in ASD, we searched genome sequence data, and the literature, for recurrent predicted damaging sequence-level variants affecting single genes. We identified 18 individuals from 16 unrelated families carrying a heterozygous guanine duplication (c.3679dup; p.Ala1227Glyfs*69) occurring within a string of 8 guanines (genomic location [hg38]g.50,721,512dup) affecting SHANK3, a prototypical ASD gene (0.08% of ASD-affected individuals carried the predicted p.Ala1227Glyfs*69 frameshift variant). Most probands carried de novo mutations, but five individuals in three families inherited it through somatic mosaicism. We scrutinized the phenotype of p.Ala1227Glyfs*69 carriers, and while everyone (17/17) formally tested for ASD carried a diagnosis, there was the variable expression of core ASD features both within and between families. Defining such recurrent mutational mechanisms underlying an ASD outcome is important for genetic counseling and early intervention.

Abstract In more than one-third of primary immunodeficiency (PID) patients, extensive genetic analysis including whole-exome sequencing (WES) fails to identify the genetic defect. Whole-genome sequencing (WGS) is able to detect variants missed by other genomics platforms, enabling the molecular diagnosis of otherwise unresolved cases. Here, we report two siblings, offspring of consanguineous parents, who experienced similar severe events encompassing early onset of colitis, lymphoproliferation, and hypogammaglobulinemia, typical of lipopolysaccharide-responsive and beige-like anchor (LRBA) or cytotoxic T lymphocyte antigen 4 (CTLA4) deficiencies. Gene-panel sequencing, comparative genomic hybridization (CGH) array, and WES failed to reveal a genetic aberration in relevant genes. WGS of these patients detected a 12.3 kb homozygous tandem duplication that was absent in control cohorts and is predicted to disrupt the reading frame of the LRBA gene. The variant was validated by PCR and Sanger sequencing, demonstrating the presence of the junction between the reference and the tandem-duplicated sequence. Droplet digital PCR (ddPCR) further confirmed the copy number in the unaffected parents (CN = 3, heterozygous) and affected siblings (CN = 4, homozygous), confirming the expected segregation pattern. In cases of suspected inherited immunodeficiency, WGS may reveal a mutation when other methods such as microarray and WES analysis failed to detect an aberration.

Abstract Background In recent years, several hundred autism spectrum disorder (ASD) implicated genes have been discovered impacting a wide range of molecular pathways. However, the molecular underpinning of ASD, particularly from the point of view of ‘brain to behaviour’ pathogenic mechanisms, remains largely unknown. Methods We undertook a study to investigate patterns of spatiotemporal and cell type expression of ASD-implicated genes by integrating large-scale brain single-cell transcriptomes (> million cells) and de novo loss-of-function (LOF) ASD variants (impacting 852 genes from 40,122 cases). Results We identified multiple single-cell clusters from three distinct developmental human brain regions (anterior cingulate cortex, middle temporal gyrus and primary visual cortex) that evidenced high evolutionary constraint through enrichment for brain critical exons and high pLI genes. These clusters also showed significant enrichment with ASD loss-of-function variant genes (p

ObjectiveChromosome 22q11.2 deletion syndrome (22q11.2DS) is associated with a more than 20-fold increased risk for developing schizophrenia. The aim of this study was to identify additional genetic factors (i.e., "second hits") that may contribute to schizophrenia expression.MethodThrough an international consortium, the authors obtained DNA samples from 329 psychiatrically phenotyped subjects with 22q11.2DS. Using a high-resolution microarray platform and established methods to assess copy number variation (CNV), the authors compared the genome-wide burden of rare autosomal CNV, outside of the 22q11.2 deletion region, between two groups: a schizophrenia group and those with no psychotic disorder at age ≥25 years. The authors assessed whether genes overlapped by rare CNVs were overrepresented in functional pathways relevant to schizophrenia.ResultsRare CNVs overlapping one or more protein-coding genes revealed significant between-group differences. For rare exonic duplications, six of 19 gene sets tested were enriched in the schizophrenia group; genes associated with abnormal nervous system phenotypes remained significant in a stepwise logistic regression model and showed significant interactions with 22q11.2 deletion region genes in a connectivity analysis. For rare exonic deletions, the schizophrenia group had, on average, more genes overlapped. The additional rare CNVs implicated known (e.g., GRM7, 15q13.3, 16p12.2) and novel schizophrenia risk genes and loci.ConclusionsThe results suggest that additional rare CNVs overlapping genes outside of the 22q11.2 deletion region contribute to schizophrenia risk in 22q11.2DS, supporting a multigenic hypothesis for schizophrenia. The findings have implications for understanding expression of psychotic illness and herald the importance of whole-genome sequencing to appreciate the overall genomic architecture of schizophrenia.

We systematically analyzed postzygotic mutations (PZMs) in whole-exome sequences from the largest collection of trios (5,947) with autism spectrum disorder (ASD) available, including 282 unpublished trios, and performed resequencing using multiple independent technologies. We identified 7.5% of de novo mutations as PZMs, 83.3% of which were not described in previous studies. Damaging, nonsynonymous PZMs within critical exons of prenatally expressed genes were more common in ASD probands than controls (P < 1 × 10-6), and genes carrying these PZMs were enriched for expression in the amygdala (P = 5.4 × 10-3). Two genes (KLF16 and MSANTD2) were significantly enriched for PZMs genome-wide, and other PZMs involved genes (SCN2A, HNRNPU and SMARCA4) whose mutation is known to cause ASD or other neurodevelopmental disorders. PZMs constitute a significant proportion of de novo mutations and contribute importantly to ASD risk.

ObjectiveWe aimed to generate a review and description of the phenotypic and genotypic spectra of ARHGEF9 mutations.MethodsPatients with mutations or chromosomal disruptions affecting ARHGEF9 were identified through our clinics and review of the literature. Detailed medical history and examination findings were obtained via a standardized questionnaire, or if this was not possible by reviewing the published phenotypic features.ResultsA total of 18 patients (including 5 females) were identified. Six had de novo, 5 had maternally inherited mutations, and 7 had chromosomal disruptions. All females had strongly skewed X-inactivation in favor of the abnormal X-chromosome. Symptoms presented in early childhood with delayed motor development alone or in combination with seizures. Intellectual disability was severe in most and moderate in patients with milder mutations. Males with severe intellectual disability had severe, often intractable, epilepsy and exhibited a particular facial dysmorphism. Patients with mutations in exon 9 affecting the protein's PH domain did not develop epilepsy.ConclusionsARHGEF9 encodes a crucial neuronal synaptic protein; loss of function of which results in severe intellectual disability, epilepsy, and a particular facial dysmorphism. Loss of only the protein's PH domain function is associated with the absence of epilepsy.

ObjectiveWe aimed to generate a review and description of the phenotypic and genotypic spectra of ARHGEF9 mutations.MethodsPatients with mutations or chromosomal disruptions affecting ARHGEF9 were identified through our clinics and review of the literature. Detailed medical history and examination findings were obtained via a standardized questionnaire, or if this was not possible by reviewing the published phenotypic features.ResultsA total of 18 patients (including 5 females) were identified. Six had de novo, 5 had maternally inherited mutations, and 7 had chromosomal disruptions. All females had strongly skewed X-inactivation in favor of the abnormal X-chromosome. Symptoms presented in early childhood with delayed motor development alone or in combination with seizures. Intellectual disability was severe in most and moderate in patients with milder mutations. Males with severe intellectual disability had severe, often intractable, epilepsy and exhibited a particular facial dysmorphism. Patients with mutations in exon 9 affecting the protein's PH domain did not develop epilepsy.ConclusionsARHGEF9 encodes a crucial neuronal synaptic protein; loss of function of which results in severe intellectual disability, epilepsy, and a particular facial dysmorphism. Loss of only the protein's PH domain function is associated with the absence of epilepsy.

Variation in the expression level and activity of genes involved in drug disposition and action ('pharmacogenes') can affect drug response and toxicity, especially when in tissues of pharmacological importance. Previous studies have relied primarily on microarrays to understand gene expression differences, or have focused on a single tissue or small number of samples. The goal of this study was to use RNA-sequencing (RNA-seq) to determine the expression levels and alternative splicing of 389 Pharmacogenomics Research Network pharmacogenes across four tissues (liver, kidney, heart and adipose) and lymphoblastoid cell lines, which are used widely in pharmacogenomics studies. Analysis of RNA-seq data from 139 different individuals across the 5 tissues (20-45 individuals per tissue type) revealed substantial variation in both expression levels and splicing across samples and tissue types. Comparison with GTEx data yielded a consistent picture. This in-depth exploration also revealed 183 splicing events in pharmacogenes that were previously not annotated. Overall, this study serves as a rich resource for the research community to inform biomarker and drug discovery and use.

Abstract Collectively, rare genetic diseases affect a significant number of individuals worldwide. In this study, we have conducted whole-exome sequencing (WES) and identified underlying pathogenic or likely pathogenic variants in five children with rare genetic diseases. We present evidence for disease-causing autosomal recessive variants in a range of disease-associated genes such as DHH-associated 46,XY gonadal dysgenesis (GD) or 46,XY sex reversal 7, GNPTAB-associated mucolipidosis II alpha/beta (ML II), BBS1-associated Bardet–Biedl Syndrome (BBS), SURF1-associated Leigh Syndrome (LS) and AP4B1-associated spastic paraplegia-47 (SPG47) in unrelated affected members from Bangladesh. Our analysis pipeline detected three homozygous mutations, including a novel c. 863 G > C (p.Pro288Arg) variant in DHH, and two compound heterozygous variants, including two novel variants: c.2972dupT (p.Met991Ilefs*) in GNPTAB and c.229 G > C (p.Gly77Arg) in SURF1. All mutations were validated by Sanger sequencing. Collectively, this study adds to the genetic heterogeneity of rare genetic diseases and is the first report elucidating the genetic profile of (consanguineous and nonconsanguineous) rare genetic diseases in the Bangladesh population.

Osteopontin (OPN), encoded by SPP1, is a phosphorylated glycoprotein predominantly synthesized in kidney tissue. Increased OPN mRNA and protein expression correlates with proteinuria, reduced creatinine clearance, and kidney fibrosis in animal models of kidney disease. But its genetic underpinnings are incompletely understood. We therefore conducted a genome-wide association study (GWAS) of OPN in a European chronic kidney disease (CKD) population. Using data from participants of the German Chronic Kidney Disease (GCKD) study (N = 4,897), a GWAS (minor allele frequency [MAF]≥1%) and aggregated variant testing (AVT, MAF

A strong motivation for undertaking psychiatric gene discovery studies is to provide novel insights into unknown biology. Although attention-deficit hyperactivity disorder (ADHD) is highly heritable, and large, rare copy number variants (CNVs) contribute to risk, little is known about its pathogenesis and it remains commonly misunderstood. We assembled and pooled five ADHD and control CNV data sets from the United Kingdom, Ireland, United States of America, Northern Europe and Canada. Our aim was to test for enrichment of neurodevelopmental gene sets, implicated by recent exome-sequencing studies of (a) schizophrenia and (b) autism as a means of testing the hypothesis that common pathogenic mechanisms underlie ADHD and these other neurodevelopmental disorders. We also undertook hypothesis-free testing of all biological pathways. We observed significant enrichment of individual genes previously found to harbour schizophrenia de novo non-synonymous single-nucleotide variants (SNVs; P=5.4 × 10(-4)) and targets of the Fragile X mental retardation protein (P=0.0018). No enrichment was observed for activity-regulated cytoskeleton-associated protein (P=0.23) or N-methyl-D-aspartate receptor (P=0.74) post-synaptic signalling gene sets previously implicated in schizophrenia. Enrichment of ADHD CNV hits for genes impacted by autism de novo SNVs (P=0.019 for non-synonymous SNV genes) did not survive Bonferroni correction. Hypothesis-free testing yielded several highly significantly enriched biological pathways, including ion channel pathways. Enrichment findings were robust to multiple testing corrections and to sensitivity analyses that excluded the most significant sample. The findings reveal that CNVs in ADHD converge on biologically meaningful gene clusters, including ones now established as conferring risk of other neurodevelopmental disorders.

The new concept of mammalian sex maintenance establishes that particular key genes must remain active in the differentiated gonads to avoid genetic sex reprogramming, as described in adult ovaries after Foxl2 ablation. Dmrt1 plays a similar role in postnatal testes, but the mechanism of adult testis maintenance remains mostly unknown. Sox9 and Sox8 are required for postnatal male fertility, but their role in the adult testis has not been investigated. Here we show that after ablation of Sox9 in Sertoli cells of adult, fertile Sox8(-/-) mice, testis-to-ovary genetic reprogramming occurs and Sertoli cells transdifferentiate into granulosa-like cells. The process of testis regression culminates in complete degeneration of the seminiferous tubules, which become acellular, empty spaces among the extant Leydig cells. DMRT1 protein only remains in non-mutant cells, showing that SOX9/8 maintain Dmrt1 expression in the adult testis. Also, Sox9/8 warrant testis integrity by controlling the expression of structural proteins and protecting Sertoli cells from early apoptosis. Concluding, this study shows that, in addition to its crucial role in testis development, Sox9, together with Sox8 and coordinately with Dmrt1, also controls adult testis maintenance.

Abstract Repeat expansions are responsible for over 40 monogenic disorders, and undoubtedly more pathogenic repeat expansions remain to be discovered. Existing methods for detecting repeat expansions in short-read sequencing data require predefined repeat catalogs. Recent discoveries emphasize the need for methods that do not require pre-specified candidate repeats. To address this need, we introduce ExpansionHunter Denovo, an efficient catalog-free method for genome-wide repeat expansion detection. Analysis of real and simulated data shows that our method can identify large expansions of 41 out of 44 pathogenic repeats, including nine recently reported non-reference repeat expansions not discoverable via existing methods.

Abstract Background Recent advances in genomic technologies have facilitated genome-wide investigation of human genetic variations. However, most efforts have focused on the major populations, yet trio genomes of indigenous populations from Southeast Asia have been under-investigated. Results We analyzed the whole-genome deep sequencing data (~ 30×) of five native trios from Peninsular Malaysia and North Borneo, and characterized the genomic variants, including single nucleotide variants (SNVs), small insertions and deletions (indels) and copy number variants (CNVs). We discovered approximately 6.9 million SNVs, 1.2 million indels, and 9000 CNVs in the 15 samples, of which 2.7% SNVs, 2.3% indels and 22% CNVs were novel, implying the insufficient coverage of population diversity in existing databases. We identified a higher proportion of novel variants in the Orang Asli (OA) samples, i.e., the indigenous people from Peninsular Malaysia, than that of the North Bornean (NB) samples, likely due to more complex demographic history and long-time isolation of the OA groups. We used the pedigree information to identify de novo variants and estimated the autosomal mutation rates to be 0.81 × 10− 8 – 1.33 × 10− 8, 1.0 × 10− 9 – 2.9 × 10− 9, and ~ 0.001 per site per generation for SNVs, indels, and CNVs, respectively. The trio-genomes also allowed for haplotype phasing with high accuracy, which serves as references to the future genomic studies of OA and NB populations. In addition, high-frequency inherited CNVs specific to OA or NB were identified. One example is a 50-kb duplication in DEFA1B detected only in the Negrito trios, implying plausible effects on host defense against the exposure of diverse microbial in tropical rainforest environment of these hunter-gatherers. The CNVs shared between OA and NB groups were much fewer than those specific to each group. Nevertheless, we identified a 142-kb duplication in AMY1A in all the 15 samples, and this gene is associated with the high-starch diet. Moreover, novel insertions shared with archaic hominids were identified in our samples. Conclusion Our study presents a full catalogue of the genome variants of the native Malaysian populations, which is a complement of the genome diversity in Southeast Asians. It implies specific population history of the native inhabitants, and demonstrated the necessity of more genome sequencing efforts on the multi-ethnic native groups of Malaysia and Southeast Asia.

Abstract Background There is limited empirical data quantifying the utility of genetic testing for families of children with autism spectrum disorder (ASD) or related neurodevelopmental disorders (NDD). We assessed the utility of clinical chromosomal microarray analysis (CMA), defined by diagnostic yield and parental empowerment, in population‐based sample of parents of affected children; and explored child, family, and health services factors predictive of empowerment. Methods Participants were families of children undergoing diagnostic assessments, between 2016 and 2019. Diagnostic yield of CMA in affected children was determined. Parental empowerment was measured through adapted version of the Genetics Counseling Outcome Scale‐24. Parents completed questionnaires to capture child, family, and health service factors. Results The diagnostic yield of CMA was 2.8% for pathogenic variants. Parental empowerment was significantly correlated with family functioning and aspects of perceived family‐centeredness of care. The model accounted for 49.8% of the variation in parental empowerment, F (10,37) = 3.67, p = 0.002. After accounting for other predictors, parental perception of the provision of general information remained significantly associated with empowerment. Conclusion The informational needs of families play an important role in their empowerment during genetic testing. Meeting these needs and monitoring empowerment can aid genomic technologies integration in personalized healthcare for ASD/NDD.

PurposeRecurrent 15q13.3 deletions are enriched in multiple neurodevelopmental conditions including intellectual disability, autism, epilepsy, and schizophrenia. However, the 15q13.3 microdeletion syndrome remains ill-defined.MethodsWe systematically compiled all cases of 15q13.3 deletion published before 2014. We also examined three locally available cohorts to identify new adults with 15q13.3 deletions.ResultsWe identified a total of 246 cases (133 children, 113 adults) with deletions overlapping or within the 15q13.3 (breakpoint (BP)4-BP5) region, including seven novel adult cases from local cohorts. No BP4-BP5 deletions were identified in 23,838 adult controls. Where known, 15q13.3 deletions were typically inherited (85.4%) and disproportionately of maternal origin (P < 0.0001). Overall, 198 cases (121 children, 77 adults; 80.5%) had at least one neuropsychiatric diagnosis. Accounting for ascertainment, developmental disability/intellectual disability was present in 57.7%, epilepsy/seizures in 28.0%, speech problems in 15.9%, autism spectrum disorder in 10.9%, schizophrenia in 10.2%, mood disorder in 10.2%, and attention deficit hyperactivity disorder in 6.5%. By contrast, major congenital malformations, including congenital heart disease (2.4%), were uncommon. Placenta previa occurred in the pregnancies of four cases.ConclusionThe 15q13.3 microdeletion syndrome is predominantly characterized by neuropsychiatric expression. There are implications for pre- and postnatal detection, genetic counseling, and anticipatory care.

BackgroundThe chromodomain helicase DNA binding domain (CHD) proteins modulate gene expression via their ability to remodel chromatin structure and influence histone acetylation. Recent studies have shown that CHD2 protein plays a critical role in embryonic development, tumor suppression and survival. Like other genes encoding members of the CHD family, pathogenic mutations in the CHD2 gene are expected to be implicated in human disease. In fact, there is emerging evidence suggesting that CHD2 might contribute to a broad spectrum of neurodevelopmental disorders. Despite growing evidence, a description of the full phenotypic spectrum of this condition is lacking.MethodsWe conducted a multicentre study to identify and characterise the clinical features associated with haploinsufficiency of CHD2. Patients with deletions of this gene were identified from among broadly ascertained clinical cohorts undergoing genomic microarray analysis for developmental delay, congenital anomalies and/or autism spectrum disorder.ResultsDetailed clinical assessments by clinical geneticists showed recurrent clinical symptoms, including developmental delay, intellectual disability, epilepsy, behavioural problems and autism-like features without characteristic facial gestalt or brain malformations observed on magnetic resonance imaging scans. Parental analysis showed that the deletions affecting CHD2 were de novo in all four patients, and analysis of high-resolution microarray data derived from 26,826 unaffected controls showed no deletions of this gene.ConclusionsThe results of this study, in addition to our review of the literature, support a causative role of CHD2 haploinsufficiency in developmental delay, intellectual disability, epilepsy and behavioural problems, with phenotypic variability between individuals.

BackgroundValidating the potential pathogenicity of copy number variants (CNVs) identified in genome-wide studies of autism spectrum disorders (ASD) requires detailed assessment of case/control frequencies, inheritance patterns, clinical correlations, and functional impact. Here, we characterize a small recurrent duplication in the annexin A1 (ANXA1) gene, identified by the Autism Genome Project (AGP) study.MethodsFrom the AGP CNV genomic screen in 2,147 ASD individuals, we selected for characterization an ANXA1 gene duplication that was absent in 4,964 population-based controls. We further screened the duplication in a follow-up sample including 1,496 patients and 410 controls, and evaluated clinical correlations and family segregation. Sequencing of exonic/downstream ANXA1 regions was performed in 490 ASD patients for identification of additional variants.ResultsThe ANXA1 duplication, overlapping the last four exons and 3'UTR region, had an overall prevalence of 11/3,643 (0.30%) in unrelated ASD patients but was not identified in 5,374 controls. Duplication carriers presented no distinctive clinical phenotype. Family analysis showed neuropsychiatric deficits and ASD traits in multiple relatives carrying the duplication, suggestive of a complex genetic inheritance. Sequencing of exonic regions and the 3'UTR identified 11 novel changes, but no obvious variants with clinical significance.ConclusionsWe provide multilevel evidence for a role of ANXA1 in ASD etiology. Given its important role as mediator of glucocorticoid function in a wide variety of brain processes, including neuroprotection, apoptosis, and control of the neuroendocrine system, the results add ANXA1 to the growing list of rare candidate genetic etiological factors for ASD.

Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community.