St. Jude Survivorship Portal: sharing and analyzing large clinical and genomic datasets from pediatric cancer survivors.

Childhood cancer survivorship studies generate comprehensive datasets comprising demographic, diagnosis, treatment, outcome, and genomic data from survivors. To broadly share this data, we created the St. Jude Survivorship Portal (https://survivorship.stjude.cloud), the first data portal for sharing, analyzing, and visualizing pediatric cancer survivorship data. Over 1,600 phenotypic variables and 400 million genetic variants from over 7,700 childhood cancer survivors can be explored on this free, open-access portal. Summary statistics of variables are computed on-the-fly and visualized through interactive and customizable charts. Survivor cohorts can be customized and/or divided into groups for comparative analysis. Users can also seamlessly perform cumulative incidence and regression analyses on the stored survivorship data. Using the portal, we explored the ototoxic effects of platinum-based chemotherapy, uncovered a novel association between mental health, age, and limb amputation, and discovered a novel haplotype in MAGI3 strongly associated with cardiomyopathy specifically in survivors of African ancestry.

Diphthamide deficiency promotes association of eEF2 with p53 to induce p21 expression and neural crest defects.

Diphthamide is a modified histidine residue unique for eukaryotic translation elongation factor 2 (eEF2), a key ribosomal protein. Loss of this evolutionarily conserved modification causes developmental defects through unknown mechanisms. In a patient with compound heterozygous mutations in Diphthamide Biosynthesis 1 (DPH1) and impaired eEF2 diphthamide modification, we observe multiple defects in neural crest (NC)-derived tissues. Knockin mice harboring the patient's mutations and Xenopus embryos with Dph1 depleted also display NC defects, which can be attributed to reduced proliferation in the neuroepithelium. DPH1 depletion facilitates dissociation of eEF2 from ribosomes and association with p53 to promote transcription of the cell cycle inhibitor p21, resulting in inhibited proliferation. Knockout of one p21 allele rescues the NC phenotypes in the knockin mice carrying the patient's mutations. These findings uncover an unexpected role for eEF2 as a transcriptional coactivator for p53 to induce p21 expression and NC defects, which is regulated by diphthamide modification.

High-grade glioma in infants and young children is histologically, molecularly, and clinically diverse: Results from the SJYC07 trial and institutional experience.

BACKGROUND: High-grade gliomas (HGG) in young children pose a challenge due to favorable but unpredictable outcomes. While retrospective studies broadened our understanding of tumor biology, prospective data is lacking. METHODS: A cohort of children with histologically diagnosed HGG from the SJYC07 trial was augmented with nonprotocol patients with HGG treated at St. Jude Children's Research Hospital from November 2007 to December 2020. DNA methylome profiling and whole genome, whole exome, and RNA sequencing were performed. These data were integrated with histopathology to yield an integrated diagnosis. Clinical characteristics and preoperative imaging were analyzed. RESULTS: Fifty-six children (0.0-4.4 years) were identified. Integrated analysis split the cohort into four categories: infant-type hemispheric glioma (IHG), HGG, low-grade glioma (LGG), and other-central nervous system (CNS) tumors. IHG was the most prevalent (n = 22), occurred in the youngest patients (median age = 0.4 years), and commonly harbored receptor tyrosine kinase gene fusions (7 ALK, 2 ROS1, 3 NTRK1/2/3, 4 MET). The 5-year event-free (EFS) and overall survival (OS) for IHG was 53.13% (95%CI: 35.52-79.47) and 90.91% (95%CI: 79.66-100.00) vs. 0.0% and 16.67% (95%CI: 2.78-99.74%) for HGG (p = 0.0043, p = 0.00013). EFS and OS were not different between IHG and LGG (p = 0.95, p = 0.43). Imaging review showed IHGs are associated with circumscribed margins (p = 0.0047), hemispheric location (p = 0.0010), and intratumoral hemorrhage (p = 0.0149). CONCLUSIONS: HGG in young children is heterogeneous and best defined by integrating histopathological and molecular features. Patients with IHG have relatively good outcomes, yet they endure significant deficits, making them good candidates for therapy de-escalation and trials of molecular targeted therapy.

Prediction of mammalian virus cross-species transmission based on host proteins.

Most emerging viruses are spilled over from mammals. Understanding the mechanism of virus cross-species transmission and identifying zoonotic viruses before their emergence are critical for the prevention and control of newly emerging viruses. This study systematically investigated the host proteins associated with the cross-species transmission of mammalian viruses based on 1,271 pairs of virus-mammal interactions including 382 viruses from 33 viral families and 73 mammal species from 11 orders. Numerous host proteins were found to contribute to the cross-species transmission of mammalian viruses. Host proteins potentially contributing to virus cross-species transmission are specific to viral families, and few overlaps of such host proteins are observed in different viral families. Based on these host proteins, the random-forest (RF) models were built to predict the cross-species transmission potential of mammalian viruses. Moderate performance was obtained when using all viruses together. However, when modeling by viral family, the performance of the RF models varied much among viral families. In 13 viral families such as Flaviviridae, Retroviridae, and Poxviridae, the AUC of the RF model was greater than 0.8. Finally, the contribution of virus receptors to cross-species transmission was evaluated, and the virus receptor was found to have a minor effect in predicting the cross-species transmission of mammalian viruses. The study deepens our understanding of the mechanism of virus cross-species transmission and provides a framework for predicting the cross-species transmission of mammalian viruses. IMPORTANCE Emerging viruses pose serious threats to humans. Understanding the mechanism of virus cross-species transmission and identifying zoonotic viruses before their emergence are critical for the prevention and control of emerging viruses. This study systematically identified host factors associated with cross-species transmission of mammalian viruses and further built machine-learning models for predicting cross-species transmission of the viruses based on host factors including virus receptors. The study not only deepens our understanding of the mechanism of virus cross-species transmission but also provides a framework for predicting the cross-species transmission of mammalian viruses based on host factors.

Systematic identification and characterization of repeat sequences in African swine fever virus genomes.

African swine fever virus (ASFV) is a large DNA virus that infects domestic pigs with high morbidity and mortality rates. Repeat sequences, which are DNA sequence elements that are repeated more than twice in the genome, play an important role in the ASFV genome. The majority of repeat sequences, however, have not been identified and characterized in a systematic manner. In this study, three types of repeat sequences, including microsatellites, minisatellites and short interspersed nuclear elements (SINEs), were identified in the ASFV genome, and their distribution, structure, function, and evolutionary history were investigated. Most repeat sequences were observed in noncoding regions and at the 5' end of the genome. Noncoding repeat sequences tended to form enhancers, whereas coding repeat sequences had a lower ratio of alpha-helix and beta-sheet and a higher ratio of loop structure and surface amino acids than nonrepeat sequences. In addition, the repeat sequences tended to encode penetrating and antimicrobial peptides. Further analysis of the evolution of repeat sequences revealed that the pan-repeat sequences presented an open state, showing the diversity of repeat sequences. Finally, CpG islands were observed to be negatively correlated with repeat sequence occurrences, suggesting that they may affect the generation of repeat sequences. Overall, this study emphasizes the importance of repeat sequences in ASFVs, and these results can aid in understanding the virus's function and evolution.

Metalloprotease ADAM9 cleaves ephrin-B ligands and differentially regulates Wnt and mTOR signaling downstream of Akt kinase in colorectal cancer cells.

Ephrin-B signaling has been implicated in many normal and pathological processes, including neural crest development and tumor metastasis. We showed previously that proteolysis of ephrin-B ligands by the disintegrin metalloprotease ADAM13 is necessary for canonical Wnt signal activation and neural crest induction in Xenopus, but it was unclear if these mechanisms are conserved in mammals. Here, we report that mammalian ADAM9 cleaves ephrin-B1 and ephrin-B2 and can substitute for Xenopus ADAM13 to induce the neural crest. We found that ADAM9 expression is elevated in human colorectal cancer (CRC) tissues and that knockdown (KD) of ADAM9 inhibits the migration and invasion of SW620 and HCT116 CRC cells by reducing the activity of Akt kinase, which is antagonized by ephrin-Bs. Akt is a signaling node that activates multiple downstream pathways, including the Wnt and mTOR pathways, both of which can promote CRC cell migration/invasion. Surprisingly, we also found that KD of ADAM9 downregulates Wnt signaling but has negligible effects on mTOR signaling in SW620 cells; in contrast, mTOR activity is suppressed while Wnt signaling remains unaffected by ADAM9 KD in HCT116 cells. These results suggest that mammalian ADAM9 cleaves ephrin-Bs to derepress Akt and promote CRC migration and invasion; however, the signaling pathways downstream of Akt are differentially regulated by ADAM9 in different CRC cell lines, reflecting the heterogeneity of CRC cells in responding to manipulations of upstream Akt regulators.

An atlas of human viruses provides new insights into diversity and tissue tropism of human viruses.

MOTIVATION: Viruses continue to threaten human health. Yet, the complete viral species carried by humans and their infection characteristics have not been fully revealed. RESULTS: This study curated an atlas of human viruses from public databases and literature, and built the Human Virus Database (HVD). The HVD contains 1131 virus species of 54 viral families which were more than twice the number of the human-infecting virus species reported in previous studies. These viruses were identified in human samples including 68 human tissues, the excreta and body fluid. The viral diversity in humans was age-dependent with a peak in the infant and a valley in the teenager. The tissue tropism of viruses was found to be associated with several factors including the viral group (DNA, RNA or reverse-transcribing viruses), enveloped or not, viral genome length and GC content, viral receptors and the virus-interacting proteins. Finally, the tissue tropism of DNA viruses was predicted using a random-forest algorithm with a middle performance. Overall, the study not only provides a valuable resource for further studies of human viruses but also deepens our understanding toward the diversity and tissue tropism of human viruses. AVAILABILITY AND IMPLEMENTATION: The HVD is available at http://computationalbiology.cn/humanVirusBase/#/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Whole-genome sequencing identifies I-SceI-mediated transgene integration sites in Xenopus tropicalis snai2:eGFP line.

Transgenesis with the meganuclease I-SceI is a safe and efficient method, but the underlying mechanisms remain unclear due to the lack of information on transgene localization. Using I-SceI, we previously developed a transgenic Xenopus tropicalis line expressing enhanced green fluorescent protein driven by the neural crest-specific snai2 promoter/enhancer, which is a powerful tool for studying neural crest development and craniofacial morphogenesis. Here, we carried out whole-genome shotgun sequencing for the snai2:eGFP embryos to identify the transgene integration sites. With a 19x sequencing coverage, we estimated that 6 copies of the transgene were inserted into the Xenopus tropicalis genome in the hemizygous transgenic embryos. Two transgene integration loci adjacent to each other were identified in a noncoding region on chromosome 1, possibly as a result of duplication after a single transgene insertion. Interestingly, genomic DNA at the boundaries of the transgene integration loci contains short sequences homologous to the I-SceI recognition site, suggesting that the integration was not random but probably mediated by sequence homology. To our knowledge, our work represents the first genome-wide sequencing study on a transgenic organism generated with I-SceI, which is useful for evaluating the potential genetic effects of I-SceI-mediated transgenesis and further understanding the mechanisms underlying this transgenic method.

Identification and characterization of circRNAs encoded by MERS-CoV, SARS-CoV-1 and SARS-CoV-2.

The life-threatening coronaviruses MERS-CoV, SARS-CoV-1 and SARS-CoV-2 (SARS-CoV-1/2) have caused and will continue to cause enormous morbidity and mortality to humans. Virus-encoded noncoding RNAs are poorly understood in coronaviruses. Data mining of viral-infection-related RNA-sequencing data has resulted in the identification of 28 754, 720 and 3437 circRNAs encoded by MERS-CoV, SARS-CoV-1 and SARS-CoV-2, respectively. MERS-CoV exhibits much more prominent ability to encode circRNAs in all genomic regions than those of SARS-CoV-1/2. Viral circRNAs typically exhibit low expression levels. Moreover, majority of the viral circRNAs exhibit expressions only in the late stage of viral infection. Analysis of the competitive interactions of viral circRNAs, human miRNAs and mRNAs in MERS-CoV infections reveals that viral circRNAs up-regulated genes related to mRNA splicing and processing in the early stage of viral infection, and regulated genes involved in diverse functions including cancer, metabolism, autophagy, viral infection in the late stage of viral infection. Similar analysis in SARS-CoV-2 infections reveals that its viral circRNAs down-regulated genes associated with metabolic processes of cholesterol, alcohol, fatty acid and up-regulated genes associated with cellular responses to oxidative stress in the late stage of viral infection. A few genes regulated by viral circRNAs from both MERS-CoV and SARS-CoV-2 were enriched in several biological processes such as response to reactive oxygen and centrosome localization. This study provides the first glimpse into viral circRNAs in three deadly coronaviruses and would serve as a valuable resource for further studies of circRNAs in coronaviruses.

Prokaryotic virus host predictor: a Gaussian model for host prediction of prokaryotic viruses in metagenomics.

BACKGROUND: Viruses are ubiquitous biological entities, estimated to be the largest reservoirs of unexplored genetic diversity on Earth. Full functional characterization and annotation of newly discovered viruses requires tools to enable taxonomic assignment, the range of hosts, and biological properties of the virus. Here we focus on prokaryotic viruses, which include phages and archaeal viruses, and for which identifying the viral host is an essential step in characterizing the virus, as the virus relies on the host for survival. Currently, the method for determining the viral host is either to culture the virus, which is low-throughput, time-consuming, and expensive, or to computationally predict the viral hosts, which needs improvements at both accuracy and usability. Here we develop a Gaussian model to predict hosts for prokaryotic viruses with better performances than previous computational methods. RESULTS: We present here Prokaryotic virus Host Predictor (PHP), a software tool using a Gaussian model, to predict hosts for prokaryotic viruses using the differences of k-mer frequencies between viral and host genomic sequences as features. PHP gave a host prediction accuracy of 34% (genus level) on the VirHostMatcher benchmark dataset and a host prediction accuracy of 35% (genus level) on a new dataset containing 671 viruses and 60,105 prokaryotic genomes. The prediction accuracy exceeded that of two alignment-free methods (VirHostMatcher and WIsH, 28-34%, genus level). PHP also outperformed these two alignment-free methods much (24-38% vs 18-20%, genus level) when predicting hosts for prokaryotic viruses which cannot be predicted by the BLAST-based or the CRISPR-spacer-based methods alone. Requiring a minimal score for making predictions (thresholding) and taking the consensus of the top 30 predictions further improved the host prediction accuracy of PHP. CONCLUSIONS: The Prokaryotic virus Host Predictor software tool provides an intuitive and user-friendly API for the Gaussian model described herein. This work will facilitate the rapid identification of hosts for newly identified prokaryotic viruses in metagenomic studies.

VirusCircBase: a database of virus circular RNAs.

Circular RNAs (circRNAs) are covalently closed long noncoding RNAs critical in diverse cellular activities and multiple human diseases. Several cancer-related viral circRNAs have been identified in double-stranded DNA viruses (dsDNA), yet no systematic study about the viral circRNAs has been reported. Herein, we have performed a systematic survey of 11 924 circRNAs from 23 viral species by computational prediction of viral circRNAs from viral-infection-related RNA sequencing data. Besides the dsDNA viruses, our study has also revealed lots of circRNAs in single-stranded RNA viruses and retro-transcribing viruses, such as the Zika virus, the Influenza A virus, the Zaire ebolavirus, and the Human immunodeficiency virus 1. Most viral circRNAs had reverse complementary sequences or repeated sequences at the flanking sequences of the back-splice sites. Most viral circRNAs only expressed in a specific cell line or tissue in a specific species. Functional enrichment analysis indicated that the viral circRNAs from dsDNA viruses were involved in KEGG pathways associated with cancer. All viral circRNAs presented in the current study were stored and organized in VirusCircBase, which is freely available at http://www.computationalbiology.cn/ViruscircBase/home.html and is the first virus circRNA database. VirusCircBase forms the fundamental atlas for the further exploration and investigation of viral circRNAs in the context of public health.

The RNA helicase DDX3 induces neural crest by promoting AKT activity.

Mutations in the RNA helicase DDX3 have emerged as a frequent cause of intellectual disability in humans. Because many individuals carrying DDX3 mutations have additional defects in craniofacial structures and other tissues containing neural crest (NC)-derived cells, we hypothesized that DDX3 is also important for NC development. Using Xenopus tropicalis as a model, we show that DDX3 is required for normal NC induction and craniofacial morphogenesis by regulating AKT kinase activity. Depletion of DDX3 decreases AKT activity and AKT-dependent inhibitory phosphorylation of GSK3ß, leading to reduced levels of ß-catenin and Snai1: two GSK3ß substrates that are crucial for NC induction. DDX3 function in regulating these downstream signaling events during NC induction is likely mediated by RAC1, a small GTPase whose translation depends on the RNA helicase activity of DDX3. These results suggest an evolutionarily conserved role of DDX3 in NC development by promoting AKT activity, and provide a potential mechanism for the NC-related birth defects displayed by individuals harboring mutations in DDX3 and its downstream effectors in this signaling cascade.

FluPhenotype-a one-stop platform for early warnings of the influenza A virus.

MOTIVATION: Newly emerging influenza viruses keep challenging global public health. To evaluate the potential risk of the viruses, it is critical to rapidly determine the phenotypes of the viruses, including the antigenicity, host, virulence and drug resistance. RESULTS: Here, we built FluPhenotype, a one-stop platform to rapidly determinate the phenotypes of the influenza A viruses. The input of FluPhenotype is the complete or partial genomic/protein sequences of the influenza A viruses. The output presents five types of information about the viruses: (i) sequence annotation including the gene and protein names as well as the open reading frames, (ii) potential hosts and human-adaptation-associated amino acid markers, (iii) antigenic and genetic relationships with the vaccine strains of different HA subtypes, (iv) mammalian virulence-related amino acid markers and (v) drug resistance-related amino acid markers. FluPhenotype will be a useful bioinformatic tool for surveillance and early warnings of the newly emerging influenza A viruses. AVAILABILITY AND IMPLEMENTATION: It is publicly available from: http://www.computationalbiology.cn : 18888/IVEW. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

iTRAQ-based proteomic analysis reveals several key metabolic pathways associated with male sterility in Salvia miltiorrhiza.

Male sterility is a common phenomenon in flowering plants, and it has been widely used in hybrid seed production in a number of economically important crops. In 2002, our team discovered a natural male sterile mutant of Salvia miltiorrhiza. It provided us with the possibility of obtaining stable and controllable quality. To study the molecular mechanism of male sterility in S. miltiorrhiza, we generated proteomic profiles comparing the male sterile mutant type (MT) and wild type (WT) using iTRAQ sequencing. We found a total of 639 differential abundant proteins (DAPs) between MT and WT buds. The DAPs associated with male sterility were mainly involved in (1) carbohydrate and energy metabolism, and (2) protein synthesis and degradation. Based on a comparison between the protein expression profiles of MT and WT, we elucidated a potential protein interaction network involved in male sterility. These results provide new potential biomarkers and insights into the molecular mechanism of male sterility in S. miltiorrhiza.

Identification of a Type-Specific Promoter Element That Differentiates between Influenza A and B Viruses.

Type A and type B influenza viruses (FluA and FluB viruses) are two major human pathogens that share common structural and functional features. FluA and FluB viruses can reassort within each type but never between the types. Here, we bioinformatically analyzed all promoter sequences of FluA and FluB viruses and confirmed the presence of the type-specific promoter elements. We then studied the promoter elements with cell-based in vivo assays and an in vitro replication initiation assay. Our results identified, for the first time, a type-specific promoter element-the nucleotide at position 5 in the 3' end of the viral RNA (vRNA)-that plays a key role(s) in modulating polymerase activity in a type-specific manner. Interestingly, swapping the promoter element between FluA and FluB recombinant viruses showed different tolerances: the replacement of FluA virus-specific U5 with FluB virus-specific C5 in influenza virus A/WSN/33 (H1N1) could be reverted to U5 after 2 to 3 passages, while the replacement of FluB virus-specific C5 with FluA virus-specific U5 in influenza virus B/Yamagata/88 could be maintained, but with significantly reduced replication efficiency. Therefore, our findings indicate that the nucleotide variation at position 5 in the 3' end of the vRNA promoter between FluA and FluB viruses contributes to their RNP incompatibility, which may shed new light on the mechanisms of intertypic exclusion of reassortment between FluA and FluB viruses.IMPORTANCE Genetic reassortment of influenza virus plays a key role in virus evolution and the emergence of pandemic strains. The reassortment occurs extensively within either FluA or FluB viruses but never between them. Here, we bioinformatically compared available promoter sequences of FluA and FluB viruses and confirmed the presence of the type-specific promoter elements. Our in vivo and in vitro mutagenesis studies showed that a type-specific promoter element-the nucleotide at position 5 in the 3' end of vRNA promoters-plays key roles in modulating polymerase activity. Interestingly, FluA and FluB viruses showed different tolerances upon key promoter element swapping in the context of virus infections. We concluded that the nucleotide at position 5 in the 3' end of the vRNA promoters of FluA and FluB viruses is a critical type-specific determinant. This work has implications for further elucidating the mechanisms of the intertypic exclusion of reassortment between FluA and FluB viruses.

Homologous recombination shapes the genetic diversity of African swine fever viruses.

The African swine fever virus (ASFV) has severely influenced the swine industry of the world. Currently, there is no effective vaccine or drugs against the ASFV. How to effectively control the virus is challenging. In this study, we have analyzed all the publicly available ASFV genomes and demonstrated that there was a large genetic diversity of ASFV genomes. Interestingly, the genetic diversity was mainly caused by extensive genomic insertions and/or deletions (indels) instead of the point mutations. Further analyses showed that the indels may be attributed much to the homologous recombination, as supported by significant associations between the occurrence of extensive recombination events and the indels in the ASFV genomes. Besides, the homologous recombination also led to changes of gene content of ASFVs. Finally, repeated elements of dozens of nucleotides in length were observed to widely distribute and cluster in the adjacent positions of ASFV genomes, which may facilitate the occurrence of homologous recombination. This work highlighted the importance of homologous recombination in shaping the genetic diversity of the ASFVs, and could help understand the evolution of the virus.

Rapid identification of human-infecting viruses.

Viruses have caused much mortality and morbidity to humans and pose a serious threat to global public health. The virome with the potential of human infection is still far from complete. Novel viruses have been discovered at an unprecedented pace as the rapid development of viral metagenomics. However, there is still a lack of methodology for rapidly identifying novel viruses with the potential of human infection. This study built several machine learning models to discriminate human-infecting viruses from other viruses based on the frequency of k-mers in the viral genomic sequences. The k-nearest neighbor (KNN) model can predict the human-infecting viruses with an accuracy of over 90%. The performance of this KNN model built on the short contigs (≥1 kb) is comparable to those built on the viral genomes. We used a reported human blood virome to further validate this KNN model with an accuracy of over 80% based on very short raw reads (150 bp). Our work demonstrates a conceptual and generic protocol for the discovery of novel human-infecting viruses in viral metagenomics studies.

[Screening of different AFLP fragments between near-isogenic lines of male sterile and fertile Salvia miltiorrhiza and their comparison analysis].

There is distinctive advantage of using male sterile lines to breed new cultivar and produce hybrids, when compared with general breeding method on yield and quality. In our previous work, near-isogenic lines (NILs) of male sterile and fertile Salvia miltiorrhiza have been obtained through continuous hybridization in many years. In this investigation, 378 primer combination were screened by using AFLP and BSA technique, in which 26 markers amplified from seven primers were found to tightly link to male sterile gene. Based on these markers, two linkage genetic maps were constructed. A 2 027,2 028 bp fragment was amplifed from NILs of fertile and sterile S. miltiorrhiza, respectively, using genome walking technique and previous E11/M4-208 marker as template. Four base mutations were found in intron when comparing both fragments. Among all different markers between NILs of male sterile and fertile S. miltiorrhiza, four was found to have 100% identities to chromosome 1, 3 and 5 of Arabidopsis, namely, E01/M09-418, E05/M13-308, E05/M04-750 and E01/M01-204. The E01/M09-418 marker was very close to male sterile gene of S. miltiorrhiza with distance of 2.1 cM, which also had 100% identities to male sterile gene MS2 in Arabidopsis. Both were distributed in chromosome 3 of Arabidopsis. The 2 028 bp fragment also had 100% identities to MS2 gene. Another E05/M04-750 marker that had 100% identities to chromosome 5 of Arabidopsis was found to have high identities to POP085-M05 gene of poplars and low affinity calcium antiporter CAX2 of Arabidopsis with very low E-value. The constructed genetic map and differential fragments with potential functions found in this study provide a solid foundation to lock male sterile genes in S. miltiorrhiza genome and to discover their functions.