The association of GNB5 with Alzheimer disease revealed by genomic analysis restricted to variants impacting gene function.

Disease-associated variants identified from genome-wide association studies (GWASs) frequently map to non-coding areas of the genome such as introns and intergenic regions. An exclusive reliance on gene-agnostic methods of genomic investigation could limit the identification of relevant genes associated with polygenic diseases such as Alzheimer disease (AD). To overcome such potential restriction, we developed a gene-constrained analytical method that considers only moderate- and high-risk variants that affect gene coding sequences. We report here the application of this approach to publicly available datasets containing 181,388 individuals without and with AD and the resulting identification of 660 genes potentially linked to the higher AD prevalence among Africans/African Americans. By integration with transcriptome analysis of 23 brain regions from 2,728 AD case-control samples, we concentrated on nine genes that potentially enhance the risk of AD: AACS, GNB5, GNS, HIPK3, MED13, SHC2, SLC22A5, VPS35, and ZNF398. GNB5, the fifth member of the heterotrimeric G protein beta family encoding Gß5, is primarily expressed in neurons and is essential for normal neuronal development in mouse brain. Homozygous or compound heterozygous loss of function of GNB5 in humans has previously been associated with a syndrome of developmental delay, cognitive impairment, and cardiac arrhythmia. In validation experiments, we confirmed that Gnb5 heterozygosity enhanced the formation of both amyloid plaques and neurofibrillary tangles in the brains of AD model mice. These results suggest that gene-constrained analysis can complement the power of GWASs in the identification of AD-associated genes and may be more broadly applicable to other polygenic diseases.

NIAGADS Alzheimer's GenomicsDB: A resource for exploring Alzheimer's disease genetic and genomic knowledge.

INTRODUCTION: The National Institute on Aging Genetics of Alzheimer's Disease Data Storage Site Alzheimer's Genomics Database (GenomicsDB) is a public knowledge base of Alzheimer's disease (AD) genetic datasets and genomic annotations. METHODS: GenomicsDB uses a custom systems architecture to adopt and enforce rigorous standards that facilitate harmonization of AD-relevant genome-wide association study summary statistics datasets with functional annotations, including over 230 million annotated variants from the AD Sequencing Project. RESULTS: GenomicsDB generates interactive reports compiled from the harmonized datasets and annotations. These reports contextualize AD-risk associations in a broader functional genomic setting and summarize them in the context of functionally annotated genes and variants. DISCUSSION: Created to make AD-genetics knowledge more accessible to AD researchers, the GenomicsDB is designed to guide users unfamiliar with genetic data in not only exploring but also interpreting this ever-growing volume of data. Scalable and interoperable with other genomics resources using data technology standards, the GenomicsDB can serve as a central hub for research and data analysis on AD and related dementias. HIGHLIGHTS: The National Institute on Aging Genetics of Alzheimer's Disease Data Storage Site (NIAGADS) offers to the public a unique, disease-centric collection of AD-relevant GWAS summary statistics datasets. Interpreting these data is challenging and requires significant bioinformatics expertise to standardize datasets and harmonize them with functional annotations on genome-wide scales. The NIAGADS Alzheimer's GenomicsDB helps overcome these challenges by providing a user-friendly public knowledge base for AD-relevant genetics that shares harmonized, annotated summary statistics datasets from the NIAGADS repository in an interpretable, easily searchable format.

Scalable approaches for functional analyses of whole-genome sequencing non-coding variants.

Non-coding genetic variants outside of protein-coding genome regions play an important role in genetic and epigenetic regulation. It has become increasingly important to understand their roles, as non-coding variants often make up the majority of top findings of genome-wide association studies (GWAS). In addition, the growing popularity of disease-specific whole-genome sequencing (WGS) efforts expands the library of and offers unique opportunities for investigating both common and rare non-coding variants, which are typically not detected in more limited GWAS approaches. However, the sheer size and breadth of WGS data introduce additional challenges to predicting functional impacts in terms of data analysis and interpretation. This review focuses on the recent approaches developed for efficient, at-scale annotation and prioritization of non-coding variants uncovered in WGS analyses. In particular, we review the latest scalable annotation tools, databases and functional genomic resources for interpreting the variant findings from WGS based on both experimental data and in silico predictive annotations. We also review machine learning-based predictive models for variant scoring and prioritization. We conclude with a discussion of future research directions which will enhance the data and tools necessary for the effective functional analyses of variants identified by WGS to improve our understanding of disease etiology.

FILER: a framework for harmonizing and querying large-scale functional genomics knowledge.

Querying massive functional genomic and annotation data collections, linking and summarizing the query results across data sources/data types are important steps in high-throughput genomic and genetic analytical workflows. However, these steps are made difficult by the heterogeneity and breadth of data sources, experimental assays, biological conditions/tissues/cell types and file formats. FILER (FunctIonaL gEnomics Repository) is a framework for querying large-scale genomics knowledge with a large, curated integrated catalog of harmonized functional genomic and annotation data coupled with a scalable genomic search and querying interface. FILER uniquely provides: (i) streamlined access to >50 000 harmonized, annotated genomic datasets across >20 integrated data sources, >1100 tissues/cell types and >20 experimental assays; (ii) a scalable genomic querying interface; and (iii) ability to analyze and annotate user's experimental data. This rich resource spans >17 billion GRCh37/hg19 and GRCh38/hg38 genomic records. Our benchmark querying 7 × 109 hg19 FILER records shows FILER is highly scalable, with a sub-linear 32-fold increase in querying time when increasing the number of queries 1000-fold from 1000 to 1 000 000 intervals. Together, these features facilitate reproducible research and streamline integrating/querying large-scale genomic data within analyses/workflows. FILER can be deployed on cloud or local servers (https://bitbucket.org/wanglab-upenn/FILER) for integration with custom pipelines and is freely available (https://lisanwanglab.org/FILER).

A developmental lineage-based gene co-expression network for mouse pancreatic ß-cells reveals a role for Zfp800 in pancreas development.

To gain a deeper understanding of pancreatic ß-cell development, we used iterative weighted gene correlation network analysis to calculate a gene co-expression network (GCN) from 11 temporally and genetically defined murine cell populations. The GCN, which contained 91 distinct modules, was then used to gain three new biological insights. First, we found that the clustered protocadherin genes are differentially expressed during pancreas development. Pcdhγ genes are preferentially expressed in pancreatic endoderm, Pcdhß genes in nascent islets, and Pcdhα genes in mature ß-cells. Second, after extracting sub-networks of transcriptional regulators for each developmental stage, we identified 81 zinc finger protein (ZFP) genes that are preferentially expressed during endocrine specification and ß-cell maturation. Third, we used the GCN to select three ZFPs for further analysis by CRISPR mutagenesis of mice. Zfp800 null mice exhibited early postnatal lethality, and at E18.5 their pancreata exhibited a reduced number of pancreatic endocrine cells, alterations in exocrine cell morphology, and marked changes in expression of genes involved in protein translation, hormone secretion and developmental pathways in the pancreas. Together, our results suggest that developmentally oriented GCNs have utility for gaining new insights into gene regulation during organogenesis.

SparkINFERNO: a scalable high-throughput pipeline for inferring molecular mechanisms of non-coding genetic variants.

SUMMARY: We report Spark-based INFERence of the molecular mechanisms of NOn-coding genetic variants (SparkINFERNO), a scalable bioinformatics pipeline characterizing non-coding genome-wide association study (GWAS) association findings. SparkINFERNO prioritizes causal variants underlying GWAS association signals and reports relevant regulatory elements, tissue contexts and plausible target genes they affect. To achieve this, the SparkINFERNO algorithm integrates GWAS summary statistics with large-scale collection of functional genomics datasets spanning enhancer activity, transcription factor binding, expression quantitative trait loci and other functional datasets across more than 400 tissues and cell types. Scalability is achieved by an underlying API implemented using Apache Spark and Giggle-based genomic indexing. We evaluated SparkINFERNO on large GWASs and show that SparkINFERNO is more than 60 times efficient and scales with data size and amount of computational resources. AVAILABILITY AND IMPLEMENTATION: SparkINFERNO runs on clusters or a single server with Apache Spark environment, and is available at https://bitbucket.org/wanglab-upenn/SparkINFERNO or https://hub.docker.com/r/wanglab/spark-inferno. CONTACT: lswang@pennmedicine.upenn.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Phospho-proteomic discovery of novel signal transducers including thioredoxin-interacting protein as mediators of erythropoietin-dependent human erythropoiesis.

Erythroid cell formation critically depends on signals transduced via erythropoietin (EPO)/EPO receptor (EPOR)/JAK2 complexes. This includes not only core response modules (e.g., JAK2/STAT5, RAS/MEK/ERK), but also specialized effectors (e.g., erythroferrone, ASCT2 glutamine transport, Spi2A). By using phospho-proteomics and a human erythroblastic cell model, we identify 121 new EPO target proteins, together with their EPO-modulated domains and phosphosites. Gene ontology (GO) enrichment for "Molecular Function" identified adaptor proteins as one top EPO target category. This includes a novel EPOR/JAK2-coupled network of actin assemblage modifiers, with adaptors DLG-1, DLG-3, WAS, WASL, and CD2AP as prime components. "Cellular Component" GO analysis further identified 19 new EPO-modulated cytoskeletal targets including the erythroid cytoskeletal targets spectrin A, spectrin B, adducin 2, and glycophorin C. In each, EPO-induced phosphorylation occurred at pY sites and subdomains, which suggests coordinated regulation by EPO of the erythroid cytoskeleton. GO analysis of "Biological Processes" further revealed metabolic regulators as a likewise unexpected EPO target set. Targets included aldolase A, pyruvate dehydrogenase α1, and thioredoxin-interacting protein (TXNIP), with EPO-modulated p-Y sites in each occurring within functional subdomains. In TXNIP, EPO-induced phosphorylation occurred at novel p-T349 and p-S358 sites, and was paralleled by rapid increases in TXNIP levels. In UT7epo-E and primary human stem cell (HSC)-derived erythroid progenitor cells, lentivirus-mediated short hairpin RNA knockdown studies revealed novel pro-erythropoietic roles for TXNIP. Specifically, TXNIP's knockdown sharply inhibited c-KIT expression; compromised EPO dose-dependent erythroblast proliferation and survival; and delayed late-stage erythroblast formation. Overall, new insight is provided into EPO's diverse action mechanisms and TXNIP's contributions to EPO-dependent human erythropoiesis.

Phospho-PTM proteomic discovery of novel EPO- modulated kinases and phosphatases, including PTPN18 as a positive regulator of EPOR/JAK2 Signaling.

The formation of erythroid progenitor cells depends sharply upon erythropoietin (EPO), its cell surface receptor (erythropoietin receptor, EPOR), and Janus kinase 2 (JAK2). Clinically, recombinant human EPO (rhEPO) additionally is an important anti-anemia agent for chronic kidney disease (CKD), myelodysplastic syndrome (MDS) and chemotherapy, but induces hypertension, and can exert certain pro-tumorigenic effects. Cellular signals transduced by EPOR/JAK2 complexes, and the nature of EPO-modulated signal transduction factors, therefore are of significant interest. By employing phospho-tyrosine post-translational modification (p-Y PTM) proteomics and human EPO- dependent UT7epo cells, we have identified 22 novel kinases and phosphatases as novel EPO targets, together with their specific sites of p-Y modification. New kinases modified due to EPO include membrane palmitoylated protein 1 (MPP1) and guanylate kinase 1 (GUK1) guanylate kinases, together with the cytoskeleton remodeling kinases, pseudopodium enriched atypical kinase 1 (PEAK1) and AP2 associated kinase 1 (AAK1). Novel EPO- modified phosphatases include protein tyrosine phosphatase receptor type A (PTPRA), phosphohistidine phosphatase 1 (PHPT1), tensin 2 (TENC1), ubiquitin associated and SH3 domain containing B (UBASH3B) and protein tyrosine phosphatase non-receptor type 18 (PTPN18). Based on PTPN18's high expression in hematopoietic progenitors, its novel connection to JAK kinase signaling, and a unique EPO- regulated PTPN18-pY389 motif which is modulated by JAK2 inhibitors, PTPN18's actions in UT7epo cells were investigated. Upon ectopic expression, wt-PTPN18 promoted EPO dose-dependent cell proliferation, and survival. Mechanistically, PTPN18 sustained the EPO- induced activation of not only mitogen-activated protein kinases 1 and 3 (ERK1/2), AKT serine/threonine kinase 1-3 (AKT), and signal transducer and activator of transcription 5A and 5B (STAT5), but also JAK2. Each effect further proved to depend upon PTPN18's EPO- modulated (p)Y389 site. In analyses of the EPOR and the associated adaptor protein RHEX (regulator of hemoglobinization and erythroid cell expansion), wt-PTPN18 increased high molecular weight EPOR forms, while sharply inhibiting the EPO-induced phosphorylation of RHEX-pY141. Each effect likewise depended upon PTPN18-Y389. PTPN18 thus promotes signals for EPO-dependent hematopoietic cell growth, and may represent a new druggable target for myeloproliferative neoplasms.

Demographic disequilibrium in living nautiloids (Nautilus and Allonautilus): Canary in the coal mines?

Averaged demographic data from previously unfished populations of Nautilus and Allonautilus (Cephalopoda) provide a baseline to determine if a population is undisturbed and in "equilibrium" or is in "disequilibrium" as a result of fishery pressure. Data are available for previously undisturbed local nautiloid populations in Papua New Guinea, Australia, Indonesia, Fiji, Palau, American Samoa, New Caledonia and Vanuatu (total n = 2,669 live-caught, tagged and released animals). The data show that unfished populations average ~75% males and ~74% mature animals. By contrast, unpublished, anecdotal and historical records since 1900 from the heavily fished central Philippines have shown a persistent decline in trap yields and a change in demographics of N. pompilius. By 1979, a sample of fished live-caught animals (n = 353) comprised only ~28% males and ~27% mature animals. Continued uncontrolled trapping caused collapse of the fishery and the shell industry has moved elsewhere, including Indonesia. In addition, we show that estimated rates of population decline are offered by unpublished tag-release records in unfished Palau. These data show that patterns of trap yields and demographic differences between fished and unfished populations in relative age class and sex ratios can indicate disequilibria wrought by fisheries pressure that can render local populations inviable. Given adequate samples (n ≥100 live-caught animals), a threshold of <50% males and mature animals in fished populations should signal the need to initiate curative conservation initiatives. The current trajectory of uncontrolled nautiloid fisheries can only mean trouble and possibly extinction of local populations of this ancient, iconic molluscan lineage.

Governing roles for Trib3 pseudokinase during stress erythropoiesis.

In response to anemia, the heightened production of erythropoietin (EPO) can sharply promote erythroid progenitor cell (EPC) formation. Specific mediators of such EPO- accelerated erythropoiesis, however, are not well understood. Presently, we first report that the expression of Trib3 in adult bone marrow EPCs in vivo is nominal at steady state, but strongly activated on EPO challenge. In a knockout mouse model, Trib3 disruption modestly increased steady-state erythrocyte numbers and decreased mean corpuscular volume. Following 5-fluorouracil myeloablation, however, rebound red blood cell production and hemoglobin levels were substantially (and selectively) compromised in Trib3-/- mice versus Trib3+/+ congenic controls. Erythrocytes from 5-fluorouracil-treated Trib3-/- mice additionally were more prone to lysis and exhibited elevated peroxide-induced reactive oxygen species. Ex vivo, the development of CD71posTer119pos erythroblasts from Trib3-/- bone marrow progenitors was attenuated, and this was associated with heightened EPO-dependent Erk1/2 activation and moderately increased Akt activation. For developmentally staged EPCs, gene profiling provided further initial insight into candidate mediators of EPO-induced Trib3 gene expression, including Cebp-beta, Atf4, Egr-1, and Nab1. Overall, Trib3 is indicated to act as a novel EPC-intrinsic governor of stress erythropoiesis.

Stat and interferon genes identified by network analysis differentially regulate primitive and definitive erythropoiesis.

BACKGROUND: Hematopoietic ontogeny is characterized by overlapping waves of primitive, fetal definitive, and adult definitive erythroid lineages. Our aim is to identify differences in the transcriptional control of these distinct erythroid cell maturation pathways by inferring and analyzing gene-interaction networks from lineage-specific expression datasets. Inferred networks are strongly connected and do not fit a scale-free model, making it difficult to identify essential regulators using the hub-essentiality standard. RESULTS: We employed a semi-supervised machine learning approach to integrate measures of network topology with expression data to score gene essentiality. The algorithm was trained and tested on the adult and fetal definitive erythroid lineages. When applied to the primitive erythroid lineage, 144 high scoring transcription factors were found to be differentially expressed between the primitive and adult definitive erythroid lineages, including all expressed STAT-family members. Differential responses of primitive and definitive erythroblasts to a Stat3 inhibitor and IFNγ in vitro supported the results of the computational analysis. Further investigation of the original expression data revealed a striking signature of Stat1-related genes in the adult definitive erythroid network. Among the potential pathways known to utilize Stat1, interferon (IFN) signaling-related genes were expressed almost exclusively within the adult definitive erythroid network. CONCLUSIONS: In vitro results support the computational prediction that differential regulation and downstream effectors of STAT signaling are key factors that distinguish the transcriptional control of primitive and definitive erythroid cell maturation.

Ontogeny of erythroid gene expression.

Erythroid ontogeny is characterized by overlapping waves of primitive and definitive erythroid lineages that share many morphologic features during terminal maturation but have marked differences in cell size and globin expression. In the present study, we compared global gene expression in primitive, fetal definitive, and adult definitive erythroid cells at morphologically equivalent stages of maturation purified from embryonic, fetal, and adult mice. Surprisingly, most transcriptional complexity in erythroid precursors is already present by the proerythroblast stage. Transcript levels are markedly modulated during terminal erythroid maturation, but housekeeping genes are not preferentially lost. Although primitive and definitive erythroid lineages share a large set of nonhousekeeping genes, annotation of lineage-restricted genes shows that alternate gene usage occurs within shared functional categories, as exemplified by the selective expression of aquaporins 3 and 8 in primitive erythroblasts and aquaporins 1 and 9 in adult definitive erythroblasts. Consistent with the known functions of Aqp3 and Aqp8 as H2O2 transporters, primitive, but not definitive, erythroblasts preferentially accumulate reactive oxygen species after exogenous H2O2 exposure. We have created a user-friendly Web site (http://www.cbil.upenn.edu/ErythronDB) to make these global expression data readily accessible and amenable to complex search strategies by the scientific community.