Skip to main content

MicroRNA profiling of primary pulmonary enteric adenocarcinoma in members from the same family reveals some similarities to pancreatic adenocarcinoma—a step towards personalized therapy

Abstract

Background

Primary pulmonary enteric adenocarcinoma (PEAC) is defined as a pulmonary adenocarcinoma with a predominant component of intestinal differentiation and tumor cells positive for at least one intestinal marker. The aim of the present study was the molecular and histological characterization of a PEAC from a patient with two other family members affected by similar lung tumors, which has never been reported before.

Findings

We evaluated the molecular characteristics of the proband’s PEAC by using a previously validated 47-microRNA (miRNA) cancer-specific array and a predictive method to estimate tissue-of-origin probabilities. Immunohistochemical (IHC) staining for thyroid transcription factor (TTF-1), napsin A, caudal-related homeobox 2 (CDX2), cytokeratins, and mucins, as well as mutational analyses for epidermal growth factor receptor (EGFR), Kirsten rat sarcoma viral oncogene homolog (KRAS), and anaplastic lymphoma kinase (ALK) were performed on formalin-fixed, paraffin-embedded (FFPE) tissues.

The occurrence of PEAC in two family members was associated with similar clinicopathological features (age at diagnosis, smoking habit, tumor localization, multiple colonic polyps), histologic findings (TTF-1 negativity and CDX2 positivity), and genetic findings (KRAS (Gly12Asp) mutation, but no EGFR/ALK aberrations). miRNA profiling revealed similarities with non-small cell lung cancer (NSCLC; 75.98 %) and some overlap with pancreatic ductal adenocarcinoma (PDAC; 23.34 %), but not with colorectal cancer (CRC; less than 0.5 %). Notably, these PEACs share key PDAC-associated miRNAs associated with tumor aggressiveness (miR-31*/-126*/-506/-508-3p/-514).

Conclusions

We describe for the first time PEAC in members from the same family, associated with similar clinical and genetic features. miRNA profiling of the PEAC resembled a NSCLC signature, with partial overlap to a PDAC pattern. This could explain its aggressive behavior and therefore help to guide future tailored-therapeutic approaches.

Findings

Introduction

Primary pulmonary enteric adenocarcinoma (PEAC) is defined as a pulmonary adenocarcinoma with a predominant component (>50 %) of intestinal differentiation and with tumor cells positive for at least one intestinal marker such as caudal-related homeobox 2 (CDX2), CK20, or MUC2 [1]. It was first described in 1991 [2], and a total of 31 cases have been reported in the literature [213]. In the International Multidisciplinary Classification of Lung Adenocarcinoma (2011), PEAC was classified as a rare variant of invasive adenocarcinoma [1]. The overlap between PEAC and metastases from other tumors (such as colorectal cancer (CRC)) represents an important diagnostic problem [213]. Several studies have demonstrated the ability of microRNA (miRNA) profiling for the identification of the tissue of origin for ambiguous cancers, assisting in the diagnosis of carcinomas of unknown origin [14, 15].

In the present study, we performed miRNA profiling in conjunction with immunohistochemical (IHC) and genetic analyses on a formalin-fixed, paraffin-embedded (FFPE) tissue from a patient affected by PEAC. Notably, this patient had a sister and a brother who were also affected by PEAC, and we describe the main clinicopathological findings in these family members.

Materials and methods

Patients

The proband was a 68-year-old male, ex-smoker, with no history of malignancy or relevant comorbidities. Computed tomography (CT) and positron emission tomography ([18F]FDG-PET) showed a 3.5-cm nodule in the right lower lobe with high glucose uptake. Lobectomy revealed a stage-IIA non-small cell lung cancer (NSCLC) with PEAC histology. Colonoscopy showed multiple benign polyps, but excluded primary CRC. Post-operative CT and [18F]FDG-PET revealed a small (8 mm) osteolytic area in the left hip bone. Considering the predominant component of intestinal differentiation, the patient was treated with XELOX (capecitabine + oxaliplatin) and a bisphosphonate, though with a very short CT scan re-evaluation, which after 2 cycles of chemotherapy showed progression of the bone lesion (2 cm). Consequently, the patient underwent surgical removal of this unique bone lesion, which was confirmed histologically as a PEAC metastasis. After 4 months, the patient relapsed with multiple bone metastases and underwent 4 cycles of carboplatin + pemetrexed which did not arrest tumor progression. Subsequent chemotherapy with docetaxel resulted in disease stability after 2 cycles.

Similarly, the 71-year-old proband’s sister underwent lung lobectomy revealing a stage-IB PEAC (Table 1), which progressed in 12 months with lung and adrenal gland metastases, but after 6 cycles of carboplatin + pemetrexed, followed by pemetrexed alone, became stable. Another brother (72 years old) presented with an advanced stage lung tumor (pleural effusion and multiple lytic bone metastases) and succumbed 1 month after diagnosis.

Table 1 Demographic and disease characteristics of two patients affected by PEAC within one family. The familial aggregation of PEAC was associated with similar clinicopathological features (age at diagnosis, smoking habit, tumor localization, multiple colon polyps), histologic findings (IHC staining negative for TTF-1 and positive for CDX2), and genetic findings (KRAS(Gly12Asp) mutation but no EGFR/ALK aberrations)

IHC and genetic and microarray analyses

The IHC analysis evaluated thyroid transcription factor (TTF-1), napsin A, cytokeratins (CK7/CK20), CDX-2, and mucins (MUC1/MUC2/MUC5AC/MUC6). For the genetic analyses, the neoplastic area was macro-dissected using a parallel hematoxylin and eosin (H&E) slide as a reference, with subsequent DNA/RNA isolation from samples with >80 % tumor infiltration. Epidermal growth factor receptor (EGFR) and Kirsten rat sarcoma viral oncogene homolog (KRAS) analyses were performed by PCR and direct sequencing, respectively. Anaplastic lymphoma kinase (ALK) status was studied by fluorescence in situ hybridization. Microsatellite instability (MSI) was evaluated by multiplex amplification and sequencing.

miRNA profiling was performed using the Agilent Human miRNA microarray v.2 (#G4470B; Agilent Technologies, Santa Clara, CA, USA). The Feature Extraction 10.7 and GeneSpring GX13 software were used to analyze the microarray raw, and we applied a previously validated method to estimate tissue-of-origin probabilities [15]. Manhattan correlation was used as a measure of similarity.

No samples for further pathological and genetic analyses were available from the proband’s sister, while cytological specimens of the proband’s brother were available to evaluate TTF-1, CDX-2, CK7, CK20, and for PCR analysis of selected miRNAs.

Further details on the analyses of EGFR, KRAS, ALK, MSI, and miRNA-PCR are in Additional file 1.

Results

Clinicopathological characteristics

The occurrence of PEAC in the proband and his sister was associated with similar clinicopathological features, summarized in Table 1, such as comparable age at diagnosis, smoking history, and identical localization of the tumor. Colonoscopy in both cases revealed multiple benign polyps and diverticulosis of the sigmoid colon. Notably, the proband’s brother had also comparable age and smoking history, but the dramatic history of his disease precluded additional analyses.

Histology, genetics, and miRNA profiling

Macroscopically, a well-circumscribed lung nodule with a maximum diameter of 3.5 cm was found in both cases. Microscopically, the nodule was composed of medium-to-large complex glands, and tumor cells were cuboidal to tall columnar with a brush border and eosinophilic cytoplasm. Areas of irregular necrosis were present, as previously described [2]. H&E revealed two different phenotypes in the lung lesion of the proband: NSCLC (30 %) and PEAC (70 %; Fig. 1). Expression of TTF-1 and CDX2 are associated with NSCLC and PEAC, respectively. Similar results were observed in the cytological specimens of the proband’s brother (Additional file 2: Figure S1). Moreover, we observed strong coexpression of MUC1 and MUC5AC in the areas with PEAC histology, while in the same areas, we did not detect MUC2. In particular, the panel 4 of the Fig. 1 summarizes the MUCs expression compared to the subtypes of pancreatic intraductal papillary mucinous neoplasms (IPMNs). Considering the intestinal-type pulmonary histology and multiple colonic polyps, we evaluated on the tissue of the proband the MSI markers, which suggested the absence of microsatellite instability (MSS phenotype). KRAS (Gly12Asp) mutations but no EGFR/ALK aberrations were found.

Fig. 1
figure 1

P1 Histological analyses (hematoxylin and eosin) revealed the presence of two different phenotypes in the lung lesion of the proband. We found non-small cell lung cancer (NSCLC, 30 %; A magnification ×10; B magnification ×40) and pulmonary enteric adenocarcinoma (PEAC, 70 %; C magnification ×10; D magnification ×40). Table 1 summarizes the marker expression associated with these two different phenotypes. The expression of TTF1 and CDX2 was associated with NSCLC and PEAC, respectively. Strong expression of CK7 was observed in both phenotypes (100 %). P2 Immunohistochemical signatures in NSCLC (left panels) and PEAC (right panels). TTF1 (A, B, C, and D); CK7 (E , F , G, and H); CDX2 (I , L , M, and N); and CK20 (O , P , Q, and R) staining showing co-localizations of TTF1/CK7 and CDX2/CK7 in NSCLC and PEAC, respectively. Columns 1 and 3 (magnification ×10); columns 2 and 4 (magnification ×40). P3 Mucins (MUCs) expression in PEAC and pancreatic intraductal papillary mucinous neoplasms (IPMNs). Strong coexpression of MUC1 and MUC5AC was observed in the area with PEAC histology. Conversely, in the same area, we did no detect MUC2 expression. Examples of the positive expression of MUC1 (D), MUC2 (E), and MUC5AC (F) associated with different subtypes of IPMNs. Columns 1 and 2 (magnification ×40). P4 Summary of mucins (MUCs) expression in the different subtypes of IPMNs and in our PEAC samples showing similar molecular markers expression. In particular, the negative expression for MUC6 (data not shown) suggests a pancreatobiliary differentiation of areas within our PEAC sample

The molecular characteristics of the proband’s PEAC were further evaluated using a miRNA profiling panel to predict tumor origin. Figure 2 shows the heat map representation of unsupervised clustering of 47 miRNAs in the proband’s sample (“sample BO01”) compared to primary and metastatic tumors of NSCLC, CRC and pancreatic ductal adenocarcinoma (PDAC), and primary gastric, kidney, prostate, skin, and liver cancers, as described previously [15]. This profiling revealed striking similarities with NSCLC (75.98 %) and a partial overlap with PDAC (23.34 %), but not with CRC (<0.5 %). The microarray data were confirmed by PCR analyses, focusing on five miRNAs, which showed a similar expression pattern in the proband’s brother (Additional file 3: Figure S2).

Fig. 2
figure 2

Unsupervised cluster analysis of 47 miRNAs in the proband’s sample of PEAC (“sample BO01”) compared to primary and metastatic tumors of lung cancer, CRC and PDAC, and also primary tumors of gastric, kidney, prostate, skin, and liver cancer. The colors represented on the heat-map correspond to the normalized miRNA expression value: green indicates downregulation while red indicates upregulation referred to the miRNA average value across all samples

Discussion

To our knowledge, this is the first study describing PEACs in members of the same family, associated with clinicopathological characteristics, histological, and genetic findings. Recently, El Hammoumi et al. described a patient with PEAC localized in the right lower lobe and associated with primary CRC [13]. However, miRNA profiling of the proband’s PEAC resembled that of lung cancer, with partial overlap of PDAC, but not CRC. This might at least in part explain the shorter disease-free survival (DFS) in the proband and his sister (1.5 and 12 months, respectively) compared to radically resected NSCLC patients, such as the 171 stage IA–IIIA NSCLC patients evaluated by Park et al., with a median DFS of 20 months [16].

Notably, we observed two different phenotypes (NSCLC and PEAC) in the lung lesions of the proband, which differed in IHC staining for TTF-1 and CDX2 markers. The NSCLC phenotype was characterized by TTF-1 positivity and CDX2 negativity. Conversely, the PEAC phenotype did not express the “lung cancer marker” TTF-1 but was instead characterized by the “intestinal marker” CDX2. Similar findings were observed in the proband’s sister and brother, underlying the potential role of these two markers (TTF-1, CDX2). Conversely, the proband and his sister showed different immunostaining for CK7 and CK20, commonly utilized to separate colorectal from other primary adenocarcinomas, especially the lung [12]. However, a previous description of 15 PEACs reported positivity for CK7 and CK20 in 93 and 27 % of cases, respectively [6].

Genetic analyses revealed KRAS (Gly12Asp) mutations, while no EGFR mutations or ALK rearrangements were described, in agreement with previous studies [213].

Using a 47-miRNA diagnostic array, we found that the miRNA profile of the proband’s PEAC was similar to that of lung carcinomas, but not CRC. However, we found a partial overlapping miRNA profile with PDAC, sharing some key miRNAs associated with tumor aggressiveness (miR-31*/-126*/-506/-508-3p/-514). In particular, miR-126*, mir-506, mir-508-3p, and mir-514 act as tumor suppressor genes, and their downregulation has been associated with proliferation, invasiveness, and poor survival outcomes [1722]. Furthermore, miR-31* acts as oncogene and its’ overexpression was an independent prognostic for shorter survival [23]. The MUCs expression in our PEACs compared to the different subtypes of IPMNs suggest pancreatobiliary differentiation and again, not intestinal, despite histological similarity [24].

Qureshi et al. [8] described a patient affected by metastatic PEAC, showing good response after 4 cycles with carboplatin and pemetrexed. The same treatment was associated with disease progression in our patient, while his sister obtained disease stability after 6 cycles, suggesting that further studies are needed to define if the combination of platinum-based drugs and pemetrexed can be used as an option for first-line treatment in PEAC. Conversely, we do not recommend the CRC regimens, since we observed disease progression after XELOX, and the miRNA profiling showed no concordance with CRC. Finally, we observed disease stabilization after docetaxel, which is interesting, as taxanes showed efficacy in both NSCLC and PDAC. On the basis of these observations, as well as on the shorter DFS and miRNA profiling, we presume that PEAC represents an aggressive histological subtype of lung adenocarcinoma. Therefore, building upon previous evidence on the key role of histology and genetics to improve the management of advanced NSCLC, especially lung adenocarcinoma, we hypothesize that a better molecular understanding of PEAC will boost future therapeutic options.

Ethics and consent

The research reported in the manuscript was performed with the approval of the appropriate ethics committee. Written informed consent was obtained for the molecular analyses and publication of this study.

Abbreviations

ALK:

anaplastic lymphoma kinase

CDX2:

caudal-related homeobox 2

CK:

cytokeratins

CRC:

colorectal cancer

CT:

computed tomography

DFS:

disease-free survival

EGFR:

epidermal growth factor receptor

FFPE:

formalin-fixed, paraffin-embedded

H&E:

hematoxylin & eosin

IHC:

immunohistochemistry

IPMN:

intraductal papillary mucinous neoplasm

KRAS:

Kirsten rat sarcoma viral oncogene homolog

miRNA:

microRNA

MSI:

microsatellite instability

MSS:

microsatellite stability

NEG:

negative

NSCLC:

non-small cell lung cancer

PDAC:

pancreatic ductal adenocarcinoma

PEAC:

primary pulmonary enteric adenocarcinoma

PET:

positron emission tomography

POS:

positive

TTF-1:

thyroid transcription factor

References

  1. Travis WD, Brambilla E, Noguchi M, Nicholson AG, Geisinger K, Yatabe Y, et al. American Thoracic Society. International Association for the Study of Lung Cancer/American Thoracic Society/European Respiratory Society international multidisciplinary classification of lung adenocarcinoma. J Thorac Oncol. 2011;6:244–85.

    Article  PubMed Central  PubMed  Google Scholar 

  2. Tsao MS, Fraser RS. Primary pulmonary adenocarcinoma with enteric differentiation. Cancer. 1991;68(8):1754–7.

    Article  CAS  PubMed  Google Scholar 

  3. Inamura K, Satoh Y, Okumura S, Nakagawa K, Tsuchiya E, Fukayama M, et al. Pulmonary adenocarcinomas with entric differentiation: histologic and immunohistochemical characteristics compared with metastatic colorectal cancers and usual pulmonary adenocarcinomas. Am J Surg Pathol. 2005;29(5):660–5.

    Article  PubMed  Google Scholar 

  4. Maeda R, Isowa N, Onuma H, Miura H. Pulmonary intestinal-type adenocarcinoma. Interact Cardiovasc Thorac Surg. 2008;7(2):349–51.

    Article  PubMed  Google Scholar 

  5. Li HC, Schmidt L, Greenson JK, Chang AC, Myers JL. Primary pulmonary adenocarcinoma with intestinal differentiation mimicking metastatic colorectal carcinoma: case report and review of literature. Am J Clin Pathol. 2009;13(1):129–33.

    Article  CAS  Google Scholar 

  6. Hatanaka K, Tsuta K, Watanabe K, Sugino K, Uekusa T. Primary pulmonary adenocarcinoma with enteric differentiation resembling metastatic colorectal carcinoma: a report of the second case negative for cytokeratin 7. Pathol Res Pract. 2011;207(3):188–91.

    Article  CAS  PubMed  Google Scholar 

  7. Lin D, Zhao Y, Li H, Xing X. Pulmonary enteric adenocarcinoma with villin brush border immunoreactivity: a case report and literature review. J Thorac Dis. 2013;5(1):E17–20.

    Article  PubMed Central  PubMed  Google Scholar 

  8. Qureshi A, Furrukh M. Enteric adenocarcinoma lung: a rare presentation in an Omani woman. BMJ Case Rep 2013

  9. László T, Lacza A, Tóth D, Molnár TF, Kálmán E. Pulmonary enteric adenocarcinoma indistinguishable from metastatic colorectal carcinoma morphologically and immunohistologically. Histopathology. 2014;65(2):283–7.

    Article  PubMed  Google Scholar 

  10. Wang CX, Liu B, Wang YF, Zhang RS, Yu B, Lu ZF, et al. Pulmonary enteric adenocarcinoma: a study of the clinicopathologic and molecular status of nine cases. Int J Clin Exp Pathol. 2014;7(3):1266–74.

    PubMed Central  PubMed  Google Scholar 

  11. Stojsic J, Kontic M, Subotic D, Popovic M, Tomasevic D, Lukic J. Intestinal type of lung adenocarcinoma in younger adults. Case Rep Pulmonol. 2014;2014:282196.

    PubMed Central  PubMed  Google Scholar 

  12. Yousem SA. Pulmonary intestinal-type adenocarcinoma does not show enteric differentiation by immunohistochemical study. Mod Pathol. 2005;18(6):816–21.

    Article  PubMed  Google Scholar 

  13. El Hammoumi MM, El Ochi R, Kabiri EH. Primary lung adenocarcinoma with enteric morphology associated with primary colon adenocarcinoma. Arch Bronconeumol, 2015 [Epub ahead of print].

  14. Rosenfeld N, Aharonov R, Meiri E, Rosenwald S, Spector Y, Zepeniuk M, et al. MicroRNAs accurately identify cancer tissue origin. Nat Biotechnol. 2008;26(4):462–9.

    Article  CAS  PubMed  Google Scholar 

  15. Ferracin M, Pedriali M, Veronese A, Zagatti B, Gafà R, Magri E, et al. MicroRNA profiling for the identification of cancers with unknown primary tissue-of-origin. J Pathol. 2011;225(1):43–53.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Park C, Lee IJ, Jang SH, Lee JW. Factors affecting tumor recurrence after curative surgery for NSCLC: impacts of lymphovascular invasion on early tumor recurrence. J Thorac Dis. 2014;6(10):1420–8.

    PubMed Central  PubMed  Google Scholar 

  17. Permuth-Wey J, Chen YA, Fisher K, McCarthy S, Qu X, Lloyd MC, et al. A genome-wide investigation of microRNA expression identifies biologically-meaningful microRNAs that distinguish between high-risk and low-risk intraductal papillary mucinous neoplasms of the pancreas. PLoS One. 2015;10(1):e0116869.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Frampton AE, Krell J, Jacob J, Stebbing J, Castellano L, Jiao LR. Loss of miR-126 is crucial to pancreatic cancer progression. Expert Rev Anticancer Ther. 2012;12(7):881–4.

    Article  CAS  PubMed  Google Scholar 

  19. Yin M, Ren X, Zhang X, Luo Y, Wang G, Huang K, et al. Selective killing of lung cancer cells by miRNA-506 molecule through inhibiting NF-κB p65 to evoke reactive oxygen species generation and p53 activation. Oncogene. 2014. doi:10.1038/onc.2013.597.

    Google Scholar 

  20. Sun Y, Hu L, Zheng H, Bagnoli M, Guo Y, Rupaimoole R, et al. MiR-506 inhibits multiple targets in the epithelial-to-mesenchymal transition network and is associated with good prognosis in epithelial ovarian cancer. J Pathol. 2015;235(1):25–36.

    Article  CAS  PubMed  Google Scholar 

  21. Zhai Q, Zhou L, Zhao C, Wan J, Yu Z, Guo X, et al. Identification of miR-508-3p and miR-509-3p that are associated with cell invasion and migration and involved in the apoptosis of renal cell carcinoma. Biochem Biophys Res Commun. 2012;419(4):621–6.

    Article  CAS  PubMed  Google Scholar 

  22. Garajová I, Le Large TY, Frampton AE, Rolfo C, Voortman J, Giovannetti E. Molecular mechanisms underlying the role of microRNAs in the chemoresistance of pancreatic cancer. Biomed Res Int. 2014;2014:678401.

    Article  PubMed Central  PubMed  Google Scholar 

  23. Frampton AE, Giovannetti E, Jamieson NB, Krell J, Gall TM, Stebbing J, et al. A microRNA meta-signature for pancreatic ductal adenocarcinoma. Expert Rev Mol Diagn. 2014;14(3):267–71.

    Article  CAS  PubMed  Google Scholar 

  24. Castellano-Megias VM, Ibarrola-de Andres C, Lopez-Alfonso G, Colina-Ruizdelgado F. Pathological features and diagnosis of intraductal papillary mucinous neoplasm of the pancreas. World J Gastrointest Oncol. 2014;6(9):311–24.

    Article  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

This work was partially supported by grants from Istituto Toscano Tumori (ITT) and AIRC Start-Up grant (to Niccola Funel and Elisa Giovannetti).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Elisa Giovannetti.

Additional information

Competing interests

The authors declare no conflicts of interest.

Authors’ contributions

All of the authors contributed to drafting or revising the article for important intellectual content. All authors gave final approval of the manuscript to be published.

Guido Biasco and Elisa Giovannetti contributed equally to this work.

Additional files

Additional file 1:

Analyses of EGFR, KRAS, ALK and MSI. Detailed methods used for the analyses of EGFR, KRAS mutation, ALK rearrangement, MSI, and miRNA-PCR. (DOC 108 kb)

Additional file 2: Supplemental Figure S1.

Immunohistochemical analyses in the cytological sample of the proband’s brother. The stainings for TTF1, CK7, CDX-2, and CK20, summarized in the table, showed patterns of protein expression similar to the expression levels observed in the PEAC component of the proband’s sample. Left column (original magnification ×20); right column (original magnification ×40). (PPTX 32,747 kb)

Additional file 3: Supplemental Figure S2.

PCR analysis of miR-31*/-126*/-506/-508-3p/-514. Quantitative “miRNA-targeted” real-time PCR of five miRNA aberrantly expressed in the proband’s PEAC according to the microarray results showed that the microarray data were validated by PCR data as well as that the pattern of modulation of expression of these miRNA was similar in the proband’s brother. Columns and bars, average levels ± SD, compared using the 2-ΔΔCt method to a pool of normal lung tissues, dashed line. (PPTX 158 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Garajová, I., Funel, N., Fiorentino, M. et al. MicroRNA profiling of primary pulmonary enteric adenocarcinoma in members from the same family reveals some similarities to pancreatic adenocarcinoma—a step towards personalized therapy. Clin Epigenet 7, 129 (2015). https://doi.org/10.1186/s13148-015-0162-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13148-015-0162-5

Keywords