Identification of sample-specific sequences in mammalian cDNA and genomic DNA by the novel ligation-mediated subtraction (Limes)

The representational difference analysis (RDA) and other subtraction techniques are used to enrich sample-specific sequences by elimination of ubiquitous sequences existing in both the sample of interest (tester) and the subtraction partner (driver). While applying the RDA to genomic DNA of cutaneous lymphoma cells in order to identify tumor relevant alterations, we predominantly isolated repetitive sequences and artificial repeat-mediated fusion products of otherwise independent PCR fragments (PCR hybrids). Since these products severely interfered with the isolation of tester-specific fragments, we developed a considerably more robust and efficient approach, termed ligation-mediated subtraction (Limes). In first applications of Limes, genomic sequences and/or transcripts of genes involved in the regulation of transcription, such as transforming growth factor β stimulated clone 22 related gene (TSC-22R), cell death and cytokine production (caspase-1) or antigen presentation (HLA class II sequences), were found to be completely absent in a cutaneous lymphoma line. On the assumption that mutations in tumor-relevant genes can affect their transcription pattern, a protocol was developed and successfully applied that allows the identification of such sequences. Due to these results, Limes may substitute/supplement other subtraction/comparison techniques such as RDA or DNA microarray techniques in a variety of different research fields.

A genomic approach to gene fusion technology

Gene expression profiling provides powerful analyses of transcriptional responses to cellular perturbation. In contrast to DNA array-based methods, reporter gene technology has been underused for this application. Here we describe a genomewide, genome-registered collection of Escherichia coli bioluminescent reporter gene fusions. DNA sequences from plasmid-borne, random fusions of E. coli chromosomal DNA to a Photorhabdus luminescens luxCDABE reporter allowed precise mapping of each fusion. The utility of this collection covering about 30% of the transcriptional units was tested by analyzing individual fusions representative of heat shock, SOS, OxyR, SoxRS, and cya/crp stress-responsive regulons. Each fusion strain responded as anticipated to environmental conditions known to activate the corresponding regulatory circuit. Thus, the collection mirrors E. coli's transcriptional wiring diagram. This genomewide collection of gene fusions provides an independent test of results from other gene expression analyses. Accordingly, a DNA microarray-based analysis of mitomycin C-treated E. coli indicated elevated expression of expected and unanticipated genes. Selected luxCDABE fusions corresponding to these up-regulated genes were used to confirm or contradict the DNA microarray results. The power of partnering gene fusion and DNA microarray technology to discover promoters and define operons was demonstrated when data from both suggested that a cluster of 20 genes encoding production of type I extracellular polysaccharide in E. coli form a single operon.

GenEST, a powerful bidirectional link between cDNA sequence data and gene expression profiles generated by cDNA-AFLP

The release of vast quantities of DNA sequence data by large-scale genome and expressed sequence tag (EST) projects underlines the necessity for the development of efficient and inexpensive ways to link sequence databases with temporal and spatial expression profiles. Here we demonstrate the power of linking cDNA sequence data (including EST sequences) with transcript profiles revealed by cDNA-AFLP, a highly reproducible differential display method based on restriction enzyme digests and selective amplification under high stringency conditions. We have developed a computer program (GenEST) that predicts the sizes of virtual transcript-derived fragments (TDFs) of in silico-digested cDNA sequences retrieved from databases. The vast majority of the resulting virtual TDFs could be traced back among the thousands of TDFs displayed on cDNA-AFLP gels. Sequencing of the corresponding bands excised from cDNA-AFLP gels revealed no inconsistencies. As a consequence, cDNA sequence databases can be screened very efficiently to identify genes with relevant expression profiles. The other way round, it is possible to switch from cDNA-AFLP gels to sequences in the databases. Using the restriction enzyme recognition sites, the primer extensions and the estimated TDF size as identifiers, the DNA sequence(s) corresponding to a TDF with an interesting expression pattern can be identified. In this paper we show examples in both directions by analyzing the plant parasitic nematode Globodera rostochiensis. Various novel pathogenicity factors were identified by combining ESTs from the infective stage juveniles with expression profiles of ∼4000 genes in five developmental stages produced by cDNA-AFLP.

Amygdala-enriched genes identified by microarray technology are restricted to specific amygdaloid subnuclei

Microarray technology represents a potentially powerful method for identifying cell type- and regionally restricted genes expressed in the brain. Here we have combined a microarray analysis of differential gene expression among five selected brain regions, including the amygdala, cerebellum, hippocampus, olfactory bulb, and periaqueductal gray, with in situ hybridization. On average, 0.3% of the 34,000 genes interrogated were highly enriched in each of the five regions, relative to the others. In situ hybridization performed on a subset of amygdala-enriched genes confirmed in most cases the overall region-specificity predicted by the microarray data and identified additional sites of brain expression not examined on the microarrays. Strikingly, the majority of these genes exhibited boundaries of expression within the amygdala corresponding to cytoarchitectonically defined subnuclei. These results define a unique set of molecular markers for amygdaloid subnuclei and provide tools to genetically dissect their functional roles in different emotional behaviors.

Exploring the DNA-binding specificities of zinc fingers with DNA microarrays

A key step in the regulation of networks that control gene expression is the sequence-specific binding of transcription factors to their DNA recognition sites. A more complete understanding of these DNA–protein interactions will permit a more comprehensive and quantitative mapping of the regulatory pathways within cells, as well as a deeper understanding of the potential functions of individual genes regulated by newly identified DNA-binding sites. Here we describe a DNA microarray-based method to characterize sequence-specific DNA recognition by zinc-finger proteins. A phage display library, prepared by randomizing critical amino acid residues in the second of three fingers of the mouse Zif268 domain, provided a rich source of zinc-finger proteins with variant DNA-binding specificities. Microarrays containing all possible 3-bp binding sites for the variable zinc fingers permitted the quantitation of the binding site preferences of the entire library, pools of zinc fingers corresponding to different rounds of selection from this library, as well as individual Zif268 variants that were isolated from the library by using specific DNA sequences. The results demonstrate the feasibility of using DNA microarrays for genome-wide identification of putative transcription factor-binding sites.

Screening insertion libraries for mutations in many genes simultaneously using DNA microarrays

We describe a method to screen pools of DNA from multiple transposon lines for insertions in many genes simultaneously. We use thermal asymmetric interlaced–PCR, a hemispecific PCR amplification protocol that combines nested, insertion-specific primers with degenerate primers, to amplify DNA flanking the transposons. In reconstruction experiments with previously characterized Arabidopsis lines carrying insertions of the maize Dissociation (Ds) transposon, we show that fluorescently labeled, transposon-flanking fragments overlapping ORFs hybridize to cognate expressed sequence tags (ESTs) on a DNA microarray. We further show that insertions can be detected in DNA pools from as many as 100 plants representing different transposon lines and that all of the tested, transposon-disrupted genes whose flanking fragments can be amplified individually also can be detected when amplified from the pool. The ability of a transposon-flanking fragment to hybridize declines rapidly with decreasing homology to the spotted DNA fragment, so that only ESTs with >90% homology to the transposon-disrupted gene exhibit significant cross-hybridization. Because thermal asymmetric interlaced–PCR fragments tend to be short, use of the present method favors recovery of insertions in and near genes. We apply the technique to screening pools of new Ds lines using cDNA microarrays containing ESTs for ≈1,000 stress-induced and -repressed Arabidopsis genes.

p53-dependent growth arrest of REF52 cells containing newly amplified DNA.

The rat cell line REF52 is not permissive for gene amplification. Simian virus 40 tumor (T) antigen converts these cells to a permissive state, as do dominant negative mutants of p53, suggesting that the effect of T antigen is due mainly to its ability to bind to p53. To manipulate permissivity, we introduced a temperature-sensitive mutant of T antigen (tsA58) into REF52 cells and selected for resistance to N-(phosphonacetyl)-L-aspartate (PALA). Most freshly isolated PALA-resistant colonies, each of approximately 200 cells, selected at a permissive temperature, arrested when shifted to a nonpermissive temperature. Growth arrest was stable, with no evidence of apoptosis, as long as T antigen was absent but was reversed when T antigen was restored. In contrast, PALA-resistant clones grown to approximately 10(7) cells at a permissive temperature did not arrest when shifted to a nonpermissive temperature. All PALA-resistant clones examined had amplified carbamoyl-phosphate synthetase-aspartate transcarbamoylase-dihydroorotase (CAD) genes, present in structures consistent with a mechanism involving bridge-breakage-fusion (BBF) cycles. We propose that p53-mediated growth arrest operates only early during the complex process of gene amplification, when newly formed PALA-resistant cells contain broken DNA, generated in BBF cycles. During propagation under permissive conditions, the broken DNA ends are healed, and, even though the p53-mediated pathway is still intact at a nonpermissive temperature and the cells contain amplified DNA, they are not arrested in the absence of broken DNA. The data support the hypothesis that BBF cycles are an important mechanism of amplification and that the broken DNA generated in each cycle is a key signal that regulates permissivity for gene amplification.

Análise de expressão gênica de Xylella fastidiosa cultivada em diferentes concentrações de glicose pela técnica de microarrays de DNA

Xylella fastidiosa (Xf) é o agente etiológico causador de doenças em uma grande variedade de cultivos de grande importância econômica, causando a c1orose variegada dos citros, uma das doenças mais danosas à indústria de citros no Brasil. Os genomas de algumas cepas deste fitopatógeno foram completamente seqüenciados promovendo assim estudos funcionais do genoma em larga escala. Neste trabalho nós nos propusemos a investigar o perfil de transcrição de Xf através da técnica microarranjos (no título da dissertação microarrays, mas a partir de agora usaremos microaarranjos) de DNA usando todo o genoma do fitopatógeno e cultivando-a sob meio definido variando a concentração de glicose. O objetivo inicial deste trabalho era observar se Xf comportava-se da mesma forma que Xac, que tem a expressão de goma aumentada devido ao aumento da concentração de glicose do meio. Nossas análises revelaram que enquanto os transcritos relacionados à goma não se mostraram afetados com a concentração de glicose, genes que codificam para análogos a Colicina-V e precursores de fimbria foram induzidos em alta concentração de glicose. Baseados nestes resultados, nós propusemos um modelo de mecanismo de produção e defesa contra a Colicina em Xf.

Transcriptional profiles of Haloferax mediterranei based on nitrogen availability

The haloarchaeon Haloferax mediterranei is able to grow in the presence of different inorganic and organic nitrogen sources by means of the assimilatory pathway under aerobic conditions. In order to identify genes of potential importance in nitrogen metabolism and its regulation in the halophilic microorganism, we have analysed its global gene expression in three culture media with different nitrogen sources: (a) cells were grown stationary and exponentially in ammonium, (b) cells were grown exponentially in nitrate, and (c) cells were shifted to nitrogen starvation conditions. The main differences in the transcriptional profiles have been identified between the cultures with ammonium as nitrogen source and the cultures with nitrate or nitrogen starvation, supporting previous results which indicate the absence of ammonium as the factor responsible for the expression of genes involved in nitrate assimilation pathway. The results have also permitted the identification of transcriptional regulators and changes in metabolic pathways related to the catabolism and anabolism of amino acids or nucleotides. The microarray data was validated by real-time quantitative PCR on 4 selected genes involved in nitrogen metabolism. This work represents the first transcriptional profiles study related to nitrogen assimilation metabolism in extreme halophilic microorganisms using microarray technology.

Tissue-specific gene expression in soybean (Glycine max) detected by cDNA microarray analysis

We have constructed cDNA microarrays for soybean (Glycine max L. Merrill), containing approximately 4,100 Unigene ESTs derived from axenic roots, to evaluate their application and utility for functional genomics of organ differentiation in legumes. We assessed microarray technology by conducting studies to evaluate the accuracy of microarray data and have found them to be both reliable and reproducible in repeat hybridisations. Several ESTs showed high levels (>50 fold) of differential expression in either root or shoot tissue of soybean. A small number of physiologically interesting, and differentially expressed sequences found by microarray analysis were verified by both quantitative real-time RT-PCR and Northern blot analysis. There was a linear correlation (r(2) = 0.99, over 5 orders of magnitude) between microarray and quantitative real-time RT-PCR data. Microarray analysis of soybean has enormous potential not only for the discovery of new genes involved in tissue differentiation and function, but also to study the expression of previously characterised genes, gene networks and gene interactions in wild-type, mutant or transgenic; plants.

Genetic study of alcoholism and novel gene expression in the alcoholic brain

Alcohol dependence may result from neuroadaptation involving alteration of gene expression after long-term alcohol exposure. The systematic study of gene expression profiles of the human alcoholic brain was initiated using the method of polymerase chain reaction (PCR)-differential display and was followed by DNA microarray. To date, more than 100 alcohol-responsive genes have been identified from the frontal cortex, motor cortex and nucleus accumbens of the human brain. These genes have a wide range of functions in the brain and indicate diverse actions of alcohol on neuronal function. This review discusses the current information on the genetic basis of alcoholism and the induction and characterization of these alcohol-responsive genes.

Uncertainty and topographic visualisations

This thesis is a study of low-dimensional visualisation methods for data visualisation under certainty of the input data. It focuses on the two main feed-forward neural network algorithms which are NeuroScale and Generative Topographic Mapping (GTM) by trying to make both algorithms able to accommodate the uncertainty. The two models are shown not to work well under high levels of noise within the data and need to be modified. The modification of both models, NeuroScale and GTM, are verified by using synthetic data to show their ability to accommodate the noise. The thesis is interested in the controversy surrounding the non-uniqueness of predictive gene lists (PGL) of predicting prognosis outcome of breast cancer patients as available in DNA microarray experiments. Many of these studies have ignored the uncertainty issue resulting in random correlations of sparse model selection in high dimensional spaces. The visualisation techniques are used to confirm that the patients involved in such medical studies are intrinsically unclassifiable on the basis of provided PGL evidence. This additional category of ‘unclassifiable’ should be accommodated within medical decision support systems if serious errors and unnecessary adjuvant therapy are to be avoided.

Amino acids in oral drug delivery:salts, ion-pairs and transcriptomics

Oral drug delivery is considered the most popular route of delivery because of the ease of administration, availability of a wide range of dosage forms and the large surface area for drug absorption via the intestinal membrane. However, besides the unfavourable biopharmaceutical properties of the therapeutic agents, efflux transporters such as Pglycoprotein (P-gp) and multiple resistance proteins (MRP) decrease the overall drug uptake by extruding the drug from the cells. Although, prodrugs have been investigated to improve drug partitioning by masking the polar groups covalently with pre-moieties promoting increased uptake, they present significant challenges including reduced solubility and increased toxicity. The current work investigates the use of amino acids as ion-pairs for three model drugs: indomethacin (weak acid), trimethoprim (weak base) and ciprofloxacin (zwitter ion) in an attempt to improve both solubility and uptake. Solubility was studied by salt formation while creating new routes for uptake across the membranes via amino acids transporter proteins or dipeptidyl transporters was the rationale to enhance absorption. New salts were prepared for the model drugs and the oppositely charged amino acids by freeze drying and they were characterised using FTIR, 1HNMR, DSC, SEM, pH solubility profile, solubility and dissolution. Permeability profiles were assessed using an in vitro cell based method; Caco-2 cells and the genetic changes occurring across the transporter genes and various pathways involved in the cellular activities were studied using DNA microarrays. Solubility data showed a significant increase in drug solubility upon preparing the new salts with the oppositely charged counter ions (ciprofloxacin glutamate salt exhibiting 2.9x103 fold enhancement when compared to the free drug). Moreover, permeability studies showed a 3 fold increase in trimethoprim and indomethacin permeabilities upon ion-pairing with amino acids and more than 10 fold when the zwitter ionic drug was paired with glutamic acid. Microarray data revealed that trimethoprim was absorbed actively via OCTN1 transporters while MRP7 is the main transporter gene that mediates its efflux. The absorption of trimethoprim from trimethoprim glutamic acid ion-paired formulations was affected by the ratio of glutamic acid in the formulation which was inversely proportional to the degree of expression of OCTN1. Interestingly, ciprofloxacin glutamic acid ion-pairs were found to decrease the up-regulation of ciprofloxacin efflux proteins (P-gp and MRP4) and over-express two solute carrier transporters; (PEPT2 and SLCO1A2) suggesting that a high aqueous binding constant (K11aq) enables the ion-paired formulations to be absorbed as one entity. In conclusion, formation of ion-pairs with amino acids can influence in a positive way solubility, transfer and gene expression effects of drugs.

The application of reduced-size short tandem repeat primer sets in the analysis of human DNA from degraded and compromised samples

The purpose of this research was to demonstrate the applicability of reduced-size STR (Miniplex) primer sets to challenging samples and to provide the forensic community with new information regarding the analysis of degraded and inhibited DNA. The Miniplex primer sets were validated in accordance with guidelines set forth by the Scientific Working Group on DNA Analysis Methods (SWGDAM) in order to demonstrate the scientific validity of the kits. The Miniplex sets were also used in the analysis of DNA extracted from human skeletal remains and telogen hair. In addition, a method for evaluating the mechanism of PCR inhibition was developed using qPCR. The Miniplexes were demonstrated to be a robust and sensitive tool for the analysis of DNA with as low as 100 pg of template DNA. They also proved to be better than commercial kits in the analysis of DNA from human skeletal remains, with 64% of samples tested producing full profiles, compared to 16% for a commercial kit. The Miniplexes also produced amplification of nuclear DNA from human telogen hairs, with partial profiles obtained from as low as 60 pg of template DNA. These data suggest smaller PCR amplicons may provide a useful alternative to mitochondrial DNA for forensic analysis of degraded DNA from human skeletal remains, telogen hairs, and other challenging samples. In the evaluation of inhibition by qPCR, the effect of amplicon length and primer melting temperature was evaluated in order to determine the binding mechanisms of different PCR inhibitors. Several mechanisms were indicated by the inhibitors tested, including binding of the polymerase, binding to the DNA, and effects on the processivity of the polymerase during primer extension. The data obtained from qPCR illustrated a method by which the type of inhibitor could be inferred in forensic samples, and some methods of reducing inhibition for specific inhibitors were demonstrated. An understanding of the mechanism of the inhibitors found in forensic samples will allow analysts to select the proper methods for inhibition removal or the type of analysis that can be performed, and will increase the information that can be obtained from inhibited samples.

New Advancements of Scalable Statistical Methods for Learning Latent Structures in Big Data

Constant technology advances have caused data explosion in recent years. Accord- ingly modern statistical and machine learning methods must be adapted to deal with complex and heterogeneous data types. This phenomenon is particularly true for an- alyzing biological data. For example DNA sequence data can be viewed as categorical variables with each nucleotide taking four different categories. The gene expression data, depending on the quantitative technology, could be continuous numbers or counts. With the advancement of high-throughput technology, the abundance of such data becomes unprecedentedly rich. Therefore efficient statistical approaches are crucial in this big data era.

Previous statistical methods for big data often aim to find low dimensional struc- tures in the observed data. For example in a factor analysis model a latent Gaussian distributed multivariate vector is assumed. With this assumption a factor model produces a low rank estimation of the covariance of the observed variables. Another example is the latent Dirichlet allocation model for documents. The mixture pro- portions of topics, represented by a Dirichlet distributed variable, is assumed. This dissertation proposes several novel extensions to the previous statistical methods that are developed to address challenges in big data. Those novel methods are applied in multiple real world applications including construction of condition specific gene co-expression networks, estimating shared topics among newsgroups, analysis of pro- moter sequences, analysis of political-economics risk data and estimating population structure from genotype data.