REMOVING TECHNICAL VARIABILITY IN RNA-SEQ DATA USING CONDITIONAL QUANTILE NORMALIZATION

The ability to measure gene expression on a genome-wide scale is one of the most promising accomplishments in molecular biology. Microarrays, the technology that first permitted this, were riddled with problems due to unwanted sources of variability. Many of these problems are now mitigated, after a decade’s worth of statistical methodology development. The recently developed RNA sequencing (RNA-seq) technology has generated much excitement in part due to claims of reduced variability in comparison to microarrays. However, we show RNA-seq data demonstrates unwanted and obscuring variability similar to what was first observed in microarrays. In particular, we find GC-content has a strong sample specific effect on gene expression measurements that, if left uncorrected, leads to false positives in downstream results. We also report on commonly observed data distortions that demonstrate the need for data normalization. Here we describe statistical methodology that improves precision by 42% without loss of accuracy. Our resulting conditional quantile normalization (CQN) algorithm combines robust generalized regression to remove systematic bias introduced by deterministic features such as GC-content, and quantile normalization to correct for global distortions.

Sendai virus RNA polymerase scanning for mRNA start sites at gene junctions

Mini-genomes expressing two reporter genes and a variable gene junction were used to study Sendai virus RNA polymerase (RdRp) scanning for the mRNA start signal of the downstream gene (gs2). We found that RdRp could scan the template efficiently as long as the initiating uridylate of gs2 (3' UCCCnnUUUC) was preceded by the conserved intergenic region (3' GAA) and the last 3 uridylates of the upstream gene end signal (ge1; 3' AUUCUUUUU). The end of the leader sequence (3' CUAAAA, which precedes gs1) could also be used for gene2 expression, but this sequence was considerably less efficient. Increasing the distance between ge1 and gs2 (up to 200 nt) led to the progressive loss of gene2 expression, in which half of gene2 expression was lost for each 70 nucleotides of intervening sequence. Beyond 200 nt, gene2 expression was lost more slowly. Our results suggest that there may be two populations of RdRp that scan at gene junctions, which can be distinguished by the efficiency with which they can scan the genome template for gs.

Photochemically induced RNA and DNA abasic sites

Two phosphoramidite building blocks were synthesized that can easily be deprotected by UV light to reveal natural abasic sites in oligoribonucleotides as well as in oligodeoxyribonucleotides. Another building block which releases a 2 ′-O-methylated abasic site upon UV radiation is also described.

Messenger RNA levels and binding sites of muscarinic acetylcholine receptors in gastrointestinal muscle layers from healthy dairy cows

Acetylcholine interacts with muscarinic receptors (M) to mediate gastrointestinal (GI) smooth muscle contractions. We have compared mRNA levels and binding sites of M(1)to M(5) in muscle tissues from fundus abomasi, pylorus, ileum, cecum, proximal loop of the ascending colon (PLAC), and external loop of the spiral colon (ELSC) of healthy dairy cows. The mRNA levels were measured by quantitative RT-PCR. The inhibition of [(3)H]-QNB (1-quinuclidinyl-[phenyl-4-(3)H]-benzilate) binding by M antagonists [atropine (M(1 - 5)), pirenzepine (M(1)), methoctramine (M(2)), 4-DAMP (M(3)), and tropicamide (M(4))] was used to identify receptors at the functional level. Maximal binding (B(max)) was determined through saturation binding with atropine as a competitor. The mRNA levels of M(1), M(2), M(3), and M(5) represented 0.2, 48, 50, and 1.8%, respectively, of the total M population, whereas mRNA of M(4) was undetectable. The mRNA levels of M(2) and of M(3) in the ileum were lower (P < 0.05) than in other GI locations, which were similar among each other. Atropine, pirenzepine, methoctramine, and 4-DAMP inhibited [(3)H]-QNB binding according to an either low- or high-affinity receptor pattern, whereas tropicamide had no effect on [(3)H]-QNB binding. The [(3)H]-QNB binding was dose-dependent and saturable. B(max) in fundus, pylorus, and PLAC was lower (P < 0.05) than in the ELSC, and in the pylorus lower (P < 0.05) than in the ileum. B(max) and mRNA levels were negatively correlated (r = -0.3; P < 0.05). In conclusion, densities of M are different among GI locations, suggesting variable importance of M for digestive functions along the GI tract.

Influence of beveling and ultrasound application on marginal adaptation of box-only Class II (slot) resin composite restorations

A laboratory study was performed to assess the influence of beveling the margins of cavities and the effects on marginal adaptation of the application of ultrasound during setting and initial light curing. After minimal access cavities had been prepared with an 80 microm diamond bur, 80 box-only Class II cavities were prepared mesially and distally in 40 extracted human molars using four different oscillating diamond coated instruments: (A) a U-shaped PCS insert as the non-beveled control (EMS), (B) Bevelshape (Intensiv), (C) SonicSys (KaVo) and (D) SuperPrep (KaVo). In groups B-D, the time taken for additional bevel finishing was measured. The cavities were filled with a hybrid composite material in three increments. Ultrasound was also applied to one cavity per tooth before and during initial light curing (10 seconds). The specimens were subjected to thermomechanical stress in a computer-controlled masticator device. Marginal quality was assessed by scanning electron microscopy and the results were compared statistically. The additional time required for finishing was B > D > C (p < or = 0.05). In all groups, thermomechanical loading resulted in a decrease in marginal quality. Beveling resulted in higher values for "continuous" margins compared with that of the unbeveled controls. The latter showed better marginal quality at the axial walls when ultrasound was used. Beveling seems essential for good marginal adaptation but requires more preparation time. The use of ultrasonic vibrations may improve the marginal quality of unbeveled fillings and warrants further investigation.