High resolution mapping of Cia3: a common arthritis quantitative trait loci in different species

Identification of the susceptible genes and mechanism underlying the comorbid presence of coronary artery disease and rheumatoid arthritis: a network modularization analysis

OBJECTIVE

The comorbidities of coronary artery disease (CAD) and rheumatoid arthritis (RA) are mutual risk factors, which lead to higher mortality, but the biological mechanisms connecting the two remain unclear. Here, we aimed to identify the risk genes for the comorbid presence of these two complex diseases using a network modularization approach, to offer insights into clinical therapy and drug development for these diseases.

METHOD

The expression profile data of patients CAD with and without RA were obtained from the GEO database (GSE110008). Based on the differentially expressed genes (DEGs), weighted gene co-expression network analysis (WGCNA) was used to construct a gene network, detect co-expression modules, and explore their relation to clinical traits. The Zsummary index, gene significance (GS), and module membership (MM) were utilized to screen the important differentiated modules and hub genes. The GO and KEGG pathway enrichment analysis were applied to analyze potential mechanisms.

RESULT

Based on the 278 DEGs obtained, 41 modules were identified, of which 17 and 24 modules were positively and negatively correlated with the comorbid occurrence of CAD and RA (CAD&RA), respectively. Thirteen modules with Zsummary < 2 were found to be the underlying modules, which may be related to CAD&RA. With GS ≥ 0.5 and MM ≥ 0.8, 49 hub genes were identified, such as ADO, ABCA11P, POT1, ZNF141, GPATCH8, ATF6 and MIA3, etc. The area under the curve values of the representative seven hub genes under the three models (LR, KNN, SVM) were greater than 0.88. Enrichment analysis revealed that the biological functions of the targeted modules were mainly involved in cAMP-dependent protein kinase activity, demethylase activity, regulation of calcium ion import, positive regulation of tyrosine, phosphorylation of STAT protein, and tissue migration, etc. CONCLUSION: Thirteen characteristic modules and 49 susceptibility hub genes were identified, and their corresponding molecular functions may reflect the underlying mechanism of CAD&RA, hence providing insights into the development of clinical therapies against these diseases.

Decreased TMEM40 expression is associated with malignant behavior of cutaneous squamous cell carcinoma and inhibits tumor progression

Cutaneous squamous cell carcinoma (CSCC) is one of the most common types of skin cancer in humans worldwide. The identification and characterization of cancer-associated transmembrane proteins are important for understanding the molecular biology of CSCC. The aim of the present study was to evaluate the expression pattern of transmembrane protein 40 (TMEM40) in CSCC and its clinical significance. The underlying mechanisms were also examined. Reverse transcription-quantitative PCR, western blot and immunohistochemistry analysis were used to determine the relative expression of TMEM40 in CSCC cell lines and clinical tissue samples. The effect of TMEM40 gene silencing on cell proliferation was also evaluated using Cell Counting Kit-8 assays. Wound healing assays, flow cytometry and Transwell assays were used to explore the migration, cell cycle distribution/apoptosis and invasion of CSCC cells following TMEM40 silencing, respectively. In the present study, increased TMEM40 expression was observed in CSCC tissue samples, compared with normal skin, and TMEM40 expression was associated with large tumor size in patients with CSCC. In vitro functional assays indicated that TMEM40 was involved in the regulation of A431 and SCL1 cell growth through its effects on the cell cycle and apoptosis. Silencing TMEM40 in A431 and SCL1 cells resulted in cell cycle arrest at the G₀/G₁ phase and promoted apoptosis. In addition, migration and invasion were significantly inhibited following silencing of TMEM40 expression in CSCC cells. Taken together, the results of the present study indicated that reduced TMEM40 expression could inhibit CSCC development and that TMEM40 may represent a therapeutic target in CSCC.

Base Excision Repair in the Immune System: Small DNA Lesions With Big Consequences

The integrity of the genome is under constant threat of environmental and endogenous agents that cause DNA damage. Endogenous damage is particularly pervasive, occurring at an estimated rate of 10,000–30,000 per cell/per day, and mostly involves chemical DNA base lesions caused by oxidation, depurination, alkylation, and deamination. The base excision repair (BER) pathway is primary responsible for removing and repairing these small base lesions that would otherwise lead to mutations or DNA breaks during replication. Next to preventing DNA mutations and damage, the BER pathway is also involved in mutagenic processes in B cells during immunoglobulin (Ig) class switch recombination (CSR) and somatic hypermutation (SHM), which are instigated by uracil (U) lesions derived from activation-induced cytidine deaminase (AID) activity. BER is required for the processing of AID-induced lesions into DNA double strand breaks (DSB) that are required for CSR, and is of pivotal importance for determining the mutagenic outcome of uracil lesions during SHM. Although uracils are generally efficiently repaired by error-free BER, this process is surprisingly error-prone at the Ig loci in proliferating B cells. Breakdown of this high-fidelity process outside of the Ig loci has been linked to mutations observed in B-cell tumors and DNA breaks and chromosomal translocations in activated B cells. Next to its role in preventing cancer, BER has also been implicated in immune tolerance. Several defects in BER components have been associated with autoimmune diseases, and animal models have shown that BER defects can cause autoimmunity in a B-cell intrinsic and extrinsic fashion. In this review we discuss the contribution of BER to genomic integrity in the context of immune receptor diversification, cancer and autoimmune diseases.

Downregulation of Transmembrane protein 40 by miR-138-5p Suppresses Cell Proliferation and Mobility in Clear Cell Renal Cell Carcinoma

Background

Clear cell renal cell carcinoma (ccRCC) represents approximately 70% of RCC,as the most frequent histological subtype of RCC. MiR-138-5p, a tumor-related microRNA (miRNA), has been reported to be implicated in the diverse types of human malignancies, but its role in ccRCCremains unclear.

Objective

The study was designed to investigate the functional behaviors and regulatory mechanisms of miR-138-5p in ccRCC.

Materials and Methods

Quantitative real-time PCR and western blotting analyses were performed to determine the expression of miR-138-5p and TMEM40 in ccRCC tissues. Pearson's correlation coefficient was utilized to evaluate the correlation between miR-138-5p and TMEM40 expression. The function of miR-138-5p and TMEM40 in the cell proliferation, migration and invasion of ccRCC cells (786-O and ACHN) was assessed by CCK-8, colony formation, wound healing and transwell assay, respectively. A luciferase reporter assay was performed to confirm the direct binding of miR-138-5p to the target gene TMEM40.

Results

We found the expression of miR-138-5p was significantly down-regulated, while TMEM40 was remarkably up-regulated in ccRCC tissues. TMEM40 expression was discovered to be inversely correlated with miR-138-5p expression in ccRCC tissues. Functional studies demonstrated that miR-138-5p overexpression or TMEM40 knockdown significantly suppressed ccRCC cell proliferation, migration and invasion in vitro. Notably, we experimentally confirmed that miR-138-5p directly recognizes the 3'-UTR of the TMEM40 transcript and down-regulated its expression in ccRCC cells.

Conclusions

Taken together, our findings provide the first clues regarding the role of miR-138-5p as a tumor suppressor in ccRCC by directly targeting of TMEM40.

High expression of TMEM40 is associated with the malignant behavior and tumorigenesis in bladder cancer

BACKGROUND

Bladder cancer (BCa) is one of the most common cancers in the urinary system among the world. Previous studies suggested that TMEM40 expression level was significantly associated with clinicopathological parameters including histological grade, clinical stage and pT status of bladder cancer. However, the molecular mechanism of TMEM40 in BCa remains poorly understood.

METHODS

Real-time quantitative RT-PCR (qRT-PCR) and western blot (WB) were used to examine the expression levels of TMEM40 in BCa tissues, paired non-cancer tissues and cell lines. A series of experiments, including CCK-8, wound healing, flow cytometry, transwell and EdU assays were performed to assess the effects of TMEM40 on cell proliferation, cell cycle and apoptosis, migration and invasion. In addition, tumor growth was evaluated in vivo using a xenogenous subcutaneously implant model. All statistical analyses were executed by using the SPSS 20.0 software. All experimental data from three independent experiments were analyzed by Student's t test and results were expressed as mean ± standard deviation.

RESULTS

In this study, we identified the role of TMEM40 in the tumorigenesis of bladder cancer and found that it was upregulated in bladder cancer tissues and cell lines, compared with their normal counterparts. The results demonstrated that effective silence of TMEM40 expression suppressed cell proliferation, blocked G1-to-S cell cycle transition, and inhibited cell migration and invasion in human bladder 5637 and EJ cell lines. Consistently, in vivo data showed that TMEM40 silencing could dramatically decreased tumor growth. Further study revealed that TMEM40 knockdown resulted in accumulation of p53 and p21 protein and decrease of c-MYC and cyclin D1 protein.

CONCLUSION

These data suggest that TMEM40 represents a potential oncogene, which exert a crucial role in the proliferation and apoptosis via the p53 signaling pathway in BCa, thus probably serve as a novel candidate biomarker and a potential therapeutic target for patients with BCa.

Comparative analysis of TCDD-induced AhR-mediated gene expression in human, mouse and rat primary B cells

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a persistent environmental pollutant that activates the aryl hydrocarbon receptor (AhR) resulting in altered gene expression. In vivo, in vitro, and ex vivo studies have demonstrated that B cells are directly impaired by TCDD, and are a sensitive target as evidenced by suppression of antibody responses. The window of sensitivity to TCDD-induced suppression of IgM secretion among mouse, rat and human B cells is similar. Specifically, TCDD must be present within the initial 12h post B cell stimulation, indicating that TCDD disrupts early signaling network(s) necessary for B lymphocyte activation and differentiation. Therefore, we hypothesized that TCDD treatment across three different species (mouse, rat and human) triggers a conserved, B cell-specific mechanism that is involved in TCDD-induced immunosuppression. RNA sequencing (RNA-Seq) was used to identify B cell-specific orthologous genes that are differentially expressed in response to TCDD in primary mouse, rat and human B cells. Time course studies identified TCDD-elicited differential expression of 515 human, 2371 mouse and 712 rat orthologous genes over the 24-h period. 28 orthologs were differentially expressed in response to TCDD in all three species. Overrepresented pathways enriched in all three species included cytokine-cytokine receptor interaction, ECM-receptor interaction, focal adhesion, regulation of actin cytoskeleton and pathways in cancer. Differentially expressed genes functionally associated with cell-cell signaling in humans, immune response in mice, and oxidation reduction in rats. Overall, these results suggest that despite the conservation of the AhR and its signaling mechanism, TCDD elicits species-specific gene expression changes.

Fine-mapping QTLs in advanced intercross lines and other outbred populations

Quantitative genetic studies in model organisms, particularly in mice, have been extremely successful in identifying chromosomal regions that are associated with a wide variety of behavioral and other traits. However, it is now widely understood that identification of the underlying genes will be far more challenging. In the last few years, a variety of populations have been utilized in an effort to more finely map these chromosomal regions with the goal of identifying specific genes. The common property of these newer populations is that linkage disequilibrium spans relatively short distances, which permits fine-scale mapping resolution. This review focuses on advanced intercross lines (AILs) which are the simplest such population. As originally proposed in 1995 by Darvasi and Soller, an AIL is the product of intercrossing two inbred strains beyond the F2 generation. Unlike recombinant inbred strains, AILs are maintained as outbred populations; brother-sister matings are specifically avoided. Each generation of intercrossing beyond the F2 further degrades linkage disequilibrium between adjacent makers, which allows for fine-scale mapping of quantitative trait loci (QTLs). Advances in genotyping technology and techniques for the statistical analysis of AILs have permitted rapid advances in the application of AILs. We review some of the analytical issues and available software, including QTLRel, EMMA, EMMAX, GEMMA, TASSEL, GRAMMAR, WOMBAT, Mendel, and others.

Molecular signature of disease onset in granulin mutation carriers: a gene expression analysis study

Mutations within Granulin (GRN) gene are causative of autosomal dominant frontotemporal lobar degeneration (FTLD). Though GRN mutations are inherited at birth, the disease onset usually occurs in the sixth decade of life. The objective of this study was to identify new genetic pathways linked to inherited GRN disease and involved in the shift from asymptomatic to symptomatic stages. Microarray gene expression analysis on leukocytes was carried out on 15 patients carrying GRN T272SfsX10 mutation, and their asymptomatic siblings with (n = 14) or without (n = 11) GRN mutation. The results were then validated by real-time polymerase chain reaction, and compared with those obtained in a cohort of FTLD without GRN mutation (n = 16). The association between candidate genes and damage of specific brain areas was investigated by voxel-based morphometry on magnetic resonance imaging scans (family-wise error-corrected). Leukocytes mRNA levels of TMEM40 and LY6G6F and other genes mainly involved in inflammation were significantly higher in patients carrying GRN mutations compared with asymptomatic carriers and other FTLD. The higher the levels of TMEM40 the greater is the damage of parietal lobule; the higher the LY6G6F gene expression the greater is the atrophy in superior frontal gyrus. Enhanced inflammation associated with the onset of GRN disease might be either related to disease pathogenetic mechanism leading to neurodegeneration or to a compensatory pathway that counteracts disease progression. The identification of specific molecular targets of GRN-FTLD disease is essential when considering future disease-modifying therapies.

Non-MHC risk alleles in rheumatoid arthritis and in the syntenic chromosome regions of corresponding animal models

Rheumatoid arthritis (RA) is a polygenic autoimmune disease primarily affecting the synovial joints. Numerous animal models show similarities to RA in humans; some of them not only mimic the clinical phenotypes but also demonstrate the involvement of homologous genomic regions in RA. This paper compares corresponding non-MHC genomic regions identified in rodent and human genome-wide association studies (GWAS). To date, over 30 non-MHC RA-associated loci have been identified in humans, and over 100 arthritis-associated loci have been identified in rodent models of RA. The genomic regions associated with the disease are designated by the name(s) of the gene having the most frequent and consistent RA-associated SNPs or a function suggesting their involvement in inflammatory or autoimmune processes. Animal studies on rats and mice preferentially have used single sequence length polymorphism (SSLP) markers to identify disease-associated qualitative and quantitative trait loci (QTLs) in the genome of F2 hybrids of arthritis-susceptible and arthritis-resistant rodent strains. Mouse GWAS appear to be far ahead of rat studies, and significantly more mouse QTLs correspond to human RA risk alleles.

High-resolution mapping of a complex disease, a model for rheumatoid arthritis, using heterogeneous stock mice

Resolving the genetic basis of complex diseases like rheumatoid arthritis will require knowledge of the corresponding diseases in experimental animals to enable translational functional studies. Mapping of quantitative trait loci in mouse models of arthritis, such as collagen-induced arthritis (CIA), using F(2) crosses has been successful, but can resolve loci only to large chromosomal regions. Using an inbred-outbred cross design, we identified and fine-mapped CIA loci on a genome-wide scale. Heterogeneous stock mice were first intercrossed with an inbred strain, B10.Q, to introduce an arthritis permitting MHCII haplotype. Homozygous H2(q) mice were then selected to set up an F(3) generation with fixed major histocompatibility complex that was used for arthritis experiments. We identified 26 loci, 18 of which are novel, controlling arthritis traits such as incidence of disease, severity and time of onset and fine-mapped a number of previously mapped loci.

Refined mapping of the renal failure RF-3 quantitative trait locus

Rf-3, a quantitative trait locus (QTL) on rat chromosome 3, affects the development of CKD in Fawn-Hooded Hypertensive (FHH) rats. This QTL spans 110 Mb and approximately 1400 genes; therefore, narrowing the position of this locus is necessary to elucidate potential candidate genes. Here, we used congenic models and comparative genomics to refine the Rf-3 candidate region. We generated congenic lines carrying smaller intervals (subcongenics) of the Rf-3 region and used these lines to reduce the Rf-3 candidate region by 94% (to 7.1 Mb). We used comparative genomics to identify QTL for both nephropathy and albuminuria in the syntenic region of this interval for both human and mouse. We also used the overlapping homologous regions to reduce the number of likely positional candidate genes to 13 known or predicted genes. By combining congenic models and cross-species studies, we narrowed the list of candidate genes to a level that we could sequence the whole interval to further identify the causative gene in future studies.

Genome-wide association studies and the problem of relatedness among advanced intercross lines and other highly recombinant populations

Model organisms offer many advantages for the genetic analysis of complex traits. However, identification of specific genes is often hampered by a lack of recombination between the genomes of inbred progenitors. Recently, genome-wide association studies (GWAS) in humans have offered gene-level mapping resolution that is possible because of the large number of accumulated recombinations among unrelated human subjects. To obtain analogous improvements in mapping resolution in mice, we used a 34th generation advanced intercross line (AIL) derived from two inbred strains (SM/J and LG/J). We used simulations to show that familial relationships among subjects must be accounted for when analyzing these data; we then used a mixed model that included polygenic effects to address this problem in our own analysis. Using a combination of F(2) and AIL mice derived from the same inbred progenitors, we identified genome-wide significant, subcentimorgan loci that were associated with methamphetamine sensitivity, (e.g., chromosome 18; LOD = 10.5) and non-drug-induced locomotor activity (e.g., chromosome 8; LOD = 18.9). The 2-LOD support interval for the former locus contains no known genes while the latter contains only one gene (Csmd1). This approach is broadly applicable in terms of phenotypes and model organisms and allows GWAS to be performed in multigenerational crosses between and among inbred strains where familial relatedness is often unavoidable.