Positive selection in cilia-related genes may facilitate deep-sea adaptation of Thermocollonia jamsteci

Deep-sea hydrothermal vents are characterized by high hydrostatic pressure, hypoxia, darkness and toxic substances. However, how organisms adapt to such extreme marine ecosystems remain poorly understood. We hypothesize that adaptive evolution plays an essential role in generating novelty for evolutionary adaptation to the deep-sea environment because adaptive evolution has been found to be critical for species origin and evolution. In this project, the chromosome-level genome of the deep-sea hydrothermal vent gastropod T. jamsteci was constructed for the first time to examine molecular mechanisms of its adaptation to the deep-sea environment. The genome size was large (2.54 Gb), ranking at the top of all species in the Vetigastropoda subclass, driven primarily by the bursts of transposable elements (TEs). The transposition of TEs may also trigger chromosomal changes including both inter-chromosomal fusions and intra-chromosomal activities involving chromosome inversions, rearrangements and fusions, as revealed by comparing the genomes of T. jamsteci and its closely related shallow-sea species Gibbula magus. Innovative changes including the expansion of the ABC transporter gene family that may facilitate detoxification, duplication of genes related to endocytosis, immunity, apoptosis, and anti-apoptotic domains that may help T. jamsteci fight against microbial pathogens, were identified. Furthermore, comparative analysis identified positive selection signals in a large number of genes including the hypoxia up-regulated protein 1, which is a chaperone that may promote adaptation of the T. jamsteci to hypoxic deepsea environments, hox2, Rx2, Pax6 and cilia-related genes BBS1, BBS2, BBS9 and RFX4. Notably, because of the critical importance of cilia and IFT in development, positive selection in cilia-related genes may play a critical role in facilitating T. jamsteci to adapt to the high-pressure deep-sea ecosystem. Results from this study thus revealed important molecular clues that may facilitate further research on the adaptation of molluscs to deep-sea hydrothermal vents.

Myzus persicae resistance to neonicotinoids-unravelling the contribution of different mechanisms to phenotype

BACKGROUND

Deciphering the mechanisms underlying insecticide resistance is key to devising appropriate strategies against this economically important trait. Myzus persicae, the green peach-potato aphid, is a major pest that has evolved resistance to many insecticide classes, including neonicotinoids. M. persicae resistance to neonicotinoids has previously been shown to result from two main mechanisms: metabolic resistance resulting from P450 overexpression and a targetsite mutation, R81T. However, their respective contribution to resistant phenotypes remains unclear.

RESULTS

By combining extensive insecticide bioassays with and without addition of the synergist PBO, and gene copy number and expression quantification of two key P450 enzymes (CYP6CY3 and CYP6CY4) in a 23 clone collection, we, (i) confirmed that metabolic resistance is correlated with P450 expression level, up to a threshold, (ii) demonstrated that the R81T mutation, in the homozygous state and in combination with P450 overexpression, leads to high levels of resistance to neonicotinoids, and, (iii) showed that there is a synergistic interaction between the P450 and R81T mechanisms, and that this interaction has the strongest impact on the strength of resistance phenotypes. However, even though the R81T mutation has a great effect on the resistance phenotype, different R81T genotypes can exhibit variation in the level of resistance, explained only partially by P450 overexpression.

CONCLUSION

To comprehend resistance phenotypes, it is important to take into account every mechanism at play, as well as the way these mechanisms interact. © 2024 The Author(s). Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Plasmid-encoded insertion sequences promote rapid adaptation in clinical enterobacteria

Plasmids are extrachromosomal genetic elements commonly found in bacteria. They are known to fuel bacterial evolution through horizontal gene transfer, and recent analyses indicate that they can also promote intragenomic adaptations. However, the role of plasmids as catalysts of bacterial evolution beyond horizontal gene transfer is poorly explored. In this study, we investigated the impact of a widespread conjugative plasmid, pOXA-48, on the evolution of several multidrug-resistant clinical enterobacteria. Combining experimental and within-patient evolution analyses, we unveiled that plasmid pOXA-48 promotes bacterial evolution through the transposition of plasmid-encoded insertion sequence 1 (IS1) elements. Specifically, IS1-mediated gene inactivation expedites the adaptation rate of clinical strains in vitro and fosters within-patient adaptation in the gut microbiota. We deciphered the mechanism underlying the plasmid-mediated surge in IS1 transposition, revealing a negative feedback loop regulated by the genomic copy number of IS1. Given the overrepresentation of IS elements in bacterial plasmids, our findings suggest that plasmid-mediated IS1 transposition represents a crucial mechanism for swift bacterial adaptation.

Genetic Variation in Jamaican Populations of the Coffee Berry Borer, Hypothenemus hampei

The coffee berry borer (CBB) Hypothenemus hampei was first described in Africa in 1867 and has spread to all major coffee-producing regions worldwide, including Jamaica. Using long-read sequencing, we produced a new high-quality reference genome (172.7 Mb) for the Jamaican strain of the CBB, with 93% of the genome assembled into 14 scaffolds. Whole genome sequencing of pooled samples from different populations across Jamaica showed that the CBB harbors low levels of genetic diversity alongside an excess of low-frequency alleles, indicative of a recent genetic bottleneck. The analyses also showed a recent surge in the activity of transposable elements (TEs), particularly LINE/R1 and LTR/Gypsy elements, within CBB populations. Our findings offer first insights into the evolutionary genomics of CBB populations in Jamaica, highlighting the potential role of TEs in shaping the genome of this important pest species.

Tyzzerella nexilis strains enriched in mobile genetic elements are involved in progressive multiple sclerosis

Multiple sclerosis (MS) is an autoimmune-demyelinating disease with an inflammatory pathology formed by self-reactive lymphocytes with activated glial cells. Progressive MS, characterized by resistance to medications, significantly differs from the non-progressive form in gut microbiome profiles. After confirming an increased abundance of "Tyzzerella nexilis" in various cohorts of progressive MS, we identified a distinct cluster of T. nexilis strains enriched in progressive MS based on long-read metagenomics. The distinct T. nexilis cluster is characterized by a large number of mobile genetic elements (MGEs) and a lack of defense systems against MGEs. Microbial genes for sulfate reduction and flagella formation with pathogenic implications are specific to this cluster. Moreover, these flagellar genes are encoded on MGEs. Mono-colonization with MGE-enriched T. nexilis made germ-free mice more susceptible to experimental autoimmune encephalomyelitis. These results indicate that the progression of MS may be promoted by MGE-enriched T. nexilis with potentially pathogenic properties.

Comparative analysis of amphibian genomes: An emerging resource for basic and applied research

Amphibians are the most threatened group of vertebrates and are in dire need of conservation intervention to ensure their continued survival. They exhibit unique features including a high diversity of reproductive strategies, permeable and specialized skin capable of producing toxins and antimicrobial compounds, multiple genetic mechanisms of sex determination and in some lineages, the ability to regenerate limbs and organs. Although genomic approaches would shed light on these unique traits and aid conservation, sequencing and assembly of amphibian genomes has lagged behind other taxa due to their comparatively large genome sizes. Fortunately, the development of long-read sequencing technologies and initiatives has led to a recent burst of new amphibian genome assemblies. Although growing, the field of amphibian genomics suffers from the lack of annotation resources, tools for working with challenging genomes and lack of high-quality assemblies in multiple clades of amphibians. Here, we analyse 51 publicly available amphibian genomes to evaluate their usefulness for functional genomics research. We report considerable variation in genome assembly quality and completeness and report some of the highest transposable element and repeat contents of any vertebrate. Additionally, we detected an association between transposable element content and climatic variables. Our analysis provides evidence of conserved genome synteny despite the long divergence times of this group, but we also highlight inconsistencies in chromosome naming and orientation across genome assemblies. We discuss sequencing gaps in the phylogeny and suggest key targets for future sequencing endeavours. Finally, we propose increased investment in amphibian genomics research to promote their conservation.

Mobilome impacts on physiology in the widely used non-toxic mutant Microcystis aeruginosa PCC 7806 ΔmcyB and toxic wildtype

The Microcystis mobilome is a well-known but understudied component of this bloom-forming cyanobacterium. Through genomic and transcriptomic comparisons, we found five families of transposases that altered the expression of genes in the well-studied toxigenic type-strain, Microcystis aeruginosa PCC 7086, and a non-toxigenic genetic mutant, Microcystis aeruginosa PCC 7806 ΔmcyB. Since its creation in 1997, the ΔmcyB strain has been used in comparative physiology studies against the wildtype strain by research labs throughout the world. Some differences in gene expression between what were thought to be otherwise genetically identical strains have appeared due to insertion events in both intra- and intergenic regions. In our ΔmcyB isolate, a sulfate transporter gene cluster (sbp-cysTWA) showed differential expression from the wildtype, which may have been caused by the insertion of a miniature inverted repeat transposable element (MITE) in the sulfate-binding protein gene (sbp). Differences in growth in sulfate-limited media also were also observed between the two isolates. This paper highlights how Microcystis strains continue to "evolve" in lab conditions and illustrates the importance of insertion sequences / transposable elements in shaping genomic and physiological differences between Microcystis strains thought otherwise identical. This study forces the necessity of knowing the complete genetic background of isolates in comparative physiological experiments, to facilitate the correct conclusions (and caveats) from experiments.

Agricultural insect pests as models for studying stress-induced evolutionary processes

Agricultural insect pests (AIPs) are widely successful in adapting to natural and anthropogenic stressors, repeatedly overcoming population bottlenecks and acquiring resistance to intensive management practices. Although they have been largely overlooked in evolutionary studies, AIPs are ideal systems for understanding rapid adaptation under novel environmental conditions. Researchers have identified several genomic mechanisms that likely contribute to adaptive stress responses, including positive selection on de novo mutations, polygenic selection on standing allelic variation and phenotypic plasticity (e.g., hormesis). However, new theory suggests that stress itself may induce epigenetic modifications, which may confer heritable physiological changes (i.e., stress-resistant phenotypes). In this perspective, we discuss how environmental stress from agricultural management generates the epigenetic and genetic modifications that are associated with rapid adaptation in AIPs. We summarise existing evidence for stress-induced evolutionary processes in the context of insecticide resistance. Ultimately, we propose that studying AIPs offers new opportunities and resources for advancing our knowledge of stress-induced evolution.

Transposable elements in Drosophila montana from harsh cold environments

BACKGROUND

Substantial discoveries during the past century have revealed that transposable elements (TEs) can play a crucial role in genome evolution by affecting gene expression and inducing genetic rearrangements, among other molecular and structural effects. Yet, our knowledge on the role of TEs in adaptation to extreme climates is still at its infancy. The availability of long-read sequencing has opened up the possibility to identify and study potential functional effects of TEs with higher precision. In this work, we used Drosophila montana as a model for cold-adapted organisms to study the association between TEs and adaptation to harsh climates.

RESULTS

Using the PacBio long-read sequencing technique, we de novo identified and manually curated TE sequences in five Drosophila montana genomes from eco-geographically distinct populations. We identified 489 new TE consensus sequences which represented 92% of the total TE consensus in D. montana. Overall, 11-13% of the D. montana genome is occupied by TEs, which as expected are non-randomly distributed across the genome. We identified five potentially active TE families, most of them from the retrotransposon class of TEs. Additionally, we found TEs present in the five analyzed genomes that were located nearby previously identified cold tolerant genes. Some of these TEs contain promoter elements and transcription binding sites. Finally, we detected TEs nearby fixed and polymorphic inversion breakpoints.

CONCLUSIONS

Our research revealed a significant number of newly identified TE consensus sequences in the genome of D. montana, suggesting that non-model species should be studied to get a comprehensive view of the TE repertoire in Drosophila species and beyond. Genome annotations with the new D. montana library allowed us to identify TEs located nearby cold tolerant genes, and present at high population frequencies, that contain regulatory regions and are thus good candidates to play a role in D. montana cold stress response. Finally, our annotations also allow us to identify for the first time TEs present in the breakpoints of three D. montana inversions.

Retrotransposon-driven environmental regulation of FLC leads to adaptive response to herbicide

The mobilization of transposable elements is a potent source of mutations. In plants, several stransposable elements respond to external cues, fuelling the hypothesis that natural transposition can create environmentally sensitive alleles for adaptation. Here we report on the detailed characterization of a retrotransposon insertion within the first intron of the Arabidopsis floral-repressor gene FLOWERING LOCUS C (FLC) and the discovery of its role for adaptation. The insertion mutation augments the environmental sensitivity of FLC by affecting the balance between coding and non-coding transcripts in response to stress, thus expediting flowering. This balance is modulated by DNA methylation and orchestrated by IBM2, a factor involved in the processing of intronic heterochromatic sequences. The stress-sensitive allele of FLC has spread across populations subjected to recurrent chemical weeding, and we show that retrotransposon-driven acceleration of the life cycle represents a rapid response to herbicide application. Our work provides a compelling example of a transposable element-driven environmentally sensitive allele that confers an adaptive response in nature.

Metagenomics of high-altitude groundwater reveal different health risks associated with antibiotic-resistant pathogens and bacterial resistome in the latitudinal gradient

Groundwater on the Tibetan Plateau is a critical water resource to people in Asia. However, its prevalence of antibiotic-resistant pathogens (ARPs), bacterial resistome and their driving factors remain unknown. Using metagenomics analysis, a hotspot of antibiotic-resistance genes (ARGs) and last-resort ARGs (LARGs) with a total of 639 subtypes was identified in the groundwater. Importantly, 164 metagenome-assembled genomes (MAGs) which possessed both ARGs and virulence factors (VFs) were assigned as potential ARPs, with the most abundant species being Acinetobacter johnsonii and Acinetobacter pittii. A total of 157 potential ARPs, involving Escherichia coli, were predicted as "natural" ARGs supercarriers. Thirty-six ARPs dominated by the genus Acinetobacter and Pseudomonas were found to harbour LARGs. Co-localizations of the ARG-mobile genetic elements (MGEs) showed that MGEs were significantly associated with ARGs in the ARPs, which suggests ARPs play a prominent role in ARG dissemination. Notably, latitudinal gradient is a driving factor in the occurrence of ARGs and ARPs. The average abundances of ARGs and ARP decreased as the latitude increased, with the highest abundance occurring in the region between 28.6◦N and 29.5◦N. MetaCompare further revealed health risks associated with the resistome decreased as the latitudes increased. These findings indicated different health risks associated with ARPs and bacterial resistome in latitudinal gradient groundwater. They raise the concerns of mitigating ARPs risk in groundwater on the Tibetan Plateau.

Long-term methylome changes after experimental seed demethylation and their interaction with recurrent water stress in Erodium cicutarium (Geraniaceae)

The frequencies and lengths of drought periods are increasing in subtropical and temperate regions worldwide. Epigenetic responses to water stress could be key for plant resilience to these largely unpredictable challenges. Experimental DNA demethylation, together with application of a stress factor is an appropriate strategy to reveal the contribution of epigenetics to plant responses to stress. We analysed leaf cytosine methylation changes in adult plants of the annual Mediterranean herb, Erodium cicutarium, in a greenhouse, after seed demethylation with 5-Azacytidine and/or recurrent water stress. We used bisulfite RADseq (BsRADseq) and a newly reported reference genome for E. cicutarium to characterize methylation changes in a 2 × 2 factorial design, controlling for plant relatedness. In the long term, 5-Azacytidine treatment alone caused both hypo- and hyper-methylation at individual cytosines, with substantial hypomethylation in CG contexts. In control conditions, drought resulted in a decrease in methylation in all but CHH contexts. In contrast, the genome of plants that experienced recurrent water stress and had been treated with 5-Azacytidine increased DNA methylation level by ca. 5%. Seed demethylation and recurrent drought produced a highly significant interaction in terms of global and context-specific cytosine methylation. Most methylation changes occurred around genic regions and within Transposable Elements. The annotation of these Differentially Methylated Regions associated with genes included several with a potential role in stress responses (e.g., PAL, CDKC, and ABCF), confirming an epigenetic contribution in response to stress at the molecular level.

Fundamental Patterns of Structural Evolution Revealed by Chromosome-Length Genomes of Cactophilic Drosophila

A thorough understanding of adaptation and speciation requires model organisms with both a history of ecological and phenotypic study as well as a complete set of genomic resources. In particular, high-quality genome assemblies of ecological model organisms are needed to assess the evolution of genome structure and its role in adaptation and speciation. Here, we generate new genomes of cactophilic Drosophila, a crucial model clade for understanding speciation and ecological adaptation in xeric environments. We generated chromosome-level genome assemblies and complete annotations for seven populations across Drosophila mojavensis, Drosophila arizonae, and Drosophila navojoa. We use these data first to establish the most robust phylogeny for this clade to date, and to assess patterns of molecular evolution across the phylogeny, showing concordance with a priori hypotheses regarding adaptive genes in this system. We then show that structural evolution occurs at constant rate across the phylogeny, varies by chromosome, and is correlated with molecular evolution. These results advance the understanding of the D. mojavensis clade by demonstrating core evolutionary genetic patterns and integrating those patterns to generate new gene-level hypotheses regarding adaptation. Our data are presented in a new public database (cactusflybase.arizona.edu), providing one of the most in-depth resources for the analysis of inter- and intraspecific evolutionary genomic data. Furthermore, we anticipate that the patterns of structural evolution identified here will serve as a baseline for future comparative studies to identify the factors that influence the evolution of genome structure across taxa.

Chromosome-Level Genome Assembly of the Viviparous Eelpout Zoarces viviparus

The viviparous eelpout Zoarces viviparus is a common fish across the North Atlantic and has successfully colonized habitats across environmental gradients. Due to its wide distribution and predictable phenotypic responses to pollution, Z. viviparus is used as an ideal marine bioindicator organism and has been routinely sampled over decades by several countries to monitor marine environmental health. Additionally, this species is a promising model to study adaptive processes related to environmental change, specifically global warming. Here, we report the chromosome-level genome assembly of Z. viviparus, which has a size of 663 Mb and consists of 607 scaffolds (N50 = 26 Mb). The 24 largest represent the 24 chromosomes of the haploid Z. viviparus genome, which harbors 98% of the complete Benchmarking Universal Single-Copy Orthologues defined for ray-finned fish, indicating that the assembly is highly contiguous and complete. Comparative analyses between the Z. viviparus assembly and the chromosome-level genomes of two other eelpout species revealed a high synteny, but also an accumulation of repetitive elements in the Z. viviparus genome. Our reference genome will be an important resource enabling future in-depth genomic analyses of the effects of environmental change on this important bioindicator species.

Causes and Consequences of Varying Transposable Element Activity: An Evolutionary Perspective

Transposable elements (TEs) are genomic parasites found in nearly all eukaryotes, including humans. This evolutionary success of TEs is due to their replicative activity, involving insertion into new genomic locations. TE activity varies at multiple levels, from between taxa to within individuals. The rapidly accumulating evidence of the influence of TE activity on human health, as well as the rapid growth of new tools to study it, motivated an evaluation of what we know about TE activity thus far. Here, we discuss why TE activity varies, and the consequences of this variation, from an evolutionary perspective. By studying TE activity in nonhuman organisms in the context of evolutionary theories, we can shed light on the factors that affect TE activity. While the consequences of TE activity are usually deleterious, some have lasting evolutionary impacts by conferring benefits on the host or affecting other evolutionary processes.

Genome sequencing of Coryphaenoides yaquinae reveals convergent and lineage-specific molecular evolution in deep-sea adaptation

Abyssal (3501-6500 m) and hadal (>6500 m) fauna evolve under harsh abiotic stresses, characterized by high hydrostatic pressure, darkness and food shortage, providing unique opportunities to investigate mechanisms underlying environmental adaptation. Genomes of several hadal species have recently been reported. However, the genetic adaptation of deep sea species across a broad spectrum of ocean depths has yet to be thoroughly investigated, due to the challenges imposed by collecting the deep sea species. To elucidate the correlation between genetic innovation and vertical distribution, we generated a chromosome-level genome assembly of the macrourids Coryphaenoides yaquinae, which is widely distributed in the abyssal/hadal zone ranging from 3655 to 7259 m in depth. Genomic comparisons among shallow, abyssal and hadal-living species identified idiosyncratic and convergent genetic alterations underlying the extraordinary adaptations of deep-sea species including light perception, circadian regulation, hydrostatic pressure and hunger tolerance. The deep-sea fishes (Coryphaenoides Sp. and Pseudoliparis swirei) venturing into various ocean depths independently have undergone convergent amino acid substitutions in multiple proteins such as rhodopsin 1, pancreatic and duodenal homeobox 1 and melanocortin 4 receptor which are known or verified in zebrafish to be related with vision adaptation and energy expenditure. Convergent evolution events were also identified in heat shock protein 90 beta family member 1 and valosin-containing protein genes known to be related to hydrostatic pressure adaptation specifically in fishes found around the hadal range. The uncovering of the molecular convergence among the deep-sea species shed new light on the common genetic innovations required for deep-sea adaptation by the fishes.

Assessing genotoxic effects of plastic leachates in Drosophila melanogaster

Plastic polymers were largely added with chemical substances to be utilized in the items and product manufacturing. The leachability of these substances is a matter of concern given the wide amount of plastic waste, particularly in terrestrial environments, where soil represents a sink for these novel contaminants and a possible pathway of human health risk. In this study, we integrated genetic, molecular, and behavioral approaches to comparatively evaluate toxicological effects of plastic leachates, virgin and oxodegradable polypropylene (PP) and polyethylene (PE), in Drosophila melanogaster, a novel in vivo model organism for environmental monitoring studies and (eco)toxicological research. The results of this study revealed that while conventional toxicological endpoints such as developmental times and longevity remain largely unaffected, exposure to plastic leachates induces chromosomal abnormalities and transposable element (TE) activation in neural tissues. The combined effects of DNA damage and TE mobilization contribute to genome instability and increase the likelihood of LOH events, thus potentiating tumor growth and metastatic behavior ofRasV12 clones. Collectively, these findings indicate that plastic leachates exert genotoxic effects in Drosophila thus highlighting potential risks associated with leachate-related plastic pollution and their implications for ecosystems and human health.

Extreme overall mushroom genome expansion in Mycena s.s. irrespective of plant hosts or substrate specializations

Mycena s.s. is a ubiquitous mushroom genus whose members degrade multiple dead plant substrates and opportunistically invade living plant roots. Having sequenced the nuclear genomes of 24 Mycena species, we find them to defy the expected patterns for fungi based on both their traditionally perceived saprotrophic ecology and substrate specializations. Mycena displayed massive genome expansions overall affecting all gene families, driven by novel gene family emergence, gene duplications, enlarged secretomes encoding polysaccharide degradation enzymes, transposable element (TE) proliferation, and horizontal gene transfers. Mainly due to TE proliferation, Arctic Mycena species display genomes of up to 502 Mbp (2-8× the temperate Mycena), the largest among mushroom-forming Agaricomycetes, indicating a possible evolutionary convergence to genomic expansions sometimes seen in Arctic plants. Overall, Mycena show highly unusual, varied mosaic-like genomic structures adaptable to multiple lifestyles, providing genomic illustration for the growing realization that fungal niche adaptations can be far more fluid than traditionally believed.

Genetic response of a perennial grass to warm and wet environments interacts and is associated with trait means as well as plasticity

The potential for rapid evolution is an important mechanism allowing species to adapt to changing climatic conditions. Although such potential has been largely studied in various short-lived organisms, to what extent we can observe similar patterns in long-lived plant species, which often dominate natural systems, is largely unexplored. We explored the potential for rapid evolution in Festuca rubra, a long-lived grass with extensive clonal growth dominating in alpine grasslands. We used a field sowing experiment simulating expected climate change in our model region. Specifically, we exposed seeds from five independent seed sources to novel climatic conditions by shifting them along a natural climatic grid and explored the genetic profiles of established seedlings after 3 years. Data on genetic profiles of plants selected under different novel conditions indicate that different climate shifts select significantly different pools of genotypes from common seed pools. Increasing soil moisture was more important than increasing temperature or the interaction of the two climatic factors in selecting pressure. This can indicate negative genetic interaction in response to the combined effects or that the effects of different climates are interactive rather than additive. The selected alleles were found in genomic regions, likely affecting the function of specific genes or their expression. Many of these were also linked to morphological traits (mainly to trait plasticity), suggesting these changes may have a consequence on plant performance. Overall, these data indicate that even long-lived plant species may experience strong selection by climate, and their populations thus have the potential to rapidly adapt to these novel conditions.

Genome Assembly of the Dyeing Poison Frog Provides Insights into the Dynamics of Transposable Element and Genome-Size Evolution

Genome size varies greatly across the tree of life and transposable elements are an important contributor to this variation. Among vertebrates, amphibians display the greatest variation in genome size, making them ideal models to explore the causes and consequences of genome size variation. However, high-quality genome assemblies for amphibians have, until recently, been rare. Here, we generate a high-quality genome assembly for the dyeing poison frog, Dendrobates tinctorius. We compare this assembly to publicly available frog genomes and find evidence for both large-scale conserved synteny and widespread rearrangements between frog lineages. Comparing conserved orthologs annotated in these genomes revealed a strong correlation between genome size and gene size. To explore the cause of gene-size variation, we quantified the location of transposable elements relative to gene features and find that the accumulation of transposable elements in introns has played an important role in the evolution of gene size in D. tinctorius, while estimates of insertion times suggest that many insertion events are recent and species-specific. Finally, we carry out population-scale mobile-element sequencing and show that the diversity and abundance of transposable elements in poison frog genomes can complicate genotyping from repetitive element sequence anchors. Our results show that transposable elements have clearly played an important role in the evolution of large genome size in D. tinctorius. Future studies are needed to fully understand the dynamics of transposable element evolution and to optimize primer or bait design for cost-effective population-level genotyping in species with large, repetitive genomes.

FishTEDB 2.0: an update fish transposable element (TE) database with new functions to facilitate TE research

We launched the initial version of FishTEDB in 2018, which aimed to establish an open-source, user-friendly, data-rich transposable element (TE) database. Over the past 5 years, FishTEDB 1.0 has gained approximately 10 000 users, accumulating more than 450 000 interactions. With the unveiling of extensive fish genome data and the increasing emphasis on TE research, FishTEDB needs to extend the richness of data and functions. To achieve the above goals, we introduced 33 new fish species to FishTEDB 2.0, encompassing a wide array of fish belonging to 48 orders. To make the updated database more functional, we added a genome browser to visualize the positional relationship between TEs and genes and the estimated TE insertion time in different species. In conclusion, we released a new version of the fish TE database, FishTEDB 2.0, designed to assist researchers in the future study of TE functions and promote the progress of biological theories related to TEs. Database URL: https://www.fishtedb.com/.

Mapping and functional characterization of structural variation in 1060 pig genomes

BACKGROUND

Structural variations (SVs) have significant impacts on complex phenotypes by rearranging large amounts of DNA sequence.

RESULTS

We present a comprehensive SV catalog based on the whole-genome sequence of 1060 pigs (Sus scrofa) representing 101 breeds, covering 9.6% of the pig genome. This catalog includes 42,487 deletions, 37,913 mobile element insertions, 3308 duplications, 1664 inversions, and 45,184 break ends. Estimates of breed ancestry and hybridization using genotyped SVs align well with those from single nucleotide polymorphisms. Geographically stratified deletions are observed, along with known duplications of the KIT gene, responsible for white coat color in European pigs. Additionally, we identify a recent SINE element insertion in MYO5A transcripts of European pigs, potentially influencing alternative splicing patterns and coat color alterations. Furthermore, a Yorkshire-specific copy number gain within ABCG2 is found, impacting chromatin interactions and gene expression across multiple tissues over a stretch of genomic region of ~200 kb. Preliminary investigations into SV's impact on gene expression and traits using the Pig Genotype-Tissue Expression (PigGTEx) data reveal SV associations with regulatory variants and gene-trait pairs. For instance, a 51-bp deletion is linked to the lead eQTL of the lipid metabolism regulating gene FADS3, whose expression in embryo may affect loin muscle area, as revealed by our transcriptome-wide association studies.

CONCLUSIONS

This SV catalog serves as a valuable resource for studying diversity, evolutionary history, and functional shaping of the pig genome by processes like domestication, trait-based breeding, and adaptive evolution.

Effects of Arboviral Infections on Transposable Element Transcript Levels in Aedes aegypti

Transposable elements are mobile repeated sequences found in all genomes. Transposable elements are controlled by RNA interference pathways in most organisms, and this control involves the PIWI-interacting RNA pathway and the small interfering RNA pathway, which is also known to be the first line of antiviral defense in invertebrates. Using Drosophila, we recently showed that viral infections result in the modulation of transposable element transcript levels through modulation of the small RNA repertoire. The Aedes aegypti mosquito is of particular interest because almost half of its genome is made of transposable elements, and it is described as a major vector of viruses (such as the dengue [DENV], Zika [ZIKV], and chikungunya [CHIKV] arboviruses). Moreover, Aedes mosquitoes are unique among insects in that the PIWI-interacting RNA pathway is also involved in the somatic antiviral response, in addition to the transposable element control and PIWI-interacting RNA pathway genes expanded in the mosquito genome. For these reasons, we studied the impacts of viral infections on transposable element transcript levels in A. aegypti samples. We retrieved public datasets corresponding to RNA-seq data obtained from viral infections by DENV, ZIKV, and CHIKV in various tissues. We found that transposable element transcripts are moderately modulated following viral infection and that the direction of the modulation varies greatly across tissues and viruses. These results highlight the need for an in-depth investigation of the tightly intertwined interactions between transposable elements and viruses.

The evolution of transposable elements in Brachypodium distachyon is governed by purifying selection, while neutral and adaptive processes play a minor role

Understanding how plants adapt to changing environments and the potential contribution of transposable elements (TEs) to this process is a key question in evolutionary genomics. While TEs have recently been put forward as active players in the context of adaptation, few studies have thoroughly investigated their precise role in plant evolution. Here, we used the wild Mediterranean grass Brachypodium distachyon as a model species to identify and quantify the forces acting on TEs during the adaptation of this species to various conditions, across its entire geographic range. Using sequencing data from more than 320 natural B. distachyon accessions and a suite of population genomics approaches, we reveal that putatively adaptive TE polymorphisms are rare in wild B. distachyon populations. After accounting for changes in past TE activity, we show that only a small proportion of TE polymorphisms evolved neutrally (<10%), while the vast majority of them are under moderate purifying selection regardless of their distance to genes. TE polymorphisms should not be ignored when conducting evolutionary studies, as they can be linked to adaptation. However, our study clearly shows that while they have a large potential to cause phenotypic variation in B. distachyon, they are not favored during evolution and adaptation over other types of mutations (such as point mutations) in this species.

Genomes of historical specimens reveal multiple invasions of LTR retrotransposons in Drosophila melanogaster during the 19th century

Transposable element invasions have a profound impact on the evolution of genomes and phenotypes. It is thus an important open question how often such TE invasions occur. To address this question, we utilize the genomes of historical specimens, sampled about 200 y ago. We found that the LTR retrotransposons Blood, Opus, and 412 spread in Drosophila melanogaster in the 19th century. These invasions constitute second waves, as degraded fragments were found for all three TEs. The composition of Opus and 412, but not of Blood, shows a pronounced geographic heterogeneity, likely due to founder effects during the invasions. Finally, we identified species from the Drosophila simulans complex as the likely origin of the TEs. We show that in total, seven TE families invaded D. melanogaster during the last 200y, thereby increasing the genome size by up to 1.2Mbp. We suggest that this high rate of TE invasions was likely triggered by human activity. Based on the analysis of strains and specimens sampled at different times, we provide a detailed timeline of TE invasions, making D. melanogaster the first organism where the invasion history of TEs during the last two centuries could be inferred.

The global distribution of angiosperm genome size is shaped by climate

Angiosperms, which inhabit diverse environments across all continents, exhibit significant variation in genome sizes, making them an excellent model system for examining hypotheses about the global distribution of genome size. These include the previously proposed large genome constraint, mutational hazard, polyploidy-mediated, and climate-mediated hypotheses. We compiled the largest genome size dataset to date, encompassing 16 017 (> 5% of known) angiosperm species, and analyzed genome size distribution using a comprehensive geographic distribution dataset for all angiosperms. We observed that angiosperms with large range sizes generally had small genomes, supporting the large genome constraint hypothesis. Climate was shown to exert a strong influence on genome size distribution along the global latitudinal gradient, while the frequency of polyploidy and the type of growth form had negligible effects. In contrast to the unimodal patterns along the global latitudinal gradient shown by plant size traits and polyploid proportions, the increase in angiosperm genome size from the equator to 40-50°N/S is probably mediated by different (mostly climatic) mechanisms than the decrease in genome sizes observed from 40 to 50°N northward. Our analysis suggests that the global distribution of genome sizes in angiosperms is mainly shaped by climatically mediated purifying selection, genetic drift, relaxed selection, and environmental filtering.

Transcriptional responses to salinity-induced changes in cell wall morphology of the euryhaline diatom Pleurosira laevis

Diatoms are unicellular algae with morphologically diverse silica cell walls, which are called frustules. The mechanism of frustule morphogenesis has attracted attention in biology and nanomaterials engineering. However, the genetic regulation of the morphology remains unclear. We therefore used transcriptome sequencing to search for genes involved in frustule morphology in the centric diatom Pleurosira laevis, which exhibits morphological plasticity between flat and domed valve faces in salinity 2 and 7, respectively. We observed differential expression of transposable elements (TEs) and transporters, likely due to osmotic response. Up-regulation of mechanosensitive ion channels and down-regulation of Ca2+-ATPases in cells with flat valves suggested that cytosolic Ca2+ levels were changed between the morphologies. Calcium signaling could be a mechanism for detecting osmotic pressure changes and triggering morphological shifts. We also observed an up-regulation of ARPC1 and annexin, involved in the regulation of actin filament dynamics known to affect frustule morphology, as well as the up-regulation of genes encoding frustule-related proteins such as BacSETs and frustulin. Taken together, we propose a model in which salinity-induced morphogenetic changes are driven by upstream responses, such as the regulation of cytosolic Ca2+ levels, and downstream responses, such as Ca2+-dependent regulation of actin dynamics and frustule-related proteins. This study highlights the sensitivity of euryhaline diatoms to environmental salinity and the role of active cellular processes in controlling gross valve morphology under different osmotic pressures.

Increasing transposase abundance with ocean depth correlates with a particle-associated lifestyle

Transposases are mobile genetic elements that move within and between genomes, promoting genomic plasticity in microorganisms. In marine microbial communities, the abundance of transposases increases with depth, but the reasons behind this trend remain unclear. Our analysis of metagenomes from the Tara Oceans and Malaspina Expeditions suggests that a particle-associated lifestyle is the main covariate for the high occurrence of transposases in the deep ocean, and this trend holds true for individual genomes as well as in a community-wide sense. We observed a strong and depth-independent correlation between transposase abundance and the presence of biofilm-associated genes, as well as the prevalence of secretory enzymes. This suggests that mobile genetic elements readily propagate among microbial communities within crowded biofilms. Furthermore, we show that particle association positively correlates with larger genome size, which is in turn associated with higher transposase abundance. Cassette sequences associated with transposons are enriched with genes related to defense mechanisms, which are more highly expressed in the deep sea. Thus, while transposons spread at the expense of their microbial hosts, they also introduce novel genes and potentially benefit the hosts in helping to compete for limited resources. Overall, our results suggest a new understanding of deep ocean particles as highways for gene sharing among defensively oriented microbial genomes.IMPORTANCEGenes can move within and between microbial genomes via mobile genetic elements, which include transposases and transposons. In the oceans, there is a puzzling increase in transposase abundance in microbial genomes as depth increases. To gain insight into this trend, we conducted an extensive analysis of marine microbial metagenomes and metatranscriptomes. We found a significant correlation between transposase abundance and a particle-associated lifestyle among marine microbes at both the metagenome and genome-resolved levels. We also observed a link between transposase abundance and genes related to defense mechanisms. These results suggest that as microbes become densely packed into crowded particles, mobile genes are more likely to spread and carry genetic material that provides a competitive advantage in crowded habitats. This may enable deep sea microbes to effectively compete in such environments.

Plant adaptation to ionizing radiation: Mechanisms and patterns

Adaptation to environmental stressors is an essential property of plants that allows them, despite an immobile lifestyle, to survive in a changeable environment. The chain of successive events culminating in the final radiobiological reaction begins with the absorption of energy of ionizing radiation in the cell. Starting from stochastic acts of molecular injury formation, radiation damage gradually acquires deterministic features, which are expressed in a limited number of phenomena that complete plant radiation damage. As plants undergo specialization, the differences between plants and animals become more pronounced, leading to distinct responses to radiation. Chronic radiation exposure may activate biological mechanisms resulting in increased radioresistance of the population. The higher the level of radiation exposure and the sensitivity of plants to radiation, the more intensive the selection. Depending on the circumstances, enhanced radioresistance of a population can be achieved in different ways or has not evolved at all. High dose rates of chronic irradiation leаd to selection for the efficiency of repair systems, while low dose rates activate epigenetic mechanisms that lead to the maintenance of oxidative balance, additional synthesis of chaperones, and control of TEs transposition. Due to huge differences in the radiosensitivity of organisms that make up the ecosystem, irradiation can result in disruption of connections between components of ecosystems which may lead to consequences that can differ drastically from those expected at the organismal and population levels. Therefore, the use of ecological knowledge is essential for understanding the responses of populations and ecosystems to radiation exposure.

Barley AGO4 proteins show overlapping functionality with distinct small RNA-binding properties in heterologous complementation

KEY MESSAGE

Barley AGO4 proteins complement expressional changes of epigenetically regulated genes in Arabidopsis ago4-3 mutant and show a distinct affinity for the 5' terminal nucleotide of small RNAs, demonstrating functional conservation and divergence. The function of Argonaute 4 (AGO4) in Arabidopsis thaliana has been extensively characterized; however, its role in monocots, which have large genomes abundantly supplemented with transposable elements (TEs), remains elusive. The study of barley AGO4 proteins can provide insights into the conserved aspects of RNA-directed DNA methylation (RdDM) and could also have further applications in the field of epigenetics or crop improvement. Bioinformatic analysis of RNA sequencing data identified two active AGO4 genes in barley, HvAGO4a and HvAGO4b. These genes function similar to AtAGO4 in an Arabidopsis heterologous complementation system, primarily binding to 24-nucleotide long small RNAs (sRNAs) and triggering methylation at specific target loci. Like AtAGO4, HvAGO4B exhibits a preference for binding sRNAs with 5' adenine residue, while also accepting 5' guanine, uracil, and cytosine residues. In contrast, HvAGO4A selectively binds only sRNAs with a 5' adenine residue. The diverse binding capacity of barley AGO4 proteins is reflected in TE-derived sRNAs and in their varying abundance. Both barley AGO4 proteins effectively restore the levels of extrachromosomal DNA and transcript abundancy of the heat-activated ONSEN retrotransposon to those observed in wild-type Arabidopsis plants. Our study provides insight into the distinct binding specificities and involvement in TE regulation of barley AGO4 proteins in Arabidopsis by heterologous complementation.

Dissection of figured wood trait in curly birch (Betula pendula Roth var. carelica (Mercklin) Hämet-Ahti) using high-throughput genotyping

Curly (Karelian) birch is a special variety of Betula pendula Roth distributed in the northwestern part of Europe. Karelian birch is well-known for its valuable figured curly wood also known as "wooden marble". The genetic basis underlying curly wood formation has been debated since last century, however, there was no data about loci responsible for the curly wood trait. In the present study, we analyzed two full-sibs populations derived from experimental crosses of curly birches and segregating for the trait. RADseq genotyping was applied to reveal how many loci are involved in 'curliness' formation and to search for genetic variants associated with this trait. One single interval on chromosome 10 was detected containing possible candidate genes. InDel marker BpCW1 was suggested for the first time for marker-assisted selection of trees with curly wood at their earliest stages of development.

Transposable Element Insertions Are Associated with Batesian Mimicry in the Pantropical Butterfly Hypolimnas misippus

Hypolimnas misippus is a Batesian mimic of the toxic African Queen butterfly (Danaus chrysippus). Female H. misippus butterflies use two major wing patterning loci (M and A) to imitate three color morphs of D. chrysippus found in different regions of Africa. In this study, we examine the evolution of the M locus and identify it as an example of adaptive atavism. This phenomenon involves a morphological reversion to an ancestral character that results in an adaptive phenotype. We show that H. misippus has re-evolved an ancestral wing pattern present in other Hypolimnas species, repurposing it for Batesian mimicry of a D. chrysippus morph. Using haplotagging, a linked-read sequencing technology, and our new analytical tool, Wrath, we discover two large transposable element insertions located at the M locus and establish that these insertions are present in the dominant allele responsible for producing mimetic phenotype. By conducting a comparative analysis involving additional Hypolimnas species, we demonstrate that the dominant allele is derived. This suggests that, in the derived allele, the transposable elements disrupt a cis-regulatory element, leading to the reversion to an ancestral phenotype that is then utilized for Batesian mimicry of a distinct model, a different morph of D. chrysippus. Our findings present a compelling instance of convergent evolution and adaptive atavism, in which the same pattern element has independently evolved multiple times in Hypolimnas butterflies, repeatedly playing a role in Batesian mimicry of diverse model species.

Insights into the genome of the 'Loco' Concholepas concholepas (Gastropoda: Muricidae) from low-coverage short-read sequencing: genome size, ploidy, transposable elements, nuclear RNA gene operon, mitochondrial genome, and phylogenetic placement in the family Muricidae

BACKGROUND

The Peruvian 'chanque' or Chilean 'loco' Concholepas concholepas is an economically, ecologically, and culturally important muricid gastropod heavily exploited by artisanal fisheries in the temperate southeastern Pacific Ocean. In this study, we have profited from a set of bioinformatics tools to recover important biological information of C. concholepas from low-coverage short-read NGS datasets. Specifically, we calculated the size of the nuclear genome, ploidy, and estimated transposable elements content using an in silico k-mer approach, we discovered, annotated, and quantified those transposable elements, we assembled and annotated the 45S rDNA RNA operon and mitochondrial genome, and we confirmed the phylogenetic position of C. concholepas within the muricid subfamily Rapaninae based on translated protein coding genes.

RESULTS

Using a k-mer approach, the haploid genome size estimated for the predicted diploid genome of C. concholepas varied between 1.83 Gbp (with kmer = 24) and 2.32 Gbp (with kmer = 36). Between half and two thirds of the nuclear genome of C. concholepas was composed of transposable elements. The most common transposable elements were classified as Long Interspersed Nuclear Elements and Short Interspersed Nuclear Elements, which were more abundant than DNA transposons, simple repeats, and Long Terminal Repeats. Less abundant repeat elements included Helitron mobile elements, 45S rRNA DNA, and Satellite DNA, among a few others.The 45S rRNA DNA operon of C. concholepas that encodes for the ssrRNA, 5.8S rRNA, and lsrRNA genes was assembled into a single contig 8,090 bp long. The assembled mitochondrial genome of C. concholepas is 15,449 bp long and encodes 13 protein coding genes, two ribosomal genes, and 22 transfer RNAs.

CONCLUSION

The information gained by this study will inform the assembly of a high quality nuclear genome for C. concholepas and will support bioprospecting and biomonitoring using environmental DNA to advance development of conservation and management plans in this overexploited marine snail.

Mechanisms of Rapid Karyotype Evolution in Mammals

Chromosome reshuffling events are often a foundational mechanism by which speciation can occur, giving rise to highly derivative karyotypes even amongst closely related species. Yet, the features that distinguish lineages prone to such rapid chromosome evolution from those that maintain stable karyotypes across evolutionary time are still to be defined. In this review, we summarize lineages prone to rapid karyotypic evolution in the context of Simpson's rates of evolution-tachytelic, horotelic, and bradytelic-and outline the mechanisms proposed to contribute to chromosome rearrangements, their fixation, and their potential impact on speciation events. Furthermore, we discuss relevant genomic features that underpin chromosome variation, including patterns of fusions/fissions, centromere positioning, and epigenetic marks such as DNA methylation. Finally, in the era of telomere-to-telomere genomics, we discuss the value of gapless genome resources to the future of research focused on the plasticity of highly rearranged karyotypes.

The impact of paramutations on the invasion dynamics of transposable elements

According to the prevailing view, the trap model, the activity of invading transposable elements (TEs) is greatly reduced when a TE copy jumps into a piRNA cluster, which triggers the emergence of piRNAs that silence the TE. One crucial component in the host defence are paramutations. Mediated by maternally deposited piRNAs, paramutations convert TE insertions into piRNA producing loci, thereby transforming selfish TEs into agents of the host defence. Despite this significant effect, the impact of paramutations on the dynamics of TE invasions remains unknown. To address this issue, we performed extensive forward simulations of TE invasions with piRNA clusters and paramutations. We found that paramutations significantly affect TE dynamics, by accelerating the silencing of TE invasions, reducing the number of insertions accumulating during the invasions and mitigating the fitness cost of TEs. We also demonstrate that piRNA production induced by paramutations, an epigenetically inherited trait, may be positively selected. Finally, we show that paramutations may account for three important open problems with the trap model. Firstly, paramutated TE insertions may compensate for the insufficient number of insertions in piRNA clusters observed in previous studies. Secondly, paramutations may explain the discrepancy between the observed and the expected abundance of different TE families in Drosophila melanogaster. Thirdly, piRNA clusters may be crucial to trigger the host defence, but paramutations render the clusters dispensable once the defence has been established. This could account for the lack of TE activation when three major piRNA clusters were deleted in a previous study.

Quantifying the role of genome size and repeat content in adaptive variation and the architecture of flowering time in Amaranthus tuberculatus

Genome size variation, largely driven by repeat content, is poorly understood within and among populations, limiting our understanding of its significance for adaptation. Here we characterize intraspecific variation in genome size and repeat content across 186 individuals of Amaranthus tuberculatus, a ubiquitous native weed that shows flowering time adaptation to climate across its range and in response to agriculture. Sequence-based genome size estimates vary by up to 20% across individuals, consistent with the considerable variability in the abundance of transposable elements, unknown repeats, and rDNAs across individuals. The additive effect of this variation has important phenotypic consequences-individuals with more repeats, and thus larger genomes, show slower flowering times and growth rates. However, compared to newly-characterized gene copy number and polygenic nucleotide changes underlying variation in flowering time, we show that genome size is a marginal contributor. Differences in flowering time are reflected by genome size variation across sexes and marginally, habitats, while polygenic variation and a gene copy number variant within the ATP synthesis pathway show consistently stronger environmental clines than genome size. Repeat content nonetheless shows non-neutral distributions across the genome, and across latitudinal and environmental gradients, demonstrating the numerous governing processes that in turn influence quantitative genetic variation for phenotypes key to plant adaptation.

Stress Induced Activation of LTR Retrotransposons in the Drosophila melanogaster Genome

Background: Retrotransposons with long terminal repeats (LTR retrotransposons) are widespread in all groups of eukaryotes and are often both the cause of new mutations and the source of new sequences. Apart from their high activity in generative and differentiation-stage tissues, LTR retrotransposons also become more active in response to different stressors. The precise causes of LTR retrotransposons' activation in response to stress, however, have not yet been thoroughly investigated. Methods: We used RT-PCR to investigate the transcriptional profile of LTR retrotransposons and piRNA clusters in response to oxidative and chronic heat stresses. We used Oxford Nanopore sequencing to investigate the genomic environment of new insertions of the retrotransposons. We used bioinformatics methods to find the stress-induced transcription factor binding sites in LTR retrotransposons. Results: We studied the transposition activity and transcription level of LTR retrotransposons in response to oxidative and chronic heat stress and assessed the contribution of various factors that can affect the increase in their expression under stress conditions: the state of the piRNA-interference system, the influence of the genomic environment on individual copies, and the presence of the stress-induced transcription factor binding sites in retrotransposon sequences. Conclusions: The main reason for the activation of LTR retrotransposons under stress conditions is the presence of transcription factor binding sites in their regulatory sequences, which are triggered in response to stress and are necessary for tissue regeneration processes. Stress-induced transposable element activation can function as a trigger mechanism, triggering multiple signal pathways and resulting in a polyvariant cell response.

Draft genome sequencing of Tilletia caries inciting common bunt of wheat provides pathogenicity-related genes

Common bunt of wheat caused by Tilletia caries is an important disease worldwide. The T. caries TC1_MSG genome was sequenced using the Illumina HiSeq 2500 and Nanopore ONT platforms. The Nanopore library was prepared using the ligation sequencing kit SQK-LSK110 to generate approximately 24 GB for sequencing. The assembly size of 38.18 Mb was generated with a GC content of 56.10%. The whole genome shotgun project was deposited at DDBJ/ENA/GenBank under the accession number JALUTQ000000000. Forty-six contigs were obtained with N50 of 1,798,756 bp. In total, 10,698 genes were predicted in the assembled genome. Out of 10,698 genes, 10,255 genes were predicted significantly in the genome. The repeat sequences made up approximately 1.57% of the genome. Molecular function, cellular components, and biological processes for predicted genes were mapped into the genome. In addition, repeat elements in the genome were assessed. In all, 0.89% of retroelements were observed, followed by long terminal repeat elements (0.86%) in the genome. In simple sequence repeat (SSR) analysis, 8,582 SSRs were found in the genome assembly. The trinucleotide SSR type (3,703) was the most abundant. Few putative secretory signal peptides and pathogenicity-related genes were predicted. The genomic information of T. caries will be valuable in understanding the pathogenesis mechanism as well as developing new methods for the management of the common bunt disease of wheat.

The First Genome Survey of the Snail Provanna glabra Inhabiting Deep-Sea Hydrothermal Vents

The snail P. glabra is an endemic species in deep-sea chemosynthetic ecosystems of the Northwest Pacific Ocean. To obtain more genetic information on this species and provide the basis for subsequent whole-genome map construction, a genome survey was performed on this snail from the hydrothermal vent of Okinawa Trough. The genomic size of P. glabra was estimated to be 1.44 Gb, with a heterozygosity of 1.91% and a repeated sequence content of 69.80%. Based on the sequencing data, a draft genome of 1.32 Gb was assembled. Transposal elements (TEs) accounted for 40.17% of the entire genome, with DNA transposons taking the highest proportion. It was found that most TEs were inserted in the genome recently. In the simple sequence repeats, the dinucleotide motif was the most enriched microsatellite type, accounting for 53% of microsatellites. A complete mitochondrial genome of P. glabra with a total length of 16,268 bp was assembled from the sequencing data. After comparison with the published mitochondrial genome of Provanna sp. from a methane seep, 331 potential single nucleotide polymorphism (SNP) sites were identified in protein-coding genes (PCGs). Except for the cox1 gene, nad2, nad4, nad5, and cob genes are expected to be candidate markers for population genetic and phylogenetic studies of P. glabra and other deep-sea snails. Compared with shallow-water species, three mitochondrial genes of deep-sea gastropods exhibited a higher evolutionary rate, indicating strong selection operating on mitochondria of deep-sea species. This study provides insights into the genome characteristics of P. glabra and supplies genomic resources for further studies on the adaptive evolution of the snail in extreme deep-sea chemosynthetic environments.

The composition of piRNA clusters in Drosophila melanogaster deviates from expectations under the trap model

BACKGROUND

It is widely assumed that the invasion of a transposable element (TE) in mammals and invertebrates is stopped when a copy of the TE jumps into a piRNA cluster (i.e., the trap model). However, recent works, which for example showed that deletion of three major piRNA clusters has no effect on TE activity, cast doubt on the trap model.

RESULTS

Here, we test the trap model from a population genetics perspective. Our simulations show that the composition of regions that act as transposon traps (i.e., potentially piRNA clusters) ought to deviate from regions that have no effect on TE activity. We investigated TEs in five Drosophila melanogaster strains using three complementary approaches to test whether the composition of piRNA clusters matches these expectations. We found that the abundance of TE families inside and outside of piRNA clusters is highly correlated, although this is not expected under the trap model. Furthermore, the distribution of the number of TE insertions in piRNA clusters is also much broader than expected.

CONCLUSIONS

We found that the observed composition of piRNA clusters is not in agreement with expectations under the simple trap model. Dispersed piRNA producing TE insertions and temporal as well as spatial heterogeneity of piRNA clusters may account for these deviations.

Genomic Plasticity of Acid-Tolerant Phenotypic Evolution in Acetobacter pasteurianus

Acetic acid bacteria have a remarkable capacity to cope with elevated concentrations of cytotoxic acetic acid in their fermentation environment. In particular, the high-level acetate tolerance of Acetobacter pasteurianus that occurs in vinegar industrial settings must be constantly selected for. However, the improved acetic acid tolerance is rapidly lost without a selection pressure. To understand genetic and molecular biology of this acquired acetic acid tolerance in A. pasteurianus, we evolved three strains A. pasteurianus CICIM B7003, CICIM B7003-02, and ATCC 33,445 over 960 generations (4 months) in two initial acetic acids of 20 g·L-1 and 30 g·L-1, respectively. An acetic acid-adapted strain M20 with significantly improved specific growth rate of 0.159 h-1 and acid productivity of 1.61 g·L-1·h-1 was obtained. Comparative genome analysis of six evolved strains revealed that the genetic variations of adaptation were mainly focused on lactate metabolism, membrane proteins, transcriptional regulators, transposases, replication, and repair system. Among of these, lactate dehydrogenase, acetolactate synthase, glycosyltransferase, ABC transporter ATP-binding protein, two-component regulatory systems, the type II toxin-antitoxin system (RelE/RelB/StbE), exodeoxyribonuclease III, type I restriction endonuclease, tRNA-uridine 2-sulfurtransferase, and transposase might collaboratively contribute to the improved acetic acid tolerance in A. pasteurianus strains. The balance between repair factors and transposition variations might be the basis for genomic plasticity of A. pasteurianus strains, allowing the survival of populations and their offspring in acetic acid stress fluctuations. These observations provide important insights into the nature of acquired acetic acid tolerance phenotype and lay a foundation for future genetic manipulation of these strains.

From genome size to trait evolution during angiosperm radiation

Angiosperm diversity arises from trait flexibility and repeated evolutionary radiations, but the role of genomic characters in these radiations remains unclear. In this opinion article, we discuss how genome size can influence angiosperm diversification via its intricate link with cell size, tissue packing, and physiological processes which, in turn, influence the macroevolution of functional traits. We propose that integrating genome size, functional traits, and phylogenetic data across a wide range of lineages allows us to test whether genome size decrease consistently leads to increased trait flexibility, while genome size increase constrains trait evolution. Combining theories from molecular biology, functional ecology and macroevolution, we provide a framework to better understand the role of genome size in trait evolution, evolutionary radiations, and the global distribution of angiosperms.

That's no moon, it's a Starship: Giant transposons driving fungal horizontal gene transfer

To date, most reports of horizontal gene transfer (HGT) in fungi rely on genome sequence data and are therefore an indirect measure of HGT after the event has occurred. However, a novel group of class II-like transposons known as Starships may soon alter this status quo. Starships are giant transposable elements that carry dozens of genes, some of which are host-beneficial, and are linked to many recent HGT events in the fungal kingdom. These transposons remain active and mobile in many fungal genomes and their transposition has recently been shown to be driven by a conserved tyrosine-recombinase called 'Captain'. This perspective explores some of the remaining unanswered questions about how these Starship transposons move, both within a genome and between different species. We seek to outline several experimental approaches that can be used to identify the genes essential for Starship-mediated HGT and draw links to other recently discovered giant transposons outside of the fungal kingdom.

Impact of the antibiotic doxycycline on the D. magna reproduction, associated microbiome and antibiotic resistance genes in treated wastewater conditions

Wastewater Treatment Plant (WWTP) effluents are important sources of antibiotics, antibiotic resistance genes (ARGs) and resistant bacteria that threaten aquatic biota and human heath. Antibiotic effects on host-associated microbiomes, spread of ARGs and the consequences for host health are still poorly described. This study investigated changes of the Daphnia magna associated microbiome exposed to the recalcitrant antibiotic doxycycline under artificial reconstituted lab water media (lab water) and treated wastewater media. D. magna individual juveniles were exposed for 10 days to treated wastewater with and without doxycycline, and similarly in lab water. We analysed 16 S rRNA gene sequences to assess changes in community structure, monitored Daphnia offspring production and quantified ARGs abundances by qPCR from both Daphnia and water (before and after the exposure). Results showed that doxycycline and media (lab water or wastewater) had a significant effect modulating Daphnia-associated microbiome composition and one of the most discriminant taxa was Enterococcus spp. Moreover, in lab water, doxycycline reduced the presence of Limnohabitans sp., which are dominant bacteria of the D. magna-associated microbiome and impaired Daphnia reproduction. Contrarily, treated wastewater increased diversity and richness of Daphnia-associated microbiome and promoted fecundity. In addition, the detected ARG genes in both lab water and treated wastewater medium included the qnrS1, sul1, and blaTEM, and the integron-related intI1 gene. The treated wastewater contained about 10 times more ARGs than lab water alone. Furthermore, there was an increase of sul1 in Daphnia cultured in treated wastewater compared to lab water. In addition, there were signs of a higher biodegradation of doxycycline by microbiomes of treated wastewater in comparison to lab water. Thus, results suggest that Daphnia-associated microbiomes are influenced by their environment, and that bacterial communities present in treated wastewater are better suited to cope with the effects of antibiotics.

Effects of soil solarization combined with manure-amended on soil ARGs and microbial communities during summer fallow

Soil solarization (SS) is a technique for managing pathogens and weeds, which involves covering with transparent plastic to increase soil temperature during summer fallow (SF). However, SS also alters the diversity of bacterial communities. Therefore, during SF, various organic modifiers are used in combination with SS to improve its efficacy. Organic amendments may contain antibiotic resistance genes (ARGs). Greenhouse vegetable production (GVP) soils are vital to ensure food security and ecological balance. However, comprehensive study on the effects of SS combined with different types of manure on ARGs in GVP soils during SF remains unclear. Therefore, this study employed high-throughput qPCR to explore the effects of different organic amendments combined with SS on the abundance changes of ARGs and mobile genetic elements (MGEs) in GVP soils during SF. The abundance and diversity of ARGs and MGEs in GVP soils with different manure fertilization and SS decreased during SF. Horizontal gene transfer via MGEs (especially integrases 45.80%) induced by changes in environmental factors (NO3--N 14.7% and NH4+-N) was the main factor responsible for the changes in ARGs. Proteobacteria (14.3%) and Firmicutes were the main potential hosts of ARGs. Network analysis suggested that Ornithinimicrobium, Idiomarina and Corynebacterium had positive correlations with aminoglycosides, MLSB, and tetracycline resistance genes. These results provide new insights to understand the fate of ARGs in the GVP soils by manure-amended combined with SS during SF, which may help to reduce the spread of ARGs.

Transposons contribute to the functional diversification of the head, gut, and ovary transcriptomes across Drosophila natural strains

Transcriptomes are dynamic, with cells, tissues, and body parts expressing particular sets of transcripts. Transposable elements (TEs) are a known source of transcriptome diversity; however, studies often focus on a particular type of chimeric transcript, analyze single body parts or cell types, or are based on incomplete TE annotations from a single reference genome. In this work, we have implemented a method based on de novo transcriptome assembly that minimizes the potential sources of errors while identifying a comprehensive set of gene-TE chimeras. We applied this method to the head, gut, and ovary dissected from five Drosophila melanogaster natural strains, with individual reference genomes available. We found that ∼19% of body part-specific transcripts are gene-TE chimeras. Overall, chimeric transcripts contribute a mean of 43% to the total gene expression, and they provide protein domains for DNA binding, catalytic activity, and DNA polymerase activity. Our comprehensive data set is a rich resource for follow-up analysis. Moreover, because TEs are present in virtually all species sequenced to date, their role in spatially restricted transcript expression is likely not exclusive to the species analyzed in this work.

Exploring Prokaryotic Communities in the Guts and Mucus of Nudibranchs, and Their Similarity to Sediment and Seawater Microbiomes

In the present study, we compared mucus and gut-associated prokaryotic communities from seven nudibranch species with sediment and seawater from Thai coral reefs using high-throughput 16S rRNA gene sequencing. The nudibranch species were identified as Doriprismatica atromarginata (family Chromodorididae), Jorunna funebris (family Discodorididae), Phyllidiella nigra, Phyllidiella pustulosa, Phyllidia carlsonhoffi, Phyllidia elegans, and Phyllidia picta (all family Phyllidiidae). The most abundant bacterial phyla in the dataset were Proteobacteria, Tenericutes, Chloroflexi, Thaumarchaeota, and Cyanobacteria. Mucus and gut-associated communities differed from one another and from sediment and seawater communities. Host phylogeny was, furthermore, a significant predictor of differences in mucus and gut-associated prokaryotic community composition. With respect to higher taxon abundance, the order Rhizobiales (Proteobacteria) was more abundant in Phyllidia species (mucus and gut), whereas the order Mycoplasmatales (Tenericutes) was more abundant in D. atromarginata and J. funebris. Mucus samples were, furthermore, associated with greater abundances of certain phyla including Chloroflexi, Poribacteria, and Gemmatimonadetes, taxa considered to be indicators for high microbial abundance (HMA) sponge species. Overall, our results indicated that nudibranch microbiomes consisted of a number of abundant prokaryotic members with high sequence similarities to organisms previously detected in sponges.

maT and mosquito transposons in cnidarians: evolutionary history and intraspecific differences

Transposable elements exert a significant effect on the size and structure of eukaryotic genomes. Tc1/mariner superfamily elements represent the widely distributed and highly variable group of DNA transposons. Tc1/mariner elements include TLE/DD34-38E, MLE/DD34D, maT/DD37D, Visitor/DD41D, Guest/DD39D, mosquito/DD37E, and L18/DD37E families, all of which are well or less scarcely studied. However, more detailed research into the patterns of prevalence and diversity of Tc1/mariner transposons enables one to better understand the coevolution of the TEs and the eukaryotic genomes. We performed a detailed analysis of the maT/DD37D family in Cnidaria. The study of 77 genomic assemblies demonstrated that maT transposons are found in a limited number of cnidarian species belonging to classes Cubozoa (1 species), Hydrozoa (3 species) и Scyphozoa (5 species) only. The identified TEs were classified into 5 clades, with the representatives from Pelagiidae (class Scyphozoa) forming a separate clade of maT transposons, which has never been described previously. The potentially functional copies of maT transposons were identified in the hydrae. The phylogenetic analysis and the studies of distribution among the taxons and the evolutionary dynamics of the elements suggest that maT transposons of the cnidarians are the descendants of several independent invasion events occurring at different periods of time. We also established that the TEs of mosquito/DD37E family are found in Hydridae (class Hydrozoa) only. A comparison of maT and mosquito prevalence in two genomic assemblies of Hydra viridissima revealed obvious differences, thus demonstrating that each individual organism might carry a unique mobilome pattern. The results of the presented research make us better understand the diversity and evolution of Tc1/mariner transposons and their effect on the eukaryotic genomes.

Long-read direct RNA sequencing reveals epigenetic regulation of chimeric gene-transposon transcripts in Arabidopsis thaliana

Transposable elements (TEs) are accumulated in both intergenic and intragenic regions in plant genomes. Intragenic TEs often act as regulatory elements of associated genes and are also co-transcribed with genes, generating chimeric TE-gene transcripts. Despite the potential impact on mRNA regulation and gene function, the prevalence and transcriptional regulation of TE-gene transcripts are poorly understood. By long-read direct RNA sequencing and a dedicated bioinformatics pipeline, ParasiTE, we investigated the transcription and RNA processing of TE-gene transcripts in Arabidopsis thaliana. We identified a global production of TE-gene transcripts in thousands of A. thaliana gene loci, with TE sequences often being associated with alternative transcription start sites or transcription termination sites. The epigenetic state of intragenic TEs affects RNAPII elongation and usage of alternative poly(A) signals within TE sequences, regulating alternative TE-gene isoform production. Co-transcription and inclusion of TE-derived sequences into gene transcripts impact regulation of RNA stability and environmental responses of some loci. Our study provides insights into TE-gene interactions that contributes to mRNA regulation, transcriptome diversity, and environmental responses in plants.

Woronichinia naegeliana: A common nontoxigenic component of temperate freshwater cyanobacterial blooms with 30% of its genome in transposons

Monitoring in the U.S. state of Washington across the period 2007-2019 showed that Woronichinia has been present in many lakes state-wide. This cyanobacterium was commonly dominant or sub-dominant in cyanobacterial blooms in the wet temperate region west of the Cascade Mountains. In these lakes, Woronichinia often co-existed with Microcystis, Dolichospermum and Aphanizomenon flos-aquae and the cyanotoxin microcystin has often been present in those blooms, although it has not been known whether Woronichinia is a toxin producer. We report the first complete genome of Woronichinia naegeliana WA131, assembled from the metagenome of a sample collected from Wiser Lake, Washington, in 2018. The genome contains no genes for cyanotoxin biosynthesis or taste-and-odor compounds, but there are biosynthetic gene clusters for other bioactive peptides, including anabaenopeptins, cyanopeptolins, microginins and ribosomally produced, post-translationally modified peptides. Genes for photosynthesis, nutrient acquisition, vitamin synthesis and buoyancy that are typical of bloom-forming cyanobacteria are present, although nitrate and nitrite reductase genes are conspicuously absent. However, the 7.9 Mbp genome is 3-4 Mbp larger than those of the above-mentioned frequently co-existing cyanobacteria. The increased genome size is largely due to an extraordinary number of insertion sequence elements (transposons), which account for 30.3% of the genome and many of which are present in multiple copies. The genome contains a relatively large number of pseudogenes, 97% of which are transposase genes. W. naegeliana WA131 thus seems to be able to limit the potentially deleterious effects of high rates of recombination and transposition to the mobilome fraction of its genome.