Generalized reporter score-based enrichment analysis for omics data.

Enrichment analysis contextualizes biological features in pathways to facilitate a systematic understanding of high-dimensional data and is widely used in biomedical research. The emerging reporter score-based analysis (RSA) method shows more promising sensitivity, as it relies on P-values instead of raw values of features. However, RSA cannot be directly applied to multi-group and longitudinal experimental designs and is often misused due to the lack of a proper tool. Here, we propose the Generalized Reporter Score-based Analysis (GRSA) method for multi-group and longitudinal omics data. A comparison with other popular enrichment analysis methods demonstrated that GRSA had increased sensitivity across multiple benchmark datasets. We applied GRSA to microbiome, transcriptome and metabolome data and discovered new biological insights in omics studies. Finally, we demonstrated the application of GRSA beyond functional enrichment using a taxonomy database. We implemented GRSA in an R package, ReporterScore, integrating with a powerful visualization module and updatable pathway databases, which is available on the Comprehensive R Archive Network (https://cran.r-project.org/web/packages/ReporterScore). We believe that the ReporterScore package will be a valuable asset for broad biomedical research fields.

Validation of key sponge symbiont pathways using genome-centric metatranscriptomics.

The sponge microbiome underpins host function through provision and recycling of essential nutrients in a nutrient poor environment. Genomic data suggest that carbohydrate degradation, carbon fixation, nitrogen metabolism, sulphur metabolism and supplementation of B-vitamins are central microbial functions. However, validation beyond the genomic potential of sponge symbiont pathways is rarely explored. To evaluate metagenomic predictions, we sequenced the metagenomes and metatranscriptomes of three common coral reef sponges: Ircinia ramosa, Ircinia microconulosa and Phyllospongia foliascens. Multiple carbohydrate active enzymes were expressed by Poribacteria, Bacteroidota and Cyanobacteria symbionts, suggesting these lineages have a central role in assimilating dissolved organic matter. Expression of entire pathways for carbon fixation and multiple sulphur compound transformations were observed in all sponges. Gene expression for anaerobic nitrogen metabolism (denitrification and nitrate reduction) were more common than aerobic metabolism (nitrification), where only the I. ramosa microbiome expressed the nitrification pathway. Finally, while expression of the biosynthetic pathways for B-vitamins was common, the expression of additional transporter genes was far more limited. Overall, we highlight consistencies and disparities between metagenomic and metatranscriptomic results when inferring microbial activity, while uncovering new microbial taxa that contribute to the health of their sponge host via nutrient exchange.

The complete and fully-phased diploid genome of a male Han Chinese.

Since the release of the complete human genome, the priority of human genomic study has now been shifting towards closing gaps in ethnic diversity. Here, we present a fully phased and well-annotated diploid human genome from a Han Chinese male individual (CN1), in which the assemblies of both haploids achieve the telomere-to-telomere (T2T) level. Comparison of this diploid genome with the CHM13 haploid T2T genome revealed significant variations in the centromere. Outside the centromere, we discovered 11,413 structural variations, including numerous novel ones. We also detected thousands of CN1 alleles that have accumulated high substitution rates and a few that have been under positive selection in the East Asian population. Further, we found that CN1 outperforms CHM13 as a reference genome in mapping and variant calling for the East Asian population owing to the distinct structural variants of the two references. Comparison of SNP calling for a large cohort of 8869 Chinese genomes using CN1 and CHM13 as reference respectively showed that the reference bias profoundly impacts rare SNP calling, with nearly 2 million rare SNPs miss-called with different reference genomes. Finally, applying the CN1 as a reference, we discovered 5.80 Mb and 4.21 Mb putative introgression sequences from Neanderthal and Denisovan, respectively, including many East Asian specific ones undetected using CHM13 as the reference. Our analyses reveal the advances of using CN1 as a reference for population genomic studies and paleo-genomic studies. This complete genome will serve as an alternative reference for future genomic studies on the East Asian population.

Testing cophylogeny between coral reef invertebrates and their bacterial and archaeal symbionts.

Marine invertebrates harbour a complex suite of bacterial and archaeal symbionts, a subset of which are probably linked to host health and homeostasis. Within a complex microbiome it can be difficult to tease apart beneficial or parasitic symbionts from nonessential commensal or transient microorganisms; however, one approach is to detect strong cophylogenetic patterns between microbial lineages and their respective hosts. We employed the Procrustean approach to cophylogeny (PACo) on 16S rRNA gene derived microbial community profiles paired with COI, 18S rRNA and ITS1 host phylogenies. Second, we undertook a network analysis to identify groups of microbes that were co-occurring within our host species. Across 12 coral, 10 octocoral and five sponge species, each host group and their core microbiota (50% prevalence within host species replicates) had a significant fit to the cophylogenetic model. Independent assessment of each microbial genus and family found that bacteria and archaea affiliated to Endozoicomonadaceae, Spirochaetaceae and Nitrosopumilaceae have the strongest cophylogenetic signals. Further, local Moran's I measure of spatial autocorrelation identified 14 ASVs, including Endozoicomonadaceae and Spirochaetaceae, whose distributions were significantly clustered by host phylogeny. Four co-occurring subnetworks were identified, each of which was dominant in a different host group. Endozoicomonadaceae and Spirochaetaceae ASVs were abundant among the subnetworks, particularly one subnetwork that was exclusively comprised of these two bacterial families and dominated the octocoral microbiota. Our results disentangle key microbial interactions that occur within complex microbiomes and reveal long-standing, essential microbial symbioses in coral reef invertebrates.

Evolutionary and biomedical insights from a marmoset diploid genome assembly.

The accurate and complete assembly of both haplotype sequences of a diploid organism is essential to understanding the role of variation in genome functions, phenotypes and diseases1. Here, using a trio-binning approach, we present a high-quality, diploid reference genome, with both haplotypes assembled independently at the chromosome level, for the common marmoset (Callithrix jacchus), an primate model system that is widely used in biomedical research2,3. The full spectrum of heterozygosity between the two haplotypes involves 1.36% of the genome-much higher than the 0.13% indicated by the standard estimation based on single-nucleotide heterozygosity alone. The de novo mutation rate is 0.43 × 10-8 per site per generation, and the paternal inherited genome acquired twice as many mutations as the maternal. Our diploid assembly enabled us to discover a recent expansion of the sex-differentiation region and unique evolutionary changes in the marmoset Y chromosome. In addition, we identified many genes with signatures of positive selection that might have contributed to the evolution of Callithrix biological features. Brain-related genes were highly conserved between marmosets and humans, although several genes experienced lineage-specific copy number variations or diversifying selection, with implications for the use of marmosets as a model system.

Efficient COI barcoding using high throughput single-end 400 bp sequencing.

BACKGROUND: Over the last decade, the rapid development of high-throughput sequencing platforms has accelerated species description and assisted morphological classification through DNA barcoding. However, the current high-throughput DNA barcoding methods cannot obtain full-length barcode sequences due to read length limitations (e.g. a maximum read length of 300 bp for the Illumina's MiSeq system), or are hindered by a relatively high cost or low sequencing output (e.g. a maximum number of eight million reads per cell for the PacBio's SEQUEL II system). RESULTS: Pooled cytochrome c oxidase subunit I (COI) barcodes from individual specimens were sequenced on the MGISEQ-2000 platform using the single-end 400 bp (SE400) module. We present a bioinformatic pipeline, HIFI-SE, that takes reads generated from the 5' and 3' ends of the COI barcode region and assembles them into full-length barcodes. HIFI-SE is written in Python and includes four function modules of filter, assign, assembly and taxonomy. We applied the HIFI-SE to a set of 845 samples (30 marine invertebrates, 815 insects) and delivered a total of 747 fully assembled COI barcodes as well as 70 Wolbachia and fungi symbionts. Compared to their corresponding Sanger sequences (72 sequences available), nearly all samples (71/72) were correctly and accurately assembled, including 46 samples that had a similarity score of 100% and 25 of ca. 99%. CONCLUSIONS: The HIFI-SE pipeline represents an efficient way to produce standard full-length barcodes, while the reasonable cost and high sensitivity of our method can contribute considerably more DNA barcodes under the same budget. Our method thereby advances DNA-based species identification from diverse ecosystems and increases the number of relevant applications.

Diverse coral reef invertebrates exhibit patterns of phylosymbiosis.

Microbiome assemblages of plants and animals often show a degree of correlation with host phylogeny; an eco-evolutionary pattern known as phylosymbiosis. Using 16S rRNA gene sequencing to profile the microbiome, paired with COI, 18S rRNA and ITS1 host phylogenies, phylosymbiosis was investigated in four groups of coral reef invertebrates (scleractinian corals, octocorals, sponges and ascidians). We tested three commonly used metrics to evaluate the extent of phylosymbiosis: (a) intraspecific versus interspecific microbiome variation, (b) topological comparisons between host phylogeny and hierarchical clustering (dendrogram) of host-associated microbial communities, and (c) correlation of host phylogenetic distance with microbial community dissimilarity. In all instances, intraspecific variation in microbiome composition was significantly lower than interspecific variation. Similarly, topological congruency between host phylogeny and the associated microbial dendrogram was more significant than would be expected by chance across all groups, except when using unweighted UniFrac distance (compared with weighted UniFrac and Bray-Curtis dissimilarity). Interestingly, all but the ascidians showed a significant positive correlation between host phylogenetic distance and associated microbial dissimilarity. Our findings provide new perspectives on the diverse nature of marine phylosymbioses and the complex roles of the microbiome in the evolution of marine invertebrates.

Metatranscriptomics of the bacterial community in response to atmospheric deposition in the Western North Pacific Ocean.

Atmospheric deposition represents a major vector of both macro- and micro-nutrients to the oligotrophic open oceans, potentially imposing a profound impact on the functioning of the microbial community. Whereas bacterial responses to atmospheric deposition are being studied at the community level, corresponding functional changes are essentially unknown. Here we conducted a microcosm experiment coupled with metatranscriptomic analyses to elucidate taxonomic and functional profiles of the bacterial community in response to East Asian aerosols in the Western North Pacific Ocean (WNPO). While the abundance of heterotrophic bacteria showed a minor change, cyanobacterial cell count number decreased dramatically, with Prochlorococcus and Synechococcus counts reduced by 83.2% and 21.5% in the aerosol treatment in relation to the control. Expression of transcripts related with Prochlorococcus, Synechococcus, Trichodesmium and Crocosphaera both were lower in the treatment (5.7%, 2.3%, 0.5% and 0.02%, respectively) than in the control (18.6%, 2.7%, 9.8% and 0.14%, respectively). Aerosol addition led to an increase in transcripts involved in iron metabolism (tonB, feoB, irr, exbB), indicating Fe limitation. Heavy metal toxicity was evidenced by an elevated expression of resistance genes, such as czcC, czcB, czcA and a probable Co/Zn/Cd efflux protein, and a range of genes functioning against oxidative stress. Our findings provide insights into an inhibitory effect of high-flux East Asian aerosols on cyanobacteria in the WNPO likely due to Fe scavenging and heavy metal toxicity.

Growth and Toxin Production of Gambierdiscus spp. Can Be Regulated by Quorum-Sensing Bacteria.

Gambierdiscus spp. are the major culprit responsible for global ciguatera fish poisoning (CFP). At present, the effects of microbiological factors on algal proliferation and toxin production are poorly understood. To evaluate the regulatory roles of quorum-sensing (QS) bacteria in the physiology of Gambierdiscus, co-culture experiments with screened QS strains were conducted in this study. Except for the growth-inhibiting effect from the strain Marinobacter hydrocarbonoclasticus, the algal host generally displayed much higher growth potential and toxin production ability with the existence of QS strains. In addition, Bacillus anthracis particularly exhibited a broad-spectrum growth enhancement effect on various Gambierdiscus types, as well as a remarkable influence on algal toxicity. The variations of algal physiological status, including growth rate, chlorophyll content, and responsive behaviors, are potential reasons for the observed positive or negative affection. This study suggests that QS bacteria regulate the algal growth and toxin production. Based on the evidence, we further speculate that QS bacteria may contribute to the site-specific distribution of CFP risk through regulating the algal host biomass and toxicity.

Unravel the hidden protistan diversity: application of blocking primers to suppress PCR amplification of metazoan DNA.

Planktonic protists, including both autotroph and heterotroph, have been recognized as a major contributor to primary production and consumers of bacteria, archaea and picophytoplankton. However, the understanding of protistan diversity is typically impeded by the large amount of metazoans when employing universal primers to environmental samples. In this study, we developed 3 blocking primers to inhibit the amplification of metazoan DNA in PCR. First, we optimized the design and concentration of blocking primers by using 2 metazoans and 3 protists in quantitative PCR. No significant difference (q > 0.05) was found in protistan community structure at phylum, family or OTU levels and genetic diversity between samples amplified with and without adding blocking primers, indicating that the blocking primer does not alter the composition of protistan community. An application to samples with a high abundance of metazoans showed that the blocking primer can reduce 42.1-72.4% of metazoan sequences, resulting the retrieval of a higher protistan richness and diversity. The use of the blocking primer helped to identify protistan community composition from a large size fraction, which is usually dominated by metazoans. Protistan community from the small (0.2-10 µm) and large (> 10 µm) size fractions exhibited a low similarity of 36.6% and shared 14.7% OTUs. About 63.8% of the OTUs were unique to the large fraction. Species from groups, such as Lepidodinium, Warnowia, Kareniaceae, Torodiniales, Phaeocystis, Chrysochromulina and Chrysophyceae, were enriched in the large fraction, indicating that they could be largely underestimated in studies that exclude the large-sized cells. Blocking primers are a promising tool to increase the efficiency for the characterization of protistan diversity in aquatic planktonic ecosystems.

Synechococcus Assemblages across the Salinity Gradient in a Salt Wedge Estuary.

Synechococcus are the most abundant and widely distributed picocyanobacteria in the ocean. The salt-wedge type of estuary possesses the complete horizontal and vertical gradient of salinity together with other physical and chemical parameters. In order to reveal whether such a complex environmental gradient harbors a high diversity of Synechococcus, we investigated the abundance, taxonomic composition and pigment genetic diversity of Synechococcus in surface and bottom waters across the salinity gradient in a salt-wedge estuary by flow cytometric analysis and pyrosequencing of the rpoC1 gene and cpcBA operon (encoding phycocyanin). Synechococcus were ubiquitously distributed in the studied region, with clear spatial variations both horizontally and vertically. The abundance and diversity of Synechococcus were low in the freshwater-dominated low salinity waters. By pyrosequencing of the rpoC1 gene, we have shown that with the increase of salinity, the dominant Synechococcus shifted from the freshwater Synechococcus to the combination of phylogenetic subcluster 5.2 and freshwater Synechococcus, and then the strictly marine subcluster 5.1 clade III. Besides, the composition of Synechococcus assemblage in the deep layer was markedly different from the surface in the stratified waters (dissimilarities: 40.32%-95.97%, SIMPER analysis). High abundance of clade III Synechococcus found in the brackish waters may revise our previous understanding that strains of this clade prefers oligotrophic environment. Our data also suggested that both the phylogenetic subcluster 5.3 Synechococcus, a lineage that was not well understood, and subcluster 5.1 clade I, a typical cold water lineage, were widely distributed in the bottom layer of the estuary. Clade I detected in the studied region was mainly contributed by subclade IG. Analysis of the cpcBA operon sequences revealed niche partitioning between type 1 and type 3 Synechococcus, with type 2 distributed broadly across the whole environmental gradients. Our results suggest that the salt wedge estuary provides various niches for different lineages of Synechococcus, making it an environment with high Synechococcus diversity compared with adjacent freshwater and shelf sea environments.

Characterization and expression of MHC class II alpha and II beta genes in mangrove red snapper (Lutjanus argentimaculatus).

The major histocompatibility complex (MHC) class II plays a key role in adaptive immunity by presenting foreign peptides to CD4(+) T cells and by triggering the adaptive immune response. While the structure and function of MHC class II have been well characterized in mammalian, limited research has been done on fishes. In this study, we characterized the gene structure and expression of MHC class II α (Lunar-DAA) and II ß (Lunar-DAB) of mangrove red snapper (Lutjanus argentimaculatus). Both genes shared, respectively, a high similarity and typical features with other vertebrate MHC class II α and II ß. The phylogenetic analysis of the deduced peptides revealed that both Lunar-DAA and Lunar-DAB were located in the teleost subclass. Western blotting analyses indicated that both MHC class II α and II ß were expressed ubiquitously in immune-related cells, tissues and organs, and that MHC class II α and II ß chains existed mainly as heterodimers. While it was highly expressed in gills, thymus, head kidney (HK), spleen, head kidney macrophage and spleen leucocytes, MHC class II ß chain was expressed with a low abundance in skin, intestine, stomach and heart. The highest expression of MHC class II ß in thymus confirmed the conclusion that thymus is one of the primary lymphoid organs in fishes. The detection of MHC class II αß dimers in HK macrophages and spleen leucocytes indicated that HK macrophages and spleen leucocytes play a critical role in the adaptive immunity in fishes. All these results provide valuable information for understanding the structure of MHC class II α and II ß and their function in immune responses.

An association network analysis among microeukaryotes and bacterioplankton reveals algal bloom dynamics.

Algal blooms are a worldwide phenomenon and the biological interactions that underlie their regulation are only just beginning to be understood. It is established that algal microorganisms associate with many other ubiquitous, oceanic organisms, but the interactions that lead to the dynamics of bloom formation are currently unknown. To address this gap, we used network approaches to investigate the association patterns among microeukaryotes and bacterioplankton in response to a natural Scrippsiella trochoidea bloom. This is the first study to apply network approaches to bloom dynamics. To this end, terminal restriction fragment (T-RF) length polymorphism analysis showed dramatic changes in community compositions of microeukaryotes and bacterioplankton over the blooming period. A variance ratio test revealed significant positive overall associations both within and between microeukaryotic and bacterioplankton communities. An association network generated from significant correlations between T-RFs revealed that S. trochoidea had few connections to other microeukaryotes and bacterioplankton and was placed on the edge. This lack of connectivity allowed for the S. trochoidea sub-network to break off from the overall network. These results allowed us to propose a conceptual model for explaining how changes in microbial associations regulate the dynamics of an algal bloom. In addition, key T-RFs were screened by principal components analysis, correlation coefficients, and network analysis. Dominant T-RFs were then identified through 18S and 16S rRNA gene clone libraries. Results showed that microeukaryotes clustered predominantly with Dinophyceae and Perkinsea while the majority of bacterioplankton identified were Alphaproteobacteria, Gammaproteobacteria, and Bacteroidetes. The ecologi-cal roles of both were discussed in the context of these findings.

Identification and functional characterization of a putative 17ß-hydroxysteroid dehydrogenase 12 in abalone (Haliotis diversicolor supertexta).

The 17ß-hydroxysteroid dehydrogenases (17ß-HSDs) are key enzymes in the downstream process of steroid hormone biosynthesis. To date, relatively little is known about the role of 17ß-HSDs in marine gastropods. In the present study, a putative cDNA sequence encoding type 12 17ß-HSD (17ß-HSD-12) was identified in abalone (Haliotis diversicolor supertexta). The full-length cDNA was 1,978 bp, including an open reading frame (ORF) of 963 bp that encoded a protein of 321 amino acids. Comparative structural analysis revealed that abalone 17ß-HSD-12 shared 39.8-42.8% amino acid identity with other 17ß-HSD-12 homologues and that the functional domains were well conserved. Phylogenetic analysis revealed that abalone 17ß-HSD-12 belonged to the short-chain dehydrogenases/reductases (SDRs) family. Functional analysis following transient transfection of the ORF in human embryonic kidney-293 (HEK-293) cells indicated that abalone 17ß-HSD-12 had the ability to convert estrone (E1) into estradiol (E2). Expression analysis in vivo demonstrated that abalone 17ß-HSD-12 was differentially expressed during the three reproductive stages (pre-spawning, spawning, and post-spawning). These results indicate that abalone 17ß-HSD-12 is an SDR family member with a key role in steroidogenesis during the reproductive period.