DNA nucleotide excision repair gene single nucleotide polymorphisms and hereditary nonpolyposis colorectal cancer

Hereditary nonpolyposis colorectal cancer (HNPCC) is an autosomal dominant disease caused by germline mutations in DNA mismatch repair(MMR) genes. The nucleotide excision repair(NER) pathway plays a very important role in cancer development. We systematically studied interactions between NER and MMR genes to identify NER gene single nucleotide polymorphism (SNP) risk factors that modify the effect of MMR mutations on risk for cancer in HNPCC. We analyzed data from polymorphisms in 10 NER genes that had been genotyped in HNPCC patients that carry MSH2 and MLH1 gene mutations. The influence of the NER gene SNPs on time to onset of colorectal cancer (CRC) was assessed using survival analysis and a semiparametric proportional hazard model. We found the median age of onset for CRC among MMR mutation carriers with the ERCC1 mutation was 3.9 years earlier than patients with wildtype ERCC1(median 47.7 vs 51.6, log-rank test p=0.035). The influence of Rad23B A249V SNP on age of onset of HNPCC is age dependent (likelihood ratio test p=0.0056). Interestingly, using the likelihood ratio test, we also found evidence of genetic interactions between the MMR gene mutations and SNPs in ERCC1 gene(C8092A) and XPG/ERCC5 gene(D1104H) with p-values of 0.004 and 0.042, respectively. An assessment using tree structured survival analysis (TSSA) showed distinct gene interactions in MLH1 mutation carriers and MSH2 mutation carriers. ERCC1 SNP genotypes greatly modified the age onset of HNPCC in MSH2 mutation carriers, while no effect was detected in MLH1 mutation carriers. Given the NER genes in this study play different roles in NER pathway, they may have distinct influences on the development of HNPCC. The findings of this study are very important for elucidation of the molecular mechanism of colon cancer development and for understanding why some mutation carriers of the MSH2 and MLH1 gene develop CRC early and others never develop CRC. Overall, the findings also have important implications for the development of early detection strategies and prevention as well as understanding the mechanism of colorectal carcinogenesis in HNPCC. ^

The effect of exposure to antibiotics on incidence and spontaneous clearance of childhood Helicobacter pylori infection

It is claimed often in the H. pylori literature that spontaneous clearance (infection loss without attempts to treat) is uncommon, though little evidence supports this claim. Emerging evidence suggests that spontaneous clearance may be frequent in young children; however, factors that determine persistence of untreated H. pylori infection in childhood are not well understood. The author hypothesized that antibiotics taken for common infections cause spontaneous clearance of H. pylori infection in children. The Pasitos Cohort Study (19982005) investigated predictors of acquisition and persistence of H. pylori infection in children from El Paso, Texas, and Juarez, Mexico, enrolled prenatally at maternal-child clinics. Children were screened for infection at target intervals of 6 months from 6-84 months of age by the 13C-urea breath test corrected for body-size-dependent variation in CO2 production. This dissertation aimed to estimate the risk of spontaneous clearance at the next test following an initial detected H. pylori infection (first detected clearance), estimate the effect of antibiotic exposure on the risk of first detected clearance (risk difference), and estimate the effect of antibiotic exposure on the rate of first detected infection (rate ratio). Data on infection status and medication history were available for 608 children followed for a mean of 3.5 years. Among 265 subjects with a first detected infection, 218 had a subsequent test, and among them, the risk of first detected clearance was 68% (95% CI: 61-74%). Children who took antibiotics during the interval between first detected infection and next test had an increased probability (risk difference of 10 percentage points) of a first detected clearance. However, there was also a similar effect of average antibiotic use >0 courses across all intervals preceding the next test. Average antibiotic exposure across all intervals preceding the first detected infection appeared to have a much stronger protective effect than interval/specific exposure when estimating incidence rate ratios (0.45 vs. 1.0). Incidental antibiotic exposure appears to influence the acquisition and duration of childhood H. pylori infection, however, given that many exposed children acquired the infection and many unexposed children cleared the infection, antibiotic exposure does not explain all infection events. ^

Characterization of glucocorticoid receptor (NR3C1) gene polymorphisms and a family-based association study with hypertension

Hypertension is usually defined as having values of systolic blood pressure ≥140 mmHg, diastolic blood pressure ≥90 mmHg. Hypertension is one of the main adverse effects of glucocorticoid on the cardiovascular system. Glucocorticoids are essential hormones, secreted from adrenal glands in circadian fashion. Glucocorticoid's effect on blood pressure is conveyed by the glucocorticoid receptor (NR3C1), an omnipresent nuclear transcription factor. Although polymorphisms in this gene have long been implicated to be a causal factor for cardiovascular diseases such as hypertension, no study has yet thoroughly interrogated the gene's polymorphisms for their effect on blood pressure levels. Therefore, I have first resequenced ∼30 kb of the gene, encompassing all exons, promoter regions, 5'/3' UTRs as well as at least 1.5 kb of the gene's flanking regions from 114 chromosome 5 monosomic cell lines, comprised of three major American ethnic groups—European American, African American and Mexican American. I observed 115 polymorphisms and 14 common molecularly phased haplotypes. A subset of markers was chosen for genotyping study populations of GENOA (Genetic Epidemiology Network of Atherosclerosis; 1022 non-Hispanic whites, 1228 African Americans and 954 Mexican Americans). Since these study populations include sibships, the family-based association test was performed on 4 blood pressure-related quantitative variables—pulse, systolic blood pressure, diastolic blood pressure and mean arterial pressure. Using these analyses, multiple correlated SNPs are significantly protective against high systolic blood pressure in non-Hispanic whites, which includes rsb198, a SNP formerly associated with beneficial body compositions. Haplotype association analysis also supports this finding and all p-values remained significant after permutation tests. I therefore conclude that multiple correlated SNPs on the gene may confer protection against high blood pressure in non-Hispanic whites. ^

The effect of DNA methylation in carcinogenesis

CpG island methylation within single gene promoters can silence expression of associated genes. We first extended these studies to bidirectional gene pairs controlled by single promoters. We showed that hypermethylation of bidirectional promoter-associated CpG island silences gene pairs (WNT9A/CD558500, CTDSPL/BC040563, and KCNK15/BF 195580) simultaneously. Hypomethylation of these promoters by 5-aza-2'-deoxycytidine treatment reactivated or enhanced gene expression bidirectionally. These results were further confirmed by luciferase assays. Methylation of WNT9A/CD558500 and CTDSPL/BC040563 promoters occurs frequently in primary colon cancers and acute lymphoid leukemia, respectively. ^ Next we sought to understand the origins of hypermethylation in cancer. CpG islands associated with tumor suppressor genes are normally free from methylation, but can be hypermethylated in cancer. It remains poorly understood how these genes are protected from methylation in normal tissues. In our studies, we aimed to determine if cis-acting elements in these genes are responsible for this protection, using the tumor suppressor gene p16 as a model. We found that Alu repeats located both upstream and downstream of the p16 promoter become hypermethylated with age. In colon cancer samples, the methylation level is particularly high, and the promoter can also be affected. Therefore, the protection in the promoter against methylation spreading could fail during tumorigenesis. This methylation pattern in p16 was also observed in cell lines of different tissue origins, and their methylation levels were found to be inversely correlated with that of active histone modification markers (H3K4-3me and H3K9-Ac). To identify the mechanism of protection against methylation spreading, we constructed serial deletions of the p16 protected region and used silencing of a neomycin reporter gene to evaluate the protective effects of these fragments. A 126 bp element was identified within the region which exerts bidirectional protection against DNA methylation, independently of its transcriptional activity. The protective strength of this element is comparable to that of the HS4 insulator. During long-term culture, the presence of this element significantly slowed methylation spreading. In conclusion, we have found that an element located in the p16 promoter is responsible for protection against DNA methylation spreading in normal tissues. The failure of protective cis-elements may be a general feature of tumor-suppressor gene silencing during tumorigenesis. ^

The molecular basis for beta-lactamase gene expression in B. anthracis, B. cereus and B. thuringiensis

The susceptibility of most Bacillus anthracis strains to β-lactam antibiotics is intriguing considering that the B. anthracis genome harbors two β-lactamase genes, bla1 and bla2, and closely-related species, Bacillus cereus and Bacillus thuringiensis, typically produce β-lactamases. This work demonstrates that B. anthracis bla expression is affected by two genes, sigP and rsp, predicted to encode an extracytoplasmic function sigma factor and an antisigma factor, respectively. Deletion of the sigP/rsp locus abolished bla expression in a penicillin-resistant clinical isolate and had no effect on bla expression in a prototypical penicillin-susceptible strain. Complementation with sigP/rsp from the penicillin-resistant strain, but not the penicillin-susceptible strain, conferred β-lactamase activity upon both mutants. These results are attributed to a nucleotide deletion near the 5' end of rsp in the penicillin-resistant strain that is predicted to result in a nonfunctional protein. B. cereus and B. thuringiensis sigP and rsp homologues are required for inducible penicillin resistance in those species. Expression of the B. cereus or B. thuringiensis sigP and rsp genes in a B. anthracis sigP/rsp-null mutant confers resistance to β-lactam antibiotics, suggesting that while B. anthracis contains the genes necessary for sensing β-lactam antibiotics, the B. anthracis sigP/rsp gene products are insufficient for bla induction. ^ Because alternative sigma factors recognize unique promoter sequence, direct targets can be elucidated by comparing transcriptional profiling results with an in silico search using the sigma factor binding sequence. Potential σP -10 and -35 promoter elements were identified upstream from bla1 bla2 and sigP. Results obtained from searching the B. anthracis genome with the conserved sequences were evaluated against transcriptional profiling results comparing B. anthracis 32 and an isogenic sigP/rsp -null strain. Results from these analyses indicate that while the absence of the sigP gene significantly affects the transcript levels of 16 genes, only bla1, bla2 and sigP are directly regulated by σP. The genomes of B. cereus and B. thuringiensis strains were also analyzed for the potential σP binding elements. The sequence was located upstream from the sigP and bla genes, and previously unidentified genes predicted to encode a penicillin-binding protein (PBP) and a D-alanyl-D-alanine carboxypeptidase, indicating that the σ P regulon in these species responds to cell-wall stress caused by β-lactam antibiotics. ^ β-lactam antibiotics prevent attachment of new peptidoglycan to the cell wall by blocking the active site of PBPs. A B. cereus and B. thuringiensis pbp-encoding gene located near bla1 contains a potential σP recognition sequence upstream from the annotated translational start. Deletion of this gene abolished β-lactam resistance in both strains. Mutations in the active site of the PBP were detrimental to β-lactam resistance in B. cereus, but not B. thuringiensis, indicating that the transpeptidase activity is only important in B. cereus. I also found that transcript levels of the PBP-encoding gene are not significantly affected by the presence of β-lactam antibiotic. Based on these data I hypothesize that the gene product acts a sensor of β-lactam antibiotic. ^

The effect of genotype-by-environment interaction on longitudinal triglyceride profiles in the Bogalusa Heart Study

Triglyceride levels are a component of plasma lipids that are thought to be an important risk factor for coronary heart disease and are influenced by genetic and environmental factors, such as single nucleotide polymorphisms (SNPs), alcohol intake, and smoking. This study used longitudinal data from the Bogalusa Heart Study, a biracial community-based survey of cardiovascular disease risk factors. A sample of 1191 individuals, 4 to 38 years of age, was measured multiple times from 1973 to 2000. The study sample consisted of 730 white and 461 African American participants. Individual growth models were developed in order to assess gene-environment interactions affecting plasma triglycerides over time. After testing for inclusion of significant covariates and interactions, final models, each accounting for the effects of a different SNP, were assessed for fit and normality. After adjustment for all other covariates and interactions, LIPC -514C/T was found to interact with age3, age2, and age and a non-significant interaction of CETP -971G/A genotype with smoking status was found (p = 0.0812). Ever-smokers had higher triglyceride levels than never smokers, but persons heterozygous at this locus, about half of both races, had higher triglyceride levels after smoking cessation compared to current smokers. Since tobacco products increase free fatty acids circulating in the bloodstream, smoking cessation programs have the potential to ultimately reduce triglyceride levels for many persons. However, due to the effect of smoking cessation on the triglyceride levels of CETP -971G/A heterozygotes, the need for smoking prevention programs is also demonstrated. Both smoking cessation and prevention programs would have a great public health impact on minimizing triglyceride levels and ultimately reducing heart disease. ^

Mechanism of methylation gene reactivation in human cancer cells

Epigenetic silencing of tumor suppressor genes by DNA hypermethylation at promoter regions is a common event in carcinogenesis and tumor progression. Abrogation of methylation and reversal of epigenetic silencing is a very potent way in cancer treatment. However, the reactivation mechanisms are poorly understood. In this study, we first developed a cell line model system named YB5, derived from SW48 cancer cell line, which bears one copy of stably integrated EGFP gene on Chromosome 1p31.1 region. The GFP gene expression is transcriptionally silenced due to the hypermethylated promoter CMV. However, the GFP expression can be restored using demethylating agent 5-aza-2' deoxycytidine (DAC), and detected by FACS and fluorescent microscopy. Using this system, we observed the heterogeneous reactivation induced by DAC treatment. After flow sorting, GFP negative cells exhibited similar level of incomplete demethylation compared to GFP positive cells on repetitive LINE1 element, tumor suppressor genes such as P16, CDH13, and RASSF1a, and CMV promoter as well. However, the local chromatin of CMV-GFP locus altered to an open structure marked by high H3 lysine 9 acetylation and low H3 lysine 27 tri-methylation in GFP positive cells, while the GFP negative cells retained mostly the original repressive marks. Thus, we concluded that DAC induced DNA hypomethylation alone does not directly determine the level of re-expression, and the resetting of the local chromatin structure under hypomethylation environment is required for gene reactivation. Besides, a lentivirus vector-based shRNA screening was performed using the YB5 system. Although it is the rare chance that vector lands in the neighboring region of GFP, we found that the exogenous vector DNA inserted into the upstream region of GFP gene locus led to the promoter demethylation and reactivated the silenced GFP gene. Thus, epigenetic state can be affected by changing of the adjacent nucleic acid sequences. Further, this hypermethylation silenced system was utilized for epigenetic drug screening. We have found that DAC combined with carboplatin would enhance the GFP% yield and increase expression of other tumor suppressor genes than DAC alone, and this synergistic effect may be related to DNA repair process. In summary, these studies reveal that reversing of methylation silencing requires coordinated alterations of DNA methylation, chromatin structure, and local microenvironment. ^

Is maternal separation anxiety associated with school readiness?

Background. Increasing rates of maternal employment highlight the need for non-maternal child care for infants at an earlier age. Several studies have shown that employment induced maternal depression or psychological distress is associated with the child's socio-emotional and cognitive development. However, separation anxiety, a common phenomenon observed among employed mothers during early years, has seldom been studied. Therefore, the purpose of this study was to evaluate the role of maternal separation anxiety in the child's cognitive development.^ Methods. Data were obtained from Phase I (birth to 36 months) of the National Institute of Child Health and Human Development, Study of Early Child Care and Youth Development (NICHD SECCYD). Bivariate and multivariate analyses were performed to determine the association between separation anxiety groups and child outcomes. Multivariate analysis was also used to examine the mediating and/or moderating effect of sensitivity and moderating effect of difficult temperament.^ Results. Separation anxiety showed a negative association with the Bracken, attachment security, maternal sensitivity and psychological state. Children whose mothers never reported high levels of separation anxiety showed higher levels of school readiness and attachment security compared to those whose mothers experienced high levels of separation anxiety at least once. There was a significant interaction between separation anxiety and maternal sensitivity for the Bracken and attachment security indicating the moderating effect of sensitivity. Maternal sensitivity was also found to partially mediate the association between high levels of separation anxiety and school readiness or attachment security. However, the interaction between difficult temperament and separation anxiety was not significant for any of the child outcomes. ^ Conclusions. High levels of separation anxiety have a negative impact on school readiness, attachment security, maternal sensitivity and psychological state. In addition, mothers who experience high levels of separation anxiety but are sensitive during the mother-child interaction have children with high school readiness and attachment security compared to those who are less sensitive.^ Keywords. Maternal separation anxiety, School readiness. ^

Association of maternal sphinganine:sphingosine ratio and folate levels with neural tube defects along the South Texas border

Neural tube defects (NTDs) remain elevated in Hispanic women along the South Texas Border, despite folate supplementation and folate fortification of cereal products. Missmer et al. examined the relationships between fumonisins, a class of corn mycotoxin, and NTDs in Hispanic women who ate corn tortillas and found increased odds ratios with increasing exposure, as measured by serum sphinganine:sphingosine (sa:so) ratios. This study examined the interactions between categorized maternal serum folate levels and stratified sa:so ratios and the resultant odds ratios of NTDs, stratified by type (anencephaly and spina bifida). The hypothesis was that the above normal folate category would have lower odds ratios of NTDs at given sa:so ratio categories and that there would be a difference in odds ratio patterns for anencephaly and spina bifida. Methods. Data for 406 Hispanic women were obtained from the Missmer case-control study. Sa:so ratios were calculated and subjects were stratified into “below normal,” “normal,” and above normal range for folate. A logistic regression model was applied, controlling for BMI, serum B12, lab batch, and conception date. Results. While OR’s of NTDs increased for increasing sa:so ratios, OR’s for “above normal” folate were not decreased at any sa:so ratio and there was no statistically significant difference between OR’s of anencephaly and spina bifida. Conclusion. Folate does not appear to be protective against the potential teratogenic effect of fumonisins and did not differ in effect on OR’s of NTD by type. More research is necessary to determine the extent of fumonisin exposure in Hispanic women along the South Texas Border.^

Genetic epidemiology of gallbladder disease in Mexican Americans and cholesterol 7$\alpha$-hydroxylase gene sequence variation

Among Mexican Americans, the second largest minority group in the United States, the prevalence of gallbladder disease is markedly elevated. Previous data from both genetic admixture and family studies indicate that there is a genetic component to the occurrence of gallbladder disease in Mexican Americans. However, prior to this thesis no formal genetic analysis of gallbladder disease had been carried out nor had any contributing genes been identified.^ The results of complex segregation analysis in a sample of 232 Mexican American pedigrees documented the existence of a major gene having two alleles with age- and gender-specific effects influencing the occurrence of gallbladder disease. The estimated frequency of the allele increasing susceptibility was 0.39. The lifetime probabilities that an individual will be affected by gallbladder disease were 1.0, 0.54, and 0.00 for females of genotypes "AA", "Aa", and "aa", respectively, and 0.68, 0.30, and 0.00 for males, respectively. This analysis provided the first conclusive evidence for the existence of a common single gene having a large effect on the occurrence of gallbladder disease.^ Human cholesterol 7$\alpha$-hydroxylase is the rate-limiting enzyme in bile acid synthesis. The results of an association study in both a random sample and a matched case/control sample showed that there is a significant association between cholesterol 7$\alpha$-hydroxylase gene variation and the occurrence of gallbladder disease in Mexican Americans males but not in females. These data have implicated a specific gene, 7$\alpha$-hydroxylase, in the etiology of gallbladder disease in this population.^ Finally, I asked whether the inferred major gene from complex segregation analysis is genetically linked to the cholesterol 7$\alpha$-hydroxylase gene. Three pedigrees predicted to be informative for linkage analysis by virtue of supporting the major gene hypothesis and having parents with informative genotypes and multiple offspring were selected for this linkage analysis. In each of these pedigrees, the recombination fractions maximized at 0 with a positive, albeit low, LOD score. The results of this linkage analysis provide preliminary and suggestive evidence that the cholesterol 7$\alpha$-hydroxylase gene and the inferred gallbladder disease susceptibility gene are genetically linked. ^

The epidemiology of left ventricular outflow tract obstruction defects in Texas, 2002-2008: An update and assessment of effect heterogeneity

Left ventricular outflow tract (LVOT) defects are an important group of congenital heart defects (CHDs) because of their associated mortality and long-term complications. LVOT defects include aortic valve stenosis (AVS), coarctation of aorta (CoA), and hypoplastic left heart syndrome (HLHS). Despite their clinical significance, their etiology is not completely understood. Even though the individual component phenotypes (AVS, CoA, and HLHS) may have different etiologies, they are often "lumped" together in epidemiological studies. Though "lumping" of component phenotypes may improve the power to detect associations, it may also lead to ambiguous findings if these defects are etiologically distinct. This is due to potential for effect heterogeneity across component phenotypes. ^ This study had two aims: (1) to identify the association between various risk factors and both the component (i.e., split) and composite (i.e., lumped) LVOT phenotypes, and (2) to assess the effect heterogeneity of risk factors across component phenotypes of LVOT defects. ^ This study was a secondary data analysis. Primary data were obtained from the Texas Birth Defect Registry (TBDR). TBDR uses an active surveillance method to ascertain birth defects in Texas. All cases of non complex LVOT defects which met our inclusion criteria during the period of 2002–2008 were included in the study. The comparison groups included all unaffected live births for the same period (2002–2008). Data from vital statistics were used to evaluate associations. Statistical associations between selected risk factors and LVOT defects was determined by calculating crude and adjusted prevalence ratio using Poisson regression analysis. Effect heterogeneity was evaluated using polytomous logistic regression. ^ There were a total of 2,353 cases of LVOT defects among 2,730,035 live births during the study period. There were a total of 1,311 definite cases of non-complex LVOT defects for analysis after excluding "complex" cardiac cases and cases associated with syndromes (n=168). Among infant characteristics, males were at a significantly higher risk of developing LVOT defects compared to females. Among maternal characteristics, significant associations were seen with maternal age > 40 years (compared to maternal age 20–24 years) and maternal residence in Texas-Mexico border (compared to non-border residence). Among birth characteristics, significant associations were seen with preterm birth and small for gestation age LVOT defects. ^ When evaluating effect heterogeneity, the following variables had significantly different effects among the component LVOT defect phenotypes: infant sex, plurality, maternal age, maternal race/ethnicity, and Texas-Mexico border residence. ^ This study found significant associations between various demographic factors and LVOT defects. While many findings from this study were consistent with results from previous studies, we also identified new factors associated with LVOT defects. Additionally, this study was the first to assess effect heterogeneity across LVOT defect component phenotypes. These findings contribute to a growing body of literature on characteristics associated with LVOT defects. ^

Gene-gene interaction and gene-pathway analysis of genome-wide association for schizophrenia

Schizophrenia (SZ) is a complex disorder with high heritability and variable phenotypes that has limited success in finding causal genes associated with the disease development. Pathway-based analysis is an effective approach in investigating the molecular mechanism of susceptible genes associated with complex diseases. The etiology of complex diseases could be a network of genetic factors and within the genes, interaction may occur. In this work we argue that some genes might be of small effect that by itself are neither sufficient nor necessary to cause the disease however, their effect may induce slight changes to the gene expression or affect the protein function, therefore, analyzing the gene-gene interaction mechanism within the disease pathway would play crucial role in dissecting the genetic architecture of complex diseases, making the pathway-based analysis a complementary approach to GWAS technique. ^ In this study, we implemented three novel linkage disequilibrium based statistics, the linear combination, the quadratic, and the decorrelation test statistics, to investigate the interaction between linked and unlinked genes in two independent case-control GWAS datasets for SZ including participants of European (EA) and African (AA) ancestries. The EA population included 1,173 cases and 1,378 controls with 729,454 genotyped SNPs, while the AA population included 219 cases and 288 controls with 845,814 genotyped SNPs. We identified 17,186 interacting gene-sets at significant level in EA dataset, and 12,691 gene-sets in AA dataset using the gene-gene interaction method. We also identified 18,846 genes in EA dataset and 19,431 genes in AA dataset that were in the disease pathways. However, few genes were reported of significant association to SZ. ^ Our research determined the pathways characteristics for schizophrenia through the gene-gene interaction and gene-pathway based approaches. Our findings suggest insightful inferences of our methods in studying the molecular mechanisms of common complex diseases.^

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.

Patenting genes and gene sequences: Analysis of public health implications and the law

At issue is whether or not isolated DNA is patent eligible under the U.S. Patent Law and the implications of that determination on public health. The U.S. Patent and Trademark Office has issued patents on DNA since the 1980s, and scientists and researchers have proceeded under that milieu since that time. Today, genetic research and testing related to the human breast cancer genes BRCA1 and BRCA2 is conducted within the framework of seven patents that were issued to Myriad Genetics and the University of Utah Research Foundation between 1997 and 2000. In 2009, suit was filed on behalf of multiple researchers, professional associations and others to invalidate fifteen of the claims underlying those patents. The Court of Appeals for the Federal Circuit, which hears patent cases, has invalidated claims for analyzing and comparing isolated DNA but has upheld claims to isolated DNA. The specific issue of whether isolated DNA is patent eligible is now before the Supreme Court, which is expected to decide the case by year's end. In this work, a systematic review was performed to determine the effects of DNA patents on various stakeholders and, ultimately, on public health; and to provide a legal analysis of the patent eligibility of isolated DNA and the likely outcome of the Supreme Court's decision. ^ A literature review was conducted to: first, identify principle stakeholders with an interest in patent eligibility of the isolated DNA sequences BRCA1 and BRCA2; and second, determine the effect of the case on those stakeholders. Published reports that addressed gene patents, the Myriad litigation, and implications of gene patents on stakeholders were included. Next, an in-depth legal analysis of the patent eligibility of isolated DNA and methods for analyzing it was performed pursuant to accepted methods of legal research and analysis based on legal briefs, federal law and jurisprudence, scholarly works and standard practice legal analysis. ^ Biotechnology, biomedical and clinical research, access to health care, and personalized medicine were identified as the principle stakeholders and interests herein. Many experts believe that the patent eligibility of isolated DNA will not greatly affect the biotechnology industry insofar as genetic testing is concerned; unlike for therapeutics, genetic testing does not require tremendous resources or lead time. The actual impact on biomedical researchers is uncertain, with greater impact expected for researchers whose work is intended for commercial purposes (versus basic science). The impact on access to health care has been surprisingly difficult to assess; while invalidating gene patents might be expected to decrease the cost of genetic testing and improve access to more laboratories and physicians' offices that provide the test, a 2010 study on the actual impact was inconclusive. As for personalized medicine, many experts believe that the availability of personalized medicine is ultimately a public policy issue for Congress, not the courts. ^ Based on the legal analysis performed in this work, this writer believes the Supreme Court is likely to invalidate patents on isolated DNA whose sequences are found in nature, because these gene sequences are a basic tool of scientific and technologic work and patents on isolated DNA would unduly inhibit their future use. Patents on complementary DNA (cDNA) are expected to stand, however, based on the human intervention required to craft cDNA and the product's distinction from the DNA found in nature. ^ In the end, the solution as to how to address gene patents may lie not in jurisprudence but in a fundamental change in business practices to provide expanded licenses to better address the interests of the several stakeholders. ^

Analysis of the mechanism of replication inhibition of the tumor suppressor gene p53

p53 plays a role in cell cycle arrest and apoptosis. p53 has also been shown to be involved in DNA replication. To study the effect of p53 on DNA replication, we utilized a SV40 based shuttle vector system. The pZ402 shuttle vector, was constructed with a mutated T-antigen unable to interact with p53 but able to support replication of the shuttle vector. When a transcriptional activation domain p53 mutant was tested for its ability to inhibit DNA replication no inhibition was observed. Competition assays with the DNA binding domain of p53 was also able to block the inhibition of DNA replication by p53 suggesting that p53 can inhibit DNA replication through the transcriptional activation of a target gene. One likely target gene, p21$\sp{\rm cip/waf}$ was tested to determine whether p53 inhibited DNA replication by transcriptionally activating p21$\sp{\rm cip/waf}$. Two independent approaches utilizing p21$\sp{\rm cip/waf}$ null cells or the expression of an anti-sense p21$\sp{\rm cip/waf}$ expression vector were utilized. p53 was able to inhibit pZ402 replication independently of p21$\sp{\rm cip/waf}$. p53 was also able to inhibit DNA replication independent of the p53 target genes Gadd45 and the replication processivity factor PCNA. The inhibition of DNA replication by p53 was also independent of direct DNA binding to a consensus site on the replicating plasmid. p53 mutants can be classified into two categories: conformational and DNA contact mutants. The two types of p53 mutants were tested for their effects on DNA replication. While all conformational mutants were unable to inhibit DNA replication three out of three DNA contact mutants tested were able to inhibit DNA replication. The work here studies the effect wild-type and mutant p53 has on DNA replication and demonstrated a possible mechanism by which wild-type p53 could inhibit DNA replication through the transcriptional activation of a target gene. ^