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

Identification of pathways to high-level vancomycin resistance in Clostridioides difficile that incur high fitness costs in key pathogenicity traits

Clostridioides difficile is an important human pathogen, for which there are very limited treatment options, primarily the glycopeptide antibiotic vancomycin. In recent years, vancomycin resistance has emerged as a serious problem in several gram-positive pathogens, but high-level resistance has yet to be reported for C. difficile, although it is not known if this is due to constraints upon resistance evolution in this species. Here, we show that resistance to vancomycin can evolve rapidly under ramping selection but is accompanied by fitness costs and pleiotropic trade-offs, including sporulation defects that would be expected to severely impact transmission. We identified 2 distinct pathways to resistance, both of which are predicted to result in changes to the muropeptide terminal D-Ala-D-Ala that is the primary target of vancomycin. One of these pathways involves a previously uncharacterised D,D-carboxypeptidase, expression of which is controlled by a dedicated two-component signal transduction system. Our findings suggest that while C. difficile is capable of evolving high-level vancomycin resistance, this outcome may be limited clinically due to pleiotropic effects on key pathogenicity traits. Moreover, our data identify potential mutational routes to resistance that should be considered in genomic surveillance.

Plasma ctDNA enhances the tissue-based detection of oncodriver mutations in colorectal cancer

PURPOSE

The advent of circulating tumor DNA (ctDNA) technology has provided a convenient and noninvasive means to continuously monitor cancer genomic data, facilitating personalized cancer treatment. This study aimed to evaluate the supplementary benefits of plasma ctDNA alongside traditional tissue-based next-generation sequencing (NGS) in identifying targetable mutations and tumor mutational burden (TMB) in colorectal cancers (CRC).

METHODS

Our study involved 76 CRC patients, collecting both tissue and plasma samples for NGS. We assessed the concordance of gene mutational status between ctDNA and tissue, focusing on actionable genes such as KRAS, NRAS, PIK3CA, BRAF, and ERBB2. Logistic regression analysis was used to explore variables associated with discordance and positive mutation rates.

RESULTS

In total, 26 cancer-related genes were identified. The most common variants in tumor tissues and plasma samples were in APC (57.9% vs 19.7%), TP53 (55.3% vs 22.4%) and KRAS (47.4% vs 43.4%). Tissue and ctDNA showed an overall concordance of 73.53% in detecting actionable gene mutations. Notably, plasma ctDNA improved detection for certain genes and gene pools. Variables significantly associated with discordance included gender and peritoneal metastases. TMB analysis revealed a higher detection rate in tissues compared to plasma, but combining both increased detection.

CONCLUSIONS

Our study highlights the importance of analyzing both tissue and plasma for detecting actionable mutations in CRC, with plasma ctDNA offering added value. Discordance is associated with gender and peritoneal metastases, and TMB analysis can benefit from a combination of tissue and plasma data. This approach provides valuable insights for personalized CRC treatment.

Rolling down that mountain: microgeographical adaptive divergence during a fast population expansion along a steep environmental gradient in European beech

Forest tree populations harbour high genetic diversity thanks to large effective population sizes and strong gene flow, allowing them to diversify through adaptation to local environmental pressures within dispersal distance. Many tree populations also experienced historical demographic fluctuations, including spatial population contraction or expansions at various temporal scales, which may constrain their ability to adapt to environmental variations. Our aim is to investigate how recent contraction and expansion events interfere with local adaptation, by studying patterns of adaptive divergence between closely related stands undergoing environmentally contrasted conditions, and having or not recently expanded. To investigate genome-wide signatures of local adaptation while accounting for demography, we analysed divergence in a European beech population by testing pairwise differentiation among four tree stands at ~35k Single Nucleotide Polymorphisms from ~9k genomic regions. We applied three divergence outlier search methods resting on different assumptions and targeting either single SNPs or contiguous genomic regions, while accounting for the effect of population size variations on genetic divergence. We found 27 signals of selective signatures in 19 target regions. Putatively adaptive divergence involved all stand pairs. We retrieved signals both when comparing old-growth stands and recently colonised areas and when comparing stands within the old-growth area. Therefore, adaptive divergence processes have taken place both over short time spans, under strong environmental contrasts, and over short ecological gradients, in populations that have been stable in the long term. This suggests that standing genetic variation supports local, microgeographic divergence processes, which can maintain genetic diversity at the landscape level.

GNA14 and GNAQ somatic mutations cause spinal and intracranial extra-axial cavernous hemangiomas

Extra-axial cavernous hemangiomas (ECHs) are complex vascular lesions mainly found in the spine and cavernous sinus. Their removal poses significant risk due to their vascularity and diffuse nature, and their genetic underpinnings remain incompletely understood. Our approach involved genetic analyses on 31 tissue samples of ECHs employing whole-exome sequencing and targeted deep sequencing. We explored downstream signaling pathways, gene expression changes, and resultant phenotypic shifts induced by these mutations, both in vitro and in vivo. In our cohort, 77.4% of samples had somatic missense variants in GNA14, GNAQ, or GJA4. Transcriptomic analysis highlighted significant pathway upregulation, with the GNAQ c.626A>G (p.Gln209Arg) mutation elevating PI3K-AKT-mTOR and angiogenesis-related pathways, while GNA14 c.614A>T (p.Gln205Leu) mutation led to MAPK and angiogenesis-related pathway upregulation. Using a mouse xenograft model, we observed enlarged vessels from these mutations. Additionally, we initiated rapamycin treatment in a 14-year-old individual harboring the GNAQ c.626A>G (p.Gln209Arg) variant, resulting in gradual regression of cutaneous cavernous hemangiomas and improved motor strength, with minimal side effects. Understanding these mutations and their pathways provides a foundation for developing therapies for ECHs resistant to current therapies. Indeed, the administration of rapamycin in an individual within this study highlights the promise of targeted treatments in treating these complex lesions.

Circulating tumor extracellular vesicles to monitor metastatic prostate cancer genomics and transcriptomic evolution

Extracellular vesicles (EVs) secreted by tumors are abundant in plasma, but their potential for interrogating the molecular features of tumors through multi-omic profiling remains widely unexplored. Genomic and transcriptomic profiling of circulating EV-DNA and EV-RNA isolated from in vitro and in vivo models of metastatic prostate cancer (mPC) reveal a high contribution of tumor material to EV-loaded DNA/RNA, validating the findings in two cohorts of longitudinal plasma samples collected from patients during androgen receptor signaling inhibitor (ARSI) or taxane-based therapy. EV-DNA genomic features recapitulate matched-patient biopsies and circulating tumor DNA (ctDNA) and associate with clinical progression. We develop a novel approach to enable transcriptomic profiling of EV-RNA (RExCuE). We report how the transcriptome of circulating EVs is enriched for tumor-associated transcripts, captures certain patient and tumor features, and reflects on-therapy tumor adaptation changes. Altogether, we show that EV profiling enables longitudinal transcriptomic and genomic profiling of mPC in liquid biopsy.

Genetic and epigenetic signatures for improved breeding of cultivated blueberry

Blueberry belongs to the Vaccinium genus and is a highly popular fruit crop with significant economic importance. It was not until the early twentieth century that they began to be domesticated through extensive interspecific hybridization. Here, we collected 220 Vaccinium accessions from various geographical locations, including 154 from the United States, 14 from China, eight from Australia, and 29 from Europe and other countries, comprising 164 Vaccinium corymbosum, 15 Vaccinium ashei, 10 lowbush blueberries, seven half-high blueberries, and others. We present the whole-genome variation map of 220 accessions and reconstructed the hundred-year molecular history of interspecific hybridization of blueberry. We focused on the two major blueberry subgroups, the northern highbush blueberry (NHB) and southern highbush blueberry (SHB) and identified candidate genes that contribute to their distinct traits in climate adaptability and fruit quality. Our analysis unveiled the role of gene introgression from Vaccinium darrowii and V. ashei into SHB in driving the differentiation between SHB and NHB, potentially facilitating SHB's adaptation to subtropical environments. Assisted by genome-wide association studies, our analysis suggested VcTBL44 as a pivotal gene regulator governing fruit firmness in SHB. Additionally, we conducted whole-genome bisulfite sequencing on nine NHB and 12 SHB cultivars, and characterized regions that are differentially methylated between the two subgroups. In particular, we discovered that the β-alanine metabolic pathway genes were enriched for DNA methylation changes. Our study provides high-quality genetic and epigenetic variation maps for blueberry, which offer valuable insights and resources for future blueberry breeding.

Comparative transcriptomics analysis on Senecavirus A-infected and non-infected cells

Senecavirus A (SVA) is an emerging virus that causes the vesicular disease in pigs, clinically indistinguishable from other high consequence vesicular diseases. This virus belongs to the genus Senecavirus in the family Picornaviridae. Its genome is a positive-sense, single-stranded RNA, approximately 7,300 nt in length, with a 3' poly(A) tail but without 5'-end capped structure. SVA can efficiently propagate in different cells, including some non-pig-derived cell lines. A wild-type SVA was previously rescued from its cDNA clone using reverse genetics in our laboratory. In the present study, the BSR-T7/5 cell line was inoculated with the passage-5 SVA. At 12 h post-inoculation, SVA-infected and non-infected cells were independently collected for the analysis on comparative transcriptomics. The results totally showed 628 differentially expressed genes, including 565 upregulated and 63 downregulated ones, suggesting that SVA infection significantly stimulated the transcription initiation in cells. GO and KEGG enrichment analyses demonstrated that SVA exerted multiple effects on immunity-related pathways in cells. Furthermore, the RNA sequencing data were subjected to other in-depth analyses, such as the single-nucleotide polymorphism, transcription factors, and protein-protein interactions. The present study, along with our previous proteomics and metabolomics researches, provides a multi-omics insight into the interaction between SVA and its hosts.

Dominance of African racial ancestry in honey bee colonies of Mexico 30 years after the migration of hybrids from South America

The Africanized honey bee, a hybrid of Apis mellifera scutellata from Africa with European subspecies, has been considered an invasive species and a problem for beekeeping. Africanized bees arrived in Mexico in 1986, 30 years after their accidental release in Brazil. Although government programs were implemented for its eradication, Africanized populations persist in Mexico, but precise information on the patterns of genetic introgression and racial ancestry is scarce. We determined maternal and parental racial ancestry of managed and feral honey bees across the five beekeeping regions of Mexico, using mitochondrial (mtDNA, COI-COII intergenic region) and nuclear markers (94 ancestrally informative SNPs), to assess the relationship between beekeeping management, beekeeping region, altitude, and latitude with the distribution of maternal and parental racial ancestry. Results revealed a predominantly African ancestry in the Mexican honey bees, but the proportion varied according to management, beekeeping regions, and latitude. The Mexican honey bees showed 31 haplotypes of four evolutionary lineages (A, M, C, and O). Managed honey bees had mitochondrial and nuclear higher proportions of European ancestry than feral honey bees, which had a higher proportion of African ancestry. Beekeeping regions of lower latitudes had higher proportions of African nuclear ancestry. Managed and feral honey bees showed differences in the proportion of maternal and nuclear racial ancestry. Managed honey bees from the Yucatan Peninsula and feral honey bees had a higher mtDNA than nuclear proportions of African ancestry. Managed honey bees, except those on the Yucatan Peninsula, had a higher nuclear than mtDNA proportion of African ancestry. Our study demonstrates that Africanized honey bee populations are genetically diverse and well established in Mexico, which highlights the limitations of management and government programs to contain the Africanization process and demands the incorporation of this lineage in any breeding program for sustainable beekeeping.

Artificial intelligence and neoantigens: paving the path for precision cancer immunotherapy

Cancer immunotherapy has witnessed rapid advancement in recent years, with a particular focus on neoantigens as promising targets for personalized treatments. The convergence of immunogenomics, bioinformatics, and artificial intelligence (AI) has propelled the development of innovative neoantigen discovery tools and pipelines. These tools have revolutionized our ability to identify tumor-specific antigens, providing the foundation for precision cancer immunotherapy. AI-driven algorithms can process extensive amounts of data, identify patterns, and make predictions that were once challenging to achieve. However, the integration of AI comes with its own set of challenges, leaving space for further research. With particular focus on the computational approaches, in this article we have explored the current landscape of neoantigen prediction, the fundamental concepts behind, the challenges and their potential solutions providing a comprehensive overview of this rapidly evolving field.

Segregation between an ornamental and a disease driver gene provides insights into pigment cell regulation

Genetic interactions are adaptive within a species. Hybridization can disrupt such species-specific genetic interactions and creates novel interactions that alter the hybrid progeny overall fitness. Hybrid incompatibility, which refers to degenerative genetic interactions that decrease the overall hybrid survival, is one of the results from combining two diverged genomes in hybrids. The discovery of spontaneous lethal tumorigenesis and underlying genetic interactions in select hybrids between diverged Xiphophorus species showed that lethal pathological process can result from degenerative genetic interactions. Such genetic interactions leading to lethal phenotype are thought to shield gene flow between diverged species. However, hybrids between certain Xiphophorus species do not develop such tumors. Here we report the identification of a locus residing in the genome of one Xiphophorus species that represses an oncogene from a different species. Our finding provides insights into normal and pathological pigment cell development, regulation and molecular mechanism in hybrid incompatibility.

Significance

The Dobzhansky-Muller model states epistatic interactions occurred between genes in diverged species underlies hybrid incompatibility. There are a few vertebrate interspecies hybrid cases that support the Dobzhansky-Muller model. This study reports a fish hybrid system where incompatible genetic interactions are involved in neuronal regulation of pigment cell biology, and also identified a novel point of regulation for pigment cells.

Genomic Context of SARS-CoV-2 Outbreaks in Farmed Mink in Spain during Pandemic: Unveiling Host Adaptation Mechanisms

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infects various mammalian species, with farmed minks experiencing the highest number of outbreaks. In Spain, we analyzed 67 whole genome sequences and eight spike sequences from 18 outbreaks, identifying four distinct lineages: B.1, B.1.177, B.1.1.7, and AY.98.1. The potential risk of transmission to humans raises crucial questions about mutation accumulation and its impact on viral fitness. Sequencing revealed numerous not-lineage-defining mutations, suggesting a cumulative mutation process during the outbreaks. We observed that the outbreaks were predominantly associated with different groups of mutations rather than specific lineages. This clustering pattern by the outbreaks could be attributed to the rapid accumulation of mutations, particularly in the ORF1a polyprotein and in the spike protein. Notably, the mutations G37E in NSP9, a potential host marker, and S486L in NSP13 were detected. Spike protein mutations may enhance SARS-CoV-2 adaptability by influencing trimer stability and binding to mink receptors. These findings provide valuable insights into mink coronavirus genetics, highlighting both host markers and viral transmission dynamics within communities.

What is new in FungiDB: a web-based bioinformatics platform for omics-scale data analysis for fungal and oomycete species

FungiDB (https://fungidb.org) serves as a valuable online resource that seamlessly integrates genomic and related large-scale data for a wide range of fungal and oomycete species. As an integral part of the VEuPathDB Bioinformatics Resource Center (https://veupathdb.org), FungiDB continually integrates both published and unpublished data addressing various aspects of fungal biology. Established in early 2011, the database has evolved to support 674 datasets. The datasets include over 300 genomes spanning various taxa (e.g. Ascomycota, Basidiomycota, Blastocladiomycota, Chytridiomycota, Mucoromycota, as well as Albuginales, Peronosporales, Pythiales, and Saprolegniales). In addition to genomic assemblies and annotation, over 300 extra datasets encompassing diverse information, such as expression and variation data, are also available. The resource also provides an intuitive web-based interface, facilitating comprehensive approaches to data mining and visualization. Users can test their hypotheses and navigate through omics-scale datasets using a built-in search strategy system. Moreover, FungiDB offers capabilities for private data analysis via the integrated VEuPathDB Galaxy platform. FungiDB also permits genome improvements by capturing expert knowledge through the User Comments system and the Apollo genome annotation editor for structural and functional gene curation. FungiDB facilitates data exploration and analysis and contributes to advancing research efforts by capturing expert knowledge for fungal and oomycete species.

Molnupiravir increases SARS-CoV-2 genome diversity and complexity: A case-control cohort study

Molnupiravir, an oral direct-acting antiviral effective in vitro against SARS-CoV-2, has been largely employed during the COVID-19 pandemic, since December 2021. After marketing and widespread usage, a progressive increase in SARS-CoV-2 lineages characterized by a higher transition/transversion ratio, a characteristic signature of molnupiravir action, appeared in the Global Initiative on Sharing All Influenza Data (GISAID) and International Nucleotide Sequence Database Collaboration (INSDC) databases. Here, we assessed the drug effects by SARS-CoV-2 whole-genome sequencing on 38 molnupiravir-treated persistently positive COVID-19 outpatients tested before and after treatment. Seventeen tixagevimab/cilgavimab-treated outpatients served as controls. Mutational analyses confirmed that SARS-CoV-2 exhibits an increased transition/transversion ratio seven days after initiation of molnupiravir. Moreover we observed an increased G->A ratio compared to controls, which was not related to apolipoprotein B mRNAediting enzyme, catalytic polypeptide-like (APOBEC) activity. In addition, we demonstrated for the first time an increased diversity and complexity of the viral quasispecies.

Molecular profiling of a bladder cancer with very high tumour mutational burden

The increasing incidence of urothelial bladder cancer is a notable global concern, as evidenced by the epidemiological data in terms of frequency, distribution, as well as mortality rates. Although numerous molecular alterations have been linked to the occurrence and progression of bladder cancer, currently there is a limited knowledge on the molecular signature able of accurately predicting clinical outcomes. In this report, we present a case of a pT3b high-grade infiltrating urothelial carcinoma with areas of squamous differentiation characterized by very high tumor mutational burden (TMB), with up-regulations of immune checkpoints. The high TMB, along with elevated expressions of PD-L1, PD-L2, and PD1, underscores the rationale for developing a personalized immunotherapy focused on the use of immune-checkpoint inhibitors. Additionally, molecular analysis revealed somatic mutations in several other cancer-related genes, including TP53, TP63 and NOTCH3. Mutations of TP53 and TP63 genes provide mechanistic insights on the molecular mechanisms underlying disease development and progression. Notably, the above-mentioned mutations and the elevated hypoxia score make the targeting of p53 and/or hypoxia related pathways a plausible personalized medicine option for this bladder cancer, particularly in combination with immunotherapy. Our data suggest a requirement for molecular profiling in bladder cancer to possibly select appropriate immune-checkpoint therapy.

A diet of oxidative stress-adapted bacteria improves stress resistance and lifespan in C. elegans via p38-MAPK

Organisms across taxa face stresses including variable temperature, redox imbalance, and xenobiotics. Successfully responding to stress and restoring homeostasis are crucial for survival. Aging is associated with a decreased stress response and alterations in the microbiome, which contribute to disease development. Animals and their microbiota share their environment; however, microbes have short generation time and can rapidly evolve and potentially affect host physiology during stress. Here, we leverage Caenorhabditis elegans and its simplified bacterial diet to demonstrate how microbial adaptation to oxidative stress affects the host's lifespan and stress response. We find that worms fed stress-evolved bacteria exhibit enhanced stress resistance and an extended lifespan. Through comprehensive genetic and metabolic analysis, we find that iron in stress-evolved bacteria enhances worm stress resistance and lifespan via activation of the mitogen-activated protein kinase pathway. In conclusion, our study provides evidence that understanding microbial stress-mediated adaptations could be used to slow aging and alleviate age-related health decline.

Mutational Signatures in Wild Type Escherichia coli Strains Reveal Predominance of DNA Polymerase Errors

While mutational processes operating in the Escherichia coli genome have been revealed by multiple laboratory experiments, the contribution of these processes to accumulation of bacterial polymorphism and evolution in natural environments is unknown. To address this question, we reconstruct signatures of distinct mutational processes from experimental data on E. coli hypermutators, and ask how these processes contribute to differences between naturally occurring E. coli strains. We show that both mutations accumulated in the course of evolution of wild-type strains in nature and in the lab-grown nonmutator laboratory strains are explained predominantly by the low fidelity of DNA polymerases II and III. By contrast, contributions specific to disruption of DNA repair systems cannot be detected, suggesting that temporary accelerations of mutagenesis associated with such disruptions are unimportant for within-species evolution. These observations demonstrate that accumulation of diversity in bacterial strains in nature is predominantly associated with errors of DNA polymerases.

How much can ChatGPT really help computational biologists in programming?

ChatGPT, a recently developed product by openAI, is successfully leaving its mark as a multi-purpose natural language based chatbot. In this paper, we are more interested in analyzing its potential in the field of computational biology. A major share of work done by computational biologists these days involve coding up bioinformatics algorithms, analyzing data, creating pipelining scripts and even machine learning modeling and feature extraction. This paper focuses on the potential influence (both positive and negative) of ChatGPT in the mentioned aspects with illustrative examples from different perspectives. Compared to other fields of computer science, computational biology has (1) less coding resources, (2) more sensitivity and bias issues (deals with medical data), and (3) more necessity of coding assistance (people from diverse background come to this field). Keeping such issues in mind, we cover use cases such as code writing, reviewing, debugging, converting, refactoring, and pipelining using ChatGPT from the perspective of computational biologists in this paper.

In silico search for and analysis of R gene variation in primitive cultivated potato species

Pathogen recognition receptors encoded by R genes play a key role in plant protection. Nowadays, R genes are a basis for breeding many crops, including potato. Many potato R genes have been discovered and found suitable for breeding thanks to the studies of a wide variety of wild potato species. The use of primitive cultivated potato species (PCPS) as representatives of the primary gene pool can also be promising in this respect. PCPS are the closest to the early domesticated forms of potato; therefore, their investigation could help understand the evolution of R genes. The present study was aimed at identifying and analyzing R genes in PCPS listed in the open database of NCBI and Solomics DB. In total, the study involved 27 accessions belonging to three species: Solanum phureja Juz. & Bukasov, S. stenotomum Juz. & Bukasov and S. goniocalyx Juz. & Bukasov Materials for the analysis were the sequencing data for the said three species from the PRJNA394943 and PRJCA006011 projects. An in silico search was carried out for sequences homologous to 26 R genes identified in potato species differing in phylogenetic distance from PCPS, namely nightshade (S. americanum), North- (S. bulbocastanum, S. demissum) and South-American (S. venturii, S. berthaultii) wild potato species, as well as the cultivated potato species S. tuberosum and S. andigenum. Homologs of all investigated protein-coding sequences were discovered in PCPS with a relatively high degree of similarity (85-100 %). Homologs of the Rpi-R3b, Rpi-amr3 and Rpi-ber1 genes have been identified in PCPS for the first time. An analysis of polymorphism of nucleotide and amino acid sequences has been carried out for 15 R genes. The differences in frequencies of substitutions in PCPS have been demonstrated by analysis of R genes, the reference sequences of which have been identified in different species. For all the studied NBS-LRR genes, the proportion of substituted amino acids in the LRR domain exceeds this figure for the NBS domain. The potential prospects of using PCPS as sources of resistance to Verticillium wilt have been shown.

Discrepancy of flowering time between genetically close sublineages of Aegilops umbellulata Zhuk

Aegilops umbellulata Zhuk., a wild diploid wheat-related species, has been used as a genetic resource for several important agronomic traits. However, its genetic variations have not been comprehensively studied. We sequenced RNA from 114 accessions of Ae. umbellulata to evaluate DNA polymorphisms and phenotypic variations. Bayesian clustering and phylogenetic analysis based on SNPs detected by RNA sequencing revealed two divergent lineages, UmbL1 and UmbL2. The main differences between them were in the sizes of spikes and spikelets, and culm diameter. UmbL1 is divided into two sublineages, UmbL1e and UmbL1w. These genetic differences corresponded to geographic distributions. UmbL1e, UmbL1w, and UmbL2 are found in Turkey, Iran/Iraq, and Greece, respectively. Although UmbL1e and UmbL1w were genetically similar, flowering time and other morphological traits were more distinct between these sublineages than those between the lineages. This discrepancy can be explained by the latitudinal and longitudinal differences in habitats. Specifically, latitudinal clines of flowering time were clearly observed in Ae. umbellulata, strongly correlated with solar radiation in the winter season. This observation implies that latitudinal differences are a factor in differences in the flowering times of Ae. umbellulata. Differences in flowering time could influence other morphological differences and promote genetic divergence between sublineages.

Genetically driven predisposition leads to an unusually genomic unstable renal cell carcinoma

Renal cell carcinoma originates from the lining of the proximal convoluted renal tubule and represents the most common type of kidney cancer. Risk factors and comorbidities might be associated to renal cell carcinoma, while a small fraction of 2-3% emerges from patients with predisposing cancer syndromes, typically associated to hereditary mutations in VHL, folliculin, fumarate hydratase or MET genes. Here, we report a case of renal cell carcinoma in patient with concurrent germline mutations in BRCA1 and RAD51 genes. This case displays an unusual high mutational burden and chromosomal aberrations compared to the typical profile of renal cell carcinoma. Mutational analysis on whole genome sequencing revealed an enrichment of the MMR2 mutational signature, which is indicative of impaired DNA repair capacity. Overall, the tumor displayed a profile of unusual high genomic instability which suggests a possible origin from germline predisposing mutations in the DNA repair genes BRCA1 and RAD51. While BRCA1 and RAD51 germline mutations are well-characterised in breast and ovarian cancer, their role in renal cell carcinoma is still largely unexplored. The genomic instability detected in this case of renal cell carcinoma, along with the presence of unusual mutations, might offer support to clinicians for the development of patient-tailored therapies.

Two complete 1918 influenza A/H1N1 pandemic virus genomes characterized by next-generation sequencing using RNA isolated from formalin-fixed, paraffin-embedded autopsy lung tissue samples along with evidence of secondary bacterial co-infection

The 1918 influenza pandemic was the most devastating respiratory pandemic in modern human history, with 50-100 million deaths worldwide. Here, we characterized the complete genomes of influenza A virus (IAV) from two fatal cases during the fall wave of 1918 influenza A (H1N1) pandemic in the United States, one from Walter Reed Army Hospital in Washington, DC, and the other from Camp Jackson, SC. The two complete IAV genomes were obtained by combining Illumina deep sequencing data from both total RNA and influenza viral genome-enriched libraries along with Sanger sequencing data from PCR across the sequencing gaps. This study confirms the previously reported 1918 IAV genomes and increases the total number of available complete or near-complete influenza viral genomes of the 1918 pandemic from four to six. Sequence comparisons among them confirm that the genomes of the 1918 pandemic virus were highly conserved during the main wave of the pandemic with geographic separation in North America and Europe. Metagenomic analyses revealed bacterial co-infections in both cases. Interestingly, in the Washington, DC, case, evidence is presented of the first reported Rhodococcus-influenza virus co-infection.

IMPORTANCE

This study applied modern molecular biotechnology and high-throughput sequencing to formalin-fixed, paraffin-embedded autopsy lung samples from two fatal cases during the fall wave of the 1918 influenza A (H1N1) pandemic in the United States. Complete influenza genomes were obtained from both cases, which increases the total number of available complete or near-complete influenza genomes of the 1918 pandemic virus from four to six. Sequence analysis confirms that the 1918 pandemic virus was highly conserved during the main wave of the pandemic with geographic separation in North America and Europe. Metagenomic analyses revealed bacterial co-infections in both cases, including the first reported evidence of Rhodococcus-influenza co-infection. Overall, this study offers a detailed view at the molecular level of the very limited samples from the most devastating influenza pandemic in modern human history.

Founder pathogenic variants in colorectal neoplasia susceptibility genes in Ashkenazi Jews undergoing colonoscopy

BACKGROUND

Colorectal neoplasia is one of the most common tumors affecting Western populations.

METHODS

In this study we used a custom amplicon sequencing platform and an in-house bioinformatic pipeline to study constitutional DNA from two different case series of Ashkenazi Jews undergoing colonoscopy (n = 765). The first series all had pathologically confirmed colorectal adenomas and/or carcinoma. The second series consisted of persons who had undergone a colonoscopy within the five years prior to ascertainment, regardless of findings. Ninety-one percent of all patients were asymptomatic at the time of colonoscopy.

RESULTS

In the first group (n = 438), we identified 65 founder variants (56 in APC, 2 in GREM1, 3 in MSH2 and 4 in BLM). In the second group (n = 327), the findings were 30, nothing, 1 and 1, respectively, as well as 2 MSH6 variants.

CONCLUSIONS

Overall, we found that 10 to 15% of Ashkenazi Jewish persons undergoing colonoscopy harbor variants of interest in colorectal and/or polyposis predisposition. This includes pathogenic variants in MSH6, which is associated with colorectal cancer but not with polyposis. We identified no pathogenic variants in more recently discovered polyposis predisposition genes (POLE, POLD1 or NTHL1), rendering the presence of such founder variants rare.

Detecting Somatic Mutations Without Matched Normal Samples Using Long Reads

DNA sequencing of tumours to identify somatic mutations has become a critical tool to guide the type of treatment given to cancer patients. The gold standard for mutation calling is comparing sequencing data from the tumour to a matched normal sample to avoid mis-classifying inherited SNPs as mutations. This procedure works extremely well, but in certain situations only a tumour sample is available. While approaches have been developed to find mutations without a matched normal, they have limited accuracy or require specific types of input data (e.g. ultra-deep sequencing). Here we explore the application of single molecule long read sequencing to calling somatic mutations without matched normal samples. We develop a simple theoretical framework to show how haplotype phasing is an important source of information for determining whether a variant is a somatic mutation. We then use simulations to assess the range of experimental parameters (tumour purity, sequencing depth) where this approach is effective. These ideas are developed into a prototype somatic mutation caller, smrest, and its use is demonstrated on two highly mutated cancer cell lines. Finally, we argue that this approach has potential to measure clinically important biomarkers that are based on the genome-wide distribution of mutations: tumour mutation burden and mutation signatures.

Genetic mutation in Escherichia coli genome during adaptation to the murine intestine is optimized for the host diet

Mammalian gut microbes colonize the intestinal tract of their host and adapt to establish a microbial ecosystem. The host diet changes the nutrient profile of the intestine and has a high impact on microbiota composition. Genetic mutations in Escherichia coli, a prevalent species in the human gut, allow for adaptation to the mammalian intestine, as reported in previous studies. However, the extent of colonization fitness in the intestine elevated by genetic mutation and the effects of diet change on these mutations in E. coli are still poorly known. Here, we show that notable mutations in sugar metabolism-related genes (gatC, araC, and malI) were detected in the E. coli K-12 genome just 2 weeks after colonization in the germ-free mouse intestine. In addition to elevated fitness by deletion of gatC, as previously reported, deletion of araC and malI also elevated E. coli fitness in the murine intestine in a host diet-dependent manner. In vitro cultures of medium containing nutrients abundant in the intestine (e.g., galactose, N-acetylglucosamine, and asparagine) also showed increased E. coli fitness after deletion of the genes-of-interest associated with their metabolism. Furthermore, the host diet was found to influence the developmental trajectory of gene mutations in E. coli. Taken together, we suggest that genetic mutations in E. coli are selected in response to the intestinal environment, which facilitates efficient utilization of nutrients abundant in the intestine under laboratory conditions. Our study offers some insight into the possible adaptation mechanisms of gut microbes.IMPORTANCEThe gut microbiota is closely associated with human health and is greatly impacted by the host diet. Bacteria such as Escherichia coli live in the gut all throughout the life of a human host and adapt to the intestinal environment. Adaptive mutations in E. coli are reported to enhance fitness in the mammalian intestine, but to what extent is still poorly known. It is also unknown whether the host diet affects what genes are mutated and to what extent fitness is affected. This study suggests that genetic mutations in the E. coli K-12 strain are selected in response to the intestinal environment and facilitate efficient utilization of abundant nutrients in the germ-free mouse intestine. Our study provides a better understanding of these intestinal adaptation mechanisms of gut microbes.

Computational investigations into structure and function impact of novel mutations identified in targeted exons from ovarian cancer cell lines

The lack of sensitive and specific biomarkers for ovarian cancer leads to late stage diagnosis of the disease in a majority of the cases. Mutation accumulation is the basis for cancer progression, thus identifying mutations is an important step in the disease diagnosis. In the present study, a comprehensive analysis of fifteen Next Generation Sequencing samples from thirteen ovarian cancer cell lines was carried out for the identification of new mutations. The study revealed eight clinically significant novel mutations in six ovarian cancer oncogenes, viz. SMARCA4, ARID1A, PPP2R1A, CTNNB1, DICER1 and PIK3CA. In-depth computational analysis revealed that the mutations affected the structure of the proteins in terms of stability, solvent accessible surface area and molecular dynamics. Moreover, the mutations were present in functionally significant domains of the proteins, thereby adversely affecting the protein functionality. PPI network for SMARCA4, CTNNB1, DICER1, PIK3CA, PPP2R1A and ARID1A showed that these genes were involved in certain significant pathways affecting various hallmarks of cancer. For further validation, in vitro studies were performed that revealed hypermutability of the CTNNB1 gene. Through this study we have identified some key mutations and have analysed their structural and functional impact. The study establishes some key mutations, which can be potentially explored as biomarker and drug target.

Prophage enhances the ability of deep-sea bacterium Shewanella psychrophila WP2 to utilize D-amino acid

Prophages are prevalent in the marine bacterial genomes and reshape the physiology and metabolism of their hosts. However, whether and how prophages influence the microbial degradation of D-amino acids (D-AAs), which is one of the widely distributed recalcitrant dissolved organic matters (RDOMs) in the ocean, remain to be explored. In this study, we addressed this issue in a representative marine bacterium, Shewanella psychrophila WP2 (WP2), and its integrated prophage SP1. Notably, compared to the WP2 wild-type strain, the SP1 deletion mutant of WP2 (WP2ΔSP1) exhibited a significantly lower D-glutamate (D-Glu) consumption rate and longer lag phase when D-Glu was used as the sole nitrogen source. The subsequent transcriptome analysis identified 1,523 differentially expressed genes involved in diverse cellular processes, especially that multiple genes related to inorganic nitrogen metabolism were highly upregulated. In addition, the dynamic profiles of ammonium, nitrate, and nitrite were distinct between the culture media of WP2 and WP2ΔSP1. Finally, we provide evidence that SP1 conferred a competitive advantage to WP2 when D-Glu was used as the sole nitrogen source and SP1-like phages may be widely distributed in the global ocean. Taken together, these findings offer novel insight into the influences of prophages on host metabolism and RDOM cycling in marine environments.IMPORTANCEThis work represents the first exploration of the impact of prophages on the D-amino acid (D-AA) metabolism of deep-sea bacteria. By using S. psychrophila WP2 and its integrated prophage SP1 as a representative system, we found that SP1 can significantly increase the catabolism rate of WP2 to D-glutamate and produce higher concentrations of ammonium, resulting in faster growth and competitive advantages. Our findings not only deepen our understanding of the interaction between deep-sea prophages and hosts but also provide new insights into the ecological role of prophages in refractory dissolved organic matter and the nitrogen cycle in deep oceans.

Disruption of the TP53 locus in osteosarcoma leads to TP53 promoter gene fusions and restoration of parts of the TP53 signalling pathway

TP53 is the most frequently mutated gene in human cancer. This gene shows not only loss-of-function mutations but also recurrent missense mutations with gain-of-function activity. We have studied the primary bone malignancy osteosarcoma, which harbours one of the most rearranged genomes of all cancers. This is odd since it primarily affects children and adolescents who have not lived the long life thought necessary to accumulate massive numbers of mutations. In osteosarcoma, TP53 is often disrupted by structural variants. Here, we show through combined whole-genome and transcriptome analyses of 148 osteosarcomas that TP53 structural variants commonly result in loss of coding parts of the gene while simultaneously preserving and relocating the promoter region. The transferred TP53 promoter region is fused to genes previously implicated in cancer development. Paradoxically, these erroneously upregulated genes are significantly associated with the TP53 signalling pathway itself. This suggests that while the classical tumour suppressor activities of TP53 are lost, certain parts of the TP53 signalling pathway that are necessary for cancer cell survival and proliferation are retained. In line with this, our data suggest that transposition of the TP53 promoter is an early event that allows for a new normal state of genome-wide rearrangements in osteosarcoma. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Bioinformatics of germline variant discovery for rare disease diagnostics: current approaches and remaining challenges

Next-generation sequencing (NGS) has revolutionized the field of rare disease diagnostics. Whole exome and whole genome sequencing are now routinely used for diagnostic purposes; however, the overall diagnosis rate remains lower than expected. In this work, we review current approaches used for calling and interpretation of germline genetic variants in the human genome, and discuss the most important challenges that persist in the bioinformatic analysis of NGS data in medical genetics. We describe and attempt to quantitatively assess the remaining problems, such as the quality of the reference genome sequence, reproducible coverage biases, or variant calling accuracy in complex regions of the genome. We also discuss the prospects of switching to the complete human genome assembly or the human pan-genome and important caveats associated with such a switch. We touch on arguably the hardest problem of NGS data analysis for medical genomics, namely, the annotation of genetic variants and their subsequent interpretation. We highlight the most challenging aspects of annotation and prioritization of both coding and non-coding variants. Finally, we demonstrate the persistent prevalence of pathogenic variants in the coding genome, and outline research directions that may enhance the efficiency of NGS-based disease diagnostics.

Empirical Bayes single nucleotide variant-calling for next-generation sequencing data

One of the fundamental computational problems in cancer genomics is the identification of single nucleotide variants (SNVs) from DNA sequencing data. Many statistical models and software implementations for SNV calling have been developed in the literature, yet, they still disagree widely on real datasets. Based on an empirical Bayesian approach, we introduce a local false discovery rate (LFDR) estimator for germline SNV calling. Our approach learns model parameters without prior information, and simultaneously accounts for information across all sites in the genomic regions of interest. We also propose another LFDR-based algorithm that reliably prioritizes a given list of mutations called by any other variant-calling algorithm. We use a suite of gold-standard cell line data to compare our LFDR approach against a collection of widely used, state of the art programs. We find that our LFDR approach approximately matches or exceeds the performance of all of these programs, despite some very large differences among them. Furthermore, when prioritizing other algorithms' calls by our LFDR score, we find that by manipulating the type I-type II tradeoff we can select subsets of variant calls with minimal loss of sensitivity but dramatic increases in precision.

Genomic and transcriptomic profiling of hepatocellular carcinoma reveals a rare molecular subtype

Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related death worldwide, occurring predominantly in patients with underlying chronic liver disease and cirrhosis. Here, we describe a case of a 62-year-old man that was admitted to our hospital and diagnosed with HCC where the cancer has already metastasized to the retroperitoneum and peritoneum. In order to better characterize the HCC, both the cancerous liver tissue and the adjacent normal liver tissue of the patient were collected and subjected to a genomic, transcriptomic and proteomic analysis. Our patient carries a highly mutated HCC, which is characterized by both somatic mutation in the following genes ALK, CDK6, TP53, PGR. In addition, we observe several molecular alterations that are associated with potential therapy resistance, for example the expression of the organic-anion-transporting polypeptide (OATP) family members B1 and B3, that mediate the transport of the anticancer drugs, has been found decreased. Overall, our molecular profiling potentially classify the patient with poor prognosis and possibly displaying resistance to pharmacological therapy.

Aluminium bioaccumulation in colon cancer, impinging on epithelial-mesenchymal-transition and cell death

We investigated the presence of aluminium (Al) in human colon cancer samples and its potential association with biological processes involved in cancer progression, such as epithelial to mesenchymal transition (EMT) and cell death. 25 consecutive colon samples were collected from patients undergoing colonic resection. Both neoplastic and normal mucosa were collected from each patient and subjected to histological, ultrastructural and immunohistochemical analyses. Moreover, colon samples from two Al-positive patients underwent multi-omic analyses, including whole genome sequencing and RNA sequencing (RNAseq). Morin staining, used to identify in situ aluminium bioaccumulation, showed the presence of Al in tumor areas of 24 % of patients. Transmission electron microscopy and energy-dispersive X-ray microanalysis confirmed the presence of Al specifically in intracytoplasmic electrondense nanodeposits adjacent to mitochondria of colon cancer cells. Immunohistochemical analyses for vimentin and nuclear β-catenin were performed to highlight the occurrence of the EMT phenomenon in association to Al bioaccumulation. Al-positive samples showed a significant increase in both the number of vimentin-positive and nuclear β-catenin-positive cancer cells compared to Al-negative samples. Moreover, Al-positive samples exhibited a significant decrease in the number of apoptotic cells, as well as the expression of the anti-apoptotic molecule BCL-2. Multi-omic analyses revealed a higher tumor mutational burden (TMB) in Al-positive colon cancers (n = 2) compared to a control cohort (n = 100). Additionally, somatic mutations in genes associated with EMT (GATA3) and apoptosis (TP53) were observed in Al-positive colon cancers. In conclusion, this study provides the first evidence of Al bioaccumulation in colon cancer and its potential role in modulating molecular pathways involved in cancer progression, such as EMT and apoptosis. Understanding the molecular mechanisms underlying Al toxicity might contribute to improve strategies for prevention, early detection, and targeted therapies for the management of colon cancer patients.

Multiple DNA repair pathways prevent acetaldehyde-induced mutagenesis in yeast

Acetaldehyde is the primary metabolite of alcohol and is present in many environmental sources including tobacco smoke. Acetaldehyde is genotoxic, whereby it can form DNA adducts and lead to mutagenesis. Individuals with defects in acetaldehyde clearance pathways have increased susceptibility to alcohol-associated cancers. Moreover, a mutation signature specific to acetaldehyde exposure is widespread in alcohol and smoking-associated cancers. However, the pathways that repair acetaldehyde-induced DNA damage and thus prevent mutagenesis are vaguely understood. Here, we used Saccharomyces cerevisiae to systematically delete genes in each of the major DNA repair pathways to identify those that alter acetaldehyde-induced mutagenesis. We found that deletion of the nucleotide excision repair (NER) genes, RAD1 or RAD14, led to an increase in mutagenesis upon acetaldehyde exposure. Acetaldehyde-induced mutations were dependent on translesion synthesis as well as DNA inter-strand crosslink (ICL) repair in Δrad1 strains. Moreover, whole genome sequencing of the mutated isolates demonstrated an increase in C→A changes coupled with an enrichment of gCn→A changes in the acetaldehyde-treated Δrad1 isolates. The gCn→A mutation signature has been shown to be diagnostic of acetaldehyde exposure in yeast and in human cancers. We also demonstrated that the deletion of the two DNA-protein crosslink (DPC) repair proteases, WSS1 and DDI1, also led to increased acetaldehyde-induced mutagenesis. Defects in base excision repair (BER) led to a mild increase in mutagenesis, while defects in mismatch repair (MMR), homologous recombination repair (HR) and post replicative repair pathways did not impact mutagenesis upon acetaldehyde exposure. Our results in yeast were further corroborated upon analysis of whole exome sequenced liver cancers, wherein, tumors with defects in ERCC1 and ERCC4 (NER), FANCD2 (ICL repair) or SPRTN (DPC repair) carried a higher gCn→A mutation load than tumors with no deleterious mutations in these genes. Our findings demonstrate that multiple DNA repair pathways protect against acetaldehyde-induced mutagenesis.

Application of Pathogen Genomics to Outbreak Investigation

Outbreaks are a risk to public health particularly when pathogenic, hypervirulent, and/or multidrug-resistant organisms (MDROs) are involved. In a hospital setting, vulnerable populations such as the immunosuppressed, intensive care patients, and neonates are most at risk. Rapid and accurate outbreak detection is essential to implement effective interventions in clinical areas to control and stop further transmission. Advances in the field of whole genome sequencing (WGS) have resulted in lowered costs, increased capacity, and improved reproducibility of results. WGS now has the potential to revolutionize the investigation and management of outbreaks replacing conventional genotyping and other discrimination systems. Here, we outline specific procedures and protocols to implement WGS into investigation of outbreaks in healthcare settings.

Clostridioides difficile canonical L,D-transpeptidases catalyze a novel type of peptidoglycan cross-links and are not required for beta-lactam resistance

Clostridioides difficile is the leading cause of antibiotic-associated diarrhea worldwide with significant morbidity and mortality. This organism is naturally resistant to several beta-lactam antibiotics that inhibit the polymerization of peptidoglycan, an essential component of the bacteria cell envelope. Previous work has revealed that C. difficile peptidoglycan has an unusual composition. It mostly contains 3-3 cross-links, catalyzed by enzymes called L,D-transpeptidases (Ldts) that are poorly inhibited by beta-lactams. It was therefore hypothesized that peptidoglycan polymerization by these enzymes could underpin antibiotic resistance. Here, we investigated the catalytic activity of the three canonical Ldts encoded by C. difficile (LdtCd1, LdtCd2, and LdtCd3) in vitro and explored their contribution to growth and antibiotic resistance. We show that two of these enzymes catalyze the formation of novel types of peptidoglycan cross-links using meso-diaminopimelic acid both as a donor and an acceptor, also observed in peptidoglycan sacculi. We demonstrate that the simultaneous deletion of these three genes only has a minor impact on both peptidoglycan structure and resistance to beta-lactams. This unexpected result therefore implies that the formation of 3-3 peptidoglycan cross-links in C. difficile is catalyzed by as yet unidentified noncanonical Ldt enzymes.

A life-history allele of large effect shortens developmental time in a wild insect population

Developmental time is a key life-history trait with large effects on Darwinian fitness. In many insects, developmental time is currently under strong selection to minimize ecological mismatches in seasonal timing induced by climate change. The genetic basis of responses to such selection, however, is poorly understood. To address this problem, we set up a long-term evolve-and-resequence experiment in the beetle Tribolium castaneum and selected replicate, outbred populations for fast or slow embryonic development. The response to this selection was substantial and embryonic developmental timing of the selection lines started to diverge during dorsal closure. Pooled whole-genome resequencing, gene expression analysis and an RNAi screen pinpoint a 222 bp deletion containing binding sites for Broad and Tramtrack upstream of the ecdysone degrading enzyme Cyp18a1 as a main target of selection. Using CRISPR/Cas9 to reconstruct this allele in the homogenous genetic background of a laboratory strain, we unravel how this single deletion advances the embryonic ecdysone peak inducing dorsal closure and show that this allele accelerates larval development but causes a trade-off with fecundity. Our study uncovers a life-history allele of large effect and reveals the evolvability of developmental time in a natural insect population.

Osteosarcomas With Few Chromosomal Alterations or Adult Onset Are Genetically Heterogeneous

Osteosarcoma is the most common primary bone malignancy, often detected in children and adolescents and commonly associated with TP53 alterations along with a high number of chromosomal rearrangements. However, osteosarcoma can affect patients of any age, and some tumors display less genetic complexity. Besides TP53 variants, data on key driving mutations are lacking for many osteosarcomas, particularly those affecting adults. To detect osteosarcoma-specific alterations, we screened transcriptomic and genomic sequencing and copy number data from 150 bone tumors originally diagnosed as osteosarcomas. To increase the precision in gene fusion detection, we developed a bioinformatic tool denoted as NAFuse, which extracts gene fusions that are verified at both the genomic and transcriptomic levels. Apart from the already reported genetic subgroups of osteosarcoma with TP53 structural variants, or MDM2 and/or CDK4 amplification, we did not identify any recurrent genetic driver that signifies the remaining cases. Among the plethora of mutations identified, we found genetic alterations characteristic of, or similar to, those of other bone and soft tissue tumors in 8 cases. These mutations were found in tumors with relatively few other genetic alterations or in adults. Due to the lack of clinical context and available tissue, we can question the diagnosis only on a genetic basis. However, our findings support the notion that osteosarcomas with few chromosomal alterations or adult onset seem genetically distinct from conventional osteosarcomas of children and adolescents.

Zinc effects on bacteria: insights from Escherichia coli by multi-omics approach

IMPORTANCE

A long-term exposure of bacteria to zinc oxide and zinc oxide nanoparticles leads to major alterations in bacterial morphology and physiology. These included biochemical and physiological processes promoting the emergence of strains with multi-drug resistance and virulence traits. After the removal of zinc pressure, bacterial phenotype reversed back to the original state; however, certain changes at the genomic, transcriptomic, and proteomic level remained. Why is this important? The extensive and intensive use of supplements in animal feed effects the intestinal microbiota of livestock and this may negatively impact the health of animals and people. Therefore, it is crucial to understand and monitor the impact of feed supplements on intestinal microorganisms in order to adequately assess and prevent potential health risks.

Molecular characterization of Listeria monocytogenes strains isolated from imported food in China from 14 countries/regions, 2003-2018

Listeria monocytogenes (Lm) is associated with severe foodborne infections and ubiquitous in the nature. Identification of characteristics of Lm transmission through trading of food products is essential for rapidly tracking Lm sources and controlling dissemination of listeriosis. In this study, a total of 44 Lm strains were isolated from food products originating from 14 countries/regions during 2003-2018 at the Shanghai port. The genomes of these Lm strains were sequenced by high-throughput sequencing. Multilocus sequence typing (MLST) analysis showed that 43 isolates were divided into 17 sequence types (STs). The distribution of STs was decentralized, with the dominant ST2 accounting for only 18.18% of the strains. The LM63 strain did not match with any of the existing STs. Core-genome MLST (cgMLST) analysis based on 1748 core genes categorized the 44 strains into 30 cgMLST types (CTs), with CT10153 and CT7892 as the most predominant CTs. Notably, LM63 and LM67 shared the same CT in the cgMLST analysis. The phylogenetic analysis based on single-copy homologous genes revealed that the 44 Lm strains were primarily classified into two lineages. The SNP analysis also indicated that these strains were roughly divided into two clades, with strains in the first clade mainly collected earlier than those in the second clade, which were predominantly collected from 2010 onwards. The analysis using the virulence factor database (VFDB) indicated that the virulence gene inlJ was the most prevalent among these 44 strains. Notably, ddrA, msbA, and sugC were enriched in this dataset, requiring further clarification of their roles in Listeria through future studies. These results might provide a clue for understanding of the global epidemiology and surveillance of Lm and present insights for implementing effective measures to reduce or prevent Listeria contamination outbreaks in imported food products.

Comparative pathogenomic analysis reveals a highly tetanus toxin-producing clade of Clostridium tetani isolates in Japan

IMPORTANCE

C. tetani is a spore-forming, anaerobic bacterium that produces a toxin causing muscle stiffness and paralysis. Tetanus is preventable with the toxoid vaccine, but it remains a significant public health threat in regions with low vaccine coverage. However, there are relatively few isolates and limited genomic information available worldwide. In Japan, about 100 cases are reported each year, but there have been no nationwide surveys of isolates, and no genomic information from Japanese isolates has been published. In our study, we analyzed the genomes of 151 strains from a limited survey of soil in Kumamoto, Japan. Our findings revealed a high degree of genetic diversity, and we also identified a subset of strains that produced significantly more toxin, which provides new insights into the pathogenesis of tetanus. Our findings lay the foundation for future studies to investigate the distribution and evolution of C. tetani in Japan and neighboring countries.

De novo identification of expressed cancer somatic mutations from single-cell RNA sequencing data

Identifying expressed somatic mutations from single-cell RNA sequencing data de novo is challenging but highly valuable. We propose RESA - Recurrently Expressed SNV Analysis, a computational framework to identify expressed somatic mutations from scRNA-seq data. RESA achieves an average precision of 0.77 on three in silico spike-in datasets. In extensive benchmarking against existing methods using 19 datasets, RESA consistently outperforms them. Furthermore, we applied RESA to analyze intratumor mutational heterogeneity in a melanoma drug resistance dataset. By enabling high precision detection of expressed somatic mutations, RESA substantially enhances the reliability of mutational analysis in scRNA-seq. RESA is available at https://github.com/ShenLab-Genomics/RESA .

Randomized phase II study of preoperative afatinib in untreated head and neck cancers: predictive and pharmacodynamic biomarkers of activity

There is no strong and reliable predictive biomarker in head and neck squamous cell carcinoma (HNSCC) for EGFR inhibitors. We aimed to identify predictive and pharmacodynamic biomarkers of efficacy of afatinib, a pan-HER tyrosine kinase inhibitor, in a window-of-opportunity trial (NCT01415674). Multi-omics analyses were carried out on pre-treatment biopsy and surgical specimen for biological assessment of afatinib activity. Sixty-one treatment-naïve and operable HNSCC patients were randomised to afatinib 40 mg/day for 21-28 days versus no treatment. Afatinib produced a high rate of metabolic response. Responders had a higher expression of pERK1/2 (P = 0.02) and lower expressions of pHER4 (P = 0.03) and pRB1 (P = 0.002) in pre-treatment biopsy compared to non-responders. At the cellular level, responders displayed an enrichment of tumor-infiltrating B cells under afatinib (P = 0.02). At the molecular level, NF-kappa B signaling was over-represented among upregulated genes in non-responders (P < 0.001; FDR = 0.01). Although exploratory, phosphoproteomics-based biomarkers deserve further investigations as predictors of afatinib efficacy.

Scalable noninvasive amplicon-based precision sequencing (SNAPseq) for genetic diagnosis and screening of β-thalassemia and sickle cell disease using a next-generation sequencing platform

β-hemoglobinopathies such as β-thalassemia (BT) and Sickle cell disease (SCD) are inherited monogenic blood disorders with significant global burden. Hence, early and affordable diagnosis can alleviate morbidity and reduce mortality given the lack of effective cure. Currently, Sanger sequencing is considered to be the gold standard genetic test for BT and SCD, but it has a very low throughput requiring multiple amplicons and more sequencing reactions to cover the entire HBB gene. To address this, we have demonstrated an extraction-free single amplicon-based approach for screening the entire β-globin gene with clinical samples using Scalable noninvasive amplicon-based precision sequencing (SNAPseq) assay catalyzing with next-generation sequencing (NGS). We optimized the assay using noninvasive buccal swab samples and simple finger prick blood for direct amplification with crude lysates. SNAPseq demonstrates high sensitivity and specificity, having a 100% agreement with Sanger sequencing. Furthermore, to facilitate seamless reporting, we have created a much simpler automated pipeline with comprehensive resources for pathogenic mutations in BT and SCD through data integration after systematic classification of variants according to ACMG and AMP guidelines. To the best of our knowledge, this is the first report of the NGS-based high throughput SNAPseq approach for the detection of both BT and SCD in a single assay with high sensitivity in an automated pipeline.

GenoPipe: identifying the genotype of origin within (epi)genomic datasets

Confidence in experimental results is critical for discovery. As the scale of data generation in genomics has grown exponentially, experimental error has likely kept pace despite the best efforts of many laboratories. Technical mistakes can and do occur at nearly every stage of a genomics assay (i.e. cell line contamination, reagent swapping, tube mislabelling, etc.) and are often difficult to identify post-execution. However, the DNA sequenced in genomic experiments contains certain markers (e.g. indels) encoded within and can often be ascertained forensically from experimental datasets. We developed the Genotype validation Pipeline (GenoPipe), a suite of heuristic tools that operate together directly on raw and aligned sequencing data from individual high-throughput sequencing experiments to characterize the underlying genome of the source material. We demonstrate how GenoPipe validates and rescues erroneously annotated experiments by identifying unique markers inherent to an organism's genome (i.e. epitope insertions, gene deletions and SNPs).

Genomic and epigenomic analysis of plasma cell-free DNA identifies stemness features associated with worse survival in AR -altered lethal prostate cancer

Metastatic castration-resistant prostate cancer (mCRPC) resistant to androgen receptor (AR)-targeted agents is often lethal. Unfortunately, biomarkers for this deadly disease remain under investigation, and underpinning mechanisms are ill-understood. Here, we applied deep sequencing to ∼100 mCRPC patients prior to the initiation of first-line AR-targeted therapy, which detected AR /enhancer alterations in over a third of patients, which correlated with lethality. To delve into the mechanism underlying why these patients with cell-free AR /enhancer alterations developed more lethal prostate cancer, we next performed genome-wide cell-free DNA epigenomics. Strikingly, we found that binding sites for transcription factors associated with developmental stemness were nucleosomally more accessible. These results were corroborated using cell-free DNA methylation data, as well as tumor RNA sequencing from a held-out cohort of mCRPC patients. Thus, we validated the importance of AR /enhancer alterations as a prognostic biomarker in lethal mCRPC, and showed that the underlying mechanism for lethality involves reprogramming developmental states toward increased stemness.

Utility of wastewater genomic surveillance compared to clinical surveillance to track the spread of the SARS-CoV-2 Omicron variant across England

The world has moved into a new stage of managing the SARS-CoV-2 pandemic with minimal restrictions and reduced testing in the population, leading to reduced genomic surveillance of virus variants in individuals. Wastewater-based epidemiology (WBE) can provide an alternative means of tracking virus variants in the population but decision-makers require confidence that it can be applied to a national scale and is comparable to individual testing data. We analysed 19,911 samples from 524 wastewater sites across England at least twice a week between November 2021 and February 2022, capturing sewage from >70% of the English population. We used amplicon-based sequencing and the phylogeny based de-mixing tool Freyja to estimate SARS-CoV-2 variant frequencies and compared these to the variant dynamics observed in individual testing data from clinical and community settings. We show that wastewater data can reconstruct the spread of the Omicron variant across England since November 2021 in close detail and aligns closely with epidemiological estimates from individual testing data. We also show the temporal and spatial spread of Omicron within London. Our wastewater data further reliably track the transition between Omicron subvariants BA1 and BA2 in February 2022 at regional and national levels. Our demonstration that WBE can track the fast-paced dynamics of SARS-CoV-2 variant frequencies at a national scale and closely match individual testing data in time shows that WBE can reliably fill the monitoring gap left by reduced individual testing in a more affordable way.

Global hypermethylation of the N6-methyladenosine RNA modification associated with apple heterografting

Grafting can facilitate better scion performance and is widely used in plants. Numerous studies have studied the involvement of mRNAs, small RNAs, and epigenetic regulations in the grafting process. However, it remains unclear whether the mRNA N6-methyladenosine (m6A) modification participates in the apple (Malus x domestica Borkh.) grafting process. Here, we decoded the landscape of m6A modification profiles in 'Golden delicious' (a cultivar, Gd) and Malus prunifolia 'Fupingqiuzi' (a unique rootstock with resistance to environmental stresses, Mp), as well as their heterografted and self-grafted plants. Interestingly, global hypermethylation of m6A occurred in both heterografted scion and rootstock compared with their self-grafting controls. Gene Ontology (GO) term enrichment analysis showed that grafting-induced differentially m6A-modified genes were mainly involved in RNA processing, epigenetic regulation, stress response, and development. Differentially m6A-modified genes harboring expression alterations were mainly involved in various stress responses and fatty acid metabolism. Furthermore, grafting-induced mobile mRNAs with m6A and gene expression alterations mainly participated in ABA synthesis and transport (e.g. carotenoid cleavage dioxygenase 1 [CCD1] and ATP-binding cassette G22 [ABCG22]) and abiotic and biotic stress responses, which might contribute to the better performance of heterografted plants. Additionally, the DNA methylome analysis also demonstrated the DNA methylation alterations during grafting. Downregulated expression of m6A methyltransferase gene MdMTA (ortholog of METTL3) in apples induced the global m6A hypomethylation and distinctly activated the expression level of DNA demethylase gene MdROS1 (REPRESSOR OF SILENCING 1) showing the possible association between m6A and 5mC methylation in apples. Our results reveal the m6A modification profiles in the apple grafting process and enhance our understanding of the m6A regulatory mechanism in plant biological processes.

Genomic profiling of subcutaneous patient-derived xenografts reveals immune constraints on tumor evolution in childhood solid cancer

Subcutaneous patient-derived xenografts (PDXs) are an important tool for childhood cancer research. Here, we describe a resource of 68 early passage PDXs established from 65 pediatric solid tumor patients. Through genomic profiling of paired PDXs and patient tumors (PTs), we observe low mutational similarity in about 30% of the PT/PDX pairs. Clonal analysis in these pairs show an aggressive PT minor subclone seeds the major clone in the PDX. We show evidence that this subclone is more immunogenic and is likely suppressed by immune responses in the PT. These results suggest interplay between intratumoral heterogeneity and antitumor immunity may underlie the genetic disparity between PTs and PDXs. We further show that PDXs generally recapitulate PTs in copy number and transcriptomic profiles. Finally, we report a gene fusion LRPAP1-PDGFRA. In summary, we report a childhood cancer PDX resource and our study highlights the role of immune constraints on tumor evolution.

Mitogenome sequences of domestic cats demonstrate lineage expansions and dynamic mutation processes in a mitochondrial minisatellite

BACKGROUND

As a population genetic tool, mitochondrial DNA is commonly divided into the ~ 1-kb control region (CR), in which single nucleotide variant (SNV) diversity is relatively high, and the coding region, in which selective constraint is greater and diversity lower, but which provides an informative phylogeny. In some species, the CR contains variable tandemly repeated sequences that are understudied due to heteroplasmy. Domestic cats (Felis catus) have a recent origin and therefore traditional CR-based analysis of populations yields only a small number of haplotypes.

RESULTS

To increase resolution we used Nanopore sequencing to analyse 119 cat mitogenomes via a long-amplicon approach. This greatly improves discrimination (from 15 to 87 distinct haplotypes in our dataset) and defines a phylogeny showing similar starlike topologies within all major clades (haplogroups), likely reflecting post-domestication expansion. We sequenced RS2, a CR tandem array of 80-bp repeat units, placing RS2 array structures within the phylogeny and increasing overall haplotype diversity. Repeat number varies between 3 and 12 (median: 4) with over 30 different repeat unit types differing largely by SNVs. Five SNVs show evidence of independent recurrence within the phylogeny, and seven are involved in at least 11 instances of rapid spread along repeat arrays within haplogroups.

CONCLUSIONS

In defining mitogenome variation our study provides key information for the forensic genetic analysis of cat hair evidence, and for the first time a phylogenetically informed picture of tandem repeat variation that reveals remarkably dynamic mutation processes at work in the mitochondrion.

Whole mitochondrial genome analysis in highland Tibetans: further matrilineal genetic structure exploration

Introduction: The Qinghai-Tibet Plateau is one of the last terrestrial environments conquered by modern humans. Tibetans are among the few high-altitude settlers in the world, and understanding the genetic profile of Tibetans plays a pivotal role in studies of anthropology, genetics, and archaeology. Methods: In this study, we investigated the maternal genetic landscape of Tibetans based on the whole mitochondrial genome collected from 145 unrelated native Lhasa Tibetans. Molecular diversity indices, haplotype diversity (HD), Tajima's D and Fu's Fs were calculated and the Bayesian Skyline Plot was obtained to determining the genetic profile and population fluctuation of Lhasa Tibetans. To further explore the genetic structure of Lhasa Tibetans, we collected 107 East Asian reference populations to perform principal component analysis (PCA), multidimensional scaling (MDS), calculated Fst values and constructed phylogenetic tree. Results: The maternal genetic landscape of Tibetans showed obvious East Asian characteristics, M9a (28.28%), R (11.03%), F1 (12.41%), D4 (9.66%), N (6.21%), and M62 (4.14%) were the dominant haplogroups. The results of PCA, MDS, Fst and phylogenetic tree were consistent: Lhasa Tibetans clustered with other highland Tibeto-Burman speakers, there was obvious genetic homogeneity of Tibetans in Xizang, and genetic similarity between Tibetans and northern Han people and geographically adjacent populations was found. In addition, specific maternal lineages of Tibetans also be determined in this study. Discussion: In general, this study further shed light on long-time matrilineal continuity on the Tibetan Plateau and the genetic connection between Tibetans and millet famers in the Yellow River Basin, and further revealed that multiple waves of population interaction and admixture during different historical periods between lowland and highland populations shaped the maternal genetic profile of Tibetans.