Biostatistical Aspects of Whole Genome Sequencing Studies: Preprocessing and Quality Control.

Rapid advances in high-throughput DNA sequencing technologies have enabled large-scale whole genome sequencing (WGS) studies. Before performing association analysis between phenotypes and genotypes, preprocessing and quality control (QC) of the raw sequence data need to be performed. Because many biostatisticians have not been working with WGS data so far, we first sketch Illumina's short-read sequencing technology. Second, we explain the general preprocessing pipeline for WGS studies. Third, we provide an overview of important QC metrics, which are applied to WGS data: on the raw data, after mapping and alignment, after variant calling, and after multisample variant calling. Fourth, we illustrate the QC with the data from the GENEtic SequencIng Study Hamburg-Davos (GENESIS-HD), a study involving more than 9000 human whole genomes. All samples were sequenced on an Illumina NovaSeq 6000 with an average coverage of 35× using a PCR-free protocol. For QC, one genome in a bottle (GIAB) trio was sequenced in four replicates, and one GIAB sample was successfully sequenced 70 times in different runs. Fifth, we provide empirical data on the compression of raw data using the DRAGEN original read archive (ORA). The most important quality metrics in the application were genetic similarity, sample cross-contamination, deviations from the expected Het/Hom ratio, relatedness, and coverage. The compression ratio of the raw files using DRAGEN ORA was 5.6:1, and compression time was linear by genome coverage. In summary, the preprocessing, joint calling, and QC of large WGS studies are feasible within a reasonable time, and efficient QC procedures are readily available.

How advancements in molecular biology impact education and training.

Molecular biology, broadly defined as the investigation of complex biomolecules in the laboratory, is a rapidly advancing field and as such the technologies available to investigators are constantly evolving. This constant advancement has obvious advantages because it allows students and researchers to perform more complex experiments in shorter periods of time. One challenge with such a rapidly advancing field is that techniques that had been vital for students to learn how to perform are now not essential for a laboratory scientist. For example, while cloning a gene in the past could have led to a publication and form the bulk of a PhD thesis project, technology has now made this process only a step toward one of these larger goals and can, in many cases, be performed by a company or core facility. As teachers and mentors, it is imperative that we understand that the technologies we teach in the lab and classroom must also evolve to match these advancements. In this perspective, we discuss how the rapid advances in gene synthesis technologies are affecting curriculum and how our classrooms should evolve to ensure our lessons prepare students for the world in which they will do science.

Recent Innovations and Technical Advances in High-Throughput Parallel Single-Cell Whole-Genome Sequencing Methods.

Single-cell whole-genome sequencing (scWGS) detects cell heterogeneity at the aspect of genomic variations, which are inheritable and play an important role in life processes such as aging and cancer progression. The recent explosive development of high-throughput single-cell sequencing methods has enabled high-performance heterogeneity detection through a vast number of novel strategies. Despite the limitation on total cost, technical advances in high-throughput single-cell whole-genome sequencing methods are made for higher genome coverage, parallel throughput, and level of integration. This review highlights the technical advancements in high-throughput scWGS in the aspects of strategies design, data efficiency, parallel handling platforms, and their applications on human genome. The experimental innovations, remaining challenges, and perspectives are summarized and discussed.

A novel variant in the SLCO2A1 gene in a Chinese patient with chronic gastroenteropathy and primary hypertrophic osteoarthropathy.

BACKGROUND: Chronic enteropathy associated with SLCO2A1 gene (CEAS) results from loss-of-function variants in SLCO2A1, which encodes the prostaglandin transporter (PGT). CEAS follows an autosomal recessive inheritance pattern. To date, approximate 30 pathogenic variants have been reported in CEAS. METHODS: We performed whole exome sequencing (WES) to screen for potential pathogenic variants in a patient suspected of having CEAS, and confirmed a variant in SLCO2A1 using Sanger sequencing. We established an in vitro minigene model to compare splicing between wild type (WT) and mutant transcripts. Quantitative polymerase chain reaction (qPCR) was used to evaluate SLCO2A1 transcription in the stomach and colon tissues from the patient and a healthy control (HC). The transcripts were further cloned and sequenced. RESULTS: The patient had a novel, homozygous, recessive c.929A > G variant in exon 7 of SLCO2A1, which has not been previously reported in CEAS or PHO. This variant altered splicing, resulting in an exon 7-truncated transcript lacking 16 bases. No normal transcript was detected in the patient's stomach or colon tissue. qPCR also showed significantly decreased SLCO2A1 transcription compared to HC. CONCLUSION: A previously unreported variant caused defective SLCO2A1 splicing and reduced mRNA levels in a patient with CEAS and PHO. This research enhances understanding of CEAS and PHO pathophysiology and aids genetic counseling and diagnosis.

Complete genome sequence of a haloalkaliphilic heterotrophic bacterium Halomonas salifodinae IM328.

The genome of a halophilic bacterium Halomonas salifodinae IM328 was completely sequenced in order to offer convenience for the research such as the synthesis of compatible solutes. The genome contains a circular chromosome which was sequenced by PacBio system.

Spontaneous Remission of Epileptic Seizures Following Norovirus Infection in a Patient With DNM1 Encephalopathy.

Epileptic seizures can be worsened by infections; however, they sometimes disappear or decrease after an acute viral infection, although this is rare. We report the spontaneous remission of epileptic seizures following norovirus-induced viral gastroenteritis in a boy with DNM1 encephalopathy. He had clonic seizures daily from the age of two months and developed epileptic spasms at 14 months of age; he was admitted to the hospital at this time. A physical examination revealed hypotonia, strabismus, tongue protrusion with drooping, and widely spaced teeth. Although brain magnetic resonance imaging was unremarkable, electroencephalography revealed frequent occipital spikes. Three days after admission, the patient developed frequent diarrhea without a fever. A rapid immunochromatographic test of norovirus in a stool sample was positive. Immediately after the appearance of diarrhea, the epileptic seizures disappeared. Currently, at the age of five years, the patient has a profound psychomotor developmental delay; he has no verbal expression and is unable to walk. He has experienced involuntary movements of the myoclonus since 10 months of age. Whole-exome sequencing of the patient's DNA revealed the presence of a heterozygous de novo variant of DNM1: c.709C>T (p.Arg237Trp). Although the findings from our patient suggest that underlying neural network abnormalities were ameliorated by immunological mechanisms as a result of the viral infection, further research is needed to clarify the mechanisms behind this spontaneous remission of seizures.

To explore the mechanism of acupoint application in the treatment of primary dysmenorrhea by 16S rDNA sequencing and metabolomics.

Graphene-based warm uterus acupoint paste (GWUAP) is an emerging non-drug alternative therapy for the treatment of primary dysmenorrhea (PD), but the underlying mechanism is still unclear. SD female rats were randomly divided into control group, model group and treatment group to explore the mechanism of GWUAP in the treatment of PD. Combined with 16S rDNA and fecal metabolomics, the diversity of microbiota and metabolites in each group was comprehensively evaluated. In this study, GWUAP reduced the torsion score of PD model rats, improved the pathological morphology of uterine tissue, reduced the pathological damage score of uterine tissue, and reversed the expression levels of inflammatory factors, pain factors and sex hormones. The 16 S rDNA sequencing of fecal samples showed that the abundance of Lactobacillus in the intestinal flora of the model group decreased and the abundance of Romboutsia increased, while the abundance of Lactobacillus in the intestinal flora of the treatment group increased and the abundance of Romboutsia decreased, which improved the imbalance of flora diversity in PD rats. In addition, 32 metabolites related to therapeutic effects were identified by metabolomics of fecal samples. Moreover, there is a close correlation between fecal microbiota and metabolites. Therefore, the mechanism of GWUAP in the treatment of PD remains to be further studied.

Identification of Novel Hb Guiyang [HBA2: c.151C > A α2 50 (CE8) His- Asn] and Phenotype- Genotype Correlation of Abnormal Hemoglobins in Guizhou, Southwest China.

Purpose: To analyze the composition of abnormal hemoglobin and the relationship between genotype and phenotype by screening abnormal hemoglobin in a subpopulation of Guizhou, China. Patients and Methods: Routine blood evaluation, capillary electrophoresis of hemoglobin, and mutation of α - and ß - thalassemia genes were evaluated in 19,976 individuals for thalassemia screening in Guizhou. Sanger sequencing of HBA1, HBA2 and HBB genes was performed in samples with abnormal bands or unexplained increases of normal bands. The types of abnormal hemoglobin were obtained by sequence analysis. Results: Abnormal hemoglobin was detected in 84 individuals (detection rate, 0.42%). Ten types each of α and ß globin chain variants were detected, including most commonly Hb E, Hb New York and Hb Port Phillip. In this study, the abnormal Hb Mizuho was identified for the first time in a Chinese population, and a novel abnormal hemoglobin Hb Guiyang (HBA2: c.151C > A) was detected for the first time. Except for Hb Mizuho, other abnormal hemoglobin heterozygotes without thalassemia or iron deficiency had no significant hematological changes. Conclusion: This study enriched the molecular epidemiological data of abnormal hemoglobin in Guizhou, China and provided reference data for genetic counseling and prenatal diagnosis of abnormal hemoglobin.

[Malignant transformation of Ollier disease-related multiple glioma with IDH1 p.R132C mutation].

A 21-year-old man who was diagnosed with Ollier disease at the age of 1 year developed incidental multiple gliomas at the age of 15 years. Subsequently, the multiple gliomas enlarged and the patient underwent three surgical removals. Genetic analysis revealed the IDH1 p.R132C mutation in the gliomas, and histopathology showed malignant transformation. Despite multimodality treatment, the gliomas could not be controlled, and the patient died at the age of 23 years. Ollier disease is a rare disease with IDH1/2 mutations and is often associated with gliomas. However, there are very few reports on genetic analysis of IDH1/2 mutations and long-term follow-up in Ollier disease-related gliomas. Genetic analysis of IDH mutations may contribute to the elucidation of its pathogenesis. The cross-departmental collaboration is required for long-term follow-up of Ollier disease-related gliomas.

Integrated approach to generate artificial samples with low tumor fraction for somatic variant calling benchmarking.

BACKGROUND: High-throughput sequencing (HTS) has become the gold standard approach for variant analysis in cancer research. However, somatic variants may occur at low fractions due to contamination from normal cells or tumor heterogeneity; this poses a significant challenge for standard HTS analysis pipelines. The problem is exacerbated in scenarios with minimal tumor DNA, such as circulating tumor DNA in plasma. Assessing sensitivity and detection of HTS approaches in such cases is paramount, but time-consuming and expensive: specialized experimental protocols and a sufficient quantity of samples are required for processing and analysis. To overcome these limitations, we propose a new computational approach specifically designed for the generation of artificial datasets suitable for this task, simulating ultra-deep targeted sequencing data with low-fraction variants and demonstrating their effectiveness in benchmarking low-fraction variant calling. RESULTS: Our approach enables the generation of artificial raw reads that mimic real data without relying on pre-existing data by using NEAT, a fine-grained read simulator that generates artificial datasets using models learned from multiple different datasets. Then, it incorporates low-fraction variants to simulate somatic mutations in samples with minimal tumor DNA content. To prove the suitability of the created artificial datasets for low-fraction variant calling benchmarking, we used them as ground truth to evaluate the performance of widely-used variant calling algorithms: they allowed us to define tuned parameter values of major variant callers, considerably improving their detection of very low-fraction variants. CONCLUSIONS: Our findings highlight both the pivotal role of our approach in creating adequate artificial datasets with low tumor fraction, facilitating rapid prototyping and benchmarking of algorithms for such dataset type, as well as the important need of advancing low-fraction variant calling techniques.

Sequence confirmation of synthetic DNA exceeding 100 nucleotides using restriction enzyme mediated digestion combined with high-resolution tandem mass spectrometry.

Oligonucleotides have emerged as important therapeutic options for inherited diseases. In recent years, RNA therapeutics, especially mRNA, have been pushed to the market. Analytical methods for these molecules have been published extensively in the last few years. Notably, mass spectrometry has proven as a state-of-the-art quality control method. For RNA based therapeutics, numerous methods are available, while DNA therapeutics lack of suitable MS-based methods when it comes to molecules exceeding approximately 60 nucleotides. We present a method which combines the use of common restriction enzymes and short enzyme-directing oligonucleotides to generate DNA digestion products with the advantages of high-resolution tandem mass spectrometry. The instrumentation includes ion pair reverse phase chromatography coupled to a time-of-flight mass spectrometer with a collision induced dissociation (CID) for sequence analysis. Utilizing this approach, we increased the sequence coverage from 23.3% for a direct CID-MS/MS experiment of a 100 nucleotide DNA molecule to 100% sequence coverage using the restriction enzyme mediated approach presented in this work. This approach is suitable for research and development and quality control purposes in a regulated environment, which makes it a versatile tool for drug development.

High-throughput nanopore DNA sequencing of large insert fosmid clones directly from bacterial colonies.

Fosmids and cosmids are vectors frequently used in functional metagenomic studies. With a large insert capacity (around 30 kb) they can encode dozens of cloned genes or in some cases, entire biochemical pathways. Fosmids with cloned inserts can be transferred to heterologous hosts and propagated to enable screening for new enzymes and metabolites. After screening, fosmids from clones with an activity of interest must be de novo sequenced, a critical step toward the identification of the gene(s) of interest. In this work, we present a new approach for rapid and high-throughput fosmid sequencing directly from Escherichia coli colonies without liquid culturing or fosmid purification. Our sample preparation involves fosmid amplification with phi29 polymerase and then direct nanopore sequencing using the Oxford Nanopore Technologies system. We also present a bioinformatics pipeline termed "phiXXer" that facilitates both de novo read assembly and vector trimming to generate a linear sequence of the fosmid insert. Finally, we demonstrate the accurate sequencing of 96 fosmids in a single run and validate the method using two fosmid libraries that contain cloned large insert (~30-40 kb) genomic or metagenomic DNA.IMPORTANCELarge-insert clone (fosmids or cosmids) sequencing is challenging and arguably the most limiting step of functional metagenomic screening workflows. Our study establishes a new method for high-throughput nanopore sequencing of fosmid clones directly from lysed Escherichia coli cells. It also describes a companion bioinformatic pipeline that enables de novo assembly of fosmid DNA insert sequences. The devised method widens the potential of functional metagenomic screening by providing a simple, high-throughput approach to fosmid clone sequencing that dramatically speeds the pace of discovery.

HATCHet2: clone- and haplotype-specific copy number inference from bulk tumor sequencing data.

Bulk DNA sequencing of multiple samples from the same tumor is becoming common, yet most methods to infer copy-number aberrations (CNAs) from this data analyze individual samples independently. We introduce HATCHet2, an algorithm to identify haplotype- and clone-specific CNAs simultaneously from multiple bulk samples. HATCHet2 extends the earlier HATCHet method by improving identification of focal CNAs and introducing a novel statistic, the minor haplotype B-allele frequency (mhBAF), that enables identification of mirrored-subclonal CNAs. We demonstrate HATCHet2's improved accuracy using simulations and a single-cell sequencing dataset. HATCHet2 analysis of 10 prostate cancer patients reveals previously unreported mirrored-subclonal CNAs affecting cancer genes.

Viral genome sequencing methods: benefits and pitfalls of current approaches.

Whole genome sequencing of viruses provides high-resolution molecular insights, enhancing our understanding of viral genome function and phylogeny. Beyond fundamental research, viral sequencing is increasingly vital for pathogen surveillance, epidemiology, and clinical applications. As sequencing methods rapidly evolve, the diversity of viral genomics applications and catalogued genomes continues to expand. Advances in long-read, single molecule, real-time sequencing methodologies present opportunities to sequence contiguous, haplotype resolved viral genomes in a range of research and applied settings. Here we present an overview of nucleic acid sequencing methods and their applications in studying viral genomes. We emphasise the advantages of different viral sequencing approaches, with a particular focus on the benefits of third-generation sequencing technologies in elucidating viral evolution, transmission networks, and pathogenesis.

Modified lentiviral globin gene therapy for pediatric ß0/ß0 transfusion-dependent ß-thalassemia: A single-center, single-arm pilot trial.

ß0/ß0 thalassemia is the most severe type of transfusion-dependent ß-thalassemia (TDT) and is still a challenge facing lentiviral gene therapy. Here, we report the interim analysis of a single-center, single-arm pilot trial (NCT05015920) evaluating the safety and efficacy of a ß-globin expression-optimized and insulator-engineered lentivirus-modified cell product (BD211) in ß0/ß0 TDT. Two female children were enrolled, infused with BD211, and followed up for an average of 25.5 months. Engraftment of genetically modified hematopoietic stem and progenitor cells was successful and sustained in both patients. No unexpected safety issues occurred during conditioning or after infusion. Both patients achieved transfusion independence for over 22 months. The treatment extended the lifespan of red blood cells by over 42 days. Single-cell DNA/RNA-sequencing analysis of the dynamic changes of gene-modified cells, transgene expression, and oncogene activation showed no notable adverse effects. Optimized lentiviral gene therapy may safely and effectively treat all ß-thalassemia.

Development of an in vitro biofilm model of the human supra-gingival microbiome for Oral microbiome transplantation.

The high prevalence of dental caries and periodontal disease place a significant burden on society, both socially and economically. Recent advances in genomic technologies have linked both diseases to shifts in the oral microbiota - a community of >700 bacterial species that live within the mouth. The development of oral microbiome transplantation draws on the success of fecal microbiome transplantation for the treatment of gut pathologies associated with disease. Many current in vitro oral biofilm models have been developed but do not fully capture the complexity of the oral microbiome which is required for successful OMT. To address this, we developed an in vitro biofilm system that maintained an oral microbiome with 252 species on average over 14 days. Six human plaque samples were grown in 3D printed flow cells on hydroxyapatite discs using artificial saliva medium (ASM). Biofilm composition and growth were monitored by high throughput sequencing and confocal microscopy/SEM, respectively. While a significant drop in bacterial diversity occurred, up to 291 species were maintained in some flow cells over 14 days with 70% viability grown with ASM. This novel in vitro biofilm model represents a marked improvement on existing oral biofilm systems and provides new opportunities to develop oral microbiome transplant therapies.

Residual disinfectant effectively suppresses Legionella species in drinking water distribution systems supplied by surface water in Minnesota, USA.

Seven public water systems in Minnesota, USA were analyzed from one to five times over a two-year period to assess temporal changes in the concentrations of total bacteria, Legionella spp., and Legionella pneumophila from source (i.e., raw water) through the water treatment process to the end water user. Bacterial biomass was collected by filtering large volumes of raw water (12 to 425 L, median: 38 L) or finished and tap water (27 to 1205 L, median: 448 L) using ultrafiltration membrane modules. Quantitative PCR (qPCR) was then used to enumerate all bacteria (16S rRNA gene fragments), all Legionella spp. (ssrA), and Legionella pneumophila (mip). Total coliforms, Escherichia coli, and L. pneumophila also were quantified in the water samples via cultivation. Median concentrations of total bacteria and Legionella spp. (ssrA) in raw water (8.5 and 4.3 log copies/L, respectively) decreased by about 2 log units during water treatment. The concentration of Legionella spp. (ssrA) in water collected from distribution systems inversely correlated with the total chlorine concentration for chloraminated systems significantly (p = 0.03). Although only 8 samples were collected from drinking water distribution systems using free chlorine as a residual disinfectant, these samples had significantly lower concentrations of Legionella spp. (ssrA) than samples collected from the chloraminated systems (p = 5 × 10-4). There was considerable incongruity between the results obtained via cultivation-independent (qPCR) and cultivation-dependent assays. Numerous samples were positive for L. pneumophila via cultivation, none of which tested positive for L. pneumophilia (mip) via qPCR. Conversely, a single sample tested positive for L. pneumophilia (mip) via qPCR, but this sample tested negative for L. pneumophilia via cultivation. Overall, the results suggest that conventional treatment is effective at reducing, but not eliminating, Legionella spp. from surface water supplies and that residual disinfection is effective at suppressing these organisms within drinking water distribution systems.

Precision Basecalling of Single DNA Nucleotide from Overlapped Transmission Readouts with Machine Learning Aided Solid-State Nanogap.

DNA sequencing with the quantum tunneling technique heralds a paradigm shift in genetic analysis, promising rapid and accurate identification for diverging applications ranging from personalized medicine to security issues. However, the widespread distribution of molecular conductance, conduction orbital alignment for resonant transport, and decoding crisscrossing conductance signals of isomorphic nucleotides have been persistent experimental hurdles for swift and precise identification. Herein, we have reported a machine learning (ML)-driven quantum tunneling study with solid-state model nanogap to determine nucleotides at single-base resolution. The optimized ML basecaller has demonstrated a high predictive basecalling accuracy of all four nucleotides from seven distinct data pools, each containing complex transmission readouts of their different dynamic conformations. ML classification of quaternary, ternary, and binary nucleotide combinations is also performed with high precision, sensitivity, and F1 score. ML explainability unravels the evidence of how extracted normalized features within overlapped nucleotide signals contribute to classification improvement. Moreover, electronic fingerprints, conductance sensitivity, and current readout analysis of nucleotides have promised practical applicability with significant sensitivity and distinguishability. Through this ML approach, our study pushes the boundaries of quantum sequencing by highlighting the effectiveness of single nucleotide basecalling with promising implications for advancing genomics and molecular diagnostics.

Genetic Screening-Emerging Issues.

In many countries, some form of genetic screening is offered to all or part of the population, either in the form of well-organized screening programs or in a less formalized way. Screening can be offered at different phases of life, such as preconception, prenatal, neonatal and later in life. Screening should only be offered if the advantages outweigh the disadvantages. Technical innovations in testing and treatment are driving changes in the field of prenatal and neonatal screening, where many jurisdictions have organized population-based screening programs. As a result, a greater number and wider range of conditions are being added to the programs, which can benefit couples' reproductive autonomy (preconception and prenatal screening) and improve early diagnosis to prevent irreversible health damage in children (neonatal screening) and in adults (cancer and cascade screening). While many developments in screening are technology-driven, citizens may also express a demand for innovation in screening, as was the case with non-invasive prenatal testing. Relatively new emerging issues for genetic screening, especially if testing is performed using DNA sequencing, relate to organization, data storage and interpretation, benefit-harm ratio and distributive justice, information provision and follow-up, all connected to acceptability in current healthcare systems.