A copula-based method for synthetic microarray data generation

In this work, we propose a copula-based method to generate synthetic gene expression data that account for marginal and joint probability distributions features captured from real data. Our method allows us to implant significant genes in the synthetic dataset in a controlled manner, giving the possibility of testing new detection algorithms under more realistic environments.

Glyoxysomal malate-dehydrogenase and malate synthase from Soybean cotyledons (Glycine max. L.): Enzyme association, antibody-production and cDNA cloning

In order to investigate a possible association between soybean malate synthase (MS; L-malate glyoxylate-lyase, CoA-acetylating, EC 4.1.3.2) and glyoxysomal malate dehydrogenase (gMDH; (S)-malate: NAD(+) oxidoreductase, EC 1.1.1.37), two consecutive enzymes in the glyoxylate cycle, their elution profiles were analyzed on Superdex 200 HR fast protein liquid chromatography columns equilibrated in low- and high-ionic-strength buffers. Starting with soluble proteins extracted from the cotyledons of 5-d-old soybean seedlings and a 45% ammonium sulfate precipitation, MS and gMDH coeluted on Superdex 200 HR (low-ionic-strength buffer) as a complex with an approximate relative molecular mass (M(r)) of 670000. Dissociation was achieved in the presence of 50 mM KCl and 5 mM MgCl2, with the elution of MS as an octamer of M, 510 000 and of gMDH as a dimer of M, 73 000. Polyclonal antibodies raised to the native copurified enzymes recognized both denatured MS and gMDH on immunoblots, and their native forms after gel filtration. When these antibodies were used to screen a lambda ZAP II expression library containing cDNA from 3-d-old soybean cotyledons, they identified seven clones encoding gMDH, whereas ten clones encoding MS were identified using an antibody to SDS-PAGE-purified MS. Of these cDNA clones a 1.8 kb clone for MS and a 1.3-kb clone for gMDH were fully sequenced. While 88% identity was found between mature soybean gMDH and watermelon gMDH, the N-terminal transit peptides showed only 37% identity. Despite this low identity, the soybean gMDH transit peptide conserves the consensus R(X(6))HL motif also found in plant and mammalian thiolases.

Dual activation of pathways regulated by steroid receptors and peptide growth factors in primary prostate cancer revealed by Factor Analysis of microarray data

Background: We use an approach based on Factor Analysis to analyze datasets generated for transcriptional profiling. The method groups samples into biologically relevant categories, and enables the identification of genes and pathways most significantly associated to each phenotypic group, while allowing for the participation of a given gene in more than one cluster. Genes assigned to each cluster are used for the detection of pathways predominantly activated in that cluster by finding statistically significant associated GO terms. We tested the approach with a published dataset of microarray experiments in yeast. Upon validation with the yeast dataset, we applied the technique to a prostate cancer dataset. Results: Two major pathways are shown to be activated in organ-confined, non-metastatic prostate cancer: those regulated by the androgen receptor and by receptor tyrosine kinases. A number of gene markers (HER3, IQGAP2 and POR1) highlighted by the software and related to the later pathway have been validated experimentally a posteriori on independent samples. Conclusion: Using a new microarray analysis tool followed by a posteriori experimental validation of the results, we have confirmed several putative markers of malignancy associated with peptide growth factor signalling in prostate cancer and revealed others, most notably ERRB3 (HER3). Our study suggest that, in primary prostate cancer, HER3, together or not with HER4, rather than in receptor complexes involving HER2, could play an important role in the biology of these tumors. These results provide new evidence for the role of receptor tyrosine kinases in the establishment and progression of prostate cancer.

A cell spot microarray method for high-throughput biology

High-throughput screening of cellular effects of RNA interference (RNAi) libraries is now being increasingly applied to explore the role of genes in specific cell biological processes and disease states. However, the technology is still limited to specialty laboratories, due to the requirements for robotic infrastructure, access to expensive reagent libraries, expertise in high-throughput screening assay development, standardization, data analysis and applications. In the future, alternative screening platforms will be required to expand functional large-scale experiments to include more RNAi constructs, allow combinatorial loss-of-function analyses (e.g. genegene or gene-drug interaction), gain-of-function screens, multi-parametric phenotypic readouts or comparative analysis of many different cell types. Such comprehensive perturbation of gene networks in cells will require a major increase in the flexibility of the screening platforms, throughput and reduction of costs. As an alternative for the conventional multi-well based high-throughput screening -platforms, here the development of a novel cell spot microarray method for production of high density siRNA reverse transfection arrays is described. The cell spot microarray platform is distinguished from the majority of other transfection cell microarray techniques by the spatially confined array layout that allow highly parallel screening of large-scale RNAi reagent libraries with assays otherwise difficult or not applicable to high-throughput screening. This study depicts the development of the cell spot microarray method along with biological application examples of high-content immunofluorescence and phenotype based cancer cell biological analyses focusing on the regulation of prostate cancer cell growth, maintenance of genomic integrity in breast cancer cells, and functional analysis of integrin protein-protein interactions in situ.

Assessment of the degree of contamination of rat germ cell preparations using specific cDNA probes

Recent reports showing a decrease in sperm count in men have brought new concerns about male infertility. Animal models have been widely used to provide some relevant information about the human male gamete, and extrapolations are made to men and to the clinical context. The present study assesses one of the methods used for separation of germ cells of the adult rat testis, namely centrifugal elutriation followed by density gradients (Percoll®). This method was chosen since it presents the best results for cell purity in separating germ cells from the rat testis. A comparison between continuous and discontinuous Percoll® gradients was performed in order to identify the best type of gradient to separate the cells. Maximal cell purity was obtained for spermatocytes (81 ± 8.2%, mean ± SEM) and spermatids (84 ± 2.6%) using centrifugal elutriation followed by continuous Percoll® gradients. A significant difference in purity was observed between elongating spermatids harvested from continuous Percoll® gradients and from discontinuous gradients. Molecular analysis was used to assess cell contamination by employing specific probes, namely transition protein 2 (TP2), mitochondrial cytochrome C oxidase II (COX II), and sulfated glycoprotein 1 (SGP1). Molecular analysis of the samples demonstrated that morphological criteria are efficient in characterizing the main composition of the cell suspension, but are not reliable for identifying minimal contamination from other cells. Reliable cell purity data should be established using molecular analysis

A tobacco cDNA reveals two different transcription patterns in vegetative and reproductive organs

In order to identify genes expressed in the pistil that may have a role in the reproduction process, we have established an expressed sequence tags project to randomly sequence clones from a Nicotiana tabacum stigma/style cDNA library. A cDNA clone (MTL-8) showing high sequence similarity to genes encoding glycine-rich RNA-binding proteins was chosen for further characterization. Based on the extensive identity of MTL-8 to the RGP-1a sequence of N. sylvestris, a primer was defined to extend the 5' sequence of MTL-8 by RT-PCR from stigma/style RNAs. The amplification product was sequenced and it was confirmed that MTL-8 corresponds to an mRNA encoding a glycine-rich RNA-binding protein. Two transcripts of different sizes and expression patterns were identified when the MTL-8 cDNA insert was used as a probe in RNA blots. The largest is 1,100 nucleotides (nt) long and markedly predominant in ovaries. The smaller transcript, with 600 nt, is ubiquitous to the vegetative and reproductive organs analyzed (roots, stems, leaves, sepals, petals, stamens, stigmas/styles and ovaries). Plants submitted to stress (wounding, virus infection and ethylene treatment) presented an increased level of the 600-nt transcript in leaves, especially after tobacco necrosis virus infection. In contrast, the level of the 1,100-nt transcript seems to be unaffected by the stress conditions tested. Results of Southern blot experiments have suggested that MTL-8 is present in one or two copies in the tobacco genome. Our results suggest that the shorter transcript is related to stress while the larger one is a flower predominant and nonstress-inducible messenger.

Pipeline for macro- and microarray analyses

The pipeline for macro- and microarray analyses (PMmA) is a set of scripts with a web interface developed to analyze DNA array data generated by array image quantification software. PMmA is designed for use with single- or double-color array data and to work as a pipeline in five classes (data format, normalization, data analysis, clustering, and array maps). It can also be used as a plugin in the BioArray Software Environment, an open-source database for array analysis, or used in a local version of the web service. All scripts in PMmA were developed in the PERL programming language and statistical analysis functions were implemented in the R statistical language. Consequently, our package is a platform-independent software. Our algorithms can correctly select almost 90% of the differentially expressed genes, showing a superior performance compared to other methods of analysis. The pipeline software has been applied to 1536 expressed sequence tags macroarray public data of sugarcane exposed to cold for 3 to 48 h. PMmA identified thirty cold-responsive genes previously unidentified in this public dataset. Fourteen genes were up-regulated, two had a variable expression and the other fourteen were down-regulated in the treatments. These new findings certainly were a consequence of using a superior statistical analysis approach, since the original study did not take into account the dependence of data variability on the average signal intensity of each gene. The web interface, supplementary information, and the package source code are available, free, to non-commercial users at http://ipe.cbmeg.unicamp.br/pub/PMmA.

Cloning, sequence analysis, and expression of cDNA coding for the major house dust mite allergen, Der f 1, in Escherichia coli

Our objective was to clone, express and characterize adult Dermatophagoides farinae group 1 (Der f 1) allergens to further produce recombinant allergens for future clinical applications in order to eliminate side reactions from crude extracts of mites. Based on GenBank data, we designed primers and amplified the cDNA fragment coding for Der f 1 by nested-PCR. After purification and recovery, the cDNA fragment was cloned into the pMD19-T vector. The fragment was then sequenced, subcloned into the plasmid pET28a(+), expressed in Escherichia coli BL21 and identified by Western blotting. The cDNA coding for Der f 1 was cloned, sequenced and expressed successfully. Sequence analysis showed the presence of an open reading frame containing 966 bp that encodes a protein of 321 amino acids. Interestingly, homology analysis showed that the Der p 1 shared more than 87% identity in amino acid sequence with Eur m 1 but only 80% with Der f 1. Furthermore, phylogenetic analyses suggested that D. pteronyssinus was evolutionarily closer to Euroglyphus maynei than to D. farinae, even though D. pteronyssinus and D. farinae belong to the same Dermatophagoides genus. A total of three cysteine peptidase active sites were found in the predicted amino acid sequence, including 127-138 (QGGCGSCWAFSG), 267-277 (NYHAVNIVGYG) and 284-303 (YWIVRNSWDTTWGDSGYGYF). Moreover, secondary structure analysis revealed that Der f 1 contained an a helix (33.96%), an extended strand (17.13%), a ß turn (5.61%), and a random coil (43.30%). A simple three-dimensional model of this protein was constructed using a Swiss-model server. The cDNA coding for Der f 1 was cloned, sequenced and expressed successfully. Alignment and phylogenetic analysis suggests that D. pteronyssinus is evolutionarily more similar to E. maynei than to D. farinae.

cDNA and deduced primary structure of basic phospholipase A2 with neurotoxic activity from the venom secretion of the Crotalus durissus collilineatus rattlesnake

To illustrate the construction of precursor complementary DNAs, we isolated mRNAs from whole venom samples. After reverse transcription polymerase chain reaction (RT-PCR), we amplified the cDNA coding for a neurotoxic protein, phospholipase A2 D49 (PLA2 D49), from the venom of Crotalus durissus collilineatus (Cdc PLA2). The cDNA encoding Cdc PLA2 from whole venom was sequenced. The deduced amino acid sequence of this cDNA has high overall sequence identity with the group II PLA2 protein family. Cdc PLA2 has 14 cysteine residues capable of forming seven disulfide bonds that characterize this group of PLA2 enzymes. Cdc PLA2 was isolated using conventional Sephadex G75 column chromatography and reverse-phase high performance liquid chromatography (RP-HPLC). The molecular mass was estimated using matrix-assisted laser desorption ionization-time-of-flight (MALDI-TOF) mass spectrometry. We tested the neuromuscular blocking activities on chick biventer cervicis neuromuscular tissue. Phylogenetic analysis of Cdc PLA2 showed the existence of two lines of N6-PLA2, denominated F24 and S24. Apparently, the sequences of the New World’s N6-F24-PLA2 are similar to those of the agkistrodotoxin from the Asian genus Gloydius. The sequences of N6-S24-PLA2 are similar to the sequence of trimucrotoxin from the genus Protobothrops, found in the Old World.

Ropinirole alters gene expression profiles in SH-SY5Y cells: a whole genome microarray study

Ropinirole (ROP) is a dopamine agonist that has been used as therapy for Parkinson's disease. In the present study, we aimed to detect whether gene expression was modulated by ROP in SH-SY5Y cells. SH-SY5Y cell lines were treated with 10 µM ROP for 2 h, after which total RNA was extracted for whole genome analysis. Gene expression profiling revealed that 113 genes were differentially expressed after ROP treatment compared with control cells. Further pathway analysis revealed modulation of the phosphatidylinositol 3-kinase (PI3K) signaling pathway, with prominent upregulation of PIK3C2B. Moreover, batches of regulated genes, including PIK3C2B, were found to be located on chromosome 1. These findings were validated by quantitative RT-PCR and Western blot analysis. Our study, therefore, revealed that ROP altered gene expression in SH-SY5Y cells, and future investigation of PIK3C2B and other loci on chromosome 1 may provide long-term implications for identifying novel target genes of Parkinson's disease.

Data analysis tools and methods for DNA microarray and high-throughput sequencing data

The recent rapid development of biotechnological approaches has enabled the production of large whole genome level biological data sets. In order to handle thesedata sets, reliable and efficient automated tools and methods for data processingand result interpretation are required. Bioinformatics, as the field of studying andprocessing biological data, tries to answer this need by combining methods and approaches across computer science, statistics, mathematics and engineering to studyand process biological data. The need is also increasing for tools that can be used by the biological researchers themselves who may not have a strong statistical or computational background, which requires creating tools and pipelines with intuitive user interfaces, robust analysis workflows and strong emphasis on result reportingand visualization. Within this thesis, several data analysis tools and methods have been developed for analyzing high-throughput biological data sets. These approaches, coveringseveral aspects of high-throughput data analysis, are specifically aimed for gene expression and genotyping data although in principle they are suitable for analyzing other data types as well. Coherent handling of the data across the various data analysis steps is highly important in order to ensure robust and reliable results. Thus,robust data analysis workflows are also described, putting the developed tools andmethods into a wider context. The choice of the correct analysis method may also depend on the properties of the specific data setandthereforeguidelinesforchoosing an optimal method are given. The data analysis tools, methods and workflows developed within this thesis have been applied to several research studies, of which two representative examplesare included in the thesis. The first study focuses on spermatogenesis in murinetestis and the second one examines cell lineage specification in mouse embryonicstem cells.

cDNA-AFLP na identificação de genes relacionados a qualidade fisiológica

Alguns trabalhos têm evidenciado a existência de genótipos de soja contrastantes para qualidade fisiológica de semente. Tais diferenças existem em virtude da presença de sementes com total ou parcial impermeabilidade do tegumento à água, o que as tornam menos susceptíveis aos danos mecânicos, as adversidades climáticas e a deterioração por umidade. A característica de tegumento semi-permeável pode ser incorporada às cultivares de alta produção por meio dos programas de melhoramento. No entanto, há necessidade de caracterizar os genes ou as regiões genômicas envolvidas com esta característica. Nesse contexto, ferramentas da biologia molecular, como a técnica de cDNA-AFLP, podem auxiliar a identificação de genes relacionados a qualidade fisiológica de sementes. O objetivo desse estudo foi verificar a eficácia da técnica de cDNA-AFLP, na obtenção de fragmentos de genes diferencialmente expressos entre tegumentos de sementes de soja com permeabilidade contrastante. Sementes provenientes dos genótipos CD-202 (tegumento amarelo, permeável e susceptível a deterioração) e TP (tegumento preto, semi-permeável e resistente a deterioração) foram cultivadas em casa-de-vegetação, sob condições homogêneas. Realizou-se a coleta das sementes em diferentes estágios de desenvolvimento (25, 30, 35, 40 e 55 dias após a antese). Procedeu-se à extração do RNA total utilizando-se três métodos, sendo o reagente Pure Link Plant RNA o mais eficiente. Para obtenção do cDNA dupla fita utilizou-se o kit Double Stranded cDNA Synthesis. A partir do cDNA obtido, aplicou-se a técnica de AFLP testando-se um total de 64 combinações de primers. Foram obtidos 47 fragmentos de cDNA diferencialmente expressos entre os tegumentos de sementes dos dois genótipos, os quais poderão ter suas funções reveladas pelo sequenciamento e análise in sílico. De acordo com os resultados obtidos, a técnica de cDNA-AFLP demonstra ser uma alternativa promissora para estudos que visem à identificação de genes relacionados a qualidade de sementes.

Isolation and analysis of a corprinus cinereus DNA fragment showing homology to a fimbrial cDNA from ustilago violacea

Surface proteinaceous fibrils, termed fimbriae, were first identified on gram negative bacteria in the 1940s. Fungal fimbriae, discovered some 25 years later, are found on members of all fungal classes. In the present study, polyclonal antiserum raised against the fimbrial proteins of U. vio/acea were used in order to identify antigenically related proteins from Coprinus cinereus and Schizophy//um commune. Two polypeptides with molecular masses of 37 and 39 kDa from C. cinereus were observed and confirm earlier results. A single previously unidentified 50 kDa polypeptide in S. commune crossreacted with the antiserum. The 50 kDa protein was found to consist of 3 isoforms with isoelectric points ranging from 5.6 to 5.8. A fimbrial cDNA derived from U. vio/acea was used to identify DNA restriction fragments from C. cinereus and S. commune showing homology to the fimbrial transcript of U. vio/acea. Heterologous hybridization with this cDNA was used in order to screen a C. cinereus genomic DNA library. A single clone, A2-3A, with a 14 kbp insert showed strong homology to the pfim3-1 cDNA. The region of homology, a 700 bp Xba I fragment, was subcloned into pUG19. This plasmid was refered to as pXX8. DNA sequence determinations of pXX8 and adjacent fragments from A2-3A suggested that the cloned DNA was a portion of the rONA repeat encoding the small subunit rRNA. DNA sequence analysis of pfim3-1 yielded an incomplete open reading frame. The predicted amino acid sequence codes for a 206 amino acid, 22 kDa polypeptide which contains a domain similar to a transmembrane domain from rat leukocyte antigen, GDS3. As well, an untranslated 576 nucleotide domain showed 81 % homology to pXX8 and 830/0 homology to the 188 rRNA sequence of Ustilago maydis. This sequence was found adjacent to a region of adenine-thymine base pairs presumed to represent the polyadenylation sequence of the fimbrial transcript. The size and extent of homology is sufficient to account for the hybridization of pfim3-1 to rDNA. It is suggested that this domain represents a completely novel regulatory domain within eukaryotes that may enable the observed rapid regeneration of fimbriae in U. violacea.