We present a discriminative learning method for pattern discovery of binding sites in nucleic acid sequences based on hidden Markov models. Sets of positive and negative example sequences are mined for sequence motifs whose occurrence frequency varies between the sets. The method offers several objective functions, but we concentrate on mutual information of condition and motif occurrence. We perform a systematic comparison of our method and numerous published motif-finding tools. Our method achieves the highest motif discovery performance, while being faster than most published methods. We present case studies of data from various technologies, including ChIP-Seq, RIP-Chip and PAR-CLIP, of embryonic stem cell transcription factors and of RNA-binding proteins, demonstrating practicality and utility of the method. For the alternative splicing factor RBM10, our analysis finds motifs known to be splicing-relevant. The motif discovery method is implemented in the free software package Discrover. It is applicable to genome- and transcriptome-scale data, makes use of available repeat experiments and aside from binary contrasts also more complex data configurations can be utilized.

Exosomes participate in intercellular communication by transporting functionally active molecules. Such cargo from the original cells comprising proteins, micro-RNA, mRNA, single-stranded (ssDNA) and double-stranded DNA (dsDNA) molecules pleiotropically transforms the target cells. Although cancer cells secrete exosomes carrying a significant level of DNA capable of modulating oncogene expression in a recipient cell, the regulatory mechanism is unknown. We have previously reported that cancer cells produce exosomes containing NANOGP8 DNA. NANOGP8 is an oncogenic paralog of embryonic stem cell transcription factor NANOG and does not express in cells since it is a pseudogene. However, in this study, we evaluated NANOGP8 expression in glioblastoma multiforme (GBM) tissue from a surgically removed brain tumor of a patient. Significantly higher NANOGP8 transcription was observed in GBM cancer stem cells (CSCs) than in GBM cancer cells or neural stem cells (NSCs), despite identical sequences of NANOGP8-upstream genomic region in all the cell lines. This finding suggests that upstream genomic sequences of NANOGP8 may have environment-dependent promoter activity. We also found that the regulatory sequences upstream of exosomal NANOGP8 GBM DNA contain multiple core promoter elements, transcription factor binding sites, and segments of human viruses known for their oncogenic role. The exosomal sequence of NANOGP8-upstream GBM DNA is different from corresponding genomic sequences in CSCs, cancer cells, and NSCs as well as from the sequences reported by NCBI. These sequence dissimilarities suggest that exosomal NANOGP8 GBM DNA may not be a part of the genomic DNA. Exosomes possibly acquire this DNA from other sources where it is synthesized by an unknown mechanism. The significance of exosome-bestowed regulatory elements in the transcription of promoter-less retrogene such as NANOGP8 remains to be determined.

NestedMICA is a new, scalable, pattern-discovery system for finding transcription factor binding sites and similar motifs in biological sequences. Like several previous methods, NestedMICA tackles this problem by optimizing a probabilistic mixture model to fit a set of sequences. However, the use of a newly developed inference strategy called Nested Sampling means NestedMICA is able to find optimal solutions without the need for a problematic initialization or seeding step. We investigate the performance of NestedMICA in a range scenario, on synthetic data and a well-characterized set of muscle regulatory regions, and compare it with the popular MEME program. We show that the new method is significantly more sensitive than MEME: in one case, it successfully extracted a target motif from background sequence four times longer than could be handled by the existing program. It also performs robustly on synthetic sequences containing multiple significant motifs. When tested on a real set of regulatory sequences, NestedMICA produced motifs which were good predictors for all five abundant classes of annotated binding sites.

The emergence of SARS-CoV-2 variants is cause for concern, because these may become resistant to current vaccines and antiviral drugs in development. Current drugs target viral proteins, resulting in a critical need for RNA-targeted nanomedicines. To address this, a comparative analysis of SARS-CoV-2 variants was performed. Several highly conserved sites were identified, of which the most noteworthy is a partial homopurine palindrome site with >99% conservation within the coding region. This sequence was compared among recently emerged, highly infectious SARS-CoV-2 variants. Conservation of the site was maintained among these emerging variants, further contributing to its potential as a regulatory target site for SARS-CoV-2. RNAfold was used to predict the structures of the highly conserved sites, with some resulting structures being common among coronaviridae. An RNA-level regulatory map of the conserved regions of SARS-CoV-2 was produced based on the predicted structures, with each representing potential target sites for antisense oligonucleotides, triplex-forming oligomers, and aptamers. Additionally, homopurine/homopyrimidine sequences within the viral genome were identified. These sequences also demonstrate appropriate target sites for antisense oligonucleotides and triplex-forming oligonucleotides. An experimental strategy to investigate these is summarized along with potential nanoparticle types for delivery, and the advantages and disadvantages of each are discussed.

Clustered regularly interspaced short palindromic repeats (CRISPR)-associated 9 (Cas9) systems have been effectively harnessed to engineer the genomes of organisms from across the tree of life. Nearly all currently characterized Cas9 proteins are derived from mesophilic bacteria, and canonical Cas9 systems are challenged by applications requiring enhanced stability or elevated temperatures. We discovered IgnaviCas9, a Cas9 protein from a hyperthermophilic Ignavibacterium identified through mini-metagenomic sequencing of samples from a hot spring. IgnaviCas9 is active at temperatures up to 100 °C in vitro, which enables DNA cleavage beyond the 44 °C limit of Streptococcus pyogenes Cas9 (SpyCas9) and the 70 °C limit of both Geobacillus stearothermophilus Cas9 (GeoCas9) and Geobacillus thermodenitrificans T12 Cas9 (ThermoCas9). As a potential application of this enzyme, we demonstrate that IgnaviCas9 can be used in bacterial RNA-seq library preparation to remove unwanted cDNA from 16s ribosomal rRNA without increasing the number of steps, thus underscoring the benefits provided by its exceptional thermostability in improving molecular biology and genomic workflows. IgnaviCas9 is an exciting addition to the CRISPR-Cas9 toolbox and expands its temperature range.

Conserved nucleic acid sequences play an essential role in transcriptional regulation. The motifs/templates derived from nucleic acid sequence datasets are usually used as biomarkers to predict biochemical properties such as protein binding sites or to identify specific non-coding RNAs. In many cases, template-based nucleic acid sequence classification performs better than some feature extraction methods, such as N-gram and k-spaced pairs classification. The availability of large-scale experimental data provides an unprecedented opportunity to improve motif extraction methods. The process for pattern extraction from large-scale data is crucial for the creation of predictive models.

Computational approaches for finding DNA regulatory motifs in promoter sequences are useful to biologists in terms of reducing the experimental costs and speeding up the discovery process of de novo binding sites. It is important for rule-based or clustering-based motif searching schemes to effectively and efficiently evaluate the similarity between a k-mer (a k-length subsequence) and a motif model, without assuming the independence of nucleotides in motif models or without employing computationally expensive Markov chain models to estimate the background probabilities of k-mers. Also, it is interesting and beneficial to use a priori knowledge in developing advanced searching tools.

Eukaryotic nucleic acid binding protein database (ENPD, http://qinlab.sls.cuhk.edu.hk/ENPD/) is a library of nucleic acid binding proteins (NBPs) and their functional information. NBPs such as DNA binding proteins (DBPs), RNA binding proteins (RBPs), and DNA and RNA binding proteins (DRBPs) are involved in every stage of gene regulation through their interactions with DNA and RNA. Due to the importance of NBPs, the database was constructed based on manual curation and a newly developed pipeline utilizing both sequenced transcriptomes and genomes. In total the database has recorded 2.8 million of NBPs and their binding motifs from 662 NBP families and 2423 species, constituting the largest NBP database. ENPD covers evolutionarily important lineages which have never been included in the previous NBP databases, while lineage-specific NBP family expansions were also found. ENPD also focuses on the involvements of DBPs, RBPs and DRBPs in non-coding RNA (ncRNA) mediated gene regulation. The predicted and experimentally validated targets of NBPs have both been recorded and manually curated in ENPD, linking the interactions between ncRNAs, DNA regulatory elements and NBPs in gene regulation. This database provides key resources for the scientific community, laying a solid foundation for future gene regulatory studies from both functional and evolutionary perspectives.

The male reproductive system has a standard immune response regulatory mechanism, However, a variety of external stimuli, including viruses, bacteria, heat, and medications can damage the testicles and cause orchitis and epididymitis. It has been shown that various RNA viruses are more likely to infect the testis than DNA viruses, inducing orchitis and impairing testicular function. It was found that local injection of the viral RNA analog poly(I:C) into the testes markedly disrupted the structure of the seminiferous tubules, accompanied by apoptosis and inflammation. Poly(I:C) mainly inhibited the expression of testosterone synthesis-associated proteins, STAR and MGARP, and affected the synthesis and metabolism of amino acids and lipids in the testis. This led to the disruption of the metabolite levels in the testis of mice, thus affecting the normal spermatogenesis process. The present study analyzed the acute inflammatory response of the testis to viral infection using a multi-omics approach. It provides insights into how RNA virus infection impairs testicular function and offers a theoretical basis for future studies on immune homeostasis and responses under stress conditions in male reproduction.

Limited data are available on the incidence of variations in nucleotide sequences of long terminal repeat (LTR) regions of Bovine Leukemia Virus (BLV). Consequently, the possible impact of SNPs on BLV LTR function are poorly elucidated. Thus, a detailed and representative study of full-length LTR sequences obtained from sixty-four BLV isolates from different geographical regions of Poland, Moldova, Croatia, Ukraine and Russia were analyzed for their genetic variability.

Regulatory sequences are not solely defined by their nucleic acid sequence but also by their relative distances to genomic landmarks such as transcription start site, exon boundaries or polyadenylation site. Deep learning has become the approach of choice for modeling regulatory sequences because of its strength to learn complex sequence features. However, modeling relative distances to genomic landmarks in deep neural networks has not been addressed.

Nucleic acid-based adjuvants such as CpG oligonucleotides (CpG ODNs) and poly(I:C) are potential vaccine adjuvants for infectious diseases and cancers. However, the mechanism by which their cell surface receptors promote their uptake into dendritic cells (DCs) and shuttle them to intracellular Toll-like receptors remains to be further investigated. Here, we demonstrated a role for nucleolin, a multifunctional DNA- and RNA-binding protein and a major constituent of the nucleolus, as one of the cell-surface receptors for nucleic acid-based adjuvants. Nucleolin on mouse DC surface bound directly to A-type CpG ODN, B-type CpG ODN, and poly(I:C) and promoted their internalization into cells following DC maturation in vitro. In human DCs, nucleolin also contributed to the binding and internalization of both types of CpG ODNs and subsequent cytokine production. Furthermore, nucleolin played a crucial role in cytokine production and activating antigen-specific antibodies and T cell responses induced by B-type CpG ODN in vivo in mice. Our findings provide valuable information that can help improve the efficacy and safety of these adjuvants.

Hedgehog (Hh) genes play major roles in animal development and studies of their evolution, expression and function point to major differences among chordates. Here we focused on Hh genes in lampreys in order to characterize the evolution of Hh signalling at the emergence of vertebrates. Screening of a cosmid library of the river lamprey Lampetra fluviatilis and searching the preliminary genome assembly of the sea lamprey Petromyzon marinus indicate that lampreys have two Hh genes, named Hha and Hhb. Phylogenetic analyses suggest that Hha and Hhb are lamprey-specific paralogs closely related to Sonic/Indian Hh genes. Expression analysis indicates that Hha and Hhb are expressed in a Sonic Hh-like pattern. The two transcripts are expressed in largely overlapping but not identical domains in the lamprey embryonic brain, including a newly-described expression domain in the nasohypophyseal placode. Global alignments of genomic sequences and local alignment with known gnathostome regulatory motifs show that lamprey Hhs share conserved non-coding elements (CNE) with gnathostome Hhs albeit with sequences that have significantly diverged and dispersed. Functional assays using zebrafish embryos demonstrate gnathostome-like midline enhancer activity for CNEs contained in intron2. We conclude that lamprey Hh genes are gnathostome Shh-like in terms of expression and regulation. In addition, they show some lamprey-specific features, including duplication and structural (but not functional) changes in the intronic/regulatory sequences.

The JC virus (JCV) control region contains AGGGAAGGGA, the tandem pentanucleotide repeat element (Pnt2). Several proteins specifically interacted via Pnt2 to regulate the expression of JCV early promoter-enhancer (JCV(E)) or late promoter-enhancer (JCV(L)). In this study, a JCV Pnt2 oligonucleotide probe was used to screen a cDNA expression library from glial P19 mouse embryonal carcinoma cells. A cDNA clone was isolated by Southwestern blot assay and it produced a protein that reproducibly and specifically bound to Pnt2. This cDNA had 100% homology to one of three previously identified mouse cDNAs called cellular nucleic acid binding proteins (Cnbps). Cnbps are a highly homologous family of eukaryotic genes implicated in functional interactions with cytoplasmic RNA and regulatory DNA elements. An mRNA of 2.2 kb of Pnt2-interacting Cnbp (PCnbp) was seen in undifferentiated, muscle or glial P19 cells. When expressed from a cDNA expression vector as a fusion protein that also contained 115 kDa from beta-galactosidase, a Pnt2 binding protein (PCNBP) specifically bound to Pnt2 in Southwestern blots as a 30 kDa component of the 145 kDa fusion protein. Furthermore, JCV(E) expression was negatively regulated by PCnbp produced in vivo from the cDNA expression vector. Regulation of JCV(L) was unaffected. We suggest a novel role for CNBP as a PCNBP that interacts with Pnt2 in the negative transcriptional regulation of JCV(E).

Alternative splicing of the bcl-x gene generates two transcripts: the anti-apoptotic bcl-xL isoform and the pro-apoptotic bcl-xS isoform. The ratio between the two isoforms is a key factor in development and in cancer progression. Here, we show that a short antisense chimeric peptide nucleic acid (PNA) oligonucleotide conjugated to a polypeptide containing eight Ser-Arg repeats (SR)(8) can modulate splicing of bcl-x both in vitro and in vivo and induces apoptosis in HeLa cells. The PNA-SR oligo was targeted to a region of bcl-x that does not contain splicing regulatory sequences and was able to override the complex network of splicing enhancers and silencers that regulates the ratio between the two bcl-x isoforms. Thus, PNA-SR oligos are powerful tools that can potentially modulate splice site choice in endogenous genes independent of the presence of other splicing regulatory mechanisms on the target gene.

NASA has made great strides in the past five years to develop a suite of instruments for the International Space Station in order to perform molecular biology in space. However, a key piece of equipment that has been lacking is an instrument that can extract nucleic acids from an array of complex human and environmental samples. The Omics in Space team has developed the μTitan (simulated micro(μ) gravity tested instrument for automated nucleic acid) system capable of automated, streamlined, nucleic acid extraction that is adapted for use under microgravity. The μTitan system was validated using a whole cell microbial reference (WCMR) standard comprised of a suspension of nine bacterial strains, titrated to concentrations that would challenge the performance of the instrument, as well as to determine the detection limits for isolating DNA. Quantitative assessment of system performance was measured by comparing instrument input challenge dose vs recovery by Qubit spectrofluorometry, qPCR, Bioanalyzer, and Next Generation Sequencing. Overall, results indicate that the μTitan system performs equal to or greater than a similar commercially available, earth-based, automated nucleic acid extraction device. The μTitan system was also tested in Yellowstone National Park (YNP) with the WCMR, to mimic a remote setting, with limited resources. The performance of the device at YNP was comparable to that in a laboratory setting. Such a portable, field-deployable, nucleic extraction system will be valuable for environmental microbiology, as well as in health care diagnostics.

The interplay of different virus species in a host cell after infection can affect the adaptation of each virus. Endogenous viral elements, such as endogenous pararetroviruses (PRVs), have arisen from vertical inheritance of viral sequences integrated into host germline genomes. As viral genomic fossils, these sequences can thus serve as valuable paleogenomic data to study the long-term evolutionary dynamics of virus-virus interactions, but they have rarely been applied for this purpose. All extant PRVs have been considered autonomous species in their parasitic life cycle in host cells. Here, we provide evidence for multiple non-autonomous PRV species with structural defects in viral activity that have frequently infected ancient grass hosts and adapted through interplay between viruses. Our paleogenomic analyses using endogenous PRVs in grass genomes revealed that these non-autonomous PRV species have participated in interplay with autonomous PRVs in a possible commensal partnership, or, alternatively, with one another in a possible mutualistic partnership. These partnerships, which have been established by the sharing of noncoding regulatory sequences (NRSs) in intergenic regions between two partner viruses, have been further maintained and altered by the sequence homogenization of NRSs between partners. Strikingly, we found that frequent region-specific recombination, rather than mutation selection, is the main causative mechanism of NRS homogenization. Our results, obtained from ancient DNA records of viruses, suggest that adaptation of PRVs has occurred by concerted evolution of NRSs between different virus species in the same host. Our findings further imply that evaluation of within-host NRS interactions within and between populations of viral pathogens may be important.

The TWIST1 bHLH transcription factor controls embryonic development and cancer processes. Although molecular and genetic analyses have provided a wealth of data on the role of bHLH transcription factors, very little is known on the molecular mechanisms underlying their binding affinity to the E-box sequence of the promoter. Here, we used an in silico model of the TWIST1/E12 (TE) heterocomplex and performed molecular dynamics (MD) simulations of its binding to specific (TE-box) and modified E-box sequences. We focused on (i) active E-box and inactive E-box sequences, on (ii) modified active E-box sequences, as well as on (iii) two box sequences with modified adjacent bases the AT- and TA-boxes. Our in silico models were supported by functional in vitro binding assays. This exploration highlighted the predominant role of protein side-chain residues, close to the heart of the complex, at anchoring the dimer to DNA sequences, and unveiled a shift towards adjacent ((-1) and (-1*)) bases and conserved bases of modified E-box sequences. In conclusion, our study provides proof of the predictive value of these MD simulations, which may contribute to the characterization of specific inhibitors by docking approaches, and their use in pharmacological therapies by blocking the tumoral TWIST1/E12 function in cancers.

Hydroxyproline-rich glycoproteins (HRGPs) are a superfamily of plant cell wall structural proteins that function in various aspects of plant growth and development, including pollen tube growth. We have previously characterized protein sequence signatures for three family members in the HRGP superfamily: the hyperglycosylated arabinogalactan-proteins (AGPs), the moderately glycosylated extensins (EXTs), and the lightly glycosylated proline-rich proteins (PRPs). However, the mechanism of pollen-specific HRGP gene expression remains unexplored. To this end, we developed an integrative analysis pipeline combining RNA-seq gene expression and promoter sequences to identify cis-regulatory motifs responsible for pollen-specific expression of HRGP genes in Arabidopsis thaliana. Specifically, we mined the public RNA-seq datasets and identified 13 pollen-specific HRGP genes. Ensemble motif discovery identified 15 conserved promoter elements between A.thaliana and A. lyrata. Motif scanning revealed two pollen related transcription factors: GATA12 and brassinosteroid (BR) signaling pathway regulator BZR1. Finally, we performed a regression analysis and demonstrated that the 15 motifs provided a good model of HRGP gene expression in pollen (R = 0.61). In conclusion, we performed the first integrative analysis of cis-regulatory motifs in pollen-specific HRGP genes, revealing important insights into transcriptional regulation in pollen tissue.

Chronic hepatitis B infection (CHB) is an area of high unmet medical need. Current standard-of-care therapies only rarely lead to a functional cure, defined as durable hepatitis B surface antigen (HBsAg) loss following treatment. The goal for next generation CHB therapies is to achieve a higher rate of functional cure with finite treatment duration. To address this urgent need, we are developing liver-targeted single-stranded oligonucleotide (SSO) therapeutics for CHB based on the locked nucleic acid (LNA) platform. These LNA-SSOs target hepatitis B virus (HBV) transcripts for RNase-H-mediated degradation. Here, we describe a HBV-specific LNA-SSO that effectively reduces intracellular viral mRNAs and viral antigens (HBsAg and HBeAg) over an extended time period in cultured human hepatoma cell lines that were infected with HBV with mean 50% effective concentration (EC50) values ranging from 1.19 to 1.66 μM. To achieve liver-specific targeting and minimize kidney exposure, this LNA-SSO was conjugated to a cluster of three N-acetylgalactosamine (GalNAc) moieties that direct specific binding to the asialoglycoprotein receptor (ASGPR) expressed specifically on the surface of hepatocytes. The GalNAc-conjugated LNA-SSO showed a strikingly higher level of potency when tested in the AAV-HBV mouse model as compared with its non-conjugated counterpart. Remarkably, higher doses of GalNAc-conjugated LNA-SSO resulted in a rapid and long-lasting reduction of HBsAg to below the detection limit for quantification, i.e., by 3 log10 (p < 0.0003). This antiviral effect depended on a close match between the sequences of the LNA-SSO and its HBV target, indicating that the antiviral effect is not due to non-specific oligonucleotide-driven immune activation. These data support the development of LNA-SSO therapeutics for the treatment of CHB infection.