Oviposition by pierid butterflies triggers defense responses in Arabidopsis.

Insect eggs represent a threat for the plant as hatching larvae rapidly start with their feeding activity. Using a whole-genome microarray, we studied the expression profile of Arabidopsis (Arabidopsis thaliana) leaves after oviposition by two pierid butterflies. For Pieris brassicae, the deposition of egg batches changed the expression of hundreds of genes over a period of 3 d after oviposition. The transcript signature was similar to that observed during a hypersensitive response or in lesion-mimic mutants, including the induction of defense and stress-related genes and the repression of genes involved in growth and photosynthesis. Deposition of single eggs by Pieris rapae caused a similar although much weaker transcriptional response. Analysis of the jasmonic acid and salicylic acid mutants coi1-1 and sid2-1 indicated that the response to egg deposition is mostly independent of these signaling pathways. Histochemical analyses showed that egg deposition is causing a localized cell death, accompanied by the accumulation of callose, and the production of reactive oxygen species. In addition, activation of the pathogenesis-related1::beta-glucuronidase reporter gene correlated precisely with the site of egg deposition and was also triggered by crude egg extract. This study provides molecular evidence for the detection of egg deposition by Arabidopsis plants and suggests that oviposition causes a localized response with strong similarity to a hypersensitive response.

The transcriptome of the arbuscular mycorrhizal fungus Glomus intraradices (DAOM 197198) reveals functional tradeoffs in an obligate symbiont.

? The arbuscular mycorrhizal symbiosis is arguably the most ecologically important eukaryotic symbiosis, yet it is poorly understood at the molecular level. To provide novel insights into the molecular basis of symbiosis-associated traits, we report the first genome-wide analysis of the transcriptome from Glomus intraradices DAOM 197198. ? We generated a set of 25,906 nonredundant virtual transcripts (NRVTs) transcribed in germinated spores, extraradical mycelium and symbiotic roots using Sanger and 454 sequencing. NRVTs were used to construct an oligoarray for investigating gene expression. ? We identified transcripts coding for the meiotic recombination machinery, as well as meiosis-specific proteins, suggesting that the lack of a known sexual cycle in G. intraradices is not a result of major deletions of genes essential for sexual reproduction and meiosis. Induced expression of genes encoding membrane transporters and small secreted proteins in intraradical mycelium, together with the lack of expression of hydrolytic enzymes acting on plant cell wall polysaccharides, are all features of G. intraradices that are shared with ectomycorrhizal symbionts and obligate biotrophic pathogens. ? Our results illuminate the genetic basis of symbiosis-related traits of the most ancient lineage of plant biotrophs, advancing future research on these agriculturally and ecologically important symbionts.

Increased EBV Reactivation in the Cerebrospinal Fluid of Patients with Early Multiple Sclerosis

Epstein-Barr virus (EBV) has been consistently associated with multiple sclerosis (MS), but whether this virus is a trigger of MS remains undetermined. Recently, EBV-infected B cells recognized by activated CD8_ T cells have been detected in the meninges of autopsied MS patients. In addition, a strong EBV-specific CD8_ T cell response in the blood of patients with MS of recent onset was reported. Here, to further explore the putative relationship between MS and EBV, we assessed the EBV-specific cellular and humoral immune responses in the blood and the cerebrospinal fluid (CSF) of patients with early MS or other neurological diseases, separated into inflammatory (IOND) and non-inflammatory (NIOND) groups. The MS non-associated neurotropic herpesvirus cytomegalovirus (CMV) served as a control. Fifty-eight study subjects were enrolled, including 44 patients (13 with early MS (onset of MS less than one year prior to the assay), 15 with IOND and 16 with NIOND) in the immunological arm of the study. The cellular immune response was investigated using a functional CFSE cytotoxic T lymphocyte (CTL) assay performed with short-term cultured EBV- or CMVspecific effector T cells from the CSF and the blood. The humoral immune response specific for these two viruses was also examined in both the blood and the CSF. The recruitment of a given virusspecific antibody in the CSF as compared to the blood was expressed as antibody indexes (AI). We found that, in the CSF of early MS patients, there was an enrichment in EBV-, but not CMV-specific, CD8_ CTL as compared to the CSF of IOND (P_ 0.003) and NIOND patients (P_0.0009), as well as compared to paired blood samples (P_0.005). Additionally, relative viral capsid antigen (VCA)-, but not EBV encoded nuclear antigen 1 (EBNA1)- or CMV-specific, AI were increased in the CSF of early MS as compared to IOND (P_0.002) or NIOND patients (P_0.008) and correlated with the EBVspecific CD8_ CTL responses in the CSF (rs_0.54, P_0.001). Fourteen additional patients were enrolled in the virological arm of the study: using semi-nested PCR, EBV-encoded nuclear RNA1 (EBER1)-a transcript expressed during all stages of EBV infection-was detected in the CSF of 2/4 early MS, but only 1/6 IOND and 0/4 NIOND patients. Altogether, our data suggest that a reactivation of EBV, but not CMV, is taking place in the central nervous system of patients with MS of recent onset. These data significantly strengthen the link between EBV and MS and may indicate a triggering role of EBV in this disease. This work was supported by grants from the Swiss National Foundation and from the Swiss Society for Multiple Sclerosis.

Fusion of the human gene for the polyubiquitination co-effector UEV-1 with Kua, a newly identified gene

UEV proteins are enzymatically inactive variants of the E2 ubiquitin-conjugating enzymes that regulate noncanonical elongation of ubiquitin chains. In Saccharomyces cerevisiae, UEV is part of the RAD6-mediated error-free DNA repair pathway. In mammalian cells, UEV proteins can modulate c-FOS transcription and the G2-M transition of the cell cycle. Here we show that the UEV genes from phylogenetically distant organisms present a remarkable conservation in their exon–intron structure. We also show that the human UEV1 gene is fused with the previously unknown gene Kua. In Caenorhabditis elegans and Drosophila melanogaster, Kua and UEV are in separated loci, and are expressed as independent transcripts and proteins. In humans, Kua and UEV1 are adjacent genes, expressed either as separate transcripts encoding independent Kua and UEV1 proteins, or as a hybrid Kua–UEV transcript, encoding a two-domain protein. Kua proteins represent a novel class of conserved proteins with juxtamembrane histidine-rich motifs. Experiments with epitope-tagged proteins show that UEV1A is a nuclear protein, whereas both Kua and Kua–UEV localize to cytoplasmic structures, indicating that the Kua domain determines the cytoplasmic localization of Kua–UEV. Therefore, the addition of a Kua domain to UEV in the fused Kua–UEV protein confers new biological properties to this regulator of variant polyubiquitination.[Kua cDNAs isolated by RT-PCR and described in this paper have been deposited in the GenBank data library under accession nos. AF1155120 (H. sapiens) and AF152361 (D. melanogaster). Genomic clones containing UEV genes: S. cerevisiae, YGL087c (accession no. Z72609); S. pombe, c338 (accession no. AL023781); P. falciparum, MAL3P2 (accession no. AL034558); A. thaliana, F26F24 (accession no. AC005292); C. elegans, F39B2 (accession no. Z92834); D. melanogaster, AC014908; and H. sapiens, 1185N5 (accession no. AL034423). Accession numbers for Kua cDNAs in GenBank dbEST: M. musculus, AA7853; T. cruzi, AI612534. Other Kua-containing sequences: A. thaliana genomic clones F10M23 (accession no. AL035440), F19K23 (accession no. AC000375), and T20K9 (accession no. AC004786).

Islands of euchromatin-like sequence and expressed polymorphic sequences within the short arm of human chromosome 21

The goals of the human genome project did not include sequencing of the heterochromatic regions. We describe here an initial sequence of 1.1 Mb of the short arm of human chromosome 21 (HSA21p), estimated to be 10% of 21p. This region contains extensive euchromatic-like sequence and includes on average one transcript every 100 kb. These transcripts show multiple inter- and intrachromosomal copies, and extensive copy number and sequence variability. The sequencing of the "heterochromatic" regions of the human genome is likely to reveal many additional functional elements and provide important evolutionary information.

Prominent use of distal 5’ transcription start sites and discovery of a large number of additional exons in ENCODE regions

This report presents systematic empirical annotation of transcript products from 399 annotated protein-coding loci across the 1% of the human genome targeted by the Encyclopedia of DNA elements (ENCODE) pilot project using a combination of 5' rapid amplification of cDNA ends (RACE) and high-density resolution tiling arrays. We identified previously unannotated and often tissue- or cell-line-specific transcribed fragments (RACEfrags), both 5' distal to the annotated 5' terminus and internal to the annotated gene bounds for the vast majority (81.5%) of the tested genes. Half of the distal RACEfrags span large segments of genomic sequences away from the main portion of the coding transcript and often overlap with the upstream-annotated gene(s). Notably, at least 20% of the resultant novel transcripts have changes in their open reading frames (ORFs), most of them fusing ORFs of adjacent transcripts. A significant fraction of distal RACEfrags show expression levels comparable to those of known exons of the same locus, suggesting that they are not part of very minority splice forms. These results have significant implications concerning (1) our current understanding of the architecture of protein-coding genes; (2) our views on locations of regulatory regions in the genome; and (3) the interpretation of sequence polymorphisms mapping to regions hitherto considered to be "noncoding," ultimately relating to the identification of disease-related sequence alterations.

A general definition and nomenclature for alternative splicing events

Understanding the molecular mechanisms responsible for the regulation of the transcriptome present in eukaryotic cells isone of the most challenging tasks in the postgenomic era. In this regard, alternative splicing (AS) is a key phenomenoncontributing to the production of different mature transcripts from the same primary RNA sequence. As a plethora ofdifferent transcript forms is available in databases, a first step to uncover the biology that drives AS is to identify thedifferent types of reflected splicing variation. In this work, we present a general definition of the AS event along with anotation system that involves the relative positions of the splice sites. This nomenclature univocally and dynamically assignsa specific ‘‘AS code’’ to every possible pattern of splicing variation. On the basis of this definition and the correspondingcodes, we have developed a computational tool (AStalavista) that automatically characterizes the complete landscape of ASevents in a given transcript annotation of a genome, thus providing a platform to investigate the transcriptome diversityacross genes, chromosomes, and species. Our analysis reveals that a substantial part—in human more than a quarter—ofthe observed splicing variations are ignored in common classification pipelines. We have used AStalavista to investigate andto compare the AS landscape of different reference annotation sets in human and in other metazoan species and found thatproportions of AS events change substantially depending on the annotation protocol, species-specific attributes, andcoding constraints acting on the transcripts. The AStalavista system therefore provides a general framework to conductspecific studies investigating the occurrence, impact, and regulation of AS.

EGASP: the human ENCODE Genome Annotation Assessment Project

Background: We present the results of EGASP, a community experiment to assess the state-ofthe-art in genome annotation within the ENCODE regions, which span 1% of the human genomesequence. The experiment had two major goals: the assessment of the accuracy of computationalmethods to predict protein coding genes; and the overall assessment of the completeness of thecurrent human genome annotations as represented in the ENCODE regions. For thecomputational prediction assessment, eighteen groups contributed gene predictions. Weevaluated these submissions against each other based on a ‘reference set’ of annotationsgenerated as part of the GENCODE project. These annotations were not available to theprediction groups prior to the submission deadline, so that their predictions were blind and anexternal advisory committee could perform a fair assessment.Results: The best methods had at least one gene transcript correctly predicted for close to 70%of the annotated genes. Nevertheless, the multiple transcript accuracy, taking into accountalternative splicing, reached only approximately 40% to 50% accuracy. At the coding nucleotidelevel, the best programs reached an accuracy of 90% in both sensitivity and specificity. Programsrelying on mRNA and protein sequences were the most accurate in reproducing the manuallycurated annotations. Experimental validation shows that only a very small percentage (3.2%) of the selected 221 computationally predicted exons outside of the existing annotation could beverified.Conclusions: This is the first such experiment in human DNA, and we have followed thestandards established in a similar experiment, GASP1, in Drosophila melanogaster. We believe theresults presented here contribute to the value of ongoing large-scale annotation projects and shouldguide further experimental methods when being scaled up to the entire human genome sequence.

Differentiated evolutionary rates in alternative exons and the implications for splicing regulation

Background: Alternatively spliced exons play an important role in the diversification of gene function in most metazoans and are highly regulated by conserved motifs in exons and introns. Two contradicting properties have been associated to evolutionary conserved alternative exons: higher sequence conservation and higher rate of non-synonymous substitutions, relative to constitutive exons. In order to clarify this issue, we have performed an analysis of the evolution of alternative and constitutive exons, using a large set of protein coding exons conserved between human and mouse and taking into account the conservation of the transcript exonic structure. Further, we have also defined a measure of the variation of the arrangement of exonic splicing enhancers (ESE-conservation score) to study the evolution of splicing regulatory sequences. We have used this measure to correlate the changes in the arrangement of ESEs with the divergence of exon and intron sequences. Results: We find evidence for a relation between the lack of conservation of the exonic structure and the weakening of the sequence evolutionary constraints in alternative and constitutive exons. Exons in transcripts with non-conserved exonic structures have higher synonymous (dS) and non-synonymous (dN) substitution rates than exons in conserved structures. Moreover, alternative exons in transcripts with non-conserved exonic structure are the least constrained in sequence evolution, and at high EST-inclusion levels they are found to be very similar to constitutive exons, whereas alternative exons in transcripts with conserved exonic structure have a dS significantly lower than average at all EST-inclusion levels. We also find higher conservation in the arrangement of ESEs in constitutive exons compared to alternative ones. Additionally, the sequence conservation at flanking introns remains constant for constitutive exons at all ESE-conservation values, but increases for alternative exons at high ESE-conservation values. Conclusion: We conclude that most of the differences in dN observed between alternative and constitutive exons can be explained by the conservation of the transcript exonic structure. Low dS values are more characteristic of alternative exons with conserved exonic structure, but not of those with non-conserved exonic structure. Additionally, constitutive exons are characterized by a higher conservation in the arrangement of ESEs, and alternative exons with an ESE-conservation similar to that of constitutive exons are characterized by a conservation of the flanking intron sequences higher than average, indicating the presence of more intronic regulatory signals.

Evaluation of the chicken transcriptome by SAGE of B cells and the DT40 cell line

Background: The understanding of whole genome sequences in higher eukaryotes depends to a large degree on the reliable definition of transcription units including exon/intron structures, translated open reading frames (ORFs) and flanking untranslated regions. The best currently available chicken transcript catalog is the Ensembl build based on the mappings of a relatively small number of full length cDNAs and ESTs to the genome as well as genome sequence derived in silico gene predictions.Results: We use Long Serial Analysis of Gene Expression (LongSAGE) in bursal lymphocytes and the DT40 cell line to verify the quality and completeness of the annotated transcripts. 53.6% of the more than 38,000 unique SAGE tags (unitags) match to full length bursal cDNAs, the Ensembl transcript build or the genome sequence. The majority of all matching unitags show single matches to the genome, but no matches to the genome derived Ensembl transcript build. Nevertheless, most of these tags map close to the 3' boundaries of annotated Ensembl transcripts.Conclusions: These results suggests that rather few genes are missing in the current Ensembl chicken transcript build, but that the 3' ends of many transcripts may not have been accurately predicted. The tags with no match in the transcript sequences can now be used to improve gene predictions, pinpoint the genomic location of entirely missed transcripts and optimize the accuracy of gene finder software.

Transcriptional regulation of MDR1, a gene involved in antifungal drug resistance in Candida albicans

ABSTRACT Upregulation of the Major Facilitator transporter gene MDR1 (Multi_drug Resistance 1) is one of the mechanisms observed in Candida albicans clinical isolates developing resistance to azole antifungal agents. To better understand this phenomenon, the cis-acting regulatory elements present in a modulatable reporter system under the control of the MDR1 promoter were characterized. In an azole-susceptible strain, transcription of this reporter is transiently upregulated in response to either benomyl or H2O2, whereas its expression is constitutively high in an azole-resistant strain (FR2). Two cis-acting regulatory elements, that are necessary and sufficient to convey the same transcriptional responses to a heterologous promoter (CDR2), were identified within the MDR1promoter. The first element, called BRE (for Benomyl Response Element, -296 to -260 with respect to the ATG start codon), is required for benomyl-dependent MDR1 upregulation and for constitutive high expression of MDR1 in FR2. The second element, termed HRE (for H2O2 Response Element, -561 to -520), is required for H2O2-dependent MDR1 upregulation, but is dispensable for constitutive high expression. Two potential binding sites (TTAG/CTAA) for the blip transcription factor Cap1p lie within the HRE. Moreover, inactivation of CAP1 abolished the transient response to H2O2 and diminished significantly the transient response to benomyl. Cap1p, which has been previously implicated in cellular responses to oxidative stress, may thus play a transacting and positive regulatory role in benomyl- and H2O2-dependent transcription of MDR1. However, it is not the only transcription factor involved in the response of MDR1 to benomyl. A minimal BRE element (-290 to -273) that is sufficient to detect in vitro sequence-specific binding of protein complexes in crude extracts prepared from C. albicans was also delimited. Genome-wide transcript profiling analyses undertaken with a matched pair of clinical isolates, one of which being azole-resistant and upregulating MDR1, and with an azole-susceptible strain exposed to benomyl, revealed that genes specifically upregulated by benomyl harbour in their promoters Cap1p binding site(s). This strengthened the idea that Cap1p plays a role in benomyl-dependent upregulation of MDR1. BRE-like sequences were also identified in several genes co-regulated with MDR1 in both conditions, which was consistent with the involvement of the BRE in both processes. A set of 147 mutants lacking a single transcription factor gene was next screened for loss of MDR1response to benomyl. Unfortunately, none of the tested mutants showed a loss of benomyl-dependent MDR1 upregulation. Nevertheless, a significant diminution of the response was observed in the mutants in which the MADS-box transcription factor Mcm1p and the C2H2 zinc finger transcription factor orf19.13374p were inactivated, suggesting that Mcm1p and orf19.13374p are involved in MDR1response to benomyl. Interestingly, the BRE contains a perfect match to the binding consensus of Mcm1p, raising the possibility that MDR1may be a direct target of this transcriptional activator. In conclusion, while the identity of the trans-acting factors that bind to the BRE and HRE remains to be confirmed, the tools we have developed during characterization of the cis-acting elements of the MDR1promoter should now serve to elucidate the nature of the components that modulate its activity. RESUME La surexpression du gène MDR1 (pour Résistance Multidrogue 1), qui code pour un transporteur de la famille des Major Facilitators, est l'un des mécanismes observés dans les isolats cliniques de la levure Candida albicans développant une résistance aux agents antifongiques appelés azoles. Pour mieux comprendre ce phénomène, les éléments de régulation agissant en cis dans un système rapporteur modulable sous le contrôle du promoteur MDR1 ont été caractérisés. Dans une souche sensible aux azoles, la transcription de ce rapporteur est transitoirement surélevée en réponse soit au bénomyl soit à l'agent oxydant H2O2, alors que son expression est constitutivement élevée dans une souche résistante aux azoles (souche FR2). Deux éléments de régulation agissant en cis, nécessaires et suffisants pour transmettre les mêmes réponses transcriptionnelles à un promoteur hétérologue (CDR2), ont été identifiés dans le promoteur MDR1. Le premier élément, appelé BRE (pour Elément de Réponse au Bénomyl, de -296 à -260 par rapport au codon d'initiation ATG) est requis pour la surexpression de MDR1dépendante du bénomyl et pour l'expression constitutive de MDR1 dans FR2. Le deuxième élément, appelé HRE (pour Elément de Réponse à l'H2O2, de -561 à -520), est requis pour la surexpression de MDR1 dépendante de l'H2O2, mais n'est pas impliqué dans l'expression constitutive du gène MDR1. Deux sites de fixation potentiels (TTAG/CTAA) pour le facteur de transcription Cap1p ont été identifiés dans l'élément HRE. De plus, l'inactivation de CAP1 abolit la réponse transitoire à l'H2O2 et diminua significativement la réponse transitoire au bénomyl. Cap1p, qui est impliqué dans les réponses de la cellule au stress oxydatif, doit donc jouer un rôle positif en trans dans la surexpression de MDR1 dépendante du bénomyl et de l'H2O2. Cependant, ce n'est pas le seul facteur de transcription impliqué dans la réponse au bénomyl. Un élément BRE d'une longueur minimale (de -290 à -273) a également été défini et est suffisant pour détecter une interaction spécifique in vitro avec des protéines provenant d'extraits bruts de C. albicans. L'analyse du profil de transcription d'une paire d'isolats cliniques comprenant une souche résistante aux azoles surexprimant MDR1, et d'une souche sensible aux azoles exposée au bénomyl, a révélé que les gènes spécifiquement surexprimés par le bénomyl contiennent dans leurs promoteurs un ou plusieurs sites de fixation pour Cap1p. Ceci renforce l'idée que Cap1p joue un rôle dans la surexpression de MDR1dépendante du bénomyl. Une ou deux séquences ressemblant à l'élément BRE ont également été identifiées dans la plupart des gènes corégulés avec MDR1 dans ces deux conditions, ce qui était attendu compte-tenu du rôle joué par cet élément dans les deux processus. Une collection de 147 mutants dans lesquels un seul facteur de transcription est inactivé a été testée pour la perte de réponse au bénomyl de MDR1. Malheureusement, la surexpression de MDR1 dépendante du bénomyl n'a été perdue dans aucun des mutants testés. Néanmoins, une diminution significative de la réponse a été observée chez des mutants dans lesquels le facteur de transcription à MADS-box Mcm1p et le facteur de transcription à doigts de zinc de type C2H2 orf19.13374p ont été inactivés, suggérant que Mcm1p et orf19.13374p sont impliqués dans la réponse de MDR1au bénomyl. Il est intéressant de noter que la BRE contient une séquence qui s'aligne parfaitement avec la séquence consensus du site de fixation de Mcm1p, ce qui soulève la possibilité que MDR1 pourrait être une cible directe de cet activateur transcriptionnel. En conclusion, alors que l'identité des facteurs agissant en trans en se fixant à la BRE et à la HRE reste à être confirmée, les outils que nous avons développés au cours de la caractérisation des éléments agissant en cis sur le promoteur MDR1 peut maintenant servir à élucider la nature des composants modulant son activité.

GENCODE: the reference human genome annotation for The ENCODE Project.

The GENCODE Consortium aims to identify all gene features in the human genome using a combination of computational analysis, manual annotation, and experimental validation. Since the first public release of this annotation data set, few new protein-coding loci have been added, yet the number of alternative splicing transcripts annotated has steadily increased. The GENCODE 7 release contains 20,687 protein-coding and 9640 long noncoding RNA loci and has 33,977 coding transcripts not represented in UCSC genes and RefSeq. It also has the most comprehensive annotation of long noncoding RNA (lncRNA) loci publicly available with the predominant transcript form consisting of two exons. We have examined the completeness of the transcript annotation and found that 35% of transcriptional start sites are supported by CAGE clusters and 62% of protein-coding genes have annotated polyA sites. Over one-third of GENCODE protein-coding genes are supported by peptide hits derived from mass spectrometry spectra submitted to Peptide Atlas. New models derived from the Illumina Body Map 2.0 RNA-seq data identify 3689 new loci not currently in GENCODE, of which 3127 consist of two exon models indicating that they are possibly unannotated long noncoding loci. GENCODE 7 is publicly available from gencodegenes.org and via the Ensembl and UCSC Genome Browsers.

HER2 shedding and serum HER2 extracellular domain: Biology and clinical utility in breast cancer.

The transmembrane protein HER2 is over-expressed in approximately 15% of invasive breast cancers as a result of HER2 gene amplification. HER2 proteolytic cleavage (HER2 shedding) generates soluble truncated HER2 molecules that include only the extracellular domain and the concentration of which can be measured in the serum fraction of blood. HER2 shedding also generates a constitutively active truncated intracellular receptor of 95kDa (p95(HER2)). Another soluble truncated HER2 protein (Herstatin), which can also be found in serum, is the product of an alternatively spliced HER2 transcript. Recent preclinical findings may provide crucial insights into the biological and clinical relevance of increased sHER2 concentrations for the outcome of HER2-positive breast cancer and sensitivity to trastuzumab and lapatinib treatment. We present here the most recent findings about the role and biology of sHER2 based on data obtained using a standardized test, which has been cleared by FDA in 2000, for measuring sHER2. This test includes quality control assessments and has been already widely used to evaluate the clinical utility of sHER2 as a biomarker in breast cancer. We will describe in detail data concerning the assessment of sHER2 as a surrogate maker to optimize the evaluation of the HER2 status of a primary tumor and as a prognosis and predictive marker of response to therapies, both in early and metastatic breast cancer.