Response to Langmore and Spottiswoode: “Visual Trickery in Avian Brood Parasites”

Organismic and Evolutionary Biolog

Next-Generation QTL Mapping: Crowdsourcing SNPs, Without Pedigrees

For many molecular ecologists, the mantra and mission of the field of ecological genomics could be encapsulated by the phrase ‘to find the genes that matter’ (Mitchell-Olds 2001; Rockman 2012). This phrase of course refers to the early hope and current increasing success in the search for genes whose variation underlies phenotypic variation and fitness in natural populations. In the years since the modern incarnation of the field of ecological genomics, many would agree that the low-hanging fruit has, at least in principle, been plucked: we now have several elegant examples of genes whose variation influences key adaptive traits in natural populations, and these examples have revealed important insights into the architecture of adaptive variation (Hoekstra et al. 2006; Shapiro et al. 2009; Chan et al. 2010). But how well will these early examples, often involving single genes of large effect on discrete or near-discrete phenotypes, represent the dynamics of adaptive change for the totality of phenotypes in nature? Will traits exhibiting continuous rather than discrete variation in natural populations have as simple a genetic basis as these early examples suggest (Prasad et al. 2012; Rockman 2012)? Two papers in this issue (Robinson et al. 2013; Santure et al. 2013) not only suggest answers to these questions but also provide useful extensions of statistical approaches for ecological geneticists to study the genetics of continuous variation in nature. Together these papers, by the same research groups studying evolution in a natural population of Great Tits (Parus major), provide a glimpse of what we should expect as the field begins to dissect the genetic basis of what is arguably the most common type of variation in nature, and how genome-wide surveys of variation can be applied to natural populations without pedigrees.Organismic and Evolutionary Biolog

The Evolution of Intron Size in Amniotes: A Role for Powered Flight?

Intronic DNA is a major component of eukaryotic genes and genomes and can be subject to selective constraint and have functions in gene regulation. Intron size is of particular interest given that it is thought to be the target of a variety of evolutionary forces and has been suggested to be linked ultimately to various phenotypic traits, such as powered flight. Using whole-genome analyses and comparative approaches that account for phylogenetic nonindependence, we examined interspecific variation in intron size variation in three data sets encompassing from 12 to 30 amniotes genomes and allowing for different levels of genome coverage. In addition to confirming that intron size is negatively associated with intron position and correlates with genome size, we found that on average mammals have longer introns than birds and nonavian reptiles, a trend that is correlated with the proliferation of repetitive elements in mammals. Two independent comparisons between flying and nonflying sister groups both showed a reduction of intron size in volant species, supporting an association between powered flight, or possibly the high metabolic rates associated with flight, and reduced intron/genome size. Small intron size in volant lineages is less easily explained as a neutral consequence of large effective population size. In conclusion, we found that the evolution of intron size in amniotes appears to be non-neutral, is correlated with genome size, and is likely influenced by powered flight and associated high metabolic rates.Human Evolutionary BiologyOrganismic and Evolutionary Biolog

A Species Tree for the Australo-Papuan Fairy-Wrens and Allies (Aves: Maluridae)

We explored the efficacy of species tree methods at the family level in birds, using the Australo-Papuan Fairy-wrens (Passeriformes: Maluridae) as a model system. Fairy-wrens of the genus Malurus are known for high intensities of sexual selection, resulting in some cases in rapid speciation. This history suggests that incomplete lineage sorting (ILS) of neutrally evolving loci could be substantial, a situation that could compromise traditional methods of combining loci in phylogenetic analysis. Using 18 molecular markers (5 anonymous loci, 7 exons, 5 introns, and 1 mitochondrial DNA locus), we show that gene tree monophyly across species could be rejected for 16 of 18 loci, suggesting substantial ILS at the family level in these birds. Using the software Concaterpillar, we also detect three statistically distinct clusters of gene trees among the 18 loci. Despite substantial variation in gene trees, species trees constructed using four different species tree estimation methods (BEST, BUCKy, and STAR) were generally well supported and similar to each other and to the concatenation tree, with a few mild discordances at nodes that could be explained by rapid and recent speciation events. By contrast, minimizing deep coalescences produced a species tree that was topologically more divergent from those of the other methods as measured by multidimensional scaling of trees. Additionally, gene and species trees were topologically more similar in the BEST analysis, presumably because of the species tree prior employed in BEST which appropriately assumes that gene trees are correlated with each other and with the species tree. Among the 18 loci, we also discovered 102 independent indel markers, which also proved phylogenetically informative, primarily among genera, and displayed a $\sim4$-fold bias towards deletions. As suggested in earlier work, the grasswrens (Amytornis) are sister to the rest of the family and the emu-wrens (Stipiturus) are sister to fairy-wrens (Malurus, Clytomyias). Our study shows that ILS is common at the family level in birds yet, despite this, species tree methods converge on broadly similar results for this family.Organismic and Evolutionary Biolog

A Maximum Pseudo-Likelihood Approach for Estimating Species Trees under the Coalescent Model

Background: Several phylogenetic approaches have been developed to estimate species trees from collections of gene trees. However, maximum likelihood approaches for estimating species trees under the coalescent model are limited. Although the likelihood of a species tree under the multispecies coalescent model has already been derived by Rannala and Yang, it can be shown that the maximum likelihood estimate (MLE) of the species tree (topology, branch lengths, and population sizes) from gene trees under this formula does not exist. In this paper, we develop a pseudo-likelihood function of the species tree to obtain maximum pseudo-likelihood estimates (MPE) of species trees, with branch lengths of the species tree in coalescent units. Results: We show that the MPE of the species tree is statistically consistent as the number $M$ of genes goes to infinity. In addition, the probability that the MPE of the species tree matches the true species tree converges to 1 at rate $O(M^{-1})$. The simulation results confirm that the maximum pseudo-likelihood approach is statistically consistent even when the species tree is in the anomaly zone. We applied our method, Maximum Pseudo-likelihood for Estimating Species Trees (MP-EST) to a mammal dataset. The four major clades found in the MP-EST tree are consistent with those in the Bayesian concatenation tree. The bootstrap supports for the species tree estimated by the MP-EST method are more reasonable than the posterior probability supports given by the Bayesian concatenation method in reflecting the level of uncertainty in gene trees and controversies over the relationship of four major groups of placental mammals. Conclusions: MP-EST can consistently estimate the topology and branch lengths (in coalescent units) of the species tree. Although the pseudo-likelihood is derived from coalescent theory, and assumes no gene flow or horizontal gene transfer (HGT), the MP-EST method is robust to a small amount of HGT in the dataset. In addition, increasing the number of genes does not increase the computational time substantially. The MP-EST method is fast for analyzing datasets that involve a large number of genes but a moderate number of species.Organismic and Evolutionary Biolog

A house finch (Haemorhous mexicanus) spleen transcriptome reveals intra- and interspecific patterns of gene expression, alternative splicing and genetic diversity in passerines

Background: With its plumage color dimorphism and unique history in North America, including a recent population expansion and an epizootic of Mycoplasma gallisepticum (MG), the house finch (Haemorhous mexicanus) is a model species for studying sexual selection, plumage coloration and host-parasite interactions. As part of our ongoing efforts to make available genomic resources for this species, here we report a transcriptome assembly derived from genes expressed in spleen. Results: We characterize transcriptomes from two populations with different histories of demography and disease exposure: a recently founded population in the eastern US that has been exposed to MG for over a decade and a native population from the western range that has never been exposed to MG. We utilize this resource to quantify conservation in gene expression in passerine birds over approximately 50 MY by comparing splenic expression profiles for 9,646 house finch transcripts and those from zebra finch and find that less than half of all genes expressed in spleen in either species are expressed in both species. Comparative gene annotations from several vertebrate species suggest that the house finch transcriptomes contain ~15 genes not yet found in previously sequenced vertebrate genomes. The house finch transcriptomes harbour ~85,000 SNPs, ~20,000 of which are non-synonymous. Although not yet validated by biological or technical replication, we identify a set of genes exhibiting differences between populations in gene expression (n = 182; 2% of all transcripts), allele frequencies (76 FST ouliers) and alternative splicing as well as genes with several fixed non-synonymous substitutions; this set includes genes with functions related to double-strand break repair and immune response. Conclusions: The two house finch spleen transcriptome profiles will add to the increasing data on genome and transcriptome sequence information from natural populations. Differences in splenic expression between house finch and zebra finch imply either significant evolutionary turnover of splenic expression patterns or different physiological states of the individuals examined. The transcriptome resource will enhance the potential to annotate an eventual house finch genome, and the set of gene-based high-quality SNPs will help clarify the genetic underpinnings of host-pathogen interactions and sexual selection

Major histocompatibility complex class I evolution in songbirds: universal primers, rapid evolution and base compositional shifts in exon 3

Genes of the Major Histocompatibility Complex (MHC) have become an important marker for the investigation of adaptive genetic variation in vertebrates because of their critical role in pathogen resistance. However, despite significant advances in the last few years the characterization of MHC variation in non-model species still remains a challenging task due to the redundancy and high variation of this gene complex. Here we report the utility of a single pair of primers for the cross-amplification of the third exon of MHC class I genes, which encodes the more polymorphic half of the peptide-binding region (PBR), in oscine passerines (songbirds; Aves: Passeriformes), a group especially challenging for MHC characterization due to the presence of large and complex MHC multigene families. In our survey, although the primers failed to amplify exon 3 from two suboscine passerine birds, they amplified exon 3 of multiple MHC class I genes in all 16 species of oscine songbirds tested, yielding a total of 120 sequences. The 16 songbird species belong to 14 different families, primarily within the Passerida, but also in the Corvida. Using a conservative approach based on the analysis of cloned amplicons (n = 16) from each species, we found between 3 and 10 MHC sequences per individual. Each allele repertoire was highly divergent, with the overall number of polymorphic sites per species ranging from 33 to 108 (out of 264 sites) and the average number of nucleotide differences between alleles ranging from 14.67 to 43.67. Our survey in songbirds allowed us to compare macroevolutionary dynamics of exon 3 between songbirds and non-passerine birds. We found compelling evidence of positive selection acting specifically upon peptide-binding codons across birds, and we estimate the strength of diversifying selection in songbirds to be about twice that in non-passerines. Analysis using comparative methods suggest weaker evidence for a higher GC content in the 3rd codon position of exon 3 in non-passerine birds, a pattern that contrasts with among-clade GC patterns found in other avian studies and may suggests different mutational mechanisms. Our primers represent a useful tool for the characterization of functional and evolutionarily relevant MHC variation across the hyperdiverse songbirds

The Maluridae: Inferring Avian Biology and Evolutionary History from DNA Sequences

The Australo-Papuan fairy-wrens, emu-wrens and grasswrens comprise the passerine family Maluridae. They have long been known for their spectacular plumages, remarkable behavioural ecology and intriguing biogeography. The family has provided an ideal model with which to explore how phylogenetic and phylogeographic analyses of DNA-sequence data can inform understanding of evolutionary history and present-day biology. We review what has been learned of the phylogeny of the group and the phylogeographic history of individual species. We conclude that there is now a strong framework within which to pursue the remaining species-level taxonomic issues, and to extend ecological and behavioural studies into a new era of more detailed genetic questions such as the role of gene–environment interactions in adaptation. We highlight some remaining examples of such questions and discuss how they might be addressed.Organismic and Evolutionary Biolog

Feather development genes and associated regulatory innovation predate the origin of Dinosauria

The evolution of avian feathers have recently been illuminated by fossils and the identification of genes involved in feather patterning and morphogenesis. However, molecular studies have focused mainly on protein-coding genes. Using comparative genomics and more than 600,000 conserved regulatory elements, we show that patterns of genome evolution in the vicinity of feather genes are consistent with a major role for regulatory innovation in the evolution of feathers. Rates of innovation at feather regulatory elements exhibit an extended period of innovation with peaks in the ancestors of amniotes and archosaurs. We estimate that 86% of such regulatory elements were present prior to the origin of Dinosauria. On the branch leading to modern birds, we detect a strong signal of regulatory innovation near IGFBP2 and IGFBP5, which have roles in body size reduction, and may represent a genomic signature for the miniaturization of dinosaurian body size preceding the origin of flight.Organismic and Evolutionary Biolog

Dense sphere packings from optimized correlation functions

Elementary smooth functions (beyond contact) are employed to construct pair correlation functions that mimic jammed disordered sphere packings. Using the g2-invariant optimization method of Torquato and Stillinger [J. Phys. Chem. B 106, 8354, 2002], parameters in these functions are optimized under necessary realizability conditions to maximize the packing fraction phi and average number of contacts per sphere Z. A pair correlation function that incorporates the salient features of a disordered packing and that is smooth beyond contact is shown to permit a phi of 0.6850: this value represents a 45% reduction in the difference between the maximum for congruent hard spheres in three dimensions, pi/sqrt{18} ~ 0.7405, and 0.64, the approximate fraction associated with maximally random jammed (MRJ) packings in three dimensions. We show that, surprisingly, the continued addition of elementary functions consisting of smooth sinusoids decaying as r^{-4} permits packing fractions approaching pi/sqrt{18}. A translational order metric is used to discriminate between degrees of order in the packings presented. We find that to achieve higher packing fractions, the degree of order must increase, which is consistent with the results of a previous study [Torquato et al., Phys. Rev. Lett. 84, 2064, 2000].Comment: 26 pages, 9 figures, 1 table; added references, fixed typos, simplified argument and discussion in Section IV