Department of Biostatistics, Columbia University, New York NY 10032, USA
Received date: July 08, 2011; Accepted date: November 23, 2011; Published date: November 30, 2011
Citation: RoyChoudhury A (2011) Loss of Rare Recessive Lethal Alleles due to Consanguineous Marriages in Inbred Populations. J Biom Biostat 2:121. doi:10.4172/2155-6180.1000121
Copyright: © 2011 RoyChoudhury A. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Visit for more related articles at Journal of Biometrics & Biostatistics
In a population practicing consanguineous marriage, rare recessive lethal alleles (RRLA) have higher chances of affecting phenotypes. As inbreeding causes more homozygosity and subsequently more deaths, the loss of individuals with RRLA decreases the frequency of these alleles. Although this phenomenon is well studied in general, here some hitherto unstudied cases are presented. An analytical formula for the RRLA frequency is presented for infinite monoecious population practicing several different types of inbreeding. In finite diecious populations, it is found that more severe inbreeding leads to quicker RRLA losses, making the upcoming generations healthier. A population of size 10,000 practicing 30% half-sib marriages loses more than 95% of its RRLA in 100 generations; a population practicing 30% cousin marriages loses about 75% of its RRLA. Our findings also suggest that given enough resources to grow, a small inbred population will be able to rebound while losing the RRLA.
Inbreeding; Consanguineous marriage; lethal allele; Deleterious allele; Inbreeding depression; Sib-marriage
Consanguineous marriage (CM), the marriage between two related individuals, is commonly practiced in many populations. At least 20% of human populations live in cultures with preference for such marriages . A birth study in Birmingham, UK found that 57% of all couples in British Pakistani population were first cousins . Half-sib marriages were allowed in ancient Greece, whereas full sib marriages were fairly common in Roman Egypt . About one in six marriages in first-third century CE Roman Egypt are recorded to be between full sibs . On the other hand, non-deliberate inbreeding may occur in a population with small size. Inbreeding resulting from CM is known to increase the chances of still-births and infant deaths [4,5,6].
It is usually held that inbreeding does not cause changes in the allele frequency. This is indeed true when one compares the parent population and the gamete population before any deaths have occurred in the latter. However, inbreeding increases number of homozygous individuals (with two copies of a recessive lethal allele), and as a result the homozygous individuals die due to selection. As a result, the frequency of such alleles reduces among the child generation (see, for example, [7,8,9]). It is known that such alleles are present in many (human or wild) populations . Here we study the amount of loss of rare recessive lethal alleles (RRLA) under different types of inbreeding, in an infinite monoecious population and in finite diecious populations with various parameters.
Infinite monoecious population
Let us consider an infinite monoecious population with discrete generations. Also assume that at a certain locus there are two possible alleles “A” and “a” and the effect of mutation is negligible. Suppose that “a” is a rare recessive lethal allele. Let pt (1 > pt > 0) be the “a”-frequency (frequency of “a”) at Generation t after any deaths have occurred due to selection. Also assume that the population practices k different types of inbreeding I1, I2,..., Ik with inbreeding coefficients γ1, γ2,..., γk respectively. A fraction αi (> 0) of the matings are of type Ii. (The remaining fraction of matings are random.) In inbreeding type , the most recent common ancestor is ti generations before the generation of the child. Thus, the probability that a child (at Generation 0) of Ii will be "aa" is
If all the children in Generation 0 had survived, then the frequency of "aa" individuals would have been:
b0 = Pr ("aa") = Pr ("aa" | parents' mating was random) Pr(Parents' mating was random)
Pr (“aa” | parents’ mating was of type li) Pr (Parents’ mating was oftype li)
= Pr (“aa” | parents’ mating was random)
(“aa”| parents’ mating was of typeli)αi
Noting that αi > 0, 1 > γi > 0 and , it follows that and subsequently both sides of Eq. (2) is positive. The "a"-frequency before any death (i.e. in the gamete population) is P-1. This is because there is no change in allele frequency between the parent and gamete populations. As "aa" individuals die, each taking away two "a" alleles
As b0 > 0, it follows that p-1 > p0, and thus the allele-frequency has decreased as a result of inbreeding. Note that even without any deliberate inbreeding there will be some homozygosity, resulting in death due to selection.
If inbreeding continues, the frequency of lethal alleles keeps decreasing. To see the long-term effect of this in finite diecious populations with mutation, we use computer simulations.
Simulation: finite diecious population
We simulated from a version of Wright-Fisher model (). The populations have discrete generations with fixed generation size: N (even number) individuals with N/2 couples (N/2 females and N/2 males). Each individual has two alleles at L independently segregating loci. The rare recessive lethal alleles (RRLA) are present among these L loci with a predetermined initial frequency. (The allele frequency is allowed to vary across the L loci, but the fraction among the total L loci is fixed.) At each new generation, each child is produced by first randomly selecting a parental couple from the previous generation, and then taking one (randomly selected) allele from each locus of each parent. We also force a pre-determined fraction of parental couples to be randomly selected from a set of related couples (sibs, half-sibs or cousins, according to predetermined parameters) so that their children are inbred. However, if at least one parent has two copies of the lethal allele (at any locus) then we do not generate a child from them and instead select another couple with same degree of relatedness. Once the children are generated, we mutate a very small fraction (1 in 108) of their healthy alleles to lethal alleles to simulate the effect of mutation. This is done by changing a healthy allele of a random individual to lethal.
We first compare the recessive lethal allele-frequencies of four simulated populations with no deliberate inbreeding (Figure 1). One of the populations has a 2% growth rate per generation. Rest have constant sizes. The RRLA frequencies in populations with N=10,000 and 1,000 go down by 50% and 70% respectively. There is more rapid decline in RRLA frequency in the populations with N=100. The decline is more in the population with constant size as it loses all its RRLA within 70 generations. However, even the growing population lost more than 97% percent of its RRLA in 100 generations.
Then we compare the frequencies in populations with constant size of 10,000 and different type of inbreeding and different initial frequency of RRLA (Figure 2). For the initial frequency of 0.01%, the populations with sib and half-sib matings lost the RRLA completely within 100 generations. Also, in the population with 1% initial frequency and 50% half-sib matings, the alleles were lost within 100 generations. The rest, although did not lose RRLA completely, showed significant decline. The population with 30% half-sib marriages loses more than 95% of its RRLA; the population with 30% cousin marriages loses about 75% of its RRLA; however, the population with 10% half-sib marriages retains about two-third of its RRLA.
Next we compare the proportion of individual with lethal traits across populations with constant size of 10,000, initial allele frequency of 1% and different type of inbreeding (Figure 3). The proportion of individuals with lethal traits (that is, "aa" genotype before any deaths) declines most rapidly for sib-matings, and little less rapidly for half-sib matings. The population with 30% full-sib marriages stops having thelethal traits within 95 generations. This is because the RRLA frequency in this population is so low at that point, that it is unlikely to have a homozygous individual.
Inbreeding alters the genotype frequencies, which results in removal of lethal alleles in homozygous individuals due to selection. Based on this basic fact, we presented three major findings. First, even a very small population has chances of surviving the effects of inbreeding and rebounding. As Figure 1 shows, both the populations that started with size N = 100 (one constant size and the other growing) lost the RRLA rapidly. In the first few generations there were many individuals with lethal traits, which facilitated the loss of the RRLA, and made the upcoming generations healthier. The high level of RRLA did not stop the population from growing in size either. Thus, small endangered wild populations that are at high risk of inbreeding has the possibility of rebounding.
The second major result focuses on populations that actively practices inbreeding. As shown in (Figures 2, 3), more severe the inbreeding, quicker the loss of RRLA. That is, the populations with higher rate of inbreeding and with more acute type of inbreeding will lose the RRLA more rapidly. We have also produced the rate of RRLA loss for different scenarios. One of our simulations shows that a population of size 10,000 practicing 30% cousin marriages loses about 75% of its RRLA in 100 generations. This is significant as there are many populations practicing 30% or more cousin marriages. (Figure 3) shows that the proportion of individuals with lethal traits declines (or disappears) after a steady period of inbreeding. A similar result is also shown with populations practicing half-sib marriages.
In addition to the above we have also provided a formula for RRLA-frequency in an infinite monoecious population practicing several different types of inbreeding. The formula is given as a function of the RRLA-frequencies of the previous generations. The formula also depends on the inbreeding coefficients of the types of inbreeding practiced, as well as on their frequencies (which are assumed to be constant across generations). Note that this formula may be used to approximate the RRLA frequency in large diecious populations; this formula is especially useful in population where several different types of inbreeding are practiced.
The author thanks "Statistical Methods in Genetic Epidemiology Journal Club" group at the Department of Biostatistics at Columbia University for helpful discussions.