Search

GLI1 encodes a transcription factor that targets cell cycle regulators affecting stem cell proliferation. GLI1 gene fusions were initially described in pericytomas with a t[7;12] translocation and more recently in gastric plexiform fibromyxomas and gastroblastomas. This study describes the clinicopathologic, immunohistochemical, and molecular features of three intestinal-based neoplasms harboring GLI1 gene fusions. We studied three unique mesenchymal small bowel tumors. Paraffin embedded tumor tissues from these cases and 62 additional tumor samples that included a plexiform fibromyxoma were sequenced using a targeted RNAseq method to detect fusion events. The study patients included two women and one man who were 52, 80, and 22 years of age at the time of diagnosis. The tumors involved the submucosa and muscularis propria of the duodenum, jejunum, and ileum. All 3 tumors contained a proliferation of monotonous oval or spindle cells with scattered, somewhat dilated vessels. Two cases showed epithelioid structures such as glands, tubules, or nests. Immunohistochemical analysis revealed cytokeratin expression in the epithelioid components of both tumors displaying these features, and variable numbers of mesenchymal cells. Diffuse CD56 positivity was seen in the mesenchymal component of 2 tumors and desmin and smooth muscle actin staining in the other tumor. Immunostains for S-100 protein, DOG-1, and CD117 were negative in all cases. GLI1 fusions with different partner genes were detected in all tumors, and in the plexiform fibromyxoma, used as a control. Validation by fluorescence in situ hybridization was performed. None of the tumors have recurred or metastasize after surgery. We describe novel GLI1 fusions in 3 mesenchymal neoplasms of the small intestine, including 2 with biphenotypic features. Thus far, all cases have pursued indolent clinical courses. We propose the term " GLI1 -rearranged enteric tumor" to encompass this group of unique neoplasms of the small intestine that harbor GLI1 gene fusions and expand the spectrum of gastrointestinal neoplasms with these alterations.

Salivary gland neoplasms may pose diagnostic difficulties due to overlapping morphologic features. Recently, specific gene fusions have been discovered that correspond to particular tumor types, and can aid in accurate diagnosis. Gene rearrangements are commonly assessed by fluorescence in situ hybridization (FISH), although use of next-generation sequencing is increasing. However, there is no "gold standard" for fusion detection. We determined the concordance between FISH and a targeted RNA sequencing panel in gene fusion detection across twenty-two salivary gland tumors, including five mucoepidermoid carcinomas, four acinic cell carcinomas, four pleomorphic adenomas, two adenoid cystic carcinomas, two NUT carcinomas, and one each of basal cell adenoma, salivary duct carcinoma ex-pleomorphic adenoma, salivary duct carcinoma, clear cell carcinoma, and secretory carcinoma. Directed FISH testing based on the diagnosis was performed on cases that did not already have FISH conducted during clinical workup. Targeted RNA sequencing of 507 genes and their partners (using the Illumina TruSight Fusion Panel) was completed. Six of twenty-two (27.3%) cases had discordant results. In three cases, FISH results were negative while RNA sequencing results found fusion transcripts, which were all confirmed with RT-PCR and Sanger sequencing. In three cases, RNA sequencing results were negative while FISH results were positive for a gene rearrangement. Thus, if fusion analysis results are conflicting with the morphologic impression, a second mode of fusion detection may be warranted. Although both methods have advantages and drawbacks, RNA sequencing provides additional information about novel fusion partners and fusions that may not have been originally considered.

Glioblastoma (GBM) is the most common primary malignant brain tumor in adults. Prior studies examining the mutational landscape of GBM revealed recurrent alterations in genes that regulate the same growth control pathways. To this regard, ~ 40% of GBM harbor EGFR alterations, whereas BRAF variants are rare. Existing data suggests that gain-of-function mutations in these genes are mutually exclusive. This study was designed to explore the clinical, pathological, and molecular differences between EGFR- and BRAF-mutated GBM. We reviewed retrospective clinical data from 89 GBM patients referred for molecular testing between November 2012 and December 2015. Differences in tumor mutational profile, location, histology, and survival outcomes were compared in patients with EGFR- versus BRAF-mutated tumors, and microarray data from The Cancer Genome Atlas was used to assess differential gene expression between the groups. Individuals with BRAF-mutant tumors were typically younger and survived longer relative to those with EGFR-mutant tumors, even in the absence of targeted treatments. BRAF-mutant tumors lacked distinct histomorphology but exhibited unique localization in the brain, typically arising adjacent to the lateral ventricles. Compared to EGFR- and IDH1-mutant tumors, BRAF-mutant tumors showed increased expression of genes related to a trophoblast-like phenotype, specifically HLA-G and pregnancy specific glycoproteins, that have been implicated in invasion and immune evasion. Taken together, these observations suggest a distinct clinical presentation, brain location, and gene expression profile for BRAF-mutant tumors. Pending further study, this may prove useful in the stratification and management of GBM.

Severe congenital neutropenia (SCN) is a collection of diverse disorders characterized by chronically low absolute neutrophil count in the peripheral blood, increased susceptibility to infection, and a significant predisposition to the development of myeloid malignancies. SCN can be acquired or inherited. Inherited forms have been linked to variants in a group of diverse genes involved in the neutrophil-differentiation process. Variants that promote resistance to treatment have also been identified. Thus, genetic testing is important for the diagnosis, prognosis, and management of SCN. Herein we describe clinically validated assay developed for assessing patients with suspected SCN. The assay is performed from a whole-exome backbone. Variants are called across all coding exons, and results are filtered to focus on 48 genes that are clinically relevant to SCN. Validation results indicated 100% analytical sensitivity and specificity for the detection of constitutional variants among the 48 reportable genes. To date, 34 individuals have been referred for testing (age range: birth to 67 years). Several pathogenic and likely pathogenic variants have been identified, including one in a patient with late-onset disease. The pattern of cases referred for testing suggests that this assay has clinical utility in a broader spectrum of patients beyond those in the pediatric population who have classic early-onset symptoms characteristic of SCN.

Ependymomas are neuroepithelial tumors that differentiate along the ependymal cell lineage, a lining of the ventricles of the brain and the central canal of the spinal cord. They are rare in adults, but account for around 9% of brain tumors in children, where they usually have an aggressive course. Efficient stratification could lead to improved care but remains a challenge even in the genomic era. Recent studies proposed a multivariate classification system based on tumor location, age, and broad genomic findings like global patterns of methylation and copy number variants (CNVs). This system shows improved prognostic utility, but is relatively impractical in the routine clinical setting because it necessitates multiple diagnostic tests. We analyzed 13 intracranial grade II and III ependymoma specimens on a DNA microarray to identify discrete CNVs that could support the existing classification. The loss of chr22 and the gain of 5p15.31 were common throughout our cohort (6 and 11 cases, respectively). Other CNVs correlated well with the previously proposed classification system. For example, gains of chr20 were unique to PF-EPN-B tumors of the posterior fossa and may differentiate them from PF-EPN-A. Given the ease of detecting CNVs using multiple, clinically validated methods, these CNVs should be further studied to confirm their diagnostic and prognostic utility, for incorporation into clinical testing algorithms.

Disorders of somatic mosaicism (DoSM) are a diverse group of syndromic and non-syndromic conditions caused by mosaic variants in genes that regulate cell survival and proliferation. Despite overlap in gene space and technical requirements, few clinical labs specialize in DoSM compared to oncology. We adapted a high-sensitivity next-generation sequencing cancer assay for DoSM in 2014. Some 343 individuals have been tested over the past 5 years, 58% of which had pathogenic and likely pathogenic (P/LP) findings, for a total of 206 P/LP variants in 22 genes. Parameters associated with the high diagnostic yield were: (1) deep sequencing (∼2,000× coverage), (2) a broad gene set, and (3) testing affected tissues. Fresh and formalin-fixed paraffin embedded tissues performed equivalently for identification of P/LP variants (62% and 71% of individuals, respectively). Comparing cultured fibroblasts to skin biopsies suggested that culturing might boost the allelic fraction of variants that confer a growth advantage, specifically gain-of-function variants in PIK3CA. Buccal swabs showed high diagnostic sensitivity in case subjects where disease phenotypes manifested in the head or brain. Peripheral blood was useful as an unaffected comparator tissue to determine somatic versus constitutional origin but had poor diagnostic sensitivity. Descriptions of all tested individuals, specimens, and P/LP variants included in this cohort are available to further the study of the DoSM population.

Meningiomas are among the most common tumors of the adult central nervous system (CNS). They are classified by the World Health Organization into three pathologic grades with increasing severity: grade I are benign with favorable treatment outcomes and low recurrence rates while grade III display malignant behavior and poor progression-free survival. Previous studies have shown that inactivation of NF-2 is the most common genetic event in high-grade meningioma; however, there is dearth of molecular data to distinguish grade II (AM-II) from the even more aggressive grade III (AM-III). As part of a routine diagnostic workup, 19 AM-II and 5 AM-III were submitted for targeted sequencing on a panel of twenty-four genes relevant to CNS tumors. The data generated during the course of clinical care was collected and re-analyzed with the aim of identifying molecular features to distinguish AM-II and AM-III. Our cases contained several well-characterized, potentially actionable mutations, but we did not find any novel, recurrent sequence variants. Copy number variations were common in both AM-II and AM-III; chr22q loss was the most prevalent followed in decreasing frequency by losses of chr1p, chr14q, and chr10. In particular, chr10 loss was noted in 4 of 5 AM-III cases but none of the AM-II cases. This suggests that chr10 loss may serve as a diagnostic and perhaps a prognostic marker to differentiate AM-II from AM-III. If confirmed in larger studies, our finding could further aid the classification of meningioma.

The detection of recurrent gene fusions can help confirm diagnoses in solid tumors, particularly when the morphology and staining are unusual or nonspecific, and can guide therapeutic decisions. Although fluorescence in situ hybridization and PCR are often used to identify fusions, the rearrangement must be suspected, with only a few prioritized probes run. It was hypothesized that the Illumina TruSight RNA Fusion Panel, which detects fusions of 507 genes and their partners, would uncover fusions with greater sensitivity than other approaches, leading to changes in diagnosis, prognosis, or therapy. Targeted RNA sequencing was performed on formalin-fixed, paraffin-embedded sarcoma and carcinoma cases in which fluorescence in situ hybridization, RT-PCR, or DNA-based sequencing was conducted during the diagnostic workup. Of the 153 cases, 138 (90%) were sequenced with adequate quality control metrics. A total of 101 of 138 (73%) cases were concordant by RNA sequencing and the prior test method. RNA sequencing identified an additional 30 cases (22%) with fusions that were not detected by conventional methods. In seven cases (5%), the additional fusion information provided by RNA sequencing would have altered diagnosis and management. A total of 19 novel fusion pairs (not previously described in the literature) were discovered (14%). Overall, the findings show that a targeted RNA-sequencing method can detect gene fusions in formalin-fixed, paraffin-embedded specimens with high sensitivity.

Background Paragonimus spp. (lung flukes) are among the most injurious foodborne helminths, infecting ∼23 million people and subjecting ∼292 million to infection risk. Paragonimiasis is acquired from infected undercooked crustaceans and primarily affects the lungs but often causes lesions elsewhere including the brain. The disease is easily mistaken for tuberculosis owing to similar pulmonary symptoms, and accordingly, diagnostics are in demand. Results We assembled, annotated, and compared draft genomes of 4 prevalent and distinct Paragonimus species: Paragonimus miyazakii, Paragonimus westermani, Paragonimus kellicotti, and Paragonimus heterotremus. Genomes ranged from 697 to 923 Mb, included 12,072–12,853 genes, and were 71.6–90.1% complete according to BUSCO. Orthologous group analysis spanning 21 species (lung, liver, and blood flukes, additional platyhelminths, and hosts) provided insights into lung fluke biology. We identified 256 lung fluke–specific and conserved orthologous groups with consistent transcriptional adult-stage Paragonimus expression profiles and enriched for iron acquisition, immune modulation, and other parasite functions. Previously identified Paragonimus diagnostic antigens were matched to genes, providing an opportunity to optimize and ensure pan-Paragonimus reactivity for diagnostic assays. Conclusions This report provides advances in molecular understanding of Paragonimus and underpins future studies into the biology, evolution, and pathogenesis of Paragonimus and related foodborne flukes. We anticipate that these novel genomic and transcriptomic resources will be invaluable for future lung fluke research.

Severe congenital neutropenia (SCN) is a collection of diverse disorders characterized by chronically low absolute neutrophil count in the peripheral blood, increased susceptibility to infection, and a significant predisposition to the development of myeloid malignancies. SCN can be acquired or inherited. Inherited forms have been linked to variants in a group of diverse genes involved in the neutrophil-differentiation process. Variants that promote resistance to treatment have also been identified. Thus, genetic testing is important for the diagnosis, prognosis, and management of SCN. Herein we describe clinically validated assay developed for assessing patients with suspected SCN. The assay is performed from a whole-exome backbone. Variants are called across all coding exons, and results are filtered to focus on 48 genes that are clinically relevant to SCN. Validation results indicated 100% analytical sensitivity and specificity for the detection of constitutional variants among the 48 reportable genes. To date, 34 individuals have been referred for testing (age range: birth to 67 years). Several pathogenic and likely pathogenic variants have been identified, including one in a patient with late-onset disease. The pattern of cases referred for testing suggests that this assay has clinical utility in a broader spectrum of patients beyond those in the pediatric population who have classic early-onset symptoms characteristic of SCN.

Meningiomas are common in adults (~35% of brain tumors) but rare in children, where they exhibit unique clinical, pathological and molecular features compared to adult counterparts. Thus, data generated from adult cohorts may be imperfectly suited to guiding diagnostic, prognostic and treatment decisions for children. We studied 50 meningioma patients ≤18 years with available clinical and pathological data to address the need for data obtained in the pediatric setting. As previously described, we noted a slight bias toward male patients and a higher proportion of spinal tumors compared to adults. Thirty-eight of 50 specimens were further analyzed by next generation sequencing. Loss-of-function mutations in NF2 and chromosome 22 losses were common, but pathogenic variants in other genes (SMARCB1, FUBP1, BRAF, TERT promoter, CHEK2, SMAD and GATA3) were identified in a minority of cases. Copy number variants outside of chromosomes 22 and 1 were infrequent. H3K27 hypomethylation, a useful biomarker in adult tumors, was not found in our cohort. In exploring the correlation between mitotic count and recurrence-free survival, we found a threshold of six mitoses per 10 high powered fields as the optimal cutoff in predicting recurrence-free survival. If independently validated in larger studies, adjusted grading thresholds could enhance the clinical management of pediatric meningiomas.

Background Next-generation sequencing (NGS) is a useful tool for evaluating patients with suspected genetic kidney disease. Clinical practice relies on the use of targeted gene panels that are ordered based on patient presentation. We compare the diagnostic yield of clinical panel-based testing to exome analysis. Methods In total, 324 consecutive patients underwent physician-ordered, panel-based NGS testing between December 2014 and October 2018. Gene panels were available for four clinical phenotypes, including atypical hemolytic uremic syndrome (n=224), nephrotic syndrome (n=56), cystic kidney disease (n=26), and Alport syndrome (n=13). Variants were analyzed and clinical reports were signed out by a pathologist or clinical geneticist at the time of testing. Subsequently, all patients underwent retrospective exome analysis to detect additional clinically significant variants in kidney disease genes that were not analyzed as part of the initial clinical gene panel. Resulting variants were classified according to the American College of Medical Genetics and Genomics 2015 guidelines. Results In the initial physician-ordered gene panels, we identified clinically significant pathogenic or likely pathogenic variants in 13% of patients (n=42/324). CFHR3-CFHR1 homozygous deletion was detected in an additional 13 patients with aHUS without a pathogenic or likely pathogenic variant. Diagnostic yield of the initial physician-ordered gene panel was 20% and varied between groups. Retrospective exome analysis identified 18 patients with a previously unknown pathogenic or likely pathogenic variant in a kidney disease gene and eight patients with a high-risk APOL1 genotype. Overall, retrospective exome analysis increased the diagnostic yield of panel-based testing from 20% to 30%. Conclusions These results highlight the importance of a broad and collaborative approach between the clinical laboratory and their physician clients that employs additional analysis when a targeted panel of kidney disease–causing genes does not return a clinically meaningful result.

A subset of endometrial adenocarcinomas (EACs) exhibit microsatellite instability and have deficient DNA mismatch repair (dMMR). The overall aim of the study was to compare the spectrum of mutations in endometrioid-type EAC with and without dMMR by using a clinically validated next-generation sequencing assay. We retrospectively identified 19 EACs with known mismatch repair status that had undergone targeted sequencing of a panel of cancer-related genes. The mismatch repair status was ascertained by immunohistochemistry against MLH1, PMS2, MSH2, and MSH6 mismatch proteins. Somatic mutations in EAC with dMMR were compared against those in cases with proficient MMR (pMMR). The dMMR EAC showed a normalized mean of 66.6 mutations/Mb per case compared with pMMR EAC with a mean of 26.2 (P

The 2007 World Health Organization Classification of Tumours of the Central Nervous System classifies lower-grade gliomas [LGGs (grades II to III diffuse gliomas)] morphologically as astrocytomas or oligodendrogliomas, and tumors with unclear ambiguous morphology as oligoastrocytomas. The World Health Organization's newly released (2016) classification incorporates molecular data. A single, targeted next-generation sequencing (NGS) panel was used for detecting single-nucleotide variation and copy number variation in 50 LGG cases originally classified using the 2007 criteria, including 36 oligoastrocytomas, 11 oligodendrogliomas, 2 astrocytomas, and 1 LGG not otherwise specified. NGS results were compared with those from IHC analysis and fluorescence in situ hybridization to assess concordance and to categorize the tumors according to the 2016 criteria. NGS results were concordant with those from IHC analysis in all cases. In 3 cases, NGS was superior to fluorescence in situ hybridization in distinguishing segmental chromosomal losses from whole-arm deletions. The NGS approach was effective in reclassifying 36 oligoastrocytomas as 30 astrocytomas (20 IDH1/2 mutant and 10 IDH1/2 wild type) and 6 oligodendrogliomas, and 1 oligodendroglioma as an astrocytoma (IDH1/2 mutant). Here we show that a single, targeted NGS assay can serve as the sole testing modality for categorizing LGG according to the World Health Organization's 2016 diagnostic scheme. This modality affords greater accuracy and efficiency while reducing specimen tissue requirements compared with multimodal approaches.

PCR amplification, a key step in next-generation sequencing (NGS) library construction, can generate an unlimited amount of product from limited input; however, it cannot create more information than was present in the original template. Thus, NGS libraries can be made from very little DNA, but reducing the input may compromise assay sensitivity in ways that are difficult to ascertain unless library complexity (ie, the number of unique DNA molecules represented in the library) and depth of coverage with unique sequence reads (those derived from input DNA molecules) versus duplicate sequence reads (those resulting from overamplification of particular molecules) are discretely measured. A series of experiments was performed to explore the impact of low DNA input on an amplicon-based NGS assay using unique molecular identifiers to track unique versus duplicate reads. At high sequencing depths, unique and total (unique plus duplicate) read coverage are not well correlated, so increasing the number of sequenced reads does not necessarily improve sensitivity. Unique coverage depth tends to improve with more input, but improvements are not consistent. Fluctuations in library complexity complicated variant detection using both standardized and clinical specimens, often resulting in technical replicates with vastly different estimates of variant allelic fraction. In conclusion, depth of coverage with unique reads must be tracked in clinical NGS to ensure that sensitivity and accuracy are maintained.

Improved diagnostic methods are needed to support ongoing efforts to eliminate onchocerciasis (river blindness). This study used an integrated approach to identify adult female Onchocerca volvulus antigens that can be explored for developing serodiagnostic tests. The first step was to develop a detailed multi-omics database of all O. volvulus proteins deduced from the genome, gene transcription data for different stages of the parasite including eight individual female worms (providing gene expression information for 94.8% of all protein coding genes), and the adult female worm proteome (detecting 2126 proteins). Next, female worm proteins were purified with IgG antibodies from onchocerciasis patients and identified using LC-MS with a high-resolution hybrid quadrupole-time-of-flight mass spectrometer. A total of 241 immunoreactive proteins were identified among those bound by IgG from infected individuals but not IgG from uninfected controls. These included most of the major diagnostic antigens described over the past 25 years plus many new candidates. Proteins of interest were prioritized for further study based on a lack of conservation with orthologs in the human host and other helminthes, their expression pattern across the life cycle, and their consistent expression among individual female worms. Based on these criteria, we selected 33 proteins that should be carried forward for testing as serodiagnostic antigens to supplement existing diagnostic tools. These candidates, together with the extensive pan-omics dataset generated in this study are available to the community (http://nematode.net) to facilitate basic and translational research on onchocerciasis.

Glioblastoma is the most common primary malignancy of the adult central nervous system. Gliomagenesis involves a complex range of alterations, including sequence changes, copy number variations (CNVs), and epigenetic modifications, that have clinical implications for disease classification and prognosis. Thus, multiple testing modalities are required to support a complete diagnostic workup. The goal of this study was to streamline the multipart workflow by predicting both sequence changes and CNVs (specifically EGFR amplifications) from a single next-generation sequencing (NGS) test. Eighty-six primary and secondary glioblastomas were submitted for clinical NGS to report sequence variants from a concise panel of cancer-relevant genes. Most specimens underwent concomitant testing by methylation-specific polymerase chain reaction, immunohistochemistry, and fluorescence in situ hybridization. Using data generated during the course of clinical testing, we found that NGS-based variant predictions were concordant with immunohistochemistry and fluorescence in situ hybridization for IDH mutation and EGFR amplification status, respectively. We also noted that EGFR amplifications correlated with polysomy of chromosome 7, 19, and 20, and loss of PTEN and CDKN2A. EGFR-unamplified cases had lower rates of chromosome 7 polysomy, and PTEN and CDKN2A loss, but more CNVs overall. TP53, NF1, ATRX, and PDGFRA mutations were nearly exclusive to specimens without EGFR amplification. EGFR amplification was not associated with longer progression-free survival in this cohort, but amplifications were enriched in a group with slightly longer overall survival despite radiographic evidence of disease progression. Further study is needed to explore the mechanisms responsible for noted patterns of co-occurring variants and to correlate them with specific clinical outcomes.

In lung adenocarcinoma, canonical EML4-ALK inversion results in a fusion protein with a constitutively active ALK kinase domain. Evidence of ALK rearrangement occurs in a minority (2-7%) of lung adenocarcinoma, and only ~60% of these patients will respond to targeted ALK inhibition by drugs such as crizotinib and ceritinib. Clinically, targeted anti-ALK therapy is often initiated based on evidence of an ALK genomic rearrangement detected by fluorescence in situ hybridization (FISH) of interphase cells in formalin-fixed, paraffin-embedded tissue sections. At the genomic level, however, ALK rearrangements are heterogeneous, with multiple potential breakpoints in EML4, and alternate fusion partners. Using next-generation sequencing of DNA and RNA together with ALK immunohistochemistry, we comprehensively characterized genomic breakpoints in 33 FISH-positive lung adenocarcinomas. Of these 33 cases, 29 (88%) had detectable DNA level ALK rearrangements involving EML4, KIF5B, or non-canonical partners including ASXL2, ATP6V1B1, PRKAR1A, and SPDYA. A subset of 12 cases had material available for RNA-Seq. Of these, eight of eight (100%) cases with DNA rearrangements showed ALK fusion transcripts from RNA-Seq; three of four cases (75%) without detectable DNA rearrangements were similarly negative by RNA-Seq, and one case was positive by RNA-Seq but negative by DNA next-generation sequencing. By immunohistochemistry, 17 of 19 (89%) tested cases were clearly positive for ALK protein expression; the remaining cases had no detectable DNA level rearrangement or had a non-canonical rearrangement not predicted to form a fusion protein. Survival analysis of patients treated with targeted ALK inhibitors demonstrates a significant difference in mean survival between patients with next-generation sequencing confirmed EML4-ALK rearrangements, and those without (20.6 months vs 5.4 months, P

Ascaris suum provides a powerful model for studying parasitic nematodes, including individual tissues such as the intestine, an established target for anthelmintic treatments. Here, we add a valuable experimental component to our existing functional, proteomic, transcriptomic and phylogenomic studies of the Ascaris suum intestine, by developing a method to manipulate intestinal cell functions via direct delivery of experimental treatments (in this case, double-stranded (ds)RNA) to the apical intestinal membrane. We developed an intestinal perfusion method for direct, controlled delivery of dsRNA/heterogeneous small interfering (hsi) RNA into the intestinal lumen for experimentation. RNA-Seq (22 samples) was used to assess influences of the method on global intestinal gene expression. Successful mRNA-specific knockdown in intestinal cells of adult A. suum was accomplished with this new experimental method. Global transcriptional profiling confirmed that targeted transcripts were knocked down more significantly than any others, with only 12 (0.07% of all genes) or 238 (1.3%) off-target gene transcripts consistently differentially regulated by dsRNA treatment or the perfusion experimental design, respectively (after 24 h). The system supports controlled, effective delivery of treatments (dsRNA/hsiRNA) to the apical intestinal membrane with relatively minor off-target effects, and builds on our experimental model to dissect A. suum intestinal cell functions with broad relevance to parasitic nematodes.

Helminth.net (http://www.helminth.net) is the new moniker for a collection of databases: Nematode.net and Trematode.net. Within this collection we provide services and resources for parasitic roundworms (nematodes) and flatworms (trematodes), collectively known as helminths. For over a decade we have provided resources for studying nematodes via our veteran site Nematode.net (http://nematode.net). In this article, (i) we provide an update on the expansions of Nematode.net that hosts omics data from 84 species and provides advanced search tools to the broad scientific community so that data can be mined in a useful and user-friendly manner and (ii) we introduce Trematode.net, a site dedicated to the dissemination of data from flukes, flatworm parasites of the class Trematoda, phylum Platyhelminthes. Trematode.net is an independent component of Helminth.net and currently hosts data from 16 species, with information ranging from genomic, functional genomic data, enzymatic pathway utilization to microbiome changes associated with helminth infections. The databases' interface, with a sophisticated query engine as a backbone, is intended to allow users to search for multi-factorial combinations of species' omics properties. This report describes updates to Nematode.net since its last description in NAR, 2012, and also introduces and presents its new sibling site, Trematode.net.

Ongoing elimination efforts have altered the global distribution of Onchocerca volvulus, the agent of river blindness, and further population restructuring is expected as efforts continue. Therefore, a better understanding of population genetic processes and their effect on biogeography is needed to support elimination goals. We describe O. volvulus genome variation in 27 isolates from the early 1990s (before widespread mass treatment) from four distinct locales: Ecuador, Uganda, the West African forest and the West African savanna. We observed genetic substructuring between Ecuador and West Africa and between the West African forest and savanna bioclimes, with evidence of unidirectional gene flow from savanna to forest strains. We identified forest:savanna-discriminatory genomic regions and report a set of ancestry informative loci that can be used to differentiate between forest, savanna and admixed isolates, which has not previously been possible. We observed mito-nuclear discordance possibly stemming from incomplete lineage sorting. The catalogue of the nuclear, mitochondrial and endosymbiont DNA variants generated in this study will support future basic and translational onchocerciasis research, with particular relevance for ongoing control programmes, and boost efforts to characterize drug, vaccine and diagnostic targets.

Background Paragonimiasis is an important and widespread neglected tropical disease. Fifteen Paragonimus species are human pathogens, but two of these, Paragonimus westermani and P. skrjabini, are responsible for the bulk of human disease. Despite their medical and economic significance, there is limited information on the gene content and expression of Paragonimus lung flukes. Results The transcriptomes of adult P. westermani and P. skrjabini were studied with deep sequencing technology. Approximately 30 million reads per species were assembled into 21,586 and 25,825 unigenes for P. westermani and P. skrjabini, respectively. Many unigenes showed homology with sequences from other food-borne trematodes, but 1,217 high-confidence Paragonimus-specific unigenes were identified. Analyses indicated that both species have the potential for aerobic and anaerobic metabolism but not de novo fatty acid biosynthesis and that they may interact with host signaling pathways. Some 12,432 P. westermani and P. skrjabini unigenes showed a clear correspondence in bi-directional sequence similarity matches. The expression of shared unigenes was mostly well correlated, but differentially expressed unigenes were identified and shown to be enriched for functions related to proteolysis for P. westermani and microtubule based motility for P. skrjabini. Conclusions The assembled transcriptomes of P. westermani and P. skrjabini, inferred proteins, and extensive functional annotations generated for this project (including identified primary sequence similarities to various species, protein domains, biological pathways, predicted proteases, molecular mimics and secreted proteins, etc.) represent a valuable resource for hypothesis driven research on these medically and economically important species. Electronic supplementary material The online version of this article (doi:10.1186/s13071-016-1785-x) contains supplementary material, which is available to authorized users.

The 2007 World Health Organization Classification of Tumours of the Central Nervous System classifies lower-grade gliomas [LGGs (grades II to III diffuse gliomas)] morphologically as astrocytomas or oligodendrogliomas, and tumors with unclear ambiguous morphology as oligoastrocytomas. The World Health Organization's newly released (2016) classification incorporates molecular data. A single, targeted next-generation sequencing (NGS) panel was used for detecting single-nucleotide variation and copy number variation in 50 LGG cases originally classified using the 2007 criteria, including 36 oligoastrocytomas, 11 oligodendrogliomas, 2 astrocytomas, and 1 LGG not otherwise specified. NGS results were compared with those from IHC analysis and fluorescence in situ hybridization to assess concordance and to categorize the tumors according to the 2016 criteria. NGS results were concordant with those from IHC analysis in all cases. In 3 cases, NGS was superior to fluorescence in situ hybridization in distinguishing segmental chromosomal losses from whole-arm deletions. The NGS approach was effective in reclassifying 36 oligoastrocytomas as 30 astrocytomas (20 IDH1/2 mutant and 10 IDH1/2 wild type) and 6 oligodendrogliomas, and 1 oligodendroglioma as an astrocytoma (IDH1/2 mutant). Here we show that a single, targeted NGS assay can serve as the sole testing modality for categorizing LGG according to the World Health Organization's 2016 diagnostic scheme. This modality affords greater accuracy and efficiency while reducing specimen tissue requirements compared with multimodal approaches.

The bovine lungworm, Dictyocaulus viviparus (order Strongylida), is an important parasite of livestock that causes substantial economic and production losses worldwide. Here we report the draft genome, variome and developmental transcriptome of D. viviparus. The genome (161 Mb) is smaller than those of related bursate nematodes and encodes fewer proteins (14,171 total). In the first genome-wide assessment of genomic variation in any parasitic nematode, we found a high degree of sequence variability in proteins predicted to be involved host-parasite interactions. Next, we used extensive RNA sequence data to track gene transcription across the life cycle of D. viviparus and identified genes that might be important in nematode development and parasitism. Finally, we predicted genes that could be vital in host-parasite interactions, genes that could serve as drug targets and putative RNAi effectors with a view to developing functional genomic tools. This extensive, well-curated dataset should provide a basis for developing new anthelmintics, vaccines and improved diagnostic tests and serve as a platform for future investigations of drug resistance and epidemiology of the bovine lungworm and related nematodes.

Chandlerella quiscali is a filarial nematode parasite of the common grackle (Quiscalus quiscula), a widespread bird species found throughout most of North America. Worms collected from wild-caught birds were morphologically identified as C. quiscali and tested for the presence of Wolbachia, an alphaproteobacterial endosymbiont required for reproduction and maturation by many filarial species. Although various methods, including polymerase chain reaction, in situ hybridization, and immunohistology, were used, we were unable to detect evidence of colonization with Wolbachia. Due to the widespread distribution of the grackle host, localization within the host, and high prevalence, C. quiscali may be among the most easily obtainable of Wolbachia-free filarial species. Further studies of C. quiscali and other Wolbachia-free filarial species may help to clarify the reason(s) that some filarial species require Wolbachia but others do not.

Genetic characterization of field isolates and clinical specimens of filarial nematodes is often limited by a shortage of DNA; therefore, we evaluated a multiple displacement amplification (MDA) based whole genome amplification method. The quality of amplified DNA was examined by conventional PCR, real-time PCR, and DNA hybridization. MDA of 5.0 ng of adult Brugia malayi DNA and one-fifteenth of the DNA isolated from a single microfilaria resulted in 6.3 and 4.2 microg of amplified DNA, respectively. Amplified DNA was equivalent to native genomic DNA for hybridization to B. malayi BAC library clones or to an oligonucleotide microarray with approximately 18,000 filarial DNA sequences. MDA is useful for whole genome amplification of filarial DNA from very small amounts of starting material. This technology will permit detailed studies of genetic diversity that were not previously feasible.

FGFR2 is recurrently amplified in 5% of gastric cancers and 1%–4% of breast cancers; however, this molecular alteration has never been reported in a primary colorectal cancer specimen. Preclinical studies indicate that several FGFR tyrosine-kinase inhibitors (TKIs), such as AZD4547, have in vitro activity against the FGFR2-amplified colorectal cell line, NCI-H716. The efficacy of these inhibitors is currently under investigation in clinical trials for breast and gastric cancer. Thus, better characterizing colorectal tumors for FGFR2 amplification could identify a subset of patients who may benefit from FGFR TKI therapies. Here, we describe a novel FGFR2 amplification identified by clinical next-generation sequencing in a primary colorectal cancer. Further characterization of the tumor by immunohistochemistry showed neuroendocrine differentiation, similar to the reported properties of the NCI-H716 cell line. These findings demonstrate that the spectrum of potentially clinically actionable mutations detected by targeted clinical sequencing panels is not limited to only single-nucleotide polymorphisms and insertions/deletions but also to copy-number alterations.

Food borne trematodes (FBTs) are an assemblage of platyhelminth parasites transmitted through the food chain, four of which are recognized as neglected tropical diseases (NTDs). Fascioliasis stands out among the other NTDs due to its broad and significant impact on both human and animal health, as Fasciola sp., are also considered major pathogens of domesticated ruminants. Here we present a reference genome sequence of the common liver fluke, Fasciola hepatica isolated from sheep, complementing previously reported isolate from cattle. A total of 14,642 genes were predicted from the 1.14 GB genome of the liver fluke. Comparative genomics indicated that F. hepatica Oregon and related food-borne trematodes are metabolically less constrained than schistosomes and cestodes, taking advantage of the richer millieux offered by the hepatobiliary organs. Protease families differentially expanded between diverse trematodes may facilitate migration and survival within the heterogeneous environments and niches within the mammalian host. Surprisingly, the sequencing of Oregon and Uruguay F. hepatica isolates led to the first discovery of an endobacteria in this species. Two contigs from the F. hepatica Oregon assembly were joined to complete the 859,205 bp genome of a novel Neorickettsia endobacterium (nFh) closely related to the etiological agents of human Sennetsu and Potomac horse fevers. Immunohistochemical studies targeting a Neorickettsia surface protein found nFh in specific organs and tissues of the adult trematode including the female reproductive tract, eggs, the Mehlis’ gland, seminal vesicle, and oral suckers, suggesting putative routes for fluke-to-fluke and fluke-to-host transmission. The genomes of F. hepatica and nFh will serve as a resource for further exploration of the biology of F. hepatica, and specifically its newly discovered trans-kingdom interaction with nFh and the impact of both species on disease in ruminants and humans., Author Summary This report presents novel findings revealing (a) the genome sequence of the food-borne trematode Fasciola hepatica (the liver fluke) isolated from sheep, which stands out among neglected tropical diseases due to its zoonotic impact on both human and animal health and (b) the first instance (and the genome) of the rickettsial endobacterium of the genus Neorickettsia in F. hepatica. Using stage-specific gene expression data, we identified liver fluke proteins likely involved in host-parasite interactions, and using immunolocalization, we confirmed Neorickettsia in organs and tissues of the adult trematode. The presence of the bacteria in fluke reproductive tissues and eggs suggests a possible mechanism for vertical transmission, and the presence of bacteria in the oral sucker used to anchor flukes to the lining of the biliary tract suggests a potential mechanism for horizontal transmission to the mammalian host. This is of interest because related Neorickettsia cause severe, even deadly, illness in a variety of species, including humans. This is the first report to localize Neorickettsia endobacteria within the tissues of adult F. hepatica. The discoveries in our manuscript have wide impact for the fields of both the pathophysiology and evolution of Fasciola and related FBTs, and the transmission strategies of Neorickettsia.

Background Paragonimiasis is a food-borne trematode infection acquired by eating raw or undercooked crustaceans. It is a major public health problem in the far East, but it also occurs in South Asia, Africa, and in the Americas. Paragonimus worms cause chronic lung disease with cough, fever and hemoptysis that can be confused with tuberculosis or other non-parasitic diseases. Treatment is straightforward, but diagnosis is often delayed due to a lack of reliable parasitological or serodiagnostic tests. Hence, the purpose of this study was to use a systems biology approach to identify key parasite proteins that may be useful for development of improved diagnostic tests. Methodology/Principal Findings The transcriptome of adult Paragonimus kellicotti was sequenced with Illumina technology. Raw reads were pre-processed and assembled into 78,674 unique transcripts derived from 54,622 genetic loci, and 77,123 unique protein translations were predicted. A total of 2,555 predicted proteins (from 1,863 genetic loci) were verified by mass spectrometric analysis of total worm homogenate, including 63 proteins lacking homology to previously characterized sequences. Parasite proteins encoded by 321 transcripts (227 genetic loci) were reactive with antibodies from infected patients, as demonstrated by immunoaffinity purification and high-resolution liquid chromatography-mass spectrometry. Serodiagnostic candidates were prioritized based on several criteria, especially low conservation with proteins in other trematodes. Cysteine proteases, MFP6 proteins and myoglobins were abundant among the immunoreactive proteins, and these warrant further study as diagnostic candidates. Conclusions The transcriptome, proteome and immunolome of adult P. kellicotti represent a major advance in the study of Paragonimus species. These data provide a powerful foundation for translational research to develop improved diagnostic tests. Similar integrated approaches may be useful for identifying novel targets for drugs and vaccines in the future., Author Summary Paragonimiasis is a food-borne trematode infection that people acquire when they eat raw or undercooked crustaceans. Disease symptoms (including cough, fever, blood in sputum, etc.) can be similar to those observed in patients with tuberculosis or bacterial pneumonia, frequently resulting in misdiagnosis. Although the infection is relatively easy to treat, diagnosis is complicated. Available diagnostic assays rely on total parasite homogenate to facilitate the detection of Paragonimus-specific antibodies in patients. Though these blot-based assays have shown high sensitivity and specificity, they are inconvenient because total parasite homogenate is not readily available. This study used next generation genomic and proteomic methods to identify transcripts and proteins expressed in adult Paragonimus flukes. We then used sera from patients infected with P. kellicotti to isolate immunoreactive proteins, and these were analyzed by mass spectrometry. The annotated transcriptome and the associated proteome of the antibody immune response represent a significant advance in research on Paragonimus. This information will be a valuable resource for further research on Paragonimus and paragonimiasis. Thus this project illustrates the potential power of employing systems biology for translational research in parasitology.

Background Alternative splicing (AS) of mRNA is a vital mechanism for enhancing genomic complexity in eukaryotes. Spliced isoforms of the same gene can have diverse molecular and biological functions and are often differentially expressed across various tissues, times, and conditions. Thus, AS has important implications in the study of parasitic nematodes with complex life cycles. Transcriptomic datasets are available from many species, but data must be revisited with splice-aware assembly protocols to facilitate the study of AS in helminthes. Methods We sequenced cDNA from the model worm Caenorhabditis elegans using 454/Roche technology for use as an experimental dataset. Reads were assembled with Newbler software, invoking the cDNA option. Several combinations of parameters were tested and assembled transcripts were verified by comparison with previously reported C. elegans genes and transcript isoforms and with Illumina RNAseq data. Results Thoughtful adjustment of program parameters increased the percentage of assembled transcripts that matched known C. elegans sequences, decreased mis-assembly rates (i.e., cis- and trans-chimeras), and improved the coverage of the geneset. The optimized protocol was used to update de novo transcriptome assemblies from nine parasitic nematode species, including important pathogens of humans and domestic animals. Our assemblies indicated AS rates in the range of 20-30%, typically with 2-3 transcripts per AS locus, depending on the species. Transcript isoforms from the nine species were translated and searched for similarity to known proteins and functional domains. Some 21 InterPro domains, including several involved in nucleotide and chromatin binding, were statistically correlated with AS genetic loci. In most cases, the Roche/454 data explored in this study are the only sequences available from the species in question; however, the recently published genome of the human hookworm Necator americanus provided an additional opportunity to validate our results. Conclusions Our optimized assembly parameters facilitated the first survey of AS among parasitic nematodes. The nine transcriptome assemblies, their protein translations, and basic annotations are available from Nematode.net as a resource for the research community. These should be useful for studies of specific genes and gene families of interest as well as for curating draft genome assemblies as they become available.

Lymphatic filariasis is caused by three closely related nematode parasites: Wuchereria bancrofti, Brugia malayi and Brugia timori. These species have many ecological variants that differ in several aspects of their biology such as mosquito vector species, host range, periodicity, and morphology. Although the genome of B. malayi (the first genome sequenced from a parasitic nematode) has been available for more than five years, very little is known about genetic variability among the lymphatic dwelling filariae. The genetic diversity among these worms is not only interesting from a biological perspective, but it may have important practical implications for the Global Program to Eliminate Lymphatic Filariasis, as the parasites may respond differently to diagnostic tests and/or medical interventions. Therefore, better information on their genetic variability is urgently needed. With improved methods for nucleic acid extraction and recent advances in sequencing chemistry and instrumentation, this gap can be filled relatively inexpensively. Improved information on filarial genetic diversity may increase the chances of success for lymphatic filariasis elimination programs.

Background Many species of filarial nematodes depend on Wolbachia endobacteria to carry out their life cycle. Other species are naturally Wolbachia-free. The biological mechanisms underpinning Wolbachia-dependence and independence in filarial nematodes are not known. Previous studies have indicated that Wolbachia have an impact on mitochondrial gene expression, which may suggest a role in energy metabolism. If Wolbachia can supplement host energy metabolism, reduced mitochondrial function in infected filarial species may account for Wolbachia-dependence. Wolbachia also have a strong influence on mitochondrial evolution due to vertical co-transmission. This could drive alterations in mitochondrial genome sequence in infected species. Comparisons between the mitochondrial genome sequences of Wolbachia-dependent and independent filarial worms may reveal differences indicative of altered mitochondrial function. Results The mitochondrial genomes of 5 species of filarial nematodes, Acanthocheilonema viteae, Chandlerella quiscali, Loa loa, Onchocerca flexuosa, and Wuchereria bancrofti, were sequenced, annotated and compared with available mitochondrial genome sequences from Brugia malayi, Dirofilaria immitis, Onchocerca volvulus and Setaria digitata. B. malayi, D. immitis, O. volvulus and W. bancrofti are Wolbachia-dependent while A. viteae, C. quiscali, L. loa, O. flexuosa and S. digitata are Wolbachia-free. The 9 mitochondrial genomes were similar in size and AT content and encoded the same 12 protein-coding genes, 22 tRNAs and 2 rRNAs. Synteny was perfectly preserved in all species except C. quiscali, which had a different order for 5 tRNA genes. Protein-coding genes were expressed at the RNA level in all examined species. In phylogenetic trees based on mitochondrial protein-coding sequences, species did not cluster according to Wolbachia dependence. Conclusions Thus far, no discernable differences were detected between the mitochondrial genome sequences of Wolbachia-dependent and independent species. Additional research will be needed to determine whether mitochondria from Wolbachia-dependent filarial species show reduced function in comparison to the mitochondria of Wolbachia-independent species despite their sequence-level similarities.

Most filarial parasites in the subfamilies Onchocercinae and Dirofilariinae depend on Wolbachia endobacteria to successfully carry out their life cycle. Recently published data indicate that the few Wolbachia-free species in these subfamilies were infected in the distant past and have subsequently shed their endosymbionts. We used an integrated transcriptomic and proteomic analysis of Onchocerca flexuosa to explore the molecular mechanisms that allow worms of this species to survive without a bacterial partner. Roche/454 sequencing of the adult transcriptome produced 16,814 isogroup and 47,252 singleton sequences that are estimated to represent approximately 41% of the complete gene set. Sequences similar to 97 Wolbachia genes were identified from the transcriptome, some of which appear on the same transcripts as sequences similar to nematode genes. Computationally predicted peptides, including those with similarity to Wolbachia proteins, were classified at the domain and pathway levels in order to assess the metabolic capabilities of O. flexuosa and compare against the Wolbachia-dependent model filaria, Brugia malayi. Transcript data further facilitated a shotgun proteomic analysis of O. flexuosa adult worm lysate, resulting in the identification of 1,803 proteins. Three of the peptides detected by mass spectroscopy map to two ABC transport-related proteins from Wolbachia. Antibodies raised to one of the Wolbachia-like peptides labeled a single 38 kDa band on Western blots of O. flexuosa lysate and stained specific worm tissues by immunohistology. Future studies will be required to determine the exact functions of Wolbachia-like peptides and proteins in O. flexuosa and to assess their roles in worm biology.

Background: Wolbachia are among the most abundant symbiotic microbes on earth; they are present in about 66% of all insect species, some spiders, mites and crustaceans, and most filarial nematode species. Infected filarial nematodes, including many pathogens of medical and veterinary importance, depend on Wolbachia for proper development and survival. The mechanisms behind this interdependence are not understood. Interestingly, a minority of filarial species examined to date are naturally Wolbachia-free. Methodology/PrincipalFindings:We used 454 pyrosequencing to survey the genomes of two distantly related Wolbachia- free filarial species, Acanthocheilonema viteae and Onchocerca flexuosa. This screen identified 49 Wolbachia-like DNA sequences in A. viteae and 114 in O. flexuosa. qRT-PCR reactions detected expression of 30 Wolbachia-like sequences in A. viteae and 56 in O. flexuosa. Approximately half of these appear to be transcribed from pseudogenes. In situ hybridization showed that two of these pseudogene transcripts were specifically expressed in developing embryos and testes of both species. Conclusions/Significance: These results strongly suggest that the last common ancestor of extant filarial nematodes was infected with Wolbachia and that this former endosymbiont contributed to their genome evolution. Horizontally transferred Wolbachia DNA may explain the ability of some filarial species to live and reproduce without the endosymbiont while other species cannot.

Background Most filarial species in the genus Onchocerca depend on Wolbachia endobacteria to successfully carry out their life cycle. O. flexuosa is a Wolbachia-free species, but its genome contains Wolbachia-like sequences presumably obtained from Wolbachia via horizontal gene transfer. Proteogenomic studies have shown that many of these Wolbachia-like sequences are expressed in adult worms. Methods Six Wolbachia-like sequences in O. flexuosa were chosen for further study based on their sequence conservation with Wolbachia genes, length of predicted open reading frames, and expression at the RNA and/or protein levels. In situ hybridization and immunohistochemical labeling were used to localize Wolbachia-like transcripts and peptides in adult worm tissues. Results RNA probes representing three of the six target sequences produced hybridization signals in worm tissues. These probes bound to transcripts in the intestine and lateral chords of both sexes, in the hypodermis, median chords and uteri in females, and in sperm precursor cells in males. Antibodies raised to three peptides corresponding to these transcripts bound to specific bands in a soluble extract of adult O. flexuosa by Western blot that were not labeled by control antibodies in pre-immune serum. Two of the three antibodies produced labeling patterns in adult worm sections that were similar to those of the RNA probes, while the third produced a different pattern. Conclusions A subset of the Wolbachia-like sequences present in the genome of the Wolbachia-free filarial species O. flexuosa are transcribed in tissues where Wolbachia reside in infected filarial species. Some of the peptides and/or proteins derived from these transcripts appear to be concentrated in the same tissues while others may be exported to other regions of the worm. These results suggest that horizontally transferred Wolbachia genes and gene products may replicate important Wolbachia functions in uninfected filarial worms.