PubMed Health⌕ Search

Biomedical subjects

Michael Hubank

Publications and source records attributed to Michael Hubank.

10 recordsLinked to original sources

Evaluating the analytical validity of circulating tumor DNA sequencing assays for precision oncology.

Circulating tumor DNA (ctDNA) sequencing is being rapidly adopted in precision oncology, but the accuracy, sensitivity and reproducibility of ctDNA assays is poorly understood. Here we report the findings of a multi-site, cross-platform evaluation of the analytical performance of five industry-leading ctDNA assays. We evaluated each stage of the ctDNA sequencing workflow with simulations, synthetic DNA spike-in experiments and proficiency testing on standardized, cell-line-derived reference samples. Above 0.5% variant allele frequency, ctDNA mutations were detected with high sensitivity, precision and reproducibility by all five assays, whereas, below this limit, detection became unreliable and varied widely between assays, especially when input material was limited. Missed mutations (false negatives) were more common than erroneous candidates (false positives), indicating that the reliable sampling of rare ctDNA fragments is the key challenge for ctDNA assays. This comprehensive evaluation of the analytical performance of ctDNA assays serves to inform best practice guidelines and provides a resource for precision oncology.

Circulating Tumor DNA↗

Correction of scaling mismatches in oligonucleotide microarray data.

BACKGROUND: Gene expression microarray data is notoriously subject to high signal variability. Moreover, unavoidable variation in the concentration of transcripts applied to microarrays may result in poor scaling of the summarized data which can hamper analytical interpretations. This is especially relevant in a systems biology context, where systematic biases in the signals of particular genes can have severe effects on subsequent analyses. Conventionally it would be necessary to replace the mismatched arrays, but individual time points cannot be rerun and inserted because of experimental variability. It would therefore be necessary to repeat the whole time series experiment, which is both impractical and expensive. RESULTS: We explain how scaling mismatches occur in data summarized by the popular MAS5 (GCOS; Affymetrix) algorithm, and propose a simple recursive algorithm to correct them. Its principle is to identify a set of constant genes and to use this set to rescale the microarray signals. We study the properties of the algorithm using artificially generated data and apply it to experimental data. We show that the set of constant genes it generates can be used to rescale data from other experiments, provided that the underlying system is similar to the original. We also demonstrate, using a simple example, that the method can successfully correct existing imbalances in the data. CONCLUSION: The set of constant genes obtained for a given experiment can be applied to other experiments, provided the systems studied are sufficiently similar. This type of rescaling is especially relevant in systems biology applications using microarray data.

Algorithms↗

Ranked prediction of p53 targets using hidden variable dynamic modeling.

Full exploitation of microarray data requires hidden information that cannot be extracted using current analysis methodologies. We present a new approach, hidden variable dynamic modeling (HVDM), which derives the hidden profile of a transcription factor from time series microarray data, and generates a ranked list of predicted targets. We applied HVDM to the p53 network, validating predictions experimentally using small interfering RNA. HVDM can be applied in many systems biology contexts to predict regulation of gene activity quantitatively.

Cell Line, Tumor↗

Radiation-induced delayed cell death in a hypomorphic Artemis cell line.

Null mutations in Artemis confer a condition described as RS-SCID, in which patients display radiosensitivity combined with severe combined immunodeficiency. Here, we characterize the defect in Artemis in a patient who displayed progressive combined immunodeficiency (CID) and elevated lymphocyte apoptosis. The patient is a compound heterozygote with novel mutations in both alleles, resulting in Artemis proteins with either L70 deletion or G126D substitution. Both mutational changes impact upon Artemis function and a fibroblast cell line derived from the patient (F96-224) has greatly reduced Artemis protein. In contrast to Artemis null cell lines, which fail to repair a subset of DNA double strand breaks (DSBs) induced by ionizing radiation, F96-224 cells show slow but residual DSB rejoining. Despite showing intermediate cellular and clinical features, F96-224 cells are as radiosensitive as Artemis null cell lines. We developed a FACS-based assay to examine cell division and cellular characteristics for 10 days following exposure to ionizing radiation (2 and 4 Gy). This analysis demonstrated that F96-224 cells show delayed cell death when compared with rapid growth arrest of an Artemis null cell line, and the emergence of a cycling population shown by a control line. F96-224 cells also display elevated chromosome aberrations when compared with control cells. F96-224 therefore represents a novel phenotype for a hypomorphic cell line. We suggest that delayed cell death contributes to the progressive CID phenotype of the Artemis patient.

Cell Death↗

Robustness of the p53 network and biological hackers.

The p53 protein interaction network is crucial in regulating the metazoan cell cycle and apoptosis. Here, the robustness of the p53 network is studied by analyzing its degeneration under two modes of attack. Linear Programming is used to calculate average path lengths among proteins and the network diameter as measures of functionality. The p53 network is found to be robust to random loss of nodes, but vulnerable to a targeted attack against its hubs, as a result of its architecture. The significance of the results is considered with respect to mutational knockouts of proteins and the directed attacks mounted by tumour inducing viruses.

Computer Simulation↗

STAT-1 facilitates the ATM activated checkpoint pathway following DNA damage.

STAT-1 plays a role in mediating stress responses to various stimuli and has also been implied to be a tumour suppressor. Here, we report that STAT-1-deficient cells have defects both in intra-S-phase and G2-M checkpoints in response to DNA damage. Interestingly, STAT-1-deficient cells showed reduced Chk2 phosphorylation on threonine 68 (Chk2(-T68)) following DNA damage, suggesting that STAT-1 might function in the ATM-Chk2 pathway. Moreover, the defects in Chk2(-T68) phosphorylation in STAT-1-deficient cells also correlated with reduced degradation of Cdc25A compared with STAT-1-expressing cells after DNA damage. We also show that STAT-1 is required for ATM-dependent phosphorylation of NBS1 and p53 but not for BRCA1 or H2AX phosphorylation following DNA damage. Expression levels of BRCT mediator/adaptor proteins MDC1 and 53BP1, which are required for ATM-mediated pathways, are reduced in cells lacking STAT-1. Enforced expression of MDC1 into STAT-1-deficient cells restored ATM-mediated phosphorylation of downstream substrates. These results imply that STAT-1 plays a crucial role in the DNA-damage-response by regulating the expression of 53BP1 and MDC1, factors known to be important for mediating ATM-dependent checkpoint pathways.

Adaptor Proteins, Signal Transducing↗

Sperm-induced modification of the oviductal gene expression profile after natural insemination in mice.

In mammals, the physiological interaction between spermatozoa and oviductal epithelia involves intimate and specific contact between the two cell types. Spermatozoa may undergo stringent selection processes within the female reproductive tract before they meet and fertilize oocytes. The physiological basis of the sperm selection process is largely unknown. Here we tested the hypothesis that the oviduct has a recognition system for spermatozoa that can detect the arrival of spermatozoa in the oviduct after insemination, resulting in alterations of the oviductal transcriptome. We initially performed a global screening of the oviductal transcriptome in mice 1) at the time of estrus (mating) and 2) 6 h after mating. Transcriptional alterations in the oviduct after mating were attributed to the presence of spermatozoa in the oviduct after mating and also to changes in the hormonal environment as female mice underwent the transition from estrus to diestrus. To distinguish these possibilities, female mice were then mated with T145H mutant mice, which because of spermatogenic arrest, produce seminal plasma but no spermatozoa. Focusing on two molecules that in the first experiment were upregulated after mating, it was found that adrenomedullin and prostaglandin endoperoxidase synthase 2 transcripts were upregulated in the oviducts of mice only after mating with fertile males; those mated with T145H infertile males showed significantly less response. These results indicate that it is the arrival of spermatozoa in the oviduct that activates one or more signal transduction pathways and leads to changes in the oviductal transcriptome profiles.

Adrenomedullin↗

Cloning of apoptosis-related genes by representational difference analysis of cDNA.

Apoptosis is frequently triggered by events that alter the expression of key target genes. Under these circumstances, the genes involved can be identified by techniques that analyze gene expression. Researchers now have a choice of reliable and effective methods for differential gene expression analysis. Comparative approaches, including gene microarray analysis, serial analysis of gene expression, and differential display provide global information about expression levels. Subtractive approaches like complementary DNA representational difference analysis (cDNA RDA) and suppression subtraction polymerase chain reaction identify a focused set of differentially expressed genes. The most suitable technique to apply depends on individual circumstances. cDNA RDA is particularly useful in nonstandard model organisms for which comprehensive gene microarrays are not available and is best used for the identification of genes with a large difference in expression levels between two populations. The technique involves the generation of amplified mixtures of cDNA fragments that are typically smaller than 1000 base pairs and represent >86% of mRNA species from each starting population. Transcriptional differences between two populations can then be identified by subtraction of cDNA amplicons followed by further polymerase chain reaction amplification. The technique is capable of detecting differences for genes expressed at less than one copy per cell and is achievable using standard laboratory apparatus. cDNA RDA can identify genes not previously described in the database, can detect low abundance transcripts (e.g., from mixed cell populations), and is best applied in experiments where relatively few differentially expressed genes are expected. Here, we describe the application of cDNA RDA to the identification of apoptosis-related genes.

Animals↗

From the top down: towards a predictive biology of signalling networks.

True quantitative and predictive biology requires the ability to interpret increasingly complex datasets in new ways to reveal underlying functional interactions. This has come a step closer with two recent articles that describe a 'top-down' modelling approach to reconstructing functional networks from microarray data.

Computer Simulation↗

Discovery of genes with highly restricted expression patterns in the Drosophila wing disc using DNA oligonucleotide microarrays.

The Drosophila wing disc is divided along the proximal-distal axis into regions giving rise to the body wall (proximal), wing hinge (central) and wing blade (distal). We applied DNA microarray analysis to discover genes with potential roles in the development of these regions. We identified a set of 94 transcripts enriched (two fold or greater) in the body wall and 56 transcripts enriched in the wing/hinge region. Transcripts that are known to have highly restricted expression patterns, such as pannier, twist and Bar-H1 (body wall) and knot, nubbin and Distal-less (wing/hinge), showed strong differential expression on the arrays. In situ hybridization for 50 previously uncharacterized genes similarly revealed that transcript enrichment identified by the array analysis was consistent with the observed spatial expression. There was a broad spectrum of patterns, in some cases suggesting that the genes could be targets of known signaling pathways. We show that three of these genes respond to wingless signaling. We also discovered genes likely to play specific roles in tracheal and myoblast cell types, as these cells are part of the body wall fragment. In summary, the identification of genes with restricted expression patterns using whole genome profiling suggests that many genes with potential roles in wing disc development remain to be characterized.

Animals↗