Using the Basic Local Alignment Search Tool (BLAST)

Chromosome-level genome assembly of tetraploid Chinese cherry (Prunus pseudocerasus)

Chinese cherry belongs to the family Rosaceae, genus Prunus, and has high nutritional and economic value. 'Duiying' is a Chinese cherry variety local to Beijing, and has better performance than sweet cherry in terms of disease resistance. However, disease resistance resources of 'Duiying' have not been fully exploited partially due to the lack of a high-quality genome. In this study, we report a high-quality chromosome-scale genome assembly for Chinese cherry 'Duiying', by combining PacBio HiFi, Bionano and Hi-C sequencing data. The assembled genome has a size of 1035.19 Mb, with a scaffold N50 of 28.99 Mb, and 978.61 Mb (94.54%) assembled into 32 pseudochromosomes. A total of 547.16 Mb (52.86%) sequences were identified as repetitive sequences, and 114,451 protein-coding genes were annotated. Moreover, a total of 1635 microRNA (miRNA), 6637 transfer RNA (tRNA), 38,258 ribosomal RNA (rRNA), and 169 small nuclear RNAs (snRNA) genes were identified. The genome assembly presented here provides valuable genomic resources to enhance our understanding of genetic and molecular basis of Chinese cherry.

A haplotype-resolved and chromosome-scale genome assembly of Hainan muntjac (Muntiacus nigripes)

The Hainan muntjac (Muntiacus nigripes) is a wild animal endemic to Hainan, China. Its species distribution and the diversity of muntjac karyotypes have attracted much attention. Although genomic resources have increased in recent years, relevant genome assembly data of Hainan muntjac are still lacking. Meanwhile, molecular evidence for the taxonomic units of this species remains lacking. In this study, we successfully assembled the Hainan muntjac haplotype genome at the chromosome level using Pacbio long read and long sequencing technologies and Hi-C data. The final assembly size was 2.66 Gb, with allele and scaffold N50 values of 29.27 and 700.27 Mb, respectively, and we scaffolded the genome sequence onto three chromosomes. The genome contains a total of 21,451 genes and 10,056 gene families. Phylogenetic analysis using single-copy gene families revealed that Hainan muntjac is most closely related to red muntjac, with a divergence time of 8-11 Ma. This new genomic resource of Hainan muntjac will be crucial for future comparative genomic analyses and genetic evolutionary studies.

Chromosome-scale genome assembly and gene annotation of the Alligator Gar (Atractosteus spatula)

Given the aggressive nature and robust survival capabilities of the alligator gar (Atractosteus spatula), if it was to exist in a new environment as an invasive species, it could cause significant disruption to the invaded ecosystem. Building on the continuity and completeness of the existing draft genome were not optimal, this study has updated a high-quality genome of the alligator gar at the chromosome level, which was assembled using Oxford Nanopore Technology and chromatin interaction mapping (Hi-C) sequencing techniques. In summary, the alligator gar genome in this study was 1.05 Gb in size with a contig N50 of 15.7 Mb and scaffold N50 of 56.8 Mb. We captured 98.26% of assembled bases in 28 pseudochromosomes. The completeness of the final chromosome-level genome reached 96.7%. Meanwhile, a total of 19,103 protein-coding genes were predicted, of which 99.83% could be predicted with functions. Taken together, the present high-quality alligator gar chromosome-level genome provides a valuable resource for exploring the underlying genomic basis to comprehend the functional genomics, chromosome evolution, and population management of this species.

The whole chromosome-level genome provides resources and insights into the endangered fish Percocypris pingi evolution and conservation

BACKGROUND

Percocypris pingi (Tchang) was classified as Endangered on the Red List of China's Vertebrates in 2015 and is widely distributed in the Upper Yangtze River. Although breeding and release into wild habitats have been performed for this commercially important fish in recent years, low genetic diversity has been found in wild populations. Genomic resources are strongly recommended before formulating and carrying out conservation strategies for P. pingi. Thus, there is an urgent need to conserve germplasm resources and improve the population diversity of P. pingi. To date, the whole genome of P. pingi has not been reported.

RESULTS

In our study, we constructed the first chromosome-level genome of P. pingi by high-throughput chromosome conformation capture (Hi-C) technology and PacBio long-read sequencing. The assembled genome was 1.7 Gb in size, with an N50 of 17,692 bp and a GC content from circular consensus sequencing of 37.67%. The Hi-C results again demonstrated that P. pingi was tetraploid (n = 98), with the genome consisting of 24-type and 25-type chromosomes. Chr.19 of the 24-type chromosomes in P. pingi resulted from the fusion of chr.19 and chr.22 in zebrafish. The divergence times between 24-type and 25-type chromosomes was around 6.1 million years ago. A total of 25,198 and 25,291 protein-coding genes were obtained from the 24-type and 25-type chromosomes, respectively. The ploidy of P. pingi is an allotetraploid. A total of 8,741 genes of P. pingi were clustered into 4,378 gene families that were shared with 14 other species, and the P. pingi genome had 68 unique gene families. Phylogenetic analyses indicated that P. pingi was most closely related to Schizothorax oconnori, and the genes were clustered on one branch. We identified 166 significantly expanded gene families and 173 significantly contracted gene families in P. pingi. The most enriched positive protein-coding genes, such as Bmp-4, Etfdh, homeobox protein HB9, and ATG3, were screened.

CONCLUSION

Our study provides a high-quality chromosome-anchored reference genome for P. pingi and provides sufficient information on the chromosomes, which will lead to valuable resources for genetic, genomic, and biological studies of P. pingi and for improving the genetic diversity, population size, and scientific conservation of endangered fish and other key cyprinid species in aquaculture.

Identification of a novel monocyclic carotenoid and prediction of its biosynthetic genes in Algoriphagus sp. oki45

Bacteria belonging to the genus Algoriphagus have been isolated from various sources, such as Antarctic sea ice, seawater, and sediment, and some strains are known to produce orange to red pigments. However, the pigment composition and biosynthetic genes have not been fully elucidated. A new red-pigmented Algoriphagus sp. strain, oki45, was isolated from the surface of seaweed collected from Senaga-Jima Island, Okinawa, Japan. Genome comparison revealed oki45's average nucleotide identity of less than 95% to its closely related species, Algoriphagus confluentis NBRC 111222 T and Algoriphagus taiwanensis JCM 19755 T. Comprehensive chemical analyses of oki45's pigments, including 1H and 13C nuclear magnetic resonance and circular dichroism spectroscopy, revealed that the pigments were mixtures of monocyclic carotenoids, (3S)-flexixanthin ((3S)-3,1'-dihydroxy-3',4'-didehydro-1',2'-dihydro-β,ψ-caroten-4-one) and (2R,3S)-2-hydroxyflexixanthin ((2R,3S)-2,3,1'-trihydroxy-3',4'-didehydro-1',2'-dihydro-β,ψ-caroten-4-one); in particular, the latter compound was new and not previously reported. Both monocyclic carotenoids were also found in A. confluentis NBRC 111222 T and A. taiwanensis JCM 19755 T. Further genome comparisons of carotenoid biosynthetic genes revealed the presence of eight genes (crtE, crtB, crtI, cruF, crtD, crtYcd, crtW, and crtZ) for flexixanthin biosynthesis. In addition, a crtG homolog gene encoding 2,2'-β-hydroxylase was found in the genome of the strains oki45, A. confluentis NBRC 111222 T, and A. taiwanensis JCM 19755 T, suggesting that the gene is involved in 2-hydroxyflexixanthin synthesis via 2-hydroxylation of flexixanthin. These findings expand our knowledge of monocyclic carotenoid biosynthesis in Algoriphagus bacteria. KEY POINTS: • Algoriphagus sp. strain oki45 was isolated from seaweed collected in Okinawa, Japan. • A novel monocyclic carotenoid 2-hydroxyflexixanthin was identified from strain oki45. • Nine genes for 2-hydroxyflexixanthin biosynthesis were found in strain oki45 genome.

Salmonella Phage vB_SpuM_X5: A Novel Approach to Reducing Salmonella Biofilms with Implications for Food Safety

Salmonella, a prevalent foodborne pathogen, poses a significant social and economic strain on both food safety and public health. The application of phages in the control of foodborne pathogens represents an emerging research area. In this study, Salmonella pullorum phage vB_SpuM_X5 (phage X5) was isolated from chicken farm sewage samples. The results revealed that phage X5 is a novel Myoviridae phage. Phage X5 has adequate temperature tolerance (28 °C-60 °C), pH stability (4-12), and a broad host range of Salmonella bacteria (87.50% of tested strains). The addition of phage X5 (MOI of 100 and 1000) to milk inoculated with Salmonella reduced the number of Salmonella by 0.72 to 0.93 log10 CFU/mL and 0.66 to 1.06 log10 CFU/mL at 4 °C and 25 °C, respectively. The addition of phage X5 (MOI of 100 and 1000) to chicken breast inoculated with Salmonella reduced bacterial numbers by 1.13 to 2.42 log10 CFU/mL and 0.81 to 1.25 log10 CFU/mL at 4 °C and 25 °C, respectively. Phage X5 has bactericidal activity against Salmonella and can be used as a potential biological bacteriostatic agent to remove mature biofilms of Salmonella or for the prevention and control of Salmonella.

Chromosome-level genome assembly of the smallscale yellowfin (Plagiognathops microlepis)

The small-scale yellowfin (Plagiognathops microlepis) is a highly valued species in East Asian aquaculture due to its adaptability and high yield. However, the lack of genomic data has impeded genetic research and breeding efforts. In this study, we utilize PacBio Hifi long-read sequencing and Hi-C technologies to construct a highly detailed genome of P. microlepis at the chromosomal level. The assembly encompasses 976.41 Mb, with an exceptional 99.84% distribution across 24 chromosomes. Notably, the contig N50 was 34.41 Mb and scaffold N50 was 38.38 Mb. The completeness of the P. microlepis genome assembly is underscored by a BUSCO score of 98.08%. A total of 25,389 protein-coding genes were identified, with a BUSCO score of 96.98%, and 99.85% of these genes were functionally annotated. Synteny relationships at the chromosome level with Danio rerio and Chanodichthys erythropterus genomes uncover small-scale chromosomal rearrangements. This high-fidelity genome assembly serves as a pivotal resource for forthcoming endeavors such as the genome structure, functional elements, comparative genomics, and evolutionary characteristics of P. microlepis and its relative species.

Chromosome-level genome assembly and annotation of the crested gecko, Correlophus ciliatus, a lizard incapable of tail regeneration

The vast majority of gecko species are capable of tail regeneration, but singular geckos of Correlophus, Uroplatus, and Nephrurus genera are unable to regrow lost tails. Of these non-regenerative geckos, the crested gecko (Correlophus ciliatus) is distinguished by ready availability, ease of care, high productivity, and hybridization potential. These features make C. ciliatus particularly suited as a model for studying the genetic, molecular, and cellular mechanisms underlying loss of tail regeneration capabilities. We report a contiguous genome of C. ciliatus with a total size of 1.65 Gb, 152 scaffolds, L50 of 6, and N50 of 109 Mb. Repetitive content consists of 40.41% of the genome, and a total of 30,780 genes were annotated. Our assembly of the crested gecko genome provides a valuable resource for future comparative genomic studies between non-regenerative and regenerative geckos and other squamate reptiles.

Findings

We report genome sequencing, assembly, and annotation for the crested gecko, Correlophus ciliatus.

High-quality genome of a modern soybean cultivar and resequencing of 547 accessions provide insights into the role of structural variation

Soybean provides protein, oil and multiple health-related compounds. Understanding the effects of structural variations (SVs) on economic traits in modern breeding is important for soybean improvement. Here we assembled the high-quality genome of modern cultivar Nongdadou2 (NDD2) and identified 25,814 SV-gene pairs compared to 29 reported genomes, with 13 NDD2-private SVs validated in 547 deep-resequencing (average = 18.05-fold) accessions, which advances our understanding of genomic variation biology. We found some insertions/deletions involved in seed protein and weight formation, an inversion related to adaptation to drought and a large intertranslocation implicated in a key divergence event in soybean. Of 749,714 SVs from 547 accessions, 6,013 were significantly associated with 22 yield-related and seed-quality-related traits determined in ten location × year environments. We uncovered 1,761 associated SVs that hit genes or regulatory regions, with 12 in GmMQT influencing oil and isoflavone contents. Our work provides resources and insights into SV roles in soybean improvement.

Chromosome-level genome assembly of the forest pest Achelura yunnanensis (Lepidoptera: Zygaenidae)

Achelura yunnanensis is a destructive pest of forests, causing substantial damage on tree growth and severe economic losses. Additionally, as a daytime-active moth, this species also holds important scientific value for investigating the genetic mechanisms governing day-night activity patterns of Lepidoptera. To facilitate effective pest control and deepen our understanding of the diurnal behavior's genetic basis of moths, genomic data for this species are crucial. In this study, we present a chromosome-level reference genome of A. yunnanensis (368.15 Mb in 32 chromosomes; scaffold N50 = 12.61 Mb; BUSCO completeness = 98.0%). Genome annotation shows that the new assembly comprises 37.10% (136.55 Mb) repetitive elements, 1,828 non-coding RNAs, and 15,523 protein-coding genes. Genes involved in lipid metabolism and xenobiotics biodegradation and metabolism, such as cytochrome P450 families, experienced significant expansion in the A. yunnanensis genome. The chromosome-level genome of A. yunnanensis provides a valuable genomic resource for devising novel pest control strategies, and will also help to study the genetic mechanism of the shift of diurnal behavior in Lepidoptera.

Metagenome-Assembled Genome of "Candidatus Aramenus sp. CH1" from the Chichon volcano, Mexico

The Chichon volcano contains several thermal manifestations including an acidic crater lake. Here we report a metagenome-assembled genome of "Candidatus Aramenus sp. CH1," a Sulfolobales archaeon inhabiting the crater lake from the Chichon volcano. In this study, we generated a novel Aramenus genome sequence from a thermal area in Southern Mexico.

Biostimulant and antagonistic potential of endophytic fungi against fusarium wilt pathogen of tomato Fusarium oxysporum f. sp. lycopersici

Endophytic fungal-based biopesticides are sustainable and ecologically-friendly biocontrol agents of several pests and diseases. However, their potential in managing tomato fusarium wilt disease (FWD) remains unexploited. This study therefore evaluated effectiveness of nine fungal isolates against tomato fusarium wilt pathogen, Fusarium oxysporum f. sp. lycopersici (FOL) in vitro using dual culture and co-culture assays. The efficacy of three potent endophytes that inhibited the pathogen in vitro was assessed against FWD incidence, severity, and ability to enhance growth and yield of tomatoes in planta. The ability of endophytically-colonized tomato (Solanum lycopersicum L.) plants to systemically defend themselves upon exposure to FOL were also assessed through defence genes expression using qPCR. In vitro assays showed that endophytes inhibited and suppressed FOL mycelial growth better than entomopathogenic fungi (EPF). Endophytes Trichoderma asperellum M2RT4, Hypocrea lixii F3ST1, Trichoderma harzianum KF2R41, and Trichoderma atroviride ICIPE 710 had the highest (68.84-99.61%) suppression and FOL radial growth inhibition rates compared to EPF which exhibited lowest (27.05-40.63%) inhibition rates. Endophytes T. asperellum M2RT4, H. lixii F3ST1 and T. harzianum KF2R41 colonized all tomato plant parts. During the in planta experiment, endophytically-colonized and FOL-infected tomato plants showed significant reduction of FWD incidence and severity compared to non-inoculated plants. In addition, these endophytes contributed to improved growth promotion parameters and yield. Moreover, there was significantly higher expression of tomato defence genes in T. asperellum M2RT4 colonized than in un-inoculated tomato plants. These findings demonstrated that H. lixii F3ST1 and T. asperellum M2RT4 are effective biocontrol agents against FWD and could sustainably mitigate tomato yield losses associated with fusarium wilt.

Functional exploration of copy number alterations in a Drosophila model of triple-negative breast cancer

Accounting for 10-20% of breast cancer cases, triple-negative breast cancer (TNBC) is associated with a disproportionate number of breast cancer deaths. One challenge in studying TNBC is its genomic profile: with the exception of TP53 loss, most breast cancer tumors are characterized by a high number of copy number alterations (CNAs), making modeling the disease in whole animals challenging. We computationally analyzed 186 CNA regions previously identified in breast cancer tumors to rank genes within each region by likelihood of acting as a tumor driver. We then used a Drosophila p53-Myc TNBC model to identify 48 genes as functional drivers. To demonstrate the utility of this functional database, we established six 3-hit models; altering candidate genes led to increased aspects of transformation as well as resistance to the chemotherapeutic drug fluorouracil. Our work provides a functional database of CNA-associated TNBC drivers, and a template for an integrated computational/whole-animal approach to identify functional drivers of transformation and drug resistance within CNAs in other tumor types.

Pharmacophore modelling based virtual screening and molecular dynamics identified the novel inhibitors and drug targets against Waddlia chondrophila

Waddlia chondrophila is a possible cause of fetal death in humans. This Chlamydia-related bacterium is an emergent pathogen that causes human miscarriages and ruminant abortions, which results in financial losses. Despite the years of efforts, the underlying mechanism behind the pathogenesis of W. chondrophila is little known which hindered the development of novel treatment options. In the framework of current study, computational approaches were used to identify novel inhibitors (phytocompounds) and drug targets against W. chondrophila. At first, RNA polymerase sigma factor SigA and 3-deoxy-D-manno-octulosonic acid transferase were identified through subtractive proteomics pipeline. Afterwards, extensive docking and simulation analyses were conducted to optimize potentially novel phytocompounds by assessing their binding affinity to target proteins. A 100ns molecular dynamics simulation well complimented the compound's binding affinity and indicated strong stability of predicted compounds at the docked site. The calculation of binding free energies with MMGBSA corroborated the significant binding affinity between phytocompounds and target protein binding sites. The proposed phytocompounds may be a viable treatment option for patients infected with W. chondrophila; however, further research is required to ensure their safety.

Chromosome-level genome assembly of the freshwater mussel Sinosolenaia oleivora (Heude, 1877)

Sinosolenaia oleivora (Bivalve, Unionida, Unionidae), is a near-endangered edible mussel. In 2022, it was selected by the Ministry of Agriculture and Rural Affairs as a top-ten aquatic germplasm resource, with potential for industrial development. Using Illumina, PacBio, and Hi-C technology, a high-quality chromosome-level genome of S. oleivora was assembled. The assembled S. oleivora genome spanned 2052.29 Mb with a contig N50 size of 20.36 Mb and a scaffold N50 size of 103.57 Mb. The 302 contigs, accounting for 98.41% of the total assembled genome, were anchored into 19 chromosomes using Hi-C scaffolding. A total of 1171.78 Mb repeat sequences were annotated and 22,971 protein-coding genes were predicted. Compared with the nearest ancestor, a total of 603 expanded and 1767 contracted gene families were found. This study provides important genomic resources for conservation, evolutionary research, and genetic improvements of many economic traits like growth performance.

Gap-free chromosome-level genomes of male and female spotted longbarbel catfish Hemibagrus guttatus

Hemibagrus guttatus, also named as spotted longbarbel catfish, is an economical fish in China. However, their gender cannot be easily distinguished from their appearance, which largely impedes their artificial breeding. Therefore, we provided two gap-free chromosome-level genomes of male and female spotted longbarbel catfish by combining wtdbg2, LR_Gapcloser and TGS-GapCloser assembly approaches with Hi-C data and accurate Pacbio HiFi long-reads. We assembled 30 chromosomes without any gap. Their genome sizes are approximately 749.1 Mb and 747.8 Mb of male and female individuals. The completeness results of BUSCO evaluation show about 94.2% and 95.0%, representing a high-level of completeness of both genomes. We also obtained 35,277 and 34,571 protein-coding gene sets from male and female individuals. Both available gap-free chromosome-level genomes of H. guttatus will provide excellent references for resequencing of male and female individuals to identify accurate markers for distinguishing gender of this fish.

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.

Phylogenetic relationships, selective pressure and molecular markers development of six species in subfamily Polygonoideae based on complete chloroplast genomes

The subfamily Polygonoideae encompasses a diverse array of medicinal and horticultural plants that hold significant economic value. However, due to the lack of a robust taxonomy based on phylogenetic relationships, the classification within this family is perplexing, and there is also a scarcity of reports on the chloroplast genomes of many plants falling under this classification. In this study, we conducted a comprehensive analysis by sequencing and characterizing the complete chloroplast genomes of six Polygonoideae plants, namely Pteroxygonum denticulatum, Pleuropterus multiflorus, Pleuropterus ciliinervis, Fallopia aubertii, Fallopia dentatoalata, and Fallopia convolvulus. Our findings revealed that these six plants possess chloroplast genomes with a typical quadripartite structure, averaging 162,931 bp in length. Comparative chloroplast analysis, codon usage analysis, and repetitive sequence analysis demonstrated a high level of conservation within the chloroplast genomes of these plants. Furthermore, phylogenetic analysis unveiled a distinct clade occupied by P. denticulatum, while P. ciliinrvis displayed a closer relationship to the three plants belonging to the Fallopia genus. Selective pressure analysis based on maximum likelihood trees showed that a total of 14 protein-coding genes exhibited positive selection, with psbB and ycf1 having the highest number of positive amino acid sites. Additionally, we identified four molecular markers, namely petN-psbM, psal-ycf4, ycf3-trnS-GGA, and trnL-UAG-ccsA, which exhibit high variability and can be utilized for the identification of these six plants.

Molecular characterization of the complete genome of a novel ormycovirus infecting the ectomycorrhizal fungus Hortiboletus rubellus

Advancements in high-throughput sequencing and the development of new bioinformatics tools for large-scale data analysis play a crucial role in uncovering virus diversity and enhancing our understanding of virus evolution. The discovery of the ormycovirus clades, a group of RNA viruses that are phylogenetically distinct from all known Riboviria members and are found in fungi, highlights the value of these tools for the discovery of novel viruses. The aim of this study was to examine viral populations in fungal hosts to gain insights into the diversity, evolution, and classification of these viruses. Here, we report the molecular characterization of a newly discovered ormycovirus, which we have named "Hortiboletus rubellus ormycovirus 1" (HrOMV1), that was found in the ectomycorrhizal fungus Hortiboletus rubellus. The bipartite genome of HrOMV1, whose nucleotide sequence was determined by HTS and RLM-RACE, consists of two RNA segments (RNA1 and RNA2) that exhibit similarity to those of previously studied ormycoviruses in their organization and the proteins they encode. The presence of upstream, in-frame AUG triplets in the 5' termini of both RNA segments suggests that HrOMV1, like certain other ormycoviruses, employs a non-canonical translation initiation strategy. Phylogenetic analysis showed that HrOMV1 is positioned within the gammaormycovirus clade. Its putative RNA-dependent RNA polymerase (RdRp) exhibits sequence similarity to those of other gammaormycovirus members, the most similarity to that of Termitomyces ormycovirus 1, with 33.05% sequence identity. This protein was found to contain conserved motifs that are crucial for RNA replication, including the distinctive GDQ catalytic triad observed in gammaormycovirus RdRps. The results of this study underscore the significance of investigating the ecological role of mycoviruses in mycorrhizal fungi. This is the first report of an ormycovirus infecting a member of the ectomycorrhizal genus Hortiboletus.

Biocontrol of Alternaria alternata in cold-stored table grapes using psychrotrophic yeasts and bioactive compounds of natural sources

Alternaria alternata is a common fungal pathogen causing postharvest decay in table grapes. This study addressed the potential of autochthonous yeasts and bioactive compounds of natural sources to act as biocontrol agents (BCAs) against A. alternata in cold-stored table grapes. With this purpose, 19 yeast capable of growing at 0-1 °C were isolated from the surface of Red Globe table grapes. These isolates, along with the pre-isolated strain Metschnikowia pulcherrima RCM2, were evaluated as BCAs in wounded berries. From these results, six yeast isolates were pre-selected to be combined with bioactive compounds of natural sources, like phenolic compounds (PCs) of side streams of wine industry, including bunch stem extract (BSE) (5-25 %), and cane extract (CE) (5-25 %), and functional polysaccharides from shrimp waste such as chitosan (CH) (0.5 %). Then, the biocontrol efficacy of combined treatments beyond individual ones was compared. The results revealed that 4 yeast isolates, namely M. pulcherrima RCM2 and ULA146, and Aureobasidium pullulans FUL14 and FUL18, were the most effective. However, when combined with the natural bioactive compounds, their efficacy against A. alternata did not increase significantly. Notably, ULA146 and FUL18 demonstrated a biocontrol efficacy of 36-37 %, comparable to that of the treatment with commercial doses of SO2, which only showed a 27 % reduction in the lesion diameter. These findings highlight the potential of using psychrotrophic yeasts as BCAs against A. alternata in cold-stored table grapes. Combining these yeast strains with BSE, CE and CH did not increase BCAs efficacy against this pathogen at the concentrations tested. The development of effective biocontrol strategies for A. alternata could contribute to reducing reliance on chemically synthesized fungicides, promoting sustainable practices, aiming to improve the quality and safety of cold-stored table grapes.

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.

CLAUDIO: automated structural analysis of cross-linking data

MOTIVATION

Cross-linking mass spectrometry has made remarkable advancements in the high-throughput characterization of protein structures and interactions. The resulting pairs of cross-linked peptides typically require geometric assessment and validation, given the availability of their corresponding structures.

RESULTS

CLAUDIO (Cross-linking Analysis Using Distances and Overlaps) is an open-source software tool designed for the automated analysis and validation of different varieties of large-scale cross-linking experiments. Many of the otherwise manual processes for structural validation (i.e. structure retrieval and mapping) are performed fully automatically to simplify and accelerate the data interpretation process. In addition, CLAUDIO has the ability to remap intra-protein links as inter-protein links and discover evidence for homo-multimers.

AVAILABILITY AND IMPLEMENTATION

CLAUDIO is available as open-source software under the MIT license at https://github.com/KohlbacherLab/CLAUDIO.

Discovery of novel inhibitor via molecular dynamics simulations against D-alanyl-D-alanine carboxypeptidase of Enterobacter cloacae

Antibiotics resistance by bacterial pathogens is a major concern to public health worldwide resulting in high health care costs and rising mortality. Subtractive proteomics prioritized D-alanyl-D-alanine carboxypeptidas (DacB) enzyme from Enterobacter cloacae ATCC 13047 as a potential candidate for drugs designing to block pathogen cell wall biosynthesis. Virtual screening of an antibacterial library against the target unraveled a hit compound (2-[(1-methylsulfonylpiperidin-3-yl)methyl]-6-(1H-pyrazol-4-yl) pyrazine) showing high affinity and stability with the target. The N-methyl-N-propyl-methanesulfonamide of the compound is seen as a closed affinity towards domain involving strong hydrogen bonds with Ser41, Lys44, Ser285, and Asn287. The 4-methyl-1H-pyrazole is posed towards the open cavity of domain I and II and formed hydrophobic and hydrophilic contacts. The system is highly stable with average carbon-alpha deviations of 1.69 Å over trajectories of 400-ns. Three vital residues projected: Arg437, Arg438 and Leu400 from enzyme pocket via Radial distribution function (RDF) assay, which actively engaged the inhibitor. Further confirmation is done by estimating binding free energies, which confirms the very low delta energy of -7.24 kcal/mol in Generalized Born (GB) method and -7.4363 kcal/mol in Poisson-Boltzmann (PB) method. WaterSwap calculations were performed that revealed the energies highly converged, an agreement on good system stability. Lastly, three DacB mutants were created to investigate the role of functional active residues and a decline in binding affinity of the residues was noticed. These computational results provide a gateway for experimentalists to further confirm their efficacy both in-vitro and in-vivo.Communicated by Ramaswamy H. Sarma.

The genome assembly and annotation of the white-lipped tree pit viper Trimeresurus albolabris

Trimeresurus albolabris, also known as the white-lipped pit viper or white-lipped tree viper, is a highly venomous snake distributed across Southeast Asia and the cause of many snakebite cases. In this study, we report the first whole genome assembly of T. albolabris obtained with next-generation sequencing from a specimen collected in Mengzi, Yunnan, China. After genome sequencing and assembly, the genome of this male T. albolabris individual was 1.51 Gb in length and included 38.42% repeat-element content. Using this genome, 21,695 genes were identified, and 99.17% of genes could be annotated using gene functional databases. Our genome assembly and annotation process was validated using a phylogenetic tree, which included six species and focused on single-copy genes of nuclear genomes. This research will contribute to future studies on Trimeresurus biology and the genetic basis of snake venom.

Haplotype-resolved chromosome-scale genomes of the Asian and African Savannah Elephants

The Proboscidea, which includes modern elephants, were once the largest terrestrial animals among extant species. They suffered mass extinction during the Ice Age. As a unique branch on the evolutionary tree, the Proboscidea are of great significance for the study of living animals. In this study, we generate chromosome-scale and haplotype-resolved genome assemblies for two extant Proboscidea species (Asian Elephant, Elephas maximus and African Savannah Elephant, Loxodonta africana) using Pacbio, Hi-C, and DNBSEQ technologies. The assembled genome sizes of the Asian and African Savannah Elephant are 3.38 Gb and 3.31 Gb, with scaffold N50 values of 130 Mb and 122 Mb, respectively. Using Hi-C technology ~97% of the scaffolds are anchored to 29 pseudochromosomes. Additionally, we identify ~9 Mb Y-linked sequences for each species. The high-quality genome assemblies in this study provide a valuable resource for future research on ecology, evolution, biology and conservation of Proboscidea species.

Genomic exploration of the endangered oriental stork, Ciconia boyciana, sheds light on migration adaptation and future conservation

BACKGROUND

The oriental stork, Ciconia boyciana, is an endangered migratory bird listed on the International Union for Conservation of Nature's Red List. The bird population has experienced a rapid decline in the past decades, with nest locations and stop-over sites largely degraded due to human-bird conflicts. Multipronged conservation efforts are required to secure the future of oriental storks. We propose that a thorough understanding of the genome-wide genetic background of this threatened bird species is critical to make future conservation strategies.

FINDINGS

In this study, the first chromosome-scale reference genome was presented for the oriental stork with high quality, contiguity, and accuracy. The assembled genome size was 1.24 Gb with a scaffold N50 of 103 Mb, and 1.23 Gb contigs (99.32%) were anchored to 35 chromosomes. Population genomic analysis did not show a genetic structure in the wild population. Genome-wide genetic diversity (π = 0.0012) of the oriental stork was at a moderate to high level among threatened bird species, and the inbreeding risk was also not significant (FROH = 5.56% ± 5.30%). Reconstruction of demographic history indicated a rapid recent population decline likely driven by human activities. Genes that were under positive selection associated with the migratory trait were identified in relation to the long-term potentiation, photoreceptor cell organization, circadian rhythm, muscle development, and energy metabolism, indicating the essential interplay between genetic and ecological adaptation.

CONCLUSIONS

Our study presents the first chromosome-scale genome assembly of the oriental stork and provides a genomic basis for understanding a genetic background of the oriental stork, the population's extinction risks, and the migratory characteristics, which will facilitate the decision of future conservation plans for this species.

Functionality of the Na+-translocating NADH:quinone oxidoreductase and quinol:fumarate reductase from Prevotella bryantii inferred from homology modeling

Members of the family Prevotellaceae are Gram-negative, obligate anaerobic bacteria found in animal and human microbiota. In Prevotella bryantii, the Na+-translocating NADH:quinone oxidoreductase (NQR) and quinol:fumarate reductase (QFR) interact using menaquinone as electron carrier, catalyzing NADH:fumarate oxidoreduction. P. bryantii NQR establishes a sodium-motive force, whereas P. bryantii QFR does not contribute to membrane energization. To elucidate the possible mode of function, we present 3D structural models of NQR and QFR from P. bryantii to predict cofactor-binding sites, electron transfer routes and interaction with substrates. Molecular docking reveals the proposed mode of menaquinone binding to the quinone site of subunit NqrB of P. bryantii NQR. A comparison of the 3D model of P. bryantii QFR with experimentally determined structures suggests alternative pathways for transmembrane proton transport in this type of QFR. Our findings are relevant for NADH-dependent succinate formation in anaerobic bacteria which operate both NQR and QFR.

Identification of an anti-CRISPR protein that inhibits the CRISPR-Cas type I-B system in Clostridioides difficile

IMPORTANCE

Clostridioides difficile is the widespread anaerobic spore-forming bacterium that is a major cause of potentially lethal nosocomial infections associated with antibiotic therapy worldwide. Due to the increase in severe forms associated with a strong inflammatory response and higher recurrence rates, a current imperative is to develop synergistic and alternative treatments for C. difficile infections. In particular, phage therapy is regarded as a potential substitute for existing antimicrobial treatments. However, it faces challenges because C. difficile has highly active CRISPR-Cas immunity, which may be a specific adaptation to phage-rich and highly crowded gut environment. To overcome this defense, C. difficile phages must employ anti-CRISPR mechanisms. Here, we present the first anti-CRISPR protein that inhibits the CRISPR-Cas defense system in this pathogen. Our work offers insights into the interactions between C. difficile and its phages, paving the way for future CRISPR-based applications and development of effective phage therapy strategies combined with the engineering of virulent C. difficile infecting phages.

Genome assembly and annotation of the king ratsnake, Elaphe carinata

The king ratsnake (Elaphe carinata) of the genus Elaphe is a common large, non-venomous snake widely distributed in Southeast and East Asia. It is an economically important farmed species. As a non-venomous snake, the king ratsnake predates venomous snakes, such as cobras and pit vipers. However, the immune and digestive mechanisms of the king ratsnake remain unclear. Despite their economic and research importance, we lack genomic resources that would benefit toxicology, phylogeography, and immunogenetics studies. Here, we used single-tube long fragment read sequencing to generate the first contiguous genome of a king ratsnake from Huangshan City, Anhui province, China. The genome size is 1.56 GB with a scaffold N50 of 6.53M. The total length of the genome is approximately 621 Mb, while the repeat content is 42.26%. Additionally, we predicted 22,339 protein-coding genes, including 22,065 with functional annotations. Our genome is a potentially useful addition to those available for snakes.

Maintenance of UK bread baking quality: Trends in wheat quality traits over 50 years of breeding and potential for future application of genomic-assisted selection

Improved selection of wheat varieties with high end-use quality contributes to sustainable food systems by ensuring productive crops are suitable for human consumption end-uses. Here, we investigated the genetic control and genomic prediction of milling and baking quality traits in a panel of 379 historic and elite, high-quality UK bread wheat (Triticum eastivum L.) varieties and breeding lines. Analysis of the panel showed that genetic diversity has not declined over recent decades of selective breeding while phenotypic analysis found a clear trend of increased loaf baking quality of modern milling wheats despite declining grain protein content. Genome-wide association analysis identified 24 quantitative trait loci (QTL) across all quality traits, many of which had pleiotropic effects. Changes in the frequency of positive alleles of QTL over recent decades reflected trends in trait variation and reveal where progress has historically been made for improved baking quality traits. It also demonstrates opportunities for marker-assisted selection for traits such as Hagberg falling number and specific weight that do not appear to have been improved by recent decades of phenotypic selection. We demonstrate that applying genomic prediction in a commercial wheat breeding program for expensive late-stage loaf baking quality traits outperforms phenotypic selection based on early-stage predictive quality traits. Finally, trait-assisted genomic prediction combining both phenotypic and genomic selection enabled slightly higher prediction accuracy, but genomic prediction alone was the most cost-effective selection strategy considering genotyping and phenotyping costs per sample.

Genome assembly and annotation of the Brown-Spotted Pit viper Protobothrops mucrosquamatus

The Brown-Spotted Pit viper (Protobothrops mucrosquamatus), also known as the Chinese habu, is a widespread and highly venomous snake distributed from Northeastern India to Eastern China. Genomics research can contribute to our understanding of venom components and natural selection in vipers. Here, we collected, sequenced and assembled the genome of a male P. mucrosquamatus individual from China. We generated a highly continuous reference genome, with a length of 1.53 Gb and 41.18% of repeat elements content. Using this genome, we identified 24,799 genes, 97.97% of which could be annotated. We verified the validity of our genome assembly and annotation process by generating a phylogenetic tree based on the nuclear genome single-copy genes of six other reptile species. The results of our research will contribute to future studies on Protobothrops biology and the genetic basis of snake venom.

Genome sequence analysis of Malayan pangolin (Manis javanica) forensic samples reveals the presence of Paraburkholderia fungorum sequences

Background

The Malayan pangolin (Manis javanica) is a placental mammal and is listed as Critically Endangered on the IUCN Red List of Threatened Species. Most previous attempts to breed pangolins in captivity have met with little success because of dietary issues, infections, and other complications, although a previous study reported breeding pangolins in captivity to the third generation. In our previous pangolin genome sequencing data analysis, we obtained a considerable amount of bacterial DNA from a pregnant female Malayan pangolin (named "UM3"), which was likely infected by Paraburkholderia fungorum-an agent of biodegradation and bioremediation in agriculture.

Methodology

Here, we further confirmed and characterized this bacterial species using PCR, histological staining, whole-genome sequencing, and bioinformatics approaches. PCR assays with in-house designed primer sets and 16S universal primers showed clear positive bands in the cerebrum, cerebellum, lung, and blood of UM3 suggesting that UM3 might have developed septicaemia. Histological staining showed the presence of Gram-negative rod-shaped bacteria in the pangolin brain and lungs, indicating the colonization of the bacteria in these two organs. In addition, PCR screening of UM3's fetal tissues revealed the presence of P. fungorum in the gastrocnemius muscle, but not in other tissues that we examined. We also sequenced and reconstructed the genome of pangolin P. fungorum, which has a genome size of 7.7 Mbps.

Conclusion

Our study is the first to present detailed evidence of the presence of P. fungorum in a pangolin and her fetus (although preliminary results were presented in our previous article). Here, we raise the concern that P. fungorum may potentially infect humans, especially YOPI (young, old, pregnant, and immunocompromised) people. Therefore, caution should be exercised when using this bacterial species as biodegradation or bioremediation agents in agriculture.

Different techniques reveal the difference of community structure and function of fungi from root and rhizosphere of Salvia miltiorrhiza Bunge

Fungi have essential functions in plant health and performance. However, the plant-associated functions of many cultured fungi have not been established in detail. Here, the fungal species diversity in Salvia miltiorrhiza roots and rhizosphere was assessed for the first time using culturomics and high-throughput sequencing. We present a comprehensive functional metagenomic analysis of these fungi and verified activity of cellulase and chitinase predicted in the metagenomic analysis. We first collected and cultured fungi from the root and rhizosphere of S. miltiorrhiza. We found 92 species across 37 families and five phyla, with Ascomycota being dominant. Many rDNA internal transcribed spacer sequences could not be assigned to lower taxonomic levels. There were 19 genera of endophytic fungi and 37 genera of rhizosphere fungi. The culturomics approach had lower taxonomic diversity than high-throughput sequencing, but some fungi were only found in cultures. Structural analyses indicated that the dominant species differed in cultured and non-cultured samples at other levels, apart from the phylum level. Functional analysis mapped 223 carbohydrate enzyme families and 393 pathways in the CAZy and KEGG databases, respectively. The most abundant families were glycoside hydrolases and those involved in carbohydrate metabolism. As predicted by metagenomics, we experimentally verified cellulase and chitinase activity for 29 and 74 fungi, respectively. We provide the first evidence of biomass recycling by fungi that are associated with plants. Culturing is essential to reveal the hidden microbial community and critical functions in plant-microbe interactions.

Pseudocercospora fijiensis Conidial Germination Is Dominated by Pathogenicity Factors and Effectors

Conidia play a vital role in the survival and rapid spread of fungi. Many biological processes of conidia, such as adhesion, signal transduction, the regulation of oxidative stress, and autophagy, have been well studied. In contrast, the contribution of pathogenicity factors during the development of conidia in fungal phytopathogens has been poorly investigated. To date, few reports have centered on the pathogenicity functions of fungal phytopathogen conidia. Pseudocercospora fijiensis is a hemibiotrophic fungus and the causal agent of the black Sigatoka disease in bananas and plantains. Here, a conidial transcriptome of P. fijiensis was characterized computationally. Carbohydrates, amino acids, and lipid metabolisms presented the highest number of annotations in Gene Ontology. Common conidial functions were found, but interestingly, pathogenicity factors and effectors were also identified. Upon analysis of the resulting proteins against the Pathogen-Host Interaction (PHI) database, 754 hits were identified. WideEffHunter and EffHunter effector predictors identified 618 effectors, 265 of them were shared with the PHI database. A total of 1107 conidial functions devoted to pathogenesis were found after our analysis. Regarding the conidial effectorome, it was found to comprise 40 canonical and 578 non-canonical effectors. Effectorome characterization revealed that RXLR, LysM, and Y/F/WxC are the largest effector families in the P. fijiensis conidial effectorome. Gene Ontology classification suggests that they are involved in many biological processes and metabolisms, expanding our current knowledge of fungal effectors.

Short-term triphenyltin exposure alters microbial homeostasis in the silkworm (Bombyx mori) midgut

Triphenyltin (TPT) is a widespread synthetic chemical used in many fields and its potential risk to organisms has been comprehensively investigated using different animal models and species. Currently, little is known about the effects of TPT exposure on microbial midgut diversity, therefore we explored these effects in the lepidopterous silkworm model using 16S rDNA sequencing. In total, 5273 and 5065 operational taxonomic units (OTUs) were identified in control and TPT-exposure group samples, ranging from 424 to 728 OTUs/sample. Alpha-diversity analyses revealed that TPT exposure induced the fluctuations of gut microbial diversity and abundance while beta-diversity analyses identified a distinct impact on major gut microbiota components. In our microbiome analyses, 23 phyla and 353 genera were recognized in the control group, while 20 phyla and 358 genera were recognized in the TPT exposure group. At the genus level, midgut microbiota were composed of several predominant bacterial genera, including Muribaculaceae, Lactobacillus, and UCG-010. In the TPT exposure group, o__Bacillales, f__Bacillaceae, and f__Caldicoprobacteraceae abundance was relatively high, while f__Oscillospiraceae, f__Fusobacteriaceae, and f__SC_I_84 abundance was relatively high in the control group. Gene function analyses in silkworm microbiota after TPT exposure showed that biosynthesis of ansamycins, fructose and mannose metabolism, glycerolipid metabolism, type II diabetes mellitus, glycolysis/gluconeogenesis, lipid metabolism, translation proteins, atrazine degradation, DNA repair and recombination proteins, nicotinate and nicotinamide metabolism were significantly increased. Collectively, our silkworm model identified gut microbial diversity risks and the adverse effects from TPT exposure, which were similar to other aquatic animals. Therefore, TPT levels in environmental samples must be monitored to prevent ecological harm.

Chromosome level genome assembly of colored calla lily (Zantedeschia elliottiana)

The colored calla lily is an ornamental floral plant native to southern Africa, belonging to the Zantedeschia genus of the Araceae family. We generated a high-quality chromosome-level genome of the colored calla lily, with a size of 1,154 Mb and a contig N50 of 42 Mb. We anchored 98.5% of the contigs (1,137 Mb) into 16 pseudo-chromosomes, and identified 60.18% of the sequences (694 Mb) as repetitive sequences. Functional annotations were assigned to 95.1% of the predicted protein-coding genes (36,165). Additionally, we annotated 469 miRNAs, 1,652 tRNAs, 10,033 rRNAs, and 1,677 snRNAs. Furthermore, Gypsy-type LTR retrotransposons insertions in the genome are the primary factor causing significant genome size variation in Araceae species. This high-quality genome assembly provides valuable resources for understanding genome size differences within the Araceae family and advancing genomic research on colored calla lily.

Real-time optical detection of endodontic infection using bacterial autofluorescence

OBJECTIVES

For successful root canal treatment (RCT), it is essential to objectively assess the presence and activity of bacteria in the root canal system. However, current methods rely on subjective observations of root canal exudates. This study aimed to confirm whether real-time optical detection using bacterial autofluorescence can evaluate endodontic infection status by assessing the red fluorescence (RF) detected from root canal exudates.

METHODS

During RCT, endodontic paper points were used to collect root canal exudates scored using conventional organoleptic tests to assess the severity of root canal infections. RF on the paper points was assessed using quantitative light-induced fluorescence (QLF) technology. RF intensity and area from the paper points were quantified, and their correlations with infection severity were assessed using their organoleptic scores. The oral microbiome composition of RF samples was compared with non-red fluorescent (non-RF) samples.

RESULTS

The RF detection rate was nil and >98% in the non-infectious and severe groups. The RF intensity and area significantly increased with infection severity (p<0.001) and showed strong correlations with organoleptic scores (r=0.72, 0.82, respectively). The diagnostic accuracy for detecting root canal infection using RF intensity was good to excellent (AUC = 0.81-0.95) and increased with infection severity. The microbial diversity of the RF samples was significantly lower than that of the non-RF samples. Gram-negative anaerobic bacteria such as Prevotella and Porphyromonas were more predominant in RF samples.

CONCLUSIONS

Optical detection using bacterial autofluorescence can objectively evaluate endodontic infection status in real-time by assessing the RF of endodontic root canal exudates.

CLINICAL SIGNIFICANCE

This real-time optical technology can be utilised to detect endodontic bacterial infection without conventional incubation, allowing clinicians to determine the endpoint of chemomechanical debridement and increase the positive outcomes of RCTs.

The sequence and de novo assembly of the genome of the Indian oil sardine, Sardinella longiceps

The Indian oil sardine, Sardinella longiceps, is a widely distributed and commercially important small pelagic fish of the Northern Indian Ocean. The genome of the Indian oil sardine has been characterized using Illumina and Nanopore platforms. The assembly is 1.077 Gb (31.86 Mb Scaffold N50) in size with a repeat content of 23.24%. The BUSCO (Benchmarking Universal Single Copy Orthologues) completeness of the assembly is 93.5% when compared with Actinopterygii (ray finned fishes) data set. A total of 46316 protein coding genes were predicted. Sardinella longiceps is nutritionally rich with high levels of omega-3 polyunsaturated fatty acids (PUFA). The core genes for omega-3 PUFA biosynthesis, such as Elovl 1a and 1b,Elovl 2, Elovl 4a and 4b,Elovl 8a and 8b,and Fads 2, were observed in Sardinella longiceps. The presence of these genes may indicate the PUFA biosynthetic capability of Indian oil sardine, which needs to be confirmed functionally.

Comparison, Analysis, and Molecular Dynamics Simulations of Structures of a Viral Protein Modeled Using Various Computational Tools

The structural analysis of proteins is a major domain of biomedical research. Such analysis requires resolved three-dimensional structures of proteins. Advancements in computer technology have led to progress in biomedical research. In silico prediction and modeling approaches have facilitated the construction of protein structures, with or without structural templates. In this study, we used three neural network-based de novo modeling approaches-AlphaFold2 (AF2), Robetta-RoseTTAFold (Robetta), and transform-restrained Rosetta (trRosetta)-and two template-based tools-the Molecular Operating Environment (MOE) and iterative threading assembly refinement (I-TASSER)-to construct the structure of a viral capsid protein, hepatitis C virus core protein (HCVcp), whose structure have not been fully resolved by laboratory techniques. Templates with sufficient sequence identity for the homology modeling of complete HCVcp are currently unavailable. Therefore, we performed domain-based homology modeling for MOE simulations. The templates for each domain were obtained through sequence-based searches on NCBI and the Protein Data Bank. Then, the modeled domains were assembled to construct the complete structure of HCVcp. The full-length structure and two truncated forms modeled using various computational tools were compared. Molecular dynamics (MD) simulations were performed to refine the structures. The root mean square deviation of backbone atoms, root mean square fluctuation of Cα atoms, and radius of gyration were calculated to monitor structural changes and convergence in the simulations. The model quality was evaluated through ERRAT and phi-psi plot analysis. In terms of the initial prediction for protein modeling, Robetta and trRosetta outperformed AF2. Regarding template-based tools, MOE outperformed I-TASSER. MD simulations resulted in compactly folded protein structures, which were of good quality and theoretically accurate. Thus, the predicted structures of certain proteins must be refined to obtain reliable structural models. MD simulation is a promising tool for this purpose.

Using macromolecular electron densities to improve the enrichment of active compounds in virtual screening

The quest for effective virtual screening algorithms is hindered by the scarcity of training data, calling for innovative approaches. This study presents the use of experimental electron density (ED) data for improving active compound enrichment in virtual screening, supported by ED's ability to reflect the time-averaged behavior of ligands and solvents in the binding pocket. Experimental ED-based grid matching score (ExptGMS) was developed to score compounds by measuring the degree of matching between their binding conformations and a series of multi-resolution experimental ED grids. The efficiency of ExptGMS was validated using both in silico tests with the Directory of Useful Decoys-Enhanced dataset and wet-lab tests on Covid-19 3CLpro-inhibitors. ExptGMS improved the active compound enrichment in top-ranked molecules by approximately 20%. Furthermore, ExptGMS identified four active inhibitors of 3CLpro, with the most effective showing an IC50 value of 1.9 µM. We also developed an online database containing experimental ED grids for over 17,000 proteins to facilitate the use of ExptGMS for academic users.

Genome Sequencing Provides Novel Insights into Mudflat Burrowing Adaptations in Eel Goby Taenioides sp. (Teleost: Amblyopinae)

Amblyopinae is one of the lineage of bony fish that preserves amphibious traits living in tidal mudflat habitats. In contrast to other active amphibious fish, Amblyopinae species adopt a seemly more passive lifestyle by living in deep burrows of mudflat to circumvent the typical negative effects associated with terrestriality. However, little is known about the genetic origin of these mudflat deep-burrowing adaptations in Amblyopinae. Here we sequenced the first genome of Amblyopinae species, Taenioides sp., to elucidate their mudflat deep-burrowing adaptations. Our results revealed an assembled genome size of 774.06 Mb with 23 pseudochromosomes anchored, which predicted 22,399 protein-coding genes. Phylogenetic analyses indicated that Taenioides sp. diverged from the active amphibious fish of mudskipper approximately 28.3 Ma ago. In addition, 185 and 977 putative gene families were identified to be under expansion, contraction and 172 genes were undergone positive selection in Taenioides sp., respectively. Enrichment categories of top candidate genes under significant expansion and selection were mainly associated with hematopoiesis or angiogenesis, DNA repairs and the immune response, possibly suggesting their involvement in the adaptation to the hypoxia and diverse pathogens typically observed in mudflat burrowing environments. Some carbohydrate/lipid metabolism, and insulin signaling genes were also remarkably alterated, illustrating physiological remolding associated with nutrient-limited subterranean environments. Interestingly, several genes related to visual perception (e.g., crystallins) have undergone apparent gene losses, pointing to their role in the small vestigial eyes development in Taenioides sp. Our work provide valuable resources for understanding the molecular mechanisms underlying mudflat deep-burrowing adaptations in Amblyopinae, as well as in other tidal burrowing teleosts.

High subtelomeric GC content in the genome of a zoonotic Cryptosporidium species

Cryptosporidium canis is a zoonotic species causing cryptosporidiosis in humans in addition to its natural hosts dogs and other fur animals. To understand the genetic basis for host adaptation, we sequenced the genomes of C. canis from dogs, minks, and foxes and conducted a comparative genomics analysis. While the genomes of C. canis have similar gene contents and organisations, they (~41.0 %) and C. felis (39.6 %) have GC content much higher than other Cryptosporidium spp. (24.3-32.9 %) sequenced to date. The high GC content is mostly restricted to subtelomeric regions of the eight chromosomes. Most of these GC-balanced genes encode Cryptosporidium-specific proteins that have intrinsically disordered regions and are involved in host-parasite interactions. Natural selection appears to play a more important role in the evolution of codon usage in GC-balanced C. canis, and most of the GC-balanced genes have undergone positive selection. While the identity in whole genome sequences between the mink- and dog-derived isolates is 99.9 % (9365 SNVs), it is only 96.0 % (362 894 SNVs) between them and the fox-derived isolate. In agreement with this, the fox-derived isolate possesses more subtelomeric genes encoding invasion-related protein families. Therefore, the change in subtelomeric GC content appears to be responsible for the more GC-balanced C. canis genomes, and the fox-derived isolate could represent a new Cryptosporidium species.

Complete Genome Sequence of Bacteriophage EO1, Which Infects Both Escherichia coli O157:H7 and Shigella

The lytic bacteriophage EO1 has been newly isolated. This phage infects Escherichia coli O157:H7 and has a broad antibacterial spectrum, including against Shigella. The complete genome sequence of phage EO1 was determined; its full length is 166,941 bp, and it has a G+C content of 35.46%.

Heavy Metal Pollution Impacts Soil Bacterial Community Structure and Antimicrobial Resistance at the Birmingham 35th Avenue Superfund Site

Heavy metals (HMs) are known to modify bacterial communities both in the laboratory and in situ. Consequently, soils in HM-contaminated sites such as the U.S. Environmental Protection Agency (EPA) Superfund sites are predicted to have altered ecosystem functioning, with potential ramifications for the health of organisms, including humans, that live nearby. Further, several studies have shown that heavy metal-resistant (HMR) bacteria often also display antimicrobial resistance (AMR), and therefore HM-contaminated soils could potentially act as reservoirs that could disseminate AMR genes into human-associated pathogenic bacteria. To explore this possibility, topsoil samples were collected from six public locations in the zip code 35207 (the home of the North Birmingham 35th Avenue Superfund Site) and in six public areas in the neighboring zip code, 35214. 35027 soils had significantly elevated levels of the HMs As, Mn, Pb, and Zn, and sequencing of the V4 region of the bacterial 16S rRNA gene revealed that elevated HM concentrations correlated with reduced microbial diversity and altered community structure. While there was no difference between zip codes in the proportion of total culturable HMR bacteria, bacterial isolates with HMR almost always also exhibited AMR. Metagenomes inferred using PICRUSt2 also predicted significantly higher mean relative frequencies in 35207 for several AMR genes related to both specific and broad-spectrum AMR phenotypes. Together, these results support the hypothesis that chronic HM pollution alters the soil bacterial community structure in ecologically meaningful ways and may also select for bacteria with increased potential to contribute to AMR in human disease. IMPORTANCE Heavy metals cross-select for antimicrobial resistance in laboratory experiments, but few studies have documented this effect in polluted soils. Moreover, despite decades of awareness of heavy metal contamination at the EPA Superfund site in North Birmingham, Alabama, this is the first analysis of the impact of this pollution on the soil microbiome. Specifically, this work advances the understanding of the relationship between heavy metals, microbial diversity, and patterns of antibiotic resistance in North Birmingham soils. Our results suggest that polluted soils carry a risk of increased exposure to antibiotic-resistant infections in addition to the direct health consequences of heavy metals. Our work provides important information relevant to both political and scientific efforts to advance environmental justice for the communities that call Superfund neighborhoods home.

Tandem mass tag-based quantitative proteomics analyses of a chicken-original virulent and its attenuated Histomonas meleagridis strain in China

Introduction

Histomonas meleagridis can cause histomonosis in poultry. Due to the prohibition of effective drugs, the prevention and treatment of the disease requires new strategies. Questions about its pathogenic mechanisms and virulence factors remain puzzling.

Methods

To address these issues, a tandem mass tag (TMT) comparative proteomic analysis of a virulent strain and its attenuated strain of Chinese chicken-origin was performed.

Results

A total of 3,494 proteins were identified in the experiment, of which 745 proteins were differentially expressed (fold change ≥1.2 or ≤0.83 and p &lt; 0.05), with 192 up-regulated proteins and 553 down-regulated proteins in the virulent strain relative to the attenuated strain.

Discussion

Surface protein BspA like, digestive cysteine proteinase, actin, and GH family 25 lysozyme were noted among the proteins up regulated in virulent strains, and these several proteins may be directly related to the pathogenic capacity of the histomonad. Ferredoxin, 60S ribosomal protein L6, 40S ribosomal protein S3, and NADP-dependent malic enzyme which associated with biosynthesis and metabolism were also noted, which have the potential to be new drug targets. The up-regulation of alpha-amylase, ras-like protein 1, ras-like protein 2, and involucrin in attenuated strains helps to understand how it is adapted to the long-term in vitro culture environment. The above results provide some candidate protein-coding genes for further functional verification, which will help to understand the molecular mechanism of pathogenicity and attenuation of H. meleagridis more comprehensively.

Fire frequency, as well as stress response and developmental gene control serotiny level variation in a widespread pioneer Mediterranean conifer, Pinus halepensis

Many plants undergo adaptation to fire. Yet, as global change is increasing fire frequency worldwide, our understanding of the genetics of adaptation to fire is still limited. We studied the genetic basis of serotiny (the ability to disseminate seeds exclusively after fire) in the widespread, pioneer Mediterranean conifer Pinus halepensis Mill., by linking individual variation in serotiny presence and level to fire frequency and to genetic polymorphism in natural populations. After filtering steps, 885 single nucleotide polymorphisms (SNPs) out of 8000 SNPs used for genotyping were implemented to perform an in situ association study between genotypes and serotiny presence and level. To identify serotiny-associated loci, we performed random forest analyses of the effect of SNPs on serotiny levels, while controlling for tree size, frequency of wildfires, and background environmental parameters. Serotiny showed a bimodal distribution, with serotinous trees more frequent in populations exposed to fire in their recent history. Twenty-two SNPs found in genes involved in stress tolerance were associated with the presence-absence of serotiny while 37 found in genes controlling for flowering were associated with continuous serotiny variation. This study shows the high potential of P. halepensis to adapt to changing fire regimes, benefiting from a large and flexible genetic basis of trait variation.

Chromosome-scale genomes reveal genomic consequences of inbreeding in the South China tiger: A comparative study with the Amur tiger

The South China tiger (Panthera tigris amoyensis, SCT) is the most critically endangered subspecies of tiger due to functional extinction in the wild. Inbreeding depression is observed among the captive population descended from six wild ancestors, resulting in high juvenile mortality and low reproduction. We assembled and characterized the first SCT genome and an improved Amur tiger (P. t. altaica, AT) genome named AmyTig1.0 and PanTig2.0. The two genomes are the most continuous and comprehensive among any tiger genomes yet reported at the chromosomal level. By using the two genomes and resequencing data of 15 SCT and 13 AT individuals, we investigated the genomic signature of inbreeding depression of the SCT. The results indicated that the effective population size of SCT experienced three phases of decline, ~5.0-1.0 thousand years ago, 100 years ago, and since captive breeding in 1963. We found 43 long runs of homozygosity fragments that were shared by all individuals in the SCT population and covered a total length of 20.63% in the SCT genome. We also detected a large proportion of identical-by-descent segments across the genome in the SCT population, especially on ChrB4. Deleterious nonsynonymous single nucleotide polymorphic sites and loss-of-function mutations were found across genomes with extensive potential influences, despite a proportion of these loads having been purged by inbreeding depression. Our research provides an invaluable resource for the formulation of genetic management policies for the South China tiger such as developing genome-based breeding and genetic rescue strategy.

Storage Fungi and Mycotoxins Associated with Rice Samples Commercialized in Thailand

The study focused on the examination of the different fungal species isolated from commercial rice samples, applying conventional culture techniques, as well as different molecular and phylogenic analyses to confirm phenotypic identification. Additionally, the mycotoxin production and contamination were analyzed using validated liquid chromatography-tandem mass spectrometry (LC-MS/MS). In total, 40 rice samples were obtained covering rice berry, red jasmine rice, brown rice, germinated brown rice, and white rice. The blotting paper technique applied on the 5 different types of rice samples detected 4285 seed-borne fungal infections (26.8%) for 16,000 rice grains. Gross morphological data revealed that 19 fungal isolates belonged to the genera Penicillium/Talaromyces (18 of 90 isolates; 20%) and Aspergillus (72 of 90 isolates; 80%). To check their morphologies, molecular data (fungal sequence-based BLAST results and a phylogenetic tree of the combined ITS, BenA, CaM, and RPB2 datasets) confirmed the initial classification. The phylogenic analysis revealed that eight isolates belonged to P. citrinum and, additionally, one isolate each belonged to P. chermesinum, A. niger, A. fumigatus, and A. tubingensis. Furthermore, four isolates of T. pinophilus and one isolate of each taxon were identified as Talaromyces (T. radicus, T. purpureogenum, and T. islandicus). The results showed that A. niger and T. pinophilus were two commonly occurring fungal species in rice samples. After subculturing, ochratoxin A (OTA), generated by T. pinophilus code W3-04, was discovered using LC-MS/MS. In addition, the Fusarium toxin beauvericin was detected in one of the samples. Aflatoxin B1 or other mycotoxins, such as citrinin, trichothecenes, and fumonisins, were detected. These preliminary findings should provide valuable guidance for hazard analysis critical control point concepts used by commercial food suppliers, including the analysis of multiple mycotoxins. Based on the current findings, mycotoxin analyses should focus on A. niger toxins, including OTA and metabolites of T. pinophilus (recently considered a producer of emerging mycotoxins) to exclude health hazards related to the traditionally high consumption of rice by Thai people.