Structural variation across 138,134 samples in the TOPMed consortium.

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hemotologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

Structural variation across 138,134 samples in the TOPMed consortium.

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hematologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

xAtlas: scalable small variant calling across heterogeneous next-generation sequencing experiments.

BACKGROUND: The growing volume and heterogeneity of next-generation sequencing (NGS) data complicate the further optimization of identifying DNA variation, especially considering that curated high-confidence variant call sets frequently used to validate these methods are generally developed from the analysis of comparatively small and homogeneous sample sets. FINDINGS: We have developed xAtlas, a single-sample variant caller for single-nucleotide variants (SNVs) and small insertions and deletions (indels) in NGS data. xAtlas features rapid runtimes, support for CRAM and gVCF file formats, and retraining capabilities. xAtlas reports SNVs with 99.11% recall and 98.43% precision across a reference HG002 sample at 60× whole-genome coverage in less than 2 CPU hours. Applying xAtlas to 3,202 samples at 30× whole-genome coverage from the 1000 Genomes Project achieves an average runtime of 1.7 hours per sample and a clear separation of the individual populations in principal component analysis across called SNVs. CONCLUSIONS: xAtlas is a fast, lightweight, and accurate SNV and small indel calling method. Source code for xAtlas is available under a BSD 3-clause license at https://github.com/jfarek/xatlas.

Advancing human genetics research and drug discovery through exome sequencing of the UK Biobank.

The UK Biobank Exome Sequencing Consortium (UKB-ESC) is a private-public partnership between the UK Biobank (UKB) and eight biopharmaceutical companies that will complete the sequencing of exomes for all ~500,000 UKB participants. Here, we describe the early results from ~200,000 UKB participants and the features of this project that enabled its success. The biopharmaceutical industry has increasingly used human genetics to improve success in drug discovery. Recognizing the need for large-scale human genetics data, as well as the unique value of the data access and contribution terms of the UKB, the UKB-ESC was formed. As a result, exome data from 200,643 UKB enrollees are now available. These data include ~10 million exonic variants-a rich resource of rare coding variation that is particularly valuable for drug discovery. The UKB-ESC precompetitive collaboration has further strengthened academic and industry ties and has provided teams with an opportunity to interact with and learn from the wider research community.

Amine-to-Azide Conversion on Native RNA via Metal-Free Diazotransfer Opens New Avenues for RNA Manipulations.

A major challenge in the field of RNA chemistry is the identification of selective and quantitative conversion reactions on RNA that can be used for tagging and any other RNA tool development. Here, we introduce metal-free diazotransfer on native RNA containing an aliphatic primary amino group using the diazotizing reagent fluorosulfuryl azide (FSO2 N3 ). The reaction provides the corresponding azide-modified RNA in nearly quantitatively yields without affecting the nucleobase amino groups. The obtained azido-RNA can then be further processed utilizing well-established bioortho-gonal reactions, such as azide-alkyne cycloadditions (Click) or Staudinger ligations. We exemplify the robustness of this approach for the synthesis of peptidyl-tRNA mimics and for the pull-down of 3-(3-amino-3-carboxypropyl)uridine (acp3 U)- and lysidine (k2 C)-containing tRNAs of an Escherichia coli tRNA pool isolated from cellular extracts. Our approach therefore adds a new dimension to the targeted chemical manipulation of diverse RNA species.

Amine-to-Azide Conversion on Native RNA via Metal-Free Diazotransfer Opens New Avenues for RNA Manipulations.

A major challenge in the field of RNA chemistry is the identification of selective and quantitative conversion reactions on RNA that can be used for tagging and any other RNA tool development. Here, we introduce metal-free diazotransfer on native RNA containing an aliphatic primary amino group using the diazotizing reagent fluorosulfuryl azide (FSO2N3). The reaction provides the corresponding azide-modified RNA in nearly quantitatively yields without affecting the nucleobase amino groups. The obtained azido-RNA can then be further processed utilizing well-established bioortho-gonal reactions, such as azide-alkyne cycloadditions (Click) or Staudinger ligations. We exemplify the robustness of this approach for the synthesis of peptidyl-tRNA mimics and for the pull-down of 3-(3-amino-3-carboxypropyl)uridine (acp3U)- and lysidine (k2C)-containing tRNAs of an Escherichia coli tRNA pool isolated from cellular extracts. Our approach therefore adds a new dimension to the targeted chemical manipulation of diverse RNA species.

Parliament2: Accurate structural variant calling at scale.

BACKGROUND: Structural variants (SVs) are critical contributors to genetic diversity and genomic disease. To predict the phenotypic impact of SVs, there is a need for better estimates of both the occurrence and frequency of SVs, preferably from large, ethnically diverse cohorts. Thus, the current standard approach requires the use of short paired-end reads, which remain challenging to detect, especially at the scale of hundreds to thousands of samples. FINDINGS: We present Parliament2, a consensus SV framework that leverages multiple best-in-class methods to identify high-quality SVs from short-read DNA sequence data at scale. Parliament2 incorporates pre-installed SV callers that are optimized for efficient execution in parallel to reduce the overall runtime and costs. We demonstrate the accuracy of Parliament2 when applied to data from NovaSeq and HiSeq X platforms with the Genome in a Bottle (GIAB) SV call set across all size classes. The reported quality score per SV is calibrated across different SV types and size classes. Parliament2 has the highest F1 score (74.27%) measured across the independent gold standard from GIAB. We illustrate the compute performance by processing all 1000 Genomes samples (2,691 samples) in <1 day on GRCH38. Parliament2 improves the runtime performance of individual methods and is open source (https://github.com/slzarate/parliament2), and a Docker image, as well as a WDL implementation, is available. CONCLUSION: Parliament2 provides both a highly accurate single-sample SV call set from short-read DNA sequence data and enables cost-efficient application over cloud or cluster environments, processing thousands of samples.

Structural distinctions between NAD+ riboswitch domains 1 and 2 determine differential folding and ligand binding.

Riboswitches are important gene regulatory elements frequently encountered in bacterial mRNAs. The recently discovered nadA riboswitch contains two similar, tandemly arrayed aptamer domains, with the first domain possessing high affinity for nicotinamide adenine dinucleotide (NAD+). The second domain which comprises the ribosomal binding site in a putative regulatory helix, however, has withdrawn from detection of ligand-induced structural modulation thus far, and therefore, the identity of the cognate ligand and the regulation mechanism have remained unclear. Here, we report crystal structures of both riboswitch domains, each bound to NAD+. Furthermore, we demonstrate that ligand binding to domain 2 requires significantly higher concentrations of NAD+ (or ADP retaining analogs) compared to domain 1. Using a fluorescence spectroscopic approach, we further shed light on the structural features which are responsible for the different ligand affinities, and describe the Mg2+-dependent, distinct folding and pre-organization of their binding pockets. Finally, we speculate about possible scenarios for nadA RNA gene regulation as a putative two-concentration sensor module for a time-controlled signal that is primed and stalled by the gene regulation machinery at low ligand concentrations (domain 1), and finally triggers repression of translation as soon as high ligand concentrations are reached in the cell (domain 2).

Postsynthetic On-Column 2' Functionalization of RNA by Convenient Versatile Method.

We report a universal straightforward strategy for the chemical synthesis of modified oligoribonucleotides containing functional groups of different structures at the 2' position of ribose. The on-column synthetic concept is based on the incorporation of two types of commercial nucleotide phosphoramidites containing orthogonal 2'-O-protecting groups, namely 2'-O-thiomorpholine-carbothioate (TC, as "permanent") and 2'-O-tert-butyl(dimethyl)silyl (tBDMS, as "temporary"), to RNA during solid-phase synthesis. Subsequently, the support-bound RNA undergoes selective deprotection and follows postsynthetic 2' functionalization of the naked hydroxyl group. This convenient method to tailor RNA, utilizing the advantages of solid phase approaches, gives an opportunity to introduce site-specifically a wide range of linkers and functional groups. By this strategy, a series of RNAs containing diverse 2' functionalities were synthesized and studied with respect to their physicochemical properties.

Crucial Roles of Two Hydrated Mg2+ Ions in Reaction Catalysis of the Pistol Ribozyme.

Pistol ribozymes constitute a new class of small self-cleaving RNAs. Crystal structures have been solved, providing three-dimensional snapshots along the reaction coordinate of pistol phosphodiester cleavage, corresponding to the pre-catalytic state, a vanadate mimic of the transition state, and the product. The results led to the proposed underlying chemical mechanism. Importantly, a hydrated Mg2+ ion remains innersphere-coordinated to N7 of G33 in all three states, and is consistent with its likely role as acid in general acid base catalysis (δ and ß catalysis). Strikingly, the new structures shed light on a second hydrated Mg2+ ion that approaches the scissile phosphate from its binding site in the pre-cleavage state to reach out for water-mediated hydrogen bonding in the cyclophosphate product. The major role of the second Mg2+ ion appears to be the stabilization of product conformation. This study delivers a mechanistic understanding of ribozyme-catalyzed backbone cleavage.

SAM-VI riboswitch structure and signature for ligand discrimination.

Riboswitches are metabolite-sensing, conserved domains located in non-coding regions of mRNA that are central to regulation of gene expression. Here we report the first three-dimensional structure of the recently discovered S-adenosyl-L-methionine responsive SAM-VI riboswitch. SAM-VI adopts a unique fold and ligand pocket that are distinct from all other known SAM riboswitch classes. The ligand binds to the junctional region with its adenine tightly intercalated and Hoogsteen base-paired. Furthermore, we reveal the ligand discrimination mode of SAM-VI by additional X-ray structures of this riboswitch bound to S-adenosyl-L-homocysteine and a synthetic ligand mimic, in combination with isothermal titration calorimetry and fluorescence spectroscopy to explore binding thermodynamics and kinetics. The structure is further evaluated by analysis of ligand binding to SAM-VI mutants. It thus provides a thorough basis for developing synthetic SAM cofactors for applications in chemical and synthetic RNA biology.

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.

Complexes of Pro-Apoptotic siRNAs and Carbosilane Dendrimers: Formation and Effect on Cancer Cells.

This paper examines the complexation of anti-cancer small interfering RNAs (siRNAs) by cationic carbosilane dendrimers, and the interaction of the formed complexes with HeLa and HL-60 cancer cells. Stepwise formation of the complexes accompanied by the evolution of their properties has been observed through the increase of the charge ratio (dendrimer/siRNA). The complexes decrease the viability of both "easy-to-transfect" cells (HeLa) and "hard-to transfect" ones (HL-60), indicating a high potential of the cationic carbosilane dendrimers for siRNA delivery into tumor cells.

Functional equivalence of genome sequencing analysis pipelines enables harmonized variant calling across human genetics projects.

Hundreds of thousands of human whole genome sequencing (WGS) datasets will be generated over the next few years. These data are more valuable in aggregate: joint analysis of genomes from many sources increases sample size and statistical power. A central challenge for joint analysis is that different WGS data processing pipelines cause substantial differences in variant calling in combined datasets, necessitating computationally expensive reprocessing. This approach is no longer tenable given the scale of current studies and data volumes. Here, we define WGS data processing standards that allow different groups to produce functionally equivalent (FE) results, yet still innovate on data processing pipelines. We present initial FE pipelines developed at five genome centers and show that they yield similar variant calling results and produce significantly less variability than sequencing replicates. This work alleviates a key technical bottleneck for genome aggregation and helps lay the foundation for community-wide human genetics studies.

Structural features of the interaction of the 3'-untranslated region of mRNA containing exosomal RNA-specific motifs with YB-1, a potential mediator of mRNA sorting.

We have previously shown that YB-1 is the only protein of the HEK293 cell cytoplasmic (S100) extract that specifically interacts with RNA hairpins each containing one of the motifs ACCAGCCU (1), CAGUGAGC (2) and UAAUCCCA (3), which had been identified as often found in exosomal RNA and proposed as potential cis-acting elements targeting RNAs into exosomes. Here we explored the interactions of YB-1 with a fragment of the 3'-untranslated region (UTR) of septin 14 mRNA (SEPT14 RNA), which contains all three motifs. We demonstrated the occurrence of YB-1 among proteins pulled down from the HEK293 S100 extract using biotinylated SEPT14 RNA. With recombinant YB-1, it was found that SEPT14 RNA can bind up to 5 moles of protein per mole of RNA in a cooperative manner, which was shown to be mainly facilitated by the presence of the above motifs. RNA hairpins with motifs 1 and 2 competed with SEPT14 RNA for binding to the protein, whereas that with motif 3 was less competitive, in accordance with the affinity of YB-1 for these RNA hairpins. With YB-1-bound RNA, nucleotides protected from attack by hydroxyl radicals were revealed in all three motifs, although hairpins with motif 2 and especially with motif 1 contained many protected nucleotides outside the motifs, suggesting that the specific environments of these motifs contribute significantly to the YB-1 binding. An analysis of the environments of motifs 1-3 in the HEK293 cell mRNA 3' UTRs gained from RNA-seq data led us to conclude that the primary binding sites of YB-1 in the 3' UTRs are hairpins containing some part of the motif along with its specific surroundings; the consensus sequences of these hairpins were derived. Thus, our findings provide a new understanding of the structural basis of the interactions between YB-1 and mRNAs carrying the aforementioned motifs.

Recent Advances in Nucleic Acid Targeting Probes and Supramolecular Constructs Based on Pyrene-Modified Oligonucleotides.

In this review, we summarize the recent advances in the use of pyrene-modified oligonucleotides as a platform for functional nucleic acid-based constructs. Pyrene is of special interest for the development of nucleic acid-based tools due to its unique fluorescent properties (sensitivity of fluorescence to the microenvironment, ability to form excimers and exciplexes, long fluorescence lifetime, high quantum yield), ability to intercalate into the nucleic acid duplex, to act as a π-π-stacking (including anchoring) moiety, and others. These properties of pyrene have been used to construct novel sensitive fluorescent probes for the sequence-specific detection of nucleic acids and the discrimination of single nucleotide polymorphisms (SNPs), aptamer-based biosensors, agents for binding of double-stranded DNAs, and building blocks for supramolecular complexes. Special attention is paid to the influence of the design of pyrene-modified oligonucleotides on their properties, i.e., the structure-function relationships. The perspectives for the applications of pyrene-modified oligonucleotides in biomolecular studies, diagnostics, and nanotechnology are discussed.

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.

Anticancer siRNA cocktails as a novel tool to treat cancer cells. Part (A). Mechanisms of interaction.

This paper examines a perspective on the use of newly engineered nanomaterials as effective and safe carriers of genes for the therapy of cancer. Three different groups of cationic dendrimers (PAMAM, phosphorus and carbosilane) were complexed with anticancer siRNA and their biophysical properties of the dendriplexes analyzed. The potential of the dendrimers as nanocarriers for anticancer siBcl-xl, siBcl-2, siMcl-1 siRNAs and a siScrambled sequence was explored. Dendrimer/siRNA complexes were characterized by methods including fluorescence, zeta potential, dynamic light scattering, circular dichroism, gel electrophoresis and transmission electron microscopy. Some of the experiments were done with heparin to check if siRNA can be easily disassociated from the complexes, and whether released siRNA maintains its structure after interaction with the dendrimer. The results indicate that siRNAs form complexes with all the dendrimers tested. Oligoribonucleotide duplexes can be released from dendriplexes after heparin treatment and the structure of siRNA is maintained in the case of PAMAM or carbosilane dendrimers. The dendrimers were also effective in protecting siRNA from RNase A activity. The selection of the best siRNA carrier will be made based on cell culture studies (Part B).

2'-Bispyrene-modified 2'-O-methyl RNA probes as useful tools for the detection of RNA: synthesis, fluorescent properties, and duplex stability.

The synthesis and properties two series of new 2'-O-methyl RNA probes, each containing a single insertion of a 2'-bispyrenylmethylphosphorodiamidate derivative of a nucleotide (U, C, A, and G), are described. As demonstrated by UV melting studies, the probes form stable complexes with model RNAs and DNAs. Significant increases (up to 21-fold) in pyrene excimer fluorescence intensity were observed upon binding of most of the probes with complementary RNAs, but not with DNAs. The fluorescence spectra are independent of the nature of the modified nucleotides. The nucleotides on the 5'-side of the modified nucleotide have no effect on the fluorescence spectra, whereas the natures of the two nucleotides on the 3'-side are important: CC, CG, and UC dinucleotide units on the 3'-side of the modified nucleotide provide the maximum increases in excimer fluorescence intensity. This study suggests that these 2'-bispyrene-labeled 2'-O-methyl RNA probes might be useful tools for detection of RNAs.

Fluorescent probes for nucleic Acid visualization in fixed and live cells.

This review analyses the literature concerning non-fluorescent and fluorescent probes for nucleic acid imaging in fixed and living cells from the point of view of their suitability for imaging intracellular native RNA and DNA. Attention is mainly paid to fluorescent probes for fluorescence microscopy imaging. Requirements for the target-binding part and the fluorophore making up the probe are formulated. In the case of native double-stranded DNA, structure-specific and sequence-specific probes are discussed. Among the latest, three classes of dsDNA-targeting molecules are described: (i) sequence-specific peptides and proteins; (ii) triplex-forming oligonucleotides and (iii) polyamide oligo(N-methylpyrrole/N-methylimidazole) minor groove binders. Polyamides seem to be the most promising targeting agents for fluorescent probe design, however, some technical problems remain to be solved, such as the relatively low sequence specificity and the high background fluorescence inside the cells. Several examples of fluorescent probe applications for DNA imaging in fixed and living cells are cited. In the case of intracellular RNA, only modified oligonucleotides can provide such sequence-specific imaging. Several approaches for designing fluorescent probes are considered: linear fluorescent probes based on modified oligonucleotide analogs, molecular beacons, binary fluorescent probes and template-directed reactions with fluorescence probe formation, FRET donor-acceptor pairs, pyrene excimers, aptamers and others. The suitability of all these methods for living cell applications is discussed.