Published papers

Splicing quantitative trait loci (sQTLs) are one of the major causal mechanisms in genome-wide association study (GWAS) loci, but their role in disease pathogenesis is poorly understood. One reason is the complexity of alternative splicing events producing many unknown isoforms. Here, we propose two approaches, namely integration and selection, for this complexity by focusing on protein-structure of isoforms. First, we integrate isoforms with the same coding sequence (CDS) and identify 369-601 integrated-isoform ratio QTLs (i2-rQTLs), which altered protein-structure, in six immune subsets. Second, we select CDS incomplete isoforms annotated in GENCODE and identify 175-337 isoform-ratio QTL (i-rQTL). By comprehensive long-read capture RNA-sequencing among these incomplete isoforms, we reveal 29 full-length isoforms with unannotated CDSs associated with GWAS traits. Furthermore, we show that disease-causal sQTL genes can be identified by evaluating their trans-eQTL effects. Our approaches highlight the understudied role of protein-altering sQTLs and are broadly applicable to other tissues and diseases.

In genetic testing of Mendelian diseases, it is a bioinformatics challenge to effectively prioritize disease-causing candidate genes listed from massively parallel sequencing. Tissue specificity of the gene expression levels may give a clue because it may reflect tissue-specific disease manifestation. However, considering poor correlations between mRNA and protein expression in some genes, it is not clear whether transcriptomics- or proteomics-based tissue specificity should be used to prioritize candidate genes. Therefore, we compared the efficiency of tissue-specific scores (TS scores) obtained from transcriptome and proteome data in prioritizing candidate genes for whole exome sequencing (WES) analysis of Mendelian disease patients. We show that both Protein and RNA TS scores are useful in prioritizing candidate genes in WES analysis, although diseases like coagulopathies get more benefit from Protein TS score. This study may provide useful evidence in developing new methods to effectively identify novel disease-causing genes.

Pathophysiological analysis and drug discovery targeting human diseases require disease models that suitably recapitulate patient pathology. Disease-specific human induced pluripotent stem cells (hiPSCs) differentiated into affected cell types can potentially recapitulate disease pathology more accurately than existing disease models. Such successful modeling of muscular diseases requires efficient differentiation of hiPSCs into skeletal muscles. hiPSCs transduced with doxycycline-inducible MYOD1 (MYOD1-hiPSCs) have been widely used; however, they require time- and labor-consuming clonal selection, and clonal variations must be overcome. Moreover, their functionality should be carefully examined. Here, we demonstrated that bulk MYOD1-hiPSCs established with puromycin selection rather than G418 selection showed rapid and highly efficient differentiation. Interestingly, bulk MYOD1-hiPSCs exhibited average differentiation properties of clonally established MYOD1-hiPSCs, suggesting that it is possible to minimize clonal variations. Moreover, disease-specific hiPSCs of spinal bulbar muscular atrophy (SBMA) could be efficiently differentiated via this method into skeletal muscle that showed disease phenotypes, suggesting the applicability of this method for disease analysis. Finally, three-dimensional muscle tissues were fabricated from bulk MYOD1-hiPSCs, which exhibited contractile force upon electrical stimulation, indicating their functionality. Thus, our bulk differentiation requires less time and labor than existing methods, efficiently generates contractible skeletal muscles, and may facilitate the generation of muscular disease models.

Pathogenic variants in WDR45 on chromosome Xp11 cause neurodegenerative disorder beta-propeller protein-associated neurodegeneration (BPAN). Currently, there is no effective therapy for BPAN. Here we report a 17-year-old female patient with BPAN and show that antisense oligonucleotide (ASO) was effective in vitro. The patient had developmental delay and later showed extrapyramidal signs since the age of 15 years. MRI findings showed iron deposition in the globus pallidus and substantia nigra on T2 MRI. Whole genome sequencing and RNA sequencing revealed generation of pseudoexon due to inclusion of intronic sequences triggered by an intronic variant that is remote from the exon–intron junction: WDR45 (OMIM #300526) chrX(GRCh37):g.48935143G > C, (NM_007075.4:c.235 + 159C > G). We recapitulated the exonization of intron sequences by a mini-gene assay and further sought antisense oligonucleotide that induce pseudoexon skipping using our recently developed, a dual fluorescent splicing reporter system that encodes two fluorescent proteins, mCherry, a transfection marker designed to facilitate evaluation of exon skipping and split eGFP, a splicing reaction marker. The results showed that the 24-base ASO was the strongest inducer of pseudoexon skipping. Our data presented here have provided supportive evidence for in vivo preclinical studies.

The decision of whether cells are activated or not is controlled through dynamic intracellular molecular networks. However, the low population of cells during the transition state of activation renders the analysis of the transcriptome of this state technically challenging. To address this issue, we have developed the Time-Dependent Cell-State Selection (TDCSS) technique, which employs live-cell imaging of secretion activity to detect an index of the transition state, followed by the simultaneous recovery of indexed cells for subsequent transcriptome analysis. In this study, we used the TDCSS technique to investigate the transition state of group 2 innate lymphoid cells (ILC2s) activation, which is indexed by the onset of interleukin (IL)-13 secretion. The TDCSS approach allowed us to identify time-dependent genes, including transiently induced genes (TIGs). Our findings of IL4 and MIR155HG as TIGs have shown a regulatory function in ILC2s activation.

Accumulating evidence indicates that long intergenic non-coding RNAs (lincRNAs) show more tissue-specific expression patterns than protein-coding genes (PCGs). However, although lincRNAs are subject to canonical transcriptional regulation like PCGs, the molecular basis for the specificity of their expression patterns remains unclear. Here, using expression data and coordinates of topologically associating domains (TADs) in human tissues, we show that lincRNA loci are significantly enriched in the more internal region of TADs compared to PCGs and that lincRNAs within TADs have higher tissue specificity than those outside TADs. Based on these, we propose an analytical framework to interpret transcriptional status using lincRNA as an indicator. We applied it to hypertrophic cardiomyopathy data and found disease-specific transcriptional regulation: ectopic expression of keratin at the TAD level and derepression of myocyte differentiation-related genes by E2F1 with down-regulation of LINC00881. Our results provide understanding of the function and regulation of lincRNAs according to genomic structure.

Sonic hedgehog signaling molecule (SHH) is a key molecule in the cilia-mediated signaling pathway and a critical morphogen in embryogenesis. The association between loss-of-function variants of SHH and holoprosencephaly is well established. In mice experiments, reduced or increased signaling of SHH have been shown to be associated with narrowing or excessive expansion of the facial midline, respectively. Herein, we report two unrelated patients with de novo truncating variants of SHH presenting with hypertelorism rather than hypotelorism. The first patient was a 13-year-old girl. Her facial features included hypertelorism, strabismus, telecanthus, malocclusion, frontal bossing, and wide widow's peak. She had borderline developmental delay and agenesis of the corpus callosum. She had a nonsense variant of SHH: Chr7(GRCh38):g.155802987C > T, NM_000193.4:c.1302G > A, p.(Trp434*). The second patient was a 25-year-old girl. Her facial features included hypertelorism and wide widow's peak. She had developmental delay and agenesis of the corpus callosum. She had a frameshift variant of SHH: Chr7(GRCh38):g.155803072_155803074delCGGinsT, NM_000193.4:c.1215_1217delCCGinsA, p.(Asp405Glufs*92). The hypertelorism phenotype contrasts sharply with the prototypical hypotelorism-holoprosencephaly phenotype associated with loss-of-function of SHH. We concluded that a subset of truncating variants of SHH could be associated with hypertelorism rather than hypotelorism.

＜文献概要＞Point　●コロナウイルスの1つである新型コロナウイルス(SARS-CoV-2)はプラス鎖のRNAウイルスであり,5'側からメッセンジャーRNA(mRNA)としてそのまま翻訳されるのと並行して,3'側から不連続に転写されたマイナス鎖を鋳型として複数のmRNAを合成することで,多種のタンパク質を作り出す複雑な仕組みをもつ.●ゲノム解析技術の発展により,SARS-CoV-2であれば検体採取から数時間以内に全ウイルスゲノム配列を決定できるようになった.そしてその配列データはデータベース化され,全世界で共有されている.●ウイルスは自身のゲノムの変異により感染力や病原性を変化させながら生き残りを図っている.そのウイルスゲノムを解析することで,ウイルス株の伝播動態や変異株の発生などを即時的に把握し,院内感染対策や行政施策,創薬研究などに活用することができる.

次世代シーケンサーが登場したことで、極少数の家系検体からでも低コストで短時間に疾患原因変異の同定が可能となり、近年の疾患変異解析研究の世界には大革新が起こった。しかし、次世代シーケンサーから出力される膨大なビッグデータから、真の1つの疾患原因遺伝子変異の特定に至る過程の解析は決して容易とは言えない。本稿では、遺伝性が疑われる疾患に対し、次世代シーケンサーを用いて疾患原因遺伝子変異を同定する際に用いられる手法とその解析方法および課題を現時点での最新情報を交えて解説する。(著者抄録)