Missing Nucleotides Research Articles

Изоформы микроРНК miR-148a и miR-203a предположительно играют роль супрессоров колоректального рака

MicroRNAs are short non-coding molecules which regulate translation in a gene-specific manner. MicroRNA isoforms that differ by few extra or missing nucleotides at the 5'-terminus (5'-isomiR) show strikingly different target specificity. This study aimed to identify functional roles of 5′-isomiR in colorectal cancers. Transcriptomic targets of microRNA isoforms were predicted using bioinformatics tools miRDB and TargetScan. The sets of putative targets identified for 5′-isomiR were integrated with mRNA and microRNA sequencing data for primary colorectal tumors retrieved from The Cancer Genome Atlas Colon Adenocarcinoma (TCGA-COAD) database. The network of interactions among miRNA, their targets and transcription factors was built using the miRGTF-net algorithm. The results indicate that microRNA isoforms highly expressed in colorectal cancer and differing by a single nucleotide position at the 5'-terminus have ≤ 30% common targets. The regulatory network of interactions enables identification of the most engaged microRNA isoforms. Anti-correlated expression levels of canonical microRNA hsa-miR-148a-3p and its putative targets including CSF1, ETS1, FLT1, ITGA5, MEIS1, MITF and RUNX2 proliferation regulators suggest an anti-tumor role for this molecule. The canonical microRNA hsa-miR-203a-3p|0 and its 5′-isoform bind different sets of anti-correlated putative targets, although both of them interact with genes involved in the epithelial-mesenchymal transition: SNAI2 and TNC.

Two linked resistance genes function divergently in defence against Verticillium Wilt in Alfalfa.

Two linked resistance genes function divergently in defence against Verticillium Wilt in Alfalfa.

A globally diverse reference alignment and panel for imputation of mitochondrial DNA variants

BackgroundVariation in mitochondrial DNA (mtDNA) identified by genotyping microarrays or by sequencing only the hypervariable regions of the genome may be insufficient to reliably assign mitochondrial genomes to phylogenetic lineages or haplogroups. This lack of resolution can limit functional and clinical interpretation of a substantial body of existing mtDNA data. To address this limitation, we developed and evaluated a large, curated reference alignment of complete mtDNA sequences as part of a pipeline for imputing missing mtDNA single nucleotide variants (mtSNVs). We call our reference alignment and pipeline MitoImpute.ResultsWe aligned the sequences of 36,960 complete human mitochondrial genomes downloaded from GenBank, filtered and controlled for quality. These sequences were reformatted for use in imputation software, IMPUTE2. We assessed the imputation accuracy of MitoImpute by measuring haplogroup and genotype concordance in data from the 1000 Genomes Project and the Alzheimer’s Disease Neuroimaging Initiative (ADNI). The mean improvement of haplogroup assignment in the 1000 Genomes samples was 42.7% (Matthew’s correlation coefficient = 0.64). In the ADNI cohort, we imputed missing single nucleotide variants.ConclusionThese results show that our reference alignment and panel can be used to impute missing mtSNVs in existing data obtained from using microarrays, thereby broadening the scope of functional and clinical investigation of mtDNA. This improvement may be particularly useful in studies where participants have been recruited over time and mtDNA data obtained using different methods, enabling better integration of early data collected using less accurate methods with more recent sequence data.

Enhancing PCR Capacity To Detect 'Candidatus Liberibacter asiaticus' Utilizing Whole Genome Sequence Information.

'Candidatus Liberibacter asiaticus' (CLas) is an unculturable α-proteobacterium associated with citrus Huanglongbing (HLB; yellow shoot disease). PCR procedures that accurately confirm or exclude CLas infection in citrus tissue/Asian citrus psyllid (ACP) samples are critical for HLB management. When CLas was described in 1994, a 23-bp signature oligonucleotide sequence (OI1) in the 16S rRNA gene (rrs, three genomic copies) was identified based on Sanger sequencing. OI1 contains single nucleotide polymorphisms (SNPs) distinguishing CLas from non-CLas species. The SNPs were used to design the primer HLBas, a key primer for a commonly used TaqMan PCR system (HLBas-PCR) for CLas detection. Recent developments in next-generation sequencing technology have led to the identification of 15 CLas whole genome sequence strains (WGSs). Analyses of CLas WGSs have generated a significant amount of biological information that could help to improve CLas detection. Utilizing the WGS information, this study re-evaluated the sequence integrity of OI1/HLBas and identified and/or confirmed a missing nucleotide G in the two primers. Replacement primers for OI1 and HLBas are proposed. At low cycle threshold (Ct) values (e.g., <30), HLBas-PCR remained reliable in CLas determination. At high Ct values (e.g., >30), HLBas-PCR alone was unreliable in differentiating whether samples contain low CLas titers or whether CLas is not present. The availability of ribonucleotide reductase (RNR)-PCR derived from the five-copy nrdB gene helped to resolve this problem. To further enhance low CLas titer detection, a 4CP-PCR system, based on a four-copy genomic locus, was developed. Evaluation of 107 HLB samples (94 citrus and 13 ACP) showed that 4CP-PCR was more sensitive than HLBas-PCR and shared similar sensitivity with RNR-PCR.

MitoIMP: A Computational Framework for Imputation of Missing Data in Low-Coverage Human Mitochondrial Genome.

The incompleteness of partial human mitochondrial genome sequences makes it difficult to perform relevant comparisons among multiple resources. To deal with this issue, we propose a computational framework for deducing missing nucleotides in the human mitochondrial genome. We applied it to worldwide mitochondrial haplogroup lineages and assessed its performance. Our approach can deduce the missing nucleotides with a precision of 0.99 or higher in most human mitochondrial DNA lineages. Furthermore, although low-coverage mitochondrial genome sequences often lead to a blurred relationship in the multidimensional scaling analysis, our approach can correct this positional arrangement according to the corresponding mitochondrial DNA lineages. Therefore, our framework will provide a practical solution to compensate for the lack of genome coverage in partial and fragmented human mitochondrial genome sequences. In this study, we developed an open-source computer program, MitoIMP, implementing our imputation procedure. MitoIMP is freely available from https://github.com/omics-tools/mitoimp.

Trinucleotide Rolling Circle Amplification: A Novel Method for the Detection of RNA and DNA

Most natural DNA and RNA are devoid of long trinucleotide (TN) sequences that lack one specific nucleotide (missing nucleotide (MN)). Here we developed a novel method that is based on rolling circle amplification (RCA), in which the TN-information of short TN stretches is sequence-specifically recognized, transferred, extended, amplified and detected by padlock probes that consist entirely of nucleotides complementary to the three nucleotides present in the target sequence (complementary TN-information). Upon specific head-to-tail annealing and ligation to the TN-target sequence, these padlock probes represent extended complementary TN versions of the target sequence that can be further amplified by trinucleotide rolling circle amplification (TN-RCA). Since during TN-RCA the MN (as dNTP) is not added, background amplification is minimized with endogenous RNA/DNA (which mostly would require all four dNTP). Therefore, various labelled dNTP can be added to the TN-RCA reaction that enables the separation, isolation and detection of the amplified single-stranded DNA (ssDNA). Here the TN-RCA method is exemplified with RNA/DNA from Zika virus and from human papilloma virus (HPV). TN-RCA is a novel isothermal amplification technique that can be used for sensitive sequence-specific detection and diagnosis of natural and synthetic DNA or RNA containing TN stretches with low background in short time.

Imputation approach for deducing a complete mitogenome sequence from low-depth-coverage next-generation sequencing data: application to ancient remains from the Moon Pyramid, Mexico.

It is considered that more than 15 depths of coverage are necessary for next-generation sequencing (NGS) data to obtain reliable complete nucleotide sequences of the mitogenome. However, it is difficult to satisfy this requirement for all nucleotide positions because of problems obtaining a uniform depth of coverage for poorly preserved materials. Thus, we propose an imputation approach that allows a complete mitogenome sequence to be deduced from low-depth-coverage NGS data. We used different types of mitogenome data files as panels for imputation: a worldwide panel comprising all the major haplogroups, a worldwide panel comprising sequences belonging to the estimated haplogroup alone, a panel comprising sequences from the population most closely related to an individual under investigation, and a panel comprising sequences belonging to the estimated haplogroup from the population most closely related to an individual under investigation. The number of missing nucleotides was drastically reduced in all the panels, but the contents obtained by imputation were quite different among the panels. The efficiency of the imputation method differed according to the panels used. The missing nucleotides were most credibly imputed using sequences of the estimated haplogroup from the population most closely related to the individual under investigation as a panel.

TSGΔ154-1054 splice variant increases TSG101 oncogenicity by inhibiting its E3-ligase-mediated proteasomal degradation.

Tumor susceptibility gene 101 (TSG101) elicits an array of cellular functions, including promoting cytokinesis, cell cycle progression and proliferation, as well as facilitating endosomal trafficking and viral budding. TSG101 protein is highly and aberrantly expressed in various human cancers. Specifically, a TSG101 splicing variant missing nucleotides 154 to 1054 (TSGΔ154-1054), which is linked to progressive tumor-stage and metastasis, has puzzled investigators for more than a decade. TSG101-associated E3 ligase (Tal)- and MDM2-mediated proteasomal degradation are the two major routes for posttranslational regulation of the total amount of TSG101. We reveal that overabundance of TSG101 results from TSGΔ154-1054 stabilizing the TSG101 protein by competitively binding to Tal, but not MDM2, thereby perturbing the Tal interaction with TSG101 and impeding subsequent polyubiquitination and proteasomal degradation of TSG101. TSGΔ154-1054 therefore specifically enhances TSG101-stimulated cell proliferation, clonogenicity, and tumor growth in nude mice. This finding shows the functional significance of TSGΔ154-1054 in preventing the ubiquitin-proteasome proteolysis of TSG101, which increases tumor malignancy and hints at its potential as a therapeutic target in cancer treatment.

Lentiviral vector platform for improved erythropoietin expression concomitant with shRNA mediated host cell elastase down regulation

Lentiviral vector (LV) mediated gene transfer holds great promise to develop stable cell lines for sustained transgene expression providing a valuable alternative to the conventional plasmid transfection based recombinant protein production methods. We report here making a third generation HIV-2 derived LV containing erythropoietin (EPO) gene expression cassette to generate a stable HEK293 cell line secreting EPO constitutively. A high producer cell clone was obtained by limiting dilution and was adapted to serum free medium. The suspension adapted cell clone stably produced milligram per liter quantities of EPO. Subsequent host metabolic engineering using lentiviral RNAi targeted to block an endogenous candidate protease elastase, identified through an in silico approach, resulted in appreciable augmentation of EPO expression above the original level. This study of LV based improved glycoprotein expression with host cell metabolic engineering for stable production of protein therapeutics thus exemplifies the versatility of LV and is of significant future biopharmaceutical importance.

Bunyamwera virus can repair both insertions and deletions during RNA replication.

The genomic termini of RNA viruses contain essential cis-acting signals for such diverse functions as packaging, genome translation, mRNA transcription, and RNA replication, and thus preservation of their sequence integrity is critical for virus viability. Sequence alteration can arise due to cellular mechanisms that add or remove nucleotides from terminal regions, or, alternatively, from introduction of sequence errors through nucleotide misincorporation by the error-prone viral RNA-dependent RNA polymerase (RdRp). To preserve template function, many RNA viruses utilize repair mechanisms to prevent accumulation of terminal alterations. Here we show that Bunyamwera virus (BUNV), the prototype of the Bunyaviridae family of segmented negative-sense RNA viruses, also can repair its genomic termini. When an intact nontranslated region (NTR) was added to the anti-genomic 3' end, it was precisely removed, to restore both length and RNA synthesis function of the wild-type template. Furthermore, when nucleotides were removed from the anti-genome 3' end, and replaced with a duplicate and intact NTR, both the external NTR were removed, and the missing nucleotides were restored, thus, indicating that the BUNV RdRp can both remove and add nucleotides to the template. We show that the mechanism for repair of terminal extensions is likely that of internal entry of the viral RdRp during genome synthesis. Possible mechanisms for repair of terminal deletions are discussed.

Long Patch Base Excision Repair Proceeds via Coordinated Stimulation of the Multienzyme DNA Repair Complex

Base excision repair, a major repair pathway in mammalian cells, is responsible for correcting DNA base damage and maintaining genomic integrity. Recent reports show that the Rad9-Rad1-Hus1 complex (9-1-1) stimulates enzymes proposed to perform a long patch-base excision repair sub-pathway (LP-BER), including DNA glycosylases, apurinic/apyrimidinic endonuclease 1 (APE1), DNA polymerase beta (pol beta), flap endonuclease 1 (FEN1), and DNA ligase I (LigI). However, 9-1-1 was found to produce minimal stimulation of FEN1 and LigI in the context of a complete reconstitution of LP-BER. We show here that pol beta is a robust stimulator of FEN1 and a moderate stimulator of LigI. Apparently, there is a maximum possible stimulation of these two proteins such that after responding to pol beta or another protein in the repair complex, only a small additional response to 9-1-1 is allowed. The 9-1-1 sliding clamp structure must serve primarily to coordinate enzyme actions rather than enhancing rate. Significantly, stimulation by the polymerase involves interaction of primer terminus-bound pol beta with FEN1 and LigI. This observation provides compelling evidence that the proposed LP-BER pathway is actually employed in cells. Moreover, this pathway has been proposed to function by sequential enzyme actions in a "hit and run" mechanism. Our results imply that this mechanism is still carried out, but in the context of a multienzyme complex that remains structurally intact during the repair process.

The Ryegrass mottle virus genome codes for a sobemovirus 3C-like serine protease and RNA-dependent RNA polymerase translated via −1 ribosomal frameshifting

In the course of sobemovirus gene cloning the complete genome of Ryegrass mottle virus (RGMoV) was sequenced. Sequence analysis revealed differences including missing and extraneous nucleotides in comparison to the previously published sequence (Zhang, Toriyama, Takanashi, J. Gen. Plant Pathol. 67, 63 (2001)). A gene coding for a typical sobemovirus 3C-like serine protease was identified in ORF2a after multiple sequence alignment analysis. The newly identified 57-amino-acid stretch in ORF2a showed similarities ranging from 38.5 to 50.9% among sequenced genes of sobemovirus proteases. ORF analysis of the RGMoV polyprotein coding sequence demonstrated the arrangement of ORF2b coding for RNA-dependent RNA polymerase (RdRP) in the -1 frame in regard to ORF2a. The localization of conserved among sobemoviruses slippery sequence (UUUAAAC) at the 3'-end of ORF2a suggests the translation of RdRP via a -1 ribosomal frameshifting mechanism, allowing to include the RGMoV in the sobemovirus group with a Cocksfoot mottle virus-like (CfMV-like) genome organization.

Determination of the 5′- and 3′-Terminal Sequences Completes the Sequences of the Two Double-stranded RNAs of Penicillium stoloniferum Virus S

The two genomic segments of Penicillium Stoloniferum virus S (PsV-S), a member of the Partitiviridae, were recently sequenced and published. We independantly sequenced PsV-S and showed that the original sequence was missing nucleotides at both the 5' and 3' termini of both segments. We determined the correct sequence in three independent experiments and found the segments to be 1753 bp (encoding the RNA-dependant RNA polymerase) and 1581 bp (encoding the Capsid Protein). Homology was shown between the 5' and 3' ends of PsV-S and other members of the Partitiviridae.

A New Effective Method for Estimating Missing Values in the Sequence Data Prior to Phylogenetic Analysis

In this article we address the problem of phylogenetic inference from nucleic acid data containing missing bases. We introduce a new effective approach, called “Probabilistic estimation of missing values” (PEMV), allowing one to estimate unknown nucleotides prior to computing the evolutionary distances between them. We show that the new method improves the accuracy of phylogenetic inference compared to the existing methods “Ignoring Missing Sites” (IMS), “Proportional Distribution of Missing and Ambiguous Bases” (PDMAB) included in the PAUP software [ 26 ]. The proposed strategy for estimating missing nucleotides is based on probabilistic formulae developed in the framework of the Jukes-Cantor [ 10 ] and Kimura 2-parameter [ 11 ] models. The relative performances of the new method were assessed through simulations carried out with the SeqGen program [ 20 ], for data generation, and the BioNJ method [ 7 ], for inferring phylogenies. We also compared the new method to the DNAML program [ 5 ] and “Matrix Representation using Parsimony” (MRP) [ 13 , 19 ] considering an example of 66 eutherian mammals originally analyzed in [ 17 ].

Efficient Rescue of Infectious Bursal Disease Virus from Cloned cDNA: Evidence for Involvement of the 3′-Terminal Sequence in Genome Replication

To study the mechanism of replication of infectious bursal disease virus (IBDV), and to determine factors on the IBDV RNA which are involved in viral replication, we used cloned full-length cDNA of both the A- and B-segments to generate infectious IBDV. Infectious IBDV was rescued from plasmids that contained full-length IBDV cDNA behind a T7 promoter, by transfecting these plasmids into cells which were infected with a recombinant Fowlpox virus that expressed T7 RNA polymerase. By using the cDNA transfection system we evaluated the effect of the length of the 3′ terminus of the A-segment plus strand of IBDV. Although wild-type IBDV predominantly contains four cytosines at the 3′ terminus, no difference in virus yield was found when virus was rescued from cDNAs containing three to six adjacent cytosines. When the 3′ terminus was shorter than three cytosines the efficiency to generate infectious IBDV from cDNA was reduced, but IBDV could still be recovered reproducibly. The rescued viruses from cDNAs containing 3′-terminal deletions appeared to have a restored 3′-terminal sequence. The missing nucleotides are probably restored by using complementary bases of a stem-loop structure as template.

Negative electrospray ionization mass spectrometry of synthetic and chemically modified oligonucleotides.

We report here on the analysis of synthetic oligonucleotides by electrospray ionization mass spectrometry (ESI-MS). After intensive removal of salt ions (especially sodium cations), negative ion mass spectra, allowing mass measurement with an accuracy of 0.01%, were obtained on several oligonucleotides up to 80 nucleotides. In most cases, the resolution was sufficient to observe n-1 and n-2 forms due to internal deletions during automated synthesis, and to identify the missing nucleotides. A 132-mer, whose size is close to the limit of automated chemical synthesis, was also successfully mass measured. A quantitative study showed that ESI-MS can provide quantitative data on oligonucleotides of similar size and structure. The described methodology is used to characterize oligonucleotide analogues such as phosphorothioate oligonucleotides designed for antisense applications. Finally, analyses in the positive ion mode on a trimer TpTpT in the presence of different amine bases were performed and allowed a better understanding of the influence of these bases on the ions formation.

Characterization of two isozymic forms of heart fructose 6-phosphate, 2-kinase:fructose 2,6-bisphosphatase

Occurrence of two isozymic forms of fructose 6-P, 2-kinase: fructose 2,6-bisphosphatase in bovine heart was investigated by transcribing mRNAs and amplifying the cDNAs with polymerase chain reactions. Analysis of the PCR products revealed 1.7 Kb and 1.5 Kb DNAs, and the determination of their nucleotide sequences showed that these DNAs are identical except for the lack of 180 base pairs near the 3' of the bovine heart enzyme DNA previously reported (6). This missing nucleotide sequence encodes Asn451-Gln510 and contains the phosphorylation sites for cAMP dependent protein kinase and protein kinase C.

Deletion of repeated sequences from tobacco mosaic virus mutants with two coat protein genes

Tobacco mosaic virus (TMV) hybrids were constructed in vitro with an extra coat protein gene (gene order 5′- 126K/183K gene-30K gene-first coat protein gene-subgenomic mRNA promoter:second coat protein gene-3′). One coat protein gene was wild type and the other was derived from one of three coat protein mutants that elicit distinctive phenotypic responses. The hybrids contained large sequence repeats (up to 747 nucleotides). RNA transcripts of each of seven hybrid constructs initiated infections in host plants, but the progeny virus contained only one coat protein gene with no additional or missing nucleotides. Progeny of mutants with the two functional coat protein genes in different order within the genome retained coat protein gene sequence predicted by the number of repeated nucleotides 3′ and 5′ of the differential nucleotides. In contrast, progeny that resulted from virus hybrids with a wild-type and a nonfunctional coat protein gene possessed only the wild-type coat protein gene, even when the predicted progeny was a free-RNA virus. This demonstrated that selection strongly affected which progeny became predominant in the resulting populations

Nucleotide and deduced amino acid sequence of the alpha subunit of yeast pyruvate dehydrogenase

A 413-base cDNA insert encoding a portion of the alpha subunit of pyruvate dehydrogenase (E 1α; EC 1.2.4.1) from Saccharomyces cerevisiae was isolated from a λgtll cDNA library by immunoscreening and by hybridization with an oligonucleotide probe which corresponded to the amino acid sequence around the phosphorylation site of E 1α. This cDNA was subcloned, sequenced and used as a probe to isolate two additional cDNA inserts which were subcloned and sequenced. These overlapping clones comprised the carboxyl-terminal part of E 1α. To identify the missing nucleotide sequence, the polymerase chain reaction was used to amplify yeast genomic DNA with synthetic oligonucleotide primers based on the amino-terminal sequence of E 1α and the 5′ end of one of the cDNA clones. Three DNA fragments were isolated and sequenced. The composite nucleotide sequence has an open reading frame of 1260 nucleotides encoding a putative presequence of 33 amino acids and a mature protein of 387 amino acids ( M r = 42,703). Hybridization analysis showed that the size of the mRNA is about I.4 kilobases.

Nicks 3' or 5' to AP sites or to mispaired bases, and one-nucleotide gaps can be sealed by T4 DNA ligase.

Using synthetic oligodeoxynucleotides with 3'-OH ends and 32P-labelled 5'-phosphate ends and the technique of polyacrylamide gel electrophoresis, it is shown that, in the presence of the complementary polynucleotide, an AP (apurinic or apyrimidinic) site at the 3' or the 5' end of the labelled oligodeoxynucleotides does not prevent their ligation by T4 DNA ligase, although the reaction rate is decreased. This decrease is more severe when the AP site is at the 3' end; the activated intermediates accumulate showing that it is the efficiency of the adenyl-5'-phosphate attack by the 3'-OH of the base-free deoxyribose which is mostly perturbed. Using the same technique, it is shown that a mispaired base at the 3' or 5' end of oligodeoxynucleotides does not prevent their ligation. A one-nucleotide gap, limited by 3'-OH and 5'-phosphate, can also be closed by T4 DNA ligase although with difficulty; here again the activation of the 5'-phosphate end does not seem to be slowed down, but rather the 3'-OH attack of the adenyl-5'-phosphate. All these anomalous ligations take place with the nick or the gap in front of a continuous complementary strand. Blunt ends ligation of correct duplexes occurs readily; however an AP site or a mispaired base at the 3' or 5' end of one strand of the duplexes prevents ligation between these strands. But a missing nucleotide (responsible for one unpaired nucleotide protruding at the 3' or 5' end of the complementary strand) does not stop ligation of the shorter oligodeoxynucleotides between independent duplexes.