Incomplete lineage sorting and phenotypic evolution in marsupials.

Incomplete lineage sorting (ILS) makes ancestral genetic polymorphisms persist during rapid speciation events, inducing incongruences between gene trees and species trees. ILS has complicated phylogenetic inference in many lineages, including hominids. However, we lack empirical evidence that ILS leads to incongruent phenotypic variation. Here, we performed phylogenomic analyses to show that the South American monito del monte is the sister lineage of all Australian marsupials, although over 31% of its genome is closer to the Diprotodontia than to other Australian groups due to ILS during ancient radiation. Pervasive conflicting phylogenetic signals across the whole genome are consistent with some of the morphological variation among extant marsupials. We detected hundreds of genes that experienced stochastic fixation during ILS, encoding the same amino acids in non-sister species. Using functional experiments, we confirm how ILS may have directly contributed to hemiplasy in morphological traits that were established during rapid marsupial speciation ca. 60 mya.

Comparative genomics uncovers the evolutionary history, demography, and molecular adaptations of South American canids.

The remarkable radiation of South American (SA) canids produced 10 extant species distributed across diverse habitats, including disparate forms such as the short-legged, hypercarnivorous bush dog and the long-legged, largely frugivorous maned wolf. Despite considerable research spanning nearly two centuries, many aspects of their evolutionary history remain unknown. Here, we analyzed 31 whole genomes encompassing all extant SA canid species to assess phylogenetic relationships, interspecific hybridization, historical demography, current genetic diversity, and the molecular bases of adaptations in the bush dog and maned wolf. We found that SA canids originated from a single ancestor that colonized South America 3.9 to 3.5 Mya, followed by diversification east of the Andes and then a single colonization event and radiation of Lycalopex species west of the Andes. We detected extensive historical gene flow between recently diverged lineages and observed distinct patterns of genomic diversity and demographic history in SA canids, likely induced by past climatic cycles compounded by human-induced population declines. Genome-wide scans of selection showed that disparate limb proportions in the bush dog and maned wolf may derive from mutations in genes regulating chondrocyte proliferation and enlargement. Further, frugivory in the maned wolf may have been enabled by variants in genes associated with energy intake from short-chain fatty acids. In contrast, unique genetic variants detected in the bush dog may underlie interdigital webbing and dental adaptations for hypercarnivory. Our analyses shed light on the evolution of a unique carnivoran radiation and how it was shaped by South American topography and climate change.

Standards recommendations for the Earth BioGenome Project.

A global international initiative, such as the Earth BioGenome Project (EBP), requires both agreement and coordination on standards to ensure that the collective effort generates rapid progress toward its goals. To this end, the EBP initiated five technical standards committees comprising volunteer members from the global genomics scientific community: Sample Collection and Processing, Sequencing and Assembly, Annotation, Analysis, and IT and Informatics. The current versions of the resulting standards documents are available on the EBP website, with the recognition that opportunities, technologies, and challenges may improve or change in the future, requiring flexibility for the EBP to meet its goals. Here, we describe some highlights from the proposed standards, and areas where additional challenges will need to be met.

Genomic Signatures of Divergent Ecological Strategies in a Recent Radiation of Neotropical Wild Cats.

Ecological differentiation among diverging species is an important component of the evolutionary process and can be investigated in rapid and recent radiations. Here, we use whole genome sequences of five species from the genus Leopardus, a recently diversified Neotropical lineage with species bearing distinctive morphological, ecological, and behavioral features, to investigate genome-wide diversity, comparative demographic history and signatures of positive selection. Our results show that divergent ecological strategies are reflected in genomic features, for example a generalist species shows historically larger effective population size and higher heterozygosity than habitat specialists. The demographic history of these cats seems to have been jointly driven by climate fluctuations and habitat specialization, with different ecological adaptations leading to distinct trajectories. Finally, a gene involved in vertebrate retinal neurogenesis (POU4F2) was found to be under positive selection in the margay, a cat with notoriously large eyes that are likely associated with its nocturnal and arboreal specializations.

Chromosome-length genome assembly and karyotype of the endangered black-footed ferret (Mustela nigripes).

The black-footed ferret (Mustela nigripes) narrowly avoided extinction to become an oft-cited example of the benefits of intensive management, research, and collaboration to save a species through ex situ conservation breeding and reintroduction into its former range. However, the species remains at risk due to possible inbreeding, disease susceptibility, and multiple fertility challenges. Here, we report the de novo genome assembly of a male black-footed ferret generated through a combination of linked-read sequencing, optical mapping, and Hi-C proximity ligation. In addition, we report the karyotype for this species, which was used to anchor and assign chromosome numbers to the chromosome-length scaffolds. The draft assembly was ~2.5 Gb in length, with 95.6% of it anchored to 19 chromosome-length scaffolds, corresponding to the 2n = 38 chromosomes revealed by the karyotype. The assembly has contig and scaffold N50 values of 148.8 kbp and 145.4 Mbp, respectively, and is up to 96% complete based on BUSCO analyses. Annotation of the assembly, including evidence from RNA-seq data, identified 21,406 protein-coding genes and a repeat content of 37.35%. Phylogenomic analyses indicated that the black-footed ferret diverged from the European polecat/domestic ferret lineage 1.6 million yr ago. This assembly will enable research on the conservation genomics of black-footed ferrets and thereby aid in the further restoration of this endangered species.

The evolutionary history of extinct and living lions.

Lions are one of the world's most iconic megafauna, yet little is known about their temporal and spatial demographic history and population differentiation. We analyzed a genomic dataset of 20 specimens: two ca. 30,000-y-old cave lions (Panthera leo spelaea), 12 historic lions (Panthera leo leo/Panthera leo melanochaita) that lived between the 15th and 20th centuries outside the current geographic distribution of lions, and 6 present-day lions from Africa and India. We found that cave and modern lions shared an ancestor ca. 500,000 y ago and that the 2 lineages likely did not hybridize following their divergence. Within modern lions, we found 2 main lineages that diverged ca. 70,000 y ago, with clear evidence of subsequent gene flow. Our data also reveal a nearly complete absence of genetic diversity within Indian lions, probably due to well-documented extremely low effective population sizes in the recent past. Our results contribute toward the understanding of the evolutionary history of lions and complement conservation efforts to protect the diversity of this vulnerable species.

Evolution of gene regulation in ruminants differs between evolutionary breakpoint regions and homologous synteny blocks.

The role of chromosome rearrangements in driving evolution has been a long-standing question of evolutionary biology. Here we focused on ruminants as a model to assess how rearrangements may have contributed to the evolution of gene regulation. Using reconstructed ancestral karyotypes of Cetartiodactyls, Ruminants, Pecorans, and Bovids, we traced patterns of gross chromosome changes. We found that the lineage leading to the ruminant ancestor after the split from other cetartiodactyls was characterized by mostly intrachromosomal changes, whereas the lineage leading to the pecoran ancestor (including all livestock ruminants) included multiple interchromosomal changes. We observed that the liver cell putative enhancers in the ruminant evolutionary breakpoint regions are highly enriched for DNA sequences under selective constraint acting on lineage-specific transposable elements (TEs) and a set of 25 specific transcription factor (TF) binding motifs associated with recently active TEs. Coupled with gene expression data, we found that genes near ruminant breakpoint regions exhibit more divergent expression profiles among species, particularly in cattle, which is consistent with the phylogenetic origin of these breakpoint regions. This divergence was significantly greater in genes with enhancers that contain at least one of the 25 specific TF binding motifs and located near bovidae-to-cattle lineage breakpoint regions. Taken together, by combining ancestral karyotype reconstructions with analysis of cis regulatory element and gene expression evolution, our work demonstrated that lineage-specific regulatory elements colocalized with gross chromosome rearrangements may have provided valuable functional modifications that helped to shape ruminant evolution.

Earth BioGenome Project: Sequencing life for the future of life.

Increasing our understanding of Earth's biodiversity and responsibly stewarding its resources are among the most crucial scientific and social challenges of the new millennium. These challenges require fundamental new knowledge of the organization, evolution, functions, and interactions among millions of the planet's organisms. Herein, we present a perspective on the Earth BioGenome Project (EBP), a moonshot for biology that aims to sequence, catalog, and characterize the genomes of all of Earth's eukaryotic biodiversity over a period of 10 years. The outcomes of the EBP will inform a broad range of major issues facing humanity, such as the impact of climate change on biodiversity, the conservation of endangered species and ecosystems, and the preservation and enhancement of ecosystem services. We describe hurdles that the project faces, including data-sharing policies that ensure a permanent, freely available resource for future scientific discovery while respecting access and benefit sharing guidelines of the Nagoya Protocol. We also describe scientific and organizational challenges in executing such an ambitious project, and the structure proposed to achieve the project's goals. The far-reaching potential benefits of creating an open digital repository of genomic information for life on Earth can be realized only by a coordinated international effort.

Pangolin genomes and the evolution of mammalian scales and immunity.

Pangolins, unique mammals with scales over most of their body, no teeth, poor vision, and an acute olfactory system, comprise the only placental order (Pholidota) without a whole-genome map. To investigate pangolin biology and evolution, we developed genome assemblies of the Malayan (Manis javanica) and Chinese (M. pentadactyla) pangolins. Strikingly, we found that interferon epsilon (IFNE), exclusively expressed in epithelial cells and important in skin and mucosal immunity, is pseudogenized in all African and Asian pangolin species that we examined, perhaps impacting resistance to infection. We propose that scale development was an innovation that provided protection against injuries or stress and reduced pangolin vulnerability to infection. Further evidence of specialized adaptations was evident from positively selected genes involving immunity-related pathways, inflammation, energy storage and metabolism, muscular and nervous systems, and scale/hair development. Olfactory receptor gene families are significantly expanded in pangolins, reflecting their well-developed olfaction system. This study provides insights into mammalian adaptation and functional diversification, new research tools and questions, and perhaps a new natural IFNE-deficient animal model for studying mammalian immunity.

Adaptive genomic evolution of opsins reveals that early mammals flourished in nocturnal environments.

BACKGROUND: Based on evolutionary patterns of the vertebrate eye, Walls (1942) hypothesized that early placental mammals evolved primarily in nocturnal habitats. However, not only Eutheria, but all mammals show photic characteristics (i.e. dichromatic vision, rod-dominated retina) suggestive of a scotopic eye design. RESULTS: Here, we used integrative comparative genomic and phylogenetic methodologies employing the photoreceptive opsin gene family in 154 mammals to test the likelihood of a nocturnal period in the emergence of all mammals. We showed that mammals possess genomic patterns concordant with a nocturnal ancestry. The loss of the RH2, VA, PARA, PARIE and OPN4x opsins in all mammals led us to advance a probable and most-parsimonious hypothesis of a global nocturnal bottleneck that explains the loss of these genes in the emerging lineage (> > 215.5 million years ago). In addition, ancestral character reconstruction analyses provided strong evidence that ancestral mammals possessed a nocturnal lifestyle, ultra-violet-sensitive vision, low visual acuity and low orbit convergence (i.e. panoramic vision). CONCLUSIONS: Overall, this study provides insight into the evolutionary history of the mammalian eye while discussing important ecological aspects of the photic paleo-environments ancestral mammals have occupied.

Landscape genomics: natural selection drives the evolution of mitogenome in penguins.

BACKGROUND: Mitochondria play a key role in the balance of energy and heat production, and therefore the mitochondrial genome is under natural selection by environmental temperature and food availability, since starvation can generate more efficient coupling of energy production. However, selection over mitochondrial DNA (mtDNA) genes has usually been evaluated at the population level. We sequenced by NGS 12 mitogenomes and with four published genomes, assessed genetic variation in ten penguin species distributed from the equator to Antarctica. Signatures of selection of 13 mitochondrial protein-coding genes were evaluated by comparing among species within and among genera (Spheniscus, Pygoscelis, Eudyptula, Eudyptes and Aptenodytes). The genetic data were correlated with environmental data obtained through remote sensing (sea surface temperature [SST], chlorophyll levels [Chl] and a combination of SST and Chl [COM]) through the distribution of these species. RESULTS: We identified the complete mtDNA genomes of several penguin species, including ND6 and 8 tRNAs on the light strand and 12 protein coding genes, 14 tRNAs and two rRNAs positioned on the heavy strand. The highest diversity was found in NADH dehydrogenase genes and the lowest in COX genes. The lowest evolutionary divergence among species was between Humboldt (Spheniscus humboldti) and Galapagos (S. mendiculus) penguins (0.004), while the highest was observed between little penguin (Eudyptula minor) and Adélie penguin (Pygoscelis adeliae) (0.097). We identified a signature of purifying selection (Ka/Ks < 1) across the mitochondrial genome, which is consistent with the hypothesis that purifying selection is constraining mitogenome evolution to maintain Oxidative phosphorylation (OXPHOS) proteins and functionality. Pairwise species maximum-likelihood analyses of selection at codon sites suggest positive selection has occurred on ATP8 (Fixed-Effects Likelihood, FEL) and ND4 (Single Likelihood Ancestral Counting, SLAC) in all penguins. In contrast, COX1 had a signature of strong negative selection. ND4 Ka/Ks ratios were highly correlated with SST (Mantel, p-value: 0.0001; GLM, p-value: 0.00001) and thus may be related to climate adaptation throughout penguin speciation. CONCLUSIONS: These results identify mtDNA candidate genes under selection which could be involved in broad-scale adaptations of penguins to their environment. Such knowledge may be particularly useful for developing predictive models of how these species may respond to severe climatic changes in the future.

The Complete Phylogeny of Pangolins: Scaling Up Resources for the Molecular Tracing of the Most Trafficked Mammals on Earth.

Pangolins, considered the most-trafficked mammals on Earth, are rapidly heading to extinction. Eight extant species of these African and Asian scale-bodied anteaters are commonly recognized, but their evolutionary relationships remain largely unexplored. Here, we present the most comprehensive phylogenetic assessment of pangolins, based on genetic variation of complete mitogenomes and 9 nuclear genes. We confirm deep divergence among Asian and African pangolins occurring not later than the Oligocene-Miocene boundary ca. 23 million years ago (Ma) (95% HPD = 18.7-27.2), limited fossil evidence suggesting dispersals from Europe. We recognize 3 genera including Manis (Asian pangolins), Smutsia (large African pangolins), and Phataginus (small African pangolins), which first diversified in the Middle-Upper Miocene (9.8-13.3 Ma) through a period of gradual cooling coinciding with a worldwide taxonomic diversification among mammals. Based on large mitogenomic distances among the 3 genera (18.3-22.8%) and numerous (18) morphological traits unique to Phataginus, we propose the subfamily Phatagininae subfam. nov. to designate small African pangolins. In contrast with the morphological-based literature, our results establish that the thick-tailed pangolin (Manis crassicaudata) is sister-species of the Sunda (Manis javanica) and Palawan (Manis culionensis) pangolins. Mitogenomic phylogenetic delineations supported additional pangolin species subdivisions (n = 13), including 6 African common pangolin (Phataginus tricuspis) lineages, but these patterns were not fully supported by our multi-locus approach. Finally, we identified more than 5000 informative mitogenomic sites and diagnostic variation from 5 nuclear genes among all species and lineages of pangolins, providing an important resource for further research and for effectively tracing the worldwide pangolin trade.

Conservation Genetics of the Cheetah: Lessons Learned and New Opportunities.

The dwindling wildlife species of our planet have become a cause célèbre for conservation groups, governments, and concerned citizens throughout the world. The application of powerful new genetic technologies to surviving populations of threatened mammals has revolutionized our ability to recognize hidden perils that afflict them. We have learned new lessons of survival, adaptation, and evolution from viewing the natural history of genomes in hundreds of detailed studies. A single case history of one species, the African cheetah, Acinonyx jubatus, is here reviewed to reveal a long-term story of conservation challenges and action informed by genetic discoveries and insights. A synthesis of 3 decades of data, interpretation, and controversy, capped by whole genome sequence analysis of cheetahs, provides a compelling tale of conservation relevance and action to protect this species and other threatened wildlife.

Bone-associated gene evolution and the origin of flight in birds.

BACKGROUND: Bones have been subjected to considerable selective pressure throughout vertebrate evolution, such as occurred during the adaptations associated with the development of powered flight. Powered flight evolved independently in two extant clades of vertebrates, birds and bats. While this trait provided advantages such as in aerial foraging habits, escape from predators or long-distance travels, it also imposed great challenges, namely in the bone structure. RESULTS: We performed comparative genomic analyses of 89 bone-associated genes from 47 avian genomes (including 45 new), 39 mammalian, and 20 reptilian genomes, and demonstrate that birds, after correcting for multiple testing, have an almost two-fold increase in the number of bone-associated genes with evidence of positive selection (~52.8 %) compared with mammals (~30.3 %). Most of the positive-selected genes in birds are linked with bone regulation and remodeling and thirteen have been linked with functional pathways relevant to powered flight, including bone metabolism, bone fusion, muscle development and hyperglycemia levels. Genes encoding proteins involved in bone resorption, such as TPP1, had a high number of sites under Darwinian selection in birds. CONCLUSIONS: Patterns of positive selection observed in bird ossification genes suggest that there was a period of intense selective pressure to improve flight efficiency that was closely linked with constraints on body size.

Olfactory Receptor Subgenomes Linked with Broad Ecological Adaptations in Sauropsida.

Olfactory receptors (ORs) govern a prime sensory function. Extant birds have distinct olfactory abilities, but the molecular mechanisms underlining diversification and specialization remain mostly unknown. We explored OR diversity in 48 phylogenetic and ecologically diverse birds and 2 reptiles (alligator and green sea turtle). OR subgenomes showed species- and lineage-specific variation related with ecological requirements. Overall 1,953 OR genes were identified in reptiles and 16,503 in birds. The two reptiles had larger OR gene repertoires (989 and 964 genes, respectively) than birds (182-688 genes). Overall, birds had more pseudogenes (7,855) than intact genes (1,944). The alligator had significantly more functional genes than sea turtle, likely because of distinct foraging habits. We found rapid species-specific expansion and positive selection in OR14 (detects hydrophobic compounds) in birds and in OR51 and OR52 (detect hydrophilic compounds) in sea turtle, suggestive of terrestrial and aquatic adaptations, respectively. Ecological partitioning among birds of prey, water birds, land birds, and vocal learners showed that diverse ecological factors determined olfactory ability and influenced corresponding olfactory-receptor subgenome. OR5/8/9 was expanded in predatory birds and alligator, suggesting adaptive specialization for carnivory. OR families 2/13, 51, and 52 were correlated with aquatic adaptations (water birds), OR families 6 and 10 were more pronounced in vocal-learning birds, whereas most specialized land birds had an expanded OR family 14. Olfactory bulb ratio (OBR) and OR gene repertoire were correlated. Birds that forage for prey (carnivores/piscivores) had relatively complex OBR and OR gene repertoires compared with modern birds, including passerines, perhaps due to highly developed cognitive capacities facilitating foraging innovations.

Molecular assessment of the phylogeny and biogeography of a recently diversified endemic group of South American canids (Mammalia: Carnivora: Canidae).

To investigate the evolution and biogeography of an endemic group of South American foxes, we examined mitochondrial DNA control region sequences for 118 individuals belonging to all six extant species of the genus Lycalopex. Phylogenetic and molecular dating analyses supported the inference that this genus has undergone a very recent and rapid radiation, stemming from a common ancestor that lived ca. 1 million years ago. The Brazilian endemic L. vetulus was supported as the most basal species in this genus, whereas the most internal group is comprised by the recently diverged (ca. 350,000 years ago) Andean/Patagonian species L. griseus and L. culpaeus. We discuss the inferred phylogenetic relationships and divergence times in the context of the current geographic distributions of these species, and the likely effects of Pleistocene climatic changes on the biogeography of this group. Furthermore, a remarkable finding was the identification of multiple individuals classified as L. gymnocercus bearing mtDNA haplotypes clearly belonging to L. griseus, sampled in regions where the latter is not known to occur. At a minimum, this result implies the need to clarify the present-day geographic distribution of each of these fox species, while it may also indicate an ongoing hybridization process between them. Future testing of this hypothesis with in-depth analyses of these populations is thus a priority for understanding the history, evolutionary dynamics and present-day composition of this endemic Neotropical genus.

Gene loss, adaptive evolution and the co-evolution of plumage coloration genes with opsins in birds.

BACKGROUND: The wide range of complex photic systems observed in birds exemplifies one of their key evolutionary adaptions, a well-developed visual system. However, genomic approaches have yet to be used to disentangle the evolutionary mechanisms that govern evolution of avian visual systems. RESULTS: We performed comparative genomic analyses across 48 avian genomes that span extant bird phylogenetic diversity to assess evolutionary changes in the 17 representatives of the opsin gene family and five plumage coloration genes. Our analyses suggest modern birds have maintained a repertoire of up to 15 opsins. Synteny analyses indicate that PARA and PARIE pineal opsins were lost, probably in conjunction with the degeneration of the parietal organ. Eleven of the 15 avian opsins evolved in a non-neutral pattern, confirming the adaptive importance of vision in birds. Visual conopsins sw1, sw2 and lw evolved under negative selection, while the dim-light RH1 photopigment diversified. The evolutionary patterns of sw1 and of violet/ultraviolet sensitivity in birds suggest that avian ancestors had violet-sensitive vision. Additionally, we demonstrate an adaptive association between the RH2 opsin and the MC1R plumage color gene, suggesting that plumage coloration has been photic mediated. At the intra-avian level we observed some unique adaptive patterns. For example, barn owl showed early signs of pseudogenization in RH2, perhaps in response to nocturnal behavior, and penguins had amino acid deletions in RH2 sites responsible for the red shift and retinal binding. These patterns in the barn owl and penguins were convergent with adaptive strategies in nocturnal and aquatic mammals, respectively. CONCLUSIONS: We conclude that birds have evolved diverse opsin adaptations through gene loss, adaptive selection and coevolution with plumage coloration, and that differentiated selective patterns at the species level suggest novel photic pressures to influence evolutionary patterns of more-recent lineages.

Genetic diversity of Rhizobium from nodulating beans grown in a variety of Mediterranean climate soils of Chile.

In spite of potentially being an important source of rhizobial diversity and a key determinant of common bean productivity, there is a paucity of data on Rhizobium genetic variation and species composition in the important bean producing area of Chile and only one species has been documented (Rhizobium leguminosarum). In this study, 240 Rhizobium isolates from Torcaza bean (Phaseolus vulgaris L.) nodules established in the highest bean producing area in Chile (33°34'S-70°38'W and 37°36'S-71°47'W) were characterized by PCR-RFLP markers for nodC gene, revealing eight banding patterns with the polymorphic enzyme Hinf I. The locality of San Agustín de Aurora in Central Chile (35°32'S-71°29'W) had the highest level of diversity. Isolates were classified by species using PCR-RFLP markers for 16S rDNA gene and were confirmed by sequencing an internal fragment of the 16S rDNA gene. The results confirmed the presence of R. leguminosarum and three other species of rhizobia nodulating beans in South Central Chile (R. etli, R. tropici and R. leucaenae). R. tropici and R. leucaenae showed the least genetic variation and were most commonly identified in acid soils, while R. etli was the most common species in slightly acidic to moderately alkaline soils, with higher levels of organic matter content. R. leguminosarum was identified in almost all soils, was the most genetically diverse, and was the most common, being documented in soils with pH that ranged between 5.3 and 8.2, and with organic matter content between 2.1 and 4 %.

Reduced Genetic Diversity and Increased Dispersal in Guigna (Leopardus guigna) in Chilean Fragmented Landscapes.

Landscape fragmentation is often a major cause of species extinction as it can affect a wide variety of ecological processes. The impact of fragmentation varies among species depending on many factors, including their life-history traits and dispersal abilities. Felids are one of the groups most threatened by fragmented landscapes because of their large home ranges, territorial behavior, and low population densities. Here, we model the impacts of habitat fragmentation on patterns of genetic diversity in the guigna (Leopardus guigna), a small felid that is closely associated with the heavily human-impacted temperate rainforests of southern South America. We assessed genetic variation in 1798 base pairs of mitochondrial DNA sequences, 15 microsatellite loci, and 2 sex chromosome genes and estimated genetic diversity, kinship, inbreeding, and dispersal in 38 individuals from landscapes with differing degrees of fragmentation on Chiloé Island in southern Chile. Increased fragmentation was associated with reduced genetic diversity, but not with increased kinship or inbreeding. However, in fragmented landscapes, there was a weaker negative correlation between pairwise kinship and geographic distance, suggesting increased dispersal distances. These results highlight the importance of biological corridors to maximize connectivity in fragmented landscapes and contribute to our understanding of the broader genetic consequences of habitat fragmentation, especially for forest-specialist carnivores.