Comparison of amplification methods for transcriptomic analyses of low abundance prokaryotic RNA sources

Microarrays have established as instrumental for bacterial detection, identification, and genotyping as well as for transcriptomic studies. For gene expression analyses using limited numbers of bacteria (derived from in vivo or ex vivo origin, for example), RNA amplification is often required prior to labeling and hybridization onto microarrays. Evaluation of the fidelity of the amplification methods is crucial for the robustness and reproducibility of microarray results. We report here the first utilization of random primers and the highly processive Phi29 phage polymerase to amplify material for transcription profiling analyses. We compared two commercial amplification methods (GenomiPhi and MessageAmp kits) with direct reverse-transcription as the reference method, focusing on the robustness of mRNA quantification using either microarrays or quantitative RT-PCR. Both amplification methods using either poly-A tailing followed by in vitro transcription, or direct strand displacement polymerase, showed appreciable linearity. Strand displacement technique was particularly affordable compared to in vitro transcription-based (IVT) amplification methods and consisted in a single tube reaction leading to high amplification yields. Real-time measurements using low-, medium-, and highly expressed genes revealed that this simple method provided linear amplification with equivalent results in terms of relative messenger abundance as those obtained by conventional direct reverse-transcription.

Effect of echo time pair selection on quantitative analysis for adrenal tumor characterization with in-phase and opposed-phase MR imaging: initial experience

PURPOSE: To determine the effect of two pairs of echo times (TEs) for in-phase (IP) and opposed-phase (OP) 3.0-T magnetic resonance (MR) imaging on (a) quantitative analysis prospectively in a phantom study and (b) diagnostic accuracy retrospectively in a clinical study of adrenal tumors, with use of various reference standards in the clinical study. MATERIALS AND METHODS: A fat-saline phantom was used to perform IP and OP 3.0-T MR imaging for various fat fractions. The institutional review board approved this HIPAA-compliant study, with waiver of informed consent. Single-breath-hold IP and OP 3.0-T MR images in 21 patients (14 women, seven men; mean age, 63 years) with 23 adrenal tumors (16 adenomas, six metastases, one adrenocortical carcinoma) were reviewed. The MR protocol involved two acquisition schemes: In scheme A, the first OP echo (approximately 1.5-msec TE) and the second IP echo (approximately 4.9-msec TE) were acquired. In scheme B, the first IP echo (approximately 2.4-msec TE) and the third OP echo (approximately 5.8-msec TE) were acquired. Quantitative analysis was performed, and analysis of variance was used to test for differences between adenomas and nonadenomas. RESULTS: In the phantom study, scheme B did not enable discrimination among voxels that had small amounts of fat. In the clinical study, no overlap in signal intensity (SI) index values between adenomas and nonadenomas was seen (P < .05) with scheme A. However, with scheme B, no overlap in the adrenal gland SI-to-liver SI ratio between adenomas and nonadenomas was seen (P < .05). With scheme B, no overlap in adrenal gland SI index-to-liver SI index ratio between adenomas and nonadenomas was seen (P < .05). CONCLUSION: This initial experience indicates SI index is the most reliable parameter for characterization of adrenal tumors with 3.0-T MR imaging when obtaining OP echo before IP echo. When acquiring IP echo before OP echo, however, nonadenomas can be mistaken as adenomas with use of the SI index value.

Quantitative stress echocardiography in coronary artery disease using contrast-based myocardial blood flow measurements: prospective comparison with coronary angiography

AIM: To test whether quantitative stress echocardiography using contrast-based myocardial blood flow (MBF, ml x min(-1) x g(-1)) measurements can detect coronary artery disease in humans. METHODS: 48 patients eligible for pharmacological stress testing by myocardial contrast echocardiography (MCE) and willing to undergo subsequent coronary angiography were prospectively enrolled in the study. Baseline and adenosine-induced (140 microg x kg(-1) x min(-1)) hyperaemic MBF was analysed according to a three-coronary-artery-territory model. Vascular territories were categorised into three groups with increasing stenosis severity defined as percentage diameter reduction by quantitative coronary angiography. RESULTS: Myocardial blood flow reserve (MBFR)-that is, the ratio of hyperaemic to baseline MBF, was obtained in 128 (89%) territories. Mean (SD) baseline MBF was 1.073 (0.395) ml x min(-1) x g(-1) and did not differ between territories supplied by coronary arteries with mild (<50% stenosis), moderate (50%-74% stenosis) or severe (>or=75% stenosis) disease. Mean (SD) hyperaemic MBF and MBFR were 2.509 (1.078) ml x min(-1) x g(-1) and 2.54 (1.03), respectively, and decreased linearly (r2 = 0.21 and r2 = 0.39) with stenosis severity. ROC analysis revealed that a territorial MBFR <1.94 detected >or=50% stenosis with 89% sensitivity and 92% specificity. CONCLUSION: Quantitative stress testing based on MBF measurements derived from contrast echocardiography is a new method for the non-invasive and reliable assessment of coronary artery disease in humans.

STATISTICAL METHODS FOR RARE AND COMMON VARIANT ASSOCIATION STUDIES

Complex human diseases are a major challenge for biological research. The goal of my research is to develop effective methods for biostatistics in order to create more opportunities for the prevention and cure of human diseases. This dissertation proposes statistical technologies that have the ability of being adapted to sequencing data in family-based designs, and that account for joint effects as well as gene-gene and gene-environment interactions in the GWA studies. The framework includes statistical methods for rare and common variant association studies. Although next-generation DNA sequencing technologies have made rare variant association studies feasible, the development of powerful statistical methods for rare variant association studies is still underway. Chapter 2 demonstrates two adaptive weighting methods for rare variant association studies based on family data for quantitative traits. The results show that both proposed methods are robust to population stratification, robust to the direction and magnitude of the effects of causal variants, and more powerful than the methods using weights suggested by Madsen and Browning [2009]. In Chapter 3, I extended the previously proposed test for Testing the effect of an Optimally Weighted combination of variants (TOW) [Sha et al., 2012] for unrelated individuals to TOW &ndash F, TOW for Family &ndash based design. Simulation results show that TOW &ndash F can control for population stratification in wide range of population structures including spatially structured populations, is robust to the directions of effect of causal variants, and is relatively robust to percentage of neutral variants. In GWA studies, this dissertation consists of a two &ndash locus joint effect analysis and a two-stage approach accounting for gene &ndash gene and gene &ndash environment interaction. Chapter 4 proposes a novel two &ndash stage approach, which is promising to identify joint effects, especially for monotonic models. The proposed approach outperforms a single &ndash marker method and a regular two &ndash stage analysis based on the two &ndash locus genotypic test. In Chapter 5, I proposed a gene &ndash based two &ndash stage approach to identify gene &ndash gene and gene &ndash environment interactions in GWA studies which can include rare variants. The two &ndash stage approach is applied to the GAW 17 dataset to identify the interaction between KDR gene and smoking status.

A new method to monitor visual field defects caused by photoreceptor degeneration by quantitative optical coherence tomography

PURPOSE: To correlate the dimension of the visual field (VF) tested by Goldman kinetic perimetry with the extent of visibility of the highly reflective layer between inner and outer segments of photoreceptors (IOS) seen in optical coherence tomography (OCT) images in patients with retinitis pigmentosa (RP). METHODS: In a retrospectively designed cross-sectional study, 18 eyes of 18 patients with RP were examined with OCT and Goldmann perimetry using test target I4e and compared with 18 eyes of 18 control subjects. A-scans of raw scan data of Stratus OCT images (Carl Zeiss Meditec, AG, Oberkochen, Germany) were quantitatively analyzed for the presence of the signal generated by the highly reflective layer between the IOS in OCT images. Starting in the fovea, the distance to which this signal was detectable was measured. Visual fields were analyzed by measuring the distance from the center point to isopter I4e. OCT and visual field data were analyzed in a clockwise fashion every 30 degrees , and corresponding measures were correlated. RESULTS: In corresponding alignments, the distance from the center point to isopter I4e and the distance to which the highly reflective signal from the IOS can be detected correlate significantly (r = 0.75, P < 0.0001). The greater the distance in VF, the greater the distance measured in OCT. CONCLUSIONS: The authors hypothesize that the retinal structure from which the highly reflective layer between the IOS emanates is of critical importance for visual and photoreceptor function. Further research is warranted to determine whether this may be useful as an objective marker of progression of retinal degeneration in patients with RP.

Quantitative OCT analysis of idiopathic perifoveal telangiectasia

PURPOSE: To identify and quantitate specific changes in optical coherence tomography (OCT) images of patients with type 2 idiopathic perifoveal telangiectasia (IPT). METHODS: In a prospectively designed, observational, case-control study, 28 eyes of 14 consecutive patients with IPT were examined with OCT and compared with eyes of 14 unaffected control subjects. Light reflectivity profiles of raw scan data of OCT images were quantitatively analyzed for differences in distance between different retinal reflectivity layers and their respective reflectivities. Maculae were examined in four separate regions: (1) central fovea, (2) nasal perifovea, (3) temporal perifovea, and (4) outside the fovea. RESULTS: Retinal thinning, shortening of the photoreceptor outer segments and loss of reflectivity of the photoreceptor ellipsoid region were found in the central foveal region as well as the nasal and temporal perifoveal regions in eyes with IPT. In addition, increased reflectivity of the outer nuclear layer was found in a sharply demarcated area of the inferotemporal perifoveal region in all affected eyes. Retinal tissue located more than 2000 mum away from the foveola was indistinguishable from that in normal eyes. CONCLUSIONS: Quantitative OCT analysis shows unique and specific changes in the photoreceptors of the central macula in IPT which can be detected from first clinical presentation. These changes may be of use as an additional diagnostic tool. Correlation of the findings in the outer nuclear layer with histologic studies may help identify the nature of the reflectivity increase and define more clearly the type of damage sustained by the photoreceptors in this condition.

Quantitative isolation of mouse and human intestinal intraepithelial lymphocytes by elutriation centrifugation

Intestinal intraepithelial lymphocytes (IEL) are specialized subsets of T cells with distinct functional capacities. While some IEL subsets are circulating, others such as CD8alphaalpha TCRalphabeta IEL are believed to represent non-circulating resident T cell subsets [Sim, G.K., Intraepithelial lymphocytes and the immune system. Adv. Immunol., 1995. 58: 297-343.]. Current methods to obtain enriched preparations of intraepithelial lymphocytes are mostly based on Percoll density gradient or magnetic bead-based technologies [Lundqvist, C., et al., Isolation of functionally active intraepithelial lymphocytes and enterocytes from human small and large intestine. J. Immunol. Methods, 1992. 152(2): 253-263.]. However, these techniques are hampered by a generally low yield of isolated cells, and potential artifacts due to the interference of the isolation procedure with subsequent functional assays, in particular, when antibodies against cell surface markers are required. Here we describe a new method for obtaining relatively pure populations of intestinal IEL (55-75%) at a high yield (>85%) by elutriation centrifugation. This technique is equally suited for the isolation and enrichment of intraepithelial lymphocytes of both mouse and human origin. Time requirements for fractionating cell suspensions by elutriation centrifugation are comparable to Percoll-, or MACS-based isolation procedures. Hence, the substantially higher yield and the consistent robust enrichment for intraepithelial lymphocytes, together with the gentle treatment of the cells during elutriation that does not interfere with subsequent functional assays, are important aspects that are in favor of using this elegant technology to obtain unmanipulated, unbiased populations of intestinal intraepithelial lymphocytes, and, if desired, also of pure epithelial cells.

Simultaneous quantitative detection of relevant biomarkers in breast cancer by quantitative real-time PCR

The assessment of ERa, PgR and HER2 status is routinely performed today to determine the endocrine responsiveness of breast cancer samples. Such determination is usually accomplished by means of immunohistochemistry and in case of HER2 amplification by means of fluorescent in situ hybridization (FISH). The analysis of these markers can be improved by simultaneous measurements using quantitative real-time PCR (Qrt-PCR). In this study we compared Qrt-PCR results for the assessment of mRNA levels of ERa, PgR, and the members of the human epidermal growth factor receptor family, HER1, HER2, HER3 and HER4. The results were obtained in two independent laboratories using two different methods, SYBR Green I and TaqMan probes, and different primers. By linear regression we demonstrated a good concordance for all six markers. The quantitative mRNA expression levels of ERa, PgR and HER2 also strongly correlated with the respective quantitative protein expression levels prospectively detected by EIA in both laboratories. In addition, HER2 mRNA expression levels correlated well with gene amplification detected by FISH in the same biopsies. Our results indicate that both Qrt-PCR methods were robust and sensitive tools for routine diagnostics and consistent with standard methodologies. The developed simultaneous assessment of several biomarkers is fast and labor effective and allows optimization of the clinical decision-making process in breast cancer tissue and/or core biopsies.

Age modeling of young non-varved lake sediments: methods and limits. Examples from two lakes in Central Chile

High-resolution and highly precise age models for recent lake sediments (last 100–150 years) are essential for quantitative paleoclimate research. These are particularly important for sedimentological and geochemical proxies, where transfer functions cannot be established and calibration must be based upon the relation of sedimentary records to instrumental data. High-precision dating for the calibration period is most critical as it determines directly the quality of the calibration statistics. Here, as an example, we compare radionuclide age models obtained on two high-elevation glacial lakes in the Central Chilean Andes (Laguna Negra: 33°38′S/70°08′W, 2,680 m a.s.l. and Laguna El Ocho: 34°02′S/70°19′W, 3,250 m a.s.l.). We show the different numerical models that produce accurate age-depth chronologies based on 210Pb profiles, and we explain how to obtain reduced age-error bars at the bottom part of the profiles, i.e., typically around the end of the 19th century. In order to constrain the age models, we propose a method with five steps: (i) sampling at irregularly-spaced intervals for 226Ra, 210Pb and 137Cs depending on the stratigraphy and microfacies, (ii) a systematic comparison of numerical models for the calculation of 210Pb-based age models: constant flux constant sedimentation (CFCS), constant initial concentration (CIC), constant rate of supply (CRS) and sediment isotope tomography (SIT), (iii) numerical constraining of the CRS and SIT models with the 137Cs chronomarker of AD 1964 and, (iv) step-wise cross-validation with independent diagnostic environmental stratigraphic markers of known age (e.g., volcanic ash layer, historical flood and earthquakes). In both examples, we also use airborne pollutants such as spheroidal carbonaceous particles (reflecting the history of fossil fuel emissions), excess atmospheric Cu deposition (reflecting the production history of a large local Cu mine), and turbidites related to historical earthquakes. Our results show that the SIT model constrained with the 137Cs AD 1964 peak performs best over the entire chronological profile (last 100–150 years) and yields the smallest standard deviations for the sediment ages. Such precision is critical for the calibration statistics, and ultimately, for the quality of the quantitative paleoclimate reconstruction. The systematic comparison of CRS and SIT models also helps to validate the robustness of the chronologies in different sections of the profile. Although surprisingly poorly known and under-explored in paleolimnological research, the SIT model has a great potential in paleoclimatological reconstructions based on lake sediments

Quality controls for two heel bone ultrasounds used in the Swiss Evaluation of the Methods of Measurement of Osteoporotic Fracture Risk Study

Because of the important morbidity and mortality associated with osteoporosis, it is essential to detect subjects at risk by screening methods, such as bone quantitative ultrasounds (QUSs). Several studies showed that QUS could predict fractures. None, however, compared prospectively different QUS devices, and few data of quality controls (QCs) have been published. The Swiss Evaluation of the Methods of Measurement of Osteoporotic Fracture Risk is a prospective multicenter study that compared three QUSs for the assessment of hip fracture risk in a population of 7609 women age >/=70 yr. Because the inclusion phase lasted 20 mo, and because 10 centers participated in this study, QC became a major issue. We therefore developed a QC procedure to assess the stability and precision of the devices, and for their cross-calibration. Our study focuses on the two heel QUSs. The water bath system (Achilles+) had a higher precision than the dry system (Sahara). The QC results were highly dependent on temperature. QUS stability was acceptable, but Sahara must be calibrated regularly. A sufficient homogeneity among all the Sahara devices could be demonstrated, whereas significant differences were found among the Achilles+ devices. For speed of sound, 52% of the differences among the Achilles+ was explained by the water s temperature. However, for broadband ultrasound attenuation, a maximal difference of 23% persisted after adjustment for temperature. Because such differences could influence measurements in vivo, it is crucial to develop standardized phantoms to be used in prospective multicenter studies.

Quantitative comparison of 3 enamel-stripping devices in vitro: how precisely can we strip teeth?

INTRODUCTION In this in-vitro study, we aimed to investigate the predictability of the expected amount of stripping using 3 common stripping devices on premolars. METHODS One hundred eighty extracted premolars were mounted and aligned in silicone. Tooth mobility was tested with Periotest (Medizintechnik Gulden, Modautal, Germany) (8.3 ± 2.8 units). The selected methods for interproximal enamel reduction were hand-pulled strips (Horico, Hapf Ringleb & Company, Berlin, Germany), oscillating segmental disks (O-drive-OD 30; KaVo Dental, Biberach, Germany), and motor-driven abrasive strips (Orthofile; SDC Switzerland, Lugano-Grancia, Switzerland). With each device, the operator intended to strip 0.1, 0.2, 0.3, or 0.4 mm on the mesial side of 15 teeth. The teeth were scanned before and after stripping with a 3-dimensional laser scanner. Superposition and measurement of stripped enamel on the most mesial point of the tooth were conducted with Viewbox software (dHal Software, Kifissia, Greece). The Wilcoxon signed rank test and the Kruskal-Wallis test were applied; statistical significance was set at alpha ≤ 0.05. RESULTS Large variations between the intended and the actual amounts of stripped enamel, and between stripping procedures, were observed. Significant differences were found at 0.1 mm of intended stripping (P ≤ 0.05) for the hand-pulled method and at 0.4 mm of intended stripping (P ≤ 0.001 to P = 0.05) for all methods. For all scenarios of enamel reduction, the actual amount of stripping was less than the predetermined and expected amount of stripping. The Kruskal-Wallis analysis showed no significant differences between the 3 methods. CONCLUSIONS There were variations in the stripped amounts of enamel, and the stripping technique did not appear to be a significant predictor of the actual amount of enamel reduction. In most cases, actual stripping was less than the intended amount of enamel reduction.

Quantitative contrast echocardiographic assessment of collateral derived myocardial perfusion during elective coronary angioplasty

OBJECTIVE To determine whether myocardial contrast echocardiography can be used to quantify collateral derived myocardial flow in humans. METHODS In 25 patients undergoing coronary angioplasty, a collateral flow index (CFI) was determined using intracoronary wedge pressure distal to the stenosis to be dilated, with simultaneous mean aortic pressure measurements. During balloon occlusion, echo contrast was injected into both main coronary arteries simultaneously. Echocardiography of the collateral receiving myocardial area was performed. The time course of myocardial contrast enhancement in images acquired at end diastole was quantified by measuring pixel intensities (256 grey units) within a region of interest. Perfusion variables, such as background subtracted peak pixel intensity and contrast transit rate, were obtained from a fitted gamma variate curve. RESULTS 16 patients had a left anterior descending coronary artery stenosis, four had a left circumflex coronary artery stenosis, and five had a right coronary artery stenosis. The mean (SD) CFI was 19 (12)% (range 0-47%). Mean contrast transit rate was 11 (8) seconds. In 17 patients, a significant collateral contrast effect was observed (defined as peak pixel intensity more than the mean + 2 SD of background). Peak pixel intensity was linearly related to CFI in patients with a significant contrast effect (p = 0.002, r = 0.69) as well as in all patients (p = 0.0003, r = 0.66). CONCLUSIONS Collateral derived perfusion of myocardial areas at risk can be demonstrated using intracoronary echo contrast injections. The peak echo contrast effect is directly related to the magnitude of collateral flow.

Time courses and quantitative analysis of atrial fibrillation episode number and duration after circular plus linear left atrial lesions: trigger elimination or substrate modification: early or delayed cure?

OBJECTIVES We sought to analyze the time course of atrial fibrillation (AF) episodes before and after circular plus linear left atrial ablation and the percentage of patients with complete freedom from AF after ablation by using serial seven-day electrocardiograms (ECGs). BACKGROUND The curative treatment of AF targets the pathophysiological corner stones of AF (i.e., the initiating triggers and/or the perpetuation of AF). The pathophysiological complexity of both may not result in an "all-or-nothing" response but may modify number and duration of AF episodes. METHODS In patients with highly symptomatic AF, circular plus linear ablation lesions were placed around the left and right pulmonary veins, between the two circles, and from the left circle to the mitral annulus using the electroanatomic mapping system. Repetitive continuous 7-day ECGs administered before and after catheter ablation were used for rhythm follow-up. RESULTS In 100 patients with paroxysmal (n = 80) and persistent (n = 20) AF, relative duration of time spent in AF significantly decreased over time (35 +/- 37% before ablation, 26 +/- 41% directly after ablation, and 10 +/- 22% after 12 months). Freedom from AF stepwise increased in patients with paroxysmal AF and after 12 months measured at 88% or 74% depending on whether 24-h ECG or 7-day ECG was used. Complete pulmonary vein isolation was demonstrated in <20% of the circular lesions. CONCLUSIONS The results obtained in patients with AF treated with circular plus linear left atrial lesions strongly indicate that substrate modification is the main underlying pathophysiologic mechanism and that it results in a delayed cure instead of an immediate cure.

Quantitative 1H-magnetic resonance spectroscopy of human brain: Influence of composition and parameterization of the basis set in linear combination model-fitting

Localized short-echo-time (1)H-MR spectra of human brain contain contributions of many low-molecular-weight metabolites and baseline contributions of macromolecules. Two approaches to model such spectra are compared and the data acquisition sequence, optimized for reproducibility, is presented. Modeling relies on prior knowledge constraints and linear combination of metabolite spectra. Investigated was what can be gained by basis parameterization, i.e., description of basis spectra as sums of parametric lineshapes. Effects of basis composition and addition of experimentally measured macromolecular baselines were investigated also. Both fitting methods yielded quantitatively similar values, model deviations, error estimates, and reproducibility in the evaluation of 64 spectra of human gray and white matter from 40 subjects. Major advantages of parameterized basis functions are the possibilities to evaluate fitting parameters separately, to treat subgroup spectra as independent moieties, and to incorporate deviations from straightforward metabolite models. It was found that most of the 22 basis metabolites used may provide meaningful data when comparing patient cohorts. In individual spectra, sums of closely related metabolites are often more meaningful. Inclusion of a macromolecular basis component leads to relatively small, but significantly different tissue content for most metabolites. It provides a means to quantitate baseline contributions that may contain crucial clinical information.

Absolute quantitative real-time polymerase chain reaction for the measurement of human papillomavirus E7 mRNA in cervical cytobrush specimens.

BACKGROUND: Few reports of the utilization of an accurate, cost-effective means for measuring HPV oncogene transcripts have been published. Several papers have reported the use of relative quantitation or more expensive Taqman methods. Here, we report a method of absolute quantitative real-time PCR utilizing SYBR-green fluorescence for the measurement of HPV E7 expression in cervical cytobrush specimens. RESULTS: The construction of a standard curve based on the serial dilution of an E7-containing plasmid was the key for being able to accurately compare measurements between cervical samples. The assay was highly reproducible with an overall coefficient of variation of 10.4%. CONCLUSION: The use of highly reproducible and accurate SYBR-based real-time polymerase chain reaction (PCR) assays instead of performing Taqman-type assays allows low-cost, high-throughput analysis of viral mRNA expression. The development of such assays will help in refining the current screening programs for HPV-related carcinomas.