ABSTRACT
Circular RNA (circRNA) is a class of RNA molecules that forms a closed loop with their 5' and 3' ends covalently bonded. CircRNAs are known to be more stable than linear RNAs, have distinct properties and functions, and are promising biomarkers. Existing methods for assembling circRNAs heavily rely on the annotated transcriptomes, hence exhibiting unsatisfactory accuracy without a high-quality transcriptome. We present TERRACE, a new algorithm for full-length assembly of circRNAs from paired-end total RNA-seq data. TERRACE uses the splice graph as the underlying data structure that organizes the splicing and coverage information. We transform the problem of assembling circRNAs into finding paths that "bridge" the three fragments in the splice graph induced by back-spliced reads. We adopt a definition for optimal bridging paths and a dynamic programming algorithm to calculate such optimal paths. TERRACE features an efficient algorithm to detect back-spliced reads missed by RNA-seq aligners, contributing to its much-improved sensitivity. It also incorporates a new machine-learning approach trained to assign a confidence score to each assembled circRNA, which is shown to be superior to using abundance for scoring. On both simulations and biological data sets, TERRACE consistently outperforms existing methods by a large margin in sensitivity while achieving better or comparable precision. In particular, when the annotations are not provided, TERRACE assembles 123%-413% more correct circRNAs than state-of-the-art methods. TERRACE presents a significant advance in assembling full-length circRNAs from RNA-seq data, and we expect it to be widely used in future research on circRNAs.
Subject(s)
Algorithms , RNA, Circular , RNA, Circular/genetics , Humans , Sequence Analysis, RNA/methods , RNA/genetics , RNA/metabolism , Machine Learning , Software , Computational Biology/methods , Transcriptome , RNA SplicingABSTRACT
BACKGROUND: Pelvic radiation plus sensitizing chemotherapy with a fluoropyrimidine (chemoradiotherapy) before surgery is standard care for locally advanced rectal cancer in North America. Whether neoadjuvant chemotherapy with fluorouracil, leucovorin, and oxaliplatin (FOLFOX) can be used in lieu of chemoradiotherapy is uncertain. METHODS: We conducted a multicenter, unblinded, noninferiority, randomized trial of neoadjuvant FOLFOX (with chemoradiotherapy given only if the primary tumor decreased in size by <20% or if FOLFOX was discontinued because of side effects) as compared with chemoradiotherapy. Adults with rectal cancer that had been clinically staged as T2 node-positive, T3 node-negative, or T3 node-positive who were candidates for sphincter-sparing surgery were eligible to participate. The primary end point was disease-free survival. Noninferiority would be claimed if the upper limit of the two-sided 90.2% confidence interval of the hazard ratio for disease recurrence or death did not exceed 1.29. Secondary end points included overall survival, local recurrence (in a time-to-event analysis), complete pathological resection, complete response, and toxic effects. RESULTS: From June 2012 through December 2018, a total of 1194 patients underwent randomization and 1128 started treatment; among those who started treatment, 585 were in the FOLFOX group and 543 in the chemoradiotherapy group. At a median follow-up of 58 months, FOLFOX was noninferior to chemoradiotherapy for disease-free survival (hazard ratio for disease recurrence or death, 0.92; 90.2% confidence interval [CI], 0.74 to 1.14; P = 0.005 for noninferiority). Five-year disease-free survival was 80.8% (95% CI, 77.9 to 83.7) in the FOLFOX group and 78.6% (95% CI, 75.4 to 81.8) in the chemoradiotherapy group. The groups were similar with respect to overall survival (hazard ratio for death, 1.04; 95% CI, 0.74 to 1.44) and local recurrence (hazard ratio, 1.18; 95% CI, 0.44 to 3.16). In the FOLFOX group, 53 patients (9.1%) received preoperative chemoradiotherapy and 8 (1.4%) received postoperative chemoradiotherapy. CONCLUSIONS: In patients with locally advanced rectal cancer who were eligible for sphincter-sparing surgery, preoperative FOLFOX was noninferior to preoperative chemoradiotherapy with respect to disease-free survival. (Funded by the National Cancer Institute; PROSPECT ClinicalTrials.gov number, NCT01515787.).
Subject(s)
Rectal Neoplasms , Adult , Humans , Anal Canal/surgery , Antineoplastic Combined Chemotherapy Protocols/adverse effects , Antineoplastic Combined Chemotherapy Protocols/therapeutic use , Chemoradiotherapy/adverse effects , Chemoradiotherapy/methods , Chemotherapy, Adjuvant , Disease-Free Survival , Fluorouracil/administration & dosage , Fluorouracil/adverse effects , Leucovorin/administration & dosage , Leucovorin/adverse effects , Neoadjuvant Therapy , Neoplasm Recurrence, Local/drug therapy , Neoplasm Staging , Organ Sparing Treatments , Oxaliplatin/administration & dosage , Oxaliplatin/adverse effects , Rectal Neoplasms/mortality , Rectal Neoplasms/pathology , Rectal Neoplasms/surgery , Preoperative Care , Preoperative PeriodABSTRACT
Accurate prediction of molecular properties is fundamental in drug discovery and development, providing crucial guidance for effective drug design. A critical factor in achieving accurate molecular property prediction lies in the appropriate representation of molecular structures. Presently, prevalent deep learning-based molecular representations rely on 2D structure information as the primary molecular representation, often overlooking essential three-dimensional (3D) conformational information due to the inherent limitations of 2D structures in conveying atomic spatial relationships. In this study, we propose employing the Gram matrix as a condensed representation of 3D molecular structures and for efficient pretraining objectives. Subsequently, we leverage this matrix to construct a novel molecular representation model, Pre-GTM, which inherently encapsulates 3D information. The model accurately predicts the 3D structure of a molecule by estimating the Gram matrix. Our findings demonstrate that Pre-GTM model outperforms the baseline Graphormer model and other pretrained models in the QM9 and MoleculeNet quantitative property prediction task. The integration of the Gram matrix as a condensed representation of 3D molecular structure, incorporated into the Pre-GTM model, opens up promising avenues for its potential application across various domains of molecular research, including drug design, materials science, and chemical engineering.
Subject(s)
Molecular Conformation , Models, Molecular , Drug Design , Deep Learning , Drug Discovery , AlgorithmsABSTRACT
BACKGROUND: Excitation-contraction (E-C) coupling processes become disrupted in heart failure (HF), resulting in abnormal Ca2+ homeostasis, maladaptive structural and transcriptional remodeling, and cardiac dysfunction. Junctophilin-2 (JP2) is an essential component of the E-C coupling apparatus but becomes site-specifically cleaved by calpain, leading to disruption of E-C coupling, plasmalemmal transverse tubule degeneration, abnormal Ca2+ homeostasis, and HF. However, it is not clear whether preventing site-specific calpain cleavage of JP2 is sufficient to protect the heart against stress-induced pathological cardiac remodeling in vivo. METHODS: Calpain-resistant JP2 knock-in mice (JP2CR) were generated by deleting the primary JP2 calpain cleavage site. Stress-dependent JP2 cleavage was assessed through in vitro cleavage assays and in isolated cardiomyocytes treated with 1 µmol/L isoproterenol by immunofluorescence. Cardiac outcomes were assessed in wild-type and JP2CR mice 5 weeks after transverse aortic constriction compared with sham surgery using echocardiography, histology, and RNA-sequencing methods. E-C coupling efficiency was measured by in situ confocal microscopy. E-C coupling proteins were evaluated by calpain assays and Western blotting. The effectiveness of adeno-associated virus gene therapy with JP2CR, JP2, or green fluorescent protein to slow HF progression was evaluated in mice with established cardiac dysfunction. RESULTS: JP2 proteolysis by calpain and in response to transverse aortic constriction and isoproterenol was blocked in JP2CR cardiomyocytes. JP2CR hearts are more resistant to pressure-overload stress, having significantly improved Ca2+ homeostasis and transverse tubule organization with significantly attenuated cardiac dysfunction, hypertrophy, lung edema, fibrosis, and gene expression changes relative to wild-type mice. JP2CR preserves the integrity of calpain-sensitive E-C coupling-related proteins, including ryanodine receptor 2, CaV1.2, and sarcoplasmic reticulum calcium ATPase 2a, by attenuating transverse aortic constriction-induced increases in calpain activity. Furthermore, JP2CR gene therapy after the onset of cardiac dysfunction was found to be effective at slowing the progression of HF and superior to wild-type JP2. CONCLUSIONS: The data presented here demonstrate that preserving JP2-dependent E-C coupling by prohibiting the site-specific calpain cleavage of JP2 offers multifaceted beneficial effects, conferring cardiac protection against stress-induced proteolysis, hypertrophy, and HF. Our data also indicate that specifically targeting the primary calpain cleavage site of JP2 by gene therapy approaches holds great therapeutic potential as a novel precision medicine for treating HF.
ABSTRACT
BACKGROUND: The docking protein IRS2 (insulin receptor substrate protein-2) is an important mediator of insulin signaling and may also regulate other signaling pathways. Murine hearts with cardiomyocyte-restricted deletion of IRS2 (cIRS2-KO) are more susceptible to pressure overload-induced cardiac dysfunction, implying a critical protective role of IRS2 in cardiac adaptation to stress through mechanisms that are not fully understood. There is limited evidence regarding the function of IRS2 beyond metabolic homeostasis regulation, particularly in the context of cardiac disease. METHODS: A retrospective analysis of an electronic medical record database was conducted to identify patients with IRS2 variants and assess their risk of cardiac arrhythmias. Arrhythmia susceptibility was examined in cIRS2-KO mice. The underlying mechanisms were investigated using confocal calcium imaging of ex vivo whole hearts and isolated cardiomyocytes to assess calcium handling, Western blotting to analyze the involved signaling pathways, and pharmacological and genetic interventions to rescue arrhythmias in cIRS2-KO mice. RESULTS: The retrospective analysis identified patients with IRS2 variants of uncertain significance with a potential association to an increased risk of cardiac arrhythmias compared with matched controls. cIRS2-KO hearts were found to be prone to catecholamine-sensitive ventricular tachycardia and reperfusion ventricular tachycardia. Confocal calcium imaging of ex vivo whole hearts and single isolated cardiomyocytes from cIRS2-KO hearts revealed decreased Ca²+ transient amplitudes, increased spontaneous Ca²+ sparks, and reduced sarcoplasmic reticulum Ca²+ content during sympathetic stress, indicating sarcoplasmic reticulum dysfunction. We identified that overactivation of the AKT1/NOS3 (nitric oxide synthase 3)/CaMKII (Ca2+/calmodulin-dependent protein kinase II)/RyR2 (type 2 ryanodine receptor) signaling pathway led to calcium mishandling and catecholamine-sensitive ventricular tachycardia in cIRS2-KO hearts. Pharmacological AKT inhibition or genetic stabilization of RyR2 rescued catecholamine-sensitive ventricular tachycardia in cIRS2-KO mice. CONCLUSIONS: Cardiac IRS2 inhibits sympathetic stress-induced AKT/NOS3/CaMKII/RyR2 overactivation and calcium-dependent arrhythmogenesis. This novel IRS2 signaling axis, essential for maintaining cardiac calcium homeostasis under stress, presents a promising target for developing new antiarrhythmic therapies.
ABSTRACT
MOTIVATION: High-throughput RNA sequencing has become indispensable for decoding gene activities, yet the challenge of reconstructing full-length transcripts persists. Traditional single-sample assemblers frequently produce fragmented transcripts, especially in single-cell RNA-seq data. While algorithms designed for assembling multiple samples exist, they encounter various limitations. RESULTS: We present Aletsch, a new assembler for multiple bulk or single-cell RNA-seq samples. Aletsch incorporates several algorithmic innovations, including a "bridging" system that can effectively integrate multiple samples to restore missed junctions in individual samples, and a new graph-decomposition algorithm that leverages "supporting" information across multiple samples to guide the decomposition of complex vertices. A standout feature of Aletsch is its application of a random forest model with 50 well-designed features for scoring transcripts. We demonstrate its robust adaptability across different chromosomes, datasets, and species. Our experiments, conducted on RNA-seq data from several protocols, firmly demonstrate Aletsch's significant outperformance over existing meta-assemblers. As an example, when measured with the partial area under the precision-recall curve (pAUC, constrained by precision), Aletsch surpasses the leading assemblers TransMeta by 22.9%-62.1% and PsiCLASS by 23.0%-175.5% on human datasets. AVAILABILITY AND IMPLEMENTATION: Aletsch is freely available at https://github.com/Shao-Group/aletsch. Scripts that reproduce the experimental results of this manuscript is available at https://github.com/Shao-Group/aletsch-test.
Subject(s)
Algorithms , RNA-Seq , Software , RNA-Seq/methods , Humans , High-Throughput Nucleotide Sequencing/methods , Sequence Analysis, RNA/methodsABSTRACT
MOTIVATION: Many tasks in sequence analysis ask to identify biologically related sequences in a large set. The edit distance, being a sensible model for both evolution and sequencing error, is widely used in these tasks as a measure. The resulting computational problem-to recognize all pairs of sequences within a small edit distance-turns out to be exceedingly difficult, since the edit distance is known to be notoriously expensive to compute and that all-versus-all comparison is simply not acceptable with millions or billions of sequences. Among many attempts, we recently proposed the locality-sensitive bucketing (LSB) functions to meet this challenge. Formally, a (d1,d2)-LSB function sends sequences into multiple buckets with the guarantee that pairs of sequences of edit distance at most d1 can be found within a same bucket while those of edit distance at least d2 do not share any. LSB functions generalize the locality-sensitive hashing (LSH) functions and admit favorable properties, with a notable highlight being that optimal LSB functions for certain (d1,d2) exist. LSB functions hold the potential of solving above problems optimally, but the existence of LSB functions for more general (d1,d2) remains unclear, let alone constructing them for practical use. RESULTS: In this work, we aim to utilize machine learning techniques to train LSB functions. With the development of a novel loss function and insights in the neural network structures that can potentially extend beyond this specific task, we obtained LSB functions that exhibit nearly perfect accuracy for certain (d1,d2), matching our theoretical results, and high accuracy for many others. Comparing to the state-of-the-art LSH method Order Min Hash, the trained LSB functions achieve a 2- to 5-fold improvement on the sensitivity of recognizing similar sequences. An experiment on analyzing erroneous cell barcode data is also included to demonstrate the application of the trained LSB functions. AVAILABILITY AND IMPLEMENTATION: The code for the training process and the structure of trained models are freely available at https://github.com/Shao-Group/lsb-learn.
Subject(s)
Algorithms , Computational Biology/methods , Machine LearningABSTRACT
Although evolutionary fates and expression patterns of duplicated genes have been extensively investigated, how duplicated genes co-regulate a biological process in polyploids remains largely unknown. Here, we identified two gsdf (gonadal somatic cell-derived factor) homeologous genes (gsdf-A and gsdf-B) in hexaploid gibel carp (Carassius gibelio), wherein each homeolog contained three highly conserved alleles. Interestingly, gsdf-A and gsdf-B transcription were mainly activated by dmrt1-A (dsx- and mab-3-related transcription factor 1) and dmrt1-B, respectively. Loss of either gsdf-A or gsdf-B alone resulted in partial male-to-female sex reversal and loss of both caused complete sex reversal, which could be rescued by a nonsteroidal aromatase inhibitor. Compensatory expression of gsdf-A and gsdf-B was observed in gsdf-B and gsdf-A mutants, respectively. Subsequently, we determined that in tissue culture cells, Gsdf-A and Gsdf-B both interacted with Ncoa5 (nuclear receptor coactivator 5) and blocked Ncoa5 interaction with Rora (retinoic acid-related orphan receptor-alpha) to repress Rora/Ncoa5-induced activation of cyp19a1a (cytochrome P450, family 19, subfamily A, polypeptide 1a). These findings illustrate that Gsdf-A and Gsdf-B can regulate male differentiation by inhibiting cyp19a1a transcription in hexaploid gibel carp and also reveal that Gsdf-A and Gsdf-B can interact with Ncoa5 to suppress cyp19a1a transcription in vitro. This study provides a typical case of cooperative mechanism of duplicated genes in polyploids and also sheds light on the conserved evolution of sex differentiation.
Subject(s)
Gonads , Sex Differentiation , Animals , Cell Differentiation/genetics , Female , Fish Proteins/genetics , Fishes/genetics , Gene Expression Regulation, Developmental , Gonads/metabolism , Male , Polyploidy , Sex Differentiation/genetics , Transcription Factors/genetics , Transcription Factors/metabolismABSTRACT
Streptococcus pneumoniae (pneumococcus) causes a wide range of important human infectious diseases, including pneumonia, pneumonia-derived sepsis, otitis media, and meningitis. Pneumococcus produces numerous secreted proteins that are critical for normal physiology and pathogenesis. The membrane targeting and translocation of these secreted proteins are partly mediated by the signal recognition particle (SRP) complex, which consists of 4.5S small cytoplasmic RNA (ScRNA), and the Ffh, and FtsY proteins. Here, we report that pneumococcal ∆scRNA, ∆ffh, and ∆ftsY mutants were significantly impaired in competence induction, competence pili production, exogenous DNA uptake, and genetic transformation. Also, the ∆scRNA mutant was significantly attenuated in the mouse models of bacteremia and pneumonia. Interestingly, unlike the ∆scRNA, both ∆ffh and ∆ftsY mutants had growth defects on Todd-Hewitt Agar, which were alleviated by the provision of free amino acids or serum. Differences in nutritional requirements between ∆ffh and ∆ftsY vs ∆scRNA suggest that Ffh and FtsY may be partially functional in the absence of ScRNA. Finally, the insertase YidC2, which could functionally rescue some SRP mutations in other streptococcal species, was not essential for pneumococcal genetic transformation. Collectively, these results indicate that ScRNA is crucial for the successful development of genetic competence and virulence in pneumococcus. IMPORTANCE: Streptococcus pneumoniae (pneumococcus) causes multiple important infectious diseases in humans. The signal recognition particle (SRP) complex, which comprised 4.5S small cytoplasmic RNA (ScRNA), and the Ffh and FtsY proteins, mediates membrane targeting and translocation of secreted proteins in all organisms. However, the role of SRP and ScRNA has not been characterized during the induction of the competence system for genetic transformation and virulence in pneumococcus. By using a combination of genetic, biochemical, proteomic, and imaging approaches, we demonstrated that the SRP complex plays a significant role in membrane targeting of competence system-regulated effectors important for genetic transformation, virulence during bacteremia and pneumonia infections, and nutritional acquisition.
Subject(s)
Bacterial Proteins , Streptococcus pneumoniae , Streptococcus pneumoniae/genetics , Streptococcus pneumoniae/pathogenicity , Streptococcus pneumoniae/metabolism , Mice , Bacterial Proteins/genetics , Bacterial Proteins/metabolism , Virulence , Animals , Signal Recognition Particle/genetics , Signal Recognition Particle/metabolism , Pneumococcal Infections/microbiology , Gene Expression Regulation, Bacterial , RNA, Bacterial/genetics , RNA, Bacterial/metabolism , DNA Transformation Competence , Bacteremia/microbiologyABSTRACT
Surface modification could enhance the cell internalization efficiency of nanovehicles for targeted gene or drug delivery. However, the influence of surface modification parameters, including recognition manners, valences, and patterns, is often clouded, especially for the endocytosis of DNA nanostructures in customized shapes. Focusing on an icosahedral DNA framework, we systematically programmed three distinct types of ligands with diverse valence and spatial distribution on their outer surface to study the internalization efficiency, endocytic pathways, and postinternalization fate. The comparison in different aspects of parameters deepens our understanding of the intricate relationship between surface modification and cell entry behavior, offering insights crucial for designing and optimizing DNA framework nanostructures for potent cell-targeted purposes.
Subject(s)
DNA , Endocytosis , Surface Properties , DNA/chemistry , Humans , Nanostructures/chemistry , Ligands , HeLa CellsABSTRACT
BACKGROUND: Dihydropyrimidine dehydrogenase (DPD) deficiency is the main known cause of life-threatening fluoropyrimidine (FP)-induced toxicities. We conducted a meta-analysis on individual patient data to assess the contribution of deleterious DPYD variants *2A/D949V/*13/HapB3 (recommended by EMA) and clinical factors, for predicting G4-5 toxicity. METHODS: Study eligibility criteria included recruitment of Caucasian patients without DPD-based FP-dose adjustment. Main endpoint was 12-week haematological or digestive G4-5 toxicity. The value of DPYD variants *2A/p.D949V/*13 merged, HapB3, and MIR27A rs895819 was evaluated using multivariable logistic models (AUC). RESULTS: Among 25 eligible studies, complete clinical variables and primary endpoint were available in 15 studies (8733 patients). Twelve-week G4-5 toxicity prevalence was 7.3% (641 events). The clinical model included age, sex, body mass index, schedule of FP-administration, concomitant anticancer drugs. Adding *2A/p.D949V/*13 variants (at least one allele, prevalence 2.2%, OR 9.5 [95%CI 6.7-13.5]) significantly improved the model (p < 0.0001). The addition of HapB3 (prevalence 4.0%, 98.6% heterozygous), in spite of significant association with toxicity (OR 1.8 [95%CI 1.2-2.7]), did not improve the model. MIR27A rs895819 was not associated with toxicity, irrespective of DPYD variants. CONCLUSIONS: FUSAFE meta-analysis highlights the major relevance of DPYD *2A/p.D949V/*13 combined with clinical variables to identify patients at risk of very severe FP-related toxicity.
Subject(s)
Antineoplastic Agents , Dihydropyrimidine Dehydrogenase Deficiency , Humans , Fluorouracil/adverse effects , Dihydrouracil Dehydrogenase (NADP)/genetics , Heterozygote , Genotype , Capecitabine/adverse effectsABSTRACT
Exosomes have emerged as a revolutionary tool for liquid biopsy (LB), as they carry specific cargo from cells. Profiling the metabolites of exosomes is crucial for cancer diagnosis and biomarker discovery. Herein, we propose a versatile platform for exosomal metabolite assay of endometrial cancer (EC). The platform is based on a nanostructured composite material comprising gold nanoparticle-coated magnetic COF with aptamer modification (Fe3O4@COF@Au-Apt). The unique design and novel synthesis strategy of Fe3O4@COF@Au-Apt provide the material with a large specific surface area, enabling the efficient and specific isolation of exosomes. The exosomes captured Fe3O4@COF@Au-Apt can be directly used as the laser desorption/ionization mass spectrometry (LDI-MS) matrix for rapid exosomal metabolic patterns. By integrating these functionalities into a single platform, the analytical process is simplified, eliminating the need for additional elution steps and minimizing potential sample loss, resulting in large-scale exosomal metabolic fingerprints. Combining with machine learning algorithms on the metabolic patterns, accurate discrimination between endometrial patients (EGs) and benign controls (CGs) was achieved, and the area under the receiver operating characteristic curve of the blind test cohort was 0.924. Confusion matrix analysis of important metabolic fingerprint features further demonstrates the high accuracy of the proposed approach toward EC diagnosis, with an overall accuracy of 94.1%. Moreover, four metabolites, namely, hydroxychalcone, l-acetylcarnitine, elaidic acid, and glutathione, have been identified as potential biomarkers of EC. These results highlight the great value of the integrated exosome metabolic fingerprint platform in facilitating low-cost and high-throughput characterization of exosomal metabolites for cancer diagnosis and biomarker discovery.
ABSTRACT
Stimulus-responsive mode is highly desirable for improving the precise monitoring and physiological efficacy of endogenous biomarkers (EB). However, its integrated application for visual detection and therapy is limited by inappropriate use of responsive triggers and poor delivery of EB signal-transducing agents, which remain challenging in simultaneous monitoring and noninvasive therapy of EB and EB-mediated pathological events. Target microRNA (miRNA) as controllable reaction triggers and DNAzyme as signal-transducing agent are proposed to develop target-stimulated multifunctional nanocabinets (MFNCs) for the visual tracking of both miRNA and miRNA-mediated anticancer events. The MFNCs, equipped with a target-discriminating sequence-incorporated DNAzyme motif, can specifically release therapeutic molecules through target-triggered conformational switches, accompanied by transduction signal output. Target detection and molecule release performance are recorded in parallel via reverse dual-signal feedback at the single-molecule level. In addition, the intrinsic thermal-replenishing of the MFNCs leads to tumor ablation without invasive exogenous aids. The system achieves visual target quantification, anticancer molecule real-time tracking, and tumor suppression in vivo and in vitro. This work proposes a new paradigm for precise visual exploration of EB or EB-mediated bio-events and provides a demonstration of efficacious all-in-one detection and therapy based on the target-triggered multifunctional nanosystem.
Subject(s)
DNA, Catalytic , MicroRNAs , Neoplasms , Humans , Feedback , MicroRNAs/genetics , Neoplasms/drug therapyABSTRACT
MOTIVATION: Modern methods for computation-intensive tasks in sequence analysis (e.g. read mapping, sequence alignment, genome assembly, etc.) often first transform each sequence into a list of short, regular-length seeds so that compact data structures and efficient algorithms can be employed to handle the ever-growing large-scale data. Seeding methods using kmers (substrings of length k) have gained tremendous success in processing sequencing data with low mutation/error rates. However, they are much less effective for sequencing data with high error rates as kmers cannot tolerate errors. RESULTS: We propose SubseqHash, a strategy that uses subsequences, rather than substrings, as seeds. Formally, SubseqHash maps a string of length n to its smallest subsequence of length k, k < n, according to a given order overall length-k strings. Finding the smallest subsequence of a string by enumeration is impractical as the number of subsequences grows exponentially. To overcome this barrier, we propose a novel algorithmic framework that consists of a specifically designed order (termed ABC order) and an algorithm that computes the minimized subsequence under an ABC order in polynomial time. We first show that the ABC order exhibits the desired property and the probability of hash collision using the ABC order is close to the Jaccard index. We then show that SubseqHash overwhelmingly outperforms the substring-based seeding methods in producing high-quality seed-matches for three critical applications: read mapping, sequence alignment, and overlap detection. SubseqHash presents a major algorithmic breakthrough for tackling the high error rates and we expect it to be widely adapted for long-reads analysis. AVAILABILITY AND IMPLEMENTATION: SubseqHash is freely available at https://github.com/Shao-Group/subseqhash.
Subject(s)
Algorithms , Mutation Rate , Probability , Sequence AlignmentABSTRACT
BACKGROUND: Postoperative adverse events (AEs) in patients with borderline resectable pancreatic ductal adenocarcinoma (BR-PC) treated with neoadjuvant therapy and pancreatectomy in the national cooperative group setting have not been previously characterized. We conducted a preplanned secondary analysis of patients enrolled on the Alliance A021501 clinical trial to quantify perioperative AE rates. METHODS: The A021501 phase 2 trial randomized patients with BR-PC to receive 8 doses of mFOLFIRINOX (Arm 1) or 7 doses of mFOLFIRINOX and hypofractionated radiotherapy (Arm 2), followed by pancreatectomy (December 31, 2016 to May 31, 2019). Adverse events were assessed 90 days after pancreatectomy. RESULTS: Of 126 enrolled patients, 51 (40%) underwent pancreatectomy (n = 32, Arm 1; n = 19, Arm 2) at 28 institutions. Five (10%) patients required reoperation within 90 days; 56% of patients (n = 27/48) experienced at least one grade 3 or higher AE (50% vs. 67%, p = 0.37). Ninety-day mortality was 2.0%. Readmission was less frequent in Arm 1 (16% vs. 42%, p = 0.05), but there were no differences between study arms in rates of reoperation (13% vs. 5%), pancreatic fistula or intra-abdominal abscess requiring drainage (9% vs. 16%), or wound infection (6% vs. 16%). Pancreatic fistula or intra-abdominal abscess requiring drainage was associated with receipt of adjuvant therapy (p = 0.012). No difference in overall survival was observed based on occurrence of postoperative AEs (hazard ratio = 1.1; 95% confidence interval 0.5-2.6). CONCLUSIONS: In this multicenter study, rates of postoperative AEs were consistent with those previously reported. Multimodality trials of preoperative therapy for BR-PC may be performed in the cooperative group setting with careful quality assurance and safety monitoring. TRIAL REGISTRATION: Clinicaltrials.gov identifier: NCT02839343.
Subject(s)
Antineoplastic Combined Chemotherapy Protocols , Carcinoma, Pancreatic Ductal , Fluorouracil , Neoadjuvant Therapy , Oxaliplatin , Pancreatectomy , Pancreatic Neoplasms , Postoperative Complications , Humans , Pancreatic Neoplasms/pathology , Pancreatic Neoplasms/therapy , Pancreatic Neoplasms/surgery , Female , Male , Middle Aged , Aged , Postoperative Complications/etiology , Antineoplastic Combined Chemotherapy Protocols/therapeutic use , Carcinoma, Pancreatic Ductal/surgery , Carcinoma, Pancreatic Ductal/pathology , Carcinoma, Pancreatic Ductal/therapy , Oxaliplatin/administration & dosage , Survival Rate , Fluorouracil/administration & dosage , Irinotecan/administration & dosage , Follow-Up Studies , Leucovorin/administration & dosage , Prognosis , Adult , Combined Modality TherapyABSTRACT
Observational studies and stand-alone trials indicate that patients with follicular lymphoma (FL) who experience disease progression within 24 months of front-line chemoimmunotherapy (POD24), have poor outcomes. We performed a pooled analysis of 13 randomized clinical trials of patients with FL in the pre- and postrituximab eras to identify clinical factors that predict POD24. Logistic regression models evaluated the association between clinical factors and POD24. Cox regression evaluated the association between POD24 as a time-dependent factor and subsequent overall survival (OS). A landmark analysis evaluated the association of POD24 with OS for the subset of patients who were alive at 24 months after trial registration. Patients without progression at 24 months at baseline had favorable performance status (PS), limited-stage (I/II) disease, low-risk FL International Prognostic Index (FLIPI) score, normal baseline hemoglobin, and normal baseline ß2 microglobulin (B2M) level. In a multivariable logistic regression model, male sex (odds ratio [OR], 1.30), PS ≥2 (OR, 1.63), B2M (≥3 mg/L; OR, 1.43), and high-risk FLIPI score (3-5; OR, 3.14) were associated with increased risk of progression before 24 months. In the time-dependent Cox model and the 24-month landmark analysis, POD24 was associated with poor subsequent OS (hazard ratio, 4.85 and 3.06, respectively). This is the largest pooled analysis of clinical trials data validating POD24 as a robust indicator of poor FL survival and identified clinical predictors of early death and progression that can aid in building comprehensive prognostic models incorporating clinical and molecular predictors of POD24.
Subject(s)
Lymphoma, Follicular , Antineoplastic Combined Chemotherapy Protocols , Disease Progression , Humans , Immunotherapy , Lymphoma, Follicular/diagnosis , Lymphoma, Follicular/drug therapy , Male , Prognosis , Risk FactorsABSTRACT
BACKGROUND: For patients with liver-confined metastatic colorectal cancer (mCRC), local therapy of isolated metastases has been associated with long-term progression-free and overall survival (OS). However, for patients with more advanced mCRC, including those with extrahepatic disease, the efficacy of local therapy is less clear although increasingly being used in clinical practice. Prospective studies to clarify the role of metastatic-directed therapies in patients with mCRC are needed. METHODS: The Evaluating Radiation, Ablation, and Surgery (ERASur) A022101/NRG-GI009 trial is a randomized, National Cancer Institute-sponsored phase III study evaluating if the addition of metastatic-directed therapy to standard of care systemic therapy improves OS in patients with newly diagnosed limited mCRC. Eligible patients require a pathologic diagnosis of CRC, have BRAF wild-type and microsatellite stable disease, and have 4 or fewer sites of metastatic disease identified on baseline imaging. Liver-only metastatic disease is not permitted. All metastatic lesions must be amenable to total ablative therapy (TAT), which includes surgical resection, microwave ablation, and/or stereotactic ablative body radiotherapy (SABR) with SABR required for at least one lesion. Patients without overt disease progression after 16-26 weeks of first-line systemic therapy will be randomized 1:1 to continuation of systemic therapy with or without TAT. The trial activated through the Cancer Trials Support Unit on January 10, 2023. The primary endpoint is OS. Secondary endpoints include event-free survival, adverse events profile, and time to local recurrence with exploratory biomarker analyses. This study requires a total of 346 evaluable patients to provide 80% power with a one-sided alpha of 0.05 to detect an improvement in OS from a median of 26 months in the control arm to 37 months in the experimental arm with a hazard ratio of 0.7. The trial uses a group sequential design with two interim analyses for futility. DISCUSSION: The ERASur trial employs a pragmatic interventional design to test the efficacy and safety of adding multimodality TAT to standard of care systemic therapy in patients with limited mCRC. TRIAL REGISTRATION: ClinicalTrials.gov: NCT05673148, registered December 21, 2022.
Subject(s)
Colonic Neoplasms , Liver Neoplasms , Radiosurgery , Rectal Neoplasms , Humans , Prospective Studies , Radiosurgery/methods , Liver Neoplasms/therapyABSTRACT
BACKGROUND: Recent data have demonstrated that in locally advanced rectal cancer (LARC), a total neoadjuvant therapy (TNT) approach improves compliance with chemotherapy and increases rates of tumor response compared to neoadjuvant chemoradiation (CRT) alone. They further indicate that the optimal sequencing of TNT involves consolidation (rather than induction) chemotherapy to optimize complete response rates. Data, largely from retrospective studies, have also shown that patients with clinical complete response (cCR) after TNT may be managed safely with the watch and wait approach (WW) instead of preemptive total mesorectal resection (TME). However, the optimal consolidation chemotherapy regimen to achieve cCR has not been established, and a randomized clinical trial has not robustly evaluated cCR as a primary endpoint. Collaborating with a multidisciplinary oncology team and patient groups, we designed this NCI-sponsored study of chemotherapy intensification to address these issues and to drive up cCR rates, to provide opportunity for organ preservation, improve quality of life for patients and improve survival outcomes. METHODS: In this NCI-sponsored multi-group randomized, seamless phase II/III trial (1:1), up to 760 patients with LARC, T4N0, any T with node positive disease (any T, N +) or T3N0 requiring abdominoperineal resection or coloanal anastomosis and distal margin within 12 cm of anal verge will be enrolled. Stratification factors include tumor stage (T4 vs T1-3), nodal stage (N + vs N0) and distance from anal verge (0-4; 4-8; 8-12 cm). Patients will be randomized to receive neoadjuvant long-course chemoradiation (LCRT) followed by consolidation doublet (mFOLFOX6 or CAPOX) or triplet chemotherapy (mFOLFIRINOX) for 3-4 months. LCRT in both arms involves 4500 cGy in 25 fractions over 5 weeks + 900 cGy boost in 5 fractions with a fluoropyrimidine (capecitabine preferred). Patients will undergo assessment 8-12 (± 4) weeks post-TNT completion. The primary endpoint for the phase II portion will compare cCR between treatment arms. A total number of 312 evaluable patients (156 per arm) will provide statistical power of 90.5% to detect a 17% increase in cCR rate, at a one-sided alpha = 0.048. The primary endpoint for the phase III portion will compare disease-free survival (DFS) between treatment arms. A total of 285 DFS events will provide 85% power to detect an effect size of hazard ratio 0.70 at a one-sided alpha of 0.025, requiring enrollment of 760 patients (380 per arm). Secondary objectives include time-to event outcomes (overall survival, organ preservation time and time to distant metastasis) and adverse event rates. Biospecimens including archival tumor tissue, plasma and buffy coat, and serial rectal MRIs will be collected for exploratory correlative research. This study, activated in late 2022, is open across the NCTN and had accrued 330 patients as of May 2024. Study support: U10CA180821, U10CA180882, U24 CA196171; https://acknowledgments.alliancefound.org . DISCUSSION: Building on data from modern day rectal cancer trials and patient input from national advocacy groups, we have designed The Janus Rectal Cancer Trial studying chemotherapy intensification via a consolidation chemotherapy approach with the intent to enhance cCR and DFS rates, increase organ preservation rates, and improve quality of life for patients with rectal cancer. TRIAL REGISTRATION: Clinicaltrials.gov ID: NCT05610163; Support includes U10CA180868 (NRG) and U10CA180888 (SWOG).
Subject(s)
Antineoplastic Combined Chemotherapy Protocols , Fluorouracil , Neoadjuvant Therapy , Rectal Neoplasms , Humans , Rectal Neoplasms/therapy , Rectal Neoplasms/pathology , Rectal Neoplasms/mortality , Rectal Neoplasms/drug therapy , Antineoplastic Combined Chemotherapy Protocols/therapeutic use , Neoadjuvant Therapy/methods , Fluorouracil/administration & dosage , Fluorouracil/therapeutic use , Male , Female , Disease-Free Survival , Leucovorin/administration & dosage , Leucovorin/therapeutic use , Oxaliplatin/administration & dosage , Oxaliplatin/therapeutic use , Capecitabine/administration & dosage , Capecitabine/therapeutic use , Irinotecan/administration & dosage , Irinotecan/therapeutic use , Middle Aged , Treatment Outcome , Quality of Life , Neoplasm Staging , Organoplatinum CompoundsABSTRACT
Fuchs endothelial corneal dystrophy (FECD) is the leading cause of endothelial keratoplasty without efficacious drug treatment. Recent studies have emphasized the involvement of epigenetic regulation in FECD development. Long non-coding RNAs (lncRNAs) are recognized as crucial epigenetic regulators in diverse cellular processes and ocular diseases. In this study, we revealed the expression patterns of lncRNAs using high-throughput sequencing technology in FECD mouse model, and identified 979 significantly dysregulated lncRNAs. By comparing the data from FECD human cell model, we obtained a series of homologous lncRNAs with similar expression patterns, and revealed that these homologous lncRNAs were enriched in FECD related biological functions, with apoptosis (mmu04210) showing the highest enrichment score. In addition, we investigated the role of lncRNA zinc finger antisense 1 (ZFAS1) in apoptotic process. This study would broaden our understanding of epigenetic regulation in FECD development, and provide potential anti-apoptotic targets for FECD therapy.
Subject(s)
Fuchs' Endothelial Dystrophy , RNA, Long Noncoding , Animals , Humans , Mice , Endothelium, Corneal/metabolism , Epigenesis, Genetic , Fuchs' Endothelial Dystrophy/genetics , Fuchs' Endothelial Dystrophy/metabolism , RNA, Long Noncoding/genetics , Zinc/metabolismABSTRACT
BACKGROUND: Transcriptional remodeling is known to contribute to heart failure (HF). Targeting stress-dependent gene expression mechanisms may represent a clinically relevant gene therapy option. We recently uncovered a salutary mechanism in the heart whereby JP2 (junctophilin-2), an essential component of the excitation-contraction coupling apparatus, is site-specifically cleaved and releases an N-terminal fragment (JP2NT [N-terminal fragment of JP2]) that translocates into the nucleus and functions as a transcriptional repressor of HF-related genes. This study aims to determine whether JP2NT can be leveraged by gene therapy techniques for attenuating HF progression in a preclinical pressure overload model. METHODS: We intraventricularly injected adeno-associated virus (AAV) (2/9) vectors expressing eGFP (enhanced green fluorescent protein), JP2NT, or DNA-binding deficient JP2NT (JP2NTΔbNLS/ARR) into neonatal mice and induced cardiac stress by transaortic constriction (TAC) 9 weeks later. We also treated mice with established moderate HF from TAC stress with either AAV-JP2NT or AAV-eGFP. RNA-sequencing analysis was used to reveal changes in hypertrophic and HF-related gene transcription by JP2NT gene therapy after TAC. Echocardiography, confocal imaging, and histology were performed to evaluate heart function and pathological myocardial remodeling following stress. RESULTS: Mice preinjected with AAV-JP2NT exhibited ameliorated cardiac remodeling following TAC. The JP2NT DNA-binding domain is required for cardioprotection as its deletion within the AAV-JP2NT vector prevented improvement in TAC-induced cardiac dysfunction. Functional and histological data suggest that JP2NT gene therapy after the onset of cardiac dysfunction is effective at slowing the progression of HF. RNA-sequencing analysis further revealed a broad reversal of hypertrophic and HF-related gene transcription by JP2NT overexpression after TAC. CONCLUSIONS: Our prevention- and intervention-based approaches here demonstrated that AAV-mediated delivery of JP2NT into the myocardium can attenuate stress-induced transcriptional remodeling and the development of HF when administered either before or after cardiac stress initiation. Our data indicate that JP2NT gene therapy holds great potential as a novel therapeutic for treating hypertrophy and HF.