Genetic and multi-omic risk assessment of Alzheimer's disease implicates core associated biological domains.

INTRODUCTION: Alzheimer's disease (AD) is the predominant dementia globally, with heterogeneous presentation and penetrance of clinical symptoms, variable presence of mixed pathologies, potential disease subtypes, and numerous associated endophenotypes. Beyond the difficulty of designing treatments that address the core pathological characteristics of the disease, therapeutic development is challenged by the uncertainty of which endophenotypic areas and specific targets implicated by those endophenotypes to prioritize for further translational research. However, publicly funded consortia driving large-scale open science efforts have produced multiple omic analyses that address both disease risk relevance and biological process involvement of genes across the genome. METHODS: Here we report the development of an informatic pipeline that draws from genetic association studies, predicted variant impact, and linkage with dementia associated phenotypes to create a genetic risk score. This is paired with a multi-omic risk score utilizing extensive sets of both transcriptomic and proteomic studies to identify system-level changes in expression associated with AD. These two elements combined constitute our target risk score that ranks AD risk genome-wide. The ranked genes are organized into endophenotypic space through the development of 19 biological domains associated with AD in the described genetics and genomics studies and accompanying literature. The biological domains are constructed from exhaustive Gene Ontology (GO) term compilations, allowing automated assignment of genes into objectively defined disease-associated biology. This rank-and-organize approach, performed genome-wide, allows the characterization of aggregations of AD risk across biological domains. RESULTS: The top AD-risk-associated biological domains are Synapse, Immune Response, Lipid Metabolism, Mitochondrial Metabolism, Structural Stabilization, and Proteostasis, with slightly lower levels of risk enrichment present within the other 13 biological domains. DISCUSSION: This provides an objective methodology to localize risk within specific biological endophenotypes and drill down into the most significantly associated sets of GO terms and annotated genes for potential therapeutic targets.

Genome-wide analysis of common and rare variants via multiple knockoffs at biobank scale, with an application to Alzheimer disease genetics.

Knockoff-based methods have become increasingly popular due to their enhanced power for locus discovery and their ability to prioritize putative causal variants in a genome-wide analysis. However, because of the substantial computational cost for generating knockoffs, existing knockoff approaches cannot analyze millions of rare genetic variants in biobank-scale whole-genome sequencing and whole-genome imputed datasets. We propose a scalable knockoff-based method for the analysis of common and rare variants across the genome, KnockoffScreen-AL, that is applicable to biobank-scale studies with hundreds of thousands of samples and millions of genetic variants. The application of KnockoffScreen-AL to the analysis of Alzheimer disease (AD) in 388,051 WG-imputed samples from the UK Biobank resulted in 31 significant loci, including 14 loci that are missed by conventional association tests on these data. We perform replication studies in an independent meta-analysis of clinically diagnosed AD with 94,437 samples, and additionally leverage single-cell RNA-sequencing data with 143,793 single-nucleus transcriptomes from 17 control subjects and AD-affected individuals, and proteomics data from 735 control subjects and affected indviduals with AD and related disorders to validate the genes at these significant loci. These multi-omics analyses show that 79.1% of the proximal genes at these loci and 76.2% of the genes at loci identified only by KnockoffScreen-AL exhibit at least suggestive signal (p < 0.05) in the scRNA-seq or proteomics analyses. We highlight a potentially causal gene in AD progression, EGFR, that shows significant differences in expression and protein levels between AD-affected individuals and healthy control subjects.

Sex-specific nicotine sensitization and imprinting of self-administration in rats inform GWAS findings on human addiction phenotypes.

Repeated nicotine exposure leads to sensitization (SST) and enhances self-administration (SA) in rodents. However, the molecular basis of nicotine SST and SA and their biological relevance to the mounting genome-wide association study (GWAS) loci of human addictive behaviors are poorly understood. Considering a gateway drug role of nicotine, we modeled nicotine SST and SA in F1 progeny of inbred rats (F344/BN) and conducted integrative genomics analyses. We unexpectedly observed male-specific nicotine SST and a parental effect of SA only present in paternal F344 crosses. Transcriptional profiling in the ventral tegmental area (VTA) and nucleus accumbens (NAc) core and shell further revealed sex- and brain region-specific transcriptomic signatures of SST and SA. We found that genes associated with SST and SA were enriched for those related to synaptic processes, myelin sheath, and tobacco use disorder or chemdependency. Interestingly, SST-associated genes were often downregulated in male VTA but upregulated in female VTA, and strongly enriched for smoking GWAS risk variants, possibly explaining the male-specific SST. For SA, we found widespread region-specific allelic imbalance of expression (AIE), of which genes showing AIE bias toward paternal F344 alleles in NAc core were strongly enriched for SA-associated genes and for GWAS risk variants of smoking initiation, likely contributing to the parental effect of SA. Our study suggests a mechanistic link between transcriptional changes underlying the NIC SST and SA and human nicotine addiction, providing a resource for understanding the neurobiology basis of the GWAS findings on human smoking and other addictive phenotypes.

Complete Genome Sequences of Two Human Oral Microbiome Commensals, Streptococcus salivarius ATCC 25975 and S. salivarius ATCC 27945.

Streptococcus salivarius strains are significant contributors to the human oral microbiome. Some possess unique fimbriae that give them the ability to coaggregate and colonize particular oral structures. We present here the complete genomes of Streptococcus salivarius Lancefield K-/K+ strains ATCC 25975 and ATCC 27945, which can and cannot, respectively, produce fimbriae.

Genetic Characterization of the Exceptionally High Heat Resistance of the Non-toxic Surrogate Clostridium sporogenes PA 3679.

Clostridium sporogenes PA 3679 is a non-toxic endospore former that is widely used as a surrogate for Clostridium botulinum by the food processing industry to validate thermal processing strategies. PA 3679 produces spores of exceptionally high heat resistance without botulinum neurotoxins, permitting the use of PA 3679 in inoculated pack studies while ensuring the safety of food processing facilities. To identify genes associated with this heat resistance, the genomes of C. sporogenes PA 3679 isolates were compared to several other C. sporogenes strains. The most significant difference was the acquisition of a second spoVA operon, spoVA2, which is responsible for transport of dipicolinic acid into the spore core during sporulation. Interestingly, spoVA2 was also found in some C. botulinum species which phylogenetically cluster with PA 3679. Most other C. sporogenes strains examined both lack the spoVA2 locus and are phylogenetically distant within the group I Clostridium, adding to the understanding that C. sporogenes are dispersed C. botulinum strains which lack toxin genes. C. sporogenes strains are thus a very eclectic group, and few strains possess the characteristic heat resistance of PA 3679.

Enhancement of Microbial Biodesulfurization via Genetic Engineering and Adaptive Evolution.

In previous work from our laboratories a synthetic gene encoding a peptide ("Sulpeptide 1" or "S1") with a high proportion of methionine and cysteine residues had been designed to act as a sulfur sink and was inserted into the dsz (desulfurization) operon of Rhodococcus erythropolis IGTS8. In the work described here this construct (dszAS1BC) and the intact dsz operon (dszABC) cloned into vector pRESX under control of the (Rhodococcus) kstD promoter were transformed into the desulfurization-negative strain CW25 of Rhodococcus qingshengii. The resulting strains (CW25[pRESX-dszABC] and CW25[pRESX-dszAS1BC]) were subjected to adaptive selection by repeated passages at log phase (up to 100 times) in minimal medium with dibenzothiophene (DBT) as sole sulfur source. For both strains DBT metabolism peaked early in the selection process and then decreased, eventually averaging four times that of the initial transformed cells; the maximum specific activity achieved by CW25[pRESX-dszAS1BC] exceeded that of CW25[pRESX-dszABC]. Growth rates increased by 7-fold (CW25[pRESX-dszABC]) and 13-fold (CW25[pRESX-dszAS1BC]) and these increases were stable. The adaptations of CW25[pRESX-dszAS1BC] were correlated with a 3-5X increase in plasmid copy numbers from those of the initial transformed cells; whole genome sequencing indicated that during its selection processes no mutations occurred to any of the dsz, S1, or other genes and promoters involved in sulfur metabolism, stress response, or DNA methylation, and that the effect of the sulfur sink produced by S1 is likely very small compared to the cells' overall cysteine and methionine requirements. Nevertheless, a combination of genetic engineering using sulfur sinks and increasing Dsz capability with adaptive selection may be a viable strategy to increase biodesulfurization ability.

Complete Genome Sequences of Two Interactive Moderate Thermophiles, Paenibacillus napthalenovorans 32O-Y and Paenibacillus sp. 32O-W.

Microorganisms with the capability to desulfurize petroleum are in high demand with escalating restrictions currently placed on fuel purity. Thermophilic desulfurizers are particularly valuable in high-temperature industrial applications. We report the whole-genome sequences of Paenibacillus napthalenovorans 32O-Y and Paenibacillus sp. 32O-W, which can and cannot, respectively, metabolize dibenzothiophene.

Genetic Diversity of Clostridium sporogenes PA 3679 Isolates Obtained from Different Sources as Resolved by Pulsed-Field Gel Electrophoresis and High-Throughput Sequencing.

Clostridium sporogenes PA 3679 is a nonpathogenic, nontoxic model organism for proteolytic Clostridium botulinum used in the validation of conventional thermal food processes due to its ability to produce highly heat-resistant endospores. Because of its public safety importance, the uncertain taxonomic classification and genetic diversity of PA 3679 are concerns. Therefore, isolates of C. sporogenes PA 3679 were obtained from various sources and characterized using pulsed-field gel electrophoresis (PFGE) and whole-genome sequencing. The phylogenetic relatedness and genetic variability were assessed based on 16S rRNA gene sequencing and whole-genome single nucleotide polymorphism (SNP) analysis. All C. sporogenes PA 3679 isolates were categorized into two clades (clade I containing ATCC 7955 NCA3679 isolates 1961-2, 1990, and 2007 and clade II containing PA 3679 isolates NFL, UW, FDA, and Campbell and ATCC 7955 NCA3679 isolate 1961-4). The 16S maximum likelihood (ML) tree clustered both clades within proteolytic C. botulinum strains, with clade I forming a distinct cluster with other C. sporogenes non-PA 3679 strains. SNP analysis revealed that clade I isolates were more similar to the genomic reference PA 3679 (NCTC8594) genome (GenBank accession number AGAH00000000.1) than clade II isolates were. The genomic reference C. sporogenes PA 3679 (NCTC8594) genome and clade I C. sporogenes isolates were genetically distinct from those obtained from other sources (University of Wisconsin, National Food Laboratory, U.S. Food and Drug Administration, and Campbell's Soup Company). Thermal destruction studies revealed that clade I isolates were more sensitive to high temperature than clade II isolates were. Considering the widespread use of C. sporogenes PA 3679 and its genetic information in numerous studies, the accurate identification and genetic characterization of C. sporogenes PA 3679 are of critical importance.

Isolation and characterization of an interactive culture of two Paenibacillus species with moderately thermophilic desulfurization ability.

OBJECTIVE: To isolate and characterize novel thermophilic bacteria capable of biodesulfurization of petroleum. RESULTS: A culture containing two Paenibacillus spp. (denoted "32O-W" and "32O-Y") was isolated by repeated passage of a soil sample at up to 55 °C in medium containing dibenzothiophene (DBT) as sulfur source. Only 32O-Y metabolized DBT, apparently via the 4S pathway; maximum activity occurred from 40 to 45 °C, with some activity up to at least 50 °C. 32O-W enhanced DBT metabolism by 32O-Y (by 22-74 % at 40-50 °C). With sulfate as sulfur source, 32O-Y and 32O-W grew well up to 58 and 63 °C, respectively. Selection of a mixed culture of 32O-Y and 32O-W at 54 °C increased DBT metabolism 36-42 % from 40 to 45 °C. Genome sequencing identified desulfurization gene homologs in the strains consistent with their desulfurization properties. CONCLUSION: The 32O-Y/32O-W culture may be a useful starting point for development of an improved thermophilic petroleum biodesulfurization process.

Comparative Analysis of Codon Usage Bias Patterns in Microsporidian Genomes.

The sub-3 Mbp genomes from microsporidian species of the Encephalitozoon genus are the smallest known among eukaryotes and paragons of genomic reduction and compaction in parasites. However, their diminutive stature is not characteristic of all Microsporidia, whose genome sizes vary by an order of magnitude. This large variability suggests that different evolutionary forces are applied on the group as a whole. In this study, we have compared the codon usage bias (CUB) between eight taxonomically distinct microsporidian genomes: Encephalitozoon intestinalis, Encephalitozoon cuniculi, Spraguea lophii, Trachipleistophora hominis, Enterocytozoon bieneusi, Nematocida parisii, Nosema bombycis and Nosema ceranae. While the CUB was found to be weak in all eight Microsporidia, nearly all (98%) of the optimal codons in S. lophii, T. hominis, E. bieneusi, N. parisii, N. bombycis and N. ceranae are fond of A/U in third position whereas most (64.6%) optimal codons in the Encephalitozoon species E. intestinalis and E. cuniculi are biased towards G/C. Although nucleotide composition biases are likely the main factor driving the CUB in Microsporidia according to correlation analyses, directed mutational pressure also likely affects the CUB as suggested by ENc-plots, correspondence and neutrality analyses. Overall, the Encephalitozoon genomes were found to be markedly different from the other microsporidians and, despite being the first sequenced representatives of this lineage, are uncharacteristic of the group as a whole. The disparities observed cannot be attributed solely to differences in host specificity and we hypothesize that other forces are at play in the lineage leading to Encephalitozoon species.

GATA-3 dose-dependent checkpoints in early T cell commitment.

GATA-3 expression is crucial for T cell development and peaks during commitment to the T cell lineage, midway through the CD4(-)CD8(-) (double-negative [DN]) stages 1-3. We used RNA interference and conditional deletion to reduce GATA-3 protein acutely at specific points during T cell differentiation in vitro. Even moderate GATA-3 reduction killed DN1 cells, delayed progression to the DN2 stage, skewed DN2 gene regulation, and blocked appearance of the DN3 phenotype. Although a Bcl-2 transgene rescued DN1 survival and improved DN2 cell generation, it did not restore DN3 differentiation. Gene expression analyses (quantitative PCR, RNA sequencing) showed that GATA-3-deficient DN2 cells quickly upregulated genes, including Spi1 (PU.1) and Bcl11a, and downregulated genes, including Cpa3, Ets1, Zfpm1, Bcl11b, Il9r, and Il17rb with gene-specific kinetics and dose dependencies. These targets could mediate two distinct roles played by GATA-3 in lineage commitment, as revealed by removing wild-type or GATA-3-deficient early T lineage cells from environmental Notch signals. GATA-3 worked as a potent repressor of B cell potential even at low expression levels, so that only full deletion of GATA-3 enabled pro-T cells to reveal B cell potential. The ability of GATA-3 to block B cell development did not require T lineage commitment factor Bcl11b. In prethymic multipotent precursors, however, titration of GATA-3 activity using tamoxifen-inducible GATA-3 showed that GATA-3 inhibits B and myeloid developmental alternatives at different threshold doses. Furthermore, differential impacts of a GATA-3 obligate repressor construct imply that B and myeloid development are inhibited through distinct transcriptional mechanisms. Thus, the pattern of GATA-3 expression sequentially produces B lineage exclusion, T lineage progression, and myeloid-lineage exclusion for commitment.