Plant nucleotide binding/leucine-rich repeat (NLR) immune receptors are activated by pathogen effectors to trigger host defenses and cell death. Toll-interleukin 1 receptor domain NLRs (TNLs) converge on the ENHANCED DISEASE SUSCEPTIBILITY1 (EDS1) family of lipase-like proteins for all resistance outputs. In Arabidopsis (Arabidopsis thaliana) TNL-mediated immunity, AtEDS1 heterodimers with PHYTOALEXIN DEFICIENT4 (AtPAD4) transcriptionally induced basal defenses. AtEDS1 uses the same surface to interact with PAD4-related SENESCENCE-ASSOCIATED GENE101 (AtSAG101), but the role of AtEDS1-AtSAG101 heterodimers remains unclear. We show that AtEDS1-AtSAG101 functions together with N REQUIRED GENE1 (AtNRG1) coiled-coil domain helper NLRs as a coevolved TNL cell death-signaling module. AtEDS1-AtSAG101-AtNRG1 cell death activity is transferable to the Solanaceous species Nicotiana benthamiana and cannot be substituted by AtEDS1-AtPAD4 with AtNRG1 or AtEDS1-AtSAG101 with endogenous NbNRG1. Analysis of EDS1-family evolutionary rate variation and heterodimer structure-guided phenotyping of AtEDS1 variants and AtPAD4-AtSAG101 chimeras identify closely aligned ɑ-helical coil surfaces in the AtEDS1-AtSAG101 partner C-terminal domains that are necessary for reconstituted TNL cell death signaling. Our data suggest that TNL-triggered cell death and pathogen growth restriction are determined by distinctive features of EDS1-SAG101 and EDS1-PAD4 complexes and that these signaling machineries coevolved with other components within plant species or clades to regulate downstream pathways in TNL immunity.

~200 words) 28 Plant intracellular nucleotide-binding/leucine-rich repeat (NLR) immune receptors are 29 activated by pathogen effectors to trigger host defenses and cell death. Toll-30 Interleukin1-receptor (TIR)-domain NLRs (TNLs) converge on the Enhanced Disease 31 Susceptibility1 (EDS1) family of lipase-like proteins for all resistance outputs. In 32 Arabidopsis TNL immunity, AtEDS1 heterodimers with Phytoalexin Deficient4 33 (AtPAD4) transcriptionally boost basal defense pathways. AtEDS1 uses the same 34 surface to interact with PAD4-related Senescence-Associated Gene101 (AtSAG101), 35 but the role of AtEDS1-AtSAG101 heterodimers was unclear. We show that AtEDS1-36 AtSAG101 function together with AtNRG1 coiled-coil domain helper NLRs as a 37 coevolved TNL cell death signaling module. AtEDS1-AtSAG101-AtNRG1 cell death 38 activity is transferable to the solanaceous species, Nicotiana benthamiana, and 39 cannot be substituted by AtEDS1-AtPAD4 with AtNRG1 or AtEDS1-AtSAG101 with 40 endogenous NbNRG1. Analysis of EDS1-family evolutionary rate variation and 41 heterodimer structure-guided phenotyping of AtEDS1 variants or AtPAD4-AtSAG101 42 chimeras identify closely aligned ɑ-helical coil surfaces in the AtEDS1-AtSAG101 43 partner C-terminal domains that are necessary for TNL cell death signaling. Our data 44 suggest that TNL-triggered cell death and pathogen growth restriction are determined 45 by distinctive features of EDS1-SAG101 and EDS1-PAD4 complexes and that these 46 signaling machineries coevolved with further components within plant species or 47 clades to regulate downstream pathways in TNL immunity. 48 49

Background The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus , a seed feeder of the family Lygaeidae. Results The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding. Conclusions With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus ’s strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes. Electronic supplementary material The online version of this article (10.1186/s13059-019-1660-0) contains supplementary material, which is available to authorized users.

Our understanding of how biological diversity arises is limited, especially in the case of speciation in the face of gene flow. Here we investigate the genomic basis of adaptive traits, focusing on a sympatrically diverging species pair of crater lake cichlid fishes. We identify the main quantitative trait loci (QTL) for two eco-morphological traits: body shape and pharyngeal jaw morphology. These traits diverge in parallel between benthic and limnetic species in the repeated adaptive radiations of this and other fish lineages. Remarkably, a single chromosomal region contains the highest effect size QTL for both traits. Transcriptomic data show that the QTL regions contain genes putatively under selection. Independent population genomic data corroborate QTL regions as areas of high differentiation between the sympatric sister species. Our results provide empirical support for current theoretical models that emphasize the importance of genetic linkage and pleiotropy in facilitating rapid divergence in sympatry.

Gene expression levels are important quantitative traits that link genotypes to molecular functions and fitness. In Drosophila, population-genetic studies have revealed substantial adaptive evolution at the genomic level, but the evolutionary modes of gene expression remain controversial. Here, we present evidence that adaptation dominates the evolution of gene expression levels in flies. We show that 64% of the observed expression divergence across seven Drosophila species are adaptive changes driven by directional selection. Our results are derived from time-resolved data of gene expression divergence across a family of related species, using a probabilistic inference method for gene-specific selection. Adaptive gene expression is stronger in specific functional classes, including regulation, sensory perception, sexual behavior, and morphology. Moreover, we identify a large group of genes with sex-specific adaptation of expression, which predominantly occurs in males. Our analysis opens an avenue to map system-wide selection on molecular quantitative traits independently of their genetic basis.

BackgroundThe Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.ResultsThe 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid-nutrition feeding.ConclusionsWith the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

The plant genus Silene has become a model for evolutionary studies of sex chromosomes and sex-determining mechanisms. A recent study performed in Silene colpophylla showed that dioecy and the sex chromosomes in this species evolved independently from those in Silene latifolia, the most widely studied dioecious Silene species. The results of this study show that the sex-determining system in Silene otites, a species related to S. colpophylla, is based on female heterogamety, a sex determination system that is unique among the Silene species studied to date. Our phylogenetic data support the placing of S. otites and S. colpophylla in the subsection Otites and the analysis of ancestral states suggests that the most recent common ancestor of S. otites and S. colpophylla was most probably dioecious. These observations imply that a switch from XX/XY sex determination to a ZZ/ZW system (or vice versa) occurred in the subsection Otites. This is the first report of two different types of heterogamety within one plant genus of this mostly nondioecious plant family. K E Y W O R D S :Evolution, sex chromosomes, sex determination, XY, ZW.

Horizontal gene transfer (HGT) is an important factor in bacterial evolution that can act across species boundaries. Yet, we know little about rate and genomic targets of cross-lineage gene transfer and about its effects on the recipient organism's physiology and fitness. Here, we address these questions in a parallel evolution experiment with two Bacillus subtilis lineages of 7% sequence divergence. We observe rapid evolution of hybrid organisms: gene transfer swaps ∼12% of the core genome in just 200 generations, and 60% of core genes are replaced in at least one population. By genomics, transcriptomics, fitness assays, and statistical modeling, we show that transfer generates adaptive evolution and functional alterations in hybrids. Specifically, our experiments reveal a strong, repeatable fitness increase of evolved populations in the stationary growth phase. By genomic analysis of the transfer statistics across replicate populations, we infer that selection on HGT has a broad genetic basis: 40% of the observed transfers are adaptive. At the level of functional gene networks, we find signatures of negative, positive, and epistatic selection, consistent with hybrid incompatibilities and adaptive evolution of network functions. Our results suggest that gene transfer navigates a complex cross-lineage fitness landscape, bridging epistatic barriers along multiple high-fitness paths.