Relaxed risk of predation drives parallel evolution of stickleback behaviour

The occurrence of similar phenotypes in multiple independent populations (viz. parallel evolution) is a testimony of evolution by natural selection. Parallel evolution implies that populations share a common phenotypic response to a common selection pressure associated with habitat similarity. Examples of parallel evolution at the genetic and phenotypic levels are fairly common, but the driving selective agents often remain elusive. Similarly, the role of phenotypic plasticity in facilitating early stages of parallel evolution is unclear. We investigated whether the relaxation of predation pressure associated with the colonization of freshwater ponds by nine-spined sticklebacks (Pungitius pungitius) likely explains the divergence in complex behaviours between marine and pond populations, and whether this divergence is parallel. Using laboratory-raised individuals exposed to different levels of perceived predation risk, we calculated vectors of phenotypic divergence for four behavioural traits between habitats and predation risk treatments. We found a significant correlation between the directions of evolutionary divergence and phenotypic plasticity, suggesting that habitat divergence in behaviour is aligned with the response to relaxation of predation pressure. Finally, we show that this alignment is found across multiple pairs of populations, and that the relaxation of predation pressure has likely driven parallel evolution of behaviour in this species.

Patterns in the Population History of Northern Eurasia from the Mesolithic to the Early Bronze Age, Based on Craniometry and Genetics

This study examines the craniometric differentiation of Northern Eurasian groups with reference to genetic and partly linguistic facts. Measurements of 66 series of male crania from that territory, dating to various periods from the Mesolithic to the Early Bronze Age, were subjected to statistical methods especially destined for detecting spatial patterns, specifi cally gradients. Using the nonmetric multidimensional scaling of the matrix of D2 distances corrected for sample size, a two-dimensional projection of group constellation was generated, and a minimum spanning tree, showing the shortest path between group centroids in the multivariate space, was constructed. East-west clines in Northern Eurasia, detected by both genetic and craniometric traits, likely indicate not so much gene fl ow as isolation by distance, resulting from an incomplete evolutionary divergence of various fi lial groups constituting the Boreal meta-population. The western fi lial component, which, in Siberia and Eastern Central Asia, is mostly represented by Afanasyevans, has evidently made little contribution to the genetic makeup of later populations. The eastern fi lial component, which had appeared in the Cis-Baikal region from across Lake Baikal no later than the Neolithic, admixed with the autochthonous Paleosiberian component. The latter’s principal marker—the ANE autosomal component—had been present in Siberia since the Upper Paleolithic. Likewise autochthonous were both Eurasian formations—Northern and Southern; statis tical analysis has made it possible to make these more inclusive, whereby the former has been expanded in the eastern direction to include the Kuznetsk Basin, and the latter westwards, to the Middle Irtysh. Nothing suggests that Eastern European groups had taken part in the origin of either the Northern Eurasian formation or the proto-Uralic groups.

Mammalian splicing divergence is shaped by drift, buffering in trans, and a scaling law

Alternative splicing is ubiquitous, but the mechanisms underlying its pattern of evolutionary divergence across mammalian tissues are still underexplored. Here, we investigated the cis-regulatory divergences and their relationship with tissue-dependent trans-regulation in multiple tissues of an F1 hybrid between two mouse species. Large splicing changes between tissues are highly conserved and likely reflect functional tissue-dependent regulation. In particular, micro-exons frequently exhibit this pattern with high inclusion levels in the brain. Cis-divergence of splicing appears to be largely non-adaptive. Although divergence is in general associated with higher densities of sequence variants in regulatory regions, events with high usage of the dominant isoform apparently tolerate more mutations, explaining why their exon sequences are highly conserved but their intronic splicing site flanking regions are not. Moreover, we demonstrate that non-adaptive mutations are often masked in tissues where accurate splicing likely is more important, and experimentally attribute such buffering effect to trans-regulatory splicing efficiency.

Sorghum Association Panel Whole-Genome Sequencing Establishes Pivotal Resource for Dissecting Genomic Diversity

Association mapping panels represent foundational resources for understanding the genetic basis of phenotypic diversity and serve to advance plant breeding by exploring genetic variation across diverse accessions with distinct histories of evolutionary divergence and local adaptation. We report the whole-genome sequencing (WGS) of 400 sorghum [Sorghum bicolor (L.) Moench] accessions from the Sorghum Association Panel (SAP) at an average coverage of 38X (25X-72X), enabling the development of a high-density genomic-marker set of 43,983,694 variants including SNPs (~38 million), indels (~5 million), and CNVs (~170,000). We observe slightly more deletions among indels and a much higher prevalence of deletions among copy number variants compared to insertions. This new marker set enabled the identification of several putatively novel genomic associations for plant height and tannin content, which were not identified when using previous lower-density marker sets. WGS identified and scored variants in 5 kb bins where available genotyping-by-sequencing (GBS) data captured no variants, with half of all bins in the genome falling into this category. The predictive ability of genomic best unbiased linear predictor (GBLUP) models was increased by an average of 30% by using WGS markers rather than GBS markers. We identified 18 selection peaks across subpopulations that formed due to evolutionary divergence during domestication, and we found six Fst peaks resulting from comparisons between converted lines and breeding lines within the SAP that were distinct from the peaks associated with historic selection. This population has been and continues to serve as a significant public resource for sorghum research and demonstrates the value of improving upon existing genomic resources.

Alignment of genetic differentiation across trophic levels in a fig community

Ecological interactions can generate close associations among species, which can in turn generate a high degree of overlap in their spatial distributions. Co-occurrence is likely to be particularly intense when species exhibit obligate comigration, in which they not only overlap in spatial distributions but also travel together from patch to patch. In theory, this pattern of ecological co-occurrence should leave a distinct signature in the pattern of genetic differentiation within and among species. Perhaps the most famous mutual co-isolation partners are fig trees and their co-evolved wasp pollinators. Here, we add another tropic level to this system by examining patterns of genomic diversity in the nematode Caenorhabditis inopinata, a close relative of the C. elegans model system that thrives in figs and obligately disperses on fig wasps. We performed RADseq on individual worms isolated from the field across three Okinawan island populations. The male/female C. inopinata is about five times more diverse than the hermaphroditic C. elegans, and polymorphism is enriched on chromosome arms relative to chromosome centers. FST is low among island population pairs, and clear population structure could not be easily detected among figs, trees, and islands, suggesting frequent migration of wasps between islands. Moreover, inbreeding coefficients are elevated in C. inopinata, consistent with field observations suggesting small C. inopinata founding populations in individual figs. These genetic patterns in C. inopinata overlap with those previously reported in its specific fig wasp vector and are consistent with C. inopinata population dynamics being driven by wasp dispersal. Thus, interspecific interactions can align patterns of genetic diversity across species separated by hundreds of millions of years of evolutionary divergence.

SHED-Dependent Oncogenic Signaling of the PEAK3 Pseudo-Kinase

The PEAK1 and Pragmin/PEAK2 pseudo-kinases have emerged as important components of the protein tyrosine kinase pathway implicated in cancer progression. They can signal using a scaffolding mechanism that involves a conserved split helical dimerization (SHED) module. We recently identified PEAK3 as a novel member of this family based on structural homology; however, its signaling mechanism remains unclear. In this study, we found that, although it can self-associate, PEAK3 shows higher evolutionary divergence than PEAK1/2. Moreover, the PEAK3 protein is strongly expressed in human hematopoietic cells and is upregulated in acute myeloid leukemia. Functionally, PEAK3 overexpression in U2OS sarcoma cells enhanced their growth and migratory properties, while its silencing in THP1 leukemic cells reduced these effects. Importantly, an intact SHED module was required for these PEAK3 oncogenic activities. Mechanistically, through a phosphokinase survey, we identified PEAK3 as a novel inducer of AKT signaling, independent of growth-factor stimulation. Then, proteomic analyses revealed that PEAK3 interacts with the signaling proteins GRB2 and ASAP1/2 and the protein kinase PYK2, and that these interactions require the SHED domain. Moreover, PEAK3 activated PYK2, which promoted PEAK3 tyrosine phosphorylation, its association with GRB2 and ASAP1, and AKT signaling. Thus, the PEAK1-3 pseudo-kinases may use a conserved SHED-dependent mechanism to activate specific signaling proteins to promote oncogenesis.

Outcrossing rates in an experimentally admixed population of self-compatible and self-incompatible Arabidopsis lyrata

The transition to self-compatibility from self-incompatibility is often associated with high rates of self-fertilization, which can restrict gene flow among populations and cause reproductive isolation of self-compatible (SC) lineages. Secondary contact between SC and self-incompatible (SI) lineages might re-establish gene flow if SC lineages remain capable of outcrossing. By contrast, intrinsic features of SC plants that reinforce high rates of self-fertilization could maintain evolutionary divergence between lineages. Arabidopsis lyrata subsp. lyrata is characterized by multiple origins of self-compatibility and high rates of self-fertilization in SC-dominated populations. It is unclear whether these high rates of selfing by SC plants have intrinsic or extrinsic causes. We estimated outcrossing rates and examined patterns of pollinator movement for 38 SC and 40 SI maternal parents sampled from an admixed array of 1509 plants sourced from six SC and six SI populations grown under uniform density. Although plants from SI populations had higher outcrossing rates (mean tm = 0.78 ± 0.05 SE) than plants from SC populations (mean tm = 0.56 ± 0.06 SE), outcrossing rates among SC plants were substantially higher than previous estimates from natural populations. Patterns of pollinator movement appeared to contribute to lower outcrossing rates for SC plants; we estimated that 40% of floral visits were geitonogamous (between flowers of the same plant). The relatively high rates of outcrossing for SC plants under standardized conditions indicate that selfing rates in natural SC populations of A. lyrata are facultative and driven by extrinsic features of A. lyrata, including patterns of pollinator movement.

Micro- and macroevolutionary change of colony-level traits in bryozoans

The evolution of trait variation among populations of animals is difficult to study due to the many overlapping genetic and environmental influences that control phenotypic expression. In a group of animals, bryozoans, it is possible to isolate genetic contributions to phenotypic variation, due to the modular nature of bryozoan colonies. Each bryozoan colony represents a snapshot of the phenotypes that correspond to a single genotype, which can be summarized as a phenotypic distribution. We test whether these phenotypic distributions are heritable across generations of colonies in two sister species of the bryozoan Stylopoma, grown and bred in a common garden breeding experiment. We find that components of phenotypic distributions, specifically median trait values of colony members, are heritable between generations of colonies. Furthermore, this heredity has macroevolutionary importance because it correlates with the morphological distance between these two species. Because parts of phenotypic distributions are heritable, and this heritability corresponds to evolutionary divergence between species, we infer that these distributions have the potential to evolve. The evolutionary potential of these phenotypic distributions may underpin the emergence of colony-level traits, like division of labor in colonies.

Complex polyploid and hybrid species in an apomictic and sexual tropical forage grass group: genomic composition and evolution in Urochloa (Brachiaria) species

Background and Aims Diploid and polyploid Urochloa (including Brachiaria, Panicum and Megathyrsus species) C4 tropical forage grasses originating from Africa are important for food security and the environment​, often being planted in marginal lands worldwide. We aimed to characterize the nature of their genomes, the repetitive DNA, and the genome composition of polyploids, leading to a model of the evolutionary pathways within the group including many apomictic species. Methods Some 362 forage grass accessions from international germplasm collections were studied, and ploidy determined using an optimized flow cytometry method. Whole-genome survey sequencing and molecular cytogenetic analysis were used to identify chromosomes and genomes in Urochloa accessions belonging to the 'brizantha' and 'humidicola' agamic complexes and U. maxima. Key Results Genome structures are complex and variable, with multiple ploidies and genome compositions within the species, and no clear geographical patterns. Sequence analysis of nine diploid and polyploid accessions enabled identification of abundant genome-specific repetitive DNA motifs​. In situ hybridization with a combination of repetitive DNA and genomic DNA probes, identified evolutionary divergence and allowed us to discriminate the different genomes present in polyploids. Conclusions We suggest a new coherent nomenclature for the genomes present​. We develop a model of evolution at the whole-genome level in diploid and polyploid accessions showing processes of grass evolution. We support the retention of narrow species concepts for U. brizantha, U. decumbens, and U. ruziziensis, and do not consider diploids and polyploids of single species as cytotypes. The results and model will be valuable in making rational choices of parents for new hybrids, assist in use of the germplasm for breeding and selection of Urochloa with improved sustainability and agronomic potential, and will assist in measuring and conserving biodiversity in grasslands.

Transcriptomic Response to Perkinsus marinus in Two Crassostrea Oysters Reveals Evolutionary Dynamics of Host-Parasite Interactions

Infectious disease outbreaks are causing widespread declines of marine invertebrates including corals, sea stars, shrimps, and molluscs. Dermo is a lethal infectious disease of the eastern oyster Crassostrea virginica caused by the protist Perkinsus marinus. The Pacific oyster Crassostrea gigas is resistant to Dermo due to differences in the host-parasite interaction that is not well understood. We compared transcriptomic responses to P. marinus challenge in the two oysters at early and late infection stages. Dynamic and orchestrated regulation of large sets of innate immune response genes were observed in both species with remarkably similar patterns for most orthologs, although responses in C. virginica were stronger, suggesting strong or over-reacting immune response could be a cause of host mortality. Between the two species, several key immune response gene families differed in their expansion, sequence variation and/or transcriptional response to P. marinus, reflecting evolutionary divergence in host-parasite interaction. Of note, significant upregulation of inhibitors of apoptosis (IAPs) was observed in resistant C. gigas but not in susceptible C. virginica, suggesting upregulation of IAPs is an active defense mechanism, not a passive response orchestrated by P. marinus. Compared with C. gigas, C. virginica exhibited greater expansion of toll-like receptors (TLRs) and positive selection in P. marinus responsive TLRs. The C1q domain containing proteins (C1qDCs) with the galactose-binding lectin domain that is involved in P. marinus recognition, were only present and significantly upregulated in C. virginica. These results point to previously undescribed differences in host defense genes between the two oyster species that may account for the difference in susceptibility, providing an expanded portrait of the evolutionary dynamics of host-parasite interaction in lophotrochozoans that lack adaptive immunity. Our findings suggest that C. virginica and P. marinus have a history of coevolution and the recent outbreaks may be due to increased virulence of the parasite.