Optimising the analysis of transcript data using high density oligonucleotide arrays and genomic DNA-based probe selection

Background: Affymetrix GeneChip arrays are widely used for transcriptomic studies in a diverse range of species. Each gene is represented on a GeneChip array by a probe- set, consisting of up to 16 probe-pairs. Signal intensities across probe- pairs within a probe-set vary in part due to different physical hybridisation characteristics of individual probes with their target labelled transcripts. We have previously developed a technique to study the transcriptomes of heterologous species based on hybridising genomic DNA (gDNA) to a GeneChip array designed for a different species, and subsequently using only those probes with good homology. Results: Here we have investigated the effects of hybridising homologous species gDNA to study the transcriptomes of species for which the arrays have been designed. Genomic DNA from Arabidopsis thaliana and rice (Oryza sativa) were hybridised to the Affymetrix Arabidopsis ATH1 and Rice Genome GeneChip arrays respectively. Probe selection based on gDNA hybridisation intensity increased the number of genes identified as significantly differentially expressed in two published studies of Arabidopsis development, and optimised the analysis of technical replicates obtained from pooled samples of RNA from rice. Conclusion: This mixed physical and bioinformatics approach can be used to optimise estimates of gene expression when using GeneChip arrays.

A Brassica exon array for whole-transcript gene expression profiling

Affymetrix GeneChip (R) arrays are used widely to study transcriptional changes in response to developmental and environmental stimuli. GeneChip (R) arrays comprise multiple 25-mer oligonucleotide probes per gene and retain certain advantages over direct sequencing. For plants, there are several public GeneChip (R) arrays whose probes are localised primarily in 39 exons. Plant whole-transcript (WT) GeneChip (R) arrays are not yet publicly available, although WT resolution is needed to study complex crop genomes such as Brassica, which are typified by segmental duplications containing paralogous genes and/or allopolyploidy. Available sequence data were sampled from the Brassica A and C genomes, and 142,997 gene models identified. The assembled gene models were then used to establish a comprehensive public WT exon array for transcriptomics studies. The Affymetrix GeneChip (R) Brassica Exon 1.0 ST Array is a 5 mu M feature size array, containing 2.4 million 25-base oligonucleotide probes representing 135,201 gene models, with 15 probes per gene distributed among exons. Discrimination of the gene models was based on an E-value cut-off of 1E(-5), with <= 98 sequence identity. The 135 k Brassica Exon Array was validated by quantifying transcriptome differences between leaf and root tissue from a reference Brassica rapa line (R-o-18), and categorisation by Gene Ontologies (GO) based on gene orthology with Arabidopsis thaliana. Technical validation involved comparison of the exon array with a 60-mer array platform using the same starting RNA samples. The 135 k Brassica Exon Array is a robust platform. All data relating to the array design and probe identities are available in the public domain and are curated within the BrassEnsembl genome viewer at http://www.brassica.info/BrassEnsembl/index.html.

A genome wide association study of mathematical ability reveals an association at chromosome 3q29, a locus associated with autism and learning difficulties: a preliminary study

Mathematical ability is heritable, but few studies have directly investigated its molecular genetic basis. Here we aimed to identify specific genetic contributions to variation in mathematical ability. We carried out a genome wide association scan using pooled DNA in two groups of U.K. samples, based on end of secondary/high school national academic exam achievement: high (n = 419) versus low (n = 183) mathematical ability while controlling for their verbal ability. Significant differences in allele frequencies between these groups were searched for in 906,600 SNPs using the Affymetrix GeneChip Human Mapping version 6.0 array. After meeting a threshold of p<1.5×10-5, 12 SNPs from the pooled association analysis were individually genotyped in 542 of the participants and analyzed to validate the initial associations (lowest p-value 1.14 ×10-6). In this analysis, one of the SNPs (rs789859) showed significant association after Bonferroni correction, and four (rs10873824, rs4144887, rs12130910 rs2809115) were nominally significant (lowest p-value 3.278 × 10-4). Three of the SNPs of interest are located within, or near to, known genes (FAM43A, SFT2D1, C14orf64). The SNP that showed the strongest association, rs789859, is located in a region on chromosome 3q29 that has been previously linked to learning difficulties and autism. rs789859 lies 1.3 kbp downstream of LSG1, and 700 bp upstream of FAM43A, mapping within the potential promoter/regulatory region of the latter. To our knowledge, this is only the second study to investigate the association of genetic variants with mathematical ability, and it highlights a number of interesting markers for future study.

A pooled genome-wide association study of Asperger Syndrome

Asperger Syndrome (AS) is a neurodevelopmental condition characterized by impairments in social interaction and communication, alongside the presence of unusually repetitive, restricted interests and stereotyped behaviour. Individuals with AS have no delay in cognitive and language development. It is a subset of Autism Spectrum Conditions (ASC), which are highly heritable and has a population prevalence of approximately 1%. Few studies have investigated the genetic basis of AS. To address this gap in the literature, we performed a genome-wide pooled DNA association study to identify candidate loci in 612 individuals (294 cases and 318 controls) of Caucasian ancestry, using the Affymetrix GeneChip Human Mapping version 6.0 array. We identified 11 SNPs that had a p-value below 1x10-5. These SNPs were independently genotyped in the same sample. Three of the SNPs (rs1268055, rs7785891 and rs2782448) were nominally significant, though none remained significant after Bonferroni correction. Two of our top three SNPs (rs7785891 and rs2782448) lie in loci previously implicated in ASC. However, investigation of the three SNPs in the ASC genome-wide association dataset from the Psychiatric Genomics Consortium indicated that these three SNPs were not significantly associated with ASC. The effect sizes of the variants were modest, indicating that our study was not sufficiently powered to identify causal variants with precision.

Transcriptome analysis of grain development in hexaploid wheat

Background: Hexaploid wheat is one of the most important cereal crops for human nutrition. Molecular understanding of the biology of the developing grain will assist the improvement of yield and quality traits for different environments. High quality transcriptomics is a powerful method to increase this understanding. Results: The transcriptome of developing caryopses from hexaploid wheat ( Triticum aestivum, cv. Hereward) was determined using Affymetrix wheat GeneChip (R) oligonucleotide arrays which have probes for 55,052 transcripts. Of these, 14,550 showed significant differential regulation in the period between 6 and 42 days after anthesis ( daa). Large changes in transcript abundance were observed which were categorised into distinct phases of differentiation ( 6 - 10 daa), grain fill ( 12 - 21 daa) and desiccation/maturation ( 28 - 42 daa) and were associated with specific tissues and processes. A similar experiment on developing caryopses grown with dry and/or hot environmental treatments was also analysed, using the profiles established in the first experiment to show that most environmental treatment effects on transcription were due to acceleration of development, but that a few transcripts were specifically affected. Transcript abundance profiles in both experiments for nine selected known and putative wheat transcription factors were independently confirmed by real time RT-PCR. These expression profiles confirm or extend our knowledge of the roles of the known transcription factors and suggest roles for the unknown ones. Conclusion: This transcriptome data will provide a valuable resource for molecular studies on wheat grain. It has been demonstrated how it can be used to distinguish general developmental shifts from specific effects of treatments on gene expression and to diagnose the probable tissue specificity and role of transcription factors.

Viral nucleolar localisation signals determine dynamic trafficking within the nucleolus

Localisation of both viral and cellular proteins to the nucleolus is determined by a variety of factors including nucleolar localisation signals (NoLSs), but how these signals operate is not clearly understood. The nucleolar trafficking of wild type viral proteins and chimeric proteins, which contain altered NoLSs, were compared to investigate the role of NoLSs in dynamic nucleolar trafficking. Three viral proteins from diverse viruses were selected which localised to the nucleolus; the coronavirus infectious bronchitis virus nucleocapsid (N) protein, the herpesvirus saimiri ORF57 protein and the HIV-1 Rev protein. The chimeric proteins were N protein and ORF57 protein which had their own NoLS replaced with those from ORF57 and Rev proteins, respectively. By analysing the sub-cellular localisation and trafficking of these viral proteins and their chimeras within and between nucleoli using confocal microscopy and photo-bleaching we show that NoLSs are responsible for different nucleolar localisations and trafficking rates.

Endothelin-converting enzyme-1 regulates trafficking and signalling of the neurokinin 1 receptor in endosomes of myenteric neurones

Neuropeptide signalling at the plasma membrane is terminated by neuropeptide degradation by cell-surface peptidases, and by beta-arrestin-dependent receptor desensitization and endocytosis. However, receptors continue to signal from endosomes by beta-arrestin-dependent processes, and endosomal sorting mediates recycling and resensitization of plasma membrane signalling. The mechanisms that control signalling and trafficking of receptors in endosomes are poorly defined. We report a major role for endothelin-converting enzyme-1 (ECE-1) in controlling substance P (SP) and the neurokinin 1 receptor (NK(1)R) in endosomes of myenteric neurones. ECE-1 mRNA and protein were expressed by myenteric neurones of rat and mouse intestine. SP (10 nM, 10 min) induced interaction of NK(1)R and beta-arrestin at the plasma membrane, and the SP-NK(1)R-beta-arrestin signalosome complex trafficked by a dynamin-mediated mechanism to ECE-1-containing early endosomes, where ECE-1 can degrade SP. After 120 min, NK(1)R recycled from endosomes to the plasma membrane. ECE-1 inhibitors (SM-19712, PD-069185) and the vacuolar H(+)ATPase inhibitor bafilomycin A(1), which prevent endosomal SP degradation, suppressed NK(1)R recycling by >50%. Preincubation of neurones with SP (10 nM, 5 min) desensitized Ca(2+) transients to a second SP challenge after 10 min, and SP signals resensitized after 60 min. SM-19712 inhibited NK(1)R resensitization by >90%. ECE-1 inhibitors also caused sustained SP-induced activation of extracellular signal-regulated kinases, consistent with stabilization of the SP-NK(1)R-beta-arrestin signalosome. By degrading SP and destabilizing endosomal signalosomes, ECE-1 has a dual role in controlling endocytic signalling and trafficking of the NK(1)R: promoting resensitization of G protein-mediated plasma membrane signalling, and terminating beta-arrestin-mediated endosomal signalling.

Trafficking and signaling of G protein-coupled receptors in the nervous system: implications for disease and therapy

G protein-coupled receptors (GPCRs) are expressed throughout the nervous system where they regulate multiple physiological processes, participate in neurological diseases, and are major targets for therapy. Given that many GPCRs respond to neurotransmitters and hormones that are present in the extracellular fluid and which do not readily cross the plasma membrane, receptor trafficking to and from the plasma membrane is a critically important determinant of cellular responsiveness. Moreover, trafficking of GPCRs throughout the endosomal system can initiate signaling events that are mechanistically and functionally distinct from those operating at the plasma membrane. This review discusses recent advances in the relationship between signaling and trafficking of GPCRs in the nervous system. It summarizes how receptor modifications influence trafficking, discusses mechanisms that regulate GPCR trafficking to and from the plasma membrane, reviews the relationship between trafficking and signaling, and considers the implications of GPCR trafficking to drug development.

Hepatocyte growth factor-regulated tyrosine kinase substrate (HRS) mediates post-endocytic trafficking of protease-activated receptor 2 and calcitonin receptor-like receptor

The E3 ligase c-Cbl ubiquitinates protease-activated receptor 2 (PAR(2)), which is required for post-endocytic sorting of PAR(2) to lysosomes, where degradation arrests signaling. The mechanisms of post-endocytic sorting of ubiquitinated receptors are incompletely understood. Here, we investigated the role of hepatocyte growth factor-regulated tyrosine kinase substrate (HRS), in post-endocytic sorting and signaling of PAR(2). In HEK-PAR(2) cells, PAR(2) activating peptide (PAR(2)-AP) induced PAR(2) trafficking from the cell surface to early endosomes containing endogenous HRS, and then to lysosomes. HRS overexpression or knockdown with small interfering RNA caused formation of enlarged HRS-positive endosomes, where activated PAR(2) and c-Cbl accumulated, and PAR(2) failed to traffic to lysosomes. Overexpression of HRS prevented PAR(2)-AP-induced degradation of PAR(2), as determined by Western blotting. Overexpression of HRS mutant lacking an ubiquitin-binding motif similarly caused retention of PAR(2) in enlarged endosomes. Moreover, HRS overexpression or knockdown caused retention of ubiquitin-resistant PAR(2)Delta14K/R in enlarged HRS-containing endosomes, preventing recycling and resensitization of PAR(2)Delta14K/R. HRS overexpression or knockdown similarly prevented lysosomal trafficking and recycling of calcitonin receptor-like receptor, a non-ubiquitinated receptor that traffics to lysosomes after sustained activation and recycles after transient activation. Thus, HRS plays a critically important role in the post-endocytic sorting of single receptors, PAR(2) and CLR, to both degradative and recycling pathways. This sorting role for HRS is independent of its ubiquitin-interacting motif, and it can regulate trafficking of both ubiquitinated and non-ubiquitinated PAR(2) and non-ubiquitinated CLR. The ultimate sorting decision to degradative or recycling pathways appears to occur downstream from HRS.

A comparison of the Thlaspi caerulescens and Thlaspi arvense shoot transcriptomes

Whole-genome transcriptome profiling is revealing how biological systems are regulated at the transcriptional level. This study reports the development of a robust method to profile and compare the transcriptomes of two nonmodel plant species, Thlaspi caerulescens, a zinc (Zn) hyperaccumulator, and Thlaspi arvense, a nonhyperaccumulator, using Affymetrix Arabidopsis thaliana ATH1-121501 GeneChip (R) arrays (Affymetrix, Santa Clara, CA, USA). Transcript abundance was quantified in the shoots of agar- and compost-grown plants of both species. Analyses were optimized using a genomic DNA (gDNA)-based probe-selection strategy based on the hybridization efficiency of Thlaspi gDNA with corresponding A. thaliana probes. In silico alignments of GeneChip (R) probes with Thlaspi gene sequences, and quantitative real-time PCR, confirmed the validity of this approach. Approximately 5000 genes were differentially expressed in the shoots of T. caerulescens compared with T. arvense, including genes involved in Zn transport and compartmentalization. Future functional analyses of genes identified as differentially expressed in the shoots of these closely related species will improve our understanding of the molecular mechanisms of Zn hyperaccumulation.

Using genomic DNA-based probe-selection to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species

High-density oligonucleotide (oligo) arrays are a powerful tool for transcript profiling. Arrays based on GeneChip® technology are amongst the most widely used, although GeneChip® arrays are currently available for only a small number of plant and animal species. Thus, we have developed a method to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species and tested the method by analysing the transcriptome of Brassica oleracea L., a species for which no GeneChip® array is available, using a GeneChip® array designed for Arabidopsis thaliana (L.) Heynh. Genomic DNA from B. oleracea was labelled and hybridised to the ATH1-121501 GeneChip® array. Arabidopsis thaliana probe-pairs that hybridised to the B. oleracea genomic DNA on the basis of the perfect-match (PM) probe signal were then selected for subsequent B. oleracea transcriptome analysis using a .cel file parser script to generate probe mask files. The transcriptional response of B. oleracea to a mineral nutrient (phosphorus; P) stress was quantified using probe mask files generated for a wide range of gDNA hybridisation intensity thresholds. An example probe mask file generated with a gDNA hybridisation intensity threshold of 400 removed > 68 % of the available PM probes from the analysis but retained >96 % of available A. thaliana probe-sets. Ninety-nine of these genes were then identified as significantly regulated under P stress in B. oleracea, including the homologues of P stress responsive genes in A. thaliana. Increasing the gDNA hybridisation intensity thresholds up to 500 for probe-selection increased the sensitivity of the GeneChip® array to detect regulation of gene expression in B. oleracea under P stress by up to 13-fold. Our open-source software to create probe mask files is freely available http://affymetrix.arabidopsis.info/xspecies/ webcite and may be used to facilitate transcriptomic analyses of a wide range of plant and animal species in the absence of custom arrays.

Trafficking of storage proteins in developing grain of wheat

The processing properties of the wheat flour are largely determined by the structures and interactions of the grain storage proteins (also called gluten proteins) which form a continuous visco-elastic network in dough. Wheat gluten proteins are classically divided into two groups, the monomeric gliadins and the polymeric glutenins, with the latter being further classified into low molecular weight (LMW) and high molecular weight (HMW) subunits. The synthesis, folding and deposition of the gluten proteins take place within the endomembrane system of the plant cell. However, determination of the precise routes of trafficking and deposition of individual gluten proteins in developing wheat grain has been limited in the past by the difficulty of developing monospecific antibodies. To overcome this limitation, a single gluten protein (a LMW subunit) was expressed in transgenic wheat with a C-terminal epitope tag, allowing the protein to be located in the cells of the developing grain using highly specific antibodies. This approach was also combined with the use of wider specificity antibodies to compare the trafficking and deposition of different gluten protein groups within the same endosperm cells. These studies are in agreement with previous suggestions that two trafficking pathways occur in wheat, with the proteins either being transported via the Golgi apparatus into the vacuole or accumulating directly within the lumen of the ER. They also suggest that the same individual protein could be trafficked by either pathway, possibly depending on the stage of development, and that segregation of gluten proteins both between and within protein bodies may occur.

Targeting stem cell niches and trafficking for cardiovascular therapy.

Regenerative cardiovascular medicine is the frontline of 21st-century health care. Cell therapy trials using bone marrow progenitor cells documented that the approach is feasible, safe and potentially beneficial in patients with ischemic disease. However, cardiovascular prevention and rehabilitation strategies should aim to conserve the pristine healing capacity of a healthy organism as well as reactivate it under disease conditions. This requires an increased understanding of stem cell microenvironment and trafficking mechanisms. Engagement and disengagement of stem cells of the osteoblastic niche is a dynamic process, finely tuned to allow low amounts of cells move out of the bone marrow and into the circulation on a regular basis. The balance is altered under stress situations, like tissue injury or ischemia, leading to remarkably increased cell egression. Individual populations of circulating progenitor cells could give rise to mature tissue cells (e.g. endothelial cells or cardiomyocytes), while the majority may differentiate to leukocytes, affecting the environment of homing sites in a paracrine way, e.g. promoting endothelial survival, proliferation and function, as well as attenuating or enhancing inflammation. This review focuses on the dynamics of the stem cell niche in healthy and disease conditions and on therapeutic means to direct stem cell/progenitor cell mobilization and recruitment into improved tissue repair.

The trafficking pathway of a wheat storage protein in transgenic rice endosperm

Background and Aims The trafficking of proteins in the endoplasmic reticulum (ER) of plant cells is a topic of considerable interest since this organelle serves as an entry point for proteins destined for other organelles, as well as for the ER itself. In the current work, transgenic rice was used to study the pattern and pathway of deposition of the wheat high molecular weight (HMW) glutenin sub-unit (GS) 1Dx5 within the rice endosperm using specific antibodies to determine whether it is deposited in the same or different protein bodies from the rice storage proteins, and whether it is located in the same or separate phases within these. Methods The protein distribution and the expression pattern of HMW sub-unit 1Dx5 in transgenic rice endosperm at different stages of development were determined using light and electron microscopy after labelling with antibodies. Key results The use of HMW-GS-specific antibodies showed that sub-unit 1Dx5 was expressed mainly in the sub-aleurone cells of the endosperm and that it was deposited in both types of protein body present in the rice endosperm: derived from the ER and containing prolamins, and derived from the vacuole and containing glutelins. In addition, new types of protein bodies were also formed within the endosperm cells. Conclusions The results suggest that the HMW 1Dx5 protein could be trafficked by either the ER or vacuolar pathway, possibly depending on the stage of development, and that its accumulation in the rice endosperm could compromise the structural integrity of protein bodies and their segregation into two distinct populations in the mature endosperm.

Intracellular trafficking, localization, and mobilization of platelet-borne thiol isomerases

OBJECTIVE: Thiol isomerases facilitate protein folding in the endoplasmic reticulum, and several of these enzymes, including protein disulfide isomerase and ERp57, are mobilized to the surface of activated platelets, where they influence platelet aggregation, blood coagulation, and thrombus formation. In this study, we examined the synthesis and trafficking of thiol isomerases in megakaryocytes, determined their subcellular localization in platelets, and identified the cellular events responsible for their movement to the platelet surface on activation. APPROACH AND RESULTS: Immunofluorescence microscopy imaging was used to localize protein disulfide isomerase and ERp57 in murine and human megakaryocytes at various developmental stages. Immunofluorescence microscopy and subcellular fractionation analysis were used to localize these proteins in platelets to a compartment distinct from known secretory vesicles that overlaps with an inner cell-surface membrane region defined by the endoplasmic/sarcoplasmic reticulum proteins calnexin and sarco/endoplasmic reticulum calcium ATPase 3. Immunofluorescence microscopy and flow cytometry were used to monitor thiol isomerase mobilization in activated platelets in the presence and absence of actin polymerization (inhibited by latrunculin) and in the presence or absence of membrane fusion mediated by Munc13-4 (absent in platelets from Unc13dJinx mice). CONCLUSIONS: Platelet-borne thiol isomerases are trafficked independently of secretory granule contents in megakaryocytes and become concentrated in a subcellular compartment near the inner surface of the platelet outer membrane corresponding to the sarco/endoplasmic reticulum of these cells. Thiol isomerases are mobilized to the surface of activated platelets via a process that requires actin polymerization but not soluble N-ethylmaleimide-sensitive fusion protein attachment receptor/Munc13-4-dependent vesicular-plasma membrane fusion.