Analysis of mRNA transcripts improves the success rate of molecular genetic testing in OTC deficiency

BACKGROUND: Ornithine transcarbamylase (OTC) deficiency is the most common inborn error of urea metabolism that can lead to hyperammonemic crises and orotic aciduria. To date, a total of 341 causative mutations within the OTC gene have been described. However, in about 20% of the patients with enzymatically confirmed OTC deficiency no mutation can be detected when sequencing of genomic DNA analyzing exons and adjacent intronic segments of the OTC gene is performed. METHODS: Standard genomic DNA analysis of the OTC gene in five consecutive patients from five families revealed no mutation. Hence, liver tissue was obtained by needle sampling or open biopsy and RNA extracted from liver was analyzed. RESULTS: Complex rearrangements of the OTC transcript (three insertions and two deletions) were found in all five patients. CONCLUSION: In patients with a strong suspicion of OTC deficiency despite normal results of sequencing exonic regions of the OTC gene, characterization of liver OTC mRNA is highly effective in resolving the genotype. Liver tissue sampling by needle aspiration allows for both enzymatic analysis and RNA based diagnostics of OTC deficiency.

Quantitative sequence analysis of FBN1 premature termination codons provides evidence for incomplete NMD in leukocytes

We improved, evaluated, and used Sanger sequencing for quantification of single nucleotide polymorphism (SNP) variants in transcripts and gDNA samples. This improved assay resulted in highly reproducible relative allele frequencies (e.g., for a heterozygous gDNA 50.0+/-1.4%, and for a missense mutation-bearing transcript 46.9+/-3.7%) with a lower detection limit of 3-9%. It provided excellent accuracy and linear correlation between expected and observed relative allele frequencies. This sequencing assay, which can also be used for the quantification of copy number variations (CNVs), methylations, mosaicisms, and DNA pools, enabled us to analyze transcripts of the FBN1 gene in fibroblasts and blood samples of patients with suspected Marfan syndrome not only qualitatively but also quantitatively. We report a total of 18 novel and 19 known FBN1 sequence variants leading to a premature termination codon (PTC), 26 of which we analyzed by quantitative sequencing both at gDNA and cDNA levels. The relative amounts of PTC-containing FBN1 transcripts in fresh and PAXgene-stabilized blood samples were significantly higher (33.0+/-3.9% to 80.0+/-7.2%) than those detected in affected fibroblasts with inhibition of nonsense-mediated mRNA decay (NMD) (11.0+/-2.1% to 25.0+/-1.8%), whereas in fibroblasts without NMD inhibition no mutant alleles could be detected. These results provide evidence for incomplete NMD in leukocytes and have particular importance for RNA-based analyses not only in FBN1 but also in other genes.

Systematic Functional Analysis of BicD Serine Phosphorylation and Intragenic Suppression of a Female Sterile Allele of BicD

Protein phosphorylation is involved in posttranslational control of essentially all biological processes. Using mass spectrometry, recent analyses of whole phosphoproteomes led to the identification of numerous new phosphorylation sites. However, the function of most of these sites remained unknown. We chose the Drosophila Bicaudal-D protein to estimate the importance of individual phosphorylation events. Being involved in different cellular processes, BicD is required for oocyte determination, for RNA transport during oogenesis and embryogenesis, and for photoreceptor nuclei migration in the developing eye. The numerous roles of BicD and the available evidence for functional importance of BicD phosphorylation led us to identify eight phosphorylation sites of BicD, and we tested a total of 14 identified and suspected phosphoserine residues for their functional importance in vivo in flies. Surprisingly, all these serines turned out to be dispensable for providing sufficient basal BicD activity for normal growth and development. However, in a genetically sensitized background where the BicD(A40V) protein variant provides only partial activity, serine 103 substitutions are not neutral anymore, but show surprising differences. The S103D substitution completely inactivates the protein, whereas S103A behaves neutral, and the S103F substitution, isolated in a genetic screen, restores BicD(A40V) function. Our results suggest that many BicD phosphorylation events may either be fortuitous or play a modulating function as shown for Ser(103). Remarkably, amongst the Drosophila serines we found phosphorylated, Ser(103) is the only one that is fully conserved in mammalian BicD.

RNA interference-mediated c-MYC inhibition prevents cell growth and decreases sensitivity to radio- and chemotherapy in childhood medulloblastoma cells

BACKGROUND: With current treatment strategies, nearly half of all medulloblastoma (MB) patients die from progressive tumors. Accordingly, the identification of novel therapeutic strategies remains a major goal. Deregulation of c-MYC is evident in numerous human cancers. In MB, over-expression of c-MYC has been shown to cause anaplasia and correlate with unfavorable prognosis. METHODS: To study the role of c-MYC in MB biology, we down-regulated c-MYC expression by using small interfering RNA (siRNA) and investigated changes in cellular proliferation, cell cycle analysis, apoptosis, telomere maintenance, and response to ionizing radiation (IR) and chemotherapeutics in a representative panel of human MB cell lines expressing different levels of c-MYC (DAOY wild-type, DAOY transfected with the empty vector, DAOY transfected with c-MYC, D341, and D425). RESULTS: siRNA-mediated c-MYC down-regulation resulted in an inhibition of cellular proliferation and clonogenic growth, inhibition of G1-S phase cell cycle progression, and a decrease in human telomerase reverse transcriptase (hTERT) expression and telomerase activity. On the other hand, down-regulation of c-MYC reduced apoptosis and decreased the sensitivity of human MB cells to IR, cisplatin, and etoposide. This effect was more pronounced in DAOY cells expressing high levels of c-MYC when compared with DAOY wild-type or DAOY cells transfected with the empty vector. CONCLUSION: In human MB cells, in addition to its roles in growth and proliferation, c-MYC is also a potent inducer of apoptosis. Therefore, targeting c-MYC might be of therapeutic benefit when used sequentially with chemo- and radiotherapy rather than concomitantly.

CD4 count slope and mortality in HIV-infected patients on antiretroviral therapy: multicohort analysis from South Africa

BACKGROUND In many resource-limited settings monitoring of combination antiretroviral therapy (cART) is based on the current CD4 count, with limited access to HIV RNA tests or laboratory diagnostics. We examined whether the CD4 count slope over 6 months could provide additional prognostic information. METHODS We analyzed data from a large multicohort study in South Africa, where HIV RNA is routinely monitored. Adult HIV-positive patients initiating cART between 2003 and 2010 were included. Mortality was analyzed in Cox models; CD4 count slope by HIV RNA level was assessed using linear mixed models. RESULTS About 44,829 patients (median age: 35 years, 58% female, median CD4 count at cART initiation: 116 cells/mm) were followed up for a median of 1.9 years, with 3706 deaths. Mean CD4 count slopes per week ranged from 1.4 [95% confidence interval (CI): 1.2 to 1.6] cells per cubic millimeter when HIV RNA was <400 copies per milliliter to -0.32 (95% CI: -0.47 to -0.18) cells per cubic millimeter with >100,000 copies per milliliter. The association of CD4 slope with mortality depended on current CD4 count: the adjusted hazard ratio (aHRs) comparing a >25% increase over 6 months with a >25% decrease was 0.68 (95% CI: 0.58 to 0.79) at <100 cells per cubic millimeter but 1.11 (95% CI: 0.78 to 1.58) at 201-350 cells per cubic millimeter. In contrast, the aHR for current CD4 count, comparing >350 with <100 cells per cubic millimeter, was 0.10 (95% CI: 0.05 to 0.20). CONCLUSIONS Absolute CD4 count remains a strong risk for mortality with a stable effect size over the first 4 years of cART. However, CD4 count slope and HIV RNA provide independently added to the model.

OMA and OPA—Software-Supported Mass Spectra Analysis of Native and Modified Nucleic Acids

The platform-independent software package consisting of the oligonucleotide mass assembler (OMA) and the oligonucleotide peak analyzer (OPA) was created to support the analysis of oligonucleotide mass spectra. It calculates all theoretically possible fragments of a given input sequence and annotates it to an experimental spectrum, thus, saving a large amount of manual processing time. The software performs analysis of precursor and product ion spectra of oligonucleotides and their analogues comprising user-defined modifications of the backbone, the nucleobases, or the sugar moiety, as well as adducts with metal ions or drugs. The ability to expand the library of building blocks and to implement individual structural variations makes it extremely useful for supporting the analysis of therapeutically active compounds. The functionality of the software tool is demonstrated on the examples of a platinated doublestranded oligonucleotide and a modified RNA sequence. Experiments also reveal the unique dissociation behavior of platinated higher-order DNA structures.

RNA-dependent genome processing during nuclear differentiation: the model systems of stichotrichous ciliates

We introduce ciliated protozoa, and more specifically the stichotrichous ciliates Oxytricha and Stylonychia, as biological model systems for the analysis of programmed DNA-reorganization processes during nuclear differentiation. These include DNA excision, DNA elimination, reordering of gene segments and specific gene amplification. We show that small nuclear RNAs specify DNA sequences to be excised or retained, but also discuss the need for a RNA template molecule derived from the parental nucleus for these processes. This RNA template guides reordering of gene segments to become functional genes and determines gene copy number in the differentiated nucleus. Since the template is derived from the parental macronucleus, gene reordering and DNA amplification are inherited in a non-Mendelian epigenetic manner.

Prediction of outcome in patients with low-grade squamous intraepithelial lesions by fluorescence in situ hybridization analysis of human papillomavirus, TERC, and MYC

BACKGROUND Cytology is an excellent method with which to diagnose preinvasive lesions of the uterine cervix, but it suffers from limited specificity for clinically significant lesions. Supplementary methods might predict the natural course of the detected lesions. The objective of the current study was to test whether a multicolor fluorescence in situ hybridization (FISH) assay might help to stratify abnormal results of Papanicolaou tests. METHODS A total of 219 liquid-based cytology specimens of low-grade squamous intraepithelial lesions (LSIL), 49 atypical squamous cells of undetermined significance (ASCUS) specimens, 52 high-grade squamous intraepithelial lesion (HSIL) specimens, and 50 normal samples were assessed by FISH with probes for the human papillomavirus (HPV), MYC, and telomerase RNA component (TERC). Subtyping of HPV by polymerase chain reaction (PCR) was performed in a subset of cases (n=206). RESULTS There was a significant correlation found between HPV detection by FISH and PCR (P<.0001). In patients with LSILs, the presence of HPV detected by FISH was significantly associated with disease progression (P<.0001). An increased MYC and/or TERC gene copy number (>2 signals in>10% of cells) prevailed in 43% of ASCUS specimens and was more frequent in HSIL (85%) than in LSIL (33%) (HSIL vs LSIL: P<.0001). Increased TERC gene copy number was significantly correlated with progression of LSIL (P<.01; odds ratio, 7.44; area under the receiver operating characteristic curve, 0.73; positive predictive value, 0.30; negative predictive value, 0.94) CONCLUSIONS: The detection of HPV by FISH analysis is feasible in liquid-based cytology and is significantly correlated with HPV analysis by PCR. The analysis of TERC gene copy number may be useful for risk stratification in patients with LSIL.

Next-generation sequencing of HIV-1 RNA genomes: determination of error rates and minimizing artificial recombination.

Next-generation sequencing (NGS) is a valuable tool for the detection and quantification of HIV-1 variants in vivo. However, these technologies require detailed characterization and control of artificially induced errors to be applicable for accurate haplotype reconstruction. To investigate the occurrence of substitutions, insertions, and deletions at the individual steps of RT-PCR and NGS, 454 pyrosequencing was performed on amplified and non-amplified HIV-1 genomes. Artificial recombination was explored by mixing five different HIV-1 clonal strains (5-virus-mix) and applying different RT-PCR conditions followed by 454 pyrosequencing. Error rates ranged from 0.04-0.66% and were similar in amplified and non-amplified samples. Discrepancies were observed between forward and reverse reads, indicating that most errors were introduced during the pyrosequencing step. Using the 5-virus-mix, non-optimized, standard RT-PCR conditions introduced artificial recombinants in a fraction of at least 30% of the reads that subsequently led to an underestimation of true haplotype frequencies. We minimized the fraction of recombinants down to 0.9-2.6% by optimized, artifact-reducing RT-PCR conditions. This approach enabled correct haplotype reconstruction and frequency estimations consistent with reference data obtained by single genome amplification. RT-PCR conditions are crucial for correct frequency estimation and analysis of haplotypes in heterogeneous virus populations. We developed an RT-PCR procedure to generate NGS data useful for reliable haplotype reconstruction and quantification.

Genome-wide analysis of genetic and epigenetic control of programmed DNA deletion

During the development of the somatic genome from the Paramecium germline genome the bulk of the copies of ∼45 000 unique, internal eliminated sequences (IESs) are deleted. IES targeting is facilitated by two small RNA (sRNA) classes: scnRNAs, which relay epigenetic information from the parental nucleus to the developing nucleus, and iesRNAs, which are produced and used in the developing nucleus. Why only certain IESs require sRNAs for their removal has been enigmatic. By analyzing the silencing effects of three genes: PGM (responsible for DNA excision), DCL2/3 (scnRNA production) and DCL5 (iesRNA production), we identify key properties required for IES elimination. Based on these results, we propose that, depending on the exact combination of their lengths and end bases, some IESs are less efficiently recognized or excised and have a greater requirement for targeting by scnRNAs and iesRNAs. We suggest that the variation in IES retention following silencing of DCL2/3 is not primarily due to scnRNA density, which is comparatively uniform relative to IES retention, but rather the genetic properties of IESs. Taken together, our analyses demonstrate that in Paramecium the underlying genetic properties of developmentally deleted DNA sequences are essential in determining the sensitivity of these sequences to epigenetic control.

Characterization of the potential role for RNA-binding protein FUS in DNA damage response: A quantitative proteomic approach

FUS/TLS (fused in sarcoma/translocated in liposarcoma) is a ubiquitously expressed protein of the hnRNP family, that has been discovered as fused to transcription factors in several human sarcomas and found in protein aggregates in neurons of patients with an inherited form of Amyotrophic Lateral Sclerosis [Vance C. et al., 2009]. FUS is a 53 kDa nuclear protein that contains structural domains, such as a RNA Recognition Motif (RRM) and a zinc finger motif, that give to FUS the ability to bind to both RNA and DNA sequences. It has been implicated in a variety of cellular processes, such as pre-mRNA splicing, miRNA processing, gene expression control and transcriptional regulation [Fiesel FC. and Kahle PJ., 2011]. Moreover, some evidences link FUS to genome stability control and DNA damage response: mice lacking FUS are hypersensitive to ionizing radiation (IR) and show high levels of chromosome instability and, in response to double-strand breaks, FUS is phosphorylated by the protein kinase ATM [Kuroda M. et al., 2000; Hicks GG. et al., 2000; Gardiner M. et al., 2008]. Furthermore, preliminary results of mass spectrometric identification of FUS interacting proteins in HEK293 cells, expressing a recombinant flag-tagged FUS protein, highlighted the interactions with proteins involved in DNA damage response, such as DNA-PK, XRCC-5/-6, and ERCC-6, raising the possibilities that FUS is involved in this pathway, even though its role still needs to be clarified. This study aims to investigate the biological roles of FUS in human cells and in particular the putative role in DNA damage response through the characterization of the proteomic profile of the neuroblastoma cell line SH-SY5Y upon FUS inducible depletion, by a quantitative proteomic approach. The SH-SY5Y cell line that will be used in this study expresses, in presence of tetracycline, a shRNA that targets FUS mRNA, leading to FUS protein depletion (SH-SY5Y FUS iKD cells). To quantify changes in proteins expression levels a SILAC strategy (Stable Isotope Labeling by Amino acids in Cell culture) will be conducted on SH-SY5Y FUS iKD cells and a control SH-SY5Y cell line (that expresses a mock shRNA) and the relative changes in proteins levels will be evaluated after five and seven days upon FUS depletion, by nanoliquid chromatography coupled to tandem mass spectrometry (nLC-MS/MS) and bioinformatics analysis. Preliminary experiments demonstrated that the SH-SY5Y FUS iKD cells, when subjected to genotoxic stress (high dose of IR), upon inducible depletion of FUS, showed a increased phosphorylation of gH2AX with respect to control cells, suggesting an higher activation of the DNA damage response.

Expression of the vault RNA protects cells from undergoing apoptosis

Non-protein-coding RNAs are a functionally versatile class of transcripts found in all domains of life exerting their biological role at the RNA level. Recently, we demonstrated that the vault-associated RNAs (vtRNAs) were significantly up-regulated in human B cells upon Epstein-Barr virus (EBV) infection [1,2]. vtRNAs are an integral part of the vault complex, a huge and evolutionarily conserved cytoplasmic ribonucleoprotein complex. The major vault protein (MVP) is the main structural component of the complex while vtRNA accounts for only 5% of its mass. Very little is known about the function(s) of the vtRNAs or the vault complex. In particular the role and significance of the previously observed vtRNA up-regulation upon EBV infection remained unclear. We individually expressed EBV-encoded genes in B cells and found the latent membrane protein 1 (LMP1) as trigger for vtRNA up-regulation. To unravel a putative functional interconnection between vtRNA expression and EBV infection, we ectopically expressed vtRNA1-1 in human B cells and observed an improved viral establishment. Furthermore, expression of vtRNA1-1 but not of the other vtRNA paralogs protected cells from undergoing apoptosis. Knock-down of MVP had no effect on these phenotypes thus revealing the vtRNA and not the vault complex to contribute to the enhanced EBV establishment and apoptosis resistance. Mutational analysis highlighted the central domain of the vtRNA to be involved in the anti-apoptotic effect. Ongoing research aims at characterizing the target of vtRNA1-1 in the apoptotic pathway. In summary, our data reveal a crucial cellular function for the so far elusive RNA biology of the vtRNAs.

Stability and kinetics of nucleic acid triplexes with chimaeric DNA/RNA third strands

A series of chimaeric DNA/RNA triplex-forming oligonucleotides (TFOs) with identical base-sequence but varying sequential composition of the sugar residues were prepared. The structural, kinetic and thermodynamic properties of triplex formation with their corresponding double-helical DNA target were investigated by spectroscopic methods. Kinetic and thermodynamic data were obtained from analysis of non-equilibrium UV-melting- and annealing curves in the range of pH 5.1 to 6.7 in a 10 mM citrate/phosphate buffer containing 0.1M NaCl and 1 mM EDTA. It was found that already single substitutions of ribo- for deoxyribonucleotides in the TFOs greatly affect stability and kinetics of triplex formation in a strongly sequence dependent manner. Within the sequence context investigated, triplex stability was found to increase when deoxyribonucleotides were present at the 5'-side and ribonucleotides in the center of the TFO. Especially the substitution of thymidines for uridines in the TFO was found to accelerate both, the association and dissociation process, in a strongly position-dependent way. Differential structural information on triplexes and TFO single-strands was obtained from CD-spectroscopy and gel mobility experiments. Only minor changes were observed in the CD spectra of the triplexes at all pH values investigated, and the electrophoretic mobility was nearly identical in all cases, indicating a high degree of structural similarity. In contrast, the single-stranded TFOs showed high structural variability as determined in the same way. The results are discussed in the context of the design of TFOs for therapeutic or biochemical applications.

Syncytiotrophoblast vesicles show altered micro-RNA and haemoglobin content after ex-vivo perfusion of placentas with haemoglobin to mimic preeclampsia.

BACKGROUND Cell-free foetal haemoglobin (HbF) has been shown to play a role in the pathology of preeclampsia (PE). In the present study, we aimed to further characterize the harmful effects of extracellular free haemoglobin (Hb) on the placenta. In particular, we investigated whether cell-free Hb affects the release of placental syncytiotrophoblast vesicles (STBMs) and their micro-RNA content. METHODS The dual ex-vivo perfusion system was used to perfuse isolated cotyledons from human placenta, with medium alone (control) or supplemented with cell-free Hb. Perfusion medium from the maternal side of the placenta was collected at the end of all perfusion phases. The STBMs were isolated using ultra-centrifugation, at 10,000×g and 150,000×g (referred to as 10K and 150K STBMs). The STBMs were characterized using the nanoparticle tracking analysis, identification of surface markers and transmission electron microscopy. RNA was extracted and nine different micro-RNAs, related to hypoxia, PE and Hb synthesis, were selected for analysis by quantitative PCR. RESULTS All micro-RNAs investigated were present in the STBMs. Mir-517a, mir-141 and mir-517b were down regulated after Hb perfusion in the 10K STBMs. Furthermore, Hb was shown to be carried by the STBMs. CONCLUSION This study showed that Hb perfusion can alter the micro-RNA content of released STBMs. Of particular interest is the alteration of two placenta specific micro-RNAs; mir-517a and mir-517b. We have also seen that STBMs may function as carriers of Hb into the maternal circulation.

Simultaneous analysis of nuclear and mitochondrial DNA, mRNA and miRNA from backspatter from inside parts of firearms generated by shots at "triple contrast" doped ballistic models

When a firearm projectile hits a biological target a spray of biological material (e.g., blood and tissue fragments) can be propelled from the entrance wound back towards the firearm. This phenomenon has become known as "backspatter" and if caused by contact shots or shots from short distances traces of backspatter may reach, consolidate on, and be recovered from, the inside surfaces of the firearm. Thus, a comprehensive investigation of firearm-related crimes must not only comprise of wound ballistic assessment but also backspatter analysis, and may even take into account potential correlations between these emergences. The aim of the present study was to evaluate and expand the applicability of the "triple contrast" method by probing its compatibility with forensic analysis of nuclear and mitochondrial DNA and the simultaneous investigation of co-extracted mRNA and miRNA from backspatter collected from internal components of different types of firearms after experimental shootings. We demonstrate that "triple contrast" stained biological samples collected from the inside surfaces of firearms are amenable to forensic co-analysis of DNA and RNA and permit sequence analysis of the entire mtDNA displacement-loop, even for "low template" DNA amounts that preclude standard short tandem repeat DNA analysis. Our findings underscore the "triple contrast" method's usefulness as a research tool in experimental forensic ballistics.