Comprehensive Cohort Analysis of Mutational Spectrum in Early Onset Breast Cancer Patients.

Early onset breast cancer (EOBC), diagnosed at age ~40 or younger, is associated with a poorer prognosis and higher mortality rate compared to breast cancer diagnosed at age 50 or older. EOBC poses a serious threat to public health and requires in-depth investigation. We studied a cohort comprising 90 Taiwanese female patients, aiming to unravel the underlying mechanisms of EOBC etiopathogenesis. Sequence data generated by whole-exome sequencing (WES) and whole-genome sequencing (WGS) from white blood cell (WBC)-tumor pairs were analyzed to identify somatic missense mutations, copy number variations (CNVs) and germline missense mutations. Similar to regular breast cancer, the key somatic mutation-susceptibility genes of EOBC include TP53 (40% prevalence), PIK3CA (37%), GATA3 (17%) and KMT2C (17%), which are frequently reported in breast cancer; however, the structural protein-coding genes MUC17 (19%), FLG (16%) and NEBL (11%) show a significantly higher prevalence in EOBC. Furthermore, the top 2 genes harboring EOBC germline mutations, MUC16 (19%) and KRT18 (19%), encode structural proteins. Compared to conventional breast cancer, an unexpectedly higher number of EOBC susceptibility genes encode structural proteins. We suspect that mutations in structural proteins may increase physical permeability to environmental hormones and carcinogens and cause breast cancer to occur at a young age.

Application of cell-free DNA sequencing in characterization of bloodborne microbes and the study of microbe-disease interactions.

It is an important issue whether microorganisms can live harmoniously with normal cells in the cardiovascular system. The answer to the question will have enormous impact on medical microbiology. To address the issue, it is essential to identify and characterize the bloodborne microbes in an efficient and comprehensive manner. Due to microbial sequence complexity and the composition of significant number of unknown microbial species in the circulatory system, traditional approaches using cell culture, PCR, or microarray are not suitable for the purpose. Recent reports indicate that cell-free DNA (cfDNA) sequencing using next-generation sequencing (NGS) or single-molecule sequencing (SMS), together with bioinformatics approaches, possesses a strong potential enabling us to distinguish microbial species at the nucleotide level. Multiple studies using microbial cfDNA sequencing to identify microbes for septic patients have shown strong agreement with cell culture. Similar approaches have also been applied to reveal previously unidentified microorganisms or to demonstrate the feasibility of comprehensive assessment of bloodborne microorganisms for healthy and/or diseased individuals. SMS using either SMRT (single-molecule real-time) sequencing or Nanopore sequencing are providing new momentum to reinforce this line of investigation. Taken together, microbial cfDNA sequencing provides a novel opportunity allowing us to further understand the involvement of bloodborne microbes in development of diseases. Similar approaches should also be applicable to the study of metagenomics for sufficient and comprehensive analysis of microbial species living in various environments. This article reviews this line of research and discuss the methodological approaches that have been developed, or are likely to be developed in the future, which may have strong potential to facilitate cfDNA- and cfRNA-based studies of cancer and acute/chronic diseases, in the hope that a better understanding of the hidden microbes in the circulatory system will improve diagnosis, prevention and treatment of problematic diseases.

Genomic sequencing of Troides aeacus nucleopolyhedrovirus (TraeNPV) from golden birdwing larvae (Troides aeacus formosanus) to reveal defective Autographa californica NPV genomic features.

BACKGROUND: The golden birdwing butterfly (Troides aeacus formosanus) is a rarely observed species in Taiwan. Recently, a typical symptom of nuclear polyhedrosis was found in reared T. aeacus larvae. From the previous Kimura-2 parameter (K-2-P) analysis based on the nucleotide sequence of three genes in this isolate, polh, lef-8 and lef-9, the underlying virus did not belong to any known nucleopolyhedrovirus (NPV) species. Therefore, this NPV was provisionally named "TraeNPV". To understand this NPV, the nucleotide sequence of the whole TraeNPV genome was determined using next-generation sequencing (NGS) technology. RESULTS: The genome of TraeNPV is 125,477 bp in length with 144 putative open reading frames (ORFs) and its GC content is 40.45%. A phylogenetic analysis based on the 37 baculoviral core genes suggested that TraeNPV is a Group I NPV that is closely related to Autographa californica nucleopolyhedrovirus (AcMNPV). A genome-wide analysis showed that TraeNPV has some different features in its genome compared with other NPVs. Two novel ORFs (Ta75 and Ta139), three truncated ORFs (pcna, he65 and bro) and one duplicated ORF (38.7 K) were found in the TraeNPV genome; moreover, there are fewer homologous regions (hrs) than there are in AcMNPV, which shares eight hrs within the TraeNPV genome. TraeNPV shares similar genomic features with AcMNPV, including the gene content, gene arrangement and gene/genome identity, but TraeNPV lacks 15 homologous ORFs from AcMNPV in its genome, such as ctx, host cell-specific factor 1 (hcf-1), PNK/PNL, vp15, and apsup, which are involved in the auxiliary functions of alphabaculoviruses. CONCLUSIONS: Based on these data, TraeNPV would be clarified as a new NPV species with defective AcMNPV genomic features. The precise relationship between TraeNPV and other closely related NPV species were further investigated. This report could provide comprehensive information on TraeNPV for evolutionary insights into butterfly-infected NPV.

Single cell transcriptome analysis of MCF-7 reveals consistently and inconsistently expressed gene groups each associated with distinct cellular localization and functions.

Single cell transcriptome (SCT) analysis provides superior resolution to illustrate tumor cell heterogeneity for clinical implications. We characterized four SCTs of MCF-7 using 143 housekeeping genes (HKGs) as control, of which lactate dehydrogenase B (LDHB) expression is silenced. These SCT libraries mapped to 11,423, 11,486, 10,380, and 11,306 RefSeq genes (UCSC), respectively. High consistency in HKG expression levels across all four SCTs, along with transcriptional silencing of LDHB, was observed, suggesting a high sensitivity and reproducibility of the SCT analysis. Cross-library comparison on expression levels by scatter plotting revealed a linear correlation and an 83-94% overlap in transcript isoforms and expressed genes were also observed. To gain insight of transcriptional diversity among the SCTs, expressed genes were split into consistently expressed (CE) (expressed in all SCTs) and inconsistently expressed (IE) (expressed in some but not all SCTs) genes for further characterization, along with the 142 expressed HKGs as a reference. Distinct transcriptional strengths were found among these groups, with averages of 1,612.0, 88.0 and 1.2 FPKM for HKGs, CE and IE, respectively. Comparison between CE and IE groups further indicated that expressions of CE genes vary more significantly than that of IE genes. Gene Ontology analysis indicated that proteins encoded by CE genes are mainly involved in fundamental intracellular activities, while proteins encoded by IE genes are mainly for extracellular activities, especially acting as receptors or ion channels. The diversified gene expressions, especially for those encoded by IE genes, may contribute to cancer drug resistance.

Analysis of microbial sequences in plasma cell-free DNA for early-onset breast cancer patients and healthy females.

BACKGROUND: Cell-free circulating DNA (cfDNA) is becoming a useful biopsy for noninvasive diagnosis of diseases. Microbial sequences in plasma cfDNA may provide important information to improve prognosis and treatment. We have developed a stringent method to identify microbial species via microbial cfDNA in the blood plasma of early-onset breast cancer (EOBC) patients and healthy females. Empirically, microbe-originated sequence reads were identified by mapping non-human PE reads in cfDNA libraries to microbial databases. Those mapped concordantly to unique microbial species were assembled into contigs, which were subsequently aligned to the same databases. Microbial species uniquely aligned were identified and compared across all individuals on MCRPM (Microbial CfDNA Reads Per Million quality PE reads) basis. RESULTS: The predominant microbial cfDNAs in all plasma samples examined are originated from bacteria and these bacteria were limited to only a few genera. Among those, Acinetobacter johnsonii XBB1 and low levels of Mycobacterium spp. were commonly found in all healthy females, but also present in an EOBC patient. Compared to those in healthy counterparts, bacterial species in EOBC patients are more diverse and more likely to present at high levels. Among these three EOBC patients tested, a patient who has record high titer (2,724 MCRPM) of Pseudomonas mendocina together with 8.82 MCRPM of Pannonibacter phragmitetus has passed away; another patient infected by multiple Sphingomonas species remains alive; while the third patient who has similar microbial species (Acinetobacter johnsonii XBB1) commonly seen in normal controls is having a normal life. CONCLUSIONS: Our preliminary data on the profiles of microbial cfDNA sequences suggested that it may have some prognostic value in cancer patients. Validation in larger number of patients is warranted.

Antibody variable domain interface and framework sequence requirements for stability and function by high-throughput experiments.

Protein structural stability and biological functionality are dictated by the formation of intradomain cores and interdomain interfaces, but the intricate sequence-structure-function interrelationships in the packing of protein cores and interfaces remain difficult to elucidate due to the intractability of enumerating all packing possibilities and assessing the consequences of all the variations. In this work, groups of ß strand residues of model antibody variable domains were randomized with saturated mutagenesis and the functional variants were selected for high-throughput sequencing and high-throughput thermal stability measurements. The results show that the sequence preferences of the intradomain hydrophobic core residues are strikingly flexible among hydrophobic residues, implying that these residues are coupled indirectly with antigen binding through energetic stabilization of the protein structures. By contrast, the interdomain interface residues are directly coupled with antigen binding. The interdomain interface should be treated as an integral part of the antigen-binding site.

Loop-sequence features and stability determinants in antibody variable domains by high-throughput experiments.

Protein loops are frequently considered as critical determinants in protein structure and function. Recent advances in high-throughput methods for DNA sequencing and thermal stability measurement have enabled effective exploration of sequence-structure-function relationships in local protein regions. Using these data-intensive technologies, we investigated the sequence-structure-function relationships of six complementarity-determining regions (CDRs) and ten non-CDR loops in the variable domains of a model vascular endothelial growth factor (VEGF)-binding single-chain antibody variable fragment (scFv) whose sequence had been optimized via a consensus-sequence approach. The results show that only a handful of residues involving long-range tertiary interactions distant from the antigen-binding site are strongly coupled with antigen binding. This implies that the loops are passive regions in protein folding; the essential sequences of these regions are dictated by conserved tertiary interactions and the consensus local loop-sequence features contribute little to protein stability and function.

Global assessment of Antrodia cinnamomea-induced microRNA alterations in hepatocarcinoma cells.

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.

Concordant and discordant regulation of target genes by miR-31 and its isoforms.

It has been shown that imprecise cleavage of a primary or precursor RNA by Drosha or Dicer, respectively, may yield a group of microRNA (miRNA) variants designated as "isomiR". Variations in the relative abundance of isoforms for a given miRNA among different species and different cell types beg the question whether these isomiRs might regulate target genes differentially. We compared the capacity of three miR-31 isoforms (miR-31-H, miR-31-P, and miR-31-M), which differ only slightly in their 5'- and/or 3'-end sequences, to regulate several known targets and a predicted target, Dicer. Notably, we found isomiR-31s displayed concordant and discordant regulation of 6 known target genes. Furthermore, we validated a predicted target gene, Dicer, to be a novel target of miR-31 but only miR-31-P could directly repress Dicer expression in both MCF-7 breast cancer cells and A549 lung cancer cells, resulting in their enhanced sensitivity to cisplatin, a known attribute of Dicer knockdown. This was further supported by reporter assay using full length 3'-untranslated region (UTR) of Dicer. Our findings not only revealed Dicer to be a direct target of miR-31, but also demonstrated that isomiRs displayed similar and disparate regulation of target genes in cell-based systems. Coupled with the variations in the distribution of isomiRs among different cells or conditions, our findings support the possibility of fine-tuning gene expression by miRNAs.

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.

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.

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.

Whole-genome cartography of estrogen receptor alpha binding sites.

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.

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.

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.

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.