Data-driven analysis approach for biomarker discovery using molecular-profiling technologies.

High-throughput molecular-profiling technologies provide rapid, efficient and systematic approaches to search for biomarkers. Supervised learning algorithms are naturally suited to analyse a large amount of data generated using these technologies in biomarker discovery efforts. The study demonstrates with two examples a data-driven analysis approach to analysis of large complicated datasets collected in high-throughput technologies in the context of biomarker discovery. The approach consists of two analytic steps: an initial unsupervised analysis to obtain accurate knowledge about sample clustering, followed by a second supervised analysis to identify a small set of putative biomarkers for further experimental characterization. By comparing the most widely applied clustering algorithms using a leukaemia DNA microarray dataset, it was established that principal component analysis-assisted projections of samples from a high-dimensional molecular feature space into a few low dimensional subspaces provides a more effective and accurate way to explore visually and identify data structures that confirm intended experimental effects based on expected group membership. A supervised analysis method, shrunken centroid algorithm, was chosen to take knowledge of sample clustering gained or confirmed by the first step of the analysis to identify a small set of molecules as candidate biomarkers for further experimentation. The approach was applied to two molecular-profiling studies. In the first study, PCA-assisted analysis of DNA microarray data revealed that discrete data structures exist in rat liver gene expression and correlated with blood clinical chemistry and liver pathological damage in response to a chemical toxicant diethylhexylphthalate, a peroxisome-proliferator-activator receptor agonist. Sixteen genes were then identified by shrunken centroid algorithm as the best candidate biomarkers for liver damage. Functional annotations of these genes revealed roles in acute phase response, lipid and fatty acid metabolism and they are functionally relevant to the observed toxicities. In the second study, 26 urine ions identified from a GC/MS spectrum, two of which were glucose fragment ions included as positive controls, showed robust changes with the development of diabetes in Zucker diabetic fatty rats. Further experiments are needed to define their chemical identities and establish functional relevancy to disease development.

Differential effects of estrogen and raloxifene on messenger RNA and matrix metalloproteinase 2 activity in the rat uterus.

A detailed analysis of the differential effects of estrogen (E) compared to raloxifene (Ral), a selective estrogen receptor modulator (SERM), following estrogen receptor (ER) binding in gynecological tissues was conducted using gene microarrays, Northern blot analysis, and matrix metalloproteinase (MMP) 2 activity studies. We profiled gene expression in the uterus following acute (1 day) and prolonged daily (5 wk) treatment of E and Ral in ovariectomized rats. Estrogen regulated twice as many genes as Ral, largely those associated with catalysis and metabolism, whereas Ral induced genes associated with cell death and negative cell regulation. Follow-up studies confirmed that genes associated with matrix integrity were differentially regulated by Ral and E at various time points in uterine and vaginal tissues. Additional experiments were conducted to determine the levels of MMP2 activity in uterus explants from ovariectomized rats following 2 wk of treatment with E, Ral, or one of two additional SERMs: lasofoxifene, and levormeloxifene. Both E and lasofoxifene stimulated uterine MMP2 activity to a level twofold that of Ral, whereas levormeloxifene elevated MMP2 activity to a level 12-fold that of Ral. These data show that one of the significant differences between E and Ral signaling in the uterus is the regulation of genes and proteins associated with matrix integrity. This may be a potential key difference between the action of SERMs in the uterus of postmenopausal women.

Promoter recognition and discrimination by EsigmaS RNA polymerase.

Although more than 30 Escherichia coli promoters utilize the RNA polymerase holoenzyme containing sigmaS (EsigmaS), and it is known that there is some overlap between the promoters recognized by EsigmaS and by the major E. coli holoenzyme (Esigma70), the sequence elements responsible for promoter recognition by EsigmaS are not well understood. To define the DNA sequences recognized best by EsigmaS in vitro, we started with random DNA and enriched for EsigmaS promoter sequences by multiple cycles of binding and selection. Surprisingly, the sequences selected by EsigmaS contained the known consensus elements (-10 and -35 hexamers) for recognition by Esigma70. Using genetic and biochemical approaches, we show that EsigmaS and Esigma70 do not achieve specificity through 'best fit' to different consensus promoter hexamers, the way that other forms of holoenzyme limit transcription to discrete sets of promoters. Rather, we suggest that EsigmaS-specific promoters have sequences that differ significantly from the consensus in at least one of the recognition hexamers, and that promoter discrimination against Esigma70 is achieved, at least in part, by the two enzymes tolerating different deviations from consensus. DNA recognition by EsigmaS versus Esigma70 thus presents an alternative solution to the problem of promoter selectivity.

Genome of the bacterium Streptococcus pneumoniae strain R6.

Streptococcus pneumoniae is among the most significant causes of bacterial disease in humans. Here we report the 2,038,615-bp genomic sequence of the gram-positive bacterium S. pneumoniae R6. Because the R6 strain is avirulent and, more importantly, because it is readily transformed with DNA from homologous species and many heterologous species, it is the principal platform for investigation of the biology of this important pathogen. It is also used as a primary vehicle for genomics-based development of antibiotics for gram-positive bacteria. In our analysis of the genome, we identified a large number of new uncharacterized genes predicted to encode proteins that either reside on the surface of the cell or are secreted. Among those proteins there may be new targets for vaccine and antibiotic development.

Bacterial promoter architecture: subsite structure of UP elements and interactions with the carboxy-terminal domain of the RNA polymerase alpha subunit.

We demonstrate here that the previously described bacterial promoter upstream element (UP element) consists of two distinct subsites, each of which, by itself, can bind the RNA polymerase holoenzyme alpha subunit carboxy-terminal domain (RNAP alphaCTD) and stimulate transcription. Using binding-site-selection experiments, we identify the consensus sequence for each subsite. The selected proximal subsites (positions -46 to -38; consensus 5'-AAAAAARNR-3') stimulate transcription up to 170-fold, and the selected distal subsites (positions -57 to -47; consensus 5'-AWWWWWTTTTT-3') stimulate transcription up to 16-fold. RNAP has subunit composition alpha(2)betabeta'sigma and thus contains two copies of alphaCTD. Experiments with RNAP derivatives containing only one copy of alphaCTD indicate, in contrast to a previous report, that the two alphaCTDs function interchangeably with respect to UP element recognition. Furthermore, function of the consensus proximal subsite requires only one copy of alphaCTD, whereas function of the consensus distal subsite requires both copies of alphaCTD. We propose that each subsite constitutes a binding site for a copy of alphaCTD, and that binding of an alphaCTD to the proximal subsite region (through specific interactions with a consensus proximal subsite or through nonspecific interactions with a nonconsensus proximal subsite) is a prerequisite for binding of the other alphaCTD to the distal subsite.

Identification of an UP element consensus sequence for bacterial promoters.

The UP element, a component of bacterial promoters located upstream of the -35 hexamer, increases transcription by interacting with the RNA polymerase alpha-subunit. By using a modification of the SELEX procedure for identification of protein-binding sites, we selected in vitro and subsequently screened in vivo for sequences that greatly increased promoter activity when situated upstream of the Escherichia coli rrnB P1 core promoter. A set of 31 of these upstream sequences increased transcription from 136- to 326-fold in vivo, considerably more than the natural rrnB P1 UP element, and was used to derive a consensus sequence: -59 nnAAA(A/T)(A/T)T(A/T)TTTTnnAAAAnnn -38. The most active selected sequence contained the derived consensus, displayed all of the properties of an UP element, and the interaction of this sequence with the alpha C-terminal domain was similar to that of previously characterized UP elements. The identification of the UP element consensus should facilitate a detailed understanding of the alpha-DNA interaction. Based on the evolutionary conservation of the residues in alpha responsible for interaction with UP elements, we suggest that the UP element consensus sequence should be applicable throughout eubacteria, should generally facilitate promoter prediction, and may be of use for biotechnological applications.

The chicken malic enzyme gene: structural organization and identification of triiodothyronine response elements in the 5'-flanking DNA.

In vivo, feeding stimulates and starvation inhibits transcription of the malic enzyme gene. In chick-embryo hepatocytes in culture, triiodothyronine (T3) stimulates and glucagon inhibits transcription of this gene. As a first step in the characterization of the involved regulatory mechanisms, fragments of genomic DNA spanning the structural and 5'-flanking regions of the chicken malic enzyme gene were cloned. The coding region of the gene is organized into 14 exons and 13 introns and is greater than 106 kb in length. The size of the gene, the number and lengths of the exons, and positions at which introns are inserted into the coding regions are virtually identical in the chicken and rat genes. When transiently transfected into chick-embryo hepatocytes, 5800 bp of 5'-flanking DNA conferred T3 responsiveness to a linked chloramphenicol acetyltransferase (CAT) reporter gene. Using deletion and site-specific mutations of 5'-flanking DNA, we identified a complex T3 response unit that contains one major T3 response element (T3RE) and several minor ones. The major element contains two degenerate copies of the hexamer, RGGWMA, separated by 4 bp and was a strong repressor in the absence of ligand. Endogenous levels of T3 receptor are sufficient to allow the T3 response elements in the upstream region of the malic enzyme gene to confer responsiveness to T3, suggesting that they are physiologically relevant.

Localization of the intrinsically bent DNA region upstream of the E.coli rrnB P1 promoter.

DNA sequences upstream of the rrnB P1 core promoter (-10, -35 region) increase transcription more than 300-fold in vivo and in vitro. This stimulation results from a cis-acting DNA sequence, the UP element, which interacts directly with the alpha subunit of RNA polymerase, increasing transcription about 30-fold, and from a positively acting transcription factor, FIS, which increases expression another 10-fold. A DNA region exhibiting a high degree of intrinsic curvature has been observed upstream of the rrnB P1 core promoter and has thus been often cited as an example of the effect of bending on transcription. However, the precise position of the curvature has not been determined. We address here whether this bend is in fact related to activation of rRNA transcription. Electrophoretic analyses were used to localize the major bend in the rrnB P1 upstream region to position approximately -100 with respect to the transcription initiation site. Since most of the effect of upstream sequences on transcription results from DNA between the -35 hexamer and position -88, i.e. downstream of the bend center, these studies indicate that the curvature leading to the unusual electrophoretic behavior of the upstream region does not play a major role in activation of rRNA transcription. Minor deviations from normal electrophoretic behavior were associated with the region just upstream of the -35 hexamer and could conceivably influence interactions between the UP element and the alpha subunit of RNA polymerase.