A comparison of methods for differential expression analysis of RNA-seq data.

BACKGROUND: Finding genes that are differentially expressed between conditions is an integral part of understanding the molecular basis of phenotypic variation. In the past decades, DNA microarrays have been used extensively to quantify the abundance of mRNA corresponding to different genes, and more recently high-throughput sequencing of cDNA (RNA-seq) has emerged as a powerful competitor. As the cost of sequencing decreases, it is conceivable that the use of RNA-seq for differential expression analysis will increase rapidly. To exploit the possibilities and address the challenges posed by this relatively new type of data, a number of software packages have been developed especially for differential expression analysis of RNA-seq data. RESULTS: We conducted an extensive comparison of eleven methods for differential expression analysis of RNA-seq data. All methods are freely available within the R framework and take as input a matrix of counts, i.e. the number of reads mapping to each genomic feature of interest in each of a number of samples. We evaluate the methods based on both simulated data and real RNA-seq data. CONCLUSIONS: Very small sample sizes, which are still common in RNA-seq experiments, impose problems for all evaluated methods and any results obtained under such conditions should be interpreted with caution. For larger sample sizes, the methods combining a variance-stabilizing transformation with the 'limma' method for differential expression analysis perform well under many different conditions, as does the nonparametric SAMseq method.

Comparative modular analysis of gene expression in vertebrate organs.

ABSTRACT: BACKGROUND: The degree of conservation of gene expression between homologous organs largely remains an open question. Several recent studies reported some evidence in favor of such conservation. Most studies compute organs' similarity across all orthologous genes, whereas the expression level of many genes are not informative about organ specificity. RESULTS: Here, we use a modularization algorithm to overcome this limitation through the identification of inter-species co-modules of organs and genes. We identify such co-modules using mouse and human microarray expression data. They are functionally coherent both in terms of genes and of organs from both organisms. We show that a large proportion of genes belonging to the same co-module are orthologous between mouse and human. Moreover, their zebrafish orthologs also tend to be expressed in the corresponding homologous organs. Notable exceptions to the general pattern of conservation are the testis and the olfactory bulb. Interestingly, some co-modules consist of single organs, while others combine several functionally related organs. For instance, amygdala, cerebral cortex, hypothalamus and spinal cord form a clearly discernible unit of expression, both in mouse and human. CONCLUSIONS: Our study provides a new framework for comparative analysis which will be applicable also to other sets of large-scale phenotypic data collected across different species.

Proof of concept for microarray-based detection of DNA-binding oncogenes in cell extracts.

The function of DNA-binding proteins is controlled not just by their abundance, but mainly at the level of their activity in terms of their interactions with DNA and protein targets. Moreover, the affinity of such transcription factors to their target sequences is often controlled by co-factors and/or modifications that are not easily assessed from biological samples. Here, we describe a scalable method for monitoring protein-DNA interactions on a microarray surface. This approach was designed to determine the DNA-binding activity of proteins in crude cell extracts, complementing conventional expression profiling arrays. Enzymatic labeling of DNA enables direct normalization of the protein binding to the microarray, allowing the estimation of relative binding affinities. Using DNA sequences covering a range of affinities, we show that the new microarray-based method yields binding strength estimates similar to low-throughput gel mobility-shift assays. The microarray is also of high sensitivity, as it allows the detection of a rare DNA-binding protein from breast cancer cells, the human tumor suppressor AP-2. This approach thus mediates precise and robust assessment of the activity of DNA-binding proteins and takes present DNA-binding assays to a high throughput level.

Monitoring of malaria parasite resistance to chloroquine and sulphadoxine-pyrimethamine in the Solomon Islands by DNA microarray technology.

BACKGROUND: Little information is available on resistance to anti-malarial drugs in the Solomon Islands (SI). The analysis of single nucleotide polymorphisms (SNPs) in drug resistance associated parasite genes is a potential alternative to classical time- and resource-consuming in vivo studies to monitor drug resistance. Mutations in pfmdr1 and pfcrt were shown to indicate chloroquine (CQ) resistance, mutations in pfdhfr and pfdhps indicate sulphadoxine-pyrimethamine (SP) resistance, and mutations in pfATPase6 indicate resistance to artemisinin derivatives. METHODS: The relationship between the rate of treatment failure among 25 symptomatic Plasmodium falciparum-infected patients presenting at the clinic and the pattern of resistance-associated SNPs in P. falciparum infecting 76 asymptomatic individuals from the surrounding population was investigated. The study was conducted in the SI in 2004. Patients presenting at a local clinic with microscopically confirmed P. falciparum malaria were recruited and treated with CQ+SP. Rates of treatment failure were estimated during a 28-day follow-up period. In parallel, a DNA microarray technology was used to analyse mutations associated with CQ, SP, and artemisinin derivative resistance among samples from the asymptomatic community. Mutation and haplotype frequencies were determined, as well as the multiplicity of infection. RESULTS: The in vivo study showed an efficacy of 88% for CQ+SP to treat P. falciparum infections. DNA microarray analyses indicated a low diversity in the parasite population with one major haplotype present in 98.7% of the cases. It was composed of fixed mutations at position 86 in pfmdr1, positions 72, 75, 76, 220, 326 and 356 in pfcrt, and positions 59 and 108 in pfdhfr. No mutation was observed in pfdhps or in pfATPase6. The mean multiplicity of infection was 1.39. CONCLUSION: This work provides the first insight into drug resistance markers of P. falciparum in the SI. The obtained results indicated the presence of a very homogenous P. falciparum population circulating in the community. Although CQ+SP could still clear most infections, seven fixed mutations associated with CQ resistance and two fixed mutations related to SP resistance were observed. Whether the absence of mutations in pfATPase6 indicates the efficacy of artemisinin derivatives remains to be proven.

Meta-analysis of heat- and chemically upregulated chaperone genes in plant and human cells.

Molecular chaperones are central to cellular protein homeostasis. In mammals, protein misfolding diseases and aging cause inflammation and progressive tissue loss, in correlation with the accumulation of toxic protein aggregates and the defective expression of chaperone genes. Bacteria and non-diseased, non-aged eukaryotic cells effectively respond to heat shock by inducing the accumulation of heat-shock proteins (HSPs), many of which molecular chaperones involved in protein homeostasis, in reducing stress damages and promoting cellular recovery and thermotolerance. We performed a meta-analysis of published microarray data and compared expression profiles of HSP genes from mammalian and plant cells in response to heat or isothermal treatments with drugs. The differences and overlaps between HSP and chaperone genes were analyzed, and expression patterns were clustered and organized in a network. HSPs and chaperones only partly overlapped. Heat-shock induced a subset of chaperones primarily targeted to the cytoplasm and organelles but not to the endoplasmic reticulum, which organized into a network with a central core of Hsp90s, Hsp70s, and sHSPs. Heat was best mimicked by isothermal treatments with Hsp90 inhibitors, whereas less toxic drugs, some of which non-steroidal anti-inflammatory drugs, weakly expressed different subsets of Hsp chaperones. This type of analysis may uncover new HSP-inducing drugs to improve protein homeostasis in misfolding and aging diseases.

Integrative molecular bioinformatics study of human adrenocortical tumors : microRNA, tissue-specific target prediction, and pathway analysis.

MicroRNAs (miRs) are involved in the pathogenesis of several neoplasms; however, there are no data on their expression patterns and possible roles in adrenocortical tumors. Our objective was to study adrenocortical tumors by an integrative bioinformatics analysis involving miR and transcriptomics profiling, pathway analysis, and a novel, tissue-specific miR target prediction approach. Thirty-six tissue samples including normal adrenocortical tissues, benign adenomas, and adrenocortical carcinomas (ACC) were studied by simultaneous miR and mRNA profiling. A novel data-processing software was used to identify all predicted miR-mRNA interactions retrieved from PicTar, TargetScan, and miRBase. Tissue-specific target prediction was achieved by filtering out mRNAs with undetectable expression and searching for mRNA targets with inverse expression alterations as their regulatory miRs. Target sets and significant microarray data were subjected to Ingenuity Pathway Analysis. Six miRs with significantly different expression were found. miR-184 and miR-503 showed significantly higher, whereas miR-511 and miR-214 showed significantly lower expression in ACCs than in other groups. Expression of miR-210 was significantly lower in cortisol-secreting adenomas than in ACCs. By calculating the difference between dCT(miR-511) and dCT(miR-503) (delta cycle threshold), ACCs could be distinguished from benign adenomas with high sensitivity and specificity. Pathway analysis revealed the possible involvement of G2/M checkpoint damage in ACC pathogenesis. To our knowledge, this is the first report describing miR expression patterns and pathway analysis in sporadic adrenocortical tumors. miR biomarkers may be helpful for the diagnosis of adrenocortical malignancy. This tissue-specific target prediction approach may be used in other tumors too.

Genotyping microarray for CSNB-associated genes.

PURPOSE: Congenital stationary night blindness (CSNB) is a clinically and genetically heterogeneous retinal disease. Although electroretinographic (ERG) measurements can discriminate clinical subgroups, the identification of the underlying genetic defects has been complicated for CSNB because of genetic heterogeneity, the uncertainty about the mode of inheritance, and time-consuming and costly mutation scanning and direct sequencing approaches. METHODS: To overcome these challenges and to generate a time- and cost-efficient mutation screening tool, the authors developed a CSNB genotyping microarray with arrayed primer extension (APEX) technology. To cover as many mutations as possible, a comprehensive literature search was performed, and DNA samples from a cohort of patients with CSNB were first sequenced directly in known CSNB genes. Subsequently, oligonucleotides were designed representing 126 sequence variations in RHO, CABP4, CACNA1F, CACNA2D4, GNAT1, GRM6, NYX, PDE6B, and SAG and spotted on the chip. RESULTS: Direct sequencing of genes known to be associated with CSNB in the study cohort revealed 21 mutations (12 novel and 9 previously reported). The resultant microarray containing oligonucleotides, which allow to detect 126 known and novel mutations, was 100% effective in determining the expected sequence changes in all known samples assessed. In addition, investigation of 34 patients with CSNB who were previously not genotyped revealed sequence variants in 18%, of which 15% are thought to be disease-causing mutations. CONCLUSIONS: This relatively inexpensive first-pass genetic testing device for patients with a diagnosis of CSNB will improve molecular diagnostics and genetic counseling of patients and their families and gives the opportunity to analyze whether, for example, more progressive disorders such as cone or cone-rod dystrophies underlie the same gene defects.

Gene expression in cultured endometrium from women with different outcomes following IVF.

Estradiol and progesterone are crucial for the acquisition of receptivity and the change in transcriptional activity of target genes in the implantation window. The aim of this study was to differentiate the regulation of genes in the endometrium of patients with recurrent implantation failure (IF) versus those who became pregnant after in vitro fertilization (IVF) treatment. Moreover, the effect of embryo-derived factors on endometrial transcriptional activity was studied. Nine women with known IVF outcome (IF, M, miscarriage, OP, ongoing pregnancy) and undergoing hysteroscopy with endometrial biopsy were enrolled. Biopsies were taken during the midluteal phase. After culture in the presence of embryo-conditioned IVF media, total RNA was extracted and submitted to reverse transcription, target cDNA synthesis, biotin labelling, fragmentation and hybridization using the Affymetrix Human Genome U133A 2.0 Chip. Differential expression of selected genes was re-analysed by quantitative PCR, in which the results were calculated as threshold cycle differences between the groups and normalized to Glyceraldehyde phosphate dehydrogenase and beta-actin. Differences were seen for several genes from endometrial tissue between the IF and the pregnancy groups, and when comparing OP with M, 1875 up- and 1807 down-regulated genes were returned. Real-time PCR analysis confirmed up-regulation for somatostatin, PLAP-2, mucin 4 and CD163, and down-regulation of glycodelin, IL-24, CD69, leukaemia inhibitory factor and prolactin receptor between Op and M. When the different embryo-conditioned media were compared, no significant differential regulation could be demonstrated. Although microarray profiling may currently not be sensitive enough for studying the effects of embryo-derived factors on the endometrium, the observed differences in gene expression between M and OP suggest that it will become an interesting tool for the identification of fertility-relevant markers produced by the endometrium.

Odorant binding proteins of the red imported fire ant, Solenopsis invicta: an example of the problems facing the analysis of widely divergent proteins.

We describe the odorant binding proteins (OBPs) of the red imported fire ant, Solenopsis invicta, obtained from analyses of an EST library and separate 454 sequencing runs of two normalized cDNA libraries. We identified a total of 18 putative functional OBPs in this ant. A third of the fire ant OBPs are orthologs to honey bee OBPs. Another third of the OBPs belong to a lineage-specific expansion, which is a common feature of insect OBP evolution. Like other OBPs, the different fire ant OBPs share little sequence similarity (∼ 20%), rendering evolutionary analyses difficult. We discuss the resulting problems with sequence alignment, phylogenetic analysis, and tests of selection. As previously suggested, our results underscore the importance for careful exploration of the sensitivity to the effects of alignment methods for data comprising widely divergent sequences.

Interactome analysis reveals that FAM161A, deficient in recessive retinitis pigmentosa, is a component of the Golgi-centrosomal network.

Defects in FAM161A, a protein of unknown function localized at the cilium of retinal photoreceptor cells, cause retinitis pigmentosa, a form of hereditary blindness. By using different fragments of this protein as baits to screen cDNA libraries of human and bovine retinas, we defined a yeast two-hybrid-based FAM161A interactome, identifying 53 bona fide partners. In addition to statistically significant enrichment in ciliary proteins, as expected, this interactome revealed a substantial bias towards proteins from the Golgi apparatus, the centrosome and the microtubule network. Validation of interaction with key partners by co-immunoprecipitation and proximity ligation assay confirmed that FAM161A is a member of the recently recognized Golgi-centrosomal interactome, a network of proteins interconnecting Golgi maintenance, intracellular transport and centrosome organization. Notable FAM161A interactors included AKAP9, FIP3, GOLGA3, KIFC3, KLC2, PDE4DIP, NIN and TRIP11. Furthermore, analysis of FAM161A localization during the cell cycle revealed that this protein followed the centrosome during all stages of mitosis, likely reflecting a specific compartmentalization related to its role at the ciliary basal body during the G0 phase. Altogether, these findings suggest that FAM161A's activities are probably not limited to ciliary tasks but also extend to more general cellular functions, highlighting possible novel mechanisms for the molecular pathology of retinal disease.

Analysis of the nucleotide sequence of the coat protein and 3'-untranslated region of two Brazilian Potato virus Y isolates

Two Brazilian Potato virus Y (PVY) isolates were biologically characterized as necrotic (PVY-NBR) and common (PVY-OBR) based upon symptoms on test plants. Additional characterization was performed by sequencing a cDNA corresponding to the 3' terminal region of the viral genome. The sequence consisted of 195 nucleotides (nt) coding part of the nuclear inclusion body b (NIb) gene, 804 nt of the coat protein (CP) gene, and 328 nt (PVY-OBR) or 326 nt (PVY-NBR) of the 3'-untranslated region (UTR). Translation of the sequence resulted in one single open reading frame with part of the NIb and a CP of 267 amino acids. The two isolates shared 95.1% similarity in the CP amino acid sequence. The CP and the 3'-UTR sequence of the Brazilian isolates were compared to those of other PVY isolates previously reported and unrooted phylogenetic trees were constructed. The trees revealed a separation of two distinct clusters, one comprising most of the common strains and the other comprising the necrotic strains. PVY-OBR was clustered in the common group and PVY-NBR in the necrotic one.

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

Functional analysis of newly discovered growth control genes: experimental approaches

A large number of DNA sequences corresponding to human and animal transcripts have been filed in data banks, as cDNAs or ESTs (expression sequence tags). However, the actual function of their corresponding gene products is still largely unknown. Several of these genes may play a role in regulation of important biological processes such as cell division, differentiation, malignant transformation and oncogenesis. Elucidation of gene function is based on 2 main approaches, namely, overexpression and expression interference, which respectively mimick or suppress a given phenotype. The currently available tools and experimental approaches to gene functional analysis and the most recent advances in mass cDNA screening by functional analysis are discussed.

Comparative sequence analysis of the P-, M- and L-coding region of the measles virus CAM-70 live attenuated vaccine strain

Measles virus is a highly contagious agent which causes a major health problem in developing countries. The viral genomic RNA is single-stranded, nonsegmented and of negative polarity. Many live attenuated vaccines for measles virus have been developed using either the prototype Edmonston strain or other locally isolated measles strains. Despite the diverse geographic origins of the vaccine viruses and the different attenuation methods used, there was remarkable sequence similarity of H, F and N genes among all vaccine strains. CAM-70 is a Japanese measles attenuated vaccine strain widely used in Brazilian children and produced by Bio-Manguinhos since 1982. Previous studies have characterized this vaccine biologically and genomically. Nevertheless, only the F, H and N genes have been sequenced. In the present study we have sequenced the remaining P, M and L genes (approximately 1.6, 1.4 and 6.5 kb, respectively) to complete the genomic characterization of CAM-70 and to assess the extent of genetic relationship between CAM-70 and other current vaccines. These genes were amplified using long-range or standard RT-PCR techniques, and the cDNA was cloned and automatically sequenced using the dideoxy chain-termination method. The sequence analysis comparing previously sequenced genotype A strains with the CAM-70 Bio-Manguinhos strain showed a low divergence among them. However, the CAM-70 strains (CAM-70 Bio-Manguinhos and a recently sequenced CAM-70 submaster seed strain) were assigned to a specific group by phylogenetic analysis using the neighbor-joining method. Information about our product at the genomic level is important for monitoring vaccination campaigns and for future studies of measles virus attenuation.