Complete and circularized genome sequences of five nitrogen-fixing Bradyrhizobium sp. strains isolated from root nodules of peanut, Arachis hypogaea, cultivated in Tunisia.

This manuscript reports the complete and circularized Oxford Nanopore Technologies (ONT) long read-based genome sequences of five nitrogen-fixing symbionts belonging to the genus Bradyrhizobium, isolated from root nodules of peanut (Arachis hypogaea) grown on soil samples collected from Tunisia.

Faecalibacterium duncaniae A2-165 regulates the expression of butyrate synthesis, ferrous iron uptake, and stress-response genes based on acetate consumption.

The promising next-generation probiotic Faecalibacterium prausnitzii is one of the most abundant acetate-consuming, butyrate-producing bacteria in the healthy human gut. Yet, little is known about how acetate availability affects this bacterium's gene expression strategies. Here, we investigated the effect of acetate on temporal changes in the transcriptome of F. duncaniae A2-165 cultures using RNA sequencing. We compared gene expression patterns between two growth phases (early stationary vs. late exponential) and two acetate levels (low: 3 mM vs. high: 23 mM). Only in low-acetate conditions, a general stress response was activated. In high-acetate conditions, there was greater expression of genes related to butyrate synthesis and to the importation of B vitamins and iron. Specifically, expression was strongly activated in the case of the feoAABC operon, which encodes a FeoB ferrous iron transporter, but not in the case of the feoAB gene, which encodes a second putative FeoAB transporter. Moreover, excess ferrous iron repressed feoB expression but not feoAB. Lastly, FeoB but not FeoAB peptides from strain A2-165 were found in abundance in a healthy human fecal metaproteome. In conclusion, we characterized two early-stationary transcriptomes based on acetate consumption and this work highlights the regulation of feoB expression in F. duncaniae A2-165.

Genomics in the long-read sequencing era.

Long-read sequencing (LRS) technologies have provided extremely powerful tools to explore genomes. While in the early years these methods suffered technical limitations, they have recently made significant progress in terms of read length, throughput, and accuracy and bioinformatics tools have strongly improved. Here, we aim to review the current status of LRS technologies, the development of novel methods, and the impact on genomics research. We will explore the most impactful recent findings made possible by these technologies focusing on high-resolution sequencing of genomes and transcriptomes and the direct detection of DNA and RNA modifications. We will also discuss how LRS methods promise a more comprehensive understanding of human genetic variation, transcriptomics, and epigenetics for the coming years.

Mistimed origin licensing and activation stabilize common fragile sites under tight DNA-replication checkpoint activation.

Genome integrity requires replication to be completed before chromosome segregation. The DNA-replication checkpoint (DRC) contributes to this coordination by inhibiting CDK1, which delays mitotic onset. Under-replication of common fragile sites (CFSs), however, escapes surveillance, resulting in mitotic chromosome breaks. Here we asked whether loose DRC activation induced by modest stresses commonly used to destabilize CFSs could explain this leakage. We found that tightening DRC activation or CDK1 inhibition stabilizes CFSs in human cells. Repli-Seq and molecular combing analyses showed a burst of replication initiations implemented in mid S-phase across a subset of late-replicating sequences, including CFSs, while the bulk genome was unaffected. CFS rescue and extra-initiations required CDC6 and CDT1 availability in S-phase, implying that CDK1 inhibition permits mistimed origin licensing and firing. In addition to delaying mitotic onset, tight DRC activation therefore supports replication completion of late origin-poor domains at risk of under-replication, two complementary roles preserving genome stability.

The Drosophila Fab-7 boundary modulates Abd-B gene activity by guiding an inversion of collinear chromatin organization and alternate promoter use.

Hox genes encode transcription factors that specify segmental identities along the anteroposterior body axis. These genes are organized in clusters, where their order corresponds to their activity along the body axis, a feature known as collinearity. In Drosophila, the BX-C cluster contains the three most posterior Hox genes, where their collinear activation incorporates progressive changes in histone modifications, chromatin architecture, and use of boundary elements and cis-regulatory regions. To dissect functional hierarchies, we compare chromatin organization in cell lines and larvae, with a focus on the Abd-B gene. Our work establishes the importance of the Fab-7 boundary for insulation between 3D domains carrying different histone modifications. Interestingly, we detect a non-canonical inversion of collinear chromatin dynamics at Abd-B, with the domain of active histone modifications progressively decreasing in size. This dynamic chromatin organization differentially activates the alternative promoters of the Abd-B gene, thereby expanding the possibilities for fine-tuning of transcriptional output.

The mineralocorticoid receptor modulates timing and location of genomic binding by glucocorticoid receptor in response to synthetic glucocorticoids in keratinocytes.

Glucocorticoids (GCs) exert potent antiproliferative and anti-inflammatory properties, explaining their therapeutic efficacy for skin diseases. GCs act by binding to the GC receptor (GR) and the mineralocorticoid receptor (MR), co-expressed in classical and non-classical targets including keratinocytes. Using knockout mice, we previously demonstrated that GR and MR exert essential nonoverlapping functions in skin homeostasis. These closely related receptors may homo- or heterodimerize to regulate transcription, and theoretically bind identical GC-response elements (GRE). We assessed the contribution of MR to GR genomic binding and the transcriptional response to the synthetic GC dexamethasone (Dex) using control (CO) and MR knockout (MREKO ) keratinocytes. GR chromatin immunoprecipitation (ChIP)-seq identified peaks common and unique to both genotypes upon Dex treatment (1 h). GREs, AP-1, TEAD, and p53 motifs were enriched in CO and MREKO peaks. However, GR genomic binding was 35% reduced in MREKO , with significantly decreased GRE enrichment, and reduced nuclear GR. Surface plasmon resonance determined steady state affinity constants, suggesting preferred dimer formation as MR-MR > GR-MR ~ GR-GR; however, kinetic studies demonstrated that GR-containing dimers had the longest lifetimes. Despite GR-binding differences, RNA-seq identified largely similar subsets of differentially expressed genes in both genotypes upon Dex treatment (3 h). However, time-course experiments showed gene-dependent differences in the magnitude of expression, which correlated with earlier and more pronounced GR binding to GRE sites unique to CO including near Nr3c1. Our data show that endogenous MR has an impact on the kinetics and differential genomic binding of GR, affecting the time-course, specificity, and magnitude of GC transcriptional responses in keratinocytes.

DNA methylation and hydroxymethylation characterize the identity of D1 and D2 striatal projection neurons.

Neuronal DNA modifications differ from those in other cells, including methylation outside CpG context and abundant 5-hydroxymethylation whose relevance for neuronal identities are unclear. Striatal projection neurons expressing D1 or D2 dopamine receptors allow addressing this question, as they share many characteristics but differ in their gene expression profiles, connections, and functional roles. We compare translating mRNAs and DNA modifications in these two populations. DNA methylation differences occur predominantly in large genomic clusters including differentially expressed genes, potentially important for D1 and D2 neurons. Decreased gene body methylation is associated with higher gene expression. Hydroxymethylation differences are more scattered and affect transcription factor binding sites, which can influence gene expression. We also find a strong genome-wide hydroxymethylation asymmetry between the two DNA strands, particularly pronounced at expressed genes and retrotransposons. These results identify novel properties of neuronal DNA modifications and unveil epigenetic characteristics of striatal projection neurons heterogeneity.

Identification of the Inner Cell Mass and the Trophectoderm Responses after an In Vitro Exposure to Glucose and Insulin during the Preimplantation Period in the Rabbit Embryo.

The prevalence of metabolic diseases is increasing, leading to more women entering pregnancy with alterations in the glucose-insulin axis. The aim of this work was to investigate the effect of a hyperglycemic and/or hyperinsulinemic environment on the development of the preimplantation embryo. In rabbit embryos developed in vitro in the presence of high insulin (HI), high glucose (HG), or both (HGI), we determined the transcriptomes of the inner cell mass (ICM) and the trophectoderm (TE). HI induced 10 differentially expressed genes (DEG) in ICM and 1 in TE. HG ICM exhibited 41 DEGs involved in oxidative phosphorylation (OXPHOS) and cell number regulation. In HG ICM, proliferation was decreased (p < 0.01) and apoptosis increased (p < 0.001). HG TE displayed 132 DEG linked to mTOR signaling and regulation of cell number. In HG TE, proliferation was increased (p < 0.001) and apoptosis decreased (p < 0.001). HGI ICM presented 39 DEG involved in OXPHOS and no differences in proliferation and apoptosis. HGI TE showed 16 DEG linked to OXPHOS and cell number regulation and exhibited increased proliferation (p < 0.001). Exposure to HG and HGI during preimplantation development results in common and specific ICM and TE responses that could compromise the development of the future individual and placenta.

Developmental timing of programmed DNA elimination in Paramecium tetraurelia recapitulates germline transposon evolutionary dynamics.

With its nuclear dualism, the ciliate Paramecium constitutes a unique model to study how host genomes cope with transposable elements (TEs). P. tetraurelia harbors two germline micronuclei (MICs) and a polyploid somatic macronucleus (MAC) that develops from one MIC at each sexual cycle. Throughout evolution, the MIC genome has been continuously colonized by TEs and related sequences that are removed from the somatic genome during MAC development. Whereas TE elimination is generally imprecise, excision of approximately 45,000 TE-derived internal eliminated sequences (IESs) is precise, allowing for functional gene assembly. Programmed DNA elimination is concomitant with genome amplification. It is guided by noncoding RNAs and repressive chromatin marks. A subset of IESs is excised independently of this epigenetic control, raising the question of how IESs are targeted for elimination. To gain insight into the determinants of IES excision, we established the developmental timing of DNA elimination genome-wide by combining fluorescence-assisted nuclear sorting with high-throughput sequencing. Essentially all IESs are excised within only one endoreplication round (32C to 64C), whereas TEs are eliminated at a later stage. We show that DNA elimination proceeds independently of replication. We defined four IES classes according to excision timing. The earliest excised IESs tend to be independent of epigenetic factors, display strong sequence signals at their ends, and originate from the most ancient integration events. We conclude that old IESs have been optimized during evolution for early and accurate excision by acquiring stronger sequence determinants and escaping epigenetic control.

Major transcriptomic, epigenetic and metabolic changes underlie the pluripotency continuum in rabbit preimplantation embryos.

Despite the growing interest in the rabbit model for developmental and stem cell biology, the characterization of embryos at the molecular level is still poorly documented. We conducted a transcriptome analysis of rabbit preimplantation embryos from E2.7 (morula stage) to E6.6 (early primitive streak stage) using bulk and single-cell RNA-sequencing. In parallel, we studied oxidative phosphorylation and glycolysis, and analysed active and repressive epigenetic modifications during blastocyst formation and expansion. We generated a transcriptomic, epigenetic and metabolic map of the pluripotency continuum in rabbit preimplantation embryos, and identified novel markers of naive pluripotency that might be instrumental for deriving naive pluripotent stem cell lines. Although the rabbit is evolutionarily closer to mice than to primates, we found that the transcriptome of rabbit epiblast cells shares common features with those of humans and non-human primates.

Molecular and Cellular Analysis of the Repair of Zebrafish Optic Tectum Meninges Following Laser Injury.

We studied cell recruitment following optic tectum (OT) injury in zebrafish (Danio rerio), which has a remarkable ability to regenerate many of its organs, including the brain. The OT is the largest dorsal layered structure in the zebrafish brain. In juveniles, it is an ideal structure for imaging and dissection. We investigated the recruited cells within the juvenile OT during regeneration in a Pdgfrß-Gal4:UAS-EGFP line in which pericytes, vascular, circulating, and meningeal cells are labeled, together with neurons and progenitors. We first performed high-resolution confocal microscopy and single-cell RNA-sequencing (scRNAseq) on EGFP-positive cells. We then tested three types of injury with very different outcomes (needle (mean depth in the OT of 200 µm); deep-laser (depth: 100 to 200 µm depth); surface-laser (depth: 0 to 100 µm)). Laser had the additional advantage of better mimicking of ischemic cerebral accidents. No massive recruitment of EGFP-positive cells was observed following laser injury deep in the OT. This type of injury does not perturb the meninx/brain-blood barrier (BBB). We also performed laser injuries at the surface of the OT, which in contrast create a breach in the meninges. Surprisingly, one day after such injury, we observed the migration to the injury site of various EGFP-positive cell types at the surface of the OT. The migrating cells included midline roof cells, which activated the PI3K-AKT pathway; fibroblast-like cells expressing numerous collagen genes and most prominently in 3D imaging; and a large number of arachnoid cells that probably migrate to the injury site through the activation of cilia motility genes, most likely being direct targets of the FOXJ1a gene. This study, combining high-content imaging and scRNAseq in physiological and pathological conditions, sheds light on meninges repair mechanisms in zebrafish that probably also operate in mammalian meninges.

Maternal age affects equine day 8 embryo gene expression both in trophoblast and inner cell mass.

BACKGROUND: Breeding a mare until she is not fertile or even until her death is common in equine industry but the fertility decreases as the mare age increases. Embryo loss due to reduced embryo quality is partly accountable for this observation. Here, the effect of mare's age on blastocysts' gene expression was explored. Day 8 post-ovulation embryos were collected from multiparous young (YM, 6-year-old, N = 5) and older (OM, > 10-year-old, N = 6) non-nursing Saddlebred mares, inseminated with the semen of one stallion. Pure or inner cell mass (ICM) enriched trophoblast, obtained by embryo bisection, were RNA sequenced. Deconvolution algorithm was used to discriminate gene expression in the ICM from that in the trophoblast. Differential expression was analyzed with embryo sex and diameter as cofactors. Functional annotation and classification of differentially expressed genes and gene set enrichment analysis were also performed. RESULTS: Maternal aging did not affect embryo recovery rate, embryo diameter nor total RNA quantity. In both compartments, the expression of genes involved in mitochondria and protein metabolism were disturbed by maternal age, although more genes were affected in the ICM. Mitosis, signaling and adhesion pathways and embryo development were decreased in the ICM of embryos from old mares. In trophoblast, ion movement pathways were affected. CONCLUSIONS: This is the first study showing that maternal age affects gene expression in the equine blastocyst, demonstrating significant effects as early as 10 years of age. These perturbations may affect further embryo development and contribute to decreased fertility due to aging.

Paramecium Polycomb repressive complex 2 physically interacts with the small RNA-binding PIWI protein to repress transposable elements.

Polycomb repressive complex 2 (PRC2) maintains transcriptionally silent genes in a repressed state via deposition of histone H3K27-trimethyl (me3) marks. PRC2 has also been implicated in silencing transposable elements (TEs), yet how PRC2 is targeted to TEs remains unclear. To address this question, we identified proteins that physically interact with the Paramecium enhancer-of-zeste Ezl1 enzyme, which catalyzes H3K9me3 and H3K27me3 deposition at TEs. We show that the Paramecium PRC2 core complex comprises four subunits, each required in vivo for catalytic activity. We also identify PRC2 cofactors, including the RNA interference (RNAi) effector Ptiwi09, which are necessary to target H3K9me3 and H3K27me3 to TEs. We find that the physical interaction between PRC2 and the RNAi pathway is mediated by a RING finger protein and that small RNA recruitment of PRC2 to TEs is analogous to the small RNA recruitment of H3K9 methylation SU(VAR)3-9 enzymes.

Dynamics of the compartmentalized Streptomyces chromosome during metabolic differentiation.

Bacteria of the genus Streptomyces are prolific producers of specialized metabolites, including antibiotics. The linear chromosome includes a central region harboring core genes, as well as extremities enriched in specialized metabolite biosynthetic gene clusters. Here, we show that chromosome structure in Streptomyces ambofaciens correlates with genetic compartmentalization during exponential phase. Conserved, large and highly transcribed genes form boundaries that segment the central part of the chromosome into domains, whereas the terminal ends tend to be transcriptionally quiescent compartments with different structural features. The onset of metabolic differentiation is accompanied by a rearrangement of chromosome architecture, from a rather 'open' to a 'closed' conformation, in which highly expressed specialized metabolite biosynthetic genes form new boundaries. Thus, our results indicate that the linear chromosome of S. ambofaciens is partitioned into structurally distinct entities, suggesting a link between chromosome folding, gene expression and genome evolution.

Data from crosslinking and analysis of cDNAs (CRAC) of Nab3 in yeast cells expressing a circular ncRNA decoy.

Pervasive transcription originating from the ubiquitous activity of RNA Polymerase II (RNAPII) generates a vast mass of non-coding RNAs (ncRNAs) that represent a potential harm to gene expression. In the compact genome of the yeast Saccharomyces cerevisiae, the main genomewide safeguard against pervasive ncRNAs is the Nrd1-Nab3-Sen1 (NNS) complex, composed of two RNA-binding proteins (Nrd1 and Nab3) and the helicase Sen1. The NNS complex directs transcription termination of ncRNA genes and promotes the rapid degradation of pervasive transcripts from yeast nuclei through its physical and functional coupling to the nuclear RNA exosome. We have recently shown that inhibition of the exosome in yeast cells leads to the accumulation of ncRNAs complexed with Nab3 and Nrd1, decreasing recycling of these termination factors to sites of transcription and inducing global termination defects at NNS targets. Consistent with the notion that ncRNAs out-titrate Nab3 and Nrd1 termination factors, we have shown that a similar genomewide termination impairment could be achieved by expressing a circular RNA decoy containing a Nab3 binding target [1]. In relation to this previous research article, here we expand our observations on the effect of the circular RNA decoy on NNS termination. We aimed at verifying that the Nab3 binding sequence present on the decoy is indeed efficiently sequestering Nab3 as intended by design, leading to the expected decrease of Nab3 binding on NNS targets. We employed the crosslinking and cDNA analysis protocol (CRAC) on yeast cells expressing the circular ncRNA decoy or a control construct. We present data from high-resolution genomewide RNA binding of Nab3 in three independent biological replicates of these S.cerevisiae cells, normalized by spiked-in S.pombe lysates. These data allow the useful assessment of the extent of co-transcriptional binding decrease of Nab3 by decoy ncRNA titration and will be valuable for further analyses of NNS targeting mechanisms.

DNA molecular combing-based replication fork directionality profiling.

The replication strategy of metazoan genomes is still unclear, mainly because definitive maps of replication origins are missing. High-throughput methods are based on population average and thus may exclusively identify efficient initiation sites, whereas inefficient origins go undetected. Single-molecule analyses of specific loci can detect both common and rare initiation events along the targeted regions. However, these usually concentrate on positioning individual events, which only gives an overview of the replication dynamics. Here, we computed the replication fork directionality (RFD) profiles of two large genes in different transcriptional states in chicken DT40 cells, namely untranscribed and transcribed DMD and CCSER1 expressed at WT levels or overexpressed, by aggregating hundreds of oriented replication tracks detected on individual DNA fibres stretched by molecular combing. These profiles reconstituted RFD domains composed of zones of initiation flanking a zone of termination originally observed in mammalian genomes and were highly consistent with independent population-averaging profiles generated by Okazaki fragment sequencing. Importantly, we demonstrate that inefficient origins do not appear as detectable RFD shifts, explaining why dispersed initiation has remained invisible to population-based assays. Our method can both generate quantitative profiles and identify discrete events, thereby constituting a comprehensive approach to study metazoan genome replication.

The C-terminal Domain of piggyBac Transposase Is Not Required for DNA Transposition.

PiggyBac(PB)-like elements (pble) are members of a eukaryotic DNA transposon family. This family is of interest to evolutionary genomics because pble transposases have been domesticated at least 9 times in vertebrates. The amino acid sequence of pble transposases can be split into three regions: an acidic N-terminal domain (~100 aa), a central domain (~400 aa) containing a DD[D/E] catalytic triad, and a cysteine-rich domain (CRD; ~90 aa). Two recent reports suggested that a functional CRD is required for pble transposase activity. Here we found that two CRD-deficient pble transposases, a PB variant and an isoform encoded by the domesticated PB-derived vertebrate transposase gene 5 (pgbd5) trigger transposition of the Ifp2 pble. When overexpressed in HeLa cells, these CRD-deficient transposases can insert Ifp2 elements with proper and improper transposon ends, associated with deleterious effects on cells. Finally, we found that mouse CRD-deficient transposase Pgbd5, as well as PB, do not insert pbles at random into chromosomes. Transposition events occurred more often in genic regions, in the neighbourhood of the transcription start sites and were often found in genes predominantly expressed in the human central nervous system.

The Rhinobiome of Exacerbated Wheezers and Asthmatics: Insights From a German Pediatric Exacerbation Network.

Although the nose, as a gateway for organism-environment interactions, may have a key role in asthmatic exacerbation, the rhinobiome of exacerbated children with asthma was widely neglected to date. The aim of this study is to understand the microbiome, the microbial immunology, and the proteome of exacerbated children and adolescents with wheeze and asthma. Considering that a certain proportion of wheezers may show a progression to asthma, the comparison of both groups provides important information regarding clinical and phenotype stratification. Thus, deep nasopharyngeal swab specimens, nasal epithelial spheroid (NAEsp) cultures, and blood samples of acute exacerbated wheezers (WH), asthmatics (AB), and healthy controls (HC) were used for culture (n = 146), 16 S-rRNA gene amplicon sequencing (n = 64), and proteomic and cytokine analyses. Interestingly, Proteobacteria were over-represented in WH, whereas Firmicutes and Bacteroidetes were associated with AB. In contrast, Actinobacteria commonly colonized HCs. Moreover, Staphylococcaceae, Enterobacteriaceae, Burkholderiaceae, Xanthobacteraceae, and Sphingomonadaceae were significantly more abundant in AB compared to WH and HC. The α-diversity analyses demonstrated an increase of bacterial abundance levels in atopic AB and a decrease in WH samples. Microbiome profiles of atopic WH differed significantly from atopic AB, whereby atopic samples of WH were more homogeneous than those of non-atopic subjects. The NAEsp bacterial exposure experiments provided a disrupted epithelial cell integrity, a cytokine release, and cohort-specific proteomic differences especially for Moraxella catarrhalis cultures. This comprehensive dataset contributes to a deeper insight into the poorly understood plasticity of the nasal microbiota, and, in particular, may enforce our understanding in the pathogenesis of asthma exacerbation in childhood.

Pattern and causes of the establishment of the invasive bacterial potato pathogen Dickeya solani and of the maintenance of the resident pathogen D. dianthicola.

Invasive pathogens can be a threat when they affect human health, food production or ecosystem services, by displacing resident species, and we need to understand the cause of their establishment. We studied the patterns and causes of the establishment of the pathogen Dickeya solani that recently invaded potato agrosystems in Europe by assessing its invasion dynamics and its competitive ability against the closely related resident D. dianthicola species. Epidemiological records over one decade in France revealed the establishment of D. solani and the maintenance of the resident D. dianthicola in potato fields exhibiting blackleg symptoms. Using experimentations, we showed that D. dianthicola caused a higher symptom incidence on aerial parts of potato plants than D. solani, while D. solani was more aggressive on tubers (i.e. with more severe symptoms). In co-infection assays, D. dianthicola outcompeted D. solani in aerial parts, while the two species co-existed in tubers. A comparison of 76 D. solani genomes (56 of which have been sequenced here) revealed balanced frequencies of two previously uncharacterized alleles, VfmBPro and VfmBSer , at the vfmB virulence gene. Experimental inoculations showed that the VfmBSer population was more aggressive on tubers, while the VfmBPro population outcompeted the VfmBSer population in stem lesions, suggesting an important role of the vfmB virulence gene in the ecology of the pathogens. This study thus brings novel insights allowing a better understanding of the pattern and causes of the D.solani invasion into potato production agrosystems, and the reasons why the endemic D. dianthicola nevertheless persisted.