Search

Genome assemblies are growing at an exponential rate and have proved indispensable for studying evolution but the effort has been biased toward vertebrates and arthropods with a particular focus on insects. Onychophora or velvet worms are an ancient group of cryptic, soil dwelling worms noted for their unique mode of prey capture, biogeographic patterns, and diversity of reproductive strategies. They constitute a poorly understood phylum of exclusively terrestrial animals that is sister group to arthropods. Due to this phylogenetic position, they are crucial in understanding the origin of the largest phylum of animals. Despite their significance, there is a paucity of genomic resources for the phylum with only one highly fragmented and incomplete genome publicly available. Initial attempts at sequencing an onychophoran genome proved difficult due to its large genome size and high repeat content. However, leveraging recent advances in long-read sequencing technology, we present here the first annotated draft genome for the phylum. With a total size of 5.6Gb, the gigantism of the Epiperipatus broadwayi genome arises from having high repeat content, intron size inflation, and extensive gene family expansion. Additionally, we report a previously unknown diversity of onychophoran hemocyanins that suggests the diversification of copper-mediated oxygen carriers occurred independently in Onychophora after its split from Arthropoda, parallel to the independent diversification of hemocyanins in each of the main arthropod lineages.

From savory to sweet Seeing a bird eat nectar from a flower is a common sight in our world. The ability to detect sugars, however, is not ancestral in the bird lineage, where most species were carnivorous. Toda et al. looked at receptors within the largest group of birds, the passerines or songbirds, and found that the emergence of sweet detection involved a single shift in a receptor for umami (see the Perspective by Barker). This ancient change facilitated sugar detection not just in nectar feeding birds, but also across the songbird group, and in a way that was different from, though convergent with, that in hummingbirds. Science , abf6505, this issue p. 226 ; see also abj6746, p. 154

Occurring independently at seven separate origins across the avian tree of life, obligate brood parasitism is a unique suite of traits observed in only approximately 1% of all bird species. Obligate brood parasites exhibit varied physiological, morphological, and behavioural traits across lineages, but common among all obligate brood parasites is that the females lay their eggs in the nest of other species. Unique among these species is the black-headed duck (Heteronetta atricapilla), a generalist brood parasite that is the only obligate brood parasite among waterfowl. This provides an opportunity to assess evolutionary changes in traits associated with brood parasitism, notably the loss of parental care behaviours, with an unspecialized brood parasite. We generated new high-quality genome assemblies and genome annotations of the black-headed duck and three related non-parasitic species (freckled duck, African pygmy-goose, and ruddy duck). With these assemblies and existing public genome assemblies, we produced a whole genome alignment across Galloanserae to identify conserved non-coding regions with atypical accelerations in the black-headed duck and coding genes with evidence of positive selection, as well as to resolve uncertainties in the duck phylogeny. To complement these data, we sequenced a population sample of black-headed ducks, allowing us to conduct McDonald-Kreitman tests of lineage-specific selection. We resolve the existing polytomy between our focal taxa with concordance from coding and non-coding sequences, and we observe stronger signals of evolution in non-coding regions of the genome than in coding regions. Collectively, the new high-quality genomes, comparative genome alignment, and population genomics provide a detailed picture of genome evolution in the only brood parasitic duck.

A major goal of comparative genomics research is modeling changes in DNA sequences between species to understand the evolutionary forces acting on species differences. Application of these models to a number of species over the past decade has revealed some commonalities across organisms, most notably a consistent role of positive selection in shaping the molecular evolution of the immune system. However, models of DNA sequence evolution also have important limitations that are increasingly being recognized, including issues with data quality and biases caused by simplifying assumptions. While new approaches have begun to address these challenges, ultimately, additional data, such as resequencing data from populations, will provide more power to fully understand the unique evolutionary forces acting on different species. In this review, I summarize the conclusions of recent genome-wide studies of selection, highlight some important challenges to applying these methods to large data sets, and discuss ways forward for the field.

The standard neutral model of molecular evolution has traditionally been used as the null model for population genomics. We gathered a collection of 45 genome-wide site frequency spectra from a diverse set of species, most of which display an excess of low and high frequency variants compared to the expectation of the standard neutral model, resulting in U-shaped spectra. We show that multiple merger coalescent models often provide a better fit to these observations than the standard Kingman coalescent. Hence, in many circumstances these under-utilized models may serve as the more appropriate reference for genomic analyses. We further discuss the underlying evolutionary processes that may result in the widespread U-shape of frequency spectra.

Nearly all eukaryotes reproduce sexually, either constitutively or facultatively, and nearly all that are thought to be entirely asexual diverged recently from sexuals, implying that loss of sex leads to early extinction and therefore that sexual reproduction is essential for evolutionary success. Nevertheless, there are several groups that have been thought to have evolved asexually for millions of years. Of these, the most extensively studied are the rotifers of Class Bdelloidea. Yet the evidence for their asexuality is entirely negative -- the failure to establish the existence of males or hermaphrodites. Opposed to this is a growing body of evidence that bdelloids do reproduce sexually, albeit rarely, retaining meiosis-associated genes and, in a limited study of allele sharing among individuals of the bdelloid Macrotrachela quadricornifera, displaying a pattern of genetic exchange indicating recent sexual reproduction. Here we present a much larger study of allele sharing among these individuals clearly showing the occurrence of sexual reproduction, thereby removing what had appeared to be a serious challenge to hypotheses for the evolutionary advantage of sex and confirming that sexual reproduction is essential for long term evolutionary success in eukaryotes. We also discuss the relation between bdelloid life history and population structure and the possible benefit of outcrossing in restoring advantageous genome-wide epistatic interactions disrupted by conversion and deletion.

Palaeognathae represent one of the two basal lineages in modern birds, and comprise the volant (flighted) tinamous and the flightless ratites. Resolving palaeognath phylogenetic relationships has historically proved difficult, and short internal branches separating major palaeognath lineages in previous molecular phylogenies suggest that extensive incomplete lineage sorting (ILS) might have accompanied a rapid ancient divergence. Here, we investigate palaeognath relationships using genome-wide data sets of three types of noncoding nuclear markers, together totaling 20,850 loci and over 41 million base pairs of aligned sequence data. We recover a fully resolved topology placing rheas as the sister to kiwi and emu + cassowary that is congruent across marker types for two species tree methods (MP-EST and ASTRAL-II). This topology is corroborated by patterns of insertions for 4274 CR1 retroelements identified from multispecies whole-genome screening, and is robustly supported by phylogenomic subsampling analyses, with MP-EST demonstrating particularly consistent performance across subsampling replicates as compared to ASTRAL. In contrast, analyses of concatenated data supermatrices recover rheas as the sister to all other nonostrich palaeognaths, an alternative that lacks retroelement support and shows inconsistent behavior under subsampling approaches. While statistically supporting the species tree topology, conflicting patterns of retroelement insertions also occur and imply high amounts of ILS across short successive internal branches, consistent with observed patterns of gene tree heterogeneity. Coalescent simulations and topology tests indicate that the majority of observed topological incongruence among gene trees is consistent with coalescent variation rather than arising from gene tree estimation error alone, and estimated branch lengths for short successive internodes in the inferred species tree fall within the theoretical range encompassing the anomaly zone. Distributions of empirical gene trees confirm that the most common gene tree topology for each marker type differs from the species tree, signifying the existence of an empirical anomaly zone in palaeognaths.

Classical evolutionary studies of protein-coding genes have established that genes in the canonical immune system are often among the most rapidly evolving within and between species. As more genomes and transcriptomes across insects are sequenced, it is becoming clear that duplications and losses of immune genes are also a likely consequence of host-pathogen interactions. Furthermore, particular species respond to diverse pathogenic challenges with a wide range of challenge-specific responses that are still poorly understood. Transcriptional studies, using RNA-seq to characterize the infection-regulated transcriptome of diverse insects, are crucial for additional progress in understanding the ecology and evolution of the full complexity of the host response.

Background The stable fly, Stomoxys calcitrans, is a major blood-feeding pest of livestock that has near worldwide distribution, causing an annual cost of over $2 billion for control and product loss in the USA alone. Control of these flies has been limited to increased sanitary management practices and insecticide application for suppressing larval stages. Few genetic and molecular resources are available to help in developing novel methods for controlling stable flies. Results This study examines stable fly biology by utilizing a combination of high-quality genome sequencing and RNA-Seq analyses targeting multiple developmental stages and tissues. In conjunction, 1600 genes were manually curated to characterize genetic features related to stable fly reproduction, vector host interactions, host-microbe dynamics, and putative targets for control. Most notable was characterization of genes associated with reproduction and identification of expanded gene families with functional associations to vision, chemosensation, immunity, and metabolic detoxification pathways. Conclusions The combined sequencing, assembly, and curation of the male stable fly genome followed by RNA-Seq and downstream analyses provide insights necessary to understand the biology of this important pest. These resources and new data will provide the groundwork for expanding the tools available to control stable fly infestations. The close relationship of Stomoxys to other blood-feeding (horn flies and Glossina) and non-blood-feeding flies (house flies, medflies, Drosophila) will facilitate understanding of the evolutionary processes associated with development of blood feeding among the Cyclorrhapha.

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity1–4. Sparse taxon sampling has previously been proposed to confound phylogenetic inference5, and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families—including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species., A dataset of the genomes of 363 species from the Bird 10,000 Genomes Project shows increased power to detect shared and lineage-specific variation, demonstrating the importance of phylogenetically diverse taxon sampling in whole-genome sequencing.

Despite more than 2,000-fold variation in genome size, key features of genome architecture are largely conserved across angiosperms. Parasitic plants have elucidated the many ways in which genomes can be modified, yet we still lack comprehensive genome data for species that represent the most extreme form of parasitism. Here, we present the highly modified genome of the iconic endophytic parasite Sapria himalayana Griff. (Rafflesiaceae), which lacks a typical plant body. First, 44% of the genes conserved in eurosids are lost in Sapria, dwarfing previously reported levels of gene loss in vascular plants. These losses demonstrate remarkable functional convergence with other parasitic plants, suggesting a common genetic roadmap underlying the evolution of plant parasitism. Second, we identified extreme disparity in intron size among retained genes. This includes a category of genes with introns longer than any so far observed in angiosperms, nearing 100 kb in some cases, and a second category of genes with exceptionally short or absent introns. Finally, at least 1.2% of the Sapria genome, including both genic and intergenic content, is inferred to be derived from host-to-parasite horizontal gene transfers (HGTs) and includes genes potentially adaptive for parasitism. Focused phylogenomic reconstruction of HGTs reveals a hidden history of former host-parasite associations involving close relatives of Sapria's modern hosts in the grapevine family. Our findings offer a unique perspective into how deeply angiosperm genomes can be altered to fit an extreme form of plant parasitism and demonstrate the value of HGTs as DNA fossils to investigate extinct symbioses.

The emergence of single-cell RNA sequencing (scRNA-seq) has led to an explosion in novel methods to study biological variation among individual cells, and to classify cells into functional and biologically meaningful categories. Here, we present a new cell type projection tool, HieRFIT (Hierarchical Random Forest for Information Transfer), based on hierarchical random forests. HieRFIT uses a priori information about cell type relationships to improve classification accuracy, taking as input a hierarchical tree structure representing the class relationships, along with the reference data. We use an ensemble approach combining multiple random forest models, organized in a hierarchical decision tree structure. We show that our hierarchical classification approach improves accuracy and reduces incorrect predictions especially for inter-dataset tasks which reflect real life applications. We use a scoring scheme that adjusts probability distributions for candidate class labels and resolves uncertainties while avoiding the assignment of cells to incorrect types by labeling cells at internal nodes of the hierarchy when necessary. Using HieRFIT, we re-analyzed publicly available scRNA-seq datasets showing its effectiveness in cell type cross-projections with inter/intra-species examples. HieRFIT is implemented as an R package and it is available at (https://github.com/yasinkaymaz/HieRFIT/releases/tag/v1.0.0)

Motivation One major goal of single-cell RNA sequencing (scRNAseq) experiments is to identify novel cell types. With increasingly large scRNAseq datasets, unsupervised clustering methods can now produce detailed catalogues of transcriptionally distinct groups of cells in a sample. However, the interpretation of these clusters is challenging for both technical and biological reasons. Popular clustering algorithms are sensitive to parameter choices, and can produce different clustering solutions with even small changes in the number of principal components used, the k nearest neighbor and the resolution parameters, among others. Results Here, we present a set of tools to evaluate cluster stability by subsampling, which can guide parameter choice and aid in biological interpretation. The R package scclusteval and the accompanying Snakemake workflow implement all steps of the pipeline: subsampling the cells, repeating the clustering with Seurat and estimation of cluster stability using the Jaccard similarity index and providing rich visualizations. Availabilityand implementation R package scclusteval: https://github.com/crazyhottommy/scclusteval Snakemake workflow: https://github.com/crazyhottommy/pyflow_seuratv3_parameter Tutorial: https://crazyhottommy.github.io/EvaluateSingleCellClustering/.

MotivationOne major goal of single-cell RNA sequencing (scRNAseq) experiments is to identify novel cell types. With increasingly large scRNAseq datasets, unsupervised clustering methods can now produce detailed catalogues of transcriptionally distinct groups of cells in a sample. However, the interpretation of these clusters is challenging for both technical and biological reasons. Popular clustering algorithms are sensitive to parameter choices, and can produce different clustering solutions with even small changes in the number of principal components used, the k nearest neighbor, and the resolution parameters, among others.ResultsHere, we present a set of tools to evaluate cluster stability by subsampling, which can guide parameter choice and aid in biological interpretation. The R package scclusteval and the accompanying Snakemake workflow implement all steps of the pipeline: subsampling the cells, repeating the clustering with Seurat, and estimation of cluster stability using the Jaccard similarity index. The Snakemake workflow takes advantage of high-performance computing clusters and dispatches jobs in parallel to available CPUs to speed up the analysis. The scclusteval package provides functions to facilitate the analysis of the output, including a series of rich visualizations.AvailabilityR package scclusteval: https://github.com/crazyhottommy/scclusteval Snakemake workflow: https://github.com/crazyhottommy/pyflow_seuratv3_parameterContacttsackton@g.harvard.edu, tangming2005@gmail.comSupplementary informationSupplementary data are available at Bioinformatics online.

The exact genetic mechanism regulating avian gonadal sex differentiation has not been completely resolved. The most likely scenario involves a dosage mechanism, whereby the Z-linked DMRT1 gene triggers testis development. However, the possibility still exists that the female-specific W chromosome may harbor an ovarian determining factor. In this study, we provide evidence that the universal gene regulating gonadal sex differentiation in birds is Z-linked DMRT1 and not a W-linked (ovarian) factor. Three candidate W-linked ovarian determinants are HINTW, female-expressed transcript 1 (FET1), and female-associated factor (FAF). To test the association of these genes with ovarian differentiation in the chicken, we examined their expression following experimentally induced female-to-male sex reversal using the aromatase inhibitor fadrozole (FAD). Administration of FAD on day 3 of embryogenesis induced a significant loss of aromatase enzyme activity in female gonads and masculinization. However, expression levels of HINTW, FAF, and FET1 were unaltered after experimental masculinization. Furthermore, comparative analysis showed that FAF and FET1 expression could not be detected in zebra finch gonads. Additionally, an antibody raised against the predicted HINTW protein failed to detect it endogenously. These data do not support a universal role for these genes or for the W sex chromosome in ovarian development in birds. We found that DMRT1 (but not the recently identified Z-linked HEMGN gene) is male upregulated in embryonic zebra finch and emu gonads, as in the chicken. As chicken, zebra finch, and emu exemplify the major evolutionary clades of birds, we propose that Z-linked DMRT1, and not the W sex chromosome, regulates gonadal sex differentiation in birds.

Evolutionary convergence has been long considered primary evidence of adaptation driven by natural selection and provides opportunities to explore evolutionary repeatability and predictability. In recent years, there has been increased interest in exploring the genetic mechanisms underlying convergent evolution, in part, owing to the advent of genomic techniques. However, the current ‘genomics gold rush’ in studies of convergence has overshadowed the reality that most trait classifications are quite broadly defined, resulting in incomplete or potentially biased interpretations of results. Genomic studies of convergence would be greatly improved by integrating deep ‘vertical’, natural history knowledge with ‘horizontal’ knowledge focusing on the breadth of taxonomic diversity. Natural history collections have and continue to be best positioned for increasing our comprehensive understanding of phenotypic diversity, with modern practices of digitization and databasing of morphological traits providing exciting improvements in our ability to evaluate the degree of morphological convergence. Combining more detailed phenotypic data with the well-established field of genomics will enable scientists to make progress on an important goal in biology: to understand the degree to which genetic or molecular convergence is associated with phenotypic convergence. Although the fields of comparative biology or comparative genomics alone can separately reveal important insights into convergent evolution, here we suggest that the synergistic and complementary roles of natural history collection-derived phenomic data and comparative genomics methods can be particularly powerful in together elucidating the genomic basis of convergent evolution among higher taxa. This article is part of the theme issue ‘Convergent evolution in the genomics era: new insights and directions’.

BackgroundThe stable fly,Stomoxys calcitrans, is a major blood-feeding pest of livestock that has near worldwide distribution, causing an annual cost of over $2 billion for control and product loss in the United States alone. Control of these flies has been limited to increased sanitary management practices and insecticide application for suppressing larval stages. Few genetic and molecular resources are available to help in developing novel methods for controlling stable flies.ResultsThis study examines stable fly biology by utilizing a combination of high-quality genome sequencing, microbiome analyses, and RNA-seq analyses targeting multiple developmental stages and tissues. In conjunction, manual curation of over 1600 genes was used to examine gene content related to stable fly reproduction, interactions with their host, host-microbe dynamics, and putative routes for control. Most notable was establishment of reproduction-associated genes and identification of expanded vision, chemosensation, immune repertoire, and metabolic detoxification pathway gene families.ConclusionsThe combined sequencing, assembly, and curation of the male stable fly genome followed by RNA-seq and downstream analyses provide insights necessary to understand the biology of this important pest. These resources and knowledge will provide the groundwork for expanding the tools available to control stable fly infestations. The close relationship ofStomoxysto other blood-feeding (Glossina) and non-blood-feeding flies (medflies,Drosophila, house flies) will allow for understanding the evolution of blood feeding among Cyclorrhapha flies.

De novo transcriptome assembly is a powerful tool, widely used over the last decade for making evolutionary inferences. However, it relies on two implicit assumptions: that the assembled transcriptome is an unbiased representation of the underlying expressed transcriptome, and that expression estimates from the assembly are good, if noisy approximations of the relative abundance of expressed transcripts. Using publicly available data for model organisms, we demonstrate that, across assembly algorithms and data sets, these assumptions are consistently violated. Bias exists at the nucleotide level, with genotyping error rates ranging from 30-83%. As a result, diversity is underestimated in transcriptome assemblies, with consistent under-estimation of heterozygosity in all but the most inbred samples. Even at the gene level, expression estimates show wide deviations from map-to-reference estimates, and positive bias at lower expression levels. Standard filtering of transcriptome assemblies improves the robustness of gene expression estimates but leads to the loss of a meaningful number of protein-coding genes, including many that are highly expressed. We demonstrate a computational method, length-rescaled CPM, to partly alleviate noise and bias in expression estimates. Researchers should consider ways to minimize the impact of bias in transcriptome assemblies.

Evolutionary convergence has been long considered primary evidence of adaptation driven by natural selection and provides opportunities to explore evolutionary repeatability and predictability. In recent years, there has been increased interest in exploring the genetic mechanisms underlying convergent evolution, in part due to the advent of genomic techniques. However, the current ‘genomics gold rush’ in studies of convergence has overshadowed the reality that most trait classifications are quite broadly defined, resulting in incomplete or potentially biased interpretations of results. Genomic studies of convergence would be greatly improved by integrating deep ‘vertical’, natural history knowledge with ‘horizontal’ knowledge focusing on the breadth of taxonomic diversity. Natural history collections have and continue to be best positioned for increasing our comprehensive understanding of phenotypic diversity, with modern practices of digitization and databasing of morphological traits providing exciting improvements in our ability to evaluate the degree of morphological convergence. Combining more detailed phenotypic data with the well-established field of genomics will enable scientists to make progress on an important goal in biology: to understand the degree to which genetic or molecular convergence is associated with phenotypic convergence. Although the fields of comparative biology or comparative genomics alone can separately reveal important insights into convergent evolution, here we suggest that the synergistic and complementary roles of natural history collection-derived phenomic data and comparative genomics methods can be particularly powerful in together elucidating the genomic basis of convergent evolution among higher taxa.

In Supplementary Table 1 of this Article, 23 samples (B10K-DU-029-32, B10K-DU-029-33, B10K-DU-029-36 to B10K-DU-029-44, B10K-DU- 029-46, B10K-DU-029-47, B10K-DU-029-49 to B10K-DU-029-53, B10K-DU- 029-75 to B10K-DU-029-77, B10K-DU-029-80, and B10K-DU-030-03; styled in boldface in the revised table) were assigned to the incorrect institution. Supplementary Table 1 has been amended to reflect the correct source institution for these samples, and associated data (tissue, museum ID/source specimen ID, site, state/province, latitude, longitude, date collected and sex) have been updated accordingly. The original table is provided as Supplementary Information to this Amendment, and the original Article has been corrected online.

Genes encode components of coevolved and interconnected networks. The effect of genotype on phenotype therefore depends on genotypic context through gene interactions known as epistasis. Epistasis is important in predicting phenotype from genotype for an individual. It is also examined in population studies to identify genetic risk factors in complex traits and to predict evolution under selection. Paradoxically, the effects of genotypic context in individuals and populations are distinct and sometimes contradictory. We argue that predicting genotype from phenotype for individuals based on population studies is difficult and, especially in human genetics, likely to result in underestimating the effects of genotypic context.

H3K9me3 (histone H3 modified with tri-methylation at lysine 9) is a hallmark of transcriptional silencing and heterochromatin. However, its global effects on the genome, including euchromatin, are less well understood. Here we develop Formaldehyde-Assisted Identification of Regulatory Elements (FAIRE) forC. elegansto examine the chromatin configuration of mutants that lack virtually all H3K9me3, while leaving H3K9me1 and H3K9me2 intact. We find that nucleosomes are mildly disrupted, and levels of H3K9me2 and H3K27me3 rise in mutant embryos. In addition to these expected changes, the most dramatic change occurs in euchromatin: many regions encompassing transcription start sites (TSSs) gain an average of two nucleosomes in mutants. The affected regions normally lack H3K9me3, revealing a locus non-autonomous role for H3K9me3. Affected TSSs are associated with genes that are active in epithelia and muscles, and implicated in development, locomotion, morphogenesis and transcription. Mutant embryos develop normally under ideal laboratory conditions but die when challenged, with defects in morphogenesis and development. Our findings reveal that H3K9me3 protects transcription start sites within euchromatin from nucleosome deposition. These results may be relevant to mammals, where diseases that disrupt the nuclear lamina and heterochromatin can alter epithelial and muscle gene expression.

Standard models of sex chromosome evolution propose that recombination suppression leads to the degeneration of the heterogametic chromosome, as is seen for the Y chromosome in mammals and the W chromosome in most birds. Unlike other birds, palaeognaths (ratites and tinamous) possess large non-degenerate regions on their sex chromosomes (PARs or pseudoautosomal regions), despite sharing the same sex determination region as neognaths (all other birds). It remains unclear why the large PARs of palaeognaths are retained over more than 100 MY of evolution, and the impact of these large PARs on sex chromosome evolution. To address this puzzle, we analysed Z chromosome evolution and gene expression across 12 palaeognaths, several of whose genomes have recently been sequenced. We confirm at the genomic levels that most palaeognaths (excepting some species of tinamous) retain large PARs. As in neognaths, we find that all palaeognaths have incomplete dosage compensation on the regions of the Z chromosome homologous to degenerated portions of the W (differentiated regions or DRs), but we find no evidence for enrichments of male-biased genes in PARs. We find limited evidence for increased evolutionary rates (faster-Z) either across the chromosome or in DRs for most palaeognaths with large PARs, but do recover signals of faster-Z evolution similar to neognaths in tinamou species with mostly degenerated W chromosomes (small PARs). Unexpectedly, in some species PAR-linked genes evolve faster on average than genes on autosomes. Increased TE density and longer introns in PARs of most palaeognaths compared to autosomes suggest that the efficacy of selection may be reduced in palaeognath PARs, contributing to the faster-Z evolution we observe. Our analysis shows that palaeognath Z chromosomes are atypical at the genomic level, but the evolutionary forces maintaining largely homomorphic sex chromosomes in these species remain elusive.

A core question in evolutionary biology is whether convergent phenotypic evolution is driven by convergent molecular changes in proteins or regulatory regions. We combined phylogenomic, developmental, and epigenomic analysis of 11 new genomes of paleognathous birds, including an extinct moa, to show that convergent evolution of regulatory regions, more so than protein-coding genes, is prevalent among developmental pathways associated with independent losses of flight. A Bayesian analysis of 284,001 conserved noncoding elements, 60,665 of which are corroborated as enhancers by open chromatin states during development, identified 2355 independent accelerations along lineages of flightless paleognaths, with functional consequences for driving gene expression in the developing forelimb. Our results suggest that the genomic landscape associated with morphological convergence in ratites has a substantial shared regulatory component.

Conservation of DNA sequence over evolutionary time is a strong indicator of function, and gain or loss of sequence conservation can be used to infer changes in function across a phylogeny. Changes in evolutionary rates on particular lineages in a phylogeny can indicate shared functional shifts, and thus can be used to detect genomic correlates of phenotypic convergence. However, existing methods do not allow easy detection of patterns of rate variation, which causes challenges for detecting convergent rate shifts or other complex evolutionary scenarios. Here we introduce PhyloAcc, a new Bayesian method to model substitution rate changes in conserved elements across a phylogeny. The method assumes several categories of substitution rate for each branch on the phylogenetic tree, estimates substitution rates per category, and detects changes of substitution rate as the posterior probability of a category switch. Simulations show that PhyloAcc can detect genomic regions with rate shifts in multiple target species better than previous methods and has a higher accuracy of reconstructing complex patterns of substitution rate changes than prevalent Bayesian relaxed clock models. We demonstrate the utility of PhyloAcc in two classic examples of convergent phenotypes: loss of flight in birds and the transition to marine life in mammals. In each case, our approach reveals numerous examples of conserved non-exonic elements with accelerations specific to the phenotypically convergent lineages. Our method is widely applicable to any set of conserved elements where multiple rate changes are expected on a phylogeny.

Palaeognathae represent one of the two basal lineages in modern birds, and comprise the volant (flighted) tinamous and the flightless ratites. Resolving palaeognath phylogenetic relationships has historically proved difficult, and short internal branches separating major palaeognath lineages in previous molecular phylogenies suggest that extensive incomplete lineage sorting (ILS) might have accompanied a rapid ancient divergence. Here, we investigate palaeognath relationships using genome-wide data sets of three types of noncoding nuclear markers, together totalling 20,850 loci and over 41 million base pairs of aligned sequence data. We recover a fully resolved topology placing rheas as the sister to kiwi and emu + cassowary that is congruent across marker types for two species tree methods (MP-EST and ASTRAL-II). This topology is corroborated by patterns of insertions for 4,274 CR1 retroelements identified from multi-species whole genome screening, and is robustly supported by phylogenomic subsampling analyses, with MP-EST demonstrating particularly consistent performance across subsampling replicates as compared to ASTRAL. In contrast, analyses of concatenated data supermatrices recover rheas as the sister to all other non-ostrich palaeognaths, an alternative that lacks retroelement support and shows inconsistent behavior under subsampling approaches. While statistically supporting the species tree topology, conflicting patterns of retroelement insertions also occur and imply high amounts of ILS across short successive internal branches, consistent with observed patterns of gene tree heterogeneity. Coalescent simulations indicate that the majority of observed topological incongruence among gene trees is consistent with coalescent variation rather than arising from gene tree estimation error alone, and estimated branch lengths for short successive internodes in the inferred species tree fall within the theoretical range encompassing the anomaly zone. Distributions of empirical gene trees confirm that the most common gene tree topology for each marker type differs from the species tree, signifying the existence of an empirical anomaly zone in palaeognaths.

High throughput sequencing (HTS) has revolutionized the field of ancient DNA (aDNA) by facilitating recovery of nuclear DNA for greater inference of evolutionary processes in extinct species than is possible from mitochondrial DNA alone. We used HTS to obtain ancient DNA from the little bush moa (Anomalopteryx didiformis), one of the iconic species of large, flightless birds that became extinct following human settlement of New Zealand in the 13 th century. In addition to a complete mitochondrial genome at 249.9X depth of coverage, we recover almost 900 Mb of the moa nuclear genome by mapping reads to a high quality reference genome for the emu (Dromaius novaehollandiae). This first nuclear genome assembly for moa covers approximately 75% of the 1.2 Gb emu reference with sequence contiguity sufficient to identify more than 85% of bird universal single-copy orthologs. From this assembly, we isolate 40 polymorphic microsatellites to serve as a community resource for future population-level studies in moa. We also compile data for a suite of candidate genes associated with vertebrate limb development and show that the wingless moa phenotype is likely not attributable to gene loss or pseudogenization among this candidate set. We also identify potential function-altering coding sequence variants in moa for future experimental assays.

The relative roles of regulatory and protein evolution in the origin and loss of convergent phenotypic traits is a core question in evolutionary biology. Here we combine phylogenomic, epigenomic and developmental data to show that convergent evolution of regulatory regions, but not protein-coding genes, is associated with flightlessness in palaeognathous birds, a classic example of a convergent phenotype. Eleven new genomes, including a draft genome from an extinct moa, resolve palaeognath phylogeny and show that the incidence of independent, convergent accelerations among 284,000 conserved non-exonic elements is significantly more frequent in ratites than other bird lineages. Ratite-specific acceleration of conserved regions and measures of open chromatin across eight tissues in the developing chick identify candidate regulatory regions that may have modified or lost function in ratites. Enhancer activity assays conducted in the early developing chicken forelimb confirm that volant versions of a conserved element in the first intron of the TEAD1 gene display conserved enhancer activity, whereas an accelerated flightless version fails to drive reporter gene expression. Our results show that convergent molecular changes associated with loss of flight are largely regulatory in nature.

Developmental genomics is a rapidly growing field, and high-quality genomes are a useful foundation for comparative developmental studies. A high-quality genome forms an essential reference onto which the data from numerous assays and experiments, including ChIP-seq, ATAC-seq, and RNA-seq, can be mapped. A genome also streamlines and simplifies the development of primers used to amplify putative regulatory regions for enhancer screens, cDNA probes for in situ hybridization, microRNAs (miRNAs) or short hairpin RNAs (shRNA) for RNA interference (RNAi) knockdowns, mRNAs for misexpression studies, and even guide RNAs (gRNAs) for CRISPR knockouts. Finally, much can be gleaned from comparative genomics alone, including the identification of highly conserved putative regulatory regions. This chapter provides an overview of laboratory and bioinformatics protocols for DNA extraction, library preparation, library quantification, and genome assembly, from fresh or frozen tissue to a draft avian genome. Generating a high-quality draft genome can provide a developmental research group with excellent resources for their study organism, opening the doors to many additional assays and experiments.

Native to Asia, the soft-skinned fruit pest Drosophila suzukii has recently invaded the United States and Europe. The eastern United States represents the most recent expansion of their range, and presents an opportunity to test alternative models of colonization history. Here, we investigate the genetic population structure of this invasive fruit fly, with a focus on the eastern United States. We sequenced six X-linked gene fragments from 246 individuals collected from a total of 12 populations. We examine patterns of genetic diversity within and between populations and explore alternative colonization scenarios using approximate Bayesian computation. Our results indicate high levels of nucleotide diversity in this species and suggest that the recent invasions of Europe and the continental United States are independent demographic events. More broadly speaking, our results highlight the importance of integrating population structure into demographic models, particularly when attempting to reconstruct invasion histories. Finally, our simulation results illustrate the general challenge in reconstructing invasion histories using genetic data and suggest that genome-level data are often required to distinguish among alternative demographic scenarios.

Thanks to a dramatic reduction in sequencing costs followed by a rapid development of bioinformatics tools, genome assembly and annotation have become accessible to many researchers in recent years. Among tetrapods, birds have genomes that display many features that facilitate their assembly and annotation, such as small genome size, low number of repeats and highly conserved genomic structure. However, we found that high genomic heterozygosity could have a great impact on the quality of the genome assembly of the thick-billed murre (Uria lomvia), an arctic colonial seabird. In this study, we tested the performance of three genome assemblers, ray/sscape, soapdenovo2 and platanus, in assembling the highly heterozygous genome of the thick-billed murre. Our results show that platanus, an assembler specifically designed for heterozygous genomes, outperforms the other two approaches and produces a highly contiguous (N50 = 15.8 Mb) and complete genome assembly (93% presence of genes from the Benchmarking Universal Single Copy Ortholog [BUSCO] gene set). Additionally, we annotated the thick-billed murre genome using a homology-based approach that takes advantage of the genomic resources available for birds and other taxa. Our study will be useful for those researchers who are approaching assembly and annotation of highly heterozygous genomes, or genomes of species of conservation concern, and/or who have limited financial resources.

With the sequencing of 12 complete euchromatic Drosophila genomes, the genus Drosophila is a leading model for comparative genomics. In this review, we discuss the novel insights into evolutionary processes afforded by the newly available genomic sequences when placed in the context of the phylogeny. We focus on three levels: insights into whole-genome content, such as changes in genome size and content across the phylogeny; insights into large-scale patterns of divergence and conservation, such as selective constraints on genes and chromosome-level evolution of sex chromosomes; and insights into finer-scale processes in individual lineages and genes, such as lineage-specific evolution in response to ecological context. As the field of comparative genomics is still young, we also discuss current challenges, such as the development of more sophisticated evolutionary models to capture nonequilibrium processes and the improvement of assembly and alignment algorithms to better capture uncertainty in the data.

Short-read sequencing techniques provide the opportunity to capture genome-wide sequence data in a single experiment. A current challenge is to identify questions that shallow-depth genomic data can address successfully and to develop corresponding analytical methods that are statistically sound. Here, we apply the Roche/454 platform to survey natural variation in strains of Drosophila melanogaster from an African (n = 3) and a North American (n = 6) population. Reads were aligned to the reference D. melanogaster genomic assembly, single nucleotide polymorphisms were identified, and nucleotide variation was quantified genome wide. Simulations and empirical results suggest that nucleotide diversity can be accurately estimated from sparse data with as little as 0.2× coverage per line. The unbiased genomic sampling provided by random short-read sequencing also allows insight into distributions of transposable elements and copy number polymorphisms found within populations and demonstrates that short-read sequencing methods provide an efficient means to quantify variation in genome organization and content. Continued development of methods for statistical inference of shallow-depth genome-wide sequencing data will allow such sparse, partial data sets to become the norm in the emerging field of population genomics.

The availability of complete genome sequence from 12 Drosophila species presents the opportunity to examine how natural selection has affected patterns of gene family evolution and sequence divergence among different components of the innate immune system. We have identified orthologs and paralogs of 245 Drosophila melanogaster immune-related genes in these recently sequenced genomes. Genes encoding effector proteins, and to a lesser extent genes encoding recognition proteins, are much more likely to vary in copy number across species than genes encoding signaling proteins. Furthermore, we can trace the apparent recent origination of several evolutionarily novel immune-related genes and gene families. Using codon-based likelihood methods, we show that immune-system genes, and especially those encoding recognition proteins, evolve under positive darwinian selection. Positively selected sites within recognition proteins cluster in domains involved in recognition of microorganisms, suggesting that molecular interactions between hosts and pathogens may drive adaptive evolution in the Drosophila immune system.