Sex differences in brain protein expression and disease.

Most complex human traits differ by sex, but we have limited insight into the underlying mechanisms. Here, we investigated the influence of biological sex on protein expression and its genetic regulation in 1,277 human brain proteomes. We found that 13.2% (1,354) of brain proteins had sex-differentiated abundance and 1.5% (150) of proteins had sex-biased protein quantitative trait loci (sb-pQTLs). Among genes with sex-biased expression, we found 67% concordance between sex-differentiated protein and transcript levels; however, sex effects on the genetic regulation of expression were more evident at the protein level. Considering 24 psychiatric, neurologic and brain morphologic traits, we found that an average of 25% of their putatively causal genes had sex-differentiated protein abundance and 12 putatively causal proteins had sb-pQTLs. Furthermore, integrating sex-specific pQTLs with sex-stratified genome-wide association studies of six psychiatric and neurologic conditions, we uncovered another 23 proteins contributing to these traits in one sex but not the other. Together, these findings begin to provide insights into mechanisms underlying sex differences in brain protein expression and disease.

Role of TET1-mediated epigenetic modulation in Alzheimer's disease.

Alzheimer's disease (AD) is a neurodegenerative disorder influenced by a complex interplay of environmental, epigenetic, and genetic factors. DNA methylation (5mC) and hydroxymethylation (5hmC) are DNA modifications that serve as tissue-specific and temporal regulators of gene expression. TET family enzymes dynamically regulate these epigenetic modifications in response to environmental conditions, connecting environmental factors with gene expression. Previous epigenetic studies have identified 5mC and 5hmC changes associated with AD. In this study, we performed targeted resequencing of TET1 on a cohort of early-onset AD (EOAD) and control samples. Through gene-wise burden analysis, we observed significant enrichment of rare TET1 variants associated with AD (p = 0.04). We also profiled 5hmC in human postmortem brain tissues from AD and control groups. Our analysis identified differentially hydroxymethylated regions (DhMRs) in key genes responsible for regulating the methylome: TET3, DNMT3L, DNMT3A, and MECP2. To further investigate the role of Tet1 in AD pathogenesis, we used the 5xFAD mouse model with a Tet1 KO allele to examine how Tet1 loss influences AD pathogenesis. We observed significant changes in neuropathology, 5hmC, and RNA expression associated with Tet1 loss, while the behavioral alterations were not significant. The loss of Tet1 significantly increased amyloid plaque burden in the 5xFAD mouse (p = 0.044) and lead to a non-significant trend towards exacerbated AD-associated stress response in 5xFAD mice. At the molecular level, we found significant DhMRs enriched in genes involved in pathways responsible for neuronal projection organization, dendritic spine development and organization, and myelin assembly. RNA-Seq analysis revealed a significant increase in the expression of AD-associated genes such as Mpeg1, Ctsd, and Trem2. In conclusion, our results suggest that TET enzymes, particularly TET1, which regulate the methylome, may contribute to AD pathogenesis, as the loss of TET function increases AD-associated pathology.

Shared mechanisms across the major psychiatric and neurodegenerative diseases.

Several common psychiatric and neurodegenerative diseases share epidemiologic risk; however, whether they share pathophysiology is unclear and is the focus of our investigation. Using 25 GWAS results and LD score regression, we find eight significant genetic correlations between psychiatric and neurodegenerative diseases. We integrate the GWAS results with human brain transcriptomes (n = 888) and proteomes (n = 722) to identify cis- and trans- transcripts and proteins that are consistent with a pleiotropic or causal role in each disease, referred to as causal proteins for brevity. Within each disease group, we find many distinct and shared causal proteins. Remarkably, 30% (13 of 42) of the neurodegenerative disease causal proteins are shared with psychiatric disorders. Furthermore, we find 2.6-fold more protein-protein interactions among the psychiatric and neurodegenerative causal proteins than expected by chance. Together, our findings suggest these psychiatric and neurodegenerative diseases have shared genetic and molecular pathophysiology, which has important ramifications for early treatment and therapeutic development.

LDL cholesterol is associated with higher AD neuropathology burden independent of APOE.

OBJECTIVE: APOE is a strong risk factor for Alzheimer's disease (AD) and associated with higher low-density lipoprotein cholesterol (LDL-C) levels. Moreover, LDL-C is associated with the development of clinically ascertained AD; however, whether this association is present with the underlying neuropathological manifestations of AD or whether it is independent of the effect of APOE is unknown and is the focus of this paper. METHODS: Individuals in the Religious Orders Study/Memory and Ageing Project cohorts with longitudinal measures of blood lipids and detailed autopsies were studied. We modelled the relationship between blood lipids and 12 age-related brain pathologies using a linear mixed model adjusted for potential confounding factors and stratified by APOE genotype with overall significance determined by meta-analysis. Blood lipids considered were LDL-C, high-density lipoprotein cholesterol and triglycerides. Brain pathologies included AD pathology measured by silver staining (Braak stage, a modified Consortium to Establish a Registry for Alzheimer's Disease [CERAD] score and global AD pathology) and immunohistochemistry (beta-amyloid and neurofibrillary tangles) as well as cerebral microinfarct, cerebral macroinfarct, cerebral amyloid angiopathy, cerebral atherosclerosis, hippocampal sclerosis, TDP-43 cytoplasmic inclusions and Lewy bodies. RESULTS: 559 participants (69.1% female) had complete data for analysis. They were followed for a median of 7 years and a median of 3 years prior to dementia onset. LDL-C was associated with all measures of AD neuropathology (neurofibrillary tangles, beta-amyloid, Braak stage, modified CERAD score and global AD pathology) and cerebral amyloid angiopathy independent of APOE after adjusting for age, sex, cholesterol-lowering medication use, body mass index, smoking and education at false discovery rate (FDR) p-value <0.05. CONCLUSIONS: These findings implicate LDL-C in the pathophysiology of AD independent of APOE and suggest LDL-C is a modifiable risk factor for AD.

Integrating human brain proteomes with genome-wide association data implicates novel proteins in post-traumatic stress disorder.

Genome-wide association studies (GWAS) have identified several risk loci for post-traumatic stress disorder (PTSD); however, how they confer PTSD risk remains unclear. We aimed to identify genes that confer PTSD risk through their effects on brain protein abundance to provide new insights into PTSD pathogenesis. To that end, we integrated human brain proteomes with PTSD GWAS results to perform a proteome-wide association study (PWAS) of PTSD, followed by Mendelian randomization, using a discovery and confirmatory study design. Brain proteomes (N = 525) were profiled from the dorsolateral prefrontal cortex using mass spectrometry. The Million Veteran Program (MVP) PTSD GWAS (n = 186,689) was used for the discovery PWAS, and the Psychiatric Genomics Consortium PTSD GWAS (n = 174,659) was used for the confirmatory PWAS. To understand whether genes identified at the protein-level were also evident at the transcript-level, we performed a transcriptome-wide association study (TWAS) using human brain transcriptomes (N = 888) and the MVP PTSD GWAS results. We identified 11 genes that contribute to PTSD pathogenesis via their respective cis-regulated brain protein abundance. Seven of 11 genes (64%) replicated in the confirmatory PWAS and 4 of 11 also had their cis-regulated brain mRNA levels associated with PTSD. High confidence level was assigned to 9 of 11 genes after considering evidence from the confirmatory PWAS and TWAS. Most of the identified genes are expressed in other PTSD-relevant brain regions and several are preferentially expressed in excitatory neurons, astrocytes, and oligodendrocyte precursor cells. These genes are novel, promising targets for mechanistic and therapeutic studies to find new treatments for PTSD.

A Genetic Study of Cerebral Atherosclerosis Reveals Novel Associations with NTNG1 and CNOT3.

Cerebral atherosclerosis is a leading cause of stroke and an important contributor to dementia. Yet little is known about its genetic basis. To examine the association of common single nucleotide polymorphisms with cerebral atherosclerosis severity, we conducted a genomewide association study (GWAS) using data collected as part of two community-based cohort studies in the United States, the Religious Orders Study (ROS) and Rush Memory and Aging Project (MAP). Both studies enroll older individuals and exclude participants with signs of dementia at baseline. From our analysis of 1325 participants of European ancestry who had genotype and neuropathologically assessed cerebral atherosclerosis measures available, we found a novel locus for cerebral atherosclerosis in NTNG1. The locus comprises eight SNPs, including two independent significant SNPs: rs6664221 (ß = -0.27, 95% CI = (-0.35, -0.19), p = 1.29 × 10-10) and rs10881463 (ß = -0.20, 95% CI = (-0.27, -0.13), p = 3.40 × 10-8). We further found that the SNPs may influence cerebral atherosclerosis by regulating brain protein expression of CNOT3. CNOT3 is a subunit of CCR4-NOT, which has been shown to be a master regulator of mRNA stability and translation and an important complex for cholesterol homeostasis. In summary, we identify a novel genetic locus for cerebral atherosclerosis and a potential mechanism linking this variation to cerebral atherosclerosis progression. These findings offer insights into the genetic effects on cerebral atherosclerosis.

Important Correlates of Purpose in Life Identified Through a Machine Learning Approach.

OBJECTIVE: A wealth of evidence has linked purpose in life (PiL) to better mental and physical health and healthy aging. Here, the authors aimed to determine important correlates of PiL using a machine learning approach. METHODS: Participants were recruited from retirement communities by the Rush Memory and Aging Project and assessed for childhood experience, adulthood sociodemographic factors (e.g., education, income, marital status), lifestyle and health behavior (e.g., cognitively stimulating activities, exercise, social activities, social network size), psychological factors (e.g., depression, loneliness, perceived discrimination, perceived social support), personality traits (e.g., PiL, harm avoidance), and medical conditions. Elastic Net was implemented to identify important correlates of PiL. RESULTS: A total of 1,839 participants were included in our analysis. Among the 23 variables provided to Elastic Net, 10 were identified as important correlates of PiL. In order of decreasing effect size, factors associated with lower PiL were loneliness, harm avoidance, older age, and depressive symptoms, while those associated with greater PiL were perceived social support, more social activities, more years of education, higher income, intact late-life cognitive performance, and more middle-age cognitive activities. CONCLUSION: Our findings identify potentially important modifiable factors as targets for intervention strategies to enhance PiL.

DNA isolation protocol effects on nuclear DNA analysis by microarrays, droplet digital PCR, and whole genome sequencing, and on mitochondrial DNA copy number estimation.

Potential bias introduced during DNA isolation is inadequately explored, although it could have significant impact on downstream analysis. To investigate this in human brain, we isolated DNA from cerebellum and frontal cortex using spin columns under different conditions, and salting-out. We first analysed DNA using array CGH, which revealed a striking wave pattern suggesting primarily GC-rich cerebellar losses, even against matched frontal cortex DNA, with a similar pattern on a SNP array. The aCGH changes varied with the isolation protocol. Droplet digital PCR of two genes also showed protocol-dependent losses. Whole genome sequencing showed GC-dependent variation in coverage with spin column isolation from cerebellum. We also extracted and sequenced DNA from substantia nigra using salting-out and phenol / chloroform. The mtDNA copy number, assessed by reads mapping to the mitochondrial genome, was higher in substantia nigra when using phenol / chloroform. We thus provide evidence for significant method-dependent bias in DNA isolation from human brain, as reported in rat tissues. This may contribute to array "waves", and could affect copy number determination, particularly if mosaicism is being sought, and sequencing coverage. Variations in isolation protocol may also affect apparent mtDNA abundance.

Role of CDKN2C Copy Number in Sporadic Medullary Thyroid Carcinoma.

BACKGROUND: The cyclin-dependent-kinase inhibitors (CDKN)/retinoblastoma (RB1) pathway has been implicated as having a role in medullary thyroid carcinoma (MTC) tumorigenesis. CDKN2C loss has been associated with RET-mediated MTC in humans but with minimal phenotypic correlation provided. The objective of this study was to evaluate the association between tumor RET mutation status, CDKN2C loss, and aggressiveness of MTC in a cohort of patients with sporadic disease. METHODS: Tumors from patients with sporadic MTC treated at a single institution were evaluated for somatic RETM918T mutation and CDKN2C copy number loss. These variables were compared to patient demographics, pathology detail, clinical course, and disease-specific and overall survival. RESULTS: Sixty-two MTC cases with an initial surgery date ranging from 1983 to 2009 met the inclusion criteria, of whom 36 (58%) were male. The median age at initial surgery was 53 years (range 22-81 years). The median tumor size was 30 mm (range 6-145 mm) with 29 (57%) possessing extrathyroidal extension. Nodal and/or distant metastasis at presentation was found in 47/60 (78%) and 12/61 (20%) patients, respectively. Median follow-up time was 10.5 years (range 1.1-27.8 years) for the censored observations. The presence of CDKN2C loss was associated with worse M stage and overall AJCC stage. Median overall survival of patients with versus without CDKN2C loss was 4.14 [confidence interval (CI) 1.93-NA] versus 18.27 [CI 17.24-NA] years (p < 0.0001). Median overall survival of patients with a combined somatic RETM918T mutation and CDKN2C loss versus no somatic RETM918T mutation and CDKN2C loss versus somatic RETM918T mutation and CDKN2C 2N versus no somatic RETM918T mutation and CDKN2C 2N was 2.38 [CI 1.67-NA] years versus 10.81 [CI 2.46-NA] versus 17.24 [CI 9.82-NA] versus not reached [CI 13.46-NA] years (p < 0.0001). CONCLUSIONS: The detection of somatic CDKN2C loss is associated with the presence of distant metastasis at presentation as well decreased overall survival, a relationship enhanced by concomitant RETM918T mutation. Further defining the genes involved in the progression of metastatic MTC will be an important step toward identifying pathways of disease progression and new therapeutic targets.

Rapid and powerful detection of subtle allelic imbalance from exome sequencing data with hapLOHseq.

MOTIVATION: The detection of subtle genomic allelic imbalance events has many potential applications. For example, identifying cancer-associated allelic imbalanced regions in low tumor-cellularity samples or in low-proportion tumor subclones can be used for early cancer detection, prognostic assessment and therapeutic selection in cancer patients. We developed hapLOHseq for the detection of subtle allelic imbalance events from next-generation sequencing data. RESULTS: Our method identified events of 10 megabases or greater occurring in as little as 16% of the sample in exome sequencing data (at 80×) and 4% in whole genome sequencing data (at 30×), far exceeding the capabilities of existing software. We also found hapLOHseq to be superior at detecting large chromosomal changes across a series of pancreatic samples from TCGA. AVAILABILITY AND IMPLEMENTATION: hapLOHseq is available at scheet.org/software, distributed under an open source MIT license. CONTACT: pscheet@alum.wustl.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genomic Landscape Established by Allelic Imbalance in the Cancerization Field of a Normal Appearing Airway.

Visually normal cells adjacent to, and extending from, tumors of the lung may carry molecular alterations characteristics of the tumor itself, an effect referred to as airway field of cancerization. This airway field has been postulated as a model for early events in lung cancer pathogenesis. Yet the genomic landscape of somatically acquired molecular alterations in airway epithelia of lung cancer patients has remained unknown. To begin to fill this void, we sought to comprehensively characterize the genomic architecture of chromosomal alterations inducing allelic imbalance (AI) in the airway field of the most common type of lung tumors, non-small cell lung cancer (NSCLC). To do so, we conducted a genome-wide survey of multiple spatially distributed normal-appearing airways, multiregion tumor specimens, and uninvolved normal tissues or blood from 45 patients with early-stage NSCLC. We detected alterations in airway epithelia from 22 patients, with an increased frequency in NSCLCs of squamous histology. Our data also indicated a spatial gradient of AI in samples at closer proximity to the NSCLC. Chromosome 9 displayed the highest levels of AI and comprised recurrent independent events. Furthermore, the airway field AI included oncogenic gains and tumor suppressor losses in known NSCLC drivers. Our results demonstrate that genome-wide AI is common in the airway field of cancerization, providing insights into early events in the pathogenesis of NSCLC that may comprise targets for early treatment and chemoprevention. Cancer Res; 76(13); 3676-83. ©2016 AACR.

Extensive Hidden Genomic Mosaicism Revealed in Normal Tissue.

Genomic mosaicism arising from post-zygotic mutation has recently been demonstrated to occur in normal tissue of individuals ascertained with varied phenotypes, indicating that detectable mosaicism may be less an exception than a rule in the general population. A challenge to comprehensive cataloging of mosaic mutations and their consequences is the presence of heterogeneous mixtures of cells, rendering low-frequency clones difficult to discern. Here we applied a computational method using estimated haplotypes to characterize mosaic megabase-scale structural mutations in 31,100 GWA study subjects. We provide in silico validation of 293 previously identified somatic mutations and identify an additional 794 novel mutations, most of which exist at lower aberrant cell fractions than have been demonstrated in previous surveys. These mutations occurred across the genome but in a nonrandom manner, and several chromosomes and loci showed unusual levels of mutation. Our analysis supports recent findings about the relationship between clonal mosaicism and old age. Finally, our results, in which we demonstrate a nearly 3-fold higher rate of clonal mosaicism, suggest that SNP-based population surveys of mosaic structural mutations should be conducted with haplotypes for optimal discovery.

Small Amounts of Archaic Admixture Provide Big Insights into Human History.

Modern humans overlapped in time and space with other hominins, such as Neanderthals and Denisovans, and limited amounts of hybridization occurred. Here, we review recent work that has identified archaic hominin sequence that survives in modern human genomes and what these genomic excavations reveal about human evolutionary history.

Clonal evolution in breast cancer revealed by single nucleus genome sequencing.

Sequencing studies of breast tumour cohorts have identified many prevalent mutations, but provide limited insight into the genomic diversity within tumours. Here we developed a whole-genome and exome single cell sequencing approach called nuc-seq that uses G2/M nuclei to achieve 91% mean coverage breadth. We applied this method to sequence single normal and tumour nuclei from an oestrogen-receptor-positive (ER(+)) breast cancer and a triple-negative ductal carcinoma. In parallel, we performed single nuclei copy number profiling. Our data show that aneuploid rearrangements occurred early in tumour evolution and remained highly stable as the tumour masses clonally expanded. In contrast, point mutations evolved gradually, generating extensive clonal diversity. Using targeted single-molecule sequencing, many of the diverse mutations were shown to occur at low frequencies (<10%) in the tumour mass. Using mathematical modelling we found that the triple-negative tumour cells had an increased mutation rate (13.3×), whereas the ER(+) tumour cells did not. These findings have important implications for the diagnosis, therapeutic treatment and evolution of chemoresistance in breast cancer.

Identification of allelic imbalance with a statistical model for subtle genomic mosaicism.

Genetic heterogeneity in a mixed sample of tumor and normal DNA can confound characterization of the tumor genome. Numerous computational methods have been proposed to detect aberrations in DNA samples from tumor and normal tissue mixtures. Most of these require tumor purities to be at least 10-15%. Here, we present a statistical model to capture information, contained in the individual's germline haplotypes, about expected patterns in the B allele frequencies from SNP microarrays while fully modeling their magnitude, the first such model for SNP microarray data. Our model consists of a pair of hidden Markov models--one for the germline and one for the tumor genome--which, conditional on the observed array data and patterns of population haplotype variation, have a dependence structure induced by the relative imbalance of an individual's inherited haplotypes. Together, these hidden Markov models offer a powerful approach for dealing with mixtures of DNA where the main component represents the germline, thus suggesting natural applications for the characterization of primary clones when stromal contamination is extremely high, and for identifying lesions in rare subclones of a tumor when tumor purity is sufficient to characterize the primary lesions. Our joint model for germline haplotypes and acquired DNA aberration is flexible, allowing a large number of chromosomal alterations, including balanced and imbalanced losses and gains, copy-neutral loss-of-heterozygosity (LOH) and tetraploidy. We found our model (which we term J-LOH) to be superior for localizing rare aberrations in a simulated 3% mixture sample. More generally, our model provides a framework for full integration of the germline and tumor genomes to deal more effectively with missing or uncertain features, and thus extract maximal information from difficult scenarios where existing methods fail.

Haplotype-based profiling of subtle allelic imbalance with SNP arrays.

Due to limitations of surgical dissection and tumor heterogeneity, tumor samples collected for cancer genomics studies are often heavily diluted with normal tissue or contain subpopulations of cells harboring important aberrations. Methods for profiling tumor-associated allelic imbalance in such scenarios break down at aberrant cell proportions of 10%-15% and below. Here, we present an approach that offers a vast improvement for detection of subtle allelic imbalance, or low proportions of cells harboring aberrant allelic ratio among nonaberrant cells, in unpaired tumor samples using SNP microarrays. We leverage the expected pattern of allele-specific intensity ratios determined by an individual's germline haplotypes, information that has been ignored in existing approaches. We demonstrate our method on real and simulated data from the CRL-2324 breast cancer cell line genotyped on the Illumina 370K array. Assuming a 5 million SNP array, we can detect the presence of aberrant cells in proportions lower than 0.25% in the breast cancer sample, approaching the sensitivity of some minimal residual disease assays. Further, we apply a hidden Markov model to identify copy-neutral LOH (loss of heterozygosity) events as short as 11 Mb in mixtures of only 4% tumor using 370K data. We anticipate our approach will offer a new paradigm for genomic profiling of heterogeneous samples.

Rare versus common variants in pharmacogenetics: SLCO1B1 variation and methotrexate disposition.

Methotrexate is used to treat autoimmune diseases and malignancies, including acute lymphoblastic leukemia (ALL). Inter-individual variation in clearance of methotrexate results in heterogeneous systemic exposure, clinical efficacy, and toxicity. In a genome-wide association study of children with ALL, we identified SLCO1B1 as harboring multiple common polymorphisms associated with methotrexate clearance. The extent of influence of rare versus common variants on pharmacogenomic phenotypes remains largely unexplored. We tested the hypothesis that rare variants in SLCO1B1 could affect methotrexate clearance and compared the influence of common versus rare variants in addition to clinical covariates on clearance. From deep resequencing of SLCO1B1 exons in 699 children, we identified 93 SNPs, 15 of which were non-synonymous (NS). Three of these NS SNPs were common, with a minor allele frequency (MAF) >5%, one had low frequency (MAF 1%-5%), and 11 were rare (MAF <1%). NS SNPs (common or rare) predicted to be functionally damaging were more likely to be found among patients with the lowest methotrexate clearance than patients with high clearance. We verified lower function in vitro of four SLCO1B1 haplotypes that were associated with reduced methotrexate clearance. In a multivariate stepwise regression analysis adjusting for other genetic and non-genetic covariates, SLCO1B1 variants accounted for 10.7% of the population variability in clearance. Of that variability, common NS variants accounted for the majority, but rare damaging NS variants constituted 17.8% of SLCO1B1's effects (1.9% of total variation) and had larger effect sizes than common NS variants. Our results show that rare variants are likely to have an important effect on pharmacogenetic phenotypes.

Genome-wide association study identifies novel loci predisposing to cutaneous melanoma.

We performed a multistage genome-wide association study of melanoma. In a discovery cohort of 1804 melanoma cases and 1026 controls, we identified loci at chromosomes 15q13.1 (HERC2/OCA2 region) and 16q24.3 (MC1R) regions that reached genome-wide significance within this study and also found strong evidence for genetic effects on susceptibility to melanoma from markers on chromosome 9p21.3 in the p16/ARF region and on chromosome 1q21.3 (ARNT/LASS2/ANXA9 region). The most significant single-nucleotide polymorphisms (SNPs) in the 15q13.1 locus (rs1129038 and rs12913832) lie within a genomic region that has profound effects on eye and skin color; notably, 50% of variability in eye color is associated with variation in the SNP rs12913832. Because eye and skin colors vary across European populations, we further evaluated the associations of the significant SNPs after carefully adjusting for European substructure. We also evaluated the top 10 most significant SNPs by using data from three other genome-wide scans. Additional in silico data provided replication of the findings from the most significant region on chromosome 1q21.3 rs7412746 (P = 6 × 10(-10)). Together, these data identified several candidate genes for additional studies to identify causal variants predisposing to increased risk for developing melanoma.