High-quality Genome Research Articles

Whole-genome sequencing and assessment of a novel protein- and gossypol-degrading Bacillus subtilis strain isolated from intestinal digesta of Tibetan Pigs

BackgroundWith the rapid development of animal husbandry, the demand for protein feed resources is increasing. Cottonseed meal (CSM) and soybean meal (SBM) are rich sources of protein. However, their application is limited due to the existence of anti-nutrients, which can be harmful to the digestion and absorption. A strain of Bacillus subtilis (Mafic-Y7) was isolated from digesta of intestines of Tibetan pigs. The strain showed high protease activity, which helps in degrading proteinic anti-nutritional factors in grain meal and in vitro degradation of free gossypol. In order to better understand this isolated strain, whole genome of Mafic-Y7 strain was sequenced and analyzed. Different effects on various grain meals were identified.ResultThe GC-depth Poisson distributions showed no bias suggesting high-quality genome assembly of Mafic-Y7. The whole genome sequencing showed that one chromosome with 4,248,845 base pairs(bp)and the genes total length with 3,736,524 bp was predicted in Mafic-Y7. Additionally, Mafic-Y7 possessed 4,254 protein-coding genes, and several protease genes were annotated by aligning them with databases. There are 55 protease genes, one phytase gene and one laccase gene were annotated in the gene sequence of Mafic-Y7. The average nucleotide identity between Mafic-Y7 and the GCA-000009045.1 homologous genome was 0.9938, suggesting a close genetic relationship between them at the species level. Compared with the closest four whole genomes, Mafic-Y7 was annotated the most abundant of protease genes (55 genes). The fermentation supernatant of Mafic-Y7 could increase the content of small peptides, water-soluble proteins, and acid-soluble proteins in vitro by 411%, 281% and 317% in SBM and 420%, 257% and 338% in CSM. After fermentation in grain meal by Mafic-Y7, the degradation rate of anti-nutritional factors in SBM, such as trypsin inhibitor, glycinin, and β-conglycinin was greater than 70%, and lectin was greater than 30%. The degradation rates of anti-nutritional factors in CSM, such as gossypol and phytic acid, were 82% and 26%, respectively.

Long-read genome assembly of the insect model organism Tribolium castaneum reveals spread of satellite DNA in gene-rich regions by recurrent burst events.

Eukaryotic genomes are replete with satellite DNAs (satDNAs), large stretches of tandemly repeated sequences that are mostly underrepresented in genome assemblies. Here we combined nanopore long-read sequencing with a reference-guided assembly approach to generate an improved, high-quality genome assembly, TcasONT, of the model beetle Tribolium castaneum Enriched by 45 Mb in repetitive regions, the new assembly comprises almost the entire genome sequence. We use the enhanced assembly to conduct global and in-depth analyses of abundant euchromatic satDNAs. Unexpectedly, we show the extensive spread of satDNAs in gene-rich regions, including long arrays. The sequence similarity relationships between satDNA monomers and arrays indicate a recent exchange of satDNA arrays between different chromosomes. We propose a scenario of their genome dynamics characterized by repeated bursts of satDNAs spreading through euchromatin, followed by a process of elongation and homogenization of arrays. We find that suppressed recombination on the X Chromosome has no significant effect on the spread of satDNAs but the X rather tolerates the amplification of satDNAs into longer arrays. Analyses of arrays' neighboring regions show a tendency of one satDNA to be associated with transposable-like elements. Using 2D electrophoresis followed by Southern blotting, we prove Cast satDNAs' presence in the fraction of extrachromosomal circular DNA (eccDNA). We point to two mechanisms that enable this satDNA spread to occur: transposition by transposable elements and insertion mediated by eccDNA. The presence of such a large proportion of satDNA in gene-rich regions inevitably gives rise to speculation about their possible influence on gene expression.

Haemoglobin Gene Repertoire in Teleost and Cichlid Fishes Shaped by Gene Duplications and Genome Rearrangements.

Haemoglobin is a key molecule for oxygen transport in vertebrates. It exhibits remarkable gene diversity in teleost fishes, reflecting adaptation to various aquatic environments. In this study, we present the dynamic evolution of haemoglobin subunit genes based on a comparison of high-quality genome assemblies of 24 vertebrate species, including 17 teleosts (of which six are cichlids). Our findings indicate that teleost genomes contain a range of haemoglobin genes, from as few as five in fugu to as many as 43 in salmon, with the latter being the largest repertoire found in vertebrates. We find evidence that the teleost ancestor had at least four Hbα and three or four Hbβ subunit genes, and that the current gene diversity emerged during teleost radiation, driven primarily by (tandem) gene duplications, genome compaction, and rearrangement dynamics. We provide insights into the genomic organisation of haemoglobin clusters in different teleost species. We further show that the evolution of paralogous rhbdf1 genes flanking both teleost clusters (LA and MN) supports the hypothesis for the origin of the LA cluster by rearrangement within teleosts, rather than by the teleost specific whole-genome duplication. We specifically focus on cichlid fishes, where adaptation to low oxygen environment plays role in species diversification. Our analysis of six cichlid genomes, including Pungu maclareni from the Barombi Mbo crater lake, for which we sequenced a representative genome, reveals 18-32 copies of the Hb genes, and elevated rates of non-synonymous substitutions compared to other teleosts. Overall, this work facilitates a deeper understanding of how haemoglobin genes contribute to the adaptive potential of teleosts.

A chromosome-level genome assembly for the desert scorpion Mesobuthus przewalskii from Asian drylands.

Scorpions are an ancient and charismatic group of arthropods with medical importance, but a high-quality reference genome for this group is still lacking. Here we perform whole-genome sequencing of Mesobuthus przewalskii, a desert scorpion endemic to the Taklimakan Desert. We combine PacBio HiFi sequencing and Hi-C chromosome conformation capturing to generate chromosomal-level, haplotype-resolved, and fully annotated genome assembly for this medically important scorpion. The assembly consists of two haplotypes (1052.01 Mbp and 1055.19 Mbp, respectively) reaching chromosome-level contiguity and >98% BUSCO completeness. Sequences were anchored in 13 chromosomes with a contig N50 of 34.44 Mbp and scaffold N50 of 81.43 Mbp. Several key genome features and the mitochondrial genome assembly were also provided. This genome represents the fifth but the most complete assembly for the order Scorpiones.

Seasonal variations of microbial communities and viral diversity in fishery-enhanced marine ranching sediments: insights into metabolic potentials and ecological interactions

BackgroundThe ecosystems of marine ranching have enhanced marine biodiversity and ecological balance and have promoted the natural recovery and enhancement of fishery resources. The microbial communities of these ecosystems, including bacteria, fungi, protists, and viruses, are the drivers of biogeochemical cycles. Although seasonal changes in microbial communities are critical for ecosystem functioning, the current understanding of microbial-driven metabolic properties and their viral communities in marine sediments remains limited. Here, we employed amplicon (16S and 18S) and metagenomic approaches aiming to reveal the seasonal patterns of microbial communities, bacterial-eukaryotic interactions, whole metabolic potential, and their coupling mechanisms with carbon (C), nitrogen (N), and sulfur (S) cycling in marine ranching sediments. Additionally, the characterization and diversity of viral communities in different seasons were explored in marine ranching sediments.ResultsThe current study demonstrated that seasonal variations dramatically affected the diversity of microbial communities in marine ranching sediments and the bacterial-eukaryotic interkingdom co-occurrence networks. Metabolic reconstruction of the 113 medium to high-quality metagenome-assembled genomes (MAGs) was conducted, and a total of 8 MAGs involved in key metabolic genes and pathways (methane oxidation - denitrification - S oxidation), suggesting a possible coupling effect between the C, N, and S cycles. In total, 338 viral operational taxonomic units (vOTUs) were identified, all possessing specific ecological characteristics in different seasons and primarily belonging to Caudoviricetes, revealing their widespread distribution and variety in marine sediment ecosystems. In addition, predicted virus-host linkages showed that high host specificity was observed, with few viruses associated with specific hosts.ConclusionsThis finding deepens our knowledge of element cycling and viral diversity in fisheries enrichment ecosystems, providing insights into microbial-virus interactions in marine sediments and their effects on biogeochemical cycling. These findings have potential applications in marine ranching management and ecological conservation.DRyDF7ykR-XALqAtu34uZtVideo

Complete genome sequence of Kozakia baliensis LTH 7215, a newly isolated strain from grape must.

We report the complete genome sequence of the novel Kozakia baliensis strain LTH 7215, isolated from grape must. This strain exhibits a mucoid phenotype on solid medium, associated with exopolysaccharide production. The high-quality and full-length genome sequence was obtained using PacBio Sequencing technology, refined with Illumina sequencing.

PneumoBrowse 2: an integrated visual platform for curated genome annotation and multiomics data analysis of Streptococcus pneumoniae.

Streptococcus pneumoniae is an opportunistic human pathogen responsible for high morbidity and mortality rates. Extensive genome sequencing revealed its large pangenome, serotype diversity, and provided insight into genome dynamics. However, functional genome analysis has lagged behind, as that requires detailed and time-consuming manual curation of genome annotations and integration of genomic and phenotypic data. To remedy this, PneumoBrowse was presented in 2018, a user-friendly interactive online platform, which provided the detailed annotation of the S. pneumoniae D39V genome, alongside transcriptomic data. Since 2018, many new studies on S. pneumoniae genome biology and protein functioning have been performed. Here, we present PneumoBrowse 2 (https://veeninglab.com/pneumobrowse), fully rebuilt in JBrowse 2. We updated annotations for transcribed and transcriptional regulatory features in the D39V genome. We added genome-wide data tracks for high-resolution chromosome conformation capture (Hi-C) data, chromatin immunoprecipitation coupled to high-throughput sequencing (ChIP-Seq), ribosome profiling, CRISPRi-seq gene essentiality data and more. Additionally, we included 18 phylogenetically diverse S. pneumoniae genomes and their annotations. By providing easy access to diverse high-quality genome annotations and links to other databases (including UniProt and AlphaFold), PneumoBrowse 2 will further accelerate research and development into preventive and treatment strategies, through increased understanding of the pneumococcal genome.

GCI: a continuity inspector for complete genome assembly.

Recent advances in long-read sequencing technologies have significantly facilitated the production of high-quality genome assembly. The telomere-to-telomere (T2T) gapless assembly has become the new golden standard of genome assembly efforts. Several recent efforts have claimed to produce T2T-level reference genomes. However, a universal standard is still missing to qualify a genome assembly to be at T2T standard. Traditional genome assembly assessment metrics (N50 and its derivatives) have no capacity in differentiating between nearly T2T assembly and the truly T2T assembly in continuity either globally or locally. Additionally, these metrics are independent of raw reads, making them inflated easily by artificial operations. Therefore, a gaplessness evaluation tool at single-nucleotide resolution to reflect true completeness is urgently needed in the era of complete genomes. Here, we present a tool called Genome Continuity Inspector (GCI), designed to assess genome assembly continuity at single-base resolution, and evaluate how close an assembly is to the T2T level. GCI utilizes multiple aligners to map long reads from various sequencing platforms back to the assembly. By incorporating curated mapping coverage of high-confidence read alignments, GCI identifies potential assembly issues. Meanwhile, it provides GCI scores that quantify overall assembly continuity on the whole genome or chromosome scales. The open-source GCI code is freely available on Github (https://github.com/yeeus/GCI) under the MIT license.

Chromosome-level genome assembly of the cosmopolitan diatom Skeletonema costatum provides insights into ecological adaptation

The cosmopolitan diatom species Skeletonema costatum is an ecologically important dominant phytoplankton frequently found in the coastal estuarine and marine waters, and often causes harmful algae blooms. Despite of its critical ecological importance, chromosome-level genome assemble is still unavailable, hindering in-depth understanding of their evolution and environmental adaption. Here, we report a chromosome-level genome assembly for the marine diatom species S. costatum. The assembled genome size was 136.49 Mb, with a contig N50 of 302 Kb and 95.30 % of the reads anchored into 23 pseudo-chromosomes with a scaffold N50 of 6.19 Mb. A total of 28,321 protein-coding genes were predicted, with 86.03 % being functional annotated. The BUSCO assessment of genome assembly and genome annotation were both above 90 %. Phylogenetic analysis showed the expected topology, with S. costatum and its closely related species S. marinoi diverged from their common ancestor around 22.6 million years ago. The genome size of S. costatum is comparatively larger than those of its closely related diatoms, due mostly to its higher transposable element contents and larger number of proteincoding genes. Collinearity analysis revealed strong collinearity between S. costatum and other Skeletonema with most chromosomes showing clear one-to-one correspondences. A larger family of nine copies of the cryptochrome genes that function as blue light photoreceptors were identified in S. costatum, which could contribute its ecological success. The availability of the high-quality chromosome-level genome assembly for S. costatum represents a valuable resource that may facilitate comparative genomics for revealing important ecological clues and gene families, and future genetics and environmental studies among Skeletonema species.

Population-wide DNA methylation polymorphisms at single-nucleotide resolution in 207 cotton accessions reveal epigenomic contributions to complex traits.

DNA methylation plays multiple regulatory roles in crop development. However, the relationships of methylation polymorphisms with genetic polymorphisms, gene expression, and phenotypic variation in natural crop populations remain largely unknown. Here, we surveyed high-quality methylomes, transcriptomes, and genomes obtained from the 20-days-post-anthesis (DPA) cotton fibers of 207 accessions and extended the classical framework of population genetics to epigenetics. Over 287 million single methylation polymorphisms (SMPs) were identified, 100 times more than the number of single nucleotide polymorphisms (SNPs). These SMPs were significantly enriched in intragenic regions while depleted in transposable elements. Association analysis further identified a total of 5,426,782 cis-methylation quantitative trait loci (cis-meQTLs), 5078 cis-expression quantitative trait methylation (cis-eQTMs), and 9157 expression quantitative trait loci (eQTLs). Notably, 36.39% of cis-eQTM genes were not associated with genetic variation, indicating that a large number of SMPs associated with gene expression variation are independent of SNPs. In addition, out of the 1715 epigenetic loci associated with yield and fiber quality traits, only 36 (2.10%) were shared with genome-wide association study (GWAS) loci. The construction of multi-omics regulatory networks revealed 43 cis-eQTM genes potentially involved in fiber development, which cannot be identified by GWAS alone. Among these genes, the role of one encoding CBL-interacting protein kinase 10 in fiber length regulation was successfully validated through gene editing. Taken together, our findings prove that DNA methylation data can serve as an additional resource for breeding purposes and can offer opportunities to enhance and expedite the crop improvement process.

Activity of novel virus families infecting soil nitrifiers is concomitant with host niche differentiation.

Chemolithoautotrophic nitrifiers are model groups for linking phylogeny, evolution, and ecophysiology. Ammonia-oxidising bacteria (AOB) typically dominate the first step of ammonia oxidation at high ammonium supply rates, ammonia-oxidising archaea (AOA) and complete ammonia-oxidising Nitrospira (comammox) are often active at lower supply rates or during AOB inactivity, and nitrite-oxidising bacteria (NOB) complete canonical nitrification. Soil virus communities are dynamic but contributions to functional processes are largely undetermined. In addition, characterising viruses infecting hosts with low relative abundance, such as nitrifiers, may be constrained by vast viral diversity, partial genome recovery, and difficulties in host linkage. Here, we describe a targeted incubation study that aimed to determine whether growth of different nitrifier groups in soil is associated with active virus populations and if process-focussed analyses facilitate characterisation of high-quality virus genomes. dsDNA viruses infecting different nitrifier groups were enriched in situ via differential host inhibition. Growth of each nitrifier group was consistent with predicted inhibition profiles and concomitant with the abundance of their viruses. These included 61 high-quality/complete virus genomes 35-173kb in length with minimal similarity to validated families. AOA viruses lacked ammonia monooxygenase sub-unit C (amoC) genes found in marine AOA viruses but some encoded AOA-specific multicopper oxidase type 1 (MCO1), previously implicated in copper acquisition, and suggesting a role in supporting energy metabolism of soil AOA. Findings demonstrate focussed incubation studies facilitate characterisation of active host-virus interactions associated with specific processes and viruses of soil AOA, AOB and NOB are dynamic and potentially influence nitrogen cycling processes.

The genome and comparative transcriptome of the euryhaline model ciliate Paramecium duboscqui reveal adaptations to environmental salinity

BackgroundAs a potential model organism for studies of environmental and cell biology, Paramecium duboscqui is a special euryhaline species of Paramecium that can be found in fresh, brackish, or marine water in natural salinity ranges between 0‰ and 33‰. However, the genome information as well as molecular mechanisms that account for its remarkable halotolerant traits remain extremely unknown. To characterize its genome feature, we combined PacBio and Illumina sequencing to assemble the first high-quality and near-complete macronuclear genome of P. duboscqui. Meanwhile, comparative transcriptomic profiles under different salinities gave underlying insight into the molecular mechanism of its adaptations to environmental salinity.ResultsThe results showed that the MAC genome of P. duboscqui comprises 160 contigs, with 113 of them possessing telomere (~ 28.82 Mb haploid genome size). Through comparative genomic analyses with the other ciliate, we found that gene families encoding transmembrane transporter proteins have been expanded in P. duboscqui, showing enormous potential in salinity adaptation. Like other Paramecium, P. duboscqui utilizes TGA as its only termination codon and has reassigned TAA and TAG to encode glutamine. P. duboscqui showed different growth rates under different salinities, with an optimum growth rate in 5‰ salinity. A comparison of the transcriptomic profiles among P. duboscqui grown under different concentrations showed that genes involved in protein folding, oxygen respiration, and glutathione-dependent detoxification were upregulated in the high-salt group, whereas genes encoding DNA-binding proteins and transcription factors were upregulated in the low-salt group, suggesting distinct mechanisms for responding to low and high salinity. Weighted gene coexpression network analysis (WGCNA) linked the hub genes expressed at 30‰ salinity to cysteine-type peptidase activity, lipid transfer, sodium hydrogen exchange, and cell division, with the hub genes expressed at 0‰ salinity involved in transmembrane transport and protein localization.ConclusionsThis study characterizes a new euryhaline model Paramecium, provides novel insights into Paramecium evolution, and describes the molecular mechanisms that have allow P. duboscqui to adapt to different osmotic environments.

Global characterization of somatic mutations and DNA methylation changes during vegetative propagation in strawberries.

Somatic mutations arise and accumulate during tissue culture and vegetative propagation, potentially affecting various traits in horticultural crops, but their characteristics are still unclear. Here, somatic mutations in regenerated woodland strawberry derived from tissue culture of shoot tips under different conditions and 12 cultivated strawberry individuals are analyzed by whole genome sequencing. The mutation frequency of single nucleotide variants is significantly increased with increased hormone levels or prolonged culture time in the range of 3.3 × 10-8-3.0 × 10-6 mutations per site. CG methylation shows a stable reduction (0.71%-8.03%) in regenerated plants, and hypoCG-DMRs are more heritable after sexual reproduction. A high-quality haplotype-resolved genome is assembled for the strawberry cultivar "Beni hoppe." The 12 "Beni hoppe" individuals randomly selected from different locations show 4731-6005 mutations relative to the reference genome, and the mutation frequency varies among the subgenomes. Our study has systematically characterized the genetic and epigenetic variants in regenerated woodland strawberry plants and different individuals of the same strawberry cultivar, providing an accurate assessment of somatic mutations at the genomic scale and nucleotide resolution in plants.

Genetic analysis of a Yayoi individual from the Doigahama site provides insights into the origins of immigrants to the Japanese Archipelago.

Mainland Japanese have been recognized as having dual ancestry, originating from indigenous Jomon people and immigrants from continental East Eurasia. Although migration from the continent to the Japanese Archipelago continued from the Yayoi to the Kofun period, our understanding of these immigrants, particularly their origins, remains insufficient due to the lack of high-quality genome samples from the Yayoi period, complicating predictions about the admixture process. To address this, we sequenced the whole nuclear genome of a Yayoi individual from the Doigahama site in Yamaguchi prefecture, Japan. A comprehensive population genetic analysis of the Doigahama Yayoi individual, along with ancient and modern populations in East Asia and Northeastern Eurasia, revealed that the Doigahama Yayoi individual, similar to Kofun individuals and modern Mainland Japanese, had three distinct genetic ancestries: Jomon-related, East Asian-related, and Northeastern Siberian-related. Among non-Japanese populations, the Korean population, possessing both East Asian-related and Northeastern Siberian-related ancestries, exhibited the highest degree of genetic similarity to the Doigahama Yayoi individual. The analysis of admixture modeling for Yayoi individuals, Kofun individuals, and modern Japanese respectively supported a two-way admixture model assuming Jomon-related and Korean-related ancestries. These results suggest that between the Yayoi and Kofun periods, the majority of immigrants to the Japanese Archipelago originated primarily from the Korean Peninsula.

What lies behind the large genome of Colletotrichum lindemuthianum.

Colletotrichum lindemuthianum is the etiological agent of anthracnose disease in common bean (Phaseolus vulgaris L.), noted for its ability to cause serious damage and significant pathogenic variability. This study reveals the features of the high-quality genome of C. lindemuthianum. Analysis showed improvements over the first assembly, with the refined genome having 119 scaffolds, ten times fewer than the first, and a 19% increase in gene number. The effector candidates increased by nearly 1.5 times. More than 40% of the amino acid sequences with homologs in the Pathogen-Host Interactions (PHI-base) database are linked to pathogenicity. Of 18 putative proteins identified as Chitinase-like Protein, six have a mutation in the enzyme catalytic motif, and three showed gene expression in the biotrophic phase, indicating they can act as effectors. Comparative genomic analyses with 30 other fungal species revealed that C. lindemuthianum is among the top three fungi encoding transport proteins. Seven Necrosis and Ethylene-Inducing Peptide 1 (Nep1)-Like Proteins (NLPs) are present in the C. lindemuthianum genome, but none had complete identity with the GHRHDWE conserved motif of NLPs; two were grouped with proteins that induce necrosis and may retain the capability to induce host necrosis. Colletotrichum species show a high number of secondary metabolite (SM) clusters, with C. lindemuthianum having 47 SM clusters. Approximately 60% of the C. lindemuthianum genome is composed of repetitive elements, a significantly higher proportion than in other fungi. These differences in transposable element (TE) numbers may explain why C. lindemuthianum has one of the largest genomes among the fungi analyzed. A significant portion of its genome comprises retroelements, particularly the Ty1/Copia superfamily, which accounts for 22% of the genome and represents 40% of the repetitive elements. The genomic profile features a remarkably high RIP-affected genomic proportion of 54.77%, indicating substantial RIP activity within this species. This high-quality genome of C. lindemuthianum, a significant pathogen in common bean cultivation, will support future research into this pathosystem, fostering a deeper understanding of the interaction between the fungus and its host.

Population Genomics of Adaptive Radiations: Exceptionally High Levels of Genetic Diversity and Recombination in an Endemic Spider From the Canary Islands.

The spider genus Dysdera has undergone a remarkable diversification in the oceanic archipelago of the Canary Islands, with ~60 endemic species having originated during the 20 million years since the origin of the archipelago. This evolutionary radiation has been accompanied by substantial dietary shifts, often characterised by phenotypic modifications encompassing morphological, metabolic and behavioural changes. Hence, these endemic spiders represent an excellent model for understanding the evolutionary drivers and to pinpoint the genomic determinants underlying adaptive radiations. Recently, we achieved the first chromosome-level genome assembly of one of the endemic species, D. silvatica, providing a high-quality reference sequence for evolutionary genomics studies. Here, we conducted a low coverage-based resequencing study of a natural population of D. silvatica from La Gomera island. Taking advantage of the new high-quality genome, we characterised genome-wide levels of nucleotide polymorphism, divergence and linkage disequilibrium, and inferred the demographic history of this population. We also performed comprehensive genome-wide scans for recent positive selection. Our findings uncovered exceptionally high levels of nucleotide diversity and recombination in this geographically restricted endemic species, indicative of large historical effective population sizes. We also identified several candidate genomic regions that are potentially under positive selection, highlighting relevant biological processes, such as vision and nitrogen extraction as potential adaptation targets. These processes may ultimately drive species diversification in this genus. This pioneering study of spiders that are endemic to an oceanic archipelago lays the groundwork for broader population genomics analyses aimed at understanding the genetic mechanisms driving adaptive radiation in island ecosystems.

High-quality draft genome of Lactiplantibacillus plantarum strain APC2688 isolated from human feces.

Lactiplantibacillus plantarum is a commensal bacterial species often found within the human gut. To expand our knowledge about this human-associated taxon with potential probiotic properties, here we present the draft genome sequence of Lactiplantibacillus plantarum strain APC2688, isolated from a human fecal sample.

Genomic and functional characterization of the Atlantic salmon gut microbiome in relation to nutrition and health.

To ensure sustainable aquaculture, it is essential to understand the path 'from feed to fish', whereby the gut microbiome plays an important role in digestion and metabolism, ultimately influencing host health and growth. Previous work has reported the taxonomic composition of the Atlantic salmon (Salmo salar) gut microbiome; however, functional insights are lacking. Here we present the Salmon Microbial Genome Atlas consisting of 211 high-quality bacterial genomes, recovered by cultivation (n = 131) and gut metagenomics (n = 80) from wild and farmed fish both in freshwater and seawater. Bacterial genomes were taxonomically assigned to 14 different orders, including 35 distinctive genera and 29 previously undescribed species. Using metatranscriptomics, we functionally characterized key bacterial populations, across five phyla, in the salmon gut. This included the ability to degrade diet-derived fibres and release vitamins and other exometabolites with known beneficial effects, which was supported by genome-scale metabolic modelling and in vitro cultivation of selected bacterial species coupled with untargeted metabolomic studies. Together, the Salmon Microbial Genome Atlas provides a genomic and functional resource to enable future studies on salmon nutrition and health.

Genome Survey of Male Rana dybowskii to Further Understand the Sex Determination Mechanism.

Rana dybowskii is one of the important aquaculture species in Northeast China. The fallopian tubes of female R. dybowskii are used to prepare oviductus ranae (an important traditional Chinese medicine). Therefore, R. dybowskii females have higher economical value than males. An increasing female R. dybowskii population can increase the benefits from R. dybowskii culture. However, the genome of amphibians is complex, making it difficult to investigate their sex determination mechanism. In this study, we analyzed the genome of male R. dybowskii using next-generation sequencing technology. A total of 200,046,452,400 bp of clean data were obtained, and the K-mer analysis indicated that the depth was 50×. The genome size of R. dybowskii was approximately 3585.05 M, with a heterozygosity rate, repeat sequence ratio, and genome GC content of 1.15%, 68.96%, and approximately 43.0%, respectively. In total, 270,785 contigs and 498 scaffolds were generated. The size of the contigs and scaffolds was 3,748,543,415 and 3,765,862,278 bp, respectively, with the N50 length of 31,988 and 336,385,783. The longest contig and scaffold were of the size 137,967,485 and 1,808,367,828 bp, respectively. The number of contigs and scaffolds > 10K nt was 99,620 and 451, respectively. Through annotation, 40,913 genes were obtained, including 156,609 CDS (i.e., 3.83 CDS per gene). Sequence alignment was performed with the assembled scaffolding genome in this study. Two and one fragment had high homology with two male-specific DNA molecular markers of R. dybowskii discovered previously (namely, MSM-222 and MSM-261, respectively). In addition, the Dmrt1 gene of R. dybowskii was obtained with a length of 18,893 bp by comparison and splicing. The forward primers amplifying MSM-222 and MSM-261 were located at 322-343 and 14,501-14,526 bp of Dmrt1, respectively. However, sequence alignment revealed that MSM-222 and MSM-261 were not located on Dmrt1, and only some homologous parts were observed. This indicated that in addition to Dmrt1, other important genes may play a crucial role in the sex determination mechanism of R. dybowskii. Our study provided a foundation for the subsequent high-quality genome construction and provided important genomic resources for future studies on R. dybowskii.

Chromosome-level genome of diamondback terrapin provides insight into the genetic basis of salinity adaptation.

Diamondback terrapins (Malaclemys terrapin centrata) exhibit strong environmental adaptability and live in both freshwater and saltwater. However, the genetic basis of this adaptability has not been the focus of research. In this study, we successfully constructed a ∼2.21-Gb chromosome-level genome assembly for M. t. centrata using high-coverage and high-depth genomic sequencing data generated on multiple platforms. The M. t. centrata genome contains 25 chromosomes and the scaffold N50 of ∼143.75 Mb, demonstrating high continuity and accuracy. In total, 53.82% of the genome assembly was composed of repetitive sequences, and 22435 protein-coding genes were predicted. Our phylogenetic analysis indicated that M. t. centrata was closely related to the red-eared slider turtle (Trachemys scripta elegans), with divergence approximately ∼23.6 million years ago (Mya) during the early Neogene period of the Cenozoic era. The population size of M. t. centrata decreased significantly over the past ∼14 Mya during the Cenozoic era. Comparative genomic analysis indicated that 36 gene families related to ion transport were expanded and several genes (AQP3, solute carrier subfamily, and potassium channel genes) underwent specific amino acid site mutations in the M. t. centrata genome. Changes to these ion transport-related genes may have contributed to the remarkable salinity adaptability of diamondback terrapin. The results of this study not only provide a high-quality reference genome for M. t. centrata but also elucidate the possible genetic basis for salinity adaptation in this species.