Search

Complex traits are characterized by multiple genes and variants acting simultaneously on a phenotype. However, studying the contribution of individual pairs of genes to complex traits has been challenging since human genetics necessitates very large population sizes, while findings from model systems do not always translate to humans. Here, we combine genetics with combinatorial RNAi (coRNAi) to systematically test for pairwise additive effects (AEs) and genetic interactions (GIs) between 30 lipid genome-wide association studies (GWAS) genes. Gene-based burden tests from 240,970 exomes show that in carriers with truncating mutations in both, APOB and either PCSK9 or LPL (“human double knock-outs”) plasma lipid levels change additively. Genetics and coRNAi identify overlapping AEs for 12 additional gene pairs. Overlapping GIs are observed for TOMM40/APOE with SORT1 and NCAN. Our study identifies distinct gene pairs that modulate plasma and cellular lipid levels primarily via AEs and nominates putative drug target pairs for improved lipid-lowering combination therapies., Studying the contribution of pairs of genes to complex traits has been challenging. Here, the authors combine exome and genotype data with RNAi to screen for genetic interactions between 30 genes identified in lipid GWAS to hint at pairs whose joint modulation may improve lipid-lowering therapies.

Recent high-throughput transcription factor (TF) binding assays revealed that TF cooperativity is a widespread phenomenon. However, a global mechanistic and functional understanding of TF cooperativity is still lacking. To address this, here we introduce a statistical learning framework that provides structural insight into TF cooperativity and its functional consequences based on next generation sequencing data. We identify DNA shape as driver for cooperativity, with a particularly strong effect for Forkhead-Ets pairs. Follow-up experiments reveal a local shape preference at the Ets-DNA-Forkhead interface and decreased cooperativity upon loss of the interaction. Additionally, we discover many functional associations for cooperatively bound TFs. Examination of the link between FOXO1:ETV6 and lymphomas reveals that their joint expression levels improve patient clinical outcome stratification. Altogether, our results demonstrate that inter-family cooperative TF binding is driven by position-specific DNA readout mechanisms, which provides an additional regulatory layer for downstream biological functions., Although transcription factor (TF) cooperativity is widespread, a global mechanistic understanding of the role of TF cooperativity is still lacking. Here the authors introduce a statistical learning framework that provides structural insight into TF cooperativity and its functional consequences based on next generation sequencing data and provide mechanistic insights into TF cooperativity and its impact on protein-phenotype interactions.

Tau is a microtubule-binding protein that can receive various post-translational modifications (PTMs) including phosphorylation, methylation, acetylation, glycosylation, nitration, sumoylation and truncation. Hyperphosphorylation of tau is linked to its aggregation and the formation of neurofibrillary tangles (NFTs), which are a hallmark of Alzheimer’s disease (AD). While more than 70 phosphorylation sites have been detected previously on NFT tau, studies of oligomeric and detergent-soluble tau in human brains during the early stages of AD are lacking. Here we apply a comprehensive electrochemiluminescence ELISA assay to analyze twenty-five different PTM sites as well as tau oligomerization in control and sporadic AD brain. The samples were classified as Braak stages 0–I, II or III–IV, corresponding to the progression of microscopically detectable tau pathology throughout different brain regions. We found that soluble tau multimers are strongly increased at Braak stages III–IV in all brain regions under investigation, including the temporal cortex, which does not contain NFTs or misfolded oligomers at this stage of pathology. We additionally identified five phosphorylation sites that are specifically and consistently increased across the entorhinal cortex, hippocampus and temporal cortex in the same donors. Three of these sites correlate with tau multimerization in all three brain regions, but do not overlap with the epitopes of phospho-sensitive antibodies commonly used for the immunohistochemical detection of NFTs. Our results thus suggest that soluble multimers are characterized by a small set of specific phosphorylation events that differ from those dominating in mature NFTs. These findings shed light on early PTM changes of tau during AD pathogenesis in human brains.

Tumor recurrence is the leading cause of breast cancer-related death. Recurrences are largely driven by cancer cells that survive therapeutic intervention. This poorly understood population is referred to as minimal residual disease. Here, using mouse models that faithfully recapitulate human disease together with organoid cultures, we have demonstrated that residual cells acquire a transcriptionally distinct state from normal epithelium and primary tumors. Gene expression changes and functional characterization revealed altered lipid metabolism and elevated ROS as hallmarks of the cells that survive tumor regression. These residual cells exhibited increased oxidative DNA damage, potentiating the acquisition of somatic mutations during hormonal-induced expansion of the mammary cell population. Inhibition of either cellular fatty acid synthesis or fatty acid transport into mitochondria reduced cellular ROS levels and DNA damage, linking these features to lipid metabolism. Direct perturbation of these hallmarks in vivo, either by scavenging ROS or by halting the cyclic mammary cell population expansion, attenuated tumor recurrence. Finally, these observations were mirrored in transcriptomic and histological signatures of residual cancer cells from neoadjuvant-treated breast cancer patients. These results highlight the potential of lipid metabolism and ROS as therapeutic targets for reducing tumor recurrence in breast cancer patients.

Background While the optimal use and timing of secondary therapy after radical prostatectomy (RP) remain controversial, there are limited data on patient-reported outcomes following multimodal therapy. Objective To assess the impact of additional radiation therapy (RT) and/or androgen deprivation therapy (ADT) on urinary continence, potency, and quality of life (QoL) after RP. Design, setting, and participants Among 13150 men who underwent RP from 1992 to 2013, 905 received RP + RT, 407 RP + ADT and 688 RP + RT + ADT. Outcome measurements and statistical analyses Urinary function, sexual function, and overall QoL were evaluated annually using self-administered validated questionnaires. Propensity score–matched and bootstrap analyses were performed, and the distributions for all functional outcomes were analyzed as a function of time after RP. Results and limitations Patients who received RP + RT had a 4% higher overall incontinence rate 3 yr after surgery, and 1% higher rate for severe incontinence (>3 pads/24h) compared to matched RP-only patients. ADT further increased the overall and severe incontinence rates by 4% and 3%, respectively, compared to matched RP + RT patients. RP + RT was associated with an 18% lower rate of potency compared to RP alone, while RP + RT + ADT was associated with a further 17% reduction compared to RP + RT. Additional RT reduced QoL by 10% and additional ADT by a further 12% compared to RP only and RP + RT, respectively. The timing of RT after RP had no influence on continence, but adjuvant compared to salvage RT was associated with significantly lower potency (37% vs 45%), but higher QoL (60% vs 56%). Limitations of our study include the observational study design and potential for selection bias in the treatments received. Conclusions Secondary RT and ADT after RP have an additive negative influence on urinary function, potency, and QoL. Patients with high-risk disease should be counseled before RP on the potential net impairment of functional outcomes due to multimodal treatment. Patient summary Men with high-risk disease choosing surgery upfront should be counseled on the potential need for additional radiation and or androgen deprivation, and the potential net impairment of functional outcomes arising from multimodal treatment.

For multiple hypothesis testing in genomics and other large-scale data analyses, the independent hypothesis weighting (IHW) approach uses data-driven P-value weight assignment to improve power while controlling the false discovery rate. Hypothesis weighting improves the power of large-scale multiple testing. We describe independent hypothesis weighting (IHW), a method that assigns weights using covariates independent of the P-values under the null hypothesis but informative of each test's power or prior probability of the null hypothesis ( http://www.bioconductor.org/packages/IHW ). IHW increases power while controlling the false discovery rate and is a practical approach to discovering associations in genomics, high-throughput biology and other large data sets.

SummaryA comprehensive map of transcription start sites (TSSs) across the highly AT-rich genome of P. falciparum would aid progress toward deciphering the molecular mechanisms that underlie the timely regulation of gene expression in this malaria parasite. Using high-throughput sequencing technologies, we generated a comprehensive atlas of transcription initiation events at single-nucleotide resolution during the parasite intra-erythrocytic developmental cycle. This detailed analysis of TSS usage enabled us to define architectural features of plasmodial promoters. We demonstrate that TSS selection and strength are constrained by local nucleotide composition. Furthermore, we provide evidence for coordinate and stage-specific TSS usage from distinct sites within the same transcription unit, thereby producing transcript isoforms, a subset of which are developmentally regulated. This work offers a framework for further investigations into the interactions between genomic sequences and regulatory factors governing the complex transcriptional program of this major human pathogen.

Nuclear transport receptors (NTRs) recognize localization signals of cargos to facilitate their passage across the central channel of nuclear pore complexes (NPCs). About 30 different NTRs constitute different transport pathways in humans and bind to a multitude of different cargos. The exact cargo spectrum of the majority of NTRs, their specificity and even the extent to which active nucleocytoplasmic transport contributes to protein localization remains understudied because of the transient nature of these interactions and the wide dynamic range of cargo concentrations. To systematically map cargo–NTR relationships in situ , we used proximity ligation coupled to mass spectrometry (BioID). We systematically fused the engineered biotin ligase BirA* to 16 NTRs. We estimate that a considerable fraction of the human proteome is subject to active nuclear transport. We quantified the specificity and redundancy in NTR interactions and identified transport pathways for cargos. We extended the BioID method by the direct identification of biotinylation sites. This approach enabled us to identify interaction interfaces and to discriminate direct versus piggyback transport mechanisms. Data are available via ProteomeXchange with identifier PXD007976.

Objective Micro-RNAs (miRNAs) play a crucial role in controlling intestinal epithelial barrier function partly by modulating the expression of tight junction (TJ) proteins. We have previously shown differential messenger RNA (mRNA) expression correlated with ultrastructural abnormalities of the epithelial barrier in patients with diarrhoea-predominant IBS (IBS-D). However, the participation of miRNAs in these differential mRNA-associated findings remains to be established. Our aims were (1) to identify miRNAs differentially expressed in the small bowel mucosa of patients with IBS-D and (2) to explore putative target genes specifically involved in epithelial barrier function that are controlled by specific dysregulated IBS-D miRNAs. Design Healthy controls and patients meeting Rome III IBS-D criteria were studied. Intestinal tissue samples were analysed to identify potential candidates by: (a) miRNA-mRNA profiling; (b) miRNA-mRNA pairing analysis to assess the co-expression profile of miRNA-mRNA pairs; (c) pathway analysis and upstream regulator identification; (d) miRNA and target mRNA validation. Candidate miRNA-mRNA pairs were functionally assessed in intestinal epithelial cells. Results IBS-D samples showed distinct miRNA and mRNA profiles compared with healthy controls. TJ signalling was associated with the IBS-D transcriptional profile. Further validation of selected genes showed consistent upregulation in 75% of genes involved in epithelial barrier function. Bioinformatic analysis of putative miRNA binding sites identified hsa-miR-125b-5p and hsa-miR-16 as regulating expression of the TJ genes CGN (cingulin) and CLDN2 (claudin-2), respectively. Consistently, protein expression of CGN and CLDN2 was upregulated in IBS-D, while the respective targeting miRNAs were downregulated. In addition, bowel dysfunction, perceived stress and depression and number of mast cells correlated with the expression of hsa-miR-125b-5p and hsa-miR-16 and their respective target proteins. Conclusions Modulation of the intestinal epithelial barrier function in IBS-D involves both transcriptional and post-transcriptional mechanisms. These molecular mechanisms include miRNAs as master regulators in controlling the expression of TJ proteins and are associated with major clinical symptoms.

Although subtypes of pancreatic ductal adenocarcinoma (PDAC) have been described, this malignancy is clinically still treated as a single disease. Here we present patient-derived models representing the full spectrum of previously identified quasi-mesenchymal (QM-PDA), classical and exocrine-like PDAC subtypes, and identify two markers--HNF1A and KRT81--that enable stratification of tumors into different subtypes by using immunohistochemistry. Individuals with tumors of these subtypes showed substantial differences in overall survival, and their tumors differed in drug sensitivity, with the exocrine-like subtype being resistant to tyrosine kinase inhibitors and paclitaxel. Cytochrome P450 3A5 (CYP3A5) metabolizes these compounds in tumors of the exocrine-like subtype, and pharmacological or short hairpin RNA (shRNA)-mediated CYP3A5 inhibition sensitizes tumor cells to these drugs. Whereas hepatocyte nuclear factor 4, alpha (HNF4A) controls basal expression of CYP3A5, drug-induced CYP3A5 upregulation is mediated by the nuclear receptor NR1I2. CYP3A5 also contributes to acquired drug resistance in QM-PDA and classical PDAC, and it is highly expressed in several additional malignancies. These findings designate CYP3A5 as a predictor of therapy response and as a tumor cell-autonomous detoxification mechanism that must be overcome to prevent drug resistance.

Signal identification in large-dimensional settings is a challenging problem in biostatistics. Recently, the method of higher criticism (HC) was shown to be an effective means for determining appropriate decision thresholds. Here, we study HC from a false discovery rate (FDR) perspective. We show that the HC threshold may be viewed as an approximation to a natural class boundary (CB) in two-class discriminant analysis which in turn is expressible as FDR threshold. We demonstrate that in a rare-weak setting in the region of the phase space where signal identification is possible both thresholds are practicably indistinguishable, and thus HC thresholding is identical to using a simple local FDR cutoff. The relationship of the HC and CB thresholds and their properties are investigated both analytically and by simulations, and are further compared by application to four cancer gene expression data sets., Comment: 19 pages, 3 figures, 2 tables