Interactive effects of environmental stress and inbreeding on reproductive traits in a wild bird population

1. Conservation biologists are concerned about the interactive effects of environmental stress and inbreeding because such interactions could affect the dynamics and extinction risk of small and isolated populations, but few studies have tested for these interactions in nature. 2. We used data from the long-term population study of song sparrows Melospiza melodia on Mandarte Island to examine the joint effects of inbreeding and environmental stress on four fitness traits that are known to be affected by the inbreeding level of adult birds: hatching success, laying date, male mating success and fledgling survival. 3. We found that inbreeding depression interacted with environmental stress to reduce hatching success in the nests of inbred females during periods of rain. 4. For laying date, we found equivocal support for an interaction between parental inbreeding and environmental stress. In this case, however, inbred females experienced less inbreeding depression in more stressful, cooler years. 5. For two other traits, we found no evidence that the strength of inbreeding depression varied with environmental stress. First, mated males fathered fewer nests per season if inbred or if the ratio of males to females in the population was high, but inbreeding depression did not depend on sex ratio. Second, fledglings survived poorly during rainy periods and if their father was inbred, but the effects of paternal inbreeding and rain did not interact. 6. Thus, even for a single species, interactions between the inbreeding level and environmental stress may not occur in all traits affected by inbreeding depression, and interactions that do occur will not always act synergistically to further decrease fitness.     

Indirect genetic effects and inbreeding: consequences of BLUP selection for socially affected traits on rate of inbreeding

Background

Social interactions often occur among living organisms, including aquatic animals. There is empirical evidence showing that social interactions may genetically affect phenotypes of individuals and their group mates. In this context, the heritable effect of an individual on the phenotype of another individual is known as an Indirect Genetic Effect (IGE). Selection for socially affected traits may increase response to artificial selection, but also affect rate of inbreeding.

Methods

A simulation study was conducted to examine the effect of Best Linear Unbiased Prediction (BLUP) selection for socially affected traits on the rate of inbreeding. A base scenario without IGE and three alternative scenarios with different magnitudes of IGE were simulated. In each generation, 25 sires and 50 dams were mated, producing eight progeny per dam. The population was selected for 20 generations using BLUP. Individuals were randomly assigned to groups of eight members in each generation, with two families per group, each contributing four individuals. “Heritabilities” (for both direct and indirect genetic effects) were equal to 0.1, 0.3 or 0.5, and direct–indirect genetic correlations were −0.8, −0.4, 0, 0.4, or 0.8. The rate of inbreeding was calculated from generation 10 to 20.

Results

For the base scenario, the rates of inbreeding were 4.09, 2.80 and 1.95% for “heritabilities” of 0.1, 0.3 and 0.5, respectively. Overall, rates of inbreeding for the three scenarios with IGE ranged from 2.21 to 5.76% and were greater than for the base scenarios. The results show that social interaction within groups of two families increases the resemblance between estimated breeding values of relatives, which, in turn, increases the rate of inbreeding.

Conclusion

BLUP selection for socially affected traits increased the rate of inbreeding. To maintain inbreeding at an acceptable rate, a selection algorithm that restricts the increase in mean kinship, such as optimum contribution selection, is required.     

Optimizing sampling design and sequencing strategy for the genomic analysis of quantitative traits in natural populations

Mapping the genes underlying ecologically relevant traits in natural populations is fundamental to develop a molecular understanding of species adaptation. Current sequencing technologies enable the characterization of a species’ genetic diversity across the landscape or even over its whole range. The relevant capture of the genetic diversity across the landscape is critical for a successful genetic mapping of traits and there are no clear guidelines on how to achieve an optimal sampling and which sequencing strategy to implement. Here we determine, through simulation, the sampling scheme that maximizes the power to map the genetic basis of a complex trait in an outbreeding species across an idealized landscape and draw genomic predictions for the trait, comparing individual and pool sequencing strategies. Our results show that quantitative trait locus detection power and prediction accuracy are higher when more populations over the landscape are sampled and this is more cost-effectively done with pool sequencing than with individual sequencing. Additionally, we recommend sampling populations from areas of high genetic diversity. As progress in sequencing enables the integration of trait-based functional ecology into landscape genomics studies, these findings will guide study designs allowing direct measures of genetic effects in natural populations across the environment.     

Density effects on life-history traits of an island lizard population

Long-term monitoring of life-history traits and the effects of density upon them were studied in an island population of the lizardEumeces okadae. Although life-history traits such as clutch size, egg size and the proportion of mature reproductive females varied little over 7 years in the intact population, manipulation of density to simulate decreased population density enhanced juvenile growth rate, age at first reproduction, frequency of female reproduction and size-specific clutch mass. In particular, the proportion of mature females reproducing annually increased almost 10 times from 5.6% to 53.8% after the removal of some lizards. However, body size at first reproduction and egg size were almost identical under both high and low density conditions. This study suggests that there were strong density-dependent effects on several life-history traits and thatE. okadae attained a density close to the carrying capacity of the environment.     

A comparative analysis of population structure of a forest tree, Eucalyptus globulus (Myrtaceae), using microsatellite markers and quantitative traits

Eucalyptus globulus (Myrtaceae) is a forest tree native to southeastern Australia, but is grown globally for pulpwood and timber. Eight microsatellite loci were used to determine the degree of selectively neutral differentiation between native populations of the geographic races of E. globulus that are used in a national breeding programme. Spatial differentiation was detected among 340 samples from across the species range (F _ST=0.09±0.02). Analysis of molecular variance showed that there was significant variation between the races, and an unweighted pair group method with arithmetic mean analysis of Nei’s genetic distance between races showed that geographically proximal races tended to be more closely related than geographically distant races. This contrasted markedly with analyses based on quantitative genetic data, where some races appeared to be highly divergent from their geographically closest neighbours. Comparison of racial differentiation based on quantitative (Q _ST) and molecular (F _ST) data suggested that at least five of the quantitative traits used for defining races of E. globulus have been influenced by natural selection, resulting in cases of both phenotypic divergence of parapatric races and phenotypic convergence of allopatric races. We conclude that selectively neutral molecular markers are more useful than quantitative genetic data for identifying the evolutionary affinities and lineages within E. globulus. However, both sources of information should be used in defining evolutionarily important units for conservation. The population structure observed in E. globulus has important consequences for future association studies and may also affect breeding strategies if significant genome co-adaptation has occurred.     

Estimation and testing of parent-of-origin effects for quantitative traits

Recent progress in developing family-based association methods has extended their use to the analysis of quantitative traits in the offspring and to the estimation, for dichotomous traits, of the relative contribution of genetic and environmental mechanisms for parent-of-origin effects. However, many traits of interest are not naturally measured on a binary scale yet are suspected or known to be influenced by imprinted genes, and there is consequent interest in seeking evidence for parent-of-origin effects at these loci. Here we show how simple linear models can be used to estimate these parent-of-origin effects for a broad class of phenotypes; in particular, normally distributed quantitative traits are easily dealt with.     

Estimating effects of a single gene and polygenes on quantitative traits from a diallel design

A genetic model is developed with additive and dominance effects of a single gene and polygenes as well as general and specific reciprocal effects for the progeny from a diallel mating design. The methods of ANOVA, minimum norm quadratic unbiased estimation (MINQUE), restricted maximum likelihood estimation (REML), and maximum likelihood estimation (ML) are suggested for estimating variance components, and the methods of generalized least squares (GLS) and ordinary least squares (OLS) for fixed effects, while best linear unbiased prediction, linear unbiased prediction (LUP), and adjusted unbiased prediction are suggested for analyzing random effects. Monte Carlo simulations were conducted to evaluate the unbiasedness and efficiency of statistical methods involving two diallel designs with commonly used sample sizes, 6 and 8 parents, with no and missing crosses, respectively. Simulation results show that GLS and OLS are almost equally efficient for estimation of fixed effects, while MINQUE (1) and REML are better estimators of the variance components and LUP is most practical method for prediction of random effects. Data from a Drosophila melanogaster experiment (Gilbert 1985a, Theor appl Genet 69:625–629) were used as a working example to demonstrate the statistical analysis. The new methodology is also applicable to screening candidate gene(s) and to other mating designs with multiple parents, such as nested (NC Design I) and factorial (NC Design II) designs. Moreover, this methodology can serve as a guide to develop new methods for detecting indiscernible major genes and mapping quantitative trait loci based on mixture distribution theory. The computer program for the methods suggested in this article is freely available from the authors.     

Methods for the study of cytoplasmic effects on quantitative traits

Summary The methods used to study cytoplasmic effects in quantitative traits often do not measure quantitative genetic parameters, while those that do are either complicated or do not take into account situations where the expression of cytoplasmic effects does not persist, but decreases in advanced generations. We present two simple models that take cytoplasmic effects and the quantitative genetic parameters into account. One of the models (A) is for cases where cytoplasmic effects remain constant through successive generations, and the second model (B) is for traits where cytoplasm-genotype interactions are present. This model also takes into account the decreasing persistence of cytoplasmic effects with advancing generations, which is often reported in the literature.     

Distinguishing adaptive from nonadaptive genetic differentiation: comparison of Q(ST) and F(ST) at two spatial scales

Genetic differentiation in 20 hierarchically sampled populations of wild barley was analyzed with quantitative traits, allozymes and Random Amplified Polymorphic DNAs (RAPDs), and compared for three marker types at two hierarchical levels. Regional subdivision for both molecular markers was much lower than for quantitative traits. For both allozymes and RAPDs, most loci exhibited minor or no regional differentiation, and the relatively high overall estimates of the latter were due to several loci with exceptionally high regional differentiation. The allozyme- and RAPD-specific patterns of differentiation were concordant in general with one another, but not with quantitative trait differentiation. Divergent selection on quantitative traits inferred from very high regional Q(ST) was in full agreement with our previous results obtained from a test of local adaptation and multilevel selection analysis. In contrast, most variation in allozyme and RAPD variation was neutral, although several allozyme loci and RAPD markers were exceptional in their levels of regional differentiation. However, it is not possible to answer the question whether these exceptional loci are directly involved in the response to selection pressure or merely linked to the selected loci. The fact that Q(ST) and F(ST) did not differ at the population scale, that is, within regions, but differed at the regional scale, for which local adaptation has been previously shown, implies that comparison of the level of subdivision in quantitative traits, as compared with molecular markers, is indicative of adaptive population differentiation only when sampling is carried out at the appropriate scale.     

Joint effects of self-fertilization and population structure on mutation load, inbreeding depression and heterosis

Both the spatial distribution of organisms and their mode of reproduction have important effects on the change in allele frequencies within populations. In this article, we study the combined effect of population structure and the rate of partial selfing of organisms on the efficiency of selection against recurrent deleterious mutations. Assuming an island model of population structure and weak selection, we express the mutation load, the within- and between-deme inbreeding depression, and heterosis as functions of the frequency of deleterious mutants in the metapopulation; we then use a diffusion model to calculate an expression for the equilibrium probability distribution of this frequency of deleterious mutants. This allows us to derive approximations for the average mutant frequency, mutation load, inbreeding depression, and heterosis, the simplest ones being Equations 35-39 in the text. We find that population structure can help to purge recessive deleterious mutations and reduce the load for some parameter values (in particular when the dominance coefficient of these mutations is <0.2-0.3), but that this effect is reversed when the selfing rate is above a given value. Conversely, within-deme inbreeding depression always decreases, while heterosis always increases, with the degree of population subdivision, for all selfing rates.     

Sorting out the effects of Wolbachia, genotype and inbreeding on life-history traits of a spider mite

Wolbachia bacteria manipulate host reproduction by inducing cytoplasmic incompatibility (CI) and sex ratio distortion. Wolbachia are transmitted from mother to offspring through the cytoplasm of the egg. Therefore, reproduction of Wolbachia is tightly coupled to reproduction of its host. Mathematical analysis predicts that in the course of evolution, traits that reduce the physiological costs of the infection will be selectively favored. For a Wolbachia-host system to evolve, traits under selection must have some genetic component and variation must be present in the population. We have previously established that highly inbred isofemale lines of the two-spotted spider mite Tetranychus urticae may differ regarding the effects of infection by Wolbachia, and that at least some of the traits affected had a genetic component. However, the effects measured could have been affected by the fact that the lines were severely inbred prior to the experiments. In this paper we attempt to distinguish between the effects of Wolbachia, isofemale line, and inbreeding. We show that Wolbachia did not affect longevity but infected females produced smaller clutch sizes, more daughter-biased sex ratios and had decreased F1 mortality; between-line variation was found for clutch size, F1 mortality and sex ratio; finally, inbreeding resulted in an overall reduction of clutch sizes, and a change in survival curves and mean longevity. This revised version was published online in July 2006 with corrections to the Cover Date.     

Adaptive genetic differentiation in a predominantly self-pollinating species analyzed by transplanting into natural environment, crossbreeding and Q(ST)-F(ST) test

? Both genetic drift and natural selection result in genetic/phenotypic differentiation over space. I analyzed the role of local adaptation in the genetic differentiation of populations of the annual grass Hordeum spontaneum sampled along an aridity gradient. ? The study included the introduction of plants having desert vs nondesert origin into natural (desert) environment, analysis of population differentiation in allozymes and random amplified polymorphic DNA (RAPD) markers vs phenotypic traits (Q(ST) -F(ST) comparison), and planting interpopulation hybrids under simulated desert conditions in a glasshouse. ? The results of the home advantage test, Q(ST) -F(ST) comparison and crossbreeding were consistent with local adaptation; that is, that differentiation of the desert plants from plants of nondesert origin in phenotypic traits was adaptive, giving them home advantage. Each method used provided additional, otherwise unavailable, information, meaning that they should be viewed as complementary rather than alternative approaches. ? Gene flow from adjacent populations (i.e. populations experiencing the desert environment) via seeds (but not pollen) had a positive effect on fitness by enhancing natural selection and counteracting drift. At the same time, the effect of genes from the species distributional core (nondesert plants) by either seed or pollen had a negative fitness effect despite its enriching effect on neutral diversity. The pattern of outbreeding depression observed in interpopulation hybrids (F(1) ) and their segregating progeny (F(2) ) was inconsistent with underdominance, but indicated the presence of additive, dominance and epistatic effects.     

A test of quantitative genetic theory using Drosophila- effects of inbreeding and rate of inbreeding on heritabilities and variance components

Inbreeding is expected to decrease the heritability within populations. However, results from empirical studies are inconclusive. In this study, we investigated the effects of three breeding treatments (fast and slow rate of inbreeding - inbred to the same absolute level - and a control) on heritability, phenotypic, genetic and environmental variances of sternopleural bristle number in Drosophila melanogaster. Heritability, and phenotypic, genetic and environmental variances were estimated in 10 replicate lines within each of the three treatments. Standard least squares regression models and Bayesian methods were used to analyse the data. Heritability and additive genetic variance within lines were higher in the control compared with both inbreeding treatments. Heritabilities and additive genetic variances within lines were higher in slow compared with fast inbred lines, indicating that slow inbred lines retain more evolutionary potential despite the same expected absolute level of inbreeding. The between line variance was larger with inbreeding and more than twice as large in the fast than in the slow inbred lines. The different pattern of redistribution of genetic variance within and between lines in the two inbred treatments cannot be explained invoking the standard model based on selective neutrality and additive gene action. Environmental variances were higher with inbreeding, and more so with fast inbreeding, indicating that inbreeding and the rate of inbreeding affect environmental sensitivity. The phenotypic variance decreased with inbreeding, but was not affected by the rate of inbreeding. No inbreeding depression for mean sternopleural bristle number was observed in this study. Considerable variance between lines in additive genetic variance within lines was observed, illustrating between line variation in evolutionary potential.     

The effect of self-fertilization, inbreeding depression, and population size on autopolyploid establishment

The minority cytotype exclusion principle describes how random mating between diploid and autotetraploid cytotypes hinders establishment of the rare cytotype. We present deterministic and stochastic models to ascertain how selfing, inbreeding depression, unreduced gamete production, and finite population size affect minority cytotype exclusion and the establishment of autotetraploids. Results demonstrate that higher selfing rates and lower inbreeding depression in autotetraploids facilitate establishment of autotetraploid populations. Stochastic effects due to finite population size increase the probability of polyploid establishment and decrease the mean time to tetraploid fixation. Our results extend the minority cytotype exclusion principle to include important features of plant reproduction and demonstrate that variation in mating system parameters significantly influences the conditions necessary for polyploid establishment.     

Dispersal,kinship and inbreeding in an island population of the Great Tit

We investigated the kinship structure of an island population of the Great Tit (Parus major). Kinship of birds could be inferred by comparing their family trees. Dispersal was also studied to explain the observed pattern of kinship. On the island of Vlieland the tits breed in several wooded areas. Both males and females preferred to breed in their natal area; males did so more strongly than females. Hence gene flow between the areas is restricted. However, within the largest wooded area females showed random dispersal, while males showed a slight tendency to breed near their natal site. The degree of kinship of neighbouring birds is a suitable control group for the relatedness of partners that takes into account the effects of dispersal. In the largest wooded area, birds were on average equally related to their partner and to their neighbours. Moreover, the mean coefficient of kinship between male and female neighbours was equal to the average kinship in this part of the population. We conclude that mating is random with respect to kinship. There is no evidence for avoidance of inbreeding. It is unlikely that kin recognition plays an important role in the process of mate choice in this population of Great Tits. We suggest that ecological factors are the main causes for the observed patterns of dispersal and mating. On the island more female than male immigrants enter the population each year. Incidental data indicate an exchange of birds between the population studied and surrounding populations. Ancestries of immigrants are not known, and indeed a first analysis of all birds, including immigrants, showed that males were more closely related than females. However, differential immigration could not fully explain the observed difference in kinship. The presence of local adaptation in males is suggested as a possible additional cause.     

Low effective population size and evidence for inbreeding in an overexploited flatfish, plaice (Pleuronectes platessa L.)

Overexploitation and subsequent collapses of major worldwide fisheries has made it clear that marine stocks are no inexhaustible. Unfortunately, the perception remains that marine fished are resilient to large population reductions, as even a commercially 'collapsed' stock will still consist of millions of individuals. Coupled with this notion is the idea that fisheries can, therefore, have little effect on the genetic diversity of stocks. We used DNA from archived otoliths collected between 1924 and 1972 together with 2002 juvenile;s tissue to estimate effective population size (Ne) in plaice (Pleuronrctes platessa). Ne was estimated at 20,000 in the North Sea and 2000 in Iceland. These values are five orders of magnitude smaller than the estimated census size foe the two locations. Populations examined between 1924 and 1960 were in Hardy-Weinberg equilibrium, whereas populations examined after 1970 were not. Extensive testing was performed to rule out genotyping artefacts and Wahlund effects. The significant heterozygote deficiencies found from 1970 onward were attributed to inbreeding. The emergence of inbreeding between 1905 and 19070 coincides with the increase in fishing mortality after World War II. Although the biological mechanisms remain speculative, our demonstration of inbreeding signals the need for understanding the social and mating behaviour in commercially important fishes.     

Analysis of founders and performance test effects on an autochthonous horse population through pedigree analysis: structure,genetic variability and inbreeding

The Maremmano is an autochthonous Italian horse breed, which probably descended from the native horses of the Etruscans (VI century B.C.); the Studbook was acknowledged in 1980, and it includes 12 368 horses born from that year up to 2015. The aim of this study was to evaluate the effect of the selection program on the genetic variability of the Maremmano population; the analysis was performed using both the ‘Endog v 4.8’ program available at http://webs.ucm.es/info/prodanim/html/JP_Web.htm and in-house software on official pedigree data. Four Reference Populations were considered, and the most important one was the population of the 12 368 Maremmano horses officially registered in the National Studbook. The pedigree completeness of this population was very good because it was more than 90% at the third parental generation and more than 70% at the fifth generation; the pedigree traced back to a maximum of 10.50 generations with an average of 3.30 complete generations and 5.70 equivalent complete generations. The average generation interval was 10.65±4.72 years, with stallions used for longer periods than mares. The intervals ranged from 10.15±4.45 (mother–daughter) to 10.99±4.93 (father–daughter). The effective number of founders (f_e) was 74 and the effective number of ancestors (f_a) was 30 so that the ratio f_e/f_a was 2.47. The founder genome equivalents (f_g) was 13.72 with a ratio f_g/f_e equal to 0.18. The mean of the genetic conservation index was 5.55±3.37, and it ranged from 0.81 to 21.32. The average inbreeding coefficient was 2.94%, with an increase of 0.1%/year, and the average relatedness coefficient was 5.52%. The effective population size (N_e) computed by an individual increase in inbreeding was 68.1±13.00; the N_e on equivalent generations was 42.00, and this value slightly increased to 42.20 when computed by Log regression on equivalent generations. The analysis confirmed the presence of seven traditional male lines. The percentage of Thoroughbred blood in the foals born in 2015 was 20.30% and has increased 0.21%/year since 1980; in particular, it increased more than twice (0.51%/year) until 1993 and afterwards slightly fluctuated. The pedigree analysis confirmed the completeness of genealogical information and the traditional importance that breeders gave to the male lines; although the genetic diversity of Maremmano seemed to be not endangered by the selection program, some effects on the population structure were found and a more scientific approach to genetic conservation should be incorporated in the selection plans.