VarScan: variant detection in massively parallel sequencing of individual and pooled samples

Novel insights into chloroplast genome evolution in the green macroalgal genus Ulva (Ulvophyceae, Chlorophyta)

To understand the evolutionary driving forces of chloroplast (or plastid) genomes (plastomes) in the green macroalgal genus Ulva (Ulvophyceae, Chlorophyta), in this study, we sequenced and constructed seven complete chloroplast genomes from five Ulva species, and conducted comparative genomic analysis of Ulva plastomes in Ulvophyceae. Ulva plastome evolution reflects the strong selection pressure driving the compactness of genome organization and the decrease of overall GC composition. The overall plastome sequences including canonical genes, introns, derived foreign sequences and non-coding regions show a synergetic decrease in GC content at varying degrees. Fast degeneration of plastome sequences including non-core genes (minD and trnR3), derived foreign sequences, and noncoding spacer regions was accompanied by the marked decrease of their GC composition. Plastome introns preferentially resided in conserved housekeeping genes with high GC content and long length, as might be related to high GC content of target site sequences recognized by intron-encoded proteins (IEPs), and to more target sites contained by long GC-rich genes. Many foreign DNA sequences integrated into different intergenic regions contain some homologous specific orfs with high similarity, indicating that they could have been derived from the same origin. The invasion of foreign sequences seems to be an important driving force for plastome rearrangement in these IR-lacking Ulva cpDNAs. Gene partitioning pattern has changed and distribution range of gene clusters has expanded after the loss of IR, indicating that genome rearrangement was more extensive and more frequent in Ulva plastomes, which was markedly different from that in IR-containing ulvophycean plastomes. These new insights greatly enhance our understanding of plastome evolution in ecologically important Ulva seaweeds.

Lactylome analysis suggests lactylation-dependent mechanisms of metabolic adaptation in hepatocellular carcinoma

Enhanced glycolysis and accumulation of lactate is a common feature in various types of cancer. Intracellular lactate drives a recently described type of posttranslational modification, lysine lactylation (Kla), on core histones. However, the impact of lactylation on biological processes of tumour cells remains largely unknown. Here we show a global lactylome profiling on a prospectively collected hepatitis B virus-related hepatocellular carcinoma (HCC) cohort. Integrative lactylome and proteome analysis of the tumours and adjacent livers identifies 9,275 Kla sites, with 9,256 sites on non-histone proteins, indicating that Kla is a prevalent modification beyond histone proteins and transcriptional regulation. Notably, Kla preferentially affects enzymes involved in metabolic pathways, including the tricarboxylic acid cycle, and carbohydrate, amino acid, fatty acid and nucleotide metabolism. We further verify that lactylation at K28 inhibits the function of adenylate kinase 2, facilitating the proliferation and metastasis of HCC cells. Our study therefore reveals that Kla plays an important role in regulating cellular metabolism and may contribute to HCC progression.

Computational studies reveal co-occurrence of two mutations in IL7R gene of high-grade serous carcinoma patients

Major cause of mortality in ovarian cancer can be attributed to a lack of specific and sensitive biomarkers for diagnosis and prognosis of the disease. Uncovering the mutations in genes involved in crucial oncogenic pathways is a key step in discovery and development of novel biomarkers. Whole exome sequencing (WES) is a powerful method for the detection of cancer driver mutations. The present work focuses on identifying functionally damaging mutations in patients with high-grade serous ovarian carcinoma (HGSC) through computational analysis of WES. In this study, WES data of HGSC patients was retrieved from the genomic literature available in sequence read archive, the variants were identified and comprehensive structural and functional analysis was performed. Interestingly, I66T and V138I mutations were found to be co-occurring in the IL7R gene in four out of five HGSC patient samples investigated in this study. The V138I mutation was located in the fibronectin type-3 domain and computationally assessed to be causing disruptive effects on the structure and dynamics of IL7R protein. This mutation was found to be co-occurring with the neutral I66T mutation in the same domain which compensated the disruptive effects of V138I variant. These comprehensive studies point to a hitherto unexplored significant role of the IL7R gene in ovarian carcinoma. It is envisaged that the work will lay the foundation for the development of a novel biomarker with potential application in molecular profiling and in estimation of the disease prognosis.Communicated by Ramaswamy H. Sarma.

XCVATR: detection and characterization of variant impact on the Embeddings of single -cell and bulk RNA-sequencing samples

BACKGROUND

RNA-sequencing has become a standard tool for analyzing gene activity in bulk samples and at the single-cell level. By increasing sample sizes and cell counts, this technique can uncover substantial information about cellular transcriptional states. Beyond quantification of gene expression, RNA-seq can be used for detecting variants, including single nucleotide polymorphisms, small insertions/deletions, and larger variants, such as copy number variants. Notably, joint analysis of variants with cellular transcriptional states may provide insights into the impact of mutations, especially for complex and heterogeneous samples. However, this analysis is often challenging due to a prohibitively high number of variants and cells, which are difficult to summarize and visualize. Further, there is a dearth of methods that assess and summarize the association between detected variants and cellular transcriptional states.

RESULTS

Here, we introduce XCVATR (eXpressed Clusters of Variant Alleles in Transcriptome pRofiles), a method that identifies variants and detects local enrichment of expressed variants within embedding of samples and cells in single-cell and bulk RNA-seq datasets. XCVATR visualizes local "clumps" of small and large-scale variants and searches for patterns of association between each variant and cellular states, as described by the coordinates of cell embedding, which can be computed independently using any type of distance metrics, such as principal component analysis or t-distributed stochastic neighbor embedding. Through simulations and analysis of real datasets, we demonstrate that XCVATR can detect enrichment of expressed variants and provide insight into the transcriptional states of cells and samples. We next sequenced 2 new single cell RNA-seq tumor samples and applied XCVATR. XCVATR revealed subtle differences in CNV impact on tumors.

CONCLUSIONS

XCVATR is publicly available to download from https://github.com/harmancilab/XCVATR .

A Persistent Tuberculosis Outbreak in the UK Is Characterized by Hydrophobic fadB4-Deficient Mycobacterium tuberculosis That Replicates Rapidly in Macrophages

The genetic diversity of Mycobacterium tuberculosis can influence disease severity and transmissibility. To better understand how this diversity influences individuals and communities, we phenotyped M. tuberculosis that was causing a persistent outbreak in the East Midlands, United Kingdom. Compared to nonoutbreak isolates, bacilli had higher lipid contents and more hydrophobic cell surfaces. In macrophage infection models, the bacteria increased more rapidly, provoked the enhanced accumulation of macrophage lipid droplets and enhanced the secretion of IL-1β. Natural deletions in fadB4, nrdB, and plcC distinguished the outbreak isolates from other lineage 3 isolates in the region. fadB4 is annotated with a putative role in cell envelope biosynthesis, so the loss of this gene has the potential to alter the interactions of bacteria with immune cells. Reintroduction of fadB4 to the outbreak strain led to a phenotype that more closely resembled those of nonoutbreak strains. The improved understanding of the microbiological characteristics and the corresponding genetic polymorphisms that associate with outbreaks have the potential to inform tuberculosis control. IMPORTANCE Tuberculosis (TB) killed 1.5 million people in 2020 and affects every country. The extent to which the natural genetic diversity of Mycobacterium tuberculosis influences disease manifestation at both the individual and epidemiological levels remains poorly understood. Insights into how pathogen polymorphisms affect patterns of TB have the potential to translate into clinical and public health practice. Two distinct lineage 3 strains isolated from local TB outbreaks, one of which (CH) was rapidly terminated and the other of which (Lro) persistently transmitted for over a decade, provided us with an opportunity to study these issues. We compared genome sequences, microbiological characteristics, and early immune responses that were evoked upon infection. Our results indicate that the natural lack of fadB4 in the Lro strain contributes to its unique features.

Systems Biology of Aromatic Compound Catabolism in Facultative Anaerobic Aromatoleum aromaticum EbN1T

Members of the genus Aromatoleum thrive in diverse habitats and use a broad range of recalcitrant organic molecules coupled to denitrification or O₂ respiration. To gain a holistic understanding of the model organism A. aromaticum EbN1^T, we studied its catabolic network dynamics in response to 3-(4-hydroxyphenyl)propanoate, phenylalanine, 3-hydroxybenzoate, benzoate, and acetate utilized under nitrate-reducing versus oxic conditions. Integrated multi-omics (transcriptome, proteome, and metabolome) covered most of the catabolic network (199 genes) and allowed for the refining of knowledge of the degradation modules studied. Their substrate-dependent regulation showed differing degrees of specificity, ranging from high with 3-(4-hydroxyphenyl)propanoate to mostly relaxed with benzoate. For benzoate, the transcript and protein formation were essentially constitutive, contrasted by that of anoxia-specific versus oxia-specific metabolite profiles. The matrix factorization of transcriptomic data revealed that the anaerobic modules accounted for most of the variance across the degradation network. The respiration network appeared to be constitutive, both on the transcript and protein levels, except for nitrate reductase (with narGHI expression occurring only under nitrate-reducing conditions). The anoxia/nitrate-dependent transcription of denitrification genes is apparently controlled by three FNR-type regulators as well as by NarXL (all constitutively formed). The resequencing and functional reannotation of the genome fostered a genome-scale metabolic model, which is comprised of 655 enzyme-catalyzed reactions and 731 distinct metabolites. The model predictions for growth rates and biomass yields agreed well with experimental stoichiometric data, except for 3-(4-hydroxyphenyl)propanoate, with which 4-hydroxybenzoate was exported. Taken together, the combination of multi-omics, growth physiology, and a metabolic model advanced our knowledge of an environmentally relevant microorganism that differs significantly from other bacterial model strains. IMPORTANCE Aromatic compounds are abundant constituents not only of natural organic matter but also of bulk industrial chemicals and fuel components of environmental concern. Considering the widespread occurrence of redox gradients in the biosphere, facultative anaerobic degradation specialists can be assumed to play a prominent role in the natural mineralization of organic matter and in bioremediation at contaminated sites. Surprisingly, differential multi-omics profiling of the A. aromaticum EbN1^T studied here revealed relaxed regulatory stringency across its four main physiological modi operandi (i.e., O₂-independent and O₂-dependent degradation reactions versus denitrification and O₂ respiration). Combining multi-omics analyses with a genome-scale metabolic model aligned with measured growth performances establishes A. aromaticum EbN1^T as a systems-biology model organism and provides unprecedented insights into how this bacterium functions on a holistic level. Moreover, this experimental platform invites future studies on eco-systems and synthetic biology of the environmentally relevant betaproteobacterial Aromatoleum/Azoarcus/Thauera cluster.

Integrative Meta-Analysis of Huntington's Disease Transcriptome Landscape

Huntington's disease (HD) is a neurodegenerative disorder with autosomal dominant inheritance caused by glutamine expansion in the Huntingtin gene (HTT). Striatal projection neurons (SPNs) in HD are more vulnerable to cell death. The executive striatal population is directly connected with the Brodmann Area (BA9), which is mainly involved in motor functions. Analyzing the disease samples from BA9 from the SRA database provides insights related to neuron degeneration, which helps to identify a promising therapeutic strategy. Most gene expression studies examine the changes in expression and associated biological functions. In this study, we elucidate the relationship between variants and their effect on gene/downstream transcript expression. We computed gene and transcript abundance and identified variants from RNA-seq data using various pipelines. We predicted the effect of genome-wide association studies (GWAS)/novel variants on regulatory functions. We found that many variants affect the histone acetylation pattern in HD, thereby perturbing the transcription factor networks. Interestingly, some variants affect miRNA binding as well as their downstream gene expression. Tissue-specific network analysis showed that mitochondrial, neuroinflammation, vasculature, and angiogenesis-related genes are disrupted in HD. From this integrative omics analysis, we propose that abnormal neuroinflammation acts as a two-edged sword that indirectly affects the vasculature and associated energy metabolism. Rehabilitation of blood-brain barrier functionality and energy metabolism may secure the neuron from cell death.

Chloroplast genome assemblies and comparative analyses of commercially important Vaccinium berry crops

Vaccinium is a large genus of shrubs that includes a handful of economically important berry crops. Given the numerous hybridizations and polyploidization events, the taxonomy of this genus has remained the subject of long debate. In addition, berries and berry-based products are liable to adulteration, either fraudulent or unintentional due to misidentification of species. The availability of more genomic information could help achieve higher phylogenetic resolution for the genus, provide molecular markers for berry crops identification, and a framework for efficient genetic engineering of chloroplasts. Therefore, in this study we assembled five Vaccinium chloroplast sequences representing the economically relevant berry types: northern highbush blueberry (V. corymbosum), southern highbush blueberry (V. corymbosum hybrids), rabbiteye blueberry (V. virgatum), lowbush blueberry (V. angustifolium), and bilberry (V. myrtillus). Comparative analyses showed that the Vaccinium chloroplast genomes exhibited an overall highly conserved synteny and sequence identity among them. Polymorphic regions included the expansion/contraction of inverted repeats, gene copy number variation, simple sequence repeats, indels, and single nucleotide polymorphisms. Based on their in silico discrimination power, we suggested variants that could be developed into molecular markers for berry crops identification. Phylogenetic analysis revealed multiple origins of highbush blueberry plastomes, likely due to the hybridization events that occurred during northern and southern highbush blueberry domestication.

Anaeropeptidivorans aminofermentans gen. nov., sp. nov., a mesophilic proteolytic salt-tolerant bacterium isolated from a laboratory-scale biogas fermenter, and emended description of Clostridium colinum

An anaerobic bacterial strain, designated strain M3/9^T, was isolated from a laboratory-scale biogas fermenter fed with maize silage supplemented with 5 % wheat straw. Cells were straight, non-motile rods, which stained Gram-negative. Optimal growth occurred between 30 and 40°C, at pH 7.5-8.5, and up to 3.9 % (w/v) NaCl was tolerated. When grown on peptone from casein and soymeal, strain M3/9^T produced mainly acetic acid, ethanol, and isobutyric acid. The major cellular fatty acids of the novel strain were C_16 : 0 and C_16 : 0 DMA. The genome of strain M3/9^T is 3757  330 bp in size with a G+C content of 38.45 mol%. Phylogenetic analysis allocated strain M3/9^T within the family Lachnospiraceae with Clostridium colinum DSM 6011^T and Anaerotignum lactatifermentans DSM 14214^T being the most closely related species sharing 57.86 and 56.99% average amino acid identity and 16S rRNA gene sequence similarities of 91.58 and 91.26 %, respectively. Based on physiological, chemotaxonomic and genetic data, we propose the description of a novel species and genus Anaeropeptidivorans aminofermentans gen. nov., sp. nov., represented by the type strain M3/9^T (=DSM 100058^T=LMG 29527^T). In addition, an emended description of Clostridium colinum is provided.

An ancestral mycobacterial effector promotes dissemination of infection

The human pathogen Mycobacterium tuberculosis typically causes lung disease but can also disseminate to other tissues. We identified a M. tuberculosis (Mtb) outbreak presenting with unusually high rates of extrapulmonary dissemination and bone disease. We found that the causal strain carried an ancestral full-length version of the type VII-secreted effector EsxM rather than the truncated version present in other modern Mtb lineages. The ancestral EsxM variant exacerbated dissemination through enhancement of macrophage motility, increased egress of macrophages from established granulomas, and alterations in macrophage actin dynamics. Reconstitution of the ancestral version of EsxM in an attenuated modern strain of Mtb altered the migratory mode of infected macrophages, enhancing their motility. In a zebrafish model, full-length EsxM promoted bone disease. The presence of a derived nonsense variant in EsxM throughout the major Mtb lineages 2, 3, and 4 is consistent with a role for EsxM in regulating the extent of dissemination.

Inferring bona fide Differentially Expressed Genes and Their Variants Associated with Vitamin K Deficiency Using a Systems Genetics Approach

Systems genetics is key for integrating a large number of variants associated with diseases. Vitamin K (VK) is one of the scarcely studied disease conditions. In this work, we ascertained the differentially expressed genes (DEGs) and variants associated with individual subpopulations of VK disease phenotypes, viz., myocardial infarction, renal failure and prostate cancer. We sought to ask whether or not any DEGs harbor pathogenic variants common in these conditions, attempt to bridge the gap in finding characteristic biomarkers and discuss the role of long noncoding RNAs (lncRNAs) in the biogenesis of VK deficiencies.

InDelGT: An integrated pipeline for extracting indel genotypes for genetic mapping in a hybrid population using next-generation sequencing data

Premise

Although several software packages are available for genotyping insertion/deletion (indel) polymorphisms in genomes using next-generation sequencing data, simultaneously calling indel genotypes across many individuals for use in genetic mapping remains challenging.

Methods and Results

We present an integrated pipeline, InDelGT, for the extraction of indel genotypes from a segregating population such as backcross or F2 lines, or from an F1 cross between outbred species. The InDelGT algorithm is implemented in three steps: generating an indel catalog, calling indel genotypes, and analyzing indel segregation. We demonstrated the use of the pipeline with an example data set from an F1 hybrid population of Populus and successfully constructed the two parental genetic linkage maps.

Conclusions

InDelGT is a practical tool that can quickly genotype a large number of indel markers within a population following Mendelian segregation. The InDelGT pipeline is freely available on GitHub (https://github.com/tongchf/InDelGT).

Discovery and characterization of a new genotype of Salmonella enterica serovar Bareilly isolated from diarrhea patients of food-borne outbreaks

Since the first food-borne outbreak of Salmonella enterica serovar Bareilly in the UK (2010), it has been recognized as a new type of food-borne pathogen in S. enterica. To detect and characterize this new serovar pathogen in South Korea, a total of 175 Salmonella strains was isolated and 31 isolates were identified as S. Bareilly from various food-borne outbreaks between 2014 and 2018. While pulsed-field gel electrophoresis (PFGE) analysis using XbaI revealed two major groups (A and B) each with two subgroups (A1, A2/B1, B2), average nucleotide identity (ANI), single nucleotide polymorphism (SNP), and in silico multilocus sequence typing (MLST) analyses confirmed only two major groups. Interestingly, extended SNP analysis with 67 S. Bareilly strains from outbreaks in other countries revealed that A group strains between 2014 and 2016 shared a close evolutionary relationship with the strains from outside of South Korea; however, the B group strains in 2018 were located in a separate SNP tree branch. These findings suggest that the A group may share common ancestor with the strains of previous outbreaks in the UK or other countries, while the B group is a new genotype. Comparative virulence factor (VF) analysis between the A and B group strains showed that S. Bareilly in the B group has more various than that of the A group. A comparative biofilm formation assay supports for this, which B group strain GG-21 has higher biofilm formation activity than A group strain GG-07. Antibiotic susceptibility test of 31 S. Bareilly strains revealed high susceptibility to 17 tested antibiotics, suggesting that S. Bareilly can be easily treated by antibiotics.

Group II intron-like reverse transcriptases function in double-strand break repair

Bacteria encode reverse transcriptases (RTs) of unknown function that are closely related to group II intron-encoded RTs. We found that a Pseudomonas aeruginosa group II intron-like RT (G2L4 RT) with YIDD instead of YADD at its active site functions in DNA repair in its native host and when expressed in Escherichia coli. G2L4 RT has biochemical activities strikingly similar to those of human DNA repair polymerase θ and uses them for translesion DNA synthesis and double-strand break repair (DSBR) via microhomology-mediated end-joining (MMEJ). We also found that a group II intron RT can function similarly in DNA repair, with reciprocal active-site substitutions showing isoleucine favors MMEJ and alanine favors primer extension in both enzymes. These DNA repair functions utilize conserved structural features of non-LTR-retroelement RTs, including human LINE-1 and other eukaryotic non-LTR-retrotransposon RTs, suggesting such enzymes may have inherent ability to function in DSBR in a wide range of organisms.

Synbiotic Intervention with Lactobacilli, Bifidobacteria, and Inulin in Healthy Volunteers Increases the Abundance of Bifidobacteria but Does Not Alter Microbial Diversity

Synbiotics combine probiotics and prebiotics and are being investigated for potential health benefits. In this single-group-design trial, we analyzed changes in the gut microbiome, stool quality, and gastrointestinal well-being in 15 healthy volunteers after a synbiotic intervention comprising Lacticaseibacillus rhamnosus (LGG), Lactobacillus acidophilus (LA-5), Lacticaseibacillus paracasei subsp. paracasei (L. CASEI 431), and Bifidobacterium animalis subsp. lactis BB-12 and 20 g of chicory-derived inulin powder consumed daily for 4 weeks. Fecal samples were collected at baseline and at completion of the intervention, and all participants completed a fecal diary based on the Bristol Stool Scale and recorded their gastrointestinal well-being. No adverse effects were observed after consumption of the synbiotic product, and stool consistency and frequency remained almost unchanged during the trial. Microbiome analysis of the fecal samples was achieved using shotgun sequencing followed by taxonomic profiling. No changes in alpha and beta diversity were seen after the intervention. Greater relative abundances of Bifidobacteriaceae were observed in 12 subjects, with indigenous bifidobacteria species constituting the main increase. All four probiotic organisms increased in abundance, and L. rhamnosus, B. animalis, and L. acidophilus were differentially abundant, compared to baseline. Comparison of the fecal strains to the B. animalis subsp. lactis BB-12 reference genome and the sequenced symbiotic product revealed only a few single-nucleotide polymorphisms differentiating the probiotic B. animalis subsp. lactis BB-12 from the fecal strains identified, indicating that this probiotic strain was detectable after the intervention. IMPORTANCE The effects of probiotics/synbiotics are seldom investigated in healthy volunteers; therefore, this study is important, especially considering the safety aspects of multiple probiotics together with prebiotic fiber in consumption by humans. The study explores at the potential of a synbiotic intervention with lactobacilli, bifidobacteria, and inulin in healthy volunteers and tracks the ingested probiotic strain B. animalis subsp. lactis.

Comprehensive Genome and Plasmidome Analysis of Antimicrobial Resistant Bacteria in Wastewater Treatment Plant Effluent of Tokyo

To characterize environmental antimicrobial resistance (AMR) in urban areas, extended-spectrum β-lactamase- (ESBL)/carbapenemase-producing bacteria (EPB/CPB, respectively) from urban wastewater treatment plant effluents in Tokyo were isolated on CHROMagar ESBL plate. Complete genome sequence analysis, including plasmids, indicated that 126 CTX-M-positive isolates (31%) were identified among the 404 obtained isolates. The CTX-M-9 group was predominant (n = 65, 52%), followed by the CTX-M-1 group (n = 44, 35%). Comparative genome analysis revealed that CTX-M-27-positive E. coli O16:H5-ST131-fimH41 exhibited a stable genome structure and clonal-global dissemination. Plasmidome network analysis revealed that 304 complete plasmid sequences among 85 isolates were grouped into 14 incompatibility (Inc) network communities (Co1 to Co14). Co10 consisted of primarily IncFIA/IncFIB plasmids harboring bla_CTX-M in E. coli, whereas Co12 consisted primarily of IncFIA(HI1)/Inc FIB(K) plasmids harboring bla_CTX-M, bla_KPC, and bla_GES in Klebsiella spp. Co11 was markedly located around Co10 and Co12. Co11 exhibited bla_CTX-M, bla_KPC, and bla_NDM, and was mainly detected in E. coli and Klebsiella spp. from human and animal sources, suggesting a mutual role of Co11 in horizontal gene transfer between E. coli and Klebsiella spp. This comprehensive resistome analysis uncovers the mode of relational transfer among bacterial species, highlighting the potential source of AMR burden on public health in urban communities.

FVC as an adaptive and accurate method for filtering variants from popular NGS analysis pipelines

The quality control of variants from whole-genome sequencing data is vital in clinical diagnosis and human genetics research. However, current filtering methods (Frequency, Hard-Filter, VQSR, GARFIELD, and VEF) were developed to be utilized on particular variant callers and have certain limitations. Especially, the number of eliminated true variants far exceeds the number of removed false variants using these methods. Here, we present an adaptive method for quality control on genetic variants from different analysis pipelines, and validate it on the variants generated from four popular variant callers (GATK HaplotypeCaller, Mutect2, Varscan2, and DeepVariant). FVC consistently exhibited the best performance. It removed far more false variants than the current state-of-the-art filtering methods and recalled ~51-99% true variants filtered out by the other methods. Once trained, FVC can be conveniently integrated into a user-specific variant calling pipeline.

FVC is a method for calling specific gene variants from whole genome data, for potential use in clinical diagnosis and human genetics research.

Vaccination decreases the risk of influenza A virus reassortment but not genetic variation in pigs

Although vaccination is broadly used in North American swine breeding herds, managing swine influenza is challenging primarily due to the continuous evolution of influenza A virus (IAV) and the ability of the virus to transmit among vaccinated pigs. Studies that have simultaneously assessed the impact of vaccination on the emergence of IAV reassortment and genetic variation in pigs are limited. Here, we directly sequenced 28 bronchoalveolar lavage fluid (BALF) samples collected from vaccinated and unvaccinated pigs co-infected with H1N1 and H3N2 IAV strains, and characterized 202 individual viral plaques recovered from 13 BALF samples. We identified 54 reassortant viruses that were grouped in 17 single and 16 mixed genotypes. Notably, we found that prime-boost vaccinated pigs had less reassortant viruses than nonvaccinated pigs, likely due to a reduction in the number of days pigs were co-infected with both challenge viruses. However, direct sequencing from BALF samples revealed limited impact of vaccination on viral variant frequency, evolutionary rates, and nucleotide diversity in any IAV coding regions. Overall, our results highlight the value of IAV vaccination not only at limiting virus replication in pigs but also at protecting public health by restricting the generation of novel reassortants with zoonotic and/or pandemic potential.

Swine influenza A viruses cause severe illness among pigs and financial losses on pig farms worldwide. These viruses can also infect humans and have caused deadly human pandemics in the past. Influenza A viruses are dangerous because viruses can be transferred between humans, birds and pigs. These co-infections can allow the viruses to swap genetic material. Viral genetic exchanges can result in new virus strains that are more dangerous or that can infect other types of animals more easily.

Farmers vaccinate their pigs to control the swine influenza A virus. The vaccines are regularly updated to match circulating virus strains. But the virus evolves rapidly to escape vaccine-induced immunity, and infections are common even in vaccinated pigs. Learning about how vaccination affects the evolution of influenza A viruses in pigs could help scientists prevent outbreaks on pig farms and avoid spillover pandemics in humans.

Li et al. show that influenza A viruses are less likely to swap genetic material in vaccinated and boosted pigs than in unvaccinated animals. In the experiments, Li et al. collected swine influenza A samples from the lungs of pigs that had received different vaccination protocols. Next, Li et al. used next-generation sequencing to identify new mutations in the virus or genetic swaps among different strains. In pigs infected with both the H1N1 and H3N2 strains of influenza, the two viruses began trading genes within a week. But less genetic mixing occurred in vaccinated and boosted pigs because they spent less time infected with both viruses than in unvaccinated pigs. The vaccination status of the pig did not have much effect on how many new mutations occurred in the viruses.

The experiments show that vaccinating and boosting pigs against influenza A viruses may protect against genetic swapping among influenza viruses. If future studies on pig farms confirm the results, the information gleaned from the study could help scientists improve farm vaccine protocols to further reduce influenza risks to animals and people.

Personalized neoantigen vaccine combined with PD-1 blockade increases CD8+ tissue-resident memory T-cell infiltration in preclinical hepatocellular carcinoma models

BACKGROUND

Personalized neoantigen vaccine could induce a robust antitumor immune response in multiple cancers, whose efficacy could be further enhanced by combining with programmed cell death 1 blockade (α-PD-1). However, the corresponding immune response and synergistic mechanisms remain largely unclear. Here, we aimed to develop clinically available combinational therapeutic strategy and further explore its potential antitumor mechanisms in hepatocellular carcinoma (HCC).

METHODS

Neoantigen peptide vaccine (NeoVAC) for murine HCC cell line Hepa1-6 was developed and optimized by neoantigen screening and adjuvant optimization. Then the synergistic efficacy and related molecular mechanisms of NeoVAC combined with α-PD-1 in HCC were evaluated by orthotopic HCC mouse model, single-cell RNA sequencing, tetramer flow cytometry, immunofluorescence, etc. The tumor-killing capacity of CD8+ tissue-resident memory T cells (CD8+ TRMs) was assessed by orthotopic HCC mouse model, and autologous patient-derived cells.

RESULTS

NeoVAC, which consisted of seven high immunogenic neoantigen peptides and clinical-grade Poly(I:C), could generate a strong antitumor immune response in HCC mouse models. Significantly, its efficacy could be further improved by combining with α-PD-1, with 80% of durable tumor regression and long-term immune memory in orthotopic HCC models. Moreover, in-depth analysis of the tumor immune microenvironment showed that the percentage of CD8+ TRMs was remarkedly increased in NeoVAC plus α-PD-1 treatment group, and positively associated with the antitumor efficacy. In vitro and in vivo T-cell cytotoxicity assay further confirmed the strong tumor-killing capacity of CD8+ TRMs sorting from orthotopic mouse HCC or patient's HCC tissue.

CONCLUSIONS

This study showed that NeoVAC plus α-PD-1 could induce a strong antitumor response and long-term tumor-specific immune memory in HCC by increasing CD8+ TRMs infiltration, which might serve as a potential immune-therapeutic target for HCC.

A public antibody class recognizes an S2 epitope exposed on open conformations of SARS-CoV-2 spike

Delineating the origins and properties of antibodies elicited by SARS-CoV-2 infection and vaccination is critical for understanding their benefits and potential shortcomings. Therefore, we investigate the SARS-CoV-2 spike (S)-reactive B cell repertoire in unexposed individuals by flow cytometry and single-cell sequencing. We show that ∼82% of SARS-CoV-2 S-reactive B cells harbor a naive phenotype, which represents an unusually high fraction of total human naive B cells (∼0.1%). Approximately 10% of these naive S-reactive B cells share an IGHV1-69/IGKV3-11 B cell receptor pairing, an enrichment of 18-fold compared to the complete naive repertoire. Following SARS-CoV-2 infection, we report an average 37-fold enrichment of IGHV1-69/IGKV3-11 B cell receptor pairing in the S-reactive memory B cells compared to the unselected memory repertoire. This class of B cells targets a previously undefined non-neutralizing epitope on the S2 subunit that becomes exposed on S proteins used in approved vaccines when they transition away from the native pre-fusion state because of instability. These findings can help guide the improvement of SARS-CoV-2 vaccines.

UV and Chemically Induced Halomonas smyrnensis Mutants for Enhanced Levan Productivity

Halomonas smyrnensis AAD6^T is a moderately halophilic bacterium proven to be a powerful biotechnological tool with its ability to accumulate valuable biopolymers such as levan and poly(3-hydroxybutyrate) (PHB). Levan is a fructose homopolymer with β-2,6 fructofuranosidic linkages on the polymer backbone, and its distinctive applications in various industries such as food, pharmaceutical, medical, and chemical have been well-defined. On the other hand, PHB is a promising raw material to produce biodegradable plastics. Although it was shown in our previous studies that H. smyrnensis AAD6^T exhibits one of the highest conversion yields of sucrose to levan reported to date, novel strategies are required to overcome high costs of levan production. In this study, we aimed at increasing levan productivity of H. smyrnensis AAD6^T cultures using random mutagenesis techniques combined (i.e., ethyl methanesulfate treatment and/or ultraviolet irradiation). After several consecutive treatments, mutant strains BAE2, BAE5 and BAE6 were selected as efficient levan producers, as BAE2 standing out as the most efficient one not only in sucrose utilization and levan production rates, but also in final PHB concentrations. The mutants' whole genome sequences were analysed to determine the mutations occurred. Several mutations in genes related to central carbon metabolism and osmoregulation were found. Our results suggest that random mutagenesis can be a facile and efficient strategy to enhance the performance of extremophiles in adverse conditions.

DNA methylation-based epigenetic signatures predict somatic genomic alterations in gliomas

Molecular classification has improved diagnosis and treatment for patients with malignant gliomas. However, classification has relied on individual assays that are both costly and slow, leading to frequent delays in treatment. Here, we propose the use of DNA methylation, as an emerging clinical diagnostic platform, to classify gliomas based on major genomic alterations and provide insight into subtype characteristics. We show that using machine learning models, DNA methylation signatures can accurately predict somatic alterations and show improvement over existing classifiers. The established Unified Diagnostic Pipeline (UniD) we develop is rapid and cost-effective for genomic alterations and gene expression subtypes diagnostic at early clinical phase and improves over individual assays currently in clinical use. The significant relationship between genetic alteration and epigenetic signature indicates broad applicability of our approach to other malignancies.

No clinical assay currently exists to classify glioma tumours based on gene expression. Here, the authors develop a DNA methylation-based classifier, Unified Diagnostic Pipeline (UniD) that identifies genomic alterations and gene expression subtypes of gliomas.

SECEDO: SNV-based subclone detection using ultra-low coverage single-cell DNA sequencing

MOTIVATION

Several recently developed single-cell DNA sequencing technologies enable whole-genome sequencing of thousands of cells. However, the ultra-low coverage of the sequenced data (<0.05× per cell) mostly limits their usage to the identification of copy number alterations in multi-megabase segments. Many tumors are not copy number-driven, and thus single-nucleotide variant (SNV)-based subclone detection may contribute to a more comprehensive view on intra-tumor heterogeneity. Due to the low coverage of the data, the identification of SNVs is only possible when superimposing the sequenced genomes of hundreds of genetically similar cells. Thus, we have developed a new approach to efficiently cluster tumor cells based on a Bayesian filtering approach of relevant loci and exploiting read overlap and phasing.

RESULTS

We developed Single Cell Data Tumor Clusterer (SECEDO, lat. 'to separate'), a new method to cluster tumor cells based solely on SNVs, inferred on ultra-low coverage single-cell DNA sequencing data. We applied SECEDO to a synthetic dataset simulating 7250 cells and eight tumor subclones from a single patient and were able to accurately reconstruct the clonal composition, detecting 92.11% of the somatic SNVs, with the smallest clusters representing only 6.9% of the total population. When applied to five real single-cell sequencing datasets from a breast cancer patient, each consisting of ≈2000 cells, SECEDO was able to recover the major clonal composition in each dataset at the original coverage of 0.03×, achieving an Adjusted Rand Index (ARI) score of ≈0.6. The current state-of-the-art SNV-based clustering method achieved an ARI score of ≈0, even after merging cells to create higher coverage data (factor 10 increase), and was only able to match SECEDOs performance when pooling data from all five datasets, in addition to artificially increasing the sequencing coverage by a factor of 7. Variant calling on the resulting clusters recovered more than twice as many SNVs as would have been detected if calling on all cells together. Further, the allelic ratio of the called SNVs on each subcluster was more than double relative to the allelic ratio of the SNVs called without clustering, thus demonstrating that calling variants on subclones, in addition to both increasing sensitivity of SNV detection and attaching SNVs to subclones, significantly increases the confidence of the called variants.

AVAILABILITY AND IMPLEMENTATION

SECEDO is implemented in C++ and is publicly available at https://github.com/ratschlab/secedo. Instructions to download the data and the evaluation code to reproduce the findings in this paper are available at: https://github.com/ratschlab/secedo-evaluation. The code and data of the submitted version are archived at: https://doi.org/10.5281/zenodo.6516955.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Non-Invasive Detection of a De Novo Frameshift Variant of STAG2 in a Female Fetus: Escape Genes Influence the Manifestation of X-Linked Diseases in Females

Background: We report on a 20-week-old female fetus with a diaphragmatic hernia and other malformations, all of which appeared after the first-trimester ultrasound. Methods and Results: Whole trio exome sequencing (WES) on cell-free fetal DNA (cff-DNA) revealed a de novo frameshift variant of the X-linked STAG2 gene. Loss-of-function (LoF) STAG2 variants cause either holoprosencephaly (HPE) or Mullegama–Klein–Martinez syndrome (MKMS), are de novo, and only affect females, indicating male lethality. In contrast, missense mutations associate with milder forms of MKMS and follow the classic X-linked recessive inheritance transmitted from healthy mothers to male offspring. STAG2 has been reported to escape X-inactivation, suggesting that disease onset in LoF females is dependent on inadequate dosing for at least some of the transcripts, as is the case with a part of the autosomal dominant diseases. Missense STAG2 variants produce a quantity of transcripts, which, while resulting in a different protein, leads to disease only in hemizygous males. Similar inheritance patterns are described for other escapee genes. Conclusions: This study confirms the advantage of WES on cff-DNA and emphasizes the role of the type of the variant in X-linked disorders.

Dioecy and chromosomal sex determination are maintained through allopolyploid speciation in the plant genus Mercurialis

Polyploidization may precipitate dramatic changes to the genome, including chromosome rearrangements, gene loss, and changes in gene expression. In dioecious plants, the sex-determining mechanism may also be disrupted by polyploidization, with the potential evolution of hermaphroditism. However, while dioecy appears to have persisted through a ploidy transition in some species, it is unknown whether the newly formed polyploid maintained its sex-determining system uninterrupted, or whether dioecy re-evolved after a period of hermaphroditism. Here, we develop a bioinformatic pipeline using RNA-sequencing data from natural populations to demonstrate that the allopolyploid plant Mercurialis canariensis directly inherited its sex-determining region from one of its diploid progenitor species, M. annua, and likely remained dioecious through the transition. The sex-determining region of M. canariensis is smaller than that of its diploid progenitor, suggesting that the non-recombining region of M. annua expanded subsequent to the polyploid origin of M. canariensis. Homeologous pairs show partial sexual subfunctionalization. We discuss the possibility that gene duplicates created by polyploidization might contribute to resolving sexual antagonism.

Functional landscapes of POLE and POLD1 mutations in checkpoint blockade-dependent antitumor immunity

Defects in pathways governing genomic fidelity have been linked to improved response to immune checkpoint blockade therapy (ICB). Pathogenic POLE/POLD1 mutations can cause hypermutation, yet how diverse mutations in POLE/POLD1 influence antitumor immunity following ICB is unclear. Here, we comprehensively determined the effect of POLE/POLD1 mutations in ICB and elucidated the mechanistic impact of these mutations on tumor immunity. Murine syngeneic tumors harboring Pole/Pold1 functional mutations displayed enhanced antitumor immunity and were sensitive to ICB. Patients with POLE/POLD1 mutated tumors harboring telltale mutational signatures respond better to ICB than patients harboring wild-type or signature-negative tumors. A mutant POLE/D1 function-associated signature-based model outperformed several traditional approaches for identifying POLE/POLD1 mutated patients that benefit from ICB. Strikingly, the spectrum of mutational signatures correlates with the biochemical features of neoantigens. Alterations that cause POLE/POLD1 function-associated signatures generate T cell receptor (TCR)-contact residues with increased hydrophobicity, potentially facilitating T cell recognition. Altogether, the functional landscapes of POLE/POLD1 mutations shape immunotherapy efficacy.

Acetaldehyde makes a distinct mutation signature in single-stranded DNA

Acetaldehyde (AA), a by-product of ethanol metabolism, is acutely toxic due to its ability to react with various biological molecules including DNA and proteins, which can greatly impede key processes such as replication and transcription and lead to DNA damage. As such AA is classified as a group 1 carcinogen by the International Agency for Research on Cancer (IARC). Previous in vitro studies have shown that AA generates bulky adducts on DNA, with signature guanine-centered (GG→TT) mutations. However, due to its weak mutagenicity, short chemical half-life, and the absence of powerful genetic assays, there is considerable variability in reporting the mutagenic effects of AA in vivo. Here, we used an established yeast genetic reporter system and demonstrate that AA treatment is highly mutagenic to cells and leads to strand-biased mutations on guanines (G→T) at a high frequency on single stranded DNA (ssDNA). We further demonstrate that AA-derived mutations occur through lesion bypass on ssDNA by the translesion polymerase Polζ. Finally, we describe a unique mutation signature for AA, which we then identify in several whole-genome and -exome sequenced cancers, particularly those associated with alcohol consumption. Our study proposes a key mechanism underlying carcinogenesis by acetaldehyde—mutagenesis of single-stranded DNA.

AtNusG, a chloroplast nucleoid protein of bacterial origin linking chloroplast transcriptional and translational machineries, is required for proper chloroplast gene expression in Arabidopsis thaliana

In Escherichia coli, transcription-translation coupling is mediated by NusG. Although chloroplasts are descendants of endosymbiotic prokaryotes, the mechanism underlying this coupling in chloroplasts remains unclear. Here, we report transcription-translation coupling through AtNusG in chloroplasts. AtNusG is localized in chloroplast nucleoids and is closely associated with the chloroplast PEP complex by interacting with its essential component PAP9. It also comigrates with chloroplast ribosomes and interacts with their two components PRPS5 (uS5c) and PRPS10 (uS10c). These data suggest that the transcription and translation machineries are coupled in chloroplasts. In the atnusg mutant, the accumulation of chloroplast-encoded photosynthetic gene transcripts, such as psbA, psbB, psbC and psbD, was not obviously changed, but that of their proteins was clearly decreased. Chloroplast polysomic analysis indicated that the decrease in these proteins was due to the reduced efficiency of their translation in this mutant, leading to reduced photosynthetic efficiency and enhanced sensitivity to cold stress. These data indicate that AtNusG-mediated coupling between transcription and translation in chloroplasts ensures the rapid establishment of photosynthetic capacity for plant growth and the response to environmental changes. Therefore, our study reveals a conserved mechanism of transcription-translation coupling between chloroplasts and E. coli, which perhaps represents a regulatory mechanism of chloroplast gene expression. This study provides insights into the underlying mechanisms of chloroplast gene expression in higher plants.

Gene Conversion Explains Elevated Diversity in the Immunity Modulating APL1 Gene of the Malaria Vector Anopheles funestus

Leucine-rich repeat proteins and antimicrobial peptides are the key components of the innate immune response to Plasmodium and other microbial pathogens in Anopheles mosquitoes. The APL1 gene of the malaria vector Anopheles funestus has exceptional levels of non-synonymous polymorphism across the range of An. funestus, with an average πn of 0.027 versus a genome-wide average of 0.002, and πn is consistently high in populations across Africa. Elevated APL1 diversity was consistent between the independent pooled-template and target-enrichment datasets, however no link between APL1 diversity and insecticide resistance was observed. Although lacking the diversity of APL1, two further mosquito innate-immunity genes of the gambicin anti-microbial peptide family had πn/πs ratios greater than one, possibly driven by either positive or balancing selection. The cecropin antimicrobial peptides were expressed much more highly than other anti-microbial peptide genes, a result discordant with current models of anti-microbial peptide activity. The observed APL1 diversity likely results from gene conversion between paralogues, as evidenced by shared polymorphisms, overlapping read mappings, and recombination events among paralogues. In conclusion, we hypothesize that higher gene expression of APL1 than its paralogues is correlated with a more open chromatin formation, which enhances gene conversion and elevated diversity at this locus.

PhenGenVar: A User-Friendly Genetic Variant Detection and Visualization Tool for Precision Medicine

Precision medicine has been revolutionized by the advent of high-throughput next-generation sequencing (NGS) technology and development of various bioinformatic analysis tools for large-scale NGS big data. At the population level, biomedical studies have identified human diseases and phenotype-associated genetic variations using NGS technology, such as whole-genome sequencing, exome sequencing, and gene panel sequencing. Furthermore, patients’ genetic variations related to a specific phenotype can also be identified by analyzing their genomic information. These breakthroughs paved the way for the clinical diagnosis and precise treatment of patients’ diseases. Although many bioinformatics tools have been developed to analyze the genetic variations from the individual patient’s NGS data, it is still challenging to develop user-friendly programs for clinical physicians who do not have bioinformatics programing skills to diagnose a patient’s disease using the genomic data. In response to this demand, we developed a Phenotype to Genotype Variation program (PhenGenVar), which is a user-friendly interface for monitoring the variations in a gene of interest for molecular diagnosis. This allows for flexible filtering and browsing of variants of the disease and phenotype-associated genes. To test this program, we analyzed the whole-genome sequencing data of an anonymous person from the 1000 human genome project data. As a result, we were able to identify several genomic variations, including single-nucleotide polymorphism, insertions, and deletions in specific gene regions. Therefore, PhenGenVar can be used to diagnose a patient’s disease. PhenGenVar is freely accessible and is available at our website.

Genomic characterization of lymphomas in patients with inborn errors of immunity

Inborn errors of immunity-associated lymphomas are characterized by distinct clinical features and genetic signatures.

Both germline and somatic alterations contribute to lymphomagenesis in patients with inborn errors of immunity.

Patients with inborn errors of immunity (IEI) have a higher risk of developing cancer, especially lymphoma. However, the molecular basis for IEI-related lymphoma is complex and remains elusive. Here, we perform an in-depth analysis of lymphoma genomes derived from 23 IEI patients. We identified and validated disease-causing or -associated germline mutations in 14 of 23 patients involving ATM, BACH2, BLM, CD70, G6PD, NBN, PIK3CD, PTEN, and TNFRSF13B. Furthermore, we profiled somatic mutations in the lymphoma genome and identified 8 genes that were mutated at a significantly higher level in IEI-associated diffuse large B-cell lymphomas (DLBCLs) than in non-IEI DLBCLs, such as BRCA2, NCOR1, KLF2, FAS, CCND3, and BRWD3. The latter, BRWD3, is furthermore preferentially mutated in tumors of a subgroup of activated phosphoinositide 3-kinase δ syndrome patients. We also identified 5 genomic mutational signatures, including 2 DNA repair deficiency-related signatures, in IEI-associated lymphomas and a strikingly high number of inter- and intrachromosomal structural variants in the tumor genome of a Bloom syndrome patient. In summary, our comprehensive genomic characterization of lymphomas derived from patients with rare genetic disorders expands our understanding of lymphomagenesis and provides new insights for targeted therapy.

Analysis of recurrent molecular alterations in phyllodes tumour of breast: insights into prognosis and pathogenesis

Phyllodes tumour (PT) of breast is a rare biphasic neoplasm. Recent next generation sequencing analyses had revealed novel genetic alterations in PT but lacked a further characterisation of their relationship to different PT features and outcome. Here, using targeted sequencing, we examined a panel of 90 recurrently altered or cancer related genes in 88 PT samples (including 49 benign, 25 borderline and 14 malignant PT). Twenty-three genes showed alterations in at least 8.0% of cases. Alterations were significantly higher with an increasing grade of PT (p=0.033), particularly for copy number alterations. The top ten alterations were TERT promoter (58.0%), MED12 (53.4%), RARA (22.8%), FLNA (19.3%), SETD2 (15.9%), SYNE1 (18.2%), PCLO (15.9%), KMT2D (14.3%), CDKN2A (15.9%) and DNAH11 (14.8%). Alterations in CDKN2A/B, EGFR, TP53, PIK3CA, PTEN and ARID1B (p≤0.039) were associated with a higher grade. Analysing alterations based on common pathways indicated a significant correlation of cell cycle pathway and epigenetic alterations with a higher PT grade (p=0.036 and 0.075 respectively). Interestingly, recurrences were not correlated with tumour grade, but related to the presence of RARA mutation (p=0.011) and the absence of alterations in epigenetic pathway (p=0.031). Analysis of synchronous pair of PT showed more differences in gene mutations with divergent MED12 mutation. By contrast, the recurrent samples showed similar genetic alterations as the primary tumours. In summary, we characterised genetic alterations in PTs of different grades and confirmed the recurrent alterations observed in earlier studies. In addition, current data implicated the roles of cell cycle, epigenetic and RARA changes in PT recurrence and tumourogenesis.

Whole Genome Variant Dataset for Enriching Studies across 18 Different Cancers

Whole genome sequencing (WGS) has helped to revolutionize biology, but the computational challenge remains for extracting valuable inferences from this information. Here, we present the cancer-associated variants from the Cancer Genome Atlas (TCGA) WGS dataset. This set of data will allow cancer researchers to further expand their analysis beyond the exomic regions of the genome to the entire genome. A total of 1342 WGS alignments available from the consortium were processed with VarScan2 and deposited to the NCI Cancer Cloud. The sample set covers 18 different cancers and reveals 157,313,519 pooled (non-unique) cancer-associated single-nucleotide variations (SNVs) across all samples. There was an average of 117,223 SNVs per sample, with a range from 1111 to 775,470 and a standard deviation of 163,273. The dataset was incorporated into BigQuery, which allows for fast access and cross-mapping, which will allow researchers to enrich their current studies with a plethora of newly available genomic data.

Dysregulated mitochondrial and cytosolic tRNA m1A methylation in Alzheimer's disease

RNA modifications affect many aspects of RNA metabolism and are involved in the regulation of many different biological processes. Mono-methylation of adenosine in the N1 position, N1-methyladensoine (m1A), is a reversible modification that is known to target rRNAs and tRNAs. m1A has been shown to increase tRNA structural stability and induce correct tRNA folding. Recent studies have begun to associate the dysregulation of epitranscriptomic control with age-related disorders such as Alzheimer's disease. Here, we applied the newly developed m1A-quant-seq approach to map the brain abundant m1A RNA modification in the cortex of an Alzheimer's disease mouse model, 5XFAD. We observed hypomethylation in both mitochondrial and cytosolic tRNAs in 5XFAD mice compared with wild type. Furthermore, the main enzymes responsible for the addition of m1A in mitochondrial (TRMT10C, HSD17B10) and cytosolic tRNAs (TRMT61A) displayed decreased expression in 5XFAD compared with wild-type mice. Knockdown of these enzymes results in a more severe phenotype in a Drosophila tau model, and differential m1A methylation is correlated with differences in mature mitochondrial tRNA expression. Collectively, this work suggests that hypo m1A modification in tRNAs may play a role in Alzheimer's disease pathogenesis.

Whole mitochondrial genome analysis of the Daur ethnic minority from Hulunbuir in the Inner Mongolia Autonomous Region of China

BACKGROUND

Mitochondrial DNA (mtDNA) variations are often associated with bioenergetics, disease, and speciation and can be used to track the history of women. Although advances in massively parallel sequencing (MPS) technology have greatly promoted our understanding of the population's history (especially genome-wide data and whole Y chromosome sequencing), the whole mtDNA sequence of many important groups has not been fully studied. In this study, we employed whole mitogenomes of 209 healthy and unrelated individuals from the Daur group, a Mongolic-speaking representative population of the indigenous groups in the Heilongjiang River basin (also known as the Amur River basin).

RESULTS

The dataset presented 127 distinct mtDNA haplotypes, resulting in a haplotype diversity of 0.9933. Most of haplotypes were assigned to eastern Eurasian-specific lineages, such as D4 (19.62%), B4 (9.09%), D5 (7.66%) and M7 (4.78%). Population comparisons showed that the Daurians do have certain connections with the ancient populations in the Heilongjiang River basin but the matrilineal genetic composition of the Daur group was also greatly influenced by other non-Mongolic groups from neighboring areas.

CONCLUSIONS

Collectively, the whole mtDNA data generated in the present study will augment the existing mtDNA database. Our study provides genetic links between the Daur population and the aborigine peoples from Siberia and the Amur-Ussuri Region. But on the whole, compared with other Mongolic-speaking groups, the modern Daur population is closer to the East Asian ancestry group.

Analysis of the Plastid Genome Sequence During Maize Seedling Development

Shoot development in maize progresses from small, non-pigmented meristematic cells to expanded cells in the green leaf. During this transition, large plastid DNA (ptDNA) molecules in proplastids become fragmented in the photosynthetically-active chloroplasts. The genome sequences were determined for ptDNA obtained from Zea mays B73 plastids isolated from four tissues: base of the stalk (the meristem region); fully-developed first green leaf; first three leaves from light-grown seedlings; and first three leaves from dark-grown (etiolated) seedlings. These genome sequences were then compared to the Z. mays B73 plastid reference genome sequence that was previously obtained from green leaves. The assembled plastid genome was identical among these four tissues to the reference genome. Furthermore, there was no difference among these tissues in the sequence at and around the previously documented 27 RNA editing sites. There were, however, more sequence variants (insertions/deletions and single-nucleotide polymorphisms) for leaves grown in the dark than in the light. These variants were tightly clustered into two areas within the inverted repeat regions of the plastid genome. We propose a model for how these variant clusters could be generated by replication-transcription conflict.

Rapid Molecular Diagnosis of Genetically Inherited Neuromuscular Disorders Using Next-Generation Sequencing Technologies

Neuromuscular diseases are genetically highly heterogeneous, and differential diagnosis can be challenging. Over a 3-year period, we prospectively analyzed 268 pediatric and adult patients with a suspected diagnosis of inherited neuromuscular disorder (INMD) using comprehensive gene-panel analysis and next-generation sequencing. The rate of diagnosis increased exponentially with the addition of genes to successive versions of the INMD panel, from 31% for the first iteration (278 genes) to 40% for the last (324 genes). The global mean diagnostic rate was 36% (97/268 patients), with a diagnostic turnaround time of 4–6 weeks. Most diagnoses corresponded to muscular dystrophies/myopathies (68.37%) and peripheral nerve diseases (22.45%). The most common causative genes, TTN, RYR1, and ANO5, accounted for almost 30% of the diagnosed cases. Finally, we evaluated the utility of the differential diagnosis tool Phenomizer, which established a correlation between the phenotype and molecular findings in 21% of the diagnosed patients. In summary, comprehensive gene-panel analysis of all genes implicated in neuromuscular diseases facilitates a rapid diagnosis and provides a high diagnostic yield.

FAT1 and MSH2 Are Predictive Prognostic Markers for Chinese Osteosarcoma Patients Following Chemotherapeutic Treatment

Osteosarcoma is characterized by diverse genetic mutations, including single-nucleotide variants (SNVs), which can complicate clinical outcomes of the treatment. This study identified key mutations or polymorphisms in genes that correlate with osteosarcoma prognoses. A total of 110 patients with osteosarcoma were assigned to "good" or "poor" cohorts depending on their 5-year disease-free survival (DFS) after surgery and chemotherapeutic treatment. We performed next-generation sequencing analysis of tumor tissues for prognosis-associated SNVs in 315 tumorigenesis-related genes, followed by modeling of clinical outcomes for these patients using random forest classification via a support vector machine (SVM). Data from the Chinese Millionome Database were used to compare SNV frequency in osteosarcoma patients and healthy people. SVM screening identified 17 nonsynonymous SNVs located in 15 genes, of which rs17224367 and rs3733406 (located in MSH2 and FAT1, respectively) were strongly correlated with osteosarcoma prognosis. These results were verified in a 26-patient validation cohort, confirming that these SNVs could be used to predict prognosis. These results demonstrated that two SNVs located in MSH2 and FAT1 are associated with prognosis of osteosarcoma patients. © 2022 American Society for Bone and Mineral Research (ASBMR).

SWAN pathway-network identification of common aneuploidy-based oncogenic drivers

Haploinsufficiency drives Darwinian evolution. Siblings, while alike in many aspects, differ due to monoallelic differences inherited from each parent. In cancer, solid tumors exhibit aneuploid genetics resulting in hundreds to thousands of monoallelic gene-level copy-number alterations (CNAs) in each tumor. Aneuploidy patterns are heterogeneous, posing a challenge to identify drivers in this high-noise genetic environment. Here, we developed Shifted Weighted Annotation Network (SWAN) analysis to assess biology impacted by cumulative monoallelic changes. SWAN enables an integrated pathway-network analysis of CNAs, RNA expression, and mutations via a simple web platform. SWAN is optimized to best prioritize known and novel tumor suppressors and oncogenes, thereby identifying drivers and potential druggable vulnerabilities within cancer CNAs. Protein homeostasis, phospholipid dephosphorylation, and ion transport pathways are commonly suppressed. An atlas of CNA pathways altered in each cancer type is released. These CNA network shifts highlight new, attractive targets to exploit in solid tumors.

Gene Mutations Associated With Clinical Characteristics in the Tumors of Patients With Breast Cancer

Background

Clinical characteristics including estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor 2 (HER2) are important biomarkers in the treatment of breast cancer, but how genomic mutations affect their status is rarely studied. This study aimed at finding genomic mutations associated with these clinical characteristics.

Methods

There were 160 patients with breast cancer enrolled in this study. Samples from those patients were used for next-generation sequencing, targeting a panel of 624 pan-cancer genes. Short nucleotide mutations, copy number variations, and gene fusions were identified for each sample. Fisher’s exact test compared each pair of genes. A similarity score was constructed with the resulting P-values. Genes were clustered with the similarity scores. The identified gene clusters were compared to the status of clinical characteristics including ER, PR, HER2, and a family history of cancer (FH) in terms of the mutations in patients.

Results

Gene-by-gene analysis found that CCND1 mutations were positively correlated with ER status while ERBB2 and CDK12 mutations were positively correlated with HER2 status. Mutation-based clustering identified four gene clusters. Gene cluster 1 (ADGRA2, ZNF703, FGFR1, KAT6A, and POLB) was significantly associated with PR status; gene cluster 2 (COL1A1, AXIN2, ZNF217, GNAS, and BRIP1) and gene cluster 3 (FGF3, FGF4, FGF19, and CCND1) were significantly associated with ER status; gene cluster 2 was also negatively associated with a family history of cancer; and gene cluster 4 was significantly negatively associated with age. Patients were classified into four corresponding groups. Patient groups 1, 2, 3, and 4 had 24.1%, 36.5%, 38.7%, and 41.3% of patients with an FDA-recognized biomarker predictive of response to an FDA-approved drug, respectively.

Conclusion

This study identified genomic mutations positively associated with ER and PR status. These findings not only revealed candidate genes in ER and PR status maintenance but also provided potential treatment targets for patients with endocrine therapy resistance.

Endosomal cargo recycling mediated by Gpa1 and phosphatidylinositol 3-kinase is inhibited by glucose starvation

Cell surface protein trafficking is regulated in response to nutrient availability, with multiple pathways directing surface membrane proteins to the lysosome for degradation in response to suboptimal extracellular nutrients. Internalized protein and lipid cargoes recycle back to the surface efficiently in glucose-replete conditions, but this trafficking is attenuated following glucose starvation. We find that cells with either reduced or hyperactive phosphatidylinositol 3-kinase (PI3K) activity are defective for endosome to surface recycling. Furthermore, we find that the yeast Gα subunit Gpa1, an endosomal PI3K effector, is required for surface recycling of cargoes. Following glucose starvation, mRNA and protein levels of a distinct Gα subunit Gpa2 are elevated following nuclear translocation of Mig1, which inhibits recycling of various cargoes. As Gpa1 and Gpa2 interact at the surface where Gpa2 concentrates during glucose starvation, we propose that this disrupts PI3K activity required for recycling, potentially diverting Gpa1 to the surface and interfering with its endosomal role in recycling. In support of this model, glucose starvation and overexpression of Gpa2 alter PI3K endosomal phosphoinositide production. Glucose deprivation therefore triggers a survival mechanism to increase retention of surface cargoes in endosomes and promote their lysosomal degradation.

An Update of Cutaneous Melanoma Patients Treated in Adjuvancy With the Allogeneic Melanoma Vaccine VACCIMEL and Presentation of a Selected Case Report With In-Transit Metastases

The CSF-470 vaccine (VACCIMEL) plus BCG and GM-CSF as adjuvants has been assayed in cutaneous melanoma patients. In the adjuvant randomized Phase II study CASVAC-0401, vaccinated patients had longer distant metastasis-free survival (DMFS) than those treated with IFNα2b. Five years after locking the data, an actualization was performed. The benefit in DMFS was maintained in the vaccinated group versus the IFNα2b-treated group (p = 0.035), with a median DMFS of 96 months for VACCIMEL and 13 months for IFNα2b. The favorable risk-benefit ratio was maintained. DMFS was also analyzed as a single cohort in all the IIB, IIC, and III patients (n = 30) who had been treated with VACCIMEL. The median DMFS was 169 months, and at 48 months follow-up, it was 71.4%, which was not statistically different from DMFS of previously published results obtained in adjuvancy with ipilimumab, pembrolizumab, nivolumab, or dabrafenib/trametinib. The possible toxicity of combining VACCIMEL with anti-immune checkpoint inhibitors (ICKi) was analyzed, especially since VACCIMEL was co-adjuvated with BCG in every vaccination. A patient with in-transit metastases was studied to produce a proof of concept. During treatment with VACCIMEL, the patient developed T-cell clones reactive towards tumor-associated antigens. Three years after ending the VACCIMEL study, the patient progressed and was treated with ICKi. During ICKi treatment, the patient did not reveal any toxicity due to previous BCG treatment. When she recurred after a 4-year treatment with nivolumab, a biopsy was obtained and immunohistochemistry and RNA-seq were performed. The tumor maintained expression of tumor-associated antigens and HLA-I and immune infiltration, with immunoreactive and immunosuppressive features. VACCIMEL plus BCG and GM-CSF is an effective treatment in adjuvancy for stages IIB, IIC, and III cutaneous melanoma patients, and it is compatible with subsequent treatments with ICKi.

Consistency and reproducibility of large panel next-generation sequencing: Multi-laboratory assessment of somatic mutation detection on reference materials with mismatch repair and proofreading deficiency

INTRODUCTION

Clinical precision oncology increasingly relies on accurate genome-wide profiling using large panel next generation sequencing; however, difficulties in accurate and consistent detection of somatic mutation from individual platforms and pipelines remain an open question.

OBJECTIVES

To obtain paired tumor-normal reference materials that can be effectively constructed and interchangeable with clinical samples, and evaluate the performance of 56 panels under routine testing conditions based on the reference samples.

METHODS

Genes involved in mismatch repair and DNA proofreading were knocked down using the CRISPR-Cas9 technology to accumulate somatic mutations in a defined GM12878 cell line. They were used as reference materials to comprehensively evaluate the reproducibility and accuracy of detection results of oncopanels and explore the potential influencing factors.

RESULTS

In total, 14 paired tumor-normal reference DNA samples from engineered cell lines were prepared, and a reference dataset comprising 168 somatic mutations in a high-confidence region of 1.8 Mb were generated. For mutations with an allele frequency (AF) of more than 5% in reference samples, 56 panels collectively reported 1306 errors, including 729 false negatives (FNs), 179 false positives (FPs) and 398 reproducibility errors. The performance metric varied among panels with precision and recall ranging from 0.773 to 1 and 0.683 to 1, respectively. Incorrect and inadequate filtering accounted for a large proportion of false discovery (including FNs and FPs), while low-quality detection, cross-contamination and other sequencing errors during the wet bench process were other sources of FNs and FPs. In addition, low AF (<5%) considerably influenced the reproducibility and comparability among panels.

CONCLUSIONS

This study provided an integrated practice for developing reference standard to assess oncopanels in detecting somatic mutations and quantitatively revealed the source of detection errors. It will promote optimization, validation, and quality control among laboratories with potential applicability in clinical use.

Evaluating Antibody Mediated Protection against Alpha, Beta, and Delta SARS-CoV-2 Variants of Concern in K18-hACE2 Transgenic Mice

SARS-CoV-2 variants of concern (VoC) are impacting responses to the COVID-19 pandemic. Here, we utilized passive immunization using human convalescent plasma (HCP) obtained from a critically ill COVID-19 patient in the early pandemic to study the efficacy of polyclonal antibodies generated to ancestral SARS-CoV-2 against the Alpha, Beta, and Delta VoC in the K18 human angiotensin converting enzyme 2 (hACE2) transgenic mouse model. HCP protected mice from challenge with the original WA-1 SARS-CoV-2 strain; however, only partially protected mice challenged with the Alpha VoC (60% survival) and failed to save Beta challenged mice from succumbing to disease. HCP treatment groups had elevated receptor binding domain (RBD) and nucleocapsid IgG titers in the serum; however, Beta VoC viral RNA burden in the lung and brain was not decreased due to HCP treatment. While mice could be protected from WA-1 or Alpha challenge with a single dose of HCP, six doses of HCP could not decrease mortality of Delta challenged mice. Overall, these data demonstrate that VoC have enhanced immune evasion and this work underscores the need for in vivo models to evaluate future emerging strains. IMPORTANCE Emerging SARS-CoV-2 VoC are posing new problems regarding vaccine and monoclonal antibody efficacy. To better understand immune evasion tactics of the VoC, we utilized passive immunization to study the effect of early-pandemic SARS-CoV-2 HCP against, Alpha, Beta, and Delta VoC. We observed that HCP from a human infected with the original SARS-CoV-2 was unable to control lethality of Alpha, Beta, or Delta VoC in the K18-hACE2 transgenic mouse model of SARS-CoV-2 infection. Our findings demonstrate that passive immunization can be used as a model to evaluate immune evasion of emerging VoC strains.

A Systematic Review of Common and Brain-Disease-Specific RNA Editing Alterations Providing Novel Insights into Neurological and Neurodegenerative Disease Manifestations

RNA editing contributes to transcriptome diversification through RNA modifications in relation to genome-encoded information (RNA–DNA differences, RDDs). The deamination of Adenosine (A) to Inosine (I) or Cytidine (C) to Uridine (U) is the most common type of mammalian RNA editing. It occurs as a nuclear co- and/or post-transcriptional event catalyzed by ADARs (Adenosine deaminases acting on RNA) and APOBECs (apolipoprotein B mRNA editing enzyme catalytic polypeptide-like genes). RNA editing may modify the structure, stability, and processing of a transcript. This review focuses on RNA editing in psychiatric, neurological, neurodegenerative (NDs), and autoimmune brain disorders in humans and rodent models. We discuss targeted studies that focus on RNA editing in specific neuron-enriched transcripts with well-established functions in neuronal activity, and transcriptome-wide studies, enabled by recent technological advances. We provide comparative editome analyses between human disease and corresponding animal models. Data suggest RNA editing to be an emerging mechanism in disease development, displaying common and disease-specific patterns. Commonly edited RNAs represent potential disease-associated targets for therapeutic and diagnostic values. Currently available data are primarily descriptive, calling for additional research to expand global editing profiles and to provide disease mechanistic insights. The potential use of RNA editing events as disease biomarkers and available tools for RNA editing identification, classification, ranking, and functional characterization that are being developed will enable comprehensive analyses for a better understanding of disease(s) pathogenesis and potential cures.

Platinum-based chemotherapy for pancreatic cancer: impact of mutations in the homologous recombination repair and Fanconi anemia genes

Background

Mutations in homologous recombination (HR) and Fanconi anemia (FA) genes may predispose to pancreatic cancer (PC) and enable the prediction of sensitivity to platinum-based chemotherapy. FOLFIRINOX is a standard treatment option for non-selected PC patients and could be effective due to undiagnosed DNA repair deficiency. Here, we aimed to determine the frequency of mutations in genes involved in the HR and FA pathways, evaluate their clinical implications, and determine the objective response rate (ORR), progression-free survival (PFS), and overall survival (OS) of PC patients treated with platinum.

Methods

We performed targeted DNA sequencing of 30 genes (ABRAXAS1, ATM, ATR, BARD1, BLM, BRCA1, BRCA2, BRIP1, CDKN2A, CHEK1, CHEK2, FANCC, FANCF, FANCG, FANCI, FANCL, FANCM, MRE11A, NBN, PALB2, PTEN, RAD50, RAD51C, RAD51D, RAD52, RAD54B, RBBP8, RINT1, SLX4, and XRCC2) for 543 PC patients.

Results

In BRCA/PALB2-mutated patients with advanced PC (33 patients, 6.1%), the PFS and OS were higher for first-line platinum therapy than for non-platinum therapy [PFS: HR = 0.28, 95% confidence interval (CI) = 0.10-0.81, p = 0.02; OS: HR = 0.31, 95% CI = 0.08-1.16, p = 0.08]. Among 93 patients (17.1%) with mutations in other HR/FA genes, no statistically significant difference in PFS and OS was observed between first-line platinum therapy and non-platinum therapy (PFS: HR = 0.83, 95% CI = 0.43-1.62, p = 0.59; OS: HR = 0.58, 95% CI = 0.28-1.22, p = 0.15). For patients with early PC, no prognostic value was observed for BRCA1/2, PALB2, or other HR/FA genes mutations. Moreover, a personal history of breast, ovarian, pancreatic, or prostate cancer was identified as the only independent predictor of the risk of BRCA/PALB2 mutations (HR = 5.83, 95% CI = 2.16-15.73, p < 0.01).

Conclusion

Mutations in the BRCA1/2 and PALB2 genes increase the sensitivity of PC to platinum agents. Thus, alterations in these genes in PC patients must be determined prior to anticancer therapy.

A Phase II Study of Pembrolizumab in Combination with Capecitabine and Oxaliplatin with Molecular Profiling in Patients with Advanced Biliary Tract Carcinoma

BACKGROUND

We conducted a phase II study of the combination of pembrolizumab with capecitabine and oxaliplatin (CAPOX) in patients with advanced biliary tract carcinoma (BTC) to assess response rate and clinical efficacy. Exploratory objectives included correlative studies of immune marker expression, tumor evolution, and immune infiltration in response to treatment.

PATIENTS AND METHODS

Adult patients with histologically confirmed BTC were enrolled and received oxaliplatin and pembrolizumab on day 1 of cycles 1-6. Capecitabine was administered orally twice daily as intermittent treatment, with the first dose on day 1 and the last dose on day 14 of cycles 1-6. Starting on cycle 7, pembrolizumab monotherapy was continued until disease progression. The primary endpoint was progression-free survival (PFS). Secondary endpoints were safety, tolerability, feasibility, and response rate. Immunohistochemistry (IHC) for PD-L1 and immune infiltrates was analyzed in paired tumor biopsies, as well as bulk transcriptome and exome profiling for five patients and single-cell RNA sequencing for one partial responder.

RESULTS

Eleven patients enrolled, three of whom had received no prior systemic therapy. Treatment was well tolerated, and the most common treatment-related grade 3 or 4 adverse events were lymphocytopenia, anemia, and decreased platelet count. Three patients (27.3%) achieved a partial response, and six (54%) had stable disease. The disease control rate was 81.8%. The median PFS was 4.1 months with a 6-month PFS rate of 45.5%. Molecular profiling suggests qualitative differences in immune infiltration and clonal evolution based on response.

CONCLUSION

Capecitabine and oxaliplatin in combination with pembrolizumab is tolerable and a potentially effective treatment for refractory advanced BTC. This study highlights a design framework for the precise characterization of individual BTC tumors.

TRIAL REGISTRATION

This study was registered in ClinicalTrials.gov (NCT03111732).

Molecular signatures of antitumor neoantigen-reactive T cells from metastatic human cancers

The accurate identification of antitumor T cell receptors (TCRs) represents a major challenge for the engineering of cell-based cancer immunotherapies. By mapping 55 neoantigen-specific TCR clonotypes (NeoTCRs) from 10 metastatic human tumors to their single-cell transcriptomes, we identified signatures of CD8+ and CD4+ neoantigen-reactive tumor-infiltrating lymphocytes (TILs). Neoantigen-specific TILs exhibited tumor-specific expansion with dysfunctional phenotypes, distinct from blood-emigrant bystanders and regulatory TILs. Prospective prediction and testing of 73 NeoTCR signature-derived clonotypes demonstrated that half of the tested TCRs recognized tumor antigens or autologous tumors. NeoTCR signatures identified TCRs that target driver neoantigens and nonmutated viral or tumor-associated antigens, suggesting a common metastatic TIL exhaustion program. NeoTCR signatures delineate the landscape of TILs across metastatic tumors, enabling successful TCR prediction based purely on TIL transcriptomic states for use in cancer immunotherapy.

Tool evaluation for the detection of variably sized indels from next generation whole genome and targeted sequencing data

Insertions and deletions (indels) in human genomes are associated with a wide range of phenotypes, including various clinical disorders. High-throughput, next generation sequencing (NGS) technologies enable the detection of short genetic variants, such as single nucleotide variants (SNVs) and indels. However, the variant calling accuracy for indels remains considerably lower than for SNVs. Here we present a comparative study of the performance of variant calling tools for indel calling, evaluated with a wide repertoire of NGS datasets. While there is no single optimal tool to suit all circumstances, our results demonstrate that the choice of variant calling tool greatly impacts the precision and recall of indel calling. Furthermore, to reliably detect indels, it is essential to choose NGS technologies that offer a long read length and high coverage coupled with specific variant calling tools.

Tissue architecture delineates field cancerization in BRAFV600E-induced tumor development

Cancer cells hijack developmental growth mechanisms but whether tissue morphogenesis and architecture modify tumorigenesis is unknown. Here, we characterized a new mouse model of sporadic thyroid carcinogenesis based on inducible expression of BRAF carrying a Val600 Glu (V600E) point mutation (BRAFV600E) from the thyroglobulin promoter (TgCreERT2). Spontaneous activation of this Braf-mutant allele due to leaky activity of the Cre recombinase revealed that intrinsic properties of thyroid follicles determined BRAF-mutant cell fate. Papillary thyroid carcinomas developed multicentrically within a normal microenvironment. Each tumor originated from a single follicle that provided a confined space for growth of a distinct tumor phenotype. Lineage tracing revealed oligoclonal tumor development in infancy and early selection of BRAFV600E kinase inhibitor-resistant clones. Somatic mutations were few, non-recurrent and limited to advanced tumors. Female mice developed larger tumors than males, reproducing the gender difference of human thyroid cancer. These data indicate that BRAFV600E-induced tumorigenesis is spatiotemporally regulated depending on the maturity and heterogeneity of follicles. Moreover, thyroid tissue organization seems to determine whether a BRAF-mutant lineage becomes a cancerized lineage. The TgCreERT2;BrafCA/+ sporadic thyroid cancer mouse model provides a new tool to evaluate drug therapy at different stages of tumor evolution.