Search

Justina P Tavana,1,* Matthew Rosene,2,* Nick O Jensen,2 Perry G Ridge,1 John SK Kauwe,1,3 Celeste M Karch2,4 1Department of Biology, Brigham Young University, Provo, UT 84602, USA; 2Department of Psychiatry, Washington University in St Louis, St Louis, MO, USA; 3Department of Neuroscience, Brigham Young University, Provo, UT 84602, USA; 4Hope Center for Neurological Disorders, Washington University in St Louis, St Louis, MO, USA *These authors contributed equally to this work Abstract: Alzheimer’s disease (AD) is mainly a late-onset neurodegenerative disorder. Substantial efforts have been made to solve the complex genetic architecture of AD as a means to identify therapeutic targets. Unfortunately, to date, no disease-altering therapeutics have been developed. As therapeutics are likely to be most effective in the early stages of disease (ie, before the onset of symptoms), a recent focus of AD research has been the identification of protective factors that prevent disease. One example is the discovery of a rare variant in the 3'-UTR of RAB10 that is protective for AD. Here, we review the possible genetic, molecular, and functional role of RAB10 in AD and potential therapeutic approaches to target RAB10. Keywords: Alzheimer’s disease, RAB10, retromer, APP, resilience, GTPase

Individuals of Pacific ancestry suffer some of the highest rates of health disparities yet remain vastly underrepresented in genomic research, including currently available linear and pangenome references. To begin addressing this, we developed the first Pacific ancestry pangenome reference using 23 individuals with diverse Pacific ancestry. We assembled 46 haploid genomes from these 23 individuals, resulting in highly accurate and contiguous genome assemblies with an average quality value of 55.0 and an average N50 of 40.7 Mb, marking the first de novo assembly of highly accurate Pacific ancestry genomes. We combined these assemblies to create a pangenome reference, which added 30.6 Mb of novel sequence missing from the Human Pangenome Reference Consortium (HPRC) reference. Mapping short reads to this pangenome reduced variant call errors and yielded more true-positive variants compared to the HPRC and T2T-CHM13 references. This Pacific ancestry pangenome reference serves as a resource to enhance genetic analyses for this underserved population.

Computer programming is a fundamental tool for life scientists, allowing them to carry out essential research tasks. However, despite various educational efforts, learning to write code can be a challenging endeavor for students and researchers in life-sciences disciplines. Recent advances in artificial intelligence have made it possible to translate human-language prompts to functional code, raising questions about whether these technologies can aid (or replace) life scientists' efforts to write code. Using 184 programming exercises from an introductory-bioinformatics course, we evaluated the extent to which one such tool-OpenAI's ChatGPT-could successfully complete programming tasks. ChatGPT solved 139 (75.5%) of the exercises on its first attempt. For the remaining exercises, we provided natural-language feedback to the model, prompting it to try different approaches. Within 7 or fewer attempts, ChatGPT solved 179 (97.3%) of the exercises. These findings have implications for life-sciences education and research. Instructors may need to adapt their pedagogical approaches and assessment techniques to account for these new capabilities that are available to the general public. For some programming tasks, researchers may be able to work in collaboration with machine-learning models to produce functional code., Competing Interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests: After completing the analyses for this paper, OpenAI contacted us and offered credits to use their application programming interface for additional research. We accepted those credits but did not use the credits in this study, nor did this arrangement influence our conclusions in this study., (Copyright: © 2023 Piccolo et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Much of today's molecular science revolves around next-generation sequencing. Frequently, the first step in analyzing such data is aligning sequencing reads to a reference genome. This step is often taken for granted, but any analysis downstream of the alignment will be affected by the aligner's ability to correctly map sequences. In most cases, for research into chromatin structure and nucleosome positioning, ATAC-seq, ChIP-seq, and MNase-seq experiments use short read lengths. How well aligners manage these reads is critical. Most aligner programs will output mapped reads and unmapped reads. However, from a biological point of view, reads will fall into one of three categories: correctly mapped, incorrectly mapped, and unmapped. While increased sequencing depth can often compensate for unmapped reads, incorrectly and correctly mapped reads appear algorithmically identical but can produce biologically significant alterations in the results. For this reason, we are benchmarking various alignment programs to determine their propensity to incorrectly map short reads. As short-read alignment is an important step in ATAC-seq, ChIP-seq, and MNase-seq experiments, caution should be taken in mapping reads to ensure that the most accurate conclusions can be made from the data generated. Our analysis is intended to help investigators new to the field pick the alignment program best suited for their experimental conditions. In general, the aligners we tested performed well. BWA, Bowtie2, and Chromap were all exceptionally accurate, and we recommend using them. Furthermore, we show that longer read lengths do in fact lead to more accurate mappings.

Caranx ignobilis , commonly known as giant kingfish or giant trevally, is a large, reef-associated apex predator. It is a prized sportfish, targeted throughout its tropical and subtropical range in the Indian and Pacific Oceans. It also gained significant interest in aquaculture due to its unusual freshwater tolerance. Here, we present a draft assembly of the estimated 625.92 Mbp nuclear genome of a C. ignobilis individual from Hawaiian waters, which host a genetically distinct population. Our 97.4% BUSCO-complete assembly has a contig NG50 of 7.3 Mbp and a scaffold NG50 of 46.3 Mbp. Twenty-five of the 203 scaffolds contain 90% of the genome. We also present noisy, long-read DNA, Hi-C, and RNA-seq datasets, the latter containing eight distinct tissues and can help with annotations and studies of freshwater tolerance. Our genome assembly and its supporting data are valuable tools for ecological and comparative genomics studies of kingfishes and other carangoid fishes., Competing Interests: The authors declare no competing interests., (© The Author(s) 2022.)

Ramp sequences occur when the average translational efficiency of codons near the 5' end of highly expressed genes is significantly lower than the rest of the gene sequence, which counterintuitively increases translational efficiency by decreasing downstream ribosomal collisions. Here, we show that the relative codon adaptiveness within different tissues changes the existence of a ramp sequence without altering the underlying genetic code. We present the first comprehensive analysis of tissue and cell type-specific ramp sequences and report 3108 genes with ramp sequences that change between tissues and cell types, which corresponds with increased gene expression within those tissues and cells. The Ramp Atlas (https://ramps.byu.edu/) allows researchers to query precomputed ramp sequences in 18 388 genes across 62 tissues and 66 cell types and calculate tissue-specific ramp sequences from user-uploaded FASTA files through an intuitive web interface. We used The Ramp Atlas to identify seven SARS-CoV-2 genes and seven human SARS-CoV-2 entry factor genes with tissue-specific ramp sequences that may help explain viral proliferation within those tissues. We anticipate that The Ramp Atlas will facilitate personalized and creative tissue-specific ramp sequence analyses for both human and viral genes that will increase our ability to utilize this often-overlooked regulatory region., (© The Author(s) 2022. Published by Oxford University Press on behalf of NAR Genomics and Bioinformatics.)

The roundjaw bonefish, Albula glossodonta , is the most widespread albulid in the Indo-Pacific and is vulnerable to extinction. We assembled the genome of a roundjaw bonefish from Hawai'i, USA, which will be instrumental for effective transboundary management and conservation when paired with population genomics datasets. The 1.05 gigabase pair (Gbp) contig-level assembly had a 4.75 megabase pair (Mbp) NG50 and a maximum contig length of 28.2 Mbp. Scaffolding yielded an LG50 of 20 and an NG50 of 14.49 Mbp, with the longest scaffold reaching 42.29 Mbp. The genome comprised 6.5% repetitive elements and was annotated with 28.3 K protein-coding genes. We then evaluated population genetic connectivity between six atolls in the Western Indian Ocean with 38,355 SNP loci across 66 A. glossodonta individuals. We discerned shallow population structure and observed genetic homogeneity between atolls in Seychelles and reduced gene flow between Seychelles and Mauritius. The South Equatorial Current might be the limiting mechanism of this reduced gene flow. The genome assembly will be useful for addressing taxonomic uncertainties of bonefishes globally., Competing Interests: The authors declare that they have no competing interests., (© The Author(s) 2022.)

Introduction: Recent clinical trials are considering inclusion of more than just apolipoprotein E ( APOE ) ε4 genotype as a way of reducing variability in analysis of outcomes., Methods: Case-control data were used to compare the capacity of age, sex, and 58 Alzheimer's disease (AD)-associated single nucleotide polymorphisms (SNPs) to predict AD status using several statistical models. Model performance was assessed with Brier scores and tenfold cross-validation. Genotype and sex × age estimates from the best performing model were combined with age and intercept estimates from the general population to develop a personalized genetic risk score, termed age, and sex-adjusted GenoRisk., Results: The elastic net model that included age, age x sex interaction, allelic APOE terms, and 29 additional SNPs performed the best. This model explained an additional 19% of the heritable risk compared to APOE genotype alone and achieved an area under the curve of 0.747., Discussion: GenoRisk could improve the risk assessment of individuals identified for prevention studies., Competing Interests: Suzanne Hendrix, as the owner of Pentara, received support from Affirmative Diagnostics for the manuscript and received contract work from AC Immmune, Acumen, ADCS, ADDF, Affirmativ Dignpostics, Alector, Alkahest, Allergan, Alzheon, Amylyx, Apodemus, Athira, Avanir, Banner, Biogen, Biohaven, Cadent, Capricor, Cerecin, Cognito, Cognoptix, Cortexyme, Elusis, Green Valley, Grifols, Ionis, Janssen, Kyowa‐Kirin, Lexeo, LuMind, Lundbeck, NCRI, Nilvad, NovoNordisk, Nutricia, Photopharmics, PTCBio, Regenera, Retrotope, Richard Isaacson, Samus, Soleno, Spaulding, Suven, Takeda, Telocyte, Tetra Discovery, Toyama, United Neuroscience, Vaccinex, Vacinity, vTv Therapeutics, and WGU. Suzanne Hendrix has also received royalties from Wiley in the last 36 months; received honoraria for lectures, presentations, etc. from Biogen; received support in the last 36 months from Grifols and Nutricia for presentations and conference attendance; claims one patent pending Brown BL, Hedges DW, Hendrix SB (2014) Extracting A periodic Components from a Time Series Wave Data Set—Patent Pending with a filing date of October 16, 2013, Application Number PCT/US2013/065327; and participated in Alzheon, Cortexyme, and Janssen DSMB or Advisory Board in the last 36 months. Samuel P. Dickson received support from Pentara, which in turn received support from Affirmativ Diagnostics for the manuscript. S. B. Booth received support from Affirmative Diagnostics for the manuscript; received travel support from her employer for work events; ADX provides profit sharing incentives with additional salary which were broadly applied to any work performed for the company (including manuscript generation). J. S. K. Kauwe received support for the manuscript from Brigham Young university and has received honoraria from UH Hilo. B. L. Brown received support through the following grants: Measuring the Interactive Effects of COVID‐19 and Latent Infections on Patterns of cognitive Timing: A Middle‐Aged to Elderly Sample from the Utah Valley Community. Translational Medicine Award from the BYU Simmons Research Endowment ($31,300). This is an intramural grant from a research endowment within our university to my departmental research fund. Brown also claims one patent pending Brown BL, Hedges DW, Hendrix SB (2014) Extracting A periodic Components from a Time Series Wave Data Set—Patent Pending with a filing date of October 16, 2013, Application Number PCT/US2013/065327. Data used in this manuscript is publically available. Participant consent is not needed, and was obtained from investigators who ran each study. M. L. Hardy has received paid consulting fees to the LLC solely owned by ML Hardy from ixLayer, NW Pathology & Labs, Mt. Baker Imaging, PlumCare, ADx Healthcare, Architectural Elements; received travel support from her employer for work events; and was the past Volunteer Secretary for the Lets Pool Together and Hope Philantropies non‐profit boards. P. G. Ridge had received honoraria from University of Kansas Alzheimer's Disease Center for a seminar presentation. R. R. Fortna is a board member of Northwest Pathology and ADx Healthcare and is co‐owner and holds stock for ADx healthcare. J. Nicodemus Johnson has nothing to disclose., (© 2021 The Authors. Alzheimer's & Dementia: Diagnosis, Assessment & Disease Monitoring published by Wiley Periodicals, LLC on behalf of Alzheimer's Association.)

The bluefin trevally, Caranx melampygus, also known as the bluefin kingfish or bluefin jack, is known for its remarkable, bright-blue fins. This marine teleost is a widely prized sportfish, but few resources have been devoted to the genomics and conservation of this species because it is not targeted by large-scale commercial fisheries. Population declines from recreational and artisanal overfishing have been observed in Hawai'i, USA, resulting in both an interest in aquaculture and concerns about the long-term conservation of this species. Most research to-date has been performed in Hawai'i, raising questions about the status of bluefin trevally populations across its Indo-Pacific range. Genomic resources allow for expanded research on stock status, genetic diversity, and population demography. We present a high quality, 711 Mb nuclear genome assembly of a Hawaiian bluefin trevally from noisy long-reads with a contig NG50 of 1.2 Mb and longest contig length of 8.9 Mb. As measured by single-copy orthologs, the assembly was 95% complete, and the genome is comprised of 16.9% repetitive elements. The assembly was annotated with 33.1 K protein-coding genes, 71.4% of which were assigned putative functions, using RNA-seq data from eight tissues from the same individual. This is the first whole-genome assembly published for the carangoid genus Caranx. Using this assembled genome, a multiple sequentially Markovian coalescent model was implemented to assess population demography. Estimates of effective population size suggest population expansion has occurred since the Late Pleistocene. This genome will be a valuable resource for comparative phylogenomic studies of carangoid fishes and will help elucidate demographic history and delineate stock structure for bluefin trevally populations throughout the Indo-Pacific., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America.)

Ramp sequences increase translational speed and accuracy when rare, slowly-translated codons are found at the beginnings of genes. Here, the results of the first analysis of ramp sequences in a phylogenetic construct are presented. Ramp sequences were compared from 247 vertebrates (114 Mammalian and 133 non-mammalian), where the presence and absence of ramp sequences was analyzed as a binary character in a parsimony and maximum likelihood framework. Additionally, ramp sequences were mapped to the Open Tree of Life synthetic tree to determine the number of parallelisms and reversals that occurred, and those results were compared to random permutations. Parsimony and maximum likelihood analyses of the presence and absence of ramp sequences recovered phylogenies that are highly congruent with established phylogenies. Additionally, 81% of vertebrate mammalian ramps and 81.2% of other vertebrate ramps had less parallelisms and reversals than the mean from 1000 randomly permuted trees. A chi-square analysis of completely orthologous ramp sequences resulted in a p-value < 0.001 as compared to random chance. Ramp sequences recover comparable phylogenies as other phylogenomic methods. Although not all ramp sequences appear to have a phylogenetic signal, more ramp sequences track speciation than expected by random chance. Therefore, ramp sequences may be used in conjunction with other phylogenomic approaches if many orthologs are taken into account. However, phylogenomic methods utilizing few orthologs should be cautious in incorporating ramp sequences because individual ramp sequences may provide conflicting signals.

A recent genome-wide association study (GWAS) of 59 cerebrospinal fluid (CSF) proteins with a connection to Alzheimer's disease (AD) demonstrated an association between increased levels of chemokine ligand 2 (CCL2) with an atypical chemokine receptor chemokine-binding protein 2 variant V41A (ACKR2-V41A; rs2228467). High levels of CCL2 are associated with increased risk of AD development as well as other inflammatory diseases. In this study we characterized the biological function of the ACKR2-V41A receptor compared to the wild type allele by measuring its ligand binding affinity, CCL2 scavenging efficiency, and cell activation sensitivity. We transfected Chinese hamster ovary cells with plasmids carrying wild type ACKR2 (ACKR2-WT) or the mutant ACKR2-V41A receptor. Binding affinity assays showed that ACKR2-V41A has a lower binding affinity for CCL2 and CCL4 than ACKR2-WT. CCL2 scavenging results aligned with binding affinity assays, with ACKR2-V41A cells scavenging CCL2 with a lower efficiency than ACKR2-WT. Cell activation assays also showed that ACKR2-V41A cells had significantly lower receptor upregulation (β-Arrestin-dependent signaling pathway) upon stimulation compared to ACKR2-WT cells. These findings provide molecular and biological mechanistic insights into the GWAS association of ACKR2-V41A with increased levels of CCL2 in CSF and possibly other chemokine ligands. Increased CCL2 levels are associated with accelerated cognitive decline and increased risk of AD. Understanding how this atypical chemokine receptor allele increases serum markers of inflammation could lead to novel therapeutic solutions for AD.

Identical codon pairing and co-tRNA codon pairing increase translational efficiency within genes when two codons that encode the same amino acid are translated by the same tRNA before it diffuses from the ribosome. We examine the phylogenetic signal in both identical and co-tRNA codon pairing across 23 428 species using alignment-free and parsimony methods. We determined that conserved codon pairing typically has a smaller window size than the length of a ribosome, and codon pairing tracks phylogenies across various taxonomic groups. We report a comprehensive analysis of codon pairing, including the extent to which each codon pairs. Our parsimony method generally recovers phylogenies that are more congruent with the established phylogenies than our alignment-free method. However, four of the ten taxonomic groups did not have sufficient orthologous codon pairings and were therefore analyzed using only the alignment-free methods. Since the recovered phylogenies using only codon pairing largely match phylogenies from the Open Tree of Life and the NCBI taxonomy, and are comparable to trees recovered by other algorithms, we propose that codon pairing biases are phylogenetically conserved and should be considered in conjunction with other phylogenomic techniques., Competing Interests: The authors have declared that no competing interests exist.

Using parsimony, we analyzed codon usages across 12,337 species and 25,727 orthologous genes to rank specific genes and codons according to their phylogenetic signal. We examined each codon within each ortholog to determine the codon usage for each species. In total, 890,814 codons were parsimony informative. Next, we compared species that used a codon with species that did not use the codon. We assessed each codon's congruence with species relationships provided in the Open Tree of Life (OTL) and determined the statistical probability of observing these results by random chance. We determined that 25,771 codons had no parallelisms or reversals when mapped to the OTL. Codon usages from orthologous genes spanning many species were 1109× more likely to be congruent with species relationships in the OTL than would be expected by random chance. Using the OTL as a reference, we show that codon usage is phylogenetically conserved within orthologous genes in archaea, bacteria, plants, mammals, and other vertebrates. We also show how to use our provided framework to test different tree hypotheses by confirming the placement of turtles as sister taxa to archosaurs., Competing Interests: Declaration of Competing Interest The authors declared that there is no conflict of interest., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

Carbapenem-resistant bacteria have quickly become a worldwide concern in nosocomial infections. Of the seven known carbapenemases, four have been shown to be particularly problematic: KPC, NDM, IMP, and VIM. To date, many local and species- or carbapenemase-specific epidemiological studies have been performed, which often focus on the organism itself. This report attempts to perform an inclusive (encompass both species and carbapenemase) epidemiologic study using publicly available plasmid sequences from NCBI. In this report, the gene content of these various plasmids has been characterized, replicon types of the plasmids identified, and the global spread and species promiscuity of the plasmids analyzed. Additionally, support to several groups targeting plasmid maintenance and transfer mechanisms to slow the spread of resistance plasmids is given.

Background: Plant chloroplasts and mitochondria utilize nuclear encoded proteins to replicate their DNA. These proteins are purposely built for replication in the organelle environment and are distinct from those involved in replication of the nuclear genome. These organelle-localized proteins have ancestral roots in bacterial and bacteriophage genes, supporting the endosymbiotic theory of their origin. We examined the interactions between three of these proteins from Arabidopsis thaliana: a DNA helicase-primase similar to bacteriophage T7 gp4 protein and animal mitochondrial Twinkle, and two DNA polymerases, Pol1A and Pol1B. We used a three-pronged approach to analyze the interactions, including Yeast-two-hybrid analysis, Direct Coupling Analysis (DCA), and thermophoresis., Results: Yeast-two-hybrid analysis reveals residues 120-295 of Twinkle as the minimal region that can still interact with Pol1A or Pol1B. This region is a part of the primase domain of the protein and slightly overlaps the zinc-finger and RNA polymerase subdomains located within. Additionally, we observed that Arabidopsis Twinkle interacts much more strongly with Pol1A versus Pol1B. Thermophoresis also confirms that the primase domain of Twinkle has higher binding affinity than any other region of the protein. Direct-Coupling-Analysis identified specific residues in Twinkle and the DNA polymerases critical to positive interaction between the two proteins., Conclusions: The interaction of Twinkle with Pol1A or Pol1B mimics the minimal DNA replisomes of T7 phage and those present in mammalian mitochondria. However, while T7 and mammals absolutely require their homolog of Twinkle DNA helicase-primase, Arabidopsis Twinkle mutants are seemingly unaffected by this loss. This implies that while Arabidopsis mitochondria mimic minimal replisomes from T7 and mammalian mitochondria, there is an extra level of redundancy specific to loss of Twinkle function.

Background: Common phylogenomic approaches for recovering phylogenies are often time-consuming and require annotations for orthologous gene relationships that are not always available. In contrast, alignment-free phylogenomic approaches typically use structure and oligomer frequencies to calculate pairwise distances between species. We have developed an approach to quickly calculate distances between species based on codon aversion., Methods: Utilizing a novel alignment-free character state, we present CAM, an alignment-free approach to recover phylogenies by comparing differences in codon aversion motifs (i.e., the set of unused codons within each gene) across all genes within a species. Synonymous codon usage is non-random and differs between organisms, between genes, and even within a single gene, and many genes do not use all possible codons. We report a comprehensive analysis of codon aversion within 229,742,339 genes from 23,428 species across all kingdoms of life, and we provide an alignment-free framework for its use in a phylogenetic construct. For each species, we first construct a set of codon aversion motifs spanning all genes within that species. We define the pairwise distance between two species, A and B, as one minus the number of shared codon aversion motifs divided by the total codon aversion motifs of the species, A or B, containing the fewest motifs. This approach allows us to calculate pairwise distances even when substantial differences in the number of genes or a high rate of divergence between species exists. Finally, we use neighbor-joining to recover phylogenies., Results: Using the Open Tree of Life and NCBI Taxonomy Database as expected phylogenies, our approach compares well, recovering phylogenies that largely match expected trees and are comparable to trees recovered using maximum likelihood and other alignment-free approaches. Our technique is much faster than maximum likelihood and similar in accuracy to other alignment-free approaches. Therefore, we propose that codon aversion be considered a phylogenetically conserved character that may be used in future phylogenomic studies., Availability: CAM, documentation, and test files are freely available on GitHub at https://github.com/ridgelab/cam., Competing Interests: The authors declare there are no competing interests.

Background: The human genome contains "dark" gene regions that cannot be adequately assembled or aligned using standard short-read sequencing technologies, preventing researchers from identifying mutations within these gene regions that may be relevant to human disease. Here, we identify regions with few mappable reads that we call dark by depth, and others that have ambiguous alignment, called camouflaged. We assess how well long-read or linked-read technologies resolve these regions., Results: Based on standard whole-genome Illumina sequencing data, we identify 36,794 dark regions in 6054 gene bodies from pathways important to human health, development, and reproduction. Of these gene bodies, 8.7% are completely dark and 35.2% are ≥ 5% dark. We identify dark regions that are present in protein-coding exons across 748 genes. Linked-read or long-read sequencing technologies from 10x Genomics, PacBio, and Oxford Nanopore Technologies reduce dark protein-coding regions to approximately 50.5%, 35.6%, and 9.6%, respectively. We present an algorithm to resolve most camouflaged regions and apply it to the Alzheimer's Disease Sequencing Project. We rescue a rare ten-nucleotide frameshift deletion in CR1, a top Alzheimer's disease gene, found in disease cases but not in controls., Conclusions: While we could not formally assess the association of the CR1 frameshift mutation with Alzheimer's disease due to insufficient sample-size, we believe it merits investigating in a larger cohort. There remain thousands of potentially important genomic regions overlooked by short-read sequencing that are largely resolved by long-read technologies.

Different species, genes, and locations within genes use different codons to fine-tune gene expression. Within genes, the ramp sequence assists in ribosome spacing and decreases downstream collisions by incorporating slowly-translated codons at the beginning of a gene. Although previously reported as occurring in some species, no previous attempt at extracting the ramp sequence from specific genes has been published. We present ExtRamp, a software package that quickly extracts ramp sequences from any species using the tRNA adaptation index or relative codon adaptiveness. Different filters facilitate the analysis of codon efficiency and enable identification of genes with a ramp sequence. We validate the existence of a ramp sequence in most species by running ExtRamp on 229 742 339 genes across 23 428 species. We evaluate differences in reported ramp sequences when we use different parameters. Using the strictest ramp sequence cut-off, we show that across most taxonomic groups, ramp sequences are approximately 20-40 codons long and occur in about 10% of gene sequences. We also show that in Drosophila melanogaster as gene expression increases, a higher proportion of genes have ramp sequences. We provide a framework for performing this analysis on other species. ExtRamp is freely available at https://github.com/ridgelab/ExtRamp., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Motivation: Orthologous gene identification is fundamental to all aspects of biology. For example, ortholog identification between species can provide functional insights for genes of unknown function and is a necessary step in phylogenetic inference. Currently, most ortholog identification algorithms require all-versus-all BLAST comparisons, which are time-consuming and memory intensive., Results: In contrast to existing approaches, JustOrthologs exploits the conservation of gene structure by using the lengths of coding sequence regions and dinucleotide percentages to identify orthologs. In comparison to OrthoMCL, OMA and OrthoFinder, JustOrthologs decreases ortholog identification runtime by more than 96% and achieves comparable precision and recall scores. The computational speedup allowed us to conduct pairwise comparisons of 1197 complete genomes (780 eukaryotes and 417 archaea). We confirmed gene annotations for 384 120 genes, grouped 1 675 415 genes in previously unreported ortholog groups, and identified 51 429 potentially mislabeled genes across 622 843 ortholog groups., Availability and Implementation: JustOrthologs is an open source collaborative software package available in the GitHub repository: https://github.com/ridgelab/JustOrthologs/. All test FASTA files used for comparisons are freely available at https://github.com/ridgelab/JustOrthologs/comparisonFastaFiles/. Reference genomes used in this work are available for download from the NCBI repository: ftp://ftp.ncbi.nih.gov/genomes/., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2018. Published by Oxford University Press.)

Polygenic scores (or genetic risk scores) quantify the aggregate of small effects from many common genetic loci that have been associated with a trait through genome-wide association. Polygenic scores were first used successfully in schizophrenia and have since been applied to multiple phenotypes including multiple sclerosis, rheumatoid arthritis, and height. Because human height is an easily-measured and complex polygenic trait, polygenic height scores provide exciting insights into the predictability of aggregate common variant effect on the phenotype. Shawn Bradley is an extremely tall former professional basketball player from Brigham Young University and the National Basketball Association (NBA), measuring 2.29 meters (7'6 ″ , 99.99999th percentile for height) tall, with no known medical conditions. Here, we present a case where a rare combination of common SNPs in one individual results in an extremely high polygenic height score that is correlated with an extreme phenotype. While polygenic scores are not clinically significant in the average case, our findings suggest that for extreme phenotypes, polygenic scores may be more successful for the prediction of individuals.

Mendelian adult-onset leukodystrophies are a spectrum of rare inherited progressive neurodegenerative disorders affecting the white matter of the central nervous system. Among these, cerebral autosomal dominant and recessive arteriopathy with subcortical infarcts and leukoencephalopathy, cerebroretinal vasculopathy, metachromatic leukodystrophy, hereditary diffuse leukoencephalopathy with spheroids, and vanishing white matter disease present with rapidly progressive dementia as dominant feature and are caused by mutations in NOTCH3, HTRA1, TREX1, ARSA, CSF1R, EIF2B1, EIF2B2, EIF2B3, EIF2B4, and EIF2B5, respectively. Given the rare incidence of these disorders and the lack of unequivocally diagnostic features, leukodystrophies are frequently misdiagnosed with common sporadic dementing diseases such as Alzheimer's disease (AD), raising the question of whether these overlapping phenotypes may be explained by shared genetic risk factors. To investigate this intriguing hypothesis, we have combined gene expression analysis (1) in 6 different AD mouse strains (APPPS1, HOTASTPM, HETASTPM, TPM, TAS10, and TAU) at 5 different developmental stages (embryo [E15], 2, 4, 8, and 18 months), (2) in APPPS1 primary cortical neurons under stress conditions (oxygen-glucose deprivation) and single-variant-based and single-gene-based (c-alpha test and sequence kernel association test (SKAT)) genetic screening in a cohort composed of 332 Caucasian late-onset AD patients and 676 Caucasian elderly controls. Csf1r was significantly overexpressed (log2FC > 1, adj. p-value < 0.05) in the cortex and hippocampus of aged HOTASTPM mice with extensive Aβ dense-core plaque pathology. We identified 3 likely pathogenic mutations in CSF1R TK domain (p.L868R, p.Q691H, and p.H703Y) in our discovery and validation cohort, composed of 465 AD and mild cognitive impairment (MCI) Caucasian patients from the United Kingdom. Moreover, NOTCH3 was a significant hit in the c-alpha test (adj p-value = 0.01). Adult-onset Mendelian leukodystrophy genes are not common factors implicated in AD. Nevertheless, our study suggests a potential pathogenic link between NOTCH3, CSF1R, and sporadic late-onset AD, which warrants further investigation., (Copyright © 2018 The Authors. Published by Elsevier Inc. All rights reserved.)

Introduction: Mitochondrial genetics are an important but largely neglected area of research in Alzheimer's disease. A major impediment is the lack of data sets., Methods: We used an innovative, rigorous approach, combining several existing tools with our own, to accurately assemble and call variants in 809 whole mitochondrial genomes., Results: To help address this impediment, we prepared a data set that consists of 809 complete and annotated mitochondrial genomes with samples from the Alzheimer's Disease Neuroimaging Initiative. These whole mitochondrial genomes include rich phenotyping, such as clinical, fluid biomarker, and imaging data, all of which is available through the Alzheimer's Disease Neuroimaging Initiative website. Genomes are cleaned, annotated, and prepared for analysis., Discussion: These data provide an important resource for investigating the impact of mitochondrial genetic variation on risk for Alzheimer's disease and other phenotypes that have been measured in the Alzheimer's Disease Neuroimaging Initiative samples., (Copyright © 2017 The Authors. Published by Elsevier Inc. All rights reserved.)

Correction: The original version of this article [1] unfortunately contained a typographical error. The 'Alzheimer's Disease Neuroimaging Initiative' was erroneously included as 'Alzheimer's Disease Neuroimaging Initative' in the author list of the article.

Purpose of Review: Alzheimer's disease (AD) is the most common form of dementia, affects an increasing number of people worldwide, has a rapidly increasing incidence, and is fatal. In the past several years, significant progress has been made towards solving the genetic architecture of AD, but our understanding remains incomplete and has not led to treatments that either cure or slow disease. There is substantial evidence that mitochondria are involved in AD: mitochondrial functional declines in AD, mitochondrial encoded gene expression changes, mitochondria are morphologically different, and mitochondrial fusion/fission are modified. While a majority of mitochondrial proteins are nuclear encoded and could lead to malfunction in mitochondria, the mitochondrial genome encodes numerous proteins important for the electron transport chain, which if damaged could possibly lead to mitochondrial changes observed in AD. Here, we review publications that describe a relationship between the mitochondrial genome and AD and make suggestions for analysis approaches and data acquisition, from existing datasets, to study the mitochondrial genetics of AD., Recent Findings: Numerous mitochondrial haplogroups and SNPs have been reported to influence risk for AD, but the majority of these have not been replicated, nor experimentally validated., Summary: The role of the mitochondrial genome in AD remains elusive, and several impediments exist to fully understand the relationship between the mitochondrial genome and AD. Yet, by leveraging existing datasets and implementing appropriate analysis approaches, determining the role of mitochondrial genetics in risk for AD is possible., Competing Interests: Compliance with Ethical StandardsThe authors declare that they have no conflict of interest.This article does not contain any studies with human or animal subjects performed by any of the authors.

Motivation: One of the main challenges with bioinformatics software is that the size and complexity of datasets necessitate trading speed for accuracy, or completeness. To combat this problem of computational complexity, a plethora of heuristic algorithms have arisen that report a 'good enough' solution to biological questions. However, in instances such as Simple Sequence Repeats (SSRs), a 'good enough' solution may not accurately portray results in population genetics, phylogenetics and forensics, which require accurate SSRs to calculate intra- and inter-species interactions., Results: We present Kmer-SSR, which finds all SSRs faster than most heuristic SSR identification algorithms in a parallelized, easy-to-use manner. The exhaustive Kmer-SSR option has 100% precision and 100% recall and accurately identifies every SSR of any specified length. To identify more biologically pertinent SSRs, we also developed several filters that allow users to easily view a subset of SSRs based on user input. Kmer-SSR, coupled with the filter options, accurately and intuitively identifies SSRs quickly and in a more user-friendly manner than any other SSR identification algorithm., Availability and Implementation: The source code is freely available on GitHub at https://github.com/ridgelab/Kmer-SSR., Contact: perry.ridge@byu.edu., (© The Author(s) 2017. Published by Oxford University Press.)

Background: While age and the APOE ε4 allele are major risk factors for Alzheimer's disease (AD), a small percentage of individuals with these risk factors exhibit AD resilience by living well beyond 75 years of age without any clinical symptoms of cognitive decline., Methods: We used over 200 "AD resilient" individuals and an innovative, pedigree-based approach to identify genetic variants that segregate with AD resilience. First, we performed linkage analyses in pedigrees with resilient individuals and a statistical excess of AD deaths. Second, we used whole genome sequences to identify candidate SNPs in significant linkage regions. Third, we replicated SNPs from the linkage peaks that reduced risk for AD in an independent dataset and in a gene-based test. Finally, we experimentally characterized replicated SNPs., Results: Rs142787485 in RAB10 confers significant protection against AD (p value = 0.0184, odds ratio = 0.5853). Moreover, we replicated this association in an independent series of unrelated individuals (p value = 0.028, odds ratio = 0.69) and used a gene-based test to confirm a role for RAB10 variants in modifying AD risk (p value = 0.002). Experimentally, we demonstrated that knockdown of RAB10 resulted in a significant decrease in Aβ42 (p value = 0.0003) and in the Aβ42/Aβ40 ratio (p value = 0.0001) in neuroblastoma cells. We also found that RAB10 expression is significantly elevated in human AD brains (p value = 0.04)., Conclusions: Our results suggest that RAB10 could be a promising therapeutic target for AD prevention. In addition, our gene discovery approach can be expanded and adapted to other phenotypes, thus serving as a model for future efforts to identify rare variants for AD and other complex human diseases.

Although many studies have documented codon usage bias in different species, the importance of codon usage in a phylogenetic framework remains largely unknown. We demonstrate that a phylogenetic signal is present in the codon usage and non-usage biases of 17 717 orthologues evaluated across 72 tetrapod species using a simple parsimony analysis of a binary matrix of codon characters. Phylogenies estimated using stop codons were more congruent with previous hypotheses than phylogenies based on any other single codon or a combination of codons. Although each codon is present in every species, specific genes have different codon preferences and may or may not use every possible codon. This observation allowed us to map the pattern of codon usage and non-usage across the topology. These results suggest that codon usage is phylogenetically conserved across shallow and deep levels within tetrapods., (© The Willi Hennig Society 2017.)

Genome-wide association studies (GWASs) have been effective approaches to dissect common genetic variability underlying complex diseases in a systematic and unbiased way. Recently, GWASs have led to the discovery of over 20 susceptibility loci for Alzheimer's disease (AD). Despite the evidence showing the contribution of these loci to AD pathogenesis, their genetic architecture has not been extensively investigated, leaving the possibility that low frequency and rare coding variants may also occur and contribute to the risk of disease. We have used exome and genome sequencing data to analyze the single independent and joint effect of rare and low-frequency protein coding variants in 9 AD GWAS loci with the strongest effect sizes after APOE (BIN1, CLU, CR1, PICALM, MS4A6A, ABCA7, EPHA1, CD33, and CD2AP) in a cohort of 332 sporadic AD cases and 676 elderly controls of British and North-American ancestry. We identified coding variability in ABCA7 as contributing to AD risk. This locus harbors a low-frequency coding variant (p.G215S, rs72973581, minor allele frequency = 4.3%) conferring a modest but statistically significant protection against AD (p-value = 0.024, odds ratio = 0.57, 95% confidence interval = 0.41-0.80). Notably, our results are not driven by an enrichment of loss of function variants in ABCA7, recently reported as main pathogenic factor underlying AD risk at this locus. In summary, our study confirms the role of ABCA7 in AD and provides new insights that should address functional studies., (Copyright © 2016 The Authors. Published by Elsevier Inc. All rights reserved.)

Unlabelled: Simple Sequence Repeats (SSRs) are used to address a variety of research questions in a variety of fields (e.g. population genetics, phylogenetics, forensics, etc.), due to their high mutability within and between species. Here, we present an innovative algorithm, SA-SSR, based on suffix and longest common prefix arrays for efficiently detecting SSRs in large sets of sequences. Existing SSR detection applications are hampered by one or more limitations (i.e. speed, accuracy, ease-of-use, etc.). Our algorithm addresses these challenges while being the most comprehensive and correct SSR detection software available. SA-SSR is 100% accurate and detected >1000 more SSRs than the second best algorithm, while offering greater control to the user than any existing software., Availability and Implementation: SA-SSR is freely available at http://github.com/ridgelab/SA-SSR CONTACT: perry.ridge@byu.edu, Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author 2016. Published by Oxford University Press.)

Background: Genome-wide association studies (GWAS) have effectively identified genetic factors for many diseases. Many diseases, including Alzheimer's disease (AD), have epistatic causes, requiring more sophisticated analyses to identify groups of variants which together affect phenotype., Results: Based on the GWAS statistical model, we developed a multi-SNP GWAS analysis to identify pairs of variants whose common occurrence signaled the Alzheimer's disease phenotype., Conclusions: Despite not having sufficient data to demonstrate significance, our preliminary experimentation identified a high correlation between GRIA3 and HLA-DRB5 (an AD gene). GRIA3 has not been previously reported in association with AD, but is known to play a role in learning and memory.

Background: Analyzing next-generation sequencing data is difficult because datasets are large, second generation sequencing platforms have high error rates, and because each position in the target genome (exome, transcriptome, etc.) is sequenced multiple times. Given these challenges, numerous bioinformatic algorithms have been developed to analyze these data. These algorithms aim to find an appropriate balance between data loss, errors, analysis time, and memory footprint. Typical analysis pipelines require multiple steps. If one or more of these steps is unnecessary, it would significantly decrease compute time and data manipulation to remove the step. One step in many pipelines is PCR duplicate removal, where PCR duplicates arise from multiple PCR products from the same template molecule binding on the flowcell. These are often removed because there is concern they can lead to false positive variant calls. Picard (MarkDuplicates) and SAMTools (rmdup) are the two main softwares used for PCR duplicate removal., Results: Approximately 92 % of the 17+ million variants called were called whether we removed duplicates with Picard or SAMTools, or left the PCR duplicates in the dataset. There were no significant differences between the unique variant sets when comparing the transition/transversion ratios (p = 1.0), percentage of novel variants (p = 0.99), average population frequencies (p = 0.99), and the percentage of protein-changing variants (p = 1.0). Results were similar for variants in the American College of Medical Genetics genes. Genotype concordance between NGS and SNP chips was above 99 % for all genotype groups (e.g., homozygous reference)., Conclusions: Our results suggest that PCR duplicate removal has minimal effect on the accuracy of subsequent variant calls.

Background: Prostatic Acid Phosphatase (PAP) is an enzyme that is produced primarily in the prostate and functions as a cell growth regulator and potential tumor suppressor. Understanding the genetic regulation of this enzyme is important because PAP plays an important role in prostate cancer and is expressed in other tissues such as the brain., Methods: We tested association between 5.8 M SNPs and PAP levels in cerebrospinal fluid across 543 individuals in two datasets using linear regression. We then performed meta-analyses using METAL =with a significance threshold of p < 5 × 10(-8) and removed SNPs where the direction of the effect was different between the two datasets, identifying 289 candidate SNPs that affect PAP cerebrospinal fluid levels. We analyzed each of these SNPs individually and prioritized SNPs that had biologically meaningful functional annotations in wANNOVAR (e.g. non-synonymous, stop gain, 3' UTR, etc.) or had a RegulomeDB score less than 3., Results: Thirteen SNPs met our criteria, suggesting they are candidate causal alleles that underlie ACPP regulation and expression., Conclusions: Given PAP's expression in the brain and its role as a cell-growth regulator and tumor suppressor, our results have important implications in brain health such as cancer and other brain diseases including neurodegenerative diseases (e.g., Alzheimer's disease and Parkinson's disease) and mental health (e.g., anxiety, depression, and schizophrenia).

Background: Prolactin is a polypeptide hormone secreted by the anterior pituitary gland that plays an essential role in lactation, tissue growth, and suppressing apoptosis to increase cell survival. Prolactin serves as a key player in many life-critical processes, including immune system and reproduction. Prolactin is also found in multiple fluids throughout the body, including plasma and cerebrospinal fluid (CSF)., Methods: In this study, we measured prolactin levels in both plasma and CSF, and performed a genome-wide association study. We then performed meta-analyses using METAL with a significance threshold of p < 5 × 10(-8) and removed SNPs where the direction of the effect was different between the two datasets., Results: We identified 12 SNPs associated with increased prolactin levels in both biological fluids., Conclusions: Our efforts will help researchers understand how prolactin is regulated in both CSF and plasma, which could be beneficial in research for the immune system and reproduction.

Background: CCL16 is a chemokine predominantly expressed in the liver, but is also found in the blood and brain, and is known to play important roles in immune response and angiogenesis. Little is known about the gene's regulation., Methods: Here, we test for potential causal SNPs that affect CCL16 protein levels in both blood plasma and cerebrospinal fluid in a genome-wide association study across two datasets. We then use METAL to performed meta-analyses with a significance threshold of p < 5x10(-8). We removed SNPs where the direction of the effect was different between the two datasets., Results: We identify 10 SNPs associated with increased CCL16 protein levels in both biological fluids., Conclusions: Our results will help understand CCL16's regulation, allowing researchers to better understand the gene's effects on human health.

The cerebral deposition of Aβ42, a neurotoxic proteolytic derivate of amyloid precursor protein (APP), is a central event in Alzheimer's disease (AD)(Amyloid hypothesis). Given the key role of APP-Aβ metabolism in AD pathogenesis, we selected 29 genes involved in APP processing, Aβ degradation and clearance. We then used exome and genome sequencing to investigate the single independent (single-variant association test) and cumulative (gene-based association test) effect of coding variants in these genes as potential susceptibility factors for AD, in a cohort composed of 332 sporadic and mainly late-onset AD cases and 676 elderly controls from North America and the UK. Our study shows that common coding variability in these genes does not play a major role for the disease development. In the single-variant association analysis, the main hits, none of which statistically significant after multiple testing correction (1.9e-4

Alzheimer's disease (AD) is a complex genetic disorder with no effective treatments. More than 20 common markers have been identified, which are associated with AD. Recently, several rare variants have been identified in Amyloid Precursor Protein (APP), Triggering Receptor Expressed On Myeloid Cells 2 (TREM2) and Unc-5 Netrin Receptor C (UNC5C) that affect risk for AD. Despite the many successes, the genetic architecture of AD remains unsolved. We used Genome-wide Complex Trait Analysis to (1) estimate phenotypic variance explained by genetics; (2) calculate genetic variance explained by known AD single nucleotide polymorphisms (SNPs); and (3) identify the genomic locations of variation that explain the remaining unexplained genetic variance. In total, 53.24% of phenotypic variance is explained by genetics, but known AD SNPs only explain 30.62% of the genetic variance. Of the unexplained genetic variance, approximately 41% is explained by unknown SNPs in regions adjacent to known AD SNPs, and the remaining unexplained genetic variance outside these regions., (Copyright © 2016 The Authors. Published by Elsevier Inc. All rights reserved.)

Introduction: Ebbert et al. reported gene-gene interactions between rs11136000-rs670139 (CLU-MS4A4E) and rs3865444-rs670139 (CD33-MS4A4E). We evaluate these interactions in the largest data set for an epistasis study., Methods: We tested interactions using 3837 cases and 4145 controls from Alzheimer's Disease Genetics Consortium using meta-analyses and permutation analyses. We repeated meta-analyses stratified by apolipoprotein E (APOE) ε4 status, estimated combined odds ratio (OR) and population attributable fraction (cPAF), and explored causal variants., Results: Results support the CLU-MS4A4E interaction and a dominant effect. An association between CLU-MS4A4E and APOE ε4 negative status exists. The estimated synergy factor, OR, and cPAF for rs11136000-rs670139 are 2.23, 2.45, and 8.0, respectively. We identified potential causal variants., Discussion: We replicated the CLU-MS4A4E interaction in a large case-control series and observed APOE ε4 and possible dominant effect. The CLU-MS4A4E OR is higher than any Alzheimer's disease locus except APOE ε4, APP, and TREM2. We estimated an 8% decrease in Alzheimer's disease incidence without CLU-MS4A4E risk alleles and identified potential causal variants., (Copyright © 2016 The Alzheimer's Association. Published by Elsevier Inc. All rights reserved.)

Genome-wide association studies of 146 plasma protein levels in 818 individuals revealed 56 genome-wide significant associations (28 novel) with 47 analytes. Loci associated with plasma levels of 39 proteins tested have been previously associated with various complex traits such as heart disease, inflammatory bowel disease, Type 2 diabetes and multiple sclerosis. These data suggest that these plasma protein levels may constitute informative endophenotypes for these complex traits. We found three potential pleiotropic genes: ABO for plasma SELE and ACE levels, FUT2 for CA19-9 and CEA plasma levels and APOE for ApoE and CRP levels. We also found multiple independent signals in loci associated with plasma levels of ApoH, CA19-9, FetuinA, IL6r and LPa. Our study highlights the power of biological traits for genetic studies to identify genetic variants influencing clinically relevant traits, potential pleiotropic effects and complex disease associations in the same locus., (© 2016. The Author(s).)

Cholesterol has been implicated in the pathogenesis of late-onset Alzheimer's disease (LOAD) and the cholesteryl ester transfer protein (CETP) is critical to cholesterol regulation within the cell, making CETP an Alzheimer's disease candidate gene. Several studies have suggested that CETP I405V (rs5882) is associated with cognitive function and LOAD risk, but findings vary and most studies have been conducted using relatively small numbers of samples. To test whether this variant is involved in cognitive function and LOAD progression, we genotyped 4486 subjects with up to 12 years of longitudinal cognitive assessment. Analyses revealed an average 0.6-point decrease per year in the rate of cognitive decline for each additional valine (p < 0.011). We failed to detect the association between CETP I405V and LOAD status (p < 0.28). We conclude that CETP I405V is associated with preserved cognition over time but is not associated with LOAD status., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Background: Proper phosphate signaling is essential for robust growth of Escherichia coli and many other bacteria. The phosphate signal is mediated by a classic two component signal system composed of PhoR and PhoB. The PhoR histidine kinase is responsible for phosphorylating/dephosphorylating the response regulator, PhoB, which controls the expression of genes that aid growth in low phosphate conditions. The mechanism by which PhoR receives a signal of environmental phosphate levels has remained elusive. A transporter complex composed of the PstS, PstC, PstA, and PstB proteins as well as a negative regulator, PhoU, have been implicated in signaling environmental phosphate to PhoR., Results: This work confirms that PhoU and the PstSCAB complex are necessary for proper signaling of high environmental phosphate. Also, we identify residues important in PhoU/PhoR interaction with genetic analysis. Using protein modeling and docking methods, we show an interaction model that points to a potential mechanism for PhoU mediated signaling to PhoR to modify its activity. This model is tested with direct coupling analysis., Conclusions: These bioinformatics tools, in combination with genetic and biochemical analysis, help to identify and test a model for phosphate signaling and may be applicable to several other systems.

Alzheimer's disease affects millions of people worldwide and incidence is expected to rise as the population ages, but no effective therapies exist despite decades of research and more than 20 known disease markers. Research has shown that Alzheimer's disease's missing heritability remains extensive with an estimated 25% of phenotypic variance unexplained by known variants. The missing heritability may be explained by missing variants or by epistasis. Researchers often focus on individual loci rather than epistatic interactions, which is likely an oversimplification of the underlying biology since most phenotypes are affected by multiple genes. Focusing research efforts on epistasis will be critical to resolving Alzheimer's disease etiology, and a major key to identifying and properly interpreting key epistatic interactions will be bridging the gap between statistical and biological epistasis. This review covers the current state of epistasis research in Alzheimer's disease and how researchers can bridge the gap between statistical and biological epistasis to help resolve Alzheimer's disease etiology.

Cerebrospinal fluid (CSF) 42 amino acid species of amyloid beta (Aβ42) and tau levels are strongly correlated with the presence of Alzheimer's disease (AD) neuropathology including amyloid plaques and neurodegeneration and have been successfully used as endophenotypes for genetic studies of AD. Additional CSF analytes may also serve as useful endophenotypes that capture other aspects of AD pathophysiology. Here we have conducted a genome-wide association study of CSF levels of 59 AD-related analytes. All analytes were measured using the Rules Based Medicine Human DiscoveryMAP Panel, which includes analytes relevant to several disease-related processes. Data from two independently collected and measured datasets, the Knight Alzheimer's Disease Research Center (ADRC) and Alzheimer's Disease Neuroimaging Initiative (ADNI), were analyzed separately, and combined results were obtained using meta-analysis. We identified genetic associations with CSF levels of 5 proteins (Angiotensin-converting enzyme (ACE), Chemokine (C-C motif) ligand 2 (CCL2), Chemokine (C-C motif) ligand 4 (CCL4), Interleukin 6 receptor (IL6R) and Matrix metalloproteinase-3 (MMP3)) with study-wide significant p-values (p<1.46×10-10) and significant, consistent evidence for association in both the Knight ADRC and the ADNI samples. These proteins are involved in amyloid processing and pro-inflammatory signaling. SNPs associated with ACE, IL6R and MMP3 protein levels are located within the coding regions of the corresponding structural gene. The SNPs associated with CSF levels of CCL4 and CCL2 are located in known chemokine binding proteins. The genetic associations reported here are novel and suggest mechanisms for genetic control of CSF and plasma levels of these disease-related proteins. Significant SNPs in ACE and MMP3 also showed association with AD risk. Our findings suggest that these proteins/pathways may be valuable therapeutic targets for AD. Robust associations in cognitively normal individuals suggest that these SNPs also influence regulation of these proteins more generally and may therefore be relevant to other diseases.

Background: Reported odds ratios and population attributable fractions (PAF) for late-onset Alzheimer's disease (LOAD) risk loci (BIN1, ABCA7, CR1, MS4A4E, CD2AP, PICALM, MS4A6A, CD33, and CLU) come from clinically ascertained samples. Little is known about the combined PAF for these LOAD risk alleles and the utility of these combined markers for case-control prediction. Here we evaluate these loci in a large population-based sample to estimate PAF and explore the effects of additive and nonadditive interactions on LOAD status prediction performance., Methods: 2419 samples from the Cache County Memory Study were genotyped for APOE and nine LOAD risk loci from AlzGene.org. We used logistic regression and receiver operator characteristic analysis to assess the LOAD status prediction performance of these loci using additive and nonadditive models and compared odds ratios and PAFs between AlzGene.org and Cache County., Results: Odds ratios were comparable between Cache County and AlzGene.org when identical single nucleotide polymorphisms were genotyped. PAFs from AlzGene.org ranged from 2.25% to 37%; those from Cache County ranged from .05% to 20%. Including non-APOE alleles significantly improved LOAD status prediction performance (area under the curve = .80) over APOE alone (area under the curve = .78) when not constrained to an additive relationship (p < .03). We identified potential allelic interactions (p values uncorrected): CD33-MS4A4E (synergy factor = 5.31; p < .003) and CLU-MS4A4E (synergy factor = 3.81; p < .016)., Conclusions: Although nonadditive interactions between loci significantly improve diagnostic ability, the improvement does not reach the desired sensitivity or specificity for clinical use. Nevertheless, these results suggest that understanding gene-gene interactions may be important in resolving Alzheimer's disease etiology., (Copyright © 2014 Society of Biological Psychiatry. Published by Elsevier Inc. All rights reserved.)

Genome-wide association studies (GWAS) have identified several risk variants for late-onset Alzheimer's disease (LOAD). These common variants have replicable but small effects on LOAD risk and generally do not have obvious functional effects. Low-frequency coding variants, not detected by GWAS, are predicted to include functional variants with larger effects on risk. To identify low-frequency coding variants with large effects on LOAD risk, we carried out whole-exome sequencing (WES) in 14 large LOAD families and follow-up analyses of the candidate variants in several large LOAD case-control data sets. A rare variant in PLD3 (phospholipase D3; Val232Met) segregated with disease status in two independent families and doubled risk for Alzheimer's disease in seven independent case-control series with a total of more than 11,000 cases and controls of European descent. Gene-based burden analyses in 4,387 cases and controls of European descent and 302 African American cases and controls, with complete sequence data for PLD3, reveal that several variants in this gene increase risk for Alzheimer's disease in both populations. PLD3 is highly expressed in brain regions that are vulnerable to Alzheimer's disease pathology, including hippocampus and cortex, and is expressed at significantly lower levels in neurons from Alzheimer's disease brains compared to control brains. Overexpression of PLD3 leads to a significant decrease in intracellular amyloid-β precursor protein (APP) and extracellular Aβ42 and Aβ40 (the 42- and 40-residue isoforms of the amyloid-β peptide), and knockdown of PLD3 leads to a significant increase in extracellular Aβ42 and Aβ40. Together, our genetic and functional data indicate that carriers of PLD3 coding variants have a twofold increased risk for LOAD and that PLD3 influences APP processing. This study provides an example of how densely affected families may help to identify rare variants with large effects on risk for disease or other complex traits.

Background: Since the advent of next-generation sequencing many previously untestable hypotheses have been realized. Next-generation sequencing has been used for a wide range of studies in diverse fields such as population and medical genetics, phylogenetics, microbiology, and others. However, this novel technology has created unanticipated challenges such as the large numbers of genetic variants. Each caucasian genome has more than four million single nucleotide variants, insertions and deletions, copy number variants, and structural variants. Several formats have been suggested for storing these variants; however, the variant call format (VCF) has become the community standard., Results: We developed new software called the Variant Tool Chest (VTC) to provide much needed tools to work with VCF files. VTC provides a variety of tools for manipulating, comparing, and analyzing VCF files beyond the functionality of existing tools. In addition, VTC was written to be easily extended with new tools., Conclusions: Variant Tool Chest brings new and important functionality that complements and integrates well with existing software. VTC is available at https://github.com/mebbert/VariantToolChest.

Background: The mitochondria are essential organelles and are the location of cellular respiration, which is responsible for the majority of ATP production. Each cell contains multiple mitochondria, and each mitochondrion contains multiple copies of its own circular genome. The ratio of mitochondrial genomes to nuclear genomes is referred to as mitochondrial copy number. Decreases in mitochondrial copy number are known to occur in many tissues as people age, and in certain diseases. The regulation of mitochondrial copy number by nuclear genes has been studied extensively. While mitochondrial variation has been associated with longevity and some of the diseases known to have reduced mitochondrial copy number, the role that the mitochondrial genome itself has in regulating mitochondrial copy number remains poorly understood., Results: We analyzed the complete mitochondrial genomes from 1007 individuals randomly selected from the Cache County Study on Memory Health and Aging utilizing the inferred evolutionary history of the mitochondrial haplotypes present in our dataset to identify sequence variation and mitochondrial haplotypes associated with changes in mitochondrial copy number. Three variants belonging to mitochondrial haplogroups U5A1 and T2 were significantly associated with higher mitochondrial copy number in our dataset., Conclusions: We identified three variants associated with higher mitochondrial copy number and suggest several hypotheses for how these variants influence mitochondrial copy number by interacting with known regulators of mitochondrial copy number. Our results are the first to report sequence variation in the mitochondrial genome that causes changes in mitochondrial copy number. The identification of these variants that increase mtDNA copy number has important implications in understanding the pathological processes that underlie these phenotypes.

Background: Population stratification is a key concern for genetic association analyses. In addition, extreme homogeneity of ethnic origins of a population can make it difficult to interpret how genetic associations in that population may translate into other populations. Here we have evaluated the genetic substructure of samples from the Cache County study relative to the HapMap Reference populations and data from the Alzheimer's Disease Neuroimaging Initiative (ADNI)., Results: Our findings show that the Cache County study is similar in ethnic diversity to the self-reported "Whites" in the ADNI sample and less homogenous than the HapMap CEU population., Conclusions: We conclude that the Cache County study is genetically representative of the general European American population in the USA and is an appropriate population for conducting broadly applicable genetic studies.

Common variable immunodeficiency (CVID) is a heterogeneous disorder characterized by antibody deficiency, poor humoral response to antigens, and recurrent infections. To investigate the molecular cause of CVID, we carried out exome sequence analysis of a family diagnosed with CVID and identified a heterozygous frameshift mutation, c.2564delA (p.Lys855Serfs(∗)7), in NFKB2 affecting the C terminus of NF-κB2 (also known as p100/p52 or p100/p49). Subsequent screening of NFKB2 in 33 unrelated CVID-affected individuals uncovered a second heterozygous nonsense mutation, c.2557C>T (p.Arg853(∗)), in one simplex case. Affected individuals in both families presented with an unusual combination of childhood-onset hypogammaglobulinemia with recurrent infections, autoimmune features, and adrenal insufficiency. NF-κB2 is the principal protein involved in the noncanonical NF-κB pathway, is evolutionarily conserved, and functions in peripheral lymphoid organ development, B cell development, and antibody production. In addition, Nfkb2 mouse models demonstrate a CVID-like phenotype with hypogammaglobulinemia and poor humoral response to antigens. Immunoblot analysis and immunofluorescence microscopy of transformed B cells from affected individuals show that the NFKB2 mutations affect phosphorylation and proteasomal processing of p100 and, ultimately, p52 nuclear translocation. These findings describe germline mutations in NFKB2 and establish the noncanonical NF-κB signaling pathway as a genetic etiology for this primary immunodeficiency syndrome., (Copyright © 2013 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)