Genetic regulation of human brain proteome reveals proteins implicated in psychiatric disorders.

Psychiatric disorders are highly heritable yet polygenic, potentially involving hundreds of risk genes. Genome-wide association studies have identified hundreds of genomic susceptibility loci with susceptibility to psychiatric disorders; however, the contribution of these loci to the underlying psychopathology and etiology remains elusive. Here we generated deep human brain proteomics data by quantifying 11,608 proteins across 268 subjects using 11-plex tandem mass tag coupled with two-dimensional liquid chromatography-tandem mass spectrometry. Our analysis revealed 788 cis-acting protein quantitative trait loci associated with the expression of 883 proteins at a genome-wide false discovery rate <5%. In contrast to expression at the transcript level and complex diseases that are found to be mainly influenced by noncoding variants, we found protein expression level tends to be regulated by non-synonymous variants. We also provided evidence of 76 shared regulatory signals between gene expression and protein abundance. Mediation analysis revealed that for most (88%) of the colocalized genes, the expression levels of their corresponding proteins are regulated by cis-pQTLs via gene transcription. Using summary data-based Mendelian randomization analysis, we identified 4 proteins and 19 genes that are causally associated with schizophrenia. We further integrated multiple omics data with network analysis to prioritize candidate genes for schizophrenia risk loci. Collectively, our findings underscore the potential of proteome-wide linkage analysis in gaining mechanistic insights into the pathogenesis of psychiatric disorders.

Direct reprogramming induces vascular regeneration post muscle ischemic injury.

Reprogramming non-cardiomyocytes (non-CMs) into cardiomyocyte (CM)-like cells is a promising strategy for cardiac regeneration in conditions such as ischemic heart disease. Here, we used a modified mRNA (modRNA) gene delivery platform to deliver a cocktail, termed 7G-modRNA, of four cardiac-reprogramming genes-Gata4 (G), Mef2c (M), Tbx5 (T), and Hand2 (H)-together with three reprogramming-helper genes-dominant-negative (DN)-TGFß, DN-Wnt8a, and acid ceramidase (AC)-to induce CM-like cells. We showed that 7G-modRNA reprogrammed 57% of CM-like cells in vitro. Through a lineage-tracing model, we determined that delivering the 7G-modRNA cocktail at the time of myocardial infarction reprogrammed ∼25% of CM-like cells in the scar area and significantly improved cardiac function, scar size, long-term survival, and capillary density. Mechanistically, we determined that while 7G-modRNA cannot create de novo beating CMs in vitro or in vivo, it can significantly upregulate pro-angiogenic mesenchymal stromal cells markers and transcription factors. We also demonstrated that our 7G-modRNA cocktail leads to neovascularization in ischemic-limb injury, indicating CM-like cells importance in other organs besides the heart. modRNA is currently being used around the globe for vaccination against COVID-19, and this study proves this is a safe, highly efficient gene delivery approach with therapeutic potential to treat ischemic diseases.

Landscape of Conditional eQTL in Dorsolateral Prefrontal Cortex and Co-localization with Schizophrenia GWAS.

Causal genes and variants within genome-wide association study (GWAS) loci can be identified by integrating GWAS statistics with expression quantitative trait loci (eQTL) and determining which variants underlie both GWAS and eQTL signals. Most analyses, however, consider only the marginal eQTL signal, rather than dissect this signal into multiple conditionally independent signals for each gene. Here we show that analyzing conditional eQTL signatures, which could be important under specific cellular or temporal contexts, leads to improved fine mapping of GWAS associations. Using genotypes and gene expression levels from post-mortem human brain samples (n = 467) reported by the CommonMind Consortium (CMC), we find that conditional eQTL are widespread; 63% of genes with primary eQTL also have conditional eQTL. In addition, genomic features associated with conditional eQTL are consistent with context-specific (e.g., tissue-, cell type-, or developmental time point-specific) regulation of gene expression. Integrating the 2014 Psychiatric Genomics Consortium schizophrenia (SCZ) GWAS and CMC primary and conditional eQTL data reveals 40 loci with strong evidence for co-localization (posterior probability > 0.8), including six loci with co-localization of conditional eQTL. Our co-localization analyses support previously reported genes, identify novel genes associated with schizophrenia risk, and provide specific hypotheses for their functional follow-up.

HLA-DRB1*11 and variants of the MHC class II locus are strong risk factors for systemic juvenile idiopathic arthritis.

Systemic juvenile idiopathic arthritis (sJIA) is an often severe, potentially life-threatening childhood inflammatory disease, the pathophysiology of which is poorly understood. To determine whether genetic variation within the MHC locus on chromosome 6 influences sJIA susceptibility, we performed an association study of 982 children with sJIA and 8,010 healthy control subjects from nine countries. Using meta-analysis of directly observed and imputed SNP genotypes and imputed classic HLA types, we identified the MHC locus as a bona fide susceptibility locus with effects on sJIA risk that transcended geographically defined strata. The strongest sJIA-associated SNP, rs151043342 [P = 2.8 × 10(-17), odds ratio (OR) 2.6 (2.1, 3.3)], was part of a cluster of 482 sJIA-associated SNPs that spanned a 400-kb region and included the class II HLA region. Conditional analysis controlling for the effect of rs151043342 found that rs12722051 independently influenced sJIA risk [P = 1.0 × 10(-5), OR 0.7 (0.6, 0.8)]. Meta-analysis of imputed classic HLA-type associations in six study populations of Western European ancestry revealed that HLA-DRB1*11 and its defining amino acid residue, glutamate 58, were strongly associated with sJIA [P = 2.7 × 10(-16), OR 2.3 (1.9, 2.8)], as was the HLA-DRB1*11-HLA-DQA1*05-HLA-DQB1*03 haplotype [6.4 × 10(-17), OR 2.3 (1.9, 2.9)]. By examining the MHC locus in the largest collection of sJIA patients assembled to date, this study solidifies the relationship between the class II HLA region and sJIA, implicating adaptive immune molecules in the pathogenesis of sJIA.

Convergence of genes and cellular pathways dysregulated in autism spectrum disorders.

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.

Genetic architecture distinguishes systemic juvenile idiopathic arthritis from other forms of juvenile idiopathic arthritis: clinical and therapeutic implications.

OBJECTIVES: Juvenile idiopathic arthritis (JIA) is a heterogeneous group of conditions unified by the presence of chronic childhood arthritis without an identifiable cause. Systemic JIA (sJIA) is a rare form of JIA characterised by systemic inflammation. sJIA is distinguished from other forms of JIA by unique clinical features and treatment responses that are similar to autoinflammatory diseases. However, approximately half of children with sJIA develop destructive, long-standing arthritis that appears similar to other forms of JIA. Using genomic approaches, we sought to gain novel insights into the pathophysiology of sJIA and its relationship with other forms of JIA. METHODS: We performed a genome-wide association study of 770 children with sJIA collected in nine countries by the International Childhood Arthritis Genetics Consortium. Single nucleotide polymorphisms were tested for association with sJIA. Weighted genetic risk scores were used to compare the genetic architecture of sJIA with other JIA subtypes. RESULTS: The major histocompatibility complex locus and a locus on chromosome 1 each showed association with sJIA exceeding the threshold for genome-wide significance, while 23 other novel loci were suggestive of association with sJIA. Using a combination of genetic and statistical approaches, we found no evidence of shared genetic architecture between sJIA and other common JIA subtypes. CONCLUSIONS: The lack of shared genetic risk factors between sJIA and other JIA subtypes supports the hypothesis that sJIA is a unique disease process and argues for a different classification framework. Research to improve sJIA therapy should target its unique genetics and specific pathophysiological pathways.

De novo unbalanced translocation (4p duplication/8p deletion) in a patient with autism, OCD, and overgrowth syndrome.

Chromosomal abnormalities, such as unbalanced translocations and copy number variants (CNVs), are found in autism spectrum disorders (ASDs) [Sanders et al. (2011) Neuron 70: 863-885]. Many chromosomal abnormalities, including sub microscopic genomic deletions and duplications, are missed by G-banded karyotyping or Fragile X screening alone and are picked up by chromosomal microarrays [Shen et al. (2010) Pediatrics 125: e727-735]. Translocations involving chromosomes 4 and 8 are possibly the second most frequent translocation in humans and are often undetected in routine cytogenetics [Giglio et al. (2002) Circulation 102: 432-437]. Deletions of 4p16 have been associated with Wolf-Hirschhorn syndrome while 4p16 duplications have been associated with an overgrowth syndrome and mild to moderate mental retardation [Partington et al. (1997) Journal of Medical Genetics 34: 719-728]. The 8p23.3 region contains the autism candidate gene DLGAP2, which can contribute to autism when disrupted [Marshall et al. (2008) The American Journal of Human Genetics 82: 477-488] . There has been a case report of a family with autism spectrum disorder (ASD), prominent obsessional behavior, and overgrowth in patients with der (8) t (4;8) p (16;23) [Partington et al. (1997)]. This is an independent report of a male patient with autism, obsessive compulsive disorder (OCD), attention-deficit hyperactivity disorder (ADHD), and an overgrowth syndrome, whose de novo unbalanced translocation der (8) t (4;8) p (16.1→ter; 23.1→ter) was initially missed by routine cytogenetics but detected with SNP microarray, allowing higher resolution of translocation breakpoints.

Identification of novel genetic causes of Rett syndrome-like phenotypes.

BACKGROUND: The aim of this work was to identify new genetic causes of Rett-like phenotypes using array comparative genomic hybridisation and a whole exome sequencing approach. METHODS AND RESULTS: We studied a cohort of 19 Portuguese patients (16 girls, 3 boys) with a clinical presentation significantly overlapping Rett syndrome (RTT). Genetic analysis included filtering of the single nucleotide variants and indels with preference for de novo, homozygous/compound heterozygous, or maternally inherited X linked variants. Examination by MRI and muscle biopsies was also performed. Pathogenic genomic imbalances were found in two patients (10.5%): an 18q21.2 deletion encompassing four exons of the TCF4 gene and a mosaic UPD of chromosome 3. Variants in genes previously implicated in neurodevelopmental disorders (NDD) were identified in six patients (32%): de novo variants in EEF1A2, STXBP1 and ZNF238 were found in three patients, maternally inherited X linked variants in SLC35A2, ZFX and SHROOM4 were detected in two male patients and one homozygous variant in EIF2B2 was detected in one patient. Variants were also detected in five novel NDD candidate genes (26%): we identified de novo variants in the RHOBTB2, SMARCA1 and GABBR2 genes; a homozygous variant in EIF4G1; compound heterozygous variant in HTT. CONCLUSIONS: Network analysis reveals that these genes interact by means of protein interactions with each other and with the known RTT genes. These findings expand the phenotypical spectrum of previously known NDD genes to encompass RTT-like clinical presentations and identify new candidate genes for RTT-like phenotypes.

Meta-analysis of SHANK Mutations in Autism Spectrum Disorders: a gradient of severity in cognitive impairments.

SHANK genes code for scaffold proteins located at the post-synaptic density of glutamatergic synapses. In neurons, SHANK2 and SHANK3 have a positive effect on the induction and maturation of dendritic spines, whereas SHANK1 induces the enlargement of spine heads. Mutations in SHANK genes have been associated with autism spectrum disorders (ASD), but their prevalence and clinical relevance remain to be determined. Here, we performed a new screen and a meta-analysis of SHANK copy-number and coding-sequence variants in ASD. Copy-number variants were analyzed in 5,657 patients and 19,163 controls, coding-sequence variants were ascertained in 760 to 2,147 patients and 492 to 1,090 controls (depending on the gene), and, individuals carrying de novo or truncating SHANK mutations underwent an extensive clinical investigation. Copy-number variants and truncating mutations in SHANK genes were present in ∼1% of patients with ASD: mutations in SHANK1 were rare (0.04%) and present in males with normal IQ and autism; mutations in SHANK2 were present in 0.17% of patients with ASD and mild intellectual disability; mutations in SHANK3 were present in 0.69% of patients with ASD and up to 2.12% of the cases with moderate to profound intellectual disability. In summary, mutations of the SHANK genes were detected in the whole spectrum of autism with a gradient of severity in cognitive impairment. Given the rare frequency of SHANK1 and SHANK2 deleterious mutations, the clinical relevance of these genes remains to be ascertained. In contrast, the frequency and the penetrance of SHANK3 mutations in individuals with ASD and intellectual disability-more than 1 in 50-warrant its consideration for mutation screening in clinical practice.

Origins and functional impact of copy number variation in the human genome.

Structural variations of DNA greater than 1 kilobase in size account for most bases that vary among human genomes, but are still relatively under-ascertained. Here we use tiling oligonucleotide microarrays, comprising 42 million probes, to generate a comprehensive map of 11,700 copy number variations (CNVs) greater than 443 base pairs, of which most (8,599) have been validated independently. For 4,978 of these CNVs, we generated reference genotypes from 450 individuals of European, African or East Asian ancestry. The predominant mutational mechanisms differ among CNV size classes. Retrotransposition has duplicated and inserted some coding and non-coding DNA segments randomly around the genome. Furthermore, by correlation with known trait-associated single nucleotide polymorphisms (SNPs), we identified 30 loci with CNVs that are candidates for influencing disease susceptibility. Despite this, having assessed the completeness of our map and the patterns of linkage disequilibrium between CNVs and SNPs, we conclude that, for complex traits, the heritability void left by genome-wide association studies will not be accounted for by common CNVs.

Functional impact of global rare copy number variation in autism spectrum disorders.

The autism spectrum disorders (ASDs) are a group of conditions characterized by impairments in reciprocal social interaction and communication, and the presence of restricted and repetitive behaviours. Individuals with an ASD vary greatly in cognitive development, which can range from above average to intellectual disability. Although ASDs are known to be highly heritable ( approximately 90%), the underlying genetic determinants are still largely unknown. Here we analysed the genome-wide characteristics of rare (<1% frequency) copy number variation in ASD using dense genotyping arrays. When comparing 996 ASD individuals of European ancestry to 1,287 matched controls, cases were found to carry a higher global burden of rare, genic copy number variants (CNVs) (1.19 fold, P = 0.012), especially so for loci previously implicated in either ASD and/or intellectual disability (1.69 fold, P = 3.4 x 10(-4)). Among the CNVs there were numerous de novo and inherited events, sometimes in combination in a given family, implicating many novel ASD genes such as SHANK2, SYNGAP1, DLGAP2 and the X-linked DDX53-PTCHD1 locus. We also discovered an enrichment of CNVs disrupting functional gene sets involved in cellular proliferation, projection and motility, and GTPase/Ras signalling. Our results reveal many new genetic and functional targets in ASD that may lead to final connected pathways.

Bone marrow failure and developmental delay caused by mutations in poly(A)-specific ribonuclease (PARN).

BACKGROUND: Deadenylation regulates RNA function and fate. Poly(A)-specific ribonuclease (PARN) is a deadenylase that processes mRNAs and non-coding RNA. Little is known about the biological significance of germline mutations in PARN. METHODS: We identified mutations in PARN in patients with haematological and neurological manifestations. Genomic, biochemical and knockdown experiments in human marrow cells and in zebrafish have been performed to clarify the role of PARN in the human disease. RESULTS: We identified large monoallelic deletions in PARN in four patients with developmental delay or mental illness. One patient in particular had a severe neurological phenotype, central hypomyelination and bone marrow failure. This patient had an additional missense mutation on the non-deleted allele and severely reduced PARN protein and deadenylation activity. Cells from this patient had impaired oligoadenylation of specific H/ACA box small nucleolar RNAs. Importantly, PARN-deficient patient cells manifested short telomeres and an aberrant ribosome profile similar to those described in some variants of dyskeratosis congenita. Knocking down PARN in human marrow cells and zebrafish impaired haematopoiesis, providing further evidence for a causal link with the human disease. CONCLUSIONS: Large monoallelic mutations of PARN can cause developmental/mental illness. Biallelic PARN mutations cause severe bone marrow failure and central hypomyelination.

Network topologies and convergent aetiologies arising from deletions and duplications observed in individuals with autism.

Autism Spectrum Disorders (ASD) are highly heritable and characterised by impairments in social interaction and communication, and restricted and repetitive behaviours. Considering four sets of de novo copy number variants (CNVs) identified in 181 individuals with autism and exploiting mouse functional genomics and known protein-protein interactions, we identified a large and significantly interconnected interaction network. This network contains 187 genes affected by CNVs drawn from 45% of the patients we considered and 22 genes previously implicated in ASD, of which 192 form a single interconnected cluster. On average, those patients with copy number changed genes from this network possess changes in 3 network genes, suggesting that epistasis mediated through the network is extensive. Correspondingly, genes that are highly connected within the network, and thus whose copy number change is predicted by the network to be more phenotypically consequential, are significantly enriched among patients that possess only a single ASD-associated network copy number changed gene (p = 0.002). Strikingly, deleted or disrupted genes from the network are significantly enriched in GO-annotated positive regulators (2.3-fold enrichment, corrected p = 2×10(-5)), whereas duplicated genes are significantly enriched in GO-annotated negative regulators (2.2-fold enrichment, corrected p = 0.005). The direction of copy change is highly informative in the context of the network, providing the means through which perturbations arising from distinct deletions or duplications can yield a common outcome. These findings reveal an extensive ASD-associated molecular network, whose topology indicates ASD-relevant mutational deleteriousness and that mechanistically details how convergent aetiologies can result extensively from CNVs affecting pathways causally implicated in ASD.

Microdeletions of ELP4 Are Associated with Language Impairment, Autism Spectrum Disorder, and Mental Retardation.

Copy-number variations (CNVs) are important in the aetiology of neurodevelopmental disorders and show broad phenotypic manifestations. We compared the presence of small CNVs disrupting the ELP4-PAX6 locus in 4,092 UK individuals with a range of neurodevelopmental conditions, clinically referred for array comparative genomic hybridization, with WTCCC controls (n = 4,783). The phenotypic analysis was then extended using the DECIPHER database. We followed up association using an autism patient cohort (n = 3,143) compared with six additional control groups (n = 6,469). In the clinical discovery series, we identified eight cases with ELP4 deletions, and one with a partial duplication of ELP4 and PAX6. These cases were referred for neurological phenotypes including language impairment, developmental delay, autism, and epilepsy. Six further cases with a primary diagnosis of autism spectrum disorder (ASD) and similar secondary phenotypes were identified with ELP4 deletions, as well as another six (out of nine) with neurodevelopmental phenotypes from DECIPHER. CNVs at ELP4 were only present in 1/11,252 controls. We found a significant excess of CNVs in discovery cases compared with controls, P = 7.5 × 10(-3) , as well as for autism, P = 2.7 × 10(-3) . Our results suggest that ELP4 deletions are highly likely to be pathogenic, predisposing to a range of neurodevelopmental phenotypes from ASD to language impairment and epilepsy.

SHANK1 Deletions in Males with Autism Spectrum Disorder.

Recent studies have highlighted the involvement of rare (<1% frequency) copy-number variations and point mutations in the genetic etiology of autism spectrum disorder (ASD); these variants particularly affect genes involved in the neuronal synaptic complex. The SHANK gene family consists of three members (SHANK1, SHANK2, and SHANK3), which encode scaffolding proteins required for the proper formation and function of neuronal synapses. Although SHANK2 and SHANK3 mutations have been implicated in ASD and intellectual disability, the involvement of SHANK1 is unknown. Here, we assess microarray data from 1,158 Canadian and 456 European individuals with ASD to discover microdeletions at the SHANK1 locus on chromosome 19. We identify a hemizygous SHANK1 deletion that segregates in a four-generation family in which male carriers--but not female carriers--have ASD with higher functioning. A de novo SHANK1 deletion was also detected in an unrelated male individual with ASD with higher functioning, and no equivalent SHANK1 mutations were found in >15,000 controls (p = 0.009). The discovery of apparent reduced penetrance of ASD in females bearing inherited autosomal SHANK1 deletions provides a possible contributory model for the male gender bias in autism. The data are also informative for clinical-genetics interpretations of both inherited and sporadic forms of ASD involving SHANK1.

Rare deletions at the neurexin 3 locus in autism spectrum disorder.

The three members of the human neurexin gene family, neurexin 1 (NRXN1), neurexin 2 (NRXN2), and neurexin 3 (NRXN3), encode neuronal adhesion proteins that have important roles in synapse development and function. In autism spectrum disorder (ASD), as well as in other neurodevelopmental conditions, rare exonic copy-number variants and/or point mutations have been identified in the NRXN1 and NRXN2 loci. We present clinical characterization of four index cases who have been diagnosed with ASD and who possess rare inherited or de novo microdeletions at 14q24.3-31.1, a region that overlaps exons of the alpha and/or beta isoforms of NRXN3. NRXN3 deletions were found in one father with subclinical autism and in a carrier mother and father without formal ASD diagnoses, indicating issues of penetrance and expressivity at this locus. Notwithstanding these clinical complexities, this report on ASD-affected individuals who harbor NRXN3 exonic deletions advances the understanding of the genetic etiology of autism, further enabling molecular diagnoses.

Genetic and functional analyses of SHANK2 mutations suggest a multiple hit model of autism spectrum disorders.

Autism spectrum disorders (ASD) are a heterogeneous group of neurodevelopmental disorders with a complex inheritance pattern. While many rare variants in synaptic proteins have been identified in patients with ASD, little is known about their effects at the synapse and their interactions with other genetic variations. Here, following the discovery of two de novo SHANK2 deletions by the Autism Genome Project, we identified a novel 421 kb de novo SHANK2 deletion in a patient with autism. We then sequenced SHANK2 in 455 patients with ASD and 431 controls and integrated these results with those reported by Berkel et al. 2010 (n = 396 patients and n = 659 controls). We observed a significant enrichment of variants affecting conserved amino acids in 29 of 851 (3.4%) patients and in 16 of 1,090 (1.5%) controls (P = 0.004, OR = 2.37, 95% CI = 1.23-4.70). In neuronal cell cultures, the variants identified in patients were associated with a reduced synaptic density at dendrites compared to the variants only detected in controls (P = 0.0013). Interestingly, the three patients with de novo SHANK2 deletions also carried inherited CNVs at 15q11-q13 previously associated with neuropsychiatric disorders. In two cases, the nicotinic receptor CHRNA7 was duplicated and in one case the synaptic translation repressor CYFIP1 was deleted. These results strengthen the role of synaptic gene dysfunction in ASD but also highlight the presence of putative modifier genes, which is in keeping with the "multiple hit model" for ASD. A better knowledge of these genetic interactions will be necessary to understand the complex inheritance pattern of ASD.

Individual common variants exert weak effects on the risk for autism spectrum disorders.

While it is apparent that rare variation can play an important role in the genetic architecture of autism spectrum disorders (ASDs), the contribution of common variation to the risk of developing ASD is less clear. To produce a more comprehensive picture, we report Stage 2 of the Autism Genome Project genome-wide association study, adding 1301 ASD families and bringing the total to 2705 families analysed (Stages 1 and 2). In addition to evaluating the association of individual single nucleotide polymorphisms (SNPs), we also sought evidence that common variants, en masse, might affect the risk. Despite genotyping over a million SNPs covering the genome, no single SNP shows significant association with ASD or selected phenotypes at a genome-wide level. The SNP that achieves the smallest P-value from secondary analyses is rs1718101. It falls in CNTNAP2, a gene previously implicated in susceptibility for ASD. This SNP also shows modest association with age of word/phrase acquisition in ASD subjects, of interest because features of language development are also associated with other variation in CNTNAP2. In contrast, allele scores derived from the transmission of common alleles to Stage 1 cases significantly predict case status in the independent Stage 2 sample. Despite being significant, the variance explained by these allele scores was small (Vm< 1%). Based on results from individual SNPs and their en masse effect on risk, as inferred from the allele score results, it is reasonable to conclude that common variants affect the risk for ASD but their individual effects are modest.

A novel approach identifies new differentially methylated regions (DMRs) associated with imprinted genes.

Imprinted genes are critical for normal human growth and neurodevelopment. They are characterized by differentially methylated regions (DMRs) of DNA that confer parent of origin-specific transcription. We developed a new strategy to identify imprinted gene-associated DMRs. Using genome-wide methylation profiling of sodium bisulfite modified DNA from normal human tissues of biparental origin, candidate DMRs were identified by selecting CpGs with methylation levels consistent with putative allelic differential methylation. In parallel, the methylation profiles of tissues of uniparental origin, i.e., paternally-derived androgenetic complete hydatidiform moles (AnCHMs), and maternally-derived mature cystic ovarian teratoma (MCT), were examined and then used to identify CpGs with parent of origin-specific DNA methylation. With this approach, we found known DMRs associated with imprinted genomic regions as well as new DMRs for known imprinted genes, NAP1L5 and ZNF597, and novel candidate imprinted genes. The paternally methylated DMR for one candidate, AXL, a receptor tyrosine kinase, was also validated in experiments with mouse embryos that demonstrated Axl was expressed preferentially from the maternal allele in a DNA methylation-dependent manner.

Transcriptomic landscape of human induced pluripotent stem cell-derived osteogenic differentiation identifies a regulatory role of KLF16.

Osteogenic differentiation is essential for bone development and metabolism, but the underlying gene regulatory networks have not been well investigated. We differentiated mesenchymal stem cells, derived from 20 human induced pluripotent stem cell lines, into preosteoblasts and osteoblasts, and performed systematic RNA-seq analyses of 60 samples for differential gene expression. We noted a highly significant correlation in expression patterns and genomic proximity among transcription factor (TF) and long noncoding RNA (lncRNA) genes. We identified TF-TF regulatory networks, regulatory roles of lncRNAs on their neighboring coding genes for TFs and splicing factors, and differential splicing of TF, lncRNA, and splicing factor genes. TF-TF regulatory and gene co-expression network analyses suggested an inhibitory role of TF KLF16 in osteogenic differentiation. We demonstrate that in vitro overexpression of human KLF16 inhibits osteogenic differentiation and mineralization, and in vivo Klf16+/- mice exhibit increased bone mineral density, trabecular number, and cortical bone area. Thus, our model system highlights the regulatory complexity of osteogenic differentiation and identifies novel osteogenic genes.