Bisulfite genomic sequencing: systematic investigation of critical experimental parameters

Bisulfite genomic sequencing is the method of choice for the generation of methylation maps with single-base resolution. The method is based on the selective deamination of cytosine to uracil by treatment with bisulfite and the sequencing of subsequently generated PCR products. In contrast to cytosine, 5-methylcytosine does not react with bisulfite and can therefore be distinguished. In order to investigate the potential for optimization of the method and to determine the critical experimental parameters, we determined the influence of incubation time and incubation temperature on the deamination efficiency and measured the degree of DNA degradation during the bisulfite treatment. We found that maximum conversion rates of cytosine occurred at 55 degrees C (4-18 h) and 95 degrees C (1 h). Under these conditions at least 84-96% of the DNA is degraded. To study the impact of primer selection, homologous DNA templates were constructed possessing cytosine-containing and cytosine-free primer binding sites, respectively. The recognition rates for cytosine (>/=97%) and 5-methylcytosine (>/=94%) were found to be identical for both templates.

Identification and localization of a new human myotubularin-related protein gene, mtmr8, on 8p22-p23

Myotubularin and myotubularin-related proteins are dual-specificity phosphatases. Several myotubularin-related proteins have been identified in humans and mice. The members of the myotubularin protein family are highly conserved, from humans to yeast. Mutations in the human myotubularin gene (MTM1) lead to X-linked myotubular myopathy. Here we isolate and localize a novel putative myotubularin-related protein gene (MTMR8) on chromosome 8p22--p23,between the markers D8S550 and D8S265, by exon-trapping experiments and RT-PCR. Genomic sequencing revealed that the gene consists of 10 exons and spans approximately 43 kb. The corresponding cDNA is 7081 bp. The open reading frame predicts a protein of 549 amino acids and a calculated molecular mass of 63 kDa. Like myotubularin-related protein-5, MTMR8 has no dual-specificity phosphatase domain. It contains a double-helical motif similar to the SET interaction domain, which is thought to have a role in the control of cell proliferation.

Comparison of the ocular hypotensive lipid AGN 192024 with timolol: dosing, efficacy, and safety evaluation of a novel compound for glaucoma management

OBJECTIVE

To compare the safety and efficacy of the ocular hypotensive lipid AGN 192024 (Lumigan) with those of timolol.

METHODS

A 30-day, randomized, investigator-masked, clinical trial involving 100 patients with elevated intraocular pressure (IOP). Study medications were instilled topically. Doses of 0.003%, 0.01%, or 0.03% AGN 192024 were given once daily for 3 weeks then twice daily for 1 week, and vehicle control or 0.5% timolol was given twice daily for 4 weeks. Mean change in IOP from baseline was the primary efficacy variable. Safety parameters included adverse events, hyperemia grading, laser flare meter analysis, heart rate, and blood pressure.

RESULTS

Timolol and all 3 concentrations of AGN 192024 lowered IOP from baseline (P < .001). A dosage of 0.03% AGN 192024 once daily lowered IOP significantly more than timolol (P < or = .02) at every study visit except day 21 (P = .053) and provided better diurnal IOP control. Twice-daily dosing of AGN 192024 provided no clinically significant benefit over once-daily dosing. All treatment regimens were safe and well tolerated, with no clinically significant effects on heart rate or blood pressure and no between-group differences in the incidence of adverse events. The only significant ocular safety finding with AGN 192024 was a dose-related mild increase in conjunctival hyperemia.

CONCLUSIONS

Of the 3 concentrations tested, 0.03% AGN 192024 once daily had the best therapeutic profile. AGN 192024 was safe and well tolerated, and it provided superior ocular hypotensive efficacy and diurnal IOP control compared with timolol in patients with ocular hypertension and glaucoma.

Human-mouse comparative sequence analysis of the NEMO gene reveals an alternative promoter within the neighboring G6PD gene

NEMO (NFkappaB essential modulator) is a non-catalytic subunit of the cytokine-dependent IkappaB kinase complex that is involved in activation of the transcription factor NFkappaB. The human NEMO gene maps to Xq28 and is arranged head to head with the proximal G6PD gene. Mutations in NEMO have recently been associated with Incontinentia Pigmenti (Smahi et al., Nature 405 (2000) 466), an X-linked dominant disorder. Three alternative transcripts with different non-coding 5' exons (1a, 1b and 1c) of NEMO have been described. In order to identify regulatory elements that control alternative transcription we have established the complete genomic sequence of the murine orthologs Nemo and G6pdx. Sequence comparison suggests the presence of two alternative promoters for NEMO/Nemo. First, a CpG island is shared by both genes driving expression of the NEMO/Nemo transcripts containing exons 1b and 1c in one direction and the housekeeping gene G6PD/G6pdx in the opposite direction. In contrast to human, an additional variant of exon 1c, named 1c+, was identified in several tissues of the mouse. This larger exon utilizes an alternative donor site located 1594 bp within intron 1c. The putative second promoter for NEMO/Nemo transcripts starting with exon 1a is unidirectional, and not associated with a CpG island. Surprisingly, this promoter is located in the second intron of G6PD/G6pdx. It shows very low basal activity and may be involved in stress/time- and/or tissue-dependent expression of NEMO. To our knowledge, an overlapping gene order similar to the G6PD/NEMO complex has not been described before.

Combined effect of operational variables and enzyme activity on aqueous enzymatic extraction of oil and protein from soybean

The individual effect of two different enzymes-protease and cellulase-on oil and protein extraction yields combined with other process parameters-enzyme concentration, time of hydrolysis, particle size and solid-to-liquid ratio-was evaluated by Response Surface Methodology. The selection of the enzymes for the study was based on preliminary experiments that showed higher increments in the extraction yield with the use of the two enzymes when compared to hemicellulase and pectinase. The levels of the quantitative parameters studied were: i) enzyme concentration: 0.1, 0.45, 2 w/w %; ii) liquid-to-solid ratio: 0.05, 0.125, 0.2; iii) mean particle size: 212.5, 449.5, 855 µm; iv) time of hydrolysis: 30; 60; 120 min. Experimental data for both oil and protein extraction yields obtained with and without enzymes correlated very well with process parameters (P < 0.0001), resulting in models with high coefficient of determination for oil and protein extraction yields (r(2) = 0.9570 and r(2) = 0.9807, respectively). The use of protease resulted in significantly higher yields over the control (protein yield increased from 27.8 to 66.2%, oil yield increased from 41.8 to 58.7%) only when heat treated flour was used, or when non-heat treated flour with large particle sizes was used in the extraction. The yields of protein and oil from non-heat treated material in general decreased slightly with the use of enzymes.

Elevated DNA sequence diversity in the genomic region of the phosphatase PPP2R3L gene in the human pseudoautosomal region

The evolution, inheritance and recombination rate of genes located in the pseudoautosomal region 1 (PAR1) is exceptional within the human genome. Pseudoautosomal genes are identical on X and Y chromosomes and are not inherited in a sex linked manner. Due to an obligatory recombination event in male meiosis, pseudoautosomal genes are exchanged frequently between X and Y chromosomes. During the isolation, characterization and sequencing of a novel gene PPP2R3L, which was classified by sequence homology as a novel member of the protein phosphatase regulatory subunit families, it became apparent that cosmids of different origin harboring this gene are highly polymorphic between individuals, both at the nucleotide level and in the number.

The complex repeats of Dictyostelium discoideum

In the course of determining the sequence of the Dictyostelium discoideum genome we have characterized in detail the quantity and nature of interspersed repetitive elements present in this species. Several of the most abundant small complex repeats and transposons (DIRS-1; TRE3-A,B; TRE5-A; skipper; Tdd-4; H3R) have been described previously. In our analysis we have identified additional elements. Thus, we can now present a complete list of complex repetitive elements in D. discoideum. All elements add up to 10% of the genome. Some of the newly described elements belong to established classes (TRE3-C, D; TRE5-B,C; DGLT-A,P; Tdd-5). However, we have also defined two new classes of DNA transposable elements (DDT and thug) that have not been described thus far. Based on the nucleotide amount, we calculated the least copy number in each family. These vary between <10 up to >200 copies. Unique sequences adjacent to the element ends and truncation points in elements gave a measure for the fragmentation of the elements. Furthermore, we describe the diversity of single elements with regard to polymorphisms and conserved structures. All elements show insertion preference into loci in which other elements of the same family reside. The analysis of the complex repeats is a valuable data resource for the ongoing assembly of whole D. discoideum chromosomes.

Complex Arrangement of Genes within a 220-kb Region of Double-Duplicated DNA on Human 2q37.1

Gene duplication events are followed by divergence of initially identical gene copies, due to the subsequent accumulation of mutations. These mutations tend to be degenerative and may lead to either nonfunctionalization or subfunctionalization of the gene copies. Here we report the molecular characterization of a 220-kb genomic DNA fragment from human 2q37.1, in which a double duplication and a partial triplication event has taken place. As a result, this region contains four copies of alkaline phosphatase (P), four copies of the ECEL1 gene (X), two copies of a newly identified gene (N), and two copies of a cholinergic receptor subunit (R), in the order N-P-X-P-X-P-X-N-P-X-R-R. While three of the four ECEL1 copies, one copy of the phosphatase gene and one copy of the newly identified gene have lost their function, three phosphatase gene copies and the two receptor subunits are still functionally active and thus may provide an example for subfunctionalization of duplicated genes.

Molecular cloning of Xp11 breakpoints in two unrelated mentally retarded females with X;autosome translocations

Mental retardation is a very common and extremely heterogeneous disorder that affects about 3% of the human population. Its molecular basis is largely unknown, but many loci have been mapped to the X chromosome. We report on two mentally retarded females with X;autosome translocations and breakpoints in Xp11, viz., t(X;17)(p11;p13) and t(X;20)(p11;q13). (Fiber-) FISH analysis assigned the breakpoints to different subbands, Xp11.4 and Xp11.23, separated by approximately 8 Mb. High-resolution mapping of the X- chromosome breakpoints using Southern blot hybridization resulted in the isolation of breakpoint-spanning genomic subclones of 3 kb and 0. 5 kb. The Xp11.4 breakpoint is contained within a single copy sequence, whereas the Xp11.23 breakpoint sequence resembles an L1 repetitive element. Several expressed sequences map close to the breakpoints, but none was found to be inactivated. Therefore, mechanisms other than disruption of X-chromosome genes likely cause the phenotypes.

Initial sequencing and analysis of the human genome

The human genome holds an extraordinary trove of information about human development, physiology, medicine and evolution. Here we report the results of an international collaboration to produce and make freely available a draft sequence of the human genome. We also present an initial analysis of the data, describing some of the insights that can be gleaned from the sequence.

A physical map of the human genome

The human genome is by far the largest genome to be sequenced, and its size and complexity present many challenges for sequence assembly. The International Human Genome Sequencing Consortium constructed a map of the whole genome to enable the selection of clones for sequencing and for the accurate assembly of the genome sequence. Here we report the construction of the whole-genome bacterial artificial chromosome (BAC) map and its integration with previous landmark maps and information from mapping efforts focused on specific chromosomal regions. We also describe the integration of sequence data with the map.

The DAPIN family: a novel domain links apoptotic and interferon response proteins

We report the discovery of a protein domain, hereafter referred to as DAPIN, in diverse vertebrate and viral proteins that is associated with tumor biology, apoptosis and inflammation. Based on a secondary structure prediction, we suggest an all-alpha fold for DAPIN, which is also adopted by apoptotic protein domains of the CARD, death domain and death effector domain type.

Topical ketorolac tromethamine 0.5% ophthalmic solution in ocular inflammation after cataract surgery

PURPOSE

To compare the efficacy and safety of ketorolac 0.5% ophthalmic solution with its vehicle in the treatment of ocular inflammation after cataract surgery and intraocular lens implantation.

DESIGN

Multicenter clinical study.

PARTICIPANTS

One hundred four patients were prospectively randomized, 52 patients in treatment group, 52 patients in control group.

METHODS

Patients received either ketorolac or vehicle four times daily in the operated eye for 14 days starting the day after surgery in a prospective, double-masked, randomized, parallel group study. Only patients with moderate or greater postoperative inflammation the day after surgery were enrolled.

MAIN OUTCOME MEASURES

The main outcome measures include inflammation (cell, flare, ciliary flush), intraocular pressure and visual acuity.

RESULTS

Ketorolac was significantly more effective than vehicle in reducing the manifestations of postoperative ocular inflammation, including: anterior chamber cells (P: = 0.002) and flare (P: = 0.009), conjunctival erythema (P: = 0.010), ciliary flush (P: = 0.022), photophobia (P: = 0.027), and pain (P: = 0.043). Five times as many patients were dropped from the study for lack of efficacy from the vehicle group (22/52) than from the ketorolac group (4/52; P: = 0.001). Ketorolac was found to be equally as safe as vehicle in terms of adverse events, changes in visual acuity, intraocular pressure, and biomicroscopic and ophthalmoscopic variables.

CONCLUSIONS

Ketorolac tromethamine 0.5% ophthalmic solution was significantly more effective than vehicle in the treatment of moderate or greater ocular inflammation following routine cataract surgery, while being as safe as vehicle.

Genome sequence of bovine herpesvirus 4, a bovine Rhadinovirus, and identification of an origin of DNA replication

Bovine herpesvirus 4 (BoHV-4) is a gammaherpesvirus of cattle. The complete long unique coding region (LUR) of BoHV-4 strain 66-p-347 was determined by a shotgun approach. Together with the previously published noncoding terminal repeats, the entire genome sequence of BoHV-4 is now available. The LUR consists of 108,873 bp with an overall G+C content of 41.4%. At least 79 open reading frames (ORFs) are present in this coding region, 17 of them unique to BoHV-4. In contrast to herpesvirus saimiri and human herpesvirus 8, BoHV-4 has a reduced set of ORFs homologous to cellular genes. Gene arrangement as well as phylogenetic analysis confirmed that BoHV-4 is a member of the genus Rhadinovirus. In addition, an origin of replication (ori) in the genome of BoHV-4 was identified by DpnI assays. A minimum of 1.69 kbp located between ORFs 69 and 71 was sufficient to act as a cis signal for replication.

Discovery of new human beta-defensins using a genomics-based approach

Epithelial beta-defensins are broad-spectrum cationic antimicrobial peptides that also act as chemokines for adaptive immune cells. In the human genome, all known defensin genes cluster to a <1 Mb region of chromosome 8p22-p23. To identify new defensin genes, the DNA sequence from a contig of large-insert genomic clones from the region containing human beta-defensin-2 (HBD-2) was analyzed for the presence of defensin genes. This sequence survey identified a novel beta-defensin, termed HBD-3. The HBD-3 gene contains two exons, is located 13 kb upstream from the HBD-2 gene, and it is transcribed in the same direction. A partial HBD-3 cDNA clone was amplified from cDNA derived from IL-1beta induced fetal lung tissue. The cDNA sequence encodes for a 67 amino acid peptide that is approximately 43% identical to HBD-2 and shares the beta-defensin six cysteine motif. By PCR analysis of two commercial cDNA panels, HBD-3 expression was detected in adult heart, skeletal muscle, placenta and in fetal thymus. From RT-PCR experiments, HBD-3 expression was observed in skin, esophagus, gingival keratinocytes, placenta and trachea. Furthermore, in fetal lung explants and gingival keratinocytes, HBD-3 mRNA expression was induced by IL-1beta. Additional sequence analysis identified the HE2 (human epididymis secretory protein) gene 17 kb upstream from the HBD-3 gene. One splice variant of this gene (HE2beta1) encodes a beta-defensin consensus cysteine motif, suggesting it represents a defensin gene product. HE2beta1 mRNA expression was detected in gingival keratinocytes and bronchial epithelia using RT-PCR analysis. The discovery of these novel beta-defensin genes may allow further understanding of the role of defensins in host immunity at mucosal surfaces.

Architecture and anatomy of the genomic locus encoding the human leukemia-associated transcription factor RUNX1/AML1

The RUNX1 gene on human chromosome 21q22.12 belongs to the 'runt domain' gene family of transcription factors (also known as AML/CBFA/PEBP2alpha). RUNX1 is a key regulator of hematopoiesis and a frequent target of leukemia associated chromosomal translocations. Here we present a detailed analysis of the RUNX1 locus based on its complete genomic sequence. RUNX1 spans 260 kb and its expression is regulated through two distinct promoter regions, that are 160 kb apart. A very large CpG island complex marks the proximal promoter (promoter-2), and an additional CpG island is located at the 3' end of the gene. Hitherto, 12 different alternatively spliced RUNX1 cDNAs have been identified. Genomic sequence analysis of intron/exon boundaries of these cDNAs has shown that all consist of properly spliced authentic coding regions. This indicates that the large repertoire of RUNX1 proteins, ranging in size between 20-52 kDa, are generated through usage of alternatively spliced exons some of which contain in frame stop codons. The gene's introns are largely depleted of repetitive sequences, especially of the LINE1 family. The RUNX1 locus marks the transition from a ~1 Mb of gene-poor region containing only pseudogenes, to a gene-rich region containing several functional genes. A search for RUNX1 sequences that may be involved in the high frequency of chromosomal translocations revealed that a 555 bp long segment originating in chromosome 11 FLI1 gene was transposed into RUNX1 intron 4.1. This intron harbors the t(8;21) and t(3;21) chromosomal breakpoints involved in acute myeloid leukemia. Interestingly, the FLI1 homologous sequence contains a breakpoint of the t(11;22) translocation associated with Ewing's tumors, and may have a similar function in RUNX1.

MethDB--a public database for DNA methylation data

Methylation of cytosine in the 5 position of the pyrimidine ring is a major modification of the DNA in most organisms. In eukaryotes, the distribution and number of 5-methylcytosines (5mC) along the DNA is heritable but can also change with the developmental state of the cell and as a response to modifications of the environment. While DNA methylation probably has a number of functions, scientific interest has recently focused on the gene silencing effect methylation can have in eukaryotic cells. In particular, the discovery of changes in the methylation level during cancer development has increased the interest in this field. In the past, a vast amount of data has been generated with different levels of resolution ranging from 5mC content of total DNA to the methylation status of single nucleotides. We present here a database for DNA methylation data that attempts to unify these results in a common resource. The database is accessible via WWW (http://www.methdb.de). It stores information about the origin of the investigated sample and the experimental procedure, and contains the DNA methylation data. Query masks allow for searching for 5mC content, species, tissue, gene, sex, phenotype, sequence ID and DNA type. The output lists all available information including the relative gene expression level. DNA methylation patterns and methylation profiles are shown both as a graphical representation and as G/A/T/C/5mC-sequences or tables with sequence positions and methylation levels, respectively.

Downregulation of Hedgehog signaling is required for organogenesis of the small intestine in Xenopus

Hedgehog ligands interact with receptor complexes containing Patched (PTC) and Smoothened (SMO) proteins to regulate many aspects of development. The mutation W535L (SmoM2) in human Smo is associated with basal cell skin cancers, causes constitutive, ligand-independent signaling through the Hedgehog pathway, and provides a powerful means to test effects of unregulated Hedgehog signaling. Expression of SmoM2 in Xenopus embryos leads to developmental anomalies that are consistent with known requirements for regulated Hedgehog signaling in the eye and pancreas. Additionally, it results in failure of midgut epithelial cytodifferentiation and of the intestine to lengthen and coil. The midgut mesenchyme shows increased cell numbers and attenuated expression of the differentiation marker smooth muscle actin. With the exception of the pancreas, differentiation of foregut and hindgut derivatives is unaffected. The intestinal epithelial abnormalities are reproduced in embryos or organ explants treated directly with active recombinant hedgehog protein. Ptc mRNA, a principal target of Hedgehog signaling, is maximally expressed at stages corresponding to the onset of the intestinal defects. In advanced embryos expressing SmoM2, Ptc expression is remarkably confined to the intestinal wall. Considered together, these findings suggest that the splanchnic mesoderm responds to endodermal Hedgehog signals by inhibiting the transition of midgut endoderm into intestinal epithelium and that attenuation of this feedback is required for normal development of the vertebrate intestine.

Drinks like a fish: zebra fish (Danio rerio) as a behavior genetic model to study alcohol effects

Zebra fish may be an ideal vertebrate model system for numerous human diseases with which the genetics and biological mechanisms of the disease may be studied. Zebra fish has been successfully used in developmental genetics, and recently, neurobiologists have also started to study this species. A potentially interesting target disease amenable for analysis with zebra fish is drug addiction, e.g. alcoholism. Although genetic tools to manipulate the genome of zebra fish are available, appropriate phenotypical testing methods are often lacking. In this paper, we describe basic behavioral tests to investigate the acute effects of alcohol on zebra fish. These behavioral paradigms will be useful for the genetic and biological analysis of acute and chronic drug effects as well as addiction. In addition to presenting findings for the acute effects of alcohol, we briefly describe our strategy for generating and screening mutants. We hope that our pilot work will facilitate the future development of behavioral tests and the use of zebra fish in the genetic analysis of the biological effects of drugs of abuse.

A regulator of transcriptional elongation controls vertebrate neuronal development

The development of distinct vertebrate neurons is defined by the unique profiles of genes that neurons express. It is accepted that neural genes are regulated at the point of transcription initiation, but the role of messenger RNA elongation in neural gene regulation has not been examined. Here we describe the mutant foggy, identified in a genetic screen for mutations that affect neuronal development in zebrafish, that displayed a reduction of dopamine-containing neurons and a corresponding surplus of serotonin-containing neurons in the hypothalamus. Positional cloning disclosed that Foggy is a brain-enriched nuclear protein that is structurally related to the transcription elongation factor Spt5 (refs 5-12). Foggy is not part of the basic transcription apparatus but a phosphorylation-dependent, dual regulator of transcription elongation. The mutation disrupts its repressive but not its stimulatory activity. Our results provide molecular, genetic and biochemical evidence that negative regulators of transcription elongation control key aspects of neuronal development.

Large-scale methylation analysis of human genomic DNA reveals tissue-specific differences between the methylation profiles of genes and pseudogenes

Cytosine in CpG dinucleotides is frequently found to be methylated in the DNA of higher eukaryotes and differential methylation has been proposed to be a key element in the organization of gene expression in man. To address this question systematically, we used bisulfite genomic sequencing to study the methylation patterns of three X-linked genes and one autosomal pseudogene in two adult individuals and across nine different tissues. Two of the genes, SLC6A8 and MSSK1, are tissue-specifically expressed. CDM is expressed ubiquitously. The pseudogene, psi SLC6A8, is exclusively expressed in the testis. The promoter regions of the SLC6A8, MSSK1 and CDM genes were found to be essentially unmethylated in all tissues, regardless of their relative expression level. In contrast, the pseudogene psi SLC6A8 shows high methylation of the CpG islands in all somatic tissues but complete demethylation in testis. Methylation profiles in different tissues are similar in shape but not identical. The data for the two investigated individuals suggest that methylation profiles of individual genes are tissue specific. Taken together, our findings support a model in which the bodies of the genes are predominantly methylated and thus insulated from the interaction with DNA-binding proteins. Only unmethylated promoter regions are accessible for binding and interaction. Based on this model we propose to use DNA methylation studies in conjunction with large-scale sequencing approaches as a tool for the prediction of cis-acting genomic regions, for the identification of cryptic and potentially active CpG islands and for the preliminary distinction of genes and pseudogenes.