Search

Standard single-cell (sc) proteomics of disease states inferred from multicellular organs or organoids cannot currently be related to single-cell physiology. Here, a scPatch-Clamp/Proteomics platform is developed on single neurons generated from hiPSCs bearing an Alzheimer's disease (AD) genetic mutation and compares them to isogenic wild-type controls. This approach provides both current and voltage electrophysiological data plus detailed proteomics information on single-cells. With this new method, the authors are able to observe hyperelectrical activity in the AD hiPSC-neurons, similar to that observed in the human AD brain, and correlate it to ≈1400 proteins detected at the single neuron level. Using linear regression and mediation analyses to explore the relationship between the abundance of individual proteins and the neuron's mutational and electrophysiological status, this approach yields new information on therapeutic targets in excitatory neurons not attainable by traditional methods. This combined patch-proteomics technique creates a new proteogenetic-therapeutic strategy to correlate genotypic alterations to physiology with protein expression in single-cells.

For the past few years, researchers in the Human Pangenome Reference Consortium (HPRC) have been working to catalog almost all human genomic diversity. Frazer and Schork preview an article recently published in Nature, "A draft human pangenome reference,"1 which represents the initial release of 47 fully phased diploid assemblies of genomes of individuals with diverse ancestries.

Background and aimsNAFLD is the most common chronic liver disease in children. Large pediatric studies identifying single nucleotide polymorphisms (SNPs) associated with risk and histologic severity of NAFLD are limited. Study aims included investigating SNPs associated with risk for NAFLD using family trios and association of candidate alleles with histologic severity.Approach and resultsChildren with biopsy-confirmed NAFLD were enrolled from the NASH Clinical Research Network. The Expert Pathology Committee reviewed liver histology. Genotyping was conducted with allele-specific primers for 60 candidate SNPs. Parents were enrolled for trio analysis. To assess risk for NAFLD, the transmission disequilibrium test was conducted in trios. Among cases, regression analysis assessed associations with histologic severity. A total of 822 children with NAFLD had mean age 13.2 years (SD 2.7) and mean ALT 101 U/L (SD 90). PNPLA3 (rs738409) demonstrated the strongest risk ( p = 2.24 × 10 -14 ) for NAFLD. Among children with NAFLD, stratifying by PNPLA3 s738409 genotype, the variant genotype associated with steatosis ( p = 0.005), lobular ( p = 0.03) and portal inflammation ( p = 0.002). Steatosis grade associated with TM6SF2 ( p = 0.0009), GCKR ( p = 0.0032), PNPLA3 rs738409 ( p = 0.0053), and MTTP ( p = 0.0051). Fibrosis stage associated with PARVB rs6006473 ( p = 0.0001), NR1I2 ( p = 0.0021), ADIPOR2 ( p = 0.0038), and OXTR ( p = 0.0065). PNPLA3 rs738409 ( p = 0.0002) associated with borderline zone 1 NASH.ConclusionsThis study demonstrated disease-associated SNPs in children with NAFLD. In particular, rs6006473 was highly associated with severity of fibrosis. These hypothesis-generating results support future mechanistic studies of development of adverse outcomes such as fibrosis and generation of therapeutic targets for NAFLD in children.

BackgroundAPOE is the largest genetic risk factor for Alzheimer's disease (AD), but there is a substantial polygenic component. Polygenic risk scores (PRS) can summarize small effects across the genome but may obscure differential risk across molecular processes and pathways that contribute to heterogeneity of disease presentation.ObjectiveWe examined polygenic risk impacting specific AD-associated pathways and its relationship with clinical status and biomarkers of amyloid, tau, and neurodegeneration (A/T/N).MethodsWe analyzed data from 1,411 participants from the Alzheimer's Disease Neuroimaging Initiative (ADNI). We applied pathway analysis and clustering to identify AD-associated "pathway clusters" and construct pathway-specific PRSs (excluding the APOE region). We tested associations with diagnostic status, abnormal levels of amyloid and ptau, and hippocampal volume.ResultsThirteen pathway clusters were identified, and eight pathway-specific PRSs were significantly associated with AD diagnosis. Amyloid-positivity was associated with endocytosis and fibril formation, response misfolded protein, and regulation protein tyrosine PRSs. Ptau positivity and hippocampal volume were both related to protein localization and mitophagy PRS, and ptau-positivity was also associated with an immune signaling PRS. A global AD PRS showed stronger associations with diagnosis and all biomarkers compared to pathway PRSs.ConclusionsPathway PRS may contribute to understanding separable disease processes, but do not add significant power for predictive purposes. These findings demonstrate that AD-phenotypes may be preferentially associated with risk in specific pathways, and defining genetic risk along multiple dimensions may clarify etiological heterogeneity in AD. This approach to delineate pathway-specific PRS can be used to study other complex diseases.

Gene variants associated with longevity are also associated with protection against cognitive decline, dementia and Alzheimer's disease, suggesting that common physiologic pathways act at the interface of longevity and cognitive function. To test the hypothesis that variants in genes implicated in cognitive function may promote exceptional longevity, we performed a comprehensive 3-stage study to identify functional longevity-associated variants in ~700 candidate genes in up to 450 centenarians and 500 controls by target capture sequencing analysis. We found an enrichment of longevity-associated genes in the nPKC and NF-κB signaling pathways by gene-based association analyses. Functional analysis of the top three gene variants (NFKBIA, CLU, PRKCH) suggests that non-coding variants modulate the expression of cognate genes, thereby reducing signaling through the nPKC and NF-κB. This matches genetic studies in multiple model organisms, suggesting that the evolutionary conservation of reduced PKC and NF-κB signaling pathways in exceptional longevity may include humans.

Parkinson's disease is characterized by accumulation of α-synuclein (αSyn). Release of oligomeric/fibrillar αSyn from damaged neurons may potentiate neuronal death in part via microglial activation. Heretofore, it remained unknown if oligomeric/fibrillar αSyn could activate the nucleotide-binding oligomerization domain (NOD)-like receptor (NLR) family pyrin domain-containing 3 (NLRP3) inflammasome in human microglia and whether anti-αSyn antibodies could prevent this effect. Here, we show that αSyn activates the NLRP3 inflammasome in human induced pluripotent stem cell (hiPSC)-derived microglia (hiMG) via dual stimulation involving Toll-like receptor 2 (TLR2) engagement and mitochondrial damage. In vitro, hiMG can be activated by mutant (A53T) αSyn secreted from hiPSC-derived A9-dopaminergic neurons. Surprisingly, αSyn-antibody complexes enhanced rather than suppressed inflammasome-mediated interleukin-1β (IL-1β) secretion, indicating these complexes are neuroinflammatory in a human context. A further increase in inflammation was observed with addition of oligomerized amyloid-β peptide (Aβ) and its cognate antibody. In vivo, engraftment of hiMG with αSyn in humanized mouse brain resulted in caspase-1 activation and neurotoxicity, which was exacerbated by αSyn antibody. These findings may have important implications for antibody therapies aimed at depleting misfolded/aggregated proteins from the human brain, as they may paradoxically trigger inflammation in human microglia.

Consumer uptake of direct-to-consumer (DTC) DNA ancestry testing is accelerating, yet few empirical studies have examined test impacts on recipients despite the DTC ancestry industry being two decades old. Participants in a longitudinal cohort study of response to health-related DTC genomic testing also received personal DNA ancestry testing at no additional cost. Baseline survey data from the primary study were analyzed together with responses to an additional follow-up survey focused on the response to ancestry results. Ancestry results were generated for 3466 individuals. Of those, 1317 accessed their results, and 322 individuals completed an ancestry response survey, in other words, approximately one in ten who received ancestry testing responded to the survey. Self-reported race/ethnicity was predictive of those most likely to view their results. While 46% of survey responders (N = 147) reported their ancestry results as surprising or unexpected, less than 1% (N = 3) were distressed by them. Importantly, however, 21% (N = 67) reported that their results reshaped their personal identity. Most (81%; N = 260) planned to share results with family, and 12% (N = 39) intended to share results with a healthcare provider. Many (61%; N = 196) reported test benefits (e.g., health insights), while 12% (N = 38) reported negative aspects (e.g., lack of utility). Over half (N = 162) reported being more likely to have other genetic tests in the future. DNA ancestry testing affected individuals with respect to personal identity, intentions to share genetic information with family and healthcare providers, and the likelihood to engage with other genetic tests in the future. These findings have implications for medical care and research, specifically, provider readiness to engage with genetic ancestry information.

signatureSearch is an R/Bioconductor package that integrates a suite of existing and novel algorithms into an analysis environment for gene expression signature (GES) searching combined with functional enrichment analysis (FEA) and visualization methods to facilitate the interpretation of the search results. In a typical GES search (GESS), a query GES is searched against a database of GESs obtained from large numbers of measurements, such as different genetic backgrounds, disease states and drug perturbations. Database matches sharing correlated signatures with the query indicate related cellular responses frequently governed by connected mechanisms, such as drugs mimicking the expression responses of a disease. To identify which processes are predominantly modulated in the GESS results, we developed specialized FEA methods combined with drug-target network visualization tools. The provided analysis tools are useful for studying the effects of genetic, chemical and environmental perturbations on biological systems, as well as searching single cell GES databases to identify novel network connections or cell types. The signatureSearch software is unique in that it provides access to an integrated environment for GESS/FEA routines that includes several novel search and enrichment methods, efficient data structures, and access to pre-built GES databases, and allowing users to work with custom databases.

von Economo neurons (VENs) are bipolar, spindle-shaped neurons restricted to layer 5 of human frontoinsula and anterior cingulate cortex that appear to be selectively vulnerable to neuropsychiatric and neurodegenerative diseases, although little is known about other VEN cellular phenotypes. Single nucleus RNA-sequencing of frontoinsula layer 5 identifies a transcriptomically-defined cell cluster that contained VENs, but also fork cells and a subset of pyramidal neurons. Cross-species alignment of this cell cluster with a well-annotated mouse classification shows strong homology to extratelencephalic (ET) excitatory neurons that project to subcerebral targets. This cluster also shows strong homology to a putative ET cluster in human temporal cortex, but with a strikingly specific regional signature. Together these results suggest that VENs are a regionally distinctive type of ET neuron. Additionally, we describe the first patch clamp recordings of VENs from neurosurgically-resected tissue that show distinctive intrinsic membrane properties relative to neighboring pyramidal neurons.

BackgroundIn standard high throughput sequencing analysis, genetic variants are not assigned to a homologous chromosome of origin. This process, called haplotype phasing, can reveal information important for understanding the relationship between genetic variants and biological phenotypes. For example, in genes that carry multiple heterozygous missense variants, phasing resolves whether one or both gene copies are altered. Here, we present a novel approach to phasing variants that takes advantage of unique properties of paired tumor:normal sequencing data from cancer studies.ResultsVAF phasing uses changes in variant allele frequency (VAF) between tumor and normal samples in regions of somatic chromosomal gain or loss to phase germline variants. We apply VAF phasing to 6180 samples from the Cancer Genome Atlas (TCGA) and demonstrate that our method is highly concordant with other standard phasing methods, and can phase an average of 33% more variants than other read-backed phasing methods. Using variant annotation tools designed to score gene haplotypes, we find a suggestive association between carrying multiple missense variants in a single copy of a cancer predisposition gene and earlier age of cancer diagnosis.ConclusionsVAF phasing exploits unique properties of tumor genomes to increase the number of germline variants that can be phased over standard read-backed methods in paired tumor:normal samples. Our phase-informed association testing results call attention to the need to develop more tools for assessing the joint effect of multiple genetic variants.

Dilated cardiomyopathy (DCM) is a leading cause of morbidity and mortality worldwide; yet how genetic variation and environmental factors impact DCM heritability remains unclear. Here, we report that compound genetic interactions between DNA sequence variants contribute to the complex heritability of DCM. By using genetic data from a large family with a history of DCM, we discovered that heterozygous sequence variants in the TROPOMYOSIN 1 (TPM1) and VINCULIN (VCL) genes cose-gregate in individuals affected by DCM. In vitro studies of patient-derived and isogenic human-pluripotent-stem-cell-derived cardio-myocytes that were genome-edited via CRISPR to create an allelic series of TPM1 and VCL variants revealed that cardiomyocytes with both TPM1 and VCL variants display reduced contractility and sarcomeres that are less organized. Analyses of mice genetically engineered to harbour these human TPM1 and VCL variants show that stress on the heart may also influence the variable penetrance and expressivity of DCM-associated genetic variants in vivo. We conclude that compound genetic variants can interact combinatorially to induce DCM, particularly when influenced by other disease-provoking stressors.

BackgroundCancer susceptibility germline variants generally require somatic alteration of the remaining allele to drive oncogenesis and, in some cases, tumor mutational profiles. Whether combined germline and somatic bi-allelic alterations are universally required for germline variation to influence tumor mutational profile is unclear. Here, we performed an exome-wide analysis of the frequency and functional effect of bi-allelic alterations in The Cancer Genome Atlas (TCGA).MethodsWe integrated germline variant, somatic mutation, somatic methylation, and somatic copy number loss data from 7790 individuals from TCGA to identify germline and somatic bi-allelic alterations in all coding genes. We used linear models to test for association between mono- and bi-allelic alterations and somatic microsatellite instability (MSI) and somatic mutational signatures.ResultsWe discovered significant enrichment of bi-allelic alterations in mismatch repair (MMR) genes and identified six bi-allelic carriers with elevated MSI, consistent with Lynch syndrome. In contrast, we find little evidence of an effect of mono-allelic germline variation on MSI. Using MSI burden and bi-allelic alteration status, we reclassify two variants of unknown significance in MSH6 as potentially pathogenic for Lynch syndrome. Extending our analysis of MSI to a set of 127 DNA damage repair (DDR) genes, we identified a novel association between methylation of SHPRH and MSI burden.ConclusionsWe find that bi-allelic alterations are infrequent in TCGA but most frequently occur in BRCA1/2 and MMR genes. Our results support the idea that bi-allelic alteration is required for germline variation to influence tumor mutational profile. Overall, we demonstrate that integrating germline, somatic, and epigenetic alterations provides new understanding of somatic mutational profiles.

Recurrent acute pancreatitis (RAP) is a complex clinical syndrome with significant morbidity, unpredictable outcomes, and limited treatment options. The National Institute of Diabetes and Digestive and Kidney Disease sponsored a workshop on July 25, 2018, in Pittsburgh, Pennsylvania, to address research gaps impeding development of effective therapies for pancreatitis. The RAP working group identified challenges to clinical progress using existing definitions, risk assessment, diagnostic and severity criteria, disease trajectories, outcomes, and research methods. Recurrent acute pancreatitis includes all the risk of acute pancreatitis and often progresses to chronic pancreatitis with variable complications of chronic pain, exocrine insufficiency, diabetes, and pancreatic cancer. However, the great variability among individuals with RAP requires better precision in defining the risks, individual episodes, as well as their frequency, pathogenic pathways, and specific outcome measures for each of the systems affected by pancreatic inflammation. Because of disease complexity, few patients are similar enough for traditional studies and methods to conduct clinical trials with small sample sizes are required. The need for genetic testing, biomarker development, and better imaging methods was highlighted. Adaptive and N-of-one study designs, better endpoints, and outcome measures including patient-reported outcomes should considered early in developing future therapeutic trial design and include all stakeholders.

We describe convergent evidence from transcriptomics, morphology, and physiology for a specialized GABAergic neuron subtype in human cortex. Using unbiased single-nucleus RNA sequencing, we identify ten GABAergic interneuron subtypes with combinatorial gene signatures in human cortical layer 1 and characterize a group of human interneurons with anatomical features never described in rodents, having large 'rosehip'-like axonal boutons and compact arborization. These rosehip cells show an immunohistochemical profile (GAD1+CCK+, CNR1-SST-CALB2-PVALB-) matching a single transcriptomically defined cell type whose specific molecular marker signature is not seen in mouse cortex. Rosehip cells in layer 1 make homotypic gap junctions, predominantly target apical dendritic shafts of layer 3 pyramidal neurons, and inhibit backpropagating pyramidal action potentials in microdomains of the dendritic tuft. These cells are therefore positioned for potent local control of distal dendritic computation in cortical pyramidal neurons.

Identification of all genetic variants associated with complex traits is one of the most important goals in modern human genetics. Genome-wide association studies (GWAS) have been successfully applied to identify common variants, which thus far explain only small portion of heritability. Interests in rare variants have been increasingly growing as an answer for this missing heritability. While next-generation sequencing allows detection of rare variants, its cost is still prohibitively high to sequence a large number of human DNA samples required for rare variant association studies. In this study, we evaluated the sensitivity and specificity of sequencing for pooled DNA samples of multiple individuals (Pool-seq) as a cost-effective and robust approach for rare variant discovery. We comparatively analyzed Pool-seq vs. individual-seq of indexed target capture of up to 960 genes in ∼1000 individuals, followed by independent genotyping validation studies. We found that Pool-seq was as effective and accurate as individual-seq in detecting rare variants and accurately estimating their minor allele frequencies (MAFs). Our results suggest that Pool-seq can be used as an efficient and cost-effective method for discovery of rare variants for population-based sequencing studies in individual laboratories.

BackgroundA fundamental characteristic of multicellular organisms is the specialization of functional cell types through the process of differentiation. These specialized cell types not only characterize the normal functioning of different organs and tissues, they can also be used as cellular biomarkers of a variety of different disease states and therapeutic/vaccine responses. In order to serve as a reference for cell type representation, the Cell Ontology has been developed to provide a standard nomenclature of defined cell types for comparative analysis and biomarker discovery. Historically, these cell types have been defined based on unique cellular shapes and structures, anatomic locations, and marker protein expression. However, we are now experiencing a revolution in cellular characterization resulting from the application of new high-throughput, high-content cytometry and sequencing technologies. The resulting explosion in the number of distinct cell types being identified is challenging the current paradigm for cell type definition in the Cell Ontology.ResultsIn this paper, we provide examples of state-of-the-art cellular biomarker characterization using high-content cytometry and single cell RNA sequencing, and present strategies for standardized cell type representations based on the data outputs from these cutting-edge technologies, including "context annotations" in the form of standardized experiment metadata about the specimen source analyzed and marker genes that serve as the most useful features in machine learning-based cell type classification models. We also propose a statistical strategy for comparing new experiment data to these standardized cell type representations.ConclusionThe advent of high-throughput/high-content single cell technologies is leading to an explosion in the number of distinct cell types being identified. It will be critical for the bioinformatics community to develop and adopt data standard conventions that will be compatible with these new technologies and support the data representation needs of the research community. The proposals enumerated here will serve as a useful starting point to address these challenges.

BackgroundCancer research to date has largely focused on somatically acquired genetic aberrations. In contrast, the degree to which germline, or inherited, variation contributes to tumorigenesis remains unclear, possibly due to a lack of accessible germline variant data. Here we called germline variants on 9618 cases from The Cancer Genome Atlas (TCGA) database representing 31 cancer types.ResultsWe identified batch effects affecting loss of function (LOF) variant calls that can be traced back to differences in the way the sequence data were generated both within and across cancer types. Overall, LOF indel calls were more sensitive to technical artifacts than LOF Single Nucleotide Variant (SNV) calls. In particular, whole genome amplification of DNA prior to sequencing led to an artificially increased burden of LOF indel calls, which confounded association analyses relating germline variants to tumor type despite stringent indel filtering strategies. The samples affected by these technical artifacts include all acute myeloid leukemia and practically all ovarian cancer samples.ConclusionsWe demonstrate how technical artifacts induced by whole genome amplification of DNA can lead to false positive germline-tumor type associations and suggest TCGA whole genome amplified samples be used with caution. This study draws attention to the need to be sensitive to problems associated with a lack of uniformity in data generation in TCGA data.

Chromogranins are pro-hormone secretory proteins released from neuroendocrine cells, with effects on control of blood pressure. We conducted a genome-wide association study for plasma catestatin, the catecholamine release inhibitory peptide derived from chromogranin A (CHGA), and other CHGA- or chromogranin B (CHGB)-related peptides, in 545 US and 1252 Australian subjects. This identified loci on chromosomes 4q35 and 5q34 affecting catestatin concentration (P = 3.40 × 10-30 for rs4253311 and 1.85 × 10-19 for rs2731672, respectively). Genes in these regions include the proteolytic enzymes kallikrein (KLKB1) and Factor XII (F12). In chromaffin cells, CHGA and KLKB1 proteins co-localized in catecholamine storage granules. In vitro, kallikrein cleaved recombinant human CHGA to catestatin, verified by mass spectrometry. The peptide identified from this digestion (CHGA360-373) selectively inhibited nicotinic cholinergic stimulated catecholamine release from chromaffin cells. A proteolytic cascade involving kallikrein and Factor XII cleaves chromogranins to active compounds both in vivo and in vitro.

Background: Alzheimer's disease (AD) exhibits considerable phenotypic heterogeneity, suggesting the potential existence of subtypes. AD is under substantial genetic influence, thus identifying systematic variation in genetic risk may provide insights into disease origins. Objective: We investigated genetic heterogeneity in AD risk through a multi-step analysis. Methods: We performed principal component analysis (PCA) on AD-associated variants in the UK Biobank (AD cases = 2,739, controls = 5,478) to assess structured genetic heterogeneity. Subsequently, a biclustering algorithm searched for distinct disease-specific genetic signatures among subsets of cases. Replication tests were conducted using the Alzheimer's Disease Neuroimaging Initiative (ADNI) dataset (AD cases = 500, controls = 470). We categorized a separate set of ADNI individuals with mild cognitive impairment (MCI; n = 399) into genetic subtypes and examined cognitive, amyloid, and tau trajectories. Results: PCA revealed three distinct clusters ("constellations") driven primarily by different correlation patterns in a region of strong LD surrounding the MAPT locus. Constellations contained a mixture of cases and controls, reflecting disease-relevant but not disease-specific structure. We found two disease-specific biclusters among AD cases. Pathway analysis linked bicluster-associated variants to neuron morphogenesis and outgrowth. Disease-relevant and disease-specific structure replicated in ADNI, and bicluster 2 exhibited increased cerebrospinal fluid p-tau and cognitive decline over time. Conclusions: This study unveils a hierarchical structure of AD genetic risk. Disease-relevant constellations may represent haplotype structure that does not increase risk directly but may alter the relative importance of other genetic risk factors. Biclusters may represent distinct AD genetic subtypes. This structure is replicable and relates to differential pathological accumulation and cognitive decline over time. [ABSTRACT FROM AUTHOR]

The electrocardiographic QRS duration, a measure of ventricular depolarization and conduction, is associated with cardiovascular mortality. While single nucleotide polymorphisms (SNPs) associated with QRS duration have been identified at 22 loci in populations of European descent, the genetic architecture of QRS duration in non-European populations is largely unknown. We therefore performed a genome-wide association study (GWAS) meta-analysis of QRS duration in 13,031 African Americans from ten cohorts and a transethnic GWAS meta-analysis with additional results from populations of European descent. In the African American GWAS, a single genome-wide significant SNP association was identified (rs3922844, P = 4 × 10-14) in intron 16 of SCN5A, a voltage-gated cardiac sodium channel gene. The QRS-prolonging rs3922844 C allele was also associated with decreased SCN5A RNA expression in human atrial tissue (P = 1.1 × 10-4). High density genotyping revealed that the SCN5A association region in African Americans was confined to intron 16. Transethnic GWAS meta-analysis identified novel SNP associations on chromosome 18 in MYL12A (rs1662342, P = 4.9 × 10-8) and chromosome 1 near CD1E and SPTA1 (rs7547997, P = 7.9 × 10-9). The 22 QRS loci previously identified in populations of European descent were enriched for significant SNP associations with QRS duration in African Americans (P = 9.9 × 10-7), and index SNP associations in or near SCN5A, SCN10A, CDKN1A, NFIA, HAND1, TBX5 and SETBP1 replicated in African Americans. In summary, rs3922844 was associated with QRS duration and SCN5A expression, two novel QRS loci were identified using transethnic meta-analysis, and a significant proportion of QRS-SNP associations discovered in populations of European descent were transferable to African Americans when adequate power was achieved.

The many subcomponents of the human cortex are known to follow an anatomical pattern and functional relationship that appears to be highly conserved between individuals. This suggests that this pattern and the relationship among cortical regions are important for cortical function and likely shaped by genetic factors, although the degree to which genetic factors contribute to this pattern is unknown. We assessed the genetic relationships among 12 cortical surface areas using brain images and genotype information on 2,364 unrelated individuals, brain images on 466 twin pairs, and transcriptome data on 6 postmortem brains in order to determine whether a consistent and biologically meaningful pattern could be identified from these very different data sets. We find that the patterns revealed by each data set are highly consistent (p

The epigenome orchestrates genome accessibility, functionality, and three-dimensional structure. Because epigenetic variation can impact transcription and thus phenotypes, it may contribute to adaptation. Here, we report 1,107 high-quality single-base resolution methylomes and 1,203 transcriptomes from the 1001 Genomes collection of Arabidopsis thaliana. Although the genetic basis of methylation variation is highly complex, geographic origin is a major predictor of genome-wide DNA methylation levels and of altered gene expression caused by epialleles. Comparison to cistrome and epicistrome datasets identifies associations between transcription factor binding sites, methylation, nucleotide variation, and co-expression modules. Physical maps for nine of the most diverse genomes reveal how transposons and other structural variants shape the epigenome, with dramatic effects on immunity genes. The 1001 Epigenomes Project provides a comprehensive resource for understanding how variation in DNA methylation contributes to molecular and non-molecular phenotypes in natural populations of the most studied model plant.

Anxiety disorders peak in incidence during adolescence, a developmental window that is marked by dynamic changes in gene expression, endocannabinoid signaling, and frontolimbic circuitry. We tested whether genetic alterations in endocannabinoid signaling related to a common polymorphism in fatty acid amide hydrolase (FAAH), which alters endocannabinoid anandamide (AEA) levels, would impact the development of frontolimbic circuitry implicated in anxiety disorders. In a pediatric imaging sample of over 1,000 3- to 21-y-olds, we show effects of the FAAH genotype specific to frontolimbic connectivity that emerge by ∼12 y of age and are paralleled by changes in anxiety-related behavior. Using a knock-in mouse model of the FAAH polymorphism that controls for genetic and environmental backgrounds, we confirm phenotypic differences in frontoamygdala circuitry and anxiety-related behavior by postnatal day 45 (P45), when AEA levels begin to decrease, and also, at P75 but not before. These results, which converge across species and level of analysis, highlight the importance of underlying developmental neurobiology in the emergence of genetic effects on brain circuitry and function. Moreover, the results have important implications for the identification of risk for disease and precise targeting of treatments to the biological state of the developing brain as a function of developmental changes in gene expression and neural circuit maturation.

Dyslexia and language impairment (LI) are complex traits with substantial genetic components. We recently completed an association scan of the DYX2 locus, where we observed associations of markers in DCDC2, KIAA0319, ACOT13, and FAM65B with reading-, language-, and IQ-related traits. Additionally, the effects of reading-associated DYX3 markers were recently characterized using structural neuroimaging techniques. Here, we assessed the neuroimaging implications of associated DYX2 and DYX3 markers, using cortical volume, cortical thickness, and fractional anisotropy. To accomplish this, we examined eight DYX2 and three DYX3 markers in 332 subjects in the Pediatrics Imaging Neurocognition Genetics study. Imaging-genetic associations were examined by multiple linear regression, testing for influence of genotype on neuroimaging. Markers in DYX2 genes KIAA0319 and FAM65B were associated with cortical thickness in the left orbitofrontal region and global fractional anisotropy, respectively. KIAA0319 and ACOT13 were suggestively associated with overall fractional anisotropy and left pars opercularis cortical thickness, respectively. DYX3 markers showed suggestive associations with cortical thickness and volume measures in temporal regions. Notably, we did not replicate association of DYX3 markers with hippocampal measures. In summary, we performed a neuroimaging follow-up of reading-, language-, and IQ-associated DYX2 and DYX3 markers. DYX2 associations with cortical thickness may reflect variations in their role in neuronal migration. Furthermore, our findings complement gene expression and imaging studies implicating DYX3 markers in temporal regions. These studies offer insight into where and how DYX2 and DYX3 risk variants may influence neuroimaging traits. Future studies should further connect the pathways to risk variants associated with neuroimaging/neurocognitive outcomes.

Neurophysiological measurements of the response to pre-pulse and startle stimuli have been suggested to represent an important endophenotype for both substance dependence and other select psychiatric disorders. We have previously shown, in young adult Mexican Americans (MA), that presentation of a short delay acoustic pre-pulse, prior to the startle stimuli can elicit a late negative component at about 400 msec (N4S), in the event-related potential (ERP), recorded from frontal cortical areas. In the present study, we investigated whether genetic factors associated with this endophenotype could be identified. The study included 420 (age 18-30 years) MA men (n = 170), and women (n = 250). DNA was genotyped using an Affymetrix Axiom Exome1A chip. An association analysis revealed that the CCKAR and CCKBR (cholecystokinin A and B receptor) genes each had a nearby variant that showed suggestive significance with the amplitude of the N4S component to pre-pulse stimuli. The neurotransmitter cholecystokinin (CCK), along with its receptors, CCKAR and CCKBR, have been previously associated with psychiatric disorders, suggesting that variants near these genes may play a role in the pre-pulse/startle response in this cohort.

PurposeThe Scripps Idiopathic Diseases of Man (IDIOM) study aims to discover novel gene-disease relationships and provide molecular genetic diagnosis and treatment guidance for individuals with novel diseases using genome sequencing integrated with clinical assessment and multidisciplinary case review. Here we describe the operational protocol and initial results of the IDIOM study.MethodsA total of 121 cases underwent first-tier review by the principal investigators to determine whether the primary inclusion criteria were satisfied, 59 (48.8%) underwent second-tier review by our clinician-scientist review panel, and 17 patients (14.0%) and their family members were enrolled.Results60% of cases resulted in a plausible molecular diagnosis, and 18% of cases resulted in a confirmed molecular diagnosis. Two of three confirmed cases led to the identification of novel gene-disease relationships. In the third confirmed case a previously described but unrecognized disease was revealed. In all three confirmed cases a new clinical management strategy was initiated based on the genetic findings.ConclusionGenome sequencing provides tangible clinical benefit for individuals with idiopathic genetic disease, not only in the context of molecular genetic diagnosis of known rare conditions but also in cases where prior clinical information regarding a new genetic disorder is lacking.

BackgroundMitochondrial DNA (mtDNA) heteroplasmy is a mixture of normal and mutated mtDNA molecules in a cell. High levels of heteroplasmy at specific mtDNA sites lead to inherited mitochondrial diseases with neurological, sensory, and movement impairments. Here we test the hypothesis that heteroplasmy levels in elderly adults are associated with impaired function resembling mild forms of mitochondrial disease.MethodsWe examined platelet mtDNA heteroplasmy at 20 disease-causing sites for associations with neurosensory and mobility function among 137 participants from the community-based Health, Aging, and Body Composition Study.ResultsElevated mtDNA heteroplasmy at four mtDNA sites in complex I and tRNA genes was nominally associated with reduced cognition, vision, hearing, and mobility: m.10158T>C with Modified Mini-Mental State Examination score (p = .009); m.11778G>A with contrast sensitivity (p = .02); m.7445A>G with high-frequency hearing (p = .047); and m.5703G>A with 400 m walking speed (p = .007).ConclusionsThese results indicate that increased mtDNA heteroplasmy at disease-causing sites is associated with neurosensory and mobility function in older persons. We propose the novel use of mtDNA heteroplasmy as a simple, noninvasive predictor of age-related neurologic, sensory, and movement impairments.

Knowing how the human brain is shaped by migration and admixture is a critical step in studying human evolution [1, 2], as well as in preventing the bias of hidden population structure in brain research [3, 4]. Yet, the neuroanatomical differences engendered by population history are still poorly understood. Most of the inference relies on craniometric measurements, because morphology of the brain is presumed to be the neurocranium's main shaping force before bones are fused and ossified [5]. Although studies have shown that the shape variations of cranial bones are consistent with population history [6-8], it is unknown how much human ancestry information is retained by the human cortical surface. In our group's previous study, we found that area measures of cortical surface and total brain volumes of individuals of European descent in the United States correlate significantly with their ancestral geographic locations in Europe [9]. Here, we demonstrate that the three-dimensional geometry of cortical surface is highly predictive of individuals' genetic ancestry in West Africa, Europe, East Asia, and America, even though their genetic background has been shaped by multiple waves of migratory and admixture events. The geometry of the cortical surface contains richer information about ancestry than the areal variability of the cortical surface, independent of total brain volumes. Besides explaining more ancestry variance than other brain imaging measurements, the 3D geometry of the cortical surface further characterizes distinct regional patterns in the folding and gyrification of the human brain associated with each ancestral lineage.

Socioeconomic disparities are associated with differences in cognitive development. The extent to which this translates to disparities in brain structure is unclear. We investigated relationships between socioeconomic factors and brain morphometry, independently of genetic ancestry, among a cohort of 1,099 typically developing individuals between 3 and 20 years of age. Income was logarithmically associated with brain surface area. Among children from lower income families, small differences in income were associated with relatively large differences in surface area, whereas, among children from higher income families, similar income increments were associated with smaller differences in surface area. These relationships were most prominent in regions supporting language, reading, executive functions and spatial skills; surface area mediated socioeconomic differences in certain neurocognitive abilities. These data imply that income relates most strongly to brain structure among the most disadvantaged children.

We sequenced the genomes of 200 individuals from 41 families multiply affected with bipolar disorder (BD) to identify contributions of rare variants to genetic risk. We initially focused on 3,087 candidate genes with known synaptic functions or prior evidence from genome-wide association studies. BD pedigrees had an increased burden of rare variants in genes encoding neuronal ion channels, including subunits of GABAA receptors and voltage-gated calcium channels. Four uncommon coding and regulatory variants also showed significant association, including a missense variant in GABRA6. Targeted sequencing of 26 of these candidate genes in an additional 3,014 cases and 1,717 controls confirmed rare variant associations in ANK3, CACNA1B, CACNA1C, CACNA1D, CACNG2, CAMK2A, and NGF. Variants in promoters and 5' and 3' UTRs contributed more strongly than coding variants to risk for BD, both in pedigrees and in the case-control cohort. The genes and pathways identified in this study regulate diverse aspects of neuronal excitability. We conclude that rare variants in neuronal excitability genes contribute to risk for BD.

BackgroundMigraine is a common comorbidity among individuals with bipolar disorder, but the underlying mechanisms for this co-occurrence are poorly understood. The aim of this study was to investigate the genetic background of bipolar patients with and without migraine.MethodsWe performed a genome-wide association analysis contrasting 460 bipolar migraneurs with 914 bipolar patients without migraine from the Bipolar Genome Study (BiGS).ResultsWe identified one genome-wide significant association between migraine in bipolar disorder patients and rs1160720, an intronic single nucleotide polymorphism (SNP) in the NBEA gene (P=2.97 × 10(-8), OR: 1.82, 95% CI: 1.47-2.25), although this was not replicated in a smaller sample of 289 migraine cases.LimitationsOur study is based on self-reported migraine.ConclusionsNBEA encodes neurobeachin, a scaffolding protein primarily expressed in the brain and involved in trafficking of vesicles containing neurotransmitter receptors. This locus has not previously been implicated in migraine per se. We found no evidence of association in data from the GWAS migraine meta-analysis consortium (n=118,710 participants) suggesting that the association might be specific to migraine co-morbid with bipolar disorder.

BackgroundResearch on the etiology of post-traumatic stress disorder (PTSD) has rapidly matured, moving from candidate gene studies to interrogation of the entire human genome in genome-wide association studies (GWAS). Here we present the results of a GWAS performed on samples from combat-exposed U.S. Marines and Sailors from the Marine Resiliency Study (MRS) scheduled for deployment to Iraq and/or Afghanistan. The MRS is a large, prospective study with longitudinal follow-up designed to identify risk and resiliency factors for combat-induced stress-related symptoms. Previously implicated PTSD risk loci from the literature and polygenic risk scores across psychiatric disorders were also evaluated in the MRS cohort.MethodsParticipants (N=3494) were assessed using the Clinician-Administered PTSD Scale and diagnosed using the DSM-IV diagnostic criterion. Subjects with partial and/or full PTSD diagnosis were called cases, all other subjects were designated controls, and study-wide maximum CAPS scores were used for longitudinal assessments. Genomic DNA was genotyped on the Illumina HumanOmniExpressExome array. Individual genetic ancestry was determined by supervised cluster analysis for subjects of European, African, Hispanic/Native American, and other descent. To test for association of SNPs with PTSD, logistic regressions were performed within each ancestry group and results were combined in meta-analyses. Measures of childhood and adult trauma were included to test for gene-by-environment (GxE) interactions. Polygenic risk scores from the Psychiatric Genomic Consortium were used for major depressive disorder (MDD), bipolar disorder (BPD), and schizophrenia (SCZ).ResultsThe array produced >800K directly genotyped and >21M imputed markers in 3494 unrelated, trauma-exposed males, of which 940 were diagnosed with partial or full PTSD. The GWAS meta-analysis identified the phosphoribosyl transferase domain containing 1 gene (PRTFDC1) as a genome-wide significant PTSD locus (rs6482463; OR=1.47, SE=0.06, p=2.04×10(-9)), with a similar effect across ancestry groups. Association of PRTFDC1 with PTSD in an independent military cohort showed some evidence for replication. Loci with suggestive evidence of association (n=25 genes, p

BackgroundMexican Americans, particularly those born in the United States, are at greater risk for alcohol associated morbidity and mortality. The present study sought to investigate whether specific genetic variants may be associated with alcohol use disorder phenotypes in a select population of Mexican American young adults.MethodsThe study evaluated a cohort of 427 (age 18 - 30 years) Mexican American men (n = 171) and women (n = 256). Information on alcohol dependence was obtained through interview using the Semi-Structured Assessment for the Genetics of Alcoholism (SSAGA). For all subjects, DNA was extracted from blood samples, followed by genotyping using an Affymetrix Axiom Exome1A chip.ResultsA protective variant (rs991316) located downstream from the ADH7 (alcohol dehydrogenase 7) gene showed suggestive significance in association with alcohol dependence symptom counts derived from DSM-III-R and DSM-IV criteria, as well as to clustered alcohol dependence symptoms. Additional linkage analysis suggested that nearby variants in linkage disequilibrium with rs991316 were not responsible for the observed association with the alcohol dependence phenotypes in this study.ConclusionsADH7 has been shown to have a protective role against alcohol dependence in previous studies involving other ethnicities, but has not been reported for Mexican Americans. These results suggest that variants near ADH7 may play a role in protection from alcohol dependence in this Mexican American cohort.