Using GOstats to test gene lists for GO term association

Developmental Behavioural Plasticity and DNA Methylation Patterns in Response to Predation Stress in Trinidadian Guppies

Early-life experiences can predict the environments experienced later in life, giving individuals an opportunity to develop adaptive behaviour appropriate to a likely future environment. Epigenetic mechanisms such as DNA methylation (DNAm) have been implicated in developmental behavioural plasticity; however, studies investigating this possibility are limited in taxonomic breadth and ecological relevance. We investigated the impact of early-life exposure to predation stress on behaviour and DNAm in the brains of Trinidadian guppies (Poecilia reticulata). We exposed guppies throughout development to either an alarm cue (conspecific skin extract), inducing predation stress, or a control cue (water) for 8 weeks and then raised them to adulthood under identical conditions. Then, we conducted two behavioural assays, an open-field and a grouping test, before performing whole-genome bisulfite sequencing on whole brains. Guppies exposed to the alarm cue during development exhibited increased grouping (shoaling) in adulthood compared to those exposed to the control treatment, but there were no detectable impacts on activity, boldness, or exploratory behaviour. We also identified stable shifts in brain DNAm in response to developmental alarm cue exposure in genes involved in behavioural regulation. Some differentially methylated sites were significantly associated with shoaling propensity in both males and females. Additionally, males and females differed in the magnitude of DNAm responses and the genes impacted, suggesting distinct roles for DNAm between the sexes. This study shows how early-life predation stress can induce behavioural changes in adulthood and that shifts in neural DNAm could be an underlying mechanism responsible for these changes.

Differences in Colour Pattern, Behaviour and Gene Expression in the Brain Suggest Divergent Camouflage Strategies in Sympatric Reef Fish Species

Camouflage is a critical survival strategy that helps to evade predation and increase hunting success. Background matching and disruptive colouration are different camouflage strategies that are subject to different selective pressures and can drive divergence in their associated traits such as colour pattern and behaviour. This study tested whether two closely related reef fish species (Hypoplectrus spp.) with distinct colour patterns exhibit different predator escape responses and differential gene expression in the brain indicative of divergent camouflage strategies. Combining field and laboratory experiments, we show that barred hamlets, characterised by disruptive colouration, are dynamic in their escape responses, while black hamlets, with their darker colouration, had a preference for hiding. The behavioural differences between these species seem to be limited to divergent predator escape responses since other behaviours such as activity or sociability did not differ. Importantly, the observed behavioural differences were accompanied by transcriptomic differences in their brains, particularly in regions associated with the perception of looming threats and less so in the region involved in conditioning. Differential expression in the diencephalon suggests enhanced neuronal plasticity in barred hamlets, which might allow for rapid adjustments in their escape response, while black hamlets exhibited upregulation in genes linked to immune response and oxygen transport in the optic tectum. Overall, our findings suggest that the two species utilise different camouflage strategies, which might contribute to the maintenance of colour pattern differences and thereby influence the speciation and diversification of these closely related sympatric reef fishes.

Combined analysis of somatic mutations and gene expression reveals nuclear speckles-associated enhanced stemness in gingivobuccal carcinoma under DNA damage response

Smokeless tobacco chewing habits in India lead to a high prevalence of Gingivobuccal oral squamous cell carcinoma (OSCC-GB). Cancer stem cells (CSCs) are a sub-population of cancer cells within a tumor with stem-like properties and are believed to contribute to tumor initiation, progression, increased resistance to drug therapy, and promote post-therapeutic cancer relapse. An RNA-sequencing data-based combined analysis of somatic mutations and gene expression was performed to explore the role of CSCs in disease progression using the novel Indian-origin OSCC-GB cell line 'IIOC019' from a patient with tobacco-chewing habit. The identified DNA damage-related known mutational signature (1 bp T/(A) nucleotide insertions and C>T mutations) indicates the impact of smokeless tobacco-related carcinogens in the IIOC019 cell line. The differentially expressed somatic variants, functional impact predictions, and survival analysis reveal the role of DNA damage response (DDR)-related genes in OSCC-GB, with the SON gene as a significant player. The study suggests that the loss-of-function in a somatic variant of the SON gene is linked to nuclear speckles-associated enhanced stemness and increased risk of disease progression in OSCC-GB under DDR conditions. The newly identified CSC-associated somatic variants appear to promote cancer spread, local recurrence, and resistance to chemotherapy or radiotherapy, contributing to the high mortality rates among Indian OSCC-GB patients.

Transcriptome remodeling drives acclimation to iron availability in the marine N2-fixing cyanobacterium Trichodesmium erythraeum IMS101

While enhanced phytoplankton growth as a result of iron (Fe) fertilization has been extensively characterized, our understanding of the underlying mechanisms remains incomplete. Here, we show in a laboratory setup mimicking Fe fertilization in the field that transcriptome remodeling is a primary driver of acclimation to Fe availability in the marine diazotrophic cyanobacterium Trichodesmium erythraeum IMS101. Fe supplementation promoted cell growth, photosynthesis and N2 fixation, and concomitant expression of the photosynthesis and N2 fixation genes. The expression of genes encoding major Fe-binding metalloproteins is tightly linked to cellular carbon and nitrogen metabolism and appears to be controlled by the ferric uptake regulator FurA, which is involved in regulating Fe uptake and homeostasis. This feedback loop is reinforced by substitutive expression of functionally equivalent or competitive genes depending on Fe availability, as well as co-expression of multiple Fe stress inducible isiA genes, an adaptive strategy evolved to elicit the Fe-responsive cascade. The study provides a genome-wide perspective on the acclimation of a prominent marine diazotroph to Fe availability, reveals an upgraded portfolio of indicator genes that can be used to better assess Fe status in the environment, and predicts scenarios of how marine diazotrophs may be affected in the future ocean.IMPORTANCEThe scarcity of trace metal iron (Fe) in global oceans has a great impact on phytoplankton growth. While enhanced primary productivity as a result of Fe fertilization has been extensively characterized, the underlying molecular mechanisms remain poorly understood. By subjecting the model marine diazotroph Trichodesmium erythraeum IMS101 to increasing concentrations of supplemented Fe, we demonstrate in it a comprehensively remodeled transcriptome that drives the mobilization of cellular Fe for coordinated carbon and nitrogen metabolism and reallocation of energy and resources. Our data provide broad genomic insight into marine diazotrophs acclimation to Fe availability, enabling the versatility and flexibility in choice of indicator genes for monitoring Fe status in the environment and having implications on how marine diazotrophs persist into the future ocean.

The gene regulatory mechanisms shaping the heterogeneity of venom production in the Cape coral snake

BACKGROUND

Venoms and their associated glands and delivery structures have evolved numerous times among animals. Within these venom systems, the molecular, cellular, and morphological components interact and co-evolve to generate distinct, venom phenotypes that are increasingly recognized as models for studying adaptive evolution. However, toxins are often unevenly distributed across venom-producing tissues in patterns that are not necessarily adaptive but instead likely result from constraints associated with protein secretion.

RESULTS

We generate a high-quality draft genome of the Cape coral snake (Aspidelaps lubricus) and combine analyses of venom gland single-cell RNA-seq data with spatial venom gland in situ toxin distributions. Our results reveal that while different toxin families are produced by distinct populations of cells, toxin expression is fine-tuned by regulatory modules that result in further specialization of toxin production within each cell population. We also find that the evolution of regulatory elements closely mirrors the evolution of their associated toxin genes, resulting in spatial association of closely related and functionally similar toxins in the venom gland. While this compartmentalization is non-adaptive, the modularity of the underlying regulatory network likely facilitated the repeated evolution of defensive venom in spitting cobras.

CONCLUSIONS

Our results provide new insight into the variability of toxin regulation across snakes, reveal the molecular mechanisms underlying the heterogeneous toxin production in snake venom glands, and provide an example of how constraints can result in non-adaptive character states that appear to be adaptive, which may nevertheless facilitate evolutionary innovation and novelty.

The immune NIK1/RPL10/LIMYB signaling module regulates photosynthesis and translation under biotic and abiotic stresses

Photosynthesis and translation are targets of metabolic control and development in plants, yet how stress signals coordinately regulate these opposing energy-producing and consuming processes remains enigmatic. Here, we unravel a growth control circuit that ties photosynthesis to translational control in response to biotic and abiotic signals. Our findings reveal that the L10-INTERACTING MYB DOMAIN-CONTAINING PROTEIN (LIMYB), a key player of the NUCLEAR SHUTTLE PROTEIN-INTERACTING KINASE 1 (NIK1)/ RIBOSOMAL PROTEIN L10 (RPL10) antiviral signaling pathway, not only downregulates translation genes, but also represses photosynthesis-related genes and photosynthesis itself. LIMYB repressor activity, regulated by phosphorylation, is crucial for the decline in photosynthesis under stress. NIK1 activation by PAMPs or the phosphomimetic NIK1-T474D represses photosynthesis-related genes and photosynthesis in control but not in limyb lines. Furthermore, heat and osmotic stress also activate the NIK1/RPL10/LIMYB signaling circuit in wild type. These stresses induce NIK1 phosphorylation, but not marker gene repression, in limyb, indicating that LIMYB connects NIK1 activation to stress-mediated downregulation of translation- and photosynthesis-related genes. This coordinated repression via the NIK1/RPL10/LIMYB module may help plants adapt to changing environments.

Transcriptomic Landscape of Polypide Development in the Freshwater Bryozoan Cristatella mucedo: From Budding to Degeneration

Colonial invertebrates consist of iterative semi-autonomous modules (usually termed zooids) whose lifespan is significantly shorter than that of the entire colony. Typically, module development begins with budding and ends with degeneration. Most studies on the developmental biology of colonial invertebrates have focused on blastogenesis, whereas the changes occurring throughout the entire zooidal life were examined only for a few tunicates. Here we provide the first description of transcriptomic changes during polypide development in the freshwater bryozoan Cristatella mucedo. For the first time for Bryozoa, we performed bulk RNA sequencing of six polypide stages in C. mucedo (buds, juvenile polypides, three mature stages, and degeneration stage) and generated a high-quality de novo reference transcriptome. Based on these data, we analyzed clusters of differentially expressed genes for enriched pathways and biological processes that may be involved in polypide budding, growth, active functioning, and degradation. Although stem cells have never been described in Bryozoa, our analysis revealed the expression of conservative "stemness" markers in developing buds and juvenile polypides. Our data also indicate that polypide degeneration is a complex regulated process involving autophagy and other types of programmed cell death. We hypothesize that the mTOR signaling pathway plays an important role in regulating the polypide lifespan.

Rapid Neural DNA Methylation Responses to Predation Stress in Trinidadian Guppies

DNA methylation (DNAm) is a well-studied epigenetic mechanism implicated in environmentally induced phenotypes and phenotypic plasticity. However, few studies investigate the timescale of DNAm shifts. Thus, it is uncertain whether DNAm can change on timescales relevant for rapid phenotypic shifts, such as during the expression of short-term behavioural plasticity. DNAm could be especially reactive in the brain, potentially increasing its relevance for behavioural plasticity. Most research investigating neural changes in methylation has been conducted in mammalian systems, on isolated individuals, and using stressors that are less ecologically relevant, reducing their generalisability to other natural systems. We exposed pairs of male and female Trinidadian guppies (Poecilia reticulata) to alarm cue, conspecific skin extract that reliably induces anti-predator behaviour, or a control cue. Whole-genome bisulphite sequencing on whole brains at various time points following cue exposure (0.5, 1, 4, 24, and 72 h) allowed us to uncover the timescale of neural DNAm responses. Males and females both showed rapid shifts in DNAm in as little as 0.5 h. However, males and females differed in the time course of their responses: both sexes showed a peak in the number of loci showing significant responses at 4 h, but males showed an additional peak at 72 h. We suggest that this finding could be due to the differing longer-term plastic responses between the sexes. This study shows that DNAm can be rapidly induced by an ecologically relevant stressor in fish and suggests that DNAm could be involved in short-term behavioural plasticity.

5G-exposed human skin cells do not respond with altered gene expression and methylation profiles

Due to the ever-increasing wirelessly transmitted data, the development of new transmission standards and higher frequencies in the 5G band is required. Despite basic biophysical considerations that argue against health effects, there is public concern about this technology. Because the skin penetration depth at these frequencies is only 1 mm or less, we exposed fibroblasts and keratinocytes to electromagnetic fields up to ten times the permissible limits, for 2 and 48 h in a fully blinded experimental design. Sham-exposed cells served as negative, and UV-exposed cells as positive controls. Differences in gene expression and methylation due to exposure were small and not higher than expected by chance. These data strongly support the assessment that there is no evidence for exposure-induced damage to human skin cells.

Social context prevents heat hormetic effects against mutagens during fish development

Since stress can be transmitted to congeners via social metabolites, it is paramount to understand how the social context of abiotic stress influences aquatic organisms' responses to global changes. Here, we integrated the transcriptomic and phenotypic responses of zebrafish embryos to a UV damage/repair assay following scenarios of heat stress, its social context and their combination. Heat stress preceding UV exposure had a hormetic effect through the cellular stress response and DNA repair, rescuing and/or protecting embryos from UV damage. However, experiencing heat stress within a social context negated this molecular hormetic effect and lowered larval fitness. We discuss the molecular basis of interindividual chemical transmission within animal groups as another layer of complexity to organisms' responses to environmental stressors.

Time-course dual RNA-seq analyses and gene identification during early stages of plant-Phytophthora infestans interactions

Late blight caused by Phytophthora infestans is a major threat to global potato and tomato production. Sustainable management of late blight requires the development of resistant crop cultivars. This process can be facilitated by high-throughput identification of functional genes involved in late blight pathogenesis. In this study, we generated a high-quality transcriptomic time-course dataset focusing on the initial 24 h of contact between P. infestans and 3 solanaceous plant species, tobacco (Nicotiana benthamiana), tomato (Solanum lycopersicum), and potato (Solanum tuberosum). Our results demonstrate species-specific transcriptional regulation in early stages of the infection. Transient silencing of putative RIBOSE-5-PHOSPHATE ISOMERASE and HMG-CoA REDUCTASE genes in N. benthamiana lowered plant resistance against P. infestans. Furthermore, heterologous expression of a putative tomato Golgi-localized nucleosugar transporter-encoding gene exacerbated P. infestans infection of N. benthamiana. In comparison, bioassays using transgenic tomato lines showed that the quantitative disease resistance genes were required but insufficient for late blight resistance; genetic knock-out of the susceptibility gene enhanced resistance. The same RNA-seq dataset was exploited to examine the transcriptional landscape of P. infestans and revealed host-specific gene expression patterns in the pathogen. This temporal transcriptomic diversity, in combination with genomic distribution features, identified the P. infestans IPI-B family GLYCINE-RICH PROTEINs as putative virulence factors that promoted disease severity or induced plant tissue necrosis when transiently expressed in N. benthamiana. These functional genes underline the effectiveness of functional gene-mining through a time-course dual RNA-seq approach and provide insight into the molecular interactions between solanaceous plants and P. infestans.

Independent origins and non-parallel selection signatures of triclabendazole resistance in Fasciola hepatica

Triclabendazole (TCBZ) is the primary treatment for fascioliasis, a global foodborne zoonosis caused by Fasciola hepatica. Widespread resistance to TCBZ (TCBZ-R) in livestock and a rapid rise in resistant human infections are significant concerns. To understand the genetic basis of TCBZ-R, we sequenced the genomes of 99 TCBZ-sensitive (TCBZ-S) and 210 TCBZ-R adult flukes from 146 bovine livers in Cusco, Peru. We identify genomic regions of high differentiation (FST outliers above the 99.9th percentile) that encod genes involved in the EGFR-PI3K-mTOR-S6K pathway and microtubule function. Transcript expression differences are observed in microtubule-related genes between TCBZ-S and -R flukes, both without drug treatment and in response to treatment. Using only 30 SNPs, it is possible to differentiate between TCBZ-S and -R parasites with ≥75% accuracy. Our outlier loci are distinct from the previously reported TCBZ-R-associated QTLs in the UK, suggesting an independent evolution of resistance alleles. Effective genetics-based TCBZ-R surveillance must consider the heterogeneity of loci under selection across diverse geographical populations.

Survey of gene, lncRNA and transposon transcription patterns in four mouse organs highlights shared and organ-specific sex-biased regulation

BACKGROUND

Sex-biased gene regulation is the basis of sexual dimorphism in phenotypes and has been studied across different cell types and different developmental stages. However, sex-biased expression of transposable elements (TEs), which represent nearly half of the mammalian genome and have the potential of influencing genome integrity and regulation, remains underexplored.

RESULTS

We report a survey of gene, lncRNA, and TE expression in four organs from mice with different combinations of gonadal and genetic sex. The data show remarkable variability among organs with respect to the impact of gonadal sex on transcription with the strongest effects observed in the liver. In contrast, the X-chromosome dosage alone had a modest influence on sex-biased transcription across organs, albeit interaction between X-dosage and gonadal sex cannot be ruled out. The presence of the Y-chromosome influences TE, but not gene or lncRNA, expression in the liver. Notably, 90% of sex-biased TEs (sDETEs) reside in clusters. Moreover, 54% of these clusters overlap or reside less than 100 kb from sex-biased genes or lncRNAs, share the same sex bias, and also have higher expression levels than sDETE clusters that do not co-localize with other types of sex-biased transcripts. We test the heterochromatic sink hypothesis that predicts higher expression of TEs in XX individuals finding no evidence to support it.

CONCLUSIONS

Our data show that sex-biased expression of TEs varies among organs with the highest numbers of sDETEs found in the liver following trends observed for genes and lncRNAs. It is enhanced by proximity to other types of sex-biased transcripts.

Differential Gene Expression and Unbalanced Homeolog Expression Bias in 4 Million-Year-Old Allopolyploids of Nicotiana Section Repandae

Allopolyploidy, a phenomenon prevalent in angiosperms involving hybridization and whole-genome duplication, results in species with multiple subgenomes, altering genome structure and gene expression, leading to novel phenotypes. Allopolyploids often experience unbalanced homeolog expression bias, the preferential expression of homeologs from one of the two progenitor genomes. To explore the consequences of allopolyploidy and unbalanced homeolog expression bias, we investigate global gene expression and the fate of homeologs in Nicotiana (Solanaceae). We focus on Nicotiana section Repandae, including three allotetraploid species, Nicotiana nudicaulis, N. repanda, and N. stocktonii, derived from diploid progenitors N. sylvestris and N. obtusifolia ∼4.3 Ma. We identify genes with differential expression and investigate expression of candidate genes for flower size variation. Our results show expression differences with the allopolyploids intermediate between the two progenitor species, with a slight bias toward N. obtusifolia. Moreover, we demonstrate unbalanced homeolog expression bias toward the N. obtusifolia subgenome across developmental stages in the allopolyploids, with a stronger bias in N. nudicaulis. In contrast, unbalanced homeolog expression bias shifts toward N. sylvestris for flower size genes in N. nudicaulis, showing that genes involved in particular phenotypes can display different patterns of unbalanced homeolog expression than the overall transcriptome. We also see differential expression of several known flower size genes across corolla developmental stages. Our results highlight the role of unbalanced homeolog expression bias in shaping the evolutionary trajectory of Nicotiana species and provide a foundation for future research into the ecological and evolutionary implications of allopolyploidy in flowering plants.

CRISPR-Cas9-mediated deletions of FvMYB46 in Fragaria vesca reveal its role in regulation of fruit set and phenylpropanoid biosynthesis

The phenylpropanoid pathway, regulated by transcription factors of the MYB family, produces secondary metabolites that play important roles in fertilization and early phase of fruit development. The MYB46 transcription factor is a key regulator of secondary cell wall structure, lignin and flavonoid biosynthesis in many plants, but little is known about its activity in flowers and berries in F. vesca. For functional analysis of FvMYB46, we designed a CRISPR-Cas9 construct with an endogenous F. vesca-specific U6 promoter for efficient and specific expression of two gRNAs targeting the first exon of FvMYB46. This generated mutants with an in-frame 81-bp deletion of the first conserved MYB domain or an out-of-frame 82-bp deletion potentially knocking out gene function. In both types of mutant plants, pollen germination and fruit set were significantly reduced compared to wild type. Transcriptomic analysis of flowers revealed that FvMYB46 positively regulates the expression of genes involved in processes like xylan biosynthesis and metabolism, homeostasis of reactive oxygen species (ROS) and the phenylpropanoid pathway, including secondary cell wall biosynthesis and flavonoid biosynthesis. Genes regulating carbohydrate metabolism and signalling were also deregulated, suggesting that FvMYB46 might regulate the crosstalk between carbohydrate metabolism and phenylpropanoid biosynthesis. In the FvMYB46-mutant flowers, the flavanol and flavan-3-ol contents, especially epicatechin, quercetin-glucoside and kaempferol-3-coumaroylhexoside, were reduced, and we observed a local reduction in the lignin content in the anthers. Together, these results suggest that FvMYB46 controls fertility and efficient fruit set by regulating the cell wall structure, flavonoid biosynthesis, carbohydrate metabolism, and sugar and ROS signalling in flowers and early fruit development in F. vesca.

Evolutionary innovations in germline biology of placental mammals identified by transcriptomics of first-wave spermatogenesis in opossum

Mammalian spermatogenesis is a highly stereotyped and conserved developmental process that is essential for fitness. At the same time, gene expression in spermatogenic cells is rapidly evolving. This combination of features has been suggested to drive rapid fixation of new gene expression patterns. Using a high-resolution dataset comprising bulk and single-cell data from juvenile and adult testes of the opossum Monodelphis domestica, a model marsupial, we define the developmental timing of the spermatogenic first wave in opossum and delineate conserved and divergent gene expression programs across the placental-marsupial split by comparison to equivalent data from mouse, a model placental mammal. Epigenomic data confirmed divergent regulation at the level of transcription, and comparison to data from four additional amniote species identified hundreds of genes with evidence of rapid fixation of expression. This gene set encompasses known and previously undescribed regulators of spermatogenic development.

Blue light restores functional circadian clocks in eyeless cave spiders

Evolution in profound darkness often leads to predictable, convergent traits, such as the loss of vision. Yet, the consequences of such repeated evolutionary experiments remain obscure, especially regarding fundamental regulatory behaviors like circadian rhythms. We studied circadian clocks of blind cave spiders and their sighted relatives. In the field, cave spiders exhibit low per expression and maintain constant activity levels. Curiously, their clocks are not permanently lost; exposure to monochromatic blue light restores both circadian gene expression and behavioral rhythms. Conversely, blocking blue light in sighted relatives induces an arrhythmic "cave phenotype." Our RNA interference experiments suggest that clock genes regulate the rhythmicity of the huddle response, establishing a link between circadian gene networks and this behavioral rhythm. We demonstrate that circadian regulation is readily toggled and may play a latent role, even in constant darkness. Overall, our study expands understanding of circadian clock variations and paves the way for future research on the maintenance of silent phenotypes.

Non-Invasive Quality Control of Organoid Cultures Using Mesofluidic CSTR Bioreactors and High-Content Imaging

Human brain organoids produce anatomically relevant cellular structures and recapitulate key aspects of in vivo brain function, which holds great potential to model neurological diseases and screen therapeutics. However, the long growth time of 3D systems complicates the culturing of brain organoids and results in heterogeneity across samples hampering their applications. We developed an integrated platform to enable robust and long-term culturing of 3D brain organoids. We designed a mesofluidic bioreactor device based on a reaction-diffusion scaling theory, which achieves robust media exchange for sufficient nutrient delivery in long-term culture. We integrated this device with longitudinal tracking and machine learning-based classification tools to enable non-invasive quality control of live organoids. This integrated platform allows for sample pre-selection for downstream molecular analysis. Transcriptome analyses of organoids revealed that our mesofluidic bioreactor promoted organoid development while reducing cell death. Our platform thus offers a generalizable tool to establish reproducible culture standards for 3D cellular systems for a variety of applications beyond brain organoids.

Glycogenin is dispensable for normal liver glycogen metabolism and body glucose homeostasis

Glycogen is a glucose-storage polysaccharide molecule present in animals, fungi and bacteria. The enzyme glycogenin can self-glycosylate, forming an oligosaccharide chain that primes glycogen synthesis. This priming role of glycogenin was first believed to be essential for glycogen synthesis, but glycogen was then found in the skeletal muscle, heart, liver and brain of glycogenin-knockout mice (Gyg KO), thereby showing that glycogen can be synthesized without glycogenin. Within the liver, glycogen is present in the form of individual glycogen particles, called β particles, and larger composite aggregates of linked β particles, called α particles. Previous studies suggested that liver glycogenin plays a role in linking β particles into α particles and thus participating in glucose homeostasis, which implies that α particles would be absent in Gyg KO mice liver. Here we test this through targeted characterization of glycogen structure and through proteomic and metabolic studies on Gyg KO mice. The results show that, contrary to what had been believed, glycogenin is not necessary for normal liver-glycogen metabolism.

IL-1 signaling enrichment in inflammatory skin disease loci with higher-risk allele frequencies in African ancestry

Inflammatory skin diseases (ISDs) exhibit varying prevalence across different ancestry background and geographical regions. Genetic research for complex ISDs has predominantly centered on European Ancestry (EurA) populations and genetic effects on immune cell responses but generally failed to consider contributions from other cell types in skin. Here, we utilized 273 genetic signals from seven different ISDs: acne, alopecia areata (AA), atopic dermatitis (AD), psoriasis, systemic lupus erythematosus (SLE), systemic sclerosis (SSc), and vitiligo, to demonstrate enriched IL1 signaling in keratinocytes, particularly in signals with higher risk allele frequencies in the African ancestry. Using a combination of ATAC-seq, Bru-seq, and promoter capture Hi-C, we revealed potential regulatory mechanisms of the acne locus on chromosome 2q13. We further demonstrated differential responses in keratinocytes upon IL1β stimulation, including the pro-inflammatory mediators CCL5, IL36G, and CXCL8. Taken together, our findings highlight IL1 signaling in epidermal keratinocytes as a contributor to ancestry-related differences in ISDs.

Inferring disease progression stages in single-cell transcriptomics using a weakly supervised deep learning approach

Application of single-cell/nucleus genomic sequencing to patient-derived tissues offers potential solutions to delineate disease mechanisms in humans. However, individual cells in patient-derived tissues are in different pathological stages, and hence, such cellular variability impedes subsequent differential gene expression analyses. To overcome such a heterogeneity issue, we present a novel deep learning approach, scIDST, that infers disease progression levels of individual cells with weak supervision framework. The disease progression-inferred cells display significant differential expression of disease-relevant genes, which cannot be detected by comparative analysis between patients and healthy donors. In addition, we demonstrate that pretrained models by scIDST are applicable to multiple independent data resources and are advantageous to infer cells related to certain disease risks and comorbidities. Taken together, scIDST offers a new strategy of single-cell sequencing analysis to identify bona fide disease-associated molecular features.

Life-history adaptation under climate warming magnifies the agricultural footprint of a cosmopolitan insect pest

Climate change is affecting population growth rates of ectothermic pests with potentially dire consequences for agriculture and global food security. However, current projection models of pest impact typically overlook the potential for rapid genetic adaptation, making current forecasts uncertain. Here, we predict how climate change adaptation in life-history traits of insect pests affects their growth rates and impact on agricultural yields by unifying thermodynamics with classic theory on resource acquisition and allocation trade-offs between foraging, reproduction, and maintenance. Our model predicts that warming temperatures will favour resource allocation towards maintenance coupled with increased resource acquisition through larval foraging, and the evolution of this life-history strategy results in both increased population growth rates and per capita host consumption, causing a double-blow on agricultural yields. We find support for these predictions by studying thermal adaptation in life-history traits and gene expression in the wide-spread insect pest, Callosobruchus maculatus; with 5 years of evolution under experimental warming causing an almost two-fold increase in its predicted agricultural footprint. These results show that pest adaptation can offset current projections of agricultural impact and emphasize the need for integrating a mechanistic understanding of life-history evolution into forecasts of pest impact under climate change.

Clinical-grade extracellular vesicles derived from umbilical cord mesenchymal stromal cells: preclinical development and first-in-human intra-articular validation as therapeutics for knee osteoarthritis

Osteoarthritis (OA) is a joint disease characterized by articular cartilage degradation. Persistent low-grade inflammation defines OA pathogenesis, with crucial involvement of pro-inflammatory M1-like macrophages. While mesenchymal stromal cells (MSC) and their small extracellular vesicles (sEV) hold promise for OA treatment, achieving consistent clinical-grade sEV products remains a significant challenge. This study aims to develop fully characterized, reproducible, clinical-grade batches of sEV derived from umbilical cord (UC)-MSC for the treatment of OA while assessing its efficacy and safety. Initially, a standardized, research-grade manufacturing protocol was established to ensure consistent sEV production. UC-MSC-sEV characterization under non-cGMP conditions showed consistent miRNA and protein profiles, suggesting their potential for standardized manufacturing. In vitro studies evaluated the efficacy, safety, and potency of sEV; animal studies confirmed their effectiveness and safety. In vitro, UC-MSC-sEV polarized macrophages to an anti-inflammatory M2b-like phenotype, through STAT1 modulation, indicating their potential to create an anti-inflammatory environment in the affected joints. In silico studies confirmed sEV's immunosuppressive signature through miRNA and proteome analysis. In an OA mouse model, sEV injected intra-articularly (IA) induced hyaline cartilage regeneration, validated by histological and μCT analyses. The unique detection of sEV signals within the knee joint over time highlights its safety profile by confirming the retention of sEV in the joint. The product development of UC-MSC-sEV involved refining, standardizing, and validating processes in compliance with GMP standards. The initial assessment of the safety of the clinical-grade product via IA administration in a first-in-human study showed no adverse effects after a 12 month follow-up period. These results support the progress of this sEV-based therapy in an early-phase clinical trial, the details of which are presented and discussed in this work. This study provides data on using UC-MSC-sEV as local therapy for OA, highlighting their regenerative and anti-inflammatory properties and safety in preclinical and a proof-of-principle clinical application.

The B Chromosome of Pseudococcus viburni: A Selfish Chromosome that Exploits Whole-Genome Meiotic Drive

Meiosis is generally a fair process: each chromosome has a 50% chance of being included into each gamete. However, meiosis can become aberrant with some chromosomes having a higher chance of making it into gametes than others. Yet, why and how such systems evolve remains unclear. Here, we study the unusual reproductive genetics of mealybugs, where only maternal-origin chromosomes are included in gametes during male meiosis, while paternal chromosomes are eliminated. One species-Pseudococcus viburni-has a segregating B chromosome that drives by escaping paternal genome elimination. We present whole genome and gene expression data from lines with and without B chromosomes. We identify B-linked sequences including 204 protein-coding genes and a satellite repeat that makes up a significant proportion of the chromosome. The few paralogs between the B and the core genome are distributed throughout the genome, arguing against a simple, or at least recent, chromosomal duplication of one of the autosomes to create the B. We do, however, find one 373 kb region containing 146 genes that appears to be a recent translocation. Finally, we show that while many B-linked genes are expressed during meiosis, most of these are encoded on the recently translocated region. Only a small number of B-exclusive genes are expressed during meiosis. Of these, only one was overexpressed during male meiosis, which is when the drive occurs: an acetyltransferase involved in H3K56Ac, which has a putative role in meiosis and is, therefore, a promising candidate for further studies.

LncRNA XLOC-040580 targeted by TPRA1 coordinate zygotic genome activation during porcine embryonic development

Long noncoding RNAs (lncRNAs) are crucial in porcine preimplantation embryonic development, yet their regulatory role during zygote genome activation (ZGA) is poorly understood. We analyzed transcriptome data from porcine fetal fibroblasts (PEF), induced pluripotent stem cells (iPS), and preimplantation embryos, identifying ZGA-specific lncRNAs like XLOC-040580, and further predicted its potentially interacting genes TPRA1 and BCL2L1 via co-expression network. XLOC-040580 was knocked down by siRNA microinjection and the expression of ZGA-related genes was detected by qRT-PCR. After microinjecting siRNA targeting TPRA1 and BCL2L1 at the one-cell stage, we counted the blastocyst development rate. The blastocyst development rate was consistent with the results from si-XLOC-040580 after si-TPRA1. Through dual-luciferase reporter assays, we found that XLOC-040580 was a downstream target of TPRA1. To further elucidate the mechanism of XLOC-040580, Single-cell mRNA sequencing after XLOC-040580 knockdown revealed its regulatory network involved in embryonic developmental defects. Transcriptome analysis revealed that XLOC-040580 was specifically expressed during zygote activation. Knockdown of XLOC-040580 decreased the blastocyst development rate and reduced both the total blastocyst cell number and TE cell number. TPRA1 and BCL2L1 were specifically co-expressed with XLOC-040580 during ZGA stage, and TPRA1 could interact with the promoter region of XLOC-040580 and regulate its expression. Knockdown of TPRA1 or XLOC-040580 blocked porcine embryonic development by affecting the expression of ZGA-related genes. We found and validated that lncRNA XLOC-040580 played a key role in the ZGA process, which was regulated by TPRA1. These results implied that the functional axis of TPRA1-XLOC-040580-downstream genes involved in ZGA-related functions also coordinated early embryonic development in porcine.

Proteomic characterization and comparison of the infective and adult life stage secretomes from Necator americanus and Ancylostoma ceylanicum

More than 470 million people globally are infected with the hookworms Ancylostoma ceylanicum and Necator americanus, resulting in an annual loss of 2.1 to 4 million disability-adjusted-life-years. Current infection management approaches are limited by modest drug efficacy, the costs associated with frequent mass drug administration campaigns, and the risk of reinfection and burgeoning drug resistance. Subunit vaccines based on proteins excreted and secreted (ES) by hookworms that reduce worm numbers and associated disease burden are a promising management strategy to overcome these limitations. However, studies on the ES proteomes of hookworms have mainly described proteins from the adult life stage which may preclude the opportunity to target the infective larva. Here, we employed high resolution mass spectrometry to identify 103 and 57 ES proteins from the infective third larvae stage (L3) as well as 106 and 512 ES proteins from the adult N. americanus and A. ceylanicum respectively. Comparisons between these developmental stages identified 91 and 41 proteins uniquely expressed in the L3 ES products of N. americanus and A. ceylanicum, respectively. We characterized these proteins based on functional annotation, KEGG pathway analysis, InterProScan signature and gene ontology. We also performed reciprocal BLAST analysis to identify orthologs across species for both the L3 and adult stages and identified five orthologous proteins in both life stages and 15 proteins that could be detected only in the L3 stage of both species. Last, we performed a three-way reciprocal BLAST on the L3 proteomes from both hookworm species together with a previously reported L3 proteome from the rodent hookworm Nippostrongylus brasiliensis, and identified eight L3 proteins that could be readily deployed for testing using well established rodent models. This novel characterization of L3 proteins and taxonomic conservation across hookworm species provides a raft of potential candidates for vaccine discovery for prevention of hookworm infection and disease.

Symbiotic T6SS affects horizontal transmission of Paraburkholderia bonniea among Dictyostelium discoideum amoeba hosts

Three species of Paraburkholderia are able to form facultative symbiotic relationships with the amoeba, Dictyostelium discoideum. These symbiotic Paraburkholderia share a type VI secretion system (T6SS) that is absent in other close relatives. We tested the phenotypic and transcriptional effect of tssH ATPase gene disruption in P. bonniea on its symbiosis with D. discoideum. We hypothesized that the ∆tssH mutant would have a significantly reduced ability to affect host fitness or transmit itself from host to host. We found that the T6SS does not directly affect host fitness. Instead, wildtype P. bonniea had significantly higher rates of horizontal transmission compared to ∆tssH. In addition, we observed significant differences in the range of infection prevalence achieved by wildtype vs. ∆tssH symbionts over multiple host social stages in the absence of opportunities for environmental symbiont acquisition. Successful symbiont transmission significantly contributes to sustained symbiotic association. Therefore, the shared T6SS appears necessary for a long-term evolutionary relationship between D. discoideum and its Paraburkholderia symbionts. The lack of difference in host fitness outcomes was confirmed by indistinguishable host gene expression patterns between hosts infected by wildtype or ∆tssH P. bonniea in an RNA-seq time series. These data also provided insight into how Paraburkholderia symbionts may evade phagocytosis by its amoeba host. Most significantly, cellular oxidant detoxification and lysosomal hydrolase delivery appear to be subject to the push and pull of host-symbiont crosstalk.

Human-specific epigenomic states in spermatogenesis

Infertility is becoming increasingly common, affecting one in six people globally. Half of these cases can be attributed to male factors, many driven by abnormalities in the process of sperm development. Emerging evidence from genome-wide association studies, genetic screening of patient cohorts, and animal models highlights an important genetic contribution to spermatogenic defects, but comprehensive identification and characterization of the genes critical for male fertility remain lacking. High divergence of gene regulation in spermatogenic cells across species poses challenges for delineating the genetic pathways required for human spermatogenesis using common model organisms. In this study, we leveraged post-translational histone modification and gene transcription data for 15,491 genes in four mammalian species (human, rhesus macaque, mouse, and opossum), to identify human-specific patterns of gene regulation during spermatogenesis. We combined H3K27me3 ChIP-seq, H3K4me3 ChIP-seq, and RNA-seq data to define epigenetic states for each gene at two stages of spermatogenesis, pachytene spermatocytes and round spermatids, in each species. We identified 239 genes that are uniquely active, poised, or dynamically regulated in human spermatogenic cells distinct from the other three species. While some of these genes have been implicated in reproductive functions, many more have not yet been associated with human infertility and may be candidates for further molecular and epidemiologic studies. Our analysis offers an example of the opportunities provided by evolutionary and epigenomic data for broadly screening candidate genes implicated in reproduction, which might lead to discoveries of novel genetic targets for diagnosis and management of male infertility and male contraception.

Genomic insights into variation in thermotolerance between hybridizing swordtail fishes

Understanding how organisms adapt to changing environments is a core focus of research in evolutionary biology. One common mechanism is adaptive introgression, which has received increasing attention as a potential route to rapid adaptation in populations struggling in the face of ecological change, particularly global climate change. However, hybridization can also result in deleterious genetic interactions that may limit the benefits of adaptive introgression. Here, we used a combination of genome-wide quantitative trait locus mapping and differential gene expression analyses between the swordtail fish species Xiphophorus malinche and X. birchmanni to study the consequences of hybridization on thermotolerance. While these two species are adapted to different thermal environments, we document a complicated architecture of thermotolerance in hybrids. We identify a region of the genome that contributes to reduced thermotolerance in individuals heterozygous for X. malinche and X. birchmanni ancestry, as well as widespread misexpression in hybrids of genes that respond to thermal stress in the parental species, particularly in the circadian clock pathway. We also show that a previously mapped hybrid incompatibility between X. malinche and X. birchmanni contributes to reduced thermotolerance in hybrids. Together, our results highlight the challenges of understanding the impact of hybridization on complex ecological traits and its potential impact on adaptive introgression.

Correlating gene expression levels with transcription factor binding sites facilitates identification of key transcription factors from transcriptome data

Identification of key transcription factors from transcriptome data by correlating gene expression levels with transcription factor binding sites is important for transcriptome data analysis. In a typical scenario, we always set a threshold to filter the top ranked differentially expressed genes and top ranked transcription factor binding sites. However, correlation analysis of filtered data can often result in spurious correlations. In this study, we tested four methods for creating the gene expression inputs (ranked gene list) in the correlation analysis: star coordinate map transformation (START), expression differential score (ED), preferential expression measure (PEM), and the specificity measure (SPM). Then, Kendall's tau correlation statistical algorithms implementing the standard (STD), LINEAR, MIX-LINEAR, DENSITY-CURVE, and MIXED-DENSITY-CURVE weighting methods were used to identify key transcription factors. ED was identified as the optimal method for creating a ranked gene list from filtered expression data, which can address the "unable to detect negative correlation" fallacy presented by other methods. The MIXED-DENSITY-CURVE was the most sensitive for identifying transcription factors from the gene set and list in which only the top proportion was correlated. Ultimately, 644 transcription factor candidates were identified from the transcriptome data of 1,206 cell lines, six of which were validated by wet lab experiments. The Jinzer and Flaver software implementing these methods can be obtained from http://www.thua45/cn/flaver under a free academic license.

Using clusterProfiler to characterize multiomics data

With the advent of multiomics, software capable of multidimensional enrichment analysis has become increasingly crucial for uncovering gene set variations in biological processes and disease pathways. This is essential for elucidating disease mechanisms and identifying potential therapeutic targets. clusterProfiler stands out for its comprehensive utilization of databases and advanced visualization features. Importantly, clusterProfiler supports various biological knowledge, including Gene Ontology and Kyoto Encyclopedia of Genes and Genomes, through performing over-representation and gene set enrichment analyses. A key feature is that clusterProfiler allows users to choose from various graphical outputs to visualize results, enhancing interpretability. This protocol describes innovative ways in which clusterProfiler has been used for integrating metabolomics and metagenomics analyses, identifying and characterizing transcription factors under stress conditions, and annotating cells in single-cell studies. In all cases, the computational steps can be completed within ~2 min. clusterProfiler is released through the Bioconductor project and can be accessed via https://bioconductor.org/packages/clusterProfiler/ .

Population Genomics of Adaptive Radiations: Exceptionally High Levels of Genetic Diversity and Recombination in an Endemic Spider From the Canary Islands

The spider genus Dysdera has undergone a remarkable diversification in the oceanic archipelago of the Canary Islands, with ~60 endemic species having originated during the 20 million years since the origin of the archipelago. This evolutionary radiation has been accompanied by substantial dietary shifts, often characterised by phenotypic modifications encompassing morphological, metabolic and behavioural changes. Hence, these endemic spiders represent an excellent model for understanding the evolutionary drivers and to pinpoint the genomic determinants underlying adaptive radiations. Recently, we achieved the first chromosome-level genome assembly of one of the endemic species, D. silvatica, providing a high-quality reference sequence for evolutionary genomics studies. Here, we conducted a low coverage-based resequencing study of a natural population of D. silvatica from La Gomera island. Taking advantage of the new high-quality genome, we characterised genome-wide levels of nucleotide polymorphism, divergence and linkage disequilibrium, and inferred the demographic history of this population. We also performed comprehensive genome-wide scans for recent positive selection. Our findings uncovered exceptionally high levels of nucleotide diversity and recombination in this geographically restricted endemic species, indicative of large historical effective population sizes. We also identified several candidate genomic regions that are potentially under positive selection, highlighting relevant biological processes, such as vision and nitrogen extraction as potential adaptation targets. These processes may ultimately drive species diversification in this genus. This pioneering study of spiders that are endemic to an oceanic archipelago lays the groundwork for broader population genomics analyses aimed at understanding the genetic mechanisms driving adaptive radiation in island ecosystems.

Identification of the Regulatory Elements and Protein Substrates of Lysine Acetoacetylation

Short chain fatty acylations establish connections between cell metabolism and regulatory pathways. Lysine acetoacetylation (Kacac) was recently identified as a new histone mark. However, regulatory elements, substrate proteins, and epigenetic functions of Kacac remain unknown, hindering further in-depth understanding of acetoacetate modulated (patho)physiological processes. Here, we created a chemo-immunological approach for reliable detection of Kacac, and demonstrated that acetoacetate serves as the primary precursor for histone Kacac. We report the enzymatic addition of the Kacac mark by the acyltransferases GCN5, p300, and PCAF, and its removal by deacetylase HDAC3. Furthermore, we establish acetoacetyl-CoA synthetase (AACS) as a key regulator of cellular Kacac levels. A comprehensive proteomic analysis has identified 139 Kacac sites on 85 human proteins. Bioinformatics analysis of Kacac substrates and RNA-seq data reveal the broad impacts of Kacac on multifaceted cellular processes. These findings unveil pivotal regulatory mechanisms for the acetoacetate-mediated Kacac pathway, opening a new avenue for further investigation into ketone body functions in various pathophysiological states.

mulea: An R package for enrichment analysis using multiple ontologies and empirical false discovery rate

Traditional gene set enrichment analyses are typically limited to a few ontologies and do not account for the interdependence of gene sets or terms, resulting in overcorrected p-values. To address these challenges, we introduce mulea, an R package offering comprehensive overrepresentation and functional enrichment analysis. mulea employs a progressive empirical false discovery rate (eFDR) method, specifically designed for interconnected biological data, to accurately identify significant terms within diverse ontologies. mulea expands beyond traditional tools by incorporating a wide range of ontologies, encompassing Gene Ontology, pathways, regulatory elements, genomic locations, and protein domains. This flexibility enables researchers to tailor enrichment analysis to their specific questions, such as identifying enriched transcriptional regulators in gene expression data or overrepresented protein domains in protein sets. To facilitate seamless analysis, mulea provides gene sets (in standardised GMT format) for 27 model organisms, covering 22 ontology types from 16 databases and various identifiers resulting in almost 900 files. Additionally, the muleaData ExperimentData Bioconductor package simplifies access to these pre-defined ontologies. Finally, mulea's architecture allows for easy integration of user-defined ontologies, or GMT files from external sources (e.g., MSigDB or Enrichr), expanding its applicability across diverse research areas. mulea is distributed as a CRAN R package downloadable from https://cran.r-project.org/web/packages/mulea/ and https://github.com/ELTEbioinformatics/mulea . It offers researchers a powerful and flexible toolkit for functional enrichment analysis, addressing limitations of traditional tools with its progressive eFDR and by supporting a variety of ontologies. Overall, mulea fosters the exploration of diverse biological questions across various model organisms.

Sirtuin 1/sirtuin 3 are robust lysine delactylases and sirtuin 1-mediated delactylation regulates glycolysis

Lysine lactylation (Kla), an epigenetic mark triggered by lactate during glycolysis, including the Warburg effect, bridges metabolism and gene regulation. Enzymes such as p300 and HDAC1/3 have been pivotal in deciphering the regulatory dynamics of Kla, though questions about additional regulatory enzymes, their specific Kla substrates, and the underlying functional mechanisms persist. Here, we identify SIRT1 and SIRT3 as key "erasers" of Kla, shedding light on their selective regulation of both histone and non-histone proteins. Proteomic analysis in SIRT1/SIRT3 knockout HepG2 cells reveals distinct substrate specificities toward Kla, highlighting their unique roles in cellular signaling. Notably, we highlight the role of specific Kla modifications, such as those on the M2 splice isoform of pyruvate kinase (PKM2), in modulating metabolic pathways and cell proliferation, thereby expanding Kla's recognized functions beyond epigenetics. Therefore, this study deepens our understanding of Kla's functional mechanisms and broadens its biological significance.

ISWI chromatin remodeler SMARCA5 is essential for meiotic gene expression and male fertility in mammals

Regulation of the transcriptome to promote meiosis is important for sperm development and fertility. However, how chromatin remodeling directs the transcriptome during meiosis in male germ cells is largely unknown. Here, we demonstrate that the ISWI family ATP-dependent chromatin remodeling factor SMARCA5 (SNF2H) plays a critical role in regulating meiotic prophase progression during spermatogenesis. Males with germ cell-specific depletion of SMARCA5 are infertile and unable to form sperm. Loss of Smarca5 results in failure of meiotic progression with abnormal spermatocytes beginning at the pachytene stage and an aberrant global increase in chromatin accessibility, especially at genes important for meiotic prophase.

Novel Sex-Specific Genes and Diverse Interspecific Expression in the Antennal Transcriptomes of Ithomiine Butterflies

The olfactory sense is crucial for organisms, facilitating environmental recognition and interindividual communication. Ithomiini butterflies exemplify this importance not only because they rely strongly on olfactory cues for both inter- and intra-sexual behaviors, but also because they show convergent evolution of specialized structures within the antennal lobe, called macroglomerular complexes (MGCs). These structures, widely absent in butterflies, are present in moths where they enable heightened sensitivity to, and integration of, information from various types of pheromones. In this study, we investigate chemosensory evolution across six Ithomiini species and identify possible links between expression profiles and neuroanatomical. To enable this, we sequenced four new high-quality genome assemblies and six sex-specific antennal transcriptomes for three of these species with different MGC morphologies. With extensive genomic analyses, we found that the expression of antennal transcriptomes across species exhibit profound divergence, and identified highly expressed ORs, which we hypothesize may be associated to MGCs, as highly expressed ORs are absent in Methona, an Ithomiini lineage which also lacks MGCs. More broadly, we show how antennal sexual dimorphism is prevalent in both chemosensory genes and non-chemosensory genes, with possible relevance for behavior. As an example, we show how lipid-related genes exhibit consistent sexual dimorphism, potentially linked to lipid transport or host selection. In this study, we investigate the antennal chemosensory adaptations, suggesting a link between genetic diversity, ecological specialization, and sensory perception with the convergent evolution of MCGs. Insights into chemosensory gene evolution, expression patterns, and potential functional implications enhance our knowledge of sensory adaptations and sexual dimorphisms in butterflies, laying the foundation for future investigations into the genetic drivers of insect behavior, adaptation, and speciation.

GAN-WGCNA: Calculating gene modules to identify key intermediate regulators in cocaine addiction

Understanding time-series interplay of genes is essential for diagnosis and treatment of disease. Spatio-temporally enriched NGS data contain important underlying regulatory mechanisms of biological processes. Generative adversarial networks (GANs) have been used to augment biological data to describe hidden intermediate time-series gene expression profiles during specific biological processes. Developing a pipeline that uses augmented time-series gene expression profiles is needed to provide an unbiased systemic-level map of biological processes and test for the statistical significance of the generated dataset, leading to the discovery of hidden intermediate regulators. Two analytical methods, GAN-WGCNA (weighted gene co-expression network analysis) and rDEG (rescued differentially expressed gene), interpreted spatiotemporal information and screened intermediate genes during cocaine addiction. GAN-WGCNA enables correlation calculations between phenotype and gene expression profiles and visualizes time-series gene module interplay. We analyzed a transcriptome dataset of two weeks of cocaine self-administration in C57BL/6J mice. Utilizing GAN-WGCNA, two genes (Alcam and Celf4) were selected as missed intermediate significant genes that showed high correlation with addiction behavior. Their correlation with addictive behavior was observed to be notably significant in aspect of statistics, and their expression and co-regulation were comprehensively mapped in terms of time, brain region, and biological process.

An anciently diverged family of RNA binding proteins maintain correct splicing of a class of ultra-long exons through cryptic splice site repression

Previously, we showed that the germ cell-specific nuclear protein RBMXL2 represses cryptic splicing patterns during meiosis and is required for male fertility (Ehrmann et al., 2019). Here, we show that in somatic cells the similar yet ubiquitously expressed RBMX protein has similar functions. RBMX regulates a distinct class of exons that exceed the median human exon size. RBMX protein-RNA interactions are enriched within ultra-long exons, particularly within genes involved in genome stability, and repress the selection of cryptic splice sites that would compromise gene function. The RBMX gene is silenced during male meiosis due to sex chromosome inactivation. To test whether RBMXL2 might replace the function of RBMX during meiosis we induced expression of RBMXL2 and the more distantly related RBMY protein in somatic cells, finding each could rescue aberrant patterns of RNA processing caused by RBMX depletion. The C-terminal disordered domain of RBMXL2 is sufficient to rescue proper splicing control after RBMX depletion. Our data indicate that RBMX and RBMXL2 have parallel roles in somatic tissues and the germline that must have been conserved for at least 200 million years of mammalian evolution. We propose RBMX family proteins are particularly important for the splicing inclusion of some ultra-long exons with increased intrinsic susceptibility to cryptic splice site selection.

The eukaryotic translation initiation factor eIF4E regulates flowering and circadian rhythm in Arabidopsis

Translation initiation is a critical, rate-limiting step in protein synthesis. The eukaryotic translation initiation factor 4E (eIF4E) plays an essential role in this process. However, the mechanisms by which eIF4E-dependent translation initiation regulates plant growth and development remain not fully understood. In this study, we found that Arabidopsis eIF4E proteins are distributed in both the nucleus and cytoplasm, with only the cytoplasmic eIF4E being involved in the control of photoperiodic flowering. Genome-wide translation profiling using Ribo-tag sequencing reveals that eIF4E may regulate plant flowering by maintaining the homeostatic translation of components in the photoperiodic flowering pathway. eIF4E not only regulates the translation of flowering genes such as FLOWERING LOCUS T (FT) and FLOWERING LOCUS D (FLD) but also influences the translation of circadian genes like CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) and PSEUDO-RESPONSE REGULATOR 9 (PRR9). Consistently, our results show that the eIF4E modulates the rhythmic oscillation of the circadian clock. Together, our study provides mechanistic insights into how the protein translation regulates multiple developmental processes in Arabidopsis, including the circadian clock and photoperiodic flowering.

Regulation of alternative polyadenylation isoforms of Timp2 is an effector event of RAS signaling in cell transformation

Alternative polyadenylation (APA) generates mRNA isoforms with different lengths of the 3' untranslated region (3' UTR). The tissue inhibitor of metalloproteinase 2 (TIMP2) plays a key role in extracellular matrix remodeling under various developmental and disease conditions. Both human and mouse genes encoding TIMP2 contain two highly conserved 3'UTR APA sites, leading to mRNA isoforms that differ substantially in 3'UTR size. APA of Timp2 is one of the most significantly regulated events in multiple cell differentiation lineages. Here we show that Timp2 APA is highly regulated in transformation of NIH3T3 cells by the oncogene HRAS G12V . Perturbations of isoform expression with long 3'UTR isoform-specific knockdown or genomic removal of the alternative UTR (aUTR) region indicate that the long 3'UTR isoform contributes to the secreted Timp2 protein much more than the short 3'UTR isoform. The short and long 3'UTR isoforms differ in subcellular localization to endoplasmic reticulum (ER). Strikingly, Timp2 aUTR enhances secreted protein expression but no effect on intracellular proteins in reporter assays. Furthermore, downregulation of Timp2 long isoform mitigates gene expression changes elicited by HRAS G12V . Together, our data indicate that regulation of Timp2 protein expression through APA isoform changes is an integral part of RAS-mediated cell transformation and 3'UTR isoforms of Timp2 can have distinct impacts on expression of secreted vs. intracellular proteins.

Transcriptional profiling in microglia across physiological and pathological states identifies a transcriptional module associated with neurodegeneration

Microglia are the resident immune cells of the central nervous system and are involved in brain development, homeostasis, and disease. New imaging and genomics technologies are revealing microglial complexity across developmental and functional states, brain regions, and diseases. We curated a set of publicly available gene expression datasets from human microglia spanning disease and health to identify sets of genes reflecting physiological and pathological microglial states. We also integrated multiple human microglial single-cell RNA-seq datasets in Alzheimer's disease (AD), multiple sclerosis (MS), and Parkinson's disease, and identified a distinct microglial transcriptional signature shared across diseases. Analysis of germ-line DNA identified genes with variants associated with AD and MS that are overrepresented in microglial gene sets, including the disease-associated transcriptional signature. This work points to genes that are dysregulated in disease states and provides a resource for the analysis of diseases in which microglia are implicated by genetic evidence.

Mouse blood cells types and aging prediction using penalized Latent Dirichlet Allocation

BACKGROUND

Aging is a complex, heterogeneous process that has multiple causes. Knowledge on genomic, epigenomic and transcriptomic changes during the aging process shed light on understanding the aging mechanism. A recent breakthrough in biotechnology, single cell RNAseq, is revolutionizing aging study by providing gene expression profile of the entire transcriptome of individual cells. Many interesting information could be inferred from this new type of data with the help of novel computational methods.

RESULTS

In this manuscript a novel statistical method, penalized Latent Dirichlet Allocation (pLDA), is applied to an aging mouse blood scRNA-seq data set. A pipeline is built for cell type and aging prediction. The sequence of models in the pipeline take scRNA-seq expression counts as input, preprocess the data using pLDA and predict the cell type and aging status.

CONCLUSIONS

pLDA learns a dimension reduced representation of the expression profile. This representation allows identification of cell types and has predictability of the age of cells.

GOTermViewer: Visualization of Gene Ontology Enrichment in Multiple Differential Gene Expression Analyses

Gene ontology phrases are a widely used set of hierarchical terms that describe the biological properties of genes. These terms are then used to annotate individual genes, making it possible to determine the likely physiological properties of groups of genes such as a list of differentially expressed genes. Consequently, their ability to predict changes in biological features and functions based on alterations in gene expression has made gene ontology terms popular in the wide range of bioinformatic fields, such as differential gene expression and evolutionary biology. However, while they make the analysis easier, it is seldom easy to convey the results in a readily understandable manner. A number of applications have been developed to visualize gene ontology (GO) term enrichment; however, these solutions tend to focus on the display of aggregated results from a single analysis, making them unsuitable for the analysis of a series of experiments such as a time course or response to different drug treatments. As multiple pair wise comparisons are becoming a common feature of RNA profiling experiments, the absence of a mechanism to easily compare them is a significant problem. Consequently, to overcome this obstacle, we have developed GOTermViewer, an application that displays GO term enrichment data as determined by GOstats such that changes in physiological response across a number of individual analyses across a time course or range of drug treatments can be visualized.

Atomic vacancies of molybdenum disulfide nanoparticles stimulate mitochondrial biogenesis

Diminished mitochondrial function underlies many rare inborn errors of energy metabolism and contributes to more common age-associated metabolic and neurodegenerative disorders. Thus, boosting mitochondrial biogenesis has been proposed as a potential therapeutic approach for these diseases; however, currently we have a limited arsenal of compounds that can stimulate mitochondrial function. In this study, we designed molybdenum disulfide (MoS2) nanoflowers with predefined atomic vacancies that are fabricated by self-assembly of individual two-dimensional MoS2 nanosheets. Treatment of mammalian cells with MoS2 nanoflowers increased mitochondrial biogenesis by induction of PGC-1α and TFAM, which resulted in increased mitochondrial DNA copy number, enhanced expression of nuclear and mitochondrial-DNA encoded genes, and increased levels of mitochondrial respiratory chain proteins. Consistent with increased mitochondrial biogenesis, treatment with MoS2 nanoflowers enhanced mitochondrial respiratory capacity and adenosine triphosphate production in multiple mammalian cell types. Taken together, this study reveals that predefined atomic vacancies in MoS2 nanoflowers stimulate mitochondrial function by upregulating the expression of genes required for mitochondrial biogenesis.

Clusters of lineage-specific genes are anchored by ZNF274 in repressive perinucleolar compartments

Long known as the site of ribosome biogenesis, the nucleolus is increasingly recognized for its role in shaping three-dimensional (3D) genome organization. Still, the mechanisms governing the targeting of selected regions of the genome to nucleolus-associated domains (NADs) remain enigmatic. Here, we reveal the essential role of ZNF274, a SCAN-bearing member of the Krüppel-associated box (KRAB)-containing zinc finger protein (KZFP) family, in sequestering lineage-specific gene clusters within NADs. Ablation of ZNF274 triggers transcriptional activation across entire genomic neighborhoods-encompassing, among others, protocadherin and KZFP-encoding genes-with loss of repressive chromatin marks, altered the 3D genome architecture and de novo CTCF binding. Mechanistically, ZNF274 anchors target DNA sequences at the nucleolus and facilitates their compartmentalization via a previously uncharted function of the SCAN domain. Our findings illuminate the mechanisms underlying NAD organization and suggest that perinucleolar entrapment into repressive hubs constrains the activation of tandemly arrayed genes to enable selective expression and modulate cell differentiation programs during development.

AAGGG repeat expansions trigger RFC1-independent synaptic dysregulation in human CANVAS neurons

Cerebellar ataxia with neuropathy and vestibular areflexia syndrome (CANVAS) is a recessively inherited neurodegenerative disorder caused by intronic biallelic, nonreference CCCTT/AAGGG repeat expansions within RFC1. To investigate how these repeats cause disease, we generated patient induced pluripotent stem cell-derived neurons (iNeurons). CCCTT/AAGGG repeat expansions do not alter neuronal RFC1 splicing, expression, or DNA repair pathway function. In reporter assays, AAGGG repeats are translated into pentapeptide repeat proteins. However, these proteins and repeat RNA foci were not detected in iNeurons, and overexpression of these repeats failed to induce neuronal toxicity. CANVAS iNeurons exhibit defects in neuronal development and diminished synaptic connectivity that is rescued by CRISPR deletion of a single expanded AAGGG allele. These deficits were neither replicated by RFC1 knockdown in control iNeurons nor rescued by RFC1 reprovision in CANVAS iNeurons. These findings support a repeat-dependent but RFC1 protein-independent cause of neuronal dysfunction in CANVAS, with implications for therapeutic development in this currently untreatable condition.

Dissecting the genetic basis of resistance to Soil-borne cereal mosaic virus (SBCMV) in durum wheat by bi-parental mapping and GWAS

Soil-borne cereal mosaic virus (SBCMV), the causative agent of wheat mosaic, is a Furovirus challenging wheat production all over Europe. Differently from bread wheat, durum wheat shows greater susceptibility and stronger yield penalties, so identification and genetic characterization of resistance sources are major targets for durum genetics and breeding. The Sbm1 locus providing high level of resistance to SBCMV was mapped in bread wheat to the 5DL chromosome arm (Bass in Genome 49:1140-1148, 2006). This excluded the direct use of Sbm1 for durum wheat improvement. Only one major QTL has been mapped in durum wheat, namely QSbm.ubo-2B, on the 2BS chromosome region coincident with Sbm2, already known in bread wheat as reported (Bayles in HGCA Project Report, 2007). Therefore, QSbm.ubo-2B = Sbm2 is considered a pillar for growing durum in SBCMV-affected areas. Herein, we report the fine mapping of Sbm2 based on bi-parental mapping and GWAS, using the Infinium 90 K SNP array and high-throughput KASP®. Fine mapping pointed out a critical haploblock of 3.2 Mb defined by concatenated SNPs successfully converted to high-throughput KASP® markers coded as KUBO. The combination of KUBO-27, wPt-2106-ASO/HRM, KUBO-29, and KUBO-1 allows unequivocal tracing of the Sbm2-resistant haplotype. The interval harbors 52 high- and 41 low-confidence genes, encoding 17 cytochrome p450, three receptor kinases, two defensins, and three NBS-LRR genes. These results pave the way for Sbm2 positional cloning. Importantly, the development of Sbm2 haplotype tagging KASP® provides a valuable case study for improving efficacy of the European variety testing system and, ultimately, the decision-making process related to varietal characterization and choice.

Integrated morphological, physiological and transcriptomic analyses reveal the responses of Toona sinensis seedlings to low-nitrogen stress

Nitrogen is one of the most essential elements for plant growth and development. In this study, the growth, physiology, and transcriptome of Toona sinensis (A. Juss) Roem seedlings were compared between low-nitrogen (LN) and normal-nitrogen (NN) conditions. These results indicate that LN stress adversely influences T. sinensis seedling growth. The activities of key enzymes related to nitrogen assimilation and phytohormone contents were altered by LN stress. A total of 2828 differentially expressed genes (DEGs) in roots and 1547 in leaves were identified between the LN and NN treatments. A differential enrichment analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways indicated that nitrogen and sugar metabolism, flavonoid biosynthesis, plant hormone signal transduction, and ABC transporters, were strongly affected by LN stress. In summary, this research provides information for further understanding the response of T. sinensis to LN stress.

The transcription factor PagLBD4 represses cell differentiation and secondary cell wall biosynthesis in Populus

LBD (LATERAL ORGAN BOUNDARIES DOMAIN) transcription factors are key regulators of plant growth and development. In this study, we functionally characterized the PagLBD4 gene in Populus (Populus alba × Populus glandulosa). Overexpression of PagLBD4 (PagLBD4OE) significantly repressed secondary xylem differentiation and secondary cell wall (SCW) deposition, while CRISPR/Cas9-mediated PagLBD4 knockout (PagLBD4KO) significantly increased secondary xylem differentiation and SCW deposition. Consistent with the functional analysis, gene expression analysis revealed that SCW biosynthesis pathways were significantly down-regulated in PagLBD4OE plants but up-regulated in PagLBD4KO plants. We also performed DNA affinity purification followed by sequencing (DAP-seq) to identify genes bound by PagLBD4. Integration of RNA sequencing (RNA-seq) and DAP-seq data identified 263 putative direct target genes (DTGs) of PagLBD4, including important regulatory genes for SCW biosynthesis, such as PagMYB103 and PagIRX12. Together, our results demonstrated that PagLBD4 is a repressor of secondary xylem differentiation and SCW biosynthesis in Populus, which possibly lead to the dramatic growth repression in PagLBD4OE plants.