Genome Evolution Research Articles

Dissolved organic matter, calcium ion and extracellular polymeric substances on living associated bacteria of Microcystis colony are crucial for unicellular Microcystis to efficiently form colonies

Microcystis typically forms colonies under natural conditions, which contributes to occurrence and prevalence of algal blooms. The colonies consist of Microcystis and associated bacteria (AB), embedded in extracellular polymeric substances (EPS). Previous studies indicate that AB can induce Microcystis to form colonies, however the efficiency is generally low and results in a uniform morphotype. In this study, by using filtrated natural water, several AB strains induced unicellular M. aeruginosa to form colonies resembling several Microcystis morphotypes. The mechanisms were investigated with Methylobacterium sp. Z5. Ca2+ was necessary for Z5 to induce Microcystis to form colonies, while dissolved organic matters (DOM) facilitated AB to agglomerate Microcystis to form large colonies. EPS of living Z5, mainly the aromatic protein components, played a key role in colony induction. Z5 initially aggregated Microcystis via the bridging effects of Ca2+ and DOM, followed by the induction of EPS synthesis and secretion in Microcystis. In this process, the colony forming mode shifted from cell adhesion to a combination of cell adhesion and cell division. Intriguingly, Z5 drove the genomic rearrangement of Microcystis by upregulating some transposase genes. This study unveiled a novel mechanism about Microcystis colony formation and identified a new driver of Microcystis genomic evolution.

Moving toward the Inclusion of Epigenomics in Bacterial Genome Evolution: Perspectives and Challenges.

The universality of DNA methylation as an epigenetic regulatory mechanism belongs to all biological kingdoms. However, while eukaryotic systems have been the primary focus of DNA methylation studies, the molecular mechanisms in prokaryotes are less known. Nevertheless, DNA methylation in prokaryotes plays a pivotal role in many cellular processes such as defense systems against exogenous DNA, cell cycle dynamics, and gene expression, including virulence. Thanks to single-molecule DNA sequencing technologies, genome-wide identification of methylated DNA is becoming feasible on a large scale, providing the possibility to investigate more deeply the presence, variability, and roles of DNA methylation. Here, we present an overview of the multifaceted roles of DNA methylation in prokaryotes and suggest research directions and tools which can enable us to better understand the contribution of DNA methylation to prokaryotic genome evolution and adaptation. In particular, we emphasize the need to understand the presence and role of transgenerational inheritance, as well as the impact of epigenomic signatures on adaptation and genome evolution. Research directions and the importance of novel computational tools are underlined.

Revealing the Satellite DNA Content in Ancistrus sp. (Siluriformes: Loricariidae) by Genomic and Bioinformatic Analysis

Introduction: Eukaryotic genomes are composed of simple, repetitive sequences, including satellite DNAs (satDNA), which are noncoding sequences arranged in tandem arrays. These sequences play a crucial role in genomic functions and innovations, influencing processes such as the maintenance of nuclear material, the formation of heterochromatin and the differentiation of sex chromosomes. In this genomic era, advances in next-generation sequencing and bioinformatics tools have facilitated the exhaustive cataloging of repetitive elements in genomes, particularly in non-model species. This study focuses on the satDNA content of Ancistrus sp., a diverse species of fish from the Loricariidae family. The genus Ancistrus shows significant karyotypic evolution, with extensive variability from the ancestral diploid number. Methods: By means of bioinformatic approaches, 40 satDNA families in Ancistrus sp., constituting 5.19% of the genome were identified. Analysis of the abundance and divergence landscape revealed diverse profiles, indicating recent amplification and homogenization of these satDNA sequences. Results: The most abundant satellite, AnSat1-142, constitutes 2.1% of the genome, while the least abundant, AnSat40-52, represents 0.0034%. The length of the monomer repeat varies from 16 to 142 base pairs, with an average length of 61 bp. These results contribute to understanding the genomic dynamics and evolution of satDNAs in Ancistrus sp. Conclusion: The study underscores the variability of satDNAs between fish species and provides valuable information on chromosome organization and the evolution of repetitive elements in non-model organisms.

Genome sequences of the first Autographiviridae phages infecting marine Roseobacter.

The ubiquitous and abundant marine phages play critical roles in shaping the composition and function of bacterial communities, impacting biogeochemical cycling in marine ecosystems. Autographiviridae is among the most abundant and ubiquitous phage families in the ocean. However, studies on the diversity and ecology of Autographiviridae phages in marine environments are restricted to isolates that infect SAR11 bacteria and cyanobacteria. In this study, ten new roseophages that infect marine Roseobacter strains were isolated from coastal waters. These new roseophages have a genome size ranging from 38 917 to 42 634 bp and G+C content of 44.6-50 %. Comparative genomics showed that they are similar to known Autographiviridae phages regarding gene content and architecture, thus representing the first Autographiviridae roseophages. Phylogenomic analysis based on concatenated conserved genes showed that the ten roseophages form three distinct subgroups within the Autographiviridae, and sequence analysis revealed that they belong to eight new genera. Finally, viromic read-mapping showed that these new Autographiviridae phages are widely distributed in global oceans, mostly inhabiting polar and estuarine locations. This study has expanded the current understanding of the genomic diversity, evolution and ecology of Autographiviridae phages and roseophages. We suggest that Autographiviridae phages play important roles in the mortality and community structure of roseobacters, and have broad ecological applications.

The Net of Life, a short story: Intricate patterns of gene flows across hundreds of extant genomes, all the way to LUCA

Over the past quarter-century, the field of evolutionary biology has been transformed by the emergence of complete genome sequences and the conceptual framework known as the 'Net of Life.' This paradigm shift challenges traditional notions of evolution as a tree-like process, emphasizing the complex, interconnected network of gene flow that may blur the boundaries between distinct lineages. In this context, gene loss, rather than horizontal gene transfer, is the primary driver of gene content, with vertical inheritance playing a principal role. The 'Net of Life' not only impacts our understanding of genome evolution but also has profound implications for classification systems, the rapid appearance of new traits, and the spread of diseases. Here, we explore the core tenets of the 'Net of Life' and its implications for genome-scale phylogenetic divergence, providing a comprehensive framework for further investigations in evolutionary biology.

The identification of retro-DNAs in primate genomes as DNA transposons mobilizing via retrotransposition

Background Mobile elements (MEs) constitute a major portion of the genome in primates and other higher eukaryotes, and they play important role in genome evolution and gene function. MEs can be divided into two fundamentally different classes: DNA transposons which transpose in the genome in a “cut-and-paste” style, and retrotransposons which propagate in a “copy-and-paste” fashion via a process involving transcription and reverse-transcription. In primate genomes, DNA transposons are mostly dead, while many retrotransposons are still highly active. We report here the identification of a new type of MEs, which we call “retro-DNAs”, for their combined characteristics of these two fundamentally different ME classes. Methods A comparative computational genomic approach was used to analyze the reference genome sequences of 10 primate species consisting of five apes, four monkeys, and marmoset. Results From our analysis, we identified a total of 1,750 retro-DNAs, representing 748 unique insertion events in the genomes of ten primate species including human. These retro-DNAs contain sequences of DNA transposons but lack the terminal inverted repeats (TIRs), the hallmark of DNA transposons. Instead, they show characteristics of retrotransposons, such as polyA tails, longer target-site duplications (TSDs), and the “TT/AAAA” insertion site motif, suggesting the use of the L1-based target-primed reverse transcription (TPRT) mechanism. At least 40% of these retro-DNAs locate into genic regions, presenting potentials for impacting gene function. More interestingly, some retro-DNAs, as well as their parent sites, show certain levels of expression, suggesting that they have the potential to create more retro-DNA copies in the present primate genomes. Conclusions Although small in number, the identification of these retro-DNAs reveals a new mechanism for propagating DNA transposons in primate genomes without active canonical DNA transposon activity. Our data also suggest that the TPRT machinery may transpose a wider variety of DNA sequences in the genomes.

Single-mitosis dissection of acute and chronic DNA mutagenesis and repair.

How chronic mutational processes and punctuated bursts of DNA damage drive evolution of the cancer genome is poorly understood. Here, we demonstrate a strategy to disentangle and quantify distinct mechanisms underlying genome evolution in single cells, during single mitoses and at single-strand resolution. To distinguish between chronic (reactive oxygen species (ROS)) and acute (ultraviolet light (UV)) mutagenesis, we microfluidically separate pairs of sister cells from the first mitosis following burst UV damage. Strikingly, UV mutations manifest as sister-specific events, revealing mirror-image mutation phasing genome-wide. In contrast, ROS mutagenesis in transcribed regions is reduced strand agnostically. Successive rounds of genome replication over persisting UV damage drives multiallelic variation at CC dinucleotides. Finally, we show that mutation phasing can be resolved to single strands across the entire genome of liver tumors from F1 mice. This strategy can be broadly used to distinguish the contributions of overlapping cancer relevant mutational processes.

Using the pan-genomic framework for the discovery of genomic islands in the haloarchaeon Halorubrum ezzemoulense.

In this study, we use pan-genomics to characterize the genomic variability of the widely dispersed halophilic archaeal species Halorubrum ezzemoulense (Hez). We include a multi-regional sampling of newly sequenced, high-quality draft genomes. The pan-genome graph of the species reveals 50 genomic islands that represent rare accessory genetic capabilities available to members. Most notably, we observe rearrangements that have led to the insertion/recombination/replacement of mutually exclusive genomic islands in equivalent genome positions ("homeocassettes"). These conflicting islands encode for similar functions, but homologs from islands located between the same core genes exhibit high divergence on the amino acid level, while the neighboring core genes are nearly identical. Both islands of a homeocassette often coexist in the same geographic location, suggesting that either island may be beyond the reach of selective sweeps and that these loci of divergence between Hez members are maintained and persist long term. This implies that subsections of the population have different niche preferences and rare metabolic capabilities. After an evaluation of the gene content in the homeocassettes, we speculate that these islands may play a role in the speciation, niche adaptability, and group selection dynamics in Hez. Though homeocassettes are first described in this study, similar replacements and divergence of genes on genomic islands have been previously reported in other Haloarchaea and distantly related Archaea, suggesting that homeocassettes may be a feature in a wide range of organisms outside of Hez.IMPORTANCEThis study catalogs the rare genes discovered in strains of the species Halorubrum ezzemoulense (Hez), an obligate halophilic archaeon, through the perspective of its pan-genome. These rare genes are often found to be arranged on islands that confer metabolic and transport functions and contain genes that have eluded previous studies. The discovery of divergent, but homologous islands occupying equivalent genome positions ("homeocassettes") in different genomes, reveals significant new information on genome evolution in Hez. Homeocassette pairs encode for similar functions, but their dissimilarity and distribution imply high rates of recombination, different specializations, and niche preferences in Hez. The coexistence of both islands of a homeocassette pair in multiple environments demonstrates that both islands are beyond the reach of selective sweeps and that these genome content differences between strains persist long term. The switch between islands through recombination under different environmental conditions may lead to a greater range of niche adaptability in Hez.

First telomere-to-telomere gapless assembly of the rice blast fungus Pyricularia oryzae

Rice blast caused by Pyricularia oryzae (syn., Magnaporthe oryzae) was one of the most destructive diseases of rice throughout the world. Genome assembly was fundamental to genetic variation identification and critically impacted the understanding of its ability to overcome host resistance. Here, we report a gapless genome assembly of rice blast fungus P. oryzae strain P131 using PacBio, Illumina and high throughput chromatin conformation capture (Hi-C) sequencing data. This assembly contained seven complete chromosomes (43,237,743 bp) and a circular mitochondrial genome (34,866 bp). Approximately 14.31% of this assembly carried repeat sequences, significantly greater than its previous assembled version. This assembly had a 99.9% complement in BUSCO evaluation. A total of 14,982 genes protein-coding genes were predicted. In summary, we assembled the first telomere-to-telomere gapless genome of P. oryzae, which would be a valuable genome resource for future research on the genome evolution and host adaptation.

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.

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 (&lt;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.

Mitotic gene regulation by the N-MYC-WDR5-PDPK1 nexus

BackgroundDuring mitosis the cell depends on proper attachment and segregation of replicated chromosomes to generate two identical progeny. In cancers defined by overexpression or dysregulation of the MYC oncogene this process becomes impaired, leading to genomic instability and tumor evolution. Recently it was discovered that the chromatin regulator WDR5—a critical MYC cofactor—regulates expression of genes needed in mitosis through a direct interaction with the master kinase PDPK1. However, whether PDPK1 and WDR5 contribute to similar mitotic gene regulation in MYC-overexpressing cancers remains unclear. Therefore, to characterize the influence of WDR5 and PDPK1 on mitotic gene expression in cells with high MYC levels, we performed a comparative transcriptomic analysis in neuroblastoma cell lines defined by MYCN-amplification, which results in high cellular levels of the N-MYC protein.ResultsUsing RNA-seq analysis, we identify the genes regulated by N-MYC and PDPK1 in multiple engineered CHP-134 neuroblastoma cell lines and compare them to previously published gene expression data collected in CHP-134 cells following inhibition of WDR5. We find that as expected N-MYC regulates a multitude of genes, including those related to mitosis, but that PDPK1 regulates specific sets of genes involved in development, signaling, and mitosis. Analysis of N-MYC- and PDPK1-regulated genes reveals a small group of commonly controlled genes associated with spindle pole formation and chromosome segregation, which overlap with genes that are also regulated by WDR5. We also find that N-MYC physically interacts with PDPK1 through the WDR5-PDPK1 interaction suggesting regulation of mitotic gene expression may be achieved through a N-MYC-WDR5-PDPK1 nexus.ConclusionsOverall, we identify a small group of genes highly enriched within functional gene categories related to mitotic processes that are commonly regulated by N-MYC, WDR5, and PDPK1 and suggest that a tripartite interaction between the three regulators may be responsible for setting the level of mitotic gene regulation in N-MYC amplified cell lines. This study provides a foundation for future studies to determine the exact mechanism by which N-MYC, WDR5, and PDPK1 converge on cell cycle related processes.

A chromosome-level genome reveals genome evolution and molecular basis of anthraquinone biosynthesis in Rheum palmatum

BackgroundRhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism.ResultsWe assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening.ConclusionOverall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies.

Mitochondrial genome variation and intergenomic sequence transfers in Hevea species.

Among the Hevea species, rubber tree (Hevea brasiliensis) is the most important source of natural rubber. In previous studies, we sequenced the complete nuclear and chloroplast genomes of Hevea species, providing an invaluable resource for studying their phylogeny, disease resistance, and breeding. However, given that plant mitochondrial genomes are more complex and more difficult to assemble than that of the other organelles, little is known about their mitochondrial genome, which limits the comprehensive understanding of Hevea genomic evolution. In this study, we sequenced and assembled the mitochondrial genomes of four Hevea species. The four mitochondrial genomes had consistent GC contents, codon usages and AT skews. However, there were significant differences in the genome lengths and sequence repeats. Specifically, the circular mitochondrial genomes of the four Hevea species ranged from 935,732 to 1,402,206 bp, with 34-35 unique protein-coding genes, 35-38 tRNA genes, and 6-13 rRNA genes. In addition, there were 17,294-46,552 bp intergenomic transfer fragments between the chloroplast and mitochondrial genomes, consisting of eight intact genes (psaA, rrn16S, tRNA-Val, rrn5S, rrn4.5S, tRNA-Arg, tRNA-Asp, and tRNA-Asn), intergenic spacer regions and partial gene sequences. The evolutionary position of Hevea species, crucial for understanding its adaptive strategies and relation to other species, was verified by phylogenetic analysis based on the protein-coding genes in the mitochondrial genomes of 21 Malpighiales species. The findings from this study not only provide valuable insights into the structure and evolution of the Hevea mitochondrial genome but also lay the foundation for further molecular, evolutionary studies, and genomic breeding studies on rubber tree and other Hevea species, thereby potentially informing conservation and utilization strategies.

Assembly and comparative analysis of the complete multichromosomal mitochondrial genome of Cymbidium ensifolium, an orchid of high economic and ornamental value

BackgroundOrchidaceae is one of the largest groups of angiosperms, and most species have high economic value and scientific research value due to their ornamental and medicinal properties. In China, Chinese Cymbidium is a popular ornamental orchid with high economic value and a long history. However, to date, no detailed information on the mitochondrial genome of any species of Chinese Cymbidium has been published.ResultsHere, we present the complete assembly and annotation of the mitochondrial genome of Cymbidium ensifolium (L.) Sw. The mitogenome of C. ensifolium was 560,647 bp in length and consisted of 19 circular subgenomes ranging in size from 21,995 bp to 48,212 bp. The genome encoded 35 protein-coding genes, 36 tRNAs, 3 rRNAs, and 3405 ORFs. Repeat sequence analysis and prediction of RNA editing sites revealed a total of 915 dispersed repeats, 162 simple repeats, 45 tandem repeats, and 530 RNA editing sites. Analysis of codon usage showed a preference for codons ending in A/T. Interorganellar DNA transfer was identified in 13 of the 19 chromosomes, with plastid-derived DNA fragments representing 6.81% of the C. ensifolium mitochondrial genome. The homologous fragments of the mitochondrial genome and nuclear genome were also analysed. Comparative analysis showed that the GC content was conserved, but the size, structure, and gene content of the mitogenomes varied greatly among plants with multichromosomal mitogenome structure. Phylogenetic analysis based on the mitogenomes reflected the evolutionary and taxonomic statuses of C. ensifolium. Interestingly, compared with the mitogenomes of Cymbidium lancifolium Hook. and Cymbidium macrorhizon Lindl., the mitogenome of C. ensifolium lost 8 ribosomal protein-coding genes.ConclusionIn this study, we assembled and annotated the mitogenome of C. ensifolium and compared it with the mitogenomes of other Liliidae and plants with multichromosomal mitogenome structures. Our findings enrich the mitochondrial genome database of orchid plants and reveal the rapid structural evolution of Cymbidium mitochondrial genomes, highlighting the potential for mitochondrial genes to help decipher plant evolutionary history.

Chromosome-level genome assembly and characterization of the Calophaca sinica genome.

Calophaca sinica is a rare plant endemic to northern China which belongs to the Fabaceae family and possesses rich nutritional value. To support the preservation of the genetic resources of this plant, we have successfully generated a high-quality genome of C. sinica (1.06 Gb). Notably, transposable elements (TEs) constituted ~73% of the genome, with long terminal repeat retrotransposons (LTR-RTs) dominating this group of elements (~54% of the genome). The average intron length of the C. sinica genome was noticeably longer than what has been observed for closely related species. The expansion of LTR-RTs and elongated introns emerged had the largest influence on the enlarged genome size of C. sinica in comparison to other Fabaceae species. The proliferation of TEs could be explained by certain modes of gene duplication, namely, whole genome duplication (WGD) and dispersed duplication (DSD). Gene family expansion, which was found to enhance genes associated with metabolism, genetic maintenance, and environmental stress resistance, was a result of transposed duplicated genes (TRD) and WGD. The presented genomic analysis sheds light on the genetic architecture of C. sinica, as well as provides a starting point for future evolutionary biology, ecology, and functional genomics studies centred around C. sinica and closely related species.

Plant genome evolution in the genus Eucalyptus is driven by structural rearrangements that promote sequence divergence.

Genomes have a highly organized architecture (nonrandom organization of functional and nonfunctional genetic elements within chromosomes) that is essential for many biological functions, particularly gene expression and reproduction. Despite the need to conserve genome architecture, a high level of structural variation has been observed within species. As species separate and diverge, genome architecture also diverges, becoming increasingly poorly conserved as divergence time increases. However, within plant genomes, the processes of genome architecture divergence are not well described. Here we use long-read sequencing and de novo assembly of 33 phylogenetically diverse, wild and naturally evolving Eucalyptus species, covering 1-50 million years of diverging genome evolution to measure genome architectural conservation and describe architectural divergence. The investigation of these genomes revealed that following lineage divergence, genome architecture is highly fragmented by rearrangements. As genomes continue to diverge, the accumulation of mutations and the subsequent divergence beyond recognition of rearrangements become the primary driver of genome divergence. The loss of syntenic regions also contribute to genome divergence but at a slower pace than that of rearrangements. We hypothesize that duplications and translocations are potentially the greatest contributors to Eucalyptus genome divergence.

Genomic evolution of human respiratory syncytial virus during a decade (2013–2023): bridging the path to monoclonal antibody surveillance

This study investigated the prevalence, genetic diversity, and evolution of human respiratory syncytial virus (HRSV) in Barcelona from 2013 to 2023. Respiratory specimens from patients with RTI suspicion at Hospital Universitari Vall d'Hebron were collected from October 2013 to May 2023 for laboratory-confirmation of respiratory viruses. Next-generation sequencing was performed in randomly-selected samples with Illumina technology. Phylogenetic analyses of whole genome sequences were performed with BEAST v1.10.4. Signals of selection and evolutionary pressures were inferred by population dynamics and evolutionary analyses. Mutations in major surface proteins were genetic and structurally characterised, emphasizing those within antigenic epitopes. Analyzing 139,625 samples, 5.3% were HRSV-positive (3008 HRSV-A, 3882 HRSV-B, 56 HRSV-A and -B, and 495 unsubtyped HRSV), with a higher prevalence observed in the paediatric population. Pandemic-related shifts in seasonal patterns returned to normal in 2022-2023. A total of 198 whole-genome sequences were obtained for HRSV-A (6.6% of the HRSV-A positive samples) belonging to GA2.3.5 lineage. For HRSV-B, 167 samples were sequenced (4.3% of the HRSV-B positive samples), belonging to GB5.0.2, GB5.0.4a and GB5.0.5a. HRSV-B exhibited a higher evolution rate. Post-SARS-CoV-2 pandemic, both subtypes showed increased evolutionary rates and decreased effective population size initially, followed by a sharp increase. Analyses indicated negative selective pressure on HRSV. Mutations in antigenic epitopes, including S276N and M274I in palivizumab-targeted site II, and I206M, Q209R, and S211N in nirsevimab-targeted site Ø, were identified. Particularly in the context of the large-scale use in 2023-2024 season of nirsevimab, continuous epidemiological and genomic surveillance is crucial.

Abstract NG02: Distinct genomic and immunologic tumor evolution in germline TP53-driven breast cancers

Abstract NG02: Distinct genomic and immunologic tumor evolution in germline TP53-driven breast cancers

Ecological Predictors of Organelle Genome Evolution: Phylogenetic Correlations with Taxonomically Broad, Sparse, Unsystematized Data.

Comparative analysis of variables across phylogenetically linked observations can reveal mechanisms and insights in evolutionary biology. As the taxonomic breadth of the sample of interest increases, challenges of data sparsity, poor phylogenetic resolution, and complicated evolutionary dynamics emerge. Here, we investigate a cross-eukaryotic question where all these problems exist: which organismal ecology features are correlated with gene retention in mitochondrial and chloroplast DNA (organelle DNA or oDNA). Through a wide palette of synthetic control studies, we first characterize the specificity and sensitivity of a collection of parametric and non-parametric phylogenetic comparative approaches to identify relationships in the face of such sparse and awkward datasets. This analysis is not directly focused on oDNA, and so provides generalizable insights into comparative approaches with challenging data. We then combine and curate ecological data coupled to oDNA genome information across eukaryotes, including a new semi-automated approach for gathering data on organismal traits from less systematized open-access resources including encyclopedia articles on species and taxa. The curation process also involved resolving several issues with existing datasets, including enforcing the clade-specificity of several ecological features and fixing incorrect annotations. Combining this unique dataset with our benchmarked comparative approaches, we confirm support for several known links between organismal ecology and organelle gene retention, identify several previously unidentified relationships constituting possible ecological contributors to oDNA genome evolution, and provide support for a recently hypothesized link between environmental demand and oDNA retention. We, with caution, discuss the implications of these findings for organelle evolution and of this pipeline for broad comparative analyses in other fields.