Application of the MAFFT sequence alignment program to large data-reexamination of the usefulness of chained guide trees

Engineering Saccharomyces Cerevisiae With Novel Functional Xylose Isomerases From Rumen Microbiota for Enhanced Biofuel Production

Xylose metabolism in Saccharomyces cerevisiae remains a significant bottleneck due to the difficulty in identifying functional and efficient xylose isomerases (XI). In the present study, publicly available metagenomic and metatranscriptomic datasets of rumen microbiota from different herbivorous mammals were used to prospect novel XIs sequences. Seven putative XIs from moose, camel, cow, and sheep were cloned into a strain modified for xylose metabolism. Out of those, five XIs demonstrated activity and efficiently converted xylose into xylulose, resulting in ethanol as the final product. A XI from camel rumen microbiota exhibited a KM of 16.25 mM, indicating high substrate affinity. The strains expressing enzymes XI11 and XI12, obtained from sheep rumen microbiota, were able to deplete 40 g/L of xylose within 72 and 96 h, achieving theoretical ethanol yields of 90% and 88%, respectively. These results are comparable to those obtained with Orpinomyces sp. ukk1 XI, a benchmark enzyme previously reported as highly efficient in S. cerevisiae. This study also provides the first report on the successful expression of XIs mined from the ruminal microbiotas of sheep and camels in S. cerevisiae, expanding the perspectives for the optimization of fermentation processes and the production of lignocellulosic biofuels from xylose.

Identification of two conserved linear antigenic epitopes on the 2C protein of Senecavirus A

Senecavirus A (SVA) is a newly emerging picornavirus threatening the global swine industry, causing vesicular disease and neonatal mortality in pigs. The non-structural protein 2C of SVA is a multifunctional virulence factor. To provide robust tools for a comprehensive study of this protein's function, we successfully generated two monoclonal antibodies (mAbs; 1F9 and 6B4) by immunizing BALB/c mice with the prokaryotically expressed 2C protein as the immunogen. Indirect immunofluorescence assays confirmed that these mAbs specifically recognized the native 2C protein. Western blot analysis further substantiated their reactivity, revealing that the recognized epitopes are linear. Both 1F9 and 6B4 were characterized as IgG1/κ isotypes. Sequence analysis of the heavy and light chain variable regions showed that the framework and complementarity-determining region (CDR) sequences were entirely distinct between the two mAbs. The antigenic epitopes recognized by 1F9 and 6B4 were precisely mapped to amino acids 162DGYKGQF168 and 34LQAWINKE41, respectively, through the expression of a series of truncated forms of 2C protein. Amino acid sequence alignment of the 2C protein from global SVA strains in the GenBank database indicated that these epitopes are highly conserved. Molecular docking revealed that mAbs 1F9 and 6B4 bind to SVA 2C via hydrophobic interactions, hydrogen bonds, and salt bridges involving specific residues in their heavy and light chain CDRs. The successful development of these mAbs provides a powerful tool for the functional investigation of SVA 2C protein.

Analysis of RNA Transcribed by RNA Polymerase III from B2 SINEs in Mouse Cells

Background/Objectives: SINEs (short interspersed elements) are eukaryotic non-autonomous retrotransposons. They are transcribed by RNA polymerase III (pol III) and generate non-coding RNAs. The 3' end of many mammalian SINEs contains a polyadenylation signal (AATAAA), a pol III transcription terminator, and an A-rich tail. Studies have shown that, in human HeLa cells that have been transiently transfected with such SINEs, short pol III-generated SINE transcripts undergo polyadenylation, resulting in the addition of a long poly(A)-tail. Notably, this AAUAAA-dependent polyadenylation is not characteristic of any other transcripts synthesized by pol III. B2 SINEs, found in the genomes of mouse-like rodents, exemplify all these features. Methods: In this study, we implemented a novel approach to sequencing pol III-generated B2 transcripts from mouse cell cultures (L929 and 4T1) and organs (brain and testis). Results: Transcription occurred in 16,000-20,000 B2 copies in each cell type, 51-62% of which were transcribed in all four cell types. Effective transcription terminators (e.g., TCT>3 and T≥4) were found in approximately 40% of the transcribed B2 copies. The transcripts of these B2 copies contained a truncated terminator sequence, as pol III transcriptional arrest is known to occur within the terminator, with a poly(A)-tail immediately downstream. Such a tail could only have formed through RNA polyadenylation. Conclusions: These results demonstrate that B2 transcripts synthesized by pol III are capable of polyadenylation in mouse cells. We discuss the transcription of B2 copies with and without moderately efficient pol III terminators (TCTTT) and provide examples of the polyadenylation of such transcripts.

Genomics comparisons provide new insights into the evolution of karyotype and body patterns in Anguilliformes species

Anguilliformes species not only possess distinctive appearance (such as body elongation and absence of pectoral and/or pelvic fins), but also display diversity in chromosome number, supporting them as a suitable model for studying karyotype evolution and related molecular mechanisms of evolutionary body patterns. However, the ancestral chromosomes and evolutionary chromosomal reorganization in various eels have not been reported yet. The most regulatory or related genes of their distinctive appearance are still unknown. Here, we predicted an eel-based ancestral chromosome karyotype for the first time, and revealed multiple chromosomal fusion and fission events that reduced the ancestral chromosome number from haploid n = 21 to the commonly extant n = 19 within the Anguilla lineage. Moreover, we carried out a genome-wide comparison of two significant gene families including homeobox (Hox) and T-box (tbx), revealing genomic loss of some Hox genes (such as HoxB9β and HoxD13α) and variation of certain tbx gene (i.e., tbx5) may be responsible for the evolutionary development of pectoral fins. Interestingly, loss of certain secretory calcium-binding phosphoprotein (SCPP) genes was identified in various eel genomes, which possibly contribute to the common reduction of scales. Overall, our current findings provide new insights into evolutionary karyotype and body pattern evolution across diverse Anguilliformes species.

Analysis of genetic diversity of Zymoseptoria tritici populations in central and south-eastern Ethiopia

Septoria tritici blotch (STB), caused by the hemibiotrophic fungus Zymoseptoria tritici, is a serious threat to global wheat production, and a major bottleneck to wheat production in Ethiopia. Accurate identification and analysis of the pathogen's genetic structure helps to develop robust STB management strategies. This study aimed at molecular identification and genetic structure analysis of 200 isolates of Z. tritici representing six populations in central and south-eastern regions of Ethiopia. A total of 165 isolates were confirmed by Sanger sequencing of the internal transcribed spacer (ITS) region of nuclear DNA (rDNA) region. The pathogen's genetic structure was further examined using 12 simple sequence repeat (SSR) markers. The microsatellite markers were highly polymorphic and informative, with mean number of alleles (Na), effective alleles (Ne), Nei's gene diversity of 6.23, 2.90, and 0.59, respectively. Analysis of molecular variance (AMOVA) confirmed the presence of low population differentiation (FST = 0.02), high gene flow (Nm = 14.7), with 95% of the total genetic variation residing within populations, and leaving only 5% for the among populations. The highest genetic diversity (Number of allele = 9.33, Effective number of allele = 3.4 and Nei's gene diversity = 0.68) was observed in the Oromia special zone surrounding Finfinnee (OSZ) Z. tritici populations, followed by Arsi and North Shewa populations, indicating that these areas are ideal for multi-location wheat germplasm resistance screening, and also the pathogen genetic and genomic analyses. Cluster analyses did not clearly divide the populations into genetically separate clusters according to their geographic areas of sampling, probably due to high gene flow. The analysis revealed existence of high genetic admixture, and all the individuals shared genomic backgrounds from two subgroups (K=2). Overall, the SSR markers are highly informative and effective genetic tools for unlocking the pathogen's genetic structure. The Z. tritici populations of central and southeast Ethiopia exhibit high genetic diversity, indicating the need to deploy durable and diverse disease management strategies. North Shewa, OSZ, Arsi and West Arsi administrative zones represent hotspots for genetic and genomic analyses of Z. tritici and excellent locations for host-pathogen interaction studies, and wheat germplasm screening for STB resistance.

Gelling and reducing agents are potential carbon and energy sources in culturing of anaerobic microorganisms

The majority of microorganisms in the environment have yet to be isolated in pure cultures, and the reasons behind this phenomenon remain elusive. In this study, we investigated the possibility of the commonly used gelling agent including agar and gellan gum as a source of carbon and energy in anaerobic roll-tube cultivation from one mangrove sediment sample and two high-temperature oilfield samples. Based on growth tests and genomic evidence, anaerobic gellan degraders were retrieved from genera of Clostridium, Lacrimispora, and lineages from the rarely cultivated phylum Atribacterota. Anaerobic agarolytic microorganisms were found to be members of Bacillus and Clostridium. We also proved the role of carbon and energy sources of L-cysteine, a routine agent used to make culture media anoxic/anaerobic in both enrichment cultures and isolated strains representing Acetomicrobium, Thermodesulfovibrio, Lacrimispora, Clostridium, Bacillus, Coprothermobacter, Citrobacter, and Enterobacter. Furthermore, the isolates and enriched microbial communities utilizing L-cysteine under anaerobic conditions were mainly through L-cysteine desulfuration to pyruvate, ammonia, and sulfide. This study demonstrates that the widely used gelling and reducing agents in the basal medium can serve as carbon and energy sources for anaerobic microorganisms and thus may bias the enrichment and isolation.

IMPORTANCE

Most microbial species inhabiting natural environments have not been isolated in pure cultures using conventional media and laboratory conditions, but the reason behind this is unclear. Here, we provided a new explanation for the phenomenon, in that both the gelling agents, like agar and gellan gum, and reducing agent L-cysteine-HCl in the media provide extra carbon and energy sources to microorganisms and therefore decrease the chance in isolation specifically for the supplemented substrate which is supposed to be the sole source of carbon and energy. This result demonstrated that further improvement in the effectiveness of isolation of targeted microorganisms will be facilitated by subtracting the overlooked organic ingredients in the medium and more innovations.

Next-generation sequencing-based population genetics unravels the evolutionary history of Rhodomyrtus tomentosa in China

BACKGROUND

Rhodomyrtus tomentosa (Ait.) Hassk. is useful for its ornamental, medicinal, and ecological characteristics, and is considered a "Neglected and Underutilized Crop Species". However, our understanding of the geographic structure and evolutionary history of its wild populations is limited. To address this gap, we investigated genomic data from 284 samples of R. tomentosa from 28 wild populations in southern China.

RESULTS

The genetic diversity of populations in different regions revealed the similar trends using whole-genome and RAD-seq data, and Hainan Island having a higher genetic diversity than other regions. The 28 populations clustered into three distinct groups: (a) GROUP1 on the eastern mainland within Guangdong, Fujian, and Hunan Provinces; (b) GROUP2 on the western mainland within Guangxi and Yunnan Provinces; and (c) GROUP3 on Hainan Island. Mantel tests and redundancy analyses revealed population differentiation was affected by distance and environmental factors such as annual average radiation. Demographic history and gene flow analyses indicated the mainland populations and the Hainan Island populations diverged around 0.93 MYA, with gene flow primarily occurring from Hainan Island and the coastal regions (such as Zhanjiang in Guangdong and Fangchenggang in Guangxi) towards the mainland, reflecting an expansion trend within the species. PSMC' analyses indicated that the populations of the three groups underwent a bottleneck during the Pleistocene due to glacial-interglacial cycles and geological events. Niche analysis revealed that the ice ages caused habitat contraction for the species, and populations with higher genetic diversity are generally distributed in areas with more suitable habitats.

CONCLUSIONS

This study elucidates the current genetic distribution of the species within China and suggests that drastic Pleistocene climate change and geographical events caused population divergence and fluctuations in effective population size, shaping the current genetic distribution of R. tomentosa. These findings provide a theoretical basis for the genetic conservation and improvement of R. tomentosa.

Genomic evidence on the distribution and ecological function of Pseudomonas in hadal zone

The hadal zone is the deepest region on Earth. It serves as a depositional zone for the sinking matter from surface ocean and continental margin, aided by its unique V-shaped structure. Due to extreme depth (over 6000 m), normally only organic matter with low degradability typically reaches the bottom of the trench. Concurrently, reports have indicated highly active carbon turnover and dense bacterial cells in the Mariana Trench. There remains a cognitive gap in understanding the connection between this phenomenon and the microbial taxa along with their metabolic activities. Here, we surveyed the Pseudomonas, one of the most widely distributed bacterial genera on Earth. The result revealed widespread distribution of Pseudomonas in the hadal zones. We obtained 21 metagenome-assembled genomes (MAGs) from seawater and sediment samples of the Mariana Trench, including three novel species. Comparative genomic analysis showed that hadal Pseudomonas possess more unique ortholog groups of genes related to energy generation and substances transport, distinct from those in other marine zones. These bacteria exhibit the ability to utilize diverse electron acceptors and accumulate compatible solutes, indicating two key strategies for adaptation for high hydrostatic pressure conditions. Furthermore, predicted genomic capabilities suggest that Pseudomonas could decompose various components of organic matter, particularly aromatics, as supported by metatranscriptomic datasets. These findings significantly enhance our understanding of Pseudomonas diversity and metabolic potential, providing valuable insights into the carbon and nitrogen cycles in hadal trench ecosystems.

Design of ancestral mammalian alkaline phosphatase bearing high stability and productivity

Mammalian alkaline phosphatase (AP) is widely used in diagnostics and molecular biology but its widespread use is impaired because it is difficult to express in Escherichia coli and has low thermostability. To overcome these challenges, we employed sequence-based protein redesign methods, specifically full consensus design (FCD) and ancestral sequence reconstruction (ASR), to create APs with enhanced properties. Biochemical analyses revealed that these newly designed APs exhibited improved levels of expression in their active form and increased thermostability compared to bovine intestinal AP isozyme II (bIAPII), without impeding enzymatic activity. Notably, the activity in culture broth of the designed APs was an order of magnitude higher than that of bIAPII, and their thermal stability increased by 13°C-17°C (measured as T50). We also assessed 28 single-point mutants of bIAPII to identify regions influencing thermostability and expression level; these mutations were common in the engineered APs but not in bIAPII. Specific mutations, such as T413E and G402S, enhanced thermostability, whereas increasing the expression level required multiple mutations. This suggests that a synergistic effect is required to enhance the expression level. Mutations enhancing thermostability were located in the crown domain, while those improving expression were close to the dimer interface, indicating distinct mechanisms underpinning these enhancements.

IMPORTANCE

Sequence-based protein redesign methods, such as full consensus design (FCD) and ancestral sequence reconstruction (ASR), have the potential to construct new enzymes utilizing protein sequence data registered in a rapidly expanding sequence database. The high thermostability of these enzymes would expand their application in diagnostics and molecular biology. These enzymes have also demonstrated a high level of active expression by Escherichia coli. These characteristics make these APs attractive candidates for industrial application. In addition, different amino acid residues are primarily responsible for thermal stability and active expression, suggesting important implications for strategies for designing enzymes by FCD and ASR.

Evolutionary history of calcium-sensing receptors unveils hyper/hypocalcemia-causing mutations

Despite advancements in understanding the structure and functions of the Calcium Sensing Receptor (CaSR), gaps persist in our knowledge of the specific functions of its residues. In this study, we used phylogeny-based techniques to identify functionally equivalent orthologs of CaSR, predict residue significance, and compute specificity-determining position (SDP) scores to understand its evolutionary basis. The analysis revealed exceptional conservation of the CaSR subfamily, emphasizing the critical role of residues with high SDP scores in receptor activation and pathogenicity. To further enhance the findings, gradient-boosting trees were applied to differentiate between gain- and loss-of-function mutations responsible for hypocalcemia and hypercalcemia. Lastly, we investigated the importance of these mutations in the context of receptor activation dynamics. In summary, through comprehensive exploration of the evolutionary history of the CaSR subfamily, coupled with innovative phylogenetic methodologies, we identified activating and inactivating residues, providing valuable insights into the regulation of calcium homeostasis and its connections to associated disorders.

Are there conserved biosynthetic genes in lichens? Genome-wide assessment of terpene biosynthetic genes suggests ubiquitous distribution of the squalene synthase cluster

Lichen-forming fungi (LFF) are prolific producers of functionally and structurally diverse secondary metabolites, most of which are taxonomically exclusive and play lineage-specific roles. To date, widely distributed, evolutionarily conserved biosynthetic pathways in LFF are not known. However, this idea stems from polyketide derivatives, since most biochemical research on lichens has concentrated on polyketide synthases (PKSs). Here, we present the first systematic identification and comparison of terpene biosynthetic genes of LFF using all the available Lecanoromycete reference genomes and 22 de novo sequenced ones (111 in total, representing 60 genera and 23 families). We implemented genome mining and gene networking approaches to identify and group the biosynthetic gene clusters (BGCs) into networks of similar BGCs. Our large-scale analysis led to the identification of 724 terpene BGCs with varying degrees of pairwise similarity. Most BGCs in the dataset were unique with no similarity to a previously known fungal or bacterial BGC or among each other. Remarkably, we found two BGCs that were widely distributed in LFF. Interestingly, both conserved BGCs contain the same core gene, i.e., putatively a squalene/phytoene synthase (SQS), involved in sterol biosynthesis. This indicates that early gene duplications, followed by gene losses/gains and gene rearrangement are the major evolutionary factors shaping the composition of these widely distributed SQS BGCs across LFF. We provide an in-depth overview of these BGCs, including the transmembrane, conserved, variable and LFF-specific regions. Our study revealed that lichenized fungi do have a highly conserved BGC, providing the first evidence that a biosynthetic gene may constitute essential genes in lichens.

Telomere-to-telomere assemblies of cattle and sheep Y-chromosomes uncover divergent structure and gene content

Reference genomes of cattle and sheep have lacked contiguous assemblies of the sex-determining Y chromosome. Here, we assemble complete and gapless telomere to telomere (T2T) Y chromosomes for these species. We find that the pseudo-autosomal regions are similar in length, but the total chromosome size is substantially different, with the cattle Y more than twice the length of the sheep Y. The length disparity is accounted for by expanded ampliconic region in cattle. The genic amplification in cattle contrasts with pseudogenization in sheep suggesting opposite evolutionary mechanisms since their divergence 19MYA. The centromeres also differ dramatically despite the close relationship between these species at the overall genome sequence level. These Y chromosomes have been added to the current reference assemblies in GenBank opening new opportunities for the study of evolution and variation while supporting efforts to improve sustainability in these important livestock species that generally use sire-driven genetic improvement strategies.

Biophysical insights into sugar-dependent medium acidification promoting YfaL protein-mediated Escherichia coli self-aggregation, biofilm formation and acid stress resistance

The ability of bacteria to interact with their environment is crucial to form aggregates and biofilms, and develop a collective stress resistance behavior. Despite its environmental and medical importance, bacterial aggregation is poorly understood and mediated by few known adhesion structures. Here, we identified a new role for a surface-exposed Escherichia coli protein, YfaL, which can self-recognize and induce bacterial autoaggregation. This process occurs only under acidic conditions generated during E. coli growth in the presence of fermentable sugars. These findings were supported by electrokinetic and atomic force spectroscopy measurements, which revealed changes in the electrostatic, hydrophobic, and structural properties of YfaL-decorated cell surface upon sugar consumption. Furthermore, YfaL-mediated autoaggregation promotes biofilm formation and enhances E. coli resistance to acid stress. The prevalence and conservation of YfaL in environmental and clinical E. coli suggest strong evolutionary selection for its function inside or outside the host. Overall, our results emphasize the importance of environmental parameters such as low pH as physicochemical cues influencing bacterial adhesion and aggregation, affecting E. coli and potentially other bacteria's resistance to environmental stress.

Detection and Molecular Characterization of GI-1 and GI-23 Avian Infectious Bronchitis Virus in Broilers Indicate the Emergence of New Genotypes in Bolivia

Infectious Bronchitis Virus (IBV) is a major threat to the poultry industry worldwide, causing significant economic losses. While the virus's genetic structure is well understood, the specific strains circulating in Bolivia have remained uncharacterized until now. This study aimed to identify and characterize new IBV strains in Bolivia. Tissue samples from broilers exhibiting clinical signs of Infectious Bronchitis were screened to detect IBV using real-time RT-PCR (RT-qPCR). Positive samples with low cycle threshold (Ct) values were selected for sequencing the full S1 gene. Of the 12 samples analyzed, 10 were determined to be positive for IBV. However, only four samples yielded sufficient genetic material for sequencing and subsequent phylogenetic analysis. The results revealed the presence of GI-1 and GI-23 lineages, both belonging to genotype I (GI). The GI-1 lineage showed >99% sequence identity to the H120 and Massachusetts vaccine strains, suggesting a close relationship. In contrast, the GI-23 lineage clustered with other IBV strains, showing a distinct subclade that is genetically distant from Brazilian strains. No evidence of recombination was found. Furthermore, amino acid substitution analysis identified specific mutations in the S1 subunit, particularly in the hypervariable regions 1, 2, and 3. These mutations could potentially alter the virus's antigenicity, leading to reduced vaccine efficacy. The findings of this study highlight the importance of continued and broad genomic surveillance of circulating IBV strains and the need to improve vaccination strategies in Bolivia.

learnMSA2: deep protein multiple alignments with large language and hidden Markov models

MOTIVATION

For the alignment of large numbers of protein sequences, tools are predominant that decide to align two residues using only simple prior knowledge, e.g. amino acid substitution matrices, and using only part of the available data. The accuracy of state-of-the-art programs declines with decreasing sequence identity and when increasingly large numbers of sequences are aligned. Recently, transformer-based deep-learning models started to harness the vast amount of protein sequence data, resulting in powerful pretrained language models with the main purpose of generating high-dimensional numerical representations, embeddings, for individual sites that agglomerate evolutionary, structural, and biophysical information.

RESULTS

We extend the traditional profile hidden Markov model so that it takes as inputs unaligned protein sequences and the corresponding embeddings. We fit the model with gradient descent using our existing differentiable hidden Markov layer. All sequences and their embeddings are jointly aligned to a model of the protein family. We report that our upgraded HMM-based aligner, learnMSA2, combined with the ProtT5-XL protein language model aligns on average almost 6% points more columns correctly than the best amino acid-based competitor and scales well with sequence number. The relative advantage of learnMSA2 over other programs tends to be greater when the sequence identity is lower and when the number of sequences is larger. Our results strengthen the evidence on the rich information contained in protein language models' embeddings and their potential downstream impact on the field of bioinformatics. Availability and implementation:  https://github.com/Gaius-Augustus/learnMSA, PyPI and Bioconda, evaluation: https://github.com/felbecker/snakeMSA.

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

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

Sacbrood viruses and select Lake Sinai virus variants dominated Apis mellifera colonies symptomatic for European foulbrood

European foulbrood (EFB) is a prevalent disease in the European honey bee (Apis mellifera) in the United States, which can lead to colony decline and collapse. The bacterial components of EFB are well-studied, but the diversity of viral infections within infected colonies has not been explored. In this study, we use meta-transcriptomics sequencing of 12 honey bee hives, symptomatic (+, n = 6) and asymptomatic (-, n = 6) for EFB, to investigate viral infection associated with the disease. We assembled 41 viral genomes, belonging to three families (Iflaviridae, Dicistroviridae, and Sinhaliviridae), all previously reported in honey bees, including Lake Sinai virus, deformed wing virus, sacbrood virus, Black queen cell virus, and Israeli acute paralysis virus. In colonies with severe EFB, we observed a higher occurrence of viral genomes (34 genomes) in contrast to fewer recovered from healthy colonies (seven genomes) and a complete absence of Dicistroviridae genomes.We observed specific Lake Sinai virus clades associated exclusively with EFB + or EFB - colonies, in addition to EFB-afflicted colonies that exhibited an increase in relative abundance of sacbrood viruses. Multivariate analyses highlighted that a combination of site and EFB disease status influenced RNA virome composition, while EFB status alone did not significantly impact it, presenting a challenge for comparisons between colonies kept in different yards. These findings contribute to the understanding of viral dynamics in honey bee colonies compromised by EFB and underscore the need for future investigations to consider viral composition when investigating EFB.IMPORTANCEThis study on the viromes of honey bee colonies affected by European foulbrood (EFB) sheds light on the dynamics of viral populations in bee colonies in the context of a prevalent bacterial brood disease. The identification of distinct Lake Sinai virus and sacbrood virus clades associated with colonies affected by severe EFB suggests a potential connection between viral composition and disease status, emphasizing the need for further investigation into the role of viruses during EFB infection. The observed increase in sacbrood viruses during EFB infection suggests a potential viral dysbiosis, with potential implications for honey bee brood health. These findings contribute valuable insights related to beekeeping practices, offering a foundation for future research aimed at understanding and mitigating the impact of bacterial and viral infection in commercial honey bee operations and the management of EFB.

A Naturally Active Spy Transposon Discovered from the Insect Genome of Colletes gigas as a Promising Novel Gene Transfer Tool

Novel active DNA transposons, such as Spy transposons from the PHIS superfamily, are identified through bioinformatics in this study. The native transposases cgSpy and cvSpy displayed transposition activities of approximately 85% and 35% compared to the hyperactive piggyBac transposase (hyPB). The cgSpy transposon showed unique characteristics, including a lack of overproduction inhibition and reduced efficiency for insertion sizes between 3.1 to 8.5 kb. Integration preferences of cgSpy are found in genes and regulatory regions, making it suitable for genetic manipulation. Evaluation in T-cell engineering demonstrated that cgSpy-mediated chimeric antigen receptor (CAR) modification is comparable to the PB system, indicating its potential utility in cell therapy. This study unveils the promising application of the active native transposase, Spy, from Colletes gigas, as a valuable tool for genetic engineering, particularly in T-cell manipulation.

Endogenous Bornavirus-like Elements in Bats: Evolutionary Insights from the Conserved Riboviral L-Gene in Microbats and Its Antisense Transcription in Myotis daubentonii

Bats are ecologically diverse vertebrates characterized by their ability to host a wide range of viruses without apparent illness and the presence of numerous endogenous viral elements (EVEs). EVEs are well preserved, expressed, and may affect host biology and immunity, but their role in bat immune system evolution remains unclear. Among EVEs, endogenous bornavirus-like elements (EBLs) are bornavirus sequences integrated into animal genomes. Here, we identified a novel EBL in the microbat Myotis daubentonii, EBLL-Cultervirus.10-MyoDau (short name is CV.10-MyoDau) that shows protein-level conservation with the L-protein of a Cultervirus (Wuhan sharpbelly bornavirus). Surprisingly, we discovered a transcript on the antisense strand comprising three exons, which we named AMCR-MyoDau. The active transcription in Myotis daubentonii tissues of AMCR-MyoDau, confirmed by RNA-Seq analysis and RT-PCR, highlights its potential role during viral infections. Using comparative genomics comprising 63 bat genomes, we demonstrate nucleotide-level conservation of CV.10-MyoDau and AMCR-MyoDau across various bat species and its detection in 22 Yangochiropera and 12 Yinpterochiroptera species. To the best of our knowledge, this marks the first occurrence of a conserved EVE shared among diverse bat species, which is accompanied by a conserved antisense transcript. This highlights the need for future research to explore the role of EVEs in shaping the evolution of bat immunity.

From dusk till dawn: cell cycle progression in the red seaweed Gracilariopsis chorda (Rhodophyta)

The conserved eukaryotic functions of cell cycle genes have primarily been studied using animal/plant models and unicellular algae. Cell cycle progression and its regulatory components in red (Rhodophyta) seaweeds are poorly understood. We analyzed diurnal gene expression data to investigate the cell cycle in the red seaweed Gracilariopsis chorda. We identified cell cycle progression and transitions in G. chorda which are induced by interactions of key regulators such as E2F/DP, RBR, cyclin-dependent kinases, and cyclins from dusk to dawn. However, several typical CDK inhibitor proteins are absent in red seaweeds. Interestingly, the G1-S transition in G. chorda is controlled by delayed transcription of GINS subunit 3. We propose that the delayed S phase entry in this seaweed may have evolved to minimize DNA damage (e.g., due to UV radiation) during replication. Our results provide important insights into cell cycle-associated physiology and its molecular mechanisms in red seaweeds.

Evolutionary and Expression Analysis of the Pig MAGE Gene Family

The melanoma-associated antigen (MAGE) family found in eukaryotes plays a crucial role in cell proliferation and differentiation, spermatogenesis, neural development, etc. This study explored the validation and evolution of MAGE genes in eukaryotic genomes and their distribution and expression patterns in pigs. In total, 249 MAGE genes were found on 13 eukaryotic species. In total, 33, 25, and 18 genes were located on human, mouse, and pig genomes, respectively. We found eight, four, and three tandemly duplicated gene clusters on the human, mouse, and pig genomes, respectively. The majority of MAGE genes in mammals are located on the X chromosome. According to the phylogenetic analysis, the MAGE family genes were classified into 11 subfamilies. The NDN gene in zebrafish (DreNDN) was the root of this evolutionary tree. In total, 10 and 11 MAGE genes on human and mouse genomes, respectively, exhibited a collinearity relationship with the MAGE genes on pig genomes. Taking the MAGE family genes in pigs, the MAGE subfamilies had similar gene structures, protein motifs, and biochemical attributes. Using the RNA-seq data of Duroc pigs and Rongchang pigs, we detected that the expression of type I MAGE genes was higher in reproductive tissues, but type II MAGE genes were predominantly expressed in the brain tissue. These findings are a valuable resource for gaining insight into the evolution and expression of the MAGE family genes.

Heterologous survey of 130 DNA transposons in human cells highlights their functional divergence and expands the genome engineering toolbox

Experimental studies on DNA transposable elements (TEs) have been limited in scale, leading to a lack of understanding of the factors influencing transposition activity, evolutionary dynamics, and application potential as genome engineering tools. We predicted 130 active DNA TEs from 102 metazoan genomes and evaluated their activity in human cells. We identified 40 active (integration-competent) TEs, surpassing the cumulative number (20) of TEs found previously. With this unified comparative data, we found that the Tc1/mariner superfamily exhibits elevated activity, potentially explaining their pervasive horizontal transfers. Further functional characterization of TEs revealed additional divergence in features such as insertion bias. Remarkably, in CAR-T therapy for hematological and solid tumors, Mariner2_AG (MAG), the most active DNA TE identified, largely outperformed two widely used vectors, the lentiviral vector and the TE-based vector SB100X. Overall, this study highlights the varied transposition features and evolutionary dynamics of DNA TEs and increases the TE toolbox diversity.

Ere, a Family of Short Interspersed Elements in the Genomes of Odd-Toed Ungulates (Perissodactyla)

Short Interspersed Elements (SINEs) are eukaryotic retrotransposons transcribed by RNA polymerase III (pol III). Many mammalian SINEs (T+ SINEs) contain a polyadenylation signal (AATAAA), a pol III transcription terminator, and an A-rich tail in their 3'-end. The RNAs of such SINEs have the capacity for AAUAAA-dependent polyadenylation, which is unique to pol III-generated transcripts. The structure, evolution, and polyadenylation of the Ere SINE of ungulates (horses, rhinos, and tapirs) were investigated in this study. A bioinformatics analysis revealed the presence of up to ~4 × 105 Ere copies in representatives of all three families. These copies can be classified into two large subfamilies, EreA and EreB, the former distinguished by an additional 60 bp sequence. The 3'-end of numerous EreA and all EreB copies exhibit a 50 bp sequence designated as a terminal domain (TD). The Ere family can be further subdivided into subfamilies EreA_0TD, EreA_1TD, EreB_1TD, and EreB_2TD, depending on the presence and number of terminal domains (TDs). Only EreA_0TD copies can be assigned to T+ SINEs as they contain the AATAAA signal and the TCTTT transcription terminator. The analysis of young Ere copies identified by comparison with related perissodactyl genomes revealed that EreA_0TD and, to a much lesser extent, EreB_2TD have retained retrotranspositional activity in the recent evolution of equids and rhinoceroses. The targeted mutagenesis and transfection of HeLa cells were used to identify sequences in equine EreA_0TD that are critical for the polyadenylation of its pol III transcripts. In addition to AATAAA and the transcription terminator, two sites in the 3' half of EreA, termed the β and τ signals, were found to be essential for this process. The evolution of Ere, with a particular focus on the emergence of T+ SINEs, as well as the polyadenylation signals are discussed in comparison with other T+ SINEs.

Evidence for cross-species transmission of human coronavirus OC43 through bioinformatics and modeling infections in porcine intestinal organoids

Cross-species transmission of coronaviruses has been continuously posing a major challenge to public health. Pigs, as the major animal reservoirs for many zoonotic viruses, frequently mediate viral transmission to humans. This study comprehensively mapped the relationship between human and porcine coronaviruses through in-depth bioinformatics analysis. We found that human coronavirus OC43 and porcine coronavirus PHEV share a close phylogenetic relationship, evidenced by high genomic homology, similar codon usage patterns and comparable tertiary structure in spike proteins. Inoculation of infectious OC43 viruses in organoids derived from porcine small and large intestine demonstrated that porcine intestinal organoids (pIOs) are highly susceptible to human coronavirus OC43 infection and support infectious virus production. Using transmission electron microscopy, we visualized OC43 viral particles in both intracellular and extracellular compartments, and observed abnormalities of multiple organelles in infected organoid cells. Robust OC43 infections in pIOs result in a significant reduction of organoids viability and widespread cell death. This study bears essential implications for better understanding the evolutionary origin of human coronavirus OC43, and provides a proof-of-concept for using pIOs as a model to investigate cross-species transmission of human coronavirus.

SARS-CoV-2 variant introduction following spring break travel and transmission mitigation strategies

BACKGROUND

University spring break carries a two-pronged SARS-CoV-2 variant transmission risk. Circulating variants from universities can spread to spring break destinations, and variants from spring break destinations can spread to universities and surrounding communities. Therefore, it is critical to implement SARS-CoV-2 variant surveillance and testing strategies to limit community spread before and after spring break to mitigate virus transmission and facilitate universities safely returning to in-person teaching.

METHODS

We examined the SARS-CoV-2 positivity rate and changes in variant lineages before and after the university spring break for two consecutive years. 155 samples were sequenced across four time periods: pre- and post-spring break 2021 and pre- and post-spring break 2022; following whole genome sequencing, samples were assigned clades. The clades were then paired with positivity and testing data from over 50,000 samples.

RESULTS

In 2021, the number of variants in the observed population increased from four to nine over spring break, with variants of concern being responsible for most of the cases; Alpha percent composition increased from 22.2% to 56.4%. In 2022, the number of clades in the population increased only from two to three, all of which were Omicron or a sub-lineage of Omicron. However, phylogenetic analysis showed the emergence of distantly related sub-lineages. 2022 saw a greater increase in positivity than 2021, which coincided with a milder mitigation strategy. Analysis of social media data provided insight into student travel destinations and how those travel events may have impacted spread.

CONCLUSIONS

We show the role that repetitive testing can play in transmission mitigation, reducing community spread, and maintaining in-person education. We identified that distantly related lineages were brought to the area after spring break travel regardless of the presence of a dominant variant of concern.

Discovery and structural mechanism of DNA endonucleases guided by RAGATH-18-derived RNAs

CRISPR-Cas systems and IS200/IS605 transposon-associated TnpBs have been utilized for the development of genome editing technologies. Using bioinformatics analysis and biochemical experiments, here we present a new family of RNA-guided DNA endonucleases. Our bioinformatics analysis initially identifies the stable co-occurrence of conserved RAGATH-18-derived RNAs (reRNAs) and their upstream IS607 TnpBs with an average length of 390 amino acids. IS607 TnpBs form programmable DNases through interaction with reRNAs. We discover the robust dsDNA interference activity of IS607 TnpB systems in bacteria and human cells. Further characterization of the Firmicutes bacteria IS607 TnpB system (ISFba1 TnpB) reveals that its dsDNA cleavage activity is remarkably sensitive to single mismatches between the guide and target sequences in human cells. Our findings demonstrate that a length of 20 nt in the guide sequence of reRNA achieves the highest DNA cleavage activity for ISFba1 TnpB. A cryo-EM structure of the ISFba1 TnpB effector protein bound by its cognate RAGATH-18 motif-containing reRNA and a dsDNA target reveals the mechanisms underlying reRNA recognition by ISFba1 TnpB, reRNA-guided dsDNA targeting, and the sensitivity of the ISFba1 TnpB system to base mismatches between the guide and target DNA. Collectively, this study identifies the IS607 TnpB family of compact and specific RNA-guided DNases with great potential for application in gene editing.

The complete mitochondrial genome of Hypsipetes amaurotis (Passeriformes: Pycnonotidae)

The brown-eared bulbul (Hypsipetes amaurotis) is a medium-sized songbird native to East Asia and characterized by its prominent reddish-brown ear-coverts. Previous studies on it have primarily been from the taxonomic and morphological aspects, with limited research in the realm of molecular biology. In this study, we sequenced and annotated the complete mitochondrial genome of H. amaurotis, which was the first reported complete mitogenome of the genus Hypsipetes. The mitogenome of H. amaurotis is 17,871 bp in length and was predicted to encode 37 typical mitochondrial genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs). Specifically, this mitogenome contains two D-loop control regions that are of similar length and sequencing pattern. A total of 8 Pycnonotidae and six outgroup taxa were used to determine the phylogenetic placement with two methods: Maximum Likelihood Approximation (IQ-TREE) and Bayesian inference (MrBayes). Our findings reveal that H. amaurotis is phylogenetically closely related to Ixos mcclellandii. The outcomes are generally consistent with the phylogenetic trees constructed in previous studies. The data gathered from this research provides valuable insights for future genomic investigations into the evolution, ecology, and conservation of this species.

Phylogenetic and transcriptomic characterization of insulin and growth factor receptor tyrosine kinases in crustaceans

Receptor tyrosine kinases (RTKs) mediate the actions of growth factors in metazoans. In decapod crustaceans, RTKs are implicated in various physiological processes, such molting and growth, limb regeneration, reproduction and sexual differentiation, and innate immunity. RTKs are organized into two main types: insulin receptors (InsRs) and growth factor receptors, which include epidermal growth factor receptor (EGFR), fibroblast growth factor receptor (FGFR), vascular endothelial growth factor receptor (VEGFR), and platelet-derived growth factor receptor (PDGFR). The identities of crustacean RTK genes are incomplete. A phylogenetic analysis of the CrusTome transcriptome database, which included all major crustacean taxa, showed that RTK sequences segregated into receptor clades representing InsR (72 sequences), EGFR (228 sequences), FGFR (129 sequences), and PDGFR/VEGFR (PVR; 235 sequences). These four receptor families were distinguished by the domain organization of the extracellular N-terminal region and motif sequences in the protein kinase catalytic domain in the C-terminus or the ligand-binding domain in the N-terminus. EGFR1 formed a single monophyletic group, while the other RTK sequences were divided into subclades, designated InsR1-3, FGFR1-3, and PVR1-2. In decapods, isoforms within the RTK subclades were common. InsRs were characterized by leucine-rich repeat, furin-like cysteine-rich, and fibronectin type 3 domains in the N-terminus. EGFRs had leucine-rich repeat, furin-like cysteine-rich, and growth factor IV domains. N-terminal regions of FGFR1 had one to three immunoglobulin-like domains, whereas FGFR2 had a cadherin tandem repeat domain. PVRs had between two and five immunoglobulin-like domains. A classification nomenclature of the four RTK classes, based on phylogenetic analysis and multiple sequence alignments, is proposed.

The first complete T2T Assemblies of Cattle and Sheep Y-Chromosomes uncover remarkable divergence in structure and gene content

Reference genomes of cattle and sheep have lacked contiguous assemblies of the sex-determining Y chromosome. We assembled complete and gapless telomere to telomere (T2T) Y chromosomes for these species. The pseudo-autosomal regions were similar in length, but the total chromosome size was substantially different, with the cattle Y more than twice the length of the sheep Y. The length disparity was accounted for by expanded ampliconic region in cattle. The genic amplification in cattle contrasts with pseudogenization in sheep suggesting opposite evolutionary mechanisms since their divergence 18MYA. The centromeres also differed dramatically despite the close relationship between these species at the overall genome sequence level. These Y chromosome have been added to the current reference assemblies in GenBank opening new opportunities for the study of evolution and variation while supporting efforts to improve sustainability in these important livestock species that generally use sire-driven genetic improvement strategies.

Multiple horizontal gene transfer events have shaped plant glycosyl hydrolase diversity and function

Plant glycosyl hydrolases (GHs) play a crucial role in selectively breaking down carbohydrates and glycoconjugates during various cellular processes, such as reserve mobilization, pathogen defense, and modification/disassembly of the cell wall. In this study, we examined the distribution of GH genes in the Archaeplastida supergroup, which encompasses red algae, glaucophytes, and green plants. We identified that the GH repertoire expanded from a few tens of genes in early archaeplastidians to over 400 genes in modern angiosperms, spanning 40 GH families in land plants. Our findings reveal that major evolutionary transitions were accompanied by significant changes in the GH repertoire. Specifically, we identified at least 23 GH families acquired by green plants through multiple horizontal gene transfer events, primarily from bacteria and fungi. We found a significant shift in the subcellular localization of GH activity during green plant evolution, with a marked increase in extracellular-targeted GH proteins associated with the diversification of plant cell wall polysaccharides and defense mechanisms against pathogens. In conclusion, our study sheds light on the macroevolutionary processes that have shaped the GH repertoire in plants, highlighting the acquisition of GH families through horizontal transfer and the role of GHs in plant adaptation and defense mechanisms.

Molecular characterization of human bocavirus in municipal wastewaters using amplicon target sequencing

Human bocavirus (HBoV) is an emerging health concern worldwide, associated with range of clinical manifestations, including gastroenteritis and respiratory infections. Therefore, it is crucial to comprehend and minimize their prevalence in different systems. In this study, we conducted regular sampling throughout the year in two different sizes and work processes of wastewater treatment plants (WWTPs) in Tianjin, China. Our objective was to investigate the occurrence, prevalence, and endurance of HBoV in wastewater, while also evaluating the efficacy of amplicon target sequencing in directly detecting HBoV in wastewater. At two WWTPs, HBoV2 (45.51 %-45.67 %) and HBoV3 (38.30 %-40.25 %) were the most common genotypes identified, and the mean concentration range of HBoV was 2.54-7.40 log10 equivalent copies/l as determined by multiplex real-time quantitative PCR assay. A positive rate of HBoV was found in 96.6 % (29/30) samples of A-WWTP, and 96.6 % (26/27) samples of B-WWTP. The phylogenetic analysis indicated that the nucleotide similarity between the HBoV DNA sequences to the reference HBoV sequences published globally ranged from 90.14 %-100 %. A significant variation in the read abundance of HBoV2 and HBoV3 in two wastewater treatment plants (p < 0.05) was detected, specifically in the Winter and Summer seasons. The findings revealed a strong correlation between the genotypes detected in wastewater and the clinical data across various regions in China. In addition, it is worth mentioning that HBoV4 was exclusively detected in wastewater and not found in the clinical samples from patients. This study highlights the high prevalence of human bocavirus in municipal wastewater. This finding illustrates that amplicon target sequencing can amplify a wide variety of viruses, enabling the identification of newly discovered viruses.

Evolution of Endogenous Retroviruses in the Subfamily of Caprinae

The interest in endogenous retroviruses (ERVs) has been fueled by their impact on the evolution of the host genome. In this study, we used multiple pipelines to conduct a de novo exploration and annotation of ERVs in 13 species of the Caprinae subfamily. Through analyses of sequence identity, structural organization, and phylogeny, we defined 28 ERV groups within Caprinae, including 19 gamma retrovirus groups and 9 beta retrovirus groups. Notably, we identified four recent and potentially active groups prevalent in the Caprinae genomes. Additionally, our investigation revealed that most long noncoding genes (lncRNA) and protein-coding genes (PC) contain ERV-derived sequences. Specifically, we observed that ERV-derived sequences were present in approximately 75% of protein-coding genes and 81% of lncRNA genes in sheep. Similarly, in goats, ERV-derived sequences were found in approximately 74% of protein-coding genes and 75% of lncRNA genes. Our findings lead to the conclusion that the majority of ERVs in the Caprinae genomes can be categorized as fossils, representing remnants of past retroviral infections that have become permanently integrated into the genomes. Nevertheless, the identification of the Cap_ERV_20, Cap_ERV_21, Cap_ERV_24, and Cap_ERV_25 groups indicates the presence of relatively recent and potentially active ERVs in these genomes. These particular groups may contribute to the ongoing evolution of the Caprinae genome. The identification of putatively active ERVs in the Caprinae genomes raises the possibility of harnessing them for future genetic marker development.

Diurnal Rhythms in the Red Seaweed Gracilariopsis chorda are Characterized by Unique Regulatory Networks of Carbon Metabolism

Cellular and physiological cycles are driven by endogenous pacemakers, the diurnal and circadian rhythms. Key functions such as cell cycle progression and cellular metabolism are under rhythmic regulation, thereby maintaining physiological homeostasis. The photoreceptors phytochrome and cryptochrome, in response to light cues, are central input pathways for physiological cycles in most photosynthetic organisms. However, among Archaeplastida, red algae are the only taxa that lack phytochromes. Current knowledge about oscillatory rhythms is primarily derived from model species such as Arabidopsis thaliana and Chlamydomonas reinhardtii in the Viridiplantae, whereas little is known about these processes in other clades of the Archaeplastida, such as the red algae (Rhodophyta). We used genome-wide expression profiling of the red seaweed Gracilariopsis chorda and identified 3,098 rhythmic genes. Here, we characterized possible cryptochrome-based regulation and photosynthetic/cytosolic carbon metabolism in this species. We found a large family of cryptochrome genes in G. chorda that display rhythmic expression over the diurnal cycle and may compensate for the lack of phytochromes in this species. The input pathway gates regulatory networks of carbon metabolism which results in a compact and efficient energy metabolism during daylight hours. The system in G. chorda is distinct from energy metabolism in most plants, which activates in the dark. The green lineage, in particular, land plants, balance water loss and CO2 capture in terrestrial environments. In contrast, red seaweeds maintain a reduced set of photoreceptors and a compact cytosolic carbon metabolism to thrive in the harsh abiotic conditions typical of intertidal zones.

George K, Kajon AE. A Fulminant Case of Adenovirus Genotype C108 Infection in a Pediatric Stem Cell Transplant Recipient with x-Linked Lymphoproliferative Syndrome Type 1

A 3-year-old male with X-linked lymphoproliferative syndrome type 1 underwent an unrelated umbilical cord blood transplant (UUCBT). The week prior to transplant the patient tested positive for adenovirus (HAdV) with a viral load of <190 copies/mL and was started on cidofovir. UUCBT proceeded as scheduled, and the patient engrafted on day +19. The patient's HAdV load in serum continued to rise with resulting hepatic dysfunction, despite ongoing therapy with cidofovir and HAdV specific T-cell infusions. The patient died 6 months after transplantation having never cleared the virus. Next generation whole genome sequencing and sequence data analyses identified an intertypic recombinant HAdV-C P1H2F2 closely related (99.6% similarity) to genotype C108 in the isolates from three blood specimens obtained during the last week of life. Incidentally, the de novo assembly strategy enabled the detection of an adeno-associated virus type 2 (AAV2) genome in the DNA purified from the plasma isolates. Proteotyping analysis revealed minor differences in the predicted amino acid sequences for E1A, E1B 19K, E1B 55K, DNA polymerase, penton base, and fiber. None of the mutations previously described for HAdV-C5 variants resistant to cidofovir were identified. In silico restriction enzyme analysis revealed a distinct Sac I profile for the identified virus, supporting its designation as a C108 variant.

In silico analysis of crustacean hyperglycemic hormone family G protein-coupled receptor candidates

Ecdysteroid molting hormone synthesis is directed by a pair of molting glands or Y-organs (YOs), and this synthesis is inhibited by molt-inhibiting hormone (MIH). MIH is a member of the crustacean hyperglycemic hormone (CHH) neuropeptide superfamily, which includes CHH and insect ion transport peptide (ITP). It is hypothesized that the MIH receptor is a Class A (Rhodopsin-like) G protein-coupled receptor (GPCR). The YO of the blackback land crab, Gecarcinus lateralis, expresses 49 Class A GPCRs, three of which (Gl-CHHR-A9, -A10, and -A12) were provisionally assigned as CHH-like receptors. CrusTome, a transcriptome database assembled from 189 crustaceans and 12 ecdysozoan outgroups, was used to deorphanize candidate MIH/CHH GPCRs, relying on sequence homology to three functionally characterized ITP receptors (BNGR-A2, BNGR-A24, and BNGR-A34) in the silk moth, Bombyx mori. Phylogenetic analysis and multiple sequence alignments across major taxonomic groups revealed extensive expansion and diversification of crustacean A2, A24, and A34 receptors, designated CHH Family Receptor Candidates (CFRCs). The A2 clade was divided into three subclades; A24 clade was divided into five subclades; and A34 was divided into six subclades. The subclades were distinguished by conserved motifs in extracellular loop (ECL) 2 and ECL3 in the ligand-binding region. Eleven of the 14 subclades occurred in decapod crustaceans. In G. lateralis, seven CFRC sequences, designated Gl-CFRC-A2α1, -A24α, -A24β1, -A24β2, -A34α2, -A34β1, and -A34β2, were identified; the three A34 sequences corresponded to Gl-GPCR-A12, -A9, and A10, respectively. ECL2 in all the CFRC sequences had a two-stranded β-sheet structure similar to human Class A GPCRs, whereas the ECL2 of decapod CFRC-A34β1/β2 had an additional two-stranded β-sheet. We hypothesize that this second β-sheet on ECL2 plays a role in MIH/CHH binding and activation, which will be investigated further with functional assays.

PiggyBac Transposon Mining in the Small Genomes of Animals

TEs, including DNA transposons, are major contributors of genome expansions, and have played a very significant role in shaping the evolution of animal genomes, due to their capacity to jump from one genomic position to the other. In this study, we investigated the evolution landscapes of PB transposons, including their distribution, diversity, activity and structure organization in 79 species of small (compact) genomes of animals comprising both vertebrate and invertebrates. Overall, 212 PB transposon types were detected from almost half (37) of the total number of the small genome species (79) investigated. The detected PB transposon types, which were unevenly distributed in various genera and phyla, have been classified into seven distinct clades or families with good bootstrap support (>80%). The PB transposon types that were identified have a length ranging from 1.23 kb to 9.51 kb. They encode transposases of approximately ≥500 amino acids in length, and possess terminal inverted repeats (TIRs) ranging from 4 bp to 24 bp. Though some of the transposon types have long TIRs (528 bp), they still maintain the consistent and reliable 4 bp target site duplication (TSD) of TTAA. However, PiggyBac-2_Cvir transposon originating from the Crassostrea virginica species exhibits a unique TSD of TATG. The TIRs of the transposons in all the seven families display high divergence, with a highly conserved 5' end motif. The core transposase domains (DDD) were better conserved among the seven different families compared to the other protein domains, which were less prevalent in the vertebrate genome. The divergent evolution dynamics analysis also indicated that the majority of the PB transposon types identified in this study are either relatively young or old, with some being active. Additionally, numerous invasions of PB transposons were found in the genomes of both vertebrate and invertebrate animals. The data reveals that the PB superfamily is widely distributed in these species. PB transposons exhibit high diversity and activity in the small genomes of animals, and might play a crucial role in shaping the evolution of these small genomes of animals.

Phylogenomics Reveals High Levels of Incomplete Lineage Sorting at the Ancestral Nodes of the Macaque Radiation

The genus Macaca includes 23 species assigned into 4 to 7 groups. It exhibits the largest geographic range and represents the most successful example of adaptive radiation of nonhuman primates. However, intrageneric phylogenetic relationships among species remain controversial and have not been resolved so far. In this study, we conducted a phylogenomic analysis on 16 newly generated and 8 published macaque genomes. We found strong evidence supporting the division of this genus into 7 species groups. Incomplete lineage sorting (ILS) was the primary factor contributing to the discordance observed among gene trees; however, we also found evidence of hybridization events, specifically between the ancestral arctoides/sinica and silenus/nigra lineages that resulted in the hybrid formation of the fascicularis/mulatta group. Combined with fossil data, our phylogenomic data were used to establish a scenario for macaque radiation. These findings provide insights into ILS and potential ancient introgression events that were involved in the radiation of macaques, which will lead to a better understanding of the rapid speciation occurring in nonhuman primates.

Chromosome-level genome assembly of the largefin longbarbel catfish (Hemibagrus macropterus)

The largefin longbarbel catfish, Hemibagrus macropterus, is an economically important fish species in southwestern China, with males growing faster than females. This study presents a high-quality chromosome-level genome assembly of the largefin longbarbel catfish, generated by integrating Illumina short reads, PacBio HiFi long reads, and Hi-C data. The assembled genome size was 858.5 Mb, with a contig and scaffold N50 of 5.8 Mb and 28.4 Mb, respectively. A total of 656 contigs were successfully anchored to 30 pseudochromosomes with a BUSCO score of 97.7%, consistent with the number of chromosomes analyzed by karyotype. The genome contained 29.5% repeat sequences, and a predicted total of 26,613 protein-coding genes, of which 25,769 (96.8%) were functionally annotated in different databases. Evolutionary analysis showed that H. macropterus was most closely related to H. wyckioides, with a divergence time of approximately 16.3 million years. Chromosomal syntenic relationships among H. macropterus, H. wyckioides, and Pelteobagrus fulvidraco revealed a one-to-one relationship for most chromosomes, except for break, fission, and inversion of some chromosomes. The first high-quality reference genome will not only provide a valuable genetic resource for the study of sex determination mechanisms and genetic breeding of largefin longbarbel catfish, but also contribute to comparative analyses of genome and chromosome evolution within Siluriformes.

Structural and functional analysis of hyper-thermostable ancestral L-amino acid oxidase that can convert Trp derivatives to D-forms by chemoenzymatic reaction

Production of D-amino acids (D-AAs) on a large-scale enables to provide precursors of peptide therapeutics. In this study, we designed a novel L-amino acid oxidase, HTAncLAAO2, by ancestral sequence reconstruction, exhibiting high thermostability and long-term stability. The crystal structure of HTAncLAAO2 was determined at 2.2 Å by X-ray crystallography, revealing that the enzyme has an octameric form like a "ninja-star" feature. Enzymatic property analysis demonstrated that HTAncLAAO2 exhibits three-order larger kcat/Km values towards four L-AAs (L-Phe, L-Leu, L-Met, and L-Ile) than that of L-Trp. Through screening the variants, we obtained the HTAncLAAO2(W220A) variant, which shows a > 6-fold increase in kcat value toward L-Trp compared to the original enzyme. This variant applies to synthesizing enantio-pure D-Trp derivatives from L- or rac-forms at a preparative scale. Given its excellent properties, HTAncLAAO2 would be a starting point for designing novel oxidases with high activity toward various amines and AAs.

Navigating the Landscape: A Comprehensive Review of Current Virus Databases

Viruses are abundant and diverse entities that have important roles in public health, ecology, and agriculture. The identification and surveillance of viruses rely on an understanding of their genome organization, sequences, and replication strategy. Despite technological advancements in sequencing methods, our current understanding of virus diversity remains incomplete, highlighting the need to explore undiscovered viruses. Virus databases play a crucial role in providing access to sequences, annotations and other metadata, and analysis tools for studying viruses. However, there has not been a comprehensive review of virus databases in the last five years. This study aimed to fill this gap by identifying 24 active virus databases and included an extensive evaluation of their content, functionality and compliance with the FAIR principles. In this study, we thoroughly assessed the search capabilities of five database catalogs, which serve as comprehensive repositories housing a diverse array of databases and offering essential metadata. Moreover, we conducted a comprehensive review of different types of errors, encompassing taxonomy, names, missing information, sequences, sequence orientation, and chimeric sequences, with the intention of empowering users to effectively tackle these challenges. We expect this review to aid users in selecting suitable virus databases and other resources, and to help databases in error management and improve their adherence to the FAIR principles. The databases listed here represent the current knowledge of viruses and will help aid users find databases of interest based on content, functionality, and scope. The use of virus databases is integral to gaining new insights into the biology, evolution, and transmission of viruses, and developing new strategies to manage virus outbreaks and preserve global health.

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

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

The universal Sua5/TsaC family evolved different mechanisms for the synthesis of a key tRNA modification

TsaC/Sua5 family of enzymes catalyzes the first step in the synthesis of N6-threonyl-carbamoyl adenosine (t⁶A) one of few truly ubiquitous tRNA modifications important for translation accuracy. TsaC is a single domain protein while Sua5 proteins contains a TsaC-like domain and an additional SUA5 domain of unknown function. The emergence of these two proteins and their respective mechanisms for t⁶A synthesis remain poorly understood. Here, we performed phylogenetic and comparative sequence and structure analysis of TsaC and Sua5 proteins. We confirm that this family is ubiquitous but the co-occurrence of both variants in the same organism is rare and unstable. We further find that obligate symbionts are the only organisms lacking sua5 or tsaC genes. The data suggest that Sua5 was the ancestral version of the enzyme while TsaC arose via loss of the SUA5 domain that occurred multiple times in course of evolution. Multiple losses of one of the two variants in combination with horizontal gene transfers along a large range of phylogenetic distances explains the present day patchy distribution of Sua5 and TsaC. The loss of the SUA5 domain triggered adaptive mutations affecting the substrate binding in TsaC proteins. Finally, we identified atypical Sua5 proteins in Archaeoglobi archaea that seem to be in the process of losing the SUA5 domain through progressive gene erosion. Together, our study uncovers the evolutionary path for emergence of these homologous isofunctional enzymes and lays the groundwork for future experimental studies on the function of TsaC/Sua5 proteins in maintaining faithful translation.

A complete telomere-to-telomere assembly of the maize genome

A complete telomere-to-telomere (T2T) finished genome has been the long pursuit of genomic research. Through generating deep coverage ultralong Oxford Nanopore Technology (ONT) and PacBio HiFi reads, we report here a complete genome assembly of maize with each chromosome entirely traversed in a single contig. The 2,178.6 Mb T2T Mo17 genome with a base accuracy of over 99.99% unveiled the structural features of all repetitive regions of the genome. There were several super-long simple-sequence-repeat arrays having consecutive thymine-adenine-guanine (TAG) tri-nucleotide repeats up to 235 kb. The assembly of the entire nucleolar organizer region of the 26.8 Mb array with 2,974 45S rDNA copies revealed the enormously complex patterns of rDNA duplications and transposon insertions. Additionally, complete assemblies of all ten centromeres enabled us to precisely dissect the repeat compositions of both CentC-rich and CentC-poor centromeres. The complete Mo17 genome represents a major step forward in understanding the complexity of the highly recalcitrant repetitive regions of higher plant genomes.

Towards the accurate alignment of over a million protein sequences: Current state of the art

Large-scale genomics requires highly scalable and accurate multiple sequence alignment methods. Results collected over this last decade suggest accuracy loss when scaling up over a few thousand sequences. This issue has been actively addressed with a number of innovative algorithmic solutions that combine low-level hardware optimization with novel higher-level heuristics. This review provides an extensive critical overview of these recent methods. Using established reference datasets we conclude that albeit significant progress has been achieved, a unified framework able to consistently and efficiently produce high-accuracy large-scale multiple alignments is still lacking.

Evolution of piggyBac Transposons in Apoidea

In this study, we investigated the presence of piggyBac (PB) transposons in 44 bee genomes from the Apoidea order, which is a superfamily within the Hymenoptera, which includes a large number of bee species crucial for pollination. We annotated the PB transposons in these 44 bee genomes and examined their evolution profiles, including structural characteristics, distribution, diversity, activity, and abundance. The mined PB transposons were divided into three clades, with uneven distribution in each genus of PB transposons in Apoidea. The complete PB transposons we discovered are around 2.23-3.52 kb in length and encode transposases of approximately 580 aa, with terminal inverted repeats (TIRs) of about 14 bp and 4 bp (TTAA) target-site duplications. Long TIRs (200 bp, 201 bp, and 493 bp) were also detected in some species of bees. The DDD domains of the three transposon types were more conserved, while the other protein domains were less conserved. Generally, most PB transposons showed low abundance in the genomes of Apoidea. Divergent evolution dynamics of PB were observed in the genomes of Apoidea. PB transposons in some identified species were relatively young, whiles others were older and with some either active or inactive. In addition, multiple invasions of PB were also detected in some genomes of Apoidea. Our findings highlight the contribution of PB transposons to genomic variation in these species and suggest their potential as candidates for future gene transfer tools.

Genomic basis of Y-linked dwarfism in cichlids pursuing alternative reproductive tactics

Sexually antagonistic selection, which favours different optima in males and females, is predicted to play an important role in the evolution of sex chromosomes. Body size is a sexually antagonistic trait in the shell-brooding cichlid fish Lamprologous callipterus, as "bourgeois" males must be large enough to carry empty snail shells to build nests whereas females must be small enough to fit into shells for breeding. In this species, there is also a second male morph: smaller "dwarf" males employ an alternative reproductive strategy by wriggling past spawning females into shells to fertilize eggs. L. callipterus male morphology is passed strictly from father to son, suggesting Y-linkage. However, sex chromosomes had not been previously identified in this species, and the genomic basis of size dimorphism was unknown. Here we used whole-genome sequencing to identify a 2.4-Mb sex-linked region on scaffold_23 with reduced coverage and single nucleotide polymorphism density in both male morphs compared to females. Within this sex region, distinct Y-haplotypes delineate the two male morphs, and candidate genes for body size (GHRHR, a known dwarfism gene) and sex determination (ADCYAP1R1) are in high linkage disequilibrium. Because differences in body size between females and males are under strong selection in L. callipterus, we hypothesize that sexual antagonism over body size initiated early events in sex chromosome evolution, followed by Y divergence to give rise to bourgeois and dwarf male reproductive strategies. Our results are consistent with the hypothesis that sexually antagonistic traits should be linked to young sex chromosomes.

Evidence for a Putative Isoprene Reductase in Acetobacterium wieringae

Recent discoveries of isoprene-metabolizing microorganisms suggest they might play an important role in the global isoprene budget. Under anoxic conditions, isoprene can be used as an electron acceptor and is reduced to methylbutene. This study describes the proteogenomic profiling of an isoprene-reducing bacterial culture to identify organisms and genes responsible for the isoprene hydrogenation reaction. A metagenome-assembled genome (MAG) of the most abundant (89% relative abundance) lineage in the enrichment, Acetobacterium wieringae, was obtained. Comparative proteogenomics and reverse transcription-PCR (RT-PCR) identified a putative five-gene operon from the A. wieringae MAG upregulated during isoprene reduction. The operon encodes a putative oxidoreductase, three pleiotropic nickel chaperones (2 × HypA, HypB), and one 4Fe-4S ferredoxin. The oxidoreductase is proposed as the putative isoprene reductase with a binding site for NADH, flavin adenine dinucleotide (FAD), two pairs of canonical [4Fe-4S] clusters, and a putative iron-sulfur cluster site in a Cys₆-bonding environment. Well-studied Acetobacterium strains, such as A. woodii DSM 1030, A. wieringae DSM 1911, or A. malicum DSM 4132, do not encode the isoprene-regulated operon but encode, like many other bacteria, a homolog of the putative isoprene reductase (~47 to 49% amino acid sequence identity). Uncharacterized homologs of the putative isoprene reductase are observed across the Firmicutes, Spirochaetes, Tenericutes, Actinobacteria, Chloroflexi, Bacteroidetes, and Proteobacteria, suggesting the ability of biohydrogenation of unfunctionalized conjugated doubled bonds in other unsaturated hydrocarbons. IMPORTANCE Isoprene was recently shown to act as an electron acceptor for a homoacetogenic bacterium. The focus of this study is the molecular basis for isoprene reduction. By comparing a genome from our isoprene-reducing enrichment culture, dominated by Acetobacterium wieringae, with genomes of other Acetobacterium lineages that do not reduce isoprene, we shortlisted candidate genes for isoprene reduction. Using comparative proteogenomics and reverse transcription-PCR we have identified a putative five-gene operon encoding an oxidoreductase referred to as putative isoprene reductase.

Intermittent antibiotic treatment of bacterial biofilms favors the rapid evolution of resistance

Bacterial antibiotic resistance is a global health concern of increasing importance and intensive study. Although biofilms are a common source of infections in clinical settings, little is known about the development of antibiotic resistance within biofilms. Here, we use experimental evolution to compare selection of resistance mutations in planktonic and biofilm Escherichia coli populations exposed to clinically relevant cycles of lethal treatment with the aminoglycoside amikacin. Consistently, mutations in sbmA, encoding an inner membrane peptide transporter, and fusA, encoding the essential elongation factor G, are rapidly selected in biofilms, but not in planktonic cells. This is due to a combination of enhanced mutation rate, increased adhesion capacity and protective biofilm-associated tolerance. These results show that the biofilm environment favors rapid evolution of resistance and provide new insights into the dynamic evolution of antibiotic resistance in biofilms.

Genome Resources for Four Clarireedia Species Causing Dollar Spot on Diverse Turfgrasses

Dollar spot (DS) is a destructive fungal disease impacting almost all warm- and cool-season turfgrasses worldwide. Multiple fungal species in the genus Clarireedia are causal agents of DS. Here, we present whole-genome assemblies of nine fungal isolates in the genus Clarireedia, including four species (C. paspali, C. hainanense, C. jacksonii, and C. monteithiana) causing DS on seashore paspalum (Paspalum vaginatum Sw.), creeping bentgrass (Agrostis stolonifera L.), and Kentucky bluegrass (Poa pratensis L.) in China. This work provides valuable baseline genomic data to support further research and management of DS pathogens on turfgrasses.

Passer, a highly active transposon from a fish genome, as a potential new robust genetic manipulation tool

The discovery of new, active DNA transposons can expand the range of genetic tools and provide more options for genomic manipulation. In this study, a bioinformatics analysis suggested that Passer (PS) transposons, which are members of the pogo superfamily, show signs of recent and current activity in animals and may be active in some species. Cell-based transposition assays revealed that the native PS transposases from Gasterosteus aculeatus and Danio rerio displayed very high activity in human cells relative to the Sleeping Beauty transposon. A typical overproduction inhibition phenomenon was observed for PS, and transposition capacity was decreased by ∼12% with each kilobase increase in the insertion size. Furthermore, PS exhibited a pronounced integration preference for genes and their transcriptional regulatory regions. We further show that two domesticated human proteins derived from PS transposases have lost their transposition activity. Overall, PS may represent an alternative with a potentially efficient genetic manipulation tool for transgenesis and mutagenesis applications.