The interplay between the transforming growth factor β (TGF-β) signaling proteins, SMAD family member 2 (SMAD2) and 3 (SMAD3), and the TGF-β-inhibiting SMAD, SMAD7, seems to play a vital role in proper pancreatic endocrine development and also in normal β-cell function in adult pancreatic islets. Here, we generated conditional SMAD7 knockout mice by crossing insulin1Cre mice with SMAD7fx/fx mice. We also created a β cell-specific SMAD7-overexpressing mouse line by crossing insulin1Dre mice with HPRT-SMAD7/RosaGFP mice. We analyzed β-cell function in adult islets when SMAD7 was either absent or overexpressed in β cells. Loss of SMAD7 in β cells inhibited proliferation, and SMAD7 overexpression enhanced cell proliferation. However, alterations in basic glucose homeostasis were not detectable following either SMAD7 deletion or overexpression in β cells. Our results show that both the absence and overexpression of SMAD7 affect TGF-β signaling and modulates β-cell proliferation but does not appear to alter β-cell function. Reversible SMAD7 overexpression may represent an attractive therapeutic option to enhance β-cell proliferation without negative effects on β-cell function.

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.

Many immune diseases occur at different rates among people with schizophrenia compared to the general population. Here, we evaluated whether this phenomenon might be explained by shared genetic risk factors. We used data from large genome-wide association studies to compare the genetic architecture of schizophrenia to 19 immune diseases. First, we evaluated the association with schizophrenia of 581 variants previously reported to be associated with immune diseases at genome-wide significance. We identified five variants with potentially pleiotropic effects. While colocalization analyses were inconclusive, functional characterization of these variants provided the strongest evidence for a model in which genetic variation at rs1734907 modulates risk of schizophrenia and Crohn's disease via altered methylation and expression of EPHB4-a gene whose protein product guides the migration of neuronal axons in the brain and the migration of lymphocytes towards infected cells in the immune system. Next, we investigated genome-wide sharing of common variants between schizophrenia and immune diseases using cross-trait LD score regression. Of the 11 immune diseases with available genome-wide summary statistics, we observed genetic correlation between six immune diseases and schizophrenia: inflammatory bowel disease (rg = 0.12 ± 0.03, P = 2.49 × 10-4), Crohn's disease (rg = 0.097 ± 0.06, P = 3.27 × 10-3), ulcerative colitis (rg = 0.11 ± 0.04, P = 4.05 × 10-3), primary biliary cirrhosis (rg = 0.13 ± 0.05, P = 3.98 × 10-3), psoriasis (rg = 0.18 ± 0.07, P = 7.78 × 10-3) and systemic lupus erythematosus (rg = 0.13 ± 0.05, P = 3.76 × 10-3). With the exception of ulcerative colitis, the degree and direction of these genetic correlations were consistent with the expected phenotypic correlation based on epidemiological data. Our findings suggest shared genetic risk factors contribute to the epidemiological association of certain immune diseases and schizophrenia.

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.

Congenital heart defects (CHDs) affect approximately 1% of newborns. Epidemiological studies have identified several genetically-mediated maternal phenotypes (e.g., pregestational diabetes, chronic hypertension) that are associated with the risk of CHDs in offspring. However, the role of the maternal genome in determining CHD risk has not been defined. We present findings from gene-level, genome-wide studies that link CHDs to maternal effect genes as well as to maternal genes related to hypertension and proteostasis. Maternal effect genes, which provide the mRNAs and proteins in the oocyte that guide early embryonic development before zygotic gene activation, have not previously been implicated in CHD risk. Our findings support a role for and suggest new pathways by which the maternal genome may contribute to the development of CHDs in offspring.

The genetic background of childhood body mass index (BMI), and the extent to which the well-known associations of childhood BMI with adult diseases are explained by shared genetic factors, are largely unknown. We performed a genome-wide association study meta-analysis of BMI in 61,111 children aged between 2 and 10 years. Twenty-five independent loci reached genome-wide significance in the combined discovery and replication analyses. Two of these, located near NEDD4L and SLC45A3, have not previously been reported in relation to either childhood or adult BMI. Positive genetic correlations of childhood BMI with birth weight and adult BMI, waist-to-hip ratio, diastolic blood pressure and type 2 diabetes were detected (Rg ranging from 0.11 to 0.76, P-values <0.002). A negative genetic correlation of childhood BMI with age at menarche was observed. Our results suggest that the biological processes underlying childhood BMI largely, but not completely, overlap with those underlying adult BMI. The well-known observational associations of BMI in childhood with cardio-metabolic diseases in adulthood may reflect partial genetic overlap, but in light of previous evidence, it is also likely that they are explained through phenotypic continuity of BMI from childhood into adulthood.

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 new lipid-associated genes and study the relationships among lipids, genotypes, gene expression and hundreds of complex human diseases from the Electronic Medical Records and Genomics (347 traits) and the UK Biobank (549 traits). Aside from 67 new 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.

Although many loci have been associated with height in European ancestry populations, very few have been identified in African ancestry individuals. Furthermore, many of the known loci have yet to be generalized to and fine-mapped within a large-scale African ancestry sample. We performed sex-combined and sex-stratified meta-analyses in up to 52,764 individuals with height and genome-wide genotyping data from the African Ancestry Anthropometry Genetics Consortium (AAAGC). We additionally combined our African ancestry meta-analysis results with published European genome-wide association study (GWAS) data. In the African ancestry analyses, we identified three novel loci (SLC4A3, NCOA2, ECD/FAM149B1) in sex-combined results and two loci (CRB1, KLF6) in women only. In the African plus European sex-combined GWAS, we identified an additional three novel loci (RCCD1, G6PC3, CEP95) which were equally driven by AAAGC and European results. Among 39 genome-wide significant signals at known loci, conditioning index SNPs from European studies identified 20 secondary signals. Two of the 20 new secondary signals and none of the 8 novel loci had minor allele frequencies (MAF) < 5%. Of 802 known European height signals, 643 displayed directionally consistent associations with height, of which 205 were nominally significant (p < 0.05) in the African ancestry sex-combined sample. Furthermore, 148 of 241 loci contained ≤20 variants in the credible sets that jointly account for 99% of the posterior probability of driving the associations. In summary, trans-ethnic meta-analyses revealed novel signals and further improved fine-mapping of putative causal variants in loci shared between African and European ancestry populations.

Hereditary cystatin C amyloid angiopathy is a dominantly inherited disease caused by a leucine to glutamine variant of human cystatin C (hCC). L68Q-hCC forms amyloid deposits in brain arteries associated with micro-infarcts, leading ultimately to paralysis, dementia and death in young adults. To evaluate the ability of molecules to interfere with aggregation of hCC while informing about cellular toxicity, we generated cells that produce and secrete WT and L68Q-hCC and have detected high-molecular weight complexes formed from the mutant protein. Incubations of either lysate or supernatant containing L68Q-hCC with reducing agents glutathione or N-acetyl-cysteine (NAC) breaks oligomers into monomers. Six L68Q-hCC carriers taking NAC had skin biopsies obtained to determine if hCC deposits were reduced following NAC treatment. Remarkably, ~50-90% reduction of L68Q-hCC staining was observed in five of the treated carriers suggesting that L68Q-hCC is a clinical target for reducing agents.

Clustering is a critical step in single cell-based studies. Most existing methods support unsupervised clustering without the a priori exploitation of any domain knowledge. When confronted by the high dimensionality and pervasive dropout events of scRNA-Seq data, purely unsupervised clustering methods may not produce biologically interpretable clusters, which complicates cell type assignment. In such cases, the only recourse is for the user to manually and repeatedly tweak clustering parameters until acceptable clusters are found. Consequently, the path to obtaining biologically meaningful clusters can be ad hoc and laborious. Here we report a principled clustering method named scDCC, that integrates domain knowledge into the clustering step. Experiments on various scRNA-seq datasets from thousands to tens of thousands of cells show that scDCC can significantly improve clustering performance, facilitating the interpretability of clusters and downstream analyses, such as cell type assignment.

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.

In contrast to the skin and the gut, where somatic stem cells and their niche are well characterized, a definitive pancreatic multipotent cell population in the adult pancreas has yet to be revealed. Of particular interest is whether such cells may be endogenous in patients with diabetes, and if so, can they be used for therapeutic purposes? In the current study, we used two separate reporter lines to target Cre-recombinase expression to the Lgr5- or glucagon-expressing cells in the pancreas. We provide evidence for the existence of a population of cells within and in the proximity of the ducts that transiently express the stem-cell marker Lgr5 during late gestational stages. Careful timing of tamoxifen treatment in Lgr5EGFP-IRES-CreERT2 ;R26 Tomato mice allowed us to show that these Lgr5-expressing progenitor cells can differentiate into α-cells during pregnancy. Furthermore, we report on a spontaneous lineage conversion of α- to β-cells specifically after parturition. The contribution of Lgr5 progeny to the β-cell compartment through an α-cell intermediate phase early after pregnancy appears to be part of a novel mechanism that would counterbalance against excessive β-cell mass reduction during β-cell involution.

Suicide is a leading cause of death worldwide, and nonfatal suicide attempts, which occur far more frequently, are a major source of disability and social and economic burden. Both have substantial genetic etiology, which is partially shared and partially distinct from that of related psychiatric disorders.

With the advances in single-cell sequencing techniques, numerous analytical methods have been developed for delineating cell development. However, most are based on Euclidean space, which would distort the complex hierarchical structure of cell differentiation. Recently, methods acting on hyperbolic space have been proposed to visualize hierarchical structures in single-cell RNA-seq (scRNA-seq) data and have been proven to be superior to methods acting on Euclidean space. However, these methods have fundamental limitations and are not optimized for the highly sparse single-cell count data. To address these limitations, we propose scDHMap, a model-based deep learning approach to visualize the complex hierarchical structures of scRNA-seq data in low-dimensional hyperbolic space. The evaluations on extensive simulation and real experiments show that scDHMap outperforms existing dimensionality-reduction methods in various common analytical tasks as needed for scRNA-seq data, including revealing trajectory branches, batch correction, and denoising the count matrix with high dropout rates. In addition, we extend scDHMap to visualize single-cell ATAC-seq data.

How abnormal neurodevelopment relates to the tumour aggressiveness of medulloblastoma (MB), the most common type of embryonal tumour, remains elusive. Here we uncover a neurodevelopmental epigenomic programme that is hijacked to induce MB metastatic dissemination. Unsupervised analyses of integrated publicly available datasets with our newly generated data reveal that SMARCD3 (also known as BAF60C) regulates Disabled 1 (DAB1)-mediated Reelin signalling in Purkinje cell migration and MB metastasis by orchestrating cis-regulatory elements at the DAB1 locus. We further identify that a core set of transcription factors, enhancer of zeste homologue 2 (EZH2) and nuclear factor I X (NFIX), coordinates with the cis-regulatory elements at the SMARCD3 locus to form a chromatin hub to control SMARCD3 expression in the developing cerebellum and in metastatic MB. Increased SMARCD3 expression activates Reelin-DAB1-mediated Src kinase signalling, which results in a MB response to Src inhibition. These data deepen our understanding of how neurodevelopmental programming influences disease progression and provide a potential therapeutic option for patients with MB.

Understanding signaling pathways that regulate pancreatic β-cell function to produce, store, and release insulin, as well as pathways that control β-cell proliferation, is vital to find new treatments for diabetes mellitus. Transforming growth factor-beta (TGF-β) signaling is involved in a broad range of β-cell functions. The canonical TGF-β signaling pathway functions through intracellular smads, including smad2 and smad3, to regulate cell development, proliferation, differentiation, and function in many organs. Here, we demonstrate the role of TGF-β/smad2 signaling in regulating mature β-cell proliferation and function using β-cell-specific smad2 null mutant mice. β-cell-specific smad2-deficient mice exhibited improved glucose clearance as demonstrated by glucose tolerance testing, enhanced in vivo and ex vivo glucose-stimulated insulin secretion, and increased β-cell mass and proliferation. Furthermore, when these mice were fed a high-fat diet to induce hyperglycemia, they again showed improved glucose tolerance, insulin secretion, and insulin sensitivity. In addition, ex vivo analysis of smad2-deficient islets showed that they displayed increased glucose-stimulated insulin secretion and upregulation of genes involved in insulin synthesis and insulin secretion. Thus, we conclude that smad2 could represent an attractive therapeutic target for type 2 diabetes mellitus.

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.

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 up-regulated in the CNS, while IL-34 remains constitutively expressed. We found that targeting CSF-1 with neutralizing antibody 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. Both anti–CSF-1 and BLZ945 treatment greatly reduced the number 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 treatment reduced the size of demyelinated lesions and microglial activation in the gray 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 ablation of CSF-1R–dependent cells.

SYNJ1 encodes Synaptojanin-1, a dual-function poly-phosphoinositide phosphatase that is expressed in the brain to regulate neuronal synaptic vesicle dynamics. Biallelic SYNJ1 variants cause a spectrum of clinical manifestations, from early onset parkinsonism to developmental and epileptic encephalopathy.