Genomic insights into the secondary aquatic transition of penguins.

Penguins lost the ability to fly more than 60 million years ago, subsequently evolving a hyper-specialized marine body plan. Within the framework of a genome-scale, fossil-inclusive phylogeny, we identify key geological events that shaped penguin diversification and genomic signatures consistent with widespread refugia/recolonization during major climate oscillations. We further identify a suite of genes potentially underpinning adaptations related to thermoregulation, oxygenation, diving, vision, diet, immunity and body size, which might have facilitated their remarkable secondary transition to an aquatic ecology. Our analyses indicate that penguins and their sister group (Procellariiformes) have the lowest evolutionary rates yet detected in birds. Together, these findings help improve our understanding of how penguins have transitioned to the marine environment, successfully colonizing some of the most extreme environments on Earth.

Genetic evidence of hybridization between Magellanic (Sphensicus magellanicus) and Humboldt (Spheniscus humboldti) penguins in the wild.

The process of hybridization between closely related species plays an important role in defining the genetic integrity and overall genetic diversity of species. The distribution range of Magellanic (Spheniscus magellanicus) and Humboldt (Spheniscus humboldti) penguins is predominantly allopatric; however, the species share a region of sympatry where they may hybridize. We analyzed four types of genetic markers (including nuclear and mitochondrial markers) to assess their utility in detecting hybridization events between Magellanic and Humboldt penguins. Genetic assessment of non-introgressed reference samples allowed us to identify three types of informative markers (microsatellites, major histocompatibility complex, and mitochondrial DNA) and detect positive evidence of introgressive hybridization in the wild. Four out of six putative hybrids showed positive evidence of hybridization, revealed by the detection of Humboldt mitochondrial DNA and Magellanic species-specific alleles from nuclear markers. Bayesian Structure analysis, including samples from the sympatric region of the species in the southern Pacific Ocean, confirmed the use of nuclear markers for detecting hybridization and genetic admixture of putative hybrids, but revealed relatively low levels of genetic introgression at the population level. These findings provide insights into the role of hybridization in regions of species sympatry and its potential consequences on the levels of genetic introgression, genetic diversity, and conservation of these penguin species.

The role of demographic history and selection in shaping genetic diversity of the Galápagos penguin (Spheniscus mendiculus).

Although many studies have documented the effects of demographic bottlenecks on the genetic diversity of natural populations, there is conflicting evidence of the roles that genetic drift and selection may play in driving changes in genetic variation at adaptive loci. We analyzed genetic variation at microsatellite and mitochondrial loci in conjunction with an adaptive MHC class II locus in the Galápagos penguin (Spheniscus mendiculus), a species that has undergone serial demographic bottlenecks associated with El Niño events through its evolutionary history. We compared levels of variation in the Galápagos penguin to those of its congener, the Magellanic penguin (Spheniscus magellanicus), which has consistently maintained a large population size and thus was used as a non-bottlenecked control. The comparison of neutral and adaptive markers in these two demographically distinct species allowed assessment of the potential role of balancing selection in maintaining levels of MHC variation during bottleneck events. Our analysis suggests that the lack of genetic diversity at both neutral and adaptive loci in the Galápagos penguin likely resulted from its restricted range, relatively low abundance, and history of demographic bottlenecks. The Galápagos penguin revealed two MHC alleles, one mitochondrial haplotype, and six alleles across five microsatellite loci, which represents only a small fraction of the diversity detected in Magellanic penguins. Despite the decreased genetic diversity in the Galápagos penguin, results revealed signals of balancing selection at the MHC, which suggest that selection can mitigate some of the effects of genetic drift during bottleneck events. Although Galápagos penguin populations have persisted for a long time, increased frequency of El Niño events due to global climate change, as well as the low diversity exhibited at immunological loci, may put this species at further risk of extinction.

High-coverage genomes to elucidate the evolution of penguins.

BACKGROUND: Penguins (Sphenisciformes) are a remarkable order of flightless wing-propelled diving seabirds distributed widely across the southern hemisphere. They share a volant common ancestor with Procellariiformes close to the Cretaceous-Paleogene boundary (66 million years ago) and subsequently lost the ability to fly but enhanced their diving capabilities. With ∼20 species among 6 genera, penguins range from the tropical Galápagos Islands to the oceanic temperate forests of New Zealand, the rocky coastlines of the sub-Antarctic islands, and the sea ice around Antarctica. To inhabit such diverse and extreme environments, penguins evolved many physiological and morphological adaptations. However, they are also highly sensitive to climate change. Therefore, penguins provide an exciting target system for understanding the evolutionary processes of speciation, adaptation, and demography. Genomic data are an emerging resource for addressing questions about such processes. RESULTS: Here we present a novel dataset of 19 high-coverage genomes that, together with 2 previously published genomes, encompass all extant penguin species. We also present a well-supported phylogeny to clarify the relationships among penguins. In contrast to recent studies, our results demonstrate that the genus Aptenodytes is basal and sister to all other extant penguin genera, providing intriguing new insights into the adaptation of penguins to Antarctica. As such, our dataset provides a novel resource for understanding the evolutionary history of penguins as a clade, as well as the fine-scale relationships of individual penguin lineages. Against this background, we introduce a major consortium of international scientists dedicated to studying these genomes. Moreover, we highlight emerging issues regarding ensuring legal and respectful indigenous consultation, particularly for genomic data originating from New Zealand Taonga species. CONCLUSIONS: We believe that our dataset and project will be important for understanding evolution, increasing cultural heritage and guiding the conservation of this iconic southern hemisphere species assemblage.

Mitogenomes Uncover Extinct Penguin Taxa and Reveal Island Formation as a Key Driver of Speciation.

The emergence of islands has been linked to spectacular radiations of diverse organisms. Although penguins spend much of their lives at sea, they rely on land for nesting, and a high proportion of extant species are endemic to geologically young islands. Islands may thus have been crucial to the evolutionary diversification of penguins. We test this hypothesis using a fossil-calibrated phylogeny of mitochondrial genomes (mitogenomes) from all extant and recently extinct penguin taxa. Our temporal analysis demonstrates that numerous recent island-endemic penguin taxa diverged following the formation of their islands during the Plio-Pleistocene, including the Galápagos (Galápagos Islands), northern rockhopper (Gough Island), erect-crested (Antipodes Islands), Snares crested (Snares) and royal (Macquarie Island) penguins. Our analysis also reveals two new recently extinct island-endemic penguin taxa from New Zealand's Chatham Islands: Eudyptes warhami sp. nov. and a dwarf subspecies of the yellow-eyed penguin, Megadyptes antipodes richdalei ssp. nov. Eudyptes warhami diverged from the Antipodes Islands erect-crested penguin between 1.1 and 2.5 Ma, shortly after the emergence of the Chatham Islands (∼3 Ma). This new finding of recently evolved taxa on this young archipelago provides further evidence that the radiation of penguins over the last 5 Ma has been linked to island emergence. Mitogenomic analyses of all penguin species, and the discovery of two new extinct penguin taxa, highlight the importance of island formation in the diversification of penguins, as well as the extent to which anthropogenic extinctions have affected island-endemic taxa across the Southern Hemisphere's isolated archipelagos.

Spatial Patterns of bphA Gene Diversity Reveal Local Adaptation of Microbial Communities to PCB and PAH Contaminants.

Biphenyl dioxygenases, encoded by the bphA gene, initiate the oxidation of polychlorinated biphenyls (PCBs) and specify the substrate range of PCB congeners metabolized by bacteria. Increased bphA gene diversity within microbial communities may allow a broader range of PCB congeners to be catabolized, thus resulting in greater PCB degradation. To assess the role of PCBs in modulating bphA gene diversity, 16S ribosomal RNA (rRNA) gene and bphA environmental DNA libraries were generated from bacterial communities in sediments with a steep gradient of PCB contamination. Multiple measures of sequence diversity revealed greater heterogeneity of bphA sequences in polluted compared to unpolluted locations. Codon-based signatures of selection in bphA sequences provided evidence of purifying selection. Unifrac analysis of 16S rRNA sequences revealed independent taxonomic lineages from polluted and unpolluted locations, consistent with the presence of locally adapted bacterial communities. Phylogenetic analysis of bphA sequences indicated that dioxygenases from sediments were closely related to previously characterized dioxygenases that metabolize PCBs and polynuclear aromatic hydrocarbons (PAHs), consistent with high levels of these contaminants within the studied sediments. Structural analyses indicated that the BphA protein of Rhodococcus jostii, capable of metabolizing both PCBs and PAHs, provided a more optimal modeling template for bphA sequences reported in this study than a BphA homologue with more restricted substrate specificity. Results from this study suggest that PCBs and PAHs may drive local adaptation of microbial communities by acting as strong selective agents for biphenyl dioxygenases capable of metabolizing a wide range of congeners.

A perspective on the evolution of germ-cell development and germinal mosaics of deleterious mutations.

In many animals a small number of primordial germ cells (PGCs) are set aside early in development, mitosis and mitochondrial DNA syntheses are arrested, transcription is stopped or reduced, and the PGCs migrate later to the emerging gonads and become germ cells. What could be the evolutionary advantage of sequestering non-dividing PGCs early in development? A commonly cited advantage is a reduction in the number of new deleterious mutations that would occur if there were additional divisions in PGCs early in development. We would like to add to this advantage the fact that these additional mutations in PGCs give rise to germinal mosaics (i.e., premeiotic clusters of mutation) in multiple progeny of the same individual, thus having a larger detrimental effect on the evolutionary fitness of their carriers. Here, we reviewed published studies providing evidence that germinal mosaics of deleterious mutant alleles are not rare, occur for all types of genetic damage, and have been observed in all tested organisms and in nature. We propose the hypothesis that PGC sequestration during early animal development may have evolved in part in response to selection for preventing the occurrence of premeiotic clusters of deleterious mutant alleles, and describe a series of predictions that would allow the assessment of the potential role of germinal mosaics on the evolution of PGC sequestration.

Characterization of major histocompatibility complex class I loci of the lark sparrow (Chondestes grammacus) and insights into avian MHC evolution.

The major histocompatibilty complex (MHC) has become increasingly important in the study of the immunocapabilities of non-model vertebrates due to its direct involvement in the immune response. The characterization of MHC class I loci in the lark sparrow (Chondestes grammacus) revealed multiple MHC class I loci with elevated genetic diversity at exon 3, evidence of differential selection between the peptide binding region (PBR) and non-PBR, and the presence of multiple pseudogenes with limited divergence. The minimum number of functional MHC class I loci was estimated at four. Sequence analysis revealed d N /d S ratios significantly less than one at non-PBR sites, indicative of negative selection, whereas PBR sites associated with antigen recognition showed ratios greater than 1 but non-significant. GenBank surveys and phylogenetic analyses of previously reported avian MHC class I sequences revealed variable signatures of evolutionary processes acting upon this gene family, including gene duplication and potential concerted evolution. An increase in the number of class I loci across species coincided with an increase in pseudogene prevalence, revealing the importance of gene duplication in the expansion of multigene families and the creation of pseudogenes.

Distribution, prevalence and host specificity of avian malaria parasites across the breeding range of the migratory lark sparrow (Chondestes grammacus).

The lark sparrow (Chondestes grammacus) is a ground-nesting passerine that breeds across much of the central North American steppe and sand barrens. Through genotyping and sequencing of avian malaria parasites we examined levels of malaria prevalence and determined the distribution of Haemoproteus and Plasmodium lineages across the breeding range of the lark sparrow. Analysis of 365 birds collected from five breeding locations revealed relatively high levels of malaria prevalence in adults (80 %) and juveniles (46 %), with infections being primarily of Haemoproteus (91 % of sequenced samples). Levels of genetic diversity and genetic structure of malaria parasites with respect to the avian host populations revealed distinct patterns for Haemoproteus and Plasmodium, most likely as a result of their distinct life histories, host specificity, and transmission vectors. With the exception of one common Haemoproteus haplotype detected in all populations, all other haplotypes were either population-specific or shared by two to three populations. A hierarchical analysis of molecular variance of Haemoproteus sequences revealed that 15-18 % of the genetic variation can be explained by differences among host populations/locations (p < 0.001). In contrast to the regional patterns of genetic differentiation detected for the lark sparrow populations, Haemoproteus parasites showed high levels of population-specific variation and no significant differences among regions, which suggests that the population dynamics of the parasites may be driven by evolutionary processes operating at small spatial scales (e.g., at the level of host populations). These results highlight the potential effects of host population structure on the demographic and evolutionary dynamics of parasites.

Darwin's diagram of divergence of taxa as a causal model for the origin of species.

On the basis that Darwin's theory of evolution encompasses two logically independent processes (common descent and natural selection), the only figure in On the Origin of Species (the Diagram of Divergence of Taxa) is often interpreted as illustrative of only one of these processes: the branching patterns representing common ancestry. Here, I argue that Darwin's Diagram of Divergence of Taxa represents a broad conceptual model of Darwin's theory, illustrating the causal efficacy of natural selection in producing well-defined varieties and ultimately species. The Tree Diagram encompasses the idea that natural selection explains common descent and the origin of organic diversity, thus representing a comprehensive model of Darwin's theory on the origin of species. I describe Darwin's Tree Diagram in relation to his argumentative strategy under the vera causa principle, and suggest that the testing of his theory based on the evidence from the geological record, the geographical distribution of organisms, and the mutual affinities of organic beings can be framed under the hypothetico-deductive method. Darwin's Diagram of Divergence of Taxa therefore represents a broad conceptual model that helps understanding the causal construction of Darwin's theory of evolution, the structure of his argumentative strategy, and the nature of his scientific methodology.

Evolutionary analysis and lateral gene transfer of two-component regulatory systems associated with heavy-metal tolerance in bacteria.

Microorganisms have adapted intricate signal transduction mechanisms to coordinate tolerance to toxic levels of metals, including two-component regulatory systems (TCRS). In particular, both cop and czc operons are regulated by TCRS; the cop operon plays a key role in bacterial tolerance to copper, whereas the czc operon is involved in the efflux of cadmium, zinc, and cobalt from the cell. Although the molecular physiology of heavy metal tolerance genes has been extensively studied, their evolutionary relationships are not well-understood. Phylogenetic relationships among heavy-metal efflux proteins and their corresponding two-component regulatory proteins revealed orthologous and paralogous relationships from species divergences and ancient gene duplications. The presence of heavy metal tolerance genes on bacterial plasmids suggests these genes may be prone to spread through horizontal gene transfer. Phylogenetic inferences revealed nine potential examples of lateral gene transfer associated with metal efflux proteins and two examples for regulatory proteins. Notably, four of the examples suggest lateral transfer across major evolutionary domains. In most cases, differences in GC content in metal tolerance genes and their corresponding host genomes confirmed lateral gene transfer events. Three-dimensional protein structures predicted for the response regulators encoded by cop and czc operons showed a high degree of structural similarity with other known proteins involved in TCRS signal transduction, which suggests common evolutionary origins of functional phenotypes and similar mechanisms of action for these response regulators.

MHC diversity and mate choice in the magellanic penguin, Spheniscus magellanicus.

We estimated levels of diversity at the major histocompatibility complex (MHC) class II DRß1 gene in 50 breeding pairs of the Magellanic penguin and compared those to estimates from Humboldt and Galapagos penguins. We tested for positive selection and 2 conditions required for the evolution of MHC-based disassortative mating: 1) greater MHC diversity between breeding pairs compared to random mating, and 2) associations between MHC genotype and fitness. Cloning and sequencing of the DRß1 gene showed that Magellanic penguins had higher levels of genetic variation than Galapagos and Humboldt penguins. Sequence analysis revealed 45 alleles with 3.6% average proportion of nucleotide differences, nucleotide diversity of 0.030, and observed heterozygosity of 0.770. A gene phylogeny showed 9 allelic lineages with interspersed DRß1 sequences from Humboldt and Galapagos penguins, indicating ancestral polymorphisms. d (N)/d (S) ratios revealed evidence for positive selection. Analysis of breeding pairs showed no disassortative mating preferences. Significant MHC genotype/fitness associations in females suggest, however, that selection for pathogen resistance plays a more important role than mate choice in maintaining diversity at the MHC in the Magellanic penguin. The differential effect of MHC heterozygosity on fitness between the sexes is likely associated with the relative role of hatching and fledging rates as reliable indicators of overall fitness in males and females.

Local adaptation of microbial communities to heavy metal stress in polluted sediments of Lake Erie.

Microbial communities must balance the assimilation of biologically necessary metals with resistance to toxic metal concentrations. To investigate the impact of heavy metal contaminants on microbial communities, we performed two experiments measuring extracellular enzyme activities (EEA) in polluted and unpolluted sediments of Lake Erie. In the first experiment, inoculations with moderate concentrations of copper and zinc appreciably diminished EEA from uncontaminated sites, whereas EEA from contaminated sediments increased or were only negligibly affected. In the second experiment, we compared the effects of three separate metals (i.e. copper, arsenic, and cadmium) on microbial community metabolism in polluted and unpolluted locations. Although copper and arsenic elicited differential effects by inhibiting EEA only in unpolluted sediments, cadmium inhibited EEA in both polluted and unpolluted sediments. Multivariate analyses of EEA from polluted sediments revealed direct associations among hydrolytic enzymes and inverse or absent associations between hydrolases and oxidases; these associations demonstrated resilience to heavy metal stress. In contrast, addition of heavy metals to unpolluted sediments appeared to have a higher impact on the multivariate pattern of EEA associations as revealed by an increase in the number of associations, more inverse relationships, and potential enzymatic trade-offs. The results of this study suggest community-level adaptations through the development of resistance mechanisms to the types and local levels of heavy metals in the environment.

The modulating role of dissolved organic matter on spatial patterns of microbial metabolism in Lake Erie sediments.

To evaluate the role of dissolved organic matter (DOM) on microbial community metabolism, we established extracellular enzyme activity (EEA) and substrate-induced respiration (SIR) profiles of sediment samples collected from littoral and profundal regions of the western, central, and eastern basins of Lake Erie. Lake Erie is spatially structured such that the central and western basins receive relatively major inputs of allochthonous DOM in comparison to the eastern basin. Overall, spatial patterns of EEA and SIR profiles suggest both greater metabolic diversity and activity in the littoral regions of the central and western basins. In contrast, the eastern basin demonstrated much less structuring between littoral and profundal areas. To evaluate whether the observed spatial patterns are the result of microbial community adaptations to local DOM availability, we performed three experimental treatments by inoculating sediment samples with polyvinylpyrrolidone, which sequesters large polyphenols, or with either vanillin or catechol, two small phenolic compounds. Our results revealed that esterase and glycosidase EEA from the eastern basin were induced by small phenolics and inhibited by large polyphenols. In contrast, the addition of small phenolics decreased esterase and glycosidase activities from the central basin, while polyphenols had a negligible effect. These results suggest that the source and composition of DOM play a significant role in the local adaptation of microbial communities, determining large-scale spatial patterns of microbial functional diversity in Lake Erie sediments.

Effects of recent population bottlenecks on reconstructing the demographic history of prairie-chickens.

Current methods of DNA sequence analysis attempt to reconstruct historical patterns of population structure and growth from contemporary samples. However, these techniques may be influenced by recent population bottlenecks, which have the potential to eliminate lineages that reveal past changes in demography. One way to examine the performance of these demographic methods is to compare samples from populations before and after recent bottlenecks. We compared estimates of demographic history from populations of greater prairie-chickens (Tympanuchus cupido) before and after recent bottlenecks using four common methods (nested clade analysis [NCA], Tajima's D, mismatch distribution, and MDIV). We found that NCA did not perform well in the presence of bottleneck events, although it did recover some genetic signals associated with increased isolation and the extinction of intermediate populations. The majority of estimates for Tajima's D, including those from bottlenecked populations, were not significantly different from zero, suggesting our data conformed to neutral expectations. In contrast, mismatch distributions including the raggedness index were more likely to identify recently bottlenecked populations with this data set. Estimates of population mutation rate (theta), population divergence time (t), and time to the most recent common ancestor (TMRCA) from MDIV were similar before and after bottlenecks; however, estimates of gene flow (M) were significantly lower in a few cases following a bottleneck. These results suggest that caution should be used when assessing demographic history from contemporary data sets, as recently fragmented and bottlenecked populations may have lost lineages that affect inferences of their demographic history.

Selection and inbreeding depression: effects of inbreeding rate and inbreeding environment.

The magnitude of inbreeding depression in small populations may depend on the effectiveness with which natural selection purges deleterious recessive alleles from populations during inbreeding. The effectiveness of this purging process, however, may be influenced by the rate of inbreeding and the environment in which inbreeding occurs. Although some experimental studies have examined these factors individually, no study has examined their joint effect or potential interaction. In the present study, therefore, we performed an experiment in which 180 lineages of Drosophila melanogaster were inbred at slow and fast inbreeding rates within each of three inbreeding environments (benign, high temperature, and competitive). The fitness of all lineages was then measured in a common benign environment. Although slow inbreeding reduced inbreeding depression in lineages inbred under high temperature stress, a similar reduction was not observed with respect to the benign or competitive treatments. Overall, therefore, the effect of inbreeding rate was nonsignificant. The inbreeding environment, in contrast, had a larger and more consistent effect on inbreeding depression. Under both slow and fast rates of inbreeding, inbreeding depression was significantly reduced in lineages inbred in the presence of a competitor D. melanogaster strain. A similar reduction of inbreeding depression occurred in lineages inbred under high temperature stress at a slow inbreeding rate. Overall, our findings show that inbreeding depression is reduced when inbreeding takes place in a stressful environment, possibly due to more effective purging under such conditions.

Associations between environmental stress, selection history, and quantitative genetic variation in Drosophila melanogaster.

Stressful environments may increase quantitative genetic variation in populations by promoting the expression of genetic variation that has not previously been eliminated or canalized by natural selection. This "selection history" hypothesis predicts that novel stressors will increase quantitative genetic variation, and that the magnitude of this effect will decrease following continued stress exposure. We tested these predictions using Drosophila melanogaster and sternopleural bristle number as a model system. In particular, we examined the effect of high temperature stress (31 degrees Celsius) on quantitative genetic variation before and after our study population had been reared at 31 degrees Celsius for 15 generations. High temperature stress was found to increase both additive genetic variance and heritability, but contrary to the selection history hypothesis prediction, the magnitude of this effect significantly increased after the study population had been reared for 15 generations under high temperature stress. These results demonstrate that high temperature stress increases quantitative genetic variation for bristle number, but do not support the selection history hypothesis as an explanation for this effect.

Ancestral inbreeding reduces the magnitude of inbreeding depression in Drosophila melanogaster.

The influence of natural selection on the magnitude of inbreeding depression is an important issue in conservation biology and the study of evolution. It is generally expected that the magnitude of inbreeding depression in small populations will depend upon the average homozygosity of individuals, as measured by the coefficient of inbreeding (F). However, if deleterious recessive alleles are selectively purged from populations during inbreeding, then inbreeding depression may differ among populations in which individuals have the same inbreeding coefficient. In such cases, the magnitude of inbreeding depression will partly depend on the ancestral inbreeding coefficient (fa), which measures the cumulative proportion of loci that have historically been homozygous and therefore exposed to natural selection. We examined the inbreeding depression that occurred in lineages of Drosophila melanogaster maintained under pedigrees that led to the same inbreeding coefficient (F = 0.375) but different levels of ancestral inbreeding (fa = 0.250 or 0.531). Although inbreeding depression varied substantially among individual lineages, we observed a significant 40% decrease in the median level of inbreeding depression in the treatment with higher ancestral inbreeding. Our results demonstrate that high levels of ancestral inbreeding are associated with greater purging effects, which reduces the inbreeding depression that occurs in isolated populations of small size.

Inbreeding depression and male survivorship in Drosophila: implications for senescence theory.

The extent to which inbreeding depression affects longevity and patterns of survivorship is an important issue from several research perspectives, including evolutionary biology, conservation biology, and the genetic analysis of quantitative traits. However, few previous inbreeding depression studies have considered longevity as a focal life-history trait. We maintained laboratory populations of Drosophila melanogaster at census population sizes of 2 and 10 male-female pairs for up to 66 generations and performed repeated assays of male survivorship throughout this time period. On average, significant levels of inbreeding depression were observed for median life span and age-specific mortality. For age-specific mortality, the severity of inbreeding depression increased over the life span. We found that a baseline inbreeding load of 0.307 lethal equivalents per gamete affected age-specific mortality, and that this value increased at a rate of 0.046 per day of the life span. With respect to some survivorship parameters, the differentiation of lineages was nonlinear with respect to the inbreeding coefficient, which suggested that nonadditive genetic variation contributed to variation among lineages. These findings provide insights into the genetic basis of longevity as a quantitative trait and have implications regarding the mutation-accumulation evolutionary explanation of senescence.