Genome-Wide DNA Methylation Analysis of a Cohort of 41 Patients Affected by Oculo-Auriculo-Vertebral Spectrum (OAVS)

Can some metabolic one-carbon cycle linked diseases be prevented? The impact of treating hypo-fertile couples carrying MTHFR SNPs with folic acid and 5-MTHF on outcomes in the offspring: a case retrospective series

PURPOSE

In our practice, testing hypo-fertile patients for circulating homocysteine (Hcy) and the two principal MTHFR SNPs (677C > T and 1298A > C) has been routine for the past 7 years. Couples carrying a genetic background known to be associated with the disease were proposed treatment regimens consisting of 5-methyl tetrahydrofolate (5-MTHF) together with nutritional support of the one-carbon cycle (1-CC). Some patients preferred to continue with folic acid (FA) as prescribed by their referring gynecologist/obstetrician: this gave us the opportunity to compare outcomes between the two groups of patients.

METHODS

After successful live birth deliveries, we compared health characteristics and circulating Hcy in the offspring from ages 2 to 6 years, i.e., after cessation of breastfeeding and before puberty. Follow-up included children of 21 couples who were treated with FA vs 36 couples treated with 5-MTHF.

RESULTS

In the FA-treated group, we found two children with autism spectrum disorder (ASD) syndrome, one child with significantly elevated circulating Hcy (19 µM at the age of 2 years), and one child affected by oculo-auriculo-vertebral spectrum (OAVS), a syndrome known to be linked to DNA methylation. No pathology of any kind was detected in children of the 5-MTHF treatment group.

CONCLUSION

Treatment with 5-MTHF is safe and effective for both males and females. It should be implemented in order to avoid disruption of methylation linked to folate metabolism during oocyte maturation and pregnancy, and subsequently in the offspring. This type of treatment should be considered to avoid metabolic diseases linked to elevated homocysteine.

Epigenetic patterns, accelerated biological aging, and enhanced epigenetic drift detected 6 months following COVID-19 infection: insights from a genome-wide DNA methylation study

BACKGROUND

The epigenetic status of patients 6-month post-COVID-19 infection remains largely unexplored. The existence of long-COVID, or post-acute sequelae of SARS-CoV-2 infection (PASC), suggests potential long-term changes. Long-COVID includes symptoms like fatigue, neurological issues, and organ-related problems, regardless of initial infection severity. The mechanisms behind long-COVID are unclear, but virus-induced epigenetic changes could play a role.

METHODS AND RESULTS

Our study explores the lasting epigenetic impacts of SARS-CoV-2 infection. We analyzed genome-wide DNA methylation patterns in an Italian cohort of 96 patients 6 months after COVID-19 exposure, comparing them to 191 healthy controls. We identified 42 CpG sites with significant methylation differences (FDR < 0.05), primarily within CpG islands and gene promoters. Dysregulated genes highlighted potential links to glutamate/glutamine metabolism, which may be relevant to PASC symptoms. Key genes with potential significance to COVID-19 infection and long-term effects include GLUD1, ATP1A3, and ARRB2. Furthermore, Horvath's epigenetic clock showed a slight but significant age acceleration in post-COVID-19 patients. We also observed a substantial increase in stochastic epigenetic mutations (SEMs) in the post-COVID-19 group, implying potential epigenetic drift. SEM analysis identified 790 affected genes, indicating dysregulation in pathways related to insulin resistance, VEGF signaling, apoptosis, hypoxia response, T-cell activation, and endothelin signaling.

CONCLUSIONS

Our study provides valuable insights into the epigenetic consequences of COVID-19. Results suggest possible associations with accelerated aging, epigenetic drift, and the disruption of critical biological pathways linked to insulin resistance, immune response, and vascular health. Understanding these epigenetic changes could be crucial for elucidating the complex mechanisms behind long-COVID and developing targeted therapeutic interventions.

Whole genome sequencing of a family with autosomal dominant features within the oculoauriculovertebral spectrum

Background

Oculoauriculovertebral Spectrum (OAVS) encompasses abnormalities on derivatives from the first and second pharyngeal arches including macrostomia, hemifacial microsomia, micrognathia, preauricular tags, ocular and vertebral anomalies. We present genetic findings on a three-generation family affected with macrostomia, preauricular tags and uni- or bilateral ptosis following an autosomal dominant pattern.

Methods

We generated whole genome sequencing data for the proband, affected parent and unaffected paternal grandparent followed by Sanger sequencing on 23 family members for the top 10 candidate genes: KCND2, PDGFRA, CASP9, NCOA3, WNT10A, SIX1, MTF1, KDR/VEGFR2, LRRK1, and TRIM2 We performed parent and sibling-based transmission disequilibrium tests and burden analysis via a penalized linear mixed model, for segregation and mutation burden respectively. Next, via bioinformatic tools we predicted protein function, mutation pathogenicity and pathway enrichment to investigate the biological relevance of mutations identified.

Results

Rare missense mutations in SIX1, KDR/VEGFR2, and PDGFRA showed the best segregation with the OAV phenotypes in this family. When considering any of the 3 OAVS phenotypes as an outcome, SIX1 had the strongest associations in parent-TDTs and sib-TDTs (p=0.025, p=0.052) (unadjusted p-values). Burden analysis identified SIX1 (RC=0.87) and PDGFRA (RC=0.98) strongly associated with OAVS severity. Using phenotype-specific outcomes, sib-TDTs identified SIX1 with uni- or bilateral ptosis (p=0.049) and ear tags (p=0.01), and PDGFRA and KDR/VEGFR2 with ear tags (both p<0.01).

Conclusion

SIX1, PDGFRA, and KDR/VEGFR2 are strongly associated to OAVS phenotypes. SIX1 has been previously associated with OAVS ear malformations and is co-expressed with EYA1 during ear development. Efforts to strengthen the genotype-phenotype co-relation underlying the OAVS are key to discover etiology, family counseling and prevention.

Identification of potential molecular mechanism related to craniofacial dysmorphism caused by FOXI3 deficiency

BACKGROUND

Hemifacial macrosomia (HFM, OMIM 164210) is a complex and highly heterogeneous disease. FORKHEAD BOX I3 (FOXI3) is a susceptibility gene for HFM, and mice with loss of function of Foxi3 did exhibit a phenotype similar to craniofacial dysmorphism. However, the specific pathogenesis of HFM caused by FOXI3 deficiency remains unclear till now.

METHOD

In this study, we first constructed a Foxi3 deficiency (Foxi3-/- ) mouse model to verify the craniofacial phenotype of Foxi3-/- mice, and then used RNAseq data for gene differential expression analysis to screen candidate pathogenic genes, and conducted gene expression verification analysis using quantitative real-time PCR.

RESULTS

By observing the phenotype of Foxi3-/- mice, we found that craniofacial dysmorphism was present. The results of comprehensive bioinformatics analysis suggested that the craniofacial dysmorphism caused by Foxi3 deficiency may be involved in the PI3K-Akt signaling pathway. Quantitative real-time PCR results showed that the expression of PI3K-Akt signaling pathway-related gene Akt2 was significantly increased in Foxi3-/- mice.

CONCLUSION

The craniofacial dysmorphism caused by the deficiency of Foxi3 may be related to the expression of Akt2 and PI3K-Akt signaling pathway. This study laid a foundation for understanding the function of FOXI3 and the pathogenesis and treatment of related craniofacial dysmorphism caused by FOXI3 dysfunction.

Genetics and Epigenetics in the Genesis and Development of Microtia

Microtia is a congenital malformation of the external and middle ear associated with varying degrees of severity that range from mild structural abnormalities to the absence of the external ear and auditory canal. Globally, it is the second most common congenital craniofacial malformation and is typically caused by inherited defects, external factors, or the interaction between genes and external factors. Epigenetics notably represents a bridge between genetics and the environment. This review has devoted attention to the current proceedings of the genetics and epigenetics of microtia and related syndromes.

Oculoauriculofrontonasal syndrome: Refining the phenotype through a new case series and literature review

The oculoauriculofrontonasal syndrome (OAFNS) is a rare condition, with unknown etiology, characterized by the association of frontonasal dysplasia (FND) and oculoauriculovertebral spectrum (OAVS). Main clinical findings include widely spaced eyes, epibulbar dermoid, broad nose, mandibular hypoplasia, and preauricular tags. Here, we describe a case series of 32 Brazilian individuals with OAFNS and review the literature ascertaining individuals presenting phenotypes compatible with the diagnosis of OAFNS, aiming to refine the phenotype. This series emphasizes the phenotypic variability of the OAFNS and highlights the occurrence of rare craniofacial clefts as a part of the phenotype. The ectopic nasal bone, a hallmark of OAFNS, was frequent in our series, reinforcing the clinical diagnosis. The absence of recurrence, consanguinity, chromosomal, and genetic abnormalities reinforces the hypothesis of a nontraditional inheritance model. The phenotypic refinement provided by this series contributes to an investigation regarding the etiology of OAFNS.

Exploration of Novel Genetic Evidence and Clinical Significance Into Hemifacial Microsomia Pathogenesis

The authors browsed through past genetic findings in hemifacial microsomia along with our previously identified mutations in ITGB4 and PDE4DIP from whole genome sequencing of hemifacial microsomia patients. Wondering whether these genes influence mandibular bone modeling by regulation on osteogenesis, the authors approached mechanisms of hemifacial microsomia through this investigation into gene knockdown effects in vitro. MC3T3E1 cells were divided into 5 groups: the negative control group without osteogenesis induction or siRNA, the positive control group with only osteogenesis induction, and 3 gene silenced groups with both osteogenesis induction and siRNA. Validation of transfection was through fluorescence microscopy and quantitative real-time Polymerase chain reaction on knockdown efficiency. Changes in expression levels of the 3 genes during osteogenesis and impact of Itgb4 and Pde4dip knockdown on osteogenesis were examined by quantitative real-time Polymerase chain reaction, alkaline phosphatase, and alizarin red staining. Elevation of osteogenic genes Alpl, Col1a1, Bglap, Spp1, and Runx2 verified successful osteogenesis. Both genes were upregulated under osteogenic induction, while they had different trends over time. Intracellular fluorophores under microscope validated successful transfection and si-m-Itgb4_003, si-m-Pde4dip_002 had satisfactory knockdown effects. During osteogenesis, Pde4dip knockdown enhanced Spp1 expression (1.95±0.13 folds, P =0.045). The authors speculated that these genes may have different involvements in osteogenesis. Stimulated expression of Spp1 by Pde4dip knockdown may suggest that Pde4dip inhibits osteogenesis.

Role of epigenetics in the clinical evolution of COVID-19 disease. Epigenome-wide association study identifies markers of severe outcome

BACKGROUND

COVID-19 has a wide spectrum of clinical manifestations and given its impact on morbidity and mortality, there is an unmet medical need to discover endogenous cellular and molecular biomarkers that predict the expected clinical course of the disease. Recently, epigenetics and especially DNA methylation have been pointed out as a promising tool for outcome prediction in several diseases.

METHODS AND RESULTS

Using the Illumina Infinium Methylation EPIC BeadChip850K, we investigated genome-wide differences in DNA methylation in an Italian Cohort of patients with comorbidities and compared severe (n = 64) and mild (123) prognosis. Results showed that the epigenetic signature, already present at the time of Hospital admission, can significantly predict risk of severe outcomes. Further analyses provided evidence of an association between age acceleration and a severe prognosis after COVID-19 infection. The burden of Stochastic Epigenetic Mutation (SEMs) has been significantly increased in patients with poor prognosis. Results have been replicated in silico considering COVID-19 negative subjects and available previously published datasets.

CONCLUSIONS

Using original methylation data and taking advantage of already published datasets, we confirmed in the blood that epigenetics is actively involved in immune response after COVID-19 infection, allowing the identification of a specific signature able to discriminate the disease evolution. Furthermore, the study showed that epigenetic drift and age acceleration are associated with severe prognosis. All these findings prove that host epigenetics undergoes notable and specific rearrangements to respond to COVID-19 infection which can be used for a personalized, timely, and targeted management of COVID-19 patients during the first stages of hospitalization.

Stochastic epigenetic mutations as possible explanation for phenotypical discordance among twins with congenital hypothyroidism

PURPOSE

The elevated frequency of discordance for congenital hypothyroidism (CH) phenotype between monozygotic twins suggests the involvement of non-mendelian mechanisms. The aim of the study was to investigate the role of epigenetics in CH pathogenesis.

METHODS

A genome-wide DNA methylation analysis was performed on the peripheral blood of 23 twin pairs (10 monozygotic and 13 dizygotic), 4 concordant and 19 discordant pairs for CH at birth.

RESULTS

Differential methylation analysis did not show significant differences in methylation levels between CH cases and controls, but a different methylation status of several genes may explain the CH discordance of a monozygotic twin couple carrying a monoallelic nonsense mutation of DUOX2. In addition, the median number of hypo-methylated Stochastic Epigenetic Mutations (SEMs) resulted significantly increased in cases compared to controls. The prioritization analysis for CH performed on the genes epimutated exclusively in the cases identified SLC26A4, FOXI1, NKX2-5 and TSHB as the genes with the highest score. The analysis of significantly SEMs-enriched regions led to the identification of two genes (FAM50B and MEG8) that resulted epigenetically dysregulated in cases.

CONCLUSION

Epigenetic modifications may potentially account for CH pathogenesis and explain discordance among monozygotic twins.

Epigenetics of Autism Spectrum Disorders: A Multi-level Analysis Combining Epi-signature, Age Acceleration, Epigenetic Drift and Rare Epivariations Using Public Datasets

BACKGROUND

Epigenetics of Autism Spectrum Disorders (ASD) is still an understudied field. The majority of the studies on the topic used an approach based on mere classification of cases and controls.

OBJECTIVE

The present study aimed at providing a multi-level approach in which different types of epigenetic analysis (epigenetic drift, age acceleration) are combined.

METHODS

We used publicly available datasets from blood (n = 3) and brain tissues (n = 3), separately. Firstly, we evaluated for each dataset and meta-analyzed the differential methylation profile between cases and controls. Secondly, we analyzed age acceleration, epigenetic drift and rare epigenetic variations.

RESULTS

We observed a significant epi-signature of ASD in blood but not in brain specimens. We did not observe significant age acceleration in ASD, while epigenetic drift was significantly higher compared to controls. We reported the presence of significant rare epigenetic variations in 41 genes, 35 of which were never associated with ASD. Almost all genes were involved in pathways linked to ASD etiopathogenesis (i.e., neuronal development, mitochondrial metabolism, lipid biosynthesis and antigen presentation).

CONCLUSION

Our data support the hypothesis of the use of blood epi-signature as a potential tool for diagnosis and prognosis of ASD. The presence of an enhanced epigenetic drift, especially in brain, which is linked to cellular replication, may suggest that alteration in epigenetics may occur at a very early developmental stage (i.e., fetal) when neuronal replication is still high.

Mutations in the TBX15-ADAMTS2 pathway associate with a novel soft palate dysplasia

We reported de novo variants in specific exons of the TBX15 and ADAMTS2 genes in a hitherto undescribed class of patients with unique craniofacial developmental defects. The nine unrelated patients represent unilateral soft palate hypoplasia, lost part of the sphenoid bone in the pterygoid process, but the uvula developed completely. Interestingly, these clinical features are contrary to the palate's anterior-posterior (A-P) developmental direction. Based on developmental characteristics, we suggested that these cases correspond to a novel craniofacial birth defect different from cleft palate, and we named it soft palate dysplasia (SPD). However, little is known about the molecular mechanism of the ADAMTS2 and TBX15 genes in the regulation of soft palate development. Phylogenetic analysis showed that the sequences around these de novo mutation sites are conserved between species. Through cellular co-transfections and chromatin immunoprecipitation assays, we demonstrate that TBX15 binds to the promoter regions of the ADAMTS2 gene and activates the promoter activity. Furthermore, we show that TBX15 and ADAMTS2 are colocalization in the posterior palatal mesenchymal cells during soft palate development in E13.5 mice embryos. Based on these data, we propose that the disruption of the TBX15-ADAMTS2 signaling pathway during embryogenesis leads to a novel SPD.

OTX2 duplications: a recurrent cause of oculo-auriculo-vertebral spectrum

Background

Oculo-auriculo-vertebral spectrum (OAVS) is the second most common cause of head and neck malformations in children after orofacial clefts. OAVS is clinically heterogeneous and characterised by a broad range of clinical features including ear anomalies with or without hearing loss, hemifacial microsomia, orofacial clefts, ocular defects and vertebral abnormalities. Various genetic causes were associated with OAVS and copy number variations represent a recurrent cause of OAVS, but the responsible gene often remains elusive.

Methods

We described an international cohort of 17 patients, including 10 probands and 7 affected relatives, presenting with OAVS and carrying a 14q22.3 microduplication detected using chromosomal microarray analysis. For each patient, clinical data were collected using a detailed questionnaire addressed to the referring clinicians. We subsequently studied the effects ofOTX2overexpression in a zebrafish model.

Results

We defined a 272 kb minimal common region that only overlaps with theOTX2gene. Head and face defects with a predominance of ear malformations were present in 100% of patients. The variability in expressivity was significant, ranging from simple chondromas to severe microtia, even between intrafamilial cases. Heterologous overexpression ofOTX2in zebrafish embryos showed significant effects on early development with alterations in craniofacial development.

Conclusions

Our results indicate that properOTX2dosage seems to be critical for the normal development of the first and second branchial arches. Overall, we demonstrated thatOTX2genomic duplications are a recurrent cause of OAVS marked by auricular malformations of variable severity.

Adhesion molecule Amigo2 is involved in the fasciculation process of the fasciculus retroflexus

BACKGROUND

The fasciculus retroflexus is the prominent efferent pathway from the habenular complex. Medial habenular axons form a core packet whereas lateral habenular axons course in a surrounding shell. Both groups of fibers share the same initial pathway but differ in the final segment of the tract, supposedly regulated by surface molecules. The gene Amigo2 codes for a membrane adhesion molecule with an immunoglobulin-like domain 2 and is selectively expressed in the medial habenula. We present it as a candidate for controlling the fasciculation behavior of medial habenula axons.

RESULTS

First, we studied the development of the habenular efferents in an Amigo2 lack of function mouse model. The fasciculus retroflexus showed a variable defasciculation phenotype. Gain of function experiments allowed us to generate a more condensed tract and rescued the Amigo2 knock-out phenotype. Changes in Amigo2 function did not alter the course of habenular fibers.

CONCLUSION

We have demonstrated that Amigo2 plays a subtle role in the fasciculation of the fasciculus retroflexus.

The Enigmatic Etiology of Oculo-Auriculo-Vertebral Spectrum (OAVS): An Exploratory Gene Variant Interaction Approach in Candidate Genes

The clinical diagnosis of oculo-auriculo-vertebral spectrum (OAVS) is established when microtia is present in association with hemifacial hypoplasia (HH) and/or ocular, vertebral, and/or renal malformations. Genetic and non-genetic factors have been associated with microtia/OAVS. Although the etiology remains unknown in most patients, some cases may have an autosomal dominant, autosomal recessive, or multifactorial inheritance. Among the possible genetic factors, gene−gene interactions may play important roles in the etiology of complex diseases, but the literature lacks related reports in OAVS patients. Therefore, we performed a gene−variant interaction analysis within five microtia/OAVS candidate genes (HOXA2, TCOF1, SALL1, EYA1 and TBX1) in 49 unrelated OAVS Mexican patients (25 familial and 24 sporadic cases). A statistically significant intergenic interaction (p-value < 0.001) was identified between variants p.(Pro1099Arg) TCOF1 (rs1136103) and p.(Leu858=) SALL1 (rs1965024). This intergenic interaction may suggest that the products of these genes could participate in pathways related to craniofacial alterations, such as the retinoic acid (RA) pathway. The absence of clearly pathogenic variants in any of the analyzed genes does not support a monogenic etiology for microtia/OAVS involving these genes in our patients. Our findings could suggest that in addition to high-throughput genomic approaches, future gene−gene interaction analyses could contribute to improving our understanding of the etiology of microtia/OAVS.

Oculo-auriculo-vertebral spectrum: new genes and literature review on a complex disease

Oculo-auriculo-vertebral spectrum (OAVS) or Goldenhar syndrome is due to an abnormal development of first and second branchial arches derivatives during embryogenesis and is characterised by hemifacial microsomia associated with auricular, ocular and vertebral malformations. The clinical and genetic heterogeneity of this spectrum with incomplete penetrance and variable expressivity, render its molecular diagnosis difficult. Only a few recurrent CNVs and genes have been identified as causatives in this complex disorder so far. Prenatal environmental causal factors have also been hypothesised. However, most of the patients remain without aetiology. In this review, we aim at updating clinical diagnostic criteria and describing genetic and non-genetic aetiologies, animal models as well as novel diagnostic tools and surgical management, in order to help and improve clinical care and genetic counselling of these patients and their families.

Copy number variation analysis implicates novel pathways in patients with oculo-auriculo-vertebral-spectrum and congenital heart defects

Oculo-auriculo-vertebral spectrum (OAVS) is a developmental disorder of craniofacial morphogenesis. Its etiology is unclear, but assumed to be complex and heterogeneous, with contribution of both genetic and environmental factors. We assessed the occurrence of copy number variants (CNVs) in a cohort of 19 unrelated OAVS individuals with congenital heart defect. Chromosomal microarray analysis identified pathogenic CNVs in 2/19 (10.5%) individuals, and CNVs classified as variants of uncertain significance in 7/19 (36.9%) individuals. Remarkably, two subjects had small intragenic CNVs involving DACH1 and DACH2, two paralogs coding for key components of the PAX-SIX-EYA-DACH network, a transcriptional regulatory pathway controlling developmental processes relevant to OAVS and causally associated with syndromes characterized by craniofacial involvement. Moreover, a third patient showed a large duplication encompassing DMBX1/OTX3, encoding a transcriptional repressor of OTX2, another transcription factor functionally connected to the DACH-EYA-PAX network. Among the other relevant CNVs, a deletion encompassing HSD17B6, a gene connected with the retinoic acid signaling pathway, whose dysregulation has been implicated in craniofacial malformations, was also identified. Our findings suggest that CNVs affecting gene dosage likely contribute to the genetic heterogeneity of OAVS, and implicate the PAX-SIX-EYA-DACH network as novel pathway involved in the etiology of this developmental trait.