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.

GC content, but not nucleosome positioning, directly contributes to intron splicing efficiency in Paramecium.

Eukaryotic genes are interrupted by introns that must be accurately spliced from mRNA precursors. With an average length of 25 nt, the more than 90,000 introns of Paramecium tetraurelia stand among the shortest introns reported in eukaryotes. The mechanisms specifying the correct recognition of these tiny introns remain poorly understood. Splicing can occur cotranscriptionally, and it has been proposed that chromatin structure might influence splice site recognition. To investigate the roles of nucleosome positioning in intron recognition, we determined the nucleosome occupancy along the P. tetraurelia genome. We show that P. tetraurelia displays a regular nucleosome array with a nucleosome repeat length of ∼151 bp, among the smallest periodicities reported. Our analysis has revealed that introns are frequently associated with inter-nucleosomal DNA, pointing to an evolutionary constraint favoring introns at the AT-rich nucleosome edge sequences. Using accurate splicing efficiency data from cells depleted for nonsense-mediated decay effectors, we show that introns located at the edge of nucleosomes display higher splicing efficiency than those at the center. However, multiple regression analysis indicates that the low GC content of introns, rather than nucleosome positioning, is associated with high splicing efficiency. Our data reveal a complex link between GC content, nucleosome positioning, and intron evolution in Paramecium.

OKseqHMM: a genome-wide replication fork directionality analysis toolkit.

During each cell division, tens of thousands of DNA replication origins are co-ordinately activated to ensure the complete duplication of the human genome. However, replication fork progression can be challenged by many factors, including co-directional and head-on transcription-replication conflicts (TRC). Head-on TRCs are more dangerous for genome integrity. To study the direction of replication fork movement and TRCs, we developed a bioinformatics toolkit called OKseqHMM (https://github.com/CL-CHEN-Lab/OK-Seq, https://doi.org/10.5281/zenodo.7428883). Then, we used OKseqHMM to analyse a large number of datasets obtained by Okazaki fragment sequencing to directly measure the genome-wide replication fork directionality (RFD) and to accurately predict replication initiation and termination at a fine resolution in organisms including yeast, mouse and human. We also successfully applied our analysis to other genome-wide sequencing techniques that also contain RFD information (e.g. eSPAN, TrAEL-seq). Our toolkit can be used to predict replication initiation and fork progression direction genome-wide in a wide range of cell models and growth conditions. Comparing the replication and transcription directions allows identifying loci at risk of TRCs, particularly head-on TRCs, and investigating their role in genome instability by checking DNA damage data, which is of prime importance for human health.

Contrasting Gene Decay in Subterranean Vertebrates: Insights from Cavefishes and Fossorial Mammals.

Evolution sometimes proceeds by loss, especially when structures and genes become dispensable after an environmental shift relaxes functional constraints. Subterranean vertebrates are outstanding models to analyze this process, and gene decay can serve as a readout. We sought to understand some general principles on the extent and tempo of the decay of genes involved in vision, circadian clock, and pigmentation in cavefishes. The analysis of the genomes of two Cuban species belonging to the genus Lucifuga provided evidence for the largest loss of eye-specific genes and nonvisual opsin genes reported so far in cavefishes. Comparisons with a recently evolved cave population of Astyanax mexicanus and three species belonging to the Chinese tetraploid genus Sinocyclocheilus revealed the combined effects of the level of eye regression, time, and genome ploidy on eye-specific gene pseudogenization. The limited extent of gene decay in all these cavefishes and the very small number of loss-of-function mutations per pseudogene suggest that their eye degeneration may not be very ancient, ranging from early to late Pleistocene. This is in sharp contrast with the identification of several vision genes carrying many loss-of-function mutations in ancient fossorial mammals, further suggesting that blind fishes cannot thrive more than a few million years in cave ecosystems.

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.

Tight protein-DNA interactions favor gene silencing.

The heterochromatin-like structure formed by the yeast silent information regulator complex (SIR) represses transcription at the silent mating type loci and telomeres. Here, we report that tight protein-DNA complexes induce ectopic recruitment of the SIR complex, promoting gene silencing and changes in subnuclear localization when cis-acting elements are nearby. Importantly, lack of the replication fork-associated helicase Rrm3 enhances this induced gene repression. Additionally, Sir3 and Sir4 are enriched genome-wide at natural replication pause sites, including tRNA genes. Consistently, inserting a tRNA gene promotes SIR-mediated silencing of a nearby gene. These results reveal that replication stress arising from tight DNA-protein interactions favors heterochromatin formation.

Evidence for late Pleistocene origin of Astyanax mexicanus cavefish.

BACKGROUND: Cavefish populations belonging to the Mexican tetra species Astyanax mexicanus are outstanding models to study the tempo and mode of adaptation to a radical environmental change. They are currently assigned to two main groups, the so-called "old" and "new" lineages, which would have populated several caves independently and at different times. However, we do not have yet accurate estimations of the time frames of evolution of these populations. RESULTS: We reanalyzed the geographic distribution of mitochondrial and nuclear DNA polymorphisms and we found that these data do not support the existence of two cavefish lineages. Using IMa2, a program that allows dating population divergence in addition to demographic parameters, we found that microsatellite polymorphism strongly supports a very recent origin of cave populations (< 20,000 years). We identified a large number of single-nucleotide polymorphisms (SNPs) in transcript sequences of pools of embryos (Pool-seq) belonging to Pachón cave population and a surface population from Texas. Based on summary statistics that can be computed with this SNP data set together with simulations of evolution of SNP polymorphisms in two recently isolated populations, we looked for sets of demographic parameters that allow the computation of summary statistics with simulated populations that are similar to the ones with the sampled populations. In most simulations for which we could find a good fit between the summary statistics of observed and simulated data, the best fit occurred when the divergence between simulated populations was less than 30,000 years. CONCLUSIONS: Although it is often assumed that some cave populations have a very ancient origin, a recent origin of these populations is strongly supported by our analyses of independent sets of nuclear DNA polymorphism. Moreover, the observation of two divergent haplogroups of mitochondrial and nuclear genes with different geographic distributions support a recent admixture of two divergent surface populations, before the isolation of cave populations. If cave populations are indeed only several thousand years old, many phenotypic changes observed in cavefish would thus have mainly involved the fixation of genetic variants present in surface fish populations and within a very short period of time.

Systematic comparison of small RNA library preparation protocols for next-generation sequencing.

BACKGROUND: Next-generation sequencing technologies have revolutionized the study of small RNAs (sRNAs) on a genome-wide scale. However, classical sRNA library preparation methods introduce serious bias, mainly during adapter ligation steps. Several types of sRNA including plant microRNAs (miRNA), piwi-interacting RNAs (piRNA) in insects, nematodes and mammals, and small interfering RNAs (siRNA) in insects and plants contain a 2'-O-methyl (2'-OMe) modification at their 3' terminal nucleotide. This inhibits 3' adapter ligation and makes library preparation particularly challenging. To reduce bias, the NEBNext kit (New England Biolabs) uses polyethylene glycol (PEG), the NEXTflex V2 kit (BIOO Scientific) uses both randomised adapters and PEG, and the novel SMARTer (Clontech) and CATS (Diagenode) kits avoid ligation altogether. Here we compared these methods with Illumina's classical TruSeq protocol regarding the detection of normal and 2' OMe RNAs. In addition, we modified the TruSeq and NEXTflex protocols to identify conditions that improve performance. RESULTS: Among the five kits tested with their respective standard protocols, the SMARTer and CATS kits had the lowest levels of bias but also had a strong formation of side products, and as a result performed relatively poorly with biological samples; NEXTflex detected the largest numbers of different miRNAs. The use of a novel type of randomised adapters called MidRand-Like (MRL) adapters and PEG improved the detection of 2' OMe RNAs both in the TruSeq as well as in the NEXTflex protocol. CONCLUSIONS: While it is commonly accepted that biases in sRNA library preparation protocols are mainly due to adapter ligation steps, the ligation-free protocols were not the best performing methods. Our modified versions of the TruSeq and NEXTflex protocols provide an improved tool for the study of 2' OMe RNAs.

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.

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.

Nucleoid organization in the radioresistant bacterium Deinococcus radiodurans.

Processes favoring the exceptional resistance to genotoxic stress of Deinococcus radiodurans are not yet completely characterized. It was postulated that its nucleoid and chromosome(s) organization could participate in the DNA double strand break repair process. Here, we investigated the organization of chromosome 1 by localization of three chromosomal loci including oriC, Ter and a locus located in its left arm. For this purpose, we used a ParB-parS system to visualize the position of the loci before and after exposure to γ-rays. By comparing the number of fluorescent foci with the number of copies of the studied loci present in the cells measured by quantitative polymerase chain reaction (qPCR), we demonstrated that the 4-10 copies of chromosome 1 per cell are dispersed within the nucleoid before irradiation, indicating that the chromosome copies are not prealigned. Chromosome segregation is progressive but not co-ordinated, allowing each locus to be paired with its sister during part of the cell cycle. After irradiation, the nucleoid organization is modified, involving a transient alignment of the loci in the late stage of DNA repair and a delay of segregation of the Ter locus. We discuss how these events can influence DNA double strand break repair.

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.

CIRCUS: a package for Circos display of structural genome variations from paired-end and mate-pair sequencing data.

BACKGROUND: Detection of large genomic rearrangements, such as large indels, duplications or translocations is now commonly achieved by next generation sequencing (NGS) approaches. Recently, several tools have been developed to analyze NGS data but the resulting files are difficult to interpret without an additional visualization step. Circos (Genome Res, 19:1639-1645, 2009), a Perl script, is a powerful visualization software that requires setting up numerous configuration files with a large number of parameters to handle. R packages like RCircos (BMC Bioinformatics, 14:244, 2013) or ggbio (Genome Biol, 13:R77, 2012) provide functions to display genomic data as circular Circos-like plots. However, these tools are very general and lack the functions needed to filter, format and adjust specific input genomic data. RESULTS: We implemented an R package called CIRCUS to analyze genomic structural variations. It generates both data and configuration files necessary for Circos, to produce graphs. Only few R pre-requisites are necessary. Options are available to deal with heterogeneous data, various chromosome numbers and multi-scale analysis. CONCLUSION: CIRCUS allows fast and versatile analysis of genomic structural variants with Circos plots for users with limited coding skills.

Zinc-mediated RNA fragmentation allows robust transcript reassembly upon whole transcriptome RNA-Seq.

Whole transcriptome RNA-Seq has emerged as a powerful tool in transcriptomics, enabling genome-wide quantitative analysis of gene expression and qualitative identification of novel coding or non-coding RNA species through transcriptome reassembly. Common protocols for preparation of RNA-Seq libraries include an RNA fragmentation step for which several RNA sizing techniques are commercially available. To date, there is no global information about their putative bias on transcriptome analysis. Here we compared the effects of RNase III- and zinc-mediated RNA fragmentation on transcript expression measurement and transcriptome reassembly in the budding yeast Saccharomyces cerevisiae. We observed that RNA cleavage by RNase III is heterogeneous along transcripts with a striking decrease of autocorrelation between adjacent nucleotides along the transcriptome. This had little impact on mRNA expression measurement, but specific classes of transcripts such as abundant non-coding RNAs were underrepresented in the libraries constructed using RNase III. Furthermore, zinc-mediated fragmentation allows proper reassembly of more transcripts, with more precise 5' and 3' ends. Together, our results show that transcriptome reassembly from RNA-Seq data is very sensitive to the RNA fragmentation technique, and that zinc-mediated fragmentation provides more robust and accurate transcript identification than cleavage by RNase III.

3D chromatin conformation correlates with replication timing and is conserved in resting cells.

Although chromatin folding is known to be of functional importance to control the gene expression program, less is known regarding its interplay with DNA replication. Here, using Circular Chromatin Conformation Capture combined with high-throughput sequencing, we identified megabase-sized self-interacting domains in the nucleus of a human lymphoblastoid cell line, as well as in cycling and resting peripheral blood mononuclear cells (PBMC). Strikingly, the boundaries of those domains coincide with early-initiation zones in every cell types. Preferential interactions have been observed between the consecutive early-initiation zones, but also between those separated by several tens of megabases. Thus, the 3D conformation of chromatin is strongly correlated with the replication timing along the whole chromosome. We furthermore provide direct clues that, in addition to the timing value per se, the shape of the timing profile at a given locus defines its set of genomic contacts. As this timing-related scheme of chromatin organization exists in lymphoblastoid cells, resting and cycling PBMC, this indicates that it is maintained several weeks or months after the previous S-phase. Lastly, our work highlights that the major chromatin changes accompanying PBMC entry into cell cycle occur while keeping largely unchanged the long-range chromatin contacts.

Megabase replication domains along the human genome: relation to chromatin structure and genome organisation.

In higher eukaryotes, the absence of specific sequence motifs, marking the origins of replication has been a serious hindrance to the understanding of (i) the mechanisms that regulate the spatio-temporal replication program, and (ii) the links between origins activation, chromatin structure and transcription. In this chapter, we review the partitioning of the human genome into megabased-size replication domains delineated as N-shaped motifs in the strand compositional asymmetry profiles. They collectively span 28.3% of the genome and are bordered by more than 1,000 putative replication origins. We recapitulate the comparison of this partition of the human genome with high-resolution experimental data that confirms that replication domain borders are likely to be preferential replication initiation zones in the germline. In addition, we highlight the specific distribution of experimental and numerical chromatin marks along replication domains. Domain borders correspond to particular open chromatin regions, possibly encoded in the DNA sequence, and around which replication and transcription are highly coordinated. These regions also present a high evolutionary breakpoint density, suggesting that susceptibility to breakage might be linked to local open chromatin fiber state. Altogether, this chapter presents a compartmentalization of the human genome into replication domains that are landmarks of the human genome organization and are likely to play a key role in genome dynamics during evolution and in pathological situations.

Maternal inheritance of functional centrioles in two parthenogenetic nematodes.

Centrioles are the core constituent of centrosomes, microtubule-organizing centers involved in directing mitotic spindle assembly and chromosome segregation in animal cells. In sexually reproducing species, centrioles degenerate during oogenesis and female meiosis is usually acentrosomal. Centrioles are retained during male meiosis and, in most species, are reintroduced with the sperm during fertilization, restoring centriole numbers in embryos. In contrast, the presence, origin, and function of centrioles in parthenogenetic species is unknown. We found that centrioles are maternally inherited in two species of asexual parthenogenetic nematodes and identified two different strategies for maternal inheritance evolved in the two species. In Rhabditophanes diutinus, centrioles organize the poles of the meiotic spindle and are inherited by both the polar body and embryo. In Disploscapter pachys, the two pairs of centrioles remain close together and are inherited by the embryo only. Our results suggest that maternally-inherited centrioles organize the embryonic spindle poles and act as a symmetry-breaking cue to induce embryo polarization. Thus, in these parthenogenetic nematodes, centrioles are maternally-inherited and functionally replace their sperm-inherited counterparts in sexually reproducing species.

Expression of Dystrophin Dp71 Splice Variants Is Temporally Regulated During Rodent Brain Development.

Dystrophin Dp71 is the major product of the Duchenne muscular dystrophy (DMD) gene in the brain, and its loss in DMD patients and mouse models leads to cognitive impairments. Dp71 is expressed as a range of proteins generated by alternative splicing of exons 71 to 74 and 78, classified in the main Dp71d and Dp71f groups that contain specific C-terminal ends. However, it is unknown whether each isoform has a specific role in distinct cell types, brain regions, and/or stages of brain development. In the present study, we characterized the expression of Dp71 isoforms during fetal (E10.5, E15.5) and postnatal (P1, P7, P14, P21 and P60) mouse and rat brain development. We finely quantified the expression of several Dp71 transcripts by RT-PCR and cloning assays in samples from whole-brain and distinct brain structures. The following Dp71 transcripts were detected: Dp71d, Dp71d∆71, Dp71d∆74, Dp71d∆71,74, Dp71d∆71-74, Dp71f, Dp71f∆71, Dp71f∆74, Dp71f∆71,74, and Dp71fΔ71-74. We found that the Dp71f isoform is the main transcript expressed at E10.5 (> 80%), while its expression is then progressively reduced and replaced by the expression of isoforms of the Dp71d group from E15.5 to postnatal and adult ages. This major finding was confirmed by third-generation nanopore sequencing. In addition, we found that the level of expression of specific Dp71 isoforms varies as a function of postnatal stages and brain structure. Our results suggest that Dp71 isoforms have different and complementary roles during embryonic and postnatal brain development, likely taking part in a variety of maturation processes in distinct cell types.

A novel strategy of transcription regulation by intragenic nucleosome ordering.

Numerous studies of chromatin structure showed that nucleosome free regions (NFRs) located at 5' gene ends contribute to transcription initiation regulation. Here, we determine the role of intragenic chromatin structure on gene expression regulation. We show that, along Saccharomyces cerevisiae genes, nucleosomes are highly organized following two types of architecture that depend only on the distance between the NFRs located at the 5' and 3' gene ends. In the first type, this distance constrains in vivo the positioning of n nucleosomes regularly organized in a "crystal-like" array. In the second type, this distance is such that the corresponding genes can accommodate either n or (n + 1) nucleosomes, thereby displaying two possible crystal-like arrays of n weakly compacted or n + 1 highly compacted nucleosomes. This adaptability confers "bi-stable" properties to chromatin and is a key to its dynamics. Compared to crystal-like genes, bi-stable genes present higher transcriptional plasticity, higher sensitivity to chromatin regulators, higher H3 turnover rate, and lower H2A.Z enrichment. The results strongly suggest that transcription elongation is facilitated by higher chromatin compaction. The data allow us to propose a new paradigm of transcriptional control mediated by the stability and the level of compaction of the intragenic chromatin architecture and open new ways for investigating eukaryotic gene expression regulation.

Impact of replication timing on non-CpG and CpG substitution rates in mammalian genomes.

Neutral nucleotide substitutions occur at varying rates along genomes, and it remains a major issue to unravel the mechanisms that cause these variations and to analyze their evolutionary consequences. Here, we study the role of replication in the neutral substitution pattern. We obtained a high-resolution replication timing profile of the whole human genome by massively parallel sequencing of nascent BrdU-labeled replicating DNA. These data were compared to the neutral substitution rates along the human genome, obtained by aligning human and chimpanzee genomes using macaque and orangutan as outgroups. All substitution rates increase monotonously with replication timing even after controlling for local or regional nucleotide composition, crossover rate, distance to telomeres, and chromatin compaction. The increase in non-CpG substitution rates might result from several mechanisms including the increase in mutation-prone activities or the decrease in efficiency of DNA repair during the S phase. In contrast, the rate of C --> T transitions in CpG dinucleotides increases in later-replicating regions due to increasing DNA methylation level that reflects a negative correlation between timing and gene expression. Similar results are observed in the mouse, which indicates that replication timing is a main factor affecting nucleotide substitution dynamics at non-CpG sites and constitutes a major neutral process driving mammalian genome evolution.