GenPipes: an open-source framework for distributed and scalable genomic analyses

Swim training induces distinct osseous gene expression patterns in anosteocytic and osteocytic teleost fish

The traditional understanding of bone mechanosensation implicates osteocytes, canaliculi, and the lacunocanalicular network in biomechanical adaptation. However, recent findings challenge this notion, as shown in advanced teleost fish where anosteocytic bone lacking osteocytes are nevertheless responsive to mechanical load. To investigate specific molecular mechanisms involved in bone mechanoadaptation in osteocytic and anosteocytic fish bone, we conducted a 5-min single swim-training experiment with zebrafish and ricefish, respectively. Through RNASeq analysis of fish spines, analyzed at various time points following swim training, we uncovered distinct gene expression patterns in osteocytic and anosteocytic fish bones. Notably, osteocytic fish bone exhibited an early response to mechanical load, contrasting to a delayed response observed in anosteocytic fish bones, both within 8 h following stimulation. We identified an increase in osteoblast differentiation in anosteocytic bone following training, while chordoblast activity was delayed. This temporal response suggests a time-dependent adaptation in anosteocytic bone, indicating the presence of intricate feedback networks within bone that lacks osteocytes.

TTF1 control of LncRNA synthesis delineates a tumor suppressor pathway directly regulating the ribosomal RNA genes

The tumor suppressor p14/19ARF regulates ribosomal RNA (rRNA) synthesis by controlling the nucleolar localization of Transcription Termination Factor 1 (TTF1). However, the role played by TTF1 in regulating the rRNA genes and in potentially controlling growth has remained unclear. We now show that TTF1 expression regulates cell growth by determining the cellular complement of ribosomes. Unexpectedly, it achieves this by acting as a "roadblock" to synthesis of the noncoding LncRNA and pRNA that we show are generated from the "Spacer Promoter" duplications present upstream of the 47S pre-rRNA promoter on the mouse and human ribosomal RNA genes. Unexpectedly, the endogenous generation of these noncoding RNAs does not induce CpG methylation or gene silencing. Rather, it acts in cis to suppress 47S preinitiation complex formation and hence de novo pre-rRNA synthesis by a mechanism reminiscent of promoter interference or occlusion. Taken together, our data delineate a pathway from p19ARF to cell growth suppression via the regulation of ribosome biogenesis by noncoding RNAs and validate a key cellular growth law in mammalian cells.

Genomic analysis of severe COVID-19 considering or not asthma comorbidity: GWAS insights from the BQC19 cohort

BACKGROUND

The severity of COVID-19 is influenced by various factors including the presence of respiratory diseases. Studies have indicated a potential relationship between asthma and COVID-19 severity.

OBJECTIVE

This study aimed to conduct a genome-wide association study (GWAS) to identify genetic and clinical variants associated with the severity of COVID-19, both among patients with and without asthma.

METHODS

We analyzed data from 2131 samples sourced from the Biobanque québécoise de la COVID-19 (BQC19), with 1499 samples from patients who tested positive for COVID-19. Among these, 1110 exhibited mild-to-moderate symptoms, 389 had severe symptoms, and 58 had asthma. We conducted a comparative analysis of clinical data from individuals in these three groups and GWAS using a logistic regression model. Phenotypic data analysis resulted in the refined covariates integrated into logistic models for genetic studies.

RESULTS

Considering a significance threshold of 1 × 10-6, seven genetic variants were associated with severe COVID-19. These variants were located proximal to five genes: sodium voltage-gated channel alpha subunit 1 (SCN10A), desmoplakin (DSP), RP1 axonemal microtubule associated (RP1), IGF like family member 1 (IGFL1), and docking protein 5 (DOK5). The GWAS comparing individuals with severe COVID-19 with asthma to those without asthma revealed four genetic variants in transmembrane protein with EGF like and two follistatin like domains 2 (TMEFF2) and huntingtin interacting protein-1 (HIP1) genes.

CONCLUSION

This study provides significant insights into the genetic profiles of patients with severe forms of the disease, whether accompanied by asthma or not. These findings enhance our comprehension of the genetic factors that affect COVID-19 severity.

KEY MESSAGES

Seven genetic variants were associated with the severe form of COVID-19; Four genetic variants were associated with the severe form of COVID-19 in individuals with comorbid asthma; These findings help define the genetic component of the severe form of COVID-19 in relation to asthma as a comorbidity.

Evolution of chromosome-arm aberrations in breast cancer through genetic network rewiring

The basal breast cancer subtype is enriched for triple-negative breast cancer (TNBC) and displays consistent large chromosomal deletions. Here, we characterize evolution and maintenance of chromosome 4p (chr4p) loss in basal breast cancer. Analysis of The Cancer Genome Atlas data shows recurrent deletion of chr4p in basal breast cancer. Phylogenetic analysis of a panel of 23 primary tumor/patient-derived xenograft basal breast cancers reveals early evolution of chr4p deletion. Mechanistically we show that chr4p loss is associated with enhanced proliferation. Gene function studies identify an unknown gene, C4orf19, within chr4p, which suppresses proliferation when overexpressed-a member of the PDCD10-GCKIII kinase module we name PGCKA1. Genome-wide pooled overexpression screens using a barcoded library of human open reading frames identify chromosomal regions, including chr4p, that suppress proliferation when overexpressed in a context-dependent manner, implicating network interactions. Together, these results shed light on the early emergence of complex aneuploid karyotypes involving chr4p and adaptive landscapes shaping breast cancer genomes.

Transcriptome analyses in juvenile yellow perch (Perca flavescens) exposed in vivo to clothianidin and chlorantraniliprole: Possible sampling bias

The St. Lawrence River is an important North American waterway that is subject to anthropogenic pressures including intensive urbanization, and agricultural development. Pesticides are widely used for agricultural activities in fields surrounding the yellow perch (Perca flavescens) habitat in Lake St. Pierre (Quebec, Canada), a fluvial lake of the river where the perch population has collapsed. Clothianidin and chlorantraniliprole were two of the most detected insecticides in surface waters near perch spawning areas. The objectives of the present study were to evaluate the transcriptional and biochemical effects of these two pesticides on juvenile yellow perch exposed for 28d to environmental doses of each compound alone and in a mixture under laboratory/aquaria conditions. Hepatic mRNA-sequencing revealed an effect of chlorantraniliprole alone (37 genes) and combined with clothianidin (251 genes), but no effects of clothianidin alone were observed in perch. Dysregulated genes were mostly related to circadian rhythms and to Ca2+ signaling, the latter effect has been previously associated with chlorantraniliprole mode of action in insects. Moreover, chronic exposure to clothianidin increased the activity of acetylcholinesterase in the brain of exposed fish, suggesting a potential non-target effect of this insecticide. Further analyses of three clock genes by qRT-PCR suggested that part of the observed effects of chlorantraniliprole on the circadian gene regulation of juvenile perch could be the result of time-of-day of sacrifice. These results provide insight into biological effects of insecticides in juvenile perch and highlight the importance of considering the circadian rhythm in experimental design and results analyses.

Host Immune Response Modulation in Avian Coronavirus Infection: Tracheal Transcriptome Profiling In Vitro and In Vivo

Infectious bronchitis virus (IBV) is a highly contagious Gammacoronavirus causing moderate to severe respiratory infection in chickens. Understanding the initial antiviral response in the respiratory mucosa is crucial for controlling viral spread. We aimed to characterize the impact of IBV Delmarva (DMV)/1639 and IBV Massachusetts (Mass) 41 at the primary site of infection, namely, in chicken tracheal epithelial cells (cTECs) in vitro and the trachea in vivo. We hypothesized that some elements of the induced antiviral responses are distinct in both infection models. We inoculated cTECs and infected young specific pathogen-free (SPF) chickens with IBV DMV/1639 or IBV Mass41, along with mock-inoculated controls, and studied the transcriptome using RNA-sequencing (RNA-seq) at 3 and 18 h post-infection (hpi) for cTECs and at 4 and 11 days post-infection (dpi) in the trachea. We showed that IBV DMV/1639 and IBV Mass41 replicate in cTECs in vitro and the trachea in vivo, inducing host mRNA expression profiles that are strain- and time-dependent. We demonstrated the different gene expression patterns between in vitro and in vivo tracheal IBV infection. Ultimately, characterizing host-pathogen interactions with various IBV strains reveals potential mechanisms for inducing and modulating the immune response during IBV infection in the chicken trachea.

An adapted protocol to derive microglia from stem cells and its application in the study of CSF1R-related disorders

BACKGROUND

Induced pluripotent stem cell-derived microglia (iMGL) represent an excellent tool in studying microglial function in health and disease. Yet, since differentiation and survival of iMGL are highly reliant on colony-stimulating factor 1 receptor (CSF1R) signaling, it is difficult to use iMGL to study microglial dysfunction associated with pathogenic defects in CSF1R.

METHODS

Serial modifications to an existing iMGL protocol were made, including but not limited to changes in growth factor combination to drive microglial differentiation, until successful derivation of microglia-like cells from an adult-onset leukoencephalopathy with axonal spheroids and pigmented glia (ALSP) patient carrying a c.2350G > A (p.V784M) CSF1R variant. Using healthy control lines, the quality of the new iMGL protocol was validated through cell yield assessment, measurement of microglia marker expression, transcriptomic comparison to primary microglia, and evaluation of inflammatory and phagocytic activities. Similarly, molecular and functional characterization of the ALSP patient-derived iMGL was carried out in comparison to healthy control iMGL.

RESULTS

The newly devised protocol allowed the generation of iMGL with enhanced transcriptomic similarity to cultured primary human microglia and with higher scavenging and inflammatory competence at ~ threefold greater yield compared to the original protocol. Using this protocol, decreased CSF1R autophosphorylation and cell surface expression was observed in iMGL derived from the ALSP patient compared to those derived from healthy controls. Additionally, ALSP patient-derived iMGL presented a migratory defect accompanying a temporal reduction in purinergic receptor P2Y12 (P2RY12) expression, a heightened capacity to internalize myelin, as well as heightened inflammatory response to Pam3CSK4. Poor P2RY12 expression was confirmed to be a consequence of CSF1R haploinsufficiency, as this feature was also observed following CSF1R knockdown or inhibition in mature control iMGL, and in CSF1RWT/KO and CSF1RWT/E633K iMGL compared to their respective isogenic controls.

CONCLUSIONS

We optimized a pre-existing iMGL protocol, generating a powerful tool to study microglial involvement in human neurological diseases. Using the optimized protocol, we have generated for the first time iMGL from an ALSP patient carrying a pathogenic CSF1R variant, with preliminary characterization pointing toward functional alterations in migratory, phagocytic and inflammatory activities.

Silencing ANGPTL8 reduces mouse preadipocyte differentiation and insulin signaling

ANGPTL8, expressed mainly in the liver and adipose tissue, regulates the activity of lipoprotein lipase (LPL) present in the extracellular space and triglyceride (TG) metabolism through its interaction with ANGPTL3 and ANGPTL4. Whether intracellular ANGPTL8 can also exert effects in tissues where it is expressed is uncertain. ANGPTL8 expression was low in preadipocytes and much increased during differentiation. To better understand the role of intracellular ANGPTL8 in adipocytes and assess whether it may play a role in adipocyte differentiation, we knocked down its expression in normal mouse subcutaneous preadipocytes. ANGPTL8 knockdown reduced adipocyte differentiation, cellular TG accumulation and also isoproterenol-stimulated lipolysis at day 7 of differentiation. RNA-Seq analysis of ANGPTL8 siRNA or control siRNA transfected SC preadipocytes on days 0, 2, 4 and 7 of differentiation showed that ANGPTL8 knockdown impeded the early (day 2) expression of adipogenic and insulin signaling genes, PPARγ, as well as genes related to extracellular matrix and NF-κB signaling. Insulin mediated Akt phosphorylation was reduced at an early stage during adipocyte differentiation. This study based on normal primary cells shows that ANGPTL8 has intracellular actions in addition to effects in the extracellular space, like modulating LPL activity. Preadipocyte ANGPTL8 expression modulates their differentiation possibly via changes in insulin signaling gene expression.

A host-directed oxadiazole compound potentiates antituberculosis treatment via zinc poisoning in human macrophages and in a mouse model of infection

Antituberculosis drugs, mostly developed over 60 years ago, combined with a poorly effective vaccine, have failed to eradicate tuberculosis. More worryingly, multiresistant strains of Mycobacterium tuberculosis (MTB) are constantly emerging. Innovative strategies are thus urgently needed to improve tuberculosis treatment. Recently, host-directed therapy has emerged as a promising strategy to be used in adjunct with existing or future antibiotics, by improving innate immunity or limiting immunopathology. Here, using high-content imaging, we identified novel 1,2,4-oxadiazole-based compounds, which allow human macrophages to control MTB replication. Genome-wide gene expression analysis revealed that these molecules induced zinc remobilization inside cells, resulting in bacterial zinc intoxication. More importantly, we also demonstrated that, upon treatment with these novel compounds, MTB became even more sensitive to antituberculosis drugs, in vitro and in vivo, in a mouse model of tuberculosis. Manipulation of heavy metal homeostasis holds thus great promise to be exploited to develop host-directed therapeutic interventions.

Sperm DNA methylation defects in a new mouse model of the 5,10-methylenetetrahydrofolate reductase 677C>T variant and correction with moderate dose folic acid supplementation

5,10-Methylenetetrahydrofolate reductase (MTHFR) is an enzyme that plays a key role in providing methyl groups for DNA methylation, including during spermatogenesis. A common genetic variant in humans (MTHFR 677C>T) results in reduced enzyme activity and has been linked to various disorders, including male infertility. A new animal model has been created by reproducing the human equivalent of the polymorphism in mice using CRISPR/Cas9. Biochemical parameters in the Mthfr 677TT mice recapitulate alterations found in MTHFR 677TT men. Our aims were to characterize the sperm DNA methylome of the Mthfr 677CC and TT mice on a control diet (2 mg folic acid/kg diet) and assess the effects of folic acid supplementation (10 mg/kg diet) on the sperm DNA methylome. Body and reproductive organ weights, testicular sperm counts, and histology were examined. DNA methylation in sperm was assessed using bisulfite pyrosequencing and whole-genome bisulfite sequencing (WGBS). Reproductive parameters and locus-specific imprinted gene methylation were unaffected by genotype or diet. Using WGBS, sperm from 677TT mice had 360 differentially methylated tiles as compared to 677CC mice, predominantly hypomethylation (60% of tiles). Folic acid supplementation mostly caused hypermethylation in sperm of males of both genotypes and was found to partially correct the DNA methylation alterations in sperm associated with the TT genotype. The new mouse model will be useful in understanding the role of MTHFR deficiency in male fertility and in designing folate supplementation regimens for the clinic.

Epigenetic variation impacts individual differences in the transcriptional response to influenza infection

Humans display remarkable interindividual variation in their immune response to identical challenges. Yet, our understanding of the genetic and epigenetic factors contributing to such variation remains limited. Here we performed in-depth genetic, epigenetic and transcriptional profiling on primary macrophages derived from individuals of European and African ancestry before and after infection with influenza A virus. We show that baseline epigenetic profiles are strongly predictive of the transcriptional response to influenza A virus across individuals. Quantitative trait locus (QTL) mapping revealed highly coordinated genetic effects on gene regulation, with many cis-acting genetic variants impacting concomitantly gene expression and multiple epigenetic marks. These data reveal that ancestry-associated differences in the epigenetic landscape can be genetically controlled, even more than gene expression. Lastly, among QTL variants that colocalized with immune-disease loci, only 7% were gene expression QTL, while the remaining genetic variants impact epigenetic marks, stressing the importance of considering molecular phenotypes beyond gene expression in disease-focused studies.

Intra-Host Evolution Analyses in an Immunosuppressed Patient Supports SARS-CoV-2 Viral Reservoir Hypothesis

Throughout the SARS-CoV-2 pandemic, several variants of concern (VOCs) have been identified, many of which share recurrent mutations in the spike glycoprotein's receptor-binding domain (RBD). This region coincides with known epitopes and can therefore have an impact on immune escape. Protracted infections in immunosuppressed patients have been hypothesized to lead to an enrichment of such mutations and therefore drive evolution towards VOCs. Here, we present the case of an immunosuppressed patient that developed distinct populations with immune escape mutations throughout the course of their infection. Notably, by investigating the co-occurrence of substitutions on individual sequencing reads in the RBD, we found quasispecies harboring mutations that confer resistance to known monoclonal antibodies (mAbs) such as S:E484K and S:E484A. These mutations were acquired without the patient being treated with mAbs nor convalescent sera and without them developing a detectable immune response to the virus. We also provide additional evidence for a viral reservoir based on intra-host phylogenetics, which led to a viral substrain that evolved elsewhere in the patient's body, colonizing their upper respiratory tract (URT). The presence of SARS-CoV-2 viral reservoirs can shed light on protracted infections interspersed with periods where the virus is undetectable, and potential explanations for long-COVID cases.

Adaptive laboratory evolution reveals regulators involved in repressing biofilm development as key players in Bacillus subtilis root colonization

Root-associated microorganisms play an important role in plant health, such as plant growth-promoting rhizobacteria (PGPR) from the Bacillus and Pseudomonas genera. Although bacterial consortia including these two genera would represent a promising avenue to efficient biofertilizer formulation, we observed that Bacillus subtilis root colonization is decreased by the presence of Pseudomonas fluorescens and Pseudomonas protegens. To determine if B. subtilis can adapt to the inhibitory effect of Pseudomonas on roots, we conducted adaptative laboratory evolution experiments with B. subtilis in mono-association or co-cultured with P. fluorescens on tomato plant roots. Evolved isolates with various colony morphology and stronger colonization capacity of both tomato plant and Arabidopsis thaliana roots emerged rapidly from the two evolution experiments. Certain evolved isolates also had better fitness on the root in the presence of other Pseudomonas species. In all independent lineages, whole-genome resequencing revealed non-synonymous mutations in genes ywcC or sinR encoding regulators involved in repressing biofilm development, suggesting their involvement in enhanced root colonization. These findings provide insights into the molecular mechanisms underlying B. subtilis adaptation to root colonization and highlight the potential of directed evolution to enhance the beneficial traits of PGPR.IMPORTANCEIn this study, we aimed to enhance the abilities of the plant-beneficial bacterium Bacillus subtilis to colonize plant roots in the presence of competing Pseudomonas bacteria. To achieve this, we conducted adaptive laboratory experiments, allowing Bacillus to evolve in a defined environment. We successfully obtained strains of Bacillus that were more effective at colonizing plant roots than the ancestor strain. To identify the genetic changes driving this improvement, we sequenced the genomes of these evolved strains. Interestingly, mutations that facilitated the formation of robust biofilms on roots were predominant. Many of these evolved Bacillus isolates also displayed the remarkable ability to outcompete Pseudomonas species. Our research sheds light on the mutational paths selected in Bacillus subtilis to thrive in root environments and offers exciting prospects for improving beneficial traits in plant growth-promoting microorganisms. Ultimately, this could pave the way for the development of more effective biofertilizers and sustainable agricultural practices.

MerTK is a mediator of alpha-synuclein fibril uptake by human microglia

Mer tyrosine kinase (MerTK) is a receptor tyrosine kinase that mediates non-inflammatory, homeostatic phagocytosis of diverse types of cellular debris. Highly expressed on the surface of microglial cells, MerTK is of importance in brain development, homeostasis, plasticity and disease. Yet, involvement of this receptor in the clearance of protein aggregates that accumulate with ageing and in neurodegenerative diseases has yet to be defined. The current study explored the function of MerTK in the microglial uptake of alpha-synuclein fibrils which play a causative role in the pathobiology of synucleinopathies. Using human primary and induced pluripotent stem cell-derived microglia, the MerTK-dependence of alpha-synuclein fibril internalization was investigated in vitro. Relevance of this pathway in synucleinopathies was assessed through burden analysis of MERTK variants and analysis of MerTK expression in patient-derived cells and tissues. Pharmacological inhibition of MerTK and siRNA-mediated MERTK knockdown both caused a decreased rate of alpha-synuclein fibril internalization by human microglia. Consistent with the non-inflammatory nature of MerTK-mediated phagocytosis, alpha-synuclein fibril internalization was not observed to induce secretion of pro-inflammatory cytokines such as IL-6 or TNF, and downmodulated IL-1β secretion from microglia. Burden analysis in two independent patient cohorts revealed a significant association between rare functionally deleterious MERTK variants and Parkinson's disease in one of the cohorts (P = 0.002). Despite a small upregulation in MERTK mRNA expression in nigral microglia from Parkinson's disease/Lewy body dementia patients compared to those from non-neurological control donors in a single-nuclei RNA-sequencing dataset (P = 5.08 × 10-21), no significant upregulation in MerTK protein expression was observed in human cortex and substantia nigra lysates from Lewy body dementia patients compared to controls. Taken together, our findings define a novel role for MerTK in mediating the uptake of alpha-synuclein fibrils by human microglia, with possible involvement in limiting alpha-synuclein spread in synucleinopathies such as Parkinson's disease. Upregulation of this pathway in synucleinopathies could have therapeutic values in enhancing alpha-synuclein fibril clearance in the brain.

Sex-based disparities in DNA methylation and gene expression in late-gestation mouse placentas

BACKGROUND

The placenta is vital for fetal development and its contributions to various developmental issues, such as pregnancy complications, fetal growth restriction, and maternal exposure, have been extensively studied in mice. The placenta forms mainly from fetal tissue and therefore has the same biological sex as the fetus it supports. Extensive research has delved into the placenta's involvement in pregnancy complications and future offspring development, with a notable emphasis on exploring sex-specific disparities. However, despite these investigations, sex-based disparities in epigenetic (e.g., DNA methylation) and transcriptomic features of the late-gestation mouse placenta remain largely unknown.

METHODS

We collected male and female mouse placentas at late gestation (E18.5, n = 3/sex) and performed next-generation sequencing to identify genome-wide sex differences in transcription and DNA methylation.

RESULTS

Our comparison between male and female revealed 358 differentially expressed genes (DEGs) on autosomes, which were associated with signaling pathways involved in transmembrane transport and the responses to viruses and external stimuli. X chromosome DEGs (n = 39) were associated with different pathways, including those regulating chromatin modification and small GTPase-mediated signal transduction. Differentially methylated regions (DMRs) were more common on the X chromosomes (n = 3756) than on autosomes (n = 1705). Interestingly, while most X chromosome DMRs had higher DNA methylation levels in female placentas and tended to be included in CpG dinucleotide-rich regions, 73% of autosomal DMRs had higher methylation levels in male placentas and were distant from CpG-rich regions. Several DEGs were correlated with DMRs. A subset of the DMRs present in late-stage placentas were already established in mid-gestation (E10.5) placentas (n = 348 DMRs on X chromosome and 19 DMRs on autosomes), while others were acquired later in placental development.

CONCLUSION

Our study provides comprehensive lists of DEGs and DMRs between male and female that collectively cause profound differences in the DNA methylation and gene expression profiles of late-gestation mouse placentas. Our results demonstrate the importance of incorporating sex-specific analyses into epigenetic and transcription studies to enhance the accuracy and comprehensiveness of their conclusions and help address the significant knowledge gap regarding how sex differences influence placental function.

The Association between Long-Term DDT or DDE Exposures and an Altered Sperm Epigenome-a Cross-Sectional Study of Greenlandic Inuit and South African VhaVenda Men

BACKGROUND

The organochlorine dichlorodiphenyltrichloroethane (DDT) is banned worldwide owing to its negative health effects. It is exceptionally used as an insecticide for malaria control. Exposure occurs in regions where DDT is applied, as well as in the Arctic, where its endocrine disrupting metabolite, p , p ' -dichlorodiphenyldichloroethylene (p , p ' -DDE) accumulates in marine mammals and fish. DDT and p , p ' -DDE exposures are linked to birth defects, infertility, cancer, and neurodevelopmental delays. Of particular concern is the potential of DDT use to impact the health of generations to come via the heritable sperm epigenome.

OBJECTIVES

The objective of this study was to assess the sperm epigenome in relation to p , p ' -DDE serum levels between geographically diverse populations.

METHODS

In the Limpopo Province of South Africa, we recruited 247 VhaVenda South African men and selected 50 paired blood serum and semen samples, and 47 Greenlandic Inuit blood and semen paired samples were selected from a total of 193 samples from the biobank of the INUENDO cohort, an EU Fifth Framework Programme Research and Development project. Sample selection was based on obtaining a range of p , p ' -DDE serum levels (mean = 870.734 ± 134.030  ng / mL ). We assessed the sperm epigenome in relation to serum p , p ' -DDE levels using MethylC-Capture-sequencing (MCC-seq) and chromatin immunoprecipitation followed by sequencing (ChIP-seq). We identified genomic regions with altered DNA methylation (DNAme) and differential enrichment of histone H3 lysine 4 trimethylation (H3K4me3) in sperm.

RESULTS

Differences in DNAme and H3K4me3 enrichment were identified at transposable elements and regulatory regions involved in fertility, disease, development, and neurofunction. A subset of regions with sperm DNAme and H3K4me3 that differed between exposure groups was predicted to persist in the preimplantation embryo and to be associated with embryonic gene expression.

DISCUSSION

These findings suggest that DDT and p , p ' -DDE exposure impacts the sperm epigenome in a dose-response-like manner and may negatively impact the health of future generations through epigenetic mechanisms. Confounding factors, such as other environmental exposures, genetic diversity, and selection bias, cannot be ruled out. https://doi.org/10.1289/EHP12013.

Draft genomes of 3 cyanobacteria strains and 17 co-habiting proteobacteria assembled from metagenomes

Cyanobium and Synechococcus are prominent, globally distributed cyanobacteria genera with ecological significance. Here, we report the genomes of the marine Synechococcus sp. CCMP836 and two strains of Cyanobium (CZS25K and CZS48M) along with the genomes of 17 co-occurring proteobacteria. These genomes will improve the strain-specific ecological positions.

Clinical Implementation of MetaFusion for Accurate Cancer-Driving Fusion Detection from RNA Sequencing

Oncogenic fusion genes may be identified from next-generation sequencing data, typically RNA-sequencing. However, in a clinical setting, identifying these alterations is challenging against a background of nonrelevant fusion calls that reduce workflow precision and specificity. Furthermore, although numerous algorithms have been developed to detect fusions in RNA-sequencing, there are variations in their individual sensitivities. Here this problem was addressed by introducing MetaFusion into clinical use. Its utility was illustrated when applied to both whole-transcriptome and targeted sequencing data sets. MetaFusion combines ensemble fusion calls from eight individual fusion-calling algorithms with practice-informed identification of gene fusions that are known to be clinically relevant. In doing so, it allows oncogenic fusions to be identified with near-perfect sensitivity and high precision and specificity, significantly outperforming the individual fusion callers it uses as well as existing clinical-grade software. MetaFusion enhances clinical yield over existing methods and is able to identify fusions that have patient relevance for the purposes of diagnosis, prognosis, and treatment.

Genome-wide association study of placental weight identifies distinct and shared genetic influences between placental and fetal growth

A well-functioning placenta is essential for fetal and maternal health throughout pregnancy. Using placental weight as a proxy for placental growth, we report genome-wide association analyses in the fetal (n = 65,405), maternal (n = 61,228) and paternal (n = 52,392) genomes, yielding 40 independent association signals. Twenty-six signals are classified as fetal, four maternal and three fetal and maternal. A maternal parent-of-origin effect is seen near KCNQ1. Genetic correlation and colocalization analyses reveal overlap with birth weight genetics, but 12 loci are classified as predominantly or only affecting placental weight, with connections to placental development and morphology, and transport of antibodies and amino acids. Mendelian randomization analyses indicate that fetal genetically mediated higher placental weight is causally associated with preeclampsia risk and shorter gestational duration. Moreover, these analyses support the role of fetal insulin in regulating placental weight, providing a key link between fetal and placental growth.

Differential gene expression provides leads to environmentally regulated soybean seed protein content

Soybean is an important global source of plant-based protein. A persistent trend has been observed over the past two decades that soybeans grown in western Canada have lower seed protein content than soybeans grown in eastern Canada. In this study, 10 soybean genotypes ranging in average seed protein content were grown in an eastern location (control) and three western locations (experimental) in Canada. Seed protein and oil contents were measured for all lines in each location. RNA-sequencing and differential gene expression analysis were used to identify differentially expressed genes that may account for relatively low protein content in western-grown soybeans. Differentially expressed genes were enriched for ontologies and pathways that included amino acid biosynthesis, circadian rhythm, starch metabolism, and lipid biosynthesis. Gene ontology, pathway mapping, and quantitative trait locus (QTL) mapping collectively provide a close inspection of mechanisms influencing nitrogen assimilation and amino acid biosynthesis between soybeans grown in the East and West. It was found that western-grown soybeans had persistent upregulation of asparaginase (an asparagine hydrolase) and persistent downregulation of asparagine synthetase across 30 individual differential expression datasets. This specific difference in asparagine metabolism between growing environments is almost certainly related to the observed differences in seed protein content because of the positive correlation between seed protein content at maturity and free asparagine in the developing seed. These results provided pointed information on seed protein-related genes influenced by environment. This information is valuable for breeding programs and genetic engineering of geographically optimized soybeans.

Establishment of paternal methylation imprint at the H19/Igf2 imprinting control region

The insulator model explains the workings of the H19 and Igf2 imprinted domain in the soma, where insulation of the Igf2 promoter from its enhancers occurs by CTCF in the maternally inherited unmethylated chromosome but not the paternally inherited methylated allele. The molecular mechanism that targets paternal methylation imprint establishment to the imprinting control region (ICR) in the male germline is unknown. We tested the function of prospermatogonia-specific broad low-level transcription in this process using mouse genetics. Paternal imprint establishment was abnormal when transcription was stopped at the entry point to the ICR. The germline epimutation persisted into the paternal allele of the soma, resulting in reduced Igf2 in fetal organs and reduced fetal growth, consistent with the insulator model and insulin-like growth factor 2 (IGF2)'s role as fetal growth factor. These results collectively support the role of broad low-level transcription through the H19/Igf2 ICR in the establishment of its paternal methylation imprint in the male germ line, with implications for Silver-Russell syndrome.

SHP-1 phosphatase acts as a coactivator of PCK1 transcription to control gluconeogenesis

We previously reported that the protein-tyrosine phosphatase SHP-1 (PTPN6) negatively regulates insulin signaling, but its impact on hepatic glucose metabolism and systemic glucose control remains poorly understood. Here, we use co-immunoprecipitation assays, chromatin immunoprecipitation sequencing, in silico methods, and gluconeogenesis assay, and found a new mechanism whereby SHP-1 acts as a coactivator for transcription of the phosphoenolpyruvate carboxykinase 1 (PCK1) gene to increase liver gluconeogenesis. SHP-1 is recruited to the regulatory regions of the PCK1 gene and interacts with RNA polymerase II. The recruitment of SHP-1 to chromatin is dependent on its association with the transcription factor signal transducer and activator of transcription 5 (STAT5). Loss of SHP-1 as well as STAT5 decrease RNA polymerase II recruitment to the PCK1 promoter and consequently PCK1 mRNA levels leading to blunted gluconeogenesis. This work highlights a novel nuclear role of SHP-1 as a key transcriptional regulator of hepatic gluconeogenesis adding a new mechanism to the repertoire of SHP-1 functions in metabolic control.

CdGAP is a talin-binding protein and a target of TGF-β signaling that promotes HER2-positive breast cancer growth and metastasis

Epithelial-to-mesenchymal transition (EMT) plays a crucial role in metastasis, which is the leading cause of death in breast cancer patients. Here, we show that Cdc42 GTPase-activating protein (CdGAP) promotes tumor formation and metastasis to lungs in the HER2-positive (HER2+) murine breast cancer model. CdGAP facilitates intravasation, extravasation, and growth at metastatic sites. CdGAP depletion in HER2+ murine primary tumors mediates crosstalk with a Dlc1-RhoA pathway and is associated with a transforming growth factor β (TGF-β)-induced EMT transcriptional signature. CdGAP is positively regulated by TGF-β signaling during EMT and interacts with the adaptor talin to modulate focal adhesion dynamics and integrin activation. Moreover, HER2+ breast cancer patients with high CdGAP mRNA expression combined with a high TGF-β-EMT signature are more likely to present lymph node invasion. Our results suggest CdGAP as a candidate therapeutic target for HER2+ metastatic breast cancer by inhibiting TGF-β and integrin/talin signaling pathways.

Unlocking Cowpea's Defense Responses: Conserved Transcriptional Signatures in the Battle against CABMV and CPSMV Viruses

Cowpea aphid-borne mosaic virus (CABMV) and Cowpea severe mosaic virus (CPSMV) threaten cowpea commercial production. This study aimed to analyze Conserved Transcriptional Signatures (CTS) in cowpea's genotypes that are resistant to these viruses. CTS covered up- (UR) or down-regulated (DR) cowpea transcripts in response to CABMV and CPSMV mechanical inoculations. The conservation of cowpea's UR defense response was primarily observed with the one hpi treatments, with decreased CTS representatives as time elapsed. This suggests that cowpea utilizes generic mechanisms during its early interaction with the studied viruses, and subsequently employs more specialized strategies for each viral agent. The potential action of the CTS-UR emphasizes the importance of redox balance, ethylene and jasmonic acid pathways. Additionally, the CTS-UR provides evidence for the involvement of R genes, PR proteins, and PRRs receptors-extensively investigated in combating bacterial and fungal pathogens-in the defense against viral inoculation. AP2-ERF, WRKY, and MYB transcription factors, as well as PIP aquaporins and MAPK cascades, also emerged as significant molecular players. The presented work represents the first study investigating conserved mechanisms in the cowpea defense response to viral inoculations, highlighting relevant processes for initial defense responses.

Gpr18 agonist dampens inflammation, enhances myogenesis, and restores muscle function in models of Duchenne muscular dystrophy

Introduction: Muscle wasting in Duchenne Muscular Dystrophy is caused by myofiber fragility and poor regeneration that lead to chronic inflammation and muscle replacement by fibrofatty tissue. Our recent findings demonstrated that Resolvin-D2, a bioactive lipid derived from omega-3 fatty acids, has the capacity to dampen inflammation and stimulate muscle regeneration to alleviate disease progression. This therapeutic avenue has many advantages compared to glucocorticoids, the current gold-standard treatment for Duchenne Muscular Dystrophy. However, the use of bioactive lipids as therapeutic drugs also faces many technical challenges such as their instability and poor oral bioavailability. Methods: Here, we explored the potential of PSB-KD107, a synthetic agonist of the resolvin-D2 receptor Gpr18, as a therapeutic alternative for Duchenne Muscular Dystrophy. Results and discussion: We showed that PSB-KD107 can stimulate the myogenic capacity of patient iPSC-derived myoblasts in vitro. RNAseq analysis revealed an enrichment in biological processes related to fatty acid metabolism, lipid biosynthesis, small molecule biosynthesis, and steroid-related processes in PSB-KD107-treated mdx myoblasts, as well as signaling pathways such as Peroxisome proliferator-activated receptors, AMP-activated protein kinase, mammalian target of rapamycin, and sphingolipid signaling pathways. In vivo, the treatment of dystrophic mdx mice with PSB-KD107 resulted in reduced inflammation, enhanced myogenesis, and improved muscle function. The positive impact of PSB-KD107 on muscle function is similar to the one of Resolvin-D2. Overall, our findings provide a proof-of concept that synthetic analogs of bioactive lipid receptors hold therapeutic potential for the treatment of Duchenne Muscular Dystrophy.

Steroidogenic factor 1 (SF-1; Nr5a1) regulates the formation of the ovarian reserve

The ovarian follicle reserve, formed pre- or perinatally, comprises all oocytes for lifetime reproduction. Depletion of this reserve results in infertility. Steroidogenic factor 1 (SF-1; Nr5a1) and liver receptor homolog 1 (LRH-1; Nr5a2) are two orphan nuclear receptors that regulate adult endocrine function, but their role in follicle formation is unknown. We developed models of conditional depletion of SF-1 or LRH-1 from prenatal ovaries. Depletion of SF-1, but not LRH-1, resulted in dramatically smaller ovaries and fewer primordial follicles. This was mediated by increased oocyte death, resulting from increased ovarian inflammation and increased Notch signaling. Major dysregulated genes were Iroquois homeobox 3 and 5 and their downstream targets involved in the establishment of the ovarian laminin matrix and oocyte-granulosa cell gap junctions. Disruptions of these pathways resulted in follicles with impaired basement membrane formation and compromised oocyte-granulosa communication networks, believed to render them more prone to atresia. This study identifies SF-1 as a key regulator of the formation of the ovarian reserve.

Linked-read based analysis of the medulloblastoma genome

Introduction

Medulloblastoma is the most common type of malignant pediatric brain tumor with group 4 medulloblastomas (G4 MBs) accounting for 40% of cases. However, the molecular mechanisms that underlie this subgroup are still poorly understood. Point mutations are detected in a large number of genes at low incidence per gene while the detection of complex structural variants in recurrently affected genes typically requires the application of long-read technologies.

Methods

Here, we applied linked-read sequencing, which combines the long-range genome information of long-read sequencing with the high base pair accuracy of short read sequencing and very low sample input requirements.

Results

We demonstrate the detection of complex structural variants and point mutations in these tumors, and, for the first time, the detection of extrachromosomal DNA (ecDNA) with linked-reads. We provide further evidence for the high heterogeneity of somatic mutations in G4 MBs and add new complex events associated with it.

Discussion

We detected several enhancer-hijacking events, an ecDNA containing the MYCN gene, and rare structural rearrangements, such a chromothripsis in a G4 medulloblastoma, chromoplexy involving 8 different chromosomes, a TERT gene rearrangement, and a PRDM6 duplication.

Glucocorticoid stimulation induces regionalized gene responses within topologically associating domains

Transcription-factor binding to cis-regulatory regions regulates the gene expression program of a cell, but occupancy is often a poor predictor of the gene response. Here, we show that glucocorticoid stimulation led to the reorganization of transcriptional coregulators MED1 and BRD4 within topologically associating domains (TADs), resulting in active or repressive gene environments. Indeed, we observed a bias toward the activation or repression of a TAD when their activities were defined by the number of regions gaining and losing MED1 and BRD4 following dexamethasone (Dex) stimulation. Variations in Dex-responsive genes at the RNA levels were consistent with the redistribution of MED1 and BRD4 at the associated cis-regulatory regions. Interestingly, Dex-responsive genes without the differential recruitment of MED1 and BRD4 or binding by the glucocorticoid receptor were found within TADs, which gained or lost MED1 and BRD4, suggesting a role of the surrounding environment in gene regulation. However, the amplitude of the response of Dex-regulated genes was higher when the differential recruitment of the glucocorticoid receptor and transcriptional coregulators was observed, reaffirming the role of transcription factor-driven gene regulation and attributing a lesser role to the TAD environment. These results support a model where a signal-induced transcription factor induces a regionalized effect throughout the TAD, redefining the notion of direct and indirect effects of transcription factors on target genes.

Chronic treatment with ATR and CHK1 inhibitors does not substantially increase the mutational burden of human cells

DNA replication stress (RS) entails the frequent slow down and arrest of replication forks by a variety of conditions that hinder accurate and processive genome duplication. Elevated RS leads to genome instability, replication catastrophe and eventually cell death. RS is particularly prevalent in cancer cells and its exacerbation to unsustainable levels by chemotherapeutic agents remains a cornerstone of cancer treatments. The adverse consequences of RS are normally prevented by the ATR and CHK1 checkpoint kinases that stabilize stressed forks, suppress origin firing and promote cell cycle arrest when replication is perturbed. Specific inhibitors of these kinases have been developed and shown to potentiate RS and cell death in multiple in vitro cancer settings. Ongoing clinical trials are now probing their efficacy against various cancer types, either as single agents or in combination with mainstay chemotherapeutics. Despite their promise as valuable additions to the anti-cancer pharmacopoeia, we still lack a genome-wide view of the potential mutagenicity of these new drugs. To investigate this question, we performed chronic long-term treatments of TP53-depleted human cancer cells with ATR and CHK1 inhibitors (ATRi, AZD6738/ceralasertib and CHK1i, MK8776/SCH-900776). ATR or CHK1 inhibition did not significantly increase the mutational burden of cells, nor generate specific mutational signatures. Indeed, no notable changes in the numbers of base substitutions, short insertions/deletions and larger scale rearrangements were observed despite induction of replication-associated DNA breaks during treatments. Interestingly, ATR inhibition did induce a slight increase in closely-spaced mutations, a feature previously attributed to translesion synthesis DNA polymerases. The results suggest that ATRi and CHK1i do not have substantial mutagenic effects in vitro when used as standalone agents.

Integration of environmental signatures and omics-based approaches on the European flounder to assist with health assessment of estuarine ecosystems in Brittany, France

This study aimed to develop a multidisciplinary approach to assess the ecological status of six moderate-sized French estuaries. For each estuary, we gathered geographical information, hydrobiological data, chemistry of pollutants and fish biology, including integration of proteomics and transcriptomics data. This integrative study covered the entire hydrological system studied, from the watershed to the estuary, and considered all the anthropogenic factors that can impact this environment. To reach this goal, European flounder (Platichthys flesus) were collected from six estuaries in September, which ensures a minimum residence time of five months within an estuary. Geographical metrics are used to characterize land use in each watershed. The concentrations of nitrite, nitrate, organic pollutants, and trace elements were measured in water, sediments and biota. All of these environmental parameters allowed to set up a typology of estuaries. Classical fish biomarkers, coupled with molecular data from transcriptomics and shotgun proteomics, highlighted the flounder's responses to stressors in its environment. We analysed the protein abundances and gene expression levels in the liver of fish from the different estuaries. We showed clear positive deregulation of proteins associated with xenobiotic detoxification in a system characterized by a large population density and industrial activity, as well as in a predominantly agricultural catchment area (mostly cultures of vegetables and pig breeding) mainly impacted by pesticides. Fish from the latter estuary also displayed strong deregulation of the urea cycle, most probably related to high nitrogen load. Proteomic and transcriptomic data also revealed a deregulation of proteins and genes related to the response to hypoxia, and a probable endocrine disruption in some estuaries. Coupling these data allowed the precise identification of the main stressors interacting within each hydrosystem.

Antigen-driven colonic inflammation is associated with development of dysplasia in primary sclerosing cholangitis

Primary sclerosing cholangitis (PSC) is an immune-mediated disease of the bile ducts that co-occurs with inflammatory bowel disease (IBD) in almost 90% of cases. Colorectal cancer is a major complication of patients with PSC and IBD, and these patients are at a much greater risk compared to patients with IBD without concomitant PSC. Combining flow cytometry, bulk and single-cell transcriptomics, and T and B cell receptor repertoire analysis of right colon tissue from 65 patients with PSC, 108 patients with IBD and 48 healthy individuals we identified a unique adaptive inflammatory transcriptional signature associated with greater risk and shorter time to dysplasia in patients with PSC. This inflammatory signature is characterized by antigen-driven interleukin-17A (IL-17A)+ forkhead box P3 (FOXP3)+ CD4 T cells that express a pathogenic IL-17 signature, as well as an expansion of IgG-secreting plasma cells. These results suggest that the mechanisms that drive the emergence of dysplasia in PSC and IBD are distinct and provide molecular insights that could guide prevention of colorectal cancer in individuals with PSC.

Capturing sex-specific and hypofertility-linked effects of assisted reproductive technologies on the cord blood DNA methylome

BACKGROUND

Children conceived through assisted reproduction are at an increased risk for growth and genomic imprinting disorders, often linked to DNA methylation defects. It has been suggested that assisted reproductive technology (ART) and underlying parental infertility can induce epigenetic instability, specifically interfering with DNA methylation reprogramming events during germ cell and preimplantation development. To date, human studies exploring the association between ART and DNA methylation defects have reported inconsistent or inconclusive results, likely due to population heterogeneity and the use of technologies with limited coverage of the epigenome. In our study, we explored the epigenetic risk of ART by comprehensively profiling the DNA methylome of 73 human cord blood samples of singleton pregnancies (n = 36 control group, n = 37 ART/hypofertile group) from a human prospective longitudinal birth cohort, the 3D (Design, Develop, Discover) Study, using a high-resolution sequencing-based custom capture panel that examines over 2.4 million autosomal CpGs in the genome.

RESULTS

We identified evidence of sex-specific effects of ART/hypofertility on cord blood DNA methylation patterns. Our genome-wide analyses identified ~ 46% more CpGs affected by ART/hypofertility in female than in male infant cord blood. We performed a detailed analysis of three imprinted genes which have been associated with altered DNA methylation following ART (KCNQ1OT1, H19/IGF2 and GNAS) and found that female infant cord blood was associated with DNA hypomethylation. When compared to less invasive procedures such as intrauterine insemination, more invasive ARTs (in vitro fertilization, intracytoplasmic sperm injection, embryo culture) resulted in more marked and distinct effects on the cord blood DNA methylome. In the in vitro group, we found a close to fourfold higher proportion of significantly enriched Gene Ontology terms involved in development than in the in vivo group.

CONCLUSIONS

Our study highlights the ability of a sensitive, targeted, sequencing-based approach to uncover DNA methylation perturbations in cord blood associated with hypofertility and ART and influenced by offspring sex and ART technique invasiveness.

Systematic comparison of culture media uncovers phenotypic shift of primary human microglia defined by reduced reliance to CSF1R signaling

Efforts to understand microglia function in health and diseases have been hindered by the lack of culture models that recapitulate in situ cellular properties. In recent years, the use of serum-free media with brain-derived growth factors (colony stimulating factor 1 receptor [CSF1R] ligands and TGF-β1/2) have been favored for the maintenance of rodent microglia as they promote morphological features observed in situ. Here we study the functional and transcriptomic impacts of such media on human microglia (hMGL). Media formulation had little impact on microglia transcriptome assessed by RNA sequencing which was sufficient to significantly alter microglia capacity to phagocytose myelin debris and to elicit an inflammatory response to lipopolysaccharide. When compared to immediately ex vivo microglia from the same donors, the addition of fetal bovine serum to culture media, but not growth factors, was found to aid in the maintenance of key signature genes including those involved in phagocytic processes. A phenotypic shift characterized by CSF1R downregulation in culture correlated with a lack of reliance on CSF1R signaling for survival. Consequently, no improvement in cell survival was observed following culture supplementation with CSF1R ligands. Our study provides better understanding of hMGL in culture, with observations that diverge from those previously made in rodent microglia.

Application of Genomic Sequencing to Refine Patient Stratification for Adjuvant Therapy in Renal Cell Carcinoma

PURPOSE

Patients with resected localized clear-cell renal cell carcinoma (ccRCC) remain at variable risk of recurrence. Incorporation of biomarkers may refine risk prediction and inform adjuvant treatment decisions. We explored the role of tumor genomics in this setting, leveraging the largest cohort to date of localized ccRCC tissues subjected to targeted gene sequencing.

EXPERIMENTAL DESIGN

The somatic mutation status of 12 genes was determined in 943 ccRCC cases from a multinational cohort of patients, and associations to outcomes were examined in a Discovery (n = 469) and Validation (n = 474) framework.

RESULTS

Tumors containing a von-Hippel Lindau (VHL) mutation alone were associated with significantly improved outcomes in comparison with tumors containing a VHL plus additional mutations. Within the Discovery cohort, those with VHL+0, VHL+1, VHL+2, and VHL+≥3 tumors had disease-free survival (DFS) rates of 90.8%, 80.1%, 68.2%, and 50.7% respectively, at 5 years. This trend was replicated in the Validation cohort. Notably, these genomically defined groups were independent of tumor mutational burden. Amongst patients eligible for adjuvant therapy, those with a VHL+0 tumor (29%) had a 5-year DFS rate of 79.3% and could, therefore, potentially be spared further treatment. Conversely, patients with VHL+2 and VHL+≥3 tumors (32%) had equivalent DFS rates of 45.6% and 35.3%, respectively, and should be prioritized for adjuvant therapy.

CONCLUSIONS

Genomic characterization of ccRCC identified biologically distinct groups of patients with divergent relapse rates. These groups account for the ∼80% of cases with VHL mutations and could be used to personalize adjuvant treatment discussions with patients as well as inform future adjuvant trial design.

Single substitution in H3.3G34 alters DNMT3A recruitment to cause progressive neurodegeneration

Germline histone H3.3 amino acid substitutions, including H3.3G34R/V, cause severe neurodevelopmental syndromes. To understand how these mutations impact brain development, we generated H3.3G34R/V/W knock-in mice and identified strikingly distinct developmental defects for each mutation. H3.3G34R-mutants exhibited progressive microcephaly and neurodegeneration, with abnormal accumulation of disease-associated microglia and concurrent neuronal depletion. G34R severely decreased H3K36me2 on the mutant H3.3 tail, impairing recruitment of DNA methyltransferase DNMT3A and its redistribution on chromatin. These changes were concurrent with sustained expression of complement and other innate immune genes possibly through loss of non-CG (CH) methylation and silencing of neuronal gene promoters through aberrant CG methylation. Complement expression in G34R brains may lead to neuroinflammation possibly accounting for progressive neurodegeneration. Our study reveals that H3.3G34-substitutions have differential impact on the epigenome, which underlie the diverse phenotypes observed, and uncovers potential roles for H3K36me2 and DNMT3A-dependent CH-methylation in modulating synaptic pruning and neuroinflammation in post-natal brains.

The Histone Variant H2A.Z C-Terminal Domain Has Locus-Specific Differential Effects on H2A.Z Occupancy and Nucleosome Localization

The incorporation of histone variant H2A.Z into nucleosomes creates specialized chromatin domains that regulate DNA-templated processes, such as gene transcription. In Saccharomyces cerevisiae, the diverging H2A.Z C terminus is thought to provide the H2A.Z exclusive functions. To elucidate the roles of this H2A.Z C terminus genome-wide, we used derivatives in which the C terminus was replaced with the corresponding region of H2A (ZA protein), or the H2A region plus a transcriptional activating peptide (ZA-rII'), with the intent of regenerating the H2A.Z-dependent regulation globally. The distribution of these H2A.Z derivatives indicates that the H2A.Z C-terminal region is crucial for both maintaining the occupation level of H2A.Z and the proper positioning of targeted nucleosomes. Interestingly, the specific contribution on incorporation efficiency versus nucleosome positioning varies enormously depending on the locus analyzed. Specifically, the role of H2A.Z in global transcription regulation relies on its C-terminal region. Remarkably, however, this mostly involves genes without a H2A.Z nucleosome in the promoter. Lastly, we demonstrate that the main chaperone complex which deposits H2A.Z to gene regulatory region (SWR1-C) is necessary to localize all H2A.Z derivatives at their specific loci, indicating that the differential association of these derivatives is not due to impaired interaction with SWR1-C. IMPORTANCE We provide evidence that the Saccharomyces cerevisiae C-terminal region of histone variant H2A.Z can mediate its special function in performing gene regulation by interacting with effector proteins and chaperones. These functional interactions allow H2A.Z not only to incorporate to very specific gene regulatory regions, but also to facilitate the gene expression process. To achieve this, we used a chimeric protein which lacks the native H2A.Z C-terminal region but contains an acidic activating region, a module that is known to interact with components of chromatin-remodeling entities and/or transcription modulators. We reasoned that because this activating region can fulfill the role of the H2A.Z C-terminal region, at least in part, the role of the latter would be to interact with these activating region targets.

The mutational impact of Illudin S on human cells

Illudin S (ILS) is a fungal sesquiterpene secondary metabolite with potent genotoxic and cytotoxic properties. Early genetic studies and more recent genome-wide CRISPR screens showed that Illudin-induced lesions are preferentially repaired by transcription-coupled nucleotide excision repair (TC-NER) with some contribution from post-replication repair pathways. In line with these results, Irofulven, a semi-synthetic ILS analog was recently shown to be particularly effective on cell lines and patient-derived xenografts with impaired NER (e.g. ERCC2/3 mutations), raising hope that ILS-derived molecules may soon enter the clinic. Despite the therapeutic potential of ILS and its analogs, we still lack a global understanding of their mutagenic potential. Here, we characterize the mutational signatures associated with chronic exposure to ILS in human cells. ILS treatment rapidly stalls DNA replication and transcription, leading to the activation of the replication stress response and the accumulation of DNA damage. Novel single and double base substitution signatures as well as a characteristic indel signature indicate that ILS treatment preferentially alkylates purine residues and induces oxidative stress, confirming prior in vitro data. Many mutation contexts exhibit a strong transcriptional strand bias, highlighting the contribution of TC-NER to the repair of ILS lesions. Finally, collateral mutations are also observed in response to ILS, suggesting a contribution of translesion synthesis pathways to ILS tolerance. Accordingly, ILS treatment led to the rapid recruitment of the Y-family DNA polymerase kappa onto chromatin, supporting its preferential use for ILS lesion bypass. Altogether, our work provides the first global assessment of the genomic impact of ILS, demonstrating the contribution of multiple DNA repair pathways to ILS resistance and mutagenicity.

GATA4 controls regionalization of tissue immunity and commensal-driven immunopathology

There is growing recognition that regionalization of bacterial colonization and immunity along the intestinal tract has an important role in health and disease. Yet, the mechanisms underlying intestinal regionalization and its dysregulation in disease are not well understood. This study found that regional epithelial expression of the transcription factor GATA4 controls bacterial colonization and inflammatory tissue immunity in the proximal small intestine by regulating retinol metabolism and luminal IgA. Furthermore, in mice without jejunal GATA4 expression, the commensal segmented filamentous bacteria promoted pathogenic inflammatory immune responses that disrupted barrier function and increased mortality upon Citrobacter rodentium infection. In celiac disease patients, low GATA4 expression was associated with metabolic alterations, mucosal Actinobacillus, and increased IL-17 immunity. Taken together, these results reveal broad impacts of GATA4-regulated intestinal regionalization on bacterial colonization and tissue immunity, highlighting an elaborate interdependence of intestinal metabolism, immunity, and microbiota in homeostasis and disease.

An experimental oil spill at a tidal freshwater wetland along the St. Lawrence River re-visited after 21 years

In 1999, a tidal wetland located along the St. Lawrence River close to Ste. Croix de Lotbinière (Quebec, Eastern Canada) was the site of an experimental oil spill. Test plots were established and subjected to an experimental crude oil spill to evaluate natural attenuation, nutrient amendment and vegetation cropping as countermeasures. In 2020, this study re-visited the test plots to investigate residual oil and habitat recovery. Only concentrations of mid-chain length n-alkanes (C10-C36), but not of polycyclic aromatic hydrocarbons (PAHs), were significantly above detection limit, and were detected in both test plot and control sediments. Hydrocarbon, total organic carbon, nitrogen and phosphate contents did not differ significantly between test plot and control sediments. Microbial analyses did not detect significant differences in microbial load, microbial diversity or microbial community composition between test plot and control sediments. Key genes for the aerobic and anaerobic degradation of n-alkanes as well as for the aerobic degradation of PAHs were detected in all sediment samples. Associated gene abundances did not differ significantly between test plot and control sediments. This study shows that oil-exposed test plot sediments of the Ste. Croix wetland can be considered completely recovered after 21 years irrespective of the performed countermeasure.

Identification of R-Spondin Gene Signature Predictive of Metastatic Progression in BRAFV600E-Positive Papillary Thyroid Cancer

Papillary thyroid carcinoma (PTC) is the most common malignancy of the thyroid gland and early stages are curable. However, a subset of PTCs shows an unusually aggressive phenotype with extensive lymph node metastasis and higher incidence of locoregional recurrence. In this study, we investigated a large cohort of PTC cases with an unusual aggressive phenotype using a high-throughput RNA sequencing (RNA-Seq) to identify differentially regulated genes associated with metastatic PTC. All metastatic PTC with mutated BRAF (V600E) but not BRAF wild-type expressed an up-regulation of R-Spondin Protein 4 (RSPO4) concomitant with an upregulation of genes involved in focal adhesion and cell-extracellular matrix signaling. Further immunohistochemistry validation confirmed the upregulation of these target genes in metastatic PTC cases. Preclinical studies using established PTC cell lines support that RSPO4 overexpression is associated with BRAF V600E mutation and is a critical upstream event that promote activation of kinases of focal adhesion signaling known to drive cancer cell locomotion and invasion. This finding opens up the potential of co-targeting B-Raf, RSPO and focal adhesion proteins as a pharmacological approach for aggressive BRAF V600E PTC.

A Multi-Year, Multi-Cultivar Approach to Differential Expression Analysis of High- and Low-Protein Soybean (Glycine max)

Soybean (Glycine max (L.) Merr.) is among the most valuable crops based on its nutritious seed protein and oil. Protein quality, evaluated as the ratio of glycinin (11S) to β-conglycinin (7S), can play a role in food and feed quality. To help uncover the underlying differences between high and low protein soybean varieties, we performed differential expression analysis on high and low total protein soybean varieties and high and low 11S soybean varieties grown in four locations across Eastern and Western Canada over three years (2018-2020). Simultaneously, ten individual differential expression datasets for high vs. low total protein soybeans and ten individual differential expression datasets for high vs. low 11S soybeans were assessed, for a total of 20 datasets. The top 15 most upregulated and the 15 most downregulated genes were extracted from each differential expression dataset and cross-examination was conducted to create shortlists of the most consistently differentially expressed genes. Shortlisted genes were assessed for gene ontology to gain a global appreciation of the commonly differentially expressed genes. Genes with roles in the lipid metabolic pathway and carbohydrate metabolic pathway were differentially expressed in high total protein and high 11S soybeans in comparison to their low total protein and low 11S counterparts. Expression differences were consistent between East and West locations with the exception of one, Glyma.03G054100. These data are important for uncovering the genes and biological pathways responsible for the difference in seed protein between high and low total protein or 11S cultivars.

K27M in canonical and noncanonical H3 variants occurs in distinct oligodendroglial cell lineages in brain midline gliomas

Canonical (H3.1/H3.2) and noncanonical (H3.3) histone 3 K27M-mutant gliomas have unique spatiotemporal distributions, partner alterations and molecular profiles. The contribution of the cell of origin to these differences has been challenging to uncouple from the oncogenic reprogramming induced by the mutation. Here, we perform an integrated analysis of 116 tumors, including single-cell transcriptome and chromatin accessibility, 3D chromatin architecture and epigenomic profiles, and show that K27M-mutant gliomas faithfully maintain chromatin configuration at developmental genes consistent with anatomically distinct oligodendrocyte precursor cells (OPCs). H3.3K27M thalamic gliomas map to prosomere 2-derived lineages. In turn, H3.1K27M ACVR1-mutant pontine gliomas uniformly mirror early ventral NKX6-1⁺/SHH-dependent brainstem OPCs, whereas H3.3K27M gliomas frequently resemble dorsal PAX3⁺/BMP-dependent progenitors. Our data suggest a context-specific vulnerability in H3.1K27M-mutant SHH-dependent ventral OPCs, which rely on acquisition of ACVR1 mutations to drive aberrant BMP signaling required for oncogenesis. The unifying action of K27M mutations is to restrict H3K27me3 at PRC2 landing sites, whereas other epigenetic changes are mainly contingent on the cell of origin chromatin state and cycling rate.

High-resolution shotgun metagenomics: the more data, the better?

In shotgun metagenomics (SM), the state-of-the-art bioinformatic workflows are referred to as high-resolution shotgun metagenomics (HRSM) and require intensive computing and disk storage resources. While the increase in data output of the latest iteration of high-throughput DNA sequencing systems can allow for unprecedented sequencing depth at a minimal cost, adjustments in HRSM workflows will be needed to properly process these ever-increasing sequence datasets. One potential adaptation is to generate so-called shallow SM datasets that contain fewer sequencing data per sample as compared with the more classic high coverage sequencing. While shallow sequencing is a promising avenue for SM data analysis, detailed benchmarks using real-data are lacking. In this case study, we took four public SM datasets, one massive and the others moderate in size and subsampled each dataset at various levels to mimic shallow sequencing datasets of various sequencing depths. Our results suggest that shallow SM sequencing is a viable avenue to obtain sound results regarding microbial community structures and that high-depth sequencing does not bring additional elements for ecological interpretation. More specifically, results obtained by subsampling as little as 0.5 M sequencing clusters per sample were similar to the results obtained with the largest subsampled dataset for human gut and agricultural soil datasets. For an Antarctic dataset, which contained only a few samples, 4 M sequencing clusters per sample was found to generate comparable results to the full dataset. One area where ultra-deep sequencing and maximizing the usage of all data was undeniably beneficial was in the generation of metagenome-assembled genomes.

Dehydration response in Stylosanthes scabra: Transcriptional, biochemical, and physiological modulations

Stylosanthes scabra, popularly known as stylo, is native to the Brazilian Caatinga semiarid region and stands out as a drought-tolerant shrub forage crop. This work provides information about the plant response during the first 48 h of water deficit, followed by a rehydration treatment. Besides root transcriptomics data, 13 physiological or biochemical parameters were scrutinized. Additionally, RNA-Seq annotated transcripts not associated with the "Viridiplantae" clade were taxonomically categorized. It was found that S. scabra quickly perceives and recovers from the oscillations of the imposed water regime. Physiologically, mechanisms that minimize evapotranspiration or protect the photosynthetic apparatus stood out. Biochemically, it was found that the root tissue invests in synthesizing compounds that can act as osmolytes (proline and sugars), emphasizing the importance of osmoregulation to water deficit acclimation. Consistently, transcriptome and qPCR analyses showed that a set of enriched biological processes with upregulated (UR) transcripts were involved in protective functions against reactive oxygen species or encoding enzymes of important metabolic pathways, which might contribute to S. scabra response to water deficit. Additionally, several UR kinases and transcription factors were identified. Finally, in an innovative approach, some naturally occurring microbial groups (such as Schizosaccharomyces, Bradyrhizobium, etc.) were identified in the S. scabra roots. This study reveals insights into the physiological, biochemical, and molecular mechanisms underlying the S. scabra response to water deficit and provides candidate genes that may be useful in developing drought-tolerant crop varieties through biotechnological applications.

Liquid biopsy-based targeted gene screening highlights tumor cell subtypes in patients with advanced prostate cancer

Prostate cancer (PCa) clinical heterogeneity underscores tumor heterogeneity, which may be best defined by cell subtypes. To test if cell subtypes contributing to progression can be assessed noninvasively, we investigated whether 14 genes representing luminal, neuroendocrine, and stem cells are detectable in whole blood RNA of patients with advanced PCa. For each gene, reverse transcription quantitative polymerase chain reaction assays were first validated using RNA from PCa cell lines, and their traceability in blood was assessed in cell spiking experiments. These were next tested in blood RNA of 40 advanced PCa cases and 40 healthy controls. Expression in controls, which was low or negative, was used to define stringent thresholds for gene overexpression in patients to account for normal variation in white blood cells. Thirty-five of 40 patients overexpressed at least one gene. Patients with more genes overexpressed had a higher risk of death (hazard ratio 1.42, range 1.12-1.77). Progression on androgen receptor inhibitors was associated with overexpression of stem (odds ratio [OR] 7.74, range 1.68-35.61) and neuroendocrine (OR 13.10, range 1.24-142.34) genes, while luminal genes were associated with taxanes (OR 2.7, range 1.07-6.82). Analyses in PCa transcriptomic datasets revealed that this gene panel was most prominent in metastases of advanced disease, with diversity among patients. Collectively, these findings support the contribution of the prostate cell subtypes to disease progression. Cell-subtype specific genes are traceable in blood RNA of patients with advanced PCa and are associated with clinically relevant end points. This opens the door to minimally invasive liquid biopsies for better management of this deadly disease.

Transcriptome Analysis of Duck and Chicken Brains Infected with Aquatic Bird Bornavirus-1 (ABBV-1)

Aquatic bird bornavirus 1 (ABBV-1) is a neurotropic virus that infects waterfowls, resulting in persistent infection. Experimental infection showed that both Muscovy ducks and chickens support persistent ABBV-1 infection in the central nervous system (CNS), up to 12 weeks post-infection (wpi), without the development of clinical disease. The aim of the present study was to describe the transcriptomic profiles in the brains of experimentally infected Muscovy ducks and chickens infected with ABBV-1 at 4 and 12 wpi. Transcribed RNA was sequenced by next-generation sequencing and analyzed by principal component analysis (PCA) and differential gene expression. The functional annotation of differentially expressed genes was evaluated by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The PCA showed that the infected ducks sampled at both 4 and 12 wpi clustered separately from the controls, while only the samples from the chickens at 12 wpi, but not at 4 wpi, formed a separate cluster. In the ducks, more genes were differentially expressed at 4 wpi than 12 wpi, and the majority of the highly differentially expressed genes (DEG) were upregulated. On the other hand, the infected chickens had fewer DEGs at 4 wpi than at 12 wpi, and the majority of those with high numbers of DEGs were downregulated at 4 wpi and upregulated at 12 wpi. The functional annotation showed that the most enriched GO terms were immune-associated in both species; however, the terms associated with the innate immune response were predominantly enriched in the ducks, whereas the chickens had enrichment of both the innate and adaptive immune response. Immune-associated pathways were also enriched according to the KEGG pathway analysis in both species. Overall, the transcriptomic analysis of the duck and chicken brains showed that the main biological responses to ABBV-1 infection were immune-associated and corresponded with the levels of inflammation in the CNS.

GmSWEET29 and Paralog GmSWEET34 Are Differentially Expressed between Soybeans Grown in Eastern and Western Canada

Over the past two decades soybeans grown in western Canada have persistently had lower seed protein than those grown in eastern Canada. To understand the discrepancy in seed protein content between eastern- and western-grown soybeans, RNA-seq and differential expression analysis have been investigated. Ten soybean genotypes, ranging from low to high in seed protein content, were grown in four locations across eastern (Ottawa) and western (Morden, Brandon, and Saskatoon) Canada. Differential expression analysis revealed 34 differentially expressed genes encoding Glycine max Sugars Will Eventually be Exported Transporters (GmSWEETs), including paralogs GmSWEET29 and GmSWEET34 (AtSWEET2 homologs) that were consistently upregulated across all ten genotypes in each of the western locations over three years. GmSWEET29 and GmSWEET34 are likely candidates underlying the lower seed protein content of western soybeans. GmSWEET20 (AtSWEET12 homolog) was downregulated in the western locations and may also play a role in lower seed protein content. These findings are valuable for improving soybean agriculture in western growing regions, establishing more strategic and efficient agricultural practices.

NF-κB p65 Attenuates Cardiomyocyte PGC-1α Expression in Hypoxia

Hypoxia exerts broad effects on cardiomyocyte function and viability, ranging from altered metabolism and mitochondrial physiology to apoptotic or necrotic cell death. The transcriptional coactivator peroxisome proliferator-activated receptor gamma coactivator 1-alpha (PGC-1α) is a key regulator of cardiomyocyte metabolism and mitochondrial function and is down-regulated in hypoxia; however, the underlying mechanism is incompletely resolved. Using primary rat cardiomyocytes coupled with electrophoretic mobility shift and luciferase assays, we report that hypoxia impaired mitochondrial energetics and resulted in an increase in nuclear localization of the Nuclear Factor-κB (NF-κB) p65 subunit, and the association of p65 with the PGC-1α proximal promoter. Tumor necrosis factor α (TNFα), an activator of NF-κB signaling, similarly reduced PGC-1α expression and p65 binding to the PGC-1α promoter in a dose-dependent manner, and TNFα-mediated down-regulation of PGC-1α expression could be reversed by the NF-κB inhibitor parthenolide. RNA-seq analysis revealed that cardiomyocytes isolated from p65 knockout mice exhibited alterations in genes associated with chromatin remodeling. Decreased PGC-1α promoter transactivation by p65 could be partially reversed by the histone deacetylase inhibitor trichostatin A. These results implicate NF-κB signaling, and specifically p65, as a potent inhibitor of PGC-1α expression in cardiac myocyte hypoxia.

iCOMIC: a graphical interface-driven bioinformatics pipeline for analyzing cancer omics data

Despite the tremendous increase in omics data generated by modern sequencing technologies, their analysis can be tricky and often requires substantial expertise in bioinformatics. To address this concern, we have developed a user-friendly pipeline to analyze (cancer) genomic data that takes in raw sequencing data (FASTQ format) as input and outputs insightful statistics. Our iCOMIC toolkit pipeline featuring many independent workflows is embedded in the popular Snakemake workflow management system. It can analyze whole-genome and transcriptome data and is characterized by a user-friendly GUI that offers several advantages, including minimal execution steps and eliminating the need for complex command-line arguments. Notably, we have integrated algorithms developed in-house to predict pathogenicity among cancer-causing mutations and differentiate between tumor suppressor genes and oncogenes from somatic mutation data. We benchmarked our tool against Genome In A Bottle benchmark dataset (NA12878) and got the highest F1 score of 0.971 and 0.988 for indels and SNPs, respectively, using the BWA MEM—GATK HC DNA-Seq pipeline. Similarly, we achieved a correlation coefficient of r = 0.85 using the HISAT2-StringTie-ballgown and STAR-StringTie-ballgown RNA-Seq pipelines on the human monocyte dataset (SRP082682). Overall, our tool enables easy analyses of omics datasets, significantly ameliorating complex data analysis pipelines.

ARMC5 is part of an RPB1-specific ubiquitin ligase implicated in adrenal hyperplasia

ARMC5 is implicated in several pathological conditions, but its function remains unknown. We have previously identified CUL3 and RPB1 (the largest subunit of RNA polymerase II (Pol II) as potential ARMC5-interacting proteins. Here, we show that ARMC5, CUL3 and RBX1 form an active E3 ligase complex specific for RPB1. ARMC5, CUL3, and RBX1 formed an active E3 specific for RPB1. Armc5 deletion caused a significant reduction in RPB1 ubiquitination and an increase in an accumulation of RPB1, and hence an enlarged Pol II pool in normal tissues and organs. The compromised RPB1 degradation did not cause generalized Pol II stalling nor depressed transcription in the adrenal glands but did result in dysregulation of a subset of genes, with most upregulated. We found RPB1 to be highly expressed in the adrenal nodules from patients with primary bilateral macronodular adrenal hyperplasia (PBMAH) harboring germline ARMC5 mutations. Mutant ARMC5 had altered binding with RPB1. In summary, we discovered that wildtype ARMC5 was part of a novel RPB1-specific E3. ARMC5 mutations resulted in an enlarged Pol II pool, which dysregulated a subset of effector genes. Such an enlarged Pol II pool and gene dysregulation was correlated to adrenal hyperplasia in humans and KO mice.