Most Pathways Can Be Related to the Pathogenesis of Alzheimer's Disease.

Alzheimer's disease (AD) is a complex neurodegenerative disorder. The relative contribution of the numerous underlying functional mechanisms is poorly understood. To comprehensively understand the context and distribution of pathways that contribute to AD, we performed text-mining to generate an exhaustive, systematic assessment of the breadth and diversity of biological pathways within a corpus of 206,324 dementia publication abstracts. A total of 91% (325/335) of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways have publications containing an association via at least 5 studies, while 63% of pathway terms have at least 50 studies providing a clear association with AD. Despite major technological advances, the same set of top-ranked pathways have been consistently related to AD for 30 years, including AD, immune system, metabolic pathways, cholinergic synapse, long-term depression, proteasome, diabetes, cancer, and chemokine signaling. AD pathways studied appear biased: animal model and human subject studies prioritize different AD pathways. Surprisingly, human genetic discoveries and drug targeting are not enriched in the most frequently studied pathways. Our findings suggest that not only is this disorder incredibly complex, but that its functional reach is also nearly global. As a consequence of our study, research results can now be assessed in the context of the wider AD literature, supporting the design of drug therapies that target a broader range of mechanisms. The results of this study can be explored at www.adpathways.org.

Region-based analysis of rare genomic variants in whole-genome sequencing datasets reveal two novel Alzheimer's disease-associated genes: DTNB and DLG2.

Alzheimer's disease (AD) is a genetically complex disease for which nearly 40 loci have now been identified via genome-wide association studies (GWAS). We attempted to identify groups of rare variants (alternate allele frequency <0.01) associated with AD in a region-based, whole-genome sequencing (WGS) association study (rvGWAS) of two independent AD family datasets (NIMH/NIA; 2247 individuals; 605 families). Employing a sliding window approach across the genome, we identified several regions that achieved association p values <10-6, using the burden test or the SKAT statistic. The genomic region around the dystobrevin beta (DTNB) gene was identified with the burden and SKAT test and replicated in case/control samples from the ADSP study reaching genome-wide significance after meta-analysis (pmeta = 4.74 × 10-8). SKAT analysis also revealed region-based association around the Discs large homolog 2 (DLG2) gene and replicated in case/control samples from the ADSP study (pmeta = 1 × 10-6). In conclusion, in a region-based rvGWAS of AD we identified two novel AD genes, DLG2 and DTNB, based on association with rare variants.

Immunotherapy for breast cancer using EpCAM aptamer tumor-targeted gene knockdown.

New strategies for cancer immunotherapy are needed since most solid tumors do not respond to current approaches. Here we used epithelial cell adhesion molecule EpCAM (a tumor-associated antigen highly expressed on common epithelial cancers and their tumor-initiating cells) aptamer-linked small-interfering RNA chimeras (AsiCs) to knock down genes selectively in EpCAM+ tumors with the goal of making cancers more visible to the immune system. Knockdown of genes that function in multiple steps of cancer immunity was evaluated in aggressive triple-negative and HER2+ orthotopic, metastatic, and genetically engineered mouse breast cancer models. Gene targets were chosen whose knockdown was predicted to promote tumor neoantigen expression (Upf2, Parp1, Apex1), phagocytosis, and antigen presentation (Cd47), reduce checkpoint inhibition (Cd274), or cause tumor cell death (Mcl1). Four of the six AsiC (Upf2, Parp1, Cd47, and Mcl1) potently inhibited tumor growth and boosted tumor-infiltrating immune cell functions. AsiC mixtures were more effective than individual AsiC and could synergize with anti-PD-1 checkpoint inhibition.

Aß-accelerated neurodegeneration caused by Alzheimer's-associated ACE variant R1279Q is rescued by angiotensin system inhibition in mice.

Recent genome-wide association studies identified the angiotensin-converting enzyme gene (ACE) as an Alzheimer's disease (AD) risk locus. However, the pathogenic mechanism by which ACE causes AD is unknown. Using whole-genome sequencing, we identified rare ACE coding variants in AD families and investigated one, ACE1 R1279Q, in knockin (KI) mice. Similar to AD, ACE1 was increased in neurons, but not microglia or astrocytes, of KI brains, which became elevated further with age. Angiotensin II (angII) and angII receptor AT1R signaling were also increased in KI brains. Autosomal dominant neurodegeneration and neuroinflammation occurred with aging in KI hippocampus, which were absent in the cortex and cerebellum. Female KI mice exhibited greater hippocampal electroencephalograph disruption and memory impairment compared to males. ACE variant effects were more pronounced in female KI mice, suggesting a mechanism for higher AD risk in women. Hippocampal neurodegeneration was completely rescued by treatment with brain-penetrant drugs that inhibit ACE1 and AT1R. Although ACE variant-induced neurodegeneration did not depend on ß-amyloid (Aß) pathology, amyloidosis in 5XFAD mice crossed to KI mice accelerated neurodegeneration and neuroinflammation, whereas Aß deposition was unchanged. KI mice had normal blood pressure and cerebrovascular functions. Our findings strongly suggest that increased ACE1/angII signaling causes aging-dependent, Aß-accelerated selective hippocampal neuron vulnerability and female susceptibility, hallmarks of AD that have hitherto been enigmatic. We conclude that repurposed brain-penetrant ACE inhibitors and AT1R blockers may protect against AD.

Silencing of the Drosophila ortholog of SOX5 leads to abnormal neuronal development and behavioral impairment.

SOX5 encodes a transcription factor that is expressed in multiple tissues including heart, lung and brain. Mutations in SOX5 have been previously found in patients with amyotrophic lateral sclerosis (ALS) and developmental delay, intellectual disability and dysmorphic features. To characterize the neuronal role of SOX5, we silenced the Drosophila ortholog of SOX5, Sox102F, by RNAi in various neuronal subtypes in Drosophila. Silencing of Sox102F led to misorientated and disorganized michrochaetes, neurons with shorter dendritic arborization (DA) and reduced complexity, diminished larval peristaltic contractions, loss of neuromuscular junction bouton structures, impaired olfactory perception, and severe neurodegeneration in brain. Silencing of SOX5 in human SH-SY5Y neuroblastoma cells resulted in a significant repression of WNT signaling activity and altered expression of WNT-related genes. Genetic association and meta-analyses of the results in several large family-based and case-control late-onset familial Alzheimer's disease (LOAD) samples of SOX5 variants revealed several variants that show significant association with AD disease status. In addition, analysis for rare and highly penetrate functional variants revealed four novel variants/mutations in SOX5, which taken together with functional prediction analysis, suggests a strong role of SOX5 causing AD in the carrier families. Collectively, these findings indicate that SOX5 is a novel candidate gene for LOAD with an important role in neuronal function. The genetic findings warrant further studies to identify and characterize SOX5 variants that confer risk for AD, ALS and intellectual disability.

Iterative sorting reveals CD133+ and CD133- melanoma cells as phenotypically distinct populations.

BACKGROUND: The heterogeneity and tumourigenicity of metastatic melanoma is attributed to a cancer stem cell model, with CD133 considered to be a cancer stem cell marker in melanoma as well as other tumours, but its role has remained controversial. METHODS: We iteratively sorted CD133+ and CD133- cells from 3 metastatic melanoma cell lines, and observed tumourigenicity and phenotypic characteristics over 7 generations of serial xeno-transplantation in NOD/SCID mice. RESULTS: We demonstrate that iterative sorting is required to make highly pure populations of CD133+ and CD133- cells from metastatic melanoma, and that these two populations have distinct characteristics not related to the cancer stem cell phenotype. In vitro, gene set enrichment analysis indicated CD133+ cells were related to a proliferative phenotype, whereas CD133- cells were of an invasive phenotype. However, in vivo, serial transplantation of CD133+ and CD133- tumours over 7 generations showed that both populations were equally able to initiate and propagate tumours. Despite this, both populations remained phenotypically distinct, with CD133- cells only able to express CD133 in vivo and not in vitro. Loss of CD133 from the surface of a CD133+ cell was observed in vitro and in vivo, however CD133- cells derived from CD133+ retained the CD133+ phenotype, even in the presence of signals from the tumour microenvironment. CONCLUSION: We show for the first time the necessity of iterative sorting to isolate pure marker-positive and marker-negative populations for comparative studies, and present evidence that despite CD133+ and CD133- cells being equally tumourigenic, they display distinct phenotypic differences, suggesting CD133 may define a distinct lineage in melanoma.

Integrated Genomic Analysis of Diverse Induced Pluripotent Stem Cells from the Progenitor Cell Biology Consortium.

The rigorous characterization of distinct induced pluripotent stem cells (iPSC) derived from multiple reprogramming technologies, somatic sources, and donors is required to understand potential sources of variability and downstream potential. To achieve this goal, the Progenitor Cell Biology Consortium performed comprehensive experimental and genomic analyses of 58 iPSC from ten laboratories generated using a variety of reprogramming genes, vectors, and cells. Associated global molecular characterization studies identified functionally informative correlations in gene expression, DNA methylation, and/or copy-number variation among key developmental and oncogenic regulators as a result of donor, sex, line stability, reprogramming technology, and cell of origin. Furthermore, X-chromosome inactivation in PSC produced highly correlated differences in teratoma-lineage staining and regulator expression upon differentiation. All experimental results, and raw, processed, and metadata from these analyses, including powerful tools, are interactively accessible from a new online portal at https://www.synapse.org to serve as a reusable resource for the stem cell community.

Exome sequencing reveals recurrent germ line variants in patients with familial Waldenström macroglobulinemia.

Familial aggregation of Waldenström macroglobulinemia (WM) cases, and the clustering of B-cell lymphoproliferative disorders among first-degree relatives of WM patients, has been reported. Nevertheless, the possible contribution of inherited susceptibility to familial WM remains unrevealed. We performed whole exome sequencing on germ line DNA obtained from 4 family members in which coinheritance for WM was documented in 3 of them, and screened additional independent 246 cases by using gene-specific mutation sequencing. Among the shared germ line variants, LAPTM5(c403t) and HCLS1(g496a) were the most recurrent, being present in 3/3 affected members of the index family, detected in 8% of the unrelated familial cases, and present in 0.5% of the nonfamilial cases and in <0.05 of a control population. LAPTM5 and HCLS1 appeared as relevant WM candidate genes that characterized familial WM individuals and were also functionally relevant to the tumor clone. These findings highlight potentially novel contributors for the genetic predisposition to familial WM and indicate that LAPTM5(c403t) and HCLS1(g496a) may represent predisposition alleles in patients with familial WM.

Mycoplasma Infection Alters Cancer Stem Cell Properties in Vitro.

Cancer cell lines can be useful to model cancer stem cells. Infection with Mycoplasma species is an insidious problem in mammalian cell culture. While investigating stem-like properties in early passage melanoma cell lines, we noted poorly reproducible results from an aliquot of a cell line that was later found to be infected with Mycoplasma hyorhinis. Deliberate infection of other early passage melanoma cell lines aliquots induced variable and unpredictable effects on expression of putative cancer stem cell markers, clonogenicity, proliferation and global gene expression. Cell lines established in stem cell media (SCM) were equally susceptible. Mycoplasma status is rarely reported in publications using cultured cells to study the cancer stem cell hypothesis. Our work highlights the importance of surveillance for Mycoplasma infection while using any cultured cells to interrogate tumor heterogeneity.

Inferring an Evolutionary Tree of Uveal Melanoma From Genomic Copy Number Aberrations.

PURPOSE: The purpose of this study is to study the genomic evolution of primary uveal melanoma. METHODS: Primary uveal melanoma genomic DNA was assayed on the Illumina Human660W-Quad v1.0 DNA Analysis BeadChip. Raw signal intensity data were quantile normalized to estimate copy number aberration with the Genome Alteration Print algorithm. Distance between samples was calculated as the Manhattan distance between the copy number profiles of the tumors. From the distance matrix, a phylogenetic network (evolutionary relationship inference) was estimated using SplitsTree4. RESULTS: Of the 57 tumors, one (1.8%) was discarded because of a failed assay, and seven (12.3%) were revealed to be mixtures of several cell populations that could not be resolved. Three clades of tumor were identified (A [59.2%], B [32.7%], and C [6.1%]), each following a distinct evolutionary path and each associated with metastatic status (P = 0.01). One tumor (2.0%) did not fit into any clade. From a normal diploid melanocyte, a few tumors (clade C) lose a large portion of chromosome 6q, but do not develop any mutations on 8q. In an alternate path, the vast majority of tumors (clade A and clade B [91.9%]) gain a copy of the telomeric half of 8q. A majority of these tumors (clade A) subsequently lose a copy of chromosome 3, as well as gain the centromeric half of 8q. The other tumors (clade B) gain copies of 6p, as well as regions on 11p and 22q. CONCLUSIONS: Our data suggest that there is little overlap in the subtypes of uveal melanoma after divergence (identified as clades A and B) and that these distinct subtypes are not likely to crossover or transform from one major clade to another.

edgeRun: an R package for sensitive, functionally relevant differential expression discovery using an unconditional exact test.

Next-generation sequencing platforms for measuring digital expression such as RNA-Seq are displacing traditional microarray-based methods in biological experiments. The detection of differentially expressed genes between groups of biological conditions has led to the development of numerous bioinformatics tools, but so far, few exploit the expanded dynamic range afforded by the new technologies. We present edgeRun, an R package that implements an unconditional exact test that is a more powerful version of the exact test in edgeR. This increase in power is especially pronounced for experiments with as few as two replicates per condition, for genes with low total expression and with large biological coefficient of variation. In comparison with a panel of other tools, edgeRun consistently captures functionally similar differentially expressed genes.

Sequencing of captive target transcripts identifies the network of regulated genes and functions of primate-specific miR-522.

Identifying microRNA (miRNA)-regulated genes is key to understanding miRNA function. However, many miRNA recognition elements (MREs) do not follow canonical "seed" base-pairing rules, making identification of bona fide targets challenging. Here, we apply an unbiased sequencing-based systems approach to characterize miR-522, a member of the oncogenic primate-specific chromosome 19 miRNA cluster, highly expressed in poorly differentiated cancers. To identify miRNA targets, we sequenced full-length transcripts captured by a biotinylated miRNA mimic. Within these targets, mostly noncanonical MREs were identified by sequencing RNase-resistant fragments. miR-522 overexpression reduced mRNA, protein levels, and luciferase activity of >70% of a random list of candidate target genes and MREs. Bioinformatic analysis suggested that miR-522 regulates cell proliferation, detachment, migration, and epithelial-mesenchymal transition. miR-522 induces G1 cell-cycle arrest and causes cells to detach without anoikis, become invasive, and express mesenchymal genes. Thus, our method provides a simple but effective technique for identifying miRNA-regulated genes and biological function.

A network of epigenetic regulators guides developmental haematopoiesis in vivo.

The initiation of cellular programs is orchestrated by key transcription factors and chromatin regulators that activate or inhibit target gene expression. To generate a compendium of chromatin factors that establish the epigenetic code during developmental haematopoiesis, a large-scale reverse genetic screen was conducted targeting orthologues of 425 human chromatin factors in zebrafish. A set of chromatin regulators was identified that target different stages of primitive and definitive blood formation, including factors not previously implicated in haematopoiesis. We identified 15 factors that regulate development of primitive erythroid progenitors and 29 factors that regulate development of definitive haematopoietic stem and progenitor cells. These chromatin factors are associated with SWI/SNF and ISWI chromatin remodelling, SET1 methyltransferase, CBP-p300-HBO1-NuA4 acetyltransferase, HDAC-NuRD deacetylase, and Polycomb repressive complexes. Our work provides a comprehensive view of how specific chromatin factors and their associated complexes play a major role in the establishment of haematopoietic cells in vivo.

A comprehensive promoter landscape identifies a novel promoter for CD133 in restricted tissues, cancers, and stem cells.

PROM1 is the gene encoding prominin-1 or CD133, an important cell surface marker for the isolation of both normal and cancer stem cells. PROM1 transcripts initiate at a range of transcription start sites (TSS) associated with distinct tissue and cancer expression profiles. Using high resolution Cap Analysis of Gene Expression (CAGE) sequencing we characterize TSS utilization across a broad range of normal and developmental tissues. We identify a novel proximal promoter (P6) within CD133(+) melanoma cell lines and stem cells. Additional exon array sampling finds P6 to be active in populations enriched for mesenchyme, neural stem cells and within CD133(+) enriched Ewing sarcomas. The P6 promoter is enriched with respect to previously characterized PROM1 promoters for a HMGI/Y (HMGA1) family transcription factor binding site motif and exhibits different epigenetic modifications relative to the canonical promoter region of PROM1.

A genome-wide siRNA screen identifies proteasome addiction as a vulnerability of basal-like triple-negative breast cancer cells.

Basal-like triple-negative breast cancers (TNBCs) have poor prognosis. To identify basal-like TNBC dependencies, a genome-wide siRNA lethality screen compared two human breast epithelial cell lines transformed with the same genes: basal-like BPLER and myoepithelial HMLER. Expression of the screen's 154 BPLER dependency genes correlated with poor prognosis in breast, but not lung or colon, cancer. Proteasome genes were overrepresented hits. Basal-like TNBC lines were selectively sensitive to proteasome inhibitor drugs relative to normal epithelial, luminal, and mesenchymal TNBC lines. Proteasome inhibition reduced growth of established basal-like TNBC tumors in mice and blocked tumor-initiating cell function and macrometastasis. Proteasome addiction in basal-like TNBCs was mediated by NOXA and linked to MCL-1 dependence.

Pathprinting: An integrative approach to understand the functional basis of disease.

New strategies to combat complex human disease require systems approaches to biology that integrate experiments from cell lines, primary tissues and model organisms. We have developed Pathprint, a functional approach that compares gene expression profiles in a set of pathways, networks and transcriptionally regulated targets. It can be applied universally to gene expression profiles across species. Integration of large-scale profiling methods and curation of the public repository overcomes platform, species and batch effects to yield a standard measure of functional distance between experiments. We show that pathprints combine mouse and human blood developmental lineage, and can be used to identify new prognostic indicators in acute myeloid leukemia. The code and resources are available at http://compbio.sph.harvard.edu/hidelab/pathprint.

Comparison of glioma stem cells to neural stem cells from the adult human brain identifies dysregulated Wnt- signaling and a fingerprint associated with clinical outcome.

Glioblastoma is the most common brain tumor. Median survival in unselected patients is <10 months. The tumor harbors stem-like cells that self-renew and propagate upon serial transplantation in mice, although the clinical relevance of these cells has not been well documented. We have performed the first genome-wide analysis that directly relates the gene expression profile of nine enriched populations of glioblastoma stem cells (GSCs) to five identically isolated and cultivated populations of stem cells from the normal adult human brain. Although the two cell types share common stem- and lineage-related markers, GSCs show a more heterogeneous gene expression. We identified a number of pathways that are dysregulated in GSCs. A subset of these pathways has previously been identified in leukemic stem cells, suggesting that cancer stem cells of different origin may have common features. Genes upregulated in GSCs were also highly expressed in embryonic and induced pluripotent stem cells. We found that canonical Wnt-signaling plays an important role in GSCs, but not in adult human neural stem cells. As well we identified a 30-gene signature highly overexpressed in GSCs. The expression of these signature genes correlates with clinical outcome and demonstrates the clinical relevance of GSCs.

Integrating murine gene expression studies to understand obstructive lung disease due to chronic inhaled endotoxin.

RATIONALE: Endotoxin is a near ubiquitous environmental exposure that that has been associated with both asthma and chronic obstructive pulmonary disease (COPD). These obstructive lung diseases have a complex pathophysiology, making them difficult to study comprehensively in the context of endotoxin. Genome-wide gene expression studies have been used to identify a molecular snapshot of the response to environmental exposures. Identification of differentially expressed genes shared across all published murine models of chronic inhaled endotoxin will provide insight into the biology underlying endotoxin-associated lung disease. METHODS: We identified three published murine models with gene expression profiling after repeated low-dose inhaled endotoxin. All array data from these experiments were re-analyzed, annotated consistently, and tested for shared genes found to be differentially expressed. Additional functional comparison was conducted by testing for significant enrichment of differentially expressed genes in known pathways. The importance of this gene signature in smoking-related lung disease was assessed using hierarchical clustering in an independent experiment where mice were exposed to endotoxin, smoke, and endotoxin plus smoke. RESULTS: A 101-gene signature was detected in three murine models, more than expected by chance. The three model systems exhibit additional similarity beyond shared genes when compared at the pathway level, with increasing enrichment of inflammatory pathways associated with longer duration of endotoxin exposure. Genes and pathways important in both asthma and COPD were shared across all endotoxin models. Mice exposed to endotoxin, smoke, and smoke plus endotoxin were accurately classified with the endotoxin gene signature. CONCLUSIONS: Despite the differences in laboratory, duration of exposure, and strain of mouse used in three experimental models of chronic inhaled endotoxin, surprising similarities in gene expression were observed. The endotoxin component of tobacco smoke may play an important role in disease development.

Population differences in transcript-regulator expression quantitative trait loci.

Gene expression quantitative trait loci (eQTL) are useful for identifying single nucleotide polymorphisms (SNPs) associated with diseases. At times, a genetic variant may be associated with a master regulator involved in the manifestation of a disease. The downstream target genes of the master regulator are typically co-expressed and share biological function. Therefore, it is practical to screen for eQTLs by identifying SNPs associated with the targets of a transcript-regulator (TR). We used a multivariate regression with the gene expression of known targets of TRs and SNPs to identify TReQTLs in European (CEU) and African (YRI) HapMap populations. A nominal p-value of <1×10(-6) revealed 234 SNPs in CEU and 154 in YRI as TReQTLs. These represent 36 independent (tag) SNPs in CEU and 39 in YRI affecting the downstream targets of 25 and 36 TRs respectively. At a false discovery rate (FDR) = 45%, one cis-acting tag SNP (within 1 kb of a gene) in each population was identified as a TReQTL. In CEU, the SNP (rs16858621) in Pcnxl2 was found to be associated with the genes regulated by CREM whereas in YRI, the SNP (rs16909324) was linked to the targets of miRNA hsa-miR-125a. To infer the pathways that regulate expression, we ranked TReQTLs by connectivity within the structure of biological process subtrees. One TReQTL SNP (rs3790904) in CEU maps to Lphn2 and is associated (nominal p-value = 8.1×10(-7)) with the targets of the X-linked breast cancer suppressor Foxp3. The structure of the biological process subtree and a gene interaction network of the TReQTL revealed that tumor necrosis factor, NF-kappaB and variants in G-protein coupled receptors signaling may play a central role as communicators in Foxp3 functional regulation. The potential pleiotropic effect of the Foxp3 TReQTLs was gleaned from integrating mRNA-Seq data and SNP-set enrichment into the analysis.

Chronic endotoxin exposure produces airflow obstruction and lung dendritic cell expansion.

Little is known about the mechanisms of persistent airflow obstruction that result from chronic occupational endotoxin exposure. We sought to analyze the inflammatory response underlying persistent airflow obstruction as a result of chronic occupational endotoxin exposure. We developed a murine model of daily inhaled endotoxin for periods of 5 days to 8 weeks. We analyzed physiologic lung dysfunction, lung histology, bronchoalveolar lavage fluid and total lung homogenate inflammatory cell and cytokine profiles, and pulmonary gene expression profiles. We observed an increase in airway hyperresponsiveness as a result of chronic endotoxin exposure. After 8 weeks, the mice exhibited an increase in bronchoalveolar lavage and lung neutrophils that correlated with an increase in proinflammatory cytokines. Detailed analyses of inflammatory cell subsets revealed an expansion of dendritic cells (DCs), and in particular, proinflammatory DCs, with a reduced percentage of macrophages. Gene expression profiling revealed the up-regulation of a panel of genes that was consistent with DC recruitment, and lung histology revealed an accumulation of DCs in inflammatory aggregates around the airways in 8-week-exposed animals. Repeated, low-dose LPS inhalation, which mirrors occupational exposure, resulted in airway hyperresponsiveness, associated with a failure to resolve the proinflammatory response, an inverted macrophage to DC ratio, and a significant rise in the inflammatory DC population. These findings point to a novel underlying mechanism of airflow obstruction as a result of occupational LPS exposure, and suggest molecular and cellular targets for therapeutic development.