Systems genetics analyses reveals key genes related to behavioral traits in the striatum of CFW mice.

The striatum plays a central role in directing many complex behaviors ranging from motor control to action choice and reward learning. In our study, we used 55 CFW mice with rapid decay linkage disequilibrium to systematically mine the striatum-related behavioral functional genes by analyzing their striatal transcriptomes and 79 measured behavioral phenotypic data. By constructing a gene co-expression network, we clustered the genes into 13 modules, with most of them being positively correlated with motor traits. Based on functional annotations as well as Fisher's exact and hypergeometric distribution tests, brown and magenta modules were identified as core modules. They were significantly enriched for striatal-related functional genes. Subsequent Mendelian randomization analysis verified the causal relationship between the core modules and dyskinesia. Through the intra-modular gene connectivity analysis, Adcy5 and Kcnma1 were identified as brown and magenta module hub genes, respectively. Knockouts of both Adcy5 and Kcnma1 lead to motor dysfunction in mice, and KCNMA1 acts as a risk gene for schizophrenia and smoking addiction in humans. We also evaluated the cellular composition of each module and identified oligodendrocytes in the striatum to have a positive role in motor regulation.Significance Statement The striatum plays a central role in guiding many complex behaviours from motor control to action selection and reward learning. Clinically, striatal dysfunction contributes to a variety of neurodegenerative diseases, including the well-known Alzheimer's disease and Huntington's disease. In our study, we systematically mined striatum-associated behavioural function genes using 55 CFW mice. And we validated our findings in multiple ways. We found that Adcy5 and Kcnma1 knockouts in mice lead to motor dysfunction in mice and that Kcnma1 is associated with schizophrenia, a finding that holds true in humans. Finally, we also assessed the role of different cells in the regulation of striatal behaviour and found that oligodendrocytes in the striatum play an active role in movement regulation.

Systems genetics of influenza A virus-infected mice identifies TRIM21 as a critical regulator of pulmonary innate immune response.

Tripartite motif 21 (TRIM21) is a cytosolic Fc receptor that targets antibody-bound, internalized pathogens for destruction. Apart from this intrinsic defense role, TRIM21 is implicated in autoimmune diseases, inflammation, and autophagy. Whether TRIM21 participates in host interactions with influenza A virus (IAV), however, is unknown. By computational modeling of body weight and lung transcriptome data from the BXD parents (C57BL/6 J (B6) and DBA/2 J (D2)) and 41 BXD mouse strains challenged by IAV, we reveal that a Trim21-associated gene network modulates the early host responses to IAV infection. Trim21 transcripts were significantly upregulated in infected mice of both B6 and D2 backgrounds. Its expression was significantly higher in infected D2 than in infected B6 early after infection and significantly correlated with body weight loss. We identified significant trans-eQTL on chromosome 14 that regulates Trim21 expression. Nr1d2 and Il3ra were among the strongest candidate genes. Pathway analysis found Trim21 to be involved in inflammation and immunity related pathways, such as inflammation signaling pathways (TNF, IL-17, and NF-κB), viral detection signaling pathways (NOD-like and RIG-I-like), influenza, and other respiratory viral infections. Knockdown of TRIM21 in human lung epithelial A549 cells significantly augmented IAV-induced expression of IFNB1, IFNL1, CCL5, CXCL10, and IFN-stimulated genes including DDX58 and IFIH1, among others. Our data suggest that a TRIM21-associated gene network is involved in several aspects of inflammation and viral detection mechanisms during IAV infection. We identify and validate TRIM21 as a critical regulator of innate immune responses to IAV in human lung epithelial cells.

Exploring Early-Stage Retinal Neurodegeneration in Murine Pigmentary Glaucoma: Insights From Gene Networks and miRNA Regulation Analyses.

Purpose: Glaucoma is a group of heterogeneous optic neuropathies characterized by the progressive degeneration of retinal ganglion cells. However, the underlying mechanisms have not been understood completely. We aimed to elucidate the genetic network associated with the development of pigmentary glaucoma with DBA/2J (D2) mouse model of glaucoma and corresponding genetic control D2-Gpnmb (D2G) mice carrying the wild type (WT) Gpnmb allele. Methods: Retinas isolated from 13 D2 and 12 D2G mice were subdivided into 2 age groups: pre-onset (1-6 months: samples were collected at approximately 1-2, 2-4, and 5-6 months) and post-onset (7-15 months: samples were collected at approximately 7-9, 10-12, and 13-15 months) glaucoma were compared. Differential gene expression (DEG) analysis and gene-set enrichment analyses were performed. To identify micro-RNAs (miRNAs) that target Gpnmb, miRNA expression levels were correlated with time point matched mRNA expression levels. A weighted gene co-expression network analysis (WGCNA) was performed using the reference BXD mouse population. Quantitative real-time PCR (qRT-PCR) was used to validate Gpnmb and miRNA expression levels. Results: A total of 314 and 86 DEGs were identified in the pre-onset and post-onset glaucoma groups, respectively. DEGs in the pre-onset glaucoma group were associated with the crystallin gene family, whereas those in the post-onset group were related to innate immune system response. Of 1329 miRNAs predicted to target Gpnmb, 3 miRNAs (miR-125a-3p, miR-3076-5p, and miR-214-5p) were selected. A total of 47 genes demonstrated overlapping with the identified DEGs between D2 and D2G, segregated into their time-relevant stages. Gpnmb was significantly downregulated, whereas 2 out of 3 miRNAs were significantly upregulated (P < 0.05) in D2 mice at both 3-and 10-month time points. Conclusions: These findings suggest distinct gene-sets involved in pre-and post-glaucoma in the D2 mouse. We identified three miRNAs regulating Gpnmb in the development of murine pigmentary glaucoma.

Expression Levels of the Tnni3k Gene in the Heart Are Highly Associated with Cardiac and Glucose Metabolism-Related Phenotypes and Functional Pathways.

BACKGROUND: Troponin-I interacting kinase encoded by the TNNI3K gene is expressed in nuclei and Z-discs of cardiomyocytes. Mutations in TNNI3K were identified in patients with cardiac conduction diseases, arrhythmias, and cardiomyopathy. METHODS: We performed cardiac gene expression, whole genome sequencing (WGS), and cardiac function analysis in 40 strains of BXD recombinant inbred mice derived from C57BL/6J (B6) and DBA/2J (D2) strains. Expression quantitative trait loci (eQTLs) mapping and gene enrichment analysis was performed, followed by validation of candidate Tnni3k-regulatory genes. RESULTS: WGS identified compound splicing and missense T659I Tnni3k variants in the D2 parent and some BXD strains (D allele) and these strains had significantly lower Tnni3k expression than those carrying wild-type Tnni3k (B allele). Expression levels of Tnni3k significantly correlated with multiple cardiac (heart rate, wall thickness, PR duration, and T amplitude) and metabolic (glucose levels and insulin resistance) phenotypes in BXDs. A significant cis-eQTL on chromosome 3 was identified for the regulation of Tnni3k expression. Furthermore, Tnni3k-correlated genes were primarily involved in cardiac and glucose metabolism-related functions and pathways. Genes Nodal, Gnas, Nfkb1, Bmpr2, Bmp7, Smad7, Acvr1b, Acvr2b, Chrd, Tgfb3, Irs1, and Ppp1cb were differentially expressed between the B and D alleles. CONCLUSIONS: Compound splicing and T659I Tnni3k variants reduce cardiac Tnni3k expression and Tnni3k levels are associated with cardiac and glucose metabolism-related phenotypes.

Differential Expression of MicroRNAs and Predicted Drug Target in Amyotrophic Lateral Sclerosis.

ALS (Amyotrophic Lateral Sclerosis) is a rare type of neurodegenerative disease. It shows progressive degradation of motor neurons in the brain and spinal cord. At present, there is no treatment available that can completely cure ALS. The available treatments can only increase a patient's life span by a few months. Recently, microRNAs (miRNAs), a sub-class of small non-coding RNAs have been shown to play an essential role in the diagnosis, prognosis, and therapy of ALS. Our study focuses on analyzing differential miRNA profiles and predicting drug targets in ALS using bioinformatics and computational approach. The study identifies eight highly differentially expressed miRNAs in ALS patients, four of which are novel. We identified 42 hub genes for these eight highly expressed miRNAs with Amyloid Precursor Protein (APP) as a candidate gene among them for highly expressed down-regulated miRNA, hsa-miR-455-3p using protein-protein interaction network and Cytoscape analysis. A novel association has been found between hsa-miR-455-3p/APP/serotonergic pathway using KEGG pathway analysis. Also, molecular docking studies have revealed curcumin as a potential drug target that may be used for the treatment of ALS. Thus, the present study has identified four novel miRNA biomarkers: hsa-miR-3613-5p, hsa-miR-24, hsa-miR-3064-5p, and hsa-miR-4455. There is a formation of a novel axis, hsa-miR-455-3p/APP/serotonergic pathway, and curcumin is predicted as a potential drug target for ALS.

Transcriptomic analysis of the Non-Obstructive Azoospermia (NOA) to address gene expression regulation in human testis.

There is a need to understand the molecular basis of testes under Non-Obstructive Azoospermia (NOA), a state of failed spermatogenesis. There has been a lack of attention to the transcriptome at the level of alternatively spliced mRNAs (iso-mRNAs) and the mechanism of gene expression regulation. Hence, we aimed to establish a reliable iso-mRNA profile of NOA-testes, and explore molecular mechanisms - especially those related to gene expression regulation. We sequenced mRNAs from testicular samples of donors with complete spermatogenesis (control samples) and a failure of spermatogenesis (NOA samples). We identified differentially expressed genes and their iso-mRNAs via standard NGS data analyses. We then listed these iso-mRNAs hierarchically based on the extent of consistency of differential quantities across samples and groups, and validated the lists via RT-qPCRs (for 80 iso-mRNAs). In addition, we performed extensive bioinformatic analysis of the splicing features, domains, interactions, and functions of differentially expressed genes and iso-mRNAs. Many top-ranking down-regulated genes and iso-mRNAs, i.e., those down-regulated more consistently across the NOA samples, are associated with mitosis, replication, meiosis, cilium, RNA regulation, and post-translational modifications such as ubiquitination and phosphorylation. Most down-regulated iso-mRNAs correspond to full-length proteins that include all expected domains. The predominance of alternative promoters and termination sites in these iso-mRNAs indicate their gene expression regulation via promoters and UTRs. We compiled a new, comprehensive list of human transcription factors (TFs) and used it to identify TF-'TF gene' interactions with potential significance in down-regulating genes under the NOA condition. The results indicate that RAD51 suppression by HSF4 prevents SP1-activation, and SP1, in turn, could regulate multiple TF genes. This potential regulatory axis and other TF interactions identified in this study could explain the down-regulation of multiple genes in NOA-testes. Such molecular interactions may also have key regulatory roles during normal human spermatogenesis.

Cardiac copper content and its relationship with heart physiology: Insights based on quantitative genetic and functional analyses using BXD family mice.

Background: Copper (Cu) is essential for the functioning of various enzymes involved in important cellular and physiological processes. Although critical for normal cardiac function, excessive accumulation, or deficiency of Cu in the myocardium is detrimental to the heart. Fluctuations in cardiac Cu content have been shown to cause cardiac pathologies and imbalance in systemic Cu metabolism. However, the genetic basis underlying cardiac Cu levels and their effects on heart traits remain to be understood. Representing the largest murine genetic reference population, BXD strains have been widely used to explore genotype-phenotype associations and identify quantitative trait loci (QTL) and candidate genes. Methods: Cardiac Cu concentration and heart function in BXD strains were measured, followed by QTL mapping. The candidate genes modulating Cu homeostasis in mice hearts were identified using a multi-criteria scoring/filtering approach. Results: Significant correlations were identified between cardiac Cu concentration and left ventricular (LV) internal diameter and volumes at end-diastole and end-systole, demonstrating that the BXDs with higher cardiac Cu levels have larger LV chamber. Conversely, cardiac Cu levels negatively correlated with LV posterior wall thickness, suggesting that lower Cu concentration in the heart is associated with LV hypertrophy. Genetic mapping identified six QTLs containing a total of 217 genes, which were further narrowed down to 21 genes that showed a significant association with cardiac Cu content in mice. Among those, Prex1 and Irx3 are the strongest candidates involved in cardiac Cu modulation. Conclusion: Cardiac Cu level is significantly correlated with heart chamber size and hypertrophy phenotypes in BXD mice, while being regulated by multiple genes in several QTLs. Prex1 and Irx3 may be involved in modulating Cu metabolism and its downstream effects and warrant further experimental and functional validations.

Echocardiography phenotyping in murine genetic reference population of BXD strains reveals significant QTLs associated with cardiac function and morphology.

The genetic reference population of recombinant inbred BXD mice has been derived from crosses between C57BL/6J and DBA/2J strains. The DBA/2J parent exhibits cardiomyopathy phenotypes, whereas C57BL/6J has normal heart. BXD mice are sequenced for studying genetic interactions in cardiomyopathies. The study aimed to assess cardiomyopathy traits in BXDs and investigate the quantitative genetic architecture of those traits. Echocardiography, blood pressure, and cardiomyocyte size parameters obtained from 44 strains of BXD family (n > 5/sex) at 4-5 mo of age were associated with heart transcriptomes and expression quantitative trait loci (eQTL) mapping was performed. More than twofold variance in ejection fraction (EF%), fractional shortening (FS%), left ventricular volumes (LVVols), internal dimensions (LVIDs), mass (LVM), and posterior wall (LVPW) thickness was found among BXDs. In male BXDs, eQTL mapping identified Ndrg4 on chromosome 8 QTL to be positively correlated with LVVol and LVID and negatively associated with cardiomyocyte diameter. In female BXDs, significant QTLs were found on chromosomes 7 and 3 to be associated with LVPW and EF% and FS%, respectively, and Josd2, Dap3, and Tpm3 were predicted as strong candidate genes. Our study found variable cardiovascular traits among BXD strains and identified multiple associated QTLs, suggesting an influence of genetic background on expression of echocardiographic and cardiomyocyte diameter traits. Increased LVVol and reduced EF% and FS% represented dilated cardiomyopathy, whereas increased LV mass and wall thickness indicated hypertrophic cardiomyopathy traits. The BXD family is ideal for identifying candidate genes, causal and modifier, that influence cardiovascular phenotypes.NEW & NOTEWORTHY This study aimed to establish a cardiac phenotype-genotype correlation in murine genetic reference population of BXD RI strains by phenotyping the echocardiography, blood pressure, and cardiomyocyte diameter traits and associating each collected phenotype with genetic background. Our study identified several QTLs and candidate genes that have significant association with cardiac hypertrophy, ventricular dilation, and function including systolic hyperfunction and dysfunction.

A systems genetics approach delineates the role of Bcl2 in leukemia pathogenesis.

Leukemia is a group of malignancies of the blood forming tissues, and is characterized by the uncontrolled proliferation of blood cells. In the United States, it accounts for approximately 3.5% and 4% of all cancer-related incidences and mortalities, respectively. The current study aimed to explore the role of Bcl2 and associated genes in leukemia pathogenesis using a systems genetics approach. The transcriptome data from BXD Recombinant Inbred (RI) mice was analyzed to identify the expression of Bcl2 in myeloid cells. eQTL mapping was performed to select the potential chromosomal region and subsequently identify the candidate gene modulating the expression of Bcl2. Furthermore, gene enrichment and protein-protein interaction (PPI) analyses of the Bcl2-coexpressed genes were performed to demonstrate the role of Bcl2 in leukemia pathogenesis. The Bcl2-coexpressed genes were found to be enriched in various hematopoietic system related functions, and multiple pathways related to signaling, immune response, and cancer. The PPI network analysis demonstrated direct interaction of hematopoietic function related genes, such as Bag3, Bak1, Bcl2l11, Bmf, Mapk9, Myc, Ppp2r5c, and Ppp3ca with Bcl2. The eQTL mapping identified a 4.5 Mb genomic region on chromosome 11, potentially regulating the expression of Bcl2. A multi-criteria filtering process identified Top2a, among the genes located in the mapped locus, as the best candidate upstream regulator for Bcl2 expression variation. Hence, the current study provides better insights into the role of Bcl2 in leukemia pathogenesis and demonstrates the significance of our approach in gaining new knowledge on leukemia. Furthermore, our findings from the PPI network analysis and eQTL mapping provide supporting evidence of leukemia-associated genes, which can be further explored for their functional importance in leukemia. DATA AVAILABILITY: The myeloid cell transcriptomic data of the BXD mice used in this study can be accessed through our GeneNetwork (http://www.genenetwork.org) with the accession number of GN144.

Ace2 and Tmprss2 Expressions Are Regulated by Dhx32 and Influence the Gastrointestinal Symptoms Caused by SARS-CoV-2.

Studies showed that the gastrointestinal (GI) tract is one of the most important pathways for SARS-CoV-2 infection and coronavirus disease 2019 (COVID-19). As SARS-CoV-2 cellular entry depends on the ACE2 receptor and TMPRSS2 priming of the spike protein, it is important to understand the molecular mechanisms through which these two proteins and their cognate transcripts interact and influence the pathogenesis of COVID-19. In this study, we quantified the expression, associations, genetic modulators, and molecular pathways for Tmprss2 and Ace2 mRNA expressions in GI tissues using a systems genetics approach and the expanded family of highly diverse BXD mouse strains. The results showed that both Tmprss2 and Ace2 are highly expressed in GI tissues with significant covariation. We identified a significant expression quantitative trait locus on chromosome 7 that controls the expression of both Tmprss2 and Ace2. Dhx32 was found to be the strongest candidate in this interval. Co-expression network analysis demonstrated that both Tmprss2 and Ace2 were located at the same module that is significantly associated with other GI-related traits. Protein-protein interaction analysis indicated that hub genes in this module are linked to circadian rhythms. Collectively, our data suggested that genes with circadian rhythms of expression may have an impact on COVID-19 disease, with implications related to the timing and treatment of COVID-19.

Active suppression of leaflet emergence as a mechanism of simple leaf development.

Angiosperm leaves show extensive shape diversity and are broadly divided into two forms; simple leaves with intact lamina and compound leaves with lamina dissected into leaflets. The mechanistic basis of margin dissection and leaflet initiation has been inferred primarily by analysing compound-leaf architecture, and thus whether the intact lamina of simple leaves has the potential to initiate leaflets upon endogenous gene inactivation remains unclear. Here, we show that the CINCINNATA-like TEOSINTE BRANCHED1, CYCLOIDEA, PROLIFERATING CELL FACTORS (CIN-TCP) transcription factors activate the class II KNOTTED1-LIKE (KNOX-II) genes and the CIN-TCP and KNOX-II proteins together redundantly suppress leaflet initiation in simple leaves. Simultaneous downregulation of CIN-TCP and KNOX-II in Arabidopsis leads to the reactivation of the stemness genes KNOX-I and CUPSHAPED COTYLEDON (CUC) and triggers ectopic organogenesis, eventually converting the simple lamina to a super-compound form that appears to initiate leaflets indefinitely. Thus, a conserved developmental mechanism promotes simple leaf architecture in which CIN-TCP-KNOX-II forms a strong differentiation module that suppresses the KNOX-I-CUC network and leaflet initiation.

The molecular basis of gender disparities in smoking lung cancer patients.

AIMS: Gender disparities exist in smoking-related lung cancer epidemiology, but the molecular basis has not been explored so far. We aimed at identifying genes with gender-bias expression pattern in smoking lung cancer patients for understanding the molecular basis of gender bias in smokers using meta-analysis of microarray gene expression data. MATERIALS AND METHODS: Transcriptome of around 1100 samples from 13 studies were used in the meta-analysis to identify 'Lung Cancer genes specific to Female-Smokers' (LCFS) and 'Lung Cancer genes specific to Male-Smokers' (LCMS). The expression profiles of these genes were validated with an independent microarray report and TCGA-RNA-sequencing data. The molecular interactions, pathway, and other functional annotations were portrayed for the key genes identified. KEY FINDINGS: We identified 1159 gender-biased genes in smoking lung cancer patients. Of these, 400 and 474 genes showed differential expression in cancerous compared to normal lung of women (LCFS) and men (LCMS), respectively. While many up-regulated LCFS were involved in 'immune responses' including T-cell activation, leukocyte cell-cell adhesion, the LCMS were mainly involved in 'positive regulation of gene expression', signaling pathways including RAS, VEGF, insulin-receptor signaling, and 'cell cycle'. SIGNIFICANCE: The strategic-method identified genes, particularly, SNX20, GIMAP6, MTMR2, FAM171B, IDH1, MOBP, FBXO17, LPXN and WIPF1, which were consistently differentially expressed in at least 4 studies, and in agreement with RNA-Seq data. Exploring their functions could be beneficial to the gender-based diagnosis, prognosis, and treatment of lung cancer in smokers. The current meta-analysis supports existing knowledge of sexual-dimorphism of immune responses in cancer.

Therapeutically actionable PAK4 is amplified, overexpressed, and involved in bladder cancer progression.

Muscle-invasive bladder carcinomas (MIBCs) are aggressive genitourinary malignancies. Metastatic urothelial carcinoma of the bladder is generally incurable by current chemotherapy and leads to early mortality. Recent studies have identified molecular subtypes of MIBCs with different sensitivities to frontline therapy, suggesting tumor heterogeneity. We have performed multi-omic profiling of the kinome in bladder cancer patients with the goal of identify therapeutic targets. Our analyses revealed amplification, overexpression, and elevated kinase activity of P21 (RAC1) activated kinase 4 (PAK4) in a subset of Bladder cancer (BLCA). Using bladder cancer cells, we confirmed the role of PAK4 in BLCA cell proliferation and invasion. Furthermore, we observed that a PAK4 inhibitor was effective in curtailing growth of BLCA cells. Transcriptomic analyses identified elevated expression of another kinase, protein tyrosine kinase 6 (PTK6), upon treatment with a PAK4 inhibitor and RNA interference of PAK4. Treatment with a combination of kinase inhibitors (vandetanib and dasatinib) showed enhanced sensitivity compared with either drug alone. Thus, PAK4 may be therapeutically actionable for a subset of MIBC patients with amplified and/or overexpressed PAK4 in their tumors. Our results also indicate that combined inhibition of PAK4 and PTK6 may overcome resistance to PAK4. These observations warrant clinical investigations with selected BLCA patients.

Systematic comparison of the protein-protein interaction databases from a user's perspective.

In absence of periodic systematic comparisons, biologists/bioinformaticians may be forced to make a subjective selection among the many protein-protein interaction (PPI) databases and tools. We conducted a comprehensive compilation and comparison of such resources. We compiled 375 PPI resources, short-listed 125 important ones (both lists are available at startbioinfo.com), and compared the features and coverage of 16 carefully-selected databases related to human PPIs. We quantitatively compared the coverage of 'experimentally verified' as well as 'total' (experimentally verified and predicted) PPIs for these 16 databases. Coverage was compared in two ways: (a) PPIs obtained in response to gene queries using the web interfaces were compared. As a query set, 108 genes expressed differently across tissues (specific to kidney, testis, and uterus, and ubiquitous - i.e., expressed in 43 human normal tissues) or associated with certain diseases (breast cancer, lung cancer, Alzheimer's, cystic fibrosis, diabetes, and cardiomyopathy) were chosen. The coverage was also compared for the well-studied genes versus the less-studied ones. The coverage of the databases for high-quality interactions was separately assessed using a set of literature curated experimentally-proven PPIs (gold standard PPI-set); (b) the back-end-data from 15 PPI databases was downloaded and compared. Combined results from STRING and UniHI covered around 84% of 'experimentally verified' PPIs. Approximately 94% of the 'total' PPIs available across the databases were retrieved by the combined use of hPRINT, STRING, and IID. Among the experimentally verified PPIs found exclusively in each database, STRING contributed around 71% of the hits. The coverage of certain databases was skewed for some gene-types. Analysis with the gold-standard PPI-set revealed that GPS-Prot, STRING, APID, and HIPPIE, each covered ~70% of the curated interactions. The database usage frequencies did not always correlate with their respective advantages, thereby justifying the need for more frequent studies of this nature.

Pro-oncogenic, intra host viral quasispecies in Diffuse large B cell lymphoma patients with occult Hepatitis B Virus infection.

Non Hodgkin lymphoma, predominantly Diffuse Large B-cell Lymphoma (DLBCL) has been reported to have a significant association with Hepatitis B virus (HBV). We investigated the presence of different gene segments of HBV in plasma, B-cells and tumor tissues from DLBCL patients and explored the genetic variability of HBV within and across different compartments in a host using Next Generation Sequencing. Despite all 40 patients being HBV seronegative, 68% showed evidence of occult HBV. Sequencing of these gene segments revealed inter-compartment viral variants in 26% of them, each with at least one non-synonymous mutation. Between compartments, core gene variants revealed Arg94Leu, Glu86Arg and Ser41Thr while X gene variants revealed Phe73Val, Ala44Val, Ser146Ala and Ser147Pro. In tumor compartments per se, several mis-sense mutations were detected, notably the classic T1762A/A1764G mutation in the basal core promoter. In addition, a virus surface antigen mis-sense mutation resulting in M125T was detected in all the samples and could account for surface antigen negativity and occult HBV status. It would be interesting to further explore if a temporal accumulation of viral variants within a favored niche, like patients' lymphocytes, could bestow survival advantage to the virus, and if certain pro-oncogenic HBV variants could drive lymphomagenesis in DLBCL.

Role of estrogen receptor alpha in human cervical cancer-associated fibroblasts: a transcriptomic study.

Cancer-Associated Fibroblasts (CAFs) are crucial in genesis and progression of tumors; however, cervical CAFs (C-CAFs) are not well characterized. Estradiol (E2) has been implicated as a cofactor in human papillomavirus (HPV)-mediated cervical cancer (CxCa), both in animal models and in women using oral contraceptives; however, the exact role of the hormone is unclear. Human C-CAFs have recently been shown to express estrogen receptor alpha (ER-α). We investigated gene expression patterns in ex vivo cultured early and late stage C-CAFs in the context of E2. CAFs were isolated from four patients with early and two patients with late stage CxCa. ER-α expression in CxCa tissues was localized to stromal fibroblast-like cells and confirmed in ex vivo cultured C-CAFs. Two ER antagonists (ICI 182,780 and Methyl Piperidino Pyrazole) were used to unravel ER signaling in CAFs. Microarray technology was used for expression profiling and validated by quantitative reverse transcription PCR. The transcriptomes of C-CAFs across stages indicated their activated state. C-CAFs had gene expression patterns associated with both pro-tumorigenic and pro-inflammatory signaling. Late-stage C-CAFs compared to those of early stage appeared to be more actively metabolizing and cycling but expressed fewer genes related to immune function. We report differential expression profiles between C-CAFs: early vs. late stage and in the presence of ER antagonists. Both ER antagonists seemed to modulate C-CAF function by down regulating genes associated with cell cycle and metabolism, affecting angiogenesis and cancer progression. This study characterized C-CAFs from early and late stage disease, and experiments with ER inhibitors emphasized the probable importance of canonical ER-α signaling. Interfering with paracrine signaling through fibroblast ER-α is worth exploiting as a targeted therapy in CxCa management.