A Mixture model with random-effects components for clustering correlated gene-expression profiles

Motivation: The clustering of gene profiles across some experimental conditions of interest contributes significantly to the elucidation of unknown gene function, the validation of gene discoveries and the interpretation of biological processes. However, this clustering problem is not straightforward as the profiles of the genes are not all independently distributed and the expression levels may have been obtained from an experimental design involving replicated arrays. Ignoring the dependence between the gene profiles and the structure of the replicated data can result in important sources of variability in the experiments being overlooked in the analysis, with the consequent possibility of misleading inferences being made. We propose a random-effects model that provides a unified approach to the clustering of genes with correlated expression levels measured in a wide variety of experimental situations. Our model is an extension of the normal mixture model to account for the correlations between the gene profiles and to enable covariate information to be incorporated into the clustering process. Hence the model is applicable to longitudinal studies with or without replication, for example, time-course experiments by using time as a covariate, and to cross-sectional experiments by using categorical covariates to represent the different experimental classes. Results: We show that our random-effects model can be fitted by maximum likelihood via the EM algorithm for which the E(expectation) and M(maximization) steps can be implemented in closed form. Hence our model can be fitted deterministically without the need for time-consuming Monte Carlo approximations. The effectiveness of our model-based procedure for the clustering of correlated gene profiles is demonstrated on three real datasets, representing typical microarray experimental designs, covering time-course, repeated-measurement and cross-sectional data. In these examples, relevant clusters of the genes are obtained, which are supported by existing gene-function annotation. A synthetic dataset is considered too.

Identification of novel "pathologs" (human disease-related gene candidates) from the RIKEN full-length mouse cDNA data set

The majority of common diseases such as cancer, allergy, diabetes, or heart disease are characterized by complex genetic traits, in which genetic and environmental components contribute to disease susceptibility. Our knowledge of the genetic factors underlying most of such diseases is limited. A major goal in the post-genomic era is to identify and characterize disease susceptibility genes and to use this knowledge for disease treatment and prevention. More than 500 genes are conserved across the invertebrate and vertebrate genomes. Because of gene conservation, various organisms including yeast, fruitfly, zebrafish, rat, and mouse have been used as genetic models.

Growth factors and cytokines modulate gene expression of cell-surface proteoglycans in human periodontal ligament cells

Cell-surface proteoglycans have been known to be involved in many functions including interactions with components of the extracellular microenvironment, and act as co-receptors which bind and modify the action of various growth factors and cytokines. The purpose of this study was to determine the regulation by growth factors and cytokines on cell-surface proteoglycan gene expression in cultured human periodontal ligament (PDL) cells. Subconfluent, quiescent PDL cells were treated with various concentrations of serum, bFGF, PDGF-BB, TGF-beta1, IL-1 beta, and IFN-gamma. RT-PCR technique was used, complemented with Northern blot for syndecan-1, to examine the effects of these agents on the mRNA expression of five cell-surface proteoglycans (syndecan-1, syndecan-2, syndecan-4, glypican and betaglycan). Syndecan-1 mRNA levels increased in response to serum, bFGF and PDGF-BB, but decreased in response to TGF-beta1, IL-1 beta and IFN-gamma. In contrast, syndecan-2 mRNA levels were upregulated by TCF-beta1 and IL-1 beta stimulation, but remained unchanged with the other agents. Betaglycan gene expression decreased in response to serum, but was upregulated by TCF-beta1 and unchanged by the other stimulants. Additionally, syndecan-4 and glypican were not significantly altered in response to the regulator molecules studied, with the exception that glypican is decreased in response to IFN-gamma. These data demonstrate that the gene expression of the five cell-surface proteoglycans studied is differentially regulated in PDL cells lending support to the nation of distinct functions for these cell-surface proteoglycans. (C) 2001 Wiley-Liss, inc.

The evolution of controlled multitasked gene networks: The role of introns and other noncoding RNAs in the development of complex organisms

Eukaryotic phenotypic diversity arises from multitasking of a core proteome of limited size. Multitasking is routine in computers, as well as in other sophisticated information systems, and requires multiple inputs and outputs to control and integrate network activity. Higher eukaryotes have a mosaic gene structure with a dual output, mRNA (protein-coding) sequences and introns, which are released from the pre-mRNA by posttranscriptional processing. Introns have been enormously successful as a class of sequences and comprise up to 95% of the primary transcripts of protein-coding genes in mammals. In addition, many other transcripts (perhaps more than half) do not encode proteins at all, but appear both to be developmentally regulated and to have genetic function. We suggest that these RNAs (eRNAs) have evolved to function as endogenous network control molecules which enable direct gene-gene communication and multitasking of eukaryotic genomes. Analysis of a range of complex genetic phenomena in which RNA is involved or implicated, including co-suppression, transgene silencing, RNA interference, imprinting, methylation, and transvection, suggests that a higher-order regulatory system based on RNA signals operates in the higher eukaryotes and involves chromatin remodeling as well as other RNA-DNA, RNA-RNA, and RNA-protein interactions. The evolution of densely connected gene networks would be expected to result in a relatively stable core proteome due to the multiple reuse of components, implying,that cellular differentiation and phenotypic variation in the higher eukaryotes results primarily from variation in the control architecture. Thus, network integration and multitasking using trans-acting RNA molecules produced in parallel with protein-coding sequences may underpin both the evolution of developmentally sophisticated multicellular organisms and the rapid expansion of phenotypic complexity into uncontested environments such as those initiated in the Cambrian radiation and those seen after major extinction events.

Using clusters of computers for large QU-GENE simulation experiments

The QU-GENE Computing Cluster (QCC) is a hardware and software solution to the automation and speedup of large QU-GENE (QUantitative GENEtics) simulation experiments that are designed to examine the properties of genetic models, particularly those that involve factorial combinations of treatment levels. QCC automates the management of the distribution of components of the simulation experiments among the networked single-processor computers to achieve the speedup.

SOX9 enhances aggrecan gene promoter/enhancer activity and is up-regulated by retinoic acid in a cartilage-derived cell line, TC6

SOX9 is a transcription factor that plays a key role in chondrogenesis, Aggrecan is one of the major structural components in cartilage; however, the molecular mechanism of aggrecan gene regulation has not yet been fully elucidated, TC6 is a clonal chondrocytic cell line derived from articular cartilage, The purpose of this study was to examine whether SOX9 modulates aggrecan gene expression and to further identify molecules that regulate Sox9 expression in TC6 cells. SOX9 overexpression in TC6 cells enhanced by similar to 3-fold the transcriptional activity of the AgCAT-8 construct containing S-kilobase (kb) promoter/first exon/first intron fragments of the aggrecan gene. SOX9 enhancement of aggrecan promoter activity was lost when we deleted a 4.5-kb fragment from the 3'-end of the 8-kb fragment corresponding to the region including the first intron, In TC6 cells, SOX9 enhanced the transcriptional activity of a reporter construct containing the Sry/Sox consensus sequence >10-fold. SOX9 enhancement of aggrecan gene promoter activity and SOX9 transactivation through the Sry/Sox consensus sequence were not observed in osteoblastic osteosarcoma cells (ROS17/2.8), indicating the dependence on the cellular background. Northern blot analysis indicated that TC6 cells constitutively express Sox9 mRNA at relatively low levels. To examine regulation of Sox9 gene expression, we investigated the effects of calciotropic hormones and cytokines, Among these, retinoic acid (RA) specifically enhanced Sox9 mRNA expression in TC6 cells. The basal levels of Sox9 expression and its enhancement by RA were observed similarly at both permissive (33 degrees C) and nonpermissive (39 degrees C) temperatures. Furthermore, RA treatment enhanced the transcriptional activity of a reporter construct containing the Sry/Sox consensus sequence in TC6 cells. Moreover, RA treatment also enhanced the transcriptional activity of another reporter construct containing the enhancer region of the type II procollagen gene in TC6 cells. These observations indicate that SOX9 enhances aggrecan promoter activity and that its expression is up-regulated by RA in TC6 cells.

Development of disabled, replication-defective gene transfer vectors from the Jembrana disease virus, a new infectious agent of cattle

Jembrana disease virus (JDV) is a newly isolated and characterised bovine lentivirus. It causes an acute disease in Ball cattle (Bos javanicus). which can be readily transmitted to susceptible cattle with 17% mortality. There is as yet no treatment or preventive vaccine. We have developed a gene transfer vector system based on JDV that has three components. The first of the components is a bicistronic transfer vector plasmid that was constructed to contain cis-sequences from the JDV genome, including 5 '- and 3 ' -long terminal repeats (LTRs), 0.4 kb of truncated gag and 1.1 kb of 3 ' -env, a multiple cloning site to accommodate the gene(s) of interest for transfer, and an internal ribosome entry site plus the neomycin phosphotransferase (Neo) gene cassette for antibiotic selection. The second element is a packaging plasmid that contains trans-sequences. including gag, pol. vif, tar and rev: but without the env and packaging signals. The third is a plasmid encoding the G glycoprotein of vesicular stomatitis virus (VSV-G) to supply the vector an envelope for pseudotyping. Cotransfection of 293T cells with these three plasmid components produced VSV-G pseudotyped. disabled, replication defective, bicistronic JDV vectors encoding the green fluorescent protein (EGFP) and the Neo resistance selection maker simultaneously with a titre range of (0.4-1.2) x 10(6) CFU/ml. Transduction of several replicating primary and transformed cells from cattle, primate and human sources and importantly growth-arrested cells with the JDV vectors showed high efficiency of EGFP gene transfer at 35-75%, which was stable and the expression of EGFP was long term. Furthermore, these JDV vectors were designed to suit the inclusion and expression of genes corresponding to JDV specific proteins, such as gag or env, for the development of vaccines for Jembrana disease. This strategy should also be applicable to other bovine diseases as wall. The design and construction of the JDV vector system should facilitate the study of the lentivirology and pathogenesis of the diseases associated with JDV or other bovine virus infections. To our knowledge, this is the first such vector system developed from a cattle virus. (C) 2001 Elsevier Science B.V. All rights reserved.

Cloning and functional analysis of the Sry-related HMG box gene, Sox18

The Sox gene family (Sry like HMG box gene) is characterised by a conserved DNA sequence encoding a domain of approximately 80 amino acids which is responsible for sequence specific DNA binding. We initially published the identification and partial cDNA sequence of murine Sox18, a new member of this gene family, isolated from a cardiac cDNA library. This sequence allowed us to classify Sox18 into the F sub-group of Sox proteins, along with Sox7 and Sox17. Recently, we demonstrated that mutations in the Sox18 activation domain underlie cardiovascular and hair follicle defects in the mouse mutation, ragged (Ra) (Pennisi et al., 2000. Mutations in Sox18 underlie cardiovascular and hair follicle defecs in ragged mice. Nat. Genet. 24, 434-437). Ra homozygotes lack vibrissae and coat hairs, have generalised oedema and an accumulation of chyle in the peritoneum. Here we have investigated the genomic sequences encoding Sox18. Screening of a mouse genomic phage library identified four overlapping clones, we sequenced a 3.25 kb XbaI fragment that defined the entire coding region and approximately 1.5 kb of 5' flanking sequences. This identified (i) an additional 91 amino acids upstream of the previously designated methionine start codon in the original cDNA, and (ii);ln intron encoded within the HMG box/DNA binding domain in exactly the same position as that found in the Sox5, -13 and -17 genes. The Sox18 gene encodes a protein of 468 aa. We present evidence that suggests HAF-2, the human HMG-box activating factor-2 protein, is the orthologue of murine Sox18. HAF-2 has been implicated in the regulation of the Human IgH enhancer in a B cell context. Random mutagenesis coupled with GAL4 hybrid analysis in the activation domain between amino acids 252 and 346, of Sox18, implicated the phosphorylation motif, SARS, and the region between amino acid residues 313 and 346 as critical components of Sox18 mediated transactivation. Finally, we examined the expression of Sox18 in multiple adult mouse tissues using RT-PCR. Low-moderate expression was observed in spleen, stomach, kidney, intestine, skeletal muscle and heart. Very abundant expression was detected in lung tissue. (C) 2001 Elsevier Science B.V. All rights reserved.

Gene-Environment Interaction Effects on Behavioral Variation and Risk of Complex Disorders: The Example of Alcoholism and Other Psychiatric Disorders

There have been few replicated examples of genotype x environment interaction effects on behavioral variation or risk of psychiatric disorder. We review some of the factors that have made detection of genotype x environment interaction effects difficult, and show how genotype x shared environment interaction (GxSE) effects are commonly confounded with genetic parameters in data from twin pairs reared together. Historic data on twin pairs reared apart can in principle be used to estimate such GxSE effects, but have rarely been used for this purpose. We illustrate this using previously published data from the Swedish Adoption Twin Study of Aging (SATSA), which suggest that GxSE effects could account for as much as 25% of the total variance in risk of becoming a regular smoker. Since few separated twin pairs will be available for study in the future, we also consider methods for modifying variance components linkage analysis to allow for environmental interactions with linked loci.

Analysis of Melanoma Onset: Assessing Familial Aggregation by Using Estimating Equations and Fitting Variance Components via Bayesian Random Effects Models

We investigate whether relative contributions of genetic and shared environmental factors are associated with an increased risk in melanoma. Data from the Queensland Familial Melanoma Project comprising 15,907 subjects arising from 1912 families were analyzed to estimate the additive genetic, common and unique environmental contributions to variation in the age at onset of melanoma. Two complementary approaches for analyzing correlated time-to-onset family data were considered: the generalized estimating equations (GEE) method in which one can estimate relationship-specific dependence simultaneously with regression coefficients that describe the average population response to changing covariates; and a subject-specific Bayesian mixed model in which heterogeneity in regression parameters is explicitly modeled and the different components of variation may be estimated directly. The proportional hazards and Weibull models were utilized, as both produce natural frameworks for estimating relative risks while adjusting for simultaneous effects of other covariates. A simple Markov Chain Monte Carlo method for covariate imputation of missing data was used and the actual implementation of the Bayesian model was based on Gibbs sampling using the free ware package BUGS. In addition, we also used a Bayesian model to investigate the relative contribution of genetic and environmental effects on the expression of naevi and freckles, which are known risk factors for melanoma.

Antagonistic interaction between abscisic acid and jasmonate-ethylene signaling pathways modulates defense gene expression and disease resistance in Arabidopsis

The plant hormones abscisic acid (ABA), jasmonic acid (JA), and ethylene are involved in diverse plant processes, including the regulation of gene expression during adaptive responses to abiotic and biotic stresses. Previously, ABA has been implicated in enhancing disease susceptibility in various plant species, but currently very little is known about the molecular mechanisms underlying this phenomenon. In this study, we obtained evidence that a complex interplay between ABA and JA-ethylene signaling pathways regulate plant defense gene expression and disease resistance. First, we showed that exogenous ABA suppressed both basal and JA-ethylene-activated transcription from defense genes. By contrast, ABA deficiency as conditioned by the mutations in the ABA1 and ABA2 genes, which encode enzymes involved in ABA biosynthesis, resulted in upregulation of basal and induced transcription from JA-ethylene responsive defense genes. Second, we found that disruption of AtMYC2 (allelic to JASMONATE INSENSITIVE1 [JIN1]), encoding a basic helix-loop-helix Leu zipper transcription factor, which is a positive regulator of ABA signaling, results in elevated levels of basal and activated transcription from JA-ethylene responsive defense genes. Furthermore, the jin1/myc2 and aba2-1 mutants showed increased resistance to the necrotrophic fungal pathogen Fusarium oxysporum. Finally, using ethylene and ABA signaling mutants, we showed that interaction between ABA and ethylene signaling is mutually antagonistic in vegetative tissues. Collectively, our results indicate that the antagonistic interactions between multiple components of ABA and the JA-ethylene signaling pathways modulate defense and stress responsive gene expression in response to biotic and abiotic stresses.

Respiratory gene clusters of Metallosphaera sedula - differential expression and transcriptional organization

Metallosphaera sedula is a thermoacidophilic Crenarchaeon which is capable of leaching metals from sulfidic ores. The authors have investigated the presence and expression of genes encoding respiratory complexes in this organism when grown heterotrophically or chemolithotrophically on either sulfur or pyrite. The presence of three gene clusters, encoding two terminal oxidase complexes, the quinol oxidase SoxABCD and the SoxM oxidase supercomplex, and a gene cluster encoding a high-potential cytochrome b and components of a bc(1) complex analogue (cbsBA-soxL2N gene cluster) was established. Expression studies showed that the soxM gene was expressed to high levels during heterotrophic growth of M. sedula on yeast extract, while the soxABCD mRNA was most abundant in cells grown on sulfur. Reduced-minus-oxidized difference spectra of cell membranes showed cytochrome-related peaks that correspond to published spectra of Sulfolobus-type terminal oxidase complexes. In pyrite-grown cells, expression levels of the two monitored oxidase gene clusters were reduced by a factor of 10-12 relative to maximal expression levels, although spectra of membranes clearly contained oxidase-associated haems, suggesting the presence of additional gene clusters encoding terminal oxidases in M. sedula. Pyrite- and sulfur-grown cells contained high levels of the cbsA transcript, which encodes a membrane-bound cytochrome b with a possible role in iron oxidation or chemolithotrophy. The cbsA gene is not co-transcribed with the soxL2N genes, and therefore does not appear to be an integral part of this bc(1) complex analogue. The data show for the first time the differential expression of the Sulfolobus-type terminal oxidase gene clusters in a Crenarchaeon in response to changing growth modes.

Butyrylcholinesterase: Association with the metabolic syndrome and identification of 2 gene loci affecting activity

Background: Plasma cholinesterase activity is known to be correlated with plasma triglycerides, HDL- and LDL-cholesterol, and other features of the metabolic syndrome. A role in triglyceride metabolism has been proposed. Genetic variants that decrease activity have been studied extensively, but the factors contributing to overall variation in the population are poorly understood. We studied plasma cholinesterase activity in a sample of 2200 adult twins to assess covariation with cardiovascular risk factors and components of the metabolic syndrome, to determine the degree of genetic effects on enzyme activity, and to search for quantitative trait loci affecting activity. Methods and Results: Cholinesterase activity was lower in women than in men before the age of 50, but increased to activity values similar to those in males after that age. There were highly significant correlations with variables associated with the metabolic syndrome: plasma triglyceride, HDL- and LDL-cholesterol, apolipoprotein B and E, urate, and insulin concentrations; gamma-glutamyltransferase and aspartate and alanine aminotransferase activities; body mass index; and blood pressure. The heritability of plasma cholinesterase activity was 65%. Linkage analysis with data from the dizygotic twin pairs showed suggestive linkage on chromosome 3 at the location of the cholinesterase WHO gene and also on chromosome 5. Conclusions: Our results confirm and extend the connection between cholinesterase, cardiovascular risk factors, and metabolic syndrome. They establish a substantial heritability for plasma cholinesterase activity that might be attributable to variation near the structural gene and at an independent locus. (c) 2006 American Association for Clinical Chemistry.

Variance components analyses of multiple asthma traits in a large sample of Australian families ascertained through a twin proband

Background: Intermediate phenotypes are often measured as a proxy for asthma. It is largely unclear to what extent the same set of environmental or genetic factors regulate these traits. Objective: Estimate the environmental and genetic correlations between self-reported and clinical asthma traits. Methods: A total of 3073 subjects from 802 families were ascertained through a twin proband. Traits measured included self-reported asthma, airway histamine responsiveness (AHR), skin prick response to common allergens including house dust mite (Dermatophagoides pteronyssinus [D. pter]), baseline lung function, total serum immunoglobulin E (IgE) and eosinophilia. Bivariate and multivariate analyses of eight traits were performed with adjustment for ascertainment and significant covariates. Results: Overall 2716 participants completed an asthma questionnaire and 2087 were clinically tested, including 1289 self-reported asthmatics (92% previously diagnosed by a doctor). Asthma, AHR, markers of allergic sensitization and eosinophilia had significant environmental correlations with each other (range: 0.23-0.89). Baseline forced expiratory volume in 1 s (FEV1) showed low environmental correlations with most traits. Fewer genetic correlations were significantly different from zero. Phenotypes with greatest genetic similarity were asthma and atopy (0.46), IgE and eosinophilia (0.44), AHR and D. pter (0.43) and AHR and airway obstruction (-0.43). Traits with greatest genetic dissimilarity were FEV1 and atopy (0.05), airway obstruction and IgE (0.07) and FEV1 and D. pter (0.11). Conclusion: These results suggest that the same set of environmental factors regulates the variation of many asthma traits. In addition, although most traits are regulated to great extent by specific genetic factors, there is still some degree of genetic overlap that could be exploited by multivariate linkage approaches.

Fate of hairpin transcript components during RNA silencing and its suppression in transgenic virus-resistant tobacco

Transgenic tobacco plants, carrying a Potato virus Y (PVY)-NIa hairpin sequence separated by a unique unrelated spacer sequence were specifically silenced and highly resistant to PVY infection. In such plants neither PVY-NIa nor spacer transgene transcripts were detectable by specific quantitative real time reverse transcriptase PCR (RT-qPCR) assays of similar relative efficiencies developed for direct comparative analysis. However, small interfering RNAs (siRNAs) specific for the PVY sequence of the transgene and none specific for the LNYV spacer sequence were detected. Following infection with Cucumber mosaic virus (CMV), which suppresses dsRNA-induced RNA silencing, transcript levels of PVY-NIa as well as spacer sequence increased manifold with the same time course. The cellular abundance of the single-stranded (ss) spacer sequence was consistently higher than that of PVY dsRNA in all cases. The results show that during RNA silencing and its suppression of a hairpin transcript in transgenic tobacco, the ssRNA spacer sequence is affected differently than the dsRNA. In PVY-silenced plants. the spacer is efficiently degraded by a mechanism not involving the accumulation of siRNAs, while following suppression of RNA silencing by CMV, the spacer appears protected from degradation. Crown Copyright (c) 2006 Published by Elsevier B.V. All rights reserved.