Autosomal allelic inactivation at loci with variable replication timing and dosage sensitivity

  1. Michael B Heskett
  2. Athanasios E Vouzas
  3. Brian Johnstone
  4. Krister P Freese
  5. Phillip A Yates
  6. Philip F Copenhaver
  7. Paul T Spellman
  8. David M Gilbert
  9. Mathew J Thayer  Is a corresponding author
  1. Division of Hematology Oncology, Departments of Medicine and Human Genetics, David Geffen School of Medicine, University of California, Los Angeles, United States
  2. San Diego Biomedical Research Institute, United States
  3. Department of Orthopedics and Rehabilitation, Oregon Health and Science University, United States
  4. Shriners Hospitals for Children - Portland, United States
  5. Department of Chemical Physiology and Biochemistry, Oregon Health and Science University, United States
  6. Department of Cell, Developmental and Cancer Biology, Oregon Health and Science University, United States
7 figures, 3 tables and 1 additional file

Figures

Figure 1 with 1 supplement
Genome-wide variable replication-timing and expression imbalance.

(a) Single-cell-derived ACP clones (color coded) were isolated, expanded, and processed for Repli-seq and RNA-seq. The maternal and paternal haplotypes were used to analyze allele-specific replication timing and expression. (b) The top panel illustrates the expected replication timing profile in multiple clones across an autosomal region. The regions that display VERT are highlighted by shading. The bottom panel illustrates the expected expression variability in different clones from the same individual. (c) We define VERT regions in Repli-seq data across each allele for each 250-kb genomic window. Windows with an SD value ≥2.5 standard deviations above the genome-wide median were classified as outliers (top panels). We define AEI as an allelic bias that is ≥80% allelic imbalance (AEI ≥0.80 or ≤0.20). The bottom panel shows both autosomal and X-linked expression data. (d) Illustration of an Inactivation/Stability Center (I/SC) showing the clonal and allele-restricted variability in replication timing and expression of both protein-coding and noncoding genes. (e) The standard deviation (Std. Dev.) in 250-kb windows (circles) across human chromosome 1. Outlier windows from all four sets of clones (ACP6, ACP7, EB3-2, and GM12878) are highlighted in different colors as shown. The top panel illustrates the G-banding (blue shading) pattern for human chromosome 1, with 23 VERT regions highlighted in red. The location of an imprinted region with asynchronous replication, containing the maternally expressed TP73 gene (https://geneimprint.com/), is marked with an +. The arrow marks the VERT region shown in panels f–h. (f–h) I/SC located between 239 and 242 Mb of chromosome 1 (see panel e). (f) VERT (shaded region) and AEI of genes (rectangles) detected in the ACP6 clones; each clone was color coded as shown. For replication, the solid lines represent the paternal (Hap1) allele, and the dotted line represents the maternal (Hap2) allele. (g) VERT (shaded region) and AEI of genes (rectangles) detected in the ACP7 clones; each clone was color coded as shown. (h) UCSC Genome Browser view of the genomic region in panels f and g. The shaded area highlights the VERT region. The asterisks mark the genes highlighted in panels f and g. (i) Venn diagram illustrating the number of VERT regions that overlap between the different clonal populations. (j) Gene Ontology (GO) enrichment analysis of the genes located within the VERT regions; the most significant hits for each set of clones, organized by cell type, are shown. GO analysis of the mouse genes located within the pre-B cell clone VERT regions is also shown.

Figure 1—source data 1

This file lists the location of the VERT regions detected in the human genome.

https://cdn.elifesciences.org/articles/109938/elife-109938-fig1-data1-v1.xlsx
Figure 1—source data 2

This file shows the number of haplotype phased sequence reads detected by RNAseq within coding genes.

https://cdn.elifesciences.org/articles/109938/elife-109938-fig1-data2-v1.xlsx
Figure 1—source data 3

This file shows the number of haplotype phased sequence reads detected by RNAseq within non-coding genes.

https://cdn.elifesciences.org/articles/109938/elife-109938-fig1-data3-v1.xlsx
Figure 1—source data 4

This file contains the number of transcripts derived from the early or late replicating alleles for Figures 1, 2 and 6 and Figure 1—figure supplement 1.

https://cdn.elifesciences.org/articles/109938/elife-109938-fig1-data4-v1.xlsx
Figure 1—source data 5

This file contains the list of human protein coding genes located within VERT regions.

https://cdn.elifesciences.org/articles/109938/elife-109938-fig1-data5-v1.xlsx
Figure 1—source data 6

This file contains the Gene Ontology analysis for the genes located in VERT regions.

https://cdn.elifesciences.org/articles/109938/elife-109938-fig1-data6-v1.xlsx
Figure 1—figure supplement 1
Imprinted regions display asynchronous replication.

(a, d) VERT regions on human chromosomes 14 and 20. The standard deviation in 250-kb windows (circles) is shown. Outlier windows from all four sets of clones (ACP6, ACP7, EB3-2, and GM12878) are highlighted in different colors as shown. The regions highlighted by the red boxes mark known imprinted loci. (b, e) For replication timing, solid lines represent the paternal allele (Hap1), and the dotted line represents the maternal allele (Hap2). The regions with asynchrony are shaded, and the gene expression is shown for each clone in a different color. (c) UCSC Genome Browser view of the asynchronous region in panel a. (b, c) The single asterisk marks the paternally expressed DLK gene, and the double asterisks mark the maternally expressed MEG8 gene. (f) UCSC Genome Browser view of the VERT region in panel d. (e, f) The # marks the paternally expressed L3MBTL1 gene.

Examples of Inactivation/Stability Centers (I/SCs).

(a, d, g) VERT regions on human chromosomes 18, 6, and 9, respectively. The top panels illustrate the G-banding (blue shading) pattern for each human chromosome, with VERT regions highlighted in red. The standard deviation (Std. Dev.) in 250-kb windows (circles) across each human chromosome is shown. Outlier windows from all four sets of clones (ACP6, ACP7, EB3-2, and GM12878) are highlighted in different colors as shown. The arrows and red boxes mark the VERT regions shown in panels b, e, h. (b, e, h) I/SCs highlighted in panels a, d, and g. For replication, the solid lines represent the paternal (Hap1) allele, and the dotted lines represent the maternal (Hap2) allele. The regions with VERT are shaded, and the AEI is shown for each clone in different colors. The location of an imprinted region with asynchronous replication, containing the paternally expressed LIN28B gene (https://geneimprint.com/), is marked with an + in panel d. (c, f, i) UCSC Genome Browser views of the genomic regions in panels b, e, and h. The shaded areas highlight the VERT regions. Genes with AEI are marked with +, *, or #. (b, c) The + marks MOCOS, which shows AEI in ACP cells (Figure 1—source data 2), and ++ marks FHOD3, which shows AEI in LCLs (Heskett et al., 2022). (e, f) The * marks transcripts from DNAH8, ** marks transcripts from KCNK5, and *** marks non-coding transcripts from ENSG00000293066 (also named: TL:6–39.9_1 Heskett et al., 2022). All three of these genes show AEI in LCL clones (Heskett et al., 2022). (h, i) The # marks transcripts from DAPK1, ## marks transcripts from SPIN1. Both genes show AEI in LCL clones (Heskett et al., 2022).

Figure 2—source data 1

This file contains the list of VERT regions that contain genes with AEI.

https://cdn.elifesciences.org/articles/109938/elife-109938-fig2-data1-v1.xlsx
Human gene clusters with AEI map to VERT regions.

(a, c) VERT regions on human chromosomes 19 and 6, respectively. The top panels illustrate the G-banding (blue shading) pattern for each human chromosome, with VERT regions highlighted in red. The standard deviation (Std. Dev.) in 250-kb windows (circles) across each human chromosome is shown. Outlier windows from all four sets of clones (ACP6, ACP7, EB3-2, and GM12878) are highlighted in different colors as shown. The red boxes and arrows identify the Inactivation/Stability Centers (I/SCs) in panels b and d. The location of an imprinted region with asynchronous replication, containing the paternally expressed LIN28B gene (https://geneimprint.com/), is marked with an + in panel c. (b) UCSC genome browser view illustrating the genomic location highlighted in panel a. The shaded area highlights the VERT region. The red underlines mark 11 olfactory receptor genes, and the green underlines mark 6 cytochrome P450 family (CYP4F) genes. (d) UCSC Genome Browser view illustrating the genomic location highlighted in panel c. The shaded area highlights the VERT region. The red underlines mark 16 HLA genes, and the green underlines mark 6 Lymphocyte antigen 6 family (LY6G) genes.

Figure 4 with 1 supplement
Mouse gene clusters with AEI map to VERT regions.

(a, d) VERT regions on mouse chromosomes 2 and 18, respectively. The standard deviation (Std. Dev.) in 50-kb windows (circles) across each mouse chromosome are shown. Outlier windows from pre-B cell clones are highlighted in orange. The red boxes and arrows identify the VERT regions in panels b (i) and c (ii). (b, c) UCSC genome browser views illustrating the genomic locations represented by i and ii in panel a. (b) Genome Browser view of the VERT region showing the location of 29 olfactory receptor genes (Or genes). The shaded area represents the VERT region. (c) Genome Browser view of the VERT region in ii above, showing the location of 242 olfactory receptor genes (Or genes). The shaded area represents the VERT region. (e) Genome Browser view of the VERT region in panel d, showing the location of clustered Protocadherin genes (15 Pcdha, 23 Pcdhb, and 22 Pcdhg genes). The shaded area represents the VERT region.

Figure 4—source data 1

This file contains the list of mouse genes located within VERT regions in the mouse genome.

https://cdn.elifesciences.org/articles/109938/elife-109938-fig4-data1-v1.xlsx
Figure 4—figure supplement 1
VERT detected at the mouse T cell receptor alpha locus.

UCSC Genome Browser view of the mouse VERT region (shaded) detected at the T cell receptor alpha (Tra) locus on mouse chromosome 14. The location of 37 Or genes, marked by red underlines, and the vomeronasal receptor gene Vmn2r88, marked by a green underline.

Figure 5 with 1 supplement
Human and mouse syntenic regions display VERT.

The short arm of human chromosome 5, with 11 VERT regions (panel a) and the centromeric region of mouse chromosome 15, with 4 VERT regions (panel b) are shown. (a) Two human VERT regions are highlighted, hi and hii, and are syntenic with locations on mouse chromosome 15, mi and mii, that also show VERT (panel b). Note that the gene order in both syntenic regions are inverted with respect to the telomeres when comparing the human and mouse loci. (c) Illustration of the I/SC on human chromosome 5 between 7.75 and 10.5 mb (hi). The shaded area marks the VERT region, and AEI of the coding gene SEMA5A in different clones is highlighted. For replication, the paternal (Hap1) and maternal (Hap2) alleles are indicated. (d) UCSC Genome Browser view of the human VERT region (shaded) in hi, showing the location of SEMA5A and a taste receptor gene (TAS2R1). (e) UCSC Genome Browser view of the mouse VERT region (shaded) in mi, showing the location of Sema5a and a taste receptor gene (Tas2r119). (f) UCSC Genome Browser view of the human VERT region (shaded) in hii above, highlighting the location of the AEI gene RETREG1 (Gimelbrant et al., 2007). (g) UCSC Genome Browser view of the mouse VERT region (shaded) in mii, showing the synteny between human and mouse for the protein-coding genes in panel f, including Retreg1. (h) VERT regions on human chromosome 16. The standard deviation (Std. Dev.) in 250 kb windows (circles) across human chromosome 16 is shown. Outlier windows are highlighted in different colors as shown. The red box and arrow mark the genomic location in panel i. (i) UCSC Genome Browser view illustrating the genomic location represented in panel h, showing the location of GRIN2A (red box). (j) VERT regions on mouse chromosome 16. The standard deviation in 50-kb windows (circles) is shown. Outlier windows from pre-B cell clones are highlighted in orange. (k) UCSC Genome Browser view illustrating the genomic location of the VERT region (shaded area) represented in panel j above, with the location of Grin2a highlighted with a red box.

Figure 5—source data 1

This file contains the list of mouse genes within VERT regions that have corresponding human genetic diseases.

https://cdn.elifesciences.org/articles/109938/elife-109938-fig5-data1-v1.xlsx
Figure 5—figure supplement 1
Mouse syntenic regions display VERT.

(a, c, e) VERT regions on mouse chromosomes 16, 2, and 18. The standard deviation in 50-kb windows (circles) is shown. Outlier windows from pre-B cell clones are highlighted in orange. The VERT regions highlighted by red boxes are expanded in panels b, d, and f. (b) UCSC Genome Browser view of the VERT region in a above, highlighting the location of the mouse App gene (red box). The shaded area represents the VERT region. (d) UCSC Genome Browser view of the VERT region in panel c, highlighting the location of the mouse Mocos and Fhod3 genes (red boxes). The shaded area represents the VERT region. (f) UCSC Genome Browser view of the VERT region in panel c, highlighting the location of five mouse Scn genes (red boxes). The shaded area represents the VERT region.

Epilepsy genes map to Inactivation/Stability Centers (I/SCs).

(a, d, g) VERT regions on human chromosomes 2, 5, and 8, respectively. The top panels illustrate the G-banding (blue shading) pattern for each human chromosome, with VERT regions highlighted in red. The standard deviation (Std. Dev.) in 250-kb windows (circles) across each human chromosome is shown. Outlier windows from all four sets of clones (ACP6, ACP7, EB3-2, and GM12878) are highlighted in different colors as shown. The location of an imprinted region with asynchronous replication, containing the paternally expressed DLGAP2 gene (https://geneimprint.com/), is marked with an + in panel g. The arrows and red boxes mark the VERT regions shown in panels b, c, e, f, and h–k. (b, e, h, j) Illustrations of I/SCs on chromosomes 2, 5, and 8, respectively. The shaded areas mark the VERT regions. Each clone was color coded to illustrate the VERT and AEI in different clones. The paternal (Hap1) and maternal (Hap2) alleles are indicated for replication and expression. (c) UCSC Genome Browser view illustrating the genomic location of the VERT region (shaded area) represented in panels a and b, with the location of five sodium channel genes (SCN3A, SCN2A, SCN1A, SCN9A, and SCN7A) highlighted with red underlines. (e) Illustration of an I/SC on chromosome 5. The shaded area marks the VERT region. AEI of GABRA1 is shown. (f) UCSC Genome Browser view illustrating the genomic location of the VERT region (shaded area) represented in panels d and e above, with the location of four GABA receptor genes (GABRB2, GABRBA6, GABRA1, and GABRG2) highlighted with red underlines. (h) Illustration of an I/SC (A and red box) on chromosome 8. Each ACP6 clone was color coded to illustrate the VERT and AEI in different clones. The paternal (Hap1) and maternal (Hap2) alleles are indicated for replication and expression. The shaded area marks the VERT region. The location and expression of OXR1 is indicated. (i) UCSC Genome Browser view illustrating the genomic location of the VERT region (shaded area) represented in panels g (A) and h, with the location of the OXR1 gene highlighted with a red underline. We note that ANGPT1 (marked with *), located within this I/SC displays AEI (see Figure 1—source data 2). (j) Illustration of an I/SC (B and red box) on chromosome 8. AEI of SAMD12 is indicated. (k) UCSC Genome Browser view illustrating the genomic location of the VERT region (shaded area) represented in panels g (B) and j, with the location of the SAMD12 highlighted with a red underline.

Figure 6—source data 1

This file lists the human disease genes located within VERT regions.

https://cdn.elifesciences.org/articles/109938/elife-109938-fig6-data1-v1.xlsx
Parkinson disease genes map to Inactivation/Stability Centers (I/SCs).

(a, c, e, g) VERT regions on human chromosomes 4, 12, 1, and 6, respectively. The top panels illustrate the G-banding (blue shading) pattern for each human chromosome, with VERT regions highlighted in red. The standard deviation (Std. Dev.) in 250-kb windows (circles) across each human chromosome is shown. Outlier windows from all four sets of clones (ACP6, ACP7, EB3-2, and GM12878) are highlighted in different colors as shown. The location of an imprinted region with asynchronous replication, containing the paternally expressed LIN28B gene (https://geneimprint.com/), is marked with an + in panel g. The location of an imprinted region with asynchronous replication, containing the maternally expressed TP73 gene (https://geneimprint.com/), is marked with an + in panel e. The arrows and red boxes mark the VERT regions shown in panels b, d, f, and h, respectively. (b) UCSC Genome Browser view illustrating the genomic location of the VERT region (shaded area) represented in panel a, with the location of the SNCA gene highlighted with a red underline. (d) UCSC Genome Browser view illustrating the genomic location of the VERT region (shaded area) highlighted in panel c, with the location of the LRRK2 gene highlighted with a red underline. (f) UCSC Genome Browser view illustrating the genomic location of the VERT region (shaded area) represented in panel e, with the location of the DNAJC6 gene highlighted with a red underline; also highlighted with a green underline is the JAK1 gene. (h) UCSC Genome Browser view illustrating the genomic location of the VERT region (shaded area) represented in panel g, with the location of the PRKN gene highlighted with a red underline.

Tables

Table 1
VERT regions at gene clusters in the human and mouse genomes.
Human VERT regionGene cluster
ChrStartStopGenes# of genesStartStop
113,250,00015,000,000PRAMEF2312,768,93913,422,488
2101,633,538102,383,540IL1R5101,987,899102,401,443
2165,250,000166,000,000SCN5165,081,204166,499,088
2176,250,000176,500,000HOXD9176,091,062176,191,493
446,750,00047,500,000GABR446,035,58347,431,569
468,384,28268,634,282UGT968,525,28369,660,413
5651,700750,000SLC5468,8881,448,927
5148,000,000148,250,000SPINK7147,822,992148,343,066
5162,250,000162,750,000GABR4161,282,809162,159,566
631,250,00032,750,000HLA1530,489,50933,093,613
LY6631,669,25931,722,476
811,250,00012,392,491DEFB611,970,74712,320,078
880,587,76583,000,000FABP481,277,94281,536,065
8142,500,000142,750,000LY68142,695,707143,162,219
11104,750,000105,000,000CASP5104,883,181105,141,441
1523,357,40023,500,000GOLGA522,451,99523,451,545
1528,750,00029,000,000GOLGA628,344,60428,855,059
179,750,00010,500,000MYH610,297,55210,658,865
1915,000,00015,250,000OR1114,792,47615,952,060
CYP615,507,10515,936,041
1920,750,00021,000,000ZNF3719,667,77324,119,420
1936,313,70036,500,000ZNF2936,171,63237,791,237
1943,000,00043,250,000PSG1042,713,52143,272,861
2024,250,00024,750,000CST1123,429,90224,957,423
2130,627,68130,877,681KRTAP3330,276,20131,041,124
Mouse VERT regionGene cluster
ChrStartStopGenes# of genesStartStop
236,322,85736,622,857Or2836,225,62437,193,679
264,622,86666,522,868Scn565,287,04266,637,709
288,950,00089,950,000Or23985,160,85990,147,389
459,850,00059,950,000Mup2259,950,00061,824,646
591,150,00091,350,000Cxcl591,186,93891,335,955
595,950,00096,250,000Pramel2393,891,82296,243,866
640,650,00041,150,000Tas2r840,386,43043,196,843
Or2440,386,43043,196,843
666,650,00066,850,000Vmn1866,498,91166,758,688
7113,809,945113,909,945Or24113,790,926114,306,041
1173,450,00074,050,000Or3173,160,40374,180,092
1837,350,00037,750,000Pcdh2337,298,05237,684,088
Table 2
VERT regions showing synteny between human and mouse genomes.
Mouse VERTHuman syntenic regionHuman VERT
ChrStartStopSynteny (%)ChrStartStopChrStartStop
173,042,58573,542,58554.02216,836,864217,400,7212217,000,000217,250,000
188,681,14789,081,14733.02234,565,154235,053,8572234,000,000235,000,000
1189558318189,758,31849.81214,144,079214,410,9561214,326,657216,000,000
264,842,28766,742,28746.52164,610,983166,685,3532165,250,000166,000,000
318,204,16418,404,16433.7864,684,74264,876,026864,337,44365,087,765
357,253,49958,553,49932.63149,407,811150,704,6413149,282,213150,032,213
378,253,38578,653,38526.14159,655,594159,952,3254159,328,848159,828,848
382,453,38582,853,38532.64154,688,934155,210,9404154,500,000154,750,000
4105034592105,234,59262.8156,331,72456,556,339155,500,00056,750,000
4142286667142,686,66740.6114,022,91514,523,879113,250,00015,000,000
810,500,00011,200,00043.213108,978,718110,065,35713109,250,000109,500,000
824,029,54424,929,54437.3840,152,84441,090,451840,250,00040,750,000
839,817,68740,717,68728.6815,936,54516,964,053815,750,00017,000,000
9116969950117,469,95042.3328,837,64729,402,683329,000,00029,500,000
1146,927,32548,127,32535.25155,465,678156,698,4305155,500,000156,250,000
1233,715,13433,815,13436.1719,330,76719,416,599719,250,00019,500,000
12105278049105,578,04941.61495,853,63096,269,3871496,000,00096,250,000
1311,095,96411,672,61929.91237,711,340238,470,9221237,604,200238,750,000
1525,120,33125,920,33142.9516,534,41217,724,462516,250,00017,000,000
1532,220,39132,420,39160.759,377,0599,575,50358,749,8889,499,888
169,067,7719,667,77141.6169,476,92410,028,307169,750,00010,000,000
Table 3
Neurodevelopmental genes within VERT regions.
VERT geneEpilepsy disease genes (inheritance: autosomal recessive—AR; autosomal dominant—AD; loss of function—LoF; gain of function—GoF)Reference
ALG14Myasthenic syndrome (AR); intellectual developmental disorder with epilepsy (AR); myopathy, epilepsy, and progressive cerebral atrophy (AR)McKusick, 2025
CACNA2D1Developmental and epileptic encephalopathy (AR)McKusick, 2025
GABRB1Developmental and epileptic encephalopathy (AD, LoF, GoF)Hernandez et al., 2023
GRIN2AEpilepsy (AD, LoF)Reutlinger et al., 2010
NTRK2Developmental and epileptic encephalopathy (AD); obesity (AD)Hamdan et al., 2017
OXR1Cerebellar hypoplasia/atrophy, epilepsy, and global developmental delay (AR)McKusick, 2025
RAPGEF2Epilepsy, familial adult myoclonic (AD, TTTCA repeat expansion)Depienne et al., 2023
RNF13Developmental and epileptic encephalopathy (AD, GoF)Edvardson et al., 2019
SAMD12Epilepsy (AD, TTTCA repeat expansion)Depienne et al., 2023
SCN1ADevelopmental and epileptic encephalopathy (AD, LoF); Dravet syndrome (AD, LoF); generalized epilepsy with febrile seizures (AD, LoF)Ding et al., 2021
SCN2ADevelopmental and epileptic encephalopathy (AD, LoF, GoF)Wolff et al., 2017
VERT geneIntellectual developmental genes (inheritance: autosomal recessive—AR; autosomal dominant—AD; loss of function—LoF; gain of function—GoF)Reference
ALG14Intellectual developmental disorder with epilepsy (AR); myasthenic syndrome (AR); myopathy, epilepsy, and progressive cerebral atrophy (AR)McKusick, 2025
CEP104Intellectual developmental disorder (AR); Joubert syndrome (AR)McKusick, 2025
CNOT2Intellectual developmental disorder (AD, LoF)Uehara et al., 2019
DPP6Intellectual developmental disorder (AD, LoF)Liao et al., 2013
ELP2Intellectual developmental disorder (AR)McKusick, 2025
FMN2Intellectual developmental disorder (AR)McKusick, 2025
GRIA1Intellectual developmental disorder (AD, GoF; AR, LoF)Ismail et al., 2022
IMPA1Intellectual developmental disorder (AR)McKusick, 2025
IQSEC1Intellectual developmental disorder (AR)McKusick, 2025
MYT1LIntellectual developmental disorder (AD, LoF)Blanchet et al., 2017
PDZD8Intellectual developmental disorder (AR)McKusick, 2025
RSRC1Intellectual developmental disorder (AR)McKusick, 2025
TAF4Intellectual developmental disorder (AD, LoF)Bertoli-Avella et al., 2021
TNIKIntellectual developmental disorder (AR)McKusick, 2025
TUSC3Intellectual developmental disorder (AR)McKusick, 2025
ZBTB18Intellectual developmental disorder (AD, LoF)Cohen et al., 2017
VERT geneNeurodevelopment disease genes (inheritance: autosomal recessive—AR; autosomal dominant—AD; loss of function—LoF; gain of function—GoF)Reference
ATP8A2Cerebellar ataxia (AR)McKusick, 2025
CSNK2BPoirier–Bienvenu neurodevelopmental syndrome (AD, LoF)Sakaguchi et al., 2017
DPH5Neurodevelopmental disorder (AR)McKusick, 2025
KCNA1Episodic ataxia/myokymia syndrome (AD, LoF, GoF)Paulhus et al., 2020; Müller et al., 2023
LMBRD2Developmental delay with variable neurologic and brain abnormalities (AD)McKusick, 2025
NBEANeurodevelopmental disorder (AD, LoF)Boulin et al., 2021
RETREG1Neuropathy (AR)McKusick, 2025
ROBO1Nystagmus (AR); neurooculorenal syndrome (AR); pituitary hormone deficiency (AD)McKusick, 2025
TIAM1Neurodevelopmental disorder (AR)McKusick, 2025
TRPM3Neurodevelopmental disorder (AD, GoF)Zhao et al., 2020
VARS1Neurodevelopmental disorder with microcephaly (AR)McKusick, 2025
ZSWIM6Acromelic frontonasal dysostosis (AD); neurodevelopmental disorder (AD, GoF)Palmer et al., 2017; Jaillard et al., 2011
VERT geneParkinson disease gene (inheritance: autosomal recessive—AR; autosomal dominant—AD; loss of function—LoF; gain of function—GoF)Reference
DNAJC6Parkinson disease 19 (AR)McKusick, 2025
LRRK2Parkinson disease (AD, GoF)Ohta et al., 2013; Taymans et al., 2023
PRKNParkinson disease (AR)McKusick, 2025
SLC18A2Parkinsonism-dystonia (AR)McKusick, 2025
SNCAParkinson disease 1 and 4 (AD, LoF/GoF)Kobayashi et al., 2003; Konno et al., 2016; Blauwendraat et al., 2021

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Michael B Heskett
  2. Athanasios E Vouzas
  3. Brian Johnstone
  4. Krister P Freese
  5. Phillip A Yates
  6. Philip F Copenhaver
  7. Paul T Spellman
  8. David M Gilbert
  9. Mathew J Thayer
(2026)
Autosomal allelic inactivation at loci with variable replication timing and dosage sensitivity
eLife 15:RP109938.
https://doi.org/10.7554/eLife.109938.4