Search

Symbiosis often occurs between partners with distinct life history characteristics and dispersal mechanisms. Many bacterial symbionts have genomes comprising multiple replicons with distinct rates of evolution and horizontal transmission. Such differences might drive differences in population structure between hosts and symbionts and among the elements of the divided genomes of bacterial symbionts. These differences might, in turn, shape the evolution of symbiotic interactions and bacterial evolution. Here we use whole genome resequencing of a hierarchically structured sample of 191 strains of Sinorhizobium meliloti collected from 21 locations in southern Europe to characterize population structures of this bacterial symbiont, which forms a root nodule symbiosis with the host plant Medicago truncatula. S. meliloti genomes showed high local (within-site) variation and little isolation by distance. This was particularly true for the two symbiosis elements, pSymA and pSymB, which have population structures that are similar to each other, but distinct from both the bacterial chromosome and the host plant. Given limited recombination on the chromosome, compared to the symbiosis elements, distinct population structures may result from differences in effective gene flow. Alternatively, positive or purifying selection, with little recombination, may explain distinct geographical patterns at the chromosome. Discordant population structure between hosts and symbionts indicates that geographically and genetically distinct host populations in different parts of the range might interact with genetically similar symbionts, potentially minimizing local specialization.

The mutualism between legumes and rhizobia is clearly the product of past coevolution. However, the nature of ongoing evolution between these partners is less clear. To characterize the nature of recent coevolution between legumes and rhizobia, we used population genomic analysis to characterize selection on functionally annotated symbiosis genes as well as on symbiosis gene candidates identified through a two-species association analysis. For the association analysis, we inoculated each of 202 accessions of the legume host Medicago truncatula with a community of 88 Sinorhizobia (Ensifer) meliloti strains. Multistrain inoculation, which better reflects the ecological reality of rhizobial selection in nature than single-strain inoculation, allows strains to compete for nodulation opportunities and host resources and for hosts to preferentially form nodules and provide resources to some strains. We found extensive host by symbiont, that is, genotype-by-genotype, effects on rhizobial fitness and some annotated rhizobial genes bear signatures of recent positive selection. However, neither genes responsible for this variation nor annotated host symbiosis genes are enriched for signatures of either positive or balancing selection. This result suggests that stabilizing selection dominates selection acting on symbiotic traits and that variation in these traits is under mutation-selection balance. Consistent with the lack of positive selection acting on host genes, we found that among-host variation in growth was similar whether plants were grown with rhizobia or N-fertilizer, suggesting that the symbiosis may not be a major driver of variation in plant growth in multistrain contexts.

Rhizobium strains that nodulate Egyptian winter clover belong to two symbiotic species, Rhizobium aegyptiacum sv. trifolii and Rhizobium bangladeshense sv. trifolii. Rhizobia from Egyptian clover are often tolerant of extreme environmental conditions and are generally poorly characterized, despite the importance of this crop species. Here, we report on the comparative genomic analysis of four stress-tolerant, highly effective nitrogen-fixing, Egyptian-clover-nodulating rhizobia, R. aegyptiacum sv. trifolii strain Rhiz950 and Rhizobium bangladeshense sv. trifolii strains Rhiz1002, Rhiz1017 and Rhiz1024, as well as other diverse Rhizobium strains. While Rhizobium bangladeshense sv. trifolii strains Rhiz1002, Rhiz1017, and Rhiz1024 had similar genome sizes of around 6.6 Mb, the salt tolerant strain Rhiz950 had a slightly larger genome size of 7.4 Mb. All four strains lacked genes responsible for methylation and sulfation of nod factors (nodS and nodHPQ). They also lacked nodZ, and instead had the functionally-equivalent nod factor fucosylation gene nolK. No strains contained homocitrate synthase (encoded by nifV) that is essential for free-living nitrogen fixing by diazotrophs. A number of genes unique to the Egyptian clover strains and to other rhizobia are potential candidates for variation in stress tolerance and host-range.

Recent policies in emerging economies (EMEs) have sought to increase the shares of households and firms that participate in the domestic banking system amid the backdrop that, compared to advanced economies (AEs), EMEs have drastically lower shares of both firm and household financial participation. We build a model with extensive margins of firm and household financial participation and equilibrium unemployment to study the impact on EME labor-market and business cycle dynamics of increasing the shares of firm and household financial participation in EMEs to AE levels. We find that a joint increase in household and firm domestic financial participation considerably narrow the differences in some, but not all, business cycle moments between EMEs and AEs. Critically, the impact of increasing only firm or only household participation can differ substantially from the impact of a joint increase. We stress the relevance of the extensive margin of domestic financial participation for labor-market and business cycle dynamics and the fact that promoting a joint increase in the shares of household and firm financial participation in EMEs can lead to smoother cyclical fluctuations, particularly in labor markets.

Symbiosis often occurs between partners with distinct life history characteristics and dispersal mechanisms. Bacterial symbionts often have genomes comprised of multiple replicons with distinct rates of evolution and horizontal transmission. Such differences might drive differences in population genetic structure between hosts and symbionts and among the elements of the divided genomes of bacterial symbionts. These differences might, in turn, shape the evolution of symbiotic interactions and bacterial evolution. Here we use whole genome resequencing of a hierarchically-structured sample of 191 strains of Ensifer meliloti collected from 21 locations in the native range to characterize population structure of this bacterial symbiont and its host plant Medicago truncatula. E. meliloti genomes showed high local (within-site) variation and little isolation by distance. This was particularly true for the two symbiosis elements pSymA and pSymB, which have population structures that are similar to each other, but distinct from both the bacterial chromosome and the host plant. The differences in population structure may result from differences in mobility or selection driving bacterial adaptation to life in the soil versus in association with plants. Discordant population structure between hosts and symbionts indicates that geographically and genetically distinct host populations in different parts of the range might interact with genetically similar symbionts, potentially minimizing the potential for local specialization.

Changing climate and land‐use practices have the potential to bring previously isolated populations of pest insects into new sympatry. This heightens the need to better understand how differing patterns of host–plant association, and unique endosymbionts, serve to promote genetic isolation or integration. We addressed these factors in populations of potato psyllid, Bactericera cockerelli (Šulc), a generalist herbivore that vectors a bacterial pathogen (Candidatus Liberibacter solanacearum, causal pathogen of zebra chip disease) of potato (Solanum tuberosum L.). Genome‐wide SNP data revealed two major genetic clusters—psyllids collected from potato crops were genetically similar to psyllids found on a common weed, Lycium spp., but dissimilar from those found on another common non‐crop host, Solanum dulcamara L. Most psyllids found on Lycium spp. and potato represented a single mitochondrial cytochrome oxidase I (COI) haplotype that has been suggested to not be native to the region, and whose arrival may have been concurrent with zebra chip disease first emerging. The putatively introduced COI haplotype usually co‐occurred with endosymbiotic Wolbachia, while the putatively resident COI haplotype generally did not. Genetic intermediates between the two genetic populations of insects were rare, consistent with recent sympatry or reproductive isolation, although admixture patterns of apparent hybrids were consistent with introgression of genes from introduced into resident populations. Our results suggest that both host–plant associations and endosymbionts are shaping the population genetic structure of sympatric psyllid populations associated with different non‐crop hosts. It is of future interest to explicitly examine vectorial capacity of the two populations and their potential hybrids, as population structure and hybridization might alter regional vector capacity and disease outbreaks.

Coevolution is predicted to depend on how the genetic diversity of interacting species is geographically structured. Plant-microbe symbioses such as the legume-rhizobium mutualism are ecologically and economically important, but distinct life history and dispersal mechanisms for these host and microbial partners, plus dynamic genome composition in bacteria, present challenges for understanding spatial genetic processes in these systems. Here we study the model rhizobium Ensifer meliloti using a hierarchically-structured sample of 191 strains from 21 sites in the native range and compare its population structure to that of its host plant Medicago truncatula. We find high local genomic variation and minimal isolation by distance across the rhizobium genome, particularly at the two symbiosis elements pSymA and pSymB, which have evolutionary histories and population structures that are similar to each other but distinct from both the chromosome and the host. While the chromosome displays weak isolation by distance, it is uncorrelated with hosts. Patterns of discordant population structure among elements with the bacterial genome has implications for bacterial adaptation to life in the soil versus symbiosis, while discordant population genetic structure of hosts and microbes might restrict local adaptation of species to each other and give rise to phenotypic mismatches in coevolutionary traits.

Coevolution is predicted to depend on how the genetic diversity of interacting species is geographically structured. Plant-microbe symbioses such as the legume-rhizobium mutualism are ecologically and economically important, but distinct life history and dispersal mechanisms for these host and microbial partners, plus dynamic genome composition in bacteria, present challenges for understanding spatial genetic processes in these systems. Here we study the model rhizobium Ensifer meliloti using a hierarchically-structured sample of 191 strains from 21 sites in the native range and compare its population structure to that of its host plant Medicago truncatula. We find high local genomic variation and minimal isolation by distance across the rhizobium genome, particularly at the two symbiosis elements pSymA and pSymB, which have evolutionary histories and population structures that are similar to each other but distinct from both the chromosome and the host. While the chromosome displays weak isolation by distance, it is uncorrelated with hosts. Patterns of discordant population structure among elements with the bacterial genome has implications for bacterial adaptation to life in the soil versus symbiosis, while discordant population genetic structure of hosts and microbes might restrict local adaptation of species to each other and give rise to phenotypic mismatches in coevolutionary traits.

We develop a theory that focuses on the general equilibrium and long-run macroeconomic consequences of trends in job utility. Given secular increases in job utility, work hours per capita can remain approximately constant over time even if the income effect of higher wages on labor supply exceeds the substitution effect. In addition, secular improvements in job utility can be substantial relative to welfare gains from ordinary technological progress. These two implications are connected by an equation ‡flowing from optimal hours choices: improvements in job utility that have a signi cant effect on labor supply tend to have large welfare effects.

Genetic studies of legume symbiosis with nitrogen-fixing rhizobial bacteria have traditionally focused on nodule and nitrogen-fixation phenotypes when hosts are inoculated with a single rhizobial strain. These approaches overlook the potential effect of host genes on rhizobial fitness (i.e. how many rhizobia are released from host nodules) and strain-specific effects of host genes (i.e. genome × genome interactions). Using Medicago truncatula mutants in the recently described nodule-specific PLAT domain (NPD) gene family, we show how inoculating plants with a mixed inoculum of 68 rhizobial strains (Ensifer meliloti) via a select-and-resequence approach can be used to efficiently assay host mutants for strain-specific effects of late-acting host genes on interacting bacteria. The deletion of a single NPD gene (npd2) or all five members of the NPD gene family (npd1-5) differentially altered the frequency of rhizobial strains in nodules even though npd2 mutants had no visible nodule morphology or N-fixation phenotype. Also, npd1-5 nodules were less diverse and had larger populations of colony-forming rhizobia despite their smaller size. Lastly, NPD mutations disrupt a positive correlation between strain fitness and wild-type host biomass. These changes indicate that the effects of NPD proteins are strain dependent and that NPD family members are not redundant with regard to their effects on rhizobial strains. Association analyses of the rhizobial strains in the mixed inoculation indicate that rhizobial genes involved in chromosome segregation, cell division, GABA metabolism, efflux systems, and stress tolerance play an important role in the strain-specific effects of NPD genes.

We use a panel structural vector autoregression methodology to study the impact of global financial risk shocks in emerging market economies (EMEs) versus small open advanced economies (SOAEs). The short-term responses of GDP, investment, and un- employment are similar in EMEs compared to SOAEs, but their medium-term responses differ considerably. In EMEs the shock propagates less across GDP and investment, and the response of unemployment is prominently more subdued. Further analysis suggests that while in EMEs these dynamics can be traced back to the response of country interest rates to global financial risk shocks and to the average degree of firm informality, in SOAEs these factors are not at play. Instead, for SOAEs the dynamics of bank credit can matter to some degree, but more broadly a global financial risk shock passes through more directly into the economy.

We document a negative and significant relationship between domestic financial development and unemployment volatility in developing and emerging economies (DEMEs) and the absence of such relationship in advanced economies (AEs). A business cycle labor search model with firm heterogeneity, collateral constraints, and interfirm input credit capital can quantitatively rationalize these facts. Greater financial development is associated with lower usage of input credit capital, greater bank credit, and greater capital accumulation, all of which make firms more resilient in the presence of financial shocks. Firms’ increased shock resiliency stabilizes employment decisions, ultimately leading to smoother unemployment fluctuations. Then, by establishing explicit linkages between firms, interfirm input credit acts as an important mechanism that fosters lower volatility across firms under greater financial development. As such, increases in financial development have a smaller impact in stabilizing unemployment in economies with higher average bank credit-GDP ratios and less input credit usage, which is the case of AEs, compared to economies with lower average bank-credit GDP ratios and more input credit usage, which is the case of DEMEs.

Understanding the genomic basis of adaptation is critical for understanding evolutionary processes and predicting how species will respond to environmental change. Spinner dolphins in the eastern tropical Pacific (ETP) present a unique system for studying adaptation. Within this large geographical region are four spinner dolphin ecotypes with weak neutral genetic divergence and no obvious barriers to gene flow, but strong spatial variation in morphology, behaviour and habitat. These ecotypes have large population sizes, which could reduce the effects of drift and facilitate selection. To identify genomic regions putatively under divergent selective pressures between ecotypes, we used genome scans with 8994 RADseq single nucleotide polymorphisms (SNPs) to identify population differentiation outliers and genotype–environment association outliers. Gene ontology enrichment analyses indicated that outlier SNPs from both types of analyses were associated with multiple genes involved in social behaviour and hippocampus development, including 15 genes associated with the human social disorder autism. Evidence for divergent selection on social behaviour is supported by previous evidence that these spinner dolphin ecotypes differ in mating systems and associated social behaviours. In particular, three of the ETP ecotypes probably have a polygynous mating system characterized by strong premating competition among males, whereas the fourth ecotype probably has a polygynandrous mating system characterized by strong postmating competition such as sperm competition. Our results provide evidence that selection for social behaviour may be an evolutionary force driving diversification of spinner dolphins in the ETP, potentially as a result of divergent sexual selection associated with different mating systems. Future studies should further investigate the potential adaptive role of the candidate genes identified here, and could probably find further signatures of selection using whole genome sequence data.

Horizontal transfer (HT) alters the repertoire of symbiosis genes in rhizobial genomes and may play an important role in the on-going evolution of the rhizobia–legume symbiosis. To gain insight into the extent of HT of symbiosis genes with different functional roles (nodulation, N-fixation, host benefit and rhizobial fitness), we conducted comparative genomic and selection analyses of the full-genome sequences from 27 rhizobial genomes. We find that symbiosis genes experience high rates of HT among rhizobial lineages but also bear signatures of purifying selection (low Ka : Ks). HT and purifying selection appear to be particularly strong in genes involved in initiating the symbiosis (e.g. nodulation) and in genome-wide association candidates for mediating benefits provided to the host. These patterns are consistent with rhizobia adapting to the host environment through the loss and gain of symbiosis genes, but not with host-imposed positive selection driving divergence of symbiosis genes through recurring bouts of positive selection.

Reforms aimed at reducing the regulatory barriers to firm entry and bolstering greater firm creation have gained significant prominence in economies beyond the U.S. Given these facts and amid growing work on the link between firm creation and cyclical housing-market dynamics, we document a significant, positive, and robust crosscountry relationship between the level of new firm creation and the cyclical volatility of house prices. Using a business cycle model with endogenous firm entry, housing, and housing finance constraints, we show how greater average firm entry can be a powerful amplification mechanism of housing-market shocks amid housing finance constraints, with these shocks and frictions playing a key role in quantitatively rationalizing the link between firm creation and house-price volatility across countries.

The environmental context of the nitrogen-fixing mutualism between leguminous plants and rhizobial bacteria varies over space and time. Variation in resource availability, population density, and composition likely affect the ecology and evolution of rhizobia and their symbiotic interactions with hosts. We examined how host genotype, nitrogen addition, rhizobial density, and community complexity affected selection on 68 rhizobia strains in the Ensifer meliloti - Medicago truncatula mutualism. As expected, the host genotype had the most substantial effect on the size, number, and strain composition of root nodules (the symbiotic organ). The understudied environmental variable of rhizobial density had a more significant effect on strain frequency in nodules than the addition of low nitrogen levels. Higher inoculum density resulted in a nodule community that was less diverse and more beneficial but only in the context of the more selective host genotype. Higher density resulted in more diverse and less beneficial nodule communities with the less selective host. Density effects on strain composition deserve additional scrutiny as they can create eco-evolutionary feedback. Lastly, we found that relative strain rankings were stable across increasing community complexity (community complexity (2, 3, 8, or 68 strains). This unexpected result suggests that higher-order interactions between strains are rare in the context of host nodule formation and development. Taken together, our empirical work highlights the importance of developing new theoretical predictions that incorporate density dependence. Further, it has translational relevance for overcoming establishment barriers in bio-inoculants and motivating host breeding programs that maintain beneficial plant-microbe interactions across diverse agro-ecological contexts.IMPORTANCELegume cash, forage, and cover crops establish beneficial associations with rhizobial bacteria who perform biological nitrogen fixation (BNF)—providing Nitrogen (N) fertilizer to plants without the economic and greenhouse gas emission costs of chemical N inputs. Here, for the first time, we examine the relative influence of three environmental factors that vary in agricultural fields on strain relative fitness in nodules when scores rhizobial strains compete. In addition to manipulating Nitrogen, we also use two biotic variables that have rarely been examined: the rhizobial community’s density and complexity. Taken together, our results suggest 1) breeding legume varieties that select beneficial strains despite environmental variation are possible, 2) changes in rhizobial population densities that occur routinely in agricultural fields could drive evolutionary changes in rhizobia populations, and 3) the lack of higher-order interactions between strains will allow the high-throughput assessments of rhizobia winners and losers during plant interactions.

Tasmanian devils (Sarcophilus harrisii) are evolving in response to a unique transmissible cancer, devil facial tumour disease (DFTD), first described in 1996. Persistence of wild populations and the recent emergence of a second independently evolved transmissible cancer suggest that transmissible cancers may be a recurrent feature in devils. Here we compared signatures of selection across temporal scales to determine whether genes or gene pathways under contemporary selection (6-8 generations) have also been subject to historical selection (65-85 million years), and test for recurrent selection in devils. First, we used a targeted sequencing approach, RAD-capture, to identify genomic regions subject to rapid evolution in approximately 2,500 devils in six populations as DFTD spread across the species range. We documented genome-wide contemporary evolution, including 186 candidate genes related to cell cycling and immune response. Then we used a molecular evolution approach to identify historical positive selection in devils compared to other marsupials and found evidence of selection in 1,773 genes. However, we found limited overlap across time scales, with historical selection detected in only 16 contemporary candidate genes, and no overlap in enriched functional gene sets. Our results are consistent with a novel, multi-locus evolutionary response of devils to DFTD. Our results can inform management actions to conserve adaptive potential of devils by identifying high priority targets for genetic monitoring and maintenance of functional diversity in managed populations.

Taking inefficiencies from taxation as given, a well-known public finance literature shows that the elasticity of taxable income (ETI) is a sufficient statistic for assessing the deadweight loss (DWL) from taxing labor income in a static neoclassical framework. Using a theoretical approach, we revisit this result from the vantage point of a general equilibrium macroeconomic model with labor search frictions. We show that, in this context, and against the backdrop of inefficient taxation, DWL can be up to 38 percent higher than the ETI under a range of reasonable parametric assumptions. Externalities arising from market participants not taking into account the impact of changes in their search- and vacancy-posting activities on other market participants can amplify this divergence substantially. However, with theoretical precision, we show how the wedge between the ETI and DWL can be controlled for, using readily observable variables.

Changing climate and land‐use practices have the potential to bring previously isolated populations of pest insects into new sympatry. This heightens the need to better understand how differing patterns of host–plant association, and unique endosymbionts, serve to promote genetic isolation or integration. We addressed these factors in populations of potato psyllid, Bactericera cockerelli (Šulc), a generalist herbivore that vectors a bacterial pathogen (Candidatus Liberibacter solanacearum, causal pathogen of zebra chip disease) of potato (Solanum tuberosum L.). Genome‐wide SNP data revealed two major genetic clusters—psyllids collected from potato crops were genetically similar to psyllids found on a common weed, Lycium spp., but dissimilar from those found on another common non‐crop host, Solanum dulcamara L. Most psyllids found on Lycium spp. and potato represented a single mitochondrial cytochrome oxidase I (COI) haplotype that has been suggested to not be native to the region, and whose arrival may have been concurrent with zebra chip disease first emerging. The putatively introduced COI haplotype usually co‐occurred with endosymbiotic Wolbachia, while the putatively resident COI haplotype generally did not. Genetic intermediates between the two genetic populations of insects were rare, consistent with recent sympatry or reproductive isolation, although admixture patterns of apparent hybrids were consistent with introgression of genes from introduced into resident populations. Our results suggest that both host–plant associations and endosymbionts are shaping the population genetic structure of sympatric psyllid populations associated with different non‐crop hosts. It is of future interest to explicitly examine vectorial capacity of the two populations and their potential hybrids, as population structure and hybridization might alter regional vector capacity and disease outbreaks.

Natural enemies often move among habitats to track prey and resources. Indeed, biocontrol often depends on natural enemies dispersing into crops after disturbances such as tillage and pesticide applications. However, the small size of many natural enemies makes it difficult to observe such movements. Here we used genetic relatedness among entomopathogenic nematodes, tiny, soil-dwelling, and thus cryptic natural enemies, to infer their movement across an agricultural landscape. We collected strains of two nematode species, Heterorhabditis bacteriophora and Steinernema feltiae, by placing sentinel hosts into eight irrigated Solanum tuberosum fields across arid central Washington State. We then used restriction site associated DNA sequencing to generate single nucleotide polymorphisms and infer relatedness among strains across our study sites. We identified 3,367 and 138,286 polymorphic loci for H. bacteriophora and S. feltiae, respectively. Genetic differentiation for both species increased with greater distance between sites, although there was considerable variation. While strains collected from the same field were generally more closely related than those from different sites, for both species, strains from different fields were sometimes quite closely related. Altogether, our results suggest a surprising amount of genetic similarity among nematodes from distant sites, despite the presumably limited distances that can be traversed by individual worms. This is consistent with the nematodes moving, at least occasionally, between far-apart locations. More generally, we suggest that recent advances in population genomics are providing powerful new tools for mapping natural enemy movement across broad landscapes.

Emerging infectious diseases are increasingly recognized in species' declines and extinctions. Landscape genetics can be used as a tool to predict disease emergence and spread. The Tasmanian devil is threatened with extinction by a nearly 100% fatal transmissible cancer, which has spread across 95% of the species' geographic range in 20 years. Here, we present a landscape genetic analysis in the last remaining uninfected parts of the Tasmanian devil's geographic range to: describe population genetic structure, characterize genetic diversity, and test the influence of landscape variables on Tasmanian devil gene flow to assess the potential for disease spread. In contrast to previous genetic studies on Tasmanian devils showing evidence for two genetic populations island-wide, our genetic based assignment tests and spatial principal components analyses suggest at least two, and possibly three, populations in a study area that is approximately 15% of the size of the overall species' geographic range. Positive spatial autocorrelation declined at about 40 km, in contrast to 80 km in eastern populations, highlighting the need for range-wide genetic studies. Strong genetic structure was found between devils in the northern part of the study area and those found south of Macquarie Harbour, with weaker structure found between the northeastern and northwestern portion of our study area. Consistent with previous work, we found low overall genetic diversity, likely owing to a combination of founder effects and extreme weather events thousands of years ago that likely caused large-scale population declines. We also found possible signs of recent bottlenecks, perhaps resulting from forest clearing for dairy farming in the central part of the study area. This human disturbance also may have contributed to weak genetic structuring detected between the northeastern and northwestern part of the study area. Individual-based least cost path modeling showed limited influence of landscape variables on gene flow, with weak effects of variation in elevation in the northeast. In the northwest, however, landscape genetic models did not perform better than the null isolation-by-distance model. At the larger spatial scale of the northern part of the study area, elevation and temperatures were negatively correlated with gene flow, consistent with low dispersal suitability of higher elevation habitats that have lower temperatures and dense, wet vegetation. Overall, Tasmanian devils are a highly vagile species for which dispersal and gene flow appear to be influenced little by landscape features, and spread of devil facial tumor disease to the remaining portion of the devil's geographic range seems imminent. Nonetheless, strong genetic structure found between the northern and southern portions of our study area, combined with low densities and limited possible colonization of DFTD from the east suggest there is some time for implementation of management strategies.

Tasmanian devils face a combination of threats to persistence, including Devil Facial Tumor Disease (DFTD), an epidemic transmissible cancer. We used RAD sequencing to investigate genome-wide patterns of genetic diversity and geographic population structure. Consistent with previous results, we found very low genetic diversity in the species as a whole, and we detected two broad genetic clusters occupying the northwestern portion of the range, and the central and eastern portions. However, these two groups overlap across a broad geographic area, and differentiation between them is modest (FST = 0.1081). Our results refine the geographic extent of the zone of mixed ancestry and substructure within it, potentially informing management of genetic variation that existed in pre-diseased populations of the species. DFTD has spread across both genetic clusters, but recent evidence points to a genomic response to selection imposed by DFTD. Any allelic variation for resistance to DFTD may be able to spread across the devil population under selection by DFTD, and/or be present as standing variation in both genetic regions.

Landscape genomics studies focus on identifying candidate genes under selection via spatial variation in abiotic environmental variables, but rarely by biotic factors such as disease. The Tasmanian devil (Sarcophilus harrisii) is found only on the environmentally heterogeneous island of Tasmania and is threatened with extinction by a nearly 100% fatal, transmissible cancer, devil facial tumor disease (DFTD). Devils persist in regions of long-term infection despite epidemiological model predictions of species’ extinction, suggesting possible adaptation to DFTD. Here, we test the extent to which spatial variation and genetic diversity are associated with the abiotic environment and/or by DFTD. We employ genetic-environment association (GEAs) analyses using a RAD-capture panel consisting of 6,886 SNPs from 3,286 individuals sampled pre- and post-disease arrival. Pre-disease, we find significant correlations of allele frequencies with environmental variables, including 349 unique loci linked to 64 genes, suggesting local adaptation to abiotic environment. Following DFTD arrival, however, we detected few of the pre-disease candidate loci, but instead frequencies of candidate loci linked to 14 genes correlated with disease prevalence. Loss of apparent signal of abiotic local adaptation following disease arrival suggests swamping by the strong selection imposed by DFTD. Further support for this result comes from the fact that post-disease candidate loci are in linkage disequilibrium with genes putatively involved in immune response, tumor suppression and apoptosis. This suggests the strength GEA associations of loci with the abiotic environment are swamped resulting from the rapid onset of a biotic selective pressure.

Whether hybridization generates or erodes species diversity has long been debated, but to date most studies have been conducted at small taxonomic scales. Salamanders (order Caudata) represent a taxonomic order in which hybridization plays a prevalent ecological and evolutionary role. We employed a recently developed model of trait-dependent diversification to test the hypothesis that hybridization impacts the diversification dynamics of species that are currently hybridizing. We find strong evidence supporting this hypothesis, showing that hybridizing salamander lineages have significantly greater net-diversification rates than non-hybridizing lineages. This pattern is driven by concurrently increased speciation rates and decreased extinction rates in hybridizing lineages. Our results support the hypothesis that hybridization can act as a generative force in macroevolutionary diversification.

Genetic studies of legume symbiosis with nitrogen-fixing rhizobial bacteria have traditionally focused on nodule and nitrogen-fixation phenotypes when hosts are inoculated with a single rhizobial strain. These approaches overlook the potential effect of host genes on rhizobial fitness (i.e. how many rhizobia are released from host nodules) and strain-specific effects of host genes (i.e. genome × genome interactions). Using Medicago truncatula mutants in the recently described nodule-specific PLAT domain (NPD) gene family, we show how inoculating plants with a mixed inoculum of 68 rhizobial strains (Ensifer meliloti) via a select-and-resequence approach can be used to efficiently assay host mutants for strain-specific effects of late-acting host genes on interacting bacteria. The deletion of a single NPD gene (npd2) or all five members of the NPD gene family (npd1–5) differentially altered the frequency of rhizobial strains in nodules even though npd2 mutants had no visible nodule morphology or N-fixation phenotype. Also, npd1–5 nodules were less diverse and had larger populations of colony-forming rhizobia despite their smaller size. Lastly, NPD mutations disrupt a positive correlation between strain fitness and wild-type host biomass. These changes indicate that the effects of NPD proteins are strain dependent and that NPD family members are not redundant with regard to their effects on rhizobial strains. Association analyses of the rhizobial strains in the mixed inoculation indicate that rhizobial genes involved in chromosome segregation, cell division, GABA metabolism, efflux systems, and stress tolerance play an important role in the strain-specific effects of NPD genes.

Understanding factors that cause species' geographic range limits is a major focus in ecology and evolution. The central marginal hypothesis (CMH) predicts that species cannot adapt to conditions beyond current geographic range edges because genetic diversity decreases from core to edge due to smaller, more isolated edge populations. We employed a population genomics framework using 24 235-33 112 SNP loci to test major predictions of the CMH in the ongoing invasion of the cane toad (Rhinella marina) in Australia. Cane toad tissue samples were collected along broad-scale, core-to-edge transects across their invasive range. Geographic and ecological core areas were identified using GIS and habitat suitability indices from ecological niche modelling. Bayesian clustering analyses revealed three genetic clusters, in the northwest invasion-front region, northeast precipitation-limited region and southeast cold temperature-limited region. Core-to-edge patterns of genetic diversity and differentiation were consistent with the CMH in the southeast, but were not supported in the northeast and showed mixed support in the northwest. Results suggest cold temperatures are a likely contributor to southeastern range limits, consistent with CMH predictions. In the northeast and northwest, ecological processes consisting of a steep physiological barrier and ongoing invasion dynamics, respectively, are more likely explanations for population genomic patterns than the CMH.

Genome-wide association analyses are a powerful approach for identifying gene function. These analyses are becoming commonplace in studies of humans, domesticated animals, and crop plants but have rarely been conducted in bacteria. We applied association analyses to 20 traits measured in Ensifer meliloti, an agriculturally and ecologically important bacterium because it fixes nitrogen when in symbiosis with leguminous plants. We identified candidate alleles and gene presence-absence variants underlying variation in symbiosis traits, antibiotic resistance, and use of various carbon sources; some of these candidates are in genes previously known to affect these traits whereas others were in genes that have not been well characterized. Our results point to the potential power of association analyses in bacteria, but also to the need to carefully evaluate the potential for false associations., Genome-wide association studies (GWAS) can identify genetic variants responsible for naturally occurring and quantitative phenotypic variation. Association studies therefore provide a powerful complement to approaches that rely on de novo mutations for characterizing gene function. Although bacteria should be amenable to GWAS, few GWAS have been conducted on bacteria, and the extent to which nonindependence among genomic variants (e.g., linkage disequilibrium [LD]) and the genetic architecture of phenotypic traits will affect GWAS performance is unclear. We apply association analyses to identify candidate genes underlying variation in 20 biochemical, growth, and symbiotic phenotypes among 153 strains of Ensifer meliloti. For 11 traits, we find genotype-phenotype associations that are stronger than expected by chance, with the candidates in relatively small linkage groups, indicating that LD does not preclude resolving association candidates to relatively small genomic regions. The significant candidates show an enrichment for nucleotide polymorphisms (SNPs) over gene presence-absence variation (PAV), and for five traits, candidates are enriched in large linkage groups, a possible signature of epistasis. Many of the variants most strongly associated with symbiosis phenotypes were in genes previously identified as being involved in nitrogen fixation or nodulation. For other traits, apparently strong associations were not stronger than the range of associations detected in permuted data. In sum, our data show that GWAS in bacteria may be a powerful tool for characterizing genetic architecture and identifying genes responsible for phenotypic variation. However, careful evaluation of candidates is necessary to avoid false signals of association. IMPORTANCE Genome-wide association analyses are a powerful approach for identifying gene function. These analyses are becoming commonplace in studies of humans, domesticated animals, and crop plants but have rarely been conducted in bacteria. We applied association analyses to 20 traits measured in Ensifer meliloti, an agriculturally and ecologically important bacterium because it fixes nitrogen when in symbiosis with leguminous plants. We identified candidate alleles and gene presence-absence variants underlying variation in symbiosis traits, antibiotic resistance, and use of various carbon sources; some of these candidates are in genes previously known to affect these traits whereas others were in genes that have not been well characterized. Our results point to the potential power of association analyses in bacteria, but also to the need to carefully evaluate the potential for false associations.

This paper presents a model in which mismatch employment arises in a constrained efficient equilibrium. In the decentralized economy, however, mismatch gives rise to a congestion externality whereby heterogeneous job seekers fail to internalize how their individual actions affect the labor market outcomes of competitors in a common unemployment pool. We provide an analytic characterization of this distortion, assess the distributional nature of the associated welfare effects, and relate it to the relative productivity of low- and high-skilled workers competing for similar jobs.

Ensifer meliloti (formerly Rhizobium meliloti and Sinorhizobium meliloti) is a model bacterium for understanding legume–rhizobial symbioses. The tripartite genome of E. meliloti consists of a chromosome, pSymA and pSymB, and in some instances strain-specific accessory plasmids. The majority of previous sequencing studies have relied on the use of assemblies generated from short read sequencing, which leads to gaps and assembly errors. Here we used PacBio-based, long-read assemblies and were able to assemble, de novo, complete circular replicons. In this study, we sequenced, de novo-assembled and analysed 10 E. meliloti strains. Sequence comparisons were also done with data from six previously published genomes. We identified genome differences between the replicons, including mol% G+C and gene content, nucleotide repeats, and transposon-associated loci. Additionally, genomic rearrangements both within and between replicons were identified, providing insight into evolutionary processes at the structural level. There were few cases of inter-replicon gene transfer of core genes between the main replicons. Accessory plasmids were more similar to pSymA than to either pSymB or the chromosome, with respect to gene content, transposon content and G+C content. In our population, the accessory plasmids appeared to share an open genome with pSymA, which contains many nodulation- and nitrogen fixation-related genes. This may explain previous observations that horizontal gene transfer has a greater effect on the content of pSymA than pSymB, or the chromosome, and why some rhizobia show unstable nodulation phenotypes on legume hosts.

Identifying the genetic architecture of complex phenotypes is a central goal of modern biology, particularly for disease-related traits. Genome-wide association methods are a classical approach for identifying the genomic basis of variation in disease phenotypes, but such analyses are particularly challenging in natural populations due to sample size difficulties. Extensive mark-recapture data, strong linkage disequilibrium, and a lethal transmissible cancer make the Tasmanian devil (Sarcophilus harrisii) an ideal model for such an association study. We used a RAD-capture approach to genotype 624 devils at ~16,000 loci and then used association analyses to assess the heritability of three cancer-related phenotypes: infection case-control (where cases were infected devils and controls were devils that were never infected), age of first infection, and survival following infection. The SNP array explained much of the phenotypic variance for female survival (>80%) and female case-control (>61%). We found that a few large-effect SNPs explained much of the variance for female survival (~5 SNPs explained >61% of the total variance) whereas more SNPs (~56) of smaller effect explained less of the variance for female case-control (~23% of the total variance). By contrast, these same SNPs did not account for a significant proportion of phenotypic variance in males, suggesting that the genetic bases of these traits and/or selection differ across sexes. Loci involved with cell adhesion and cell-cycle regulation underlay trait variation, suggesting that the devil immune system is rapidly evolving to recognize and potentially suppress cancer growth through these pathways. Overall, our study provided necessary data for genomics-based conservation and management in Tasmanian devils.

Tasmanian devils face a combination of threats to persistence, including Devil Facial Tumor Disease (DFTD), an epidemic transmissible cancer. We used RAD sequencing to investigate genome-wide patterns of genetic diversity and geographic population structure. Consistent with previous results, we found very low genetic diversity in the species as a whole, and we detected two broad genetic clusters occupying the northwestern portion of the range, and the central and eastern portions. However, these two groups overlap across a broad geographic area, and differentiation between them is modest (FST = 0.1081). Our results refine the geographic extent of the zone of mixed ancestry and substructure within it, potentially informing management of genetic variation that existed in pre-diseased populations of the species. DFTD has spread across both genetic clusters, but recent evidence points to a genomic response to selection imposed by DFTD. Any allelic variation for resistance to DFTD may be able to spread across the devil population under selection by DFTD, and/or be present as standing variation in both genetic regions.

Amid total factor productivity (TFP) shocks job-to-job flows amplify the volatility of unemployment, but the aggregate implications of job-to-job flows amid financial shocks are less understood. To develop such understanding we model a general equilibrium labor-search framework that incorporates on-the-job (OTJ) search and distinctly accounts for the differential impact of TFP and financial shocks. Surprisingly, we find that the interaction of OTJ search with financial shocks is sufficiently different from its interaction with TFP shocks so that, under standard calibrations, our model generates aggregate dynamics exceedingly in line with the behavior of key U.S. macro data across several decades and in the wake of the Global Financial Crisis (GFC) as well. Importantly, as in the data, the model yields relatively high volatilities of consumption, labor income, and unemployment. As such, our work contributes to resolving two limitations of current general equilibrium labor-search theory: under standard calibrations and amid TFP shocks, only, models without OTJ search generate implausibly low unemployment volatility, while models with OTJ search generate unemployment volatility closer to the data but at the expense of implausibly low consumption and labor-income volatility.

Many countries have large employment shares in micro and small firms that have limited access to formal financing and therefore rely on input credit. Such countries are mainly emerging and developing economies, whose business cycle dynamics are increasingly important for the global economy in light of the dramatic rise in international linkages and spillovers that have occurred over the last several decades. Emerging and developing economies implemented a host of countercyclical labor market policies amid the global financial crisis, but data limitations on high-frequency labor and job flows prevent a detailed empirical assessment of the effectiveness of these policies. To address this problem, we develop a business cycle model with frictional labor markets that is novel in light of its consistency with the employment and firm structure of emerging and developing economies. We use the model to assess the aggregate impact of key countercyclical labor market policies. We find that hiring subsidies and job intermediation services for large firms are particularly effective in aiding recoveries. Policies targeting smaller firms yield limited aggregate benefits and may even be detrimental to the recovery process. The labor market structure shapes sectoral allocation and explains the economy's differential response to policy.

Herbivores often move among spatially interspersed host plants, tracking high-quality resources through space and time. This dispersal is of particular interest for vectors of plant pathogens. Existing molecular tools to track such movement have yielded important insights, but often provide insufficient genetic resolution to infer spread at finer spatiotemporal scales. Here, we explore the use of Nextera-tagmented reductively-amplified DNA (NextRAD) sequencing to infer movement of a highly-mobile winged insect, the potato psyllid (Bactericera cockerelli), among host plants. The psyllid vectors the pathogen that causes zebra chip disease in potato (Solanum tuberosum), but understanding and managing the spread of this pathogen is limited by uncertainty about the insect’s host plant(s) outside of the growing season. We identified 8,443 polymorphic loci among psyllids separated spatiotemporally on potato or in patches of bittersweet nightshade (S. dulcumara), a weedy plant proposed to be the source of potato-colonizing psyllids. A subset of the psyllids on potato exhibited close genetic similarity to insects on nightshade, consistent with regular movement between these two host plants. However, a second subset of potato-collected psyllids was genetically distinct from those collected on bittersweet nightshade; this suggests that a currently unrecognized host-plant species could be contributing to psyllid populations in potato. Oftentimes, dispersal of vectors of plant or animal pathogens must be tracked at a relatively fine scale in order to understand, predict, and manage disease spread. We demonstrate that emerging sequencing technologies that detect SNPs across a vector’s entire genome can be used to infer such localized movement.

Ranaviruses, a genus of the Iridoviridae, are large double-stranded DNA viruses that infect cold-blooded vertebrates worldwide. Ranaviruses have caused severe epizootics in commercial frog and fish populations, and are currently classified as notifiable pathogens in international trade. Previous work shows that a ranavirus that infects tiger salamanders throughout Western North America (Ambystoma tigrinum virus, or ATV) is in high prevalence among salamanders in the fishing bait trade. Bait ATV strains have elevated virulence and are transported long distances by humans, providing widespread opportunities for pathogen pollution. We sequenced the genomes of 15 strains of ATV collected from tiger salamanders across western North America and performed phylogenetic and population genomic analyses and tests for recombination. We find that ATV forms a monophyletic clade within the rest of the Ranaviruses and that it likely emerged within the last several thousand years, before human activities influenced its spread. We also identify several genes under strong positive selection, some of which appear to be involved in viral virulence and/or host immune evasion. In addition, we provide support for the pathogen pollution hypothesis with evidence of recombination among ATV strains, and potential bait-endemic strain recombination.

We extend a standard model of international risk sharing to include an empirically plausible distortion: Taxes. The tax-inclusive theory implies, even under full risk sharing, a predictable relationship between consumption growth and the consumption and capital income tax rates, both within and across countries. We find strong empirical evidence in favor of this relationship. While idiosyncratic output fluctuations account for substantially more of cross-country consumption growth variability than do taxes, trends in tax differentials are found to be informative about the dynamic evolution of international risk sharing. In particular, adjusting for capital taxes reveals a marked improvement in risk sharing over the last three decades that is absent in baseline measures. This improvement has been driven by the convergence of average tax rates on capital income across OECD countries toward the United States average capital tax rate.

Skill-mismatch employment occurs when high-skilled individuals accept employment in jobs for which they are over-qualified. These employment relationships can be beneficial because they allow high-skilled individuals to more rapidly transition out of unemployment. They come at the cost, however, in the form of lower wage compensation. Moreover, an externality arises as high-skilled individuals do not take into account the effect that their search activity in the market for low-tech jobs has on low-skilled individuals. This paper presents a tractable general equilibrium model featuring mismatch employment and on-the-job search to articulate these tradeoffs. We derive a set of efficiency conditions that describe the labor market distortions associated with these two model features and illustrate how they alter the standard notion of the labor wedges inherent in general equilibrium search models. Finally, we calibrate the model to U.S. data and show that the distortions associated with mismatch employment are largely distributional and can be quantitatively large.

Structural variation, including variation in gene copy number and presence or absence of genes, is a widespread and important source of genomic variation. We used whole-genome DNA sequences from 48 strains of Sinorhizobium (recently renamed Ensifer), including 20 strains of Sinorhizobium meliloti and 12 strains of S. medicae that were the focus of the analyses, to study the fitness effects of new structural variants created by duplication and horizontal gene transfer. We find that derived duplicated and horizontally transferred (HT) genes segregate at lower frequency than synonymous and nonsynonymous nucleotide variants in S. meliloti and S. medicae. Furthermore, the relative frequencies of different types of variants are more similar in S. medicae than in S. meliloti, the species with the larger effective population size. These results are consistent with the hypothesis that most duplications and HT genes have deleterious effects. Diversity of duplications, as measured by segregating duplicated genes per gene, is greater than nucleotide diversity, consistent with a high rate of duplication. Our results suggest that the vast majority of structural variants found among closely related bacterial strains are short-lived and unlikely to be involved in species-wide adaptation.

We study the role of labor market mismatch in the adjustment to a trade liberalization that results in the offshoring of high-tech production. Our model features two-sided heterogeneity in the labor market: high- and low-skilled workers are matched in a frictional labor market with high- and low-tech firms. Mismatch employment occurs when high-skilled workers choose to accept a less desirable job in the low-tech industry. The main result is that perhaps counter-intuitively this type of job displacement is actually beneficial for the labor market in the country doing the offshoring. Mismatch allows the economy to reallocate domestic high-skilled labor across both high- and low-tech industries. In doing so, mismatch dampens both the increase in the aggregate unemployment rate and the decline in aggregate wages that come as a consequence of shifting domestic production abroad.