BACKGROUND: During gestation, stressors to the fetus, including viral exposure or maternal psychological distress, can fundamentally alter the neonatal epigenome, and may be associated with long-term impaired developmental outcomes. The impact of in utero exposure to the COVID-19 pandemic on the newborn epigenome has yet to be described. METHODS: This study aimed to determine whether there are unique epigenetic signatures in newborns who experienced otherwise healthy pregnancies that occurred during the COVID-19 pandemic (Project RESCUE). The pre-pandemic control and pandemic cohorts (Project RESCUE) included in this study are part of a prospective observational and longitudinal cohort study that evaluates the impact of elevated prenatal maternal stress during the COVID-19 pandemic on early childhood neurodevelopment. Using buccal swabs collected at birth, differential DNA methylation analysis was performed using the Infinium MethylationEPIC arrays and linear regression analysis. Pathway analysis and gene ontology enrichment were performed on resultant gene lists. RESULTS: Widespread differential methylation was found between neonates exposed in utero to the pandemic and pre-pandemic neonates. In contrast, there were no apparent epigenetic differences associated with maternal COVID-19 infection during pregnancy. Differential methylation was observed among genomic sites that underpin important neurological pathways that have been previously reported in the literature to be differentially methylated because of prenatal stress, such as NR3C1. CONCLUSIONS: The present study reveals potential associations between exposure to the COVID-19 pandemic during pregnancy and subsequent changes in the newborn epigenome. While this finding warrants further investigation, it is a point that should be considered in any study assessing newborn DNA methylation studies obtained during this period, even in otherwise healthy pregnancies.

BACKGROUND:Disorders of sex development (DSD) have an estimated frequency of 0.5% of live births encompassing a variety of urogenital anomalies ranging from mild hypospadias to a discrepancy between sex chromosomes and external genitalia. In order to identify the underlying genetic etiology, we had performed exome sequencing in a subset of DSD cases with 46,XY karyotype and were able to identify the causative genetic variant in 35% of cases. While the genetic etiology was not ascertained in more than half of the cases, a large number of variants of unknown clinical significance (VUS) were identified in those exomes. METHODS:To investigate the relevance of these VUS in regards to the patient's phenotype, we utilized a mouse model in which the presence of a Y chromosome from the poschiavinus strain (Y POS ) on a C57BL/6J (B6) background results in XY undervirilization and sex reversal, a phenotype characteristic to a large subset of human 46,XY DSD cases. We assessed gene expression differences between B6-Y B6 and undervirilized B6-Y POS gonads at E11.5 and identified 515 differentially expressed genes (308 underexpressed and 207 overexpressed in B6-Y POS males). RESULTS:We identified 15 novel candidate genes potentially involved in 46,XY DSD pathogenesis by filtering the list of human VUS-carrying genes provided by exome sequencing with the list of differentially expressed genes from B6-Y POS mouse model. Additionally, we identified that 7 of the 15 candidate genes were significantly underexpressed in the XY gonads of mice with suppressed Sox9 expression in Sertoli cells suggesting that some of the candidate genes may be downstream of a well-known sex determining gene, Sox9. CONCLUSION:The use of a DSD-specific animal model improves variant interpretation by correlating human sequence variants with transcriptome variation.

Background

Massively parallel DNA sequencing, such as exome sequencing, has become a routine clinical procedure to identify pathogenic variants responsible for a patient's phenotype. Exome sequencing has the capability of reliably identifying inherited and de novo single-nucleotide variants, small insertions, and deletions. However, due to the use of 100-300-bp fragment reads, this platform is not well powered to sensitively identify moderate to large structural variants (SV), such as insertions, deletions, inversions, and translocations.

Methods

To overcome these limitations, we used next-generation mapping (NGM) to image high molecular weight double-stranded DNA molecules (megabase size) with fluorescent tags in nanochannel arrays for de novo genome assembly. We investigated the capacity of this NGM platform to identify pathogenic SV in a series of patients diagnosed with Duchenne muscular dystrophy (DMD), due to large deletions, insertion, and inversion involving the DMD gene.

Results

We identified deletion, duplication, and inversion breakpoints within DMD. The sizes of deletions were in the range of 45-250 Kbp, whereas the one identified insertion was approximately 13 Kbp in size. This method refined the location of the break points within introns for cases with deletions compared to current polymerase chain reaction (PCR)-based clinical techniques. Heterozygous SV were detected in the known carrier mothers of the DMD patients, demonstrating the ability of the method to ascertain carrier status for large SV. The method was also able to identify a 5.1-Mbp inversion involving the DMD gene, previously identified by RNA sequencing.

Conclusions

We showed the ability of NGM technology to detect pathogenic structural variants otherwise missed by PCR-based techniques or chromosomal microarrays. NGM is poised to become a new tool in the clinical genetic diagnostic strategy and research due to its ability to sensitively identify large genomic variations.

We report on the development of a methylation analysis workflow for optical detection of fluorescent methylation profiles along chromosomal DNA molecules. In combination with Bionano Genomics genome mapping technology, these profiles provide a hybrid genetic/epigenetic genome-wide map composed of DNA molecules spanning hundreds of kilobase pairs. The method provides kilobase pair-scale genomic methylation patterns comparable to whole-genome bisulfite sequencing (WGBS) along genes and regulatory elements. These long single-molecule reads allow for methylation variation calling and analysis of large structural aberrations such as pathogenic macrosatellite arrays not accessible to single-cell second-generation sequencing. The method is applied here to study facioscapulohumeral muscular dystrophy (FSHD), simultaneously recording the haplotype, copy number, and methylation status of the disease-associated, highly repetitive locus on Chromosome 4q.