Single-cell immune profiling reveals markers of emergency myelopoiesis that distinguish severe from mild respiratory syncytial virus disease in infants

Whereas most infants infected with respiratory syncytial virus (RSV) show no or only mild symptoms, an estimated 3 million children under five are hospitalized annually due to RSV disease. This study aimed to investigate biological mechanisms and associated biomarkers underlying RSV disease heterogeneity in young infants, enabling the potential to objectively categorize RSV-infected infants according to their medical needs. Immunophenotypic and functional profiling demonstrated the emergence of immature and progenitor-like neutrophils, proliferative monocytes (HLA-DRLow , Ki67+), impaired antigen-presenting function, downregulation of T cell response and low abundance of HLA-DRLow B cells in severe RSV disease. HLA-DRLow monocytes were found as a hallmark of RSV-infected infants requiring hospitalization. Complementary transcriptomics identified genes associated with disease severity and pointed to the emergency myelopoiesis response. These results shed new light on mechanisms underlying the pathogenesis and development of severe RSV disease and identified potential new candidate biomarkers for patient stratification.

Blocking NS3-NS4B interaction inhibits dengue virus in non-human primates

Dengue is a major health threat and the number of symptomatic infections caused by the four dengue serotypes is estimated to be 96 million1 with annually around 10,000 deaths2. However, no antiviral drugs are available for the treatment or prophylaxis of dengue. We recently described the interaction between non-structural proteins NS3 and NS4B as a promising target for the development of pan-serotype dengue virus (DENV) inhibitors3. Here we present JNJ-1802-a highly potent DENV inhibitor that blocks the NS3-NS4B interaction within the viral replication complex. JNJ-1802 exerts picomolar to low nanomolar in vitro antiviral activity, a high barrier to resistance and potent in vivo efficacy in mice against infection with any of the four DENV serotypes. Finally, we demonstrate that the small-molecule inhibitor JNJ-1802 is highly effective against viral infection with DENV-1 or DENV-2 in non-human primates. JNJ-1802 has successfully completed a phase I first-in-human clinical study in healthy volunteers and was found to be safe and well tolerated4. These findings support the further clinical development of JNJ-1802, a first-in-class antiviral agent against dengue, which is now progressing in clinical studies for the prevention and treatment of dengue.

Virology analysis in HCV genotype 1-infected patients treated with the combination of simeprevir and TMC647055/ritonavir, with and without ribavirin, and JNJ-56914845

BACKGROUND

In study TMC647055HPC2001, a 3-direct-acting-antiviral (DAA) regimen combining NS3/4A protease inhibitor simeprevir (SMV), non-nucleoside NS5B inhibitor TMC647055/ritonavir (RTV) and NS5A inhibitor JNJ-56914845 resulted in high sustained virologic response 12 weeks after actual end of treatment (SVR12) in chronic hepatitis C virus (HCV) genotype 1-infected patients. SVR12 rates were generally lower in the 2-DAA regimen SMV + TMC647055/RTV with or without ribavirin. The objective of this study was to identify and characterise pre-existing and emerging resistance-associated variants (RAVs) in patients enrolled in study TMC647055HPC2001.

METHODS

HCV population sequencing analyses were performed on baseline isolates from all patients (n = 90) and post-baseline isolates from patients with virologic failure (n = 22). In addition, deep sequencing and phenotypic analyses were performed on selected baseline and post-baseline isolates.

RESULTS

The majority of patients with virologic failure had emerging RAVs to all study drugs at the time of failure: in all 22 patients SMV RAVs emerged at NS3 positions 80, 155, 156 and/or 168, consistent with the known SMV resistance profile. Emerging TMC647055 RAVs at NS5B position 495 were detected in the majority of patients (16/22), and all 5 patients who failed the 3-DAA regimen had emerging JNJ-56914845 RAVs at NS5A positions 30 and/or 31. While at the end of study emerging SMV and TMC647055 RAVs were no longer observed by population sequencing in 40% (8/20) and 62.5% (10/16) of patients with follow-up data available, respectively, emerging JNJ-56914845 RAVs were still detected in all (5/5) patients.

CONCLUSIONS

Virologic failure in the 2- and 3-DAA combinations was, in the majority of patients, associated with the emergence of RAVs to all study drugs. While emerging SMV and TMC647055 RAVs became undetectable during follow-up, JNJ-56914845 RAVs in NS5A were still observed at end of study.

TRIAL REGISTRATION NUMBER

NCT01724086 (date of registration: September 26, 2012).

Living Organ Donation by Minors: An Analysis of the Regulations in European Union Member States

Living organ donation (LD) is an increasingly established practice. Whereas in the United States and Canada LD by minors has occasionally been reported, LD by minors seems to be largely absent in the European Union (EU). It is currently unclear whether this is the result of a different legal approach. This study is the first to systematically analyze the regulations of EU member states, Norway, and Iceland toward LD by minors. Relevant regulations were identified by searching government websites, translated, compared, and sent for verification to national legal experts. We identified five countries where LD by minors is allowed. In two of these (Belgium and the United Kingdom), some minors may be deemed sufficiently mature to make an autonomous decision regarding LD. In contrast, in the three other countries (Luxembourg, Norway, and Sweden), LD by minors is only allowed subject to parental permission and the assent (or absence of objection) of the donor. Where allowed, regulations differ significantly with regard to the substantive and procedural safeguards in place. In view of the controversial nature of the procedure, as illustrated by recent reports and surveys, we argue for a very cautious approach and greater harmonization in countries where LD by minors is allowed.

HCMV Displays a Unique Transcriptome of Immunomodulatory Genes in Primary Monocyte-Derived Cell Types

Human cytomegalovirus (HCMV) is a betaherpesvirus which rarely presents problems in healthy individuals, yet may result in severe morbidity in immunocompromised patients and in immune-naïve neonates. HCMV has a large 235 kb genome with a coding capacity of at least 165 open reading frames (ORFs). This large genome allows complex gene regulation resulting in different sets of transcripts during lytic and latent infection. While latent virus mainly resides within monocytes and CD34+ progenitor cells, reactivation to lytic infection is driven by differentiation towards terminally differentiated myeloid dendritic cells and macrophages. Consequently, it has been suggested that macrophages and dendritic cells contribute to viral spread in vivo. Thus far only limited knowledge is available on the expression of HCMV genes in terminally differentiated myeloid primary cells and whether or not the virus exhibits a different set of lytic genes in primary cells compared with lytic infection in NHDF fibroblasts. To address these questions, we used Illumina next generation sequencing to determine the HCMV transcriptome in macrophages and dendritic cells during lytic infection and compared it to the transcriptome in NHDF fibroblasts. Here, we demonstrate unique expression profiles in macrophages and dendritic cells which significantly differ from the transcriptome in fibroblasts mainly by modulating the expression of viral transcripts involved in immune modulation, cell tropism and viral spread. In a head to head comparison between macrophages and dendritic cells, we observed that factors involved in viral spread and virion composition are differentially regulated suggesting that the plasticity of the virion facilitates the infection of surrounding cells. Taken together, this study provides the full transcript expression analysis of lytic HCMV genes in monocyte-derived type 1 and type 2 macrophages as well as in monocyte-derived dendritic cells. Thereby underlining the potential of HCMV to adapt to or influence different cellular environments to promote its own survival.

Pre-existence and Persistence of Resistant Minority Hepatitis C Virus Variants in Genotype 1-Infected Patients Treated With Simeprevir/Peginterferon/Ribavirin

Deep sequencing analyses in HCV genotype 1-infected patients treated with simeprevir/peginterferon/ribavirin showed that pre-existing minority HCV variants did not impact treatment outcome and emerging resistant variants in patients failing treatment were not persisting at minority levels

Background. The pre-existence of minority hepatitis C virus (HCV) variants and their impact on treatment outcome, as well as the persistence of emerging resistant variants posttreatment in patients failing treatment with simeprevir/peginterferon/ribavirin (SMV/PR), were assessed by deep sequencing (DS).

Methods. Population sequencing (PS) and Illumina DS were performed on HCV genotype 1 isolates from patients treated with SMV/PR in Phase 2b (PILLAR [NCT00882908] and ASPIRE [NCT00980330]) and Phase 3 (QUEST-1 [NCT01289782], QUEST-2 [NCT01290679], and PROMISE [NCT01281839]) trials.

Results. Minority polymorphisms (ie, detected pretreatment by DS only) reducing SMV activity in vitro were uncommon (3.6%, 19 of 534 patients). These SMV-resistant minority polymorphisms were detected in similar proportions of patients achieving (3.7%) and not achieving (3.3%) sustained virologic response with SMV/PR and generally did not emerge as major variants at time of failure. SMV-resistant variants emerging at time of failure were no longer detected at end of study in 69.3% and 52.0% of the patients by PS and DS, respectively.

Conclusions. Minority polymorphisms did not impact outcome of SMV/PR treatment. The majority of emerging variants that became undetectable at end of study by PS were also undetectable by DS. These results suggest no added value of DS for clinical usage of SMV.

Assessment of etravirine resistance in HIV-1-infected paediatric patients using population and deep sequencing: final results of the PIANO study

BACKGROUND

We assessed etravirine resistance in treatment-experienced, HIV-1-infected children (n=41)/adolescents (n=60) who received twice-daily etravirine 5.2 mg/kg and a background regimen (boosted protease inhibitor plus nucleoside/nucleotide reverse transcriptase inhibitors, optional enfuvirtide/raltegravir) in a Phase II, open-label, multicentre trial (PIANO).

METHODS

In addition to phenotypes, viral genotypes were assessed by population and deep sequencing (PS and DS) in virological failures (VFs; baseline and end point) and responders (baseline). Minority resistance-associated mutations (RAMs) were defined as those with frequencies above 1% and not detected with PS.

RESULTS

By week 48, 41/101 (40.6%) patients experienced VF; 17/41 (41.5%) VFs and 22/54 (40.8%) responders had ≥1 baseline etravirine RAM by PS, mainly A98G, K101E, V106I and G190A. Baseline minority etravirine RAMs (n) were detected in 8/40 VFs (V90I [2], A98G [1], L100I [1], V106I [1], E138G [1] and Y181C [2]) and 5/38 responders (V90I [3], A98G [1], V106I [1] and E138G [1]). The most frequent emerging non-nucleoside reverse transcriptase inhibitor RAMs detected by PS (≥3 VFs; n) were the etravirine RAMs Y181C (8), V90I (3), L100I (3) and E138A (3). In 15 of 29 (51.7%) VFs with baseline DS/PS and end point PS data, ≥1 emerging etravirine RAM was detected by PS, which was not detected at baseline by DS in most cases (12/15 [80.0%]). In 10/26 (38.5%) VFs with baseline/end point DS data, ≥1 additional emerging minority etravirine RAM was detected.

CONCLUSIONS

Patterns of etravirine resistance in adults, adolescents and children experiencing VF are similar. The presence of minority etravirine RAMs at baseline was not consistently associated with treatment failure. ClinicalTrials.gov: NCT00665847.

QQ-SNV: single nucleotide variant detection at low frequency by comparing the quality quantiles

Background

Next generation sequencing enables studying heterogeneous populations of viral infections. When the sequencing is done at high coverage depth (“deep sequencing”), low frequency variants can be detected. Here we present QQ-SNV (http://sourceforge.net/projects/qqsnv), a logistic regression classifier model developed for the Illumina sequencing platforms that uses the quantiles of the quality scores, to distinguish true single nucleotide variants from sequencing errors based on the estimated SNV probability. To train the model, we created a dataset of an in silico mixture of five HIV-1 plasmids. Testing of our method in comparison to the existing methods LoFreq, ShoRAH, and V-Phaser 2 was performed on two HIV and four HCV plasmid mixture datasets and one influenza H1N1 clinical dataset.

Results

For default application of QQ-SNV, variants were called using a SNV probability cutoff of 0.5 (QQ-SNV_D). To improve the sensitivity we used a SNV probability cutoff of 0.0001 (QQ-SNV_HS). To also increase specificity, SNVs called were overruled when their frequency was below the 80^th percentile calculated on the distribution of error frequencies (QQ-SNV_HS-P80).

When comparing QQ-SNV versus the other methods on the plasmid mixture test sets, QQ-SNV_D performed similarly to the existing approaches. QQ-SNV_HS was more sensitive on all test sets but with more false positives. QQ-SNV_HS-P80 was found to be the most accurate method over all test sets by balancing sensitivity and specificity. When applied to a paired-end HCV sequencing study, with lowest spiked-in true frequency of 0.5 %, QQ-SNV_HS-P80 revealed a sensitivity of 100 % (vs. 40–60 % for the existing methods) and a specificity of 100 % (vs. 98.0–99.7 % for the existing methods). In addition, QQ-SNV required the least overall computation time to process the test sets.

Finally, when testing on a clinical sample, four putative true variants with frequency below 0.5 % were consistently detected by QQ-SNV_HS-P80 from different generations of Illumina sequencers.

Conclusions

We developed and successfully evaluated a novel method, called QQ-SNV, for highly efficient single nucleotide variant calling on Illumina deep sequencing virology data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0812-9) contains supplementary material, which is available to authorized users.

Deep sequencing analysis of HIV-1 reverse transcriptase at baseline and time of failure in patients receiving rilpivirine in the phase III studies ECHO and THRIVE

Minority variants (1.0-25.0%) were evaluated by deep sequencing (DS) at baseline and virological failure (VF) in a selection of antiretroviral treatment-naïve, HIV-1-infected patients from the rilpivirine ECHO/THRIVE phase III studies. Linkage between frequently emerging resistance-associated mutations (RAMs) was determined. DS (llIumina®) and population sequencing (PS) results were available at baseline for 47 VFs and time of failure for 48 VFs; and at baseline for 49 responders matched for baseline characteristics. Minority mutations were accurately detected at frequencies down to 1.2% of the HIV-1 quasispecies. No baseline minority rilpivirine RAMs were detected in VFs; one responder carried 1.9% F227C. Baseline minority mutations associated with resistance to other non-nucleoside reverse transcriptase inhibitors (NNRTIs) were detected in 8/47 VFs (17.0%) and 7/49 responders (14.3%). Baseline minority nucleoside/nucleotide reverse transcriptase inhibitor (NRTI) RAMs M184V and L210W were each detected in one VF (none in responders). At failure, two patients without NNRTI RAMs by PS carried minority rilpivirine RAMs K101E and/or E138K; and five additional patients carried other minority NNRTI RAMs V90I, V106I, V179I, V189I, and Y188H. Overall at failure, minority NNRTI RAMs and NRTI RAMs were found in 29/48 (60.4%) and 16/48 VFs (33.3%), respectively. Linkage analysis showed that E138K and K101E were usually not observed on the same viral genome. In conclusion, baseline minority rilpivirine RAMs and other NNRTI/NRTI RAMs were uncommon in the rilpivirine arm of the ECHO and THRIVE studies. DS at failure showed emerging NNRTI resistant minority variants in seven rilpivirine VFs who had no detectable NNRTI RAMs by PS.

Performance assessment of the Illumina massively parallel sequencing platform for deep sequencing analysis of viral minority variants

Massively parallel sequencing (MPS) technology has opened new avenues to study viral dynamics and treatment-induced resistance mechanisms of infections such as human immunodeficiency virus (HIV) and hepatitis C virus (HCV). Whereas the Roche/454 platform has been used widely for the detection of low-frequent drug resistant variants, more recently developed short-read MPS technologies have the advantage of delivering a higher sequencing depth at a lower cost per sequenced base. This study assesses the performance characteristics of Illumina MPS technology for the characterization of genetic variability in viral populations by deep sequencing. The reported results from MPS experiments comprising HIV and HCV plasmids demonstrate that a 0.5-1% lower limit of detection can be achieved readily with Illumina MPS while retaining good accuracy also at low frequencies. Deep sequencing of a set of clinical samples (12 HIV and 9 HCV patients), designed at a similar budget for both MPS platforms, reveals a comparable lower limit of detection for Illumina and Roche/454. Finally, this study shows the possibility to apply Illumina's paired-end sequencing as a strategy to assess linkage between different mutations identified in individual viral subspecies. These results support the use of Illumina as another MPS platform of choice for deep sequencing of viral minority species.

ViVaMBC: estimating viral sequence variation in complex populations from illumina deep-sequencing data using model-based clustering

Background

Deep-sequencing allows for an in-depth characterization of sequence variation in complex populations. However, technology associated errors may impede a powerful assessment of low-frequency mutations. Fortunately, base calls are complemented with quality scores which are derived from a quadruplet of intensities, one channel for each nucleotide type for Illumina sequencing. The highest intensity of the four channels determines the base that is called. Mismatch bases can often be corrected by the second best base, i.e. the base with the second highest intensity in the quadruplet. A virus variant model-based clustering method, ViVaMBC, is presented that explores quality scores and second best base calls for identifying and quantifying viral variants. ViVaMBC is optimized to call variants at the codon level (nucleotide triplets) which enables immediate biological interpretation of the variants with respect to their antiviral drug responses.

Results

Using mixtures of HCV plasmids we show that our method accurately estimates frequencies down to 0.5%. The estimates are unbiased when average coverages of 25,000 are reached. A comparison with the SNP-callers V-Phaser2, ShoRAH, and LoFreq shows that ViVaMBC has a superb sensitivity and specificity for variants with frequencies above 0.4%. Unlike the competitors, ViVaMBC reports a higher number of false-positive findings with frequencies below 0.4% which might partially originate from picking up artificial variants introduced by errors in the sample and library preparation step.

Conclusions

ViVaMBC is the first method to call viral variants directly at the codon level. The strength of the approach lies in modeling the error probabilities based on the quality scores. Although the use of second best base calls appeared very promising in our data exploration phase, their utility was limited. They provided a slight increase in sensitivity, which however does not warrant the additional computational cost of running the offline base caller. Apparently a lot of information is already contained in the quality scores enabling the model based clustering procedure to adjust the majority of the sequencing errors. Overall the sensitivity of ViVaMBC is such that technical constraints like PCR errors start to form the bottleneck for low frequency variant detection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0458-7) contains supplementary material, which is available to authorized users.

VirVarSeq: a low-frequency virus variant detection pipeline for Illumina sequencing using adaptive base-calling accuracy filtering

MOTIVATION

In virology, massively parallel sequencing (MPS) opens many opportunities for studying viral quasi-species, e.g. in HIV-1- and HCV-infected patients. This is essential for understanding pathways to resistance, which can substantially improve treatment. Although MPS platforms allow in-depth characterization of sequence variation, their measurements still involve substantial technical noise. For Illumina sequencing, single base substitutions are the main error source and impede powerful assessment of low-frequency mutations. Fortunately, base calls are complemented with quality scores (Qs) that are useful for differentiating errors from the real low-frequency mutations.

RESULTS

A variant calling tool, Q-cpileup, is proposed, which exploits the Qs of nucleotides in a filtering strategy to increase specificity. The tool is imbedded in an open-source pipeline, VirVarSeq, which allows variant calling starting from fastq files. Using both plasmid mixtures and clinical samples, we show that Q-cpileup is able to reduce the number of false-positive findings. The filtering strategy is adaptive and provides an optimized threshold for individual samples in each sequencing run. Additionally, linkage information is kept between single-nucleotide polymorphisms as variants are called at the codon level. This enables virologists to have an immediate biological interpretation of the reported variants with respect to their antiviral drug responses. A comparison with existing SNP caller tools reveals that calling variants at the codon level with Q-cpileup results in an outstanding sensitivity while maintaining a good specificity for variants with frequencies down to 0.5%.

AVAILABILITY

The VirVarSeq is available, together with a user's guide and test data, at sourceforge: http://sourceforge.net/projects/virtools/?source=directory.

Acquired resistance of Mycobacterium tuberculosis to bedaquiline

Bedaquiline (BDQ), an ATP synthase inhibitor, is the first drug to be approved for treatment of multi-drug resistant tuberculosis in decades. In vitro resistance to BDQ was previously shown to be due to target-based mutations. Here we report that non-target based resistance to BDQ, and cross-resistance to clofazimine (CFZ), is due to mutations in Rv0678, a transcriptional repressor of the genes encoding the MmpS5-MmpL5 efflux pump. Efflux-based resistance was identified in paired isolates from patients treated with BDQ, as well as in mice, in which it was confirmed to decrease bactericidal efficacy. The efflux inhibitors verapamil and reserpine decreased the minimum inhibitory concentrations of BDQ and CFZ in vitro, but verapamil failed to increase the bactericidal effect of BDQ in mice and was unable to reverse efflux-based resistance in vivo. Cross-resistance between BDQ and CFZ may have important clinical implications.

Deep-sequencing analysis of the gene encoding the hepatitis C virus nonstructural 3-4A protease confirms a low prevalence of telaprevir-resistant variants at baseline and the end of the REALIZE study

BACKGROUND

Population sequencing (PS) has shown that telaprevir-resistant variants are not typically detectable at baseline (prevalence, ≤5% of patients), and most variants present at the time of treatment failure are no longer detectable at the end of the study.

METHODS

To gain insight into the evolution of telaprevir-resistant variants, their baseline prevalence and persistence after treatment was investigated using a more sensitive, deep-sequencing (DS) technique in a large number of treatment-experienced patients from the REALIZE study who were infected with hepatitis C virus genotype 1.

RESULTS

Before treatment initiation, telaprevir-resistant variants (T54A, T54S, or R155K in 1%-2% of the viral population) were detected by DS in a fraction (2%) of patients for whom PS failed to detect resistance; these variants were not necessarily detected at the time of treatment failure. Of 49 patients in whom telaprevir-resistant variants were detected by PS at the time of treatment failure but not at the end of the study, DS revealed the presence of variants (V36A/L/M, T54S, or R155K in 1%-36% of the viral population) in 16 patients (33%) at the end of the study.

CONCLUSIONS

Similar to PS findings, DS analysis revealed that the frequency of telaprevir-resistant variants before treatment was also low, and variants detected at the time of treatment failure were no longer detectable in the majority of patients during follow-up.

A method enabling high-throughput sequencing of human cytomegalovirus complete genomes from clinical isolates

Human cytomegalovirus (HCMV) is a ubiquitous virus that can cause serious sequelae in immunocompromised patients and in the developing fetus. The coding capacity of the 235 kbp genome is still incompletely understood, and there is a pressing need to characterize genomic contents in clinical isolates. In this study, a procedure for the high-throughput generation of full genome consensus sequences from clinical HCMV isolates is presented. This method relies on low number passaging of clinical isolates on human fibroblasts, followed by digestion of cellular DNA and purification of viral DNA. After multiple displacement amplification, highly pure viral DNA is generated. These extracts are suitable for high-throughput next-generation sequencing and assembly of consensus sequences. Throughout a series of validation experiments, we showed that the workflow reproducibly generated consensus sequences representative for the virus population present in the original clinical material. Additionally, the performance of 454 GS FLX and/or Illumina Genome Analyzer datasets in consensus sequence deduction was evaluated. Based on assembly performance data, the Illumina Genome Analyzer was the platform of choice in the presented workflow. Analysis of the consensus sequences derived in this study confirmed the presence of gene-disrupting mutations in clinical HCMV isolates independent from in vitro passaging. These mutations were identified in genes RL5A, UL1, UL9, UL111A and UL150. In conclusion, the presented workflow provides opportunities for high-throughput characterization of complete HCMV genomes that could deliver new insights into HCMV coding capacity and genetic determinants of viral tropism and pathogenicity.

Efficiency of whole genome amplification of single circulating tumor cells enriched by CellSearch and sorted by FACS

BACKGROUND

Tumor cells in the blood of patients with metastatic carcinomas are associated with poor survival. Knowledge of the cells' genetic make-up can help to guide targeted therapy. We evaluated the efficiency and quality of isolation and amplification of DNA from single circulating tumor cells (CTC).

METHODS

The efficiency of the procedure was determined by spiking blood with SKBR-3 cells, enrichment with the CellSearch system, followed by single cell sorting by fluorescence-activated cell sorting (FACS) and whole genome amplification. A selection of single cell DNA from fixed and unfixed SKBR-3 cells was exome sequenced and the DNA quality analyzed. Single CTC from patients with lung cancer were used to demonstrate the potential of single CTC molecular characterization.

RESULTS

The overall efficiency of the procedure from spiked cell to amplified DNA was approximately 20%. Losses attributed to the CellSearch system were around 20%, transfer to FACS around 25%, sorting around 5% and DNA amplification around 25%. Exome sequencing revealed that the quality of the DNA was affected by the fixation of the cells, amplification, and the low starting quantity of DNA. A single fixed cell had an average coverage at 20× depth of 30% when sequencing to an average of 40× depth, whereas a single unfixed cell had 45% coverage. GenomiPhi-amplified genomic DNA had a coverage of 72% versus a coverage of 87% of genomic DNA. Twenty-one percent of the CTC from patients with lung cancer identified by the CellSearch system could be isolated individually and amplified.

CONCLUSIONS

CTC enriched by the CellSearch system were sorted by FACS, and DNA retrieved and amplified with an overall efficiency of 20%. Analysis of the sequencing data showed that this DNA could be used for variant calling, but not for quantitative measurements such as copy number detection. Close to 55% of the exome of single SKBR-3 cells were successfully sequenced to 20× depth making it possible to call 72% of the variants. The overall coverage was reduced to 30% at 20× depth, making it possible to call 56% of the variants in CellSave-fixed cells.

Ultra-deep sequencing of HIV-1 reverse transcriptase before start of an NNRTI-based regimen in treatment-naive patients

There are conflicting data on the impact of low frequency HIV-1 drug-resistant mutants on the response of first-line highly active antiretroviral therapy (HAART), more specifically containing a NNRTI. As population sequencing does not detect resistant viruses representing less than 15-25% of the viral population, more sensitive techniques have been developed but still need clinical validation. We evaluated ultra-deep sequencing (UDPS), recently more available and affordable, as a tool for the detection of HIV-1 minority species carrying drug resistant mutation (DRM) in a clinical setting. A retrospective analysis of the reverse transcriptase (RT) gene of plasma HIV-1 from 70 patients starting a NNRTI based regimen was performed. Minority populations were defined as representing > 1% and < 20% of the total viral population. Using UDPS, we could not confirm an association between the presence of low minority variants harbouring RT mutations at the start of therapy and primary or secondary therapeutic failure.

HIV-1 nucleotide mixture detection in the virco(®)TYPE HIV-1 genotyping assay: a comparison between Sanger sequencing and 454 pyrosequencing

HIV-1 Protease (PR) and Reverse Transcriptase (RT) genotyping is well established for the management of antiretroviral (ARV) drug therapy, as it is able to detect gene mutations encoding resistance to ARV compounds or drug classes, that are associated with reduced drug susceptibility (i.e. phenotype). A correct phenotypic interpretation from the derived PR-RT genotype (i.e. virtual phenotype), requires a well characterized geno-phenotype correlative database and appropriate statistical predictive models. The applicability of the virtual phenotype for the patient, will, however, not only depend on the accuracy of the statistical models and the database they rely on, but also depend largely on the sequence information that is provided. Since HIV-1 evolves as a complex of closely related but non-identical viral genomes (i.e. quasispecies) it is crucial that the sequencing method used, is able to characterize most of the genetic mixtures that make up the different quasispecies within a single patient. US regulatory agencies require that developers of HIV-1 genotyping assays, determine and report the HIV-1 mixture detection level of their assay. Hence, the mixture scoring sensitivity of the population-based Sanger sequencing method, along with the defined mixture scoring rules, used to drive the virco(®)TYPE HIV-1 virtual phenotype, was investigated by comparing it to the 454 pyrosequencing technique, which is able to generate the complete viral population sequence. To this end the PR-RT coding sequence of 20 clinical isolates was determined by both sequencing methodologies. The genotyping assay which feeds the virco(®)TYPE HIV-1 virtual phenotype was able to call automatically 97.5% (i.e. 268 mixtures) and 95.3% (i.e. 326 mixtures) of the mixtures that were present between 25 and 75% and between 20 and 80% in the viral population, as detected by 454. From the not called mixtures, all but one did present a mixture sequence in the Sanger DNA chromatograms, however, with a peak surface area for the second peak that was below the threshold setting for automatic mixture calling in the basecaller software (i.e. 25%). Viral loads ranged from 470 to 629,000 copies/mL and exerted no effect on the mixture calling relationship between both sequencing methodologies (R(2)=0.92). In some occasions (i.e. 55 mixtures) the genotyping assay would detect automatically mixtures that were present below 20% in the viral population, when measured by 454. Hence, the mixture scoring sensitivity of the automated high throughput virco(®)TYPE HIV-1 genotyping assay is currently set at 97.5% and 95.3%, for mixtures present at 25 and 20% in the viral population and may identify occasionally mutations that are present at lower frequencies. These findings were not influenced by the viral load of the examined samples.

HIV-1 V3 envelope deep sequencing for clinical plasma specimens failing in phenotypic tropism assays

Background

HIV-1 infected patients for whom standard gp160 phenotypic tropism testing failed are currently excluded from co-receptor antagonist treatment. To provide patients with maximal treatment options, massively parallel sequencing of the envelope V3 domain, in combination with tropism prediction tools, was evaluated as an alternative tropism determination strategy. Plasma samples from twelve HIV-1 infected individuals with failing phenotyping results were available. The samples were submitted to massive parallel sequencing and to confirmatory recombinant phenotyping using a fraction of the gp120 domain.

Results

A cut-off for sequence reads interpretation of 5 to10 times the sequencing error rate (0.2%) was implemented. On average, each sample contained 7 different V3 haplotypes. V3 haplotypes were submitted to tropism prediction algorithms, and 4/14 samples returned with presence of a dual/mixed (D/M) tropic virus, respectively at 3%, 10%, 11%, and 95% of the viral quasispecies. V3 tropism prediction was confirmed by gp120 phenotyping, except for two out of 4 D/M predicted viruses (with 3 and 95%) which were phenotypically R5-tropic. In the first case, the result was discordant due to the limit of detection for the phenotyping technology, while in the latter case the prediction algorithms were not computing the viral tropism correctly.

Conclusions

Although only demonstrated on a limited set of samples, the potential of the combined use of "deep sequencing + prediction algorithms" in cases where routine gp160 phenotype testing cannot be employed was illustrated. While good concordance was observed between gp120 phenotyping and prediction of R5-tropic virus, the results suggest that accurate prediction of X4-tropic virus would require further algorithm development.

Adenoviral vectors expressing siRNAs for discovery and validation of gene function

RNA interference is a powerful tool for studying gene function and for drug target discovery in diverse organisms and cell types. In mammalian systems, small interfering RNAs (siRNAs), or DNA plasmids expressing these siRNAs, have been used to down-modulate gene expression. However, inefficient transfection protocols, in particular, for primary cell types, have hampered the use of these tools in disease-relevant cellular assays. To be able to use this technology for genome-wide function screening, a more robust transduction protocol, resulting in a longer duration of the knock-down effect, is required. Here, we describe the validation of adenoviral vectors that express hairpin RNAs that are further processed to siRNAs. Infection of cell lines, or primary human cells, with these viruses leads to an efficient, sequence-specific, and prolonged reduction of the corresponding target mRNA, resulting in a reduction of the encoded protein level in the cell. For knock-down of one of the targets, GalphaS, we have measured inhibition of ligand-dependent, G-protein-coupled signaling. It is expected that this technology will prove to be of great value in target validation and target discovery efforts.

Homeobox gene expression profile in human hematopoietic multipotent stem cells and T-cell progenitors: implications for human T-cell development

Class I homeobox (HOX) genes comprise a large family of transcription factors that have been implicated in normal and malignant hematopoiesis. However, data on their expression or function during T-cell development is limited. Using degenerated RT-PCR and Affymetrix microarray analysis, we analyzed the expression pattern of this gene family in human multipotent stem cells from fetal liver (FL) and adult bone marrow (ABM), and in T-cell progenitors from child thymus. We show that FL and ABM stem cells are similar in terms of HOX gene expression, but significant differences were observed between these two cell types and child thymocytes. As the most immature thymocytes are derived from immigrated FL and ABM stem cells, this indicates a drastic change in HOX gene expression upon entry into the thymus. Further analysis of HOX-A7, HOX-A9, HOX-A10, and HOX-A11 expression with specific RT-PCR in all thymocyte differentiation stages showed a sequential loss of 3' region HOX-A cluster genes during intrathymic T-cell development and an unexpected expression of HOX-A11, previously not recognized to play a role in hematopoiesis. Also HOX-B3 and HOX-C4 were expressed throughout thymocyte development. Overall, these data provide novel evidence for an important role of certain HOX genes in human T-cell development.