Direct Sequencing Research Articles

Capturing clinically relevant Campylobacter attributes through direct whole genome sequencing of stool.

Campylobacter is the leading bacterial cause of infectious intestinal disease, but the pathogen typically accounts for a very small proportion of the overall stool microbiome in each patient. Diagnosis is even more difficult due to the fastidious nature of Campylobacter in the laboratory setting. This has, in part, driven a change in recent years, from culture-based to rapid PCR-based diagnostic assays which have improved diagnostic detection, whilst creating a knowledge gap in our clinical and epidemiological understanding of Campylobacter genotypes - no isolates to sequence. In this study, direct metagenomic sequencing approaches were used to assess the possibility of replacing genome sequences with metagenome sequences; metagenomic sequencing outputs were used to describe clinically relevant attributes of Campylobacter genotypes. A total of 37 diarrhoeal stool samples with Campylobacter and five samples with an unknown pathogen result were collected and processed with and without filtration, DNA was extracted, and metagenomes were sequenced by short-read sequencing. Culture-based methods were used to validate Campylobacter metagenome-derived genome (MDG) results. Sequence output metrics were assessed for Campylobacter genome quality and accuracy of characterization. Of the 42 samples passing quality checks for analysis, identification of Campylobacter to the genus and species level was dependent on Campylobacter genome read count, coverage and genome completeness. A total of 65% (24/37) of samples were reliably identified to the genus level through Campylobacter MDG, 73% (27/37) by culture and 97% (36/37) by qPCR. The Campylobacter genomes with a genome completeness of over 60% (n=21) were all accurately identified at the species level (100%). Of those, 72% (15/21) were identified to sequence types (STs), and 95% (20/21) accurately identified antimicrobial resistance (AMR) gene determinants. Filtration of stool samples enhanced Campylobacter MDG recovery and genome quality metrics compared to the corresponding unfiltered samples, which improved the identification of STs and AMR profiles. The phylogenetic analysis in this study demonstrated the clustering of the metagenome-derived with culture-derived genomes and revealed the reliability of genomes from direct stool sequencing. Furthermore, Campylobacter genome spiking percentages ranging from 0 to 2% total metagenome abundance in the ONT MinION sequencer, configured to adaptive sequencing, exhibited better assembly quality and accurate identification of STs, particularly in the analysis of metagenomes containing 2 and 1% of Campylobacter jejuni genomes. Direct sequencing of Campylobacter from stool samples provides clinically relevant and epidemiologically important genomic information without the reliance on cultured genomes.

Symmetrical acral keratoderma associated with new variants in the filaggrin gene.

Symmetrical acral keratoderma (SAK) is a rare skin disorder with symmetric hyperkeratotic patches on the acral regions. Variants in the filaggrin gene (FLG) have been associated with SAK since 2020. To explore the clinical and genetic basis in six patients with SAK. Whole-exome sequencing, direct sequencing, and prediction of protein structure and function were performed. In this study, we identified two novel variants, c.3320del and c.4909del, and seven previously reported variants, c.3099C>G, c.4544C>A, c.6950_6957del, c.7264G>T, c.7945del, c.8117C>G, c.12064A>T. The findings of this study bolster the existing evidence implicating FLG variants in SAK, introducing two novel variants to the database of FLG variants associated with the condition.

Single Nucleotide Polymorphisms in APE1, hOGG1, RAD51 Genes and their Association with Radiotherapy Induced Toxicity among Head and Neck Cancer Patients.

Radiotherapy (RT) is a crucial treatment for head and neck cancer however, it causes adverse reactions to the normal tissue and organs adjacent to target tumor. The present study was carried out to investigate possible association of single nucleotide polymorphism in DNA repair genes with toxicity effects of radiotherapy on normal tissue. Three hundred and fifty head and neck cancer patients receiving radiotherapy treatment were enrolled in this study. The adverse after effects of radiotherapy on the normal tissue in the form of skin reactions were recorded. Single nucleotide polymorphisms of APE1 (rs1130409), hOGG1 (rs1052133) and Rad51 (rs1801320, rs1801321) genes were studied by polymerase chain reaction-Restriction fragment length polymorphism (PCR-RFLP) and direct DNA sequencing methods and their association with development of severe radio-toxicity effects was evaluated logistic regression analysis. The 172G/T polymorphism of Rad51 was 2.85 times higher and significantly associated with skin reactions (OR=2.85, 95% CI: 1.50-5.41; p=0.001) and severe oral mucositis (OR=4.96, 95% CI: 2.40-10.25; p<0.0001). These results suggested that the polymorphic nature of Rad51 is responsible for risk of radiotherapy adverse effects in HNC patients. The variant 326Cys and heterozygous 326Ser/Cys genotype of hOGG1 was significantly associated with high tumor grade (OR=3.16 95% CI: 1.66-5.99; p=0.0004, and OR=3.97 95% CI: 2.15-7.34; p=<0.0001 respectively). The homozygous variant 172TT genotype of Rad51 showed positive association with poor response of both tumor and nodes towards radiotherapy treatment (p=0.007 and p=0.022). Interpretation of our results revealed significant association of rs1801321 SNP of Rad51 with development of adverse toxicity reactions in normal tissue of head and neck cancer patients treated with radiotherapy.

Genetic Heterogeneity in Cellular Angiofibromas.

Cellular angiofibroma, a rare benign mesenchymal neoplasm, is classified within the 13q/RB1 family of tumors due to morphological, immunohistochemical, and genetic similarities with spindle cell lipoma. Here, genetic data reveal pathogenetic heterogeneity in cellular angiofibroma. Three cellular angiofibromas were studied using G-banding/Karyotyping, array comparative genomic hybridization, RNA sequencing, and direct cycling sequencing. The first tumor carried a del(13)(q12) together with heterozygous loss and minimal expression of the RB1 gene. Tumors two and three displayed chromosome 8 abnormalities associated with chimeras of the pleomorphic adenoma gene 1 (PLAG1). In tumor 2, the cathepsin B (CTSB) fused to PLAG1 (CTSB::PLAG1) while in tumor 3, the mir-99a-let-7c cluster host gene (MIR99AHG) fused to PLAG1 (MIR99AHG::PLAG1), both leading to elevated expression of PLAG1 and insulin growth factor 2. This study uncovers two genetic pathways contributing to the pathogenetic heterogeneity within cellular angiofibromas. The first aligns with the 13q/RB1 family of tumors and the second involves PLAG1-chimeras. These findings highlight the diverse genetic landscape of cellular angiofibromas, providing insights into potential diagnostic strategies.

Linking Survivor Stories to Forensic Engineering: How an Interscience Approach Reveals Opportunities for Reducing Tornado Vulnerability in Residential Structures

Abstract When a tornado strikes a permanent or mobile/manufactured home, occupants are at risk of injury and death from blunt force trauma caused by debris-loaded winds and failure of the structure. Mechanisms for these failures have been studied for the past few decades and identified common weaknesses in the structural load path. Also under study in recent decades, much has been learned about how people receive and understand warnings and determine how, when, and if they will shelter in advance. Recent research, for example, shows most people do not shelter until close to impact, after seeing, hearing, or feeling the approaching tornado. To advance beyond these innovations, a new, multidisciplinary approach was fielded in nine southeast U.S. tornadoes between 2019 and 2022. For each tornado, 1) wind engineering assessments documented near-surface wind fields, 2) structural engineering assessments documented the primary wind load path for each structure, and 3) social science interviews captured the survivor’s narrative and asked several follow-up questions to assure key items of interest were addressed in each interview. When possible, the team was multidisciplinary during the interview, enabling survivors to ask questions and better understand their experiences. Most survivors became aware of the approaching tornado with at least a few minutes of lead time, and most were able to reach a place of refuge. Most survivors recalled sensory experiences during the tornado, and about half could describe the direction or temporal sequences of damage. A case study of the Cookeville, Tennessee, tornado of 3 March 2020 illustrates the power of the integrated data assessment.

Long-range and real-time PCR identification of a large SERPINC1 deletion in a patient with antithrombin deficiency.

Congenital antithrombin (AT) or serpin C1 deficiency, caused by a SERPINC1 abnormality, is a high-risk factor for venous thrombosis. SERPINC1 is prone to genetic rearrangement, because it contains numerous Alu elements. In this study, a Japanese patient who developed deep vein thrombosis during pregnancy and exhibited low AT activity underwent SERPINC1 gene analysis using routine methods: long-range polymerase chain reaction (PCR) and real-time PCR. Sequencing using long-range PCR products revealed no pathological variants in SERPINC1 exons or exon-intron junctions, and all the identified variants were homozygous, suggesting a deletion in one SERPINC1 allele. Copy number quantification for each SERPINC1 exon using real-time PCR revealed half the number of exon 1 and 2 copies compared with controls. Moreover, a deletion region was deduced by quantifying the 5'-upstream region copy number of SERPINC1 for each constant region. Direct long-range PCR sequencing with primers for the 5'-end of each presumed deletion region revealed a large Alu-mediated deletion (∼13kb) involving SERPINC1 exons 1 and 2. Thus, a large deletion was identified in SERPINC1 using conventional PCR methods.

A Rare Combination of Compound Heterozygous Mutations in the PAH Gene in Three Unrelated Consanguineous Iranian Families with Classical Phenylketonuria.

The PAH gene mutations have been linked to the development of phenylketonuria (PKU), which is recognized as the most common inborn metabolic disorder, and is caused by a deficiency in the phenylalanine hydroxylase (PAH) enzyme. The Iranian population, known for its diversity and high consanguinity, offers a valuable sample for studying autosomal recessive disorders. Our study investigated three unrelated families with PKU from Iran, utilizing clinical, laboratory, and computational approaches. We performed direct PCR sequencing for 13 exons of the PAH gene on three Iranian patients who were diagnosed with PKU. Then, Sanger sequencing confirmed the segregation of the mutations from parents to probands. Pathogenicity predictor tools, including ACMG, CADD, SIFT, Polyphen-2, and Mutation Taster, were utilized to analyze the identified genetic variants. The three-dimensional structure of the mutant forms of the protein was predicted. We also analyzed the protein-protein interactions of PAH using the STRING database. All three patients exhibited rare compound heterozygosity rearrangements in the PAH gene (NM_000277.3). These included three missense variants: [c.533A>G/c.1222C>T], [c.526C>T/c.1222C>T], and [c.533A> / c.526C>T]. This study adds to the body of evidence establishing the association between PAH mutations and the development of PKU. We speculated that the conjunction of a high consanguinity rate in populations such as Iran, coupled with the founder effect, can give rise to atypical genetic profiles, as observed in the rearrangement of compound heterozygosity in this study Moreover, our research underscores the significance of genetic testing in the precise diagnosis of individuals affected by inborn errors of metabolism.

Dengue virus type 3 infection in a traveler returning from Costa Rica to Japan in 2023

The number of dengue cases has increased dramatically in recent years. In Latin America, the number of cases and deaths in 2023 was the highest ever recorded. We report on a patient who had been infected with dengue virus during his stay in Costa Rica in September 2023, and developed the disease after returning to Japan. Plasma obtained from the patient was used for diagnosis and dengue virus serotyping by real-time PCR. The nucleotide sequence of the envelope region of dengue virus was then determined by the direct sequencing method, and this sequence was used for phylogenetic analyses. The patient was found to be infected with dengue virus type 3 genotype III. The sequence from the present case was more homologous with sequences registered in Florida, USA, associated with travel to Cuba in 2022 than with sequences registered in Costa Rica 10 years ago. The Pan American Health Organization reported that only dengue virus type 1 and 2 cases were reported in Costa Rica in 2019–2021, whereas dengue virus type 3 and 4 cases started being reported in 2022. In 2023, the reported numbers of cases with dengue virus types 3 and 4 exceeded those of dengue virus types 1 and 2. In addition, regional differences in endemic strains have been observed in Costa Rica. Our findings suggest that the dengue virus type 3 that infected the patient was more likely an influx of a strain that had been circulating in Caribbean countries such as Cuba in recent years, rather than a re-emergence of an indigenous virus in Costa Rica. The serotypes of dengue virus prevalent in Costa Rica have been changing since 2022. All four serotypes were prevalent in 2023, with a particularly sharp increase in the number of cases of dengue virus types 3 and 4. Future monitoring and surveillance are essential because changes in endemic serotypes can cause antibody-dependent enhancement, which can lead to severe dengue disease presentations.

Detection of HPV DNA in Saliva of Patients with HPV-Associated Oropharyngeal Cancer Treated with Radiotherapy.

To investigate the technical feasibility of RT-PCR and direct sequencing to quantify HPV DNA in the saliva of patients with Human-Papilloma-Virus related oropharyngeal cancer (HPV-OPC), the level of which is known to predict prognosis after treatment. Nine patients with locally advanced HPV-OPC treated with definitive radiotherapy with chemotherapy or cetuximab, or radiotherapy alone between April 2016 and September 2017, were enrolled, two of whom also received induction chemotherapy. Saliva was collected before (baseline), during (mid-RT) and after (post-RT) radiotherapy. HPV-16 DNAs (E6 and E7) in saliva were quantified by RT-PCR and sequencing, the latter using a custom cancer panel. Correlations between HPV DNA levels and clinical outcomes were assessed. Compared to the baseline, the relative cycle threshold (Ct) value of E6 and E7 reduced at the point of mid-RT in the majority of the patients (100% and 75% for E6 and E7, respectively). Similarly, the relative Ct value from the baseline to post-RT reduced in 86% and 100% of the patients for E6 and E7, respectively. During the follow-up period, three patients (33%) experienced disease progression. The relative baseline Ct values of these three patients were in the top 4 of all the patients. The sequences of HPV DNA were detected in five (83%) of six samples of the baseline saliva that underwent DNA sequencing, along with several gene mutations, such as TP53,CDKN2A and PIK3CA. This study demonstrates that, in addition to detection and quantification of HPV DNA by RT-PCR, detection by sequencing of HPV-DNA using a customized cancer panel is technically possible.

Clinical Characteristics and Prognostic Relevance of Co-Mutated Genes in Acute Myeloid Leukemia Patients with FLT3 Mutations

To investigate the clinical characteristics and influence of co-mutated gene on acute myeloid leukemia patients (AML) with FMS-like tyrosine kinase-3 (FLT3) mutations. A total of 273 FLT3+ AML patients were enrolled, and the co-mutation gene data of the patients were collected to further analyze the prognosis of the patients. FLT3 and other common mutations were quantified by PCR amplification products direct sequencing and second-generation sequencing (NGS). When patients were divided into FLT3- ITD +, FLT3- TKD +, FLT3- ITD ++TKD + and FLT3- ITD -+TKD - group according to the type of FLT3 mutations, it was found that the frequencies of TET2, GATA2, NRAS and ASXL1 mutation were significantly different among the 4 groups (all P < 0.05). When patients were divided into allelic ratio (AR) ≥0.5 and <0.5 group, it was found that the frequencies of FLT3- ITD +, FLT3 -ITD - +TKD -, NPM1, NRAS and C-kit were significantly different between the two groups (all P < 0.05). When patients were divided into normal and abnormal karyotype group, it was found that the frequencies of FLT3- ITD +, FLT3- TKD +, NPM1, GATA2 and C-kit were significantly different between the two groups (all P < 0.05). The median overall survival (OS) of AML patients with FLT3 -TKD + (including FLT3- ITD ++TKD +) was longer than that of patients with FLT3- ITD + alone (P < 0.05). The OS and relapse-free survival (RFS) of AML patients with FLT3++TET2+ were both shorter than those of patients with FLT3++TET2- (both P < 0.05). The mutation frequencies of co-mutated genes are correlated with subtypes of FLT3, karyotype and AR. AML patients with FLT3 -TKD + have longer OS than patients with FLT3- ITD + alone, and patients with co-mutation of TET2 have shorter median OS and RFS.

Dysfibrinogenemia and hypofibrinogenemia - Spectrum of pathogenic variants in Slovak patients.

Congenital hypofibrinogenemia (CH) and congenital dysfibrinogenemia (CD) are rare coagulation disorders caused by quantitative or qualitative defects in the fibrinogen gene. The aim of this study was to characterize the genetic background and the clinical manifestations of congenital fibrinogen disorders in the patients from Slovakia registered at the National Haemophilia Centre. Results of genetic analysis of the fibrinogen genes FGA, FGB and FGG using polymerase chain reaction followed by direct sequencing were evaluated in 36 patients. Molecular-genetic analysis revealed six novel variants - FGA c.923_968dup p.(Gly324Lysfs*44) and FGG c.1105C>T p.(His369Tyr) were identified in CD patients. In CH patients, in the FGG gene c.8G>A p.(Trp3*), c.823G>T p.(Glu275*) and c.323C>A p.(Ala108Asp) variants were detected. In the FGB gene c.1427C>T p.(Ser476Leu) was identified. This study is a positive contribution towards expanding knowledge about genetic variants in patients with congenital fibrinogen disorders.

Association of a rapidly selected 4.3kb transposon-containing structural variation with a P450-based resistance to pyrethroids in the African malaria vector Anopheles funestus.

Deciphering the evolutionary forces controlling insecticide resistance in malaria vectors remains a prerequisite to designing molecular tools to detect and assess resistance impact on control tools. Here, we demonstrate that a 4.3kb transposon-containing structural variation is associated with pyrethroid resistance in central/eastern African populations of the malaria vector Anopheles funestus. In this study, we analysed Pooled template sequencing data and direct sequencing to identify an insertion of 4.3kb containing a putative retro-transposon in the intergenic region of two P450s CYP6P5-CYP6P9b in mosquitoes of the malaria vector Anopheles funestus from Uganda. We then designed a PCR assay to track its spread temporally and regionally and decipher its role in insecticide resistance. The insertion originates in or near Uganda in East Africa, where it is fixed and has spread to high frequencies in the Central African nation of Cameroon but is still at low frequency in West Africa and absent in Southern Africa. A marked and rapid selection was observed with the 4.3kb-SV frequency increasing from 3% in 2014 to 98% in 2021 in Cameroon. A strong association was established between this SV and pyrethroid resistance in field populations and is reducing pyrethroid-only nets' efficacy. Genetic crosses and qRT-PCR revealed that this SV enhances the expression of CYP6P9a/b but not CYP6P5. Within this structural variant (SV), we identified putative binding sites for transcription factors associated with the regulation of detoxification genes. An inverse correlation was observed between the 4.3kb SV and malaria parasite infection, indicating that mosquitoes lacking the 4.3kb SV were more frequently infected compared to those possessing it. Our findings highlight the underexplored role and rapid spread of SVs in the evolution of insecticide resistance and provide additional tools for molecular surveillance of insecticide resistance.

Использование алгоритмов Process Mining для выявления паттернов поведения студентов

This work is devoted to studying the possibilities of using process analytics methods (Process Mining) to analyze student activity based on digital traces that students leave in learning management systems (LMS). This work examines the specifics of process mining algorithms that can be used to analyze educational processes, namely, heuristic and inductive algorithms are considered as the most effective for building models and suitable for use for the purpose of analyzing educational data. The work involved creating a way to use process analytics algorithms to identify clusters of students with similar behavior patterns. The development of a process analysis algorithm was carried out on the basis of the event log of the distance learning system of Kostroma State University. As a result of the work, models of student behavior were built and visualized, including the identification and clustering of students with similar behavior, the construction of heuristic networks, Petri nets, a direct sequence graph, a BPMN model and a decision tree. An analysis of the resulting models was carried out, which showed that the developed method makes it possible to study the behavioral patterns of students. The proposed method of using intellectual analysis of educational processes can be used to solve issues of increasing the productivity of the educational process, early detection of problems, especially in the context of changing student behavior in the system, as well as the development and optimization of educational programs. In addition, the limitations of this system have been identified, which may hinder its implementation and application in the educational environment of universities.

Comparative analysis of 43 distinct RNA modifications by nanopore tRNA sequencing.

Transfer RNAs are the fundamental adapter molecules of protein synthesis and the most abundant and heterogeneous class of noncoding RNA molecules in cells. The study of tRNA repertoires remains challenging, complicated by the presence of dozens of post transcriptional modifications. Nanopore sequencing is an emerging technology with promise for both tRNA sequencing and the detection of RNA modifications; however, such studies have been limited by the throughput and accuracy of direct RNA sequencing methods. Moreover, detection of the complete set of tRNA modifications by nanopore sequencing remains challenging. Here we show that recent updates to nanopore direct RNA sequencing chemistry (RNA004) combined with our own optimizations to tRNA sequencing protocols and analysis workflows enable high throughput coverage of tRNA molecules and characterization of nanopore signals produced by 43 distinct RNA modifications. We share best practices and protocols for nanopore sequencing of tRNA and further report successful detection of low abundance mitochondrial and viral tRNAs, providing proof of concept for use of nanopore sequencing to study tRNA populations in the context of infection and organelle biology. This work provides a roadmap to guide future efforts towards de novo detection of RNA modifications across multiple organisms using nanopore sequencing.

Characteristics and Driving Mechanisms of Coastal Wind Speed during the Typhoon Season: A Case Study of Typhoon Lekima

The development and utilization of wind energy is of great significance to the sustainable development of China’s economy and the realization of the “dual carbon” goal. Under typhoon conditions, the randomness and volatility of wind speed significantly impact the energy efficiency and design of wind turbines. This paper analyzed the changes in wind speed and direction using the BFAST method and Hurst index based on data collected at 10 m, 30 m, 50 m, and 70 m heights from a wind power tower in Yancheng, Jiangsu Province. Furthermore, the paper examined the causes of wind speed and direction changes using wind speed near the typhoon center, distance from the typhoon center to the wind tower, topographic data, and mesoscale system wind direction data. The conclusions drawn are as follows: (i) Using the BEAST method, change points were identified at 10 m, 30 m, 50 m, and 70 m heights, with 5, 5, 6, and 6 change points respectively. The change points at 10 m, 30 m, and 50 m occurred around node 325, while the change time at 70 m was inconsistent with other heights. Hurst index results indicated stronger inconsistency at 70 m altitude compared to other altitudes. (ii) By analyzing the wind direction sequence at 10 m, 30 m, 50 m, and 70 m, it was found that the wind direction changes follow the sequence Southeast (SE)—East (E)—Southeast (SE)—Southwest (SW)—West (W)—Northwest (NW). Notably, the trend of wind direction at 70 m significantly differed from other altitudes during the wind speed strengthening and weakening stages. (iii) Wind speed at 10 m and 70 m altitudes responded differently to the distance from the typhoon center and the wind near the typhoon center. The correlation between wind speed and the distance to the typhoon center was stronger at 10 m than at 70 m. The surface type and the mesoscale system’s wind direction also influenced the wind speed and direction. This study provides methods and theoretical support for analyzing short-term wind speed changes during typhoons, offering reliable support for selecting wind power forecast indicators and designing wind turbines under extreme gale weather conditions.

The Advantage of Targeted Next-Generation Sequencing over qPCR in Testing for Druggable EGFR Variants in Non-Small-Cell Lung Cancer.

The emergence of targeted therapies in non-small-cell lung cancer (NSCLC), including inhibitors of epidermal growth factor receptor (EGFR) tyrosine kinase, has increased the need for robust companion diagnostic tests. Nowadays, detection of actionable variants in exons 18-21 of the EGFR gene by qPCR and direct DNA sequencing is often replaced by next-generation sequencing (NGS). In this study, we evaluated the diagnostic usefulness of targeted NGS for druggable EGFR variants testing in clinical NSCLC material previously analyzed by the IVD-certified qPCR test with respect to DNA reference material. We tested 59 NSCLC tissue and cytology specimens for EGFR variants using the NGS 'TruSight Tumor 15' assay (Illumina) and the qPCR 'cobas EGFR mutation test v2' (Roche Diagnostics). The sensitivity and specificity of targeted NGS assay were evaluated using the biosynthetic and biological DNA reference material with known allelic frequencies (VAF) of EGFR variants. NGS demonstrated a sufficient lower detection limit for diagnostic applications (VAF < 5%) in DNA reference material; all EGFR variants were correctly identified. NGS showed high repeatability of VAF assessment between runs (CV% from 0.02 to 3.98). In clinical material, the overall concordance between NGS and qPCR was 76.14% (Cohen's Kappa = 0.5933). The majority of discordant results concerned false-positive detection of EGFR exon 20 insertions by qPCR. A total of 9 out of 59 (15%) clinical samples showed discordant results for one or more EGFR variants in both assays. Additionally, we observed TP53 to be a frequently co-mutated gene in EGFR-positive NSCLC patients. In conclusion, targeted NGS showed a number of superior features over qPCR in EGFR variant detection (exact identification of variants, calculation of allelic frequency, high analytical sensitivity), which might enhance the basic diagnostic report.

High level of genomic divergence in orf-I p12 and hbz genes of HTLV-1 subtype-C in Central Australia

BackgroundHuman T cell lymphotropic virus type 1 (HTLV-1) infection remains a largely neglected public health problem, particularly in resource-poor areas with high burden of communicable and non-communicable diseases, such as some remote populations in Central Australia where an estimated 37% of adults are infected with HTLV-1. Most of our understanding of HTLV-1 infection comes from studies of the globally spread subtype-A (HTLV-1a), with few molecular studies reported with the Austral-Melanesian subtype-C (HTLV-1c) predominant in the Indo-Pacific and Oceania regions.ResultsUsing a primer walking strategy and direct sequencing, we constructed HTLV-1c genomic consensus sequences from 22 First Nations participants living with HTLV-1c in Central Australia. Phylogenetic and pairwise analysis of this subtype-C proviral gDNA showed higher levels of genomic divergence in comparison to previously published HTLV-1a genomes. While the overall genomic homology between subtypes was 92.5%, the lowest nucleotide and amino acid sequence identity occurred near the 3′ end of the proviral genome coding regulatory genes, especially overlapping hbz (85.37%, 77.46%, respectively) and orf-I product p12 (82.00%, 70.30%, respectively). Strikingly, the HTLV-1c genomic consensus sequences uniformly showed a defective translation start codon for the immune regulatory proteins p12/p8 encoded by the HTLV-1A orf-I. Deletions in the proviral genome were detected in many subjects, particularly in the structural gag, pol and env genes. Similarly, using a droplet digital PCR assay measuring the copies of gag and tax per reference host genome, we quantitatively confirmed that provirus retains the tax gene region at higher levels than gag.ConclusionsOur genomic analysis of HTLV-1c in Central Australia in conjunction with earlier Melanesian HTLV-1c sequences, elucidate substantial differences with respect to the globally spread HTLV-1a. Future studies should address the impact these genomic differences have on infection and the regionally distinctive frequency of associated pulmonary disease. Understanding the host and virus subtype factors which contribute to the differential morbidity observed, is crucial for the development of much needed therapeutics and vaccine strategies against this highly endemic infection in remote First Nations communities in Central Australia.

Methodological approaches to gene identification of tea raw materials and raw material composition of tea-based soft drinks

Tea or tea shrub is a plant of the Camellia sinensis species, the leaves of which, previously prepared in a special way, are the traditional raw material for the production of tea products. Varietal gene identification of tea allows us to increase the level of assessment of the authenticity of tea raw materials and tea products. It is predominantly based on DNA technologies for the detection and interpretation of SNP markers (Single Nucleotide Polymorphism), represented by a wide arsenal of both expensive high-tech methods and publicly available laboratory approaches. Species gene identification of the raw material composition of tea-based soft drinks is an equally important area of research due to the risk of falsification of this type of product. The purpose of this study was to find methodological approaches to the varietal gene identification of tea raw materials and to the species gene identification of the raw material composition of tea-based soft drinks. As a result of a bioinformatics study to identify polymorphic restriction sites in the nucleotide sequences of Camellia sinensis genome loci, diagnostically significant restriction enzymes were selected that were capable of detecting SNPs and identifying tea genotypes using the analyzed markers. At the same time, 16 loci had potential for practical application, of which 11 belonged to the group of the most informative SNP markers. A post-analytical assessment of tea varieties was carried out with them regarding their genotypic affiliation and identifiability as part of solving the first task of the study. To achieve the second task, a molecular genetic approach to the species identification of the raw composition of soft drinks based on green tea was tested. The study included the analysis of experimental drinks (with natural flavoring “Lemon” and synthetic flavoring “Peach 716”), as well as commercial concentrates “TIAKVA” (based on extracts from the coarse stems of green or black tea). The methods used in the work were PCR (Polymerase Chain Reaction), RFLP (Restriction Fragment Length Polymorphism) and direct sequencing of the amplified chloroplast DNA locus. The combination of two methods (PCR and sequencing) showed its effectiveness in establishing the belonging of the analyzed nucleic acid samples to the Camellia sinensis species, the raw material base of the studied drinks and concentrates. However, to unlock the authentication potential of PCR with primers #1 and #2 combined with RFLP analysis, it will be necessary to select diagnostically significant restriction enzymes suitable for generating species-specific combinations of PCR-RFLP profiles of marker sequence.

Analysis of bacterial transcriptome and epitranscriptome using nanopore direct RNA sequencing.

Bacterial gene expression is a complex process involving extensive regulatory mechanisms. Along with growing interests in this field, Nanopore Direct RNA Sequencing (DRS) provides a promising platform for rapid and comprehensive characterization of bacterial RNA biology. However, the DRS of bacterial RNA is currently deficient in the yield of mRNA-mapping reads and has yet to be exploited for transcriptome-wide RNA modification mapping. Here, we showed that pre-processing of bacterial total RNA (size selection followed by ribosomal RNA depletion and polyadenylation) guaranteed high throughputs of sequencing data and considerably increased the amount of mRNA reads. This way, complex transcriptome architectures were reconstructed for Escherichia coliand Staphylococcus aureusand extended the boundaries of 225 known E. coli operons and 89 defined S. aureus operons. Utilizing unmodified in vitro-transcribed (IVT) RNA libraries as a negative control, several Nanopore-based computational tools globally detected putative modification sites in the E. coli and S. aureus transcriptomes. Combined with Next-Generation Sequencing-based N6-methyladenosine (m6A) detection methods, 75 high-confidence m6A candidates were identified in the E. coli protein-coding transcripts, while none were detected in S. aureus. Altogether, we demonstrated the potential of Nanopore DRS in systematic and convenient transcriptome and epitranscriptome analysis.

Improved characterization of single-cell RNA-seq libraries with paired-end avidity sequencing.

Prevailing poly(dT)-primed 3' single-cell RNA-seq protocols generate barcoded cDNA fragments containing the reverse transcriptase priming site, which is expected to be the poly(A) tail or a genomic adenine homopolymer. Direct sequencing across this priming site was historically difficult because of DNA sequencing errors induced by the homopolymeric primer at the 'barcode' end. Here, we evaluate the capability of "avidity base chemistry" DNA sequencing from Element Biosciences to sequence through this homopolymer accurately, and the impact of the additional cDNA sequence on read alignment and precise quantification of polyadenylation site usage. We find that the Element Aviti instrument sequences through the thymine homopolymer into the subsequent cDNA sequence without detectable loss of accuracy. The resulting paired-end alignments enable direct and independent assignment of reads to polyadenylation sites, which bypasses complexities and limitations of conventional approaches but does not consistently improve read mapping rates compared to single-end alignment. We also characterize low-level artifacts and arrive at an adjusted adapter trimming and alignment workflow that significantly improves the alignment of sequence data from Element and Illumina, particularly in the context of extended read lengths. Our analyses confirm that Element avidity sequencing is an effective alternative to Illumina sequencing for standard single-cell RNA-seq, particularly for polyadenylation site analyses but do not rule out the potential for similar performance from other emerging platforms.