Liquid drop of DNA libraries reveals total genome information.

Conventional "bulk" PCR often yields inefficient and nonuniform amplification of complex templates in DNA libraries, introducing unwanted biases. Amplification of single DNA molecules encapsulated in a myriad of emulsion droplets (emulsion PCR, ePCR) allows the mitigation of this problem. Different ePCR regimes were experimentally analyzed to identify the most robust techniques for enhanced amplification of DNA libraries. A phenomenological mathematical model that forms an essential basis for optimal use of ePCR for library amplification was developed. A detailed description by high-throughput sequencing of amplified DNA-encoded libraries highlights the principal advantages of ePCR over bulk PCR. ePCR outperforms PCR, reduces gross DNA errors, and provides a more uniform distribution of the amplified sequences. The quasi single-molecule amplification achieved via ePCR represents the fundamental requirement in case of complex DNA templates being prone to diversity degeneration and provides a way to preserve the quality of DNA libraries.

Mutation at the Site of Hydroxylation in the Ribosomal Protein uL15 (RPL27a) Causes Specific Changes in the Repertoire of mRNAs Translated in Mammalian Cells.

Ribosomal protein uL15 (RPL27a) carries a specific modification, hydroxylation, at the His39 residue, which neighbors the CCA terminus of the E-site-bound tRNA at the mammalian ribosome. Under hypoxia, the level of hydroxylation of this protein decreases. We transiently transfected HEK293T cells with constructs expressing wild-type uL15 or mutated uL15 (His39Ala) incapable of hydroxylation, and demonstrated that ribosomes containing both proteins are competent in translation. By applying RNA-seq to the total cellular and polysome-associated mRNAs, we identified differentially expressed genes (DEGs) in cells containing exogenous uL15 or its mutant form. Analyzing mRNA features of up- and down-regulated DEGs, we found an increase in the level of more abundant mRNAs and shorter CDSs in cells with uL15 mutant for both translated and total cellular mRNAs. The level of longer and rarer mRNAs, on the contrary, decreased. Our data show how ribosome heterogeneity can change the composition of the translatome and transcriptome, depending on the properties of the translated mRNAs.

Multi-Omics Profiling of Human Endothelial Cells from the Coronary Artery and Internal Thoracic Artery Reveals Molecular but Not Functional Heterogeneity.

Major adverse cardiovascular events occurring upon coronary artery bypass graft surgery are typically accompanied by endothelial dysfunction. Total arterial revascularisation, which employs both left and right internal thoracic arteries instead of the saphenous vein to create a bypass, is associated with better mid- and long-term outcomes. We suggested that molecular profiles of human coronary artery endothelial cells (HCAECs) and human internal mammary artery endothelial cells (HITAECs) are coherent in terms of transcriptomic and proteomic signatures, which were then investigated by RNA sequencing and ultra-high performance liquid chromatography-mass spectrometry, respectively. Both HCAECs and HITAECs overexpressed molecules responsible for the synthesis of extracellular matrix (ECM) components, basement membrane assembly, cell-ECM adhesion, organisation of intercellular junctions, and secretion of extracellular vesicles. HCAECs were characterised by higher enrichment with molecular signatures of basement membrane construction, collagen biosynthesis and folding, and formation of intercellular junctions, whilst HITAECs were notable for augmented pro-inflammatory signaling, intensive synthesis of proteins and nitrogen compounds, and enhanced ribosome biogenesis. Despite HCAECs and HITAECs showing a certain degree of molecular heterogeneity, no specific markers at the protein level have been identified. Coherence of differentially expressed molecular categories in HCAECs and HITAECs suggests synergistic interactions between these ECs in a bypass surgery scenario.

Transcriptomic Analysis of CRISPR/Cas9-Mediated PARP1-Knockout Cells under the Influence of Topotecan and TDP1 Inhibitor.

Topoisomerase 1 (TOP1) is an enzyme that regulates DNA topology and is essential for replication, recombination, and other processes. The normal TOP1 catalytic cycle involves the formation of a short-lived covalent complex with the 3' end of DNA (TOP1 cleavage complex, TOP1cc), which can be stabilized, resulting in cell death. This fact substantiates the effectiveness of anticancer drugs-TOP1 poisons, such as topotecan, that block the relegation of DNA and fix TOP1cc. Tyrosyl-DNA phosphodiesterase 1 (TDP1) is able to eliminate TOP1cc. Thus, TDP1 interferes with the action of topotecan. Poly(ADP-ribose) polymerase 1 (PARP1) is a key regulator of many processes in the cell, such as maintaining the integrity of the genome, regulation of the cell cycle, cell death, and others. PARP1 also controls the repair of TOP1cc. We performed a transcriptomic analysis of wild type and PARP1 knockout HEK293A cells treated with topotecan and TDP1 inhibitor OL9-119 alone and in combination. The largest number of differentially expressed genes (DEGs, about 4000 both up- and down-regulated genes) was found in knockout cells. Topotecan and OL9-119 treatment elicited significantly fewer DEGs in WT cells and negligible DEGs in PARP1-KO cells. A significant part of the changes caused by PARP1-KO affected the synthesis and processing of proteins. Differences under the action of treatment with TOP1 or TDP1 inhibitors alone were found in the signaling pathways for the development of cancer, DNA repair, and the proteasome. The drug combination resulted in DEGs in the ribosome, proteasome, spliceosome, and oxidative phosphorylation pathways.

mRNA regions where 80S ribosomes pause during translation elongation in vivo interact with protein uS19, a component of the decoding site.

In eukaryotic ribosomes, the conserved protein uS19, formerly known as S15, extends with its C-terminal tail to the decoding site. The cross-linking of uS19 to the A site codon has been detected using synthetic mRNAs bearing 4-thiouridine (s4U) residues. Here, we showed that the A-site tRNA prevents this cross-linking and that the P site codon does not contact uS19. Next, we focused on determining uS19-mRNA interactions in vivo by applying the photoactivatable-ribonucleoside enhancing cross-linking and immunoprecipitation method to a stable HEK293 cell line producing FLAG-tagged uS19 and grown in a medium containing s4U. We found that when translation was stopped by cycloheximide, uS19 was efficiently cross-linked to mRNA regions with a high frequency of Glu, Lys and, more rarely, Arg codons. The results indicate that the complexes, in which the A site codon is not involved in the formation of the mRNA-tRNA duplex, are present among the cycloheximide-arrested 80S complexes, which implies pausing of elongating ribosomes at the above mRNA regions. Thus, our findings demonstrate that the human ribosomal protein uS19 interacts with mRNAs during translation elongation and highlight the regions of mRNAs where ribosome pausing occurs, bringing new structural and functional insights into eukaryotic translation in vivo.

Size and Methylation Index of Cell-Free and Cell-Surface-Bound DNA in Blood of Breast Cancer Patients in the Contest of Liquid Biopsy.

Aberrantly methylated circulating DNA (cirDNA) has proven to be a good cancer marker, but its detection is limited by low concentrations, fragmentation, and insufficiency. Since the methylated cirDNA was shown to be more stable in circulation than the unmethylated one and was shown to bind with the blood cell surface, we studied the concentration, representation, and fragmentation of tumor-derived methylated DNA in cell-free and cell-surface-associated DNA. We found that long DNA fragments (more than 10 kb) are mainly associated with the surface of blood cells. However, in plasma short DNA fragments (100-1000 bp) were also found along with long DNA fragments. Isolation of short fragments after separation of cirDNA in 6% PAGE followed by quantitative PCR (L1 element) has shown that short DNA fragments in healthy females represent 22% versus 0.5-4.4% in breast cancer patients. The methylated form of the RARß2 gene was detected only in long DNA fragments by Real-time TaqMan PCR of bisulfite-converted DNA. The methylation index of cirDNA from healthy women was estimated at 0%, 9%, and 7% in plasma, PBS-EDTA, and trypsin eluates from the surface of blood cells, respectively. The methylation index of breast cancer patients' DNA was found to be 33%, 15%, and 61% in the same fractions confirming the overrepresentation of methylated DNA in csbDNA.

Changes in the Transcriptome Caused by Mutations in the Ribosomal Protein uS10 Associated with a Predisposition to Colorectal Cancer.

A number of mutations in the RPS20 gene encoding the ribosomal protein uS10 have been found to be associated with a predisposition to hereditary non-polyposis colorectal carcinoma (CRC). We transfected HEK293T cells with constructs carrying the uS10 minigene with mutations identical to those mentioned above and examined the effects of the produced proteins on the cellular transcriptome. We showed that uS10 with mutations p.V50SfsX23 or p.L61EfsX11 cannot be incorporated into 40S ribosomal subunits, while the protein with the missense mutation p.V54L functionally replaces the respective endogenous protein in the 40S subunit assembly and the translation process. The comparison of RNA-seq data obtained from cells producing aberrant forms of uS10 with data for those producing the wild-type protein revealed overlapping sets of upregulated and downregulated differently expressed genes (DEGs) related to several pathways. Among the limited number of upregulated DEGs, there were genes directly associated with the progression of CRC, e.g., PPM1D and PIGN. Our findings indicate that the accumulation of the mutant forms of uS10 triggers a cascade of cellular events, similar to that which is triggered when the cell responds to a large number of erroneous proteins, suggesting that this may increase the risk of cancer.

Comparative Analysis of Tumor-Associated microRNAs and Tetraspanines from Exosomes of Plasma and Ascitic Fluids of Ovarian Cancer Patients.

Ovarian cancer (OC) is one of the most common and fatal types of gynecological cancer. In the early phase of OC detection, the current treatment and diagnostic methods are not efficient and sensitive enough. Therefore, it is crucial to explore the mechanisms of OC metastasis and discover valuable factors for early diagnosis of female cancers and novel therapeutic strategies for metastasis. Exosomes are known to be involved in the development, migration, and invasion of cancer cells, and their cargo could be useful for the non-invasive biopsy development. CD151- and Tspan8-positive exosomes are known to support the degradation of the extracellular matrix, and are involved in stroma remodeling, angiogenesis and cell motility, as well as the association of miR-24 and miR-101 with these processes. The objective of this study was to explore the relationship of these components of exosomal cargo, in patients with OC, to clarify the clinical significance of these markers in liquid biopsies. The levels of tetraspanins Tspan8+ and CD151+ exosomes were significantly higher in plasma exosomes of OC patients compared with healthy females (HFs). The relative levels of miR-24 and miR-101 in plasma exosomes of HFs were significantly higher than in plasma exosomes of OC patients, while the levels of these microRNAs in exosomes from plasma and ascites of ill females showed no difference. Our study revealed a strong direct correlation between the change in the ascites exosomes CD151+Tspan8+ subpopulation level and the expression levels of the ascites (R = 0.81, p < 0.05) and plasma exosomal miR-24 (R = 0.74, p < 0.05) in OC patients, which confirms the assumption that exosomal cargo act synergistically to increase cellular motility, affecting cellular processes and signaling. Bioinformatics analysis confirmed the involvement of CD151 and Tspan8 tetraspanins and genes controlled by miR-24-3p and miR-101 in signaling pathways, which are crucial for carcinogenesis, demonstrating that these tetraspanins and microRNAs are potential biomarkers for OC screening, and predictors of poor clinicopathological behavior in tumors.

Reorganization of the Landscape of Translated mRNAs in NSUN2-Deficient Cells and Specific Features of NSUN2 Target mRNAs.

The RNA cytosine C5 methyltransferase NSUN2 has a variety of RNA substrates and plays an important role in mRNA metabolism. NSUN2 binds to specific sequences enriched in exosomal mRNAs, suggesting its possible involvement in the sorting of mRNAs into exosomes. We applied the photoactivatable.4-thiouridine-enhanced cross-linking and immunoprecipitation assay involving high-throughput RNA sequencing (RNA-seq) to HEK293T cells to determine NSUN2 mRNA targets. NSUN2 cross-linking sites were found in more than one hundred relatively abundant mRNAs with a high GC content and a pronounced secondary structure. Then, utilizing RNA-seq for the total and polysome-associated mRNA from HEK293T cells with and without the knockdown of NSUN2, we identified differentially expressed genes, as well as genes with altered translational efficiency (GATEs). It turned out that the up-regulated GATE mRNAs were much shorter on average than the down-regulated ones, and their GC content was higher; moreover, they contained motifs with C residues located in GC-rich environments. Our findings reveal the specific features of mRNAs that make them potential targets for NSUN2 and expand our understanding of the role of NSUN2 in controlling translation and, possibly, in mRNA sorting into exosomes implemented through the methylation of cytosine residues.

The Cytogenetic Map of the Nile Crocodile (Crocodylus niloticus, Crocodylidae, Reptilia) with Fluorescence In Situ Localization of Major Repetitive DNAs.

Tandemly arranged and dispersed repetitive DNA sequences are important structural and functional elements that make up a significant portion of vertebrate genomes. Using high throughput, low coverage whole genome sequencing followed by bioinformatics analysis, we have identified seven major tandem repetitive DNAs and two fragments of LTR retrotransposons in the genome of the Nile crocodile (Crocodylus niloticus, 2n = 32). The repeats showed great variability in structure, genomic organization, and chromosomal distribution as revealed by fluorescence in situ hybridization (FISH). We found that centromeric and pericentromeric heterochromatin of C. niloticus is composed of previously described in Crocodylus siamensis CSI-HindIII and CSI-DraI repetitive sequence families, a satellite revealed in Crocodylus porosus, and additionally contains at least three previously unannotated tandem repeats. Both LTR sequences identified here belong to the ERV1 family of endogenous retroviruses. Each pericentromeric region was characterized by a diverse set of repeats, with the exception of chromosome pair 4, in which we found only one type of satellite. Only a few repeats showed non-centromeric signals in addition to their centromeric localization. Mapping of 18S-28S ribosomal RNA genes and telomeric sequences (TTAGGG)n did not demonstrate any co-localization of these sequences with revealed centromeric and pericentromeric heterochromatic blocks.

Flow-Seq Evaluation of Translation Driven by a Set of Natural Escherichia coli 5'-UTR of Variable Length.

Flow-seq is a method that combines fluorescently activated cell sorting and next-generation sequencing to deduce a large amount of data about translation efficiency from a single experiment. Here, we constructed a library of fluorescent protein-based reporters preceded by a set of 648 natural 5'-untranslated regions (5'-UTRs) of Escherichia coli genes. Usually, Flow-seq libraries are constructed using uniform-length sequence elements, in contrast to natural situations, where functional elements are of heterogenous lengths. Here, we demonstrated that a 5'-UTR library of variable length could be created and analyzed with Flow-seq. In line with previous Flow-seq experiments with randomized 5'-UTRs, we observed the influence of an RNA secondary structure and Shine-Dalgarno sequences on translation efficiency; however, the variability of these parameters for natural 5'-UTRs in our library was smaller in comparison with randomized libraries. In line with this, we only observed a 30-fold difference in translation efficiency between the best and worst bins sorted with this factor. The results correlated with those obtained with ribosome profiling.

Deficiency of the Ribosomal Protein uL5 Leads to Significant Rearrangements of the Transcriptional and Translational Landscapes in Mammalian Cells.

Protein uL5 (formerly called L11) is an integral component of the large (60S) subunit of the human ribosome, and its deficiency in cells leads to the impaired biogenesis of 60S subunits. Using RNA interference, we reduced the level of uL5 in HEK293T cells by three times, which caused an almost proportional decrease in the content of the fraction corresponding to 80S ribosomes, without a noticeable diminution in the level of polysomes. By RNA sequencing of uL5-deficient and control cell samples, which were those of total mRNA and mRNA from the polysome fraction, we identified hundreds of differentially expressed genes (DEGs) at the transcriptome and translatome levels and revealed dozens of genes with altered translational efficiency (GATEs). Transcriptionally up-regulated DEGs were mainly associated with rRNA processing, pre-mRNA splicing, translation and DNA repair, while down-regulated DEGs were genes of membrane proteins; the type of regulation depended on the GC content in the 3' untranslated regions of DEG mRNAs. The belonging of GATEs to up-regulated and down-regulated ones was determined by the coding sequence length of their mRNAs. Our findings suggest that the effects observed in uL5-deficient cells result from an insufficiency of translationally active ribosomes caused by a deficiency of 60S subunits.

Knockdown of the Ribosomal Protein eL38 in HEK293 Cells Changes the Translational Efficiency of Specific Genes.

The protein eL38 is one of the smallest proteins of the mammalian ribosome, which is a component of its large (60S) subunit. The haploinsufficiency of eL38 in mice leads to the Tail-short mutant phenotype characterized by defects in the development of the axial skeleton caused by the poor translation of mRNA subsets of Hox genes. Using the ribosome profiling assay applied to HEK293 cells knocked down of eL38, we examined the effects of the lack of eL38 in 60S subunits on gene expression at the level of translation. A four-fold decrease in the cell content of eL38 was shown to result in significant changes in the translational efficiencies of 150 genes. Among the genes, whose expression at the level of translation was enhanced, there were mainly those associated with basic metabolic processes; namely, translation, protein folding, chromosome organization, splicing, and others. The set of genes with reduced translation efficiencies contained those that are mostly involved in the processes related to the regulation of transcription, including the activation of Hox genes. Thus, we demonstrated that eL38 insufficiency significantly affects the expression of certain genes at the translational level. Our findings facilitate understanding the possible causes of some anomalies in eL38-deficient animals.

Human Gut Microbiome Response to Short-Term Bifidobacterium-Based Probiotic Treatment.

Gut microbiota is believed to play a crucial role in modulating obesity in humans, and probiotics affecting gut microbiota can alleviate some of the obesity-related health complications. The study was aimed to investigate changes in the composition of the gut microbiome in obese humans due to short-term (2 weeks) treatment of obese patients with a probiotic preparation containing Bifidobacterium longum. Faecal microbiome diversity was studied using the 16S amplicon sequencing by Illumina MiSeq. Bioinformatic analysis showed distribution across 14 phyla (with Firmicutes and Bacteroidetes dominating), 21 class, 125 genera and 973 OTUs. The probiotic treatment decreased relative abundance of Bacteroidetes (Prevotellaceae and Bacteroidaceae), while increasing that of Actinobacteria (Bifidobacteriaceae and Coriobacteriaceae), and Firmicutes (Negativicutes: Veillonellaceae and Clostridia: Peptostreptococcaceae). The probiotic treatment decreased total blood sugar and increased patients' assessment of their physical and mental health. Thus even the short-term Bifidobacterium-based probiotic treatment brought significant compositional changes in the 16S rRNA gene diversity in faecal bacterial assemblages by increasing beneficial and decreasing pathogenic or opportunistic bacteria; the related shifts in life quality assessment necessitate further research into the causal relationships involved.

Development and characterization of novel 2'-F-RNA aptamers specific to human total and glycated hemoglobins.

Aptamers are short DNA and RNA fragments which bind their molecular targets with affinity and specificity comparable to those of antibodies. Here, we describe the selection of novel 2'-F-RNA aptamers against total human hemoglobin or its glycated form HbA1c. After SELEX and high-throughput sequencing of the enriched libraries, affinities and specificities of candidate aptamers and their truncated variants were examined by the solid-phase bioluminescent assay. As a result, we identified aptamers specific to both hemoglobins or only glycated HbA1c. The developed 2'-F-RNA aptamers have shown their applicability for detection of total and glycated hemoglobin in one sample.

Human ribosomal protein eS1 is engaged in cellular events related to processing and functioning of U11 snRNA.

Ribosomal proteins are involved in many cellular processes through interactions with various RNAs. Here, applying the photoactivatable-ribonucleoside-enhanced cross-linking and immunoprecipitation approach to HEK293 cells overproducing ribosomal protein (rp) eS1, we determined the products of RNU5A-1 and RNU11 genes encoding U5 and U11 snRNAs as the RNA partners of ribosome-unbound rp eS1. U11 pre-snRNA-associated rp eS1 was revealed in the cytoplasm and nucleus where rp eS1-bound U11/U12 di-snRNP was also found. Utilizing recombinant rp eS1 and 4-thiouridine-containing U11 snRNA transcript, we identified an N-terminal peptide contacting the U-rich sequence in the Sm site-containing RNA region. We also showed that the rp eS1 binding site on U11 snRNA is located in the cleft between stem-loops I and III and that its structure mimics the respective site on the 18S rRNA. It was found that cell depletion of rp eS1 leads to a decrease in the splicing efficiency of minor introns and to an increase in the level of U11 pre-snRNA with the unprocessed 3' terminus. Our findings demonstrate the engagement of human rp eS1 in events related to the U11 snRNA processing and to minor-class splicing. Contacts of rp eS1 with U5 snRNA in the minor pre-catalytic spliceosome are discussed.

Recognition but no repair of abasic site in single-stranded DNA by human ribosomal uS3 protein residing within intact 40S subunit.

Isolated human ribosomal protein uS3 has extra-ribosomal functions including those related to base excision DNA repair, e.g. AP lyase activity that nicks double-stranded (ds) DNA 3΄ to the abasic (AP) site. However, the ability of uS3 residing within ribosome to recognize and cleave damaged DNA has never been addressed. Here, we compare interactions of single-stranded (ss) DNA and dsDNA bearing AP site with human ribosome-bound uS3 and with the isolated protein, whose interactions with ssDNA were not yet studied. The AP lyase activity of free uS3 was much higher with ssDNA than with dsDNA, whereas ribosome-bound uS3 was completely deprived of this activity. Nevertheless, an exposed peptide of ribosome-bound uS3 located far away from the putative catalytic center previously suggested for isolated uS3 cross-linked to full-length uncleaved ssDNA, but not to dsDNA. In contrast, free uS3 cross-linked mainly to the 5΄-part of the damaged DNA strand after its cleavage at the AP site. ChIP-seq analysis showed preferential uS3 binding to nucleolus-associated chromatin domains. We conclude that free and ribosome-bound uS3 proteins interact with AP sites differently, exhibiting their non-translational functions in DNA repair in and around the nucleolus and in regulation of DNA damage response in looped DNA structures, respectively.

Chromosomal-Level Assembly of the Asian Seabass Genome Using Long Sequence Reads and Multi-layered Scaffolding.

We report here the ~670 Mb genome assembly of the Asian seabass (Lates calcarifer), a tropical marine teleost. We used long-read sequencing augmented by transcriptomics, optical and genetic mapping along with shared synteny from closely related fish species to derive a chromosome-level assembly with a contig N50 size over 1 Mb and scaffold N50 size over 25 Mb that span ~90% of the genome. The population structure of L. calcarifer species complex was analyzed by re-sequencing 61 individuals representing various regions across the species' native range. SNP analyses identified high levels of genetic diversity and confirmed earlier indications of a population stratification comprising three clades with signs of admixture apparent in the South-East Asian population. The quality of the Asian seabass genome assembly far exceeds that of any other fish species, and will serve as a new standard for fish genomics.

Correction: Chromosomal-Level Assembly of the Asian Seabass Genome Using Long Sequence Reads and Multi-layered Scaffolding.

[This corrects the article DOI: 10.1371/journal.pgen.1005954.].

Cytosolic YB-1 and NSUN2 are the only proteins recognizing specific motifs present in mRNAs enriched in exosomes.

Exosomes, membranous vesicles secreted by various cells, are involved in intercellular communication and carry vast repertoires of RNAs and proteins. Processes mediating RNA sorting into exosomes are currently poorly understood. Using bioinformatics approaches, three structural motifs ACCAGCCU, CAGUGAGC and UAAUCCCA have been discovered as enriched in exosomal mRNAs and long noncoding RNAs. Here, utilizing short RNA hairpins, each containing one of the motifs, in a pull-down assay of cytosolic extract of human embryonic kidney 293 (HEK293) cells, we prove that multifunctional RNA-binding protein YB-1 specifically interacts with all three motifs, whereas methyltransferase NSUN2 recognizes only the motif CAGUGAGC. RNA hairpins other than those mentioned above pull out neither YB-1 nor NSUN2. Both these proteins are found in exosomes secreted by HEK293 cells. YB-1 for all that is detected as a form having a slightly higher electrophoretic mobility than that of YB-1 associated with the above RNA hairpins, assuming changes in posttranslational modifications of the protein during its transfer from cytoplasm into exosomes. Next generation sequencing of total exosomal RNA (eRNA) reveals a large representative set of RNA species, including mRNAs containing the above-mentioned motifs. The degree of enrichment in exosomes with this kind of mRNAs strongly depends on the locations of eRNA-specific motifs within the mRNA sequences. Altogether, our findings point to YB-1 and NSUN2 as possible mediators of the process of transfer of specific mRNAs into exosomes, allowing us to speculate on an involvement of these proteins in the mRNA sorting via the recognition of the above motifs.