A rapid assessment of wastewater for genomic surveillance of SARS-CoV-2 variants at sewershed scale in Louisville, KY.

In this communication, we report on the genomic surveillance of SARS-CoV-2 using wastewater samples in Jefferson County, KY. In February 2021, we analyzed seven wastewater samples for SARS-CoV-2 genomic surveillance. Variants observed in smaller catchment areas, such as neighborhood manhole locations, were not necessarily consistent when compared to associated variant results in downstream treatment plants, suggesting catchment size or population could impact the ability to detect diversity.

Gene alteration in zebrafish exposed to a mixture of substances of abuse.

A recent surge in the use and abuse of diverse prescribed psychotic and illicit drugs necessitates the surveillance of drug residues in source water and the associated ecological impacts of chronic exposure to the aquatic organism. Thirty-six psychotic and illicit drug residues were determined in discharged wastewater from two centralized municipal wastewater treatment facilities and two wastewater receiving creeks for seven consecutive days in Kentucky. Zebrafish (Danio rerio) larvae were exposed to the environmental relevant mixtures of all drug residues, all illicit drugs, and all prescribed psychotic drugs. The extracted RNA from fish homogenates was sequenced, and differentially expressed sequences were analyzed for known or predicted nervous system expression, and screened annotated protein-coding genes to the true environmental cocktail mixture. Illicit stimulant (cocaine and one metabolite), opioids (methadone, methadone metabolite, and oxycodone), hallucinogen (MDA), benzodiazepine (oxazepam and temazepam), carbamazepine, and all target selective serotonin reuptake inhibitors including sertraline, fluoxetine, venlafaxine, and citalopram were quantified in 100% of collected samples from both creeks. The high dose cocktail mixture exposure group revealed the largest group of differentially expressed genes: 100 upregulated and 77 downregulated (p ≤ 0.05; q ≤ 0.05). The top 20 differentially expressed sequences in each exposure group comprise 82 unique transcripts corresponding to 74% annotated genes, 7% non-coding sequences, and 19% uncharacterized sequences. Among 61 differentially expressed sequences that corresponded to annotated protein-coding genes, 23 (38%) genes or their homologs are known to be expressed in the nervous system of fish or other organisms. Several of the differentially expressed sequences are associated primarily with the immune system, including several major histocompatibility complex class I and interferon-induced proteins. Interleukin-1 beta (downregulated in this study) abnormalities are considered a risk factor for psychosis. This is the first study to assess the contributions of multiple classes of psychotic and illicit drugs in combination with developmental gene expression.

Gene structure prediction and alternative splicing analysis using genomically aligned ESTs.

With the availability of a nearly complete sequence of the human genome, aligning expressed sequence tags (EST) to the genomic sequence has become a practical and powerful strategy for gene prediction. Elucidating gene structure is a complex problem requiring the identification of splice junctions, gene boundaries, and alternative splicing variants. We have developed a software tool, Transcript Assembly Program (TAP), to delineate gene structures using genomically aligned EST sequences. TAP assembles the joint gene structure of the entire genomic region from individual splice junction pairs, using a novel algorithm that uses the EST-encoded connectivity and redundancy information to sort out the complex alternative splicing patterns. A method called polyadenylation site scan (PASS) has been developed to detect poly-A sites in the genome. TAP uses these predictions to identify gene boundaries by segmenting the joint gene structure at polyadenylated terminal exons. Reconstructing 1007 known transcripts, TAP scored a sensitivity (Sn) of 60% and a specificity (Sp) of 92% at the exon level. The gene boundary identification process was found to be accurate 78% of the time. also reports alternative splicing patterns in EST alignments. An analysis of alternative splicing in 1124 genic regions suggested that more than half of human genes undergo alternative splicing. Surprisingly, we saw an absolute majority of the detected alternative splicing events affect the coding region. Furthermore, the evolutionary conservation of alternative splicing between human and mouse was analyzed using an EST-based approach. (See http://stl.wustl.edu/~zkan/TAP/)

Sequence assembly validation by multiple restriction digest fragment coverage analysis.

DNA sequence analysis depends on the accurate assembly of fragment reads for the determination of a consensus sequence. This report examines the possibility of analyzing multiple, independent restriction digests as a method for testing the fidelity of sequence assembly. A dynamic programming algorithm to determine the maximum likelihood alignment of error prone electrophoretic mobility data to the expected fragment mobilities given the consensus sequence and restriction enzymes is derived and used to assess the likelihood of detecting rearrangements in genomic sequencing projects. The method is shown to reliably detect errors in sequence fragment assembly without the necessity of making reference to an overlying physical map. An html form-based interface is available at http:/(/)www.ibc.wustl.edu/services/validate. html.