Clinical relevance of polymorphic DNA copy number variation (CNV) in African American women with stage I-II breast cancer

Objectives. The chief goal of this study was to analyze copy number variation (CNV) in breast cancer tumors from 25 African American women with early stage breast cancer (BC) using molecular inversion probes (MIP) in order to: (1) compare the degree of CNV in tumors compared to normal lymph nodes, and (2) determine whether gains and/or losses of genes in specific chromosomes differ between pathologic subtypes of breast cancer defined by known prognostic markers, (3) determine whether gains/losses in CN are associated with known oncogenes or tumor suppressor genes, and (4) determine whether increased gains/losses in CN for specific chromosomes were associated with differences in breast cancer recurrence. ^ Methods. Twenty to 37 nanograms of DNA extracted from 25 formalin-fixed paraffin embedded (FFPE) tumor samples and matched normal lymph nodes were added to individual tubes. Oligonucleotide probes with recognition sequences at each terminus were hybridized with a genomic target sequence to form a circular structure. Probes are released from genomic DNA obtained from FFPE samples, and those which have been correctly "circularized" in the proper allele/nucleotide reaction combination are amplified using polymerase chain reaction (PCR) primers. Amplicons were fluorescently labeled and the tag sequences released from the genome homology regions by treatment with uracil-N-glycosylase to cleave the probe at the site where uracils are present, and detected using a complementary tag array developed by Affymetrix. ^ Results. Analysis of CN gains and losses from tumors and normal tissues showed marked differences in tumors with numerous chromosomes affected. Similar changes were not observed in normal lymph nodes. When tumors were stratified into four groups based on expression or lack of expression of the estrogen receptor and HER2/neu, distinct patterns of CNV for different chromosomes were observed. Gains or losses in CN for specific chromosomes correlated with amplifications/deletions of particular oncogenes or tumor suppressor genes (i.e. such as found on chromosome 17) known to be associated with aggressive tumor phenotype and poor prognosis. There was a trend for increases in CN observed for chromosome 17 to correlate inversely with time to recurrence of BC (p=0.14 for trend). CNV was also observed for chromosomes 5, 8, 10, 11, and 16, which are known sites for several breast cancer susceptibility alleles. ^ Conclusions. This study is the first to validate the MIP technique, to correlate differences in gene expression with known prognostic tumor markers, and to correlate significant increases/decreases in CN with known tumor markers associated with prognosis. The results of this study may have far reaching public health implications towards identifying new high-risk groups based on genomic differences in CNP, both with respect to prognosis and response to therapy, and to eventually identify new therapeutic targets for prevention and treatment of this disease. ^

Genome-wide algorithm for detecting CNV associations with diseases

SNP genotyping arrays have been developed to characterize single-nucleotide polymorphisms (SNPs) and DNA copy number variations (CNVs). The quality of the inferences about copy number can be affected by many factors including batch effects, DNA sample preparation, signal processing, and analytical approach. Nonparametric and model-based statistical algorithms have been developed to detect CNVs from SNP genotyping data. However, these algorithms lack specificity to detect small CNVs due to the high false positive rate when calling CNVs based on the intensity values. Association tests based on detected CNVs therefore lack power even if the CNVs affecting disease risk are common. In this research, by combining an existing Hidden Markov Model (HMM) and the logistic regression model, a new genome-wide logistic regression algorithm was developed to detect CNV associations with diseases. We showed that the new algorithm is more sensitive and can be more powerful in detecting CNV associations with diseases than an existing popular algorithm, especially when the CNV association signal is weak and a limited number of SNPs are located in the CNV.^

KEAP1 E3 ligase-mediated downregulation of NF-kappaB signaling by targeting IKKbeta.

IkappaB kinase beta (IKKbeta) is involved in tumor development and progression through activation of the nuclear factor (NF)-kappaB pathway. However, the molecular mechanism that regulates IKKbeta degradation remains largely unknown. Here, we show that a Cullin 3 (CUL3)-based ubiquitin ligase, Kelch-like ECH-associated protein 1 (KEAP1), is responsible for IKKbeta ubiquitination. Depletion of KEAP1 led to the accumulation and stabilization of IKKbeta and to upregulation of NF-kappaB-derived tumor angiogenic factors. A systematic analysis of the CUL3, KEAP1, and RBX1 genomic loci revealed a high percentage of genome loss and missense mutations in human cancers that failed to facilitate IKKbeta degradation. Our results suggest that the dysregulation of KEAP1-mediated IKKbeta ubiquitination may contribute to tumorigenesis.

Role of obesity, weight gain, physical activity, and polymorphisms in the mTOR pathway and the risk of renal cell carcinoma

The interplay between obesity, physical activity, weight gain and genetic variants in mTOR pathway have not been studied in renal cell carcinoma (RCC). We examined the associations between obesity, weight gain, physical activity and RCC risk. We also analyzed whether genetic variants in the mTOR pathway could modify the association. Incident renal cell carcinoma cases and healthy controls were recruited from the University of Texas MD Anderson Cancer Center in Houston, Texas. Cases and controls were frequency-matched by age (±5 years), ethnicity, sex, and county of residence. Epidemiologic data were collected via in-person interview. A total of 577 cases and 593 healthy controls (all white) were included. One hundred ninety-two (192) SNPs from 22 genes were available and their genotyping data were extracted from previous genome-wide association studies. Logistic regression and regression spline were performed to obtain odds ratios. Obesity at age 20, 40, and 3 years prior to diagnosis/recruitment, and moderate and large weight gain from age 20 to 40 were each significantly associated with increased RCC risk. Low physical activity was associated with a 4.08-fold (95% CI: 2.92-5.70) increased risk. Five single nucleotide polymorphisms (SNPs) were significantly associated with RCC risk and their cumulative effect increased the risk by up to 72% (95% CI: 1.20-2.46). Strata specific effects for weight change and genotyping cumulative groups were observed. However, no interaction was suggested by our study. In conclusion, energy balance related risk factors and genetic variants in the mTOR pathway may jointly influence susceptibility to RCC. ^

INTEGRATIVE BIOMARKER IDENTIFICATION AND CLASSIFICATION USING HIGH THROUGHPUT ASSAYS

It is well accepted that tumorigenesis is a multi-step procedure involving aberrant functioning of genes regulating cell proliferation, differentiation, apoptosis, genome stability, angiogenesis and motility. To obtain a full understanding of tumorigenesis, it is necessary to collect information on all aspects of cell activity. Recent advances in high throughput technologies allow biologists to generate massive amounts of data, more than might have been imagined decades ago. These advances have made it possible to launch comprehensive projects such as (TCGA) and (ICGC) which systematically characterize the molecular fingerprints of cancer cells using gene expression, methylation, copy number, microRNA and SNP microarrays as well as next generation sequencing assays interrogating somatic mutation, insertion, deletion, translocation and structural rearrangements. Given the massive amount of data, a major challenge is to integrate information from multiple sources and formulate testable hypotheses. This thesis focuses on developing methodologies for integrative analyses of genomic assays profiled on the same set of samples. We have developed several novel methods for integrative biomarker identification and cancer classification. We introduce a regression-based approach to identify biomarkers predictive to therapy response or survival by integrating multiple assays including gene expression, methylation and copy number data through penalized regression. To identify key cancer-specific genes accounting for multiple mechanisms of regulation, we have developed the integIRTy software that provides robust and reliable inferences about gene alteration by automatically adjusting for sample heterogeneity as well as technical artifacts using Item Response Theory. To cope with the increasing need for accurate cancer diagnosis and individualized therapy, we have developed a robust and powerful algorithm called SIBER to systematically identify bimodally expressed genes using next generation RNAseq data. We have shown that prediction models built from these bimodal genes have the same accuracy as models built from all genes. Further, prediction models with dichotomized gene expression measurements based on their bimodal shapes still perform well. The effectiveness of outcome prediction using discretized signals paves the road for more accurate and interpretable cancer classification by integrating signals from multiple sources.

The Role of Acidic Fibroblast Growth Factor and Fibroblast Growth Factor Receptor 4 in High Grade Serous Carcinoma

Background: High grade serous carcinoma whether ovarian, tubal or primary peritoneal, continues to be the most lethal gynecologic malignancy in the USA. Although combination chemotherapy and aggressive surgical resection has improved survival in the past decade the majority of patients still succumb to chemo-resistant disease recurrence. It has recently been reported that amplification of 5q31-5q35.3 is associated with poor prognosis in patients with high grade serous ovarian carcinoma. Although the amplicon contains over 50 genes, it is notable for the presence of several members of the fibroblast growth factor signaling axis. In particular acidic fibroblast growth factor (FGF1) has been demonstrated to be one of the driving genes in mediating the observed prognostic effect of the amplicon in ovarian cancer patients. This study seeks to further validate the prognostic value of fibroblast growth receptor 4 (FGFR4), another candidate gene of the FGF/FGFR axis located in the same amplicon. The emphasis will be delineating the role the FGF1/FGFR4 signaling axis plays in high grade serous ovarian carcinoma; and test the feasibility of targeting the FGF1/FGFR4 axis therapeutically. Materials and Methods: Spearman and Pearson correlation studies on data generated from array CGH and transcriptome profiling analyses on 51 microdissected tumor samples were used to identify genes located on chromosome 5q31-35.3 that showed significant correlation between DNA and mRNA copy numbers. Significant correlation between FGF1 and FGFR4 DNA copy numbers was further validated by qPCR analysis on DNA isolated from 51 microdissected tumor samples. Immunolocalization and quantification of FGFR4 expression were performed on paraffin embedded tissue samples from 183 cases of high-grade serous ovarian carcinoma. The expression was then correlated with clinical data to assess impact on survival. The expression of FGF1 and FGFR4 in vitro was quantified by real-time PCR and western blotting in six high-grade serous ovarian carcinoma cell lines and compared to those in human ovarian surface epithelial cells to identify overexpression. The effect of FGF1 on these cell lines after serum starvation was quantified for in vitro cellular proliferation, migration/invasion, chemoresistance and survival utilizing a combination of commercially available colorimetric, fluorometric and electrical impedance assays. FGFR4 expression was then transiently silenced via siRNA transfection and the effects on response to FGF1, cellular proliferation, and migration were quantified. To identify relevant cellular pathways involved, responsive cell lines were transduced with different transcription response elements using the Cignal-Lenti reporter system and treated with FGF1 with and without transient FGFR4 knock down. This was followed by western blot confirmation for the relevant phosphoproteins. Anti-FGF1 antibodies and FGFR trap proteins were used to attempt inhibition of FGF mediated phenotypic changes and relevant signaling in vitro. Orthotopic intraperitoneal tumors were established in nude mice using serous cell lines that have been previously transfected with luciferase expressing constructs. The mice were then treated with FGFR trap protein. Tumor progression was then followed via bioluminescent imaging. The FGFR4 gene from 52 clinical samples was sequenced to screen for mutations. Results: FGFR4 DNA and mRNA copy numbers were significantly correlated and FGFR4 DNA copy number was significantly correlated with that of FGF1. Survival of patients with high FGFR4 expressing tumors was significantly shorter that those with low expression(median survival 28 vs 55 month p< 0.001) In a multivariate cox regression model FGFR expression significantly increased risk of death (HR 2.1, p<0.001). FGFR4 expression was significantly higher in all cell lines tested compared to HOSE, OVCA432 cell line in particular had very high expression suggesting amplification. FGF1 was also particularly overexpressed in OVCA432. FGF1 significantly increased cell survival after serum deprivation in all cell lines. Transient knock down of FGFR4 caused significant reduction in cell migration and proliferation in vitro and significantly decreased the proliferative effects of FGF1 in vitro. FGFR1, FGFR4 traps and anti-FGF1 antibodies did not show activity in vitro. OVCA432 transfected with the cignal lenti reporter system revealed significant activation of MAPK, NFkB and WNT pathways, western blotting confirmed the results. Reverse phase protein array (RPPA) analysis also showed activation of MAPK, AKT, WNT pathways and down regulation of E Cadherin. FGFR trap protein significantly reduced tumor growth in vivo in an orthotopic mouse model. Conclusions: Overexpression and amplification of several members of the FGF signaling axis present on the amplicon 5q31-35.3 is a negative prognostic indicator in high grade serous ovarian carcinoma and may drive poor survival associated with that amplicon. Activation of The FGF signaling pathway leads to downstream activation of MAPK, AKT, WNT and NFkB pathways leading to a more aggressive cancer phenotype with increased tumor growth, evasion of apoptosis and increased migration and invasion. Inhibition of FGF pathway in vivo via FGFR trap protein leads to significantly decreased tumor growth in an orthotopic mouse model.

Contribution of Ectodomain Mutations in Epidermal Growth Factor Receptor to Signaling in Glioblastoma Multiforme

CONTRIBUTION OF ECTODOMAIN MUTATIONS IN EPIDERMAL GROWTH FACTOR RECEPTOR TO SIGNALING IN GLIOBLASTOMA MULTIFORME Publication No._________ Marta Rojas, M.S. Supervisory Professor: Oliver Bögler, Ph.D. The Cancer Genome Atlas (TCGA) has conducted a comprehensive analysis of a large tumor cohort and has cataloged genetic alterations involving primary sequence variations and copy number aberrations of genes involved in key signaling pathways in glioblastoma (GBM). This dataset revealed missense ectodomain point mutations in epidermal growth factor receptor (EGFR), but the biological and clinical significance of these mutations is not well defined in the context of gliomas. In our study, we focused on understanding and defining the molecular mechanisms underlying the functions of EGFR ectodomain mutants. Using proteomic approaches to broadly analyze cell signaling, including antibody array and mass spectrometry-based methods, we found a differential spectrum of tyrosine phosphorylation across the EGFR ectodomain mutations that enabled us to stratify them into three main groups that correlate with either wild type EGFR (EGFR) or the long-studied mutant, EGFRvIII. Interestingly, one mutant shared characteristics of both groups suggesting a continuum of behaviors along which different mutants fall. Surprisingly, no substantial differences were seen in activation of classical downstream signaling pathways such as Akt and S6 pathways between these classes of mutants. Importantly, we demonstrated that ectodomain mutations lead to differential tumor growth capabilities in both in vitro (anchorage independent colony formation) and in vivo conditions (xenografts). Our data from the biological characterization allowed us to categorize the mutants into three main groups: the first group typified by EGFRvIII are mutations with a more aggressive phenotype including R108K and A289T; a second group characterized by a less aggressive phenotype exemplified by EGFR and the T263P mutation; and a third group which shared characteristics from both groups and is exemplified by the mutation A289D. In addition, we treated cells overexpressing the mutants with various agents employed in the clinic including temozolomide, cisplatin and tarceva. We found that cells overexpressing the mutants in general displayed resistance to the treatments. Our findings yield insights that help with the molecular characterization of these mutants. In addition, our results from the drug studies might be valuable in explaining differential responses to specific treatments in GBM patients.

Molecular alterations in nucleotide excision repair genes in malignant glioma and their relevance to malignant progression and drug resistance

Recently, it has become apparent that DNA repair mechanisms are involved in the malignant progression and resistance to therapy of gliomas. Many investigators have shown that increased levels of O6-methyl guanine DNA alkyltransferase, a DNA monoalkyl adduct repair enzyme, are correlated with resistance of malignant glioma cell lines to nitrosourea-based chemotherapy. Three important DNA excision repair genes ERCC1 (excision repair cross complementation group 1), ERCC2 (excision repair cross complementation group 2), and ERCC6 (excision repair cross complementation group 6) have been studied in human tumors. Gene copy number variation of ERCC1 and ERCC2 has been observed in primary glioma tissues. A number of reports describing a relationship between ERCC1 gene alterations and resistance to anti-cancer drugs have been also described. The levels of ERCC1 gene expression, however, have not been correlated with drug resistance in gliomas. The expression of ERCC6 gene transcribes has been shown to vary with tissue types and to be highest in the brain. There have been no comprehensive studies so far, however, of ERCC6 gene expression and molecular alterations in malignant glioma. This project examined the ERCC1 expression levels and correlated them with cisplatin resistance in malignant glioma cell lines. We also examined the molecular alterations of ERCC6 gene in primary glioma tissues and cells and analyzed whether these alterations are related to tumor progression and chemotherapy resistance. Our results indicate the presence of mutations and/or deletions in exons II and V of the ERCC6 gene, and these alterations are more frequent in exon II. Furthermore, the mutations and/or deletions in exon II were shown to be associated with increased malignant grade of gliomas. The results on the Levels of ERCC1 gene transcripts showed that expression levels correlate with cisplatin resistance. The increase in ERCC1 mRNA induced by cisplatin could be down-regulated by cyclosporin A and herbimycin A. The results of this study are likely to provide useful information for clinical treatment of human gliomas. ^

RANGE ADAPTIVE PROTON THERAPY FOR PROSTATE CANCER

Purpose: The rapid distal falloff of a proton beam allows for sparing of normal tissues distal to the target. However proton beams that aim directly towards critical structures are avoided due to concerns of range uncertainties, such as CT number conversion and anatomy variations. We propose to eliminate range uncertainty and enable prostate treatment with a single anterior beam by detecting the proton’s range at the prostate-rectal interface and adaptively adjusting the range in vivo and in real-time. Materials and Methods: A prototype device, consisting of an endorectal liquid scintillation detector and dual-inverted Lucite wedges for range compensation, was designed to test the feasibility and accuracy of the technique. Liquid scintillation filled volume was fitted with optical fiber and placed inside the rectum of an anthropomorphic pelvic phantom. Photodiode-generated current signal was generated as a function of proton beam distal depth, and the spatial resolution of this technique was calculated by relating the variance in detecting proton spills to its maximum penetration depth. The relative water-equivalent thickness of the wedges was measured in a water phantom and prospectively tested to determine the accuracy of range corrections. Treatment simulation studies were performed to test the potential dosimetric benefit in sparing the rectum. Results: The spatial resolution of the detector in phantom measurement was 0.5 mm. The precision of the range correction was 0.04 mm. The residual margin to ensure CTV coverage was 1.1 mm. The composite distal margin for 95% treatment confidence was 2.4 mm. Planning studies based on a previously estimated 2mm margin (90% treatment confidence) for 27 patients showed a rectal sparing up to 51% at 70 Gy and 57% at 40 Gy relative to IMRT and bilateral proton treatment. Conclusion: We demonstrated the feasibility of our design. Use of this technique allows for proton treatment using a single anterior beam, significantly reducing the rectal dose.

Rtr1 is the Saccharomyces cerevisiae homolog of a novel family of RNA polymerase II-binding proteins.

Cells must rapidly sense and respond to a wide variety of potentially cytotoxic external stressors to survive in a constantly changing environment. In a search for novel genes required for stress tolerance in Saccharomyces cerevisiae, we identified the uncharacterized open reading frame YER139C as a gene required for growth at 37 degrees C in the presence of the heat shock mimetic formamide. YER139C encodes the closest yeast homolog of the human RPAP2 protein, recently identified as a novel RNA polymerase II (RNAPII)-associated factor. Multiple lines of evidence support a role for this gene family in transcription, prompting us to rename YER139C RTR1 (regulator of transcription). The core RNAPII subunits RPB5, RPB7, and RPB9 were isolated as potent high-copy-number suppressors of the rtr1Delta temperature-sensitive growth phenotype, and deletion of the nonessential subunits RPB4 and RPB9 hypersensitized cells to RTR1 overexpression. Disruption of RTR1 resulted in mycophenolic acid sensitivity and synthetic genetic interactions with a number of genes involved in multiple phases of transcription. Consistently, rtr1Delta cells are defective in inducible transcription from the GAL1 promoter. Rtr1 constitutively shuttles between the cytoplasm and nucleus, where it physically associates with an active RNAPII transcriptional complex. Taken together, our data reveal a role for members of the RTR1/RPAP2 family as regulators of core RNAPII function.

Interlinked dual-time feedback loops can enhance robustness to stochasticity and persistence of memory.

Multiple interlinked positive feedback loops shape the stimulus responses of various biochemical systems, such as the cell cycle or intracellular Ca2+ release. Recent studies with simplified models have identified two advantages of coupling fast and slow feedback loops. This dual-time structure enables a fast response while enhancing resistances of responses and bistability to stimulus noise. We now find that (1) the dual-time structure similarly confers resistance to internal noise due to molecule number fluctuations, and (2) model variants with altered coupling, which better represent some specific biochemical systems, share all the above advantages. We also develop a similar bistable model with coupling of a fast autoactivation loop to a slow loop. This model's topology was suggested by positive feedback proposed to play a role in long-term synaptic potentiation (LTP). The advantages of fast response and noise resistance are also present in this autoactivation model. Empirically, LTP develops resistance to reversal over approximately 1h . The model suggests this resistance may result from increased amounts of synaptic kinases involved in positive feedback.

Gene by BMI Interactions Influencing C-Reactive Protein Levels in European-Americans

C-Reactive Protein (CRP) is a biomarker indicating tissue damage, inflammation, and infection. High-sensitivity CRP (hsCRP) is an emerging biomarker often used to estimate an individual’s risk for future coronary heart disease (CHD). hsCRP levels falling below 1.00 mg/l indicate a low risk for developing CHD, levels ranging between 1.00 mg/l and 3.00 mg/l indicate an elevated risk, and levels exceeding 3.00 mg/l indicate high risk. Multiple Genome-Wide Association Studies (GWAS) have identified a number of genetic polymorphisms which influence CRP levels. SNPs implicated in such studies have been found in or near genes of interest including: CRP, APOE, APOC, IL-6, HNF1A, LEPR, and GCKR. A strong positive correlation has also been found to exist between CRP levels and BMI, a known risk factor for CHD and a state of chronic inflammation. We conducted a series of analyses designed to identify loci which interact with BMI to influence CRP levels in a subsample of European-Americans in the ARIC cohort. In a stratified GWA analysis, 15 genetic regions were identified as having significantly (p-value < 2.00*10-3) distinct effects on hsCRP levels between the two obesity strata: lean (18.50 kg/m2 < BMI < 24.99 kg/m2) and obese (BMI ≥ 30.00 kg/m2). A GWA analysis performed on all individuals combined (i.e. not a priori stratified for obesity status) with the inclusion of an additional parameter for BMI by gene interaction, identified 11 regions which interact with BMI to influence hsCRP levels. Two regions containing the genes GJA5 and GJA8 (on chromosome 1) and FBXO11 (on chromosome 2) were identified in both methods of analysis suggesting that these genes possibly interact with BMI to influence hsCRP levels. We speculate that atrial fibrillation (AF), age-related cataracts and the TGF-β pathway may be the biological processes influenced by the interaction of GJA5, GJA8 and FBXO11, respectively, with BMI to cause changes in hsCRP levels. Future studies should focus on the influence of gene x bmi interaction on AF, age-related cataracts and TGF-β.

DUAL ENERGY COMPUTED TOMOGRAPHY FOR PROTON THERAPY TREATMENT PLANNING

Proton radiation therapy is gaining popularity because of the unique characteristics of its dose distribution, e.g., high dose-gradient at the distal end of the percentage-depth-dose curve (known as the Bragg peak). The high dose-gradient offers the possibility of delivering high dose to the target while still sparing critical organs distal to the target. However, the high dose-gradient is a double-edged sword: a small shift of the highly conformal high-dose area can cause the target to be substantially under-dosed or the critical organs to be substantially over-dosed. Because of that, large margins are required in treatment planning to ensure adequate dose coverage of the target, which prevents us from realizing the full potential of proton beams. Therefore, it is critical to reduce uncertainties in the proton radiation therapy. One major uncertainty in a proton treatment is the range uncertainty related to the estimation of proton stopping power ratio (SPR) distribution inside a patient. The SPR distribution inside a patient is required to account for tissue heterogeneities when calculating dose distribution inside the patient. In current clinical practice, the SPR distribution inside a patient is estimated from the patient’s treatment planning computed tomography (CT) images based on the CT number-to-SPR calibration curve. The SPR derived from a single CT number carries large uncertainties in the presence of human tissue composition variations, which is the major drawback of the current SPR estimation method. We propose to solve this problem by using dual energy CT (DECT) and hypothesize that the range uncertainty can be reduced by a factor of two from currently used value of 3.5%. A MATLAB program was developed to calculate the electron density ratio (EDR) and effective atomic number (EAN) from two CT measurements of the same object. An empirical relationship was discovered between mean excitation energies and EANs existing in human body tissues. With the MATLAB program and the empirical relationship, a DECT-based method was successfully developed to derive SPRs for human body tissues (the DECT method). The DECT method is more robust against the uncertainties in human tissues compositions than the current single-CT-based method, because the DECT method incorporated both density and elemental composition information in the SPR estimation. Furthermore, we studied practical limitations of the DECT method. We found that the accuracy of the DECT method using conventional kV-kV x-ray pair is susceptible to CT number variations, which compromises the theoretical advantage of the DECT method. Our solution to this problem is to use a different x-ray pair for the DECT. The accuracy of the DECT method using different combinations of x-ray energies, i.e., the kV-kV, kV-MV and MV-MV pair, was compared using the measured imaging uncertainties for each case. The kV-MV DECT was found to be the most robust against CT number variations. In addition, we studied how uncertainties propagate through the DECT calculation, and found general principles of selecting x-ray pairs for the DECT method to minimize its sensitivity to CT number variations. The uncertainties in SPRs estimated using the kV-MV DECT were analyzed further and compared to those using the stoichiometric method. The uncertainties in SPR estimation can be divided into five categories according to their origins: the inherent uncertainty, the DECT modeling uncertainty, the CT imaging uncertainty, the uncertainty in the mean excitation energy, and SPR variation with proton energy. Additionally, human body tissues were divided into three tissue groups – low density (lung) tissues, soft tissues and bone tissues. The uncertainties were estimated separately because their uncertainties were different under each condition. An estimate of the composite range uncertainty (2s) was determined for three tumor sites – prostate, lung, and head-and-neck, by combining the uncertainty estimates of all three tissue groups, weighted by their proportions along typical beam path for each treatment site. In conclusion, the DECT method holds theoretical advantages in estimating SPRs for human tissues over the current single-CT-based method. Using existing imaging techniques, the kV-MV DECT approach was capable of reducing the range uncertainty from the currently used value of 3.5% to 1.9%-2.3%, but it is short to reach our original goal of reducing the range uncertainty by a factor of two. The dominant source of uncertainties in the kV-MV DECT was the uncertainties in CT imaging, especially in MV CT imaging. Further reduction in beam hardening effect, the impact of scatter, out-of-field object etc. would reduce the Hounsfeld Unit variations in CT imaging. The kV-MV DECT still has the potential to reduce the range uncertainty further.

CHROMOSOMAL MAPPING, LINKAGE RELATIONSHIPS, AND EVOLUTIONARY STUDIES OF THE HUMAN ANONYMOUS DNA CLONE D1S1 (IN SITU HYBRIDIZATION, PRIMATES, GENE MAPPING, AUTOSOMAL DOMINANT RETINITIS PIGMENTOSA, GENE DUPLICATION)

D1S1, an anonymous human DNA clone originally called (lamda)Ch4-H3 or (lamda)H3, was the first single copy mapped to a human chromosome (1p36) by in situ hybridization. The chromosomal assignment has been confirmed in other laboratories by repeating the in situ hybridization but not by another method. In the present study, hybridization to a panel of hamster-human somatic cell hybrids revealed copies of D1S1 on both chromosomes 1 and 3. Subcloning D1S1 showed that the D1S1 clone itself is from chromosome 3, and the sequence detected by in situ hybridization is at least two copies of part of the chromosome 3 copy. This finding demonstrates the importance of verifying gene mapping with two methods and questions the accuracy of in situ hybridization mapping.^ Non-human mammals have only one copy of D1S1, and the non-human primate D1S1 map closely resembles the human chromosome 3 copy. Thus, the human chromosome 1 copies appear to be part of a very recent duplication that occurred after the divergence between humans and the other great apes.^ A moderately informative HindIII D1S1 RFLP was mapped to chromosome 3. This marker and 12 protein markers were applied to a linkage study of autosomal dominant retinitis pigmentosa (ADRP). None of the markers proved linkage, but adding the three families examined to previously published data raises the ADRP:Rh lod score to 1.92 at (THETA) = 0.30. ^

EXPRESSION OF THE HUMAN PLACENTAL LACTOGEN GENES

The placenta is the site of synthesis of various peptide and steroid hormones related to pregnancy. Human placental lactogen (hPL) is the predominant peptide hormone secreted by term placenta and its synthesis is tissue-specific and coupled to placenta development. The objective of this work was to study the structure and expression of the hPL.^ Poly(A('+))RNA from human term placenta was translated in a mouse-derived cell-free system. A major band corresponding to pre-hPL and a minor band comigrating with mature hPL, represent (TURN)15% of the total radioactively labeled proteins. Analysis of the poly(A('+))RNA showed a prominent band at approximately 860 nucleotides. A corresponding band was observed in Northern blots of total RNA, hybridized with {('32)P}-labeled recombinant plasmid containing a portion of hPL cDNA. Similar analyses of nuclear RNA showed at least four additional bands at 990, 1200, 1460 and 1760 nucleotides, respectively, which are likely precursors of hPL mRNA. Poly(A('+))RNA was used to construct a cDNA library, of which approximately 5% of the clones were found to hybridize to hPL DNA sequences. Heteroduplexes constructed between a clone containing a 815 bp hPL cDNA insert and a hPL genomic DNA clone revealed four small intervening sequences which can account for the lengths observed in hnRNA molecules.^ Recombinant plasmid HCS-pBR322 containing a 550 bp insert of a cDNA transcript of human placental lactogen (hPL) mRNA was ('3)H-labeled an hybridized in situ to human chromosome preparations. These experiments allowed assignment of the hPL and growth hormone (hGH) genes, which have over 90% nucleotide homology in their coding sequences, to band q22-24 of chromosome 17. A gene copy number experiment showed that both genes are present in (TURN)3 copies per haploid genome.^ Experiments were designed to determine if all members of the hPL gene cluster, consisting of four non-allelic genes, are transcribed in term placenta. Advantage was taken of differences in restriction endonuclease sites in the coding portions of the different hPL genes, to distinguish the putative cDNAs of the transcriptionally active genes. Two genes were found to be represented in the cDNA library and their cDNA transcripts were isolated and characterized. Three independent methods showed that their corresponding mRNAs are about equally represented in the hPL mRNA population. The two cDNAs code for prehPL proteins which differ at a single amino acid position. However the secreted hPLs have identical amino acid sequences. A tetramer insertion duplication was found in a palindrome area of the 3' untranslated region of one of the hPL mRNAs. ^