Combining genomics and epidemiology to track mumps virus transmission in the United States.

Unusually large outbreaks of mumps across the United States in 2016 and 2017 raised questions about the extent of mumps circulation and the relationship between these and prior outbreaks. We paired epidemiological data from public health investigations with analysis of mumps virus whole genome sequences from 201 infected individuals, focusing on Massachusetts university communities. Our analysis suggests continuous, undetected circulation of mumps locally and nationally, including multiple independent introductions into Massachusetts and into individual communities. Despite the presence of these multiple mumps virus lineages, the genomic data show that one lineage has dominated in the US since at least 2006. Widespread transmission was surprising given high vaccination rates, but we found no genetic evidence that variants arising during this outbreak contributed to vaccine escape. Viral genomic data allowed us to reconstruct mumps transmission links not evident from epidemiological data or standard single-gene surveillance efforts and also revealed connections between apparently unrelated mumps outbreaks

Fluorescent amplified fragment length polymorphism analysis of Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis isolates

DNA from over 300 Bacillus thuringiensis, Bacillus cereus, and Bacillus anthracis isolates was analyzed by fluorescent amplified fragment length polymorphism (AFLP). B. thuringiensis and B. cereus isolates were from diverse sources and locations, including soil, clinical isolates and food products causing diarrheal and emetic outbreaks, and type strains from the American Type Culture Collection, and over 200 B. thuringiensis isolates representing 36 serovars or subspecies were from the U.S. Department of Agriculture collection. Twenty-four diverse B. anthracis isolates were also included. Phylogenetic analysis of AFLP data revealed extensive diversity within B. thuringiensis and B. cereus compared to the monomorphic nature of B. anthracis. All of the B. anthracis strains were more closely related to each other than to any other Bacillus isolate, while B. cereus and B. thuringiensis strains populated the entire tree. Ten distinct branches were defined, with many branches containing both B. cereus and B. thuringiensis isolates. A single branch contained all the B. anthracis isolates plus an unusual B. thuringiensis isolate that is pathogenic in mice. In contrast, B. thuringiensis subsp. kurstaki (ATCC 33679) and other isolates used to prepare insecticides mapped distal to the B. anthracis isolates. The interspersion of B. cereus and B. thuringiensis isolates within the phylogenetic tree suggests that phenotypic traits used to distinguish between these two species do not reflect the genomic content of the different isolates and that horizontal gene transfer plays an important role in establishing the phenotype of each of these microbes. B. thuringiensis isolates of a particular subspecies tended to cluster together

Risk Factors Associated with Positive QuantiFERON-TB Gold In-Tube and Tuberculin Skin Tests Results in Zambia and South Africa

INTRODUCTION: The utility of T-cell based interferon-gamma release assays for the diagnosis of latent tuberculosis infection remains unclear in settings with a high burden of tuberculosis. OBJECTIVES: To determine risk factors associated with positive QuantiFERON-TB Gold In-Tube (QFT-GIT) and tuberculin skin test (TST) results and the level of agreement between the tests; to explore the hypotheses that positivity in QFT-GIT is more related to recent infection and less affected by HIV than the TST. METHODS: Adult household contacts of tuberculosis patients were invited to participate in a cross-sectional study across 24 communities in Zambia and South Africa. HIV, QFT-GIT and TST tests were done. A questionnaire was used to assess risk factors. RESULTS: A total of 2,220 contacts were seen. 1,803 individuals had interpretable results for both tests, 1,147 (63.6%) were QFT-GIT positive while 725 (40.2%) were TST positive. Agreement between the tests was low (kappa = 0.24). QFT-GIT and TST results were associated with increasing age (adjusted OR [aOR] for each 10 year increase for QFT-GIT 1.15; 95% CI: 1.06-1.25, and for TST aOR: 1.10; 95% CI 1.01-1.20). HIV positivity was less common among those with positive results on QFT-GIT (aOR: 0.51; 95% CI: 0.39-0.67) and TST (aOR: 0.61; 95% CI: 0.46-0.82). Smear positivity of the index case was associated with QFT-GIT (aOR: 1.25; 95% CI: 0.90-1.74) and TST (aOR: 1.39; 95% CI: 0.98-1.98) results. We found little evidence in our data to support our hypotheses. CONCLUSION: QFT-GIT may not be more sensitive than the TST to detect risk factors associated with tuberculous infection. We found little evidence to support the hypotheses that positivity in QFT-GIT is more related to recent infection and less affected by HIV than the TST

Pan-cancer Alterations of the MYC Oncogene and Its Proximal Network across the Cancer Genome Atlas

Although theMYConcogene has been implicated incancer, a systematic assessment of alterations ofMYC, related transcription factors, and co-regulatoryproteins, forming the proximal MYC network (PMN),across human cancers is lacking. Using computa-tional approaches, we define genomic and proteo-mic features associated with MYC and the PMNacross the 33 cancers of The Cancer Genome Atlas.Pan-cancer, 28% of all samples had at least one ofthe MYC paralogs amplified. In contrast, the MYCantagonists MGA and MNT were the most frequentlymutated or deleted members, proposing a roleas tumor suppressors.MYCalterations were mutu-ally exclusive withPIK3CA,PTEN,APC,orBRAFalterations, suggesting that MYC is a distinct onco-genic driver. Expression analysis revealed MYC-associated pathways in tumor subtypes, such asimmune response and growth factor signaling; chro-matin, translation, and DNA replication/repair wereconserved pan-cancer. This analysis reveals insightsinto MYC biology and is a reference for biomarkersand therapeutics for cancers with alterations ofMYC or the PMN

Genomic, Pathway Network, and Immunologic Features Distinguishing Squamous Carcinomas

This integrated, multiplatform PanCancer Atlas study co-mapped and identified distinguishing molecular features of squamous cell carcinomas (SCCs) from five sites associated with smokin

Rare variation at the TNFAIP3 locus and susceptibility to rheumatoid arthritis

Genome-wide association studies (GWAS) conducted using commercial single nucleotide polymorphisms (SNP) arrays have proven to be a powerful tool for the detection of common disease susceptibility variants. However, their utility for the detection of lower frequency variants is yet to be practically investigated. Here we describe the application of a rare variant collapsing method to a large genome-wide SNP dataset, the Wellcome Trust Case Control Consortium rheumatoid arthritis (RA) GWAS. We partitioned the data into gene-centric bins and collapsed genotypes of low frequency variants (defined here as MAF ≤0.05) into a single count coupled with univariate analysis. We then prioritised gene regions for further investigation in an independent cohort of 3,355 cases and 2,427 controls based on rare variant signal p value and prior evidence to support involvement in RA. A total of 14,536 gene bins were investigated in the primary analysis and signals mapping to the TNFAIP3 and chr17q24 loci were selected for further investigation. We detected replicating association to low frequency variants in the TNFAIP3 gene (combined p = 6.6 × 10−6). Even though rare variants are not well-represented and can be difficult to genotype in GWAS, our study supports the application of low frequency variant collapsing methods to genome-wide SNP datasets as a means of exploiting data that are routinely ignored

Pan-Cancer Analysis of lncRNA Regulation Supports Their Targeting of Cancer Genes in Each Tumor Context

Long noncoding RNAs (lncRNAs) are commonly dys-regulated in tumors, but only a handful are known toplay pathophysiological roles in cancer. We inferredlncRNAs that dysregulate cancer pathways, onco-genes, and tumor suppressors (cancer genes) bymodeling their effects on the activity of transcriptionfactors, RNA-binding proteins, and microRNAs in5,185 TCGA tumors and 1,019 ENCODE assays.Our predictions included hundreds of candidateonco- and tumor-suppressor lncRNAs (cancerlncRNAs) whose somatic alterations account for thedysregulation of dozens of cancer genes and path-ways in each of 14 tumor contexts. To demonstrateproof of concept, we showed that perturbations tar-geting OIP5-AS1 (an inferred tumor suppressor) andTUG1 and WT1-AS (inferred onco-lncRNAs) dysre-gulated cancer genes and altered proliferation ofbreast and gynecologic cancer cells. Our analysis in-dicates that, although most lncRNAs are dysregu-lated in a tumor-specific manner, some, includingOIP5-AS1, TUG1, NEAT1, MEG3, and TSIX, synergis-tically dysregulate cancer pathways in multiple tumorcontexts

Spatial Organization and Molecular Correlation of Tumor-Infiltrating Lymphocytes Using Deep Learning on Pathology Images

Beyond sample curation and basic pathologic characterization, the digitized H&E-stained images of TCGA samples remain underutilized. To highlight this resource, we present mappings of tumorinfiltrating lymphocytes (TILs) based on H&E images from 13 TCGA tumor types. These TIL maps are derived through computational staining using a convolutional neural network trained to classify patches of images. Affinity propagation revealed local spatial structure in TIL patterns and correlation with overall survival. TIL map structural patterns were grouped using standard histopathological parameters. These patterns are enriched in particular T cell subpopulations derived from molecular measures. TIL densities and spatial structure were differentially enriched among tumor types, immune subtypes, and tumor molecular subtypes, implying that spatial infiltrate state could reflect particular tumor cell aberration states. Obtaining spatial lymphocytic patterns linked to the rich genomic characterization of TCGA samples demonstrates one use for the TCGA image archives with insights into the tumor-immune microenvironment

Genome-wide linkage screen for testicular germ cell tumour susceptibility loci

A family history of disease is a strong risk factor for testicular germ cell tumour (TGCT). In order to identify the location of putative TGCT susceptibility gene(s) we conducted a linkage search in 237 pedigrees with two or more cases of TGCT. One hundred and seventy-nine pedigrees were evaluated genome-wide with an average inter-marker distance of 10 cM. An additional 58 pedigrees were used to more intensively investigate several genomic regions of interest. Genetic linkage analysis was performed with the ALLEGRO software using two model-based parametric analyses and a non-parametric analysis. Six genomic regions on chromosomes 2p23, 3p12, 3q26, 12p13-q21, 18q21-q23 and Xq27 showed heterogeneity LOD (HLOD) scores of greater than 1, with a maximum HLOD of 1.94 at 3q26. Genome-wide simulation studies indicate that the observed number of HLOD peaks greater than one does not differ significantly from that expected by chance. A TGCT locus at Xq27 has been previously reported. Of the 237 pedigrees examined in this study, 66 were previously unstudied at Xq27, no evidence for linkage to this region was observed in this new pedigree set. Overall, the results indicate that no single major locus can account for the majority of the familial aggregation of TGCT, and suggests that multiple susceptibility loci with weak effects contribute to the diseas

Annual Risk of Tuberculous Infection Using Different Methods in Communities with a High Prevalence of TB and HIV in Zambia and South Africa

BACKGROUND: The annual risk of tuberculous infection (ARTI) is a key epidemiological indicator of the extent of transmission in a community. Several methods have been suggested to estimate the prevalence of tuberculous infection using tuberculin skin test data. This paper explores the implications of using different methods to estimate prevalence of infection and ARTI. The effect of BCG vaccination on these estimates is also investigated. METHODOLOGY/PRINCIPAL FINDINGS: Tuberculin surveys among school children in 16 communities in Zambia and 8 in South Africa (SA) were performed in 2005, as part of baseline data collection and for randomisation purposes of the ZAMSTAR study. Infection prevalence and ARTI estimates were calculated using five methods: different cut-offs with or without adjustments for sensitivity, the mirror method, and mixture analysis. A total of 49,835 children were registered for the surveys, of which 25,048 (50%) had skin tests done and 22,563 (90%) of those tested were read. Infection prevalence was higher in the combined SA than Zambian communities. The mirror method resulted in the least difference of 7.8%, whereas that estimated by the cut-off methods varied from 12.2% to 17.3%. The ARTI in the Zambian and SA communities was between 0.8% and 2.8% and 2.5% and 4.2% respectively, depending on the method used. In the SA communities, the ARTI was higher among the younger children. BCG vaccination had little effect on these estimates. CONCLUSIONS/SIGNIFICANCE: ARTI estimates are dependent on the calculation method used. All methods agreed that there were substantial differences in infection prevalence across the communities, with higher rates in SA. Although TB notification rates have increased over the past decades, the difference in cumulative exposure between younger and older children is less dramatic and a rise in risk of infection in parallel with the estimated incidence of active tuberculosis cannot be excluded