Next generation sequencing has a potential to revolutionize the management of cancer patients within the framework of precision oncology. Nevertheless, lack of standardization decelerated entering of the technology into the clinical testing space. Here we dissected a number of common problems of NGS diagnostics in oncology and introduced ways they can be resolved.

Progression of RNA polymerase II (RNAPII) transcription relies on the appropriately positioned activities of elongation factors. The resulting profile of factors and chromatin signatures along transcription units provides a "positional information system" for transcribing RNAPII. Here, we investigate a chromatin-based mechanism that suppresses intragenic initiation of RNAPII transcription. We demonstrate that RNAPII transcription across gene promoters represses their function in plants. This repression is characterized by reduced promoter-specific molecular signatures and increased molecular signatures associated with RNAPII elongation. The conserved FACT histone chaperone complex is required for this repression mechanism. Genome-wide Transcription Start Site (TSS) mapping reveals thousands of discrete intragenic TSS positions in fact mutants, including downstream promoters that initiate alternative transcript isoforms. We find that histone H3 lysine 4 mono-methylation (H3K4me1), an Arabidopsis RNAPII elongation signature, is enriched at FACT-repressed intragenic TSSs. Our analyses suggest that FACT is required to repress intragenic TSSs at positions that are in part characterized by elevated H3K4me1 levels. In sum, conserved and plant-specific chromatin features correlate with the co-transcriptional repression of intragenic TSSs. Our insights into TSS repression by RNAPII transcription promise to inform the regulation of alternative transcript isoforms and the characterization of gene regulation through the act of pervasive transcription across eukaryotic genomes.

Most DNA in the genomes of higher organisms does not code for proteins. RNA Polymerase II (Pol II) transcribes non-coding DNA into long non-coding RNAs (lncRNAs), but biological roles of lncRNA are unclear. We find that mutations in the yeast lncRNA CUT60 result in poor growth. Defective termination of CUT60 transcription causes read-through transcription across the ATP16 gene promoter. Read-through transcription localizes chromatin signatures associated with Pol II elongation to the ATP16 promoter. The act of Pol II elongation across this promoter represses functional ATP16 expression by a Transcriptional Interference (TI) mechanism. Atp16p function in the mitochondrial ATP-synthase complex promotes mitochondrial DNA stability. ATP16 repression by TI through inefficient termination of CUT60 therefore triggers mitochondrial genome loss. Our results expand the functional and mechanistic implications of non-coding DNA in eukaryotes by highlighting termination of nuclear lncRNA transcription as mechanism to stabilize an organellar genome.

As the use of next-generation sequencing (NGS) for the Mendelian diseases diagnosis is expanding, the performance of this method has to be improved in order to achieve higher quality. Typically, performance measures are considered to be designed in the context of each application and, therefore, account for a spectrum of clinically relevant variants. We present EphaGen, a new computational methodology for bioinformatics quality control (QC). Given a single NGS dataset in BAM format and a pre-compiled VCF-file of targeted clinically relevant variants it associates this dataset with a single arbiter parameter. Intrinsically, EphaGen estimates the probability to miss any variant from the defined spectrum within a particular NGS dataset. Such performance measure virtually resembles the diagnostic sensitivity of given NGS dataset. Here we present case studies of the use of EphaGen in context of BRCA1/2 and CFTR sequencing in a series of 14 runs across 43 blood samples and 504 publically available NGS datasets. EphaGen is superior to conventional bioinformatics metrics such as coverage depth and coverage uniformity. We recommend using this software as a QC step in NGS studies in the clinical context. Availability: https://github.com/m4merg/EphaGen or https://hub.docker.com/r/m4merg/ephagen.

Bone injuries represent a major social and financial impairment, commonly requiring surgical intervention due to a limited healing capacity of the tissue, particularly regarding critical-sized defects and non-union fractures. Regenerative medicine with the application of bone implants has been developing in the past decades towards the manufacturing of appropriate devices. This work intended to evaluate medical 316L stainless steel (SS)-based devices covered by a polymer poly (L-lactic acid) (PLLA) coating for bone lesion mechanical and functional support. SS316L devices were subjected to a previously described silanization process, following a three-layer PLLA film coating. Devices were further characterized and evaluated towards their cytocompatibility and osteogenic potential using human dental pulp stem cells, and biocompatibility via subcutaneous implantation in a rat animal model. Results demonstrated PLLA-SS316L devices to present superior in vitro and in vivo outcomes and suggested the PLLA coating to provide osteo-inductive properties to the device. Overall, this work represents a preliminary study on PLLA-SS316L devices' potential towards bone tissue regenerative techniques, showing promising outcomes for bone lesion support.

The quality of gene annotation determines the interpretation of results obtained in transcriptomic studies. The growing number of genome sequence information calls for experimental and computational pipelines for de novo transcriptome annotation. Ideally, gene and transcript models should be called from a limited set of key experimental data.

The liver plays a central role in the maintenance of homeostasis and health in general. However, there is substantial inter-individual variation in hepatic gene expression, and although numerous genetic factors have been identified, less is known about the epigenetic factors.

Genotyping of cell-free DNA (cfDNA) in plasma samples has the potential to allow for a noninvasive assessment of tumor biology, avoiding the inherent shortcomings of tissue biopsy. Next generation sequencing (NGS), a leading technology for liquid biopsy analysis, continues to be hurdled with several major issues with cfDNA samples, including low cfDNA concentration and high fragmentation. In this study, by employing Ion Torrent PGM semiconductor technology, we performed a comparison between two multi-biomarker amplicon-based NGS panels characterized by a substantial difference in average amplicon length. In course of the analysis of the peripheral blood from 13 diagnostic non-small cell lung cancer patients, equivalence of two panels, in terms of overall diagnostic sensitivity and specificity was shown. A pairwise comparison of the allele frequencies for the same somatic variants obtained from the pairs of panel-specific amplicons, demonstrated an identical analytical sensitivity in range of 140 to 170 bp amplicons in size. Further regression analysis between amplicon length and its coverage, illustrated that NGS sequencing of plasma cfDNA equally tolerates amplicons with lengths in the range of 120 to 170 bp. To increase the sensitivity of mutation detection in cfDNA, we performed a computational analysis of the features associated with genome-wide nucleosome maps, evident from the data on the prevalence of cfDNA fragments of certain sizes and their fragmentation patterns. By leveraging the support vector machine-based machine learning approach, we showed that a combination of nucleosome map associated features with GC content, results in the increased accuracy of prediction of high inter-sample sequencing coverage variation (areas under the receiver operating curve: 0.75, 95% CI: 0.750-0.752 vs. 0.65, 95% CI: 0.63-0.67). Thus, nucleosome-guided fragmentation should be utilized as a guide to design amplicon-based NGS panels for the genotyping of cfDNA samples.

RNA polymerase II (RNAPII) transcription is crucial for gene expression. RNAPII density peaks at gene boundaries, associating these key regions for gene expression control with limited RNAPII movement. The connections between RNAPII transcription speed and gene regulation in multicellular organisms are poorly understood. Here, we directly modulate RNAPII transcription speed by point mutations in the second largest subunit of RNAPII in Arabidopsis thaliana. A RNAPII mutation predicted to decelerate transcription is inviable, while accelerating RNAPII transcription confers phenotypes resembling auto-immunity. Nascent transcription profiling revealed that RNAPII complexes with accelerated transcription clear stalling sites at both gene ends, resulting in read-through transcription. The accelerated transcription mutant NRPB2-Y732F exhibits increased association with 5' splice site (5'SS) intermediates and enhanced splicing efficiency. Our findings highlight potential advantages of RNAPII stalling through local reduction in transcription speed to optimize gene expression for the development of multicellular organisms.

Cystic fibrosis (CF) is one of the most common life-threatening genetic disorders. Around 2000 variants in the CFTR gene have been identified, with some proportion known to be pathogenic and 300 disease-causing mutations have been characterized in detail by CFTR2 database, which complicates its analysis with conventional methods.

Small interfering RNAs (siRNAs) with N-acetylgalactosamine (GalNAc) conjugation for improved liver uptake represent an emerging class of drugs to treat liver diseases. Understanding how pharmacokinetics and pharmacodynamics translate is pivotal for in vivo study design and human dose prediction. However, the literature is sparse on translational data for this modality, and pharmacokinetics in the liver is seldom measured. To overcome these difficulties, we collected time-course biomarker data for 11 GalNAc-siRNAs in various species and applied the kinetic-pharmacodynamic modeling approach to estimate the biophase (liver) half-life and the potency. Our analysis indicates that the biophase half-life is 0.6-3 weeks in mouse, 1-8 weeks in monkey, and 1.5-14 weeks in human. For individual siRNAs, the biophase half-life is 1-8 times longer in human than in mouse, and generally 1-3 times longer in human than in monkey. The analysis indicates that the siRNAs are more potent in human than in mouse and monkey.

To improve the epigenomic analysis of tissues rich in 5-hydroxymethylcytosine (hmC), we developed a novel protocol called TAB-Methyl-SEQ, which allows for single base resolution profiling of both hmC and 5-methylcytosine by targeted next-generation sequencing. TAB-Methyl-SEQ data were extensively validated by a set of five methodologically different protocols. Importantly, these extensive cross-comparisons revealed that protocols based on Tet1-assisted bisulfite conversion provided more precise hmC values than TrueMethyl-based methods. A total of 109 454 CpG sites were analyzed by TAB-Methyl-SEQ for mC and hmC in 188 genes from 20 different adult human livers. We describe three types of variability of hepatic hmC profiles: (i) sample-specific variability at 40.8% of CpG sites analyzed, where the local hmC values correlate to the global hmC content of livers (measured by LC-MS), (ii) gene-specific variability, where hmC levels in the coding regions positively correlate to expression of the respective gene and (iii) site-specific variability, where prominent hmC peaks span only 1 to 3 neighboring CpG sites. Our data suggest that both the gene- and site-specific components of hmC variability might contribute to the epigenetic control of hepatic genes. The protocol described here should be useful for targeted DNA analysis in a variety of applications.

Nucleosome-depleted regions (NDRs) at gene promoters support initiation of RNA polymerase II transcription. Interestingly, transcription often initiates in both directions, resulting in an mRNA and a divergent non-coding (DNC) transcript of unclear purpose. Here, we characterized the genetic architecture and molecular mechanism of DNC transcription in budding yeast. Using high-throughput reverse genetic screens based on quantitative single-cell fluorescence measurements, we identified the Hda1 histone deacetylase complex (Hda1C) as a repressor of DNC transcription. Nascent transcription profiling showed a genome-wide role of Hda1C in repression of DNC transcription. Live-cell imaging of transcription revealed that mutations in the Hda3 subunit increased the frequency of DNC transcription. Hda1C contributed to decreased acetylation of histone H3 in DNC transcription regions, supporting DNC transcription repression by histone deacetylation. Our data support the interpretation that DNC transcription results as a consequence of the NDR-based architecture of eukaryotic promoters, but that it is governed by locus-specific repression to maintain genome fidelity.

Most DNA in the genomes of higher organisms does not encode proteins, yet much is transcribed by RNA polymerase II (RNAPII) into long non-coding RNAs (lncRNAs). The biological significance of most lncRNAs is largely unclear. Here, we identify a lncRNA (SVALKA) in a cold-sensitive region of the Arabidopsis genome. Mutations in SVALKA affect CBF1 expression and freezing tolerance. RNAPII read-through transcription of SVALKA results in a cryptic lncRNA overlapping CBF1 on the antisense strand, termed asCBF1. Our molecular dissection reveals that CBF1 is suppressed by RNAPII collision stemming from the SVALKA-asCBF1 lncRNA cascade. The SVALKA-asCBF1 cascade provides a mechanism to tightly control CBF1 expression and timing that could be exploited to maximize freezing tolerance with mitigated fitness costs. Our results provide a compelling example of local gene regulation by lncRNA transcription having a profound impact on the ability of plants to appropriately acclimate to challenging environmental conditions.

DNA methylation is one of the most important epigenetic alterations involved in the control of gene expression. Bisulfite sequencing of genomic DNA is currently the only method to study DNA methylation patterns at single-nucleotide resolution. Hence, next-generation sequencing of bisulfite-converted DNA is the method of choice to investigate DNA methylation profiles at the genome-wide scale. Nevertheless, whole genome sequencing for analysis of human methylomes is expensive, and a method for targeted gene analysis would provide a good alternative in many cases where the primary interest is restricted to a set of genes. Here, we report the successful use of a custom Agilent SureSelect Target Enrichment system for the hybrid capture of bisulfite-converted DNA. We prepared bisulfite-converted next-generation sequencing libraries, which are enriched for the coding and regulatory regions of 174 ADME genes (i.e. genes involved in the metabolism and distribution of drugs). Sequencing of these libraries on Illumina's HiSeq2000 revealed that the method allows a reliable quantification of methylation levels of CpG sites in the selected genes, and validation of the method using pyrosequencing and the Illumina 450K methylation BeadChips revealed good concordance.

Temperature profoundly affects the kinetics of biochemical reactions, yet how large molecular complexes such as the transcription machinery accommodate changing temperatures to maintain cellular function is poorly understood. Here, we developed plant native elongating transcripts sequencing (plaNET-seq) to profile genome-wide nascent RNA polymerase II (RNAPII) transcription during the cold-response of Arabidopsis thaliana with single-nucleotide resolution. Combined with temporal resolution, these data revealed transient genome-wide reprogramming of nascent RNAPII transcription during cold, including characteristics of RNAPII elongation and thousands of non-coding transcripts connected to gene expression. Our results suggest a role for promoter-proximal RNAPII stalling in predisposing genes for transcriptional activation during plant-environment interactions. At gene 3'-ends, cold initially facilitated transcriptional termination by limiting the distance of read-through transcription. Within gene bodies, cold reduced the kinetics of co-transcriptional splicing leading to increased intragenic stalling. Our data resolved multiple distinct mechanisms by which temperature transiently altered the dynamics of nascent RNAPII transcription and associated RNA processing, illustrating potential biotechnological solutions and future focus areas to promote food security in the context of a changing climate.

How the noncoding genome affects cellular functions is a key biological question. A particular challenge is to distinguish the effects of noncoding DNA elements from long noncoding RNAs (lncRNAs) that coincide at the same loci. Here, we identified the flowering-associated intergenic lncRNA (FLAIL) in Arabidopsis through early flowering flail mutants. Expression of FLAIL RNA from a different chromosomal location in combination with strand-specific RNA knockdown characterized FLAIL as a trans-acting RNA molecule. FLAIL directly binds to differentially expressed target genes that control flowering via RNA-DNA interactions through conserved sequence motifs. FLAIL interacts with protein and RNA components of the spliceosome to affect target mRNA expression through co-transcriptional alternative splicing (AS) and linked chromatin regulation. In the absence of FLAIL, splicing defects at the direct FLAIL target flowering gene LACCASE 8 (LAC8) correlated with reduced mRNA expression. Double mutant analyses support a model where FLAIL-mediated splicing of LAC8 promotes its mRNA expression and represses flowering. Our study suggests lncRNAs as accessory components of the spliceosome that regulate AS and gene expression to impact organismal development.

A fraction of patients referred for complex molecular profiling of biopsied tumors may harbor germline variants in genes associated with the development of hereditary cancer syndromes (HCS). Neither the bioinformatic analysis nor the reporting of such incidental germline findings are standardized.