Genome Of Line Research Articles

Climate-driven succession in marine microbiome biodiversity and biogeochemical function.

Seasonal and El Niño-Southern Oscillation (ENSO) warming result in similar ocean changes as predicted with climate change. Climate-driven environmental cycles have strong impacts on microbiome diversity, but impacts on microbiome function are poorly understood. We quantified changes in microbial genomic diversity and functioning over 11 years covering seasonal and ENSO cycles at a coastal site in the southern California Current. We observed seasonal oscillations between large genome lineages during cold, nutrient rich conditions in winter and spring versus small genome lineages, including Prochlorococcus and Pelagibacter , in summer and fall. Parallel interannual changes separated communities depending on ENSO condition. Biodiversity shifts translated into clear oscillations in microbiome functional potential. Ocean warming induced an ecosystem with less iron but more macronutrient stress genes, depressed organic carbon degradation potential and biomass, and elevated carbon-to-nutrient biomass ratios. The consistent microbial response observed across time-scales points towards large climate-driven changes in marine ecosystems and biogeochemical cycles.

Vancomycin-resistant Staphylococcus aureus (VRSA) can overcome the cost of antibiotic resistance and may threaten vancomycin's clinical durability.

Vancomycin has proven remarkably durable to resistance evolution by Staphylococcus aureus despite widespread treatment with vancomycin in the clinic. Only 16 cases of vancomycin-resistant S. aureus (VRSA) have been documented in the United States. It is thought that the failure of VRSA to spread is partly due to the fitness cost imposed by the vanA operon, which is the only known means of high-level resistance. Here, we show that the fitness cost of vanA-mediated resistance can be overcome through laboratory evolution of VRSA in the presence of vancomycin. Adaptation to vancomycin imposed a tradeoff such that fitness in the presence of vancomycin increased, while fitness in its absence decreased in evolved lineages. Comparing the genomes of vancomycin-exposed and vancomycin-unexposed lineages pinpointed the D-alanine:D-alanine ligase gene (ddl) as the target of loss-of-function mutations, which were associated with the observed fitness tradeoff. Vancomycin-exposed lineages exhibited vancomycin dependence and abnormal colony morphology in the absence of drug, which were associated with mutations in ddl. However, further evolution of vancomycin-exposed lineages in the absence of vancomycin enabled some evolved lineages to escape this fitness tradeoff. Many vancomycin-exposed lineages maintained resistance in the absence of vancomycin, unlike their ancestral VRSA strains. These results indicate that VRSA might be able to compensate for the fitness deficit associated with vanA-mediated resistance, which may pose a threat to the prolonged durability of vancomycin in the clinic. Our results also suggest vancomycin treatment should be immediately discontinued in patients after VRSA is identified to mitigate potential adaptations.

Investigation of SARS-CoV-2 Infection among Companion Animals in Households with Confirmed Human COVID-19 Cases.

We aimed to characterize SARS-CoV-2 infection in companion animals living in households with COVID-19-positive people and understand the dynamics surrounding how these animals become infected. Public health investigators contacted households with at least one confirmed, symptomatic person with COVID-19 for study recruitment. Blood, nasal, and rectal swab specimens were collected from pet dogs and cats and a questionnaire was completed. Specimens were tested for SARS-CoV-2 by RT-PCR, and for neutralizing antibodies; genomic sequencing was performed on viral-positive samples. A total of 36.4% of 110 pets enrolled had evidence of infection with SARS-CoV-2. Pets were more likely to test positive if the pet was immunocompromised, and if more than one person in the home was positive for COVID-19. Among 12 multi-pet households where at least one pet was positive, 10 had at least one other pet test positive. Whole-genome sequencing revealed the genomes of viral lineages circulating in the community during the time of sample collection. Our findings suggest a high likelihood of viral transmission in households with multiple pets and when pets had very close interactions with symptomatic humans. Further surveillance studies are needed to characterize how new variants impact animals and to understand opportunities for infection and spillover in susceptible species.

Is the persistence of Listeria monocytogenes in food processing facilities and its resistance to pathogen intervention linked to its phylogeny?

The foodborne pathogen Listeria monocytogenes is differentiated into four distinct lineages which differ in their virulence. It remains unknown, however, whether the four lineages also differ with respect to their ability to persist in food processing facilities, their resistance to high pressure, a preservation method that is used commercially for Listeria control on ready-to-eat meats, and their ability to form biofilms. This study aimed to determine differences in the pressure resistance and biofilm formation of 59 isolates of L. monocytogenes representing lineages I and II. Furthermore, the genetic similarity of 9 isolates of L. monocytogenes that were obtained from a meat processing facility over a period of 1 year and of 20 isolates of L. monocytogenes from food processing facilities was analyzed to assess whether the ability of the lineages of L. monocytogenes to persist in these facilities differs. Analysis of 386 genomes with respect to the source of isolation revealed that genomes of lineage II are over-represented in meat isolates when compared with clinical isolates. Of the 38 strains of Lm. monocytogenes that persisted in food processing facilities (this study or published studies), 31 were assigned to lineage II. Isolates of lineage I were more resistant to treatments at 400 to 600 MPa. The thickness of biofilms did not differ between lineages. In conclusion, strains of lineage II are more likely to persist in food processing facilities while strains of lineage I are more resistant to high pressure.IMPORTANCEListeria monocytogenes substantially contributes to the mortality of foodborne disease in developed countries. The virulence of strains of four lineages of L. monocytogenes differs, indicating that risks associated with the presence of L. monocytogenes are lineage specific. Our study extends the current knowledge by documentation that the lineage-level phylogeny of L. monocytogenes plays a role in the source of isolation, in the persistence in food processing facilities, and in the resistance to pathogen intervention technologies. In short, the control of risks associated with the presence of L. monocytogenes in food is also lineage specific. Understanding the route of contamination L. monocytogenes is an important factor to consider when designing improved control measures.

Circular genetic structure of the Abies nephrolepis species complex shaped by the circular landform of Northeast Asia

AbstractAimPinaceae have been noted for their tendency towards reticulation as a result of interspecific hybridization. Here, we demonstrated the phylogeographic dynamics of the native Northeast Asian sub‐alpine conifer, Abies nephrolepis species complex (ANSC), evolving in circular overlaps along the Northeast Asian landform, which functions as a corridor with the sea, thus acting as a geographic barrier.LocationNortheast Asia: the Korean Peninsula, Japanese Archipelago, Russian Far East and northeastern China.TaxonAbies nephrolepis species complex (Family Pinaceae, Genus Abies, Section Balsamea, Abies nephrolepis, A. koreana, A. veitchii and A. sachalinensis).MethodsA total of 728 individuals from 38 ANSC populations were analysed using multiplexed inter‐simple sequence repeat genotyping by sequencing to capture variations in biparental nuclear genomes. Eight mitochondrial regions and eight chloroplast regions of each individual were sequenced using the MiSeq platform and Sanger sequencing. Species distribution models were generated.ResultsBayesian clustering with 507 nuclear single nucleotide polymorphisms and the discrepancy between cytoplasmic and nuclear genome lineages implied contemporary and ancient connections between neighbouring species in the form of circular overlap. This genetic connectivity was supported by principal component analysis. Strong correlations between genetic distance and geographic distance were observed, suggesting that gene flow occurs through a continuous chain around the sea. We also found that gene flow direction and intensity changed over time, with support from palaeodistribution modelling.ConclusionPast hybridization events were captured in cytoplasmic genomes, generating heterogeneity across maternal ancestries. This intensive phylogeographic study demonstrates speciation with incomplete reproductive isolation (continuous gene flow) among neighbouring species with an alteration of the direction and intensity of gene flow due to climate change. The divergence of ANSC due to repeated isolation and reconnection caused by heterogenous physiological environments and climate fluctuation provides a model to solve evolutionary scenarios for reticulate evolution in Pinaceae and other plants.

Molecular epidemiology of Omicron's CH.1.1 lineage: evidence of recombination event

The genome of the Omicron (B.1.1.529) variant is highly modified and contagious because of natural selection, which has resulted in rapid genomic evolution and the emergence of several lineages. The spike (S) protein's high number of amino acid variations in the receptor binding domain (RBD) have increased its binding affinity to the host angiotensin-converting enzyme 2 (ACE2) receptor and thus the transmissibility. However, recombination is a critical event in the evolution of RNA viruses, especially the positive single-stranded RNA viruses such as the viruses of the Coronaviridae family. The recombination event is caused by co-infection, which occurs when the genomes of lineages combine within a host cell, giving rise to a lineage with new genomic characteristics. This molecular epidemiological study demonstrates that CH.1.1, a designated lineage primarily found in Europe, is a possible recombinant derived from the likely parental BM.4 and BA.2.75.2 Omicron lineages. Regarding the genomic analysis for the possible recombination, the genomes of Wuhan/WIV04, CH.1.1, BM.4, and BA.2.75.2 were uploaded as FASTA file from the GISAID platform and aligned by UGEN software. Dissimilarity percentages were counted based on the Distance Matrix Hamming dissimilarity count. Through the genomic analysis, we identified significant amino acid variations from the hypothetical parental lineages and CH.1.1. However, the hypothesis of the possibility that CH.1.1 arises from these parental lineages is supported by their circulation simultaneously in Austria in July 2022. The dissimilarity percentage between CH.1.1 and the parental lineages, in addition to the shared inherited amino acid variations in CH.1.1.

Macroevolutionary dynamics of gene family gain and loss along multicellular eukaryotic lineages

The gain and loss of genes fluctuate over evolutionary time in major eukaryotic clades. However, the full profile of these macroevolutionary trajectories is still missing. To give a more inclusive view on the changes in genome complexity across the tree of life, here we recovered the evolutionary dynamics of gene family gain and loss ranging from the ancestor of cellular organisms to 352 eukaryotic species. We show that in all considered lineages the gene family content follows a common evolutionary pattern, where the number of gene families reaches the highest value at a major evolutionary and ecological transition, and then gradually decreases towards extant organisms. This supports theoretical predictions and suggests that the genome complexity is often decoupled from commonly perceived organismal complexity. We conclude that simplification by gene family loss is a dominant force in Phanerozoic genomes of various lineages, probably underpinned by intense ecological specializations and functional outsourcing.

Evolutionary dynamics of the accessory genomes of Staphylococcus aureus.

Staphylococcus aureus is a ubiquitous commensal and opportunistic bacterial pathogen that can cause a wide gamut of infections, which are exacerbated by the presence of multidrug-resistant and methicillin-resistant S. aureus. S. aureus is genetically heterogeneous and consists of numerous distinct lineages. Using 558 complete genomes of S. aureus, we aim to determine how the accessory genome content among phylogenetic lineages of S. aureus is structured and has evolved. Bayesian hierarchical clustering identified 10 sequence clusters, of which seven contained major sequence types (ST 1, 5, 8, 30, 59, 239, and 398). The seven sequence clusters differed in their accessory gene content, including genes associated with antimicrobial resistance and virulence. Focusing on the two largest clusters, BAPS8 and BAPS10, and each consisting mostly of ST5 and ST8, respectively, we found that the structure and connected components in the co-occurrence networks of accessory genomes varied between them. These differences are explained, in part, by the variation in the rates at which the two sequence clusters gained and lost accessory genes, with the highest rate of gene accumulation occurring recently in their evolutionary histories. We also identified a divergent group within BAPS10 that has experienced high gene gain and loss early in its history. Together, our results show highly variable and dynamic accessory genomes in S. aureus that are structured by the history of the specific lineages that carry them.IMPORTANCEStaphylococcus aureus is an opportunistic, multi-host pathogen that can cause a variety of benign and life-threatening infections. Our results revealed considerable differences in the structure and evolution of the accessory genomes of major lineages within S. aureus. Such genomic variation within a species can have important implications on disease epidemiology, pathogenesis of infection, and interactions with the vertebrate host. Our findings provide important insights into the underlying genetic basis for the success of S. aureus as a highly adaptable and resistant pathogen, which will inform current efforts to control and treat staphylococcal diseases.

Progress in the development and use of monoclonal antibodies to study the evolution and function of the immune systems in the extant lineages of ungulates

Details on the origin and function of the immune system are beginning to emerge from genomic studies tracing the origin of B and T cells and the major histocompatibility complex. This is being accomplished through identification of DNA sequences of ancestral genes present in the genomes of lineages of vertebrates that have evolved from a common primordial ancestor. Information on the evolution of the composition and function of the immune system is being obtained through development of monoclonal antibodies (mAbs) specific for the MHC class I and II molecules and differentially expressed on leukocytes differentiation molecules (LDM). The mAbs have provided the tools needed to compare the similarities and differences in the phenotype and function of immune systems that have evolved during speciation. The majority of information currently available on evolution of the composition and function of the immune system is derived from study of the immune systems in humans and mice. As described in the present review, further information is beginning to emerge from comparative studies of the immune systems in the extant lineages of species present in the two orders of ungulates, Perissodactyla and Artiodactyla. Methods have been developed to facilitate comparative research across species on pathogens affecting animal and human health.

ViTAL: Vision TrAnsformer based Low coverage SARS-CoV-2 lineage assignment.

Rapid spread of viral diseases such as Coronavirus disease 2019 (COVID-19) highlights an urgent need for efficient surveillance of virus mutation and transmission dynamics, which requires fast, inexpensive and accurate viral lineage assignment. The first two goals might be achieved through low-coverage whole-genome sequencing (LC-WGS) which enables rapid genome sequencing at scale and at reduced costs. Unfortunately, LC-WGS significantly diminishes the genomic details, rendering accurate lineage assignment very challenging. We present ViTAL, a novel deep learning algorithm specifically designed to perform lineage assignment of low coverage-sequenced genomes. ViTAL utilizes a combination of MinHash for genomic feature extraction and Vision Transformer for fine-grain genome classification and lineage assignment. We show that ViTAL outperforms state-of-the-art tools across diverse coverage levels, reaching up to 87.7% lineage assignment accuracy at 1× coverage where state-of-the-art tools such as UShER and Kraken2 achieve the accuracy of 5.4% and 27.4% respectively. ViTAL achieves comparable accuracy results with up to 8× lower coverage than state-of-the-art tools. We explore ViTAL's ability to identify the lineages of novel genomes, i.e. genomes the Vision Transformer was not trained on. We show how ViTAL can be applied to preliminary phylogenetic placement of novel variants. The data underlying this article are available in https://github.com/zuherJahshan/vital and can be accessed with 10.5281/zenodo.10688110.

ToxCodAn-Genome: an automated pipeline for toxin-gene annotation in genome assembly of venomous lineages.

The rapid development of sequencing technologies resulted in a wide expansion of genomics studies using venomous lineages. This facilitated research focusing on understanding the evolution of adaptive traits and the search for novel compounds that can be applied in agriculture and medicine. However, the toxin annotation of genomes is a laborious and time-consuming task, and no consensus pipeline is currently available. No computational tool currently exists to address the challenges specific to toxin annotation and to ensure the reproducibility of the process. Here, we present ToxCodAn-Genome, the first software designed to perform automated toxin annotation in genomes of venomous lineages. This pipeline was designed to retrieve the full-length coding sequences of toxins and to allow the detection of novel truncated paralogs and pseudogenes. We tested ToxCodAn-Genome using 12 genomes of venomous lineages and achieved high performance on recovering their current toxin annotations. This tool can be easily customized to allow improvements in the final toxin annotation set and can be expanded to virtually any venomous lineage. ToxCodAn-Genome is fast, allowing it to run on any personal computer, but it can also be executed in multicore mode, taking advantage of large high-performance servers. In addition, we provide a guide to direct future research in the venomics field to ensure a confident toxin annotation in the genome being studied. As a case study, we sequenced and annotated the toxin repertoire of Bothrops alternatus, which may facilitate future evolutionary and biomedical studies using vipers as models. ToxCodAn-Genome is suitable to perform toxin annotation in the genome of venomous species and may help to improve the reproducibility of further studies. ToxCodAn-Genome and the guide are freely available at https://github.com/pedronachtigall/ToxCodAn-Genome.

Genomic landscape of introgression from the ghost lineage in a gobiid fish uncovers the generality of forces shaping hybrid genomes.

Extinct lineages can leave legacies in the genomes of extant lineages through ancient introgressive hybridization. The patterns of genomic survival of these extinct lineages provide insight into the role of extinct lineages in current biodiversity. However, our understanding on the genomic landscape of introgression from extinct lineages remains limited due to challenges associated with locating the traces of unsampled 'ghost' extinct lineages without ancient genomes. Herein, we conducted population genomic analyses on the East China Sea (ECS) lineage of Chaenogobius annularis, which was suspected to have originated from ghost introgression, with the aim of elucidating its genomic origins and characterizing its landscape of introgression. By combining phylogeographic analysis and demographic modelling, we demonstrated that the ECS lineage originated from ancient hybridization with an extinct ghost lineage. Forward simulations based on the estimated demography indicated that the statistic γ of the HyDe analysis can be used to distinguish the differences in local introgression rates in our data. Consistent with introgression between extant organisms, we found reduced introgression from extinct lineage in regions with low recombination rates and with functional importance, thereby suggesting a role of linked selection that has eliminated the extinct lineage in shaping the hybrid genome. Moreover, we identified enrichment of repetitive elements in regions associated with ghost introgression, which was hitherto little known but was also observed in the re-analysis of published data on introgression between extant organisms. Overall, our findings underscore the unexpected similarities in the characteristics of introgression landscapes across different taxa, even in cases of ghost introgression.

Ct Value from RT-qPCR Can Predict SARS-CoV-2 Virus Assembly and Lineage Assignment Success

During the recent pandemics of COVID-19, sequencing technics became a powerful tool for gaining information about the SARS-CoV-2 virus and using this knowledge to our advantage. Thanks to this advantage, scientists all over the world were able to search for emerging variations, watching the virus evolve in real time. Assembly of the virus genomes is a crucial part of obtaining this kind of useful information. In our study, we sequenced 79 samples from nasopharyngeal swabs of COVID-19 patients. Positivity to disease was evaluated using RT-qPCR. In this work, we described the relationship between RT-qPCR Ct value and genome construction success (plus genome lineage assignment). The specific value of this study is that this relationship was described for data from metatranscriptomic sequencing of human tissue, while there was no step of viral RNA isolation (usual for genome assembly procedure). RT-qPCR Ct value and assembly quality metric NG50 were correlated. We observed that the RT-qPCR Ct value threshold of the certain success of genome assembly (Ct value < 25) and certain failure (Ct value > 30) could be drawn, while results varied for values between (with completed, completed with lower quality, and failed assemblies).

Phenotypic and genotypic characterization of Neisseria gonorrhoeae isolates from Yaoundé, Cameroon, 2019 to 2020.

This study investigated antimicrobial resistance (AMR) phenotypes and genotypes exhibited by Neisseria gonorrhoeae from Yaoundé, Cameroon. AMR to tetracycline, penicillin and ciprofloxacin was observed although none of the isolates had reduced susceptibility to azithromycin, cefixime or ceftriaxone. Whole genome sequence (WGS) data were obtained and, using a threshold of 300 or fewer locus differences in the N. gonorrhoeae core gene multilocus sequence typing (cgMLST) scheme, four distinct core genome lineages were identified. Publicly available WGS data from 1355 gonococci belonging to these four lineages were retrieved from the PubMLST database, allowing the Cameroonian isolates to be examined in the context of existing data and compared with related gonococci. Examination of AMR genotypes in this dataset found an association between the core genome and AMR with, for example, isolates belonging to the core genome group, Ng_cgc_300 : 21, possessing GyrA and ParC alleles with amino acid substitutions conferring high-level resistance to ciprofloxacin while lineages Ng_cgc_300 : 41 and Ng_cgc_300 : 243 were predicted to be susceptible to several antimicrobials. A core genome lineage, Ng_cgc_300 : 498, was observed which largely consisted of gonococci originating from Africa. Analyses from this study demonstrate the advantages of using the N. gonorrhoeae cgMLST scheme to find related gonococci to carry out genomic analyses that enhance our understanding of the population biology of this important pathogen.

No Transcriptional Compensation for Extreme Gene Dosage Imbalance in Fragmented Bacterial Endosymbionts of Cicadas.

Bacteria that form long-term intracellular associations with host cells lose many genes, a process that often results in tiny, gene-dense, and stable genomes. Paradoxically, the some of the same evolutionary processes that drive genome reduction and simplification may also cause genome expansion and complexification. A bacterial endosymbiont of cicadas, Hodgkinia cicadicola, exemplifies this paradox. In many cicada species, a single Hodgkinia lineage with a tiny, gene-dense genome has split into several interdependent cell and genome lineages. Each new Hodgkinia lineage encodes a unique subset of the ancestral unsplit genome in a complementary way, such that the collective gene contents of all lineages match the total found in the ancestral single genome. This splitting creates genetically distinct Hodgkinia cells that must function together to carry out basic cellular processes. It also creates a gene dosage problem where some genes are encoded by only a small fraction of cells while others are much more abundant. Here, by sequencing DNA and RNA of Hodgkinia from different cicada species with different amounts of splitting-along with its structurally stable, unsplit partner endosymbiont Sulcia muelleri-we show that Hodgkinia does not transcriptionally compensate to rescue the wildly unbalanced gene and genome ratios that result from lineage splitting. We also find that Hodgkinia has a reduced capacity for basic transcriptional control independent of the splitting process. Our findings reveal another layer of degeneration further pushing the limits of canonical molecular and cell biology in Hodgkinia and may partially explain its propensity to go extinct through symbiont replacement.

Sequencing of monkeypox virus from infected patients reveals viral genomes with APOBEC3-like editing, gene inactivation, and bacterial agents of skin superinfection.

A large outbreak of Monkeypox virus (MPXV) infections has arisen in May 2022 in nonendemic countries. Here, we performed DNA metagenomics using next-generation sequencing with Illumina or Nanopore technologies for clinical samples from MPXV-infected patients diagnosed between June and July 2022. Classification of the MPXV genomes and determination of their mutational patterns were performed using Nextclade. Twenty-five samples from 25 patients were studied. A MPXV genome was obtained for 18 patients, essentially from skin lesions and rectal swabbing. All 18 genomes were classified in clade IIb, lineage B.1, and we identified four B.1 sublineages (B.1.1, B.1.10, B.1.12, B.1.14). We detected a high number of mutations (range, 64-73) relatively to a 2018 Nigerian genome (genome GenBank Accession no. NC_063383.1), which were harbored by a large part of a set of 3184 MPXV genomes of lineage B.1 recovered from GenBank and Nextstrain; and we detected 35 mutations relatively to genome ON563414.3 (a B.1 lineage reference genome). Nonsynonymous mutations occurred in genes encoding central proteins, among which transcription factors and core and envelope proteins, and included two mutations that would truncate a RNA polymerase subunit and a phospholipase d-like protein, suggesting an alternative start codon and gene inactivation, respectively. A large majority (94%) of nucleotide substitutions were G > A or C > U, suggesting the action of human APOBEC3 enzymes. Finally, >1000 reads were identified as from Staphylococcus aureus and Streptococcus pyogenes for 3 and 6 samples, respectively. These findings warrant a close genomic monitoring of MPXV to get a better picture of thegenetic micro-evolution and mutational patterns of this virus, and a close clinical monitoring of skin bacterial superinfection in monkeypox patients.

Lineage skewing and genome instability underlie marrow failure in a zebrafish model of GATA2 deficiency

Inherited bone marrow failure associated with heterozygous mutations in GATA2 predisposes toward hematological malignancies, but the mechanisms remain poorly understood. Here, we investigate the mechanistic basis of marrow failure in a zebrafish model of GATA2 deficiency. Single-cell transcriptomics and chromatin accessibility assays reveal that loss of gata2a leads to skewing toward the erythroid lineage at the expense of myeloid cells, associated with loss of cebpa expression and decreased PU.1 and CEBPA transcription factor accessibility in hematopoietic stem and progenitor cells (HSPCs). Furthermore, gata2a mutants show impaired expression of npm1a, the zebrafish NPM1 ortholog. Progressive loss of npm1a in HSPCs is associated with elevated levels of DNA damage in gata2a mutants. Thus, Gata2a maintains myeloid lineage priming through cebpa and protects against genome instability and marrow failure by maintaining expression of npm1a. Our results establish a potential mechanism underlying bone marrow failure in GATA2 deficiency.

Using Chaos-Game-Representation for Analysing the SARS-CoV-2 Lineages, Newly Emerging Strains and Recombinants.

Viruses have high mutation rates, facilitating rapid evolution and the emergence of new species, subspecies, strains and recombinant forms. Accurate classification of these forms is crucial for understanding viral evolution and developing therapeutic applications. Phylogenetic classification is typically performed by analyzing molecular differences at the genomic and sub-genomic levels. This involves aligning homologous proteins or genes. However, there is growing interest in developing alignment-free methods for whole-genome comparisons that are computationally efficient. Here we elaborate on the Chaos Game Representation (CGR) method, based on concepts of statistical physics and free of sequence alignment assumptions. We adopt the CGR method for classification of the closely related clades/lineages A and B of the SARS-Corona virus 2019 (SARS-CoV-2), which is one of the fastest evolving viruses. Our study shows that the CGR approach can easily yield the SARS-CoV-2 phylogeny from the available whole genomes of lineage A and lineage B sequences. It also shows an accurate classification of eight different strains and the newly evolved XBB variant from its parental strains. Compared to alignment-based methods (Neighbour-Joining and Maximum Likelihood), the CGR method requires low computational resources, is fast and accurate for long sequences, and, being a K-mer based approach, allows simultaneous comparison of a large number of closely-related sequences of different sizes. Further, we developed an R pipeline CGRphylo, available on GitHub, which integrates the CGR module with various other R packages to create phylogenetic trees and visualize them. Our findings demonstrate the efficacy of the CGR method for accurate classification and tracking of rapidly evolving viruses, offering valuable insights into the evolution and emergence of new SARS-CoV-2 strains and recombinants.

Two Complete Genomes of Male-Killing Wolbachia Infecting Ostrinia Moth Species Illuminate Their Evolutionary Dynamics and Association with Hosts.

Wolbachia is an extremely widespread intracellular symbiont which causes reproductive manipulation on various arthropod hosts. Male progenies are killed in Wolbachia-infected lineages of the Japanese Ostrinia moth population. While the mechanism of male killing and the evolutionary interaction between host and symbiont are significant concerns for this system, the absence of Wolbachia genomic information has limited approaches to these issues. We determined the complete genome sequences of wFur and wSca, the male-killing Wolbachia of Ostrinia furnacalis and Ostrinia scapulalis. The two genomes shared an extremely high degree of homology, with over 95% of the predicted protein sequences being identical. A comparison of these two genomes revealed nearly minimal genome evolution, with a strong emphasis on the frequent genome rearrangements and the rapid evolution of ankyrin repeat-containing proteins. Additionally, we determined the mitochondrial genomes of both species' infected lineages and performed phylogenetic analyses to deduce the evolutionary dynamics of Wolbachia infection in the Ostrinia clade. According to the inferred phylogenetic relationship, two possible scenarios were proposed: (1) Wolbachia infection was established in the Ostrinia clade prior to the speciation of related species such as O. furnacalis and O. scapulalis, or (2) Wolbachia infection in these species was introgressively transferred from a currently unidentified relative. Simultaneously, the relatively high homology of mitochondrial genomes suggested recent Wolbachia introgression between infected Ostrinia species. The findings of this study collectively shed light on the host-symbiont interaction from an evolutionary standpoint.

Ferrous iron uptake via IRT1/ZIP evolved at least twice in green plants.

Iron (Fe) is essential for virtually all organisms, being irreplaceable because of its electrochemical properties that enable many biochemical processes, including photosynthesis. Besides its abundance, Fe is generally found in the poorly soluble form of ferric iron (Fe3+ ), while most plants uptake the soluble form Fe2+ . The model angiosperm Arabidopsis thaliana, for example, captures Fe through a mechanism that lowers rhizosphere pH through proton pumping that increases Fe3+ solubility, which is then reduced by a membrane-bound reductase and transported into the cell by the zinc-regulated, iron-regulated transporter-like protein (ZIP) family protein AtIRT1. ZIP proteins are transmembrane transporters of divalent metals such as Fe2+ , Zn2+ , Mn2+ , and Cd2+ . In this work, we investigated the evolution of functional homologs of IRON-REGULATED TRANSPORTER 1/ZIP in the supergroup Archaeplastida (Viridiplantae + Rhodophyta + Glaucophyta) using 51 genomes of diverse lineages. Our analyses suggest that Fe is acquired through deeply divergent ZIP proteins in land plants and chlorophyte green algae, indicating that Fe2+ uptake by ZIP proteins evolved independently at least twice throughout green plant evolution. Our results indicate that the archetypical IRON-REGULATED TRANSPORTER (IRT) proteins from angiosperms likely emerged before the origin of land plants during early streptophyte algae terrestrialization, a process that required the evolution of Fe acquisition in terrestrial subaerial settings.