Adaptive evolution of hybrid bacteria by horizontal gene transfer.

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.

A Coevolved EDS1-SAG101-NRG1 Module Mediates Cell Death Signaling by TIR-Domain Immune Receptors.

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.

Possible mechanisms responsible for absence of a retrotransposon family on a plant Y chromosome.

Some transposable elements (TEs) show extraordinary variance in abundance along sex chromosomes but the mechanisms responsible for this variance are unknown. Here, we studied Ogre long terminal repeat (LTR) retrotransposons in Silene latifolia, a dioecious plant with evolutionarily young heteromorphic sex chromosomes. Ogre elements are ubiquitous in the S. latifolia genome but surprisingly absent on the Y chromosome. Bacterial artificial chromosome (BAC) library analysis and fluorescence in situ hybridization (FISH) were used to determine Ogre structure and chromosomal localization. Next generation sequencing (NGS) data were analysed to assess the transcription level and abundance of small RNAs. Methylation of Ogres was determined by bisulphite sequencing. Phylogenetic analysis was used to determine mobilization time and selection forces acting on Ogre elements. We characterized three Ogre families ubiquitous in the S. latifolia genome. One family is nearly absent on the Y chromosome despite all the families having similar structures and spreading mechanisms. We showed that Ogre retrotransposons evolved before sex chromosomes appeared but were mobilized after formation of the Y chromosome. Our data suggest that the absence of one Ogre family on the Y chromosome may be caused by 24-nucleotide (24-nt) small RNA-mediated silencing leading to female-specific spreading. Our findings highlight epigenetic silencing mechanisms as potentially crucial factors in sex-specific spreading of some TEs, but other possible mechanisms are also discussed.

Vivid COVID-19 LAMP is an ultrasensitive, quadruplexed test using LNA-modified primers and a zinc ion and 5-Br-PAPS colorimetric detection system.

Sensitive and rapid point-of-care assays have been crucial in the global response to SARS-CoV-2. Loop-mediated isothermal amplification (LAMP) has emerged as an important diagnostic tool given its simplicity and minimal equipment requirements, although limitations exist regarding sensitivity and the methods used to detect reaction products. We describe the development of Vivid COVID-19 LAMP, which leverages a metallochromic detection system utilizing zinc ions and a zinc sensor, 5-Br-PAPS, to circumvent the limitations of classic detection systems dependent on pH indicators or magnesium chelators. We make important strides in improving RT-LAMP sensitivity by establishing principles for using LNA-modified LAMP primers, multiplexing, and conducting extensive optimizations of reaction parameters. To enable point-of-care testing, we introduce a rapid sample inactivation procedure without RNA extraction that is compatible with self-collected, non-invasive gargle samples. Our quadruplexed assay (targeting E, N, ORF1a, and RdRP) reliably detects 1 RNA copy/µl of sample (=8 copies/reaction) from extracted RNA and 2 RNA copies/µl of sample (=16 copies/reaction) directly from gargle samples, making it one of the most sensitive RT-LAMP tests and even comparable to RT-qPCR. Additionally, we demonstrate a self-contained, mobile version of our assay in a variety of high-throughput field testing scenarios on nearly 9,000 crude gargle samples. Vivid COVID-19 LAMP can be an important asset for the endemic phase of COVID-19 as well as preparing for future pandemics.

Bisprimer--a program for the design of primers for bisulfite-based genomic sequencing of both plant and Mammalian DNA samples.

Plants and animals differ in the sequence context of the methylated sites in DNA. Plants exhibit cytosine methylation in CG, CHG, and CHH sites, whereas CG methylation is the only form present in mammals (with an exception of the early embryonic development). This fact must be taken into account in the design of primers for bisulfite-based genomic sequencing because CHG and CHH sites can remain unmodified. Surprisingly, no user-friendly primer design program is publicly available that could be used to design primers in plants and to simultaneously check the properties of primers such as the potential for primer-dimer formation. For studies concentrating on particular DNA loci, the correct design of primers is crucial. The program, called BisPrimer, includes 2 different subprograms for the primer design, the first one for mammals and the second one for angiosperm plants. Each subprogram is divided into 2 variants. The first variant serves to design primers that preferentially bind to the bisulfite-modified primer-binding sites (C to U conversion). This type of primer preferentially amplifies the bisulfite-converted DNA strands. This feature can help to avoid problems connected with an incomplete bisulfite modification that can sometimes occur for technical reasons. The second variant is intended for the analysis of samples that are supposed to consist of a mixture of DNA molecules that have different levels of cytosine methylation (e.g., pollen DNA). In this case, the aim is to minimize the selection in favor of either less methylated or more methylated molecules.

Sequential development of several RT-qPCR tests using LNA nucleotides and dual probe technology to differentiate SARS-CoV-2 from influenza A and B.

Sensitive and accurate RT-qPCR tests are the primary diagnostic tools to identify SARS-CoV-2-infected patients. While many SARS-CoV-2 RT-qPCR tests are available, there are significant differences in test sensitivity, workflow (e.g. hands-on-time), gene targets and other functionalities that users must consider. Several publicly available protocols shared by reference labs and public health authorities provide useful tools for SARS-CoV-2 diagnosis, but many have shortcomings related to sensitivity and laborious workflows. Here, we describe a series of SARS-CoV-2 RT-qPCR tests that are originally based on the protocol targeting regions of the RNA-dependent RNA polymerase (RdRp) and envelope (E) coding genes developed by the Charité Berlin. We redesigned the primers/probes, utilized locked nucleic acid nucleotides, incorporated dual probe technology and conducted extensive optimizations of reaction conditions to enhance the sensitivity and specificity of these tests. By incorporating an RNase P internal control and developing multiplexed assays for distinguishing SARS-CoV-2 and influenza A and B, we streamlined the workflow to provide quicker results and reduced consumable costs. Some of these tests use modified enzymes enabling the formulation of a room temperature-stable master mix and lyophilized positive control, thus increasing the functionality of the test and eliminating cold chain shipping and storage. Moreover, a rapid, RNA extraction-free version enables high sensitivity detection of SARS-CoV-2 in about an hour using minimally invasive, self-collected gargle samples. These RT-qPCR assays can easily be implemented in any diagnostic laboratory and can provide a powerful tool to detect SARS-CoV-2 and the most common seasonal influenzas during the vaccination phase of the pandemic.

Epigenetic Modulating Chemicals Significantly Affect the Virulence and Genetic Characteristics of the Bacterial Plant Pathogen Xanthomonas campestris pv. campestris.

Epigenetics is the study of heritable alterations in phenotypes that are not caused by changes in DNA sequence. In the present study, we characterized the genetic and phenotypic alterations of the bacterial plant pathogen Xanthomonas campestris pv. campestris (Xcc) under different treatments with several epigenetic modulating chemicals. The use of DNA demethylating chemicals unambiguously caused a durable decrease in Xcc bacterial virulence, even after its reisolation from infected plants. The first-time use of chemicals to modify the activity of sirtuins also showed some noticeable results in terms of increasing bacterial virulence, but this effect was not typically stable. Changes in treated strains were also confirmed by using methylation sensitive amplification (MSAP), but with respect to registered SNPs induction, it was necessary to consider their contribution to the observed polymorphism. The molecular basis of the altered virulence was deciphered by using dualRNA-seq analysis of treated Xcc strains infecting Brassica rapa plants. The results of the present study should promote more intensive research in the generally understudied field of bacterial epigenetics, where artificially induced modification by epigenetic modulating chemicals can significantly increase the diversity of bacterial properties and potentially contribute to the further development of the fields, such as bacterial ecology and adaptation.

Surveillance of SARS-CoV-2 lineage B.1.1.7 in Slovakia using a novel, multiplexed RT-qPCR assay.

The emergence of a novel SARS-CoV-2 B.1.1.7 variant sparked global alarm due to increased transmissibility, mortality, and uncertainty about vaccine efficacy, thus accelerating efforts to detect and track the variant. Current approaches to detect B.1.1.7 include sequencing and RT-qPCR tests containing a target assay that fails or results in reduced sensitivity towards the B.1.1.7 variant. Since many countries lack genomic surveillance programs and failed assays detect unrelated variants containing similar mutations as B.1.1.7, we used allele-specific PCR, and judicious placement of LNA-modified nucleotides to develop an RT-qPCR test that accurately and rapidly differentiates B.1.1.7 from other SARS-CoV-2 variants. We validated the test on 106 clinical samples with lineage status confirmed by sequencing and conducted a country-wide surveillance study of B.1.1.7 prevalence in Slovakia. Our multiplexed RT-qPCR test showed 97% clinical sensitivity and retesting 6,886 SARS-CoV-2 positive samples obtained during three campaigns performed within one month, revealed pervasive spread of B.1.1.7 with an average prevalence of 82%. Labs can easily implement this test to rapidly scale B.1.1.7 surveillance efforts and it is particularly useful in countries with high prevalence of variants possessing only the ΔH69/ΔV70 deletion because current strategies using target failure assays incorrectly identify these as putative B.1.1.7 variants.

Discovery of a Family of Mixed Lineage Kinase Domain-like Proteins in Plants and Their Role in Innate Immune Signaling.

HeLo domain-containing mixed lineage kinase domain-like protein (MLKL), a pseudokinase, mediates necroptotic cell death in animals. Here, we report the discovery of a conserved protein family across seed plants that structurally resembles vertebrate MLKL. The Arabidopsis genome encodes three MLKLs (AtMLKLs) with overlapping functions in disease resistance mediated by Toll-interleukin 1-receptor domain intracellular immune receptors (TNLs). The HeLo domain of AtMLKLs confers cell death activity but is dispensable for immunity. Cryo-EM structures reveal a tetrameric configuration, in which the HeLo domain is buried, suggestive of an auto-repressed complex. The mobility of AtMLKL1 along microtubules is reduced by chitin, a fungal immunity-triggering molecule. An AtMLKL1 phosphomimetic variant exhibiting reduced mobility enhances immunity. Coupled with the predicted presence of HeLo domains in plant helper NLRs, our data reveal the importance of HeLo domain proteins for TNL-dependent immunity and argue for a cell death-independent immune mechanism mediated by MLKLs.

Evolution of sex determination and heterogamety changes in section Otites of the genus Silene.

Switches in heterogamety are known to occur in both animals and plants. Although plant sex determination systems probably often evolved more recently than those in several well-studied animals, including mammals, and have had less time for switches to occur, we previously detected a switch in heterogamety in the plant genus Silene: section Otites has both female and male heterogamety, whereas S. latifolia and its close relatives, in a different section of the genus, Melandrium (subgenus Behenantha), all have male heterogamety. Here we analyse the evolution of sex chromosomes in section Otites, which is estimated to have evolved only about 0.55 MYA. Our study confirms female heterogamety in S. otites and newly reveals female heterogamety in S. borysthenica. Sequence analyses and genetic mapping show that the sex-linked regions of these two species are the same, but the region in S. colpophylla, a close relative with male heterogamety, is different. The sex chromosome pairs of S. colpophylla and S. otites each correspond to an autosome of the other species, and both differ from the XY pair in S. latifolia. Silene section Otites species are suitable for detailed studies of the events involved in such changes, and our phylogenetic analysis suggests a possible change from female to male heterogamety within this section. Our analyses suggest a possibility that has so far not been considered, change in heterogamety through hybridization, in which a male-determining chromosome from one species is introgressed into another one, and over-rides its previous sex-determining system.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome.

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.

Adaptive Evolution of Gene Expression in Drosophila.

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.

Genetic linkage of distinct adaptive traits in sympatrically speciating crater lake cichlid fish.

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.

The evolutionary fate of the horizontally transferred agrobacterial mikimopine synthase gene in the genera Nicotiana and Linaria.

Few cases of spontaneously horizontally transferred bacterial genes into plant genomes have been described to date. The occurrence of horizontally transferred genes from the T-DNA of Agrobacterium rhizogenes into the plant genome has been reported in the genus Nicotiana and in the species Linaria vulgaris. Here we compare patterns of evolution in one of these genes (a gene encoding mikimopine synthase, mis) following three different events of horizontal gene transfer (HGT). As this gene plays an important role in Agrobacterium, and there are known cases showing that genes from pathogens can acquire plant protection function, we hypothesised that in at least some of the studied species we will find signs of selective pressures influencing mis sequence. The mikimopine synthase (mis) gene evolved in a different manner in the branch leading to Nicotiana tabacum and N. tomentosiformis, in the branch leading to N. glauca and in the genus Linaria. Our analyses of the genus Linaria suggest that the mis gene began to degenerate soon after the HGT. In contrast, in the case of N. glauca, the mis gene evolved under significant selective pressures. This suggests a possible role of mikimopine synthase in current N. glauca and its ancestor(s). In N. tabacum and N. tomentosiformis, the mis gene has a common frameshift mutation that disrupted its open reading frame. Interestingly, our results suggest that in spite of the frameshift, the mis gene could evolve under selective pressures. This sequence may still have some regulatory role at the RNA level as suggested by coverage of this sequence by small RNAs in N. tabacum.

Characterization of the HMA7 gene and transcriptomic analysis of candidate genes for copper tolerance in two Silene vulgaris ecotypes.

Silene vulgaris possesses ecotype-specific tolerance to high levels of copper in the soil. Although this was reported a few decades ago, little is known about this trait on a molecular level. The aim of this study was to analyze the transcription response to elevated copper concentrations in two S. vulgaris ecotypes originating from copper-contrasting soil types - copper-tolerant Lubietova and copper-sensitive Stranska skala. To reveal if plants are transcriptionally affected, we first analyzed the HMA7 gene, a known key player in copper metabolism. Based on BAC library screening, we identified a BAC clone containing a SvHMA7 sequence with all the structural properties specific for plant copper-transporting ATPases. The functionality of the gene was tested using heterologous complementation in yeast mutants. Analyses of SvHMA7 transcription patterns showed that both ecotypes studied up-regulated SvHMA7 transcription after the copper treatment. Our data are supported by analysis of appropriate reference genes based on RNA-Seq databases. To identify genes specifically involved in copper response in the studied ecotypes, we analyzed transcription profiles of genes coding Cu-transporting proteins and genes involved in the prevention of copper-induced oxidative stress in both ecotypes. Our data show that three genes (APx, POD and COPT5) differ in their transcription pattern between the ecotypes with constitutively increased transcription in Lubietova. Taken together, we have identified transcription differences between metallifferous and non-metalliferous ecotypes of S. vulgaris, and we have suggested candidate genes participating in metal tolerance in this species.

Evolution of sex determination systems with heterogametic males and females in silene.

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.