Vocal learning-associated convergent evolution in mammalian proteins and regulatory elements

Vocal production learning ("vocal learning") is a convergently evolved trait in vertebrates. To identify brain genomic elements associated with mammalian vocal learning, we integrated genomic, anatomical, and neurophysiological data from the Egyptian fruit bat (Rousettus aegyptiacus) with analyses of the genomes of 215 placental mammals. First, we identified a set of proteins evolving more slowly in vocal learners. Then, we discovered a vocal motor cortical region in the Egyptian fruit bat, an emergent vocal learner, and leveraged that knowledge to identify active cis-regulatory elements in the motor cortex of vocal learners. Machine learning methods applied to motor cortex open chromatin revealed 50 enhancers robustly associated with vocal learning whose activity tended to be lower in vocal learners. Our research implicates convergent losses of motor cortex regulatory elements in mammalian vocal learning evolution.

Sensommatic: an efficient pipeline to mine and predict sensory receptor genes in the era of reference-quality genomes

SUMMARY

Sensory receptor gene families have undergone extensive expansion and loss across vertebrate evolution, leading to significant variation in receptor counts between species. However, due to their species-specific nature, conventional reference-based annotation tools often underestimate the true number of sensory receptors in a given species. While there has been an exponential increase in the taxonomic diversity of publicly available genome assemblies in recent years, only ∼30% of vertebrate species on the NCBI database are currently annotated. To overcome these limitations, we developed 'Sensommatic', an automated and accessible sensory receptor annotation pipeline. Sensommatic implements BLAST and AUGUSTUS to mine and predict sensory receptor genes from whole genome assemblies, adopting a one-to-many gene mapping approach. While designed for vertebrates, Sensommatic can be extended to run on non-vertebrate species by generating customized reference files, making it a scalable and generalizable tool.

AVAILABILITY AND IMPLEMENTATION

Source code and associated files are available at: https://github.com/GMHughes/Sensommatic.

Single-haplotype comparative genomics provides insights into lineage-specific structural variation during cat evolution

The role of structurally dynamic genomic regions in speciation is poorly understood due to challenges inherent in diploid genome assembly. Here we reconstructed the evolutionary dynamics of structural variation in five cat species by phasing the genomes of three interspecies F1 hybrids to generate near-gapless single-haplotype assemblies. We discerned that cat genomes have a paucity of segmental duplications relative to great apes, explaining their remarkable karyotypic stability. X chromosomes were hotspots of structural variation, including enrichment with inversions in a large recombination desert with characteristics of a supergene. The X-linked macrosatellite DXZ4 evolves more rapidly than 99.5% of the genome clarifying its role in felid hybrid incompatibility. Resolved sensory gene repertoires revealed functional copy number changes associated with ecomorphological adaptations, sociality and domestication. This study highlights the value of gapless genomes to reveal structural mechanisms underpinning karyotypic evolution, reproductive isolation and ecological niche adaptation.

The evolution of antimicrobial peptides in Chiroptera

High viral tolerance coupled with an extraordinary regulation of the immune response makes bats a great model to study host-pathogen evolution. Although many immune-related gene gains and losses have been previously reported in bats, important gene families such as antimicrobial peptides (AMPs) remain understudied. We built an exhaustive bioinformatic pipeline targeting the major gene families of defensins and cathelicidins to explore AMP diversity and analyze their evolution and distribution across six bat families. A combination of manual and automated procedures identified 29 AMP families across queried species, with α-, β-defensins, and cathelicidins representing around 10% of AMP diversity. Gene duplications were inferred in both α-defensins, which were absent in five species, and three β-defensin gene subfamilies, but cathelicidins did not show significant shifts in gene family size and were absent in Anoura caudifer and the pteropodids. Based on lineage-specific gains and losses, we propose diet and diet-related microbiome evolution may determine the evolution of α- and β-defensins gene families and subfamilies. These results highlight the importance of building species-specific libraries for genome annotation in non-model organisms and shed light on possible drivers responsible for the rapid evolution of AMPs. By focusing on these understudied defenses, we provide a robust framework for explaining bat responses to pathogens.

A novel nematode species from the Siberian permafrost shares adaptive mechanisms for cryptobiotic survival with C. elegans dauer larva

Some organisms in nature have developed the ability to enter a state of suspended metabolism called cryptobiosis when environmental conditions are unfavorable. This state-transition requires execution of a combination of genetic and biochemical pathways that enable the organism to survive for prolonged periods. Recently, nematode individuals have been reanimated from Siberian permafrost after remaining in cryptobiosis. Preliminary analysis indicates that these nematodes belong to the genera Panagrolaimus and Plectus. Here, we present precise radiocarbon dating indicating that the Panagrolaimus individuals have remained in cryptobiosis since the late Pleistocene (~46,000 years). Phylogenetic inference based on our genome assembly and a detailed morphological analysis demonstrate that they belong to an undescribed species, which we named Panagrolaimus kolymaensis. Comparative genome analysis revealed that the molecular toolkit for cryptobiosis in P. kolymaensis and in C. elegans is partly orthologous. We show that biochemical mechanisms employed by these two species to survive desiccation and freezing under laboratory conditions are similar. Our experimental evidence also reveals that C. elegans dauer larvae can remain viable for longer periods in suspended animation than previously reported. Altogether, our findings demonstrate that nematodes evolved mechanisms potentially allowing them to suspend life over geological time scales.

Evolutionary constraint and innovation across hundreds of placental mammals

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Relating enhancer genetic variation across mammals to complex phenotypes using machine learning

Protein-coding differences between species often fail to explain phenotypic diversity, suggesting the involvement of genomic elements that regulate gene expression such as enhancers. Identifying associations between enhancers and phenotypes is challenging because enhancer activity can be tissue-dependent and functionally conserved despite low sequence conservation. We developed the Tissue-Aware Conservation Inference Toolkit (TACIT) to associate candidate enhancers with species' phenotypes using predictions from machine learning models trained on specific tissues. Applying TACIT to associate motor cortex and parvalbumin-positive interneuron enhancers with neurological phenotypes revealed dozens of enhancer-phenotype associations, including brain size-associated enhancers that interact with genes implicated in microcephaly or macrocephaly. TACIT provides a foundation for identifying enhancers associated with the evolution of any convergently evolved phenotype in any large group of species with aligned genomes.

Evolutionary constraint and innovation across hundreds of placental mammals

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Integrating gene annotation with orthology inference at scale

Annotating coding genes and inferring orthologs are two classical challenges in genomics and evolutionary biology that have traditionally been approached separately, limiting scalability. We present TOGA (Tool to infer Orthologs from Genome Alignments), a method that integrates structural gene annotation and orthology inference. TOGA implements a different paradigm to infer orthologous loci, improves ortholog detection and annotation of conserved genes compared with state-of-the-art methods, and handles even highly fragmented assemblies. TOGA scales to hundreds of genomes, which we demonstrate by applying it to 488 placental mammal and 501 bird assemblies, creating the largest comparative gene resources so far. Additionally, TOGA detects gene losses, enables selection screens, and automatically provides a superior measure of mammalian genome quality. TOGA is a powerful and scalable method to annotate and compare genes in the genomic era.

Leveraging base-pair mammalian constraint to understand genetic variation and human disease

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

The contribution of historical processes to contemporary extinction risk in placental mammals

Species persistence can be influenced by the amount, type, and distribution of diversity across the genome, suggesting a potential relationship between historical demography and resilience. In this study, we surveyed genetic variation across single genomes of 240 mammals that compose the Zoonomia alignment to evaluate how historical effective population size (N_e) affects heterozygosity and deleterious genetic load and how these factors may contribute to extinction risk. We find that species with smaller historical N_e carry a proportionally larger burden of deleterious alleles owing to long-term accumulation and fixation of genetic load and have a higher risk of extinction. This suggests that historical demography can inform contemporary resilience. Models that included genomic data were predictive of species' conservation status, suggesting that, in the absence of adequate census or ecological data, genomic information may provide an initial risk assessment.

Mammalian evolution of human cis-regulatory elements and transcription factor binding sites

Understanding the regulatory landscape of the human genome is a long-standing objective of modern biology. Using the reference-free alignment across 241 mammalian genomes produced by the Zoonomia Consortium, we charted evolutionary trajectories for 0.92 million human candidate cis-regulatory elements (cCREs) and 15.6 million human transcription factor binding sites (TFBSs). We identified 439,461 cCREs and 2,024,062 TFBSs under evolutionary constraint. Genes near constrained elements perform fundamental cellular processes, whereas genes near primate-specific elements are involved in environmental interaction, including odor perception and immune response. About 20% of TFBSs are transposable element-derived and exhibit intricate patterns of gains and losses during primate evolution whereas sequence variants associated with complex traits are enriched in constrained TFBSs. Our annotations illuminate the regulatory functions of the human genome.

A genomic timescale for placental mammal evolution

The precise pattern and timing of speciation events that gave rise to all living placental mammals remain controversial. We provide a comprehensive phylogenetic analysis of genetic variation across an alignment of 241 placental mammal genome assemblies, addressing prior concerns regarding limited genomic sampling across species. We compared neutral genome-wide phylogenomic signals using concatenation and coalescent-based approaches, interrogated phylogenetic variation across chromosomes, and analyzed extensive catalogs of structural variants. Interordinal relationships exhibit relatively low rates of phylogenomic conflict across diverse datasets and analytical methods. Conversely, X-chromosome versus autosome conflicts characterize multiple independent clades that radiated during the Cenozoic. Genomic time trees reveal an accumulation of cladogenic events before and immediately after the Cretaceous-Paleogene (K-Pg) boundary, implying important roles for Cretaceous continental vicariance and the K-Pg extinction in the placental radiation.

The functional and evolutionary impacts of human-specific deletions in conserved elements

Conserved genomic sequences disrupted in humans may underlie uniquely human phenotypic traits. We identified and characterized 10,032 human-specific conserved deletions (hCONDELs). These short (average 2.56 base pairs) deletions are enriched for human brain functions across genetic, epigenomic, and transcriptomic datasets. Using massively parallel reporter assays in six cell types, we discovered 800 hCONDELs conferring significant differences in regulatory activity, half of which enhance rather than disrupt regulatory function. We highlight several hCONDELs with putative human-specific effects on brain development, including HDAC5, CPEB4, and PPP2CA. Reverting an hCONDEL to the ancestral sequence alters the expression of LOXL2 and developmental genes involved in myelination and synaptic function. Our data provide a rich resource to investigate the evolutionary mechanisms driving new traits in humans and other species.

Leveraging Base Pair Mammalian Constraint to Understand Genetic Variation and Human Disease

Although thousands of genomic regions have been associated with heritable human diseases, attempts to elucidate biological mechanisms are impeded by a general inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function that is agnostic to cell type or disease mechanism. Here, single base phyloP scores from the whole genome alignment of 240 placental mammals identified 3.5% of the human genome as significantly constrained, and likely functional. We compared these scores to large-scale genome annotation, genome-wide association studies (GWAS), copy number variation, clinical genetics findings, and cancer data sets. Evolutionarily constrained positions are enriched for variants explaining common disease heritability (more than any other functional annotation). Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Bat pluripotent stem cells reveal unusual entanglement between host and viruses

Bats are distinctive among mammals due to their ability to fly, use laryngeal echolocation, and tolerate viruses. However, there are currently no reliable cellular models for studying bat biology or their response to viral infections. Here, we created induced pluripotent stem cells (iPSCs) from two species of bats: the wild greater horseshoe bat (Rhinolophus ferrumequinum) and the greater mouse-eared bat (Myotis myotis). The iPSCs from both bat species showed similar characteristics and had a gene expression profile resembling that of cells attacked by viruses. They also had a high number of endogenous viral sequences, particularly retroviruses. These results suggest that bats have evolved mechanisms to tolerate a large load of viral sequences and may have a more intertwined relationship with viruses than previously thought. Further study of bat iPSCs and their differentiated progeny will provide insights into bat biology, virus host relationships, and the molecular basis of bats' special traits.

The pale spear-nosed bat: A neuromolecular and transgenic model for vocal learning

Vocal learning, the ability to produce modified vocalizations via learning from acoustic signals, is a key trait in the evolution of speech. While extensively studied in songbirds, mammalian models for vocal learning are rare. Bats present a promising study system given their gregarious natures, small size, and the ability of some species to be maintained in captive colonies. We utilize the pale spear-nosed bat (Phyllostomus discolor) and report advances in establishing this species as a tractable model for understanding vocal learning. We have taken an interdisciplinary approach, aiming to provide an integrated understanding across genomics (Part I), neurobiology (Part II), and transgenics (Part III). In Part I, we generated new, high-quality genome annotations of coding genes and noncoding microRNAs to facilitate functional and evolutionary studies. In Part II, we traced connections between auditory-related brain regions and reported neuroimaging to explore the structure of the brain and gene expression patterns to highlight brain regions. In Part III, we created the first successful transgenic bats by manipulating the expression of FoxP2, a speech-related gene. These interdisciplinary approaches are facilitating a mechanistic and evolutionary understanding of mammalian vocal learning and can also contribute to other areas of investigation that utilize P. discolor or bats as study species.

Genome analysis of the metabolically versatile Pseudomonas umsongensis GO16: the genetic basis for PET monomer upcycling into polyhydroxyalkanoates

The throwaway culture related to the single-use materials such as polyethylene terephthalate (PET) has created a major environmental concern. Recycling of PET waste into biodegradable plastic polyhydroxyalkanoate (PHA) creates an opportunity to improve resource efficiency and contribute to a circular economy. We sequenced the genome of Pseudomonas umsongensis GO16 previously shown to convert PET-derived terephthalic acid (TA) into PHA and performed an in-depth genome analysis. GO16 can degrade a range of aromatic substrates in addition to TA, due to the presence of a catabolic plasmid pENK22. The genetic complement required for the degradation of TA via protocatechuate was identified and its functionality was confirmed by transferring the tph operon into Pseudomonas putida KT2440, which is unable to utilize TA naturally. We also identified the genes involved in ethylene glycol (EG) metabolism, the second PET monomer, and validated the capacity of GO16 to use EG as a sole source of carbon and energy. Moreover, GO16 possesses genes for the synthesis of both medium and short chain length PHA and we have demonstrated the capacity of the strain to convert mixed TA and EG into PHA. The metabolic versatility of GO16 highlights the potential of this organism for biotransformations using PET waste as a feedstock.

Large-scale genome sampling reveals unique immunity and metabolic adaptations in bats

Comprising more than 1,400 species, bats possess adaptations unique among mammals including powered flight, unexpected longevity, and extraordinary immunity. Some of the molecular mechanisms underlying these unique adaptations includes DNA repair, metabolism and immunity. However, analyses have been limited to a few divergent lineages, reducing the scope of inferences on gene family evolution across the Order Chiroptera. We conducted an exhaustive comparative genomic study of 37 bat species, one generated in this study, encompassing a large number of lineages, with a particular emphasis on multi-gene family evolution across immune and metabolic genes. In agreement with previous analyses, we found lineage-specific expansions of the APOBEC3 and MHC-I gene families, and loss of the proinflammatory PYHIN gene family. We inferred more than 1,000 gene losses unique to bats, including genes involved in the regulation of inflammasome pathways such as epithelial defence receptors, the natural killer gene complex and the interferon-gamma induced pathway. Gene set enrichment analyses revealed genes lost in bats are involved in defence response against pathogen-associated molecular patterns and damage-associated molecular patterns. Gene family evolution and selection analyses indicate bats have evolved fundamental functional differences compared to other mammals in both innate and adaptive immune system, with the potential to enhance antiviral immune response while dampening inflammatory signalling. In addition, metabolic genes have experienced repeated expansions related to convergent shifts to plant-based diets. Our analyses support the hypothesis that, in tandem with flight, ancestral bats had evolved a unique set of immune adaptations whose functional implications remain to be explored.

Evolution of mammalian longevity: age-related increase in autophagy in bats compared to other mammals

Autophagy maintains cellular homeostasis and its dysfunction has been implicated in aging. Bats are the longest-lived mammals for their size, but the molecular mechanisms underlying their extended healthspan are not well understood. Here, drawing on >8 years of mark-recapture field studies, we report the first longitudinal analysis of autophagy regulation in bats. Mining of published population level aging blood transcriptomes (M. myotis, mouse and human) highlighted a unique increase of autophagy related transcripts with age in bats, but not in other mammals. This bat-specific increase in autophagy transcripts was recapitulated by the western blot determination of the autophagy marker, LC3II/I ratio, in skin primary fibroblasts (Myotis myotis,Pipistrellus kuhlii, mouse), that also showed an increase with age in both bat species. Further phylogenomic selection pressure analyses across eutherian mammals (n=70 taxa; 274 genes) uncovered 10 autophagy-associated genes under selective pressure in bat lineages. These molecular adaptations potentially mediate the exceptional age-related increase of autophagy signalling in bats, which may contribute to their longer healthspans.

Broad host range of SARS-CoV-2 predicted by comparative and structural analysis of ACE2 in vertebrates

The novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of COVID-19. The main receptor of SARS-CoV-2, angiotensin I converting enzyme 2 (ACE2), is now undergoing extensive scrutiny to understand the routes of transmission and sensitivity in different species. Here, we utilized a unique dataset of ACE2 sequences from 410 vertebrate species, including 252 mammals, to study the conservation of ACE2 and its potential to be used as a receptor by SARS-CoV-2. We designed a five-category binding score based on the conservation properties of 25 amino acids important for the binding between ACE2 and the SARS-CoV-2 spike protein. Only mammals fell into the medium to very high categories and only catarrhine primates into the very high category, suggesting that they are at high risk for SARS-CoV-2 infection. We employed a protein structural analysis to qualitatively assess whether amino acid changes at variable residues would be likely to disrupt ACE2/SARS-CoV-2 spike protein binding and found the number of predicted unfavorable changes significantly correlated with the binding score. Extending this analysis to human population data, we found only rare (frequency <0.001) variants in 10/25 binding sites. In addition, we found significant signals of selection and accelerated evolution in the ACE2 coding sequence across all mammals, and specific to the bat lineage. Our results, if confirmed by additional experimental data, may lead to the identification of intermediate host species for SARS-CoV-2, guide the selection of animal models of COVID-19, and assist the conservation of animals both in native habitats and in human care.

Six reference-quality genomes reveal evolution of bat adaptations

Bats possess extraordinary adaptations, including flight, echolocation, extreme longevity and unique immunity. High-quality genomes are crucial for understanding the molecular basis and evolution of these traits. Here we incorporated long-read sequencing and state-of-the-art scaffolding protocols1 to generate, to our knowledge, the first reference-quality genomes of six bat species (Rhinolophus ferrumequinum, Rousettus aegyptiacus, Phyllostomus discolor, Myotis myotis, Pipistrellus kuhlii and Molossus molossus). We integrated gene projections from our 'Tool to infer Orthologs from Genome Alignments' (TOGA) software with de novo and homology gene predictions as well as short- and long-read transcriptomics to generate highly complete gene annotations. To resolve the phylogenetic position of bats within Laurasiatheria, we applied several phylogenetic methods to comprehensive sets of orthologous protein-coding and noncoding regions of the genome, and identified a basal origin for bats within Scrotifera. Our genome-wide screens revealed positive selection on hearing-related genes in the ancestral branch of bats, which is indicative of laryngeal echolocation being an ancestral trait in this clade. We found selection and loss of immunity-related genes (including pro-inflammatory NF-κB regulators) and expansions of anti-viral APOBEC3 genes, which highlights molecular mechanisms that may contribute to the exceptional immunity of bats. Genomic integrations of diverse viruses provide a genomic record of historical tolerance to viral infection in bats. Finally, we found and experimentally validated bat-specific variation in microRNAs, which may regulate bat-specific gene-expression programs. Our reference-quality bat genomes provide the resources required to uncover and validate the genomic basis of adaptations of bats, and stimulate new avenues of research that are directly relevant to human health and disease1.

Broad Host Range of SARS-CoV-2 Predicted by Comparative and Structural Analysis of ACE2 in Vertebrates

The novel coronavirus SARS-CoV-2 is the cause of Coronavirus Disease-2019 (COVID-19). The main receptor of SARS-CoV-2, angiotensin I converting enzyme 2 (ACE2), is now undergoing extensive scrutiny to understand the routes of transmission and sensitivity in different species. Here, we utilized a unique dataset of 410 vertebrates, including 252 mammals, to study cross-species conservation of ACE2 and its likelihood to function as a SARS-CoV-2 receptor. We designed a five-category ranking score based on the conservation properties of 25 amino acids important for the binding between receptor and virus, classifying all species from very high to very low. Only mammals fell into the medium to very high categories, and only catarrhine primates in the very high category, suggesting that they are at high risk for SARS-CoV-2 infection. We employed a protein structural analysis to qualitatively assess whether amino acid changes at variable residues would be likely to disrupt ACE2/SARS-CoV-2 binding, and found the number of predicted unfavorable changes significantly correlated with the binding score. Extending this analysis to human population data, we found only rare (<0.1%) variants in 10/25 binding sites. In addition, we observed evidence of positive selection in ACE2 in multiple species, including bats. Utilized appropriately, our results may lead to the identification of intermediate host species for SARS-CoV-2, justify the selection of animal models of COVID-19, and assist the conservation of animals both in native habitats and in human care.

Olfactory receptor repertoire size in dinosaurs

The olfactory bulb (OB) ratio is the size of the OB relative to the cerebral hemisphere, and is used to estimate the proportion of the forebrain devoted to smell. In birds, OB ratio correlates with the number of olfactory receptor (OR) genes and therefore has been used as a proxy for olfactory acuity. By coupling OB ratios with known OR gene repertoires in birds, we infer minimum repertoire sizes for extinct taxa, including non-avian dinosaurs, using phylogenetic modelling, ancestral state reconstruction and comparative genomics. We highlight a shift in the scaling of OB ratio to body size along the lineage leading to modern birds, demonstrating variable OR repertoires present in different dinosaur and crown-bird lineages, with varying factors potentially influencing sensory evolution in theropods. We investigate the ancestral sensory space available to extinct taxa, highlighting potential adaptations to ecological niches. Through combining morphological and genomic data, we show that, while genetic information for extinct taxa is forever lost, it is potentially feasible to investigate evolutionary trajectories in extinct genomes.

As Blind as a Bat? Opsin Phylogenetics Illuminates the Evolution of Color Vision in Bats

Through their unique use of sophisticated laryngeal echolocation bats are considered sensory specialists amongst mammals and represent an excellent model in which to explore sensory perception. Although several studies have shown that the evolution of vision is linked to ecological niche adaptation in other mammalian lineages, this has not yet been fully explored in bats. Recent molecular analysis of the opsin genes, which encode the photosensitive pigments underpinning color vision, have implicated high-duty cycle (HDC) echolocation and the adoption of cave roosting habits in the degeneration of color vision in bats. However, insufficient sampling of relevant taxa has hindered definitive testing of these hypotheses. To address this, novel sequence data was generated for the SWS1 and MWS/LWS opsin genes and combined with existing data to comprehensively sample species representing diverse echolocation types and niches (SWS1 n = 115; MWS/LWS n = 45). A combination of phylogenetic analysis, ancestral state reconstruction, and selective pressure analyses were used to reconstruct the evolution of these visual pigments in bats and revealed that although both genes are evolving under purifying selection in bats, MWS/LWS is highly conserved but SWS1 is highly variable. Spectral tuning analyses revealed that MWS/LWS opsin is tuned to a long wavelength, 555-560 nm in the bat ancestor and the majority of extant taxa. The presence of UV vision in bats is supported by our spectral tuning analysis, but phylogenetic analyses demonstrated that the SWS1 opsin gene has undergone pseudogenization in several lineages. We do not find support for a link between the evolution of HDC echolocation and the pseudogenization of the SWS1 gene in bats, instead we show the SWS1 opsin is functional in the HDC echolocator, Pteronotus parnellii. Pseudogenization of the SWS1 is correlated with cave roosting habits in the majority of pteropodid species. Together these results demonstrate that the loss of UV vision in bats is more widespread than was previously considered and further elucidate the role of ecological niche specialization in the evolution of vision in bats.

The Birth and Death of Olfactory Receptor Gene Families in Mammalian Niche Adaptation

The olfactory receptor (OR) gene families, which govern mammalian olfaction, have undergone extensive expansion and contraction through duplication and pseudogenization. Previous studies have shown that broadly defined environmental adaptations (e.g., terrestrial vs. aquatic) are correlated with the number of functional and non-functional OR genes retained. However, to date, no study has examined species-specific gene duplications in multiple phylogenetically divergent mammals to elucidate OR evolution and adaptation. Here, we identify the OR gene families driving adaptation to different ecological niches by mapping the fate of species-specific gene duplications in the OR repertoire of 94 diverse mammalian taxa, using molecular phylogenomic methods. We analyze >70,000 OR gene sequences mined from whole genomes, generated from novel amplicon sequencing data, and collated with data from previous studies, comprising one of the largest OR studies to date. For the first time, we demonstrate statistically significant patterns of OR species-specific gene duplications associated with the presence of a functioning vomeronasal organ. With respect to dietary niche, we uncover a novel link between a large number of duplications in OR family 5/8/9 and herbivory. Our results also highlight differences between social and solitary niches, indicating that a greater OR repertoire expansion may be associated with a solitary lifestyle. This study demonstrates the utility of species-specific duplications in elucidating gene family evolution, revealing how the OR repertoire has undergone expansion and contraction with respect to a number of ecological adaptations in mammals.

Growing old, yet staying young: The role of telomeres in bats' exceptional longevity

Understanding aging is a grand challenge in biology. Exceptionally long-lived animals have mechanisms that underpin extreme longevity. Telomeres are protective nucleotide repeats on chromosome tips that shorten with cell division, potentially limiting life span. Bats are the longest-lived mammals for their size, but it is unknown whether their telomeres shorten. Using >60 years of cumulative mark-recapture field data, we show that telomeres shorten with age in Rhinolophus ferrumequinum and Miniopterus schreibersii, but not in the bat genus with greatest longevity, Myotis. As in humans, telomerase is not expressed in Myotis myotis blood or fibroblasts. Selection tests on telomere maintenance genes show that ATM and SETX, which repair and prevent DNA damage, potentially mediate telomere dynamics in Myotis bats. Twenty-one telomere maintenance genes are differentially expressed in Myotis, of which 14 are enriched for DNA repair, and 5 for alternative telomere-lengthening mechanisms. We demonstrate how telomeres, telomerase, and DNA repair genes have contributed to the evolution of exceptional longevity in Myotis bats, advancing our understanding of healthy aging.

Is there a link between aging and microbiome diversity in exceptional mammalian longevity?

A changing microbiome has been linked to biological aging in mice and humans, suggesting a possible role of gut flora in pathogenic aging phenotypes. Many bat species have exceptional longevity given their body size and some can live up to ten times longer than expected with little signs of aging. This study explores the anal microbiome of the exceptionally long-lived Myotis myotis bat, investigating bacterial composition in both adult and juvenile bats to determine if the microbiome changes with age in a wild, long-lived non-model organism, using non-lethal sampling. The anal microbiome was sequenced using metabarcoding in more than 50 individuals, finding no significant difference between the composition of juvenile and adult bats, suggesting that age-related microbial shifts previously observed in other mammals may not be present in Myotis myotis. Functional gene categories, inferred from metabarcoding data, expressed in the M. myotis microbiome were categorized identifying pathways involved in metabolism, DNA repair and oxidative phosphorylation. We highlight an abundance of ‘Proteobacteria’ relative to other mammals, with similar patterns compared to other bat microbiomes. Our results suggest that M. myotis may have a relatively stable, unchanging microbiome playing a role in their extended ‘health spans’ with the advancement of age, and suggest a potential link between microbiome and sustained, powered flight.

Genome-wide signatures of complex introgression and adaptive evolution in the big cats

The great cats of the genus Panthera comprise a recent radiation whose evolutionary history is poorly understood. Their rapid diversification poses challenges to resolving their phylogeny while offering opportunities to investigate the historical dynamics of adaptive divergence. We report the sequence, de novo assembly, and annotation of the jaguar (Panthera onca) genome, a novel genome sequence for the leopard (Panthera pardus), and comparative analyses encompassing all living Panthera species. Demographic reconstructions indicated that all of these species have experienced variable episodes of population decline during the Pleistocene, ultimately leading to small effective sizes in present-day genomes. We observed pervasive genealogical discordance across Panthera genomes, caused by both incomplete lineage sorting and complex patterns of historical interspecific hybridization. We identified multiple signatures of species-specific positive selection, affecting genes involved in craniofacial and limb development, protein metabolism, hypoxia, reproduction, pigmentation, and sensory perception. There was remarkable concordance in pathways enriched in genomic segments implicated in interspecies introgression and in positive selection, suggesting that these processes were connected. We tested this hypothesis by developing exome capture probes targeting ~19,000 Panthera genes and applying them to 30 wild-caught jaguars. We found at least two genes (DOCK3 and COL4A5, both related to optic nerve development) bearing significant signatures of interspecies introgression and within-species positive selection. These findings indicate that post-speciation admixture has contributed genetic material that facilitated the adaptive evolution of big cat lineages.

Nanoparticle Trapping and Characterization Using Open Microcavities

Characterization and trapping of nanoparticles in solution is of great importance for lab-on-a-chip applications in biomedical, environmental, and materials sciences. Devices are now starting to emerge allowing such manipulations and investigations in real-time. Better insights into the interaction between the nanoparticle and the optical trap is therefore necessary in order to move forward in this field. In this work, we present a new kind of nanotweezers based on open microcavities. We show that by monitoring the cavity mode wavelength shift as the particle diffuses through the cavity, it is possible to establish both the nanoparticle polarizability and its coefficient of friction. Additionally, our experiment provides a deep insight in the interaction between the nanoparticle and the cavity mode. The technique has built-in calibration of the trap strength and spring constant, making it attractive for practical applications. This work illustrates the potential of such optical microcavities for future developments in nanoparticle sensors and lab-on-a-chip devices.

A novel method of microsatellite genotyping-by-sequencing using individual combinatorial barcoding

This study examines the potential of next-generation sequencing based 'genotyping-by-sequencing' (GBS) of microsatellite loci for rapid and cost-effective genotyping in large-scale population genetic studies. The recovery of individual genotypes from large sequence pools was achieved by PCR-incorporated combinatorial barcoding using universal primers. Three experimental conditions were employed to explore the possibility of using this approach with existing and novel multiplex marker panels and weighted amplicon mixture. The GBS approach was validated against microsatellite data generated by capillary electrophoresis. GBS allows access to the underlying nucleotide sequences that can reveal homoplasy, even in large datasets and facilitates cross laboratory transfer. GBS of microsatellites, using individual combinatorial barcoding, is potentially faster and cheaper than current microsatellite approaches and offers better and more data.

Open-access optical microcavities for lab-on-a-chip refractive index sensing

Open-access optical microcavities provide a novel approach to label-free lab-on-a-chip optofluidic sensing. They offer direct access to a highly confined electromagnetic field, and yield a femtoliter detection volume. This article describes the characteristics of these devices for refractive index sensing. We show that most of the ambient noise can be removed from the refractive index data by simultaneous tracking of resonances across an array of cavities. A sensitivity of 3.5 × 10(-4) RIU is demonstrated which corresponds to detecting the refractive index change caused by the presence of 500,000 glucose molecules in aqueous solution.

Stable isotope imaging of biological samples with high resolution secondary ion mass spectrometry and complementary techniques

Stable isotopes are ideal labels for studying biological processes because they have little or no effect on the biochemical properties of target molecules. The NanoSIMS is a tool that can image the distribution of stable isotope labels with up to 50 nm spatial resolution and with good quantitation. This combination of features has enabled several groups to undertake significant experiments on biological problems in the last decade. Combining the NanoSIMS with other imaging techniques also enables us to obtain not only chemical information but also the structural information needed to understand biological processes. This article describes the methodologies that we have developed to correlate atomic force microscopy and backscattered electron imaging with NanoSIMS experiments to illustrate the imaging of stable isotopes at molecular, cellular, and tissue scales. Our studies make it possible to address 3 biological problems: (1) the interaction of antimicrobial peptides with membranes; (2) glutamine metabolism in cancer cells; and (3) lipoprotein interactions in different tissues.

Loss of olfactory receptor function in hominin evolution

The mammalian sense of smell is governed by the largest gene family, which encodes the olfactory receptors (ORs). The gain and loss of OR genes is typically correlated with adaptations to various ecological niches. Modern humans have 853 OR genes but 55% of these have lost their function. Here we show evidence of additional OR loss of function in the Neanderthal and Denisovan hominin genomes using comparative genomic methodologies. Ten Neanderthal and 8 Denisovan ORs show evidence of loss of function that differ from the reference modern human OR genome. Some of these losses are also present in a subset of modern humans, while some are unique to each lineage. Morphological changes in the cranium of Neanderthals suggest different sensory arrangements to that of modern humans. We identify differences in functional olfactory receptor genes among modern humans, Neanderthals and Denisovans, suggesting varied loss of function across all three taxa and we highlight the utility of using genomic information to elucidate the sensory niches of extinct species.

Using Illumina next generation sequencing technologies to sequence multigene families in de novo species

The advent of Next Generation Sequencing Technology (NGST) has revolutionized molecular biology research, allowing for rapid gene/genome sequencing from a multitude of diverse species. As high throughput sequencing becomes more accessible, more efficient workflows must be developed to deal with the amounts of data produced and better assemble the genomes of de novo lineages. We combine traditional laboratory methods with Illumina NGST to amplify and sequence the largest mammalian multigene family, the Olfactory Receptor gene family, for species with and without a reference genome. We develop novel assembly methods to annotate and filter these data, which can be utilized for any gene family or any species. We find no significant difference between the ratio of genes within their respective gene families of our data compared with available genomic data. Using simulated data we explore the limitations of short-read sequence data and our assembly in recovering this gene family. We highlight the benefits and shortcomings of these methods. Compared with data generated from traditional polymerase chain reaction, cloning and Sanger sequencing methodologies, sequence data generated using our pipeline increases yield and sequencing efficiency without reducing the number of unique genes amplified. A cloning step is not required, therefore shortening data generation time. The novel downstream methodologies and workflows described provide a tool to be utilized by many fields of biology, to access and analyze the vast quantities of data generated. By combining laboratory and in silico methods, we provide a means of extracting genomic information for multigene families without complete genome sequencing.

Acute effects of a herb extract formulation and inulin fibre on appetite, energy intake and food choice

The impact of two commercially available products, a patented herb extract Yerbe Maté, Guarana and Damiana (YGD) formulation and an inulin-based soluble fermentable fibre (SFF), alone or in combination, on appetite and food intake were studied for the first time in a double blind, placebo-controlled, cross-over design. 58 normal to slightly overweight women consumed a fixed-load breakfast followed 4h later by an ad libitum lunch. They were administered YGD (3 tablets) and SFF (5g in 100ml water), YGD and water (100ml), SFF and placebo (3 tablets) or water and placebo 15min before meals. Appetite was assessed using visual analogue scales, and energy intake was measured at lunch. Significant reductions in food intake and energy intake were observed when YGD was present (59.5g, 16.3%; 112.4kcal, 17.3%) and when SFF was present (31.9g, 9.1%; 80kcal, 11.7%) compared with conditions were products were absent. The lowest intake (gram and kcal) was in the YGD+SFF condition. Significant reductions in AUC hunger and AUC desire to eat were also observed after YGD+SFF combination. The data demonstrate that YGD produces a robust short-term effect on caloric intake, an effect augmented by SFF. Caloric compensation for SFF indicates independent effects on appetite regulation.

Intra-operative 3-T MRI for paediatric brain tumours: challenges and perspectives

MRI is the ideal modality for imaging intracranial tumours. Intraoperative MRI (ioMRI) makes it possible to obtain scans during a neurosurgical operation that can aid complete macroscopic tumour resection—a major prognostic factor in the majority of brain tumours in children. Intraoperative MRI can also help limit damage to normal brain tissue. It therefore has the potential to improve the survival of children with brain tumours and to minimise morbidity, including neurological deficits. The use of ioMRI is also likely to reduce the need for second look surgery, and may reduce the need for chemotherapy and radiotherapy. Highfield MRI systems provide better anatomical information and also enable effective utilisation of advanced MRI techniques such as perfusion imaging, diffusion tensor imaging, and magnetic resonance spectroscopy. However, high-field ioMRI facilities require substantial capital investment, and careful planning is required for optimal benefit. Safe ioMRI requires meticulous attention to detail and rigorous application of magnetic field safety precautions. Interpretation of ioMRI can be challenging and requires experience and understanding of artefacts that are common in the intra-operative setting.

Preparation of location-specific thin foils from Fe-3% Si bi- and tri-crystals for examination in a FEG-STEM

Bi-crystals and tri-crystals of a nominal Fe-3% Si (wt%) of well-defined orientations have been grown using a floating-zone technique with optical heating. The manufacture of these unique crystals and the preparation technique involved in harvesting thin foils from specific locations for transmission electron microscopy are described in detail. In particular, the grain boundary triple junction has been extracted from the tri-crystal and examined in high-resolution aberration-corrected FEG-STEM instruments. To achieve the necessary resolution, the foils have to be uniformly thin, in the range 50-100 nm over large areas of the specimen. For ferromagnetic materials, there are further challenges arising from the magnetic field interaction, with the electron beam placing significant demands on the aberration correction system. One way to minimise this interaction is to reduce the total mass of magnetic material. To achieve this, an in situ focused ion beam lift-out technique has been combined with an additional precision ion-polishing stage to reproducibly provide thin-foil specimens suitable for high-resolution EELS and EDX analysis. Examination of the foils reveals that the final precision ion-polishing stage removes residual damage arising from the use of focused ion beam milling procedures.

Coping, depressive feelings and gender differences in late life widowhood

The study investigated the relationship between depressive feelings and coping amongst older widowed men and women. Participants were interviewed about their affective experiences of widowhood and completed two depression questionnaire assessments, the Symptoms of Anxiety and Depression Scale (SAD) and the Hospital Anxiety and Depression Scale (HADS). Participants were assessed as either coping or not coping. The results showed that both measures were effective at differentiating those who coped (Copers) from those who did not (Non-Copers) in the sample as a whole. Amongst the widows the HADS significantly differentiated the two groups. Amongst men, neither measure significantly distinguished Copers from Non-Copers. However, an examination of the interviews suggested that widowers reported depressive feelings significantly more often than widows. The results suggest that depressive feelings are associated with non-coping in older widowed people. There is also evidence to suggest that widows and widowers respond differentially to assessment measures.

Colposcopy information leaflets: what women want to know and when they want to receive this information

OBJECTIVES

To evaluate whether the information leaflets produced by UK colposcopy clinics provide women with the information they desire and to determine when they would like to receive this information.

DESIGN

Questionnaire study and structured evaluation.

SETTING

The colposcopy clinic of a UK cancer centre.

PARTICIPANTS

Forty-two women attending a pre-colposcopy counselling session and 100 consecutive women attending the colposcopy clinic.

METHODS

Thirty-eight standards derived from the concerns/questions asked by women attending a pre-colposcopy counselling session were used to assess locally produced colposcopy clinic leaflets from UK colposcopy clinics, the leaflets produced by the Royal College of Obstetricians and Gynaecologists and the National Health Service Cervical Screening Programme (NHSCSP), and two "leaflets" obtained from internet sites. The Gunning fog test was used to assess the leaflets' readability. A questionnaire survey of 100 women attending the colposcopy clinic was used to determine when women wanted to receive information about colposcopy.

MAIN OUTCOME MEASURES

Percentage of questions answered by a given leaflet and Gunning fog scores for readability.

RESULTS

The information leaflets of 128 colposcopy clinics were received and assessed. Thirty-two clinics only sent women the NHSCSP leaflet. No leaflet answered all 38 questions. Less than half (36/100) of the leaflets answered more than 50% of the questions. In addition to the lack of advice given, different leaflets frequently gave conflicting advice. The average Gunning fog score was 9.7 (range 5.5-15.5). The majority of women (70%) wanted to receive information about colposcopy at or prior to the time of receiving their abnormal smear test result, although only 42% of women actually received information at this time.

CONCLUSIONS

Many UK colposcopy clinics do not appear to be providing women with the information they require to understand their condition and the procedure that they are about to undergo. Furthermore, this information is often not provided at the appropriate time in the screening process.

Can pre-colposcopy sessions reduce anxiety at the time of colposcopy? A prospective randomised study

The main objective of this prospective randomised study was to evaluate whether offering pre-colposcopy group sessions reduces anxiety at the time of colposcopy. We also examined whether this strategy improved knowledge about abnormal smears and colposcopy and improved satisfaction with the colposcopy service provided. One hundred and forty-seven women undergoing colposcopy for the first time were randomised into two groups. The control group (n = 75) received conventional management. The study group (n = 72), in addition to conventional management, were invited to attend a pre-colposcopy group session led by a trained colposcopy nurse. Questionnaires were used to determine state anxiety inventory scores and knowledge scores at the time of randomisation, immediately before colposcopy and 6 weeks after the clinic visit. Satisfaction questionnaires were completed 6 weeks after the clinic visit. We found that women attending colposcopy clinics are anxious. Those women who attended the pre-colposcopy session had improved knowledge scores (P = 0.039) at the time of colposcopy and satisfaction (P = 0.037). However, the intervention failed to significantly reduce anxiety at the time of colposcopy (P > 0.05).

Richard Pattle and the early days of comparative respiratory biology

After a distinguished academic career Richard Pattle was well qualified as a Physicist and Physical Chemist. Following a brief spell in the army (REME) he was appointed as a biophysicist in the Chemical Defence Experimental Laboratory at Porton Down. Such centres here and in USA initiated work which led to important discoveries relating to lung surfactants. His 1955 paper led to an expansion in studies of surfactant and the lung lining and appreciation of its importance in Respiratory Distress Syndrome of the newborn. This period (late 1950s and 1960s) was also when comparative studies in respiration were beginning to expand. Because of the relevance of his work to RDS it was inevitable that Pattle became closely associated with Human Physiologists. But he was also interested in other mammals and later the lungs of primates and lower vertebrates and this was extended to fish and the transition from water to air breathing. It was a time notable for the getting-together of medically-trained and zoologically-based scientists. I was always fascinated by the simplicity of his bubble technique and encouraged him to summarise it and his findings at Bhagalpur in 1974 when he introduced the term 'surpellic'. He also made important studies on bird lung surfactants.

Study of the effects of brief exposure to an organophosphorous insecticide (methidathion) on blood characteristics of carp (Cyprinus carpio)

1. Carp with the dorsal aorta cannulated have been used to study the effects of a 5 hour exposure to methidathion at concentrations of 2 or 6 mg/litre. Blood samples taken during control, exposed and recovery periods were used to determine plasma acetylcholinesterase (AChE) activity, haematocrit value, mean cell volume, total plasma protein, and filtration time through Nucleopore filters containing pores of 8 microns. 2. A drastic inhibition of AChE activity was observed and this continued during the recovery for at least 5 days. Changes in other blood parameters were less marked and recovered soon after fish were returned to non-polluted water. The increased haematocrit and decrease in filtration time suggests some impairment of gas transfer during exposure but the depletion of AChE activity and associated muscular and neural disturbances are more serious results of pollution with this organophosphate insecticide.

A comparison of blood composition in the common carp (Cyprinus carpio L.) using samples obtained from caudal vessels and dorsal aorta cannulae

Blood samples have been taken from chronic cannulae in the dorsal aorta of Carp under anaesthesia and at intervals of 24, 29 and 101 hours later. A few minutes after the 29 hours and 101 hours sampling was completed similar quantities (1 ml) were withdrawn from a caudal vessel. All samples were treated in the same way and determinations made of haematocrit value, total protein, glucose concentration and activity of the following plasma enzymes-acetylcholinesterase (AChE), glutamic oxaloacetic transaminase (GOT), glutamic pyruvic transaminase (GPT) and lactic dehydrogenase (LDH). Comparisons between dorsal aorta and caudal samples at 29 hours and 101 hours showed no significant differences for any of the enzyme activities that were tested. However significant (paired "t"-test, P < 0.05) differences were observed for haematocrit and glucose concentration. It is suggested that the higher haematocrit and lower glucose of caudal samples are consistent with the expected differences between arterial and venous blood. Apart from these differences which reflect normal physiological function, values obtained for samples from the two sites were identical. It is concluded that blood sampling from caudal vessels is a satisfactory method to obtain blood for biomonitoring purposes at least for those enzymes commonly used in recent surveys.

The relationship of sustained exercise training and bone mineral density in aging male runners

Fifty-six men aged 42-73 years (50.2 +/- 10.0 years), who were competitive distance runners 20-25 years previously, were examined for bone mineral density (BMD) to determine the relationship between sustained distance running and BMD. Subjects were classified as being highly trained (HT, n = 17), moderately trained (MT, n = 29) or untrained (UT, n = 10) according to their training in recent years. Subjects in each group were of similar age (HT 46.5 +/- 2.01, MT 53.0 +/- 1.51, UT 46.7 +/- 2.44 years) and lean body mass. Total body weight (kg) and percentage fat, however, were significantly greater (P < 0.05) in the UT group than in either the MT or HT groups (UT 80.6 +/- 2.44 kg, 22.0 +/- 1.16%; MT 74.9 +/- 1.51 kg, 17.5 +/- 0.61%; HT 70.5 +/- 1.71 kg, 13.5 +/- 0.59%). Lumbar vertebrae and hip region BMD (g.cm-2) was determined via dual energy X-ray absorptiometry (DEXA). No differences in BMD were found among the three groups in either the lumbar (HT 1.00 +/- 0.02, MT 1.02 +/- 0.03, UT 1.07 +/- 0.04 g.cm-2) or the hip regions (HT 0.99 +/- 0.03, MT 0.98 +/- 0.02, UT 1.06 +/- 0.04 g.cm-2). Furthermore, none of the groups had BMD that was significantly different from age-matched normative values taken from a reference database. A moderate correlation was found between body weight and BMD when combining all subjects (r = 0.38 for lumbar and r = 0.41 for hip). These results indicate that middle-aged to older males who have sustained exercise training in the form of running do not have significantly different lumbar vertebrae or hip region BMD compared to individuals who run less or not at all.