Identification of SARS-CoV-2 variants using viral sequencing for the Centers for Disease Control and Prevention genomic surveillance program.

BACKGROUND: The Centers for Disease Control and Prevention contracted with laboratories to sequence the SARS-CoV-2 genome from positive samples across the United States to enable public health officials to investigate the impact of variants on disease severity as well as the effectiveness of vaccines and treatment. Herein we present the initial results correlating RT-PCR quality control metrics with sample collection and sequencing methods from full SARS-CoV-2 viral genomic sequencing of 24,441 positive patient samples between April and June 2021. METHODS: RT-PCR confirmed (N Gene Ct value < 30) positive patient samples, with nucleic acid extracted from saliva, nasopharyngeal and oropharyngeal swabs were selected for viral whole genome SARS-CoV-2 sequencing. Sequencing was performed using Illumina COVIDSeq™ protocol on either the NextSeq550 or NovaSeq6000 systems. Informatic variant calling, and lineage analysis were performed using DRAGEN COVID Lineage applications on Illumina's Basespace cloud analytical system. All sequence data and variant calls were uploaded to NCBI and GISAID. RESULTS: An association was observed between higher sequencing coverage, quality, and samples with a lower Ct value, with < 27 being optimal, across both sequencing platforms and sample collection methods. Both nasopharyngeal swabs and saliva samples were found to be optimal samples of choice for SARS-CoV-2 surveillance sequencing studies, both in terms of strain identification and sequencing depth of coverage, with NovaSeq 6000 providing higher coverage than the NextSeq 550. The most frequent variants identified were the B.1.617.2 Delta (India) and P.1 Gamma (Brazil) variants in the samples sequenced between April 2021 and June 2021. At the time of submission, the most common variant > 99% of positives sequenced was Omicron. CONCLUSION: These initial analyses highlight the importance of sequencing platform, sample collection methods, and RT-PCR Ct values in guiding surveillance efforts. These surveillance studies evaluating genetic changes of SARS-CoV-2 have been identified as critical by the CDC that can affect many aspects of public health including transmission, disease severity, diagnostics, therapeutics, and vaccines.

Splicing factor ESRP1 controls ER-positive breast cancer by altering metabolic pathways.

The epithelial splicing regulatory proteins 1 and 2 (ESRP1 and ESRP2) control the epithelial-to-mesenchymal transition (EMT) splicing program in cancer. However, their role in breast cancer recurrence is unclear. In this study, we report that high levels of ESRP1, but not ESRP2, are associated with poor prognosis in estrogen receptor positive (ER+) breast tumors. Knockdown of ESRP1 in endocrine-resistant breast cancer models decreases growth significantly and alters the EMT splicing signature, which we confirm using TCGA SpliceSeq data of ER+ BRCA tumors. However, these changes are not accompanied by the development of a mesenchymal phenotype or a change in key EMT-transcription factors. In tamoxifen-resistant cells, knockdown of ESRP1 affects lipid metabolism and oxidoreductase processes, resulting in the decreased expression of fatty acid synthase (FASN), stearoyl-CoA desaturase 1 (SCD1), and phosphoglycerate dehydrogenase (PHGDH) at both the mRNA and protein levels. Furthermore, ESRP1 knockdown increases the basal respiration and spare respiration capacity. This study reports a novel role for ESRP1 that could form the basis for the prevention of tamoxifen resistance in ER+ breast cancer.

Altered expression of telomere-associated genes in leukocytes among BRCA1 and BRCA2 carriers.

Telomere dysfunction resulting from telomere shortening and deregulation of shelterin components has been linked to the pathogenesis of age-related disorders, including cancer. Recent evidence suggests that BRCA1/2 (BRCA1 and BRCA2) tumor suppressor gene products play an important role in telomere maintenance. Although telomere shortening has been reported in BRCA1/2 carriers, the direct effects of BRCA1/2 haploinsufficiency on telomere maintenance and predisposition to cancer development are not completely understood. In this study, we assessed the telomere-associated and telomere-proximal gene expression profiles in peripheral blood leukocytes from patients with a BRCA1 or BRCA2 mutation, compared to samples from sporadic and familial breast cancer individuals. We found that 25 genes, including TINF2 gene (a negative regulator of telomere length), were significantly differentially expressed in BRCA1 carriers. Leukocyte telomere length analysis revealed that BRCA1/2 carriers had relatively shorter telomeres than healthy controls. Further, affected BRCA1/2 carriers were well differentiated from unaffected BRCA1/2 carriers by the expression of telomere-proximal genes. Our results link BRCA1/2 haploinsufficiency to changes in telomere length, telomere-associated as well as telomere-proximal gene expression. Thus, this work supports the effect of BRCA1/2 haploinsufficiency in the biology underlying telomere dysfunction in cancer development. Future studies evaluating these findings will require a large study population.

Transcriptional Networks in Single Perivascular Cells Sorted from Human Adipose Tissue Reveal a Hierarchy of Mesenchymal Stem Cells.

Adipose tissue is a rich source of multipotent mesenchymal stem-like cells, located in the perivascular niche. Based on their surface markers, these have been assigned to two main categories: CD31- /CD45- /CD34+ /CD146- cells (adventitial stromal/stem cells [ASCs]) and CD31- /CD45- /CD34- /CD146+ cells (pericytes [PCs]). These populations display heterogeneity of unknown significance. We hypothesized that aldehyde dehydrogenase (ALDH) activity, a functional marker of primitivity, could help to better define ASC and PC subclasses. To this end, the stromal vascular fraction from a human lipoaspirate was simultaneously stained with fluorescent antibodies to CD31, CD45, CD34, and CD146 antigens and the ALDH substrate Aldefluor, then sorted by fluorescence-activated cell sorting. Individual ASCs (n = 67) and PCs (n = 73) selected from the extremities of the ALDH-staining spectrum were transcriptionally profiled by Fluidigm single-cell quantitative polymerase chain reaction for a predefined set (n = 429) of marker genes. To these single-cell data, we applied differential expression and principal component and clustering analysis, as well as an original gene coexpression network reconstruction algorithm. Despite the stochasticity at the single-cell level, covariation of gene expression analysis yielded multiple network connectivity parameters suggesting that these perivascular progenitor cell subclasses possess the following order of maturity: (a) ALDHbr ASC (most primitive); (b) ALDHdim ASC; (c) ALDHbr PC; (d) ALDHdim PC (least primitive). This order was independently supported by specific combinations of class-specific expressed genes and further confirmed by the analysis of associated signaling pathways. In conclusion, single-cell transcriptional analysis of four populations isolated from fat by surface markers and enzyme activity suggests a developmental hierarchy among perivascular mesenchymal stem cells supported by markers and coexpression networks. Stem Cells 2017;35:1273-1289.

Expression levels of SF3B3 correlate with prognosis and endocrine resistance in estrogen receptor-positive breast cancer.

De novo or acquired resistance to endocrine therapy limits its utility in a significant number of estrogen receptor-positive (ER-positive) breast cancers. It is crucial to identify novel targets for therapeutic intervention and improve the success of endocrine therapies. Splicing factor 3b, subunit 1 (SF3B1) mutations are described in luminal breast cancer albeit in low frequency. In this study, we evaluated the role of SF3B1 and SF3B3, critical parts of the SF3b splicing complex, in ER-positive endocrine resistance. To ascertain the role of SF3B1/SF3B3 in endocrine resistance, their expression levels were evaluated in ER-positive/endocrine-resistant cell lines (MCF-7/LCC2 and MCF-7/LCC9) using a real-time quantitative reverse transcription PCR (qRT-PCR). To further determine their clinical relevance, expression analysis was performed in a cohort of 60 paraffin-embedded ER-positive, node-negative breast carcinomas with low, intermediate, and high Oncotype DX recurrence scores. Expression levels of SF3B1 and SF3B3 and their prognostic value were validated in large cohorts using publicly available gene expression data sets including The Cancer Genome Atlas. SF3B1 and SF3B3 levels were significantly increased in ERα-positive cells with acquired tamoxifen (MCF-7/LCC2; both P<0.0002) and fulvestrant/tamoxifen resistance (MCF-7/LCC9; P=0.008 for SF3B1 and P=0.0006 for SF3B3). Expression levels of both MCF-7/LCC2 and MCF-7/LCC9 were not affected by additional treatments with E2 and/or tamoxifen. Furthermore, qRT-PCR analysis confirmed that SF3B3 expression is significantly upregulated in Oncotype DX high-risk groups when compared with low risk (P=0.019). Similarly, in publicly available breast cancer gene expression data sets, overexpression of SF3B3, but not SF3B1, was significantly correlated with overall survival. Furthermore, the correlation was significant in ER-positive, but not in ER-negative tumors.This is the first study to document the role of SF3B3 in endocrine resistance and prognosis in ER-positive breast cancer. Potential strategies for therapeutic targeting of the splicing mechanism(s) need to be evaluated.

Predicting early brain metastases based on clinicopathological factors and gene expression analysis in advanced HER2-positive breast cancer patients.

The overexpression or amplification of the human epidermal growth factor receptor 2 gene (HER2/neu) is associated with high risk of brain metastasis (BM). The identification of patients at highest immediate risk of BM could optimize screening and facilitate interventional trials. We performed gene expression analysis using complementary deoxyribonucleic acid-mediated annealing, selection, extension and ligation and real-time quantitative reverse transcription PCR (qRT-PCR) in primary tumor samples from two independent cohorts of advanced HER2 positive breast cancer patients. Additionally, we analyzed predictive relevance of clinicopathological factors in this series. Study group included discovery Cohort A (84 patients) and validation Cohort B (75 patients). The only independent variables associated with the development of early BM in both cohorts were the visceral location of first distant relapse [Cohort A: hazard ratio (HR) 7.4, 95 % CI 2.4-22.3; p < 0.001; Cohort B: HR 6.1, 95 % CI 1.5-25.6; p = 0.01] and the lack of trastuzumab administration in the metastatic setting (Cohort A: HR 5.0, 95 % CI 1.4-10.0; p = 0.009; Cohort B: HR 10.0, 95 % CI 2.0-100.0; p = 0.008). A profile including 13 genes was associated with early (≤36 months) symptomatic BM in the discovery cohort. This was refined by qRT-PCR to a 3-gene classifier (RAD51, HDGF, TPR) highly predictive of early BM (HR 5.3, 95 % CI 1.6-16.7; p = 0.005; multivariate analysis). However, predictive value of the classifier was not confirmed in the independent validation Cohort B. The presence of visceral metastases and the lack of trastuzumab administration in the metastatic setting apparently increase the likelihood of early BM in advanced HER2-positive breast cancer.

PROGgeneV2: enhancements on the existing database.

BACKGROUND: We recently published PROGgene, a tool that can be used to study prognostic implications of genes in various cancers. The first version of the tool had several areas for improvement. In this paper we present some major enhancements we have made on the existing tool in the new version, PROGgeneV2. RESULTS: In PROGgeneV2, we have made several modifications to enhance survival analysis capability of the tool. First, we have increased the repository of public studies catalogued in our tool by almost two folds. We have also added additional functionalities to perform survival analysis in a variety of new ways. Survival analysis can now be performed on a) single genes b) multiple genes as a signature, c) ratio of expression of two genes, and d) curated/published gene signatures in new version. Users can now also adjust the survival analysis models for available covariates. Users can study prognostic implications of entire gene signatures in different cancer types, which are searchable by keywords. Also, unique to our tool, in the new version, users will be able to upload and use their own datasets to perform survival analysis on genes of interest. CONCLUSIONS: We believe, like its predecessor, PROGGeneV2 will continue to be useful for the scientific community for formulating research hypotheses and designing mechanistic studies. With added datasets PROGgeneV2 is the most comprehensive survival analysis tool available. PROGgeneV2 is available at http://www.compbio.iupui.edu/proggene.

Comprehensive genomic profiling in diabetic nephropathy reveals the predominance of proinflammatory pathways.

Despite advances in the treatment of diabetic nephropathy (DN), currently available therapies have not prevented the epidemic of progressive chronic kidney disease (CKD). The morbidity of CKD, and the inexorable increase in the prevalence of end-stage renal disease, demands more effective approaches to prevent and treat progressive CKD. We undertook next-generation sequencing in a rat model of diabetic nephropathy to study in depth the pathogenic alterations involved in DN with progressive CKD. We employed the obese, diabetic ZS rat, a model that develops diabetic nephropathy, characterized by progressive CKD, inflammation, and fibrosis, the hallmarks of human disease. We then used RNA-seq to examine the combined effects of renal cells and infiltrating inflammatory cells acting as a pathophysiological unit. The comprehensive systems biology analysis of progressive CKD revealed multiple interactions of altered genes that were integrated into morbid networks. These pathological gene assemblies lead to renal inflammation and promote apoptosis and cell cycle arrest in progressive CKD. Moreover, in what is clearly a major therapeutic challenge, multiple and redundant pathways were found to be linked to renal fibrosis, a major cause of kidney loss. We conclude that systems biology applied to progressive CKD in DN can be used to develop novel therapeutic strategies directed to restore critical anomalies in affected gene networks.

Induced Mist1 expression promotes remodeling of mouse pancreatic acinar cells.

BACKGROUND & AIMS: Early embryogenesis involves cell fate decisions that define the body axes and establish pools of progenitor cells. Development does not stop once lineages are specified; cells continue to undergo specific maturation events, and changes in gene expression patterns lead to their unique physiological functions. Secretory pancreatic acinar cells mature postnatally to synthesize large amounts of protein, polarize, and communicate with other cells. The transcription factor MIST1 is expressed by only secretory cells and regulates maturation events. MIST1-deficient acinar cells in mice do not establish apical-basal polarity, properly position zymogen granules, or communicate with adjacent cells, disrupting pancreatic function. We investigated whether MIST1 directly induces and maintains the mature phenotype of acinar cells. METHODS: We analyzed the effects of Cre-mediated expression of Mist1 in adult Mist1-deficient (Mist1(KO)) mice. Pancreatic tissues were collected and analyzed by light and electron microscopy, immunohistochemistry, real-time polymerase chain reaction analysis, and chromatin immunoprecipitation. Primary acini were isolated from mice and analyzed in amylase secretion assays. RESULTS: Induced expression of Mist1 in adult Mist1(KO) mice restored wild-type gene expression patterns in acinar cells. The acinar cells changed phenotypes, establishing apical-basal polarity, increasing the size of zymogen granules, reorganizing the cytoskeletal network, communicating intercellularly (by synthesizing gap junctions), and undergoing exocytosis. CONCLUSIONS: The exocrine pancreas of adult mice can be remodeled by re-expression of the transcription factor MIST1. MIST1 regulates acinar cell maturation and might be used to repair damaged pancreata in patients with pancreatic disorders.

Regulation of microRNA expression by rifampin in human hepatocytes.

Rifampin causes drug interactions by altering hepatic drug metabolism. Because microRNAs (miRNAs) have been shown to regulate genes involved in drug metabolism, we determined the effect of rifampin on the expression of hepatic miRNAs. Primary human hepatocytes from seven subjects were treated with rifampin, and the expression of miRNA and cytochrome P450 (P450) mRNAs was measured by TaqMan assays and RNA-seq, respectively. Rifampin induced the expression of 10 clinically important and 13 additional P450 genes and repressed the expression of 9 other P450 genes (P < 0.05). Rifampin induced the expression of 33 miRNAs and repressed the expression of 35 miRNAs (P < 0.05). Several of these changes were highly negatively correlated with the rifampin-induced changes in the expression of their predicted target P450 mRNAs, supporting the possibility of miRNA-induced regulation of P450 mRNA expression. In addition, several other miRNA changes were positively correlated with the changes in P450 mRNA expression, suggesting similar regulatory mechanisms. Despite the interindividual variability in the rifampin effects on miRNA expression, principal components analysis clearly separated the rifampin-treated samples from the controls. In conclusion, rifampin treatment alters miRNA expression patterns in human hepatocytes, and some of the changes were correlated with the rifampin-induced changes in expression of the P450 mRNAs they are predicted to target.