A blood RNA transcriptome signature for COVID-19.

BACKGROUND: COVID-19 is a respiratory viral infection with unique features including a more chronic course and systemic disease manifestations including multiple organ involvement; and there are differences in disease severity between ethnic groups. The immunological basis for disease has not been fully characterised. Analysis of whole-blood RNA expression may provide valuable information on disease pathogenesis. METHODS: We studied 45 patients with confirmed COVID-19 infection within 10 days from onset of illness and a control group of 19 asymptomatic healthy volunteers with no known exposure to COVID-19 in the previous 14 days. Relevant demographic and clinical information was collected and a blood sample was drawn from all participants for whole-blood RNA sequencing. We evaluated differentially-expressed genes in COVID-19 patients (log2 fold change ≥ 1 versus healthy controls; false-discovery rate < 0.05) and associated protein pathways and compared these to published whole-blood signatures for respiratory syncytial virus (RSV) and influenza. We developed a disease score reflecting the overall magnitude of expression of internally-validated genes and assessed the relationship between the disease score and clinical disease parameters. RESULTS: We found 135 differentially-expressed genes in the patients with COVID-19 (median age 35 years; 82% male; 36% Chinese, 53% South Asian ethnicity). Of the 117 induced genes, 14 were found in datasets from RSV and 40 from influenza; 95 genes were unique to COVID-19. Protein pathways were mostly generic responses to viral infections, including apoptosis by P53-associated pathway, but also included some unique pathways such as viral carcinogenesis. There were no major qualitative differences in pathways between ethnic groups. The composite gene-expression score was correlated with the time from onset of symptoms and nasal swab qPCR CT values (both p < 0.01) but was not related to participant age, gender, ethnicity or the presence or absence of chest X-ray abnormalities (all p > 0.05). CONCLUSIONS: The whole-blood transcriptome of COVID-19 has overall similarity with other respiratory infections but there are some unique pathways that merit further exploration to determine clinical relevance. The approach to a disease score may be of value, but needs further validation in a population with a greater range of disease severity.

Distinct genetic architectures and environmental factors associate with host response to the γ2-herpesvirus infections.

Kaposi's sarcoma-associated herpesvirus (KSHV) and Epstein-Barr Virus (EBV) establish life-long infections and are associated with malignancies. Striking geographic variation in incidence and the fact that virus alone is insufficient to cause disease, suggests other co-factors are involved. Here we present epidemiological analysis and genome-wide association study (GWAS) in 4365 individuals from an African population cohort, to assess the influence of host genetic and non-genetic factors on virus antibody responses. EBV/KSHV co-infection (OR = 5.71(1.58-7.12)), HIV positivity (OR = 2.22(1.32-3.73)) and living in a more rural area (OR = 1.38(1.01-1.89)) are strongly associated with immunogenicity. GWAS reveals associations with KSHV antibody response in the HLA-B/C region (p = 6.64 × 10-09). For EBV, associations are identified for VCA (rs71542439, p = 1.15 × 10-12). Human leucocyte antigen (HLA) and trans-ancestry fine-mapping substantiate that distinct variants in HLA-DQA1 (p = 5.24 × 10-44) are driving associations for EBNA-1 in Africa. This study highlights complex interactions between KSHV and EBV, in addition to distinct genetic architectures resulting in important differences in pathogenesis and transmission.

A blood RNA transcript signature for TB exposure in household contacts.

BACKGROUND: Current tools for diagnosing latent TB infection (LTBI) detect immunological memory of past exposure but are unable to determine whether exposure is recent. We sought to identify a whole-blood transcriptome signature of recent TB exposure. METHODS: We studied household contacts of TB patients; healthy volunteers without recent history of TB exposure; and patients with active TB. We performed whole-blood RNA sequencing (in all), an interferon gamma release assay (IGRA; in contacts and healthy controls) and PET/MRI lung scans (in contacts only). We evaluated differentially-expressed genes in household contacts (log2 fold change ≥1 versus healthy controls; false-discovery rate < 0.05); compared these to differentially-expressed genes seen in the active TB group; and assessed the association of a composite gene expression score to independent exposure/treatment/immunological variables. RESULTS: There were 186 differentially-expressed genes in household contacts (n = 26, age 22-66, 46% male) compared with healthy controls (n = 5, age 29-38, 100% male). Of these genes, 141 (76%) were also differentially expressed in active TB (n = 14, age 27-69, 71% male). The exposure signature included genes from inflammatory response, type I interferon signalling and neutrophil-mediated immunity pathways; and genes such as BATF2 and SCARF1 known to be associated with incipient TB. The composite gene-expression score was higher in IGRA-positive contacts (P = 0.04) but not related to time from exposure, isoniazid prophylaxis, or abnormalities on PET/MRI (all P > 0.19). CONCLUSIONS: Transcriptomics can detect TB exposure and, with further development, may be an approach of value for epidemiological research and targeting public health interventions.

Sustained Low-Level Transmission of Zika and Chikungunya Viruses after Emergence in the Fiji Islands.

Zika and chikungunya viruses were first detected in Fiji in 2015. Examining surveillance and phylogenetic and serologic data, we found evidence of low-level transmission of Zika and chikungunya viruses during 2013-2017, in contrast to the major outbreaks caused by closely related virus strains in other Pacific Island countries.

HLA-C variants associated with amino acid substitutions in the peptide binding groove influence susceptibility to Kawasaki disease.

Kawasaki disease (KD) is a pediatric vasculitis caused by an unknown trigger in genetically susceptible children. The incidence varies widely across genetically diverse populations. Several associations with HLA Class I alleles have been reported in single cohort studies. Using a genetic approach, from the nine single nucleotide variants (SNVs) associated with KD susceptibility in children of European descent, we identified SNVs near the HLA-C (rs6906846) and HLA-B genes (rs2254556) whose association was replicated in a Japanese descent cohort (rs6906846 p = 0.01, rs2254556 p = 0.005). The risk allele (A at rs6906846) was also associated with HLA-C*07:02 and HLA-C*04:01 in both US multi-ethnic and Japanese cohorts and HLA-C*12:02 only in the Japanese cohort. The risk A-allele was associated with eight non-conservative amino acid substitutions (amino acid positions); Asp or Ser (9), Arg (14), Ala (49), Ala (73), Ala (90), Arg (97), Phe or Ser (99), and Phe or Ser (116) in the HLA-C peptide binding groove that binds peptides for presentation to cytotoxic T cells (CTL). This raises the possibility of increased affinity to a "KD peptide" that contributes to the vasculitis of KD in genetically susceptible children.

Identification of regulatory variants associated with genetic susceptibility to meningococcal disease.

Non-coding genetic variants play an important role in driving susceptibility to complex diseases but their characterization remains challenging. Here, we employed a novel approach to interrogate the genetic risk of such polymorphisms in a more systematic way by targeting specific regulatory regions relevant for the phenotype studied. We applied this method to meningococcal disease susceptibility, using the DNA binding pattern of RELA - a NF-kB subunit, master regulator of the response to infection - under bacterial stimuli in nasopharyngeal epithelial cells. We designed a custom panel to cover these RELA binding sites and used it for targeted sequencing in cases and controls. Variant calling and association analysis were performed followed by validation of candidate polymorphisms by genotyping in three independent cohorts. We identified two new polymorphisms, rs4823231 and rs11913168, showing signs of association with meningococcal disease susceptibility. In addition, using our genomic data as well as publicly available resources, we found evidences for these SNPs to have potential regulatory effects on ATXN10 and LIF genes respectively. The variants and related candidate genes are relevant for infectious diseases and may have important contribution for meningococcal disease pathology. Finally, we described a novel genetic association approach that could be applied to other phenotypes.

Primary macrophages and J774 cells respond differently to infection with Mycobacterium tuberculosis.

Macrophages play an essential role in the early immune response to Mycobacterium tuberculosis and are the cell type preferentially infected in vivo. Primary macrophages and macrophage-like cell lines are commonly used as infection models, although the physiological relevance of cell lines, particularly for host-pathogen interaction studies, is debatable. Here we use high-throughput RNA-sequencing to analyse transcriptome dynamics of two macrophage models in response to M. tuberculosis infection. Specifically, we study the early response of bone marrow-derived mouse macrophages and cell line J774 to infection with live and γ-irradiated (killed) M. tuberculosis. We show that infection with live bacilli specifically alters the expression of host genes such as Rsad2, Ifit1/2/3 and Rig-I, whose potential roles in resistance to M. tuberculosis infection have not yet been investigated. In addition, the response of primary macrophages is faster and more intense than that of J774 cells in terms of number of differentially expressed genes and magnitude of induction/repression. Our results point to potentially novel processes leading to immune containment early during M. tuberculosis infection, and support the idea that important differences exist between primary macrophages and cell lines, which should be taken into account when choosing a macrophage model to study host-pathogen interactions.

Inhibition of megakaryocyte development in the bone marrow underlies dengue virus-induced thrombocytopenia in humanized mice.

A characteristic clinical feature of dengue virus infection is thrombocytopenia, though its underlying mechanism is not definitively determined. By adoptive transfer of human CD34(+) fetal liver cells into immunodeficient mice, we have constructed humanized mice with significant levels of human platelets, monocytes/macrophages, and hepatocytes. Infection of these mice with both lab-adapted and clinical strains of dengue virus induces characteristic human hematological changes, including transient leukopenia and thrombocytopenia. We show that the specific depletion of human platelets is not mediated by antibodies in the periphery or reduced production of human thrombopoietin in the liver but reduction of human megakaryocytes and megakaryocyte progenitors in the bone marrow of the infected mice. These findings identify inhibition of platelet production in the bone marrow as a key mechanism underlying dengue-induced thrombocytopenia and suggest the utility of the improved humanized mouse model in studying dengue virus infection and pathogenesis in a human cell context.

Genetic variants of MICB and PLCE1 and associations with non-severe dengue.

BACKGROUND: A recent genome-wide association study (GWAS) identified susceptibility loci for dengue shock syndrome (DSS) at MICB rs3132468 and PLCE1 rs3740360. The aim of this study was to define the extent to which MICB (rs3132468) and PLCE1 (rs3740360) were associated with less severe clinical phenotypes of pediatric and adult dengue. METHODS: 3961 laboratory-confirmed dengue cases and 5968 controls were genotyped at MICB rs3132468 and PLCE1 rs3740360. Per-allele odds ratios (OR) with 95% confidence intervals (CI) were calculated for each patient cohort. Pooled analyses were performed for adults and paediatrics respectively using a fixed effects model. RESULTS: Pooled analysis of the paediatric and adult cohorts indicated a significant association between MICB rs3132468 and dengue cases without shock (OR  =  1.15; 95%CI: 1.07 - 1.24; P  =  0.0012). Similarly, pooled analysis of pediatric and adult cohorts indicated a significant association between dengue cases without shock and PLCE1 rs3740360 (OR  =  0.92; 95%CI: 0.85 - 0.99; P  =  0.018). We also note significant association between both SNPs (OR  =  1.48; P  =  0.0075 for MICB rs3132468 and OR  =  0.75, P  =  0.041 for PLCE1 rs3740360) and dengue in infants. DISCUSSION: This study confirms that the MICB rs3132468 and PLCE1 rs3740360 risk genotypes are not only associated with DSS, but are also associated with less severe clinical phenotypes of dengue, as well as with dengue in infants. These findings have implications for our understanding of dengue pathogenesis.

Whole-genome reconstruction and mutational signatures in gastric cancer.

BACKGROUND: Gastric cancer is the second highest cause of global cancer mortality. To explore the complete repertoire of somatic alterations in gastric cancer, we combined massively parallel short read and DNA paired-end tag sequencing to present the first whole-genome analysis of two gastric adenocarcinomas, one with chromosomal instability and the other with microsatellite instability. RESULTS: Integrative analysis and de novo assemblies revealed the architecture of a wild-type KRAS amplification, a common driver event in gastric cancer. We discovered three distinct mutational signatures in gastric cancer--against a genome-wide backdrop of oxidative and microsatellite instability-related mutational signatures, we identified the first exome-specific mutational signature. Further characterization of the impact of these signatures by combining sequencing data from 40 complete gastric cancer exomes and targeted screening of an additional 94 independent gastric tumors uncovered ACVR2A, RPL22 and LMAN1 as recurrently mutated genes in microsatellite instability-positive gastric cancer and PAPPA as a recurrently mutated gene in TP53 wild-type gastric cancer. CONCLUSIONS: These results highlight how whole-genome cancer sequencing can uncover information relevant to tissue-specific carcinogenesis that would otherwise be missed from exome-sequencing data.

Disruption of vascular homeostasis in patients with Kawasaki disease: involvement of vascular endothelial growth factor and angiopoietins.

OBJECTIVE: In Kawasaki disease (KD), a pediatric vasculitis of medium-sized arteries, the coronary arteries are most commonly affected. Angiopoietins and vascular endothelial growth factor (VEGF) play an important role in maintaining vascular homeostasis. Recently, we identified ANGPT1 and VEGFA as susceptibility loci for KD. This study was undertaken to fine-map these associations and to gain further insight into their role in this vasculitis of unknown etiology to further the search for improved diagnostic and therapeutic options. METHODS: A total of 292 single-nucleotide polymorphisms (SNPs) located in VEGF and ANGPT and their receptors were genotyped in 574 families, including 462 trios. For replication, 123 cases and 171 controls were genotyped. RESULTS: A significant association with KD susceptibility was observed with 5 SNPs in the ANGPT1 gene (most significantly associated SNP +265037 C>T; Pcombined=2.3×10(-7) ) and 2 SNPs in VEGFA (most significantly associated SNP rs3025039; Pcombined=2.5×10(-4) ). Both ANGPT1 +265037 C>T and VEGFA rs3025039 are located in 3' regulatory regions at putative transcription factor binding sites. We observed significantly down-regulated transcript levels of angiopoietin 1 (Ang-1) in patients with acute KD compared to patients with convalescent KD. In patients with acute KD, high serum protein levels of VEGF and Ang-2 were observed compared to patients with convalescent KD and to both controls with and controls without fever. Immunohistochemistry demonstrated VEGF and angiopoietin expression in the coronary artery wall in autopsy tissue. CONCLUSION: Our data support the hypothesis that dysregulation of VEGF and angiopoietins contributes to the disruption of vascular homeostasis in KD.

Genome-wide association study identifies susceptibility loci for dengue shock syndrome at MICB and PLCE1.

Hypovolemic shock (dengue shock syndrome (DSS)) is the most common life-threatening complication of dengue. We conducted a genome-wide association study of 2,008 pediatric cases treated for DSS and 2,018 controls from Vietnam. Replication of the most significantly associated markers was carried out in an independent Vietnamese sample of 1,737 cases and 2,934 controls. SNPs at two loci showed genome-wide significant association with DSS. We identified a susceptibility locus at MICB (major histocompatibility complex (MHC) class I polypeptide-related sequence B), which was within the broad MHC region on chromosome 6 but outside the class I and class II HLA loci (rs3132468, P(meta) = 4.41 × 10(-11), per-allele odds ratio (OR) = 1.34 (95% confidence interval: 1.23-1.46)). We identified associated variants within PLCE1 (phospholipase C, epsilon 1) on chromosome 10 (rs3765524, P(meta) = 3.08 × 10(-10), per-allele OR = 0.80 (95% confidence interval: 0.75-0.86)). We identify two loci associated with susceptibility to DSS in people with dengue, suggesting possible mechanisms for this severe complication of dengue.

Genome wide expression profiling reveals suppression of host defence responses during colonisation by Neisseria meningitides but not N. lactamica.

Both Neisseria meningitidis and the closely related bacterium Neisseria lactamica colonise human nasopharyngeal mucosal surface, but only N. meningitidis invades the bloodstream to cause potentially life-threatening meningitis and septicaemia. We have hypothesised that the two neisserial species differentially modulate host respiratory epithelial cell gene expression reflecting their disease potential. Confluent monolayers of 16HBE14 human bronchial epithelial cells were exposed to live and/or dead N. meningitidis (including capsule and pili mutants) and N. lactamica, and their transcriptomes were compared using whole genome microarrays. Changes in expression of selected genes were subsequently validated using Q-RT-PCR and ELISAs. Live N. meningitidis and N. lactamica induced genes involved in host energy production processes suggesting that both bacterial species utilise host resources. N. meningitidis infection was associated with down-regulation of host defence genes. N. lactamica, relative to N. meningitidis, initiates up-regulation of proinflammatory genes. Bacterial secreted proteins alone induced some of the changes observed. The results suggest N. meningitidis and N. lactamica differentially regulate host respiratory epithelial cell gene expression through colonisation and/or protein secretion, and that this may contribute to subsequent clinical outcomes associated with these bacteria.

Transforming growth factor-beta signaling pathway in patients with Kawasaki disease.

BACKGROUND: Transforming growth factor (TGF)-ß is a multifunctional peptide that is important in T-cell activation and cardiovascular remodeling, both of which are important features of Kawasaki disease (KD). We postulated that variation in TGF-ß signaling might be important in KD susceptibility and disease outcome. METHODS AND RESULTS: We investigated genetic variation in 15 genes belonging to the TGF-ß pathway in a total of 771 KD subjects of mainly European descent from the United States, the United Kingdom, Australia, and the Netherlands. We analyzed transcript abundance patterns using microarray and reverse transcriptase-polymerase chain reaction for these same genes, and measured TGF-ß2 protein levels in plasma. Genetic variants in TGFB2, TGFBR2, and SMAD3 and their haplotypes were consistently and reproducibly associated with KD susceptibility, coronary artery aneurysm formation, aortic root dilatation, and intravenous immunoglobulin treatment response in different cohorts. A SMAD3 haplotype associated with KD susceptibility replicated in 2 independent cohorts and an intronic single nucleotide polymorphism in a separate haplotype block was also strongly associated (A/G, rs4776338) (P=0.000022; odds ratio, 1.50; 95% confidence interval, 1.25 to 1.81). Pathway analysis using all 15 genes further confirmed the importance of the TGF-ß pathway in KD pathogenesis. Whole-blood transcript abundance for these genes and TGF-ß2 plasma protein levels changed dynamically over the course of the illness. CONCLUSIONS: These studies suggest that genetic variation in the TGF-ß pathway influences KD susceptibility, disease outcome, and response to therapy, and that aortic root and coronary artery Z scores can be used for phenotype/genotype analyses. Analysis of transcript abundance and protein levels further support the importance of this pathway in KD pathogenesis.

Dengue virus regulates type I interferon signalling in a strain-dependent manner in human cell lines.

Outbreaks of dengue disease are constant threats to tropical and subtropical populations but range widely in severity, from mild to haemorrhagic fevers, for reasons that are still elusive. We investigated the interferon (IFN) response in infected human cell lines A549 and HepG2, using two strains (NGC and TSV01) of dengue serotype 2 (DEN2) and found that the two viruses exhibited a marked difference in inducing type I IFN response. While TSV01 infection led to activation of type I antiviral genes such as EIF2AK2 (PKR), OAS, ADAR and MX, these responses were absent in NGC-infected cells. Biochemical analysis revealed that NGC but not TSV01 suppressed STAT-1 and STAT-2 activation in response to type I IFN (alpha and beta). However, these two strains did not differ in their response to type II IFN (gamma). Although unable to suppress IFN signalling, TSV01 infection caused a weaker IFN-beta induction compared with NGC, suggesting an alternative mechanism of innate immune escape. We extended our study to clinical isolates of various serotypes and found that while MY10245 (DEN2) and MY22713 (DEN4) could suppress the IFN response in a similar fashion to NGC, three other strains of dengue [EDEN167 (DEN1), MY02569 (DEN1) and MY10340 (DEN2)] were unable to suppress the IFN response, suggesting that this difference is strain-dependent but not serotype-specific. Our report indicates the existence of a strain-specific virulence factor that may impact on disease severity.

Whole-genome expression profiling reveals that inhibition of host innate immune response pathways by Ebola virus can be reversed by a single amino acid change in the VP35 protein.

Ebola hemorrhagic fever is a rapidly progressing acute febrile illness characterized by high virus replication, severe immunosuppression, and case fatalities of ca. 80%. Inhibition of phosphorylation of interferon regulatory factor 3 (IRF-3) by the Ebola VP35 protein may block the host innate immune response and play an important role in the severity of disease. We used two precisely defined reverse genetics-generated Ebola viruses to investigate global host cell responses resulting from the inhibition of IRF-3 phosphorylation. The two viruses encoded either wild-type (WT) VP35 protein (recEbo-VP35/WT) or VP35 with an arginine (R)-to-alanine (A) amino acid substitution at position 312 (recEbo-VP35/R312A) within a previously defined IRF-3 inhibitory domain. When sucrose-gradient purified virus was used for infection, host cell whole-genome expression profiling revealed striking differences in human liver cell responses to these viruses differing by a single amino acid. The inhibition of host innate immune responses by WT Ebola virus was so potent that little difference in interferon and antiviral gene expression could be discerned between cells infected with purified WT, inactivated virus, or mock-infected cells. However, infection with recEbo-VP35/R312A virus resulted in a strong innate immune response including increased expression of MDA-5, RIG-I, RANTES, MCP-1, ISG-15, ISG-54, ISG-56, ISG-60, STAT1, IRF-9, OAS, and Mx1. The clear gene expression differences were obscured if unpurified virus stocks were used to initiate infection, presumably due to soluble factors present in virus-infected cell supernatant preparations. Ebola virus VP35 protein clearly plays a pivotal role in the potent inhibition of the host innate immune responses, and the present study indicates that VP35 has a wider effect on host cell responses than previously shown. The ability to eliminate this inhibitory effect with a single amino acid change in VP35 demonstrates the critical role this protein must play in the severe aspects this highly fatal disease.

Host gene expression profiling of dengue virus infection in cell lines and patients.

BACKGROUND: Despite the seriousness of dengue-related disease, with an estimated 50-100 million cases of dengue fever and 250,000-500,000 cases of dengue hemorrhagic fever/dengue shock syndrome each year, a clear understanding of dengue pathogenesis remains elusive. Because of the lack of a disease model in animals and the complex immune interaction in dengue infection, the study of host response and immunopathogenesis is difficult. The development of genomics technology, microarray and high throughput quantitative PCR have allowed researchers to study gene expression changes on a much broader scale. We therefore used this approach to investigate the host response in dengue virus-infected cell lines and in patients developing dengue fever. METHODOLOGY/PRINCIPAL FINDINGS: Using microarray and high throughput quantitative PCR method to monitor the host response to dengue viral replication in cell line infection models and in dengue patient blood samples, we identified differentially expressed genes along three major pathways; NF-kappaB initiated immune responses, type I interferon (IFN) and the ubiquitin proteasome pathway. Among the most highly upregulated genes were the chemokines IP-10 and I-TAC, both ligands of the CXCR3 receptor. Increased expression of IP-10 and I-TAC in the peripheral blood of ten patients at the early onset of fever was confirmed by ELISA. A highly upregulated gene in the IFN pathway, viperin, was overexpressed in A549 cells resulting in a significant reduction in viral replication. The upregulation of genes in the ubiquitin-proteasome pathway prompted the testing of proteasome inhibitors MG-132 and ALLN, both of which reduced viral replication. CONCLUSION/SIGNIFICANCE: Unbiased gene expression analysis has identified new host genes associated with dengue infection, which we have validated in functional studies. We showed that some parts of the host response can be used as potential biomarkers for the disease while others can be used to control dengue viral replication, thus representing viable targets for drug therapy.