The effective family size of immigrant founders predicts their long-term demographic outcome: From Québec settlers to their 20th-century descendants.

Population history reconstruction, using extant genetic diversity data, routinely relies on simple demographic models to project the past through ascending genealogical-tree branches. Because genealogy and genetics are intimately related, we traced descending genealogies of the Québec founders to pursue their fate and to assess their contribution to the present-day population. Focusing on the female and male founder lines, we observed important sex-biased immigration in the early colony years and documented a remarkable impact of these early immigrants on the genetic make-up of 20th-century Québec. We estimated the immigrants' survival ratio as a proportion of lineages found in the 1931-60 Québec to their number introduced within the immigration period. We assessed the effective family size, EFS, of all immigrant parents and their Québec-born descendants. The survival ratio of the earliest immigrants was the highest and declined over centuries in association with the immigrants' EFS. Parents with high EFS left plentiful married descendants, putting EFS as the most important variable determining the parental demographic success throughout time for generations ahead. EFS of immigrant founders appears to predict their long-term demographic and, consequently, their genetic outcome. Genealogically inferred immigrants' "autosomal" genetic contribution to 1931-60 Québec from consecutive immigration periods follow the same yearly pattern as the corresponding maternal and paternal lines. Québec genealogical data offer much broader information on the ancestral diversity distribution than genetic scrutiny of a limited population sample. Genealogically inferred population history could assist studies of evolutionary factors shaping population structure and provide tools to target specific health interventions.

Historical human remains identification through maternal and paternal genetic signatures in a founder population with extensive genealogical record.

OBJECTIVES: We describe a method to identify human remains excavated from unmarked graves in historical Québec cemeteries by combining parental-lineage genetic markers with the whole-population genealogy of Québec contained in the BALSAC database. MATERIALS AND METHODS: The remains of six men were exhumed from four historical cemeteries in the province of Québec, Canada. DNA was extracted from the remains and genotyped to reveal their mitochondrial and Y-chromosome haplotypes, which were compared to a collection of haplotypes of genealogically-anchored modern volunteers. Maternal and paternal genealogies were searched in the BALSAC genealogical record for parental couples matching the mitochondrial and the Y-chromosome haplotypic signatures, to identify candidate sons from whom the remains could have originated. RESULTS: Analysis of the matching genealogies identified the parents of one man inhumed in the cemetery of the investigated parish during its operating time. The candidate individual died in 1833 at the age of 58, a plausible age at death in light of osteological analysis of the remains. DISCUSSION: This study demonstrates the promising potential of coupling genetic information from living individuals to genealogical data in BALSAC to identify historical human remains. If genetic coverage is increased, the genealogical information in BALSAC could enable the identification of 87% of the men (n = 178,435) married in Québec before 1850, with high discriminatory power in most cases since >75% of the parental couples have unique biparental signatures in most regions. Genotyping and identifying Québec's historical human remains are a key to reconstructing the genomes of the founders of Québec and reinhuming archeological remains with a marked grave.

Molecular characterization of a large group of Mucopolysaccharidosis type IIIC patients reveals the evolutionary history of the disease.

Mucopolysaccharidosis type IIIC (MPSIIIC) is a severe, rare autosomal recessive disorder caused by variants in the heparan-α-glucosaminide N-acetyltransferase (HGSNAT) gene which result in lysosomal accumulation of heparan sulfate. We analyzed clinical presentation, molecular defects and their haplotype context in 78 (27 novel) MPSIIIC cases from 22 countries, the largest group studied so far. We describe for the first time disease-causing variants in the patients from Brazil, Algeria, Azerbaijan, and Iran, and extend their spectrum within Canada, Colombia, Turkey, and the USA. Six variants are novel: two missense, c.773A>T/p.N258I and c.1267G>T/p.G423W, a nonsense c.164T>A/p.L55*, a splice-site mutation c.494-1G>A/p.[P165_L187delinsQSCYVTQAGVRWHHLGSLQALPPGFTPFSYLSLLSSWNC,P165fs], a deletion c.1348delG/p.(D450fs) and an insertion c.1479dupA/p.(Leu494fs). The missense HGSNAT variants lacked lysosomal targeting, enzymatic activity, and likely the correct folding. The haplotype analysis identified founder mutations, p.N258I, c.525dupT, and p.L55* in the Brazilian state of Paraiba, c.493+1G>A in Eastern Canada/Quebec, p.A489E in the USA, p.R384* in Poland, p.R344C and p.S518F in the Netherlands and suggested that variants c.525dupT, c.372-2G>A, and c.234+1G>A present in cis with c.564-98T>C and c.710C>A rare single-nucleotide polymorphisms, have been introduced by Portuguese settlers in Brazil. Altogether, our results provide insights into the origin, migration roots and founder effects of HGSNAT disease-causing variants, and reveal the evolutionary history of MPSIIIC.

Inferring Transmission Histories of Rare Alleles in Population-Scale Genealogies.

Learning the transmission history of alleles through a family or population plays an important role in evolutionary, demographic, and medical genetic studies. Most classical models of population genetics have attempted to do so under the assumption that the genealogy of a population is unavailable and that its idiosyncrasies can be described by a small number of parameters describing population size and mate choice dynamics. Large genetic samples have increased sensitivity to such modeling assumptions, and large-scale genealogical datasets become a useful tool to investigate realistic genealogies. However, analyses in such large datasets are often intractable using conventional methods. We present an efficient method to infer transmission paths of rare alleles through population-scale genealogies. Based on backward-time Monte Carlo simulations of genetic inheritance, we use an importance sampling scheme to dramatically speed up convergence. The approach can take advantage of available genotypes of subsets of individuals in the genealogy including haplotype structure as well as information about the mode of inheritance and general prevalence of a mutation or disease in the population. Using a high-quality genealogical dataset of more than three million married individuals in the Quebec founder population, we apply the method to reconstruct the transmission history of chronic atrial and intestinal dysrhythmia (CAID), a rare recessive disease. We identify the most likely early carriers of the mutation and geographically map the expected carrier rate in the present-day French-Canadian population of Quebec.

Genome diversity in the Neolithic Globular Amphorae culture and the spread of Indo-European languages.

It is unclear whether Indo-European languages in Europe spread from the Pontic steppes in the late Neolithic, or from Anatolia in the Early Neolithic. Under the former hypothesis, people of the Globular Amphorae culture (GAC) would be descended from Eastern ancestors, likely representing the Yamnaya culture. However, nuclear (six individuals typed for 597 573 SNPs) and mitochondrial (11 complete sequences) DNA from the GAC appear closer to those of earlier Neolithic groups than to the DNA of all other populations related to the Pontic steppe migration. Explicit comparisons of alternative demographic models via approximate Bayesian computation confirmed this pattern. These results are not in contrast to Late Neolithic gene flow from the Pontic steppes into Central Europe. However, they add nuance to this model, showing that the eastern affinities of the GAC in the archaeological record reflect cultural influences from other groups from the East, rather than the movement of people.

Mother's curse neutralizes natural selection against a human genetic disease over three centuries.

According to evolutionary theory, mitochondria could be poisoned gifts that mothers transmit to their sons. This is because mutations harmful to males are expected to accumulate in the mitochondrial genome, the so-called 'mother's curse'. However, the contribution of the mother's curse to the mutation load in nature remains largely unknown and hard to predict, because compensatory mechanisms could impede the spread of deleterious mitochondria. Here we provide evidence for the mother's curse in action over 290 years in a human population. We studied a mutation causing Leber's hereditary optical neuropathy, a disease with male-biased prevalence and which has long been suspected to be maintained in populations by the mother's curse. Male carriers showed a low fitness relative to non-carriers and to females, mostly explained by their high rate of infant mortality. Despite poor male fitness, selection analysis predicted a slight (albeit non-significant) increase in frequency, which sharply contrasts with the 35.5% per-generation decrease predicted if mitochondrial DNA transmission had been through males instead of females. Our results are therefore even suggestive of positive selection through the female line that may exacerbate effects of the mother's curse. This study supports a contribution of the mother's curse to the reduction of male lifespan, uncovering a large fitness effect associated with a single mitochondrial variant.

Genetic Structure of First Nation Communities in the Pacific Northwest.

This study presents genetic data for nine Native American populations from northern North America. Analyses of genetic variation focus on the Pacific Northwest (PNW). Using mitochondrial, Y chromosomal, and autosomal DNA variants, we aimed to more closely address the relationships of geography and language with present genetic diversity among the regional PNW Native American populations. Patterns of genetic diversity exhibited by the three genetic systems were consistent with our hypotheses: genetic variation was more strongly explained by geographic proximity than by linguistic structure. Our findings were corroborated through a variety on analytic approaches, with the unrooted trees for the three genetic systems consistently separating inland from coastal PNW populations. Furthermore, analyses of molecular variance support the trends exhibited by the unrooted trees, with geographic partitioning of PNW populations (FCT = 19.43%, p = 0.010 ± 0.009) accounting for over twice as much of the observed genetic variation as linguistic partitioning of the same populations (FCT = 9.15%, p = 0.193 ± 0.013). These findings demonstrate a consensus with previous PNW population studies examining the relationships of genome-wide variation, mitochondrial haplogroup frequencies, and skeletal morphology with geography and language.

Joubert Syndrome in French Canadians and Identification of Mutations in CEP104.

Joubert syndrome (JBTS) is a primarily autosomal-recessive disorder characterized by a distinctive mid-hindbrain and cerebellar malformation, oculomotor apraxia, irregular breathing, developmental delay, and ataxia. JBTS is a genetically heterogeneous ciliopathy. We sought to characterize the genetic landscape associated with JBTS in the French Canadian (FC) population. We studied 43 FC JBTS subjects from 35 families by combining targeted and exome sequencing. We identified pathogenic (n = 32 families) or possibly pathogenic (n = 2 families) variants in genes previously associated with JBTS in all of these subjects, except for one. In the latter case, we found a homozygous splice-site mutation (c.735+2T>C) in CEP104. Interestingly, we identified two additional non-FC JBTS subjects with mutations in CEP104; one of these subjects harbors a maternally inherited nonsense mutation (c.496C>T [p.Arg166*]) and a de novo splice-site mutation (c.2572-2A>G), whereas the other bears a homozygous frameshift mutation (c.1328_1329insT [p.Tyr444fs*3]) in CEP104. Previous studies have shown that CEP104 moves from the mother centriole to the tip of the primary cilium during ciliogenesis. Knockdown of CEP104 in retinal pigment epithelial (RPE1) cells resulted in severe defects in ciliogenesis. These observations suggest that CEP104 acts early during cilia formation by regulating the conversion of the mother centriole into the cilia basal body. We conclude that disruption of CEP104 causes JBTS. Our study also reveals that the cause of JBTS has been elucidated in the great majority of our FC subjects (33/35 [94%] families), even though JBTS shows substantial locus and allelic heterogeneity in this population.

Global diversity, population stratification, and selection of human copy-number variation.

In order to explore the diversity and selective signatures of duplication and deletion human copy-number variants (CNVs), we sequenced 236 individuals from 125 distinct human populations. We observed that duplications exhibit fundamentally different population genetic and selective signatures than deletions and are more likely to be stratified between human populations. Through reconstruction of the ancestral human genome, we identify megabases of DNA lost in different human lineages and pinpoint large duplications that introgressed from the extinct Denisova lineage now found at high frequency exclusively in Oceanic populations. We find that the proportion of CNV base pairs to single-nucleotide-variant base pairs is greater among non-Africans than it is among African populations, but we conclude that this difference is likely due to unique aspects of non-African population history as opposed to differences in CNV load.

GENLIB: an R package for the analysis of genealogical data.

BACKGROUND: Founder populations have an important role in the study of genetic diseases. Access to detailed genealogical records is often one of their advantages. These genealogical data provide unique information for researchers in evolutionary and population genetics, demography and genetic epidemiology. However, analyzing large genealogical datasets requires specialized methods and software. The GENLIB software was developed to study the large genealogies of the French Canadian population of Quebec, Canada. These genealogies are accessible through the BALSAC database, which contains over 3 million records covering the whole province of Quebec over four centuries. Using this resource, extended pedigrees of up to 17 generations can be constructed from a sample of present-day individuals. RESULTS: We have extended and implemented GENLIB as a package in the R environment for statistical computing and graphics, thus allowing optimal flexibility for users. The GENLIB package includes basic functions to manage genealogical data allowing, for example, extraction of a part of a genealogy or selection of specific individuals. There are also many functions providing information to describe the size and complexity of genealogies as well as functions to compute standard measures such as kinship, inbreeding and genetic contribution. GENLIB also includes functions for gene-dropping simulations. The goal of this paper is to present the full functionalities of GENLIB. We used a sample of 140 individuals from the province of Quebec (Canada) to demonstrate GENLIB's functions. Ascending genealogies for these individuals were reconstructed using BALSAC, yielding a large pedigree of 41,523 individuals. Using GENLIB's functions, we provide a detailed description of these genealogical data in terms of completeness, genetic contribution of founders, relatedness, inbreeding and the overall complexity of the genealogical tree. We also present gene-dropping simulations based on the whole genealogy to investigate identical-by-descent sharing of alleles and chromosomal segments of different lengths and estimate probabilities of identical-by-descent sharing. CONCLUSIONS: The R package GENLIB provides a user friendly and flexible environment to analyze extensive genealogical data, allowing an efficient and easy integration of different types of data, analytical methods and additional developments and making this tool ideal for genealogical analysis.

Disruption of CLPB is associated with congenital microcephaly, severe encephalopathy and 3-methylglutaconic aciduria.

BACKGROUND: The heterogeneous group of 3-methylglutaconic aciduria disorders includes several inborn errors of metabolism that affect mitochondrial function through poorly understood mechanisms. We describe four newborn siblings, from a consanguineous family, who showed microcephaly, small birth weight, severe encephalopathy and 3-methylglutaconic aciduria. Their neurological examination was characterised by severe hypertonia and the induction of prolonged clonic movements of the four limbs upon minimal tactile stimulation. METHODS AND RESULTS: Using homozygosity mapping and exome sequencing, we identified a homozygous truncating mutation (p.I562Tfs*23) in CLPB segregating with the disease in this family. CLPB codes for a member of the family of ATPases associated with various cellular activities (AAA(+) proteins) whose function remains unknown. We found that CLPB expression is abolished in fibroblasts from the patients. To investigate the function of this gene, we interfered with the translation of the zebrafish clpb orthologue using an antisense morpholino. The clpb morphants showed an abnormal touch-evoked response with increased swim velocity and tail beat frequency. This motor phenotype is reminiscent of that observed in the patients and is suggestive of increased excitability in neuronal circuits. Interestingly, knocking down clpb reduced the number of inhibitory glycinergic interneurons and increased a population of excitatory glutamatergic neurons in the spinal cord. CONCLUSIONS: Altogether, our study suggests that disruption of CLPB causes a novel form of neonatal encephalopathy associated with 3-methylglutaconic aciduria.

Nonsense mutations in BAG3 are associated with early-onset dilated cardiomyopathy in French Canadians.

BACKGROUND: Dilated cardiomyopathy (DCM) is a major cause of heart failure that may require heart transplantation. Approximately one third of DCM cases are familial. Next-generation DNA sequencing of large panels of candidate genes (ie, targeted sequencing) or of the whole exome can rapidly and economically identify pathogenic mutations in familial DCM. METHODS: We recruited 64 individuals from 26 DCM families followed at the Montreal Heart Institute Cardiovascular Genetic Center and sequenced the whole exome of 44 patients and 2 controls. Both affected and unaffected family members underwent genotyping for segregation analysis. RESULTS: We found 2 truncating mutations in BAG3 in 4 DCM families (15%) and confirmed segregation with disease status by linkage (log of the odds [LOD] score = 3.8). BAG3 nonsense mutations conferred a worse prognosis as evidenced by a younger age of clinical onset (37 vs 48 years for carriers and noncarriers respectively; P = 0.037). We also found truncating mutations in TTN in 5 families (19%). Finally, we identified potential pathogenic mutations for 9 DCM families in 6 candidate genes (DSP, LMNA, MYH7, MYPN, RBM20, and TNNT2). We still need to confirm several of these mutations by segregation analysis. CONCLUSIONS: Screening an extended panel of 41 candidate genes allowed us to identify probable pathogenic mutations in 69% of families with DCM in our cohort of mostly French-Canadian patients. We confirmed the prevalence of TTN nonsense mutations in DCM. Furthermore, to our knowledge, we are the first to present an association between nonsense mutations in BAG3 and early-onset DCM.

Mutations in SGOL1 cause a novel cohesinopathy affecting heart and gut rhythm.

The pacemaking activity of specialized tissues in the heart and gut results in lifelong rhythmic contractions. Here we describe a new syndrome characterized by Chronic Atrial and Intestinal Dysrhythmia, termed CAID syndrome, in 16 French Canadians and 1 Swede. We show that a single shared homozygous founder mutation in SGOL1, a component of the cohesin complex, causes CAID syndrome. Cultured dermal fibroblasts from affected individuals showed accelerated cell cycle progression, a higher rate of senescence and enhanced activation of TGF-ß signaling. Karyotypes showed the typical railroad appearance of a centromeric cohesion defect. Tissues derived from affected individuals displayed pathological changes in both the enteric nervous system and smooth muscle. Morpholino-induced knockdown of sgol1 in zebrafish recapitulated the abnormalities seen in humans with CAID syndrome. Our findings identify CAID syndrome as a novel generalized dysrhythmia, suggesting a new role for SGOL1 and the cohesin complex in mediating the integrity of human cardiac and gut rhythm.

The catalytic function of hormone-sensitive lipase is essential for fertility in male mice.

In male mice, deficiency of hormone sensitive lipase (HSL, Lipe gene, E.C.3.1.1.3) causes deficient spermatogenesis, azoospermia, and infertility. Postmeiotic germ cells express a specific HSL isoform that includes a 313 amino acid N-terminus encoded by a testis-specific exon (exon T1). The remainder of testicular HSL is identical to adipocyte HSL. The amino acid sequence of the testis-specific exon is poorly conserved, showing only a 46% amino acid identity with orthologous human and rat sequences, compared with 87% over the remainder of the HSL coding sequence, providing no evidence in favor of a vital functional role for the testis-specific N-terminus of HSL. However, exon T1 is important for Lipe transcription; in mouse testicular mRNA, we identified 3 major Lipe transcription start sites, finding numerous testicular transcription factor binding motifs upstream of the transcription start site. We directly explored two possible mechanisms for the infertility of HSL-deficient mice, using mice that expressed mutant HSL transgenes only in postmeiotic germ cells on a HSL-deficient background. One transgene expressed human HSL lacking enzyme activity but containing the testis-specific N-terminus (HSL-/-muttg mice). The other transgene expressed catalytically inactive HSL with the testis-specific N-terminal peptide (HSL-/-atg mice). HSL-/-muttg mice were infertile, with abnormal histology of the seminiferous epithelium and absence of spermatozoa in the epididymal lumen. In contrast, HSL-/-atg mice had normal fertility and normal testicular morphology. In conclusion, whereas the catalytic function of HSL is necessary for spermatogenesis in mice, the presence of the N-terminal testis-specific fragment is not essential.

Genome-wide patterns of identity-by-descent sharing in the French Canadian founder population.

In genetics the ability to accurately describe the familial relationships among a group of individuals can be very useful. Recent statistical tools succeeded in assessing the degree of relatedness up to 6-7 generations with good power using dense genome-wide single-nucleotide polymorphism data to estimate the extent of identity-by-descent (IBD) sharing. It is therefore important to describe genome-wide patterns of IBD sharing for more remote and complex relatedness between individuals, such as that observed in a founder population like Quebec, Canada. Taking advantage of the extended genealogical records of the French Canadian founder population, we first compared different tools to identify regions of IBD in order to best describe genome-wide IBD sharing and its correlation with genealogical characteristics. Results showed that the extent of IBD sharing identified with FastIBD correlates best with relatedness measured using genealogical data. Total length of IBD sharing explained 85% of the genealogical kinship's variance. In addition, we observed significantly higher sharing in pairs of individuals with at least one inbred ancestor compared with those without any. Furthermore, patterns of IBD sharing and average sharing were different across regional populations, consistent with the settlement history of Quebec. Our results suggest that, as expected, the complex relatedness present in founder populations is reflected in patterns of IBD sharing. Using these patterns, it is thus possible to gain insight on the types of distant relationships in a sample from a founder population like Quebec.

X-linked MTMR8 diversity and evolutionary history of sub-Saharan populations.

The genetic diversity within an 11 kb segment of the MTMR8 gene in a sample of 111 sub-Saharan and 49 non-African X chromosomes was investigated to assess the early evolutionary history of sub-Saharan Africans and the out-of-Africa expansion. The analyses revealed a complex genetic structure of the Africans that contributed to the emergence of modern humans. We observed partitioning of two thirds of old lineages among southern, west/central and east African populations indicating ancient population stratification predating the out of Africa migration. Age estimates of these lineages, older than coalescence times of uniparentally inherited markers, raise the question whether contemporary humans originated from a single population or as an amalgamation of different populations separated by years of independent evolution, thus suggesting a greater antiquity of our species than generally assumed. While the oldest sub-Saharan lineages, ~500 thousand years, are found among Khoe-San from southern-Africa, a distinct haplotype found among Biaka is likely due to admixture from an even older population. An East African population that gave rise to non-Africans underwent a selective sweep affecting the subcentromeric region where MTMR8 is located. This and similar sweeps in four other regions of the X chromosome, documented in the literature, effectively reduced genetic diversity of non-African chromosomes and therefore may have exacerbated the effect of the demographic bottleneck usually ascribed to the out of Africa migration. Our data is suggestive, however, that a bottleneck, occurred in Africa before range expansion.

Reconciling migration models to the Americas with the variation of North American native mitogenomes.

In this study we evaluated migration models to the Americas by using the information contained in native mitochondrial genomes (mitogenomes) from North America. Molecular and phylogeographic analyses of B2a mitogenomes, which are absent in Eskimo-Aleut and northern Na-Dene speakers, revealed that this haplogroup arose in North America ∼11-13 ka from one of the founder Paleo-Indian B2 mitogenomes. In contrast, haplogroup A2a, which is typical of Eskimo-Aleuts and Na-Dene, but also present in the easternmost Siberian groups, originated only 4-7 ka in Alaska, led to the first Paleo-Eskimo settlement of northern Canada and Greenland, and contributed to the formation of the Na-Dene gene pool. However, mitogenomes also show that Amerindians from northern North America, without any distinction between Na-Dene and non-Na-Dene, were heavily affected by an additional and distinctive Beringian genetic input. In conclusion, most mtDNA variation (along the double-continent) stems from the first wave from Beringia, which followed the Pacific coastal route. This was accompanied or followed by a second inland migratory event, marked by haplogroups X2a and C4c, which affected all Amerindian groups of Northern North America. Much later, the ancestral A2a carriers spread from Alaska, undertaking both a westward migration to Asia and an eastward expansion into the circumpolar regions of Canada. Thus, the first American founders left the greatest genetic mark but the original maternal makeup of North American Natives was subsequently reshaped by additional streams of gene flow and local population dynamics, making a three-wave view too simplistic.

Native American admixture in the Quebec founder population.

For years, studies of founder populations and genetic isolates represented the mainstream of genetic mapping in the effort to target genetic defects causing Mendelian disorders. The genetic homogeneity of such populations as well as relatively homogeneous environmental exposures were also seen as primary advantages in studies of genetic susceptibility loci that underlie complex diseases. European colonization of the St-Lawrence Valley by a small number of settlers, mainly from France, resulted in a founder effect reflected by the appearance of a number of population-specific disease-causing mutations in Quebec. The purported genetic homogeneity of this population was recently challenged by genealogical and genetic analyses. We studied one of the contributing factors to genetic heterogeneity, early Native American admixture that was never investigated in this population before. Consistent admixture estimates, in the order of one per cent, were obtained from genome-wide autosomal data using the ADMIXTURE and HAPMIX software, as well as with the fastIBD software evaluating the degree of the identity-by-descent between Quebec individuals and Native American populations. These genomic results correlated well with the genealogical estimates. Correlations are imperfect most likely because of incomplete records of Native founders' origin in genealogical data. Although the overall degree of admixture is modest, it contributed to the enrichment of the population diversity and to its demographic stratification. Because admixture greatly varies among regions of Quebec and among individuals, it could have significantly affected the homogeneity of the population, which is of importance in mapping studies, especially when rare genetic susceptibility variants are in play.

Mutations in TMEM231 cause Joubert syndrome in French Canadians.

BACKGROUND: Joubert syndrome (JBTS) is a predominantly autosomal recessive disorder characterised by a distinctive midhindbrain malformation, oculomotor apraxia, breathing abnormalities and developmental delay. JBTS is genetically heterogeneous, involving genes required for formation and function of non-motile cilia. Here we investigate the genetic basis of JBTS in 12 French-Canadian (FC) individuals. METHODS AND RESULTS: Exome sequencing in all subjects showed that six of them carried rare compound heterozygous mutations in CC2D2A or C5ORF42, known JBTS genes. In addition, three individuals (two families) were compound heterozygous for the same rare mutations in TMEM231(c.12T>A[p.Tyr4*]; c.625G>A[p.Asp209Asn]). All three subjects showed a severe neurological phenotype and variable presence of polydactyly, retinopathy and renal cysts. These mutations were not detected among 385 FC controls. TMEM231 has been previously shown to localise to the ciliary transition zone, and to interact with several JBTS gene products in a complex involved in the formation of the diffusion barrier between the cilia and plasma membrane. siRNA knockdown of TMEM231 was also shown to affect barrier integrity, resulting in a reduction of cilia formation and ciliary localisation of signalling receptors. CONCLUSIONS: Our data suggest that mutations in TMEM231 cause JBTS, reinforcing the relationship between this condition and the disruption of the barrier at the ciliary transition zone.

Genotype-based test in mapping cis-regulatory variants from allele-specific expression data.

Identifying and understanding the impact of gene regulatory variation is of considerable importance in evolutionary and medical genetics; such variants are thought to be responsible for human-specific adaptation and to have an important role in genetic disease. Regulatory variation in cis is readily detected in individuals showing uneven expression of a transcript from its two allelic copies, an observation referred to as allelic imbalance (AI). Identifying individuals exhibiting AI allows mapping of regulatory DNA regions and the potential to identify the underlying causal genetic variant(s). However, existing mapping methods require knowledge of the haplotypes, which make them sensitive to phasing errors. In this study, we introduce a genotype-based mapping test that does not require haplotype-phase inference to locate regulatory regions. The test relies on partitioning genotypes of individuals exhibiting AI and those not expressing AI in a 2×3 contingency table. The performance of this test to detect linkage disequilibrium (LD) between a potential regulatory site and a SNP located in this region was examined by analyzing the simulated and the empirical AI datasets. In simulation experiments, the genotype-based test outperforms the haplotype-based tests with the increasing distance separating the regulatory region from its regulated transcript. The genotype-based test performed equally well with the experimental AI datasets, either from genome-wide cDNA hybridization arrays or from RNA sequencing. By avoiding the need of haplotype inference, the genotype-based test will suit AI analyses in population samples of unknown haplotype structure and will additionally facilitate the identification of cis-regulatory variants that are located far away from the regulated transcript.