Influence of Gene Flow and Breeding Tactics on Gene Diversity within Populations

Expressions describing the accumulation of gene correlations within and among lineages and individuals of a population are derived. The model permits different migration rates by males and females and accounts for various breeding tactics within lineages. The resultant equations enable calculation of the probabilistic quantities for the fixation indices, rates of loss of genetic variation, accumulation of inbreeding, and coefficients of relationship for the population at any generation. All fixation indices were found to attain asymptotic values rapidly despite the consistent loss of genetic variation and accumulation of inbreeding within the population. The time required to attain asymptotic values, however, was prolonged when gene flow among lineages was relatively low (less than 20%). The degree of genetic differentiation among breeding groups, inbreeding coefficients, and gene correlations within lineages were found to be primarily functions of breeding tactics within groups rather than gene flow among groups. Thus, the asymptotic value of S. Wright's island model is not appropriate for describing genetic differences among groups within populations. An alternative solution is provided that under limited conditions will reduce to the original island model. The evolution of polygynous breeding tactics appears to be more favorable for promoting intragroup gene correlations than modification of migration rates. Inbreeding and variance effective sizes are derived for populations that are structured by different migration and breeding tactics. Processes that reduce the inbreeding effective population size result in a concomitant increase in variance effective population size.     

Genetic Diversity of a Tropical Tree Species,Shorea leprosula Miq. (Dipterocarpaceae), in Malaysia: Implications for Conservation of Genetic Resources and Tree Improvement1

Genetic diversity and population genetic structure of Shorea leprosula was investigated using seven natural populations distributed throughout Peninsular Malaysia and one natural population from Borneo. The mean population and species level genetic diversity were exceptionally high (H_e= 0.369 ± 0.025 and 0.406 ± 0.070, respectively). Heterozygosity varied among populations, ranging from 0.326 to 0.400, with the highest values found in the populations from central Peninsular Malaysia. Correlations among ecological factors (longitude, latitude, and annual rainfall) were not significant (P > 0.05), indicating that these ecological variables were not responsible for the observed genetic differences among populations. The Bangi adult population exhibited a higher level of observed heterozygosity but lower fixation indices in comparison to its seedling population. All other seedling populations also showed positive fixation indices (f), indicating a general excess of homozygotes. This also may suggest selection against homozygotes between the seedling and adult stages. A low level of population differentiation was detected (G_ST= 0.117 with the Lambir population and GST= 0.085 without the Lambir population). Furthermore, gene flow (N_m) between populations was not significantly correlated with geographical distances for the populations within Peninsular Malaysia. Cluster analysis also did not reflect geographical proximity and gave little insight into the genetic relatedness of the populations. This may indicate that the populations sampled are part of a continuous population with fragmentation having occurred in the recent past.     

Bayesian analysis of genetic differentiation between populations

We introduce a Bayesian method for estimating hidden population substructure using multilocus molecular markers and geographical information provided by the sampling design. The joint posterior distribution of the substructure and allele frequencies of the respective populations is available in an analytical form when the number of populations is small, whereas an approximation based on a Markov chain Monte Carlo simulation approach can be obtained for a moderate or large number of populations. Using the joint posterior distribution, posteriors can also be derived for any evolutionary population parameters, such as the traditional fixation indices. A major advantage compared to most earlier methods is that the number of populations is treated here as an unknown parameter. What is traditionally considered as two genetically distinct populations, either recently founded or connected by considerable gene flow, is here considered as one panmictic population with a certain probability based on marker data and prior information. Analyses of previously published data on the Moroccan argan tree (Argania spinosa) and of simulated data sets suggest that our method is capable of estimating a population substructure, while not artificially enforcing a substructure when it does not exist. The software (BAPS) used for the computations is freely available from http://www.rni.helsinki.fi/~mjs.     

Genetic structure among populations in the endemic Hawaiian Plantago lineage: insights from microsatellite variation

Endemic Hawaiian species in the genus Plantago show considerable morphological and ecological diversity. Despite their variation, a recent phylogenetic analysis based on DNA sequence data showed that the group is monophyletic and that sequence variation among species and morphotypes is low. This lack of sequence polymorphisms resulted in an inability to resolve species and population affinities within the most recently derived clade of this lineage. To assess species boundaries, population genetic structure and interpopulation connectivity among the morphologically and ecologically distinct populations within this clade, genetic variation was examined using eight microsatellite loci. Within‐population genetic diversity was found to be lowest in the Maunaiu, Hawai'i population of the endangered P. hawaiensis, and highest in the large P. pachyphylla population from 'Eke, West Maui. Isolation by distance across the range of populations was detected and indicated restricted dispersal. This result is likely to be attributable to few interisland dispersal events in the evolutionary history of this lineage. Genetic differentiation within islands tended to be higher among populations occurring in contrasting bog and woodland habitats, suggesting ecological barriers to gene flow and the potential role of ecological divergence in population diversification. Overall, these results are consistent with findings from phylogenetic analysis of the entire lineage. Our data bring new insights regarding patterns of dispersal and population genetic structure to this endemic and endangered group of island taxa. As island environments become increasingly fragmented, information of this type has important implications for the successful management of these fragile populations and habitats.     

Habitat geometry, population viscosity and the rate of genetic drift

In all natural populations, individuals located close to one another tend to interact more than those further apart. The extent of population viscosity can have important implications for ecological and evolutionary processes. Here we develop a spatially explicit population model to examine how the rate of genetic drift depends upon both spatial population structure and habitat geometry. The results show that the time to fixation for a new and selectively neutral mutation is dramatically increased in viscous populations. Furthermore, in viscous populations the time to fixation depends critically on habitat geometry. Fixation time for populations of identical size increases markedly as landscape width decreases and length increases. We suggest that similar effects will also be important in metapopulations, with the spatial arrangement of subpopulations and their connectivity likely to determine the rate of drift. We argue that the recent increases in computer power should facilitate major advances in our understanding of evolutionary landscape ecology over the next few years, and suggest that the time is ripe for a unification of spatial population dynamics theory, landscape ecology and population genetics.     

Age structure: an indicator to monitor populations of large herbivores

Most populations of large herbivores are hunted over much of their distribution. The size of these populations over time must be monitored to detect responses to management actions and environmental changes. Here, I first review the limits in precision and accuracy of most methods used to estimate population size in wild herbivores. I then show how density-dependence of age structure can be used to monitor harvested populations, using simulated and empirical data. Accurate estimates of population size are difficult, costly in time and effort, and mostly impracticable at a large scale. Indices are widely used to monitor population size at a large scale, but are rarely tested, and assumptions of constant detection rates over years are questionable. In a recent approach, the effects of population density on habitat quality and individual performance were used as ecological indicators to infer changes in population size over time. Populations of large mammals are strongly age structured. In harvested populations, density-dependence and harvest affect age structure through changes in age-specific vital rates. Because age- and sex-specific vital rates are density-dependent, changes over time in age structure should reflect the response of population to changes in environmental conditions such as density and harvest rate. Time series of age-structure may be used as an ecological indicator to monitor harvested populations, and to infer their demographic trajectory, especially when used within a set of ecological indicators. However, attention needs to be paid in the use of ecological indicators because of the potential bias introduced in data collection by hunter selectivity.     

Toward ecologically scaled landscape indices

Nature conservation is increasingly based on a landscape approach rather than a species approach. Landscape planning that includes nature conservation goals requires integrated ecological tools. However, species differ widely in their response to landscape change. We propose a framework of ecologically scaled landscape indices that takes into account this variation. Our approach is based on a combination of field studies of spatially structured populations (metapopulations) and model simulations in artificial landscapes. From these, we seek generalities in the relationship among species features, landscape indices, and metapopulation viability. The concept of ecological species profiles is used to group species according to characteristics that are important in metapopulations' response to landscape change: individual area requirements as the dominant characteristic of extinction risk in landscape patches and dispersal distance as the main determinant of the ability to colonize patches. The ecological profiles and landscape indices are then integrated into two ecologically scaled landscape indices (ESLI): average patch carrying capacity and average patch connectivity. The field data show that the fraction of occupied habitat patches is correlated with the two ESLI. To put the ESLI into a perspective of metapopulation persistence, we determine the viability for six ecological profiles at different degrees of habitat fragmentation using a metapopulation model and computer-generated landscapes. The model results show that the fraction of occupied patches is a good indicator for metapopulation viability. We discuss how ecological profiles, ESLI, and the viability threshold can be applied for landscape planning and design in nature conservation.     

Census data aggregation decisions can affect population‐level inference in heterogeneous populations

1. Conservation and population management decisions often rely on population models parameterized using census data. However, the sampling regime, precision, sample size, and methods used to collect census data are usually heterogeneous in time and space. Decisions about how to derive population‐wide estimates from this patchwork of data are complicated and may bias estimated population dynamics, with important implications for subsequent management decisions.
2. Here, we explore the impact of site selection and data aggregation decisions on pup survival estimates, and downstream estimates derived from parameterized matrix population models (MPMs), using a long‐term dataset on grey seal (Halichoerus grypus) pup survival from southwestern Wales. The spatiotemporal and methodological heterogeneity of the data are fairly typical for ecological census data and it is, therefore, a good model to address this topic.
3. Data were collected from 46 sampling locations (sites) over 25 years, and we explore the impact of data handling decisions by varying how years and sampling locations are combined to parameterize pup survival in population‐level MPMs. We focus on pup survival because abundant high‐quality data are available on this developmental stage.
4. We found that survival probability was highly variable with most variation being at the site level, and poorly correlated among sampling sites. This variation could generate marked differences in predicted population dynamics depending on sampling strategy. The sample size required for a confident survival estimate also varied markedly geographically.
5. We conclude that for populations with highly variable vital rates among sub‐populations, site selection and data aggregation methods are important. In particular, including peripheral or less frequently used areas can introduce substantial variation into population estimates. This is likely to be context‐dependent, but these choices, including the use of appropriate weights when summarizing across sampling areas, should be explored to ensure that management actions are successful.

     

Root architecture,plant size and soil nutrient variation in natural populations of <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Phenotypic variation in ecologically important traits may vary at large and small geographic scales, and may be shaped by natural selection. Here our explicit aim is to evaluate phenotypic differentiation among local populations and examine its relationship with ecological edaphic and climatic features that could lead to local adaptation. We characterized six populations of the model plant Arabidopsis thaliana over 3 years in the field in its native range. At each site, we measured edaphic conditions and aboveground and belowground phenotypes. In addition, we grew plants from the six characterized populations in a common greenhouse along with an additional fifteen populations from the Iberian Peninsula to examine evolutionary and ecological differentiation among populations, and relationships between geographic and ecological distance to phenotypic differences among populations. Significant differences in aboveground and belowground traits, population density, and micro- and macro-nutrient soil concentrations were found among the field populations. In particular, root architectural traits differed significantly among field populations. Complex patterns of ecological differences among population and plant phenotypes emerged when examining edaphic conditions in the Extremadura region, and geographic and climate variables at a broader scale of the Iberian Peninsula. We report levels of phenotypic variation at the local scale comparable to those found at broad geographic scales and report that local edaphic conditions contribute to population-level phenotypic variation in root and shoot traits. To our knowledge, these are the first reports of among population root architectural variation from natural field populations for this model organism. We demonstrate how ecological features, such as soil nutrients, can be associated with the phenotypic variation of A. thaliana measured in natural populations and may contribute to adaptive differentiation at a local scale.     

Heterozygosity-fitness correlations among wild populations of European tree frog (Hyla arborea) detect fixation load

Quantifying the impacts of inbreeding and genetic drift on fitness traits in fragmented populations is becoming a major goal in conservation biology. Such impacts occur at different levels and involve different sets of loci. Genetic drift randomly fixes slightly deleterious alleles leading to different fixation load among populations. By contrast, inbreeding depression arises from highly deleterious alleles in segregation within a population and creates variation among individuals. A popular approach is to measure correlations between molecular variation and phenotypic performances. This approach has been mainly used at the individual level to detect inbreeding depression within populations and sometimes at the population level but without consideration about the genetic processes measured. For the first time, we used in this study a molecular approach considering both the interpopulation and intrapopulation level to discriminate the relative importance of inbreeding depression vs. fixation load in isolated and non-fragmented populations of European tree frog (Hyla arborea), complemented with interpopulational crosses. We demonstrated that the positive correlations observed between genetic heterozygosity and larval performances on merged data were mainly caused by co-variations in genetic diversity and fixation load among populations rather than by inbreeding depression and segregating deleterious alleles within populations. Such a method is highly relevant in a conservation perspective because, depending on how populations lose fitness (inbreeding vs. fixation load), specific management actions may be designed to improve the persistence of populations.     

Subpopulation genetic structure of a plant panmictic population of Castanea sequinii as revealed by microsatellite markers

Castanea squinii Dode,an endemic tree widely distributed in China,plays an important role both in chestnut breeding and forest ecosystem function.The spatial genetic structure within and among populations is an important part of the evolutionary and ecological genetic dynamics of natural populations,and can provide insights into effective conservation of genetic resources.In the present study,the spatial genetic structure of a panmictic natural population of C.sequinii in the Dabie Mountain region was investigated using microsatellite markers.Nine prescreened microsatellite loci generated 29-33 alleles each,and were used for spatial autocorrelation analysis.Based on Moran's I coefficient,a panmictic population of C.sequinii in the Dabie Mountain region was found to be lacking a spatial genetic structure.These results suggest that a high pollen-mediated gene flow among subpopulations counteract genetic drift and/or genetic differentiation and plays an important role in maintaining a random and panmictic population structure in C.sequinii populations.Further,a spatial genetic structure was detected in each subpopulation's scale (0.228 km),with all three subpopulations showing significant fine-scale structure.The genetic variation was found to be nonrandomly distributed within 61 m in each subpopulation (Moran's I positive values).Although Moran's I values varied among the different subpopulations,Moran's I in all the three subpopulations reached the expected values with an increase in distances,suggesting a generally patchy distribution in the subpopulations.The fine-scale structure seems to reflect restricted seed dispersal and microenvironment selection in C.sequinii.These results have important implications for understanding the evolutionary history and ecological process of the natural population of C.sequinii and provide baseline data for formulating a conservation strategy of Castanea species.     

Variation among populations of Galapagos land iguanas (Conolophus): contrasts of phylogeny and ecology

A phylogenetic scheme derived via multivariate analyses of adaptively neutral scale characteristics is compared to patterns of ecological adaptation in body size and shape, hatchling size, clutch size, and reproductive seasonality, in extant populations of Galapagos land iguanas (genus Conolophus). Three groups of land iguana populations are identified, the oldest being the population of Isla Santa Fe, followed by the populations of the central islands (Santa Cruz, Plaza Sur and Baltra), the youngest populations are those of the western islands (Fernandina and Isabela). Patterns of ecological similarity among these populations are not concordant with phylogenetic lineage. Populations most similar in ecological characteristics are often phylogenetically divergent. Adaptation to local conditions by iguana populations is apparently more important than phylogenetic constraint in explaining variation in ecological characteristics. The assumption that phylogenetically closely-related organisms are also ecologically more similar than less closely-related organisms is not supported by this evidence. Some previous studies may have been misled by using ecological characteristics to derive phylogenetic lineages, resulting in circular support of the assumption.     

Genetic relationships of American alligator populations distributed across different ecological and geographic scales

Although much work has been conducted on coastal populations of the American alligator (Alligator mississippiensis), less is known about the population dynamics and genetic structure of populations of alligators confined to inland habitats. DNA microsatellite loci, derived from the American alligator, were used to investigate patterns of genetic variation within and between populations of alligators distributed at coastal and inland localities in Texas. These data were used to evaluate the genetic discreteness of different alligator stocks relative to their basic ecology at these sites. Observed mean heterozygosities across seven loci for both coastal and inland populations ranged from 0.50-0.61, with both inland and coastal populations revealing similar patterns of variation. Measures of F(st) revealed significant population differentiation among all populations; however, analyses of molecular variance (AMOVAs) failed to demonstrate any apparent geographic pattern relative to the population differentiation indicated by F(st) values. Each population contained unique alleles for at least one locus. Additionally, assignment tests based on the distribution of genotypes placed 76% of individuals to their source population. These genetic data suggest considerable subdivision among alligator populations, possibly influenced by demographic and life history differences as well as barriers to dispersal. These results have clear implications for management. Rather than managing alligators in Texas as a single panmictic population, translocation programs and harvest quotas should consider the ecological and genetic distinctiveness of local alligator populations.     

Genetic structure of Lake Magadi tilapia populations

The control region of the mitochondrial DNA haplotype frequencies were significantly different among the two separate lagoon populations of Oreochromis alcalicus grahami in Lake Magadi and of O. a. alcalicus from lake Natron, and DNA fingerprint similarity indices were significantly higher for intra-population comparisons of the two Magadi lagoon populations and the Lake Natron population than the inter-population similarity indices among these populations. A modified F _st measure indicated population sub-division and the phylogeographic partitioning of the VNTR fragments observed were unique to specific populations further indicating substantial genetic differentiation. The lagoon populations within Lake Magadi demonstrated the same degree of genetic differentiation as either of these populations did to the outgroup (the Lake Natron population). There appears to be limited gene flow between Lake Magadi tilapia populations and this population structure has important implications for protecting locally adapted populations within this unique ecosystem.     

Estimation and adjustment of microsatellite null alleles in nonequilibrium populations

Nonamplified (null) alleles are a common feature of microsatellite genotyping and can bias estimates of allele and genotype frequencies, thereby hindering population genetic analyses. The frequency of microsatellite null alleles in diploid populations can be estimated for populations that are in Hardy–Weinberg equilibrium. However, many microsatellite data sets are from nonequilibrium populations, often with known inbreeding coefficients (F) or fixation indices (F_IS or F_ST). Here, we propose a novel null allele estimator that can be used to estimate the null allele frequency and adjust visible allele frequencies in populations for which independent estimates of F, F_IS or F_ST are available. The algorithm is currently available as an Excel macro that can be downloaded at no cost from http://www.microchecker.hull.ac.uk/ and will be incorporated into the software micro ‐checker .     

Genetic diversity of the endangered and narrow endemic Piperia yadonii (Orchidaceae) assessed with ISSR polymorphisms

Highly endangered plants that are also narrow endemics are generally found to be genetically depauperate and thus are exceedingly susceptible to ecological and anthropological threats that can lead to their extinction. Piperia yadonii is restricted to a single California county within a biodiversity hotspot. We used nine primers to generate intersimple sequence repeat (ISSR) data to assess its genetic diversity and structure. Within each population, 99% of the loci were polymorphic, expected heterozygosity was low, and a majority of the loci were shared with few other populations. Forty percent of the total variation could be attributed to population differentiation while the rest (60%) resides within populations, and the genetic distances between populations were independent of the corresponding geographical distances. High divergence among populations is likely due to fragmentation and limited gene flow. Each population contains several private loci, and ideally, each should be protected to preserve the overall diversity of the species. Because P. yadonii currently retains a modest amount of genetic variation among individuals within populations, preserving and expanding the habitat at each site to allow natural expansion of populations would be additional strategies for its conservation before populations become too small to persist naturally.     

Estimating Regions of Oceanographic Importance for Seabirds Using A-Spatial Data

Advances in GPS tracking technologies have allowed for rapid assessment of important oceanographic regions for seabirds. This allows us to understand seabird distributions, and the characteristics which determine the success of populations. In many cases, quality GPS tracking data may not be available; however, long term population monitoring data may exist. In this study, a method to infer important oceanographic regions for seabirds will be presented using breeding sooty shearwaters as a case study. This method combines a popular machine learning algorithm (generalized boosted regression modeling), geographic information systems, long-term ecological data and open access oceanographic datasets. Time series of chick size and harvest index data derived from a long term dataset of Maori ‘muttonbirder’ diaries were obtained and used as response variables in a gridded spatial model. It was found that areas of the sub-Antarctic water region best capture the variation in the chick size data. Oceanographic features including wind speed and charnock (a derived variable representing ocean surface roughness) came out as top predictor variables in these models. Previously collected GPS data demonstrates that these regions are used as “flyways” by sooty shearwaters during the breeding season. It is therefore likely that wind speeds in these flyways affect the ability of sooty shearwaters to provision for their chicks due to changes in flight dynamics. This approach was designed to utilize machine learning methodology but can also be implemented with other statistical algorithms. Furthermore, these methods can be applied to any long term time series of population data to identify important regions for a species of interest.     

Modelling the effects of genetics and habitat on the demography of a grassland herb

There is growing evidence that genetic and ecological factors interact in determining population persistence. The demographic effects of inbreeding depression can largely depend on the ecological milieu. We used demographic data of the perennial herb Succisa pratensis from six populations in grazed and ungrazed sites with different soil moisture. We built an individual-based model assessing the demographic consequences of inbreeding depression in populations with different management and habitat. Today this plant has to cope with severe landscape fragmentation, deteriorating habitat conditions in terms of decreasing grazing intensity, and the effects of inbreeding depression. For each population we performed simulations testing two inbreeding depression hypotheses (partial dominance and overdominance) and three epistatic functions among loci. The results indicated stronger inbreeding depression effects for populations in unfavourable sites without grazing or in xeric habitats compared to populations in favourable mesic sites with grazing. Overall, we found stronger effects with overdominance, a result that emphasizes the importance of understanding the genetic mechanisms of inbreeding depression. Hence, management practices can interact with the genetic consequences of inbreeding depression in population dynamics, which may have important implications for plant population ecology and evolutionary dynamics of inbreeding depression.     

Measuring the contribution of evolution to community trait structure in freshwater zooplankton

There are currently few predictions about when evolutionary processes are likely to play an important role in structuring community features. Determining predictors that indicate when evolution is expected to impact ecological processes in natural landscapes can help researchers identify eco-evolutionary ‘hotspots', where eco-evolutionary interactions are more likely to occur. Using data collected from a survey in freshwater cladoceran communities, landscape population genetic data and phenotypic trait data measured in a common garden, we applied a Bayesian linear model to assess whether the impact of local trait evolution in the keystone species Daphnia magna on cladoceran community trait values could be predicted by population genetic properties (within-population genetic diversity, genetic distance among populations), ecological properties (Simpson's diversity, phenotypic divergence) or environmental divergence. We found that the impact of local trait evolution varied among communities. Moreover, community diversity and phenotypic divergence were found to be better predictors of the contribution of evolution to community trait values than environmental features or genetic properties of the evolving species. Our results thus indicate the importance of ecological context for the impact of evolution on community features. Our study also demonstrates one way to detect signatures of eco-evolutionary interactions in communities inhabiting heterogeneous landscapes using survey data of contemporary ecological and evolutionary structure.     

Effect of Habitat Type and Population Density on the Stress Level of Midday Gerbils (Meriones meridianus Pall.) in Free-Living Populations

New methods of noninvasive evaluation of the endocrine status of animals by the content of hormones in their feces were used to study the relationship between the stress level and social (population density) and ecological (habitat integrity) indices in natural populations of midday gerbil in the Southern Kalmykia in 2000. Stress level proved to increase with habitat disturbance but did not depend on population density of animals. The obtained data are discussed in the context of species-specific ecology and social behavior of midday gerbils.