Correction: A novel power-amplified jumping behavior in larval beetles (Coleoptera: Laemophloeidae).

[This corrects the article DOI: 10.1371/journal.pone.0256509.].

A novel method for the establishment of autologous skin cell suspensions: characterisation of cellular sub-populations, epidermal stem cell content and wound response-enhancing biological properties.

Introduction: Autologous cell suspension (ACS)-based therapy represents a highly promising approach for burns and chronic wounds. However, existing technologies have not achieved the desired clinical success due to several limitations. To overcome practical and cost-associated obstacles of existing ACS methods, we have established a novel methodology for rapid, enzymatic disaggregation of human skin cells and their isolation using a procedure that requires no specialist laboratory instrumentation and is performed at room temperature. Methods: Cells were isolated using enzymatic disaggregation of split-thickness human skin followed by several filtration steps for isolation of cell populations, and cell viability was determined. Individual population recovery was confirmed in appropriate culture medium types, and the presence of epidermal stem cells (EpSCs) within keratinocyte sub-populations was defined by flow cytometry via detection of CD49 and CD71. Positive mediators of wound healing secreted by ACS-derived cultures established on a collagen-based wound-bed mimic were detected by proteome arrays and quantified by ELISA, and the role of such mediators was determined by cell proliferation assays. The effect of ACS-derived conditioned-medium on myofibroblasts was investigated using an in-vitro model of myofibroblast differentiation via detection of α-SMA using immunoblotting and immunofluorescence microscopy. Results: Our methodology permitted efficient recovery of keratinocytes, fibroblasts and melanocytes, which remained viable upon long-term culture. ACS-derivatives comprised sub-populations with the CD49-high/CD71-low expression profile known to demarcate EpSCs. Via secretion of mitogenic factors and wound healing-enhancing mediators, the ACS secretome accelerated keratinocyte proliferation and markedly curtailed cytodifferentiation of myofibroblasts, the latter being key mediators of fibrosis and scarring. Discussion: The systematic characterisation of the cell types within our ACS isolates provided evidence for their superior cell viability and the presence of EpSCs that are critical drivers of wound healing. We defined the biological properties of ACS-derived keratinocytes, which include ability to secrete positive mediators of wound healing as well as suppression of myofibroblast cytodifferentiation. Thus, our study provides several lines of evidence that the established ACS isolates comprise highly-viable cell populations which can physically support wound healing and possess biological properties that have the potential to enhance not only the speed but also the quality of wound healing.

A morphofunctional study of the jumping apparatus in globular springtails.

Springtails are notable for their jumping apparatus and latch-mediated spring mechanism. The challenge, in the light of the tiny size and rapid movement of these organisms, has been to understand the morphological intricacies of this spring system. This study takes an approach that integrates SEM, MicroCT, cLSM and high-speed video recordings to understand the composition and functionality of the jumping apparatus in Megalothorax minimus (Neelipleona), Dicyrtomina ornata and Dicyrtomina minuta (Symphypleona). We focus on reconstructing, describing, and understanding the functioning of structures such as basal plates, musculature and furca. The dimensions of the jumping apparatus in Dicyrtomina and Megalothorax differ significantly from those in elongated springtails. A hypothesis of functional coherence between taxa, based on muscle connections and basal plates, is postulated. High-speed video recordings provide information on: 1) furca release timing and function during jumping and self-righting; 2) performance properties of manubrium, dens and mucro in interaction with the ground and in take-off; 3) possible pre-release furca moves. The study underscores the need for further research employing a variety of visualization methods in order to explore additional aspects such as retinaculum unlatching and furca flexion/extension muscles.

Temporal trends and transmission dynamics of pre-treatment HIV-1 drug resistance within and between risk groups in Kenya, 1986-2020.

BACKGROUND: Evidence on the distribution of pre-treatment HIV-1 drug resistance (HIVDR) among risk groups is limited in Africa. We assessed the prevalence, trends and transmission dynamics of pre-treatment HIVDR within and between MSM, people who inject drugs (PWID), female sex workers (FSWs), heterosexuals (HETs) and perinatally infected children in Kenya. METHODS: HIV-1 partial pol sequences from antiretroviral-naive individuals collected from multiple sources between 1986 and 2020 were used. Pre-treatment reverse transcriptase inhibitor (RTI), PI and integrase inhibitor (INSTI) mutations were assessed using the Stanford HIVDR database. Phylogenetic methods were used to determine and date transmission clusters. RESULTS: Of 3567 sequences analysed, 550 (15.4%, 95% CI: 14.2-16.6) had at least one pre-treatment HIVDR mutation, which was most prevalent amongst children (41.3%), followed by PWID (31.0%), MSM (19.9%), FSWs (15.1%) and HETs (13.9%). Overall, pre-treatment HIVDR increased consistently, from 6.9% (before 2005) to 24.2% (2016-20). Among HETs, pre-treatment HIVDR increased from 6.6% (before 2005) to 20.2% (2011-15), but dropped to 6.5% (2016-20). Additionally, 32 clusters with shared pre-treatment HIVDR mutations were identified. The majority of clusters had R0 ≥ 1.0, indicating ongoing transmissions. The largest was a K103N cluster involving 16 MSM sequences sampled between 2010 and 2017, with an estimated time to the most recent common ancestor (tMRCA) of 2005 [95% higher posterior density (HPD), 2000-08], indicating propagation over 12 years. CONCLUSIONS: Compared to HETs, children and key populations had higher levels of pre-treatment HIVDR. Introduction of INSTIs after 2017 may have abrogated the increase in pre-treatment RTI mutations, albeit in the HET population only. Taken together, our findings underscore the need for targeted efforts towards equitable access to ART for children and key populations in Kenya.

Ancient chicken remains reveal the origins of virulence in Marek's disease virus.

The pronounced growth in livestock populations since the 1950s has altered the epidemiological and evolutionary trajectory of their associated pathogens. For example, Marek's disease virus (MDV), which causes lymphoid tumors in chickens, has experienced a marked increase in virulence over the past century. Today, MDV infections kill >90% of unvaccinated birds, and controlling it costs more than US$1 billion annually. By sequencing MDV genomes derived from archeological chickens, we demonstrate that it has been circulating for at least 1000 years. We functionally tested the Meq oncogene, one of 49 viral genes positively selected in modern strains, demonstrating that ancient MDV was likely incapable of driving tumor formation. Our results demonstrate the power of ancient DNA approaches to trace the molecular basis of virulence in economically relevant pathogens.

Comparative effectiveness and duration of protection of ChAdOx1, CoronaVac, BNT162b2, mRNA-1273, and Ad26.COV2.S COVID-19 vaccines for symptomatic and hospitalized Mu, Delta, and Omicron: A test-negative case-control study.

OBJECTIVE: We carried out a study to estimate the vaccine effectiveness (VE) of homologous vaccination schedules against COVID-19, using data from mandatory information systems from Bogota, Colombia. METHODS: A test-negative case-control study in adults from Bogota (Colombia), between March 1st of 2021 and February 25th of 2022. We assess VE among symptomatic COVID-19 cases during the Mul, Delta, and Omicron predominance periods in Bogota, with controls matched by sex, age (±5 years), and date of testing (±7 days), using a case:control ratio of 1:1. We selected homologous vaccination schedules with ChAdOx1, CoronaVac, BNT162b2, mRNA-1273, and Ad26.COV2.S. VE was reported as one minus the odds ratio in adjusted conditional logistic regressions, with their 95% confidence intervals (CI). A p-value < 0.05 was considered statistically significant. RESULTS: 52,913 cases were matched to controls, 16,722 for Mu, 14,094 for Delta, and 22,097 for Omicron. VE was high against COVID-19 during Mu weeks with full vaccination using the monovalent BNT162b2 (VE: 69; 95% CI, 65 to 72) vaccine and ChAdOx1 (VE: 64; 95% CI, 31 to 81) and significantly lower with CoronaVac (P < 0.001) and Ad26.COV2.S (P = 0.005). During Delta, VE against COVID-19 was higher with BNT162b2 (VE: 55; 95% CI, 51 to 58). The VE for COVID-19 cases during Omicron was higher with a booster dose of monovalent BNT162b2 (VE: 45; 95% CI, 34 to 54). The VE of primary series and booster for ChAdOx1, Ad26.COV2.S, and CoronaVac did not show protection for Omicron. CONCLUSION: Our study provides further evidence on the protective effect of mRNA vaccines for Omicron, and warrant that the duration of protection against symptomatic infection may last for only a few months.

Advances in Animal Models and Cutting-Edge Research in Alternatives: Proceedings of the Third International Conference on 3Rs Research and Progress, Vishakhapatnam, 2022.

Animal experimentation has been integral to drug discovery and development and safety assessment for many years, since it provides insights into the mechanisms of drug efficacy and toxicity (e.g. pharmacology, pharmacokinetics and pharmacodynamics). However, due to species differences in physiology, metabolism and sensitivity to drugs, the animal models can often fail to replicate the effects of drugs and chemicals in human patients, workers and consumers. Researchers across the globe are increasingly applying the Three Rs principles by employing innovative methods in research and testing. The Three Rs concept focuses on: the replacement of animal models (e.g. with in vitro and in silico models or human studies), on the reduction of the number of animals required to achieve research objectives, and on the refinement of existing experimental practices (e.g. eliminating distress and enhancing animal wellbeing). For the last two years, Oncoseek Bio-Acasta Health, a 3-D cell culture-based cutting-edge translational biotechnology company, has organised an annual International Conference on 3Rs Research and Progress. This series of global conferences aims to bring together researchers with diverse expertise and interests, and provides a platform where they can share and discuss their research to promote practices according to the Three Rs principles. In November 2022, the 3rd international conference, Advances in Animal Models and Cutting-Edge Research in Alternatives, took place at the GITAM University in Vishakhapatnam (AP, India) in a hybrid format (i.e. online and in-person). These conference proceedings provide details of the presentations, which were categorised under five different topic sessions. It also describes a special interactive session on in silico strategies for preclinical research in oncology, which was held at the end of the first day.

Understanding the evolution of immune genes in jawed vertebrates.

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

Reductions in Childhood Pneumonia Mortality After Vaccination in the United States.

BACKGROUND: We aim to estimate the magnitude of the reduction in pneumococcal pneumonia and meningitis mortality after the mass introduction of pneumococcal conjugate vaccine (PCV)7 and PCV13 in children in the United States. METHODS: We assessed the trends in mortality rates from pneumococcal pneumonia and meningitis, in the United States between 1994 and 2017. We fitted an interrupted time-series negative binomial regression model (adjusted by trend, seasonality, PCV7/PCV13 coverage, and H. influenzae type b vaccine coverage) to estimate the counterfactual rates without vaccination. We reported a percent reduction in mortality estimates relative to the projected no-vaccination scenario, using the formula 1 minus the incidence risk ratio, with 95% confidence intervals (CIs). RESULTS: Between 1994 and 1999 (the prevaccination period), the all-cause pneumonia mortality rate for 0-1-month-old children was 2.55 per 100,00 pop., whereas for 2-11 months-old children, this rate was 0.82 deaths per 100,000 pop. During the PCV7-period in 0-59-month-old children in the United States, the adjusted reduction of all-cause pneumonia was 13% (95% CI: 4-21) and 19% (95% CI: 0-33) of all-cause meningitis For PCV13, the reductions in this age group were 21% (95% CI: 4-35) for all-cause pneumonia mortality and 22% (95% CI: -19 to 48) for all-cause meningitis mortality. PCV13 had greater reductions of all-cause pneumonia than PCV13 in 6-11-month-old infants. CONCLUSIONS: The universal introduction of PCV7, and later PCV13, for children 0-59 months old in the United States was associated with decreases in mortality due to all-cause pneumonia.

Diagnostic performance of two rapid tests for syphilis screening in people living with HIV in Cali, Colombia.

INTRODUCTION: There is insufficient evidence supporting the use of rapid diagnostic tests (RDTs) for syphilis in people living with HIV (PLWH). We evaluated the diagnostic performance of two commercially available RDTs (Bioline and Determine) in PLWH in Cali, Colombia. METHODS: A cross-sectional field validation study on consecutive adults with confirmed HIV diagnosis attending three outpatient clinics. Both RDTs were performed on capillary blood (CB), obtained by finger prick, and sera, by venipuncture. A combination of treponemal enzyme linked immunosorbent assay (ELISA) and Treponema pallidum haemagglutination assay (TPHA) on serum samples was the reference standard. Rapid plasma reagin (RPR) and clinical criteria were added to define active syphilis. Sensitivity and specificity, predictive values and likelihood ratios (LR) of RDTs were estimated with their corresponding 95% confidence interval (95% CI). Stratified analyses by sample type, patient characteristics, non-treponemal titers, operator and re-training were performed. RESULTS: 244 PLWH were enrolled, of whom 112 (46%) had positive treponemal reference tests and 26/234 (11.1%) had active syphilis. The sensitivities of Bioline on CB and sera were similar (96.4% vs 94.6%, p = 0.6). In contrast, Determine had a lower sensitivity on CB than sera (87.5% vs 99.1%, p<0.001). Sensitivities were lower in PLWH not receiving ART (Bioline 87.1% and Determine 64.5%, p<0.001) and for one of the operators (Bioline 85% and Determine 60%, p<0.001). Specificities of the RDTs were > 95% in most analyses. Predictive values were 90% or higher. For active syphilis, the RDTs showed a similar performance pattern but with decreased specificities. CONCLUSION: The studied RDTs have an excellent performance in PLWH to screen for syphilis and potentially for active syphilis, yet Determine performs better on sera than CB. Patient characteristics and potential difficulties operators may face in acquiring enough blood volume from finger pricks should be considered for the implementation and the interpretation of RDTs.

Norecopa: A global knowledge base of resources for improving animal research and testing.

There are good ethical, legal and scientific reasons for ensuring that our use of animals in research and testing is limited to the lowest number of animals, and that those which are used are treated as humanely as possible, while at the same time providing reliable, reproducible and translatable data which is adequately reported. Unfortunately, there is widespread evidence that there is room for improvement in all these areas. This paper describes the Norecopa website, which offers links to global resources which can be used to resolve these issues. Much of the website content is linked to the PREPARE guidelines for planning any research or testing which appears to need animals. Attention to detail on all steps of the pathway from early planning to manuscript submission should lead to better science, improved animal welfare, and fewer health and safety accidents. This will also minimize the chances of manuscript rejection due to inadequate planning, avoiding a waste of human resources and animal lives.

Life-threatening arrhythmogenic CaM mutations disrupt CaM binding to a distinct RyR2 CaM-binding pocket.

Calmodulin (CaM) modulates the activity of several proteins that play a key role in excitation-contraction coupling (ECC). In cardiac muscle, the major binding partner of CaM is the type-2 ryanodine receptor (RyR2) and altered CaM binding contributes to defects in sarcoplasmic reticulum (SR) calcium (Ca2+) release. Many genetic studies have reported a series of CaM missense mutations in patients with a history of severe arrhythmogenic cardiac disorders. In the present study, we generated four missense CaM mutants (CaMN98I, CaMD132E, CaMD134H and CaMQ136P) and we used a CaM-RyR2 co-immunoprecipitation and a [3H]ryanodine binding assay to directly compare the relative RyR2-binding of wild type and mutant CaM proteins and to investigate the functional effects of these CaM mutations on RyR2 activity. Furthermore, isothermal titration calorimetry (ITC) experiments were performed to investigate and compare the interactions of the wild-type and mutant CaM proteins with various synthetic peptides located in the well-established RyR2 CaM-binding region (3584-3602aa), as well as another CaM-binding region (4255-4271aa) of human RyR2. Our data revealed that all four CaM mutants displayed dramatically reduced RyR2 interaction and defective modulation of [3H]ryanodine binding to RyR2, regardless of LQTS or CPVT association. Moreover, our isothermal titration calorimetry ITC data suggest that RyR2 3584-3602aa and 4255-4271aa regions interact with significant affinity with wild-type CaM, in the presence and absence of Ca2+, two regions that might contribute to a putative intra-subunit CaM-binding pocket. In contrast, screening the interaction of the four arrhythmogenic CaM mutants with two synthetic peptides that correspond to these RyR2 regions, revealed disparate binding properties and signifying differential mechanisms that contribute to reduced RyR2 association.

Impact of host age on viral and bacterial communities in a waterbird population.

Wildlife harbour pathogens that can harm human or livestock health and are the source of most emerging infectious diseases. It is rarely considered how changes in wildlife population age-structures or how age-stratified behaviours might alter the level of pathogen detection within a species, or risk of spillover to other species. Micro-organisms that occur in healthy animals can be an important model for understanding and predicting the dynamics of pathogens of greater health concern, which are hard to study in wild populations due to their relative rarity. We therefore used a metagenomic approach to jointly characterise viral and prokaryotic carriage in faeces collected from a healthy wild bird population (Cygnus olor; mute swan) that has been subject to long-term study. Using 223 samples from known individuals allowed us to compare differences in prokaryotic and eukaryotic viral carriage between adults and juveniles at an unprecedented level of detail. We discovered and characterised 77 novel virus species, of which 21% belong putatively to bird-infecting families, and described the core prokaryotic microbiome of C. olor. Whilst no difference in microbiota diversity was observed between juveniles and adult individuals, 50% (4/8) of bird-infecting virus families (picornaviruses, astroviruses, adenoviruses and bornaviruses) and 3.4% (9/267) of prokaryotic families (including Helicobacteraceae, Spirochaetaceae and Flavobacteriaceae families) were differentially abundant and/or prevalent between juveniles and adults. This indicates that perturbations that affect population age-structures of wildlife could alter circulation dynamics and spillover risk of microbes, potentially including pathogens.

High resolution parallel sequencing reveals multistrain Campylobacter in broiler chicken flocks testing 'negative' by conventional culture methods: implications for control of Campylobacter infection.

Contaminated chicken meat is a major source of human Campylobacteriosis and rates of infection remain high, despite efforts to limit the colonisation of broiler (meat) chicken flocks on farms. Using conventional testing methods of culture or qPCR, Campylobacter is typically detected amongst broiler flocks from 3 wk of age, leading to the assumption that infection is introduced horizontally into chicken rearing houses at this time. In this study, we use parallel sequencing of a fragment of the Campylobacter outer membrane protein, encoded by the porA gene, to test for presence of Campylobacter DNA amongst fresh fecal samples collected from broiler flocks aged 23 to 28 d. Campylobacter DNA was detected in all of the 290 samples tested using the porA target, and in 48% of samples using 16S bacterial profiling, irrespective of whether or not Campylobacter could be detected using conventional qPCR thresholds. A single porAf2 variant was predominant among flocks that would be determined to be Campylobacter 'positive' by conventional means, but a diverse pattern was seen among flocks that were Campylobacter 'negative'. The ability to routinely detect low levels of Campylobacter amongst broiler flocks at a much earlier age than would conventionally be identified requires a re-examination of how and when biosecurity measures are best applied for live birds. In addition, it may be useful to investigate why single Campylobacter variants proliferate in some broiler flocks and not others.

Automated quantification of airway wall thickness on chest CT using retina U-Nets - Performance evaluation and application to a large cohort of chest CTs of COPD patients.

PURPOSE: Airway wall thickening is a consequence of chronic inflammatory processes and usually only qualitatively described in CT radiology reports. The purpose of this study is to automatically quantify airway wall thickness in multiple airway generations and assess the diagnostic potential of this parameter in a large cohort of patients with Chronic Obstructive Pulmonary Disease (COPD). MATERIALS AND METHODS: This retrospective, single-center study included a series of unenhanced chest CTs. Inclusion criteria were the mentioning of an explicit COPD GOLD stage in the written radiology report and time period (01/2019-12/2021). A control group included chest CTs with completely unremarkable lungs according to the report. The DICOM images of all cases (axial orientation; slice-thickness: 1 mm; soft-tissue kernel) were processed by an AI algorithm pipeline consisting of (A) a 3D-U-Net for det detection and tracing of the bronchial tree centerlines (B) extraction of image patches perpendicular to the centerlines of the bronchi, and (C) a 2D U-Net for segmentation of airway walls on those patches. The performance of centerline detection and wall segmentation was assessed. The imaging parameter average wall thickness was calculated for bronchus generations 3-8 (AWT3-8) across the lungs. Mean AWT3-8 was compared between five groups (control, COPD Gold I-IV) using non-parametric statistics. Furthermore, the established emphysema score %LAV-950 was calculated and used to classify scans (normal vs. COPD) alone and in combination with AWT3-8. RESULTS: A total of 575 chest CTs were processed. Algorithm performance was very good (airway centerline detection sensitivity: 86.9%; airway wall segmentation Dice score: 0.86). AWT3-8 was statistically significantly greater in COPD patients compared to controls (2.03 vs. 1.87 mm, p < 0.001) and increased with COPD stage. The classifier that combined %LAV-950 and AWT3-8 was superior to the classifier using only %LAV-950 (AUC = 0.92 vs. 0.79). CONCLUSION: Airway wall thickness increases in patients suffering from COPD and is automatically quantifiable. AWT3-8 could become a CT imaging parameter in COPD complementing the established emphysema biomarker %LAV-950. CLINICAL RELEVANCE STATEMENT: Quantitative measurements considering the complete visible bronchial tree instead of qualitative description could enhance radiology reports, allow for precise monitoring of disease progression and diagnosis of early stages of disease.

Effectiveness of CoronaVac and BNT162b2 COVID-19 mass vaccination in Colombia: A population-based cohort study.

Background: In February 2021, Colombia began mass vaccination against COVID-19 using mainly BNT162b2 and CoronaVac vaccines. We aimed to estimate vaccine effectiveness (VE) to prevent COVID-19 symptomatic cases, hospitalization, critical care admission, and deaths in a cohort of 796,072 insured subjects older than 40 years in northern Colombia, a setting with a high SARS-CoV-2 transmission. Methods: We identified individuals vaccinated between March 1st of 2021 and August 15th of 2021. We included symptomatic cases, hospitalizations, critical care admissions, and deaths in patients with confirmed COVID-19 as main outcomes. We calculated VE for each outcome from the hazard ratio in Cox proportionally hazards regressions (adjusted by age, sex, place of residence, diabetes, human immunodeficiency virus, cancer, hypertension, tuberculosis, neurological diseases, and chronic renal disease), with 95% confidence intervals (CI). Findings: A total of 719,735 insured participants of 40 and more years were followed. We found 21,545 laboratory-confirmed symptomatic COVID-19 among unvaccinated population, along with 2874 hospitalizations, 1061 critical care admissions, and 1329 deaths, for a rate of 207.2 per million person-days, 27.1 per million person-days, 10.0 per million person-days, and 12.5 per million person-days, respectively. We found CoronaVac was not effective for any outcome in subjects above 80 years old; but for people 40-79 years of age, we found two doses of CoronaVac reduced hospitalization (33.1%; 95% CI, 14.5-47.7), critical care admission (47.2%; 95% CI, 18.5-65.8), and death (55.7%; 95% CI, 32.5-70.0). We found BNT162b2 was effective for all outcomes in the entire population of subjects above 40 years of age, significantly declining for subjects ≥80 years. Interpretation: Two doses of either CoronaVac in population between 40 and 79 years of age, or BNT162b2 among vaccinated above 40 years old significantly reduced deaths of confirmed COVID-19 in a cohort of individuals from Colombia. Vaccine effectiveness for CoronaVac and BNT162b2 declined with increasing age. Funding: UK National Institute for Health Research, the European Union's Horizon 2020 research and innovation programme, and the Bill & Melinda Gates Foundation.

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.