Esr1 but Not CYP19A1 Overexpression in Mammary Epithelial Cells during Reproductive Senescence Induces Pregnancy-Like Proliferative Mammary Disease Responsive to Anti-Hormonals.

Molecular-level analyses of breast carcinogenesis benefit from vivo disease models. Estrogen receptor 1 (Esr1) and cytochrome P450 family 19 subfamily A member 1 (CYP19A1) overexpression targeted to mammary epithelial cells in genetically engineered mouse models induces largely similar rates of proliferative mammary disease in prereproductive senescent mice. Herein, with natural reproductive senescence, Esr1 overexpression compared with CYP19A1 overexpression resulted in significantly higher rates of preneoplasia and cancer. Before reproductive senescence, Esr1, but not CYP19A1, overexpressing mice are tamoxifen resistant. However, during reproductive senescence, Esr1 mice exhibited responsiveness. Both Esr1 and CYP19A1 are responsive to letrozole before and after reproductive senescence. Gene Set Enrichment Analyses of RNA-sequencing data sets showed that higher disease rates in Esr1 mice were accompanied by significantly higher expression of cell proliferation genes, including members of prognostic platforms for women with early-stage hormone receptor-positive disease. Tamoxifen and letrozole exposure induced down-regulation of these genes and resolved differences between the two models. Both Esr1 and CYP19A1 overexpression induced abnormal developmental patterns of pregnancy-like gene expression. This resolved with progression through reproductive senescence in CYP19A1 mice, but was more persistent in Esr1 mice, resolving only with tamoxifen and letrozole exposure. In summary, genetically engineered mouse models of Esr1 and CYP19A1 overexpression revealed a diversion of disease processes resulting from the two distinct molecular pathophysiological mammary gland-targeted intrusions into estrogen signaling during reproductive senescence.

Overexpression of Estrogen Receptor α in Mammary Glands of Aging Mice Is Associated with a Proliferative Risk Signature and Generation of Estrogen Receptor α-Positive Mammary Adenocarcinomas.

Age is a risk factor for human estrogen receptor-positive breast cancer, with highest prevalence following menopause. While transcriptome risk profiling is available for human breast cancers, it is not yet developed for prognostication for primary or secondary breast cancer development utilizing at-risk breast tissue. Both estrogen receptor α (ER) and aromatase overexpression have been linked to human breast cancer. Herein, conditional genetically engineered mouse models of estrogen receptor 1 (Esr1) and cytochrome P450 family 19 subfamily A member 1 (CYP19A1) were used to show that induction of Esr1 overexpression just before or with reproductive senescence and maintained through age 30 months resulted in significantly higher prevalence of estrogen receptor-positive adenocarcinomas than CYP19A1 overexpression. All adenocarcinomas tested showed high percentages of ER+ cells. Mammary cancer development was preceded by a persistent proliferative transcriptome risk signature initiated within 1 week of transgene induction that showed parallels to the Prosigna/Prediction Analysis of Microarray 50 human prognostic signature for early-stage human ER+ breast cancer. CYP19A1 mice also developed ER+ mammary cancers, but histology was more divided between adenocarcinoma and adenosquamous, with one ER- adenocarcinoma. Results demonstrate that, like humans, generation of ER+ adenocarcinoma in mice was facilitated by aging mice past the age of reproductive senescence. Esr1 overexpression was associated with a proliferative estrogen pathway-linked signature that preceded appearance of ER+ mammary adenocarcinomas.

Mouse Mammary Gland Whole Mount Density Assessment across Different Morphologies Using a Bifurcated Program for Image Processing.

Mammographic density is associated with increased breast cancer risk. Conventional visual assessment of murine mouse models does not include quantified total density analysis. A bifurcated method was sufficient to obtain relative density scores on a broad range of two-dimensional whole mount images that contained both normal and abnormal findings. Image processing techniques, including a ridge operator and a gaussian denoising method, were used to isolate background away from mammary epithelium and use mean pixel intensity to represent mammary density on genetically engineered mouse models for breast cancer in mice 4 to 29 months of age. The bifurcated method allowed for application of an optimal image processing approach for the structural elements present in the whole mount images. Gaussian denoising was the optimal approach when more dense lobular growth and tertiary branching dominate and a ridge operator when epithelial growth was more sparse and secondary branching was the more dominant structural feature. The two processing approaches were combined in a single experimental flow program using an initial image density measurement as the decision point between the two approaches. Higher density was associated with lobular growth, tertiary branching, fibrotic stroma, and presence of cancer. The significance of the study is development of a readily accessible program for digital assessment of mammary gland whole mount density across a range of mammary gland morphologies.

The PPARγ agonist efatutazone delays invasive progression and induces differentiation of ductal carcinoma in situ.

PURPOSE: Ductal carcinoma in situ (DCIS) is a pre-invasive lesion of the breast considered a precursor of invasive ductal carcinoma. This study aimed to determine whether activated PPARγ acts as a tumor suppressor in human DCIS progression. METHODS: We utilized the high-affinity PPARγ agonist, efatutazone, to activate endogenous PPARγ in a well-defined model for the progression of basal (triple negative) DCIS, MCFDCIS cells, cultured under 2D and 3D conditions. We studied the effects of activated PPARγ on DCIS progression in MCFDCIS xenograft and C3(1)/Tag transgenic mice treated with 30 mg/kg of efatutazone. RESULTS: In vitro, efatutazone did not alter the MCFDCIS cell proliferation but induced phenotypic and gene expression changes, indicating that activated PPARγ is able to differentiate MCFDCIS cells into more luminal and lactational-like cells. In addition, MCFDCIS tumorsphere formation in 3D was reduced by PPARγ activation. In vivo, efatutazone-treated MCFDCIS tumors exhibited fat deposition along with upregulation of PPARγ responsive genes in both epithelial and stromal compartments, suggesting features of milk-producing mammary epithelial cell differentiation. The efatutazone-treated lesions were less invasive with fewer CD44+/p63+ basal progenitor cells. PPARγ activation downregulated Akt phosphorylation in these tumors, although the ERK pathway remained unchanged. Similar trends in gene expression changes consistent with lactational and luminal cell differentiation were observed in the C3(1)/Tag mouse model after efatutazone treatment. CONCLUSIONS: Our data suggest that activation of the PPARγ pathway differentiates DCIS lesions and may be a useful approach to delay DCIS progression.

Comparison of tamoxifen and letrozole response in mammary preneoplasia of ER and aromatase overexpressing mice defines an immune-associated gene signature linked to tamoxifen resistance.

Response to breast cancer chemoprevention can depend upon host genetic makeup and initiating events leading up to preneoplasia. Increased expression of aromatase and estrogen receptor (ER) is found in conjunction with breast cancer. To investigate response or resistance to endocrine therapy, mice with targeted overexpression of Esr1 or CYP19A1 to mammary epithelial cells were employed, representing two direct pathophysiological interventions in estrogen pathway signaling. Both Esr1 and CYP19A1 overexpressing mice responded to letrozole with reduced hyperplastic alveolar nodule prevalence and decreased mammary epithelial cell proliferation. CYP19A1 overexpressing mice were tamoxifen sensitive but Esr1 overexpressing mice were tamoxifen resistant. Increased ER expression occurred with tamoxifen resistance but no consistent changes in progesterone receptor, pSTAT3, pSTAT5, cyclin D1 or cyclin E levels in association with response or resistance were found. RNA-sequencing (RNA-seq) was employed to seek a transcriptome predictive of tamoxifen resistance using these models and a second tamoxifen-resistant model, BRCA1 deficient/Trp53 haploinsufficient mice. Sixty-eight genes associated with immune system processing were upregulated in tamoxifen-resistant Esr1- and Brca1-deficient mice, whereas genes related to aromatic compound metabolic process were upregulated in tamoxifen-sensitive CYP19A1 mice. Interferon regulatory factor 7 was identified as a key transcription factor regulating these 68 immune processing genes. Two loci encoding novel transcripts with high homology to human immunoglobulin lambda-like polypeptide 1 were uniquely upregulated in the tamoxifen-resistant models. Letrozole proved to be a successful alternative to tamoxifen. Further study of transcriptional changes associated with tamoxifen resistance including immune-related genes could expand our mechanistic understanding and lead to biomarkers predictive of escape or response to endocrine therapies.

Oestrogen receptor-alpha regulates non-canonical Hedgehog-signalling in the mammary gland.

Mesenchymal dysplasia (mes) mice harbour a truncation in the C-terminal region of the Hh-ligand receptor, Patched-1 (mPtch1). While the mes variant of mPtch1 binds to Hh-ligands with an affinity similar to that of wild type mPtch1 and appears to normally regulate canonical Hh-signalling via smoothened, the mes mutation causes, among other non-lethal defects, a block to mammary ductal elongation at puberty. We demonstrated previously Hh-signalling induces the activation of Erk1/2 and c-src independently of its control of smo activity. Furthermore, mammary epithelial cell-directed expression of an activated allele of c-src rescued the block to ductal elongation in mes mice, albeit with delayed kinetics. Given that this rescue was accompanied by an induction in estrogen receptor-alpha (ERα) expression and that complex regulatory interactions between ERα and c-src are required for normal mammary gland development, it was hypothesized that expression of ERα would also overcome the block to mammary ductal elongation at puberty in the mes mouse. We demonstrate here that conditional expression of ERα in luminal mammary epithelial cells on the mes background facilitates ductal morphogenesis with kinetics similar to that of the MMTV-c-src(Act) mice. We demonstrate further that Erk1/2 is activated in primary mammary epithelial cells by Shh-ligand and that this activation is blocked by the inhibitor of c-src, PP2, is partially blocked by the ERα inhibitor, ICI 182780 but is not blocked by the smo-inhibitor, SANT-1. These data reveal an apparent Hh-signalling cascade operating through c-src and ERα that is required for mammary gland morphogenesis at puberty.

The methyltransferase EZH2 is not required for mammary cancer development, although high EZH2 and low H3K27me3 correlate with poor prognosis of ER-positive breast cancers.

Enhancer of zeste homolog 2 (EZH2) catalyzes trimethylation of histone H3 lysine 27 (H3K27me3) and its demethylation is catalyzed by UTX. EZH2 levels are frequently elevated in breast cancer and have been proposed to control gene expression through regulating repressive H3K27me3 marks. However, it is not fully established whether breast cancers with different levels of H3K27me3, EZH2 and UTX exhibit different biological behaviors. Levels of H3K27me3, EZH2 and UTX and their prognostic significance were evaluated in 146 cases of breast cancer. H3K27me3 levels were higher in HER2-negative samples. EZH2 expression was higher in cancers that were LN+, size > 20mm, and with higher tumor grade and stage. Using a Cox regression model, H3K27me3 levels and EZH2 expression were identified as independent prognostic factors for overall survival for all the breast cancers studied as well as the ER-positive subgroup. The combination of low H3K27me3 and high EZH2 expression levels were significantly associated with shorter survival. UTX expression was not significantly associated with prognosis and there were no correlations between H3K27me3 levels and EZH2/UTX expression. To determine if EZH2 is required to establish H3K27me3 marks in mammary cancer, Brca1 and Ezh2 were deleted in mammary stem cells in mice. Brca1-deficient mammary cancers with unaltered H3K27me3 levels developed in the absence of EZH2, demonstrating that EZH2 is not a mandatory H3K27 methyltransferase in mammary neoplasia and providing genetic evidence for biological independence between H3K27me3 and EZH2 in this tissue.

The PPARγ agonist efatutazone increases the spectrum of well-differentiated mammary cancer subtypes initiated by loss of full-length BRCA1 in association with TP53 haploinsufficiency.

Peroxisome proliferator-activated receptor gamma (PPARγ) agonists have anticancer activity and influence cell differentiation. We examined the impact of the selective PPARγ agonist efatutazone on mammary cancer pathogenesis in a mouse model of BRCA1 mutation. Mice with conditional loss of full-length BRCA1 targeted to mammary epithelial cells in association with germline TP53 insufficiency were treated with efatutazone through the diet starting at age 4 months and were euthanized at age 12 months or when palpable tumor reached 1 cm(3). Although treatment did not reduce percentage of mice developing invasive cancer, it significantly reduced prevalence of noninvasive cancer and total number of cancers per mouse and increased prevalence of well-differentiated cancer subtypes not usually seen in this mouse model. Invasive cancers from controls were uniformly estrogen receptor α negative and undifferentiated, whereas well-differentiated estrogen receptor α-positive papillary invasive cancers appeared in efatutazone-treated mice. Expression levels of phosphorylated AKT and CDK6 were significantly reduced in the cancers developing in efatutazone-treated mice. Efatutazone treatment reduced rates of mammary epithelial cell proliferation and development of hyperplastic alveolar nodules and increased expression levels of the PPARγ target genes Adfp, Fabp4, and Pdhk4 in preneoplastic mammary tissue. Intervention efatutazone treatment in mice with BRCA1 deficiency altered mammary cancer development by promoting development of differentiated invasive cancer and reducing prevalence of noninvasive cancer and preneoplastic disease.

Variant- and vaccination-specific alternative splicing profiles in SARS-CoV-2 infections.

The COVID-19 pandemic, driven by the SARS-CoV-2 virus and its variants, highlights the important role of understanding host-viral molecular interactions influencing infection outcomes. Alternative splicing post-infection can impact both host responses and viral replication. We analyzed RNA splicing patterns in immune cells across various SARS-CoV-2 variants, considering immunization status. Using a dataset of 190 RNA-seq samples from our prior studies, we observed a substantial deactivation of alternative splicing and RNA splicing-related genes in COVID-19 patients. The alterations varied significantly depending on the infecting variant and immunization history. Notably, Alpha or Beta-infected patients differed from controls, while Omicron-infected patients displayed a splicing profile closer to controls. Particularly, vaccinated Omicron-infected individuals showed a distinct dynamic in alternative splicing patterns not widely shared among other groups. Our findings underscore the intricate interplay between SARS-CoV-2 variants, vaccination-induced immunity, and alternative splicing, emphasizing the need for further investigations to deepen understanding and guide therapeutic development.

Blood transcriptomics analysis offers insights into variant-specific immune response to SARS-CoV-2.

Bulk RNA sequencing (RNA-seq) of blood is typically used for gene expression analysis in biomedical research but is still rarely used in clinical practice. In this study, we propose that RNA-seq should be considered a diagnostic tool, as it offers not only insights into aberrant gene expression and splicing but also delivers additional readouts on immune cell type composition as well as B-cell and T-cell receptor (BCR/TCR) repertoires. We demonstrate that RNA-seq offers insights into a patient's immune status via integrative analysis of RNA-seq data from patients infected with various SARS-CoV-2 variants (in total 196 samples with up to 200 million reads sequencing depth). We compare the results of computational cell-type deconvolution methods (e.g., MCP-counter, xCell, EPIC, quanTIseq) to complete blood count data, the current gold standard in clinical practice. We observe varying levels of lymphocyte depletion and significant differences in neutrophil levels between SARS-CoV-2 variants. Additionally, we identify B and T cell receptor (BCR/TCR) sequences using the tools MiXCR and TRUST4 to show that-combined with sequence alignments and BLASTp-they could be used to classify a patient's disease. Finally, we investigated the sequencing depth required for such analyses and concluded that 10 million reads per sample is sufficient. In conclusion, our study reveals that computational cell-type deconvolution and BCR/TCR methods using bulk RNA-seq analyses can supplement missing CBC data and offer insights into immune responses, disease severity, and pathogen-specific immunity, all achievable with a sequencing depth of 10 million reads per sample.

Variant- and Vaccination-Specific Alternative Splicing Profiles in SARS-CoV-2 Infections.

The COVID-19 pandemic, caused by the coronavirus SARS-CoV-2, and its subsequent variants has underscored the importance of understanding the host-viral molecular interactions to devise effective therapeutic strategies. A significant aspect of these interactions is the role of alternative splicing in modulating host responses and viral replication mechanisms. Our study sought to delineate the patterns of alternative splicing of RNAs from immune cells across different SARS-CoV-2 variants and vaccination statuses, utilizing a robust dataset of 190 RNA-seq samples from our previous studies, encompassing an average of 212 million reads per sample. We identified a dynamic alteration in alternative splicing and genes related to RNA splicing were highly deactivated in COVID-19 patients and showed variant- and vaccination-specific expression profiles. Overall, Omicron-infected patients exhibited a gene expression profile akin to healthy controls, unlike the Alpha or Beta variants. However, significantly, we found identified a subset of infected individuals, most pronounced in vaccinated patients infected with Omicron variant, that exhibited a specific dynamic in their alternative splicing patterns that was not widely shared amongst the other groups. Our findings underscore the complex interplay between SARS-CoV-2 variants, vaccination-induced immune responses, and alternative splicing, emphasizing the necessity for further investigations into these molecular cross-talks to foster deeper understanding and guide strategic therapeutic development.

Blood transcriptomics analysis offers insights into variant-specific immune response to SARS-CoV-2.

Bulk RNA sequencing (RNA-seq) of blood is typically used for gene expression analysis in biomedical research but is still rarely used in clinical practice. In this study, we argue that RNA-seq should be considered a routine diagnostic tool, as it offers not only insights into aberrant gene expression and splicing but also delivers additional readouts on immune cell type composition as well as B-cell and T-cell receptor (BCR/TCR) repertoires. We demonstrate that RNA-seq offers vital insights into a patient's immune status via integrative analysis of RNA-seq data from patients infected with various SARS-CoV-2 variants (in total 240 samples with up to 200 million reads sequencing depth). We compare the results of computational cell-type deconvolution methods (e.g., MCP-counter, xCell, EPIC, quanTIseq) to complete blood count data, the current gold standard in clinical practice. We observe varying levels of lymphocyte depletion and significant differences in neutrophil levels between SARS-CoV-2 variants. Additionally, we identify B and T cell receptor (BCR/TCR) sequences using the tools MiXCR and TRUST4 to show that - combined with sequence alignments and pBLAST - they could be used to classify a patient's disease. Finally, we investigated the sequencing depth required for such analyses and concluded that 10 million reads per sample is sufficient. In conclusion, our study reveals that computational cell-type deconvolution and BCR/TCR methods using bulk RNA-seq analyses can supplement missing CBC data and offer insights into immune responses, disease severity, and pathogen-specific immunity, all achievable with a sequencing depth of 10 million reads per sample.

Analysis of immune responses in patients with CLL after heterologous COVID-19 vaccination.

Patients with chronic lymphocytic leukemia (CLL) treated with B-cell pathway inhibitors and anti-CD20 antibodies exhibit low humoral response rates following SARS-CoV-2 vaccination. To investigate this observation, a prospective single-institution study was conducted comparing peripheral blood mononuclear cell transcriptional response with antibody and T-cell response rates following heterologous BNT162b2/ChAdOx1 vaccination of 15 patients with CLL/small lymphocytic lymphoma (SLL). Two-dose antibody response rate was 40%, increasing to 53% after booster. Patients on Bruton tyrosine kinase inhibitor (BTKi) and venetoclax ± anti-CD20 antibody within 12 months of vaccination responded inferiorly compared with those under BTKi alone. The 2-dose-T-cell response rate was 80%, which increased to 93% after the booster dose. Key transcriptional findings were that interferon-mediated signaling activation including activation of the JAK-STAT pathway generally occurred within days of vaccination, but was independent from the magnitude of the antibody response. Increasing counts of IGHV genes were associated with B-cell reconstitution and improved humoral response rate in the vaccinated patients. T-cell responses in patients with CLL appeared independent of treatment status, whereas higher humoral response rate was associated with BTKi treatment and B-cell reconstitution. Boosting was particularly effective when intrinsic immune status was improved by CLL treatment. Limitations included studying a relatively small cohort, with different treatments and vaccination schedules.

circRNA-sponging: a pipeline for extensive analysis of circRNA expression and their role in miRNA sponging.

Motivation: Circular RNAs (circRNAs) are long noncoding RNAs (lncRNAs) often associated with diseases and considered potential biomarkers for diagnosis and treatment. Among other functions, circRNAs have been shown to act as microRNA (miRNA) sponges, preventing the role of miRNAs that repress their targets. However, there is no pipeline to systematically assess the sponging potential of circRNAs. Results: We developed circRNA-sponging, a nextflow pipeline that (i) identifies circRNAs via backsplicing junctions detected in RNA-seq data, (ii) quantifies their expression values in relation to their linear counterparts spliced from the same gene, (iii) performs differential expression analysis, (iv) identifies and quantifies miRNA expression from miRNA-sequencing (miRNA-seq) data, (v) predicts miRNA binding sites on circRNAs, (vi) systematically investigates potential circRNA-miRNA sponging events, (vii) creates a network of competing endogenous RNAs and (viii) identifies potential circRNA biomarkers. We showed the functionality of the circRNA-sponging pipeline using RNA sequencing data from brain tissues, where we identified two distinct types of circRNAs characterized by a specific ratio of the number of the binding site to the length of the transcript. The circRNA-sponging pipeline is the first end-to-end pipeline to identify circRNAs and their sponging systematically with raw total RNA-seq and miRNA-seq files, allowing us to better indicate the functional impact of circRNAs as a routine aspect in transcriptomic research. Availability and implementation: https://github.com/biomedbigdata/circRNA-sponging. Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Unlocking Translational Potential: Conditionally Reprogrammed Cells in Advancing Breast Cancer Research.

Preclinical in vitro models play an important role in studying cancer cell biology and facilitating translational research, especially in the identification of drug targets and drug discovery studies. This is particularly relevant in breast cancer, where the global burden of disease is quite high based on prevalence and a relatively high rate of lethality. Predictive tools to select patients who will be responsive to invasive or morbid therapies (radiotherapy, chemotherapy, immunotherapy, and/or surgery) are relatively lacking. To be clinically relevant, a model must accurately replicate the biology and cellular heterogeneity of the primary tumor. Addressing these requirements and overcoming the limitations of most existing cancer cell lines, which are typically derived from a single clone, we have recently developed conditional reprogramming (CR) technology. The CR technology refers to a co-culture system of primary human normal or tumor cells with irradiated murine fibroblasts in the presence of a Rho-associated kinase inhibitor to allow the primary cells to acquire stem cell properties and the ability to proliferate indefinitely in vitro without any exogenous gene or viral transfection. This innovative approach fulfills many of these needs and offers an alternative that surpasses the deficiencies associated with traditional cancer cell lines. These CR cells (CRCs) can be reprogrammed to maintain a highly proliferative state and reproduce the genomic and histological characteristics of the parental tissue. Therefore, CR technology may be a clinically relevant model to test and predict drug sensitivity, conduct gene profile analysis and xenograft research, and undertake personalized medicine. This review discusses studies that have applied CR technology to conduct breast cancer research.

Network medicine-based epistasis detection in complex diseases: ready for quantum computing.

Most heritable diseases are polygenic. To comprehend the underlying genetic architecture, it is crucial to discover the clinically relevant epistatic interactions (EIs) between genomic single nucleotide polymorphisms (SNPs)1-3. Existing statistical computational methods for EI detection are mostly limited to pairs of SNPs due to the combinatorial explosion of higher-order EIs. With NeEDL (network-based epistasis detection via local search), we leverage network medicine to inform the selection of EIs that are an order of magnitude more statistically significant compared to existing tools and consist, on average, of five SNPs. We further show that this computationally demanding task can be substantially accelerated once quantum computing hardware becomes available. We apply NeEDL to eight different diseases and discover genes (affected by EIs of SNPs) that are partly known to affect the disease, additionally, these results are reproducible across independent cohorts. EIs for these eight diseases can be interactively explored in the Epistasis Disease Atlas (https://epistasis-disease-atlas.com). In summary, NeEDL is the first application that demonstrates the potential of seamlessly integrated quantum computing techniques to accelerate biomedical research. Our network medicine approach detects higher-order EIs with unprecedented statistical and biological evidence, yielding unique insights into polygenic diseases and providing a basis for the development of improved risk scores and combination therapies.

Limited cross-variant immune response from SARS-CoV-2 Omicron BA.2 in naïve but not previously infected outpatients.

Omicron is currently the dominant SARS-CoV-2 variant and several sublineages have emerged. Questions remain about the impact of previous SARS-CoV-2 exposure on cross-variant immune responses elicited by BA.2 infection compared to BA.1. Here we show that without previous history of COVID-19, BA.2 infection induces a reduced immune response against all variants of concern (VOC) compared to BA.1 infection. The absence of ACE2 binding in sera of previously naïve BA.1 and BA.2 patients indicates a lack of meaningful neutralization. In contrast, anti-spike antibody levels and neutralizing activity greatly increased in the BA.1 and BA.2 patients with a previous history of COVID-19. Transcriptome analyses of peripheral immune cells showed significant differences in immune response and specific antibody generation between BA.1 and BA.2 patients as well as significant differences in expression of specific immune genes. In summary, prior infection status significantly impacts the innate and adaptive immune response against VOC following BA.2 infection.

Limited cross-variant immune response from SARS-CoV-2 Omicron BA.2 in naïve but not previously infected outpatients.

Omicron is currently the dominant SARS-CoV-2 variant and several sublineages have emerged. Questions remain about the impact of previous SARS-CoV-2 exposure on cross-variant immune responses elicited by the SARS-CoV-2 Omicron sublineage BA.2 compared to BA.1. Here we show that without previous history of COVID-19, BA.2 infection induces a reduced immune response against all variants of concern (VOC) compared to BA.1 infection. The absence of ACE2 binding in sera of previously naïve BA.1 and BA.2 patients indicates a lack of meaningful neutralization. In contrast, anti-spike antibody levels and neutralizing activity greatly increased in the BA.1 and BA.2 patients with a previous history of COVID-19. Transcriptome analyses of peripheral immune cells showed significant differences in immune response and specific antibody generation between BA.1 and BA.2 patients as well as significant differences in the expression of specific immune genes. In summary, prior infection status significantly impacts the innate and adaptive immune response against VOC following BA.2 infection.

Heterologous ChAdOx1-BNT162b2 vaccination in Korean cohort induces robust immune and antibody responses that includes Omicron.

Heterologous ChAdOx1-BNT162b2 vaccination induces a stronger immune response than BNT162b2-BNT162b2. Here, we investigated the molecular transcriptome, germline allelic variants of immunoglobulin loci, and anti-Omicron antibody levels in 46 office and lab workers from the Republic of Korea following ChAdOx1-BNT162b2 vaccination. Anti-spike-specific IgG antibody levels against the ancestral SARS-CoV-2 strain increased from 70 AU/ml to 14,000 AU/ml to 142,000 AU/ml one, three and seven days following the second vaccination. Titers against VOC, including Omicron, were two-fold to three-fold lower, yet higher than those measured following BNT162b2-BNT162b2 vaccination. RNA-seq of peripheral immune cells demonstrated activation of interferon pathways with increased IGHV clonal transcripts encoding neutralizing antibodies. scRNA-seq revealed enriched B cell and CD4+ T cell responses in both ChAdOx1-BNT162b2 and BNT162b2-BNT162b2 recipients, but a stronger clonal expansion of memory B cells with ChAdOx1-BNT162b2. In summary, heterologous ChAdOx1-BNT162b2 provides an innate and adaptive immune response that exceeds homologous BNT162b2 vaccination.