Comparison of the performance of two targeted metagenomic virus capture probe-based methods using reference control materials and clinical samples.

Viral enrichment by probe hybridization has been reported to significantly increase the sensitivity of viral metagenomics. This study compares the analytical performance of two targeted metagenomic virus capture probe-based methods: (i) SeqCap EZ HyperCap by Roche (ViroCap) and (ii) Twist Comprehensive Viral Research Panel workflow, for diagnostic use. Sensitivity, specificity, and limit of detection were analyzed using 25 synthetic viral sequences spiked in increasing proportions of human background DNA, eight clinical samples, and American Type Culture Collection (ATCC) Virome Virus Mix. Sensitivity and specificity were 95% and higher for both methods using the synthetic and reference controls as gold standard. Combining thresholds for viral sequence read counts and genome coverage [respectively 500 reads per million (RPM) and 10% coverage] resulted in optimal prediction of true positive results. Limits of detection were approximately 50-500 copies/mL for both methods as determined by ddPCR. Increasing proportions of spike-in cell-free human background sequences up to 99.999% (50 ng/mL) did not negatively affect viral detection, suggesting effective capture of viral sequences. These data show analytical performances in ranges applicable to clinical samples, for both probe hybridization metagenomic approaches. This study supports further steps toward more widespread use of viral metagenomics for pathogen detection, in clinical and surveillance settings using low biomass samples. IMPORTANCE: Viral metagenomics has been gradually applied for broad-spectrum pathogen detection of infectious diseases, surveillance of emerging diseases, and pathogen discovery. Viral enrichment by probe hybridization methods has been reported to significantly increase the sensitivity of viral metagenomics. During the past years, a specific hybridization panel distributed by Roche has been adopted in a broad range of different clinical and zoonotic settings. Recently, Twist Bioscience has released a new hybridization panel targeting human and animal viruses. This is the first report comparing the performance of viral metagenomic hybridization panels.

A comparison of five Illumina, Ion Torrent, and nanopore sequencing technology-based approaches for whole genome sequencing of SARS-CoV-2.

Rapid identification of the rise and spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern remains critical for monitoring of the efficacy of diagnostics, therapeutics, vaccines, and control strategies. A wide range of SARS-CoV-2 next-generation sequencing (NGS) methods have been developed over the last years, but cross-sequence technology benchmarking studies have been scarce. In the current study, 26 clinical samples were sequenced using five protocols: AmpliSeq SARS-CoV-2 (Illumina), EasySeq RC-PCR SARS-CoV-2 (Illumina/NimaGen), Ion AmpliSeq SARS-CoV-2 (Thermo Fisher), custom primer sets (Oxford Nanopore Technologies (ONT)), and capture probe-based viral metagenomics (Roche/Illumina). Studied parameters included genome coverage, depth of coverage, amplicon distribution, and variant calling. The median SARS-CoV-2 genome coverage of samples with cycle threshold (Ct) values of 30 and lower ranged from 81.6 to 99.8% for, respectively, the ONT protocol and Illumina AmpliSeq protocol. Correlation of coverage with PCR Ct values varied per protocol. Amplicon distribution signatures differed across the methods, with peak differences of up to 4 log10 at disbalanced positions in samples with high viral loads (Ct values ≤ 23). Phylogenetic analyses of consensus sequences showed clustering independent of the workflow used. The proportion of SARS-CoV-2 reads in relation to background sequences, as a (cost-)efficiency metric, was the highest for the EasySeq protocol. The hands-on time was the lowest when using EasySeq and ONT protocols, with the latter additionally having the shortest sequence runtime. In conclusion, the studied protocols differed on a variety of the studied metrics. This study provides data that assist laboratories when selecting protocols for their specific setting.

Fusion transcripts and their genomic breakpoints in polyadenylated and ribosomal RNA-minus RNA sequencing data.

BACKGROUND: Fusion genes are typically identified by RNA sequencing (RNA-seq) without elucidating the causal genomic breakpoints. However, non-poly(A)-enriched RNA-seq contains large proportions of intronic reads that also span genomic breakpoints. RESULTS: We have developed an algorithm, Dr. Disco, that searches for fusion transcripts by taking an entire reference genome into account as search space. This includes exons but also introns, intergenic regions, and sequences that do not meet splice junction motifs. Using 1,275 RNA-seq samples, we investigated to what extent genomic breakpoints can be extracted from RNA-seq data and their implications regarding poly(A)-enriched and ribosomal RNA-minus RNA-seq data. Comparison with whole-genome sequencing data revealed that most genomic breakpoints are not, or minimally, transcribed while, in contrast, the genomic breakpoints of all 32 TMPRSS2-ERG-positive tumours were present at RNA level. We also revealed tumours in which the ERG breakpoint was located before ERG, which co-existed with additional deletions and messenger RNA that incorporated intergenic cryptic exons. In breast cancer we identified rearrangement hot spots near CCND1 and in glioma near CDK4 and MDM2 and could directly associate this with increased expression. Furthermore, in all datasets we find fusions to intergenic regions, often spanning multiple cryptic exons that potentially encode neo-antigens. Thus, fusion transcripts other than classical gene-to-gene fusions are prominently present and can be identified using RNA-seq. CONCLUSION: By using the full potential of non-poly(A)-enriched RNA-seq data, sophisticated analysis can reliably identify expressed genomic breakpoints and their transcriptional effects.

Profiling of promoter occupancy by PPARalpha in human hepatoma cells via ChIP-chip analysis.

The transcription factor peroxisome proliferator-activated receptor alpha (PPARalpha) is an important regulator of hepatic lipid metabolism. While PPARalpha is known to activate transcription of numerous genes, no comprehensive picture of PPARalpha binding to endogenous genes has yet been reported. To fill this gap, we performed Chromatin immunoprecipitation (ChIP)-chip in combination with transcriptional profiling on HepG2 human hepatoma cells treated with the PPARalpha agonist GW7647. We found that GW7647 increased PPARalpha binding to 4220 binding regions. GW7647-induced binding regions showed a bias around the transcription start site and most contained a predicted PPAR binding motif. Several genes known to be regulated by PPARalpha, such as ACOX1, SULT2A1, ACADL, CD36, IGFBP1 and G0S2, showed GW7647-induced PPARalpha binding to their promoter. A GW7647-induced PPARalpha-binding region was also assigned to SREBP-targets HMGCS1, HMGCR, FDFT1, SC4MOL, and LPIN1, expression of which was induced by GW7647, suggesting cross-talk between PPARalpha and SREBP signaling. Our data furthermore demonstrate interaction between PPARalpha and STAT transcription factors in PPARalpha-mediated transcriptional repression, and suggest interaction between PPARalpha and TBP, and PPARalpha and C/EBPalpha in PPARalpha-mediated transcriptional activation. Overall, our analysis leads to important new insights into the mechanisms and impact of transcriptional regulation by PPARalpha in human liver and highlight the importance of cross-talk with other transcription factors.

Tip-1 induces filopodia growth and is important for gastrulation movements during zebrafish development.

Wnt signaling is essential during animal development and also plays important roles in pathological conditions. Two mayor pathways have been described: the beta-catenin-dependent canonical (or classical) pathway and the beta-catenin-independent non-canonical Wnt pathway. Recent binding studies suggest links between the small PDZ protein TIP-1 (Tax-1 interacting protein) to components of both Wnt pathways. We have cloned and characterized the zebrafish tip-1 gene. Whole mount in situ hybridization and semiquantitative reverse transcriptase-polymerase chain reaction (RT-PCR) indicated that zebrafish tip-1 is present as a maternal RNA and is ubiquitously expressed during early development. After 24 h of development, tip-1 expression was high in the central nervous system (CNS) whereas only weak expression was detected in the caudal regions of the zebrafish embryo. Tip-1 knockdown using antisense morpholino oligonucleotides, as well as ectopic tip-1 expression, led to elongation defects in zebrafish embryos and larvae. Both knockdown and overexpression of tip-1 resulted in a widened goosecoid (gsc) expression domain in shield stage embryos, led to an abbreviated prechordal plate, and to reduced convergent extension movements during gastrulation. We constructed a green fluorescence protein (GFP)/TIP-1 fusion protein which, when expressed in cultured fibroblasts (ZF4-cells), induced filopodia growth. Our observations indicate a role for TIP-1 in gastrulation movements and in filopodia growth induction.

Zebrafish cypher is important for somite formation and heart development.

Mammalian CYPHER (Oracle, KIA0613), a member of the PDZ-LIM family of proteins (Enigma/LMP-1, ENH, ZASP/Cypher, RIL, ALP, and CLP-36), has been associated with cardiac and muscular myopathies. Targeted deletion of Cypher in mice is neonatal lethal possibly caused by myopathies. To further investigate the role of cypher in development, we have cloned the zebrafish orthologue. We present here the gene, domain structure, and expression pattern of zebrafish cypher during development. Cypher was not present as a maternal mRNA and was absent during early development. Cypher mRNA was first detected at the 3-somite stage in adaxial somites, and as somites matured, cypher expression gradually enveloped the whole somite. Later, cypher expression was also found in the heart, in head and jaw musculature, and in the brain. We further identified 13 alternative spliced forms of cypher from zebrafish heart and skeletal muscle tissue, among them a very short form containing the PDZ domain but lacking the ZM (ZASP-like) motif and the LIM domains. Targeted gene knock-down experiments using cypher antisense morpholinos led to severe defects, including truncation of the embryo, deformation of somites, dilatation of the pericardium, and thinning of the ventricular wall. The phenotype could be rescued by a cypher form, which contains the PDZ domain and the ZM motif, but lacks all three LIM domains. These findings indicate that a PDZ domain protein is important for normal somite formation and in normal heart development. Treatment of zebrafish embryos with cyclopamine, which disrupts hedgehog signaling, abolished cypher expression in 9 somite and 15-somite stage embryos. Taken together, our data suggest that cypher may play a role downstream of sonic hedgehog, in a late stage of somite development, when slow muscle fibers differentiate and migrate from the adaxial cells.

Gene expression profiling of the long-term adaptive response to hypoxia in the gills of adult zebrafish.

Low oxygen levels (hypoxia) play a role in clinical conditions such as stroke, chronic ischemia, and cancer. To better understand these diseases, it is crucial to study the responses of vertebrates to hypoxia. Among vertebrates, some teleosts have developed the ability to adapt to extremely low oxygen levels. We have studied long-term adaptive responses to hypoxia in adult zebrafish. We used zebrafish that survived severe hypoxic conditions for 3 wk and showed adaptive behavioral and phenotypic changes. We used cDNA microarrays to investigate hypoxia-induced changes in expression of 15,532 genes in the respiratory organs (the gills). We have identified 367 differentially expressed genes of which 117 showed hypoxia-induced and 250 hypoxia-reduced expressions. Metabolic depression was indicated by repression of genes in the TCA cycle in the electron transport chain and of genes involved in protein biosynthesis. We observed enhanced expression of the monocarboxylate transporter and of the oxygen transporter myoglobin. The hypoxia-induced group further included the genes for Niemann-Pick C disease and for Wolman disease [lysosomal acid lipase (LAL)]. Both diseases lead to a similar intra- and extracellular accumulation of cholesterol and glycolipids. The Niemann-Pick C protein binds to cholesterol from internal lysosomal membranes and is involved in cholesterol trafficking. LAL is responsible for lysosomal cholesterol degradation. Our data suggest a novel adaptive mechanism to hypoxia, the induction of genes for lysosomal lipid trafficking and degradation. Studying physiological responses to hypoxia in species tolerant for extremely low oxygen levels can help identify novel regulatory genes, which may have important clinical implications.