Evolution of reduced post-copulatory molecular interactions in Drosophila populations lacking sperm competition.

In many species with internal fertilization, molecules transferred in the male ejaculate trigger and interact with physiological changes in females. It is controversial to what extent these interactions between the sexes act synergistically to mediate the female switch to a reproductive state or instead reflect sexual antagonism evolved as a by product of sexual selection on males. To address this question, we eliminated sexual selection by enforcing monogamy in populations of Drosophila melanogaster for 65 generations and then measured the expression of male seminal fluid protein genes and genes involved in the female response to mating. In the absence of sperm competition, male and female reproductive interests are perfectly aligned and any antagonism should be reduced by natural selection. Consistent with this idea, males from monogamous populations showed reduced expression of seminal fluid protein genes, 16% less on average than in polygamous males. Further, we identified 428 genes that responded to mating in females. After mating, females with an evolutionary history of monogamy exhibited lower relative expression of genes that were up regulated in response to mating and higher expression of genes that were down-regulated--in other words, their post-mating transcriptome appeared more virgin-like. Surprisingly, these genes showed a similar pattern even before mating, suggesting that monogamous females evolved to be less poised for mating and the accompanying receipt of male seminal fluid proteins. This reduced investment by both monogamous males and females in molecules involved in post-copulatory interactions points to a pervasive role of sexual conflict in shaping these interactions.

Estimating sampling error of evolutionary statistics based on genetic covariance matrices using maximum likelihood.

We explore the estimation of uncertainty in evolutionary parameters using a recently devised approach for resampling entire additive genetic variance-covariance matrices (G). Large-sample theory shows that maximum-likelihood estimates (including restricted maximum likelihood, REML) asymptotically have a multivariate normal distribution, with covariance matrix derived from the inverse of the information matrix, and mean equal to the estimated G. This suggests that sampling estimates of G from this distribution can be used to assess the variability of estimates of G, and of functions of G. We refer to this as the REML-MVN method. This has been implemented in the mixed-model program WOMBAT. Estimates of sampling variances from REML-MVN were compared to those from the parametric bootstrap and from a Bayesian Markov chain Monte Carlo (MCMC) approach (implemented in the R package MCMCglmm). We apply each approach to evolvability statistics previously estimated for a large, 20-dimensional data set for Drosophila wings. REML-MVN and MCMC sampling variances are close to those estimated with the parametric bootstrap. Both slightly underestimate the error in the best-estimated aspects of the G matrix. REML analysis supports the previous conclusion that the G matrix for this population is full rank. REML-MVN is computationally very efficient, making it an attractive alternative to both data resampling and MCMC approaches to assessing confidence in parameters of evolutionary interest.

Response of canopy nitrogen uptake to a rapid decrease in bulk nitrate deposition in two eastern Canadian boreal forests.

A few studies have reported a recent and rapid decline in NO3(-) deposition in eastern North America. Whether this trend can be observed at remote boreal sites with low rates of N deposition and how it could impact canopy uptake (CU) of N remain unknown. Here we report trends between 1997/1999 and 2012 for precipitation, throughfall N deposition as well as inorganic N CU for two boreal forest sites of Quebec, Canada, with contrasted N deposition rates and tree species composition. NO3(-) bulk deposition declined by approximately 50% at both sites over the studied period while no change was observed for NH4(+). As a result, the contribution of NH4(+) to inorganic N deposition changed from ~33% to more than 50% during the study period. On average, 52-59% of N deposition was intercepted by the canopy, the retention being higher for NH4(+) (60-67%) than for NO3(-) (45-54%). The decrease in NO3(-) bulk deposition and the increase in the NH4(+):NO3(-) ratio had important impacts on N-canopy interactions. The contribution of NH4(+) CU to that of total inorganic N CU increased at both sites but the trend was significant only at Tirasse (lowest N deposition). At this site, absolute NO3(-) CU significantly decreased (as did total N CU) during the study period, a consequence of the strong relationship (r(2) = 0.88) between NO3(-) bulk deposition and NO3(-) CU. Our data suggest that N interactions with forest canopies may change rapidly with changes in N deposition as well as with tree species composition.

Tree species partition N uptake by soil depth in boreal forests.

It is recognized that the coexistence of herbaceous species in N-depleted habitats can be facilitated by N partitioning; however, the existence of such a phenomenon for trees has not yet been demonstrated. Here, we show from both foliage and soil 15N natural abundance values and from a 12-year in situ 15N addition experiment, that black spruce (Picea mariana) and jack pine (Pinus banksiana), two widespread species of the Canadian boreal forest, take up N at different depths. While black spruce takes up N from the organic soil, jack pine acquires it deeper within the highly N-depleted mineral soil. Systematic difference in foliar 15N natural abundance between the two species across seven sites distributed throughout the eastern Canadian boreal forest shows that N spatial partitioning is a widespread phenomenon. Distinct relationships between delta15N and N concentration in leaves of both species further emphasize their difference in N acquisition strategies. This result suggests that such complementary mechanisms of N acquisition could facilitate tree species coexistence in such N-depleted habitats and could contribute to the positive biodiversity-productivity relationship recently revealed for the eastern Canadian boreal forest, where jack pine is present. It also has implications for forest management and provides new insights to interpret boreal forest regeneration following natural or anthropogenic perturbations.

The effects of CO2 pressure and pH on the Suzuki coupling of basic nitrogen containing substrates.

The Suzuki coupling reaction of basic nitrogen containing substrates (2-bromo- and 2-chloro-4-aminopyridine, and 2-bromo and 2-chloropyridine) with phenylboronic acid using Pd(TPP)2Cl2/K3PO4 in acetonitrile-water biphasic solvent systems under a CO2 or a N2 atmosphere is discussed. It was observed that 2-halo-4-aminopyridine produced quantitative yields of coupled products under a CO2 atmosphere while the yields for the 2-halopyridines were poor. In contrast, the yields of coupled products for the 2-halopyridines substrates were quantitative under a N2 atmosphere while only poor yields were realized for the 2-halo-4-aminopyridines under the same conditions. Evidence is presented which suggests that the presence of CO2 alters the pH of the aqueous phase of the reaction system and the accompanying efficiency of the coupling process. Using a series of buffers to adjust the pH of the aqueous phase, the pH dependence associated with the efficiency of the coupling process is illustrated.

Key factors influencing Hg levels and trends in unperturbed oligotrophic temperate and boreal lakes.

Mercury (Hg) is a toxic metal that presents a major risk to ecosystems, biota, human health, and remains a priority concern. In temperate and boreal lakes Hg and methylmercury (MMHg) are expected to vary as a function of atmospheric Hg deposition, lake water chemistry, catchment characteristics and climate variables. The aim of this study was to quantify Hg and MMHg in unperturbed oligotrophic lakes and to identify the factors controlling their distribution. We first hypothesized that lake Hg (and MMHg to lesser extent) spatial variations are linked to atmospheric deposition, catchment characteristics, and terrestrial exportation of dissolved organic carbon (DOC). We secondly examined if lake Hg concentrations have followed the decrease in atmospheric Hg emission observed between the mid-1990s to the end-2010s. We found that overall, atmospheric Hg has little impact on lake Hg and MMHg concentrations, which are both primarily influenced by DOC input originating from the forest catchment. The relationship between DOC and Hg differed between the spring and the fall, with a Hg-to-DOC ratio twice as high in spring. This seems related to snowmelt input of Hg (with a relatively reduced input of DOC) or the internal lake build-up of Hg during the ice-covered period. Of the 10 lakes intensively visited over a 20-year period, only 3 showed significant lake Hg decreases despite significant negative trends in atmospheric Hg concentrations, suggesting a lag between atmospheric and surface water temporal trends. Overall, terrestrial catchments retain around 80% of atmospheric Hg implying that large Hg pools have been built up in soils in the last decades. As such, the reduction of atmospheric Hg alone will not necessarily result in Hg decreases in lakes, since the Hg concentrations may be modulated by DOC export trends and catchment characteristics. This stresses the need to improve our understanding of the processes governing Hg transfers from catchments into lakes.

Populations with elevated mutation load do not benefit from the operation of sexual selection.

Theory predicts that if most mutations are deleterious to both overall fitness and condition-dependent traits affecting mating success, sexual selection will purge mutation load and increase nonsexual fitness. We explored this possibility with populations of mutagenized Drosophila melanogaster exhibiting elevated levels of deleterious variation and evolving in the presence or absence of male-male competition and female choice. After 60 generations of experimental evolution, monogamous populations exhibited higher total reproductive output than polygamous populations. Parental environment also affected fitness measures - flies that evolved in the presence of sexual conflict showed reduced nonsexual fitness when their parents experienced a polygamous environment, indicating trans-generational effects of male harassment and highlighting the importance of a common garden design. This cost of parental promiscuity was nearly absent in monogamous lines, providing evidence for the evolution of reduced sexual antagonism. There was no overall difference in egg-to-adult viability between selection regimes. If mutation load was reduced by the action of sexual selection in this experiment, the resultant gain in fitness was not sufficient to overcome the costs of sexual antagonism.

The effect of inbreeding on fluctuating asymmetry of wing veins in two laboratory strains of Drosophila melanogaster.

Many authors have proposed that inbreeding destabilizes developmental processes. This destabilization may be reflected by increased fluctuating asymmetry (FA) in inbred compared to relatively outbred populations, but many studies have failed to find such differences. We measured the left and right wings of a large number of individual Drosophila melanogaster from two genetically distinct populations to estimate changes in FA caused by inbreeding. The large sample size and experimental design allowed removal of potentially confounding directional asymmetry (DA) and measurement error terms. Trait means in the two populations were essentially unchanged by inbreeding (less than 0.5% smaller in both populations). Inbred lines showed higher signed FA variances (16 and 38% higher, significantly so in one population) and higher unsigned FA means (3.7 and 13.2%, significantly increased in one population). Significant DA was found in both populations, although the pattern differed between populations. DA did not change due to inbreeding.

Measuring and comparing evolvability and constraint in multivariate characters.

The Lande equation forms the basis for our understanding of the short-term evolution of quantitative traits in a multivariate context. It predicts the response to selection as the product of an additive genetic variance matrix and a selection gradient. The selection gradient approximates the force and direction of selection, and the genetic variance matrix quantifies the role of the genetic system in evolution. Attempts to understand the evolutionary significance of the genetic variance matrix are hampered by the fact that the majority of the methods used to characterize and compare variance matrices have not been derived in an explicit theoretical context. We use the Lande equation to derive new measures of the ability of a variance matrix to allow or constrain evolution in any direction in phenotype space. Evolvability captures the ability of a population to evolve in the direction of selection when stabilizing selection is absent. Conditional evolvability captures the ability of a population to respond to directional selection in the presence of stabilizing selection on other trait combinations. We then derive measures of character autonomy and integration from these evolvabilities. We study the properties of these measures and show how they can be used to interpret and compare variance matrices. As an illustration, we show that divergence of wing shape in the dipteran family Drosophilidae has proceeded in directions that have relatively high evolvabilities.

Northeastern North America as a potential refugium for boreal forests in a warming climate.

High precipitation in boreal northeastern North America could help forests withstand the expected temperature-driven increase in evaporative demand, but definitive evidence is lacking. Using a network of tree-ring collections from 16,450 stands across 583,000 km(2) of boreal forests in Québec, Canada, we observe a latitudinal shift in the correlation of black spruce growth with temperature and reduced precipitation, from negative south of 49°N to largely positive to the north of that latitude. Our results suggest that the positive effect of a warmer climate on growth rates and growing season length north of 49°N outweighs the potential negative effect of lower water availability. Unlike the central and western portions of the continent's boreal forest, northeastern North America may act as a climatic refugium in a warmer climate.

Allozyme-associated heterosis in Drosophila melanogaster.

Two large experiments designed to detect allozyme-associated heterosis for growth rate in Drosophila melanogaster were performed. Heterosis associated with allozyme genotypes may be explained either by functional overdominance at the allozyme loci, or closely linked loci; or by genotypic correlations between allozyme loci and loci at which deleterious recessive alleles segregate. Such genotypic correlations would be favored by consanguineous mating, small effective population size, population mixing and strong natural or artificial selection. D. melanogaster is outbred, has large effective population size and there is little evidence for genotypic disequilibria. Therefore it would be unlikely to show allozyme heterosis due to genotypic correlations. In the first experiment I estimated the genotypic values of 97 replicated genotypes. In the second experiment, 500 individuals were raised in a fluctuating, stressful environment. In neither experiment was there any consistent evidence for allozyme heterosis in size or development rate, fluctuating asymmetry for size or in tendency to deviate from the population mean. In the first experiment, heterosis explained less than 5.6% of the genetic variance in growth characters. In the second, heterosis explained less than 0.1% of the phenotypic variance in growth characters. Outside of the molluscs, species which show allozyme heterosis have population structures or histories which tend to promote genotypic correlations. There is little evidence that functional overdominance is responsible for observations of allozyme-associated heterosis.

Comparing evolvability and variability of quantitative traits.

There are two distinct reasons for making comparisons of genetic variation for quantitative characters. The first is to compare evolvabilities, or ability to respond to selection, and the second is to make inferences about the forces that maintain genetic variability. Measures of variation that are standardized by the trait mean, such as the additive genetic coefficient of variation, are appropriate for both purposes. Variation has usually been compared as narrow sense heritabilities, but this is almost always an inappropriate comparative measure of evolvability and variability. Coefficients of variation were calculated from 842 estimates of trait means, variances and heritabilities in the literature. Traits closely related to fitness have higher additive genetic and nongenetic variability by the coefficient of variation criterion than characters under weak selection. This is the reverse of the accepted conclusion based on comparisons of heritability. The low heritability of fitness components is best explained by their high residual variation. The high additive genetic and residual variability of fitness traits might be explained by the great number of genetic and environmental events they are affected by, or by a lack of stabilizing selection to reduce their phenotypic variance. Over one-third of the quantitative genetics papers reviewed did not report trait means or variances. Researchers should always report these statistics, so that measures of variation appropriate to a variety of situations may be calculated.

Comparing mutational variabilities.

We have reviewed the available data on VM, the amount of genetic variation in phenotypic traits produced each generation by mutation. We use these data to make several qualitative tests of the mutation-selection balance hypothesis for the maintenance of genetic variance (MSB). To compare VM values, we use three dimensionless quantities: mutational heritability, VM/VE, the mutational coefficient of variation, CVM; and the ratio of the standing genetic variance to VM, VC/VM. Since genetic coefficients of variation for life history traits are larger than those for morphological traits, we predict that under MSB, life history traits should also have larger CVM. This is confirmed; life history traits have a median CVM value more than six times higher than that for morphological traits. VC/VM approximates the persistence time of mutations under MSB in an infinite population. In order for MSB to hold, VC/VM must be small, substantially less than 1000, and life history traits should have smaller values than morphological traits. VC/VM averages about 50 generations for life history traits and 100 generations for morphological traits. These observations are all consistent with the predictions of a mutation-selection balance model.

Age-specific properties of spontaneous mutations affecting mortality in Drosophila melanogaster.

An analysis of the effects of spontaneous mutations affecting age-specific mortality was conducted using 29 lines of Drosophila melanogaster that had accumulated spontaneous mutations for 19 generations. Divergence among the lines was used to estimate the mutational variance for weekly mortality rates and the covariance between weekly mortality rates at different ages. Significant mutational variance was observed in both males and females early in life (up to approximately 30 days of age). Mutational variance was not significantly different from zero for mortality rates at older ages. Mutational correlations between ages separated by 1 or 2 wk were generally positive, but they declined monotonically with increasing separation such that mutational effects on early-age mortality were uncorrelated with effects at later ages. Analyses of individual lines revealed several instances of mutation-induced changes in mortality over a limited range of ages. Significant age-specific effects of mutations were identified in early and middle ages, but surprisingly, mortality rates at older ages were essentially unaffected by the accumulation procedure. Our results provide strong evidence for the existence of a class of polygenic mutations that affect mortality rates on an age-specific basis. The patterns of mutational effects measured here relate directly to recently published estimates of standing genetic variance for mortality in Drosophila, and they support mutation accumulation as a viable mechanism for the evolution of senescence.

The evolution of age-specific mortality rates in Drosophila melanogaster: genetic divergence among unselected lines.

Age-specific effects of spontaneous mutations on mortality rates in Drosophila are inferred from three large demographic experiments. Data were collected from inbred lines that were allowed to accumulate spontaneous mutations for 10, 19, and 47 generations. Estimates of age-specific mutational variance for mortality were based on data from all three experiments, totalling approximately 225,000 flies, using a model developed for genetic analysis of age-dependent traits (the character process model). Both within- and among-generation analyses suggest that the input of genetic variance is greater for early life mortality rates than for mortality at older ages. In females, age-specific mutational variances ranged over an order of magnitude from 5.96 x 10(-3) at 2 wk posteclosion to 0.02 x 10(-3) at 7 wk. The male data show a similar pattern. Age-specific genetic variances were substantially less at generation 47 than at generation 19-an unexplained observation that is likely due to block effects. Mutational correlations among mortality rates at different ages tend to increase with the accumulation of new mutations. Comparison of the mutation-accumulation lines at generations 19 and 47 with their respective control lines suggests little age-specific mutational bias.

How informative is Wright's estimator of the number of genes affecting a quantitative character?

S. Wright suggested an estimator, m, of the number of loci, m, contributing to the difference in a quantitative character between two differentiated populations, which is calculated from the phenotypic means and variances in the two parental populations and their F1 and F2 hybrids. The same method can also be used to estimate m contributing to the genetic variance within a single population, by using divergent selection to create differentiated lines from the base population. In this paper we systematically examine the utility and problems of this technique under the influences of unequal allelic effects and initial allele frequencies, and linkage, which are known to lead m to underestimate m. In addition, we examine the effects of population size and selection intensity during the generations of selection. During selection, the estimator m rapidly approaches its expected value at the selection limit. With reasonable assumptions about unequal allelic effects and initial allele frequencies, the expected value of m without linkage is likely to be on the order of one-third of the number of genes. The estimates suffer most seriously from linkage. The practical maximum expectation of m is just about the number of chromosomes, considerably less than the "recombination index" which has been assumed to be the upper limit. The estimates are also associated with large sampling variances. An estimator of the variance of m derived by R. Lande substantially underestimates the actual variance. Modifications to the method can ameliorate some of the problems. These include using F3 or later generation variances or the genetic variance in the base population, and replicating the experiments and estimation procedure. However, even in the best of circumstances, information from m is very limited and can be misleading.

Homozygous and Hemizygous Viability Variation on the X Chromosome of DROSOPHILA MELANOGASTER.

We report here a study of viability inbreeding depression associated with the X chromosome of Drosophila melanogaster. Fifty wild chromosomes from Mt. Sinai, New York, and 90 wild chromosomes from Death Valley, California, were extracted using the marked FM6 balancer chromosome and viabilities measured for homozygous and heterozygous females, and for hemizygous males, relative to FM6 males as a standard genotype. No statistically significant female genetic load was observed for either chromosome set, although a 95% confidence limit estimated the total load <0.046 for the samples pooled. About 10% of the Death Valley chromosomes appear to be "supervital" as homozygotes. There is little evidence for a pervasive sex-limited detrimental load on the X chromosome; the evidence indicates nearly identical viability effects in males and homozygous females excluding the supervital chromosomes. The average degree of dominance for viability polygenes is estimated between 0.23 to 0.36, which is consistent with autosomal variation and implies near additivity. We conclude that there is little genetic load associated with viability variation on the X chromosome and that the substantial reduction in total fitness observed for chromosome homozygosity in an earlier study may be due largely to sex-limited fertility in females.

Targeted selection experiments and enzyme polymorphism: negative evidence for octanoate selection at the G6PD locus in Drosophila melanogaster.

Published studies have reported significant selection with respect to the G6pd locus for Drosophila melanogaster reared on Na-octanoate food. We have reexamined the selective effects of Na-octanoate on egg to adult viability with respect to the G6pd polymorphism using specially constructed X chromosomes. Four experiments were carried out using different 6Pgd backgrounds in two recombinant sets of chromosomes segregating for the G6pd locus but constructed so as to minimize variation over most of the X chromosome. In addition, two measures of viability were used, and the size of the experiments and their associated degrees of freedom are approximately double those reported in the former studies. Our results find no evidence for differential selection on G6pd genotypes (males and females) by Na-octanoate and, therefore, do not corroborate the positive results of selection reported by other investigators. The reasons for our different results are discussed.

The effects of spontaneous mutation on quantitative traits. I. Variances and covariances of life history traits.

We have accumulated spontaneous mutations in the absence of natural selection in Drosophila melanogaster by backcrossing 200 heterozygous replicates of a single high fitness second chromosome to a balancer stock for 44 generations. At generations 33 and 44 of accumulation, we extracted samples of chromosomes and assayed their homozygous performance for female fecundity early and late in adult life, male and female longevity, male mating ability early and late in adult life, productivity (a measure of fecundity times viability) and body weight. The variance among lines increased significantly for all traits except male mating ability and weight. The rate of increase in variance was similar to that found in previous studies of egg-to-adult viability, when calculated relative to trait means. The mutational correlations among traits were all strongly positive. Many correlations were significantly different from 0, while none was significantly different from 1. These data suggest that the mutation-accumulation hypothesis is not a sufficient explanation for the evolution of senescence in D. melanogaster. Mutation-selection balance does seem adequate to explain a substantial proportion of the additive genetic variance for fecundity and longevity.

Random depletion of T cells that bear specific T cell receptor V beta sequences in AIDS patients.

A cross-sectional PCR analysis of the TCR V beta repertoires in HIV-1 seronegative controls and HIV-1 infected individuals with either clinically or immunologically defined AIDS [1] was performed to examine the proposed superantigen model for HIV-1 pathogenesis. In contrast to previous reports, we find neither uniform specific losses nor uniform clonal expansions of particular TCR V beta gene families in subjects with AIDS. Instead our study, which was designed specifically to qualitatively determine the presence or absence of TCR V beta families in both subject populations, indicates an overall diminution in the expression of TCR V beta gene families in HIV-1 infected individuals with AIDS compared with controls. This is commensurate with the decrease in CD4 T cells in the AIDS population. Our data are therefore not directly suggestive of a common superantigen model of HIV-1 induced T cell clonal depletion or anergy, but instead emphasize a broad decrease in signals throughout the TCR V beta repertoire in AIDS versus control groups. This random depletion in the TCR V beta repertoire is most likely caused by aspects of HIV-1 pathogenesis other than virus-encoded superantigens.