Immune response stability to the SARS-CoV-2 mRNA vaccine booster is influenced by differential splicing of HLA genes.

Many molecular mechanisms that lead to the host antibody response to COVID-19 vaccines remain largely unknown. In this study, we used serum antibody detection combined with whole blood RNA-based transcriptome analysis to investigate variability in vaccine response in healthy recipients of a booster (third) dose schedule of the mRNA BNT162b2 vaccine against COVID-19. The cohort was divided into two groups: (1) low-stable individuals, with antibody concentration anti-SARS-CoV IgG S1 below 0.4 percentile at 180 days after boosting vaccination; and (2) high-stable individuals, with antibody values greater than 0.6 percentile of the range in the same period (median 9525 [185-80,000] AU/mL). Differential gene expression, expressed single nucleotide variants and insertions/deletions, differential splicing events, and allelic imbalance were explored to broaden our understanding of the immune response sustenance. Our analysis revealed a differential expression of genes with immunological functions in individuals with low antibody titers, compared to those with higher antibody titers, underscoring the fundamental importance of the innate immune response for boosting immunity. Our findings also provide new insights into the determinants of the immune response variability to the SARS-CoV-2 mRNA vaccine booster, highlighting the significance of differential splicing regulatory mechanisms, mainly concerning HLA alleles, in delineating vaccine immunogenicity.

Differential haplotype expression in class I MHC genes during SARS-CoV-2 infection of human lung cell lines.

Introduction: Cell entry of SARS-CoV-2 causes genome-wide disruption of the transcriptional profiles of genes and biological pathways involved in the pathogenesis of COVID-19. Expression allelic imbalance is characterized by a deviation from the Mendelian expected 1:1 expression ratio and is an important source of allele-specific heterogeneity. Expression allelic imbalance can be measured by allele-specific expression analysis (ASE) across heterozygous informative expressed single nucleotide variants (eSNVs). ASE reflects many regulatory biological phenomena that can be assessed by combining genome and transcriptome information. ASE contributes to the interindividual variability associated with the disease. We aim to estimate the transcriptome-wide impact of SARS-CoV-2 infection by analyzing eSNVs. Methods: We compared ASE profiles in the human lung cell lines Calu-3, A459, and H522 before and after infection with SARS-CoV-2 using RNA-Seq experiments. Results: We identified 34 differential ASE (DASE) sites in 13 genes (HLA-A, HLA-B, HLA-C, BRD2, EHD2, GFM2, GSPT1, HAVCR1, MAT2A, NQO2, SUPT6H, TNFRSF11A, UMPS), all of which are enriched in protein binding functions and play a role in COVID-19. Most DASE sites were assigned to the MHC class I locus and were predominantly upregulated upon infection. DASE sites in the MHC class I locus also occur in iPSC-derived airway epithelium basal cells infected with SARS-CoV-2. Using an RNA-Seq haplotype reconstruction approach, we found DASE sites and adjacent eSNVs in phase (i.e., predicted on the same DNA strand), demonstrating differential haplotype expression upon infection. We found a bias towards the expression of the HLA alleles with a higher binding affinity to SARS-CoV-2 epitopes. Discussion: Independent of gene expression compensation, SARS-CoV-2 infection of human lung cell lines induces transcriptional allelic switching at the MHC loci. This suggests a response mechanism to SARS-CoV-2 infection that swaps HLA alleles with poor epitope binding affinity, an expectation supported by publicly available proteome data.

Associations Between the Purinergic Receptor P2X7 and Leprosy Disease.

Leprosy is an infectious disease still highly prevalent in Brazil, having been detected around 27,863 new cases in 2019. Exposure to Mycobacterium leprae may not be sufficient to trigger the disease, which seems to be influenced by host immunogenetics to determine resistance or susceptibility. The purinergic receptor P2X7 plays a crucial role in immunity, inflammation, neurological function, bone homeostasis, and neoplasia and is associated with several infectious and non-infectious diseases. Here, we first compare the P2RX7 expression in RNA-seq experiments from 16 leprosy cases and 16 healthy controls to establish the magnitude of allele-specific expression for single-nucleotide polymorphisms of the gene P2RX7 and to determine the level of gene expression in healthy and diseased skin. In addition, we also evaluated the association of two P2RX7 single-nucleotide polymorphisms (c.1513A>C/rs3751143 and c.1068A>G/rs1718119) with leprosy risk. The expression of P2RX7 was found significantly upregulated at macrophage cells from leprosy patients compared with healthy controls, mainly in macrophages from lepromatous patients. Significant risk for leprosy disease was associated with loss function of rs3751143 homozygous mutant CC [CC vs. AA: p = 0.001; odds ratio (OR) = 1.676, 95% CI = 1.251-2.247] but not with heterozygous AC (AC vs. AA: p = 0.001; OR = 1.429, 95% CI = 1.260-1.621). Contrary, the polymorphic A allele from the gain function of rs1718119 was associated with protection for the development of leprosy, as observed in the dominant model (AA + AG × GG p = 0.0028; OR = 0.03516; CI = 0.1801-0.6864). So, our results suggest that the functional P2X7 purinergic receptor may exert a key role in the Mycobacterium death inside macrophages and inflammatory response, which is necessary to control the disease.

Exome-Wide Search for Genes Associated With Central Nervous System Inflammatory Demyelinating Diseases Following CHIKV Infection: The Tip of the Iceberg.

Chikungunya virus (CHIKV) is a re-emergent arbovirus that causes a disease characterized primarily by fever, rash and severe persistent polyarthralgia, although <1% of cases develop severe neurological manifestations such as inflammatory demyelinating diseases (IDD) of the central nervous system (CNS) like acute disseminated encephalomyelitis (ADEM) and extensive transverse myelitis. Genetic factors associated with host response and disease severity are still poorly understood. In this study, we performed whole-exome sequencing (WES) to identify HLA alleles, genes and cellular pathways associated with CNS IDD clinical phenotype outcomes following CHIKV infection. The cohort includes 345 patients of which 160 were confirmed for CHIKV. Six cases presented neurological manifestation mimetizing CNS IDD. WES data analysis was performed for 12 patients, including the CNS IDD cases and 6 CHIKV patients without any neurological manifestation. We identified 29 candidate genes harboring rare, pathogenic, or probably pathogenic variants in all exomes analyzed. HLA alleles were also determined and patients who developed CNS IDD shared a common signature with diseases such as Multiple sclerosis (MS) and Neuromyelitis Optica Spectrum Disorders (NMOSD). When these genes were included in Gene Ontology analyses, pathways associated with CNS IDD syndromes were retrieved, suggesting that CHIKV-induced CNS outcomesmay share a genetic background with other neurological disorders. To our knowledge, this study was the first genome-wide investigation of genetic risk factors for CNS phenotypes in CHIKV infection. Our data suggest that HLA-DRB1 alleles associated with demyelinating diseases may also confer risk of CNS IDD outcomes in patients with CHIKV infection.