Modeling phenotypic heterogeneity towards evolutionarily inspired osteosarcoma therapy.

Osteosarcoma is the most common bone sarcoma in children and young adults. While universally delivered, chemotherapy only benefits roughly half of patients with localized disease. Increasingly, intratumoral heterogeneity is recognized as a source of therapeutic resistance. In this study, we develop and evaluate an in vitro model of osteosarcoma heterogeneity based on phenotype and genotype. Cancer cell populations vary in their environment-specific growth rates and in their sensitivity to chemotherapy. We present the genotypic and phenotypic characterization of an osteosarcoma cell line panel with a focus on co-cultures of the most phenotypically divergent cell lines, 143B and SAOS2. Modest environmental (pH, glutamine) or chemical perturbations dramatically shift the success and composition of cell lines. We demonstrate that in nutrient rich culture conditions 143B outcompetes SAOS2. But, under nutrient deprivation or conventional chemotherapy, SAOS2 growth can be favored in spheroids. Importantly, when the simplest heterogeneity state is evaluated, a two-cell line coculture, perturbations that affect the faster growing cell line have only a modest effect on final spheroid size. Thus the only evaluated therapies to eliminate the spheroids were by switching therapies from a first strike to a second strike. This extensively characterized, widely available system, can be modeled and scaled to allow for improved strategies to anticipate resistance in osteosarcoma due to heterogeneity.

Managing a Large-Scale Multiomics Project: A Team Science Case Study in Proteogenomics.

Highly collaborative scientists are often called on to extend their expertise to different types of projects and to expand the scope and scale of projects well beyond their previous experience. For a large-scale project involving "big data" to be successful, several different aspects of the research plan need to be developed and tested, which include but are not limited to the experimental design, sample collection, sample preparation, metadata recording, technical capability, data acquisition, approaches for data analysis, methods for integration of different data types, recruitment of additional expertise as needed to guide the project, and strategies for clear communication throughout the project. To capture this process, we describe an example project in proteogenomics that built on our collective expertise and experience. Key steps included definition of hypotheses, identification of an appropriate clinical cohort, pilot projects to assess feasibility, refinement of experimental designs, and extensive discussions involving the research team throughout the process. The goal of this chapter is to provide the reader with a set of guidelines to support development of other large-scale multiomics projects.

Fabrication of All-Solid Organic Electrochromic Devices on Absorptive Paper Substrates Utilizing a Simplified Lateral Architecture.

Poly(3,4-ethylenedioxythiophene) doped with the polymer anion poly(styrenesulfonate), PEDOT:PSS, is a common electrochromic material used in the preparation of electrochromic devices (ECDs). In this paper, the PEDOT:PSS doped with a solvent was used both as the electrode and the electrochromic functional layer for fabrication of ECDs on absorptive paper surfaces. The doped PEDOT:PSS dispersion was assessed for the film-forming evenness, sheet resistance and conductivity, and the performance of prepared ECDs for their color contrast and switching dynamics. The ECD performance is discussed in relation to the absorptive characteristics of the substrates. The results indicate that it is feasible to prepare ECDs onto absorptive substrates, despite the partial polymer material imbibition into them. The extent of polymer absorption influences the ECD performance: an increased absorption reduces the color contrast but speeds up the color switching. The electrochemical properties of the used solid electrolyte were found to be crucial for functioning of the ECDs. Insufficient ion transport and associated high resistance led to failure of the devices.

The genomic landscape of undifferentiated embryonal sarcoma of the liver is typified by C19MC structural rearrangement and overexpression combined with TP53 mutation or loss.

Undifferentiated embryonal sarcoma of the liver (UESL) is a rare and aggressive malignancy. Though the molecular underpinnings of this cancer have been largely unexplored, recurrent chromosomal breakpoints affecting a noncoding region on chr19q13, which includes the chromosome 19 microRNA cluster (C19MC), have been reported in several cases. We performed comprehensive molecular profiling on samples from 14 patients diagnosed with UESL. Congruent with prior reports, we identified structural variants in chr19q13 in 10 of 13 evaluable tumors. From whole transcriptome sequencing, we observed striking expressional activity of the entire C19MC region. Concordantly, in 7 of 7 samples undergoing miRNAseq, we observed hyperexpression of the miRNAs within this cluster to levels >100 fold compared to matched normal tissue or a non-C19MC amplified cancer cell line. Concurrent TP53 mutation or copy number loss was identified in all evaluable tumors with evidence of C19MC overexpression. We find that C19MC miRNAs exhibit significant negative correlation to TP53 regulatory miRNAs and K-Ras regulatory miRNAs. Using RNA-seq we identified that pathways relevant to cellular differentiation as well as mRNA translation machinery are transcriptionally enriched in UESL. In summary, utilizing a combination of next-generation sequencing and high-density arrays we identify the combination of C19MC hyperexpression via chromosomal structural event with TP53 mutation or loss as highly recurrent genomic features of UESL.

Proteogenomic landscape of squamous cell lung cancer.

How genomic and transcriptomic alterations affect the functional proteome in lung cancer is not fully understood. Here, we integrate DNA copy number, somatic mutations, RNA-sequencing, and expression proteomics in a cohort of 108 squamous cell lung cancer (SCC) patients. We identify three proteomic subtypes, two of which (Inflamed, Redox) comprise 87% of tumors. The Inflamed subtype is enriched with neutrophils, B-cells, and monocytes and expresses more PD-1. Redox tumours are enriched for oxidation-reduction and glutathione pathways and harbor more NFE2L2/KEAP1 alterations and copy gain in the 3q2 locus. Proteomic subtypes are not associated with patient survival. However, B-cell-rich tertiary lymph node structures, more common in Inflamed, are associated with better survival. We identify metabolic vulnerabilities (TP63, PSAT1, and TFRC) in Redox. Our work provides a powerful resource for lung SCC biology and suggests therapeutic opportunities based on redox metabolism and immune cell infiltrates.

Application of Bayesian predictive probability for interim futility analysis in single-arm phase II trial.

BACKGROUND: Bayesian predictive probability design, with a binary endpoint, is gaining attention for the phase II trial due to its innovative strategy. To make the Bayesian design more accessible, we elucidate this Bayesian approach with a R package to streamline a statistical plan, so biostatisticians and clinicians can easily integrate the design into clinical trial. METHODS: We utilize a Bayesian framework using Bayesian posterior probability and predictive probability to build a R package and develop a statistical plan for the trial design. With pre-defined sample sizes, the approach employs the posterior probability with a threshold to calculate the minimum number of responders needed at end of the study to claim efficacy. Then the predictive probability is applied to evaluate future success at interim stages and form stopping rule at each stage. RESULTS: An R package, 'BayesianPredictiveFutility', with associated graphical interface is developed for easy utilization of the trial design. The statistical tool generates a professional statistical plan with comprehensive results including a summary, details of study design, a series of tables and figures from stopping boundary for futility, Bayesian predictive probability, performance [probability of early termination (PET), type I error, and power], PET at each interim analysis, sensitivity analysis for predictive probability, posterior probability, sample size, and beta prior distribution. The statistical plan presents the methodology in a readable language fashion while preserving rigorous statistical arguments. The output formats (Word or PDF) are available to communicate with physicians or to be incorporated in the trial protocol. Two clinical trials in lung cancer are used to demonstrate its usefulness. CONCLUSIONS: Bayesian predictive probability method presents a flexible design in clinical trial. The statistical tool brings an added value to broaden the application.

Redox-Neutral Dual Functionalization of Electron-Deficient Alkenes.

Visible-light photoredox catalysis has been utilized in a new multicomponent reaction forming ß-functionalized δ-diketones under mild conditions in an operationally convenient manner. Single-electron reduction of in situ generated carboxylic acid derivatives forms acyl radicals that react further via 1,2-acylalkylation of olefins in an intermolecular, three-components cascade reaction, giving valuable synthetic entities from readily available starting materials. A diverse set of substrates has been used, demonstrating robust methodology with broad substrate scope.

Environmentally friendly transistors and circuits on paper.

We created environmentally friendly low-voltage, ion-modulated transistors (IMTs) that can be fabricated successfully on a paper substrate. A range of ionic liquids (ILs) based on choline chloride (ChoCl) were used as the electrolytic layer in the IMTs. Different organic compounds were mixed with ChoCl to create solution-processable deep eutectic mixtures that are liquid or semiliquid at room temperature. In the final, solid version of the IMT, the ILs are also solidified by using a commercial binder to create printable transistor structures The semiconductor layer in the IMT is also substituted with a blend of the original semiconductor and a biodegradable polymer insulator. This reduces the amount of expensive and potentially harmful semiconductor used, and it also provides increased transistor performance, especially increasing the device switching speed. These environmentally friendly IMTs are then used to create ring oscillators, logic gates, and memories on paper.

Impact of humidity on functionality of on-paper printed electronics.

A multilayer coated paper substrate, combining barrier and printability properties was manufactured utilizing a pilot-scale slide curtain coating technique. The coating structure consists of a thin mineral pigment layer coated on top of a barrier layer. The surface properties, i.e. smoothness and surface porosity, were adjusted by the choice of calendering parameters. The influence of surface properties on the fine line printability and conductivity of inkjet-printed silver lines was studied. Surface roughness played a significant role when printing narrow lines, increasing the risk of defects and discontinuities, whereas for wider lines the influence of surface roughness was less critical. A smooth, calendered surface resulted in finer line definition, i.e. less edge raggedness. Dimensional stability and its influence on substrate surface properties as well as on the functionality of conductive tracks and transistors were studied by exposure to high/low humidity cycles. The barrier layer of the multilayer coated paper reduced the dimensional changes and surface roughness increase caused by humidity and helped maintain the conductivity of the printed tracks. Functionality of a printed transistor during a short, one hour humidity cycle was maintained, but a longer exposure to humidity destroyed the non-encapsulated transistor.

An impedimetric study of DNA hybridization on paper-supported inkjet-printed gold electrodes.

In this study, two different supramolecular recognition architectures for impedimetric detection of DNA hybridization have been formed on disposable paper-supported inkjet-printed gold electrodes. The gold electrodes were fabricated using a gold nanoparticle based ink. The first recognition architecture consists of subsequent layers of biotinylated self-assembly monolayer (SAM), streptavidin and biotinylated DNA probe. The other recognition architecture is constructed by immobilization of thiol-functionalized DNA probe (HS-DNA) and subsequent backfill with 11-mercapto-1-undecanol (MUOH) SAM. The binding capacity and selectivity of the recognition architectures were examined by surface plasmon resonance (SPR) measurements. SPR results showed that the HS-DNA/MUOH system had a higher binding capacity for the complementary DNA target. Electrochemical impedance spectroscopy (EIS) measurements showed that the hybridization can be detected with impedimetric spectroscopy in picomol range for both systems. EIS signal indicated a good selectivity for both recognition architectures, whereas SPR showed very high unspecific binding for the HS-DNA/MUOH system. The factors affecting the impedance signal were interpreted in terms of the complexity of the supramolecular architecture. The more complex architecture acts as a less ideal capacitive sensor and the impedance signal is dominated by the resistive elements.

Synthesis, pharmacological evaluation and QSAR modeling of mono-substituted 4-phenylpiperidines and 4-phenylpiperazines.

A series of mono-substituted 4-phenylpiperidines and -piperazines have been synthesized and their effects on the dopaminergic system tested in vivo. The structure activity relationship (SAR) revealed that the position and physicochemical character of the aromatic substituent proved to be critical for the levels of 3,4-dihydroxyphenylacetic acid (DOPAC) in the brain of freely moving rats. In order to investigate how the structural properties of these compounds affect the response, a set of tabulated and calculated physicochemical descriptors were modeled against the in vivo effects using partial least square (PLS) regression. Furthermore, the binding affinities to the dopamine D2 (DA D2) receptor and monoamine oxidase A (MAO A) enzyme were determined for a chosen subset and QSAR models using the same descriptors as in the in vivo model were produced to investigate the mechanisms leading to the observed DOPAC response. These models, in combination with a strong correlation between the levels of striatal DOPAC and the affinities to DA D2 and MAO A, provides a comprehensive understanding of the biological response for compounds in this class.

Sequence signatures involved in targeting the Male-Specific Lethal complex to X-chromosomal genes in Drosophila melanogaster.

BACKGROUND: In Drosophila melanogaster, the dosage-compensation system that equalizes X-linked gene expression between males and females, thereby assuring that an appropriate balance is maintained between the expression of genes on the X chromosome(s) and the autosomes, is at least partially mediated by the Male-Specific Lethal (MSL) complex. This complex binds to genes with a preference for exons on the male X chromosome with a 3' bias, and it targets most expressed genes on the X chromosome. However, a number of genes are expressed but not targeted by the complex. High affinity sites seem to be responsible for initial recruitment of the complex to the X chromosome, but the targeting to and within individual genes is poorly understood. RESULTS: We have extensively examined X chromosome sequence variation within five types of gene features (promoters, 5' UTRs, coding sequences, introns, 3' UTRs) and intergenic sequences, and assessed its potential involvement in dosage compensation. Presented results show that: the X chromosome has a distinct sequence composition within its gene features; some of the detected variation correlates with genes targeted by the MSL-complex; the insulator protein BEAF-32 preferentially binds upstream of MSL-bound genes; BEAF-32 and MOF co-localizes in promoters; and that bound genes have a distinct sequence composition that shows a 3' bias within coding sequence. CONCLUSIONS: Although, many strongly bound genes are close to a high affinity site neither our promoter motif nor our coding sequence signatures show any correlation to HAS. Based on the results presented here, we believe that there are sequences in the promoters and coding sequences of targeted genes that have the potential to direct the secondary spreading of the MSL-complex to nearby genes.

Coexpression network analysis in abdominal and gluteal adipose tissue reveals regulatory genetic loci for metabolic syndrome and related phenotypes.

Metabolic Syndrome (MetS) is highly prevalent and has considerable public health impact, but its underlying genetic factors remain elusive. To identify gene networks involved in MetS, we conducted whole-genome expression and genotype profiling on abdominal (ABD) and gluteal (GLU) adipose tissue, and whole blood (WB), from 29 MetS cases and 44 controls. Co-expression network analysis for each tissue independently identified nine, six, and zero MetS-associated modules of coexpressed genes in ABD, GLU, and WB, respectively. Of 8,992 probesets expressed in ABD or GLU, 685 (7.6%) were expressed in ABD and 51 (0.6%) in GLU only. Differential eigengene network analysis of 8,256 shared probesets detected 22 shared modules with high preservation across adipose depots (D(ABD-GLU) = 0.89), seven of which were associated with MetS (FDR P<0.01). The strongest associated module, significantly enriched for immune response-related processes, contained 94/620 (15%) genes with inter-depot differences. In an independent cohort of 145/141 twins with ABD and WB longitudinal expression data, median variability in ABD due to familiality was greater for MetS-associated versus un-associated modules (ABD: 0.48 versus 0.18, P = 0.08; GLU: 0.54 versus 0.20, P = 7.8×10(-4)). Cis-eQTL analysis of probesets associated with MetS (FDR P<0.01) and/or inter-depot differences (FDR P<0.01) provided evidence for 32 eQTLs. Corresponding eSNPs were tested for association with MetS-related phenotypes in two GWAS of >100,000 individuals; rs10282458, affecting expression of RARRES2 (encoding chemerin), was associated with body mass index (BMI) (P = 6.0×10(-4)); and rs2395185, affecting inter-depot differences of HLA-DRB1 expression, was associated with high-density lipoprotein (P = 8.7×10(-4)) and BMI-adjusted waist-to-hip ratio (P = 2.4×10(-4)). Since many genes and their interactions influence complex traits such as MetS, integrated analysis of genotypes and coexpression networks across multiple tissues relevant to clinical traits is an efficient strategy to identify novel associations.

Synthesis and evaluation of a set of para-substituted 4-phenylpiperidines and 4-phenylpiperazines as monoamine oxidase (MAO) inhibitors.

A series of para-substituted 4-phenylpiperidines/piperazines have been synthesized and their affinity to recombinant rat cerebral cortex monoamine oxidases A (MAO A) and B (MAO B) determined. Para-substituents with low dipole moment increased the affinity to MAO A, whereas groups with high dipole moment yielded compounds with no or weak affinity. In contrast, the properties affecting MAO B affinity were the polarity and bulk of the para-substituent, with large hydrophobic substituents producing compounds with high MAO B affinity. In addition, these compounds were tested in freely moving rats and the effect on the post-mortem neurochemistry was measured. A linear correlation was demonstrated between the affinity for MAO A, but not MAO B, and the levels of 3,4-dihydroxyphenylacetic acid (DOPAC) and 3-methoxytyramine (3-MT) in the striatum.

A method for bidirectional solution exchange--"liquid bullet" applications of acetylcholine to α7 nicotinic receptors.

Fast solution exchange techniques have revolutionized the study of synaptic transmission and promise to remain an important neuroscience research tool. Here we provide evidence for the hypothesis that using continuous, rapid transitions through an agonist solution can significantly increase the exchange rate around a cell by reducing the diffusion boundary at the membrane. This novel approach of rapid solution exchange during whole-cell recordings--described as a "liquid bullet" (LB) application--takes advantage of a bidirectional solution flow around the cell, allowing for a full solution exchange within a range of several milliseconds. An exchange rate (10-90% rise time) of about 2 ms could be achieved during both agonist application and washout. We recorded whole-cell currents from cells expressing the rapidly desensitizing α7 neuronal nicotinic receptor (NNR) subtype that exhibited very fast rise times of around 4-5 ms. We further demonstrated the advantages of a LB application over conventional methods by the ability of this method to elicit concentration-dependent responses for rapidly desensitizing compounds that were not measurable with conventional agonist applications. In addition, we illustrate the utility of this approach for frequency-based assays through fast, repeated agonist applications at frequencies of 1 Hz and 30 Hz. This approach could therefore be useful for the study of rapid agonist-receptor interactions that closely mimic the physiological conditions in the synaptic cleft during bursts of neuronal activity.

The use of genome-wide eQTL associations in lymphoblastoid cell lines to identify novel genetic pathways involved in complex traits.

The integrated analysis of genotypic and expression data for association with complex traits could identify novel genetic pathways involved in complex traits. We profiled 19,573 expression probes in Epstein-Barr virus-transformed lymphoblastoid cell lines (LCLs) from 299 twins and correlated these with 44 quantitative traits (QTs). For 939 expressed probes correlating with more than one QT, we investigated the presence of eQTL associations in three datasets of 57 CEU HapMap founders and 86 unrelated twins. Genome-wide association analysis of these probes with 2.2 m SNPs revealed 131 potential eQTLs (1,989 eQTL SNPs) overlapping between the HapMap datasets, five of which were in cis (58 eQTL SNPs). We then tested 535 SNPs tagging the eQTL SNPs, for association with the relevant QT in 2,905 twins. We identified nine potential SNP-QT associations (P<0.01) but none significantly replicated in five large consortia of 1,097-16,129 subjects. We also failed to replicate previous reported eQTL associations with body mass index, plasma low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides levels derived from lymphocytes, adipose and liver tissue. Our results and additional power calculations suggest that proponents may have been overoptimistic in the power of LCLs in eQTL approaches to elucidate regulatory genetic effects on complex traits using the small datasets generated to date. Nevertheless, larger tissue-specific expression data sets relevant to specific traits are becoming available, and should enable the adoption of similar integrated analyses in the near future.

Basic statistical analysis in genetic case-control studies.

This protocol describes how to perform basic statistical analysis in a population-based genetic association case-control study. The steps described involve the (i) appropriate selection of measures of association and relevance of disease models; (ii) appropriate selection of tests of association; (iii) visualization and interpretation of results; (iv) consideration of appropriate methods to control for multiple testing; and (v) replication strategies. Assuming no previous experience with software such as PLINK, R or Haploview, we describe how to use these popular tools for handling single-nucleotide polymorphism data in order to carry out tests of association and visualize and interpret results. This protocol assumes that data quality assessment and control has been performed, as described in a previous protocol, so that samples and markers deemed to have the potential to introduce bias to the study have been identified and removed. Study design, marker selection and quality control of case-control studies have also been discussed in earlier protocols. The protocol should take ~1 h to complete.

Enhanced performance of printed organic diodes using a thin interfacial barrier layer.

Printed, organic diodes with a thin organic interfacial layer forming a Schottky barrier were fabricated and characterized. Experiments indicated that the thickness of the barrier layer is <10 nm. The interfacial layer reduces the reverse current of the diode by 2 orders of magnitude without significantly affecting the forward characteristics above 1 V. As a result, printed organic diodes with a rectification ratio of 5 orders of magnitude were fabricated. The diodes enable applications where low reverse currents are needed.

Data quality control in genetic case-control association studies.

This protocol details the steps for data quality assessment and control that are typically carried out during case-control association studies. The steps described involve the identification and removal of DNA samples and markers that introduce bias. These critical steps are paramount to the success of a case-control study and are necessary before statistically testing for association. We describe how to use PLINK, a tool for handling SNP data, to perform assessments of failure rate per individual and per SNP and to assess the degree of relatedness between individuals. We also detail other quality-control procedures, including the use of SMARTPCA software for the identification of ancestral outliers. These platforms were selected because they are user-friendly, widely used and computationally efficient. Steps needed to detect and establish a disease association using case-control data are not discussed here. Issues concerning study design and marker selection in case-control studies have been discussed in our earlier protocols. This protocol, which is routinely used in our labs, should take approximately 8 h to complete.

Leprosy and the adaptation of human toll-like receptor 1.

Leprosy is an infectious disease caused by the obligate intracellular pathogen Mycobacterium leprae and remains endemic in many parts of the world. Despite several major studies on susceptibility to leprosy, few genomic loci have been replicated independently. We have conducted an association analysis of more than 1,500 individuals from different case-control and family studies, and observed consistent associations between genetic variants in both TLR1 and the HLA-DRB1/DQA1 regions with susceptibility to leprosy (TLR1 I602S, case-control P = 5.7 x 10(-8), OR = 0.31, 95% CI = 0.20-0.48, and HLA-DQA1 rs1071630, case-control P = 4.9 x 10(-14), OR = 0.43, 95% CI = 0.35-0.54). The effect sizes of these associations suggest that TLR1 and HLA-DRB1/DQA1 are major susceptibility genes in susceptibility to leprosy. Further population differentiation analysis shows that the TLR1 locus is extremely differentiated. The protective dysfunctional 602S allele is rare in Africa but expands to become the dominant allele among individuals of European descent. This supports the hypothesis that this locus may be under selection from mycobacteria or other pathogens that are recognized by TLR1 and its co-receptors. These observations provide insight into the long standing host-pathogen relationship between human and mycobacteria and highlight the key role of the TLR pathway in infectious diseases.