Construction of a trio-based structural variation panel utilizing activated T lymphocytes and long-read sequencing technology

Comprehensive and deep evaluation of structural variation detection pipelines with third-generation sequencing data

BACKGROUND

Structural variation (SV) detection methods using third-generation sequencing data are widely employed, yet accurately detecting SVs remains challenging. Different methods often yield inconsistent results for certain SV types, complicating tool selection and revealing biases in detection.

RESULTS

This study comprehensively evaluates 53 SV detection pipelines using simulated and real data from PacBio (CLR: Continuous Long Read, CCS: Circular Consensus Sequencing) and Nanopore (ONT) platforms. We assess their performance in detecting various sizes and types of SVs, breakpoint biases, and genotyping accuracy with various sequencing depths. Notably, pipelines such as Minimap2-cuteSV2, NGMLR-SVIM, PBMM2-pbsv, Winnowmap-Sniffles2, and Winnowmap-SVision exhibit comparatively higher recall and precision. Our findings also show that combining multiple pipelines with the same aligner, like pbmm2 or winnowmap, can significantly enhance performance. The individual pipelines' detailed ranking and performance metrics can be viewed in a dynamic table: http://pmglab.top/SVPipelinesRanking .

CONCLUSIONS

This study comprehensively characterizes the strengths and weaknesses of numerous pipelines, providing valuable insights that can improve SV detection in third-generation sequencing data and inform SV annotation and function prediction.

Whole blood transcriptome analysis for age- and gender-specific gene expression profiling in Japanese individuals

Whole blood transcriptome analysis is a valuable approachin medical research, primarily due to the ease of sample collection and the richness of the information obtained. Since the expression profile of individual genes in the analysis is influenced by medical traits and demographic attributes such as age and gender, there has been a growing demand for a comprehensive database for blood transcriptome analysis. Here, we performed whole blood RNA sequencing (RNA-seq) analysis on 576 participants stratified by age (20-30s and 60-70s) and gender from cohorts of the Tohoku Medical Megabank (TMM). A part of female segment included pregnant women. We did not exclude the globin gene family in our RNA-seq study, which enabled us to identify instances of hereditary persistence of fetal hemoglobin based on the HBG1 and HBG2 expression information. Comparing stratified populations allowed us to identify groups of genes associated with age-related changes and gender differences. We also found that the immune response status, particularly measured by neutrophil-to-lymphocyte ratio (NLR), strongly influences the diversity of individual gene expression profiles in whole blood transcriptome analysis. This stratification has resulted in a data set that will be highly beneficial for future whole blood transcriptome analysis in the Japanese population.

A homozygous structural variant of RPGRIP1 is frequently associated with achromatopsia in Japanese patients with IRD

Purpose

Achromatopsia (ACHM) is an early-onset cone dysfunction caused by 5 genes with cone-specific functions (CNGA3, CNGB3, GNAT2, PDE6C, and PDE6H) and by ATF6, a transcription factor with ubiquitous expression. To improve the relatively low variant detection ratio in these genes in a cohort of exome-sequenced Japanese patients with inherited retinal diseases (IRD), we performed genome sequencing to detect structural variants and intronic variants in patients with ACHM.

Methods

Genome sequencing of 10 ACHM pedigrees was performed after exome sequencing. Structural, non-coding, and coding variants were filtered based on segregation between the affected and unaffected in each pedigree. Variant frequency and predicted damage scores were considered in identifying pathogenic variants.

Results

A homozygous deletion involving exon 18 of RPGRIP1 was detected in 5 of 10 ACHM probands, and variant inheritance from each parent was confirmed. This deletion was relatively frequent (minor allele frequency = 0.0023) in the Japanese population but was only homozygous in patients with ACHM among the 199 Japanese IRD probands analyzed by the same genome sequencing pipeline.

Conclusion

The deletion involving exon 18 of RPGRIP1 is a prevalent cause of ACHM in Japanese patients and contributes to the wide spectrum of RPGRIP1-associated IRD phenotypes, from Leber congenital amaurosis to ACHM.

The landscape of genomic structural variation in Indigenous Australians

Indigenous Australians harbour rich and unique genomic diversity. However, Aboriginal and Torres Strait Islander ancestries are historically under-represented in genomics research and almost completely missing from reference datasets1-3. Addressing this representation gap is critical, both to advance our understanding of global human genomic diversity and as a prerequisite for ensuring equitable outcomes in genomic medicine. Here we apply population-scale whole-genome long-read sequencing4 to profile genomic structural variation across four remote Indigenous communities. We uncover an abundance of large insertion-deletion variants (20-49 bp; n = 136,797), structural variants (50  b-50 kb; n = 159,912) and regions of variable copy number (>50 kb; n = 156). The majority of variants are composed of tandem repeat or interspersed mobile element sequences (up to 90%) and have not been previously annotated (up to 62%). A large fraction of structural variants appear to be exclusive to Indigenous Australians (12% lower-bound estimate) and most of these are found in only a single community, underscoring the need for broad and deep sampling to achieve a comprehensive catalogue of genomic structural variation across the Australian continent. Finally, we explore short tandem repeats throughout the genome to characterize allelic diversity at 50 known disease loci5, uncover hundreds of novel repeat expansion sites within protein-coding genes, and identify unique patterns of diversity and constraint among short tandem repeat sequences. Our study sheds new light on the dimensions and dynamics of genomic structural variation within and beyond Australia.

Tohoku Medical Megabank Brain Magnetic Resonance Imaging Study: Rationale, Design, and Background

The Tohoku Medical Megabank Brain Magnetic Resonance Imaging Study (TMM Brain MRI Study) was established to collect multimodal information through neuroimaging and neuropsychological assessments to evaluate the cognitive function and mental health of residents who experienced the Great East Japan Earthquake (GEJE) and associated tsunami. The study also aimed to promote advances in personalized healthcare and medicine related to mental health and cognitive function among the general population. We recruited participants for the first (baseline) survey starting in July 2014, enrolling individuals who were participating in either the TMM Community-Based Cohort Study (TMM CommCohort Study) or the TMM Birth and Three-Generation Cohort Study (TMM BirThree Cohort Study). We collected multiple magnetic resonance imaging (MRI) sequences, including 3D T1-weighted sequences, magnetic resonance angiography (MRA), diffusion tensor imaging (DTI), pseudo-continuous arterial spin labeling (pCASL), and three-dimensional fluid-attenuated inversion recovery (FLAIR) sequences. To assess neuropsychological status, we used both questionnaire- and interview-based rating scales. The former assessments included the Tri-axial Coping Scale, Impact of Event Scale in Japanese, Profile of Mood States, and 15-item Depression, Anxiety, and Stress Scale, whereas the latter assessments included the Mini-Mental State Examination, Japanese version. A total of 12,164 individuals were recruited for the first (baseline) survey, including those unable to complete all assessments. In parallel, we returned the MRI results to the participants and subsequently shared the MRI data through the TMM Biobank. At present, the second (first follow-up) survey of the study started in October 2019 is underway. In this study, we established a large and comprehensive database that included robust neuroimaging data as well as psychological and cognitive assessment data. In combination with genomic and omics data already contained in the TMM Biobank database, these data could provide new insights into the relationships of pathological processes with neuropsychological disorders, including age-related cognitive impairment.

Accelerated nanopore basecalling with SLOW5 data format

MOTIVATION

Nanopore sequencing is emerging as a key pillar in the genomic technology landscape but computational constraints limiting its scalability remain to be overcome. The translation of raw current signal data into DNA or RNA sequence reads, known as 'basecalling', is a major friction in any nanopore sequencing workflow. Here, we exploit the advantages of the recently developed signal data format 'SLOW5' to streamline and accelerate nanopore basecalling on high-performance computing (HPC) and cloud environments.

RESULTS

SLOW5 permits highly efficient sequential data access, eliminating a potential analysis bottleneck. To take advantage of this, we introduce Buttery-eel, an open-source wrapper for Oxford Nanopore's Guppy basecaller that enables SLOW5 data access, resulting in performance improvements that are essential for scalable, affordable basecalling.

AVAILABILITY AND IMPLEMENTATION

Buttery-eel is available at https://github.com/Psy-Fer/buttery-eel.