KAAS: an automatic genome annotation and pathway reconstruction server

Identification and expression profiling of microRNAs in leaf tissues of Foeniculum vulgare Mill. under salinity stress

Foeniculum vulgare Mill. commonly known as fennel, is a globally recognized aromatic medicinal plant and culinary herb with widespread popularity due to its antimicrobial, antioxidant, carminative, and diuretic properties, among others. Although the phenotypic effects of salinity stress have been previously explored in fennel, the molecular mechanisms underlying responses to elevated salinity in this plant remain elusive. MicroRNAs (miRNAs) are tiny, endogenous, and extensively conserved non-coding RNAs (ncRNAs) typically ranging from 20 to 24 nucleotides (nt) in length that play a major role in a myriad of biological functions. In fact, a number of miRNAs have been extensively associated with responses to abiotic stress in plants. Consequently, employing computational methodologies and rigorous filtering criteria, 40 putative miRNAs belonging to 25 different families were characterized from fennel in this study. Subsequently, employing the psRNATarget tool, a total of 67 different candidate target transcripts for the characterized fennel miRNAs were predicted. Additionally, the expression patterns of six selected fennel miRNAs (i.e. fvu-miR156a, fvu-miR162a-3p, fvu-miR166a-3p, fvu-miR167a-5p, fvu-miR171a-3p, and fvu-miR408-3p) were analyzed under salinity stress conditions via qPCR. This article holds notable significance as it identifies not only 40 putative miRNAs in fennel, a non-model plant, but also pioneers the analysis of their expression under salinity stress conditions.

Groundwater chromate removal by autotrophic sulfur disproportionation

Chromate [Cr(VI)] contamination in groundwater is a global environmental challenge. Traditional elemental sulfur-based biotechnologies for Cr(VI) removal depend heavily on the synthesis of dissolved organic carbon to fuel heterotrophic Cr(VI) reduction, a bottleneck in the remediation process. Here we show an alternative approach by leveraging sulfur-disproportionating bacteria (SDB) inherent to groundwater ecosystems, offering a novel and efficient Cr(VI) removal strategy. We implemented SDB within a sulfur-packed bed reactor for treating Cr(VI)-contaminated groundwater, achieving a notable removal rate of 6.19 mg L-1 h-1 under oligotrophic conditions. We identified the chemical reduction of Cr(VI) via sulfide, produced through sulfur disproportionation, as a key mechanism, alongside microbial Cr(VI) reduction within the sulfur-based biosystem. Genome-centric metagenomic analysis revealed a symbiotic relationship among SDB, sulfur-oxidizing, and chromate-reducing bacteria within the reactor, suggesting that Cr(VI) detoxification by these microbial communities enhances the sulfur-disproportionation process. This research highlights the significance of sulfur disproportionation in the cryptic sulfur cycle in Cr(VI)-contaminated groundwater and proposes its practical application in groundwater remediation efforts.

Coping with salinity extremes: Gill transcriptome profiling in the black-chinned tilapia (Sarotherodon melanotheron)

Steeper and sometimes extreme salinity gradients increasingly affect aquatic organisms because of climate change. Hypersalinity habitats demand powerful physiological adaptive strategies. Few teleost species have the capacity to spend their whole life cycle in salinities way over seawater levels. Focusing on the multifunctional gill, we unraveled the tilapia S. melanotheron key strategies to cope with different environmental conditions, ranging from freshwater up to hypersaline habitats. De novo transcriptome assembly based on RNAseq allowed for the analysis of 40,967 annotated transcripts among samples collected in three wild populations at 0, 40 and 80 ‰. A trend analysis of the expression patterns revealed responses across the salinity gradient with different gene pathways involved. Genes linked to ion transport, pH regulation and cell surface receptor signaling were mainly upregulated in the high salinity habitat. We identified tight junction proteins that were critical in high salinity habitats and that were different from the well-known tightening junctional proteins identified and expressed in fresh water. Expression profiles also suggest a change in the vascular tone that could be linked to an osmorespiratory compromise not only in fresh water, but also in high salinity environments. A striking downregulation of genes linked to the immune system and to the heat shock response was observed suggesting an energetic trade-off between immunity and acclimation/adaptation in the hypersaline habitat. The high expression of transcripts coding for immune and heat shock response in the freshwater habitat suggests the establishment of powerful mechanisms to protect gills from environmental threats and to maintain protein integrity. Non-directional expression trends were also detected with an upregulation of genes only in the hypersaline habitat (80 ‰) or only in the marine habitat (40 ‰). Unravel physiological strategies in S. melanotheron populations will help to better understand the molecular basis of fish euryhalinity in salinity-contrasted environments.

Construction of protein-protein interaction network in sulfate-reducing bacteria: Unveiling of global response to Hg

Sulfate-reducing bacteria (SRB) play pivotal roles in the biotransformation of mercury (Hg). However, unrevealed global responses of SRB to Hg have restricted our understanding of details of Hg biotransformation processes. The absence of protein-protein interaction (PPI) network under Hg stimuli has been a bottleneck of proteomic analysis for molecular mechanisms of Hg transformation. This study constructed the first comprehensive PPI network of SRB in response to Hg, encompassing 67 connected nodes, 26 independent nodes, and 121 edges, covering 93% of differentially expressed proteins from both previous studies and this study. The network suggested that proteomic changes of SRB in response to Hg occurred globally, including microbial metabolism in diverse environments, carbon metabolism, nucleic acid metabolism and translation, nucleic acid repair, transport systems, nitrogen metabolism, and methyltransferase activity, partial of which could cover the known knowledge. Antibiotic resistance was the original response revealed by this network, providing insights into of Hg biotransformation mechanisms. This study firstly provided the foundational network for a comprehensive understanding of SRB's responses to Hg, convenient for exploration of potential targets for Hg biotransformation. Furthermore, the network indicated that Hg enhances the metabolic activities and modification pathways of SRB to maintain cellular activities, shedding light on the influences of Hg on the carbon, nitrogen, and sulfur cycles at the cellular level.

Comparative transcriptomics elucidates the cellular responses of an aeroterrestrial zygnematophyte to UV radiation

The zygnematophytes are the closest relatives of land plants and comprise several lineages that adapted to a life on land. Species of the genus Serritaenia form colorful, mucilaginous capsules, which surround the cells and block harmful solar radiation, one of the major terrestrial stressors. In eukaryotic algae, this 'sunscreen mucilage' represents a unique photoprotective strategy, whose induction and chemical background are unknown. We generated a de novo transcriptome of Serritaenia testaceovaginata and studied its gene regulation under moderate UV radiation (UVR) that triggers sunscreen mucilage under experimental conditions. UVR induced the repair of DNA and the photosynthetic apparatus as well as the synthesis of aromatic specialized metabolites. Specifically, we observed pronounced expressional changes in the production of aromatic amino acids, phenylpropanoid biosynthesis genes, potential cross-membrane transporters of phenolics, and extracellular, oxidative enzymes. Interestingly, the most up-regulated enzyme was a secreted class III peroxidase, whose embryophyte homologs are involved in apoplastic lignin formation. Overall, our findings reveal a conserved, plant-like UVR perception system (UVR8 and downstream factors) in zygnematophyte algae and point to a polyphenolic origin of the sunscreen pigment of Serritaenia, whose synthesis might be extracellular and oxidative, resembling that of plant lignins.

Discovery of a novel symbiotic lineage associated with a hematophagous leech from the genus Haementeria

Similarly to other strict blood feeders, leeches from the Haementeria genus (Hirudinida: Glossiphoniidae) have established a symbiotic association with bacteria harbored intracellularly in esophageal bacteriomes. Previous genome sequence analyses of these endosymbionts revealed co-divergence with their hosts, a strong genome reduction, and a simplified metabolism largely dedicated to the production of B vitamins, which are nutrients lacking from a blood diet. 'Candidatus Providencia siddallii' has been identified as the obligate nutritional endosymbiont of a monophyletic clade of Mexican and South American Haementeria spp. However, the Haementeria genus includes a sister clade of congeners from Central and South America, where the presence or absence of the aforementioned symbiont taxon remains unknown. In this work, we report on a novel bacterial endosymbiont found in a representative from this Haementeria clade. We found that this symbiont lineage has evolved from within the Pluralibacter genus, known mainly from clinical but also environmental strains. Similarly to Ca. Providencia siddallii, the Haementeria-associated Pluralibacter symbiont displays clear signs of genome reduction, accompanied by an A+T-biased sequence composition. Genomic analysis of its metabolic potential revealed a retention of pathways related to B vitamin biosynthesis, supporting its role as a nutritional endosymbiont. Finally, comparative genomics of both Haementeria symbiont lineages suggests that an ancient Providencia symbiont was likely replaced by the novel Pluralibacter one, thus constituting the first reported case of nutritional symbiont replacement in a leech without morphological changes in the bacteriome.

IMPORTANCE

Obligate symbiotic associations with a nutritional base have likely evolved more than once in strict blood-feeding leeches. Unlike those symbioses found in hematophagous arthropods, the nature, identity, and evolutionary history of these remains poorly studied. In this work, we further explored obligate nutritional associations between Haementeria leeches and their microbial symbionts, which led to the unexpected discovery of a novel symbiosis with a member of the Pluralibacter genus. When compared to Providencia siddallii, an obligate nutritional symbiont of other Haementeria leeches, this novel bacterial symbiont shows convergent retention of the metabolic pathways involved in B vitamin biosynthesis. Moreover, the genomic characteristics of this Pluralibacter symbiont suggest a more recent association than that of Pr. siddallii and Haementeria. We conclude that the once-thought stable associations between blood-feeding Glossiphoniidae and their symbionts (i.e., one bacteriome structure, one symbiont lineage) can break down, mirroring symbiont turnover observed in various arthropod lineages.

Convergent reductive evolution of cyanobacteria in symbiosis with Dinophysiales dinoflagellates

The diversity of marine cyanobacteria has been extensively studied due to their vital roles in ocean primary production. However, little is understood about the diversity of cyanobacterial species involved in symbiotic relationships. In this study, we successfully sequenced the complete genome of a cyanobacterium in symbiosis with Citharistes regius, a dinoflagellate species thriving in the open ocean. A phylogenomic analysis revealed that the cyanobacterium (CregCyn) belongs to the marine picocyanobacterial lineage, akin to another cyanobacterial symbiont (OmCyn) of a different dinoflagellate closely related to Citharistes. Nevertheless, these two symbionts are representing distinct lineages, suggesting independent origins of their symbiotic lifestyles. Despite the distinct origins, the genome analyses of CregCyn revealed shared characteristics with OmCyn, including an obligate symbiotic relationship with the host dinoflagellates and a degree of genome reduction. In contrast, a detailed analysis of genome subregions unveiled that the CregCyn genome carries genomic islands that are not found in the OmCyn genome. The presence of the genomic islands implies that exogenous genes have been integrated into the CregCyn genome at some point in its evolution. This study contributes to our understanding of the complex history of the symbiosis between dinoflagellates and cyanobacteria, as well as the genomic diversity of marine picocyanobacteria.

Characterization of Sf9 cell clones with differential susceptibilities to Sf-rhabdovirus X+3.7 and Sf-rhabdovirus X- replication

Our laboratory previously discovered a novel rhabdovirus in the Spodoptera frugiperda Sf9 insect cell line that was designated as Sf-rhabdovirus. Using limiting dilution, this cell line was found to be a mixed population of cells infected by Sf-rhabdovirus variants containing either the full length X accessory gene with a 3.7 kb internal duplication (designated as Sf-rhabdovirus X+3.7) or lacking the duplication and part of the X gene (designated as Sf-rhabdovirus X-), and cells that were negative for Sf-rhabdovirus. In this paper, we found that the Sf-rhabdovirus negative cell clones had sub-populations with different susceptibilities to the replication of Sf-rhabdovirus X+3.7 and X- variants: cell clone Sf9-13F12 was more sensitive to replication by both virus variants compared to Sf9-3003; moreover, Sf9-3003 showed more resistance to X+3.7 replication than to X- replication. RNA-Seq analysis indicated significant differentially expressed genes in the Sf9-13F12 and Sf9-3003 cell clones further supporting that distinct sub-populations of virus-negative cells co-exist in the parent Sf9 cell line.

Whole genome resequencing unveils low-temperature stress tolerance specific genomic variations in jute (Corchorus sp.)

Jute (Corchorus sp.), a commercially important and eco-friendly crop, is widely cultivated in Bangladesh, India, and China. Some varieties of this tropical plant such as the Corchorus olitorius. Variety accession no. 2015 (acc. 2015) has been found to be low-temperature tolerant. The current study was designed to explore the genome-wide variations present in the tolerant plant acc. 2015 in comparison to the sensitive farmer popular variety Corchorus olitorius var. O9897 using the whole genome resequencing technique. Among different variations, intergenic Single Nucleotide Polymorphism (SNPs) and Insertion-Deletion (InDels) were found in the highest percentage whereas approximately 3% SNPs and 2% InDels were found in exonic regions in both plants. Gene enrichment analysis indicated the presence of acc. 2015 specific SNPs in the genes encoding peroxidase, ER lumen protein retaining receptor, and hexosyltransferase involved in stress response (GO:0006950) which were not present in sensitive variety O9897. Besides, distinctive copy number variation regions (CNVRs) comprising 120 gene loci were found in acc. 2015 with a gain of function from multiple copy numbers but absent in O9897. Gene ontology analysis revealed these gene loci to possess different receptors like kinases, helicases, phosphatases, transcription factors especially Myb transcription factors, regulatory proteins containing different binding domains, annexin, laccase, acyl carrier protein, potassium transporter, and vesicular transporter proteins that are responsible for low temperature induced adaptation pathways in plants. This work of identifying genomic variations linked to cold stress tolerance traits will help to develop successful markers that will pave the way to develop genetically modified cold-resistant jute lines for year-round cultivation to meet the demand for a sustainable fiber crop economy.

Phosphorylation of 399S at CsHsp70 of Cymbidium sinense is essential to maintain chlorophyll stability

The Chinese orchids symbolise nobility and gentility in China, and the variation of leaf color makes Cymbidium sinense more diversified and valuable. However, its color variations especially at the protein level still remain largely unexplored. In this study, the proteomics and phosphoproteomics of Cymbidium sinense leaf color variation mutants were studied. A total of 1059 differentially abundant proteins (DAPs) and 1127 differentially abundant phosphorylation sites belonging to 644 phosphoproteins (DAPPs) were identified in the yellow section of leaf variegation mutant of Cymbidium sinense (MY) compared with the green section (MG). Moreover, 349 co-expressing proteins were found in both omics' datasets, while only 26 proteins showed the same expression patterns in the two omics. The interaction network analysis of kinases and phosphatases showed that DAPs and DAPPs in photosynthesis, response to hormones, pigment metabolic process, phosphorylation, glucose metabolic process, and dephosphorylation might contribute to leaf color variation. The abundance of 28 Hsps and 28 phosphorylation sites belonging to 10 Hsps showed significant differences between MG and MY. CsHsp70 was selected to explore the function in Cymbidium sinense leaf variegation. The results showed CsHsp70 is essential for maintaining photosynthetic pigment content and the 399S phosphorylation site is crucial to the function of CsHsp70. Collectively, our findings construct a comprehensive coverage of protein and protein phosphorylation in leaf variegation of C. sinense, providing valuable insights into its formation mechanisms.

Detection of serum major histocompatibility complex I (HLA-1) and β2-microglobulin (β2M) in pre-eclampsia using isobaric tags for relative and absolute quantitation (iTRAQ)

OBJECTIVE

The purpose of this preliminary investigation into the pathogenesis of pre-eclampsia was to screen the differential proteins in the serum of pregnant women with normal pregnancy and early-onset pre-eclampsia using isobaric tags for relative and absolute quantitation (iTRAQ), so as to identify serum biomarkers for the early diagnosis of pre-eclampsia.

METHODS

We examined the peripheral serum of 58 normal pregnant women and 42 pregnant women with early-onset pre-eclampsia using iTRAQ; the differentially expressed proteins were screened for bioinformatics analysis; and the expression of candidate proteins human leukocyte antigen-1 (HLA-1) and β2-microglobulin (β2M) in placental tissues was detected using western blot.

RESULTS

We identified a total of 63 differential proteins in the serum of patients from the normal control group and the pre-eclampsia group, and this included 24 up-regulated proteins and 39 down-regulated proteins. The western blot results of placental tissue showed reduced HLA-1 expression (1.12 ± 0.23) in the placenta in the pre-eclampsia group as compared with the normal control group (1.34 ± 0.22). Consistent with the results observed in the serum, β2M in the placenta in the pre-eclampsia group was significantly elevated (1.05 ± 0.47) in comparison with the normal group (0.75 ± 0.33) (P < 0.05).

CONCLUSION

In this study, we found that iTRAQ technology was useful for identifying differentially expressed proteins in the peripheral serum of pregnant women with pre-eclampsia, and that HLA-1 and β2M, which may be involved in the occurrence of pre-eclampsia, show promise as predictive markers of pre-eclampsia.

Different oxygen affinities of methanotrophs and Comammox Nitrospira inform an electrically induced symbiosis for nitrogen loss

Aerobic methanotrophs establish a symbiotic association with denitrifiers to facilitate the process of aerobic methane oxidation coupled with denitrification (AME-D). However, the symbiosis has been frequently observed in hypoxic conditions continuing to pose an enigma. The present study has firstly characterized an electrically induced symbiosis primarily governed by Methylosarcina and Hyphomicrobium for the AME-D process in a hypoxic niche caused by Comammox Nitrospira. The kinetic analysis revealed that Comammox Nitrospira exhibited a higher apparent oxygen affinity compared to Methylosarcina. While the coexistence of comammox and AME-D resulted in an increase in methane oxidation and nitrogen loss rates, from 0.82 ± 0.10 to 1.72 ± 0.09 mmol CH4 d-1 and from 0.59 ± 0.04 to 1.30 ± 0.15 mmol N2 d-1, respectively. Furthermore, the constructed microbial fuel cells demonstrated a pronounced dependence of the biocurrents on AME-D due to oxygen competition, suggesting the involvement of direct interspecies electron transfer in the AME-D process under hypoxic conditions. Metagenomic and metatranscriptomic analysis revealed that Methylosarcina efficiently oxidized methane to formaldehyde, subsequently generating abundant NAD(P)H for nitrate reduction by Hyphomicrobium through the dissimilatory RuMP pathway, leading to CO2 production. This study challenges the conventional understanding of survival mechanism employed by AME-D symbionts, thereby contributing to the characterization responsible for limiting methane emissions and promoting nitrogen removal in hypoxic regions.

Genetic dissection of Brassica napus seed vigor after aging

KEY MESSAGE

Stable and novel QTLs that affect seed vigor under different storage durations were discovered, and BnaOLE4, located in the interval of cqSW-C2-3, increased seed vigor after aging. Seed vigor is an important trait in crop breeding; however, the underlying molecular regulatory mechanisms governing this trait in rapeseed remain largely unknown. In the present study, vigor-related traits were analyzed in seeds from a doubled haploid (DH) rapeseed (Brassica napus) population grown in 2 different environments using seeds stored for 7, 5, and 3 years under natural storage conditions. A total of 229 quantitative trait loci (QTLs) were identified and were found to explain 3.78%-17.22% of the phenotypic variance for seed vigor-related traits after aging. We further demonstrated that seed vigor-related traits were positively correlated with oil content (OC) but negatively correlated with unsaturated fatty acids (FAs). Some pleiotropic QTLs that collectively regulate OC, FAs, and seed vigor, such as uq.A8, uq.A3-2, uq.A9-2, and uq.C3-1, were identified. The transcriptomic results from extreme pools of DH lines with distinct seed vigor phenotypes during accelerated aging revealed that various biological pathways and metabolic processes (such as glutathione metabolism and reactive oxygen species) were involved in seed vigor. Through integration of QTL analysis and RNA-Seq, a regulatory network for the control of seed vigor was constructed. Importantly, a candidate (BnaOLE4) from cqSW-C2-3 was selected for functional analysis, and transgenic lines overexpressing BnaOLE4 showed increased seed vigor after artificial aging. Collectively, these results provide novel information on QTL and potential candidate genes for molecular breeding for improved seed storability.

Foliar mycobiome remains unaltered under urban air-pollution but differentially express stress-related genes

Air pollution caused by tropospheric ozone contributes to the decline of forest ecosystems; for instance, sacred fir, Abies religiosa (Kunth) Schltdl. & Cham. forests in the peri-urban region of Mexico City. Individual trees within these forests exhibit variation in their response to ozone exposure, including the severity of visible symptoms in needles. Using RNA-Seq metatranscriptomic data and ITS2 metabarcoding, we investigated whether symptom variation correlates with the taxonomic and functional composition of fungal mycobiomes from needles collected in this highly polluted area in the surroundings of Mexico City. Our findings indicate that ozone-related symptoms do not significantly correlate with changes in the taxonomic composition of fungal mycobiomes. However, genes coding for 30 putative proteins were differentially expressed in the mycobiome of asymptomatic needles, including eight genes previously associated with resistance to oxidative stress. These results suggest that fungal communities likely play a role in mitigating the oxidative burst caused by tropospheric ozone in sacred fir. Our study illustrates the feasibility of using RNA-Seq data, accessible from global sequence repositories, for the characterization of fungal communities associated with plant tissues, including their gene expression.

Subtractive Proteomics and Reverse-Vaccinology Approaches for Novel Drug Target Identification and Chimeric Vaccine Development against Bartonella henselae Strain Houston-1

Bartonella henselae is a Gram-negative bacterium causing a variety of clinical symptoms, ranging from cat-scratch disease to severe systemic infections, and it is primarily transmitted by infected fleas. Its status as an emerging zoonotic pathogen and its capacity to persist within host erythrocytes and endothelial cells emphasize its clinical significance. Despite progress in understanding its pathogenesis, limited knowledge exists about the virulence factors and regulatory mechanisms specific to the B. henselae strain Houston-1. Exploring these aspects is crucial for targeted therapeutic strategies against this versatile pathogen. Using reverse-vaccinology-based subtractive proteomics, this research aimed to identify the most antigenic proteins for formulating a multi-epitope vaccine against the B. henselae strain Houston-1. One crucial virulent and antigenic protein, the PAS domain-containing sensor histidine kinase protein, was identified. Subsequently, the identification of B-cell and T-cell epitopes for the specified protein was carried out and the evaluated epitopes were checked for their antigenicity, allergenicity, solubility, MHC binding capability, and toxicity. The filtered epitopes were merged using linkers and an adjuvant to create a multi-epitope vaccine construct. The structure was then refined, with 92.3% of amino acids falling within the allowed regions. Docking of the human receptor (TLR4) with the vaccine construct was performed and demonstrated a binding energy of -1047.2 Kcal/mol with more interactions. Molecular dynamic simulations confirmed the stability of this docked complex, emphasizing the conformation and interactions between the molecules. Further experimental validation is necessary to evaluate its effectiveness against B. henselae.

Subtractive genomics study for the identification of therapeutic targets against Cronobacter sakazakii: A threat to infants

Cronobacter sakazakii is an opportunistic pathogen that has been associated with severe infection in neonates such as necrotizing enterocolitis (NEC), neonatal meningitis, and bacteremia. This pathogen can survive in a relatively dry environment, especially in powdered infant formula (PIF). Unfortunately, conventional drugs that were once effective against C. sakazakii are gradually losing their efficacy due to rising antibiotic resistance. In this study, a subtractive genomic approach was followed in order to identify potential therapeutic targets in the pathogen. The whole proteome of the pathogen was filtered through a step-by-step process, which involved removing paralogous proteins, human homologs, sequences that are less essential for survival, proteins with shared metabolic pathways, and proteins that are located in cells other than the cytoplasmic membrane. As a result, nine novel drug targets were identified. Further, the analysis also unveiled that the FDA-approved drug Terbinafine can be repurposed against the Glutathione/l-cysteine transport system ATP-binding/permease protein CydC of C. sakazakii. Moreover, molecular docking and dynamics studies of Terbinafine and CydC suggested that this drug can be used to treat C. sakazakii infection in neonates. However, for clinical purposes further in vitro and in vivo studies are necessary.

Exploring optimal drug targets through subtractive proteomics analysis and pangenomic insights for tailored drug design in tuberculosis

Tuberculosis (TB), caused by Mycobacterium tuberculosis, ranks among the top causes of global human mortality, as reported by the World Health Organization's 2022 TB report. The prevalence of M. tuberculosis strains that are multiple and extensive-drug resistant represents a significant barrier to TB eradication. Fortunately, having many completely sequenced M. tuberculosis genomes available has made it possible to investigate the species pangenome, conduct a pan-phylogenetic investigation, and find potential new drug targets. The 442 complete genome dataset was used to estimate the pangenome of M. tuberculosis. This study involved phylogenomic classification and in-depth analyses. Sequential filters were applied to the conserved core genome containing 2754 proteins. These filters assessed non-human homology, virulence, essentiality, physiochemical properties, and pathway analysis. Through these intensive filtering approaches, promising broad-spectrum therapeutic targets were identified. These targets were docked with FDA-approved compounds readily available on the ZINC database. Selected highly ranked ligands with inhibitory potential include dihydroergotamine and abiraterone acetate. The effectiveness of the ligands has been supported by molecular dynamics simulation of the ligand-protein complexes, instilling optimism that the identified lead compounds may serve as a robust basis for the development of safe and efficient drugs for TB treatment, subject to further lead optimization and subsequent experimental validation.

Colonization of vancomycin-resistant Enterococcus faecium in human-derived colonic epithelium: unraveling the transcriptional dynamics of host-enterococcal interactions

Enterococcus faecium is an opportunistic pathogen able to colonize the intestines of hospitalized patients. This initial colonization is an important step in the downstream pathogenesis, which includes outgrowth of the intestinal microbiota and potential infection of the host. The impact of intestinal overgrowth on host-enterococcal interactions is not well understood. We therefore applied a RNAseq approach in order to unravel the transcriptional dynamics of E. faecium upon co-culturing with human derived colonic epithelium. Co-cultures of colonic epithelium with a hospital-associated vancomycin resistant (vanA-type) E. faecium (VRE) showed that VRE resided on top of the colonic epithelium when analyzed by microscopy. RNAseq revealed that exposure to the colonic epithelium resulted in upregulation of 238 VRE genes compared to the control condition, including genes implicated in pili expression, conjugation (plasmid_2), genes related to sugar uptake, and biofilm formation (chromosome). In total, 260 were downregulated, including the vanA operon located on plasmid_3. Pathway analysis revealed an overall switch in metabolism to amino acid scavenging and reduction. In summary, our study demonstrates that co-culturing of VRE with human colonic epithelium promotes an elaborate gene response in VRE, enhancing our insight in host-E. faecium interactions, which might facilitate the design of novel anti-infectivity strategies.

Integrated genome-transcriptome analysis unveiled the mechanism of Debaryomyces hansenii-mediated arsenic stress amelioration in rice

Globally, rice is becoming more vulnerable to arsenic (As) pollution, posing a serious threat to public food safety. Previously Debaryomyces hansenii was found to reduce grain As content of rice. To better understand the underlying mechanism, we performed a genome analysis to identify the key genes in D. hansenii responsible for As tolerance and plant growth promotion. Notably, genes related to As resistance (ARR, Ycf1, and Yap) were observed in the genome of D. hansenii. The presence of auxin pathway and glutathione metabolism-related genes may explain the plant growth-promoting potential and As tolerance mechanism of this novel yeast strain. The genome annotation of D. hansenii indicated that it contains a repertoire of genes encoding antioxidants, well corroborated with the in vitro studies of GST, GR, and glutathione content. In addition, the effect of D. hansenii on gene expression profiling of rice plants under As stress was also examined. The Kyoto Encyclopedia of Genes and Genomes (KEGG) database revealed 307 genes, annotated in D. hansenii-treated rice, related to metabolic pathways (184), photosynthesis (12), glutathione (10), tryptophan (4), and biosynthesis of secondary metabolite (117). Higher expression of regulatory elements like AUX/IAA and WRKY transcription factors (TFs), and defense-responsive genes dismutases, catalases, peroxiredoxin, and glutaredoxins during D. hansenii+As exposure was also observed. Combined analysis revealed that D. hansenii genes are contributing to stress mitigation in rice by supporting plant growth and As-tolerance. The study lays the foundation to develop yeast as a beneficial biofertilizer for As-prone areas.

PacBio Hi-Fi genome assembly of Sipha maydis, a model for the study of multipartite mutualism in insects

Dependence on multiple nutritional endosymbionts has evolved repeatedly in insects feeding on unbalanced diets. However, reference genomes for species hosting multi-symbiotic nutritional systems are lacking, even though they are essential for deciphering the processes governing cooperative life between insects and anatomically integrated symbionts. The cereal aphid Sipha maydis is a promising model for addressing these issues, as it has evolved a nutritional dependence on two bacterial endosymbionts that complement each other. In this study, we used PacBio High fidelity (HiFi) long-read sequencing to generate a highly contiguous genome assembly of S. maydis with a length of 410 Mb, 3,570 contigs with a contig N50 length of 187 kb, and BUSCO completeness of 95.5%. We identified 117 Mb of repetitive sequences, accounting for 29% of the genome assembly, and predicted 24,453 protein-coding genes, of which 2,541 were predicted enzymes included in an integrated metabolic network with the two aphid-associated endosymbionts. These resources provide valuable genetic and metabolic information for understanding the evolution and functioning of multi-symbiotic systems in insects.

Genomic and Transcriptomic Diversification of Flagellin Genes Provides Insight into Environmental Adaptation and Phylogeographic Characteristics in Aeromonas hydrophila

Aeromonas hydrophila is an opportunistic motile pathogen with a broad host range, infecting both terrestrial and aquatic animals. Environmental and geographical conditions exert selective pressure on both geno- and phenotypes of pathogens. Flagellin, directly exposed to external environments and containing important immunogenic epitopes, may display significant variability in response to external conditions. In this study, we conducted a comparative analysis of ~  150 A. hydrophila genomes, leading to the identification of six subunits of the flagellin gene (fla-1 to fla-4, flaA, and flaB). Individual strains harbored different composition of flagellin subunits and copies. The composition of subunits showed distinct patterns depending on environmental sources. Strains from aquatic environments were mainly comprised of fla-1 to fla-4 subunits, while terrestrial strains predominated in groups harboring flaA and flaB subunits. Each flagellin showed varying levels of expression, with flaA and flaB demonstrating significantly higher expression compared to others. One of the chemotaxis pathways that control flagellin movement through a two-component system was significantly upregulated in flaA(+ 1)/flaB(+ 1) group, whereas flaA and flaB showed different transcriptomic expressions. The genes positively correlated with flaA expression were relevant to biofilm formation and bacterial chemotaxis, but flaB showed a negative correlation with the genes in ABC transporters and quorum sensing pathway. However, the expression patterns of fla-2 to fla-4 were identical. This suggests various types of flagellin subunits may have different biological functions. The composition and expression levels of flagellin subunits could provide valuable insights into the adaptation of A. hydrophila and the differences among strains in response to various external environments.

Deinococcus arenicola sp. nov., a novel radiation-resistant bacterium isolated from sandy soil in Antarctica

A Gram-stain-positive, aerobic, non-mobile and spherical strain, designated ZS9-10T, belonging to the genus Deinococcus was isolated from soil sampled at the Chinese Zhong Shan Station, Antarctica. Growth was observed in the presence of 0-4 % (w/v) NaCl, at pH 7.0-8.0 and at 4-25 °C. Phylogenetic analysis based on 16S rRNA gene sequences showed that strain ZS9-10T formed a lineage in the genus Deinococcus. It exhibited highest sequence similarity (97.4 %) to Deinococcus marmoris DSM 12784T. The major phospholipids of ZS9-10T were unidentified phosphoglycolipid, unidentified glycolipids and unidentified lipids. The major fatty acids were summed feature 3 (C16 : 1  ω7c and/or C16 : 1  ω6c), C16 : 0 and C16 : 1  ω7c. MK-8 was the predominant respiratory quinone. The digital DNA-DNA hybridization and average nucleotide identity values between strain ZS9-10T and its close relative D. marmoris DSM 12784T were 27.4 and 83.9 %, respectively. Based on phenotypic, phylogenetic and genotypic data, a novel species, named Deinococcus arenicola sp. nov., is proposed. The type strain iis ZS9-10T (=CCTCC AB 2019392T=KCTC43192T).

Metabolic relationships between marine red algae and algae-associated bacteria

Mutualistic interactions between marine phototrophs and associated bacteria are an important strategy for their successful survival in the ocean, but little is known about their metabolic relationships. Here, bacterial communities in the algal sphere (AS) and bulk solution (BS) of nine marine red algal cultures were analyzed, and Roseibium and Phycisphaera were identified significantly more abundantly in AS than in BS. The metabolic features of Roseibium RMAR6-6 (isolated and genome-sequenced), Phycisphaera MAG 12 (obtained by metagenomic sequencing), and a marine red alga, Porphyridium purpureum CCMP1328 (from GenBank), were analyzed bioinformatically. RMAR6-6 has the genetic capability to fix nitrogen and produce B vitamins (B1, B2, B5, B6, B9, and B12), bacterioferritin, dimethylsulfoniopropionate (DMSP), and phenylacetate that may enhance algal growth, whereas MAG 12 may have a limited metabolic capability, not producing vitamins B9 and B12, DMSP, phenylacetate, and siderophores, but with the ability to produce bacitracin, possibly modulating algal microbiome. P. purpureum CCMP1328 lacks the genetic capability to fix nitrogen and produce vitamin B12, DMSP, phenylacetate, and siderophore. It was shown that the nitrogen-fixing ability of RMAR6-6 promoted the growth of P. purpureum, and DMSP reduced the oxidative stress of P. purpureum. The metabolic interactions between strain RMAR6-6 and P. purpureum CCMP1328 were also investigated by the transcriptomic analyses of their monoculture and co-culture. Taken together, potential metabolic relationships between Roseibium and P. purpureum were proposed. This study provides a better understanding of the metabolic relationships between marine algae and algae-associated bacteria for successful growth.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-024-00227-z.

Phylogenomic analysis of metagenome-assembled genomes indicates new taxa in the order Spirochaetales and proposal of Thalassospirochaeta sargassi gen. nov. sp. nov. from seaweeds

Six metagenome-assembled genomes (JB008Ts, JB007, JB015, JB003, JB004, and JB002) belonging to the order Spirochaetales were generated from seaweed samples collected from the Gulf of Mannar, India. The binned genomes JB008Ts and JB007 shared highest 16S rRNA gene identity of 94.9 % and 92.2-93.4 %, respectively with uncultivated Spirochaetaceae family members, and < 90 % identity with Marispirochaeta aestuari JC444T. While, the bin JB015 showed 99.1 % identity with Pleomorphochaeta naphthae SEBR 4209T. The phylogenomic and 16S rRNA gene-based phylogenetic analysis of the binned genomes JB007 and JB008Ts confirmed that these members belong to the family Spirochaetaceae and bins JB015, JB002, JB003, and JB004 belong to the genus Pleomorphochaeta within the family Sphaerochaetaceae. The AAI values of the binned genomes JB007 and JB008Ts compared to other members of the Spirochaetaceae family were between 53.9- 56.8 % and 53.8-57.1 %, respectively. Furthermore, the comparison of ANI, dDDH, and POCP metrics of the binned genomes JB007 and JB008Ts, both among themselves and with the members of Spirochaetaceae, was also below the suggested thresholds for genera delineation. Consequently, the binned genome JB008Ts is proposed as a new genus according to the guidelines of code of nomenclature of prokaryotes described from sequence data (SeqCode) with the name Thalassospirochaeta sargassi gen. nov. sp. nov., in the family Spirochaetaceae while the bin JB007 could not be proposed as novel taxa due to low-quality estimates. The bin JB015 and its additional genomes form a distinct clade, but their taxonomic status remains ambiguous due to the absence of genomic evidence from other Pleomorphochaeta members.

De novo genome assembly and population genomics of a shrub tree Barthea barthei (Hance) krass provide insights into the adaptive color variations

Flower color is a classic example of an ecologically important trait under selection in plants. Understanding the genetic mechanisms underlying shifts in flower color can provide key insights into ecological speciation. In this study, we investigated the genetic basis of flower color divergence in Barthea barthei, a shrub tree species exhibiting natural variation in flower color. We assembled a high-quality genome assembly for B. barthei with a contig N50 of 2.39 Mb and a scaffold N50 of 16.21 Mb. The assembly was annotated with 46,430 protein-coding genes and 1,560 non-coding RNAs. Genome synteny analysis revealed two recent tetraploidization events in B. barthei, estimated to have occurred at approximately 17 and 63 million years ago. These tetraploidization events resulted in massive duplicated gene content, with over 70% of genes retained in collinear blocks. Gene family members of the core regulators of the MBW complex were significantly expanded in B. barthei compared to Arabidopsis, suggesting that these duplications may have provided raw genetic material for the evolution of novel regulatory interactions and the diversification of anthocyanin pigmentation. Transcriptome profiling of B. barthei flowers revealed differential expression of 9 transcription factors related to anthocyanin biosynthesis between the two ecotypes. Six of these differentially expressed transcription factors were identified as high-confidence candidates for adaptive evolution based on positive selection signals. This study provides insights into the genetic basis of flower color divergence and the evolutionary mechanisms underlying ecological adaptation in plants.

Microbacterium horticulturae sp. nov., a novel actinobacterium isolated from flowerpot soil

Strain CJN36-1NT, a Gram-stain-positive, non-flagellated, strictly aerobic and short rod-shaped bacterium, was isolated from flowerpot soil sampled in the Jeonju region of the Republic of Korea. Based on 16S rRNA gene sequences and the resulting phylogenetic tree, the strain belonged to the genus Microbacterium. Strain CJN36-1NT contained a chromosome of 3.6 Mbp with a G+C content of 68.5 mol%. The strain grew at 10-37 °C (optimally at 28 °C), at pH 5.0-8.0 (optimally at pH 8.0), and in the presence of 0-7 % NaCl (w/v; optimally with 0 % NaCl). Digital DNA-DNA hybridization, average nucleotide identity and average amino acid identity values between strain CJN36-1NT and its closest related species, Microbacterium protaetiae DFW100M-13T, were 82.0, 81.2, and 23.2 %, respectively. We propose naming this novel species Microbacterium horticulturae sp. nov., with CJN36-1NT (=KACC 23027T=NBRC 116065T) as the type strain.

Identification and clinicopathological analysis of potential p73-regulated biomarkers in colorectal cancer via integrative bioinformatics

This study aims to decipher crucial biomarkers regulated by p73 for the early detection of colorectal cancer (CRC) by employing a combination of integrative bioinformatics and expression profiling techniques. The transcriptome profile of HCT116 cell line p53- / - p73+ / + and p53- / - p73 knockdown was performed to identify differentially expressed genes (DEGs). This was corroborated with three CRC tissue expression datasets available in Gene Expression Omnibus. Further analysis involved KEGG and Gene ontology to elucidate the functional roles of DEGs. The protein-protein interaction (PPI) network was constructed using Cytoscape to identify hub genes. Kaplan-Meier (KM) plots along with GEPIA and UALCAN database analysis provided the insights into the prognostic and diagnostic significance of these hub genes. Machine/deep learning algorithms were employed to perform TNM-stage classification. Transcriptome profiling revealed 1289 upregulated and 1897 downregulated genes. When intersected with employed CRC datasets, 284 DEGs were obtained. Comprehensive analysis using gene ontology and KEGG revealed enrichment of the DEGs in metabolic process, fatty acid biosynthesis, etc. The PPI network constructed using these 284 genes assisted in identifying 20 hub genes. Kaplan-Meier, GEPIA, and UALCAN analyses uncovered the clinicopathological relevance of these hub genes. Conclusively, the deep learning model achieved TNM-stage classification accuracy of 0.78 and 0.75 using 284 DEGs and 20 hub genes, respectively. The study represents a pioneer endeavor amalgamating transcriptomics, publicly available tissue datasets, and machine learning to unveil key CRC-associated genes. These genes are found relevant regarding the patients' prognosis and diagnosis. The unveiled biomarkers exhibit robustness in TNM-stage prediction, thereby laying the foundation for future clinical applications and therapeutic interventions in CRC management.

Phytoplankton-derived polysaccharides and microbial peptidoglycans are key nutrients for deep-sea microbes in the Mariana Trench

BACKGROUND

The deep sea represents the largest marine ecosystem, driving global-scale biogeochemical cycles. Microorganisms are the most abundant biological entities and play a vital role in the cycling of organic matter in such ecosystems. The primary food source for abyssal biota is the sedimentation of particulate organic polymers. However, our knowledge of the specific biopolymers available to deep-sea microbes remains largely incomplete. One crucial rate-limiting step in organic matter cycling is the depolymerization of particulate organic polymers facilitated by extracellular enzymes (EEs). Therefore, the investigation of active EEs and the microbes responsible for their production is a top priority to better understand the key nutrient sources for deep-sea microbes.

RESULTS

In this study, we conducted analyses of extracellular enzymatic activities (EEAs), metagenomics, and metatranscriptomics from seawater samples of 50-9305 m from the Mariana Trench. While a diverse array of microbial groups was identified throughout the water column, only a few exhibited high levels of transcriptional activities. Notably, microbial populations actively transcribing EE genes involved in biopolymer processing in the abyssopelagic (4700 m) and hadopelagic zones (9305 m) were primarily associated with the class Actinobacteria. These microbes actively transcribed genes coding for enzymes such as cutinase, laccase, and xyloglucanase which are capable of degrading phytoplankton polysaccharides as well as GH23 peptidoglycan lyases and M23 peptidases which have the capacity to break down peptidoglycan. Consequently, corresponding enzyme activities including glycosidases, esterase, and peptidases can be detected in the deep ocean. Furthermore, cell-specific EEAs increased at 9305 m compared to 4700 m, indicating extracellular enzymes play a more significant role in nutrient cycling in the deeper regions of the Mariana Trench.

CONCLUSIONS

Transcriptomic analyses have shed light on the predominant microbial population actively participating in organic matter cycling in the deep-sea environment of the Mariana Trench. The categories of active EEs suggest that the complex phytoplankton polysaccharides (e.g., cutin, lignin, and hemicellulose) and microbial peptidoglycans serve as the primary nutrient sources available to deep-sea microbes. The high cell-specific EEA observed in the hadal zone underscores the robust polymer-degrading capacities of hadal microbes even in the face of the challenging conditions they encounter in this extreme environment. These findings provide valuable new insights into the sources of nutrition, the key microbes, and the EEs crucial for biopolymer degradation in the deep seawater of the Mariana Trench. Video Abstract.

Genomic characterization and probiotic potential assessment of an exopolysaccharide-producing strain Pediococcus pentosaceus LL-07 isolated from fermented meat

BACKGROUND

The genomic information available for Pediococcus pentosaceus is primarily derived from fermented fruits and vegetables, with less information available from fermented meat. P. pentosaceus LL-07, a strain isolated from fermented meat, has the capability of producing exopolysaccharides (EPS). To assess the probiotic attributes of P. pentosaceus LL-07, we conducted whole-genome sequencing (WGS) using the PacBio SequelIIe and Illumina MiSeq platforms, followed by in vitro experiments to explore its probiotic potential.

RESULTS

The genome size of P. pentosaceus LL-07 is 1,782,685 bp, comprising a circular chromosome and a circular plasmid. Our investigation revealed the absence of a CRISPR/Cas system. Sugar fermentation experiments demonstrated the characteristics of carbohydrate metabolism. P. pentosaceus LL-07 contains an EPS synthesis gene cluster consisting of 13 genes, which is different from the currently known gene cluster structure. NO genes associated with hemolysis or toxin synthesis were detected. Additionally, eighty-six genes related to antibiotic resistance were identified but not present in the prophage, transposon or plasmid. In vitro experiments demonstrated that P. pentosaceus LL-07 was comparable to the reference strain P. pentosaceus ATCC25745 in terms of tolerance to artificial digestive juice and bile, autoaggregation and antioxidation, and provided corresponding genomic evidence.

CONCLUSION

This study confirmed the safety and probiotic properties of P. pentosaceus LL-07 via complete genome and phenotype analysis, supporting its characterization as a potential probiotic candidate.

Anaerostipes hadrus, a butyrate-producing bacterium capable of metabolizing 5-fluorouracil

Anaerostipes hadrus (A. hadrus) is a dominant species in the human gut microbiota and considered a beneficial bacterium for producing probiotic butyrate. However, recent studies have suggested that A. hadrus may negatively affect the host through synthesizing fatty acid and metabolizing the anticancer drug 5-fluorouracil, indicating that the impact of A. hadrus is complex and unclear. Therefore, comprehensive genomic studies on A. hadrus need to be performed. We integrated 527 high-quality public A. hadrus genomes and five distinct metagenomic cohorts. We analyzed these data using the approaches of comparative genomics, metagenomics, and protein structure prediction. We also performed validations with culture-based in vitro assays. We constructed the first large-scale pan-genome of A. hadrus (n = 527) and identified 5-fluorouracil metabolism genes as ubiquitous in A. hadrus genomes as butyrate-producing genes. Metagenomic analysis revealed the wide and stable distribution of A. hadrus in healthy individuals, patients with inflammatory bowel disease, and patients with colorectal cancer, with healthy individuals carrying more A. hadrus. The predicted high-quality protein structure indicated that A. hadrus might metabolize 5-fluorouracil by producing bacterial dihydropyrimidine dehydrogenase (encoded by the preTA operon). Through in vitro assays, we validated the short-chain fatty acid production and 5-fluorouracil metabolism abilities of A. hadrus. We observed for the first time that A. hadrus can convert 5-fluorouracil to α-fluoro-β-ureidopropionic acid, which may result from the combined action of the preTA operon and adjacent hydA (encoding bacterial dihydropyrimidinase). Our results offer novel understandings of A. hadrus, exceptionally functional features, and potential applications.

IMPORTANCE

This work provides new insights into the evolutionary relationships, functional characteristics, prevalence, and potential applications of Anaerostipes hadrus.

Genome-based taxonomic classification of the genus Sulfitobacter along with the proposal of a new genus Parasulfitobacter gen. nov. and exploring the gene clusters associated with sulfur oxidation

BACKGROUND

The genus Sulfitobacter, a member of the family Roseobacteraceae, is widely distributed in the ocean and is believed to play crucial roles in the global sulfur cycle. However, gene clusters associated with sulfur oxidation in genomes of the type strains of this genus have been poorly studied. Furthermore, taxonomic errors have been identified in this genus, potentially leading to significant confusion in ecological and evolutionary interpretations in subsequent studies of the genus Sulfitobacter. This study aims to investigate the taxonomic status of this genus and explore the metabolism associated with sulfur oxidation.

RESULTS

This study suggests that Sulfitobacter algicola does not belong to Sulfitobacter and should be reclassified into a novel genus, for which we propose the name Parasulfitobacter gen. nov., with Parasulfitobacter algicola comb. nov. as the type species. Additionally, enzymes involved in the sulfur oxidation process, such as the sulfur oxidization (Sox) system, the disulfide reductase protein family, and the sulfite dehydrogenase (SoeABC), were identified in almost all Sulfitobacter species. This finding implies that the majority of Sulfitobacter species can oxidize reduced sulfur compounds. Differences in the modular organization of sox gene clusters among Sulfitobacter species were identified, along with the presence of five genes with unknown function located in some of the sox gene clusters. Lastly, this study revealed the presence of the demethylation pathway and the cleavage pathway used by many Sulfitobacter species to degrade dimethylsulfoniopropionate (DMSP). These pathways enable these bacteria to utilize DMSP as important source of sulfur and carbon or as a defence strategy.

CONCLUSIONS

Our findings contribute to interpreting the mechanism by which Sulfitobacter species participate in the global sulfur cycle. The taxonomic rearrangement of S. algicola into the novel genus Parasulfitobacter will prevent confusion in ecological and evolutionary interpretations in future studies of the genus Sulfitobacter.

Statistical Genomics Analysis of Simple Sequence Repeats from the Paphiopedilum Malipoense Transcriptome Reveals Control Knob Motifs Modulating Gene Expression

Simple sequence repeats (SSRs) are found in nonrandom distributions in genomes and are thought to impact gene expression. The distribution patterns of 48 295 SSRs of Paphiopedilum malipoense are mined and characterized based on the first full-length transcriptome and comprehensive transcriptome dataset from 12 organs. Statistical genomics analyses are used to investigate how SSRs in transcripts affect gene expression. The results demonstrate the correlations between SSR distributions, characteristics, and expression level. Nine expression-modulating motifs (expMotifs) are identified and a model is proposed to explain the effect of their key features, potency, and gene function on an intra-transcribed region scale. The expMotif-transcribed region combination is the most predominant contributor to the expression-modulating effect of SSRs, and some intra-transcribed regions are critical for this effect. Genes containing the same type of expMotif-SSR elements in the same transcribed region are likely linked in function, regulation, or evolution aspects. This study offers novel evidence to understand how SSRs regulate gene expression and provides potential regulatory elements for plant genetic engineering.

Genome-wide characterization of post-transcriptional processes related to wood formation in Dalbergia odorifera

BACKGROUND

Alternative polyadenylation (APA), alternative splicing (AS), and long non-coding RNAs (lncRNAs) play regulatory roles in post-transcriptional processes in plants. However, little is known about their involvement in xylem development in Dalbergia odorifera, a valuable rosewood species with medicinal and commercial significance. We addressed this by conducting Isoform Sequencing (Iso-Seq) using PacBio's SMRT technology and combined it with RNA-seq analysis (RNA sequencing on Illumina platform) after collecting xylem samples from the transition zone and the sapwood of D. odorifera.

RESULTS

We identified 14,938 full-length transcripts, including 9,830 novel isoforms, which has updated the D. odorifera genome annotation. Our analysis has revealed that 4,164 genes undergo APA, whereas 3,084 genes encounter AS. We have also annotated 118 lncRNAs. Furthermore, RNA-seq analysis identified 170 differential alternative splicing (DAS) events, 344 genes with differential APA site usage (DE-APA), and 6 differentially expressed lncRNAs in the transition zone when compared to the sapwood. AS, APA, and lncRNAs are differentially regulated during xylem development. Differentially expressed APA genes were enriched for terpenoid and flavonoid metabolism, indicating their role in the heartwood formation. Additionally, DE-APA genes were associated with cell wall biosynthesis and terpenoid metabolism, implying an APA's role in wood formation. A DAS gene (involved in chalcone accumulation) with a significantly greater inclusion of the last exon in the transition zone than in the sapwood was identified. We also found that differentially expressed lncRNAs targeted the genes related to terpene synthesis.

CONCLUSIONS

This study enhances our understanding of the molecular regulatory mechanisms underlying wood formation in D. odorifera, and provides valuable genetic resources and insights for its molecular-assisted breeding.

KEGG orthology prediction of bacterial proteins using natural language processing

BACKGROUND

The advent of high-throughput technologies has led to an exponential increase in uncharacterized bacterial protein sequences, surpassing the capacity of manual curation. A large number of bacterial protein sequences remain unannotated by Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology, making it necessary to use auto annotation tools. These tools are now indispensable in the biological research landscape, bridging the gap between the vastness of unannotated sequences and meaningful biological insights.

RESULTS

In this work, we propose a novel pipeline for KEGG orthology annotation of bacterial protein sequences that uses natural language processing and deep learning. To assess the effectiveness of our pipeline, we conducted evaluations using the genomes of two randomly selected species from the KEGG database. In our evaluation, we obtain competitive results on precision, recall, and F1 score, with values of 0.948, 0.947, and 0.947, respectively.

CONCLUSIONS

Our experimental results suggest that our pipeline demonstrates performance comparable to traditional methods and excels in identifying distant relatives with low sequence identity. This demonstrates the potential of our pipeline to significantly improve the accuracy and comprehensiveness of KEGG orthology annotation, thereby advancing our understanding of functional relationships within biological systems.

The ecological role of microbiome at community-, taxonomic - and genome-levels in black-odorous waters

Black-odorous waters (BOWs) are heavily polluted waters where microbial information remains elusive mechanistically. Based on gene amplicon and metagenomics sequencing, a comprehensive study was conducted to investigate the microbial communities in urban and rural BOWs. The results revealed that microbial communities' assembly in urban and rural BOWs was predominantly governed by stochastic factors at the community level. At the taxonomic level, there were 62 core species (58.48%) in water and 207 core species (44.56%) in sediment across urban and rural areas. Notably, significant differences were observed in the functional genetic composition of BOWs between urban and rural areas. Specifically, rural areas exhibited an enhanced abundance of genes involved in nitrogen fixation, Fe2+ transport, and sulfate reduction. Conversely, urban areas showed higher abundances of some genes associated with carbon fixation, nitrification and denitrification. A sulfur-centered ecological model of microbial communities was constructed by integrating data from the three levels of analysis, and 14 near-complete draft genomes were generated, representing a substantial portion of the microbial community (35.04% in rural BOWs and 29.97% in urban BOWs). This research provides significant insights into the sustainable management and preservation of aquatic ecosystems affected by BOWs.

The massive 340 megabase genome of Anisogramma anomala, a biotrophic ascomycete that causes eastern filbert blight of hazelnut

BACKGROUND

The ascomycete fungus Anisogramma anomala causes Eastern Filbert Blight (EFB) on hazelnut (Corylus spp.) trees. It is a minor disease on its native host, the American hazelnut (C. americana), but is highly destructive on the commercially important European hazelnut (C. avellana). In North America, EFB has historically limited commercial production of hazelnut to west of the Rocky Mountains. A. anomala is an obligately biotrophic fungus that has not been grown in continuous culture, rendering its study challenging. There is a 15-month latency before symptoms appear on infected hazelnut trees, and only a sexual reproductive stage has been observed. Here we report the sequencing, annotation, and characterization of its genome.

RESULTS

The genome of A. anomala was assembled into 108 scaffolds totaling 342,498,352 nt with a GC content of 34.46%. Scaffold N50 was 33.3 Mb and L50 was 5. Nineteen scaffolds with lengths over 1 Mb constituted 99% of the assembly. Telomere sequences were identified on both ends of two scaffolds and on one end of another 10 scaffolds. Flow cytometry estimated the genome size of A. anomala at 370 Mb. The genome exhibits two-speed evolution, with 93% of the assembly as AT-rich regions (32.9% GC) and the other 7% as GC-rich (57.1% GC). The AT-rich regions consist predominantly of repeats with low gene content, while 90% of predicted protein coding genes were identified in GC-rich regions. Copia-like retrotransposons accounted for more than half of the genome. Evidence of repeat-induced point mutation (RIP) was identified throughout the AT-rich regions, and two copies of the rid gene and one of dim-2, the key genes in the RIP mutation pathway, were identified in the genome. Consistent with its homothallic sexual reproduction cycle, both MAT1-1 and MAT1-2 idiomorphs were found. We identified a large suite of genes likely involved in pathogenicity, including 614 carbohydrate active enzymes, 762 secreted proteins and 165 effectors.

CONCLUSIONS

This study reveals the genomic structure, composition, and putative gene function of the important pathogen A. anomala. It provides insight into the molecular basis of the pathogen's life cycle and a solid foundation for studying EFB.

Karyotype and LTR-RTs analysis provide insights into oak genomic evolution

BACKGROUND

Whole-genome duplication and long terminal repeat retrotransposons (LTR-RTs) amplification in organisms are essential factors that affect speciation, local adaptation, and diversification of organisms. Understanding the karyotype projection and LTR-RTs amplification could contribute to untangling evolutionary history. This study compared the karyotype and LTR-RTs evolution in the genomes of eight oaks, a dominant lineage in Northern Hemisphere forests.

RESULTS

Karyotype projections showed that chromosomal evolution was relatively conservative in oaks, especially on chromosomes 1 and 7. Modern oak chromosomes formed through multiple fusions, fissions, and rearrangements after an ancestral triplication event. Species-specific chromosomal rearrangements revealed fragments preserved through natural selection and adaptive evolution. A total of 441,449 full-length LTR-RTs were identified from eight oak genomes, and the number of LTR-RTs for oaks from section Cyclobalanopsis was larger than in other sections. Recent amplification of the species-specific LTR-RTs lineages resulted in significant variation in the abundance and composition of LTR-RTs among oaks. The LTR-RTs insertion suppresses gene expression, and the suppressed intensity in gene regions was larger than in promoter regions. Some centromere and rearrangement regions indicated high-density peaks of LTR/Copia and LTR/Gypsy. Different centromeric regional repeat units (32, 78, 79 bp) were detected on different Q. glauca chromosomes.

CONCLUSION

Chromosome fusions and arm exchanges contribute to the formation of oak karyotypes. The composition and abundance of LTR-RTs are affected by its recent amplification. LTR-RTs random retrotransposition suppresses gene expression and is enriched in centromere and chromosomal rearrangement regions. This study provides novel insights into the evolutionary history of oak karyotypes and the organization, amplification, and function of LTR-RTs.

The Cissus quadrangularis genome reveals its adaptive features in an arid habitat

Cissus quadrangularis is a tetraploid species belonging to the Vitaceae family and is known for the Crassulacean acid metabolism (CAM) pathway in the succulent stem, while the leaves perform C3 photosynthesis. Here, we report a high-quality genome of C. quadrangularis comprising a total size of 679.2 Mb which was phased into two subgenomes. Genome annotation identified 51 857 protein-coding genes, while approximately 47.75% of the genome was composed of repetitive sequences. Gene expression ratios of two subgenomes demonstrated that the sub-A genome as the dominant subgenome played a vital role during the drought tolerance. Genome divergence analysis suggests that the tetraploidization event occurred around 8.9 million years ago. Transcriptome data revealed that pathways related to cutin, suberine, and wax metabolism were enriched in the stem during drought treatment, suggesting that these genes contributed to the drought adaption. Additionally, a subset of CAM-related genes displayed diurnal expression patterns in the succulent stems but not in leaves, indicating that stem-biased expression of existing genes contributed to the CAM evolution. Our findings provide insights into the mechanisms of drought adaptation and photosynthesis transition in plants.

Transcriptome profiling and characterization of genes associated with tuberization under high temperature in aeroponics in potato cv. Kufri Anand

BACKGROUND

High temperature stress is an important abiotic factor, which affects tuberization and ultimately causes heavy yield reduction in potato.

OBJECTIVES

Identification and characterization of genes associated with tuberization under high temperature stress is essential for future management through biotechnology.

METHODOLOGY

Two contrasting potato varieties Kufri Anand (profuse tuber-bearing) versus Kufri Frysona (very less/scanty tuber-bearing, control) were cultivated in aeroponics under high temperature stress, and transcriptomes were analyzed.

RESULTS

Potato cv. Kufri Anand was found superior over control (Kufri Frysona) for tuber yield and its component traits along with root morphology under aeroponics. Transcriptomes of tuber and leaf tissues were analyzed. Statistically significant (p < 0.05) differentially expressed genes (DEGs) were categorised into up-regulated (> 2 log2 fold change, FC) and down-regulated (< -2 log2 FC) genes. DEGs were annotated by gene ontology and KEGG pathways. A few selected up-regulated genes of both tissues were identified, and phylogeny tree and motif analysis were analysed based on 36 peptide sequences representing 15 selected DEGs in this study. Further, gene expression markers were developed and validated by real time qPCR analysis for the identification of high temperature tolerant genotypes.

CONCLUSION

A few key genes associated in tuberization under high temperature conditions were heat shock proteins (e.g. 18.5 kDa class I heat shock protein), sugar metabolism (e.g. glucosyltransferase), transcription factor (e.g. WRKY), and phytohormones (e.g. auxin-induced beta-glucosidase). Our study provides an overview of key genes involved in tuberization under high temperature stress in potato cv. Kufri Anand under aeroponics.

New plastids, old proteins: repeated endosymbiotic acquisitions in kareniacean dinoflagellates

Dinoflagellates are a diverse group of ecologically significant micro-eukaryotes that can serve as a model system for plastid symbiogenesis due to their susceptibility to plastid loss and replacement via serial endosymbiosis. Kareniaceae harbor fucoxanthin-pigmented plastids instead of the ancestral peridinin-pigmented ones and support them with a diverse range of nucleus-encoded plastid-targeted proteins originating from the haptophyte endosymbiont, dinoflagellate host, and/or lateral gene transfers (LGT). Here, we present predicted plastid proteomes from seven distantly related kareniaceans in three genera (Karenia, Karlodinium, and Takayama) and analyze their evolutionary patterns using automated tree building and sorting. We project a relatively limited ( ~ 10%) haptophyte signal pointing towards a shared origin in the family Chrysochromulinaceae. Our data establish significant variations in the functional distributions of these signals, emphasizing the importance of micro-evolutionary processes in shaping the chimeric proteomes. Analysis of plastid genome sequences recontextualizes these results by a striking finding the extant kareniacean plastids are in fact not all of the same origin, as two of the studied species (Karlodinium armiger, Takayama helix) possess plastids from different haptophyte orders than the rest.

Gills de novo assembly reveals oxidative stress, unfolded protein, and immune response on red cusk-eel (Genypterus chilensis) under thermal stress

The heat waves on the South Pacific coast could lead to thermal stress in native fish. The red cusk-eel (Genypterus chilensis) is relevant for Chilean artisanal fisheries and aquaculture diversification. This study examined the effect of high-temperature stress in the gills of G. chilensis in control (14 °C) and high-temperature stress (19 °C) conditions. High-temperature stress induces a significant increase in gills cortisol levels. Additionally, oxidative damage was observed in gills (protein carbonylation and lipoperoxidation). RNA-seq data was used to build the first transcriptome assembly of gills in this species (23,656 annotated transcripts). A total of 1138 down-regulated and 1531 up-regulated transcripts were observed in response to high-temperature stress in gills. The enrichment analysis showed immune response and replication enriched processes (on down-regulated transcripts), and processes related to the folding of proteins, endoplasmic reticulum, and transporter activity (on up-regulated transcripts). The present study showed how gills could be affected by high-temperature stress.

Near-gapless genome and transcriptome analyses provide insights into fruiting body development in Lentinula edodes

Fruiting body development in macrofungi is an intensive research subject. In this study, high-quality genomes were assembled for two sexually compatible monokaryons from a heterokaryotic Lentinula edodes strain WX1, and variations in L. edodes genomes were analyzed. Specifically, differential gene expression and allele-specific expression (ASE) were analyzed using the two monokaryotic genomes and transcriptome data from four different stages of fruiting body development in WX1. Results revealed that after aeration, mycelia sensed cell wall stress, pheromones, and a decrease in CO2 concentration, leading to up-regulated expression in genes related to cell adhesion, cell wall remodeling, proteolysis, and lipid metabolism, which may promote primordium differentiation. Aquaporin genes and those related to proteolysis, mitosis, lipid, and carbohydrate metabolism may play important roles in primordium development, while genes related to tissue differentiation and sexual reproduction were active in fruiting body. Several essential genes for fruiting body development were allele-specifically expressed and the two nuclear types could synergistically regulate fruiting body development by dominantly expressing genes with different functions. ASE was probably induced by long terminal repeat-retrotransposons. Findings here contribute to the further understanding of the mechanism of fruiting body development in macrofungi.

Genome architecture and genetic diversity of allopolyploid okra (Abelmoschus esculentus)

The allopolyploid okra (Abelmoschus esculentus) unveiled telomeric repeats flanking distal gene-rich regions and short interstitial TTTAGGG telomeric repeats, possibly representing hallmarks of chromosomal speciation. Ribosomal RNA (rRNA) genes organize into 5S clusters, distinct from the 18S-5.8S-28S units, indicating an S-type rRNA gene arrangement. The assembly, in line with cytogenetic and cytometry observations, identifies 65 chromosomes and a 1.45 Gb genome size estimate in a haploid sibling. The lack of aberrant meiotic configurations implies limited to no recombination among sub-genomes. k-mer distribution analysis reveals 75% has a diploid nature and 15% heterozygosity. The configurations of Benchmarking Universal Single-Copy Ortholog (BUSCO), k-mer, and repeat clustering point to the presence of at least two sub-genomes one with 30 and the other with 35 chromosomes, indicating the allopolyploid nature of the okra genome. Over 130 000 putative genes, derived from mapped IsoSeq data and transcriptome data from public okra accessions, exhibit a low genetic diversity of one single nucleotide polymorphisms per 2.1 kbp. The genes are predominantly located at the distal chromosome ends, declining toward central scaffold domains. Long terminal repeat retrotransposons prevail in central domains, consistent with the observed pericentromeric heterochromatin and distal euchromatin. Disparities in paralogous gene counts suggest potential sub-genome differentiation implying possible sub-genome dominance. Amino acid query sequences of putative genes facilitated phenol biosynthesis pathway annotation. Comparison with manually curated reference KEGG pathways from related Malvaceae species reveals the genetic basis for putative enzyme coding genes that likely enable metabolic reactions involved in the biosynthesis of dietary and therapeutic compounds in okra.

Genome assembly of a symbiotic balantidia (Balantidium ctenopharyngodoni) in fish hindgut

Balantidium ctenopharyngodoni is identified as the sole ciliate species that exclusively resides within the hindgut of grass carp with high prevalence and intensity. In this study, the successful cultivation of B. ctenopharyngodoni enabled us to collect enough cells for genome sequencing. Consequently, we acquired a high-quality genome assembly spanning 68.66 Mb, encompassing a total of 22,334 nanochromosomes. Furthermore, we predicted 29,348 protein-coding genes, and 95.5% of them was supported by the RNA-seq data. The trend of GC content in the subtelomeric regions of single-gene chromosomes was similar to other ciliates containing nanochromosomes. A large number of genes encoding carbohydrate-binding modules with affinities for starch and peptidoglycans was identified. The identification of mitochondrion-related organelles (MROs) within genome indicates its well-suited adaptation to the anaerobic conditions in the hindgut environment. In summary, our results will offer resources for understanding the genetic basis and molecular adaptations of balantidia to hindgut of herbivorous fish.

Salinity Tolerance and Osmoadaptation Strategies in Four Genera of Anammox Bacteria: Brocadia, Jettenia, Kuenenia, and Scalindua

The salinity tolerance and osmoadaptation strategies in four phylogenetically distant anammox species, Brocadia, Jettenia, Kuenenia, and Scalindua, were investigated by using highly enriched cell cultures. The first-emerged "Ca. Scalindua sp." showed optimum growth at 1.5-3% salinity and was tolerant to ∼10% salinity (a slight halophile). The second-emerged "Ca. Kuenenia stuttgartiensis" was tolerant to ∼6% salinity with optimum growth at 0.25-1.5% (a halotolerant). These early-emerged "Ca. Scalindua sp." and ″Ca. K. stuttgartiensis" rapidly accumulated K+ ions and simultaneously synthesized glutamate as a counterion. Subsequently, part of the glutamate was replaced by trehalose. In contrast, the late-emerged "Ca. B. sinica" and "Ca. J. caeni" were unable to accumulate sufficient amounts of K+─glutamate and trehalose, resulting in a significant decrease in activity even at 1-2% salinity (nonhalophiles). In addition, the external addition of glutamate may increase anammox activity at high salinity. The species-dependent salinity tolerance and osmoadaptation strategies were consistent with the genetic potential required for the biosynthesis and transport of these osmolytes and the evolutionary history of anammox bacteria: Scalindua first emerged in marine environments and then Kuenenia and other two species gradually expanded their habitat to estuaries, freshwater, and terrestrial environments, while Brocadia and Jettenia likely lost their ability to accumulate K+─glutamate.

RNA-seq analysis reveals changes in mRNA expression during development in Daphnia mitsukuri

Temporal transcriptional variation is a major contributor to functional evolution and the developmental process. Parthenogenetic water fleas of the genus Daphnia (Cladocera) provide an ideal model to characterize gene expression patterns across distinct developmental stages. Herein, we report RNA-seq data for female Daphnia mitsukuri at three developmental stages: the embryo, juvenile (three timepoints) and adult. Comparisons of gene expression patterns among these three developmental stages and weighted gene co-expression network analysis based on expression data across developmental stages identified sets of genes underpinning each of the developmental stages of D. mitsukuri. Specifically, highly expressed genes (HEGs) at the embryonic developmental stage were associated with cell proliferation, ensuring the necessary foundation for subsequent development; HEGs at the juvenile stages were associated with chemosensory perception, visual perception and neurotransmission, allowing individuals to enhance detection of potential environmental risks; HEGs at the adult stage were associated with antioxidative defensive systems, enabling adults to mount an efficient response to perceived environmental risks. Additionally, we found a significant overlap between expanded gene families of Daphnia species and HEGs at the juvenile stages, and these genes were associated with visual perception and neurotransmission. Our work provides a resource of developmental transcriptomes, and comparative analyses that characterize gene expression dynamics throughout development of Daphnia.

Comparative genomic analysis of symbiotic and free-living Fluviibacter phosphoraccumulans strains provides insights into the evolutionary origins of obligate Euplotes-bacterial endosymbioses

Endosymbiosis is a widespread and important phenomenon requiring diverse model systems. Ciliates are a widespread group of protists that often form symbioses with diverse microorganisms. Endosymbioses between the ciliate Euplotes and heritable bacterial symbionts are common in nature, and four essential symbionts were described: Polynucleobacter necessarius, "Candidatus Protistobacter heckmanni," "Ca. Devosia symbiotica," and "Ca. Devosia euplotis." Among them, only the genus Polynucleobacter comprises very close free-living and symbiotic representatives, which makes it an excellent model for investigating symbiont replacements and recent symbioses. In this article, we characterized a novel endosymbiont inhabiting the cytoplasm of Euplotes octocarinatus and found that it is a close relative of the free-living bacterium Fluviibacter phosphoraccumulans (Betaproteobacteria and Rhodocyclales). We present the complete genome sequence and annotation of the symbiotic Fluviibacter. Comparative analyses indicate that the genome of symbiotic Fluviibacter is small in size and rich in pseudogenes when compared with free-living strains, which seems to fit the prediction for recently established endosymbionts undergoing genome erosion. Further comparative analysis revealed reduced metabolic capacities in symbiotic Fluviibacter, which implies that the symbiont relies on the host Euplotes for carbon sources, organic nitrogen and sulfur, and some cofactors. We also estimated substitution rates between symbiotic and free-living Fluviibacter pairs for 233 genes; the results showed that symbiotic Fluviibacter displays higher dN/dS mean value than free-living relatives, which suggested that genetic drift is the main driving force behind molecular evolution in endosymbionts.

IMPORTANCE

In the long history of symbiosis research, most studies focused mainly on organelles or bacteria within multicellular hosts. The single-celled protists receive little attention despite harboring an immense diversity of symbiotic associations with bacteria and archaea. One subgroup of the ciliate Euplotes species is strictly dependent on essential symbionts for survival and has emerged as a valuable model for understanding symbiont replacements and recent symbioses. However, almost all of our knowledge about the evolution and functions of Euplotes symbioses comes from the Euplotes-Polynucleobacter system. In this article, we report a novel essential symbiont, which also has very close free-living relatives. Genome analysis indicated that it is a recently established endosymbiont undergoing genome erosion and relies on the Euplotes host for many essential molecules. Our results provide support for the notion that essential symbionts of the ciliate Euplotes evolve from free-living progenitors in the natural water environment.

Extremely oligotrophic and complex-carbon-degrading microaerobic bacteria from Arabian Sea oxygen minimum zone sediments

Sediments underlying marine hypoxic zones are huge sinks of unreacted complex organic matter, where despite acute O2 limitation, obligately aerobic bacteria thrive, and steady depletion of organic carbon takes place within a few meters below the seafloor. However, little knowledge exists about the sustenance and complex carbon degradation potentials of aerobic chemoorganotrophs in these sulfidic ecosystems. We isolated and characterized a number of aerobic bacterial chemoorganoheterotrophs from across a ~ 3 m sediment horizon underlying the perennial hypoxic zone of the eastern Arabian Sea. High levels of sequence correspondence between the isolates' genomes and the habitat's metagenomes and metatranscriptomes illustrated that the strains were widespread and active across the sediment cores explored. The isolates catabolized several complex organic compounds of marine and terrestrial origins in the presence of high or low, but not zero, O2. Some of them could also grow anaerobically on yeast extract or acetate by reducing nitrate and/or nitrite. Fermentation did not support growth, but enabled all the strains to maintain a fraction of their cell populations over prolonged anoxia. Under extreme oligotrophy, limited growth followed by protracted stationary phase was observed for all the isolates at low cell density, amid high or low, but not zero, O2 concentration. While population control and maintenance could be particularly useful for the strains' survival in the critically carbon-depleted layers below the explored sediment depths (core-bottom organic carbon: 0.5-1.0% w/w), metagenomic data suggested that in situ anoxia could be surmounted via potential supplies of cryptic O2 from previously reported sources such as Nitrosopumilus species.

Molecular adaptations underlying high-frequency hearing in the brain of CF bats species

BACKGROUND

The majority of bat species have developed remarkable echolocation ability, especially for the laryngeally echolocating bats along with high-frequency hearing. Adaptive evolution has been widely detected for the cochleae in the laryngeally echolocating bats, however, limited understanding for the brain which is the central to echolocation signal processing in the auditory perception system, the laryngeally echolocating bats brain may also undergo adaptive changes.

RESULT

In order to uncover the molecular adaptations related with high-frequency hearing in the brain of laryngeally echolocating bats, the genes expressed in the brain of Rhinolophus ferrumequinum (CF bat) and Myotis pilosus (FM bat) were both detected and also compared. A total of 346,891 genes were detected and the signal transduction mechanisms were annotated by the most abundant genes, followed by the transcription. In hence, there were 3,088 DEGs were found between the two bat brains, with 1,426 highly expressed in the brain of R. ferrumequinum, which were significantly enriched in the neuron and neurodevelopmental processes. Moreover, we found a key candidate hearing gene, ADCY1, playing an important role in the R. ferrumequinum brain and undergoing adaptive evolution in CF bats.

CONCLUSIONS

Our study provides a new insight to the molecular bases of high-frequency hearing in two laryngeally echolocating bats brain and revealed different nervous system activities during auditory perception in the brain of CF bats.

Comparative genomics analysis reveals sequence characteristics potentially related to host preference in Cryptosporidium xiaoi

Cryptosporidium spp. are important diarrhea-associated pathogens in humans and livestock. Among the known species, Cryptosporidium xiaoi, which causes cryptosporidiosis in sheep and goats, was previously recognized as a genotype of the bovine-specific Cryptosporidium bovis based on their high sequence identity in the ssrRNA gene. However, the lack of genomic data has limited characterization of the genetic differences between the two closely related species. In this study, we sequenced the genomes of two C. xiaoi isolates and performed comparative genomic analysis to identify the sequence uniqueness of this ovine-adapted species compared with other Cryptosporidium spp. Our results showed that C. xiaoi is genetically related to C. bovis as shown by their 95.8% genomic identity and similar gene content. Consistent with this, both C. xiaoi and C. bovis appear to have fewer genes encoding mitochondrial metabolic enzymes and invasion-related protein families. However, they appear to possess several species-specific genes. Further analysis indicates that the sequence differences between these two Cryptosporidium spp. are mainly in 24 highly polymorphic genes, half of which are located in the subtelomeric regions. Some of these subtelomeric genes encode secretory proteins that have undergone positive selection. In addition, the genomes of two C. xiaoi isolates, identified as subtypes XXIIIf and XXIIIh, share 99.9% nucleotide sequence identity, with six highly divergent genes encoding putative secretory proteins. Therefore, these species-specific genes and sequence polymorphism in subtelomeric genes probably contribute to the different host preference of C. xiaoi and C. bovis.