Parthenogenetic mosaicism: generation via second polar body retention and unmasking of a likely causative PER2 variant for hypersomnia
Clinical Epigenetics volume 13, Article number: 73 (2021)
Parthenogenetic mosaicism is an extremely rare condition identified only in five subjects to date. The previous studies indicate that this condition is mediated by parthenogenetic activation and is free from a specific phenotype ascribed to unmaking of a maternally inherited recessive variant in the parthenogenetic cell lineage.
We examined a 28-year-old Japanese 46,XX female with Silver-Russell syndrome and idiopathic hypersomnia. The results revealed (1) predominance of maternally derived alleles for all the differentially methylated regions examined; (2) no disease-related copy-number variant; (3) two types of regions for all chromosomes, i.e., four BAF (B-allele frequency) band regions with single major microsatellite peaks of maternal origin and single minor microsatellite peaks of non-maternal (paternal) origin, and six BAF band regions with single major microsatellite peaks of maternal origin and two minor microsatellite peaks of maternal and non-maternal (paternal) origin; (4) an unmasked extremely rare PER2 variant (c.1403G>A:p.(Arg468Gln)) with high predicted pathogenicity; (5) mildly affected local structure with altered hydrogen bonds of the p.Arg468Gln-PER2 protein; and (6) nucleus-dominant subcellular distribution of the p.Arg468Gln-PER2 protein.
The above findings imply that the second polar body retention occurred around fertilization, resulting in the generation of the parthenogenetic cell lineage by endoreplication of a female pronucleus and the normal cell lineage by fusion of male and female pronuclei, and that the homozygous PER2 variant in the parthenogenetic cells is the likely causative factor for idiopathic hypersomnia.
Normal mammalian development requires both maternal and paternal genomes, because of the presence of multiple imprinted genes . Thus, parthenogenetic mammals are incompatible with life because of absent expression of paternally expressed genes and excessive expression of maternally expressed genes, and androgenetic mammals are also lethal because of absent expression of maternally expressed genes and excessive expression of paternally expressed genes. However, such mammals could be viable in the presence of a co-existing biparental cell lineage. Indeed, parthenogenetic and androgenetic mosaicisms have been identified in the human [2,3,4,5,6,7].
To our knowledge, parthenogenetic mosaicism has been identified in five patients [2,3,4,5,6]. Molecular data informative for the assessment of the generation mechanism invariably imply that parthenogenetic activation has taken place around the time of fertilization, resulting in the generation of the 46,XX parthenogenetic cell lineage by endoreplication of a female pronucleus and the biparental cell lineage by fusion of male and female pronuclei [2, 4,5,6]. Thus, although other mechanisms such as retention of second polar body and failure of the paternal pronucleus to duplicate, followed by endoreplication of a female pronucleus and fusion of biparentally derived pronuclei, have been assumed , there is no report describing such generation mechanisms. Furthermore, while such patients frequently exhibit growth failure with some Silver-Russell syndrome (SRS) features [2,3,4,5,6], no specific phenotype ascribed to unmasking of a recessive pathogenic variant(s) of maternal origin in the parthenogenetic cell lineage has been described, probably due to the buffering effect of a co-existing paternally derived wildtype (WT) allele or the absence of such a variant(s) in the previously reported patients.
Here, we report parthenogenetic mosaicism mediated by second polar body retention and an unmasked PER2 (period circadian regulator 2) variant in a female with SRS and idiopathic hypersomnia (IH).
This Japanese female was conceived naturally and was born at 42 weeks of gestation after an uncomplicated pregnancy and delivery. The parents were both 29 years of age at her birth and were non-consanguineous and healthy, with the paternal height of 168 cm (− 0.5 SD) and the maternal height of 156 cm (− 0.4 SD). The elder brother was also healthy.
She exhibited SRS-compatible clinical features satisfying the Netchine-Harbison criteria  in infancy to early childhood (Table 1 and Fig. 1a). She also exhibited overt hypotonia and borderline developmental delay: she controlled her head at six months of age, sat without support at nine months of age, walked without support at 23 months of age, and spoke single words at two years of age. She received growth hormone therapy for short stature from 6 to 12 years of age, and gonadotropin releasing hormone analog therapy for central precocious puberty with menarche at 7.5 years of age (the menarchial age in Japanese girls, 12.25 ± 1.25 years) from 7 to 12 years of age. She also had hypersomnia from childhood and truncal obesity and oligomenorrhea from late teens.
At 28 years of age, she was referred to us for the genetic studies of SRS. At that time, she manifested short stature, relative macrocephaly, hemihypoplasia, truncal obesity, borderline intellectual disability, and mild hypotonia, together with several other features (Table 1 and its Footnote-6 and Fig. 1a). Biochemical and radiological findings were grossly normal, except for mild fatty liver, low blood estradiol (24 pg/mL) (reference values: follicular phase 29 – 197 pg/mL, luteal phase 44 – 492 pg/mL), and hypoplastic ovaries and uterus. Her karyotype was 46,XX in all the 100 lymphocytes examined. Estrogen and progesterone replacement therapy was started, successfully inducing menses.
We performed detailed examinations for hypersomnia which was specific to her and compromised her quality of life. She had prolonged nocturnal sleep, daytime sleepiness, unrefreshing naps, and difficulty in awaking, and had no episode suggestive of cataplexy or severe snoring. Orexin concentration was 242 pg/mL (reference value, > 110 pg/mL) in the cerebrospinal fluid, and HLA typing indicated DQB1*03:01/DQB1*06:01 genotype and excluded the narcolepsy-associated DQB1*06:02 allele . Polysomnography (PSG) revealed a total sleep time of 505 min (age-matched average sleep time, 460 min) with normal sleep structure, apnea-hypoxia index of 2.9/h (reference value, < 5/h), and 3% oxygen desaturation index of 5.3/hour (reference value, < 15/h). Multiple sleep latency test (MSLT) showed remarkably short sleep latency of 1.5 min, and no sleep onset rapid eye movement period (SOREMP). These findings, especially the sleep-related symptoms and the combination of sleep latency of < 8 min and no SOREMP on PSG/MSLT, indicated IH rather than sleep apnea syndrome, narcolepsy (both type 1 with cataplexy and type 2 without cataplexy), and circadian rhythm sleep-wake disorders such as advanced or delayed sleep-wake phase disorder [9,10,11]. After thorough consultation, she received modafinil treatment , which appeared effective for the first ~ two months but became obviously ineffective thereafter. Thus, it was discontinued after six months of therapy.
We obtained leukocyte gDNA samples of the patient, the mother, and the brother, and salivary cell gDNA sample of the patient. Unfortunately, paternal sample was not available, because he had been deceased.
Pyrosequencing analysis showed abnormal, but not extremely skewed (e.g., < 10% or > 90%), methylation indices (MIs, the ratios of methylated clones) not only for SRS-related differentially methylated regions (DMRs) but also for other imprinting disease-related DMRs (Fig. 1b). The MIs were consistent with the predominance of maternally derived DMRs.
Genomewide array-based analysis
CytoScan HD revealed no known pathogenic copy-number variant (CNV) or CNV of unknown significance absent from public databases. B allele frequency (BAF) analysis showed two types of regions for all chromosomes, i.e., one type of regions with four bands for different BAFs (0.00, ~ 0.15, ~ 0.85, and 1.00) and the other type of regions with six bands for different BAFs (0.00, ~ 0.15, ~ 0.30, ~ 0.70, ~ 0.85, and 1.00) (Fig. 1c and Additional file 1: Fig. S1). The pericentromeric portions were invariably composed of four BAF band regions.
Microsatellite analysis delineated a different pattern between four and six BAF band regions (Fig. 1d and Additional file 1: Table S1). For loci on four BAF band regions (shown in blue), informative results were obtained for several loci, showing two peaks in the patient, one major peaks consistent with maternal isodisomy and one minor peaks of non-maternal (paternal) origin. For loci on six BAF band regions (shown in red), informative results were also obtained for several loci, showing three peaks in the patient, one major peaks consistent with maternal isodisomy and two minor peaks of maternal and non-maternal (paternal) origin. The frequency of the parthenogenetic cells was calculated as 70% in leukocytes and 67% in salivary cells, using informative microsatellite data for loci on four BAF band regions (e.g., D2S2211) (Additional file 1: Note).
Whole exome sequencing
Whole exome sequencing was performed, to identify an underlying genetic factor(s) for IH specific to this patient. Consequently, a missense variant on the third nuclear export signal (NES3) of PER2 (NM_022817.2:c.1403G>A:p.(Arg468Gln)) (https://www.ncbi.nlm.nih.gov/genbank) was identified as the most likely underlying factor for IH, because of its extreme rarity in the public and in-house databases (registered in gnomeAD_exome_EAS in a heterozygous condition with the allele frequency of 0.0001604), high predicted pathogenicity, and relevance to the circadian rhythm  (Fig. 2a–c). This variant was detected in 55 of 66 reads (83.3%) in this patient and in 25 of 52 reads (48.0%) in the mother. The frequency was well consistent with the genomic position of PER2 on the four BAF band region where BAF is assumed to be ~ 85% (Fig. 2d). No other rare and pathogenic variants were identified in sleep-related genes .
Protein modelling analysis
In the mouse crystal structure of a PAS domain fragment (Protein Data bank ID code: 3GDI), p.Arg466Gln, which corresponds to p.Arg468Gln in the human, was evaluated to mildly affect the local protein structure and hydrogen bonds (Fig. 2e).
Subcellular distribution analysis
Since this variant resided at NES3, we compared subcellular distribution pattern between the GFP-labelled WT-PER2 and p.Arg468Gln-PER2 proteins. The p.Arg468Gln-PER2 protein was mildly but significantly more localized to the nucleus than the WT-PER2 protein (Fig. 2f).
Molecular studies indicate that this patient had parthenogenetic mosaicism. It is inferred that the second polar retention occurred around the time of fertilization, resulting in the generation of the parthenogenetic cell lineage by endoreplication of a female pronucleus and the normal cell lineage by fusion of male and female pronuclei (Fig. 1e)  (thus, in a strict sense, this condition is neither mosaicism caused by a mitotic error of a single zygote nor chimerism resulting from the fusion of two different zygotes). In this case, this patient is assumed to have three haploid sets of maternal origin (two in the parthenogenetic cell lineage and one in the biparental cell lineage) consisting of homozygous and heterozygous regions, with pericentromeric portions being composed of homozygous regions, together with one haploid set of paternal origin (Fig. 1e). This explains why the MIs were outside the normal range but were not extremely skewed, and why this patient had not only four BAF band regions with single major microsatellite peaks of maternal origin and single minor microsatellite peaks of non-maternal (paternal) origin but also six BAF band regions with single major microsatellite peaks of maternal origin and two minor microsatellite peaks of maternal and non-maternal (paternal) origin (Additional file 1: Fig. S2). The usefulness of SNP array in the elucidation of generation mechanisms of such rare conditions has been reported previously [4, 6]. Indeed, CytoScan HD analysis has confirmed the underlying generation mechanisms in the previously reported parthenogenetic and androgenetic mosaic patients (Additional file 1: Figs. S3 and S4) [5, 7].
Several findings are notable with regard to the parthenogenetic mosaicism of this patient. First, the mosaic ratio estimated by microsatellite analysis was well consistent with the BAFs obtained by SNP array (Additional file 1: Fig. S2). This would argue for the accuracy of the molecular studies employed in this study. Second, parthenogenetic cells accounted for the majority of leukocytes, despite the apparently global negative selective pressure of the parthenogenesis . This would primarily be due to the bias in the human study. Indeed, genetic abnormalities would not be identified if they are barely present in leukocytes. In addition, parthenogenetic cells might be well preserved in several tissues of subjects with hitherto unknown specific genetic and environmental factors. Third, the mosaic ratio was similar between leukocytes and salivary cells. This would probably be co-incidental, because it was rather different between the two tissues in the previously reported patient (case 4 in Table 1) . Indeed, since adult saliva cells primarily consist of leukocytes of mesodermal origin and epithelial cells of ectodermal origin with a roughly similar frequency , the mosaic ratio could be variable between the leukocytes and salivary cells.
To date, parthenogenetic mosaicism has been identified in six subjects including this patient (Table 1) [2,3,4,5,6]. Although clinical information remains fragmentary, they primarily exhibit growth failure with some SRS features including hemihypoplasia. This would primarily be consistent with the occurrence of such phenotypes in upd(7)mat, upd(11)mat, upd(14)mat, upd(16)mat, and upd(20)mat [7, 15]. In this regard, hemihypoplasia would probably be due to laterality in the distribution of parthenogenetic cells, and disorders of sex development in cases 1 and 2 is explained by the co-existing 46,XY cell lineage. In addition, upd(15)mat might be relevant to the development of several features such as growth failure, intellectual disability, and hypotonia observed in Prader-Willi syndrome , although such non-specific features are also reported in SRS . For our patient, it is notable that she exhibited SRS-compatible phenotype in infancy and central precocious puberty and truncal obesity in a later age. Since such phenotype is characteristic of Temple syndrome caused by upd(14)mat , this may suggest an epi-dominant clinical effect of the chromosome 14q32.2 imprinted region and a relatively high frequency of the parthenogenetic cells in tissues/organs related to the phenotype of this patient.
Notably, this patient had IH and a c.1403G>A:p.(Arg468Gln) variant in PER2. PER2 is involved in the mammalian circadian clock system, together with cryptochrome (CRY), circadian locomotor output cycles kaput (CLOCK), and brain and muscle Arnt-like protein 1 (BMAL1) [18,19,20]. PER2 and CRY are positively expressed by the CLOCK-BMAL1 heterodimer, and PER2 and CRY proteins are exported from the nucleus to the cytoplasm. Then, the PER2-CRY complex generated in the cytoplasm is imported into the nucleus to inhibit their own transcriptions by directly interacting with the CLOCK-BMAL1 complex. This feedback cycle underlies the circadian rhythm (Fig. 2g and its legend) [18,19,20]. In this context, the p.Arg468Gln-PER2 protein of this patient distributed more frequently in the nucleus than the WT-PER2 protein, consistent with the variant being present at NES3. This would have resulted in the reduced PER2 and CRY production via a prolonged repressing effect of the PER2-CRY dimer on the CLOCK-BMAL1 complex, leading to the development of IH as a circadian rhythm disorder (Fig. 2g) . In this context, since the mother heterozygous for the c.1403G>A:p.(Arg468Gln) variant had no sleep-related phenotype, this variant would have a recessive effect. The presence of this extremely rare variant in a heterozygous condition in the gnomeAD_exome_EAS database would also support the recessive effect. It is likely, therefore, that the recessive effect of the c.1403G>A:p.(Arg468Gln) variant was unmasked in the parthenogenetic cells, and that parthenogenetic cells accounted for a considerable degree of cells in the sleep control region such as the suprachiasmatic nucleus . Thus, the c.1403G>A:p.(Arg468Gln) variant would be regarded as a likely causative variant for the development of IH. In addition, since PER2 exerts biological effects on multiple systems , the c.1403G>A:p.(Arg468Gln) variant might also be involved in the development of phenotypes other than IH, such as central obesity, of this patient.
Several matters would be worth pointing out for PER2. First, the subcellular distribution was not so markedly different, though significant, between the WT-PER2 and p.Arg468Gln-PER2 proteins. This would not be unexpected, because: (1) p.Arg468Gln was assessed to just mildly affect the local protein structure and hydrogen bonds (Fig. 2e); (2) p.Arg468Gln is an exchange between hydrophilic basic and polar amino acids, while hydrophobic amino acids (e.g., leucine and isoleucine) are likely to play a crucial role in the NES ; and (3) substitution of one leucine and one isoleucine residues at NES3 with hydrophobic alanine residues has exerted a significant but mild effect on the subcellular distribution as observed in this study (replacement of two leucine residues at NES1 and that of four leucine residues at NES2 with alanine residues have markedly and mildly affected the subcellular distribution, respectively) . It appears therefore, that such a relatively mild alteration in the subcellular distribution can cause a clinically recognizable sleep phenotype. Second, a pathogenic p.Ser662Gly variant has been identified at the CKIε domain in a family with dominantly inherited advanced sleep-wake phase disorder (ASWPD) . In contrast to the p.Arg468Gln, the p.Ser662Gly has been implicated to cause premature nuclear export of the p.Ser662Gly-PER2 protein [19, 20, 24, 25] (for details, see Additional file 1: Fig. S5). It is likely, therefore, that PER2 variants lead to different sleep phenotypes (IH and ASWPD), probably depending on the subcellular distribution pattern of the variant PER2 proteins. Lastly, significant associations have been identified between multiple SNPs and various sleep-wake phenotypes (Additional file 1: Table S2). This indicates that PER2 variants function not only as causative factors but also as susceptibility factors.
We identified the first case of parthenogenetic mosaicism that was generated via second polar body retention and unmasked a recessive PER2 variant for IH. This implies the presence of different generation mechanisms for parthenogenetic mosaicism and the unmasking of a maternal recessive variant by parthenogenesis.
Primers used are shown in Additional file 1: Table S3.
Pyrosequencing was performed for 59 CpG sites at nine DMRs involved in the development of SRS and other human imprinting disorders, using bisulfite-treated leukocyte gDNA samples. The procedure was as described in the manufacturer′s instructions (Qiagen). Subsequently, MIs were obtained using PyroMark Q24 (Qiagen). The data from 50 healthy subjects were utilized as references.
Genomewide array-based analysis
CNV and BAF analyses were performed using CytoScan HD (ThermoFisher Scientific) harboring more than 2.4 million CNV markers and approximately 750,000 SNP markers. The log ratio and B-allele frequency were calculated using Rawcopy R package . We selected CNVs larger than 100-kb, using log-ratio thresholds of > +0.2 for gains and < −0.3 for deletions, as suggested previously , and these CNVs were annotated using AnnotSV (Annotation and Ranking of Human Structural Variations) (https://www.lbgi.fr/AnnotSV/runjob). We searched for known pathogenic CNVs and CNVs of unknown significance absent from multiple public databases adopted in AnnotSV. CNVs located on segmental duplications were regarded as begin CNVs. Violin plot for BAFs was made with snp.profile.Rdata by ggplot2 (https://ggplot2.tidyverse.org/reference/index.html).
PCR was performed with a fluorescently labeled forward primer and an unlabeled reverse primer flanking each microsatellite locus. Subsequently, PCR products were determined for physical sizes and area under curves on an ABI PRISM 3100xl Genetic Analyzer using Gene Mapper software version 4.1 (ThermoFisher Scientific).
Whole exome sequencing
Whole exome sequencing was carried out using SureSelect Human All Exon V6 (Agilent Technologies). Captured libraries were sequenced by NextSeq 500 (Illumina) with 150 bp paired-end reads. Exome data processing, variant calling, and variant annotation were carried out, as described previously . Human GRCh37 was utilized as the reference genome, and the PER2 variant was described according to the HGVS recommendations , using NM_022817.2 as the reference sequence (https://www.ncbi.nlm.nih.gov/genbank). We extracted rare variants with minor allele frequencies of < 0.005 in all the following public databases and in-house database (n = 218): (1) whole genome and exome data for East Asian Population in Genome Aggregation Database (gnomAD_genome_EAS and gnomAD_exome_EAS) (http://gnomad.broadinstitute.org/); (2) Human Genetic Variation Database (HGVD) (http://www.hgvd.genome.med.kyoto-u.ac.jp/); and (3) whole-genome sequences of 8380 healthy Japanese individuals and construction of the highly accurate Japanese population reference panel (8.3KJPN) (https://ijgvd.megabank.tohoku.ac.jp/). Final variants were annotated with Annovar .
In silico pathogenicity predictions were carried out for identified rare variants, using the following methods: (1) CADD (Combined Annotation–Dependent Depletion) (http://cadd.gs.washington.edu/); (2) PP2_HVAR (Polyphen-2 Hum Var) (http://genetics.bwh.harvard.edu/pph2/); (3) SIFT (Sorting Intolerant From Tolerant) (http://sift.jcvi.org/); and (4) MutationTaster (http://www.mutationtaster.org/). We further examined the involvement of genes with rare variants in the circadian rhythm .
Sanger direct sequencing was performed for PCR products containing identified variants on the ABI 3130xl Genetic Analyzer (ThermoFisher Scientific). To confirm the variant, the PCR products were subcloned with the TOPO TA Cloning Kit (ThermoFisher Scientific), and wildtype and variant alleles were sequenced separately.
Protein modelling analysis
The impact of p.Arg468Gln on protein structure was examined for mouse crystal structure of a PAS domain fragment (Protein Data bank ID code: 3GDI) using Swiss PDB Viewer. The 468th Arg codon in the human corresponds to the 466th Arg codon in the mouse.
Subcellular distribution analysis
Full-length wildtype PER2 cDNA (WT-PER2) embedded in pF1K0F vector was obtained from Kazusa DNA Research Institute (Product ID FXC06329), and variant PER2 cDNA with c.1403G>A (c.1403G>A-PER2) was created by site-directed mutagenesis using PrimeSTAR Mutagenesis Basal Kit (Takara Bio). cDNAs for the WT-PER2 and c.1403G>A-PER2 were fused to the C-terminal side of the GFP gene in pEF1α-EGFP vector. Transfection was performed for COS7 cells, which were seeded into 6-well dishes containing 3 mL of DMEM with 10% fetal bovine serum (3.0 × 105 cells/well) and were maintained for 24 h, using lipofectamine 3000 (ThermoFisher Scientific). Transfected COS7 cells were further cultured for 24 h, and fluorescent signals were obtained using ZOE Fluorescent Cell Imager (Bio-Rad) with and without the 4,6 diamidino-2-phenylindole (DAPI) staining. Subcellular distribution patterns were evaluated as a nuclear localization pattern (N), a nuclear-cytoplasmic distribution pattern (N + C), and a cytoplasmic localization pattern (C) by a technician who did not know the purpose of this experiment. The experiments were repeated six times, and a total of ~ 2000 cells were observed for the subcellular distribution of the WT-PER2 and p.Arg468Gln-PER2 proteins. Statistical analysis of the subcellular distribution pattern between the WT-PER2 and p.Arg468Gln-PER2 proteins (2 × 3 contingency table) was performed with the two-tailed Fisher's exact probability test.
Availability of data and materials
All data generated or analyzed during this study are available from the corresponding author on reasonable request.
B allele frequency
Differentially methylated region
Multiple sleep latency test
Nuclear export signal
Period circadian regulator 2
Sleep onset rapid eye movement period
Bartolomei MS, Ferguson-Smith AC. Mammalian genomic imprinting. Cold Spring Harb Perspect Biol. 2011;3:a002592.
Strain L, Warner JP, Johnston T, Bonthron DT. A human parthenogenetic chimaera. Nat Genet. 1995;11:164–9.
Horike S, Ferreira JC, Meguro-Horike M, Choufani S, Smith AC, Shuman C, et al. Screening of DNA methylation at the H19 promoter or the distal region of its ICR1 ensures efficient detection of chromosome 11p15 epimutations in Russell-Silver syndrome. Am J Med Genet A. 2009;149A:2415–23.
Conlin LK, Thiel BD, Bonnemann CG, Medne L, Ernst LM, Zackai EH, et al. Mechanisms of mosaicism, chimerism and uniparental disomy identified by single nucleotide polymorphism array analysis. Hum Mol Genet. 2010;19:1263–75.
Yamazawa K, Nakabayashi K, Kagami M, Sato T, Saitoh S, Horikawa R, et al. Parthenogenetic chimaerism/mosaicism with a Silver-Russell syndrome-like phenotype. J Med Genet. 2010;47:782–5.
Bens S, Luedeke M, Richter T, Graf M, Kolarova J, Barbi G, et al. Mosaic genome-wide maternal isodiploidy: an extreme form of imprinting disorder presenting as prenatal diagnostic challenge. Clin Epigenetics. 2017;9:111.
Yamazawa K, Nakabayashi K, Matsuoka K, Matsubara K, Hata K, Horikawa R, et al. Androgenetic/biparental mosaicism in a girl with Beckwith-Wiedemann syndrome-like and upd(14)pat-like phenotypes. J Hum Genet. 2011;56:91–3.
Wakeling EL, Brioude F, Lokulo-Sodipe O, O’Connell SM, Salem J, Bliek J, et al. Diagnosis and management of Silver-Russell syndrome: first international consensus statement. Nat Rev Endocrinol. 2017;13:105–24.
Scammell TE. Narcolepsy. N Engl J Med. 2015;373:2654–62.
American Academy of Sleep Medicine. International classification of sleep disorders. 3rd ed. Darien: Elsevier; 2014.
Trotti LM. Idiopathic hypersomnia. Sleep Med Clin. 2017;12:331–44.
Lowrey PL, Takahashi JS. Genetics of circadian rhythms in mammalian model organisms. Adv Genet. 2011;74:175–230.
Fundele R, Norris ML, Barton SC, Reik W, Surani MA. Systematic elimination of parthenogenetic cells in mouse chimeras. Development. 1989;106:29e35.
Theda C, Hwang SH, Czajko A, Loke YJ, Leong P, Craig JM. Quantitation of the cellular content of saliva and buccal swab samples. Sci Rep. 2018;8:6944.
Luk HM, Ivan Lo FM, Sano S, Matsubara K, Nakamura A, Ogata T, et al. Silver-Russell syndrome in a patient with somatic mosaicism for upd(11)mat identified by buccal cell analysis. Am J Med Genet A. 2016;170:1938–41.
Cassidy SB, Schwartz S, Miller JL, Driscoll DJ. Prader-Willi syndrome. Genet Med. 2012;14:10–26.
Kagami M, Nagasaki K, Kosaki R, Horikawa R, Naiki Y, Saitoh S, et al. Temple syndrome: comprehensive molecular and clinical findings in 32 Japanese patients. Genet Med. 2017;19:1356–66.
Mieda M. The central circadian clock of the suprachiasmatic nucleus as an ensemble of multiple oscillatory neurons. Neurosci Res. 2020;156:24–31.
Vanselow K, Vanselow JT, Westermark PO, Reischl S, Maier B, Korte T, et al. Differential effects of PER2 phosphorylation: molecular basis for the human familial advanced sleep phase syndrome (FASPS). Genes Dev. 2006;20:2660–72.
Albrecht U, Bordon A, Schmutz I, Ripperger J. The multiple facets of Per2. Cold Spring Harb Symp Quant Biol. 2007;72:95–104.
Landzberg D, Trotti LM. Is idiopathic hypersomnia a circadian rhythm disorder? Curr Sleep Med Rep. 2019;5:201–6.
Xu D, Farmer A, Collett G, Grishin NV, Chook YM. Sequence and structural analyses of nuclear export signals in the NESdb database. Mol Biol Cell. 2012;23:3677–93.
Yagita K, Tamanini F, Yasuda M, Hoeijmakers JH, van der Horst GT, Okamura H. Nucleocytoplasmic shuttling and mCRY-dependent inhibition of ubiquitylation of the mPER2 clock protein. EMBO J. 2002;21:1301–14.
Toh KL, Jones CR, He Y, Eide EJ, Hinz WA, Virshup DM, et al. An hPer2 phosphorylation site mutation in familial advanced sleep phase syndrome. Science. 2001;291:1040–3.
Xu Y, Toh KL, Jones CR, Shin JY, Fu YH, Ptácek LJ. Modeling of a human circadian mutation yields insights into clock regulation by PER2. Cell. 2007;128:59–70.
Mayrhofer M, Viklund B, Isaksson A. Rawcopy: Improved copy number analysis with Affymetrix arrays. Sci Rep. 2016;6:36158.
Miyado M, Fukami M, Takada S, Terao M, Nakabayashi K, Hata K, et al. Germline-derived gain-of-function variants of Gsα-coding GNAS gene identified in nephrogenic syndrome of inappropriate antidiuresis. J Am Soc Nephrol. 2019;30:877–89.
den Dunnen JT, Dalgleish R, Maglott DR, Hart RK, Greenblatt MS, McGowan-Jordan J, et al. HGVS recommendations for the description of sequence variants: 2016 update. Hum Mutat. 2016;37:564–9.
Wang K, Li M, Hakonarson H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38:e164.
We are grateful to the patient and her family members for their cooperation. We would like to thank Dr. Eisaku Okada for the statistical analysis, Dr. Katsushi Tokunaga for the HLA typing, Dr. Takashi Kanbayashi for the measurement of orexin, Dr. Yumiko Ohkubo for providing us with some clinical data, and Ms. Aya Kitamoto for their technical support.
This work was supported by Grants from the Japan Agency for Medical Research and Development (AMED) (JP20ek0109301, JP20ek0109373, and JP20ek0109464), the Japan Society for the Promotion of Science (17H06428), and the Takeda Science Foundation.
This study was approved by the Institutional Review Board Committee at Hamamatsu University School of Medicine, and was performed after obtaining written informed consent from the patient.
Consent for publication
We obtained written permission to publish photographs and clinical and molecular data of the patient.
The authors declare no competing of interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: The results of microsatellite analysis. Table S2: Association of sleep-wake phenotypes with SNPs in PER2. Table S3: Primers utilized in this study. Fig. S1: Array-based copy-number and BAF analyses for each chromosome, using CytoScan HD. Fig. S2: Parthenogenetic mosaicism mediated by the second polar body retention. Fig. S3: Generation mechanism of the parthenogenetic mosaicism in the patient previously described by Yamazawa et al. Fig. S4: Generation mechanism of the androgenetic mosaicism in the patient described by Yamazawa et al. Fig. S5: The p.Ser662Gly PER2 variant identified in a dominantly inherited advanced sleep-wake phase disorder (ASWPD). Note: Calculation of the parthenogenetic cell frequency in leukocytes and salivary cells.
About this article
Cite this article
Masunaga, Y., Kagami, M., Kato, F. et al. Parthenogenetic mosaicism: generation via second polar body retention and unmasking of a likely causative PER2 variant for hypersomnia. Clin Epigenet 13, 73 (2021). https://doi.org/10.1186/s13148-021-01062-0