The MicroArray Quality Control (MAQC)-II study of common practices for the development and validation of microarray-based predictive models.

Gene expression data from microarrays are being applied to predict preclinical and clinical endpoints, but the reliability of these predictions has not been established. In the MAQC-II project, 36 independent teams analyzed six microarray data sets to generate predictive models for classifying a sample with respect to one of 13 endpoints indicative of lung or liver toxicity in rodents, or of breast cancer, multiple myeloma or neuroblastoma in humans. In total, >30,000 models were built using many combinations of analytical methods. The teams generated predictive models without knowing the biological meaning of some of the endpoints and, to mimic clinical reality, tested the models on data that had not been used for training. We found that model performance depended largely on the endpoint and team proficiency and that different approaches generated models of similar performance. The conclusions and recommendations from MAQC-II should be useful for regulatory agencies, study committees and independent investigators that evaluate methods for global gene expression analysis.

A positive FGFR3/FOXN1 feedback loop underlies benign skin keratosis versus squamous cell carcinoma formation in humans.

Seborrheic keratoses (SKs) are common, benign epithelial tumors of the skin that do not, or very rarely, progress into malignancy, for reasons that are not understood. We investigated this by gene expression profiling of human SKs and cutaneous squamous cell carcinomas (SCCs) and found that several genes previously connected with keratinocyte tumor development were similarly modulated in SKs and SCCs, whereas the expression of others differed by only a few fold. In contrast, the tyrosine kinase receptor FGF receptor-3 (FGFR3) and the transcription factor forkhead box N1 (FOXN1) were highly expressed in SKs, and close to undetectable in SCCs. We also showed that increased FGFR3 activity was sufficient to induce FOXN1 expression, counteract the inhibitory effect of EGFR signaling on FOXN1 expression and differentiation, and induce differentiation in a FOXN1-dependent manner. Knockdown of FOXN1 expression in primary human keratinocytes cooperated with oncogenic RAS in the induction of SCC-like tumors, whereas increased FOXN1 expression triggered the SCC cells to shift to a benign SK-like tumor phenotype, which included increased FGFR3 expression. Thus,we have uncovered a positive regulatory loop between FGFR3 and FOXN1 that underlies a benign versus malignant skin tumor phenotype.

Integrative analysis of RUNX1 downstream pathways and target genes.

BACKGROUND: The RUNX1 transcription factor gene is frequently mutated in sporadic myeloid and lymphoid leukemia through translocation, point mutation or amplification. It is also responsible for a familial platelet disorder with predisposition to acute myeloid leukemia (FPD-AML). The disruption of the largely unknown biological pathways controlled by RUNX1 is likely to be responsible for the development of leukemia. We have used multiple microarray platforms and bioinformatic techniques to help identify these biological pathways to aid in the understanding of why RUNX1 mutations lead to leukemia. RESULTS: Here we report genes regulated either directly or indirectly by RUNX1 based on the study of gene expression profiles generated from 3 different human and mouse platforms. The platforms used were global gene expression profiling of: 1) cell lines with RUNX1 mutations from FPD-AML patients, 2) over-expression of RUNX1 and CBFbeta, and 3) Runx1 knockout mouse embryos using either cDNA or Affymetrix microarrays. We observe that our datasets (lists of differentially expressed genes) significantly correlate with published microarray data from sporadic AML patients with mutations in either RUNX1 or its cofactor, CBFbeta. A number of biological processes were identified among the differentially expressed genes and functional assays suggest that heterozygous RUNX1 point mutations in patients with FPD-AML impair cell proliferation, microtubule dynamics and possibly genetic stability. In addition, analysis of the regulatory regions of the differentially expressed genes has for the first time systematically identified numerous potential novel RUNX1 target genes. CONCLUSION: This work is the first large-scale study attempting to identify the genetic networks regulated by RUNX1, a master regulator in the development of the hematopoietic system and leukemia. The biological pathways and target genes controlled by RUNX1 will have considerable importance in disease progression in both familial and sporadic leukemia as well as therapeutic implications

Affymetrix Whole-Transcript Human Gene 1.0 ST array is highly concordant with standard 3' expression arrays.

The recently released Affymetrix Human Gene 1.0 ST array has two major differences compared with standard 3' based arrays: (i) it interrogates the entire mRNA transcript, and (ii) it uses DNA targets. To assess the impact of these differences on array performance, we performed a series of comparative hybridizations between the Human Gene 1.0 ST and the Affymetrix HG-U133 Plus 2.0 and the Illumina HumanRef-8 BeadChip arrays. Additionally, both RNA and DNA targets were hybridized on HG-U133 Plus 2.0 arrays. The results show that the overall reproducibility of the Gene 1.0 ST array is best. When looking only at the high intensity probes, the reproducibility of the Gene 1.0 ST array and the Illumina BeadChip array is equally good. Concordance of array results was assessed using different inter-platform mappings. Agreements are best between the two labeling protocols using HG-U133 Plus 2.0 array. The Gene 1.0 ST array is most concordant with the HG-U133 array hybridized with cDNA targets. This may reflect the impact of the target type. Overall, the high degree of correspondence provides strong evidence for the reliability of the Gene 1.0 ST array.

Fourteen new cases contribute to the characterization of the 7q11.23 microduplication syndrome.

Interstitial deletions of 7q11.23 cause Williams-Beuren syndrome, one of the best characterized microdeletion syndromes. The clinical phenotype associated with the reciprocal duplication however is not well defined, though speech delay is often mentioned. We present 14 new 7q11.23 patients with the reciprocal duplication of the Williams-Beuren syndrome critical region, nine familial and five de novo. These were identified by either array-based MLPA or by array-CGH/oligonucleotide analysis in a series of patients with idiopathic mental retardation with an estimated population frequency of 1:13,000-1:20,000. Variable speech delay is a constant finding in our patient group, confirming previous reports. Cognitive abilities range from normal to moderate mental retardation. The association with autism is present in five patients and in one father who also carries the duplication. There is an increased incidence of hypotonia and congenital anomalies: heart defects (PDA), diaphragmatic hernia, cryptorchidism and non-specific brain abnormalities on MRI. Specific dysmorphic features were noted in our patients, including a short philtrum, thin lips and straight eyebrows. Our patient collection demonstrates that the 7q11.23 microduplication not only causes language delay, but is also associated with congenital anomalies and a recognizable face.

Pregnancy-induced chromatin remodeling in the breast of postmenopausal women.

Early pregnancy and multiparity are known to reduce the risk of women to develop breast cancer at menopause. Based on the knowledge that the differentiation of the breast induced by the hormones of pregnancy plays a major role in this protection, this work was performed with the purpose of identifying what differentiation-associated molecular changes persist in the breast until menopause. Core needle biopsies (CNB) obtained from the breast of 42 nulliparous (NP) and 71 parous (P) postmenopausal women were analyzed in morphology, immunocytochemistry and gene expression. Whereas in the NP breast, nuclei of epithelial cells were large and euchromatic, in the P breast they were small and hyperchromatic, showing strong methylation of histone 3 at lysine 9 and 27. Transcriptomic analysis performed using Affymetrix HG_U133 oligonucleotide arrays revealed that in CNB of the P breast, there were 267 upregulated probesets that comprised genes controlling chromatin organization, transcription regulation, splicing machinery, mRNA processing and noncoding elements including XIST. We concluded that the differentiation process induced by pregnancy is centered in chromatin remodeling and in the mRNA processing reactome, both of which emerge as important regulatory pathways. These are indicative of a safeguard step that maintains the fidelity of the transcription process, becoming the ultimate mechanism mediating the protection of the breast conferred by full-term pregnancy.

An exquisite cross-control mechanism among endothelial cell fate regulators directs the plasticity and heterogeneity of lymphatic endothelial cells.

Arteriovenous-lymphatic endothelial cell fates are specified by the master regulators, namely, Notch, COUP-TFII, and Prox1. Whereas Notch is expressed in the arteries and COUP-TFII in the veins, the lymphatics express all 3 cell fate regulators. Previous studies show that lymphatic endothelial cell (LEC) fate is highly plastic and reversible, raising a new concept that all 3 endothelial cell fates may co-reside in LECs and a subtle alteration can result in a reprogramming of LEC fate. We provide a molecular basis verifying this concept by identifying a cross-control mechanism among these cell fate regulators. We found that Notch signal down-regulates Prox1 and COUP-TFII through Hey1 and Hey2 and that activated Notch receptor suppresses the lymphatic phenotypes and induces the arterial cell fate. On the contrary, Prox1 and COUP-TFII attenuate vascular endothelial growth factor signaling, known to induce Notch, by repressing vascular endothelial growth factor receptor-2 and neuropilin-1. We show that previously reported podoplanin-based LEC heterogeneity is associated with differential expression of Notch1 in human cutaneous lymphatics. We propose that the expression of the 3 cell fate regulators is controlled by an exquisite feedback mechanism working in LECs and that LEC fate is a consequence of the Prox1-directed lymphatic equilibrium among the cell fate regulators.

Cross-species and cross-platform gene expression studies with the Bioconductor-compliant R package 'annotationTools'.

Background: The variety of DNA microarray formats and datasets presently available offers an unprecedented opportunity to perform insightful comparisons of heterogeneous data. Cross-species studies, in particular, have the power of identifying conserved, functionally important molecular processes. Validation of discoveries can now often be performed in readily available public data which frequently requires cross-platform studies.Cross-platform and cross-species analyses require matching probes on different microarray formats. This can be achieved using the information in microarray annotations and additional molecular biology databases, such as orthology databases. Although annotations and other biological information are stored using modern database models ( e. g. relational), they are very often distributed and shared as tables in text files, i.e. flat file databases. This common flat database format thus provides a simple and robust solution to flexibly integrate various sources of information and a basis for the combined analysis of heterogeneous gene expression profiles.Results: We provide annotationTools, a Bioconductor-compliant R package to annotate microarray experiments and integrate heterogeneous gene expression profiles using annotation and other molecular biology information available as flat file databases. First, annotationTools contains a specialized set of functions for mining this widely used database format in a systematic manner. It thus offers a straightforward solution for annotating microarray experiments. Second, building on these basic functions and relying on the combination of information from several databases, it provides tools to easily perform cross-species analyses of gene expression data.Here, we present two example applications of annotationTools that are of direct relevance for the analysis of heterogeneous gene expression profiles, namely a cross-platform mapping of probes and a cross-species mapping of orthologous probes using different orthology databases. We also show how to perform an explorative comparison of disease-related transcriptional changes in human patients and in a genetic mouse model.Conclusion: The R package annotationTools provides a simple solution to handle microarray annotation and orthology tables, as well as other flat molecular biology databases. Thereby, it allows easy integration and analysis of heterogeneous microarray experiments across different technological platforms or species.

Involvement of connexin 43 in meiotic maturation of bovine oocytes.

In ovarian follicles, cumulus cells provide the oocyte with small molecules that permit growth and control maturation. These nutrients reach the germinal cell through gap junction channels, which are present between the cumulus cells and the oocyte, and between the cumulus cells. In this study the involvement of intercellular communication mediated by gap junction channels on oocyte maturation of in vitro cultured bovine cumulus-oocyte complexes (COCs) was investigated. The stages of oocyte maturation were determined by Hoechst 33342 staining, which showed that 90% of COCs placed in the maturation medium for 24 h progress to the metaphase II stage. Bovine COC gap junction communication was disrupted initially using n-alkanols, which inhibit any passage through gap junctions. In the presence of 1-heptanol (3 mmol l(-1)) or octanol (3.0 mmol l(-1) and 0.3 mmol l(-1)), only 29% of the COCs reached metaphase II. Removal of the uncoupling agent was associated with restoration of oocyte maturation, indicating that treatment with n-alkanols was neither cytotoxic nor irreversible. Concentrations of connexin 43 (Cx43), the major gap junction protein expressed in the COCs, were decreased specifically using a recombinant adenovirus expressing the antisense Cx43 cDNA (Ad-asCx43). The efficacy of adenoviral infection was > 95% in cumulus cells evaluated after infection with recombinant adenoviruses expressing the green fluorescence protein. RT-PCR performed on total RNA isolated from Ad-asCx43-infected COCs showed that the rat Cx43 cDNA was transcribed. Western blot analysis revealed a three-fold decrease in Cx43 expression in COCs expressing the antisense RNA for Cx43. Injection of cumulus cells with Lucifer yellow demonstrated further that the resulting lower amount of Cx43 in infected COCs is associated with a two-fold decrease in the extent of coupling between cumulus cells. In addition, oocyte maturation was decreased by 50% in the infected COC cultures. These results indicate that Cx43-mediated communication between cumulus cells plays a crucial role in maturation of bovine oocytes.

Pancreatic-specific expression of the glucose transporter type 2 gene: identification of cis-elements and islet-specific trans-acting factors.

A defect in glucose sensing of the pancreatic beta-cells has been observed in several animal models of type II diabetes and has been correlated with a reduced gene expression of the glucose transporter type 2 (Glut2). In a transgenic mouse model, expression of Glut2 antisense RNA in pancreatic beta-cells has recently been shown to be associated with an impaired glucose-induced insulin secretion and the development of diabetes. To identify factors that may be involved in the specific decrease of Glut2 in the beta-cells of the diabetic animal, an attempt was made to localize the cis-elements and trans-acting factors involved in the control of Glut2 expression in the endocrine pancreas. It was demonstrated by transient transfection studies that only 338 base pairs (bp) of the murine Glut2 proximal promoter are needed for reporter gene expression in pancreatic islet-derived cell lines, whereas no activity was detected in nonpancreatic cells. Three cis-elements, GTI, GTII, and GTIII, have been identified by DNAse I footprinting and gel retardation experiments within these 338 bp. GTI and GTIII bind distinct but ubiquitously expressed trans-acting factors. On the other hand, nuclear proteins specifically expressed in pancreatic cell lines interact with GTII, and their relative abundance correlates with endogenous Glut2 expression. These GTII-binding factors correspond to nuclear proteins of 180 and 90 kilodaltons as defined by Southwestern analysis. The 180-kilodalton factor is present in pancreatic beta-cell lines but not in an alpha-cell line. Mutation of the GTI or GTIII cis-elements decreases transcriptional activity directed by the 338-bp promoter, whereas mutation of GTII increases gene transcription. Thus negative and positive regulatory sequences are identified within the proximal 338 bp of the GLUT2 promoter and may participate in the islet-specific expression of the gene by binding beta-cell specific trans-acting factors.

The scaffold protein IB1/JIP-1 controls the activation of JNK in rat stressed urothelium.

The c-Jun N-terminal kinase (JNK) is critical for cell survival, differentiation, apoptosis and tumorigenesis. This signalling pathway requires the presence of the scaffold protein Islet-Brain1/c-Jun N-terminal kinase interacting protein-1 (IB1/JIP-1). Immunolabeling and in situ hybridisation of bladder sections showed that IB1/JIP-1 is expressed in urothelial cells. The functional role of IB1/JIP-1 in the urothelium was therefore studied in vivo in a model of complete rat bladder outlet obstruction. This parietal stress, which is due to urine retention, reduced the content of IB1/JIP-1 in urothelial cells and consequently induced a drastic increase in JNK activity and AP-1 binding activity. Using a viral gene transfer approach, the stress-induced activation of JNK was prevented by overexpressing IB1/JIP-1. Conversely, the JNK activity was increased in urothelial cells where the IB1/JIP-1 content was experimentally reduced using an antisense RNA strategy. Furthermore, JNK activation was found to be increased in non-stressed urothelial cells of heterozygous mice carrying a selective disruption of the IB1/JIP-1 gene. These data established that mechanical stress in urothelial cells in vivo induces a robust JNK activation as a consequence of regulated expression of the scaffold protein IB1/JIP-1. This result highlights a critical role for that scaffold protein in the homeostasis of the urothelium and unravels a new potential target to regulate the JNK pathway in this tissue.

Analysis of epidermis- and mesophyll-specific transcript accumulation in powdery mildew-inoculated wheat leaves.

Powdery mildew is an important disease of wheat caused by the obligate biotrophic fungus Blumeria graminis f. sp. tritici. This pathogen invades exclusively epidermal cells after penetrating directly through the cell wall. Because powdery mildew colonizes exclusively epidermal cells, it is of importance not only to identify genes which are activated, but also to monitor tissue specificity of gene activation. Acquired resistance of wheat to powdery mildew can be induced by a previous inoculation with the non-host pathogen B. graminis f. sp. hordei, the causal agent of barley powdery mildew. The establishment of the resistant state is accompanied by the activation of genes. Here we report the tissue-specific cDNA-AFLP analysis and cloning of transcripts accumulating 6 and 24 h after the resistance-inducing inoculation with B. graminis f. sp. hordei. A total of 25,000 fragments estimated to represent about 17,000 transcripts were displayed. Out of these, 141 transcripts, were found to accumulate after Bgh inoculation using microarray hybridization analysis. Forty-four accumulated predominantly in the epidermis whereas 76 transcripts accumulated mostly in mesophyll tissue.

The mitochondrial targeting function of randomly generated peptide sequences correlates with predicted helical amphiphilicity.

A pool of oligonucleotides encoding a start methionine and nine random amino acids was inserted at the 5'-end of the gene for the yeast cytochrome oxidase subunit IV lacking its own mitochondrial targeting sequence. Approximately one-quarter of the randomly generated sequences targeted subunit IV to its correct intramitochondrial location in vivo. Sequence analysis of 89 randomly generated sequences showed that their efficiencies as mitochondrial targeting signals correlated with the potential to fold into an amphiphilic alpha-helix. Functional targeting sequences were enriched in arginine and isoleucine residues but contained few aspartate, glutamate, and proline residues. Nonfunctional sequences predicted to have significant helical amphiphilicity often had at least one acidic or multiple helix-breaking residues that would be expected to interfere with targeting functioning. These results support the hypothesis that the signal for targeting a protein into the mitochondrial matrix is usually a positively charged amphiphilic helix.

Identification of molecular apocrine breast tumours by microarray analysis.

Previous microarray studies on breast cancer identified multiple tumour classes, of which the most prominent, named luminal and basal, differ in expression of the oestrogen receptor alpha gene (ER). We report here the identification of a group of breast tumours with increased androgen signalling and a 'molecular apocrine' gene expression profile. Tumour samples from 49 patients with large operable or locally advanced breast cancers were tested on Affymetrix U133A gene expression microarrays. Principal components analysis and hierarchical clustering split the tumours into three groups: basal, luminal and a group we call molecular apocrine. All of the molecular apocrine tumours have strong apocrine features on histological examination (P=0.0002). The molecular apocrine group is androgen receptor (AR) positive and contains all of the ER-negative tumours outside the basal group. Kolmogorov-Smirnov testing indicates that oestrogen signalling is most active in the luminal group, and androgen signalling is most active in the molecular apocrine group. ERBB2 amplification is commoner in the molecular apocrine than the other groups. Genes that best split the three groups were identified by Wilcoxon test. Correlation of the average expression profile of these genes in our data with the expression profile of individual tumours in four published breast cancer studies suggest that molecular apocrine tumours represent 8-14% of tumours in these studies. Our data show that it is possible with microarray data to divide mammary tumour cells into three groups based on steroid receptor activity: luminal (ER+ AR+), basal (ER- AR-) and molecular apocrine (ER- AR+).