Direct Sequencing Research Articles

The dynamics of N6-methyladenine RNA modification in resistant and susceptible rice varieties responding to rice stem borer damage.

N6-methyladenosine (m6A) is the most prevalent modification in cellular RNA which orchestrates diverse physiological and pathological processes during stress response. However, the differential m6A modifications that cope with herbivore stress in resistant and susceptible crop varieties remain unclear. Here, we found that rice stem borer (RSB) larvae grew better on indica rice (e.g., MH63, IR64, Nanjing 11) than on japonica rice varieties (e.g., Nipponbare, Zhonghua 11, Xiushui 11). Then, transcriptome-wide m6A profiling of representative resistant (Nipponbare) and susceptible (MH63) rice varieties were performed using a nanopore direct RNA sequencing approach, to reveal variety-specific m6A modifications against RSB. Upon RSB infestation, m6A methylation occurred in actively expressed genes in Nipponbare and MH63, but the number of methylation sites decreased across rice chromosomes. Integrative analysis showed that m6A methylation levels were closely associated with transcriptional regulation. Genes involved in herbivorous resistance related to mitogen-activated protein kinase, jasmonic acid (JA), and terpenoid biosynthesis pathways, as well as JA-mediated trypsin protease inhibitors, were heavily methylated by m6A, and their expression was more pronounced in RSB-infested Nipponbare than in RSB-infested MH63, which may have contributed to RSB resistance in Nipponbare. Therefore, dynamics of m6A modifications act as the main regulatory strategy for expression of genes involved in plant-insect interactions, which is attributed to differential responses of resistant and susceptible rice varieties to RSB infestation. These findings could contribute to developing molecular breeding strategies for controlling herbivorous pests.

RNA m6A detection using raw current signals and basecalling errors from Nanopore direct RNA sequencing reads.

Nanopore direct RNA sequencing (DRS) enables the detection of RNA N6-methyladenosine (m6A) without extra laboratory techniques. A number of supervised or comparative approaches have been developed to identify m6A from Nanopore DRS reads. However, existing methods typically utilize either statistical features of the current signals or basecalling-error features, ignoring the richer information of the raw signals of DRS reads. Here, we propose RedNano, a deep-learning method designed to detect m6A from Nanopore DRS reads by utilizing both raw signals and basecalling errors. RedNano processes the raw-signal feature and basecalling-error feature through residual networks. We validated the effectiveness of RedNano using synthesized, Arabidopsis, and human DRS data. The results demonstrate that RedNano surpasses existing methods by achieving higher area under the ROC curve (AUC) and area under the precision-recall curve (AUPRs) in all three datasets. Furthermore, RedNano performs better in cross-species validation, demonstrating its robustness. Additionally, when detecting m6A from an independent dataset of Populus trichocarpa, RedNano achieves the highest AUC and AUPR, which are 3.8%-9.9% and 5.5%-13.8% higher than other methods, respectively. The source code of RedNano is freely available at https://github.com/Derryxu/RedNano.

PT-seq: A method for metagenomic analysis of phosphorothioate epigenetics in complex microbial communities.

Among dozens of known epigenetic marks, naturally occurring phosphorothioate (PT) DNA modifications are unique in replacing a non-bridging phosphate oxygen with redox-active sulfur and function in prokaryotic restriction-modification and transcriptional regulation. Interest in PTs has grown due to the widespread distribution of the dnd, ssp, and brx genes among bacteria and archaea, as well as the discovery of PTs in 5-10% of gut microbes. Efforts to map PTs in complex microbiomes using existing next-generation and direct sequencing technologies have failed due to poor sensitivity. Here we developed PT-seq as a high-sensitivity method to quantitatively map PTs across genomes and metagenomically identify PT-containing microbes in complex genomic mixtures. Like other methods for mapping PTs in individual genomes, PT-seq exploits targeted DNA strand cleavage at PTs by iodine, followed by sequencing library construction using ligation or template switching approaches. However, PT-specific sequencing reads are dramatically increased by adding steps to heat denature the DNA, block pre-existing 3'-ends, fragment DNA after T-tailing, and enrich iodine-induced breaks using biotin-labeling and streptavidin beads capture. Iterative optimization of the sensitivity and specificity of PT-seq is demonstrated with individual bacteria and human fecal DNA.

Detection and Quantification of 5moU RNA Modification from Direct RNA Sequencing Data.

Chemically modified therapeutic mRNAs have gained momentum recently. In addition to commonly used modifications (e.g., pseudouridine), 5moU is considered a promising substitution for uridine in therapeutic mRNAs. Accurate identification of 5-methoxyuridine (5moU) would be crucial for the study and quality control of relevant in vitro-transcribed (IVT) mRNAs. However, current methods exhibit deficiencies in providing quantitative methodologies for detecting such modification. Utilizing the capabilities of Oxford nanopore direct RNA sequencing, in this study, we present NanoML-5moU, a machine-learning framework designed specifically for the read-level detection and quantification of 5moU modification for IVT data. Nanopore direct RNA sequencing data from both 5moU-modified and unmodified control samples were collected. Subsequently, a comprehensive analysis and modeling of signal event characteristics (mean, median current intensities, standard deviations, and dwell times) were performed. Furthermore, classical machine learning algorithms, notably the Support Vector Machine (SVM), Random Forest (RF), and XGBoost were employed to discern 5moU modifications within NNUNN (where N represents A, C, U, or G) 5-mers. Notably, the signal event attributes pertaining to each constituent base of the NNUNN 5-mers, in conjunction with the utilization of the XGBoost algorithm, exhibited remarkable performance levels (with a maximum AUROC of 0.9567 in the "AGTTC" reference 5-mer dataset and a minimum AUROC of 0.8113 in the "TGTGC" reference 5-mer dataset). This accomplishment markedly exceeded the efficacy of the prevailing background error comparison model (ELIGOs AUC 0.751 for site-level prediction). The model's performance was further validated through a series of curated datasets, which featured customized modification ratios designed to emulate broader data patterns, demonstrating its general applicability in quality control of IVT mRNA vaccines. The NanoML-5moU framework is publicly available on GitHub (https://github.com/JiayiLi21/NanoML-5moU). NanoML-5moU enables accurate read-level profiling of 5moU modification with nanopore direct RNA-sequencing, which is a powerful tool specialized in unveiling signal patterns in in vitro-transcribed (IVT) mRNAs.

A machine learning algorithm based on multi-omics biomarkers for the detection of tumor microsatellite instability.

1554 Background: Microsatellite instability (MSI) occurs via defects in mismatch repair (MMR) and is characterized by high tumor mutational burden (TMB) and a robust inflammatory response. Accurate MSI-high classification of immune-evasive tumors is vital due to the effectiveness of immune checkpoint inhibitor (ICI) therapy, such as pembrolizumab, which is FDA approved for all MSI-high tumors irrespective of anatomic origin. Detection methods of MSI-status vary widely, including immunohistochemistry (IHC), polymerase chain reaction (PCR), and next-generation sequencing (NGS). We developed a machine learning (ML) model to predict MSI status in patients with solid tumors using comprehensive genomic and immune profiling (CGIP), independent of direct sequencing data from microsatellite sites. Methods: We analyzed samples from 1,838 patients with colorectal cancer (CRC) by CGIP, which included DNA panel testing (523 genes) for pathogenic single nucleotide variants (pSNV) and determination of TMB (mut/Mb), and RNA sequencing (397 genes) for gene expression (GEX). The Boruta algorithm was used to select key genomic and GEX changes associated with MSI status. A distributed gradient boosting algorithm created a predictive model that was trained on 70% of the cohort. The model was tested on the remaining CRC cohort (Test), and in the PanCancer Atlas CRCs (TCGA), endometrial adenocarcinomas (EMCA), and CRC with indeterminate MSI testing by CGIP. Performance of the trained model was assessed using sensitivity, specificity, and positive (PPV) and negative (NPV) predictive value. Results: Feature selection identified 79 genes with pSNVs, 63 GEX changes, and TMB as informative for MSI prediction. The model showed high predictive accuracy in differentiating MSI-H tumors. Of the 39 cases that failed MSI component of CGIP (Indeterminate), 17 had MMR IHC results available, of which 1 demonstrated loss of MLH1/PMS2 and the other 16 showed intact (normal) expression. Testing of this cohort identified 3 cases as MSI-high, at least one of which was confirmed by IHC. No cases with intact MMR protein IHC were identified as MSI-high by the algorithm. Conclusions: Our ML-driven approach accurately assessed MSI status of CRC and EMCA using CGIP data. This approach also identified potential cases with MSI-high status where direct sequencing of microsatellites failed. This study highlights a method to identify patients with potential MSI-high status for orthogonal screening with the MSI component of a test fails. [Table: see text]

Congenital FⅦ Deficiency Associated with a Novel Mutation in F7 Gene

To identify the genetic mutation of coagulation factor Ⅶ ( F7) gene in a pedigree with coagulation factor Ⅶ (FⅦ) deficiency and explore the molecular pathogenesis. The prothrombin time (PT), activated partial thromboplastin time (APTT), thrombin time (TT), D-dimer (DD), fibrin degradation products (FDP) and coagulation factor Ⅶ activity (FⅦ:C) of the proband and her family members were detected by Sysmex-CS5100 analyzer. All exons and exon-intron boundaries of the F7 gene were amplified by PCR followed by direct sequencing. The detected mutation was confirmed by reverse sequencing. The ClustalW software was used to analyze the conservatism of the mutant site. Pathogenicity of the mutation was assessed with Mutation Taster and PolyPhen-2 online bioinformatics software. Structure of the mutant protein was analyzed using Swiss-PdbViewer software. The results of routine coagulation tests showed that PT of the proband was markedly extended to 42.5 s, and her FⅦ:C significantly reduced to 2%. The FⅦ:C of her grandmother, mother and sister had slightly reduced to 49%, 51%, and 42%, respectively. These coagulant parameters of her father were within the normal range. Genetic analysis reveled a heterozygous G>A change at cDNA 646 in exon 6 of F7 gene in the proband, resulting in a replacement of glycine at 156 of FⅦ catalytic region with serine (p.Gly156Ser). The sequencing results of other exons and exon-intron boundaries of her F7 gene were normal. The proband's grandmother, mother and sister were all the carriers of this missense mutation except her father. Bioinformatics analysis showed that the p.Gly156Ser mutation caused polarity change of the amino acid at this site and formation of side chains, leading to increase of protein instability, which may affect catalytic activity of structural domain. Meanwhile, both Mutation Taster and PolyPhen-2 online bioinformatics software also predicted the pathogenicity of this missense mutation with high scores. The heterozygous p.Gly156Ser mutation is the direct cause of the reduced FⅦ in this proband.

Full-length PacBio Amplicon Sequencing to Unveil RNA Editing Sites

Background: RNA editing enriches post-transcriptional sequence changes. Currently detecting RNA editing sites is mostly based on the Sanger sequencing platform and second-generation sequencing. However, detection with Sanger sequencing is limited by the disturbing background peaks using the direct sequencing method and the clone number using the clone sequencing method, while second-generation sequencing detection is constrained by its short read. Objective: We aimed to design a pipeline that can accurately detect RNA editing sites for full-length long-read amplicons to meet the requirement when focusing on a few specific genes of interest. Method: We developed a novel high-throughput RNA editing sites detection pipeline based on the PacBio circular consensus sequences sequencing which is accurate with high-throughput and long-read coverage. We tested the pipeline on cytosolic malate dehydrogenase in the hard-shelled mussel Mytilus coruscus and further validated it using direct Sanger sequencing. Results: Data generated from the PacBio circular consensus sequences (CCS) amplicons in three mussels were first filtered by quality and then selected by open reading frame. After filtering, 225-2047 sequences of the three mussels, respectively, were used to identify RNA editing sites. With corresponding genomic DNA sequences, we extracted 227-799 candidate RNA editing sites excluding heterozygous sites. We further figured out 7-11 final RESs using a new error model specially designed for RNA editing site detection. The resulting RNA editing sites all agree with the validation using the Sanger sequencing. Conclusion: We report a near-zero error rate method in identifying RNA editing sites of long-read amplicons with the use of PacBio CCS sequencing.

MRNA psi profiling using nanopore DRS reveals cell type-specific pseudouridylation.

Pseudouridine (psi) is one of the most abundant human mRNA modifications generated from the isomerization of uridine via psi synthases, including TRUB1 and PUS7. Nanopore direct RNA sequencing combined with our recent tool, Mod-p ID, enables psi mapping, transcriptome-wide, without chemical derivatization of the input RNA and/or conversion to cDNA. This method is sensitive for detecting changes in positional psi occupancies across cell types, which can inform our understanding of the impact on gene expression. We sequenced, mapped, and compared the positional psi occupancy across six immortalized human cell lines derived from diverse tissue types. We found that lung-derived cells have the highest proportion of psi, while liver-derived cells have the lowest. Further, among a list of highly conserved sites across cell types, most are TRUB1 substrates and fall within the coding sequence. We find that these conserved psi positions correspond to higher levels of protein expression than expected, suggesting translation regulation. Interestingly, we identify cell type-specific sites of psi modification in ubiquitously expressed genes. We validate these sites by ruling out single-nucleotide variants, analyzing current traces, and performing enzymatic knockdowns of psi synthases. Finally, we characterize sites with multiple psi modifications on the same transcript (hypermodification type II) and found that these can be conserved or cell type specific. Among these, we discovered examples of multiple psi modifications within the same k-mer for the first time and analyzed the effect on current distribution. Our data support the hypothesis that motif sequence and the presence of psi synthase are insufficient to drive modifications, that psi modifications contribute to regulating translation and that cell type-specific trans-acting factors play a major role in driving pseudouridylation.

Analysis of the Class 1 Integrons, Carbapenemase Genes and Biofilm Formation Genes Occurrence in Acinetobacter baumannii Clinical Isolates.

Acinetobacter baumannii is a non-fermentative Gram-negative bacterium that can cause nosocomial infections in critically ill patients. Carbapenem-resistant A. baumannii (CRAB) has spread rapidly in clinical settings and has become a key concern. The main objective of this study was to identify the distribution of integrons and biofilm-formation-related virulence genes in CRAB isolates. A total of 269 A. baumannii isolates (219 isolates of CRAB and 50 isolates of carbapenem-sensitive A. baumannii (CSAB)) were collected. Carbapenemase genes (bla KPC, bla VIM, bla IMP, bla NDM, and bla OXA-23-like) and biofilm-formation-related virulence genes (abal, bfms, bap, and cusE) were screened with PCR. Class 1 integron was screened with PCR, and common promoters and gene cassette arrays were determined with restriction pattern analysis combined with primer walking sequencing. Whole-genome sequencing was conducted, and data were analyzed for a bla OXA-23-like-negative isolate. All 219 CRAB isolates were negative for bla KPC, bla VIM, bla IMP, and bla NDM, while bla OXA-23-like was detected in 218 isolates. The detection rates for abal, bfms, bap, and cusE in 219 CRAB were 93.15%, 63.93%, 88.13%, and 77.63%, respectively. Class 1 integron was detected in 75 CRAB (34.25%) and in 3 CSAB. The single gene cassette array aacA4-catB8-aadA1 with relatively strong PcH2 promoter was detected in class 1 integrons. The bla OXA-23-like-negative CRAB isolate was revealed to be a new sequence type (Oxford 3272, Pasteur 2520) carrying bla OXA-72, bla OXA-259, and bla ADC-26. In conclusion, bla OXA-23-like was the main reason for CRAB's resistance to carbapenems. A new (Oxford 3272, Pasteur 2520) CRAB sequence type carrying the bla OXA-72, bla OXA-259, and bla ADC-26 was reported.

Multicellular, IVT-derived, unmodified human transcriptome for nanopore-direct RNA analysis.

Nanopore direct RNA sequencing (DRS) enables measurements of RNA modifications. Modification-free transcripts are a practical and targeted control for DRS, providing a baseline measurement for canonical nucleotides within a matched and biologically derived sequence context. However, these controls can be challenging to generate and carry nanopore-specific nuances that can impact analysis. We produced DRS datasets using modification-free transcripts from in vitro transcription (IVT) of cDNA from six immortalized human cell lines. We characterized variation across cell lines and demonstrated how these may be interpreted. These data will serve as a versatile control and resource to the community for RNA modification analysis of human transcripts.

Sequencing accuracy and systematic errors of nanopore direct RNA sequencing

BackgroundDirect RNA sequencing (dRNA-seq) on the Oxford Nanopore Technologies (ONT) platforms can produce reads covering up to full-length gene transcripts, while containing decipherable information about RNA base modifications and poly-A tail lengths. Although many published studies have been expanding the potential of dRNA-seq, its sequencing accuracy and error patterns remain understudied.ResultsWe present the first comprehensive evaluation of sequencing accuracy and characterisation of systematic errors in dRNA-seq data from diverse organisms and synthetic in vitro transcribed RNAs. We found that for sequencing kits SQK-RNA001 and SQK-RNA002, the median read accuracy ranged from 87% to 92% across species, and deletions significantly outnumbered mismatches and insertions. Due to their high abundance in the transcriptome, heteropolymers and short homopolymers were the major contributors to the overall sequencing errors. We also observed systematic biases across all species at the levels of single nucleotides and motifs. In general, cytosine/uracil-rich regions were more likely to be erroneous than guanines and adenines. By examining raw signal data, we identified the underlying signal-level features potentially associated with the error patterns and their dependency on sequence contexts. While read quality scores can be used to approximate error rates at base and read levels, failure to detect DNA adapters may be a source of errors and data loss. By comparing distinct basecallers, we reason that some sequencing errors are attributable to signal insufficiency rather than algorithmic (basecalling) artefacts. Lastly, we generated dRNA-seq data using the latest SQK-RNA004 sequencing kit released at the end of 2023 and found that although the overall read accuracy increased, the systematic errors remain largely identical compared to the previous kits.ConclusionsAs the first systematic investigation of dRNA-seq errors, this study offers a comprehensive overview of reproducible error patterns across diverse datasets, identifies potential signal-level insufficiency, and lays the foundation for error correction methods.

Molecular Characterization of the Actin Gene and 5′ Flanking Sequence from Brown Macroalga Saccharina japonica (Laminariales, Phaeophyta)

The brown macroalga Saccharina japonica (Laminariales, Phaeophyta) is the most productive cultured seaweed in the world. In order to improve the biosafety of transgenic kelp, it is necessary to develop endogenous constitutive promoters, replacing those of virus origin. In this study, the housekeeping actin gene from S. japonica (SjACT) was found to contain three exons and two introns, representing a unique actin gene structure pattern in brown algae. Additionally, the 5′ upstream region was obtained using genome walking, and fused to the reporter gene lacZ or EGFP to construct promoter-detective vectors. Using an established genetic transformation system, kelps in different life-cycle stages were transformed. The detection results showed that, in the diploid sporophyte stage, the transient expression from the lacZ gene could be observed in the frond, stipe, or holdfast of kelps, indicating a manner of being non-tissue-specific. And, in the haploid gametophyte stage of S. japonica, the fluorescence of the expressed EGFP were detected in vivo in gametophyte cells of both genders. These results indicate that the promoter of the SjACT gene (pSjACT) functions in a constitutive manner and is expected to be a key endogenous element in the genetic manipulation of kelps.

Studying the Pathogenic Effect of New Nucleotide Changes in the Promoter Region of the TERT Gene in Patients with Glioblastoma Brain Tumor

Introduction: Brain tumors are considered to be one of the most dangerous types of cancer, despite decades of research and treatment efforts. The activation of the telomerase enzyme is a critical factor in the immortality of these cells. Mutations in the promoter region of the TERT gene, particularly in the promoter hot spots, can influence the binding of activating transcription factors and inhibit the binding of negative regulatory factors of the TERT gene, ultimately regulating the gene's expression. This study aimed to identify and investigate the mutations of the TERT gene promoter region in glioblastoma, a malignant brain tumor, using bioinformatics. Methods: A study was conducted using the case-control method where 35 patients with glioblastoma multiform brain tumor were examined along with 40 people who were examined as control samples. In this research, the Touchdown PCR technique and direct DNA sequencing were used to detect mutations in the promoter region of the TERT gene in patients with glioblastoma. Furthermore, bioinformatic analyses were conducted to investigate the pathogenic effect of nucleotide changes in this gene region. Results: In the core region of the TERT gene promoter, three-point mutations were identified (c.*146C>T، c.*144C>G and c.*143G>C). Two of these mutations were novel and have been observed for the first time in glioblastoma multiform. The remaining mutation was located in the promoter core's hot spot of the gene. This mutation was known as c.*146C>T. Conclusion: In this research, the harmful impact of TERT promoter region mutations as a biomarker for glioblastoma was examined. These findings suggest that mutations in regulatory regions, as well as coding sequences, could be significant contributory factors in the process of carcinogenesis.

Biochemical-free enrichment or depletion of RNA classes in real-time during direct RNA sequencing with RISER

The heterogeneous composition of cellular transcriptomes poses a major challenge for detecting weakly expressed RNA classes, as they can be obscured by abundant RNAs. Although biochemical protocols can enrich or deplete specified RNAs, they are time-consuming, expensive and can compromise RNA integrity. Here we introduce RISER, a biochemical-free technology for the real-time enrichment or depletion of RNA classes. RISER performs selective rejection of molecules during direct RNA sequencing by identifying RNA classes directly from nanopore signals with deep learning and communicating with the sequencing hardware in real time. By targeting the dominant messenger and mitochondrial RNA classes for depletion, RISER reduces their respective read counts by more than 85%, resulting in an increase in sequencing depth of 47% on average for long non-coding RNAs. We also apply RISER for the depletion of globin mRNA in whole blood, achieving a decrease in globin reads by more than 90% as well as an increase in non-globin reads by 16% on average. Furthermore, using a GPU or a CPU, RISER is faster than GPU-accelerated basecalling and mapping. RISER’s modular and retrainable software and intuitive command-line interface allow easy adaptation to other RNA classes. RISER is available at https://github.com/comprna/riser.

The Novel Long QT Syndrome Type 2-associated F129I Mutation in the KCNH2 Gene Significantly Affects I Kr Through the hERG1 Homomeric and Heteromeric Potassium Channels

Objective: The long QT syndrome type 2 is caused by the loss-of-function mutations in the KCNH2 gene, which encodes hERG1, the voltage-gated potassium channel. The hERG1 channels conduct rapid delayed rectifier K+ currents (I Kr) in the human cardiac tissue. KCNH2 encodes 2 main isoforms—hERG1a and hERG1b, which assemble to form the homomeric or heteromeric hERG1 channels. However, the functional characteristics of the heteromeric hERG1 channels in long QT syndrome type 2 are not clear. In this study, a novel mutation in the N-terminus of hERG1a (F129I) was identified in a proband of long QT syndrome type 2. The purpose of this study was to identify the electrophysiological change of homomeric and heteromeric hERG1 channels with the F129I-hERG1a. Methods: Candidate genes were screened by direct sequencing. F129I-hERG1a was cloned in the pcDNA3.1 vector by site-directed mutagenesis. Then, the wild-type (WT) hERG1a and/or F129I-hERG1a were transiently expressed in the HEK293 cells with or without hERG1b co-expression. The expression levels of the transgenes, cellular distribution of hERG1a and hERG1b, and the electrophysiological features of the homomeric and the heteromeric hERG1 channels with the WT-hERG1a or F129I-hERG1a were analyzed using whole-cell patch-clamp electrophysiology, western blotting, and immunofluorescence techniques. Results: The proband was clinically diagnosed with long QT syndrome type 2 and carried a heterozygous mutation c.385T>A (F129I) in the KCNH2 gene. Electrophysiology study proved that the F129I substitution in hERG1a significantly decreased I Kr in both the homomeric and heteromeric hERG1channels by 86% and 70%, respectively (WT-hERG1a (54.88 ± 18.74) pA/pF vs. F129I-hERG1a (7.34 ± 1.90) pA/pF, P < 0.001; WT-hERG1a/hERG1b (89.92 ± 24.51) pA/pF vs. F129I-hERG1a/hERG1b (26.54 ± 9.83) pA/pF, P < 0.001). The voltage dependence of IKr activation (V½ and k) was not affected by the mutation in both the homomeric and heteromeric hERG1 channels. The peak current densities and the kinetic characteristics of IKr were comparable for both WT/F129I-hERG1a and WT-hERG1a. The channel inactivation and deactivation analysis showed that F129I substitution did not affect deactivation of the homomeric hERG1a channel, but significantly accelerated the deactivation and recovery from inactivation of the heteromeric hERG1a/hERG1b channel based on the time constants of fast and slow recovery from deactivation F129I-hERG1a/hERG1b vs. WT-hERG1a/hERG1b (P < 0.05). Western blotting and immunofluorescence labeling experiments showed that maturation and intracellular trafficking of the F129I-hERG1a protein was impaired and potentially increased the ratio of hERG1b to hERG1a in the F129I-hERG1a/hERG1b tetramer channel, thereby resulting in electrophysiological changes characteristic of the long QT syndrome type 2 pathology. Conclusions: I Kr was significantly reduced in the homomeric and heteromeric hERG1 channels with F129I-hERG1a. The F129I mutation significantly accelerated the deactivation and recovery from inactivation of the heteromeric F129I-hERG1a/hERG1b channel. F129I-hERG1a exhibited impaired maturation and intracellular trafficking, thereby potentially increasing the ratio of the hERG1b to hERG1a stoichiometry in the hERG1 tetrameric channel. These changes demonstrated the importance of the heteromeric hERG1 channel in long QT syndrome type 2 pathophysiology.

Surface material analysis for human papillomavirus detection in nail Bowen's disease caused by HPV type 58.

Over the past few years, cases of human papillomavirus (HPV) infection in nail Bowen's disease have been reported. This disease presents diagnostic challenges due to its similarity to nail malignant melanoma, particularly with respect to the clinical manifestation of black nail streaks. While skin biopsy is usually employed for diagnosis, it is an invasive procedure. We report the case of a 52-year-old healthy Japanese male with a pigmented streak on the nail of the fourth finger of his right hand, which had extended from the central to the lateral nail fold within 4 months. Dermoscopic examination revealed a dark-brown pigmented band with splinter microhemorrhage. Clinically, nail Bowen's disease was suspected. The lesion was excised in strips under local anesthesia. Histopathological examination revealed hyperkeratosis, parakeratosis, papillomatosis, and dyskeratotic cells with atypical nuclei irregularly arranged. Immunohistochemistry using anti-HPV L1 antibody detected HPV-positive cells in the upper epidermis and stratum corneum of the nail matrix. Mucosal high-risk HPV type 58 DNA was detected from brush cytology of the keratotic surface prior to surgery, which was confirmed in formalin-fixed, paraffin-embedded excised samples using polymerase chain reaction (PCR) and subsequent direct DNA sequencing. Our case highlights HPV type 58 as a potential causative agent of nail Bowen's disease and shows that brush cytology of the surface material prior to excision may be a useful and less invasive way for mucosal high-risk HPV detection. PCR analysis of the nail surface could serve as a supplementary diagnostic tool for nail Bowen's disease.

The Effects of a Direct Instruction Curriculum Sequence on the Reading Comprehension Responses of Early Readers

Reading comprehension is critical for early learners’ academic success. Derived relations are an important feature of reading comprehension within the science of behavior because reading comprehension involves experientially acquired derived stimulus relations (Helou-Care, 2008; Neves in Behavioral Interventions 38:376–400, 2023, Park, 2005; da Silva in Experimental Analysis of Human Behavior Bulletin 24:1–8, 2006). That is, derived relations pertain to a learner demonstrating comprehension by making untrained connections to previously trained stimuli (e.g., words). For example, a learner who is trained that “kitten” is a type of “cat” and that all “cats meow” should then be able to demonstrate the untrained connection that “kittens meow.” The purpose of this study was to investigate the effects of the Corrective Reading Comprehension (CR) curriculum, which employs an instructional sequence involving derived relations, on reading comprehension responses entailing derived relations. Four typically developing first-grade students were the participants in this study and the students were divided into two dyads. The principal researcher instructed each dyad on each of the five lessons of CR until mastery was obtained and then conducted postprobes on the dependent measures that were metaphors, derived relations from sentences, derived relations from letters/numbers, and implicit/explicit comprehension questions. The researchers used a delayed multiple probe design across 1st-grade dyads to assess the effect of five CR lessons on derived relations and reading comprehension responses of the participants. Posttest assessments showed gains for all four participants. The article concludes with a discussion of the main findings, implications for future research, and the limitations of this study.

#314 Systemic lupus erythematous (SLE) and isolated multiple sclerosis: analysis of BAFF—VAR variant on gene TNSF13B

Abstract Background and Aims We report a case of young a woman (33 years old) affected by SLE since 2007. In this year, a renal biopsy showed class IV lupus nephritis and the patient received corticosteroids and immunosuppressant for 3 years. The patient remained thereafter in disease remission with normal renal function, no proteinuria and normal complement levels. In 2013 the patient performed brain MRI for sudden hearing loss. The MRI showed a pattern compatible with isolated multiple sclerosis (RIS), confirmed on subsequent radiographic checks. In 2023, due to the initial onset of renal failure, a new re-staging biopsy showed microscopic aspects compatible with class I/II lupus nephritis (INS/RPS Classification 2004, 2018); the activity index calculated with SLEDAI—2K was +4. A new brain MRI carried out without onset of new symptoms showed a new hyperintense region at the level of the head of the right caudate nucleus and enlargement of the intrasellar CSF content as per the condition of a “partially empty” sella. A genetic variant of the TNFSF13B gene was therefore suspected, since it is involved in encoding the cytochine and drug target B-cell activating factor (BAFF), associated with multiple sclerosis as well as SLE. We searched for the genetic variant in order to undertake target therapy early in our patient. Method We performed PCR amplification and direct sequencing of all coding exons plus all the intron-exon junctions of the analyzed gene, by automatic sequencer with SeqStudio Genetic Analyzer (analytical sensitivity and specificity >99%). Results We found the heterozygous presence of the BAFF-var variant in the 3' UTR region of the TNFSF13B gene (Heterozygous genotype for BAFF-var). This variant involves an insertion/deletion (GCTGT>A) at the level of the 3' UTR region with consequent production of an alternative polyadenylation site (APA). This APA leads to the formation of a shorter transcript (BAFF-var mRNA) which causes greater production of the encoded protein compared to the wild-type transcript (BAFF-WT mRNA). Conclusion The molecular analysis conducted highlighted the presence of the BAFF-var variant in heterozygosis which leads to increased levels of soluble BAFF in the serum, with consequent up-regulation of humoral immunity. It is reported in the literature1 that this variant is associated with both Multiple Sclerosis and SLE and may represent a valid therapeutic target being treatable with the monoclonal antibody Belimumab.

The Relationship between Intelligence and Mitochondrial DNA Polymorphism at Site of 15925 Base Pairs

Intelligence has a huge influences and is very closely related with human’s life. Because intelligence is indefinable, scholars have specified it through testing tools and named it intelligence quotient (IQ). In measuring IQ, researchers have found that IQ has a wide distribution and tried to find out the reasons by assessing the brain's functional level and environmental and genetic factors. One study reported an association between IQ and mitochondrial DNA (mtDNA) sequence 15,925, and this study investigated the association in Koreans. The number of participants was 57 in the mental retardation group without chromosomal abnormalities and 148 in the normal group, respectively. PCR technique and direct sequencing analysis were used to know participant’s mtDNA sequences. At 15,925 base, all participants had only wild type C and no variant type of m.15,925 C > T. We further studied variants in mtDNA sequences 15,926-15,927. At 15,926 base, the variant type of 15,926 C < T was found in only one case (0.7%) in the control group. At the 15,927 base, variant type of 15,927 G < A was found in 9 cases (6.1%) in the control group and 1 case (1.8%) in the retardation group, respectively. As a result, the polymorphisms of mtDNA sequences at bases 15,925-15,927 had no statistical significance in the Korean control and mental retardation group.

PGPg_finder: A comprehensive and user-friendly pipeline for identifying plant growth-promoting genes in genomic and metagenomic data

Identifying and comparing plant growth-promoting traits (PGPT) within whole-genome and metagenomic sequencing data can significantly advance agricultural research and promote sustainable crop production. This study introduces PGPg_finder, a comprehensive pipeline designed to annotate and compare PGPT from both whole-genome and metagenome sequencing datasets. This pipeline utilizes direct sequence annotation alongside de novo assembly methods to accurately detect PGPT. By cross-referencing sequences from the PLaBAse database, it identifies and quantifies the presence of these genes within the original datasets, facilitating an intuitive comparison of the abundance and distribution of PGPT across various samples. We evaluated the performance of PGPg_finder by analyzing genomes from five rhizobacterial strains: Paenibacillus vini, Paenibacillus polymyxa, Fictibacillus sp., Brevibacillus agri, and Bacillus cereus, and also metagenomic samples from bulk soils subjected to forest-to-pasture conversion in the Amazon rainforest. The genomic workflow revealed several genes associated with substrate utilization, abiotic stress neutralization, phosphate solubilization, and iron acquisition. It also identified genes unique to specific lineages, including those associated with colonization and plant-derived substrate usage in P. polymyxa, quorum sensing response and biofilm formation in P. vini, heavy metal detoxification and nitrogen acquisition in B. agri, and spore production and neutralizing biotic stress in B. cereus. The strain Fictibacillus sp. presented several unique genes related to surface attachment, stress response, xenobiotic degradation, phosphate solubilization, and phytohormone production. The use of PGPg_finder highlights its potential to uncover novel inoculants and strains. The metagenomic workflow distinguished plant-growth promotion gene profiles between soils from the Amazon rainforest and pasture, with the latter showing a profile more aligned with simple carbohydrate consumption, abiotic stress tolerance, motility and chemotaxis, and phosphorus mineralization. Native forests exhibited a profile associated with the degradation of complex organic matter, oxidative stress tolerance, xenobiotic degradation, bactericidal activity, iron acquisition, and volatile pathways. These findings underscore the effectiveness and sensitivity of PGPg_finder in accurately identifying and comparing PGPT genes, highlighting both commonalities and variations across samples. The application of this pipeline has the potential to significantly facilitate the identification of plant growth-promoting microbes.