SpeCond: a method to detect condition-specific gene expression

Transcriptomic studies routinely measure expression levels across numerous conditions. These datasets allow identification of genes that are specifically expressed in a small number of conditions. However, there are currently no statistically robust methods for identifying such genes. Here we present SpeCond, a method to detect condition-specific genes that outperforms alternative approaches. We apply the method to a dataset of 32 human tissues to determine 2,673 specifically expressed genes. An implementation of SpeCond is freely available as a Bioconductor package at http://www.bioconductor.org/packages/release/bioc/html/SpeCond.html

Improving Precancerous Case Characterization via Transformer-based Ensemble Learning

The application of natural language processing (NLP) to cancer pathology reports has been focused on detecting cancer cases, largely ignoring precancerous cases. Improving the characterization of precancerous adenomas assists in developing diagnostic tests for early cancer detection and prevention, especially for colorectal cancer (CRC). Here we developed transformer-based deep neural network NLP models to perform the CRC phenotyping, with the goal of extracting precancerous lesion attributes and distinguishing cancer and precancerous cases. We achieved 0.914 macro-F1 scores for classifying patients into negative, non-advanced adenoma, advanced adenoma and CRC. We further improved the performance to 0.923 using an ensemble of classifiers for cancer status classification and lesion size named entity recognition (NER). Our results demonstrated the potential of using NLP to leverage real-world health record data to facilitate the development of diagnostic tests for early cancer prevention

Genome-wide allele- and strand-specific expression profiling

Recent reports have shown that most of the genome is transcribed and that transcription frequently occurs concurrently on both DNA strands. In diploid genomes, the expression level of each allele conditions the degree to which sequence polymorphisms affect the phenotype. It is thus essential to quantify expression in an allele- and strand-specific manner. Using a custom-designed tiling array and a new computational approach, we piloted measuring allele- and strand-specific expression in yeast. Confident quantitative estimates of allele-specific expression were obtained for about half of the coding and non-coding transcripts of a heterozygous yeast strain, of which 371 transcripts (13%) showed significant allelic differential expression greater than 1.5-fold. The data revealed complex allelic differential expression on opposite strands. Furthermore, combining allele-specific expression with linkage mapping enabled identifying allelic variants that act in cis and in trans to regulate allelic expression in the heterozygous strain. Our results provide the first high-resolution analysis of differential expression on all four strands of an eukaryotic genome

Dietary strategies of Pleistocene Pongo sp. and Homo erectus on Java (Indonesia)

During the Early to Middle Pleistocene, Java was inhabited by hominid taxa of great diversity. However, their seasonal dietary strategies have never been explored. We undertook geochemical analyses of orangutan (Pongo sp.), Homo erectus and other mammalian Pleistocene teeth from Sangiran. We reconstructed past dietary strategies at subweekly resolution and inferred seasonal ecological patterns. Histologically controlled spatially resolved elemental analyses by laser-based plasma mass spectrometry confirmed the preservation of authentic biogenic signals despite the effect of spatially restricted diagenetic overprint. The Sr/Ca record of faunal remains is in line with expected trophic positions, contextualizing fossil hominid diet. Pongo sp. displays marked seasonal cycles with ~3 month-long strongly elevated Sr/Ca peaks, reflecting contrasting plant food consumption presumably during the monsoon season, while lower Sr/Ca ratios suggest different food availability during the dry season. In contrast, omnivorous H. erectus shows low and less accentuated intra-annual Sr/Ca variability compared to Pongo sp., with δ13C data of one individual indicating a dietary shift from C4 to a mix of C3 and C4 plants. Our data suggest that H. erectus on Java was maximizing the resources available in more open mosaic habitats and was less dependent on variations in seasonal resource availability. While still influenced by seasonal food availability, we infer that H. erectus was affected to a lesser degree than Pongo sp., which inhabited monsoonal rain forests on Java. We suggest that H. erectus maintained a greater degree of nutritional independence by exploiting the regional diversity of food resources across the seasons

Cognitive performance at first episode of psychosis and the relationship with future treatment resistance: Evidence from an international prospective cohort study

Background: Antipsychotic treatment resistance affects up to a third of individuals with schizophrenia, with recent research finding systematic biological differences between antipsychotic resistant and responsive patients. Our aim was to determine whether cognitive impairment at first episode significantly differs between future antipsychotic responders and resistant cases. Methods: Analysis of data from seven international cohorts of first-episode psychosis (FEP) with cognitive data at baseline (N = 683) and follow-up data on antipsychotic treatment response: 605 treatment responsive and 78 treatment resistant cases. Cognitive measures were grouped into seven cognitive domains based on the pre-existing literature. We ran multiple imputation for missing data and used logistic regression to test for associations between cognitive performance at FEP and treatment resistant status at follow-up. Results: On average patients who were future classified as treatment resistant reported poorer performance across most cognitive domains at baseline. Univariate logistic regressions showed that antipsychotic treatment resistance cases had significantly poorer IQ/general cognitive functioning at FEP (OR = 0.70, p = .003). These findings remained significant after adjusting for additional variables in multivariable analyses (OR = 0.76, p = .049). Conclusions: Although replication in larger studies is required, it appears that deficits in IQ/general cognitive functioning at first episode are associated with future treatment resistance. Cognitive variables may be able to provide further insight into neurodevelopmental factors associated with treatment resistance or act as early predictors of treatment resistance, which could allow prompt identification of refractory illness and timely interventions

Activation of NF-κB and p300/CBP potentiates cancer chemoimmunotherapy through induction of MHC-I antigen presentation

Many cancers evade immune rejection by suppressing major histocompatibility class I (MHC-I) antigen processing and presentation (AgPP). Such cancers do not respond to immune checkpoint inhibitor therapies (ICIT) such as PD-1/PD-L1 [PD-(L)1] blockade. Certain chemotherapeutic drugs augment tumor control by PD-(L)1 inhibitors through potentiation of T-cell priming but whether and how chemotherapy enhances MHC-I-dependent cancer cell recognition by cytotoxic T cells (CTLs) is not entirely clear. We now show that the lysine acetyl transferases p300/CREB binding protein (CBP) control MHC-I AgPPM expression and neoantigen amounts in human cancers. Moreover, we found that two distinct DNA damaging drugs, the platinoid oxaliplatin and the topoisomerase inhibitor mitoxantrone, strongly up-regulate MHC-I AgPP in a manner dependent on activation of nuclear factor kappa B (NF-κB), p300/CBP, and other transcription factors, but independently of autocrine IFNγ signaling. Accordingly, NF-κB and p300 ablations prevent chemotherapy-induced MHC-I AgPP and abrogate rejection of low MHC-I-expressing tumors by reinvigorated CD8+ CTLs. Drugs like oxaliplatin and mitoxantrone may be used to overcome resistance to PD-(L)1 inhibitors in tumors that had "epigenetically down-regulated," but had not permanently lost MHC-I AgPP activity

Clinical predictors of antipsychotic treatment resistance: Development and internal validation of a prognostic prediction model by the STRATA-G consortium.

Our aim was to, firstly, identify characteristics at first-episode of psychosis that are associated with later antipsychotic treatment resistance (TR) and, secondly, to develop a parsimonious prediction model for TR. We combined data from ten prospective, first-episode psychosis cohorts from across Europe and categorised patients as TR or non-treatment resistant (NTR) after a mean follow up of 4.18 years (s.d. = 3.20) for secondary data analysis. We identified a list of potential predictors from clinical and demographic data recorded at first-episode. These potential predictors were entered in two models: a multivariable logistic regression to identify which were independently associated with TR and a penalised logistic regression, which performed variable selection, to produce a parsimonious prediction model. This model was internally validated using a 5-fold, 50-repeat cross-validation optimism-correction. Our sample consisted of N = 2216 participants of which 385 (17 %) developed TR. Younger age of psychosis onset and fewer years in education were independently associated with increased odds of developing TR. The prediction model selected 7 out of 17 variables that, when combined, could quantify the risk of being TR better than chance. These included age of onset, years in education, gender, BMI, relationship status, alcohol use, and positive symptoms. The optimism-corrected area under the curve was 0.59 (accuracy = 64 %, sensitivity = 48 %, and specificity = 76 %). Our findings show that treatment resistance can be predicted, at first-episode of psychosis. Pending a model update and external validation, we demonstrate the potential value of prediction models for TR

Transcription Factors Bind Thousands of Active and Inactive Regions in the Drosophila Blastoderm

Identifying the genomic regions bound by sequence-specific regulatory factors is central both to deciphering the complex DNA cis-regulatory code that controls transcription in metazoans and to determining the range of genes that shape animal morphogenesis. We used whole-genome tiling arrays to map sequences bound in Drosophila melanogaster embryos by the six maternal and gap transcription factors that initiate anterior–posterior patterning. We find that these sequence-specific DNA binding proteins bind with quantitatively different specificities to highly overlapping sets of several thousand genomic regions in blastoderm embryos. Specific high- and moderate-affinity in vitro recognition sequences for each factor are enriched in bound regions. This enrichment, however, is not sufficient to explain the pattern of binding in vivo and varies in a context-dependent manner, demonstrating that higher-order rules must govern targeting of transcription factors. The more highly bound regions include all of the over 40 well-characterized enhancers known to respond to these factors as well as several hundred putative new cis-regulatory modules clustered near developmental regulators and other genes with patterned expression at this stage of embryogenesis. The new targets include most of the microRNAs (miRNAs) transcribed in the blastoderm, as well as all major zygotically transcribed dorsal–ventral patterning genes, whose expression we show to be quantitatively modulated by anterior–posterior factors. In addition to these highly bound regions, there are several thousand regions that are reproducibly bound at lower levels. However, these poorly bound regions are, collectively, far more distant from genes transcribed in the blastoderm than highly bound regions; are preferentially found in protein-coding sequences; and are less conserved than highly bound regions. Together these observations suggest that many of these poorly bound regions are not involved in early-embryonic transcriptional regulation, and a significant proportion may be nonfunctional. Surprisingly, for five of the six factors, their recognition sites are not unambiguously more constrained evolutionarily than the immediate flanking DNA, even in more highly bound and presumably functional regions, indicating that comparative DNA sequence analysis is limited in its ability to identify functional transcription factor targets

Transcriptomic Profiling of Virus-Host Cell Interactions following Chicken Anaemia Virus (CAV) Infection in an In Vivo Model.

Chicken Anaemia Virus (CAV) is an economically important virus that targets lymphoid and erythroblastoid progenitor cells leading to immunosuppression. This study aimed to investigate the interplay between viral infection and the host's immune response to better understand the pathways that lead to CAV-induced immunosuppression. To mimic vertical transmission of CAV in the absence of maternally-derived antibody, day-old chicks were infected and their responses measured at various time-points post-infection by qRT-PCR and gene expression microarrays. The kinetics of mRNA expression levels of signature cytokines of innate and adaptive immune responses were determined by qRT-PCR. The global gene expression profiles of mock-infected (control) and CAV-infected chickens at 14 dpi were also compared using a chicken immune-related 5K microarray. Although in the thymus there was evidence of induction of an innate immune response following CAV infection, this was limited in magnitude. There was little evidence of a Th1 adaptive immune response in any lymphoid tissue, as would normally be expected in response to viral infection. Most cytokines associated with Th1, Th2 or Treg subsets were down-regulated, except IL-2, IL-13, IL-10 and IFNγ, which were all up-regulated in thymus and bone marrow. From the microarray studies, genes that exhibited significant (greater than 1.5-fold, false discovery rate <0.05) changes in expression in thymus and bone marrow on CAV infection were mainly associated with T-cell receptor signalling, immune response, transcriptional regulation, intracellular signalling and regulation of apoptosis. Expression levels of a number of adaptor proteins, such as src-like adaptor protein (SLA), a negative regulator of T-cell receptor signalling and the transcription factor Special AT-rich Binding Protein 1 (SATB1), were significantly down-regulated by CAV infection, suggesting potential roles for these genes as regulators of viral infection or cell defence. These results extend our understanding of CAV-induced immunosuppression and suggest a global immune dysregulation following CAV infection

Interaction Testing and Polygenic Risk Scoring to Estimate the Association of Common Genetic Variants With Treatment Resistance in Schizophrenia

Importance: About 20% to 30% of people with schizophrenia have psychotic symptoms that do not respond adequately to first-line antipsychotic treatment. This clinical presentation, chronic and highly disabling, is known as treatment-resistant schizophrenia (TRS). The causes of treatment resistance and their relationships with causes underlying schizophrenia are largely unknown. Adequately powered genetic studies of TRS are scarce because of the difficulty in collecting data from well-characterized TRS cohorts. Objective: To examine the genetic architecture of TRS through the reassessment of genetic data from schizophrenia studies and its validation in carefully ascertained clinical samples. Design, Setting, and Participants: Two case-control genome-wide association studies (GWASs) of schizophrenia were performed in which the case samples were defined as individuals with TRS (n = 10 501) and individuals with non-TRS (n = 20 325). The differences in effect sizes for allelic associations were then determined between both studies, the reasoning being such differences reflect treatment resistance instead of schizophrenia. Genotype data were retrieved from the CLOZUK and Psychiatric Genomics Consortium (PGC) schizophrenia studies. The output was validated using polygenic risk score (PRS) profiling of 2 independent schizophrenia cohorts with TRS and non-TRS: a prevalence sample with 817 individuals (Cardiff Cognition in Schizophrenia [CardiffCOGS]) and an incidence sample with 563 individuals (Genetics Workstream of the Schizophrenia Treatment Resistance and Therapeutic Advances [STRATA-G]). Main Outcomes and Measures: GWAS of treatment resistance in schizophrenia. The results of the GWAS were compared with complex polygenic traits through a genetic correlation approach and were used for PRS analysis on the independent validation cohorts using the same TRS definition. Results: The study included a total of 85 490 participants (48 635 [56.9%] male) in its GWAS stage and 1380 participants (859 [62.2%] male) in its PRS validation stage. Treatment resistance in schizophrenia emerged as a polygenic trait with detectable heritability (1% to 4%), and several traits related to intelligence and cognition were found to be genetically correlated with it (genetic correlation, 0.41-0.69). PRS analysis in the CardiffCOGS prevalence sample showed a positive association between TRS and a history of taking clozapine (r2 = 2.03%; P = .001), which was replicated in the STRATA-G incidence sample (r2 = 1.09%; P = .04). Conclusions and Relevance: In this GWAS, common genetic variants were differentially associated with TRS, and these associations may have been obscured through the amalgamation of large GWAS samples in previous studies of broadly defined schizophrenia. Findings of this study suggest the validity of meta-analytic approaches for studies on patient outcomes, including treatment resistance