Ribose-Map: a bioinformatics toolkit to map ribonucleotides embedded in genomic DNA

Host specificity and cophylogeny in the "animal-gut bacteria-phage" tripartite system

Cophylogeny has been identified between gut bacteria and their animal host and is highly relevant to host health, but little research has extended to gut bacteriophages. Here we use bee model to investigate host specificity and cophylogeny in the "animal-gut bacteria-phage" tripartite system. Through metagenomic sequencing upon different bee species, the gut phageome revealed a more variable composition than the gut bacteriome. Nevertheless, the bacteriome and the phageome showed a significant association of their dissimilarity matrices, indicating a reciprocal interaction between the two kinds of communities. Most of the gut phages were host generalist at the viral cluster level but host specialist at the viral OTU level. While the dominant gut bacteria Gilliamella and Snodgrassella exhibited matched phylogeny with bee hosts, most of their phages showed a diminished level of cophylogeny. The evolutionary rates of the bee, the gut bacteria and the gut phages showed a remarkably increasing trend, including synonymous and non-synonymous substitution and gene content variation. For all of the three codiversified tripartite members, however, their genes under positive selection and genes involving gain/loss during evolution simultaneously enriched the functions into metabolism of nutrients, therefore highlighting the tripartite coevolution that results in an enhanced ecological fitness for the whole holobiont.

Deciphering the allosteric regulation of mycobacterial inosine-5'-monophosphate dehydrogenase

Allosteric regulation of inosine 5'-monophosphate dehydrogenase (IMPDH), an essential enzyme of purine metabolism, contributes to the homeostasis of adenine and guanine nucleotides. However, the precise molecular mechanism of IMPDH regulation in bacteria remains unclear. Using biochemical and cryo-EM approaches, we reveal the intricate molecular mechanism of the IMPDH allosteric regulation in mycobacteria. The enzyme is inhibited by both GTP and (p)ppGpp, which bind to the regulatory CBS domains and, via interactions with basic residues in hinge regions, lock the catalytic core domains in a compressed conformation. This results in occlusion of inosine monophosphate (IMP) substrate binding to the active site and, ultimately, inhibition of the enzyme. The GTP and (p)ppGpp allosteric effectors bind to their dedicated sites but stabilize the compressed octamer by a common mechanism. Inhibition is relieved by the competitive displacement of GTP or (p)ppGpp by ATP allowing IMP-induced enzyme expansion. The structural knowledge and mechanistic understanding presented here open up new possibilities for the development of allosteric inhibitors with antibacterial potential.

Unraveling the complexity of asymmetric DNA replication: Advancements in ribonucleotide mapping techniques and beyond

DNA replication is a fundamental process for cell proliferation, governed by intricate mechanisms involving leading and lagging strand synthesis. In eukaryotes, canonical DNA replication occurs during the S phase of the cell cycle, facilitated by various components of the replicative machinery at sites known as replication origins. Leading and lagging strands exhibit distinct replication dynamics, with leading strand replication being relatively straightforward compared to the complex synthesis of lagging strands involving Okazaki fragment maturation. Central to DNA synthesis are DNA polymerases, with Polα, Polε, and Polδ playing pivotal roles, each specializing in specific tasks during replication. Notably, leading and lagging strands are replicated by different polymerases, contributing to the division of labor in DNA replication. Understanding the enzymology of asymmetric DNA replication has been challenging, with methods relying on ribonucleotide incorporation and next-generation sequencing techniques offering comprehensive insights. These methodologies, such as HydEn-seq, PU-seq, ribose-seq, and emRiboSeq, offer insights into polymerase activity and strand synthesis, aiding in understanding DNA replication dynamics. Recent advancements include novel conditional mutants for ribonucleotide excision repair, enzymatic cleavage alternatives, and unified pipelines for data analysis. Further developments in adapting techniques to different organisms, studying non-canonical polymerases, and exploring new sequencing platforms hold promise for expanding our understanding of DNA replication dynamics. Integrating strand-specific information into single-cell studies could offer novel insights into enzymology, opening avenues for future research and applications in repair and replication biology.

Neuronal and non-neuronal functions of the synaptic cell adhesion molecule neurexin in Nematostella vectensis

The evolutionary transition from diffusion-mediated cell-cell communication to faster, targeted synaptic signaling in animal nervous systems is still unclear. Genome sequencing analyses have revealed a widespread distribution of synapse-related genes among early-diverging metazoans, but how synaptic machinery evolved remains largely unknown. Here, we examine the function of neurexins (Nrxns), a family of presynaptic cell adhesion molecules with critical roles in bilaterian chemical synapses, using the cnidarian model, Nematostella vectensis. Delta-Nrxns are expressed mainly in neuronal cell clusters that exhibit both peptidergic and classical neurotransmitter signaling. Knockdown of δ-Nrxn reduces spontaneous peristalsis of N. vectensis polyps. Interestingly, gene knockdown and pharmacological studies suggest that δ-Nrxn is involved in glutamate- and glycine-mediated signaling rather than peptidergic signaling. Knockdown of the epithelial α-Nrxn reveals a major role in cell adhesion between ectodermal and endodermal epithelia. Overall, this study provides molecular, functional, and cellular insights into the pre-neural function of Nrxns, as well as key information for understanding how and why they were recruited to the synaptic machinery.

A comprehensive molecular phylogeny of Cephalotrichum and Microascus provides novel insights into their systematics and evolutionary history

The genera Cephalotrichum and Microascus contain ecologically, morphologically and lifestyle diverse fungi in Microascaceae (Microascales, Sordariomycetes) with a world-wide distribution. Despite previous studies having elucidated that Cephalotrichum and Microascus are highly polyphyletic, the DNA phylogeny of many traditionally morphology-defined species is still poorly resolved, and a comprehensive taxonomic overview of the two genera is lacking. To resolve this issue, we integrate broad taxon sampling strategies and the most comprehensive multi-gene (ITS, LSU, tef1 and tub2) datasets to date, with fossil calibrations to address the phylogenetic relationships and divergence times among major lineages of Microascaceae. Two previously recognised main clades, Cephalotrichum (24 species) and Microascus (49 species), were re-affirmed based on our phylogenetic analyses, as well as the phylogenetic position of 15 genera within Microascaceae. In this study, we provide an up-to-date overview on the taxonomy and phylogeny of species belonging to Cephalotrichum and Microascus, as well as detailed descriptions and illustrations of 21 species of which eight are newly described. Furthermore, the divergence time estimates indicate that the crown age of Microascaceae was around 210.37 Mya (95 % HPD: 177.18-246.96 Mya) in the Late Triassic, and that Cephalotrichum and Microascus began to diversify approximately 27.07 Mya (95 % HPD: 20.47-34.37 Mya) and 70.46 Mya (95 % HPD: 56.96-86.24 Mya), respectively. Our results also demonstrate that multigene sequence data coupled with broad taxon sampling can help elucidate previously unresolved clade relationships. Citation: Wei TP, Wu YM, Zhang X, et al. 2024. A comprehensive molecular phylogeny of Cephalotrichum and Microascus provides novel insights into their systematics and evolutionary history. Persoonia 52: 119-160. https://doi.org/10.3767/persoonia.2024.52.05 .

Protein semisynthesis reveals plasticity in HECT E3 ubiquitin ligase mechanisms

Lys ubiquitination is catalysed by E3 ubiquitin ligases and is central to the regulation of protein stability and cell signalling in normal and disease states. There are gaps in our understanding of E3 mechanisms, and here we use protein semisynthesis, chemical rescue, microscale thermophoresis and other biochemical approaches to dissect the role of catalytic base/acid function and conformational interconversion in HECT-domain E3 catalysis. We demonstrate that there is plasticity in the use of the terminal side chain or backbone carboxylate for proton transfer in HECT E3 ubiquitin ligase reactions, with yeast Rsp5 orthologues appearing to be possible evolutionary intermediates. We also show that the HECT-domain ubiquitin covalent intermediate appears to eject the E2 conjugating enzyme, promoting catalytic turnover. These findings provide key mechanistic insights into how protein ubiquitination occurs and provide a framework for understanding E3 functions and regulation.

Polymorphism-Aware Models in RevBayes: Species Trees, Disentangling Balancing Selection, and GC-Biased Gene Conversion

The role of balancing selection is a long-standing evolutionary puzzle. Balancing selection is a crucial evolutionary process that maintains genetic variation (polymorphism) over extended periods of time; however, detecting it poses a significant challenge. Building upon the Polymorphism-aware phylogenetic Models (PoMos) framework rooted in the Moran model, we introduce a PoMoBalance model. This novel approach is designed to disentangle the interplay of mutation, genetic drift, and directional selection (GC-biased gene conversion), along with the previously unexplored balancing selection pressures on ultra-long timescales comparable with species divergence times by analyzing multi-individual genomic and phylogenetic divergence data. Implemented in the open-source RevBayes Bayesian framework, PoMoBalance offers a versatile tool for inferring phylogenetic trees as well as quantifying various selective pressures. The novel aspect of our approach in studying balancing selection lies in polymorphism-aware phylogenetic models' ability to account for ancestral polymorphisms and incorporate parameters that measure frequency-dependent selection, allowing us to determine the strength of the effect and exact frequencies under selection. We implemented validation tests and assessed the model on the data simulated with SLiM and a custom Moran model simulator. Real sequence analysis of Drosophila populations reveals insights into the evolutionary dynamics of regions subject to frequency-dependent balancing selection, particularly in the context of sex-limited color dimorphism in Drosophila erecta.

Distinct features of ribonucleotides within genomic DNA in Aicardi-Goutières syndrome ortholog mutants of Saccharomyces cerevisiae

Ribonucleoside monophosphates (rNMPs) are abundantly found within genomic DNA of cells. The embedded rNMPs alter DNA properties and impact genome stability. Mutations in ribonuclease (RNase) H2, a key enzyme for rNMP removal, are associated with the Aicardi-Goutières syndrome (AGS), a severe neurological disorder. Here, we engineered orthologs of the human RNASEH2A-G37S and RNASEH2C-R69W AGS mutations in yeast Saccharomyces cerevisiae: rnh201-G42S and rnh203-K46W. Using the ribose-seq technique and the Ribose-Map bioinformatics toolkit, we unveiled rNMP abundance, composition, hotspots, and sequence context in these AGS-ortholog mutants. We found a high rNMP presence in the nuclear genome of rnh201-G42S-mutant cells, and an elevated rCMP content in both mutants, reflecting preferential cleavage of RNase H2 at rGMP. We discovered unique rNMP patterns in each mutant, showing differential activity of the AGS mutants on the leading or lagging replication strands. This study guides future research on rNMP characteristics in human genomes with AGS mutations.

Comparative chloroplast genomics, phylogenetic relationships and molecular markers development of Aglaonema commutatum and seven green cultivars of Aglaonema

Aglaonema commutatum is a famous species in the Aglaonema genus, which has important ornamental and economic value. However, its chloroplast genome information and phylogenetic relationships among popular green cultivars of Aglaonema in southern China have not been reported. Herein, chloroplast genomes of one variety of A. commutatum and seven green cultivars of Aglaonema, namely, A. commutatum 'San Remo', 'Kai Sa', 'Pattaya Beauty', 'Sapphire', 'Silver Queen', 'Snow White', 'White Gem', and 'White Horse Prince', were sequenced and assembled for comparative analysis and phylogeny. These eight genomes possessed a typical quadripartite structure that consisted of a LSC region (90,799-91,486 bp), an SSC region (20,508-21,137 bp) and a pair of IR regions (26,661-26,750 bp). Each genome contained 112 different genes, comprising 79 protein-coding genes, 29 tRNA genes and 4 rRNA genes. The gene orders, GC contents, codon usage frequency, and IR/SC boundaries were highly conserved among these eight genomes. Long repeats, SSRs, SNPs and indels were analyzed among these eight genomes. Comparative analysis of 15 Aglaonema chloroplast genomes identified 7 highly variable regions, including trnH-GUG-exon1-psbA, trnS-GCU-trnG-UCC-exon1, trnY-GUA-trnE-UUC, psbC-trnS-UGA, trnF-GAA-ndhJ, ccsA-ndhD, and rps15-ycf1-D2. Reconstruction of the phylogenetic trees based on chloroplast genomes, strongly supported that Aglaonema was a sister to Anchomanes, and that the Aglaonema genus was classified into two sister clades including clade I and clade II, which corresponded to two sections, Aglaonema and Chamaecaulon, respectively. One variety and five cultivars, including A. commutatum 'San Remo', 'Kai Sa', 'Pattaya Beauty', 'Silver Queen', 'Snow White', and 'White Horse Prince', were classified into clade I; and the rest of the two cultivars, including 'Sapphire' and 'White Gem', were classified into clade II. Positive selection was observed in 34 protein-coding genes at the level of the amino acid sites among 77 chloroplast genomes of the Araceae family. Based on the highly variable regions and SSRs, 4 DNA markers were developed to differentiate the clade I and clade II in Aglaonema. In conclusion, this study provided chloroplast genomic resources for Aglaonema, which were useful for its classification and phylogeny.

Identification of the mutual gliding locus as a factor for gut colonization in non-native bee hosts using the ARTP mutagenesis

BACKGROUND

The gut microbiota and their hosts profoundly affect each other's physiology and evolution. Identifying host-selected traits is crucial to understanding the processes that govern the evolving interactions between animals and symbiotic microbes. Current experimental approaches mainly focus on the model bacteria, like hypermutating Escherichia coli or the evolutionary changes of wild stains by host transmissions. A method called atmospheric and room temperature plasma (ARTP) may overcome the bottleneck of low spontaneous mutation rates while maintaining mild conditions for the gut bacteria.

RESULTS

We established an experimental symbiotic system with gnotobiotic bee models to unravel the molecular mechanisms promoting host colonization. By in vivo serial passage, we tracked the genetic changes of ARTP-treated Snodgrassella strains from Bombus terrestris in the non-native honeybee host. We observed that passaged isolates showing genetic changes in the mutual gliding locus have a competitive advantage in the non-native host. Specifically, alleles in the orphan mglB, the GTPase activating protein, promoted colonization potentially by altering the type IV pili-dependent motility of the cells. Finally, competition assays confirmed that the mutations out-competed the ancestral strain in the non-native honeybee gut but not in the native host.

CONCLUSIONS

Using the ARTP mutagenesis to generate a mutation library of gut symbionts, we explored the potential genetic mechanisms for improved gut colonization in non-native hosts. Our findings demonstrate the implication of the cell mutual-gliding motility in host association and provide an experimental system for future study on host-microbe interactions. Video Abstract.

rNMPID: a database for riboNucleoside MonoPhosphates in DNA

Motivation

Ribonucleoside monophosphates (rNMPs) are the most abundant non-standard nucleotides embedded in genomic DNA. If the presence of rNMP in DNA cannot be controlled, it can lead to genome instability. The actual regulatory functions of rNMPs in DNA remain mainly unknown. Considering the association between rNMP embedment and various diseases and cancer, the phenomenon of rNMP embedment in DNA has become a prominent area of research in recent years.

Results

We introduce the rNMPID database, which is the first database revealing rNMP-embedment characteristics, strand bias, and preferred incorporation patterns in the genomic DNA of samples from bacterial to human cells of different genetic backgrounds. The rNMPID database uses datasets generated by different rNMP-mapping techniques. It provides the researchers with a solid foundation to explore the features of rNMP embedded in the genomic DNA of multiple sources, and their association with cellular functions, and, in future, disease. It also significantly benefits researchers in the fields of genetics and genomics who aim to integrate their studies with the rNMP-embedment data.

Availability and implementation

rNMPID is freely accessible on the web at https://www.rnmpid.org.

Chloroplast genome analysis and evolutionary insights in the versatile medicinal plant Calendula officinalis L

Calendula officinalis L.is a versatile medicinal plant with numerous applications in various fields. However, its chloroplast genome structure, features, phylogeny, and patterns of evolution and mutation remain largely unexplored. This study examines the chloroplast genome, phylogeny, codon usage bias, and divergence time of C. officinalis, enhancing our understanding of its evolution and adaptation. The chloroplast genome of C. officinalis is a 150,465 bp circular molecule with a G + C content of 37.75% and comprises 131 genes. Phylogenetic analysis revealed a close relationship between C. officinalis, C. arvensis, and Osteospermum ecklonis. A key finding is the similarity in codon usage bias among these species, which, coupled with the divergence time analysis, supports their close phylogenetic proximity. This similarity in codon preference and divergence times underscores a parallel evolutionary adaptation journey for these species, highlighting the intricate interplay between genetic evolution and environmental adaptation in the Asteraceae family. Moreover unique evolutionary features in C. officinalis, possibly associated with certain genes were identified, laying a foundation for future research into the genetic diversity and medicinal value of C. officinalis.

Evolutionary analysis reveals the origin of sodium coupling in glutamate transporters

Secondary active membrane transporters harness the energy of ion gradients to concentrate their substrates. Homologous transporters evolved to couple transport to different ions in response to changing environments and needs. The bases of such diversification, and thus principles of ion coupling, are unexplored. Employing phylogenetics and ancestral protein reconstruction, we investigated sodium-coupled transport in prokaryotic glutamate transporters, a mechanism ubiquitous across life domains and critical to neurotransmitter recycling in humans. We found that the evolutionary transition from sodium-dependent to independent substrate binding to the transporter preceded changes in the coupling mechanism. Structural and functional experiments suggest that the transition entailed allosteric mutations, making sodium binding dispensable without affecting ion-binding sites. Allosteric tuning of transporters' energy landscapes might be a widespread route of their functional diversification.

Antioxidant compounds produced by endolichenic fungus Penicillium sp.-strain 1322P isolated from Pyxine subcinerea

Endolichenic fungi are expecting for new bioresources of pharmacological compounds. However, the number of investigations targeting antioxidant compounds produced by endolichenic fungi remains limited. To discover new antioxidant compounds, we analyzed the antioxidant activity of the methanol extracts derived from isolated lichen mycobionts or endolichenic fungi induced from Pyxine subcinerea. We performed this analysis using the oxygen radical absorbance capacity (ORAC) method. As a result, we isolated from an endolichenic fungus identified as Penicillium sp.-stain 1322P in Pyxine subcinerea. This fungus produced a red pigment, and its chemical structure was determined to be sclerotioramine based on the analytical data obtained from NMR, LC-MS/MS, and HPLC-PDA. Sclerotioramine exhibited high antioxidant activity, and the ORAC values (mean ± SD) of sclerotioramine and sclerotiorin were 11.4 ± 0.36 and 4.86 ± 0.70 mmol TE per gram of the respective pure compound. Thus, the antioxidant activity of sclerotioramine was greater than twice that of sclerotiorin. This work represents the first report that the antioxidant activity of sclerotioramine is higher than that of the sclerotiorin.

Light-strand bias and enriched zones of embedded ribonucleotides are associated with DNA replication and transcription in the human-mitochondrial genome

Abundant ribonucleoside-triphosphate (rNTP) incorporation into DNA by DNA polymerases in the form of ribonucleoside monophosphates (rNMPs) is a widespread phenomenon in nature, resulting in DNA-structural change and genome instability. The rNMP distribution, characteristics, hotspots and association with DNA metabolic processes in human mitochondrial DNA (hmtDNA) remain mostly unknown. Here, we utilize the ribose-seq technique to capture embedded rNMPs in hmtDNA of six different cell types. In most cell types, the rNMPs are preferentially embedded on the light strand of hmtDNA with a strong bias towards rCMPs; while in the liver-tissue cells, the rNMPs are predominately found on the heavy strand. We uncover common rNMP hotspots and conserved rNMP-enriched zones across the entire hmtDNA, including in the control region, which links the rNMP presence to the frequent hmtDNA replication-failure events. We show a strong correlation between coding-sequence size and rNMP-embedment frequency per nucleotide on the non-template, light strand in all cell types, supporting the presence of transient RNA-DNA hybrids preceding light-strand replication. Moreover, we detect rNMP-embedment patterns that are only partly conserved across the different cell types and are distinct from those found in yeast mtDNA. The study opens new research directions to understand the biology of hmtDNA and genomic rNMPs.

Seagrass genomes reveal ancient polyploidy and adaptations to the marine environment

We present chromosome-level genome assemblies from representative species of three independently evolved seagrass lineages: Posidonia oceanica, Cymodocea nodosa, Thalassia testudinum and Zostera marina. We also include a draft genome of Potamogeton acutifolius, belonging to a freshwater sister lineage to Zosteraceae. All seagrass species share an ancient whole-genome triplication, while additional whole-genome duplications were uncovered for C. nodosa, Z. marina and P. acutifolius. Comparative analysis of selected gene families suggests that the transition from submerged-freshwater to submerged-marine environments mainly involved fine-tuning of multiple processes (such as osmoregulation, salinity, light capture, carbon acquisition and temperature) that all had to happen in parallel, probably explaining why adaptation to a marine lifestyle has been exceedingly rare. Major gene losses related to stomata, volatiles, defence and lignification are probably a consequence of the return to the sea rather than the cause of it. These new genomes will accelerate functional studies and solutions, as continuing losses of the 'savannahs of the sea' are of major concern in times of climate change and loss of biodiversity.

The Vibrio cholerae CBASS phage defence system modulates resistance and killing by antifolate antibiotics

Toxic bacterial modules such as toxin-antitoxin systems hold antimicrobial potential, though successful applications are rare. Here we show that in Vibrio cholerae the cyclic-oligonucleotide-based anti-phage signalling system (CBASS), another example of a toxic module, increases sensitivity to antifolate antibiotics up to 10×, interferes with their synergy and ultimately enables bacterial lysis by these otherwise classic bacteriostatic antibiotics. Cyclic-oligonucleotide production by the CBASS nucleotidyltransferase DncV upon antifolate treatment confirms full CBASS activation under these conditions, and suggests that antifolates release DncV allosteric inhibition by folates. Consequently, the CBASS-antifolate interaction is specific to CBASS systems with closely related nucleotidyltransferases and similar folate-binding pockets. Last, antifolate resistance genes abolish the CBASS-antifolate interaction by bypassing the effects of on-target antifolate activity, thereby creating potential for their coevolution with CBASS. Altogether, our findings illustrate how toxic modules can impact antibiotic activity and ultimately confer bactericidal activity to classical bacteriostatic antibiotics.

Population-level variation in gut bifidobacterial composition and association with geography, age, ethnicity, and staple food

Bifidobacteria are key gut commensals that confer various health benefits and are commonly used as probiotics. However, little is known about the population-level variation in gut bifidobacterial composition and its affecting factors. Therefore, we analyzed Bifidobacterium species with amplicon sequencing of the groEL gene on fecal samples of 1674 healthy individuals, who belonged to eight ethnic groups and resided in 60 counties/cities of 28 provinces across China. We found that the composition of the bifidobacterial community was associated with geographical factors, demographic characteristics, staple food type, and urbanization. First, geography, which reflects a mixed effect of other variables, explained the largest variation in the bifidobacterial profile. Second, middle adolescence (age 14-17) and age 30 were two key change points in the bifidobacterial community development, and a bifidobacterial community resembling that of adults occurred in middle adolescence, which is much later than the maturation of the whole gut microbial community at approximately age 3. Third, each ethnicity showed a distinct bifidobacterial profile, and the remarkable amount of unknown Bifidobacterium species in the Tibetan gut suggested undiscovered biodiversity. Fourth, wheat as the main staple food promoted the flourish of B. adolescentis and B. longum. Fifth, alpha diversity of the bifidobacterial community decreased with urbanization. Collectively, our findings provide insight into the environmental and host factors that shape the human gut bifidobacterial community, which is fundamental for precision probiotics.

Systematic review and integrated data analysis reveal diverse pangolin-associated microbes with infection potential

There has been increasing global concern about the spillover transmission of pangolin-associated microbes. To assess the risk of these microbes for emergence as human pathogens, we integrated data from multiple sources to describe the distribution and spectrum of microbes harbored by pangolins. Wild and trafficked pangolins have been mainly recorded in Asia and Africa, while captive pangolins have been reported in European and North American countries. A total of 128 microbes, including 92 viruses, 25 bacteria, eight protists, and three uncharacterized microbes, have been identified in five pangolin species. Out of 128 pangolin-associated microbes, 31 (including 13 viruses, 15 bacteria, and three protists) have been reported in humans, and 54 are animal-associated viruses. The phylogenetic analysis of human-associated viruses carried by pangolins reveals that they are genetically close to those naturally circulating among human populations in the world. Pangolins harbor diverse microbes, many of which have been previously reported in humans and animals. Abundant viruses initially detected in pangolins might exhibit risks for spillover transmission.

Distinct features of ribonucleotides within genomic DNA in Aicardi-Goutières syndrome (AGS)-ortholog mutants of Saccharomyces cerevisiae

Ribonucleoside monophosphates (rNMPs) are abundantly found within genomic DNA of cells. The embedded rNMPs alter DNA properties and impact genome stability. Mutations in ribonuclease (RNase) H2, a key enzyme for rNMP removal, are associated with the Aicardi-Goutières syndrome (AGS), a severe neurological disorder. Here, we engineered two AGS-ortholog mutations in Saccharomyces cerevisiae: rnh201-G42S and rnh203-K46W. Using the ribose-seq technique and the Ribose-Map bioinformatics toolkit, we unveiled rNMP abundance, composition, hotspots, and sequence context in these yeast AGS-ortholog mutants. We found higher rNMP incorporation in the nuclear genome of rnh201-G42S than in wild-type and rnh203-K46W-mutant cells, and an elevated rCMP content in both mutants. Moreover, we uncovered unique rNMP patterns in each mutant, highlighting a differential activity of the AGS mutants towards rNMPs embedded on the leading or on the lagging strand of DNA replication. This study guides future research on rNMP characteristics in human genomic samples carrying AGS mutations.

Mapping Structural Distribution and Gating-Property Impacts of Disease-Associated Missense Mutations in Voltage-Gated Sodium Channels

Thousands of voltage-gated sodium (Nav) channel variants contribute to a variety of disorders, including epilepsy, autism, cardiac arrhythmia, and pain disorders. Yet variant effects of more mutations remain unclear. The conventional gain-of-function (GoF) or loss-of-function (LoF) classifications is frequently employed to interpret of variant effects on function and guide precision therapy for sodium channelopathies. Our study challenges this binary classification by analyzing 525 mutations associated with 34 diseases across 366 electrophysiology studies, revealing that diseases with similar phenotypic effects can stem from unique molecular mechanisms. Our results show a high biophysical agreement (86%) between homologous disease-associated variants in different Nav genes, significantly surpassing the 60% phenotype (GoFo/LoFo) agreement among homologous mutants, suggesting the need for more nuanced disease categorization and treatment based on specific gating-property changes. Using UniProt data, we mapped over 2,400 disease-associated missense variants across nine human Nav channels and identified three clusters of mutation hotspots. Our findings indicate that mutations near the selectivity filter generally diminish the maximal current amplitude, while those in the fast inactivation region lean towards a depolarizing shift in half-inactivation voltage in steady-state activation, and mutations in the activation gate commonly enhance persistent current. In contrast to mutations in the PD, those within the VSD exhibit diverse impacts and subtle preferences on channel activity. This study shows great potential to enhance prediction accuracy for variant effects based on the structural context, laying the groundwork for targeted drug design in precision medicine.

The diverse evolutionary histories of domesticated metaviral capsid genes in mammals

Selfish genetic elements and their remnants comprise at least half of the human genome. Active transposons duplicate by inserting copies at new sites in a host genome. Following insertion, transposons can acquire mutations that render them inactive; the accrual of additional mutations can render them unrecognizable over time. However, in rare instances, segments of transposons become useful for the host, in a process called gene domestication. Using the first complete human genome assembly and 25 additional vertebrate genomes, we analyzed the evolutionary trajectories and functional potential of genes domesticated from the capsid genes of Metaviridae, a retroviral-like retrotransposon family. Our analysis reveals four families of domesticated capsid genes in placental mammals with varied evolutionary outcomes, ranging from universal retention to lineage-specific duplications or losses and from purifying selection to lineage-specific rapid evolution. The four families of domesticated capsid genes have divergent amino-terminal domains, inherited from four distinct ancestral metaviruses. Structural predictions reveal that many domesticated genes encode a previously unrecognized RNA-binding domain retained in multiple paralogs in mammalian genomes both adjacent to and independent from the capsid domain. Collectively, our study reveals diverse outcomes of domestication of diverse metaviruses, which led to structurally and evolutionarily diverse genes that encode important, but still largely-unknown functions in placental mammals. (207).

Diversity of Bathyarchaeia viruses in metagenomes and virus-encoded CRISPR system components

Bathyarchaeia represent a class of archaea common and abundant in sedimentary ecosystems. The virome of Bathyarchaeia so far has not been characterized. Here we report 56 metagenome-assembled genomes of Bathyarchaeia viruses identified in metagenomes from different environments. Gene sharing network and phylogenomic analyses led to the proposal of four virus families, including viruses of the realms Duplodnaviria and Adnaviria, and archaea-specific spindle-shaped viruses. Genomic analyses uncovered diverse CRISPR elements in these viruses. Viruses of the proposed family 'Fuxiviridae' harbor an atypical type IV-B CRISPR-Cas system and a Cas4 protein that might interfere with host immunity. Viruses of the family 'Chiyouviridae' encode a Cas2-like endonuclease and two mini-CRISPR arrays, one with a repeat identical to that in the host CRISPR array, potentially allowing the virus to recruit the host CRISPR adaptation machinery to acquire spacers that could contribute to competition with other mobile genetic elements or to inhibition of host defenses. These findings present an outline of the Bathyarchaeia virome and offer a glimpse into their counter-defense mechanisms.

Comparative and phylogenetic analysis of complete chloroplast genome sequences of Salvia regarding its worldwide distribution

Salvia is widely used as medicine, food, and ornamental plants all over the world, with three main distribution centers, the Central and western Asia/Mediterranean (CAM), the East Aisa (EA), and the Central and South America (CASA). Along with its large number of species and world-wide distribution, Salvia is paraphyletic with multiple diversity. Chloroplast genomes (CPs) are useful tools for analyzing the phylogeny of plants at lower taxonomic levels. In this study, we reported chloroplast genomes of five species of Salvia and performed phylogenetic analysis with current available CPs of Salvia. Repeated sequence analysis and comparative analysis of Salvia CPs were also performed with representative species from different distribution centers. The results showed that the genetic characters of the CPs are related to the geographic distribution of plants. Species from CAM diverged first to form a separate group, followed by species from EA, and finally species from CASA. Larger variations of CPs were observed in species from CAM, whereas more deficient sequences and less repeated sequences in the CPs were observed in species from CASA. These results provide valuable information on the development and utilization of the worldwide genetic resources of Salvia.

One-Step Reverse-Transcription Loop-Mediated Isothermal Amplification (RT-LAMP) Assay for the Detection of Canna Yellow Streak Virus

Canna yellow streak virus (CaYSV) is a potyvirus that causes severe damage to the ornamental plant canna in the United Kingdom and Brazil. Here, we identified CaYSV in China by isolating total RNA from an infected plant, amplifying the virus genome segments, and cloning and sequencing the amplicons. After assembly, the full-length genome of the virus was obtained and uploaded to the NCBI database. Phylogenetic analysis results showed that the Guizhou isolate (OL546222) was most closely related to the KS isolate (MG545919.1). Virus detection is essential for virus disease control but the subclinical infection of CaYSV on canna in its early development increases the difficulty of CaYSV diagnosis. The goal of this study was to develop an efficient method for detection of CaYSV. We designed the primers, optimized the reaction conditions, and finally established a one-step reverse-transcription loop-mediated isothermal amplification (RT-LAMP) method. The product of RT-LAMP can be analyzed by both agarose gel electrophoresis and visible color change. The established one-step RT-LAMP assay showed high specificity and sensitivity in detecting CaYSV. This RT-LAMP method was also applied in analysis of 61 field samples collected from Guizhou and Jiangsu Provinces. The results showed that the infection rates of CaYSV on canna samples from these two provinces were very high (63 and 96% respectively).

Alternative Z-genome biosynthesis pathway shows evolutionary progression from Archaea to phage

Many bacteriophages evade bacterial immune recognition by substituting adenine with 2,6-diaminopurine (Z) in their genomes. The Z-genome biosynthetic pathway involves PurZ that belongs to the PurA (adenylosuccinate synthetase) family and bears particular similarity to archaeal PurA. However, how the transition of PurA to PurZ occurred during evolution is not clear; recapturing this process may shed light on the origin of Z-containing phages. Here we describe the computer-guided identification and biochemical characterization of a naturally existing PurZ variant, PurZ0, which uses guanosine triphosphate as the phosphate donor rather than the ATP used by PurZ. The atomic resolution structure of PurZ0 reveals a guanine nucleotide binding pocket highly analogous to that of archaeal PurA. Phylogenetic analyses suggest PurZ0 as an intermediate during the evolution of archaeal PurA to phage PurZ. Maintaining the balance of different purines necessitates further evolvement of guanosine triphosphate-using PurZ0 to ATP-using PurZ in adaptation to Z-genome life.

Chromosome-level genome assemblies from two sandalwood species provide insights into the evolution of the Santalales

Sandalwood is one of the most expensive woods in the world and is well known for its long-lasting and distinctive aroma. In our study, chromosome-level genome assemblies for two sandalwood species (Santalum album and Santalum yasi) were constructed by integrating NGS short reads, RNA-seq, and Hi-C libraries with PacBio HiFi long reads. The S. album and S. yasi genomes were both assembled into 10 pseudochromosomes with a length of 229.59 Mb and 232.64 Mb, containing 21,673 and 22,816 predicted genes and a repeat content of 28.93% and 29.54% of the total genomes, respectively. Further analyses resolved a Santalum-specific whole-genome triplication event after divergence from ancestors of the Santalales lineage Malania, yet due to dramatic differences in transposon content, the Santalum genomes were only one-sixth the size of the Malania oleifera genome. Examination of RNA-seq data revealed a suite of genes that are differentially expressed in haustoria and might be involved in host hemiparasite interactions. The two genomes presented here not only provide an important comparative dataset for studying genome evolution in early diverging eudicots and hemiparasitic plants but will also hasten the application of conservation genomics for a lineage of trees recovering from decades of overexploitation.

The complete chloroplast genome sequence of Sedum bulbiferum (Crassulaceae)

Sedum bulbiferum is a traditional medicinal plant in China, with few reports on its chloroplast genome. In this study, the chloroplast genome of Sedum bulbiferum was characterized, and its phylogenetic position among other closely related species was studied. The results showed that the full length of the chloroplast genome was 150,074 bp, containing a large single-copy (LSC) region and a small single-copy (SSC) region of 81,730 and 16,726 bp, respectively, as well as two inverted repeat regions (IRs) of 25,809 bp like other plants. A total of 128 genes were found, including 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Phylogenetic analysis showed that Sedum bulbiferum is closely related to Sedum emarginatum, Sedum alfredii, Sedum tricarpum, Sedum plumbizincicola, and Sedum sarmentosum.

The soil microbiomes of forest ecosystems in Kenya: their diversity and environmental drivers

Soil microbiomes in forest ecosystems act as both nutrient sources and sinks through a range of processes including organic matter decomposition, nutrient cycling, and humic compound incorporation into the soil. Most forest soil microbial diversity studies have been performed in the northern hemisphere, and very little has been done in forests within African continent. This study examined the composition, diversity and distribution of prokaryotes in Kenyan forests top soils using amplicon sequencing of V4-V5 hypervariable region of the 16S rRNA gene. Additionally, soil physicochemical characteristics were measured to identify abiotic drivers of prokaryotic distribution. Different forest soils were found to have statistically distinct microbiome compositions, with Proteobacteria and Crenarchaeota taxa being the most differentially abundant across regions within bacterial and archaeal phyla, respectively. Key bacterial community drivers included pH, Ca, K, Fe, and total N while archaeal diversity was shaped by Na, pH, Ca, total P and total N. To contextualize the prokaryote diversity of Kenyan forest soils on a global scale, the sample set was compared to amplicon data obtained from forest biomes across the globe; displaying them to harbor distinct microbiomes with an over-representation of uncultured taxa such as TK-10 and Ellin6067 genera.

Genomic insights into the coupling of a Chlorella-like microeukaryote and sulfur bacteria in the chemocline of permanently stratified Lake Cadagno

Meromictic Lake Cadagno is a permanently stratified system with a persistent microbial bloom within the oxic-anoxic boundary called the chemocline. The association between oxygenic and anoxygenic photosynthesis within the chemocline has been known for at least two decades. Although anoxygenic purple and green sulfur bacteria have been well studied, reports on oxygenic phytoplankton have remained sparse since their discovery in the 1920s. Nearly a century later, this study presents the first near-complete genome of a photosynthetic microbial eukaryote from the chemocline of Lake Cadagno, provisionally named Chlorella-like MAG. The 18.9 Mbp nuclear genome displays a high GC content (71.5%), and the phylogenetic placement suggests that it is a novel species of the genus Chlorella of Chlorophytes. Functional annotation of the Chlorella-like metagenome-assembled genome predicted 10,732 protein-coding genes, with an approximate 0.6% proportion potentially involved in carbon, sulfur, and nitrogen (C, N, and S) metabolism. In addition to C4 photosynthesis, this study detected genes for heat shock proteins (HSPs) in the Chlorella-like algae, consistent with the other Chlorella species. Altogether, the genomic insights in this study suggest the cooperation of photosynthetic algae with phototrophic sulfur bacteria via C, N, and S metabolism, which may aid their collective persistence in the Lake Cadagno chemocline. Furthermore, this work additionally presents the chloroplast genome of Cryptomonas-like species, which was likely to be presumed as cyanobacteria in previous studies because of the presence of phycobilisomes.

Amine recognizing domain in diverse receptors from bacteria and archaea evolved from the universal amino acid sensor

Bacteria contain many different receptor families that sense different signals permitting an optimal adaptation to the environment. A major limitation in microbiology is the lack of information on the signal molecules that activate receptors. Due to a significant sequence divergence, the signal recognized by sensor domains is only poorly reflected in overall sequence identity. Biogenic amines are of central physiological relevance for microorganisms and serve for example as substrates for aerobic and anaerobic growth, neurotransmitters or osmoprotectants. Based on protein structural information and sequence analysis, we report here the identification of a sequence motif that is specific for amine-sensing dCache sensor domains (dCache_1AM). These domains were identified in more than 13,000 proteins from 8,000 bacterial and archaeal species. dCache_1AM containing receptors were identified in all major receptor families including sensor kinases, chemoreceptors, receptors involved in second messenger homeostasis and Ser/Thr phosphatases. The screening of compound libraries and microcalorimetric titrations of selected dCache_1AM domains confirmed their capacity to specifically bind amines. Mutants in the amine binding motif or domains that contain a single mismatch in the binding motif, had either no or a largely reduced affinity for amines, illustrating the specificity of this motif. We demonstrate that the dCache_1AM domain has evolved from the universal amino acid sensing domain, providing novel insight into receptor evolution. Our approach enables precise "wet"-lab experiments to define the function of regulatory systems and thus holds a strong promise to address an important bottleneck in microbiology: the identification of signals that stimulate numerous receptors.

Characterization and phylogenetic analysis of the chloroplast genome of galium spurium

Galium spurium is a farmland weed, with strong stress resistance. However, its chloroplast genome has never been reported. In this study, the complete sequence of the chloroplast genome of G. spurium was characterized, which is a circular molecule, 153,481 bp in length, and with a large single copy region of 84,334 bp, a small single copy region of 17,057 bp, and a pair of inverted repeat regions of 26,045 bp. The whole genome contained 127 genes, including 82 protein-coding genes, 37 transfer RNA genes, and eight ribosomal RNA genes. Phylogenetic analysis shows that it relates closely to G. aparine. This study provides a basis for the further phylogenic study of Galium.

Sulfoquinovose is a widespread organosulfur substrate for Roseobacter clade bacteria in the ocean

Sulfoquinovose (SQ) is one of the most abundant organosulfur compounds in the biosphere, and its biosynthesis and degradation can represent an important contribution to the sulfur cycle. To data, in marine environments, the microorganisms capable of metabolising SQ have remained unidentified and the sources of SQ are still uncertain. Herein, the marine Roseobacter clade bacteria (RCB) Dinoroseobacter shibae DFL 12 and Roseobacter denitrificans OCh 114 were found to grow using SQ as the sole source of carbon and energy. In the presence of SQ, we identified a set of highly up-regulated proteins encoded by gene clusters in these two organisms, of which four homologues to proteins in the SQ monooxygenase pathway of Agrobacterium fabrum C58 may confer the ability to metabolise SQ to these marine bacteria. The sulfite released from SQ desulfonation by FMN-dependent SQ monooxygenase (SmoC) may provide bacteria with reduced sulfur for assimilation, while proteins associated with sulfite production via assimilatory sulfate reduction were significantly down-regulated. Such SQ catabolic genes are restricted to a limited number of phylogenetically diverse bacterial taxa with the predominate genera belonging to the Roseobacter clade (Roseobacteraceae). Moreover, transcript analysis of Tara Oceans project and coastal Bohai Sea samples provided additional evidence for SQ metabolism by RCB. SQ was found to be widely distributed in marine phytoplankton and cyanobacteria with variable intracellular concentrations ranging from micromolar to millimolar levels, and the amounts of SQ on particulate organic matter in field samples were, on average, lower than that of dimethylsulfoniopropionate (DMSP) by one order of magnitude. Together, the phototroph-derived SQ actively metabolised by RCB represents a previously unidentified link in the marine sulfur cycle.

The complete chloroplast genome sequence of Centaurea cyanus (Asteraceae)

Centaurea cyanus has been a weed in farmland for a long time. In this study, the chloroplast genome of C. cyanus was sequenced to establish the phylogenetic relationship between its genomic characteristics and other related species. The chloroplast gene structure of C. cyanus is a circular molecule with a length of 152,433 bp, including a large single-copy (LSC) region of 83,464 bp, a small single-copy (SSC) region of 18,545 bp, and a pair of inverted repeats sequences (IRs) region of 25,212 bp. The whole genome contains 130 genes, including 86 protein-coding genes, 36 tRNA genes, and eight rRNA genes. Phylogenetic analysis showed that C. cyanus is close to Carthamus. tinctorius, C. tinctorius, C. diffusa, and C. maculosa, and all of them were in one clade. This study provides genetic resource information for the further study of Centaurea.

The complete chloroplast genome of Striga asiatica (L.) Kuntze 1891 (Orobanchaceae), a hemiparasitic weed from Guangxi China

Striga asiatica (L.) Kuntze 1891 is a hemiparasitic plant native to Asia and Africa. It is invasive and causes yield losses in crops such as corn, rice and sorghum. Lack of chloroplast genomic data has limited research into its obligate parasitic lifestyle. In this study, the complete chloroplast genome of Striga asiatica was sequenced and characterized. It is a quadripartite structure with a total length of 191,085 bp and a GC content of 37.86%. It has a large single copy region (LSC) of 51,406 bp, a small single copy region (SSC) of 273 bp, and two copies of the reverse repeat sequence (IRA and IRB) of 69,703 bp. A total of 122 protein-coding genes, 8 rRNA genes, and 44 tRNA genes were annotated in the chloroplast genome. There were a lot of ndh gene deletions and pseudogenizations in this chloroplast genome. For example, ndhA, D, E, H, I, and K were all pseudogenes because they were missing the 5' end start codon. ndhB, C, and J had shorter gene lengths than their homologs, and ndhF and ndhG were missing genes. The phylogenetic tree reveals that all Striga species form a clade, and a bootstrap value of 100 indicates that S. asiatica is closely related to Striga hermonthica and Striga sepera. The comprehensive chloroplast genomic resource of S. asiatica would assist researchers in comprehending hemiparasitic mechanisms, molecular markers, and evolutionary patterns of the genus Striga.

The complete plastome sequence of Momordica cochinchinensis (Cucurbitaceae)

Momordica cochinchinensis (Lour.) Spreng. is an important medicinal plant that is used to treat various diseases in South and Southeast Asia. In this study, the complete plastome of M. cochinchinensis was sequenced and found to exhibit a total length of 158,955 bp, with a large single copy (LSC) region of 87,924 bp and a small single copy (SSC) region of 18,479 bp, as well as with two inverted repeats (IRs) that were both 26,726 bp in length. In total, 129 genes were detected, comprising 86 protein-encoding genes, 8 ribosomal RNA (rRNA) genes, and 35 transfer RNA (tRNA) genes. Furthermore, the inferred phylogenetic tree confirmed that M. cochinchinensis belongs to the genus Momordica in the Cucurbitaceae family. The research results will be used for authenticating M. cochinchinensis plant materials and for analyzing the genetic diversity and phylogenetic relationships in Momordica.

Mutational profile confers increased stability of SARS-CoV-2 spike protein in Brazilian isolates

Spike (S) protein has been recognized as a promising molecular target for diagnostic, vaccines and antiviral drugs development for COVID-19. In this study, we analyzed the most predominant mutations in the S protein of Brazilian isolates and predicted the effect of these amino acid alterations to protein conformation. A total of 25,924 sequences were obtained from GISAID for five regions of Brazilian territory (Midwest, North, Northeast, South, and Southeast), according to exclusion criteria. Most of the SARS-CoV-2 isolates belongs to the G clade and showed a large occurrence of D614G, N501Y and L18F substitutions. Prediction effects of these amino acid substitutions on the structure dynamics of the spike protein indicated a positive ΔΔG values and negative ΔΔS_Vib in most cases which is associated to structural stabilization and flexibility reduction of the S protein. Mutations E484K, N501Y and K417N belong to several SARS-CoV-2 variants of concern such as Alpha, Beta, Gamma and Delta, and showed high incidence among Brazilian isolates. These mutations have been described to increase RBD affinity to ACE-2 host and abolishment of RBD affinity to potent neutralizing ant-RBD. The increase in rates of infection and reinfection requires continuous genomic surveillance studies in order to characterize emerging mutations and monitor vaccine efficacy, and thus consideration structural data and dynamics in the observed phenotypes.Communicated by Ramaswamy H. Sarma.

Variation analysis of SARS-CoV-2 complete sequences from Iran

Aim: SARS-CoV-2 is an emerging coronavirus that was discovered in China and rapidly spread throughout the world. The authors looked at nucleotide and amino acid variations in SARS-CoV-2 genomes, as well as phylogenetic and evolutionary events in viral genomes, in Iran. Materials & methods: All SARS-CoV-2 sequences that were publicly released between the start of the pandemic and 15 October 2021 were included. Results: The majority of mutations were found in vaccine target proteins, Spike and Nucleocapsid proteins, and nonstructural proteins. The majority of the viruses that circulated in the early stages of the pandemic belonged to the B.4 lineage. Conclusion: We discovered the prevalence of viral populations in Iran. As a result, tracking the virus’s variation in Iran and comparing it with a variety of nearby neighborhoods may reveal a pattern for future variant introductions.

The complete mitochondrial genome of Ruspolia yunnana (Orthoptera: Tettigoniidae: Conocephalinae)

The complete mitochondrial genome of Ruspolia yunnana was 15,794 bp and consisted of 37 genes and a control region. The A + T content of the mitogenome was as high as 71.3%. Most protein-coding genes started with the codon ATN and ended with stop codons TAA. Phylogenetic analyses based on mitogenomes suggested that Ruspolia yunnana and Ruspolia dubia have the closest genetic relationship. In general, this study provided meaningful genetic information for Ruspolia yunnana and validated the phylogenetic relationships within the Conocephalinae.

The complete chloroplast genome of Suaeda physophora Pall. (Chenopodiaceae)

Suaeda physophora Pall. (Chenopodiaceae) is a leaf succulent shrub species with potential usefulness as fodder for the desert animal. However, the phylogeny of S. physophora is lacking. Here, we sequenced and assembled a complete chloroplast genome of S. physophora and further reconstructed the phylogeny of Chenopodiaceae. The chloroplast genome of S. physophora is 151,104 bp in length, consisting of an 18,597 bp small single-copy (SSC), an 82,845 bp large single-copy (LSC), and a pair of 24,831 bp inverted repeat (IR) regions. The genome encodes 131 genes, including 87 protein-coding genes, 36 tRNA genes, and eight rRNA genes. Phylogenetic analysis revealed that the genus Suaeda forms a monophyletic taxon, and S. physophora is closely related to S. eltonica. Chloroplast genome and phylogenetic studies provided an essential foundation for the conservation of S. physophora.

Three types of genes underlying the Gametophyte factor1 locus cause unilateral cross incompatibility in maize

Unilateral cross incompatibility (UCI) occurs between popcorn and dent corn, and represents a critical step towards speciation. It has been reported that ZmGa1P, encoding a pectin methylesterase (PME), is a male determinant of the Ga1 locus. However, the female determinant and the genetic relationship between male and female determinants at this locus are unclear. Here, we report three different types, a total of seven linked genes underlying the Ga1 locus, which control UCI phenotype by independently affecting pollen tube growth in both antagonistic and synergistic manners. These include five pollen-expressed PME genes (ZmGa1Ps-m), a silk-expressed PME gene (ZmPME3), and another silk-expressed gene (ZmPRP3), encoding a pathogenesis-related (PR) proteins. ZmGa1Ps-m confer pollen compatibility. Presence of ZmPME3 causes silk to reject incompatible pollen. ZmPRP3 promotes incompatibility pollen tube growth and thereby breaks the blocking effect of ZmPME3. In addition, evolutionary genomics analyses suggest that the divergence of the Ga1 locus existed before maize domestication and continued during breeding improvement. The knowledge gained here deepen our understanding of the complex regulation of cross incompatibility.

Unilateral cross incompatibility (UCI) is a type of prezygotic reproductive isolation, which is associated with multiple loci in maize. Here, the authors use genetic analysis to separate the Ga1 locus into two functional components and identify seven linked genes encoding three types of proteins.

The complete chloroplast genome sequence of Gomesa flexuosa

Gomesa flexuosa (Lodd) M.W.Chase & N.H.Williams is a new species of orchid revised in 2009. In this study, the chloroplast genome of G. flexuosa was sequenced to determine its genomic characteristics and phylogenetic relationship with other related species. G. flexuosa has a chloroplast genome size of 147,764 bp, comprising 25,757 bp of two inverted repeat (IR) regions, 83,579 bp of large single-copy (LSC) region, and 12,671 bp of small single-copy (SSC) region. Moreover, the whole genome contains 73 protein-coding genes, 38 tRNA genes, and eight rRNA genes. Phylogenetic analysis showed that G. flexuosa is closely related to Oncidium sphacelatum.

The complete chloroplast genome sequence of Lithocarpus longinux (Fagaceae)

Lithocarpus longinux (Hu) Chun ex Y.C.Hsu & H.Wei Jen is a Critically Endangered tree distributed in Ma-li-po county in the Southeastern Yunnan Province China. Less than ten individuals have been found since the species was established 70 years ago. In this study, we assembled and annotated the complete chloroplast genome of L. longinux. The complete chloroplast genome of the species is 161,420 bp in length and has a GC content of 36.8%, including one large single-copy region (LSC, 90,409 bp), one small single-copy region (SSC, 19,255 bp), and two copies of inverted repeat regions (IRs, 25,878 bp). A total of 112 unique genes were detected, including 81 protein-coding genes, 29 tRNA genes, and 4 rRNA genes. Phylogenetic analysis of 31 representative chloroplast genomes of the Fagales suggests Lithocarpus is monophyletic with strong bootstrap support and that L. longinux is closely related to L. balansae.

The complete chloroplast genome and phylogenetic analysis of Salvia karwinskii (Lamiaceae)

Salvia karwinskii Benth. 1835 is a perennial herb in the Lamiaceae family native in Mexico and Central America. The complete chloroplast (cp) genome of S. karwinskii was sequenced using the Illumina platform and assembled for the first time. The complete plastid genome of S. karwinskii was 150,907 bp in length including a large single-copy (LSC) region of 82,205 bp, a small single-copy (SSC) region of 17,538 bp, and a pair of inverted repeat (IR) regions of 25,582 bp. The total GC content of this genome was 38.05%, and that of LSC, SSC, and IR regions was 36.22%, 31.77%, and 43.14%, respectively. The cp genome contained 114 unique genes, including 80 protein-coding genes, 30 tRNA genes, and four rRNA genes. The maximum-likelihood phylogenetic tree was constructed with 38 complete cp genomes, supporting a close relationship between S. karwinskii and a 10 species lineage, all of which belong to the subg. Calosphace of Salvia. The cp genome of S. karwinskii provides a foundation for further studies on genetic diversity and improving the classification system of Salvia.

The complete chloroplast genome sequence of Picrasma quassioides (D. Don) Benn. 1844 (Simaroubaceae)

Picrasma quassioides is a member of the Simaroubaceae family and is widely used as a medicinal plant. In this study, we sequenced and assembled the complete chloroplast genome of P. quassioides. The chloroplast genome is 160,015 bp in length, with a large single-copy region of 87,136 bp, a small single-copy region of 18,069 bp, and a pair of inverted repeat regions of 27,405 bp. It contains a total of 110 unique genes, including 77 protein-coding genes, 29 tRNA genes, and 4 rRNA genes. Phylogenetic analysis showed that P. quassioides clustered well with Simaroubaceae plants, Eurycoma longifolia, Leitneria floridana, and Ailanthus latissimus.

The complete chloroplast genome sequence of Dictamnus albus L

Dictamnus albus L. refers to a perennial herb with both ornamental and medicinal value. In the present study, we obtained the complete chloroplast genome sequence of D. albus through high-throughput sequencing. The length of the chloroplast genome was 157,139 bp, while the large single-copy and small single-copy regions were 84,478 bp and 18,587 bp, respectively. The pair of inverted repeat sequences was 27,037 bp, and the GC content was 38.5%. A total of 132 genes were annotated, including 87 protein-coding genes (PCGs), eight ribosomal RNA (rRNA) genes, and 37 transfer RNA (tRNA) genes. The chloroplast genomes of D. albus and eight species of Rutaceae were subjected to maximum-likelihood phylogenetic tree analysis. D. albus was found to be most closely related to Orixa japonica.

Chloroplast genome structure and phylogenetic analysis of Glycosmis parviflora (Sims) Little 1948, a folk medicinal plant featured in Lingnan Region, China

Glycosmis parviflora is the most widely spread and the most morphologically varied species of Chinese Glycosmis, and its roots and leaves serve as folk medicines. We sequenced the complete chloroplast (cp) genome of G. parviflora. The cp genome obtained was a circular DNA molecule of 159,825 bp in length, containing one large and one small single copy region (LSC and SSC) of 87,517 and 18,352 bp separated by a pair of 26,978 bp inverted repeat regions (IRs). The overall GC content of the cp genome was 38.40%. The phylogenetic analysis revealed that Glycosmis was strongly supported as a monophyletic group belonging to Clauseneae, and G. parviflora was closely related to G. pentaphylla. The results will provide the basis for the further study of molecular markers and phylogeny of G. parviflora.

NGS data vectorization, clustering, and finding key codons in SARS-CoV-2 variations

The rapid global spread and dissemination of SARS-CoV-2 has provided the virus with numerous opportunities to develop several variants. Thus, it is critical to determine the degree of the variations and in which part of the virus those variations occurred. Therefore, in this study, methods that could be used to vectorize the sequence data, perform clustering analysis, and visualize the results were proposed using machine learning methods. To conduct this study, a total of 224,073 cases of SARS-CoV-2 sequence data were collected through NCBI and GISAID, and the data were visualized using dimensionality reduction and clustering analysis models such as T-SNE and DBSCAN. The SARS-CoV-2 virus, which was first detected, was distinguished from different variations, including Omicron and Delta, in the cluster results. Furthermore, it was possible to examine which codon changes in the spike protein caused the variants to be distinguished using feature importance extraction models such as Random Forest or Shapely Value. The proposed method has the advantage of being able to analyse and visualize a large amount of data at once compared to the existing tree-based sequence data analysis. The proposed method was able to identify and visualize significant changes between the SARS-CoV-2 virus, which was first detected in Wuhan, China, in December 2019, and the newly formed mutant virus group. As a result of clustering analysis using sequence data, it was possible to confirm the formation of clusters among various variants in a two-dimensional graph, and by extracting the importance of variables, it was possible to confirm which codon changes played a major role in distinguishing variants. Furthermore, since the proposed method can handle a variety of data sequences, it can be used for all kinds of diseases, including influenza and SARS-CoV-2. Therefore, the proposed method has the potential to become widely used for the effective analysis of disease variations.

The complete chloroplast genome of Psychotria asiatica, Linnaeus 1759 (Rubiaceae)

Psychotria asiatica is a typical traditional medicinal plant. Herein, we acquired and characterized the complete chloroplast (cp) genome sequence of P. asiatica to provide genomic resources for conservation genetics and phylogenetic studies of P. asiatica. The cp genome of P. asiatica is 154,652 bp in length and consists of a large single-copy (LSC) region with 85,106 bp, a small single-copy (SSC) region with 17,960 bp, and two inverted repeat regions (IRs) with 25,793 bp. The cp genome of P. asiatica comprises 127 genes, including 82 protein-coding genes, 37 tRNA genes, and eight rRNA genes. The phylogenetic result confirmed that P. asiatica was closely related to Psychotria kirkii within the Rubiaceae.