Mouse models and quantitative trait loci of susceptibility to bleomycin- and radiation-induced pulmonary fibrosis

Radiotherapy involving the thoracic cavity and chemotherapy with the drug bleomycin are both dose limited by the development of pulmonary fibrosis. From evidence that there is variation in the population in susceptibility to pulmonary fibrosis, and animal data, it was hypothesized that individual variation in susceptibility to bleomycin-induced, or radiation-induced, pulmonary fibrosis is, in part, genetically controlled. In this thesis a three generation mouse genetic model of C57BL/6J (fibrosis prone) and C3Hf/Kam (fibrosis resistant) mouse strains and F1 and F2 (F1 intercross) progeny derived from the parental strains was developed to investigate the genetic basis of susceptibility to fibrosis. In the bleomycin studies the mice received 100 mg/kg (125 for females) of bleomycin, via mini osmotic pump. The animals were sacrificed at eight weeks following treatment or when their breathing rate indicated respiratory distress. In the radiation studies the mice were given a single dose of 14 or 16 Gy (Co$\sp{60})$ to the whole thorax and were sacrificed when moribund. The phenotype was defined as the percent of fibrosis area in the left lung as quantified with image analysis of histological sections. Quantitative trait loci (QTL) mapping was used to identify the chromosomal location of genes which contribute to susceptibility to bleomycin-induced pulmonary fibrosis in C57BL/6J mice compared to C3Hf/Kam mice and to determine if the QTL's which influence susceptibility to bleomycin-induced lung fibrosis in these progenitor strains could be implicated in susceptibility to radiation-induced lung fibrosis. For bleomycin, a genome wide scan revealed QTL's on chromosome 17, at the MHC, (LOD = 11.7 for males and 7.2 for females) accounting for approximately 21% of the phenotypic variance, and on chromosome 11 (LOD = 4.9), in male mice only, adding 8% of phenotypic variance. The bleomycin QTL on chromosome 17 was also implicated for susceptibility to radiation-induced fibrosis (LOD = 5.0) and contributes 7% of the phenotypic variance in the radiation study. In conclusion, susceptibility to both bleomycin-induced and radiation-induced pulmonary fibrosis are heritable traits, and are influenced by a genetic factor which maps to a genomic region containing the MHC. ^

Cost and power analysis of linkage disequilibrium methods for locating genes influencing quantitative traits in humans

Linkage disequilibrium methods can be used to find genes influencing quantitative trait variation in humans. Linkage disequilibrium methods can require smaller sample sizes than linkage equilibrium methods, such as the variance component approach to find loci with a specific effect size. The increase in power is at the expense of requiring more markers to be typed to scan the entire genome. This thesis compares different linkage disequilibrium methods to determine which factors influence the power to detect disequilibrium. The costs of disequilibrium and equilibrium tests were compared to determine whether the savings in phenotyping costs when using disequilibrium methods outweigh the additional genotyping costs.^ Nine linkage disequilibrium tests were examined by simulation. Five tests involve selecting isolated unrelated individuals while four involved the selection of parent child trios (TDT). All nine tests were found to be able to identify disequilibrium with the correct significance level in Hardy-Weinberg populations. Increasing linked genetic variance and trait allele frequency were found to increase the power to detect disequilibrium, while increasing the number of generations and distance between marker and trait loci decreased the power to detect disequilibrium. Discordant sampling was used for several of the tests. It was found that the more stringent the sampling, the greater the power to detect disequilibrium in a sample of given size. The power to detect disequilibrium was not affected by the presence of polygenic effects.^ When the trait locus had more than two trait alleles, the power of the tests maximized to less than one. For the simulation methods used here, when there were more than two-trait alleles there was a probability equal to 1-heterozygosity of the marker locus that both trait alleles were in disequilibrium with the same marker allele, resulting in the marker being uninformative for disequilibrium.^ The five tests using isolated unrelated individuals were found to have excess error rates when there was disequilibrium due to population admixture. Increased error rates also resulted from increased unlinked major gene effects, discordant trait allele frequency, and increased disequilibrium. Polygenic effects did not affect the error rates. The TDT, Transmission Disequilibrium Test, based tests were not liable to any increase in error rates.^ For all sample ascertainment costs, for recent mutations ($<$100 generations) linkage disequilibrium tests were less expensive than the variance component test to carry out. Candidate gene scans saved even more money. The use of recently admixed populations also decreased the cost of performing a linkage disequilibrium test. ^

AN INVESTIGATION INTO THE EFFECT OF TESTOSTERONE ON THE QUANTITATIVE MAINTENANCE OF SPERMATOGENESIS IN THE HYPOPHYSECTOMIZED ADULT RAT (PHARMACOKINETICS, MORPHOMETRY, STATISTICS)

Administration of gonadotropins or testosterone (T) will maintain qualitatively normal spermatogenesis and fertility in hypophysectomized (APX) rats. However, quantitative maintenance of the spermatogenic process in APX rats treated with T alone or in combination with follicle stimulating hormone (FSH) has not been demonstrated. Studies reported here were conducted to determine whether it would be possible to increase intratesticular testosterone (ITT) levels in APX rats to those found in normal animals by administration of appropriate amounts of testosterone propionate (TP) and if under these conditions spermatogenesis can be maintained quantitatively. Quantitative analysis of spermatogenesis was performed on stages VI and VII of the spermatogenic cycle utilizing criteria of Leblond and Clermont (1952) all cell types were enumerated. In a series of experiments designed to investigate the effects of T on spermatogenesis, TP was administered to 60 day old APX rats twice daily for 30 days in doses ranging from 0.6 to 15 mg/day or from 0.6 to 6.0 mg/day in combination with FSH. The results of this study demonstrate that the efficiency of transformation of type A to type B spermatogonia and the efficacy of the meiotic prophase are related to ITT levels, and that quantitatively normal completion of the reduction division requires normal ITT levels. The ratio of spermatids to spermatocytes in the vehicle-treated APX rats was 1:1.38; in the APX rats treated with 15 mg of TP it was 1:4.0 (the theoretically expected number). This study is probably the first to demonstrate: (1) the pharmacokinetics of TP, (2) the profile and quantity of T-immunoactivity in both serum and testicular tissue of APX and IC rats as well as APX rats treated with TP alone or in combination with FSH, (3) the direct correlation of serum T and ITT levels in treated APX rats (r = 0.9, p < 0.001) as well as in the IC rats (r = 0.9, p < 0.001), (4) the significant increase in the number of Type B spermatogonia, preleptotene and pachytene spermatocytes and round spermatids in TP-treated APX rats, (5) the correlation of the number of round spermatids formed in IC rats to ITT levels (r = 0.9, p < 0.001), and (6) the correlation of the quantitative maintenance of spermatogenesis with ITT levels (r = 0.7, p < 0.001) in the testes of TP-treated APX rats. These results provide direct experimental evidence for the key role of T in the spermatogenic process. ^

Adiposity and the perilipin gene: The role of single locus and multilocus variation, obesity-related quantitative traits, and disease-related phenotypes in the Atherosclerosis Risk in Communities (ARIC) study

Obesity is a complex multifactorial disease and is a public health priority. Perilipin coats the surface of lipid droplets in adipocytes and is believed to stabilize these lipid bodies by protecting triglyceride from early lipolysis. This research project evaluated the association between genetic variation within the human perilipin (PLIN) gene and obesity-related quantitative traits and disease-related phenotypes in Non-Hispanic White (NHW) and African American (AA) participants from the Atherosclerosis Risk in Communities (ARIC) Study. ^ Multivariate linear regression, multivariate logistic regression, and Cox proportional hazards models evaluated the association between single gene variants (rs2304794, rs894160, rs8179071, and rs2304795) and multilocus variation (rs894160 and rs2304795) within the PLIN gene and both obesity-related quantitative traits (body weight, body mass index [BMI], waist girth, waist-to-hip ratio [WHR], estimated percent body fat, and plasma total triglycerides) and disease-related phenotypes (prevalent obesity, metabolic syndrome [MetS], prevalent coronary heart disease [CHD], and incident CHD). Single variant analyses were stratified by race and gender within race while multilocus analyses were stratified by race. ^ Single variant analyses revealed that rs2304794 and rs894160 were significantly related to plasma triglyceride levels in all NHWs and NHW women. Among AA women, variant rs8179071 was associated with triglyceride levels and rs2304794 was associated with risk-raising waist circumference (>0.8 in women). The multilocus effects of variants rs894160 and rs2304795 were significantly associated with body weight, waist girth, WHR, estimated percent body fat, class II obesity (BMI ≥ 35 kg/m2), class III obesity (BMI ≥ 35 kg/m2), and risk-raising WHR (>0.9 in men and >0.8 in women) in AAs. Variant rs2304795 was significantly related to prevalent MetS among AA males and prevalent CHD in NHW women; multilocus effects of the PLIN gene were associated with prevalent CHD among NHWs. Rs2304794 was associated with incident CHD in the absence of the MetS among AAs. These findings support the hypothesis that variation within the PLIN gene influences obesity-related traits and disease-related phenotypes. ^ Understanding these effects of the PLIN genotype on the development of obesity can potentially lead to tailored health promotion interventions that are more effective. ^

An integrated molecular biology and computational approach to CYP1A1 expression in human brain

The cytochromes P450 comprise a superfamily of heme-containing mono-oxygenases. These enzymes metabolize numerous xenobiotics, but also play a role in metabolism of endogenous compounds. The P450 1A1 enzyme generally metabolizes polycyclic aromatic hydrocarbons, and its expression can be induced by aryl hydrocarbon receptor (AhR) activation. CYP1A1 is an exception to the generality that the majority of CYPs demonstrate highest expression in liver; CYP1Al is present in numerous extrahepatic tissues, including brain. This P450 has been observed in two forms, wildtype (WT) and brain variant (BV), arising from alternatively spliced mRNA transcripts. The CYP1A1 BV mRNA presented an exon deletion and was detected in human brain but not liver tissue of the same individuals. ^ Quantitative PCR analyses were performed to determine CYP1A1 WT and BV transcript expression levels in normal, bipolar disorder or schizophrenic groups. In our samples, we show that CYP1A1 BV mRNA, when present, is found alongside the full-length form. Furthermore, we demonstrate a significant decrease in expression of CYP1A1 in patients with bipolar disorder or schizophrenia. The expression level was not influenced by post-mortem interval, tissue pH, age, tobacco use, or lifetime antipsychotic medication load. ^ There is no indication of increased brain CYP1A1 expression in normal smokers versus non-smokers in these samples. We observed slightly increased CYP1A1 expression only in bipolar and schizophrenic smokers versus non-smokers. This may be indicative of complex interactions between neuronal chemical environments and AhR-mediated CYP1A1 induction in brain. ^ Structural homology modeling demonstrated that P450 1A1 BV has several alterations to positions/orientations of substrate recognition site residues compared to the WT isoform. Automated substrate docking was employed to investigate the potential binding of neurological signaling molecules and neurotropic drugs, as well as to differentiate specificities of the two P450 1A1 isoforms. We consistently observed that the BV isoform produced energetically favorable substrate dockings in orientations not observed for the same substrate in the WT isoform. These results demonstrated that structural differences, namely an expanded substrate access channel and active site, confer greater capacity for unique compound docking positions suggesting a metabolic profile distinct from the wildtype form for these test compounds. ^

Identification of genes associated with quantitative traits involved in cardiovascular disease and lipoprotein metabolism

Cardiovascular disease (CVD) is a threat to public health. It has been reported to be the leading cause of death in United States. The invention of next generation sequencing (NGS) technology has revolutionized the biomedical research. To investigate NGS data of CVD related quantitative traits would contribute to address the unknown etiology and disease mechanism of CVD. NHLBI's Exome Sequencing Project (ESP) contains CVD related phenotypes and their associated NGS exomes sequence data. Initially, a subset of next generation sequencing data consisting of 13 CVD-related quantitative traits was investigated. Only 6 traits, systolic blood pressure (SBP), diastolic blood pressure (DBP), height, platelet counts, waist circumference, and weight, were analyzed by functional linear model (FLM) and 7 currently existing methods. FLM outperformed all currently existing methods by identifying the highest number of significant genes and had identified 96, 139, 756, 1162, 1106, and 298 genes associated with SBP, DBP, Height, Platelet, Waist, and Weight respectively. ^

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.^