Talita Aguiar1,2,3†
Anne Teixeira1,2†Marília O. Scliar2
Juliana Sobral de Barros1,2Renan B. Lemes1,2Silvia Souza1,2
Giovanna Tolezano1,2Fernanda Santos4
Israel Tojal5Monica Cypriano6
Silvia Regina Caminada de Toledo6
Eugênia Valadares7
Raquel Borges Pinto8
Osvaldo Afonso Pinto Artigalas9Joaquim Caetano de Aguirre Neto10Estela Novak11,12Lilian Maria Cristofani11
Sofia M. Miura Sugayama13Vicente Odone11
Isabela Werneck Cunha14
Cecilia Maria Lima da Costa4Carla Rosenberg1,2
Ana Krepischi1,2*- 1Department of Genetics and Evolutionary Biology, Institute of Biosciences, University of São Paulo, São Paulo, Brazil
- 2Human Genome and Stem Cell Research Center, Institute of Biosciences, University of São Paulo, São Paulo, Brazil
- 3Columbia University Irving Medical Center, New York, NY, United States
- 4Department of Pediatric Oncology, A. C. Camargo Cancer Center, São Paulo, Brazil
- 5International Center for Research, A. C. Camargo Cancer Center, São Paulo, Brazil
- 6GRAACC—Grupo de Apoio Ao Adolescente e Criança Com Câncer, Federal University of São Paulo, São Paulo, Brazil
- 7Benjamim Guimarães Foundation - Department of Pediatrics Hospital da Baleia, Belo Horizonte, Brazil
- 8Department of Genetics, Hospital da Criança Conceição, Hospitalar Conceição Group, Porto Alegre, Brazil
- 9Department Pediatric Gastroenterology, Hospital da Criança Conceição, Hospitalar Conceição Group, Porto Alegre, Brazil
- 10Paediatric Haemato-oncology, Hospital Santa Casa de Belo Horizonte, Belo Horizonte, Brazil
- 11Pediatric Cancer Institute (ITACI) at the Pediatric Department, São Paulo University Medical School, São Paulo, Brazil
- 12Molecular Genetics—Foundation Pro Sangue Blood Center of São Paulo, São Paulo, Brazil
- 13Department of Pediatric, Faculty of Medicine of the University of São Paulo, São Paulo, Brazil
- 14Department of Pathology, Rede D’OR-São Luiz, São Paulo, Brazil
The ultrarare hepatoblastoma (HB) is the most common pediatric liver cancer. HB risk is related to a few rare syndromes, and the molecular bases remain elusive for most cases. We investigated the burden of rare damaging germline variants in 30 Brazilian patients with HB and the presence of additional clinical signs. A high frequency of prematurity (20%) and birth defects (37%), especially craniofacial (17%, including craniosynostosis) and kidney (7%) anomalies, was observed. Putative pathogenic or likely pathogenic monoallelic germline variants mapped to 10 cancer predisposition genes (CPGs: APC, CHEK2, DROSHA, ERCC5, FAH, MSH2, MUTYH, RPS19, TGFBR2 and VHL) were detected in 33% of the patients, only 40% of them with a family history of cancer. These findings showed a predominance of CPGs with a known link to gastrointestinal/colorectal and renal cancer risk. A remarkable feature was an enrichment of rare damaging variants affecting different classes of DNA repair genes, particularly those known as Fanconi anemia genes. Moreover, several potentially deleterious variants mapped to genes impacting liver functions were disclosed. To our knowledge, this is the largest assessment of rare germline variants in HB patients to date, contributing to elucidate the genetic architecture of HB risk.
Introduction
The etiology of pediatric cancer is largely unknown (Saletta et al., 2015). Despite intensive research, gaps remain in our understanding of the genetic landscape of pediatric cancer susceptibility. It is assumed that germline mutations in cancer-predisposing genes (CPGs) in children and adolescents are rare events (Downing et al., 2012; Zhang et al., 2015; Kindler et al., 2018). However, the genetic predisposition to childhood cancer is most likely under-identified. Most investigations have focused on known CPGs and sequenced patients without parental samples, an approach that impairs a broad evaluation of the full range of genetic mechanisms underlying pediatric cancer risk, such as de novo mutations and the identification of new candidate CPGs. Recent studies of large cohorts of pediatric cancer patients have confirmed that approximately 8–18% of patients carry a germline pathogenic variant in a broad spectrum of known CPGs (Zhang et al., 2015; Ma et al., 2018a; Grobner et al., 2018; Alejandro Sweet-Cordero and Biegel, 2019; Akhavanfard et al., 2020; Capasso et al., 2020; Newman et al., 2021). These studies also highlighted that isolated factors, such as tumor type and a positive family history of cancer, have low predictive power for the presence of germline CPG mutations.
HB is the most common malignant liver tumor in the pediatric population (Heck et al., 2013), although it is considered an ultrarare disease, accounting for only 1% of all pediatric tumors (Stiller et al., 2006; Czauderna and Garnier, 2018; Feng et al., 2019). In Brazil, collected data on HB are concordant with the worldwide incidence of 0.5–1.5 cases per million (Institute Nacional de Câncer José Alencar Gomes da Silva (INCA), 2016; Aguiar et al., 2020). Most cases are diagnosed before the age of 4 years, and a male preponderance is reported (Ries et al., 1999). Nongenetic factors known to be associated with HB risk are related to very low birth weight (<1,500 g), including preterm birth (<33 weeks), small for gestational age and multiple birth pregnancies (Heck et al., 2013; Turcotte et al., 2014), and in vitro fertilization (Spector et al., 2019). A slow increase in HB incidence is observed in North America and Europe (Pateva et al., 2017), which can be partly due to the increased survival of children with low birth weight (Tanimura et al., 1998). An increased risk for HB development has been reported in association with a few specific genetic conditions, including Beckwith-Wiedemann syndrome (DeBaun and Tucker, 1998; Kim et al., 2017), familial adenomatous polyposis (APC gene) (Hirschman et al., 2005), Li-Fraumeni syndrome (TP53 gene) (Curia et al., 2008), Aicardi syndrome (Kamien and Gabbett, 2009), and trisomy 18 (Pereira et al., 2012).
Here, we investigated the germline exome of 30 children who developed HB, 13 of whom (43%) exhibited additional clinical signs. Our analysis provides a framework for investigating candidate genes involved in HB predisposition, as well as the tumor association with specific birth defects.
Patients and Methodology
Participants
Thirty children diagnosed with HB were enrolled in this study, which was performed at the Institute of Biosciences, University of São Paulo, Brazil. Patients were recruited from five different Brazilian institutions, most of them (n = 27) from three large pediatric cancer centers of the city of São Paulo, namely, the A. C. Camargo Cancer Center, Adolescent and Child with Cancer Support Group (GRAACC), and Pediatric Cancer Institute (ITACI). In addition, three patients were recruited from other institutions, including the Hospital da Baleia (n = 1), Hospital da Criança Conceição (n = 1), and Hospital São Lucas (n = 1). This study was approved by the Research Ethics Committee of the Institute of Biosciences (CAAE: 09163818.4.0000.5464), and informed consent was obtained from the parents. Patient’s clinical data were documented or recovered from medical records by the oncologists who contributed to this work. Peripheral blood samples were collected from 28 patients, and normal liver tissues were recovered from two patients. Genomic samples were also obtained from 27 available parents (19 patients).
Library Preparation and Whole-Exome Sequencing (WES)
DNA samples were extracted by phenol-chloroform followed by ethanol precipitation (Sambrook and Russell, 2006). Genomic libraries were constructed with 1 µg of genomic DNA using the following kits: Sureselect QXT V6 (Agilent Technologies), OneSeq Constitutional Research Panel (Agilent Technologies), or xGen Exome Research Panel v1.0 (IDT - Integrated DNA Technologies). The sequencing of enriched libraries was performed on the Illumina HiSeq 2,500 platform using 150 base paired end reads. The sequences were aligned to the GRCh37/hg19 human genome reference with the BWA_MEM algorithm (Li, 2013). Picard tools (v.1.8, http://broadinstitute.github.io/picard/) were used to convert the SAM file into BAM and to mark PCR duplicates. The Genome Analysis Toolkit (GATK 3.7) (McKenna et al., 2010) was used to realign indels, recalibrate the bases, and call (Unified Genotyper) and recalibrate variants (VQSR). Finally, multiallelic variants were split into different lines using the script split_multiallelic_rows.rb from Atlas2 (Challis et al., 2012) to obtain the VCF files used for analysis.
WES Data Analysis
SNV and indel variant annotation were conducted using VarSeq software version 1.5.0 (Golden Helix) and selected by the reading depth (>10), Phred score (>20), and alternative allele frequency (>0.35). Based on the public variant databases in ABraOM (http://abraom.ib.usp.br - (Naslavsky et al., 2020), GnomAD (https://gnomad.broadinstitute.org-(Karczewski et al., 2020), and 1,000 Genomes - phase three {https://www.internationalgenome.org–(Clarke et al., 2017)}, we filtered out germline variants with frequencies below 1%, as well as those mapped to hypervariable genes (Fuentes Fajardo et al., 2012) or detected in an in-house dataset including data from 19 healthy controls.
Coding variants - Coding nonsynonymous missense and loss-of-function (LoF; frameshift, stop loss/gain, essential splice site, nonsense) variants were maintained for further analysis. In silico pathogenicity prediction for missense variants was based on six algorithms provided by the database dbNSFP (version 2.4); those predicted to be damaging to protein function by at least five different tools were prioritized. The final set of genes with rare damaging coding variants was annotated using Varelect (Stelzer et al., 2016) and HPO (Köhler et al., 2019) for ranking in association with specific phenotypes. All LoF and prioritized missense variants were validated by visual inspection of the BAM files, further annotated using the Varsome tool (Kopanos et al., 2019), and classified according to the American College of Medical Genetics and Genomics (ACMG) guidelines (David Bick et al., 2015; Hampel et al., 2015).
Noncoding variants - Intronic, intergenic, 3′ and 5’ untranslated regions (UTRs), and splice region variants were annotated using SNPnexus v4 (Oscanoa et al., 2020), a web-based annotation tool for the analysis and interpretation of variants that includes databases of regulatory elements and regions, such as miRbase (ftp://mirbase.org/pub/mirbase/20/genomes/), Vista HMR (http://hgdownload.cse.ucsc.edu/goldenPath/hg19/database/), and ENCODE (ftp://ftp.ensembl.org/pub/grch37/release-95/mysql/regulation_mart_95/); phenotype and disease association (Genetic Association of Complex Diseases and Disorders (GAD) - http://hgdownload.cse.ucsc.edu/goldenPath/hg19/database/, ClinVar, and COSMIC); and noncoding scoring (Combined Annotation Dependent Depletion (CADD) - http://cadd.gs.washington.edu/ download, Fitness Consequences of Functional Annotation (fitCons)-http://compgen.cshl.edu/fitCons/0downloads/tracks/V1.01/i6/scores/; and Chromatin Effects of Sequence Alterations (DeepSEA)-http://deepsea.princeton.edu/help/).
Copy number variants (CNVs) - CNVs and region of homozygosity (ROH) events were derived from WES data using the software Nexus Copy Number 9 (Biodiscovery) with the SNP-FASST2 segmentation algorithm (threshold log2 Cy3/Cy5 ratio of |0.3| for gains and losses; minimum ROH length of 5 Mb). Common CNVs (Database of Genomic Variants, http://dgv.tcag.ca/dgv/app/home) were disregarded. For CNV validation, chromosome microarray analysis (CMA) was performed using the 180K platform (Agilent Technologies), as previously reported (Oliveira et al., 2018).
Results
Clinical Characterization of the HB Cohort
The clinical features of the 30 patients who developed HB are described in Table 1; the details of their tumors can be found in Supplementary Table S1. The mean age at HB diagnosis was 24 months, excluding one patient who was diagnosed at 17 years (P07). Twenty-one patients were male (70%), which agrees with the literature on sex bias in HB (Spector and Birch, 2012; Feng et al., 2019). Six patients (∼20%) presented with a family history of cancer (relatives of different degrees developed different tumors at varying ages).
TABLE 1. Clinical description of the 30 cases of patients who developed hepatoblastoma.
Sixteen patients were diagnosed with high-risk HB, classified according to the CHIC stratification (Czauderna et al., 2016; Meyers et al., 2017), and 10 presented pulmonary metastasis at diagnosis. Except for P23, who underwent surgery at diagnosis, all patients received neoadjuvant chemotherapy protocols followed by tumor resection or transplantation. Eight patients received a liver transplant, and two relapsed. Four patients died from the disease.
Six out of the 30 patients were born prematurely (<37 weeks), corresponding to 20% of the group. Eleven patients had birth defects (37%), and eight of them were classified as syndromic due to the presentation of two or more congenital clinical features. Among them, five patients had craniofacial anomalies, two of whom were diagnosed with craniosynostosis (P13 and P28); two patients were born with kidney anomalies (P11 and P09); and P07, who was diagnosed with HB at an advanced age (17 years), was born with mild hepatomegaly.
P11, female, had a congenital HB diagnosed at 1 month of age; in addition, this patient was born with unilateral renal agenesis. P13, a male, was born extremely premature at 27 weeks. More details about the clinical features of both patients can be found in our previous study (Aguiar et al., 2020).
Two patients were born with Hirschsprung disease and other clinical features (P17 and P18). P17, female, was the third child of no consanguineous parents. She was born at term, and her two siblings had a normal phenotype. Abdominal volume and no evacuation were detected in the first 24 h after birth, and the diagnosis of Hirschsprung disease was made. In the clinical evaluation at 5 years old, she presented with global neuropsychomotor delay, facial dysmorphisms, clinodactyly, and nail dysplasia (hypoplastic). HB was diagnosed at the age of four and classified as an epithelial subtype with a predominance of embryonal cells, PRETEXT IV, high risk. She underwent the SIOPEL six chemotherapy protocol and died before the surgical procedure. P18, male, was the third child of a consanguineous couple; his sister was born with congenital bilateral cataracts, while his brother exhibited intestinal atresia-terminal ileus. The patient was born prematurely (28 weeks), with a syndromic phenotype composed of congenital ileal atresia, bilateral cataracts, and sensorineural deafness. His mother, who had gestational risk (cardiac defect and preeclampsia), died during his birth due to congestive heart failure. HB tumors were diagnosed at 1 year of age and classified as fetal epithelial subtype, PRETEXT II, and low risk. He underwent a chemotherapy protocol with cisplatin, doxorubicin, and ifosfamide, followed by partial hepatectomy. Currently, the patient is in post treatment follow-up, and clinical details have been previously published (Pinto et al., 2016).
P24 had facial dysmorphisms and dysplastic nails of the hands and feet, in addition to developmental delay. P29 was born with congenital malformations of the VACTERL spectrum and exhibited postnatal microcephaly and developmental delay. P30 presented craniofacial dysmorphisms, turricephaly, short neck, laryngomalacia, swallowing disorder, severe bronchodysplasia, and polysyndactyly of the right fifth digit, in addition to severe malnutrition and neuropsychomotor delay.
We performed comparisons between disease risk status and other variables (gender, premature birth, death, associated clinical conditions, familial clinical history, and age at diagnosis). When only categorical variables were considered, the non-parametric chi-square test with Yates continuity correction and Fisher exact test were applied to determine if the distribution of the data differed statistically from random expectations (null hypothesis). In cases in which continuous variables were considered, we performed the non-parametric Kruskal–Wallis test. Although no significant results were found, mainly due to the small size of the sample, when we compared the distribution of individuals with craniofacial or kidney defects according to risk stratification, a marginally significant result was obtained (p-value = 0.10), which can indicate a tendency of individuals with high-risk status to be more likely to present craniofacial defects, while individuals with less severe status are more prone to have kidney defects (Supplementary Figure S1).
Germline Coding and Noncoding Variants (SNVs and Indels)
Figure 1 summarizes the analysis workflow of the WES data. The mean sequencing depth of the exomes was 92× (Supplementary Table S2 presents sequencing metrics and the type of genomic library for each sample); P03 was the only sample presenting 10x on-target coverage below 80%.
FIGURE 1. Workflow of the exome sequencing analysis. Variants were filtered according to quality (Phred score >20, read depth >10, variant allele frequency >35%); coding and noncoding variants were separately analyzed; frequency (GnomAD, ABraOM, 1K Genomes); the effect on coding variants (frameshift, stop loss/gain, missense, splice site, nonsense); for missense variants, the prediction of pathogenicity in at least five out of six algorithms; HPO annotation (hepatoblastoma, abnormalities of the liver, and cancer). The filtered variants were visually examined using Integrative Genomics Viewer (IGV) software (http//www.broadinstitute.org/igv) to further filter out possible strand bias and homopolymeric region artifacts. All the filtered variants mapped to cancer predisposition genes were classified using the ACMG guidelines. One- Intronic variants, 3′UTR, 5′UTR; two- MAF: GnomAD, ABraOM, 1K Genomes; three- Frameshift, stop loss, stop gain, missense, splice site, nonsense variants; four- Missense variants with dbNSFP Functional Prediction of pathogenicity in at least five out of six algorithms; five- Terms used for Varelect and HPO annotation: Hepatoblastoma, abnormalities of the liver, and cancer. CNV—copy number variation, ROH - region of homozygosity, MAF—maximum allele frequency, CPG—cancer predisposition gene, VUS—variant of uncertain significance, HPO—Human Phenotype Ontology, HB—hepatoblastoma.
A total of 9,467 rare (population frequency <1%) germline coding nonsynonymous variants were detected in the cohort of 30 HB patients, mapping to 6,102 genes. Details of all variants can be found in Supplementary Table S3. A total of 2,107 of these rare variants, related to 1,737 different genes and including 1,671 missense mutations and 436 LoF variants (Figure 2A), met our criteria of a read depth >10, Phred score >20, and alternative allele frequency >0.35. Pathogenic (P) or likely pathogenic (LP) variants mapped to morbid OMIM genes that could explain the syndromic phenotypes of some patients were not detected.
FIGURE 2. Distribution of the detected high-quality rare coding and noncoding variants detected in 30 HB patients (A). Sequence ontology of the rare coding variants detected after selection by the read depth (>10), Phred score (>20), alternative allele frequency (>0.35), and population frequency (<1%). A total of 2,107 variants were classified into 1,671 missense mutations and 436 LoF variants (B). Sequence ontology of the rare noncoding variants detected after selection by the read depth (>10), Phred score (>20), alternative allele frequency (>0.35), and population frequency (<0.1%). A total of 2,070 noncoding variants were distributed in intronic, intergenic, and 3′ and 5′ UTRs.
Using a list of 222 CPGs composed of 119 known CPGs (49 of them reported in OMIM), in addition to 103 candidates that were compiled by revision of recent publications (Supplementary Table S4), we investigated the presence of rare coding variants that could be related to cancer development. No homozygous or compound heterozygous pathogenic (P) or likely pathogenic (LP) variants were observed in known or candidate CPGs. Eleven heterozygous putative P/LP variants mapped to ten CPGs were detected in ten patients, comprising 33% of the group (Table 2; details provided in Supplementary Table S5); VHL variants were detected in two patients. Among the 10 patients with potentially P/LP variants in CPGs, only four presented a family history of cancer; four of them were syndromic, and three were born prematurely. One of these patients (P28) carried two variants mapped to known recessive CPGs. Eight of these variants were detected in seven autosomal-dominant CPGs (known or candidate), including an APC LoF variant (OMIM #175100 FAMILIAL ADENOMATOUS POLYPOSIS one; gastrointestinal carcinomas) and six missense variants, mapped to the CHEK2 (LI-FRAUMENI SYNDROME two; colorectal, breast and prostate cancer), DROSHA (50), MSH2 (LYNCH SYNDROME I; colorectal cancer/MISMATCH REPAIR CANCER SYNDROME two; hematologic malignancy, brain tumors, and gastrointestinal tumors), RPS19 (DIAMOND-BLACKFAN ANEMIA one; osteogenic sarcoma, myelodysplastic syndrome, colon cancer), VHL (VON HIPPEL-LINDAU SYNDROME; renal cell carcinoma, pheochromocytoma, hemangioblastoma, hypernephroma, pancreatic cancer, paraganglioma, adenocarcinoma of the ampulla of Vater), and TGFBR2 (COLORECTAL CANCER, HEREDITARY NONPOLYPOSIS, TYPE 6) genes. Three variants were detected in three CPGs associated with recessive conditions: an ERCC5 LoF (XERODERMA PIGMENTOSUM, COMPLEMENTATION GROUP G; skin cancers) and missense variants mapped to FAH (TYROSINEMIA, TYPE I; hepatocellular carcinoma) and MUTYH (FAMILIAL ADENOMATOUS POLYPOSIS two; colorectal carcinomas).
TABLE 2. Description of 11 potentially pathogenic/likely pathogenic germline heterozygous variants mapped to 10 known/candidate cancer predisposition genes and detected in 10 hepatoblastoma patients.
We investigated genomic data available of tumor tissues derived from three of these patients with relevant germline variants described in Table 2 (P03, P13, P21); however, additional somatic mutations in the same gene were not observed in these cases, and tumor tissue from patient P13 was no longer available to test microsatellite instability. In particular, Patient P27, who carries an APC pathogenic variant, had a strong familial cancer history, in which her mother and sister were diagnosed with familial adenomatous polyposis (FAP), and her maternal grandmother and uncle, already dead, had colon cancer. Sanger sequencing analysis confirmed the segregation of the detected pathogenic APC variant with the cancer phenotype in this family (Figure 3).
FIGURE 3. APC pathogenic variant segregates with the cancer phenotype in the family of patient P27. The patient P27 is indicated with the black arrow. P27, her mother, sister, and niece are carriers of the pathogenic p. Cys1249* variant in the APC gene, which was identified by exome sequencing in the patient, and validated by Sanger sequencing in their indicated relatives. Her mother and sister had colon polyps and stomach polyps, respectively. The maternal grandmother and uncle had colon cancer (both are dead).
In addition, 44 VUS mapped to 34 known/candidate CPGs were observed in 21 patients (70%) (Table 3), most of whom carried more than one variant. VUS mapped to the ATM, BRCA2, COL7A1, DHCR7, DOCK8, FANCs, and GLI3 genes were detected in more than one patient.
TABLE 3. Description of 44 germline rare heterozygous variants of uncertain significance (VUS) mapped to known/candidate cancer predisposition genes detected in 21 hepatoblastoma patients.
Thirty-two genes related to liver differentiation or function were found to be affected by rare variants (ABCB11, ABCB4, ABCC2, ABCC3, AFP, ALB, CYP1A1, CYP1A2, CYP2C19, CYP2C8, CYP2C9, CYP2D6, CYP3A7, FAH, FOXA2, KRT7, KRT8, MET, NR1I2, ONECUT1, PAH, POU5F1, PPARG, SOX17) (Supplementary Table S6).
We also investigated whether the observed sex bias in the group could be explained by an increased burden of rare damaging variants in one of the sexes; we did not detect significant differences considering the average of rare damaging variants in male and female patients (∼66.7 and 68.1, respectively), LoF variants (∼15 variants in both groups), and rare damaging CPG variants (∼8 and 6, respectively).
A total of 2,069 noncoding variants passed our filters (read depth >10, Phred score >20, alternative allele frequency >0.35, frequency in population databases >0.1%), including intronic (76%), intergenic (7%), 3′ prime UTR (7%), 5′ prime UTR (6%), and splice region (4%) variants (Supplementary Table S7; Figure 2B). These variants were annotated using SNP Nexus (Oscanoa et al., 2020), and those with CADD scores above 15 and associated with cancer (https://geneticassociationdb.nih.gov/) were prioritized for further analyses. Table 4 details the 23 prioritized rare noncoding variants. Two variants were observed in the intronic regions of the CPGs BRAF and CREBBP, but with no evidence of a functional effect. In particular, P12 carried a de novo mutation in the 5’ UTR of the TCF7 gene, an important effector protein in the Wnt pathway; this patient also carries a paternally inherited coding TCF7 VUS (c.1060C > G).
TABLE 4. Rare germline non-coding variants detected in HB patients which were prioritized (CADD score above 15 and association with cancer).
Using the exome data of HB patients and a healthy control group (n = 19, data not shown), we inspected a list of 220 DNA repair genes (distributed in 16 categories, including several bona fide CPGs; https://www.mdanderson.org/documents/Labs/Wood-Laboratory/human-dna-repair-genes.html), searching for rare LoF and missense variants with high in silico damage prediction (5 or more algorithms; (Supplementary Table S8); Figure 4 shows the frequency of rare damaging variants detected in 12 DNA repair categories in both controls and patients. Thirty-four heterozygous variants mapped to DNA repair genes were observed in 21 patients (70%), while nine heterozygous variants were detected in nine healthy controls (47%). Although not statistically significant, there was an apparent excess of damaging variants mapped to DNA repair genes in patients. Moreover, rare damaging variants affecting specific DNA repair gene categories, such as ubiquitin modification, poly (ADP-ribose) polymerase (PARP) enzymes that bind to DNA, nonhomologous end-joining, homologous recombination (BRCA1, EME2, SPIDR, and RAD54L), Fanconi anemia genes (BRCA2, FANCA, BRIP1, SLX4, FANCD2, and FAAP24), genes associated with DNA sensitivity to damaging agents, and base excision repair, were observed only in patients. Significant enrichment for rare damaging variants mapped to Fanconi anemia genes was detected in the group of patients (p value 0.0338; Fisher’s test).
FIGURE 4. Frequency of high-quality rare germline coding variants mapped to DNA repair genes in HB patients and a control group. A list of 220 DNA repair genes distributed in 16 categories was analyzed; the 12 categories with variants detected in either patients or controls are represented. PARP—poly (ADP-ribose) polymerase. *p value 0.0338; Fisher’s test.
Germline Copy Number Variation (CNV) and Regions of Homozygosity
Assessing germline CNVs using exome data, seven rare copy number changes were detected, four with clinical relevance, all of which were validated by CMA (Supplementary Figure S2). The encompassed genomic regions did not include CPGs or genes associated with the patients’ phenotypes (Table 5). A pathogenic Y chromosome aneuploidy was detected in P19 (Jacob syndrome 47, XYY). P21 and P28 paternally inherited 15q11.2 copy number gain and loss, respectively, which are classified as risk factors for neurodevelopmental disorders, mainly global delay, and intellectual disability.
TABLE 5. Rare germline CNVs containing coding regions and copy neutral events (>10 Mb) detected by whole-exome sequencing.
Regions of homozygosity (ROH) were also investigated in the germline exome data (Table 5). The individuals P05, P16, and P18 harbor several large ROH, evidencing parental consanguinity. Germline 11p15.5 loss of heterozygosity (LOH) was not observed in any of the patients by exome analysis. However, in one of the patients (P21), a somatic 11p LOH event was detected in the respective tumor sample (Supplementary Figure S3).
Discussion
We explored the clinical features and spectrum of rare germline variants possibly associated with cancer risk in 30 children who developed HB. We report here an enrichment of birth defects in HB patients (37%), mainly craniofacial (17%, including craniosynostosis) and kidney (7%) anomalies, as well as Hirschsprung disease (7%) and nail dysplasia (7%). There is a growing robust bulk of evidence emphasizing that various birth defects occur in association with a significant increase in the risk of developing childhood cancer (Krivit and Good, 1957; Agha et al., 2005; Johnson et al., 2017). Recent studies also supported an increased risk of pediatric cancer in individuals with birth defects unrelated to chromosomal abnormalities or known genetic syndromes (Norwood et al., 2017; Siegel et al., 2017). HB are reported to occur in association with a wide variety of congenital abnormalities (Narod et al., 1997; Ansell et al., 2005; Scollon et al., 2017), especially craniosynostosis and renal anomalies (de Camargo et al., 2011), which we also observed in this study. Two recent studies used large population-based linkage data to investigate the association between birth defects and childhood cancer, showing that HB are more prevalent in children with birth defects (5.0% in this group compared to 1.3% in children with cancer; Schraw et al., 2020) and frequently occurs in association with craniosynostosis (Lupo et al., 2019). It is interesting that we found in this group of patients a marginally significant result suggesting a tendency of high-risk tumors to develop in patients who also have craniofacial defects. In addition, a large case-control study confirmed the association between HB and kidney and bladder abnormalities (Venkatramani et al., 2014). In this study, we also described the cases of two patients with Hirschsprung disease and HB, a condition with extensive genetic heterogeneity (Sergi et al., 2017; Heuckeroth, 2018); most Hirschsprung disease cases are sporadic (approximately 70%), while a smaller portion of patients have other associated congenital anomalies. However, pathogenic, or likely pathogenic variants mapped to morbid OMIM genes that could explain the syndromic phenotypes of our patients were not detected, reinforcing that, despite the long-standing knowledge about the intersection between biological pathways of cancer and development, the etiology of most of these associations in specific cases remains unknown.
There is a recognized sexual dimorphism in HB, with increased prevalence in males (Spector and Birch, 2012; Williams and Spector, 2019), and this study also reflected this tendency. We investigated whether this sex bias could be explained by an increased burden of rare damaging variants in one of the sexes; however, significant differences were not detected, suggesting that most likely environmental factors modulate the sexual dimorphism in HB prevalence.
In addition to SNV/indel variants, germline CNVs were already reported as causally related to cancer (Krepischi et al., 2012), but pathogenic CNVs affecting known/candidate CPGs were not detected in this study. One patient (P19) carried a gain of an entire Y chromosome, which is associated with an increased risk of learning disabilities, behavioral disturbances, and other clinical signs (Lalatta et al., 2012; Bardsley et al., 2013). Moreover, two 15q11.2 CNVs, considered risk factors for neurodevelopmental disorders (Butler, 2017; Rafi and Butler, 2020), were identified in two HB patients (P21 and P28). Only one patient (P28) presented an unspecific clinical feature related to the condition associated with its CNV (developmental delay). However, all these detected copy number alterations are known to present incomplete penetrance and variable expressivity. Importantly, germline 11p15.5 LOH, a common molecular mechanism resulting in epigenetic alterations known to increase the HB risk (Carrillo-Reixach et al., 2020; Fiala et al., 2020) as not detected in any of the patients.
Germline mutations have been detected in 8–18% of children and adolescents with cancer (Zhang et al., 2015; Akhavanfard et al., 2020; Capasso et al., 2020; Newman et al., 2021), and the most prevalent CPGs reported to be mutated are TP53, APC, BRCA2, NF1, PMS2, RB1, and RUNX1 (Zhang et al., 2015). In a very recent work about the spectrum of germline mutations in childhood cancer, most of the mutations (55%) were mapped to genes not previously associated with the patient’s tumor type (Newman et al., 2021); this work included a small number of HB patients (3 patients) with no pathogenic or likely pathogenic variants. Our findings revealed a high burden of damaging germline variants mapped to CPGs, with only 40% of these patients presenting a family history of cancer. Most of the variants affected autosomal-dominant CPGs not previously linked to HB, since only one mutation was detected in a known gene of HB development (APC), in a familial context of relatives diagnosed with FAP/colon cancer which segregated with the APC mutation. Interestingly, our data showed a clear predominance of CPGs related to gastrointestinal/colorectal cancer risk, such as familial adenomatous polyposis one and 2 (APC and MUTYH), nonpolyposis colorectal cancer hereditary types 1 (MSH2) and 6 (TGFBR2), Li-Fraumeni syndrome 2 (CHEK2), and Diamond-Blackfan anemia 1 (RPS19). Furthermore, the CPGs DROSHA and VHL, with putative P/LP variants, are known to be associated with renal cancer. Finally, a LoF variant was identified in the DNA repair gene ERCC5, in which homozygous mutations increase the skin cancer risk, and a damaging missense variant was mapped to FAH and causes an autosomal-recessive disorder characterized by progressive liver disease with increased liver cancer risk. Compared to previous pan-cancer studies, although not in HB patients, the findings of germline variants in APC, CHEK2, ERCC5, MSH2, MUTYH, and VHL were also detected in our work (Zhang et al., 2015; Ma et al., 2018a; Grobner et al., 2018; Newman et al., 2021).
We also observed numerous rare damaging VUS in known or candidate CPGs. Variants mapped to the ATM, BRCA2, COL7A1, DHCR7, DOCK8, FANCD2, FANCM, and GLI3 genes were observed in more than one HB patient. In particular, the same DHCR7 VUS (c.988G > A) was spotted in two patients (P09, born with only one functional kidney, and P26), both diagnosed with low-risk epithelial fetal HB; this gene encodes an enzyme that catalyzes the conversion of 7-dehydrocholesterol to cholesterol (Fitzky et al., 1998) and is a candidate for familial breast cancer in BRCA1-and BRCA2-negative breast cancer families (Shahi et al., 2019), in a recessive mode of action. As in our study, VUS mapped to MET, ATM, SLX4, and FANCA were also reported in HB patients in a recently published large study (Newman et al., 2021).
Strong support exists for the hypothesis that monoallelic CPG mutations confer an increased risk of cancer for adult carriers, while biallelic carriers would have a high risk of childhood cancer (Rahman and Scott, 2007). More recently, large studies have shown that variants in heterozygosity affecting recessive genes can also increase the predisposition to pediatric cancer (Savary et al., 2020), and a second hit in the tumor, such as a loss of heterozygosity or an inactivation of the second allele, has been observed in some of these patients. A highlight has been given to the role of germline monoallelic variants in cancer predisposition in genes involved in the recognition and repair of DNA damage, such as the ATM, PALB2, and Fanconi anemia genes (Thompson et al., 2005; Renwick et al., 2006; Antoniou et al., 2014; Helgason et al., 2015; Rebbeck et al., 2015; Esteban-Jurado et al., 2016; Esai Selvan et al., 2019). Since 1971 (Swift et al., 1971), it was proposed that individuals who were heterozygous for the Fanconi anemia genes might be at increased risk of cancer, and the premise would be that a modest reduction in the DNA repair efficiency could lead to tumor development. In our work, we detected relevant monoallelic variants in CPGs associated with recessive clinical conditions (FAH, ERCC5, and MUTYH) and in DNA repair genes (MSH2, CHEK2, ERCC5, and MUTYH). Recently, it has been shown that carriers of monoallelic loss of function MUTYH germline variants are at a higher risk of developing cancer, especially in tumors with frequent loss of heterozygosity events, such as adrenal adenocarcinoma, although the overall risk is still low (Barreiro et al., 2022); unfortunately, tumor sample of Patient P28, who carry a germline heterozygous MUTYH variant, was not available to check if loss of heterozygosity was present as a second hit, as suggested by authors. In addition, we observed a general enrichment of heterozygous germline damaging VUS affecting DNA repair genes related to Fanconi anemia, nucleotide and base excision repair, and homologous recombination repair. Sequencing the tumors of the carriers could contribute to changing the classification of such variants (Newman et al., 2021) and clarifying the role of these variants in cancer predisposition.
One of the most interesting findings of this study was the disclosure of rare germline damaging variants affecting genes linked to liver differentiation and function, such as FAH, that cause a recessive disorder related to hepatocellular carcinoma risk. Several rare damaging variants were observed in genes of the cytochrome P450 (CYP) family, including one pathogenic heterozygous variant in CYP21A2 (OMIM #201910 ADRENAL HYPERPLASIA, CONGENITAL, DUE TO 21-HYDROXYLASE DEFICIENCY, a recessive condition associated with testicular neoplasia in adults—P22) and two heterozygous pathogenic CYP1B1 variants (OMIM #617315, #231300; recessive conditions) detected in different patients (P03 and P30). Ten rare variants mapped to CYP1A1 (one LoF and nine missense variants), which is associated with primary liver metabolism, were present in eight patients (26%). CYP1A1 encodes a xenobiotic-metabolizing enzyme acting in the placenta (Yan et al., 2005; Suter et al., 2010; Stejskalova and Pavek, 2011), as well as in several drugs and compounds widely used in pharmacotherapy (Willis et al., 2018; Kapelyukh et al., 2019) or present in the diet (Ito et al., 2007). CYP1A1 expression is transcriptionally regulated through the AhR receptor (Kawashiro et al., 1998; Delescluse et al., 2000; Kamenickova et al., 2013) and various exogenous AhR binders, such as nitrosamines, polycyclic aromatic hydrocarbons, polychlorinated biphenyls, and halogenated dioxins, can be found in products from combustion processes, such as chimney soot, grilled food, and cigarette smoke, or as the product of incineration waste (Larigot et al., 2018; Lu et al., 2020; Schanz et al., 2020; Shahin et al., 2020). Elevated CYP1A1 activity through AhR activation in the placentas of female smokers has been associated with pregnancy complications, including low birth weight (Stejskalova et al., 2011), a known risk factor for HB. In our previous study (Aguiar et al., 2020), in which we investigated the mutational burden of HB, a mutational signature like the COSMIC 29 signature was detected, which was initially observed only in oral gengivo-buccal squamous cell carcinoma, which develops in individuals with the habit of chewing tobacco. Currently, this mutational signature has also been detected in patients with lung, thyroid, bladder, biliary, breast, pancreas, liver, and kidney cancer and can be linked to nitrosamine exposure (Alexandrov et al., 2015). It is interesting to notice that two of the patients carrying germline CYP1A1 variants developed HBs with predominance of S29. Therefore, we can speculate that a defective response during pregnancy to xenobiotic exposure, such as nitrosamines, could be linked to CYP1A1-damaging germline variants, increasing the risk for HB development, and resulting in the mutational signature we previously found in HB.
In conclusion, our major findings in HB patients were the detection of germline variants in CPGs associated with gastrointestinal/colorectal and renal cancer risk, not always linked to a familial history of cancer; the enrichment of monoallelic variants in CPGs and DNA repair genes; and a high frequency of birth defects. One limitation of our study is the small size of the group of patients, which precludes robust statistical analysis and demands further evaluation of germline data of larger cohorts; however, HB is indeed an ultrarare condition, and this is the largest assessment of germline variants in HB patients to date. Another point that deserves to be mentioned is that we did not investigate germline 11p15.5 epimutations, which are known to be associated with increased HB risk; this limitation was only partially overcome by the investigation of 11p15 regions of homozygosity in exome data. Finally, most studies of HB tumors are from North America, Europe, and Asia, and this work is pioneering in South America and contributes to elucidating the genetic architecture of HB risk. Further validation of the genes highlighted as HB germline risk factors in other cohorts can provide new insights regarding HB development.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors. Data is submitted to the https://www.ncbi.nlm.nih.gov/clinvar/ with identifier: SUB11168694 (SCV002103098 - SCV002103150).
Ethics Statement
Thirty children diagnosed with hepatoblastoma were enrolled in this study, and their parents, when available. Patients are from five different Brazilian`s institutes: A. C. Camargo Cancer Center, Adolescent and Child with Cancer Support Group (GRAACC), Pediatric Cancer Institute (ITACI), Baleia`s Hospital and Hospital da Criança Conceição. The Research Ethics Committee of the respective Institutions approved this research using these biological samples, and all samples were collected after informed signed consent was obtained from parents or legal guardians. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by TA, AT, JS and AK. The first draft of the manuscript was written by TA, AT and AK and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.
Funding
The present study was supported by grants from FAPESP (CEPID - Human Genome and Stem Cell Research Center 2013/08028-1; 2018/21047-9; fellowships: 2018/05961-2, 2016/04785-0, 2019/17423-8) and CNPq (141625/2016-3). The funders had no roles in the study design, data collection, and analysis, decision to publish, or preparation of the manuscript.
Conflict of Interest
Author IC was employed by Rede D’OR-São Luiz.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank the patients and their families who enrolled in this study.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.858396/full#supplementary-material
Abbreviations
CNV, Copy number variant; CPGs, Cancer predisposing genes; FA, Fanconi anemia; FAP, Familial adenomatous polyposis; FFPE, Formalin-fixed paraffin-embedded; HB, Hepatoblastoma; IGV, Integrated Genomics Viewer; LoF, Loss of function; LP, Likely pathogenic; P, Pathogenic; VUS, Variant of uncertain significance; WES, Whole-exome sequencing.
References
Agha, M. M., Williams, J. I., Marrett, L., To, T., Zipursky, A., and Dodds, L. (2005). Congenital Abnormalities and Childhood Cancer. Cancer 103, 1939–1948. doi:10.1002/cncr.20985
Aguiar, T. F. M., Rivas, M. P., Costa, S., Maschietto, M., Rodrigues, T., Sobral de Barros, J., et al. (2020). Insights into the Somatic Mutation Burden of Hepatoblastomas from Brazilian Patients. Front. Oncol. 10, 556. doi:10.3389/fonc.2020.00556
Akhavanfard, S., Padmanabhan, R., Yehia, L., Cheng, F., and Eng, C. (2020). Comprehensive Germline Genomic Profiles of Children, Adolescents and Young Adults with Solid Tumors. Nat. Commun. 11, 1–13. doi:10.1038/s41467-020-16067-1
Alexandrov, L. B., Jones, P. H., Wedge, D. C., Sale, J. E., Campbell, P. J., Nik-Zainal, S., et al. (2015). Clock-like Mutational Processes in Human Somatic Cells. Nat. Genet. 47, 1402–1407. doi:10.1038/ng.3441
Ansell, P., Mitchell, C. D., Roman, E., Simpson, J., Birch, J. M., and Eden, T. O. B. (2005). Relationships between Perinatal and Maternal Characteristics and Hepatoblastoma: a Report from the UKCCS. Eur. J. Cancer 41, 741–748. doi:10.1016/j.ejca.2004.10.024
Antoniou, A. C., Casadei, S., Heikkinen, T., Barrowdale, D., Pylkäs, K., Roberts, J., et al. (2014). Breast-cancer Risk in Families with Mutations in PALB2. N. Engl. J. Med. 371, 497–506. doi:10.1056/NEJMoa1400382
Baeg, G. H., Matsumine, A., Kuroda, T., Bhattacharjee, R. N., Miyashiro, I., Toyoshima, K., et al. (1995). The Tumour Suppressor Gene Product APC Blocks Cell Cycle Progression from G0/G1 to S Phase. EMBO J. 14 (22), 5618–5625. doi:10.1002/j.1460-2075.1995.tb00249.x
Bąk, A., Janiszewska, H., Junkiert-Czarnecka, A., Heise, M., Pilarska-Deltow, M., Laskowski, R., et al. (2014). A Risk of Breast Cancer in Women - Carriers of Constitutional CHEK2 Gene Mutations, Originating from the North - Central Poland. Hered. Cancer Clin. Pract. 12, 10. doi:10.1186/1897-4287-12-10
Bardsley, M. Z., Kowal, K., Levy, C., Gosek, A., Ayari, N., Tartaglia, N., et al. (2013). 47,XYY Syndrome: Clinical Phenotype and Timing of Ascertainment. J. Pediatr. 163, 1085–1094. doi:10.1016/J.JPEDS.2013.05.037
Barreiro, R. A. S., Sabbaga, J., Rossi, B. M., Achatz, M. I. W., Bettoni, F., Camargo, A. A., et al. (2022). Monoallelic Deleterious MUTYH Germline Variants as a Driver for Tumorigenesis. J. Pathol. 256 (2), 214–222. doi:10.1002/path.5829
Butler, M. G. (2017). Clinical and Genetic Aspects of the 15q11.2 BP1-BP2 Microdeletion Disorder. J. Intellect. Disabil. Res. 61, 568–579. doi:10.1111/jir.12382
Capasso, M., Montella, A., Tirelli, M., Maiorino, T., Cantalupo, S., and Iolascon, A. (2020). Genetic Predisposition to Solid Pediatric Cancers. Front. Oncol. 10, 2083. doi:10.3389/FONC.2020.590033
Capellini, A., Williams, M., Onel, K., and Huang, K.-L. (2021). The Functional Hallmarks of Cancer Predisposition Genes. Cmar Vol. 13, 4351–4357. doi:10.2147/CMAR.S311548
Carraro, D. M., Ramalho, R. F., Ramalho, R. F., and Maschietto, M. (2016). Gene Expression in Wilms Tumor: Disturbance of the Wnt Signaling Pathway and MicroRNA Biogenesis. Wilms Tumor 10, 149–162. doi:10.15586/CODON.WT.2016.CH10
Carrillo-Reixach, J., Torrens, L., Simon-Coma, M., Royo, L., Domingo-Sàbat, M., Abril-Fornaguera, J., et al. (2020). Epigenetic Footprint Enables Molecular Risk Stratification of Hepatoblastoma with Clinical Implications. J. Hepatol. 73 (2), 328–341. doi:10.1016/j.jhep.2020.03.025
Carta, R., Del Baldo, G., Miele, E., Po, A., Besharat, Z. M., Nazio, F., et al. (2020). Cancer Predisposition Syndromes and Medulloblastoma in the Molecular Era. Front. Oncol. 10, 566822. doi:10.3389/FONC.2020.566822/BIBTEX
Challis, D., Yu, J., Evani, U. S., Jackson, A. R., Paithankar, S., Coarfa, C., et al. (2012). An Integrative Variant Analysis Suite for Whole Exome Next-Generation Sequencing Data. BMC Bioinformatics 13, 8. doi:10.1186/1471-2105-13-8
Clarke, L., Fairley, S., Zheng-Bradley, X., Streeter, I., Perry, E., Lowy, E., et al. (2017). The International Genome Sample Resource (IGSR): A Worldwide Collection of Genome Variation Incorporating the 1000 Genomes Project Data. Nucleic Acids Res. 45, D854–D859. doi:10.1093/NAR/GKW829
Curia, M. C., Zuckermann, M., De Lellis, L., Catalano, T., Lattanzio, R., Aceto, G., et al. (2008). Sporadic Childhood Hepatoblastomas Show Activation of β-catenin, Mismatch Repair Defects and P53 Mutations. Mod. Pathol. 21, 7–14. doi:10.1038/modpathol.3800977
Czauderna, P., and Garnier, H. (2018). Hepatoblastoma: Current Understanding, Recent Advances, and Controversies. F1000Res 7, 53. doi:10.12688/f1000research.12239.1
Czauderna, P., Haeberle, B., Hiyama, E., Rangaswami, A., Krailo, M., Maibach, R., et al. (2016). The Children's Hepatic Tumors International Collaboration (CHIC): Novel Global Rare Tumor Database Yields New Prognostic Factors in Hepatoblastoma and Becomes a Research Model. Eur. J. Cancer 52, 92–101. doi:10.1016/j.ejca.2015.09.023
David Bick, M., Aziz, N., Bale, S., Richards, L., Soma Das, P., Gastier-Foster, J., et al. (2015). Standards and Guidelines for the Interpretation of Sequence Variants: A Joint Consensus Recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology Sue. Genet. Med. 17, 405–424. doi:10.1038/gim.2015.30.Standards
de Camargo, B., de Oliveira Ferreira, J. M., de Souza Reis, R., Ferman, S., de Oliveira Santos, M., and Pombo-de-Oliveira, M. S. (2011). Socioeconomic Status and the Incidence of Non-central Nervous System Childhood Embryonic Tumours in Brazil. BMC Cancer 11, 160. doi:10.1186/1471-2407-11-160
DeBaun, M. R., and Tucker, M. A. (1998). Risk of Cancer during the First Four Years of Life in Children from the Beckwith-Wiedemann Syndrome Registry. J. Pediatr. 132, 398–400. doi:10.1016/s0022-3476(98)70008-3
Delescluse, C., Lemaire, G., de Sousa, G., and Rahmani, R. (2000). Is CYP1A1 Induction Always Related to AHR Signaling Pathway? Toxicology 153, 73–82. doi:10.1016/s0300-483x(00)00305-x
Downing, J. R., Wilson, R. K., Zhang, J., Mardis, E. R., Pui, C.-H., Ding, L., et al. (2012). The Pediatric Cancer Genome Project. Nat. Genet. 44, 619–622. doi:10.1038/ng.2287
Esai Selvan, M., Klein, R. J., and Gümüş, Z. H. (2019). Rare, Pathogenic Germline Variants in Fanconi Anemia Genes Increase Risk for Squamous Lung Cancer. Clin. Cancer Res. 25, 1517–1525. doi:10.1158/1078-0432.CCR-18-2660
Esteban-Jurado, C., Franch-Expósito, S., Franch-Expósito, S., Muñoz, J., Ocaña, T., Carballal, S., et al. (2016). The Fanconi Anemia DNA Damage Repair Pathway in the Spotlight for Germline Predisposition to Colorectal Cancer. Eur. J. Hum. Genet. 24, 1501–1505. doi:10.1038/ejhg.2016.44
Fassihi, H., Sethi, M., Fawcett, H., Wing, J., Chandler, N., Mohammed, S., et al. (2016). Deep Phenotyping of 89 Xeroderma Pigmentosum Patients Reveals Unexpected Heterogeneity Dependent on the Precise Molecular Defect. Proc. Natl. Acad. Sci. U.S.A. 113 (9), E1236–E1245. doi:10.1073/pnas.1519444113
Felicio, P. S., Grasel, R. S., Campacci, N., Paula, A. E., Galvão, H. C. R., Torrezan, G. T., et al. (2021). Whole‐exome Sequencing of Non‐ BRCA1/BRCA2 Mutation Carrier Cases at High‐risk for Hereditary Breast/ovarian Cancer. Hum. Mutat. 42, 290–299. doi:10.1002/HUMU.24158
Feng, J., Polychronidis, G., Heger, U., Frongia, G., Mehrabi, A., and Hoffmann, K. (2019). Incidence Trends and Survival Prediction of Hepatoblastoma in Children: A Population-Based Study. Cancer Commun. 39, 62. doi:10.1186/s40880-019-0411-7
Fiala, E. M., Ortiz, M. V., Kennedy, J. A., Glodzik, D., Fleischut, M. H., Duffy, K. A., et al. (2020). 11p15.5 Epimutations in Children with Wilms Tumor and Hepatoblastoma Detected in Peripheral Blood. Cancer 126, 3114–3121. doi:10.1002/cncr.32907
Fitzky, B. U., Witsch-Baumgartner, M., Erdel, M., Lee, J. N., Paik, Y.-K., Glossmann, H., et al. (1998). Mutations in the 7-sterol Reductase Gene in Patients with the Smith-Lemli-Opitz Syndrome. Proc. Natl. Acad. Sci. 95, 8181–8186. doi:10.1073/PNAS.95.14.8181
Fuentes Fajardo, K. V., Adams, D., Mason, C. E., Sincan, M., Tifft, C., Toro, C., et al. (2012). Detecting False-Positive Signals in Exome Sequencing. Hum. Mutat. 33, 609–613. doi:10.1002/humu.22033
Gröbner, S. N., Worst, B. C., Worst, B. C., Weischenfeldt, J., Buchhalter, I., Kleinheinz, K., et al. (2018). The Landscape of Genomic Alterations across Childhood Cancers. Nature 555, 321–327. doi:10.1038/nature25480
Hampel, H., Bennett, R. L., Buchanan, A., Pearlman, R., and Wiesner, G. L. (2015). Guideline Development Group, American College of Medical Genetics and Genomics Professional Practice and Guidelines Committee and National Society of Genetic Counselors Practice Guidelines CommitteeA Practice Guideline from the American College of Medical Genetics and Genomics and the National Society of Genetic Counselors: Referral Indications for Cancer Predisposition Assessment. Genet. Med. 17, 70–87. doi:10.1038/gim.2014.147
Heck, J. E., Meyers, T. J., Lombardi, C., Park, A. S., Cockburn, M., Reynolds, P., et al. (2013). Case-control Study of Birth Characteristics and the Risk of Hepatoblastoma. Cancer Epidemiol. 37, 390–395. doi:10.1016/j.canep.2013.03.004
Helgason, H., Rafnar, T., Olafsdottir, H. S., Jonasson, J. G., Sigurdsson, A., Stacey, S. N., et al. (2015). Loss-of-function Variants in ATM Confer Risk of Gastric Cancer. Nat. Genet. 47, 906–910. doi:10.1038/ng.3342
Heuckeroth, R. O. (2018). Hirschsprung Disease - Integrating Basic Science and Clinical Medicine to Improve Outcomes. Nat. Rev. Gastroenterol. Hepatol. 15, 152–167. doi:10.1038/nrgastro.2017.149
Hirschman, B. A., Pollock, B. H., and Tomlinson, G. E. (2005). The Spectrum of APC Mutations in Children with Hepatoblastoma from Familial Adenomatous Polyposis Kindreds. J. Pediatr. 147, 263–266. doi:10.1016/j.jpeds.2005.04.019
Institute Nacional de Câncer José Alencar Gomes da Silva (INCA) (2016). Incidência, mortalidade e morbidade hospitalar por câncer em crianças, adolescentes e adultos jovens no Brasil: informações dos registros de câncer e do sistema de mortalidade. Rio de Janeiro: Coordenação de Prevenção e VigilânciaINCA.
Ito, S., Chen, C., Satoh, J., Yim, S., and Gonzalez, F. J. (2007). Dietary Phytochemicals Regulate Whole-Body CYP1A1 Expression through an Arylhydrocarbon Receptor Nuclear Translocator-dependent System in Gut. J. Clin. Invest. 117, 1940–1950. doi:10.1172/JCI31647
Johnson, K. J., Lee, J. M., Ahsan, K., Padda, H., Feng, Q., Partap, S., et al. (2017). Pediatric Cancer Risk in Association with Birth Defects: A Systematic Review. PLoS One 12, e0181246. doi:10.1371/journal.pone.0181246
Kamenickova, A., Anzenbacherova, E., Pavek, P., Soshilov, A. A., and Denison, M. S. (2013). Effects of Anthocyanins on the AhR-Cyp1a1 Signaling Pathway in Human Hepatocytes and Human Cancer Cell Lines. Toxicol. Lett. 221, 1–8. doi:10.1016/J.TOXLET.2013.05.007
Kamien, B. A., and Gabbett, M. T. (2009). Aicardi Syndrome Associated with Hepatoblastoma and Pulmonary Sequestration. Am. J. Med. Genet. 149A, 1850–1852. doi:10.1002/ajmg.a.32985
Kapelyukh, Y., Henderson, C. J., Scheer, N., Rode, A., and Wolf, C. R. (2019). Defining the Contribution of CYP1A1 and CYP1A2 to Drug Metabolism Using Humanized CYP1A1/1A2 and Cyp1a1/Cyp1a2 Knockout Mice. Drug Metab. Dispos. 47, 907–918. doi:10.1124/DMD.119.087718
Karczewski, K. J., Francioli, L. C., Tiao, G., Cummings, B. B., Alföldi, J., Wang, Q., et al. (2020). The Mutational Constraint Spectrum Quantified from Variation in 141,456 Humans. Nature 581, 434–443. doi:10.1038/s41586-020-2308-7
Kawashiro, T., Yamashita, K., Zhao, X. J., Koyama, E., Tani, M., Chiba, K., et al. (1998). A Study on the Metabolism of Etoposide and Possible Interactions with Antitumor or Supporting Agents by Human Liver Microsomes. J. Pharmacol. Exp. Ther. 286, 1294–1300.
Kim, S. Y., Jung, S.-H., Kim, M. S., Han, M.-R., Park, H.-C., Jung, E. S., et al. (2017). Genomic Profiles of a Hepatoblastoma from a Patient with Beckwith-Wiedemann Syndrome with Uniparental Disomy on Chromosome 11p15 and Germline Mutation of APC and PALB2. Oncotarget 8, 91950–91957. doi:10.18632/oncotarget.20515
Kindler, O., Quehenberger, F., Benesch, M., and Seidel, M. G. (2018). The Iceberg Map of Germline Mutations in Childhood Cancer. Curr. Opin. Pediatr. 30, 855–863. doi:10.1097/MOP.0000000000000680
Köhler, S., Carmody, L., Vasilevsky, N., Jacobsen, J. O. B., Danis, D., Gourdine, J. P., et al. (2019). Expansion of the Human Phenotype Ontology (HPO) Knowledge Base and Resources. Nucleic Acids Res. 47, D1018–D1027. doi:10.1093/nar/gky1105
Kopanos, C., Tsiolkas, V., Kouris, A., Chapple, C. E., Albarca Aguilera, M., Meyer, R., et al. (2019). VarSome: the Human Genomic Variant Search Engine. Bioinformatics 35, 1978–1980. doi:10.1093/bioinformatics/bty897
Krepischi, A. C. V., Pearson, P. L., and Rosenberg, C. (2012). Germline Copy Number Variations and Cancer Predisposition. Future Oncol. 8, 441–450. doi:10.2217/fon.12.34
Krivit, W., and Good, R. A. (1957). Simultaneous Occurrence of Mongolism and Leukemia. AMA Am. J. Dis. Child. 94, 289. doi:10.1001/archpedi.1957.04030040075012
Lalatta, F., Folliero, E., Cavallari, U., Di Segni, M., Gentilin, R., Quagliarini, D., et al. (2012). Early Manifestations in a Cohort of Children Prenatally Diagnosed with 47,XYY. Role of Multidisciplinary Counseling for Parental Guidance and Prevention of Aggressive Behavior. Ital. J. Pediatr. 38, 52. doi:10.1186/1824-7288-38-52
Larigot, L., Juricek, L., Dairou, J., and Coumoul, X. (2018). AhR Signaling Pathways and Regulatory Functions. Biochimie Open 7, 1–9. doi:10.1016/J.BIOPEN.2018.05.001
Li, H. (2013). Aligning Sequence Reads, Clone Sequences and Assembly Contigs with BWA-MEM. Available at: http://github.com/lh3/bwa (Accessed May 18, 2020).
Lu, J., Shang, X., Zhong, W., Xu, Y., Shi, R., and Wang, X. (2020). New Insights of CYP1A in Endogenous Metabolism: a Focus on Single Nucleotide Polymorphisms and Diseases. Acta Pharmaceutica Sinica B 10, 91–104. doi:10.1016/J.APSB.2019.11.016
Luft, F. (2010). The Rise of a Ribosomopathy and Increased Cancer Risk. J. Mol. Med. 88, 1–3. doi:10.1007/s00109-009-0570-0
Lupo, P. J., Schraw, J. M., Desrosiers, T. A., Nembhard, W. N., Langlois, P. H., Canfield, M. A., et al. (2019). Association between Birth Defects and Cancer Risk Among Children and Adolescents in a Population-Based Assessment of 10 Million Live Births. JAMA Oncol. 5 (8), 1150–1158. doi:10.1001/jamaoncol.2019.1215
Ma, X., Beeghly-Fadiel, A., Lu, W., Shi, J., Xiang, Y.-B., Cai, Q., et al. (2012). Pathway Analyses Identify TGFBR2 as Potential Breast Cancer Susceptibility Gene: Results from a Consortium Study Among Asians. Cancer Epidemiol. Biomarkers Prev. 21, 1176–1184. doi:10.1158/1055-9965.EPI-12-0118
Ma, X., Liu, Y., Liu, Y., Alexandrov, L. B., Edmonson, M. N., Gawad, C., et al. (2018a). Pan-cancer Genome and Transcriptome Analyses of 1,699 Paediatric Leukaemias and Solid Tumours. Nature 555, 371–376. doi:10.1038/nature25795
McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., et al. (2010). The Genome Analysis Toolkit: A MapReduce Framework for Analyzing Next-Generation DNA Sequencing Data. Genome Res. 20, 1297–1303. doi:10.1101/gr.107524.110
Meyers, R. L., Maibach, R., Hiyama, E., Häberle, B., Krailo, M., Rangaswami, A., et al. (2017). Risk-stratified Staging in Paediatric Hepatoblastoma: a Unified Analysis from the Children's Hepatic Tumors International Collaboration. Lancet Oncol. 18, 122–131. doi:10.1016/S1470-2045(16)30598-8
Narod, S. A., Hawkins, M. M., Robertson, C. M., and Stiller, C. A. (1997). Congenital Anomalies and Childhood Cancer in Great Britain. Am. J. Hum. Genet. 60, 474–485.
Naslavsky, M. S., Scliar, M. O., Yamamoto, G. L., Wang, J. Y. T., Zverinova, S., Karp, T., et al. (2020). Whole-genome Sequencing of 1,171 Elderly Admixed Individuals from the Largest Latin American metropolis (São Paulo, Brazil). bioRxiv 10, 298026. doi:10.1101/2020.09.15.298026
Newman, S., Nakitandwe, J., Kesserwan, C. A., Azzato, E. M., Wheeler, D. A., Rusch, M., et al. (2021). Genomes for Kids: The Scope of Pathogenic Mutations in Pediatric Cancer Revealed by Comprehensive DNA and RNA Sequencing. Cancer Discov. 11, 3008–3027. doi:10.1158/2159-8290.CD-20-1631
Norwood, M. S., Lupo, P. J., Chow, E. J., Scheurer, M. E., Plon, S. E., Danysh, H. E., et al. (2017). Childhood Cancer Risk in Those with Chromosomal and Non-chromosomal Congenital Anomalies in Washington State: 1984-2013. PLoS One 12, e0179006. doi:10.1371/journal.pone.0179006
Oliveira, D., Leal, G. F., Sertié, A. L., Caires Jr, L. C., Goulart, E., Musso, C. M., et al. (2018). 10q23.31 Microduplication Encompassing PTEN Decreases mTOR Signalling Activity and Is Associated with Autosomal Dominant Primary Microcephaly. J. Med. Genet. 56, 543–547. doi:10.1136/jmedgenet-2018-105471
Oscanoa, J., Sivapalan, L., Gadaleta, E., Dayem Ullah, A. Z., Lemoine, N. R., and Chelala, C. (2020). SNPnexus: a Web Server for Functional Annotation of Human Genome Sequence Variation (2020 Update). Nucleic Acids Res. 48, W185–W192. doi:10.1093/nar/gkaa420
Pateva, I. B., Egler, R. A., and Stearns, D. S. (2017). Hepatoblastoma in an 11-Year-Old. Medicine (Baltimore) 96, e5858. doi:10.1097/MD.0000000000005858
Pereira, E. M., Marion, R., Ramesh, K. H., Kim, J. S., Ewart, M., and Ricafort, R. (2012). Hepatoblastoma in a Mosaic Trisomy 18 Patient. J. Pediatr. Hematol. Oncol. 34, e145–e148. doi:10.1097/MPH.0b013e3182459ee8
Pinto, R. B., Ramos, A. R. L., Backes, A. N., Santos, B. J. d., Provenzi, V. O., Carbonera, M. R., et al. (2016). Hirschsprung Disease and Hepatoblastoma: Case Report of a Rare Association. Sao Paulo Med. J. 134, 171–175. doi:10.1590/1516-3180.2014.9200311
Rafi, S. K., and Butler, M. G. (2020). The 15q11.2 BP1-BP2 Microdeletion (Burnside-Butler) Syndrome: In Silico Analyses of the Four Coding Genes Reveal Functional Associations with Neurodevelopmental Disorders. Ijms 21, 3296. doi:10.3390/ijms21093296
Rahman, N., and Scott, R. H. (2007). Cancer Genes Associated with Phenotypes in Monoallelic and Biallelic Mutation Carriers: New Lessons from Old Players. Hum. Mol. Genet. 16, R60–R66. doi:10.1093/hmg/ddm026
Rebbeck, T. R., Mitra, N., Wan, F., Sinilnikova, O. M., Healey, S., McGuffog, L., et al. (2015). Association of Type and Location of BRCA1 and BRCA2 Mutations with Risk of Breast and Ovarian Cancer. JAMA 313, 1347–1361. doi:10.1001/jama.2014.5985
Renwick, A., Thompson, D., Thompson, D., Seal, S., Kelly, P., Chagtai, T., et al. (2006). ATM Mutations that Cause Ataxia-Telangiectasia Are Breast Cancer Susceptibility Alleles. Nat. Genet. 38, 873–875. doi:10.1038/ng1837
Ries, L., Smith, M., Gurney, J., Linet, M., Tamra, T., Young, J., et al. (1999). Cancer Incidence and Survival Among Children and Adolescents: United States SEER Program 1975-1995. Natl. Cancer Inst. SEER Program, 113, 2575–2596.
Saletta, F., Dalla Pozza, L., and Byrne, J. A. (2015). Genetic Causes of Cancer Predisposition in Children and Adolescents. Transl. Pediatr. 4, 67–75. doi:10.3978/j.issn.2224-4336.2015.04.08
Sambrook, J., and Russell, D. W. (2006). Purification of Nucleic Acids by Extraction with Phenol:chloroform. Cold Spring Harb Protoc. 2006, pdb.prot4455. doi:10.1101/pdb.prot4455
Savary, C., Kim, A., Lespagnol, A., Gandemer, V., Pellier, I., Andrieu, C., et al. (2020). Depicting the Genetic Architecture of Pediatric Cancers through an Integrative Gene Network Approach. Sci. Rep. 10, 1–15. doi:10.1038/s41598-020-58179-0
Schanz, O., Chijiiwa, R., Cengiz, S. C., Majlesain, Y., Weighardt, H., Takeyama, H., et al. (2020). Dietary AhR Ligands Regulate AhRR Expression in Intestinal Immune Cells and Intestinal Microbiota Composition. Ijms 21, 3189. doi:10.3390/IJMS21093189
Schmidt, L. S., and Linehan, W. M. (2016). Genetic Predisposition to Kidney Cancer. Semin. Oncol. 43 (5), 566–574. doi:10.1053/j.seminoncol.2016.09.001
Schraw, J. M., Desrosiers, T. A., Nembhard, W. N., Langlois, P. H., Meyer, R. E., Canfield, M. A., et al. (2020). Cancer Diagnostic Profile in Children with Structural Birth Defects: An Assessment in 15,000 Childhood Cancer Cases. Cancer 126 (15), 3483–3492. doi:10.1002/cncr.32982
Scollon, S., Anglin, A. K., Thomas, M., Turner, J. T., and Wolfe Schneider, K. (2017). A Comprehensive Review of Pediatric Tumors and Associated Cancer Predisposition Syndromes. J. Genet. Counsel 26, 387–434. doi:10.1007/s10897-017-0077-8
Sergi, C. M., Caluseriu, O., McColl, H., and Eisenstat, D. D. (2017). Hirschsprung's Disease: Clinical Dysmorphology, Genes, Micro-RNAs, and Future Perspectives. Pediatr. Res. 81, 177–191. doi:10.1038/pr.2016.202
Shahi, R. B., De Brakeleer, S., Caljon, B., Pauwels, I., Bonduelle, M., Joris, S., et al. (2019). Identification of Candidate Cancer Predisposing Variants by Performing Whole-Exome Sequencing on index Patients from BRCA1 and BRCA2-Negative Breast Cancer Families. BMC Cancer 19, 313. doi:10.1186/s12885-019-5494-7
Shahin, N. N., Abd-Elwahab, G. T., Tawfiq, A. A., and Abdelgawad, H. M. (2020). Potential Role of Aryl Hydrocarbon Receptor Signaling in Childhood Obesity. Biochim. Biophys. Acta (Bba) - Mol. Cel Biol. Lipids 1865, 158714. doi:10.1016/J.BBALIP.2020.158714
Siegel, R. L., Miller, K. D., and Jemal, A. (2017). Cancer Statistics, 2017. CA: A Cancer J. Clinicians 67, 7–30. doi:10.3322/caac.21387
Southey, M. C., Goldgar, D. E., Winqvist, R., Pylkäs, K., Couch, F., Tischkowitz, M., et al. (2016). PALB2, CHEK2 and ATM Rare Variants and Cancer Risk: Data from COGS. J. Med. Genet. 53 (12), 800–811. doi:10.1136/jmedgenet-2016-103839
Spector, L. G., and Birch, J. (2012). The Epidemiology of Hepatoblastoma. Pediatr. Blood Cancer 59, 776–779. doi:10.1002/pbc.24215
Spector, L. G., Brown, M. B., Wantman, E., Letterie, G. S., Toner, J. P., Doody, K., et al. (2019). Association of In Vitro Fertilization with Childhood Cancer in the United States. JAMA Pediatr. 173, e190392. doi:10.1001/JAMAPEDIATRICS.2019.0392
Stejskalova, L., Dvorak, Z., and Pavek, P. (2011). Endogenous and Exogenous Ligands of Aryl Hydrocarbon Receptor: Current State of Art. Curr. Drug Metab. 12, 198–212. doi:10.2174/138920011795016818
Stejskalova, L., and Pavek, P. (2011). The Function of Cytochrome P450 1A1 Enzyme (CYP1A1) and Aryl Hydrocarbon Receptor (AhR) in the Placenta. Cpb 12, 715–730. doi:10.2174/138920111795470994
Stelzer, G., Plaschkes, I., Oz-Levi, D., Alkelai, A., Olender, T., Zimmerman, S., et al. (2016). VarElect: The Phenotype-Based Variation Prioritizer of the GeneCards Suite. BMC Genomics 17, 444. doi:10.1186/s12864-016-2722-2
Stiller, C. A., Pritchard, J., and Steliarova-Foucher, E. (2006). Liver Cancer in European Children: Incidence and Survival, 1978-1997. Report from the Automated Childhood Cancer Information System Project. Eur. J. Cancer 42, 2115–2123. doi:10.1016/j.ejca.2006.05.011
Suter, M., Abramovici, A., Showalter, L., Hu, M., Shope, C. D., Varner, M., et al. (2010). In Utero tobacco Exposure Epigenetically Modifies Placental CYP1A1 Expression. Metabolism 59, 1481–1490. doi:10.1016/J.METABOL.2010.01.013
Sweet-Cordero, E. A., and Biegel, J. A. (2019). The Genomic Landscape of Pediatric Cancers: Implications for Diagnosis and Treatment. Science 363, 1170–1175. doi:10.1126/science.aaw3535
Swift, M., Zimmerman, D., and McDonough, E. R. (1971). Squamous Cell Carcinomas in Fanconi's Anemia. JAMA 216, 325–326. doi:10.1001/jama.216.2.325
Tanimura, M., Matsui, I., Abe, J., Ikeda, H., Kobayashi, N., Ohira, M., et al. (1998). Increased Risk of Hepatoblastoma Among Immature Children with a Lower Birth Weight. Cancer Res. 58, 3032–3035.
Thompson, D., Duedal, S., Kirner, J., McGuffog, L., Last, J., Reiman, A., et al. (2005). Cancer Risks and Mortality in Heterozygous ATM Mutation Carriers. J. Natl. Cancer Inst. 97, 813–822. doi:10.1093/jnci/dji141
Turcotte, L. M., Georgieff, M. K., Ross, J. A., Feusner, J. H., Tomlinson, G. E., Malogolowkin, M. H., et al. (2014). Neonatal Medical Exposures and Characteristics of Low Birth Weight Hepatoblastoma Cases: A Report from the Children's Oncology Group. Pediatr. Blood Cancer 61, 2018–2023. doi:10.1002/pbc.25128
Vahteristo, P., Bartkova, J., Eerola, H., Syrjäkoski, K., Ojala, S., Kilpivaara, O., et al. (2002). A CHEK2 Genetic Variant Contributing to a Substantial Fraction of Familial Breast Cancer. Am. J. Hum. Genet. 71 (2), 432–438. doi:10.1002/ajmg.a.36638
Venkatramani, R., Spector, L. G., Georgieff, M., Tomlinson, G., Krailo, M., Malogolowkin, M., et al. (2014). Congenital Abnormalities and Hepatoblastoma: a Report from the Children’s Oncology Group (COG) and the Utah Population Database (UPDB). Am. J. Med. Genet. A. 164A, 2250–2255. doi:10.1002/ajmg.a.36638
Wegert, J., Ishaque, N., Vardapour, R., Geörg, C., Gu, Z., Bieg, M., et al. (2015). Mutations in the SIX1/2 Pathway and the DROSHA/DGCR8 miRNA Microprocessor Complex Underlie High-Risk Blastemal Type Wilms Tumors. Cancer Cell 27, 298–311. doi:10.101610.1016/j.ccell.2015.01.002
Williams, L. A., and Spector, L. G. (2019). Survival Differences between Males and Females Diagnosed with Childhood Cancer. JNCI Cancer Spectr. 3, pkz032. doi:10.1093/jncics/pkz032
Willis, A. J., Indra, R., Wohak, L. E., Sozeri, O., Feser, K., Mrizova, I., et al. (2018). The Impact of Chemotherapeutic Drugs on the CYP1A1-Catalysed Metabolism of the Environmental Carcinogen Benzo[a]pyrene: Effects in Human Colorectal HCT116 TP53(+/+), TP53(+/−) and TP53(−/−) Cells. Toxicology 398-399, 1–12. doi:10.1016/J.TOX.2018.02.006
Xu, Y., and Pasche, B. (2007). TGF-β Signaling Alterations and Susceptibility to Colorectal Cancer. Hum. Mol. Genet. 16, R14–R20. doi:10.1093/HMG/DDL486
Yan, Y.-E., Wang, H., and Feng, Y.-H. (2005). Alterations of Placental Cytochrome P450 1A1 and P-Glycoprotein in Tobacco-Induced Intrauterine Growth Retardation in Rats1. Acta Pharmacologica Sinica 26, 1387–1394. doi:10.1111/j.1745-7254.2005.00209.x
Keywords: hepatoblastoma, cancer predisposition, DNA repair, MSH2, VHL, ERCC5, MUTYH, CYP1A1
Citation: Aguiar T, Teixeira A, Scliar MO, Sobral de Barros J, Lemes RB, Souza S, Tolezano G, Santos F, Tojal I, Cypriano M, Caminada de Toledo SR, Valadares E, Borges Pinto R, Pinto Artigalas OA, Caetano de Aguirre Neto J, Novak E, Cristofani LM, Miura Sugayama SM, Odone V, Cunha IW, Lima da Costa CM, Rosenberg C and Krepischi A (2022) Unraveling the Genetic Architecture of Hepatoblastoma Risk: Birth Defects and Increased Burden of Germline Damaging Variants in Gastrointestinal/Renal Cancer Predisposition and DNA Repair Genes. Front. Genet. 13:858396. doi: 10.3389/fgene.2022.858396
Received: 19 January 2022; Accepted: 08 March 2022;
Published: 12 April 2022.
Edited by:
Augusto Rojas-Martinez, Tecnológico de Monterrey, MexicoReviewed by:
Päivi Peltomäki, University of Helsinki, FinlandAdam de Smith, University of Southern California, United States
Copyright © 2022 Aguiar, Teixeira, Scliar, Sobral de Barros, Lemes, Souza, Tolezano, Santos, Tojal, Cypriano, Caminada de Toledo, Valadares, Borges Pinto, Pinto Artigalas, Caetano de Aguirre Neto, Novak, Cristofani, Miura Sugayama, Odone, Cunha, Lima da Costa, Rosenberg and Krepischi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ana Krepischi, ana.krepischi@ib.usp.br
†These authors have contributed equally to this work