Search

Background Eosinophilic esophagitis (EoE) is a chronic inflammatory disorder of the esophagus marked by eosinophilic infiltration. Cumulative evidence indicates that the risk of EoE involves the complex interplay of both genetic and environmental factors. Because only a few genetic loci have been identified in EoE, the genetic underpinning of EoE remains largely elusive. Objective We sought to identify genetic loci associated with EoE. Methods Four EoE cohorts were genotyped using the Illumina single nucleotide polymorphism array platform, totaling 1,930 cases and 13,634 controls of European ancestry. Genotype imputation was performed with the Michigan Imputation Server using the Trans-Omics for Precision Medicine reference panel including whole-genome sequencing data from more than 100,000 individuals. Meta-analysis was conducted to identify potential novel genetic loci associated with EoE. Results Our study identified 11 new genome-wide significant loci, of which 6 are common variant loci, including 5q31.1 (rs2106984, P = 4.16 × 10−8; odds ratio [OR], 1.26, RAD50), 15q22.2 (rs2279293, P = 1.23 × 10−10; OR, 0.69, RORA), and 15q23 (rs56062135, P = 2.91 × 10−11; OR, 1.29, SMAD3), which have been previously associated with allergic conditions. Interestingly, a low-frequency synonymous mutation within the MATN2 gene was identified as the most significant single nucleotide polymorphism at the 8q22.1 locus. We also identified 5 sex-specific loci in the EoE cases, including an inflammatory bowel disease–associated locus at 9p24.1 (rs62541556, P = 4.4 × 10−8; OR, 1.11, JAK2). Conclusions Our findings demonstrate shared genetic underpinnings between EoE and other immune-mediated diseases and provide novel candidate genes for therapeutic target identification and prioritization.

We report a 38-year-old Saudi male with Ehlers-Danlos Syndrome (EDS). The patient presented with rare and unusual neurological manifestations, including but not limited to ophthalmoplegia and myopathic pattern on his electromyography. In addition to hand weakness, there was skin hyperextensibility, joint hyperflexibility, and frontal baldness. Next-generation sequencing was performed on target exon sequences, using whole exome sequencing and Burrows-Wheeler Aligner for alignment/base calling. Genome Analysis Toolkit and reference genome Homo sapiens (UCSC hg19) were used for sequence processing and analysis. Variant classification was done according to standard international recommendations. A novel homozygous variant, NM_019105.6: c.8488C>T p.(Gln2830*), was detected in the TNXB gene. This variant is not reported in the literature nor dbSNP or gnomAD databases. Additionally, this variant is predicted to create a premature stop codon and produce a truncated protein or nonsense-mediated mRNA decay. Hence, it is classified as a likely pathogenic variant. The same point variant was found in a heterozygous state in the patient's father and sister. Both presented with milder symptoms associated with Ehlers-Danlos syndromes and heritable connective tissue disorders. Therefore, the patient was diagnosed as a tenascin-X (TNX) deficient type of EDS known as classical-like Ehlers-Danlos syndrome. TNX deficient patients may present with clinical and electrophysiological manifestations that are unusual in EDS like frontal baldness, ophthalmoplegia, and myotonia, which mimic myotonic dystrophy type I. Clinicians should be aware of the potential overlap of symptoms among these two diseases to ensure correct diagnosis is made.

This study aimed to evaluate the effect of micro-osteoperforations (MOPs) on the gene expression profile of the periodontal ligament (PDL) of orthodontically moved teeth. Fifteen participants were randomly assigned into two groups: tooth movement only (Tr1, n = 7) and tooth movement supplemented with MOPs (Tr2, n = 8). In each subject, orthodontic tooth movement (OTM) was performed on premolar in one side, while no force was applied on contralateral premolar (Unt, n = 15). Seven days after loading, premolars were extracted for orthodontic reasons. RNA extraction from PDL and subsequent RNA-sequencing were performed. False discovery rates (Padj

Summary Tumor-associated macrophages (TAMs) play an important role in tumor immunity and comprise of subsets that have distinct phenotype, function, and ontology. Transcriptomic analyses of human medulloblastoma, the most common malignant pediatric brain cancer, showed that medulloblastomas (MBs) with activated sonic hedgehog signaling (SHH-MB) have significantly more TAMs than other MB subtypes. Therefore, we examined MB-associated TAMs by single-cell RNA sequencing of autochthonous murine SHH-MB at steady state and under two distinct treatment modalities: molecular-targeted inhibitor and radiation. Our analyses reveal significant TAM heterogeneity, identify markers of ontologically distinct TAM subsets, and show the impact of brain microenvironment on the differentiation of tumor-infiltrating monocytes. TAM composition undergoes dramatic changes with treatment and differs significantly between molecular-targeted and radiation therapy. We identify an immunosuppressive monocyte-derived TAM subset that emerges with radiation therapy and demonstrate its role in regulating T cell and neutrophil infiltration in MB.

Current understanding of the underlying molecular network and mechanism for attention-deficit hyperactivity disorder (ADHD) is lacking and incomplete. Previous studies suggest that genomic structural variations play an important role in the pathogenesis of ADHD. For effective modeling, deep learning approaches have become a method of choice, with ability to predict the impact of genetic variations involving complicated mechanisms. In this study, we examined copy number variation in whole genome sequencing from 116 African Americans ADHD children and 408 African American controls. We divided the human genome into 150 regions, and the variation intensity in each region was applied as feature vectors for deep learning modeling to classify ADHD patients. The accuracy of deep learning for predicting ADHD diagnosis is consistently around 78% in a two-fold shuffle test, compared with ∼50% by traditional k-mean clustering methods. Additional whole genome sequencing data from 351 European Americans children, including 89 ADHD cases and 262 controls, were applied as independent validation using feature vectors obtained from the African American ethnicity analysis. The accuracy of ADHD labeling was lower in this setting (∼70–75%) but still above the results from traditional methods. The regions with highest weight overlapped with the previously reported ADHD-associated copy number variation regions, including genes such as GRM1 and GRM8, key drivers of metabotropic glutamate receptor signaling. A notable discovery is that structural variations in non-coding genomic (intronic/intergenic) regions show prediction weights that can be as high as prediction weight from variations in coding regions, results that were unexpected.

Plasma lipids are known heritable risk factors for cardiovascular disease, but increasing evidence also supports shared genetics with diseases of other organ systems. We devised a comprehensive three-phase framework to identify novel lipid-associated genes and study the relationships between lipids, genotypes, gene expression and hundreds of complex human diseases from electronic Medical Records and Genomics (347 traits) and UK Biobank (549 traits) cohorts. Aside from 67 novel lipid-associated genes with strong replication, we found evidence for pleiotropic SNPs/genes between lipids and diseases across the phenome. These include discordant pleiotropy in the HLA region between lipids and multiple sclerosis and putative causal paths between triglycerides and gout, among several others. Our findings give insights into the genetic basis of the relationship between plasma lipids and diseases on a phenome-wide scale and can provide context for future prevention and treatment strategies.

CLEC16A has been shown to play a role in autophagy/mitophagy processes. Additionally, genetic variants in CLEC16A have been implicated in multiple autoimmune diseases. We generated an inducible whole-body knockout, Clec16aΔUBC mice, to investigate the loss of function of CLEC16A. The mice exhibited a neuronal phenotype including tremors and impaired gait that rapidly progressed to dystonic postures. Nerve conduction studies and pathological analysis revealed loss of sensory axons that are associated with this phenotype. Activated microglia and astrocytes were found in regions of the CNS. Several mitochondrial-related proteins were up- or down-regulated. Upregulation of interferon stimulated gene 15 (IGS15) were observed in neuronal tissues. CLEC16A expression inversely related to IGS15 expression. ISG15 may be the link between CLEC16A and downstream autoimmune, inflammatory processes. Our results demonstrate that a whole-body, inducible knockout of Clec16a in mice results in an inflammatory neurodegenerative phenotype resembling spinocerebellar ataxia.

Background Not all cells in a given individual are identical in their genomic makeup. Mosaicism describes such a phenomenon where a mixture of genotypic states in certain genomic segments exists within the same individual. Mosaicism is a prevalent and impactful class of non-integer state copy number variation (CNV). Mosaicism implies that certain cell types or subset of cells contain a CNV in a segment of the genome while other cells in the same individual do not. Several studies have investigated the impact of mosaicism in single patients or small cohorts but no comprehensive scan of mosaic CNVs has been undertaken to accurately detect such variants and interpret their impact on human health and disease. Results We developed a tool called Montage to improve the accuracy of detection of mosaic copy number variants in a high throughput fashion. Montage directly interfaces with ParseCNV2 algorithm to establish disease phenotype genome-wide association and determine which genomic ranges had more or less than expected frequency of mosaic events. We screened for mosaic events in over 350,000 samples using 1% allele frequency as the detection limit. Additionally, we uncovered disease associations of multiple phenotypes with mosaic CNVs at several genomic loci. We additionally investigated the allele imbalance observations genome-wide to define non-diploid and non-integer copy number states. Conclusions Our novel algorithm presents an efficient tool with fast computational runtime and high levels of accuracy of mosaic CNV detection. A curated mosaic CNV callset of 3716 events in 2269 samples is presented with comparability to previous reports and disease phenotype associations. The new algorithm can be freely accessed via: https://github.com/CAG-CNV/MONTAGE.

Background Elevated triglycerides (TG) are associated with, and may be causal for, cardiovascular disease (CVD), and co-morbidities such as type II diabetes and metabolic syndrome. Pathogenic variants in APOA5 and APOC3 as well as risk SNVs in other genes [APOE (rs429358, rs7412), APOA1/C3/A4/A5 gene cluster (rs964184), INSR (rs7248104), CETP (rs7205804), GCKR (rs1260326)] have been shown to affect TG levels. Knowledge of genetic causes for elevated TG may lead to early intervention and targeted treatment for CVD. We previously identified linkage and association of a rare, highly conserved missense variant in SLC25A40, rs762174003, with hypertriglyceridemia (HTG) in a single large family, and replicated this association with rare, highly conserved missense variants in a European American and African American sample. Methods Here, we analyzed a longitudinal mixed-ancestry cohort (European, African and Asian ancestry, N = 8966) from the Electronic Medical Record and Genomics (eMERGE) Network. We tested associations between median TG and the genes of interest, using linear regression, adjusting for sex, median age, median BMI, and the first two principal components of ancestry. Results We replicated the association between TG and APOC3, APOA5, and risk variation at APOE, APOA1/C3/A4/A5 gene cluster, and GCKR. We failed to replicate the association between rare, highly conserved variation at SLC25A40 and TG, as well as for risk variation at INSR and CETP. Conclusions Analysis using data from electronic health records presents challenges that need to be overcome. Although large amounts of genotype data is becoming increasingly accessible, usable phenotype data can be challenging to obtain. We were able to replicate known, strong associations, but were unable to replicate moderate associations due to the limited sample size and missing drug information.

Neuropsychiatric disorders, such as autism spectrum disorder (ASD), attention deficit hyperactivity disorder (ADHD), schizophrenia (SCZ), bipolar disorder (BIP), and major depressive disorder (MDD) share common clinical presentations, suggesting etiologic overlap. A substantial proportion of SNP-based heritability for neuropsychiatric disorders is attributable to genetic components, and genome-wide association studies (GWASs) focusing on individual diseases have identified multiple genetic loci shared between these diseases. Here, we aimed at identifying novel genetic loci associated with individual neuropsychiatric diseases and genetic loci shared by neuropsychiatric diseases. We performed multi-trait joint analyses and meta-analysis across five neuropsychiatric disorders based on their summary statistics from the Psychiatric Genomics Consortium (PGC), and further carried out a replication study of ADHD among 2726 cases and 16299 controls in an independent pediatric cohort. In the multi-trait joint analyses, we found five novel genome-wide significant loci for ADHD, one novel locus for BIP, and ten novel loci for MDD. We further achieved modest replication in our independent pediatric dataset. We conducted fine-mapping and functional annotation through an integrative multi-omics approach and identified causal variants and potential target genes at each novel locus. Gene expression profile and gene-set enrichment analysis further suggested early developmental stage expression pattern and postsynaptic membrane compartment enrichment of candidate genes at the genome-wide significant loci of these neuropsychiatric disorders. Therefore, through a multi-omics approach, we identified novel genetic loci associated with the five neuropsychiatric disorders which may help to better understand the underlying molecular mechanism of neuropsychiatric diseases.

The receptor for colony stimulating factor 1 (CSF-1R) is important for the survival and function of myeloid cells that mediate pathology during experimental autoimmune encephalomyelitis (EAE), an animal model of multiple sclerosis (MS). CSF-1 and IL-34, the ligands of CSF-1R, have similar bioactivities but distinct tissue and context-dependent expression patterns, suggesting that they have different roles. This could be the case in EAE, given that CSF-1 expression is upregulated in the CNS, while IL-34 remains constitutively expressed. We found that targeting CSF-1 with neutralizing mAb halted ongoing EAE, with efficacy superior to CSF-1R inhibitor BLZ945, whereas IL-34 neutralization had no effect, suggesting that pathogenic myeloid cells were maintained by CSF-1, not IL-34. Both anti-CSF-1- and BLZ945-treatment greatly reduced numbers of monocyte-derived cells and microglia in the CNS. However, anti-CSF-1 selectively depleted inflammatory microglia and monocytes in inflamed CNS areas, whereas BLZ945 depleted virtually all myeloid cells, including quiescent microglia, throughout the CNS. Anti-CSF-1 treatments reduced the size of demyelinated lesions, and microglial activation in the grey matter. Lastly, we found that bone marrow-derived immune cells were the major mediators of CSF-1R-dependent pathology, while microglia played a lesser role. Our findings suggest that targeting CSF-1 could be effective in ameliorating MS pathology, while preserving the homeostatic functions of myeloid cells, thereby minimizing risks associated with total ablation of CSF-1R-dependent cells.Significance StatementMultiple sclerosis (MS) and its animal model, experimental autoimmune encephalomyelitis (EAE), are autoimmune diseases characterized by accumulation of myeloid immune cells into the central nervous system (CNS). Both harmful and beneficial myeloid cells are present in EAE/MS, and a goal of MS therapy is to preferentially remove harmful myeloid cells. The receptor for CSF-1 (CSF-1R) is found on myeloid cells and it is important for their survival. CSF-1R can bind two ligands, CSF-1 and IL-34, but is unknown whether their functions in EAE/MS differ. We found that blocking CSF-1 depleted only harmful myeloid cells in the CNS and suppressed EAE, whereas blocking IL-34 had no effect. Thus, we propose that blocking CSF-1 could be a novel therapy for MS.

Background/Objectives Melanocortin-4 receptor (MC4R) plays an essential role in food intake and energy homeostasis. More than 170 MC4R variants have been described over the past two decades, with conflicting reports regarding the prevalence and phenotypic effects of these variants in diverse cohorts. To determine the frequency of MC4R variants in large cohort of different ancestries, we evaluated the MC4R coding region for 20,537 eMERGE participants with sequencing data plus additional 77,454 independent individuals with genome-wide genotyping data at this locus. Subjects/Methods The sequencing data were obtained from the eMERGE phase III study, in which multisample variant call format calls have been generated, curated, and annotated. In addition to penetrance estimation using body mass index (BMI) as a binary outcome, GWAS and PheWAS were performed using median BMI in linear regression analyses. All results were adjusted for principal components, age, sex, and sites of genotyping. Results Targeted sequencing data of MC4R revealed 125 coding variants in 1839 eMERGE participants including 30 unreported coding variants that were predicted to be functionally damaging. Highly penetrant unreported variants included (L325I, E308K, D298N, S270F, F261L, T248A, D111V, and Y80F) in which seven participants had obesity class III defined as BMI ≥ 40 kg/m2. In GWAS analysis, in addition to known risk haplotype upstream of MC4R (best variant rs6567160 (P = 5.36 × 10−25, Beta = 0.37), a novel rare haplotype was detected which was protective against obesity and encompassed the V103I variant with known gain-of-function properties (P = 6.23 × 10−08, Beta = −0.62). PheWAS analyses extended this protective effect of V103I to type 2 diabetes, diabetic nephropathy, and chronic renal failure independent of BMI. Conclusions MC4R screening in a large eMERGE cohort confirmed many previous findings, extend the MC4R pleotropic effects, and discovered additional MC4R rare alleles that probably contribute to obesity.

Previous studies of attention-deficit hyperactivity disorder (ADHD) have suggested that structural variants (SVs) play an important role but these were mainly studied in subjects of European ancestry and focused on coding regions. In this study, we sought to address the role of SVs in non-European populations and outside of coding regions. To that end, we generated whole genome sequence (WGS) data on 875 individuals, including 205 ADHD cases and 670 non-ADHD controls. The ADHD cases included 116 African Americans (AA) and 89 of European Ancestry (EA) with SVs in comparison with 408 AA and 262 controls, respectively. Multiple SVs and target genes that associated with ADHD from previous studies were identified or replicated, and novel recurrent ADHD-associated SV loci were discovered. We identified clustering of non-coding SVs around neuroactive ligand-receptor interaction pathways, which are involved in neuronal brain function, and highly relevant to ADHD pathogenesis and regulation of gene expression related to specific ADHD phenotypes. There was little overlap (around 6%) in the genes impacted by SVs between AA and EA. These results suggest that SVs within non-coding regions may play an important role in ADHD development and that WGS could be a powerful discovery tool for studying the molecular mechanisms of ADHD

Research on DNA methylation on N6-adenine (6mA) in eukaryotes has received much recent attention. Recent studies have generated a large amount of 6mA genomic data, yet the role of DNA 6mA in eukaryotes remains elusive, or even controversial. We argue that the sparsity of DNA 6mA in eukaryotes, the limitations of current biotechnologies for 6mA detection and the sophistication of the 6mA regulatory mechanism together pose great challenges for elucidation of DNA 6mA. To exploit existing 6mA genomic data and address this challenge, here we develop a deep-learning-based algorithm for predicting potential DNA 6mA sites de novo from sequence at single-nucleotide resolution, with application to three representative model organisms, Arabidopsis thaliana, Drosophila melanogaster and Escherichia coli. Extensive experiments demonstrate the accuracy of our algorithm and its superior performance compared with conventional k-mer-based approaches. Furthermore, our saliency maps-based context analysis protocol reveals interesting cis-regulatory patterns around the 6mA sites that are missed by conventional motif analysis. Our proposed analytical tools and findings will help to elucidate the regulatory mechanisms of 6mA and benefit the in-depth exploration of their functional effects. Finally, we offer a complete catalogue of potential 6mA sites based on in silico whole-genome prediction. The role of DNA methylation on N6-adenine (6mA) in eukaryotes is a challenging research problem. Tan et al. develop a deep-learning-based algorithm to predict 6mA sites from sequences at single-nucleotide resolution, and apply the method to three representative model organisms. The method is further developed to visualize regulatory patterns around 6mA sites.

The intestine is a highly dynamic environment that requires tight control of the various inputs to maintain homeostasis and allow for proper responses to injury. It was recently found that the stem cell niche and epithelium is regenerated after injury by de-differentiated adult cells, through a process that gives rise to Sca1+ fetal-like cells and is driven by a transient population of Clu+ revival stem cells (revSCs). However, the molecular mechanisms that regulate this dynamic process have not been fully defined. Here we show that TNFAIP8 (also known as TIPE0) is a regulator of intestinal homeostasis that is vital for proper regeneration. TIPE0 functions through inhibiting basal Akt activation by the commensal microbiota via modulating membrane phospholipid abundance. Loss of TIPE0 in mice results in injury-resistant enterocytes, that are hyperproliferative, yet have regenerative deficits and are shifted towards a de-differentiated state. Tipe0−/− enterocytes show basal induction of the Clu+ regenerative program and a fetal gene expression signature marked by Sca1, but upon injury are unable to generate Sca-1+/Clu+ revSCs and could not regenerate the epithelium. This work demonstrates the role of TIPE0 in regulating the dynamic signaling that determines the injury response and enables intestinal epithelial cell regenerative plasticity., The molecular mechanisms that regulate intestinal Clu+ revival stem cells (revSCs) and their niche to enable regeneration in response to injury are unclear. Here, the authors show that mice without the phospholipid transport protein, TNFAIP8, causes less revSCs to be induced following injury.

Background and Aims Among the >240 genetic loci described to date which confer susceptibility to inflammatory bowel disease, a small subset have been fine-mapped to an individual, non-coding single nucleotide polymorphism [SNP]. To illustrate a model mechanism by which a presumed-causal non-coding SNP can function, we analysed rs1887428, located in the promoter region of the Janus kinase 2 [JAK2] gene. Methods We utilized comparative affinity purification-mass spectrometry, DNA–protein binding assays, CRISPR/Cas9 genome editing, transcriptome sequencing and methylome quantitative trait locus methods to characterize the role of this SNP. Results We determined that the risk allele of rs1887428 is bound by the transcription factor [TF] RBPJ, while the protective allele is bound by the homeobox TF CUX1. While rs188748 only has a very modest influence on JAK2 expression, this effect was amplified downstream through the expression of pathway member STAT5B and epigenetic modification of the JAK2 locus. Conclusion Despite the absence of a consensus TF-binding motif or expression quantitative trait locus, we have used improved methods to characterize a putatively causal SNP to yield insight into inflammatory bowel disease mechanisms. Podcast This article has an associated podcast which can be accessed at https://academic.oup.com/ecco-jcc/pages/podcast

Baseline lung function is a standard diagnostic criterion used by clinicians to detect lung diseases. It is a complex trait significantly influenced by both genetics and environmental factors..., Baseline lung function, quantified as forced expiratory volume in the first second of exhalation (FEV1), is a standard diagnostic criterion used by clinicians to identify and classify lung diseases. Using whole-genome sequencing data from the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine project, we identified a novel genetic association with FEV1 on chromosome 12 in 867 African American children with asthma (P = 1.26 × 10−8, β = 0.302). Conditional analysis within 1 Mb of the tag signal (rs73429450) yielded one major and two other weaker independent signals within this peak. We explored statistical and functional evidence for all variants in linkage disequilibrium with the three independent signals and yielded nine variants as the most likely candidates responsible for the association with FEV1. Hi-C data and expression QTL analysis demonstrated that these variants physically interacted with KITLG (KIT ligand, also known as SCF), and their minor alleles were associated with increased expression of the KITLG gene in nasal epithelial cells. Gene-by-air-pollution interaction analysis found that the candidate variant rs58475486 interacted with past-year ambient sulfur dioxide exposure (P = 0.003, β = 0.32). This study identified a novel protective genetic association with FEV1, possibly mediated through KITLG, in African American children with asthma. This is the first study that has identified a genetic association between lung function and KITLG, which has established a role in orchestrating allergic inflammation in asthma.

BackgroundRubinstein–Taybi syndrome (RTS) is a rare, congenital, plurimalformative, and neurodevelopmental disorder. Previous studies have reported that large deletions contribute to more severe RTS phenotypes than those caused by CREBBP point mutations, suggesting a concurrent pathogenetic role of flanking genes, typical of contiguous gene syndromes, but the detailed genetics are unclear.ResultsThis study presented a rare case of Rubinstein-Taybi (RT) syndrome with serious cardiac abnormalities. Based on the clinical and genetic analysis of the patient, theADCY9gene deletion was highlighted as a plausible explanation of cardiac abnormalities. Inadcy9morphant zebrafish, cardiac malformation was observed. Immunofluorescence study disclosed increased macrophage migration and cardiac apoptosis. RNA sequencing in zebrafish model highlighted the changes of a number of genes, including increased expression of the mmp9 gene which encodes a matrix metalloproteinase with the main function to degrade and remodel extracellular matrix.ConclusionsIn this study, we identified a plausible new candidate geneADCY9of CHD through the clinical and genetic analysis of a rare case of Rubinstein-Taybi (RT) syndrome with serious cardiac abnormalities. By functional study of zebrafish, we demonstrated that deletion ofadcy9is the causation for the cardiac abnormalities. Cardiac apoptosis and increased expression of the MMP9 gene are involved in the pathogenesis.

Background Neuroblastoma is a childhood malignancy that arises from the developing sympathetic nervous system. Although mitochondrial dysfunctions have been implicated in the pathophysiology of neuroblastoma, the role of mitochondrial DNA (mtDNA) has not been extensively investigated. Methods A total of 2404 Caucasian children diagnosed with neuroblastoma and 9310 ancestry-matched controls were recruited at the Children’s Hospital of Philadelphia. The mtDNA haplogroups were identified from SNP array data of two independent cohorts. We conducted a case-control study to explore potential associations of mtDNA haplogroups with the susceptibility of neuroblastoma. The genetic effect of neuroblastoma was measured by odds ratios (ORs) of mitochondrial haplogroups. All tests were two-sided. Results Haplogroup K was statistically significantly associated with reduced risk of neuroblastoma in the discovery cohort consisting of 1474 cases and 5699 controls (OR = 0.72, 95% confidence interval [CI] = 0.57 to 0.90; P = 4.8 × 10-3). The association was replicated in an independent cohort (OR = 0.69, 95% CI = 0.53 to 0.92; P = .01) of 930 cases and 3611 controls. Pooled analysis was performed by combining the two data sets. The association remained highly statistically significant after correction for multiple testing (OR = 0.71, 95% CI = 0.59 to 0.84, P = 1.96 × 10-4, Pcorrected = .002). Further analysis focusing on neuroblastoma subtypes indicated haplogroup K was more associated with high-risk neuroblastoma (OR = 0.57, 95% CI = 0.43 to 0.76; P = 1.46 × 10–4) than low-risk and intermediate-risk neuroblastoma. Conclusions Haplogroup K is an independent genetic factor associated with reduced risk of developing neuroblastoma in European descents. These findings provide new insights into the genetic basis of neuroblastoma, implicating mitochondrial DNA encoded proteins in the etiology of neuroblastoma.

Most replicated genetic determinants for type 1 diabetes are common (minor allele frequency [MAF] >5%). We aimed to identify novel rare or low-frequency (MAF

Linked-read sequencing provides long-range information on short-read sequencing data by barcoding reads originating from the same DNA molecule, and can improve detection and breakpoint identification for structural variants (SVs). Here we present LinkedSV for SV detection on linked-read sequencing data. LinkedSV considers barcode overlapping and enriched fragment endpoints as signals to detect large SVs, while it leverages read depth, paired-end signals and local assembly to detect small SVs. Benchmarking studies demonstrate that LinkedSV outperforms existing tools, especially on exome data and on somatic SVs with low variant allele frequencies. We demonstrate clinical cases where LinkedSV identifies disease-causal SVs from linked-read exome sequencing data missed by conventional exome sequencing, and show examples where LinkedSV identifies SVs missed by high-coverage long-read sequencing. In summary, LinkedSV can detect SVs missed by conventional short-read and long-read sequencing approaches, and may resolve negative cases from clinical genome/exome sequencing studies., Compared to single nucleotide variants and short indels, structural variants (SVs) are often more challenging to detect using high-throughput sequencing based methods. Here, the authors develop LinkedSV, a computational tool for SV detection using linked-read exome and genome sequencing data.

Hepatitis is an inflammatory condition of the liver, which is frequently caused by the infection of hepatitis B virus (HBV) or hepatitis C virus (HCV). Hepatitis can lead to the development of chronic complications including cancer, making it a major public health burden. Co-infection of HBV and HCV can result in faster disease progression. Therefore, it is important to identify shared genetic susceptibility loci for HBV and HCV infection to further understand the underlying mechanism. Through a meta-analysis based on genome-wide association summary statistics of HBV and HCV infection, we found one novel locus in the Asian population and two novel loci in the European population. By functional annotation based on multi-omics data, we identified the likely target genes at each novel locus, such as HMGB1 and ATF3, which play a critical role in autophagy and immune response to virus. By re-analyzing a microarray dataset from Hmgb1–/– mice and RNA-seq data from mouse liver tissue overexpressing ATF3, we found that differential expression of autophagy and immune and metabolic gene pathways underlie these conditions. Our study reveals novel common susceptibility loci to HBV and HCV infection, supporting their role in linking autophagy signaling and immune response.

With polygenic risk score (PRS) for autoimmune type 1 diabetes (T1D), this study identified T1D cases with low T1D PRS and searched for susceptibility loci in these cases. Our hypothesis is that genetic effects (likely mediated by relatively rare genetic variants) of non-mainstream (or non-autoimmune) T1D might have been diluted in the previous studies on T1D cases in general. Two cohorts for the PRS modeling and testing respectively were included. The first cohort consisted of 3302 T1D cases and 6181 controls, and the independent second cohort consisted of 3297 T1D cases and 6169 controls. Cases with low T1D PRS were identified using PRSice-2 and compared to controls with low T1D PRS by genome-wide association (GWA) test. Thirteen novel genetic loci with high imputation quality (Quality Score r2 > 0.91) were identified of SNPs/SNVs associated with low PRS T1D at genome-wide significance (P ≤ 5.0 × E−08), in addition to 4 established T1D loci, 3 reported loci by our previous study, as well as 9 potential novel loci represented by rare SNVs, but with relatively low imputation quality (Quality Score r2

Conotruncal defects with normally related great vessels (CTD-NRGVs) occur in both patients with and without 22q11.2 deletion syndrome (22q11.2DS), but it is unclear to what extent the genetically complex etiologies of these heart defects may overlap across these two groups, potentially involving variation within and/or outside of the 22q11.2 region. To explore this potential overlap, we conducted genome-wide SNP-level, gene-level, and gene set analyses using common variants, separately in each of five cohorts, including two with 22q11.2DS (N = 1472 total cases) and three without 22q11.2DS (N = 935 total cases). Results from the SNP-level analyses were combined in meta-analyses, and summary statistics from these analyses were also used in gene and gene set analyses. Across all these analyses, no association was significant after correction for multiple comparisons. However, several SNPs, genes, and gene sets with suggestive evidence of association were identified. For common inherited variants, we did not identify strong evidence for shared genomic mechanisms for CTD-NRGVs across individuals with and without 22q11.2 deletions. Nevertheless, several of our top gene-level and gene set results have been linked to cardiogenesis and may represent candidates for future work.

Assumptions are made about the genetic model of single nucleotide polymorphisms (SNPs) when choosing a traditional genetic encoding: additive, dominant, and recessive. Furthermore, SNPs across the genome are unlikely to demonstrate identical genetic models. However, running SNP-SNP interaction analyses with every combination of encodings raises the multiple testing burden. Here, we present a novel and flexible encoding for genetic interactions, the elastic data-driven genetic encoding (EDGE), in which SNPs are assigned a heterozygous value based on the genetic model they demonstrate in a dataset prior to interaction testing. We assessed the power of EDGE to detect genetic interactions using 29 combinations of simulated genetic models and found it outperformed the traditional encoding methods across 10%, 30%, and 50% minor allele frequencies (MAFs). Further, EDGE maintained a low false-positive rate, while additive and dominant encodings demonstrated inflation. We evaluated EDGE and the traditional encodings with genetic data from the Electronic Medical Records and Genomics (eMERGE) Network for five phenotypes: age-related macular degeneration (AMD), age-related cataract, glaucoma, type 2 diabetes (T2D), and resistant hypertension. A multi-encoding genome-wide association study (GWAS) for each phenotype was performed using the traditional encodings, and the top results of the multi-encoding GWAS were considered for SNP-SNP interaction using the traditional encodings and EDGE. EDGE identified a novel SNP-SNP interaction for age-related cataract that no other method identified: rs7787286 (MAF: 0.041; intergenic region of chromosome 7)–rs4695885 (MAF: 0.34; intergenic region of chromosome 4) with a Bonferroni LRT p of 0.018. A SNP-SNP interaction was found in data from the UK Biobank within 25 kb of these SNPs using the recessive encoding: rs60374751 (MAF: 0.030) and rs6843594 (MAF: 0.34) (Bonferroni LRT p: 0.026). We recommend using EDGE to flexibly detect interactions between SNPs exhibiting diverse action., Author summary Although traditional genetic encodings are widely implemented in genetics research, including in genome-wide association studies (GWAS) and epistasis, each method makes assumptions that may not reflect the underlying etiology. Here, we introduce a novel encoding method that estimates and assigns an individualized data-driven encoding for each single nucleotide polymorphism (SNP): the elastic data-driven genetic encoding (EDGE). With simulations, we demonstrate that this novel method is more accurate and robust than traditional encoding methods in estimating heterozygous genotype values, reducing the type I error, and detecting SNP-SNP interactions. We further applied the traditional encodings and EDGE to biomedical data from the Electronic Medical Records and Genomics (eMERGE) Network for five phenotypes, and EDGE identified a novel interaction for age-related cataract not detected by traditional methods, which replicated in data from the UK Biobank. EDGE provides an alternative approach to understanding and modeling diverse SNP models and is recommended for studying complex genetics in common human phenotypes.

NKAP is a ubiquitously expressed nucleoplasmic protein that is currently known as a transcriptional regulatory molecule via its interaction with HDAC3 and spliceosomal proteins. Here, we report a disorder of transcriptional regulation due to missense mutations in the X chromosome gene, NKAP. These mutations are clustered in the C-terminal region of NKAP where NKAP interacts with HDAC3 and post-catalytic spliceosomal complex proteins. Consistent with a role for the C-terminal region of NKAP in embryogenesis, nkap mutant zebrafish with a C-terminally truncated NKAP demonstrate severe developmental defects. The clinical features of affected individuals are highly conserved and include developmental delay, hypotonia, joint contractures, behavioral abnormalities, Marfanoid habitus, and scoliosis. In affected cases, transcriptome analysis revealed the presence of a unique transcriptome signature, which is characterized by the downregulation of long genes with higher exon numbers. These observations indicate the critical role of NKAP in transcriptional regulation and demonstrate that perturbations of the C-terminal region lead to developmental defects in both humans and zebrafish.

Research suggested accumulation of tau proteins might lead to the degeneration of functional networks. Studies investigating the impact of genetic risk for Alzheimer's disease (AD) on early brain connections might shed light on mechanisms leading to AD development later in life. Here, we aim to investigate whether the polygenic risk score for Alzheimer's disease (AD-PRS) influences the connectivity among regions susceptible to tau pathology during childhood and adolescence. Participants were youth, aged 6–14 years, and recruited in Porto Alegre (discovery sample, n = 332) and Sao Paulo (replication sample, n = 304), Brazil. Subjects underwent genotyping and 6-min resting state funcional magnetic resonance imaging. Connections between the local maxima of tau pathology networks were used as dependent variables. The AD-PRS was associated with the connectivity between the right precuneus and the right superior temporal gyrus (discovery sample: β = 0.180, padjusted = 0.036; replication sample: β = 0.202, p = 0.031). This connectivity was also associated with inhibitory control (β = 0.157, padjusted = 0.035) and moderated the association between the AD-PRS and both immediate and delayed recall. These findings suggest the AD-PRS may affect brain connectivity in youth, which might impact memory performance and inhibitory control in early life.

The duration of pregnancy is influenced by fetal and maternal genetic and non-genetic factors. Here we report a fetal genome-wide association meta-analysis of gestational duration, and early preterm, preterm, and postterm birth in 84,689 infants. One locus on chromosome 2q13 is associated with gestational duration; the association is replicated in 9,291 additional infants (combined P = 3.96 × 10−14). Analysis of 15,588 mother-child pairs shows that the association is driven by fetal rather than maternal genotype. Functional experiments show that the lead SNP, rs7594852, alters the binding of the HIC1 transcriptional repressor. Genes at the locus include several interleukin 1 family members with roles in pro-inflammatory pathways that are central to the process of parturition. Further understanding of the underlying mechanisms will be of great public health importance, since giving birth either before or after the window of term gestation is associated with increased morbidity and mortality., Gestational duration depends on both maternal and fetal genetic influences. Here, the authors perform a fetal genome-wide association meta-analysis and find that a locus on 2q13 is associated with pregnancy duration and further show that the lead SNP rs7594852 changes the binding properties of transcriptional repressor HIC1.

Objective Up to 40% of patients with epilepsy become drug resistant (DRE). Genetic factors are likely to play a role. While efforts have focused on the transporter and target hypotheses, neither of them fully explains the pan-pharmacoresistance seen in DRE. Materials and methods In this study, we developed and used a phenotyping algorithm for the identification of DRE, responders, and epilepsy-free controls that were sequenced using a gene panel developed by the Pharmacogenomics Research Network (PGRN), which includes 82 genes involved in drug response. We tested the transporter hypothesis of DRE, the association between drug resistance and variants in the ATP-binding cassette family of genes previously associated with DRE, and also investigated potential new genetic factors. Results In the analysis of DRE vs controls, NTRK2 was significantly associated with DRE (rs76950094; P = 1.19 × 10-7 and gene-based P-value = 1.67 × 10-4 ). NTRK2 encodes TrkB, which is involved in the development and maturation of the central nervous system, and increased activation of TrkB signaling is suggested to promote epilepsy. Conclusion Although the role of NTRK2 in DRE needs to be elucidated, these results support alternative mechanisms underlying DRE, complementary to the existing hypotheses, that should be evaluated.

Purpose:To provide a validated method to confidently identify exon-containing copy number variants (CNVs), with a low false discovery rate (FDR), in targeted sequencing data from a clinical laboratory with particular focus on single-exon CNVs.Methods:DNA sequence coverage data are normalized within each sample and subsequently exonic CNVs are identified in a batch of samples (midpool), when the target log2 ratio of the sample to the batch median exceeds defined thresholds. The quality of exonic CNV calls is assessed by C-scores (Z-like scores) using thresholds derived from gold standard samples and simulation studies. We integrate an ExonQC threshold to lower FDR and compare performance with alternate software (VisCap).Results:Thirteen CNVs were used as a truth set to validate Atlas-CNV and compared with VisCap. We demonstrated FDR reduction in validation, simulation and 10,926 eMERGESeq samples without sensitivity loss. Sixty-four multi-exon and 29 single-exon CNVs with high C-scores were assessed by MLPA.Conclusions:Atlas-CNV is validated as a method to identify exonic CNVs in targeted sequencing data generated in the clinical laboratory. The ExonQC and C-score assignment can reduce FDR (identification of targets with high variance) and improve calling accuracy of single-exon CNVs respectively. We proposed guidelines and criteria to identify high confidence single-exon CNVs.

The study of patients with schizophrenia (SZ) at different clinical stages may help clarify what effects could be due to the disease itself, to the pharmacological treatment, or to the disease progression. We compared expression levels of targeted genes in blood from individuals in different stages of SZ: clinical high risk for psychosis (CHR), first episode of psychosis (FEP), and chronic SZ (CSZ). Then, we further verified whether single-nucleotide polymorphisms (SNPs) could be related to gene expression differences. We investigated 12 genes in 394 individuals (27 individuals with CHR, 70 antipsychotic-naive individuals with FEP, 157 CSZ patients, and 140 healthy controls (HCs)). For a subsample, genotype data were also available, and we extracted SNPs that were previously associated with the expression of selected genes in whole blood or brain tissue. We generated a mediation model in which a putative cause (SNP) is related to a presumed effect (disorder) via an intermediate variable (gene expression). MBP and NDEL1 were upregulated in FEP compared to all other groups; DGCR8 was downregulated in FEP compared to HC and CHR; DGCR2 was downregulated in CSZ compared to FEP and HCs; DISC1 was upregulated in schizophrenia compared to controls or FEP, possibly induced by the rs3738398 and rs10864693 genotypes, which were associated with DISC1 expression; and UFD1 was upregulated in CSZ and CHR compared to FEP and HC. Our results indicated changes in gene expression profiles throughout the different clinical stages of SZ, reinforcing the need for staging approaches to better capture SZ heterogeneity.

Variability in SARS-CoV-2 susceptibility and COVID-19 disease severity between individuals is partly due to genetic factors. Here, we identify 4 genomic loci with suggestive associations for SARS-CoV-2 susceptibility and 19 for COVID-19 disease severity. Four of these 23 loci likely have an ethnicity-specific component. Genome-wide association study (GWAS) signals in 11 loci colocalize with expression quantitative trait loci (eQTLs) associated with the expression of 20 genes in 62 tissues/cell types (range: 1:43 tissues/gene), including lung, brain, heart, muscle, and skin as well as the digestive system and immune system. We perform genetic fine mapping to compute 99% credible SNP sets, which identify 10 GWAS loci that have eight or fewer SNPs in the credible set, including three loci with one single likely causal SNP. Our study suggests that the diverse symptoms and disease severity of COVID-19 observed between individuals is associated with variants across the genome, affecting gene expression levels in a wide variety of tissue types., Graphical abstract, D’Antonio et al. characterize associations between GWAS signals for COVID-19 disease and eQTLs in 69 human tissues to identify causal variants and their underlying molecular mechanisms. They show that diverse symptoms and disease severity of COVID-19 are associated with variants affecting gene expression in a wide variety of tissues.

Hispanic/Latinos have been underrepresented in genome-wide association studies (GWAS) for anthropometric traits despite notable anthropometric variability with ancestry proportions, and a high burden of growth stunting and overweight/obesity in Hispanic/Latino populations. This address this knowledge gap, we analyzed densely-imputed genetic data in a sample of Hispanic/Latino adults, to identify and fine-map common genetic variants associated with body mass index (BMI), height, and BMI-adjusted waist-to-hip ratio (WHRadjBMI). We conducted a GWAS of 18 studies/consortia as part of the Hispanic/Latino Anthropometry (HISLA) Consortium (Stage 1, n=59,769) and validated our findings in 9 additional studies (HISLA Stage 2, n=9,336). We conducted a trans-ethnic GWAS with summary statistics from HISLA Stage 1 and existing consortia of European and African ancestries. In our HISLA Stage 1+2 analyses, we discovered one novel BMI locus, as well two novel BMI signals and another novel height signal, each within established anthropometric loci. In our trans-ethnic meta- analysis, we identified three additional novel BMI loci, one novel height locus, and one novel WHRadjBMI locus. We also identified three secondary signals for BMI, 28 for height, and two for WHRadjBMI. We replicated >60 established anthropometric loci in Hispanic/Latino populations at genome-wide significance—representing up to 30% of previously-reported index SNP anthropometric associations. Trans-ethnic meta-analysis of the three ancestries showed a small-to-moderate impact of uncorrected population stratification on the resulting effect size estimates. Our novel findings demonstrate that future studies may also benefit from leveraging differences in linkage disequilibrium patterns to discover novel loci and additional signals with less residual population stratification.

Exome sequencing and Drosophila modeling revealed variants in SMARCA5 as a novel cause of a neurodevelopmental disorder., Intellectual disability encompasses a wide spectrum of neurodevelopmental disorders, with many linked genetic loci. However, the underlying molecular mechanism for more than 50% of the patients remains elusive. We describe pathogenic variants in SMARCA5, encoding the ATPase motor of the ISWI chromatin remodeler, as a cause of a previously unidentified neurodevelopmental disorder, identifying 12 individuals with de novo or dominantly segregating rare heterozygous variants. Accompanying phenotypes include mild developmental delay, frequent postnatal short stature and microcephaly, and recurrent dysmorphic features. Loss of function of the SMARCA5 Drosophila ortholog Iswi led to smaller body size, reduced sensory dendrite complexity, and tiling defects in larvae. In adult flies, Iswi neural knockdown caused decreased brain size, aberrant mushroom body morphology, and abnormal locomotor function. Iswi loss of function was rescued by wild-type but not mutant SMARCA5. Our results demonstrate that SMARCA5 pathogenic variants cause a neurodevelopmental syndrome with mild facial dysmorphia.

The pioneer transcription factor (TF) PU.1 controls hematopoietic cell fate by decompacting stem cell heterochromatin and allowing nonpioneer TFs to enter otherwise inaccessible genomic sites. PU.1 deficiency fatally arrests lymphopoiesis and myelopoiesis in mice, but human congenital PU.1 disorders have not previously been described. We studied six unrelated agammaglobulinemic patients, each harboring a heterozygous mutation (four de novo, two unphased) of SPI1, the gene encoding PU.1. Affected patients lacked circulating B cells and possessed few conventional dendritic cells. Introducing disease-similar SPI1 mutations into human hematopoietic stem and progenitor cells impaired early in vitro B cell and myeloid cell differentiation. Patient SPI1 mutations encoded destabilized PU.1 proteins unable to nuclear localize or bind target DNA. In PU.1-haploinsufficient pro–B cell lines, euchromatin was less accessible to nonpioneer TFs critical for B cell development, and gene expression patterns associated with the pro– to pre–B cell transition were undermined. Our findings molecularly describe a novel form of agammaglobulinemia and underscore PU.1’s critical, dose-dependent role as a hematopoietic euchromatin gatekeeper.

The traditional view is that the occurrence and development of hallux valgus (HV) are mainly due to environmental factors. Recent studies have suggested the large contribution of genetic heritability to HV, but it remains elusive about the genetic variants underlying the development of HV. To gain knowledge about the molecular mechanisms of HV pathogenesis by genetic approach, whole exome sequencing studies were performed in 10 individuals (7 affected by HV and 3 unaffected) from three independent families. Specific mutations were found to be related to the pathogenesis of HV and conform to the laws of inheritance. A total of 36 genes with functional candidate single nucleotide variants were identified. Genetic predisposition plays an important role in the development of HV. Interestingly, some of these genes are related to chronic arthritis, such as the complement encoding gene C7, or are related to long toe or long fingers, such as TTN, COL6A3, LARS, FIG4, and CBS. This study identified rare potentially pathogenic mutations represented by genes related to digital anomalies and chronic arthritis underlying the familial types of HV, which acquired new insights into the genetic and physiological foundations of HV, thereby might improve accurate prevention and drug development for HV.

Oral‐facial‐digital syndromes (OFDS) are a heterogeneous and rare group of Mendelian disorders characterized by developmental abnormalities of the oral cavity, face, and digits caused by dysfunction of the primary cilium, a mechanosensory organelle that exists atop most cell types that facilitates organ patterning and growth. OFDS is inherited both in an X‐linked dominant, X‐linked recessive, and autosomal recessive manner. Importantly, though many of the causal genes for OFDS have been identified, up to 40% of OFD syndromes are of unknown genetic basis. Here we describe three children with classical presentations of OFDS including lingual hamartomas, polydactyly, and characteristic facial features found by exome sequencing to harbor variants in causal genes not previously associated with OFDS. We describe a female with hypothalamic hamartoma, urogenital sinus, polysyndactyly, and multiple lingual hamartomas consistent with OFDVI with biallelic pathogenic variants in CEP164, a gene associated with ciliopathy‐spectrum disease, but never before with OFDS. We additionally describe two unrelated probands with postaxial polydactyly, multiple lingual hamartomas, and dysmorphic features both found to be homozygous for an identical TOPORS missense variant, c.29 C>A; (p.Pro10Gln). Heterozygous TOPORS pathogenic gene variants are associated with autosomal dominant retinitis pigmentosa, but never before with syndromic ciliopathy. Of note, both probands are of Dominican ancestry, suggesting a possible founder allele.

Mutations affecting the transcriptional regulator ANKRD11 are mainly associated with the multisystem developmental disorder known as KBG syndrome, but have also been identified in individuals with Cornelia de Lange syndrome (CdLS) and other developmental disorders caused by variants affecting different chromatin regulators. The extensive functional overlap of these proteins results in shared phenotypical features, which complicate the assessment of the clinical diagnosis. Additionally, re-evaluation of individuals at a later age occasionally reveals that the initial phenotype has evolved towards clinical features more reminiscent of a developmental disorder different from the one that was initially diagnosed. For this reason, variants in ANKRD11 can be ascribed to a broader class of disorders that fall within the category of the so-called chromatinopathies. In this work, we report on the clinical characterization of 23 individuals with variants in ANKRD11. The subjects present primarily with developmental delay, intellectual disability and dysmorphic features, and all but two received an initial clinical diagnosis of either KBG syndrome or CdLS. The number and the severity of the clinical signs are overlapping but variable and result in a broad spectrum of phenotypes, which could be partially accounted for by the presence of additional molecular diagnoses and distinct pathogenic mechanisms. This article is protected by copyright. All rights reserved.

There is strong evidence for a genetic contribution to non-syndromic congenital heart defects (CHDs). However, exome- and genome-wide studies conducted at the variant and gene-level have identified few genome-wide significant CHD-related genes. Gene-set analyses are a useful complement to such studies and candidate gene-set analyses of rare variants have provided insight into the genetics of CHDs. However, similar analyses have not been conducted using data on common genetic variants. Consequently, we conducted common variant analyses of 15 CHD candidate gene-sets, using data from two common types of CHDs: conotruncal heart defects (1431 cases) and left ventricular outflow tract defects (509 cases). After Bonferroni correction for evaluation of multiple gene-sets, the cytoskeletal gene-set was significantly associated with conotruncal heart defects (βS = 0.09, 95% confidence interval (CI) 0.03–0.15). This association was stronger when analyses were restricted to the sub-set of cytoskeletal genes that have been observed to harbor rare damaging genotypes in at least two CHD cases (βS = 0.32, 95% CI 0.08–0.56). These findings add to the evidence linking cytoskeletal genes to CHDs and suggest that, for cytoskeletal genes, common variation may contribute to the risk of CHDs.

RNA-seq has been a powerful method to detect the differentially expressed genes/long non-coding RNAs (lncRNAs) in schizophrenia (SCZ) patients, however, due to overfitting problems differentially expressed targets (DETs) cannot be used properly as biomarkers. This study used machine learning to reduce gene/non-coding RNA features. Dorsolateral prefrontal cortex (dlpfc) RNA-seq data from 254 individuals was obtained from the CommonMind consortium. The average predictive accuracy for SCZ patients was 67% based on coding genes, and 96% based on long non-coding RNAs (lncRNAs). Machine learning is a powerful algorithm to reduce functional biomarkers in SCZ patients. The lncRNAs capture the characteristics of SCZ tissue more accurately than mRNA as the former regulate every level of gene expression, not limited to mRNA levels.

Attention-deficit hyperactivity disorder (ADHD) is a neurodevelopmental disorder with poorly understood molecular mechanisms that results in significant impairment in children. In this study, we sought to assess the role of rare recurrent variants in non-European populations and outside of coding regions. We generated whole genome sequence (WGS) data on 875 individuals, including 205 ADHD cases and 670 non-ADHD controls. The cases included 116 African Americans (AA) and 89 European Americans (EA), and the controls included 408 AA and 262 EA. Multiple novel rare recurrent variants were identified in exonic regions, functionally classified as stop-gains and frameshifts for known ADHD genes. Deletion in introns of the protocadherins families and the ncRNA HGB8P were identified in two independent EA ADHD patients. A meta-analysis of the two ethnicities for differential ADHD recurrent variants compared to controls shows a small number of overlaps. These results suggest that rare recurrent variants in noncoding regions may be involved in the pathogenesis of ADHD in children of both AA and EA ancestry, thus, WGS could be a powerful discovery tool for studying the molecular mechanisms of ADHD.

SUMMARYRare deletions and duplications of genomic segments, collectively known as rare copy number variants (rCNVs), contribute to a broad spectrum of human diseases. To date, most disease-association studies of rCNVs have focused on recognized genomic disorders or on the impact of haploinsufficiency caused by deletions. By comparison, our understanding of duplications in disease remains rudimentary as very few individual genes are known to be triplosensitive (i.e., duplication intolerant). In this study, we meta-analyzed rCNVs from 753,994 individuals across 30 primarily neurological disease phenotypes to create a genome-wide catalog of rCNV association statistics across disorders. We discovered 114 rCNV-disease associations at 52 distinct loci surpassing genome-wide significance (P=3.72×10−6), 42% of which involve duplications. Using Bayesian fine-mapping methods, we further prioritized 38 novel triplosensitive disease genes (e.g., GMEB2 in brain abnormalities), including three known haploinsufficient genes that we now reveal as bidirectionally dosage sensitive (e.g., ANKRD11 in growth abnormalities). By integrating our results with prior literature, we found that disease-associated rCNV segments were enriched for genes constrained against damaging coding variation and identified likely dominant driver genes for about one-third (32%) of rCNV segments based on de novo mutations from exome sequencing studies of developmental disorders. However, while the presence of constrained driver genes was a common feature of many pathogenic large rCNVs across disorders, most of the rCNVs showing genome-wide significant association were incompletely penetrant (mean odds ratio=11.6) and we also identified two examples of noncoding disease-associated rCNVs (e.g., intronic CADM2 deletions in behavioral disorders). Finally, we developed a statistical model to predict dosage sensitivity for all genes, which defined 3,006 haploinsufficient and 295 triplosensitive genes where the effect sizes of rCNVs were comparable to deletions of genes constrained against truncating mutations. These dosage sensitivity scores classified disease genes across molecular mechanisms, prioritized pathogenic de novo rCNVs in children with autism, and revealed features that distinguished haploinsufficient and triplosensitive genes, such as insulation from other genes and local cis-regulatory complexity. Collectively, the cross-disorder rCNV maps and metrics derived in this study provide the most comprehensive assessment of dosage sensitive genomic segments and genes in disease to date and set the foundation for future studies of dosage sensitivity throughout the human genome.

BackgroundBone accrual impacts lifelong skeletal health, but genetic discovery has been primarily limited to cross-sectional study designs and hampered by uncertainty about target effector genes. Here, we capture this dynamic phenotype by modeling longitudinal bone accrual across 11,000 bone scans in a cohort of healthy children and adolescents, followed by genome-wide association studies (GWAS) and variant-to-gene mapping with functional follow-up.ResultsWe identify 40 loci, 35 not previously reported, with various degrees of supportive evidence, half residing in topological associated domains harboring known bone genes. Of several loci potentially associated with later-life fracture risk, a candidate SNP lookup provides the most compelling evidence for rs11195210 (SMC3). Variant-to-gene mapping combining ATAC-seq to assay open chromatin with high-resolution promoter-focused Capture C identifies contacts between GWAS loci and nearby gene promoters. siRNA knockdown of gene expression supports the putative effector gene at three specific loci in two osteoblast cell models. Finally, using CRISPR-Cas9 genome editing, we confirm that the immediate genomic region harboring the putative causal SNP influences PRPF38A expression, a location which is predicted to coincide with a set of binding sites for relevant transcription factors.ConclusionsUsing a new longitudinal approach, we expand the number of genetic loci putatively associated with pediatric bone gain. Functional follow-up in appropriate cell models finds novel candidate genes impacting bone accrual. Our data also raise the possibility that the cell fate decision between osteogenic and adipogenic lineages is important in normal bone accrual.

Teebi hypertelorism syndrome (THS; OMIM 145420) is a rare craniofacial disorder characterized by hypertelorism, prominent forehead, short nose with broad or depressed nasal root. Some cases of THS have been attributed to SPECC1L variants. Homozygous variants in CDH11 truncating the transmembrane and intracellular domains have been implicated in Elsahy-Waters syndrome (EWS; OMIM 211380) with hypertelorism. We report THS due to CDH11 heterozygous missense variants on 19 subjects from 9 families. All affected residues in the extracellular region of Cadherin-11 (CHD11) are highly conserved across vertebrate species and classical cadherins. Six of the variants that cluster around the EC2-EC3 and EC3-EC4 linker regions are predicted to affect Ca2+ binding that is required for cadherin stability. Two of the additional variants [c.164G > C, p.(Trp55Ser) and c.418G > A, p.(Glu140Lys)] are also notable as they are predicted to directly affect trans-homodimer formation. Immunohistochemical study demonstrates that CDH11 is strongly expressed in human facial mesenchyme. Using multiple functional assays, we show that five variants from the EC1, EC2-EC3 linker, and EC3 regions significantly reduced the cell-substrate trans adhesion activity and one variant from EC3-EC4 linker results in changes in cell morphology, focal adhesion, and migration, suggesting dominant negative effect. Characteristic features in this cohort included depressed nasal root, cardiac and umbilical defects. These features distinguished this phenotype from that seen in SPECC1L-related hypertelorism syndrome and CDH11-related EWS. Our results demonstrate heterozygous variants in CDH11, which decrease cell-cell adhesion and increase cell migratory behavior, cause a form of THS, as termed CDH11-related THS.

EBioMedicine 72, 103588 (2021). doi:10.1016/j.ebiom.2021.103588, Published by Elsevier, Amsterdam [u.a.]

BackgroundGenetics and biology may influence the age of onset of anorexia nervosa (AN). The aims of this study were to determine whether common genetic variation contributes to age of onset of AN and to investigate the genetic associations between age of onset of AN and age at menarche.MethodsA secondary analysis of the Psychiatric Genomics Consortium genome-wide association study (GWAS) of AN was performed, which included 9335 cases and 31,981 screened controls, all from European ancestries. We conducted GWASs of age of onset, early-onset AN (ResultsTwo loci were genome-wide significant in the typical-onset AN GWAS. Heritability estimates (single nucleotide polymorphism–h2) were 0.01–0.04 for age of onset, 0.16–0.25 for early-onset AN, and 0.17–0.25 for typical-onset AN. Early- and typical-onset AN showed distinct genetic correlation patterns with putative risk factors for AN. Specifically, early-onset AN was significantly genetically correlated with younger age at menarche, and typical-onset AN was significantly negatively genetically correlated with anthropometric traits. Genetic risk scores for age of onset and early-onset AN estimated from independent GWASs significantly predicted age of onset. Mendelian randomization analysis suggested a causal link between younger age at menarche and early-onset AN.ConclusionsOur results provide evidence consistent with a common variant genetic basis for age of onset and implicate biological pathways regulating menarche and reproduction. Background: Genetics and biology may influence the age of onset of anorexia nervosa (AN). The aims of this study were to determine whether common genetic variation contributes to age of onset of AN and to investigate the genetic associations between age of onset of AN and age at menarche. Methods: A secondary analysis of the Psychiatric Genomics Consortium genome-wide association study (GWAS) of AN was performed, which included 9335 cases and 31,981 screened controls, all from European ancestries. We conducted GWASs of age of onset, early-onset AN (2) were 0.01–0.04 for age of onset, 0.16–0.25 for early-onset AN, and 0.17–0.25 for typical-onset AN. Early- and typical-onset AN showed distinct genetic correlation patterns with putative risk factors for AN. Specifically, early-onset AN was significantly genetically correlated with younger age at menarche, and typical-onset AN was significantly negatively genetically correlated with anthropometric traits. Genetic risk scores for age of onset and early-onset AN estimated from independent GWASs significantly predicted age of onset. Mendelian randomization analysis suggested a causal link between younger age at menarche and early-onset AN. Conclusions: Our results provide evidence consistent with a common variant genetic basis for age of onset and implicate biological pathways regulating menarche and reproduction.

Spermatogenesis-associated 5 like 1 (SPATA5L1) represents an orphan gene encoding a protein of unknown function. We report 28 bi-allelic variants in SPATA5L1 associated with sensorineural hearing loss in 47 individuals from 28 (26 unrelated) families. In addition, 25/47 affected individuals (53%) presented with microcephaly, developmental delay/intellectual disability, cerebral palsy, and/or epilepsy. Modeling indicated damaging effect of variants on the protein, largely via destabilizing effects on protein domains. Brain imaging revealed diminished cerebral volume, thin corpus callosum, and periventricular leukomalacia, and quantitative volumetry demonstrated significantly diminished white matter volumes in several individuals. Immunofluorescent imaging in rat hippocampal neurons revealed localization of Spata511 in neuronal and glial cell nuclei and more prominent expression in neurons. In the rodent inner ear, Spata511 is expressed in the neurosensory hair cells and inner ear supporting cells. Transcriptomic analysis performed with fibroblasts from affected individuals was able to distinguish affected from controls by principal components. Analysis of differentially expressed genes and networks suggested a role for SPATA5L1 in cell surface adhesion receptor function, intracellular focal adhesions, and DNA replication and mitosis. Collectively, our results indicate that bi-allelic SPATA5L1 variants lead to a human disease characterized by sensorineural hearing loss (SNHL) with or without a nonprogressive mixed neurodevelopmental phenotype.

Summary Both mild and severe epilepsies are influenced by variants in the same genes, yet an explanation for the resulting phenotypic variation is unknown. As part of the ongoing Epi25 Collaboration, we performed a whole-exome sequencing analysis of 13,487 epilepsy-affected individuals and 15,678 control individuals. While prior Epi25 studies focused on gene-based collapsing analyses, we asked how the pattern of variation within genes differs by epilepsy type. Specifically, we compared the genetic architectures of severe developmental and epileptic encephalopathies (DEEs) and two generally less severe epilepsies, genetic generalized epilepsy and non-acquired focal epilepsy (NAFE). Our gene-based rare variant collapsing analysis used geographic ancestry-based clustering that included broader ancestries than previously possible and revealed novel associations. Using the missense intolerance ratio (MTR), we found that variants in DEE-affected individuals are in significantly more intolerant genic sub-regions than those in NAFE-affected individuals. Only previously reported pathogenic variants absent in available genomic datasets showed a significant burden in epilepsy-affected individuals compared with control individuals, and the ultra-rare pathogenic variants associated with DEE were located in more intolerant genic sub-regions than variants associated with non-DEE epilepsies. MTR filtering improved the yield of ultra-rare pathogenic variants in affected individuals compared with control individuals. Finally, analysis of variants in genes without a disease association revealed a significant burden of loss-of-function variants in the genes most intolerant to such variation, indicating additional epilepsy-risk genes yet to be discovered. Taken together, our study suggests that genic and sub-genic intolerance are critical characteristics for interpreting the effects of variation in genes that influence epilepsy.

Contains fulltext : 234992.pdf (Publisher’s version ) (Closed access) PURPOSE: MED12 is a subunit of the Mediator multiprotein complex with a central role in RNA polymerase II transcription and regulation of cell growth, development, and differentiation. This might underlie the variable phenotypes in males carrying missense variants in MED12, including X-linked recessive Ohdo, Lujan, and FG syndromes. METHODS: By international matchmaking we assembled variant and clinical data on 18 females presenting with variable neurodevelopmental disorders (NDDs) and harboring de novo variants in MED12. RESULTS: Five nonsense variants clustered in the C-terminal region, two splice variants were found in the same exon 8 splice acceptor site, and 11 missense variants were distributed over the gene/protein. Protein truncating variants were associated with a severe, syndromic phenotype consisting of intellectual disability (ID), facial dysmorphism, short stature, skeletal abnormalities, feeding difficulties, and variable other abnormalities. De novo missense variants were associated with a less specific, but homogeneous phenotype including severe ID, autistic features, limited speech and variable other anomalies, overlapping both with females with truncating variants as well as males with missense variants. CONCLUSION: We establish de novo truncating variants in MED12 as causative for a distinct NDD and de novo missense variants as causative for a severe, less specific NDD in females.