Striking differences between knockout and wild-type mice in global gene expression variability

Comparative Proteomic Analysis of Toxoplasma gondii RH Wild-Type and Four SRS29B (SAG1) Knock-Out Clones Reveals Significant Differences between Individual Strains

In T. gondii, as well as in other model organisms, gene knock-out using CRISPR-Cas9 is a suitable tool to identify the role of specific genes. The general consensus implies that only the gene of interest is affected by the knock-out. Is this really the case? In a previous study, we generated knock-out (KO) clones of TgRH88_077450 (SRS29B; SAG1) which differed in the numbers of the integrated dihydrofolate-reductase-thymidylate-synthase (MDHFR-TS) drug-selectable marker. Clones 18 and 33 had a single insertion of MDHFR-TS within SRS29B. Clone 6 was disrupted by the insertion of a short unrelated DNA-sequence, but the marker was integrated elsewhere. In clone 30, the marker was inserted into SRS29B, and several other MDHFR-TS copies were found in the genome. KO and wild-type (WT) tachyzoites had similar shapes, dimensions, and vitality. This prompted us to investigate the impact of genetic engineering on the overall proteome patterns of the four clones as compared to the respective WT. Comparative shotgun proteomics of the five strains was performed. Overall, 3236 proteins were identified. Principal component analysis of the proteomes revealed five distinct clusters corresponding to the five strains by both iTop3 and iLFQ algorithms. Detailed analysis of the differentially expressed proteins revealed that the target of the KO, srs29B, was lacking in all KO clones. In addition to this protein, 20 other proteins were differentially expressed between KO clones and WT or between different KO clones. The protein exhibiting the highest variation between the five strains was srs36D encoded by TgRH_016110. The deregulated expression of SRS36D was further validated by quantitative PCR. Moreover, the transcript levels of three other selected SRS genes, namely SRS36B, SRS46, and SRS57, exhibited significant differences between individual strains. These results indicate that knocking out a given gene may affect the expression of other genes. Therefore, care must be taken when specific phenotypes are regarded as a direct consequence of the KO of a given gene.

Increased gene expression variability in BRCA1-associated and basal-like breast tumours

Purpose

Inherited variants in the cancer susceptibility genes, BRCA1 and BRCA2 account for up to 5% of breast cancers. Multiple gene expression studies have analysed gene expression patterns that maybe associated with BRCA12 pathogenic variant status; however, results from these studies lack consensus. These studies have focused on the differences in population means to identified genes associated with BRCA1/2-carriers with little consideration for gene expression variability, which is also under genetic control and is a feature of cellular function.

Methods

We measured differential gene expression variability in three of the largest familial breast cancer datasets and a 2116 breast cancer meta-cohort. Additionally, we used RNA in situ hybridisation to confirm expression variability of EN1 in an independent cohort of more than 500 breast tumours.

Results

BRCA1-associated breast tumours exhibited a 22.8% (95% CI 22.3–23.2) increase in transcriptome-wide gene expression variability compared to BRCAx tumours. Additionally, 40 genes were associated with BRCA1-related breast cancers that had ChIP-seq data suggestive of enriched EZH2 binding. Of these, two genes (EN1 and IGF2BP3) were significantly variable in both BRCA1-associated and basal-like breast tumours. RNA in situ analysis of EN1 supported a significant (p = 6.3 × 10⁻⁰⁴) increase in expression variability in BRCA1-associated breast tumours.

Conclusion

Our novel results describe a state of increased gene expression variability in BRCA1-related and basal-like breast tumours. Furthermore, genes with increased variability may be driven by changes in DNA occupancy of epigenetic effectors. The variation in gene expression is replicable and led to the identification of novel associations between genes and disease phenotypes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10549-021-06328-y.

Despite Genetic Iron Overload, Hfe-Hemochromatosis Mice Do Not Show Bone Loss

One of the most prevalent genetic iron overload disorders in Caucasians is caused by mutations in the HFE gene. Both HFE patients and Hfe‐mouse models develop a progressive accumulation of iron in the parenchymal cells of various tissues, eventually resulting in liver cirrhosis, hepatocellular carcinoma, cardiomyopathies, hypogonadism, and other pathologies. Clinical data and preclinical models have brought considerable attention to the correlation between iron overload and the development of osteoporosis in HFE/Hfe hemochromatosis. Our study critically challenges this concept. We show that systemic iron overload, at the degree present in Hfe^−/− mice, does not associate with the microarchitecture impairment of long bones, thus excluding a negative effect of iron overload on bone integrity. We further reveal that Hfe actions in osteoblasts and osteoclasts are dispensable for the maintenance of bone and iron homeostasis in mice under steady‐state conditions. We conclude that, despite systemic iron overload, Hfe^−/− mice present normal physiological bone homeostasis. © 2019 The Authors. JBMR Plus in published by Wiley Periodicals, Inc. on behalf of the American Society for Bone and Mineral Research.

Cancer cells copy migratory behavior and exchange signaling networks via extracellular vesicles

Recent data showed that cancer cells from different tumor subtypes with distinct metastatic potential influence each other's metastatic behavior by exchanging biomolecules through extracellular vesicles (EVs). However, it is debated how small amounts of cargo can mediate this effect, especially in tumors where all cells are from one subtype, and only subtle molecular differences drive metastatic heterogeneity. To study this, we have characterized the content of EVs shed in vivo by two clones of melanoma (B16) tumors with distinct metastatic potential. Using the Cre‐LoxP system and intravital microscopy, we show that cells from these distinct clones phenocopy their migratory behavior through EV exchange. By tandem mass spectrometry and RNA sequencing, we show that EVs shed by these clones into the tumor microenvironment contain thousands of different proteins and RNAs, and many of these biomolecules are from interconnected signaling networks involved in cellular processes such as migration. Thus, EVs contain numerous proteins and RNAs and act on recipient cells by invoking a multi‐faceted biological response including cell migration.

Characterization of free fatty acid receptors expression in an obesity rat model with high sucrose diet

INTRODUCTION/AIMS

In recent years, it has been shown that free fatty acids receptors (FFAR) of whose function in the cell surface plays a significant role in the regulation of cell function and nutrition as well are activated by various endogenous ligands, but mainly by fatty acids. Within FFAR of our interest are GPR 41, 43 and 120. The functions of these receptors are varied and dependent on the tissue where they are. The activation and signaling of these receptors, FFAR, are involved in many physiological processes, and currently the target of many drugs in metabolic disorders like obesity, diabetes and atherosclerosis.

MATERIAL AND METHODS

Obesity was induced with hypercaloric diet (HD) in male Wistar rats for 20 weeks (n = 10). At the end, adipose tissue (abdominal and subcutaneous) was taken to perform assays for relative quantification mRNA expression by end-point RT-PCR and protein level expression by Western blot.

RESULTS

These present data have shown for the first time that total mRNA isolation and protein expression from both adipose tissues (abdominal and subcutaneous) of rat in obesity condition yield significative statistical difference among the control versus obese groups, showing that the diet high in carbohydrates modifies the total presence of mRNA and protein level expression of the receptors GPR41, 43 and 120.

CONCLUSIONS

Further comparative methods are in process to clarify whether or not the obesity changes the functional receptors in these two tissues for new pharmacological approaches.

Calcium Oscillation Frequency-Sensitive Gene Regulation and Homeostatic Compensation in Pancreatic β-Cells

Pancreatic islet [Formula: see text]-cells are electrically excitable cells that secrete insulin in an oscillatory fashion when the blood glucose concentration is at a stimulatory level. Insulin oscillations are the result of cytosolic [Formula: see text] oscillations that accompany bursting electrical activity of [Formula: see text]-cells and are physiologically important. ATP-sensitive [Formula: see text] channels (K(ATP) channels) play the key role in setting the overall activity of the cell and in driving bursting, by coupling cell metabolism to the membrane potential. In humans, when there is a defect in K(ATP) channel function, [Formula: see text]-cells fail to respond appropriately to changes in the blood glucose level, and electrical and [Formula: see text] oscillations are lost. However, mice compensate for K(ATP) channel defects in islet [Formula: see text]-cells by employing alternative mechanisms to maintain electrical and [Formula: see text] oscillations. In a recent study, we showed that in mice islets in which K(ATP) channels are genetically knocked out another [Formula: see text] current, provided by inward-rectifying [Formula: see text] channels, is increased. With mathematical modeling, we demonstrated that a sufficient upregulation in these channels can account for the paradoxical electrical bursting and [Formula: see text] oscillations observed in these [Formula: see text]-cells. However, the question of determining the correct level of upregulation that is necessary for this compensation remained unanswered, and this question motivates the current study. [Formula: see text] is a well-known regulator of gene expression, and several examples have been shown of genes that are sensitive to the frequency of the [Formula: see text] signal. In this mathematical modeling study, we demonstrate that a [Formula: see text] oscillation frequency-sensitive gene transcription network can adjust the gene expression level of a compensating [Formula: see text] channel so as to rescue electrical bursting and [Formula: see text] oscillations in a model [Formula: see text]-cell in which the key K(ATP) current is removed. This is done without the prescription of a target [Formula: see text] level, but evolves naturally as a consequence of the feedback between the [Formula: see text]-dependent enzymes and the cell's electrical activity. More generally, the study indicates how [Formula: see text] can provide the link between gene expression and cellular electrical activity that promotes wild-type behavior in a cell following gene knockout.

Osteopontin Affects Insulin Vesicle Localization and Ca2+ Homeostasis in Pancreatic Beta Cells from Female Mice

Type 2 diabetic patients suffer from insulin resistance and reduced insulin secretion. Osteopontin (OPN), a versatile protein expressed in several tissues throughout the body including the islets of Langerhans, has previously been implicated in the development of insulin resistance. Here we have investigated the role of OPN in insulin secretion using an OPN knock out mouse model (OPN^-/-). Ultra-structural analyzes of islets from OPN^-/- and WT mice indicated weaker cell-cell connections between the islet cells in the OPN^-/- mouse compared to WT. Analysis of the insulin granule distribution in the beta cells showed that although OPN^-/- and WT beta cells have the same number of insulin granules OPN^-/- beta cells have significantly fewer docked granules. Both OPN^-/- and WT islets displayed synchronized Ca²⁺ oscillations indicative of an intact beta cell communication. OPN^-/- islets displayed higher intracellular Ca²⁺ concentrations when stimulated with 16.7 mM glucose than WT islets and the initial dip upon elevated glucose concentrations (which is associated with Ca²⁺ uptake into ER) was significantly lower in these islets. Glucose-induced insulin secretion was similar in OPN^-/- and WT islets. Likewise, non-fasted blood glucose levels were the same in both groups. In summary, deletion of OPN results in several minor beta-cell defects that can be compensated for in a healthy system.

Increased Variability of Genomic Transcription in Schizophrenia

Schizophrenia (SZ) is a severe chronic mental disorder with a high heritability. Current microarray analyses typically focus on identifying differentially expressed genes or enriched pathways relevant to phenotypes. Whether there is a variability change of the genomic transcription in diseases has rarely been explored. In this study, we compared coefficient of variation (CV, the ratio of the standard deviation to the mean) of genome transcription of early-onset SZ (EOS) patients with controls in a blood mRNA microarray dataset and a blood microRNA (miRNA) microarray dataset. Furthermore, we compared CV of the expression levels of 17 genes in blood of the 30 patients before and after the 12-week treatment using real-time quantitative PCR (RT-qPCR) analysis. Our results indicated a significant increase of CV of genome transcription in patients compared with controls in both the mRNA and the miRNA datasets. The 30 after-treatment patients showed a significant decrease of CV of gene expression compared with the before-treatment patients. Our study may implicate the blood gene expression variability in SZ, providing further evidence supporting the abnormality of peripheral blood transcriptome in SZ. Given that peripheral blood can be easily collected from patients and followed longitudinally, our results may indicate a new way to facilitate the identification of the signatures of clinical subtypes, their prognosis and treatment response.

Transcriptome-based reconstructions from the murine knockout suggest involvement of the urate transporter, URAT1 (slc22a12), in novel metabolic pathways

URAT1 (slc22a12) was identified as the transporter responsible for renal reabsorption of the medically important compound, uric acid. However, subsequent studies have indicated that other transporters make contributions to this process, and that URAT1 transports other organic anions besides urate (including several in common with the closely related multi-specific renal organic anion transporters, OAT1 (slc22a6) and OAT3 (slc22a8)). These findings raise the possibility that urate transport is not the sole physiological function of URAT1. We previously characterized mice null for the murine ortholog of URAT1 (mURAT1; previously cloned as RST), finding a relatively modest decrement in urate reabsorptive capacity. Nevertheless, there were shifts in the plasma and urinary concentrations of multiple small molecules, suggesting significant metabolic changes in the knockouts. Although these molecules remain unidentified, here we have computationally delineated the biochemical networks consistent with transcriptomic data from the null mice. These analyses suggest alterations in the handling of not only urate but also other putative URAT1 substrates comprising intermediates in nucleotide, carbohydrate, and steroid metabolism. Moreover, the analyses indicate changes in multiple other pathways, including those relating to the metabolism of glycosaminoglycans, methionine, and coenzyme A, possibly reflecting downstream effects of URAT1 loss. Taken together with the available substrate and metabolomic data for the other OATs, our findings suggest that the transport and biochemical functions of URAT1 overlap those of OAT1 and OAT3, and could contribute to our understanding of the relationship between uric acid and the various metabolic disorders to which it has been linked.

•

URAT1 handles multiple substrates suggesting functions beyond urate transport

•

We determined metabolic constraints of gene expression changes in URAT1 null mice

•

These suggest URAT1 involvement in multiple bioenergtic and biosynthetic pathways