Genome-wide mapping of RELA(p65) binding identifies E2F1 as a transcriptional activator recruited by NF-kappaB upon TLR4 activation

NF-kappaB is a key mediator of inflammation. Here, we mapped the genome-wide loci bound by the RELA subunit of NF-kappaB in lipopolysaccharide (LPS)-stimulated human monocytic cells, and together with global gene expression profiling, found an overrepresentation of the E2F1-binding motif among RELA-bound loci associated with NF-kappaB target genes. Knockdown of endogenous E2F1 impaired the LPS inducibility of the proinflammatory cytokines CCL3(MIP-1alpha), IL23A(p19), TNF-alpha, and IL1-beta. Upon LPS stimulation, E2F1 is rapidly recruited to the promoters of these genes along with p50/RELA heterodimer via a mechanism that is dependent on NF-kappaB activation. Together with the observation that E2F1 physically interacts with p50/RELA in LPS-stimulated cells, our findings suggest that NF-kappaB recruits E2F1 to fully activate the transcription of NF-kappaB target genes. Global gene expression profiling subsequently revealed a spectrum of NF-kappaB target genes that are positively regulated by E2F1, further demonstrating the critical role of E2F1 in the Toll-like receptor 4 pathway.

Mapping of transcription factor binding regions in mammalian cells by ChIP: comparison of array- and sequencing-based technologies

Recent progress in mapping transcription factor (TF) binding regions can largely be credited to chromatin immunoprecipitation (ChIP) technologies. We compared strategies for mapping TF binding regions in mammalian cells using two different ChIP schemes: ChIP with DNA microarray analysis (ChIP-chip) and ChIP with DNA sequencing (ChIP-PET). We first investigated parameters central to obtaining robust ChIP-chip data sets by analyzing STAT1 targets in the ENCODE regions of the human genome, and then compared ChIP-chip to ChIP-PET. We devised methods for scoring and comparing results among various tiling arrays and examined parameters such as DNA microarray format, oligonucleotide length, hybridization conditions, and the use of competitor Cot-1 DNA. The best performance was achieved with high-density oligonucleotide arrays, oligonucleotides >/=50 bases (b), the presence of competitor Cot-1 DNA and hybridizations conducted in microfluidics stations. When target identification was evaluated as a function of array number, 80%-86% of targets were identified with three or more arrays. Comparison of ChIP-chip with ChIP-PET revealed strong agreement for the highest ranked targets with less overlap for the low ranked targets. With advantages and disadvantages unique to each approach, we found that ChIP-chip and ChIP-PET are frequently complementary in their relative abilities to detect STAT1 targets for the lower ranked targets; each method detected validated targets that were missed by the other method. The most comprehensive list of STAT1 binding regions is obtained by merging results from ChIP-chip and ChIP-sequencing. Overall, this study provides information for robust identification, scoring, and validation of TF targets using ChIP-based technologies.

Fusion transcripts and transcribed retrotransposed loci discovered through comprehensive transcriptome analysis using Paired-End diTags (PETs)

Identification of unconventional functional features such as fusion transcripts is a challenging task in the effort to annotate all functional DNA elements in the human genome. Paired-End diTag (PET) analysis possesses a unique capability to accurately and efficiently characterize the two ends of DNA fragments, which may have either normal or unusual compositions. This unique nature of PET analysis makes it an ideal tool for uncovering unconventional features residing in the human genome. Using the PET approach for comprehensive transcriptome analysis, we were able to identify fusion transcripts derived from genome rearrangements and actively expressed retrotransposed pseudogenes, which would be difficult to capture by other means. Here, we demonstrate this unique capability through the analysis of 865,000 individual transcripts in two types of cancer cells. In addition to the characterization of a large number of differentially expressed alternative 5' and 3' transcript variants and novel transcriptional units, we identified 70 fusion transcript candidates in this study. One was validated as the product of a fusion gene between BCAS4 and BCAS3 resulting from an amplification followed by a translocation event between the two loci, chr20q13 and chr17q23. Through an examination of PETs that mapped to multiple genomic locations, we identified 4055 retrotransposed loci in the human genome, of which at least three were found to be transcriptionally active. The PET mapping strategy presented here promises to be a useful tool in annotating the human genome, especially aberrations in human cancer genomes.

Pseudogenes in the ENCODE regions: consensus annotation, analysis of transcription, and evolution

Arising from either retrotransposition or genomic duplication of functional genes, pseudogenes are "genomic fossils" valuable for exploring the dynamics and evolution of genes and genomes. Pseudogene identification is an important problem in computational genomics, and is also critical for obtaining an accurate picture of a genome's structure and function. However, no consensus computational scheme for defining and detecting pseudogenes has been developed thus far. As part of the ENCyclopedia Of DNA Elements (ENCODE) project, we have compared several distinct pseudogene annotation strategies and found that different approaches and parameters often resulted in rather distinct sets of pseudogenes. We subsequently developed a consensus approach for annotating pseudogenes (derived from protein coding genes) in the ENCODE regions, resulting in 201 pseudogenes, two-thirds of which originated from retrotransposition. A survey of orthologs for these pseudogenes in 28 vertebrate genomes showed that a significant fraction ( approximately 80%) of the processed pseudogenes are primate-specific sequences, highlighting the increasing retrotransposition activity in primates. Analysis of sequence conservation and variation also demonstrated that most pseudogenes evolve neutrally, and processed pseudogenes appear to have lost their coding potential immediately or soon after their emergence. In order to explore the functional implication of pseudogene prevalence, we have extensively examined the transcriptional activity of the ENCODE pseudogenes. We performed systematic series of pseudogene-specific RACE analyses. These, together with complementary evidence derived from tiling microarrays and high throughput sequencing, demonstrated that at least a fifth of the 201 pseudogenes are transcribed in one or more cell lines or tissues.

Measurement of psi2S radiative decays

Using 14 x 10(6) psi(2S) events accumulated at the BESII detector, we report first measurements of branching fractions or upper limits for psi(2S) decays into gammapp, gamma2(pi+pi-), gammaKS0K+pi-+c.c., gammaK+K-pi+pi-, gammaK*0K-pi++c.c., gammaK*0K*0, gammapi+pi-pp, gamma2(K+K-), gamma3(pi+pi-), and gamma2(pi+pi-)K+K- with the invariant mass of hadrons below 2.9 GeV/c2. We also report branching fractions of psi(2S) decays into 2(pi+pi-)pi0, omegapi+pi-, omegaf2(1270), b1+/-pi-/+, and pi02(pi+pi-)K+K-.

Pathway aberrations of murine melanoma cells observed in Paired-End diTag transcriptomes

BACKGROUND

Melanoma is the major cause of skin cancer deaths and melanoma incidence doubles every 10 to 20 years. However, little is known about melanoma pathway aberrations. Here we applied the robust Gene Identification Signature Paired End diTag (GIS-PET) approach to investigate the melanoma transcriptome and characterize the global pathway aberrations.

METHODS

GIS-PET technology directly links 5' mRNA signatures with their corresponding 3' signatures to generate, and then concatenate, PETs for efficient sequencing. We annotated PETs to pathways of KEGG database and compared the murine B16F1 melanoma transcriptome with three non-melanoma murine transcriptomes (Melan-a2 melanocytes, E14 embryonic stem cells, and E17.5 embryo). Gene expression levels as represented by PET counts were compared across melanoma and melanocyte libraries to identify the most significantly altered pathways and investigate the expression levels of crucial cancer genes.

RESULTS

Melanin biosynthesis genes were solely expressed in the cells of melanocytic origin, indicating the feasibility of using the PET approach for transcriptome comparison. The most significantly altered pathways were metabolic pathways, including upregulated pathways: purine metabolism, aminophosphonate metabolism, tyrosine metabolism, selenoamino acid metabolism, galactose utilization, nitrobenzene degradation, and bisphenol A degradation; and downregulated pathways: oxidative phosphorylation, ATPase synthesis, TCA cycle, pyruvate metabolism, and glutathione metabolism. The downregulated pathways concurrently indicated a slowdown of mitochondrial activities. Mitochondrial permeability was also significantly altered, as indicated by transcriptional activation of ATP/ADP, citrate/malate, Mg++, fatty acid and amino acid transporters, and transcriptional repression of zinc and metal ion transporters. Upregulation of cell cycle progression, MAPK, and PI3K/Akt pathways were more limited to certain region(s) of the pathway. Expression levels of c-Myc and Trp53 were also higher in melanoma. Moreover, transcriptional variants resulted from alternative transcription start sites or alternative polyadenylation sites were found in Ras and genes encoding adhesion or cytoskeleton proteins such as integrin, beta-catenin, alpha-catenin, and actin.

CONCLUSION

The highly correlated results unmistakably point to a systematic downregulation of mitochondrial activities, which we hypothesize aims to downgrade the mitochondria-mediated apoptosis and the dependency of cancer cells on angiogenesis. Our results also demonstrate the advantage of using the PET approach in conjunction with KEGG database for systematic pathway analysis.

Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project

We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.

Whole-genome cartography of estrogen receptor alpha binding sites

Using a chromatin immunoprecipitation-paired end diTag cloning and sequencing strategy, we mapped estrogen receptor α (ERα) binding sites in MCF-7 breast cancer cells. We identified 1,234 high confidence binding clusters of which 94% are projected to be bona fide ERα binding regions. Only 5% of the mapped estrogen receptor binding sites are located within 5 kb upstream of the transcriptional start sites of adjacent genes, regions containing the proximal promoters, whereas vast majority of the sites are mapped to intronic or distal locations (>5 kb from 5′ and 3′ ends of adjacent transcript), suggesting transcriptional regulatory mechanisms over significant physical distances. Of all the identified sites, 71% harbored putative full estrogen response elements (EREs), 25% bore ERE half sites, and only 4% had no recognizable ERE sequences. Genes in the vicinity of ERα binding sites were enriched for regulation by estradiol in MCF-7 cells, and their expression profiles in patient samples segregate ERα-positive from ERα-negative breast tumors. The expression dynamics of the genes adjacent to ERα binding sites suggest a direct induction of gene expression through binding to ERE-like sequences, whereas transcriptional repression by ERα appears to be through indirect mechanisms. Our analysis also indicates a number of candidate transcription factor binding sites adjacent to occupied EREs at frequencies much greater than by chance, including the previously reported FOXA1 sites, and demonstrate the potential involvement of one such putative adjacent factor, Sp1, in the global regulation of ERα target genes. Unexpectedly, we found that only 22%–24% of the bona fide human ERα binding sites were overlapping conserved regions in whole genome vertebrate alignments, which suggest limited conservation of functional binding sites. Taken together, this genome-scale analysis suggests complex but definable rules governing ERα binding and gene regulation.

Estrogen receptors (ERs) play key roles in facilitating the transcriptional effects of hormone functions in target tissues. To obtain a genome-wide view of ERα binding sites, we applied chromatin immunoprecipitation coupled with a cloning and sequencing strategy using chromatin immunoprecipitation pair end-tagging technology to map ERα binding sites in MCF-7 human breast cancer cells. We identified 1,234 high quality ERα binding sites in the human genome and demonstrated that the binding sites are frequently adjacent to genes significantly associated with breast cancer disease status and outcome. The mapping results also revealed that ERα can influence gene expression across distances of up to 100 kilobases or more, that genes that are induced or repressed utilize sites in different regions relative to the transcript (suggesting different mechanisms of action), and that ERα binding sites are only modestly conserved in evolution. Using computational approaches, we identified potential interactions with other transcription factor binding sites adjacent to the ERα binding elements. Taken together, these findings suggest complex but definable rules governing ERα binding and gene regulation and provide a valuable dataset for mapping the precise control nodes for one of the most important nuclear hormone receptors in breast cancer biology.

Measurements of the continuum R(uds) and R values in e(+)e(-) annihilation in the energy region between 3.650 and 3.872 GeV

We report measurements of the continuum R(uds) near the center-of-mass energy of 3.70 GeV, the R[uds(c)+psi(3770)](s) and the R(had)(s) values in e(+)e(-) annihilation at 68 energy points in the energy region between 3.650 and 3.872 GeV with the BES-II detector at the BEPC Collider. We obtain the R(uds) for the continuum light hadron (containing u, d, and s quarks) production near the DD threshold to be R(uds)=2.141+/-0.025+/-0.085.

Search for invisible decays of eta and eta' in J/psi --> phi eta and phi eta'

Using a data sample of 58 x 10(6) J/psi decays collected with the Beijing Spectrometer II detector at the Beijing Electron Positron Collider, searches for invisible decays of eta and eta' in J/psi to phi eta and phi eta' are performed. The phi signals, which are reconstructed in K+K- final states, are used to tag the eta and eta' decays. No signals are found for the invisible decays of either eta or eta', and upper limits at the 90% confidence level are determined to be 1.65 x 10(-3) for the ratio B(eta-->invisible)/B(eta --> gamma gamma) and 6.69 x 10(-2) for B(eta' --> invisible)/B(eta' --> gammagamma). These are the first searches for eta and eta' decays into invisible final states.

Global mapping of c-Myc binding sites and target gene networks in human B cells

The protooncogene MYC encodes the c-Myc transcription factor that regulates cell growth, cell proliferation, cell cycle, and apoptosis. Although deregulation of MYC contributes to tumorigenesis, it is still unclear what direct Myc-induced transcriptomes promote cell transformation. Here we provide a snapshot of genome-wide, unbiased characterization of direct Myc binding targets in a model of human B lymphoid tumor using ChIP coupled with pair-end ditag sequencing analysis (ChIP-PET). Myc potentially occupies > 4,000 genomic loci with the majority near proximal promoter regions associated frequently with CpG islands. Using gene expression profiles with ChIP-PET, we identified 668 direct Myc-regulated gene targets, including 48 transcription factors, indicating that Myc is a central transcriptional hub in growth and proliferation control. This first global genomic view of Myc binding sites yields insights of transcriptional circuitries and cis regulatory modules involving Myc and provides a substantial framework for our understanding of mechanisms of Myc-induced tumorigenesis.

Observation of a broad 1-- resonant structure around 1.5 GeV/c2 in the K+K- mass spectrum in J/psi-->K+K-pi0

A broad peak is observed at low K+K- invariant mass in J/psi-->K+K-pi(0) decays found in a sample of 5.8x10(7) J/psi events collected with the BESII detector. The statistical significance of the broad resonance is much larger than 5sigma. A partial wave analysis shows that the J;{PC} of this structure is 1--. Its pole position is determined to be [1576(-55)(+49)(stat)-91+98(syst)] MeV/c(2)-i/2[818(-23)(+22)(stat)-133+64(syst)] MeV/c(2). These parameters are not compatible with any known meson resonances.

Measurements of the branching fractions for psi(3770)-->D(0)D[over ](0), D+D-, DD[over ], and the resonance parameters of psi(3770) and psi(2S)

We measure the branching fractions for psi(3770)-->D(0)D[over ](0), D+D-, DD[over ], and non-DD[over ] to be (46.7+/-4.7+/-2.3)%, (36.9+/-3.7+/-2.8)%, (83.6+/-7.3+/-4.2)%, and (16.4+/-7.3+/-4.2)%, respectively. The resonance parameters of psi(3770) and psi(2S) are measured to be M_(psi(3770))=3772.2+/-0.7+/-0.3 MeV, Gamma_(psi(3770))(tot)=26.9+/-2.4+/-0.3 MeV, and Gamma_(psi(3770))(ee)=251+/-26+/-11 eV; M_(psi(2S))=3685.5+/-0.0+/-0.3 MeV, Gamma_(psi(2S))(tot)=331+/-58+/-2 keV, and Gamma_(psi(2S))(ee)=2.330+/-0.036+/-0.110 keV. We also measure the light hadron R value to be R(uds)=2.262+/-0.054+/-0.109 in the energy region from 3.660 to 3.872 GeV.

PET-Tool: a software suite for comprehensive processing and managing of Paired-End diTag (PET) sequence data

BACKGROUND

We recently developed the Paired End diTag (PET) strategy for efficient characterization of mammalian transcriptomes and genomes. The paired end nature of short PET sequences derived from long DNA fragments raised a new set of bioinformatics challenges, including how to extract PETs from raw sequence reads, and correctly yet efficiently map PETs to reference genome sequences. To accommodate and streamline data analysis of the large volume PET sequences generated from each PET experiment, an automated PET data process pipeline is desirable.

RESULTS

We designed an integrated computation program package, PET-Tool, to automatically process PET sequences and map them to the genome sequences. The Tool was implemented as a web-based application composed of four modules: the Extractor module for PET extraction; the Examiner module for analytic evaluation of PET sequence quality; the Mapper module for locating PET sequences in the genome sequences; and the Project Manager module for data organization. The performance of PET-Tool was evaluated through the analyses of 2.7 million PET sequences. It was demonstrated that PET-Tool is accurate and efficient in extracting PET sequences and removing artifacts from large volume dataset. Using optimized mapping criteria, over 70% of quality PET sequences were mapped specifically to the genome sequences. With a 2.4 GHz LINUX machine, it takes approximately six hours to process one million PETs from extraction to mapping.

CONCLUSION

The speed, accuracy, and comprehensiveness have proved that PET-Tool is an important and useful component in PET experiments, and can be extended to accommodate other related analyses of paired-end sequences. The Tool also provides user-friendly functions for data quality check and system for multi-layer data management.

Observation of two new N* peaks in J/psi-->ppi-n and ppi+n decays

The decay J/psi-->NNpi provides an effective isospin 1/2 filter for the piN system due to isospin conservation. Using 58x10(6) J/psi decays collected with the Beijing Electromagnetic Spectrometer at the Beijing Electron Positron Collider, more than 100 thousand J/psi-->ppi-n+c.c. events are obtained. Besides the two well-known N* peaks at around 1500 MeV/c2 and 1670 MeV/c2, there are two new, clear N* peaks in the ppi invariant mass spectrum around 1360 MeV/c2 and 2030 MeV/c2 with statistical significance of 11sigma and 13sigma, respectively. We identify these as the first direct observation of the N*(1440) peak and a long-sought missing N* peak above 2 GeV/c2 in the piN invariant mass spectrum.

Multiplex sequencing of paired-end ditags (MS-PET): a strategy for the ultra-high-throughput analysis of transcriptomes and genomes

The paired-end ditagging (PET) technique has been shown to be efficient and accurate for large-scale transcriptome and genome analysis. However, as with other DNA tag-based sequencing strategies, it is constrained by the current efficiency of Sanger technology. A recently developed multiplex sequencing method (454-sequencing™) using picolitre-scale reactions has achieved a remarkable advance in efficiency, but suffers from short-read lengths, and a lack of paired-end information. To further enhance the efficiency of PET analysis and at the same time overcome the drawbacks of the new sequencing method, we coupled multiplex sequencing with paired-end ditagging (MS-PET) using modified PET procedures to simultaneously sequence 200 000 to 300 000 dimerized PET (diPET) templates, with an output of nearly half-a-million PET sequences in a single 4 h machine run. We demonstrate the utility and robustness of MS-PET by analyzing the transcriptome of human breast carcinoma cells, and by mapping p53 binding sites in the genome of human colorectal carcinoma cells. This combined sequencing strategy achieved an approximate 100-fold efficiency increase over the current standard for PET analysis, and furthermore enables the short-read-length multiplex sequencing procedure to acquire paired-end information from large DNA fragments.

Observation of a near-threshold enhancement in the omega(phi) mass spectrum from the doubly OZI-suppressed decay J/psi-->gamma(omega)phi

An enhancement near threshold is observed in the omega(phi) invariant mass spectrum from the doubly Okubo-Zweig-Iizuka-suppressed decays of J/psi-->gamma(omega)phi, based on a sample of 5.8 x 10(7) J/psi events collected with the BESII detector. A partial wave analysis shows that this enhancement favors JP=0+, and its mass and width are M=1812(+19)(-26)(stat)+/-18(syst) MeV/c2 and Gamma=105+/-20(stat)+/-28(syst) MeV/c2. The product branching fraction is determined to be B(J/psi-->gammaX)B(X-->omega(phi))=[2.61+/-0.27(stat)+/-0.65(syst)]x10(-4).

The Oct4 and Nanog transcription network regulates pluripotency in mouse embryonic stem cells

Oct4 and Nanog are transcription factors required to maintain the pluripotency and self-renewal of embryonic stem (ES) cells. Using the chromatin immunoprecipitation paired-end ditags method, we mapped the binding sites of these factors in the mouse ES cell genome. We identified 1,083 and 3,006 high-confidence binding sites for Oct4 and Nanog, respectively. Comparative location analyses indicated that Oct4 and Nanog overlap substantially in their targets, and they are bound to genes in different configurations. Using de novo motif discovery algorithms, we defined the cis-acting elements mediating their respective binding to genomic sites. By integrating RNA interference-mediated depletion of Oct4 and Nanog with microarray expression profiling, we demonstrated that these factors can activate or suppress transcription. We further showed that common core downstream targets are important to keep ES cells from differentiating. The emerging picture is one in which Oct4 and Nanog control a cascade of pathways that are intricately connected to govern pluripotency, self-renewal, genome surveillance and cell fate determination.

A global map of p53 transcription-factor binding sites in the human genome

The ability to derive a whole-genome map of transcription-factor binding sites (TFBS) is crucial for elucidating gene regulatory networks. Herein, we describe a robust approach that couples chromatin immunoprecipitation (ChIP) with the paired-end ditag (PET) sequencing strategy for unbiased and precise global localization of TFBS. We have applied this strategy to map p53 targets in the human genome. From a saturated sampling of over half a million PET sequences, we characterized 65,572 unique p53 ChIP DNA fragments and established overlapping PET clusters as a readout to define p53 binding loci with remarkable specificity. Based on this information, we refined the consensus p53 binding motif, identified at least 542 binding loci with high confidence, discovered 98 previously unidentified p53 target genes that were implicated in novel aspects of p53 functions, and showed their clinical relevance to p53-dependent tumorigenesis in primary cancer samples.

Observation of a resonance in Chi(1835) in J/psi --> gammapi+ pi- eta-

The decay channel J/psi --> gamma(pi)(+)pi(-)eta is analyzed using a sample of 5.8 x 10(7) J/psi events collected with the BESII detector. A resonance, the Chi(1835), is observed in the pi(+)pi(-)eta invariant-mass spectrum with a statistical significance of 7.7 sigma. A fit with a Breit-Wigner function yields a mass M = 1833.7 +/- 6.1(stat) +/- 2.7(syst) MeV/c(2), a width Tau = 67.7 +/- 20.3(stat) +/- 7.7(syst) MeV/c(2), and a product branching fraction B(J/psi --> gammaChi) . B(Chi --> pi(+)pi(-)eta) = [2.2 +/- 0.4(stat) +/- 0.4(syst)] x 10(-4). The mass and width of the Chi(1835) are not compatible with any known meson resonance. Its properties are consistent with expectations for the state that produces the strong pp mass threshold enhancement observed in the J/psi --> gammapp process at BESII.

The transcriptional landscape of the mammalian genome

This study describes comprehensive polling of transcription start and termination sites and analysis of previously unidentified full-length complementary DNAs derived from the mouse genome. We identify the 5' and 3' boundaries of 181,047 transcripts with extensive variation in transcripts arising from alternative promoter usage, splicing, and polyadenylation. There are 16,247 new mouse protein-coding transcripts, including 5154 encoding previously unidentified proteins. Genomic mapping of the transcriptome reveals transcriptional forests, with overlapping transcription on both strands, separated by deserts in which few transcripts are observed. The data provide a comprehensive platform for the comparative analysis of mammalian transcriptional regulation in differentiation and development.

Gene identification signature (GIS) analysis for transcriptome characterization and genome annotation

We have developed a DNA tag sequencing and mapping strategy called gene identification signature (GIS) analysis, in which 5' and 3' signatures of full-length cDNAs are accurately extracted into paired-end ditags (PETs) that are concatenated for efficient sequencing and mapped to genome sequences to demarcate the transcription boundaries of every gene. GIS analysis is potentially 30-fold more efficient than standard cDNA sequencing approaches for transcriptome characterization. We demonstrated this approach with 116,252 PET sequences derived from mouse embryonic stem cells. Initial analysis of this dataset identified hundreds of previously uncharacterized transcripts, including alternative transcripts of known genes. We also uncovered several intergenically spliced and unusual fusion transcripts, one of which was confirmed as a trans-splicing event and was differentially expressed. The concept of paired-end ditagging described here for transcriptome analysis can also be applied to whole-genome analysis of cis-regulatory and other DNA elements and represents an important technological advance for genome annotation.

Observation of a threshold enhancement in the plambda invariant-mass spectrum

An enhancement near the m(p)+M(Lambda) mass threshold is observed in the combined pLambda and pLambda invariant-mass spectrum from J/psi-->pK(-)Lambda;+c.c. decays. It can be fit with an S-wave Breit-Wigner resonance with a mass m=2075+/-12(stat)+/-5(syst) MeV and a width of Gamma=90+/-35(stat)+/-9(syst) MeV; it can also be fit with a P-wave Breit-Wigner resonance. Evidence for a similar enhancement is also observed in psi(')-->pK(-)Lambda;+c.c. decays. The analysis is based on samples of 5.8x10(7)J/psi and 1.4x10(7)psi(') decays accumulated in the BES II detector at the Beijing Electron-Positron Collider.

5' Long serial analysis of gene expression (LongSAGE) and 3' LongSAGE for transcriptome characterization and genome annotation

Complete genome annotation relies on precise identification of transcription units bounded by a transcription initiation site (TIS) and a polyadenylation site (PAS). To facilitate this process, we developed a set of two complementary methods, 5' Long serial analysis of gene expression (LS) and 3'LS. These analyses are based on the original SAGE and LS methods coupled with full-length cDNA cloning, and enable the high-throughput extraction of the first and the last 20 bp of each transcript. We demonstrate that the mapping of 5'LS and 3'LS tags to the genome allows the localization of TIS and PAS. By using 537 tag pairs mapping to the region of known genes, we confirmed that >90% of the tag pairs appropriately assigned to the first and last exons. Moreover, by using tag sequences as primers for RT-PCRs, we were able to recover putative full-length transcripts in 81% of the attempts. This large-scale generation of transcript terminal tags is at least 20-40 times more efficient than full-length cDNA cloning and sequencing in the identification of complete transcription units. The apparent precision and deep coverage makes 5'LS and 3'LS an advanced approach for genome annotation through whole-transcriptome characterization.

Observation of the decay psi(2S)-->K0SK0L

The decay psi(2S)-->K(0)(S)K(0)(L) is observed using psi(2S) data collected with the Beijing Spectrometer at the Beijing Electron-Positron Collider; the branching fraction is determined to be B(psi(2S)-->K(0)(S)K(0)(L))=(5.24+/-0.47+/-0.48)x10(-5). Compared with J/psi-->K(0)(S)K(0)(L), the psi(2S) branching fraction is enhanced relative to the prediction of the perturbative QCD "12%" rule. The result, together with the branching fractions of psi(2S) decays to other pseudoscalar meson pairs (pi(+)pi(-) and K+K-), is used to investigate the relative phase between the three-gluon and the one-photon annihilation amplitudes of psi(2S) decays.