Comprehensive literature review and statistical considerations for microarray meta-analysis

With the rapid advances of various high-throughput technologies, generation of '-omics' data is commonplace in almost every biomedical field. Effective data management and analytical approaches are essential to fully decipher the biological knowledge contained in the tremendous amount of experimental data. Meta-analysis, a set of statistical tools for combining multiple studies of a related hypothesis, has become popular in genomic research. Here, we perform a systematic search from PubMed and manual collection to obtain 620 genomic meta-analysis papers, of which 333 microarray meta-analysis papers are summarized as the basis of this paper and the other 249 GWAS meta-analysis papers are discussed in the next companion paper. The review in the present paper focuses on various biological purposes of microarray meta-analysis, databases and software and related statistical procedures. Statistical considerations of such an analysis are further scrutinized and illustrated by a case study. Finally, several open questions are listed and discussed.     

Novel Adaptors of Amyloid Precursor Protein Intracellular Domain and Their Functional Implications

Amyloid precursor protein intracellular domain(AICD) is one of the potential candidates in deciphering the complexity of Alzheimer’s disease.It plays important roles in determining cell fate and neurodegeneration through its interactions with several adaptors.The presence or absence of phosphorylation at specific sites determines the choice of partners.In this study,we identified 20 novel AICDinteracting proteins by in vitro pull down experiments followed by 2D gel electrophoresis and MALDI-MS analysis.The identified proteins can be grouped into different functional classes including molecular chaperones,structural proteins,signaling and transport molecules,adaptors,motor proteins and apoptosis determinants.Interactions of nine proteins were further validated either by colocalization using confocal imaging or by co-immunoprecipitation followed by immunoblotting.The cellular functions of most of the proteins can be correlated with AD.Hence,illustration of their interactions with AICD may shed some light on the disease pathophysiology.     

Evidence for an Elevated Aspartate pKa in the Active Site of Human Aromatase

Aromatase (CYP19A1), the enzyme that converts androgens to estrogens, is of significant mechanistic and therapeutic interest. Crystal structures and computational studies of this enzyme shed light on the critical role of Asp³⁰⁹ in substrate binding and catalysis. These studies predicted an elevated pK_a for Asp³⁰⁹ and proposed that protonation of this residue was required for function. In this study, UV-visible absorption, circular dichroism, resonance Raman spectroscopy, and enzyme kinetics were used to study the impact of pH on aromatase structure and androstenedione binding. Spectroscopic studies demonstrate that androstenedione binding is pH-dependent, whereas, in contrast, the D309N mutant retains its ability to bind to androstenedione across the entire pH range studied. Neither pH nor mutation perturbed the secondary structure or heme environment. The origin of the observed pH dependence was further narrowed to the protonation equilibria of Asp³⁰⁹ with a parallel set of spectroscopic studies using exemestane and anastrozole. Because exemestane interacts with Asp³⁰⁹ based on its co-crystal structure with the enzyme, its binding is pH-dependent. Aromatase binding to anastrozole is pH-independent, consistent with the hypothesis that this ligand exploits a distinct set of interactions in the active site. In summary, we assign the apparent pK_a of 8.2 observed for androstenedione binding to the side chain of Asp³⁰⁹. To our knowledge, this work represents the first experimental assignment of a pK_a value to a residue in a cytochrome P450. This value is in agreement with theoretical calculations (7.7–8.1) despite the reliance of the computational methods on the conformational snapshots provided by crystal structures.     

The cell death regulator GRIM-19 is involved in HIV-1 induced T-cell apoptosis

One of the hallmarks of Human Immunodeficiency Virus-1 (HIV-1) infection is progressive depletion of the infected and bystander CD4+ T-cells by apoptosis. Different mitochondrial proteins have been implicated in this apoptotic process; however, the role of different subunits of mitochondrial oxidative phosphorylation (OXPHOS) complexes in apoptosis is not clearly understood. Some of the OXPHOS complex subunits seem to perform other functions in addition to their primary role in energy generating process. GRIM-19 (gene associated with retinoid-interferon-induced-mortality-19), a subunit of mitochondrial complex-I was previously implicated in Interferon-β and retionoic acid induced apoptosis in many tumor cells. In this study we report, using differential gene expression analysis, that GRIM-19 is up-regulated in HIV-1 infected apoptotic T-cells. A temporal up regulation of this subunit was observed in different HIV-1 infected T-cell lines and human PBMC and the extent of increase correlated to increasing apoptosis and virus production. Moreover, silencing GRIM-19 in HIV-1 infected cells reduced apoptosis, indicating its involvement in HIV-1 induced T-cell death.     

Protein folding by domain V of Escherichia coli 23S rRNA: specificity of RNA-protein interactions

The peptidyl transferase center, present in domain V of 23S rRNA of eubacteria and large rRNA of plants and animals, can act as a general protein folding modulator. Here we show that a few specific nucleotides in Escherichia coli domain V RNA bind to unfolded proteins and, as shown previously, bring the trapped proteins to a folding-competent state before releasing them. These nucleotides are the same for the proteins studied so far: bovine carbonic anhydrase, lactate dehydrogenase, malate dehydrogenase, and chicken egg white lysozyme. The amino acids that interact with these nucleotides are also found to be specific in the two cases tested: bovine carbonic anhydrase and lysozyme. They are either neutral or positively charged and are present in random coils on the surface of the crystal structure of both the proteins. In fact, two of these amino acid-nucleotide pairs are identical in the two cases. How these features might help the process of protein folding is discussed.     

"Omics" data and levels of evidence for biomarker discovery

With the development of new technologies for assaying biological activity on a global basis in experimental samples, various new "-omics" signatures have been developed to predict disease progression. Such signatures hold the potential to alter the nature of clinical management of human disease. In this article, we describe some necessary statistical considerations needed to take these signatures from the discovery phase to a clinically useful assay. Much of the work discussed is in the area of cancer.     

Combined analysis of murine and human microarrays and ChIP analysis reveals genes associated with the ability of MYC to maintain tumorigenesis

The MYC oncogene has been implicated in the regulation of up to thousands of genes involved in many cellular programs including proliferation, growth, differentiation, self-renewal, and apoptosis. MYC is thought to induce cancer through an exaggerated effect on these physiologic programs. Which of these genes are responsible for the ability of MYC to initiate and/or maintain tumorigenesis is not clear. Previously, we have shown that upon brief MYC inactivation, some tumors undergo sustained regression. Here we demonstrate that upon MYC inactivation there are global permanent changes in gene expression detected by microarray analysis. By applying StepMiner analysis, we identified genes whose expression most strongly correlated with the ability of MYC to induce a neoplastic state. Notably, genes were identified that exhibited permanent changes in mRNA expression upon MYC inactivation. Importantly, permanent changes in gene expression could be shown by chromatin immunoprecipitation (ChIP) to be associated with permanent changes in the ability of MYC to bind to the promoter regions. Our list of candidate genes associated with tumor maintenance was further refined by comparing our analysis with other published results to generate a gene signature associated with MYC-induced tumorigenesis in mice. To validate the role of gene signatures associated with MYC in human tumorigenesis, we examined the expression of human homologs in 273 published human lymphoma microarray datasets in Affymetrix U133A format. One large functional group of these genes included the ribosomal structural proteins. In addition, we identified a group of genes involved in a diverse array of cellular functions including: BZW2, H2AFY, SFRS3, NAP1L1, NOLA2, UBE2D2, CCNG1, LIFR, FABP3, and EDG1. Hence, through our analysis of gene expression in murine tumor models and human lymphomas, we have identified a novel gene signature correlated with the ability of MYC to maintain tumorigenesis.     

Probing the role of active site histidine residues in the catalytic activity of lacrimal gland peroxidase

The role of active site histidine residues in SCN^– oxidation by lacrimal gland peroxidase (LGP) has been probed after modification with diethylpyrocarbonate (DEPC). The enzyme is irreversibly inactivated following pseudo-first order kinetics with a second order rate constant of 0.26 M^–1 sec^–1 at 25°C. The pH dependent rate of inactivation shows an inflection point at 6.6 indicating histidine derivatization. The UV difference spectrum of the modified versus native enzyme shows a peak at 242 nm indicating formation of N-carbethoxyhistidine. Carbethoxyhistidine formation and associated inactivation are reversed by hydroxylamine indicating histidine modification. The stoichiometry of histidine modification and the extent of inactivation show that out of five histidine residues modified, modification of two residues inactivates the enzyme. Substrate protection with SCN^– during modification indicates that although one histidine is protected, it does not prevent inactivation. The spectroscopically detectable compound II formation is lost due to modification and is not evident after SCN^– protection. The data indicate that out of two histidines, one regulates compound I formation while the other one controls SCN^– binding. SCN^– protected enzyme is inactive due to loss of compound I formation. SCN^– binding studies by optical difference spectroscopy indicate that while the native enzyme binds SCN^– with the K_d of 15 mM, the modified enzyme shows very weak binding with the K_d of 660 mM. From the pH dependent binding of SCN^–, a plot of log K_d vs. pH shows a sigmoidal curve from which the involvement of an enzyme ionizable group of pKa 6.6 is ascertained and attributed to the histidine residue controlling SCN^– binding. LGP has thus two distinctly different essential histidine residues – one regulates compound I formation while the other one controls SCN^– binding.