An introduction to low-level analysis methods of DNA microarray data

This article gives an overview over the methods used in the low--level analysis of gene expression data generated using DNA microarrays. This type of experiment allows to determine relative levels of nucleic acid abundance in a set of tissues or cell populations for thousands of transcripts or loci simultaneously. Careful statistical design and analysis are essential to improve the efficiency and reliability of microarray experiments throughout the data acquisition and analysis process. This includes the design of probes, the experimental design, the image analysis of microarray scanned images, the normalization of fluorescence intensities, the assessment of the quality of microarray data and incorporation of quality information in subsequent analyses, the combination of information across arrays and across sets of experiments, the discovery and recognition of patterns in expression at the single gene and multiple gene levels, and the assessment of significance of these findings, considering the fact that there is a lot of noise and thus random features in the data. For all of these components, access to a flexible and efficient statistical computing environment is an essential aspect.

Fluorinated olefinic peptide nucleic acid: synthesis and pairing properties with complementary DNA

The fluorinated olefinic peptide nucleic acid (F-OPA) system was designed as a peptide nucleic acid (PNA) analogue in which the base carrying amide moiety was replaced by an isostructural and isoelectrostatic fluorinated C-C double bond, locking the nucleobases in one of the two possible rotameric forms. By comparison of the base-pairing properties of this analogue with its nonfluorinated analogue OPA and PNA, we aimed at a closer understanding of the role of this amide function in complementary DNA recognition. Here we present the synthesis of the F-OPA monomer building blocks containing the nucleobases A, T, and G according to the MMTr/Acyl protecting group scheme. Key steps are a selective desymmetrization of the double bond in the monomer precursor via lactonization as well as a highly regioselective Mitsunobu reaction for the introduction of the bases. PNA decamers containing single F-OPA mutations and fully modified F-OPA decamers and pentadecamers containing the bases A and T were synthesized by solid-phase peptide chemistry, and their hybridization properties with complementary parallel and antiparallel DNA were assessed by UV melting curves and CD spectroscopic methods. The stability of the duplexes formed by the decamers containing single (Z)-F-OPA modifications with parallel and antiparallel DNA was found to be strongly dependent on their position in the sequence with T(m) values ranging from +2.4 to -8.1 degrees C/modification as compared to PNA. Fully modified F-OPA decamers and pentadecamers were found to form parallel duplexes with complementary DNA with reduced stability compared to PNA or OPA. An asymmetric F-OPA pentadecamer was found to form a stable self-complex (T(m) approximately 65 degrees C) of unknown structure. The generally reduced affinity to DNA may therefore be due to an increased propensity for self-aggregation

Hybridization of single-stranded DNA targets to immobilized complementary DNA probes: comparison of hairpin versus linear capture probes

A microtiter-based assay system is described in which DNA hairpin probes with dangling ends and single-stranded, linear DNA probes were immobilized and compared based on their ability to capture single-strand target DNA. Hairpin probes consisted of a 16 bp duplex stem, linked by a T2-biotin·dT-T2 loop. The third base was a biotinylated uracil (UB) necessary for coupling to avidin coated microtiter wells. The capture region of the hairpin was a 3′ dangling end composed of either 16 or 32 bases. Fundamental parameters of the system, such as probe density and avidin adsorption capacity of the plates were characterized. The target DNA consisted of 65 bases whose 3′ end was complementary to the dangling end of the hairpin or to the linear probe sequence. The assay system was employed to measure the time dependence and thermodynamic stability of target hybridization with hairpin and linear probes. Target molecules were labeled with either a 5′-FITC, or radiolabeled with [γ-33P]ATP and captured by either linear or hairpin probes affixed to the solid support. Over the range of target concentrations from 10 to 640 pmol hybridization rates increased with increasing target concentration, but varied for the different probes examined. Hairpin probes displayed higher rates of hybridization and larger equilibrium amounts of captured targets than linear probes. At 25 and 45°C, rates of hybridization were better than twice as great for the hairpin compared with the linear capture probes. Hairpin–target complexes were also more thermodynamically stable. Binding free energies were evaluated from the observed equilibrium constants for complex formation. Results showed the order of stability of the probes to be: hairpins with 32 base dangling ends > hairpin probes with l6 base dangling ends > 16 base linear probes > 32 base linear probes. The physical characteristics of hairpins could offer substantial advantages as nucleic acid capture moieties in solid support based hybridization systems.

Complementary DNA Cloning and Characterization of Ferredoxin Localized in Bundle-Sheath Cells of Maize Leaves1

In maize (Zea mays L.) two leaf-specific ferredoxin (Fd) isoproteins, Fd I and Fd II, are distributed differentially in mesophyll and bundle-sheath cells. A novel cDNA encoding the precursor of Fd II (pFD2) was isolated by heterologous hybridization using a cDNA for Fd I (pFD1) as a probe. The assignment of the cDNAs to the Fds was verified by capillary liquid-chromatography/electrospray ionization-mass spectrometry. RNA-blot analysis demonstrated that transcripts for Fd I and Fd II accumulated specifically in mesophyll and bundle-sheath cells, respectively. The mature regions of pFD1 and pFD2 were expressed in Escherichia coli as functional Fds. Fd I and Fd II had similar redox potentials of −423 and −406 mV, respectively, but the Km value of Fd-NADP+ reductase for Fd II was about 3-fold larger than that for Fd I. Asparagine at position 65 of Fd II is a unique residue compared with Fd I and other Fds from various plants, which have aspartic acid or glutamic acid at the corresponding position as an electrostatic interaction site with Fd-NADP+ reductase. Substitution of asparagine-65 with aspartic acid increased the affinity of Fd II with Fd-NADP+ reductase to a level comparable to that of Fd I. These structural and functional differences of Fd I and Fd II may be related to their cell-specific expression in the leaves of a C4 plant.

Detailed molecular characterization of cord blood-derived endothelial progenitors

Objective. Given their involvement in pathological and physiological angiogenesis, there has been growing interest in understanding and manipulating endothellial progenitor cells (EPC) for therapeutic purposes. However, detailed molecular analysis of EPC before and during endothelial differentiation is lacking and is the subject of the present study. Materials and Methods. We report a detailed microarray gene-expression profile of freshly isolated (day 0) human cord blood (CB)-derived EPC (CD133(+)KDR(+) or CD34(+)KDR(+)), and at different time points during in vitro differentiation (early: day 13; late: day 27). Results. Data obtained reflect an EPC transcriptome enriched in genes related to stem/progenitor cells properties (chromatin remodeling, self-renewal, signaling, cytoskeleton organization and biogenesis, recruitment, and adhesion). Using a complementary DNA microarray enriched in intronic transcribed sequences, we observed, as well, that naturally transcribed intronic noncoding RNAs were specifically expressed at the EPC stage. Conclusion. Taken together, we have defined the global gene-expression profile of CB-derived EPC during the process of endothelial differentiation, which can be used to identify genes involved in different vascular pathologies. (C) 2008 ISEH - Society for Hematology and Stem Cells. Published by Elsevier Inc.

Identification of protein expression signatures in gastric carcinomas using clustering analysis

Background and Aim: The identification of gastric carcinomas (GC) has traditionally been based on histomorphology. Recently, DNA microarrays have successfully been used to identify tumors through clustering of the expression profiles. Random forest clustering is widely used for tissue microarrays and other immunohistochemical data, because it handles highly-skewed tumor marker expressions well, and weighs the contribution of each marker according to its relatedness with other tumor markers. In the present study, we e identified biologically- and clinically-meaningful groups of GC by hierarchical clustering analysis of immunohistochemical protein expression. Methods: We selected 28 proteins (p16, p27, p21, cyclin D1, cyclin A, cyclin B1, pRb, p53, c-met, c-erbB-2, vascular endothelial growth factor, transforming growth factor [TGF]-beta I, TGF-beta II, MutS homolog-2, bcl-2, bax, bak, bcl-x, adenomatous polyposis coli, clathrin, E-cadherin, beta-catenin, mucin (MUC) 1, MUC2, MUC5AC, MUC6, matrix metalloproteinase [ MMP]-2, and MMP-9) to be investigated by immunohistochemistry in 482 GC. The analyses of the data were done using a random forest-clustering method. Results: Proteins related to cell cycle, growth factor, cell motility, cell adhesion, apoptosis, and matrix remodeling were highly expressed in GC. We identified protein expressions associated with poor survival in diffuse-type GC. Conclusions: Based on the expression analysis of 28 proteins, we identified two groups of GC that could not be explained by any clinicopathological variables, and a subgroup of long-surviving diffuse-type GC patients with a distinct molecular profile. These results provide not only a new molecular basis for understanding the biological properties of GC, but also better prediction of survival than the classic pathological grouping.

Bone morphogenetic protein-7 is a MYC target with prosurvival functions in childhood medulloblastoma

Medulloblastoma (MB) is the most common malignant brain tumor in children. It is known that overexpression and/or amplification of the MYC oncogene is associated with poor clinical outcome, but the molecular mechanisms and the MYC downstream effectors in MB remain still elusive. Besides contributing to elucidate how progression of MB takes place, most importantly, the identification of novel MYC-target genes will suggest novel candidates for targeted therapy in MB. A group of 209 MYC-responsive genes was obtained from a complementary DNA microarray analysis of a MB-derived cell line, following MYC overexpression and silencing. Among the MYC-responsive genes, we identified the members of the bone morphogenetic protein (BMP) signaling pathway, which have a crucial role during the development of the cerebellum. In particular, the gene BMP7 was identified as a direct target of MYC. A positive correlation between MYC and BMP7 expression was documented by analyzing two distinct sets of primary MB samples. Functional studies in vitro using a small-molecule inhibitor of the BMP/SMAD signaling pathway reproduced the effect of the small interfering RNA-mediated silencing of BMP7. Both approaches led to a block of proliferation in a panel of MB cells and to inhibition of SMAD phosphorylation. Altogether, our findings indicate that high MYC levels drive BMP7 overexpression, promoting cell survival in MB cells. This observation suggests the potential relevance of targeting the BMP/SMAD pathway as a novel therapeutic approach for the treatment of childhood MB.

Expression profiling in spondyloarthropathy synovial biopsies highlights changes in expression of inflammatory genes in conjunction with tissue remodelling genes

Background: In the spondyloarthropathies, the underlying molecular and cellular pathways driving disease are poorly understood. By undertaking a study in knee synovial biopsies from spondyloarthropathy (SpA) and ankylosing spondylitis (AS) patients we aimed to elucidate dysregulated genes and pathways. Methods RNA was extracted from six SpA, two AS, three osteoarthritis (OA) and four normal control knee synovial biopsies. Whole genome expression profiling was undertaken using the Illumina DASL system, which assays 24000 cDNA probes. Differentially expressed candidate genes were then validated using quantitative PCR and immunohistochemistry. Results: Four hundred and sixteen differentially expressed genes were identified that clearly delineated between AS/SpA and control groups. Pathway analysis showed altered gene-expression in oxidoreductase activity, B-cell associated, matrix catabolic, and metabolic pathways. Altered «myogene» profiling was also identified. The inflammatory mediator, MMP3, was strongly upregulated (5-fold) in AS/SpA samples and the Wnt pathway inhibitors DKK3 (2.7-fold) and Kremen1 (1.5-fold) were downregulated. Conclusions: Altered expression profiling in SpA and AS samples demonstrates that disease pathogenesis is associated with both systemic inflammation as well as local tissue alterations that may underlie tissue damaging modelling and remodelling outcomes. This supports the hypothesis that initial systemic inflammation in spondyloarthropathies transfers to and persists in the local joint environment, and might subsequently mediate changes in genes directly involved in the destructive tissue remodelling.

Whole blood transcriptional profiling in ankylosing spondylitis identifies novel candidate genes that might contribute to the inflammatory and tissue-destructive disease aspects

Introduction: A number of genetic-association studies have identified genes contributing to ankylosing spondylitis (AS) susceptibility but such approaches provide little information as to the gene activity changes occurring during the disease process. Transcriptional profiling generates a 'snapshot' of the sampled cells' activity and thus can provide insights into the molecular processes driving the disease process. We undertook a whole-genome microarray approach to identify candidate genes associated with AS and validated these gene-expression changes in a larger sample cohort. Methods: A total of 18 active AS patients, classified according to the New York criteria, and 18 gender- and age-matched controls were profiled using Illumina HT-12 whole-genome expression BeadChips which carry cDNAs for 48,000 genes and transcripts. Class comparison analysis identified a number of differentially expressed candidate genes. These candidate genes were then validated in a larger cohort using qPCR-based TaqMan low density arrays (TLDAs). Results: A total of 239 probes corresponding to 221 genes were identified as being significantly different between patients and controls with a P-value <0.0005 (80% confidence level of false discovery rate). Forty-seven genes were then selected for validation studies, using the TLDAs. Thirteen of these genes were validated in the second patient cohort with 12 downregulated 1.3- to 2-fold and only 1 upregulated (1.6-fold). Among a number of identified genes with well-documented inflammatory roles we also validated genes that might be of great interest to the understanding of AS progression such as SPOCK2 (osteonectin) and EP300, which modulate cartilage and bone metabolism. Conclusions: We have validated a gene expression signature for AS from whole blood and identified strong candidate genes that may play roles in both the inflammatory and joint destruction aspects of the disease.

High-resolution SNP microarray investigation of copy number variations on chromosome 18 in a control cohort

Copy number variations (CNVs) as described in the healthy population are purported to contribute significantly to genetic heterogeneity. Recent studies have described CNVs using lymphoblastoid cell lines or by application of specifically developed algorithms to interrogate previously described data. However, the full extent of CNVs remains unclear. Using high-density SNP array, we have undertaken a comprehensive investigation of chromosome 18 for CNV discovery and characterisation of distribution and association with chromosome architecture. We identified 399 CNVs, of which loss represents 98%, 58% are less than 2.5 kb in size and 71% are intergenic. Intronic deletions account for the majority of copy number changes with gene involvement. Furthermore, one-third of CNVs do not have putative breakpoints within repetitive sequences. We conclude that replicative processes, mediated either by repetitive elements or microhomology, account for the majority of CNVs in the healthy population. Genomic instability involving the formation of a non-B structure is demonstrated in one region.