Skip to main content
Download PDF

Determining oncogenic patterns and cancer predisposition through the transcriptomic profile in Mitchell–Riley syndrome with heterotopic gastric mucosa and duodenal atresia: a case report

Orphanet Journal of Rare Diseases volume 16, Article number: 455 (2021) Cite this article

Abstract

Background

Homozygous mutations in the transcription factor RFX6 are the cause of the Mitchell–Riley syndrome (MRS) associating neonatal diabetes, congenital digestive system, such as biliary atresia, pancreatic hypoplasia, duodenal and/or jejunal atresia, intestinal malrotation, gallbladder aplasia, cholestasis. A constitutive inactivation of RFX6 leads also to gastric heterotopia. Application of RNA-seq in human diseases may help to better understand pathogenic mechanism of diseases and to predict the risk of developing chronic disorders and personalizing their prevention and treatment. We evaluated oncogenic patterns and cancer predisposition using the transcriptomic profile in a case of MRS with neonatal diabetes, duodenal atresia, and extensive intestinal tract gastric heterotopia.

Results

We signalled the interactors of RFX6 with other up and downregulated genes, that may be interested in severity of diabetic condition, in multi-organs impairment and cancer predisposition. Furthermore, several dysregulated genes are involved in biological processes that can lead to promote cancer including “Evading apoptosis” (BAD, BBC3, EGF, FGFR2, FLT3LG, HMOX1, HRAS, IFNAR2, IGF1R, IL12RB1, IL13RA1, IL15, IL2RB, IL2RG, IL6R, KEAP1, MGST1, PDGFA, PDGFRB, PIK3R3, RALB, RALGDS, RASSF1, SOS1, TGFA, TXNRD3), “Proliferation” (APC, BRAF, CCND2, CCND3, CCNE2, FGFR2, FLT3LG, FZD1, FZD6, HMOX1, HRAS, IGF1R, KEAP1, LRP6, MAPK3, MGST1, PDGFA, PDGFB, PDGFRB, RB1, SOS1, TGFA, TXNRD3, WNT10B), “Sustained angiogenesis” (BRAF, FGFR2, FLT3LG, HRAS, IGF1R, JAG1, MAPK3, NOTCH2, PDGFA, PDGFB, PDGFRB, SOS1, TGFA, TGFB1), “Genomic instability” (BAD, BBC3) and “Insensitivity to anti-growth signals” (SMAD2, TGFB1). We also inspected the signalings and their related genes in cancer, such as “PI3K signaling”, “ERK signaling”, “JAK-STAT signaling”, “Calcium signaling”, “Other RAS signaling”, “WNT signaling”.

Conclusions

In our MRS patient, we signaled the interactors of RFX6 with other up- and downregulated genes that may be related to severe diabetic condition, multi-organ impairment, and cancer predisposition. Notably, many dysregulated genes may lead to triggering carcinogenesis. The possibility of the patient developing cancer degeneration in heterotopic gastric mucosa and/or additional long-term tumoral sequelae is not excluded. Personalized prevention and treatment strategies should be proposed.

Background

Mitchell–Riley syndrome (MRS) is an autosomal recessive disorder caused by mutations in the RFX6 gene in which a combination of neonatal diabetes mellitus and congenital gastrointestinal defects—such as atresia, stenosis or malrotation of the small intestine, gallbladder hypoplasia or agenesis, intrahepatic or extrahepatic ductal atresia, or hypoplastic or annular pancreas—occur [1,2,3,4,5,6]. In a limited number of patients with MRS, heterotopic jejunal gastric mucosa, including in the duodenal and jejunal tract, has been described [2, 7, 8]. Even though the link between RFX6 and heterotopic gastric mucosa has not been extensively studied, as reported by Piccand et al. [9], RFX6 is required for the maintenance of intestinal cell identity and the constitutive inactivation of RFX6 leads to gastric heterotopia.

Heterotopic gastric mucosa is a pathological condition that has rarely been reported on and consists of ectopic gastric mucosa. It can be discovered anywhere throughout the gastrointestinal tract and may be asymptomatic or present with intussusception, obstruction, pain, bleeding, ulceration, or perforation. The association between intestinal tract gastric heterotopia and carcinogenesis is controversial [10,11,12,13,14,15]. Heterotopic gastric mucosa grow in the submucosa. Surface mucosa are vulnerable to repetitive erosion and regeneration, and these could be the cause of aberrations in the gastric mucosa and the development of cancer degeneration at the surface mucosa [16].

The pathogenesis of congenital anomalies (CAs), such as duodenal atresia associated with intestinal gastric heterotopia, has not been fully elucidated. Due to interactions between genes and the environment, during the organogenesis phase a modification of normal embryo-fetal development may occur [17]. The memory of the insult will be retained by the organism and may result in pathology later on, such as an increased risk of noncommunicable diseases (NCDs) and cancer predisposition [18, 19].

The recent application of scRNA-seq in human diseases may help us to better understand the pathogenic mechanisms of diseases [18]. The identification of gene expression patterns may be useful for predicting the risk of developing chronic diseases and help to personalize prevention and treatment methods [18,19,20,21,22].

We evaluated oncogenic patterns and cancer predisposition using the transcriptomic profile in a case of Mitchell–Riley syndrome with neonatal diabetes, duodenal atresia, and extensive intestinal tract gastric heterotopia, including the duodenum and the jejunum. Sequencing could aid in the implementation of personalized prevention and treatment strategies.

Results

The analysis of DEGs revealed that 4834 transcripts passed all the filters in MRS. Among them, 2202 transcripts had a fold change higher than 2, while 2632 had one lower than − 2. The mean of the distribution of the fold change in the upregulated genes is 4.35+/− 2.18, while that in the downregulated genes is − 4.33+/− 1.92.

MRS is known to be associated with mutations in the RFX6 gene. For this reason, we inspected the closest interactors of RFX6 deregulated in our analysis. We observed nine genes with upregulated fold changes (AGAP4, ATR, EHBP1, EIF2AK3, GLIS3, IFT88, RFX7, RPGRIP1L, SYTL4) and five genes with downregulated fold changes (GHRL, GPR68, IER3IP1, MAFB, PLAGL1), as reported in the fold change column of Table 1. Additionally, the table shows for each gene the transcripts’ mean counts for the control or MRS group obtained after DESeq2 normalization. Among these interactors, EIF2AK3, GLIS3, IER3IP1, and PLAGL1 are linked to Diabetes mellitus in Swiss-Prot. In addition, GPR68 is associated with tumor suppressor and MAFB is associated with proto-oncogene and tumor suppressor activities. In Figure 1 the 14 DEG interactors of RFX6 are plotted. Notably, EIF2AK3, IER3IP1, GLIS3, RFX7, and IFT88 were present at the highest degree (3).

Table 1. DEGs interactors of RFX6 in our analysis
Full size table
Figure 1.
figure 1

Interactors of RFX6 that are deregulated in our analysis obtained with STRING. Ignoring the construction RFX6, the genes EIF2AK3, IER3IP1, GLIS3, RFX7 and IFT88 have 3 as the degree that is the highest value in the network

Full size image

Among the DEGs highlighted in our analysis, 100 were characterized in several processes in the KEGG map “pathways in cancer” (Figure 2). We observed (Table 2) that the major implications occurred for “evading apoptosis” (BAD, BBC3, EGF, FGFR2, FLT3LG, HMOX1, HRAS, IFNAR2, IGF1R, IL12RB1, IL13RA1, IL15, IL2RB, IL2RG, IL6R, KEAP1, MGST1, PDGFA, PDGFRB, PIK3R3, RALB, RALGDS, RASSF1, SOS1, TGFA, TXNRD3), “proliferation” (APC, BRAF, CCND2, CCND3, CCNE2, FGFR2, FLT3LG, FZD1, FZD6, HMOX1, HRAS, IGF1R, KEAP1, LRP6, MAPK3, MGST1, PDGFA, PDGFB, PDGFRB, RB1, SOS1, TGFA, TXNRD3, WNT10B), “sustained angiogenesis” (BRAF, FGFR2, FLT3LG, HRAS, IGF1R, JAG1, MAPK3, NOTCH2, PDGFA, PDGFB, PDGFRB, SOS1, TGFA, TGFB1), “genomic instability” (BAD, BBC3) and “insensitivity to anti-growth signals” (SMAD2, TGFB1).

Figure 2.
figure 2

Number of DEGs in the processes identified inside the KEGG map “Pathways in cancer”. The red line represent the median value so that “Evading apoptosis”, “Proliferation”, “Sustained angiogenesis”, “Genomic instability” and “Insensitivity to anti-growth signals” have a major implication in our study

Full size image
Table 2. DEGs inside the KEGG map “Pathways in cancer” identified in the most implicated processes
Full size table

In addition to the biological implications, we also inspected the signals and their related genes in cancer (Figure 3). The signals in which the most of the genes were implicated (Table 3) are “PI3K signaling” (BAD, EGF, FGFR2, FLT3LG, HRAS, IGF1R, PDGFA, PDGFRB, PIK3R3, SOS1, TGFA), “ERK signaling” (BRAF, FGFR2, FLT3LG, HRAS, IGF1R, MAPK3, PDGFA, PDGFB, PDGFRB, SOS1, TGFA), “JAK-STAT signaling” (EGF, FLT3LG, IFNAR2, IL12RB1, IL13RA1, IL15, IL2RB, IL2RG, IL6R), “calcium signaling” (CAMK2D, EGF, PDGFA, PDGFB, PDGFRB, PRKCG, TGFA), “other RAS signaling” (EGF, HRAS, RALB, RALGDS, RASSF1, SOS1), and “WNT signaling” (APC, FZD1, FZD6, LRP6, WNT10B).

Figure 3.
figure 3

Number of DEGs in the signaling identified inside the KEGG map “Pathways in cancer”. The red line represents the median value so that “ERK signaling”, “PI3K signaling”, “JAK-STAT signaling”, “Calcium signaling”, “Other RAS signaling” and “WNT signaling” have a major implication in our study

Full size image
Table 3. DEGs inside the KEGG map “Pathways in cancer” identified in the most implicated signalings
Full size table

Discussion

Homozygous mutations in the transcription factor (TF) RFX6 are the cause of the MRS associating neonatal diabetes congenital digestive system defects, including biliary atresia, pancreatic hypoplasia, duodenal and/or jejunal atresia, intestinal malrotation, gallbladder aplasia, and cholestasis. In our case, heterotopic gastric mucosa in the small bowel tissue is also reported [2, 7, 8].

RFX6 is included in the TF regulatory network of human pancreas development. During early pancreas specification and lineage commitment, specific TFs and other critical markers are expressed at each stage [9, 23]. Mutations in RFX6 are involved in endocrine and exocrine pancreatic insufficiency and also in the altered maturation of the enteroendocrine cell subpopulation in the gastrointestinal tract, leading to diabetes and severe malabsorption [23]. Additionally, the RFX6 gene controls genetic programs, regulating Peptidergic Enteroendocrine cell differentiation and identity. The constitutive inactivation of RFX6 leads to a lack of functional compensation in some pluripotent endodermal stem cells, which consequently lose or do not acquire their proper regional identity, resulting in patchy patterns in gastric tissue [9].

Our network analysis highlighted AGAP4, ATR, EHBP1, EIF2AK3, GHRL, GLIS3, GPR68, IER3IP1, IFT88, MAFB, PLAGL1, RFX7, RPGRIP1L, and SYTL4 as interactors of RFX6 (Table 1). Among them, EIF2AK3, IER3IP1, GLIS3, RFX7, and IFT88 were present to the highest degree (Fig. 1).

The interaction of RFX6 with other upregulated genes, such as EIF2AK3, may induce a severe diabetic condition and could be related to multi-organ impairment and cancer degeneration. In fact, the EIF2AK3 gene is also involved in the mechanism by which endoplasmic reticulum (ER) stress leads to programmed cell death (PCD) [24, 25]. PCD is an essential part of the maintenance of cellular homeostasis and survival of multicellular organisms during embryologic development, after birth, and during adulthood. Cell death is suspected to be also partially responsible for degenerative diseases and the uncontrolled proliferation of cancer cells [24, 25]. In heterotopic gastric mucosa, PCD and the vulnerability of surface mucosa to repetitive erosion and regeneration could be the cause of aberrations in the gastric mucosa and cancer degeneration at the surface mucosa [16, 26].

To date, in the literature less than twenty patients with MRS have been described and long-term follow ups of these patients are limited. Considering the reported risk of NCDs and cancer in patients with chromosomal CAs [18], the analysis of the transcriptomic profile may be useful to discover genetic alteration and support links to the disease. It is noteworthy that the analysis can be particularly useful when the individual has a known cancer predisposition, as in our patient who presented with heterotopic gastric mucosa.

EIF2AK3 also interacts with IER3IP1 and GLIS3. IER3IP1 encodes the immediate early response 3 interacting protein 1 into the ER. It also participates in permanent neonatal diabetes mellitus and seems to alter the death and proliferation rate of the β-cells that potentially occur after ER stress [27]. Sun, J, and Ren, D., observed an alteration in IER3IP1 expression with apoptotic proteins of the BCL-2 family [28], while the PI3K/AKT/mTOR pathway seems likely to play a key role in β-cell growth, proliferation, and survival [29]. Interestingly, in the PI3K signaling investigated in our study, we observed the downregulation of the BCL2-associated agonist of cell death encoded by the BAD gene along with the deregulation of EGF, FGFR2, FLT3LG, HRAS, IGF1R, PDGFA, PDGFRB, PIK3R3, SOS1, and TGFA. On the other hand, GLIS3 belongs to the GLIS protein family of zinc finger and occurs in many diseases, including neonatal diabetes mellitus, aplasia, hypoplasia, hypothyroidism, growth retardation, atrial septal defects, autoimmune diseases, and neurological disorders. GLIS3 deregulation was also associated with liver, thyroid, and breast cancer. In particular, breast carcinogenesis seems to take advantage of the WNT/B-catenin pathway [30]. The WNT signaling in our study is altered by the deregulation of APC, FZD1, FZD6, LRP6, and WNT10B, which are included by KEGG in the process of cancer proliferation.

RFX6 also interacts with IFT88 and RFX7, which in our network are present to a high degree and, in turn, interact each other.

Intraflagellar Transport 88, encoded by IFT88, is involved in cilium biogenesis and the genetic mutations carried on this gene were firstly associated with kidney disease. Then, cilia dysfunction was also linked to diabetes and cancer. Recently, IFT88 was also suggested to play a role in the primary cilium in Hedgehog and WNT signaling [31]. Interestingly, its interactor RPGRIP1L, which encodes for RPGRIP1 Like, is localized in ciliated cells and seems to regulate the activity of the ciliary proteosome, which was observed to be altered in many cancers [32]. Contrary to other members of Regulatory Factor X family, RFX7 is poorly characterized but shows a high tumor suppressor potential in lymphoid cancers at least [33, 34]. Furthermore, RFX7 is upregulated in breast cancer but negatively correlates with metastatic development [35]. Additionally, RFX6 and RFX7 directly interact with ATR, encoding for ATR Serine/Threonine Kinase. ATR Serine/Threonine Kinase, along with its downstream Checkpoint Kinase 1, plays a crucial role in the regulation of the cell cycle for repairing DNA in response to damage. Similar to our analysis, other studies have shown that the activation of ATR leads to cell survival and proliferation; thus, several inhibitors of ATR/CHK1 have been proposed to treat cancer [28, 36].

Similar to GLIS3, the PLAGL1 gene encodes for a zinc finger protein associated with diabetes and has been nominated as a tumor suppressor. It can handle the cell cycle as well as apoptosis and was found to be downregulated in many cancers, such as prostate, colon, ovarian, and breast [37]. The downregulation of PLAGL1 in our analysis is in line with this observation, so its role in the cancer degeneration of heterotopic gastric areas cannot be excluded.

MAFB and GPR68 genes are downregulated interactors of RFX6 that Swiss-Prot associates with cancer. The MAF BZIP Transcription Factor B, encoded by MAFB, contributes to the differentiation of pancreatic α- and β-cells and to adult islet function. Nevertheless, it expresses proto-oncogene or tumor suppressor potential depending on the cell context. Lu et al. showed that MafB could play a pivot role in the proliferation of β-cells as well as in tumorigenesis condition [38]. The tumor suppressor gene GPR68, known as ovarian cancer G protein-coupled receptor 1 (OGR1), is a transmembrane receptor of the proton-sensing G protein-coupled receptors that is activated when extracellular levels of PH are altered. RFX6 regulates the transcription of GPR68 in adult human β-cells, in turn promoting the production of inflammatory interleukin 8 [39]. Herein, we observed the downregulation of IL15 along with the receptors IL2RB, IL2RG, IL6R, IL12RB1, and IL13RA1 in JAK-STAT signaling. In cancer, this culminates with the evading apoptosis event, as reported in KEGG. Additionally, the signal could be also implicated in the development of eosinophilic colitis. The GHRL gene, encoding ghrelin and obestatin prepropeptide, is associated with insulin secretion and downregulated. Our research group has already observed in a previous work that an infant with VACTERL and esophageal atresia carried a missense mutation on GHRL linked to metabolic syndrome [40]. In addition, the role of GHRL seems to be crucial in cell proliferation, migration, and invasion, as well as in inflammation in many cancers. Even so, its role as a promoter or inhibitor is still up for debate [41]. SYTL4 encodes for the Synaptogamin Like 4, which is localized on the microtubule cytoskeleton. Liu et al. showed that Synaptogamin Like 4 hinders microtubule polymerization, reducing the stability. Furthermore, high levels of SYTL4, as in our case study, are associated with a poor prognosis for breast cancer [42]. EHBP1 encodes EH Domain Binding Protein 1 and is poorly characterized. Nevertheless, EHBP1 has already been associated with cancer and is required for the insulin-mediated translocation of glucose transporter type 4 [43, 44]. AGAP4 encodes for ArfGAP with GTPase domain and Ankyrin Repeat and PH Domain 4, but very little is known about it. It is part of the family of centaurins as proteins with a GTPase-like domain which are known to regulate cell proliferation and vesicular trafficking. For similarity, GTPase-like centaurin γ − 1 is associated with cancer, in which it promotes cell invasion and prevents apoptosis [45].

In our analysis, we also wanted to further inspect the different hallmarks of cancer, as well as which of them could allow its development [46]. We categorized each process based on the amount of DEGs involved (Table 2). Thus, we found that “evading apoptosis”, “proliferation”, “sustained angiogenesis”, “genomic instability”, and “insensitivity to anti-growth signals” are the most relevant adopted strategies (Fig. 2). Noteworthy, research on cancer has revealed many pathways that are able to promote these processes (Fig. 3). For this reason, we inspected the key signals in MSR and we observed that “ERK signaling”, “PI3K signaling”, “JAK-STAT signaling”, “calcium signaling”, “other RAS signaling”, and “WNT signaling” show the highest involvement (Table 3). Notably, the PI3K, Wnt/β-catenin, and RAS/ERK signaling pathways can handle cellular metabolism and consequently influence signal transduction and oxidative stress potential [47]. Furthermore, cancer proliferation and invasiveness were reported to be caused by altered calcium signaling in the tumor microenvironment [48]. The regulation of the microenvironment can also be changed through the alteration of Jak-Stat signaling. This signaling is mediated by inflammatory cytokines that promote the self-renewal of cancer stem cells and differentiation [48, 49].

Conclusions

MRS is caused by the mutation of the RFX6 gene and is characterized by neonatal diabetes, pancreatic hypoplasia, intestinal atresia, and gallbladder hypoplasia or aplasia and chronic diarrhea. A constitutive inactivation of RFX6 may also lead to gastric heterotopia. In our MRS patient, we signaled the interactors of RFX6 with other up- and downregulated genes that may be related to severe diabetic condition, multi-organ impairment, and cancer predisposition. Notably, many dysregulated genes take place in the mechanisms of evading apoptosis, proliferation, sustained angiogenesis, genomic instability, and insensitivity to anti-growth signals, which may lead to triggering carcinogenesis. The possibility of the patient developing cancer degeneration in heterotopic gastric mucosa and/or additional long-term tumoral sequelae is not excluded. Personalized prevention and treatment strategies should be proposed.

Methods

Patient

The patient is a male infant born to consanguineous Pakistani parents at 37 weeks+2 days of gestation, with intrauterine growth restriction (weight 1417 g, <3rd percentile; length 41 cm, < 3rd percentile; head circumference 34 cm, 75th percentile). The mother was diagnosed with hyperthyroidism during pregnancy, with a normal glucose profile. The father is diabetic. A prenatal suspicion of duodenal atresia was posed.

He scored an APGAR of 1 and 8 at 1 and 5 minutes, respectively, requiring non-invasive ventilation during the first 24 hours. At birth, the patent foramen ovale and duodenal atresia were detected. The baby underwent surgical repair of the type III duodenal atresia. At operation, gallbladder malposition, hypoplasia, an ectopic pancreas, and duodenal heterotopic gastric mucosa were found. Histologic evaluation confirmed the duodenal presence of gastric mucosa without signs of erosion.

The infant started insulin infusion for hyperglycemia from day 1; during a hyperglycemic episode, his C-peptide level was very low (< 0.1 ng/mL, nv 0.8–4.2 ng/mL).

From day 8, watery diarrhea became evident and persistent, associated with failure to thrive and dependency on parenteral nutrition.

The infant underwent genetic analysis, which detected homozygous missense mutations not reported previously in the RFX6 [p.Ser500Gly (c.1498A>G)], and MRS was diagnosed. A heterozygous mutation was detected in the parents.

Neonatal diabetes was confirmed and insulin pump therapy and continuous glucose monitoring were started on day 40.

At the age of 2 months, bilateral hyperechogenic renal parenchyma with normal renal function was detected.

During monitoring, protracted diarrhea persisted despite several interventions, including dietary adaptations with semi-elemental, elemental, and low long-chain triglyceride formulas and the use of pancreatic enzymes. Recurrent hospitalization for sepsis was recorded. Metabolic control of diabetes was near-optimal.

Stool was intermittently positive for microscopic blood, and iron substitution was required for chronic anemia. At the age of 8 months, the anemia worsened and repeated blood transfusions were required. Endoscopic investigation confirmed duodenal ectopic gastric mucosa including the entire duodenal surface. An extensive jejunal gastric heterotopia was also recorded, and histologic evaluation confirmed the jejunal heterotopic gastric mucosa with signs of intestinal mucosa erosion. An eosinophilic colitis was also confirmed. Total parenteral nutrition was introduced. Progressively, the boy started to improve clinically and we began a reintroduction of elemental, high medium-chain triglyceride formula by the mouth, with the use of pancreatic enzymes and multiple vitamin supplements associated with parental nutrition. Diarrhea persisted and consisted of 5 to 8 watery stools per day. At 13 months, antihypertensive treatment was started.

Currently, the patient is 15 months years old, 71.5 cm in length (< 3rd percentile), and with a 7835 kg body weight (< 3rd percentile). He shows an unstable control of diabetes, with high glycemic variability and a HbA1c 7.3% (56.3 mmol/mol). In addition to personalized parental nutrition, a hypoallergenic diet (using an elemental formula in milk) with pancreatic enzymes and multiple vitamin supplementation was adopted.

Transcriptomic analysis

The use of a control group is necessary in order to inspect the genes differentially expressed in MRS. From the Sequence Read Archive (SRA) [51], we selected the runs of the healthy samples GSM2370017, GSM2370185, GSM2370217, GSM2370225, GSM2370231, GSM2370237, GSM2370251, GSM2370261, GSM2370269, and GSM2370271 that belong to the bioproject PRJNA352062 [50].

MRS and control samples were analyzed, taking advantage of the same workflow. Fastq raw data were inspected with fastQC in order to analyze their quality. We toke advantage of Trimmomatic (version 0.38, Usadel Lab, Aachen, Germany) [51] to drop the poor-quality bases and potential adapters. Then, the Spliced Transcripts Alignment to a Reference (STAR) RNA-seq aligner [52] was used to align and sort the reads against the reference homo sapiens GRCh38 genome. The python package htseq-count counted the number of transcripts in each region [53]. The package DESeq2 of Bioconductor found the genes that were differentially expressed (DEGs) between MRS and the control group using the R programming language [54]. The Benjamini–Hochberg post hoc test was used to remove the false positives to correct the p-value. All the genes whose q-value was lower than 0.05 and whose fold change was lower than − 2.0 or over 2.0 were kept.

To highlight oncogene and tumor suppressor genes in our analysis, we observed the DEGs that take place in the KEGG [55] map “pathways in cancer” from homo sapiens (hsa05200). Indeed, this map collects the different signaling pathways that are activated in cancer and the different biological processes that they trigger. Specifically, we observed which DEGs were included in the map, which biological processes they promote, and in which pathways they play a role. Furthermore, in order to study the behavior of the proteins that are known to interact with the RFX6 gene, we used the STRING [56] database. Thus, we input our DEGs along with RFX6 and kept only the genes that transcribed proteins known to be direct interactors. Finally, we took advantage of the manual curated Swiss-Prot [57] database to provide a role of each DEG identified by STRING as an RFX6 interactor. In detail, the Swiss-Prot “Keyword – Disease” section was inspected and the involvement of the interactors with unhealthy status was further explored in the literature.

Availability of data and materials

Not applicable.

Abbreviations

MRS:

Mitchell–Riley syndrome

CAs:

Congenital anomalies

NCDs:

Noncommunicable diseases

PCD:

Programmed cell death

SRA:

Sequence read archive

DEGs:

Differentially expressed genes

References

  1. Kambal MA, Al-Harbi DA, Al-Sunaid AR, Al-Atawi MS. Mitchell–Riley syndrome due to a novel mutation in RFX6. Front Pediatr. 2019;7:243.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Sansbury FH, Kirel B, Caswell R, Allen HL, Flanagan SE, Hattersley AT, et al. Biallelic RFX6 mutations can cause childhood as well as neonatal onset diabetes mellitus. Eur J Hum Genet. 2015;23(12):1744–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Khan N, Dandan W, Al Hassani N, Hadi S. A Newly-Discovered Mutation in the RFX6 Gene of the Rare Mitchell–Riley Syndrome. J Clin Res Pediatr Endocrinol. 2016;8(2):246–9.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Zegre Amorim M, Houghton JA, Carmo S, Salva I, Pita A, Pereira-da-Silva L. Mitchell–Riley syndrome: a novel mutation in RFX6 Gene. Case Rep Genet. 2015;2015:937.

    Google Scholar 

  5. Concepcion JP, Reh CS, Daniels M, Liu X, Paz VP, Ye H, et al. Neonatal diabetes, gallbladder agenesis, duodenal atresia, and intestinal malrotation caused by a novel homozygous mutation in RFX6. Pediatr Diabetes. 2014;15(1):67–72.

    Article  CAS  PubMed  Google Scholar 

  6. Pearl EJ, Jarikji Z, Horb ME. Functional analysis of Rfx6 and mutant variants associated with neonatal diabetes. Dev Biol. 2011;351(1):135–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Meng X, Gong CX. Mitchell–Riley syndrome with late onset diabetes and heterotopic gastric mucosa: a case report. Chin J Pediatr. 2020;58(1):62–4.

    CAS  Google Scholar 

  8. Skopkova M, Ciljakova M, Havlicekova Z, Vojtkova J, Valentinova L, Danis D, et al. Two novel RFX6 variants in siblings with Mitchell–Riley syndrome with later diabetes onset and heterotopic gastric mucosa. Eur J Med Genet. 2016;59(9):429–35.

    Article  PubMed  Google Scholar 

  9. Piccand J, Vagne C, Blot F, Meunier A, Beucher A, Strasser P, et al. Rfx6 promotes the differentiation of peptide-secreting enteroendocrine cells while repressing genetic programs controlling serotonin production. Mol Metabol. 2019;29:24–39.

    Article  CAS  Google Scholar 

  10. Lemaire J, Delaunoit T, Molle G. Adenocarcinoma arising in gastric heterotopic pancreas Case report and review of the literature. Acta Chir Bel. 2014;114(1):79–81.

    CAS  Google Scholar 

  11. Neumann WL, Lujan GM, Genta RM. Gastric heterotopia in the proximal oesophagus (“inlet patch”): association with adenocarcinomas arising in Barrett mucosa. Dig Liver Dis Off J Ital Soc Gastroenterol Ital Assoc Study Liver. 2012;44(4):292–6.

    Google Scholar 

  12. Kosugi S, Kanda T, Hatakeyama K. Adenocarcinoma arising from heterotopic gastric mucosa in the stomach. J Gastroenterol Hepatol. 2006;21(2):483–4.

    Article  PubMed  Google Scholar 

  13. von Rahden BH, Stein HJ, Becker K, Liebermann-Meffert D, Siewert JR. Heterotopic gastric mucosa of the esophagus: literature-review and proposal of a clinicopathologic classification. Am J Gastroenterol. 2004;99(3):543–51.

    Article  Google Scholar 

  14. Pai S, Deshpande R, Naresh KN. Adenocarcinoma of cervical esophagus arising in aberrant gastric mucosa. Indian J Gastroenterol Off J Indian Soc Gastroenterol. 1997;16(4):157–8.

    CAS  Google Scholar 

  15. Yamagiwa H, Matsuzaki O, Ishihara A, Yoshimura H. Heterotopic gastric glands in the submucosa of the stomach. Acta Pathol Jpn. 1979;29(3):347–50.

    CAS  PubMed  Google Scholar 

  16. Kim DH, Kim KM, Oh SJ, Oh JA, Choi MG, Noh JH, et al. Early gastric cancer arising from heterotopic gastric mucosa in the gastric submucosa. J Korean Surg Soc. 2011;80(Suppl 1):S6–11.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Baldacci S, Gorini F, Santoro M, Pierini A, Minichilli F, Bianchi F. Environmental and individual exposure and the risk of congenital anomalies: a review of recent epidemiological evidence. Epidemiol Prev. 2018;42(3–4):1–34.

    PubMed  Google Scholar 

  18. Norwood MS, Lupo PJ, Chow EJ, Scheurer ME, Plon SE, Danysh HE, et al. Childhood cancer risk in those with chromosomal and non-chromosomal congenital anomalies in Washington State: 1984–2013. PLoS ONE. 2017;12(6):e0179006.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. Charles MA, Delpierre C, Breant B. Developmental origin of health and adult diseases (DOHaD): evolution of a concept over three decades. Med Sci. 2016;32(1):15–20.

    Google Scholar 

  20. Kim D, Kobayashi T, Voisin B, Jo JH, Sakamoto K, Jin SP, et al. Targeted therapy guided by single-cell transcriptomic analysis in drug-induced hypersensitivity syndrome: a case report. Nat Med. 2020;26(2):236–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Sandhu C, Qureshi A, Emili A. Panomics for precision medicine. Trends Mol Med. 2018;24(1):85–101.

    Article  PubMed  Google Scholar 

  22. Manzoni C, Kia DA, Vandrovcova J, Hardy J, Wood NW, Lewis PA, et al. Genome, transcriptome and proteome: the rise of omics data and their integration in biomedical sciences. Brief Bioinf. 2018;19(2):286–302.

    Article  CAS  Google Scholar 

  23. Aigha II, Abdelalim EM. Transcription factor: a crucial regulator of pancreatic beta cell development identity and proliferation. Stem Cell Res Ther. 2020;11(1):459.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Oyadomari S, Araki E, Mori M. Endoplasmic reticulum stress-mediated apoptosis in pancreatic beta-cells. Apop Int J Program Cell Death. 2002;7(4):335–45.

    Article  CAS  Google Scholar 

  25. Ghemrawi R, Khair M. Endoplasmic reticulum stress and unfolded protein response in neurodegenerative diseases. Int J Mol Sci. 2020;21(17).

  26. Sun J, Ren D. IER3IP1 deficiency leads to increased beta-cell death and decreased beta-cell proliferation. Oncotarget. 2017;8(34):56768–79.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Jahnavi S, Poovazhagi V, Kanthimathi S, Gayathri V, Mohan V, Radha V. EIF2AK3 mutations in South Indian children with permanent neonatal diabetes mellitus associated with Wolcott-Rallison syndrome. Pediatr Diabetes. 2014;15(4):313–8.

    Article  CAS  PubMed  Google Scholar 

  28. Li Z, Pearlman AH, Hsieh P. DNA mismatch repair and the DNA damage response. DNA Repair. 2016;38:94–101.

    Article  PubMed  CAS  Google Scholar 

  29. Balcazar-Morales N, de Plata CA. Role of AKT/mTORC1 pathway in pancreatic beta-cell proliferation. Colomb Med. 2012;43(3):235–43.

    Article  Google Scholar 

  30. Pinto K, Chetty R. Gene of the month: GLIS1-3. J Clin Pathol. 2020;73(9):527–30.

    Article  CAS  PubMed  Google Scholar 

  31. Veland IR, Awan A, Pedersen LB, Yoder BK, Christensen ST. Primary cilia and signaling pathways in mammalian development, health and disease. Nephron Physiol. 2009;111(3):p39-53.

    Article  CAS  PubMed  Google Scholar 

  32. Gerhardt C, Leu T, Lier JM, Ruther U. The cilia-regulated proteasome and its role in the development of ciliopathies and cancer. Cilia. 2016;5:14.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Fischer BA, Chelbi ST, Guarda G. Regulatory Factor X 7 and its Potential Link to Lymphoid Cancers. Trends Cancer. 2020;6(1):6–9.

    Article  CAS  PubMed  Google Scholar 

  34. Weber J, de la Rosa J, Grove CS, Schick M, Rad L, Baranov O, et al. PiggyBac transposon tools for recessive screening identify B-cell lymphoma drivers in mice. Nat Commun. 2019;10(1):1415.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. Yau C, Esserman L, Moore DH, Waldman F, Sninsky J, Benz CC. A multigene predictor of metastatic outcome in early stage hormone receptor-negative and triple-negative breast cancer. Breast Cancer Res. 2010;12(5):R85.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Qiu Z, Oleinick NL, Zhang J. ATR/CHK1 inhibitors and cancer therapy. Radiother Oncol J Eur Soc Ther Radiol Oncol. 2018;126(3):450–64.

    Article  CAS  Google Scholar 

  37. Vega-Benedetti AF, Saucedo C, Zavattari P, Vanni R, Zugaza JL, Parada LA. PLAGL1: an important player in diverse pathological processes. J Appl Genet. 2017;58(1):71–8.

    Article  CAS  PubMed  Google Scholar 

  38. Lu J, Hamze Z, Bonnavion R, Herath N, Pouponnot C, Assade F, et al. Reexpression of oncoprotein MafB in proliferative beta-cells and Men1 insulinomas in mouse. Oncogene. 2012;31(31):3647–54.

    Article  CAS  PubMed  Google Scholar 

  39. Chandra V, Karamitri A, Richards P, Cormier F, Ramond C, Jockers R, et al. Extracellular acidification stimulates GPR68 mediated IL-8 production in human pancreatic beta cells. Sci Rep. 2016;6:25765.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Pelizzo G, Chiricosta L, Mazzon E, Zuccotti GV, Avanzini MA, Croce S, et al. Discovering genotype variants in an infant with VACTERL through clinical exome sequencing: a support for personalized risk assessment and disease prevention. Pediatr Rep. 2021;13(1):45–56.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Lin TC, Hsiao M. Ghrelin and cancer progression. Biochim Biophys Acta Rev Cancer. 2017;1868(1):51–7.

    Article  CAS  PubMed  Google Scholar 

  42. Liu XY, Jiang W, Ma D, Ge LP, Yang YS, Gou ZC, et al. SYTL4 downregulates microtubule stability and confers paclitaxel resistance in triple-negative breast cancer. Theranostics. 2020;10(24):10940–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Ghalali A, Wiklund F, Zheng H, Stenius U, Hogberg J. Atorvastatin prevents ATP-driven invasiveness via P2X7 and EHBP1 signaling in PTEN-expressing prostate cancer cells. Carcinogenesis. 2014;35(7):1547–55.

    Article  CAS  PubMed  Google Scholar 

  44. Guilherme A, Soriano NA, Furcinitti PS, Czech MP. Role of EHD1 and EHBP1 in perinuclear sorting and insulin-regulated GLUT4 recycling in 3T3-L1 adipocytes. J Biol Chem. 2004;279(38):40062–75.

    Article  CAS  PubMed  Google Scholar 

  45. Soundararajan M, Yang X, Elkins JM, Sobott F, Doyle DA. The centaurin gamma-1 GTPase-like domain functions as an NTPase. Biochem J. 2007;401(3):679–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.

    Article  CAS  PubMed  Google Scholar 

  47. Park JH, Pyun WY, Park HW. Cancer metabolism: phenotype, signaling and therapeutic targets. Cells. 2020;9(10).

  48. Bong AHL, Monteith GR. Calcium signaling and the therapeutic targeting of cancer cells. Biochim Biophys Acta Mol Cell Res. 2018;1865(11):1786–94.

    Article  CAS  PubMed  Google Scholar 

  49. Groner B, von Manstein V. Jak Stat signaling and cancer: opportunities, benefits and side effects of targeted inhibition. Mol Cell Endocrinol. 2017;451:1–14.

    Article  CAS  PubMed  Google Scholar 

  50. Thompson EG, Du Y, Malherbe ST, Shankar S, Braun J, Valvo J, et al. Host blood RNA signatures predict the outcome of tuberculosis treatment. Tuberculosis. 2017;107:48–58.

    Article  CAS  PubMed  Google Scholar 

  51. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21.

    Article  CAS  PubMed  Google Scholar 

  53. Anders S, Pyl PT, Huber W. HTSeq–a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9.

    Article  CAS  PubMed  Google Scholar 

  54. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  55. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.

    Article  CAS  PubMed  Google Scholar 

  57. UniProt C. UniProt: the universal protein knowledgebase in 2021. Nucleic Acids Res. 2021;49(D1):D480–9.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was supported by Current Research Fund 2020, Ministry of Health, Italy.

Author information

Author notes

    Authors and Affiliations

    1. Pediatrics and Adolescentology Unit, Department of Internal Medicine, University of Pavia, Pavia, Italy

      Valeria Calcaterra

    2. Pediatric Department, Children’s Hospital “Vittore Buzzi”, Milan, Italy

      Valeria Calcaterra, Elvira Verduci, Dario Dilillo & Gianvincenzo Zuccotti

    3. IRCCS Centro Neurolesi “Bonino-Pulejo”, Messina, Italy

      Luigi Chiricosta, Emanuela Mazzon & Agnese Gugnandolo

    4. Pediatric Surgery Department, “Spedali Civili” Children’s Hospital, Brescia, Italy

      Daniele Alberti

    5. Pediatric Surgery Department, Children’s Hospital “Vittore Buzzi”, Via Lodovico Castelvetro n.32, 20154, Milan, Italy

      Luciano Maestri, Milena Meroni, Elettra Vestri & Gloria Pelizzo

    6. Department of Health Sciences, University of Milano, Milan, Italy

      Elvira Verduci

    7. Department of Biomedical and Clinical Science “L. Sacco”, University of Milano, Milan, Italy

      Gianvincenzo Zuccotti & Gloria Pelizzo

    Authors
    1. Valeria Calcaterra
      View author publications

      You can also search for this author in PubMed Google Scholar

    2. Luigi Chiricosta
      View author publications

      You can also search for this author in PubMed Google Scholar

    3. Emanuela Mazzon
      View author publications

      You can also search for this author in PubMed Google Scholar

    4. Agnese Gugnandolo
      View author publications

      You can also search for this author in PubMed Google Scholar

    5. Daniele Alberti
      View author publications

      You can also search for this author in PubMed Google Scholar

    6. Luciano Maestri
      View author publications

      You can also search for this author in PubMed Google Scholar

    7. Milena Meroni
      View author publications

      You can also search for this author in PubMed Google Scholar

    8. Elettra Vestri
      View author publications

      You can also search for this author in PubMed Google Scholar

    9. Elvira Verduci
      View author publications

      You can also search for this author in PubMed Google Scholar

    10. Dario Dilillo
      View author publications

      You can also search for this author in PubMed Google Scholar

    11. Gianvincenzo Zuccotti
      View author publications

      You can also search for this author in PubMed Google Scholar

    12. Gloria Pelizzo
      View author publications

      You can also search for this author in PubMed Google Scholar

    Contributions

    Conceptualization, V.C., L.C., E.M., G.P.; methodology and formal analysis, L.C., G.A., E.M.; software L.C.; writing—original draft preparation V.C., L.C., E.M., G.A., L.M, M.M., E.V., E.V., D.D., and G.P.; writing—review and editing, V.C., L.C., E.M., G.Z., G.P. All authors have read and agreed to the published version of the manuscript.

    Corresponding author

    Correspondence to Gloria Pelizzo.

    Ethics declarations

    Consent to participate

    The parents of the patient provided informed consent for genetic testing and publication. No images or videos are used.

    Competing interests

    The authors declare no conflict of interest.

    Additional information

    Publisher's Note

    Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

    Valeria Calcaterra and Luigi Chiricosta are contributed equally

    Rights and permissions

    Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

    Reprints and permissions

    About this article

    Check for updates. Verify currency and authenticity via CrossMark

    Cite this article

    Calcaterra, V., Chiricosta, L., Mazzon, E. et al. Determining oncogenic patterns and cancer predisposition through the transcriptomic profile in Mitchell–Riley syndrome with heterotopic gastric mucosa and duodenal atresia: a case report. Orphanet J Rare Dis 16, 455 (2021). https://doi.org/10.1186/s13023-021-02093-9

    Download citation

    • Received:

    • Accepted:

    • Published:

    • DOI: https://doi.org/10.1186/s13023-021-02093-9

    Keywords

    Orphanet Journal of Rare Diseases

    ISSN: 1750-1172