Maser: one-stop platform for NGS big data from analysis to visualization

A unique substrate specificity of PonAAS2, an aromatic aldehyde synthase, involved in phytohormone auxin biosynthesis in a gall-inducing sawfly Euura sp. "Pontania"

The aromatic aldehyde synthase (AAS), PonAAS2, from the gall-inducing sawfly has been identified as a biosynthetic enzyme for indole-3-acetic acid (IAA), a key molecule of the plant hormone auxin, which is thought to play a role in gall induction. Unlike other insect AASs that convert Dopa, PonAAS2 uniquely converts l-tryptophan (Trp) into indole-3-acetaldehyde, a precursor of IAA. In this study, an examination of AAS enzymes from various insect species revealed that the ability to convert Trp has been acquired in only a very limited taxonomic group. Comparative analysis between PonAAS2 and DjAAS2 from a gall wasp showed that, despite having conserved substrate-recognizing amino acids, they exhibit different substrate specificities. This difference likely arises from variations in how these enzymes' monomers interact during dimer formation, as demonstrated by amino acid substitution experiments and structural predictions.

Ninein, a candidate gene for ethanol anxiolysis, shows complex exon-specific expression and alternative splicing differences between C57BL/6J and DBA/2J mice

Ethanol's anxiolytic actions contribute to increased consumption and the development of Alcohol Use Disorder (AUD). Our laboratory previously identified genetic loci contributing to the anxiolytic-like properties of ethanol in BXD recombinant inbred mice, derived from C57BL/6J (B6) and DBA/2J (D2) progenitor strains. That work identified Ninein (Nin) as a candidate gene underlying ethanol's acute anxiolytic-like properties in BXD mice. Nin has a complex exonic content with known alternative splicing events that alter cellular distribution of the NIN protein. We hypothesize that strain-specific differences in Nin alternative splicing contribute to changes in Nin gene expression and B6/D2 strain differences in ethanol anxiolysis. Using quantitative reverse-transcriptase PCR to target specific Nin splice variants, we identified isoform-specific exon expression differences between B6 and D2 mice in prefrontal cortex, nucleus accumbens and amygdala. We extended this analysis using deep RNA sequencing in B6 and D2 nucleus accumbens samples and found that total Nin expression was significantly higher in D2 mice. Furthermore, exon utilization and alternative splicing analyses identified eight differentially utilized exons and significant exon-skipping events between the strains, including three novel splicing events in the 3' end of the Nin gene that were specific to the D2 strain. Additionally, we document multiple single nucleotide polymorphisms in D2 Nin exons that are predicted to have deleterious effects on protein function. Our studies provide the first in-depth analysis of Nin alternative splicing in brain and identify a potential genetic mechanism altering Nin expression and function between B6 and D2 mice, thus possibly contributing to differences in the anxiolytic-like properties of ethanol between these strains. This work adds novel information to our understanding of genetic differences modulating ethanol actions on anxiety that may contribute to the risk for alcohol use disorder.

Characterization of unique EDTA-insensitive methylthioalkylmalate synthase from Eutrema japonicum and its potential application in synthetic biology

6-(Methylsulfinyl)hexyl isothiocyanate (6-MSITC), a derivative of glucosinolate with a six-carbon chain, is a compound found in wasabi and has diverse health-promoting properties. The biosynthesis of glucosinolates from methionine depends on a crucial step catalyzed methylthioalkylmalate synthases (MAMs), which are responsible for the generation of glucosinolates with varying chain lengths. In this study, our primary focus was the characterization of two methylthioalkyl malate synthases, MAM1-1 and MAM1-2, derived from Eutrema japonicum, commonly referred to as Japanese wasabi. Eutremajaponicum MAMs (EjMAMs) were expressed in an Escherichiacoli expression system, subsequently purified, and in vitro enzymatic activity was assayed. We explored the kinetic properties, optimal pH conditions, and cofactor preferences of EjMAMs and compared them with those of previously documented MAMs. Surprisingly, EjMAM1-2, categorized as a metallolyase family enzyme, displayed 20% of its maximum activity even in the absence of divalent metal cofactors or under high concentrations of EDTA. Additionally, we utilized AlphaFold2 to generate structural homology models of EjMAMs, and used in silico analysis and mutagenesis studies to investigate the key residues participating in catalytic activity. Moreover, we examined in vivo biosynthesis in E. coli containing Arabidopsis thaliana branched-chain amino acid transferase 3 (AtBCAT3) along with AtMAMs or EjMAMs and demonstrated that EjMAM1-2 exhibited the highest conversion rate among those MAMs, converting l-methionine to 2-(2-methylthio) ethyl malate (2-(2-MT)EM). EjMAM1-2 shows a unique property in vitro and highest activity on converting l-methionine to 2-(2-MT)EM in vivo which displays high potential for isothiocyanate biosynthesis in E. coli platform.

Genomic properties of a Bartonella quintana strain from Japanese macaque (Macaca fuscata) revealed by genome comparison with human and rhesus macaque strains

Bartonella quintana, the causative agent of trench fever, is an intracellular bacterium that infects human erythrocytes and vascular endothelial cells. For many years, humans were considered the only natural hosts for B. quintana; however, it was recently discovered that wild Japanese macaques (Macaca fuscata) also serve as hosts for B. quintana. To elucidate the genetic characteristics of the B. quintana strain MF1-1 isolated from a Japanese macaque, we determined the complete genome sequence of the strain and compared it with those of strain Toulouse from a human and strain RM-11 from a rhesus macaque. General genomic features and orthologous gene cluster profiles are similar among the three strains, and strain MF1-1 is genetically closer to strain RM-11 than strain Toulouse based on the average nucleotide identity values; however, a significant inversion of approximately 0.68 Mb was detected in the chromosome of strain MF1-1. Moreover, the Japanese macaque strains lacked the bepA gene, which is responsible for anti-apoptotic function, and the trwL2, trwL4, and trwL6 genes, which may be involved in adhesion to erythrocytes of rhesus macaque and human. These features likely represent the genomic traits acquired by Japanese macaque strains in their host-associated evolution.

Integrating mitochondrial and nuclear genomic data to decipher the evolutionary history of Eubranchipus species in Japan

Understanding the genetic diversity and evolutionary history of species is crucial for their conservation and management. In this study, we investigated the genetic diversity and phylogenetic relationships among Eubranchipus species occurring in Japan. Phylogenetic analyses revealed that nuclear and mitochondrial data yield incompatible results. In E. uchidai, nuclear data support the monophyly of the Shimokita area, while mitochondrial data indicate a clustering of Higashidori2 individuals with Hokkaido (Ishikari and Wakkanai) E. uchidai. Similar incongruences were observed in E. hatanakai, where nuclear data favor the monophyly of the Chokai area, while mitochondrial data cluster some Chokai pool 3 individuals with Aizu individuals. These incompatibilities might be caused by mitochondrial gene flow. The findings emphasize the importance of considering both nuclear and mitochondrial data during phylogenetic studies and provide valuable insights into the complex dynamics of migration and genetic exchange in Eubranchipus species.

Identification of an operon and its regulator required for autoaggregation in Tetragenococcus halophilus

IMPORTANCE

Tetragenococcus halophilus is a halophilic lactic acid bacterium generally used as a starter culture in fermenting soy and fish sauces. Aggregating strains can be useful in fermenting and obtaining clear soy sauce because cell clumps are trapped by the filter cake when the soy sauce mash is pressed. However, the genetic mechanisms of aggregation in T. halophilus are unknown. In this study, we identified genes encoding aggregation factor and its regulator. These findings may provide a foundation for developing improved T. halophilus starter cultures for soy sauce fermentation, leading to more efficient and consistent clear soy sauce production.

A newly identified enzyme from Japanese common squid Todarodes pacificus has the ability to biosynthesize d-aspartate

Amino acids exist in two chiral forms, namely L and D. Although l-amino acids are predominant in vivo, certain limited circumstances have reported the usage of d-amino acids. d-aspartate (Asp), among them, plays crucial physiological roles in living organisms and is biosynthesized from L-Asp by the enzyme named aspartate racemase (AspRase). D-Asp is known to accumulate in large amounts in the nervous system of cephalopods. To understand the function of D-Asp in nervous system in more detail, it is necessary to elucidate its metabolic pathway; however, AspRase gene has not been identified in cephalopods as in the case of mammals. In this study, we successfully identified a novel gene encoding AspRase from the optic ganglion of Japanese common squid Todarodes pacificus. Our discovery of the squid AspRase challenges the prevailing assumption that AspRases across different animals share similar structures. Surprisingly, the squid AspRase is a unique enzyme that differs significantly from known AspRases, being structurally and phylogenetically related to aspartate aminotransferase (AST) and possessing both AspRase and AST activities. The optimum pH and temperature for AspRase activity using L-Asp as a substrate are approximately 7.0 and 20 °C, respectively. Moreover, we have found that AspRase activity is enhanced in the presence of 2-oxoacids. These findings have far-reaching implications for the understanding of enzymology and suggest that yet-to-be-identified mammalian AspRases may also be phylogenetically related to AST, rather than conventional AspRases. Furthermore, our results provide valuable insights into the evolution of the D-Asp biosynthetic pathway.

Social isolation shortens lifespan through oxidative stress in ants

Social isolation negatively affects health, induces detrimental behaviors, and shortens lifespan in social species. Little is known about the mechanisms underpinning these effects because model species are typically short-lived and non-social. Using colonies of the carpenter ant Camponotus fellah, we show that social isolation induces hyperactivity, alters space-use, and reduces lifespan via changes in the expression of genes with key roles in oxidation-reduction and an associated accumulation of reactive oxygen species. These physiological effects are localized to the fat body and oenocytes, which perform liver-like functions in insects. We use pharmacological manipulations to demonstrate that the oxidation-reduction pathway causally underpins the detrimental effects of social isolation on behavior and lifespan. These findings have important implications for our understanding of how social isolation affects behavior and lifespan in general.

Symbiosis of Carpenter Bees with Uncharacterized Lactic Acid Bacteria Showing NAD Auxotrophy

Eusocial bees (such as honey bees and bumble bees) harbor core gut microbiomes that are transmitted through social interaction between nestmates. Carpenter bees are not eusocial; however, recent microbiome analyses found that Xylocopa species harbor distinctive core gut microbiomes. In this study, we analyzed the gut microbiomes of three Xylocopa species in Japan between 2016 and 2021 by V1 to V2 region-based 16S rDNA amplicon sequencing, and 14 candidate novel species were detected based on the full-length 16S rRNA gene sequences. All Xylocopa species harbor core gut microbiomes consisting of primarily lactic acid bacteria (LAB) that were phylogenetically distant from known species. Although they were difficult to cultivate, two LAB species from two different Xylocopa species were isolated by supplementing bacterial culture supernatants. Both genomes exhibited an average LAB genome size with a large set of genes for carbohydrate utilization but lacked genes to synthesize an essential coenzyme NAD, which is unique among known insect symbionts. Our findings of phylogenetically distinct core LAB of NAD auxotrophy reflected the evolution of Xylocopa-restricted bacteria retention and maintenance through vertical transmission of microbes during solitary life. We propose five candidate novel species belonging to the families Lactobacillaceae and Bifidobacteriaceae, including a novel genus, and their potential functions in carbohydrate utilization. IMPORTANCE Recent investigations found unique microbiomes in carpenter bees, but the description of individual microbes, including isolation and genomics, remains largely unknown. Here, we found that the Japanese Xylocopa species also harbor core gut microbiomes. Although most of them were difficult to isolate a pure colony, we successfully isolated several strains. We performed whole-genome sequencing of the isolated candidate novel species and found that the two Lactobacillaceae strains belonging to the Xylocopa-specific novel LAB clade lack the genes for synthesizing NAD, a coenzyme central to metabolism in all living organisms. Here, we propose a novel genus for the two LAB species based on very low 16S rRNA gene sequence similarities and genotypic characters.

Resources and tools for rare disease variant interpretation

Collectively, rare genetic disorders affect a substantial portion of the world's population. In most cases, those affected face difficulties in receiving a clinical diagnosis and genetic characterization. The understanding of the molecular mechanisms of these diseases and the development of therapeutic treatments for patients are also challenging. However, the application of recent advancements in genome sequencing/analysis technologies and computer-aided tools for predicting phenotype-genotype associations can bring significant benefits to this field. In this review, we highlight the most relevant online resources and computational tools for genome interpretation that can enhance the diagnosis, clinical management, and development of treatments for rare disorders. Our focus is on resources for interpreting single nucleotide variants. Additionally, we present use cases for interpreting genetic variants in clinical settings and review the limitations of these results and prediction tools. Finally, we have compiled a curated set of core resources and tools for analyzing rare disease genomes. Such resources and tools can be utilized to develop standardized protocols that will enhance the accuracy and effectiveness of rare disease diagnosis.

Rapid whole genome sequencing methods for RNA viruses

RNA viruses are the etiological agents of many infectious diseases. Since RNA viruses are error-prone during genome replication, rapid, accurate and economical whole RNA viral genome sequence determination is highly demanded. Next-generation sequencing (NGS) techniques perform whole viral genome sequencing due to their high-throughput sequencing capacity. However, the NGS techniques involve a significant burden for sample preparation. Since to generate complete viral genome coverage, genomic nucleic acid enrichment is required by reverse transcription PCR using virus-specific primers or by viral particle concentration. Furthermore, conventional NGS techniques cannot determine the 5' and 3' terminal sequences of the RNA viral genome. Therefore, the terminal sequences are determined one by one using rapid amplification of cDNA ends (RACE). However, since some RNA viruses have segmented genomes, the burden of the determination using RACE is proportional to the number of segments. To date, there is only one study attempting whole genome sequencing of multiple RNA viruses without using above mentioned methods, but the generated sequences' accuracy compared to the reference sequences was up to 97% and did not reach 100% due to the low read depth. Hence, we established novel methods, named PCR-NGS and RCA-NGS, that were optimized for an NGS machine, MinION. These methods do not require nucleic acid amplification with virus-specific PCR primers, physical viral particle enrichment, and RACE. These methods enable whole RNA viral genome sequencing by combining the following techniques: (1) removal of unwanted DNA and RNA other than the RNA viral genome by nuclease treatment; (2) the terminal of viral genome sequence determination by barcoded linkers ligation; (3) amplification of the viral genomic cDNA using ligated linker sequences-specific PCR or an isothermal DNA amplification technique, such as rolling circle amplification (RCA). The established method was evaluated using isolated RNA viruses with single-stranded, double-stranded, positive-stranded, negative-stranded, non-segmented or multi-segmented genomes. As a result, all the viral genome sequences could be determined with 100% accuracy, and these mean read depths were greater than 2,500×, at least using either of the methods. This method should allow for easy and economical determination of accurate RNA viral genomes.

Changes of the gene expression in silkworm larvae and Cordyceps militaris at late stages of the pathogenesis

Cordyceps militaris is an entomopathogenic fungus that forms its fruiting body. The gene expression change in C. militaris and silkworm larvae were analyzed using RNA-seq to investigate the relationship of C. militaris with the host, silkworm larvae before the death by mycosis. At 144 h after the injection of C. militaris conidia, genes encoding proteases, protease inhibitors, and cuticle proteins in the fat body of silkworm larvae were upregulated, but genes encoding lipoproteins and other proteins in hemolymph were downregulated. On the other hand, at 168 h after the injection of C. militaris conidia, genes encoding amino acid and oligopeptide transporters and permeases in C. militaris were upregulated, suggesting that C. militaris may use peptides and amino acids in silkworm larvae as a nutrient to grow in vivo. Additionally, one gene cluster composed of genes putatively involved in the degradation of phenolic substrates was also upregulated. The addition of 4,5-dichlorocatechol, an inhibitor of catechol 1,2-dioxygenase, inhibited the in vivo growth of C. militaris, Beauveria bassiana and Metarhizium anisopliae. These results also suggest that the expression of the gene cluster may be crucial for the in vivo growth of C. militaris and entomopathogenic fungi. This study will clarify how C. militaris grows in insect hosts by avoiding host's immune systems.

Estrogen-Inducible LncRNA BNAT1 Functions as a Modulator for Estrogen Receptor Signaling in Endocrine-Resistant Breast Cancer Cells

Recent advances in RNA studies have revealed that functional long noncoding RNAs (lncRNAs) contribute to the biology of cancers. In breast cancer, estrogen receptor α (ERα) is an essential transcription factor that primarily promotes the growth of luminal-type cancer, although only a small number of lncRNAs are identified as direct ERα targets and modulators for ERα signaling. In this study, we performed RNA-sequencing for ER-positive breast cancer cells and identified a novel estrogen-inducible antisense RNA in the COL18A1 promoter region, named breast cancer natural antisense transcript 1 (BNAT1). In clinicopathological study, BNAT1 may have clinical relevance as a potential diagnostic factor for prognoses of ER-positive breast cancer patients based on an in situ hybridization study for breast cancer specimens. siRNA-mediated BNAT1 silencing significantly inhibited the in vitro and in vivo growth of tamoxifen-resistant ER-positive breast cancer cells. Notably, BNAT1 silencing repressed cell cycle progression whereas it promoted apoptosis. Microarray analysis revealed that BNAT1 silencing in estrogen-sensitive breast cancer cells repressed estrogen signaling. We showed that BNAT1 knockdown decreased ERα expression and repressed ERα transactivation. RNA immunoprecipitation showed that BNAT1 physically binds to ERα protein. In summary, BNAT1 would play a critical role in the biology of ER-positive breast cancer by modulating ERα-dependent transcription regulation. We consider that BNAT1 could be a potential molecular target for diagnostic and therapeutic options targeting luminal-type and endocrine-resistant breast cancer.

Phagocytosing differentiated cell-fragments is a novel mechanism for controlling somatic stem cell differentiation within a short time frame

Stem cells undergo cytokine-driven differentiation, but this process often takes longer than several weeks to complete. A novel mechanism for somatic stem cell differentiation via phagocytosing ‘model cells’ (apoptotic differentiated cells) was found to require only a short time frame. Pluripotent-like Muse cells, multipotent mesenchymal stem cells (MSCs), and neural stem cells (NSCs) phagocytosed apoptotic differentiated cells via different phagocytic receptor subsets than macrophages. The phagocytosed-differentiated cell-derived contents (e.g., transcription factors) were quickly released into the cytoplasm, translocated into the nucleus, and bound to promoter regions of the stem cell genomes. Within 24 ~ 36 h, the cells expressed lineage-specific markers corresponding to the phagocytosed-differentiated cells, both in vitro and in vivo. At 1 week, the gene expression profiles were similar to those of the authentic differentiated cells and expressed functional markers. Differentiation was limited to the inherent potential of each cell line: triploblastic-, adipogenic-/chondrogenic-, and neural-lineages in Muse cells, MSCs, and NSCs, respectively. Disruption of phagocytosis, either by phagocytic receptor inhibition via small interfering RNA or annexin V treatment, impeded differentiation in vitro and in vivo. Together, our findings uncovered a simple mechanism by which differentiation-directing factors are directly transferred to somatic stem cells by phagocytosing apoptotic differentiated cells to trigger their rapid differentiation into the target cell lineage.

AnFiSA: An open-source computational platform for the analysis of sequencing data for rare genetic disease

Despite genomic sequencing rapidly transforming from being a bench-side tool to a routine procedure in a hospital, there is a noticeable lack of genomic analysis software that supports both clinical and research workflows as well as crowdsourcing. Furthermore, most existing software packages are not forward-compatible in regards to supporting ever-changing diagnostic rules adopted by the genetics community. Regular updates of genomics databases pose challenges for reproducible and traceable automated genetic diagnostics tools. Lastly, most of the software tools score low on explainability amongst clinicians. We have created a fully open-source variant curation tool, AnFiSA, with the intention to invite and accept contributions from clinicians, researchers, and professional software developers. The design of AnFiSA addresses the aforementioned issues via the following architectural principles: using a multidimensional database management system (DBMS) for genomic data to address reproducibility, curated decision trees adaptable to changing clinical rules, and a crowdsourcing-friendly interface to address difficult-to-diagnose cases. We discuss how we have chosen our technology stack and describe the design and implementation of the software. Finally, we show in detail how selected workflows can be implemented using the current version of AnFiSA by a medical geneticist.

Identification of Candidate Genes Responsible for Flower Colour Intensity in Gentiana triflora

Gentians cultivated in Japan (Gentiana triflora and Gentiana scabra and hybrids) have blue flowers, but flower colour intensity differs among cultivars. The molecular mechanism underlying the variation in flower colour intensity is unclear. Here, we produced F₂ progeny derived from an F₁ cross of intense- and faint-blue lines and attempted to identify the genes responsible for flower colour intensity using RNA-sequencing analyses. Comparative analysis of flower colour intensity and transcriptome data revealed differentially expressed genes (DEGs), although known flavonoid biosynthesis-related genes showed similar expression patterns. From quantitative RT-PCR (qRT-PCR) analysis, we identified two and four genes with significantly different expression levels in the intense- and faint-blue flower lines, respectively. We conducted further analyses on one of the DEGs, termed GtMIF1, which encodes a putative mini zinc-finger protein homolog, which was most differently expressed in faint-blue individuals. Functional analysis of GtMIF1 was performed by producing stable tobacco transformants. GtMIF1-overexpressing tobacco plants showed reduced flower colour intensity compared with untransformed control plants. DNA-marker analysis also confirmed that the GtMIF1 allele of the faint-blue flower line correlated well with faint flower colour in F₂ progeny. These results suggest that GtMIF1 is one of the key genes involved in determining the flower colour intensity of gentian.

Exploration of Pericyte-Derived Factors Implicated in Lung Cancer Brain Metastasis Protection: A Pilot Messenger RNA Sequencing Using the Blood-Brain Barrier In Vitro Model

Metastatic brain tumors have poor prognoses and pose unmet clinical problems for the patients. The blood-brain barrier (BBB) implication is supposed to play a major role in brain metastasis. However, the role of pericytes remains to be elucidated in the brain metastasis. This pilot study described the expression profile of interactions between pericytes, endothelial cells, and cancer cells. We applied an in vitro BBB model with rat primary cultured BBB-related cells (endothelial cells and pericytes), and performed the gene expression analyses of pericytes under the lung cancer cells coculture conditions. Pericytes demonstrated inhibition of the cancer cell proliferation significantly (p < 0.05). RNA was extracted from the pericytes, complementary DNA library was prepared, and RNA-seq was performed. The sequence read data were analyzed on the Management and Analysis System for Enormous Reads and Tag Count Comparison-Graphical User Interface platforms. No statistically or biologically significant differentially expressed genes (DEGs) were detected in the explanatory analyses. Lot-specific DEG detection demonstrated significant decreases in the expression of two genes (Wwtr1 and Acin1), and enrichment analyses using Metascape software revealed the inhibition of apoptotic processes in fibroblasts. Our results suggest that the expression profiles of brain pericytes are partially implicated in the prevention of lung cancer metastasis to the brain. Pericytes exerted an anti-metastatic effect in the BBB model, and their neurohumoral factors remain to be elucidated.

Complete mitochondrial genomes of three fairy shrimps from snowmelt pools in Japan

Background

Fairy shrimps belong to order Anostraca, class Branchiopoda, subphylum Crustacea, and phylum Arthropoda. Three fairy shrimp species (Eubranchipus uchidai, E. asanumai, and E. hatanakai) that inhabit snowmelt pools are currently known in Japan. Whole mitochondrial genomes are useful genetic information for conducting phylogenetic analyses. Mitochondrial genome sequences for Branchiopoda members are gradually being collated.

Results

Six whole mitochondrial genomes from the three Eubranchipus species are presented here. Eubranchipus species share the anostracan pattern of gene arrangement in their mitochondrial genomes. The mitochondrial genomes of the Eubranchipus species have a higher GC content than those of other anostracans. Accelerated substitution rates in the lineage of Eubranchipus species were observed.

Conclusion

This study is the first to obtain whole mitochondrial genomes for Far Eastern Eubranchipus species. We show that the nucleotide sequences of cytochrome oxidase subunit I and the 16S ribosomal RNA of E. asanumai presented in a previous study were nuclear mitochondrial DNA segments. Higher GC contents and accelerated substitution rates are specific characteristics of the mitochondrial genomes of Far Eastern Eubranchipus. The results will be useful for further investigations of the evolution of Anostraca as well as Branchiopoda.

Differential Biosynthesis and Roles of Two Ferrichrome-Type Siderophores, ASP2397/AS2488053 and Ferricrocin, in Acremonium persicinum

Ferrichromes are a family of fungal siderophores with cyclic hexapeptide structures. Most fungi produce one or two ferrichrome-type siderophores. Acremonium persicinum MF-347833 produces ferrichrome-like potent Trojan horse antifungal antibiotics ASP2397 and AS2488053, the aluminum- and iron-chelating forms of AS2488059, respectively. Here, we show by gene sequencing followed by gene deletion experiments that A. persicinum MF-347833 possesses two nonribosomal peptide synthetase genes responsible for AS2488059 and ferricrocin assembly. AS2488059 was produced under iron starvation conditions and excreted into the media to serve as a defense metabolite and probably an iron courier. In contrast, ferricrocin was produced under iron-replete conditions and retained inside the cells, likely serving as an iron-sequestering molecule. Notably, the phylogenetic analyses suggest the different evolutionary origin of AS2488059 from that of conventional ferrichrome-type siderophores. Harnessing two ferrichrome-type siderophores with distinct biological properties may give A. persicinum a competitive advantage for surviving the natural environment.

RNA Sequencing Data Analysis on the Maser Platform and the Tag-Count Comparison Graphical User Interface

The RNA sequencing (RNA-seq) process that allows for comprehensive transcriptome analysis has become increasingly simple. Analysis and interpretation of RNA-seq output data are indispensable for research, but bioinformatics experts are not always available to assist. Currently, however, even a wet-lab specialist can perform the pipeline analysis of RNA-seq described in this chapter using the Maser platform and the Tag-Count Comparison Graphical User Interface (TCC-GUI). These are free of charge for scientific use.

Hydrophilic Shell Matrix Proteins of Nautilus pompilius and the Identification of a Core Set of Conchiferan Domains

Despite being a member of the shelled mollusks (Conchiferans), most members of extant cephalopods have lost their external biomineralized shells, except for the basally diverging Nautilids. Here, we report the result of our study to identify major Shell Matrix Proteins and their domains in the Nautilid Nautilus pompilius, in order to gain a general insight into the evolution of Conchiferan Shell Matrix Proteins. In order to do so, we performed a multiomics study on the shell of N. pompilius, by conducting transcriptomics of its mantle tissue and proteomics of its shell matrix. Analyses of obtained data identified 61 distinct shell-specific sequences. Of the successfully annotated 27 sequences, protein domains were predicted in 19. Comparative analysis of Nautilus sequences with four Conchiferans for which Shell Matrix Protein data were available (the pacific oyster, the pearl oyster, the limpet and the Euhadra snail) revealed that three proteins and six protein domains were conserved in all Conchiferans. Interestingly, when the terrestrial Euhadra snail was excluded, another five proteins and six protein domains were found to be shared among the four marine Conchiferans. Phylogenetic analyses indicated that most of these proteins and domains were probably present in the ancestral Conchiferan, but employed in shell formation later and independently in most clades. Even though further studies utilizing deeper sequencing techniques to obtain genome and full-length sequences, and functional analyses, must be carried out in the future, our results here provide important pieces of information for the elucidation of the evolution of Conchiferan shells at the molecular level.

The genetic variation and relationship among the natural hybrids of Mangifera casturi Kosterm

Mangifera casturi Kosterm., a mango plant from Kalimantan Selatan, Indonesia, has limited genetic information, severely limiting the research on its genetic variation and phylogeny. We collected M. casturi's genomic information using next-generation sequencing, developed microsatellite markers and performed Sanger sequencing for DNA barcoding analysis. These markers were used to confirm parental origin and genetic diversity of M. casturi hybrids. The clean reads of the Kasturi accession were assembled de novo, producing 259 872 scaffolds (N50 = 1 445 bp). Fourteen polymorphic microsatellite markers were developed from 11 040 microsatellite motif-containing sequences. In total, 58 alleles were produced with a mean of 4.14 alleles per locus. Microsatellite marker analysis revealed broad genetic variation in M. casturi. Phylogenetic analysis was performed using internal transcribed spacers (ITS), matK, rbcL, and trnH-psbA. The phylogenetic tree of chloroplast markers placed Kasturi, Cuban, Pelipisan, Pinari, and Hambawang in one group, with M. indica as the female ancestor. Meanwhile, the phylogenetic tree of ITS markers indicated several Mangifera species as ancestors of M. casturi. Thus, M. casturi very likely originated from the cross-hybridization of multiple ancestors. Furthermore, crossing the F1 hybrids of M. indica and M. quadrifida with other Mangifera spp. may have generated much genetic variation. The genetic information for M. casturi will be a resource for breeding improvement, and conservation studies.

Tomato E8 Encodes a C-27 Hydroxylase in Metabolic Detoxification of α-Tomatine during Fruit Ripening

Tomato (Solanum lycopersicum) contains α-tomatine, a steroidal glycoalkaloid that contributes to the plant defense against pathogens and herbivores through its bitter taste and toxicity. It accumulates at high levels in all the plant tissues, especially in leaves and immature green fruits, whereas it decreases during fruit ripening through metabolic conversion to the nontoxic esculeoside A, which accumulates in the mature red fruit. This study aimed to identify the gene encoding a C-27 hydroxylase that is a key enzyme in the metabolic conversion of α-tomatine to esculeoside A. The E8 gene, encoding a 2-oxoglutalate-dependent dioxygenase, is well known as an inducible gene in response to ethylene during fruit ripening. The recombinant E8 was found to catalyze the C-27 hydroxylation of lycoperoside C to produce prosapogenin A and is designated as Sl27DOX. The ripe fruit of E8/Sl27DOX-silenced transgenic tomato plants accumulated lycoperoside C and exhibited decreased esculeoside A levels compared with the wild-type (WT) plants. Furthermore, E8/Sl27DOX deletion in tomato accessions resulted in higher lycoperoside C levels in ripe fruits than in WT plants. Thus, E8/Sl27DOX functions as a C-27 hydroxylase of lycoperoside C in the metabolic detoxification of α-tomatine during tomato fruit ripening, and the efficient detoxification by E8/27DOX may provide an advantage in the domestication of cultivated tomatoes.

Violet bioluminescent Polycirrus sp. (Annelida: Terebelliformia) discovered in the shallow coastal waters of the Noto Peninsula in Japan

Terebellidae worms have large numbers of tentacles responsible for various biological functions. Some Terebellidae worms whose tentacles emit light are found around the world, including exceptional violet-light-emitting Polycirrus spp. found in Europe and North America. However, there is no video-recorded observation of the luminous behavior of such unique species in nature, and the genetic information related to their ecology are lacking. Here, for the first time, we video-recorded the violet-light-emitting behavior of an undescribed Japanese worm in its natural habitat. The worm was designated as Polycirrus sp. ISK based on morphological observations, and the luminescence spectrum showed a peak at 444 nm, which is an exceptionally short wavelength for bioluminescence in a shallow coastal water environment. An analysis of differentially expressing genes based on separate RNA-Seq analysis for the tentacles and the rest of body revealed the specific expression of genes that are probably involved in innate immunity in the tentacles exposed to predators. We also found a Renilla luciferase homologous gene, but coelenterazine was not detected in the worm extract by analyses using a liquid chromatography and a recombinant Renilla luciferase. These results will promote an understanding of the ecology and luminescence mechanisms of luminous Polycirrus spp.

An efficient i-GONAD method for creating and maintaining lethal mutant mice using an inversion balancer identified from the C3H/HeJJcl strain

As the efficiency of the clustered regularly interspaced short palindromic repeats/Cas system is extremely high, creation and maintenance of homozygous lethal mutants are often difficult. Here, we present an efficient in vivo electroporation method called improved genome editing via oviductal nucleic acid delivery (i-GONAD), wherein one of two alleles in the lethal gene was selectively edited in the presence of a non-targeted B6.C3H-In(6)1J inversion identified from the C3H/HeJJcl strain. This method did not require isolation, culture, transfer, or other in vitro handling of mouse embryos. The edited lethal genes were stably maintained in heterozygotes, as recombination is strongly suppressed within this inversion interval. Using this strategy, we successfully generated the first Tprkb null knockout strain with an embryonic lethal mutation and showed that B6.C3H-In(6)1J can efficiently suppress recombination. As B6.C3H-In(6)1J was tagged with a gene encoding the visible coat color marker, Mitf, the Tprkb mutation could be visually recognized. We listed the stock balancer strains currently available as public bioresources to create these lethal gene knockouts. This method will allow for more efficient experiments for further analysis of lethal mutants.

Changes in Cell Wall Structure During Rhizoid Formation of Silvetia babingtonii (Fucales, Phaeophyceae) Zygotes

We examined the ultrastructure of the cell wall and immunolocalization of alginates using specific antibodies against M-rich alginates and MG blocks during rhizoid formation in fucoid zygotes, Silvetia babingtonii. The thallus region of 24-h-old zygotes had a cell wall made of three layers with different fiber distribution. In the 12-h-old zygotes, three layers in the thallus were observed before rhizoid formation, namely the inner, middle, and outer layers. During rhizoid elongation, only the inner layer was apparent close to the rhizoid tip area. Immunoelectron microscopy detected M-rich blocks of alginate on the inner half of the cell wall, irrespective of the number of layers in the thallus and rhizoid regions. The MG blocks were seen to cover a slightly wider area than M-rich alginate blocks. It was suggested that parts of M in mannuronan would be rapidly converted to G, and MG-blocks are generated. Transcriptome analysis was performed using 3 -, 10 -, and 24-h-old zygotes after fertilization to examine the relationship between gene expression and alginate synthesis over time. The expression of two mannuronan C5-epimerase homologs that convert mannuronic acid into guluronic acid in alginates was upregulated or downregulated over the course of the examination.

Identification and Characterization of Genes for the Allo A Lectins in Japanese Rhinoceros Beetle (Trypoxylus dichotomus [Allomyrina dichotoma])

The Japanese rhinoceros beetle (Trypoxylus dichotomus [Allomyrina dichotoma]) produces the lectins allo A-I and allo A-II, which have strong N-acetyllactosamine (Galβ1-4GlcNAc)-binding activity. It has been suggested that the two lectins are formed from three subunits (α, β, and γ), with allo A-I comprising α and γ subunits and allo A-II comprising β and γ subunits. Here, we determined the cDNA sequences of these subunits using both conventional polymerase chain reaction (PCR)-cloning-sequencing and transcriptome-sequencing analyses. For the α and β subunits, one gene (locus) for each was predicted, whereas for the γ subunit, two types of cDNA sequences were obtained, which we named γ1 and γ2. These two types probably have distinct loci. Average nucleotide sequence identities among the subunits ranged from 87.6% (between α and γ1) to 92.6% (between γ1 and γ2), suggesting that they form a gene family. Although no homology was found between the sequences of allo A and other known lectin proteins in a protein database search, some unknown proteins containing the DUF3421 domain were identified. Those DUF3421 domain-encoding proteins are upregulated in the insect larval midgut. Thus, we infer that allo A genes also play an important role in larvae and that their lectin activity may have been obtained collaterally.

Effects of Cordycepin in Cordyceps militaris during Its Infection to Silkworm Larvae

Cordyceps militaris produces cordycepin, a secondary metabolite that exhibits numerous bioactive properties. However, cordycepin pharmacology in vivo is not yet understood. In this study, the roles of cordycepin in C. militaris during its infection were investigated. After the injection of conidia, C. militaris NBRC100741 killed silkworm larvae more rapidly than NBRC103752. At 96 and 120 h, Cmcns genes (Cmcns1-4), which are part of the cordycepin biosynthesis gene cluster, were expressed in fat bodies and cuticles. Thus, cordycepin may be produced in the infection of silkworm larvae. Further, cordycepin enhanced pathogenicity toward silkworm larvae of Metarhizium anisopliae and Beauveria bassiana, that are also entomopathogenic fungi and do not produce cordycepin. In addition, by RNA-seq analysis, the increased expression of the gene encoding a lipoprotein 30K-8 (Bmlp20, KWMTBOMO11934) and decreased expression of genes encoding cuticular proteins (KWMTBOMO13140, KWMTBOMO13167) and a serine protease inhibitor (serpin29, KWMTBOMO08927) were observed when cordycepin was injected into silkworm larvae. This result suggests that cordycepin may aid the in vivo growth of C. militaris in silkworm larvae by the influence of the expression of some genes in silkworm larvae.

Usage of the Sea Urchin Hemicentrotus pulcherrimus Database, HpBase

HpBase ( http://cell-innovation.nig.ac.jp/Hpul/ ) is a database that provides genome and transcriptome resources of the sea urchin Hemicentrotus pulcherrimus. In addition to downloading the bulk data, several analysis tools for resource use are available: gene search, homology search, and genome browsing. HpBase also discloses the protocols for biological experiments using H. pulcherrimus that have been accumulated so far. Therefore, HpBase can assist efficient use of genome resources for researchers from various fields-evolutionary, developmental, and cell biology. In this chapter we present an overview and usage of tools in HpBase.

Advancing computer-aided drug discovery (CADD) by big data and data-driven machine learning modeling

Advancing a new drug to market requires substantial investments in time as well as financial resources. Crucial bioactivities for drug candidates, including their efficacy, pharmacokinetics (PK), and adverse effects, need to be investigated during drug development. With advancements in chemical synthesis and biological screening technologies over the past decade, a large amount of biological data points for millions of small molecules have been generated and are stored in various databases. These accumulated data, combined with new machine learning (ML) approaches, such as deep learning, have shown great potential to provide insights into relevant chemical structures to predict in vitro, in vivo, and clinical outcomes, thereby advancing drug discovery and development in the big data era.

De novo assembly of transcriptome dataset from leaves of Dryobalanops aromatica (Syn. Dryobalanops sumatrensis) seedlings grown in two contrasting potting media

Objectives

Efforts to restore tropical peat swamp forests in Indonesia face huge challenges of potential failures due to socio-economic factors and ecological dynamics attributed to lack of knowledge on the adaptive mechanisms of potential tree species such as Kapur (Dryobalanops aromatica C.F.Gaertn Syn. Dryobalanops sumatrensis J.F. Gmelin A.J.G.H Kostermans). This species is a multi-purpose tree that, commonly grows in mineral soils, but also in peat swamp as previously reported, which raised a fundamental question regarding the molecular mechanism of this adaptation. Therefore, a dataset was created aiming to detect candidates of adaptive genes in D. aromatica seedlings, cultivated in two contrasting potting media, namely mineral soil and peat media, based on RNA Sequencing Transcriptome Analysis.

Data description

The RNA transcriptome data of D. aromatica’s seedlings derived from young leaves of three one-year-old seedlings, raised in each dry mineral soil media and peat media, were generated by using Illumina HiSeq 4000 platform in NovogenAIT, Singapore. The acquired data, as the first transcriptome dataset for D. aromatica, is of a great importance in understanding molecular mechanism and responses of the involved genes of D. aromatica to the contrasting, growing potting media conditions that could also be useful to generate molecular markers.

The Dynamics, Causes, and Impacts of Mammalian Evolutionary Rates Revealed by the Analyses of Capybara Draft Genome Sequences

Capybara (Hydrochoerus hydrochaeri) is the largest species among the extant rodents. The draft genome of capybara was sequenced with the estimated genome size of 2.6 Gb. Although capybara is about 60 times larger than guinea pig, comparative analyses revealed that the neutral evolutionary rates of the two species were not substantially different. However, analyses of 39 mammalian genomes revealed very heterogeneous evolutionary rates. The highest evolutionary rate, 8.5 times higher than the human rate, was found in the Cricetidae-Muridae common ancestor after the divergence of Spalacidae. Muridae, the family with the highest number of species among mammals, emerged after the rate acceleration. Factors responsible for the evolutionary rate heterogeneity were investigated through correlations between the evolutionary rate and longevity, gestation length, litter frequency, litter size, body weight, generation interval, age at maturity, and taxonomic order. The regression analysis of these factors showed that the model with three factors (taxonomic order, generation interval, and litter size) had the highest predictive power (R2 = 0.74). These three factors determine the number of meiosis per unit time. We also conducted transcriptome analysis and found that the evolutionary rate dynamics affects the evolution of gene expression patterns.

CSI NGS Portal: An Online Platform for Automated NGS Data Analysis and Sharing

Next-generation sequencing (NGS) has been a widely-used technology in biomedical research for understanding the role of molecular genetics of cells in health and disease. A variety of computational tools have been developed to analyse the vastly growing NGS data, which often require bioinformatics skills, tedious work and a significant amount of time. To facilitate data processing steps minding the gap between biologists and bioinformaticians, we developed CSI NGS Portal, an online platform which gathers established bioinformatics pipelines to provide fully automated NGS data analysis and sharing in a user-friendly website. The portal currently provides 16 standard pipelines for analysing data from DNA, RNA, smallRNA, ChIP, RIP, 4C, SHAPE, circRNA, eCLIP, Bisulfite and scRNA sequencing, and is flexible to expand with new pipelines. The users can upload raw data in FASTQ format and submit jobs in a few clicks, and the results will be self-accessible via the portal to view/download/share in real-time. The output can be readily used as the final report or as input for other tools depending on the pipeline. Overall, CSI NGS Portal helps researchers rapidly analyse their NGS data and share results with colleagues without the aid of a bioinformatician. The portal is freely available at: https://csibioinfo.nus.edu.sg/csingsportal.

Identification of the putative N-acetylglucosaminidase CseA associated with daughter cell separation in Tetragenococcus halophilus

The lactic acid bacterium Tetragenococcus halophilus, which is used as a starter to brew soy sauce, comprises both cluster-forming strains and dispersed strains. The cluster-forming strains are industrially useful for obtaining clear soy sauce, because the cell clusters are trapped by filter cloth when the soy sauce mash is pressed. However, the molecular mechanism underlying cell cluster formation is unknown. Whole genome sequence analysis and subsequent target sequence analysis revealed that the cluster-forming strains commonly have functional defects in N-acetylglucosaminidase CseA, a peptidoglycan hydrolase. CseA is a multimodular protein that harbors a GH73 domain and six peptidoglycan-binding LysM domains. Recombinant CseA hydrolyzed peptidoglycan and promoted cell separation. Functional analysis of truncated CseA derivatives revealed that the LysM domains play an important role in efficient peptidoglycan degradation and cell separation. Taken together, the results of this study identify CseA as a factor that greatly affects the cluster formation in T. halophilus.

Bioinformatics Research Methodology of Non-coding RNAs in Cardiovascular Diseases

The transcriptional complexity generated by the human genomic output is within the core of cell and organ physiology, but also could be in the origin of pathologies. In cardiovascular diseases, the role of specific families of RNA transcripts belonging to the group of the non-coding RNAs started to be unveiled in the last two decades. The knowledge of the functional rules and roles of non-coding RNAs in the context of cardiovascular diseases is an important factor to derive new diagnostic methods, but also to design targeted therapeutic strategies. The characterization and analysis of ncRNA function requires a deep knowledge of the regulatory mechanism of these RNA species that often relies on intricated interaction networks. The use of specific bioinformatic tools to interrogate biological data and to derive functional implications is particularly relevant and needs to be extended to the general practice of translational researchers. This chapter briefly summarizes the bioinformatic tools and strategies that could be used for the characterization and functional analysis of non-coding RNAs, with special emphasis in their applications to the cardiovascular field.

Luciferase gene of a Caribbean fireworm (Syllidae) from Puerto Rico

The fireworms Odontosyllis spp. are globally distributed and well-known for their characteristic and fascinating mating behavior, with secreted mucus emitting bluish-green light. However, knowledge about the molecules involved in the light emission are still scarce. The fireworms are believed to emit light with a luciferin-luciferase reaction, but biochemical evidence of the luciferase is established for only one species living in Japan and no information is available for its luciferin structure. In this study, we identified a luciferase gene from a related Puerto Rican fireworm. We identified eight luciferase-like genes in this Puerto Rican fireworm, finding amino acid identities between Japanese and Puerto Rican luciferase-like genes to be less than 60%. We confirmed cross reactivity of extracts of the Japanese fireworm luciferin with a recombinant Puerto Rican luciferase (PR1). The emission spectrum of recombinant PR1 was similar to the crude extract of the native luciferase, suggesting that PR1 is a functional luciferase of this Puerto Rican fireworm. Our results indicate that the molecular mechanism of luminescence is widely conserved among fireworms.

Phylogenomic analyses and divergence time estimation of Elateroidea (Coleoptera) based on RNA-Seq data

Bioluminescence, the emission of visible light in a living organism, is an intriguing phenomenon observed in different species and environments. In terrestrial organisms, the bioluminescence is observed mainly in beetles of the Elateroidea superfamily (Coleoptera). Several phylogenetic studies have been used different strategies to propose a scenario for the origin and evolution of the bioluminescence within this group, however some of them showed incongruences, mainly about the relationship of the bioluminescent families. In order to increase the number of molecular markers available for Elateroidea species and to propose a more accurate phylogeny, with high supported topology, we employed the Next-Generation Sequencing (NGS) methodology to perform the RNA-Seq analysis of luminescent (Elateridae, Phengodidae, Rhagophthalmidae, and Lampyridae) and non-luminescent (Cantharidae) species of Neotropical beetles. We used the RNA-Seq data to construct a calibrated phylogeny of Elateroidea superfamily using a large number of nuclear molecular markers. The results indicate Lampyridae and Phengodidae/Rhagophthalmidae as sister-groups, suggesting that the bioluminescence evolved later in Elateridae than other families (Lampyridae, Phengodidae, and Rhagophthalmidae), and indicating the Upper Cretaceous as the period for the main diversification of Elateroidea bioluminescent species.

Experimental evolution: its principles and applications in developing stress-tolerant yeasts

Stress tolerance and resistance in industrial yeast strains are important attributes for cost-effective bioprocessing. The source of stress-tolerant yeasts ranges from extremophilic environments to laboratory engineered strains. However, industrial stress-tolerant yeasts are very rare in nature as the natural environment forces them to evolve traits that optimize survival and reproduction and not the ability to withstand harsh habitat-irrelevant industrial conditions. Experimental evolution is a frequent method used to uncover the mechanisms of evolution and microbial adaption towards environmental stresses. It optimizes biological systems by means of adaptation to environmental stresses and thus has immense power of development of robust stress-tolerant yeasts. This mini-review briefly outlines the basics and implications of evolution experiments and their applications to industrial biotechnology. This work is meant to serve as an introduction to those new to the field of experimental evolution, and as a guide to biologists working in the field of yeast stress response. Future perspectives of experimental evolution for potential biotechnological applications have also been elucidated.

Molecular basis of wax-based color change and UV reflection in dragonflies

Many animals change their body color for visual signaling and environmental adaptation. Some dragonflies show wax-based color change and ultraviolet (UV) reflection, but the biochemical properties underlying the phenomena are totally unknown. Here we investigated the UV-reflective abdominal wax of dragonflies, thereby identifying very long-chain methyl ketones and aldehydes as unique and major wax components. Little wax was detected on young adults, but dense wax secretion was found mainly on the dorsal abdomen of mature males of Orthetrum albistylum and O. melania, and pruinose wax secretion was identified on the ventral abdomen of mature females of O. albistylum and Sympetrum darwinianum. Comparative transcriptomics demonstrated drastic upregulation of the ELOVL17 gene, a member of the fatty acid elongase gene family, whose expression reflected the distribution of very long-chain methyl ketones. Synthetic 2-pentacosanone, the major component of dragonfly’s wax, spontaneously formed light-scattering scale-like fine structures with strong UV reflection, suggesting its potential utility for biomimetics.

Humans have often looked to nature for answers to problems. Living things has evolved for millions of years to deal with life’s challenges, and so engineers and inventors faced with similar challenges can also take inspiration from the natural world. Several plants and animals, for instance, reflect ultraviolet light. This ability may protect them from some of the damaging effects of sunlight; materials with similar properties would have a range of uses, including as coatings on windows that protect our homes and furniture or as cosmetics that protect ourselves in the same way.

Some dragonflies – including the white-tailed skimmer, which is particularly common in Japan – are partly coated with a wax that reflects both ultraviolet and visible light. These insects can also see ultraviolet light, which means it is likely that they also use the reflective wax to send visual signals to one another. However, the biochemistry of this wax and the genes involved in its production remained poorly understood.

Futahashi et al. have now found that the dragonfly wax consists mostly of very long-chain molecules known as methyl ketones and aldehydes; neither of which are a common components of other waxes. The wax was found in distinct patches on the bodies of adults; these patches were colored white with a hint of blue, while the rest of the dragonfly was mostly brown. Looking at gene activity in different parts of the dragonflies showed that a gene called ELOVL17 is much more active in the wax-coated areas. This gene encodes an enzyme that makes long-chain molecules, and its activity closely matched the distribution of the especially long-chain methyl ketones on the dragonflies’ surface.

Futahashi et al. then synthesized the major component of the surface wax – specifically, a chemical called 2-pentacosanone – in the laboratory, and saw that it spontaneously formed fine, scale-like structures that strongly reflected ultraviolet light. Further work is now needed to explore the potential applications of this bio-inspired wax, and to understand exactly what the dragonflies use it for in the wild.