Across diverse organisms, the strength and ecological drivers of sexual selection vary enormously. In birds, some of the families with the most elaborate plumage and display-such as birds of paradise, manakins, and cotingas-are also specialist frugivores, yet links between shifts in diet, changes in breeding system, and the evolution of elaborate traits are poorly understood. We focus on manakins, a radiation of frugivorous Neotropical birds well known for spectacular courtship rituals and colorful plumage, and present an integrative analysis of the transition in both diet and mating systems in this clade to examine the causes and consequences of strong sexual selection. In manakins, we find reduced genetic diversity on the Z sex chromosome relative to autosomes, a predicted signature of sexual selection. We also identify targets of positive selection across the manakin radiation, including genes related to muscle function, visual perception, and the transition to frugivory. Among these, we observe selection on sugar-sensing taste receptors, as well as on lactase-phlorizin hydrolase, implicated in the consumption of chemically defended fruits. For both, we confirm that selection signatures correspond to functional changes and infer the relative time of these changes, as well as of shifts in diet, breeding systems, and plumage coloration: elaborated traits evolved subsequent to changes in mating systems and after key physiological changes facilitating fruit-eating. Altogether, these results suggest that intensified frugivory set the stage for the radiation of one of the planet's most colorful and acrobatic avian lineages.

In a hybrid zone between two tropical lekking birds, yellow male plumage of one species has introgressed asymmetrically replacing white plumage of another via sexual selection. Here, we present a detailed analysis of the plumage trait to uncover its physical and genetic bases and trace its evolutionary history. We determine that the carotenoid lutein underlies the yellow phenotype and describe microstructural feather features likely to enhance color appearance. These same features reduce predicted water shedding capacity of feathers, a potential liability in the tropics. Through genome-scale DNA sequencing of hybrids and each species in the genus, we identify BCO2 as the major gene responsible for the color polymorphism. The BCO2 gene tree and genome-wide allele frequency patterns suggest that carotenoid-pigmented collars initially arose in a third species and reached the hybrid zone through historical gene flow. Complex interplay between sexual selection and hybridization has thus shaped phenotypes of these species, where conspicuous sexual traits are key to male reproductive success.

Hybrid zones can be studied by modeling clines of trait variation (e.g., morphology, genetics) over a linear transect. Yet, hybrid zones can also be spatially complex, can shift over time, and can even lead to the formation of hybrid lineages with the right combination of dispersal and vicariance. We reassessed Sibley's (1950) gradient between Collared Towhee (Pipilo ocai) and Spotted Towhee (Pipilo maculatus) in Central Mexico to test whether it conformed to a typical tension-zone cline model. By comparing historical and modern data, we found that cline centers for genetic and phenotypic traits have not shifted over the course of 70 years. This equilibrium suggests that secondary contact between these species, which originally diverged over 2 million years ago, likely dates to the Pleistocene. Given the amount of mtDNA divergence, parental ends of the cline have very low autosomal nuclear differentiation (FST = 0.12). Dramatic and coincident cline shifts in mtDNA and throat color suggest the possibility of sexual selection as a factor in differential introgression, while a contrasting cline shift in green back color hints at a role for natural selection. Supporting the idea of a continuum between clinal variation and hybrid lineage formation, the towhee gradient can be analyzed as one population under isolation-by-distance, as a two-population cline, and as three lineages experiencing divergence with gene flow. In the middle of the gradient, a hybrid lineage has become partly isolated, likely due to forested habitat shrinking and fragmenting as it moved upslope after the last glacial maximum and a stark environmental transition. This towhee system offers a window into the potential outcomes of hybridization across a dynamic landscape including the creation of novel genomic and phenotypic combinations and incipient hybrid lineages.

Abstract Damage to the adult central nervous system often leads to long-term disruptions in function due to the limited capacity for neurological recovery. The central nervous system of the Mediterranean field cricket, Gryllus bimaculatus, shows an unusual capacity for compensatory plasticity, most obviously in the auditory system and the cercal escape system. In both systems, unilateral sensory disruption leads the central circuitry to compensate by forming and/or strengthening connections with the contralateral sensory organ. While this compensatory plasticity relies on robust dendritic sprouting and novel synapse formation in the auditory system, the compensatory plasticity in the cercal escape circuitry shows little obvious dendritic sprouting and instead may rely on shifts in excitatory and inhibitory synaptic strength. In order to better understand what types of molecular pathways might underlie this compensatory shift in the cercal system, we used a multiple k-mer approach to assemble a terminal ganglion transcriptome that included ganglia collected one, three, and seven days after unilateral cercal ablation in adult, male animals. We performed differential expression analysis using EdgeR and DESeq2 and examined Gene Ontologies to identify candidates potentially involved in this plasticity. Enriched GO terms included those related to the ubiquitin-proteosome protein degradation system, chromatin-mediated transcriptional pathways, and the GTPase-related signaling system. Further exploration of these GO terms will provide a clearer picture of the processes involved in compensatory recovery of the cercal escape system in the cricket and can be compared and contrasted with the distinct pathways that have been identified upon deafferentation of the auditory system in this same animal.

Cryptic and hybridizing species may lack diagnostic taxonomic characters leaving researchers with semi-diagnostic ones. Identification based on such characters is probabilistic, the probability of correct identification depending on the species composition in a mixed population. Here we test the possibilities of applying a semi-diagnostic conchological character for distinguishing two cryptic species of blue mussels, Mytilus edulis and M. trossulus. These ecologically, stratigraphically and economically important molluscs co-occur and hybridize in many areas of the North Atlantic and the neighboring Arctic. Any cues for distinguishing them in sympatry without genotyping would save much research effort. Recently these species have been shown to statistically differ in the White Sea, where a simple character of the shell was used to distinguish two mussel morphotypes. In this paper, we analyzed the associations between morphotypes and species-specific genotypes based on an abundant material from the waters of the Kola Peninsula (White Sea, Barents Sea) and a more limited material from Norway, the Baltic Sea, Scotland and the Gulf of Maine. The performance of the “morphotype test” for species identification was formally evaluated using approaches from evidence-based medicine. Interspecific differences in the morphotype frequencies were ubiquitous and unidirectional, but their scale varied geographically (from 75% in the White Sea to 15% in the Baltic Sea). In addition, salinity-related variation of this character within M. edulis was revealed in the Arctic Barents Sea. For every studied region, we established relationships between the proportions of the morphotypes in the populations as well as between the proportions of the morphotypes in samples and the probabilities of mussels of different morphotypes being M. trossulus and M. edulis. We provide recommendations for the application of the morphotype test to mussels from unstudied contact zones and note that they may apply equally well to other taxa identified by semi-diagnostic traits.

Most adult organisms are limited in their capacity to recover from neurological damage. The auditory system of the Mediterranean field cricket, Gryllus bimaculatus, presents a compelling model for investigating neuroplasticity due to its unusual capabilities for structural reorganization into adulthood. Specifically, the dendrites of the central auditory neurons of the prothoracic ganglion sprout in response to the loss of auditory afferents. Deafferented auditory dendrites grow across the midline, a boundary they normally respect, and form functional synapses with the contralateral auditory afferents, restoring tuning-curve specificity. The molecular pathways underlying these changes are entirely unknown. Here, we used a multiple k-mer approach to re-assemble a previously reported prothoracic ganglion transcriptome that included ganglia collected one, three, and seven days after unilateral deafferentation in adult, male animals. We used EdgeR and DESeq2 to perform differential expression analysis and we examined Gene Ontologies to further understand the potential molecular basis of this compensatory anatomical plasticity. Enriched GO terms included those related to protein translation and degradation, enzymatic activity, and Toll signaling. Extracellular space GO terms were also enriched and included the upregulation of several protein yellow family members one day after deafferentation. Investigation of these regulated GO terms help to provide a broader understanding of the types of pathways that might be involved in this compensatory growth and can be used to design hypotheses around identified molecular mechanisms that may be involved in this unique example of adult structural plasticity.

AbstractMost adult organisms are limited in their capacity to recover from neurological damage. The auditory system of the Mediterranean field cricket, Gryllus bimaculatus, presents a compelling model for investigating neuroplasticity due to its unusual capabilities for structural reorganization into adulthood. Specifically, the dendrites of the central auditory neurons of the prothoracic ganglion sprout in response to the loss of auditory afferents. Deafferented auditory dendrites grow across the midline, a boundary they normally respect, and form functional synapses with the contralateral auditory afferents, restoring tuning-curve specificity. The molecular pathways underlying these changes are entirely unknown. Here, we used a multiple k-mer approach to re-assemble a previously reported prothoracic ganglion transcriptome that included ganglia collected one, three, and seven days after unilateral deafferentation in adult, male animals. We used EdgeR and DESeq2 to perform differential expression analysis and we examined Gene Ontologies to further understand the potential molecular basis of this compensatory anatomical plasticity. Enriched GO terms included those related to protein translation and degradation, enzymatic activity, and Toll signaling. Extracellular space GO terms were also enriched and included the upregulation of several protein yellow family members one day after deafferentation. Investigation of these regulated GO terms help to provide a broader understanding of the types of pathways that might be involved in this compensatory growth and can be used to design hypotheses around identified molecular mechanisms that may be involved in this unique example of adult structural plasticity.

The biogeochemistry of the Gulf of Maine (GOM) is rapidly changing in response to the changing climate, including rising temperatures, acidification, and declining primary productivity. These impacts are projected to worsen over the next 100 y and will apply selective pressure on populations of marine calcifiers. This study investigates the transcriptome expression response to these changes in ecologically and economically important marine calcifiers, blue mussels. Wild mussels (Mytilus edulis and Mytilus trossulus) were sampled from sites spanning the GOM and exposed to two different biogeochemical water conditions: (1) present-day conditions in the GOM and (2) simulated future conditions, which included elevated temperature, increased acidity, and decreased food supply. Patterns of gene expression were measured using RNA sequencing from 24 mussel samples and contrasted between ambient and future conditions. The net calcification rate, a trait predicted to be under climate-induced stress. was measured for each individual over a 2-wk exposure period and used as a covariate along with gene expression patterns. Generalized linear models, with and without the calcification rate, were used to identify differentially expressed transcripts between ambient and future conditions. The comparison revealed transcripts that likely comprise a core stress response characterized by the induction of molecular chaperones, genes involved in aerobic metabolism, and indicators of cellular stress. Furthermore, the model contrasts revealed transcripts that may be associated with individual variation in calcification rate and suggest possible biological processes that may have downstream effects on calcification phenotypes, such as zinc-ion binding and protein degradation. Overall, these findings contribute to the understanding of blue mussel adaptive responses to imminent climate change and suggest metabolic pathways are resilient in variable environments.

Captive breeding and reintroduction can be effective conservation tools, but the origin of founders for such programs is key for success. The red siskin Spinus cucullatus, a bird from northern South America, is Critically Endangered in Venezuela due to decades of trapping for the illegal wildlife trade. As a result, many red siskins are held in captivity worldwide, but several potential problems with captive birds make considering founders from the wild more desirable. A recently discovered population of red siskins in Guyana presents such an opportunity, but, due to its disjunct distribution from the main range, the possibility of genetic differentiation is a concern. We used a variety of standard and novel analyses of amplified fragment length polymorphism (AFLP) and mitochondrial (mtDNA) markers to evaluate genetic divergence of the Guyana (GU) population, using 5 individuals from GU and 13 captive birds of presumed Venezuelan (‘VE’) origin. All analyses of nuclear loci revealed 2 clusters separating GU from ‘VE’ individuals, with FST values varying from 0.15 to 0.24, depending on assumptions about individual relatedness. Furthermore, all 5 GU birds shared an mtDNA haplotype that differed by 2 or more substitutions (0.11%) from the 3 ‘VE’ haplotypes. The GU population thus appears to be differentiated from the ‘VE’ population in both nuclear and mtDNA. While further genetic evidence is needed, these data suggest that the GU population is not an optimal source of founders for recovery efforts in Venezuela, and should be treated as a separate elemental conservation unit until additional data are available.

Shellfish aquaculture (in particular oyster cultivation) has the potential to play a significant role in refreshing the coastal economy in the state of Maine (USA). Although ocean warming and acidification are often listed as the primary manifestations of climate change in the marine environment, the issue of expanding geographic ranges of parasites and pathogens is an equally serious threat to shellfish populations. Protozoan parasites of the genera Bonamia, Haplosporidium, and Perkinsus, the bacterial pathogen responsible for Roseovarius oyster disease, and the disease condition is known as disseminated neoplasia are currently recognized as some of the significant threats to natural and farmed bivalve populations. We have analyzed the peer-reviewed literature for reports of these pathogens/conditions in Maine. Most reports focus on oysters from the Damariscotta River Estuary and are the result of directed studies into the biology of a particular pathogen rather than the result of intensive monitoring programs. The sampling effort could be interpreted in several ways including; the pathogens/conditions not impacting the annual harvest to any great extent due to limited distribution of the causative agent, the oysters developing a tolerance to the diseases, or just a lack of resources directed at studying this topic. With the shellfish aquaculture industry expected to grow in the next several decades, we recommend that a rigorous and sustained survey of parasitic diseases and believe that such an effort is fundamental to the success, resilience, and well-being of Maine's shellfish aquaculture industry.

A key component to understanding the evolutionary response to a changing climate is linking underlying genetic variation to phenotypic variation in stress response. Here, we use a genome-wide association approach (GWAS) to understand the genetic architecture of calcification rates under simulated climate stress. We take advantage of the genomic gradient across the blue mussel hybrid zone (Mytilus edulis and Mytilus trossulus) in the Gulf of Maine (GOM) to link genetic variation with variance in calcification rates in response to simulated climate change. Falling calcium carbonate saturation states are predicted to negatively impact many marine organisms that build calcium carbonate shells - like blue mussels. We sampled wild mussels and measured net calcification phenotypes after exposing mussels to a climate change' common garden, where we raised temperature by 3 degrees C, decreased pH by 0.2 units and limited food supply by filtering out planktonic particles >5m, compared to ambient GOM conditions in the summer. This climate change exposure greatly increased phenotypic variation in net calcification rates compared to ambient conditions. We then used regression models to link the phenotypic variation with over 170000 single nucleotide polymorphism loci (SNPs) generated by genotype by sequencing to identify genomic locations associated with calcification phenotype, and estimate heritability and architecture of the trait. We identified at least one of potentially 2-10 genomic regions responsible for 30% of the phenotypic variation in calcification rates that are potential targets of natural selection by climate change. Our simulations suggest a power of 13.7% with our study's average effective sample size of 118 individuals and rare alleles, but a power of >90% when effective sample size is 900.

Phylogenomics, the use of large-scale data matrices in phylogenetic analyses, has been viewed as the ultimate solution to the problem of resolving difficult nodes in the tree of life. However, it has become clear that analyses of these large genomic data sets can also result in conflicting estimates of phylogeny. Here, we use the early divergences in Neoaves, the largest clade of extant birds, as a "model system" to understand the basis for incongruence among phylogenomic trees. We were motivated by the observation that trees from two recent avian phylogenomic studies exhibit conflicts. Those studies used different strategies: 1) collecting many characters [$\sim$ 42 mega base pairs (Mbp) of sequence data] from 48 birds, sometimes including only one taxon for each major clade; and 2) collecting fewer characters ($\sim$ 0.4 Mbp) from 198 birds, selected to subdivide long branches. However, the studies also used different data types: the taxon-poor data matrix comprised 68% non-coding sequences whereas coding exons dominated the taxon-rich data matrix. This difference raises the question of whether the primary reason for incongruence is the number of sites, the number of taxa, or the data type. To test among these alternative hypotheses we assembled a novel, large-scale data matrix comprising 90% non-coding sequences from 235 bird species. Although increased taxon sampling appeared to have a positive impact on phylogenetic analyses the most important variable was data type. Indeed, by analyzing different subsets of the taxa in our data matrix we found that increased taxon sampling actually resulted in increased congruence with the tree from the previous taxon-poor study (which had a majority of non-coding data) instead of the taxon-rich study (which largely used coding data). We suggest that the observed differences in the estimates of topology for these studies reflect data-type effects due to violations of the models used in phylogenetic analyses, some of which may be difficult to detect. If incongruence among trees estimated using phylogenomic methods largely reflects problems with model fit developing more "biologically-realistic" models is likely to be critical for efforts to reconstruct the tree of life. [Birds; coding exons; GTR model; model fit; Neoaves; non-coding DNA; phylogenomics; taxon sampling.].

The maintenance or breakdown of reproductive isolation is an observable outcome of secondary contact between species. In cases where hybrids beyond the F1 are formed, the representation of each species' ancestry can vary dramatically among genomic regions. This genomic heterogeneity in ancestry and introgression can offer insight into evolutionary processes, particularly if introgression is compared in multiple hybrid zones. Similarly, considerable heterogeneity exists across the genome in the extent to which populations and species have diverged, reflecting the combined effects of different evolutionary processes on genetic variation. We studied hybridization across two hybrid zones of two phenotypically well-differentiated bird species in Mexico (Pipilo maculatus and P.ocai), to investigate genomic heterogeneity in differentiation and introgression. Using genotyping-by-sequencing (GBS) and hierarchical Bayesian models, we genotyped 460 birds at over 41000 single nucleotide polymorphism (SNP) loci. We identified loci exhibiting extreme introgression relative to the genome-wide expectation using a Bayesian genomic cline model. We also estimated locus-specific F-ST and identified loci with exceptionally high genetic divergence between the parental species. We found some concordance of locus-specific introgression in the two independent hybrid zones (6-20% of extreme loci shared across zones), reflecting areas of the genome that experience similar gene flow when the species interact. Additionally, heterogeneity in introgression and divergence across the genome revealed another subset of loci under the influence of locally specific factors. These results are consistent with a history in which reproductive isolation has been influenced by a common set of loci in both hybrid zones, but where local environmental and stochastic factors also lead to genomic differentiation.

Hybridization can either reinforce or erode species boundaries; therefore, hybrid zones offer a natural experimental setting in which to assess the dynamics of reproductive isolation. Secondary contact zones, in particular, present a partial separation of the original divergence mechanisms and the subsequent genomic architecture of reproductive isolation (or lack thereof). The spatial context of secondary contact and its consequent effect on dispersal play vital roles on the contact’s outcome. In a hybrid complex between two towhee species in Mexico, Pipilo maculatus and Pipilo ocai, two major hybrid gradients provide natural replicates for comparison. However, genetic analyses demonstrate significant divergence between geographically separate parental populations of each species and divergence of populations within each hybrid zone. The two hybrid transects (Teziutlán and Transvolcanic) are distinct and evidence suggests allelic introgression both across the species boundary and between the two transects. Habitat corridors for dispersal represent functional connectivity hotspots where the two transects meet. Both habitat connectivity and genetic differentiation between geographically disparate parental populations appear to influence the dynamics of gene flow across the hybrid gradients. In southern sympatric populations (Mt. Orizaba and Oaxaca) where morphological evidence for hybridization is scarce, opposing species’ alleles appear to traverse through the hybrid zones rather than arising from cryptic local hybridization. These results illustrate the importance of environmentally mediated gene flow in the context of secondary contact as an important force influencing evolutionary trajectory.

Hybridization can either reinforce or erode species boundaries; therefore, hybrid zones offer a natural experimental setting in which to assess the dynamics of reproductive isolation. Secondary contact zones, in particular, present a partial separation of the original divergence mechanisms and the subsequent genomic architecture of reproductive isolation (or lack thereof). The spatial context of secondary contact and its consequent effect on dispersal play vital roles on the contact's outcome. In a hybrid complex between two towhee species in Mexico, Pipilo maculatus and Pipilo ocai, two major hybrid gradients provide natural replicates for comparison. However, genetic analyses demonstrate significant divergence between geographically separate parental populations of each species and divergence of populations within each hybrid zone. The two hybrid transects (Teziutlan and Transvolcanic) are distinct and evidence suggests