Tracking SARS-CoV-2 variants of concern in wastewater: an assessment of nine computational tools using simulated genomic data.

Wastewater-based surveillance (WBS) is an important epidemiological and public health tool for tracking pathogens across the scale of a building, neighbourhood, city, or region. WBS gained widespread adoption globally during the SARS-CoV-2 pandemic for estimating community infection levels by qPCR. Sequencing pathogen genes or genomes from wastewater adds information about pathogen genetic diversity, which can be used to identify viral lineages (including variants of concern) that are circulating in a local population. Capturing the genetic diversity by WBS sequencing is not trivial, as wastewater samples often contain a diverse mixture of viral lineages with real mutations and sequencing errors, which must be deconvoluted computationally from short sequencing reads. In this study we assess nine different computational tools that have recently been developed to address this challenge. We simulated 100 wastewater sequence samples consisting of SARS-CoV-2 BA.1, BA.2, and Delta lineages, in various mixtures, as well as a Delta-Omicron recombinant and a synthetic 'novel' lineage. Most tools performed well in identifying the true lineages present and estimating their relative abundances and were generally robust to variation in sequencing depth and read length. While many tools identified lineages present down to 1 % frequency, results were more reliable above a 5 % threshold. The presence of an unknown synthetic lineage, which represents an unclassified SARS-CoV-2 lineage, increases the error in relative abundance estimates of other lineages, but the magnitude of this effect was small for most tools. The tools also varied in how they labelled novel synthetic lineages and recombinants. While our simulated dataset represents just one of many possible use cases for these methods, we hope it helps users understand potential sources of error or bias in wastewater sequencing analysis and to appreciate the commonalities and differences across methods.

Characterization of genetic predisposition to molecular subtypes of breast cancer in Brazilian patients.

Introduction: BRCA1 and BRCA2 germline pathogenic variants (GPVs) account for most of the 5-10% of breast cancer (BC) that is attributable to inherited genetic variants. BRCA1 GPVs are associated with the triple negative subtype, whereas BRCA2 GPVs are likely to result in higher grade, estrogen-receptor positive BCs. The contribution of other genes of high and moderate risk for BC has not been well defined and risk estimates to specific BC subtypes is lacking, especially for an admixed population like Brazilian. Objective: The aim of this study is to evaluate the value of a multigene panel in detecting germline mutations in cancer-predisposing genes for Brazilian BC patients and its relation with molecular subtypes and the predominant molecular ancestry. Patients and methods: A total of 321 unrelated BC patients who fulfilled NCCN criteria for BRCA1/2 testing between 2016-2018 were investigated with a 94-genes panel. Molecular subtypes were retrieved from medical records and ancestry-specific variants were obtained from off-target reads obtained from the sequencing data. Results: We detected 83 GPVs in 81 patients (positivity rate of 25.2%). Among GPVs, 47% (39/83) were identified in high-risk BC genes (BRCA1/2, PALB2 and TP53) and 18% (15/83) in moderate-penetrance genes (ATM, CHEK2 and RAD51C). The remainder of the GPVs (35% - 29/83), were identified in lower-risk genes. As for the molecular subtypes, triple negative BC had a mutation frequency of 31.6% (25/79), with predominance in BRCA1 (12.6%; 10/79). Among the luminal subtypes, except Luminal B HER2-positive, 18.7% (29/155) had GPV with BRCA1/2 genes contributing 7.1% (11/155) and non-BRCA1/2 genes, 12.9% (20/155). For Luminal B HER2-positive subtype, 40% (16/40) had GPVs, with a predominance of ATM gene (15% - 6/40) and BRCA2 with only 2.5% (1/40). Finally, HER2-enriched subtype presented a mutation rate of 30.8% (4/13) with contribution of BRCA2 of 7.5% (1/13) and non-BRCA1/2 of 23% (3/13). Variants of uncertain significance (VUS) were identified in 77.6% (249/321) of the patients and the number of VUS was increased in patients with Asian and Native American ancestry. Conclusion: The multigene panel contributed to identify GPVs in genes other than BRCA1/2, increasing the positivity of the genetic test from 9.6% (BRCA1/2) to 25.2% and, considering only the most clinically relevant BC predisposing genes, to 16.2%. These results indicate that women with clinical criteria for hereditary BC may benefit from a multigene panel testing, as it allows identifying GPVs in genes that directly impact the clinical management of these patients and family members.

A mixture model for determining SARS-Cov-2 variant composition in pooled samples.

MOTIVATION: Despite of the fast development of highly effective vaccines to control the current COVID-19 pandemics, the unequal distribution and availability of these vaccines worldwide and the number of people infected in the world lead to the continuous emergence of Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2) variants of concern. Therefore, it is likely that real-time genomic surveillance will be continuously needed as an unceasing monitoring tool, necessary to follow the spread of the disease and the evolution of the virus. In this context, new genomic variants of SARS-CoV-2, including variants refractory to current vaccines, makes genomic surveillance programs tools of utmost importance. Nevertheless, the lack of appropriate analytical tools to quickly and effectively access the viral composition in meta-transcriptomic sequencing data, including environmental surveillance, represent possible challenges that may impact the fast adoption of this approach to mitigate the spread and transmission of viruses. RESULTS: We propose a statistical model for the estimation of the relative frequencies of SARS-CoV-2 variants in pooled samples. This model is built by considering a previously defined selection of genomic polymorphisms that characterize SARS-CoV-2 variants. The methods described here support both raw sequencing reads for polymorphisms-based markers calling and predefined markers in the variant call format. Results obtained using simulated data show that our method is quite effective in recovering the correct variant proportions. Further, results obtained by considering longitudinal data from wastewater samples of two locations in Switzerland agree well with those describing the epidemiological evolution of COVID-19 variants in clinical samples of these locations. Our results show that the described method can be a valuable tool for tracking the proportions of SARS-CoV-2 variants in complex mixtures such as waste water and environmental samples. AVAILABILITY AND IMPLEMENTATION: http://github.com/rvalieris/LCS. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The impact of HER2 overexpression on the miRNA and circRNA transcriptomes in two breast cell lines and their vesicles.

HER2 upregulation is related with poor outcome in many tumor types. Whereas anti-HER2 treatment is the standard approach as adjuvant therapy in HER2-overexpressing breast cancer, the frequent relapses reinforce the need for alternative treatments. Here we used next-generation sequencing (NGS) to evaluate miRNAs and circRNAs in the cell-lines HB4a and C5.2, where the latter is a HER2-overexpressing clone of the former, and also from two different populations of their secreted extracellular vesicles. Whereas circRNA-levels were stable, we found at least 16 miRNAs apparently modulated by HER2-expression. The miR223-3p, miR-421 and miR-21-5p were validated in an independent cohort of 431 breast cancer patients from The Cancer Genome Atlas (TCGA). The consistent modulation of these molecules and their possible involvement in the HER2-axis makes them promising new targets to overcome HER2-activation.

Utility of trio-based exome sequencing in the elucidation of the genetic basis of isolated syndromic intellectual disability: illustrative cases.

INTRODUCTION: Exome sequencing is recognized as a powerful tool for identifying the genetic cause of intellectual disability (ID). It is uncertain, however, whether only the exome of the proband should be sequenced or if the sequencing of parental genomes is also required, and the resulting increase in diagnostic yield justifies the increase in costs. PATIENTS AND METHODS: We sequenced the exomes of eight individuals with sporadic syndromic ID and their parents. RESULTS AND DISCUSSION: Likely pathogenic variants were detected in eight candidate genes, namely homozygous or compound heterozygous variants in three autosomal genes (ADAMTSL2, NALCN, VPS13B), one in an X-linked gene (MID1), and de novo heterozygous variants in four autosomal genes (RYR2, GABBR2, CDK13, DDX3X). Two patients harbored rare variants in two or more candidate genes, while in three other patients no candidate was identified. In five probands (62%), the detected variants explained their clinical findings. The causative recessive variants would have led to diagnosis even without parental exome sequencing, but for the heterozygous dominant ones, the exome trio-based approach was fundamental in the identification of the de novo likely pathogenic variants.

Distinct chromatin functional states correlate with HIV latency reactivation in infected primary CD4+ T cells.

Human immunodeficiency virus (HIV) infection is currently incurable, due to the persistence of latently infected cells. The 'shock and kill' approach to a cure proposes to eliminate this reservoir via transcriptional activation of latent proviruses, enabling direct or indirect killing of infected cells. Currently available latency-reversing agents (LRAs) have however proven ineffective. To understand why, we used a novel HIV reporter strain in primary CD4+ T cells and determined which latently infected cells are reactivatable by current candidate LRAs. Remarkably, none of these agents reactivated more than 5% of cells carrying a latent provirus. Sequencing analysis of reactivatable vs. non-reactivatable populations revealed that the integration sites were distinguishable in terms of chromatin functional states. Our findings challenge the feasibility of 'shock and kill', and suggest the need to explore other strategies to control the latent HIV reservoir.

The G protein-coupled receptor GPR31 promotes membrane association of KRAS.

The product of the KRAS oncogene, KRAS4B, promotes tumor growth when associated with the plasma membrane (PM). PM association is mediated, in part, by farnesylation of KRAS4B, but trafficking of nascent KRAS4B to the PM is incompletely understood. We performed a genome-wide screen to identify genes required for KRAS4B membrane association and identified a G protein-coupled receptor, GPR31. GPR31 associated with KRAS4B on cellular membranes in a farnesylation-dependent fashion, and retention of GPR31 on the endoplasmic reticulum inhibited delivery of KRAS4B to the PM. Silencing of GPR31 expression partially mislocalized KRAS4B, slowed the growth of KRAS-dependent tumor cells, and blocked KRAS-stimulated macropinocytosis. Our data suggest that GPR31 acts as a secretory pathway chaperone for KRAS4B.

Human dendritic cells (DCs) are derived from distinct circulating precursors that are precommitted to become CD1c+ or CD141+ DCs.

In humans, conventional dendritic cells (cDCs) exist as two unique populations characterized by expression of CD1c and CD141. cDCs arise from increasingly restricted but well-defined bone marrow progenitors that include the common DC progenitor that differentiates into the pre-cDC, which is the direct precursor of cDCs. In this study, we show that pre-cDCs in humans are heterogeneous, consisting of two distinct populations of precursors that are precommitted to become either CD1c+ or CD141+ cDCs. The two groups of lineage-primed precursors can be distinguished based on differential expression of CD172a. Both subpopulations of pre-cDCs arise in the adult bone marrow and can be found in cord blood and adult peripheral blood. Gene expression analysis revealed that CD172a+ and CD172a- pre-cDCs represent developmentally discrete populations that differentially express lineage-restricted transcription factors. A clinical trial of Flt3L injection revealed that this cytokine increases the number of both CD172a- and CD172a+ pre-cDCs in human peripheral blood.