Variant Location Is a Novel Risk Factor for Individuals With Arrhythmogenic Cardiomyopathy Due to a Desmoplakin (DSP) Truncating Variant.

BACKGROUND: Truncating variants in desmoplakin (DSPtv) are an important cause of arrhythmogenic cardiomyopathy; however the genetic architecture and genotype-specific risk factors are incompletely understood. We evaluated phenotype, risk factors for ventricular arrhythmias, and underlying genetics of DSPtv cardiomyopathy. METHODS: Individuals with DSPtv and any cardiac phenotype, and their gene-positive family members were included from multiple international centers. Clinical data and family history information were collected. Event-free survival from ventricular arrhythmia was assessed. Variant location was compared between cases and controls, and literature review of reported DSPtv performed. RESULTS: There were 98 probands and 72 family members (mean age at diagnosis 43±8 years, 59% women) with a DSPtv, of which 146 were considered clinically affected. Ventricular arrhythmia (sudden cardiac arrest, sustained ventricular tachycardia, appropriate implantable cardioverter defibrillator therapy) occurred in 56 (33%) individuals. DSPtv location and proband status were independent risk factors for ventricular arrhythmia. Further, gene region was important with variants in cases (cohort n=98; Clinvar n=167) more likely to occur in the regions resulting in nonsense mediated decay of both major DSP isoforms, compared with n=124 genome aggregation database control variants (148 [83.6%] versus 29 [16.4%]; P<0.0001). CONCLUSIONS: In the largest series of individuals with DSPtv, we demonstrate that variant location is a novel risk factor for ventricular arrhythmia, can inform variant interpretation, and provide critical insights to allow for precision-based clinical management.

CAG Repeat Size Influences the Progression Rate of Spinocerebellar Ataxia Type 3.

OBJECTIVE: In spinocerebellar ataxia type 3/Machado-Joseph disease (SCA3/MJD), the expanded cytosine adenine guanine (CAG) repeat in ATXN3 is the causal mutation, and its length is the main factor in determining the age at onset (AO) of clinical symptoms. However, the contribution of the expanded CAG repeat length to the rate of disease progression after onset has remained a matter of debate, even though an understanding of this factor is crucial for experimental data on disease modifiers and their translation to clinical trials and their design. METHODS: Eighty-two Dutch patients with SCA3/MJD were evaluated annually for 15 years using the International Cooperative Ataxia Rating Scale (ICARS). Using linear growth curve models, ICARS progression rates were calculated and tested for their relation to the length of the CAG repeat expansion and to the residual age at onset (RAO): The difference between the observed AO and the AO predicted on the basis of the CAG repeat length. RESULTS: On average, ICARS scores increased 2.57 points/year of disease. The length of the CAG repeat was positively correlated with a more rapid ICARS progression, explaining 30% of the differences between patients. Combining both the length of the CAG repeat and RAO as comodifiers explained up to 47% of the interpatient variation in ICARS progression. INTERPRETATION: Our data imply that the length of the expanded CAG repeat in ATXN3 is a major determinant of clinical decline, which suggests that CAG-dependent molecular mechanisms similar to those responsible for disease onset also contribute to the rate of disease progression in SCA3/MJD. ANN NEUROL 2021;89:66-73.

NIPTeR: an R package for fast and accurate trisomy prediction in non-invasive prenatal testing.

BACKGROUND: Various algorithms have been developed to predict fetal trisomies using cell-free DNA in non-invasive prenatal testing (NIPT). As basis for prediction, a control group of non-trisomy samples is needed. Prediction accuracy is dependent on the characteristics of this group and can be improved by reducing variability between samples and by ensuring the control group is representative for the sample analyzed. RESULTS: NIPTeR is an open-source R Package that enables fast NIPT analysis and simple but flexible workflow creation, including variation reduction, trisomy prediction algorithms and quality control. This broad range of functions allows users to account for variability in NIPT data, calculate control group statistics and predict the presence of trisomies. CONCLUSION: NIPTeR supports laboratories processing next-generation sequencing data for NIPT in assessing data quality and determining whether a fetal trisomy is present. NIPTeR is available under the GNU LGPL v3 license and can be freely downloaded from https://github.com/molgenis/NIPTeR or CRAN.

A "late-but-fitter revertant cell" explains the high frequency of revertant mosaicism in epidermolysis bullosa.

Revertant mosaicism, or "natural gene therapy", is the phenomenon in which germline mutations are corrected by somatic events. In recent years, revertant mosaicism has been identified in all major types of epidermolysis bullosa, the group of heritable blistering disorders caused by mutations in the genes encoding epidermal adhesion proteins. Moreover, revertant mosaicism appears to be present in all patients with a specific subtype of recessive epidermolysis bullosa. We therefore hypothesized that revertant mosaicism should be expected at least in all patients with recessive forms of epidermolysis bullosa. Naturally corrected, patient-own cells are of extreme interest for their promising therapeutic potential, and their presence in all patients would open exciting, new treatment perspectives to those patients. To test our hypothesis, we determined the probability that single nucleotide reversions occur in patients' skin using a mathematical developmental model. According to our model, reverse mutations are expected to occur frequently (estimated 216x) in each patient's skin. Reverse mutations should, however, occur early in embryogenesis to be able to drive the emergence of recognizable revertant patches, which is expected to occur in only one per ~10,000 patients. This underestimate, compared to our clinical observations, can be explained by the "late-but-fitter revertant cell" hypothesis: reverse mutations arise at later stages of development, but provide revertant cells with a selective growth advantage in vivo that drives the development of recognizable healthy skin patches. Our results can be extrapolated to any other organ with stem cell division numbers comparable to skin, which may offer novel future therapeutic options for other genetic conditions if these revertant cells can be identified and isolated.

Lamin A/C-Related Cardiac Disease: Late Onset With a Variable and Mild Phenotype in a Large Cohort of Patients With the Lamin A/C p.(Arg331Gln) Founder Mutation.

BACKGROUND: Interpretation of missense variants can be especially difficult when the variant is also found in control populations. This is what we encountered for the LMNA c.992G>A (p.(Arg331Gln)) variant. Therefore, to evaluate the effect of this variant, we combined an evaluation of clinical data with functional experiments and morphological studies. METHODS AND RESULTS: Clinical data of 23 probands and 35 family members carrying this variant were retrospectively collected. A time-to-event analysis was performed to compare the course of the disease with carriers of other LMNA mutations. Myocardial biopsies were studied with electron microscopy and by measuring force development of the sarcomeres. Morphology of the nuclear envelope was assessed with immunofluorescence on cultured fibroblasts. The phenotype in probands and family members was characterized by atrioventricular conduction disturbances (61% and 44%, respectively), supraventricular arrhythmias (69% and 52%, respectively), and dilated cardiomyopathy (74% and 14%, respectively). LMNA p.(Arg331Gln) carriers had a significantly better outcome regarding the composite end point (malignant ventricular arrhythmias, end-stage heart failure, or death) compared with carriers of other pathogenic LMNA mutations. A shared haplotype of 1 Mb around LMNA suggested a common founder. The combined logarithm of the odds score was 3.46. Force development in membrane-permeabilized cardiomyocytes was reduced because of decreased myofibril density. Structural nuclear LMNA-associated envelope abnormalities, that is, blebs, were confirmed by electron microscopy and immunofluorescence microscopy. CONCLUSIONS: Clinical, morphological, functional, haplotype, and segregation data all indicate that LMNA p.(Arg331Gln) is a pathogenic founder mutation with a phenotype reminiscent of other LMNA mutations but with a more benign course.

Sensitive Monogenic Noninvasive Prenatal Diagnosis by Targeted Haplotyping.

During pregnancy, cell-free DNA (cfDNA) in maternal blood encompasses a small percentage of cell-free fetal DNA (cffDNA), an easily accessible source for determination of fetal disease status in risk families through non-invasive procedures. In case of monogenic heritable disease, background maternal cfDNA prohibits direct observation of the maternally inherited allele. Non-invasive prenatal diagnostics (NIPD) of monogenic diseases therefore relies on parental haplotyping and statistical assessment of inherited alleles from cffDNA, techniques currently unavailable for routine clinical practice. Here, we present monogenic NIPD (MG-NIPD), which requires a blood sample from both parents, for targeted locus amplification (TLA)-based phasing of heterozygous variants selectively at a gene of interest. Capture probes-based targeted sequencing of cfDNA from the pregnant mother and a tailored statistical analysis enables predicting fetal gene inheritance. MG-NIPD was validated for 18 pregnancies, focusing on CFTR, CYP21A2, and HBB. In all cases we could predict the inherited alleles with >98% confidence, even at relatively early stages (8 weeks) of pregnancy. This prediction and the accuracy of parental haplotyping was confirmed by sequencing of fetal material obtained by parallel invasive procedures. MG-NIPD is a robust method that requires standard instrumentation and can be implemented in any clinic to provide families carrying a severe monogenic disease with a prenatal diagnostic test based on a simple blood draw.

NIPTRIC: an online tool for clinical interpretation of non-invasive prenatal testing (NIPT) results.

To properly interpret the result of a pregnant woman's non-invasive prenatal test (NIPT), her a priori risk must be taken into account in order to obtain her personalised a posteriori risk (PPR), which more accurately expresses her true likelihood of carrying a foetus with trisomy. Our aim was to develop a tool for laboratories and clinicians to calculate easily the PPR for genome-wide NIPT results, using diploid samples as a control group. The tool takes the a priori risk and Z-score into account. Foetal DNA percentage and coefficient of variation can be given default settings, but actual values should be used if known. We tested the tool on 209 samples from pregnant women undergoing NIPT. For Z-scores < 5, the PPR is considerably higher at a high a priori risk than at a low a priori risk, for NIPT results with the same Z-score, foetal DNA percentage and coefficient of variation. However, the PPR is effectively independent under all conditions for Z-scores above 6. A high PPR for low a priori risks can only be reached at Z-scores > 5. Our online tool can assist clinicians in understanding NIPT results and conveying their true clinical implication to pregnant women, because the PPR is crucial for individual counselling and decision-making.

Familial Aggregation between the 14th and 21st Century and Type 2 Diabetes Risk in an Isolated Dutch Population.

INTRODUCTION: The development of type 2 diabetes results from an interaction of hereditary factors and environmental factors. This study aimed to investigate the contribution of interrelatedness to the risk of developing type 2 diabetes in an isolated Dutch population. MATERIALS AND METHODS: A genealogical database from inhabitants living on the former island Urk between the 14th and 21st century was constructed. In a case-control study, effects of interrelatedness and the risk of type 2 diabetes were estimated with Kinship Coefficients (KCs). Relative risks in first, second, and third degree relatives and spouses of inhabitants with type 2 diabetes were compared to matched controls. RESULTS: Patients with type 2 diabetes were more interrelated, expressed by a higher KC compared to controls (7.2 vs. 5.2, p=0.001). First, second and third degree relatives had an increased risk of developing type 2 diabetes. Second degree relatives had a similar risk,1.7 (1.5-2.0) as third degree relatives,1.8 (1.5-2.2). Spouses of patients with diabetes had a 3.4 (2.7-4.4) higher risk of developing type 2 diabetes. CONCLUSIONS: Interrelatedness was higher among inhabitants with type 2 diabetes compared to controls. This differences extended beyond the nuclear family, thereby supporting the hypothesis that interrelatedness contributed to the development of type 2 diabetes on Urk. However, the size of this effect was small and the patterns of risk in first, second and third degree relatives suggested that factors other than interrelatedness were the main contributors to the development of type 2 diabetes on Urk.

Gene expression analysis identifies global gene dosage sensitivity in cancer.

Many cancer-associated somatic copy number alterations (SCNAs) are known. Currently, one of the challenges is to identify the molecular downstream effects of these variants. Although several SCNAs are known to change gene expression levels, it is not clear whether each individual SCNA affects gene expression. We reanalyzed 77,840 expression profiles and observed a limited set of 'transcriptional components' that describe well-known biology, explain the vast majority of variation in gene expression and enable us to predict the biological function of genes. On correcting expression profiles for these components, we observed that the residual expression levels (in 'functional genomic mRNA' profiling) correlated strongly with copy number. DNA copy number correlated positively with expression levels for 99% of all abundantly expressed human genes, indicating global gene dosage sensitivity. By applying this method to 16,172 patient-derived tumor samples, we replicated many loci with aberrant copy numbers and identified recurrently disrupted genes in genomically unstable cancers.

Changes in winter depression phenotype correlate with white blood cell gene expression profiles: a combined metagene and gene ontology approach.

In the present study we evaluate the feasibility of gene expression in white blood cells as a peripheral marker for winter depression. Sixteen patients with winter type seasonal affective disorder were included in the study. Blood was taken by venous puncture at three time points; in winter prior and following bright light therapy and in summer. RNA was isolated, converted into cRNA, amplified and hybridized on Illumina® gene expression arrays. The raw optical array data were quantile normalized and thereafter analyzed using a metagene approach, based on previously published Affymetrix gene array data. The raw data were also subjected to a secondary analysis focusing on circadian genes and genes involved in serotonergic neurotransmission. Differences between the conditions were analyzed, using analysis of variance on the principal components of the metagene score matrix. After correction for multiple testing no statistically significant differences were found. Another approach uses the correlation between metagene factor weights and the actual expression values, averaged over conditions. When comparing the correlations of winter vs. summer and bright light therapy vs. summer significant changes for several metagenes were found. Subsequent gene ontology analyses (DAVID and GeneTrail) of 5 major metagenes suggest an interaction between brain and white blood cells. The hypothesis driven analysis with a smaller group of genes failed to demonstrate any significant effects. The results from the combined metagene and gene ontology analyses support the idea of communication between brain and white blood cells. Future studies will need a much larger sample size to obtain information at the level of single genes.

Periconceptional folic acid associated with an increased risk of oral clefts relative to non-folate related malformations in the Northern Netherlands: a population based case-control study.

Periconceptional folic acid has been associated with a reduced risk of neural tube defects, but findings on its effect in oral clefts are largely inconclusive. This case-control study assesses the effects of periconceptional folic acid on cleft risk, using complementary data from the Dutch Oral Cleft Registry and a population-based birth defects registry (Eurocat) of children and foetuses born in the Northern Netherlands between 1997 and 2009. Cases were live-born infants with non-syndromic clefts (n = 367) and controls were infants or foetuses with chromosomal/syndromal (n = 924) or non-folate related anomalies (n = 2,021). We analyzed type/timing/duration of supplement use related to traditional cleft categories as well as to their timing (early/late embryonic periods) and underlying embryological processes (fusion/differentiation defects). Consistent supplement use during the aetiologically relevant period (weeks 0-12 postconception) was associated with an increased risk of clefts (adjusted odds ratio 1.72, 95% confidence interval 1.19-2.49), especially of cleft lip/alveolus (3.16, 1.69-5.91). Further analysis systematically showed twofold to threefold increased risks for late differentiation defects-mainly clefts of the lip/alveolus-with no significant associations for early/late fusion defects. Effects were attributable to folic acid and not to other multivitamin components, and inclusion of partial use (not covering the complete aetiologically relevant period) generally weakened associations. In conclusion, this study presents several lines of evidence indicating that periconceptional folic acid in the Northern Netherlands is associated with an increased risk of clefts, in particular of cleft lip/alveolus. This association is strengthened by the specificity, consistency, systematic pattern, and duration of exposure-response relationship of our findings, underlining the need to evaluate public health strategies regarding folic acid and to further investigate potential adverse effects.

PLS3 mutations in X-linked osteoporosis with fractures.

Plastin 3 (PLS3), a protein involved in the formation of filamentous actin (F-actin) bundles, appears to be important in human bone health, on the basis of pathogenic variants in PLS3 in five families with X-linked osteoporosis and osteoporotic fractures that we report here. The bone-regulatory properties of PLS3 were supported by in vivo analyses in zebrafish. Furthermore, in an additional five families (described in less detail) referred for diagnosis or ruling out of osteogenesis imperfecta type I, a rare variant (rs140121121) in PLS3 was found. This variant was also associated with a risk of fracture among elderly heterozygous women that was two times as high as that among noncarriers, which indicates that genetic variation in PLS3 is a novel etiologic factor involved in common, multi-factorial osteoporosis.

Contribution of rare and common variants determine complex diseases-Hirschsprung disease as a model.

Finding genes for complex diseases has been the goal of many genetic studies. Most of these studies have been successful by searching for genes and mutations in rare familial cases, by screening candidate genes and by performing genome wide association studies. However, only a small fraction of the total genetic risk for these complex genetic diseases can be explained by the identified mutations and associated genetic loci. In this review we focus on Hirschsprung disease (HSCR) as an example of a complex genetic disorder. We describe the genes identified in this congenital malformation and postulate that both common 'low penetrant' variants in combination with rare or private 'high penetrant' variants determine the risk on HSCR, and likely, on other complex diseases. We also discuss how new technological advances can be used to gain further insights in the genetic background of complex diseases. Finally, we outline a few steps to develop functional assays in order to determine the involvement of these variants in disease development.

Comparison of brain and blood gene expression in an animal model of negative symptoms in schizophrenia.

OBJECTIVES: To investigate the potential of white blood cells as probes for central processes we have measured gene expression in both the anterior cingulate cortex and white blood cells using a putative animal model of negative symptoms in schizophrenia. METHODS: The model is based on the capability of ketamine to induce psychotic symptoms in healthy volunteers and to worsen such symptoms in schizophrenic patients. Classical fear conditioning is used to assess emotional processing and cognitive function in animals exposed to sub-chronic ketamine vs. controls. Gene expression was measured using a commercially sourced whole genome rat gene array. Data analyses were performed using ANOVA (Systat 11). RESULTS: In both anterior cingulate cortex and white blood cells a significant interaction between ketamine and fear conditioning could be observed. The outcome is largely supported by our subsequent metagene analysis. Moreover, the correlation between gene expression in brain and blood is about constant when no ketamine is present (r~0.4). With ketamine, however, the correlation becomes very low (r~0.2) when there is no fear, but it increases to ~0.6 when fear and ketamine are both present. Our results show that under normal conditions ketamine lowers gene expression in the brain, but this effect is completely reversed in combination with fear conditioning, indicating a stimulatory action. CONCLUSION: This paradoxical outcome indicates that extreme care must be taken when using gene expression data from white blood cells as marker for psychiatric disorders, especially when pharmacological and environmental interactions are at play.

Unraveling the regulatory mechanisms underlying tissue-dependent genetic variation of gene expression.

It is known that genetic variants can affect gene expression, but it is not yet completely clear through what mechanisms genetic variation mediate this expression. We therefore compared the cis-effect of single nucleotide polymorphisms (SNPs) on gene expression between blood samples from 1,240 human subjects and four primary non-blood tissues (liver, subcutaneous, and visceral adipose tissue and skeletal muscle) from 85 subjects. We characterized four different mechanisms for 2,072 probes that show tissue-dependent genetic regulation between blood and non-blood tissues: on average 33.2% only showed cis-regulation in non-blood tissues; 14.5% of the eQTL probes were regulated by different, independent SNPs depending on the tissue of investigation. 47.9% showed a different effect size although they were regulated by the same SNPs. Surprisingly, we observed that 4.4% were regulated by the same SNP but with opposite allelic direction. We show here that SNPs that are located in transcriptional regulatory elements are enriched for tissue-dependent regulation, including SNPs at 3' and 5' untranslated regions (P = 1.84×10(-5) and 4.7×10(-4), respectively) and SNPs that are synonymous-coding (P = 9.9×10(-4)). SNPs that are associated with complex traits more often exert a tissue-dependent effect on gene expression (P = 2.6×10(-10)). Our study yields new insights into the genetic basis of tissue-dependent expression and suggests that complex trait associated genetic variants have even more complex regulatory effects than previously anticipated.

CLMP is required for intestinal development, and loss-of-function mutations cause congenital short-bowel syndrome.

BACKGROUND & AIMS: Short-bowel syndrome usually results from surgical resection of the small intestine for diseases such as intestinal atresias, volvulus, and necrotizing enterocolitis. Patients with congenital short-bowel syndrome (CSBS) are born with a substantial shortening of the small intestine, to a mean length of 50 cm, compared with a normal length at birth of 190-280 cm. They also are born with intestinal malrotation. Because CSBS occurs in many consanguineous families, it is considered to be an autosomal-recessive disorder. We aimed to identify and characterize the genetic factor causing CSBS. METHODS: We performed homozygosity mapping using 610,000 K single-nucleotide polymorphism arrays to analyze the genomes of 5 patients with CSBS. After identifying a gene causing the disease, we determined its expression pattern in human embryos. We also overexpressed forms of the gene product that were and were not associated with CSBS in Chinese Hamster Ovary and T84 cells and generated a zebrafish model of the disease. RESULTS: We identified loss-of-function mutations in Coxsackie- and adenovirus receptor-like membrane protein (CLMP) in CSBS patients. CLMP is a tight-junction-associated protein that is expressed in the intestine of human embryos throughout development. Mutations in CLMP prevented its normal localization to the cell membrane. Knock-down experiments in zebrafish resulted in general developmental defects, including shortening of the intestine and the absence of goblet cells. Because goblet cells are characteristic for the midintestine in zebrafish, which resembles the small intestine in human beings, the zebrafish model mimics CSBS. CONCLUSIONS: Loss-of-function mutations in CLMP cause CSBS in human beings, likely by interfering with tight-junction formation, which disrupts intestinal development. Furthermore, we developed a zebrafish model of CSBS.

Trans-eQTLs reveal that independent genetic variants associated with a complex phenotype converge on intermediate genes, with a major role for the HLA.

For many complex traits, genetic variants have been found associated. However, it is still mostly unclear through which downstream mechanism these variants cause these phenotypes. Knowledge of these intermediate steps is crucial to understand pathogenesis, while also providing leads for potential pharmacological intervention. Here we relied upon natural human genetic variation to identify effects of these variants on trans-gene expression (expression quantitative trait locus mapping, eQTL) in whole peripheral blood from 1,469 unrelated individuals. We looked at 1,167 published trait- or disease-associated SNPs and observed trans-eQTL effects on 113 different genes, of which we replicated 46 in monocytes of 1,490 different individuals and 18 in a smaller dataset that comprised subcutaneous adipose, visceral adipose, liver tissue, and muscle tissue. HLA single-nucleotide polymorphisms (SNPs) were 10-fold enriched for trans-eQTLs: 48% of the trans-acting SNPs map within the HLA, including ulcerative colitis susceptibility variants that affect plausible candidate genes AOAH and TRBV18 in trans. We identified 18 pairs of unlinked SNPs associated with the same phenotype and affecting expression of the same trans-gene (21 times more than expected, P<10(-16)). This was particularly pronounced for mean platelet volume (MPV): Two independent SNPs significantly affect the well-known blood coagulation genes GP9 and F13A1 but also C19orf33, SAMD14, VCL, and GNG11. Several of these SNPs have a substantially higher effect on the downstream trans-genes than on the eventual phenotypes, supporting the concept that the effects of these SNPs on expression seems to be much less multifactorial. Therefore, these trans-eQTLs could well represent some of the intermediate genes that connect genetic variants with their eventual complex phenotypic outcomes.

MixupMapper: correcting sample mix-ups in genome-wide datasets increases power to detect small genetic effects.

MOTIVATION: Sample mix-ups can arise during sample collection, handling, genotyping or data management. It is unclear how often sample mix-ups occur in genome-wide studies, as there currently are no post hoc methods that can identify these mix-ups in unrelated samples. We have therefore developed an algorithm (MixupMapper) that can both detect and correct sample mix-ups in genome-wide studies that study gene expression levels. RESULTS: We applied MixupMapper to five publicly available human genetical genomics datasets. On average, 3% of all analyzed samples had been assigned incorrect expression phenotypes: in one of the datasets 23% of the samples had incorrect expression phenotypes. The consequences of sample mix-ups are substantial: when we corrected these sample mix-ups, we identified on average 15% more significant cis-expression quantitative trait loci (cis-eQTLs). In one dataset, we identified three times as many significant cis-eQTLs after correction. Furthermore, we show through simulations that sample mix-ups can lead to an underestimation of the explained heritability of complex traits in genome-wide association datasets. AVAILABILITY AND IMPLEMENTATION: MixupMapper is freely available at http://www.genenetwork.nl/mixupmapper/

A haplotype sharing method for determining the relative age of SNP alleles.

There are two aspects regarding the age of alleles that are relevant as indicators of the timing of mutational events. The first is to know which alleles are species-specific; the second is about the time of origin of species-specific alleles. Both aspects can be analyzed using haplotype-sharing methods, by using the length of shared haplotypes as a measure of the speed of coalescence to common ancestors. The availability of sequence data for closely related species makes it possible to infer the original SNP allele. The allele present in more than one species is the original allele. In general, original alleles are expected to be more frequent, because the cumulative effects of genetic drift determine the maximum frequency a new mutant can reach. The human species is relatively young, and founder effects are still observable as extended linkage disequilibrium. Coalescence to a single founder takes place in human populations over a time frame that is so small that original haplotypes spanning several markers are still observable in current high-density SNP genotyping arrays. We show here that the length of shared haplotypes surrounding alleles is an indicator of the relative ages of alleles, and it is applicable to original and species-specific alleles.