The human noncoding genome defined by genetic diversity.

Understanding the significance of genetic variants in the noncoding genome is emerging as the next challenge in human genomics. We used the power of 11,257 whole-genome sequences and 16,384 heptamers (7-nt motifs) to build a map of sequence constraint for the human species. This build differed substantially from traditional maps of interspecies conservation and identified regulatory elements among the most constrained regions of the genome. Using new Hi-C experimental data, we describe a strong pattern of coordination over 2 Mb where the most constrained regulatory elements associate with the most essential genes. Constrained regions of the noncoding genome are up to 52-fold enriched for known pathogenic variants as compared to unconstrained regions (21-fold when compared to the genome average). This map of sequence constraint across thousands of individuals is an asset to help interpret noncoding elements in the human genome, prioritize variants and reconsider gene units at a larger scale.

Profiling of Short-Tandem-Repeat Disease Alleles in 12,632 Human Whole Genomes.

Short tandem repeats (STRs) are hyper-mutable sequences in the human genome. They are often used in forensics and population genetics and are also the underlying cause of many genetic diseases. There are challenges associated with accurately determining the length polymorphism of STR loci in the genome by next-generation sequencing (NGS). In particular, accurate detection of pathological STR expansion is limited by the sequence read length during whole-genome analysis. We developed TREDPARSE, a software package that incorporates various cues from read alignment and paired-end distance distribution, as well as a sequence stutter model, in a probabilistic framework to infer repeat sizes for genetic loci, and we used this software to infer repeat sizes for 30 known disease loci. Using simulated data, we show that TREDPARSE outperforms other available software. We sampled the full genome sequences of 12,632 individuals to an average read depth of approximately 30× to 40× with Illumina HiSeq X. We identified 138 individuals with risk alleles at 15 STR disease loci. We validated a representative subset of the samples (n = 19) by Sanger and by Oxford Nanopore sequencing. Additionally, we validated the STR calls against known allele sizes in a set of GeT-RM reference cell-line materials (n = 6). Several STR loci that are entirely guanine or cytosines (G or C) have insufficient read evidence for inference and therefore could not be assayed precisely by TREDPARSE. TREDPARSE extends the limit of STR size detection beyond the physical sequence read length. This extension is critical because many of the disease risk cutoffs are close to or beyond the short sequence read length of 100 to 150 bases.

Deep sequencing of 10,000 human genomes.

We report on the sequencing of 10,545 human genomes at 30×-40× coverage with an emphasis on quality metrics and novel variant and sequence discovery. We find that 84% of an individual human genome can be sequenced confidently. This high-confidence region includes 91.5% of exon sequence and 95.2% of known pathogenic variant positions. We present the distribution of over 150 million single-nucleotide variants in the coding and noncoding genome. Each newly sequenced genome contributes an average of 8,579 novel variants. In addition, each genome carries on average 0.7 Mb of sequence that is not found in the main build of the hg38 reference genome. The density of this catalog of variation allowed us to construct high-resolution profiles that define genomic sites that are highly intolerant of genetic variation. These results indicate that the data generated by deep genome sequencing is of the quality necessary for clinical use.

Library preparation methodology can influence genomic and functional predictions in human microbiome research.

Observations from human microbiome studies are often conflicting or inconclusive. Many factors likely contribute to these issues including small cohort sizes, sample collection, and handling and processing differences. The field of microbiome research is moving from 16S rDNA gene sequencing to a more comprehensive genomic and functional representation through whole-genome sequencing (WGS) of complete communities. Here we performed quantitative and qualitative analyses comparing WGS metagenomic data from human stool specimens using the Illumina Nextera XT and Illumina TruSeq DNA PCR-free kits, and the KAPA Biosystems Hyper Prep PCR and PCR-free systems. Significant differences in taxonomy are observed among the four different next-generation sequencing library preparations using a DNA mock community and a cell control of known concentration. We also revealed biases in error profiles, duplication rates, and loss of reads representing organisms that have a high %G+C content that can significantly impact results. As with all methods, the use of benchmarking controls has revealed critical differences among methods that impact sequencing results and later would impact study interpretation. We recommend that the community adopt PCR-free-based approaches to reduce PCR bias that affects calculations of abundance and to improve assemblies for accurate taxonomic assignment. Furthermore, the inclusion of a known-input cell spike-in control provides accurate quantitation of organisms in clinical samples.

Chemical structure requirements and cellular targeting of microRNA-122 by peptide nucleic acids anti-miRs.

Anti-miRs are oligonucleotide inhibitors complementary to miRNAs that have been used extensively as tools to gain understanding of specific miRNA functions and as potential therapeutics. We showed previously that peptide nucleic acid (PNA) anti-miRs containing a few attached Lys residues were potent miRNA inhibitors. Using miR-122 as an example, we report here the PNA sequence and attached amino acid requirements for efficient miRNA targeting and show that anti-miR activity is enhanced substantially by the presence of a terminal-free thiol group, such as a Cys residue, primarily due to better cellular uptake. We show that anti-miR activity of a Cys-containing PNA is achieved by cell uptake through both clathrin-dependent and independent routes. With the aid of two PNA analogues having intrinsic fluorescence, thiazole orange (TO)-PNA and [bis-o-(aminoethoxy)phenyl]pyrrolocytosine (BoPhpC)-PNA, we explored the subcellular localization of PNA anti-miRs and our data suggest that anti-miR targeting of miR-122 may take place in or associated with endosomal compartments. Our findings are valuable for further design of PNAs and other oligonucleotides as potent anti-miR agents.

Site-selective artificial ribonucleases: oligonucleotide conjugates containing multiple imidazole residues in the catalytic domain.

Design of site-selective artificial ribonucleases (aRNases) is one of the most challenging tasks in RNA targeting. Here, we designed and studied oligonucleotide-based aRNases containing multiple imidazole residues in the catalytic part and systematically varied structure of cleaving constructs. We demonstrated that the ribonuclease activity of the conjugates is strongly affected by the number of imidazole residues in the catalytic part, the length of a linker between the catalytic imidazole groups of the construct and the oligonucleotide, and the type of anchor group, connecting linker structure and the oligonucleotide. Molecular modeling of the most active aRNases showed that preferable orientation(s) of cleaving constructs strongly depend on the structure of the anchor group and length of the linker. The inclusion of deoxyribothymidine anchor group significantly reduced the probability of cleaving groups to locate near the cleavage site, presumably due to a stacking interaction with the neighbouring nucleotide residue. Altogether the obtained results show that dynamics factors play an important role in site-specific RNA cleavage. Remarkably high cleavage activity was displayed by the conjugates with the most flexible and extended cleaving construct, which presumably provides a better opportunity for imidazole residues to be correctly positioned in the vicinity of scissile phosphodiester bond.

MicroRNA fate upon targeting with anti-miRNA oligonucleotides as revealed by an improved Northern-blot-based method for miRNA detection.

MicroRNAs (miRNAs) are small non-coding RNAs involved in fine-tuning of gene regulation. Antisense oligonucleotides (ONs) are promising tools as anti-miRNA (anti-miR) agents toward therapeutic applications and to uncover miRNA function. Such anti-miR ONs include 2'-O-methyl (OMe), cationic peptide nucleic acids like K-PNA-K3, and locked nucleic acid (LNA)-based anti-miRs such as LNA/DNA or LNA/OMe. Northern blotting is a widely used and robust technique to detect miRNAs. However, miRNA quantification in the presence of anti-miR ONs has proved to be challenging, due to detection artifacts, which has led to poor understanding of miRNA fate upon anti-miR binding. Here we show that anti-miR ON bound to miR-122 can prevent the miRNA from being properly precipitated into the purified RNA fraction using the standard RNA extraction protocol (TRI-Reagent), yielding an RNA extract that does not reflect the real cellular levels of the miRNA. An increase in the numbers of equivalents of isopropanol during the precipitation step leads to full recovery of the targeted miRNA back into the purified RNA extract. Following our improved protocol, we demonstrate by Northern blotting, in conjunction with a PNA decoy strategy and use of high denaturing PAGE, that high-affinity anti-miRs (K-PNA-K3, LNA/DNA, and LNA/OMe) sequester miR-122 without causing miRNA degradation, while miR-122 targeting with a lower-affinity anti-miR (OMe) seems to promote degradation of the miRNA. The technical issues explored in this work will have relevance for other hybridization-based techniques for miRNA quantification in the presence of anti-miR ONs.

miR-132 regulates antiviral innate immunity through suppression of the p300 transcriptional co-activator.

MicroRNAs are small, non-coding RNAs that negatively regulate gene expression. It has been proposed that microRNAs could function in the regulation of innate immunity, but this has not been demonstrated for viral infection. Here we test this hypothesis using the human pathogenic virus Kaposi's sarcoma-associated herpesvirus (KSHV) and one of its putative natural cellular targets, primary lymphatic endothelial cells (LECs). We show that an early antiviral microRNA response (6 h post-infection) includes expression of microRNAs that enhance viral gene expression. In particular, the CREB-induced miR-132 microRNA is highly upregulated after infection and has a negative effect on the expression of interferon-stimulated genes, facilitating viral replication. We show a similar function for miR-132 during infection of monocytes with herpes simplex virus-1 (HSV-1) and human cytomegalovirus (HCMV). miR-132 regulates innate antiviral immunity by inhibiting expression of the p300 transcriptional co-activator. p300 is downregulated early after KSHV infection, and inhibition of miR-132 induction restores p300 expression. Furthermore, p300 regulates miR-132 levels, revealing a dynamic equilibrium between miR-132 and p300. By targeting p300, rather than a transcription factor or signalling protein, miR-132 has a broad role in the regulation of antiviral immunity.

Efficient inhibition of miR-155 function in vivo by peptide nucleic acids.

MicroRNAs (miRNAs) play an important role in diverse physiological processes and are potential therapeutic agents. Synthetic oligonucleotides (ONs) of different chemistries have proven successful for blocking miRNA expression. However, their specificity and efficiency have not been fully evaluated. Here, we show that peptide nucleic acids (PNAs) efficiently block a key inducible miRNA expressed in the haematopoietic system, miR-155, in cultured B cells as well as in mice. Remarkably, miR-155 inhibition by PNA in primary B cells was achieved in the absence of any transfection agent. In mice, the high efficiency of the treatment was demonstrated by a strong overlap in global gene expression between B cells isolated from anti-miR-155 PNA-treated and miR-155-deficient mice. Interestingly, PNA also induced additional changes in gene expression. Our analysis provides a useful platform to aid the design of efficient and specific anti-miRNA ONs for in vivo use.

PNA-peptide conjugates as intracellular gene control agents.

Serum-stabilized PNA-internalization peptides (Pip) conjugated to PNA complementary to the 705 aberrant beta-globin splice site are able to correct splicing and increase luciferase production in Hela pLuc705 cells with sub microM EC(50) in the absence of a transfection agent. Inhibition of microRNA-122 in liver cells is achieved by treatment with complementary PNA containing just a few attached Lys residues, again without need of a transfection agent.

miR-122 targeting with LNA/2'-O-methyl oligonucleotide mixmers, peptide nucleic acids (PNA), and PNA-peptide conjugates.

MicroRNAs are small noncoding RNAs that regulate many cellular processes in a post-transcriptional mode. MicroRNA knockdown by antisense oligonucleotides is a useful strategy to explore microRNA functionality and as potential therapeutics. MicroRNA-122 (miR-122) is a liver-specific microRNA, the main function of which has been linked with lipid metabolism and liver homeostasis. Here, we show that lipofection of an antisense oligonucleotide based on a Locked Nucleic Acids (LNA)/2'-O-methyl mixmer or electroporation of a Peptide Nucleic Acid (PNA) oligomer is effective at blocking miR-122 activity in human and rat liver cells. These oligonucleotide analogs, evaluated for the first time in microRNA inhibition, are more effective than standard 2'-O-methyl oligonucleotides in binding and inhibiting microRNA action. We also show that microRNA inhibition can be achieved without the need for transfection or electroporation of the human or rat cell lines, by conjugation of an antisense PNA to the cell-penetrating peptide R6-Penetratin, or merely by linkage to just four Lys residues, highlighting the potential of PNA for future therapeutic applications as well as for studying microRNA function.

RNA targeting with peptide conjugates of oligonucleotides, siRNA and PNA.

Towards the development of oligonucleotide analogues and siRNA as drugs, one potential alternative to the use of liposomal transfection agents is the covalent conjugation of a cell-penetrating peptide (CPP), with the intention of imparting on the oligonucleotide or siRNA an enhanced ability to enter mammalian cells and reach the appropriate RNA target. We have developed robust methods for the chemical synthesis of disulfide-linked conjugates of oligonucleotide analogues, siRNA and peptide nucleic acids (PNA) with a range of cationic and other CPPs. In a HeLa cell assay with integrated plasmid reporters of Tat-dependent trans-activation at the TAR RNA target in the cell nucleus, we were unable to obtain steric block inhibition of gene expression for conjugates of CPPs with a 12-mer oligonucleotide mixmer of 2'-O-methyl and locked nucleic acids units. By contrast, we were able to obtain some reductions in expression of P38alpha MAP kinase mRNA in HeLa cells using microM concentrations of Penetratin or Tat peptides conjugated to the 3'-end of the sense strand of siRNA. However, the most promising results to date have been with a 16-mer PNA conjugated to the CPP Transportan or a double CPP R(6)-Penetratin, where we have demonstrated Tat-dependent trans-activation inhibition in HeLa cells. Results to date suggest the possibility of development of CPP-PNA conjugates as anti-HIV agents as well as other potential applications involving nuclear cell delivery, such as the redirection of splicing.

Oligonucleotide analogs as antiviral agents.

The recent availability of a wider range of oligonucleotide analogs has stimulated renewed interest in their application as potential antiviral agents through a variety of mechanisms. These mechanisms include RNase H-mediated antisense, steric block antisense and small interfering RNA targeting viral RNAs, but also other mechanisms, including blockage of virus uptake by cells and the stimulation of a Toll-like receptor-9-dependent immune response by CpG oligonucleotides.

Targeting the HIV-1 RNA leader sequence with synthetic oligonucleotides and siRNA: chemistry and cell delivery.

New candidates for development as potential drugs or virucides against HIV-1 infection and AIDS continue to be needed. The HIV-1 RNA leader sequence has many essential functional sites for virus replication and regulation that includes several highly conserved sequences. The review describes the historical context of targeting the HIV-1 RNA leader sequence with antisense phosphorothioate oligonucleotides, such as GEM 91, and goes on to describe modern approaches to targeting this region with steric blocking oligonucleotide analogues having newer and more advantageous chemistries, as well as recent studies on siRNA, towards the attainment of antiviral activity. Recent attempts to obtain improved cell delivery are highlighted, including exciting new developments in the use of peptide conjugates of peptide nucleic acid (PNA) as potential virucides.

Sequence-specific artificial ribonucleases. I. Bis-imidazole-containing oligonucleotide conjugates prepared using precursor-based strategy.

Antisense oligonucleotide conjugates, bearing constructs with two imidazole residues, were synthesized using a precursor-based technique employing post-synthetic histamine functionalization of oligonucleotides bearing methoxyoxalamido precursors at the 5'-termini. The conjugates were assessed in terms of their cleavage activities using both biochemical assays and conformational analysis by molecular modelling. The oligonucleotide part of the conjugates was complementary to the T-arm of yeast tRNA(Phe) (44-60 nt) and was expected to deliver imidazole groups near the fragile sequence C61-ACA-G65 of the tRNA. The conjugates showed ribonuclease activity at neutral pH and physiological temperature resulting in complete cleavage of the target RNA, mainly at the C63-A64 phosphodiester bond. For some constructs, cleavage was completed within 1-2 h under optimal conditions. Molecular modelling was used to determine the preferred orientation(s) of the cleaving group(s) in the complexes of the conjugates with RNA target. Cleaving constructs bearing two imidazole residues were found to be conformationally highly flexible, adopting no preferred specific conformation. No interactions other than complementary base pairing between the conjugates and the target were found to be the factors stabilizing the 'active' cleaving conformation(s).