The pathophysiology of obesity and type 2 diabetes mellitus is associated with abnormalities in endocrine signaling in adipose tissue and one of the key signaling affectors operative in these disorders is the nuclear hormone transcription factor peroxisome proliferator-activated receptor-gamma (PPARgamma). PPARgamma has pleiotropic functions affecting a wide range of fundamental biological processes including the regulation of genes that modulate insulin sensitivity, adipocyte differentiation, inflammation and atherosclerosis. To date, only a limited number of direct targets for PPARgamma have been identified through research using the well established pre-adipogenic cell line, 3T3-L1. In order to obtain a genome-wide view of PPARgamma binding sites, we applied the pair end-tagging technology (ChIP-PET) to map PPARgamma binding sites in 3T3-L1 preadipocyte cells.

Humoral immunity against diverse pathogens is rapidly elicited from natural antibody repertoires of limited complexity. But the organizing principles underlying the antibody repertoires that facilitate this immunity are not well-understood. We used HER2 as a model immunogen and reverse-engineered murine antibody response through constructing an artificial antibody library encoded with rudimentary sequence and structural characteristics learned from high throughput sequencing of antibody variable domains. Antibodies selected in vitro from the phage-displayed synthetic antibody library bound to the model immunogen with high affinity and specificities, which reproduced the specificities of natural antibody responses. We conclude that natural antibody structural repertoires are shaped to allow functional antibodies to be encoded efficiently, within the complexity limit of an individual antibody repertoire, to bind to diverse protein antigens with high specificity and affinity. Phage-displayed synthetic antibody libraries, in conjunction with high-throughput sequencing, can thus be designed to replicate natural antibody responses and to generate novel antibodies against diverse antigens.

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.

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.

Recent studies have demonstrated a potent anticancer potential of medicinal fungus Antrodia cinnamomea, especially against hepatocarcinoma. These studies, however, were performed with prolonged treatments, and the early anticancer events remain missing. To probe the early anticancer mechanisms of A. cinnamomea, we treated SK-Hep-1 liver cancer cell with A. cinnamomea fruiting body extract for 2 and 4 hours, sequenced RNA samples with next-generation sequencing approach, and profiled the genome-wide miRNA and mRNA transcriptomes. Results unmistakably associated the early anticancer effect of A. cinnamomea fruiting body extract with a global downregulation of miRNAs which occurred solely in the A. cinnamomea fruiting body extract-treated SK-Hep-1 cells. Moreover, the inhibitory effect of A. cinnamomea fruiting body extract upon cancer miRNAs imposed no discrimination against any particular miRNA species, with oncomirs miR-21, miR-191 and major oncogenic clusters miR-17-92 and miR-106b-25 among the most severely downregulated. Western blotting further indicated a decrease in Drosha and Dicer proteins which play a key role in miRNA biogenesis, together with an increase of XRN2 known to participate in miRNA degradation pathway. Transcriptome profiling followed by GO and pathway analyses indicated that A. cinnamomea induced apoptosis, which was tightly associated with a downregulation of PI3K/AKT and MAPK pathways. Phosphorylation assay further suggested that JNK and c-Jun were closely involved in the apoptotic process. Taken together, our data indicated that the anticancer effect of A. cinnamomea can take place within a few hours by targeting multiple proteins and the miRNA system. A. cinnamomea indiscriminately induced a global downregulation of miRNAs by simultaneously inhibiting the key enzymes involved in miRNA maturation and activating XRN2 protein involved in miRNA degradation. Collapsing of the miRNA system together with downregulation of cell growth and survival pathways and activation of JNK signaling unleash the extrinsic and intrinsic apoptosis pathways, leading to the cancer cell death.

Epigenetic modifications are crucial for proper lineage specification and embryo development. To explore the chromatin modification landscapes in human ES cells, we profiled two histone modifications, H3K4me3 and H3K27me3, by ChIP coupled with the paired-end ditags sequencing strategy. H3K4me3 was found to be a prevalent mark and occurred in close proximity to the promoters of two-thirds of total human genes. Among the H3K27me3 loci identified, 56% are associated with promoters and the vast majority of them are comodified by H3K4me3. By deep-transcript digital counting, 80% of H3K4me3 and 36% of comodified promoters were found to be transcribed. Remarkably, we observed that different combinations of histone methylations are associated with genes from distinct functional categories. These global histone methylation maps provide an epigenetic framework that enables the discovery of novel transcriptional networks and delineation of different genetic compartments of the pluripotent cell genome.

Using a chromatin immunoprecipitation-paired end diTag cloning and sequencing strategy, we mapped estrogen receptor alpha (ERalpha) 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 ERalpha 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 ERalpha binding sites were enriched for regulation by estradiol in MCF-7 cells, and their expression profiles in patient samples segregate ERalpha-positive from ERalpha-negative breast tumors. The expression dynamics of the genes adjacent to ERalpha binding sites suggest a direct induction of gene expression through binding to ERE-like sequences, whereas transcriptional repression by ERalpha 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 ERalpha target genes. Unexpectedly, we found that only 22%-24% of the bona fide human ERalpha 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 ERalpha binding and gene regulation.

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.

Body fluid DNA sequencing is a powerful noninvasive approach for the diagnosis of genetic defects, infectious agents and diseases. The success relies on the quantity and quality of the DNA samples. However, numerous clinical samples are either at low quantity or of poor quality due to various reasons. To overcome these problems, we have developed T oligo-primed polymerase chain reaction (TOP-PCR) for full-length nonselective amplification of minute quantity of DNA fragments. TOP-PCR adopts homogeneous "half adaptor" (HA), generated by annealing P oligo (carrying a phosphate group at the 5' end) and T oligo (carrying a T-tail at the 3' end), for efficient ligation to target DNA and subsequent PCR amplification primed by the T oligo alone. Using DNA samples from body fluids, we demonstrate that TOP-PCR recovers minute DNA fragments and maintains the DNA size profile, while enhancing the major molecular populations. Our results also showed that TOP-PCR is a superior method for detecting apoptosis and outperforms the method adopted by Illumina for DNA amplification.

The cause of severe acute respiratory syndrome (SARS) has been identified as a new coronavirus. Whole genome sequence analysis of various isolates might provide an indication of potential strain differences of this new virus. Moreover, mutation analysis will help to develop effective vaccines.

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.

The Taiwanese (Formosan) macaque (Macaca cyclopis) is the only nonhuman primate endemic to Taiwan. This primate species is valuable for evolutionary studies and as subjects in medical research. However, only partial fragments of the mitochondrial genome (mitogenome) of this primate species have been sequenced, not mentioning its nuclear genome. We employed next-generation sequencing to generate 2 x 90 bp paired-end reads, followed by reference-assisted de novo assembly with multiple k-mer strategy to characterize the M. cyclopis mitogenome. We compared the assembled mitogenome with that of other macaque species for phylogenetic analysis. Our results show that, the M. cyclopis mitogenome consists of 16,563 nucleotides encoding for 13 protein-coding genes, 2 ribosomal RNAs and 22 transfer RNAs. Phylogenetic analysis indicates that M. cyclopis is most closely related to M. mulatta lasiota (Chinese rhesus macaque), supporting the notion of Asia-continental origin of M. cyclopis proposed in previous studies based on partial mitochondrial sequences. Our work presents a novel approach for assembling a mitogenome that utilizes the capabilities of de novo genome assembly with assistance of a reference genome. The availability of the complete Taiwanese macaque mitogenome will facilitate the study of primate evolution and the characterization of genetic variations for the potential usage of this species as a non-human primate model for medical research.