Biliary atresia is associated with polygenic susceptibility in ciliogenesis and planar polarity effector genes.

BACKGROUND & AIMS: Biliary atresia (BA) is poorly understood and leads to liver transplantation (LT), with the requirement for and associated risks of lifelong immunosuppression, in most children. We performed a genome-wide association study (GWAS) to determine the genetic basis of BA. METHODS: We performed a GWAS in 811 European BA cases treated with LT in US, Canadian and UK centers, and 4,654 genetically matched controls. Whole-genome sequencing of 100 cases evaluated synthetic association with rare variants. Functional studies included whole liver transcriptome analysis of 64 BA cases and perturbations in experimental models. RESULTS: A GWAS of common single nucleotide polymorphisms (SNPs), i.e. allele frequencies >1%, identified intronic SNPs rs6446628 in AFAP1 with genome-wide significance (p = 3.93E-8) and rs34599046 in TUSC3 at sub-threshold genome-wide significance (p = 1.34E-7), both supported by credible peaks of neighboring SNPs. Like other previously reported BA-associated genes, AFAP1 and TUSC3 are ciliogenesis and planar polarity effectors (CPLANE). In gene-set-based GWAS, BA was associated with 6,005 SNPs in 102 CPLANE genes (p = 5.84E-15). Compared with non-CPLANE genes, more CPLANE genes harbored rare variants (allele frequency <1%) that were assigned Human Phenotype Ontology terms related to hepatobiliary anomalies by predictive algorithms, 87% vs. 40%, p <0.0001. Rare variants were present in multiple genes distinct from those with BA-associated common variants in most BA cases. AFAP1 and TUSC3 knockdown blocked ciliogenesis in mouse tracheal cells. Inhibition of ciliogenesis caused biliary dysgenesis in zebrafish. AFAP1 and TUSC3 were expressed in fetal liver organoids, as well as fetal and BA livers, but not in normal or disease-control livers. Integrative analysis of BA-associated variants and liver transcripts revealed abnormal vasculogenesis and epithelial tube formation, explaining portal vein anomalies that co-exist with BA. CONCLUSIONS: BA is associated with polygenic susceptibility in CPLANE genes. Rare variants contribute to polygenic risk in vulnerable pathways via unique genes. IMPACT AND IMPLICATIONS: Liver transplantation is needed to cure most children born with biliary atresia, a poorly understood rare disease. Transplant immunosuppression increases the likelihood of life-threatening infections and cancers. To improve care by preventing this disease and its progression to transplantation, we examined its genetic basis. We find that this disease is associated with both common and rare mutations in highly specialized genes which maintain normal communication and movement of cells, and their organization into bile ducts and blood vessels during early development of the human embryo. Because defects in these genes also cause other birth defects, our findings could lead to preventive strategies to lower the incidence of biliary atresia and potentially other birth defects.

Evidence of epistasis in regions of long-range linkage disequilibrium across five complex diseases in the UK Biobank and eMERGE datasets.

Leveraging linkage disequilibrium (LD) patterns as representative of population substructure enables the discovery of additive association signals in genome-wide association studies (GWASs). Standard GWASs are well-powered to interrogate additive models; however, new approaches are required for invesigating other modes of inheritance such as dominance and epistasis. Epistasis, or non-additive interaction between genes, exists across the genome but often goes undetected because of a lack of statistical power. Furthermore, the adoption of LD pruning as customary in standard GWASs excludes detection of sites that are in LD but might underlie the genetic architecture of complex traits. We hypothesize that uncovering long-range interactions between loci with strong LD due to epistatic selection can elucidate genetic mechanisms underlying common diseases. To investigate this hypothesis, we tested for associations between 23 common diseases and 5,625,845 epistatic SNP-SNP pairs (determined by Ohta's D statistics) in long-range LD (>0.25 cM). Across five disease phenotypes, we identified one significant and four near-significant associations that replicated in two large genotype-phenotype datasets (UK Biobank and eMERGE). The genes that were most likely involved in the replicated associations were (1) members of highly conserved gene families with complex roles in multiple pathways, (2) essential genes, and/or (3) genes that were associated in the literature with complex traits that display variable expressivity. These results support the highly pleiotropic and conserved nature of variants in long-range LD under epistatic selection. Our work supports the hypothesis that epistatic interactions regulate diverse clinical mechanisms and might especially be driving factors in conditions with a wide range of phenotypic outcomes.

Rare recurrent copy number variations in metabotropic glutamate receptor interacting genes in children with neurodevelopmental disorders.

BACKGROUND: Neurodevelopmental disorders (NDDs), such as attention deficit hyperactivity disorder (ADHD) and autism spectrum disorder (ASD), are examples of complex and partially overlapping phenotypes that often lack definitive corroborating genetic information. ADHD and ASD have complex genetic associations implicated by rare recurrent copy number variations (CNVs). Both of these NDDs have been shown to share similar biological etiologies as well as genetic pleiotropy. METHODS: Platforms aimed at investigating genetic-based associations, such as high-density microarray technologies, have been groundbreaking techniques in the field of complex diseases, aimed at elucidating the underlying disease biology. Previous studies have uncovered CNVs associated with genes within shared candidate genomic networks, including glutamate receptor genes, across multiple different NDDs. To examine shared biological pathways across two of the most common NDDs, we investigated CNVs across 15,689 individuals with ADHD (n = 7920), ASD (n = 4318), or both (n = 3,416), as well as 19,993 controls. Cases and controls were matched by genotype array (i.e., Illumina array versions). Three case-control association studies each calculated and compared the observed vs. expected frequency of CNVs across individual genes, loci, pathways, and gene networks. Quality control measures of confidence in CNV-calling, prior to association analyses, included visual inspection of genotype and hybridization intensity. RESULTS: Here, we report results from CNV analysis in search for individual genes, loci, pathways, and gene networks. To extend our previous observations implicating a key role of the metabotropic glutamate receptor (mGluR) network in both ADHD and autism, we exhaustively queried patients with ASD and/or ADHD for CNVs associated with the 273 genomic regions of interest within the mGluR gene network (genes with one or two degrees protein-protein interaction with mGluR 1-8 genes). Among CNVs in mGluR network genes, we uncovered CNTN4 deletions enriched in NDD cases (P = 3.22E - 26, OR = 2.49). Additionally, we uncovered PRLHR deletions in 40 ADHD cases and 12 controls (P = 5.26E - 13, OR = 8.45) as well as clinically diagnostic relevant 22q11.2 duplications and 16p11.2 duplications in 23 ADHD + ASD cases and 9 controls (P = 4.08E - 13, OR = 15.05) and 22q11.2 duplications in 34 ADHD + ASD cases and 51 controls (P = 9.21E - 9, OR = 3.93); those control samples were not with previous 22qDS diagnosis in their EHR records. CONCLUSION: Together, these results suggest that disruption in neuronal cell-adhesion pathways confers significant risk to NDDs and showcase that rare recurrent CNVs in CNTN4, 22q11.2, and 16p11.2 are overrepresented in NDDs that constitute patients predominantly suffering from ADHD and ASD. TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT02286817 First Posted: 10 November 14, ClinicalTrials.gov Identifier: NCT02777931 first posted: 19 May 2016, ClinicalTrials.gov Identifier: NCT03006367 first posted: 30 December 2016, ClinicalTrials.gov Identifier: NCT02895906 first posted: 12 September 2016.

Elevated Levels of the Cytokine LIGHT in Pediatric Crohn's Disease.

LIGHT (homologous to lymphotoxins, exhibits inducible expression, and competes with HSV glycoprotein D for herpes virus entry mediator, a receptor expressed by T lymphocytes), encoded by the TNFSF14 gene, is a cytokine belonging to the TNF superfamily. On binding to its receptors, herpes virus entry mediator and lymphotoxin ß receptor, it activates inflammatory responses. We conducted this study to determine whether plasma LIGHT levels are elevated in Crohn's disease (CD) in a pediatric population with the aim of nominating this cytokine as a therapeutic target. We used a single-molecule immunoassay to determine the circulating levels of free LIGHT in plasma from pediatric patients with CD in our biobank (n = 183), a panel of healthy pediatric (n = 9) or adult (n = 22) reference samples, and pediatric biobank controls (n = 19). We performed correlational analyses between LIGHT levels and the clinical characteristics of the CD cohort, including age, Montreal classification, family history, medical/surgical therapy, and routine blood test parameters. LIGHT levels were greatly elevated in CD, with an average of 305 versus 32.4 pg/ml for controls from the biobank (p < 0.0001). The outside reference samples showed levels of 57 pg/ml in pediatric controls and 55 pg/ml in adults (p < 0.0001). We found a statistically significant correlation between white blood cell count and free LIGHT (p < 0.046). We conclude that free, soluble LIGHT is increased 5- to 10-fold in pediatric CD across an array of disease subtypes and characteristics.

ParseCNV2: efficient sequencing tool for copy number variation genome-wide association studies.

Improved copy number variation (CNV) detection remains an area of heavy emphasis for algorithm development; however, both CNV curation and disease association approaches remain in its infancy. The current practice of focusing on candidate CNVs, where researchers study specific CNVs they believe to be pathological while discarding others, refrains from considering the full spectrum of CNVs in a hypothesis-free GWAS. To address this, we present a next-generation approach to CNV association by natively supporting the popular VCF specification for sequencing-derived variants as well as SNP array calls using a PennCNV format. The code is fast and efficient, allowing for the analysis of large (>100,000 sample) cohorts without dividing up the data on a compute cluster. The scripts are condensed into a single tool to promote simplicity and best practices. CNV curation pre and post-association is rigorously supported and emphasized to yield reliable results of highest quality. We benchmarked two large datasets, including the UK Biobank (n > 450,000) and CAG Biobank (n > 350,000) both of which are genotyped at >0.5 M probes, for our input files. ParseCNV has been actively supported and developed since 2008. ParseCNV2 presents a critical addition to formalizing CNV association for inclusion with SNP associations in GWAS Catalog. Clinical CNV prioritization, interactive quality control (QC), and adjustment for covariates are revolutionary new features of ParseCNV2 vs. ParseCNV. The software is freely available at: https://github.com/CAG-CNV/ParseCNV2 .

COVID-19 in pediatrics: Genetic susceptibility.

The uptick in SARS-CoV-2 infection has resulted in a worldwide COVID-19 pandemic, which has created troublesome health and economic problems. We performed case-control meta-analyses in both African and European ethnicity COVID-19 disease cases based on laboratory test and phenotypic criteria. The cases had laboratory-confirmed SARS-CoV-2 infection. We uniquely investigated COVID infection genetics in a pediatric population. Our cohort has a large African ancestry component, also unique to our study. We tested for genetic variant association in 498 cases vs. 1,533 controls of African ancestry and 271 cases vs. 855 controls of European ancestry. We acknowledge that the sample size is relatively small, owing to the low prevalence of COVID infection among pediatric individuals. COVID-19 cases averaged 13 years of age. Pediatric genetic studies enhance the ability to detect genetic associations with a limited possible environment impact. Our findings support the notion that some genetic variants, most notably at the SEMA6D, FMN1, ACTN1, PDS5B, NFIA, ADGRL3, MMP27, TENM3, SPRY4, MNS1, and RSU1 loci, play a role in COVID-19 infection susceptibility. The pediatric cohort also shows nominal replication of previously reported adult study results: CCR9, CXCR6, FYCO1, LZTFL1, TDGF1, CCR1, CCR2, CCR3, CCR5, MAPT-AS1, and IFNAR2 gene variants. Reviewing the biological roles of genes implicated here, NFIA looks to be the most interesting as it binds to a palindromic sequence observed in both viral and cellular promoters and in the adenovirus type 2 origin of replication.

Identification of novel loci in obstructive sleep apnea in European American and African American children.

STUDY OBJECTIVES: To identify genetic susceptibility variants in pediatric obstructive sleep apnea in European American and African American children. METHODS: A phenotyping algorithm using electronic medical records was developed to recruit cases with OSA and control subjects from the Center for Applied Genomics at Children's Hospital of Philadelphia (CHOP.) Genome-wide association studies (GWAS) were performed in pediatric OSA cases and control subjects with European American (EA) and African American (AA) ancestry followed by meta-analysis and sex stratification. RESULTS: The algorithm accrued 1,486 subjects (46.3% European American, 53.7% African American). We identified genomic loci at 1p36.22 and 15q26.1 that associated with OSA risk in EA and AA, respectively. We also revealed a shared risk locus at 18p11.32 (rs114124196, P =1.72 ×10 -8) across EA and AA populations. Additionally, association at 1q43 (rs12754698) and 2p25.1 (rs72775219) was identified in the male-only analysis of EA children with OSA, while association at 8q21.11 (rs6472959), 11q24.3 (rs4370952) and 15q21.1 (rs149936782) was detected in the female-only analysis of EA children and association at 18p11.23 (rs9964029) was identified in the female-only analysis of African-American children. Moreover, the 18p11.32 locus was replicated in an EA cohort (rs114124196, P =8.8 ×10 -3). CONCLUSIONS: We report the first GWAS for pediatric OSA in European Americans and African Americans. Our results provide novel insights to the genetic underpins of pediatric OSA.

Mutation burden analysis of six common mental disorders in African Americans by whole genome sequencing.

Mental disorders present a global health concern and have limited treatment options. In today's medical practice, medications such as antidepressants are prescribed not only for depression but also for conditions such as anxiety and attention deficit hyperactivity disorder (ADHD). Therefore, identifying gene targets for specific disorders is important and offers improved precision. In this study, we performed a genetic analysis of six common mental disorders-ADHD, anxiety, depression, delays in mental development, intellectual disabilities (IDs) and speech/language disorder-in the ethnic minority of African Americans (AAs) using whole genome sequencing (WGS). WGS data were generated from blood-derived DNA from 4178 AA individuals, including 1384 patients with the diagnosis of at least one mental disorder. Mutation burden analysis was applied based on rare and deleterious mutations in the AA population between cases and controls, and further analyzed in the context of patients with single mental disorder diagnosis. Certain genes uncovered demonstrated significant P-values in mutation burden analysis. In addition, exclusive recurrences in specific type of disorder were scanned through gene-drug interaction databases to assess for availability of potential medications. We uncovered 15 genes harboring deleterious mutations, including 3-Hydroxy-3-Methylglutaryl-CoA Reductase (HMGCR) and Uronyl 2-Sulfotransferase (UST) for ADHD; Farnesyltransferase, CAAX Box, Beta (FNTB) for anxiety; Xin Actin Binding Repeat Containing 2 (XIRP2), Natriuretic Peptide C (NPPC), Serine/Threonine Kinase 33 (STK33), Pannexin 1 (PANX1) and Neurotensin (NTS) for depression; RUNX Family Transcription Factor 3 (RUNX3), Tachykinin Receptor 1 (TACR1) and NADH:Ubiquinone Oxidoreductase Core Subunit S7 (NDUFS7) for delays in mental development; Hepsin (HPN) for ID and Collagen Type VI Alpha 3 Chain (COL6A3), Damage Specific DNA Binding Protein 1 (DDB1) and NADH:Ubiquinone Oxidoreductase Subunit A11 (NDUFA11) for speech/language disorder. Taken together, we have established critical insights into the development of new precision medicine approaches for mental disorders in AAs.

An electronic health record (EHR) phenotype algorithm to identify patients with attention deficit hyperactivity disorders (ADHD) and psychiatric comorbidities.

BACKGROUND: In over half of pediatric cases, ADHD presents with comorbidities, and often, it is unclear whether the symptoms causing impairment are due to the comorbidity or the underlying ADHD. Comorbid conditions increase the likelihood for a more severe and persistent course and complicate treatment decisions. Therefore, it is highly important to establish an algorithm that identifies ADHD and comorbidities in order to improve research on ADHD using biorepository and other electronic record data. METHODS: It is feasible to accurately distinguish between ADHD in isolation from ADHD with comorbidities using an electronic algorithm designed to include other psychiatric disorders. We sought to develop an EHR phenotype algorithm to discriminate cases with ADHD in isolation from cases with ADHD with comorbidities more effectively for efficient future searches in large biorepositories. We developed a multi-source algorithm allowing for a more complete view of the patient's EHR, leveraging the biobank of the Center for Applied Genomics (CAG) at Children's Hospital of Philadelphia (CHOP). We mined EHRs from 2009 to 2016 using International Statistical Classification of Diseases and Related Health Problems (ICD) codes, medication history and keywords specific to ADHD, and comorbid psychiatric disorders to facilitate genotype-phenotype correlation efforts. Chart abstractions and behavioral surveys added evidence in support of the psychiatric diagnoses. Most notably, the algorithm did not exclude other psychiatric disorders, as is the case in many previous algorithms. Controls lacked psychiatric and other neurological disorders. Participants enrolled in various CAG studies at CHOP and completed a broad informed consent, including consent for prospective analyses of EHRs. We created and validated an EHR-based algorithm to classify ADHD and comorbid psychiatric status in a pediatric healthcare network to be used in future genetic analyses and discovery-based studies. RESULTS: In this retrospective case-control study that included data from 51,293 subjects, 5840 ADHD cases were discovered of which 46.1% had ADHD alone and 53.9% had ADHD with psychiatric comorbidities. Our primary study outcome was to examine whether the algorithm could identify and distinguish ADHD exclusive cases from ADHD comorbid cases. The results indicate ICD codes coupled with medication searches revealed the most cases. We discovered ADHD-related keywords did not increase yield. However, we found including ADHD-specific medications increased our number of cases by 21%. Positive predictive values (PPVs) were 95% for ADHD cases and 93% for controls. CONCLUSION: We established a new algorithm and demonstrated the feasibility of the electronic algorithm approach to accurately diagnose ADHD and comorbid conditions, verifying the efficiency of our large biorepository for further genetic discovery-based analyses. TRIAL REGISTRATION: ClinicalTrials.gov, NCT02286817 . First posted on 10 November 2014. CLINICALTRIALS: gov, NCT02777931 . First posted on 19 May 2016. CLINICALTRIALS: gov, NCT03006367 . First posted on 30 December 2016. CLINICALTRIALS: gov, NCT02895906 . First posted on 12 September 2016.

Saudi Arabian CML patient with a novel four-way translocation at t(9;22;5;2)(q34;q11.2;p13;q44).

BACKGROUND: The vast majority of chronic myeloid leukemia (CML) patients have a single translocation t(9;22)(q34;q11), BCR/ABL1 fusion genes, which is regarded as the hallmark of CML. However, around 5 to 10% of CML patients exhibit the involvement of a third chromosome. In some very rare cases a fourth or even fifth chromosome can be involved with the t(9;22). METHODS: This case report is based on a 40-year-old Saudi Arabian male patient, diagnosed with CML in lymphoid blast crisis, and observed to have a four-way 46 XY, t(9;22;5;2)(q34;q11.2;p13;q44) translocation. The BCR/ABL1 fusion was identified by fluorescent in situ hybridization (FISH). Additionally, the BCR/ABL1 p210 mRNA fusion transcripts were identified by a molecular test. RESULTS: The clinical and prognostic impact of additional partner chromosomes to t(9;22) remains unknown. The CML patient with this novel four-way translocation t(9;22;5;2) progressed to blast crisis and was resistant to Tyrosine Kinase Inhibitor (TKI) therapy. Therefore, this case is more in alignment with the negative impact of additional partner chromosomes to the translocation at t(9;22). CONCLUSION: Here we report for the first time a novel four-way translocation at t(9;22;5;2)(q34;q11.2;p13;q44).

Burden of rare coding variants reveals genetic heterogeneity between obese and non-obese asthma patients in the African American population.

BACKGROUND: Asthma is a complex condition largely attributed to the interactions among genes and environments as a heterogeneous phenotype. Obesity is significantly associated with asthma development, and genetic studies on obese vs. non-obese asthma are warranted. METHODS: To investigate asthma in the minority African American (AA) population with or without obesity, we performed a whole genome sequencing (WGS) study on blood-derived DNA of 4289 AA individuals, included 2226 asthma patients (1364 with obesity and 862 without obesity) and 2006 controls without asthma. The burden analysis of functional rare coding variants was performed by comparing asthma vs. controls and by stratified analysis of obese vs. non-obese asthma, respectively. RESULTS: Among the top 66 genes with P < 0.01 in the asthma vs. control analysis, stratified analysis by obesity showed inverse correlation of natural logarithm (LN) of P value between obese and non-obese asthma (r = - 0.757, P = 1.90E-13). Five genes previously reported in a genome-wide association study (GWAS) on asthma, including TSLP, SLC9A4, PSMB8, IGSF5, and IKZF4 were demonstrated association in the asthma vs. control analysis. The associations of IKZF4 and IGSF5 are only associated with obese asthma; and the association of SLC9A4 is only observed in non-obese asthma. In addition, the association of RSPH3 (the gene is related to primary ciliary dyskinesia) is observed in non-obese asthma. CONCLUSIONS: These findings highlight genetic heterogeneity between obese and non-obese asthma in patients of AA ancestry.

Association of Pathogenic Variants in Hereditary Cancer Genes With Multiple Diseases.

Importance: Knowledge about the spectrum of diseases associated with hereditary cancer syndromes may improve disease diagnosis and management for patients and help to identify high-risk individuals. Objective: To identify phenotypes associated with hereditary cancer genes through a phenome-wide association study. Design, Setting, and Participants: This phenome-wide association study used health data from participants in 3 cohorts. The Electronic Medical Records and Genomics Sequencing (eMERGEseq) data set recruited predominantly healthy individuals from 10 US medical centers from July 16, 2016, through February 18, 2018, with a mean follow-up through electronic health records (EHRs) of 12.7 (7.4) years. The UK Biobank (UKB) cohort recruited participants from March 15, 2006, through August 1, 2010, with a mean (SD) follow-up of 12.4 (1.0) years. The Hereditary Cancer Registry (HCR) recruited patients undergoing clinical genetic testing at Vanderbilt University Medical Center from May 1, 2012, through December 31, 2019, with a mean (SD) follow-up through EHRs of 8.8 (6.5) years. Exposures: Germline variants in 23 hereditary cancer genes. Pathogenic and likely pathogenic variants for each gene were aggregated for association analyses. Main Outcomes and Measures: Phenotypes in the eMERGEseq and HCR cohorts were derived from the linked EHRs. Phenotypes in UKB were from multiple sources of health-related data. Results: A total of 214 020 participants were identified, including 23 544 in eMERGEseq cohort (mean [SD] age, 47.8 [23.7] years; 12 611 women [53.6%]), 187 234 in the UKB cohort (mean [SD] age, 56.7 [8.1] years; 104 055 [55.6%] women), and 3242 in the HCR cohort (mean [SD] age, 52.5 [15.5] years; 2851 [87.9%] women). All 38 established gene-cancer associations were replicated, and 19 new associations were identified. These included the following 7 associations with neoplasms: CHEK2 with leukemia (odds ratio [OR], 3.81 [95% CI, 2.64-5.48]) and plasma cell neoplasms (OR, 3.12 [95% CI, 1.84-5.28]), ATM with gastric cancer (OR, 4.27 [95% CI, 2.35-7.44]) and pancreatic cancer (OR, 4.44 [95% CI, 2.66-7.40]), MUTYH (biallelic) with kidney cancer (OR, 32.28 [95% CI, 6.40-162.73]), MSH6 with bladder cancer (OR, 5.63 [95% CI, 2.75-11.49]), and APC with benign liver/intrahepatic bile duct tumors (OR, 52.01 [95% CI, 14.29-189.29]). The remaining 12 associations with nonneoplastic diseases included BRCA1/2 with ovarian cysts (OR, 3.15 [95% CI, 2.22-4.46] and 3.12 [95% CI, 2.36-4.12], respectively), MEN1 with acute pancreatitis (OR, 33.45 [95% CI, 9.25-121.02]), APC with gastritis and duodenitis (OR, 4.66 [95% CI, 2.61-8.33]), and PTEN with chronic gastritis (OR, 15.68 [95% CI, 6.01-40.92]). Conclusions and Relevance: The findings of this genetic association study analyzing the EHRs of 3 large cohorts suggest that these new phenotypes associated with hereditary cancer genes may facilitate early detection and better management of cancers. This study highlights the potential benefits of using EHR data in genomic medicine.

Identification of Novel Loci Shared by Juvenile Idiopathic Arthritis Subtypes Through Integrative Genetic Analysis.

OBJECTIVE: Juvenile idiopathic arthritis (JIA) is the most common chronic immune-mediated joint disease among children and encompasses a heterogeneous group of immune-mediated joint disorders classified into 7 subtypes according to clinical presentation. However, phenotype overlap and biologic evidence suggest a shared mechanistic basis between subtypes. This study was undertaken to systematically investigate shared genetic underpinnings of JIA subtypes. METHODS: We performed a heterogeneity-sensitive genome-wide association study encompassing a total of 1,245 JIA cases (classified into 7 subtypes) and 9,250 controls, followed by fine-mapping of candidate causal variants at each genome-wide significant locus, functional annotation, and pathway and network analysis. We further identified candidate drug targets and drug repurposing opportunities by in silico analyses. RESULTS: In addition to the major histocompatibility complex locus, we identified 15 genome-wide significant loci shared between at least 2 JIA subtypes, including 10 novel loci. Functional annotation indicated that candidate genes at these loci were expressed in diverse immune cell types. CONCLUSION: This study identified novel genetic loci shared by JIA subtypes. Our findings identified candidate mechanisms underlying JIA subtypes and candidate targets with drug repurposing opportunities for JIA treatment.

Application of deep learning algorithm on whole genome sequencing data uncovers structural variants associated with multiple mental disorders in African American patients.

Mental disorders present a global health concern, while the diagnosis of mental disorders can be challenging. The diagnosis is even harder for patients who have more than one type of mental disorder, especially for young toddlers who are not able to complete questionnaires or standardized rating scales for diagnosis. In the past decade, multiple genomic association signals have been reported for mental disorders, some of which present attractive drug targets. Concurrently, machine learning algorithms, especially deep learning algorithms, have been successful in the diagnosis and/or labeling of complex diseases, such as attention deficit hyperactivity disorder (ADHD) or cancer. In this study, we focused on eight common mental disorders, including ADHD, depression, anxiety, autism, intellectual disabilities, speech/language disorder, delays in developments, and oppositional defiant disorder in the ethnic minority of African Americans. Blood-derived whole genome sequencing data from 4179 individuals were generated, including 1384 patients with the diagnosis of at least one mental disorder. The burden of genomic variants in coding/non-coding regions was applied as feature vectors in the deep learning algorithm. Our model showed ~65% accuracy in differentiating patients from controls. Ability to label patients with multiple disorders was similarly successful, with a hamming loss score less than 0.3, while exact diagnostic matches are around 10%. Genes in genomic regions with the highest weights showed enrichment of biological pathways involved in immune responses, antigen/nucleic acid binding, chemokine signaling pathway, and G-protein receptor activities. A noticeable fact is that variants in non-coding regions (e.g., ncRNA, intronic, and intergenic) performed equally well as variants in coding regions; however, unlike coding region variants, variants in non-coding regions do not express genomic hotspots whereas they carry much more narrow standard deviations, indicating they probably serve as alternative markers.

Rare neurological manifestations in a Saudi Arabian patient with Ehlers-Danlos syndrome and a novel homozygous variant in the TNXB gene.

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.

Circulating LIGHT (TNFSF14) and Interleukin-18 Levels in Sepsis-Induced Multi-Organ Injuries.

The novel therapeutic target cytokine LIGHT (TNFSF14) was recently shown to play a major role in COVID19-induced acute respiratory distress syndrome (ARDS). This study aims to investigate the associations of plasma LIGHT and another potentially targetable cytokine, Interleukin-18 (IL-18), with ARDS, acute hypoxic respiratory failure (AHRF) or acute kidney injury (AKI), caused by non-COVID19 viral or bacterial sepsis. A cohort of 280 subjects diagnosed with sepsis, including 91 cases with sepsis triggered by viral infections, were investigated in this cohort study. Day 0 plasma LIGHT and IL-18, as well as 59 other biomarkers (cytokines, chemokines and acute-phase reactants) were measured by sensitive bead immunoassay and associated with symptom severity. We observed significantly increased LIGHT level in both bacterial sepsis patients (P=1.80E-05) and patients with sepsis from viral infections (P=1.78E-03). In bacterial sepsis, increased LIGHT level was associated with ARDS, AKI and higher Apache III scores, findings also supported by correlations of LIGHT with other biomarkers of organ failures. IL-18 levels were highly variable across individuals, and consistently correlated with Apache III scores, mortality, and AKI, in both bacterial and viral sepsis. There was no correlation between LIGHT and IL-18. For the first time, we demonstrate independent effects of LIGHT and IL-18 in septic organ failures. The association of plasma LIGHT with AHRF suggests that targeting the pathway warrants exploration, and ongoing trials may soon elucidate whether this is beneficial. Given the large variance of plasma IL-18 among septic subjects, targeting this pathway requires a precision application.

Rare variant analysis in eczema identifies exonic variants in DUSP1, NOTCH4 and SLC9A4.

Previous genome-wide association studies revealed multiple common variants involved in eczema but the role of rare variants remains to be elucidated. Here, we investigate the role of rare variants in eczema susceptibility. We meta-analyze 21 study populations including 20,016 eczema cases and 380,433 controls. Rare variants are imputed with high accuracy using large population-based reference panels. We identify rare exonic variants in DUSP1, NOTCH4, and SLC9A4 to be associated with eczema. In DUSP1 and NOTCH4 missense variants are predicted to impact conserved functional domains. In addition, five novel common variants at SATB1-AS1/KCNH8, TRIB1/LINC00861, ZBTB1, TBX21/OSBPL7, and CSF2RB are discovered. While genes prioritized based on rare variants are significantly up-regulated in the skin, common variants point to immune cell function. Over 20% of the single nucleotide variant-based heritability is attributable to rare and low-frequency variants. The identified rare/low-frequency variants located in functional protein domains point to promising targets for novel therapeutic approaches to eczema.

HIF-1α Pulmonary Phenotype Wide Association Study Unveils a Link to Inflammatory Airway Conditions.

Despite experimental data linking HIF-1α dysfunction to inflammatory airway conditions, the effect of single nucleotide polymorphisms within the HIF1A gene on these conditions remains poorly understood. In the current study, we complete a phenotype wide association study to assess the link between SNPs with known disease associations and respiratory phenotypes. We report two SNPs of the HIF1A gene, the intronic rs79865957 and the missense rs41508050. In these positions the A and the T allele are significantly associated with allergic rhinitis and acute bronchitis and bronchiolitis, respectively. These findings further support the role of HIF-1α in inflammatory pulmonary conditions and may serve as a basis to refine our understanding of other HIF-1α associated phenotypes.

Insights into non-autoimmune type 1 diabetes with 13 novel loci in low polygenic risk score patients.

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 < 0.90). For the 13 novel loci, 9 regions have been reported of association with obesity related traits by previous GWA studies. Three loci encoding long intergenic non-protein coding RNAs (lncRNA), and 2 loci involved in N-linked glycosylation are also highlighted in this study.