Screening the human exome: a comparison of whole genome and whole transcriptome sequencing.

BACKGROUND: There is considerable interest in the development of methods to efficiently identify all coding variants present in large sample sets of humans. There are three approaches possible: whole-genome sequencing, whole-exome sequencing using exon capture methods, and RNA-Seq. While whole-genome sequencing is the most complete, it remains sufficiently expensive that cost effective alternatives are important. RESULTS: Here we provide a systematic exploration of how well RNA-Seq can identify human coding variants by comparing variants identified through high coverage whole-genome sequencing to those identified by high coverage RNA-Seq in the same individual. This comparison allowed us to directly evaluate the sensitivity and specificity of RNA-Seq in identifying coding variants, and to evaluate how key parameters such as the degree of coverage and the expression levels of genes interact to influence performance. We find that although only 40% of exonic variants identified by whole genome sequencing were captured using RNA-Seq; this number rose to 81% when concentrating on genes known to be well-expressed in the source tissue. We also find that a high false positive rate can be problematic when working with RNA-Seq data, especially at higher levels of coverage. CONCLUSIONS: We conclude that as long as a tissue relevant to the trait under study is available and suitable quality control screens are implemented, RNA-Seq is a fast and inexpensive alternative approach for finding coding variants in genes with sufficiently high expression levels.

Metabolic programming and PDHK1 control CD4+ T cell subsets and inflammation.

Activation of CD4+ T cells results in rapid proliferation and differentiation into effector and regulatory subsets. CD4+ effector T cell (Teff) (Th1 and Th17) and Treg subsets are metabolically distinct, yet the specific metabolic differences that modify T cell populations are uncertain. Here, we evaluated CD4+ T cell populations in murine models and determined that inflammatory Teffs maintain high expression of glycolytic genes and rely on high glycolytic rates, while Tregs are oxidative and require mitochondrial electron transport to proliferate, differentiate, and survive. Metabolic profiling revealed that pyruvate dehydrogenase (PDH) is a key bifurcation point between T cell glycolytic and oxidative metabolism. PDH function is inhibited by PDH kinases (PDHKs). PDHK1 was expressed in Th17 cells, but not Th1 cells, and at low levels in Tregs, and inhibition or knockdown of PDHK1 selectively suppressed Th17 cells and increased Tregs. This alteration in the CD4+ T cell populations was mediated in part through ROS, as N-acetyl cysteine (NAC) treatment restored Th17 cell generation. Moreover, inhibition of PDHK1 modulated immunity and protected animals against experimental autoimmune encephalomyelitis, decreasing Th17 cells and increasing Tregs. Together, these data show that CD4+ subsets utilize and require distinct metabolic programs that can be targeted to control specific T cell populations in autoimmune and inflammatory diseases.

Structural and Functional Analysis of the Caspase –dependent and –independent Domains of the X-linked Inhibitor of Apoptosis Protein in Inflammatory Breast Cancer Tumor Biology

Inflammatory breast cancer (IBC) is an extremely rare but highly aggressive form of breast cancer characterized by the rapid development of therapeutic resistance leading to particularly poor survival. Our previous work focused on the elucidation of factors that mediate therapeutic resistance in IBC and identified increased expression of the anti-apoptotic protein, X-linked inhibitor of apoptosis protein (XIAP), to correlate with the development of resistance to chemotherapeutics. Although XIAP is classically thought of as an inhibitor of caspase activation, multiple studies have revealed that XIAP can also function as a signaling intermediate in numerous pathways. Based on preliminary evidence revealing high expression of XIAP in pre-treatment IBC cells rather than only subsequent to the development of resistance, we hypothesized that XIAP could play an important signaling role in IBC pathobiology outside of its heavily published apoptotic inhibition function. Further, based on our discovery of inhibition of chemotherapeutic efficacy, we postulated that XIAP overexpression might also play a role in resistance to other forms of therapy, such as immunotherapy. Finally, we posited that targeting of specific redox adaptive mechanisms, which are observed to be a significant barrier to successful treatment of IBC, could overcome therapeutic resistance and enhance the efficacy of chemo-, radio-, and immuno- therapies. To address these hypotheses our objectives were: 1. to determine a role for XIAP in IBC pathobiology and to elucidate the upstream regulators and downstream effectors of XIAP; 2. to evaluate and describe a role for XIAP in the inhibition of immunotherapy; and 3. to develop and characterize novel redox modulatory strategies that target identified mechanisms to prevent or reverse therapeutic resistance.

Using various genomic and proteomic approaches, combined with analysis of cellular viability, proliferation, and growth parameters both in vitro and in vivo, we demonstrate that XIAP plays a central role in both IBC pathobiology in a manner mostly independent of its role as a caspase-binding protein. Modulation of XIAP expression in cells derived from patients prior to any therapeutic intervention significantly altered key aspects IBC biology including, but not limited to: IBC-specific gene signatures; the tumorigenic capacity of tumor cells; and the metastatic phenotype of IBC, all of which are revealed to functionally hinge on XIAP-mediated NFκB activation, a robust molecular determinant of IBC. Identification of the mechanism of XIAP-mediated NFκB activation led to the characterization of novel peptide-based antagonist which was further used to identify that increased NFκB activation was responsible for redox adaptation previously observed in therapy-resistant IBC cells. Lastly, we describe the targeting of this XIAP-NFκB-ROS axis using a novel redox modulatory strategy both in vitro and in vivo. Together, the data presented here characterize a novel and crucial role for XIAP both in therapeutic resistance and the pathobiology of IBC; these results confirm our previous work in acquired therapeutic resistance and establish the feasibility of targeting XIAP-NFκB and the redox adaptive phenotype of IBC as a means to enhance survival of patients.

Role of ABL Family Kinases in Breast Cancer

The ABL family of non-receptor tyrosine kinases, ABL1 (also known as c-ABL) and ABL2 (also known as Arg), links diverse extracellular stimuli to signaling pathways that control cell growth, survival, adhesion, migration and invasion. ABL tyrosine kinases play an oncogenic role in human leukemias. However, the role of ABL kinases in solid tumors including breast cancer progression and metastasis is just emerging.

To evaluate whether ABL family kinases are involved in breast cancer development and metastasis, we first analyzed genomic data from large-scale screen of breast cancer patients. We found that ABL kinases are up-regulated in invasive breast cancer patients and high expression of ABL kinases correlates with poor prognosis and early metastasis. Using xenograft mouse models combined with genetic and pharmacological approaches, we demonstrated that ABL kinases are required for regulating breast cancer progression and metastasis to the bone. Using next generation sequencing and bioinformatics analysis, we uncovered a critical role for ABL kinases in promoting multiple oncogenic pathways including TAZ and STAT5 signaling networks and the epithelial to mesenchymal transition (EMT). These findings revealed a role for ABL kinases in regulating breast cancer tumorigenesis and bone metastasis and provide a rationale for targeting breast tumors with ABL-specific inhibitors.

X-linked inhibitor of apoptosis protein mediates tumor cell resistance to antibody-dependent cellular cytotoxicity.

Inflammatory breast cancer (IBC) is the deadliest, distinct subtype of breast cancer. High expression of epidermal growth factor receptors [EGFR or human epidermal growth factor receptor 2 (HER2)] in IBC tumors has prompted trials of anti-EGFR/HER2 monoclonal antibodies to inhibit oncogenic signaling; however, de novo and acquired therapeutic resistance is common. Another critical function of these antibodies is to mediate antibody-dependent cellular cytotoxicity (ADCC), which enables immune effector cells to engage tumors and deliver granzymes, activating executioner caspases. We hypothesized that high expression of anti-apoptotic molecules in tumors would render them resistant to ADCC. Herein, we demonstrate that the most potent caspase inhibitor, X-linked inhibitor of apoptosis protein (XIAP), overexpressed in IBC, drives resistance to ADCC mediated by cetuximab (anti-EGFR) and trastuzumab (anti-HER2). Overexpression of XIAP in parental IBC cell lines enhances resistance to ADCC; conversely, targeted downregulation of XIAP in ADCC-resistant IBC cells renders them sensitive. As hypothesized, this ADCC resistance is in part a result of the ability of XIAP to inhibit caspase activity; however, we also unexpectedly found that resistance was dependent on XIAP-mediated, caspase-independent suppression of reactive oxygen species (ROS) accumulation, which otherwise occurs during ADCC. Transcriptome analysis supported these observations by revealing modulation of genes involved in immunosuppression and oxidative stress response in XIAP-overexpressing, ADCC-resistant cells. We conclude that XIAP is a critical modulator of ADCC responsiveness, operating through both caspase-dependent and -independent mechanisms. These results suggest that strategies targeting the effects of XIAP on caspase activation and ROS suppression have the potential to enhance the activity of monoclonal antibody-based immunotherapy.

Morphing low-affinity ligands into high-avidity nanoparticles by thermally triggered self-assembly of a genetically encoded polymer.

Multivalency is the increase in avidity resulting from the simultaneous interaction of multiple ligands with multiple receptors. This phenomenon, seen in antibody-antigen and virus-cell membrane interactions, is useful in designing bioinspired materials for targeted delivery of drugs or imaging agents. While increased avidity offered by multivalent targeting is attractive, it can also promote nonspecific receptor interaction in nontarget tissues, reducing the effectiveness of multivalent targeting. Here, we present a thermal targeting strategy--dynamic affinity modulation (DAM)--using elastin-like polypeptide diblock copolymers (ELP(BC)s) that self-assemble from a low-affinity to high-avidity state by a tunable thermal "switch", thereby restricting activity to the desired site of action. We used an in vitro cell binding assay to investigate the effect of the thermally triggered self-assembly of these ELP(BC)s on their receptor-mediated binding and cellular uptake. The data presented herein show that (1) ligand presentation does not disrupt ELP(BC) self-assembly; (2) both multivalent ligand presentation and upregulated receptor expression are needed for receptor-mediated interaction; (3) increased size of the hydrophobic segment of the block copolymer promotes multivalent interaction with membrane receptors, potentially due to changes in the nanoscale architecture of the micelle; and (4) nanoscale presentation of the ligand is important, as presentation of the ligand by micrometer-sized aggregates of an ELP showed a low level of binding/uptake by receptor-positive cells compared to its presentation on the corona of a micelle. These data validate the concept of thermally triggered DAM and provide rational design parameters for future applications of this technology for targeted drug delivery.

SplicerAV: a tool for mining microarray expression data for changes in RNA processing.

BACKGROUND: Over the past two decades more than fifty thousand unique clinical and biological samples have been assayed using the Affymetrix HG-U133 and HG-U95 GeneChip microarray platforms. This substantial repository has been used extensively to characterize changes in gene expression between biological samples, but has not been previously mined en masse for changes in mRNA processing. We explored the possibility of using HG-U133 microarray data to identify changes in alternative mRNA processing in several available archival datasets. RESULTS: Data from these and other gene expression microarrays can now be mined for changes in transcript isoform abundance using a program described here, SplicerAV. Using in vivo and in vitro breast cancer microarray datasets, SplicerAV was able to perform both gene and isoform specific expression profiling within the same microarray dataset. Our reanalysis of Affymetrix U133 plus 2.0 data generated by in vitro over-expression of HRAS, E2F3, beta-catenin (CTNNB1), SRC, and MYC identified several hundred oncogene-induced mRNA isoform changes, one of which recognized a previously unknown mechanism of EGFR family activation. Using clinical data, SplicerAV predicted 241 isoform changes between low and high grade breast tumors; with changes enriched among genes coding for guanyl-nucleotide exchange factors, metalloprotease inhibitors, and mRNA processing factors. Isoform changes in 15 genes were associated with aggressive cancer across the three breast cancer datasets. CONCLUSIONS: Using SplicerAV, we identified several hundred previously uncharacterized isoform changes induced by in vitro oncogene over-expression and revealed a previously unknown mechanism of EGFR activation in human mammary epithelial cells. We analyzed Affymetrix GeneChip data from over 400 human breast tumors in three independent studies, making this the largest clinical dataset analyzed for en masse changes in alternative mRNA processing. The capacity to detect RNA isoform changes in archival microarray data using SplicerAV allowed us to carry out the first analysis of isoform specific mRNA changes directly associated with cancer survival.

Boosting high-intensity focused ultrasound-induced anti-tumor immunity using a sparse-scan strategy that can more effectively promote dendritic cell maturation.

BACKGROUND: The conventional treatment protocol in high-intensity focused ultrasound (HIFU) therapy utilizes a dense-scan strategy to produce closely packed thermal lesions aiming at eradicating as much tumor mass as possible. However, this strategy is not most effective in terms of inducing a systemic anti-tumor immunity so that it cannot provide efficient micro-metastatic control and long-term tumor resistance. We have previously provided evidence that HIFU may enhance systemic anti-tumor immunity by in situ activation of dendritic cells (DCs) inside HIFU-treated tumor tissue. The present study was conducted to test the feasibility of a sparse-scan strategy to boost HIFU-induced anti-tumor immune response by more effectively promoting DC maturation. METHODS: An experimental HIFU system was set up to perform tumor ablation experiments in subcutaneous implanted MC-38 and B16 tumor with dense- or sparse-scan strategy to produce closely-packed or separated thermal lesions. DCs infiltration into HIFU-treated tumor tissues was detected by immunohistochemistry and flow cytometry. DCs maturation was evaluated by IL-12/IL-10 production and CD80/CD86 expression after co-culture with tumor cells treated with different HIFU. HIFU-induced anti-tumor immune response was evaluated by detecting growth-retarding effects on distant re-challenged tumor and tumor-specific IFN-gamma-secreting cells in HIFU-treated mice. RESULTS: HIFU exposure raised temperature up to 80 degrees centigrade at beam focus within 4 s in experimental tumors and led to formation of a well-defined thermal lesion. The infiltrated DCs were recruited to the periphery of lesion, where the peak temperature was only 55 degrees centigrade during HIFU exposure. Tumor cells heated to 55 degrees centigrade in 4-s HIFU exposure were more effective to stimulate co-cultured DCs to mature. Sparse-scan HIFU, which can reserve 55 degrees-heated tumor cells surrounding the separated lesions, elicited an enhanced anti-tumor immune response than dense-scan HIFU, while their suppressive effects on the treated primary tumor were maintained at the same level. Flow cytometry analysis showed that sparse-scan HIFU was more effective than dense-scan HIFU in enhancing DC infiltration into tumor tissues and promoting their maturation in situ. CONCLUSION: Optimizing scan strategy is a feasible way to boost HIFU-induced anti-tumor immunity by more effectively promoting DC maturation.

Diagnosis of partial body radiation exposure in mice using peripheral blood gene expression profiles.

In the event of a terrorist-mediated attack in the United States using radiological or improvised nuclear weapons, it is expected that hundreds of thousands of people could be exposed to life-threatening levels of ionizing radiation. We have recently shown that genome-wide expression analysis of the peripheral blood (PB) can generate gene expression profiles that can predict radiation exposure and distinguish the dose level of exposure following total body irradiation (TBI). However, in the event a radiation-mass casualty scenario, many victims will have heterogeneous exposure due to partial shielding and it is unknown whether PB gene expression profiles would be useful in predicting the status of partially irradiated individuals. Here, we identified gene expression profiles in the PB that were characteristic of anterior hemibody-, posterior hemibody- and single limb-irradiation at 0.5 Gy, 2 Gy and 10 Gy in C57Bl6 mice. These PB signatures predicted the radiation status of partially irradiated mice with a high level of accuracy (range 79-100%) compared to non-irradiated mice. Interestingly, PB signatures of partial body irradiation were poorly predictive of radiation status by site of injury (range 16-43%), suggesting that the PB molecular response to partial body irradiation was anatomic site specific. Importantly, PB gene signatures generated from TBI-treated mice failed completely to predict the radiation status of partially irradiated animals or non-irradiated controls. These data demonstrate that partial body irradiation, even to a single limb, generates a characteristic PB signature of radiation injury and thus may necessitate the use of multiple signatures, both partial body and total body, to accurately assess the status of an individual exposed to radiation.

Specialized motor-driven dusp1 expression in the song systems of multiple lineages of vocal learning birds.

Mechanisms for the evolution of convergent behavioral traits are largely unknown. Vocal learning is one such trait that evolved multiple times and is necessary in humans for the acquisition of spoken language. Among birds, vocal learning is evolved in songbirds, parrots, and hummingbirds. Each time similar forebrain song nuclei specialized for vocal learning and production have evolved. This finding led to the hypothesis that the behavioral and neuroanatomical convergences for vocal learning could be associated with molecular convergence. We previously found that the neural activity-induced gene dual specificity phosphatase 1 (dusp1) was up-regulated in non-vocal circuits, specifically in sensory-input neurons of the thalamus and telencephalon; however, dusp1 was not up-regulated in higher order sensory neurons or motor circuits. Here we show that song motor nuclei are an exception to this pattern. The song nuclei of species from all known vocal learning avian lineages showed motor-driven up-regulation of dusp1 expression induced by singing. There was no detectable motor-driven dusp1 expression throughout the rest of the forebrain after non-vocal motor performance. This pattern contrasts with expression of the commonly studied activity-induced gene egr1, which shows motor-driven expression in song nuclei induced by singing, but also motor-driven expression in adjacent brain regions after non-vocal motor behaviors. In the vocal non-learning avian species, we found no detectable vocalizing-driven dusp1 expression in the forebrain. These findings suggest that independent evolutions of neural systems for vocal learning were accompanied by selection for specialized motor-driven expression of the dusp1 gene in those circuits. This specialized expression of dusp1 could potentially lead to differential regulation of dusp1-modulated molecular cascades in vocal learning circuits.

Sensitive and precise quantification of insulin-like mRNA expression in Caenorhabditis elegans.

Insulin-like signaling regulates developmental arrest, stress resistance and lifespan in the nematode Caenorhabditis elegans. However, the genome encodes 40 insulin-like peptides, and the regulation and function of individual peptides is largely uncharacterized. We used the nCounter platform to measure mRNA expression of all 40 insulin-like peptides as well as the insulin-like receptor daf-2, its transcriptional effector daf-16, and the daf-16 target gene sod-3. We validated the platform using 53 RNA samples previously characterized by high density oligonucleotide microarray analysis. For this set of genes and the standard nCounter protocol, sensitivity and precision were comparable between the two platforms. We optimized conditions of the nCounter assay by varying the mass of total RNA used for hybridization, thereby increasing sensitivity up to 50-fold and reducing the median coefficient of variation as much as 4-fold. We used deletion mutants to demonstrate specificity of the assay, and we used optimized conditions to assay insulin-like gene expression throughout the C. elegans life cycle. We detected expression for nearly all insulin-like genes and find that they are expressed in a variety of distinct patterns suggesting complexity of regulation and specificity of function. We identified insulin-like genes that are specifically expressed during developmental arrest, larval development, adulthood and embryogenesis. These results demonstrate that the nCounter platform provides a powerful approach to analyzing insulin-like gene expression dynamics, and they suggest hypotheses about the function of individual insulin-like genes.

A complex intronic enhancer regulates expression of the CFTR gene by direct interaction with the promoter.

Genes can maintain spatiotemporal expression patterns by long-range interactions between cis-acting elements. The cystic fibrosis transmembrane conductance regulator gene (CFTR) is expressed primarily in epithelial cells. An element located within a DNase I-hypersensitive site (DHS) 10 kb into the first intron was previously shown to augment CFTR promoter activity in a tissue-specific manner. Here, we reveal the mechanism by which this element influences CFTR transcription. We employed a high-resolution method of mapping DHS using tiled microarrays to accurately locate the intron 1 DHS. Transfection of promoter-reporter constructs demonstrated that the element displays classical tissue-specific enhancer properties and can independently recruit factors necessary for transcription initiation. In vitro DNase I footprinting analysis identified a protected region that corresponds to a conserved, predicted binding site for hepatocyte nuclear factor 1 (HNF1). We demonstrate by electromobility shift assays (EMSA) and chromatin immunoprecipitation (ChIP) that HNF1 binds to this element both in vitro and in vivo. Moreover, using chromosome conformation capture (3C) analysis, we show that this element interacts with the CFTR promoter in CFTR-expressing cells. These data provide the first insight into the three- dimensional (3D) structure of the CFTR locus and confirm the contribution of intronic cis-acting elements to the regulation of CFTR gene expression.

A noisy linear map underlies oscillations in cell size and gene expression in bacteria.

During bacterial growth, a cell approximately doubles in size before division, after which it splits into two daughter cells. This process is subjected to the inherent perturbations of cellular noise and thus requires regulation for cell-size homeostasis. The mechanisms underlying the control and dynamics of cell size remain poorly understood owing to the difficulty in sizing individual bacteria over long periods of time in a high-throughput manner. Here we measure and analyse long-term, single-cell growth and division across different Escherichia coli strains and growth conditions. We show that a subset of cells in a population exhibit transient oscillations in cell size with periods that stretch across several (more than ten) generations. Our analysis reveals that a simple law governing cell-size control-a noisy linear map-explains the origins of these cell-size oscillations across all strains. This noisy linear map implements a negative feedback on cell-size control: a cell with a larger initial size tends to divide earlier, whereas one with a smaller initial size tends to divide later. Combining simulations of cell growth and division with experimental data, we demonstrate that this noisy linear map generates transient oscillations, not just in cell size, but also in constitutive gene expression. Our work provides new insights into the dynamics of bacterial cell-size regulation with implications for the physiological processes involved.

PenPC: A two-step approach to estimate the skeletons of high-dimensional directed acyclic graphs.

Estimation of the skeleton of a directed acyclic graph (DAG) is of great importance for understanding the underlying DAG and causal effects can be assessed from the skeleton when the DAG is not identifiable. We propose a novel method named PenPC to estimate the skeleton of a high-dimensional DAG by a two-step approach. We first estimate the nonzero entries of a concentration matrix using penalized regression, and then fix the difference between the concentration matrix and the skeleton by evaluating a set of conditional independence hypotheses. For high-dimensional problems where the number of vertices p is in polynomial or exponential scale of sample size n, we study the asymptotic property of PenPC on two types of graphs: traditional random graphs where all the vertices have the same expected number of neighbors, and scale-free graphs where a few vertices may have a large number of neighbors. As illustrated by extensive simulations and applications on gene expression data of cancer patients, PenPC has higher sensitivity and specificity than the state-of-the-art method, the PC-stable algorithm.