Prenatal phenotype of PNKP-related primary microcephaly associated with variants in the FHA and Phosphatase domain

Biallelic PNKP variants cause heterogeneous disorders ranging from neurodevelopmental disorder with microcephaly/seizures to adult-onset Charcot-Marie-Tooth disease. To date, only postnatal descriptions exist. We present the first prenatal diagnosis of PNKP-related primary microcephaly. Detailed pathological examination of a male fetus revealed micrencephaly with extracerebral malformations and thus presumed syndromic microcephaly. A recessive disorder was suspected because of previous pregnancy termination for similar abnormalities in a sibling fetus. Prenatal trio exome sequencing identified compound-heterozygosity for the PNKP variants c.498G>A, p.[(=),0?] and c.302C>T, p.(Pro101Leu). Segregation confirmed both variants in the sibling fetus. Through RNA analyses, we characterized skipping of exon 4 affecting the PNKP Forkhead-associated (FHA) and Phosphatase domains (p.Leu67_Lys166del) as the predominant effect of the c.498G>A variant. We retrospectively investigated two unrelated individuals diagnosed with biallelic PNKP-variants to compare prenatal/postnatal phenotypes. Both carry the same splice-donor variant c.1029+2T>C in trans with a variant in the FHA domain (c.311T>C, p.(Leu104Pro) and c.151G>C, p.(Val51Leu), respectively). RNA-seq showed complex splicing events for c.1029+2T>C and c.151G>C. Computational modelling and structural analysis revealed significant clustering of missense variants in the FHA domain, with some variants potentially generating structural damage. Our detailed clinical description extends the PNKP-continuum to the prenatal stage. Investigating possible PNKP-variant effects using RNA and structural modelling, we highlight the mutational complexity and exemplify a framework for variant characterization in this multi-domain protein.


INTRODUCTION
Microcephaly, or rather micrencephaly (abnormally small brain) in the narrow sense, is defined as an occipital frontal circumference (OFC) below -2 SD of the mean for (gestational) age and sex and can occur in isolated form or in a syndromic context (1). If detected prenatally, it is classified as primary microcephaly (PM) in contrast to secondary microcephaly (SM) developing after birth. Infections, traumata, ischemic events, exposure to teratogens and genetic disorders are possible etiologies (1,2). As head growth -among many other factors -depends on normal neuronal tissue proliferation, requiring continuous cell division, several genetic neurodevelopmental and neurodegenerative disorders in this context are caused by variants affecting DNA repair genes, highlighting the importance of the pathways in neurogenesis (3).
Most pathogenic PNKP variants described so far are either truncating or located in the Cterminal Kinase domain (11). While genotype-phenotype correlations have been attempted and C-terminal variants have been implied to cause the milder adult-onset diseases, no clear relation could yet be established. Instead, it has even been postulated that the pathogenic variants observed present with rather mild mutational effects, due to survivorship bias, and more damaging variants would result in non-viability (11).
Here, we describe the first prenatal identification of biallelic PNKP variants affecting the region between the N-terminal Forkhead-associated (FHA) domain and the Phosphatase domain causative for severe early onset of PM. We provide detailed descriptions based on . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted September 28, 2021. ; https://doi.org/10.1101/2021.09.25.21261035 doi: medRxiv preprint prenatal imaging and syndrome-oriented fetal autopsies of two affected sibling fetuses, we compare the fetal phenotype with two individuals with PNKP-associated disorder diagnosed previously and give a review of similar cases in the literature. Additionally, we performed extensive RNA analyses to characterize aberrant splicing of identified variants and used structural modelling to investigate the effects of missense variants.

Ethics approval
The study adheres to the principles set out in the Declaration of Helsinki. The Ethical Committee of the Medical Faculty, Leipzig University approved genetic testing in a research setting for all probands within the study. Written consent of the parents to publish genetic and clinical data, as well as prenatal images and postmortem photographs (P1, P2), sonography (P1, P2), and magnetic resonance imaging (MRI) images (P1, P3 and P4) was received and archived by the authors.

Genetic analyses and review of PNKP variants
P2 and both his parents underwent trio-exome sequencing. In P3 and P4, clinical exome sequencing (CES) was performed. Segregation was confirmed through Sanger sequencing in all. Technical details and primer sequences are provided in Supplementary notes. All PNKP variants have been submitted to ClinVar (File S3 (12) sheet "PNKP_variants").

Clinical data collection
We used a questionnaire for retrospective phenotype analysis in which clinical terms were standardized using the Human Phenotype Ontology (HPO) terminology (15) based on a review of previously reported clinical associations in PNKP-disorders (6)(7)(8)(9)(10). The sheet . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted September 28, 2021. ; https://doi.org/10.1101/2021.09.25.21261035 doi: medRxiv preprint was sent for evaluation to the pediatric neurologist or pathologist, respectively and data from available clinical reports were added. Pre-and postnatal anthropometric measurements were compared to standard values according to Potter et Craig (16), or WHO child growth standards (17). Comprehensive results are shown in File S2 (12) sheet "clinical" and Figure S1.

Fetal autopsy and RNA extraction from fetal tissue
Fetal pathological examination of P1 and P2 was performed as previously described (18), detailed information is shown in File S1. Cryopreserved native skeletal muscle tissue of P2 was processed with QIAshredder (Qiagen, Hilden, Germany) and RNA was extracted according to the manufacturer's protocol (RNeasy Mini, Qiagen, Hilden, Germany).

RNA analyses
In family 1 we performed RT-PCR as described previously (19) using PAXgene RNA (Becton Dickinson, Franklin Lakes, NJ) in the parents and fetal skeletal muscle RNA derived cDNA. In family 3, we performed RNA-seq from PAXgene RNA using the TruSeq RNA Library Prep Kit v2 according to the manufacturer's instructions and pairedend sequencing on an Illumina NextSeq platform. Bioinformatic workup included an established pipeline from our institute. In brief, reads were demultiplexed, adapters trimmed, overrepresented sequences removed before we aligned the remaining reads to the hg38 reference using STAR aligner (20). We visualized and inspected the alignments for aberrant splicing as described previously (21) and applied iREAD (22) to quantify the visually observed intron retention events. Plotting of Figure 3 was performed with R using the packages "ggplot2", "Gviz" and "trackViewer" and Inkscape was used to adjust figure components. Detailed procedures including primer sequences, reagents and software versions are provided in the Supplementary notes.

Analysis of missense variant spectrum
Analysis of disease-associated missense variants in the linear protein representation, clustering analysis in 3D and structural modelling of missense variants using the crystal structure 2BRF (23) was performed as described previously (11,13,19) and is detailed in the Supplementary notes. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint

Prenatal phenotype in two sibling fetuses
In the pregnancy of a healthy non-consanguineous couple, routine sonography revealed microcephaly, abnormal skull shape and microretrognathia in the male fetus (P2). Followup ultrasound examinations displayed progression of the described anomalies.
Amniocentesis for genetic testing (trio ES) was performed. Severe fetal anomalies and supposing genetic background caused the parents to decide for termination of pregnancy.
A previous pregnancy had been terminated after prenatal imaging had confirmed multiple anomalies in a female fetus (P1). Ultrasound had shown microcephaly, asymmetric skull shape, abnormal brain development, cerebellar hypoplasia, cataract of both eyes, and facial abnormalities. Prenatal MRI confirmed microcephaly, large supratentorial defects of brain parenchyma in occipital, parietal and frontal regions, severe cerebellar hypoplasia, dilatation and fusion of both lateral ventricles ( Figure 1C). Corpus callosum and septum pellucidum were not determinable. Brainstem and thalamus appeared normal, basal ganglia were not accessible. Bulbi of the eyes differed in size and signal. Other organs were described as normal.
Both fetuses were examined after TOP. The female fetus P1 measured 24.0 cm crown to heel length (standard value: 26.2±3.6 cm (16) having a weight of 301 g (standard value: 353±125 g (16)). OFC was 15.0 cm (-4.3 SD; 5 th percentile at GA: 17.3 cm (24)). Autopsy presumed severe cerebral parenchymatous defects, profound hypoplasia of posterior cranial fossa and confluence of the side ventricles. The findings were initially interpreted as arhinencephaly/holoprosencephaly. Further, extracerebral anomalies were not assessable, also because of autolysis. Placenta appeared hypotrophic with signs of insufficiency.
The sibling fetus P2 was almost age-appropriate in terms of crown to heel length of 18.8 cm and weight (18.8 cm resp. 124.3 g) (standard values: 20 cm and 150 g), but had extremely small OFC (11.3 cm, -5.85 SD; 5 th percentile at this GA: 14.1 cm (24). Fetal autopsy confirmed severe micro-and brachycephaly, short receding forehead, narrow fontanelles and associated facial dysmorphisms including hypertelorism, anteverted nares, long philtrum, small upper lip and small outer ears. Moreover, the fetus showed . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted September 28, 2021. ; https://doi.org/10.1101/2021.09.25.21261035 doi: medRxiv preprint contractures according to early manifestation of arthrogryposis. Brain volume was reduced (about 40% of standard value, GA norm (16)) with slight enlargement of the ventricles.
The frontal lobes were hypoplastic, occipital lobes were shortened and appeared wing-like.
Temporo-parietal lobulation and corpus callosum were missing. Cerebellum was hypoplastic with a diameter of 1.2 cm (5 th percentile for GA: 1.6 cm (24)). All examined brain sections appeared histologically normal. The findings were interpreted as micrencephaly without neuronal migration disorder or structural malformations.
Examination of the eyes revealed discrete anisophthalmia and anisocoria with partial lens luxation due to dysplasia of the iris and persistent hyaloid artery of the left eye. There were no signs for infectious, toxic or hypoxic influences and no further organ abnormalities.
Fetal autopsy results of P2 were suspected as monogenic syndromic type of microcephaly.

Genetic analyses
Initial genetic investigations in P1, including conventional karyotyping, chromosomal

Postnatal phenotype in two unrelated individuals
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Analysis of missense variants in the FHA domain
The PNKP protein consists of the N-terminal FHA domain, a Linker region, connecting the According to Missense3D, the c.302C>T, p.(Pro101Leu) variant triggers a local steric clash alert ( Figure 2D). Additionally, based on the DynaMut web server predictor, the effect of this variant is stabilizing and the Δ Vibrational Entropy Energy between wildtype . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted September 28, 2021. ; https://doi.org/10.1101/2021.09.25.21261035 doi: medRxiv preprint and mutant structure is predicted to slightly decrease molecule flexibility ( Figure 2D). The c.311T>C, p.(Leu104Pro) substitution introduces a buried proline in the core of the protein domain, which tends to be particularly damaging with its restricted backbone conformation ( Figure 2E). In fact, according to Dynamut, the variant has the opposite effect of p.(Pro101Leu) and is predicted to be destabilizing with increase of the molecule flexibility ( Figure 2E). Comprehensive results of computational analyses for all variants in the FHA domain are listed in Table 2.

RNA analyses
To determine the effect of the parentally inherited c.498G>A variant in P1 and P2, we Despite the various PNKP-associated phenotypes (5-9,38-41) prenatal presentations in humans, especially noticeable concomitant brain anomalies, were unreported to date. Our compilation of prenatal diagnostic procedures and fetal pathological examination of P1 and P2 revealed neurodevelopmental and neurodegenerative brain alterations comparable to those described in mouse models with neuronal tissue-specific inactivation of PNKP (32).
These include general hypoplasia of different cerebral and cerebellar regions, without a . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. Given the OFC values observed immediately after birth in P3 and P4, the presence of prenatal microcephaly in these individuals is obvious. The diagnosis of PM relies on ultrasound (US) measurements in comparison to distributions for the respective gestational age and exclusion of exogenous causes such as infectious diseases. Until the recent Zika virus outbreak, no international standards and guidelines had been defined and research on the diagnostic performance of US measurements for fetal microcephaly was hampered by the overall rarity of the condition (33). Performance for prenatal US diagnosis seems good at the more extreme ends (< -4 SD) or when additional brain anomalies are present (33), like in the case of the two fetuses of family 1. Additionally, improved US technology and specialist training during the last years might have enabled these prenatal diagnoses.
We recommend interpreting the phenotypic presentations associated with biallelic variants in PNKP as continuous spectrum instead of the separated clinical entities MCSZ, AOA4 and CMT2B2, in agreement with previous suggestions (28). In this disease model, our . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review) preprint
The copyright holder for this this version posted September 28, 2021. ; Figure 4) clinical assessment will be needed to reach precise genotype-phenotype prediction.
Here, we extended the "PNKP-associated disorder continuum" to the prenatal period and complemented missense variant interpretation with 3D structure analysis and presented the first RNA-analyses known so far for PNKP variants. Future studies will need to combine these techniques with detailed phenotyping taking into account the PNKP-continuum and ideally add massive parallel functional tests. Due to the associated genetic heterogeneity the diagnosis of primary microcephaly requires (trio-) exome sequencing. The knowledge of distinct fetal phenotypes will be helpful for genetic variant assessment, especially those with unknown significance. Only with knowledge of variant pathogenicity and expected symptoms, we will be able to improve counselling in the prenatal setting, management in the postnatal period, prenatal diagnosis in subsequent pregnancies, and finally enable potential evaluation of treatment in the future.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

DATA AVAILABILITY
All data generated or analyzed during this study can be found in the online version of this article at the publisher's website with access to supplementary data resources or on Zenodo.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.  . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint   . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint