Large Genome Research Articles

Diatom abundance in the polar oceans is predicted by genome size.

A principal goal in ecology is to identify the determinants of species abundances in nature. Body size has emerged as a fundamental and repeatable predictor of abundance, with smaller organisms occurring in greater numbers than larger ones. A biogeographic component, known as Bergmann's rule, describes the preponderance, across taxonomic groups, of larger-bodied organisms in colder areas. Although undeniably important, the extent to which body size is the key trait underlying these patterns is unclear. We explored these questions in diatoms, unicellular algae of global importance for their roles in carbon fixation and energy flow through marine food webs. Using a phylogenomic dataset from a single lineage with worldwide distribution, we found that body size (cell volume) was strongly correlated with genome size, which varied by 50-fold across species and was driven by differences in the amount of repetitive DNA. However, directional models identified temperature and genome size, not cell size, as having the greatest influence on maximum population growth rate. A global metabarcoding dataset further identified genome size as a strong predictor of species abundance in the ocean, but only in colder regions at high and low latitudes where diatoms with large genomes dominated, a pattern consistent with Bergmann's rule. Although species abundances are shaped by myriad interacting abiotic and biotic factors, genome size alone was a remarkably strong predictor of abundance. Taken together, these results highlight the cascading cellular and ecological consequences of macroevolutionary changes in an emergent trait, genome size, one of the most fundamental and irreducible properties of an organism.

Genome size is positively correlated with extinction risk in herbaceous angiosperms.

Angiosperms with large genomes experience nuclear-, cellular-, and organism-level constraints that may limit their phenotypic plasticity and ecological niche, which could increase their risk of extinction. Therefore, we test the hypotheses that large-genomed species are more likely to be threatened with extinction than those with small genomes, and that the effect of genome size varies across three selected covariates: life form, endemism, and climatic zone. We collated genome size and extinction risk information for a representative sample of angiosperms comprising 3250 species, which we analyzed alongside life form, endemism, and climatic zone variables using a phylogenetic framework. Genome size is positively correlated with extinction risk, a pattern driven by a signal in herbaceous but not woody species, regardless of climate and endemism. The influence of genome size is stronger in endemic herbaceous species, but is relatively homogenous across different climates. Beyond its indirect link via endemism and climate, genome size is associated with extinction risk directly and significantly. Genome size may serve as a proxy for difficult-to-measure parameters associated with resilience and vulnerability in herbaceous angiosperms. Therefore, it merits further exploration as a useful biological attribute for understanding intrinsic extinction risk and augmenting plant conservation efforts.

GuidedNOMe-seq quantifies chromatin states at single allele resolution for hundreds of custom regions in parallel

Since the introduction of next generation sequencing technologies, the field of epigenomics has evolved rapidly. However, most commonly used assays are enrichment-based methods and thus only semi-quantitative. Nucleosome occupancy and methylome sequencing (NOMe-seq) allows for quantitative inference of chromatin states with single locus resolution, but this requires high sequencing depth and is therefore prohibitively expensive to routinely apply to organisms with large genomes. To overcome this limitation, we introduce guidedNOMe-seq, where we combine NOMe profiling with large scale sgRNA synthesis and Cas9-mediated region-of-interest (ROI) liberation. To facilitate quantitative comparisons between multiple samples, we additionally develop an R package to standardize differential analysis of any type of NOMe-seq data. We extensively benchmark guidedNOMe-seq in a proof-of-concept study, dissecting the interplay of ChAHP and CTCF on chromatin. In summary we present a cost-effective, scalable, and customizable target enrichment extension to the existing NOMe-seq protocol allowing genome-scale quantification of nucleosome occupancy and transcription factor binding at single allele resolution.

Potential Transcriptional Enhancers in Coronaviruses: From Infectious Bronchitis Virus to SARS-CoV-2.

Coronaviruses constitute a global threat to human and animal health. It is essential to investigate the long-distance RNA-RNA interactions that approximate remote regulatory elements in strategies, including genome circularization, discontinuous transcription, and transcriptional enhancers, aimed at the rapid replication of their large genomes, pathogenicity, and immune evasion. Based on the primary sequences and modeled RNA-RNA interactions of two experimentally defined coronaviral enhancers, we detected via an in silico primary and secondary structural analysis potential enhancers in various coronaviruses, from the phylogenetically ancient avian infectious bronchitis virus (IBV) to the recently emerged SARS-CoV-2. These potential enhancers possess a core duplex-forming region that could transition between closed and open states, as molecular switches directed by viral or host factors. The duplex open state would pair with remote sequences in the viral genome and modulate the expression of downstream crucial genes involved in viral replication and host immune evasion. Consistently, variations in the predicted IBV enhancer region or its distant targets coincide with cases of viral attenuation, possibly driven by decreased open reading frame (ORF)3a immune evasion protein expression. If validated experimentally, the annotated enhancer sequences could inform structural prediction tools and antiviral interventions.

CircHIPK3 nucleates IGF2BP2 and functions as a competing endogenous RNA

Circular RNAs represent a class of endogenous RNAs that regulate gene expression and influence cell biological decisions with implications for the pathogenesis of several diseases. Here, we disclose a novel gene-regulatory role of circHIPK3 by combining analyses of large genomics datasets and mechanistic cell biological follow-up experiments. Using time-course depletion of circHIPK3 and specific candidate RNA-binding proteins, we identify several perturbed genes by RNA sequencing analyses. Expression-coupled motif analyses identify an 11-mer motif within circHIPK3, which also becomes enriched in genes that are downregulated upon circHIPK3 depletion. By mining eCLIP datasets and combined with RNA immunoprecipitation assays, we demonstrate that the 11-mer motif constitutes a strong binding site for IGF2BP2 in bladder cancer cell lines. Our results suggest that circHIPK3 can sequester IGF2BP2 as a competing endogenous RNA (ceRNA), leading to target mRNA stabilization. As an example of a circHIPK3-regulated gene, we focus on the STAT3 mRNA as a specific substrate of IGF2BP2 and validate that manipulation of circHIPK3 regulates IGF2BP2-STAT3 mRNA binding and, thereby, STAT3 mRNA levels. Surprisingly, absolute copy number quantifications demonstrate that IGF2BP2 outnumbers circHIPK3 by orders of magnitude, which is inconsistent with a simple 1:1 ceRNA hypothesis. Instead, we show that circHIPK3 can nucleate multiple copies of IGF2BP2, potentially via phase separation, to produce IGF2BP2 condensates. Our results support a model where a few cellular circHIPK3 molecules can induce IGF2BP2 condensation, thereby regulating key factors for cell proliferation.

CircHIPK3 nucleates IGF2BP2 and functions as a competing endogenous RNA.

Circular RNAs represent a class of endogenous RNAs that regulate gene expression and influence cell biological decisions with implications for the pathogenesis of several diseases. Here, we disclose a novel gene-regulatory role of circHIPK3 by combining analyses of large genomics datasets and mechanistic cell biological follow-up experiments. Using time-course depletion of circHIPK3 and specific candidate RNA-binding proteins, we identify several perturbed genes by RNA sequencing analyses. Expression-coupled motif analyses identify an 11-mer motif within circHIPK3, which also becomes enriched in genes that are downregulated upon circHIPK3 depletion. By mining eCLIP datasets and combined with RNA immunoprecipitation assays, we demonstrate that the 11-mer motif constitutes a strong binding site for IGF2BP2 in bladder cancer cell lines. Our results suggest that circHIPK3 can sequester IGF2BP2 as a competing endogenous RNA (ceRNA), leading to target mRNA stabilization. As an example of a circHIPK3-regulated gene, we focus on the STAT3 mRNA as a specific substrate of IGF2BP2 and validate that manipulation of circHIPK3 regulates IGF2BP2-STAT3 mRNA binding and, thereby, STAT3 mRNA levels. Surprisingly, absolute copy number quantifications demonstrate that IGF2BP2 outnumbers circHIPK3 by orders of magnitude, which is inconsistent with a simple 1:1 ceRNA hypothesis. Instead, we show that circHIPK3 can nucleate multiple copies of IGF2BP2, potentially via phase separation, to produce IGF2BP2 condensates. Our results support a model where a few cellular circHIPK3 molecules can induce IGF2BP2 condensation, thereby regulating key factors for cell proliferation.

AMRomics: a scalable workflow to analyze large microbial genome collections

Whole genome analysis for microbial genomics is critical to studying and monitoring antimicrobial resistance strains. The exponential growth of microbial sequencing data necessitates a fast and scalable computational pipeline to generate the desired outputs in a timely and cost-effective manner. Recent methods have been implemented to integrate individual genomes into large collections of specific bacterial populations and are widely employed for systematic genomic surveillance. However, they do not scale well when the population expands and turnaround time remains the main issue for this type of analysis. Here, we introduce AMRomics, an optimized microbial genomics pipeline that can work efficiently with big datasets. We use different bacterial data collections to compare AMRomics against competitive tools and show that our pipeline can generate similar results of interest but with better performance. The software is open source and is publicly available at https://github.com/amromics/amromics under an MIT license.

Population genomics of emerging multidrug-resistant Salmonella derby from pork and human in Guangdong, China

Salmonella is a foodborne pathogen that poses a severe public health threat. Salmonella enterica serovar Derby (S. Derby) is one of the top five serovars isolated from humans and animals, and it carries the potential risk of transmission to humans through consumption of pork and related products. This study isolated 98 S. Derby, mainly from pork, and clinical samples in Guangdong to investigate and compare their molecular characteristics, phylogenetic relationships, transmission routes, antimicrobial resistance phenotypes, genotypes, and virulence profiles. Whole-genome sequencing (WGS) revealed that all S. Derby isolates belonged to ST40 and displayed a relatively large core genome. Phylogenomic analysis indicated the possibility of clonal transmission of S. Derby not only between pork and humans but also across geographic areas in Guangdong over an extended period of time. Discriminant analysis of principal components (DAPC) model showed that a significant proportion of S. Derby isolates retrieved from multiple hosts originated from pork. Furthermore, similar antimicrobial resistance phenotypes and genotypes were found in S. Derby strains recovered from animal products and humans. All S. Derby isolates were found to contain numerous virulence genes, primarily related to fimbrial adherence determinants and secretion system components. This study offers valuable insights into the transmission routes of S. Derby between pork and humans and provides essential data for public health management.

Genome assembly catalog for species in the Japanese Red List: unlocking endangered biodiversity through genomic inventory.

Improvements in DNA sequencing technology are allowing the dramatic increase of whole genome data for a wide variety of species. Such genome sequence data can assist the monitoring of intraspecific genetic diversity, but is often lacking for threatened species. In this project, we focused on the national Red List, a catalog of extinct and threatened species, issued by the Japanese government. We combined the data included in it with the record of genome assembly in NCBI and tabulated the assembly availability of the species in the list. The combined data shows a low percentage (2.1%) of the availability of whole genome sequence data for the taxa ranked on the Japanese Red List as well as a strong bias towards mammals and birds in Animalia and vascular plants in Plantae. Our data presentation highlights potential systematic limitations in genome sequencing (e.g., budget for sequencing large genomes of amphibians) and instructs future policies including which taxon needs more effort for genome sequencing. The resultant tables are available in the original website https://treethinkers.nig.ac.jp/redlist/ and are regularly updated.

Genomic Epidemiology of Rift Valley Fever Virus Involved in the 2018 and 2022 Outbreaks in Livestock in Rwanda.

Rift Valley fever (RVF), a mosquito-borne transboundary zoonosis, was first confirmed in Rwanda's livestock in 2012 and since then sporadic cases have been reported almost every year. In 2018, the country experienced its first large outbreak, which was followed by a second one in 2022. To determine the circulating virus lineages and their ancestral origin, two genome sequences from the 2018 outbreak, and thirty-six, forty-one, and thirty-eight sequences of small (S), medium (M), and large (L) genome segments, respectively, from the 2022 outbreak were generated. All of the samples from the 2022 outbreak were collected from slaughterhouses. Both maximum likelihood and Bayesian-based phylogenetic analyses were performed. The findings showed that RVF viruses belonging to a single lineage, C, were circulating during the two outbreaks, and shared a recent common ancestor with RVF viruses isolated in Uganda between 2016 and 2019, and were also linked to the 2006/2007 largest East Africa RVF outbreak reported in Kenya, Tanzania, and Somalia. Alongside the wild-type viruses, genetic evidence of the RVFV Clone 13 vaccine strain was found in slaughterhouse animals, demonstrating a possible occupational risk of exposure with unknown outcome for people working in meat-related industry. These results provide additional evidence of the ongoing wide spread of RVFV lineage C in Africa and emphasize the need for an effective national and international One Health-based collaborative approach in responding to RVF emergencies.

Analysis of 14,392 whole genomes reveals 3.5% of Qataris carry medically actionable variants.

Arabic populations are underrepresented in large genome projects; therefore, the frequency of clinically actionable variants among Arabs is largely unknown. Here, we investigated genetic variation in 14,392 whole genomes from the Qatar Genome Program (QGP) across the list of 78 actionable genes (v3.1) determined by the American College of Medical Genetics and Genomics (ACMG). Variants were categorized into one of the following groups: (1) Pathogenic (P), (2) Likely pathogenic (LP), and (3) Rare variants of uncertain significance with evidence of pathogenicity. For the classification, we used variant databases, effect predictors, and the disease-relevant phenotypes available for the cohort. Data on cardiovascular disease, cancer, and hypercholesterolemia allowed us to assess the disease-relevant phenotype association of rare missense variants. We identified 248 distinct variants in 50 ACMG genes that fulfilled our criteria to be included in one of the three groups affecting 1036 genotype-positive participants of the QGP cohort. The most frequent variants were in TTN, followed by RYR1 and ATP7B. The prevalence of reportable secondary findings was 3.5%. A further 46 heterozygous variants in six genes with an autosomal recessive mode of inheritance were detected in 200 individuals, accounting for an additional 1.4%. Altogether, they affect 5% of the population. Due to the high consanguinity rate in the QGP cohort (28% in spouses and 60% in parents), P and LP variants both in genes with dominant and recessive inheritance are important for developing better treatment options and preventive strategies in Qatar and the Arabic population of the Middle East.

Immunological characteristics of a recombinant alphaherpesvirus with an envelope-embedded Cap protein of circovirus.

Variant pseudorabies virus (PRV) is a newly emerged zoonotic pathogen that can cause human blindness. PRV can take advantage of its large genome and multiple non-essential genes to construct recombinant attenuated vaccines carrying foreign genes. However, a major problem is that the foreign genes in recombinant PRV are only integrated into the genome for independent expression, rather than assembled on the surface of virion. We reported a recombinant PRV with deleted gE/TK genes and an inserted porcine circovirus virus 2 (PCV2) Cap gene into the extracellular domain of the PRV gE gene using the Cre-loxP recombinant system combined with the CRISPR-Cas9 gene editing system. This recombinant PRV (PRV-Cap), with the envelope-embedded Cap protein, exhibits a similar replication ability to its parental virus. An immunogenicity assay revealed that PRV-Cap immunized mice have 100% resistance to lethal PRV and PCV2 attacks. Neutralization antibody and ELISPOT detections indicated that PRV-Cap can enhance neutralizing antibodies to PRV and produce IFN-γ secreting T cells specific for both PRV and PCV2. Immunological mechanistic investigation revealed that initial immunization with PRV-Cap stimulates significantly early activation and expansion of CD69+ T cells, promoting the activation of CD4 Tfh cell dependent germinal B cells and producing effectively specific effector memory T and B cells. Booster immunization with PRV-Cap recalled the activation of PRV-specific IFN-γ+IL-2+CD4+ T cells and IFN-γ+TNF-α+CD8+ T cells, as well as PCV2-specific IFN-γ+TNF-α+CD8+ T cells. Collectively, our data suggested an immunological mechanism in that the recombinant PRV with envelope-assembled PCV2 Cap protein can serve as an excellent vaccine candidate for combined immunity against PRV and PCV2, and provided a cost-effective method for the production of PRV- PCV2 vaccine.

GeCKO: user-friendly workflows for genotyping complex genomes using target enrichment capture. A use case on the large tetraploid durum wheat genome

BackgroundGenotyping of individuals plays a pivotal role in various biological analyses, with technology choice influenced by multiple factors including genomic constraints, number of targeted loci and individuals, cost considerations, and the ease of sample preparation and data processing. Target enrichment capture of specific polymorphic regions has emerged as a flexible and cost-effective genomic reduction method for genotyping, especially adapted to the case of very large genomes. However, this approach necessitates complex bioinformatics treatment to extract genotyping data from raw reads. Existing workflows predominantly cater to phylogenetic inference, leaving a gap in user-friendly tools for genotyping analysis based on capture methods. In response to these challenges, we have developed GeCKO (Genotyping Complexity Knocked-Out). To assess the effectiveness of combining target enrichment capture with GeCKO, we conducted a case study on durum wheat domestication history, involving sequencing, processing, and analyzing variants in four relevant durum wheat groups.ResultsGeCKO encompasses four distinct workflows, each designed for specific steps of genomic data processing: (i) read demultiplexing and trimming for data cleaning, (ii) read mapping to align sequences to a reference genome, (iii) variant calling to identify genetic variants, and (iv) variant filtering. Each workflow in GeCKO can be easily configured and is executable across diverse computational environments. The workflows generate comprehensive HTML reports including key summary statistics and illustrative graphs, ensuring traceable, reproducible results and facilitating straightforward quality assessment. A specific innovation within GeCKO is its 'targeted remapping' feature, specifically designed for efficient treatment of targeted enrichment capture data. This process consists of extracting reads mapped to the targeted regions, constructing a smaller sub-reference genome, and remapping the reads to this sub-reference, thereby enhancing the efficiency of subsequent steps.ConclusionsThe case study results showed the expected intra-group diversity and inter-group differentiation levels, confirming the method's effectiveness for genotyping and analyzing genetic diversity in species with complex genomes. GeCKO streamlined the data processing, significantly improving computational performance and efficiency. The targeted remapping enabled straightforward SNP calling in durum wheat, a task otherwise complicated by the species' large genome size. This illustrates its potential applications in various biological research contexts.

Large Language Models and Genomics for Summarizing the Role of microRNA in Regulating mRNA Expression.

microRNA (miRNA)-messenger RNA (mRNA or gene) interactions are pivotal in various biological processes, including the regulation of gene expression, cellular differentiation, proliferation, apoptosis, and development, as well as the maintenance of cellular homeostasis and pathogenesis of numerous diseases, such as cancer, cardiovascular diseases, neurological disorders, and metabolic conditions. Understanding the mechanisms of miRNA-mRNA interactions can provide insights into disease mechanisms and potential therapeutic targets. However, extracting these interactions efficiently from a huge collection of published articles in PubMed is challenging. In the current study, we annotated a miRNA-mRNA Interaction Corpus (MMIC) and used it for evaluating the performance of a variety of machine learning (ML) models, deep learning-based transformer (DLT) models, and large language models (LLMs) in extracting the miRNA-mRNA interactions mentioned in PubMed. We used the genomics approaches for validating the extracted miRNA-mRNA interactions. Among the ML, DLT, and LLM models, PubMedBERT showed the highest precision, recall, and F-score, with all equal to 0.783. Among the LLM models, the performance of Llama-2 is better when compared to others. Llama 2 achieved 0.56 precision, 0.86 recall, and 0.68 F-score in a zero-shot experiment and 0.56 precision, 0.87 recall, and 0.68 F-score in a three-shot experiment. Our study shows that Llama 2 achieves better recall than ML and DLT models and leaves space for further improvement in terms of precision and F-score.

Rapid discrimination of African swine fever virus nucleic acid and virions using BenzoNuclease

African swine fever (ASF) is an acute and severe infectious disease caused by the African Swine Fever Virus (ASFV). ASFV exhibits significant resistance and stability in the environment, which, coupled with its double-stranded DNA and large genome, predisposes it to contaminate laboratory samples. This characteristic can lead to false-positive results in swine farm settings even days after disinfection, as detectable through polymerase chain reaction (PCR) or real-time fluorescent quantitative PCR (qPCR) assays. To meet the demand for efficient clinical methods capable of discriminating between ASFV nucleic acid and ASFV virions, this study aims to ascertain the efficacy of the nuclease “BenzoNuclease” in distinguishing ASFV nucleic acid (ASFV-DNA) from ASFV virions. BenzoNuclease is a versatile nucleic acid enzyme with the capacity to degrade nearly all forms of DNA and RNA. Initially, this research established a highly sensitive general PCR detection method for ASFV. Subsequently, a positive control was constructed using the M13 bacteriophage to substitute for active ASFV, facilitating the development of an improved qPCR method. It is important to note that common disinfectants have the potential to deactivate BenzoNuclease. However, in an environment simulating actual production applications, residual disinfectants do not interfere with the enzymatic efficacy of BenzoNuclease, thus not affecting the detection capabilities of this method. Positive clinical samples from pig farms, upon testing with the improved method, revealed that three samples were positive, indicating the presence of viral particles, while the remaining samples were negative, indicating the presence of nucleic acids. This provides an additional new option for sample testing in pig farms.

Viroporin-like activity of the hairpin transmembrane domain of African swine fever virus B169L protein.

African swine fever virus (ASFV) is the causative agent of a contagious disease affecting wild and domestic swine. The function of B169L protein, as a potential integral structural membrane protein, remains to be experimentally characterized. Using state-of-the-art bioinformatics tools, we confirm here earlier predictions indicating the presence of an integral membrane helical hairpin, and further suggest anchoring of this protein to the ER membrane, with both terminal ends facing the lumen of the organelle. Our evolutionary analysis confirmed the importance of purifying selection in the preservation of the identified domains during the evolution of B169L in nature. Also, we address the possible function of this hairpin transmembrane domain (HTMD) as a class IIA viroporin. Expression of GFP fusion proteins in the absence of a signal peptide supported B169L insertion into the ER as a Type III membrane protein and the formation of oligomers therein. Overlapping peptides that spanned the B169L HTMD were reconstituted into ER-like membranes and the adopted structures analyzed by infrared spectroscopy. Consistent with the predictions, B169L transmembrane sequences adopted α-helical conformations in lipid bilayers. Moreover, single vesicle permeability assays demonstrated the assembly of lytic pores in ER-like membranes by B169L transmembrane helices, a capacity confirmed by ion-channel activity measurements in planar bilayers. Emphasizing the relevance of these observations, pore-forming activities were not observed in the case of transmembrane helices derived from EP84R, another ASFV protein predicted to anchor to membranes through a α-helical HTMD. Overall, our results support predictions of viroporin-like function for the B169L HTMD.IMPORTANCEAfrican swine fever (ASF), a devastating disease affecting domestic swine, is widely spread in Eurasia, producing significant economic problems in the pork industry. Approaches to prevent/cure the disease are mainly restricted to the limited information concerning the role of most of the genes encoded by the large (160-170 kba) virus genome. In this report, we present the experimental data on the functional characterization of the African swine fever virus (ASFV) gene B169L. Data presented here indicates that the B169L gene encodes for an essential membrane-associated protein with a viroporin function.

Deciphering Pteronia's evolution in the Cape Floristic Region: A comprehensive study disputes polyploid deficiency and affirms diploid radiation.

The Greater Cape Floristic Region (GCFR) is renowned for its exceptional biodiversity, accommodating over 11 000 plant species, notable degree of endemism, and substantial diversification within limited plant lineages, a phenomenon ascribed to historical radiation events. While both abiotic and biotic factors contribute to this diversification, comprehensive genomic alterations, recognized as pivotal in the diversification of angiosperms, are perceived as uncommon. This investigation focuses on the genus Pteronia, a prominent representative of the Asteraceae family in the GCFR. Employing NGS-based HybSeq and RADSeq methodologies, flow cytometry, karyology, and ecological modeling, we scrutinize the intricacies of its polyploid evolution. Phylogenetic reconstructions using 951 low-copy nuclear genes confirm Pteronia as a well-supported, distinct clade within the tribe Astereae. The ingroup displays a structure indicative of rapid radiation likely antedating polyploid establishment, with the two main groups demarcated by their presence or absence in the fynbos biome. Genome size analysis encompasses 1293 individuals across 347 populations, elucidating significant variation ranging from 6.1 to 34.2 pg (2C-value). Pteronia demonstrates substantially large genome sizes within Astereae and phanerophytes. Polyploidy is identified in 31% of the studied species, with four discerned ploidy levels (2x, 4x, 6x, 8x). Cytotypes exhibit marked distinctions in environmental traits, influencing their distribution across biomes and augmenting their niche differentiation. These revelations challenge the presumed scarcity of polyploidy in the Cape flora, underscoring the imperative need for detailed population studies. The intricate evolutionary history of Pteronia, characterized by recent polyploidy and genome size variation, contributes substantially to the comprehension of diversification patterns within the GCFR biodiversity hotspot.

Repeated horizontal acquisition of lagriamide-producing symbionts in Lagriinae beetles.

Microbial symbionts associate with multicellular organisms on a continuum from facultative associations to mutual codependency. In some of the oldest intracellular symbioses there is exclusive vertical symbiont transmission, and co-diversification of symbiotic partners over millions of years. Such symbionts often undergo genome reduction due to low effective population sizes, frequent population bottlenecks, and reduced purifying selection. Here, we describe multiple independent acquisition events of closely related defensive symbionts followed by genome erosion in a group of Lagriinae beetles. Previous work in Lagria villosa revealed the dominant genome-eroded symbiont of the genus Burkholderia produces the antifungal compound lagriamide and protects the beetle's eggs and larvae from antagonistic fungi. Here, we use metagenomics to assemble 11 additional genomes of lagriamide-producing symbionts from seven different host species within Lagriinae from five countries, to unravel the evolutionary history of this symbiotic relationship. In each host species, we detected one dominant genome-eroded Burkholderia symbiont encoding the lagriamide biosynthetic gene cluster (BGC). Surprisingly, however, we did not find evidence for host-symbiont co-diversification, or for a monophyly of the lagriamide-producing symbionts. Instead, our analyses support at least four independent acquisition events of lagriamide-encoding symbionts and subsequent genome erosion in each of these lineages. By contrast, a clade of plant-associated relatives retained large genomes but secondarily lost the lagriamide BGC. In conclusion, our results reveal a dynamic evolutionary history with multiple independent symbiont acquisitions characterized by high degree of specificity. They highlight the importance of the specialized metabolite lagriamide for the establishment and maintenance of this defensive symbiosis.

Copy number variation introduced by a massive mobile element facilitates global thermal adaptation in a fungal wheat pathogen

Copy number variation (CNV) can drive rapid evolution in changing environments. In microbial pathogens, such adaptation is a key factor underpinning epidemics and colonization of new niches. However, the genomic determinants of such adaptation remain poorly understood. Here, we systematically investigate CNVs in a large genome sequencing dataset spanning a worldwide collection of 1104 genomes from the major wheat pathogen Zymoseptoria tritici. We found overall strong purifying selection acting on most CNVs. Genomic defense mechanisms likely accelerated gene loss over episodes of continental colonization. Local adaptation along climatic gradients was likely facilitated by CNVs affecting secondary metabolite production and gene loss in general. One of the strongest loci for climatic adaptation is a highly conserved gene of the NAD-dependent Sirtuin family. The Sirtuin CNV locus localizes to an ~68-kb Starship mobile element unique to the species carrying genes highly expressed during plant infection. The element has likely lost the ability to transpose, demonstrating how the ongoing domestication of cargo-carrying selfish elements can contribute to selectable variation within populations. Our work highlights how standing variation in gene copy numbers at the global scale can be a major factor driving climatic and metabolic adaptation in microbial species.

MyxoPortal: a database of myxobacterial genomic features.

Myxobacteria are predatory bacteria with antimicrobial activity, utilizing complex mechanisms to kill their prey and assimilate their macromolecules. Having large genomes encoding hundreds of secondary metabolites, hydrolytic enzymes and antimicrobial peptides, these organisms are widely studied for their antibiotic potential. MyxoPortal is a comprehensive genomic database hosting 262 genomes of myxobacterial strains. Datasets included provide genome annotations with gene locations, functions, amino acids and nucleotide sequences, allowing analysis of evolutionary and taxonomical relationships between strains and genes. Biosynthetic gene clusters are identified by AntiSMASH, and dbAMP-generated antimicrobial peptide sequences are included as a resource for novel antimicrobial discoveries, while curated datasets of CRISPR/Cas genes, regulatory protein sequences, and phage associated genes give useful insights into each strain's biological properties. MyxoPortal is an intuitive open-source database that brings together application-oriented genomic features that can be used in taxonomy, evolution, predation and antimicrobial research. MyxoPortal can be accessed at http://dicsoft1.physics.iisc.ac.in/MyxoPortal/. Database URL: http://dicsoft1.physics.iisc.ac.in/MyxoPortal/. Graphical Abstract.