The lymphocyte guard-IANs: regulation of lymphocyte survival by IAN/GIMAP family proteins

Identification and validation of GIMAP family genes as immune-related prognostic biomarkers in lung adenocarcinoma

Background

The GIMAP family genes play a key role in immune function. Increasing evidence suggests that GIMAP genes were implicated in the tumorigenesis of lung adenocarcinoma (LUAD). This study aimed to investigate the clinical significance of GIMAP family genes in LUAD.

Methods

In this study, we explored the expression, mutation, prognostic value of GIMAP family genes and the correlation with immune microenvironment in LUAD. We further investigated the relationship between GIMAP family genes expression and immunotherapy response in GEO LUAD and melanoma cohorts.

Results

Among the GIMAP family genes, the expression levels of GIMAP1, GIMAP2, GIMAP4, GIMAP5, GIMAP6, GIMAP7, and GIMAP8 were significantly lower in LUAD tumor tissues than normal tissues. Most GIMAP genes were closely related to age, tumor grade and T stage, but not significantly related to sex, N stage and M stage. In the overall population, patients with high expression of GIMAP family genes had a significant longer overall survival (OS). GO and KEGG enrichment analysis showed that GIMAP family genes were highly enriched in immune-related biological process. The expression of GIMAP family genes was positively correlated with immune cell infiltration and immune checkpoint molecules. Furthermore, high expression of GIMAP family genes were correlated with therapeutic response to immunotherapy in LUAD and melanoma patients.

Conclusion

In this study, we identified that GIMAP family genes were significantly associated with immune cell infiltration and immune checkpoint molecules. They potentially play a critical role in anti-tumor immunity and serve as immunotherapy biomarkers.

Immunomodulatory effects of a probiotic combination treatment to improve the survival of Pacific oyster (Crassostrea gigas) larvae against infection by Vibrio coralliilyticus

Introduction

The culture of Pacific oysters (Crassostrea gigas) is of significant socio-economic importance in the U.S. Pacific Northwest and other temperate regions worldwide, with disease outbreaks acting as significant bottlenecks to the successful production of healthy seed larvae. Therefore, the current study aims to describe the mechanisms of a probiotic combination in improving the survival of C. gigas larvae. Specifically, we investigate changes in C. gigas larval gene expression in response to V. coralliilyticus infection with or without a pre-treatment of a novel probiotic combination.

Methods

Treatment groups consisted of replicates of Pacific oyster larvae exposed to a) a combination of four probiotic bacteria at a total concentration of 3.0 x 105 CFU/mL at 18 hours post-fertilization (hpf), b) pathogenic V. coralliilyticus RE22 at a concentration of 6.0 x 103 CFU/mL at 48 hpf, and c) the probiotic combination at 18 hpf and V. coralliilyticus RE22 at 48 hpf. RNA was extracted from washed larvae after 72 hpf, and transcriptome sequencing was used to identify significant differentially expressed genes (DEGs) within each treatment.

Results

Larvae challenged with V. coralliilyticus showed enhanced expression of genes responsible for inhibiting immune signaling (i.e., TNFAIP3, PSMD10) and inducing apoptosis (i.e., CDIP53). However, when pre-treated with the probiotic combination, these genes were no longer differentially expressed relative to untreated control larvae. Additionally, pre-treatment with the probiotic combination increased expression of immune signaling proteins and immune effectors (i.e., IL-17, MyD88). Apparent immunomodulation in response to probiotic treatment corresponds to an increase in the survival of C. gigas larvae infected with V. coralliilyticus by up to 82%.

Discussion

These results indicate that infection with V. coralliilyticus can suppress the larval immune response while also prompting cell death. Furthermore, the results suggest that the probiotic combination treatment negates the deleterious effects of V. coralliilyticus on larval gene expression while stimulating the expression of genes involved in infection defense mechanisms.

Genes differentially expressed between pathogenic and non-pathogenic Entamoeba histolytica clones influence pathogenicity-associated phenotypes by multiple mechanisms

Recently, two genes involved in amoebic liver abscess formation in a mouse model were identified by their differential expression of non-pathogenic (A1np) and pathogenic (B2p) clones of the Entamoeba histolytica isolate HM:1-IMSS. While overexpression of a gene encoding the metallopeptidase EhMP8-2 reduces the virulence of the pathogenic clone B2p, overexpression of the gene ehi_127670 (ehhp127), encoding a hypothetical protein, increases the virulence of the non-pathogenic clone A1np, while silencing this gene in the pathogenic B2p reduces virulence. To understand the role of both molecules in determining the pathogenicity of E. histolytica, silencing, and overexpression transfectants were characterized in detail. Silencing of ehmp8-2, of the homologous gene ehmp8-1, or both in non-pathogenic A1np trophozoites significantly altered the transcript levels of 347, 216, and 58 genes, respectively. This strong change in the expression profiles caused by the silencing of ehmp8-1 and ehmp8-2 implies that these peptidases regulate the expression of numerous genes. Consequently, numerous phenotypic characteristics, including cytopathic, hemolytic, and cysteine peptidase activity, were altered in response to their silencing. Silencing of ehhp127 in pathogenic B2p trophozoites did not affect the expression of other genes, whereas its overexpression in non-pathogenic A1np trophozoites results in an altered expression of approximately 140 genes. EhHP127 is important for trophozoite motility, as its silencing reduces, while its overexpression enhances movement activity. Interestingly, the specific silencing of ehhp127 also significantly affects cytopathic, cysteine peptidase, and hemolytic activities. All three molecules characterized in this study, namely EhMP8-1, EhMP8-2, and EhHP127, are present in amoeba vesicles. The results show that ehmp8-2 and ehhp127 are not only differentially expressed between pathogenic and non-pathogenic amoebae, but that they also significantly affect amoeba pathogenicity-associated phenotypes by completely different mechanisms. This observation suggests that the regulation of amoeba pathogenicity is achieved by a complex network of molecular mechanisms rather than by single factors.

Investigation and verification of GIMAP6 as a robust biomarker for prognosis and tumor immunity in lung adenocarcinoma

BACKGROUND AND AIM

According to previous reports, GTPase of immunity-associated protein 6 (GIMAP6) is essential for autophagy. However, it is unclear how GIMAP6 affects the development and tumor immunity of lung adenocarcinoma (LUAD).

METHODS

In the present study, the role of GIMAP6 in vivo and in vitro was examined using reverse transcription-quantitative PCR, western blotting, and Cell Counting Kit-8, colony formation and Transwell assays. Datasets from The Cancer Genome Atlas and Genotype-Tissue Expression databases were thoroughly analyzed using R software. A nomogram was created using GIMAP6 and prognostic characteristics. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes and Gene Set Enrichment Analysis were applied to explore the potential mechanism of GIMAP6 in lung cancer. The link between GIMAP6 and the immunological landscape was studied using single-cell RNA sequencing datasets from Tumor Immune Estimation Resource (TIMER) 2.0 and Tumor Immune Single-cell Hub.

RESULTS

Patients with high GIMAP6 expression had improved overall and disease-specific survival compared with those patients with low GIMAP6 expression. According to the receiver operating characteristic and calibration curve, the nomogram based on T stage, N stage and GIMAP6 had predictive value for prognosis. According to functional enrichment analysis, GIMAP6 was primarily involved in T-cell receptor signaling pathway, chemokine signaling pathway, cytokine and cytokine receptor interaction. GIMAP6 was shown to be favorably linked with the infiltration of immune cells and immune-related molecules, including cytotoxic T-lymphocyte associated protein 4, programmed death-ligand 1, and T cell immunoreceptor with Ig and ITIM domains, by single-cell sequencing and TIMER2.0 analysis. The role of GIMAP6 in lung cancer cell proliferation, invasion, migration and immunity was experimentally verified.

CONCLUSION

These findings confirmed that GIMAP6 was an effective prognostic molecule that was involved in the regulation of the immune microenvironment of LUAD, and may become a predictor for the efficacy of immunotherapy.

Attenuation of In Vitro and In Vivo Virulence Is Associated with Repression of Gene Expression of AIG1 Gene in Entamoeba histolytica

Entamoeba histolytica virulence results from complex host-parasite interactions implicating multiple amoebic components (e.g., Gal/GalNAc lectin, cysteine proteinases, and amoebapores) and host factors (microbiota and immune response). UG10 is a strain derived from E. histolytica virulent HM-1:IMSS strain that has lost its virulence in vitro and in vivo as determined by a decrease of hemolytic, cytopathic, and cytotoxic activities, increased susceptibility to human complement, and its inability to form liver abscesses in hamsters. We compared the transcriptome of nonvirulent UG10 and its parental HM-1:IMSS strain. No differences in gene expression of the classical virulence factors were observed. Genes downregulated in the UG10 trophozoites encode for proteins that belong to small GTPases, such as Rab and AIG1. Several protein-coding genes, including iron-sulfur flavoproteins and heat shock protein 70, were also upregulated in UG10. Overexpression of the EhAIG1 gene (EHI_180390) in nonvirulent UG10 trophozoites resulted in augmented virulence in vitro and in vivo. Cocultivation of HM-1:IMSS with E. coli O55 bacteria cells reduced virulence in vitro, and the EhAIG1 gene expression was downregulated. In contrast, virulence was increased in the monoxenic strain UG10, and the EhAIG1 gene expression was upregulated. Therefore, the EhAIG1 gene (EHI_180390) represents a novel virulence determinant in E. histolytica.

Prognostic Value of GIMAP4 and Its Role in Promoting Immune Cell Infiltration into Tumor Microenvironment of Lung Adenocarcinoma

GIMAPs are recognized as an important regulator in the carcinogenesis and development of lung cancer, but the function of GIMAP4 in the tumor microenvironment (TME) of lung cancers is unclear. In this study, we investigated the expression and variation of GIMAP4 in lung adenocarcinoma (LUAD), to explore its association with infiltration of immune cells. The Cancer Genome Atlas (TCGA) data and Gene Expression Omnibus (GEO) data were analyzed. Infiltration of immune cells was identified with TIMER (Tumor Immune Estimation Resource) and TISIDB (an integrated repository portal for tumor-immune system interactions). GIMAP4 expression declined in non-small-cell lung cancer (NSCLC), correlated with a poor overall survival (OS) in LUAD, indicating that GIMAP4 was a promising prognostic biomarker in LUAD. GIMAP4 mutation frequency was 1.76% in TCGA cohort and was relevant to the expression of immune components. TIMER and CIBERSORT analysis further confirmed that high GIMAP4 expression possibly promoted immune cell infiltration into the TME, with low GIMAP4 impairing the efficacy of immunotherapies targeting common immune check point inhibitors (ICI). GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses were performed to provide insights into biological processes involved in LUAD. GIMAP4 was expected to be a prognostic biomarker in LUAD and provides potential adjuvant or neoadjuvant therapeutic strategies for targeting ICIs.

GIMAP6 regulates autophagy, immune competence, and inflammation in mice and humans

Inborn errors of immunity (IEIs) unveil regulatory pathways of human immunity. We describe a new IEI caused by mutations in the GTPase of the immune-associated protein 6 (GIMAP6) gene in patients with infections, lymphoproliferation, autoimmunity, and multiorgan vasculitis. Patients and Gimap6^−/− mice show defects in autophagy, redox regulation, and polyunsaturated fatty acid (PUFA)–containing lipids. We find that GIMAP6 complexes with GABARAPL2 and GIMAP7 to regulate GTPase activity. Also, GIMAP6 is induced by IFN-γ and plays a critical role in antibacterial immunity. Finally, we observed that Gimap6^−/− mice died prematurely from microangiopathic glomerulosclerosis most likely due to GIMAP6 deficiency in kidney endothelial cells.

GIMAP7 as a Potential Predictive Marker for Pan-Cancer Prognosis and Immunotherapy Efficacy

Background

Methods

Results

Conclusion

Aberrant GIMAP2 expression affects oral squamous cell carcinoma progression by promoting cell cycle and inhibiting apoptosis

GTPases of immunity-associated protein 2 (GIMAP2) is a GTPase family member associated with T cell survival. However, its mechanisms of action in oral squamous cell carcinoma (OSCC) remain largely unknown. Therefore, the present study aimed to elucidate the possible role of GIMAP2 in OSCC development by investigating its expression levels and molecular mechanisms in OSCC. Reverse transcription quantitative PCR, immunoblotting and immunohistochemistry indicated that GIMAP2 expression was significantly upregulated (P<0.05) in OSCC-derived cell lines and primary OSCC specimens compared with that in their normal counterparts. GIMAP2-knockdown OSCC cells exhibited decreased cell growth, which was associated with cyclin-dependent kinase (CDK)4, CDK6 and phosphorylated Rb downregulation and p53 and p21 upregulation. In addition to cell cycle arrest, GIMAP2 affected anti-apoptotic functions in GIMAP2-knockdown cells by upregulating Bcl-2 and downregulating Bax and Bak. These findings indicated that GIMAP2 may significantly influence OSCC development and apoptosis inhibition and thus is a potential biomarker of OSCC.

Topical Application of Temperature-Sensitive Gel Containing Caerin 1.1 and 1.9 Peptides on TC-1 Tumour-Bearing Mice Induced High-Level Immune Response in the Tumour Microenvironment

The development of topical cream drugs that increase the immune activation of tumour-infiltrating lymphocytes against tumour and chronic viral infection-associated lesions is of great immunotherapeutic significance. This study demonstrates that the topical application of a temperature-sensitive gel containing caerin 1.1 and 1.9 peptides reduces nearly 50% of the tumour weight of HPV16 E6/E7-transformed TC-1 tumour-bearing mice via improving the tumour microenvironment. Confocal microscopy confirms the time-dependent penetration of caerin 1.9 through the epidermal layer of the ear skin structure of mice. Single-cell transcriptomic analysis shows that the caerin 1.1/1.9 gel expands the populations with high immune activation level and largely stimulates the pro-inflammatory activity of NK and dendritic cells. Closely associated with INFα response, Cebpb seems to play a key role in altering the function of all Arg1^hi macrophages in the caerin group. In addition, the caerin gel treatment recruits almost two-fold more activated CD8⁺ T cells to the TME, relative to the untreated tumour, which shows a synergistic effect derived from the regulation of S1pr1, Ccr7, Ms4a4b and Gimap family expression. The TMT10plex-labelling proteomic quantification further demonstrates the activation of interferon-alpha/beta secretion and response to cytokine stimulus by the caerin gel, while the protein contents of several key regulators were elevated by more than 30%, such as Cd5l, Gzma, Ifit1, Irf9 and Stat1. Computational integration of the proteome with the single-cell transcriptome consistently suggested greater activation of NK and T cells with the topical application of caerin peptide gel.

The Prognostic Value of Bone Morphogenetic Proteins and Their Receptors in Lung Adenocarcinoma

Background

Bone morphogenetic proteins (BMPs) regulate tumor progression via binding to their receptors (BMPRs). However, the expression and clinical significance of BMPs/BMPRs in lung adenocarcinoma remain unclear due to a lack of systematic studies.

Methods

This study screened differentially expressed BMPs/BMPRs (deBMPs/BMPRs) in a training dataset combining TCGA-LUAD and GTEx-LUNG and verified them in four GEO datasets. Their prognostic value was evaluated via univariate and multivariate Cox regression analyses. LASSO was performed to construct an initial risk model. Subsequently, after weighted gene co-expression network analysis (WGCNA), differential expression analysis, and univariate Cox regression analysis, hub genes co-expressed with differentially expressed BMPs/BMPRs were filtered out to improve the risk model and explore potential mechanisms. The improved risk model was re-established via LASSO combining hub genes with differentially expressed BMPs/BMPRs as the core. In the testing cohort including 93 lung adenocarcinoma patients, immunohistochemistry (IHC) was performed to verify BMP5 protein expression and its association with prognosis.

Results

BMP2, BMP5, BMP6, GDF10, and ACVRL1 were verified as downregulated in lung adenocarcinoma. Survival analysis identified BMP5 as an independent protective prognostic factor. We also found that BMP5 was significantly correlated with EGFR expression and mutations, suggesting that BMP5 may play a role in targeted therapy. The initial risk model containing only BMP5 showed a significant correlation (HR: 1.71, 95% CI: 1.28−2.28, p: 3e-04) but low prognostic accuracy (AUC of 1-year survival: 0.6, 3-year survival: 0.6, 5-year survival: 0.63). Seventy-nine hub genes co-expressed with BMP5 were identified, and their functions were enriched in cell migration and tumor metastasis. The re-established risk model showed greater prognostic correlation (HR: 2.58, 95% CI: 1.92–3.46, p: 0) and value (AUC of 1-year survival: 0.72, 3-year survival: 0.69, and 5-year survival: 0.68). IHC results revealed that BMP5 protein was also downregulated in lung adenocarcinoma and higher expression was markedly associated with better prognosis (HR: 0.44, 95% CI: 0.23–0.85, p: 0.0145).

Conclusion

BMP5 is a potential crucial target for lung adenocarcinoma treatment based on significant differential expression and superior prognostic value.

The three as: Alternative splicing, alternative polyadenylation and their impact on apoptosis in immune function

The latest advances in next-generation sequencing studies and transcriptomic profiling over the past decade have highlighted a surprising frequency of genes regulated by RNA processing mechanisms in the immune system. In particular, two control steps in mRNA maturation, namely alternative splicing and alternative polyadenylation, are now recognized to occur in the vast majority of human genes. Both have the potential to alter the identity of the encoded protein, as well as control protein abundance or even protein localization or association with other factors. In this review, we will provide a summary of the general mechanisms by which alternative splicing (AS) and alternative polyadenylation (APA) occur, their regulation within cells of the immune system, and their impact on immunobiology. In particular, we will focus on how control of apoptosis by AS and APA is used to tune cell fate during an immune response.

Taxon-Specific Proteins of the Pathogenic Entamoeba Species E. histolytica and E. nuttalli

The human protozoan parasite Entamoeba histolytica can live in the human intestine for months or years without generating any symptoms in the host. For unknown reasons, amoebae can suddenly destroy the intestinal mucosa and become invasive. This can lead to amoebic colitis or extraintestinal amoebiasis whereby the amoebae spread to other organs via the blood vessels, most commonly the liver where abscesses develop. Entamoeba nuttalli is the closest genetic relative of E. histolytica and is found in wild macaques. Another close relative is E. dispar, which asyptomatically infects the human intestine. Although all three species are closely related, only E. histolytica and E. nuttalli are able to penetrate their host’s intestinal epithelium. Lineage-specific genes and gene families may hold the key to understanding differences in virulence among species. Here we discuss those genes found in E. histolytica that have relatives in only one or neither of its sister species, with particular focus on the peptidase, AIG, Ariel, and BspA families.

Arabidopsis immune-associated nucleotide-binding genes repress heat tolerance at the reproductive stage by inhibiting the unfolded protein response and promoting cell death

Plants are vulnerable to heat stress, especially during reproductive development. The heat shock response (HSR) in the cytosol and nucleus, as well as the unfolded protein response (UPR) in the endoplasmic reticulum (ER), are two mechanisms that enable plants to survive heat stress. Excessive heat or ER stresses lead to cell death when the UPR cannot repair stress damage, but the means by which cell survival or death is determined remains unclear. In this study, we used a genome-wide association study (GWAS) to identify that a cluster of five Immune-associated nucleotide-binding protein (IAN) genes (IAN2 to IAN6) is responsible for variation in heat tolerance at the reproductive stage in Arabidopsis thaliana. These IAN genes have both unique and overlapping functions in the negative regulation of heat tolerance, and their loss of function singly or in combination confers increased heat tolerance, measured by a lower number of barren siliques and a higher seedling survival rate under heat. The loss of rice IAN1 gene function also leads to enhanced heat tolerance, suggesting a conserved function of plant IANs. Transcriptome analysis revealed enhanced expression of HSR and UPR genes, as well as reduced cell death, under heat and ER stress in the mutant of IAN6, a major effect member in Arabidopsis. Furthermore, the IAN proteins were found to promote cell death induced by heat stress, ER stress, and cell death-inducing molecules. Thus, the Arabidopsis IAN genes repress heat tolerance, probably through the HSR and UPR and by enhancing the cell death pathway. The IAN2 to IAN6 proteins are partially localized to the ER, suggesting a direct role in the UPR and UPR-mediated cell death. In addition, a natural IAN6 variant from more heat-tolerant Arabidopsis accessions confers greater heat tolerance and induces less cell death compared with the natural variant from less heat-tolerant accessions. The heat-tolerant IAN6 variant is associated with a higher maximum temperature of the warmest month at its collection sites compared with the heat-sensitive variant. Taken together, these results reveal an important role of Arabidopsis IAN2 to IAN6 genes in the regulation of the HSR, UPR, and cell death, and suggest that their natural variations have adaptive functions in heat tolerance.

A human case of GIMAP6 deficiency: a novel primary immune deficiency

The GTPase of immunity-associated proteins (GIMAPs) are a family of genes believed to contribute to lymphocyte development, signaling, and apoptosis, thus playing an important role in immune system homeostasis. While models of gene derangement have been described in both mice and immortalized cell lines, human examples of these diseases remain exceptionally rare. In this manuscript we describe the first documented human cases of a homozygous deleterious GIMAP6 variant in the GIMAP6 gene and their subsequent clinical and immunological phenotype. In order to interrogate the patients’ immune defect, we performed whole-exome sequencing, western blot, flow cytometry analysis, lymphocyte activation and proliferation studies, cytokine release assays, and apoptosis studies. We found two siblings with a predicted deleterious homozygous variant in the GIMAP6 gene with no expression of GIMAP6 protein on western blot. Patients demonstrated accelerated apoptosis, but largely normal lymphocyte subpopulations, activation and proliferation and cytokine release. There appears to be a spectrum of clinical features associated with deficiency of GIMAP6 protein, with one patient suffering lymphopenia and recurrent sinopulmonary infections, and the other clinically asymptomatic. Biallelic variants in the GIMAP6 gene have now been shown to demonstrate disease in humans. The absence of GIMAP6 protein is associated with a spectrum of clinical manifestations and much remains to be learnt about the pathogenic mechanisms underlying this disease. We suggest that biallelic variants in the gene for GIMAP6 should be considered in children with lymphopenia and recurrent sinopulmonary infections.

Differential expression of Exaiptasia pallida GIMAP genes upon induction of apoptosis and autophagy suggests a potential role in cnidarian symbiosis and disease

Coral reefs, one of the world's most productive and diverse ecosystems, are currently threatened by a variety of stressors that result in increased prevalence of both bleaching and disease. Therefore, understanding the molecular mechanisms involved in these responses is critical to mitigate future damage to the reefs. One group of genes that is potentially involved in cnidarian immunity and symbiosis is GTPases of immunity associated proteins (GIMAP). In vertebrates, this family of proteins is involved in regulating the fate of developing lymphocytes and interacts with proteins involved in apoptosis and autophagy. As apoptosis, autophagy and immunity have previously been shown to be involved in cnidarian symbiosis and disease, the goal of this research was to determine the role of cnidarian GIMAPs in these processes using the anemone Exaiptasia pallida To do so, GIMAP genes were characterized in the E. pallida genome and changes in gene expression were measured using qPCR in response to chemical induction of apoptosis, autophagy and treatment with the immune stimulant lipopolysaccharide (LPS) in both aposymbiotic and symbiotic anemones. The results revealed four GIMAP-like genes in E. pallida, referred to as Ep_GIMAPs Induction of apoptosis and autophagy resulted in a general downregulation of Ep_GIMAPs, but no significant changes were observed in response to LPS treatment. This indicates that Ep_GIMAPs may be involved in the regulation of apoptosis and autophagy, and therefore could play a role in cnidarian-dinoflagellate symbiosis. Overall, these results increase our knowledge on the function of GIMAPs in a basal metazoan.

Systemic analyses of expression patterns and clinical features for GIMAPs family members in lung adenocarcinoma

GTPase of immunity-associated proteins (GIMAPs) are frequently prescribed as important components of immune regulation complexes, which were known to play key roles in lung adenocarcinoma. However, little is known about the function of distinct GIMAPs in lung adenocarcinoma. To address this issue, this study investigated the biological function and pathway of GIMAPs in lung adenocarcinoma using multiple public databases. Absent expression of GIMAPs was found in lung adenocarcinoma at mRNA and protein levels. While a purity-corrected value uncovered that all GIMAPs were positively associated with the immune infiltration of lung adenocarcinoma. Furthermore, the expressions of GIMAPs were considered to be negatively associated with clinical cancer stages, patient’s gender and pathological tumor grades in patients with lung adenocarcinoma. Besides, higher mRNA expression of GIMAPs was significantly associated with longer overall survival of patients with lung adenocarcinoma. Taken together, these results may enable GIMAPs family members as diagnostic and survival biomarker candidates or even potential therapeutic targets for patients with lung adenocarcinoma.

Contrasting Immunomodulatory Effects of Probiotic and Pathogenic Bacteria on Eastern Oyster, Crassostrea Virginica, Larvae

Several Vibrio spp. cause acute and severe mortality events in hatcheries where larvae of bivalve mollusks are reared, potentially leading to subsequent shortage of bivalve seed for the grow-out industry. In particular, strains of Vibrio coralliilyticus have been identified as a major cause of disease in Pacific, Crassostrea gigas, and eastern, C. virginica, oyster hatcheries in the United States of America. Probiotic bacteria are an inexpensive, practical, and natural method of disease control. Previous research shows that pretreatment of larval oysters with probiotic bacteria Bacillus pumilus RI06-95 (RI) and Phaeobacter inhibens S4 (S4) significantly decreases mortality caused by experimental challenge with the bacterial pathogen V. coralliilyticus RE22 (RE22). This study aims to characterize the immune response of 6-10-day-old eastern oyster larvae to experimental challenge with pathogen V. coralliilyticus RE22 and probionts RI and S4. Treatments included (a) pathogen and probiont exposure at a concentration of 5 × 10⁴ CFU per mL (~2500 bacterial cells per larva) for a duration of 6 h, (b) probiont exposure at the same concentration for a duration of 24 h, and (c) probiont RI daily treatment of larvae in the hatchery for 4, 11, and 15 days. Differential gene expression analysis compared pathogen or probiotic-treated transcriptomes to unexposed controls. Probiotic and pathogen treatment led to upregulation of transcripts coding for several immune pattern recognition receptors (PRRs) involved in environmental sensing and detection of microbes in oyster larvae. Larval oyster responses to pathogen RE22 suggested suppression of expression of genes in immune signaling pathways (myd88, tak1, nkap), failure in upregulation of immune effector genes, high metabolic demand, and oxidative stress that potentially contributed to mortality. On the other hand, the transcriptomic response to probiotic bacteria RI and S4 suggested activation of immune signaling pathways and expression of immune effectors (e.g., Cv-spi2, mucins and perforin-2). These key features of the host immune response to probiotic bacteria were shared despite the length of probiotic exposure, probiotic species, and the type of environment in which exposures were conducted. This study suggests that pre-exposure of eastern oyster larvae to probiotics for 6-24 h prior to pathogenic challenge leads to a robust and effective immune response that may contribute to protecting larvae from subsequent challenge with V. coralliilyticus RE22. This research provides new insights into host-microbe interactions in larval oysters that could be applied in the management of vibriosis in bivalve hatcheries.

Genome-wide discovery, and computational and transcriptional characterization of an AIG gene family in the freshwater snail Biomphalaria glabrata, a vector for Schistosoma mansoni

BACKGROUND

The AIG (avrRpt2-induced gene) family of GTPases, characterized by the presence of a distinctive AIG1 domain, is mysterious in having a peculiar phylogenetic distribution, a predilection for undergoing expansion and loss, and an uncertain functional role, especially in invertebrates. AIGs are frequently represented as GIMAPs (GTPase of the immunity associated protein family), characterized by presence of the AIG1 domain along with coiled-coil domains. Here we provide an overview of the remarkably expanded AIG repertoire of the freshwater gastropod Biomphalaria glabrata, compare it with AIGs in other organisms, and detail patterns of expression in B. glabrata susceptible or resistant to infection with Schistosoma mansoni, responsible for the neglected tropical disease of intestinal schistosomiasis.

RESULTS

We define the 7 conserved motifs that comprise the AIG1 domain in B. glabrata and detail its association with at least 7 other domains, indicative of functional versatility of B. glabrata AIGs. AIG genes were usually found in tandem arrays in the B. glabrata genome, suggestive of an origin by segmental gene duplication. We found 91 genes with complete AIG1 domains, including 64 GIMAPs and 27 AIG genes without coiled-coils, more than known for any other organism except Danio (with > 100). We defined expression patterns of AIG genes in 12 different B. glabrata organs and characterized whole-body AIG responses to microbial PAMPs, and of schistosome-resistant or -susceptible strains of B. glabrata to S. mansoni exposure. Biomphalaria glabrata AIG genes clustered with expansions of AIG genes from other heterobranch gastropods yet showed unique lineage-specific subclusters. Other gastropods and bivalves had separate but also diverse expansions of AIG genes, whereas cephalopods seem to lack AIG genes.

CONCLUSIONS

The AIG genes of B. glabrata exhibit expansion in both numbers and potential functions, differ markedly in expression between strains varying in susceptibility to schistosomes, and are responsive to immune challenge. These features provide strong impetus to further explore the functional role of AIG genes in the defense responses of B. glabrata, including to suppress or support the development of medically relevant S. mansoni parasites.

The genome of a subterrestrial nematode reveals adaptations to heat

The nematode Halicephalobus mephisto was originally discovered inhabiting a deep terrestrial aquifer 1.3 km underground. H. mephisto can thrive under conditions of abiotic stress including heat and minimal oxygen, where it feeds on a community of both chemolithotrophic and heterotrophic prokaryotes in an unusual ecosystem isolated from the surface biosphere. Here we report the comprehensive genome and transcriptome of this organism, identifying a signature of adaptation: an expanded repertoire of 70 kilodalton heat-shock proteins (Hsp70) and avrRpt2 induced gene 1 (AIG1) proteins. The expanded Hsp70 genes are transcriptionally induced upon growth under heat stress, and we find that positive selection is detectable in several members of this family. We further show that AIG1 may have been acquired by horizontal gene transfer (HGT) from a rhizobial fungus. Over one-third of the genes of H. mephisto are novel, highlighting the divergence of this nematode from other sequenced organisms. This work sheds light on the genomic basis of heat tolerance in a complete subterrestrial eukaryotic genome.

Exploring the expression and preliminary function of chicken Gimap5 gene

GTPase immune-associated protein 5 (Gimap5) plays a key role in maintaining T cell homeostasis, immunological tolerance and inflammatory processes. However, there are no reports on the chicken Gimap5 gene. In this study, the Gimap5 gene was first cloned from chicken and characterized its tissue expression characteristics in different developmental stages. The transcriptional activities of the Gimap5 gene in immune response were identified. The results showed that full-length cDNA sequence of Gimap5 contained 771 bp and encoded a 256-amino acid protein. The Gimap5 gene was transcribed in various tissues and different development stages. The transcriptional activities of Gimap5 gene in the most tissues increased with the development of chicken, but significantly up to peak in liver and large intestine of 10-month-old chicken. The Gimap5 gene exhibited differential transcriptional activities in immune-related tissues in immune responses, with down-regulated in liver (P < 0.01), spleen (P < 0.05) and bursa of Fabricius (P < 0.05), and up-regulated in thymus (P < 0.01). The results show that Gimap5 may be a multifunctional gene involved in tissue function, development and immune response in chicken. These data can provide the foundation for further study of Gimap5.

Elimination of Reference Mapping Bias Reveals Robust Immune Related Allele-Specific Expression in Crossbred Sheep

Pervasive allelic variation at both gene and single nucleotide level (SNV) between individuals is commonly associated with complex traits in humans and animals. Allele-specific expression (ASE) analysis, using RNA-Seq, can provide a detailed annotation of allelic imbalance and infer the existence of cis-acting transcriptional regulation. However, variant detection in RNA-Seq data is compromised by biased mapping of reads to the reference DNA sequence. In this manuscript, we describe an unbiased standardized computational pipeline for allele-specific expression analysis using RNA-Seq data, which we have adapted and developed using tools available under open license. The analysis pipeline we present is designed to minimize reference bias while providing accurate profiling of allele-specific expression across tissues and cell types. Using this methodology, we were able to profile pervasive allelic imbalance across tissues and cell types, at both the gene and SNV level, in Texel×Scottish Blackface sheep, using the sheep gene expression atlas data set. ASE profiles were pervasive in each sheep and across all tissue types investigated. However, ASE profiles shared across tissues were limited, and instead, they tended to be highly tissue-specific. These tissue-specific ASE profiles may underlie the expression of economically important traits and could be utilized as weighted SNVs, for example, to improve the accuracy of genomic selection in breeding programs for sheep. An additional benefit of the pipeline is that it does not require parental genotypes and can therefore be applied to other RNA-Seq data sets for livestock, including those available on the Functional Annotation of Animal Genomes (FAANG) data portal. This study is the first global characterization of moderate to extreme ASE in tissues and cell types from sheep. We have applied a robust methodology for ASE profiling to provide both a novel analysis of the multi-dimensional sheep gene expression atlas data set and a foundation for identifying the regulatory and expressed elements of the genome that are driving complex traits in livestock.

Heterogeneity of human bone marrow and blood natural killer cells defined by single-cell transcriptome

Natural killer (NK) cells are critical to both innate and adaptive immunity. However, the development and heterogeneity of human NK cells are yet to be fully defined. Using single-cell RNA-sequencing technology, here we identify distinct NK populations in human bone marrow and blood, including one population expressing higher levels of immediate early genes indicative of a homeostatic activation. Functionally matured NK cells with high expression of CX3CR1, HAVCR2 (TIM-3), and ZEB2 represents terminally differentiated status with the unique transcriptional profile. Transcriptomic and pseudotime analyses identify a transitional population between CD56^bright and CD56^dim NK cells. Finally, a donor with GATA2^T354M mutation exhibits reduced percentage of CD56^bright NK cells with altered transcriptome and elevated cell death. These data expand our understanding of the heterogeneity and development of human NK cells.

Natural killer (NK) cells are important innate immune cells with diverse functions. Here the authors use single-cell RNA-sequencing of purified human bone marrow and peripheral blood NK cells to define five populations of NK cells with distinct transcriptomic profile to further our understanding of NK development and heterogeneity.

The IMMUNE-ASSOCIATED NUCLEOTIDE-BINDING 9 Protein Is a Regulator of Basal Immunity in Arabidopsis thaliana

A robust regulation of plant immune responses requires a multitude of positive and negative regulators that act in concert. The immune-associated nucleotide-binding (IAN) gene family members are associated with immunity in different organisms, although no characterization of their function has been carried out to date in plants. In this work, we analyzed the expression patterns of IAN genes and found that IAN9 is repressed upon pathogen infection or treatment with immune elicitors. IAN9 encodes a plasma membrane-localized protein that genetically behaves as a negative regulator of immunity. A novel ian9 mutant generated by CRISPR/Cas9 shows increased resistance to Pseudomonas syringae, while transgenic plants overexpressing IAN9 show a slight increase in susceptibility. In vivo immunoprecipitation of IAN9-green fluorescent protein followed by mass spectrometry analysis revealed that IAN9 associates with a previously uncharacterized C3HC4-type RING-finger domain-containing protein that we named IAN9-associated protein 1 (IAP1), which also acts as a negative regulator of basal immunity. Interestingly, neither ian9 or iap1 mutant plants show any obvious developmental phenotype, suggesting that they display enhanced inducible immunity rather than constitutive immune responses. Because both IAN9 and IAP1 have orthologs in important crop species, they could be suitable targets to generate plants more resistant to diseases caused by bacterial pathogens without yield penalty.

Targeted Quantitative Profiling of GTP-Binding Proteins in Cancer Cells Using Isotope-Coded GTP Probes

GTP-binding proteins play important roles in many essential biological processes, including cell signaling, trafficking, and protein synthesis. To assess quantitatively these proteins at the whole proteome level, we developed a high-throughput targeted proteomic method based on the use of isotope-coded GTP probes and multiple-reaction monitoring (MRM) analysis. Targeted proteins were labeled with desthiobiotin-GTP probes, digested with trypsin, and the ensuing desthiobiotin-conjugated peptides were enriched with streptavidin beads for LC-MS/MS analysis. We also established a Skyline MRM library based on shotgun proteomic data acquired for 12 different human cell lines. The library contained 605 tryptic peptides derived from 217 GTP-binding proteins, representing approximately 60% of the annotated human GTP-binding proteome. By using this library, in conjunction with isotope-coded GTP probes and scheduled LC-MRM analysis, we investigated the differential expression of GTP-binding proteins in a pair of primary/metastatic colon cancer cell lines (SW480 and SW620). We were able to quantify 97 GTP-binding proteins, and we further validated the differential expression of several GTP-binding proteins by Western blot analysis. Together, we developed a facile targeted quantitative proteomic method for the high-throughput analysis of GTP-binding proteins and applied the method for probing the altered expression of these proteins involved in colon cancer metastasis.

Oncogenic transcriptional program driven by TAL1 in T-cell acute lymphoblastic leukemia

TAL1/SCL is a prime example of an oncogenic transcription factor that is abnormally expressed in acute leukemia due to the replacement of regulator elements. This gene has also been recognized as an essential regulator of hematopoiesis. TAL1 expression is strictly regulated in a lineage- and stage-specific manner. Such precise control is crucial for the switching of the transcriptional program. The misexpression of TAL1 in immature thymocytes leads to a widespread series of orchestrated downstream events that affect several different cellular machineries, resulting in a lethal consequence, namely T-cell acute lymphoblastic leukemia (T-ALL). In this article, we will discuss the transcriptional regulatory network and downstream target genes, including protein-coding genes and non-coding RNAs, controlled by TAL1 in normal hematopoiesis and T-cell leukemogenesis.

Small GTPase Immunity-Associated Proteins Mediate Resistance to Toxoplasma gondii Infection in Lewis Rat

Rats vary in their susceptibilities to Toxoplasma gondii infection depending on the rat strain. Compared to the T. gondii-susceptible Brown Norway (BN) rat, the Lewis (LEW) rat is extremely resistant to T. gondii Thus, these two rat strains are ideal models for elucidating host mechanisms that are important for host resistance to T. gondii infection. Therefore, in our efforts to unravel molecular factors directing the protective early innate immune response in the LEW rat, we performed RNA sequencing analysis of the LEW versus BN rat with or without T. gondii infection. We identified three candidate small GTPase immunity-associated proteins (GIMAPs) that were upregulated (false discovery rate, 0.05) in the LEW rat in response to T. gondii infection. Subsequently, we engineered T. gondii-susceptible NR8383 rat macrophage cells for overexpression of LEW rat-derived candidate GIMAP 4, 5, and 6. By immunofluorescence analysis we observed that GIMAP 4, 5, and 6 in T. gondii-infected NR8383 cells each colocalized with GRA5, a parasite parasitophorous vacuole membrane (PVM) marker protein, suggesting their translocation to the PVM. Interestingly, overexpression of each candidate GIMAP in T. gondii-infected NR8383 cells induced translocation of LAMP1, a lysosome marker protein, to the T. gondii surface membrane. Importantly, overexpression of GIMAP 4, 5, or 6 individually inhibited intracellular T. gondii growth, with GIMAP 4 having the highest inhibitory effect. Together, our findings indicate that upregulation of GIMAP 4, 5, and 6 contributes to the robust refractoriness of the LEW rat to T. gondii through induction of lysosomal fusion to the otherwise nonfusogenic PVM.

AIG1 affects in vitro and in vivo virulence in clinical isolates of Entamoeba histolytica

The disease state of amebiasis, caused by Entamoeba histolytica, varies from asymptomatic to severe manifestations that include dysentery and extraintestinal abscesses. The virulence factors of the pathogen, and host defense mechanisms, contribute to the outcomes of infection; however, the underlying genetic factors, which affect clinical outcomes, remain to be fully elucidated. To identify these genetic factors in E. histolytica, we used Illumina next-generation sequencing to conduct a comparative genomic analysis of two clinical isolates obtained from diarrheal and asymptomatic patients (strains KU50 and KU27, respectively). By mapping KU50 and KU27 reads to the genome of a reference HM-1:IMSS strain, we identified two genes (EHI_089440 and EHI_176590) that were absent in strain KU27. In KU27, a single AIG1 (avrRpt2-induced gene 1) family gene (EHI_176590) was found to be deleted, from a tandem array of three AIG1 genes, by homologous recombination between the two flanking genes. Overexpression of the EHI_176590 gene, in strain HM-1:IMSS cl6, resulted in increased formation of cell-surface protrusions and enhanced adhesion to human erythrocytes. The EHI_176590 gene was detected by PCR in 56% of stool samples from symptomatic patients infected with E. histolytica, but only in 15% of stool samples from asymptomatic individuals. This suggests that the presence of the EHI_176590 gene is correlated with the outcomes of infection. Taken together, these data strongly indicate that the AIG1 family protein plays a pivotal role in E. histolytica virulence via regulation of host cell adhesion. Our in-vivo experiments, using a hamster liver abscess model, showed that overexpression or gene silencing of EHI_176590 reduced and increased liver abscess formation, respectively. This suggests that the AIG1 genes may have contrasting roles in virulence depending on the genetic background of the parasite and host environment.

Leukemia-Initiating Cells in T-Cell Acute Lymphoblastic Leukemia

T-cell acute lymphoblastic leukemia (T-ALL) is a hematological malignancy characterized by the clonal proliferation of immature T-cell precursors. T-ALL has many similar pathophysiological features to acute myeloid leukemia, which has been extensively studied in the establishment of the cancer stem cell (CSC) theory, but the CSC concept in T-ALL is still debatable. Although leukemia-initiating cells (LICs), which can generate leukemia in a xenograft setting, have been found in both human T-ALL patients and animal models, the nature and origin of LICs are largely unknown. In this review, we discuss recent studies on LICs in T-ALL and the potential mechanisms of LIC emergence in this disease. We focus on the oncogenic transcription factors TAL1, LMO2, and NOTCH1 and highlight the significance of the transcriptional regulatory programs in normal hematopoietic stem cells and T-ALL.

Functional and biochemical characterization of a T cell-associated anti-apoptotic protein, GIMAP6

GTPases of immunity-associated proteins (GIMAPs) are expressed in lymphocytes and regulate survival/death signaling and cell development within the immune system. We found that human GIMAP6 is expressed primarily in T cell lines. By sorting human peripheral blood mononuclear cells and performing quantitative RT-PCR, GIMAP6 was found to be expressed in CD3+ cells. In Jurkat cells that had been knocked down for GIMAP6, treatment with hydrogen peroxide, FasL, or okadaic acid significantly increased cell death/apoptosis. Exogenous expression of GMAP6 protected Huh-7 cells from apoptosis, suggesting that GIMAP6 is an anti-apoptotic protein. Furthermore, knockdown of GIMAP6 not only rendered Jurkat cells sensitive to apoptosis but also accelerated T cell activation under phorbol 12-myristate 13-acetate/ionomycin treatment conditions. Using this experimental system, we also observed a down-regulation of p65 phosphorylation (Ser-536) in GIMAP6 knockdown cells, indicating that GIMAP6 might display anti-apoptotic function through NF-κB activation. The conclusion from the study on cultured T cells was corroborated by the analysis of primary CD3+ T cells, showing that specific knockdown of GIMAP6 led to enhancement of phorbol 12-myristate 13-acetate/ionomycin-mediated activation signals. To characterize the biochemical properties of GIMAP6, we purified the recombinant GIMAP6 to homogeneity and revealed that GIMAP6 had ATPase as well as GTPase activity. We further demonstrated that the hydrolysis activity of GIMAP6 was not essential for its anti-apoptotic function in Huh-7 cells. Combining the expression data, biochemical properties, and cellular features, we conclude that GIMAP6 plays a role in modulating immune function and that it does this by controlling cell death and the activation of T cells.

Aberrant activation of the GIMAP enhancer by oncogenic transcription factors in T-cell acute lymphoblastic leukemia

The transcription factor TAL1/SCL is one of the most prevalent oncogenes in T-cell acute lymphoblastic leukemia (T-ALL), a malignant disorder resulting from leukemic transformation of thymus T-cell precursors. TAL1 is normally expressed in hematopoietic stem cells (HSCs) but is silenced in immature thymocytes. We hypothesize that TAL1 contributes to leukemogenesis by activating genes that are normally repressed in immature thymocytes. Herein, we identified a novel TAL1-regulated super-enhancer controlling the GIMAP locus, which resides within an insulated chromosomal locus in T-ALL cells. The GIMAP genes are expressed in HSCs and mature T-cells but are downregulated during the immature stage of thymocyte differentiation. The GIMAP enhancer is activated by TAL1, RUNX1 and GATA3 in human T-ALL cells but is repressed by E-proteins. Overexpression of human GIMAP genes in immature thymocytes alone does not induce tumorigenesis but accelerates leukemia development in zebrafish. Our results demonstrate that aberrant activation of the GIMAP enhancer contributes to T-cell leukemogenesis.

Dysregulation of GTPase IMAP family members in hepatocellular cancer

Hepatocellular carcinoma (HCC) is one of the most life‑threatening diseases in the world. Members of the GTPase of the immunity‑associated protein (GIMAP) family are important in regulating apoptosis in cancer cells. However, the basic mechanism of GIMAP in HCC remains to be fully elucidated. The present study was performed to investigate the dysregulation of GIMAP family members in HCC. The techniques of polymerase chain reaction analysis, immunohistochemistry and ELISA were used to analyze the expression of GIMAP5 and GIMAP6 in HCC tissues, in matched noncancerous tissue samples, and in blood samples obtained from patients with HCC and healthy subjects. It was found that the mRNA expression levels of GIMAP5 and GIMAP6 were significantly downregulated in the HCC tumor samples, compared with the levels of expression in the matched non‑tumor tissue samples. Similarly, the mRNA expression levels of GIMAP5 and GIMAP6 were also significantly downregulated in the blood samples from patients with HCC, compared with the expression levels in the blood from healthy subjects. At the protein level, it was found that the GIMAP5 and GIMAP6 proteins were expressed at lower levels in the tumor tissue samples, compared with the matched normal tissue samples, and their expression levels were also lower in the blood samples from patients with HCC, compared with the blood samples from the healthy subjects. These data, demonstrating the downregulation of the mRNA and protein expression levels of GIMAP5 and GIMAP6 in the tumor tissues and blood of patients with HCC, suggested the involvement of GIMAP5 and GIMAP6 in the pathogenesis of HCC, and indicate their possible use as diagnostic markers for HCC.

Innate-like functions of natural killer T cell subsets result from highly divergent gene programs

Natural killer T cells (NKT cells) have stimulatory or inhibitory effects on the immune response that can be attributed in part to the existence of functional subsets of NKT cells. These subsets have been characterized only on the basis of the differential expression of a few transcription factors and cell-surface molecules. Here we have analyzed purified populations of thymic NKT cell subsets at both the transcriptomic level and epigenomic level and by single-cell RNA sequencing. Our data indicated that despite their similar antigen specificity, the functional NKT cell subsets were highly divergent populations with many gene-expression and epigenetic differences. Therefore, the thymus 'imprints' distinct gene programs on subsets of innate-like NKT cells that probably impart differences in proliferative capacity, homing, and effector functions.

Pathogen-associated molecular patterns activate expression of genes involved in cell proliferation, immunity and detoxification in the amebocyte-producing organ of the snail Biomphalaria glabrata

The anterior pericardial wall of the snail Biomphalaria glabrata has been identified as a site of hemocyte production, hence has been named the amebocyte-producing organ (APO). A number of studies have shown that exogenous abiotic and biotic substances, including pathogen associated molecular patterns (PAMPs), are able to stimulate APO mitotic activity and/or enlarge its size, implying a role for the APO in innate immunity. The molecular mechanisms underlying such responses have not yet been explored, in part due to the difficulty in obtaining sufficient APO tissue for gene expression studies. By using a modified RNA extraction technique and microarray technology, we investigated transcriptomic responses of APOs dissected from snails at 24 h post-injection with two bacterial PAMPs, lipopolysaccharide (LPS) and peptidoglycan (PGN), or with fucoidan (FCN), which may mimic fucosyl-rich glycan PAMPs on sporocysts of Schistosoma mansoni. Based upon the number of genes differentially expressed, LPS exhibited the strongest activity, relative to saline-injected controls. A concurrent activation of genes involved in cell proliferation, immune response and detoxification metabolism was observed. A gene encoding checkpoint 1 kinase, a key regulator of mitosis, was highly expressed after stimulation by LPS. Also, seven different aminoacyl-tRNA synthetases that play an essential role in protein synthesis were found to be highly expressed. In addition to stimulating genes involved in cell proliferation, the injected substances, especially LPS, also induced expression of a number of immune-related genes including arginase, peptidoglycan recognition protein short form, tumor necrosis factor receptor, ficolin, calmodulin, bacterial permeability increasing proteins and E3 ubiquitin-protein ligase. Importantly, significant up-regulation was observed in four GiMAP (GTPase of immunity-associated protein) genes, a result which provides the first evidence suggesting an immune role of GiMAP in protostome animals. Moreover, altered expression of genes encoding cytochrome P450, glutathione-S-transferase, multiple drug resistance protein as well as a large number of genes encoding enzymes associated with degradation and detoxification metabolism was elicited in response to the injected substances.

Extensive variation between tissues in allele specific expression in an outbred mammal

BACKGROUND

Allele specific gene expression (ASE), with the paternal allele more expressed than the maternal allele or vice versa, appears to be a common phenomenon in humans and mice. In other species the extent of ASE is unknown, and even in humans and mice there are several outstanding questions. These include; to what extent is ASE tissue specific? how often does the direction of allele expression imbalance reverse between tissues? how often is only one of the two alleles expressed? is there a genome wide bias towards expression of the paternal or maternal allele; and finally do genes that are nearby on a chromosome share the same direction of ASE? Here we use gene expression data (RNASeq) from 18 tissues from a single cow to investigate each of these questions in turn, and then validate some of these findings in two tissues from 20 cows.

RESULTS

Between 40 and 100 million sequence reads were generated per tissue across three replicate samples for each of the eighteen tissues from the single cow (the discovery dataset). A bovine gene expression atlas was created (the first from RNASeq data), and differentially expressed genes in each tissue were identified. To analyse ASE, we had access to unambiguously phased genotypes for all heterozygous variants in the cow's whole genome sequence, where these variants were homozygous in the whole genome sequence of her sire, and as a result we were able to map reads to parental genomes, to determine SNP and genes showing ASE in each tissue. In total 25,251 heterozygous SNP within 7985 genes were tested for ASE in at least one tissue. ASE was pervasive, 89 % of genes tested had significant ASE in at least one tissue. This large proportion of genes displaying ASE was confirmed in the two tissues in a validation dataset. For individual tissues the proportion of genes showing significant ASE varied from as low as 8-16 % of those tested in thymus to as high as 71-82 % of those tested in lung. There were a number of cases where the direction of allele expression imbalance reversed between tissues. For example the gene SPTY2D1 showed almost complete paternal allele expression in kidney and thymus, and almost complete maternal allele expression in the brain caudal lobe and brain cerebellum. Mono allelic expression (MAE) was common, with 1349 of 4856 genes (28 %) tested with more than one heterozygous SNP showing MAE. Across all tissues, 54.17 % of all genes with ASE favoured the paternal allele. Genes that are closely linked on the chromosome were more likely to show higher expression of the same allele (paternal or maternal) than expected by chance. We identified several long runs of neighbouring genes that showed either paternal or maternal ASE, one example was five adjacent genes (GIMAP8, GIMAP7 copy1, GIMAP4, GIMAP7 copy 2 and GIMAP5) that showed almost exclusive paternal expression in brain caudal lobe.

CONCLUSIONS

Investigating the extent of ASE across 18 bovine tissues in one cow and two tissues in 20 cows demonstrated 1) ASE is pervasive in cattle, 2) the ASE is often MAE but ranges from MAE to slight overexpression of the major allele, 3) the ASE is most often tissue specific and that more than half the time displays divergent allele specific expression patterns across tissues, 4) across all genes there is a slight bias towards expression of the paternal allele and 5) genes expressing the same parental allele are clustered together more than expected by chance, and there are several runs of large numbers of genes expressing the same parental allele.

GIMAP GTPase family genes: potential modifiers in autoimmune diabetes, asthma, and allergy

GTPase of the immunity-associated protein (GIMAP) family members are differentially regulated during human Th cell differentiation and have been previously connected to immune-mediated disorders in animal studies. GIMAP4 is believed to contribute to the Th cell subtype-driven immunological balance via its role in T cell survival. GIMAP5 has a key role in BB-DR rat and NOD mouse lymphopenia. To elucidate GIMAP4 and GIMAP5 function and role in human immunity, we conducted a study combining genetic association in different immunological diseases and complementing functional analyses. Single nucleotide polymorphisms tagging the GIMAP haplotype variation were genotyped in Finnish type 1 diabetes (T1D) families and in a prospective Swedish asthma and allergic sensitization birth cohort. Initially, GIMAP5 rs6965571 was associated with risk for asthma and allergic sensitization (odds ratio [OR] 3.74, p = 0.00072, and OR 2.70, p = 0.0063, respectively) and protection from T1D (OR 0.64, p = 0.0058); GIMAP4 rs13222905 was associated with asthma (OR 1.28, p = 0.035) and allergic sensitization (OR 1.27, p = 0.0068). However, after false discovery rate correction for multiple testing, only the associations of GIMAP4 with allergic sensitization and GIMAP5 with asthma remained significant. In addition, transcription factor binding sites surrounding the associated loci were predicted. A gene-gene interaction in the T1D data were observed between the IL2RA rs2104286 and GIMAP4 rs9640279 (OR 1.52, p = 0.0064) and indicated between INS rs689 and GIMAP5 rs2286899. The follow-up functional analyses revealed lower IL-2RA expression upon GIMAP4 knockdown and an effect of GIMAP5 rs2286899 genotype on protein expression. Thus, the potential role of GIMAP4 and GIMAP5 as modifiers of immune-mediated diseases cannot be discarded.

Central role of gimap5 in maintaining peripheral tolerance and T cell homeostasis in the gut

Inflammatory bowel disease (IBD) including Crohn's disease and ulcerative colitis is often precipitated by an abnormal immune response to microbiota due to host genetic aberrancies. Recent studies highlight the importance of the host genome and microflora interactions in the pathogenesis of mucosal inflammation including IBD. Specifically, genome-wide (GWAS) and also next-generation sequencing (NGS)—including whole exome or genome sequencing—have uncovered a large number of susceptibility loci that predispose to autoimmune diseases and/or the two phenotypes of IBD. In addition, the generation of “IBD-prone” animal models using both reverse and forward genetic approaches has not only helped confirm the identification of susceptibility loci but also shed critical insight into the underlying molecular and cellular pathways that drive colitis development. In this review, we summarize recent findings derived from studies involving a novel early-onset model of colitis as it develops in GTPase of immunity-associated protein 5- (Gimap5-) deficient mice. In humans, GIMAP5 has been associated with autoimmune diseases although its function is poorly defined. Here, we discuss how defects in Gimap5 function impair immunological tolerance and lymphocyte survival and ultimately drive the development of CD4⁺ T cell-mediated early-onset colitis.

Quantitative changes in Gimap3 and Gimap5 expression modify mitochondrial DNA segregation in mice

Mammalian mitochondrial DNA (mtDNA) is a high-copy maternally inherited genome essential for aerobic energy metabolism. Mutations in mtDNA can lead to heteroplasmy, the co-occurence of two different mtDNA variants in the same cell, which can segregate in a tissue-specific manner affecting the onset and severity of mitochondrial dysfunction. To investigate mechanisms regulating mtDNA segregation we use a heteroplasmic mouse model with two polymorphic neutral mtDNA haplotypes (NZB and BALB) that displays tissue-specific and age-dependent selection for mtDNA haplotypes. In the hematopoietic compartment there is selection for the BALB mtDNA haplotype, a phenotype that can be modified by allelic variants of Gimap3. Gimap3 is a tail-anchored member of the GTPase of the immunity-associated protein (Gimap) family of protein scaffolds important for leukocyte development and survival. Here we show how the expression of two murine Gimap3 alleles from Mus musculus domesticus and M. m. castaneus differentially affect mtDNA segregation. The castaneus allele has incorporated a uORF (upstream open reading frame) in-frame with the Gimap3 mRNA that impairs translation and imparts a negative effect on the steady-state protein abundance. We found that quantitative changes in the expression of Gimap3 and the paralogue Gimap5, which encodes a lysosomal protein, affect mtDNA segregation in the mouse hematopoietic tissues. We also show that Gimap3 localizes to the endoplasmic reticulum and not mitochondria as previously reported. Collectively these data show that the abundance of protein scaffolds on the endoplasmic reticulum and lysosomes are important to the segregation of the mitochondrial genome in the mouse hematopoietic compartment.

Generation and characterisation of mice deficient in the multi-GTPase domain containing protein, GIMAP8

Background

GTPases of the immunity-associated protein family (GIMAPs) are predominantly expressed in mature lymphocytes. Studies of rodents deficient in GIMAP1, GIMAP4, or GIMAP5 have demonstrated that these GTPases regulate lymphocyte survival. In contrast to the other family members, GIMAP8 contains three potential GTP-binding domains (G-domains), a highly unusual feature suggesting a novel function for this protein. To examine a role for GIMAP8 in lymphocyte biology we examined GIMAP8 expression during lymphocyte development. We also generated a mouse deficient in GIMAP8 and examined lymphocyte development and function.

Principal Findings

We show that GIMAP8 is expressed in the very early and late stages of T cell development in the thymus, at late stages during B cell development, and peripheral T and B cells. We find no defects in T or B lymphocyte development in the absence of GIMAP8. A marginal decrease in the number of recirculating bone marrow B cells suggests that GIMAP8 is important for the survival of mature B cells within the bone marrow niche. We also show that deletion of GIMAP8 results in a delay in apoptotic death of mature T cell in vitro in response to dexamethasone or γ-irradiation. However, despite these findings we find that GIMAP8-deficient mice mount normal primary and secondary responses to a T cell dependent antigen.

Conclusions

Despite its unique structure, GIMAP8 is not required for lymphocyte development but appears to have a minor role in maintaining recirculating B cells in the bone marrow niche and a role in regulating apoptosis of mature T cells.

Tubulin- and actin-associating GIMAP4 is required for IFN-γ secretion during Th cell differentiation

Although GTPase of the immunity-associated protein (GIMAP) family are known to be most highly expressed in the cells of the immune system, their function and role remain still poorly characterized. Small GTPases in general are known to be involved in many cellular processes in a cell type-specific manner and to contribute to specific differentiation processes. Among GIMAP family, GIMAP4 is the only member reported to have true GTPase activity, and its transcription is found to be differentially regulated during early human CD4(+) T helper (Th) lymphocyte differentiation. GIMAP4 has been previously connected mainly with T- and B-cell development and survival and T-cell apoptosis. Here we show GIMAP4 to be localized into cytoskeletal elements and with the component of the trans golgi network, which suggests it to have a function in cellular transport processes. We demonstrate that depletion of GIMAP4 with RNAi results in downregulation of endoplasmic reticulum localizing chaperone VMA21. Most importantly, we discovered that GIMAP4 regulates secretion of cytokines in early differentiating human CD4(+) Th lymphocytes and in particular the secretion of interferon-γ also affecting its downstream targets.

Transcriptome of American oysters, Crassostrea virginica, in response to bacterial challenge: insights into potential mechanisms of disease resistance

The American oyster Crassostrea virginica, an ecologically and economically important estuarine organism, can suffer high mortalities in areas in the Northeast United States due to Roseovarius Oyster Disease (ROD), caused by the gram-negative bacterial pathogen Roseovarius crassostreae. The goals of this research were to provide insights into: 1) the responses of American oysters to R. crassostreae, and 2) potential mechanisms of resistance or susceptibility to ROD. The responses of oysters to bacterial challenge were characterized by exposing oysters from ROD-resistant and susceptible families to R. crassostreae, followed by high-throughput sequencing of cDNA samples from various timepoints after disease challenge. Sequence data was assembled into a reference transcriptome and analyzed through differential gene expression and functional enrichment to uncover genes and processes potentially involved in responses to ROD in the American oyster. While susceptible oysters experienced constant levels of mortality when challenged with R. crassostreae, resistant oysters showed levels of mortality similar to non-challenged oysters. Oysters exposed to R. crassostreae showed differential expression of transcripts involved in immune recognition, signaling, protease inhibition, detoxification, and apoptosis. Transcripts involved in metabolism were enriched in susceptible oysters, suggesting that bacterial infection places a large metabolic demand on these oysters. Transcripts differentially expressed in resistant oysters in response to infection included the immune modulators IL-17 and arginase, as well as several genes involved in extracellular matrix remodeling. The identification of potential genes and processes responsible for defense against R. crassostreae in the American oyster provides insights into potential mechanisms of disease resistance.

Expression profiling of genes modulated by minocycline in a rat model of neuropathic pain

Background

The molecular mechanisms underlying neuropathic pain are constantly being studied to create new opportunities to prevent or alleviate neuropathic pain. The aim of our study was to determine the gene expression changes induced by sciatic nerve chronic constriction injury (CCI) that are modulated by minocycline, which can effectively diminish neuropathic pain in animal studies. The genes associated with minocycline efficacy in neuropathic pain should provide insight into the etiology of neuropathic pain and identify novel therapeutic targets.

Results

We screened the ipsilateral dorsal part of the lumbar spinal cord of the rat CCI model for differentially expressed genes. Out of 22,500 studied transcripts, the abundance levels of 93 transcripts were altered following sciatic nerve ligation. Percentage analysis revealed that 54 transcripts were not affected by the repeated administration of minocycline (30 mg/kg, i.p.), but the levels of 39 transcripts were modulated following minocycline treatment. We then selected two gene expression patterns, B1 and B2. The first transcription pattern, B1, consisted of 10 mRNA transcripts that increased in abundance after injury, and minocycline treatment reversed or inhibited the effect of the injury; the B2 transcription pattern consisted of 7 mRNA transcripts whose abundance decreased following sciatic nerve ligation, and minocycline treatment reversed the effect of the injury. Based on the literature, we selected seven genes for further analysis: Cd40, Clec7a, Apobec3b, Slc7a7, and Fam22f from pattern B1 and Rwdd3 and Gimap5 from pattern B2. Additionally, these genes were analyzed using quantitative PCR to determine the transcriptional changes strongly related to the development of neuropathic pain; the ipsilateral DRGs (L4-L6) were also collected and analyzed in these rats using qPCR.

Conclusion

In this work, we confirmed gene expression alterations previously identified by microarray analysis in the spinal cord and analyzed the expression of selected genes in the DRG. Moreover, we reviewed the literature to illustrate the relevance of these findings for neuropathic pain development and therapy. Further studies are needed to elucidate the roles of the individual genes in neuropathic pain and to determine the therapeutic role of minocycline in the rat neuropathic pain model.

Gimap and T cells: a matter of life or death

GTPase immune-associated proteins (Gimap) genes encode evolutionarily conserved GTP-binding proteins that are preferentially expressed in immune cells. Specific members have been shown to be involved in lymphocyte development, or are associated with inflammatory and autoimmune diseases. However, the function of these proteins remains poorly understood, both at the cellular and molecular levels. A new study in this issue of the European Journal of Immunology [Eur. J. Immunol. 2014. 44: 561-572] points to the distinct but partly overlapping functions of two members of this family, Gimap3 and Gimap5, and offers new insight into their potential functions in T cells.

Variability in DNA methylation defines novel epigenetic subgroups of DLBCL associated with different clinical outcomes

Diffuse large B-cell lymphoma (DLBCL) is the most common aggressive form of non-Hodgkin lymphoma with variable biology and clinical behavior. The current classification does not fully explain the biological and clinical heterogeneity of DLBCLs. In this study, we carried out genomewide DNA methylation profiling of 140 DLBCL samples and 10 normal germinal center B cells using the HpaII tiny fragment enrichment by ligation-mediated polymerase chain reaction assay and hybridization to a custom Roche NimbleGen promoter array. We defined methylation disruption as a main epigenetic event in DLBCLs and designed a method for measuring the methylation variability of individual cases. We then used a novel approach for unsupervised hierarchical clustering based on the extent of DNA methylation variability. This approach identified 6 clusters (A-F). The extent of methylation variability was associated with survival outcomes, with significant differences in overall and progression-free survival. The novel clusters are characterized by disruption of specific biological pathways such as cytokine-mediated signaling, ephrin signaling, and pathways associated with apoptosis and cell-cycle regulation. In a subset of patients, we profiled gene expression and genomic variation to investigate their interplay with methylation changes. This study is the first to identify novel epigenetic clusters of DLBCLs and their aberrantly methylated genes, molecular associations, and survival.

Gimap3 and Gimap5 cooperate to maintain T-cell numbers in the mouse

Gimap3 (IAN4) and Gimap5 (IAN5) are highly homologous GTP-binding proteins of the Gimap family. Gimap3 and Gimap5, whose transcripts are abundant in mature lymphocytes, can associate with antiapoptotic Bcl-2 family proteins. While it is established that Gimap5 regulates T-cell survival, the in vivo role of Gimap3 is unclear. Here we report the preparation and characteristics of mouse strains lacking Gimap3 and/or Gimap5. We found that the number of T cells was markedly reduced in mice deficient in both Gimap3 and Gimap5. The defects in T-cell cellularity were more severe in mice lacking both Gimap3 and Gimap5 than in mice lacking only Gimap5. No defects in the cellularity of T cells were detected in mice lacking only Gimap3, whereas bone marrow cells from Gimap3-deficient mice showed reduced T-cell production in a competitive hematopoietic environment. Moreover, retroviral overexpression and short hairpin RNAs-mediated silencing of Gimap3 in bone marrow cells elevated and reduced, respectively, the number of T cells produced in irradiated mice. These results suggest that Gimap3 is a regulator of T-cell numbers in the mouse and that multiple Gimap family proteins cooperate to maintain T-cell survival.

The emerging role of photorespiration and non-photorespiratory peroxisomal metabolism in pathogen defence

Photorespiration represents one of the major highways of primary plant metabolism and is the most prominent example of metabolic cell organelle integration, since the pathway requires the concerted action of plastidial, peroxisomal, mitochondrial and cytosolic enzymes and organellar transport proteins. Oxygenation of ribulose-1,5-bisphosphate by Rubisco leads to the formation of large amounts of 2-phosphoglycolate, which are recycled to 3-phosphoglycerate by the photorespiratory C2 cycle, concomitant with stoichiometric production rates of H2 O2 in peroxisomes. Apart from its significance for agricultural productivity, a secondary function of photorespiration in pathogen defence has emerged only recently. Here, we summarise literature data supporting the crosstalk between photorespiration and pathogen defence and perform a meta-expression analysis of photorespiratory genes during pathogen attack. Moreover, we screened Arabidopsis proteins newly predicted using machine learning methods to be targeted to peroxisomes, the central H2 O2 -producing organelle of photorespiration, for homologues of known pathogen defence proteins and analysed their expression during pathogen infection. The analyses further support the idea that photorespiration and non-photorespiratory peroxisomal metabolism play multi-faceted roles in pathogen defence beyond metabolism of reactive oxygen species.

The acute transcriptional response of the coral Acropora millepora to immune challenge: expression of GiMAP/IAN genes links the innate immune responses of corals with those of mammals and plants

Background

As a step towards understanding coral immunity we present the first whole transcriptome analysis of the acute responses of Acropora millepora to challenge with the bacterial cell wall derivative MDP and the viral mimic poly I:C, defined immunogens provoking distinct but well characterised responses in higher animals.

Results

These experiments reveal similarities with the responses both of arthropods and mammals, as well as coral-specific effects. The most surprising finding was that MDP specifically induced three members of the GiMAP gene family, which has been implicated in immunity in mammals but is absent from Drosophila and Caenorhabditis. Like their mammalian homologs, GiMAP genes are arranged in a tandem cluster in the coral genome.

Conclusions

A phylogenomic survey of this gene family implies ancient origins, multiple independent losses and lineage-specific expansions during animal evolution. Whilst functional convergence cannot be ruled out, GiMAP expression in corals may reflect an ancestral role in immunity, perhaps in phagolysosomal processing.

Structural insights into the mechanism of GTPase activation in the GIMAP family

GTPases of immunity-associated proteins (GIMAPs) are regulators of lymphocyte survival and homeostasis. We previously determined the structural basis of GTP-dependent GIMAP2 scaffold formation on lipid droplets. To understand how its GTP hydrolysis is activated, we screened for other GIMAPs on lipid droplets and identified GIMAP7. In contrast to GIMAP2, GIMAP7 displayed dimerization-stimulated GTP hydrolysis. The crystal structure of GTP-bound GIMAP7 showed a homodimer that assembled via the G domains, with the helical extensions protruding in opposite directions. We identified a catalytic arginine that is supplied to the opposing monomer to stimulate GTP hydrolysis. GIMAP7 also stimulated GTP hydrolysis by GIMAP2 via an analogous mechanism. Finally, we found GIMAP2 and GIMAP7 expression differentially regulated in several human T cell lymphoma lines. Our findings suggest that GTPase activity in the GIMAP family is controlled by homo- and heterodimerization. This may have implications for the differential roles of some GIMAPs in lymphocyte survival.

Tumor-microenvironment interactions studied by zonal transcriptional profiling of squamous cell lung carcinoma

Invasion is a critical step in lung tumor progression. The interaction between tumor cells and their surroundings may play an important role in tumor invasion and metastasis. To better understand the mechanisms of tumor invasion and tumor-microenvironment interactions in lung tumors, total RNA was isolated from the inner tumor, tumor invasion front, adjacent lung, and distant normal lung tissue from 17 patients with primary squamous cell lung carcinoma using punch-aided laser capture microdissection. Messenger RNA expression profiles were obtained by microarray analysis, and microRNA profiles were generated from eight of these samples using TaqMan Low Density Arrays. Statistical analysis of the expression data showed extensive changes in gene expression in the inner tumor and tumor front compared with the normal lung and adjacent lung tissue. Only a few genes were differentially expressed between tumor front and the inner tumor. Several genes were validated by immunohistochemistry. Evaluation of the microRNA data revealed zonal expression differences in nearly a fourth of the microRNAs analyzed. Validation of selected microRNAs by in situ hybridization demonstrated strong expression of hsa-miR-196a in the inner tumor; moderate expression of hsa-miR-224 in the inner tumor and tumor front, and strong expression of hsa-miR-650 in the adjacent lung tissue. Pathway analysis placed the majority of genes differentially expressed between tumor and nontumor cells in intrinsic processes associated with inflammation and extrinsic processes related to lymphocyte physiology. Genes differentially expressed between the inner tumor and the adjacent lung/normal lung tissue affected pathways of arachidonic acid metabolism and eicosanoid signaling.

An insertional mutagenesis screen identifies genes that cooperate with Mll-AF9 in a murine leukemogenesis model

Patients with a t(9;11) translocation (MLL-AF9) develop acute myeloid leukemia (AML), and while in mice the expression of this fusion oncogene also results in the development of myeloid leukemia, it is with long latency. To identify mutations that cooperate with Mll-AF9, we infected neonatal wild-type (WT) or Mll-AF9 mice with a murine leukemia virus (MuLV). MuLV-infected Mll-AF9 mice succumbed to disease significantly faster than controls presenting predominantly with myeloid leukemia while infected WT animals developed predominantly lymphoid leukemia. We identified 88 candidate cancer genes near common sites of proviral insertion. Analysis of transcript levels revealed significantly elevated expression of Mn1, and a trend toward increased expression of Bcl11a and Fosb in Mll-AF9 murine leukemia samples with proviral insertions proximal to these genes. Accordingly, FOSB and BCL11A were also overexpressed in human AML harboring MLL gene translocations. FOSB was revealed to be essential for growth in mouse and human myeloid leukemia cells using shRNA lentiviral vectors in vitro. Importantly, MN1 cooperated with Mll-AF9 in leukemogenesis in an in vivo BM viral transduction and transplantation assay. Together, our data identified genes that define transcription factor networks and important genetic pathways acting during progression of leukemia induced by MLL fusion oncogenes.