Rfam 13.0: shifting to a genome-centric resource for non-coding RNA families

The chromosome-level genome assembly of the red swamp crayfish Procambarus clarkii

Red swamp crayfish, Procambarus clarkii, is the most cultured freshwater crayfish species. It attracts significant research attention due to its considerable economic importance. However, the limited availability of genome information has impeded further genetic studies and breeding programs. By utilizing Illumina, PacBio, and Hi-C sequencing technologies, we present a more comprehensive and continuous chromosome-level assembly for P. clarkii than the published one. The final genome size is 4.03 Gb, consisting of 2,358 scaffolds with a N50 of 42.87 Mb. Notably, 3.68 Gb, corresponding to 91.42% of the genome, was anchored to 94 chromosomes. The assembly comprises 70.64% repetitive sequences, including 5.21% tandem repeats and 65.40% transposable elements. Additionally, a total of 4,456 non-coding RNAs and 28,852 protein-coding genes were predicted in the P. clarkii genome, with 96.26% of the genes were annotated. This high-quality genome assembly not only represents a significant improvement for the genome of P. clarkii and provides insights into the unique genome evolution, but also offers valuable information for developing freshwater aquaculture and accelerating genetic breeding.

Genomic variation, environmental adaptation, and feralization in ramie, an ancient fiber crop

Feralization is an important evolutionary process, but the mechanisms behind it remain poorly understood. Here, we use the ancient fiber crop ramie (Boehmeria nivea (L.) Gaudich.) as a model to investigate genomic changes associated with both domestication and feralization. We first produced a chromosome-scale de novo genome assembly of feral ramie and investigated structural variations between feral and domesticated ramie genomes. Next, we gathered 915 accessions from 23 countries, comprising cultivars, major landraces, feral populations, and the wild progenitor. Based on whole-genome resequencing of these accessions, we constructed the most comprehensive ramie genomic variation map to date. Phylogenetic, demographic, and admixture signal detection analyses indicated that feral ramie is of exoferal or exo-endo origin, i.e., descended from hybridization between domesticated ramie and the wild progenitor or ancient landraces. Feral ramie has higher genetic diversity than wild or domesticated ramie, and genomic regions affected by natural selection during feralization differ from those under selection during domestication. Ecological analyses showed that feral and domesticated ramie have similar ecological niches that differ substantially from the niche of the wild progenitor, and three environmental variables are associated with habitat-specific adaptation in feral ramie. These findings advance our understanding of feralization, providing a scientific basis for the excavation of new crop germplasm resources and offering novel insights into the evolution of feralization in nature.

Characteristics of Corynespora cassiicola, the causal agent of tobacco Corynespora leaf spot, revealed by genomic and metabolic phenomic analysis

Corynespora cassiicola is a highly diverse fungal pathogen that commonly occurs in tropical, subtropical, and greenhouse environments worldwide. In this study, the isolates were identified as C. cassiicola, and the optimum growth and sporulation were studied. The phenotypic characteristics of C. cassiicola, concerning 950 different growth conditions, were tested using Biolog PM plates 1-10. In addition, the strain of C. cassiicola DWZ from tobacco hosts was sequenced for the using Illumina PE150 and Pacbio technologies. The host resistance of tobacco Yunyan 87 with different maturity levels was investigated. In addition, the resistance evaluation of 10 common tobacco varieties was investigated. The results showed that C. cassiicola metabolized 89.47% of the tested carbon source, 100% of the nitrogen source, 100% of the phosphorus source, and 97.14% of the sulfur source. It can adapt to a variety of different osmotic pressure and pH environments, and has good decarboxylase and deaminase activities. The optimum conditions for pathogen growth and sporulation were 25-30 °C, and the growth was better on AEA and OA medium. The total length of the genome was 45.9 Mbp, the GC content was 51.23%, and a total of 13,061 protein-coding genes, 202 non-coding RNAs and 2801 and repeat sequences were predicted. Mature leaves were more susceptible than proper mature and immature leaves, and the average diameter of diseased spots reached 17.74 mm at 12 days. None of the tested ten cultivars exhibited obvious resistance to Corynespora leaf spot of tobacco, whereby all disease spot diameters reached > 10 mm and > 30 mm when at 5 and 10 days after inoculation, respectively. The phenotypic characteristics, genomic analysis of C. cassiicola and the cultivar resistance assessment of this pathogen have increased our understanding of Corynespora leaf spot of tobacco.

Assessing a respiratory toxic infectious bronchitis virus (IBV) strain: isolation, identification, pathogenicity, and immunological failure insights

Infectious bronchitis virus (IBV) is caused by avian coronavirus and poses a global economic threat to the poultry industry. In 2023, a highly pathogenic IBV strain, IBV/CN/GD20230501, was isolated and identified from chickens vaccinated with IBV-M41 in Guangdong, China. This study comprehensively investigated the biological characteristics of the isolated IBV strain, including its genotype, whole genome sequence analysis of its S1 gene, pathogenicity, host immune response, and serum non-targeted metabolomics. Through the analysis of the S1 gene sequence, serum neutralization tests, and comparative genomics, it was proven that IBV/CN/GD20230501 belongs to the GI-I type of strain and is serotype II. One alanine residue in the S1 subunit of the isolated strain was mutated into serine, and some mutations were observed in the ORF1ab gene and the terminal region of the genome. Animal challenge experiments using the EID50 and TCID50 calculations showed that IBV/CN/GD20230501 possesses strong respiratory pathogenicity, with early and long-term shedding of viruses and rapid viral spread. Antibody detection indicated that chickens infected with IBV/CN/GD20230501 exhibited delayed expression of early innate immune genes, while those infected with M41 showed rapid gene induction and effective viral control. Metabolomics analysis demonstrated that this virus infection led to differential expression of 291 ions in chicken serum, mainly affecting the citric acid cycle (tricarboxylic acid cycle).IMPORTANCEThis study identified an infectious bronchitis virus (IBV) strain isolated from vaccinated chickens in an immunized population that had certain sequence differences compared to IBV-M41, resulting in significantly enhanced pathogenicity and host defense. This strain has the potential to replace M41 as a more suitable challenge model for drug research. The non-targeted metabolomics analysis highlighting the citric acid cycle provides a new avenue for studying this highly virulent strain.

Comprehensive analysis of the differential expression of mRNAs, lncRNAs, and miRNAs in Zi goose testis with high and low sperm mobility

Sperm mobility (SM) is an objective index for measuring sperm motility; however, the mechanisms underlying its regulation in geese remain unclear. The present study sought to elucidate the genetic mechanism underlying SM traits in Zi geese (Anser cygnoides L.). To this end, three successive experiments were performed. In Experiment I, SM was determined in 40 ganders; the 3 ganders with the highest mobility and three with the lowest mobility were assigned to the high and low sperm mobility rank (SMR) groups, respectively. In Experiment II, the differences in fertility between the two SMR groups were assessed within two breeding flocks comprising the selected six ganders from Experiment I and 30 females (each flock had 3 ganders and 15 females). In Experiment III, the testes of the 6 ganders were harvested for histological observation and whole-transcriptome sequencing. Results revealed better fertility, well-developed seminiferous tubules, and abundant mature sperm in the high-SMR-flock compared to those of the low-SMR-flock (89 vs. 81%) (P < 0.05). Differential expression (DE) analysis identified 76 mRNAs, 344 lncRNAs, and 17 miRNAs between the SMR groups, with LOC106049708, XPNPEP3, GNB3, ADCY8, PRKAG3, oha-miR-182-5p, and ocu-miR-10b-5p identified as key mRNAs and miRNAs contributing to SM. Enrichment analysis implicated these DE RNAs in pathways related to ATP binding, cell metabolism, apelin signaling, Wnt signaling, and Adherens junctions. Additionally, competing endogenous RNA (ceRNA) networks comprising 9 DE mRNAs, 17 DE miRNAs, and 169 DE lncRNAs were constructed. Two ceRNA network pathways (LOC106049708-oha-miR-182-5p-MSTRG.2479.6 and PRKAG3-ocu-miR-10b-5p-MSTRG.9047.14) were identified as key regulators of SM in geese. These findings offer crucial insights into the identification of key genes and ceRNA pathways influencing sperm mobility in geese.

Comprehensive analysis of insertion sequences within rRNA genes of CPR bacteria and biochemical characterization of a homing endonuclease encoded by these sequences

The Candidate Phyla Radiation (CPR) represents an extensive bacterial clade comprising primarily uncultured lineages and is distinguished from other bacteria by a significant prevalence of insertion sequences (ISs) within their rRNA genes. However, our understanding of the taxonomic distribution and characteristics of these ISs remains limited. In this study, we used a comprehensive approach to systematically determine the nature of the rRNA ISs in CPR bacteria. The analysis of hundreds of rRNA gene sequences across 65 CPR phyla revealed that ISs are present in 48% of 16S rRNA genes and 82% of 23S rRNA genes, indicating a broad distribution across the CPR clade, with exceptions in the 16S and 23S rRNA genes of Candidatus (Ca.) Saccharibacteria and the 16S rRNA genes of Ca. Peregrinibacteria. Over half the ISs display a group-I-intron-like structure, whereas specific 16S rRNA gene ISs display features reminiscent of group II introns. The ISs frequently encode proteins with homing endonuclease (HE) domains, centered around the LAGLIDADG motif. The LAGLIDADG HE (LHE) proteins encoded by the rRNA ISs of CPR bacteria predominantly have a single-domain structure, deviating from the usual single- or double-domain configuration observed in typical prokaryotic LHEs. Experimental analysis of one LHE protein, I-ShaI from Ca. Shapirobacteria, confirmed that its endonuclease activity targets the DNA sequence of its insertion site, and chemical cross-linking experiments demonstrated its capacity to form homodimers. These results provide robust evidence supporting the hypothesis that the explosive proliferation of rRNA ISs in CPR bacteria was facilitated by mechanisms involving LHEs.

IMPORTANCE

Insertion sequences (ISs) in rRNA genes are relatively limited and infrequent in most bacterial phyla. With a comprehensive bioinformatic analysis, we show that in CPR bacteria, these ISs occur in 48% of 16S rRNA genes and 82% of 23S rRNA genes. We also report the systematic and biochemical characterization of the LAGLIDADG homing endonucleases (LHEs) encoded by these ISs in the first such analysis of the CPR bacteria. This study significantly extends our understanding of the phylogenetic positions of rRNA ISs within CPR bacteria and the biochemical features of their LHEs.

Halotolerant Bacillus sp. strain RA coordinates myo-inositol metabolism to confer salt tolerance to tomato

Soil salinity is a worldwide problem threatening crop yields. Some plant growth-promoting rhizobacteria (PGPR) could survive in high salt environment and assist plant adaptation to stress. Nevertheless, the genomic and metabolic features, as well as the regulatory mechanisms promoting salt tolerance in plants by these bacteria remain largely unknown. In the current work, a novel halotolerant PGPR strain, namely, Bacillus sp. strain RA can enhance tomato tolerance to salt stress. Comparative genomic analysis of strain RA with its closely related species indicated a high level of evolutionary plasticity exhibited by strain-specific genes and evolutionary constraints driven by purifying selection, which facilitated its genomic adaptation to salt-affected soils. The transcriptome further showed that strain RA could tolerate salt stress by balancing energy metabolism via the reprogramming of biosynthetic pathways. Plants exude a plethora of metabolites that can strongly influence plant fitness. The accumulation of myo-inositol in leaves under salt stress was observed, leading to the promotion of plant growth triggered by Bacillus sp. strain RA. Importantly, myo-inositol serves as a selective force in the assembly of the phyllosphere microbiome and the recruitment of plant-beneficial species. It promotes destabilizing properties in phyllosphere bacterial co-occurrence networks, but not in fungal networks. Furthermore, interdomain interactions between bacteria and fungi were strengthened by myo-inositol in response to salt stress. This work highlights the genetic adaptation of RA to salt-affected soils and its ability to impact phyllosphere microorganisms through the adjustment of myo-inositol metabolites, thereby imparting enduring resistance against salt stress in tomato.

Promethearchaeum syntrophicum gen. nov., sp. nov., an anaerobic, obligately syntrophic archaeon, the first isolate of the lineage 'Asgard' archaea, and proposal of the new archaeal phylum Promethearchaeota phyl. nov. and kingdom Promethearchaeati regn. nov

An anaerobic, mesophilic, syntrophic, archaeon strain MK-D1T, was isolated as a pure co-culture with Methanogenium sp. strain MK-MG from deep-sea methane seep sediment. This organism is, to our knowledge, the first cultured representative of 'Asgard' archaea, an archaeal group closely related to eukaryotes. Here, we describe the detailed physiology and phylogeny of MK-D1T and propose Promethearchaeum syntrophicum gen. nov., sp. nov. to accommodate this strain. Cells were non-motile, small cocci, approximately 300-750 nm in diameter and produced membrane vesicles, chains of blebs and membrane-based protrusions. MK-D1T grew at 4-30 °C with optimum growth at 20 °C. The strain grew chemoorganotrophically with amino acids, peptides and yeast extract with obligate dependence on syntrophy with H2-/formate-utilizing organisms. MK-D1T showed the fastest growth and highest maximum cell yield when grown with yeast extract as the substrate: approximately 3 months to full growth, reaching up to 6.7×106 16S rRNA gene copies ml-1. MK-D1T had a circular 4.32 Mb chromosome with a DNA G+C content of 31.1 mol%. The results of phylogenetic analyses of the 16S rRNA gene and conserved marker proteins indicated that the strain is affiliated with 'Asgard' archaea and more specifically DHVC1/DSAG/MBG-B and 'Lokiarchaeota'/'Lokiarchaeia'. On the basis of the results of 16S rRNA gene sequence analysis, the most closely related isolated relatives were Infirmifilum lucidum 3507LTT (76.09 %) and Methanothermobacter tenebrarum RMAST (77.45 %) and the closest relative in enrichment culture was Candidatus 'Lokiarchaeum ossiferum' (95.39 %). The type strain of the type species is MK-D1T (JCM 39240T and JAMSTEC no. 115508). We propose the associated family, order, class, phylum, and kingdom as Promethearchaeaceae fam. nov., Promethearchaeales ord. nov., Promethearchaeia class. nov., Promethearchaeota phyl. nov., and Promethearchaeati regn. nov., respectively. These are in accordance with ICNP Rules 8 and 22 for nomenclature, Rule 30(3)(b) for validation and maintenance of the type strain, and Rule 31a for description as a member of an unambiguous syntrophic association.

Highly efficient production of 2-phenylethanol by wild-type Saccharomyces bayanus strain

2-Phenylethanol (2-PE) is a highly valuable aromatic alcohol utilized in fragrance, cosmetics and food industries. Due to the toxic by-products from chemical synthesis and the low productivity of the extraction method, bioproduction of 2-PE by yeast is considered promising. In this study, a wild-type Saccharomyces bayanus L1 strain producing 2-PE was isolated from soy sauce mash. Transcriptional analysis showed that 2-PE was synthesized via the Ehrlich pathway and Shikimate pathway in S. bayanus L1. By improving the fermentation conditions in shaking flasks, the maximum 2-PE titer reached 4.2 g/L with a productivity of 0.058 g/L/h within 72 h. In fed-batch fermentation, S. bayanus L1 strain produced 6.5 g/L of 2-PE within 60 h, achieving a productivity of 0.108 g/L/h. These findings suggest that S. bayanus L1 strain is an efficient 2-PE producer, paving the way for highly efficient 2-PE production.

Genome sequencing, comparative analysis, and gene expression responses of cytochrome P450 genes in Oryzias curvinotus provide insights into environmental adaptation

The mangrove fish (Oryzias curvinotus) serves as a model for researching environmental adaptation and sexual development. To further such research, we sequenced and assembled a high-quality 842 Mb reference genome for O. curvinotus. Comparative genomic analysis revealed 891 expanded gene families, including significantly expanded cytochrome P450 (CYP) detoxification genes known to be involved in xenobiotic defense. We identified 69 O. curvinotus CYPs (OcuCYPs) across 18 families and 10 clans using multiple methods. Extensive RNA-seq and qPCR analysis demonstrated diverse spatiotemporal expression patterns of OcuCYPs by developmental stage, tissue type, sex, and pollutant exposure (17β-estradiol (E2) and testosterone (MT)). Many OcuCYPs exhibited sexual dimorphism in gonads, suggesting reproductive roles in steroidogenesis, while their responsiveness to model toxicants indicates their importance in environmental adaptation through enhanced detoxification. Pathway analysis highlighted expanded CYP genes in arachidonic acid metabolism, drug metabolism, and steroid hormone biosynthesis. This chromosome-level genomic resource provides crucial biological insights to elucidate the functional roles of expanded CYPs in environmental adaptation, sexual development, early life history, and conservation in the anthropogenically impacted mangrove habitats of O. curvinotus. It also enables future ecotoxicology research leveraging O. curvinotus as a pollution sentinel species.

Changes in microRNA expression related to ischemia-reperfusion injury in the kidney of the thirteen-lined ground squirrel during torpor

During the hibernation season, the thirteen-lined ground squirrel undergoes cyclical torpor and arousal periods. The decrease and restoration of metabolic rate and oxygen delivery during torpor and arousal, respectively, may cause reperfusion-ischemia injury in the kidneys. In order to maintain the structural integrity of the kidneys necessary for renal function resumption during arousal, the thirteen-lined ground squirrel has developed adaptive methods to prevent and repair kidney injury. In this present study, computational methods were used to clean and analyze sequenced kidney RNA samples. Significantly differentially expressed microRNAs and enriched gene sets were also determined. From the gene set analysis, the results showed an increase in ubiquitin-related processes and p53 signaling pathways which suggested the occurrence of kidney damage during torpor. There was also an observed increase in cell cycle processes and the anchoring junction cellular compartment which may lend to the prevention of kidney injury. From the differentially expressed microRNAs, miR-27a (log2FC = 1.639; p-value = 0.023), miR-129 (log2FC = 2.516; p-value = 0.023), miR-let-7b (log2FC = 2.360; p-value = 0.025), miR-let-7c (log2FC = 2.291; p-value = 0.037) and miR-let-7i (log2FC = 1.564; p-value = 0.039) were found to be significantly upregulated. These biochemical adaptations may allow the thirteen-lined ground squirrel to maintain kidney structure and function during hibernation.

Phylogenomics and genetic analysis of solvent-producing Clostridium species

The genus Clostridium is a large and diverse group within the Bacillota (formerly Firmicutes), whose members can encode useful complex traits such as solvent production, gas-fermentation, and lignocellulose breakdown. We describe 270 genome sequences of solventogenic clostridia from a comprehensive industrial strain collection assembled by Professor David Jones that includes 194 C. beijerinckii, 57 C. saccharobutylicum, 4 C. saccharoperbutylacetonicum, 5 C. butyricum, 7 C. acetobutylicum, and 3 C. tetanomorphum genomes. We report methods, analyses and characterization for phylogeny, key attributes, core biosynthetic genes, secondary metabolites, plasmids, prophage/CRISPR diversity, cellulosomes and quorum sensing for the 6 species. The expanded genomic data described here will facilitate engineering of solvent-producing clostridia as well as non-model microorganisms with innately desirable traits. Sequences could be applied in conventional platform biocatalysts such as yeast or Escherichia coli for enhanced chemical production. Recently, gene sequences from this collection were used to engineer Clostridium autoethanogenum, a gas-fermenting autotrophic acetogen, for continuous acetone or isopropanol production, as well as butanol, butanoic acid, hexanol and hexanoic acid production.

Coping with alpine habitats: genomic insights into the adaptation strategies of Triplostegia glandulifera (Caprifoliaceae)

How plants find a way to thrive in alpine habitats remains largely unknown. Here we present a chromosome-level genome assembly for an alpine medicinal herb, Triplostegia glandulifera (Caprifoliaceae), and 13 transcriptomes from other species of Dipsacales. We detected a whole-genome duplication event in T. glandulifera that occurred prior to the diversification of Dipsacales. Preferential gene retention after whole-genome duplication was found to contribute to increasing cold-related genes in T. glandulifera. A series of genes putatively associated with alpine adaptation (e.g. CBFs, ERF-VIIs, and RAD51C) exhibited higher expression levels in T. glandulifera than in its low-elevation relative, Lonicera japonica. Comparative genomic analysis among five pairs of high- vs low-elevation species, including a comparison of T. glandulifera and L. japonica, indicated that the gene families related to disease resistance experienced a significantly convergent contraction in alpine plants compared with their lowland relatives. The reduction in gene repertory size was largely concentrated in clades of genes for pathogen recognition (e.g. CNLs, prRLPs, and XII RLKs), while the clades for signal transduction and development remained nearly unchanged. This finding reflects an energy-saving strategy for survival in hostile alpine areas, where there is a tradeoff with less challenge from pathogens and limited resources for growth. We also identified candidate genes for alpine adaptation (e.g. RAD1, DMC1, and MSH3) that were under convergent positive selection or that exhibited a convergent acceleration in evolutionary rate in the investigated alpine plants. Overall, our study provides novel insights into the high-elevation adaptation strategies of this and other alpine plants.

Genomic analysis reveals phylogeny of Zygophyllales and mechanism for water retention of a succulent xerophyte

Revealing the genetic basis for stress-resistant traits in extremophile plants will yield important information for crop improvement. Zygophyllum xanthoxylum, an extant species of the ancient Mediterranean, is a succulent xerophyte that can maintain a favorable water status under desert habitats; however, the genetic basis of this adaptive trait is poorly understood. Furthermore, the phylogenetic position of Zygophyllales, to which Z. xanthoxylum belongs, remains controversial. In this study, we sequenced and assembled the chromosome-level genome of Z. xanthoxylum. Phylogenetic analysis showed that Zygophyllales and Myrtales form a separated taxon as a sister to the clade comprising fabids and malvids, clarifying the phylogenetic position of Zygophyllales at whole-genome scale. Analysis of genomic and transcriptomic data revealed multiple critical mechanisms underlying the efficient osmotic adjustment using Na+ and K+ as "cheap" osmolytes that Z. xanthoxylum has evolved through the expansion and synchronized expression of genes encoding key transporters/channels and their regulators involved in Na+/K+ uptake, transport, and compartmentation. It is worth noting that ZxCNGC1;1 (cyclic nucleotide-gated channels) and ZxCNGC1;2 constituted a previously undiscovered energy-saving pathway for Na+ uptake. Meanwhile, the core genes involved in biosynthesis of cuticular wax also featured an expansion and upregulated expression, contributing to the water retention capacity of Z. xanthoxylum under desert environments. Overall, these findings boost the understanding of evolutionary relationships of eudicots, illustrate the unique water retention mechanism in the succulent xerophyte that is distinct from glycophyte, and thus provide valuable genetic resources for the improvement of stress tolerance in crops and insights into the remediation of sodic lands.

The genome of the colonial hydroid Hydractinia reveals that their stem cells use a toolkit of evolutionarily shared genes with all animals

Hydractinia is a colonial marine hydroid that shows remarkable biological properties, including the capacity to regenerate its entire body throughout its lifetime, a process made possible by its adult migratory stem cells, known as i-cells. Here, we provide an in-depth characterization of the genomic structure and gene content of two Hydractinia species, Hydractinia symbiolongicarpus and Hydractinia echinata, placing them in a comparative evolutionary framework with other cnidarian genomes. We also generated and annotated a single-cell transcriptomic atlas for adult male H. symbiolongicarpus and identified cell-type markers for all major cell types, including key i-cell markers. Orthology analyses based on the markers revealed that Hydractinia's i-cells are highly enriched in genes that are widely shared amongst animals, a striking finding given that Hydractinia has a higher proportion of phylum-specific genes than any of the other 41 animals in our orthology analysis. These results indicate that Hydractinia's stem cells and early progenitor cells may use a toolkit shared with all animals, making it a promising model organism for future exploration of stem cell biology and regenerative medicine. The genomic and transcriptomic resources for Hydractinia presented here will enable further studies of their regenerative capacity, colonial morphology, and ability to distinguish self from nonself.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Genomic analysis and experimental pathogenic characterization of Riemerella anatipestifer isolates from chickens in China

Waterfowl have a high likelihood of being infected with Riemerella anatipestifer. Although the pathogen is found in domestic ducks, turkeys, geese, and wild birds, there is little information available about the consequences of infection during egg laying and hatching in chickens. Here, we present the first report of a novel sequence type of R. anatipestifer S63 isolated from chickens in China. On the basis of pan-genome analysis, we showed S63's genome occupies a distinct branch with other R. anatipestifer isolates from other hosts. Galleria mellonella larval tests indicated that S63 is less virulent than R. anatipestifer Ra36 isolated from ducks. Ducks and hens are susceptible to S63 infection. There is no mortality rate for chickens or ducks, but adult chickens experience neurological symptoms that reduce egg production and hatching rates. In chickens, S63 might be passed vertically from parents to offspring, resulting in "jelly-like" lifeless embryos. Using quantitative PCR, S63 was detected in the brain, liver, reproductive organs, and embryos. As far as we know, this is the first report of R. anatipestifer in hens, a disease that can reduce egg productivity, lower hatching rates, and produce jelly-like lifeless embryos, and the first report to raise the possibility that hens can be infected by roosters via semen.

Mating system is associated with seed phenotypes upon loss of RNA-directed DNA methylation in Brassicaceae

In plants, de novo DNA methylation is guided by 24-nt short interfering (si)RNAs in a process called RNA-directed DNA methylation (RdDM). Primarily targeted at transposons, RdDM causes transcriptional silencing and can indirectly influence expression of neighboring genes. During reproduction, a small number of siRNA loci are dramatically upregulated in the maternally derived seed coat, suggesting that RdDM might have a special function during reproduction. However, the developmental consequence of RdDM has been difficult to dissect because disruption of RdDM does not result in overt phenotypes in Arabidopsis (Arabidopsis thaliana), where the pathway has been most thoroughly studied. In contrast, Brassica rapa mutants lacking RdDM have a severe seed production defect, which is determined by the maternal sporophytic genotype. To explore the factors that underlie the different phenotypes of these species, we produced RdDM mutations in 3 additional members of the Brassicaceae family: Camelina sativa, Capsella rubella, and Capsella grandiflora. Among these 3 species, only mutations in the obligate outcrosser, C. grandiflora, displayed a seed production defect similar to Brassica rapa mutants, suggesting that mating system is a key determinant for reproductive phenotypes in RdDM mutants.

A high-quality chromosome-level genome assembly of the Chinese medaka Oryzias sinensis

Oryzias sinensis, also known as Chinese medaka or Chinese ricefish, is a commonly used animal model for aquatic environmental assessment in the wild as well as gene function validation or toxicology research in the lab. Here, a high-quality chromosome-level genome assembly of O. sinensis was generated using single-tube long fragment read (stLFR) reads, Nanopore long-reads, and Hi-C sequencing data. The genome is 796.58 Mb, and a total of 712.17 Mb of the assembled sequences were anchored to 23 pseudo-chromosomes. A final set of 22,461 genes were annotated, with 98.67% being functionally annotated. The Benchmarking Universal Single-Copy Orthologs (BUSCO) benchmark of genome assembly and gene annotation reached 95.1% (93.3% single-copy) and 94.6% (91.7% single-copy), respectively. Furthermore, we also use ATAC-seq to uncover chromosome transposase-accessibility as well as related genome area function enrichment for Oryzias sinensis. This study offers a new improved foundation for future genomics research in Chinese medaka.

Perception of the Biocontrol Potential and Palmitic Acid Biosynthesis Pathway of Bacillus subtilis H2 through Merging Genome Mining with Chemical Analysis

Bacillus has been widely studied for its potential to protect plants from pathogens. Here, we report the whole genome sequence of Bacillus subtilis H2, which was isolated from the tea garden soil of Guiyang Forest Park. Strain H2 showed a broad spectrum of antagonistic activities against many plant fungal pathogens and bacteria pathogens, including the rice blast fungus Magnaporthe oryzae, and showed a good field control effect against rice blast. The complete genome of B. subtilis H2 contained a 4,160,635-bp circular chromosome, with an average G + C content of 43.78%. Through the genome mining of strain H2, we identified 7 known antimicrobial compound biosynthetic gene clusters (BGCs) including sporulation killing factor, surfactin, bacillaene, fengycin, bacillibactin, subtilosin A, and bacilysin. Palmitic acid (PA), a secondary metabolite, was detected and identified in the H2 strain through genome mining analysis and gas chromatography-mass spectrometry (GC-MS). Additionally, we propose, for the first time, that the type II fatty acid synthesis (FAS) pathway in Bacillus is responsible for PA biosynthesis. This finding was confirmed by studying the antimicrobial activity of PA and conducting reverse transcription-quantitative polymerase chain reaction (RT-qPCR) experiments. We also identified numerous genes associated with plant-bacteria interactions in the H2 genome, including more than 94 colonization-related genes, more than 34 antimicrobial genes, and more than 13 plant growth-promoting genes. These findings contribute to our understanding of the biocontrol mechanisms of B. subtilis H2 and have potential applications in crop disease control.

The American Cherimoya Genome Reveals Insights into the Intra-Specific Divergence, the Evolution of Magnoliales, and a Putative Gene Cluster for Acetogenin Biosynthesis

Annona cherimola (cherimoya) is a species renowned for its delectable fruit and medicinal properties. In this study, we developed a chromosome-level genome assembly for the cherimoya 'Booth' cultivar from the United States. The genome assembly has a size of 794 Mb with a N50 = 97.59 Mb. The seven longest scaffolds account for 87.6% of the total genome length, which corresponds to the seven pseudo-chromosomes. A total of 45,272 protein-coding genes (≥30 aa) were predicted with 92.9% gene content completeness. No recent whole genome duplications were identified by an intra-genome collinearity analysis. Phylogenetic analysis supports that eudicots and magnoliids are more closely related to each other than to monocots. Moreover, the Magnoliales was found to be more closely related to the Laurales than the Piperales. Genome comparison revealed that the 'Booth' cultivar has 200 Mb less repeats than the Spanish cultivar 'Fino de Jete', despite their highly similar (>99%) genome sequence identity and collinearity. These two cultivars were diverged during the early Pleistocene (1.93 Mya), which suggests a different origin and domestication of the cherimoya. Terpene/terpenoid metabolism functions were found to be enriched in Magnoliales, while TNL (Toll/Interleukin-1-NBS-LRR) disease resistance gene has been lost in Magnoliales during evolution. We have also identified a gene cluster that is potentially responsible for the biosynthesis of acetogenins, a class of natural products found exclusively in Annonaceae. The cherimoya genome provides an invaluable resource for supporting characterization, conservation, and utilization of Annona genetic resources.

A genome and gene catalog of the aquatic microbiomes of the Tibetan Plateau

The Tibetan Plateau supplies water to nearly 2 billion people in Asia, but climate change poses threats to its aquatic microbial resources. Here, we construct the Tibetan Plateau Microbial Catalog by sequencing 498 metagenomes from six water ecosystems (saline lakes, freshwater lakes, rivers, hot springs, wetlands and glaciers). Our catalog expands knowledge of regional genomic diversity by presenting 32,355 metagenome-assembled genomes that de-replicated into 10,723 representative genome-based species, of which 88% were unannotated. The catalog contains nearly 300 million non-redundant gene clusters, of which 15% novel, and 73,864 biosynthetic gene clusters, of which 50% novel, thus expanding known functional diversity. Using these data, we investigate the Tibetan Plateau aquatic microbiome's biogeography along a distance of 2,500 km and >5 km in altitude. Microbial compositional similarity and the shared gene count with the Tibetan Plateau microbiome decline along with distance and altitude difference, suggesting a dispersal pattern. The Tibetan Plateau Microbial Catalog stands as a substantial repository for high-altitude aquatic microbiome resources, providing potential for discovering novel lineages and functions, and bridging knowledge gaps in microbiome biogeography.

Diversity and modularity of tyrosine-accepting tRNA-like structures

Certain positive-sense single-stranded RNA viruses contain elements at their 3' termini that structurally mimic tRNAs. These tRNA-like structures (TLSs) are classified based on which amino acid is covalently added to the 3' end by host aminoacyl-tRNA synthetase. Recently, a cryoEM reconstruction of a representative tyrosine-accepting tRNA-like structure (TLSTyr) from brome mosaic virus (BMV) revealed a unique mode of recognition of the viral anticodon-mimicking domain by tyrosyl-tRNA synthetase. Some viruses in the hordeivirus genus of Virgaviridae are also selectively aminoacylated with tyrosine, yet these TLS RNAs have a different architecture in the 5' domain that comprises the atypical anticodon loop mimic. Herein, we present bioinformatic and biochemical data supporting a distinct secondary structure for the 5' domain of the hordeivirus TLSTyr compared to those in Bromoviridae Despite forming a different secondary structure, the 5' domain is necessary to achieve robust in vitro aminoacylation. Furthermore, a chimeric RNA containing the 5' domain from the BMV TLSTyr and the 3' domain from a hordeivirus TLSTyr are aminoacylated, illustrating modularity in these structured RNA elements. We propose that the structurally distinct 5' domain of the hordeivirus TLSTyrs performs the same role in mimicking the anticodon loop as its counterpart in the BMV TLSTyr Finally, these structurally and phylogenetically divergent types of TLSTyr provide insight into the evolutionary connections between all classes of viral tRNA-like structures.

Insights into the nuclear and mitochondrial genome of the Lemon shark Negaprion brevirostris using low-coverage sequencing: Genome size, repetitive elements, mitochondrial genome, and phylogenetic placement

The Lemon shark Negaprion brevirostris is an important species experiencing conservation issues that is in need of genomic resources. Herein, we conducted a genome survey sequencing in N. brevirostris and determined genome size, explored repetitive elements, assembled and annotated the 45S rRNA DNA operon, and assembled and described in detail the mitochondrial genome. Lastly, the phylogenetic position of N. brevirostris in the family Carcharhinidae was examined using translated protein coding genes. The estimated haploid genome size ranged between 2.29 and 2.58 Gbp using a k-mer analysis, which is slightly below the genome size estimated for other sharks belonging to the family Carcharhinidae. Using a k-mer analysis, approx. 64-71 % of the genome of N. brevirostris was composed of repetitive elements. A relatively large proportion of the 'repeatome' could not be annotated. Taking into account only annotated repetitive elements, Class I - Long Interspersed Nuclear Element (LINE) were the most abundant repetitive elements followed by Class I - Penelope and Satellite DNA. The nuclear ribosomal operon was fully assembled. The AT-rich complete mitochondrial genome was 16,703 bp long and encoded 13 protein coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Negaprion brevirostris is closely related to the genera Carcharhinus, Glyphis and Lamiopsis in the family Carcharinidae. This new genomic resources will aid with the development of conservation plans for this large coastal shark.

Two chromosome-level genomes of Smittia aterrima and Smittia pratorum (Diptera, Chironomidae)

Chironomids are one of the most abundant aquatic insects and are widely distributed in various biological communities. However, the lack of high-quality genomes has hindered our ability to study the evolution and ecology of this group. Here, we used Nanopore long reads and Hi-C data to produce two chromosome-level genomes from mixed genomic data. The genomes of Smittia aterrima (SateA) and Smittia pratorum (SateB) were assembled into three chromosomes, with sizes of 78.45 Mb and 71.56 Mb, scaffold N50 lengths of 25.73 and 23.53 Mb, and BUSCO completeness of 98.5% and 97.8% (n = 1,367), 5.68 Mb (7.24%) and 1.94 Mb (2.72%) of repetitive elements, and predicted 12,330 (97.70% BUSCO completeness) and 11,250 (97.40%) protein-coding genes, respectively. These high-quality genomes will serve as valuable resources for comprehending the evolution and environmental adaptation of chironomids.

Chromosome-level genome assembly of humpback grouper using PacBio HiFi reads and Hi-C technologies

The humpback grouper (Cromileptes altivelis), a medium-sized coral reef teleost, is a naturally rare species distributed in the tropical waters of the Indian and Pacific Oceans. It has high market value, but artificial reproduction and breeding remain limited and need to be improved. Here, we assembled the genome with 1.08 Gb, with a contig N50 of 43.78 Mb. A total of 96.59% of the assembly anchored to 24 pseudochromosomes using Hi-C technology. It contained 24,442 protein-coding sequences, of which 99.3% were functionally annotated. The completeness of the assembly was estimated to be 97.3% using BUSCO. The phylogenomic analysis suggested that humpback grouper should be classified into the genus Epinephelus rather than Cromileptes. The comparative genomic analysis revealed that the gene families related to circadian entrainment were significantly expanded. The high-quality reference genome provides useful genomic tools for exploiting the genomic resource of humpback grouper and supports the functional genomic study of this species in the future.

A chromosome-level genome assembly for Onobrychis viciifolia reveals gene copy number gain underlying enhanced proanthocyanidin biosynthesis

Sainfoin (Onobrychis viciifolia), which belongs to subfamily Papilionoideae of Leguminosae, is a vital perennial forage known as "holy hay" due to its high contents of crude proteins and proanthocyanidins (PAs, also called condensed tannins) that have various pharmacological properties in animal feed, such as alleviating rumen tympanic disease in ruminants. In this study, we select an autotetraploid common sainfoin (2n = 4x = 28) and report its high-quality chromosome-level genome assembly with 28 pseudochromosomes and four haplotypes (~1950.14 Mb, contig N50 = 10.91 Mb). The copy numbers of genes involved in PA biosynthesis in sainfoin are significantly greater than those in four selected Fabales species, namely, autotetraploid Medicago sativa and three other diploid species, Lotus japonicus, Medicago truncatula, and Glycine max. Furthermore, gene expansion is confirmed to be the key contributor to the increased expression of these genes and subsequent PA enhancement in sainfoin. Transcriptomic analyses reveal that the expression of genes involved in the PA biosynthesis pathway is significantly increased in the lines with high PA content compared to the lines with medium and low PA content. The sainfoin genome assembly will improve our understanding of leguminous genome evolution and biosynthesis of secondary metabolites in sainfoin.

The Nucleic Acid Knowledgebase: a new portal for 3D structural information about nucleic acids

The Nucleic Acid Knowledgebase (nakb.org) is a new data resource, updated weekly, for experimentally determined 3D structures containing DNA and/or RNA nucleic acid polymers and their biological assemblies. NAKB indexes nucleic acid-containing structures derived from all major structure determination methods (X-ray, NMR and EM), including all held by the Protein Data Bank (PDB). As the planned successor to the Nucleic Acid Database (NDB), NAKB's design preserves all functionality of the NDB and provides novel nucleic acid-centric content, including structural and functional annotations, as well as annotations from and links to external resources. A variety of custom interactive tools have been developed to enable rapid exploration and drill-down of NAKB's content.

Exploring the cobia (Rachycentron canadum) genome: unveiling putative male heterogametic regions and identification of sex-specific markers

BACKGROUND

Cobia (Rachycentron canadum) is the only member of the Rachycentridae family and exhibits considerable sexual dimorphism in growth rate. Sex determination in teleosts has been a long-standing basic biological question, and the molecular mechanisms of sex determination/differentiation in cobia are completely unknown.

RESULTS

Here, we reported 2 high-quality, chromosome-level annotated male and female cobia genomes with assembly sizes of 586.51 Mb (contig/scaffold N50: 86.0 kb/24.3 Mb) and 583.88 Mb (79.9 kb/22.5 Mb), respectively. Synteny inference among perciform genomes revealed that cobia and the remora Echeneis naucrates were sister groups. Further, whole-genome resequencing of 31 males and 60 females, genome-wide association study, and sequencing depth analysis identified 3 short male-specific regions within a 10.7-kb continuous genomic region on male chromosome 18, which hinted at an undifferentiated sex chromosome system with a putative XX/XY mode of sex determination in cobia. Importantly, the only 2 genes within/between the male-specific regions, epoxide hydrolase 1 (ephx1, renamed cephx1y) and transcription factor 24 (tcf24, renamed ctcf24y), showed testis-specific/biased gene expression, whereas their counterparts cephx1x and ctf24x, located in female chromosome 18, were similarly expressed in both sexes. In addition, male-specific PCR targeting the cephx1y gene revealed that this genomic feature is conserved in cobia populations from Panama, Brazil, Australia, and Japan.

CONCLUSION

The first comprehensive genomic survey presented here is a valuable resource for future studies on cobia population structure and dynamics, conservation, and evolutionary history. Furthermore, it establishes evidence of putative male heterogametic regions with 2 genes playing a potential role in the sex determination of the species, and it provides further support for the rapid evolution of sex-determining mechanisms in teleost fish.

Classification and Identification of Non-canonical Base Pairs and Structural Motifs

The 3D structures of many ribonucleic acid (RNA) loops are characterized by highly organized networks of non-canonical interactions. Multiple computational methods have been developed to annotate structures with those interactions or automatically identify recurrent interaction networks. By contrast, the reverse problem that aims to retrieve the geometry of a look from its sequence or ensemble of interactions remains much less explored. In this chapter, we will describe how to retrieve and build families of conserved structural motifs using their underlying network of non-canonical interactions. Then, we will show how to assign sequence alignments to those families and use the software BayesPairing to build statistical models of structural motifs with their associated sequence alignments. From this model, we will apply BayesPairing to identify in new sequences regions where those loop geometries can occur.

Evolutionary Structure Conservation and Covariance Scores

Effective homology search for non-coding RNAs is frequently not possible via sequence similarity alone. Current methods leverage evolutionary information like structure conservation or covariance scores to identify homologs in organisms that are phylogenetically more distant. In this chapter, we introduce the theoretical background of evolutionary structure conservation and covariance score, and we show hands-on how current methods in the field are applied on example datasets.

A chromosome-level genome assembly for the Silkie chicken resolves complete sequences for key chicken metabolic, reproductive, and immunity genes

A set of high-quality pan-genomes would help identify important genes that are still hidden/incomplete in bird reference genomes. In an attempt to address these issues, we have assembled a de novo chromosome-level reference genome of the Silkie (Gallus gallus domesticus), which is an important avian model for unique traits, like fibromelanosis, with unclear genetic foundation. This Silkie genome includes the complete genomic sequences of well-known, but unresolved, evolutionarily, endocrinologically, and immunologically important genes, including leptin, ovocleidin-17, and tumor-necrosis factor-α. The gap-less and manually annotated MHC (major histocompatibility complex) region possesses 38 recently identified genes, with differentially regulated genes recovered in response to pathogen challenges. We also provide whole-genome methylation and genetic variation maps, and resolve a complex genetic region that may contribute to fibromelanosis in these animals. Finally, we experimentally show leptin binding to the identified leptin receptor in chicken, confirming an active leptin ligand-receptor system. The Silkie genome assembly not only provides a rich data resource for avian genome studies, but also lays a foundation for further functional validation of resolved genes.

Annotating Macromolecular Complexes in the Protein Data Bank: Improving the FAIRness of Structure Data

Macromolecular complexes are essential functional units in nearly all cellular processes, and their atomic-level understanding is critical for elucidating and modulating molecular mechanisms. The Protein Data Bank (PDB) serves as the global repository for experimentally determined structures of macromolecules. Structural data in the PDB offer valuable insights into the dynamics, conformation, and functional states of biological assemblies. However, the current annotation practices lack standardised naming conventions for assemblies in the PDB, complicating the identification of instances representing the same assembly. In this study, we introduce a method leveraging resources external to PDB, such as the Complex Portal, UniProt and Gene Ontology, to describe assemblies and contextualise them within their biological settings accurately. Employing the proposed approach, we assigned standard names to over 90% of unique assemblies in the PDB and provided persistent identifiers for each assembly. This standardisation of assembly data enhances the PDB, facilitating a deeper understanding of macromolecular complexes. Furthermore, the data standardisation improves the PDB's FAIR attributes, fostering more effective basic and translational research and scientific education.

High-quality assembly and methylome of a Tibetan wild tree peony genome (Paeonia ludlowii) reveal the evolution of giant genome architecture

Tree peony belongs to one of the Saxifragales families, Paeoniaceae. It is one of the most famous ornamental plants, and is also a promising woody oil plant. Although two Paeoniaceae genomes have been released, their assembly qualities are still to be improved. Additionally, more genomes from wild peonies are needed to accelerate genomic-assisted breeding. Here we assemble a high-quality and chromosome-scale 10.3-Gb genome of a wild Tibetan tree peony, Paeonia ludlowii, which features substantial sequence divergence, including around 75% specific sequences and gene-level differentials compared with other peony genomes. Our phylogenetic analyses suggest that Saxifragales and Vitales are sister taxa and, together with rosids, they are the sister taxon to asterids. The P. ludlowii genome is characterized by frequent chromosome reductions, centromere rearrangements, broadly distributed heterochromatin, and recent continuous bursts of transposable element (TE) movement in peony, although it lacks recent whole-genome duplication. These recent TE bursts appeared during the uplift and glacial period of the Qinghai-Tibet Plateau, perhaps contributing to adaptation to rapid climate changes. Further integrated analyses with methylome data revealed that genome expansion in peony might be dynamically affected by complex interactions among TE proliferation, TE removal, and DNA methylation silencing. Such interactions also impact numerous recently duplicated genes, particularly those related to oil biosynthesis and flower traits. This genome resource will not only provide the genomic basis for tree peony breeding but also shed light on the study of the evolution of huge genome structures as well as their protein-coding genes.

Whole genome sequencing provides evidence for Bacillus velezensis SH-1471 as a beneficial rhizosphere bacterium in plants

Bacillus is widely used in agriculture due to its diverse biological activities. We isolated a Bacillus velezensis SH-1471 from the rhizosphere soil of healthy tobacco, which has broad-spectrum antagonistic activity against a variety of plant pathogenic fungi such as Fusarium oxysporum, and can be colonized in the rhizosphere of a variety of plants. This study will further explore its mechanism by combining biological and molecular biology methods. SH-1471 contains a ring chromosome of 4,181,346 bp with a mean G + C content of 46.18%. We identified 14 homologous genes related to biosynthesis of resistant secondary metabolite, and three clusters encoded potential new antibacterial substances. It also contains a large number of genes from colonizing bacteria and genes related to plant bacterial interactions. It also contains genes related to environmental stress, as well as genes related to drug resistance. We also found that there are many metabolites in the strain that can inhibit the growth of pathogens. In addition, our indoor pot test found that SH-1471 has a good control effect on tomato wilt, and could significantly improve plant height, stem circumference, root length, root weight, and fresh weight and dry weight of the aboveground part of tomato seedlings. Therefore, SH-1471 is a potential biological control strain with important application value. The results of this study will help to further study the mechanism of SH-1471 in biological control of plant diseases and promote its application.

Long noncoding RNAs emerge from transposon-derived antisense sequences and may contribute to infection stage-specific transposon regulation in a fungal phytopathogen

BACKGROUND

The genome of the obligate biotrophic phytopathogenic barley powdery mildew fungus Blumeria hordei is inflated due to highly abundant and possibly active transposable elements (TEs). In the absence of the otherwise common repeat-induced point mutation transposon defense mechanism, noncoding RNAs could be key for regulating the activity of TEs and coding genes during the pathogenic life cycle.

RESULTS

We performed time-course whole-transcriptome shotgun sequencing (RNA-seq) of total RNA derived from infected barley leaf epidermis at various stages of fungal pathogenesis and observed significant transcript accumulation and time point-dependent regulation of TEs in B. hordei. Using a manually curated consensus database of 344 TEs, we discovered phased small RNAs mapping to 104 consensus transposons, suggesting that RNA interference contributes significantly to their regulation. Further, we identified 5,127 long noncoding RNAs (lncRNAs) genome-wide in B. hordei, of which 823 originated from the antisense strand of a TE. Co-expression network analysis of lncRNAs, TEs, and coding genes throughout the asexual life cycle of B. hordei points at extensive positive and negative co-regulation of lncRNAs, subsets of TEs and coding genes.

CONCLUSIONS

Our work suggests that similar to mammals and plants, fungal lncRNAs support the dynamic modulation of transcript levels, including TEs, during pivotal stages of host infection. The lncRNAs may support transcriptional diversity and plasticity amid loss of coding genes in powdery mildew fungi and may give rise to novel regulatory elements and virulence peptides, thus representing key drivers of rapid evolutionary adaptation to promote pathogenicity and overcome host defense.

Targeting FMN, TPP, SAM-I, and glmS Riboswitches with Chimeric Antisense Oligonucleotides for Completely Rational Antibacterial Drug Development

Antimicrobial drug resistance has emerged as a significant challenge in contemporary medicine due to the proliferation of numerous bacterial strains resistant to all existing antibiotics. Meanwhile, riboswitches have emerged as promising targets for discovering antibacterial drugs. Riboswitches are regulatory elements in certain bacterial mRNAs that can bind to specific molecules and control gene expression via transcriptional termination, prevention of translation, or mRNA destabilization. By targeting riboswitches, we aim to develop innovative strategies to combat antibiotic-resistant bacteria and enhance the efficacy of antibacterial treatments. This convergence of challenges and opportunities underscores the ongoing quest to revolutionize medical approaches against evolving bacterial threats. For the first time, this innovative review describes the rational design and applications of chimeric antisense oligonucleotides as antibacterial agents targeting four riboswitches selected based on genome-wide bioinformatic analyses. The antisense oligonucleotides are coupled with the cell-penetrating oligopeptide pVEC, which penetrates Gram-positive and Gram-negative bacteria and specifically targets glmS, FMN, TPP, and SAM-I riboswitches in Staphylococcus aureus, Listeria monocytogenes, and Escherichia coli. The average antibiotic dosage of antisense oligonucleotides that inhibits 80% of bacterial growth is around 700 nM (4.5 μg/mL). Antisense oligonucleotides do not exhibit toxicity in human cell lines at this concentration. The results demonstrate that these riboswitches are suitable targets for antibacterial drug development using antisense oligonucleotide technology. The approach is fully rational because selecting suitable riboswitch targets and designing ASOs that target them are based on predefined criteria. The approach can be used to develop narrow or broad-spectrum antibiotics against multidrug-resistant bacterial strains for a short time. The approach is easily adaptive to new resistance using targeting NGS technology.

Whole-Genome Analysis of Termite-Derived Bacillus velezensis BV-10 and Its Application in King Grass Silage

Bacillus velezensis (B. velezensis) is a cellulose-degrading strain that has the potential as an additive in fermented feed. B. velezensis BV-10 was isolated and screened from the termite gut. We sequenced the whole genome of this new source of B. velezensis to reveal its potential for use in cellulose degradation. Whole-genome sequencing of B. velezensis BV-10 showed that it has a circular chromosome of 3929792 bp containing 3873 coding genes with a GC content of 45.51% and many genes related to cellulose, hemicellulose, and lignin degradation. King grass silage was inoculated with B. velezensis BV-10 and mixed with other feed additives to assess the effect of B. velezensis BV-10 on the fermentation quality of silage. Six treatment groups were established: the control, B. velezensis BV-10, molasses, cellulase, B. velezensis BV-10 plus molasses, and B. velezensis BV-10 plus cellulase groups. After 30 days of silage-fermentation testing, B. velezensis BV-10 was found to rapidly reduce the silage pH value and significantly reduce the acid-detergent fiber (ADF) content (p < 0.05). The addition of B. velezensis BV-10 plus molasses and cellulase in fermented feed significantly reduced the silage neutral-detergent fiber and ADF content and promoted organic-acid accumulation (p < 0.05). The above results demonstrate that B. velezensis BV-10 promotes the fermentation quality of silage and that this effect is greater when other silage-fermentation additives are included. In conclusion, genes involved in cellulose degradation in B. velezensis BV-10 were identified by whole-genome sequencing and further experiments explored the effects of B. velezensis BV-10 and different feed additives on the fermentation quality of king grass silage, revealing the potential of Bacillus velezensis as a new silage additive.

MicroRNA397 regulates tolerance to drought and fungal infection by regulating lignin deposition in chickpea root

Plants deposit lignin in the secondary cell wall as a common response to drought and pathogen attacks. Cell wall localised multicopper oxidase family enzymes LACCASES (LACs) catalyse the formation of monolignol radicals and facilitate lignin formation. We show an upregulation of the expression of several LAC genes and a downregulation of microRNA397 (CamiR397) in response to natural drought in chickpea roots. CamiR397 was found to target LAC4 and LAC17L out of twenty annotated LACs in chickpea. CamiR397 and its target genes are expressed in the root. Overexpression of CamiR397 reduced expression of LAC4 and LAC17L and lignin deposition in chickpea root xylem causing reduction in xylem wall thickness. Downregulation of CamiR397 activity by expressing a short tandem target mimic (STTM397) construct increased root lignin deposition in chickpea. CamiR397-overexpressing and STTM397 chickpea lines showed sensitivity and tolerance, respectively, towards natural drought. Infection with a fungal pathogen Macrophomina phaseolina, responsible for dry root rot (DRR) disease in chickpea, induced local lignin deposition and LAC gene expression. CamiR397-overexpressing and STTM397 chickpea lines showed more sensitivity and tolerance, respectively, to DRR. Our results demonstrated the regulatory role of CamiR397 in root lignification during drought and DRR in an agriculturally important crop chickpea.

A transient intermediate RNA structure underlies the regulatory function of the E. coli thiB TPP translational riboswitch

Riboswitches are cis-regulatory RNA elements that regulate gene expression in response to ligand binding through the coordinated action of a ligand-binding aptamer domain (AD) and a downstream expression platform (EP). Previous studies of transcriptional riboswitches have uncovered diverse examples that utilize structural intermediates that compete with the AD and EP folds to mediate the switching mechanism on the timescale of transcription. Here we investigate whether similar intermediates are important for riboswitches that control translation by studying the Escherichia coli thiB thiamin pyrophosphate (TPP) riboswitch. Using cellular gene expression assays, we first confirmed that the riboswitch acts at the level of translational regulation. Deletion mutagenesis showed the importance of the AD-EP linker sequence for riboswitch function. Sequence complementarity between the linker region and the AD P1 stem suggested the possibility of an intermediate nascent RNA structure called the antisequestering stem that could mediate the thiB switching mechanism. Experimentally informed secondary structure models of the thiB folding pathway generated from chemical probing of nascent thiB structures in stalled transcription elongation complexes confirmed the presence of the antisequestering stem, and showed it may form cotranscriptionally. Additional mutational analysis showed that mutations to the antisequestering stem break or bias thiB function according to whether the antisequestering stem or P1 is favored. This work provides an important example of intermediate structures that compete with AD and EP folds to implement riboswitch mechanisms.

Genome Analysis and Optimization of Caproic Acid Production of Clostridium butyricum GD1-1 Isolated from the Pit Mud of Nongxiangxing Baijiu

Caproic acid is a precursor substance for the synthesis of ethyl caproate, the main flavor substance of nongxiangxing baijiu liquor. In this study, Clostridium butyricum GD1-1, a strain with high caproic acid concentration (3.86 g/l), was isolated from the storage pit mud of nongxiangxing baijiu for sequencing and analysis. The strain's genome was 3,840,048 bp in length with 4,050 open reading frames. In addition, virulence factor annotation analysis showed C. butyricum GD1-1 to be safe at the genetic level. However, the annotation results using the Kyoto Encyclopedia of Genes and Genomes Automatic Annotation Server predicted a deficiency in the strain's synthesis of alanine, methionine, and biotin. These results were confirmed by essential nutrient factor validation experiments. Furthermore, the optimized medium conditions for caproic acid concentration by strain GD1-1 were (g/l): glucose 30, NaCl 5, yeast extract 10, peptone 10, beef paste 10, sodium acetate 11, L-cysteine 0.6, biotin 0.004, starch 2, and 2.0% ethanol. The optimized fermentation conditions for caproic acid production by C. butyricum GD1-1 on a single-factor basis were: 5% inoculum volume, 35°C, pH 7, and 90% loading volume. Under optimal conditions, the caproic acid concentration of strain GD1-1 reached 5.42 g/l, which was 1.40 times higher than the initial concentration. C. butyricum GD1-1 could be further used in caproic acid production, NXXB pit mud strengthening and maintenance, and artificial pit mud preparation.

Diverse regulatory pathways modulate bet hedging of competence induction in epigenetically-differentiated phase variants of Streptococcus pneumoniae

Despite enabling Streptococcus pneumoniae to acquire antibiotic resistance and evade vaccine-induced immunity, transformation occurs at variable rates across pneumococci. Phase variants of isolate RMV7, distinguished by altered methylation patterns driven by the translocating variable restriction-modification (tvr) locus, differed significantly in their transformation efficiencies and biofilm thicknesses. These differences were replicated when the corresponding tvr alleles were introduced into an RMV7 derivative lacking the locus. RNA-seq identified differential expression of the type 1 pilus, causing the variation in biofilm formation, and inhibition of competence induction in the less transformable variant, RMV7domi. This was partly attributable to RMV7domi's lower expression of ManLMN, which promoted competence induction through importing N-acetylglucosamine. This effect was potentiated by analogues of some proteobacterial competence regulatory machinery. Additionally, one of RMV7domi's phage-related chromosomal island was relatively active, which inhibited transformation by increasing expression of the stress response proteins ClpP and HrcA. However, HrcA increased competence induction in the other variant, with its effects depending on Ca2+ supplementation and heat shock. Hence the heterogeneity in transformation efficiency likely reflects the diverse signalling pathways by which it is affected. This regulatory complexity will modulate population-wide responses to synchronising quorum sensing signals to produce co-ordinated yet stochastic bet hedging behaviour.

Hibernation-Induced microRNA Expression Promotes Signaling Pathways and Cell Cycle Dysregulation in Ictidomys tridecemlineatus Cardiac Tissue

The thirteen-lined ground squirrel Ictidomys tridecemlineatus is a rodent that lives throughout the United States and Canada and uses metabolic rate depression to facilitate circannual hibernation which helps it survive the winter. Metabolic rate depression is the reorganization of cellular physiology and molecular biology to facilitate a global downregulation of nonessential genes and processes, which conserves endogenous fuel resources and prevents the buildup of waste byproducts. Facilitating metabolic rate depression requires a complex interplay of regulatory approaches, including post-transcriptional modes such as microRNA. MicroRNA are short, single-stranded RNA species that bind to mRNA transcripts and target them for degradation or translational suppression. Using next-generation sequencing, we analyzed euthermic vs. hibernating cardiac tissue in I. tridecemlineatus to predict seven miRNAs (let-7e-5p, miR-122-5p, miR-2355-3p, miR-6715b-3p, miR-378i, miR-9851-3p, and miR-454-3p) that may be differentially regulated during hibernation. Gene ontology and KEGG pathway analysis suggested that these miRNAs cause a strong activation of ErbB2 signaling which causes downstream effects, including the activation of MAPK and PI3K/Akt signaling and concurrent decreases in p53 signaling and cell cycle-related processes. Taken together, these results predict critical miRNAs that may change during hibernation in the hearts of I. tridecemlineatus and identify key signaling pathways that warrant further study in this species.

The SMN complex drives structural changes in human snRNAs to enable snRNP assembly

Spliceosomal snRNPs are multicomponent particles that undergo a complex maturation pathway. Human Sm-class snRNAs are generated as 3'-end extended precursors, which are exported to the cytoplasm and assembled together with Sm proteins into core RNPs by the SMN complex. Here, we provide evidence that these pre-snRNA substrates contain compact, evolutionarily conserved secondary structures that overlap with the Sm binding site. These structural motifs in pre-snRNAs are predicted to interfere with Sm core assembly. We model structural rearrangements that lead to an open pre-snRNA conformation compatible with Sm protein interaction. The predicted rearrangement pathway is conserved in Metazoa and requires an external factor that initiates snRNA remodeling. We show that the essential helicase Gemin3, which is a component of the SMN complex, is crucial for snRNA structural rearrangements during snRNP maturation. The SMN complex thus facilitates ATP-driven structural changes in snRNAs that expose the Sm site and enable Sm protein binding.

lhCLIP reveals the in vivo RNA-RNA interactions recognized by hnRNPK

RNA-RNA interactions play a crucial role in regulating gene expression and various biological processes, but identifying these interactions on a transcriptomic scale remains a challenge. To address this, we have developed a new biochemical technique called pCp-biotin labelled RNA hybrid and ultraviolet crosslinking and immunoprecipitation (lhCLIP) that enables the transcriptome-wide identification of intra- and intermolecular RNA-RNA interactions mediated by a specific RNA-binding protein (RBP). Using lhCLIP, we have uncovered a diverse landscape of intermolecular RNA interactions recognized by hnRNPK in human cells, involving all major classes of noncoding RNAs (ncRNAs) and mRNA. Notably, hnRNPK selectively binds with snRNA U4, U11, and U12, and shapes the secondary structure of these snRNAs, which may impact RNA splicing. Our study demonstrates the potential of lhCLIP as a user-friendly and widely applicable method for discovering RNA-RNA interactions mediated by a particular protein of interest and provides a valuable tool for further investigating the role of RBPs in gene expression and biological processes.

Uncultured prokaryotic genomes in the spotlight: An examination of publicly available data from metagenomics and single-cell genomics

Owing to the ineffectiveness of traditional culture techniques for the vast majority of microbial species, culture-independent analyses utilizing next-generation sequencing and bioinformatics have become essential for gaining insight into microbial ecology and function. This mini-review focuses on two essential methods for obtaining genetic information from uncultured prokaryotes, metagenomics and single-cell genomics. We analyzed the registration status of uncultured prokaryotic genome data from major public databases and assessed the advantages and limitations of both the methods. Metagenomics generates a significant quantity of sequence data and multiple prokaryotic genomes using straightforward experimental procedures. However, in ecosystems with high microbial diversity, such as soil, most genes are presented as brief, disconnected contigs, and lack association of highly conserved genes and mobile genetic elements with individual species genomes. Although technically more challenging, single-cell genomics offers valuable insights into complex ecosystems by providing strain-resolved genomes, addressing issues in metagenomics. Recent technological advancements, such as long-read sequencing, machine learning algorithms, and in silico protein structure prediction, in combination with vast genomic data, have the potential to overcome the current technical challenges and facilitate a deeper understanding of uncultured microbial ecosystems and microbial dark matter genes and proteins. In light of this, it is imperative that continued innovation in both methods and technologies take place to create high-quality reference genome databases that will support future microbial research and industrial applications.

Characterization of the male-specific region containing the candidate sex-determining gene in Amur catfish (Silurus asotus) using third-generation- and pool-sequencing data

Amur catfish (Silurus asotus) is an ecologically and economically important fish species in Asia. Here, we assembled the female and male Amur catfish genomes, with genome sizes of 757.15 and 755.44 Mb, respectively, at the chromosome level using nanopore and Hi-C technologies. Consistent with the known diploid chromosome count, both genomes contained 29 chromosome-size scaffolds covering 98.80 and 98.73 % of the complete haplotypic assembly with scaffold N50 of 28.87 and 27.29 Mb, respectively. The female (n = 40) and male (n = 40) pools were re-sequenced. Comparative analysis of sequencing and re-sequencing data from both sexes confirmed the presence of an XX/XY sex determination system in Amur catfish and revealed Chr5 as the sex chromosome containing an approximately 400 kb Y-specific region (MSY). Gene annotation revealed a male-specific duplicate of amhr2, namely amhr2y, in MSY, which is male-specific in different wild populations and expressed only in the testes. Amur catfish shared partially syntenic MSY and amhr2y genes with the southern catfish (S. meridionalis, Chr24), which were located on different chromosomes. High sequence divergence between amhr2y and amhr2 and high sequence similarity with amhr2y were observed in both species. These results indicate the common origin of the sex-determining (SD) gene and transition of amhr2y in the two Silurus species. Accumulation of repetitive elements in the MSY of both species may be the main driver of the transition of amhr2y. Overall, our study provides valuable catfish genomic resources. Moreover, determination of amhr2y as the candidate SD gene in Amur catfish provides another example of amhr2 as the SD gene in fish.

The genome of the colonial hydroid Hydractinia reveals their stem cells utilize a toolkit of evolutionarily shared genes with all animals

Hydractinia is a colonial marine hydroid that exhibits remarkable biological properties, including the capacity to regenerate its entire body throughout its lifetime, a process made possible by its adult migratory stem cells, known as i-cells. Here, we provide an in-depth characterization of the genomic structure and gene content of two Hydractinia species, H. symbiolongicarpus and H. echinata, placing them in a comparative evolutionary framework with other cnidarian genomes. We also generated and annotated a single-cell transcriptomic atlas for adult male H. symbiolongicarpus and identified cell type markers for all major cell types, including key i-cell markers. Orthology analyses based on the markers revealed that Hydractinia's i-cells are highly enriched in genes that are widely shared amongst animals, a striking finding given that Hydractinia has a higher proportion of phylum-specific genes than any of the other 41 animals in our orthology analysis. These results indicate that Hydractinia's stem cells and early progenitor cells may use a toolkit shared with all animals, making it a promising model organism for future exploration of stem cell biology and regenerative medicine. The genomic and transcriptomic resources for Hydractinia presented here will enable further studies of their regenerative capacity, colonial morphology, and ability to distinguish self from non-self.

Inter-comparison of marine microbiome sampling protocols

Research on marine microbial communities is growing, but studies are hard to compare because of variation in seawater sampling protocols. To help researchers in the inter-comparison of studies that use different seawater sampling methodologies, as well as to help them design future sampling campaigns, we developed the EuroMarine Open Science Exploration initiative (EMOSE). Within the EMOSE framework, we sampled thousands of liters of seawater from a single station in the NW Mediterranean Sea (Service d'Observation du Laboratoire Arago [SOLA], Banyuls-sur-Mer), during one single day. The resulting dataset includes multiple seawater processing approaches, encompassing different material-type kinds of filters (cartridge membrane and flat membrane), three different size fractionations (>0.22 µm, 0.22-3 µm, 3-20 µm and >20 µm), and a number of different seawater volumes ranging from 1 L up to 1000 L. We show that the volume of seawater that is filtered does not have a significant effect on prokaryotic and protist diversity, independently of the sequencing strategy. However, there was a clear difference in alpha and beta diversity between size fractions and between these and "whole water" (with no pre-fractionation). Overall, we recommend care when merging data from datasets that use filters of different pore size, but we consider that the type of filter and volume should not act as confounding variables for the tested sequencing strategies. To the best of our knowledge, this is the first time a publicly available dataset effectively allows for the clarification of the impact of marine microbiome methodological options across a wide range of protocols, including large-scale variations in sampled volume.

Terabase-Scale Coassembly of a Tropical Soil Microbiome

Petabases of environmental metagenomic data are publicly available, presenting an opportunity to characterize complex environments and discover novel lineages of life. Metagenome coassembly, in which many metagenomic samples from an environment are simultaneously analyzed to infer the underlying genomes' sequences, is an essential tool for achieving this goal. We applied MetaHipMer2, a distributed metagenome assembler that runs on supercomputing clusters, to coassemble 3.4 terabases (Tbp) of metagenome data from a tropical soil in the Luquillo Experimental Forest (LEF), Puerto Rico. The resulting coassembly yielded 39 high-quality (>90% complete, <5% contaminated, with predicted 23S, 16S, and 5S rRNA genes and ≥18 tRNAs) metagenome-assembled genomes (MAGs), including two from the candidate phylum Eremiobacterota. Another 268 medium-quality (≥50% complete, <10% contaminated) MAGs were extracted, including the candidate phyla Dependentiae, Dormibacterota, and Methylomirabilota. In total, 307 medium- or higher-quality MAGs were assigned to 23 phyla, compared to 294 MAGs assigned to nine phyla in the same samples individually assembled. The low-quality (<50% complete, <10% contaminated) MAGs from the coassembly revealed a 49% complete rare biosphere microbe from the candidate phylum FCPU426 among other low-abundance microbes, an 81% complete fungal genome from the phylum Ascomycota, and 30 partial eukaryotic MAGs with ≥10% completeness, possibly representing protist lineages. A total of 22,254 viruses, many of them low abundance, were identified. Estimation of metagenome coverage and diversity indicates that we may have characterized ≥87.5% of the sequence diversity in this humid tropical soil and indicates the value of future terabase-scale sequencing and coassembly of complex environments. IMPORTANCE Petabases of reads are being produced by environmental metagenome sequencing. An essential step in analyzing these data is metagenome assembly, the computational reconstruction of genome sequences from microbial communities. "Coassembly" of metagenomic sequence data, in which multiple samples are assembled together, enables more complete detection of microbial genomes in an environment than "multiassembly," in which samples are assembled individually. To demonstrate the potential for coassembling terabases of metagenome data to drive biological discovery, we applied MetaHipMer2, a distributed metagenome assembler that runs on supercomputing clusters, to coassemble 3.4 Tbp of reads from a humid tropical soil environment. The resulting coassembly, its functional annotation, and analysis are presented here. The coassembly yielded more, and phylogenetically more diverse, microbial, eukaryotic, and viral genomes than the multiassembly of the same data. Our resource may facilitate the discovery of novel microbial biology in tropical soils and demonstrates the value of terabase-scale metagenome sequencing.