Public platform with 39,472 exome control samples enables association studies without genotype sharing.

Acquiring a sufficiently powered cohort of control samples matched to a case sample can be time-consuming or, in some cases, impossible. Accordingly, an ability to leverage genetic data from control samples that were already collected elsewhere could dramatically improve power in genetic association studies. Sharing of control samples can pose significant challenges, since most human genetic data are subject to strict sharing regulations. Here, using the properties of singular value decomposition and subsampling algorithm, we developed a method allowing selection of the best-matching controls in an external pool of samples compliant with personal data protection and eliminating the need for genotype sharing. We provide access to a library of 39,472 exome sequencing controls at http://dnascore.net enabling association studies for case cohorts lacking control subjects. Using this approach, control sets can be selected from this online library with a prespecified matching accuracy, ensuring well-calibrated association analysis for both rare and common variants.

CR1 variants contribute to FSGS susceptibility across multiple populations.

Focal segmental glomerulosclerosis (FSGS) is a common cause of nephrotic syndrome with an annual incidence in the United States in African-Americans compared to European-Americans of 24 cases and 5 cases per million, respectively. Among glomerular diseases in Europe and Latin-America, FSGS was the second most frequent diagnosis, and in Asia the fifth. We expand previous efforts in understanding genetics of FSGS by performing a case-control study involving ethnically-diverse groups FSGS cases (726) and a pool of controls (13,994), using panel sequencing of approximately 2,500 podocyte-expressed genes. Through rare variant association tests, we replicated known risk genes - KANK1, COL4A4, and APOL1. A novel significant association was observed for the gene encoding complement receptor 1 (CR1). High-risk rare variants in CR1 in the European-American cohort were commonly observed in Latin- and African-Americans. Therefore, a combined rare and common variant analysis was used to replicate the CR1 association in non-European populations. The CR1 risk variant, rs17047661, gives rise to the Sl1/Sl2 (R1601G) allele that was previously associated with protection against cerebral malaria. Pleiotropic effects of rs17047661 may explain the difference in allele frequencies across continental ancestries and suggest a possible role for genetically-driven alterations of adaptive immunity in the pathogenesis of FSGS.

Discordant calls across genotype discovery approaches elucidate variants with systematic errors.

Large-scale high-throughput sequencing data sets have been transformative for informing clinical variant interpretation and for use as reference panels for statistical and population genetic efforts. Although such resources are often treated as ground truth, we find that in widely used reference data sets such as the Genome Aggregation Database (gnomAD), some variants pass gold-standard filters, yet are systematically different in their genotype calls across genotype discovery approaches. The inclusion of such discordant sites in study designs involving multiple genotype discovery strategies could bias results and lead to false-positive hits in association studies owing to technological artifacts rather than a true relationship to the phenotype. Here, we describe this phenomenon of discordant genotype calls across genotype discovery approaches, characterize the error mode of wrong calls, provide a list of discordant sites identified in gnomAD that should be treated with caution in analyses, and present a metric and machine learning classifier trained on gnomAD data to identify likely discordant variants in other data sets. We find that different genotype discovery approaches have different sets of variants at which this problem occurs, but there are characteristic variant features that can be used to predict discordant behavior. Discordant sites are largely shared across ancestry groups, although different populations are powered for the discovery of different variants. We find that the most common error mode is that of a variant being heterozygous for one approach and homozygous for the other, with heterozygous in the genomes and homozygous reference in the exomes making up the majority of miscalls.

Gene, cell type, and drug prioritization analysis suggest genetic basis for the utility of diuretics in treating Alzheimer disease.

We introduce a user-friendly tool for risk gene, cell type, and drug prioritization for complex traits: GCDPipe. It uses gene-level GWAS-derived data and gene expression data to train a model for the identification of disease risk genes and relevant cell types. Gene prioritization information is then coupled with known drug target data to search for applicable drug agents based on their estimated functional effects on the identified risk genes. We illustrate the utility of our approach in different settings: identification of the cell types, implicated in disease pathogenesis, was tested in inflammatory bowel disease (IBD) and Alzheimer disease (AD); gene target and drug prioritization was tested in IBD and schizophrenia. The analysis of phenotypes with known disease-affected cell types and/or existing drug candidates shows that GCDPipe is an effective tool to unify genetic risk factors with cellular context and known drug targets. Next, analysis of the AD data with GCDPipe suggested that gene targets of diuretics, as an Anatomical Therapeutic Chemical drug subgroup, are significantly enriched among the genes prioritized by GCDPipe, indicating their possible effect on the course of the disease.

The 22q11.2 region regulates presynaptic gene-products linked to schizophrenia.

It is unclear how the 22q11.2 deletion predisposes to psychiatric disease. To study this, we generated induced pluripotent stem cells from deletion carriers and controls and utilized CRISPR/Cas9 to introduce the heterozygous deletion into a control cell line. Here, we show that upon differentiation into neural progenitor cells, the deletion acted in trans to alter the abundance of transcripts associated with risk for neurodevelopmental disorders including autism. In excitatory neurons, altered transcripts encoded presynaptic factors and were associated with genetic risk for schizophrenia, including common and rare variants. To understand how the deletion contributed to these changes, we defined the minimal protein-protein interaction network that best explains gene expression alterations. We found that many genes in 22q11.2 interact in presynaptic, proteasome, and JUN/FOS transcriptional pathways. Our findings suggest that the 22q11.2 deletion impacts genes that may converge with psychiatric risk loci to influence disease manifestation in each deletion carrier.

Genotype imputation and polygenic score estimation in northwestern Russian population.

Numerous studies demonstrated the lack of transferability of polygenic score (PGS) models across populations and the problem arising from unequal presentation of ancestries across genetic studies. However, even within European ancestry there are ethnic groups that are rarely presented in genetic studies. For instance, Russians, being one of the largest, diverse, and yet understudied group in Europe. In this study, we evaluated the reliability of genotype imputation for the Russian cohort by testing several commonly used imputation reference panels (e.g. HRC, 1000G, HGDP). HRC, in comparison with two other panels, showed the most accurate results based on both imputation accuracy and allele frequency concordance between masked and imputed genotypes. We built polygenic score models based on GWAS results from the UK biobank, measured the explained phenotypic variance in the Russian cohort attributed to polygenic scores for 11 phenotypes, collected in the clinic for each participant, and finally explored the role of allele frequency discordance between the UK biobank and the study cohort in the resulting PGS performance.

Shiny GATOM: omics-based identification of regulated metabolic modules in atom transition networks.

Multiple high-throughput omics techniques provide different angles on systematically quantifying and studying metabolic regulation of cellular processes. However, an unbiased analysis of such data and, in particular, integration of multiple types of data remains a challenge. Previously, for this purpose we developed GAM web-service for integrative metabolic network analysis. Here we describe an updated pipeline GATOM and the corresponding web-service Shiny GATOM, which takes as input transcriptional and/or metabolomic data and finds a metabolic subnetwork most regulated between the two conditions of interest. GATOM features a new metabolic network topology based on atom transition, which significantly improves interpretability of the analysis results. To address computational challenges arising with the new network topology, we introduce a new variant of the maximum weight connected subgraph problem and provide a corresponding exact solver. To make the used networks up-to-date we upgraded the KEGG-based network construction pipeline and developed one based on the Rhea database, which allows analysis of lipidomics data. Finally, we simplified local installation, providing R package mwcsr for solving relevant graph optimization problems and R package gatom, which implements the GATOM pipeline. The web-service is available at https://ctlab.itmo.ru/shiny/gatom and https://artyomovlab.wustl.edu/shiny/gatom.

Immuno Tomography (IT) and Imaging Mass Cytometry (IMC) for constructing spatially resolved, multiplexed 3D IMC data sets.

PURPOSE: We have previously used Immuno Tomography (IT) to identify label-retaining stem cell populations in the cornea and meibomian gland. While this method provides the unique ability to quantify stem cell populations comprised of 1-4 cells, the number of antigens that can be sequentially used to characterize these unique cells is limited by antigen stability after antibody stripping and re-probing. To address this deficiency, we have evaluated the capability of Imaging Mass Cytometry™ (IMC™) to generate multiplexed images using metal-conjugated antibodies to label IT plastic sections and generate 3-dimensional IMC data sets (3D IMC). METHODS: K5-H2B-GFP mice, 56 days after doxycycline chase, were sacrificed and eyelid tissue processed for IT. A total of 400 serial, plastic sections, 2 µm thick, were then probed using metal-tagged antibodies specific for sox 9, collagen type I, E-cadherin, Ki67, GFP, αSMA, vimentin, and DNA intercalator. Multiplexed images were then generated using an Imaging Mass Cytometry system (Fluidigm®), and 3D reconstructions were assembled. RESULTS: All 8 metal-labeled tags were detected and their images were successfully assembled into 3D IMC data sets. GFP-labeled nuclei were identified within the meibomian glands in comparable numbers to those previously reported for slow-cycling meibomian gland stem cells. CONCLUSIONS: These findings demonstrate that IMC can be used on plastic sections to generate multiplexed, 3D data sets that can be reconstructed to show the spatial localization of meibomian gland stem cells. We propose that 3D IMC might prove valuable in more fully characterizing stem cell populations in different tissues.

TREM2 Maintains Microglial Metabolic Fitness in Alzheimer's Disease.

Elevated risk of developing Alzheimer's disease (AD) is associated with hypomorphic variants of TREM2, a surface receptor required for microglial responses to neurodegeneration, including proliferation, survival, clustering, and phagocytosis. How TREM2 promotes such diverse responses is unknown. Here, we find that microglia in AD patients carrying TREM2 risk variants and TREM2-deficient mice with AD-like pathology have abundant autophagic vesicles, as do TREM2-deficient macrophages under growth-factor limitation or endoplasmic reticulum (ER) stress. Combined metabolomics and RNA sequencing (RNA-seq) linked this anomalous autophagy to defective mammalian target of rapamycin (mTOR) signaling, which affects ATP levels and biosynthetic pathways. Metabolic derailment and autophagy were offset in vitro through Dectin-1, a receptor that elicits TREM2-like intracellular signals, and cyclocreatine, a creatine analog that can supply ATP. Dietary cyclocreatine tempered autophagy, restored microglial clustering around plaques, and decreased plaque-adjacent neuronal dystrophy in TREM2-deficient mice with amyloid-ß pathology. Thus, TREM2 enables microglial responses during AD by sustaining cellular energetic and biosynthetic metabolism.

Imaging Mass Cytometry.

Imaging Mass Cytometry (IMC) is an expansion of mass cytometry, but rather than analyzing single cells in suspension, it uses laser ablation to generate plumes of particles that are carried to the mass cytometer by a stream of inert gas. Images reconstructed from tissue sections scanned by IMC have a resolution comparable to light microscopy, with the high content of mass cytometry enabled through the use of isotopically labeled probes and ICP-MS detection. Importantly, IMC can be performed on paraffin-embedded tissue sections, so can be applied to the retrospective analysis of patient cohorts whose outcome is known, and eventually to personalized medicine. Since the original description in 2014, IMC has evolved rapidly into a commercial instrument of unprecedented power for the analysis of histological sections. In this Review, we discuss the underlying principles of this new technology, and outline emerging applications of IMC in the analysis of normal and pathological tissues. © 2017 International Society for Advancement of Cytometry.

GAM: a web-service for integrated transcriptional and metabolic network analysis.

Novel techniques for high-throughput steady-state metabolomic profiling yield information about changes of nearly thousands of metabolites. Such metabolomic profiles, when analyzed together with transcriptional profiles, can reveal novel insights about underlying biological processes. While a number of conceptual approaches have been developed for data integration, easily accessible tools for integrated analysis of mammalian steady-state metabolomic and transcriptional data are lacking. Here we present GAM ('genes and metabolites'): a web-service for integrated network analysis of transcriptional and steady-state metabolomic data focused on identification of the most changing metabolic subnetworks between two conditions of interest. In the web-service, we have pre-assembled metabolic networks for humans, mice, Arabidopsis and yeast and adapted exact solvers for an optimal subgraph search to work in the context of these metabolic networks. The output is the most regulated metabolic subnetwork of size controlled by false discovery rate parameters. The subnetworks are then visualized online and also can be downloaded in Cytoscape format for subsequent processing. The web-service is available at: https://artyomovlab.wustl.edu/shiny/gam/.

A novel ion trap that enables high duty cycle and wide m/z range on an orthogonal injection TOF mass spectrometer.

Although TOF analyzers with orthogonal ion injection provide the whole spectrum without scanning, their duty cycle is low compared with scanning analyzers in single ion monitoring mode. Typical duty cycle is in the range of 5% to 30% depending on the instrument geometry and ion m/z value. We present here a novel trapping/releasing setup, which offers the duty cycle near 100% over a wide range. Operation in the mass range from m/z 120 to almost 2000 is demonstrated. Ions are trapped in a short linear ion trap at the end of the collision cell in an axial pseudopotential well created by additional rf ("AC") voltage applied to all four rods of the trap with the same amplitude and phase. The pseudopotential created by AC field is mass dependent, and by ramping down the AC voltage, ions can be released from the trap sequentially from high m/z to low, while all ions are gaining the same kinetic energy. Upon entering the TOF accelerator, ions with lower m/z catch up with heavier ions, and the AC ramp parameters can be selected to make all ions meet in the center of the TOF extraction region, resulting in sensitivity gains from 3 to 14 without loss of mass accuracy or resolution.

Design considerations for high speed quantitative mass spectrometry with MALDI ionization.

A MALDI ion source on a triple quadrupole mass spectrometer constructed for the purpose of obtaining high speed quantitative measurements on drugs and other low molecular weight compounds is described. Particular attention is given to the ion generation and transport phenomena that affect analysis speed, throughput, and practical instrument robustness. In this regard parameters that affect desorption speed, beam spreading, ion flight times, sensitivity, signal-to-noise, ion fragmentation, sample carry-over, and instrument contamination are examined and experimental results are provided. MALDI and electrospray sensitivity is compared, to provide a practical frame of reference.

Novel ion mobility setup combined with collision cell and time-of-flight mass spectrometer.

An ion mobility cell of a novel type was coupled to an orthogonal injection time-of-flight (TOF) mass spectrometer. The mobility cell operates at low-pressure and contains a segmented RF ion guide providing an axial electric field that drives the ions towards the exit. A flow of gas is arranged inside the ion guide in such a way that the gas drag counteracts the force exerted by the axial field. Ions with different mobility coefficients can be scanned out of the ion guide by ramping the axial field strength. The ions can be analyzed intact or fragmented in a collision cell before introduction into an orthogonal TOF mass spectrometer. An ion source with matrix assisted laser desorption/ionization (MALDI) was attached to the instrument. The setup was evaluated for the analysis of peptide and protein mixture, with sequential fragmentation of multiple precursor ions from a protein digest and with mobility separation of fragment ions formed by in-source fragmentation of pure peptides. The mobility resolution for peptides was observed to be three times higher than the theoretical resolution predicted for a classical mobility setup with similar operating conditions (pressure, field strength, and length).

Clustering of pentacene and functionalized pentacene ions in a matrix-assisted laser desorption/ionization orthogonal TOF mass spectrometer.

A high-performance orthogonal time-of flight (TOF) mass spectrometer, in combination with the matrix assisted laser desorption/ionization (MALDI) source operating at elevated pressure ( approximately 1 torr in N(2)), was used to perform MALDI-TOF analyses of pentacene and some of its derivatives with and without an added matrix. These molecules are among the most interesting semiconductor materials for organic thin film transistor applications (OTFT). The observation of ion-molecule reactions between "cold" analyte ions and neutral analyte molecules in the gas phase has provided some insight into the mechanism of pentacene cluster formation and its functionalized derivatives. Furthermore, some of the matrices employed to assist the desorption/ionization process of these compounds were observed to influence the outcome via ion-molecule reactions of analyte ions and matrix molecules in the gas phase. The stability and reactivity of the compounds and their clusters in the MALDI plume during gas-phase expansion were evaluated; possible structures of the resulting clusters are discussed. The MALDI-TOF technique was also helpful in distinguishing between two isomeric forms of bis-[(triisopropylsilyl)-ethynyl]-pentacene.

Sequencing of rat liver cytosolic proteins by matrix-assisted laser desorption ionization-quadrupole time-of-flight mass spectrometry following electrophoretic separation and extraction.

A new technique is described that enables the direct determination of the complete or partial amino acid sequence of cytosolic proteins separated by gel electrophoresis and allows for the further observation of disease- or drug-induced posttranslational modifications. The procedure uses a two-phase extraction strategy (ethyl acetate/ammonium bicarbonate) for the efficient separation of proteins/peptides from an electrophoretic matrix and subsequent sequence analysis by matrix-assisted laser desorption ionization-quadrupole time-of-flight mass spectrometry. The method was tested using hepatocyte cytosolic proteins and compared to a complementary approach using direct solvent extraction from in-gel digests. Although the latter procedure identified the proteins, it did not enable complete amino acid sequence determination. In contrast, high sequence coverage was obtained using the peptide extraction procedure, without any apparent dependence on protein size. The technique minimized the chemically inconsistent modifications generated from in-gel digestion, thus aiding mass spectrometric interpretation and valid protein sequence identification.