The whole blood transcriptional regulation landscape in 465 COVID-19 infected samples from Japan COVID-19 Task Force

「コロナ制圧タスクフォース」COVID-19患者由来の血液細胞における遺伝子発現の網羅的解析 --重症度に応じた遺伝子発現の変化には、ヒトゲノム配列の個人差が影響する--. 京都大学プレスリリース. 2022-08-23.Coronavirus disease 2019 (COVID-19) is a recently-emerged infectious disease that has caused millions of deaths, where comprehensive understanding of disease mechanisms is still unestablished. In particular, studies of gene expression dynamics and regulation landscape in COVID-19 infected individuals are limited. Here, we report on a thorough analysis of whole blood RNA-seq data from 465 genotyped samples from the Japan COVID-19 Task Force, including 359 severe and 106 non-severe COVID-19 cases. We discover 1169 putative causal expression quantitative trait loci (eQTLs) including 34 possible colocalizations with biobank fine-mapping results of hematopoietic traits in a Japanese population, 1549 putative causal splice QTLs (sQTLs; e.g. two independent sQTLs at TOR1AIP1), as well as biologically interpretable trans-eQTL examples (e.g., REST and STING1), all fine-mapped at single variant resolution. We perform differential gene expression analysis to elucidate 198 genes with increased expression in severe COVID-19 cases and enriched for innate immune-related functions. Finally, we evaluate the limited but non-zero effect of COVID-19 phenotype on eQTL discovery, and highlight the presence of COVID-19 severity-interaction eQTLs (ieQTLs; e.g., CLEC4C and MYBL2). Our study provides a comprehensive catalog of whole blood regulatory variants in Japanese, as well as a reference for transcriptional landscapes in response to COVID-19 infection

DOCK2 is involved in the host genetics and biology of severe COVID-19

「コロナ制圧タスクフォース」COVID-19疾患感受性遺伝子DOCK2の重症化機序を解明 --アジア最大のバイオレポジトリーでCOVID-19の治療標的を発見--. 京都大学プレスリリース. 2022-08-10.Identifying the host genetic factors underlying severe COVID-19 is an emerging challenge. Here we conducted a genome-wide association study (GWAS) involving 2, 393 cases of COVID-19 in a cohort of Japanese individuals collected during the initial waves of the pandemic, with 3, 289 unaffected controls. We identified a variant on chromosome 5 at 5q35 (rs60200309-A), close to the dedicator of cytokinesis 2 gene (DOCK2), which was associated with severe COVID-19 in patients less than 65 years of age. This risk allele was prevalent in East Asian individuals but rare in Europeans, highlighting the value of genome-wide association studies in non-European populations. RNA-sequencing analysis of 473 bulk peripheral blood samples identified decreased expression of DOCK2 associated with the risk allele in these younger patients. DOCK2 expression was suppressed in patients with severe cases of COVID-19. Single-cell RNA-sequencing analysis (n = 61 individuals) identified cell-type-specific downregulation of DOCK2 and a COVID-19-specific decreasing effect of the risk allele on DOCK2 expression in non-classical monocytes. Immunohistochemistry of lung specimens from patients with severe COVID-19 pneumonia showed suppressed DOCK2 expression. Moreover, inhibition of DOCK2 function with CPYPP increased the severity of pneumonia in a Syrian hamster model of SARS-CoV-2 infection, characterized by weight loss, lung oedema, enhanced viral loads, impaired macrophage recruitment and dysregulated type I interferon responses. We conclude that DOCK2 has an important role in the host immune response to SARS-CoV-2 infection and the development of severe COVID-19, and could be further explored as a potential biomarker and/or therapeutic target

Retrospective evaluation of whole exome and genome mutation calls in 746 cancer samples

Funder: NCI U24CA211006Abstract: The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) curated consensus somatic mutation calls using whole exome sequencing (WES) and whole genome sequencing (WGS), respectively. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, which aggregated whole genome sequencing data from 2,658 cancers across 38 tumour types, we compare WES and WGS side-by-side from 746 TCGA samples, finding that ~80% of mutations overlap in covered exonic regions. We estimate that low variant allele fraction (VAF < 15%) and clonal heterogeneity contribute up to 68% of private WGS mutations and 71% of private WES mutations. We observe that ~30% of private WGS mutations trace to mutations identified by a single variant caller in WES consensus efforts. WGS captures both ~50% more variation in exonic regions and un-observed mutations in loci with variable GC-content. Together, our analysis highlights technological divergences between two reproducible somatic variant detection efforts

Fluctuation of Global Gene Expression by Endogenous miRNA Response to the Introduction of an Exogenous miRNA

Most of the intracellular endogenous microRNAs (endo-miRNAs) are considered to be saturated in Argonaute (Ago) proteins in the RNA-induced silencing complexes (RISCs). When exogenous miRNAs (exo-miRNAs) are introduced into cells, endo-miRNAs in the RISC may be replaced with exo-miRNAs or exo-miRNAs, and endo-miRNAs might also compete for the position in the newly synthesized RISC with each other. This would lead to the fluctuation of global gene expression not only by repression of exo-miRNA target gene expression, but also by the increase of the endo-miRNA target gene expression. In the present study, we quantified the changes in the expression levels of target genes of exo-miRNA and endo-miRNA in the cells transfected with fifteen different exo-miRNAs by microarray experiments. Different exo-miRNAs increased ratios of expression levels of target genes of a given endo-miRNA to different extents, suggesting that the replacement efficiencies might differ according to the exo-miRNA types. However, the increased ratios in the expression levels of each endo-miRNA target genes by the transfection of any particular exo-miRNA were mostly equivalent, suggesting that the endo-miRNAs present in the RISC might be replaced with excessive exo-miRNAs at similar levels, probably because they exist in single-stranded forms in the RISC

High-entropy intermetallics serve an ultrastable single-atom Pt for propane dehydrogenation

Propane dehydrogenation (PDH) has been a promising propylene production process that can compensate for the increasing global demand for propylene. However, Pt-based catalysts with high stability at ≥600°C have barely been reported because the catalysts typically result in short catalyst life owing to side reactions and coke formation. Herein, we report a new class of heterogeneous catalysts using high-entropy intermetallics (HEIs). Pt–Pt ensembles, which cause side reactions, are entirely diluted by the component inert metals in PtGe-type HEI; thereby, unfavorable side reactions are drastically inhibited. The resultant HEI: (PtCoCu)(GeGeSn)/Ca–SiO2 exhibited an outstandingly high catalytic stability, even at 600°C (kd−1 = τ = 4146 h = 173 d), and almost no deactivation of the catalyst was observed two months for the first time

Flow Characteristics of Slurry with Rare-Earth Rich Mud under Deep Seabed around Minamitorishima

Flow characteristics of slurry with rare-earth rich mud are investigated to examine lifting systems for mining rare-earth elements from deep seabed. Twelve samples, extracted from different locations and depths from seabed around the Minamitorishima with several concentration of rare-earth elements are mixed up with sea water to make slurry with specified volume concentration of mud 1.0% to 10.0%. Cone-plate type rotary viscometers are used to examine relations between shear stress and shear rate of slurry in each volume concentration of mud. As a result, it was found that shape of graph; the shear stress in the vertical axis and the shear rate in the horizontal axis was concave down and increasing during whole range of the shear rate. The shear stress increased gradually as the shear rate increased in case of low volume concentration of mud up to 3.0%. On the other hand, the shear stress changed significantly at small shear rate, and then gradually increased, then lineally in the end in case of larger volume concentration. In addition, similar characteristics under the same volume concentration, even though tested samples were extracted from different locations, depths and concentration of rare-earth elements. Further, three types of fluid model; the Power low model, the Bingham-Papanastasiou model and the Herchel-Bulkley-Papanastasiou model were fitted on the data using the least square techniques, then compared with each other. The last two models, i.e., the Bingham-Papanastasiou model and the Herchel-Bulkley-Papanastasiou model corrects deviations from the data when using “original” the Bingham model and the Herchel-Bulkley model, especially in the range of small shear rate under high volume concentration of mud. The Herchel-Bulkley-Papanastasiou model was the most appropriate model within the three models. Furthermore, correlation equations for parameters of the HerchelBulkley-Papanastasiou model were derived related to volume concentration of mud

IntroMap: a signal analysis based method for the detection of genomic introgressions

Background: Breeding programs often rely on marker-assisted tests or variant calling of next generation sequence (NGS) data to identify regions of genomic introgression arising from the hybridization of two plant species. In this paper we present IntroMap, a bioinformatics pipeline for the screening of candidate plants through the application of signal processing techniques to NGS data, using alignment to a reference genome sequence (annotation is not required) that shares homology with the recurrent parental cultivar, and without the need for de novo assembly of the read data or variant calling. Results: We show the accurate identification of introgressed genomic regions using both in silico simulated genomes, and a hybridized cultivar data set using our pipeline. Additionally we show, through targeted marker-based assays, validation of the IntroMap predicted regions for the hybrid cultivar. Conclusions: This approach can be used to automate the screening of large populations, reducing the time and labor required, and can improve the accuracy of the detection of introgressed regions in comparison to a marker-based approach. In contrast to other approaches that generally rely upon a variant calling step, our method achieves accurate identification of introgressed regions without variant calling, relying solely upon alignment