Search

Human evolution involves population splits, size fluctuations, founder effects, and admixture. Population history reconstruction based on genetic diversity data routinely relies on simple demographic models while projecting the past. No specific demographic assumptions are needed to understand the genetic structure of the founder population of Québec. Because genealogy and genetics are intimately related, we used descending genealogies of this population to pursue the fate of its founder lineages. Maternal and paternal lines reflect the transmission of mtDNA and the Y-chromosome, respectively. We followed their transmission in real-time, from the 17th century down to its 20th-century population. We counted the number of married children of immigrants (i.e., their effective family size, EFS), estimated the proportion of successful immigrants in terms of their survival ratio, and assessed net growth rates and extinction. Likewise, we evaluated the same parameters for their Québec-born descendants. The survival ratio of the first immigrants was the highest and declined over time in association with the decreasing immigrants’ EFS. Parents with high EFS left plentiful married progeny, putting EFS as the most important variable determining the parental demographic success throughout time for generations ahead. The 17th and 18th-century immigrants bear the most remarkable demographic and genetic impact on the 20th-century population of Québec. Lessons learned from Québec genealogies can teach us about the consequences of founder effects and migrations through real people’s history. The effective family size of immigrant founders predicts their long-term demographic outcome.

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.

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.

[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., European Regional Development Fund, Grant/Award Number: PGC2018‐095931‐B‐100; Fonds de Recherche du Québec ‐ Santé; Fonds de Recherche du Québec‐Société et Culture; Natural Sciences and Engineering Research Council of Canada.

Farming was first introduced to southeastern Europe in the mid-7th millennium BCE - brought by migrants from Anatolia who settled in the region before spreading throughout Europe. However, the dynamics of the interaction between the first farmers and the indigenous hunter-gatherers remain poorly understood because of the near absence of ancient DNA from the region. We report new genome-wide ancient DNA data from 204 individuals-65 Paleolithic and Mesolithic, 93 Neolithic, and 46 Copper, Bronze and Iron Age-who lived in southeastern Europe and surrounding regions between about 12,000 and 500 BCE. We document that the hunter-gatherer populations of southeastern Europe, the Baltic, and the North Pontic Steppe were distinctive from those of western Europe, with a West-East cline of ancestry. We show that the people who brought farming to Europe were not part of a single population, as early farmers from southern Greece are not descended from the Neolithic population of northwestern Anatolia that was ancestral to all other European farmers. The ancestors of the first farmers of northern and western Europe passed through southeastern Europe with limited admixture with local hunter-gatherers, but we show that some groups that remained in the region mixed extensively with local hunter-gatherers, with relatively sex-balanced admixture compared to the male-biased hunter-gatherer admixture that we show prevailed later in the North and West. After the spread of farming, southeastern Europe continued to be a nexus between East and West, with intermittent steppe ancestry, including in individuals from the Varna I cemetery and associated with the Cucuteni-Trypillian archaeological complex, up to 2,000 years before the Steppe migration that replaced much of northern Europe9s population.

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. In evolutionary biology, ‘mother’s curse’ refers to the possibility of passing on harmful mutations through mitochondria. Here evidence is presented for the mother’s curse in action over 290 years in a human population.

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.

We performed exome analysis in two affected siblings with severe intellectual disability (ID), microcephaly and spasticity from an Ashkenazi Jewish consanguineous family. We identified only one rare variant, a missense in SLC1A4 (c. 766G>A [p. E256K]), that is homozygous in both siblings but not in any of their 11 unaffected siblings or their parents (Logarithm of odds, LOD score: 2.6). This variant is predicted damaging. We genotyped 450 controls of Ashkenazi Jewish ancestry and identified only 5 individuals who are heterozygous for this variant (minor allele frequency: 0.0056). SLC1A4 (ASCT1) encodes a transporter for neutral aminoacids such as alanine, serine, cysteine and threonine. l-Serine is essential for neuronal survival and differentiation. Indeed, l-serine biosynthesis disorders affect brain development and cause severe ID. In the brain, l-serine is synthesized in astrocytes but not in neurons. It has been proposed that ASCT1 mediates the uptake of l-serine into neurons and the release of glia-borne l-serine to neighboring cells. SLC1A4 disruption may thus impair brain development and function by decreasing the levels of l-serine in neurons. The identification of additional families with mutations in SLC1A4 would be necessary to confirm its involvement in ID.

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.

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.

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.

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.

L’identite ethnique d’un individu releve d’un sentiment d’appartenance plus ou moins prononce a un groupe possedant certaines caracteristiques propres. Dans les recensements et de nombreuses enquetes quantitatives, une question sur l’origine ethnique est posee aux participants et on peut se demander dans quelle mesure les reponses fournies correspondent aux origines ancestrales des repondants. A partir d’informations genealogiques portant sur pres de 400 individus residant dans la region de la Gaspesie (Est du Quebec), cette etude compare l’origine declaree par les participants a celles des premiers immigrants de leur ascendance etablis sur le sol quebecois. Les genealogies des participants ont ete reconstruites, les ancetres immigrants identifies et leurs origines geographiques documentees. Les resultats montrent que l’origine declaree par les participants gaspesiens trouve une correspondance sur le plan des origines ancestrales. Bien que presque tous les participants aient des ancetres immigrants d’origines diverses, le groupe d’appartenance est, dans la majorite des cas, le plus represente parmi ces ancetres. En particulier, les repondants ont davantage tendance a s’identifier a leurs origines paternelles qu’a leurs origines maternelles lorsque celles-ci different. Malgre une histoire de peuplement diversifiee et un metissage important, les Gaspesiens ont donc conserve un sentiment d’appartenance ethnique qui reflete aujourd’hui encore une certaine realite demographique et genetique.

Dans le cadre d’un programme de recherche sur la diversité génétique des populations régionales du Québec, nous avons mesuré la contribution des ancêtres amérindiens au pool génique de 794 participants résidant dans les régions de Montréal, du Saguenay—Lac-Saint-Jean (SLSJ), de la Gaspésie et de la Côte-Nord. Les ancêtres amérindiens ont été identifiés à partir de reconstructions généalogiques (fichier BALSAC) et de données génétiques (ADN mitochondrial) à l’aide d’une approche qui repose sur l’exploitation conjointe de ces deux types de données. Les résultats indiquent que plus de la moitié des participants ont au moins un ancêtre amérindien dans leur généalogie. Ceci veut dire que la majorité des participants sont porteurs de gènes reçus de fondateurs amérindiens. Cependant la contribution génétique totale de ces ancêtres aux quatre groupes régionaux demeure faible. En effet, elle est de moins de 1 % au SLSJ et à Montréal alors qu’elle dépasse à peine 1 % sur la Côte-Nord et en Gaspésie., In the context of a research programme on the genetic diversity of regional populations in Quebec, we measured the contribution of Amerindian ancestors to the gene pool of 794 participants living in the Montreal, Saguenay-Lac-St-Jean (SLSJ), Gaspe Peninsula and North Shore regions. Amerindian ancestors were identified through genealogical reconstructions (BALSAC database) and genetic data (mitochondrial DNA) using an approach relying on the joint exploitation of the two types of data. Results show that more than half of the participants have at least one Amerindian ancestor in their genealogy. This means that a majority of participants carry genes received from Amerindian founders. However the total genetic contribution of these ancestors to the four regional groups is low, reaching less than 1 % in SLSJ and Montreal and slightly more than 1 % on the North Shore and in the Gaspe Peninsula.

Here we report the Simons Genome Diversity Project data set: high quality genomes from 300 individuals from 142 diverse populations. These genomes include at least 5.8 million base pairs that are not present in the human reference genome. Our analysis reveals key features of the landscape of human genome variation, including that the rate of accumulation of mutations has accelerated by about 5% in non-Africans compared to Africans since divergence. We show that the ancestors of some pairs of present-day human populations were substantially separated by 100,000 years ago, well before the archaeologically attested onset of behavioural modernity. We also demonstrate that indigenous Australians, New Guineans and Andamanese do not derive substantial ancestry from an early dispersal of modern humans; instead, their modern human ancestry is consistent with coming from the same source as that of other non-Africans., Facultad de Ciencias Naturales y Museo

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.

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.

We report an integrated analysis of nuclear (autosomal, X- and Y-chromosome) short tandem repeat (STR) data and mtDNA D-loop sequences obtained in the same set of 22 Native populations from across the Americas. A north to south gradient of decreasing population diversity was observed, in agreement with a settlement of the Americas from the extreme northwest of the continent. This correlation is stronger with "least cost distances," which consider the coasts as facilitators of migration. Continent-wide estimates of population structure are highest for the Y-chromosome and lowest for the autosomes, consistent with the effective size of the different marker systems examined. Population differentiation is highest in East South America and lowest in Meso America and the Andean region. Regional analyses suggest a deviation from mutation-drift equilibrium consistent with population expansion in Meso America and the Andes and population contraction in Northwest and East South America. These data hint at an early divergence of Andean and non-Andean South Americans and at a contrasting demographic history for populations from these regions.

Was the past genetic contribution of women and men to the current human population equal? Was polygyny (excess of breeding women) present among hominid lineages? We addressed these questions by measuring the ratio of population recombination rates between the X chromosome and the autosomes, rho(X)/rho(A). The X chromosome recombines only in female meiosis, whereas autosomes undergo crossovers in both sexes; thus, rho(X)/rho(A) reflects the female-to-male breeding ratio, beta. We estimated beta from rho(X)/rho(A) inferred from genomic diversity data and calibrated with recombination rates derived from pedigree data. For the HapMap populations, we obtained beta of 1.4 in the Yoruba from West Africa, 1.3 in Europeans, and 1.1 in East Asian samples. These values are consistent with a high prevalence of monogamy and limited polygyny in human populations. More mutations occur during male meiosis as compared to female meiosis at the rate ratio referred to as alpha. We show that at alpha not equal 1, the divergence rates and genetic diversities of the X chromosome relative to the autosomes are complex functions of both alpha and beta, making their independent estimation difficult. Because our estimator of beta does not require any knowledge of the mutation rates, our approach should allow us to dissociate the effects of alpha and beta on the genetic diversity and divergence rate ratios of the sex chromosomes to the autosomes.

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 x 10(-4) (95% credible interval [CI], 1.38 x 10(-5) - 2.02 x 10(-3)) to 7.44 x 10(-2) (95% Cl, 6.51 x 10(-2) - 9.09 x 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.

Comparative studies of linkage disequilibrium (LD) can provide insights into human demographic history. Here, we characterize LD in six Native American populations using seven microsatellite markers in Xq13, a region of the genome extensively studied in populations around the world. Native Americans show relatively low diversity and high LD, in agreement with recent genome-wide survey and a scenario of sequential founder effects accompanying human population dispersal around the globe. LD in Native Americans is similar to that observed in some recently described small population isolates and higher than in large European isolates (e.g., Finns), which have been extensively analyzed in medical genetics studies. Haplotype analyses are consistent with a colonization of the New World by a differentiated East Asian population, followed by extensive genetic drift in the Americas.

PCSK9 is a liver-secreted blood protein that promotes the degradation of low-density lipoprotein receptors, leading to reduced hepatic uptake of plasma cholesterol. Nucleotide variations in its gene have been linked to hypo- and hyper-cholesterolemia. Two nonsense mutations, Y142X and C679X, are associated to lifelong hypocholesterolemia and a remarkable protection against coronary heart disease (CHD) in African Americans. The aim of this study was to determine the frequency of these cardioprotective mutations in West Africans. Subjects (n = 520) from different ethnic groups were recruited in Burkina-Faso, Benin, and Togo. Only the C679X mutation was detected. All carriers were heterozygous. The overall heterozygosity frequency was 3.3%. It varied significantly among ethnic groups, ranging from 0% to 6.9%. The overall high frequency of the cardioprotective C679X mutation in Africa may contribute to the lower incidence of CHD on this continent. The interethnic frequency differences may reflect historical settlement and migration patterns in the region, possibly combined with positive selection for the mutation driven by yet-unknown environmental factors.

The Altaian Kazakhs, a Turkic speaking group, now reside in the southern part of the Altai Republic in south-central Russia. According to historical accounts, they are one of several ethnic and geographical subdivisions of the Kazakh nomadic group that migrated from China and Western Mongolia into the Altai region during the 19th Century. However, their population history of the Altaian Kazakhs and the genetic relationships with other Kazakh groups and neighboring Turkic-speaking populations is not well understood. To begin elucidating their genetic history, we analyzed the mtDNAs from 237 Altaian Kazakhs through a combination of SNP analysis and HVS1 sequencing. This analysis revealed that their mtDNA gene pool was comprised of roughly equal proportions of East (A-G, M7, M13, Y and Z) and West (H, HV, pre-HV, R, IK, JT, X, U) Eurasian haplogroups, with the haplotypic diversity within haplogroups C, D, H, and U being particularly high. This pattern of diversity likely reflects the complex interactions of the Kazakhs with other Turkic groups, Mongolians, and indigenous Altaians. Overall, these data have important implications for Kazakh population history, the genetic prehistory of the Altai-Sayan region, and the phylogeography of major mitochondrial lineages in Eurasia.

In this article we present a new heuristic approach (informative recombinations, InfRec) to analyze recombination density at the sequence level. InfRec is intuitive and easy and combines previously developed methods that (i) resolve genotypes into haplotypes, (ii) estimate the minimum number of recombinations, and (iii) evaluate the fraction of informative recombinations. We tested this approach in its sliding-window version on 117 genes from the SeattleSNPs program, resequenced in 24 African-Americans (AAs) and 23 European-Americans (EAs). We obtained population recombination rate estimates (ρobs) of 0.85 and 0.37 kb−1 in AAs and EAs, respectively. Coalescence simulations indicated that these values account for both the recombinations and the gene conversions in the history of the sample. The intensity of ρobs varied considerably along the sequence, revealing the presence of recombination hotspots. Overall, we observed ∼80% of recombinations in one-third and ∼50% in only 10% of the sequence. InfRec performance, tested on published simulated and additional experimental data sets, was similar to that of other hotspot detection methods. Fast, intuitive, and visual, InfRec is not constrained by sample size limitations. It facilitates understanding data and provides a simple and flexible tool to analyze recombination intensity along the sequence.

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.

Genetic variability of the compound interrupted microsatellite DXS1238, in intron 44 of the dystrophin gene, provides evidence for a complex structure of the ancestral population that led to the emergence of modern humans. We sequenced DXS1238 in 600 X-chromosomes from all over the world. Forty four percent of African-specific chromosomes belong to the ancestral lineage that did not participate in the out-of-Africa expansion and subsequent colonization of other continents. Based on the coalescence analysis these lineages separated from those that contributed to the out-of-Africa expansion 366 +/- 136 thousands years ago (Kya). Independently, the analysis of the variance in the repeat length and of the decay of the ancestral alleles of the two DXS1238 repeats, GT and GA, dates this separation at more than 200 Kya. This suggests a complex demographic history and genetic structure of the African melting pot that led to the emergence of modern humans and their out-of-Africa migration. The subsequent subdivisions of human populations among different continents appear to be preceded by even more structured population history within Africa itself, which resulted from a restricted gene flow between lineages allowing for genetic differences to accumulate. If the transition to modern humans occurred during that time, it necessarily follows that genes associated with this transformation spread between subpopulations via gene flow. Otherwise, in spite of subsequent anatomical variation, Homo sapiens as a species could have emerged in Africa already between 300 and 200 Kya, i.e. before the mitochondrial DNA and well before the Y-chromosome most recent common ancestors.

Human colorectal carcinoma (Caco-2) cells undergo in culture spontaneous enterocytic differentiation, characterized by polarization and appearance of the functional apical brush border membrane. To provide insights into the biology of differentiation, we have performed a comparative proteomic analysis of the plasma membranes from proliferating cells (PCs) and the apical membranes from differentiated cells (DCs). Proteins were resolved by SDS-PAGE, in-gel digested and analyzed by RP-LC and MS/MS. Alternatively, proteins were digested in solution, and tryptic peptides were labeled with isotopic tags and analyzed by 2-D LC followed by MS/MS. Among the 1125 proteins identified in both proteomes, 76 were found to be significantly increased in the membranes of DCs and 61 were increased in PCs. Majority of the proteins increased in the apical membranes were metabolic enzymes, proteins involved in the maintenance of cellular structure, transmembrane transporters, and proteins regulating vesicular transport. In contrast, majority of the proteins increased in the membranes of PCs were involved in gene expression, protein synthesis, and folding. Both groups contained many novel proteins with yet to be identified functions, which could provide potential new markers of the intestinal cells or of colorectal cancer.

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. Hum Mutat 28(5), 441–450, 2007. © 2007 Wiley-Liss, Inc.

Sudmant, Peter H. et al., 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., This project has been funded in part with federal funds from the National Cancer Institute, NIH, under contract HHSN26120080001E. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does the mention of trade names, commercial products, or organizations imply endorsement by the U.S. government. This research was supported in part by the Intramural Research Program of the NIH, National Cancer Institute, Center for Cancer Research. This work was also partly supported by NIH grant 2R01HG002385 and a grant (11631) from the Paul G. Allen Family Foundation to E.E.E. The sequencing for this study was supported by a grant from the Simons Foundation to D.R. (SFARI 280376) and by a HOMINID grant from the NSF to D.R. (BCS-1032255). T.K. is supported by a European Research Council Starting Investigator grant (FP7 - 26213). R.S. and S.D. received support from the Ministry of Education and Science, Russian Federation (14.Z50.31.0010). H.S., E.M., R.V., and M.M. are supported by Institutional Research Funding from the Estonian Research Council IUT24-1 and by the European Regional Development Fund (European Union) through the Centre of Excellence in Genomics to Estonian Biocentre and University of Tartu. S.A.T. is supported by NIH grants 5DP1ES022577 05, 1R01DK104339-01, and 1R01GM113657-01. C.T.-S. is supported by Wellcome Trust grant 098051. C.M.B. is supported by the NSF (award numbers 0924726 and 1153911). E.E.E. and D.R. are investigators of the Howard Hughes Medical Institute. Data are deposited into ENA (PRJEB9586 or ERP010710), and variant calls are deposited in dbVar (PRJNA285786). E.E.E. is on the scientific advisory board of DNAnexus, Incorporated, and is a consultant for Kunming University of Science and Technology (KUST) as part of the 1000 China Talent Program.

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. Electronic supplementary material The online version of this article (doi:10.1186/s12859-015-0581-5) contains supplementary material, which is available to authorized users.

The predominance of the T14484C mutation in French Canadians with Leber hereditary optic neuropathy is due to a founder effect. By use of genealogical reconstructions of maternal lineages, a woman married in Quebec City in 1669 is identified as the shared female ancestor for 11 of 13 affected individuals, who were previously not known to be related. These individuals carry identical mitochondrial haplogroups. The current geographic distribution of French Canadian cases overlaps with that of the founder's female descendants in 1800. This is the first example of genealogical reconstruction to identify the introduction of a mitochondrial mutation by a woman in a founder population.

Senataxin recently was identified as the mutated gene in ataxia-oculomotor apraxia 2, which is characterized by ataxia, oculomotor apraxia, and increased α-fetoprotein levels. In this study, we evaluated 24 ataxic patients from 10 French-Canadian families. All cases have a homogeneous phenotype consisting of a progressive ataxia appearing between 2 and 20 (mean age, 14.8) years of age with associated dysarthria, saccadic ocular pursuit, distal amyotrophy, sensory and motor neuropathy, and increased α-fetoprotein levels but absence of oculomotor apraxia. Linkage disequilibrium was observed with markers in the ataxia-oculomotor apraxia 2 locus on chromosome 9q34. We have identified four mutations in senataxin in the French-Canadian population including two novel missense mutations: the 5927TG mutation changes the leucine encoded by codon 1976 to an arginine in the helicase domain (L1976R), and the 193GA mutation changes a glutamic acid encoded by codon 65 into a lysine in the N-terminal domain of the protein (E65K). The common L1976R mutation is shared by 17 of 20 (85%) carrier chromosomes. The study of this large French-Canadian cohort better defines the phenotype of this ataxia and presents two novel mutations in senataxin including the more common founder mutation in the French-Canadian population. Ann Neurol 2005;57:408–414

The central role of methylenetetrahydrofolate reductase (MTHFR) in the folate metabolism renders MTHFR gene polymorphisms (C677T and A1298C) potential modulators of a variety of disorders whose development depends on folate/homocysteine imbalance. Here, we provide additional evidence on the protective role of these polymorphisms in acute lymphoblastic leukemia (ALL), the most common pediatric cancer. A case-control study was conducted in 270 ALL patients and 300 healthy controls of French-Canadian origin. The TT677/AA1298 and CC677/CC1298 individuals were associated with reduced risk of ALL (crude odds ratio [OR] = 0.4; 95% confidence interval [CI], 0.2-0.9; and OR = 0.3; 95% CI, 0.1-0.6; respectively). Further stratification in patients born before and after January 1996 (approximate time of Health Canada recommendation for folic acid supplement in pregnancy) revealed that the protective effect of MTHFR variants is accentuated and present only in children born before 1996. Similar results were obtained when a transmission disequilibrium test was performed on a subset of children (n = 95) in a family-based study. This finding suggests gene-environment interaction and its role in the susceptibility to childhood ALL, which is consistent with previous findings associating either folate deficiency or MTHFR polymorphisms with risk of leukemia.

The modification of cellular DNA by environmental substances is thought to be a crucial event in chemical induced carcinogenesis. Among the environmental carcinogens, aromatic amines are known for the fact that they can induce several types of cancers through the formation of so-called DNA adducts. We took advantage of the potential of the SELEX method to select for highly specific RNA ligands that recognize specific genotoxic aromatic amines. The aromatic amine 4,4'-methylenedianiline (MDA) was used as a target. Following in vitro selection, we obtained specific MDA-binding RNA molecules based on an affinity chromatography assay. These results open the possibility of using the SELEX technique to generate RNA molecules as diagnostic tools for the detection of DNA damaging compounds and ultimately DNA adducts.

Although Africa has played a central role in human evolutionary history, certain studies have suggested that not all contemporary human genetic diversity is of recent African origin. We investigated 35 simple polymorphic sites and one T(n) microsatellite in an 8-kb segment of the dystrophin gene. We found 86 haplotypes in 1,343 chromosomes from around the world. Although a classical out-of-Africa topology was observed in trees based on the variant frequencies, the tree of haplotype sequences reveals three lineages accounting for present-day diversity. The proportion of new recombinants and the diversity of the T(n) microsatellite were used to estimate the age of haplotype lineages and the time of colonization events. The lineage that underwent the great expansion originated in Africa prior to the Upper Paleolithic (27,000-56,000 years ago). A second group, of structurally distinct haplotypes that occupy a central position on the tree, has never left Africa. The third lineage is represented by the haplotype that lies closest to the root, is virtually absent in Africa, and appears older than the recent out-of-Africa expansion. We propose that this lineage could have left Africa before the expansion (as early as 160,000 years ago) and admixed, outside of Africa, with the expanding lineage. Contemporary human diversity, although dominated by the recently expanded African lineage, thus represents a mosaic of different contributions.

Within Europe, classical genetic markers, nuclear autosomal and Y-chromosome DNA polymorphisms display an east-west frequency gradient. This has been taken as evidence for the westward migration of Neolithic farmers from the Middle East. In contrast, most studies of mtDNA variation in Europe and the Middle East have not revealed clinal distributions. Here we report an analysis of dys44 haplotypes, consisting of 35 polymorphisms on an 8 kb segment of the dystrophin gene on Xp21, in a sample of 1203 Eurasian chromosomes. Our results do not show a significant genetic structure in Europe, though when Middle Eastern samples are included a very low but significant genetic structure, rooted in Middle Eastern heterogeneity, is observed. This structure was not correlated to either geography or language, indicating that neither of these factors are a barrier to gene flow within Europe and/or the Middle East. Spatial autocorrelation analysis did not show clinal variation from the Middle East to Europe, though an underlying and ancient east-west cline across the Eurasian continent was detected. Clines provide a strong signal of ancient major population migration(s), and we suggest that the observed cline likely resulted from an ancient, bifurcating migration out of Africa that influenced the colonizing of Europe, Asia and the Americas. Our study reveals that, in addition to settlements from the Near East, Europe has been influenced by other major population movements, such as expansion(s) from Asia, as well as by recent gene flow from within Europe and the Middle East.

Acute lymphoblastic leukaemia (ALL) is the most common childhood cancer. Genetic variants in the coding regions of the mismatch repair genes MLH1 (Ile-219Val) and MSH3 (Arg-940Glu and Thr-1036Ala) could contribute to an individual's susceptibility as modifiers in leukaemogenesis. To investigate this possibility, we conducted a case-control study on 287 children with ALL and 320 healthy controls both of French-Canadian origin. MLH1 and MSH3 variants, when taken independently, seemed to play little or no role in the aetiology of childhood ALL. However, when the MLH1 genotypes were combined with genotypes previously shown to influence ALL susceptibility, we found that the MLH1 variant Val-219 further increases the risk of GSTM1 null and CYP1A1*2A genotypes [combined odds ratio (OR) = 6.0, P = 0.002] as well as that of CYP2E1*5 (OR = 15.8, P = 0.001). No association was found with MSH3 variants. This study suggests an association of leukaemogenesis in children with both xenobiotic metabolism and DNA repair, and thus points to the effect of environmental exposure.

Several isoforms of estrogen receptor (ER) β (also known as NR3A2) have been reported, including variants with different N-terminal ends. In rodents, two in-frame initiation codons (ATGs) are used to produce proteins of 530 and 549 amino acids, respectively. In humans, the upstream ATG is out of frame in all clones reported, until recently, when human clones with an extra A-T base pair placing the upstream ATG in frame were reported. The authors suggested that this could represent a novel polymorphism in the ERβ gene. Because human ERβ548 (hERβ548) and hERβ530 display different functional characteristics in vitro, it is of interest to determine if this variant constitutes a polymorphism in human populations. We therefore determined the frequency of this novel isoform in several populations including African (n = 96), Caucasian (n = 100), and Asian (n = 128) subjects using denaturing HPLC. We did not detect any alleles that correspond to hERβ548 in these samples or in additional samples of heterogeneous origin. It is concluded that hERβ548 is not a common variant in Africans, Caucasians, or Asians.