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.

On the genes, genealogies, and geographies of Quebec.

Population genetic models only provide coarse representations of real-world ancestry. We used a pedigree compiled from 4 million parish records and genotype data from 2276 French and 20,451 French Canadian individuals to finely model and trace French Canadian ancestry through space and time. The loss of ancestral French population structure and the appearance of spatial and regional structure highlights a wide range of population expansion models. Geographic features shaped migrations, and we find enrichments for migration, genetic, and genealogical relatedness patterns within river networks across regions of Quebec. Finally, we provide a freely accessible simulated whole-genome sequence dataset with spatiotemporal metadata for 1,426,749 individuals reflecting intricate French Canadian population structure. Such realistic population-scale simulations provide opportunities to investigate population genetics at an unprecedented resolution.

GENLIB: new function to simulate haplotype transmission in large complex genealogies.

SUMMARY: Founder populations with deep genealogical data are well suited for investigating genetic variants contributing to diseases. Here, we present a major update of the genealogical analysis R package GENLIB, centered around a new function which can simulate the transmission of haplotypes from founders to probands along very large and complex user-specified genealogies. AVAILABILITY AND IMPLEMENTATION: The latest update of the GENLIB package (v1.1.9) contains the new gen.simuHaplo() function and is available on the CRAN repository and from https://github.com/R-GENLIB/GENLIB. Examples can be accessed at https://github.com/R-GENLIB/simuhaplo_functions.

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.

Genetic burden linked to founder effects in Saguenay-Lac-Saint-Jean illustrates the importance of genetic screening test availability.

The Saguenay-Lac-Saint-Jean (SLSJ) region located in the province of Quebec was settled in the 19th century by pioneers issued from successive migration waves starting in France in the 17th century and continuing within Quebec until the beginning of the 20th century. The genetic structure of the SLSJ population is considered to be the product of a triple founder effect and is characterised by a higher prevalence of some rare genetic diseases. Several studies were performed to elucidate the historical, demographic and genetic background of current SLSJ inhabitants to assess the origins of these rare disorders and their distribution in the population. Thanks to the development of new sequencing technologies, the genes and the variants responsible for the most prevalent conditions were identified. Combined with other resources such as the BALSAC population database, identifying the causal genes and the pathogenic variants allowed to assess the impacts of some of these founder mutations on the population health and to design precision medicine public health strategies based on carrier testing. Furthermore, it stimulated the establishment of many public programmes.We report here a review and an update of a subset of inherited disorders and founder mutations in the SLSJ region. Data were collected from published scientific sources. This work expands the knowledge about the current frequencies of these rare disorders, the frequencies of other rare genetic diseases in this population, the relevance of the carrier tests offered to the population, as well as the current available treatments and research about future therapeutic avenues for these inherited disorders.

Correspondence Between Genomic- and Genealogical/Coalescent-Based Inference of Homozygosity by Descent in Large French-Canadian Genealogies.

Research on the genetics of complex traits overwhelmingly focuses on the additive effects of genes. Yet, animal studies have shown that non-additive effects, in particular homozygosity effects, can shape complex traits. Recent investigations in human studies found some significant homozygosity effects. However, most human populations display restricted ranges of homozygosity by descent (HBD), making the identification of homozygosity effects challenging. Founder populations give rise to higher HBD levels. When deep genealogical data are available in a founder population, it is possible to gain information on the time to the most recent common ancestor (MRCA) from whom a chromosomal segment has been transmitted to both parents of an individual and in turn to that individual. This information on the time to MRCA can be combined with the time to MRCA inferred from coalescent models of gene genealogies. HBD can also be estimated from genomic data. The extent to which the genomic HBD measures correspond to the genealogical/coalescent measures has not been documented in founder populations with extensive genealogical data. In this study, we used simulations to relate genomic and genealogical/coalescent HBD measures. We based our simulations on genealogical data from two ongoing studies from the French-Canadian founder population displaying different levels of inbreeding. We simulated single-nucleotide polymorphisms (SNPs) in a 1-Mb genomic segment from a coalescent model in conjunction with the observed genealogical data. We compared genealogical/coalescent HBD to two genomic methods of HBD estimation based on hidden Markov models (HMMs). We found that genomic estimates of HBD correlated well with genealogical/coalescent HBD measures in both study genealogies. We described generation time to coalescence in terms of genomic HBD estimates and found a large variability in generation time captured by genomic HBD when considering each SNP. However, SNPs in longer segments were more likely to capture recent time to coalescence, as expected. Our study suggests that estimating the coalescent gene genealogy from the genomic data to use in conjunction with observed genealogical data could provide valuable information on HBD.

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.

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.

Relaxed Selection During a Recent Human Expansion.

Humans have colonized the planet through a series of range expansions, which deeply impacted genetic diversity in newly settled areas and potentially increased the frequency of deleterious mutations on expanding wave fronts. To test this prediction, we studied the genomic diversity of French Canadians who colonized Quebec in the 17th century. We used historical information and records from ∼4000 ascending genealogies to select individuals whose ancestors lived mostly on the colonizing wave front and individuals whose ancestors remained in the core of the settlement. Comparison of exomic diversity reveals that: (i) both new and low-frequency variants are significantly more deleterious in front than in core individuals, (ii) equally deleterious mutations are at higher frequencies in front individuals, and (iii) front individuals are two times more likely to be homozygous for rare very deleterious mutations present in Europeans. These differences have emerged in the past six to nine generations and cannot be explained by differential inbreeding, but are consistent with relaxed selection mainly due to higher rates of genetic drift on the wave front. Demographic inference and modeling of the evolution of rare variants suggest lower effective size on the front, and lead to an estimation of selection coefficients that increase with conservation scores. Even though range expansions have had a relatively limited impact on the overall fitness of French Canadians, they could explain the higher prevalence of recessive genetic diseases in recently settled regions of Quebec.

Socio-economic status and fertility decline: Insights from historical transitions in Europe and North America.

The timings of historical fertility transitions in different regions are well understood by demographers, but much less is known regarding their specific features and causes. In the study reported in this paper, we used longitudinal micro-level data for five local populations in Europe and North America to analyse the relationship between socio-economic status and fertility during the fertility transition. Using comparable analytical models and class schemes for each population, we examined the changing socio-economic differences in marital fertility and related these to common theories on fertility behaviour. Our results do not provide support for the hypothesis of universally high fertility among the upper classes in pre-transitional society, but do support the idea that the upper classes acted as forerunners by reducing their fertility before other groups. Farmers and unskilled workers were the latest to start limiting their fertility. Apart from these similarities, patterns of class differences in fertility varied significantly between populations.

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.

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.

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.

Deep human genealogies reveal a selective advantage to be on an expanding wave front.

Since their origin, human populations have colonized the whole planet, but the demographic processes governing range expansions are mostly unknown. We analyzed the genealogy of more than one million individuals resulting from a range expansion in Quebec between 1686 and 1960 and reconstructed the spatial dynamics of the expansion. We find that a majority of the present Saguenay Lac-Saint-Jean population can be traced back to ancestors having lived directly on or close to the wave front. Ancestors located on the front contributed significantly more to the current gene pool than those from the range core, likely due to a 20% larger effective fertility of women on the wave front. This fitness component is heritable on the wave front and not in the core, implying that this life-history trait evolves during range expansions.

Admixed ancestry and stratification of Quebec regional populations.

Population stratification results from unequal, nonrandom genetic contribution of ancestors and should be reflected in the underlying genealogies. In Quebec, the distribution of Mendelian diseases points to local founder effects suggesting stratification of the contemporary French Canadian gene pool. Here we characterize the population structure through the analysis of the genetic contribution of 7,798 immigrant founders identified in the genealogies of 2,221 subjects partitioned in eight regions. In all but one region, about 90% of gene pools were contributed by early French founders. In the eastern region where this contribution was 76%, we observed higher contributions of Acadians, British and American Loyalists. To detect population stratification from genealogical data, we propose an approach based on principal component analysis (PCA) of immigrant founders' genetic contributions. This analysis was compared with a multidimensional scaling of pairwise kinship coefficients. Both methods showed evidence of a distinct identity of the northeastern and eastern regions and stratification of the regional populations correlated with geographical location along the St-Lawrence River. In addition, we observed a West-East decreasing gradient of diversity. Analysis of PC-correlated founders illustrates the differential impact of early versus latter founders consistent with specific regional genetic patterns. These results highlight the importance of considering the geographic origin of samples in the design of genetic epidemiology studies conducted in Quebec. Moreover, our results demonstrate that the study of deep ascending genealogies can accurately reveal population structure.

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.

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.

Mutability of Y-chromosomal microsatellites: rates, characteristics, molecular bases, and forensic implications.

Nonrecombining Y-chromosomal microsatellites (Y-STRs) are widely used to infer population histories, discover genealogical relationships, and identify males for criminal justice purposes. Although a key requirement for their application is reliable mutability knowledge, empirical data are only available for a small number of Y-STRs thus far. To rectify this, we analyzed a large number of 186 Y-STR markers in nearly 2000 DNA-confirmed father-son pairs, covering an overall number of 352,999 meiotic transfers. Following confirmation by DNA sequence analysis, the retrieved mutation data were modeled via a Bayesian approach, resulting in mutation rates from 3.78 × 10(-4) (95% credible interval [CI], 1.38 × 10(-5) - 2.02 × 10(-3)) to 7.44 × 10(-2) (95% CI, 6.51 × 10(-2) - 9.09 × 10(-2)) per marker per generation. With the 924 mutations at 120 Y-STR markers, a nonsignificant excess of repeat losses versus gains (1.16:1), as well as a strong and significant excess of single-repeat versus multirepeat changes (25.23:1), was observed. Although the total repeat number influenced Y-STR locus mutability most strongly, repeat complexity, the length in base pairs of the repeated motif, and the father's age also contributed to Y-STR mutability. To exemplify how to practically utilize this knowledge, we analyzed the 13 most mutable Y-STRs in an independent sample set and empirically proved their suitability for distinguishing close and distantly related males. This finding is expected to revolutionize Y-chromosomal applications in forensic biology, from previous male lineage differentiation toward future male individual identification.

Genealogical analysis of maternal and paternal lineages in the Quebec population.

The Quebec population is one of the rare populations of its size for which genealogical information is available for an uninterrupted period of almost four centuries. This allows for in-depth studies on the formation and evolution of a young founder population. Using data from two major population registers, in this study we focus on the maternal and paternal lineages (i.e., strictly female or male genealogical lines) that can be traced back within the Quebec genealogies. Through the analysis of these lineages it is possible to characterize the founders who transmitted to the contemporary population their mitochondrial (for females) and Y-chromosome (for males) DNA. The basic material consists of 2,221 ascending genealogies of subjects who married in the Quebec population between 1945 and 1965. On average, more than nine generations of ancestors were identified among the lineages. Analyses of maternal and paternal lineages show that the number of paternal founders is higher and their origins and genetic contributions are more variable than that of maternal founders, leading to a larger effective population size and greater diversity of Y chromosomes than of mtDNA. This is explained for the most part by differential migratory patterns among male and female founders of the Quebec population. Comparisons of sex-specific genetic contributions with total genetic contribution showed a strong correlation between the two values, with some discrepancies related to sex ratio differences among the founders' first descendants.

Is there a trade-off between fertility and longevity? A comparative study of women from three large historical databases accounting for mortality selection.

Frontier populations provide exceptional opportunities to test the hypothesis of a trade-off between fertility and longevity. In such populations, mechanisms favoring reproduction usually find fertile ground, and if these mechanisms reduce longevity, demographers should observe higher postreproductive mortality among highly fertile women. We test this hypothesis using complete female reproductive histories from three large demographic databases: the Registre de la population du Québec ancien (Université de Montréal), which covers the first centuries of settlement in Quebec; the BALSAC database (Université du Québec à Chicoutimi), including comprehensive records for the Saguenay-Lac-St-Jean (SLSJ) in Quebec in the nineteenth and twentieth centuries; and the Utah Population Database (University of Utah), including all individuals who experienced a vital event on the Mormon Trail and their descendants. Together, the three samples allow for comparisons over time and space, and represent one of the largest set of natural fertility cohorts used to simultaneously assess reproduction and longevity. Using survival analyses, we found a negative influence of parity and a positive influence of age at last child on postreproductive survival in the three populations, as well as a significant interaction between these two variables. The effect sizes of all these parameters were remarkably similar in the three samples. However, we found little evidence that early fertility affects postreproductive survival. The use of Heckman's procedure assessing the impact of mortality selection during reproductive ages did not appreciably alter these results. We conclude our empirical investigation by discussing the advantages of comparative approaches.