GENLIB: an R package for the analysis of genealogical data

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.

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.

Polygenic risk of Alzheimer's disease in the Faroe Islands

INTRODUCTION

The Faroe Islands are a geographically isolated population in the North Atlantic with a similar prevalence Alzheimer's disease (AD) and all cause dementia as other European nations. However, the genetic risk underlying Alzheimer's disease and other dementia susceptibility has yet to be elucidated.

METHODS

Forty-nine single nucleotide polymorphisms (SNPs) were genotyped in 174 patients with AD and other dementias and 159 healthy controls. Single variant and polygenic risk score (PRS) associations, with/without APOE variability, were assessed by logistic regression. Performance was examined using receiver operating characteristics 'area under the curve' analysis (ROC AUC).

RESULTS

APOE rs429358 was associated with AD in the Faroese cohort after correction for multiple testing (OR=6.32, CI[3.98-10.05], p=6.31e^-15 ), with suggestive evidence for three other variants: NECTIN2 rs41289512 (OR 2.05, CI[1.20-3.51], p=0.01), HLA-DRB1 rs6931277 (OR 0.67, CI[0.48-0.94], p=0.02), and APOE rs7412 [ε2] (OR 0.28, CI[0.11-0.73], p=0.01). PRS were associated with AD with or without the inclusion of APOE (PRS^+APOE OR=4.5. CI[2.90-5.85, p=4.56e-¹⁵ and PRS^-APOE OR=1.53, CI[1.21-1.98], p=6.82e^-4 ). AD ROC AUC analyses demonstrated a PRS^+APOE AUC=80.3% and PRS^-APOE AUC=63.4%. However, PRS^+APOE was also significantly associated with all cause dementia (OR=3.39, CI[2.51-4.71], p= 2.50e^-14 ) with an AUC=76.9%, i.e. all cause dementia did show similar results albeit less significant.

DISCUSSION

In the Faroe Islands, SNP analyses highlighted APOE and immunogenomic variability in AD and dementia risk. PRS^+APOE , based on 25 SNPs/loci, had excellent sensitivity and specificity for Alzheimer's disease with AUC of 80.3%. High PRS were also associated with an earlier onset of late-onset AD.

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.

Amyotrophic lateral sclerosis in the Faroe Islands - a genealogical study

In the Faroe Islands, a clustering of amyotrophic lateral sclerosis (ALS) was observed on the geographically isolated island, Suðuroy. This study aims to estimate the frequency of familial ALS (fALS) in the Faroes including 43 patients diagnosed with ALS. Patients with fALS were identified through medical records and the Faroese Multi Generation Register. Firstly, fALS was recognized when occurring between first- or second-degree relatives. Secondly, families and individuals with fALS were recognized through pedigrees (≥3 cases within 3 generations). The prevalence of ALS was 3 times higher in Suðuroy compared to the nationwide prevalence. The frequency of fALS was at least 14% (n = 6) and mean survival time was 1.7 years shorter for fALS compared to sporatic ALS (p = 0.01. SD = 0.5, range 1.0-2.2). This study is suggestive of familial clustering in excess of expected for ALS and supports a genetic contribution to ALS in the Faroe Islands albeit environmental exposure within families cannot be excluded.

Genetic parameters of weeping teats in Italian Saanen and Alpine dairy goats and their relationship with milk production and somatic cell score

This paper reports a quantitative genetics analysis of weeping teats (WT), an abnormality of the mammary gland in goats. Weeping teats are characterized by milk oozing out of the teat or by the presence of multiple cysts near its base. This abnormality has been routinely recorded in Italian Alpine and Saanen goats since 2000 using a score of 0 or 1 (0 = defect not present, 1 = defect present). No information is available on the genetic background of WT or its relationship with production or other udder-related traits. Data obtained by the Italian National Sheep and Goat Breeders Association (Rome, Italy) were used to estimate both heritability of WT and its genetic correlation with milk yield, somatic cell score, and udder traits. The final data set used in the analysis included 2,178 Saanen and 2,309 Alpine primiparous goats kidding from 2009 to 2014. The pedigree data included 7,333 Saanen and 7,421 Alpines, respectively. A threshold-linear multivariate animal model was used to estimate variance and covariance components. A genealogical data analysis was also implemented, including genealogical data completeness, inbreeding, and identification of possible most recent common ancestors. On average, around 4 and 13% of primiparous Saanen and Alpine females kidding from 2009 to 2014 showed mammary gland abnormality, respectively. Weeping teats heritability was 0.27 and 0.26 for Saanen and Alpine, respectively. Genetic correlations between milk production or somatic cell score ranged from -0.16 in Saanen to 0.43 in Alpine, but the standard error of the estimates was very large. Positive genetic correlations were observed among WT and teat characteristics in both Saanen and Alpine. The average inbreeding of abnormality carriers was 2.4 and 5.1 for Saanen and Alpine, respectively. The genealogical data analysis identified 4 common ancestors of affected does in Saanen and 2 in Alpine. These results indicate that WT have a possible genetic background. A genome-wide association study might help in understanding the polygenic or monogenic determination of this abnormality.

Analysis of 100 high-coverage genomes from a pedigreed captive baboon colony

Baboons (genus Papio) are broadly studied in the wild and in captivity. They are widely used as a nonhuman primate model for biomedical studies, and the Southwest National Primate Research Center (SNPRC) at Texas Biomedical Research Institute has maintained a large captive baboon colony for more than 50 yr. Unlike other model organisms, however, the genomic resources for baboons are severely lacking. This has hindered the progress of studies using baboons as a model for basic biology or human disease. Here, we describe a data set of 100 high-coverage whole-genome sequences obtained from the mixed colony of olive (P. anubis) and yellow (P. cynocephalus) baboons housed at the SNPRC. These data provide a comprehensive catalog of common genetic variation in baboons, as well as a fine-scale genetic map. We show how the data can be used to learn about ancestry and admixture and to correct errors in the colony records. Finally, we investigated the consequences of inbreeding within the SNPRC colony and found clear evidence for increased rates of infant mortality and increased homozygosity of putatively deleterious alleles in inbred individuals.

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.

Deep genealogical analysis of a large cohort of participants in the CARTaGENE project (Quebec, Canada)

BACKGROUND

Genealogical analysis helps to better understand the genetic structure of populations. The population of Quebec (Canada) often serves as a model for this type of analysis, having one of the world's most complete genealogical databases.

AIM

The main objective of this study was to reconstruct, analyse and compare the ascending genealogies of participants to CARTaGENE, a project that aims at building a database on various aspects of public health.

SUBJECTS AND METHODS

In total, 5110 genealogies from four Quebec regions were reconstructed. Distribution of ancestors, completeness and depth of the genealogies, characteristics of immigrant ancestors and kinship and inbreeding coefficients were analysed.

RESULTS

Most genealogies go back to the 17th century, with a mean genealogical depth of 10 generations. Origins of immigrant ancestors are more diverse in the Montreal region, resulting in lower inbreeding and kinship among the participants from this region. Inbreeding and kinship values are mainly explained by remote genealogical links (from 6 to 11 generations).

CONCLUSION

Deep genealogies allowed for a precise measurement of the geographic origins of the participants' immigrant ancestors, as well as inbreeding and kinship ties in the population, which may be crucial for studies aiming to identify genetic variations associated with Mendelian or complex diseases.