Wastewater-based epidemiology applied at the building-level reveals distinct virome profiles based on the age of the contributing individuals.

BACKGROUND: Human viruses released into the environment can be detected and characterized in wastewater. The study of wastewater virome offers a consolidated perspective on the circulation of viruses within a population. Because the occurrence and severity of viral infections can vary across a person's lifetime, studying the virome in wastewater samples contributed by various demographic segments can provide valuable insights into the prevalence of viral infections within these segments. In our study, targeted enrichment sequencing was employed to characterize the human virome in wastewater at a building-level scale. This was accomplished through passive sampling of wastewater in schools, university settings, and nursing homes in two cities in Catalonia. Additionally, sewage from a large urban wastewater treatment plant was analysed to serve as a reference for examining the collective excreted human virome. RESULTS: The virome obtained from influent wastewater treatment plant samples showcased the combined viral presence from individuals of varying ages, with astroviruses and human bocaviruses being the most prevalent, followed by human adenoviruses, polyomaviruses, and papillomaviruses. Significant variations in the viral profiles were observed among the different types of buildings studied. Mamastrovirus 1 was predominant in school samples, salivirus and human polyomaviruses JC and BK in the university settings while nursing homes showed a more balanced distribution of viral families presenting papillomavirus and picornaviruses and, interestingly, some viruses linked to immunosuppression. CONCLUSIONS: This study shows the utility of building-level wastewater-based epidemiology as an effective tool for monitoring the presence of viruses circulating within specific age groups. It provides valuable insights for public health monitoring and epidemiological studies.

Mitochondrial genome structure and composition in 70 fishes: a key resource for fisheries management in the South Atlantic.

BACKGROUND: Phylogenetic gaps of public databases of reference sequences are a major obstacle for comparative genomics and management of marine resources, particularly in the Global South, where economically important fisheries and conservation flagship species often lack closely-related references. We applied target-enrichment to obtain complete mitochondrial genomes of marine ichthyofauna from the Brazilian coast selected based on economic significance, conservation status and lack of phylogenetically-close references. These included sardines (Dorosomatidae, Alosidae), mackerels (Scombridae) croakers (Sciaenidae), groupers (Epinephelidae) and snappers (Lutjanidae). RESULTS: Custom baits were designed to enrich mitochondrial DNA across a broad phylogenetic range of fishes. Sequencing generated approximately 100k reads per sample, which were assembled in a total of 70 complete mitochondrial genomes and include fifty-two new additions to GenBank, including five species with no previous mitochondrial data. Departures from the typical gene content and order occurred in only three taxa and mostly involved tRNA gene duplications. Start-codons for all genes, except Cytochrome C Oxidase subunit I (COI), were consistently ATG, whilst a wide range of stop-codons deviated from the prevailing TAA. Phylogenetic analysis confirmed assembly accuracy and revealed signs of cryptic diversification within the Mullus genus. Lineage delimitation methods using Sardinella aurita and S. brasiliensis mitochondrial genomes support a single Operational Taxonomic Unit. CONCLUSIONS: Target enrichment was highly efficient, providing complete novel mitochondrial genomes with little sequencing effort. These sequences are deposited in public databases to enable subsequent studies in population genetics and adaptation of Latin American fish species and serve as a vital resource for conservation and management programs that rely on molecular data for species and genus-level identification.

The phylogeny and global biogeography of Primulaceae based on high-throughput DNA sequence data.

The angiosperm family Primulaceae is morphologically diverse and distributed nearly worldwide. However, phylogenetic uncertainty has obstructed the identification of major morphological and biogeographic transitions within the clade. We used target capture sequencing with the Angiosperms353 probes, taxon-sampling encompassing nearly all genera of the family, tree-based sequence curation, and multiple phylogenetic approaches to investigate the major clades of Primulaceae and their relationship to other Ericales. We generated dated phylogenetic trees and conducted broad-scale biogeographic analyses as well as stochastic character mapping of growth habit. We show that Ardisia, a pantropical genus and the largest in the family, is not monophyletic, with at least 19 smaller genera nested within it. Neotropical members of Ardisia and several smaller genera form a clade, an ancestor of which arrived in the Neotropics and began diversifying about 20 Ma. This Neotropical clade is most closely related to Elingamita and Tapeinosperma, which are most diverse on islands of the Pacific. Both Androsace and Primula are non-monophyletic by the inclusion of smaller genera. Ancestral state reconstructions revealed that there have either been parallel transitions to an herbaceous habit in Primuloideae, Samolus, and at least three lineages of Myrsinoideae, or a common ancestor of nearly all Primulaceae was herbaceous. Our results provide a robust estimate of phylogenetic relationships across Primulaceae and show that a revised classification of Myrsinoideae and several other clades within the family is necessary to render all genera monophyletic.

A Large-Scale Candidate-Gene Association Mapping for Drought Tolerance and Agronomic Traits in Sugarcane.

Dissection of the genetic loci controlling drought tolerance traits with a complex genetic inheritance is important for drought-tolerant sugarcane improvement. In this study, we conducted a large-scale candidate gene association study of 649 candidate genes in a sugarcane diversity panel to identify genetic variants underlying agronomic traits and drought tolerance indices evaluated in plant cane and ratoon cane under water-stressed (WS) and non-stressed (NS) environments. We identified 197 significant marker-trait associations (MTAs) in 141 candidate genes associated with 18 evaluated traits with the Bonferroni correction threshold (α = 0.05). Out of the total, 95 MTAs in 78 candidate genes and 62 MTAs in 58 candidate genes were detected under NS and WS conditions, respectively. Most MTAs were found only in specific water regimes and crop seasons. These MTAs explained 7.93-30.52% of phenotypic variation. Association mapping results revealed that 34, 59, and 104 MTAs involved physiological and molecular adaptation, phytohormone metabolism, and drought-inducible genes. They identified 19 pleiotropic genes associated with more than one trait and many genes related to drought tolerance indices. The genetic and genomic resources identified in this study will enable the combining of yield-related traits and sugar-related traits with agronomic value to optimize the yield of sugarcane cultivars grown under drought-stressed and non-stressed environments.

Broadening the horizon of crop research: a decade of advancements in plant molecular genetics to divulge phenotype governing genes.

MAIN CONCLUSION: Advancements in sequencing, genotyping, and computational technologies during the last decade (2011-2020) enabled new forward-genetic approaches, which subdue the impediments of precise gene mapping in varied crops. The modern crop improvement programs rely heavily on two major steps-trait-associated QTL/gene/marker's identification and molecular breeding. Thus, it is vital for basic and translational crop research to identify genomic regions that govern the phenotype of interest. Until the advent of next-generation sequencing, the forward-genetic techniques were laborious and time-consuming. Over the last 10 years, advancements in the area of genome assembly, genotyping, large-scale data analysis, and statistical algorithms have led faster identification of genomic variations regulating the complex agronomic traits and pathogen resistance. In this review, we describe the latest developments in genome sequencing and genotyping along with a comprehensive evaluation of the last 10-year headways in forward-genetic techniques that have shifted the focus of plant research from model plants to diverse crops. We have classified the available molecular genetic methods under bulk-segregant analysis-based (QTL-seq, GradedPool-Seq, QTG-Seq, Exome QTL-seq, and RapMap), target sequence enrichment-based (RenSeq, AgRenSeq, and TACCA), and mutation-based groups (MutMap, NIKS algorithm, MutRenSeq, MutChromSeq), alongside improvements in classical mapping and genome-wide association analyses. Newer methods for outcrossing, heterozygous, and polyploid plant genetics have also been discussed. The use of k-mers has enriched the nature of genetic variants which can be utilized to identify the phenotype-causing genes, independent of reference genomes. We envisage that the recent methods discussed herein will expand the repertoire of useful alleles and help in developing high-yielding and climate-resilient crops.

A catalogue of resistance gene homologs and a chromosome-scale reference sequence support resistance gene mapping in winter wheat.

A resistance gene atlas is an integral component of the breeder's arsenal in the fight against evolving pathogens. Thanks to high-throughput sequencing, catalogues of resistance genes can be assembled even in crop species with large and polyploid genomes. Here, we report on capture sequencing and assembly of resistance gene homologs in a diversity panel of 907 winter wheat genotypes comprising ex situ genebank accessions and current elite cultivars. In addition, we use accurate long-read sequencing and chromosome conformation capture sequencing to construct a chromosome-scale genome sequence assembly of cv. Attraktion, an elite variety representative of European winter wheat. We illustrate the value of our resource for breeders and geneticists by (i) comparing the resistance gene complements in plant genetic resources and elite varieties and (ii) conducting genome-wide associations scans (GWAS) for the fungal diseases yellow rust and leaf rust using reference-based and reference-free GWAS approaches. The gene content under GWAS peaks was scrutinized in the assembly of cv. Attraktion.

Admixture may be extensive among hyperdominant Amazon rainforest tree species.

Admixture is a mechanism by which species of long-lived plants may acquire novel alleles. However, the potential role of admixture in the origin and maintenance of tropical plant diversity is unclear. We ask whether admixture occurs in an ecologically important clade of Eschweilera (Parvifolia clade, Lecythidaceae), which includes some of the most widespread and abundant tree species in Amazonian forests. Using target capture sequencing, we conducted a detailed phylogenomic investigation of 33 species in the Parvifolia clade and investigated specific hypotheses of admixture within a robust phylogenetic framework. We found strong evidence of admixture among three ecologically dominant species, E. coriacea, E. wachenheimii and E. parviflora, but a lack of evidence for admixture among other lineages. Accepted species were largely distinguishable from one another, as was geographic structure within species. We show that hybridization may play a role in the evolution of the most widespread and ecologically variable Amazonian tree species. While admixture occurs among some species of Eschweilera, it has not led to widespread erosion of most species' genetic or morphological identities. Therefore, current morphological based species circumscriptions appear to provide a useful characterization of the clade's lineage diversity.

Incomplete reproductive isolation and low genetic differentiation despite floral divergence across varying geographic scales in Castilleja.

PREMISE: Divergence depends on the strength of selection and frequency of gene flow between taxa, while reproductive isolation relies on mating barriers and geographic distance. Less is known about how these processes interact at early stages of speciation. Here, we compared population-level differentiation in floral phenotype and genetic sequence variation among recently diverged Castilleja to explore patterns of diversification under different scenarios of reproductive isolation. METHODS: Using target enrichment enabled by the Angiosperms353 probe set, we assessed genetic distance among 50 populations of four Castilleja species. We investigated whether patterns of genetic divergence are explained by floral trait variation or geographic distance in two focal groups: the widespread C. sessiliflora and the more restricted C. purpurea species complex. RESULTS: We document that C. sessiliflora and the C. purpurea complex are characterized by high diversity in floral color across varying geographic scales. Despite phenotypic divergence, groups were not well supported in phylogenetic analyses, and little genetic differentiation was found across targeted Angiosperms353 loci. Nonetheless, a principal coordinate analysis of single nucleotide polymorphisms revealed differentiation within C. sessiliflora across floral morphs and geography and less differentiation among species of the C. purpurea complex. CONCLUSIONS: Patterns of genetic distance in C. sessiliflora suggest species cohesion maintained over long distances despite variation in floral traits. In the C. purpurea complex, divergence in floral color across narrow geographic clines may be driven by recent selection on floral color. These contrasting patterns of floral and genetic differentiation reveal that divergence can arise via multiple eco-evolutionary paths.

RLP/K enrichment sequencing; a novel method to identify receptor-like protein (RLP) and receptor-like kinase (RLK) genes.

The identification of immune receptors in crop plants is time-consuming but important for disease control. Previously, resistance gene enrichment sequencing (RenSeq) was developed to accelerate mapping of nucleotide-binding domain and leucine-rich repeat containing (NLR) genes. However, resistances mediated by pattern recognition receptors (PRRs) remain less utilized. Here, our pipeline shows accelerated mapping of PRRs. Effectoromics leads to precise identification of plants with target PRRs, and subsequent RLP/K enrichment sequencing (RLP/KSeq) leads to detection of informative single nucleotide polymorphisms that are linked to the trait. Using Phytophthora infestans as a model, we identified Solanum microdontum plants that recognize the apoplastic effectors INF1 or SCR74. RLP/KSeq in a segregating Solanum population confirmed the localization of the INF1 receptor on chromosome 12, and led to the rapid mapping of the response to SCR74 to chromosome 9. By using markers obtained from RLP/KSeq in conjunction with additional markers, we fine-mapped the SCR74 receptor to a 43-kbp G-LecRK locus. Our findings show that RLP/KSeq enables rapid mapping of PRRs and is especially beneficial for crop plants with large and complex genomes. This work will enable the elucidation and characterization of the nonNLR plant immune receptors and ultimately facilitate informed resistance breeding.

Target enrichment sequencing of 307 germplasm accessions identified ancestry of ancient and modern hybrids and signatures of adaptation and selection in sugarcane (Saccharum spp.), a 'sweet' crop with 'bitter' genomes.

Sugarcane (Saccharum spp.) is a highly energy-efficient crop primarily for sugar and bio-ethanol production. Sugarcane genetics and cultivar improvement have been extremely challenging largely due to its complex genomes with high polyploidy levels. In this study, we deeply sequenced the coding regions of 307 sugarcane germplasm accessions. Nearly five million sequence variations were catalogued. The average of 98× sequence depth enabled different allele dosages of sequence variation to be differentiated in this polyploid collection. With selected high-quality genome-wide SNPs, we performed population genomic studies and environmental association analysis. Results illustrated that the ancient sugarcane hybrids, S. barberi and S. sinense, and modern sugarcane hybrids are significantly different in terms of genomic compositions, hybridization processes and their potential ancestry contributors. Linkage disequilibrium (LD) analysis showed a large extent of LD in sugarcane, with 962.4 Kbp, 2739.2 Kbp and 3573.6 Kbp for S. spontaneum, S. officinarum and modern S. hybrids respectively. Candidate selective sweep regions and genes were identified during domestication and historical selection processes of sugarcane in addition to genes associated with environmental variables at the original locations of the collection. This research provided an extensive amount of genomic resources for sugarcane community and the in-depth population genomic analyses shed light on the breeding and evolution history of sugarcane, a highly polyploid species.

Genome-Wide Association Studies Identified Resistance Loci to Orange Rust and Yellow Leaf Virus Diseases in Sugarcane (Saccharum spp.).

Sugarcane (Saccharum spp.) supplies globally ∼80% of table sugar and 60% of bioethanol. Sugarcane orange rust and Sugarcane yellow leaf virus (SCYLV) are major sugarcane diseases, causing up to 50 and 40% yield losses, respectively. Sugarcane cultivars resistant to these diseases are needed to sustain sugarcane production in several regions. Dissecting DNA sequence variants controlling disease resistance provides a valuable tool for fulfilling a breeding strategy to develop resistant cultivars. In this study, we evaluated disease reactions to orange rust and SCYLV of a sugarcane diversity panel in repeated trials. We conducted a genome-wide association study between high-density markers and disease resistance reactions. We identified 91 putative DNA markers and 82 candidate genes significantly associated with resistance to one of the two diseases. These provide an important genetic resource for finding genes and molecular markers for disease resistance. Our results emphasized the importance of utilizing a wide germplasm collection for breeding resistant sugarcane cultivars.

Target enrichment sequencing in cultivated peanut (Arachis hypogaea L.) using probes designed from transcript sequences.

Enabled by the next generation sequencing, target enrichment sequencing (TES) is a powerful method to enrich genomic regions of interest and to identify sequence variations. The objective of this study was to explore the feasibility of probe design from transcript sequences for TES application in calling sequence variants in peanut, an important allotetraploid crop with a large genome size. In this study, we applied an in-solution hybridization method to enrich DNA sequences of seven peanut genotypes. Our results showed that it is feasible to apply TES with probes designed from transcript sequences in polyploid peanut. Using a set of 31,123 probes, a total of 5131 and 7521 genes were targeted in peanut A and B genomes, respectively. For each genotype used in this study, the probe target capture regions were efficiently covered with high depth. The average on-target rate of sequencing reads was 42.47%, with a significant amount of off-target reads coming from genomic regions homologous to target regions. In this study, when given predefined genomic regions of interest and the same amount of sequencing data, TES provided the highest coverage of target regions when compared to whole genome sequencing, RNA sequencing, and genotyping by sequencing. Single nucleotide polymorphism (SNP) calling and subsequent validation revealed a high validation rate (85.71%) of homozygous SNPs, providing valuable markers for peanut genotyping. This study demonstrated the success of applying TES for SNP identification in peanut, which shall provide valuable suggestions for TES application in other non-model species without a genome reference available.

Local interspecies introgression is the main cause of extreme levels of intraspecific differentiation in mussels.

Structured populations, and replicated zones of contact between species, are an ideal opportunity to study regions of the genome with unusual levels of differentiation; and these can illuminate the genomic architecture of species isolation, and the spread of adaptive alleles across species ranges. Here, we investigated the effects of gene flow on divergence and adaptation in the Mytilus complex of species, including replicated parental populations in quite distant geographical locations. We used target enrichment sequencing of 1269 contigs of a few kb each, including some genes of known function, to infer gene genealogies at a small chromosomal scale. We show that geography is an important determinant of the genomewide patterns of introgression in Mytilus and that gene flow between different species, with contiguous ranges, explained up to half of the intraspecific outliers. This suggests that local introgression is both widespread and tends to affect larger chromosomal regions than purely intraspecific processes. We argue that this situation might be common, and this implies that genome scans should always consider the possibility of introgression from sister species, unsampled differentiated backgrounds, or even extinct relatives, for example Neanderthals in humans. The hypothesis that reticulate evolution over long periods of time contributes widely to adaptation, and to the spatial and genomic reorganization of genetic backgrounds, needs to be more widely considered to make better sense of genome scans.

Phylogenomic reconstruction of Cryptosporidium spp. captured directly from clinical samples reveals extensive genetic diversity.

Cryptosporidium is a leading cause of severe diarrhea and mortality in young children and infants in Africa and southern Asia. More than twenty Cryptosporidium species infect humans, of which C. parvum and C. hominis are the major agents causing moderate to severe diarrhea. Relatively few genetic markers are typically applied to genotype and/or diagnose Cryptosporidium. Most infections produce limited oocysts making it difficult to perform whole genome sequencing (WGS) directly from stool samples. Hence, there is an immediate need to apply WGS strategies to 1) develop high-resolution genetic markers to genotype these parasites more precisely, 2) to investigate endemic regions and detect the prevalence of different genotypes, and the role of mixed infections in generating genetic diversity, and 3) to investigate zoonotic transmission and evolution. To understand Cryptosporidium global population genetic structure, we applied Capture Enrichment Sequencing (CES-Seq) using 74,973 RNA-based 120 nucleotide baits that cover ~92% of the genome of C. parvum. CES-Seq is sensitive and successfully sequenced Cryptosporidium genomic DNA diluted up to 0.005% in human stool DNA. It also resolved mixed strain infections and captured new species of Cryptosporidium directly from clinical/field samples to promote genome-wide phylogenomic analyses and prospective GWAS studies.

Exploring waterborne viruses in groundwater: Quantification and Virome characterization via passive sampling and targeted enrichment sequencing.

Aquifers, which provide drinking water for nearly half the world's population, face significant challenges from microbial contamination, particularly from waterborne viruses such as human adenovirus (HAdV), norovirus (NoV) and enterovirus (EV). This study, conducted as part of the UPWATER project, investigates the sources of urban groundwater contamination using viral passive sampling (VPS) and target enrichment sequencing (TES). We assessed the abundance of eight viral pathogens (HAdV, EV, NoV genogroup I and II, rotavirus, influenza A virus, hepatitis E virus and SARS-CoV-2) and investigated the virome diversity of groundwater in the aquifer of the Besòs River Delta in Catalonia. Over a period of 7 months, we collected 114 samples from the aquifer using nylon and nitrocellulose membranes to adsorb viruses over a 10-day period. Human faecal contamination was detected in nearly 50 % of the groundwater samples, with mean HAdV total counts ranging from 1.23E+02 to 3.66E+03 GC, and occasional detections of EV and NoV GI and GII. In addition, deep sequencing revealed a diverse virome in the aquifer, with detection of human pathogens, including adenovirus, astrovirus, calicivirus, enterovirus, herpesvirus, papillomavirus and rotavirus. Time-integrated sampling using VPS increases the likelihood of virus detection and, when combined with TES, can provide a deeper understanding of virus prevalence in this important water compartment. This approach is expected to streamline long-term monitoring efforts and enable small communities or water managers with limited resources to effectively manage their groundwater reservoirs.

Assessing environmental exposure to viruses in wastewater treatment plant and swine farm scenarios with next-generation sequencing and occupational risk approaches.

Occupational exposure to pathogens can pose health risks. This study investigates the viral exposure of workers in a wastewater treatment plant (WWTP) and a swine farm by analyzing aerosol and surfaces samples. Viral contamination was evaluated using quantitative polymerase chain reaction (qPCR) assays, and target enrichment sequencing (TES) was performed to identify the vertebrate viruses to which workers might be exposed. Additionally, Quantitative Microbial Risk Assessment (QMRA) was conducted to estimate the occupational risk associated with viral exposure for WWTP workers, choosing Human Adenovirus (HAdV) as the reference pathogen. In the swine farm, QMRA was performed as an extrapolation, considering a hypothetical zoonotic virus with characteristics similar to Porcine Adenovirus (PAdV). The modelled exposure routes included aerosol inhalation and oral ingestion through contaminated surfaces and hand-to-mouth contact. HAdV and PAdV were widespread viruses in the WWTP and the swine farm, respectively, by qPCR assays. TES identified human and other vertebrate viruses WWTP samples, including viruses from families such as Adenoviridae, Circoviridae, Orthoherpesviridae, Papillomaviridae, and Parvoviridae. In the swine farm, most of the identified vertebrate viruses were porcine viruses belonging to Adenoviridae, Astroviridae, Circoviridae, Herpesviridae, Papillomaviridae, Parvoviridae, Picornaviridae, and Retroviridae. QMRA analysis revealed noteworthy risks of viral infections for WWTP workers if safety measures are not taken. The probability of illness due to HAdV inhalation was higher in summer compared to winter, while the greatest risk from oral ingestion was observed in workspaces during winter. Swine farm QMRA simulation suggested a potential occupational risk in the case of exposure to a hypothetical zoonotic virus. This study provides valuable insights into WWTP and swine farm worker's occupational exposure to human and other vertebrate viruses. QMRA and NGS analyses conducted in this study will assist managers in making evidence-based decisions, facilitating the implementation of protection measures, and risk mitigation practices for workers.

The longitudinal kinetics of AAV5 vector integration profiles and evaluation of clonal expansion in mice.

Adeno-associated virus (AAV)-based vectors are used clinically for gene transfer and persist as extrachromosomal episomes. A small fraction of vector genomes integrate into the host genome, but the theoretical risk of tumorigenesis depends on vector regulatory features. A mouse model was used to investigate integration profiles of an AAV serotype 5 (AAV5) vector produced using Sf and HEK293 cells that mimic key features of valoctocogene roxaparvovec (AAV5-hFVIII-SQ), a gene therapy for severe hemophilia A. The majority (95%) of vector genome reads were derived from episomes, and mean (± standard deviation) integration frequency was 2.70 ± 1.26 and 1.79 ± 0.86 integrations per 1,000 cells for Sf- and HEK293-produced vector. Longitudinal integration analysis suggested integrations occur primarily within 1 week, at low frequency, and their abundance was stable over time. Integration profiles were polyclonal and randomly distributed. No major differences in integration profiles were observed for either vector production platform, and no integrations were associated with clonal expansion. Integrations were enriched near transcription start sites of genes highly expressed in the liver (p = 1 × 10-4) and less enriched for genes of lower expression. We found no evidence of tumorigenesis or fibrosis caused by the vector integrations.

Comparing molecular and computational approaches for detecting viral integration of AAV gene therapy constructs.

Many current gene therapy targets use recombinant adeno-associated virus (AAV). The majority of delivered AAV therapeutics persist as episomes, separate from host DNA, yet some viral DNA can integrate into host DNA in different proportions and at genomic locations. The potential for viral integration leading to oncogenic transformation has led regulatory agencies to require investigation into AAV integration events following gene therapy in preclinical species. In the present study, tissues were collected from cynomolgus monkeys and mice 6 and 8 weeks, respectively, following administration of an AAV vector delivering transgene cargo. We compared three different next-generation sequencing approaches (shearing extension primer tag selection ligation-mediated PCR, targeted enrichment sequencing [TES], and whole-genome sequencing) to contrast the specificity, scope, and frequency of integration detected by each method. All three methods detected dose-dependent insertions with a limited number of hotspots and expanded clones. While the functional outcome was similar for all three methods, TES was the most cost-effective and comprehensive method of detecting viral integration. Our findings aim to inform the direction of molecular efforts to ensure a thorough hazard assessment of AAV viral integration in our preclinical gene therapy studies.

LINE1-Mediated Reverse Transcription and Genomic Integration of SARS-CoV-2 mRNA Detected in Virus-Infected but Not in Viral mRNA-Transfected Cells.

SARS-CoV-2 sequences can be reverse-transcribed and integrated into the genomes of virus-infected cells by a LINE1-mediated retrotransposition mechanism. Whole-genome sequencing (WGS) methods detected retrotransposed SARS-CoV-2 subgenomic sequences in virus-infected cells overexpressing LINE1, while an enrichment method (TagMap) identified retrotranspositions in cells that did not overexpress LINE1. LINE1 overexpression increased retrotranspositions about 1000-fold as compared to non-overexpressing cells. Nanopore WGS can directly recover retrotransposed viral and flanking host sequences, but its sensitivity depends on the depth of sequencing (a typical 20-fold sequencing depth would only examine 10 diploid cell equivalents). In contrast, TagMap enriches the host-virus junctions and can interrogate up to 20,000 cells and is able to detect rare viral retrotranspositions in LINE1 non-overexpressing cells. Although Nanopore WGS is 10-20-fold more sensitive per tested cell, TagMap can interrogate 1000-2000-fold more cells and, therefore, can identify infrequent retrotranspositions. When comparing SARS-CoV-2 infection and viral nucleocapsid mRNA transfection by TagMap, retrotransposed SARS-CoV-2 sequences were only detected in infected but not in transfected cells. Retrotransposition in virus-infected cells, in contrast to transfected cells, may be facilitated because virus infection, in contrast to viral RNA transfection, results in significantly higher viral RNA levels and stimulates LINE1 expression by causing cellular stress.

DNA methylation biomarkers distinguishing early-stage prostate cancer from benign prostatic hyperplasia.

Background: DNA methylation markers are considered robust diagnostic features in various cancer types, as epigenetic marks are commonly altered during cancer progression. Differentiation between benign prostatic hyperplasia (BPH) and early-stage prostate cancer (PCa) is clinically difficult, relying on the information of the patient's symptoms or levels of prostate-specific antigen. Methods: A total of 42 PCa patients and 11 BPH patients were recruited. Genomic DNA was purified from tissues and used for the library preparation of the target-enriched methylome with enzymatic conversion and a Twist 85 Mbp EM-seq panel. Paired-end sequencing (150 bp) was performed using NovaSeq 6000 or NextSeq 550. After quality control, including adapter trimming and de-duplication of raw sequencing data, differential methylation patterns were analyzed between the BPH and PCa groups. Results: We report DNA methylation patterns existing between BPH and PCa. The major finding is that broad hypermethylation occurred at genic loci in PCa tissues as compared to the BPH. Gene ontology analysis suggested that hypermethylation of genic loci involved in chromatin and transcriptional regulation is involved in cancer progression. We also compared PCa tissues with high Gleason scores to tissues with low Gleason scores. The high-Gleason PCa tissues showed hundreds of focal differentially methylated CpG sites corresponding to genes functioning in cancer cell proliferation or metastasis. This suggests that dissecting early-to-advanced-grade cancer stages requires an in-depth analysis of differential methylation at the single CpG site level. Conclusions: Our study reports that enzymatic methylome sequencing data can be used to distinguish PCa from BPH and advanced PCa from early-stage PCa. The stage-specific methylation patterns in this study will be valuable resources for diagnostic purposes as well as further development of liquid biopsy approaches for the early detection of PCa.