MetaSquare: an integrated metadatabase of 16S rRNA gene amplicon for microbiome taxonomic classification.

MOTIVATION: Taxonomic classification of 16S ribosomal RNA gene amplicon is an efficient and economic approach in microbiome analysis. 16S rRNA sequence databases like SILVA, RDP, EzBioCloud and HOMD used in downstream bioinformatic pipelines have limitations on either the sequence redundancy or the delay on new sequence recruitment. To improve the 16S rRNA gene-based taxonomic classification, we merged these widely used databases and a collection of novel sequences systemically into an integrated resource. RESULTS: MetaSquare version 1.0 is an integrated 16S rRNA sequence database. It is composed of more than 6 million sequences and improves taxonomic classification resolution on both long-read and short-read methods. AVAILABILITY AND IMPLEMENTATION: Accessible at https://hub.docker.com/r/lsbnb/metasquare_db and https://github.com/lsbnb/MetaSquare. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Intelligent De Novo Design of Novel Antimicrobial Peptides against Antibiotic-Resistant Bacteria Strains.

Because of the growing number of clinical antibiotic resistance cases in recent years, novel antimicrobial peptides (AMPs) may be ideal for next-generation antibiotics. This study trained a Wasserstein generative adversarial network with gradient penalty (WGAN-GP) based on known AMPs to generate novel AMP candidates. The quality of the GAN-designed peptides was evaluated in silico, and eight of them, named GAN-pep 1-8, were selected by an AMP Artificial Intelligence (AI) classifier and synthesized for further experiments. Disc diffusion testing and minimum inhibitory concentration (MIC) determinations were used to identify the antibacterial effects of the synthesized GAN-designed peptides. Seven of the eight synthesized GAN-designed peptides displayed antibacterial activity. Additionally, GAN-pep 3 and GAN-pep 8 presented a broad spectrum of antibacterial effects and were effective against antibiotic-resistant bacteria strains, such as methicillin-resistant Staphylococcus aureus and carbapenem-resistant Pseudomonas aeruginosa. GAN-pep 3, the most promising GAN-designed peptide candidate, had low MICs against all the tested bacteria. In brief, our approach shows an efficient way to discover AMPs effective against general and antibiotic-resistant bacteria strains. In addition, such a strategy also allows other novel functional peptides to be quickly designed, identified, and synthesized for validation on the wet bench.

Enrichment of Prevotella intermedia in human colorectal cancer and its additive effects with Fusobacterium nucleatum on the malignant transformation of colorectal adenomas.

BACKGROUND: Owing to the heterogeneity of microbiota among individuals and populations, only Fusobacterium nucleatum and Bacteroides fragilis have been reported to be enriched in colorectal cancer (CRC) in multiple studies. Thus, the discovery of additional bacteria contributing to CRC development in various populations can be expected. We aimed to identify bacteria associated with the progression of colorectal adenoma to carcinoma and determine the contribution of these bacteria to malignant transformation in patients of Han Chinese origin. METHODS: Microbiota composition was determined through 16S rRNA V3-V4 amplicon sequencing of autologous adenocarcinomas, adenomatous polyps, and non-neoplastic colon tissue samples (referred to as "tri-part samples") in patients with CRC. Enriched taxa in adenocarcinoma tissues were identified through pairwise comparison. The abundance of candidate bacteria was quantified through genomic quantitative polymerase chain reaction (qPCR) in tissue samples from 116 patients. Associations of candidate bacteria with clinicopathological features and genomic and genetic alterations were evaluated through odds ratio tests. Additionally, the effects of candidate bacteria on CRC cell proliferation, migration, and invasion were evaluated through the co-culture of CRC cells with bacterial cells or with conditioned media from bacteria. RESULTS: Prevotella intermedia was overrepresented in adenocarcinomas compared with paired adenomatous polyps. Furthermore, co-abundance of P. intermedia and F. nucleatum was observed in tumor tissues. More notably, the coexistence of these two bacteria in adenocarcinomas was associated with lymph node involvement and distant metastasis. These two bacteria also exerted additive effects on the enhancement of the migration and invasion abilities of CRC cells. Finally, conditioned media from P. intermedia promoted the migration and invasion of CRC cells. CONCLUSION: This report is the first to demonstrate that P. intermedia is enriched in colorectal adenocarcinoma tissues and enhances the migration and invasion abilities of CRC cells. Moreover, P. intermedia and F. nucleatum exert additive effects on the malignant transformation of colorectal adenomas into carcinomas. These findings can be used to identify patients at a high risk of malignant transformation of colorectal adenomas or metastasis of CRC, and they can accordingly be provided optimal clinical management.

EpiMOLAS: an intuitive web-based framework for genome-wide DNA methylation analysis.

BACKGROUND: DNA methylation is a crucial epigenomic mechanism in various biological processes. Using whole-genome bisulfite sequencing (WGBS) technology, methylated cytosine sites can be revealed at the single nucleotide level. However, the WGBS data analysis process is usually complicated and challenging. RESULTS: To alleviate the associated difficulties, we integrated the WGBS data processing steps and downstream analysis into a two-phase approach. First, we set up the required tools in Galaxy and developed workflows to calculate the methylation level from raw WGBS data and generate a methylation status summary, the mtable. This computation environment is wrapped into the Docker container image DocMethyl, which allows users to rapidly deploy an executable environment without tedious software installation and library dependency problems. Next, the mtable files were uploaded to the web server EpiMOLAS_web to link with the gene annotation databases that enable rapid data retrieval and analyses. CONCLUSION: To our knowledge, the EpiMOLAS framework, consisting of DocMethyl and EpiMOLAS_web, is the first approach to include containerization technology and a web-based system for WGBS data analysis from raw data processing to downstream analysis. EpiMOLAS will help users cope with their WGBS data and also conduct reproducible analyses of publicly available data, thereby gaining insights into the mechanisms underlying complex biological phenomenon. The Galaxy Docker image DocMethyl is available at https://hub.docker.com/r/lsbnb/docmethyl/. EpiMOLAS_web is publicly accessible at http://symbiosis.iis.sinica.edu.tw/epimolas/.

The role of the bacterial protease Prc in the uropathogenesis of extraintestinal pathogenic Escherichia coli.

BACKGROUND: Extraintestinal pathogenic E. coli (ExPEC) remains one of the most prevalent bacterial pathogens that cause extraintestinal infections, including neonatal meningitis, septicemia, and urinary tract (UT) infections (UTIs). Antibiotic therapy has been the conventional treatment for such infections, but its efficacy has decreased due to the emergence of antibiotic-resistant bacteria. Identification and characterization of bacterial factors that contribute to the severity of infection would facilitate the development of novel therapeutic strategies. The ExPEC periplasmic protease Prc contributes to the pathogen's ability to evade complement-mediated killing in the serum. Here, we further investigated the role of the Prc protease in ExPEC-induced UTIs and the underlying mechanism. METHODS: The uropathogenic role of Prc was determined in a mouse model of UTIs. Using global quantitative proteomic analyses, we revealed that the expression of FliC and other outer membrane-associated proteins was altered by Prc deficiency. Comparative transcriptome analyses identified that Prc deficiency affected expression of the flagellar regulon and genes that are regulated by five extracytoplasmic signaling systems. RESULTS: A mutant ExPEC with a prc deletion was attenuated in bladder and kidney colonization. Global quantitative proteomic analyses of the prc mutant and wild-type ExPEC strains revealed significantly reduced flagellum expression in the absence of Prc, consequently impairing bacterial motility. The prc deletion triggered downregulation of the flhDC operon encoding the master transcriptional regulator of flagellum biogenesis. Overexpressing flhDC restored the prc mutant's motility and ability to colonize the UT, suggesting that the impaired motility is responsible for attenuated UT colonization of the mutant. Further comparative transcriptome analyses revealed that Prc deficiency activated the σE and RcsCDB signaling pathways. These pathways were responsible for the diminished flhDC expression. Finally, the activation of the RcsCDB system was attributed to the intracellular accumulation of a known Prc substrate Spr in the prc mutant. Spr is a peptidoglycan hydrolase and its accumulation destabilizes the bacterial envelope. CONCLUSIONS: We demonstrated for the first time that Prc is essential for full ExPEC virulence in UTIs. Our results collectively support the idea that Prc is essential for bacterial envelope integrity, thus explaining how Prc deficiency results in an attenuated ExPEC.

Systematic identification of anti-interferon function on hepatitis C virus genome reveals p7 as an immune evasion protein.

Hepatitis C virus (HCV) encodes mechanisms to evade the multilayered antiviral actions of the host immune system. Great progress has been made in elucidating the strategies HCV employs to down-regulate interferon (IFN) production, impede IFN signaling transduction, and impair IFN-stimulated gene (ISG) expression. However, there is a limited understanding of the mechanisms governing how viral proteins counteract the antiviral functions of downstream IFN effectors due to the lack of an efficient approach to identify such interactions systematically. To study the mechanisms by which HCV antagonizes the IFN responses, we have developed a high-throughput profiling platform that enables mapping of HCV sequences critical for anti-IFN function at high resolution. Genome-wide profiling performed with a 15-nt insertion mutant library of HCV showed that mutations in the p7 region conferred high levels of IFN sensitivity, which could be alleviated by the expression of WT p7 protein. This finding suggests that p7 protein of HCV has an immune evasion function. By screening a liver-specific ISG library, we identified that IFI6-16 significantly inhibits the replication of p7 mutant viruses without affecting WT virus replication. In contrast, knockout of IFI6-16 reversed the IFN hypersensitivity of p7 mutant virus. In addition, p7 was found to be coimmunoprecipitated with IFI6-16 and to counteract the function of IFI6-16 by depolarizing the mitochondria potential. Our data suggest that p7 is a critical immune evasion protein that suppresses the antiviral IFN function by counteracting the function of IFI6-16.

SQUAT: a Sequencing Quality Assessment Tool for data quality assessments of genome assemblies.

BACKGROUND: With the rapid increase in genome sequencing projects for non-model organisms, numerous genome assemblies are currently in progress or available as drafts, but not made available as satisfactory, usable genomes. Data quality assessment of genome assemblies is gaining importance not only for people who perform the assembly/re-assembly processes, but also for those who attempt to use assemblies as maps in downstream analyses. Recent studies of the quality control, quality evaluation/ assessment of genome assemblies have focused on either quality control of reads before assemblies or evaluation of the assemblies with respect to their contiguity and correctness. However, correctness assessment depends on a reference and is not applicable for de novo assembly projects. Hence, development of methods providing both post-assembly and pre-assembly quality assessment reports for examining the quality/correctness of de novo assemblies and the input reads is worth studying. RESULTS: We present SQUAT, an efficient tool for both pre-assembly and post-assembly quality assessment of de novo genome assemblies. The pre-assembly module of SQUAT computes quality statistics of reads and presents the analysis in a well-designed interface to visualize the distribution of high- and poor-quality reads in a portable HTML report. The post-assembly module of SQUAT provides read mapping analytics in an HTML format. We categorized reads into several groups including uniquely mapped reads, multiply mapped, unmapped reads; for uniquely mapped reads, we further categorized them into perfectly matched, with substitutions, containing clips, and the others. We carefully defined the poorly mapped (PM) reads into several groups to prevent the underestimation of unmapped reads; indeed, a high PM% would be a sign of a poor assembly that requires researchers' attention for further examination or improvements before using the assembly. Finally, we evaluate SQUAT with six datasets, including the genome assemblies for eel, worm, mushroom, and three bacteria. The results show that SQUAT reports provide useful information with details for assessing the quality of assemblies and reads. AVAILABILITY: The SQUAT software with links to both its docker image and the on-line manual is freely available at https://github.com/luke831215/SQUAT .

A gene profiling deconvolution approach to estimating immune cell composition from complex tissues.

BACKGROUND: A new emerged cancer treatment utilizes intrinsic immune surveillance mechanism that is silenced by those malicious cells. Hence, studies of tumor infiltrating lymphocyte populations (TILs) are key to the success of advanced treatments. In addition to laboratory methods such as immunohistochemistry and flow cytometry, in silico gene expression deconvolution methods are available for analyses of relative proportions of immune cell types. RESULTS: Herein, we used microarray data from the public domain to profile gene expression pattern of twenty-two immune cell types. Initially, outliers were detected based on the consistency of gene profiling clustering results and the original cell phenotype notation. Subsequently, we filtered out genes that are expressed in non-hematopoietic normal tissues and cancer cells. For every pair of immune cell types, we ran t-tests for each gene, and defined differentially expressed genes (DEGs) from this comparison. Equal numbers of DEGs were then collected as candidate lists and numbers of conditions and minimal values for building signature matrixes were calculated. Finally, we used v -Support Vector Regression to construct a deconvolution model. The performance of our system was finally evaluated using blood biopsies from 20 adults, in which 9 immune cell types were identified using flow cytometry. The present computations performed better than current state-of-the-art deconvolution methods. CONCLUSIONS: Finally, we implemented the proposed method into R and tested extensibility and usability on Windows, MacOS, and Linux operating systems. The method, MySort, is wrapped as the Galaxy platform pluggable tool and usage details are available at https://testtoolshed.g2.bx.psu.edu/view/moneycat/mysort/e3afe097e80a .

TEA: the epigenome platform for Arabidopsis methylome study.

BACKGROUND: Bisulfite sequencing (BS-seq) has become a standard technology to profile genome-wide DNA methylation at single-base resolution. It allows researchers to conduct genome-wise cytosine methylation analyses on issues about genomic imprinting, transcriptional regulation, cellular development and differentiation. One single data from a BS-Seq experiment is resolved into many features according to the sequence contexts, making methylome data analysis and data visualization a complex task. RESULTS: We developed a streamlined platform, TEA, for analyzing and visualizing data from whole-genome BS-Seq (WGBS) experiments conducted in the model plant Arabidopsis thaliana. To capture the essence of the genome methylation level and to meet the efficiency for running online, we introduce a straightforward method for measuring genome methylation in each sequence context by gene. The method is scripted in Java to process BS-Seq mapping results. Through a simple data uploading process, the TEA server deploys a web-based platform for deep analysis by linking data to an updated Arabidopsis annotation database and toolkits. CONCLUSIONS: TEA is an intuitive and efficient online platform for analyzing the Arabidopsis genomic DNA methylation landscape. It provides several ways to help users exploit WGBS data. TEA is freely accessible for academic users at: http://tea.iis.sinica.edu.tw .

Characterization of FN1-FGFR1 and novel FN1-FGF1 fusion genes in a large series of phosphaturic mesenchymal tumors.

Phosphaturic mesenchymal tumors typically cause paraneoplastic osteomalacia, chiefly as a result of FGF23 secretion. In a prior study, we identified FN1-FGFR1 fusion in 9 of 15 phosphaturic mesenchymal tumors. In this study, a total of 66 phosphaturic mesenchymal tumors and 7 tumors resembling phosphaturic mesenchymal tumor but without known phosphaturia were studied. A novel FN1-FGF1 fusion gene was identified in two cases without FN1-FGFR1 fusion by RNA sequencing and cross-validated with direct sequencing and western blot. Fluorescence in situ hybridization analyses revealed FN1-FGFR1 fusion in 16 of 39 (41%) phosphaturic mesenchymal tumors and identified an additional case with FN1-FGF1 fusion. The two fusion genes were mutually exclusive. Combined with previous data, the overall prevalence of FN1-FGFR1 and FN1-FGF1 fusions was 42% (21/50) and 6% (3/50), respectively. FGFR1 immunohistochemistry was positive in 82% (45/55) of phosphaturic mesenchymal tumors regardless of fusion status. By contrast, 121 cases of potential morphologic mimics (belonging to 13 tumor types) rarely expressed FGFR1, the main exceptions being solitary fibrous tumors (positive in 40%), chondroblastomas (40%), and giant cell tumors of bone (38%), suggesting a possible role for FGFR1 immunohistochemistry in the diagnosis of phosphaturic mesenchymal tumor. With the exception of one case reported in our prior study, none of the remaining tumors resembling phosphaturic mesenchymal tumor had either fusion type or expressed significant FGFR1. Our findings provide insight into possible mechanisms underlying the pathogenesis of phosphaturic mesenchymal tumor and imply a central role of the FGF1-FGFR1 signaling pathway. The novel FN1-FGF1 protein is expected to be secreted and serves as a ligand that binds and activates FGFR1 to achieve an autocrine loop. Further study is required to determine the functions of these fusion proteins.

A quantitative high-resolution genetic profile rapidly identifies sequence determinants of hepatitis C viral fitness and drug sensitivity.

Widely used chemical genetic screens have greatly facilitated the identification of many antiviral agents. However, the regions of interaction and inhibitory mechanisms of many therapeutic candidates have yet to be elucidated. Previous chemical screens identified Daclatasvir (BMS-790052) as a potent nonstructural protein 5A (NS5A) inhibitor for Hepatitis C virus (HCV) infection with an unclear inhibitory mechanism. Here we have developed a quantitative high-resolution genetic (qHRG) approach to systematically map the drug-protein interactions between Daclatasvir and NS5A and profile genetic barriers to Daclatasvir resistance. We implemented saturation mutagenesis in combination with next-generation sequencing technology to systematically quantify the effect of every possible amino acid substitution in the drug-targeted region (domain IA of NS5A) on replication fitness and sensitivity to Daclatasvir. This enabled determination of the residues governing drug-protein interactions. The relative fitness and drug sensitivity profiles also provide a comprehensive reference of the genetic barriers for all possible single amino acid changes during viral evolution, which we utilized to predict clinical outcomes using mathematical models. We envision that this high-resolution profiling methodology will be useful for next-generation drug development to select drugs with higher fitness costs to resistance, and also for informing the rational use of drugs based on viral variant spectra from patients.

Dual-compartmental transcriptomic + proteomic analysis of a marine endosymbiosis exposed to environmental change.

As significant anthropogenic pressures are putting undue stress on the world's oceans, there has been a concerted effort to understand how marine organisms respond to environmental change. Transcriptomic approaches, in particular, have been readily employed to document the mRNA-level response of a plethora of marine invertebrates exposed to an array of simulated stress scenarios, with the tacit and untested assumption being that the respective proteins show a corresponding trend. To better understand the degree of congruency between mRNA and protein expression in an endosymbiotic marine invertebrate, mRNAs and proteins were sequenced from the same samples of the common, Indo-Pacific coral Seriatopora hystrix exposed to stable or upwelling-simulating conditions for 1 week. Of the 167 proteins downregulated at variable temperature, only two were associated with mRNAs that were also differentially expressed between treatments. Of the 378 differentially expressed genes, none were associated with a differentially expressed protein. Collectively, these results highlight the inherent risk of inferring cellular behaviour based on mRNA expression data alone and challenge the current, mRNA-focused approach taken by most marine and many molecular biologists.

Identification of a novel FN1-FGFR1 genetic fusion as a frequent event in phosphaturic mesenchymal tumour.

Phosphaturic mesenchymal tumours (PMTs) are uncommon soft tissue and bone tumours that typically cause hypophosphataemia and tumour-induced osteomalacia (TIO) through secretion of phosphatonins including fibroblast growth factor 23 (FGF23). PMT has recently been accepted by the World Health Organization as a formal tumour entity. The genetic basis and oncogenic pathways underlying its tumourigenesis remain obscure. In this study, we identified a novel FN1-FGFR1 fusion gene in three out of four PMTs by next-generation RNA sequencing. The fusion transcripts and proteins were subsequently confirmed with RT-PCR and western blotting. Fluorescence in situ hybridization analysis showed six cases with FN1-FGFR1 fusion out of an additional 11 PMTs. Overall, nine out of 15 PMTs (60%) harboured this fusion. The FN1 gene possibly provides its constitutively active promoter and the encoded protein's oligomerization domains to overexpress and facilitate the activation of the FGFR1 kinase domain. Interestingly, unlike the prototypical leukaemia-inducing FGFR1 fusion genes, which are ligand-independent, the FN1-FGFR1 chimeric protein was predicted to preserve its ligand-binding domains, suggesting an advantage of the presence of its ligands (such as FGF23 secreted at high levels by the tumour) in the activation of the chimeric receptor tyrosine kinase, thus effecting an autocrine or a paracrine mechanism of tumourigenesis.

Precise genotyping and recombination detection of Enterovirus.

Enteroviruses (EV) with different genotypes cause diverse infectious diseases in humans and mammals. A correct EV typing result is crucial for effective medical treatment and disease control; however, the emergence of novel viral strains has impaired the performance of available diagnostic tools. Here, we present a web-based tool, named EVIDENCE (EnteroVirus In DEep conception, http://symbiont.iis.sinica.edu.tw/evidence), for EV genotyping and recombination detection. We introduce the idea of using mixed-ranking scores to evaluate the fitness of prototypes based on relatedness and on the genome regions of interest. Using phylogenetic methods, the most possible genotype is determined based on the closest neighbor among the selected references. To detect possible recombination events, EVIDENCE calculates the sequence distance and phylogenetic relationship among sequences of all sliding windows scanning over the whole genome. Detected recombination events are plotted in an interactive figure for viewing of fine details. In addition, all EV sequences available in GenBank were collected and revised using the latest classification and nomenclature of EV in EVIDENCE. These sequences are built into the database and are retrieved in an indexed catalog, or can be searched for by keywords or by sequence similarity. EVIDENCE is the first web-based tool containing pipelines for genotyping and recombination detection, with updated, built-in, and complete reference sequences to improve sensitivity and specificity. The use of EVIDENCE can accelerate genotype identification, aiding clinical diagnosis and enhancing our understanding of EV evolution.

Subset selection of high-depth next generation sequencing reads for de novo genome assembly using MapReduce framework.

BACKGROUND: Recent progress in next-generation sequencing technology has afforded several improvements such as ultra-high throughput at low cost, very high read quality, and substantially increased sequencing depth. State-of-the-art high-throughput sequencers, such as the Illumina MiSeq system, can generate ~15 Gbp sequencing data per run, with >80% bases above Q30 and a sequencing depth of up to several 1000x for small genomes. Illumina HiSeq 2500 is capable of generating up to 1 Tbp per run, with >80% bases above Q30 and often >100x sequencing depth for large genomes. To speed up otherwise time-consuming genome assembly and/or to obtain a skeleton of the assembly quickly for scaffolding or progressive assembly, methods for noise removal and reduction of redundancy in the original data, with almost equal or better assembly results, are worth studying. RESULTS: We developed two subset selection methods for single-end reads and a method for paired-end reads based on base quality scores and other read analytic tools using the MapReduce framework. We proposed two strategies to select reads: MinimalQ and ProductQ. MinimalQ selects reads with minimal base-quality above a threshold. ProductQ selects reads with probability of no incorrect base above a threshold. In the single-end experiments, we used Escherichia coli and Bacillus cereus datasets of MiSeq, Velvet assembler for genome assembly, and GAGE benchmark tools for result evaluation. In the paired-end experiments, we used the giant grouper (Epinephelus lanceolatus) dataset of HiSeq, ALLPATHS-LG genome assembler, and QUAST quality assessment tool for comparing genome assemblies of the original set and the subset. The results show that subset selection not only can speed up the genome assembly but also can produce substantially longer scaffolds. AVAILABILITY: The software is freely available at https://github.com/moneycat/QReadSelector.

The effect of red light and far-red light conditions on secondary metabolism in agarwood.

BACKGROUND: Agarwood, a heartwood derived from Aquilaria trees, is a valuable commodity that has seen prevalent use among many cultures. In particular, it is widely used in herbal medicine and many compounds in agarwood are known to exhibit medicinal properties. Although there exists much research into medicinal herbs and extraction of high value compounds, few have focused on increasing the quantity of target compounds through stimulation of its related pathways in this species. RESULTS: In this study, we observed that cucurbitacin yield can be increased through the use of different light conditions to stimulate related pathways and conducted three types of high-throughput sequencing experiments in order to study the effect of light conditions on secondary metabolism in agarwood. We constructed genome-wide profiles of RNA expression, small RNA, and DNA methylation under red light and far-red light conditions. With these profiles, we identified a set of small RNA which potentially regulates gene expression via the RNA-directed DNA methylation pathway. CONCLUSIONS: We demonstrate that light conditions can be used to stimulate pathways related to secondary metabolism, increasing the yield of cucurbitacins. The genome-wide expression and methylation profiles from our study provide insight into the effect of light on gene expression for secondary metabolism in agarwood and provide compelling new candidates towards the study of functional secondary metabolic components.

Compartment-specific transcriptomics in a reef-building coral exposed to elevated temperatures.

Although rising ocean temperatures threaten scleractinian corals and the reefs they construct, certain reef corals can acclimate to elevated temperatures to which they are rarely exposed in situ. Specimens of the model Indo-Pacific reef coral Pocillopora damicornis collected from upwelling reefs of Southern Taiwan were previously found to have survived a 36-week exposure to 30°C, a temperature they encounter infrequently and one that can elicit the breakdown of the coral-dinoflagellate (genus Symbiodinium) endosymbiosis in many corals of the Pacific Ocean. To gain insight into the subcellular pathways utilized by both the coral hosts and their mutualistic Symbiodinium populations to acclimate to this temperature, mRNAs from both control (27°C) and high (30°C)-temperature samples were sequenced on an Illumina platform and assembled into a 236 435-contig transcriptome. These P. damicornis specimens were found to be ~60% anthozoan and 40% microbe (Symbiodinium, other eukaryotic microbes, and bacteria), from an mRNA-perspective. Furthermore, a significantly higher proportion of genes from the Symbiodinium compartment were differentially expressed after two weeks of exposure. Specifically, at elevated temperatures, Symbiodinium populations residing within the coral gastrodermal tissues were more likely to up-regulate the expression of genes encoding proteins involved in metabolism than their coral hosts. Collectively, these transcriptome-scale data suggest that the two members of this endosymbiosis have distinct strategies for acclimating to elevated temperatures that are expected to characterize many of Earth's coral reefs in the coming decades.

Decoding the genome of bloodsucking midge Forcipomyia taiwana (Diptera: Ceratopogonidae): Insights into odorant receptor expansion.

Biting midges, notably those within the Ceratopogonidae family, have long been recognized for their epidemiological significance, both as nuisances and vectors for disease transmission in vertebrates. Despite their impact, genomic insights into these insects, particularly beyond the Culicoides genus, remain limited. In this study, we assembled the Forcipomyia taiwana (Shiraki) genome, comprising 113 scaffolds covering 130.4 Mbps-with the longest scaffold reaching 7.6 Mbps and an N50 value of 2.6 Mbps-marking a pivotal advancement in understanding the genetic architecture of ceratopogonid biting midges. Phylogenomic analyses reveal a shared ancestry between F. taiwana and Culicoides sonorensis Wirth & Jones, dating back approximately 124 million years, and highlight a dynamic history of gene family expansions and contractions within the Ceratopogonidae family. Notably, a substantial expansion of the odorant receptor (OR) gene family was observed, which is crucial for the chemosensory capabilities that govern biting midges' interactions with their environment, including host seeking and oviposition behaviors. The distribution of OR genes across the F. taiwana genome displays notable clusters on scaffolds, indicating localized tandem gene duplication events. Additionally, several collinear regions were identified, hinting at segmental duplications, inversions, and translocations, contributing to the olfactory system's evolutionary complexity. Among the 156 ORs identified in F. taiwana, 134 are biting midge-specific ORs, distributed across three distinct clades, each exhibiting unique motif features that distinguish them from the others. Through weighted gene co-expression network analysis, we correlated distinct gene modules with sex and reproductive status, laying the groundwork for future investigations into the interplay between gene expression and adaptive behaviors in F. taiwana. In conclusion, our study not only highlights the unique olfactory repertoire of ceratopogonid biting midges but also sets the stage for future studies into the genetic underpinnings of their unique biological traits and ecological strategies.

Identification of Sexually Dimorphic Genes in Pectoral Fin as Molecular Markers for Assessing the Sex of Japanese Silver Eels (Anguilla japonica).

The Japanese eel (Anguilla japonica) is an important species in East Asian aquaculture. However, the production of seedlings for this purpose still depends on natural resources, as the commercial production of glass eels is not yet possible. Confusion about the sex of silver eels is one of the factors affecting the success rate of artificial maturation. This study sought to devise a harmless method to precisely assess the sex of silver eels. Partial pectoral fins were collected from females and males and the total RNA was extracted for transcriptomic analysis to identify sexually dimorphic genes as molecular markers for sex typing. An online database was constructed to integrate the annotations of transcripts and perform comparative transcriptome analysis. This analysis identified a total of 29 candidate sexually dimorphic genes. Ten were selected for a real-time quantitative polymerase chain reaction (RT-qPCR) to validate the transcriptomic data and evaluate their feasibility as markers. The transcriptomic analysis and RT-qPCR data implicated three potential markers (LOC111853410, kera, and dcn) in sex typing. The expression of LOC111853410 was higher in females than in males. In contrast, the expression of kera and dcn was higher in males than in females. The ΔCT values of three markers were analyzed to determine their inferred thresholds, which can be used to determine the sex of Japanese eels. The results suggested that if a silver eel had a pectoral fin with the pectoral fin having the ΔCT of LOC111853410 < 11.3, the ΔCT of kera > 11.4, or the ΔCT of dcn > 6.5 can be assessed it could be assessed as female. Males could be assessed by the ΔCT of LOC111853410 > 11.3, the ΔCT of kera < 11.4, or the ΔCT of dcn < 6.5 in their pectoral fins. The molecular functions of these markers and the biological significance of their differential expression require further exploration.