SequencErr: measuring and suppressing sequencer errors in next-generation sequencing data

Degradation of edible mushroom waste by Hermetia illucens L. and consequent adaptation of its gut microbiota

The edible fungus industry is one of the pillar industries in the Yunnan-Guizhou Plateau, China. The expansion of the planting scale has led to the release of various mushroom residues, such as mushroom feet, and other wastes, which are not treated adequately, resulting in environmental pollution. This study investigated the ability of black soldier fly (Hermetia illucens L.) larvae (BSFL) to degrade mushroom waste. Moreover, this study analyzed changes in the intestinal bacterial community and gene expression of BSFL after feeding on mushroom waste. Under identical feeding conditions, the remaining amount of mushroom waste in Pleurotus ostreatus treatment group was reduced by 18.66%, whereas that in Flammulina velutipes treatment group was increased by 31.08%. Regarding gut microbial diversity, compared with wheat bran-treated control group, Dysgonomonas, Providencia, Enterococcus, Pseudochrobactrum, Actinomyces, Morganella, Ochrobactrum, Raoultella, and Ignatzschineria were the most abundant bacteria in the midgut of BSFL in F. velutipes treatment group. Furthermore, Dysgonomonas, Campylobacter, Providencia, Ignatzschineria, Actinomyces, Enterococcus, Morganella, Raoultella, and Pseudochrobactrum were the most abundant bacteria in the midgut of BSFL in P. ostreatus treatment group. Compared with wheat bran-treated control group, 501 upregulated and 285 downregulated genes were identified in F. velutipes treatment group, whereas 211 upregulated and 43 downregulated genes were identified in P. ostreatus treatment group. Using Kyoto Encyclopedia of Genes and Genomes and Gene Ontology enrichment analyses, we identified 14 differentially expressed genes (DEGs) related to amino sugar and nucleotide sugar metabolism in F. velutipes treatment group, followed by 12 DEGs related to protein digestion and absorption. Moreover, in P. ostreatus treatment group, two DEGs were detected for fructose and mannose metabolism, and two were noted for fatty acid metabolism. These results indicate that feeding on edible mushroom waste can alter the intestinal microbial community structure of BSFL; moreover, the larval intestine can generate a corresponding feedback. These changes contribute to the degradation of edible mushroom waste by BSFL and provide a reference for treating edible mushroom waste using BSFL.

Transcriptomic atlas for hypoxia and following re-oxygenation in Ancherythroculter nigrocauda heart and brain tissues: insights into gene expression, alternative splicing, and signaling pathways

Hypoxia is a mounting problem that affects the world's freshwaters, with severe consequence for many species, including death and large economical loss. The hypoxia problem has increased recently due to the combined effects of water eutrophication and global warming. In this study, we investigated the transcriptome atlas for the bony fish Ancherythroculter nigrocauda under hypoxia for 1.5, 3, and 4.5 h and its recovery to normal oxygen levels in heart and brain tissues. We sequenced 21 samples for brain and heart tissues (a total of 42 samples) plus three control samples and obtained an average of 32.40 million raw reads per sample, and 95.24% mapping rate of the filtered clean reads. This robust transcriptome dataset facilitated the discovery of 52,428 new transcripts and 6,609 novel genes. In the heart tissue, the KEGG enrichment analysis showed that genes linked to the Vascular smooth muscle contraction and MAPK and VEGF signaling pathways were notably altered under hypoxia. Re-oxygenation introduced changes in genes associated with abiotic stimulus response and stress regulation. In the heart tissue, weighted gene co-expression network analysis pinpointed a module enriched in insulin receptor pathways that was correlated with hypoxia. Conversely, in the brain tissue, the response to hypoxia was characterized by alterations in the PPAR signaling pathway, and re-oxygenation influenced the mTOR and FoxO signaling pathways. Alternative splicing analysis identified an average of 27,226 and 28,290 events in the heart and brain tissues, respectively, with differential events between control and hypoxia-stressed groups. This study offers a holistic view of transcriptomic adaptations in A. nigrocauda heart and brain tissues under oxygen stress and emphasizes the role of gene expression and alternative splicing in the response mechanisms.

Many purported pseudogenes in bacterial genomes are bona fide genes

BACKGROUND

Microbial genomes are largely comprised of protein coding sequences, yet some genomes contain many pseudogenes caused by frameshifts or internal stop codons. These pseudogenes are believed to result from gene degradation during evolution but could also be technical artifacts of genome sequencing or assembly.

RESULTS

Using a combination of observational and experimental data, we show that many putative pseudogenes are attributable to errors that are incorporated into genomes during assembly. Within 126,564 publicly available genomes, we observed that nearly identical genomes often substantially differed in pseudogene counts. Causal inference implicated assembler, sequencing platform, and coverage as likely causative factors. Reassembly of genomes from raw reads confirmed that each variable affects the number of putative pseudogenes in an assembly. Furthermore, simulated sequencing reads corroborated our observations that the quality and quantity of raw data can significantly impact the number of pseudogenes in an assembler dependent fashion. The number of unexpected pseudogenes due to internal stops was highly correlated (R2 = 0.96) with average nucleotide identity to the ground truth genome, implying relative pseudogene counts can be used as a proxy for overall assembly correctness. Applying our method to assemblies in RefSeq resulted in rejection of 3.6% of assemblies due to significantly elevated pseudogene counts. Reassembly from real reads obtained from high coverage genomes showed considerable variability in spurious pseudogenes beyond that observed with simulated reads, reinforcing the finding that high coverage is necessary to mitigate assembly errors.

CONCLUSIONS

Collectively, these results demonstrate that many pseudogenes in microbial genome assemblies are actually genes. Our results suggest that high read coverage is required for correct assembly and indicate an inflated number of pseudogenes due to internal stops is indicative of poor overall assembly quality.

SJPedPanel: A pan-cancer gene panel for childhood malignancies

Background

Large scale genomics projects have identified driver alterations for most childhood cancers that provide reliable biomarkers for clinical diagnosis and disease monitoring using targeted sequencing. However, there is lack of a comprehensive panel that matches the list of known driver genes. Here we fill this gap by developing SJPedPanel for childhood cancers.

Results

SJPedPanel covers 5,275 coding exons of 357 driver genes, 297 introns frequently involved in rearrangements that generate fusion oncoproteins, commonly amplified/deleted regions (e.g., MYCN for neuroblastoma, CDKN2A and PAX5 for B-/T-ALL, and SMARCB1 for AT/RT), and 7,590 polymorphism sites for interrogating tumors with aneuploidy, such as hyperdiploid and hypodiploid B-ALL or 17q gain neuroblastoma. We used driver alterations reported from an established real-time clinical genomics cohort (n=253) to validate this gene panel. Among the 485 pathogenic variants reported, our panel covered 417 variants (86%). For 90 rearrangements responsible for oncogenic fusions, our panel covered 74 events (82%). We re-sequenced 113 previously characterized clinical specimens at an average depth of 2,500X using SJPedPanel and recovered 354 (91%) of the 389 reported pathogenic variants. We then investigated the power of this panel in detecting mutations from specimens with low tumor purity (as low as 0.1%) using cell line-based dilution experiments and discovered that this gene panel enabled us to detect ∼80% variants with allele fraction of 0.2%, while the detection rate decreases to ∼50% when the allele fraction is 0.1%. We finally demonstrate its utility in disease monitoring on clinical specimens collected from AML patients in morphologic remission.

Conclusions

SJPedPanel enables the detection of clinically relevant genetic alterations including rearrangements responsible for subtype-defining fusions for childhood cancers by targeted sequencing of ∼0.15% of human genome. It will enhance the analysis of specimens with low tumor burdens for cancer monitoring and early detection.

Differential Response of MYB Transcription Factor Gene Transcripts to Circadian Rhythm in Tea Plants (Camellia sinensis)

The circadian clock refers to the formation of a certain rule in the long-term evolution of an organism, which is an invisible 'clock' in the body of an organism. As one of the largest TF families in higher plants, the MYB transcription factor is involved in plant growth and development. MYB is also inextricably correlated with the circadian rhythm. In this study, the transcriptome data of the tea plant 'Baiyeyihao' were measured at a photoperiod interval of 4 h (24 h). A total of 25,306 unigenes were obtained, including 14,615 unigenes that were annotated across 20 functional categories within the GO classification. Additionally, 10,443 single-gene clusters were annotated to 11 sublevels of metabolic pathways using KEGG. Based on the results of gene annotation and differential gene transcript analysis, 22 genes encoding MYB transcription factors were identified. The G10 group in the phylogenetic tree had 13 members, of which 5 were related to the circadian rhythm, accounting for 39%. The G1, G2, G8, G9, G15, G16, G18, G19, G20, G21 and G23 groups had no members associated with the circadian rhythm. Among the 22 differentially expressed MYB transcription factors, 3 members of LHY, RVE1 and RVE8 were core circadian rhythm genes belonging to the G10, G12 and G10 groups, respectively. Real-time fluorescence quantitative PCR was used to detect and validate the expression of the gene transcripts encoding MYB transcription factors associated with the circadian rhythm.

Single duplex DNA sequencing with CODEC detects mutations with high sensitivity

Detecting mutations from single DNA molecules is crucial in many fields but challenging. Next-generation sequencing (NGS) affords tremendous throughput but cannot directly sequence double-stranded DNA molecules ('single duplexes') to discern the true mutations on both strands. Here we present Concatenating Original Duplex for Error Correction (CODEC), which confers single duplex resolution to NGS. CODEC affords 1,000-fold higher accuracy than NGS, using up to 100-fold fewer reads than duplex sequencing. CODEC revealed mutation frequencies of 2.72 × 10-8 in sperm of a 39-year-old individual, and somatic mutations acquired with age in blood cells. CODEC detected genome-wide, clonal hematopoiesis mutations from single DNA molecules, single mutated duplexes from tumor genomes and liquid biopsies, microsatellite instability with 10-fold greater sensitivity and mutational signatures, and specific tumor mutations with up to 100-fold fewer reads. CODEC enables more precise genetic testing and reveals biologically significant mutations, which are commonly obscured by NGS errors.

satmut_utils: a simulation and variant calling package for multiplexed assays of variant effect

The impact of millions of individual genetic variants on molecular phenotypes in coding sequences remains unknown. Multiplexed assays of variant effect (MAVEs) are scalable methods to annotate relevant variants, but existing software lacks standardization, requires cumbersome configuration, and does not scale to large targets. We present satmut_utils as a flexible solution for simulation and variant quantification. We then benchmark MAVE software using simulated and real MAVE data. We finally determine mRNA abundance for thousands of cystathionine beta-synthase variants using two experimental methods. The satmut_utils package enables high-performance analysis of MAVEs and reveals the capability of variants to alter mRNA abundance.

Etiology of oncogenic fusions in 5,190 childhood cancers and its clinical and therapeutic implication

Oncogenic fusions formed through chromosomal rearrangements are hallmarks of childhood cancer that define cancer subtype, predict outcome, persist through treatment, and can be ideal therapeutic targets. However, mechanistic understanding of the etiology of oncogenic fusions remains elusive. Here we report a comprehensive detection of 272 oncogenic fusion gene pairs by using tumor transcriptome sequencing data from 5190 childhood cancer patients. We identify diverse factors, including translation frame, protein domain, splicing, and gene length, that shape the formation of oncogenic fusions. Our mathematical modeling reveals a strong link between differential selection pressure and clinical outcome in CBFB-MYH11. We discover 4 oncogenic fusions, including RUNX1-RUNX1T1, TCF3-PBX1, CBFA2T3-GLIS2, and KMT2A-AFDN, with promoter-hijacking-like features that may offer alternative strategies for therapeutic targeting. We uncover extensive alternative splicing in oncogenic fusions including KMT2A-MLLT3, KMT2A-MLLT10, C11orf95-RELA, NUP98-NSD1, KMT2A-AFDN and ETV6-RUNX1. We discover neo splice sites in 18 oncogenic fusion gene pairs and demonstrate that such splice sites confer therapeutic vulnerability for etiology-based genome editing. Our study reveals general principles on the etiology of oncogenic fusions in childhood cancer and suggests profound clinical implications including etiology-based risk stratification and genome-editing-based therapeutics.

The genetic basis of adaptation to copper pollution in Drosophila melanogaster

Introduction: Heavy metal pollutants can have long lasting negative impacts on ecosystem health and can shape the evolution of species. The persistent and ubiquitous nature of heavy metal pollution provides an opportunity to characterize the genetic mechanisms that contribute to metal resistance in natural populations. Methods: We examined variation in resistance to copper, a common heavy metal contaminant, using wild collections of the model organism Drosophila melanogaster. Flies were collected from multiple sites that varied in copper contamination risk. We characterized phenotypic variation in copper resistance within and among populations using bulked segregant analysis to identify regions of the genome that contribute to copper resistance. Results and Discussion: Copper resistance varied among wild populations with a clear correspondence between resistance level and historical exposure to copper. We identified 288 SNPs distributed across the genome associated with copper resistance. Many SNPs had population-specific effects, but some had consistent effects on copper resistance in all populations. Significant SNPs map to several novel candidate genes involved in refolding disrupted proteins, energy production, and mitochondrial function. We also identified one SNP with consistent effects on copper resistance in all populations near CG11825, a gene involved in copper homeostasis and copper resistance. We compared the genetic signatures of copper resistance in the wild-derived populations to genetic control of copper resistance in the Drosophila Synthetic Population Resource (DSPR) and the Drosophila Genetic Reference Panel (DGRP), two copper-naïve laboratory populations. In addition to CG11825, which was identified as a candidate gene in the wild-derived populations and previously in the DSPR, there was modest overlap of copper-associated SNPs between the wild-derived populations and laboratory populations. Thirty-one SNPs associated with copper resistance in wild-derived populations fell within regions of the genome that were associated with copper resistance in the DSPR in a prior study. Collectively, our results demonstrate that the genetic control of copper resistance is highly polygenic, and that several loci can be clearly linked to genes involved in heavy metal toxicity response. The mixture of parallel and population-specific SNPs points to a complex interplay between genetic background and the selection regime that modifies the effects of genetic variation on copper resistance.

Advances of mRNA vaccine in tumor: a maze of opportunities and challenges

High-frequency mutations in tumor genomes could be exploited as an asset for developing tumor vaccines. In recent years, with the tremendous breakthrough in genomics, intelligence algorithm, and in-depth insight of tumor immunology, it has become possible to rapidly target genomic alterations in tumor cell and rationally select vaccine targets. Among a variety of candidate vaccine platforms, the early application of mRNA was limited by instability low efficiency and excessive immunogenicity until the successful development of mRNA vaccines against SARS-COV-2 broken of technical bottleneck in vaccine preparation, allowing tumor mRNA vaccines to be prepared rapidly in an economical way with good performance of stability and efficiency. In this review, we systematically summarized the classification and characteristics of tumor antigens, the general process and methods for screening neoantigens, the strategies of vaccine preparations and advances in clinical trials, as well as presented the main challenges in the current mRNA tumor vaccine development.

Genomic and Transcriptomic Analysis Reveal Multiple Strategies for the Cadmium Tolerance in Vibrio parahaemolyticus N10-18 Isolated from Aquatic Animal Ostrea gigas Thunberg

The waterborne Vibrio parahaemolyticus can cause acute gastroenteritis, wound infection, and septicemia in humans. Pollution of heavy metals in aquatic environments is proposed to link high incidence of the multidrug-resistant (MDR) pathogen. Nevertheless, the genome evolution and heavy metal tolerance mechanism of V. parahaemolyticus in aquatic animals remain to be largely unveiled. Here, we overcome the limitation by characterizing an MDR V. parahaemolyticus N10-18 isolate with high cadmium (Cd) tolerance using genomic and transcriptomic techniques. The draft genome sequence (4,910,080 bp) of V. parahaemolyticus N10-18 recovered from Ostrea gigas Thunberg was determined, and 722 of 4653 predicted genes had unknown function. Comparative genomic analysis revealed mobile genetic elements (n = 11) and heavy metal and antibiotic-resistance genes (n = 38 and 7). The bacterium significantly changed cell membrane structure to resist the Cd2+ (50 μg/mL) stress (p < 0.05). Comparative transcriptomic analysis revealed seven significantly altered metabolic pathways elicited by the stress. The zinc/Cd/mercury/lead transportation and efflux and the zinc ATP-binding cassette (ABC) transportation were greatly enhanced; metal and iron ABC transportation and thiamine metabolism were also up-regulated; conversely, propanoate metabolism and ribose and maltose ABC transportation were inhibited (p < 0.05). The results of this study demonstrate multiple strategies for the Cd tolerance in V. parahaemolyticus.

Next-Generation Sequencing Applications for the Study of Fungal Pathogens

Next-generation sequencing (NGS) has become a widely used technology in biological research. NGS applications for clinical pathogen detection have become vital technologies. It is increasingly common to perform fast, accurate, and specific detection of clinical specimens using NGS. Pathogenic fungi with high virulence and drug resistance cause life-threatening clinical infections. NGS has had a significant biotechnological impact on detecting bacteria and viruses but is not equally applicable to fungi. There is a particularly urgent clinical need to use NGS to help identify fungi causing infections and prevent negative impacts. This review summarizes current research on NGS applications for fungi and offers a visual method of fungal detection. With the development of NGS and solutions for overcoming sequencing limitations, we suggest clinicians test specimens as soon as possible when encountering infections of unknown cause, suspected infections in vital organs, or rapidly progressive disease.

Therapeutic and prognostic insights from the analysis of cancer mutational signatures

The somatic mutations in each cancer genome are caused by multiple mutational processes, each of which leaves a characteristic imprint (or 'signature'), potentially caused by specific etiologies or exposures. Deconvolution of these signatures offers a glimpse into the evolutionary history of individual tumors. Recent work has shown that mutational signatures may also yield therapeutic and prognostic insights, including the identification of cell-intrinsic signatures as biomarkers of drug response and prognosis. For example, mutational signatures indicating homologous recombination deficiency are associated with poly(ADP)-ribose polymerase (PARP) inhibitor sensitivity, whereas APOBEC-associated signatures are associated with ataxia telangiectasia and Rad3-related kinase (ATR) inhibitor sensitivity. Furthermore, therapy-induced mutational signatures implicated in cancer progression have also been uncovered, including the identification of thiopurine-induced TP53 mutations in leukemia. In this review, we explore the various ways mutational signatures can reveal new therapeutic and prognostic insights, thus extending their traditional role in identifying disease etiology.

Advancing NGS quality control to enable measurement of actionable mutations in circulating tumor DNA

The primary objective of the FDA-led Sequencing and Quality Control Phase 2 (SEQC2) project is to develop standard analysis protocols and quality control metrics for use in DNA testing to enhance scientific research and precision medicine. This study reports a targeted next-generation sequencing (NGS) method that will enable more accurate detection of actionable mutations in circulating tumor DNA (ctDNA) clinical specimens. To accomplish this, a synthetic internal standard spike-in was designed for each actionable mutation target, suitable for use in NGS following hybrid capture enrichment and unique molecular index (UMI) or non-UMI library preparation. When mixed with contrived ctDNA reference samples, internal standards enabled calculation of technical error rate, limit of blank, and limit of detection for each variant at each nucleotide position in each sample. True-positive mutations with variant allele fraction too low for detection by current practice were detected with this method, thereby increasing sensitivity.

CRISPR Screens in Synthetic Lethality and Combinatorial Therapies for Cancer

Cancer is a complex disease resulting from the accumulation of genetic dysfunctions. Tumor heterogeneity causes the molecular variety that divergently controls responses to chemotherapy, leading to the recurrent problem of cancer reappearance. For many decades, efforts have focused on identifying essential tumoral genes and cancer driver mutations. More recently, prompted by the clinical success of the synthetic lethality (SL)-based therapy of the PARP inhibitors in homologous recombinant deficient tumors, scientists have centered their novel research on SL interactions (SLI). The state of the art to find new genetic interactions are currently large-scale forward genetic CRISPR screens. CRISPR technology has rapidly evolved to be a common tool in the vast majority of laboratories, as tools to implement CRISPR screen protocols are available to all researchers. Taking advantage of SLI, combinatorial therapies have become the ultimate model to treat cancer with lower toxicity, and therefore better efficiency. This review explores the CRISPR screen methodology, integrates the up-to-date published findings on CRISPR screens in the cancer field and proposes future directions to uncover cancer regulation and individual responses to chemotherapy.