Functional analysis of the sea urchin {\it orthodenticle\/} -related gene, {\it SPOTX\/}, in aboral ectoderm-specific gene activation and aboral ectoderm cell fate specification

An important question in developmental biology is how embryonic cell types are derived from a fertilized egg. To address this question, this thesis investigates the mechanisms by which the aboral ectoderm-specific Spec2a gene is spatially and temporally regulated during sea urchin embryogenesis. The Spec2a gene of the sea urchin Strongylocentratus purpuratus has served as a valuable maker to understand the basis of lineage-specific gene activation and the role of transcription factors in cell fate specification. The hypothesis is that transcription factors responsible for cell type-specific gene activation are key components in the initial cell specification step. The Spec2a gene, which encodes a small cytosolic calcium-binding protein, is expressed exclusively in aboral ectoderm cell lineages. The 1516-bp control region of the Spec2a gene contains a 188-bp enhancer element required for temporal activation and aboral ectoderm/mesenchyme cell expression, while an unidentified element upstream of the enhancer represses expression in mesenchyme cells. Using an enhancer activation assay, combined with site-directed mutagenesis, I showed that three TAATCC/T sites within the enhancer are responsible for enhancer activity. Mutagenizing these sites and a fourth one just upstream abolished all activity from the Spec2a control region. A 77-bp DNA fragment from the Spec2a enhancer containing two of the TAATCC/T sites is sufficient for aboral ectoderm/mesenchyme cell expression. A cDNA encoding SpOtx, an orthodenticle-related protein, was cloned from S. purpuratus and shown to bind with high affinity to the TAATCC/T sequences within the Spec2a control region. SpOtx transcripts were found initially in all cells of the cleaving embryo, but they gradually became restricted to oral ectoderm and endoderm cells, suggesting that SpOtx might play a role in the initial temporal activation of the Spec2a gene and most likely has additional functions in the developing embryo. To reveal the broader biological functions of SpOtx, I injected SpOtx mRNA into living sea urchin eggs to determine what effects overexpressing the SpOtx protein might have on embryo development. SpOtx mRNA-injected embryos displayed dramatic alterations in development. Instead of developing into pluteus larvae with 15 different cell types, uniform epithelia balls were formed. These balls consisted of a thin layer of squamous cells with short cilia highly reminiscent of aboral ectoderm. Immunohistochemical staining and RT-PCR demonstrated that the SpOtx-injected embryoids expressed aboral ectoderm markers uniformly, but showed very weak or no expression of markers for non-aboral ectoderm cell types. These data strongly suggested that overexpression of SpOtx redirected the normal fate of non-aboral ectoderm cells to that of aboral ectoderm. These results show that SpOtx is involved in aboral ectoderm differentiation by activating aboral ectoderm-specific genes and that modulating its expression can lead to changes in cell fate. ^

Analysis of differential gene expression in nontumorigenic glioblastoma multiforme

Loss of chromosome 10 represents the most common cytogenetic abnormality in high grade gliomas (glioblastoma multiforme). To identify genes involved in the malignant progression of human gliomas, a subtractive hybridization was performed between a tumorigenic glioblastoma cell line (LG11) and a nontumorgenic hybrid cell (LG11.3) containing an introduced chromosome 10. LG11 mRNA was subtracted from LG11.3 cDNA to produce cDNA probes enriched for sequences whose expression differs quantitatively from the parental tumorigenic cells. Both known and novel sequences were identified as a result of the subtraction. Northern blot analysis was then used to confirm differential expression of several subtracted clones. One novel clone, clone 17, identified a 2.6 kb message that showed a consistent two to four fold increase in expression in the LG11.3 nontumorigenic cells. Clone 17 (340 bp) was used successfully to screen for a near full-length version, RIG (regulated in glioma), which was 2,569 bp in size. The RIG cDNA sequence showed homology to clone 17 and to an anonymous EST (IB666), but to no previously identified genes. This screening effort also identified several independent clones representing novel sequences, most of which failed to show increased expression in the nontumorigenic GBM cells. Tissue distribution studies of RIG indicated highest levels of expression in human brain with appreciably lower levels in heart and lung. In vitro transcription and translation experiments demonstrated the ability of RIG to direct the synthesis of a 13 kD protein product. However, open reading frame analysis revealed no identify with previously described motifs or any known proteins. Using a combination of somatic cell hybrid panels and in situ hybridization, the RIG gene was mapped to chromosome 11p14-11p15. Further study of RIG and related gene products may provide insight into the negative regulation of glial oncogenesis. ^

Molecular studies of human high grade B cell non -Hodgkin's lymphomas

The non-Hodgkin's B cell lymphomas are a diverse group of neoplastic diseases. The incidence rate of the malignant tumors has been rising rapidly over the past twenty years in the United States and worldwide. The lack of insight to pathogenesis of the disease poses a significant problem in the early detection and effective treatment of the human malignancies. These studies attempted to investigate the molecular basis of pathogenesis of the human high grade B cell non-Hodgkin's lymphomas with a reverse genetic approach. The specific objective was to clone gene(s) which may play roles in development and progression of human high grade B cell non-Hodgkin's lymphomas.^ The messenger RNAs from two high grade B cell lymphoma lines, CJ and RR, were used for construction of cDNA libraries. Differential screening of the derived cDNA libraries yielded a 1.4 kb cDNA clone. The gene, designated as NHL-B1.4, was shown to be highly amplified and over-expressed in the high grade B cell lymphoma lines. It was not expressed in the peripheral blood lymphoid cells from normal donors. However, it was inducible in peripheral blood T lymphocytes by a T cell mitogen, PHA, but could not be activated in normal B cells by B cell mitogen PMA. Further molecular characterization revealed that the gene may have been rearranged in the RR and some other B cell lymphoma lines. The coding capacity of the cDNA has been confirmed by a rabbit reticulocyte lysate and wheat germ protein synthesis system. A recombinant protein with a molecular weight of approximate 30 kDa was visualized in autoradiogram. Polyclonal antisera have been generated by immunization of two rabbits with the NHL-B1.4 recombinant protein produced in the E. coli JM109. The derived antibody can recognize a natural protein with molecular weight of 49 kDa in cell lysate of activated peripheral T lymphocytes of normal donors and both the cell lysate and supernatant of RR B cell lymphoma lines. The possible biologic functions of the molecule has been tested preliminarily in a B lymphocyte proliferation assay. It was found that the Q-sepharose chromatograph purified supernatant of COS cell transfection could increase tritiated thymidine uptake by B lymphocytes but not by T lymphocytes. The B cell stimulatory activity of the supernatant of COS cell tranfection could be neutralized by the polyclonal antisera, indicating that the NHL-B1.4 gene product may be a molecule with BCGF-like activity.^ The expression profiles of NHL-B1.4 in normal and neoplastic lymphoid cells were consistent with the current B lymphocyte activation model and autocrine hypothesis of high grade B cell lymphomagenesis. These results suggested that the NHL-B1.4 cDNA may be a disease-related gene of human high grade B cell lymphomas, which may codes for a postulated B cell autocrine growth factor. ^

A urokinase receptor promoter element (-152/-135) bound with an AP-2-alpha-related protein and SP1 is required for the induction of gene expression by PMA and SRC

The urokinase-type plasminogen activator receptor (u-PAR) promotes extracellular matrix degradation, invasion and metastasis. A first objective of this dissertation was to identify cis-elements and trans-acting factors activating u-PAR gene expression through a previously footprinted (–148/–124) promoter region. Mobility shifting experiments on nuclear extracts of a high u-PAR-expressing colon cancer cell line (RKO) indicated Sp1, Sp3 and a factor similar to, but distinct from, AP-2α bound to an oligonucleotide spanning –152/–135. Mutations preventing the binding of the AP-2α-related factor reduced u-PAR promoter activity. In RKO, the expression of a dominant negative AP-2 (AP-2αB) diminished u-PAR promoter activity, protein and u-PAR mediated laminin degradation. Conversely, u-PAR promoter activity in low u-PAR-expressing GEO cells was increased by AP-2αA expression. PMA treatment, which induces u-PAR expression, caused an increased amount of the AP-2α-related factor-containing complex in GEO, and mutations preventing AP-2α-like and Sp1/Sp3 binding reduced the u-PAR promoter stimulation by PMA. In resected colon cancers, u-PAR protein amounts were related to the amount of the AP-2α-related factor-containing complex. In conclusion, constitutive and PMA- inducible u-PAR gene expression and -proteolysis are mediated partly through transactivation via a promoter sequence (–152/435) bound with an AP-2α-related factor and Sp1/Sp3. ^ A second interest of this dissertation was to determine if a constitutively active Src regulates the transcription of the u-PAR gene, since c-src expression increases invasion in colon cancer. Increased u-PAR protein and laminin degradation paralleling elevated Src activity was evident in SW480 colon cancer cells stably expressing a constitutively active Src (Y- c-src527F). Nuclear run-on experiments indicated that this was due largely to transcriptional activation. While transient transfection of SW480 cells with Y-c-src527F induced a u-PAR-CAT-reporter, mutations preventing Sp1-binding to promoter region –152/435 abolished this induction. Mobility shift assays revealed increased Sp1 binding to region –152/135 with nuclear extracts of Src-transfected SW480 cells. Finally, the amounts of endogenous u-PAR in resected colon cancers significantly correlated with Src-activity. These data suggest that u-PAR gene expression and proteolysis are regulated by Src, this requiring the promoter region (–152/–135) bound with Sp1, thus, demonstrating for the first time that transcription factor Sp1 is a downstream effector of Src. ^

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

GENETIC ANALYSIS OF THE FUNCTION OF THE DROSOPHILA DOUBLESEX-RELATED FACTOR dmrt93B

DMRT (Doublesex and Mab-3 related transcription factor) proteins generally associated with sexual differentiation in many organisms share a common DNA binding domain and are often expressed in reproductive tissues. Aside from doublesex, which is a central factor in the regulation of sex determination, Drosophila possesses three different dmrt genes that are of unknown function. Because the association with sexual differentiation and reproduction is not universal and some DMRT proteins have been found to play other developmental roles we chose to further characterize one of these Drosophila genes. We carried out genetic analysis of dmrt93B, which was previously found to be expressed sex-specifically in the developing somatic gonad and to affect testis morphogenesis in RNAi knockdowns. In order to disrupt this gene, the GAL4 yeast transcriptional activator followed by a polyadenylation signal was inserted after the dmrt93B start codon and introduced into the genome by homologous recombination. Analysis of the knock-in mutation as well as a small deletion removing all dmrt93B sequence demonstrate that loss of function causes partial lethality at the late pupal stage. Surprisingly, these mutations have no significant effect on gonad formation or male fertility. Analysis of GAL4-driven GFP reporter expression indicates that the dmrt93B promoter activity is highly specific to neurons in the suboesophageal and proventricular ganglion in larva and adult of both sexes suggesting a possible role in digestive tract function. Using the Capillary Feeder (CAFÉ) assay to measure daily food intake we find that reduction in this gene’s function leads to an increase in food consumption. These results suggest dmrt93 plays an important role in the formation or maintenance of neurons that affect feeding and support the idea that dmrt genes may not be restricted to roles in sexual differentiation.

A NEW TUMOR SUPPRESSOR GENE CANDIDATE REGULATED BY THE NON-CODING RNA PCA3 IN HUMAN PROSTATE CANCER

Prostate cancer is the second leading cause of cancer-related death and the most common non-skin cancer in men in the USA. Considerable advancements in the practice of medicine have allowed a significant improvement in the diagnosis and treatment of this disease and, in recent years, both incidence and mortality rates have been slightly declining. However, it is still estimated that 1 man in 6 will be diagnosed with prostate cancer during his lifetime, and 1 man in 35 will die of the disease. In order to identify novel strategies and effective therapeutic approaches in the fight against prostate cancer, it is imperative to improve our understanding of its complex biology since many aspects of prostate cancer initiation and progression still remain elusive. The study of tumor biomarkers, due to their specific altered expression in tumor versus normal tissue, is a valid tool for elucidating key aspects of cancer biology, and may provide important insights into the molecular mechanisms underlining the tumorigenesis process of prostate cancer. PCA3, is considered the most specific prostate cancer biomarker, however its biological role, until now, remained unknown. PCA3 is a long non-coding RNA (ncRNA) expressed from chromosome 9q21 and its study led us to the discovery of a novel human gene, PC-TSGC, transcribed from the opposite strand and in an antisense orientation to PCA3. With the work presented in this thesis, we demonstrate that PCA3 exerts a negative regulatory role over PC-TSGC, and we propose PC-TSGC to be a new tumor suppressor gene that contrasts the transformation of prostate cells by inhibiting Rho-GTPases signaling pathways. Our findings provide a biological role for PCA3 in prostate cancer and suggest a new mechanism of tumor suppressor gene inactivation mediated by non-coding RNA. Also, the characterization of PCA3 and PC-TSGC led us to propose a new molecular pathway involving both genes in the transformation process of the prostate, thus providing a new piece of the jigsaw puzzle representing the complex biology of prostate cancer.

Multifactor effects and evidence of potential interaction between complement factor H Y402H and LOC387715 A69S in age-related macular degeneration.

BACKGROUND: Variants in the complement cascade genes and the LOC387715/HTRA1, have been widely reported to associate with age-related macular degeneration (AMD), the most common cause of visual impairment in industrialized countries. METHODS/PRINCIPAL FINDINGS: We investigated the association between the LOC387715 A69S and complement component C3 R102G risk alleles in the Finnish case-control material and found a significant association with both variants (OR 2.98, p = 3.75 x 10(-9); non-AMD controls and OR 2.79, p = 2.78 x 10(-19), blood donor controls and OR 1.83, p = 0.008; non-AMD controls and OR 1.39, p = 0.039; blood donor controls), respectively. Previously, we have shown a strong association between complement factor H (CFH) Y402H and AMD in the Finnish population. A carrier of at least one risk allele in each of the three susceptibility loci (LOC387715, C3, CFH) had an 18-fold risk of AMD when compared to a non-carrier homozygote in all three loci. A tentative gene-gene interaction between the two major AMD-associated loci, LOC387715 and CFH, was found in this study using a multiplicative (logistic regression) model, a synergy index (departure-from-additivity model) and the mutual information method (MI), suggesting that a common causative pathway may exist for these genes. Smoking (ever vs. never) exerted an extra risk for AMD, but somewhat surprisingly, only in connection with other factors such as sex and the C3 genotype. Population attributable risks (PAR) for the CFH, LOC387715 and C3 variants were 58.2%, 51.4% and 5.8%, respectively, the summary PAR for the three variants being 65.4%. CONCLUSIONS/SIGNIFICANCE: Evidence for gene-gene interaction between two major AMD associated loci CFH and LOC387715 was obtained using three methods, logistic regression, a synergy index and the mutual information (MI) index.

Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies.

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

SURVIVAL PREDICTION FOR BRAIN TUMOR PATIENTS USING GENE EXPRESSION DATA

Brain tumor is one of the most aggressive types of cancer in humans, with an estimated median survival time of 12 months and only 4% of the patients surviving more than 5 years after disease diagnosis. Until recently, brain tumor prognosis has been based only on clinical information such as tumor grade and patient age, but there are reports indicating that molecular profiling of gliomas can reveal subgroups of patients with distinct survival rates. We hypothesize that coupling molecular profiling of brain tumors with clinical information might improve predictions of patient survival time and, consequently, better guide future treatment decisions. In order to evaluate this hypothesis, the general goal of this research is to build models for survival prediction of glioma patients using DNA molecular profiles (U133 Affymetrix gene expression microarrays) along with clinical information. First, a predictive Random Forest model is built for binary outcomes (i.e. short vs. long-term survival) and a small subset of genes whose expression values can be used to predict survival time is selected. Following, a new statistical methodology is developed for predicting time-to-death outcomes using Bayesian ensemble trees. Due to a large heterogeneity observed within prognostic classes obtained by the Random Forest model, prediction can be improved by relating time-to-death with gene expression profile directly. We propose a Bayesian ensemble model for survival prediction which is appropriate for high-dimensional data such as gene expression data. Our approach is based on the ensemble "sum-of-trees" model which is flexible to incorporate additive and interaction effects between genes. We specify a fully Bayesian hierarchical approach and illustrate our methodology for the CPH, Weibull, and AFT survival models. We overcome the lack of conjugacy using a latent variable formulation to model the covariate effects which decreases computation time for model fitting. Also, our proposed models provides a model-free way to select important predictive prognostic markers based on controlling false discovery rates. We compare the performance of our methods with baseline reference survival methods and apply our methodology to an unpublished data set of brain tumor survival times and gene expression data, selecting genes potentially related to the development of the disease under study. A closing discussion compares results obtained by Random Forest and Bayesian ensemble methods under the biological/clinical perspectives and highlights the statistical advantages and disadvantages of the new methodology in the context of DNA microarray data analysis.

Repression of tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) but not its receptors during oral cancer progression.

BACKGROUND: TRAIL plays an important role in host immunosurveillance against tumor progression, as it induces apoptosis of tumor cells but not normal cells, and thus has great therapeutic potential for cancer treatment. TRAIL binds to two cell-death-inducing (DR4 and DR5) and two decoy (DcR1, and DcR2) receptors. Here, we compare the expression levels of TRAIL and its receptors in normal oral mucosa (NOM), oral premalignancies (OPM), and primary and metastatic oral squamous cell carcinomas (OSCC) in order to characterize the changes in their expression patterns during OSCC initiation and progression. METHODS: DNA microarray, immunoblotting and immunohistochemical analyses were used to examine the expression levels of TRAIL and its receptors in oral epithelial cell lines and in archival tissues of NOM, OPM, primary and metastatic OSCC. Apoptotic rates of tumor cells and tumor-infiltrating lymphocytes (TIL) in OSCC specimens were determined by cleaved caspase 3 immunohistochemistry. RESULTS: Normal oral epithelia constitutively expressed TRAIL, but expression was progressively lost in OPM and OSCC. Reduction in DcR2 expression levels was noted frequently in OPM and OSCC compared to respective patient-matched uninvolved oral mucosa. OSCC frequently expressed DR4, DR5 and DcR1 but less frequently DcR2. Expression levels of DR4, DR5 and DcR1 receptors were not significantly altered in OPM, primary OSCC and metastatic OSCC compared to patient-matched normal oral mucosa. Expression of proapoptotic TRAIL-receptors DR4 and DR5 in OSCC seemed to depend, at least in part, on whether or not these receptors were expressed in their parental oral epithelia. High DR5 expression in primary OSCC correlated significantly with larger tumor size. There was no significant association between TRAIL-R expression and OSSC histology grade, nodal status or apoptosis rates of tumor cells and TIL. CONCLUSION: Loss of TRAIL expression is an early event during oral carcinogenesis and may be involved in dysregulation of apoptosis and contribute to the molecular carcinogenesis of OSCC. Differential expressions of TRAIL receptors in OSCC do not appear to play a crucial role in their apoptotic rate or metastatic progression.

Identification and characterization of the herpes simplex virus type 2 gene encoding the essential capsid protein ICP32/VP19c.

We describe the characterization of the herpes simplex virus type 2 (HSV-2) gene encoding infected cell protein 32 (ICP32) and virion protein 19c (VP19c). We also demonstrate that the HSV-1 UL38/ORF.553 open reading frame (ORF), which has been shown to specify a viral protein essential for capsid formation (B. Pertuiset, M. Boccara, J. Cebrian, N. Berthelot, S. Chousterman, F. Puvian-Dutilleul, J. Sisman, and P. Sheldrick, J. Virol. 63: 2169-2179, 1989), must encode the cognate HSV type 1 (HSV-1) ICP32/VP19c protein. The region of the HSV-2 genome deduced to contain the gene specifying ICP32/VP19c was isolated and subcloned, and the nucleotide sequence of 2,158 base pairs of HSV-2 DNA mapping immediately upstream of the gene encoding the large subunit of the viral ribonucleotide reductase was determined. This region of the HSV-2 genome contains a large ORF capable of encoding two related 50,538- and 49,472-molecular-weight polypeptides. Direct evidence that this ORF encodes HSV-2 ICP32/VP19c was provided by immunoblotting experiments that utilized antisera directed against synthetic oligopeptides corresponding to internal portions of the predicted polypeptides encoded by the HSV-2 ORF or antisera directed against a TrpE/HSV-2 ORF fusion protein. The type-common immunoreactivity of the two antisera and comparison of the primary amino acid sequences of the predicted products of the HSV-2 ORF and the equivalent genomic region of HSV-1 provided evidence that the HSV-1 UL38 ORF encodes the HSV-1 ICP32/VP19c. Analysis of the expression of the HSV-1 and HSV-2 ICP32/VP19c cognate proteins indicated that there may be differences in their modes of synthesis. Comparison of the predicted structure of the HSV-2 ICP32/VP19c protein with the structures of related proteins encoded by other herpes viruses suggested that the internal capsid architecture of the herpes family of viruses varies substantially.

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

ABERRATIONS OF A PUTATIVE TUMOR SUPPRESSOR GENE SEL1L IN PANCREATIC DUCTAL ADENOCARCINOMA

Introduction: Pancreatic cancer is the fourth leading cause of cancer-related death among males and females in the United States. Sel-1-like (SEL1L) is a putative tumor suppressor gene that is downregulated in a significant proportion of human pancreatic ductal adenocarcinoma (PDAC). It was hypothesized that SEL1L expression could be down-modulated by somatic mutation, loss of heterozygosity (LOH), CpG island hypermethylation and/or aberrantly expressed microRNAs (miRNAs). Material and methods: In 42 PDAC tumors, the SEL1L coding region was amplified using reverse transcription polymerase chain reaction (RT-PCR), and analyzed by agarose gel electrophoresis and sequenced to search for mutations. Using fluorescent fragment analysis, two intragenic microsatellites in the SEL1L gene region were examined to detect LOH in a total of 73 pairs of PDAC tumors and normal-appearing adjacent tissues. Bisulfite DNA sequencing was performed to determine the methylation status of the SEL1L promoter in 41 PDAC tumors and 6 PDAC cell lines. Using real-time quantitative PCR, the expression levels of SEL1L mRNA and 7 aberrantly upregulated miRNAs that potentially target SEL1L were assessed in 42 PDAC tumor and normal pairs. Statistical methods were applied to evaluate the correlation between SEL1L mRNA and the miRNAs. Further the interaction was determined by functional analysis using a molecular biological approach. Results: No mutations were detected in the SEL1L coding region. More than 50% of the samples displayed abnormally alternate or aberrant spliced transcripts of SEL1L. About 14.5% of the tumors displayed LOH at the CAR/CAL microsatellite locus and 10.7% at the RepIN20 microsatellite locus. However, the presence of LOH did not show significant association with SEL1L downregulation. No methylation was observed in the SEL1L promoter. Statistical analysis showed that SEL1L mRNA expression levels significantly and inversely correlated with the expression of hsa-mir-143, hsa-mir-155, and hsa-mir-223. Functional analysis indicated that hsa-mir-155 acted as a suppressor of SEL1L in PL18 and MDAPanc3 PDAC cell lines. Discussion: Evidence from these studies suggested that SEL1L was possibly downregulated by aberrantly upregulated miRNAs in PDAC. Future studies should be directed towards developing a better understanding of the mechanisms for generation of aberrant SEL1L transcripts, and further analysis of miRNAs that may downregulate SEL1L.

A cladistic analysis of apolipoprotein B gene variation and plasma lipid and apolipoprotein B levels, using familial data

Any functionally important mutation is embedded in an evolutionary matrix of other mutations. Cladistic analysis, based on this, is a method of investigating gene effects using a haplotype phylogeny to define a set of tests which localize causal mutations to branches of the phylogeny. Previous implementations of cladistic analysis have not addressed the issue of analyzing data from related individuals, though in human studies, family data are usually needed to obtain unambiguous haplotypes. In this study, a method of cladistic analysis is described in which haplotype effects are parameterized in a linear model which accounts for familial correlations. The method was used to study the effect of apolipoprotein (Apo) B gene variation on total-, LDL-, and HDL-cholesterol, triglyceride, and Apo B levels in 121 French families. Five polymorphisms defined Apo B haplotypes: the signal peptide Insertion/deletion, Bsp 1286I, XbaI, MspI, and EcoRI. Eleven haplotypes were found, and a haplotype phylogeny was constructed and used to define a set of tests of haplotype effects on lipid and apo B levels.^ This new method of cladistic analysis, the parametric method, found significant effects for single haplotypes for all variables. For HDL-cholesterol, 3 clusters of evolutionarily-related haplotypes affecting levels were found. Haplotype effects accounted for about 10% of the genetic variance of triglyceride and HDL-cholesterol levels. The results of the parametric method were compared to those of a method of cladistic analysis based on permutational testing. The permutational method detected fewer haplotype effects, even when modified to account for correlations within families. Simulation studies exploring these differences found evidence of systematic errors in the permutational method due to the process by which haplotype groups were selected for testing.^ The applicability of cladistic analysis to human data was shown. The parametric method is suggested as an improvement over the permutational method. This study has identified candidate haplotypes for sequence comparisons in order to locate the functional mutations in the Apo B gene which may influence plasma lipid levels. ^