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.     

Genetic population structure of the white sifaka (Propithecus verreauxi verreauxi) at Beza Mahafaly Special Reserve,southwest Madagascar (1992-2001)

Gene flow within and between social groups is contingent on behaviourally mediated patterns of mating and dispersal. To understand how these patterns affect the genetic structure of primate populations, long-term data are required. In this study, we analyse 10 years of demographic and genetic data from a wild lemur population (Propithecus verreauxi verreauxi) at Beza Mahafaly Special Reserve, southwest Madagascar. Our goal is to specify how patterns of mating and dispersal determine kinship and genetic diversity among animals in the population. Specifically, we use microsatellite, parentage, and census data to obtain estimates of genetic subdivision (FST), within group homozygosity (FIS), and relatedness (r) within and among social groups in the population. We analyse different classes of individuals (i.e. adults, offspring, males, females) separately in order to discern which classes most strongly influence aspects of population structure. Microsatellite data reveal that, across years, offspring are consistently more heterozygous than expected within social groups (FIS mean = -0.068) while adults show both positive and negative deviations from expected genotypic frequencies within groups (FIS mean = 0.003). Offspring cohorts are more genetically subdivided than adults (FST mean = 0.108 vs. 0.052) and adult females are more genetically subdivided than adult males (FST mean = 0.098 vs. 0.046). As the proportion of females in social groups increases, the proportion of offspring sired by resident males decreases. Offspring are characterized by a heterozygote excess as resident males (vs. nonresident males) sire the majority of offspring within groups. We link these genetic data to patterns of female philopatry, male dispersal, exogamy, and offspring sex-ratio. Overall, these data reveal how mating and dispersal tactics influence the genetic population structure in this species.     

Population structure and its effects on patterns of nucleotide polymorphism in teosinte (Zea mays ssp. parviglumis)

Surveys of nucleotide diversity in the wild ancestor of maize, Zea mays ssp. parviglumis, have revealed genomewide departures from the standard neutral equilibrium (NE) model. Here we investigate the degree to which population structure may account for the excess of rare polymorphisms frequently observed in species-wide samples. On the basis of sequence data from five nuclear and two chloroplast loci, we found significant population genetic structure among seven subpopulations from two geographic regions. Comparisons of estimates of population genetic parameters from species-wide samples and subpopulation-specific samples showed that population genetic subdivision influenced observed patterns of nucleotide polymorphism. In particular, Tajima's D was significantly higher (closer to zero) in subpopulation-specific samples relative to species-wide samples, and therefore more closely corresponded to NE expectations. In spite of these overall patterns, the extent to which levels and patterns of polymorphism within subpopulations differed from species-wide samples and NE expectations depended strongly on the geographic region (Jalisco vs. Balsas) from which subpopulations were sampled. This may be due to the demographic history of subpopulations in those regions. Overall, these results suggest that explicitly accounting for population structure may be important for studies examining the genetic basis of ecologically and agronomically important traits as well as for identifying loci that have been the targets of selection.     

Genetic consequences of polygyny and social structure in an Indian fruit bat, Cynopterus sphinx. I. Inbreeding, outbreeding, and population subdivision

Population subdivision into behaviorally cohesive kin groups influences rates of inbreeding and genetic drift and has important implications for the evolution of social behavior. Here we report the results of a study designed to test the hypothesis that harem social structure promotes inbreeding and genetic subdivision in a population with overlapping generations. Genetic consequences of harem social structure were investigated in a natural population of a highly polygynous fruit bat, Cynopterus sphinx (Chiroptera: Pteropodidae), in western India. The partitioning of genetic variance within and among breeding groups was assessed using 10-locus microsatellite genotypes for 431 individually marked bats. Genetic analysis of the C. sphinx study population was integrated with field data on demography and social structure to determine the specific ways in which mating, dispersal, and new social group formation influenced population genetic structure. Microsatellite data revealed striking contrasts in genetic structure between consecutive offspring cohorts and between generations. Relative to the 1998 (dry-season) offspring cohort, the 1997 (wet-season) cohort was characterized by a more extensive degree of within-group heterozygote excess (F(IS) = -0.164 vs. -0.050), a greater degree of among-group subdivision (F(ST) = 0.123 vs. 0.008), and higher average within-group relatedness (r = 0.251 vs. 0.017). Differences in genetic structure between the two offspring cohorts were attributable to seasonal differences in the number and proportional representation of male parents. Relative to adult age-classes, offspring cohorts were characterized by more extensive departures from allelic and genotypic equilibria and a greater degree of genetic subdivision. Generational differences in F-statistics indicated that genetic structuring of offspring cohorts was randomized by natal dispersal prior to recruitment into the breeding population. Low relatedness among harem females (r = 0.002-0.005) was primarily attributable to high rates of natal dispersal and low rates of juvenile survivorship. Kin selection is therefore an unlikely explanation for the formation and maintenance of behaviorally cohesive breeding groups in this highly social mammal.     

Genetic profile of cosmopolitan populations: effects of hidden subdivision

Natural populations of many organisms exhibit excess of rare alleles in comparison with the predictions of the neutral mutation hypothesis. It has been shown before that either a population bottleneck or the presence of slightly deleterious mutations can explain this phenomenon. A third explanation is presented in this work, showing that hidden subdivision within a population can also lead to an excess of rare alleles in the total population when the expectations of the neutral model are based on the allele frequency profile of the entire population data. With two examples (mitochondrial DNA-morph distribution and isozyme allele frequency distributions), it is shown that most cosmopolitan human populations exhibit excess of rare as well as total allele counts, when these are compared with the expectations of the neutral mutation hypothesis. The mitochondrial data demonstrate that such excesses can be detected from genetic variation at a single locus as well, and this is not due to stochastic error of allele frequency distributions. Contrast of the present observations with the allele frequency profiles in agglomerated tribal populations from South and Central America shows that even when the neutral expectations hold for individual subpopulations, if all subpopulations are grouped into a single population, the pooled data exhibit an excess of total number of alleles that is mainly due to the excess of rare alleles. Therefore, a primary cause of the excess number of rare alleles could be the hidden subdivision, and the magnitude of the excess indicates the extent of substructuring. The two components of hidden subdivision are: 1) Number of subpopulations, and 2) the average genetic distance among them. The implications of this observation in estimating mutation rate are discussed indicating the difficulties of comparing mutation rates from different population surveys.     

STORM: software for testing hypotheses of relatedness and mating patterns

Storm is a software package that allows users to test a variety of hypotheses regarding patterns of relatedness and patterns of mate choice and/or mate compatibility within a population. These functions are based on four main calculations that can be conducted either independently or in the hypothesis-testing framework: internal relatedness; homozygosity by loci; pairwise relatedness; and a new metric called allele inheritance, which calculates the proportion of loci at which an offspring inherits a paternal allele different from that inherited from its mother. STORM allows users to test four hypotheses based on these calculations and Monte Carlo simulations: (i) are individuals within observed associations or groupings more/less related than expected; (ii) do observed offspring have more/less genetic variability (based on internal relatedness or homozygosity by loci) than expected from the gene pool; (iii) are observed mating pairs more/less related than expected if mating is random with respect to relatedness; and (iv) do observed offspring inherit paternal alleles different from those inherited from the mother more/less often than expected based on Mendelian inheritance.     

Analysis of marital structure in Massachusetts using repeating pairs of surnames.

Analysis of surnames from marriages is now a well-established method in the study of marital and genetic structure. Traditional methods of partitioning inbreeding into random and nonrandom components rely on the total number of isonymous marriages. Because this number is often low, standard errors of inbreeding estimates tend to be high. Lasker and Kaplan (1985) devised a method that circumvents this problem by focusing on the total number of repeating pairs (RP) of surnames among marriages. The observed value of RP can be compared with the value expected at random (RPr) to assess patterns of subdivision within a population. The RP method is applied here to data from 3431 marriages that took place from 1800 to 1849 in 4 Massachusetts towns. The level of excess RP [(RP-RPr)/RPr] is positively associated with population size and exogamy rate. These results indicate a tendency for greater relative subdivision in larger, more exogamous populations. One possible reason for increased subdivision is preferential marriage by social class, although adequate data are not available for a test of this hypothesis.     

Population genetics and relatedness in a critically endangered island raptor, Ridgway’s Hawk Buteo ridgwayi

Many island avian populations are of conservation interest because they have a higher risk of extinction than mainland populations. Susceptibility of island birds to extinction is primarily related to human induced change through habitat loss, persecution, and introduction of exotic species, in combination with genetic factors. We used microsatellite profiles from 11 loci to assess genetic diversity and relatedness in the critically endangered hawk Buteo ridgwayi endemic to the island of Hispaniola in the Caribbean. Using samples collected between 2005 and 2009, our results revealed a relatively high level of heterozygosity, evidence of a recent genetic bottleneck, and the occurrence of inbreeding within the population. Pair relatedness analysis found 4 of 7 sampled breeding pairs to be related similar to that of first cousin or greater. Pedigree estimates indicated that up to 18 % of potential pairings would be between individuals with relatedness values similar to that of half-sibling. We discuss our findings in the context of conservation genetic management suggesting both carefully managed translocations and the initiation of a captive population as a safeguard of the remaining genetic diversity.     

Repetition of the same pairs of names in marriages in Fogo Island,Newfoundland, and genetic variation

Surnames can be used to investigate the genetic structure of human populations. The repeated-pairs approach (RP) uses information on the repetition of the same pairs of names from marital data sets to indicate the influence of clanlike behavior on mate choice. RP estimates the subdivision of the population into subgroups that breed among themselves and the kind of inbreeding ascribed to the Wahlund effect. The application of this method to the Fogo Island, Newfoundland, data set indicates a large percentage excess of observed over random repetitions. Presumably at least a part of this excess in the island population is due to denominational subdivision and endogamy. The expected relationship between RP and heterozygosity is observed in this case. Given these results, it would seem worthwhile to explore the relationship between the frequency of repetitions of pairs of surnames in marriages and genetic variation in other populations as well and, at the same time, to begin to inquire into the validity of such assumptions concerning the use of surnames as the monophyletic origin of common names.     

On the robustness of regular systems of inbreeding.

Half-sib, first cousin, half nth cousin, and nth cousin mating systems are robust in that small deviations from the mating structure will not significantly alter the levels of genetic identity or effect qualitative distinctions between the models. Substitution of nearest kin in matings may either increase or decrease the level of homozygosity depending on the mating structure; the effect of a single error in the mating structure is not an accurate indicator of the equilibrium resulting from recurrent errors. Models of mixed half nth cousin or nth cousin mating show that the relative frequency of the lowest order inbred mating essentially determines the level of homozygosity. Any positive relative frequency of more distant matings will reduce the probability of identity by descent under half-sib or first cousin mating to less than 1.     

Temporal and Spatial Heterogeneity of Mtdna Polymorphisms in Natural Populations of Drosophila mercatorum

Restriction endonuclease analysis of mtDNA was used to examine the genetic relatedness of several geographically separated isolines of the Drosophila mercatorum subgroup. In addition, we examined the temporal and spatial distribution of two mtDNA restriction site polymorphisms produced by the enzymes BstEII and BstNI at a single locality--Kamuela, Hawaii. Due to small sample sizes of some collections and the undesirable dependance of the estimation of polymorphism frequency on its variance, an arcsin square root transformation of the frequency data was used. We also use an Fst estimator of our transformed frequencies to demonstrate considerable spatial and temporal differentiation within the Kamuela population. In contrast, isozyme data from the same population reveals no pattern of differentiation. The temporal and geographic heterogeneity and population subdivision detected with mtDNA analysis also is consistent with the known dispersal behavior and ecological constraints of this species. The mtDNA data in conjunction with the isozyme data show that the population structure of the Kamuela D. mercatorum is close to the boundary line separating panmixia from subdivision, a conclusion that could not be made from isozyme data alone.     

Estimating relatedness in social groups

Genetic relatedness is a vital parameter in the evolution of social behaviour by kin selection. It can be easily estimated using genetic markers and calculating the genotypic correlation or regression of group members. Spatial gene frequency differentiation, due to population subdivision or isolation by distance, boosts the relatedness estimates. In such cases it may be useful to partition the estimate into components, the operational relatedness is normally that among individuals in social groups within the same subpopulation. Although it is straightforward to estimate the average relatedness in social groups, estimating values for specific individuals with the help of genetic markers is still problematic. Current estimators tend to give biased values and the sampling error is large. In spite of these shortcomings, studies of social behaviour combining relatedness and reproductive success are sorely needed.     

Genetic analysis of the population structure of socially organized oystercatchers (Haematopus ostralegus) using microsatellites

On the island of Schiermonnikoog (The Netherlands), the breeding population of oystercatchers can be divided into two groups: 'residents' and 'leapfrogs', based on their distinct social characteristics and limited probabilities of status change between breeding seasons. In order to investigate whether this social organization has caused local genetic differentiation, leapfrogs and residents were compared at eight polymorphic microsatellite loci. No significant genetic subdivision between residents and leapfrogs was observed (theta = 0.0000; 95% confidence interval (CI), -0.0027-0.0033), indicating that the oystercatcher population on the island of Schiermonnikoog has to be considered as one panmictic unit. Investigation of three additional locations in the northern part of The Netherlands did not reveal significant genetic population subdivision either (theta = -0.0005; 95% CI, -0.0045-0.0037), despite the fact that adult osytercatchers show extreme fidelity to their breeding localities. These results indicate panmixis and considerable levels of gene flow within the northern part of The Netherlands. Thus, the results from genetical analyses do not seem to be in agreement with observational data on the dispersal behaviour of breeding individuals. It is argued that the lack of population structure, locally on Schiermonnikoog as well as across larger geographical distances, is to be attributed to high levels of gene flow through dispersal of juvenile birds.     

Model-free Estimation of Recent Genetic Relatedness

Genealogical inference from genetic data is essential for a variety of applications in human genetics. In genome-wide and sequencing association studies, for example, accurate inference on both recent genetic relatedness, such as family structure, and more distant genetic relatedness, such as population structure, is necessary for protection against spurious associations. Distinguishing familial relatedness from population structure with genotype data, however, is difficult because both manifest as genetic similarity through the sharing of alleles. Existing approaches for inference on recent genetic relatedness have limitations in the presence of population structure, where they either (1) make strong and simplifying assumptions about population structure, which are often untenable, or (2) require correct specification of and appropriate reference population panels for the ancestries in the sample, which might be unknown or not well defined. Here, we propose PC-Relate, a model-free approach for estimating commonly used measures of recent genetic relatedness, such as kinship coefficients and IBD sharing probabilities, in the presence of unspecified structure. PC-Relate uses principal components calculated from genome-screen data to partition genetic correlations among sampled individuals due to the sharing of recent ancestors and more distant common ancestry into two separate components, without requiring specification of the ancestral populations or reference population panels. In simulation studies with population structure, including admixture, we demonstrate that PC-Relate provides accurate estimates of genetic relatedness and improved relationship classification over widely used approaches. We further demonstrate the utility of PC-Relate in applications to three ancestrally diverse samples that vary in both size and genealogical complexity.     

Low genome‐wide homozygosity in 11 Spanish ovine breeds

The population of Spanish sheep has decreased from 24 to 15 million heads in the last 75 years due to multiple social and economic factors. Such a demographic reduction might have caused an increase in homozygosity and inbreeding, thus limiting the viability of local breeds with excellent adaptations to harsh ecosystems. The main goal of our study was to investigate the homozygosity patterns of 11 Spanish ovine breeds and to elucidate the relationship of these Spanish breeds with reference populations from Europe, Africa and the Near East. By using Ovine SNP50 BeadChip data retrieved from previous publications, we have found that the majority of studied Spanish ovine breeds have close genetic relatedness with other European populations; the one exception is the Canaria de Pelo breed, which is similar to North African breeds. Our analysis has also demonstrated that, with few exceptions, the genomes of Spanish sheep harbor fewer than 50 runs of homozygosity (ROH) with a total length of less than 350 Mb. Moreover, the frequencies of very long ROH (>30 Mb) are very low, and the inbreeding coefficients (F_ROH) are generally small (F_ROH < 0.10), ranging from 0.008 (Rasa Aragonesa) to 0.086 (Canaria de Pelo). The low levels of homozygosity observed in the 11 Spanish sheep under analysis might be due to their extensive management and the high number of small to medium farms.     

Balancing in- and out-breeding by the predatory mite <Emphasis Type="Italic">Phytoseiulus persimilis</Emphasis>

In- and out-breeding depressions are commonly observed phenomena in sexually reproducing organisms with a patchy distribution pattern, and spatial segmentation and/or isolation of groups. At the genetic level, inbreeding depression is due to increased homozygosity, whereas outbreeding depression is due to inferior genetic compatibility of mates. Optimal outbreeding theory suggests that intermediate levels of mate relatedness should provide for the highest fitness gains. Here, we assessed the fitness consequences of genetic relatedness between mates in plant-inhabiting predatory mites Phytoseiulus persimilis, which are obligatory sexually reproducing but haplo-diploid. Both females and males arise from fertilized eggs but males lose the paternal chromosome set during embryogenesis, dubbed pseudo-arrhenotoky. Phytoseiulus persimilis are highly efficacious in reducing crop-damaging spider mite populations and widely used in biological control. Using iso-female lines of two populations, from Sicily and Greece, we assessed the fecundity of females, and sex ratio of their offspring, that mated with either a sibling, a male from the same population or a male from the other population. Additionally, we recorded mating latency and duration. Females mating with a male from the same population produced more eggs, with a lower female bias, over a longer time than females mating with a sibling or with a male from the other population. Mating latency was unaffected by mate relatedness; mating duration was disproportionally long in sibling couples, likely indicating female reluctance to mate and sub-optimal spermatophore transfer. Our study provides a rare example of in- and out-breeding depression in a haplo-diploid arthropod, supporting the optimal outbreeding theory.     

Interpopulation variation in inbreeding is primarily driven by tolerance of mating with relatives in a spermcasting invertebrate

The degree to which individuals inbreed is a fundamental aspect of population biology shaped by both passive and active processes. Yet, the relative influences of random and non-random mating on the overall magnitude of inbreeding are not well characterized for many taxa. We quantified variation in inbreeding among qualitatively accessible and isolated populations of a sessile marine invertebrate (the colonial ascidian Lissoclinum verrilli) in which hermaphroditic colonies cast sperm into the water column for subsequent uptake and internal fertilization. We compared estimates of inbreeding to simulations predicting random mating within sites to evaluate if levels of inbreeding were (1) less than expected because of active attempts to limit inbreeding, (2) as predicted by genetic subdivision and passive inbreeding tolerance, or (3) greater than simulations due to active attempts to promote inbreeding via self-fertilization or a preference for related mates. We found evidence of restricted gene flow and significant differences in the genetic diversity of L. verrilli colonies among sites, indicating that on average colonies were weakly related in accessible locations, but their levels of relatedness matched that of first cousins or half-siblings on isolated substrates. Irrespective of population size, progeny arrays revealed variation in the magnitude of inbreeding across sites that tracked with the mean relatedness of conspecifics. Biparental reproduction was confirmed in most offspring (86%) and estimates of total inbreeding largely overlapped with simulations of random mating, suggesting that interpopulation variation in mother–offspring resemblance was primarily due to genetic subdivision and passive tolerance of related mates. Our results highlight the influence of demographic isolation on the genetic composition of populations, and support theory predicting that tolerance of biparental inbreeding, even when mates are closely related, may be favoured under a broad set of ecological and evolutionary conditions.     

Genetic Structure of the Veined Rapa Whelk (Rapana venosa) Populations Along the Coast of China

The veined rapa whelk (Rapana venosa) is a valuable and important fishery resource in China. In order to provide guidelines for fisheries management, the genetics population structure was assessed using 10 polymorphic allozyme loci from seven populations of R. venosa across the species' range on the Chinese coast. The mean allele richness ranged from 1.65 to 2.13, and the mean heterozygosity ranged from 0.086 to 0.149. Significant genetic differentiation was present, and the theta value was 0.016 across all populations. The nMDS plot of pairwise theta values, UPGMA dendrogram, and AMOVA analysis were in good agreement and identified three geographic subdivision groups. We propose that the genetic structure may be due to larval dispersal barriers and localized adaptation or a combination thereof. The implications of these results for the management of R. venosa populations are discussed.     

spip 1.0: a program for simulating pedigrees and genetic data in age‐structured populations

We present the program spip for simulating multilocus genetic data on individuals in age‐structured populations. In addition to genetic data on sampled individuals, the pedigree connecting all individuals in the population is recorded. This allows investigation of the relationship between family structure and population parameters. We foresee that spip will be useful for evaluating multilocus estimators of pairwise relatedness and population structure, and for simulating the distribution of relatedness in populations with varying demographies. It also provides a method for simulating genetic drift in complex populations.     

Comparative genetic structure reflects underlying life histories of three sympatric small mammal species in continuous forest of south-eastern Australia

There are strong links between habitat fragmentation, population size and genetic structure. However, to fully understand the long‐term effects of fragmentation on population viability, it is necessary to first understand the relationship between life history traits and genetic characteristics in un‐fragmented habitats. This is best done through comparing patterns of genetic diversity in sympatric species, since relative data may be maximally informative. We compared genetic structure and diversity among three sympatric small mammals – of which two were marsupial species – over a small spatial scale (<4 km) using microsatellite genotypic data from individuals sampled from four grids spaced in a linear fashion. Expected heterozygosity was high for all three species (average H_e range: 0.781–0.886), but the smallest species had significantly higher genetic diversity (both heterozygosity and allelic diversity) than the two larger species, possibly due to it having; 1) large effective population size and little overlap in generations, and/or 2) high mutation rates in large alleles. Genotypic isolation by distance (measured as relatedness) was detected only at the smallest distance of 750 m, for any species. In the two marsupial species F_st, relatedness and percentage of individuals assigned to site of capture were significantly lower in males than in females, corroborating previous reports of male‐biased dispersal. However, sex‐biased dispersal was not detected for the third species (the native bush rat Rattus fuscipes), and presence of significant heterozygous deficits suggested this resulted in inbreeding within kin‐structured demes. We speculate that habitat fragmentation will have differing effects on population dynamics, social organisation and mating systems for the two marsupial species compared to the native murid rodent, due to their differing population structure and divergent inbreeding avoidance mechanisms.