Efficient Test and Visualization of Multi-Set Intersections

A personalized metabolic modelling approach through integrated analysis of RNA-Seq-based genomic variants and gene expression levels in Alzheimer's disease

Generating condition-specific metabolic models by mapping gene expression data to genome-scale metabolic models (GEMs) is a routine approach to elucidate disease mechanisms from a metabolic perspective. On the other hand, integrating variants that perturb enzyme functionality from the same RNA-seq data may enhance GEM accuracy, offering insights into genome-wide metabolic pathology. Our study pioneers the extraction of both transcriptomic and genomic data from the same RNA-seq data to reconstruct personalized metabolic models. We map genes with significantly higher load of pathogenic variants in Alzheimer's disease (AD) onto a human GEM together with the gene expression data. Comparative analysis of the resulting personalized patient metabolic models with the control models shows enhanced accuracy in detecting AD-associated metabolic pathways compared to the case where only expression data is mapped on the GEM. Besides, several otherwise would-be missed pathways are annotated in AD by considering the effect of genomic variants.

The intensity of the transcriptional response varies across infection with distinct viral strains in an insect host

Organisms respond to infectious agents through diverse immune strategies, and may need to cater a specific response to distinct pathogen challenges, such as various strains of a virus, to maximize fitness. Deformed wing virus (DWV) is one of the most damaging viruses of honey bees (Apis mellifera) across the globe, with variant DWV-B currently expanding at the expense of variant DWV-A. While previous research has characterized general host transcriptomic responses to viral exposure, host responses to different DWV strains have not been fully explored. Here, we performed experimental infections with the two dominant strains of DWV, A and B, as well as a mixed infection, and conducted transcriptomic analyses to compare differences in host molecular response to infection. We confirmed canonical anti-viral response to DWV infection, including upregulation of Toll pathway genes and the antimicrobial peptides abaecin and hymenoptaecin. Furthermore, our results suggest a potential role of aerobic glycolysis during viral infection in honey bees. DWV-A and mixed infections were associated with differential expression of a much larger number of host genes than infection with DWV-B. That DWV-B potentially elicits a reduced host immune response may provide a mechanistic explanation for its higher virulence and global emergence. Overall, this study provides the first evidence for strain-specific immune responses to DWV infection, and integrates these findings into the broader domain of insect immunity and host-pathogen dynamics.

Ecological divergence of sibling allopolyploid marsh orchids is associated with species specific plasticity and distinct fungal communities

Phenotypic plasticity, the dynamic adjustment of traits to environmental variations, is crucial for enabling species to exploit broader niches and withstand suboptimal conditions. This adaptability is particularly relevant for newly formed allopolyploids, which possess redundant gene copies and must become established in diverse environments distinct from their parents and other relatives. By evaluating gene expression and root mycobiome among two ecologically divergent sibling allopolyploid marsh orchids (Dactylorhiza majalis and D. traunsteineri) in reciprocal transplants at localities where both species are native, we aimed to understand the drivers of species persistence in the face of interspecific gene flow. Despite consistent abiotic differences characterising the alternative environments at each locality, the majority of gene expression differences between the allopolyploids appears to be plastic. Ecologically relevant processes, such as photosynthesis and transmembrane transport, include some genes that are differentially expressed between the two orchids regardless of the environment, while others change their activity plastically in one species or the other. This suggests that although plasticity helps define the specific ecological range of each sibling allopolyploid, it also mediates gene flow between them, thereby preventing differentiation. Extending our investigations to the root mycobiome, we uncover more diverse fungal communities for either species when grown in the environment with nutrient-poor soils, indicating that both abiotic and biotic factors drive the distribution of sibling marsh orchids. Altogether, our results indicate that plasticity can simultaneously promote diversification and homogenisation of lineages, influencing the establishment and persistence of recurrently formed allopolyploid species.

Permanent defects in renal medullary structure and function after reversal of urinary obstruction

Urinary obstruction causes injury to the renal medulla, impairing the ability to concentrate urine and increasing the risk of progressive kidney disease. However, the regenerative capacity of the renal medulla after reversal of obstruction is poorly understood. To investigate this, we developed a mouse model of reversible urinary obstruction. Despite robust regeneration and complete histological recovery of the renal medulla, these mice exhibited a permanent defect in urinary concentrating capacity. However, there were lasting changes in the composition, organization, and transcriptional profiles of epithelial, endothelial, and interstitial cells. Persistent inflammatory responses were also seen in patients with renal stone disease, but there were also adaptive responses to the increasingly hypoxic environment of the renal medulla that occurred only after reversal of obstruction. These findings indicate that while partial repair occurs after reversal of urinary obstruction, there are lasting structural and functional changes across all major cellular compartments of the renal medulla. These changes reflect shared and distinct responses to different renal medullary injuries in humans and mice.

Single-Cell Multi-Omics Profiling of Immune Cells Isolated from Atherosclerotic Plaques in Male ApoE Knockout Mice Exposed to Arsenic

BACKGROUND

Millions worldwide are exposed to elevated levels of arsenic that significantly increase their risk of developing atherosclerosis, a pathology primarily driven by immune cells. While the impact of arsenic on immune cell populations in atherosclerotic plaques has been broadly characterized, cellular heterogeneity is a substantial barrier to in-depth examinations of the cellular dynamics for varying immune cell populations.

OBJECTIVES

This study aimed to conduct single-cell multi-omics profiling of atherosclerotic plaques in apolipoprotein E knockout (ApoE-/-) mice to elucidate transcriptomic and epigenetic changes in immune cells induced by arsenic exposure.

METHODS

The ApoE-/- mice were fed a high-fat diet and were exposed to either 200  ppb arsenic in drinking water or a tap water control, and single-cell multi-omics profiling was performed on atherosclerotic plaque-resident immune cells. Transcriptomic and epigenetic changes in immune cells were analyzed within the same cell to understand the effects of arsenic exposure.

RESULTS

Our data revealed that the transcriptional profile of macrophages from arsenic-exposed mice were significantly different from that of control mice and that differences were subtype specific and associated with cell-cell interaction and cell fates. Additionally, our data suggest that differences in arsenic-mediated changes in chromosome accessibility in arsenic-exposed mice were statistically more likely to be due to factors other than random variation compared to their effects on the transcriptome, revealing markers of arsenic exposure and potential targets for intervention.

DISCUSSION

These findings in mice provide insights into how arsenic exposure impacts immune cell types in atherosclerosis, highlighting the importance of considering cellular heterogeneity in studying such effects. The identification of subtype-specific differences and potential intervention targets underscores the significance of understanding the molecular mechanisms underlying arsenic-induced atherosclerosis. Further research is warranted to validate these findings and explore therapeutic interventions targeting immune cell dysfunction in arsenic-exposed individuals. https://doi.org/10.1289/EHP14285.

A transcriptome atlas of zygotic and somatic embryogenesis in Norway spruce

Somatic embryogenesis (SE) is a powerful model system for studying embryo development and an important method for scaling up availability of elite and climate-adapted genetic material of Norway spruce (Picea abies L. Karst). However, there are several steps during the development of the somatic embryo (Sem) that are suboptimal compared to zygotic embryo (Zem) development. These differences are poorly understood and result in substantial yield losses during plant production, which limits cost-effective large-scale production of SE plants. This study presents a comprehensive data resource profiling gene expression during zygotic and somatic embryo development to support studies aiming to advance understanding of gene regulatory programmes controlling embryo development. Transcriptome expression patterns were analysed during zygotic embryogenesis (ZE) in Norway spruce, including separated samples of the female gametophytes and Zem, and at multiple stages during SE. Expression data from eight developmental stages of SE, starting with pro-embryogenic masses (PEMs) up until germination, revealed extensive modulation of the transcriptome between the early and mid-stage maturing embryos and at the transition of desiccated embryos to germination. Comparative analysis of gene expression changes during ZE and SE identified differences in the pattern of gene expression changes and functional enrichment of these provided insight into the associated biological processes. Orthologs of transcription factors known to regulate embryo development in angiosperms were differentially regulated during Zem and Sem development and in the different zygotic embryo tissues, providing clues to the differences in development observed between Zem and Sem. This resource represents the most comprehensive dataset available for exploring embryo development in conifers.

APOE from astrocytes restores Alzheimer's Aβ-pathology and DAM-like responses in APOE deficient microglia

The major genetic risk factor for Alzheimer's disease (AD), APOE4, accelerates beta-amyloid (Aβ) plaque formation, but whether this is caused by APOE expressed in microglia or astrocytes is debated. We express here the human APOE isoforms in astrocytes in an Apoe-deficient AD mouse model. This is not only sufficient to restore the amyloid plaque pathology but also induces the characteristic transcriptional pathological responses in Apoe-deficient microglia surrounding the plaques. We find that both APOE4 and the protective APOE2 from astrocytes increase fibrillar plaque deposition, but differentially affect soluble Aβ aggregates. Microglia and astrocytes show specific alterations in function of APOE genotype expressed in astrocytes. Our experiments indicate a central role of the astrocytes in APOE mediated amyloid plaque pathology and in the induction of associated microglia responses.

Dosage-sensitive maternal siRNAs determine hybridization success in Capsella

Hybrid seed failure arising from wide crosses between plant species is a recurring obstacle in plant breeding, impeding the transfer of desirable traits. This postzygotic reproductive barrier primarily occurs in the endosperm, a tissue that nourishes the embryo and functions similarly to the placenta in mammals. We found that incompatible seeds show a loss of DNA methylation and chromatin condensation in the endosperm, similar to seeds lacking maternal RNA polymerase IV activity. This similarity is linked to a decrease in small interfering RNAs in the endosperm (sirenRNAs), maternal RNA polymerase IV-dependent short interfering RNAs that regulate DNA methylation. Several AGAMOUS-like MADS-box transcription factor genes (AGLs), key regulators of endosperm development, are targeted by sirenRNAs in cis and in trans. This finding aligns with the enrichment of AGL target genes among deregulated genes. We propose that hybrid seed failure results from reduced maternal sirenRNAs combined with increased AGL expression, leading to abnormal gene regulation in the endosperm.

Polyploids broadly generate novel haplotypes from trans-specific variation in Arabidopsis arenosa and Arabidopsis lyrata

Polyploidy, the result of whole genome duplication (WGD), is widespread across the tree of life and is often associated with speciation and adaptability. It is thought that adaptation in autopolyploids (within-species polyploids) may be facilitated by increased access to genetic variation. This variation may be sourced from gene flow with sister diploids and new access to other tetraploid lineages, as well as from increased mutational targets provided by doubled DNA content. Here, we deconstruct in detail the origins of haplotypes displaying the strongest selection signals in established, successful autopolyploids, Arabidopsis lyrata and Arabidopsis arenosa. We see strong signatures of selection in 17 genes implied in meiosis, cell cycle, and transcription across all four autotetraploid lineages present in our expanded sampling of 983 sequenced genomes. Most prominent in our results is the finding that the tetraploid-characteristic haplotypes with the most robust signals of selection were completely absent in all diploid sisters. In contrast, the fine-scaled variant 'mosaics' in the tetraploids originated from highly diverse evolutionary sources. These include widespread novel reassortments of trans-specific polymorphism from diploids, new mutations, and tetraploid-specific inter-species hybridization-a pattern that is in line with the broad-scale acquisition and reshuffling of potentially adaptive variation in tetraploids.

iDOMO: identification of drug combinations via multi-set operations for treating diseases

Combination therapy has become increasingly important for treating complex diseases which often involve multiple pathways and targets. However, experimental screening of drug combinations is costly and time-consuming. The availability of large-scale transcriptomic datasets (e.g. CMap and LINCS) from in vitro drug treatment experiments makes it possible to computationally predict drug combinations with synergistic effects. Towards this end, we developed a computational approach, termed Identification of Drug Combinations via Multi-Set Operations (iDOMO), to predict drug synergy based on multi-set operations of drug and disease gene signatures. iDOMO quantifies the synergistic effect of a pair of drugs by taking into account the combination's beneficial and detrimental effects on treating a disease. We evaluated iDOMO, in a DREAM Challenge dataset with the matched, pre- and post-treatment gene expression data and cell viability information. We further evaluated the performance of iDOMO by concordance index and Spearman correlation on predicting the Highest Single Agency (HSA) synergy scores for four most common cancer types in two large-scale drug combination databases, showing that iDOMO  significantly outperformed two existing popular drug combination approaches including the Therapeutic Score and the SynergySeq Orthogonality Score. Application of iDOMO to triple-negative breast cancer (TNBC) identified drug pairs with potential synergistic effects, with the combination of trifluridine and monobenzone being the most synergistic. Our in vitro experiments confirmed that the top predicted drug combination exerted a significant synergistic effect in inhibiting TNBC cell growth. In summary, iDOMO is an effective method for the in silico screening of synergistic drug combinations and will be a valuable tool for the development of novel therapeutics for complex diseases.

Pre- and Postsynaptic MEF2C Promotes Experience-Dependent, Input-Specific Development of Cortical Layer 4 to Layer 2/3 Excitatory Synapses and Regulates Activity-Dependent Expression of Synaptic Cell Adhesion Molecules

Experience- and activity-dependent transcription is a candidate mechanism to mediate development and refinement of specific cortical circuits. Here, we demonstrate that the activity-dependent transcription factor myocyte enhancer factor 2C (MEF2C) is required in both presynaptic layer (L) 4 and postsynaptic L2/3 mouse (male and female) somatosensory (S1) cortical neurons for development of this specific synaptic connection. While postsynaptic deletion of Mef2c weakens L4 synaptic inputs, it has no effect on inputs from local L2/3, contralateral S1, or the ipsilateral frontal/motor cortex. Similarly, homozygous or heterozygous deletion of Mef2c in presynaptic L4 neurons weakens L4 to L2/3 excitatory synaptic inputs by decreasing presynaptic release probability. Postsynaptic MEF2C is specifically required during an early postnatal, experience-dependent, period for L4 to L2/3 synapse function, and expression of transcriptionally active MEF2C (MEF2C-VP16) rescues weak L4 to L2/3 synaptic strength in sensory-deprived mice. Together, these results suggest that experience- and/or activity-dependent transcriptional activation of MEF2C promotes development of L4 to L2/3 synapses. Additionally, MEF2C regulates the expression of many pre- and postsynaptic genes in postnatal cortical neurons. Interestingly, MEF2C was necessary for activity-dependent expression of many presynaptic genes, including those that function in transsynaptic adhesion and neurotransmitter release. This work provides mechanistic insight into the experience-dependent development of specific cortical circuits.

Integrative analyses of convergent adaptation in sympatric extremophile fishes

The evolution of independent lineages along replicated environmental transitions frequently results in convergent adaptation, yet the degree to which convergence is present across multiple levels of biological organization is often unclear. Additionally, inherent biases associated with shared ancestry and variation in selective regimes across geographic replicates often pose challenges for confidently identifying patterns of convergence. We investigated a system in which three species of poeciliid fishes sympatrically occur in a toxic spring rich in hydrogen sulfide (H2S) and an adjacent nonsulfidic stream to examine patterns of adaptive evolution across levels of biological organization. We found convergence in morphological and physiological traits and genome-wide patterns of gene expression among all three species. In addition, there were shared signatures of selection on genes encoding H2S toxicity targets in the mitochondrial genomes of each species. However, analyses of nuclear genomes revealed neither evidence for substantial genomic islands of divergence around genes involved in H2S toxicity and detoxification nor substantial congruence of strongly differentiated regions across population pairs. These non-convergent, heterogeneous patterns of genomic divergence may indicate that sulfide tolerance is highly polygenic, with shared allele frequency shifts present at many loci with small effects along the genome. Alternatively, H2S tolerance may involve substantial genetic redundancy, with non-convergent, lineage-specific variation at multiple loci along the genome underpinning similar changes in phenotypes and gene expression. Overall, we demonstrate variability in the extent of convergence across organizational levels and highlight the challenges of linking patterns of convergence across scales.

Repeated Shifts in Sociality Are Associated With Fine-tuning of Highly Conserved and Lineage-Specific Enhancers in a Socially Flexible Bee

Comparative genomic studies of social insects suggest that changes in gene regulation are associated with evolutionary transitions in social behavior, but the activity of predicted regulatory regions has not been tested empirically. We used self-transcribing active regulatory region sequencing, a high-throughput enhancer discovery tool, to identify and measure the activity of enhancers in the socially variable sweat bee, Lasioglossum albipes. We identified over 36,000 enhancers in the L. albipes genome from 3 social and 3 solitary populations. Many enhancers were identified in only a subset of L. albipes populations, revealing rapid divergence in regulatory regions within this species. Population-specific enhancers were often proximal to the same genes across populations, suggesting compensatory gains and losses of regulatory regions may preserve gene activity. We also identified 1,182 enhancers with significant differences in activity between social and solitary populations, some of which are conserved regulatory regions across species of bees. These results indicate that social trait variation in L. albipes is associated with the fine-tuning of ancient enhancers as well as lineage-specific regulatory changes. Combining enhancer activity with population genetic data revealed variants associated with differences in enhancer activity and identified a subset of differential enhancers with signatures of selection associated with social behavior. Together, these results provide the first empirical map of enhancers in a socially flexible bee and highlight links between cis-regulatory variation and the evolution of social behavior.

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.

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.

Calorie restriction and rapamycin distinctly restore non-canonical ORF translation in the muscles of aging mice

Loss of protein homeostasis is one of the hallmarks of aging. As such, interventions that restore proteostasis should slow down the aging process and improve healthspan. Two of the most broadly used anti-aging interventions that are effective in organisms from yeast to mammals are calorie restriction (CR) and rapamycin (RM) treatment. To identify the regulatory mechanisms by which these interventions improve the protein homeostasis, we carried out ribosome footprinting in the muscle of mice aged under standard conditions, or under long-term treatment with CR or RM. We found that the treatments distinctly impact the non-canonical translation, RM primarily remodeling the translation of upstream open reading frames (uORFs), while CR restores stop codon readthrough and the translation of downstream ORFs. Proteomics analysis revealed the expression of numerous non-canonical ORFs at the protein level. The corresponding peptides may provide entry points for therapies aiming to maintain muscle function and extend health span.

Partial repair causes permanent defects in papillary structure and function after reversal of urinary obstruction

Urinary obstruction causes injury to the renal papilla and leads to defects in the ability to concentrate urine which predisposes to progressive kidney injury. However, the regenerative capacity of the papilla after reversal of obstruction is poorly understood. To address this, we developed a mouse model of reversible urinary obstruction which is characterized by extensive papillary injury, followed by a robust regeneration response and complete histological recovery over a 3- month period. However, these mice have a pronounced defect in urinary concentrating capacity. We now show that this is due to permanent changes in the composition, organization, and transcriptional signatures of epithelial, endothelial, and interstitial cell lineages in the papilla. There are persistent inflammatory responses that are also seen in patients with renal stone disease but are associated with cell-specific adaptive responses to the increasingly hypoxic environment of the papilla after reversal of obstruction. Taken together, our analysis of a new model of reversible urinary obstruction reveals that partial repair leads to permanent changes in the structure and function of all of the major cellular compartments in the papilla that include both shared and distinct responses to different types of renal papillary injury in humans and mice.

Summary

Partial repair after reversal of urinary obstruction leads to permanent changes in structure and function of all major cellular compartments in the renal papilla.

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.

Dear-PSM: A deep learning-based peptide search engine enables full database search for proteomics

Peptide spectrum matching is the process of linking mass spectrometry data with peptide sequences. An experimental spectrum can match thousands of candidate peptides with variable modifications leading to an exponential increase in candidates. Completing the search within a limited time is a key challenge. Traditional searches expedite the process by restricting peptide mass errors and variable modifications, but this limits interpretive capability. To address this challenge, we propose Dear-PSM, a peptide search engine that supports full database searching. Dear-PSM does not restrict peptide mass errors, matching each spectrum to all peptides in the database and increasing the number of variable modifications per peptide from the conventional 3-20. Leveraging inverted index technology, Dear-PSM creates a high-performance index table of experimental spectra and utilizes deep learning algorithms for peptide validation. Through these techniques, Dear-PSM achieves a speed breakthrough 7 times faster than mainstream search engines on a regular desktop computer, with a remarkable 240-fold reduction in memory consumption. Benchmark test results demonstrate that Dear-PSM, in full database search mode, can reproduce over 90% of the results obtained by mainstream search engines when handling complex mass spectrometry data collected from different species using various instruments. Furthermore, it uncovers a substantial number of new peptides and proteins. Dear-PSM has been publicly released on the GitHub repository https://github.com/jianweishuai/Dear-PSM.

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.

Strong Activation of ID1, ID2, and ID3 Genes Is Coupled with the Formation of Vasculogenic Mimicry Phenotype in Melanoma Cells

Gene expression patterns are very sensitive to external influences and are reflected in phenotypic changes. It was previously described that transferring melanoma cells from a plastic surface to Matrigel led to formation of de novo vascular networks-vasculogenic mimicry-that are characteristic to a stemness phenotype in aggressive tumors. Up to now there was no detailed data about the gene signature accompanying this process. Here, we show that this transfer shortly led to extremely strong epigenetic changes in gene expression in the melanoma cells. We observed that on Matrigel numerous genes controlling ribosome biogenesis were upregulated. However, most of the activated genes were inhibitors of the differentiation genes (ID1, ID2, and ID3). At the same time, the genes that control differentiation were downregulated. Both the upregulated and the downregulated genes are simultaneously targeted by different transcription factors shaping sets of co-expressed genes. The specific group of downregulated genes shaping contacts with rDNA genes are also associated with the H3K27me3 mark and with numerous lincRNAs and miRNAs. We conclude that the stemness phenotype of melanoma cells is due to the downregulation of developmental genes and formation of dedifferentiated cells.

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.

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;43:807-820. Published 2023. This article is a U.S. Government work and is in the public domain in the USA.

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.

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.