{\it In vivo\/} magnetic resonance studies of experimental liver disease: Carbon tetrachloride hepatotoxicity and alcohol-induced fatty liver in rat

Magnetic resonance imaging (MRI) and magnetic resonance spectroscopy (MRS) were used to non-invasively determine if cirrhosis induced by carbon tetrachloride (CCl$\sb4$) and phospholipase-D (PLD) could be distinguished from fatty infiltration in rat. MRS localization and water suppression methods were developed, implemented and evaluated in terms of their application to in vivo proton NMR studies of experimental liver disease. MRS studies were also performed to quantitate fatty infiltration resulting from carbon tetrachloride (CCl$\sb4$) or alcohol (ethanol) administration and the MRS results were confirmed using biochemical total lipid analysis and histology. $\rm T\sb1$ weighted MR images acquired weekly, 48 hours post administration, demonstrated only a slight increase in overall liver intensity with CCl$\sb4$ or alcohol administration, which is consistent with previously reported results. The MR images were able to detect nodules resulting from CCl$\sb4$+PLD induced cirrhosis as hypointense regions, also consistent with previous reports. Localized in vivo water and lipid proton $\rm T\sb1$ relaxation time measurements were performed and demonstrated no statistically significant trends for either agent. In vivo proton spectra were also acquired using stimulated echo techniques to quantitatively follow the changes in liver lipid content. The changes in liver lipid content observed using MRS were verified by total lipid analysis using the Folch technique and histology. The in vivo $\rm T\sb1$ and lipid quantification data str inconsistent with the previous hypothesis that the changes in $\rm T\sb1$ weighted images were the result of increased "free" water content and, therefore, increased water $\rm T\sb1$ relaxation times. These data indicate that the long term changes are more likely the result of changes in lipid content. The data are also shown to agree with the accepted hypothesis that the time course and mechanism of fatty infiltration are different for CCl$\sb4$ and alcohol. The hypothesis that the lipids resulting from either protocol are from the same lipid fraction(s), presumably triglycerides, is also supported. And lastly, on the basis of MR images and quantitative MRS lipid information, it was shown that cirrhosis could be distinguished from fatty infiltration. ^

New methods for quantification and analysis of quantitative real-time polymerase chain reaction data

Quantitative real-time polymerase chain reaction (qPCR) is a sensitive gene quantitation method that has been widely used in the biological and biomedical fields. The currently used methods for PCR data analysis, including the threshold cycle (CT) method, linear and non-linear model fitting methods, all require subtracting background fluorescence. However, the removal of background fluorescence is usually inaccurate, and therefore can distort results. Here, we propose a new method, the taking-difference linear regression method, to overcome this limitation. Briefly, for each two consecutive PCR cycles, we subtracted the fluorescence in the former cycle from that in the later cycle, transforming the n cycle raw data into n-1 cycle data. Then linear regression was applied to the natural logarithm of the transformed data. Finally, amplification efficiencies and the initial DNA molecular numbers were calculated for each PCR run. To evaluate this new method, we compared it in terms of accuracy and precision with the original linear regression method with three background corrections, being the mean of cycles 1-3, the mean of cycles 3-7, and the minimum. Three criteria, including threshold identification, max R2, and max slope, were employed to search for target data points. Considering that PCR data are time series data, we also applied linear mixed models. Collectively, when the threshold identification criterion was applied and when the linear mixed model was adopted, the taking-difference linear regression method was superior as it gave an accurate estimation of initial DNA amount and a reasonable estimation of PCR amplification efficiencies. When the criteria of max R2 and max slope were used, the original linear regression method gave an accurate estimation of initial DNA amount. Overall, the taking-difference linear regression method avoids the error in subtracting an unknown background and thus it is theoretically more accurate and reliable. This method is easy to perform and the taking-difference strategy can be extended to all current methods for qPCR data analysis.^

Effect of continuous gastric enteral feeding on tonometric measurement of gastric intramucosal carbon dioxide levels in healthy volunteers

Purpose. The aim of this research was to evaluate the effect of enteral feeding on tonometric measurement of gastric regional carbon dioxide levels (PrCO2) in normal healthy volunteers. Design and methods. The sample included 12 healthy volunteers recruited by the University Clinical Research Center (UCRC). An air tonometry system monitored PrCO2 levels using a tonometer placed in the lumen of the stomach via orogastric intubation. PrCO2 was automatically measured and recorded every 10 minutes throughout the five hour study period. An oral dose of famotidine 40 mg was self-administered the evening prior to and the morning of the study. Instillation of Isocal® High Nitrogen (HN) was used for enteral feeding in hourly escalating doses of 0, 40, 60, and 80 ml/hr with no feeding during the fifth hour. Results . PrCO2 measurements at time 0 and 10 minutes (41.4 ± 6.5 and 41.8 ± 5.7, respectively) demonstrated biologic precision (Levene's Test statistic = 0.085, p-value 0.774). Biologic precision was lost between T130 and T140 40 when compared to baseline TO (Levene's Test statistic = 1.70, p-value 0.205; and 3.205, p-value 0.042, respectively) and returned to non-significant levels between T270 and T280 (Levene's Test statistic = 3.083, p-value 0.043; and 2.307, p-value 0.143, respectively). Isocal® HN significantly affected the biologic accuracy of PrCO2 measurements (repeated measures ANOVA F 4.91, p-value <0.001). After 20 minutes of enteral feeding at 40 ml/hr, PrCO2 significantly increased (41.4 ± 6.5 to 46.6 ± 4.25, F = 5.4, p-value 0.029). Maximum variance from baseline (41.4 ± 6.5 to 61.3 ± 15.2, F = 17.22, p-value <0.001) was noted after 30 minutes of Isocal® HN at 80 ml/hr or 210 minutes from baseline. The significant elevations in PrCO2 continued throughout the study. Sixty minutes after discontinuation of enteral feeding, PrCO2 remained significantly elevated from baseline (41.4 ± 6.5 to 51.8 ± 9.2, F = 10.15, p-value 0.004). Conclusion. Enteral feeding with Isocal® HN significantly affects the precision and accuracy of PrCO2 measurements in healthy volunteers. ^

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.