Investigating shared genetic architecture between inflammatory bowel diseases and primary biliary cholangitis.

Background & Aims: Inflammatory bowel disease (IBD) is commonly associated with extraintestinal complications, including autoimmune liver disease. The co-occurrence of IBD and primary biliary cholangitis (PBC) has been increasingly observed, but the underlying relationship between these conditions remains unclear. Methods: Using summary statistics from genome-wide association studies (GWAS), we investigated the causal effects between PBC and IBD, including Crohn's disease (CD) and ulcerative colitis (UC). We also analyzed the shared genetic architecture between IBD and PBC using data from GWAS, bulk-tissue RNA sequencing, and single cell RNA sequencing, and explored potential functional genes. Result: There was a strong positive genetic correlation between PBC and IBD (linkage disequilibrium score regression: rg = 0.2249, p = 3.38 × 10-5). Cross-trait analysis yielded 10 shared-risk single nucleotide polymorphisms (SNPs), as well as nine novel SNPs, which were associated with both traits. Using Mendelian randomization, a stable causal effect was established of PBC on IBD. Genetically predicted PBC was found to have a risk effect on IBD (1.105; 95% CI: 1.058-1.15; p = 1.16 × 10-10), but not vice versa. Shared tissue-specific heritability enrichment was identified for PBC and IBD (including CD and UC) in lung, spleen, and whole-blood samples. Furthermore, shared enrichment was observed of specific cell types (T cells, B cells, and natural killer cells) and their subtypes. Nine functional genes were identified based on summary statistics-based Mendelian randomization. Conclusions: This study detected shared genetic architecture between IBD and PBC and demonstrated a stable causal relationship of genetically predicted PBC on the risk of IBD. These findings shed light on the biological basis of comorbidity between IBD and PBC, and have important implications for intervention and treatment targets of these two diseases simultaneously. Impact and Implications: The discovery of novel shared single nucleotide polymorphisms (SNPs) and functional genes provides insights into the common targets between inflammatory bowel disease (IBD) and primary biliary cholangitis (PBC), serving as a basis for new drug development and contributing to the study of disease pathogenesis. Additionally, the established significant causality and genetic correlation underscore the importance of clinical intervention in preventing the comorbidity of IBD and PBC. The enrichment of SNP heritability in specific tissues and cell types reveals the role of immune factors in the potential disease mechanisms shared between IBD and PBC. This stimulates further research on potential interventions and could lead to the development of new targets for immune-based therapies.

ARGNet: using deep neural networks for robust identification and classification of antibiotic resistance genes from sequences.

BACKGROUND: Emergence of antibiotic resistance in bacteria is an important threat to global health. Antibiotic resistance genes (ARGs) are some of the key components to define bacterial resistance and their spread in different environments. Identification of ARGs, particularly from high-throughput sequencing data of the specimens, is the state-of-the-art method for comprehensively monitoring their spread and evolution. Current computational methods to identify ARGs mainly rely on alignment-based sequence similarities with known ARGs. Such approaches are limited by choice of reference databases and may potentially miss novel ARGs. The similarity thresholds are usually simple and could not accommodate variations across different gene families and regions. It is also difficult to scale up when sequence data are increasing. RESULTS: In this study, we developed ARGNet, a deep neural network that incorporates an unsupervised learning autoencoder model to identify ARGs and a multiclass classification convolutional neural network to classify ARGs that do not depend on sequence alignment. This approach enables a more efficient discovery of both known and novel ARGs. ARGNet accepts both amino acid and nucleotide sequences of variable lengths, from partial (30-50 aa; 100-150 nt) sequences to full-length protein or genes, allowing its application in both target sequencing and metagenomic sequencing. Our performance evaluation showed that ARGNet outperformed other deep learning models including DeepARG and HMD-ARG in most of the application scenarios especially quasi-negative test and the analysis of prediction consistency with phylogenetic tree. ARGNet has a reduced inference runtime by up to 57% relative to DeepARG. CONCLUSIONS: ARGNet is flexible, efficient, and accurate at predicting a broad range of ARGs from the sequencing data. ARGNet is freely available at https://github.com/id-bioinfo/ARGNet , with an online service provided at https://ARGNet.hku.hk . Video Abstract.

Development and validation of risk prediction model for recurrent cardiovascular events among Chinese: the Personalized CARdiovascular DIsease risk Assessment for Chinese model.

Aims: Cardiovascular disease (CVD) is a leading cause of mortality, especially in developing countries. This study aimed to develop and validate a CVD risk prediction model, Personalized CARdiovascular DIsease risk Assessment for Chinese (P-CARDIAC), for recurrent cardiovascular events using machine learning technique. Methods and results: Three cohorts of Chinese patients with established CVD were included if they had used any of the public healthcare services provided by the Hong Kong Hospital Authority (HA) since 2004 and categorized by their geographical locations. The 10-year CVD outcome was a composite of diagnostic or procedure codes with specific International Classification of Diseases, Ninth Revision, Clinical Modification. Multivariate imputation with chained equations and XGBoost were applied for the model development. The comparison with Thrombolysis in Myocardial Infarction Risk Score for Secondary Prevention (TRS-2°P) and Secondary Manifestations of ARTerial disease (SMART2) used the validation cohorts with 1000 bootstrap replicates. A total of 48 799, 119 672 and 140 533 patients were included in the derivation and validation cohorts, respectively. A list of 125 risk variables were used to make predictions on CVD risk, of which 8 classes of CVD-related drugs were considered interactive covariates. Model performance in the derivation cohort showed satisfying discrimination and calibration with a C statistic of 0.69. Internal validation showed good discrimination and calibration performance with C statistic over 0.6. The P-CARDIAC also showed better performance than TRS-2°P and SMART2. Conclusion: Compared with other risk scores, the P-CARDIAC enables to identify unique patterns of Chinese patients with established CVD. We anticipate that the P-CARDIAC can be applied in various settings to prevent recurrent CVD events, thus reducing the related healthcare burden.

Sleep patterns, genetic susceptibility, and digestive diseases: A large-scale longitudinal cohort study.

BACKGROUND: Sleep problems are prevalent. However, the impact of sleep patterns on digestive diseases remains uncertain. Moreover, the interaction between sleep patterns and genetic predisposition with digestive diseases has not been comprehensively explored. METHODS: 410,586 participants from UK Biobank with complete sleep information were included in the analysis. Sleep patterns were measured by sleep scores as the primary exposure, based on five healthy sleep behaviors. Individual sleep behaviors were secondary exposures. Genetic risk of the digestive diseases was characterized by polygenic risk score. Primary outcome was incidence of 16 digestive diseases. RESULTS: Healthy sleep scores showed dose-response associations with reduced risks of digestive diseases. Compared to participants scoring 0-1, those scoring 5 showed a 28% reduced risk of any digestive disease, including a 50% decrease in irritable bowel syndrome, 37% in non-alcoholic fatty liver disease, 35% in peptic ulcer, 34% in dyspepsia, 32% in gastroesophageal reflux disease, 28% in constipation, 25% in diverticulosis, 24% in severe liver disease, and 18% in gallbladder disease, whereas no correlation was observed with inflammatory bowel disease and pancreatic disease. Participants with poor sleep and high genetic risk exhibited approximately a 60% increase in the risk of digestive diseases. A healthy sleep pattern is linked to lower digestive disease risk in participants of all genetic risk levels. CONCLUSIONS: In this large population-based cohort, a healthy sleep pattern was associated with reduced risk of digestive diseases, regardless of the genetic susceptibility. Our findings underscore the potential impact of healthy sleep traits in mitigating the risk of digestive diseases.

Adverse clinical outcomes and immunosuppressive microenvironment of RHO-GTPase activation pattern in hepatocellular carcinoma.

BACKGROUND: Emerging evidence suggests that Rho GTPases play a crucial role in tumorigenesis and metastasis, but their involvement in the tumor microenvironment (TME) and prognosis of hepatocellular carcinoma (HCC) is not well understood. METHODS: We aim to develop a tumor prognosis prediction system called the Rho GTPases-related gene score (RGPRG score) using Rho GTPase signaling genes and further bioinformatic analyses. RESULTS: Our work found that HCC patients with a high RGPRG score had significantly worse survival and increased immunosuppressive cell fractions compared to those with a low RGPRG score. Single-cell cohort analysis revealed an immune-active TME in patients with a low RGPRG score, with strengthened communication from T/NK cells to other cells through MIF signaling networks. Targeting these alterations in TME, the patients with high RGPRG score have worse immunotherapeutic outcomes and decreased survival time in the immunotherapy cohort. Moreover, the RGPRG score was found to be correlated with survival in 27 other cancers. In vitro experiments confirmed that knockdown of the key Rho GTPase-signaling biomarker SFN significantly inhibited HCC cell proliferation, invasion, and migration. CONCLUSIONS: This study provides new insight into the TME features and clinical use of Rho GTPase gene pattern at the bulk-seq and single-cell level, which may contribute to guiding personalized treatment and improving clinical outcome in HCC.

ClusterV-Web: a user-friendly tool for profiling HIV quasispecies and generating drug resistance reports from nanopore long-read data.

Summary: Third-generation long-read sequencing is an increasingly utilized technique for profiling human immunodeficiency virus (HIV) quasispecies and detecting drug resistance mutations due to its ability to cover the entire viral genome in individual reads. Recently, the ClusterV tool has demonstrated accurate detection of HIV quasispecies from Nanopore long-read sequencing data. However, the need for scripting skills and a computational environment may act as a barrier for many potential users. To address this issue, we have introduced ClusterV-Web, a user-friendly web-based application that enables easy configuration and execution of ClusterV, both remotely and locally. Our tool provides interactive tables and data visualizations to aid in the interpretation of results. This development is expected to democratize access to long-read sequencing data analysis, enabling a wider range of researchers and clinicians to efficiently profile HIV quasispecies and detect drug resistance mutations. Availability and implementation: ClusterV-Web is freely available and open source, with detailed documentation accessible at http://www.bio8.cs.hku.hk/ClusterVW/. The standalone Docker image and source code are also available at https://github.com/HKU-BAL/ClusterV-Web.

Integrated modeling framework reveals co-regulation of transcription factors, miRNAs and lncRNAs on cardiac developmental dynamics.

AIMS: Dissecting complex interactions among transcription factors (TFs), microRNAs (miRNAs) and long noncoding RNAs (lncRNAs) are central for understanding heart development and function. Although computational approaches and platforms have been described to infer relationships among regulatory factors and genes, current approaches do not adequately account for how highly diverse, interacting regulators that include noncoding RNAs (ncRNAs) control cardiac gene expression dynamics over time. METHODS: To overcome this limitation, we devised an integrated framework, cardiac gene regulatory modeling (CGRM) that integrates LogicTRN and regulatory component analysis bioinformatics modeling platforms to infer complex regulatory mechanisms. We then used CGRM to identify and compare the TF-ncRNA gene regulatory networks that govern early- and late-stage cardiomyocytes (CMs) generated by in vitro differentiation of human pluripotent stem cells (hPSC) and ventricular and atrial CMs isolated during in vivo human cardiac development. RESULTS: Comparisons of in vitro versus in vivo derived CMs revealed conserved regulatory networks among TFs and ncRNAs in early cells that significantly diverged in late staged cells. We report that cardiac genes ("heart targets") expressed in early-stage hPSC-CMs are primarily regulated by MESP1, miR-1, miR-23, lncRNAs NEAT1 and MALAT1, while GATA6, HAND2, miR-200c, NEAT1 and MALAT1 are critical for late hPSC-CMs. The inferred TF-miRNA-lncRNA networks regulating heart development and contraction were similar among early-stage CMs, among individual hPSC-CM datasets and between in vitro and in vivo samples. However, genes related to apoptosis, cell cycle and proliferation, and transmembrane transport showed a high degree of divergence between in vitro and in vivo derived late-stage CMs. Overall, late-, but not early-stage CMs diverged greatly in the expression of "heart target" transcripts and their regulatory mechanisms. CONCLUSIONS: In conclusion, we find that hPSC-CMs are regulated in a cell autonomous manner during early development that diverges significantly as a function of time when compared to in vivo derived CMs. These findings demonstrate the feasibility of using CGRM to reveal dynamic and complex transcriptional and posttranscriptional regulatory interactions that underlie cell directed versus environment-dependent CM development. These results with in vitro versus in vivo derived CMs thus establish this approach for detailed analyses of heart disease and for the analysis of cell regulatory systems in other biomedical fields.

Boosting variant-calling performance with multi-platform sequencing data using Clair3-MP.

BACKGROUND: With the continuous advances in third-generation sequencing technology and the increasing affordability of next-generation sequencing technology, sequencing data from different sequencing technology platforms is becoming more common. While numerous benchmarking studies have been conducted to compare variant-calling performance across different platforms and approaches, little attention has been paid to the potential of leveraging the strengths of different platforms to optimize overall performance, especially integrating Oxford Nanopore and Illumina sequencing data. RESULTS: We investigated the impact of multi-platform data on the performance of variant calling through carefully designed experiments with a deep learning-based variant caller named Clair3-MP (Multi-Platform). Through our research, we not only demonstrated the capability of ONT-Illumina data for improved variant calling, but also identified the optimal scenarios for utilizing ONT-Illumina data. In addition, we revealed that the improvement in variant calling using ONT-Illumina data comes from an improvement in difficult genomic regions, such as the large low-complexity regions and segmental and collapse duplication regions. Moreover, Clair3-MP can incorporate reference genome stratification information to achieve a small but measurable improvement in variant calling. Clair3-MP is accessible as an open-source project at: https://github.com/HKU-BAL/Clair3-MP . CONCLUSIONS: These insights have important implications for researchers and practitioners alike, providing valuable guidance for improving the reliability and efficiency of genomic analysis in diverse applications.

Long-Read Sequencing with Hierarchical Clustering for Antiretroviral Resistance Profiling of Mixed Human Immunodeficiency Virus Quasispecies.

BACKGROUND: HIV infections often develop drug resistance mutations (DRMs), which can increase the risk of virological failure. However, it has been difficult to determine if minor mutations occur in the same genome or in different virions using Sanger sequencing and short-read sequencing methods. Oxford Nanopore Technologies (ONT) sequencing may improve antiretroviral resistance profiling by allowing for long-read clustering. METHODS: A new ONT sequencing-based method for profiling DRMs in HIV quasispecies was developed and validated. The method used hierarchical clustering of long amplicons that cover regions associated with different types of antiretroviral drugs. A gradient series of an HIV plasmid and 2 plasma samples was prepared to validate the clustering performance. The ONT results were compared to those obtained with Sanger sequencing and Illumina sequencing in 77 HIV-positive plasma samples to evaluate the diagnostic performance. RESULTS: In the validation study, the abundance of detected quasispecies was concordant with the predicted result with the R2 of > 0.99. During the diagnostic evaluation, 59/77 samples were successfully sequenced for DRMs. Among 18 failed samples, 17 were below the limit of detection of 303.9 copies/µL. Based on the receiver operating characteristic analysis, the ONT workflow achieved an F1 score of 0.96 with a cutoff of 0.4 variant allele frequency. Four cases were found to have quasispecies with DRMs, in which 2 harbored quasispecies with more than one class of DRMs. Treatment modifications were recommended for these cases. CONCLUSIONS: Long-read sequencing coupled with hierarchical clustering could differentiate the quasispecies resistance profiles in HIV-infected samples, providing a clearer picture for medical care.

Evaluation of Mycobacterium tuberculosis enrichment in metagenomic samples using ONT adaptive sequencing and amplicon sequencing for identification and variant calling.

Sensitive detection of Mycobacterium tuberculosis (TB) in small percentages in metagenomic samples is essential for microbial classification and drug resistance prediction. However, traditional methods, such as bacterial culture and microscopy, are time-consuming and sometimes have limited TB detection sensitivity. Oxford nanopore technologies (ONT) MinION sequencing allows rapid and simple sample preparation for sequencing. Its recently developed adaptive sequencing selects reads from targets while allowing real-time base-calling to achieve sequence enrichment or depletion during sequencing. Another common enrichment method is PCR amplification of the target TB genes. In this study, we compared both methods using ONT MinION sequencing for TB detection and variant calling in metagenomic samples using both simulation runs and those with synthetic and patient samples. We found that both methods effectively enrich TB reads from a high percentage of human (95%) and other microbial DNA. Adaptive sequencing with readfish and UNCALLDE achieved a 3.9-fold and 2.2-fold enrichment compared to the control run. We provide a simple automatic analysis framework to support the detection of TB for clinical use, openly available at https://github.com/HKU-BAL/ONT-TB-NF . Depending on the patient's medical condition and sample type, we recommend users evaluate and optimize their workflow for different clinical specimens to improve the detection limit.

Ultra-low-coverage genome-wide association study-insights into gestational age using 17,844 embryo samples with preimplantation genetic testing.

BACKGROUND: Very low-coverage (0.1 to 1×) whole genome sequencing (WGS) has become a promising and affordable approach to discover genomic variants of human populations for genome-wide association study (GWAS). To support genetic screening using preimplantation genetic testing (PGT) in a large population, the sequencing coverage goes below 0.1× to an ultra-low level. However, the feasibility and effectiveness of ultra-low-coverage WGS (ulcWGS) for GWAS remains undetermined. METHODS: We built a pipeline to carry out analysis of ulcWGS data for GWAS. To examine its effectiveness, we benchmarked the accuracy of genotype imputation at the combination of different coverages below 0.1× and sample sizes from 2000 to 16,000, using 17,844 embryo PGT samples with approximately 0.04× average coverage and the standard Chinese sample HG005 with known genotypes. We then applied the imputed genotypes of 1744 transferred embryos who have gestational ages and complete follow-up records to GWAS. RESULTS: The accuracy of genotype imputation under ultra-low coverage can be improved by increasing the sample size and applying a set of filters. From 1744 born embryos, we identified 11 genomic risk loci associated with gestational ages and 166 genes mapped to these loci according to positional, expression quantitative trait locus, and chromatin interaction strategies. Among these mapped genes, CRHBP, ICAM1, and OXTR were more frequently reported as preterm birth related. By joint analysis of gene expression data from previous studies, we constructed interrelationships of mainly CRHBP, ICAM1, PLAGL1, DNMT1, CNTLN, DKK1, and EGR2 with preterm birth, infant disease, and breast cancer. CONCLUSIONS: This study not only demonstrates that ulcWGS could achieve relatively high accuracy of adequate genotype imputation and is capable of GWAS, but also provides insights into the associations between gestational age and genetic variations of the fetal embryos from Chinese population.

Duet: SNP-assisted structural variant calling and phasing using Oxford nanopore sequencing.

BACKGROUND: Whole genome sequencing using the long-read Oxford Nanopore Technologies (ONT) MinION sequencer provides a cost-effective option for structural variant (SV) detection in clinical applications. Despite the advantage of using long reads, however, accurate SV calling and phasing are still challenging. RESULTS: We introduce Duet, an SV detection tool optimized for SV calling and phasing using ONT data. The tool uses novel features integrated from both SV signatures and single-nucleotide polymorphism signatures, which can accurately distinguish SV haplotype from a false signal. Duet was benchmarked against state-of-the-art tools on multiple ONT sequencing datasets of sequencing coverage ranging from 8× to 40×. At low sequencing coverage of 8×, Duet performs better than all other tools in SV calling, SV genotyping and SV phasing. When the sequencing coverage is higher (20× to 40×), the F1-score for SV phasing is further improved in comparison to the performance of other tools, while its performance of SV genotyping and SV calling remains higher than other tools. CONCLUSION: Duet can perform accurate SV calling, SV genotyping and SV phasing using low-coverage ONT data, making it very useful for low-coverage genomes. It has great performance when scaled to high-coverage genomes, which is adaptable to various clinical applications. Duet is open source and is available at https://github.com/yekaizhou/duet .

Assembly-free discovery of human novel sequences using long reads.

DNA sequences that are absent in the human reference genome are classified as novel sequences. The discovery of these missed sequences is crucial for exploring the genomic diversity of populations and understanding the genetic basis of human diseases. However, various DNA lengths of reads generated from different sequencing technologies can significantly affect the results of novel sequences. In this work, we designed an assembly-free novel sequence (AF-NS) approach to identify novel sequences from Oxford Nanopore Technology long reads. Among the newly detected sequences using AF-NS, more than 95% were omitted from those using long-read assemblers and 85% were not present in short reads of Illumina. We identified the common novel sequences among all the samples and revealed their association with the binding motifs of transcription factors. Regarding the placements of the novel sequences, we found about 70% enriched in repeat regions and generated 430 for one specific subpopulation that might be related to their evolution. Our study demonstrates the advance of the assembly-free approach to capture more novel sequences over other assembler based methods. Combining the long-read data with powerful analytical methods can be a robust way to improve the completeness of novel sequences.

Clair3-trio: high-performance Nanopore long-read variant calling in family trios with trio-to-trio deep neural networks.

Accurate identification of genetic variants from family child-mother-father trio sequencing data is important in genomics. However, state-of-the-art approaches treat variant calling from trios as three independent tasks, which limits their calling accuracy for Nanopore long-read sequencing data. For better trio variant calling, we introduce Clair3-Trio, the first variant caller tailored for family trio data from Nanopore long-reads. Clair3-Trio employs a Trio-to-Trio deep neural network model, which allows it to input the trio sequencing information and output all of the trio's predicted variants within a single model to improve variant calling. We also present MCVLoss, a novel loss function tailor-made for variant calling in trios, leveraging the explicit encoding of the Mendelian inheritance. Clair3-Trio showed comprehensive improvement in experiments. It predicted far fewer Mendelian inheritance violation variations than current state-of-the-art methods. We also demonstrated that our Trio-to-Trio model is more accurate than competing architectures. Clair3-Trio is accessible as a free, open-source project at https://github.com/HKU-BAL/Clair3-Trio.

Streptococcus oriscaviae sp. nov. Infection Associated with Guinea Pigs.

Pet bite-related infections are commonly caused by the pet's oral flora transmitted to the animal handlers through the bite wounds. In this study, we isolated a streptococcus, HKU75T, in pure culture from the purulent discharge collected from a guinea pig bite wound in a previously healthy young patient. HKU75T was alpha-hemolytic on sheep blood agar and agglutinated with Lancefield group D and group G antisera. API 20 STREP showed that the most likely identity for HKU75T was S. suis I with 85.4% confidence while Vitek 2 showed that HKU75T was unidentifiable. MALDI-TOF MS identified HKU75T as Streptococcus suis (score of 1.86 only). 16S rRNA gene sequencing showed that HKU75T was most closely related to S. parasuis (98.3% nucleotide identity), whereas partial groEL and rpoB gene sequencing showed that it was most closely related to S. suis (81.8% and 89.8% nucleotide identity respectively). Whole genome sequencing and intergenomic distance determined by ANI revealed that there was <85% identity between the genome of HKU75T and those of all other known Streptococcus species. Genome classification using concatenated sequences of 92 bacterial core genes showed that HKU75T belonged to the Suis group. groEL gene sequences identical to that of HKU75T could be directly amplified from the oral cavities of the two guinea pigs owned by the patient. HKU75T is a novel Streptococcus species, which we propose to be named S. oriscaviae. The oral cavity of guinea pigs is presumably a reservoir of S. oriscaviae. Some of the reported S. suis strains isolated from clinical specimens may be S. oriscaviae. IMPORTANCE We reported the discovery of a novel Streptococcus species, propose to be named Streptococcus oriscaviae, from the pus collected from a guinea pig bite wound in a healthy young patient. The bacterium was initially misidentified as S. suis/S. parasuis by biochemical tests, mass spectrometry. and housekeeping genes sequencing. Its novelty was confirmed by whole genome sequencing. Comparative genomic studies showed that S. oriscaviae belongs to the Suis group. S. oriscaviae sequences were detected in the oral cavities of the two guinea pigs owned by the patient, suggesting that the oral cavity of guinea pigs could be a reservoir of S. oriscaviae. Some of the reported S. suis strains may be S. oriscaviae. Further studies are warranted to refine our knowledge on this novel Streptococcus species.

Temporal Control of the WNT Signaling Pathway During Cardiac Differentiation Impacts Upon the Maturation State of Human Pluripotent Stem Cell Derived Cardiomyocytes.

Inefficient differentiation and insufficient maturation are barriers to the application of human pluripotent stem cell (hPSC)-derived cardiomyocytes (CMs) for research and therapy. Great strides have been made to the former, and multiple groups have reported cardiac differentiation protocol that can generate hPSC-CMs at high efficiency. Although many such protocols are based on the modulation of the WNT signaling pathway, they differ in their timing and in the WNT inhibitors used. Little is currently known about whether and how conditions of differentiation affect cardiac maturation. Here we adapted multiple cardiac differentiation protocols to improve cost-effectiveness and consistency, and compared the properties of the hPSC-CMs generated. Our results showed that the schedule of differentiation, but not the choice of WNT inhibitors, was a critical determinant not only of differentiation efficiency, which was expected, but also CM maturation. Among cultures with comparable purity, hPSC-CMs generated with different differentiation schedules vary in the expression of genes which are important for developmental maturation, and in their structural, metabolic, calcium transient and proliferative properties. In summary, we demonstrated that simple changes in the schedule of cardiac differentiation could promote maturation. To this end, we have optimized a cardiac differentiation protocol that can simultaneously achieve high differentiation efficiency and enhanced developmental maturation. Our findings would advance the production of hPSC-CMs for research and therapy.

Detecting structural variations with precise breakpoints using low-depth WGS data from a single oxford nanopore MinION flowcell.

Structural variation (SV) is a major cause of genetic disorders. In this paper, we show that low-depth (specifically, 4×) whole-genome sequencing using a single Oxford Nanopore MinION flow cell suffices to support sensitive detection of SV, particularly pathogenic SV for supporting clinical diagnosis. When using 4× ONT WGS data, existing SV calling software often fails to detect pathogenic SV, especially in the form of long deletion, terminal deletion, duplication, and unbalanced translocation. Our new SV calling software SENSV can achieve high sensitivity for all types of SV and a breakpoint precision typically ± 100 bp; both features are important for clinical concerns. The improvement achieved by SENSV stems from several new algorithms. We evaluated SENSV and other software using both real and simulated data. The former was based on 24 patient samples, each diagnosed with a genetic disorder. SENSV found the pathogenic SV in 22 out of 24 cases (all heterozygous, size from hundreds of kbp to a few Mbp), reporting breakpoints within 100 bp of the true answers. On the other hand, no existing software can detect the pathogenic SV in more than 10 out of 24 cases, even when the breakpoint requirement is relaxed to ± 2000 bp.

ECNano: A cost-effective workflow for target enrichment sequencing and accurate variant calling on 4800 clinically significant genes using a single MinION flowcell.

BACKGROUND: The application of long-read sequencing using the Oxford Nanopore Technologies (ONT) MinION sequencer is getting more diverse in the medical field. Having a high sequencing error of ONT and limited throughput from a single MinION flowcell, however, limits its applicability for accurate variant detection. Medical exome sequencing (MES) targets clinically significant exon regions, allowing rapid and comprehensive screening of pathogenic variants. By applying MES with MinION sequencing, the technology can achieve a more uniform capture of the target regions, shorter turnaround time, and lower sequencing cost per sample. METHOD: We introduced a cost-effective optimized workflow, ECNano, comprising a wet-lab protocol and bioinformatics analysis, for accurate variant detection at 4800 clinically important genes and regions using a single MinION flowcell. The ECNano wet-lab protocol was optimized to perform long-read target enrichment and ONT library preparation to stably generate high-quality MES data with adequate coverage. The subsequent variant-calling workflow, Clair-ensemble, adopted a fast RNN-based variant caller, Clair, and was optimized for target enrichment data. To evaluate its performance and practicality, ECNano was tested on both reference DNA samples and patient samples. RESULTS: ECNano achieved deep on-target depth of coverage (DoC) at average > 100× and > 98% uniformity using one MinION flowcell. For accurate ONT variant calling, the generated reads sufficiently covered 98.9% of pathogenic positions listed in ClinVar, with 98.96% having at least 30× DoC. ECNano obtained an average read length of 1000 bp. The long reads of ECNano also covered the adjacent splice sites well, with 98.5% of positions having ≥ 30× DoC. Clair-ensemble achieved > 99% recall and accuracy for SNV calling. The whole workflow from wet-lab protocol to variant detection was completed within three days. CONCLUSION: We presented ECNano, an out-of-the-box workflow comprising (1) a wet-lab protocol for ONT target enrichment sequencing and (2) a downstream variant detection workflow, Clair-ensemble. The workflow is cost-effective, with a short turnaround time for high accuracy variant calling in 4800 clinically significant genes and regions using a single MinION flowcell. The long-read exon captured data has potential for further development, promoting the application of long-read sequencing in personalized disease treatment and risk prediction.

Same-Cell Co-Occurrence of RAS Hotspot and BRAF V600E Mutations in Treatment-Naive Colorectal Cancer.

PURPOSE: Mitogen-activated protein kinase pathway-activating mutations occur in the majority of colorectal cancer (CRC) cases and show mutual exclusivity. We identified 47 epidermal growth factor receptor/BRAF inhibitor-naive CRC patients with dual RAS hotspot/BRAF V600E mutations (CRC-DD) from a cohort of 4,561 CRC patients with clinical next-generation sequencing results. We aimed to define the molecular phenotypes of the CRC-DD and to test if the dual RAS hotspot/BRAF V600E mutations coexist within the same cell. MATERIALS AND METHODS: We developed a single-cell genotyping method with a mutation detection rate of 96.3% and a genotype prediction accuracy of 92.1%. Mutations in the CRC-DD cohort were analyzed for clonality, allelic imbalance, copy number, and overall survival. RESULTS: Application of single-cell genotyping to four CRC-DD revealed the co-occurrence of both mutations in the following percentages of cells per case: NRAS G13D/KRAS G12C, 95%; KRAS G12D/NRAS G12V, 48%; BRAF V600E/KRAS G12D, 44%; and KRAS G12D/NRAS G13V, 14%, respectively. Allelic imbalance favoring the oncogenic allele was less frequent in CRC-DD (24 of 76, 31.5%, somatic mutations) compared with a curated cohort of CRC with a single-driver mutation (CRC-SD; 119 of 232 mutations, 51.3%; P = .013). Microsatellite instability-high status was enriched in CRC-DD compared with CRC-SD (23% v 11.4%, P = .028). Of the seven CRC-DD cases with multiregional sequencing, five retained both driver mutations throughout all sequenced tumor sites. Both CRC-DD cases with discordant multiregional sequencing were microsatellite instability-high. CONCLUSION: Our findings indicate that dual-driver mutations occur in a rare subset of CRC, often within the same tumor cells and across multiple tumor sites. Their presence and a lower rate of allelic imbalance may be related to dose-dependent signaling within the mitogen-activated protein kinase pathway.