Proto-oncogene c-Myb potentiates cisplatin resistance of ovarian cancer cells by downregulating lncRNA NKILA and modulating cancer stemness and LIN28A-let7 axis.

Ovarian cancer is a major gynecological cancer that has poor prognosis associated mainly to its late diagnosis. Cisplatin is an FDA approved ovarian cancer therapy and even though the therapy is initially promising, the patients mostly progress to resistance against cisplatin. The underlying mechanisms are complex and not very clearly understood. Using two different paired cell lines representing cisplatin-sensitive and the cisplatin-resistant ovarian cancer cells, the ES2 and the A2780 parental and cisplatin-resistant cells, we show an elevated proto-oncogene c-Myb in resistant cells. We further show down-regulated lncRNA NKILA in resistant cells with its de-repression in resistant cells when c-Myb is silenced. NKILA negatively correlates with cancer cell and invasion but has no effect on cellular proliferation or cell cycle. C-Myb activates NF-κB signaling which is inhibited by NKILA. The cisplatin resistant cells are also marked by upregulated stem cell markers, particularly LIN28A and OCT4, and downregulated LIN28A-targeted let-7 family miRNAs. Whereas LIN28A and downregulated let-7s individually de-repress c-Myb-mediated cisplatin resistance, the ectopic expression of let-7s attenuates LIN28A effects, thus underlying a c-Myb-NKILA-LIN28A-let-7 axis in cisplatin resistance of ovarian cancer cells that needs to be further explored for therapeutic intervention.

Unveiling Promising Neuroimaging Biomarkers for Schizophrenia Through Clinical and Genetic Perspectives.

Schizophrenia is a complex and serious brain disorder. Neuroscientists have become increasingly interested in using magnetic resonance-based brain imaging-derived phenotypes (IDPs) to investigate the etiology of psychiatric disorders. IDPs capture valuable clinical advantages and hold biological significance in identifying brain abnormalities. In this review, we aim to discuss current and prospective approaches to identify potential biomarkers for schizophrenia using clinical multimodal neuroimaging and imaging genetics. We first described IDPs through their phenotypic classification and neuroimaging genomics. Secondly, we discussed the applications of multimodal neuroimaging by clinical evidence in observational studies and randomized controlled trials. Thirdly, considering the genetic evidence of IDPs, we discussed how can utilize neuroimaging data as an intermediate phenotype to make association inferences by polygenic risk scores and Mendelian randomization. Finally, we discussed machine learning as an optimum approach for validating biomarkers. Together, future research efforts focused on neuroimaging biomarkers aim to enhance our understanding of schizophrenia.

The Causal Relationships Between Gut Microbiota, Brain Volume, and Intelligence: A Two-Step Mendelian Randomization Analysis.

BACKGROUND: Growing evidence indicates that dynamic changes in gut microbiome can affect intelligence; however, whether these relationships are causal remains elusive. We aimed to disentangle the poorly understood causal relationship between gut microbiota and intelligence. METHODS: We performed a 2-sample Mendelian randomization (MR) analysis using genetic variants from the largest available genome-wide association studies of gut microbiota (N = 18,340) and intelligence (N = 269,867). The inverse-variance weighted method was used to conduct the MR analyses complemented by a range of sensitivity analyses to validate the robustness of the results. Considering the close relationship between brain volume and intelligence, we applied 2-step MR to evaluate whether the identified effect was mediated by regulating brain volume (N = 47,316). RESULTS: We found a risk effect of the genus Oxalobacter on intelligence (odds ratio = 0.968 change in intelligence per standard deviation increase in taxa; 95% CI, 0.952-0.985; p = 1.88 × 10-4) and a protective effect of the genus Fusicatenibacter on intelligence (odds ratio = 1.053; 95% CI, 1.024-1.082; p = 3.03 × 10-4). The 2-step MR analysis further showed that the effect of genus Fusicatenibacter on intelligence was partially mediated by regulating brain volume, with a mediated proportion of 33.6% (95% CI, 6.8%-60.4%; p = .014). CONCLUSIONS: Our results provide causal evidence indicating the role of the microbiome in intelligence. Our findings may help reshape our understanding of the microbiota-gut-brain axis and development of novel intervention approaches for preventing cognitive impairment.

A landscape of gene expression regulation for synovium in arthritis.

The synovium is an important component of any synovial joint and is the major target tissue of inflammatory arthritis. However, the multi-omics landscape of synovium required for functional inference is absent from large-scale resources. Here we integrate genomics with transcriptomics and chromatin accessibility features of human synovium in up to 245 arthritic patients, to characterize the landscape of genetic regulation on gene expression and the regulatory mechanisms mediating arthritic diseases predisposition. We identify 4765 independent primary and 616 secondary cis-expression quantitative trait loci (cis-eQTLs) in the synovium and find that the eQTLs with multiple independent signals have stronger effects and heritability than single independent eQTLs. Integration of genome-wide association studies (GWASs) and eQTLs identifies 84 arthritis related genes, revealing 38 novel genes which have not been reported by previous studies using eQTL data from the GTEx project or immune cells. We further develop a method called eQTac to identify variants that could affect gene expression by affecting chromatin accessibility and identify 1517 regions with potential regulatory function of chromatin accessibility. Altogether, our study provides a comprehensive synovium multi-omics resource for arthritic diseases and gains new insights into the regulation of gene expression.

Integrative high-throughput enhancer surveying and functional verification divulges a YY2-condensed regulatory axis conferring risk for osteoporosis.

The precise roles of chromatin organization at osteoporosis risk loci remain largely elusive. Here, we combined chromatin interaction conformation (Hi-C) profiling and self-transcribing active regulatory region sequencing (STARR-seq) to qualify enhancer activities of prioritized osteoporosis-associated single-nucleotide polymorphisms (SNPs). We identified 319 SNPs with biased allelic enhancer activity effect (baaSNPs) that linked to hundreds of candidate target genes through chromatin interactions across 146 loci. Functional characterizations revealed active epigenetic enrichment for baaSNPs and prevailing osteoporosis-relevant regulatory roles for their chromatin interaction genes. Further motif enrichment and network mapping prioritized several putative, key transcription factors (TFs) controlling osteoporosis binding to baaSNPs. Specifically, we selected one top-ranked TF and deciphered that an intronic baaSNP (rs11202530) could allele-preferentially bind to YY2 to augment PAPSS2 expression through chromatin interactions and promote osteoblast differentiation. Our results underline the roles of TF-mediated enhancer-promoter contacts for osteoporosis, which may help to better understand the intricate molecular regulatory mechanisms underlying osteoporosis risk loci.

Mechanism of Efferocytosis in Determining Ischaemic Stroke Resolution-Diving into Microglia/Macrophage Functions and Therapeutic Modality.

After ischaemic cerebral vascular injury, efferocytosis-a process known as the efficient clearance of apoptotic cells (ACs) by various phagocytes in both physiological and pathological states-is crucial for maintaining central nervous system (CNS) homeostasis and regaining prognosis. The mechanisms of efferocytosis in ischaemic stroke and its influence on preventing inflammation progression from secondary injury were still not fully understood, despite the fact that the fundamental process of efferocytosis has been described in a series of phases, including AC recognition, phagocyte engulfment, and subsequent degradation. The genetic reprogramming of macrophages and brain-resident microglia after an ischaemic stroke has been equated by some researchers to that of the peripheral blood and brain. Based on previous studies, some molecules, such as signal transducer and activator of transcription 6 (STAT6), peroxisome proliferator-activated receptor γ (PPARG), CD300A, and sigma non-opioid intracellular receptor 1 (SIGMAR1), were discovered to be largely associated with aspects of apoptotic cell elimination and accompanying neuroinflammation, such as inflammatory cytokine release, phenotype transformation, and suppressing of antigen presentation. Exacerbated stroke outcomes are brought on by defective efferocytosis and improper modulation of pertinent signalling pathways in blood-borne macrophages and brain microglia, which also results in subsequent tissue inflammatory damage. This review focuses on recent researches which contain a number of recently discovered mechanisms, such as studies on the relationship between benign efferocytosis and the regulation of inflammation in ischaemic stroke, the roles of some risk factors in disease progression, and current immune approaches that aim to promote efferocytosis to treat some autoimmune diseases. Understanding these pathways provides insight into novel pathophysiological processes and fresh characteristics, which can be used to build cerebral ischaemia targeting techniques.

Mitochondrial RNA editing sites affect the phylogenetic reconstruction of gymnosperms.

•RNA editing sites may contain homoplasious signals that cause artifactual inferences in phylogenetic analyses.•Excluding RNA editing sites from gymnosperm mitochondrial genes restored the sister relationship of gnetophytes and Pinaceae.•Phylogenetic analysis based on mitochondrial genomic data should carefully evaluate the impact of RNA editing sites.

High-throughput functional dissection of noncoding SNPs with biased allelic enhancer activity for insulin resistance-relevant phenotypes.

Most of the single-nucleotide polymorphisms (SNPs) associated with insulin resistance (IR)-relevant phenotypes by genome-wide association studies (GWASs) are located in noncoding regions, complicating their functional interpretation. Here, we utilized an adapted STARR-seq to evaluate the regulatory activities of 5,987 noncoding SNPs associated with IR-relevant phenotypes. We identified 876 SNPs with biased allelic enhancer activity effects (baaSNPs) across 133 loci in three IR-relevant cell lines (HepG2, preadipocyte, and A673), which showed pervasive cell specificity and significant enrichment for cell-specific open chromatin regions or enhancer-indicative markers (H3K4me1, H3K27ac). Further functional characterization suggested several transcription factors (TFs) with preferential allelic binding to baaSNPs. We also incorporated multi-omics data to prioritize 102 candidate regulatory target genes for baaSNPs and revealed prevalent long-range regulatory effects and cell-specific IR-relevant biological functional enrichment on them. Specifically, we experimentally verified the distal regulatory mechanism at IRS1 locus, in which rs952227-A reinforces IRS1 expression by long-range chromatin interaction and preferential binding to the transcription factor HOXC6 to augment the enhancer activity. Finally, based on our STARR-seq screening data, we predicted the enhancer activity of 227,343 noncoding SNPs associated with IR-relevant phenotypes (fasting insulin adjusted for BMI, HDL cholesterol, and triglycerides) from the largest available GWAS summary statistics. We further provided an open resource (http://www.bigc.online/fnSNP-IR) for better understanding genetic regulatory mechanisms of IR-relevant phenotypes.

Assessment of bidirectional relationships between brain imaging-derived phenotypes and stroke: a Mendelian randomization study.

BACKGROUND: Stroke is a major cause of mortality and long-term disability worldwide. Whether the associations between brain imaging-derived phenotypes (IDPs) and stroke are causal is uncertain. METHODS: We performed two-sample bidirectional Mendelian randomization (MR) analyses to explore the causal associations between IDPs and stroke. Summary data of 587 brain IDPs (up to 33,224 individuals) from the UK Biobank and five stroke types (sample size range from 301,663 to 446,696, case number range from 5,386 to 40,585) from the MEGASTROKE consortium were used. RESULTS: Forward MR indicated 14 IDPs belong to projection fibers or association fibers were associated with stroke. For example, higher genetically determined mean diffusivity (MD) in the right external capsule was causally associated with an increased risk of small vessel stroke (IVW OR = 2.76, 95% CI 2.07 to 3.68, P = 5.87 × 10-12). Reverse MR indicated that genetically determined higher risk of any ischemic stroke was associated with increased isotropic or free water volume fraction (ISOVF) in body of corpus callosum (IVW ß = 0.23, 95% CI 0.14 to 0.33, P = 3.22 × 10-7). This IDP is a commissural fiber and it is not included in the IDPs identified by forward MR. CONCLUSIONS: We identified 14 IDPs with statistically significant evidence of causal effects on stroke or stroke subtypes. We also identified potential causal effects of stroke on one IDP of commissural fiber. These findings might guide further work toward identifying preventative strategies at the brain imaging levels.

Enhancer variants on chromosome 2p14 regulating SPRED2 and ACTR2 act as a signal amplifier to protect against rheumatoid arthritis.

Genome-wide association studies (GWASs) have repeatedly reported multiple non-coding single-nucleotide polymorphisms (SNPs) at 2p14 associated with rheumatoid arthritis (RA), but their functional roles in the pathological mechanisms of RA remain to be explored. In this study, we integrated a series of bioinformatics and functional experiments and identified three intronic RA SNPs (rs1876518, rs268131, and rs2576923) within active enhancers that can regulate the expression of SPRED2 directly. At the same time, SPRED2 and ACTR2 influence each other as a positive feedback signal amplifier to strengthen the protective role in RA by inhibiting the migration and invasion of rheumatoid fibroblast-like synoviocytes (FLSs). In particular, the transcription factor CEBPB preferentially binds to the rs1876518-T allele to increase the expression of SPRED2 in FLSs. Our findings decipher the molecular mechanisms behind the GWAS signals at 2p14 for RA and emphasize SPRED2 as a potential candidate gene for RA, providing a potential target and direction for precise treatment of RA.

Surgical methods influence on the risk of anastomotic fistula after pancreaticoduodenectomy: a systematic review and network meta-analysis.

BACKGROUND: Pancreaticoduodenectomy is the first choice surgical intervention for the radical treatment of pancreatic tumors. However, an anastomotic fistula is a common complication after pancreaticoduodenectomy with a high mortality rate. With the development of minimally invasive surgery, open pancreaticoduodenectomy (OPD), laparoscopic pancreaticoduodenectomy (LPD), and robotic pancreaticoduodenectomy (RPD) are gaining interest. But the impact of these surgical methods on the risk of anastomosis has not been confirmed. Therefore, we aimed to integrate relevant clinical studies and explore the effects of these three surgical methods on the occurrence of anastomotic fistula after pancreaticoduodenectomy. METHODS: A systematic literature search was conducted for studies reporting the RPD, LPD, and OPD. Network meta-analysis of postoperative anastomotic fistula (Pancreatic fistula, biliary leakage, gastrointestinal fistula) was performed. RESULTS: Sixty-five studies including 10,026 patients were included in the network meta-analysis. The rank of risk probability of pancreatic fistula for RPD (0.00) was better than LPD (0.37) and OPD (0.62). Thus, the analysis suggests the rank of risk of the postoperative pancreatic fistula for RPD, LPD, and OPD. The rank of risk probability for biliary leakage was similar for RPD (0.15) and LPD (0.15), and both were better than OPD (0.68). CONCLUSIONS: This network meta-analysis provided ranking for three different types of pancreaticoduodenectomy. The RPD and LPD can effectively improve the quality of surgery and are safe as well as feasible for OPD.

Multitissue Integrative Analysis Identifies Susceptibility Genes for Atopic Dermatitis.

Atopic dermatitis (AD) is a chronic relapsing inflammatory skin disease with multiple environmental and genetic factors involved in its etiology. Although lots of genetic loci associated with AD have been reported by GWASs, only a small part of phenotypic variations can be explained. To identify additional susceptibility genes on AD, we conducted a large-scale transcriptome-wide association study using a joint-tissue imputation approach in ∼840,000 European individuals combined with six precomputed gene expression weights of four AD-relevant tissues, including skin fibroblast, lymphocyte, and whole blood. The Mendelian randomization causal inference analysis was performed to estimate the causal effect of transcriptome-wide association study‒identified genes. We identified 51 genes significantly associated with AD after Bonferroni corrections, and 19 genes showed putatively causal associations such as an established gene FLG (P = 3.98 × 10‒10) and seven genes that have not been implicated in previous transcriptome-wide association studies, such as AQP3 (P = 4.43 × 10‒7) and PDCD1 (P = 7.66 × 10‒7). Among them, four genes (AQP3, PDCD1, ADCY3, and DOLPP1) were further supported in differential expression analyses or the Mouse Genome Informatics database. Overall, our study identified susceptibility genes associated with AD, providing, to our knowledge, previously unreported clues in revealing the genetic mechanisms in AD.

Lineage-selective super enhancers mediate core regulatory circuitry during adipogenic and osteogenic differentiation of human mesenchymal stem cells.

Human mesenchymal stem cells (hMSCs) can be differentiated into osteoblasts and adipocytes. During these processes, super enhancers (SEs) play important roles. Here, we performed comprehensive characterization of the SEs changes associated with adipogenic and osteogenic differentiation of hMSCs, and revealed that SEs changed more dramatically compared with typical enhancers. We identified a set of lineage-selective SEs, whose target genes were enriched with cell type-specific functions. Functional experiments in lineage-selective SEs demonstrated their specific roles in directed differentiation of hMSCs. We also found that some key transcription factors regulated by lineage-selective SEs could form core regulatory circuitry (CRC) to regulate each other's expression and control the hMSCs fate determination. In addition, we found that GWAS SNPs of osteoporosis and obesity were significantly enriched in osteoblasts-selective SEs or adipocytes-selective SEs, respectively. Taken together, our studies unveiled important roles of lineage-selective SEs in hMSCs differentiation into osteoblasts and adipocytes.

Mendelian randomization analyses support causal relationships between brain imaging-derived phenotypes and risk of psychiatric disorders.

Observational studies have reported the correlations between brain imaging-derived phenotypes (IDPs) and psychiatric disorders; however, whether the relationships are causal is uncertain. We conducted bidirectional two-sample Mendelian randomization (MR) analyses to explore the causalities between 587 reliable IDPs (N = 33,224 individuals) and 10 psychiatric disorders (N = 9,725 to 161,405). We identified nine IDPs for which there was evidence of a causal influence on risk of schizophrenia, anorexia nervosa and bipolar disorder. For example, 1 s.d. increase in the orientation dispersion index of the forceps major was associated with 32% lower odds of schizophrenia risk. Reverse MR indicated that only genetically predicted schizophrenia was positively associated with two IDPs, the cortical surface area and the volume of the right pars orbitalis. We established the BrainMR database ( http://www.bigc.online/BrainMR/ ) to share our results. Our findings provide potential strategies for the prediction and intervention for psychiatric disorder risk at the brain-imaging level.

Bidirectional two-sample Mendelian randomization analysis identifies causal associations between relative carbohydrate intake and depression.

Growing evidence suggests that relative carbohydrate intake affects depression; however, the association between carbohydrates and depression remains controversial. To test this, we performed a two-sample bidirectional Mendelian randomization (MR) analysis using genetic variants associated with relative carbohydrate intake (N = 268,922) and major depressive disorder (N = 143,265) from the largest available genome-wide association studies. MR evidence suggested a causal relationship between higher relative carbohydrate intake and lower depression risk (odds ratio, 0.42 for depression per one-standard-deviation increment in relative carbohydrate intake; 95% confidence interval, 0.28 to 0.62; P = 1.49 × 10-5). Multivariable MR indicated that the protective effect of relative carbohydrate intake on depression persisted after conditioning on other diet compositions. The mediation analysis via two-step MR showed that this effect was partly mediated by body mass index, with a mediated proportion of 15.4% (95% confidence interval, 6.7% to 24.1%). These findings may inform prevention strategies and interventions directed towards relative carbohydrate intake and depression.

Lineage-specific rearrangement of chromatin loops and epigenomic features during adipocytes and osteoblasts commitment.

Human mesenchymal stem cells (hMSCs) can be differentiated into adipocytes and osteoblasts. The processes are driven by the rewiring of chromatin architectures and transcriptomic/epigenomic changes. Here, we induced hMSCs to adipogenic and osteogenic differentiation, and performed 2 kb resolution Hi-C experiments for chromatin loops detection. We also generated matched RNA-seq, ChIP-seq and ATAC-seq data for integrative analysis. After comprehensively comparing adipogenesis and osteogenesis, we quantitatively identified lineage-specific loops and screened out lineage-specific enhancers and open chromatin. We reveal that lineage-specific loops can activate gene expression and facilitate cell commitment through combining enhancers and accessible chromatin in a lineage-specific manner. We finally proposed loop-mediated regulatory networks and identified the controlling factors for adipocytes and osteoblasts determination. Functional experiments validated the lineage-specific regulation networks towards IRS2 and RUNX2 that are associated with adipogenesis and osteogenesis, respectively. These results are expected to help better understand the chromatin conformation determinants of hMSCs fate commitment.

ExsgRNA: reduce off-target efficiency by on-target mismatched sgRNA.

Clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 gene editing technology has been widely used to facilitate efficient genome editing. Current popular sgRNA design tools only consider the sgRNA perfectly matched to the target site and provide the results without any on-target mismatch. We suppose taking on-target gRNA-DNA mismatches into consideration might provide better sgRNA with similar binding activity and reduced off-target sites. Here, we trained a seq2seq-attention model with feedback-loop architecture, to automatically generate sgRNAs with on-target mismatches. Dual-luciferase reporter experiment showed that multiple sgRNAs with three mismatches could achieve the 80% of the relative activity of the perfect matched sgRNA. Meanwhile, it could reduce the number of off-target sites using sgRNAs with on-target mismatches. Finally, we provided a freely accessible web server sgRNA design tool named ExsgRNA. Users could submit their target sequence to this server and get optimal sgRNAs with less off-targets and similar on-target activity compared with the perfect-matched sgRNA.

An Osteoporosis Susceptibility Allele at 11p15 Regulates SOX6 Expression by Modulating TCF4 Chromatin Binding.

Osteoporosis is an age-related complex disease clinically diagnosed with bone mineral density (BMD). Although several genomewide association studies (GWASs) have discovered multiple noncoding genetic variants at 11p15 influencing osteoporosis risk, the functional mechanisms of these variants remain unknown. Through integrating bioinformatics and functional experiments, a potential functional single-nucleotide polymorphism (SNP; rs1440702) located in an enhancer element was identified and the A allele of rs1440702 acted as an allelic specificities enhancer to increase its distal target gene SOX6 (~600 Kb upstream) expression, which plays a key role in bone formation. We also validated this long-range regulation via conducting chromosome conformation capture (3C) assay. Furthermore, we demonstrated that SNP rs1440702 with a risk allele (rs1440702-A) could increase the activity of the enhancer element by altering the binding affinity of the transcription factor TCF4, resulting in the upregulation expression of SOX6 gene. Collectively, our integrated analyses revealed how the noncoding genetic variants (rs1440702) affect osteoporosis predisposition via long-range gene regulatory mechanisms and identified its target gene SOX6 for downstream biomarker and drug development. © 2022 American Society for Bone and Mineral Research (ASBMR).

Individualized pathway activity algorithm identifies oncogenic pathways in pan-cancer analysis.

BACKGROUND: Accumulative evidences have shown that dysregulation of biological pathways contributed to the initiation and progression of malignant tumours. Several methods for pathway activity measurement have been proposed, but they are restricted to making comparisons between groups or sensitive to experimental batch effects. METHODS: We introduced a novel method for individualized pathway activity measurement (IPAM) that is based on the ranking of gene expression levels in individual sample. Taking advantage of IPAM, we calculated the pathway activity of 318 pathways from KEGG database in the 10528 tumour/normal samples of 33 cancer types from TCGA to identify characteristic dysregulated pathways among different cancer types. FINDINGS: IPAM precisely quantified the level of activity of each pathway in pan-cancer analysis and exhibited better performance in cancer classification and prognosis prediction over five widely used tools. The average ROC-AUC of cancer diagnostic model using tumour-educated platelets (TEPs) reached 92.84%, suggesting the potential of our algorithm in early diagnosis of cancer. We identified several pathways significantly deregulated and associated with patient survival in a large fraction of cancer types, such as tyrosine metabolism, fatty acid degradation, cell cycle, p53 signalling pathway and DNA replication. We also confirmed the dominant role of metabolic pathways in cancer pathway dysregulation and identified the driving factors of specific pathway dysregulation, such as PPARA for branched-chain amino acid metabolism and NR1I2, NR1I3 for fatty acid metabolism. INTERPRETATION: Our study will provide novel clues for understanding the pathological mechanisms of cancer, ultimately paving the way for personalized medicine of cancer. FUNDING: A full list of funding can be found in the Acknowledgements section.

MetaDecoder: a novel method for clustering metagenomic contigs.

BACKGROUND: Clustering the metagenomic contigs into potential genomes is a key step to investigate the functional roles of microbial populations. Existing algorithms have achieved considerable success with simulated or real sequencing datasets. However, accurately classifying contigs from complex metagenomes is still a challenge. RESULTS: We introduced a novel clustering algorithm, MetaDecoder, which can classify metagenomic contigs based on the frequencies of k-mers and coverages. MetaDecoder was built as a two-layer model with the first layer being a GPU-based modified Dirichlet process Gaussian mixture model (DPGMM), which controls the weight of each DPGMM cluster to avoid over-segmentation by dynamically dissolving contigs in small clusters and reassigning them to the remaining clusters. The second layer comprises a semi-supervised k-mer frequency probabilistic model and a modified Gaussian mixture model for modeling the coverage based on single copy marker genes. Benchmarks on simulated and real-world datasets demonstrated that MetaDecoder can be served as a promising approach for effectively clustering metagenomic contigs. CONCLUSIONS: In conclusion, we developed the GPU-based MetaDecoder for effectively clustering metagenomic contigs and reconstructing microbial communities from microbial data. Applying MetaDecoder on both simulated and real-world datasets demonstrated that it could generate more complete clusters with lower contamination. Using MetaDecoder, we identified novel high-quality genomes and expanded the existing catalog of bacterial genomes. Video Abstract.