Published papers

Background: Mycorrhizal symbiosis is one of the most fundamental types of mutualistic plant-microbe interaction. Among the many classes of mycorrhizae, the arbuscular mycorrhizae have the most general symbiotic style and the longest history. However, the genomes of arbuscular mycorrhizal (AM) fungi are not well characterized due to difficulties in cultivation and genetic analysis. In this study, we sequenced the genome of the AM fungus Rhizophagus clarus HR1, compared the sequence with the genome sequence of the model species R. irregularis, and checked for missing genes that encode enzymes in metabolic pathways related to their obligate biotrophy. Results: In the genome of R. clarus, we confirmed the absence of cytosolic fatty acid synthase (FAS), whereas all mitochondrial FAS components were present. A KEGG pathway map identified the absence of genes encoding enzymes for several other metabolic pathways in the two AM fungi, including thiamine biosynthesis and the conversion of vitamin B6 derivatives. We also found that a large proportion of the genes encoding glucose-producing polysaccharide hydrolases, that are present even in ectomycorrhizal fungi, also appear to be absent in AM fungi. Conclusions: In this study, we found several new genes that are absent from the genomes of AM fungi in addition to the genes previously identified as missing. Missing genes for enzymes in primary metabolic pathways imply that AM fungi may have a higher dependency on host plants than other biotrophic fungi. These missing metabolic pathways provide a genetic basis to explore the physiological characteristics and auxotrophy of AM fungi.

Abstract Arbuscular mycorrhizal fungi (AMF) establish symbiotic relationships with most land plants, mainly for the purpose of nutrient exchange. Many studies have revealed the regulation of processes in AMF, such as nutrient absorption from soil, metabolism and exchange with host plants, and the genes involved. However, the spatial regulation of the genes within the structures comprising each developmental stage is not well understood. Here, we demonstrate the structure-specific transcriptome of the model AMF species, Rhizophagus irregularis. We performed an ultra-low input RNA-seq analysis, SMART-seq2, comparing five extraradical structures, germ tubes, runner hyphae, branched absorbing structures (BAS), immature spores and mature spores. In addition, we reanalyzed the recently reported RNA-seq data comparing intraradical mycelium and arbuscule. Our analyses captured the distinct features of each structure and revealed the structure-specific expression patterns of genes related to nutrient transport and metabolism. Of note, the transcriptional profiles suggest distinct functions of BAS in nutrient absorption. These findings provide a comprehensive dataset to advance our understanding of the transcriptional dynamics of fungal nutrition in this symbiotic system.

The signaling lymphocyte activation molecule (SLAM) is a receptor for morbilliviruses. To understand the recent host range expansion of canine distemper virus (CDV) in carnivores, we determined the nucleotide sequences of SLAMs of various carnivores and generated three-dimensional homology SLAM models. Thirty-four amino acid residues were found for the candidates binding to CDV on the interface of the carnivore SLAMs. SLAM of the domestic dog (Canis lupus familiaris) were similar to those of other members of the suborder Caniformia, indicating that the animals in this group have similar sensitivity to dog CDV. However, they were different at nine positions from those of felids. Among the nine residues, four of domestic cat (Felis catus) SLAM (72, 76, 82, and 129) and three of lion (Panthera leo persica) SLAM (72, 82, and 129) were associated with charge alterations, suggesting that the felid interfaces have lower affinities to dog CDV. Only the residue at 76 was different between domestic cat and lion SLAM interfaces. The domestic cat SLAM had threonine at 76, whereas the lion SLAM had arginine, a positively charged residue like that of the dog SLAM. The cat SLAM with threonine is likely to have lower affinity to CDV-H and to confer higher resistance against dog CDV. Thus, the four residues (72, 76, 82, and 129) on carnivore SLAMs are important for the determination of affinity and sensitivity with CDV. Additionally, the CDV-H protein of felid strains had a substitution of histidine for tyrosine at 549 of dog CDV-H and may have higher affinity to lion SLAM. Three-dimensional model construction is a new risk assessment method of morbillivirus infectivity. Because the method is applicable to animals that have no information about virus infection, it is especially useful for morbillivirus risk assessment and wildlife conservation.

Beetles, representing the majority of the insect species diversity, are characterized by thick and hard cuticle, which plays important roles for their environmental adaptation and underpins their inordinate diversity and prosperity. Here, we report a bacterial endosymbiont extremely specialized for sustaining beetle's cuticle formation. Many weevils are associated with gamma-proteobacterial endosymbiont lineage Nardonella, whose evolutionary origin is estimated as older than 100 million years, but its functional aspect has been elusive. Sequencing of Nardonella genomes from diverse weevils unveiled drastic size reduction to 0.2 Mb, in which minimal complete gene sets for bacterial replication, transcription, and translation were present but almost all of the other metabolic pathway genes were missing. Notably, the only metabolic pathway retained in the Nardonella genomes was the tyrosine synthesis pathway, identifying tyrosine provisioning as Nardonella's sole biological role. Weevils are armored with hard cuticle, tyrosine is the principal precursor for cuticle formation, and experimental suppression of Nardonella resulted in emergence of reddish and soft weevils with low tyrosine titer, confirming the importance of Nardonella-mediated tyrosine production for host's cuticle formation and hardening. Notably, Nardonella's tyrosine synthesis pathway was incomplete, lacking the final step transaminase gene. RNA sequencing identified host's aminotransferase genes up-regulated in the bacteriome. RNA interference targeting the aminotransferase genes induced reddish and soft weevils with low tyrosine titer, verifying host's final step regulation of the tyrosine synthesis pathway. Our finding highlights an impressively intimate and focused aspect of the host-symbiont metabolic integrity via streamlined evolution for a single biological function of ecological relevance.

On low intertidal and shallow subtidal shores on the west coast of Okinawa, Japan, we investigated the trophic associations of sacoglossan opisthobranchs associated with Bryopsidalean green algae. During 11 short research visits (55 days total) from 2002 to 2008, we recorded almost 500 specimens of I I species. These sacoglossans include a new record for Japan (Caliphylla A. Costa, 1867), a recent record for Japan (Placida daguilarensis Jensen, 1990), two undescribed species (Placida Trinchese, 1876 and Elysia Risso, 1818), one unnamed (but well-described) species (Placida sp. sensu Baba 1986), and six other Indo-Pacific species. Not only did we record more sacoglossan species but also we found higher slug abundances than other colleagues in Okinawa or the Indo-Pacific region. Quantitative population attributes and feeding preferences are described for these sacoglossans. In contrast to temperate geographic regions, several of these Japanese sacoglossans specialized on a single algal genus rather than two or more genera in different families. This specificity is consistent with narrower host-plant associations in high-diversity communities; yet monophagy has not yet been demonstrated in this guild of Okinawan sacoglossans. Given the broad geographic ranges, restricted host ranges, often predictable populations, and high frequency of life cycles with planktotrophic larvae, western Pacific subtropical sacoglossans should be considered "predictably rare" (sensu Rabinowitz 1981) rather than at "special risk" (sensu Clark 1994).

Maternal investment for offspring's growth and survival is widespread among diverse organisms [1-3]. Vertical symbiont transmission via maternal passage is also pivotal for offspring's growth and survival in many organisms [4-6]. Hence, it is expected that vertical symbiont transmission may coevolve with various organismal traits concerning maternal investment in offspring. Here we report a novel phenotypic syndrome entailing morphological, histological, behavioral, and ecological specializations for maternal investment and vertical symbiont transmission in stinkbugs of the family Urostylididae [7-9]. Adult females develop huge ovaries exaggerated for polysaccharide excretion, possess novel ovipositor-associated organs for vertical transmission of a bacterial symbiont ("Candidatus Tachikawaea gelatinosa"), and lay eggs covered with voluminous symbiont-supplemented jelly. Newborns hatch in midwinter, feed solely on the jelly, acquire the symbiont, and grow during winter. In spring, the insects start feeding on plant sap, wherein the symbiont localizes to a specialized midgut region and supplies essential amino acids deficient in the host's diet. The reduced symbiont genome and host-symbiont co-speciation indicate their obligate association over evolutionary time. Experimental deprivation of the jelly results in nymphal mortality, whereas restoration of the jelly leads to recovered nymphal growth, confirming that the jelly supports nymphal growth in winter. Chemical analyses demonstrate that the galactan-based jelly contains a sufficient quantity of amino acids to sustain nymphal growth to the third instar. The versatile biological roles of the symbiont-containing egg-covering jelly highlight intricate evolutionary interactions between maternal resource investment and vertical symbiont transmission, which are commonly important for offspring's growth, survival, and ecological adaptation.

Arbuscular mycorrhizal (AM) symbiosis is a mutually beneficial interaction between fungi and land plants and promotes global phosphate cycling in terrestrial ecosystems. AM fungi are recognised as obligate symbionts that require root colonisation to complete a life cycle involving the production of propagules, asexual spores. Recently, it has been shown that Rhizophagus irregularis can produce infection-competent secondary spores asymbiotically by adding a fatty acid, palmitoleic acid. Furthermore, asymbiotic growth can be supported using myristate as a carbon and energy source for their asymbiotic growth to increase fungal biomass. However, the spore production and the ability of these spores to colonise host roots were still limited compared to the co-culture of the fungus with plant roots. Here we show that a combination of two plant hormones, strigolactone and jasmonate, induces the production of a large number of infection-competent spores in asymbiotic cultures of Rhizophagus clarus HR1 in the presence of myristate and organic nitrogen. Inoculation of asymbiotically-generated spores promoted the growth of host plants, as observed for spores produced by symbiotic culture system. Our findings provide a foundation for the elucidation of hormonal control of the fungal life cycle and the development of inoculum production schemes.

Elysia chlorotica, a sacoglossan sea slug found off the East Coast of the United States, is well-known for its ability to sequester chloroplasts from its algal prey and survive by photosynthesis for up to 12 months in the absence of food supply. Here we present a draft genome assembly of E. chlorotica that was generated using a hybrid assembly strategy with Illumina short reads and PacBio long reads. The genome assembly comprised 9,989 scaffolds, with a total length of 557 Mb and a scaffold N50 of 442 kb. BUSCO assessment indicated that 93.3% of the expected metazoan genes were completely present in the genome assembly. Annotation of the E. chlorotica genome assembly identified 176 Mb (32.6%) of repetitive sequences and a total of 24,980 protein-coding genes. We anticipate that the annotated draft genome assembly of the E. chlorotica sea slug will promote the investigation of sacoglossan genetics, evolution, and particularly, the genetic signatures accounting for the long-term functioning of algal chloroplasts in an animal.

© 2018, The Author(s). Arbuscular mycorrhizal fungus (AMF) species are some of the most widespread symbionts of land plants. Our much improved reference genome assembly of a model AMF, Rhizophagus irregularis DAOM-181602 (total contigs = 210), facilitated a discovery of repetitive elements with unusual characteristics. R. irregularis has only ten or 11 copies of complete 45S rDNAs, whereas the general eukaryotic genome has tens to thousands of rDNA copies. R. irregularis rDNAs are highly heterogeneous and lack a tandem repeat structure. These findings provide evidence for the hypothesis that rDNA heterogeneity depends on the lack of tandem repeat structures. RNA-Seq analysis confirmed that all rDNA variants are actively transcribed. Observed rDNA/rRNA polymorphisms may modulate translation by using different ribosomes depending on biotic and abiotic interactions. The non-tandem repeat structure and intragenomic heterogeneity of AMF rDNA/rRNA may facilitate successful adaptation to various environmental conditions, increasing host compatibility of these symbiotic fungi.

Members of the genus Rickettsiella are bacterial pathogens of insects and other arthropods. Recently, a novel facultative endosymbiont, “Candidatus Rickettsiella viridis,” was described in the pea aphid Acyrthosiphon pisum, whose infection causes a striking host phenotype: red and green genetic color morphs exist in aphid populations, and upon infection with the symbiont, red aphids become green due to increased production of green polycyclic quinone pigments. Here we determined the complete genome sequence of the symbiont. The 1.6-Mb circular genome, harboring some 1,400 protein-coding genes, was similar to the genome of entomopathogenic Rickettsiella grylli (1.6 Mb) but was smaller than the genomes of phylogenetically allied human pathogens Coxiella burnetii (2.0 Mb) and Legionella pneumophila (3.4 Mb). The symbiont’s metabolic pathways exhibited little complementarity to those of the coexisting primary symbiont Buchnera aphidicola, reflecting the facultative nature of the symbiont. The symbiont genome harbored neither polyketide synthase genes nor the evolutionarily allied fatty acid synthase genes that are suspected to catalyze the polycyclic quinone synthesis, indicating that the green pigments are produced not by the symbiont but by the host aphid. The symbiont genome retained many type IV secretion system genes and presumable effector protein genes, whose homologues in L. pneumophila were reported to modulate a variety of the host’s cellular processes for facilitating infection and virulence. These results suggest the possibility that the symbiont is involved in the green pigment production by affecting the host’s metabolism using the secretion machineries for delivering the effector molecules into the host cells. IMPORTANCE Insect body color is relevant to a variety of biological aspects such as species recognition, sexual selection, mimicry, aposematism, and crypsis. Hence, the bacterial endosymbiont “Candidatus Rickettsiella viridis,” which alters aphid body color from red to green, is of ecological interest, given that different predators preferentially exploit either red-or green-colored aphids. Here we determined the complete 1.6-Mb genome of the symbiont and uncovered that, although the red-green color transition was ascribed to upregulated production of green polycyclic quinone pigments, the symbiont genome harbored few genes involved in the polycyclic quinone biosynthesis. Meanwhile, the symbiont genome contained type IV secretion system genes and presumable effector protein genes, whose homologues modulate eukaryotic cellular processes for facilitating infection and virulence in the pathogen Legionella pneumophila. We propose the hypothesis that the symbiont may upregu-late the host’s production of polycyclic quinone pigments via cooption of secretion machineries and effector molecules for pathogenicity.

The sea slug Plakobranchus ocellatus (Sacoglossa, Gastropoda) retains photosynthetically active chloroplasts from ingested algae (functional kleptoplasts) in the epithelial cells of its digestive gland for up to 10 months. While its feeding behavior has not been observed in natural habitats, two hypotheses have been proposed: 1) adult P. ocellatus uses kleptoplasts to obtain photosynthates and nutritionally behaves as a photoautotroph without replenishing the kleptoplasts; or 2) it behaves as a mixotroph (photoautotroph and herbivorous consumer) and replenishes kleptoplasts continually or periodically. To address the question of which hypothesis is more likely, we examined the source algae for kleptoplasts and temporal changes in kleptoplast composition and nutritional contribution. By characterizing the temporal diversity of P. ocellatus kleptoplasts using rbcL sequences, we found that P. ocellatus harvests kleptoplasts from at least 8 different siphonous green algal species, that kleptoplasts from more than one species are present in each individual sea slug, and that the kleptoplast composition differs temporally. These results suggest that wild P. ocellatus often feed on multiple species of siphonous algae from which they continually obtain fresh chloroplasts. By estimating the trophic position of wild and starved P. ocellatus using the stable nitrogen isotopic composition of amino acids, we showed that despite the abundance of kleptoplasts, their photosynthates do not contribute greatly to the nutrition of wild P. ocellatus, but that kleptoplast photosynthates form a significant source of nutrition for starved sea slugs. The herbivorous nature of wild P. ocellatus is consistent with insights from molecular analyses indicating that kleptoplasts are frequently replenished from ingested algae, leading to the conclusion that natural populations of P. ocellatus do not rely on photosynthesis but mainly on the digestion of ingested algae.

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.

Background: Two Calyptogena clam intracellular obligate symbionts, Ca. Vesicomyosocius okutanii (Vok; C. okutanii symbiont) and Ca. Ruthia magnifica (Rma; C. magnifica symbiont), have small genomes (1.02 and 1.16 Mb, respectively) with low G+C contents (31.6% and 34.0%, respectively) and are thought to be in an ongoing stage of reductive genome evolution (RGE). They lack recA and some genes for DNA repair, including mutY. The loss of recA and mutY is thought to contribute to the stabilization of their genome architectures and GC bias, respectively. To understand how these genes were lost from the symbiont genomes, we surveyed these genes in the genomes from 10 other Calyptogena clam symbionts using the polymerase chain reaction (PCR). Results: Phylogenetic trees reconstructed using concatenated 16S and 23S rRNA gene sequences showed that the symbionts formed two clades, clade I (symbionts of C. kawamurai, C. laubieri, C. kilmeri, C. okutanii and C. soyoae) and clade II (those of C. pacifica, C. fausta, C. nautilei, C. stearnsii, C. magnifica, C. fossajaponica and C. phaseoliformis). recA was detected by PCR with consensus primers for recA in the symbiont of C. phaseoliformis. A detailed homology search revealed a remnant recA in the Rma genome. Using PCR with a newly designed primer set, intact recA or its remnant was detected in clade II symbionts. In clade I symbionts, the recA coding region was found to be mostly deleted. In the Rma genome, a pseudogene of mutY was found. Using PCR with newly designed primer sets, mutY was not found in clade I symbionts but was found in clade II symbionts. The G+C content of 16S and 23S rRNA genes in symbionts lacking mutY was significantly lower than in those with mutY. Conclusions: The extant Calyptogena clam symbionts in clade II were shown to have recA and mutY or their remnants, while those in clade I did not. The present results indicate that the extant symbionts are losing these genes in RGE, and that the loss of mutY contributed to the GC bias of the genomes during their evolution.

Sacoglossan sea slugs retain the chloroplasts from food algae in their cells (kleptoplasty) and obtain the photosynthetic products, but the capability of kleptoplasty differs among species. One evolutionary hypothesis for kleptoplasty is that the family Volvatellidae is the ancestral sacoglossan lineage in which kleptoplasty arose, but this is based on ambiguous phylogenetic relationships. To test this hypothesis, we analyzed sacoglossan phylogeny for concatenated sequences (approximately 4000 bp) of nuclear (18S rRNA) and mitochondrial (CO1, trnV. and 16S rRNA) genes of 18 sacoglossan species, using neighbor-joining, maximum parsimony, maximum likelihood, and Bayesian methods. The resultant phylogenetic trees clearly revealed monophyly of the order Sacoglossa (including the debatable genus Cylindrobulla), which contained two sister clades: the shelled sacoglossan group (Oxynoacea + Cylindrobulla) and the non-shelled group (Plakobranchacea). The family Volvatellidae was in a derived position within the Oxynoacea. The most parsimonious rearrangement of the character states of kleptoplasty on the phylogenetic tree suggested that (1) kleptoplasty was acquired at the basal position of the Sacoglossa (including Cylindrobulla) as non-functional kleptoplasty, (2) functional kleptoplasty was gained in the Plakobranchacea clade, and (3) non-functional kleptoplasty was lost in the Volvatellidae lineage in the Oxynoacea clade.

Bradyrhizobium sp. strain SUTN9-2 is a symbiotic and endophytic diazotrophic bacterium found in legume and rice plants and has the potential to promote growth. The present results revealed that SUTN9-2 underwent cell enlargement, increased its DNA content, and efficiently performed nitrogen fixation in response to rice extract. Some factors in rice extract induced the expression of cell cycle and nitrogen fixation genes. According to differentially expressed genes (DEGs) from the transcriptomic analysis, SUTN9-2 was affected by rice extract and the deletion of the bclA gene. The up-regulated DEGs encoding a class of oxidoreductases, which act with oxygen atoms and may have a role in controlling oxygen at an appropriate level for nitrogenase activity, followed by GroESL chaperonins are required for the function of nitrogenase. These results indicate that following its exposure to rice extract, nitrogen fixation by SUTN9-2 is induced by the collective effects of GroESL and oxidoreductases. The expression of the sensitivity to antimicrobial peptides transporter (sapDF) was also up-regulated, resulting in cell differentiation, even when bclA (sapDF) was mutated. This result implies similarities in the production of defensin-like antimicrobial peptides (DEFs) by rice and nodule-specific cysteine-rich (NCR) peptides in legume plants, which affect bacterial cell differentiation.

Some sea slugs sequester chloroplasts from algal food in their intestinal cells and photosynthesize for months. This phenomenon, kleptoplasty, poses a question of how the chloroplast retains its activity without the algal nucleus. There have been debates on the horizontal transfer of algal genes to the animal nucleus. To settle the arguments, this study reported the genome of a kleptoplastic sea slug, Plakobranchus ocellatus, and found no evidence of photosynthetic genes encoded on the nucleus. Nevertheless, it was confirmed that light illumination prolongs the life of mollusk under starvation. These data presented a paradigm that a complex adaptive trait, as typified by photosynthesis, can be transferred between eukaryotic kingdoms by a unique organelle transmission without nuclear gene transfer. Our phylogenomic analysis showed that genes for proteolysis and immunity undergo gene expansion and are up-regulated in chloroplast-enriched tissue, suggesting that these molluskan genes are involved in the phenotype acquisition without horizontal gene transfer.