Baiji genomes reveal low genetic variability and new insights into secondary aquatic adaptations

Telomere-to-telomere gap-free genome assembly of the endangered Yangtze finless porpoise and East Asian finless porpoise

BACKGROUND

The Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis, YFP) and the East Asian finless porpoise (Neophocaena asiaeorientalis sunameri, EFP) are 2 subspecies of the narrow-ridged finless porpoise that live in freshwater and saltwater, respectively. The main objective of this study was to provide contiguous chromosome-level genome assemblies for YFP and EFP.

RESULTS

Here, we generated and upgraded the genomes of YFP and EFP at the telomere-to-telomere level through the integration of PacBio HiFi long reads, ultra-long ONT reads, and Hi-C sequencing data with a total size of 2.48 Gb and 2.50 Gb, respectively. The scaffold N50 of 2 genomes was 125.12 Mb (YFP) and 128 Mb (EFP) with 1 contig for 1 chromosome. The telomere repeat and centromere position were clearly identified in both YFP and EFP genomes. In total, 5,480 newfound genes were detected in the YFP genome, including 56 genes located in the newly identified centromere regions. Additionally, synteny blocks, structural similarities, phylogenetic relationships, gene family expansion, and inference of selection were studied in connection with the genomes of other related mammals.

CONCLUSIONS

Our research findings provide evidence for the gradual adaptation of EFP in a marine environment and the potential sensitivity of YFP to genetic damage. Compared to the 34 cetacean genomes sourced from public databases, the 2 new assemblies demonstrate superior continuity with the longest contig N50 and scaffold N50 values, as well as the lowest number of contigs. The improvement of telomere-to-telomere gap-free reference genome resources supports conservation genetics and population management for finless porpoises.

Population genomics of the muskox' resilience in the near absence of genetic variation

Genomic studies of species threatened by extinction are providing crucial information about evolutionary mechanisms and genetic consequences of population declines and bottlenecks. However, to understand how species avoid the extinction vortex, insights can be drawn by studying species that thrive despite past declines. Here, we studied the population genomics of the muskox (Ovibos moschatus), an Ice Age relict that was at the brink of extinction for thousands of years at the end of the Pleistocene yet appears to be thriving today. We analysed 108 whole genomes, including present-day individuals representing the current native range of both muskox subspecies, the white-faced and the barren-ground muskox (O. moschatus wardi and O. moschatus moschatus) and a ~21,000-year-old ancient individual from Siberia. We found that the muskox' demographic history was profoundly shaped by past climate changes and post-glacial re-colonizations. In particular, the white-faced muskox has the lowest genome-wide heterozygosity recorded in an ungulate. Yet, there is no evidence of inbreeding depression in native muskox populations. We hypothesize that this can be explained by the effect of long-term gradual population declines that allowed for purging of strongly deleterious mutations. This study provides insights into how species with a history of population bottlenecks, small population sizes and low genetic diversity survive against all odds.

New chromosome-scale genomes provide insights into marine adaptations of sea snakes (Hydrophis: Elapidae)

BACKGROUND

Sea snakes underwent a complete transition from land to sea within the last ~ 15 million years, yet they remain a conspicuous gap in molecular studies of marine adaptation in vertebrates.

RESULTS

Here, we generate four new annotated sea snake genomes, three of these at chromosome-scale (Hydrophis major, H. ornatus and H. curtus), and perform detailed comparative genomic analyses of sea snakes and their closest terrestrial relatives. Phylogenomic analyses highlight the possibility of near-simultaneous speciation at the root of Hydrophis, and synteny maps show intra-chromosomal variations that will be important targets for future adaptation and speciation genomic studies of this system. We then used a strict screen for positive selection in sea snakes (against a background of seven terrestrial snake genomes) to identify genes over-represented in hypoxia adaptation, sensory perception, immune response and morphological development.

CONCLUSIONS

We provide the best reference genomes currently available for the prolific and medically important elapid snake radiation. Our analyses highlight the phylogenetic complexity and conserved genome structure within Hydrophis. Positively selected marine-associated genes provide promising candidates for future, functional studies linking genetic signatures to the marine phenotypes of sea snakes and other vertebrates.

Genomic Divergence and the Evolution of Ecotypes in Bottlenose Dolphins (Genus Tursiops)

Climatic changes have caused major environmental restructuring throughout the world's oceans. Marine organisms have responded to novel conditions through various biological systems, including genomic adaptation. Growing accessibility of next-generation DNA sequencing methods to study nonmodel species has recently allowed genomic changes underlying environmental adaptations to be investigated. This study used double-digest restriction-site associated DNA (ddRAD) sequence data to investigate the genomic basis of ecotype formation across currently recognized species and subspecies of bottlenose dolphins (genus Tursiops) in the Southern Hemisphere. Subspecies-level genomic divergence was confirmed between the offshore common bottlenose dolphin (T. truncatus truncatus) and the inshore Lahille's bottlenose dolphin (T. t. gephyreus) from the southwestern Atlantic Ocean (SWAO). Similarly, subspecies-level divergence is suggested between inshore (eastern Australia) Indo-Pacific bottlenose dolphin (T. aduncus) and the proposed Burrunan dolphin (T. australis) from southern Australia. Inshore bottlenose dolphin lineages generally had lower genomic diversity than offshore lineages, a pattern particularly evident for T. t. gephyreus, which showed exceptionally low diversity. Genomic regions associated with cardiovascular, musculoskeletal, and energy production systems appear to have undergone repeated adaptive evolution in inshore lineages across the Southern Hemisphere. We hypothesize that comparable selective pressures in the inshore environment drove similar adaptive responses in each lineage, supporting parallel evolution of inshore bottlenose dolphins. With climate change altering marine ecosystems worldwide, it is crucial to gain an understanding of the adaptive capacity of local species and populations. Our study provides insights into key adaptive pathways that may be important for the long-term survival of cetaceans and other organisms in a changing marine environment.

When morphology does not fit the genomes: the case of rodent olfaction

Linking genes to phenotypes has been a major question in evolutionary biology for the last decades. In the genomic era, few studies attempted to link olfactory-related genes to different anatomical proxies. However, they found very inconsistent results. This study is the first to investigate a potential relation between olfactory turbinals and olfactory receptor (OR) genes. We demonstrated that despite the use of similar methodology in the acquisition of data, OR genes do not correlate with the relative and the absolute surface area of olfactory turbinals. These results challenged the interpretations of several studies based on different proxies related to olfaction and their potential relation to olfactory capabilities.

Gapless genome assembly of East Asian finless porpoise

In recent years, conservation efforts have increased for rare and endangered aquatic wildlife, especially cetaceans. However, the East Asian finless porpoise (Neophocaena asiaeorientalis sunameri), which has a wide distribution in China, has received far less attention and protection. As an endangered small cetacean, the lack of a chromosomal-level reference for the East Asian finless porpoise limits our understanding of its population genetics and conservation biology. To address this issue, we combined PacBio HiFi long reads and Hi-C sequencing data to generate a gapless genome of the East Asian finless porpoise that is approximately 2.5 Gb in size over its 21 autosomes and two sex chromosomes (X and Y). A total of 22,814 protein-coding genes were predicted where ~97.31% were functionally annotated. This high-quality genome assembly of East Asian finless porpoise will not only provide new resources for the comparative genomics of cetaceans and conservation biology of threatened species, but also lay a foundation for more speciation, ecology, and evolutionary studies. Measurement(s) Neophocaena asiaeorientalis sunameri • Gapless genome assembly • sequence annotation Technology Type(s) MGISEQ. 2000 • PacBio HiFi Sequencing • Hi-C Sample Characteristic - Organism Neophocaena asiaeorientalis sunameri Sample Characteristic - Environment seawater Sample Characteristic - Location Yellow Sea near Lianyungang City, Jiangsu Province, China.

Accelerated evolution of dim-light vision-related arrestin in deep-diving amniotes

Arrestins are key molecules involved in the signaling of light-sensation initiated by visual pigments in retinal photoreceptor cells. Vertebrate photoreceptor cells have two types of arrestins, rod arrestin, which is encoded by SAG and is expressed in both rods and cones, and cone arrestin, encoded by ARR3 in cones. The arrestins can bind to visual pigments, and thus regulate either dim-light vision via interactions with rhodopsin or bright-light vision together with cone visual pigments. After adapting to terrestrial life, several amniote lineages independently went back to the sea and evolved deep-diving habits. Interestingly, the rhodopsins in these species exhibit specialized phenotypes responding to rapidly changing dim-light environments. However, little is known about whether their rod arrestin also experienced adaptive evolution associated with rhodopsin. Here, we collected SAG coding sequences from >250 amniote species, and examined changes in selective pressure experienced by the sequences from deep-diving taxa. Divergent patterns of evolution of SAG were observed in the penguin, pinniped and cetacean clades, suggesting possible co-adaptation with rhodopsin. After verifying pseudogenes, the same analyses were performed for cone arrestin (ARR3) in deep-diving species and only sequences from cetacean species, and not pinnipeds or penguins, have experienced changed selection pressure compared to other species. Taken together, this evidence for changes in selective pressures acting upon arrestin genes strengthens the suggestion that rapid dim-light adaptation for deep-diving amniotes require SAG, but not ARR3.

Chromosomal-scale genome assembly of the near-extinction big-head schizothorcin (Aspiorhynchus laticeps)

The big-head schizothorcin (Aspiorhynchus laticeps) is an endemic and near-extinction freshwater fish in Xinjiang, China. In this study, a chromosome-scale genome assembly of A. laticeps was generated using PacBio and Hi-C techniques. The PacBio sequencing data resulted in a 1.58 Gb assembly with a contig N50 of 1.27 Mb. Using Hi-C scaffolding approach, 88.38% of the initial assembled sequences were anchored and oriented into a chromosomal-scale assembly. The final assembly consisted of 25 pseudo-chromosomes that yielded 1.37 Gb of sequence, with a scaffold N50 of 44.02 Mb. BUSCO analysis showed a completeness score of 93.7%. The genome contained 48,537 predicted protein-coding genes and 58.31% of the assembly was annotated as repetitive sequences. Whole genome duplication events were further confirmed using 4dTv analysis. The genome assembly of A. laticeps should be valuable and important to understand the genetic adaptation and endangerment process of this species, which could lead to more effective management and conservation of the big-head schizothorcin and related freshwater fish species.

Measurement(s)	Genome Assembly Sequence
Technology Type(s)	Next Generation Sequencing
Sample Characteristic - Organism	Aspiorhynchus laticeps
Sample Characteristic - Location	Xinjiang Autonomous Region

Seascape genomics of coastal bottlenose dolphins along strong gradients of temperature and salinity

Heterogeneous seascapes and strong environmental gradients in coastal waters are expected to influence adaptive divergence, particularly in species with large population sizes where selection is expected to be highly efficient. However, these influences might also extend to species characterized by strong social structure, natal philopatry and small home ranges. We implemented a seascape genomic study to test this hypothesis in Indo-Pacific bottlenose dolphins (Tursiops aduncus) distributed along the environmentally heterogeneous coast of southern Australia. The datasets included oceanographic and environmental variables thought to be good predictors of local adaptation in dolphins and 8,081 filtered single nucleotide polymorphisms (SNPs) genotyped for individuals sampled from seven different bioregions. From a neutral perspective, population structure and connectivity of the dolphins were generally influenced by habitat type and social structuring. Genotype-environment association analysis identified 241 candidate adaptive loci and revealed that sea surface temperature and salinity gradients influenced adaptive divergence in these animals at both large- (1,000s km) and fine-scales (<100 km). Enrichment analysis and annotation of candidate genes revealed functions related to sodium-activated ion transport, kidney development, adipogenesis and thermogenesis. The findings of spatial adaptive divergence and inferences of putative physiological adaptations challenge previous suggestions that marine megafauna is most likely to be affected by environmental and climatic changes via indirect, trophic effects. Our work contributes to conservation management of coastal bottlenose dolphins subjected to anthropogenic disturbance and to efforts of clarifying how seascape heterogeneity influences adaptive diversity and evolution in small cetaceans.

Ancient demographics determine the effectiveness of genetic purging in endangered lizards

The purging of deleterious alleles has been hypothesized to mitigate inbreeding depression, but its effectiveness in endangered species remains debatable. To understand how deleterious alleles are purged during population contractions, we analyzed genomes of the endangered Chinese crocodile lizard (Shinisaurus crocodilurus), which is the only surviving species of its family and currently isolated into small populations. Population genomic analyses revealed four genetically distinct conservation units and sharp declines in both effective population size and genetic diversity. By comparing the relative genetic load across populations and conducting genomic simulations, we discovered that seriously deleterious alleles were effectively purged during population contractions in this relict species, although inbreeding generally enhanced the genetic burden. However, despite with the initial purging, our simulations also predicted that seriously deleterious alleles will gradually accumulate under prolonged bottlenecking. Therefore, we emphasize the importance of maintaining a minimum population capacity and increasing the functional genetic diversity in conservation efforts to preserve populations of the crocodile lizard and other endangered species.

Vulture Genomes Reveal Molecular Adaptations Underlying Obligate Scavenging and Low Levels of Genetic Diversity

Obligate scavenging on the dead and decaying animal matter is a rare dietary specialization that in extant vertebrates is restricted to vultures. These birds perform essential ecological services, yet many vulture species have undergone recent steep population declines and are now endangered. To test for molecular adaptations underlying obligate scavenging in vultures, and to assess whether genomic features might have contributed to their population declines, we generated high-quality genomes of the Himalayan and bearded vultures, representing both independent origins of scavenging within the Accipitridae, alongside a sister taxon, the upland buzzard. By comparing our data to published sequences from other birds, we show that the evolution of obligate scavenging in vultures has been accompanied by widespread positive selection acting on genes underlying gastric acid production, and immunity. Moreover, we find evidence of parallel molecular evolution, with amino acid replacements shared among divergent lineages of these scavengers. Our genome-wide screens also reveal that both the Himalayan and bearded vultures exhibit low levels of genetic diversity, equating to around a half of the mean genetic diversity of other bird genomes examined. However, demographic reconstructions indicate that population declines began at around the Last Glacial Maximum, predating the well-documented dramatic declines of the past three decades. Taken together, our genomic analyses imply that vultures harbor unique adaptations for processing carrion, but that modern populations are genetically depauperate and thus especially vulnerable to further genetic erosion through anthropogenic activities.

The global significance of biodiversity science in China: an overview

Biodiversity science in China has seen rapid growth over recent decades, ranging from baseline biodiversity studies to understanding the processes behind evolution across dynamic regions such as the Qinghai-Tibetan Plateau. We review research, including species catalogues; biodiversity monitoring; the origins, distributions, maintenance and threats to biodiversity; biodiversity-related ecosystem function and services; and species and ecosystems' responses to global change. Next, we identify priority topics and offer suggestions and priorities for future biodiversity research in China. These priorities include (i) the ecology and biogeography of the Qinghai-Tibetan Plateau and surrounding mountains, and that of subtropical and tropical forests across China; (ii) marine and inland aquatic biodiversity; and (iii) effective conservation and management to identify and maintain synergies between biodiversity and socio-economic development to fulfil China's vision for becoming an ecological civilization. In addition, we propose three future strategies: (i) translate advanced biodiversity science into practice for biodiversity conservation; (ii) strengthen capacity building and application of advanced technologies, including high-throughput sequencing, genomics and remote sensing; and (iii) strengthen and expand international collaborations. Based on the recent rapid progress of biodiversity research, China is well positioned to become a global leader in biodiversity research in the near future.

The inflated significance of neutral genetic diversity in conservation genetics

The current rate of species extinction is rapidly approaching unprecedented highs, and life on Earth presently faces a sixth mass extinction event driven by anthropogenic activity, climate change, and ecological collapse. The field of conservation genetics aims at preserving species by using their levels of genetic diversity, usually measured as neutral genome-wide diversity, as a barometer for evaluating population health and extinction risk. A fundamental assumption is that higher levels of genetic diversity lead to an increase in fitness and long-term survival of a species. Here, we argue against the perceived importance of neutral genetic diversity for the conservation of wild populations and species. We demonstrate that no simple general relationship exists between neutral genetic diversity and the risk of species extinction. Instead, a better understanding of the properties of functional genetic diversity, demographic history, and ecological relationships is necessary for developing and implementing effective conservation genetic strategies.

Phylogenomic Resolution of the Cetacean Tree of Life Using Target Sequence Capture

The evolution of cetaceans, from their early transition to an aquatic lifestyle to their subsequent diversification, has been the subject of numerous studies. However, although the higher-level relationships among cetacean families have been largely settled, several aspects of the systematics within these groups remain unresolved. Problematic clades include the oceanic dolphins (37 spp.), which have experienced a recent rapid radiation, and the beaked whales (22 spp.), which have not been investigated in detail using nuclear loci. The combined application of high-throughput sequencing with techniques that target specific genomic sequences provide a powerful means of rapidly generating large volumes of orthologous sequence data for use in phylogenomic studies. To elucidate the phylogenetic relationships within the Cetacea, we combined sequence capture with Illumina sequencing to generate data for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\sim $\end{document}3200 protein-coding genes for 68 cetacean species and their close relatives including the pygmy hippopotamus. By combining data from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$>$\end{document}38,000 exons with existing sequences from 11 cetaceans and seven outgroup taxa, we produced the first comprehensive comparative genomic data set for cetaceans, spanning 6,527,596 aligned base pairs (bp) and 89 taxa. Phylogenetic trees reconstructed with maximum likelihood and Bayesian inference of concatenated loci, as well as with coalescence analyses of individual gene trees, produced mostly concordant and well-supported trees. Our results completely resolve the relationships among beaked whales as well as the contentious relationships among oceanic dolphins, especially the problematic subfamily Delphinidae. We carried out Bayesian estimation of species divergence times using MCMCTree and compared our complete data set to a subset of clocklike genes. Analyses using the complete data set consistently showed less variance in divergence times than the reduced data set. In addition, integration of new fossils (e.g., Mystacodon selenensis) indicates that the diversification of Crown Cetacea began before the Late Eocene and the divergence of Crown Delphinidae as early as the Middle Miocene. [Cetaceans; phylogenomics; Delphinidae; Ziphiidae; dolphins; whales.]

An Indo-Pacific Humpback Dolphin Genome Reveals Insights into Chromosome Evolution and the Demography of a Vulnerable Species

The Indo-Pacific humpback dolphin (Sousa chinensis) is a small inshore species of odontocete cetacean listed as Vulnerable on the IUCN Red List. Here, we report on the evolution of S. chinensis chromosomes from its cetruminant ancestor and elucidate the evolutionary history and population genetics of two neighboring S. chinensis populations. We found that breakpoints in ancestral chromosomes leading to S. chinensis could have affected the function of genes related to kidney filtration, body development, and immunity. Resequencing of individuals from two neighboring populations in the northwestern South China Sea, Leizhou Bay and Sanniang Bay, revealed genetic differentiation, low diversity, and small contemporary effective population sizes. Demographic analyses showed a marked decrease in the population size of the two investigated populations over the last ~4,000 years, possibly related to climatic oscillations. This study implies a high risk of extinction and strong conservation requirement for the Indo-Pacific humpback dolphin.

•

Deducing chromosome evolution from ancestral Cetruminantia and ancestral Odontoceti

•

Reconstructing the demographic history of Sousa chinensis

•

Implying high risk of extinction and strong conservation requirement for S. chinensis

Characterization of the primate TRIM gene family reveals the recent evolution in primates

The tripartite motif (TRIM) gene family encodes diverse distinct proteins that play important roles in many biological processes. However, the molecular evolution and phylogenetic relationships of TRIM genes in primates are still elusive. We performed a genomic approach to identify and characterize TRIM genes in human and other six primate genomes. In total, 537 putative functional TRIM genes were identified and TRIM members varied among primates. A neighbor joining (NJ) tree based on the protein sequences of 82 human TRIM genes indicates seven TRIM groups, which is consistent with the results based on the architectural motifs. Many TRIM gene duplication events were identified, indicating a recent expansion of TRIM family in primate lineages. Interestingly, the chimpanzee genome shows the greatest TRIM gene expansion among the primates; however, its congeneric species, bonobo, has the least number of TRIM genes and no duplication event. Moreover, we identified a ~ 200 kb deletion on chromosome 11 of bonobos that results in a loss of cluster3 TRIM genes. The loss of TRIM genes might have occurred within the last 2 mys. Analysis of positive selection recovered 9 previously reported and 21 newly identified positively selected TRIM genes. In particular, most positive selected sites are located in the B30.2 domains. Our results have provided new insight into the evolution of primate TRIM genes and will broaden our understanding on the functions of the TRIM family.

Comparative transcriptomic analysis of fireflies (Coleoptera: Lampyridae) to explore the molecular adaptations to fresh water

Aquatic insects are well adapted to freshwater environments, but the molecular basis of these adaptations remains largely unknown. Most firefly species (Coleoptera: Lampyridae) are terrestrial, but the larvae of several species are aquatic. Here, larval and adult transcriptomes from Aquatica leii (freshwater) and Lychnuris praetexta (terrestrial) were generated to test whether the genes associated with metabolic efficiency and morphology have undergone adaptive evolution to fresh water. The aquatic fireflies had a significantly lower ratio of nonsynonymous to synonymous substitutions than the terrestrial insects, indicating a genomewide evolutionary constraint in the aquatic fireflies. We identified 341 fast-evolving genes and 116 positively selected genes in the aquatic fireflies. Of these, 76 genes exhibiting both fast evolution and positive selection were primarily involved in ATP production, energy metabolism and the hypoxia response. We identified 7,271 differentially expressed genes (DEGs) in A. leii (adults versus larvae) and 8,309 DEGs in L. praetexta (adults versus larvae). DEGs specific to the aquatic firefly (n = 1,445) were screened via interspecific comparisons (A. leii versus L. praetexta) and were significantly enriched for genes involved in metabolic efficiency (e.g., ATP production, hypoxia, and immune responses) and certain aspects of morphology (e.g., cuticle chitin, tracheal and compound eye morphology). These results indicate that sequence and expression-level changes in genes associated with both metabolic efficiency and morphological attributes related to the freshwater lifestyle contributed to freshwater adaptation in fireflies. This study provides new insights into the molecular mechanisms of aquatic adaptation in insects.

Diving deep: understanding the genetic components of hypoxia tolerance in marine mammals

Marine mammals have highly specialized physiology, exhibited in many species by extreme breath-holding capabilities that allow deep dives and extended submergence. Cardiovascular control and cell-level hypoxia tolerance are key features of this phenotype. Identifying genomic signatures tied to physiology will be valuable in understanding these natural model species, which may generate translational opportunities to human diseases arising from hypoxic stress or tissue injury. Genomic analyses have now been conducted in dolphins, river dolphins, minke whales, bowhead whales, and polar bears, with multispecies studies exploring evolutionary signals across marine mammal lineages, encompassing extinct and extant divers. Single-species genome studies for sirenians do not yet exist. Extant marine mammals arose in three lineages from separate aquatic recolonizations. Their physiological specializations, along with these independent origins create an interesting case to examine convergent evolution. Although molecular mechanisms of hypoxia tolerance are not universally apparent across marine mammal genomic studies, altered evolutionary rates have been identified for genes linked to oxygen binding and transport (e.g., MB, HBA, and HBB), blood pressure control (e.g., endothelin pathway genes), and cell protection in multiple species. Despite convergent phenotypes across clades, instances of identical molecular convergence have been uncommon. Given the inherent logistical and regulatory difficulties associated with functional genetic experiments in marine mammals, several avenues of further investigation are suggested to enable validation of candidate genes for hypoxia tolerance: leveraging phylogeny to better understand convergent phenotypes; ontogenic studies to identify regulation of key genes underlying the elite, adult, hypoxia-tolerant physiology; and cell culture manipulations to understand gene function.

Genes lost during the transition from land to water in cetaceans highlight genomic changes associated with aquatic adaptations

The transition from land to water in whales and dolphins (cetaceans) was accompanied by remarkable adaptations. To reveal genomic changes that occurred during this transition, we screened for protein-coding genes that were inactivated in the ancestral cetacean lineage. We found 85 gene losses. Some of these were likely beneficial for cetaceans, for example, by reducing the risk of thrombus formation during diving (F12 and KLKB1), erroneous DNA damage repair (POLM), and oxidative stress-induced lung inflammation (MAP3K19). Additional gene losses may reflect other diving-related adaptations, such as enhanced vasoconstriction during the diving response (mediated by SLC6A18) and altered pulmonary surfactant composition (SEC14L3), while loss of SLC4A9 relates to a reduced need for saliva. Last, loss of melatonin synthesis and receptor genes (AANAT, ASMT, and MTNR1A/B) may have been a precondition for adopting unihemispheric sleep. Our findings suggest that some genes lost in ancestral cetaceans were likely involved in adapting to a fully aquatic lifestyle.

Whole Genome Sequencing of Chinese White Dolphin (Sousa chinensis) for High-Throughput Screening of Antihypertensive Peptides

Chinese white dolphin (Sousa chinensis), also known as the Indo-Pacific humpback dolphin, has been classified as “Vulnerable” on the IUCN Red List of Threatened Species. It is a special cetacean species that lives in tropical and subtropical nearshore waters, with significant differences from other cetaceans. Here, we sequenced and assembled a draft genome of the Chinese white dolphin with a total length of 2.3 Gb and annotation of 18,387 protein-coding genes. Genes from certain expanded families are potentially involved in DNA replication and repairing, suggesting that they may be related to adaptation of this marine mammal to nearshore environments. We also discovered that its historical population had undergone a remarkable bottleneck incident before the Mindel glaciation. In addition, a comparative genomic survey on antihypertensive peptides (AHTPs) among five representative mammals with various residential habitats (such as remarkable differences in exogenous ion concentrations and sea depth) revealed that these small bioactive peptides were highly conserved among these examined mammals, and they had the most abundant hits in collagen subunit proteins, especially for two putative AHTP peptides Gly-Leu-Pro (GLP) and Leu-Gly-Pro (LGP). Our genome assembly will be a valuable resource for further genetic researches on adaptive ecology and conservation biology of cetaceans, and for in-depth investigations into bioactive peptides in aquatic and terrestrial mammals for development of peptide-based drugs to treat various human cardiovascular diseases.

Complete Inactivation of Sebum-Producing Genes Parallels the Loss of Sebaceous Glands in Cetacea

Genomes are dynamic biological units, with processes of gene duplication and loss triggering evolutionary novelty. The mammalian skin provides a remarkable case study on the occurrence of adaptive morphological innovations. Skin sebaceous glands (SGs), for instance, emerged in the ancestor of mammals serving pivotal roles, such as lubrication, waterproofing, immunity, and thermoregulation, through the secretion of sebum, a complex mixture of various neutral lipids such as triacylglycerol, free fatty acids, wax esters, cholesterol, and squalene. Remarkably, SGs are absent in a few mammalian lineages, including the iconic Cetacea. We investigated the evolution of the key molecular components responsible for skin sebum production: Dgat2l6, Awat1, Awat2, Elovl3, Mogat3, and Fabp9. We show that all analyzed genes have been rendered nonfunctional in Cetacea species (toothed and baleen whales). Transcriptomic analysis, including a novel skin transcriptome from blue whale, supports gene inactivation. The conserved mutational pattern found in most analyzed genes, indicates that pseudogenization events took place prior to the diversification of modern Cetacea lineages. Genome and skin transcriptome analysis of the common hippopotamus highlighted the convergent loss of a subset of sebum-producing genes, notably Awat1 and Mogat3. Partial loss profiles were also detected in non-Cetacea aquatic mammals, such as the Florida manatee, and in terrestrial mammals displaying specialized skin phenotypes such as the African elephant, white rhinoceros and pig. Our findings reveal a unique landscape of “gene vestiges” in the Cetacea sebum-producing compartment, with limited gene loss observed in other mammalian lineages: suggestive of specific adaptations or specializations of skin lipids.

Differential Evolution of the Epidermal Keratin Cytoskeleton in Terrestrial and Aquatic Mammals

Keratins are the main intermediate filament proteins of epithelial cells. In keratinocytes of the mammalian epidermis they form a cytoskeleton that resists mechanical stress and thereby are essential for the function of the skin as a barrier against the environment. Here, we performed a comparative genomics study of epidermal keratin genes in terrestrial and fully aquatic mammals to determine adaptations of the epidermal keratin cytoskeleton to different environments. We show that keratins K5 and K14 of the innermost (basal), proliferation-competent layer of the epidermis are conserved in all mammals investigated. In contrast, K1 and K10, which form the main part of the cytoskeleton in the outer (suprabasal) layers of the epidermis of terrestrial mammals, have been lost in whales and dolphins (cetaceans) and in the manatee. Whereas in terrestrial mammalian epidermis K6 and K17 are expressed only upon stress-induced epidermal thickening, high levels of K6 and K17 are consistently present in dolphin skin, indicating constitutive expression and substitution of K1 and K10. K2 and K9, which are expressed in a body site-restricted manner in human and mouse suprabasal epidermis, have been lost not only in cetaceans and manatee but also in some terrestrial mammals. The evolution of alternative splicing of K10 and differentiation-dependent upregulation of K23 have increased the complexity of keratin expression in the epidermis of terrestrial mammals. Taken together, these results reveal evolutionary diversification of the epidermal cytoskeleton in mammals and suggest a complete replacement of the quantitatively predominant epidermal proteins of terrestrial mammals by originally stress-inducible keratins in cetaceans.

Evidence that DNA repair genes, a family of tumor suppressor genes, are associated with evolution rate and size of genomes

Adaptive radiation and evolutionary stasis are characterized by very different evolution rates. The main aim of this study was to investigate if any genes have a special role to a high or low evolution rate. The availability of animal genomes permitted comparison of gene content of genomes of 24 vertebrate species that evolved through adaptive radiation (representing high evolutionary rate) and of 20 vertebrate species that are considered as living fossils (representing a slow evolutionary rate or evolutionary stasis). Mammals, birds, reptiles, and bony fishes were included in the analysis. Pathway analysis was performed for genes found to be specific in adaptive radiation or evolutionary stasis respectively. Pathway analysis revealed that DNA repair and cellular response to DNA damage are important (false discovery rate = 8.35 × 10⁻⁵; 7.15 × 10⁻⁶, respectively) for species evolved through adaptive radiation. This was confirmed by further genetic in silico analysis (p = 5.30 × 10⁻³). Nucleotide excision repair and base excision repair were the most significant pathways. Additionally, the number of DNA repair genes was found to be linearly related to the genome size and the protein number (proteome) of the 44 animals analyzed (p < 1.00 × 10⁻⁴), this being compatible with Drake’s rule. This is the first study where radiated and living fossil species have been genetically compared. Evidence has been found that cancer-related genes have a special role in radiated species. Linear association of the number of DNA repair genes with the species genome size has also been revealed. These comparative genetics results can support the idea of punctuated equilibrium evolution.

The genome resources for conservation of Indo-Pacific humpback dolphin, Sousa chinensis

The Indo-Pacific humpback dolphin (Sousa chinensis), is a threatened marine mammal and belongs to the First Order of the National Key Protected Wild Aquatic Animals List in China. However, limited genomic information is available for studies of its population genetics and biological conservation. Here, we have assembled a genomic sequence of this species using a whole genome shotgun (WGS) sequencing strategy after a pilot low coverage genome survey. The total assembled genome size was 2.34 Gb: with a contig N50 of 67 kb and a scaffold N50 of 9 Mb (107.6-fold sequencing coverage). The S. chinensis genome contained 24,640 predicted protein-coding genes and had approximately 37% repeated sequences. The completeness of the genome assembly was evaluated by benchmarking universal single copy orthologous genes (BUSCOs): 94.3% of a total 4,104 expected mammalian genes were identified as complete, and 2.3% were identified as fragmented. This newly produced high-quality assembly and annotation of the genome will greatly promote the future studies of the genetic diversity, conservation and evolution.

The Singularity of Cetacea Behavior Parallels the Complete Inactivation of Melatonin Gene Modules

Melatonin, the hormone of darkness, is a peculiar molecule found in most living organisms. Emerging as a potent broad-spectrum antioxidant, melatonin was repurposed into extra roles such as the modulation of circadian and seasonal rhythmicity, affecting numerous aspects of physiology and behaviour, including sleep entrainment and locomotor activity. Interestingly, the pineal gland—the melatonin synthesising organ in vertebrates—was suggested to be absent or rudimentary in some mammalian lineages, including Cetacea. In Cetacea, pineal regression is paralleled by their unique bio-rhythmicity, as illustrated by the unihemispheric sleeping behaviour and long-term vigilance. Here, we examined the genes responsible for melatonin synthesis (Aanat and Asmt) and signalling (Mtnr1a and Mtnr1b) in 12 toothed and baleen whale genomes. Based on an ample genomic comparison, we deduce that melatonin-related gene modules are eroded in Cetacea.

Diversity and Evolution of Sensor Histidine Kinases in Eukaryotes

Histidine kinases (HKs) are primary sensor proteins that act in cell signaling pathways generically referred to as "two-component systems" (TCSs). TCSs are among the most widely distributed transduction systems used by both prokaryotic and eukaryotic organisms to detect and respond to a broad range of environmental cues. The structure and distribution of HK proteins are now well documented in prokaryotes, but information is still fragmentary for eukaryotes. Here, we have taken advantage of recent genomic resources to explore the structural diversity and the phylogenetic distribution of HKs in the prominent eukaryotic supergroups. Searches of the genomes of 67 eukaryotic species spread evenly throughout the phylogenetic tree of life identified 748 predicted HK proteins. Independent phylogenetic analyses of predicted HK proteins were carried out for each of the major eukaryotic supergroups. This allowed most of the compiled sequences to be categorized into previously described HK groups. Beyond the phylogenetic analysis of eukaryotic HKs, this study revealed some interesting findings: 1) characterization of some previously undescribed eukaryotic HK groups with predicted functions putatively related to physiological traits; 2) discovery of HK groups that were previously believed to be restricted to a single kingdom in additional supergroups, and 3) indications that some evolutionary paths have led to the appearance, transfer, duplication, and loss of HK genes in some phylogenetic lineages. This study provides an unprecedented overview of the structure and distribution of HKs in the Eukaryota and represents a first step toward deciphering the evolution of TCS signaling in living organisms.

Endangered Père David's deer genome provides insights into population recovering

The Milu (Père David's deer, Elaphurus davidianus) were once widely distributed in the swamps (coastal areas to inland areas) of East Asia. The dramatic recovery of the Milu population is now deemed a classic example of how highly endangered animal species can be rescued. However, the molecular mechanisms that underpinned this population recovery remain largely unknown. Here, different approaches (genome sequencing, resequencing, and salinity analysis) were utilized to elucidate the aforementioned molecular mechanisms. The comparative genomic analyses revealed that the largest recovered Milu population carries extensive genetic diversity despite an extreme population bottleneck. And the protracted inbreeding history might have facilitated the purging of deleterious recessive alleles. Seventeen genes that are putatively related to reproduction, embryonic (fatal) development, and immune response were under high selective pressure. Besides, SCNN1A, a gene involved in controlling reabsorption of sodium in the body, was positively selected. An additional 29 genes were also observed to be positively selected, which are involved in blood pressure regulation, cardiovascular development, cholesterol regulation, glycemic control, and thyroid hormone synthesis. It is possible that these genetic adaptations were required to buffer the negative effects commonly associated with a high-salt diet. The associated genetic adaptions are likely to have enabled increased breeding success and fetal survival. The future success of Milu population management might depend on the successful reintroduction of the animal to historically important distribution regions.

De-Extinction

De-extinction projects for species such as the woolly mammoth and passenger pigeon have greatly stimulated public and scientific interest, producing a large body of literature and much debate. To date, there has been little consistency in descriptions of de-extinction technologies and purposes. In 2016, a special committee of the International Union for the Conservation of Nature (IUCN) published a set of guidelines for de-extinction practice, establishing the first detailed description of de-extinction; yet incoherencies in published literature persist. There are even several problems with the IUCN definition. Here I present a comprehensive definition of de-extinction practice and rationale that expounds and reconciles the biological and ecological inconsistencies in the IUCN definition. This new definition brings together the practices of reintroduction and ecological replacement with de-extinction efforts that employ breeding strategies to recover unique extinct phenotypes into a single “de-extinction” discipline. An accurate understanding of de-extinction and biotechnology segregates the restoration of certain species into a new classification of endangerment, removing them from the purview of de-extinction and into the arena of species’ recovery. I term these species as “evolutionarily torpid species”; a term to apply to species falsely considered extinct, which in fact persist in the form of cryopreserved tissues and cultured cells. For the first time in published literature, all currently active de-extinction breeding programs are reviewed and their progress presented. Lastly, I review and scrutinize various topics pertaining to de-extinction in light of the growing body of peer-reviewed literature published since de-extinction breeding programs gained public attention in 2013.

Conservation genetics and genomics of threatened vertebrates in China

Conservation genetics and genomics are two independent disciplines that focus on using new techniques in genetics and genomics to solve problems in conservation biology. During the past two decades, conservation genetics and genomics have experienced rapid progress. Here, we summarize the research advances in the conservation genetics and genomics of threatened vertebrates (e.g., carnivorans, primates, ungulates, cetaceans, avians, amphibians and reptiles) in China. First, we introduce the concepts of conservation genetics and genomics and their development. Second, we review the recent advances in conservation genetics research, including noninvasive genetics and landscape genetics. Third, we summarize the progress in conservation genomics research, which mainly focuses on resolving genetic problems relevant to conservation such as genetic diversity, genetic structure, demographic history, and genomic evolution and adaptation. Finally, we discuss the future directions of conservation genetics and genomics.

Molecular Footprints of Aquatic Adaptation Including Bone Mass Changes in Cetaceans

Cetaceans (whales, dolphins, and porpoises) are a group of specialized mammals that evolved from terrestrial ancestors and are fully adapted to aquatic habitats. Taking advantage of the recently sequenced finless porpoise genome, we conducted comparative analyses of the genomes of seven cetaceans and related terrestrial species to provide insight into the molecular bases of adaptation of these aquatic mammals. Changes in gene sequences were identified in main lineages of cetaceans, offering an evolutionary picture of cetacean genomes that reveal new pathways that could be associated with adaptation to aquatic lifestyle. We profiled bone microanatomical structures across 28 mammals, including representatives of cetaceans, pinnipeds, and sirenians. Subsequent phylogenetic comparative analyses revealed genes (including leptin, insulin-like growth factor 1, and collagen type I alpha 2 chain) with the root-to-tip substitution rate significantly correlated with bone compactness, implicating these genes could be involved in bone mass control. Overall, this study described adjustments of the genomes of cetaceans according to lifestyle, phylogeny, and bone mass.

Ancient convergent losses of Paraoxonase 1 yield potential risks for modern marine mammals

Mammals diversified by colonizing drastically different environments, with each transition yielding numerous molecular changes, including losses of protein function. Though not initially deleterious, these losses could subsequently carry deleterious pleiotropic consequences. We have used phylogenetic methods to identify convergent functional losses across independent marine mammal lineages. In one extreme case, Paraoxonase 1 (PON1) accrued lesions in all marine lineages, while remaining intact in all terrestrial mammals. These lesions coincide with PON1 enzymatic activity loss in marine species' blood plasma. This convergent loss is likely explained by parallel shifts in marine ancestors' lipid metabolism and/or bloodstream oxidative environment affecting PON1's role in fatty acid oxidation. PON1 loss also eliminates marine mammals' main defense against neurotoxicity from specific man-made organophosphorus compounds, implying potential risks in modern environments.

High-quality whole-genome sequence of an abundant Holarctic odontocete, the harbour porpoise (Phocoena phocoena)

The harbour porpoise (Phocoena phocoena) is a highly mobile cetacean found across the Northern hemisphere. It occurs in coastal waters and inhabits basins that vary broadly in salinity, temperature and food availability. These diverse habitats could drive subtle differentiation among populations, but examination of this would be best conducted with a robust reference genome. Here, we report the first harbour porpoise genome, assembled de novo from an individual originating in the Kattegat Sea (Sweden). The genome is one of the most complete cetacean genomes currently available, with a total size of 2.39 Gb and 50% of the total length found in just 34 scaffolds. Using 122 of the longest scaffolds, we were able to show high levels of synteny with the genome of the domestic cattle (Bos taurus). Our draft annotation comprises 22,154 predicted genes, which we further annotated through matches to the NCBI nucleotide database, GO categorization and motif prediction. Within the predicted genes, we have confirmed the presence of >20 genes or gene families that have been associated with adaptive evolution in other cetaceans. Overall, this genome assembly and draft annotation represent a crucial addition to the genomic resources currently available for the study of porpoises and Phocoenidae evolution, phylogeny and conservation.

Molecular footprints of inshore aquatic adaptation in Indo-Pacific humpback dolphin (Sousa chinensis)

The Indo-Pacific humpback dolphin, Sousa chinensis, being a member of cetaceans, had fully adapted to inshore waters. As a threatened marine mammal, little molecular information available for understanding the genetic basis of ecological adaptation. We firstly sequenced and obtained the draft genome map of S. chinensis. Phylogenetic analysis in this study, based on the single copy orthologous genes of the draft genome, is consistent with traditional phylogenetic classification. The comparative genomic analysis indicated that S. chinensis had 494 species-specific gene families, which involved immune, DNA repair and sensory systems associated with the potential adaption mechanism. We also identified the expansion and positive selection genes in S. chinensis lineage to investigate the potential adaptation mechanism. Our study provided the potential insight into the molecular bases of ecological adaptation in Indo-Pacific humpback dolphin and will be also valuable for future understanding the ecological adaptation and evolution of cetaceans at the genomic level.

Genome Sequence of the Freshwater Yangtze Finless Porpoise

The Yangtze finless porpoise (Neophocaena asiaeorientalis ssp. asiaeorientalis) is a subspecies of the narrow-ridged finless porpoise (N. asiaeorientalis). In total, 714.28 gigabases (Gb) of raw reads were generated by whole-genome sequencing of the Yangtze finless porpoise, using an Illumina HiSeq 2000 platform. After filtering the low-quality and duplicated reads, we assembled a draft genome of 2.22 Gb, with contig N50 and scaffold N50 values of 46.69 kilobases (kb) and 1.71 megabases (Mb), respectively. We identified 887.63 Mb of repetitive sequences and predicted 18,479 protein-coding genes in the assembled genome. The phylogenetic tree showed a relationship between the Yangtze finless porpoise and the Yangtze River dolphin, which diverged approximately 20.84 million years ago. In comparisons with the genomes of 10 other mammals, we detected 44 species-specific gene families, 164 expanded gene families, and 313 positively selected genes in the Yangtze finless porpoise genome. The assembled genome sequence and underlying sequence data are available at the National Center for Biotechnology Information under BioProject accession number PRJNA433603.

Population genomics of finless porpoises reveal an incipient cetacean species adapted to freshwater

Cetaceans (whales, dolphins, and porpoises) are a group of mammals adapted to various aquatic habitats, from oceans to freshwater rivers. We report the sequencing, de novo assembly and analysis of a finless porpoise genome, and the re-sequencing of an additional 48 finless porpoise individuals. We use these data to reconstruct the demographic history of finless porpoises from their origin to the occupation into the Yangtze River. Analyses of selection between marine and freshwater porpoises identify genes associated with renal water homeostasis and urea cycle, such as urea transporter 2 and angiotensin I-converting enzyme 2, which are likely adaptations associated with the difference in osmotic stress between ocean and rivers. Our results strongly suggest that the critically endangered Yangtze finless porpoises are reproductively isolated from other porpoise populations and harbor unique genetic adaptations, supporting that they should be considered a unique incipient species.

Whales, dolphins and porpoises are adapted to various aquatic habitats. Here, Zhou et al. show that polymorphisms associated with renal function and the urea cycle have undergone selection in the freshwater Yangtze finless porpoise and provide genomic evidence of incipient speciation.

Beluga whale pVHL enhances HIF-2α activity via inducing HIF-2α proteasomal degradation under hypoxia

Aquatic mammals, such as cetaceans experience various depths, with accordingly diverse oxygenation, thus, cetaceans have developed adaptations for hypoxia, but mechanisms underlying this tolerance to low oxygen are unclear. Here we analyzed VHL and HIF-2α, in the hypoxia signaling pathway. Variations in VHL are greater than HIF-2α between cetaceans and terrestrial mammals, and beluga whale VHL (BW-VHL) promotes HIF-2α degradation under hypoxia. BW-VHL catalyzes BW-HIF-2α to form K48-linked poly-ubiquitin chains mainly at the lysine 429 of BW-HIF-2α (K429) and induces BW-HIF-2α for proteasomal degradation. W100 within BW-VHL is a key site for BW-VHL functionally and BW-VHL enhances transcriptional activity of BW-HIF-2α under hypoxia. Our data therefore reveal that BW-VHL has a unique function that may contribute to hypoxic adaptation.

Mitochondrial OXPHOS genes provides insights into genetics basis of hypoxia adaptation in anchialine cave shrimps

Cave shrimps from the genera Typhlatya, Stygiocaris and Typhlopatsa (TST complex) comprises twenty cave-adapted taxa, which mainly occur in the anchialine environment. Anchialine habitats may undergo drastic environmental fluctuations, including spatial and temporal changes in salinity, temperature, and dissolved oxygen content. Previous studies of crustaceans from anchialine caves suggest that they have possessed morphological, behavioral, and physiological adaptations to cope with the extreme conditions, similar to other cave-dwelling crustaceans. However, the genetic basis has not been thoroughly explored in crustaceans from anchialine habitats, which can experience hypoxic regimes. To test whether the TST shrimp-complex hypoxia adaptations matched adaptive evolution of mitochondrial OXPHOS genes. The 13 OXPHOS genes from mitochondrial genomes of 98 shrimps and 1 outgroup were examined. For each of these genes was investigated and compared to orthologous sequences using both gene (i.e. branch-site and Datamonkey) and protein (i.e. TreeSAAP) level approaches. Positive selection was detected in 11 of the 13 candidate genes, and the radical amino acid changes sites scattered throughout the entire TST complex phylogeny. Additionally, a series of parallel/convergent amino acid substitutions were identified in mitochondrial OXPHOS genes of TST complex shrimps, which reflect functional convergence or similar genetic mechanisms of cave adaptation. The extensive occurrence of positive selection is suggestive of their essential role in adaptation to hypoxic anchialine environment, and further implying that TST complex shrimps might have acquired a finely capacity for energy metabolism. These results provided some new insights into the genetic basis of anchialine hypoxia adaptation.

Their loss is our gain: regressive evolution in vertebrates provides genomic models for uncovering human disease loci

Throughout Earth's history, evolution's numerous natural 'experiments' have resulted in a diverse range of phenotypes. Though de novo phenotypes receive widespread attention, degeneration of traits inherited from an ancestor is a very common, yet frequently neglected, evolutionary path. The latter phenomenon, known as regressive evolution, often results in vertebrates with phenotypes that mimic inherited disease states in humans. Regressive evolution of anatomical and/or physiological traits is typically accompanied by inactivating mutations underlying these traits, which frequently occur at loci identical to those implicated in human diseases. Here we discuss the potential utility of examining the genomes of vertebrates that have experienced regressive evolution to inform human medical genetics. This approach is low cost and high throughput, giving it the potential to rapidly improve knowledge of disease genetics. We discuss two well-described examples, rod monochromacy (congenital achromatopsia) and amelogenesis imperfecta, to demonstrate the utility of this approach, and then suggest methods to equip non-experts with the ability to corroborate candidate genes and uncover new disease loci.

Inactivation of the olfactory marker protein (OMP) gene in river dolphins and other odontocete cetaceans

Various toothed whales (Odontoceti) are unique among mammals in lacking olfactory bulbs as adults and are thought to be anosmic (lacking the olfactory sense). At the molecular level, toothed whales have high percentages of pseudogenic olfactory receptor genes, but species that have been investigated to date retain an intact copy of the olfactory marker protein gene (OMP), which is highly expressed in olfactory receptor neurons and may regulate the temporal resolution of olfactory responses. One hypothesis for the retention of intact OMP in diverse odontocete lineages is that this gene is pleiotropic with additional functions that are unrelated to olfaction. Recent expression studies provide some support for this hypothesis. Here, we report OMP sequences for representatives of all extant cetacean families and provide the first molecular evidence for inactivation of this gene in vertebrates. Specifically, OMP exhibits independent inactivating mutations in six different odontocete lineages: four river dolphin genera (Platanista, Lipotes, Pontoporia, Inia), sperm whale (Physeter), and harbor porpoise (Phocoena). These results suggest that the only essential role of OMP that is maintained by natural selection is in olfaction, although a non-olfactory role for OMP cannot be ruled out for lineages that retain an intact copy of this gene. Available genome sequences from cetaceans and close outgroups provide evidence of inactivating mutations in two additional genes (CNGA2, CNGA4), which imply further pseudogenization events in the olfactory cascade of odontocetes. Selection analyses demonstrate that evolutionary constraints on all three genes (OMP, CNGA2, CNGA4) have been greatly reduced in Odontoceti, but retain a signature of purifying selection on the stem Cetacea branch and in Mysticeti (baleen whales). This pattern is compatible with the 'echolocation-priority' hypothesis for the evolution of OMP, which posits that negative selection was maintained in the common ancestor of Cetacea and was not relaxed significantly until the evolution of echolocation in Odontoceti.

Analysis of the FGF gene family provides insights into aquatic adaptation in cetaceans

Cetacean body structure and physiology exhibit dramatic adaptations to their aquatic environment. Fibroblast growth factors (FGFs) are a family of essential factors that regulate animal development and physiology; however, their role in cetacean evolution is not clearly understood. Here, we sequenced the fin whale genome and analysed FGFs from 8 cetaceans. FGF22, a hair follicle-enriched gene, exhibited pseudogenization, indicating that the function of this gene is no longer necessary in cetaceans that have lost most of their body hair. An evolutionary analysis revealed signatures of positive selection for FGF3 and FGF11, genes related to ear and tooth development and hypoxia, respectively. We found a D203G substitution in cetacean FGF9, which was predicted to affect FGF9 homodimerization, suggesting that this gene plays a role in the acquisition of rigid flippers for efficient manoeuvring. Cetaceans utilize low bone density as a buoyancy control mechanism, but the underlying genes are not known. We found that the expression of FGF23, a gene associated with reduced bone density, is greatly increased in the cetacean liver under hypoxic conditions, thus implicating FGF23 in low bone density in cetaceans. Altogether, our results provide novel insights into the roles of FGFs in cetacean adaptation to the aquatic environment.

Extreme genomic erosion after recurrent demographic bottlenecks in the highly endangered Iberian lynx

BACKGROUND

Genomic studies of endangered species provide insights into their evolution and demographic history, reveal patterns of genomic erosion that might limit their viability, and offer tools for their effective conservation. The Iberian lynx (Lynx pardinus) is the most endangered felid and a unique example of a species on the brink of extinction.

RESULTS

We generate the first annotated draft of the Iberian lynx genome and carry out genome-based analyses of lynx demography, evolution, and population genetics. We identify a series of severe population bottlenecks in the history of the Iberian lynx that predate its known demographic decline during the 20th century and have greatly impacted its genome evolution. We observe drastically reduced rates of weak-to-strong substitutions associated with GC-biased gene conversion and increased rates of fixation of transposable elements. We also find multiple signatures of genetic erosion in the two remnant Iberian lynx populations, including a high frequency of potentially deleterious variants and substitutions, as well as the lowest genome-wide genetic diversity reported so far in any species.

CONCLUSIONS

The genomic features observed in the Iberian lynx genome may hamper short- and long-term viability through reduced fitness and adaptive potential. The knowledge and resources developed in this study will boost the research on felid evolution and conservation genomics and will benefit the ongoing conservation and management of this emblematic species.

Barcoding the largest animals on Earth: ongoing challenges and molecular solutions in the taxonomic identification of ancient cetaceans

Over the last few centuries, many cetacean species have witnessed dramatic global declines due to industrial overharvesting and other anthropogenic influences, and thus are key targets for conservation. Whale bones recovered from archaeological and palaeontological contexts can provide essential baseline information on the past geographical distribution and abundance of species required for developing informed conservation policies. Here we review the challenges with identifying whale bones through traditional anatomical methods, as well as the opportunities provided by new molecular analyses. Through a case study focused on the North Sea, we demonstrate how the utility of this (pre)historic data is currently limited by a lack of accurate taxonomic information for the majority of ancient cetacean remains. We then discuss current opportunities presented by molecular identification methods such as DNA barcoding and collagen peptide mass fingerprinting (zooarchaeology by mass spectrometry), and highlight the importance of molecular identifications in assessing ancient species' distributions through a case study focused on the Mediterranean. We conclude by considering high-throughput molecular approaches such as hybridization capture followed by next-generation sequencing as cost-effective approaches for enhancing the ecological informativeness of these ancient sample sets.This article is part of the themed issue 'From DNA barcodes to biomes'.

Population Genomics Reveals Low Genetic Diversity and Adaptation to Hypoxia in Snub-Nosed Monkeys

Snub-nosed monkeys (genus Rhinopithecus) are a group of endangered colobines endemic to South Asia. Here, we re-sequenced the whole genomes of 38 snub-nosed monkeys representing four species within this genus. By conducting population genomic analyses, we observed a similar load of deleterious variation in snub-nosed monkeys living in both smaller and larger populations and found that genomic diversity was lower than that reported in other primates. Reconstruction of Rhinopithecus evolutionary history suggested that episodes of climatic variation over the past 2 million years, associated with glacial advances and retreats and population isolation, have shaped snub-nosed monkey demography and evolution. We further identified several hypoxia-related genes under selection in R. bieti (black snub-nosed monkey), a species that exploits habitats higher than any other nonhuman primate. These results provide the first detailed and comprehensive genomic insights into genetic diversity, demography, genetic burden, and adaptation in this radiation of endangered primates.

Genomic Methods Take the Plunge: Recent Advances in High-Throughput Sequencing of Marine Mammals

The dramatic increase in the application of genomic techniques to non-model organisms (NMOs) over the past decade has yielded numerous valuable contributions to evolutionary biology and ecology, many of which would not have been possible with traditional genetic markers. We review this recent progression with a particular focus on genomic studies of marine mammals, a group of taxa that represent key macroevolutionary transitions from terrestrial to marine environments and for which available genomic resources have recently undergone notable rapid growth. Genomic studies of NMOs utilize an expanding range of approaches, including whole genome sequencing, restriction site-associated DNA sequencing, array-based sequencing of single nucleotide polymorphisms and target sequence probes (e.g., exomes), and transcriptome sequencing. These approaches generate different types and quantities of data, and many can be applied with limited or no prior genomic resources, thus overcoming one traditional limitation of research on NMOs. Within marine mammals, such studies have thus far yielded significant contributions to the fields of phylogenomics and comparative genomics, as well as enabled investigations of fitness, demography, and population structure. Here we review the primary options for generating genomic data, introduce several emerging techniques, and discuss the suitability of each approach for different applications in the study of NMOs.