Distinct epigenetic signatures delineate transcriptional programs during virus-specific CD8(+) T cell differentiation.

The molecular mechanisms that regulate the rapid transcriptional changes that occur during cytotoxic T lymphocyte (CTL) proliferation and differentiation in response to infection are poorly understood. We have utilized ChIP-seq to assess histone H3 methylation dynamics within naive, effector, and memory virus-specific T cells isolated directly ex vivo after influenza A virus infection. Our results show that within naive T cells, codeposition of the permissive H3K4me3 and repressive H3K27me3 modifications is a signature of gene loci associated with gene transcription, replication, and cellular differentiation. Upon differentiation into effector and/or memory CTLs, the majority of these gene loci lose repressive H3K27me3 while retaining the permissive H3K4me3 modification. In contrast, immune-related effector gene promoters within naive T cells lacked the permissive H3K4me3 modification, with acquisition of this modification occurring upon differentiation into effector/memory CTLs. Thus, coordinate transcriptional regulation of CTL genes with related functions is achieved via distinct epigenetic mechanisms.

A positive feedback between p53 and miR-34 miRNAs mediates tumor suppression.

As bona fide p53 transcriptional targets, miR-34 microRNAs (miRNAs) exhibit frequent alterations in many human tumor types and elicit multiple p53 downstream effects upon overexpression. Unexpectedly, miR-34 deletion alone fails to impair multiple p53-mediated tumor suppressor effects in mice, possibly due to the considerable redundancy in the p53 pathway. Here, we demonstrate that miR-34a represses HDM4, a potent negative regulator of p53, creating a positive feedback loop acting on p53. In a Kras-induced mouse lung cancer model, miR-34a deficiency alone does not exhibit a strong oncogenic effect. However, miR-34a deficiency strongly promotes tumorigenesis when p53 is haploinsufficient, suggesting that the defective p53-miR-34 feedback loop can enhance oncogenesis in a specific context. The importance of the p53/miR-34/HDM4 feedback loop is further confirmed by an inverse correlation between miR-34 and full-length HDM4 in human lung adenocarcinomas. In addition, human lung adenocarcinomas generate an elevated level of a short HDM4 isoform through alternative polyadenylation. This short HDM4 isoform lacks miR-34-binding sites in the 3' untranslated region (UTR), thereby evading miR-34 regulation to disable the p53-miR-34 positive feedback. Taken together, our results elucidated the intricate cross-talk between p53 and miR-34 miRNAs and revealed an important tumor suppressor effect generated by this positive feedback loop.

Colon tumour secretopeptidome: insights into endogenous proteolytic cleavage events in the colon tumour microenvironment.

The secretopeptidome comprises endogenous peptides derived from proteins secreted into the tumour microenvironment through classical and non-classical secretion. This study characterised the low-Mr (<3kDa) component of the human colon tumour (LIM1215, LIM1863) secretopeptidome, as a first step towards gaining insights into extracellular proteolytic cleavage events in the tumour microenvironment. Based on two biological replicates, this secretopeptidome isolation strategy utilised differential centrifugal ultrafiltration in combination with analytical RP-HPLC and nanoLC-MS/MS. Secreted peptides were identified using a combination of Mascot and post-processing analyses including MSPro re-scoring, extended feature sets and Percolator, resulting in 474 protein identifications from 1228 peptides (≤1% q-value, ≤5% PEP) - a 36% increase in peptide identifications when compared with conventional Mascot (homology ionscore thresholding). In both colon tumour models, 122 identified peptides were derived from 41 cell surface protein ectodomains, 23 peptides (12 proteins) from regulated intramembrane proteolysis (RIP), and 12 peptides (9 proteins) generated from intracellular domain proteolysis. Further analyses using the protease/substrate database MEROPS, (http://merops.sanger.ac.uk/), revealed 335 (71%) proteins classified as originating from classical/non-classical secretion, or the cell membrane. Of these, peptides were identified from 42 substrates in MEROPS with defined protease cleavage sites, while peptides generated from a further 205 substrates were fragmented by hitherto unknown proteases. A salient finding was the identification of peptides from 88 classical/non-classical secreted substrates in MEROPS, implicated in tumour progression and angiogenesis (FGFBP1, PLXDC2), cell-cell recognition and signalling (DDR1, GPA33), and tumour invasiveness and metastasis (MACC1, SMAGP); the nature of the proteases responsible for these proteolytic events is unknown. To confirm reproducibility of peptide fragment abundance in this study, we report the identification of a specific cleaved peptide fragment in the secretopeptidome from the colon-specific GPA33 antigen in 4/14 human CRC models. This improved secretopeptidome isolation and characterisation strategy has extended our understanding of endogenous peptides generated through proteolysis of classical/non-classical secreted proteins, extracellular proteolytic processing of cell surface membrane proteins, and peptides generated through RIP. The novel peptide cleavage site information in this study provides a useful first step in detailing proteolytic cleavage associated with tumourigenesis and the extracellular environment. This article is part of a Special Issue entitled: An Updated Secretome.

G protein-linked signaling pathways in bipolar and major depressive disorders.

The G-protein linked signaling system (GPLS) comprises a large number of G-proteins, G protein-coupled receptors (GPCRs), GPCR ligands, and downstream effector molecules. G-proteins interact with both GPCRs and downstream effectors such as cyclic adenosine monophosphate (cAMP), phosphatidylinositols, and ion channels. The GPLS is implicated in the pathophysiology and pharmacology of both major depressive disorder (MDD) and bipolar disorder (BPD). This study evaluated whether GPLS is altered at the transcript level. The gene expression in the dorsolateral prefrontal (DLPFC) and anterior cingulate (ACC) were compared from MDD, BPD, and control subjects using Affymetrix Gene Chips and real time quantitative PCR. High quality brain tissue was used in the study to control for confounding effects of agonal events, tissue pH, RNA integrity, gender, and age. GPLS signaling transcripts were altered especially in the ACC of BPD and MDD subjects. Transcript levels of molecules which repress cAMP activity were increased in BPD and decreased in MDD. Two orphan GPCRs, GPRC5B and GPR37, showed significantly decreased expression levels in MDD, and significantly increased expression levels in BPD. Our results suggest opposite changes in BPD and MDD in the GPLS, "activated" cAMP signaling activity in BPD and "blunted" cAMP signaling activity in MDD. GPRC5B and GPR37 both appear to have behavioral effects, and are also candidate genes for neurodegenerative disorders. In the context of the opposite changes observed in BPD and MDD, these GPCRs warrant further study of their brain effects.

Copy number variation in patients with disorders of sex development due to 46,XY gonadal dysgenesis.

Disorders of sex development (DSD), ranging in severity from mild genital abnormalities to complete sex reversal, represent a major concern for patients and their families. DSD are often due to disruption of the genetic programs that regulate gonad development. Although some genes have been identified in these developmental pathways, the causative mutations have not been identified in more than 50% 46,XY DSD cases. We used the Affymetrix Genome-Wide Human SNP Array 6.0 to analyse copy number variation in 23 individuals with unexplained 46,XY DSD due to gonadal dysgenesis (GD). Here we describe three discrete changes in copy number that are the likely cause of the GD. Firstly, we identified a large duplication on the X chromosome that included DAX1 (NR0B1). Secondly, we identified a rearrangement that appears to affect a novel gonad-specific regulatory region in a known testis gene, SOX9. Surprisingly this patient lacked any signs of campomelic dysplasia, suggesting that the deletion affected expression of SOX9 only in the gonad. Functional analysis of potential SRY binding sites within this deleted region identified five putative enhancers, suggesting that sequences additional to the known SRY-binding TES enhancer influence human testis-specific SOX9 expression. Thirdly, we identified a small deletion immediately downstream of GATA4, supporting a role for GATA4 in gonad development in humans. These CNV analyses give new insights into the pathways involved in human gonad development and dysfunction, and suggest that rearrangements of non-coding sequences disturbing gene regulation may account for significant proportion of DSD cases.

Conference Scene: Lessons learned from the 5th Statistical Analysis Workshop of the Pharmacogenetics Research Network.

The Pharmacogenomics Research Network holds a statistical analysis workshop every other year to share novel statistical methods and study designs for pharmacogenomics research, as well as insightful analyses of substantive ongoing studies. The 5th workshop was held 15 April 2009, in Rochester (MN, USA), in conjunction with the general Pharmacogenomics Research Network meeting. This summary of the ten contributed talks highlights a variety of timely topics, including identification of functional variants, how to maximize power using various study designs, and pathway analysis approaches. We also discuss the keynote invited presentation by Terry Speed, which provided an overview of statistical issues with next-generation sequence data with an emphasis on some statistical challenges in mRNA sequence data. Novel applications of Poisson regression models demonstrated innovative, yet practical, approaches to distinguish between biological and technical sources of variation in the counts of mRNA transcript reads. Overall, the workshop emphasized the need for diverse approaches to conducting pharmacogenomics studies, as well as the evolving nature of the field.

Genomic profiling and expression studies reveal both positive and negative activities for the Drosophila Myb MuvB/dREAM complex in proliferating cells.

Myb-MuvB (MMB)/dREAM is a nine-subunit complex first described in Drosophila as a repressor of transcription, dependent on E2F2 and the RBFs. Myb, an integral member of MMB, curiously plays no role in the silencing of the test genes previously analyzed. Moreover, Myb plays an activating role in DNA replication in Drosophila egg chamber follicle cells. The essential functions for Myb are executed as part of MMB. This duality of function lead to the hypothesis that MMB, which contains both known activator and repressor proteins, might function as part of a switching mechanism that is dependent on DNA sites and developmental context. Here, we used proliferating Drosophila Kc tissue culture cells to explore both the network of genes regulated by MMB (employing RNA interference and microarray expression analysis) and the genomic locations of MMB following chromatin immunoprecipitation (ChIP) and tiling array analysis. MMB occupied 3538 chromosomal sites and was promoter-proximal to 32% of Drosophila genes. MMB contains multiple DNA-binding factors, and the data highlighted the combinatorial way by which the complex was targeted and utilized for regulation. Interestingly, only a subset of chromatin-bound complexes repressed genes normally expressed in a wide range of developmental pathways. At many of these sites, E2F2 was critical for repression, whereas at other nonoverlapping sites, Myb was critical for repression. We also found sites where MMB was a positive regulator of transcript levels that included genes required for mitotic functions (G2/M), which may explain some of the chromosome instability phenotypes attributed to loss of Myb function in myb mutants.

A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes.

Recent studies suggest that thousands of genes may contribute to breast cancer pathophysiologies when deregulated by genomic or epigenomic events. Here, we describe a model "system" to appraise the functional contributions of these genes to breast cancer subsets. In general, the recurrent genomic and transcriptional characteristics of 51 breast cancer cell lines mirror those of 145 primary breast tumors, although some significant differences are documented. The cell lines that comprise the system also exhibit the substantial genomic, transcriptional, and biological heterogeneity found in primary tumors. We show, using Trastuzumab (Herceptin) monotherapy as an example, that the system can be used to identify molecular features that predict or indicate response to targeted therapies or other physiological perturbations.

Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data.

Genome-wide expression profiling is a powerful tool for implicating novel gene ensembles in cellular mechanisms of health and disease. The most popular platform for genome-wide expression profiling is the Affymetrix GeneChip. However, its selection of probes relied on earlier genome and transcriptome annotation which is significantly different from current knowledge. The resultant informatics problems have a profound impact on analysis and interpretation the data. Here, we address these critical issues and offer a solution. We identified several classes of problems at the individual probe level in the existing annotation, under the assumption that current genome and transcriptome databases are more accurate than those used for GeneChip design. We then reorganized probes on more than a dozen popular GeneChips into gene-, transcript- and exon-specific probe sets in light of up-to-date genome, cDNA/EST clustering and single nucleotide polymorphism information. Comparing analysis results between the original and the redefined probe sets reveals approximately 30-50% discrepancy in the genes previously identified as differentially expressed, regardless of analysis method. Our results demonstrate that the original Affymetrix probe set definitions are inaccurate, and many conclusions derived from past GeneChip analyses may be significantly flawed. It will be beneficial to re-analyze existing GeneChip data with updated probe set definitions.

Colon cancer prognosis prediction by gene expression profiling.

This study assessed the possibility to build a prognosis predictor, based on microarray gene expression measures, in stage II and III colon cancer patients. Tumour (T) and non-neoplastic mucosa (NM) mRNA samples from 18 patients (nine with a recurrence, nine with no recurrence) were profiled using the Affymetrix HGU133A GeneChip. The k-nearest neighbour method was used for prognosis prediction using T and NM gene expression measures. Six-fold cross-validation was applied to select the number of neighbours and the number of informative genes to include in the predictors. Based on this information, one T-based and one NM-based predictor were proposed and their accuracies were estimated by double cross-validation. In six-fold cross-validation, the lowest numbers of informative genes giving the lowest numbers of false predictions (two out of 18) were 30 and 70 with the T and NM gene expression measures, respectively. A 30-gene T-based predictor and a 70-gene NM-based predictor were then built, with estimated accuracies of 78 and 83%, respectively. This study suggests that one can build an accurate prognosis predictor for stage II and III colon cancer patients, based on gene expression measures, and one can use either tumour or non-neoplastic mucosa for this purpose.

Evolution of the relaxin-like peptide family: from neuropeptide to reproduction.

The relaxin-like peptide family consists of relaxin-1, relaxin-2, and relaxin-3 and the insulin-like peptides (INSL)-3, INSL4, INSL5, and INSL6 (human relaxin-2 is equivalent to relaxin-1 in other species). Evolution of this family has been contentious. We therefore sought to clarify the issue by performing phylogenetic analysis of all relaxin-like peptides from the genomic databases available. Surprisingly, the phylogeny, combined with previous biologic characterizations, suggest that although relaxin's original function was likely in the brain, its reproductive role was acquired just prior to the divergence of amphibians. This phylogeny also illuminates inconsistencies in relaxin evolution in invertebrates, chickens, and cows.

Rooting a phylogenetic tree with nonreversible substitution models.

BACKGROUND: We compared two methods of rooting a phylogenetic tree: the stationary and the nonstationary substitution processes. These methods do not require an outgroup. METHODS: Given a multiple alignment and an unrooted tree, the maximum likelihood estimates of branch lengths and substitution parameters for each associated rooted tree are found; rooted trees are compared using their likelihood values. Site variation in substitution rates is handled by assigning sites into several classes before the analysis. RESULTS: In three test datasets where the trees are small and the roots are assumed known, the nonstationary process gets the correct estimate significantly more often, and fits data much better, than the stationary process. Both processes give biologically plausible root placements in a set of nine primate mitochondrial DNA sequences. CONCLUSIONS: The nonstationary process is simple to use and is much better than the stationary process at inferring the root. It could be useful for situations where an outgroup is unavailable.

Normalization of cDNA microarray data.

Normalization means to adjust microarray data for effects which arise from variation in the technology rather than from biological differences between the RNA samples or between the printed probes. This paper describes normalization methods based on the fact that dye balance typically varies with spot intensity and with spatial position on the array. Print-tip loess normalization provides a well-tested general purpose normalization method which has given good results on a wide range of arrays. The method may be refined by using quality weights for individual spots. The method is best combined with diagnostic plots of the data which display the spatial and intensity trends. When diagnostic plots show that biases still remain in the data after normalization, further normalization steps such as plate-order normalization or scale-normalization between the arrays may be undertaken. Composite normalization may be used when control spots are available which are known to be not differentially expressed. Variations on loess normalization include global loess normalization and two-dimensional normalization. Detailed commands are given to implement the normalization techniques using freely available software.

Genes for glycosylphosphatidylinositol toxin biosynthesis in Plasmodium falciparum.

About 2.5 million people die of Plasmodium falciparum malaria every year. Fatalities are associated with systemic and organ-specific inflammation initiated by a parasite toxin. Recent studies show that glycosylphosphatidylinositol (GPI) functions as the dominant parasite toxin in the context of infection. GPIs also serve as membrane anchors for several of the most important surface antigens of parasite invasive stages. GPI anchoring is a complex posttranslational modification produced through the coordinated action of a multicomponent biosynthetic pathway. Here we present eight new genes of P. falciparum selected for encoding homologs of proteins essential for GPI synthesis: PIG-A, PIG-B, PIG-M, PIG-O, GPI1, GPI8, GAA-1, and DPM1. We describe the experimentally verified mRNA and predicted amino acid sequences and in situ localization of the gene products to the parasite endoplasmic reticulum. Moreover, we show preliminary evidence for the PIG-L and PIG-C genes. The biosynthetic pathway of the malaria parasite GPI offers potential targets for drug development and may be useful for studying parasite cell biology and the molecular basis for the pathophysiology of parasitic diseases.

An HMM model for coiled-coil domains and a comparison with PSSM-based predictions.

MOTIVATION: Large-scale sequence data require methods for the automated annotation of protein domains. Many of the predictive methods are based either on a Position Specific Scoring Matrix (PSSM) of fixed length or on a window-less Hidden Markov Model (HMM). The performance of the two approaches is tested for Coiled-Coil Domains (CCDs). The prediction of CCDs is used frequently, and its optimization seems worthwhile. RESULTS: We have conceived MARCOIL, an HMM for the recognition of proteins with a CCD on a genomic scale. A cross-validated study suggests that MARCOIL improves predictions compared to the traditional PSSM algorithm, especially for some protein families and for short CCDs. The study was designed to reveal differences inherent in the two methods. Potential confounding factors such as differences in the dimension of parameter space and in the parameter values were avoided by using the same amino acid propensities and by keeping the transition probabilities of the HMM constant during cross-validation. AVAILABILTY: The prediction program and the databases are available at http://www.wehi.edu.au/bioweb/Mauro/Marcoil