Assessing different next-generation sequencing technologies for wastewater-based epidemiology.

Wastewater-based epidemiology has proven to be an important public health asset during the COVID-19 pandemic. It can provide less biassed and more cost-effective population-level monitoring of the disease burden as compared to clinical testing. An essential component of SARS-CoV-2 wastewater monitoring is next-generation sequencing, providing genomic data to identify and quantify circulating viral strains rapidly. However, the specific choice of sequencing method influences the quality and timeliness of generated data and hence its usefulness for wastewater-based pathogen surveillance. Here, we systematically benchmarked Illumina Novaseq 6000, Element Aviti, ONT R9.4.1 MinION flow cell, and ONT R9.4.1 Flongle flow cell sequencing data to facilitate the selection of sequencing technology. Using a time series of wastewater samples from influent of six wastewater treatment plants throughout Switzerland, along with spike-in experiments, we show that higher sequencing error rates of ONT Nanopore sequencing reduce the accuracy of estimates of the relative abundance of viral variants, but the overall trend is in good concordance among all technologies. We find that the sequencing runtime for ONT Nanopore flow cells can be reduced to as little as five hours without significant impact on the quality of variant estimates. Our findings suggest that SARS-CoV-2 variant tracking is readily achievable with all tested technologies, albeit with different tradeoffs in terms of cost, timeliness and accuracy.

Tlr4 Deletion Modulates Cytokine and Extracellular Matrix Expression in Chronic Spinal Cord Injury, Leading to Improved Secondary Damage and Functional Recovery.

Toll-like receptors (TLRs) play an important role in the innate immune response after CNS injury. Although TLR4 is one of the best characterized, its role in chronic stages after spinal cord injury (SCI) is not well understood. We examined the role of TLR4 signaling in injury-induced responses at 1 d, 7 d, and 8 weeks after spinal cord contusion injury in adult female TLR4 null and wild-type mice. Analyses include secondary damage, a range of transcriptome and protein analyses of inflammatory, cell death, and extracellular matrix (ECM) molecules, as well as immune cell infiltration and changes in axonal sprouting and locomotor recovery. Lack of TLR4 signaling results in reduced neuronal and myelin loss, reduced activation of NFκB, and decreased expression of inflammatory cytokines and necroptotic cell death pathway at a late time point (8 weeks) after injury. TLR4 null mice also showed reduction of scar-related ECM molecules at 8 weeks after SCI, accompanied by increase in ECM molecules associated with perineuronal nets, increased sprouting of serotonergic fibers, and improved locomotor recovery. These findings reveal novel effects of TLR4 signaling in chronic SCI. We show that TLR4 influences inflammation, cell death, and ECM deposition at late-stage post-injury when secondary injury processes are normally considered to be over. This highlights the potential for late-stage targeting of TLR4 as a potential therapy for chronic SCI.

Early detection and surveillance of SARS-CoV-2 genomic variants in wastewater using COJAC.

The continuing emergence of SARS-CoV-2 variants of concern and variants of interest emphasizes the need for early detection and epidemiological surveillance of novel variants. We used genomic sequencing of 122 wastewater samples from three locations in Switzerland to monitor the local spread of B.1.1.7 (Alpha), B.1.351 (Beta) and P.1 (Gamma) variants of SARS-CoV-2 at a population level. We devised a bioinformatics method named COJAC (Co-Occurrence adJusted Analysis and Calling) that uses read pairs carrying multiple variant-specific signature mutations as a robust indicator of low-frequency variants. Application of COJAC revealed that a local outbreak of the Alpha variant in two Swiss cities was observable in wastewater up to 13 d before being first reported in clinical samples. We further confirmed the ability of COJAC to detect emerging variants early for the Delta variant by analysing an additional 1,339 wastewater samples. While sequencing data of single wastewater samples provide limited precision for the quantification of relative prevalence of a variant, we show that replicate and close-meshed longitudinal sequencing allow for robust estimation not only of the local prevalence but also of the transmission fitness advantage of any variant. We conclude that genomic sequencing and our computational analysis can provide population-level estimates of prevalence and fitness of emerging variants from wastewater samples earlier and on the basis of substantially fewer samples than from clinical samples. Our framework is being routinely used in large national projects in Switzerland and the UK.

Quantification of the spread of SARS-CoV-2 variant B.1.1.7 in Switzerland.

BACKGROUND: In December 2020, the United Kingdom (UK) reported a SARS-CoV-2 Variant of Concern (VoC) which is now named B.1.1.7. Based on initial data from the UK and later data from other countries, this variant was estimated to have a transmission fitness advantage of around 40-80 % (Volz et al., 2021; Leung et al., 2021; Davies et al., 2021). AIM: This study aims to estimate the transmission fitness advantage and the effective reproductive number of B.1.1.7 through time based on data from Switzerland. METHODS: We generated whole genome sequences from 11.8 % of all confirmed SARS-CoV-2 cases in Switzerland between 14 December 2020 and 11 March 2021. Based on these data, we determine the daily frequency of the B.1.1.7 variant and quantify the variant's transmission fitness advantage on a national and a regional scale. RESULTS: We estimate B.1.1.7 had a transmission fitness advantage of 43-52 % compared to the other variants circulating in Switzerland during the study period. Further, we estimate B.1.1.7 had a reproductive number above 1 from 01 January 2021 until the end of the study period, compared to below 1 for the other variants. Specifically, we estimate the reproductive number for B.1.1.7 was 1.24 [1.07-1.41] from 01 January until 17 January 2021 and 1.18 [1.06-1.30] from 18 January until 01 March 2021 based on the whole genome sequencing data. From 10 March to 16 March 2021, once B.1.1.7 was dominant, we estimate the reproductive number was 1.14 [1.00-1.26] based on all confirmed cases. For reference, Switzerland applied more non-pharmaceutical interventions to combat SARS-CoV-2 on 18 January 2021 and lifted some measures again on 01 March 2021. CONCLUSION: The observed increase in B.1.1.7 frequency in Switzerland during the study period is as expected based on observations in the UK. In absolute numbers, B.1.1.7 increased exponentially with an estimated doubling time of around 2-3.5 weeks. To monitor the ongoing spread of B.1.1.7, our plots are available online.

Mapping Spatial Genetic Landscapes in Tissue Sections through Microscale Integration of Sampling Methodology into Genomic Workflows.

In cancer research, genomic profiles are often extracted from homogenized macrodissections of tissues, with the histological context lost and a large fraction of material underutilized. Pertinently, the spatial genomic landscape provides critical complementary information in deciphering disease heterogeneity and progression. Microscale sampling methods such as microdissection to obtain such information are often destructive to a sizeable fraction of the biopsy sample, thus showing limited multiplexability and adaptability to different assays. A modular microfluidic technology is here implemented to recover cells at the microscale from tumor tissue sections, with minimal disruption of unsampled areas and tailored to interface with genome profiling workflows, which is directed here toward evaluating intratumoral genomic heterogeneity. The integrated workflow-GeneScape-is used to evaluate heterogeneity in a metastatic mammary carcinoma, showing distinct single nucleotide variants and copy number variations in different tumor tissue regions, suggesting the polyclonal origin of the metastasis as well as development driven by multiple location-specific drivers.

Comparative Analysis of Single-Cell RNA Sequencing Platforms and Methods.

Single-cell RNA sequencing (scRNA-seq) offers great new opportunities for increasing our understanding of complex biological processes. In particular, development of an accurate Human Cell Atlas is largely dependent on the rapidly advancing technologies and molecular chemistries employed in scRNA-seq. These advances have already allowed an increase in throughput for scRNA-seq from 96 to 80,000 cells on a single instrument run by capturing cells within nanoliter droplets. Although this increase in throughput is critical for many experimental questions, a thorough comparison between microfluidic-based, plate-based, and droplet-based technologies or between multiple available platforms utilizing these technologies is largely lacking. Here, we report scRNA-seq data from SUM149PT cells treated with the histone deacetylase inhibitor trichostatin A versus untreated controls across several scRNA-seq platforms (Fluidigm C1, WaferGen iCell8, 10x Genomics Chromium Controller, and Illumina/BioRad ddSEQ). The primary goal of this project was to demonstrate RNA sequencing methods for profiling the ultra-low amounts of RNA present in individual cells, and this report discusses the results of the study, as well as technical challenges and lessons learned and present general guidelines for best practices in sample preparation and analysis.

De Novo Genome Assembly of the Japanese Wheat Cultivar Norin 61 Highlights Functional Variation in Flowering Time and Fusarium-Resistant Genes in East Asian Genotypes.

Bread wheat is a major crop that has long been the focus of basic and breeding research. Assembly of its genome has been difficult because of its large size and allohexaploid nature (AABBDD genome). Following the first reported assembly of the genome of the experimental strain Chinese Spring (CS), the 10+ Wheat Genomes Project was launched to produce multiple assemblies of worldwide modern cultivars. The only Asian cultivar in the project is Norin 61, a representative Japanese cultivar adapted to grow across a broad latitudinal range, mostly characterized by a wet climate and a short growing season. Here, we characterize the key aspects of its chromosome-scale genome assembly spanning 15 Gb with a raw scaffold N50 of 22 Mb. Analysis of the repetitive elements identified chromosomal regions unique to Norin 61 that encompass a tandem array of the pathogenesis-related 13 family. We report novel copy-number variations in the B homeolog of the florigen gene FT1/VRN3, pseudogenization of its D homeolog and the association of its A homeologous alleles with the spring/winter growth habit. Furthermore, the Norin 61 genome carries typical East Asian functional variants different from CS, ranging from a single nucleotide to multi-Mb scale. Examples of such variation are the Fhb1 locus, which confers Fusarium head-blight resistance, Ppd-D1a, which confers early flowering, Glu-D1f for Asian noodle quality and Rht-D1b, which introduced semi-dwarfism during the green revolution. The adoption of Norin 61 as a reference assembly for functional and evolutionary studies will enable comprehensive characterization of the underexploited Asian bread wheat diversity.

Multiple wheat genomes reveal global variation in modern breeding.

Advances in genomics have expedited the improvement of several agriculturally important crops but similar efforts in wheat (Triticum spp.) have been more challenging. This is largely owing to the size and complexity of the wheat genome1, and the lack of genome-assembly data for multiple wheat lines2,3. Here we generated ten chromosome pseudomolecule and five scaffold assemblies of hexaploid wheat to explore the genomic diversity among wheat lines from global breeding programs. Comparative analysis revealed extensive structural rearrangements, introgressions from wild relatives and differences in gene content resulting from complex breeding histories aimed at improving adaptation to diverse environments, grain yield and quality, and resistance to stresses4,5. We provide examples outlining the utility of these genomes, including a detailed multi-genome-derived nucleotide-binding leucine-rich repeat protein repertoire involved in disease resistance and the characterization of Sm16, a gene associated with insect resistance. These genome assemblies will provide a basis for functional gene discovery and breeding to deliver the next generation of modern wheat cultivars.

Contribution of epigenetic variation to adaptation in Arabidopsis.

In plants, transgenerational inheritance of some epialleles has been demonstrated but it remains controversial whether epigenetic variation is subject to selection and contributes to adaptation. Simulating selection in a rapidly changing environment, we compare phenotypic traits and epigenetic variation between Arabidopsis thaliana populations grown for five generations under selection and their genetically nearly identical ancestors. Selected populations of two distinct genotypes show significant differences in flowering time and plant architecture, which are maintained for at least 2-3 generations in the absence of selection. While we cannot detect consistent genetic changes, we observe a reduction of epigenetic diversity and changes in the methylation state of about 50,000 cytosines, some of which are associated with phenotypic changes. Thus, we propose that epigenetic variation is subject to selection and can contribute to rapid adaptive responses, although the extent to which epigenetics plays a role in adaptation is still unclear.

Correction to: Limited utility of tissue micro-arrays in detecting intra-tumoral heterogeneity in stem cell characteristics and tumor progression markers in breast cancer.

Following publication of the original article [1], a typesetting mistake is reported. For Fig. 7b, a copy of Fig. 6b has been published. The correct Fig. 7b is given in this correction and the original article has been updated.

Limited utility of tissue micro-arrays in detecting intra-tumoral heterogeneity in stem cell characteristics and tumor progression markers in breast cancer.

BACKGROUND: Intra-tumoral heterogeneity has been recently addressed in different types of cancer, including breast cancer. A concept describing the origin of intra-tumoral heterogeneity is the cancer stem-cell hypothesis, proposing the existence of cancer stem cells that can self-renew limitlessly and therefore lead to tumor progression. Clonal evolution in accumulated single cell genomic alterations is a further possible explanation in carcinogenesis. In this study, we addressed the question whether intra-tumoral heterogeneity can be reliably detected in tissue-micro-arrays in breast cancer by comparing expression levels of conventional predictive/prognostic tumor markers, tumor progression markers and stem cell markers between central and peripheral tumor areas. METHODS: We analyzed immunohistochemical expression and/or gene amplification status of conventional prognostic tumor markers (ER, PR, HER2, CK5/6), tumor progression markers (PTEN, PIK3CA, p53, Ki-67) and stem cell markers (mTOR, SOX2, SOX9, SOX10, SLUG, CD44, CD24, TWIST) in 372 tissue-micro-array samples from 72 breast cancer patients. Expression levels were compared between central and peripheral tumor tissue areas and were correlated to histopathological grading. 15 selected cases additionally underwent RNA sequencing for transcriptome analysis. RESULTS: No significant difference in any of the analyzed between central and peripheral tumor areas was seen with any of the analyzed methods/or results that showed difference. Except mTOR, PIK3CA and SOX9 (nuclear) protein expression, all markers correlated significantly (p < 0.05) with histopathological grading both in central and peripheral areas. CONCLUSION: Our results suggest that intra-tumoral heterogeneity of stem-cell and tumor-progression markers cannot be reliably addressed in tissue-micro-array samples in breast cancer. However, most markers correlated strongly with histopathological grading confirming prognostic information as expression profiles were independent on the site of the biopsy was taken.

Drop drying on surfaces determines chemical reactivity - the specific case of immobilization of oligonucleotides on microarrays.

BACKGROUND: Drop drying is a key factor in a wide range of technical applications, including spotted microarrays. The applied nL liquid volume provides specific reaction conditions for the immobilization of probe molecules to a chemically modified surface. RESULTS: We investigated the influence of nL and µL liquid drop volumes on the process of probe immobilization and compare the results obtained to the situation in liquid solution. In our data, we observe a strong relationship between drop drying effects on immobilization and surface chemistry. In this work, we present results on the immobilization of dye labeled 20mer oligonucleotides with and without an activating 5'-aminoheptyl linker onto a 2D epoxysilane and a 3D NHS activated hydrogel surface. CONCLUSIONS: Our experiments identified two basic processes determining immobilization. First, the rate of drop drying that depends on the drop volume and the ambient relative humidity. Oligonucleotides in a dried spot react unspecifically with the surface and long reaction times are needed. 3D hydrogel surfaces allow for immobilization in a liquid environment under diffusive conditions. Here, oligonucleotide immobilization is much faster and a specific reaction with the reactive linker group is observed. Second, the effect of increasing probe concentration as a result of drop drying. On a 3D hydrogel, the increasing concentration of probe molecules in nL spotting volumes accelerates immobilization dramatically. In case of µL volumes, immobilization depends on whether the drop is allowed to dry completely. At non-drying conditions, very limited immobilization is observed due to the low oligonucleotide concentration used in microarray spotting solutions. The results of our study provide a general guideline for microarray assay development. They allow for the initial definition and further optimization of reaction conditions for the immobilization of oligonucleotides and other probe molecule classes to different surfaces in dependence of the applied spotting and reaction volume.

AGRONOMICS1: a new resource for Arabidopsis transcriptome profiling.

Transcriptome profiling has become a routine tool in biology. For Arabidopsis (Arabidopsis thaliana), the Affymetrix ATH1 expression array is most commonly used, but it lacks about one-third of all annotated genes present in the reference strain. An alternative are tiling arrays, but previous designs have not allowed the simultaneous analysis of both strands on a single array. We introduce AGRONOMICS1, a new Affymetrix Arabidopsis microarray that contains the complete paths of both genome strands, with on average one 25mer probe per 35-bp genome sequence window. In addition, the new AGRONOMICS1 array contains all perfect match probes from the original ATH1 array, allowing for seamless integration of the very large existing ATH1 knowledge base. The AGRONOMICS1 array can be used for diverse functional genomics applications such as reliable expression profiling of more than 30,000 genes, detection of alternative splicing, and chromatin immunoprecipitation coupled to microarrays (ChIP-chip). Here, we describe the design of the array and compare its performance with that of the ATH1 array. We find results from both microarrays to be of similar quality, but AGRONOMICS1 arrays yield robust expression information for many more genes, as expected. Analysis of the ATH1 probes on AGRONOMICS1 arrays produces results that closely mirror those of ATH1 arrays. Finally, the AGRONOMICS1 array is shown to be useful for ChIP-chip experiments. We show that heterochromatic H3K9me2 is strongly confined to the gene body of target genes in euchromatic chromosome regions, suggesting that spreading of heterochromatin is limited outside of pericentromeric regions.

16S rRNA gene-based phylogenetic microarray for simultaneous identification of members of the genus Burkholderia.

For cultivation-independent and highly parallel analysis of members of the genus Burkholderia, an oligonucleotide microarray (phylochip) consisting of 131 hierarchically nested 16S rRNA gene-targeted oligonucleotide probes was developed. A novel primer pair was designed for selective amplification of a 1.3 kb 16S rRNA gene fragment of Burkholderia species prior to microarray analysis. The diagnostic performance of the microarray for identification and differentiation of Burkholderia species was tested with 44 reference strains of the genera Burkholderia, Pandoraea, Ralstonia and Limnobacter. Hybridization patterns based on presence/absence of probe signals were interpreted semi-automatically using the novel likelihood-based strategy of the web-tool Phylo- Detect. Eighty-eight per cent of the reference strains were correctly identified at the species level. The evaluated microarray was applied to investigate shifts in the Burkholderia community structure in acidic forest soil upon addition of cadmium, a condition that selected for Burkholderia species. The microarray results were in agreement with those obtained from phylogenetic analysis of Burkholderia 16S rRNA gene sequences recovered from the same cadmiumcontaminated soil, demonstrating the value of the Burkholderia phylochip for determinative and environmental studies.

Processing protocols for high quality glass-based microarrays: applications in DNA, peptide, antibody, and carbohydrate microarraying.

Based on the selection of a suitable surface chemistry and bearing the option for optimization using a defined workflow, standard experimental protocols for the processing of microarrays can be used as starting points for a successful experiment. In Chapters 2 and 3, general quality considerations and the selection of surface chemistry have been discussed. A workflow for the selection of slide surface architectures and the optimization of microarray processing parameters also has been described. In the present article, processing parameters for DNA, peptide, antibody, and carbohydrate microarrays are outlined that serve as a first recommended step in the iterative establishment process. For a number of popular applications of microarray technology the outlined protocols can be applied to directly generate high-quality results.

Quality considerations and selection of surface chemistry for glass-based DNA, peptide, antibody, carbohydrate, and small molecule microarrays.

The complexity of workflows for the production of high quality microarrays asks for the careful evaluation and implementation of materials and methods. As a cornerstone of the whole microarray process, the microarray substrate has to be chosen appropriately and a number of crucial considerations in respect to matching the research question with the technical requirements and possibilities have to be taken into account. In the following, how to lay the fundamental for high performance microarray experiments by evaluating basic quality requirements and the selection of suitable slide surface architectures for a variety of applications was concentrated.

Optimization workflow for the processing of high quality glass-based microarrays: applications in DNA, peptide, antibody, and carbohydrate microarraying.

As the performance of microarray experiments is directly dependent on the quality of the materials, the suitability of the protocols, and the accuracy of the work performed, optimization of existing microarray workflows is needed in almost every experiment to achieve higher quality and meaningfulness of the generated data. In the following, we describe a workflow for the optimization of microarray processing parameters, based on the previously selected surface structure. Four simple model experiments with dye-labeled compounds is used to determine crucial experimental parameters including spotting concentration, spotting solution, immobilization efficiency, and blocking conditions even in cases where recommendations from the slide manufacturer or from the literature are missing. In this article, processing parameters for DNA, peptide, antibody, and carbohydrate microarrays are outlined. The applicability of the model experiments is demonstrated and described in detail on the example of short oligonucleotides.