Clinical Characterization and Underlying Genetic Findings in Brazilian Patients with Syndromic Microcephaly Associated with Neurodevelopmental Disorders.

Microcephaly is characterized by an occipitofrontal circumference at least two standard deviations below the mean for age and sex. Neurodevelopmental disorders (NDD) are commonly associated with microcephaly, due to perturbations in brain development and functioning. Given the extensive genetic heterogeneity of microcephaly, managing patients is hindered by the broad spectrum of diagnostic possibilities that exist before conducting molecular testing. We investigated the genetic basis of syndromic microcephaly accompanied by NDD in a Brazilian cohort of 45 individuals and characterized associated clinical features, as well as evaluated the effectiveness of whole-exome sequencing (WES) as a diagnostic tool for this condition. Patients previously negative for pathogenic copy number variants underwent WES, which was performed using a trio approach for isolated index cases (n = 31), only the index in isolated cases with parental consanguinity (n = 8) or affected siblings in familial cases (n = 3). Pathogenic/likely pathogenic variants were identified in 19 families (18 genes) with a diagnostic yield of approximately 45%. Nearly 86% of the individuals had global developmental delay/intellectual disability and 51% presented with behavioral disturbances. Additional frequent clinical features included facial dysmorphisms (80%), brain malformations (67%), musculoskeletal (71%) or cardiovascular (47%) defects, and short stature (54%). Our findings unraveled the underlying genetic basis of microcephaly in half of the patients, demonstrating a high diagnostic yield of WES for microcephaly and reinforcing its genetic heterogeneity. We expanded the phenotypic spectrum associated with the condition and identified a potentially novel gene (CCDC17) for congenital microcephaly.

Germline mutations in cancer predisposition genes among pediatric patients with cancer and congenital anomalies.

BACKGROUND: Childhood cancer has a poorly known etiology, and investigating the underlying genetic background may provide novel insights. A recognized association exists between non-chromosomal birth defects and childhood cancer susceptibility. METHODS: We performed whole-exome sequencing and chromosomal microarray analysis in a cohort of childhood cancer (22 individuals, 50% with congenital anomalies) to unravel deleterious germline variants. RESULTS: A diagnostic yield of 14% was found, encompassing heterozygous variants in bona fide dominant Cancer Predisposition Genes (CPGs). Considering candidate and recessive CPGs harboring monoallelic variants, which were also deemed to play a role in the phenotype, the yield escalated to 45%. Most of the deleterious variants were mapped in genes not conventionally linked to the patient's tumor type. Relevant findings were detected in 55% of the syndromic individuals, mostly variants potentially underlying both phenotypes. CONCLUSION: We uncovered a remarkable prevalence of germline deleterious CPG variants, highlighting the significance of a comprehensive genetic analysis in pediatric cancer, especially when coupled with additional clinical signs. Moreover, our findings emphasized the potential for oligogenic inheritance, wherein multiple genes synergistically increase cancer risk. Lastly, our investigation unveiled potentially novel genotype-phenotype associations, such as SETD5 in neuroblastoma, KAT6A in gliomas, JAG1 in hepatoblastomas, and TNFRSF13B in Langerhans cell histiocytosis. IMPACT: Novel gene-phenotype associations and candidate genes for pediatric cancer were unraveled, such as KAT6A in gliomas, SETD5 in neuroblastoma, JAG1 in hepatoblastomas, and TNFRSF13B in Langerhans cell histiocytosis. Our analysis revealed a high frequency of deleterious germline variants, particularly in cases accompanied by additional clinical signs, highlighting the importance of a comprehensive genetic evaluation in childhood cancer. Our findings also underscored the potential for oligogenic inheritance in pediatric cancer risk. Understanding the cancer etiology is crucial for genetic counseling, often influencing therapeutic decisions and offering valuable insights into molecular targets for the development of oncological therapies.

Optimal population-specific HLA imputation with dimension reduction.

Human genomics has quickly evolved, powering genome-wide association studies (GWASs). SNP-based GWASs cannot capture the intense polymorphism of HLA genes, highly associated with disease susceptibility. There are methods to statistically impute HLA genotypes from SNP-genotypes data, but lack of diversity in reference panels hinders their performance. We evaluated the accuracy of the 1000 Genomes data as a reference panel for imputing HLA from admixed individuals of African and European ancestries, focusing on (a) the full dataset, (b) 10 replications from 6 populations, and (c) 19 conditions for the custom reference panels. The full dataset outperformed smaller models, with a good F1-score of 0.66 for HLA-B. However, custom models outperformed the multiethnic or population models of similar size (F1-scores up to 0.53, against up to 0.42). We demonstrated the importance of using genetically specific models for imputing populations, which are currently underrepresented in public datasets, opening the door to HLA imputation for every genetic population.

Analysis of the Mutational Landscape of Osteosarcomas Identifies Genes Related to Metastasis and Prognosis and Disrupted Biological Pathways of Immune Response and Bone Development.

Osteosarcoma (OS) is the most prevalent type of bone tumor, but slow progress has been achieved in disentangling the full set of genomic events involved in its initiation and progression. We assessed by NGS the mutational spectrum of 28 primary OSs from Brazilian patients, and identified 445 potentially deleterious SNVs/indels and 1176 copy number alterations (CNAs). TP53 was the most recurrently mutated gene, with an overall rate of ~60%, considering SNVs/indels and CNAs. The most frequent CNAs (~60%) were gains at 1q21.2q21.3, 6p21.1, and 8q13.3q24.22, and losses at 10q26 and 13q14.3q21.1. Seven cases presented CNA patterns reminiscent of complex events (chromothripsis and chromoanasynthesis). Putative RB1 and TP53 germline variants were found in five samples associated with metastasis at diagnosis along with complex genomic patterns of CNAs. PTPRQ, KNL1, ZFHX4, and DMD alterations were prevalent in metastatic or deceased patients, being potentially indicative of poor prognosis. TNFRSF11B, involved in skeletal system development and maintenance, emerged as a candidate for osteosarcomagenesis due to its biological function and a high frequency of copy number gains. A protein-protein network enrichment highlighted biological pathways involved in immunity and bone development. Our findings reinforced the high genomic OS instability and heterogeneity, and led to the identification of novel disrupted genes deserving further evaluation as biomarkers due to their association with poor outcomes.

Three generation families: Analysis of de novo variants in autism.

De novo variants (DNVs) analysis has proven to be a powerful approach to gene discovery in Autism Spectrum Disorder (ASD), which has not yet been shown in a Brazilian ASD cohort. The relevance of inherited rare variants has also been suggested, particularly in oligogenic models. We hypothesized that three-generation analyses of DNVs could provide new insights into the relevance of de novo and inherited variants across generations. To accomplish this goal, we performed whole-exome sequencing of 33 septet families composed of probands, parents, and grandparents (n = 231 individuals) and compared DNV rates (DNVr) between generations and those from two control cohorts. The DNVr in the probands (DNVr = 1.16) was marginally higher than in parents (DNVr = 0.60; p = 0.054), and in controls (DNVr = 0.68; p = 0.035, congenital heart disorder and DNVr = 0.70; p = 0.047, unaffected ASD siblings from Simons Simplex Collection). Moreover, most of the DNVs were found to have paternal origin in both generations (84.6%). Finally, we observed that 40% (6/15) of the DNVs in parents transmitted for probands are in ASD or ASD candidate genes, representing recently emerged risk variants to ASD in their families and suggest ZNF536, MSL2 and HDAC9 as ASD candidate genes. We did not observe an enrichment of risk variants nor sex bias of transmitted variants in the three generations, that can be due to sample size. These results further reinforce the relevance of de novo variants in ASD.

Characterization of pharmacogenomic variants in a Brazilian admixed cohort of elderly individuals based on whole-genome sequencing data.

Introduction: Research in the field of pharmacogenomics (PGx) aims to identify genetic variants that modulate response to drugs, through alterations in their pharmacokinetics (PK) or pharmacodynamics (PD). The distribution of PGx variants differs considerably among populations, and whole-genome sequencing (WGS) plays a major role as a comprehensive approach to detect both common and rare variants. This study evaluated the frequency of PGx markers in the context of the Brazilian population, using data from a population-based admixed cohort from Sao Paulo, Brazil, which includes variants from WGS of 1,171 unrelated, elderly individuals. Methods: The Stargazer tool was used to call star alleles and structural variants (SVs) from 38 pharmacogenes. Clinically relevant variants were investigated, and the predicted drug response phenotype was analyzed in combination with the medication record to assess individuals potentially at high-risk of gene-drug interaction. Results: In total, 352 unique star alleles or haplotypes were observed, of which 255 and 199 had a frequency < 0.05 and < 0.01, respectively. For star alleles with frequency > 5% (n = 97), decreased, loss-of-function and unknown function accounted for 13.4%, 8.2% and 27.8% of alleles or haplotypes, respectively. Structural variants (SVs) were identified in 35 genes for at least one individual, and occurred with frequencies >5% for CYP2D6, CYP2A6, GSTM1, and UGT2B17. Overall 98.0% of the individuals carried at least one high risk genotype-predicted phenotype in pharmacogenes with PharmGKB level of evidence 1A for drug interaction. The Electronic Health Record (EHR) Priority Result Notation and the cohort medication registry were combined to assess high-risk gene-drug interactions. In general, 42.0% of the cohort used at least one PharmGKB evidence level 1A drug, and 18.9% of individuals who used PharmGKB evidence level 1A drugs had a genotype-predicted phenotype of high-risk gene-drug interaction. Conclusion: This study described the applicability of next-generation sequencing (NGS) techniques for translating PGx variants into clinically relevant phenotypes on a large scale in the Brazilian population and explores the feasibility of systematic adoption of PGx testing in Brazil.

Corrigendum: Follow-up of young adult monozygotic twins after simultaneous critical coronavirus disease 2019: a case report.

[This corrects the article DOI: 10.3389/fmed.2022.1008585.].

DNA methylation patterns suggest the involvement of DNMT3B and TET1 in osteosarcoma development.

DNA methylation may be involved in the development of osteosarcomas. Osteosarcomas commonly arise during the bone growth and remodeling in puberty, making it plausible to infer the involvement of epigenetic alterations in their development. As a highly studied epigenetic mechanism, we investigated DNA methylation and related genetic variants in 28 primary osteosarcomas aiming to identify deregulated driver alterations. Methylation and genomic data were obtained using the Illumina HM450K beadchips and the TruSight One sequencing panel, respectively. Aberrant DNA methylation was spread throughout the osteosarcomas genomes. We identified 3146 differentially methylated CpGs comparing osteosarcomas and bone tissue samples, with high methylation heterogeneity, global hypomethylation and focal hypermethylation at CpG islands. Differentially methylated regions (DMR) were detected in 585 loci (319 hypomethylated and 266 hypermethylated), mapped to the promoter regions of 350 genes. These DMR genes were enriched for biological processes related to skeletal system morphogenesis, proliferation, inflammatory response, and signal transduction. Both methylation and expression data were validated in independent groups of cases. Six tumor suppressor genes harbored deletions or promoter hypermethylation (DLEC1, GJB2, HIC1, MIR149, PAX6, and WNT5A), and four oncogenes presented gains or hypomethylation (ASPSCR1, NOTCH4, PRDM16, and RUNX3). Our analysis also revealed hypomethylation at 6p22, a region that contains several histone genes. Copy-number changes in DNMT3B (gain) and TET1 (loss), as well as overexpression of DNMT3B in osteosarcomas provide a possible explanation for the observed phenotype of CpG island hypermethylation. While the detected open-sea hypomethylation likely contributes to the well-known osteosarcoma genomic instability, enriched CpG island hypermethylation suggests an underlying mechanism possibly driven by overexpression of DNMT3B likely resulting in silencing of tumor suppressors and DNA repair genes.

Skewed X-chromosome Inactivation in Women with Idiopathic Intellectual Disability is Indicative of Pathogenic Variants.

Intellectual disability (ID) is an early onset impairment in cognitive functioning and adaptive behavior, affecting approximately 1% of the population worldwide. Extreme skewing of X-chromosome inactivation (XCI) can be associated with ID phenotypes caused by pathogenic variants in the X chromosome. We analyzed the XCI pattern in blood samples of 194 women with idiopathic ID, using the androgen receptor gene (AR) methylation assay. Among the 136 patients who were informative, 11 (8%) presented with extreme or total XCI skewing (≥ 90%), which was significantly higher than expected by chance. Whole-exome data obtained from these 11 patients revealed the presence of dominant pathogenic variants in eight of them, all sporadic cases, resulting in a molecular diagnostic rate of 73% (8/11 patients). All variants were mapped to ID-related genes with dominant phenotypes: four variants in the X-linked genes DDX3X (an XCI escape gene; two cases), WDR45, and PDHA1, and four variants in the autosomal genes KCNB1, CTNNB1, YY1, and ANKRD11. Three of the autosomal genes had no obvious correlation with the observed XCI skewing. However, YY1 is a known transcriptional repressor that acts in the binding of the XIST long noncoding RNA on the inactive X chromosome, providing a mechanistic link between the pathogenic variant and the detected skewed XCI in the carrier. These data confirm that extreme XCI skewing in females with ID is highly indicative of causative X-linked pathogenic variants, and point to the possibility of identifying causative variants in autosomal genes with a XCI role.

Immunological evaluation of young unvaccinated patients with Turner syndrome after COVID-19.

OBJECTIVES: The X-chromosome contains the largest number of immune-related genes, which play a major role in COVID-19 symptomatology and susceptibility. Here, we had a unique opportunity to investigate, for the first time, COVID-19 outcomes in six unvaccinated young Brazilian patients with Turner syndrome (TS; 45, X0), including one case of critical illness in a child aged 10 years, to evaluate their immune response according to their genetic profile. METHODS: A serological analysis of humoral immune response against SARS-CoV-2, phenotypic characterization of antiviral responses in peripheral blood mononuclear cells after stimuli, and the production of cytotoxic cytokines of T lymphocytes and natural killer cells were performed in blood samples collected from the patients with TS during the convalescence period. Whole exome sequencing was also performed. RESULTS: Our volunteers with TS showed a delayed or insufficient humoral immune response to SARS-CoV-2 (particularly immunoglobulin G) and a decrease in interferon-γ production by cluster of differentiation (CD)4+ and CD8+ T lymphocytes after stimulation with toll-like receptors 7/8 agonists. In contrast, we observed a higher cytotoxic activity in the volunteers with TS than the volunteers without TS after phorbol myristate acetate/ionomycin stimulation, particularly granzyme B and perforin by CD8+ and natural killer cells. Interestingly, two volunteers with TS carry rare genetic variants in genes that regulate type I and III interferon immunity. CONCLUSION: Following previous reports in the literature for other conditions, our data showed that patients with TS may have an impaired immune response against SARS-CoV-2. Furthermore, other medical conditions associated with TS could make them more vulnerable to COVID-19.

SCAF4-related syndromic intellectual disability.

The causal link between variants in the SCAF4 gene and a syndromic form of intellectual disability (ID) was established in 2020 by Fliedner et al. Since then, no additional cases have been reported. We performed exome sequencing in a 16-year-old Brazilian male presenting with ID, epilepsy, behavioral problems, speech impairment, facial dysmorphisms, heart malformations, and obesity. A de novo pathogenic variant [SCAF4(NM_020706.2):c.374_375dup(p.Glu126LeufsTer20)] was identified. This is the second study reporting the involvement of SCAF4 in syndromic ID, and the description of the patient's clinical features contributes to defining the phenotypic spectrum of this recently described Mendelian disorder.

Challenges in selecting admixture models and marker sets to infer genetic ancestry in a Brazilian admixed population.

The inference of genetic ancestry plays an increasingly prominent role in clinical, population, and forensic genetics studies. Several genotyping strategies and analytical methodologies have been developed over the last few decades to assign individuals to specific biogeographic regions. However, despite these efforts, ancestry inference in populations with a recent history of admixture, such as those in Brazil, remains a challenge. In admixed populations, proportion and components of genetic ancestry vary on different levels: (i) between populations; (ii) between individuals of the same population, and (iii) throughout the individual's genome. The present study evaluated 1171 admixed Brazilian samples to compare the genetic ancestry inferred by tri-/tetra-hybrid admixture models and evaluated different marker sets from those with small numbers of ancestry informative markers panels (AIMs), to high-density SNPs (HDSNP) and whole-genome-sequence (WGS) data. Analyses revealed greater variation in the correlation coefficient of ancestry components within and between admixed populations, especially for minority ancestral components. We also observed positive correlation between the number of markers in the AIMs panel and HDSNP/WGS. Furthermore, the greater the number of markers, the more accurate the tri-/tetra-hybrid admixture models.

MUC22, HLA-A, and HLA-DOB variants and COVID-19 in resilient super-agers from Brazil.

Background: Although aging correlates with a worse prognosis for Covid-19, super elderly still unvaccinated individuals presenting mild or no symptoms have been reported worldwide. Most of the reported genetic variants responsible for increased disease susceptibility are associated with immune response, involving type I IFN immunity and modulation; HLA cluster genes; inflammasome activation; genes of interleukins; and chemokines receptors. On the other hand, little is known about the resistance mechanisms against SARS-CoV-2 infection. Here, we addressed polymorphisms in the MHC region associated with Covid-19 outcome in super elderly resilient patients as compared to younger patients with a severe outcome. Methods: SARS-CoV-2 infection was confirmed by RT-PCR test. Aiming to identify candidate genes associated with host resistance, we investigated 87 individuals older than 90 years who recovered from Covid-19 with mild symptoms or who remained asymptomatic following positive test for SARS-CoV-2 as compared to 55 individuals younger than 60 years who had a severe disease or died due to Covid-19, as well as to the general elderly population from the same city. Whole-exome sequencing and an in-depth analysis of the MHC region was performed. All samples were collected in early 2020 and before the local vaccination programs started. Results: We found that the resilient super elderly group displayed a higher frequency of some missense variants in the MUC22 gene (a member of the mucins' family) as one of the strongest signals in the MHC region as compared to the severe Covid-19 group and the general elderly control population. For example, the missense variant rs62399430 at MUC22 is two times more frequent among the resilient super elderly (p = 0.00002, OR = 2.24). Conclusion: Since the pro-inflammatory basal state in the elderly may enhance the susceptibility to severe Covid-19, we hypothesized that MUC22 might play an important protective role against severe Covid-19, by reducing overactive immune responses in the senior population.

The oldest unvaccinated Covid-19 survivors in South America.

BACKGROUND: Although older adults are at a high risk of severe or critical Covid-19, there are many cases of unvaccinated centenarians who had a silent infection or recovered from mild or moderate Covid-19. We studied three Brazilian supercentenarians, older than 110 years, who survived Covid-19 in 2020 before being vaccinated. RESULTS: Despite their advanced age, humoral immune response analysis showed that these individuals displayed robust levels of IgG and neutralizing antibodies (NAbs) against SARS-CoV-2. Enrichment of plasma proteins and metabolites related to innate immune response and host defense was also observed. None presented autoantibodies (auto-Abs) to type I interferon (IFN). Furthermore, these supercentenarians do not carry rare variants in genes underlying the known inborn errors of immunity, including particular inborn errors of type I IFN. CONCLUSION: These observations suggest that their Covid-19 resilience might be a combination of their genetic background and their innate and adaptive immunity.

Follow-up of young adult monozygotic twins after simultaneous critical coronavirus disease 2019: A case report.

Background: The influence of the host genome on coronavirus disease 2019 (COVID-19) susceptibility and severity is supported by reports on monozygotic (MZ) twins where both were infected simultaneously with similar disease outcomes, including several who died due to the SARS-CoV-2 infection within days apart. However, successive exposures to pathogens throughout life along with other environmental factors make the immune response unique for each individual, even among MZ twins. Case presentation and methods: Here we report a case of a young adult monozygotic twin pair, who caught attention since both presented simultaneously severe COVID-19 with the need for oxygen support despite age and good health conditions. One of the twins, who spent more time hospitalized, reported symptoms of long-COVID even 7 months after infection. Immune cell profile and specific responses to SARS-CoV-2 were evaluated as well as whole exome sequencing. Conclusion: Although the MZ twin brothers shared the same genetic mutations which may be associated with their increased risk of developing severe COVID-19, their clinical progression was different, reinforcing the role of both immune response and genetics in the COVID-19 presentation and course. Besides, post-COVID syndrome was observed in one of them, corroborating an association between the duration of hospitalization and the occurrence of long-COVID symptoms.

Global and local ancestry modulate APOE association with Alzheimer's neuropathology and cognitive outcomes in an admixed sample.

Dementia is more prevalent in Blacks than in Whites, likely due to a combination of environmental and biological factors. Paradoxically, clinical studies suggest an attenuation of APOE ε4 risk of dementia in African ancestry (AFR), but a dearth of neuropathological data preclude the interpretation of the biological factors underlying these findings, including the association between APOE ε4 risk and Alzheimer's disease (AD) pathology, the most frequent cause of dementia. We investigated the interaction between African ancestry, AD-related neuropathology, APOE genotype, and functional cognition in a postmortem sample of 400 individuals with a range of AD pathology severity and lack of comorbid neuropathology from a cohort of community-dwelling, admixed Brazilians. Increasing proportions of African ancestry (AFR) correlated with a lower burden of neuritic plaques (NP). However, for individuals with a severe burden of NP and neurofibrillary tangles (NFT), AFR proportion was associated with worse Clinical Dementia Rating sum of boxes (CDR-SOB). Among APOE ε4 carriers, the association between AFR proportion and CDR-SOB disappeared. APOE local ancestry inference of a subset of 309 individuals revealed that, in APOE ε4 noncarriers, non-European APOE background correlated with lower NP burden and, also, worse cognitive outcomes than European APOE when adjusting by NP burden. Finally, APOE ε4 was associated with worse AD neuropathological burden only in a European APOE background. APOE genotype and its association with AD neuropathology and clinical pattern are highly influenced by ancestry, with AFR associated with lower NP burden and attenuated APOE ε4 risk compared to European ancestry.

NAToRA, a relatedness-pruning method to minimize the loss of dataset size in genetic and omics analyses.

Genetic and omics analyses frequently require independent observations, which is not guaranteed in real datasets. When relatedness cannot be accounted for, solutions involve removing related individuals (or observations) and, consequently, a reduction of available data. We developed a network-based relatedness-pruning method that minimizes dataset reduction while removing unwanted relationships in a dataset. It uses node degree centrality metric to identify highly connected nodes (or individuals) and implements heuristics that approximate the minimal reduction of a dataset to allow its application to complex datasets. When compared with two other popular population genetics methodologies (PLINK and KING), NAToRA shows the best combination of removing all relatives while keeping the largest possible number of individuals in all datasets tested and also, with similar effects on the allele frequency spectrum and Principal Component Analysis than PLINK and KING. NAToRA is freely available, both as a standalone tool that can be easily incorporated as part of a pipeline, and as a graphical web tool that allows visualization of the relatedness networks. NAToRA also accepts a variety of relationship metrics as input, which facilitates its use. We also release a genealogies simulator software used for different tests performed in this study.

Genetic investigation of syndromic forms of obesity.

BACKGROUND: Syndromic obesity (SO) refers to obesity with additional phenotypes, including intellectual disability (ID)/developmental delay (DD), dysmorphic features, or organ-specific abnormalities. SO is rare, has high phenotypic variability, and frequently follows a monogenic pattern of inheritance. However, the genetic etiology of most cases of SO has not been elucidated. SUBJECTS AND METHODS: In this study, we investigated 20 SO patients by whole-exome sequencing (WES) trios to identify causal genetic variants. RESULTS: 4/20 patients had negative results for array comparative genomic hybridization (aCGH) analyses. In the remaining 15 patients, in addition to SNVs and indels, CNVs were also evaluated. Pathogenic/likely pathogenic (P/LP) SNVs/indels were detected in 6/20 patients (involving MED13L, AHDC1, EHMT1, MYT1L, GRIA3, and SETD1A), while two patients carried an inherited VUS. In addition, P/LP CNVs were observed in 3/15 patients (involving SATG2, KIAA0442, and MEIS2). CONCLUSIONS: All nine detected P/LP variants involved genes already known to lead to syndromic ID/DD; however, for only two genes (EHMT1 and MYT1L) is the link with obesity well established. This is the first study applying a comprehensive genomic investigation of an SO cohort, showing a high diagnostic yield (~47%). Additionally, our findings suggested that several known ID/DD genes may also predispose individuals to SO.

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.

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.