The draft genome sequence of the Japanese rhinoceros beetle Trypoxylus dichotomus septentrionalis towards an understanding of horn formation.

The Japanese rhinoceros beetle Trypoxylus dichotomus is a giant beetle with distinctive exaggerated horns present on the head and prothoracic regions of the male. T. dichotomus has been used as a research model in various fields such as evolutionary developmental biology, ecology, ethology, biomimetics, and drug discovery. In this study, de novo assembly of 615 Mb, representing 80% of the genome estimated by flow cytometry, was obtained using the 10 × Chromium platform. The scaffold N50 length of the genome assembly was 8.02 Mb, with repetitive elements predicted to comprise 49.5% of the assembly. In total, 23,987 protein-coding genes were predicted in the genome. In addition, de novo assembly of the mitochondrial genome yielded a contig of 20,217 bp. We also analyzed the transcriptome by generating 16 RNA-seq libraries from a variety of tissues of both sexes and developmental stages, which allowed us to identify 13 co-expressed gene modules. We focused on the genes related to horn formation and obtained new insights into the evolution of the gene repertoire and sexual dimorphism as exemplified by the sex-specific splicing pattern of the doublesex gene. This genomic information will be an excellent resource for further functional and evolutionary analyses, including the evolutionary origin and genetic regulation of beetle horns and the molecular mechanisms underlying sexual dimorphism.

The genome of Lyophyllum shimeji provides insight into the initial evolution of ectomycorrhizal fungal genomes.

Mycorrhizae are one of the most fundamental symbioses between plants and fungi, with ectomycorrhizae being the most widespread in boreal forest ecosystems. Ectomycorrhizal fungi are hypothesized to have evolved convergently from saprotrophic ancestors in several fungal clades, especially members of the subdivision Agaricomycotina. Studies on fungal genomes have identified several typical characteristics of mycorrhizal fungi, such as genome size expansion and decreases in plant cell-wall degrading enzymes (PCWDEs). However, genomic changes concerning the evolutionary transition to the ectomycorrhizal lifestyle are largely unknown. In this study, we sequenced the genome of Lyophyllum shimeji, an ectomycorrhizal fungus that is phylogenetically related to saprotrophic species and retains some saprotroph-like traits. We found that the genome of Ly. shimeji strain AT787 lacks both incremental increases in genome size and reduced numbers of PCWDEs. Our findings suggest that the previously reported common genomic traits of mycorrhizal fungi are not essential for the ectomycorrhizal lifestyle, but are a result of abolishing saprotrophic activity. Since Ly. shimeji is commercially consumed as an edible mushroom, the newly available genomic information may also impact research designed to enhance the cultivation of this mushroom.

Antiemetic effects of baclofen in a shrew model of postoperative nausea and vomiting: Whole-transcriptome analysis in the nucleus of the solitary tract.

AIMS: The molecular genetic mechanisms underlying postoperative nausea and vomiting (PONV) in the brain have not been fully elucidated. This study aimed to determine the changes in whole transcriptome in the nucleus of the solitary tract (NTS) in an animal model of PONV, to screen a drug candidate and to elucidate the molecular genetic mechanisms of PONV development. METHODS: Twenty-one female musk shrews were assigned into three groups: the Surgery group (shrew PONV model, n = 9), the Sham group (n = 6), and the Naïve group (n = 6). In behavioral studies, the main outcome was the number of emetic episodes. In genetic experiments, changes in the transcriptome in the NTS were measured. In a separate study, 12 shrews were used to verify the candidate mechanism underlying PONV. RESULTS: A median of six emetic episodes occurred in both the Sham and Surgery groups. Whole-transcriptome analysis indicated the inhibition of the GABAB receptor-mediated signaling pathway in the PONV model. Baclofen (GABAB receptor agonist) administration eliminated emetic behaviors in the shrew PONV model. CONCLUSIONS: Our findings suggest that the GABAB receptor-mediated signaling pathway is involved in emesis and that baclofen may be a novel therapeutic or prophylactic agent for PONV.

Corrigendum: Networking and Specificity-Changing DNA Methyltransferases in Helicobacter pylori.

[This corrects the article DOI: 10.3389/fmicb.2020.01628.].

Ethylene signaling mediates host invasion by parasitic plants.

Parasitic plants form a specialized organ, a haustorium, to invade host tissues and acquire water and nutrients. To understand the molecular mechanism of haustorium development, we performed a forward genetics screening to isolate mutants exhibiting haustorial defects in the model parasitic plant Phtheirospermum japonicum. We isolated two mutants that show prolonged and sometimes aberrant meristematic activity in the haustorium apex, resulting in severe defects on host invasion. Whole-genome sequencing revealed that the two mutants respectively have point mutations in homologs of ETHYLENE RESPONSE 1 (ETR1) and ETHYLENE INSENSITIVE 2 (EIN2), signaling components in response to the gaseous phytohormone ethylene. Application of the ethylene signaling inhibitors also caused similar haustorial defects, indicating that ethylene signaling regulates cell proliferation and differentiation of parasite cells. Genetic disruption of host ethylene production also perturbs parasite invasion. We propose that parasitic plants use ethylene as a signal to invade host roots.

Networking and Specificity-Changing DNA Methyltransferases in Helicobacter pylori.

Epigenetic DNA base methylation plays important roles in gene expression regulation. We here describe a gene expression regulation network consisting of many DNA methyltransferases each frequently changing its target sequence-specificity. Our object Helicobacter pylori, a bacterium responsible for most incidence of stomach cancer, carries a large and variable repertoire of sequence-specific DNA methyltransferases. By creating a dozen of single-gene knockout strains for the methyltransferases, we revealed that they form a network controlling methylome, transcriptome and adaptive phenotype sets. The methyltransferases interact with each other in a hierarchical way, sometimes regulated positively by one methyltransferase but negatively with another. Motility, oxidative stress tolerance and DNA damage repair are likewise regulated by multiple methyltransferases. Their regulation sometimes involves translation start and stop codons suggesting coupling of methylation, transcription and translation. The methyltransferases frequently change their sequence-specificity through gene conversion of their target recognition domain and switch their target sets to remodel the network. The emerging picture of a metamorphosing gene regulation network, or firework, consisting of epigenetic systems ever-changing their specificity in search for adaptation, provides a new paradigm in understanding global gene regulation and adaptive evolution.

Genomic insights of body plan transitions from bilateral to pentameral symmetry in Echinoderms.

Echinoderms are an exceptional group of bilaterians that develop pentameral adult symmetry from a bilaterally symmetric larva. However, the genetic basis in evolution and development of this unique transformation remains to be clarified. Here we report newly sequenced genomes, developmental transcriptomes, and proteomes of diverse echinoderms including the green sea urchin (L. variegatus), a sea cucumber (A. japonicus), and with particular emphasis on a sister group of the earliest-diverged echinoderms, the feather star (A. japonica). We learned that the last common ancestor of echinoderms retained a well-organized Hox cluster reminiscent of the hemichordate, and had gene sets involved in endoskeleton development. Further, unlike in other animal groups, the most conserved developmental stages were not at the body plan establishing phase, and genes normally involved in bilaterality appear to function in pentameric axis development. These results enhance our understanding of the divergence of protostomes and deuterostomes almost 500 Mya.

Gene expression analysis of three homeobox genes throughout early and late development of a feather star Anneissia japonica.

Crinoids are considered as the most basal extant echinoderms. They retain aboral nervous system with a nerve center, which has been degraded in the eleutherozoan echinoderms. To investigate the evolution of patterning of the nervous systems in crinoids, we examined temporal and spatial expression patterns of three neural patterning-related homeobox genes, six3, pax6, and otx, throughout the development of a feather star Anneissia japonica. These genes were involved in the patterning of endomesodermal tissues instead of the ectodermal neural tissues in the early planktonic stages. In the stages after larval attachment, the expression of these genes was mainly observed in the podia and the oral nervous systems instead of the aboral nerve center. Our results indicate the involvement of these three genes in the formation of oral nervous system in the common ancestor of the echinoderms and suggest that the aboral nerve center is not evolutionally related to the brain of other bilaterians.

Genomes of the Venus Flytrap and Close Relatives Unveil the Roots of Plant Carnivory.

Most plants grow and develop by taking up nutrients from the soil while continuously under threat from foraging animals. Carnivorous plants have turned the tables by capturing and consuming nutrient-rich animal prey, enabling them to thrive in nutrient-poor soil. To better understand the evolution of botanical carnivory, we compared the draft genome of the Venus flytrap (Dionaea muscipula) with that of its aquatic sister, the waterwheel plant Aldrovanda vesiculosa, and the sundew Drosera spatulata. We identified an early whole-genome duplication in the family as source for carnivory-associated genes. Recruitment of genes to the trap from the root especially was a major mechanism in the evolution of carnivory, supported by family-specific duplications. Still, these genomes belong to the gene poorest land plants sequenced thus far, suggesting reduction of selective pressure on different processes, including non-carnivorous nutrient acquisition. Our results show how non-carnivorous plants evolved into the most skillful green hunters on the planet.

Construction of JRG (Japanese reference genome) with single-molecule real-time sequencing.

In recent genome analyses, population-specific reference panels have indicated important. However, reference panels based on short-read sequencing data do not sufficiently cover long insertions. Therefore, the nature of long insertions has not been well documented. Here, we assembled a Japanese genome using single-molecule real-time sequencing data and characterized insertions found in the assembled genome. We identified 3691 insertions ranging from 100 bps to ~10,000 bps in the assembled genome relative to the international reference sequence (GRCh38). To validate and characterize these insertions, we mapped short-reads from 1070 Japanese individuals and 728 individuals from eight other populations to insertions integrated into GRCh38. With this result, we constructed JRGv1 (Japanese Reference Genome version 1) by integrating the 903 verified insertions, totaling 1,086,173 bases, shared by at least two Japanese individuals into GRCh38. We also constructed decoyJRGv1 by concatenating 3559 verified insertions, totaling 2,536,870 bases, shared by at least two Japanese individuals or by six other assemblies. This assembly improved the alignment ratio by 0.4% on average. These results demonstrate the importance of refining the reference assembly and creating a population-specific reference genome. JRGv1 and decoyJRGv1 are available at the JRG website.

Construction of full-length Japanese reference panel of class I HLA genes with single-molecule, real-time sequencing.

Human leukocyte antigen (HLA) is a gene complex known for its exceptional diversity across populations, importance in organ and blood stem cell transplantation, and associations of specific alleles with various diseases. We constructed a Japanese reference panel of class I HLA genes (ToMMo HLA panel), comprising a distinct set of HLA-A, HLA-B, HLA-C, and HLA-H alleles, by single-molecule, real-time (SMRT) sequencing of 208 individuals included in the 1070 whole-genome Japanese reference panel (1KJPN). For high-quality allele reconstruction, we developed a novel pipeline, Primer-Separation Assembly and Refinement Pipeline (PSARP), in which the SMRT sequencing and additional short-read data were used. The panel consisted of 139 alleles, which were all extended from known IPD-IMGT/HLA sequences, contained 40 with novel variants, and captured more than 96.5% of allelic diversity in 1KJPN. These newly available sequences would be important resources for research and clinical applications including high-resolution HLA typing, genetic association studies, and analyzes of cis-regulatory elements.

Genome of the pitcher plant Cephalotus reveals genetic changes associated with carnivory.

Carnivorous plants exploit animals as a nutritional source and have inspired long-standing questions about the origin and evolution of carnivory-related traits. To investigate the molecular bases of carnivory, we sequenced the genome of the heterophyllous pitcher plant Cephalotus follicularis, in which we succeeded in regulating the developmental switch between carnivorous and non-carnivorous leaves. Transcriptome comparison of the two leaf types and gene repertoire analysis identified genetic changes associated with prey attraction, capture, digestion and nutrient absorption. Analysis of digestive fluid proteins from C. follicularis and three other carnivorous plants with independent carnivorous origins revealed repeated co-options of stress-responsive protein lineages coupled with convergent amino acid substitutions to acquire digestive physiology. These results imply constraints on the available routes to evolve plant carnivory.

Comparative analysis of the brain transcriptome in a hyper-aggressive fruit fly, Drosophila prolongata.

Aggressive behavior is observed in many animals, but its intensity differs between species. In a model animal of genetics, Drosophila melanogaster, genetic basis of aggressive behavior has been studied intensively, including transcriptome analyses to identify genes whose expression level was associated with intra-species variation in aggressiveness. However, whether these genes are also involved in the evolution of aggressiveness among different species has not been examined. In this study, we performed de novo transcriptome analysis in the brain of Drosophila prolongata to identify genes associated with the evolution of aggressiveness. Males of D. prolongata were hyper-aggressive compared with closely related species. Comparison of the brain transcriptomes identified 21 differentially expressed genes in males of D. prolongata. They did not overlap with the list of aggression-related genes identified in D. melanogaster, suggesting that genes involved in the evolution of aggressiveness were independent of those associated with the intra-species variation in aggressiveness in Drosophila. Although females of D. prolongata were not aggressive as the males, expression levels of the 21 genes identified in this study were more similar between sexes than between species.

Experimental Approach Reveals the Role of alx1 in the Evolution of the Echinoderm Larval Skeleton.

Over the course of evolution, the acquisition of novel structures has ultimately led to wide variation in morphology among extant multicellular organisms. Thus, the origins of genetic systems for new morphological structures are a subject of great interest in evolutionary biology. The larval skeleton is a novel structure acquired in some echinoderm lineages via the activation of the adult skeletogenic machinery. Previously, VEGF signaling was suggested to have played an important role in the acquisition of the larval skeleton. In the present study, we compared expression patterns of Alx genes among echinoderm classes to further explore the factors involved in the acquisition of a larval skeleton. We found that the alx1 gene, originally described as crucial for sea urchin skeletogenesis, may have also played an essential role in the evolution of the larval skeleton. Unlike those echinoderms that have a larval skeleton, we found that alx1 of starfish was barely expressed in early larvae that have no skeleton. When alx1 overexpression was induced via injection of alx1 mRNA into starfish eggs, the expression patterns of certain genes, including those possibly involved in skeletogenesis, were altered. This suggested that a portion of the skeletogenic program was induced solely by alx1. However, we observed no obvious external phenotype or skeleton. We concluded that alx1 was necessary but not sufficient for the acquisition of the larval skeleton, which, in fact, requires several genetic events. Based on these results, we discuss how the larval expression of alx1 contributed to the acquisition of the larval skeleton in the putative ancestral lineage of echinoderms.

Comparative genome sequencing reveals genomic signature of extreme desiccation tolerance in the anhydrobiotic midge.

Anhydrobiosis represents an extreme example of tolerance adaptation to water loss, where an organism can survive in an ametabolic state until water returns. Here we report the first comparative analysis examining the genomic background of extreme desiccation tolerance, which is exclusively found in larvae of the only anhydrobiotic insect, Polypedilum vanderplanki. We compare the genomes of P. vanderplanki and a congeneric desiccation-sensitive midge P. nubifer. We determine that the genome of the anhydrobiotic species specifically contains clusters of multi-copy genes with products that act as molecular shields. In addition, the genome possesses several groups of genes with high similarity to known protective proteins. However, these genes are located in distinct paralogous clusters in the genome apart from the classical orthologues of the corresponding genes shared by both chironomids and other insects. The transcripts of these clustered paralogues contribute to a large majority of the mRNA pool in the desiccating larvae and most likely define successful anhydrobiosis. Comparison of expression patterns of orthologues between two chironomid species provides evidence for the existence of desiccation-specific gene expression systems in P. vanderplanki.

Whole-Genome Sequence of Burkholderia sp. Strain RPE67, a Bacterial Gut Symbiont of the Bean Bug Riptortus pedestris.

Burkholderia sp. strain RPE67 is a bacterial symbiont isolated from a field-collected bean bug, Riptortus pedestris. To understand the genetic basis of the insect-microbe symbiosis, we performed whole-genome sequencing of the Burkholderia strain, revealing an 8.69-Mb genome consisting of three chromosomes and three plasmids.

Methylome diversification through changes in DNA methyltransferase sequence specificity.

Epigenetic modifications such as DNA methylation have large effects on gene expression and genome maintenance. Helicobacter pylori, a human gastric pathogen, has a large number of DNA methyltransferase genes, with different strains having unique repertoires. Previous genome comparisons suggested that these methyltransferases often change DNA sequence specificity through domain movement--the movement between and within genes of coding sequences of target recognition domains. Using single-molecule real-time sequencing technology, which detects N6-methyladenines and N4-methylcytosines with single-base resolution, we studied methylated DNA sites throughout the H. pylori genome for several closely related strains. Overall, the methylome was highly variable among closely related strains. Hypermethylated regions were found, for example, in rpoB gene for RNA polymerase. We identified DNA sequence motifs for methylation and then assigned each of them to a specific homology group of the target recognition domains in the specificity-determining genes for Type I and other restriction-modification systems. These results supported proposed mechanisms for sequence-specificity changes in DNA methyltransferases. Knocking out one of the Type I specificity genes led to transcriptome changes, which suggested its role in gene expression. These results are consistent with the concept of evolution driven by DNA methylation, in which changes in the methylome lead to changes in the transcriptome and potentially to changes in phenotype, providing targets for natural or artificial selection.

Complete Genome Sequence of Burkholderia sp. Strain RPE64, Bacterial Symbiont of the Bean Bug Riptortus pedestris.

We isolated Burkholderia symbiont strain RPE64 from the bean bug Riptortus pedestris. Analysis of the complete 6.96-Mb genome, which consists of three chromosomes and two plasmids, will facilitate further understanding of insect-microbe symbiosis and the development of pest-control technologies.

Glycomics of a novel type-2 N-acetyllactosamine-specific lectin purified from the feather star, Oxycomanthus japonicus (Pelmatozoa: Crinoidea).

A lectin - designated OXYL for the purposes of this study that strongly recognizes complex-type oligosaccharides of serum glycoproteins - was purified from a crinoid, the feather star Oxycomanthus japonicus, the most basal group among extant echinoderms. OXYL was purified through a combination of anion-exchange and affinity chromatography using Q-sepharose and fetuin-sepharose gel, respectively. Lectin was determined to be a 14-kDa polypeptide by sodium dodecyl sulphate-polyacrylamide gel electrophoresis under reducing conditions. However, 14-kDa and 28-kDa bands appeared in the same proportion under non-reducing conditions. Gel permeation chromatography showed a 54-kDa peak, suggesting that lectin consists of four 14-kDa subunits. Divalent cations were not indicated, and stable haemagglutination activity was demonstrated at pH 4-12 and temperatures below 60°C. Surface plasmon resonance analysis of OXYL against fetuin showed k(ass) and k(diss) values of 1.4×10(-6)M(-1)s(-1) and 3.1×10(-3)s(-1), respectively, indicating that it has a strong binding affinity to the glycoprotein as lectin. Frontal affinity chromatography using 25 types of prydylamine-conjugated glycans indicated that OXYL specifically recognizes multi-antennary complex-type oligosaccharides containing type-2 N-acetyllactosamines (Galß1-4GlcNAc) if α2-3-linked sialic acid is linked at the non-reducing terminal. However, type-1 N-acetyllactosamine (Galß1-3GlcNAc) chains and α2-6-linked sialic acids were never recognized by OXYL. This profiling study showed that OXYL essentially recognizes ß1-4-linkage at C-1 position and free OH group at C-6 position of Gal in addition to the conservation of N-acetyl groups at C-2 position and free OH groups at C-3 position of GlcNAc in N-acetyllactosamine. This is the first report on glycomics on a lectin purified from an echinoderm belonging to the subphylum Pelmatozoa.