A MOLECULAR APPROACH FOR DETECTING GENETIC DAMAGE USING THE ECO R1 FAMILY OF HUMAN DNA (SEQUENCING)

A clone of the primary Eco R1 family of human DNA sequences has been used as an indicator sequence for detecting alterations induced by a toxic agent. Specific clones of this family have been examined and compared to the consensus sequence to determine the normal variability of this family. Though variations were observed, data indicated that such clones can be used to study induced DNA modifications. This DNA was exposed to the toxic agent dimethyl sulfate under various conditions and a distinct pattern of aberrations was shown to occur. It is suggested that this approach be used to characterize patterns of damage induced by various agents in the ultimate development of a system capable of monitoring human genotoxic exposure. ^

Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

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

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

Identification of the Causative Bacteria in Musculoskeletal Infections Using 16s rDNA - Denaturing Gradient Gel Electrophoresis Analysis

Musculoskeletal infections are infections of the bone and surrounding tissues. They are currently diagnosed based on culture analysis, which is the gold standard for pathogen identification. However, these clinical laboratory methods are frequently inadequate for the identification of the causative agents, because a large percentage (25-50%) of confirmed musculoskeletal infections are false negatives in which no pathogen is identified in culture. My data supports these results. The goal of this project was to use PCR amplification of a portion of the 16S rRNA gene to test an alternative approach for the identification of these pathogens and to assess the diversity of the bacteria involved. The advantages of this alternative method are that it should increase sample sensitivity and the speed of detection. In addition, bacteria that are non-culturable or in low abundance can be detected using this molecular technique. However, a complication of this approach is that the majority of musculoskeletal infections are polymicrobial, which prohibits direct identification from the infected tissue by DNA sequencing of the initial 16S rDNA amplification products. One way to solve this problem is to use denaturing gradient gel electrophoresis (DGGE) to separate the PCR products before DNA sequencing. Denaturing gradient gel electrophoresis (DGGE) separates DNA molecules based on their melting point, which is determined by their DNA sequence. This analytical technique allows a mixture of PCR products of the same length that electrophoreses through agarose gels as one band, to be separated into different bands and then used for DNA sequence analysis. In this way, the DGGE allows for the identification of individual bacterial species in polymicrobial-infected tissue, which is critical for improving clinical outcomes. By combining the 16S rDNA amplification and the DGGE techniques together, an alternative approach for identification has been used. The 16S rRNA gene PCR-DGGE method includes several critical steps: DNA extraction from tissue biopsies, amplification of the bacterial DNA, PCR product separation by DGGE, amplification of the gel-extracted DNA, and DNA sequencing and analysis. Each step of the method was optimized to increase its sensitivity and for rapid detection of the bacteria present in human tissue samples. The limit of detection for the DNA extraction from tissue was at least 20 Staphylococcus aureus cells and the limit of detection for PCR was at least 0.05 pg of template DNA. The conditions for DGGE electrophoreses were optimized by using a double gradient of acrylamide (6 – 10%) and denaturant (30-70%), which increased the separation between distinct PCR products. The use of GelRed (Biotium) improved the DNA visualization in the DGGE gel. To recover the DNA from the DGGE gels the gel slices were excised, shredded in a bead beater, and the DNA was allowed to diffuse into sterile water overnight. The use of primers containing specific linkers allowed the entire amplified PCR product to be sequenced and then analyzed. The optimized 16S rRNA gene PCR-DGGE method was used to analyze 50 tissue biopsy samples chosen randomly from our collection. The results were compared to those of the Memorial Hermann Hospital Clinical Microbiology Laboratory for the same samples. The molecular method was congruent for 10 of the 17 (59%) culture negative tissue samples. In 7 of the 17 (41%) culture negative the molecular method identified a bacterium. The molecular method was congruent with the culture identification for 7 of the 33 (21%) positive cultured tissue samples. However, in 8 of the 33 (24%) the molecular method identified more organisms. In 13 of the 15 (87%) polymicrobial cultured tissue samples the molecular method identified at least one organism that was also identified by culture techniques. Overall, the DGGE analysis of 16S rDNA is an effective method to identify bacteria not identified by culture analysis.

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.

Genetic and physical analysis of the bacteriophage P22 tail protein

A plasmid based genetic system was developed for the tail protein of the Salmonella typhimurium bacteriophage P22 and used to isolate and characterize tail protein mutants. The tail protein is a trimeric structural protein of the phage and an endorhamnosidase whose activity is essential for infection. The gene for the tail protein has previously been cloned into a plasmid expression vector and sequenced. A plate complementation assay for tail protein produced from the cloned gene was developed and used to isolate 27 tail protein mutants following mutagenesis of the cloned gene. These mutations were mapped into 12 deletion intervals using deletions which were made on plasmids in vitro and crossed onto P22. The base substitutions were determined by DNA sequencing. The majority of mutants had missense or nonsense mutations in the protein coding portion of the gene; however four of the mutants were in the putative transcription terminator. The oligomeric state of tail protein from the 15 missense mutants was investigated using SDS and nondenaturing polyacrylamide gel electrophoresis of cell lysates. Wild-type tail protein retains its trimeric structure in SDS gels at room temperature. Two of the mutant proteins also migrated as trimers in SDS gels, yet one of these had a considerably faster mobility than wild-type trimer. Its migration was the same as wild-type in a nondenaturing gel, so it is thought to be a trimer which is partially denatured by SDS. Four of the mutants produced proteins which migrate at the position of a monomer in an SDS gel but cannot be seen on a nondenaturing gel. These proteins are thought to be either monomers or soluble aggregates which cannot enter the nondenaturing gel. The remainder of mutants produce protein which is degraded. The mutant tail protein which had normal trimeric mobility on SDS and nondenaturing gels was purified. This protein has essentially wild-type ability to attach to phage capsids, but its endorhamnosidase activity is only 4% of wild-type. ^

Mutagenicity of herpes simplex virus type 1 in a shuttle vector

The shuttle vector plasmid pZ189 was used to find the kinds of mutations that are induced by herpes simplex virus type-1 (HSV-1). In cells infected by HSV-1 the frequency of mutation in supF gene, the mutagenesis marker, was increased over background by from two- to seven-fold, reaching 0.14-0.45%. No increase was induced by infection by vaccinia virus under the same conditions. Mutagenesis was an early event, showing a four-fold increase in mutation frequency at only two hours after infection, and peaking at a seven-fold increase at four hours after infection. DNA sequencing and gel electrophoresis analysis were performed on 105 HSV-1 induced mutants and 65 spontaneous mutants and provided the following information: (1) A change in plasmid size was seen in 54% of HSV-1 related mutants, compared with only 37% of spontaneous mutants. (2) Among point mutations, the predominant type was G:C to A:T transition, which accounted for 51% of point mutations in mutants isolated from cells infected with HSV-1, and 32% of point mutations in spontaneous mutants. (3) Deletions of DNA were seen in HSV-1 related mutants at a frequency of 40%, compared with 29% in spontaneous mutants. The HSV-1 related deletions were about half the length of spontaneous mutants and three contained short filler sequences. (4) Fifteen (15%) of HSV-1 induced mutants revealed the altered restriction patterns on agarose gel electrophoresis analysis and were due either to rearrangements of plasmid DNA, and/or to insertion of sequences derived from chromosomal DNA (seven plasmids). No insertions of DNA from HSV-1 were detected. Among spontaneous mutants, only 5 (7.7%) were rearrangements and none had inserted chromosomal DNA. (5) DNA sequence analysis of seven plasmids with inserted chromosomal DNA revealed that four cases had repetitive DNA sequences integrated and the other three were unidentified sequences from the GenBank database. Three repetitive DNA included $\alpha$ satellite, Alu and KpnI family sequences. The other sequence was identified as tRNA-like component. The observed mutations have implications for the mechanism of malignant transformation of cells by HSV-1. ^

Identification of cadmium-induced mutations in a mammalian cell line

Epidemiological studies have shown cadmium to induce cancer in humans, while experimental studies have proven this metal to be a potent tumor inducer in animals. However, cadmium appears nonmutagenic in most prokaryotic and eukaryotic mutagenesis assays. In this study, we present the identification of mutations in normal rat kidney cells infected with the mutant MuSVts110 retrovirus (6m2 cells) as a result of treatment with cadmium chloride. The detection of these mutations was facilitated by the use of a novel mutagenesis assay established in this laboratory. The 6m2 reversion assay is a positive selection system based on the conditional expression of the MuSVts110 v-mos gene. In MuSVts110 the gag and mos genes are fused out of frame, thus the translation of the v-mos sequence requires a frameshift in the genomic RNA. In 6m2 cells this frameshift is accomplished by the temperature-dependent splicing of the primary MuSVts110 transcript. Splicing of MuSVts110, which is mediated by cis-acting sequences, occurs when 6m2 cells are grown at 33$\sp\circ$C and below, but not at 39$\sp\circ$C. Therefore, 6m2 cells appear transformed at low growth temperatures, but take on a morphologically normal appearance when grown at high temperatures. The treatment of 6m2 cells with cadmium chloride resulted in the outgrowth of a number of cells that reverted to the transformed state at high growth temperatures. Analysis of the viral proteins expressed in these cadmium-induced 6m2 revertants suggested that they contained mutations in their MuSVts110 DNA. Sequencing of the viral DNA from three revertants that constitutively expressed the P85$\sp{gag{-}mos}$ transforming protein revealed five different mutations. The Cd-B2 revertant contained three of those mutations: an A-to-G transition 48 bases downstream of the MuSVts110 3$\sp\prime$ splice site, plus a G-to-T and an A-to-T transversion 84 and 100 bases downstream of the 5$\sp\prime$ splice site, respectively. The Cd-15-5 revertant also contained a point mutation, a T-to-C transition 46 bases downstream of the 5$\sp\prime$ splice site, while Cd-10-5 contained a three base deletion of MuSVts110 11 bases upstream of the 3$\sp\prime$ splice site. A fourth revertant, Cd-10, expressed a P100$\sp{gag{-}mos}$ transforming protein, and was found to have a two base deletion. This deletion accomplished the frameshift necessary for v-mos expression, but did not alter MuSVts110 RNA splicing and the expression of p85$\sp{gag{-}mos}.$ Lastly, sequencing of the MuSVts110 DNA from three spontaneous revertants revealed the same G to T transversion in each one. This was the same mutation that was found in the Cd-B2 revertant. These findings provide the first example of mutations resulting from exposure to cadmium and suggest, by the difference in each mutation, the complexity of the mechanism utilized by cadmium to induce DNA damage. ^

Positional candidate cloning of the RP10 form of retinitis pigmentosa

Retinitis pigmentosa (RP) is a genetically heterogeneous group of retinal degenerations that affects over one million people worldwide. To date, 11 autosomal dominant, 13 autosomal recessive, and 5 X-linked forms of retinitis pigmentosa have been identified through linkage analysis, but the disease-causing genes and mutations have been found for only half of these loci. My research uses a positional candidate cloning approach to identify the gene and mutations responsible for one type of autosomal dominant retinitis pigmentosa, RP10. The premise is that identifying the genes and mutations responsible for disease will provide insight into disease mechanisms and provide treatment options. Previous research mapped the RP10 locus to a 5cM region on chromosome 7q31 between markers D7S686 and D7S530. Linkage and fine-point haplotype analysis was used to reduce and refine the RP10 disease interval to a 4cM region located between D7S2471 and a new marker located 45,000bp telomeric of D7S461. In order to identify genes located in the RP10 interval, an extensive EST map was created of this region. Five EST clusters from this map were analyzed to determine if mutations in these genes cause the RP10 form of retinitis pigmentosa. The genomic structure of a known metabotrophic glutamate receptor, GRMS8, was determined first. DNA sequencing of GRM8 in RP10 family members did not identify any disease-causing mutations. Four other EST clusters (A170, A173, A189, and A258) were characterized and determined to be part of the same gene, UBNL1 (ubinuclein-like 1). The full-length mRNA sequence and genomic structure of UBNL1 was determined and then screened in patients. No disease-causing mutations were identified in any of the RP10 family members tested. Recent data made available with the release of the public and Celera genome assemblies indicates that UBNL1 is outside of the RP10 disease region. Despite this complication, characterization of UBNL1 is still important in the understanding of normal visual processes and it is possible that mutations in UBNL1 could cause other forms of retinopathy. The EST map and list of RP10 candidates will continue to aid others in the search for the RP10 gene and mutations. ^

Analysis of Mycobacterium tuberculosis in the state of Texas for rifampin resistance using molecular beacons

Mycobacterium tuberculosis, a bacillus known to cause disease in humans since ancient times, is the etiological agent of tuberculosis (TB). The infection is primarily pulmonary, although other organs may also be affected. The prevalence of pulmonary TB disease in the US is highest along the US-Mexico border, and of the four US states bordering Mexico, Texas had the second highest percentage of cases of TB disease among Mexico-born individuals in 1999 (CDC, 2001). Between the years of 1993 and 1998, the prevalence of drug-resistant (DR) TB was 9.1% among Mexican-born individuals and 4.4% among US-born individuals (CDC, 2001). In the same time period, the prevalence of multi-drug resistant (MDR) TB was 1.4% among Mexican-born individuals and 0.6% among US-born individuals (CDC, 2001). There is a renewed urgency in the quest for faster and more effective screening, diagnosis, and treatment methods for TB due to the resurgence of tuberculosis in the US during the mid-1980s and early 1990s (CDC, 2007a), and the emergence of drug-resistant, multidrug-resistant, and extremely drug-resistant tuberculosis worldwide. Failure to identify DR and MDR-TB quickly leads to poorer treatment outcomes (CDC, 2007b). The recent rise in TB/HIV comorbidity further complicates TB control efforts. The gold standard for identification of DR-TB requires mycobacterial growth in culture, a technique taking up to three weeks, during which time DR/MDR-TB individuals harboring resistant organisms may be receiving inappropriate treatment. The goal of this study was to determine the sensitivity and specificity of real-time quantitative polymerase chain reaction (qPCR) using molecular beacons in the Texas population. qPCR using molecular beacons is a novel approach to detect mycobacterial mutations conferring drug resistance. This technique is time-efficient and has been shown to have high sensitivity and specificity in several populations worldwide. Rifampin (RIF) susceptibility was chosen as the test parameter because strains of M. tuberculosis which are resistant to RIF are likely to also be MDR. Due to its status as a point of entry for many immigrants into the US, control efforts against TB and drug-resistant TB in Texas is a vital component of prevention efforts in the US as a whole. We show that qPCR using molecular beacons has high sensitivity and specificity when compared with culture (94% and 87%, respectively) and DNA sequencing (90% and 96%, respectively). We also used receiver operator curve analysis to calculate cutoff values for the objective determination of results obtained by qPCR using molecular beacons. ^

Molecular and genetic analysis of Myxococcus xanthus early *development

To identify more mutations that can affect the early development of Myxococcus xanthus, the synthetic transposon TnT41 was designed and constructed. By virtue of its special features, it can greatly facilitate the processes of mutation screening/selection, mapping, cloning and DNA sequencing. In addition, it allows for the systematic discovery of genes in regulatory hierarchies using their target promoters. In this study, the minimal regulatory region of the early developmentally regulated gene 4521 was used as a reporter in the TnT41 mutagenesis. Both positive (P) mutations and negative (N) mutations were isolated based on their effects on 4521 expression.^ Four of these mutations, i.e. N1, N2, P52 and P54 were analyzed in detail. Mutations N1 and N2 are insertion mutations in a gene designated sasB. The sasB gene is also identified in this study by genetic and molecular analysis of five UV-generated 4521 suppressor mutations. The sasB gene encodes a protein without meaningful homology in the databases. The sasB gene negatively regulates 4521 expression possibly through the SasS-SasR two component system. A wild-type sasB gene is required for normal M. xanthus fruiting body formation and sporulation.^ Cloning and sequencing analysis of the P52 mutation led to the identification of an operon that encodes the M. xanthus high-affinity branched-chain amino acid transporter system. This liv operon consists of five genes designated livK, livH, livM, livC, and livF, respectively. The Liv proteins are highly similar to their counterparts from other bacteria in both amino acid sequences, functional motifs and predicted secondary structures. This system is required for development since liv null mutations cause abnormality in fruiting body formation and a 100-fold decrease in sporulation efficiency.^ Mutation P54 is a TnT41 insertion in the sscM gene of the ssc chemotaxis system, which has been independently identified by Dr. Shi's lab. The sscM gene encodes a MCP (methyl-accepting chemotaxis protein) homologue. The SscM protein is predicted to contain two transmembrane domains, a signaling domain and at least one putative methylation site. Null mutations of this gene abolish the aggregation of starving cells at a very early stage, though the sporulation levels of the mutant can reach 10% that of wild-type cells. ^

POLG DNA testing as an emerging standard of care before instituting valproic acid therapy for pediatric seizure disorders.

PURPOSE: To review our clinical experience and determine if there are appropriate signs and symptoms to consider POLG sequencing prior to valproic acid (VPA) dosing in patients with seizures. METHODS: Four patients who developed VPA-induced hepatotoxicity were examined for POLG sequence variations. A subsequent chart review was used to describe clinical course prior to and after VPA dosing. RESULTS: Four patients of multiple different ethnicities, age 3-18 years, developed VPA-induced hepatotoxicity. All were given VPA due to intractable partial seizures. Three of the patients had developed epilepsia partialis continua. The time from VPA exposure to liver failure was between 2 and 3 months. Liver failure was reversible in one patient. Molecular studies revealed homozygous p.R597W or p.A467T mutations in two patients. The other two patients showed compound heterozygous mutations, p.A467T/p.Q68X and p.L83P/p.G888S. Clinical findings and POLG mutations were diagnostic of Alpers-Huttenlocher syndrome. CONCLUSION: Our cases underscore several important findings: POLG mutations have been observed in every ethnic group studied to date; early predominance of epileptiform discharges over the occipital region is common in POLG-induced epilepsy; the EEG and MRI findings varying between patients and stages of the disease; and VPA dosing at any stage of Alpers-Huttenlocher syndrome can precipitate liver failure. Our data support an emerging proposal that POLG gene testing should be considered in any child or adolescent who presents or develops intractable seizures with or without status epilepticus or epilepsia partialis continua, particularly when there is a history of psychomotor regression.

Sequencing of TGF-beta pathway genes in familial cases of intracranial aneurysm.

BACKGROUND AND PURPOSE: Familial aggregation of intracranial aneurysms (IA) strongly suggests a genetic contribution to pathogenesis. However, genetic risk factors have yet to be defined. For families affected by aortic aneurysms, specific gene variants have been identified, many affecting the receptors to transforming growth factor-beta (TGF-beta). In recent work, we found that aortic and intracranial aneurysms may share a common genetic basis in some families. We hypothesized, therefore, that mutations in TGF-beta receptors might also play a role in IA pathogenesis. METHODS: To identify genetic variants in TGF-beta and its receptors, TGFB1, TGFBR1, TGFBR2, ACVR1, TGFBR3, and ENG were directly sequenced in 44 unrelated patients with familial IA. Novel variants were confirmed by restriction digestion analyses, and allele frequencies were analyzed in cases versus individuals without known intracranial disease. Similarly, allele frequencies of a subset of known SNPs in each gene were also analyzed for association with IA. RESULTS: No mutations were found in TGFB1, TGFBR1, TGFBR2, or ACVR1. Novel variants identified in ENG (p.A60E) and TGFBR3 (p.W112R) were not detected in at least 892 reference chromosomes. ENG p.A60E showed significant association with familial IA in case-control studies (P=0.0080). No association with IA could be found for any of the known polymorphisms tested. CONCLUSIONS: Mutations in TGF-beta receptor genes are not a major cause of IA. However, we identified rare variants in ENG and TGFBR3 that may be important for IA pathogenesis in a subset of families.