The thermal environment of sperm affects offspring success: a test of the anticipatory paternal effects hypothesis in the blue mussel.

There has been an explosion of recent evidence that environments experienced by fathers or their ejaculates can influence offspring phenotypes (paternal effects). However, little is known about whether such effects are adaptive, which would have far-reaching implications for the many species facing rapidly changing environments. For example, some arguments suggest paternal effects might be a source of cross-generational plasticity, preparing offspring to face similar conditions to their father (anticipatory hypothesis). Alternatively, ejaculate-mediated effects on offspring may be non-adaptive by-products of stress. Here, we conduct an experiment to distinguish between these predictions, exposing ejaculates of the externally fertilizing mussel Mytilus galloprovincialis to ambient (19°C) and high (24°C) temperatures, then rearing offspring groups in temperatures that match and mismatch those of sperm. We find that, overall, high temperature-treated sperm induced higher rates of normal offspring development and higher success in transitioning to second-stage larvae, which may represent adaptive epigenetic changes or selection on sperm haplotypes. However, the progeny of high temperature-treated sperm did not perform better than those of ambient temperature-treated sperm when rearing temperatures were high. Overall, these findings offer little support for the anticipatory hypothesis and suggest instead that beneficial paternal effects may be eroded when offspring develop under stressful conditions.

Heritable variation in heat shock gene expression: a potential mechanism for adaptation to thermal stress in embryos of sea turtles.

The capacity of species to respond adaptively to warming temperatures will be key to their survival in the Anthropocene. The embryos of egg-laying species such as sea turtles have limited behavioural means for avoiding high nest temperatures, and responses at the physiological level may be critical to coping with predicted global temperature increases. Using the loggerhead sea turtle (Caretta caretta) as a model, we used quantitative PCR to characterise variation in the expression response of heat-shock genes (hsp60, hsp70 and hsp90; molecular chaperones involved in cellular stress response) to an acute non-lethal heat shock. We show significant variation in gene expression at the clutch and population levels for some, but not all hsp genes. Using pedigree information, we estimated heritabilities of the expression response of hsp genes to heat shock and demonstrated both maternal and additive genetic effects. This is the first evidence that the heat-shock response is heritable in sea turtles and operates at the embryonic stage in any reptile. The presence of heritable variation in the expression of key thermotolerance genes is necessary for sea turtles to adapt at a molecular level to warming incubation environments.

Increased expression of Hsp70 and Hsp90 mRNA as biomarkers of thermal stress in loggerhead turtle embryos (Caretta Caretta).

The survival and viability of sea turtle embryos is dependent upon favourable nest temperatures throughout the incubation period. Consequently, future generations of sea turtles may be at risk from increasing nest temperatures due to climate change, but little is known about how embryos respond to heat stress. Heat shock genes are likely to be important in this process because they code for proteins that prevent cellular damage in response to environmental stressors. This study provides the first evidence of an expression response in the heat shock genes of embryos of loggerhead sea turtles (Caretta caretta) exposed to realistic and near-lethal temperatures (34°C and 36°C) for 1 or 3 hours. We investigated changes in Heat shock protein 60 (Hsp60), Hsp70, and Hsp90 mRNA in heart (n=24) and brain tissue (n=29) in response to heat stress. Under the most extreme treatment (36°C, 3h), Hsp70 increased mRNA expression by a factor of 38.8 in heart tissue and 15.7 in brain tissue, while Hsp90 mRNA expression increased by a factor of 98.3 in heart tissue and 14.7 in brain tissue. Hence, both Hsp70 and Hsp90 are useful biomarkers for assessing heat stress in the late-stage embryos of sea turtles. The method we developed can be used as a platform for future studies on variation in the thermotolerance response from the clutch to population scale, and can help us anticipate the resilience of reptile embryos to extreme heating events.

Exploring Symbiodinium diversity and host specificity in Acropora corals from geographical extremes of Western Australia with 454 amplicon pyrosequencing.

Scleractinian corals have demonstrated the ability to shuffle their endosymbiotic dinoflagellate communities (genus Symbiodinium) during periods of acute environmental stress. This has been proposed as a mechanism of acclimation, which would be increased by a diverse and flexible association with Symbiodinium. Conventional molecular techniques used to evaluate Symbiodinium diversity are unable to identify genetic lineages present at background levels below 10%. Next generation sequencing (NGS) offers a solution to this problem and can resolve microorganism diversity at much finer scales. Here we apply NGS to evaluate Symbiodinium diversity and host specificity in Acropora corals from contrasting regions of Western Australia. The application of 454 pyrosequencing allowed for detection of Symbiodinium operational taxonomic units (OTUs) occurring at frequencies as low as 0.001%, offering a 10,000-fold increase in sensitivity compared to traditional methods. All coral species from both regions were overwhelmingly dominated by a single clade C OTU (accounting for 98% of all recovered sequences). Only 8.5% of colonies associated with multiple clades (clades C and D, or C and G), suggesting a high level of symbiont specificity in Acropora assemblages in Western Australia. While only 40% of the OTUs were shared between regions, the dominance of a single OTU resulted in no significant difference in Symbiodinium community structure, demonstrating that the coral-algal symbiosis can remain stable across more than 15° of latitude and a range of sea surface temperature profiles. This study validates the use of NGS platforms as tools for providing fine-scale estimates of Symbiodinium diversity and can offer critical insight into the flexibility of the coral-algal symbiosis.

Resilience to periodic disturbances and the long-term genetic stability in Acropora coral.

Climate change is restructuring natural ecosystems. The direct impacts of these events on biodiversity and community structure are widely documented, but the impacts on the genetic variation of populations remains largely unknown. We monitored populations of Acropora coral on a remote coral reef system in northwest Australia for two decades and through multiple cycles of impact and recovery. We combined these demographic data with a temporal genetic dataset of a common broadcast spawning corymbose Acropora to explore the spatial and temporal patterns of connectivity underlying recovery. Our data show that broad-scale dispersal and post-recruitment survival drive recovery from recurrent disturbances, including mass bleaching and mortality. Consequently, genetic diversity and associated patterns of connectivity are maintained through time in the broader metapopulation. The results highlight an inherent resilience in these globally threatened species of coral and showcase their ability to cope with multiple disturbances, given enough time to recover is permitted.

Estimating relatedness and inbreeding using molecular markers and pedigrees: the effect of demographic history.

Estimates of inbreeding and relatedness are commonly calculated using molecular markers, although the accuracy of such estimates has been questioned. As a further complication, in many situations, such estimates are required in populations with reduced genetic diversity, which is likely to affect their accuracy. We investigated the correlation between microsatellite- and pedigree-based coefficients of inbreeding and relatedness in laboratory populations of Drosophila melanogaster that had passed through bottlenecks to manipulate their genetic diversity. We also used simulations to predict expected correlations between marker- and pedigree-based estimates and to investigate the influence of linkage between loci and null alleles. Our empirical data showed lower correlations between marker- and pedigree-based estimates in our control (nonbottleneck) population than were predicted by our simulations or those found in similar studies. Correlations were weaker in bottleneck populations, confirming that extreme reductions in diversity can compromise the ability of molecular estimates to detect recent inbreeding events. However, this result was highly dependent on the strength of the bottleneck and we did not observe or predict any reduction in correlations in our population that went through a relatively severe bottleneck of N = 10 for one generation. Our results are therefore encouraging, as molecular estimates appeared robust to quite severe reductions in genetic diversity. It should also be remembered that pedigree-based estimates may not capture realized identity-by-decent and that marker-based estimates may actually be more useful in certain situations.

Assortative mating for relatedness in a large naturally occurring population of Drosophila melanogaster.

New theoretical work on kin selection and inclusive fitness benefits predicts that individuals will sometimes choose close or intermediate relatives as mates to maximize their fitness. However, empirical examples supporting such predictions are rare. In this study, we look for such evidence in a natural population of Drosophila melanogaster. We compared mating and nonmating individuals to test whether mating was nonrandom with respect to relatedness. Consistent with optimal inbreeding, males were more closely related to their mate than to randomly sampled females. However, all individuals collected mating showed higher relatedness and males were not significantly more related to their mate than to other mating females. We also found a negative relationship between relatedness and fecundity. Our results are consistent with the hypothesis that inclusive fitness benefits may drive inbreeding tolerance despite direct costs to fitness; however, an experimental approach is needed to investigate the link between mate preference and relatedness.

Nonclinality of molecular variation implicates selection in maintaining a morphological cline of Drosophila melanogaster.

One general approach for assessing whether phenotypic variation is due to selection is to test its correlation with presumably neutral molecular variation. Neutral variation is determined by population history, the most likely alternative explanation of spatial genetic structure, whereas phenotypic variation may be influenced by the spatial pattern of selection pressure. Several methods for comparing the spatial apportionment of molecular and morphological variation have been used. Here, we present an analysis of variance framework that compares the magnitudes of latitudinal effects for molecular and morphological variation along a body size cline in Australian Drosophila populations. Explicit incorporation of the relevant environmental gradient can result in a simple and powerful test of selection. For the Australian cline, our analysis provides strong internal evidence that the cline is due to selection.

Characterization of microsatellite loci for the littorine snail Bembicium vittatum.

We describe the isolation and development of 17 polymorphic microsatellite loci for the intertidal snail Bembicium vittatum (Gastropoda: Littorinidae). The loci were tested in 46 individuals from a single population situated near the centre of the species distribution. No evidence of linkage disequilibrium was detected between any pair of loci. However, two loci showed significant departures from Hardy-Weinberg expectations. The number of alleles per locus ranged from two to 15.

Lack of genetic structure among ecologically adapted populations of an Australian rainforest Drosophila species as indicated by microsatellite markers and mitochondrial DNA sequences.

Although fragmented rainforest environments represent hotspots for invertebrate biodiversity, few genetic studies have been conducted on rainforest invertebrates. Thus, it is not known if invertebrate species in rainforests are highly genetically fragmented, with the potential for populations to show divergent selection responses, or if there are low levels of gene flow sufficient to maintain genetic homogeneity among fragmented populations. Here we use microsatellite markers and DNA sequences from the mitochondrial ND5 locus to investigate genetic differences among Drosophila birchii populations from tropical rainforests in Queensland, Australia. As found in a previous study, mitochondrial DNA diversity was low with no evidence for population differentiation among rainforest fragments. The pattern of mitochondrial haplotype variation was consistent with D. birchii having undergone substantial past population growth. Levels of nuclear genetic variation were high in all populations while F(ST) values were very low, even for flies from geographically isolated areas of rainforest. No significant differentiation was observed between populations on either side of the Burdekin Gap (a long-term dry corridor), although there was evidence for higher gene diversity in low-latitude populations. Spatial autocorrelation coefficients were low and did not differ significantly from random, except for one locus which revealed a clinal-like pattern. Comparisons of microsatellite differentiation contrasted with previously established clinal patterns in quantitative traits in D. birchii, and indicate that the patterns in quantitative traits are likely to be due to selection. These results suggest moderate gene flow in D. birchii over large distances. Limited population structure in this species appears to be due to recent range expansions or cycles of local extinctions followed by recolonizations/expansions. Nevertheless, patterns of local adaptation have developed in D. birchii that may result in populations showing different selection responses when faced with environmental change.

The genetic bases of divergence in desiccation and starvation resistance among tropical and temperate populations of Drosophila melanogaster.

Desiccation and starvation resistance are two stress-related traits which vary geographically with climate in Drosophila melanogaster. To investigate the contribution of epistasis to population divergence for these traits, we crossed tropical and temperate populations from two different geographical regions to produce F1, F2 and first backcross generations. Line-cross analysis of generation means revealed that genetic bases of divergence for both traits were complex and remarkably similar in a number of respects. Strong additive and dominance effects were present in most of the models, whereas epistatic and maternal effects were less common. The presence of epistatic effects in approximately half of the models presented in this study is consistent with line-cross studies of diverged traits in other animals, and does not support the view that epistasis is the predominant means by which populations diverge. In addition, evidence of maternal effects in both traits adds to a growing body of recent evidence that suggests maternal contributions to population differentiation are more widespread than previously thought. This finding undermines the accuracy of studies inferring epistasis directly from the magnitude of F2 breakdown. More line-cross analysis studies of naturally diverged populations that take into account maternal effects will shed further light on the true incidence of epistasis and its importance in the evolutionary process.

Quantitative genetic analysis of natural variation in body size in Drosophila melanogaster.

Latitudinal, genetic variation in body size is a commonly observed phenomenon in many invertebrate species and is shaped by natural selection. In this study, we use a chromosome substitution and a quantitative trait locus (QTL) mapping approach to identify chromosomes and genomic regions associated with adaptive variation in body size in natural populations of Drosophila melanogaster from the extreme ends of clines in South America and Australia. Chromosome substitution revealed the largest effects on chromosome three in both continents, and minor effects on the X and second chromosome. Similarly, QTL analysis of the Australian cline identified QTL with largest effects on the third chromosome, with smaller effects on the second. However, no QTL were found on the X chromosome. We also compared the coincidence of locations of QTL with the locations of five microsatellite loci previously shown to vary clinally in Australia. Permutation tests using both the sum of the LOD scores and the sum distance to nearest QTL peak revealed there were no significant associations between locations of clinal markers and QTL's. The lack of significance may, in part, be due to broad QTL peaks identified in this study. Future studies using higher resolution QTL maps should reveal whether the degree of clinality in microsatellite allele frequencies can be used to identify QTL in traits that vary along an environmental gradient.