Efficient Test and Visualization of Multi-Set Intersections

A proteome-wide analysis unveils a core Epstein-Barr virus antibody signature of classic Hodgkin lymphoma across ethnically diverse populations

Epstein-Barr virus (EBV) is an oncogenic virus associated with various malignancies, including classical Hodgkin lymphoma (cHL). Despite its known association, the specific role of humoral immune response to EBV remains poorly characterized in cHL. To address this, we conducted a study using a custom protein microarray to measure the antibody responses in cHL patients and matched healthy controls recruited from an East-Asian hospital-based case-control study. We identified 16 IgG antibodies significantly elevated in EBV-positive cHL compared with controls, defining an "East-Asian antibody signature of EBV-positive cHL." We evaluated responses against these 16 antibodies in a distinct European population, leveraging data from our previous European cHL case-control study from the UK, Denmark, and Sweden. A subset of antibodies (14/16, 87.5%) from the "East-Asian antibody signature of EBV-positive cHL" exhibited significant associations with cHL in the European population. Conversely, we assessed the "European antibody signature of EBV-positive cHL" identified in our prior study which consisted of 18 EBV antibodies (2 IgA, 16 IgG), in the East-Asian population. A subset of these antibodies (15/18, 83.3%) maintained significant associations with cHL in the East-Asian population. This cross-comparison of antibody signatures underscores the robust generalizability of EBV antibodies across populations. Five anti-EBV IgG antibodies (LMP-1, TK, BALF2, BDLF3, and BBLF1), found in both population-specific antibody signatures, represent a "core signature of EBV-positive cHL." Our findings suggest that the antibody responses targeting these core EBV proteins reflect a specific EBV gene expression pattern, serving as potential biomarkers for EBV-positive cHL independent of population-specific factors.

Assessing the global prevalence of wild birds in trade

Trade represents a significant threat to many wild species and is often clandestine and poorly monitored. Information on which species are most prevalent in trade and potentially threatened by it therefore remains fragmentary. We used 7 global data sets on birds in trade to identify species or groups of species at particular risk and assessed the extent to which they were congruent in terms of the species recorded in trade. We used the frequency with which species were recorded in the data sets as the basis for a trade prevalence score that was applied to all bird species globally. Literature searches and questionnaire surveys were used to develop a list of species known to be heavily traded to validate the trade prevalence score. The score was modeled to identify significant predictors of trade. Although the data sets sampled different parts of the broad trade spectrum, congruence among them was statistically strong in all comparisons. Furthermore, the frequency with which species were recorded within data sets was positively correlated with their occurrence across data sets, indicating that the trade prevalence score captured information on trade volume. The trade prevalence score discriminated well between species identified from semi-independent assessments as heavily or unsustainably traded and all other species. Globally, 45.1% of all bird species and 36.7% of globally threatened bird species had trade prevalence scores ≥1. Species listed in Appendices I or II of CITES, species with large geographical distributions, and nonpasserines tended to have high trade prevalence scores. Speciose orders with high mean trade prevalence scores included Falconiformes, Psittaciformes, Accipitriformes, Anseriformes, Bucerotiformes, and Strigiformes. Despite their low mean prevalence score, Passeriformes accounted for the highest overall number of traded species of any order but had low representation in CITES appendices. Geographical hotspots where large numbers of traded species co-occur differed among passerines (Southeast Asia and Eurasia) and nonpasserines (central South America, sub-Saharan Africa, and India). This first attempt to quantify and map the relative prevalence in trade of all bird species globally can be used to identify species and groups of species that may be at particular risk of harm from trade and can inform conservation and policy interventions to reduce its adverse impacts.

Patient subtyping analysis of baseline multi-omic data reveals distinct pre-immune states associated with antibody response to seasonal influenza vaccination

Understanding the molecular mechanisms underpinning diverse vaccination responses is critical for developing efficient vaccines. Molecular subtyping can offer insights into heterogeneous nature of responses and aid in vaccine design. We analyzed multi-omic data from 62 haemagglutinin seasonal influenza vaccine recipients (2019-2020), including transcriptomics, proteomics, glycomics, and metabolomics data collected pre-vaccination. We performed a subtyping analysis on the integrated data revealing five subtypes with distinct molecular signatures. These subtypes differed in the expression of pre-existing adaptive or innate immunity signatures, which were linked to significant variation in baseline immunoglobulin A (IgA) and hemagglutination inhibition (HAI) titer levels. It is worth noting that these differences persisted through day 28 post-vaccination, indicating the effect of initial immune state on vaccination response. These findings highlight the significance of interpersonal variation in baseline immune status as a crucial factor in determining the effectiveness of seasonal vaccines. Ultimately, incorporating molecular profiling could enable personalized vaccine optimization.

Multi-omics analyses highlight molecular differences between clinical and neuropathological diagnoses in Alzheimer's disease

Both clinical diagnosis and neuropathological diagnosis are commonly used in literature to categorize individuals as Alzheimer's disease (AD) or non-AD in omics analyses. Whether these diagnostic strategies result in distinct profiles of molecular abnormalities is poorly understood. Here, we analysed one of the most commonly used AD omics datasets in the literature from the Religious Orders Study and Memory and Aging Project (ROSMAP) cohort and compared the two diagnosis strategies using brain transcriptome and metabolome by grouping individuals as non-AD and AD according to clinical or neuropathological diagnosis separately. Differentially expressed genes, associated pathways related with AD hallmarks and AD-related genes showed that the categorization based on neuropathological diagnosis more accurately reflects the disease state at the molecular level than the categorization based on clinical diagnosis. We further identified consensus biomarker candidates between the two diagnosis strategies such as 5-hydroxylysine, sphingomyelin and 1-myristoyl-2-palmitoyl-GPC as metabolite biomarkers and sphingolipid metabolism as a pathway biomarker, which could be robust AD biomarkers since they are independent of diagnosis strategies. We also used consensus AD and consensus non-AD individuals between the two diagnostic strategies to train a machine-learning based model, which we used to classify the individuals who were cognitively normal but diagnosed as AD based on neuropathological diagnosis (asymptomatic AD individuals). The majority of these individuals were classified as consensus AD patients for both omics data types. Our study provides a detailed characterization of both diagnostic strategies in terms of the association of the corresponding multi-omics profiles with AD.

Comparative Transcriptomic Reveals Greater Similarities in Response to Temperature Than to Invasive Alien Predator in the Damselfly Ischnura elegans Across Different Geographic Scales

The impact of global changes on populations may not be necessarily uniform across a species' range. Here, we aim at comparing the phenotypic and transcriptomic response to warming and an invasive predator cue in populations across different geographic scales in the damselfly Ischnura elegans. We collected adult females in two ponds in southern Poland (central latitude) and two ponds in southern Sweden (high latitude). We raised their larvae in growth chambers and exposed them to combination of temperature and a predator cue released by the crayfish Orconectes limosus. When larvae reached the prefinal larval stage, they were phenotyped for traits related to growth and size and collected for a gene expression analysis. High-latitude populations exhibited greater phenotypic and transcriptomic variation than central-latitude populations. Across latitudes and ponds, temperature generally increased growth rate and the predator cue decreased mass, but the effects of temperature were also pond-specific. Comparison of the transcriptomic profiles revealed a greater overlap in the response to temperature across latitudes and ponds, especially for pathway-related oxidative stress and sugar and lipid metabolism. The transcriptomic response to a predator cue and to the interaction temperature × predator cue was more pond-specific and overlapped only for few genes and pathways related to cuticle, development and signal transduction. We demonstrated that central- and high-latitude populations may partially respond through similar mechanisms to warming and, to a lower extent to a predator cue and to the interaction temperature × predator cue. For the predator cue and the interaction, the large fraction of ponds-specific genes suggests local adaptation. We show that high-latitude populations were generally more plastic at the phenotypic and transcriptomic level and may be more capable to cope with environmental changes than their central-latitude counterparts.

Genomics of ecological adaptation in Canary Island Descurainia (Brassicaceae) and comparisons with other Brassicaceae

Oceanic archipelagos provide striking examples of lineages that have radiated over pronounced ecological gradients. Accompanying this diversification, lineages have evolved adaptations allowing survival in extreme environments. Here, we investigate the genomic basis of ecological adaptation in Canary Island Descurainia (Brassicaceae), an island relative of Arabidopsis. The seven endemic species have diversified in situ along an elevational and ecological gradient, from low-elevation scrub to high-elevation sub-alpine desert. We first generated a reference genome for Descurainia millefolia, phylogenetic analysis of which placed it as sister to D. sophioides. Ninety-six gene families were found to be specific to D. millefolia and a further 1087 and 1469 gene families have expanded or contracted in size, respectively, along the D. millefolia branch. We then employed genome re-sequencing to sample 14 genomes across the seven species of Canary Island Descurainia and an outgroup. Phylogenomic analyses were consistent with previous reconstructions of Canary Island Descurainia in resolving low- and high-elevation clades. Using the branch-site dN/dS method, we detected positive selection for 275 genes on the branch separating the low- and high-elevation species and these positively selected genes (PSGs) were significantly enriched for functions related to reproduction and stress tolerance. Comparing PSGs to those in analyses of adaptation to elevation and/or latitude in other Brassicaceae, we found little evidence of widespread convergence and gene reuse, except for two examples, one of which was a significant overlap between Descurainia and Draba nivalis, a species restricted to high latitudes. The study of Canary Island Descurainia suggests that the transition to high-elevation environments such as that found in the high mountains of the Canary Islands involves selection on genes related to reproduction and stress tolerance but that repeated evolution across different lineages that have evolved into similar habitats is limited, indicating substantially different molecular trajectories to adaptation.

Selection against domestication alleles in introduced rabbit populations

Humans have moved domestic animals around the globe for thousands of years. These have occasionally established feral populations in nature, often with devastating ecological consequences. To understand how natural selection shapes re-adaptation into the wild, we investigated one of the most successful colonizers in history, the European rabbit. By sequencing the genomes of 297 rabbits across three continents, we show that introduced populations exhibit a mixed wild-domestic ancestry. We show that alleles that increased in frequency during domestication were preferentially selected against in novel natural environments. Interestingly, causative mutations for common domestication traits sometimes segregate at considerable frequencies if associated with less drastic phenotypes (for example, coat colour dilution), whereas mutations that are probably strongly maladaptive in nature are absent. Whereas natural selection largely targeted different genomic regions in each introduced population, some of the strongest signals of parallelism overlap genes associated with neuronal or brain function. This limited parallelism is probably explained by extensive standing genetic variation resulting from domestication together with the complex mixed ancestry of introduced populations. Our findings shed light on the selective and molecular mechanisms that enable domestic animals to re-adapt to the wild and provide important insights for the mitigation and management of invasive populations.

scCross: a deep generative model for unifying single-cell multi-omics with seamless integration, cross-modal generation, and in silico exploration

Single-cell multi-omics data reveal complex cellular states, providing significant insights into cellular dynamics and disease. Yet, integration of multi-omics data presents challenges. Some modalities have not reached the robustness or clarity of established transcriptomics. Coupled with data scarcity for less established modalities and integration intricacies, these challenges limit our ability to maximize single-cell omics benefits. We introduce scCross, a tool leveraging variational autoencoders, generative adversarial networks, and the mutual nearest neighbors (MNN) technique for modality alignment. By enabling single-cell cross-modal data generation, multi-omics data simulation, and in silico cellular perturbations, scCross enhances the utility of single-cell multi-omics studies.

Clinical and in vitro models identify distinct adaptations enhancing Staphylococcus aureus pathogenesis in human macrophages

Staphylococcus aureus is a facultative intracellular pathogen of human macrophages, which facilitates chronic infection. The genotypes, pathways, and mutations influencing that phenotype remain incompletely explored. Here, we used two distinct strategies to ascertain S. aureus gene mutations affecting pathogenesis in macrophages. First, we analyzed isolates collected serially from chronic cystic fibrosis (CF) respiratory infections. We found that S. aureus strains evolved greater macrophage invasion capacity during chronic human infection. Bacterial genome-wide association studies (GWAS) identified 127 candidate genes for which mutation was significantly associated with macrophage pathogenesis in vivo. In parallel, we passaged laboratory S. aureus strains in vitro to select for increased infection of human THP-1 derived macrophages, which identified 15 candidate genes by whole-genome sequencing. Functional validation of candidate genes using isogenic transposon mutant knockouts and CRISPR interference (CRISPRi) knockdowns confirmed virulence contributions from 37 of 39 tested genes (95%) implicated by in vivo studies and 7 of 10 genes (70%) ascertained from in vitro selection, with one gene in common to the two strategies. Validated genes included 17 known virulence factors (39%) and 27 newly identified by our study (61%), some encoding functions not previously associated with macrophage pathogenesis. Most genes (80%) positively impacted macrophage invasion when disrupted, consistent with the phenotype readily arising from loss-of-function mutations in vivo. This work reveals genes and mechanisms that contribute to S. aureus infection of macrophages, highlights differences in mutations underlying convergent phenotypes arising from in vivo and in vitro systems, and supports the relevance of S. aureus macrophage pathogenesis during chronic respiratory infection in CF. Additional studies will be needed to illuminate the exact mechanisms by which implicated mutations affect their phenotypes.

Targeted suppression of siRNA biogenesis in Arabidopsis pollen promotes triploid seed viability

In plants, small-interfering RNAs (siRNAs) mediate epigenetic silencing via the RNA-directed DNA methylation (RdDM) pathway, which is particularly prominent during reproduction and seed development. However, there is limited understanding of the origins and dynamics of reproductive siRNAs acting in different cellular and developmental contexts. Here, we used the RNaseIII-like protein RTL1 to suppress siRNA biogenesis in Arabidopsis pollen, and found distinct siRNA subsets produced during pollen development. We demonstrate that RTL1 expression in the late microspore and vegetative cell strongly impairs epigenetic silencing, and resembles RdDM mutants in their ability to bypass interploidy hybridization barriers in the seed. However, germline-specific RTL1 expression did not impact transgenerational inheritance of triploid seed lethality. These results reveal the existence of multiple siRNA subsets accumulated in mature pollen, and suggest that mobile siRNAs involved in the triploid block are produced in germline precursor cells after meiosis, or in the vegetative cell during pollen mitosis.

Meta-analysis towards FSHD reveals misregulation of neuromuscular junction, nuclear envelope, and spliceosome

Facioscapulohumeral muscular dystrophy (FSHD) is one of the most common autosomal dominant muscle disorders, yet no cure or amelioration exists. The clinical presentation is diverse, making it difficult to identify the actual driving pathomechanism among many downstream events. To unravel this complexity, we performed a meta-analysis of 13 original omics datasets (in total 171 FSHD and 129 control samples). Our approach confirmed previous findings about the disease pathology and specified them further. We confirmed increased expression of former proposed DUX4 biomarkers, and furthermore impairment of the respiratory chain. Notably, the meta-analysis provides insights about so far not reported pathways, including misregulation of neuromuscular junction protein encoding genes, downregulation of the spliceosome, and extensive alterations of nuclear envelope protein expression. Finally, we developed a publicly available shiny app to provide a platform for researchers who want to search our analysis for genes of interest in the future.

Single-cell multi-cohort dissection of the schizophrenia transcriptome

The complexity and heterogeneity of schizophrenia have hindered mechanistic elucidation and the development of more effective therapies. Here, we performed single-cell dissection of schizophrenia-associated transcriptomic changes in the human prefrontal cortex across 140 individuals in two independent cohorts. Excitatory neurons were the most affected cell group, with transcriptional changes converging on neurodevelopment and synapse-related molecular pathways. Transcriptional alterations included known genetic risk factors, suggesting convergence of rare and common genomic variants on neuronal population-specific alterations in schizophrenia. Based on the magnitude of schizophrenia-associated transcriptional change, we identified two populations of individuals with schizophrenia marked by expression of specific excitatory and inhibitory neuronal cell states. This single-cell atlas links transcriptomic changes to etiological genetic risk factors, contextualizing established knowledge within the human cortical cytoarchitecture and facilitating mechanistic understanding of schizophrenia pathophysiology and heterogeneity.

Lifespan analysis of repeat expression reveals age-dependent upregulation of HERV-K in the neurotypical human brain

DNA repetitive sequences (or repeats) comprise over 50% of the human genome and have a crucial regulatory role, specifically regulating transcription machinery. The human brain is the tissue with the highest detectable repeat expression and dysregulations on the repeat activity are related to several neurological and neurodegenerative disorders, as repeat-derived products can stimulate a pro-inflammatory response. Even so, it is unclear how repeat expression acts on the aging neurotypical brain. Here, we leverage a large postmortem transcriptome cohort spanning the human lifespan to assess global repeat expression in the neurotypical brain. We identified 21,696 differentially expressed repeats (DERs) that varied across seven age bins (Prenatal; 0-15; 16-29; 30-39; 40-49; 50-59; 60+) across the caudate nucleus (n=271), dorsolateral prefrontal cortex (n=304), and hippocampus (n=310). Interestingly, we found that long interspersed nuclear elements and long terminal repeats (LTRs) DERs were the most abundant repeat families when comparing infants to early adolescence (0-15) with older adults (60+). Of these differentially regulated LTRs, we identified 17 shared across all brain regions, including increased expression of HERV-K-int in older adult brains (60+). Co-expression analysis from each of the three brain regions also showed repeats from the HERV subfamily were intramodular hubs in its subnetworks. While we do not observe a strong global relationship between repeat expression and age, we identified HERV-K as a repeat signature associated with the aging neurotypical brain. Our study is the first global assessment of repeat expression in the neurotypical brain.

Targeted accurate RNA consensus sequencing (tARC-seq) reveals mechanisms of replication error affecting SARS-CoV-2 divergence

RNA viruses, like SARS-CoV-2, depend on their RNA-dependent RNA polymerases (RdRp) for replication, which is error prone. Monitoring replication errors is crucial for understanding the virus's evolution. Current methods lack the precision to detect rare de novo RNA mutations, particularly in low-input samples such as those from patients. Here we introduce a targeted accurate RNA consensus sequencing method (tARC-seq) to accurately determine the mutation frequency and types in SARS-CoV-2, both in cell culture and clinical samples. Our findings show an average of 2.68 × 10-5 de novo errors per cycle with a C > T bias that cannot be solely attributed to APOBEC editing. We identified hotspots and cold spots throughout the genome, correlating with high or low GC content, and pinpointed transcription regulatory sites as regions more susceptible to errors. tARC-seq captured template switching events including insertions, deletions and complex mutations. These insights shed light on the genetic diversity generation and evolutionary dynamics of SARS-CoV-2.

Comparing Transcriptomic Points of Departure to Apical Effect Concentrations For Larval Fathead Minnow Exposed to Chemicals with Four Different Modes Of Action

It is postulated that below a transcriptomic-based point of departure, adverse effects are unlikely to occur, thereby providing a chemical concentration to use in screening level hazard assessment. The present study extends previous work describing a high-throughput fathead minnow assay that can provide full transcriptomic data after exposure to a test chemical. One-day post-hatch fathead minnows were exposed to ten concentrations of three representatives of four chemical modes of action: organophosphates, ecdysone receptor agonists, plant photosystem II inhibitors, and estrogen receptor agonists for 24 h. Concentration response modeling was performed on whole body gene expression data from each exposure, using measured chemical concentrations when available. Transcriptomic points of departure in larval fathead minnow were lower than apical effect concentrations across fish species but not always lower than toxic effect concentrations in other aquatic taxa like crustaceans and insects. The point of departure was highly dependent on measured chemical concentration which were often lower than the nominal concentration. Differentially expressed genes between chemicals within modes of action were compared and often showed statistically significant overlap. In addition, reproducibility between identical exposures using a positive control chemical (CuSO4) and variability associated with the transcriptomic point of departure using in silico sampling were considered. Results extend a transcriptomic-compatible fathead minnow high-throughput assay for possible use in ecological hazard screening.

Molecular signatures of angiogenesis inhibitors: a single-embryo untargeted metabolomics approach in zebrafish

Angiogenesis is a key process in embryonic development, a disruption of this process can lead to severe developmental defects, such as limb malformations. The identification of molecular perturbations representative of antiangiogenesis in zebrafish embryo (ZFE) may guide the assessment of developmental toxicity from an endpoint- to a mechanism-based approach, thereby improving the extrapolation of findings to humans. Thus, the aim of the study was to discover molecular changes characteristic of antiangiogenesis and developmental toxicity. We exposed ZFEs to two antiangiogenic drugs (SU4312, sorafenib) and two developmental toxicants (methotrexate, rotenone) with putative antiangiogenic action. Molecular changes were measured by performing untargeted metabolomics in single embryos. The metabolome response was accompanied by the occurrence of morphological alterations. Two distinct metabolic effect patterns were observed. The first pattern comprised common effects of two specific angiogenesis inhibitors and the known teratogen methotrexate, strongly suggesting a shared mode of action of antiangiogenesis and developmental toxicity. The second pattern involved joint effects of methotrexate and rotenone, likely related to disturbances in energy metabolism. The metabolites of the first pattern, such as phosphatidylserines, pterines, retinol, or coenzyme Q precursors, represented potential links to antiangiogenesis and related developmental toxicity. The metabolic effect pattern can contribute to biomarker identification for a mechanism-based toxicological testing.

Spatially Resolved Transcriptomic Signatures of Hippocampal Subregions and Arc-Expressing Ensembles in Active Place Avoidance Memory

The rodent hippocampus is a spatially organized neuronal network that supports the formation of spatial and episodic memories. We conducted bulk RNA sequencing and spatial transcriptomics experiments to measure gene expression changes in the dorsal hippocampus following the recall of active place avoidance (APA) memory. Through bulk RNA sequencing, we examined the gene expression changes following memory recall across the functionally distinct subregions of the dorsal hippocampus. We found that recall induced differentially expressed genes (DEGs) in the CA1 and CA3 hippocampal subregions were enriched with genes involved in synaptic transmission and synaptic plasticity, while DEGs in the dentate gyrus (DG) were enriched with genes involved in energy balance and ribosomal function. Through spatial transcriptomics, we examined gene expression changes following memory recall across an array of spots encompassing putative memory-associated neuronal ensembles marked by the expression of the IEGs Arc, Egr1, and c-Jun. Within samples from both trained and untrained mice, the subpopulations of spatial transcriptomic spots marked by these IEGs were transcriptomically and spatially distinct from one another. DEGs detected between Arc+ and Arc- spots exclusively in the trained mouse were enriched in several memory-related gene ontology terms, including "regulation of synaptic plasticity" and "memory." Our results suggest that APA memory recall is supported by regionalized transcriptomic profiles separating the CA1 and CA3 from the DG, transcriptionally and spatially distinct IEG expressing spatial transcriptomic spots, and biological processes related to synaptic plasticity as a defining the difference between Arc+ and Arc- spatial transcriptomic spots.

Early results of the integrative epigenomic-transcriptomic landscape of colorectal adenoma and cancer

BACKGROUND

Aberrant methylation is common during the initiation and progression of colorectal cancer (CRC), and detecting these changes that occur during early adenoma (ADE) formation and CRC progression has clinical value.

AIM

To identify potential DNA methylation markers specific to ADE and CRC.

METHODS

Here, we performed SeqCap targeted bisulfite sequencing and RNA-seq analysis of colorectal ADE and CRC samples to profile the epigenomic-transcriptomic landscape.

RESULTS

Comparing 22 CRC and 25 ADE samples, global methylation was higher in the former, but both showed similar methylation patterns regarding differentially methylated gene positions, chromatin signatures, and repeated elements. High-grade CRC tended to exhibit elevated methylation levels in gene promoter regions compared to those in low-grade CRC. Combined with RNA-seq gene expression data, we identified 14 methylation-regulated differentially expressed genes, of which only AGTR1 and NECAB1 methylation had prognostic significance.

CONCLUSION

Our results suggest that genome-wide alterations in DNA methylation occur during the early stages of CRC and demonstrate the methylation signatures associated with colorectal ADEs and CRC, suggesting prognostic biomarkers for CRC.

Patient Subtyping Analysis of Baseline Multi-omic Data Reveals Distinct Pre-immune States Predictive of Vaccination Responses

Understanding the molecular mechanisms that underpin diverse vaccination responses is a critical step toward developing efficient vaccines. Molecular subtyping approaches can offer valuable insights into the heterogeneous nature of responses and aid in the design of more effective vaccines. In order to explore the molecular signatures associated with the vaccine response, we analyzed baseline transcriptomics data from paired samples of whole blood, proteomics and glycomics data from serum, and metabolomics data from urine, obtained from influenza vaccine recipients (2019-2020 season) prior to vaccination. After integrating the data using a network-based model, we performed a subtyping analysis. The integration of multiple data modalities from 62 samples resulted in five baseline molecular subtypes with distinct molecular signatures. These baseline subtypes differed in the expression of pre-existing adaptive or innate immunity signatures, which were linked to significant variation across subtypes in baseline immunoglobulin A (IgA) and hemagglutination inhibition (HAI) titer levels. It is worth noting that these significant differences persisted through day 28 post-vaccination, indicating the effect of initial immune state on vaccination response. These findings highlight the significance of interpersonal variation in baseline immune status as a crucial factor in determining vaccine response and efficacy. Ultimately, incorporating molecular profiling could enable personalized vaccine optimization.

Crosstalk between chromatin and Shavenbaby defines transcriptional output along the Drosophila intestinal stem cell lineage

The transcription factor Shavenbaby (Svb), the only member of the OvoL family in Drosophila, controls the fate of various epithelial embryonic cells and adult stem cells. Post-translational modification of Svb produces two protein isoforms, Svb-ACT and Svb-REP, which promote adult intestinal stem cell renewal or differentiation, respectively. To define Svb mode of action, we used engineered cell lines and develop an unbiased method to identify Svb target genes across different contexts. Within a given cell type, Svb-ACT and Svb-REP antagonistically regulate the expression of a set of target genes, binding specific enhancers whose accessibility is constrained by chromatin landscape. Reciprocally, Svb-REP can influence local chromatin marks of active enhancers to help repressing target genes. Along the intestinal lineage, the set of Svb target genes progressively changes, together with chromatin accessibility. We propose that Svb-ACT-to-REP transition promotes enterocyte differentiation of intestinal stem cells through direct gene regulation and chromatin remodeling.

Molecular profiling of human substantia nigra identifies diverse neuron types associated with vulnerability in Parkinson's disease

Parkinson's disease (PD) is characterized pathologically by the loss of dopaminergic (DA) neurons in the substantia nigra (SN). Whether cell types beyond DA neurons in the SN show vulnerability in PD remains unclear. Through transcriptomic profiling of 315,867 high-quality single nuclei in the SN from individuals with and without PD, we identified cell clusters representing various neuron types, glia, endothelial cells, pericytes, fibroblasts, and T cells and investigated cell type-dependent alterations in gene expression in PD. Notably, a unique neuron cluster marked by the expression of RIT2, a PD risk gene, also displayed vulnerability in PD. We validated RIT2-enriched neurons in midbrain organoids and the mouse SN. Our results demonstrated distinct transcriptomic signatures of the RIT2-enriched neurons in the human SN and implicated reduced RIT2 expression in the pathogenesis of PD. Our study sheds light on the diversity of cell types, including DA neurons, in the SN and the complexity of molecular and cellular changes associated with PD pathogenesis.

Serum Proteomic Signatures of Common Health Outcomes among Older Adults

INTRODUCTION

In aging populations, the coexistence of multiple health comorbidities represents a significant challenge for clinicians and researchers. Leveraging advances in omics techniques to characterize these health conditions may provide insight into disease pathogenesis as well as reveal biomarkers for monitoring, prognostication, and diagnosis. Researchers have previously established the utility of big data approaches with respect to comprehensive health outcome measurements in younger populations, identifying protein markers that may provide significant health information with a single blood sample.

METHODS

Here, we employed a similar approach in two cohorts of older adults, the Baltimore Longitudinal Study of Aging (mean age = 76.12 years) and InCHIANTI Study (mean age = 66.05 years), examining the relationship between levels of serum proteins and 5 key health outcomes: kidney function, fasting glucose, physical activity, lean body mass, and percent body fat.

RESULTS

Correlations between proteins and health outcomes were primarily shared across both older adult cohorts. We further identified that most proteins associated with health outcomes in the older adult cohorts were not associated with the same outcomes in a prior study of a younger population. A subset of proteins, adiponectin, MIC-1, and NCAM-120, were associated with at least three health outcomes in both older adult cohorts but not in the previously published younger cohort, suggesting that they may represent plausible markers of general health in older adult populations.

CONCLUSION

Taken together, these findings suggest that comprehensive protein health markers have utility in aging populations and are distinct from those identified in younger adults, indicating unique mechanisms of disease with aging.

Proteomic Signaling of Dual-Specificity Phosphatase 4 (DUSP4) in Alzheimer's Disease

DUSP4 is a member of the DUSP (dual-specificity phosphatase) subfamily that is selective to the mitogen-activated protein kinases (MAPK) and has been implicated in a range of biological processes and functions in Alzheimer's disease (AD). In this study, we utilized the stereotactic delivery of adeno-associated virus (AAV)-DUSP4 to overexpress DUSP4 in the dorsal hippocampus of 5xFAD and wildtype (WT) mice, then used mass spectrometry (MS)-based proteomics along with the label-free quantification to profile the proteome and phosphoproteome in the hippocampus. We identified protein expression and phosphorylation patterns modulated in 5xFAD mice and examined the sex-specific impact of DUSP4 overexpression on the 5xFAD proteome/phosphoproteome. In 5xFAD mice, a substantial number of proteins were up- or down-regulated in both male and female mice in comparison to age and sex-matched WT mice, many of which are involved in AD-related biological processes, such as activated immune response or suppressed synaptic activities. Many proteins in pathways, such as immune response were found to be suppressed in response to DUSP4 overexpression in male 5xFAD mice. In contrast, such a shift was absent in female mice. For the phosphoproteome, we detected an array of phosphorylation sites regulated in 5xFAD compared to WT and modulated via DUSP4 overexpression in each sex. Interestingly, 5xFAD- and DUSP4-associated phosphorylation changes occurred in opposite directions. Strikingly, both the 5xFAD- and DUSP4-associated phosphorylation changes were found to be mostly in neurons and play key roles in neuronal processes and synaptic functions. Site-centric pathway analysis revealed that both the 5xFAD- and DUSP4-associated phosphorylation sites were enriched for a number of kinase sets in females but only a limited number of sets of kinases in male mice. Taken together, our results suggest that male and female 5xFAD mice responded to DUSP4 overexpression via shared and sex-specific molecular mechanisms, which might underly similar reductions in amyloid pathology in both sexes while learning deficits were reduced in only females with DUSP4 overexpression. Finally, we validated our findings with the sex-specific AD-associated proteomes in human cohorts and further developed DUSP4-centric proteomic network models and signaling maps for each sex.

Histone acetylation in an Alzheimer's disease cell model promotes homeostatic amyloid-reducing pathways

Alzheimer's Disease (AD) is a disorder characterized by cognitive decline, neurodegeneration, and accumulation of amyloid plaques and tau neurofibrillary tangles in the brain. Dysregulation of epigenetic histone modifications may lead to expression of transcriptional programs that play a role either in protecting against disease genesis or in worsening of disease pathology. One such histone modification, acetylation of histone H3 lysine residue 27 (H3K27ac), is primarily localized to genomic enhancer regions and promotes active gene transcription. We previously discovered H3K27ac to be more abundant in AD patient brain tissue compared to the brains of age-matched non-demented controls. In this study, we use iPSC-neurons derived from familial AD patients with an amyloid precursor protein (APP) duplication (APPDup neurons) as a model to study the functional effect of lowering CBP/P300 enzymes that catalyze H3K27ac. We found that homeostatic amyloid-reducing genes were upregulated in the APPDup neurons compared to non-demented controls. We lowered CBP/P300 to reduce H3K27ac, which led to decreased expression of numerous of these homeostatic amyloid-reducing genes, along with increased extracellular secretion of a toxic amyloid-β species, Aβ(1-42). Our findings suggest that epigenomic histone acetylation, including H3K27ac, drives expression of compensatory genetic programs in response to AD-associated insults, specifically those resulting from APP duplication, and thus may play a role in mitigating AD pathology in neurons.

Effects of Age and Exposure Duration on the Sensitivity of Early Life Stage Fathead Minnow (Pimephales promelas) to Waterborne Propranolol Exposure

Propranolol is a heavily prescribed, nonspecific beta-adrenoceptor (bAR) antagonist frequently found in wastewater effluents, prompting concern over its potential to adversely affect exposed organisms. In the present study, the transcriptional responses of 4, 5, and 6 days postfertilization (dpf) ±1 h fathead minnow, exposed for 6, 24, or 48 h to 0.66 or 3.3 mg/L (nominal) propranolol were characterized using RNA sequencing. The number of differentially expressed genes (DEGs) was used as an estimate of sensitivity. A trend toward increased sensitivity with age was observed; fish >7 dpf at the end of exposure were particularly sensitive to propranolol. The DEGs largely overlapped among treatment groups, suggesting a highly consistent response that was independent of age. Cluster analysis was performed using normalized count data for unexposed and propranolol-exposed fish. Control fish clustered tightly by age, with fish ≥7 dpf clustering away from younger fish, reflecting developmental differences. When clustering was conducted using exposed fish, in cases where propranolol induced a minimal or no transcriptional response, the results mirrored those of the control fish and did not appreciably cluster by treatment. In treatment groups that displayed a more robust transcriptional response, the effects of propranolol were evident; however, fish <7 dpf clustered away from older fish, despite having similar numbers of DEGs. Increased sensitivity at 7 dpf coincided with developmental milestones with the potential to alter propranolol pharmacokinetics or pharmacodynamics, such as the onset of exogenous feeding and gill functionality as well as increased systemic expression of bAR. These results may have broader implications because toxicity testing often utilizes fish <4 dpf, prior to the onset of these potentially important developmental milestones, which may result in an underestimation of risk for some chemicals. Environ Toxicol Chem 2024;00:1-14. Published 2023. This article is a U.S. Government work and is in the public domain in the USA.

Multiscale protein networks systematically identify aberrant protein interactions and oncogenic regulators in seven cancer types

Global proteomic data generated by advanced mass spectrometry (MS) technologies can help bridge the gap between genome/transcriptome and functions and hold great potential in elucidating unbiased functional models of pro-tumorigenic pathways. To this end, we collected the high-throughput, whole-genome MS data and conducted integrative proteomic network analyses of 687 cases across 7 cancer types including breast carcinoma (115 tumor samples; 10,438 genes), clear cell renal carcinoma (100 tumor samples; 9,910 genes), colorectal cancer (91 tumor samples; 7,362 genes), hepatocellular carcinoma (101 tumor samples; 6,478 genes), lung adenocarcinoma (104 tumor samples; 10,967 genes), stomach adenocarcinoma (80 tumor samples; 9,268 genes), and uterine corpus endometrial carcinoma UCEC (96 tumor samples; 10,768 genes). Through the protein co-expression network analysis, we identified co-expressed protein modules enriched for differentially expressed proteins in tumor as disease-associated pathways. Comparison with the respective transcriptome network models revealed proteome-specific cancer subnetworks associated with heme metabolism, DNA repair, spliceosome, oxidative phosphorylation and several oncogenic signaling pathways. Cross-cancer comparison identified highly preserved protein modules showing robust pan-cancer interactions and identified endoplasmic reticulum-associated degradation (ERAD) and N-acetyltransferase activity as the central functional axes. We further utilized these network models to predict pan-cancer protein regulators of disease-associated pathways. The top predicted pan-cancer regulators including RSL1D1, DDX21 and SMC2, were experimentally validated in lung, colon, breast cancer and fetal kidney cells. In summary, this study has developed interpretable network models of cancer proteomes, showcasing their potential in unveiling novel oncogenic regulators, elucidating underlying mechanisms, and identifying new therapeutic targets.

Exploring multisite heterogeneity of human basal cell carcinoma proteome and transcriptome

Basal cell carcinoma (BCC) is the most common type of skin cancer. Due to multiple, potential underlying molecular tumor aberrations, clinical treatment protocols are not well-defined. This study presents multisite molecular heterogeneity profiles of human BCC based on RNA and proteome profiling. Three areas from lesions excised from 9 patients were analyzed. The focus was gene expression profiles based on proteome and RNA measurements of intra-tumor heterogeneity from the same patient and inter-tumor heterogeneity in nodular, infiltrative, and superficial BCC tumor subtypes from different patients. We observed significant overlap in intra- and inter-tumor variability of proteome and RNA expression profiles, showing significant multisite heterogeneity of protein expression in the BCC tumors. Inter-subtype analysis has also identified unique proteins for each BCC subtype. This profiling leads to a deeper understanding of BCC molecular heterogeneity and potentially contributes to developing new sampling tools for personalized diagnostics therapeutic approaches to BCC.

Distinct and shared endothermic strategies in the heat producing tissues of tuna and other teleosts

Although most fishes are ectothermic, some, including tuna and billfish, achieve endothermy through specialized heat producing tissues that are modified muscles. How these heat producing tissues evolved, and whether they share convergent molecular mechanisms, remain unresolved. Here, we generated a high-quality genome from the mackerel tuna (Euthynnus affinis) and investigated the heat producing tissues of this fish by single-nucleus and bulk RNA sequencing. Compared with other teleosts, tuna-specific genetic variation is strongly associated with muscle differentiation. Single-nucleus RNA-seq revealed a high proportion of specific slow skeletal muscle cell subtypes in the heat producing tissues of tuna. Marker genes of this cell subtype are associated with the relative sliding of actin and myosin, suggesting that tuna endothermy is mainly based on shivering thermogenesis. In contrast, cross-species transcriptome analysis indicated that endothermy in billfish relies mainly on non-shivering thermogenesis. Nevertheless, the heat producing tissues of the different species do share some tissue-specific genes, including vascular-related and mitochondrial genes. Overall, although tunas and billfishes differ in their thermogenic strategies, they share similar expression patterns in some respects, highlighting the complexity of convergent evolution.

Proteomic signaling of dual specificity phosphatase 4 (DUSP4) in Alzheimer's disease

DUSP4 is a member of the DUSP (Dual-Specificity Phosphatase) subfamily that is selective to the mitogen-activated protein kinases (MAPK) and has been implicated in a range of biological processes and functions in Alzheimer's disease (AD). In this study, we utilized stereotactic delivery of adeno-associated virus (AAV)-DUSP4 to overexpress DUSP4 in the dorsal hippocampus of 5xFAD and wildtype (WT) mice, then used mass spectrometry (MS)-based proteomics along with label-free quantification to profile the proteome and phosphoproteome in the hippocampus. We identified patterns of protein expression and phosphorylation that are modulated in 5xFAD mice and examined the sex-specific impact of DUSP4 overexpression on the 5xFAD proteome/phosphoproteome. In 5xFAD mice, a substantial number of proteins were up- or down-regulated in both male and female mice in comparison to age and sex-matched WT mice, many of which are involved in AD-related biological processes, such as the activated immune response or suppression of synaptic activities. Upon DUSP4 overexpression, significantly regulated proteins were found in pathways that were suppressed, such as the immune response, in male 5xFAD mice. In contrast, such a shift was absent in female mice. For the phosphoproteome, we detected an array of phosphorylation sites that are regulated in 5xFAD compared to WT, and are modulated by DUSP4 overexpression in each sex. Interestingly, the changes in 5xFAD- and DUSP4-associated phosphorylation occurred in opposite directions. Strikingly, both the 5xFAD- and DUSP4-associated phosphorylation changes were found for the most part in neurons, and play key roles in neuronal processes and synaptic function. Site-centric pathway analysis revealed that both the 5xFAD- and DUSP4-associated phosphorylation sites were enriched for a number of kinase sets in female, but only a limited number of sets of kinases in male mice. Taken together, our results suggest that male and female 5xFAD mice respond to DUSP4 overexpression via shared and sex-specific molecular mechanisms, which might underly similar reductions in amyloid pathology in both sexes, while learning deficits were reduced in only females with DUSP4 overexpression. Finally, we validated our findings with the sex-specific AD-associated proteomes in human cohorts and further developed DUSP4-centric proteomic network models and signaling maps for each sex.

Analysis of asymptomatic Drosophila models for ALS and SMA reveals convergent impact on functional protein complexes linked to neuro-muscular degeneration

BACKGROUND

Spinal Muscular Atrophy (SMA) and Amyotrophic Lateral Sclerosis (ALS) share phenotypic and molecular commonalities, including the fact that they can be caused by mutations in ubiquitous proteins involved in RNA metabolism, namely SMN, TDP-43 and FUS. Although this suggests the existence of common disease mechanisms, there is currently no model to explain the resulting motor neuron dysfunction. In this work we generated a parallel set of Drosophila models for adult-onset RNAi and tagged neuronal expression of the fly orthologues of the three human proteins, named Smn, TBPH and Caz, respectively. We profiled nuclear and cytoplasmic bound mRNAs using a RIP-seq approach and characterized the transcriptome of the RNAi models by RNA-seq. To unravel the mechanisms underlying the common functional impact of these proteins on neuronal cells, we devised a computational approach based on the construction of a tissue-specific library of protein functional modules, selected by an overall impact score measuring the estimated extent of perturbation caused by each gene knockdown.

RESULTS

Transcriptome analysis revealed that the three proteins do not bind to the same RNA molecules and that only a limited set of functionally unrelated transcripts is commonly affected by their knock-down. However, through our integrative approach we were able to identify a concerted effect on protein functional modules, albeit acting through distinct targets. Most strikingly, functional annotation revealed that these modules are involved in critical cellular pathways for motor neurons, including neuromuscular junction function. Furthermore, selected modules were found to be significantly enriched in orthologues of human neuronal disease genes.

CONCLUSIONS

The results presented here show that SMA and ALS disease-associated genes linked to RNA metabolism functionally converge on neuronal protein complexes, providing a new hypothesis to explain the common motor neuron phenotype. The functional modules identified represent promising biomarkers and therapeutic targets, namely given their alteration in asymptomatic settings.

Proteomic signaling of dual specificity phosphatase 4 (DUSP4) in Alzheimer's disease

DUSP4 is a member of the DUSP (Dual-Specificity Phosphatase) subfamily that is selective to the mitogen-activated protein kinases (MAPK) and has been implicated in a range of biological processes and functions in Alzheimer's disease (AD). In this study, we utilized stereotactic delivery of adeno-associated virus (AAV)-DUSP4 to overexpress DUSP4 in the dorsal hippocampus of 5×FAD and wildtype (WT) mice, then used mass spectrometry (MS)-based proteomics along with label-free quantification to profile the proteome and phosphoproteome in the hippocampus. We identified patterns of protein expression and phosphorylation that are modulated in 5×FAD mice and examined the sex-specific impact of DUSP4 overexpression on the 5×FAD proteome/phosphoproteome. In 5×FAD mice, a substantial number of proteins were up- or down-regulated in both male and female mice in comparison to age and sex-matched WT mice, many of which are involved in AD-related biological processes, such as the activated immune response or suppression of synaptic activities. Upon DUSP4 overexpression, significantly regulated proteins were found in pathways that were suppressed, such as the immune response, in male 5×FAD mice. In contrast, such a shift was absent in female mice. For the phosphoproteome, we detected an array of phosphorylation sites that are regulated in 5×FAD compared to WT, and are modulated by DUSP4 overexpression in each sex. Interestingly, the changes in 5×FAD- and DUSP4-associated phosphorylation occurred in opposite directions. Strikingly, both the 5×FAD- and DUSP4-associated phosphorylation changes were found for the most part in neurons, and play key roles in neuronal processes and synaptic function. Site-centric pathway analysis revealed that both the 5×FAD- and DUSP4-associated phosphorylation sites were enriched for a number of kinase sets in female, but only a limited number of sets of kinases in male mice. Taken together, our results suggest that male and female 5×FAD mice respond to DUSP4 overexpression via shared and sex-specific molecular mechanisms, which might underly similar reductions in amyloid pathology in both sexes, while learning deficits were reduced in only females with DUSP4 overexpression. Finally, we validated our findings with the sex-specific AD-associated proteomes in human cohorts and further developed DUSP4-centric proteomic network models and signaling maps for each sex.

Parallel adaptation to lower altitudes is associated with enhanced plasticity in Heliosperma pusillum (Caryophyllaceae)

High levels of phenotypic plasticity are thought to be inherently costly in stable or extreme environments, but enhanced plasticity may evolve as a response to new environments and foster novel phenotypes. Heliosperma pusillum forms glabrous alpine and pubescent montane ecotypes that diverged recurrently and polytopically (parallel evolution) and can serve as evolutionary replicates. The specific alpine and montane localities are characterized by distinct temperature conditions, available moisture, and light. Noteworthy, the ecotypes show a home-site fitness advantage in reciprocal transplantations. To disentangle the relative contribution of constitutive versus plastic gene expression to altitudinal divergence, we analyze the transcriptomic profiles of two parallely evolved ecotype pairs, grown in reciprocal transplantations at native altitudinal sites. In this incipient stage of divergence, only a minor proportion of genes appear constitutively differentially expressed between the ecotypes in both pairs, regardless of the growing environment. Both derived, montane populations bear comparatively higher plasticity of gene expression than the alpine populations. Genes that change expression plastically or constitutively underlie similar ecologically relevant pathways, related to response to drought and trichome formation. Other relevant processes, such as photosynthesis, rely mainly on plastic changes. The enhanced plasticity consistently observed in the montane ecotype likely evolved as a response to the newly colonized, drier, and warmer niche. We report a striking parallelism of directional changes in gene expression plasticity. Thus, plasticity appears to be a key mechanism shaping the initial stages of phenotypic evolution, likely fostering adaptation to novel environments.

Epigenetic signature of human immune aging in the GESTALT study

Age-associated DNA methylation in blood cells convey information on health status. However, the mechanisms that drive these changes in circulating cells and their relationships to gene regulation are unknown. We identified age-associated DNA methylation sites in six purified blood-borne immune cell types (naive B, naive CD4+ and CD8+ T cells, granulocytes, monocytes, and NK cells) collected from healthy individuals interspersed over a wide age range. Of the thousands of age-associated sites, only 350 sites were differentially methylated in the same direction in all cell types and validated in an independent longitudinal cohort. Genes close to age-associated hypomethylated sites were enriched for collagen biosynthesis and complement cascade pathways, while genes close to hypermethylated sites mapped to neuronal pathways. In silico analyses showed that in most cell types, the age-associated hypo- and hypermethylated sites were enriched for ARNT (HIF1β) and REST transcription factor (TF) motifs, respectively, which are both master regulators of hypoxia response. To conclude, despite spatial heterogeneity, there is a commonality in the putative regulatory role with respect to TF motifs and histone modifications at and around these sites. These features suggest that DNA methylation changes in healthy aging may be adaptive responses to fluctuations of oxygen availability.

Intrinsically dysregulated cellular stress signaling genes and gene networks in postpartum depression

Postpartum depression (PPD) is a leading cause of morbidity and mortality among women. Clinically, the administration and withdrawal of supraphysiologic estradiol and progesterone (E2 + P) can cause affective symptom reoccurrence in women with a history of PPD, but not matched controls. To investigate the cellular basis underlying this differential affective response, lymphoblastoid cell lines (LCLs) were derived from women with and without past PPD and compared transcriptomically in hormone conditions mimicking pregnancy and parturition: supraphysiologic E2 + P-addback; supraphysiologic E2 + P-withdrawal; and no added E2 + P (Baseline). RNA-sequencing identified unique differentially expressed genes (DEGs) in all hormone conditions, but the majority tended to be downregulated in PPD and observed in E2 + P-addback. Two of these DEGs were evolutionarily conserved cellular stress regulators: IMPACT, an integrative response protein maintaining translational homeostasis, and WWTR1, a transcriptional coactivator in the 'Hippo' pathway mediating cell proliferation and survival. Correspondingly, significant gene network modules were linked to cell cycle progression, estrogen response, and immune dysregulation, suggesting innate differences in intracellular signaling in PPD. In certain hormone conditions, PPD LCLs displayed increased GATA3 expression (an upstream regulator of IMPACT and WWTR1) and differentially phosphorylated eiF2α (the ultimate downstream target of IMPACT). Taken together, these transcriptomic data primarily implicate innately dysregulated cellular responses as potentially influencing mood and/or escalating PPD risk. Furthermore, the intrinsic downregulation of IMPACT's translation and WWTR1's transcription networks may suggest a novel link between PPD and a compromised ability to maintain homeostasis in the context of cellular stress occurring during pregnancy and parturition.

Medium-throughput image-based phenotypic siRNA screen to unveil the molecular basis of B cell polarization

Cell polarity is an essential and highly conserved process governing cell function. Cell polarization is generally triggered by an external signal that induces the relocation of the centrosome, thus defining the polarity axis of the cell. Here, we took advantage of B cells as a model to study cell polarity and perform a medium-throughput siRNA-based imaging screen to identify new molecular regulators of polarization. We first identified candidates based on a quantitative proteomic analysis of proteins differentially associated with the centrosome of resting non-polarized and stimulated polarized B cells. We then targeted 233 candidates in a siRNA screen and identified hits regulating the polarization of the centrosome and/or lysosomes in B cells upon stimulation. Our dataset of proteomics, images, and polarity indexes provides a valuable source of information for a broad community of scientists interested in the molecular mechanisms regulating cell polarity.

An Ancestral Balanced Inversion Polymorphism Confers Global Adaptation

Since the pioneering work of Dobzhansky in the 1930s and 1940s, many chromosomal inversions have been identified, but how they contribute to adaptation remains poorly understood. In Drosophila melanogaster, the widespread inversion polymorphism In(3R)Payne underpins latitudinal clines in fitness traits on multiple continents. Here, we use single-individual whole-genome sequencing, transcriptomics, and published sequencing data to study the population genomics of this inversion on four continents: in its ancestral African range and in derived populations in Europe, North America, and Australia. Our results confirm that this inversion originated in sub-Saharan Africa and subsequently became cosmopolitan; we observe marked monophyletic divergence of inverted and noninverted karyotypes, with some substructure among inverted chromosomes between continents. Despite divergent evolution of this inversion since its out-of-Africa migration, derived non-African populations exhibit similar patterns of long-range linkage disequilibrium between the inversion breakpoints and major peaks of divergence in its center, consistent with balancing selection and suggesting that the inversion harbors alleles that are maintained by selection on several continents. Using RNA-sequencing, we identify overlap between inversion-linked single-nucleotide polymorphisms and loci that are differentially expressed between inverted and noninverted chromosomes. Expression levels are higher for inverted chromosomes at low temperature, suggesting loss of buffering or compensatory plasticity and consistent with higher inversion frequency in warm climates. Our results suggest that this ancestrally tropical balanced polymorphism spread around the world and became latitudinally assorted along similar but independent climatic gradients, always being frequent in subtropical/tropical areas but rare or absent in temperate climates.

Multiscale network analysis identifies potential receptors for SARS-CoV-2 and reveals their tissue-specific and age-dependent expression

The coronavirus disease 2019 (COVID-19) pandemic has affected tens of millions of individuals and caused hundreds of thousands of deaths worldwide. Here, we present a comprehensive, multiscale network analysis of the transcriptional response to the virus. In particular, we focused on key regulators, cell receptors, and host processes that were hijacked by the virus for its advantage. ACE2-controlled processes involved CD300e (a TYROBP receptor) as a key regulator and the activation of IL-2 pro-inflammatory cytokine signaling. We further investigated the age dependency of such receptors in different tissues. In summary, this study provides novel insights into the gene regulatory organization during the SARS-CoV-2 infection and the tissue-specific, age-dependent expression of the cell receptors involved in COVID-19.

Dependency of NELF-E-SLUG-KAT2B epigenetic axis in breast cancer carcinogenesis

Cancer cells undergo transcriptional reprogramming to drive tumor progression and metastasis. Using cancer cell lines and patient-derived tumor organoids, we demonstrate that loss of the negative elongation factor (NELF) complex inhibits breast cancer development through downregulating epithelial-mesenchymal transition (EMT) and stemness-associated genes. Quantitative multiplexed Rapid Immunoprecipitation Mass spectrometry of Endogenous proteins (qPLEX-RIME) further reveals a significant rewiring of NELF-E-associated chromatin partners as a function of EMT and a co-option of NELF-E with the key EMT transcription factor SLUG. Accordingly, loss of NELF-E leads to impaired SLUG binding on chromatin. Through integrative transcriptomic and genomic analyses, we identify the histone acetyltransferase, KAT2B, as a key functional target of NELF-E-SLUG. Genetic and pharmacological inactivation of KAT2B ameliorate the expression of EMT markers, phenocopying NELF ablation. Elevated expression of NELF-E and KAT2B is associated with poorer prognosis in breast cancer patients, highlighting the clinical relevance of our findings. Taken together, we uncover a crucial role of the NELF-E-SLUG-KAT2B epigenetic axis in breast cancer carcinogenesis.

Mimicked Mixing-Induced Heterogeneities of Industrial Bioreactors Stimulate Long-Lasting Adaption Programs in Ethanol-Producing Yeasts

Commercial-scale bioreactors create an unnatural environment for microbes from an evolutionary point of view. Mixing insufficiencies expose individual cells to fluctuating nutrient concentrations on a second-to-minute scale while transcriptional and translational capacities limit the microbial adaptation time from minutes to hours. This mismatch carries the risk of inadequate adaptation effects, especially considering that nutrients are available at optimal concentrations on average. Consequently, industrial bioprocesses that strive to maintain microbes in a phenotypic sweet spot, during lab-scale development, might suffer performance losses when said adaptive misconfigurations arise during scale-up. Here, we investigated the influence of fluctuating glucose availability on the gene-expression profile in the industrial yeast Ethanol Red™. The stimulus-response experiment introduced 2 min glucose depletion phases to cells growing under glucose limitation in a chemostat. Even though Ethanol Red™ displayed robust growth and productivity, a single 2 min depletion of glucose transiently triggered the environmental stress response. Furthermore, a new growth phenotype with an increased ribosome portfolio emerged after complete adaptation to recurring glucose shortages. The results of this study serve a twofold purpose. First, it highlights the necessity to consider the large-scale environment already at the experimental development stage, even when process-related stressors are moderate. Second, it allowed the deduction of strain engineering guidelines to optimize the genetic background of large-scale production hosts.

G. de Lomana A, Gomes-Filho JV, Vêncio RZN, Moritz RL, Koide T, Baliga NS. A Genome-Scale Atlas Reveals Complex Interplay of Transcription and Translation in an Archaeon

The scale of post-transcriptional regulation and the implications of its interplay with other forms of regulation in environmental acclimation are underexplored for organisms of the domain Archaea. Here, we have investigated the scale of post-transcriptional regulation in the extremely halophilic archaeon Halobacterium salinarum NRC-1 by integrating the transcriptome-wide locations of transcript processing sites (TPSs) and SmAP1 binding, the genome-wide locations of antisense RNAs (asRNAs), and the consequences of RNase_2099C knockout on the differential expression of all genes. This integrated analysis has discovered that 54% of all protein-coding genes in the genome of this haloarchaeon are likely targeted by multiple mechanisms for putative post-transcriptional processing and regulation, with about 20% of genes likely being regulated by combinatorial schemes involving SmAP1, asRNAs, and RNase_2099C. Comparative analysis of mRNA levels (transcriptome sequencing [RNA-Seq]) and protein levels (sequential window acquisition of all theoretical fragment ion spectra mass spectrometry [SWATH-MS]) for 2,579 genes over four phases of batch culture growth in complex medium generated additional evidence for the conditional post-transcriptional regulation of 7% of all protein-coding genes. We demonstrate that post-transcriptional regulation may act to fine-tune specialized and rapid acclimation to stressful environments, e.g., as a switch to turn on gas vesicle biogenesis to promote vertical relocation under anoxic conditions and modulate the frequency of transposition by insertion sequence (IS) elements of the IS200/IS605, IS4, and ISH3 families. Findings from this study are provided as an atlas in a public Web resource (https://halodata.systemsbiology.net). IMPORTANCE While the transcriptional regulation landscape of archaea has been extensively investigated, we currently have limited knowledge about post-transcriptional regulation and its driving mechanisms in this domain of life. In this study, we collected and integrated omics data from multiple sources and technologies to infer post-transcriptionally regulated genes and the putative mechanisms modulating their expression at the protein level in Halobacterium salinarum NRC-1. The results suggest that post-transcriptional regulation may drive environmental acclimation by regulating hallmark biological processes. To foster discoveries by other research groups interested in the topic, we extended our integrated data to the public in the form of an interactive atlas (https://halodata.systemsbiology.net).

The large GTPase AtGBPL3 links nuclear envelope formation and morphogenesis to transcriptional repression

Guanylate binding proteins (GBPs) are prominent regulators of immunity not known to be required for nuclear envelope formation and morphogenesis. Here we identify the Arabidopsis GBP orthologue AtGBPL3 as a lamina component with essential functions in mitotic nuclear envelope reformation, nuclear morphogenesis and transcriptional repression during interphase. AtGBPL3 is preferentially expressed in mitotically active root tips, accumulates at the nuclear envelope and interacts with centromeric chromatin as well as with lamina components transcriptionally repressing pericentromeric chromatin. Reduced expression of AtGBPL3 or associated lamina components similarly altered nuclear morphology and caused overlapping transcriptional deregulation. Investigating the dynamics of AtGBPL3-GFP and other nuclear markers during mitosis (1) revealed that AtGBPL3 accumulation on the surface of daughter nuclei precedes nuclear envelope reformation and (2) uncovered defects in this process in roots of AtGBPL3 mutants, which cause programmed cell death and impair growth. AtGBPL3 functions established by these observations are unique among dynamin-family large GTPases.

Integrating transcriptomic datasets across neurological disease identifies unique myeloid subpopulations driving disease-specific signatures

Microglia and bone marrow-derived monocytes are key elements of central nervous system (CNS) inflammation, both capable of enhancing and dampening immune-mediated pathology. However, the study-specific focus on individual cell types, disease models or experimental approaches has limited our ability to infer common and disease-specific responses. This meta-analysis integrates bulk and single-cell transcriptomic datasets of microglia and monocytes from disease models of autoimmunity, neurodegeneration, sterile injury, and infection to build a comprehensive resource connecting myeloid responses across CNS disease. We demonstrate that the bulk microglial and monocyte program is highly contingent on the disease environment, challenging the notion of a universal microglial disease signature. Integration of six single-cell RNA-sequencing datasets revealed that these disease-specific signatures are likely driven by differing proportions of unique myeloid subpopulations that were individually expanded in different disease settings. These subsets were functionally-defined as neurodegeneration-associated, inflammatory, interferon-responsive, phagocytic, antigen-presenting, and lipopolysaccharide-responsive cellular states, revealing a core set of myeloid responses at the single-cell level that are conserved across CNS pathology. Showcasing the predictive and practical value of this resource, we performed differential expression analysis on microglia and monocytes across disease and identified Cd81 as a new neuroinflammatory-stable gene that accurately identified microglia and distinguished them from monocyte-derived cells across all experimental models at both the bulk and single-cell level. Together, this resource dissects the influence of disease environment on shared immune response programmes to build a unified perspective of myeloid behavior across CNS pathology.

Convergent and complementary selection shaped gains and losses of eusociality in sweat bees

Sweat bees have repeatedly gained and lost eusociality, a transition from individual to group reproduction. Here we generate chromosome-length genome assemblies for 17 species and identify genomic signatures of evolutionary trade-offs associated with transitions between social and solitary living. Both young genes and regulatory regions show enrichment for these molecular patterns. We also identify loci that show evidence of complementary signals of positive and relaxed selection linked specifically to the convergent gains and losses of eusociality in sweat bees. This includes two pleiotropic proteins that bind and transport juvenile hormone (JH)-a key regulator of insect development and reproduction. We find that one of these proteins is primarily expressed in subperineurial glial cells that form the insect blood-brain barrier and that brain levels of JH vary by sociality. Our findings are consistent with a role of JH in modulating social behaviour and suggest that eusocial evolution was facilitated by alteration of the proteins that bind and transport JH, revealing how an ancestral developmental hormone may have been co-opted during one of life's major transitions. More broadly, our results highlight how evolutionary trade-offs have structured the molecular basis of eusociality in these bees and demonstrate how both directional selection and release from constraint can shape trait evolution.

Polygenic routes lead to parallel altitudinal adaptation in Heliosperma pusillum (Caryophyllaceae)

Understanding how organisms adapt to the environment is a major goal of modern biology. Parallel evolution-the independent evolution of similar phenotypes in different populations-provides a powerful framework to investigate the evolutionary potential of populations, the constraints of evolution, its repeatability and therefore its predictability. Here, we quantified the degree of gene expression and functional parallelism across replicated ecotype formation in Heliosperma pusillum (Caryophyllaceae), and gained insights into the architecture of adaptive traits. Population structure analyses and demographic modelling support a previously formulated hypothesis of parallel polytopic divergence of montane and alpine ecotypes. We detect a large proportion of differentially expressed genes (DEGs) underlying divergence within each replicate ecotype pair, with a strikingly low number of shared DEGs across pairs. Functional enrichment of DEGs reveals that the traits affected by significant expression divergence are largely consistent across ecotype pairs, in strong contrast to the nonshared genetic basis. The remarkable redundancy of differential gene expression indicates a polygenic architecture for the diverged adaptive traits. We conclude that polygenic traits appear key to opening multiple routes for adaptation, widening the adaptive potential of organisms.

Prolonged partner separation erodes nucleus accumbens transcriptional signatures of pair bonding in male prairie voles

The loss of a spouse is often cited as the most traumatic event in a person's life. However, for most people, the severity of grief and its maladaptive effects subside over time via an understudied adaptive process. Like humans, socially monogamous prairie voles (Microtus ochrogaster) form opposite-sex pair bonds, and upon partner separation, show stress phenotypes that diminish over time. We test the hypothesis that extended partner separation diminishes pair bond-associated behaviors and causes pair bond transcriptional signatures to erode. Opposite-sex or same-sex paired males were cohoused for 2 weeks and then either remained paired or were separated for 48 hours or 4 weeks before collecting fresh nucleus accumbens tissue for RNAseq. In a separate cohort, we assessed partner-directed affiliation at these time points. We found that these behaviors persist despite prolonged separation in both same-sex and opposite-sex paired voles. Opposite-sex pair bonding led to changes in accumbal transcription that were stably maintained while animals remained paired but eroded following prolonged partner separation. Eroded genes are associated with gliogenesis and myelination, suggesting a previously undescribed role for glia in pair bonding and loss. Further, we pioneered neuron-specific translating ribosomal affinity purification in voles. Neuronally enriched transcriptional changes revealed dopaminergic-, mitochondrial-, and steroid hormone signaling-associated gene clusters sensitive to acute pair bond disruption and loss adaptation. Our results suggest that partner separation erodes transcriptomic signatures of pair bonding despite core behavioral features of the bond remaining intact, revealing potential molecular processes priming a vole to be able to form a new bond.

Convergent abnormalities in striatal gene networks in human cocaine use disorder and mouse cocaine administration models

Cocaine use disorder (CUD) is an intractable syndrome, and rising overdose death rates represent a substantial public health crisis that exacts tremendous personal and financial costs on patients and society. Sharp increases in cocaine use drive the urgent need for better mechanistic insight into this chronic relapsing brain disorder that currently lacks effective treatment options. To investigate the transcriptomic changes involved, we conducted RNA sequencing on two striatal brain regions that are heavily implicated in CUD, the nucleus accumbens and caudate nucleus, from men suffering from CUD and matched controls. Weighted gene coexpression analyses identified CUD-specific gene networks enriched in ionotropic receptors and linked to lowered neuroinflammation, contrasting the proinflammatory responses found in opioid use disorder. Integration of comprehensive transcriptomic datasets from mouse cocaine self-administration models revealed evolutionarily conserved gene networks in CUD that implicate especially D1 medium spiny neurons as drivers of cocaine-induced plasticity.

Genome-wide signatures of the geographic expansion and breeding of soybean

Soybean is a leguminous crop that provides oil and protein. Exploring the genomic signatures of soybean evolution is crucial for breeding varieties with improved adaptability to environmental extremes. We analyzed the genome sequences of 2,214 soybeans and proposed a soybean evolutionary route, i.e., the expansion of annual wild soybean (Glycine soja Sieb. & Zucc.) from southern China and its domestication in central China, followed by the expansion and local breeding selection of its landraces (G. max (L.) Merr.). We observed that the genetic introgression in soybean landraces was mostly derived from sympatric rather than allopatric wild populations during the geographic expansion. Soybean expansion and breeding were accompanied by the positive selection of flowering time genes, including GmSPA3c. Our study sheds light on the evolutionary history of soybean and provides valuable genetic resources for its future breeding.

Dietary stress remodels the genetic architecture of lifespan variation in outbred Drosophila

Evolutionary theory suggests that lifespan-reducing alleles should be purged from the gene pool, and yet decades of genome-wide association and model organism studies have shown that they persist. One potential explanation is that alleles that regulate lifespan do so only in certain environmental contexts. We exposed outbred Drosophila to control and high-sugar diets and genotyped more than 10,000 adult flies to track allele frequency changes over the course of a single adult lifespan. We identified thousands of lifespan-associated alleles associated with early versus late-life trade-offs, late-onset effects and genotype-by-environment interactions. Remarkably, a third of lifespan-associated genetic variation had environmentally dependent effects on lifespan. We find that lifespan-reducing alleles are often recently derived, have stronger effects on a high-sugar diet and show signatures of selection in wild Drosophila populations, consistent with the evolutionary mismatch hypothesis. Our results provide insight into the highly polygenic and context-dependent genetic architecture of lifespan variation and the evolutionary processes that shape this key trait.

Causal effect of COVID-19 on Alzheimer's disease: A Mendelian randomization study

It was reported that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection may cause brain size reduction and cognitive decline. Whether COVID-19 may contribute to the development of Alzheimer's disease (AD) is not known. We conducted genetic correlation and Mendelian randomization (MR) analyses to assess genetic relationships and potential causal associations between AD and three COVID-19 outcomes (SARS-CoV-2 infection, COVID-19 hospitalization, and critical COVID-19) by utilizing genome-wide association study datasets on these traits. A map of COVID-19-driven molecular pathways was constructed to investigate potential mechanisms underlying the COVID-19 and AD connection. Genetic correlation analyses indicated that AD had a significant positive genetic correlation with hospitalized COVID-19 (rg  = 0.271). The MR analysis from the inverse-variance-weighted model showed that genetic liabilities to hospitalized COVID-19 (odds ratio: 1.02, 95% confidence interval: 1.01-1.03) and critical COVID-19 (1.01, 1.00-1.02) were associated with an increased risk for AD. However, no causal effect of genetic liability to SARS-CoV-2 infection on AD was detected (1.03, 0.97-1.09). A total of 60 functionally interconnected genes were reported to mediate the COVID-19-AD connection, which showed functional enrichment in immunity-related pathways and tissue enrichment in the lung and brain. Our study suggests that severe COVID-19 may contribute to the development of AD, while suffering a mild case of COVID-19 may not increase the risk for AD. The influence of COVID-19 on AD may be mediated by immunity-related pathways acting predominantly in the lung and brain.

Contamination of aquatic environment with anticancer reagents influences Daphnia magna - Ecotoxicogenomics approach

Pharmaceuticals used in human medicine contaminate freshwater ecosystems. Chemotherapeutics applied in cancer treatment are found in freshwaters at low concentrations (in the range of ng L^-1) which, however, can be toxic or mutagenic to aquatic organisms. The aim of this study was to determine the impact of the alkylating/crosslinking anticancer agents, cyclophosphamide (CP) and cisplatin (CDDP), at the concentration detected in water, on Daphnia magna life history, transcriptome, and proteome. This filter feeding cladoceran is an important member of the aquatic food webs controlling algal biomass and forming basic food for planktivorous fish. Here, observations of the D. magna growth rate, age at first reproduction, and the number of eggs produced were performed in the presence of CP or CDDP. The D. magna proteins and RNA were isolated and analysed by mass spectrometry and the mRNA-seq method, respectively. Five generations of contact with the pharmaceuticals in question significantly influenced the D. magna life history parameters with the growth rate and number of laid eggs decreased, whereas age at first reproduction was increased. A decrease in survivorship was observed when daphnids were exposed to CP. These changes are the result of modifications in the gene/transcript expression followed by differences in the proteome profile in comparison to the untreated control. The proteome changes were generally in accordance with the modified transcriptome. The ecotoxicogenomics approach makes it possible to get closer to a complete picture of the influence of CP and CDDP on Daphnia. We have gathered evidence that animals in the presence of anticancer pharmaceuticals attempt to cope with permanent stress by changing their proteome and transcriptome profile. Additionally, our analyses indicate that CDDP showed a stronger effect on tested organisms than CP.