Optimized pipeline of MuTect and GATK tools to improve the detection of somatic single nucleotide polymorphisms in whole-exome sequencing data

Single-cell transcriptome dissecting the microenvironment remodeled by PD1 blockade combined with photodynamic therapy in a mouse model of oral carcinogenesis

Oral squamous cell carcinoma (OSCC) stands as a predominant and perilous malignant neoplasm globally, with the majority of cases originating from oral potential malignant disorders (OPMDs). Despite this, effective strategies to impede the progression of OPMDs to OSCC remain elusive. In this study, we established mouse models of oral carcinogenesis via 4-nitroquinoline 1-oxide induction, mirroring the sequential transformation from normal oral mucosa to OPMDs, culminating in OSCC development. By intervening during the OPMDs stage, we observed that combining PD1 blockade with photodynamic therapy (PDT) significantly mitigated oral carcinogenesis progression. Single-cell transcriptomic sequencing unveiled microenvironmental dysregulation occurring predominantly from OPMDs to OSCC stages, fostering a tumor-promoting milieu characterized by increased Treg proportion, heightened S100A8 expression, and decreased Fib_Igfbp5 (a specific fibroblast subtype) proportion, among others. Notably, intervening with PD1 blockade and PDT during the OPMDs stage hindered the formation of the tumor-promoting microenvironment, resulting in decreased Treg proportion, reduced S100A8 expression, and increased Fib_Igfbp5 proportion. Moreover, combination therapy elicited a more robust treatment-associated immune response compared with monotherapy. In essence, our findings present a novel strategy for curtailing the progression of oral carcinogenesis.

Technical strategy for monozygotic twin discrimination by single-nucleotide variants

Monozygotic (MZ) twins are theoretically genetically identical. Although they are revealed to accumulate mutations after the zygote splits, discriminating between twin genomes remains a formidable challenge in the field of forensic genetics. Single-nucleotide variants (SNVs) are responsible for a substantial portion of genetic variation, thus potentially serving as promising biomarkers for the identification of MZ twins. In this study, we sequenced the whole genome of a pair of female MZ twins when they were 27 and 33 years old to approximately 30 × coverage using peripheral blood on an Illumina NovaSeq 6000 Sequencing System. Potentially discordant SNVs supported by whole-genome sequencing were validated extensively by amplicon-based targeted deep sequencing and Sanger sequencing. In total, we found nine bona fide post-twinning SNVs, all of which were identified in the younger genomes and found in the older genomes. None of the SNVs occurred within coding exons, three of which were observed in introns, supported by whole-exome sequencing results. A double-blind test was employed, and the reliability of MZ twin discrimination by discordant SNVs was endorsed. All SNVs were successfully detected when input DNA amounts decreased to 0.25 ng, and reliable detection was limited to seven SNVs below 0.075 ng input. This comprehensive analysis confirms that SNVs could serve as cost-effective biomarkers for MZ twin discrimination.

Comparative Transcriptomic Analysis of Gossypium hirsutum Fiber Development in Mutant Materials (xin w 139) Provides New Insights into Cotton Fiber Development

Cotton is the most widely planted fiber crop in the world, and improving cotton fiber quality has long been a research hotspot. The development of cotton fibers is a complex process that includes four consecutive and overlapping stages, and although many studies on cotton fiber development have been reported, most of the studies have been based on cultivars that are promoted in production or based on lines that are used in breeding. Here, we report a phenotypic evaluation of Gossypium hirsutum based on immature fiber mutant (xin w 139) and wild-type (Xin W 139) lines and a comparative transcriptomic study at seven time points during fiber development. The results of the two-year study showed that the fiber length, fiber strength, single-boll weight and lint percentage of xin w 139 were significantly lower than those of Xin W 139, and there were no significant differences in the other traits. Principal component analysis (PCA) and cluster analysis of the RNA-sequencing (RNA-seq) data revealed that these seven time points could be clearly divided into three different groups corresponding to the initiation, elongation and secondary cell wall (SCW) synthesis stages of fiber development, and the differences in fiber development between the two lines were mainly due to developmental differences after twenty days post anthesis (DPA). Differential expression analysis revealed a total of 5131 unique differentially expressed genes (DEGs), including 290 transcription factors (TFs), between the 2 lines. These DEGs were divided into five clusters. Each cluster functional category was annotated based on the KEGG database, and different clusters could describe different stages of fiber development. In addition, we constructed a gene regulatory network by weighted correlation network analysis (WGCNA) and identified 15 key genes that determined the differences in fiber development between the 2 lines. We also screened seven candidate genes related to cotton fiber development through comparative sequence analysis and qRT-PCR; these genes included three TFs (GH_A08G1821 (bHLH), GH_D05G3074 (Dof), and GH_D13G0161 (C3H)). These results provide a theoretical basis for obtaining an in-depth understanding of the molecular mechanism of cotton fiber development and provide new genetic resources for cotton fiber research.

Meiosis in an asymmetric dikaryotic genome of Tremella fuciformis Tr01 facilitates new chromosome formation

BACKGROUND

The dikaryotic stage dominates most of the life cycle in basidiomycetes, and each cell carries two different haploid nuclei. Accurate phasing of these two nuclear genomes and their interactions have long been of interest.

RESULTS

We combine PacBio HiFi reads, Nanopore ultra-long reads, and Hi-C data to generate a complete, high-quality asymmetric dikaryotic genome of Tremella fuciformis Tr01, including Haplotypes A and B genomes. We assemble a meiotic haploid DBZ04 genome and detect three recombination events in these two haplotypes. We identify several chromosomal rearrangements that lead to differences in chromosome number, length, content, and sequence arrangement between these two haplotypes. Each nucleus contains a two-speed genome, harboring three accessory chromosomes and two accessory compartments that affect horizontal chromatin transfer between nuclei. We find few basidiospores are ejected from fruiting bodies of Tr01. Most monospore isolates sequenced belong to Tr01-Haplotype A genome architecture. More than one-third of monospore isolates carry one or two extra chromosomes including Chr12B and two new chromosomes ChrN1 and ChrN2. We hypothesize that homologous regions of seven sister chromatids pair into a large complex during meiosis, followed by inter-chromosomal recombination at physical contact sites and formation of new chromosomes.

CONCLUSION

We assemble two haplotype genomes of T. fuciformis Tr01 and provide the first overview of basidiomycetous genomes with discrete genomic architecture. Meiotic activities of asymmetric dikaryotic genomes result in formation of new chromosomes, aneuploidy of some daughter cells, and inviability of most other daughter cells. We propose a new approach for breeding of sporeless mushroom.

Transcriptome profile reveals novel candidate genes associated with bone strength in end-of-lay hens

Bone weakness causes many problems such as osteoporosis, bone fractures, and economic loss, especially at the late stage of lay, in laying hen production. However, the genetic factors and molecular mechanism affecting the bone strength is still largely unknown. To elucidate the molecular mechanism and genetic factors affecting bone strength, a total of six cDNA libraries were constructed and used to compare genetic differences between tibia with higher（Group HBS）and lower（Group LBS）breaking strength in Hyline grey layers. A comparison between Groups HBS and LBS revealed nine differentially expressed genes, of which five were upregulated and four were downregulated in the LBS relative to the HBS in tibia. Our results showed novel candidate genes concerned with bone strength in the late laying period. These include transcription factor paired box protein Pax-5 (Pax5), tissue inhibitor of Metallopoteinase-4 (TIMP4), Kelch-like protein 14 (KLHL14), predicted MAGUK p55 subfamily member 7 isoform X4 (MPP7) and Osteoclast-associated Ig-like receptor (OSCAR). Our data provide a vital resource for discovering important candidate genes associated with bone strength and will help further study the molecular mechanisms for bone remodeling.

Whole-exome and targeted gene sequencing of large-cell lung carcinoma reveals recurrent mutations in the PI3K pathway

BACKGROUND

Large cell lung carcinoma (LCLC) is an exceptionally aggressive disease with a poor prognosis. At present, little is known about the molecular pathology of LCLC.

METHODS

Ultra-deep sequencing of cancer-related genes and exome sequencing were used to detect the LCLC mutational in 118 tumor-normal pairs. The cell function test was employed to confirm the potential carcinogenic mutation of PI3K pathway.

RESULTS

The mutation pattern is determined by the predominance of A > C mutations. Genes with a significant non-silent mutation frequency (FDR) < 0.05) include TP53 (47.5%), EGFR (13.6%) and PTEN (12.1%). Moreover, PI3K signaling (including EGFR, FGRG4, ITGA1, ITGA5, and ITGA2B) is the most mutated pathway, influencing 61.9% (73/118) of the LCLC samples. The cell function test confirmed that the potential carcinogenic mutation of PI3K pathway had a more malignant cell function phenotype. Multivariate analysis further revealed that patients with the PI3K signaling pathway mutations have a poor prognosis (P = 0.007).

CONCLUSIONS

These results initially identified frequent mutation of PI3K signaling pathways in LCLC and indicate potential targets for the treatment of this fatal type of LCLC.

Imaging-Based Screening Identifies Modulators of the eIF3 Translation Initiation Factor Complex in Candida albicans

Fungal pathogens like Candida albicans can cause devastating human disease. Treatment of candidemia is complicated by the high rate of resistance to common antifungal therapies. Additionally, there is host toxicity associated with many antifungal compounds due to the conservation between essential mammalian and fungal proteins. An attractive new approach for antimicrobial development is to target virulence factors: non-essential processes that are required for the organism to cause disease in human hosts. This approach expands the potential target space while reducing the selective pressure toward resistance, as these targets are not essential for viability. In C. albicans, a key virulence factor is the ability to transition to hyphal morphology. We developed a high-throughput image analysis pipeline to distinguish between yeast and filamentous growth in C. albicans at the single cell level. Based on this phenotypic assay, we screened the FDA drug repurposing library of 2,017 compounds for their ability to inhibit filamentation and identified 33 compounds that block the hyphal transition in C. albicans with IC50 values ranging from 0.2 to 150 μM. Multiple compounds showed a phenyl sulfone chemotype, prompting further analysis. Of these phenyl sulfones, NSC 697923 displayed the most efficacy, and by selecting for resistant mutants, we identified eIF3 as the target of NSC 697923 in C. albicans.

A pilot study of lymphodepletion intensity for peripheral blood mononuclear cell-derived neoantigen-specific CD8 + T cell therapy in patients with advanced solid tumors

Currently, the optimal lymphodepletion intensity for peripheral blood mononuclear cell-derived neoantigen-specific CD8 + T cell (Neo-T) therapy has yet to be determined. We report a single-arm, open-label and non-randomized phase 1 study (NCT02959905) of Neo-T therapy with lymphodepletion at various dose intensity in patients with locally advanced or metastatic solid tumors that are refractory to standard therapies. The primary end point is safety and the secondary end points are disease control rate (DCR), progression-free survival (PFS), overall survival (OS). Results show that the treatment is well tolerated with lymphopenia being the most common adverse event in the highest-intensity lymphodepletion groups. Neo-T infusion-related adverse events are only grade 1-2 in the no lymphodepletion group. The median PFS is 7.1 months (95% CI:3.7-9.8), the median OS is 16.8 months (95% CI: 11.9-31.7), and the DCR is 66.7% (6/9) among all groups. Three patients achieve partial response, two of them are in the no lymphodepletion group. In the group without lymphodepletion pretreatment, one patient refractory to prior anti-PD1 therapy shows partial response to Neo-T therapy. Neoantigen specific TCRs are examined in two patients and show delayed expansion after lymphodepletion treatment. In summary, Neo-T therapy without lymphodepletion could be a safe and promising regimen for advanced solid tumors.

Whole-exome sequencing for screening noise-induced hearing loss susceptibility genes

BACKGROUND

High-throughput sequencing of genes indicating susceptibility to noise-induced hearing loss has not previously been reported.

AIMS/OBJECTIVES

To identify and analyze genes associated with susceptibility to noise-induced hearing loss (NIHL) and characterize differences in susceptibility to hearing loss by genotype.

MATERIAL AND METHODS

Pure tone audiometry tests were performed on 113 workers exposed to high-intensity noise. Whole-exome sequencing (WES) was conducted and NIHL susceptibility genes screened for training unsupervised and supervised machine learning models. Immunofluorescence staining of mouse cochlea was used to observe patterns of NIHL susceptibility gene expression.

RESULTS

Participants were divided into a NIHL and a control group, according to the results of audiometry tests. Seventy-three possible NIHL susceptibility genes were input into the machine learning model. Two subgroups of NIHL could be distinguished by unsupervised machine learning and the classification was evaluated by the supervised machine learning algorithm. The VWF gene had the highest mutation frequency in the NIHL group and was expressed mainly in the spiral ligament.

CONCLUSIONS AND SIGNIFICANCE

NIHL susceptibility genes were screened and NIHL subgroups could be distinguished. VWF may be a novel NIHL susceptibility gene.

Genome assembly, resequencing and genome-wide association analyses provide novel insights into the origin, evolution and flower colour variations of flowering cherry

Flowering cherry is a very popular species around the world. High-quality genome resources for different elite cultivars are needed, and the understanding of their origins and the regulation of key ornamental traits are limited for this tree. Here, a high-quality chromosome-scale genome of Prunus campanulata 'Plena' (PCP), which is a native and elite flowering cherry cultivar in China, was generated. The contig N50 of the genome was 18.31 Mb, and 99.98% of its contigs were anchored to eight chromosomes. Furthermore, a total of 306 accessions of flowering cherry germplasm and six lines of outgroups were collected. Resequencing of these 312 lines was performed, and 761 267 high-quality genomic variants were obtained. The origins of flowering cherry were predicted, and these 306 accessions could be classified into three clades, A, B and C. According to phylogenetic analysis, we predicted two origins of flowering cherry. Flowering cherry in clade A originated in southern China, such as in the Himalayan Mountains, while clades B and C originated in northeastern China. Finally, a genome-wide association study of flower colour was performed for all 312 accessions of flowering cherry germplasm. A total of seven quantitative trait loci (QTLs) were identified. One gene encoding glycosylate transferase was predicted as the candidate gene for one QTL. Taken together, our results provide a valuable genomic resource and novel insights into the origin, evolution and flower colour variations of flowering cherry.

satmut_utils: a simulation and variant calling package for multiplexed assays of variant effect

The impact of millions of individual genetic variants on molecular phenotypes in coding sequences remains unknown. Multiplexed assays of variant effect (MAVEs) are scalable methods to annotate relevant variants, but existing software lacks standardization, requires cumbersome configuration, and does not scale to large targets. We present satmut_utils as a flexible solution for simulation and variant quantification. We then benchmark MAVE software using simulated and real MAVE data. We finally determine mRNA abundance for thousands of cystathionine beta-synthase variants using two experimental methods. The satmut_utils package enables high-performance analysis of MAVEs and reveals the capability of variants to alter mRNA abundance.

Imaging-based screening identifies modulators of the eIF3 translation initiation factor complex in Candida albicans

Fungal pathogens like Candida albicans can cause devastating human disease. Treatment of candidemia is complicated by the high rate of resistance to common antifungal therapies. Additionally, there is host toxicity associated with many antifungal compounds due to the conservation between essential mammalian and fungal proteins. An attractive new approach for antimicrobial development is to target virulence factors: non-essential processes that are required for the organism to cause disease in human hosts. This approach expands the potential target space while reducing the selective pressure towards resistance, as these targets are not essential for viability. In C. albicans, a key virulence factor is the ability to transition to hyphal morphology. We developed a high-throughput image analysis pipeline to distinguish between yeast and filamentous growth in C. albicans at the single cell level. Based on this phenotypic assay, we screened the FDA drug repurposing library of 2,017 compounds for their ability to inhibit filamentation and identified 33 compounds that block the hyphal transition in C. albicans with IC 50 values ranging from 0.2 to 150 µM. Multiple compounds showed a phenyl vinyl sulfone chemotype, prompting further analysis. Of these phenyl vinyl sulfones, NSC 697923 displayed the most efficacy, and by selecting for resistant mutants, we identified eIF3 as the target of NSC 697923 in C. albicans .

Afro-TB dataset as a large scale genomic data of Mycobacterium tuberuclosis in Africa

Mycobacterium tuberculosis (MTB) is a pathogenic bacterium accountable for 10.6 million new infections with tuberculosis (TB) in 2021. The fact that the genetic sequences of M. tuberculosis vary widely provides a basis for understanding how this bacterium causes disease, how the immune system responds to it, how it has evolved over time, and how it is distributed geographically. However, despite extensive research efforts, the evolution and transmission of MTB in Africa remain poorly understood. In this study, we used 17,641 strains from 26 countries to create the first curated African Mycobacterium tuberculosis (MTB) classification and resistance dataset, containing 13,753 strains. We identified 157 mutations in 12 genes associated with resistance and additional new mutations potentially associated with resistance. The resistance profile was used to classify strains. We also performed a phylogenetic classification of each isolate and prepared the data in a format that can be used for phylogenetic and comparative analysis of tuberculosis worldwide. These genomic data will extend current information for comparative genomic studies to understand the mechanisms and evolution of MTB drug resistance.

PIK3CA is recurrently mutated in canine mammary tumors, similarly to in human mammary neoplasia

Biological features of neoplastic disease affecting mammary gland tissue are shared between canines and humans. Research performed in either species has translational value and early phase clinical trials performed in canines with spontaneous disease could be informative for human trials. The purpose of this study was to investigate the somatic genetic aberrations occurring in canine mammary neoplasia by exome capture and next generation sequencing. Based on 55 tumor-normal pairs we identified the PIK3CA gene as the most commonly mutated gene in canine mammary tumors, with 25% of samples carrying mutations in this gene. A recurrent missense mutation was identified, p.H1047R, which is homologous to the human PIK3CA hotspot mutation found in different types of breast neoplasia. Mutations homologous to other known human mutation hotspots such as the PIK3CA p.E545K and the KRAS p.G12V/D were also identified. We identified copy number aberrations affecting important tumor suppressor and oncogenic pathways including deletions affecting the PTEN tumor suppressor gene. We suggest that activation of the KRAS or PIK3CA oncogenes or loss of the PTEN suppressor gene may be important for mammary tumor development in dogs. This data endorses the conservation of cancer across species and the validity of studying cancer in non-human species.

SNP-based bulk segregant analysis revealed disease resistance QTLs associated with northern corn leaf blight in maize

Maize (Zea mays L.) is the most important food security crop worldwide. Northern corn leaf blight (NCLB), caused by Exserohilum turcicum, severely reduces production causing millions of dollars in losses worldwide. Therefore, this study aimed to identify significant QTLs associated with NCLB by utilizing next-generation sequencing-based bulked-segregant analysis (BSA). Parental lines GML71 (resistant) and Gui A10341 (susceptible) were used to develop segregating population F₂. Two bulks with 30 plants each were further selected from the segregating population for sequencing along with the parental lines. High throughput sequencing data was used for BSA. We identified 10 QTLs on Chr 1, Chr 2, Chr 3, and Chr 5 with 265 non-synonymous SNPs. Moreover, based on annotation information, we identified 27 candidate genes in the QTL regions. The candidate genes associated with disease resistance include AATP1, At4g24790, STICHEL-like 2, BI O 3-BIO1, ZAR1, SECA2, ABCG25, LECRK54, MKK7, MKK9, RLK902, and DEAD-box ATP-dependent RNA helicase. The annotation information suggested their involvement in disease resistance-related pathways, including protein phosphorylation, cytoplasmic vesicle, protein serine/threonine kinase activity, and ATP binding pathways. Our study provides a substantial addition to the available information regarding QTLs associated with NCLB, and further functional verification of identified candidate genes can broaden the scope of understanding the NCLB resistance mechanism in maize.

Aberrant Cholesterol Metabolism and Wnt/β-Catenin Signaling Coalesce via Frizzled5 in Supporting Cancer Growth

Frizzled (Fzd) proteins are Wnt receptors and play essential roles in development, homeostasis, and oncogenesis. How Wnt/Fzd signaling is coupled to physiological regulation remains unknown. Cholesterol is reported as a signaling molecule regulating morphogen such as Hedgehog signaling. Despite the elusiveness of the in-depth mechanism, it is well-established that pancreatic cancer specially requires abnormal cholesterol metabolism levels for growth. In this study, it is unexpectedly found that among ten Fzds, Fzd5 has a unique capacity to bind cholesterol specifically through its conserved extracellular linker region. Cholesterol-binding enables Fzd5 palmitoylation, which is indispensable for receptor maturation and trafficking to the plasma membrane. In Wnt-addicted pancreatic ductal adenocarcinoma (PDAC), cholesterol stimulates tumor growth via Fzd5-mediated Wnt/β-catenin signaling. A natural oxysterol, 25-hydroxylsterol competes with cholesterol and inhibits Fzd5 maturation and Wnt signaling, thereby alleviating PDAC growth. This cholesterol-receptor interaction and ensuing receptor lipidation uncover a novel mechanism by which Fzd5 acts as a cholesterol sensor and pivotal connection coupling lipid metabolism to morphogen signaling. These findings further suggest that cholesterol-targeting may provide new therapeutic opportunities for treating Wnt-dependent cancers.

Recurrent noncoding somatic and germline WT1 variants converge to disrupt MYB binding in acute promyelocytic leukemia

Genetic alternations can occur at noncoding regions, but how they contribute to cancer pathogenesis is poorly understood. Here, we established a mutational landscape of cis-regulatory regions (CREs) in acute promyelocytic leukemia (APL) based on whole-genome sequencing analysis of paired tumor and germline samples from 24 patients and epigenetic profiling of 16 patients. Mutations occurring in CREs occur preferentially in active enhancers bound by the complex of master transcription factors in APL. Among significantly enriched mutated CREs, we found a recurrently mutated region located within the third intron of WT1, an essential regulator of normal and malignant hematopoiesis. Focusing on noncoding mutations within this WT1 intron, an analysis on 169 APL patients revealed that somatic mutations were clustered into a focal hotspot region, including one site identified as a germline polymorphism contributing to APL risk. Significantly decreased WT1 expression was observed in APL patients bearing somatic and/or germline noncoding WT1 variants. Furthermore, biallelic WT1 inactivation was recurrently found in APL patients with noncoding WT1 variants, which resulted in the complete loss of WT1. The high incidence of biallelic inactivation suggested the tumor suppressor activity of WT1 in APL. Mechanistically, noncoding WT1 variants disrupted MYB binding on chromatin and suppressed the enhancer activity and WT1 expression through destroying the chromatin looping formation. Our study highlights the important role of noncoding variants in the leukemogenesis of APL.

Pulmonate slug evolution is reflected in the de novo genome of Arion vulgaris Moquin-Tandon, 1855

Stylommatophoran pulmonate land slugs and snails successfully completed the water-to-land transition from an aquatic ancestor and flourished on land. Of the 30,000 estimated species, very few genomes have so far been published. Here, we assembled and characterized a chromosome-level genome of the “Spanish” slug, Arion vulgaris Moquin-Tandon, 1855, a notorious pest land slug in Europe. Using this reference genome, we conclude that a whole-genome duplication event occurred approximately 93–109 Mya at the base of Stylommatophora and might have promoted land invasion and adaptive radiation. Comparative genomic analyses reveal that genes related to the development of kidney, blood vessels, muscle, and nervous systems had expanded in the last common ancestor of land pulmonates, likely an evolutionary response to the terrestrial challenges of gravity and water loss. Analyses of A. vulgaris gene families and positively selected genes show the slug has evolved a stronger ability to counteract the greater threats of external damage, radiation, and water loss lacking a protective shell. Furthermore, a recent burst of long interspersed elements in the genome of A. vulgaris might affect gene regulation and contribute to rapid phenotype changes in A. vulgaris, which might be conducive to its rapid adaptation and invasiveness.

Improvements in Quality Control and Library Preparation for Targeted Sequencing Allowed Detection of Potentially Pathogenic Alterations in Circulating Cell-Free DNA Derived from Plasma of Brain Tumor Patients

Malignant gliomas are the most frequent primary brain tumors in adults. They are genetically heterogenous and invariably recur due to incomplete surgery and therapy resistance. Circulating tumor DNA (ctDNA) is a component of circulating cell-free DNA (ccfDNA) and represents genetic material that originates from the primary tumor or metastasis. Brain tumors are frequently located in the eloquent brain regions, which makes biopsy difficult or impossible due to severe postoperative complications. The analysis of ccfDNA from a patient's blood presents a plausible and noninvasive alternative. In this study, freshly frozen tumors and corresponding blood samples were collected from 84 brain tumor patients and analyzed by targeted next-generation sequencing (NGS). The cohort included 80 glioma patients, 2 metastatic cancer patients, and 2 primary CNS lymphoma (PCNSL) patients. We compared the pattern of genetic alterations in the tumor DNA (tDNA) with that of ccfDNA. The implemented technical improvements in quality control and library preparation allowed for the detection of ctDNA in 8 out of 84 patients, including 5 out of 80 glioma patients. In 32 out of 84 patients, we found potentially pathogenic genetic alterations in ccfDNA that were not detectable in tDNA. While sequencing ccfDNA from plasma has a low efficacy as a diagnostic tool for glioma patients, we concluded that further improvements in sample processing and library preparation can make liquid biopsy a valuable diagnostic tool for glioma patients.

Prevalence and Molecular Characteristics of Polymyxin-Resistant Pseudomonas aeruginosa in a Chinese Tertiary Teaching Hospital

Polymyxin-resistant Pseudomonas aeruginosa is a major threat to public health globally. We investigated the prevalence of polymyxin-resistant P. aeruginosa in a Chinese teaching hospital and determined the genetic and drug-resistant phenotypes of the resistant isolates. P. aeruginosa isolates identified by MALDI-TOF MS were collected across a 3-month period in Ruijin Hospital. Antimicrobial susceptibility was determined by a Vitek-2 Compact system with broth dilution used to determine polymyxin B (PMB) susceptibility. Polymyxin-resistant isolates were further characterized by molecular typing using PCR, multi-locus sequence typing (MLST) and whole-genome sequencing. Phylogenetic relationships were analyzed using single nucleotide polymorphism (SNP) from the whole-genome sequencing. Of 362 P. aeruginosa isolates collected, 8 (2.2%) isolates from separate patients across six wards were polymyxin-resistant (MIC range, PMB 4–16 μg/mL and colistin 4–≥16 μg/mL). Four patients received PMB treatments (intravenous, aerosolized and/or topical) and all patients survived to discharge. All polymyxin-resistant isolates were genetically related and were assigned to five different clades (Isolate 150 and Isolate 211 being the same ST823 type). Genetic variations V51I, Y345H, G68S and R155H in pmrB and L71R in pmrA were identified, which might confer polymyxin resistance in these isolates. Six of the polymyxin-resistant isolates showed reduced susceptibility to imipenem and meropenem (MIC range ≥ 16 μg/mL), while two of the eight isolates were resistant to ceftazidime. We revealed a low prevalence of polymyxin-resistant P. aeruginosa in a Chinese teaching hospital with most polymyxin-resistant isolates being multidrug-resistant. Therefore, effective infection control measures are urgently needed to prevent further spread of resistance to the last-line polymyxins.

Prediction of risk-associated genes and high-risk liver cancer patients from their mutation profile: Benchmarking of mutation calling techniques

Identification of somatic mutations with high precision is one of the major challenges in the prediction of high-risk liver-cancer patients. In the past, number of mutations calling techniques have been developed that include MuTect2, MuSE, Varscan2, and SomaticSniper. In this study, an attempt has been made to benchmark the potential of these techniques in predicting the prognostic biomarkers for liver cancer. Initially, we extracted somatic mutations in liver cancer patients using Variant Call Format (VCF) and Mutation Annotation Format (MAF) files from the cancer genome atlas. In terms of size, the MAF files are 42 times smaller than VCF files and containing only high-quality somatic mutations. Further, machine learning based models have been developed for predicting high-risk cancer patients using mutations obtained from different techniques. The performance of different techniques and data files have been compared based on their potential to discriminate high and low-risk liver-cancer patients. Based on correlation analysis, we selected 80 genes having significant negative-correlation with the overall survival of liver cancer patients. The univariate survival analysis revealed the prognostic role of highly mutated genes. Single-gene based analysis showed that MuTect2 technique based MAF file has achieved maximum hazard ratio (HRLAMC3) of 9.25 with p-value 1.78E-06. Further, we developed various prediction models using risk-associated top-10 genes for each technique. Our results indicate that MuTect2 technique based VCF files outperform all other methods with maximum Area Under the Receiver-Operating Characteristic (AUROC) curve of 0.765 and HR 4.50 (p-value 3.83E-15). Eventually, VCF file generated using MuTect2 technique performs better among other mutation calling techniques for the prediction of high-risk liver cancer patients. We hope that our findings will provide a useful and comprehensive comparison of various mutation calling techniques for the prognostic analysis of cancer patients. In order to serve the scientific community, we have provided a Python-based pipeline to develop the prediction models using mutation profiles (VCF/MAF) of cancer patients. It is available on GitHub at https://github.com/raghavagps/mutation_bench.

Pitaya Genome and Multiomics Database (PGMD): A Comprehensive and Integrative Resource of Selenicereus undatus

Pitaya (Selenicereus) is a kind of novel fruit with a delicious taste and superior horticulture ornamental value. The potential economic impact of the pitaya lies in its diverse uses not only as agricultural produce and processed foods but also in industrial and medicinal products. It is also an excellent plant material for basic and applied biological research. A comprehensive database of pitaya would facilitate studies of pitaya and the other Cactaceae plant species. Here, we constructed pitaya genome and multiomics database, which is a collection of the most updated and high-quality pitaya genomic assemblies. The database contains various information such as genomic variation, gene expression, miRNA profiles, metabolite and proteomic data from various tissues and fruit developmental stages of different pitaya cultivars. In PGMD, we also uploaded videos on the flowering process and planting tutorials for practical usage of pitaya. Overall, these valuable data provided in the PGMD will significantly facilitate future studies on population genetics, molecular breeding and function research of pitaya.

Co-segregation of variant NSUN2 Lue198Arg among Iranian family with intellectual disability: a case report

Background

Intellectual disability is characterized by impairments in adaptive behavior and cognitive functioning manifested during the developmental period. Since disabilities are heterogeneous, variant analysis can help us confirm and accurately diagnose children with intellectual disabilities. Some papers reported that bi-allelic variants of the NSUN2 gene caused a group of neurological disorders, including non-syndromic autosomal recessive intellectual disability (NS-ARID), Dubowitz syndrome, and familial restrictive cardiomyopathy 1 (RCM1). We report on a consanguineous family with three siblings diagnosed with intellectual disability.

Case presentation

The 7-year-old female was referred to Ali-Asghar hospital, Zahedan, Iran, with clinical manifestations comprising moderate intellectual disability, ptosis, long face, and short stature. Chromosome banding, metabolic testing, and magnetic resonance imaging examinations revealed no abnormalities. Accordingly, other affected siblings born of the same parents were considered. Whole-exome sequencing (WES) was conducted on the sufferer to consider NS-ARID variants. Findings identified a variant with uncertain significance (NM_017755.6: c.593 T > G) in the NSUN2 gene in the proband. This variant was confirmed through Sanger sequencing of the affected and unaffected family members. Besides, the computational results showed that the L198R exchange could change the interaction between wild-type and other residues in the protein. The affected patients with NS-ARID had similar clinical characteristics and genetic abnormalities.

Conclusion

Taken together, we described the variant in three Iranian siblings; further expanding of the other variants involved in the disease will be evident by using high-throughput sequencing technologies.

Development of a 5-Gene Signature to Evaluate Lung Adenocarcinoma Prognosis Based on the Features of Cancer Stem Cells

Cancer stem cells (CSCs) can induce recurrence and chemotherapy resistance of lung adenocarcinoma (LUAD). Reliable markers identified based on CSC characteristic of LUAD may improve patients' chemotherapy response and prognosis. OCLR was used to calculate mRNA expression-based stemness index (mRNAsi) of LUAD patients' data in TCGA. Association analysis of mRNAsi was performed with clinical features, somatic mutation, and tumor immunity. A prognostic prediction model was established with LASSO Cox regression. Kaplan-Meier Plotter (KM-plotter) and time-dependent ROC were applied to assess signature performance. For LUAD, univariate and multivariate Cox analysis was performed to identify independent prognostic factors. LUAD tissues showed a noticeably higher mRNAsi in than nontumor tissues, and it showed significant differences in T, N, M, AJCC stages, and smoking history. The most frequently mutated gene was TP53, with a higher mRNAsi relating to more frequent mutation of TP53. The mRNAsi was significantly negatively correlated with immune score, stromal score, and ESTIMATE score in LUAD. The blue module was associated with mRNAsi. The 5-gene signature was confirmed as an independent indicator of LUAD prognosis that could promote personalized treatment of LUAD and accurately predict overall survival (OS) of LUAD patients.

Identification of Novel Compound Heterozygous Variants of the PNPLA6 Gene in Boucher-Neuhäuser Syndrome

Background: Boucher–Neuhäuser syndrome (BNS, MIM 215470) is a rare autosomal recessive syndrome caused by mutations in the PNPLA6 gene. Few BNS cases have been reported for functional validation at the RNA level. Herein, we report on the family of a 17-year-old girl with clinical characteristics of BNS, genetic validation, and a systematic review of PNPLA6 variants related to BNS.

Methods: Clinical data and blood samples were collected from the patient and their parents, and whole-exome sequencing was performed and confirmed by Sanger sequencing. RNA-sequencing (RNA-Seq) and quantitative RT-PCR (qRT-PCR) were performed, and the three-dimensional protein structures of the variants were predicted.

Results: We report a 17-year-old female with progressive night blindness since the age of four, primary amenorrhea, and non-development of secondary sexual characteristics. Her impaired vision was diagnosed as retinal pigmentary degeneration of the retina. She had congenital hypogonadotropic hypogonadism (CHH) but no cerebellar ataxia at present. Two novel compound heterozygous variants (c.2241del/p.Met748TrpfsTer65 and c.2986A>G/p.Thr996Ala) of the PNPLA6 gene (NM_006702.4) were identified by whole-exome sequencing. The former variant was carried from her healthy father and has not been reported previously. The latter was inherited from her healthy mother and was noted in a report without functional studies. The RT-PCR results showed that the mRNA expression of PNPLA6 was lower in this patient and her father than in the control group. She was diagnosed with BNS. Both variants (c.2241del and c.2986A>G) were likely pathogenic according to the ACMG criteria. The novel variants in the PNPLA6 gene related to Boucher–Neuhäuser syndrome were summarized in this article.

Conclusion: The possibility of Boucher–Neuhäuser syndrome should be considered when patients present with night blindness, impaired vision, and hypogonadotropic hypogonadism. Gene sequencing is currently the primary diagnostic method. Herein, novel compound heterozygous variants of PNPLA6 were identified in a BNS patient, and its function was verified at the RNA level. The PNPLA6 c.2241del variant is novel and potentially pathogenic, expanding the mutation spectrum in PNPLA6.

Fusion Gene Detection Using Whole-Exome Sequencing Data in Cancer Patients

Several fusion genes are directly involved in the initiation and progression of cancers. Numerous bioinformatics tools have been developed to detect fusion events, but they are mainly based on RNA-seq data. The whole-exome sequencing (WES) represents a powerful technology that is widely used for disease-related DNA variant detection. In this study, we build a novel analysis pipeline called Fuseq-WES to detect fusion genes at DNA level based on the WES data. The same method applies also for targeted panel sequencing data. We assess the method to real datasets of acute myeloid leukemia (AML) and prostate cancer patients. The result shows that two of the main AML fusion genes discovered in RNA-seq data, PML-RARA and CBFB-MYH11, are detected in the WES data in 36 and 63% of the available samples, respectively. For the targeted deep-sequencing of prostate cancer patients, detection of the TMPRSS2-ERG fusion, which is the most frequent chimeric alteration in prostate cancer, is 91% concordant with a manually curated procedure based on four other methods. In summary, the overall results indicate that it is challenging to detect fusion genes in WES data with a standard coverage of ∼ 15–30x, where fusion candidates discovered in the RNA-seq data are often not detected in the WES data and vice versa. A subsampling study of the prostate data suggests that a coverage of at least 75x is necessary to achieve high accuracy.

Theranostic Interpolation of Genomic Instability in Breast Cancer

Breast cancer is a diverse disease caused by mutations in multiple genes accompanying epigenetic aberrations of hazardous genes and protein pathways, which distress tumor-suppressor genes and the expression of oncogenes. Alteration in any of the several physiological mechanisms such as cell cycle checkpoints, DNA repair machinery, mitotic checkpoints, and telomere maintenance results in genomic instability. Theranostic has the potential to foretell and estimate therapy response, contributing a valuable opportunity to modify the ongoing treatments and has developed new treatment strategies in a personalized manner. “Omics” technologies play a key role while studying genomic instability in breast cancer, and broadly include various aspects of proteomics, genomics, metabolomics, and tumor grading. Certain computational techniques have been designed to facilitate the early diagnosis of cancer and predict disease-specific therapies, which can produce many effective results. Several diverse tools are used to investigate genomic instability and underlying mechanisms. The current review aimed to explore the genomic landscape, tumor heterogeneity, and possible mechanisms of genomic instability involved in initiating breast cancer. We also discuss the implications of computational biology regarding mutational and pathway analyses, identification of prognostic markers, and the development of strategies for precision medicine. We also review different technologies required for the investigation of genomic instability in breast cancer cells, including recent therapeutic and preventive advances in breast cancer.

Comparison of seven SNP calling pipelines for the next-generation sequencing data of chickens

Single nucleotide polymorphisms (SNPs) are widely used in genome-wide association studies and population genetics analyses. Next-generation sequencing (NGS) has become convenient, and many SNP-calling pipelines have been developed for human NGS data. We took advantage of a gap knowledge in selecting the appropriated SNP calling pipeline to handle with high-throughput NGS data. To fill this gap, we studied and compared seven SNP calling pipelines, which include 16GT, genome analysis toolkit (GATK), Bcftools-single (Bcftools single sample mode), Bcftools-multiple (Bcftools multiple sample mode), VarScan2-single (VarScan2 single sample mode), VarScan2-multiple (VarScan2 multiple sample mode) and Freebayes pipelines, using 96 NGS data with the different depth gradients of approximately 5X, 10X, 20X, 30X, 40X, and 50X coverage from 16 Rhode Island Red chickens. The sixteen chickens were also genotyped with a 50K SNP array, and the sensitivity and specificity of each pipeline were assessed by comparison to the results of SNP arrays. For each pipeline, except Freebayes, the number of detected SNPs increased as the input read depth increased. In comparison with other pipelines, 16GT, followed by Bcftools-multiple, obtained the most SNPs when the input coverage exceeded 10X, and Bcftools-multiple obtained the most when the input was 5X and 10X. The sensitivity and specificity of each pipeline increased with increasing input. Bcftools-multiple had the highest sensitivity numerically when the input ranged from 5X to 30X, and 16GT showed the highest sensitivity when the input was 40X and 50X. Bcftools-multiple also had the highest specificity, followed by GATK, at almost all input levels. For most calling pipelines, there were no obvious changes in SNP numbers, sensitivities or specificities beyond 20X. In conclusion, (1) if only SNPs were detected, the sequencing depth did not need to exceed 20X; (2) the Bcftools-multiple may be the best choice for detecting SNPs from chicken NGS data, but for a single sample or sequencing depth greater than 20X, 16GT was recommended. Our findings provide a reference for researchers to select suitable pipelines to obtain SNPs from the NGS data of chickens or nonhuman animals.

Novel variants underlying autosomal recessive neurodevelopmental disorders with intellectual disability in Iranian consanguineous families

Background

Intellectual disability (ID) is a heterogeneous group of neurodevelopmental disorders that is characterized by significant impairment in intellectual and adaptive functioning with onset during the developmental period. Whole‐exome sequencing (WES)‐based studies in the consanguineous families with individuals affected with ID have shown a high burden of relevant variants. So far, over 700 genes have been reported in syndromic and non‐syndromic ID. However, genetic causes in more than 50% of ID patients still remain unclear.

Methods

Whole‐exome sequencing was applied for investigation of various variants of ID, then Sanger sequencing and in silico analysis in ten patients from five Iranian consanguineous families diagnosed with autosomal recessive neurodevelopmental disorders, intellectual disability, performed for confirming the causative mutation within the probands. The most patients presented moderate‐to‐severe intellectual disability, developmental delay, seizure, speech problem, high level of lactate, and onset before 10 years.

Results

Filtering the data identified by WES, two novel homozygous missense variants in FBXO31 and TIMM50 genes and one previously reported mutation in the CEP290 gene in the probands were found. Sanger sequencing confirmed the homozygote variant's presence of TIMM50 and FBXO31 genes in six patients and two affected siblings in their respective families. Our computational results predicted that the variants are located in the conserved regions across different species and have the impacts on the protein stability.

Conclusion

Hence, we provide evidence for the pathogenicity of two novel variants in the patients which will expand our knowledge about potential mutation involved in the heterogeneous disease.

Chromosome-Level Genome Assembly and HazelOmics Database Construction Provides Insights Into Unsaturated Fatty Acid Synthesis and Cold Resistance in Hazelnut (Corylus heterophylla)

Corylus heterophylla (2n = 22) is the most widely distributed, unique, and economically important nut species in China. Chromosome-level genomes of C. avellana, C. heterophylla, and C. mandshurica have been published in 2021, but a satisfactory hazelnut genome database is absent. Northeast China is the main distribution and cultivation area of C. heterophylla, and the mechanism underlying the adaptation of C. heterophylla to extremely low temperature in this area remains unclear. Using single-molecule real-time sequencing and the chromosomal conformational capture (Hi-C) assisted genome assembly strategy, we obtained a high-quality chromosome-scale genome sequence of C. heterophylla, with a total length of 343 Mb and scaffold N50 of 32.88 Mb. A total of 94.72% of the test genes from the assembled genome could be aligned to the Embryophyta_odb9 database. In total, 22,319 protein-coding genes were predicted, and 21,056 (94.34%) were annotated in the assembled genome. A HazelOmics online database (HOD) containing the assembled genome, gene-coding sequences, protein sequences, and various types of annotation information was constructed. This database has a user-friendly and straightforward interface. In total, 439 contracted genes and 3,810 expanded genes were identified through genome evolution analysis, and 17 expanded genes were significantly enriched in the unsaturated fatty acid biosynthesis pathway (ko01040). Transcriptome analysis results showed that FAD (Cor0058010.1), SAD (Cor0141290.1), and KAT (Cor0122500.1) with high expression abundance were upregulated at the ovule maturity stage. We deduced that the expansion of these genes may promote high unsaturated fatty acid content in the kernels and improve the adaptability of C. heterophylla to the cold climate of Northeast China. The reference genome and database will be beneficial for future molecular breeding and gene function studies in this nut species, as well as for evolutionary research on species of the order Fagales.

In silico epitope prediction analyses highlight the potential for distracting antigen immunodominance with allogeneic cancer vaccines

Allogeneic cancer vaccines are designed to induce antitumor immune responses with the goal of impacting tumor growth. Typical allogeneic cancer vaccines are produced by expansion of established cancer cell lines, transfection with vectors encoding immunostimulatory cytokines, and lethal irradiation. More than 100 clinical trials have investigated the clinical benefit of allogeneic cancer vaccines in various cancer types. Results show limited therapeutic benefit in clinical trials and currently there are no FDA approved allogeneic cancer vaccines. We used recently developed bioinformatics tools including the pVAC-seq suite of software tools to analyze DNA/RNA sequencing data from the TCGA to examine the repertoire of antigens presented by a typical allogeneic cancer vaccine, and to simulate allogeneic cancer vaccine clinical trials. Specifically, for each simulated clinical trial we modeled the repertoire of antigens presented by allogeneic cancer vaccines consisting of three hypothetical cancer cell lines to 30 patients with the same cancer type. Simulations were repeated ten times for each cancer type. Each tumor sample in the vaccine and the vaccine recipient was subjected to HLA typing, differential expression analyses for tumor associated antigens (TAAs), germline variant calling, and neoantigen prediction. These analyses provided a robust, quantitative comparison between potentially beneficial TAAs and neoantigens versus distracting antigens present in the allogeneic cancer vaccines. We observe that distracting antigens greatly outnumber shared TAAs and neoantigens, providing one potential explanation for the lack of observed responses to allogeneic cancer vaccines. This analysis provides additional rationale for the redirection of efforts towards a personalized cancer vaccine approach.

Genetic alteration and clonal evolution of primary glioblastoma into secondary gliosarcoma

Aims

Secondary gliosarcoma (SGS) rarely arises post treatment of primary glioblastoma multiforme (GBM), and contains gliomatous and sarcomatous components. The origin and clonal evolution of SGS sarcomatous components remain uncharacterized. Therapeutic radiation is mutagenic and can induce sarcomas in patients with other tumor phenotypes, but possible causal relationships between radiotherapy and induction of SGS sarcomatous components remain unexplored. Herein, we investigated the clonal origin of SGS in a patient with primary GBM progressing into SGS post‐radiochemotherapy.

Methods

Somatic mutation profile in GBM and SGS was examined using whole‐genome sequencing and deep‐whole‐exome sequencing. Mutation signatures were characterized to investigate relationships between radiochemotherapy and SGS pathogenesis.

Results

A mutation cluster containing two founding mutations in tumor‐suppressor genes NF1 (variant allele frequency [VAF]: 50.0% in GBM and 51.1% in SGS) and TP53 (VAF: 26.7% in GBM and 50.8% in SGS) was shared in GBM and SGS. SGS exhibited an overpresented C>A (G>T) transversion (oxidative DNA damage signature) but no signature 11 mutations (alkylating‐agents – exposure signature). Since radiation induces DNA lesions by generating reactive oxygen species, the mutations observed in this case of SGS were likely the result of radiotherapy rather than chemotherapy.

Conclusions

Secondary gliosarcoma components likely have a monoclonal origin, and the clone possessing mutations in NF1 and TP53 was likely the founding clone in this case of SGS.

Tumor-infiltrating CD8+ T cells recognize a heterogeneously expressed functional neoantigen in clear cell renal cell carcinoma

Immune checkpoint inhibitors (ICIs) are used in cancer immunotherapy to block programmed death-1 and cytotoxic T-lymphocyte antigen 4, but the response rate for ICIs is still low and tumor cell heterogeneity is considered to be responsible for resistance to immunotherapy. Tumor-infiltrating lymphocytes (TILs) have an essential role in the anti-tumor effect of cancer immunotherapy; however, the specificity of TILs in renal cell carcinoma (RCC) is elusive. In this study, we analyzed a 58-year-old case with clear cell RCC (ccRCC) with the tumor showing macroscopic and microscopic heterogeneity. The tumor was composed of low-grade and high-grade ccRCC. A tumor cell line (1226 RCC cells) and TILs were isolated from the high-grade ccRCC lesion, and a TIL clone recognized a novel neoantigen peptide (YVVPGSPCL) encoded by a missense mutation of the tensin 1 (TNS1) gene in a human leukocyte antigen-C*03:03-restricted fashion. The TNS1 gene mutation was not detected in the low-grade ccRCC lesion and the TIL clone did not recognized low-grade ccRCC cells. The missense mutation of TNS1 encoding the S1309Y mutation was found to be related to cell migration by gene over-expression. These findings suggest that macroscopically and microscopically heterogenous tumors might show heterogenous gene mutations and reactivity to TILs.

Genomic Sequence Analysis of Methicillin- and Carbapenem-Resistant Bacteria Isolated from Raw Sewage

Antibiotic resistance is one of the largest threats facing global health. Wastewater treatment plants are well-known hot spots for interaction between diverse bacteria, genetic exchange, and antibiotic resistance. Nonpathogenic bacteria theoretically act as reservoirs of antibiotic resistance subsequently transferring antibiotic resistance genes to pathogens, indicating that evolutionary processes occur outside clinical settings and may drive patterns of drug-resistant infections. We isolated and sequenced 100 bacterial strains from five wastewater treatment plants to analyze regional dynamics of antibiotic resistance in the California Central Valley. The results demonstrate the presence of a wide diversity of pathogenic and nonpathogenic bacteria, with an arithmetic mean of 5.1 resistance genes per isolate. Forty-three percent of resistance genes were located on plasmids, suggesting that large levels of gene transfer between bacteria that otherwise may not co-occur are facilitated by wastewater treatment. One of the strains detected was a Bacillus carrying pX01 and pX02 anthrax-like plasmids and multiple drug resistance genes. A correlation between resistance genes and taxonomy indicates that taxon-specific evolutionary studies may be useful in determining and predicting patterns of antibiotic resistance. Conversely, a lack of geographic correlation may indicate that landscape genetic studies to understand the spread of antibiotic resistance genes should be carried out at broader scales. This large data set provides insights into how pathogenic and nonpathogenic bacteria interact in wastewater environments and the resistance genes which may be horizontally transferred between them. This can help in determining the mechanisms leading to the increasing prevalence of drug-resistant infections observed in clinical settings. IMPORTANCE The reasons for the increasing prevalence of antibiotic-resistant infections are complex and associated with myriad clinical and environmental processes. Wastewater treatment plants operate as nexuses of bacterial interaction and are known hot spots for genetic exchange between bacteria, including antibiotic resistance genes. We isolated and sequenced 100 drug-resistant bacteria from five wastewater treatment plants in California's Central Valley, characterizing widespread gene sharing between pathogens and nonpathogens. We identified a novel, multiresistant Bacillus carrying anthrax-like plasmids. This empirical study supports the likelihood of evolutionary and population processes in the broader environment affecting the prevalence of clinical drug-resistant infections and identifies several taxa that may operate as reservoirs and vectors of antibiotic resistance genes.

Prospects and challenges of cancer systems medicine: from genes to disease networks

It is becoming evident that holistic perspectives toward cancer are crucial in deciphering the overwhelming complexity of tumors. Single-layer analysis of genome-wide data has greatly contributed to our understanding of cellular systems and their perturbations. However, fundamental gaps in our knowledge persist and hamper the design of effective interventions. It is becoming more apparent than ever, that cancer should not only be viewed as a disease of the genome but as a disease of the cellular system. Integrative multilayer approaches are emerging as vigorous assets in our endeavors to achieve systemic views on cancer biology. Herein, we provide a comprehensive review of the approaches, methods and technologies that can serve to achieve systemic perspectives of cancer. We start with genome-wide single-layer approaches of omics analyses of cellular systems and move on to multilayer integrative approaches in which in-depth descriptions of proteogenomics and network-based data analysis are provided. Proteogenomics is a remarkable example of how the integration of multiple levels of information can reduce our blind spots and increase the accuracy and reliability of our interpretations and network-based data analysis is a major approach for data interpretation and a robust scaffold for data integration and modeling. Overall, this review aims to increase cross-field awareness of the approaches and challenges regarding the omics-based study of cancer and to facilitate the necessary shift toward holistic approaches.

Technological advances in cancer immunity: from immunogenomics to single-cell analysis and artificial intelligence

Immunotherapies play critical roles in cancer treatment. However, given that only a few patients respond to immune checkpoint blockades and other immunotherapeutic strategies, more novel technologies are needed to decipher the complicated interplay between tumor cells and the components of the tumor immune microenvironment (TIME). Tumor immunomics refers to the integrated study of the TIME using immunogenomics, immunoproteomics, immune-bioinformatics, and other multi-omics data reflecting the immune states of tumors, which has relied on the rapid development of next-generation sequencing. High-throughput genomic and transcriptomic data may be utilized for calculating the abundance of immune cells and predicting tumor antigens, referring to immunogenomics. However, as bulk sequencing represents the average characteristics of a heterogeneous cell population, it fails to distinguish distinct cell subtypes. Single-cell-based technologies enable better dissection of the TIME through precise immune cell subpopulation and spatial architecture investigations. In addition, radiomics and digital pathology-based deep learning models largely contribute to research on cancer immunity. These artificial intelligence technologies have performed well in predicting response to immunotherapy, with profound significance in cancer therapy. In this review, we briefly summarize conventional and state-of-the-art technologies in the field of immunogenomics, single-cell and artificial intelligence, and present prospects for future research.

An integrated DNA and RNA variant detector identifies a highly conserved three base exon in the MAP4K5 kinase locus

RNA variants that emerge from editing and alternative splicing form important regulatory stages in protein signalling. In this report, we apply an integrated DNA and RNA variant detection workbench to define the range of RNA variants that deviate from the reference genome in a human melanoma cell model. The RNA variants can be grouped into (i) classic ADAR-like or APOBEC-like RNA editing events and (ii) multiple-nucleotide variants (MNVs) including three and six base pair in-frame non-canonical unmapped exons. We focus on validating representative genes of these classes. First, clustered non-synonymous RNA edits (A-I) in the CDK13 gene were validated by Sanger sequencing to confirm the integrity of the RNA variant detection workbench. Second, a highly conserved RNA variant in the MAP4K5 gene was detected that results most likely from the splicing of a non-canonical three-base exon. The two RNA variants produced from the MAP4K5 locus deviate from the genomic reference sequence and produce V569E or V569del isoform variants. Low doses of splicing inhibitors demonstrated that the MAP4K5-V569E variant emerges from an SF3B1-dependent splicing event. Mass spectrometry of the recombinant SBP-tagged MAP4K5^V569E and MAP4K5^V569del proteins pull-downs in transfected cell systems was used to identify the protein-protein interactions of these two MAP4K5 isoforms and propose possible functions. Together these data highlight the utility of this integrated DNA and RNA variant detection platform to detect RNA variants in cancer cells and support future analysis of RNA variant detection in cancer tissue.

Immunogenomic profiling and pathological response results from a clinical trial of docetaxel and carboplatin in triple-negative breast cancer

PURPOSE

Patients with triple-negative breast cancer (TNBC) who do not achieve pathological complete response (pCR) following neoadjuvant chemotherapy have a high risk of recurrence and death. Molecular characterization may identify patients unlikely to achieve pCR. This neoadjuvant trial was conducted to determine the pCR rate with docetaxel and carboplatin and to identify molecular alterations and/or immune gene signatures predicting pCR.

EXPERIMENTAL DESIGN

Patients with clinical stages II/III TNBC received 6 cycles of docetaxel and carboplatin. The primary objective was to determine if neoadjuvant docetaxel and carboplatin would increase the pCR rate in TNBC compared to historical expectations. We performed whole-exome sequencing (WES) and immune profiling on pre-treatment tumor samples to identify alterations that may predict pCR. Thirteen matching on-treatment samples were also analyzed to assess changes in molecular profiles.

RESULTS

Fifty-eight of 127 (45.7%) patients achieved pCR. There was a non-significant trend toward higher mutation burden for patients with residual cancer burden (RCB) 0/I versus RCB II/III (median 80 versus 68 variants, p 0.88). TP53 was the most frequently mutated gene, observed in 85.7% of tumors. EGFR, RB1, RAD51AP2, SDK2, L1CAM, KPRP, PCDHA1, CACNA1S, CFAP58, COL22A1, and COL4A5 mutations were observed almost exclusively in pre-treatment samples from patients who achieved pCR. Seven mutations in PCDHA1 were observed in pre-treatment samples from patients who did not achieve pCR. Several immune gene signatures including IDO1, PD-L1, interferon gamma signaling, CTLA4, cytotoxicity, tumor inflammation signature, inflammatory chemokines, cytotoxic cells, lymphoid, PD-L2, exhausted CD8, Tregs, and immunoproteasome were upregulated in pre-treatment samples from patients who achieved pCR.

CONCLUSION

Neoadjuvant docetaxel and carboplatin resulted in a pCR of 45.7%. WES and immune profiling differentiated patients with and without pCR.

TRIAL REGISTRATION

Clinical trial information: NCT02124902, Registered 24 April 2014 & NCT02547987, Registered 10 September 2015.

Longitudinal single-cell profiling reveals molecular heterogeneity and tumor-immune evolution in refractory mantle cell lymphoma

The mechanisms driving therapeutic resistance and poor outcomes of mantle cell lymphoma (MCL) are incompletely understood. We characterize the cellular and molecular heterogeneity within and across patients and delineate the dynamic evolution of tumor and immune cell compartments at single cell resolution in longitudinal specimens from ibrutinib-sensitive patients and non-responders. Temporal activation of multiple cancer hallmark pathways and acquisition of 17q are observed in a refractory MCL. Multi-platform validation is performed at genomic and cellular levels in PDX models and larger patient cohorts. We demonstrate that due to 17q gain, BIRC5/survivin expression is upregulated in resistant MCL tumor cells and targeting BIRC5 results in marked tumor inhibition in preclinical models. In addition, we discover notable differences in the tumor microenvironment including progressive dampening of CD8+ T cells and aberrant cell-to-cell communication networks in refractory MCLs. This study reveals diverse and dynamic tumor and immune programs underlying therapy resistance in MCL.

Sex chromosome transformation and the origin of a male-specific X chromosome in the creeping vole

The mammalian sex chromosome system (XX female/XY male) is ancient and highly conserved. The sex chromosome karyotype of the creeping vole (Microtus oregoni) represents a long-standing anomaly, with an X chromosome that is unpaired in females (X0) and exclusively maternally transmitted. We produced a highly contiguous male genome assembly, together with short-read genomes and transcriptomes for both sexes. We show that M. oregoni has lost an independently segregating Y chromosome and that the male-specific sex chromosome is a second X chromosome that is largely homologous to the maternally transmitted X. Both maternally inherited and male-specific sex chromosomes carry fragments of the ancestral Y chromosome. Consequences of this recently transformed sex chromosome system include Y-like degeneration and gene amplification on the male-specific X, expression of ancestral Y-linked genes in females, and X inactivation of the male-specific chromosome in male somatic cells. The genome of M. oregoni elucidates the processes that shape the gene content and dosage of mammalian sex chromosomes and exemplifies a rare case of plasticity in an ancient sex chromosome system.

Detection of new drivers of frequent B-cell lymphoid neoplasms using an integrated analysis of whole genomes

B-cell lymphoproliferative disorders exhibit a diverse spectrum of diagnostic entities with heterogeneous behaviour. Multiple efforts have focused on the determination of the genomic drivers of B-cell lymphoma subtypes. In the meantime, the aggregation of diverse tumors in pan-cancer genomic studies has become a useful tool to detect new driver genes, while enabling the comparison of mutational patterns across tumors. Here we present an integrated analysis of 354 B-cell lymphoid disorders. 112 recurrently mutated genes were discovered, of which KMT2D, CREBBP, IGLL5 and BCL2 were the most frequent, and 31 genes were putative new drivers. Mutations in CREBBP, TNFRSF14 and KMT2D predominated in follicular lymphoma, whereas those in BTG2, HTA-A and PIM1 were more frequent in diffuse large B-cell lymphoma. Additionally, we discovered 31 significantly mutated protein networks, reinforcing the role of genes such as CREBBP, EEF1A1, STAT6, GNA13 and TP53, but also pointing towards a myriad of infrequent players in lymphomagenesis. Finally, we report aberrant expression of oncogenes and tumor suppressors associated with novel noncoding mutations (DTX1 and S1PR2), and new recurrent copy number aberrations affecting immune check-point regulators (CD83, PVR) and B-cell specific genes (TNFRSF13C). Our analysis expands the number of mutational drivers of B-cell lymphoid neoplasms, and identifies several differential somatic events between disease subtypes.

FGFR2 Extracellular Domain In-Frame Deletions are Therapeutically Targetable Genomic Alterations that Function as Oncogenic Drivers in Cholangiocarcinoma

We conducted next generation DNA sequencing on 335 biliary tract cancers and characterized the genomic landscape by anatomic site within the biliary tree. In addition to frequent FGFR2 fusions among patients with intrahepatic cholangiocarcinoma (IHCC), we identified FGFR2 extracellular domain in-frame deletions (EIDs) in 5 of 178 (2.8%) patients with IHCC, including two patients with FGFR2 p.H167_N173del. Expression of this FGFR2 EID in NIH3T3 cells resulted in constitutive FGFR2 activation, oncogenic transformation, and sensitivity to FGFR inhibitors. Three patients with FGFR2 EIDs were treated with Debio 1347, an oral FGFR-1/2/3 inhibitor, and all showed partial responses. One patient developed an acquired L618F FGFR2 kinase domain mutation at disease progression and experienced a further partial response for 17 months to an irreversible FGFR2 inhibitor, futibatinib. Together, these findings reveal FGFR2 EIDs as an alternative mechanism of FGFR2 activation in IHCC that predict sensitivity to FGFR inhibitors in the clinic.

Genome-wide transcriptome and translatome analyses reveal the role of protein extension and domestication in liver cancer oncogenesis

One gene could be transcribed to different RNA isoforms, and then produce various forms of protein sequences. This mechanism largely diversifies the cellular pool and allows natural selection to select from a wider range of substrates. In the cancer field, the isoform switches between tumor and normal tissues, such as the alternative splicing, stop codon read-through, or protein domestication, are significantly ignored by the traditional differential expression analyses. The intention of this work is to fill this gap. We collected public transcriptome and translatome data from ten patients with liver cancer, and performed genome-wide comparison on the stop codon read-through and protein domestication events. Both events diversify the proteome during long-term evolution. Surprisingly, we found that the tumor tissues globally have higher occurrence of stop codon read-through events as well as protein domestication events (translation signals of non-coding repetitive elements). These read-through and domestication events show limited overlapping across the ten patients, indicating the randomness of their occurrence and their deleterious nature. These tumor-specific events might have been purged by natural selection if they are not collected timely. Our work manifests the role of protein extension and domestication in liver cancer oncogenesis, adding new aspects to the cancer field.

Genomic profile of metastatic breast cancer patient-derived xenografts established using percutaneous biopsy

Background

Metastatic breast cancer (mBC) is a complex and life-threatening disease and although it is difficult to cure, patients can benefit from sequential anticancer treatment, including endocrine therapy, targeted therapy and cytotoxic chemotherapy. The patient-derived xenograft (PDX) model is suggested as a practical tool to predict the clinical outcome of this disease as well as to screen novel drugs. This study aimed to establish PDX models in Korean patients and analyze their genomic profiles and utility for translational research.

Methods

Percutaneous core needle biopsy or punch biopsy samples were used for xenotransplantation. Whole exome sequencing and transcriptome analysis were performed to assess the genomic and RNA expression profiles, respectively. Copy number variation and mutational burden were analyzed and compared with other metastatic breast cancer genomic results. Mutational signatures were also analyzed. The antitumor effect of an ATR inhibitor was tested in the relevant PDX model.

Results

Of the 151 cases studied, 40 (26%) PDX models were established. Notably, the take rate of all subtypes, including the hormone receptor-positive (HR +) subtype, exceeded 20%. The PDX model had genomic fidelity and copy number variation that represented the pattern of its donor sample. TP53, PIK3CA, ESR1, and GATA3 mutations were frequently found in our samples, with TP53 being the most frequently mutated, and the somatic mutations in these genes strengthened their frequency in the PDX model. The ESR1 mutation, CCND1 amplification, and the APOBEC signature were significant features in our HR + HER2- PDX model. Fulvestrant in combination with palbociclib showed a partial response to the relevant patient’s tumor harboring the ESR1 mutation, and CCND1 amplification was found in the PDX model. AZD6738, an ATR inhibitor, delayed tumor growth in a relevant PDX model.

Conclusions

Our PDX model was established using core needle biopsy samples from primary and metastatic tissues. Genomic profiles of the samples reflected their original tissue characteristics and could be used for the interpretation of clinical outcomes.

Risk prediction and marker selection in nonsynonymous single nucleotide polymorphisms using whole genome sequencing data

Despite the various existing studies about nonsynonymous single nucleotide polymorphisms (nsSNPs), genome-wide studies based on nsSNPs are rare. NsSNPs alter amino acid sequences, affect protein structure and function, and have deleterious effects. By predicting the deleterious effect of nsSNPs, we determined the total risk score per individual. Additionally, the machine learning technique was utilized to find an optimal nsSNP subset that best explains the complete nsSNP effect. A total of 16,100 nsSNPs were selected as the best representatives among 89,519 regressed nsSNPs. In the gene ontology analysis encompassing the 16,100 nsSNPs, DNA metabolic process, chemokine- and immune-related, and reproduction were the most enriched terms. We expect that our risk score prediction and nsSNP marker selection will contribute to future development of extant genome-wide association studies and breeding science more broadly.

Systematically Characterizing A-to-I RNA Editing Neoantigens in Cancer

A-to-I RNA editing can contribute to the transcriptomic and proteomic diversity of many diseases including cancer. It has been reported that peptides generated from RNA editing could be naturally presented by human leukocyte antigen (HLA) molecules and elicit CD8+ T cell activation. However, a systematical characterization of A-to-I RNA editing neoantigens in cancer is still lacking. Here, an integrated RNA-editing based neoantigen identification pipeline PREP (Prioritizing of RNA Editing-based Peptides) was presented. A comprehensive RNA editing neoantigen profile analysis on 12 cancer types from The Cancer Genome Atlas (TCGA) cohorts was performed. PREP was also applied to 14 ovarian tumor samples and two clinical melanoma cohorts treated with immunotherapy. We finally proposed an RNA editing neoantigen immunogenicity score scheme, i.e. REscore, which takes RNA editing level and infiltrating immune cell population into consideration. We reported variant peptide from protein IFI30 in breast cancer which was confirmed expressed and presented in two samples with mass spectrometry data support. We showed that RNA editing neoantigen could be identified from RNA-seq data and could be validated with mass spectrometry data in ovarian tumor samples. Furthermore, we characterized the RNA editing neoantigen profile of clinical melanoma cohorts treated with immunotherapy. Finally, REscore showed significant associations with improved overall survival in melanoma cohorts treated with immunotherapy. These findings provided novel insights of cancer biomarker and enhance our understanding of neoantigen derived from A-to-I RNA editing as well as more types of candidates for personalized cancer vaccines design in the context of cancer immunotherapy.

DNA sequencing of anatomy lab cadavers to provide hands-on precision medicine introduction to medical students

BACKGROUND

Medical treatment informed by Precision Medicine is becoming a standard practice for many diseases, and patients are curious about the consequences of genomic variants in their genome. However, most medical students' understanding of Precision Medicine derives from classroom lectures. This format does little to foster an understanding for the potential and limitations of Precision Medicine. To close this gap, we implemented a hands-on Precision Medicine training program utilizing exome sequencing to prepare a clinical genetic report of cadavers studied in the anatomy lab. The program reinforces Precision Medicine related learning objectives for the Genetics curriculum.

METHODS

Pre-embalmed blood samples and embalmed tissue were obtained from cadavers (donors) used in the anatomy lab. DNA was isolated and sequenced and illustrative genetic reports provided to the students. The reports were used to facilitate discussion with students on the implications of pathogenic genomic variants and the potential correlation of these variants in each "donor" with any anatomical anomalies identified during cadaver dissection.

RESULTS

In 75% of cases, analysis of whole exome sequencing data identified a variant associated with increased risk for a disease/abnormal condition noted in the donor's cause of death or in the students' anatomical findings. This provided students with real-world examples of the potential relationship between genomic variants and disease risk. Our students also noted that diseases associated with 92% of the pathogenic variants identified were not related to the anatomical findings, demonstrating the limitations of Precision Medicine.

CONCLUSION

With this study, we have established protocols and classroom procedures incorporating hands-on Precision Medicine training in the medical student curriculum and a template for other medical educators interested in enhancing their Precision Medicine training program. The program engaged students in discovering variants that were associated with the pathophysiology of the cadaver they were studying, which led to more exposure and understanding of the potential risks and benefits of genomic medicine.

Exploration into biomarker potential of region-specific brain gene co-expression networks

The human brain is a complex organ that consists of several regions each with a unique gene expression pattern. Our intent in this study was to construct a gene co-expression network (GCN) for the normal brain using RNA expression profiles from the Genotype-Tissue Expression (GTEx) project. The brain GCN contains gene correlation relationships that are broadly present in the brain or specific to thirteen brain regions, which we later combined into six overarching brain mini-GCNs based on the brain's structure. Using the expression profiles of brain region-specific GCN edges, we determined how well the brain region samples could be discriminated from each other, visually with t-SNE plots or quantitatively with the Gene Oracle deep learning classifier. Next, we tested these gene sets on their relevance to human tumors of brain and non-brain origin. Interestingly, we found that genes in the six brain mini-GCNs showed markedly higher mutation rates in tumors relative to matched sets of random genes. Further, we found that cortex genes subdivided Head and Neck Squamous Cell Carcinoma (HNSC) tumors and Pheochromocytoma and Paraganglioma (PCPG) tumors into distinct groups. The brain GCN and mini-GCNs are useful resources for the classification of brain regions and identification of biomarker genes for brain related phenotypes.