Accounting for long-range correlations in genome-wide simulations of large cohorts.

Coalescent simulations are widely used to examine the effects of evolution and demographic history on the genetic makeup of populations. Thanks to recent progress in algorithms and data structures, simulators such as the widely-used msprime now provide genome-wide simulations for millions of individuals. However, this software relies on classic coalescent theory and its assumptions that sample sizes are small and that the region being simulated is short. Here we show that coalescent simulations of long regions of the genome exhibit large biases in identity-by-descent (IBD), long-range linkage disequilibrium (LD), and ancestry patterns, particularly when the sample size is large. We present a Wright-Fisher extension to msprime, and show that it produces more realistic distributions of IBD, LD, and ancestry proportions, while also addressing more subtle biases of the coalescent. Further, these extensions are more computationally efficient than state-of-the-art coalescent simulations when simulating long regions, including whole-genome data. For shorter regions, efficiency can be maintained via a hybrid model which simulates the recent past under the Wright-Fisher model and uses coalescent simulations in the distant past.

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.

Testing association of rare genetic variants with resistance to three common antiseizure medications.

OBJECTIVE: Drug resistance is a major concern in the treatment of individuals with epilepsy. No genetic markers for resistance to individual antiseizure medication (ASM) have yet been identified. We aimed to identify the role of rare genetic variants in drug resistance for three common ASMs: levetiracetam (LEV), lamotrigine (LTG), and valproic acid (VPA). METHODS: A cohort of 1622 individuals of European descent with epilepsy was deeply phenotyped and underwent whole exome sequencing (WES), comprising 575 taking LEV, 826 LTG, and 782 VPA. We performed gene- and gene set-based collapsing analyses comparing responders and nonresponders to the three drugs to determine the burden of different categories of rare genetic variants. RESULTS: We observed a marginally significant enrichment of rare missense, truncating, and splice region variants in individuals who were resistant to VPA compared to VPA responders for genes involved in VPA pharmacokinetics. We also found a borderline significant enrichment of truncating and splice region variants in the synaptic vesicle glycoprotein (SV2) gene family in nonresponders compared to responders to LEV. We did not see any significant enrichment using a gene-based approach. SIGNIFICANCE: In our pharmacogenetic study, we identified a slightly increased burden of damaging variants in gene groups related to drug kinetics or targeting in individuals presenting with drug resistance to VPA or LEV. Such variants could thus determine a genetic contribution to drug resistance.

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.

Reconstructing Native American population history.

The peopling of the Americas has been the subject of extensive genetic, archaeological and linguistic research; however, central questions remain unresolved. One contentious issue is whether the settlement occurred by means of a single migration or multiple streams of migration from Siberia. The pattern of dispersals within the Americas is also poorly understood. To address these questions at a higher resolution than was previously possible, we assembled data from 52 Native American and 17 Siberian groups genotyped at 364,470 single nucleotide polymorphisms. Here we show that Native Americans descend from at least three streams of Asian gene flow. Most descend entirely from a single ancestral population that we call 'First American'. However, speakers of Eskimo-Aleut languages from the Arctic inherit almost half their ancestry from a second stream of Asian gene flow, and the Na-Dene-speaking Chipewyan from Canada inherit roughly one-tenth of their ancestry from a third stream. We show that the initial peopling followed a southward expansion facilitated by the coast, with sequential population splits and little gene flow after divergence, especially in South America. A major exception is in Chibchan speakers on both sides of the Panama isthmus, who have ancestry from both North and South America.

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.

Deciphering the genetic structure of the Quebec founder population using genealogies.

Using genealogy to study the demographic history of a population makes it possible to overcome the models and assumptions often used in population genetics. The Quebec founder population is one of the few populations in the world having access to the complete genealogy of the last 400 years. The goal of this study is to follow the evolution of the Quebec population structure over time from the beginning of European colonization until the present day. To do so, we calculated the kinship coefficients of all ancestors' pairs in the ascending genealogy of 665 subjects from eight regional and ethnocultural groups per 25-year period. We show that the Quebec population structure appeared progressively in the St. Lawrence valley as early as 1750 with the distinction of the Saguenay and Gaspesian groups. At that time, the ancestors of two groups, the Sagueneans and the Acadians from the Gaspé Peninsula, experienced a marked increase in kinship and inbreeding levels which have shaped the structure and led to the contemporary population structure. Interestingly, this structure arose before the colonization of the Saguenay region and at the very beginning of the Gaspé Peninsula settlement. The resulting regional founder effects in these groups led to differences in the present-day identity-by-descent sharing, the Gaspé and North Shore groups sharing more large segments and the Sagueneans more short segments. This is also reflected by the distribution of the number of most recent common ancestors at different generations and their genetic contribution to the studied subjects.

Unraveling the role of non-coding rare variants in epilepsy.

The discovery of new variants has leveled off in recent years in epilepsy studies, despite the use of very large cohorts. Consequently, most of the heritability is still unexplained. Rare non-coding variants have been largely ignored in studies on epilepsy, although non-coding single nucleotide variants can have a significant impact on gene expression. We had access to whole genome sequencing (WGS) from 247 epilepsy patients and 377 controls. To assess the functional impact of non-coding variants, ExPecto, a deep learning algorithm was used to predict expression change in brain tissues. We compared the burden of rare non-coding deleterious variants between cases and controls. Rare non-coding highly deleterious variants were significantly enriched in Genetic Generalized Epilepsy (GGE), but not in Non-Acquired Focal Epilepsy (NAFE) or all epilepsy cases when compared with controls. In this study we showed that rare non-coding deleterious variants are associated with epilepsy, specifically with GGE. Larger WGS epilepsy cohort will be needed to investigate those effects at a greater resolution. Nevertheless, we demonstrated the importance of studying non-coding regions in epilepsy, a disease where new discoveries are scarce.

An X-linked haplotype of Neandertal origin is present among all non-African populations.

Recent work on the Neandertal genome has raised the possibility of admixture between Neandertals and the expanding population of Homo sapiens who left Africa between 80 and 50 Kya (thousand years ago) to colonize the rest of the world. Here, we provide evidence of a notable presence (9% overall) of a Neandertal-derived X chromosome segment among all contemporary human populations outside Africa. Our analysis of 6,092 X-chromosomes from all inhabited continents supports earlier contentions that a mosaic of lineages of different time depths and different geographic provenance could have contributed to the genetic constitution of modern humans. It indicates a very early admixture between expanding African migrants and Neandertals prior to or very early on the route of the out-of-Africa expansion that led to the successful colonization of the planet.

Assessment of burden and segregation profiles of CNVs in patients with epilepsy.

OBJECTIVE: Microdeletions are associated with different forms of epilepsy but show incomplete penetrance, which is not well understood. We aimed to assess whether unmasked variants or double CNVs could explain incomplete penetrance. METHODS: We analyzed copy number variants (CNVs) in 603 patients with four different subgroups of epilepsy and 945 controls. CNVs were called from genotypes and validated on whole-genome (WGS) or whole-exome sequences (WES). CNV burden difference between patients and controls was obtained by fitting a logistic regression. CNV burden was assessed for small and large (>1 Mb) deletions and duplications and for deletions overlapping different gene sets. RESULTS: Large deletions were enriched in genetic generalized epilepsies (GGE) compared to controls. We also found enrichment of deletions in epilepsy genes and hotspots for GGE. We did not find truncating or functional variants that could have been unmasked by the deletions. We observed a double CNV hit in two patients. One patient also carried a de novo deletion in the 22q11.2 hotspot. INTERPRETATION: We could corroborate previous findings of an enrichment of large microdeletions and deletions in epilepsy genes in GGE. We could also replicate that microdeletions show incomplete penetrance. However, we could not validate the hypothesis of unmasked variants nor the hypothesis of double CNVs to explain the incomplete penetrance. We found a de novo CNV on 22q11.2 that could be of interest. We also observed GGE families carrying a deletion on 15q13.3 hotspot that could be investigated in the Quebec founder population.

The role of common genetic variation in presumed monogenic epilepsies.

BACKGROUND: The developmental and epileptic encephalopathies (DEEs) are the most severe group of epilepsies which co-present with developmental delay and intellectual disability (ID). DEEs usually occur in people without a family history of epilepsy and have emerged as primarily monogenic, with damaging rare mutations found in 50% of patients. Little is known about the genetic architecture of patients with DEEs in whom no pathogenic variant is identified. Polygenic risk scoring (PRS) is a method that measures a person's common genetic burden for a trait or condition. Here, we used PRS to test whether genetic burden for epilepsy is relevant in individuals with DEEs, and other forms of epilepsy with ID. METHODS: Genetic data on 2,759 cases with DEEs, or epilepsy with ID presumed to have a monogenic basis, and 447,760 population-matched controls were analysed. We compared PRS for 'all epilepsy', 'focal epilepsy', and 'genetic generalised epilepsy' (GGE) between cases and controls. We performed pairwise comparisons between cases stratified for identifiable rare deleterious genetic variants and controls. FINDINGS: Cases of presumed monogenic severe epilepsy had an increased PRS for 'all epilepsy' (p<0.0001), 'focal epilepsy' (p<0.0001), and 'GGE' (p=0.0002) relative to controls, which explain between 0.08% and 3.3% of phenotypic variance. PRS was increased in cases both with and without an identified deleterious variant of major effect, and there was no significant difference in PRS between the two groups. INTERPRETATION: We provide evidence that common genetic variation contributes to the aetiology of DEEs and other forms of epilepsy with ID, even when there is a known pathogenic variant of major effect. These results provide insight into the genetic underpinnings of the severe epilepsies and warrant a shift in our understanding of the aetiology of the DEEs as complex, rather than monogenic, disorders. FUNDING: Science foundation Ireland, Human Genome Research Institute; National Heart, Lung, and Blood Institute; German Research Foundation.

When genetics and genealogies tell different stories-maternal lineages in Gaspesia.

Data from uniparentally inherited genetic systems were used to trace evolution of human populations. Reconstruction of the past primarily relies on variation in present-day populations, limiting historical inference to lineages that are found among living subjects. Our analysis of four population groups in the Gaspé Peninsula, demonstrates how this may occasionally lead to erroneous interpretations. Mitochondrial DNA analysis of Gaspesians revealed an important admixture with Native Americans. The most likely scenario links this admixture to French-Canadians from the St. Lawrence Valley who moved to Gaspesia in the 19th century. However, in contrast to genetic data, analysis of genealogical record shows that Native American maternal lineages were brought to Gaspesia in the 18th century by Acadians who settled on the south-western coast of the peninsula. Intriguingly, within three generations, virtually all Métis Acadian families separated from their nonadmixed relatives and moved eastward mixing in with other Gaspesian groups, in which Native American maternal lines are present in relatively high frequencies. Over time, the carriers of these lines eventually lost memory of their mixed Amerindian-Acadian origin. Our results show that a reliable reconstruction of population history requires cross-verification of different data sources for consistency, thus favouring multidisciplinary approaches.

Genomic and genealogical investigation of the French Canadian founder population structure.

Characterizing the genetic structure of worldwide populations is important for understanding human history and is essential to the design and analysis of genetic epidemiological studies. In this study, we examined genetic structure and distant relatedness and their effect on the extent of linkage disequilibrium (LD) and homozygosity in the founder population of Quebec (Canada). In the French Canadian founder population, such analysis can be performed using both genomic and genealogical data. We investigated genetic differences, extent of LD, and homozygosity in 140 individuals from seven sub-populations of Quebec characterized by different demographic histories reflecting complex founder events. Genetic findings from genome-wide single nucleotide polymorphism data were correlated with genealogical information on each of these sub-populations. Our genomic data showed significant population structure and relatedness present in the contemporary Quebec population, also reflected in LD and homozygosity levels. Our extended genealogical data corroborated these findings and indicated that this structure is consistent with the settlement patterns involving several founder events. This provides an independent and complementary validation of genomic-based studies of population structure. Combined genomic and genealogical data in the Quebec founder population provide insights into the effects of the interplay of two important sources of bias in genetic epidemiological studies, unrecognized genetic structure and cryptic relatedness.

Assessing the role of rare genetic variants in drug-resistant, non-lesional focal epilepsy.

OBJECTIVE: Resistance to antiseizure medications (ASMs) is one of the major concerns in the treatment of epilepsy. Despite the increasing number of ASMs available, the proportion of individuals with drug-resistant epilepsy remains unchanged. In this study, we aimed to investigate the role of rare genetic variants in ASM resistance. METHODS: We performed exome sequencing of 1,128 individuals with non-familial non-acquired focal epilepsy (NAFE) (762 non-responders, 366 responders) and were provided with 1,734 healthy controls. We undertook replication in a cohort of 350 individuals with NAFE (165 non-responders, 185 responders). We performed gene-based and gene-set-based kernel association tests to investigate potential enrichment of rare variants in relation to drug response status and to risk for NAFE. RESULTS: We found no gene or gene set that reached genome-wide significance. Yet, we identified several prospective candidate genes - among them DEPDC5, which showed a potential association with resistance to ASMs. We found some evidence for an enrichment of truncating variants in dominant familial NAFE genes in our cohort of non-familial NAFE and in association with drug-resistant NAFE. INTERPRETATION: Our study identifies potential candidate genes for ASM resistance. Our results corroborate the role of rare variants for non-familial NAFE and imply their involvement in drug-resistant epilepsy. Future large-scale genetic research studies are needed to substantiate these findings.

Polygenic risk scores of several subtypes of epilepsies in a founder population.

OBJECTIVE: Polygenic risk scores (PRSs) are used to quantify the cumulative effects of a number of genetic variants, which may individually have a very small effect on susceptibility to a disease; we used PRSs to better understand the genetic contribution to common epilepsy and its subtypes. METHODS: We first replicated previous single associations using 373 unrelated patients. We then calculated PRSs in the same French Canadian patients with epilepsy divided into 7 epilepsy subtypes and population-based controls. We fitted a logistic mixed model to calculate the variance explained by the PRS using pseudo-R2 statistics. RESULTS: We show that the PRS explains more of the variance in idiopathic generalized epilepsy than in patients with nonacquired focal epilepsy. We also demonstrate that the variance explained is different within each epilepsy subtype. CONCLUSIONS: Globally, we support the notion that PRSs provide a reliable measure to rightfully estimate the contribution of genetic factors to the pathophysiologic mechanism of epilepsies, but further studies are needed on PRSs before they can be used clinically.

Genetic heterogeneity in regional populations of Quebec--parental lineages in the Gaspe Peninsula.

Stable colonization of the Gaspe Peninsula by Europeans started in the middle of the 18th century at the time of the British conquest of New France. The earliest settlers were Acadians, escaping British deportation policies, followed by Loyalists from the US, who preferred to remain under British rule after the Declaration of Independence. In the 19th century, the developing fishing industry attracted French Canadians from the St. Lawrence Valley and newcomers from Europe including Channel Islanders from Jersey and Guernsey. We analyzed parental lineages of the self-declared descendants of these four groups of settlers by mtDNA D-loop sequencing and Y-chromosome genotyping and compared them with French, British, and Irish samples. Their representation in terms of haplotype frequency classes reveals different signatures of founder effects, such as a loss of rare haplotypes, modification of intermediate frequency haplotypes, reduction in genetic diversity (seen in Acadians), but also enrichment by admixture. Parental lineages correlate with group identity. Descendants of early settlers, Acadians and Loyalists, preserved their identity more than those of French Canadian and Channel Islander "latecomers." Although overall genetic diversity among Gaspesians is comparable with their European source populations, F(ST) analysis indicated their greater differentiation. Distinct settlement history, a limited number of founders and relative genetic isolation contributed to the regionalization of the Quebec gene pool that appears less homogenous than usually anticipated.

X-chromosome lineages and the settlement of the Americas.

Most genetic studies on the origins of Native Americans have examined data from mtDNA and Y-chromosome DNA. To complement these studies and to broaden our understanding of the origin of Native American populations, we present an analysis of 1,873 X-chromosomes representing Native American (n = 438) and other continental populations (n = 1,435). We genotyped 36 polymorphic sites, forming an informative haplotype within an 8-kb DNA segment spanning exon 44 of the dystrophin gene. The data reveal continuity from a common Eurasian ancestry between Europeans, Siberians, and Native Americans. However, the loss of two haplotypes frequent in Eurasia (18.8 and 7%) and the rise in frequency of a third haplotype rare elsewhere, indicate a major population bottleneck in the peopling of the Americas. Although genetic drift appears to have played a greater role in the genetic differentiation of Native Americans than in the latitudinally distributed Eurasians, we also observe a signal of a differentiated ancestry of southern and northern populations that cannot be simply explained by the serial southward dilution of genetic diversity. It is possible that the distribution of X-chromosome lineages reflects the genetic structure of the population of Beringia, itself issued from founder effects and a source of subsequent southern colonization(s).

Genomic inference using diffusion models and the allele frequency spectrum.

Evolutionary, biological, and demographic processes together shape observed variation in populations. Understanding how these processes influence variation allows us to infer past demography and the nature of selection in populations. Forward in time models such as the diffusion approximation provide a powerful tool for performing inference based on the distribution of allele frequencies. Here, we discuss recent computational developments and their application to reconstructing human demographic history. Using whole-genome sequence data for 797 French Canadian individuals, we assess the neutrality of synonymous variants and show that selection can bias inferred demography, mutation rates, and distributions of fitness effects. We argue that the simple evolutionary models investigated by Kimura and Ohta still provide important insight into modern genetic research.

Patterns of variation in DNA segments upstream of transcription start sites.

It is likely that evolutionary differences among species are driven by sequence changes in regulatory regions. Likewise, polymorphisms in the promoter regions may be responsible for interindividual differences at the level of populations. We present an unbiased survey of genetic variation in 2-kb segments upstream of the transcription start sites of 28 protein-coding genes, characterized in five population groups of different geographic origin. On average, we found 9.1 polymorphisms and 8.8 haplotypes per segment with corresponding nucleotide and haplotype diversities of 0.082% and 58%, respectively. We characterized these segments through different summary statistics, Hardy-Weinberg equilibria fixation index (Fst) estimates, and neutrality tests, as well as by analyzing the distributions of haplotype allelic classes, introduced here to assess the departure from neutrality and examined by coalescent simulations under a simple population model, assuming recombinations or different demography. Our results suggest that genetic diversity in some of these regions could have been shaped by purifying selection and driven by adaptive changes in the other, thus explaining the relatively large variance in the corresponding genetic diversity indices loci. However, some of these effects could be also due to linkage with surrounding sequences, and the neutralists' explanations cannot be ruled out given uncertainty in the underlying demographic histories and the possibility of random effects due to the small size of the studied segments.

[Founder effects and genetic variability in Quebec]. / Effets fondateurs et variabilité génétique au Québec.

Knowledge of the genetic population structure lies at the heart of mapping studies aiming genes responsible for Mendelian and complex traits. The Quebec population, which is of mostly French descent, is considered an excellent model for such genetic epidemiological endeavours because it is a young founder population. Yet, the assessment of the founder effect has relied mostly on the observed distribution of monogenic diseases and on the analysis of the underlying mutations with investigations focusing on the Saguenay region. To eliminate this clinical bias and to obtain a more complete image of the genetic diversity, different regional populations of Quebec were investigated by analysing neutral markers that represent maternal, paternal and X chromosome lineages. Results indicate that Quebec does not appear more homogeneous nor significantly different from European populations. However, a series of regional founder effects, particularly visible at the level of rare variants, are observed. These effects can be explained by the successive migrations of descendants of the first immigrants from the initial sites of settlement towards the outer regions. Depending on the number of founders and their diversity, as well as on the degree of isolation and the magnitude of the interbreeding with the neighbouring or local populations, such as Amerindians or later migrants, the consequences of these regional founder effects are more or less detectable in the contemporary population.