Impact of fetal expression quantitative trait loci on transcriptome-wide association study of childhood leukemia.

Transcriptome-wide association studies increase the yield of loci associated with disease phenotypes by focusing on expression quantitative trait loci (eQTL). The major source of eQTL data for is the Gene and Tissue Expression (GTEx) project, which is comprised entirely of adults, mainly those >50 years of age at death. Since gene expression levels differ by developmental stage, it is not clear whether eQTLs derived from adult data sources are best suited for use in young-onset diseases such as pediatric cancers. To fill in this knowledge gap, we performed a large-scale eQTL mapping analysis in the GenCord study with newborn samples and compared it with GTEx. Under matched conditions, we found around 80% of the eQTLs in one study can be replicated in the other. However, among all eQTLs identified in GenCord (GTEx), 584 (1045) showed statistically significant differences in effect sizes in GTEx (GenCord). We further investigated how using fetal eQTL data can facilitate the genetic association study of acute lymphoblastic leukemia. GenCord and GTEx identified the same genetic loci with statistical significance; however, the overall association pattern was only weakly correlated. Our paper demonstrates age-differential eQTLs and shows their potential influence on childhood leukemia research.

Time-dependent Pax3-mediated chromatin remodeling and cooperation with Six4 and Tead2 specify the skeletal myogenic lineage in developing mesoderm.

The transcriptional mechanisms driving lineage specification during development are still largely unknown, as the interplay of multiple transcription factors makes it difficult to dissect these molecular events. Using a cell-based differentiation platform to probe transcription function, we investigated the role of the key paraxial mesoderm and skeletal myogenic commitment factors-mesogenin 1 (Msgn1), T-box 6 (Tbx6), forkhead box C1 (Foxc1), paired box 3 (Pax3), Paraxis, mesenchyme homeobox 1 (Meox1), sine oculis-related homeobox 1 (Six1), and myogenic factor 5 (Myf5)-in paraxial mesoderm and skeletal myogenesis. From this study, we define a genetic hierarchy, with Pax3 emerging as the gatekeeper between the presomitic mesoderm and the myogenic lineage. By assaying chromatin accessibility, genomic binding and transcription profiling in mesodermal cells from mouse and human Pax3-induced embryonic stem cells and Pax3-null embryonic day (E)9.5 mouse embryos, we identified conserved Pax3 functions in the activation of the skeletal myogenic lineage through modulation of Hedgehog, Notch, and bone morphogenetic protein (BMP) signaling pathways. In addition, we demonstrate that Pax3 molecular function involves chromatin remodeling of its bound elements through an increase in chromatin accessibility and cooperation with sine oculis-related homeobox 4 (Six4) and TEA domain family member 2 (Tead2) factors. To our knowledge, these data provide the first integrated analysis of Pax3 function, demonstrating its ability to remodel chromatin in mesodermal cells from developing embryos and proving a mechanistic footing for the transcriptional hierarchy driving myogenesis.

Reprogramming responsiveness to checkpoint blockade in dysfunctional CD8 T cells.

Established T cell dysfunction is a barrier to antitumor responses, and checkpoint blockade presumably reverses this. Many patients fail to respond to treatment and/or develop autoimmune adverse events. The underlying reason for T cell responsiveness remains elusive. Here, we show that susceptibility to checkpoint blockade is dependent on the activation status of T cells. Newly activated self-specific CD8 T cells respond to checkpoint blockade and cause autoimmunity, which is mitigated by inhibiting the mechanistic target of rapamycin. However, once tolerance is established, self-specific CD8 T cells display a gene signature comparable to tumor-specific CD8 T cells in a fixed state of dysfunction. Tolerant self-specific CD8 T cells do not respond to single or combinatorial dosing of anti-CTLA4, anti-PD-L1, anti-PD-1, anti-LAG-3, and/or anti-TIM-3. Despite this, T cell responsiveness can be induced by vaccination with cognate antigen, which alters the previously fixed transcriptional signature and increases antigen-sensing machinery. Antigenic reeducation of tolerant T cells synergizes with checkpoint blockade to generate functional CD8 T cells, which eliminate tumors without concomitant autoimmunity and are transcriptionally distinct from classic effector T cells. These data demonstrate that responses to checkpoint blockade are dependent on the activation state of a T cell and show that checkpoint blockade-insensitive CD8 T cells can be induced to respond to checkpoint blockade with robust antigenic stimulation to participate in tumor control.

Development of an exosomal gene signature to detect residual disease in dogs with osteosarcoma using a novel xenograft platform and machine learning.

Osteosarcoma has a guarded prognosis. A major hurdle in developing more effective osteosarcoma therapies is the lack of disease-specific biomarkers to predict risk, prognosis, or therapeutic response. Exosomes are secreted extracellular microvesicles emerging as powerful diagnostic tools. However, their clinical application is precluded by challenges in identifying disease-associated cargo from the vastly larger background of normal exosome cargo. We developed a method using canine osteosarcoma in mouse xenografts to distinguish tumor-derived from host-response exosomal messenger RNAs (mRNAs). The model allows for the identification of canine osteosarcoma-specific gene signatures by RNA sequencing and a species-differentiating bioinformatics pipeline. An osteosarcoma-associated signature consisting of five gene transcripts (SKA2, NEU1, PAF1, PSMG2, and NOB1) was validated in dogs with spontaneous osteosarcoma by real-time quantitative reverse transcription PCR (qRT-PCR), while a machine learning model assigned dogs into healthy or disease groups. Serum/plasma exosomes were isolated from 53 dogs in distinct clinical groups ("healthy", "osteosarcoma", "other bone tumor", or "non-neoplastic disease"). Pre-treatment samples from osteosarcoma cases were used as the training set, and a validation set from post-treatment samples was used for testing, classifying as "osteosarcoma detected" or "osteosarcoma-NOT detected". Dogs in a validation set whose post-treatment samples were classified as "osteosarcoma-NOT detected" had longer remissions, up to 15 months after treatment. In conclusion, we identified a gene signature predictive of molecular remissions with potential applications in the early detection and minimal residual disease settings. These results provide proof of concept for our discovery platform and its utilization in future studies to inform cancer risk, diagnosis, prognosis, and therapeutic response.

Extensive libraries of gene truncation variants generated by in vitro transposition.

The detailed analysis of the impact of deletions on proteins or nucleic acids can reveal important functional regions and lead to variants with improved macromolecular properties. We present a method to generate large libraries of mutants with deletions of varying length that are randomly distributed throughout a given gene. This technique facilitates the identification of crucial sequence regions in nucleic acids or proteins. The approach utilizes in vitro transposition to generate 5΄ and 3΄ fragment sub-libraries of a given gene, which are then randomly recombined to yield a final library comprising both terminal and internal deletions. The method is easy to implement and can generate libraries in three to four days. We used this approach to produce a library of >9000 random deletion mutants of an artificial RNA ligase enzyme representing 32% of all possible deletions. The quality of the library was assessed by next-generation sequencing and detailed bioinformatics analysis. Finally, we subjected this library to in vitro selection and obtained fully functional variants with deletions of up to 18 amino acids of the parental enzyme that had been 95 amino acids in length.

Adjusting scoring matrices to correct overextended alignments.

MOTIVATION: Sequence similarity searches performed with BLAST, SSEARCH and FASTA achieve high sensitivity by using scoring matrices (e.g. BLOSUM62) that target low identity (<33%) alignments. Although such scoring matrices can effectively identify distant homologs, they can also produce local alignments that extend beyond the homologous regions. RESULTS: We measured local alignment start/stop boundary accuracy using a set of queries where the correct alignment boundaries were known, and found that 7% of BLASTP and 8% of SSEARCH alignment boundaries were overextended. Overextended alignments include non-homologous sequences; they occur most frequently between sequences that are more closely related (>33% identity). Adjusting the scoring matrix to reflect the identity of the homologous sequence can correct higher identity overextended alignment boundaries. In addition, the scoring matrix that produced a correct alignment could be reliably predicted based on the sequence identity seen in the original BLOSUM62 alignment. Realigning with the predicted scoring matrix corrected 37% of all overextended alignments, resulting in more correct alignments than using BLOSUM62 alone.

Genetic analyses identify evidence for a causal relationship between Ewing sarcoma and hernias.

Knowledge of Ewing sarcoma (EWS) risk factors is exceedingly limited; however, multiple small, independent studies have suggested a possible connection between hernia and EWS. By leveraging hernia summary statistics from the UK Biobank and a recently published genome-wide association study of EWS (733 EWS cases and 1,346 controls), we conducted a genetic investigation of the relationship of 5 hernia types (diaphragmatic, inguinal, umbilical, femoral, and ventral) and EWS. We discovered a positive causal relationship between inguinal hernia and EWS (OR 1.27, 95% confidence interval [CI] 1.01-1.59, and p = 0.041) through Mendelian randomization analysis. Further analyses suggested shared pathways through three genes: HMGA2, LOX, and FBXW7. Diaphragmatic hernia showed a stronger causal relationship with EWS among all of the hernia types (OR 2.26, 95% CI 1.30-3.95, p = 0.004), but no statistically significant local correlation pattern was observed. No evidence of a causal or genetic relationship was observed between EWS and the other three hernia types, including umbilical hernia, despite a previous report indicating an OR as high as 3.3. The finding of our genetic analysis provided additional support to the hypothesis that EWS and hernias may share a common origin.

KMT2A-D pathogenicity, prevalence, and variation according to a population database.

INTRODUCTION: The KMT2 family of genes is essential epigenetic regulators promoting gene expression. The gene family contains three subgroups, each with two paralogues: KMT2A and KMT2B; KMT2C and KMT2D; KMT2F and KMT2G. KMT2A-D are among the most frequent somatically altered genes in several different cancer types. Somatic KMT2A rearrangements are well-characterized in infant leukemia (IL), and growing evidence supports the role of additional family members (KMT2B, KMT2C, and KMT2D) in leukemogenesis. Enrichment of rare heterozygous frameshift variants in KMT2A and C has been reported in acute myeloid leukemia (AML), IL, and solid tumors. Currently, the non-synonymous variation, prevalence, and penetrance of these four genes are unknown. METHODS: This study determined the prevalence of pathogenic/likely pathogenic (P/LP) germline KMT2A-D variants in a cancer-free adult population from the Genome Aggregation Database (gnomAD). Two methods of variant interpretation were utilized: a manual genomic variant interpretation and an automated ACMG pipeline. RESULTS: The ACMG pipeline identified considerably fewer P/LP variants (n = 89) compared to the manual method (n = 660) in all 4 genes. Consequently, the total P/LP prevalence and allele frequency (AF) were higher in the manual method (1:112, AF = 4.46E-03) than in ACMG (1:832, AF = 6.01E-04). Multiple ancestry-exclusive P/LP variants were identified along with an increased frequency in males compared to females. Many of these variants identified in this population database are also associated with severe juvenile conditions. CONCLUSION: These data demonstrate that putatively functional germline variation in these developmentally important genes is more common than previously appreciated and identification in cancer-free adults may indicate incomplete penetrance for many of these variants. Future research should examine a genetic predisposing role in IL and other pediatric cancers.

Genetic ancestry, differential gene expression, and survival in pediatric B-cell acute lymphoblastic leukemia.

BACKGROUND: Black children have lower incidence yet worse survival than White and Latinx children with B-cell acute lymphoblastic leukemia (B-ALL). It is unclear how reported race/ethnicity (RRE) is associated with death in B-ALL after accounting for differentially expressed genes associated with genetic ancestry. METHODS: Using Phase 1 and 2 NCI TARGET B-ALL cases (N = 273; RRE-Black = 21, RRE-White = 162, RRE-Latinx = 69, RRE-Other = 9, RRE-Unknown = 12), we estimated proportions of African (AFR), European (EUR), and Amerindian (AMR) genetic ancestry. We estimated hazard ratios (HR) and 95% confidence intervals (95% CI) between ancestry and death while adjusting for RRE and clinical measures. We identified genes associated with genetic ancestry and adjusted for them in RRE and death associations. RESULTS: Genetic ancestry varied within RRE (RRE-Black, AFR proportion: Mean: 78.5%, Range: 38.2%-93.6%; RRE-White, EUR proportion: Mean: 94%, Range: 1.6%-99.9%; RRE-Latinx, AMR proportion: Mean: 52.0%, Range: 1.2%-98.7%). We identified 10, 1, and 6 differentially expressed genes (padjusted  <0.05) associated with AFR, AMR, and EUR ancestry proportion, respectively. We found AMR and AFR ancestry were statistically significantly associated with death (AMR each 10% HR: 1.05, 95% CI: 1.03-1.17, AFR each 10% increase HR: 1.03, 95% CI:1.01-1.19). RRE differences in the risk of death were larger in magnitude upon adjustment for genes associated with genetic ancestry for RRE-Black, but not RRE-Latinx children (RRE-Black HR: 3.35, 95% CI: 1.31, 8.53; RRE-Latinx HR: 1.47, 0.88-2.45). CONCLUSIONS: Our work highlights B-ALL survival differences by RRE after adjusting for ancestry differentially expressed genes suggesting other factors impacting survival are important.

Distinct mechanisms of PTEN inactivation in dogs and humans highlight convergent molecular events that drive cell division in the pathogenesis of osteosarcoma.

A hallmark of osteosarcoma in both human and canine tumors is somatic fragmentation and rearrangement of chromosome structure which leads to recurrent increases and decreases in DNA copy number. The PTEN gene has been implicated as an important tumor suppressor in osteosarcoma via forward genetic screens. Here, we analyzed copy number changes, promoter methylation and transcriptomes to better understand the role of PTEN in canine and human osteosarcoma. Reduction in PTEN copy number was observed in 23 of 95 (25%) of the canine tumors examined leading to corresponding decreases in PTEN transcript levels from RNA-Seq samples. Unexpectedly, canine tumors with an intact PTEN locus had higher levels of PTEN transcripts than human tumors. This variation in transcript abundance was used to evaluate the role of PTEN in osteosarcoma biology. Decreased PTEN copy number and transcript level was observed in - and likely an important driver of - increases in cell cycle transcripts in four independent canine transcriptional datasets. In human osteosarcoma, homozygous copy number loss was not observed, instead increased methylation of the PTEN promoter was associated with increased cell cycle transcripts. Somatic modification of PTEN, either by homozygous deletion in dogs or by promoter methylation in humans, is clinically relevant to osteosarcoma, because the cell cycle related transcripts are associated with patient outcomes. The PTEN gene is part of a syntenic rearrangement unique to the canine genome, making it susceptible to somatic loss of both copies of distal chromosome 26 which also includes the FAS death receptor. SIGNIFICANCE STATEMENT: PTEN function is abrogated by different mechanisms in canine and human osteosarcoma tumors leading to uncontrolled cell cycling. Somatic loss of this canine specific syntenic region may help explain why the canine genome appears to be uniquely susceptible to osteosarcoma. Syntenic arrangement, in the context of copy number change, may lead to synergistic interactions that in turn modify species specific cancer risk. Comparative models of tumorigenesis may utilize different driver mechanisms.

Sex differences in methylation profiles are apparent in medulloblastoma, particularly among SHH tumors.

Background: Medulloblastoma, the most common malignant pediatric brain tumor, displays marked sex differences in prevalence of the four main molecular subgroups: SHH, WNT, Group 3 and Group 4. Males are more frequently diagnosed with SHH, Group 3 and 4 tumors, which have worse prognoses than WNT tumors. Little is known about sex differences in methylation profiles within subgroups. Methods: Using publicly available methylation data (Illumina HumanMethylation450K array), we compared beta values for males versus females. Differentially methylated positions (DMP) by sex within medulloblastoma subgroups were identified on the autosomes. DMPs were mapped to genes and Reactome pathway analysis was run by subgroup. Kaplan-Meier survival curves (Log-Rank p-values) were assessed for each sex within subgroup. MethylCIBERSORT was used to investigate the tumor microenvironment using deconvolution to estimate the abundances of immune cell types using DNA methylation data. Results: There were statistically significant differences in sex by medulloblastoma subgroups (chi-squared p-value=0.00004): Group 3 (n=144; 65% male), Group 4 (n=326; 67% male), SHH (n=223; 57% male) and WNT (n=70; 41% male). Females had worse survival than males for SHH (p-value=0.02). DMPs by sex were identified within subgroups: SHH (n=131), Group 4 (n=29), Group 3 (n=19), and WNT (n=16) and validated in an independent dataset. Unsupervised hierarchical clustering showed that sex-DMPs in SHH did not correlate with other tumor attributes. Ten genes with sex DMPs (RFTN1, C1orf103, FKBP1B, COL25A1, NPDC1, B3GNT1, FOXN3, RNASEH2C, TLE1, and PHF17) were shared across subgroups. Significant pathways (p<0.05) associated with DMPs were identified for SHH (n=22) and Group 4 (n=4) and included signaling pathways for RET proto-oncogene, advanced glycosylation end product receptor, regulation of KIT, neurotrophic receptors, NOTCH, and TGF-ß. In SHH, we identified DMPs in four genes (CDK6, COL25A1, MMP16, PRIM2) that encode proteins which are the target of therapies in clinical trials for other cancers. There were few sex differences in immune cell composition within tumor subgroups. Conclusion: There are sexually dimorphic methylation profiles for SHH medulloblastoma where survival differences were observed. Sex-specific therapies in medulloblastoma may impact outcomes.

Cas9-induced targeted integration of large DNA payloads in primary human T cells via homology-mediated end-joining DNA repair.

The reliance on viral vectors for the production of genetically engineered immune cells for adoptive cellular therapies remains a translational bottleneck. Here we report a method leveraging the DNA repair pathway homology-mediated end joining, as well as optimized reagent composition and delivery, for the Cas9-induced targeted integration of large DNA payloads into primary human T cells with low toxicity and at efficiencies nearing those of viral vectors (targeted knock-in of 1-6.7 kb payloads at rates of up to 70% at multiple targeted genomic loci and with cell viabilities of over 80%). We used the method to produce T cells with an engineered T-cell receptor or a chimaeric antigen receptor and show that the cells maintained low levels of exhaustion markers and excellent capacities for proliferation and cytokine production and that they elicited potent antitumour cytotoxicity in vitro and in mice. The method is readily adaptable to current good manufacturing practices and scale-up processes, and hence may be used as an alternative to viral vectors for the production of genetically engineered T cells for cancer immunotherapies.

Neoadjuvant systemic oncolytic vesicular stomatitis virus is safe and may enhance long-term survivorship in dogs with naturally occurring osteosarcoma.

Osteosarcoma is a devastating bone cancer that disproportionally afflicts children, adolescents, and young adults. Standard therapy includes surgical tumor resection combined with multiagent chemotherapy, but many patients still suffer from metastatic disease progression. Neoadjuvant systemic oncolytic virus (OV) therapy has the potential to improve clinical outcomes by targeting primary and metastatic tumor sites and inducing durable antitumor immune responses. Here we describe the first evaluation of neoadjuvant systemic therapy with a clinical-stage recombinant oncolytic vesicular stomatitis virus (VSV), VSV-IFNß-NIS, in naturally occurring cancer, specifically appendicular osteosarcoma in companion dogs. Canine osteosarcoma has a similar natural disease history as its human counterpart. VSV-IFNß-NIS was administered prior to standard of care surgical resection, permitting microscopic and genomic analysis of tumors. Treatment was well-tolerated and a "tail" of long-term survivors (∼35%) was apparent in the VSV-treated group, a greater proportion than observed in two contemporary control cohorts. An increase in tumor inflammation was observed in VSV-treated tumors and RNA-seq analysis showed that all the long-term responders had increased expression of a T cell anchored immune gene cluster. We conclude that neoadjuvant VSV-IFNß-NIS is safe and may increase long-term survivorship in dogs with naturally occurring osteosarcoma, particularly those that exhibit pre-existing antitumor immunity.

Neoadjuvant systemic oncolytic vesicular stomatitis virus is safe and may enhance long-term survivorship in dogs with naturally occurring osteosarcoma.

Osteosarcoma is a devastating bone cancer that disproportionally afflicts children, adolescents, and young adults. Standard therapy includes surgical tumor resection combined with multiagent chemotherapy, but many patients still suffer from metastatic disease progression. Neoadjuvant systemic oncolytic virus (OV) therapy has the potential to improve clinical outcomes by targeting primary and metastatic tumor sites and inducing durable antitumor immune responses. Here we described the first evaluation of neoadjuvant systemic therapy with a clinical-stage recombinant oncolytic Vesicular stomatitis virus (VSV), VSV-IFNß-NIS, in naturally occurring cancer, specifically appendicular osteosarcoma in companion dogs. Canine osteosarcoma has a similar natural disease history as its human counterpart. VSV-IFNß-NIS was administered prior to standard of care surgical resection, permitting microscopic and genomic analysis of tumors. Treatment was well-tolerated and a 'tail' of long-term survivors (~35%) was apparent in the VSV-treated group, a greater proportion than observed in two contemporary control cohorts. An increase in tumor inflammation was observed in VSV-treated tumors and RNAseq analysis showed that all the long-term responders had increased expression of a T-cell anchored immune gene cluster. We conclude that neoadjuvant VSV-IFNß-NIS is safe and may increase long-term survivorship in dogs with naturally occurring osteosarcoma, particularly those that exhibit pre-existing antitumor immunity.

Comparative analysis of genome-wide DNA methylation identifies patterns that associate with conserved transcriptional programs in osteosarcoma.

Osteosarcoma is an aggressive tumor of the bone that primarily affects young adults and adolescents. Osteosarcoma is characterized by genomic chaos and heterogeneity. While inactivation of tumor protein p53 (TP53) is nearly universal other high frequency mutations or structural variations have not been identified. Despite this genomic heterogeneity, key conserved transcriptional programs associated with survival have been identified across human, canine and induced murine osteosarcoma. The epigenomic landscape, including DNA methylation, plays a key role in establishing transcriptional programs in all cell types. The role of epigenetic dysregulation has been studied in a variety of cancers but has yet to be explored at scale in osteosarcoma. Here we examined genome-wide DNA methylation patterns in 24 human and 44 canine osteosarcoma samples identifying groups of highly correlated DNA methylation marks in human and canine osteosarcoma samples. We also link specific DNA methylation patterns to key transcriptional programs in both human and canine osteosarcoma. Building on previous work, we built a DNA methylation-based measure for the presence and abundance of various immune cell types in osteosarcoma. Finally, we determined that the underlying state of the tumor, and not changes in cell composition, were the main driver of differences in DNA methylation across the human and canine samples. SIGNIFICANCE: Genome wide comparison of DNA methylation patterns in osteosarcoma across two species lays the ground work for the exploration of DNA methylation programs that help establish conserved transcriptional programs in the context of varied mutational landscapes.

Sex differences in expression of immune elements emerge in children, young adults and mice with osteosarcoma.

BACKGROUND: Males < 40 years old are more likely to be diagnosed with and die from osteosarcoma (OS). The underlying mechanisms may depend on sex differences in immune response. METHODS: We used SEER data to estimate survival differences between males and females aged < 40 years at OS diagnosis. In NCI TARGET-OS cases, we determined sex differences in gene expression, conducted Gene Set Enrichment Analysis (GSEA), and applied the LM22 signature to identify biologic sex differences. We compared sex differences in gene expression profiles in TARGET-OS to those observed in Sleeping Beauty (SB) transposon mutagenesis accelerated Trp53R270H-mutant mouse-OS and healthy adult osteoblasts. RESULTS: Males had worse 17-year overall survival than females (SEER p < 0.0001). From 87 TARGET-OS cases, we observed 1018 genes and 69 pathways that differed significantly by sex (adjusted p < 0.05). Pathway and gene lists overlapped with those from mice (p = 0.03) and healthy osteoblasts (p = 0.017), respectively. Pathways that differed significantly by sex were largely immune-based and included the PD-1/PD-L1 immunotherapy pathway. We observed sex differences in M2 macrophages (LM22; p = 0.056) and M1-M2 macrophage transition (GSEA; p = 0.037) in TARGET-OS. LM22 trends were similar in mice. Twenty-four genes differentially expressed by sex in TARGET-OS had existing cancer therapies. CONCLUSIONS: Sex differences in OS gene expression were similar across species and centered on immune pathways. Identified sex-specific therapeutic targets may improve outcomes in young individuals with OS.

Development of a User-Friendly Pipeline for Mutational Analyses of HIV Using Ultra-Accurate Maximum-Depth Sequencing.

Human immunodeficiency virus type 2 (HIV-2) accumulates fewer mutations during replication than HIV type 1 (HIV-1). Advanced studies of HIV-2 mutagenesis, however, have historically been confounded by high background error rates in traditional next-generation sequencing techniques. In this study, we describe the adaptation of the previously described maximum-depth sequencing (MDS) technique to studies of both HIV-1 and HIV-2 for the ultra-accurate characterization of viral mutagenesis. We also present the development of a user-friendly Galaxy workflow for the bioinformatic analyses of sequencing data generated using the MDS technique, designed to improve replicability and accessibility to molecular virologists. This adapted MDS technique and analysis pipeline were validated by comparisons with previously published analyses of the frequency and spectra of mutations in HIV-1 and HIV-2 and is readily expandable to studies of viral mutation across the genomes of both viruses. Using this novel sequencing pipeline, we observed that the background error rate was reduced 100-fold over standard Illumina error rates, and 10-fold over traditional unique molecular identifier (UMI)-based sequencing. This technical advancement will allow for the exploration of novel and previously unrecognized sources of viral mutagenesis in both HIV-1 and HIV-2, which will expand our understanding of retroviral diversity and evolution.

Corrigendum: Methylation of FKBP5 and SLC6A4 in Relation to Treatment Response to Mindfulness Based Stress Reduction for Posttraumatic Stress Disorder.

[This corrects the article DOI: 10.3389/fpsyt.2018.00418.].

Genome-wide association study accounting for anticholinergic burden to examine cognitive dysfunction in psychotic disorders.

Identifying genetic contributors to cognitive impairments in psychosis-spectrum disorders can advance understanding of disease pathophysiology. Although CNS medications are known to affect cognitive performance, they are often not accounted for in genetic association studies. In this study, we performed a genome-wide association study (GWAS) of global cognitive performance, measured as composite z-scores from the Brief Assessment of Cognition in Schizophrenia (BACS), in persons with psychotic disorders and controls (N = 817; 682 cases and 135 controls) from the Bipolar-Schizophrenia Network on Intermediate Phenotypes (B-SNIP) study. Analyses accounting for anticholinergic exposures from both psychiatric and non-psychiatric medications revealed five significantly associated variants located at the chromosome 3p21.1 locus, with the top SNP rs1076425 in the inter-alpha-trypsin inhibitor heavy chain 1 (ITIH1) gene (P = 3.25×E-9). The inclusion of anticholinergic burden improved association models (P < 0.001) and the number of significant SNPs identified. The effect sizes and direction of effect of the top variants remained consistent when investigating findings within individuals receiving specific antipsychotic drugs and after accounting for antipsychotic dose. These associations were replicated in a separate study sample of untreated first-episode psychosis. The chromosome 3p21.1 locus was previously reported to have association with the risk for psychotic disorders and cognitive performance in healthy individuals. Our findings suggest that this region may be a psychosis risk locus that is associated with cognitive mechanisms. Our data highlight the general point that the inclusion of medication exposure information may improve the detection of gene-cognition associations in psychiatric genetic research.

Nearly Half of TP53 Germline Variants Predicted To Be Pathogenic in Patients With Osteosarcoma Are De Novo: A Report From the Children's Oncology Group.

PURPOSE: To ascertain the prevalence of recurrent de novo variants among 240 pediatric patients with osteosarcoma (OS; age < 20 years) unselected for family history of cancer. METHODS: The identification of de novo variants was implemented in 2 phases. In the first, we identified genes with a rare (minor allele frequency < 0.01) de novo variant in > 1 of the 95 case-parent trios examined by whole-exome sequencing (WES) who passed quality control measures. In phase 2, 145 additional patients with OS were evaluated by targeted sequencing to identify rare de novo variants in genes nominated from phase 1. Recurrent rare variants identified from phase 1 and 2 were verified as either de novo or inherited by Sanger sequencing of affected patients and their parents. Categorical and continuous data were analyzed using Fisher exact test and t tests, respectively. RESULTS: Among 95 case-parent trios who underwent WES, we observed 61 de novo variants in 60 genes among 47 patients, with TP53 identified as the only gene with a pathogenic or likely pathogenic (P/LP) de novo variant in more than one case-parent trio. Among all 240 patients with OS, 13 (5.4%) harbored a P/LP TP53 germline variant, of which 6 (46.2%) were confirmed to be de novo. CONCLUSION: Apart from TP53, we did not observe any other recurrent de novo P/LP variants in the case-parent trios, suggesting that new mutations in other genes are not a frequent cause of pediatric OS. That nearly half of P/LP TP53 variants in our sample were de novo suggests universal screening for germline TP53 P/LP variants among pediatric patients with OS should be considered.