Dual Roles of Poly(dA:dT) Tracts in Replication Initiation and Fork Collapse.

Replication origins, fragile sites, and rDNA have been implicated as sources of chromosomal instability. However, the defining genomic features of replication origins and fragile sites are among the least understood elements of eukaryote genomes. Here, we map sites of replication initiation and breakage in primary cells at high resolution. We find that replication initiates between transcribed genes within nucleosome-depleted structures established by long asymmetrical poly(dA:dT) tracts flanking the initiation site. Paradoxically, long (>20 bp) (dA:dT) tracts are also preferential sites of polar replication fork stalling and collapse within early-replicating fragile sites (ERFSs) and late-replicating common fragile sites (CFSs) and at the rDNA replication fork barrier. Poly(dA:dT) sequences are fragile because long single-strand poly(dA) stretches at the replication fork are unprotected by the replication protein A (RPA). We propose that the evolutionary expansion of poly(dA:dT) tracts in eukaryotic genomes promotes replication initiation, but at the cost of chromosome fragility.

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.

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.

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 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.

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.

The Third Revolution in Sequencing Technology.

Forty years ago the advent of Sanger sequencing was revolutionary as it allowed complete genome sequences to be deciphered for the first time. A second revolution came when next-generation sequencing (NGS) technologies appeared, which made genome sequencing much cheaper and faster. However, NGS methods have several drawbacks and pitfalls, most notably their short reads. Recently, third-generation/long-read methods appeared, which can produce genome assemblies of unprecedented quality. Moreover, these technologies can directly detect epigenetic modifications on native DNA and allow whole-transcript sequencing without the need for assembly. This marks the third revolution in sequencing technology. Here we review and compare the various long-read methods. We discuss their applications and their respective strengths and weaknesses and provide future perspectives.

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.

Developmental and cancer-associated plasticity of DNA replication preferentially targets GC-poor, lowly expressed and late-replicating regions.

The spatiotemporal program of metazoan DNA replication is regulated during development and altered in cancers. We have generated novel OK-seq, Repli-seq and RNA-seq data to compare the DNA replication and gene expression programs of twelve cancer and non-cancer human cell types. Changes in replication fork directionality (RFD) determined by OK-seq are widespread but more frequent within GC-poor isochores and largely disconnected from transcription changes. Cancer cell RFD profiles cluster with non-cancer cells of similar developmental origin but not with different cancer types. Importantly, recurrent RFD changes are detected in specific tumour progression pathways. Using a model for establishment and early progression of chronic myeloid leukemia (CML), we identify 1027 replication initiation zones (IZs) that progressively change efficiency during long-term expression of the BCR-ABL1 oncogene, being twice more often downregulated than upregulated. Prolonged expression of BCR-ABL1 results in targeting of new IZs and accentuation of previous efficiency changes. Targeted IZs are predominantly located in GC-poor, late replicating gene deserts and frequently silenced in late CML. Prolonged expression of BCR-ABL1 results in massive deletion of GC-poor, late replicating DNA sequences enriched in origin silencing events. We conclude that BCR-ABL1 expression progressively affects replication and stability of GC-poor, late-replicating regions during CML progression.

Postembryonic Fish Brain Proliferation Zones Exhibit Neuroepithelial-Type Gene Expression Profile.

In mammals, neuroepithelial cells play an essential role in embryonic neurogenesis, whereas glial stem cells are the principal source of neurons at postembryonic stages. By contrast, neuroepithelial-like stem/progenitor (NE) cells have been shown to be present throughout life in teleosts. We used three-dimensional (3D) reconstructions of cleared transgenic wdr12:GFP medaka brains to demonstrate that this cell type is widespread in juvenile and to identify new regions containing NE cells. We established the gene expression profile of optic tectum (OT) NE cells by cell sorting followed by RNA-seq. Our results demonstrate that most OT NE cells are indeed active stem cells and that some of them exhibit long G2 phases. We identified several novel pathways (e.g., DNA repair pathways) potentially involved in NE cell homeostasis. In situ hybridization studies showed that all NE populations in the postembryonic medaka brain have a similar molecular signature. Our findings highlight the importance of NE progenitors in medaka and improve our understanding of NE-cell biology. These cells are potentially useful not only for neural stem cell studies but also for improving the characterization of neurodevelopmental diseases, such as microcephaly. Stem Cells 2017;35:1505-1518.

Finding sRNA-associated phenotypes by competition assays: An example with Staphylococcus aureus.

Bacteria optimize their fitness in response to a changing environment by tight regulation of gene expression. Regulation can be controlled at both transcriptional and post-transcriptional levels via key players such as sigma factors, regulatory proteins and regulatory RNAs. The identification of phenotypes associated with gene deletions is the established method for finding gene functions but may require testing many conditions for each studied mutant. As regulatory RNAs often contribute to fine-tuning gene expression, phenotypes associated with their inactivation are often weak and difficult to detect. Nevertheless, minor phenotypes conferring modest advantages, may allow bacteria to emerge after some generations under selective pressure. A strategy employing DNA barcodes can be used to perform competition experiments between mutants and to monitor fitness associated with mutations in different growth conditions. We combined this strategy with deep sequencing to study regulatory RNAs in Staphylococcus aureus, a major opportunistic pathogen.

Ten years of next-generation sequencing technology.

Ten years ago next-generation sequencing (NGS) technologies appeared on the market. During the past decade, tremendous progress has been made in terms of speed, read length, and throughput, along with a sharp reduction in per-base cost. Together, these advances democratized NGS and paved the way for the development of a large number of novel NGS applications in basic science as well as in translational research areas such as clinical diagnostics, agrigenomics, and forensic science. Here we provide an overview of the evolution of NGS and discuss the most significant improvements in sequencing technologies and library preparation protocols. We also explore the current landscape of NGS applications and provide a perspective for future developments.

Library preparation methods for next-generation sequencing: tone down the bias.

Next-generation sequencing (NGS) has caused a revolution in biology. NGS requires the preparation of libraries in which (fragments of) DNA or RNA molecules are fused with adapters followed by PCR amplification and sequencing. It is evident that robust library preparation methods that produce a representative, non-biased source of nucleic acid material from the genome under investigation are of crucial importance. Nevertheless, it has become clear that NGS libraries for all types of applications contain biases that compromise the quality of NGS datasets and can lead to their erroneous interpretation. A detailed knowledge of the nature of these biases will be essential for a careful interpretation of NGS data on the one hand and will help to find ways to improve library quality or to develop bioinformatics tools to compensate for the bias on the other hand. In this review we discuss the literature on bias in the most common NGS library preparation protocols, both for DNA sequencing (DNA-seq) as well as for RNA sequencing (RNA-seq). Strikingly, almost all steps of the various protocols have been reported to introduce bias, especially in the case of RNA-seq, which is technically more challenging than DNA-seq. For each type of bias we discuss methods for improvement with a view to providing some useful advice to the researcher who wishes to convert any kind of raw nucleic acid into an NGS library.

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.

A cis-regulatory signature for chordate anterior neuroectodermal genes.

One of the striking findings of comparative developmental genetics was that expression patterns of core transcription factors are extraordinarily conserved in bilaterians. However, it remains unclear whether cis-regulatory elements of their target genes also exhibit common signatures associated with conserved embryonic fields. To address this question, we focused on genes that are active in the anterior neuroectoderm and non-neural ectoderm of the ascidian Ciona intestinalis. Following the dissection of a prototypic anterior placodal enhancer, we searched all genomic conserved non-coding elements for duplicated motifs around genes showing anterior neuroectodermal expression. Strikingly, we identified an over-represented pentamer motif corresponding to the binding site of the homeodomain protein OTX, which plays a pivotal role in the anterior development of all bilaterian species. Using an in vivo reporter gene assay, we observed that 10 of 23 candidate cis-regulatory elements containing duplicated OTX motifs are active in the anterior neuroectoderm, thus showing that this cis-regulatory signature is predictive of neuroectodermal enhancers. These results show that a common cis-regulatory signature corresponding to K50-Paired homeodomain transcription factors is found in non-coding sequences flanking anterior neuroectodermal genes in chordate embryos. Thus, field-specific selector genes impose architectural constraints in the form of combinations of short tags on their target enhancers. This could account for the strong evolutionary conservation of the regulatory elements controlling field-specific selector genes responsible for body plan formation.

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.

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.