Distribution and early development of microarray technology in Europe

A Mechanistic Beta-Binomial Probability Model for mRNA Sequencing Data

A main application for mRNA sequencing (mRNAseq) is determining lists of differentially-expressed genes (DEGs) between two or more conditions. Several software packages exist to produce DEGs from mRNAseq data, but they typically yield different DEGs, sometimes markedly so. The underlying probability model used to describe mRNAseq data is central to deriving DEGs, and not surprisingly most softwares use different models and assumptions to analyze mRNAseq data. Here, we propose a mechanistic justification to model mRNAseq as a binomial process, with data from technical replicates given by a binomial distribution, and data from biological replicates well-described by a beta-binomial distribution. We demonstrate good agreement of this model with two large datasets. We show that an emergent feature of the beta-binomial distribution, given parameter regimes typical for mRNAseq experiments, is the well-known quadratic polynomial scaling of variance with the mean. The so-called dispersion parameter controls this scaling, and our analysis suggests that the dispersion parameter is a continually decreasing function of the mean, as opposed to current approaches that impose an asymptotic value to the dispersion parameter at moderate mean read counts. We show how this leads to current approaches overestimating variance for moderately to highly expressed genes, which inflates false negative rates. Describing mRNAseq data with a beta-binomial distribution thus may be preferred since its parameters are relatable to the mechanistic underpinnings of the technique and may improve the consistency of DEG analysis across softwares, particularly for moderately to highly expressed genes.

DNA binding by estrogen receptor-alpha is essential for the transcriptional response to estrogen in the liver and the uterus

The majority of the biological effects of estrogens in the reproductive tract are mediated by estrogen receptor (ER)alpha, which regulates transcription by several mechanisms. Because the tissue-specific effects of some ERalpha ligands may be caused by tissue-specific transcriptional mechanisms of ERalpha, we aimed to identify the contribution of DNA recognition to these mechanisms in two clinically important target organs, namely uterus and liver. We used a genetic mouse model that dissects DNA binding-dependent vs. independent transcriptional regulation elicited by ERalpha. The EAAE mutant harbors amino acid exchanges at four positions of the DNA-binding domain (DBD) of ERalpha. This construct was knocked in the ERalpha gene locus to produce ERalpha((EAAE/EAAE)) mice devoid of a functional ERalpha DBD. The phenotype of the ERalpha((EAAE/EAAE)) mice resembles the general loss-of-function phenotype of alphaER knockout mutant mice with hypoplastic uteri, hemorrhagic ovaries, and impaired mammary gland development. In agreement with this phenotype, the expression pattern of the ERalpha((EAAE/EAAE)) mutant mice in liver obtained by genome-wide gene expression profiling supports the observation of a near-complete loss of estrogen-dependent gene regulation in comparison with the wild type. Further gene expression analyses to validate the results of the microarray data were performed by quantitative RT-PCR. The analyses indicate that both gene activation and repression by estrogen-bound ERalpha rely on an intact DBD in vivo.

The bifunctional flavokinase/flavin adenine dinucleotide synthetase from Streptomyces davawensis produces inactive flavin cofactors and is not involved in resistance to the antibiotic roseoflavin

Streptomyces davawensis synthesizes the antibiotic roseoflavin, one of the few known natural riboflavin analogs, and is roseoflavin resistant. It is thought that the endogenous flavokinase (EC 2.7.1.26)/flavin adenine dinucleotide (FAD) synthetase (EC 2.7.7.2) activities of roseoflavin-sensitive organisms are responsible for the antibiotic effect of roseoflavin, producing the inactive cofactors roseoflavin-5'-monophosphate (RoFMN) and roseoflavin adenine dinucleotide (RoFAD) from roseoflavin. To confirm this, the FAD-dependent Sus scrofa D-amino acid oxidase (EC 1.4.3.3) was tested with RoFAD as a cofactor and found to be inactive. It was hypothesized that a flavokinase/FAD synthetase (RibC) highly specific for riboflavin may be present in S. davawensis, which would not allow the formation of toxic RoFMN/RoFAD. The gene ribC from S. davawensis was cloned. RibC from S. davawensis was overproduced in Escherichia coli and purified. Analysis of the flavokinase activity of RibC revealed that the S. davawensis enzyme is not riboflavin specific (roseoflavin, kcat/Km = 1.7 10(-2) microM(-1) s(-1); riboflavin, kcat/Km = 7.5 10(-3) microM(-1) s(-1)). Similar results were obtained for RibC from the roseoflavin-sensitive bacterium Bacillus subtilis (roseoflavin, kcat/Km = 1.3 10(-2) microM(-1) s(-1); riboflavin, kcat/Km = 1.3 10(-2) microM(-1) s(-1)). Both RibC enzymes synthesized RoFAD and RoFMN. The functional expression of S. davawensis ribC did not confer roseoflavin resistance to a ribC-defective B. subtilis strain.

Thyrocyte-specific Gq/G11 deficiency impairs thyroid function and prevents goiter development

The function of the adult thyroid is regulated by thyroid-stimulating hormone (TSH), which acts through a G protein-coupled receptor. Overactivation of the TSH receptor results in hyperthyroidism and goiter. The Gs-mediated stimulation of adenylyl cyclase-dependent cAMP formation has been regarded as the principal intracellular signaling mechanism mediating the action of TSH. Here we show that the Gq/G11-mediated signaling pathway plays an unexpected and essential role in the regulation of thyroid function. Mice lacking the alpha subunits of Gq and G11 specifically in thyroid epithelial cells showed severely reduced iodine organification and thyroid hormone secretion in response to TSH, and many developed hypothyroidism within months after birth. In addition, thyrocyte-specific Galphaq/Galpha11-deficient mice lacked the normal proliferative thyroid response to TSH or goitrogenic diet, indicating an essential role of this pathway in the adaptive growth of the thyroid gland. Our data suggest that Gq/G11 and their downstream effectors are promising targets to interfere with increased thyroid function and growth.

Identification and characterization of two Streptomyces davawensis riboflavin biosynthesis gene clusters

In Streptomyces davawensis roseoflavin is synthesized from GTP and ribulose-5-phosphate through riboflavin. As a first step towards the molecular analysis of flavin metabolism in S. davawensis the genes involved in riboflavin biosynthesis were cloned by hybridization of heterologous probes to a genomic library on a high-density colony-array. The genes ribB (riboflavin synthase, alpha-chain; EC 2.5.1.9), ribM (putative membrane protein), ribA (bifunctional GTP cyclohydrolase II/3,4-dihydroxy-2-butanone-4-phosphate synthase; EC 3.5.4.25) and ribH (lumazine synthase; EC 2.5.1.9) are organized in an operon-like cluster. Northern blot analysis of this cluster revealed two transcripts of 1.7 and 3.1 kb, respectively. The gene ribB was overexpressed in Escherichia coli. The specific riboflavin synthase activity in a cell-free extract of a recombinant strain was 0.246 nmol mg(-1 )min(-1). Overexpression of ribM enhanced the transport of riboflavin in the corresponding recombinant E. coli strain. Furthermore, overexpression of ribM increased roseoflavin sensitivity of E. coli. On another subgenomic fragment a putative S. davawensis ribG gene coding for the missing pyrimidine deaminase/reductase (EC 3.5.4.26 and EC 1.1.1.193) of the riboflavin biosynthetic pathway and ribY coding for a second (monofunctional) GTP cyclohydrolase II were identified.

Introduction: array technology - an overview

Microarray technology has its roots in high-throughput parallel synthesis of biomacromolecules, combined with combinatorial science. In principle, the preparation of arrays can be performed either by in situ synthesis of biomacromolecules on solid substrates or by spotting of ex situ synthesized biomacromolecules onto the substrate surface. The application of microarrays includes spatial addressing with target (macro) molecules and screening for interactions between immobilized probe and target. The screening is simplified by the microarray format, which features a known structure of every immobilized library element. The area of nucleic acid arrays is best developed, because such arrays are allowed to follow the biosynthetic pathway from genes to proteins, and because nucleic acid hybridization is a most straightforward screening tool. Applications to genomics, transcriptomics, proteomics, and glycomics are currently in the foreground of interest; in this postgenomic phase they are allowed to gain new insights into the molecular basis of cellular processes and the development of disease.

Definition of estrogen receptor pathway critical for estrogen positive feedback to gonadotropin-releasing hormone neurons and fertility

The mechanisms through which estrogen regulates gonadotropin-releasing hormone (GnRH) neurons to control mammalian ovulation are unknown. We found that estrogen positive feedback to generate the preovulatory gonadotropin surge was normal in estrogen receptor beta knockout (ERbeta) mutant mice, but absent in ERalpha mutant mice. An ERalpha-selective compound was sufficient to generate positive feedback in wild-type mice. As GnRH neurons do not express ERalpha, estrogen positive feedback upon GnRH neurons must be indirect in nature. To establish the cell type responsible, we generated a neuron-specific ERalpha mutant mouse line. These mice failed to exhibit estrogen positive feedback, demonstrating that neurons expressing ERalpha are critical. We then used a GnRH neuron-specific Pseudorabies virus (PRV) tracing approach to show that the ERalpha-expressing neurons innervating GnRH neurons are located within rostral periventricular regions of the hypothalamus. These studies demonstrate that ovulation is driven by estrogen actions upon ERalpha-expressing neuronal afferents to GnRH neurons.

Differential requirements for COPI transport during vertebrate early development

The coatomer vesicular coat complex is essential for normal Golgi and secretory activities in eukaryotic cells. Through positional cloning of genes controlling zebrafish notochord development, we found that the sneezy, happy, and dopey loci encode the alpha, beta, and beta' subunits of the coatomer complex. Export from mutant endoplasmic reticulum is blocked, Golgi structure is disrupted, and mutant embryos eventually degenerate due to widespread apoptosis. The early embryonic phenotype, however, demonstrates that despite its "housekeeping" functions, coatomer activity is specifically and cell autonomously required for normal chordamesoderm differentiation, perinotochordal basement membrane formation, and melanophore pigmentation. Hence, differential requirements for coatomer activity among embryonic tissues lead to tissue-specific developmental defects. Moreover, we note that the mRNA encoding alpha coatomer is strikingly upregulated in notochord progenitors, and we present data suggesting that alpha coatomer transcription is tuned to activity- and cell type-specific secretory loads.

Construction and initial analysis of bacterial artificial chromosome (BAC) contigs from the sex-determining region of the platyfish Xiphophorus maculatus

Despite the major importance of sex determination in aquaculture, no master sex-determining gene has been identified so far in teleost fish. In the platyfish Xiphophorus maculatus, this master gene is flanked by two receptor tyrosine kinase genes, the Xmrk oncogene responsible for melanoma formation in some Xiphophorus interspecific hybrids, and its proto-oncogenic counterpart. Both Xmrk genes, which have already been characterised at the molecular level, delimit a region of about 1 Mb that contains other gene loci involved in sexual maturity, pigmentation and melanoma formation. We have constructed a genomic bacterial artificial chromosome (BAC) library of X. maculatus with a tenfold coverage of the haploid genome and walked on both X and Y sex chromosomes starting from both Xmrk genes. This led to the assembly of BAC contigs from the sex-determining region covering approximately 950 kb of the X and 750 kb of the Y chromosome. To our knowledge, these are the largest contigs reported so far for sex chromosomes in fish. Molecular analysis suggests that the sex-determining region of X. maculatus frequently undergoes retrotranspositions and other kinds of rearrangements. This genomic plasticity might be related to the high genetic variability observed in Xiphophorus for sex determination, sexual maturity, pigmentation and melanoma formation, which are encoded by gene loci located in the sex-determining region.

Estrogen-dependent growth of a rat pituitary tumor involves, but does not require, a high level of vascular endothelial growth factor

Long-term (10-week) treatment of Fischer 344 (F344) rats with the synthetic estrogen diethylstilbestrol (DES) increases the level of vascular endothelial growth factor (VEGF) in the pituitary. This is concurrent with the development of a large tumor of the pituitary of F344 rats. A role for VEGF in estrogen-dependent pituitary tumor growth is also supported by the fact that pituitary VEGF level is not increased by estrogen treatment in rats of the tumor-resistant Brown Norway (BN) strain. However, VEGF is not increased by estrogen treatment in an F(1) hybrid of F344 and BN, even though F(1) hybrid rats do form pituitary tumors in response to estrogen. Quantitative trait locus (QLT) mapping reveals that control of estrogen-dependent VEGF expression is linked to the Edpm5 QTL, which was previously identified as a QTL for estrogen-dependent pituitary tumor growth. In contrast, the QTL Edpm2-1 and Edpm9-2, which have been shown to each have a significant effect on estrogen-dependent pituitary mass of a magnitude similar to Edpm5, do not have any effect on VEGF level. Taken together, our results support the association of VEGF expression with growth of the estrogen-induced rat pituitary tumor, as has been reported by others, but they also indicate that there is significant pathways of growth regulation that are independent of high-level VEGF expression.

Zebrafish mutants identify an essential role for laminins in notochord formation

Basement membranes are thought to be essential for organ formation, providing the scaffold on which individual cells organize to form complex tissues. Laminins are integral components of basement membranes. To understand the development of a simple vertebrate organ, we have used positional cloning to characterize grumpy and sleepy, two zebrafish loci known to control notochord formation, and find that they encode laminin β1 and laminin γ1, respectively. Removal of either chain results in the dramatic loss of laminin 1 staining throughout the embryo and prevents formation of the basement membrane surrounding the notochord. Notochord cells fail to differentiate and many die by apoptosis. By transplantation, we demonstrate that, for both grumpy and sleepy, notochord differentiation can be rescued by exogenous sources of the missing laminin chain, although notochordal sources are also sufficient for rescue. These results demonstrate a clear in vivo requirement for laminin β1 and laminin γ1 in the formation of a specific vertebrate organ and show that laminin or the laminin-dependent basement membrane is essential for the differentiation of chordamesoderm to notochord.

Targeting mammary epithelial cells using a bacterial artificial chromosome

We describe the generation of transgenic mouse lines expressing Cre recombinase in epithelial cells of the lactating mammary gland. As an expression vector, we used a P1-derived bacterial artificial chromosome (PAC) which harbors the gene for the secretory milk protein, whey acidic protein (Wap). Using homologous recombination in E. coli, the PAC was modified to carry the improved coding sequence of Cre recombinase (iCre). Transgenic lines carrying the WAPiCre PAC express Cre recombinase efficiently in the majority of mammary epithelial cells upon lactation. Of only four transgenic lines produced, three express Cre recombinase to a high efficiency. LoxP-flanked DNA sequences are recombined in virtually all epithelial cells of WAPiCre transgenic mice at lactation day 3.

Genomic organisation of the approximately 1.5 Mb Smith-Magenis syndrome critical interval: transcription map, genomic contig, and candidate gene analysis

Smith-Magenis syndrome (SMS) is a multiple congenital anomalies/mental retardation syndrome associated with an interstitial deletion of chromosome 17 involving band p11.2. SMS is hypothesised to be a contiguous gene syndrome in which the phenotype arises from the haploinsufficiency of multiple, functionally-unrelated genes in close physical proximity, although the true molecular basis of SMS is not yet known. In this study, we have generated the first overlapping and contiguous transcription map of the SMS critical interval, linking the proximal 17p11.2 region near the SMS-REPM and the distal region near D17S740 in a minimum tiling path of 16 BACs and two PACs. Additional clones provide greater coverage throughout the critical region. Not including the repetitive sequences that flank the critical interval, the map is comprised of 13 known genes, 14 ESTs, and six genomic markers, and is a synthesis of Southern hybridisation and polymerase chain reaction data from gene and marker localisation to BACs and PACs and database sequence analysis from the human genome project high-throughput draft sequence. In order to identify possible candidate genes, we performed sequence analysis and determined the tissue expression pattern analysis of 10 novel ESTs that are deleted in all SMS patients. We also present a detailed review of six promising candidate genes that map to the SMS critical region.

Brain-derived neurotrophic factor gene organization and transcription in the zebrafish embryo

The gene encoding zebrafish brain-derived neurotrophic factor (BDNF) was cloned from a PAC genomic DNA library. The entire transcription unit was contained in two independently isolated clones that together encompass 120 kb of genomic DNA. The intron/exon organization of the zebrafish gene was found to be identical to that of the mammalian gene but only one promoter has so far been identified. The associated 5' exon is 67% identical to exon 1c of the rat BDNF gene. A search of the 5' flank of the cloned promoter for sequence similarities with known transcription factor binding sites revealed potential AP-1, CREB, and SP1 binding sites. Fusion constructs containing the cloned promoter and 1.7 kb of 5' flank and an enhanced green fluorescent protein reporter that becomes membrane-anchored were injected into 1-8 cell stage embryos. Expression was seen in notochord, muscle, epithelial and endothelial cells of the 1-day-old embryo in consonance with the endogenous gene. These results demonstrate that the cloned promoter mediates cell-specific expression.

Identification of a novel immunoreceptor tyrosine-based activation motif-containing molecule, STAM2, by mass spectrometry and its involvement in growth factor and cytokine receptor signaling pathways

In an effort to clone novel tyrosine-phosphorylated substrates of the epidermal growth factor receptor, we have initiated an approach coupling affinity purification using anti-phosphotyrosine antibodies to mass spectrometry-based identification. Here, we report the identification of a signaling molecule containing a Src homology 3 domain as well as an immunoreceptor tyrosine-based activation motif (ITAM). This molecule is 55% identical to a previously isolated molecule designated signal transducing adaptor molecule (STAM) that was identified as an interleukin (IL)-2-induced phosphoprotein and is therefore designated STAM2. Tyrosine phosphorylation of STAM2 is induced by growth factors such as epidermal growth factor and platelet-derived growth factor as well as by cytokines like IL-3. Several of the deletion mutants tested except the one containing only the amino-terminal region underwent tyrosine phosphorylation upon growth factor stimulation, implying that STAM2 is phosphorylated on several tyrosine residues. STAM2 is downstream of the Jak family of kinases since coexpression of STAM2 with Jak1 or Jak2 but not an unrelated Tec family kinase, Etk, resulted in its tyrosine phosphorylation. In contrast to epidermal growth factor receptor-induced phosphorylation, this required the ITAM domain since mutants lacking this region did not undergo tyrosine phosphorylation. Finally, overexpression of wild type STAM2 led to an increase in IL-2-mediated induction of c-Myc promoter activation indicating that it potentiates cytokine receptor signaling.

Rapid cloning of metK encoding methionine adenosyltransferase from Corynebacterium glutamicum by screening a genomic library on a high density colony-array

The genes SAM1 and SAM2 encoding the two different methionine adenosyltransferases (EC 2.5.1.6) in Saccharomyces cerevisiae were used as templates to generate specific DNA-probes. This heterologous mixture of DNA-probes was hybridized under low stringency hybridization conditions to a Corynebacterium glutamicum colony-array representing the complete genome. Subsequently, one genomic fragment was isolated which contained the C. glutamicum methionine adenosyltransferase gene metK (1.224 kb). When overproduced in Escherichia coli, MetK (44.2 kDa) of C. glutamicum had methionine adenosyltransferase activity. In addition, overexpression of metK in C. glutamicum led to an increased intracellular S-adenosylmethionine concentration. The metK transcript was detected by reverse transcription PCR in C. glutamicum cells in the exponential growth phase but not in the stationary phase.

Identification of a novel human sterol-sensitive ATP-binding cassette transporter (ABCA7)

We report the identification of the full-length cDNA for a novel ATP-binding cassette (ABC) transporter from human macrophages. The mRNA is of 6.8 kb size and contains an open reading frame encoding a polypeptide of 2146 amino acids with a calculated molecular weight of 220 kDa. The predicted protein product is composed of two transmembrane domains and two nucleotide binding folds indicating that it pertains to the group of full-size ABC transporters. The novel transporter shows highest protein sequence homology with the recently cloned human cholesterol and phospholipid exporter ABCA1 (54%) and the human retinal transporter ABCR (49%), both members of the ABC transporter subfamily A. In accordance with the currently proposed classification, the novel transporter was designated ABCA7. ABCA7 mRNA was detected predominantly in myelo-lymphatic tissues with highest expression in peripheral leukocytes, thymus, spleen, and bone marrow. Expression of ABCA7 is induced during in vitro differentiation of human monocytes into macrophages. In macrophages, both the ABCA7 mRNA and protein expression are upregulated in the presence of modified low density lipoprotein and downregulated by HDL(3). Our results suggest a role for ABCA7 in macrophage transmembrane lipid transport.

Mushroom body defect, a gene involved in the control of neuroblast proliferation in Drosophila, encodes a coiled-coil protein

Neurogenesis relies on the establishment of the proper number and precisely controlled proliferation of neuroblasts, the neuronal precursor cells. A role for the mushroom body defect (mud) gene in both of these aspects of neuroblast behavior, as well as possible roles in other aspects of fruit fly biology, is implied by phenotypes associated with mud mutations. We have localized mud by determining the sequence change in one point mutant, identifying a predicted ORF affected by the mutation, and showing that an appropriate segment of the genome rescues mud mutant phenotypes. An analysis of mud cDNAs and a survey of mud transcripts by Northern blotting indicate that the gene is subject to differential splicing and is expressed primarily during embryogenesis but also, at lower levels, during subsequent developmental stages in a sexually dimorphic manner. The gene is predicted to encode a polypeptide without obvious homologs but with two prominent structural features, a long coiled coil that constitutes the central core of the protein and a carboxyl-terminal transmembrane domain.

Observing the living genome

Dynamic pictures of living genomes are now beginning to emerge from systematic studies of gene expression patterns using DNA microarrays. The rich information represented in the variation in each gene's expression provides the basis for a new kind of genomic map.