Genomic Reconstruction of the Successful Establishment of a Feralized Bovine Population on the Subantarctic Island of Amsterdam.

The feral cattle of the subantarctic island of Amsterdam provide an outstanding case study of a large mammalian population that was established by a handful of founders and thrived within a few generations in a seemingly inhospitable environment. Here, we investigated the genetic history and composition of this population using genotyping and sequencing data. Our inference showed an intense but brief founding bottleneck around the late 19th century and revealed contributions from European taurine and Indian Ocean Zebu in the founder ancestry. Comparative analysis of whole-genome sequences further revealed a moderate reduction in genetic diversity despite high levels of inbreeding. The brief and intense bottleneck was associated with high levels of drift, a flattening of the site frequency spectrum and a slight relaxation of purifying selection on mildly deleterious variants. Unlike some populations that have experienced prolonged reductions in effective population size, we did not observe any significant purging of highly deleterious variants. Interestingly, the population's success in the harsh environment can be attributed to preadaptation from their European taurine ancestry, suggesting no strong bioclimatic challenge, and also contradicting evidence for insular dwarfism. Genome scan for footprints of selection uncovered a majority of candidate genes related to nervous system function, likely reflecting rapid feralization driven by behavioral changes and complex social restructuring. The Amsterdam Island cattle offers valuable insights into rapid population establishment, feralization, and genetic adaptation in challenging environments. It also sheds light on the unique genetic legacies of feral populations, raising ethical questions according to conservation efforts.

DeSUMOylation of chromatin-bound proteins limits the rapid transcriptional reprogramming induced by daunorubicin in acute myeloid leukemias.

Genotoxicants have been used for decades as front-line therapies against cancer on the basis of their DNA-damaging actions. However, some of their non-DNA-damaging effects are also instrumental for killing dividing cells. We report here that the anthracycline Daunorubicin (DNR), one of the main drugs used to treat Acute Myeloid Leukemia (AML), induces rapid (3 h) and broad transcriptional changes in AML cells. The regulated genes are particularly enriched in genes controlling cell proliferation and death, as well as inflammation and immunity. These transcriptional changes are preceded by DNR-dependent deSUMOylation of chromatin proteins, in particular at active promoters and enhancers. Surprisingly, inhibition of SUMOylation with ML-792 (SUMO E1 inhibitor), dampens DNR-induced transcriptional reprogramming. Quantitative proteomics shows that the proteins deSUMOylated in response to DNR are mostly transcription factors, transcriptional co-regulators and chromatin organizers. Among them, the CCCTC-binding factor CTCF is highly enriched at SUMO-binding sites found in cis-regulatory regions. This is notably the case at the promoter of the DNR-induced NFKB2 gene. DNR leads to a reconfiguration of chromatin loops engaging CTCF- and SUMO-bound NFKB2 promoter with a distal cis-regulatory region and inhibition of SUMOylation with ML-792 prevents these changes.

The NADPH oxidase NOX2 is a marker of adverse prognosis involved in chemoresistance of acute myeloid leukemias.

Resistance to chemotherapeutic drugs is a major cause of treatment failure in acute myeloid leukemias (AML). To better characterize the mechanisms of chemoresistance, we first identified genes whose expression is dysregulated in AML cells resistant to daunorubicin or cytarabine, the main drugs used for induction therapy. The genes found to be activated are mostly linked to immune signaling and inflammation. Among them, we identified a strong upregulation of the NOX2 NAPDH oxidase subunit genes (CYBB, CYBA, NCF1, NCF2, NCF4 and RAC2). The ensuing increase in NADPH oxidase expression and production of reactive oxygen species, which is particularly strong in daunorubicin-resistant cells, participates in the acquisition and/or maintenance of resistance to daunorubicin. Gp91phox (CYBB-encoded Nox2 catalytic subunit), was found to be more expressed and active in leukemic cells from patients with the French-American-British (FAB) M4/M5 subtypes of AML than in those from patients with the FAB M0-M2 ones. Moreover, its expression was increased at the surface of patients' chemotherapy-resistant AML cells. Finally, using a gene expression based score we demonstrated that high expression of NOX2 subunit genes is a marker of adverse prognosis in AML patients. The prognostic NOX score we defined is independent of the cytogenetic-based risk classification, FAB subtype, FLT3/NPM1 mutational status and age.

Family-effects in the epigenomic response of red blood cells to a challenge test in the European sea bass (Dicentrarchus labrax, L.).

BACKGROUND: In fish, minimally invasive blood sampling is widely used to monitor physiological stress with blood plasma biomarkers. As fish blood cells are nucleated, they might be a source a potential new markers derived from 'omics technologies. We modified the epiGBS (epiGenotyping By Sequencing) technique to explore changes in genome-wide cytosine methylation in the red blood cells (RBCs) of challenged European sea bass (Dicentrarchus labrax), a species widely studied in both natural and farmed environments. RESULTS: We retrieved 501,108,033 sequencing reads after trimming, with a mean mapping efficiency of 73.0% (unique best hits). Minor changes in RBC methylome appeared to manifest after the challenge test and a family-effect was detected. Only fifty-seven differentially methylated cytosines (DMCs) close to 51 distinct genes distributed on 17 of 24 linkage groups (LGs) were detected between RBCs of pre- and post-challenge individuals. Thirty-seven of these genes were previously reported as differentially expressed in the brain of zebrafish, most of them involved in stress coping differences. While further investigation remains necessary, few DMC-related genes associated to the Brain Derived Neurotrophic Factor, a protein that favors stress adaptation and fear memory, appear relevant to integrate a centrally produced stress response in RBCs. CONCLUSION: Our modified epiGBS protocol was powerful to analyze patterns of cytosine methylation in RBCs of D. labrax and to evaluate the impact of a challenge using minimally invasive blood samples. This study is the first approximation to identify epigenetic biomarkers of exposure to stress in fish.

Variant calling and genotyping accuracy of ddRAD-seq: Comparison with 20X WGS in layers.

Whole Genome Sequencing (WGS) remains a costly or unsuitable method for routine genotyping of laying hens. Until now, breeding companies have been using or developing SNP chips. Nevertheless, alternatives methods based on sequencing have been developed. Among these, reduced representation sequencing approaches can offer sequencing quality and cost-effectiveness by reducing the genomic regions covered by sequencing. The aim of this study was to evaluate the ability of double digested Restriction site Associated DNA sequencing (ddRAD-seq) to identify and genotype SNPs in laying hens, by comparison with a presumed reliable WGS approach. Firstly, the sensitivity and precision of variant calling and the genotyping reliability of ddRADseq were determined. Next, the SNP Call Rate (CRSNP) and mean depth of sequencing per SNP (DPSNP) were compared between both methods. Finally, the effect of multiple combinations of thresholds for these parameters on genotyping reliability and amount of remaining SNPs in ddRAD-seq was studied. In raw form, the ddRAD-seq identified 349,497 SNPs evenly distributed on the genome with a CRSNP of 0.55, a DPSNP of 11X and a mean genotyping reliability rate per SNP of 80%. Considering genomic regions covered by expected enzymatic fragments (EFs), the sensitivity of the ddRAD-seq was estimated at 32.4% and its precision at 96.4%. The low CRSNP and DPSNP values were explained by the detection of SNPs outside the EFs theoretically generated by the ddRAD-seq protocol. Indeed, SNPs outside the EFs had significantly lower CRSNP (0.25) and DPSNP (1X) values than SNPs within the EFs (0.7 and 17X, resp.). The study demonstrated the relationship between CRSNP, DPSNP, genotyping reliability and the number of SNPs retained, to provide a decision-support tool for defining filtration thresholds. Severe quality control over ddRAD-seq data allowed to retain a minimum of 40% of the SNPs with a CcR of 98%. Then, ddRAD-seq was defined as a suitable method for variant calling and genotyping in layers.

Deep population structure linked to host vernalization requirement in the barley net blotch fungal pathogen.

Invasive fungal pathogens pose a substantial threat to widely cultivated crop species, owing to their capacity to adapt to new hosts and new environmental conditions. Gaining insights into the demographic history of these pathogens and unravelling the mechanisms driving coevolutionary processes are crucial for developing durably effective disease management programmes. Pyrenophora teres is a significant fungal pathogen of barley, consisting of two lineages, Ptt and Ptm, with global distributions and demographic histories reflecting barley domestication and spread. However, the factors influencing the population structure of P. teres remain poorly understood, despite the varietal and environmental heterogeneity of barley agrosystems. Here, we report on the population genomic structure of P. teres in France and globally. We used genotyping-by-sequencing to show that Ptt and Ptm can coexist in the same area in France, with Ptt predominating. Furthermore, we showed that differences in the vernalization requirement of barley varieties were associated with population differentiation within Ptt in France and at a global scale, with one population cluster found on spring barley and another population cluster found on winter barley. Our results demonstrate how cultivation conditions, possibly associated with genetic differences between host populations, can be associated with the maintenance of divergent invasive pathogen populations coexisting over large geographic areas. This study not only advances our understanding of the coevolutionary dynamics of the Pt-barley pathosystem but also prompts further research on the relative contributions of adaptation to the host versus adaptation to abiotic conditions in shaping Ptt populations.

Draft genome sequences of Neurospora crassa clade B, isolated from burned Cytisus sp. plants in France.

Neurospora crassa clade A is a model system for genetics, biochemistry, molecular biology, and experimental evolution. Here, we present the draft genome sequences of four isolates of N. crassa clade B. These data represent a valuable resource to investigate the population biology and evolutionary history of N. crassa sensu lato.

Characterizing intergenic transcription at RNA polymerase II binding sites in normal and cancer tissues.

Intergenic transcription in normal and cancerous tissues is pervasive but incompletely understood. To investigate this, we constructed an atlas of over 180,000 consensus RNA polymerase II (RNAPII)-bound intergenic regions from 900 RNAPII chromatin immunoprecipitation sequencing (ChIP-seq) experiments in normal and cancer samples. Through unsupervised analysis, we identified 51 RNAPII consensus clusters, many of which mapped to specific biotypes and revealed tissue-specific regulatory signatures. We developed a meta-clustering methodology to integrate our RNAPII atlas with active transcription across 28,797 RNA sequencing (RNA-seq) samples from The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), and Encyclopedia of DNA Elements (ENCODE). This analysis revealed strong tissue- and disease-specific interconnections between RNAPII occupancy and transcriptional activity. We demonstrate that intergenic transcription at RNAPII-bound regions is a novel per-cancer and pan-cancer biomarker. This biomarker displays genomic and clinically relevant characteristics, distinguishing cancer subtypes and linking to overall survival. Our results demonstrate the effectiveness of coherent data integration to uncover intergenic transcriptional activity in normal and cancer tissues.

In vivo investigation of Lcr35® anti-candidiasis properties in Caenorhabditis elegans reveals the involvement of highly conserved immune pathways.

Lactic acid bacteria, including the microorganisms formerly designated as Lactobacillus, are the major representatives of Live Biotherapeutic Microorganisms (LBM) when used for therapeutic purposes. However, in most cases, the mechanisms of action remain unknown. The antifungal potential of LBM has already been demonstrated using preclinical models (cell cultures, laboratory animals). Understanding their mechanisms of action is strategic for the development of new therapeutics for humans. Here, Caenorhabditis elegans was used as an in vivo model to analyze pro-longevity, anti-aging and anti-candidiasis effects of the LBM Lacticaseibacillus rhamnosus (formerly Lactobacillus rhamnosus) Lcr35®. A high-throughput transcriptomic analysis revealed a specific response of C. elegans depending on whether it is in the presence of the LBM L. rhamnosus Lcr35® (structural response), the yeast Candida albicans (metabolic response) or both (structural and metabolic responses) in a preventive and a curative conditions. Studies on C. elegans mutants demonstrated that the p38 MAPK (sek-1, skn-1) and the insulin-like (daf-2, daf-16) signaling pathways were involved in the extended lifespan provided by L. rhamnosus Lcr35® strain whereas the JNK pathway was not involved (jnk-1). In addition, the anti C. albicans effect of the bacterium requires the daf-16 and sek-1 genes while it is independent of daf-2 and skn-1. Moreover, the anti-aging effect of Lcr35®, linked to the extension of longevity, is not due to protection against oxidative stress (H2O2). Taken together, these results formally show the involvement of the p38 MAP kinase and insulin-like signaling pathways for the longevity extension and anti-Candida albicans properties of Lcr35® with, however, differences in the genes involved. Overall, these findings provide new insight for understanding the mechanisms of action of a probiotic strain with antimicrobial potential.

MEOX2 Transcription Factor Is Involved in Survival and Adhesion of Glioma Stem-like Cells.

The high expression of MEOX2 transcription factor is closely associated with poor overall survival in glioma. MEOX2 has recently been described as an interesting prognostic biomarker, especially for lower grade glioma. MEOX2 has never been studied in glioma stem-like cells (GSC), responsible for glioma recurrence. The aim of our study was to investigate the role of MEOX2 in GSC. Loss of function approach using siRNA was used to assess the impact of MEOX2 on GSC viability and stemness phenotype. MEOX2 was localized in the nucleus and its expression was heterogeneous between GSCs. MEOX2 expression depends on the methylation state of its promoter and is strongly associated with IDH mutations. MEOX2 is involved in cell proliferation and viability regulation through ERK/MAPK and PI3K/AKT pathways. MEOX2 loss of function correlated with GSC differentiation and acquisition of neuronal lineage characteristics. Besides, inhibition of MEOX2 is correlated with increased expression of CDH10 and decreased pFAK. In this study, we unraveled, for the first time, MEOX2 contribution to cell viability and proliferation through AKT/ERK pathway and its potential involvement in phenotype and adhesion properties of GSC.

Concurrent BMP Signaling Maintenance and TGF-ß Signaling Inhibition Is a Hallmark of Natural Resistance to Muscle Atrophy in the Hibernating Bear.

Muscle atrophy arises from a multiplicity of physio-pathological situations and has very detrimental consequences for the whole body. Although knowledge of muscle atrophy mechanisms keeps growing, there is still no proven treatment to date. This study aimed at identifying new drivers for muscle atrophy resistance. We selected an innovative approach that compares muscle transcriptome between an original model of natural resistance to muscle atrophy, the hibernating brown bear, and a classical model of induced atrophy, the unloaded mouse. Using RNA sequencing, we identified 4415 differentially expressed genes, including 1746 up- and 2369 down-regulated genes, in bear muscles between the active versus hibernating period. We focused on the Transforming Growth Factor (TGF)-ß and the Bone Morphogenetic Protein (BMP) pathways, respectively, involved in muscle mass loss and maintenance. TGF-ß- and BMP-related genes were overall down- and up-regulated in the non-atrophied muscles of the hibernating bear, respectively, and the opposite occurred for the atrophied muscles of the unloaded mouse. This was further substantiated at the protein level. Our data suggest TGF-ß/BMP balance is crucial for muscle mass maintenance during long-term physical inactivity in the hibernating bear. Thus, concurrent activation of the BMP pathway may potentiate TGF-ß inhibiting therapies already targeted to prevent muscle atrophy.