Juvenile idiopathic arthritis (JIA) is the most common chronic rheumatic disease among children, the etiology of which involves a strong genetic component, but much of the underlying genetic determinants still remain unknown. Our aim was to identify novel genetic variants that predispose to JIA.

Absence of the anterior (ACL) or posterior cruciate ligament (PCL) are rare congenital malformations that result in knee joint instability, with a prevalence of 1.7 per 100,000 live births and can be associated with other lower-limb abnormalities such as ACL agnesia and absence of the menisci of the knee. While a few cases of absence of ACL/PCL are reported in the literature, a number of large familial case series of related conditions such as ACL agnesia suggest a potential underlying monogenic etiology. We performed whole exome sequencing of a family with two individuals affected by ACL/PCL.

Through an international multi-center collaboration, 13 individuals from nine unrelated families and affected by likely pathogenic biallelic variants in TBC1-domain-containing kinase (TBCK) were identified through whole-exome sequencing. All affected individuals were found to share a core phenotype of intellectual disability and hypotonia, and many had seizures and showed brain atrophy and white-matter changes on neuroimaging. Minor non-specific facial dysmorphism was also noted in some individuals, including multiple older children who developed coarse features similar to those of storage disorders. TBCK has been shown to regulate the mammalian target of rapamycin (mTOR) signaling pathway, which is also stimulated by exogenous leucine supplementation. TBCK was absent in cells from affected individuals, and decreased phosphorylation of phospho-ribosomal protein S6 was also observed, a finding suggestive of downregulation of mTOR signaling. Lastly, we demonstrated that activation of the mTOR pathway in response to L-leucine supplementation was retained, suggesting a possible avenue for directed therapies for this condition.

The standard expression quantitative trait loci (eQTL) detects polymorphisms associated with gene expression without revealing causality. We introduce a coupled Bayesian regression approach--eQTeL, which leverages epigenetic data to estimate regulatory and gene interaction potential, and identifies combination of regulatory single-nucleotide polymorphisms (SNPs) that explain the gene expression variance. On human heart data, eQTeL not only explains a significantly greater proportion of expression variance but also predicts gene expression more accurately than other methods. Based on realistic simulated data, we demonstrate that eQTeL accurately detects causal regulatory SNPs, including those with small effect sizes. Using various functional data, we show that SNPs detected by eQTeL are enriched for allele-specific protein binding and histone modifications, which potentially disrupt binding of core cardiac transcription factors and are spatially proximal to their target. eQTeL SNPs capture a substantial proportion of genetic determinants of expression variance and we estimate that 58% of these SNPs are putatively causal.

In addition to HLA genetic incompatibility, non-HLA difference between donor and recipients of transplantation leading to allograft rejection are now becoming evident. We aimed to create a unique genome-wide platform to facilitate genomic research studies in transplant-related studies. We designed a genome-wide genotyping tool based on the most recent human genomic reference datasets, and included customization for known and potentially relevant metabolic and pharmacological loci relevant to transplantation.

Genetic association studies have identified 21 loci associated with atopic dermatitis risk predominantly in populations of European ancestry. To identify further susceptibility loci for this common, complex skin disease, we performed a meta-analysis of >15 million genetic variants in 21,399 cases and 95,464 controls from populations of European, African, Japanese and Latino ancestry, followed by replication in 32,059 cases and 228,628 controls from 18 studies. We identified ten new risk loci, bringing the total number of known atopic dermatitis risk loci to 31 (with new secondary signals at four of these loci). Notably, the new loci include candidate genes with roles in the regulation of innate host defenses and T cell function, underscoring the important contribution of (auto)immune mechanisms to atopic dermatitis pathogenesis.

Eosinophilic esophagitis (EoE) is an allergic disorder characterized by infiltration of the oesophagus with eosinophils. We had previously reported association of the TSLP/WDR36 locus with EoE. Here we report genome-wide significant associations at four additional loci; c11orf30 and STAT6, which have been previously associated with both atopic and autoimmune diseases, and two EoE-specific loci, ANKRD27 that regulates the trafficking of melanogenic enzymes to epidermal melanocytes and CAPN14, that encodes a calpain whose expression is highly enriched in the oesophagus. The identification of five EoE loci, not only expands our aetiological understanding of the disease but may also represent new therapeutic targets to treat the most debilitating aspect of EoE, oesophageal inflammation and remodelling.

Calcium, magnesium, potassium, sodium, chloride and phosphorus are the major dietary minerals involved in various biological functions and are commonly measured in the blood serum. Sufficient mineral intake is especially important for children due to their rapid growth. Currently, the genetic mechanisms influencing serum mineral levels are poorly understood, especially for children. We carried out a genome-wide association (GWA) study on 5,602 European-American children and 4,706 African-American children who had mineral measures available in their electronic medical records (EMR). While no locus met the criteria for genome-wide significant association, our results demonstrated a nominal association of total serum calcium levels with a missense variant in the calcium -sensing receptor (CASR) gene on 3q13 (rs1801725, P = 1.96 × 10(-3)) in the African-American pediatric cohort, a locus previously reported in Caucasians. We also confirmed the association result in our pediatric European-American cohort (P = 1.38 × 10(-4)). We further replicated two other loci associated with serum calcium levels in the European-American cohort (rs780094, GCKR, P = 4.26 × 10(-3); rs10491003, GATA3, P = 0.02). In addition, we replicated a previously reported locus on 1q21, demonstrating association of serum magnesium levels with MUC1 (rs4072037, P = 2.04 × 10(-6)). Moreover, in an extended gene-based association analysis we uncovered evidence for association of calcium levels with the previously reported gene locus DGKD in both European-American children and African-American children. Taken together, our results support a role for CASR and DGKD mediated calcium regulation in both African-American and European-American children, and corroborate the association of calcium levels with GCKR and GATA3, and the association of magnesium levels with MUC1 in the European-American children.

Congenital anomalies of the kidney and urinary tract (CAKUT) are a major cause of pediatric kidney failure. We performed a genome-wide analysis of copy number variants (CNVs) in 2,824 cases and 21,498 controls. Affected individuals carried a significant burden of rare exonic (that is, affecting coding regions) CNVs and were enriched for known genomic disorders (GD). Kidney anomaly (KA) cases were most enriched for exonic CNVs, encompassing GD-CNVs and novel deletions; obstructive uropathy (OU) had a lower CNV burden and an intermediate prevalence of GD-CNVs; and vesicoureteral reflux (VUR) had the fewest GD-CNVs but was enriched for novel exonic CNVs, particularly duplications. Six loci (1q21, 4p16.1-p16.3, 16p11.2, 16p13.11, 17q12 and 22q11.2) accounted for 65% of patients with GD-CNVs. Deletions at 17q12, 4p16.1-p16.3 and 22q11.2 were specific for KA; the 16p11.2 locus showed extensive pleiotropy. Using a multidisciplinary approach, we identified TBX6 as a driver for the CAKUT subphenotypes in the 16p11.2 microdeletion syndrome.

The link between cardiovascular diseases and neurological disorders has been widely observed in the aging population. Disease prevention and treatment rely on understanding the potential genetic nexus of multiple diseases in these categories. In this study, we were interested in detecting pleiotropy, or the phenomenon in which a genetic variant influences more than one phenotype. Marker-phenotype association approaches can be grouped into univariate, bivariate, and multivariate categories based on the number of phenotypes considered at one time. Here we applied one statistical method per category followed by an eQTL colocalization analysis to identify potential pleiotropic variants that contribute to the link between cardiovascular and neurological diseases. We performed our analyses on ~530,000 common SNPs coupled with 65 electronic health record (EHR)-based phenotypes in 43,870 unrelated European adults from the Electronic Medical Records and Genomics (eMERGE) network. There were 31 variants identified by all three methods that showed significant associations across late onset cardiac- and neurologic- diseases. We further investigated functional implications of gene expression on the detected "lead SNPs" via colocalization analysis, providing a deeper understanding of the discovered associations. In summary, we present the framework and landscape for detecting potential pleiotropy using univariate, bivariate, multivariate, and colocalization methods. Further exploration of these potentially pleiotropic genetic variants will work toward understanding disease causing mechanisms across cardiovascular and neurological diseases and may assist in considering disease prevention as well as drug repositioning in future research.

To date, most transcriptome studies of schizophrenia focus on the analysis of protein-coding genes. Long noncoding RNAs (lncRNAs) are emerging as key tissue-specific regulators of cellular and disease processes. The amygdala brain region has been implicated in the pathophysiology of schizophrenia. We performed unbiased whole transcriptome profiling of amygdala tissues from 22 schizophrenia patients and 24 non-psychiatric controls using RNA-seq. We reconstructed amygdala transcriptome and employed systems biology approaches to annotating the functional roles of lncRNAs. As a result, we identified 839 novel lncRNAs in amygdala. We found in amygdala lncRNAs are more subtype-specific than protein-coding genes. We identified functional modules associated with "synaptic transmission", "ribosome", and "immune responses" which were related to schizophrenia pathophysiology that involved lncRNAs. Integrative functional analyses associating individual lncRNAs with specific pathways and functions further show that amygdala lncRNAs are connected with all of these pathways. Our study presents the first systematic landscape of lncRNAs in amygdala tissue from schizophrenia cases.

The goal of this paper is to review recent research on copy number variations (CNVs) and their association with complex and rare diseases. In the latter part of this paper, we focus on how large biorepositories such as the electronic medical record and genomics (eMERGE) consortium may be best leveraged to systematically mine for potentially pathogenic CNVs, and we end with a discussion of how such variants might be reported back for inclusion in electronic medical records as part of medical history.

Rare copy-number variation (CNV) is an important source of risk for autism spectrum disorders (ASDs). We analyzed 2,446 ASD-affected families and confirmed an excess of genic deletions and duplications in affected versus control groups (1.41-fold, p = 1.0 × 10(-5)) and an increase in affected subjects carrying exonic pathogenic CNVs overlapping known loci associated with dominant or X-linked ASD and intellectual disability (odds ratio = 12.62, p = 2.7 × 10(-15), ∼3% of ASD subjects). Pathogenic CNVs, often showing variable expressivity, included rare de novo and inherited events at 36 loci, implicating ASD-associated genes (CHD2, HDAC4, and GDI1) previously linked to other neurodevelopmental disorders, as well as other genes such as SETD5, MIR137, and HDAC9. Consistent with hypothesized gender-specific modulators, females with ASD were more likely to have highly penetrant CNVs (p = 0.017) and were also overrepresented among subjects with fragile X syndrome protein targets (p = 0.02). Genes affected by de novo CNVs and/or loss-of-function single-nucleotide variants converged on networks related to neuronal signaling and development, synapse function, and chromatin regulation.

Estimates from large scale genome sequencing studies indicate that each human carries up to 20 genetic variants that are predicted to results in loss of function (LOF) of protein-coding genes. While some are known disease-causing variants or common, tolerated, LOFs in non-essential genes, the majority remain of unknown consequence. We explore the possibility of using imputed GWAS data from large biorepositories such as the electronic medical record and genomics (eMERGE) consortium to determine the effects of rare LOFs. Here, we show that two hypocholesterolemia-associated LOF mutations in the PCSK9 gene can be accurately imputed into large-scale GWAS datasets which raises the possibility of assessing LOFs through genomics-linked medical records.

Candidate gene studies have revealed limited genetic bases for opioid analgesic response variability. Genome-wide association studies facilitate impartial queries of common genetic variants, allowing identification of novel genetic contributions to drug effect. Illumina (Illumina Inc, San Diego, CA, USA) single nucleotide polymorphism (SNP) arrays were used to investigate SNP associations with total morphine requirement as a quantitative trait locus and with postoperative pain in a retrospective population of opioid-naïve children ages 4-18years who had undergone day surgery tonsillectomy and adenoidectomy. In an independent replication cohort, significant genome-wide association studies-identified SNPs were assayed using TaqMan probes. Among 617 comprehensively phenotyped children, the 277 subjects of European Caucasian (EC) ancestry demonstrated nominal association between morphine dose and a series of novel SNPs (top rs795484, P=1.01 × 10(-6) and rs1277441, P=2.77 × 10(-6)) at the TAOK3 locus. Age, body mass index, and physical status were included covariates. Morphine requirement averaged 132.4 μg/kg (SD 40.9). Each minor allele at rs795484 (guanine [G]>adenine [A]) contributed +17.6 μg/kg (95% confidence interval [CI] 10.7-24.4) to dose. Effect direction and magnitude were replicated in an independent cohort of 75 EC children (P<0.05). No association with morphine dose was detected in African Americans (AA) (n=241). Postoperative pain scores ≥ 7/10 were associated with rs795484 (G>A) in the EC cohort (odds ratio 2.35, 95% CI 1.56-3.52, P<0.00005) and this association replicated in AA children (odds ratio 1.76, 95% CI 1.14-2.71, P<0.01). Variants in TAOK3 encoding the serine/threonine-protein kinase, TAO3, are associated with increased morphine requirement in children of EC ancestry and with increased acute postoperative pain in both EC and AA subjects.

Social communication difficulties represent an autistic trait that is highly heritable and persistent during the course of development. However, little is known about the underlying genetic architecture of this phenotype.

Statin treatment and variants in the gene encoding HMG-CoA reductase are associated with reductions in both the concentration of LDL cholesterol and the risk of coronary heart disease, but also with modest hyperglycaemia, increased bodyweight, and modestly increased risk of type 2 diabetes, which in no way offsets their substantial benefits. We sought to investigate the associations of LDL cholesterol-lowering PCSK9 variants with type 2 diabetes and related biomarkers to gauge the likely effects of PCSK9 inhibitors on diabetes risk.

Hepatoblastoma (HBL), the most common childhood liver cancer is cured with surgical resection after chemotherapy or with liver transplantation if local invasion and multifocality preclude resection. However, variable survival rates of 60-80% and debilitating chemotherapy sequelae argue for more informed treatment selection, which is not possible by grading the Wnt-β-catenin over activity present in most HBL tumors. A hypothesis-generating whole transcriptome analysis shows that HBL tumors removed at transplantation are enriched most for cancer signaling pathways which depend predominantly on epidermal growth factor (EGF) signaling, and to a lesser extent, on aberrant Wnt-β-catenin signaling. We therefore evaluated whether EGFR, ASAP1, ERBB2 and ERBB4, which signal downstream after ligation of EGF, and which show aberrant expression in several other invasive cancers, would also predict HBL tumor invasiveness. Immunohistochemistry of HBL tumors (n = 60), which are histologically heterogeneous, shows that compared with well-differentiated fetal cells, less differentiated embryonal and undifferentiated small cells (SCU) progressively lose EGFR and ASAP1 expression. This trend is exaggerated in unresectable, locally invasive or metastatic tumors, in which embryonal tumor cells are EGFR-negative, while SCU cells are EGFR-negative and ASAP1-negative. Loss of EGFR-ASAP1 signaling characterizes undifferentiated and invasive HBL. EGFR-expressing HBL tumors present novel therapeutic targeting opportunities.

The endocrine system is crucial for maintaining whole-body homeostasis. Little is known regarding endocrine hormones secreted by the heart other than atrial/brain natriuretic peptides discovered over 30 years ago. Here, we identify growth differentiation factor 15 (GDF15) as a heart-derived hormone that regulates body growth. We show that pediatric heart disease induces GDF15 synthesis and secretion by cardiomyocytes. Circulating GDF15 in turn acts on the liver to inhibit growth hormone (GH) signaling and body growth. We demonstrate that blocking cardiomyocyte production of GDF15 normalizes circulating GDF15 level and restores liver GH signaling, establishing GDF15 as a bona fide heart-derived hormone that regulates pediatric body growth. Importantly, plasma GDF15 is further increased in children with concomitant heart disease and failure to thrive (FTT). Together these studies reveal a new endocrine mechanism by which the heart coordinates cardiac function and body growth. Our results also provide a potential mechanism for the well-established clinical observation that children with heart diseases often develop FTT.

Neuroblastoma is a cancer of the developing sympathetic nervous system that most commonly presents in young children and accounts for approximately 12% of pediatric oncology deaths. Here, we report on a genome-wide association study (GWAS) in a discovery cohort or 2,101 cases and 4,202 controls of European ancestry. We identify two new association signals at 3q25 and 4p16 that replicated robustly in multiple independent cohorts comprising 1,163 cases and 4,396 controls (3q25: rs6441201 combined P = 1.2x10-11, Odds Ratio 1.23, 95% CI:1.16-1.31; 4p16: rs3796727 combined P = 1.26x10-12, Odds Ratio 1.30, 95% CI: 1.21-1.40). The 4p16 signal maps within the carboxypeptidase Z (CPZ) gene. The 3q25 signal resides within the arginine/serine-rich coiled-coil 1 (RSRC1) gene and upstream of the myeloid leukemia factor 1 (MLF1) gene. Increased expression of MLF1 was observed in neuroblastoma cells homozygous for the rs6441201 risk allele (P = 0.02), and significant growth inhibition was observed upon depletion of MLF1 (P < 0.0001) in neuroblastoma cells. Taken together, we show that common DNA variants within CPZ at 4p16 and upstream of MLF1 at 3q25 influence neuroblastoma susceptibility and MLF1 likely plays an important role in neuroblastoma tumorigenesis.