Isozyme polymorphisms maintained by lethal loci in inbred strains of Aedes triseriatus

We analyzed the progeny of brother-sister matings of two inbred strains of the treehole mosquito, Aedes triseriatus, at several isozyme loci to determine which factors were maintaining heterozygosity at these loci after prolonged inbreeding. The triseriatus Kramer (TK) lines were polymorphic for the Odh and Hbd loci after 12 generations of full-sib mating, and the triseriatus Vero Beach (TV) lines were polymorphic for the Odh, Hk-4, Pgm, and Hbd loci after 14 generations of full-sib mating. Genetic analysis of 22 TK F13 matings and nine TV F15 matings showed that heterozygosity was enforced by lethal loci to which all the polymorphic isozyme loci were linked. The lethal loci formed a balanced lethal system. We determined the relative positions and map distances of lethals and enzyme loci.     

Isoenzyme variation in Aedes aegypti correlated with Dirofilaria immitis infectability

From the Vero Beach strain of the mosquito Aedes (Stegomyia) aegypti (L.) (Diptera: Culicidae), substrains were selected for susceptibility (SS) and refractoriness (RR) to the dog heartworm Dirofilaria immitis (Leidy) (Filarioidea: Onchocercidae). These two lines and their reciprocal F1 hybrids were analysed for genetic variation at 14 enzyme loci, using polyacrylamide gel electrophoresis. Six of the enzyme loci showed variation (sample size 48 alleles/locus/line). Three of these were monomorphic in the refractory line but polymorphic in the susceptible, i.e. aconitase hydratase (Acoh), isocitrate dehydrogenase-1 (Idh-1) and phosphoglucomutase (Pgm). The other three loci, glucose-6-phosphate isomerase (Gpi), hexokinase-1 (Hk-1) and isocitrate dehydrogenase-2 (Idh-2), were polymorphic in both SS and RR lines and their hybrids. At two loci (Hk-1, Pgm) three alleles were detected, whereas the other polymorphic loci had only two alleles. For Hk-1, the most frequent allele was Hk-1(80) (0.563) in refractory and Hk-1(100) in the susceptible (0.521) and F1 hybrids. For Pgm the most frequent alleles were Pgm125 in the susceptible line (0.646) and Pgm100 in the F1 hybrids (0.563 and 0.604) and refractory line (1.000). The mean observed heterozygosity (Ho), the mean Hardy-Weinberg expected heterozygosity (He) and the mean number of alleles per locus in the refractory line were lower, but not significantly so, than in the susceptible line and their reciprocal F1 hybrids; the proportion of polymorphic loci was significantly lower in the refractory than in the susceptible line and their F1 hybrids. Within both lines all polymorphisms were in Hardy-Weinberg equilibrium, whereas significant departures from predicted frequencies were observed in SS x RR hybrids at four polymorphic loci (Acoh, Gpi, Hk-1, Pgm) and at three polymorphic loci (Acoh, Hk-1, Pgm) in RR x SS hybrids. The average Nei's and modified Rogers' genetic distances between the lines were 0.024 and 0.139, respectively. These electrophoretic data show that the refractory line (putatively lacking fi allele) can be distinguished from the susceptible line (fi/fi) and their hybrids (heterozygous fi) by isozyme marker frequencies, but it remains to be seen whether this difference is causal or chance linkage. In any case, this model system of Ae. aegypti/D. immitis provides opportunities to better understand and manipulate the molecular biology of filariasis transmission.     

Analysis of inbreeding in a colonizing population of Drosophila subobscura

An analysis of the effects of inbreeding on the genetic structure of a colonizing population of Drosophila subobscura has been carried out. Species of Drosophila, particularly D. subobscura, may have lethal alleles associated with chromosomal inversions and our aim was to assess the extent to which the genome is balanced in this way. The frequencies of chromosomal inversions were compared between a large population and a set of 72 lines that were maintained by brother-sister mating for 10 generations. Fisher's matrix method was used to calculate the expected homozygosity in these inbred lines for 5 allozyme loci (Aph, Hk-1, Lap, Odh, and Pept-1) used as markers of large chromosomal segments. Furthermore, the expected rates of fixation corresponding to these allozyme loci were also calculated. The results show that the amount of homozygosis observed did not differ significantly from expectations (with the corresponding loss of lines as a consequence of the reduction in viability). However, two deviations from strict neutrality were observed: there was a heterozygote excess at the Lap locus, and the frequency of the O ₅ inversion (always associated with a lethal gene in colonizing populations) was higher than expected.     

Effect of rate of inbreeding on inbreeding depression in Drosophila melanogaster

Summary This experiment was designed to study the relationship between rate of inbreeding and observed inbreeding depression of larval viability, adult fecundity and cold shock mortality in Drosophila melanogaster. Rates of inbreeding used were full-sib mating and closed lines of N=4 and N=20. Eight generations of mating in the N=20 lines, three generations in the N=4 lines and one generation of full-sib mating were synchronised to simultaneously produce individuals with an expected level of inbreeding coefficient (F) of approximately 0.25. Inbreeding depression for the three traits was significant at F=0.25. N=20 lines showed significantly less inbreeding depression than full-sib mated lines for larval viability at approximately the same level of F. A similar trend was observed for fecundity. No effect of rate of inbreeding depression was found for cold shock mortality, but this trait was measured with less precision than the other two. Natural selection acting on loci influencing larval viability and fecundity during the process of inbreeding could explain these results. Selection is expected to be more effective with slow rates of inbreeding because there are more generations and greater opportunity for selection to act before F=0.25 is reached. Selection intensities seem to have been different in the three traits measured. Selection was most intense for larval viability, less intense for fecundity and, perhaps, negligible at loci influencing cold shock mortality.     

Variable outcrossing and the genetic structure of predominantly self-pollinated species

Several studies of highly inbreeding plant species have revealed appreciable levels of heterozygosity, exceeding those anticipated from the mating system. Hence net heterozygote advantage has been postulated. An alternative proposal for the observed heterozygous excess is that it arises from fluctuating outcrossing rates. This study analyses the effects of variable outcrossing on heterozygosity and on the expected inbreeding, using simple models and simulations. Generally, estimates of observed heterozygosity are not seriously biased. However, estimates of the predicted fixation index, estimated from the current outcrossing rate, were biased upward, especially when both over-dominant selection operates and the distribution of outcrossing is markedly skewed. Such an upward bias would increase the difference between observed and expected heterozygosity. Even so, only a small portion of the observed excess heterozygosity in inbreeders is likely to be due to variable outcrossing. The bias can be effectively reduced by a few independent measurements of outcrossing in different generations.     

中华稻蝗四个种群的遗传分化

采用水平淀粉凝胶电泳技术研究中华稻蝗(Oxya chinensis)四个种群的12个基因座位,探讨其遗传分化.这四个种群分别采自内蒙古呼和浩特、山西代县、山西太原和陕西西安.Ck和Mdh-2在四个种群中均为单态,其余的基因座位至少在一个种群内有两个以上的等位基因;在Ldh和Mdh-1的等位基因频率呈现梯度分布趋势.多态基因座位百分率(P)和平均每个基因座位的等位基因数目(A)分别为58.3%-66.7%和2.2-2.8,平均杂合度为Ho=0.173-0.240.除Gpi,Hk-2,Idh,Ldh和Mdh-1在部分蝗虫种群符合Hardy-Weinberg平衡外,其余绝大多数基因座位的基因型频率显著偏离Hardy-Weinberg平衡.四个种群间FST平均值不存在显著差异(FST=0.0510,P＞0.05),结合高的Nei's遗传一致度(I＞0.97)可知,种群之间遗传分化不明显.我们认为人类的农业活动可能促进了种群间的基因交流,从而降低了分化程度[动物学报50(2)187-192,2004].     

Microsatellite analysis of genetic diversity of the Vietnamese sika deer (Cervus nippon pseudaxis)

The Vietnamese sika deer (Cervus nippon pseudaxis) is an endangered subspecies of economic and traditional value in Vietnam. Most living individuals are held in traditional farms in central Vietnam, others being found in zoos around the world. Here we study the neutral genetic diversity and population structure of this subspecies using nine microsatellite loci in order to evaluate the consequences of the limited number of individuals from which this population was initiated and of the breeding practices (i.e., possible inbreeding). Two hundred individuals were sampled from several villages. Our data show both evidence for limited local inbreeding and isolation by distance with a mean F(ST) value of 0.02 between villages. This suggests that exchange of animals occurs at a local scale, at a rate such that highly inbred mating is avoided. However, the genetic diversity, with an expected heterozygosity (H(e)) of 0.60 and mean number of alleles (k) of 5.7, was not significantly larger than that estimated from zoo populations of much smaller census size (17 animals sampled; H(e) = 0.65, k = 4.11). Our results also suggest that the Vietnamese population might have experienced a slight bottleneck. However, this population is sufficiently variable to constitute a source of individuals for reintroduction in the wild in Vietnam.     

Population genetic structure of the endangered tropical tree species Caryocar brasiliense, based on variability at microsatellite loci

We report the population genetic structure of the endangered tropical tree species Caryocar brasiliense, based on variability at 10 microsatellite loci. Additionally, we compare heterozygosity and inbreeding estimates for continuous and fragmented populations and discuss the consequences for conservation. For a total of 314 individuals over 10 populations, the number of alleles per locus ranged from 20 to 27 and expected and observed heterozygosity varied from 0.129 to 0.924 and 0.067 to 1.000, respectively. Significant values of theta and R(ST) showed important genetic differentiation among populations. theta was much lower than R(ST), suggesting that identity by state and identity by descent have diverged in these populations. Although a significant amount of inbreeding was found under the identity by descent model (f = 0.11), an estimate of inbreeding for microsatellite markers based on a more adequate stepwise mutation model showed no evidence of nonrandom mating (R(IS) = 0.04). Differentiation (pairwise F(ST)) was positively correlated with geographical distance, as expected under the isolation by distance model. No effect of fragmentation on heterozygosity or inbreeding could be detected. This is most likely due to the fact that Cerrado fragmentation is a relatively recent event (approximately 60 years) compared to the species life cycle. Also, the populations surveyed from both fragmented and disturbed areas were composed mainly of adult individuals, already present prior to ecosystem fragmentation. Adequate hypothesis testing of the effect of habitat fragmentation will require the recurrent analysis of juveniles across generations in both fragmented and nonfragmented areas.     

Analysis of the relationship between allozyme heterozygosity and fitness in the rare Gentiana pneumonanthe L.

Especially for rare species occurring in small populations, which are prone to loss of genetic variation and inbreeding, detailed knowledge of the relationship between heterozygosity and fitness is generally lacking. After reporting on allozyme variation and fitness in relation to population size in the rare plant Gentiana pneumonanthe, we present a more detailed analysis of the association between heterozygosity and individual fitness. The aim of this study was to test whether increased fitness of more heterozygous individuals is explained best by the ‘inbreeding’ hypothesis or by the ‘overdominance’ hypothesis. Individual fitness was measured during 8 months of growth in the greenhouse as the performance for six life-history parameters. PCA reduced these parameters to four main Fitness Components. Individual heterozygosity was scored for seven polymorphic allozyme loci. For some of these loci (e.g. Aat3, Pgm1 and 6Pgdh2) heterozygotes showed a significantly higher relative fitness than homozygotes. To test the inbreeding model, regression analyses were performed between each Fitness Component and the number of heterozygous loci per individual. Multiple regressions with the adaptive distance of five loci as independent variables were used to test the overdominance model. Only the inbreeding model was a statistically significant explanation for the relationship between heterozygosity and fitness in G. pneumonanthe. The number of heterozygous loci was significantly negatively correlated with the coefficients of variation of three of the six initially measured fitness parameters. This suggests a lower developmental stability among more homozygous plants and may explain the higher phenotypic variation in small populations of the species observed earlier. The importance of the results for conservation biology is discussed.     

Effects of population structures and selection strategies on the purging of inbreeding depression due to deleterious mutations

Stochastic simulations were run to compare the effects of nine breeding schemes, using full-sib mating, on the rate of purging of inbreeding depression due to mutations with equal deleterious effect on viability at unlinked loci in an outbred population. A number of full-sib mating lines were initiated from a large outbred population and maintained for 20 generations (if not extinct). Selection against deleterious mutations was allowed to occur within lines only, between lines or equal within and between lines, and surviving lines were either not crossed or crossed following every one or three generations of full-sib mating. The effectiveness of purging was indicated by the decreased number of lethal equivalents and the increased fitness of the purged population formed from crossing surviving lines after 20 generations under a given breeding scheme. The results show that the effectiveness of purging, the survival of the inbred lines and the inbreeding level attained are generally highest with between-line selection and lowest with within-line selection. Compared with no crossing, line crossing could lower the risk of extinction and the inbreeding coefficient of the purged population substantially with little loss of the effectiveness of purging. Compromising between the effectiveness of purging, and the risk of extinction and inbreeding coefficient, the breeding scheme with equal within- and between-line selection and crossing alternatively with full-sib mating is generally the most desirable scheme for purging deleterious mutations. Unless most deleterious mutations have relatively large effects on fitness in species with reproductive ability high enough to cope with the depressed fitness and thus increased risk of extinction with inbreeding, it is not justified to apply a breeding programme aimed at purging inbreeding depression by inbreeding and selection to a population of conservation concern.     

Larval growth, juvenile size and heterozygosity in laboratory reared mussels, Mytilus edulis

Studies with marine bivalve juveniles have shown a positive correlation between growth and allozyme multi-locus heterozygosity (MLH), and, in some cases, between larval growth and juvenile growth, but there has been little research on the relationship between allozyme heterozygosity and larval growth. Larvae of M. edulis from different mating systems (half-sib families with a single female, or a single male parent, a reciprocal cross of two malesxtwo females and two mass matings of 13x13 and 8x17 females and males, respectively) were reared in the laboratory and selected into fast and slow growing groups when about 10-30% were undergoing metamorphosis. Offspring were reared to the juvenile stage (>3.00 mm) and both groups of each mating were electrophoresed and genotyped at up to 12 allozyme loci. There was generally good agreement with Mendelian inheritance (half-sibs and reciprocal cross) or the Hardy-Weinberg model (mass matings). Null alleles were detected at the Odh and Lap loci but there was no evidence that null allele heterozygotes grew slower than other genotypes. Over all cohorts, juveniles from the fast growing larval group were not significantly larger, or smaller, than juveniles from the slow growing group which suggests that larval growth rate may be independent of juvenile growth rate. This observation agrees with some, but not all, earlier studies and has commercial relevance. Tests of heterozygosity and juvenile shell length indicated no association between average heterozygosity across all allozyme loci and the size of juveniles in any cohort regardless of the mating system used or their larval growth rate. The association between MLH and juvenile growth in bivalves is seldom detected in cohorts from a limited genetic background. The lack of an association between heterozygosity and size might therefore be expected in the half-sib and reciprocal cross cohorts, but not in the mass matings. The results argue against any significant association between heterozygosity and larval size in mussels.     

黄山钓桥青冈种群的交配系统与近交衰退

采用垂直板型不连续聚丙烯酰胺凝胶电泳检测了黄山钓桥青冈种群的遗传变异、交配系统及近交衰退程度。黄山钓桥青冈群维持有中等程度的遗传变异,多态位点百分比和期望杂合度分别为５０％和０．１８１４。采用ＰＯＤ－２、ＥＳＴ－１和ＥＳＴ－２等３个多态位点,利用ＭＬＴ程序估计的花粉库和母株基因频率之间存在一些差异,可能是由于成熟个体花粉产量不同和种群内个体非随机分布造成的。单位点杂交率平均为０．３９２,略高于多位     

Decline in Heterozygosity under Full-Sib and Double First-Cousin Inbreeding in Drosophila Melanogaster

The effects of inbreeding on heterozygosities and reproductive fitness were determined by carrying out full-sib and double first-cousin inbreeding in Drosophila melanogaster populations for up to 18 generations. Parents were scored each generation for five or six polymorphic enzyme loci, and progeny numbers per pair were recorded. Inbreeding depression, in the form of significant reductions in progeny numbers and significant extinction of lines, was observed. Heterozygosity decreased at a significantly slower rate than predicted, being about 80% of expected. The full-sib and double first-cousin treatments showed similar disagreement with expectations over comparable ranges of inbreeding. Natural selection was shown to favor heterozygotes in the inbred lines. Associative overdominance was the most probable explanation for the slower than expected decline in heterozygosity.     

Effects on heterozygosity and reproductive fitness of inbreeding with and without selection on fitness in Drosophila melanogaster

The effects of inbreeding, with (IS) and without selection (IO) for reproductive fitness, on inbreeding depression and heterozygosity were evaluated in 20 lines of each treatment inbred over seven generations using full-sib mating. The survival of lines was significantly greater in IS (20/20) than in IO (15/20). The competitive index measure of reproductive fitness was significantly lower in the inbred lines than in the outbred base population, but not significantly different in surviving IS and IO lines. There was a trend for higher fitness in the IS treatment as relative fitnesses were 19% higher in IS than IO for surviving lines and 59% higher for all lines. Heterozygosities were lower in the inbred lines than in the base population, and significantly higher in the IS than the IO lines. Consequently, the reduction of inbreeding depression in IS has been achieved, at least in part, by slowing the rate of fixation.     

新铁炮百合自交初代遗传分化的等位酶分析

采用等位酶分析技术,检测了新铁炮百合自交初代各株系的遗传多样性与遗传分化。结果表明,在人工辅助自交和选择作用下,随着世代的增加,多态位点的比率下降,每个位点平均等位基因数降低,观测杂合度和期望杂合度下降。由于选择作用,株系间纯合的位点和速度不同,固定的等位基因不同。F4代株系间平均遗传距离较F3、代株系有一定增大,加强了株系间的分化。     

Effects of Partial Inbreeding on Fixation Rates and Variation of Mutant Genes

Diffusion methods were used to investigate the fixation probability, average time until fixation and extinction, and cumulative heterozygosity and genetic variance for single mutant genes in finite populations with partial inbreeding. The critical parameters in the approximation are the coefficient of inbreeding due to nonrandom mating (F) and the effective population size (Ne), which also depends on F and the variance of family size. For large Ns, the fixation probability (u) is u = 2(Ne/N)s (F + h - Fh), where N is the population census, s is the coefficient of selection of the mutant homozygote and h is the coefficient of dominance. For Poisson family size (independent Poisson distributions of selfed and nonselfed offspring with partial selfing, and independent Poisson distributions of male and female numbers with partial sib mating), Ne = N/(1 + F), and the time until fixation is approximately equal to Ne/N times the time to fixation with random mating, but this relation does not hold, however, for other distributions of family size. The cumulative nonadditive variance until fixation or loss for dominant genes is reduced with increasing F while for recessive genes it is increased with intermediate values of F. The average time until extinction of deleterious mutations is reduced by increasing F. This reduction, when expressed as a proportion, is approximately independent of the initial gene frequency as well as the selective disadvantage if this is large.     

Major histocompatibility complex heterozygosity enhances reproductive success

We investigated how heterozygosity at the major histocompatibility complex (MHC) affects fitness in wild-derived (F2) house mice (Mus musculus musculus). To compare and control for potential confounding effects from close inbreeding and genome-wide heterozygosity, we used mice that were systematically outbred. We assessed how heterozygosity at MHC and background loci (using 15 microsatellite markers on 11 different chromosomes) affects individual survival and reproductive success (RS) in large, semi-natural population enclosures. We found that overall heterozygosity significantly increased RS, and this correlation was entirely explained by heterozygosity at two MHC loci. Moreover, we found that the effects of MHC heterozygosity depend on the level of background heterozygosity, and the benefits of maximal MHC heterozygosity show a curvilinear effect with increasing background heterozygosity. The enhanced RS from MHC heterozygosity was not because of increased survival, and although MHC heterozygosity was correlated with body mass, body mass did not correlate with RS when heterozygosity is controlled. Breeders were more MHC heterozygous than nonbreeders for both sexes, indicating that MHC heterozygosity enhanced fecundity, mating success or both. Our results show that (i) MHC heterozygosity enhances fitness among wild, outbred as well as congenic laboratory mice; (ii) heterozygosity-fitness correlations can potentially be explained by a few loci, such as MHC; (iii) MHC heterozygosity can increase fitness, even without affecting survival, by increasing mating and RS; and (iv) MHC effects depend on background genes, and maximal MHC heterozygosity is most beneficial at intermediate or optimal levels of background heterozygosity.     

Genetic diversity and the mating system of a rare Mexican pinon, Pinus pinceana, and a comparison with Pinus maximartinezii (Pinaceae)

Weeping pi?on (Pinus pinceana) has a restricted and fragmented range, trees are widely scattered within populations, and reproduction is limited. Nevertheless, genetic diversity was high; based on 27 isozyme loci in 18 enzyme systems, unbiased expected heterozygosity averaged 0.174. Differentiation also was high (F(ST) = 0.152), reflecting isolation between southern, central, and northern fragments of the range. Among populations in the northern fragment, F(ST) was only 0.056, and the number of migrants per generation (Nm) was 4.21, which should preclude fixation. Nm between central and southern populations or between them and populations in the northern fragment was lower, 0.99-1.66, indicating a degree of genetic isolation. Multilocus outcrossing rates (t(m)) ranged from 0.836 in the south to 0.897 in the north. Therefore, selfing is low but statistically significant. The equilibrium inbreeding coefficient (F(e)) calculated from t(m) was in good agreement with observed inbreeding coefficients, suggesting that weeping pi?on may be near equilibrium with respect to inbreeding and selection against selfed trees. Weeping pi?on was variable at all loci polymorphic in maxipi?on (Pinus maximartinezii) and, therefore, qualifies as a possible progenitor of maxipi?on. Because of the high level of diversity, reasonable levels of gene flow within the northern fragment of weeping pi?on's range, high rates of outcrossing, and, perhaps, only weak selection against inbred trees, protection in reserves would be a viable option for conservation.     

Understanding the relationship between the inbreeding coefficient and multilocus heterozygosity: theoretical expectations and empirical data

Geneticists have been interested in inbreeding and inbreeding depression since the time of Darwin. Two alternative approaches that can be used to measure how inbred an individual is involve the use of pedigree records to estimate inbreeding coefficients or molecular markers to measure multilocus heterozygosity. However, the relationship between inbreeding coefficient and heterozygosity has only rarely been investigated. In this paper, a framework to predict the relationship between the two variables is presented. In addition, microsatellite genotypes at 138 loci spanning all 26 autosomes of the sheep genome were used to investigate the relationship between inbreeding coefficient and multilocus heterozygosity. Multilocus heterozygosity was only weakly correlated with inbreeding coefficient, and heterozygosity was not positively correlated between markers more often than expected by chance. Inbreeding coefficient, but not multilocus heterozygosity, detected evidence of inbreeding depression for morphological traits. The relevance of these findings to the causes of heterozygosity--fitness correlations is discussed and predictions for other wild and captive populations are presented.     

Rhh: an R extension for estimating multilocus heterozygosity and heterozygosity-heterozygosity correlation

Individual multilocus heterozygosity estimates based on a limited number of loci are expected to correlate only weakly with the inbreeding level of an individual. Before using multilocus heterozygosity estimates in studies of inbreeding, their ability to capture information on inbreeding in the given setting should be tested. A convenient method for this is to compute the heterozygosity-heterozygosity correlation, i.e. the mean correlation between multilocus heterozygosity estimates calculated from random samples of loci, which should be positive if multilocus heterozygosity carries a signature of inbreeding. Rhh is an extension package for the statistical software r that estimates this correlation and calculates three measures of individual multilocus heterozygosity: homozygosity by loci, internal relatedness and standardized heterozygosity. The extension package is available through the CRAN (http://cran.r-project.org) and has a homepage at http://www.helsinki.fi/biosci/egru/research/software.