Association of clonal hematopoiesis and mosaic chromosomal alterations with solid malignancy incidence and mortality.

BACKGROUND: Understanding the impact of clonal hematopoiesis of indeterminate potential (CHIP) and mosaic chromosomal alterations (mCAs) on solid tumor risk and mortality can shed light on novel cancer pathways. METHODS: The authors analyzed whole genome sequencing data from the Trans-Omics for Precision Medicine Women's Health Initiative study (n = 10,866). They investigated the presence of CHIP and mCA and their association with the development and mortality of breast, lung, and colorectal cancers. RESULTS: CHIP was associated with higher risk of breast (hazard ratio [HR], 1.30; 95% confidence interval [CI], 1.03-1.64; p = .02) but not colorectal (p = .77) or lung cancer (p = .32). CHIP carriers who developed colorectal cancer also had a greater risk for advanced-stage (p = .01), but this was not seen in breast or lung cancer. CHIP was associated with increased colorectal cancer mortality both with (HR, 3.99; 95% CI, 2.41-6.62; p < .001) and without adjustment (HR, 2.50; 95% CI, 1.32-4.72; p = .004) for advanced-stage and a borderline higher breast cancer mortality (HR, 1.53; 95% CI, 0.98-2.41; p = .06). Conversely, mCA (cell fraction [CF] >3%) did not correlate with cancer risk. With higher CFs (mCA >5%), autosomal mCA was associated with increased breast cancer risk (HR, 1.39; 95% CI, 1.06-1.83; p = .01). There was no association of mCA (>3%) with breast, colorectal, or lung mortality except higher colon cancer mortality (HR, 2.19; 95% CI, 1.11-4.3; p = .02) with mCA >5%. CONCLUSIONS: CHIP and mCA (CF >5%) were associated with higher breast cancer risk and colorectal cancer mortality individually. These data could inform on novel pathways that impact cancer risk and lead to better risk stratification.

Investigating SH-SY5Y Neuroblastoma Cell Surfaceome as a Model for Neuronal-Targeted Novel Therapeutic Modalities.

The SH-SY5Y neuroblastoma cells are a widely used in vitro model approximating neurons for testing the target engagement of therapeutics designed for neurodegenerative diseases and pain disorders. However, their potential as a model for receptor-mediated delivery and uptake of novel modalities, such as antibody-drug conjugates, remains understudied. Investigation of the SH-SY5Y cell surfaceome will aid in greater in vitro to in vivo correlation of delivery and uptake, thereby accelerating drug discovery. So far, the majority of studies have focused on total cell proteomics from undifferentiated and differentiated SH-SY5Y cells. While some studies have investigated the expression of specific proteins in neuroblastoma tissue, a global approach for comparison of neuroblastoma cell surfaceome to the brain and dorsal root ganglion (DRG) neurons remains uninvestigated. Furthermore, an isoform-specific evaluation of cell surface proteins expressed on neuroblastoma cells remains unexplored. In this study, we define a bioinformatic workflow for the identification of high-confidence surface proteins expressed on brain and DRG neurons using tissue proteomic and transcriptomic data. We then delineate the SH-SY5Y cell surfaceome by surface proteomics and show that it significantly overlaps with the human brain and DRG neuronal surface proteome. We find that, for 32% of common surface proteins, SH-SY5Y-specific major isoforms are alternatively spliced, maintaining their protein-coding ability, and are predicted to localize to the cell surface. Validation of these isoforms using surface proteomics confirms a SH-SY5Y-specific alternative NRCAM (neuron-glia related cell adhesion molecule) isoform, which is absent in typical brain neurons, but present in neuroblastomas, making it a receptor of interest for neuroblastoma-specific therapeutics.

Population-based meta-analysis and gene-set enrichment identifies FXR/RXR pathway as common to fatty liver disease and serum lipids.

Nonalcoholic fatty liver disease (NAFLD) is prevalent worldwide. NAFLD is associated with elevated serum triglycerides (TG), low-density lipoprotein cholesterol (LDL), and reduced high-density lipoprotein cholesterol (HDL). Both NAFLD and blood lipid levels are genetically influenced and may share a common genetic etiology. We used genome-wide association studies (GWAS)-ranked genes and gene-set enrichment analysis to identify pathways that affect serum lipids and NAFLD. We identified credible genes in these pathways and characterized missense variants in these for effects on serum traits. We used MAGENTA to identify 58 enriched pathways from publicly available TG, LDL, and HDL GWAS (n = 99,000). Three of these pathways were also enriched for associations with European-ancestry NAFLD GWAS (n = 7176). One pathway, farnesoid X receptor (FXR)/retinoid X receptor (RXR) activation, was replicated for association in an African-ancestry NAFLD GWAS (n = 3214) and plays a role in serum lipids and NAFLD. Credible genes (proteins) in FXR/RXR activation include those associated with cholesterol/bile/bilirubin transport/absorption (ABCC2 (MRP2) [ATP binding cassette subfamily C member (multidrug resistance-associated protein 2)], ABCG5, ABCG8 [ATP-binding cassette (ABC) transporters G5 and G8], APOB (APOB) [apolipoprotein B], FABP6 (ILBP) [fatty acid binding protein 6 (ileal lipid-binding protein)], MTTP (MTP) [microsomal triglyceride transfer protein], SLC4A2 (AE2) [solute carrier family 4 member 2 (anion exchange protein 2)]), nuclear hormone-mediated control of metabolism (NR0B2 (SHP) [nuclear receptor subfamily 0 group B member 2 (small heterodimer partner)], NR1H4 (FXR) [nuclear receptor subfamily 1 group H member 4 (FXR)], PPARA (PPAR) [peroxisome proliferator activated receptor alpha], FOXO1 (FOXO1A) [forkhead box O1]), or other pathways (FETUB (FETUB) [fetuin B]). Missense variants in ABCC2 (MRP2), ABCG5 (ABCG5), ABCG8 (ABCG8), APOB (APOB), MTTP (MTP), NR0B2 (SHP), NR1H4 (FXR), and PPARA (PPAR) that associate with serum LDL levels also associate with serum liver function tests in UK Biobank. Conclusion: Genetic variants in NR1H4 (FXR) that protect against liver steatosis increase serum LDL cholesterol while variants in other members of the family have congruent effects on these traits. Human genetic pathway enrichment analysis can help guide therapeutic development by identifying effective targets for NAFLD/serum lipid manipulation while minimizing side effects. In addition, missense variants could be used in companion diagnostics to determine their influence on drug effectiveness.

Morphological cell profiling of SARS-CoV-2 infection identifies drug repurposing candidates for COVID-19.

The global spread of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), and the associated disease COVID-19, requires therapeutic interventions that can be rapidly identified and translated to clinical care. Traditional drug discovery methods have a >90% failure rate and can take 10 to 15 y from target identification to clinical use. In contrast, drug repurposing can significantly accelerate translation. We developed a quantitative high-throughput screen to identify efficacious agents against SARS-CoV-2. From a library of 1,425 US Food and Drug Administration (FDA)-approved compounds and clinical candidates, we identified 17 hits that inhibited SARS-CoV-2 infection and analyzed their antiviral activity across multiple cell lines, including lymph node carcinoma of the prostate (LNCaP) cells and a physiologically relevant model of alveolar epithelial type 2 cells (iAEC2s). Additionally, we found that inhibitors of the Ras/Raf/MEK/ERK signaling pathway exacerbate SARS-CoV-2 infection in vitro. Notably, we discovered that lactoferrin, a glycoprotein found in secretory fluids including mammalian milk, inhibits SARS-CoV-2 infection in the nanomolar range in all cell models with multiple modes of action, including blockage of virus attachment to cellular heparan sulfate and enhancement of interferon responses. Given its safety profile, lactoferrin is a readily translatable therapeutic option for the management of COVID-19.

BET Bromodomain Inhibition Blocks an AR-Repressed, E2F1-Activated Treatment-Emergent Neuroendocrine Prostate Cancer Lineage Plasticity Program.

PURPOSE: Lineage plasticity in prostate cancer-most commonly exemplified by loss of androgen receptor (AR) signaling and a switch from a luminal to alternate differentiation program-is now recognized as a treatment resistance mechanism. Lineage plasticity is a spectrum, but neuroendocrine prostate cancer (NEPC) is the most virulent example. Currently, there are limited treatments for NEPC. Moreover, the incidence of treatment-emergent NEPC (t-NEPC) is increasing in the era of novel AR inhibitors. In contradistinction to de novo NEPC, t-NEPC tumors often express the AR, but AR's functional role in t-NEPC is unknown. Furthermore, targetable factors that promote t-NEPC lineage plasticity are also unclear. EXPERIMENTAL DESIGN: Using an integrative systems biology approach, we investigated enzalutamide-resistant t-NEPC cell lines and their parental, enzalutamide-sensitive adenocarcinoma cell lines. The AR is still expressed in these t-NEPC cells, enabling us to determine the role of the AR and other key factors in regulating t-NEPC lineage plasticity. RESULTS: AR inhibition accentuates lineage plasticity in t-NEPC cells-an effect not observed in parental, enzalutamide-sensitive adenocarcinoma cells. Induction of an AR-repressed, lineage plasticity program is dependent on activation of the transcription factor E2F1 in concert with the BET bromodomain chromatin reader BRD4. BET inhibition (BETi) blocks this E2F1/BRD4-regulated program and decreases growth of t-NEPC tumor models and a subset of t-NEPC patient tumors with high activity of this program in a BETi clinical trial. CONCLUSIONS: E2F1 and BRD4 are critical for activating an AR-repressed, t-NEPC lineage plasticity program. BETi is a promising approach to block this program.

Standing Variations Modeling Captures Inter-Individual Heterogeneity in a Deterministic Model of Prostate Cancer Response to Combination Therapy.

Sipuleucel-T (Provenge) is the first live cell vaccine approved for advanced, hormonally refractive prostate cancer. However, survival benefit is modest and the optimal combination or schedule of sipuleucel-T with androgen depletion remains unknown. We employ a nonlinear dynamical systems approach to modeling the response of hormonally refractive prostate cancer to sipuleucel-T. Our mechanistic model incorporates the immune response to the cancer elicited by vaccination, and the effect of androgen depletion therapy. Because only a fraction of patients benefit from sipuleucel-T treatment, inter-individual heterogeneity is clearly crucial. Therefore, we introduce our novel approach, Standing Variations Modeling, which exploits inestimability of model parameters to capture heterogeneity in a deterministic model. We use data from mouse xenograft experiments to infer distributions on parameters critical to tumor growth and to the resultant immune response. Sampling model parameters from these distributions allows us to represent heterogeneity, both at the level of the tumor cells and the individual (mouse) being treated. Our model simulations explain the limited success of sipuleucel-T observed in practice, and predict an optimal combination regime that maximizes predicted efficacy. This approach will generalize to a range of emerging cancer immunotherapies.

A Noncoding Variant Near PPP1R3B Promotes Liver Glycogen Storage and MetS, but Protects Against Myocardial Infarction.

CONTEXT: Glycogen storage diseases are rare. Increased glycogen in the liver results in increased attenuation. OBJECTIVE: Investigate the association and function of a noncoding region associated with liver attenuation but not histologic nonalcoholic fatty liver disease. DESIGN: Genetics of Obesity-associated Liver Disease Consortium. SETTING: Population-based. MAIN OUTCOME: Computed tomography measured liver attenuation. RESULTS: Carriers of rs4841132-A (frequency 2%-19%) do not show increased hepatic steatosis; they have increased liver attenuation indicative of increased glycogen deposition. rs4841132 falls in a noncoding RNA LOC157273 ~190 kb upstream of PPP1R3B. We demonstrate that rs4841132-A increases PPP1R3B through a cis genetic effect. Using CRISPR/Cas9 we engineered a 105-bp deletion including rs4841132-A in human hepatocarcinoma cells that increases PPP1R3B, decreases LOC157273, and increases glycogen perfectly mirroring the human disease. Overexpression of PPP1R3B or knockdown of LOC157273 increased glycogen but did not result in decreased LOC157273 or increased PPP1R3B, respectively, suggesting that the effects may not all occur via affecting RNA levels. Based on electronic health record (EHR) data, rs4841132-A associates with all components of the metabolic syndrome (MetS). However, rs4841132-A associated with decreased low-density lipoprotein (LDL) cholesterol and risk for myocardial infarction (MI). A metabolic signature for rs4841132-A includes increased glycine, lactate, triglycerides, and decreased acetoacetate and beta-hydroxybutyrate. CONCLUSIONS: These results show that rs4841132-A promotes a hepatic glycogen storage disease by increasing PPP1R3B and decreasing LOC157273. rs4841132-A promotes glycogen accumulation and development of MetS but lowers LDL cholesterol and risk for MI. These results suggest that elevated hepatic glycogen is one cause of MetS that does not invariably promote MI.

Morphological Cell Profiling of SARS-CoV-2 Infection Identifies Drug Repurposing Candidates for COVID-19.

The global spread of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), and the associated disease COVID-19, requires therapeutic interventions that can be rapidly identified and translated to clinical care. Traditional drug discovery methods have a >90% failure rate and can take 10-15 years from target identification to clinical use. In contrast, drug repurposing can significantly accelerate translation. We developed a quantitative high-throughput screen to identify efficacious agents against SARS-CoV-2. From a library of 1,425 FDA-approved compounds and clinical candidates, we identified 17 dose-responsive compounds with in vitro antiviral efficacy in human liver Huh7 cells and confirmed antiviral efficacy in human colon carcinoma Caco-2, human prostate adenocarcinoma LNCaP, and in a physiologic relevant model of alveolar epithelial type 2 cells (iAEC2s). Additionally, we found that inhibitors of the Ras/Raf/MEK/ERK signaling pathway exacerbate SARS-CoV-2 infection in vitro. Notably, we discovered that lactoferrin, a glycoprotein classically found in secretory fluids, including mammalian milk, inhibits SARS-CoV-2 infection in the nanomolar range in all cell models with multiple modes of action, including blockage of virus attachment to cellular heparan sulfate and enhancement of interferon responses. Given its safety profile, lactoferrin is a readily translatable therapeutic option for the management of COVID-19.

Mitigating temozolomide resistance in glioblastoma via DNA damage-repair inhibition.

Glioblastomas are among the most lethal cancers, with a 5 year survival rate below 25%. Temozolomide is typically used in glioblastoma treatment; however, the enzymes alkylpurine-DNA-N-glycosylase (APNG) and methylguanine-DNA-methyltransferase (MGMT) efficiently mediate the repair of DNA damage caused by temozolomide, reducing treatment efficacy. Consequently, APNG and MGMT inhibition has been proposed as a way of overcoming chemotherapy resistance. Here, we develop a mechanistic mathematical model that explicitly incorporates the effects of chemotherapy on tumour cells, including the processes of DNA damage induction, cell arrest and DNA repair. Our model is carefully parametrized and validated, and then used to virtually recreate the response of heteroclonal glioblastomas to dual treatment with temozolomide and inhibitors of APNG/MGMT. Using our mechanistic model, we identify four combination treatment strategies optimized by tumour cell phenotype, and isolate the strategy most likely to succeed in a pre-clinical and clinical setting. If confirmed in clinical trials, these strategies have the potential to offset chemotherapy resistance in patients with glioblastoma and improve overall survival.

The plasma metabolome of women in early pregnancy differs from that of non-pregnant women.

BACKGROUND: In comparison to the non-pregnant state, the first trimester of pregnancy is characterized by systemic adaptation of the mother. The extent to which these adaptive processes are reflected in the maternal blood metabolome is not well characterized. OBJECTIVE: To determine the differences between the plasma metabolome of non-pregnant and pregnant women before 16 weeks gestation. STUDY DESIGN: This study included plasma samples from 21 non-pregnant women and 50 women with a normal pregnancy (8-16 weeks of gestation). Combined measurements by ultrahigh performance liquid chromatography/tandem mass spectrometry and by gas chromatography/mass spectrometry generated molecular abundance measurements for each sample. Molecular species detected in at least 10 samples were included in the analysis. Differential abundance was inferred based on false discovery adjusted p-values (FDR) from Mann-Whitney-Wilcoxon U tests <0.1 and a minimum median abundance ratio (fold change) of 1.5. Alternatively, metabolic data were quantile normalized to remove sample-to-sample differences in the overall metabolite abundance (adjusted analysis). RESULTS: Overall, 637 small molecules met the inclusion criteria and were tested for association with pregnancy; 44% (281/637) of small molecules had significantly different abundance, of which 81% (229/281) were less abundant in pregnant than in non-pregnant women. Eight percent (14/169) of the metabolites that remained significant in the adjusted analysis also changed as a function of gestational age. A pathway analysis revealed enrichment in steroid metabolites related to sex hormones, caffeine metabolites, lysolipids, dipeptides, and polypeptide bradykinin derivatives (all, FDR < 0.1). CONCLUSIONS: This high-throughput mass spectrometry study identified: 1) differences between pregnant vs. non-pregnant women in the abundance of 44% of the profiled plasma metabolites, including known and novel molecules and pathways; and 2) specific metabolites that changed with gestational age.

The focal adhesion scaffold protein Hic-5 regulates vimentin organization in fibroblasts.

Focal adhesion (FA)-stimulated reorganization of the F-actin cytoskeleton regulates cellular size, shape, and mechanical properties. However, FA cross-talk with the intermediate filament cytoskeleton is poorly understood. Genetic ablation of the FA-associated scaffold protein Hic-5 in mouse cancer-associated fibroblasts (CAFs) promoted a dramatic collapse of the vimentin network, which was rescued following EGFP-Hic-5 expression. Vimentin collapse correlated with a loss of detergent-soluble vimentin filament precursors and decreased vimentin S72/S82 phosphorylation. Additionally, fluorescence recovery after photobleaching analysis indicated impaired vimentin dynamics. Microtubule (MT)-associated EB1 tracking and Western blotting of MT posttranslational modifications indicated no change in MT dynamics that could explain the vimentin collapse. However, pharmacological inhibition of the RhoGTPase Cdc42 in Hic-5 knockout CAFs rescued the vimentin collapse, while pan-formin inhibition with SMIFH2 promoted vimentin collapse in Hic-5 heterozygous CAFs. Our results reveal novel regulation of vimentin organization/dynamics by the FA scaffold protein Hic-5 via modulation of RhoGTPases and downstream formin activity.

Identifying Interaction Clusters for MiRNA and MRNA Pairs in TCGA Network.

Existing methods often fail to recognize the conversions for the biological roles of the pairs of genes and microRNAs (miRNAs) between the tumor and normal samples. We have developed a novel cluster scoring method to identify messenger RNA (mRNA) and miRNA interaction pairs and clusters while considering tumor and normal samples jointly. Our method has identified 54 significant clusters for 15 cancer types selected from The Cancer Genome Atlas project. We also determined the shared clusters across tumor types and/or subtypes. In addition, we compared gene and miRNA overlap between lists identified in our liver hepatocellular carcinoma (LIHC) study and regulatory relationships reported from human and rat nonalcoholic fatty liver disease studies (NAFLD). Finally, we analyzed biological functions for the single significant cluster in LIHC and uncovered a significantly enriched pathway (phospholipase D signaling pathway) with six genes represented in the cluster, symbols: DGKQ, LPAR2, PDGFRB, PIK3R3, PTGFR and RAPGEF3.

17-Beta Hydroxysteroid Dehydrogenase 13 Is a Hepatic Retinol Dehydrogenase Associated With Histological Features of Nonalcoholic Fatty Liver Disease.

Nonalcoholic fatty liver disease (NAFLD) is a common cause of chronic liver disease. A single-nucleotide polymorphism (SNP), rs6834314, was associated with serum liver enzymes in the general population, presumably reflecting liver fat or injury. We studied rs6834314 and its nearest gene, 17-beta hydroxysteroid dehydrogenase 13 (HSD17B13), to identify associations with histological features of NAFLD and to characterize the functional role of HSD17B13 in NAFLD pathogenesis. The minor allele of rs6834314 was significantly associated with increased steatosis but decreased inflammation, ballooning, Mallory-Denk bodies, and liver enzyme levels in 768 adult Caucasians with biopsy-proven NAFLD and with cirrhosis in the general population. We found two plausible causative variants in the HSD17B13 gene. rs72613567, a splice-site SNP in high linkage with rs6834314 (r2 = 0.94) generates splice variants and shows a similar pattern of association with NAFLD histology. Its minor allele generates simultaneous expression of exon 6-skipping and G-nucleotide insertion variants. Another SNP, rs62305723 (encoding a P260S mutation), is significantly associated with decreased ballooning and inflammation. Hepatic expression of HSD17B13 is 5.9-fold higher (P = 0.003) in patients with NAFLD. HSD17B13 is targeted to lipid droplets, requiring the conserved amino acid 22-28 sequence and amino acid 71-106 region. The protein has retinol dehydrogenase (RDH) activity, with enzymatic activity dependent on lipid droplet targeting and cofactor binding site. The exon 6 deletion, G insertion, and naturally occurring P260S mutation all confer loss of enzymatic activity. Conclusion: We demonstrate the association of variants in HSD17B13 with specific features of NAFLD histology and identify the enzyme as a lipid droplet-associated RDH; our data suggest that HSD17B13 plays a role in NAFLD through its enzymatic activity.

AmpliSeq transcriptome analysis of human alveolar and monocyte-derived macrophages over time in response to Mycobacterium tuberculosis infection.

Human alveolar macrophages (HAM) are primary bacterial niche and immune response cells during Mycobacterium tuberculosis (M.tb) infection, and human blood monocyte-derived macrophages (MDM) are a model for investigating M.tb-macrophage interactions. Here, we use a targeted RNA-Seq method to measure transcriptome-wide changes in RNA expression patterns of freshly obtained HAM (used within 6 h) and 6 day cultured MDM upon M.tb infection over time (2, 24 and 72 h), in both uninfected and infected cells from three donors each. The Ion AmpliSeq™ Transcriptome Human Gene Expression Kit (AmpliSeq) uses primers targeting 18,574 mRNAs and 2,228 non-coding RNAs (ncRNAs) for a total of 20,802 transcripts. AmpliSeqTM yields highly precise and reproducible gene expression profiles (R2 >0.99). Taking advantage of AmpliSeq's reproducibility, we establish well-defined quantitative RNA expression patterns of HAM versus MDM, including significant M.tb-inducible genes, in networks and pathways that differ in part between MDM and HAM. A similar number of expressed genes are detected at all time-points between uninfected MDM and HAM, in common pathways including inflammatory and immune functions, but canonical pathway differences also exist. In particular, at 2 h, multiple genes relevant to the immune response are preferentially expressed in either uninfected HAM or MDM, while the HAM RNA profiles approximate MDM profiles over time in culture, highlighting the unique RNA expression profile of freshly obtained HAM. MDM demonstrate a greater transcriptional response than HAM upon M.tb infection, with 2 to >10 times more genes up- or down-regulated. The results identify key genes involved in cellular responses to M.tb in two different human macrophage types. Follow-up bioinformatics analysis indicates that approximately 30% of response genes have expression quantitative trait loci (eQTLs in GTEx), common DNA variants that can influence host gene expression susceptibility or resistance to M.tb, illustrated with the TREM1 gene cluster and IL-10.

Gene-Hormone Therapy Interaction and Fracture Risk in Postmenopausal Women.

Context: Evidence supports a protective effect of menopausal hormone therapy (HT) on bone. However, whether genetic susceptibility modifies the association of HT and fracture risk is not sufficiently explored. Objective: The objective was to test an interaction between genetic susceptibility and HT on fracture risk. Design: We constructed two weighted genetic risk scores (GRSs) based on 16 fracture-associated variants (Fx-GRSs) and 50 bone mineral density variants (BMD-GRSs). We used Cox regression to estimate the main effects of GRSs and their interactions with HT on fracture risk. We estimated the relative excess risk due to interaction (RERI) as a measure of additive interaction. We also used the case-only approach to test for a multiplicative interaction. Setting: Forty US clinical centers. Participants: A total of 9922 genotyped white postmenopausal women (age, 50 to 79) from the Women's Health Initiative HT randomized trials. Main Outcome Measures: Adjudicated fracture incidence. Results: Both GRSs were associated with fracture risk per 1-unit increment in GRS (hazard ratio, 1.04 [95% confidence interval, 1.02 to 1.06] for Fx-GRS and hazard ratio, 1.03 [95% confidence interval,1.02-1.04] for BMD-GRS). We found no evidence for multiplicative interaction for either of the GRS. However, we observed a substantial additive interaction, where the highest quartile of both GRSs and randomization to placebo have excess fracture risk: Fx-GRS P for RERI = 0.047, BMD-GRS P for RERI = 0.046. Conclusions: These results suggest that HT reduces fracture risk in postmenopausal women, especially in those at highest genetic risk of fracture and low BMD.

RNA sequencing of transcriptomes in human brain regions: protein-coding and non-coding RNAs, isoforms and alleles.

BACKGROUND: We used RNA sequencing to analyze transcript profiles of ten autopsy brain regions from ten subjects. RNA sequencing techniques were designed to detect both coding and non-coding RNA, splice isoform composition, and allelic expression. Brain regions were selected from five subjects with a documented history of smoking and five non-smokers. Paired-end RNA sequencing was performed on SOLiD instruments to a depth of >40 million reads, using linearly amplified, ribosomally depleted RNA. Sequencing libraries were prepared with both poly-dT and random hexamer primers to detect all RNA classes, including long non-coding (lncRNA), intronic and intergenic transcripts, and transcripts lacking poly-A tails, providing additional data not previously available. The study was designed to generate a database of the complete transcriptomes in brain region for gene network analyses and discovery of regulatory variants. RESULTS: Of 20,318 protein coding and 18,080 lncRNA genes annotated from GENCODE and lncipedia, 12 thousand protein coding and 2 thousand lncRNA transcripts were detectable at a conservative threshold. Of the aligned reads, 52 % were exonic, 34 % intronic and 14 % intergenic. A majority of protein coding genes (65 %) was expressed in all regions, whereas ncRNAs displayed a more restricted distribution. Profiles of RNA isoforms varied across brain regions and subjects at multiple gene loci, with neurexin 3 (NRXN3) a prominent example. Allelic RNA ratios deviating from unity were identified in > 400 genes, detectable in both protein-coding and non-coding genes, indicating the presence of cis-acting regulatory variants. Mathematical modeling was used to identify RNAs stably expressed in all brain regions (serving as potential markers for normalizing expression levels), linked to basic cellular functions. An initial analysis of differential expression analysis between smokers and nonsmokers implicated a number of genes, several previously associated with nicotine exposure. CONCLUSIONS: RNA sequencing identifies distinct and consistent differences in gene expression between brain regions, with non-coding RNA displaying greater diversity between brain regions than mRNAs. Numerous RNAs exhibit robust allele selective expression, proving a means for discovery of cis-acting regulatory factors with potential clinical relevance.

Association between Prediagnostic Allergy-Related Serum Cytokines and Glioma.

Allergy is inversely related to glioma risk. To determine whether prediagnostic allergy-related serum proteins are associated with glioma, we conducted a nested case-control study of seven cytokines (IL4, IL13, IL5, IL6, IL10, IFNG, TGFB2), two soluble cytokine receptors (sIL4RA, sIL13RA2) and three allergy-related transcription factors (FOXP3, STAT3, STAT6) using serum specimens from the Janus Serum Bank Cohort in Oslo, Norway. Blood donors subsequently diagnosed with glioma (n = 487) were matched to controls (n = 487) on age and date of blood draw and sex. We first estimated individual effects of the 12 serum proteins and then interactions between IL4 and IL13 and their receptors using conditional logistic regression. We next tested equality of case-control inter-correlations among the 12 serum proteins. We found that TGFB2 is inversely related to glioblastoma (Odds Ratio (OR) = 0.87, 95% Confidence Interval (CI)) = 0.76, 0.98). In addition, ≤ 5 years before diagnosis, we observed associations between IL4 (OR = 0.82, 95% CI = 0.66, 1.01), sIL4RA (OR = 0.80, 95% CI = 0.65, 1.00), their interaction (OR = 1.06, 95% CI = 1.01, 1.12) and glioblastoma. This interaction was apparent > 20 years before diagnosis (IL4-sIL4RA OR = 1.20, 95% CI = 1.05, 1.37). Findings for glioma were similar. Case correlations were different from control correlations stratified on time before diagnosis. Five years or less before diagnosis, correlations among case serum proteins were weaker than were those among controls. Our findings suggest that IL4 and sIL4RA reduce glioma risk long before diagnosis and early gliomagenesis affects circulating immune function proteins.

Meta-analysis of genome-wide association studies identifies two loci associated with circulating osteoprotegerin levels.

Osteoprotegerin (OPG) is involved in bone homeostasis and tumor cell survival. Circulating OPG levels are also important biomarkers of various clinical traits, such as cancers and atherosclerosis. OPG levels were measured in serum or in plasma. In a meta-analysis of genome-wide association studies in up to 10 336 individuals from European and Asian origin, we discovered that variants >100 kb upstream of the TNFRSF11B gene encoding OPG and another new locus on chromosome 17q11.2 were significantly associated with OPG variation. We also identified a suggestive locus on chromosome 14q21.2 associated with the trait. Moreover, we estimated that over half of the heritability of OPG levels could be explained by all variants examined in our study. Our findings provide further insight into the genetic regulation of circulating OPG levels.

Whole transcriptome RNA-Seq allelic expression in human brain.

BACKGROUND: Measuring allelic RNA expression ratios is a powerful approach for detecting cis-acting regulatory variants, RNA editing, loss of heterozygosity in cancer, copy number variation, and allele-specific epigenetic gene silencing. Whole transcriptome RNA sequencing (RNA-Seq) has emerged as a genome-wide tool for identifying allelic expression imbalance (AEI), but numerous factors bias allelic RNA ratio measurements. Here, we compare RNA-Seq allelic ratios measured in nine different human brain regions with a highly sensitive and accurate SNaPshot measure of allelic RNA ratios, identifying factors affecting reliable allelic ratio measurement. Accounting for these factors, we subsequently surveyed the variability of RNA editing across brain regions and across individuals. RESULTS: We find that RNA-Seq allelic ratios from standard alignment methods correlate poorly with SNaPshot, but applying alternative alignment strategies and correcting for observed biases significantly improves correlations. Deploying these methods on a transcriptome-wide basis in nine brain regions from a single individual, we identified genes with AEI across all regions (SLC1A3, NHP2L1) and many others with region-specific AEI. In dorsolateral prefrontal cortex (DLPFC) tissues from 14 individuals, we found evidence for frequent regulatory variants affecting RNA expression in tens to hundreds of genes, depending on stringency for assigning AEI. Further, we find that the extent and variability of RNA editing is similar across brain regions and across individuals. CONCLUSIONS: These results identify critical factors affecting allelic ratios measured by RNA-Seq and provide a foundation for using this technology to screen allelic RNA expression on a transcriptome-wide basis. Using this technology as a screening tool reveals tens to hundreds of genes harboring frequent functional variants affecting RNA expression in the human brain. With respect to RNA editing, the similarities within and between individuals leads us to conclude that this post-transcriptional process is under heavy regulatory influence to maintain an optimal degree of editing for normal biological function.

Identification of eukaryotic and prokaryotic methylthiotransferase for biosynthesis of 2-methylthio-N6-threonylcarbamoyladenosine in tRNA.

Bacterial and eukaryotic transfer RNAs have been shown to contain hypermodified adenosine, 2-methylthio-N(6)-threonylcarbamoyladenosine, at position 37 (A(37)) adjacent to the 3'-end of the anticodon, which is essential for efficient and highly accurate protein translation by the ribosome. Using a combination of bioinformatic sequence analysis and in vivo assay coupled to HPLC/MS technique, we have identified, from distinct sequence signatures, two methylthiotransferase (MTTase) subfamilies, designated as MtaB in bacterial cells and e-MtaB in eukaryotic and archaeal cells. Both subfamilies are responsible for the transformation of N(6)-threonylcarbamoyladenosine into 2-methylthio-N(6)-threonylcarbamoyladenosine. Recently, a variant within the human CDKAL1 gene belonging to the e-MtaB subfamily was shown to predispose for type 2 diabetes. CDKAL1 is thus the first eukaryotic MTTase identified so far. Using purified preparations of Bacillus subtilis MtaB (YqeV), a CDKAL1 bacterial homolog, we demonstrate that YqeV/CDKAL1 enzymes, as the previously studied MTTases MiaB and RimO, contain two [4Fe-4S] clusters. This work lays the foundation for elucidating the function of CDKAL1.