Severe Bottleneck Impacted the Genomic Structure of Egg-Eating Cichlids in Lake Victoria

Within 15,000 years, the explosive adaptive radiation of haplochromine cichlids in Lake Victoria, East Africa, generated 500 endemic species. In the 1980s, the upsurge of Nile perch, a carnivorous fish artificially introduced to the lake, drove the extinction of more than 200 endemic cichlids. The Nile perch predation particularly harmed piscivorous cichlids, including paedophages, cichlids eat eggs and fries, which is an example of the unique trophic adaptation seen in African cichlids. Here, aiming to investigate past demographic events possibly triggered by the invasion of Nile perch and the subsequent impacts on the genetic structure of cichlids, we conducted large-scale comparative genomics. We discovered evidence of recent bottleneck events in four species, including two paedophages, which began during the 1970s-1980s, and population size rebounded during the 1990s-2000s. The timing of the bottleneck corresponded to the historical records of endemic haplochromines' disappearance and later resurgence, which is likely associated with the introduction of Nile perch by commercial demand to Lake Victoria in the 1950s. Interestingly, among the four species that likely experienced bottleneck, Haplochromis sp. 'matumbi hunter,' a paedophagous cichlid, showed the most severe bottleneck signatures. The components of shared ancestry inferred by ADMIXTURE suggested a high genetic differentiation between matumbi hunter and other species. In contrast, our phylogenetic analyses highly supported the monophyly of the five paedophages, consistent with the results of previous studies. We conclude that high genetic differentiation of matumbi hunter occurred due to the loss of shared genetic components among haplochromines in Lake Victoria caused by the recent severe bottleneck.

Loss of one X and the Y chromosome changes the configuration of the X inactivation center in the genus Tokudaia

INTRODUCTION

X chromosome inactivation (XCI) is an essential mechanism for dosage compensation between females and males in mammals. In females, XCI is controlled by a complex, conserved locus termed the X inactivation center (Xic), in which the lncRNA Xist is the key regulator. However, little is known about the Xic in species with unusual sex chromosomes. The genus Tokudaia includes three rodent species endemic to Japan. Tokudaia osimensis (TOS) and Tokudaia tokunoshimensis (TTO) lost the Y chromosome (XO/XO), while Tokudaia muenninki (TMU) acquired a neo-X region by fusion of the X chromosome and an autosome (XX/XY). We compared the gene location and structure in the Xic among Tokudaia species.

METHODS

Gene structure of nine genes in Xic were predicted, and the gene location and genome sequences of Xic were compared between mouse and Tokudaia species. The expression level of gene was confirmed by TPM calculation using RNA-seq data.

RESULTS

Compared to mouse, the Xic gene order and location were conserved in Tokudaia species. However, remarkable structure changes were observed in lncRNA genes, Xist and Tsix, in the XO/XO species. In Xist, important functional repeats, B-, C-, D-, and E-repeats, were partially or completely lost due to deletions in these species. RNA-seq data showed that female-specific expression patterns of Xist and Tsix were confirmed in TMU, however not in the XO/XO species. Additionally, three deletions and one inversion were confirmed in the intergenic region between Jpx and Ftx in the XO/XO species.

CONCLUSION

Our findings indicate that even if the Xist and Tsix lncRNAs are expressed, they are incapable of producing a successful and lasting XCI in the XO/XO species. We hypothesized that the significant structure change in intergenic region of Jpx-Ftx resulted in the inability to perform the X chromosome inactivation, and, as a result, a lack of Xist expression. Our results collectively suggest that structural changes in the Xic occurred in the ancestral lineage of XO/XO species, likely due to the loss of one X chromosome and the Y chromosome and a consequence of the degradation of XCI system.

SDF-1/CXCR4 signal is involved in the induction of Primordial Germ Cell migration in a model marine fish, Japanese anchovy (Engraulis japonicus)

Primordial germ cells (PGCs) are pivotal for gonadal development and reproductive success. Though artificial induction of sterility by targeting PGCs are gaining popularity due to its advantages in fish surrogacy and biodiversity management, it is often skill and time intensive. In this study, we have focused on understanding the role of PGCs and the chemotactic SDF-1/CXCR4 signaling on gonad development of Japanese anchovy (JA, Engraulis japonicus), an upcoming marine model organism with eco-commercial values, with an aim to develop a novel, easy, and versatile gonad sterilization method. Our data showed that PGC migration related genes, i.e., sdf-1a, sdf-1b, cxcr4a, cxcr4b and vasa, are phylogenetically closer relatives of respective herring (Clupea harengus) and zebrafish (Danio rerio) homolog. Subsequently, PGC marking and live tracing experiments confirmed that PGC migration in JA initiates from 16 hours post fertilization (hpf) followed by PGC settlement in the gonadal ridge at 44 hpf. We found that overexpression of zebrafish sdf-1a mRNA in the germ cell suppresses cxcr4a and increases cxcr4b transcription at 8 hpf, dose dependently disrupts PGC migration at 24-48 hpf, induces PGC death and upregulates sdf-1b at 5 days post hatching. 48 h of immersion treatment with CXCR4 antagonist (AMD3100, Abcam) also accelerated PGC mismigration and pushed the PGC away from gonadal ridge in a dose responsive manner, and further when grown to adulthood caused germ cell less gonad formation in some individuals. Cumulatively, our data, for the first time, suggests that JA PGC migration is largely regulated by SDF1/CXCR4 signaling, and modulation of this signaling has strong potential for sterile, germ cell less gonad preparation at a mass scale. However, further in-depth analysis is pertinent to apply this methodology in marine fish species to successfully catapult Japanese anchovy into a true marine fish model.

Metabolic engineering with adaptive laboratory evolution for phenylalanine production by Corynebacterium glutamicum

The mutants resistant to a phenylalanine analog, 4-fluorophenylalanine (4FP), were obtained for metabolic engineering of Corynebacterium glutamicum for producing aromatic amino acids synthesized through the shikimate pathway by adaptive laboratory evolution. Culture experiments of the C. glutamicum strains which carry the mutations found in the open reading frame from the 4FP-resistant mutants revealed that the mutations in the open reading frames of aroG (NCgl2098), pheA (NCgl2799) and aroP (NCgl1062) encoding 3-deoxy-d-arabino-heptulosonate-7-phosphate, prephenate dehydratase, and aromatic amino acid transporter are responsible for 4FP resistance and higher concentration of aromatic amino acids in their culture supernatants in the 4FP-resistant strains. It was expected that aroG and pheA mutations would release feedback inhibition of the enzymes involved in the shikimate pathway by phenylalanine and that aroP mutations would prevent intracellular uptake of aromatic amino acids. Therefore, we conducted metabolic engineering of the C. glutamicum wild-type strain for aromatic amino acid production and found that phenylalanine production at 6.11 ± 0.08 g L-1 was achieved by overexpressing the mutant pheA and aroG genes from the 4FP-resistant mutants and deleting aroP gene. This study demonstrates that adaptive laboratory evolution is an effective way to obtain useful mutant genes related to production of target material and to establish metabolic engineering strategies.

Genome sequence and cell biological toolbox of the highly regenerative, coenocytic green feather alga Bryopsis

Green feather algae (Bryopsidales) undergo a unique life cycle in which a single cell repeatedly executes nuclear division without cytokinesis, resulting in the development of a thallus (>100 mm) with characteristic morphology called coenocyte. Bryopsis is a representative coenocytic alga that has exceptionally high regeneration ability: extruded cytoplasm aggregates rapidly in seawater, leading to the formation of protoplasts. However, the genetic basis of the unique cell biology of Bryopsis remains poorly understood. Here, we present a high-quality assembly and annotation of the nuclear genome of Bryopsis sp. (90.7 Mbp, 27 contigs, N50 = 6.7 Mbp, 14 034 protein-coding genes). Comparative genomic analyses indicate that the genes encoding BPL-1/Bryohealin, the aggregation-promoting lectin, are heavily duplicated in Bryopsis, whereas homologous genes are absent in other ulvophyceans, suggesting the basis of regeneration capability of Bryopsis. Bryopsis sp. possesses >30 kinesins but only a single myosin, which differs from other green algae that have multiple types of myosin genes. Consistent with this biased motor toolkit, we observed that the bidirectional motility of chloroplasts in the cytoplasm was dependent on microtubules but not actin in Bryopsis sp. Most genes required for cytokinesis in plants are present in Bryopsis, including those in the SNARE or kinesin superfamily. Nevertheless, a kinesin crucial for cytokinesis initiation in plants (NACK/Kinesin-7II) is hardly expressed in the coenocytic part of the thallus, possibly underlying the lack of cytokinesis in this portion. The present genome sequence lays the foundation for experimental biology in coenocytic macroalgae.

Modeling the SDF-1/CXCR4 protein using advanced artificial intelligence and antagonist screening for Japanese anchovy

SDF-1/CXCR4 chemokine signaling are indispensable for cell migration, especially the Primordial Germ Cell (PGC) migration towards the gonadal ridge during early development. We earlier found that this signaling is largely conserved in the Japanese anchovy (Engraulis japonicus, EJ), and a mere treatment of CXCR4 antagonist, AMD3100, leads to germ cell depletion and thereafter gonad sterilization. However, the effect of AMD3100 was limited. So, in this research, we scouted for CXCR4 antagonist with higher potency by employing advanced artificial intelligence deep learning-based computer simulations. Three potential candidates, AMD3465, WZ811, and LY2510924, were selected and in vivo validation was conducted using Japanese anchovy embryos. We found that seven transmembrane motif of EJ CXCR4a and EJ CXCR4b were extremely similar with human homolog while the CXCR4 chemokine receptor N terminal (PF12109, essential for SDF-1 binding) was missing in EJ CXCR4b. 3D protein analysis and cavity search predicted the cavity in EJ CXCR4a to be five times larger (6,307 Å³) than that in EJ CXCR4b (1,241 Å³). Docking analysis demonstrated lower binding energy of AMD3100 and AMD3465 to EJ CXCR4a (Vina score -9.6) and EJ CXCR4b (Vina score -8.8), respectively. Furthermore, we observed significant PGC mismigration in microinjected AMD3465 treated groups at 10, 100 and 1 × 105 nM concentration in 48 h post fertilized embryos. The other three antagonists showed various degrees of PGC dispersion, but no significant effect compared to their solvent control at tested concentrations was observed. Cumulatively, our results suggests that AMD3645 might be a better candidate for abnormal PGC migration in Japanese anchovy and warrants further investigation.

Chromosomal-level assembly of Tokudaia osimensis, Tokudaia tokunoshimensis, and Tokudaia muenninki genomes

Herein, we present the first high-quality long-read-based chromosome-level genome assemblies and gene annotations of the genomes of three endangered Tokudaia species: Tokudaia osimensis, Tokudaia tokunoshimensis, and Tokudaia muenninki. These species, which are endemic to different islands of the Ryukyu Islands, Japan, exhibited unique karyotypes and sex chromosomal characteristics. The genome assemblies generated using PacBio, Illumina, and Hi-C sequence data consisted of 13 (corresponded to 12 autosomes and one X chromosome), 23 (corresponded to 22 autosomes and one X chromosome), and 23 (corresponded to 21 autosomes and the neo- and ancestral X regions) chromosome-level scaffolds that contained 2,445, 2,477, and 2,661 Mbp of sequence data, respectively. Annotations of protein-coding genes were performed using RNA-Seq-based, homology-based, and Ab initio methods. BUSCO completeness values for every species exceeded 96% for genomes and 98% for genes. These data can be an important resource for contributing to our understanding of species genomes resulting from allopatric speciation and provide insights into mammalian sex-determination mechanisms and sex chromosome evolution.

Complete genomes of mutualistic bacterial co-symbionts "Candidatus Sulcia muelleri" and "Candidatus Nasuia deltocephalinicola" of the rice green leafhopper Nephotettix cincticeps

The genomes of obligate bacterial co-symbionts of the green rice leafhopper Nephotettix cincticeps, which is notorious as an agricultural pest, were determined. The streamlined genomes of "Candidatus Sulcia muelleri" and "Candidatus Nasuia deltocephalinicola" exhibited complementary metabolic pathways for synthesizing essential nutrients that contribute to host adaptation.

Chromosomal-level Genome Assembly of the Coffee Bee Hawk Moth Reveals the Evolution of Chromosomes and the Molecular Basis of Distinct Phenotypes

Cephonodes hylas, the coffee bee hawk moth is a hawk moth species with unique characteristics, such as larvae feeding on gardenia, overcoming the toxicity of its iridoid glycosides, diurnal adults, and transparent wings. Although C. hylas is a fascinating model for molecular biological research, genome sequence analysis-based genetic approaches to elucidate these peculiarities have not yet been undertaken. We successfully achieved de novo genome assembly at the chromosome level of C. hylas comparable to the Lepidoptera model organism, silkworm. Additionally, 16,854 protein-coding genes were annotated, and the constructed genome sequence and annotated genes were of the highest quality BUSCO completion compared to closely related species. Comparative genome analysis revealed the process of chromosomal evolution from the Bombycoidea ancestral (n = 31) genome and changes in turnover at the chromosome level associated with chromosomal fusion events, such as the rate of repetitive sequence insertion. These analyses were only possible because the genome was constructed at the chromosome level. Additionally, increased the nonsynonymous/synonymous rate (dN/dS) ratios were observed in multiple photoreceptor-related genes that were strongly associated with the acquisition of diurnal activity. Furthermore, tandemly duplicated expanded genes containing many digestive and other enzymes and larval midgut-specific expression were also confirmed. These genes may be involved in the metabolism of genipin, a toxin found in gardenias. Using the genome sequence of C. hylas determined at the chromosome level, we have successfully identified new insights into the chromosomal evolution of Bombycoidea, as well as the relationship between the genome sequence and its characteristic traits.

Bioinformatic and fine-scale chromosomal mapping reveal the nature and evolution of eliminated chromosomes in the Japanese hagfish, Eptatretus burgeri, through analysis of repetitive DNA families

In the Japanese hagfish, Eptatretus burgeri, approximately 21% of the genomic DNA in germ cells (2n = 52) consists of 16 chromosomes (eliminated [E]-chromosomes) that are eliminated from presumptive somatic cells (2n = 36). To uncover the eliminated genome (E-genome), we have identified 16 eliminated repetitive DNA families from eight hagfish species, with 11 of these repeats being selectively amplified in the germline genome of E. burgeri. Furthermore, we have demonstrated that six of these sequences, namely EEEb1-6, are exclusively localized on all 16 E-chromosomes. This has led to the hypothesis that the eight pairs of E-chromosomes are derived from one pair of ancestral chromosomes via multiple duplication events over a prolonged evolutionary period. NGS analysis has recently facilitated the re-assembly of two distinct draft genomes of E. burgeri, derived from the testis and liver. This advancement allows for the prediction of not only nonrepetitive eliminated sequences but also over 100 repetitive and eliminated sequences, accomplished through K-mer-based analysis. In this study, we report four novel eliminated repetitive DNA sequences (designated as EEEb7-10) and confirm the relative chromosomal localization of all eliminated repeats (EEEb1-10) by fluorescence in situ hybridization (FISH). With the exception of EEEb10, all sequences were exclusively detected on EEEb1-positive chromosomes. Surprisingly, EEEb10 was detected as an intense signal on EEEb1-positive chromosomes and as a scattered signal on other chromosomes in germ cells. The study further divided the eight pairs of E-chromosomes into six groups based on the signal distribution of each DNA family, and fiber-FISH experiments showed that the EEEb2-10 family was dispersed in the EEEb1-positive extended chromatin fiber. These findings provide new insights into the mechanisms underlying chromosome elimination and the evolution of E-chromosomes, supporting our previous hypothesis.

GINGER: an integrated method for high-accuracy prediction of gene structure in higher eukaryotes at the gene and exon level

The prediction of gene structure within the genome sequence is the starting point of genome analysis, and its accuracy has a significant impact on the quality of subsequent analyses. Gene structure prediction is roughly divided into RNA-Seq-based methods, ab initio-based methods, homology-based methods, and the integration of individual prediction methods. Integrated methods are mainstream in recent genome projects because they improve prediction accuracy by combining or taking the best individual prediction findings; however, adequate prediction accuracy for eukaryotic species has not yet been achieved. Therefore, we developed an integrated tool, GINGER, that solves various issues related to gene structure prediction in higher eukaryotes. By handling artefacts in alignments of RNA and protein sequences, reconstructing gene structures via dynamic programming with appropriately weighted and scored exon/intron/intergenic regions, and applying different prediction processes and filtering criteria to multi-exon and single-exon genes, we achieved a significant improvement in accuracy compared to the existing integration methods. The feature of GINGER is its high prediction accuracy at the gene and exon levels, which is pronounced for species with more complex gene architectures. GINGER is implemented using Nextflow, which allows for the efficient and effective use of computing resources.

Characterization and development of a plastid genome base editor, ptpTALECD

The modification of photosynthesis-related genes in plastid genomes may improve crop yields. Recently, we reported that a plastid-targeting base editor named ptpTALECD, in which a cytidine deaminase DddA functions as the catalytic domain, can homoplasmically substitute a targeted C to T in plastid genomes of Arabidopsis thaliana. However, some target Cs were not substituted. In addition, although ptpTALECD could substitute Cs on the 3' side of T and A, it was unclear whether it could also substitute Cs on the 3' side of G and C. In this study, we identified the preferential positions of the substituted Cs in ptpTALECD-targeting sequences in the Arabidopsis plastid genome. We also found that ptpTALECD could substitute Cs on the 3' side of all four bases in plastid genomes of Arabidopsis. More recently, a base editor containing an improved version of DddA (DddA11) was reported to substitute Cs more efficiently, and to substitute Cs on the 3' side of more varieties of bases in human mitochondrial genomes than a base editor containing DddA. Here, we also show that ptpTALECD_v2, in which a modified version of DddA11 functions as the catalytic domain, more frequently substituted Cs than ptpTALECD in the Arabidopsis plastid genome. We also found that ptpTALECD_v2 tended to substitute Cs at more positions than ptpTALECD. Our results reveal that ptpTALECD can cause a greater variety of codon changes and amino acid substitutions than previously thought, and that ptpTALECD and ptpTALECD_v2 are useful tools for the targeted base editing of plastid genomes.

GreenHill: a de novo chromosome-level scaffolding and phasing tool using Hi-C

Chromosome-level haplotype-resolved genome assembly is an important resource in molecular biology. However, current de novo haplotype assemblers require parental data or reference genomes and often fail to provide chromosome-level results. We present GreenHill, a novel scaffolding and phasing tool that considers various assemblers' contigs as input to reconstruct chromosome-level haplotypes using Hi-C without parental or reference data. Its unique functions include new error correction based on Hi-C contacts and the simultaneous use of Hi-C and long reads. Benchmarks reveal that GreenHill outperforms other approaches in contiguity and phasing accuracy, and the majority of chromosome arms are entirely phased.

Haplotype-resolved chromosomal-level assembly of wasabi (Eutrema japonicum) genome

In Japan, wasabi (Eutrema japonicum) is an important traditional condiment, and is recognized as an endemic species. In the present study, we generated a chromosome-level and haplotype-resolved reference genome for E. japonicum using PacBio CLR (continuous long reads), Illumina, and Hi-C sequencing data. The genome consists of 28 chromosomes that contain 1,512.1 Mb of sequence data, with a scaffold N50 length of 55.67 Mb. We also reported the subgenome and haplotype assignment of the 28 chromosomes by read-mapping and phylogenic analysis. Three validation methods (Benchmarking Universal Single-Copy Orthologs, Merqury, and Inspector) indicated that our obtained genome sequences were a high-quality and high-completeness genome assembly. Comparison of genome assemblies from previously published genomes showed that our obtained genome was of higher quality. Therefore, our genome will serve as a valuable genetic resource for both chemical ecology and evolution research of the genera Eutrema and Brassicaceae, as well as for wasabi breeding.

Alteration of a Shiga toxin-encoding phage associated with a change in toxin production level and disease severity in Escherichia coli

Among the nine clades of Shiga toxin (Stx)-producing Escherichia coli O157:H7, clade 8 is thought to be highly pathogenic, as it causes severe disease more often than other clades. Two subclades have been proposed, but there are conflicting reports on intersubclade differences in Stx2 levels, although Stx2 production is a risk factor for severe disease development. The global population structure of clade 8 has also yet to be fully elucidated. Here, we present genome analyses of a global clade 8 strain set (n=510), including 147 Japanese strains sequenced in this study. The complete genome sequences of 18 of the 147 strains were determined to perform detailed clade-wide genome analyses together with 17 publicly available closed genomes. Intraclade variations in Stx2 production level and disease severity were also re-evaluated within the phylogenetic context. Based on phylogenomic analysis, clade 8 was divided into four lineages corresponding to the previously proposed SNP genotypes (SGs): SG8_30, SG8_31A, SG8_31B and SG8_32. SG8_30 and the common ancestor of the other SGs were first separated, with SG8_31A and SG8_31B emerging from the latter and SG8_32 emerging from SG8_31B. Comparison of 35 closed genomes revealed the overall structure of chromosomes and pO157 virulence plasmids and the prophage contents to be well conserved. However, Stx2a phages exhibit notable genomic diversity, even though all are integrated into the argW locus, indicating that subtype changes in Stx2a phage occurred from the γ subtype to its variant (γ_v1) in SG8_31A and from γ to δ in SG8_31B and SG8_32 via replacement of parts or almost entire phage genomes, respectively. We further show that SG8_30 strains (all carrying γ Stx2a phages) produce significantly higher levels of Stx2 and cause severe disease more frequently than SG8_32 strains (all carrying δ Stx2a phages). Clear conclusions on SG8_31A and SG8_31B cannot be made due to the small number of strains available, but as SG8_31A (carrying γ_v1 Stx2a phages) contains strains that produce much more Stx2 than SG8_30 strains, attention should also be paid to this SG.

Gene Recruitments and Dismissals in the Argonaut Genome Provide Insights into Pelagic Lifestyle Adaptation and Shell-like Eggcase Reacquisition

The paper nautilus or greater argonaut, Argonauta argo, is a species of octopods which is characterized by its pelagic lifestyle and by the presence of a protective spiral-shaped shell-like eggcase in females. To reveal the genomic background of how the species adapted to the pelagic lifestyle and acquired its shell-like eggcase, we sequenced the draft genome of the species. The genome size was 1.1 Gb, which is the smallest among the cephalopods known to date, with the top 215 scaffolds (average length 5,064,479 bp) covering 81% (1.09 Gb) of the total assembly. A total of 26,433 protein-coding genes were predicted from 16,802 assembled scaffolds. From these, we identified nearly intact HOX, Parahox, Wnt clusters, and some gene clusters that could probably be related to the pelagic lifestyle, such as reflectin, tyrosinase, and opsin. The gene models also revealed several homologous genes related to calcified shell formation in Conchiferan mollusks, such as Pif-like, SOD, and TRX. Interestingly, comparative genomics analysis revealed that the homologous genes for such genes were also found in the genome of the shell-less octopus, as well as Nautilus, which has a true outer shell. Therefore, the draft genome sequence of Arg. argo presented here has helped us to gain further insights into the genetic background of the dynamic recruitment and dismissal of genes to form an important, converging extended phenotypic structure such as the shell and the shell-like eggcase. Additionally, it allows us to explore the evolution of from benthic to pelagic lifestyles in cephalopods and octopods.

Genomic architecture and functional unit of mimicry supergene in female limited Batesian mimic Papilio butterflies

It has long been suggested that dimorphic female-limited Batesian mimicry of two closely related Papilio butterflies, Papilio memnon and Papilio polytes, is controlled by supergenes. Whole-genome sequencing, genome-wide association studies and functional analyses have recently identified mimicry supergenes, including the doublesex (dsx) gene. Although supergenes of both the species are composed of highly divergent regions between mimetic and non-mimetic alleles and are located at the same chromosomal locus, they show critical differences in genomic architecture, particularly with or without an inversion: P. polytes has an inversion, but P. memnon does not. This review introduces and compares the detailed genomic structure of mimicry supergenes in two Papilio species, including gene composition, repetitive sequence composition, breakpoint/boundary site structure, chromosomal inversion and linkage disequilibrium. Expression patterns and functional analyses of the respective genes within or flanking the supergene suggest that dsx and other genes are involved in mimetic traits. In addition, structural comparison of the corresponding region for the mimicry supergene among further Papilio species suggests three scenarios for the evolution of the mimicry supergene between the two Papilio species. The structural features revealed in the Papilio mimicry supergene provide insight into the formation, maintenance and evolution of supergenes. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Evolutionary history of sexual differentiation mechanism in insects

Alternative splicing underpins functional diversity in proteins and the complexity and diversity of eukaryotes. An example is the doublesex gene, the key transcriptional factor in arthropod sexual differentiation. doublesex is controlled by sex-specific splicing and promotes both male and female differentiation in holometabolan insects, whereas in hemimetabolan species, doublesex has sex-specific isoforms but is not required for female differentiation. How doublesex evolved to be essential for female development remains largely unknown. Here, we investigate ancestral states of doublesex using Thermobia domestica belonging to Zygentoma, the sister group of Pterygota, that is, winged insects. We find that, in T. domestica, doublesex expresses sex-specific isoforms but is only necessary for male differentiation of sexual morphology. This result supports the hypothesis that doublesex initially promoted male differentiation during insect evolution. However, T. domestica doublesex has a short female-specific region and upregulates the expression of vitellogenin homologs in females, suggesting that doublesex may already play some role in female morphogenesis of the common ancestor of Pterygota. Reconstruction of the ancestral sequence and prediction of protein structures show that the female-specific isoform of doublesex has an extended C-terminal disordered region in holometabolan insects but not in nonholometabolan species. We propose that doublesex acquired its function in female morphogenesis through a change in the protein motif structure rather than the emergence of the female-specific exon.

Whole-genome sequencing analysis and protocol for RNA interference of the endoparasitoid wasp Asobara japonica

Asobara japonica is an endoparasitic wasp that parasitizes Drosophila flies. It synthesizes various toxic components in the venom gland and injects them into host larvae during oviposition. To identify and characterize these toxic components for enabling parasitism, we performed the whole-genome sequencing (WGS) and devised a protocol for RNA interference (RNAi) with A. japonica. Because it has a parthenogenetic lineage due to Wolbachia infection, we generated a clonal strain from a single wasp to obtain highly homogenous genomic DNA. The WGS analysis revealed that the estimated genome size was 322 Mb with a heterozygosity of 0.132%. We also performed RNA-seq analyses for gene annotation. Based on the qualified WGS platform, we cloned ebony-Aj, which encodes the enzyme N-β-alanyl dopamine synthetase, which is involved in melanin production. The microinjection of double-stranded RNA (dsRNA) targeting ebony-Aj led to body colour changes in adult wasps, phenocopying ebony-Dm mutants. Furthermore, we identified putative venom genes as a target of RNAi, confirming that dsRNA injection-based RNAi specifically suppressed the expression of the target gene in wasp adults. Taken together, our results provide a powerful genetic toolkit for studying the molecular mechanisms of parasitism.

Global population structure of the Serratia marcescens complex and identification of hospital-adapted lineages in the complex

Serratia marcescens is an important nosocomial pathogen causing various opportunistic infections, such as urinary tract infections, bacteremia and sometimes even hospital outbreaks. The recent emergence and spread of multidrug-resistant (MDR) strains further pose serious threats to global public health. This bacterium is also ubiquitously found in natural environments, but the genomic differences between clinical and environmental isolates are not clear, including those between S. marcescens and its close relatives. In this study, we performed a large-scale genome analysis of S. marcescens and closely related species (referred to as the 'S. marcescens complex'), including more than 200 clinical and environmental strains newly sequenced here. Our analysis revealed their phylogenetic relationships and complex global population structure, comprising 14 clades, which were defined based on whole-genome average nucleotide identity. Clades 10, 11, 12 and 13 corresponded to S. nematodiphila, S. marcescens sensu stricto, S. ureilytica and S. surfactantfaciens, respectively. Several clades exhibited distinct genome sizes and GC contents and a negative correlation of these genomic parameters was observed in each clade, which was associated with the acquisition of mobile genetic elements (MGEs), but different types of MGEs, plasmids or prophages (and other integrative elements), were found to contribute to the generation of these genomic variations. Importantly, clades 1 and 2 mostly comprised clinical or hospital environment isolates and accumulated a wide range of antimicrobial resistance genes, including various extended-spectrum β-lactamase and carbapenemase genes, and fluoroquinolone target site mutations, leading to a high proportion of MDR strains. This finding suggests that clades 1 and 2 represent hospital-adapted lineages in the S. marcescens complex although their potential virulence is currently unknown. These data provide an important genomic basis for reconsidering the classification of this group of bacteria and reveal novel insights into their evolution, biology and differential importance in clinical settings.

Mining RNA-seq data reveals the massive regulon of GcvB small RNA and its physiological significance in maintaining amino acid homeostasis in Escherichia coli

Bacterial small RNAs regulate the expression of multiple genes through imperfect base-pairing with target mRNAs mediated by RNA chaperone proteins such as Hfq. GcvB is the master sRNA regulator of amino acid metabolism and transport in a wide range of Gram-negative bacteria. Recently, independent RNA-seq approaches identified a plethora of transcripts interacting with GcvB in Escherichia coli. In this study, the compilation of RIL-seq, CLASH, and MAPS data sets allowed us to identify GcvB targets with high accuracy. We validated 21 new GcvB targets repressed at the posttranscriptional level, raising the number of direct targets to >50 genes in E. coli. Among its multiple seed sequences, GcvB utilizes either R1 or R3 to regulate most of these targets. Furthermore, we demonstrated that both R1 and R3 seed sequences are required to fully repress the expression of gdhA, cstA, and sucC genes. In contrast, the ilvLXGMEDA polycistronic mRNA is targeted by GcvB through at least four individual binding sites in the mRNA. Finally, we revealed that GcvB is involved in the susceptibility of peptidase-deficient E. coli strain (Δpeps) to Ala-Gln dipeptide by regulating both Dpp dipeptide importer and YdeE dipeptide exporter via R1 and R3 seed sequences, respectively.

MetaPlatanus: a metagenome assembler that combines long-range sequence links and species-specific features

De novo metagenome assembly is effective in assembling multiple draft genomes, including those of uncultured organisms. However, heterogeneity in the metagenome hinders assembly and introduces interspecies misassembly deleterious for downstream analysis. For this purpose, we developed a hybrid metagenome assembler, MetaPlatanus. First, as a characteristic function, it assembles the basic contigs from accurate short reads and then iteratively utilizes long-range sequence links, species-specific sequence compositions, and coverage depth. The binning information was also used to improve contiguity. Benchmarking using mock datasets consisting of known bacteria with long reads or mate pairs revealed the high contiguity MetaPlatanus with a few interspecies misassemblies. For published human gut data with nanopore reads from potable sequencers, MetaPlatanus assembled many biologically important elements, such as coding genes, gene clusters, viral sequences, and over-half bacterial genomes. In the benchmark with published human saliva data with high-throughput nanopore reads, the superiority of MetaPlatanus was considerably more evident. We found that some high-abundance bacterial genomes were assembled only by MetaPlatanus as near-complete. Furthermore, MetaPlatanus can circumvent the limitations of highly fragmented assemblies and frequent interspecies misassembles obtained by the other tools. Overall, the study demonstrates that MetaPlatanus could be an effective approach for exploring large-scale structures in metagenomes.

The global population structure and evolutionary history of the acquisition of major virulence factor-encoding genetic elements in Shiga toxin-producing Escherichia coli O121:H19

Shiga toxin (Stx)-producing Escherichia coli (STEC) are foodborne pathogens causing serious diseases, such as haemorrhagic colitis and haemolytic uraemic syndrome. Although O157:H7 STEC strains have been the most prevalent, incidences of STEC infections by several other serotypes have recently increased. O121:H19 STEC is one of these major non-O157 STECs, but systematic whole genome sequence (WGS) analyses have not yet been conducted on this STEC. Here, we performed a global WGS analysis of 638 O121:H19 strains, including 143 sequenced in this study, and a detailed comparison of 11 complete genomes, including four obtained in this study. By serotype-wide WGS analysis, we found that O121:H19 strains were divided into four lineages, including major and second major lineages (named L1 and L3, respectively), and that the locus of enterocyte effacement (LEE) encoding a type III secretion system (T3SS) was acquired by the common ancestor of O121:H19. Analyses of 11 complete genomes belonging to L1 or L3 revealed remarkable interlineage differences in the prophage pool and prophage-encoded T3SS effector repertoire, independent acquisition of virulence plasmids by the two lineages, and high conservation in the prophage repertoire, including that for Stx2a phages in lineage L1. Further sequence determination of complete Stx2a phage genomes of 49 strains confirmed that Stx2a phages in lineage L1 are highly conserved short-tailed phages, while those in lineage L3 are long-tailed lambda-like phages with notable genomic diversity, suggesting that an Stx2a phage was acquired by the common ancestor of L1 and has been stably maintained. Consistent with these genomic features of Stx2a phages, most lineage L1 strains produced much higher levels of Stx2a than lineage L3 strains. Altogether, this study provides a global phylogenetic overview of O121:H19 STEC and shows the interlineage genomic differences and the highly conserved genomic features of the major lineage within this serotype of STEC.

Analysis of Sex Chromosome Evolution in the Clade Palaeognathae from Phased Genome Assembly

Birds in the clade Palaeognathae, excluding Tinamiformes, have morphologically conserved karyotypes and less differentiated ZW sex chromosomes compared with those of other birds. In particular, the sex chromosomes of the ostrich and emu have exceptionally large recombining pseudoautosomal regions (PARs), whereas non-PARs are classified into two strata according to the date of their origins: stratum 0 and stratum 1 (S1). However, the construction and analysis of the genome sequences in these regions in the clade Palaeognathae can be challenging because assembling the S1 region is difficult owing to low sequence diversity between gametologs (Z-linked and W-linked sequences). We addressed this issue by applying the Platanus-allee assembler and successfully constructed the haplotype-resolved (phased) assembly for female emu, cassowary, and ostrich using only sequence read data derived from the Illumina platform. Comparative genomic and phylogenetic analyses based on assembled Z-linked and W-linked sequences confirmed that the S1 region of emu and cassowary formed in their common ancestor. Moreover, the interspersed repetitive sequence landscapes in the S1 regions of female emu showed an expansion of younger repetitive elements in the W-linked S1 region, suggesting an interruption in homologous recombination in the S1 region. These results provide novel insights into the trajectory of sex chromosome evolution in the clade Palaeognathae and suggest that the Illumina-based phased assembly method is an effective approach for elucidating the evolutionary process underlying the transition from homomorphic to differentiated sex chromosomes.

A chromosome-level genome sequence of Chrysanthemum seticuspe, a model species for hexaploid cultivated chrysanthemum

Chrysanthemums are one of the most industrially important cut flowers worldwide. However, their segmental allopolyploidy and self-incompatibility have prevented the application of genetic analysis and modern breeding strategies. We thus developed a model strain, Gojo-0 (Chrysanthemum seticuspe), which is a diploid and self-compatible pure line. Here, we present the 3.05 Gb chromosome-level reference genome sequence, which covered 97% of the C. seticuspe genome. The genome contained more than 80% interspersed repeats, of which retrotransposons accounted for 72%. We identified recent segmental duplication and retrotransposon expansion in C. seticuspe, contributing to arelatively large genome size. Furthermore, we identified a retrotransposon family, SbdRT, which was enriched in gene-dense genome regions and had experienced a very recent transposition burst. We also demonstrated that the chromosome-level genome sequence facilitates positional cloning in C. seticuspe. The genome sequence obtained here can greatly contribute as a reference for chrysanthemum in front-line breeding including genome editing.

Elucidation of the speciation history of three sister species of crown-of-thorns starfish (Acanthaster spp.) based on genomic analysis

The crown-of-thorns starfish (COTS) is a coral predator that is widely distributed in Indo-Pacific Oceans. A previous phylogenetic study using partial mitochondrial sequences suggested that COTS had diverged into four distinct species, but a nuclear genome-based analysis to confirm this was not conducted. To address this, COTS species nuclear genome sequences were analysed here, sequencing Northern Indian Ocean (NIO) and Red Sea (RS) species genomes for the first time, followed by a comparative analysis with the Pacific Ocean (PO) species. Phylogenetic analysis and ADMIXTURE analysis revealed clear divergences between the three COTS species. Furthermore, within the PO species, the phylogenetic position of the Hawaiian sample was further away from the other Pacific-derived samples than expected based on the mitochondrial data, suggesting that it may be a PO subspecies. The pairwise sequentially Markovian coalescent model showed that the trajectories of the population size diverged by region during the Mid-Pleistocene transition when the sea-level was dramatically decreased, strongly suggesting that the three COTS species experienced allopatric speciation. Analysis of the orthologues indicated that there were remarkable genes with species-specific positive selection in the genomes of the PO and RS species, which suggested that there may be local adaptations in the COTS species.

A spinach genome assembly with remarkable completeness, and its use for rapid identification of candidate genes for agronomic traits

Spinach (Spinacia oleracea) is grown as a nutritious leafy vegetable worldwide. To accelerate spinach breeding efficiency, a high-quality reference genome sequence with great completeness and continuity is needed as a basic infrastructure. Here, we used long-read and linked-read technologies to construct a de novo spinach genome assembly, designated SOL_r1.1, which was comprised of 287 scaffolds (total size: 935.7 Mb; N50 = 11.3 Mb) with a low proportion of undetermined nucleotides (Ns = 0.34%) and with high gene completeness (BUSCO complete 96.9%). A genome-wide survey of resistance gene analogues identified 695 genes encoding nucleotide-binding site domains, receptor-like protein kinases, receptor-like proteins and transmembrane-coiled coil domains. Based on a high-density double-digest restriction-site associated DNA sequencing-based linkage map, the genome assembly was anchored to six pseudomolecules representing ∼73.5% of the whole genome assembly. In addition, we used SOL_r1.1 to identify quantitative trait loci for bolting timing and fruit/seed shape, which harbour biologically plausible candidate genes, such as homologues of the FLOWERING LOCUS T and EPIDERMAL PATTERNING FACTOR-LIKE genes. The new genome assembly, SOL_r1.1, will serve as a useful resource for identifying loci associated with important agronomic traits and for developing molecular markers for spinach breeding/selection programs.

Targeted base editing in the plastid genome of Arabidopsis thaliana

Bacterial cytidine deaminase fused to the DNA binding domains of transcription activator-like effector nucleases was recently reported to transiently substitute a targeted C to a T in mitochondrial DNA of mammalian cultured cells¹. We applied this system to targeted base editing in the Arabidopsis thaliana plastid genome. The targeted Cs were homoplasmically substituted to Ts in some plantlets of the T₁ generation and the mutations were inherited by their offspring independently of their nuclear-introduced vectors.

Platanus_B: an accurate de novo assembler for bacterial genomes using an iterative error-removal process

De novo assembly of short DNA reads remains an essential technology, especially for large-scale projects and high-resolution variant analyses in epidemiology. However, the existing tools often lack sufficient accuracy required to compare closely related strains. To facilitate such studies on bacterial genomes, we developed Platanus_B, a de novo assembler that employs iterations of multiple error-removal algorithms. The benchmarks demonstrated the superior accuracy and high contiguity of Platanus_B, in addition to its ability to enhance the hybrid assembly of both short and nanopore long reads. Although the hybrid strategies for short and long reads were effective in achieving near full-length genomes, we found that short-read-only assemblies generated with Platanus_B were sufficient to obtain ≥90% of exact coding sequences in most cases. In addition, while nanopore long-read-only assemblies lacked fine-scale accuracies, inclusion of short reads was effective in improving the accuracies. Platanus_B can, therefore, be used for comprehensive genomic surveillances of bacterial pathogens and high-resolution phylogenomic analyses of a wide range of bacteria.

The full mitochondrial genome sequence of the greater argonaut Argonauta argo (Cephalopoda, Argonautoidea) and its phylogenetic position in Octopodiformes

The greater argonaut Argonauta argo is a species of the paper nautilus (Argonautidae), which is a family in Octopoda. In this paper, we report its full mitogenome sequence, which was obtained from a specimen collected in the Japan Seas near Oki Island, Shimane Prefecture, in Japan. The sequence was determined using the NGS Illumina HiSeq platform. With its 37 genes, the mitogenome shows a typical metazoan and Octopoda genomic structure, and similar to the mitogenome of the previously reported congener, A. hians. To confirm A. argo phylogenetic position in Octopoda, we conducted maximum likelihood phylogenetic analysis, using a data set including publicly available 17 Octopodiformes, five Decapodiformes, three Nautiloids and two outgroup Conchiferans. The result confirmed the affinity of Argonautidae to Tremoctopus, and the sister group position of this clade against the rest of incirrate Octopods. The mitogenome and phylogeny of A. argo reported here will be useful for future studies involving this enigmatic species, including on the reacquisition of external calcified shell structures in mollusks.

Hamster PIWI proteins bind to piRNAs with stage-specific size variations during oocyte maturation

In animal gonads, transposable elements are actively repressed to preserve genome integrity through the PIWI-interacting RNA (piRNA) pathway. In mice, piRNAs are abundantly expressed in male germ cells, and form effector complexes with three distinct PIWIs. The depletion of individual Piwi genes causes male-specific sterility with no discernible phenotype in female mice. Unlike mice, most other mammals have four PIWI genes, some of which are expressed in the ovary. Here, purification of PIWI complexes from oocytes of the golden hamster revealed that the size of the PIWIL1-associated piRNAs changed during oocyte maturation. In contrast, PIWIL3, an ovary-specific PIWI in most mammals, associates with short piRNAs only in metaphase II oocytes, which coincides with intense phosphorylation of the protein. An improved high-quality genome assembly and annotation revealed that PIWIL1- and PIWIL3-associated piRNAs appear to share the 5′-ends of common piRNA precursors and are mostly derived from unannotated sequences with a diminished contribution from TE-derived sequences, most of which correspond to endogenous retroviruses. Our findings show the complex and dynamic nature of biogenesis of piRNAs in hamster oocytes, and together with the new genome sequence generated, serve as the foundation for developing useful models to study the piRNA pathway in mammalian oocytes.

Genomic Signatures for Species-Specific Adaptation in Lake Victoria Cichlids Derived from Large-Scale Standing Genetic Variation

The cichlids of Lake Victoria are a textbook example of adaptive radiation, as >500 endemic species arose in just 14,600 years. The degree of genetic differentiation among species is very low due to the short period of time after the radiation, which allows us to ascertain highly differentiated genes that are strong candidates for driving speciation and adaptation. Previous studies have revealed the critical contribution of vision to speciation by showing the existence of highly differentiated alleles in the visual opsin gene among species with different habitat depths. In contrast, the processes of species-specific adaptation to different ecological backgrounds remain to be investigated. Here, we used genome-wide comparative analyses of three species of Lake Victoria cichlids that inhabit different environments-Haplochromis chilotes, H. sauvagei, and Lithochromis rufus-to elucidate the processes of adaptation by estimating population history and by searching for candidate genes that contribute to adaptation. The patterns of changes in population size were quite distinct among the species according to their habitats. We identified many novel adaptive candidate genes, some of which had surprisingly long divergent haplotypes between species, thus showing the footprint of selective sweep events. Molecular phylogenetic analyses revealed that a large fraction of the allelic diversity among Lake Victoria cichlids was derived from standing genetic variation that originated before the adaptive radiation. Our analyses uncovered the processes of species-specific adaptation of Lake Victoria cichlids and the complexity of the genomic substrate that facilitated this adaptation.

Genomics and Transcriptomics of the green mussel explain the durability of its byssus

Mussels, which occupy important positions in marine ecosystems, attach tightly to underwater substrates using a proteinaceous holdfast known as the byssus, which is tough, durable, and resistant to enzymatic degradation. Although various byssal proteins have been identified, the mechanisms by which it achieves such durability are unknown. Here we report comprehensive identification of genes involved in byssus formation through whole-genome and foot-specific transcriptomic analyses of the green mussel, Perna viridis. Interestingly, proteins encoded by highly expressed genes include proteinase inhibitors and defense proteins, including lysozyme and lectins, in addition to structural proteins and protein modification enzymes that probably catalyze polymerization and insolubilization. This assemblage of structural and protective molecules constitutes a multi-pronged strategy to render the byssus highly resistant to environmental insults.

Comparative genomics of Glandirana rugosa using unsupervised AI reveals a high CG frequency

Genome sequencing of a unique frog (Glandirana rugosa) having XY/ZW systems within the species and comparative genomics with other six frogs were performed using a batch-learning self-organizing map, which is unsupervised AI for oligonucleotide compositions, to clarify its genome characteristics.

The Japanese wrinkled frog (Glandirana rugosa) is unique in having both XX-XY and ZZ-ZW types of sex chromosomes within the species. The genome sequencing and comparative genomics with other frogs should be important to understand mechanisms of turnover of sex chromosomes within one species or during a short period. In this study, we analyzed the newly sequenced genome of G. rugosa using a batch-learning self-organizing map which is unsupervised artificial intelligence for oligonucleotide compositions. To clarify genome characteristics of G. rugosa, we compared its short oligonucleotide compositions in all 1-Mb genomic fragments with those of other six frog species (Pyxicephalus adspersus, Rhinella marina, Spea multiplicata, Leptobrachium leishanense, Xenopus laevis, and Xenopus tropicalis). In G. rugosa, we found an Mb-level large size of repeat sequences having a high identity with the W chromosome of the African bullfrog (P. adspersus). Our study concluded that G. rugosa has unique genome characteristics with a high CG frequency, and its genome is assumed to heterochromatinize a large size of genome via methylataion of CG.

Targeted gene disruption of ATP synthases 6-1 and 6-2 in the mitochondrial genome of Arabidopsis thaliana by mitoTALENs

We recently achieved targeted disruptions of cytoplasmic male sterility (CMS)-associated genes in the mitochondrial genomes of rice and rapeseed by using mitochondria-targeted transcription activator-like effector nucleases (mitoTALENs). It was the first report of stable and heritable targeted gene modification of plant mitochondrial genomes. Here, we attempted to use mitoTALENs to disrupt two mitochondrial genes in the model plant Arabidopsis thaliana(Arabidopsis) using three different promoters and two types of TALENs. The targets were the two isoforms of the ATP synthase subunit 6 gene, atp6-1 and atp6-2. Each of these genes was successfully deleted and the mitochondrial genomes were recovered in a homoplasmic state. The nuclear genome also has a copy of atp6-1, and we were able to confirm that it was the mitochondrial gene and not the nuclear pseudogene that was knocked out. Among the three mitoTALEN promoters tried, the RPS5A promoter was the most effective. Conventional mitoTALENs were more effective than single-molecule mito-compactTALENs. Targeted mitochondrial gene deletion was achieved by crossing as well as by floral-dip transformation to introduce the mitoTALEN constructs into the nucleus. The gene disruptions were caused by large (kb-size) deletions. The ends of the remaining sequences were connected to distant loci, mostly by illegitimate homologous recombinations between repeats.

Impact of east Japan earthquake disaster with massive tsunami for prevalence of Takotsubo syndrome – a multicenter regional registry before and after east Japan earthquake disaster

Background

Previous studies demonstrated Takotsubo syndrome (TS) was easy provoked by earthquake disaster. However, a previous other regional report demonstrated TS was not increased after 2011 east Japan earthquake disaster. The purpose of this study was to clarify incidence of TS after the earthquake disaster in Iwate prefecture during long term period.

Method

Consecutive hospitalized TS patients were registered during 8 years between 2009 and 2016 in our medical university and five Iwate prefecture hospitals. Moreover, patients were divided into two groups, i.e., those with the inland and those with tsunami-stricken area groups. Prevalence of TS were calculated by standard incidence ratio (SIR) before and after the earthquake disaster. Moreover, long-term prognosis in the both groups was compared using Kaplan-Meier analysis.

Results

A total of 112 TS (male 25 and female 87) were registered from acute coronary syndrome registry in each hospital (n=4,163). Averaged age was 75.3 year-old. A total number of TS just after the two months of the earthquake (March and April 2011) was nine and significance monthly variation was observed comparing with the other months (p=0.029). SIR before and after the disaster is as following Figure. There were no significant differences for long-term prognosis between the two groups (p=0.20).

Conclusion

Incidence of TS was increased in acute phase after east Japan earthquake disaster. However, significance increases were maintained during long-term period, although number of TS was decreased after acute phase. TS is increased not only acute but also chronic phase after the serious earthquake disaster.

Standard incidence ratio

Funding Acknowledgement

Type of funding source: None

Parallel reductive genome evolution in Desulfovibrio ectosymbionts independently acquired by Trichonympha protists in the termite gut

Several Trichonympha protist species in the termite gut have independently acquired Desulfovibrio ectosymbionts in apparently different stages of symbiosis. Here, we obtained the near-complete genome sequence of Desulfovibrio phylotype ZnDsv-02, which attaches to the surface of Trichonympha collaris cells, and compared it with a previously obtained genome sequence of 'Candidatus Desulfovibrio trichonymphae' phylotype Rs-N31, which is almost completely embedded in the cytoplasm of Trichonympha agilis. Single-nucleotide polymorphism analysis indicated that although Rs-N31 is almost clonal, the ZnDsv-02 population on a single host cell is heterogeneous. Despite these differences, the genome of ZnDsv-02 has been reduced to 1.6 Mb, which is comparable to that of Rs-N31 (1.4 Mb), but unlike other known ectosymbionts of protists with a genome similar in size to their free-living relatives. Except for the presence of a lactate utilization pathway, cell-adhesion components and anti-phage defense systems in ZnDsv-02, the overall gene-loss pattern between the two genomes is very similar, including the loss of genes responsive to environmental changes. Our study suggests that genome reduction can occur in ectosymbionts, even when they can be transmitted horizontally and obtain genes via lateral transfer, and that the symbiont genome size depends heavily on their role in the symbiotic system.

A ubiquitous subcuticular bacterial symbiont of a coral predator, the crown-of-thorns starfish, in the Indo-Pacific

BACKGROUND

Population outbreaks of the crown-of-thorns starfish (Acanthaster planci sensu lato; COTS), a primary predator of reef-building corals in the Indo-Pacific Ocean, are a major threat to coral reefs. While biological and ecological knowledge of COTS has been accumulating since the 1960s, little is known about its associated bacteria. The aim of this study was to provide fundamental information on the dominant COTS-associated bacteria through a multifaceted molecular approach.

METHODS

A total of 205 COTS individuals from 17 locations throughout the Indo-Pacific Ocean were examined for the presence of COTS-associated bacteria. We conducted 16S rRNA metabarcoding of COTS to determine the bacterial profiles of different parts of the body and generated a full-length 16S rRNA gene sequence from a single dominant bacterium, which we designated COTS27. We performed phylogenetic analysis to determine the taxonomy, screening of COTS27 across the Indo-Pacific, FISH to visualize it within the COTS tissues, and reconstruction of the bacterial genome from the hologenome sequence data.

RESULTS

We discovered that a single bacterium exists at high densities in the subcuticular space in COTS forming a biofilm-like structure between the cuticle and the epidermis. COTS27 belongs to a clade that presumably represents a distinct order (so-called marine spirochetes) in the phylum Spirochaetes and is universally present in COTS throughout the Indo-Pacific Ocean. The reconstructed genome of COTS27 includes some genetic traits that are probably linked to adaptation to marine environments and evolution as an extracellular endosymbiont in subcuticular spaces.

CONCLUSIONS

COTS27 can be found in three allopatric COTS species, ranging from the northern Red Sea to the Pacific, implying that the symbiotic relationship arose before the speciation events (approximately 2 million years ago). The universal association of COTS27 with COTS and nearly mono-specific association at least with the Indo-Pacific COTS provides a useful model system for studying symbiont-host interactions in marine invertebrates and may have applications for coral reef conservation. Video Abstract.

A genome database for a Japanese population of the larvacean Oikopleura dioica

The larvacean Oikopleura dioica is a planktonic chordate and is a tunicate that belongs to the closest relatives to vertebrates. Its simple and transparent body, invariant embryonic cell lineages, and short life cycle of 5 days make it a promising model organism for the study of developmental biology. The genome browser OikoBase was established in 2013 using Norwegian O. dioica. However, genome information for other populations is not available, even though many researchers have studied local populations. In the present study, we sequenced using Illumina and PacBio RSII technologies the genome of O. dioica from a southwestern Japanese population that was cultured in our laboratory for 3 years. The genome of Japanese O. dioica was assembled into 576 scaffold sequences with a total length and N50 length of 56.6 and 1.5 Mb, respectively. A total of 18,743 gene models (transcript models) were predicted in the genome assembly, named OSKA2016. In addition, 19,277 non-redundant transcripts were assembled using RNA-seq data. The OSKA2016 has global sequence similarity of only 86.5% when compared with the OikoBase, highlighting the sequence difference between the two far distant O. dioica populations on the globe. The genome assembly, transcript assembly, and transcript models were incorporated into ANISEED (https://www.aniseed.cnrs.fr/) for genome browsing and BLAST searches. Mapping of reads obtained from male- or female-specific genome libraries yielded male-specific scaffolds in the OSKA2016 and revealed that over 2.6 Mb of sequence were included in the male-specific Y-region. The genome and transcriptome resources from two distinct populations will be useful datasets for developmental biology, evolutionary biology, and molecular ecology using this model organism.

0001 Activation of Nociceptin/Orphanin-FQ Peptide (NOP) Receptors Produces an Increase in Non-REM Sleep in Rats and Constitutes a Novel and Attractive Target for the Treatment of Insomnia

Introduction

Treatments for insomnia have targeted GABA, histamine, serotonin, melatonin and orexin receptors. The nociceptin/orphanin-FQ peptide (NOP) receptor is widely expressed in the nervous system. High doses of NOP agonists administered systemically or locally into the CNS can result in sedation, however, the utility of targeting this receptor to treat insomnia has not been fully described.

Methods

V117957 is a recently described investigational oral, potent and selective NOP receptor partial agonist. We determined the brain Kp in whole brain and multiple sub-regions (50mg/kg) and receptor occupancy in the hypothalamus (30, 300mg/kg) via in vivo displacement using [3H]-NOP-1A. EEG/EMG were determined in rats chronically implanted with electrodes (cortex and dorsal neck muscle) and recorded via telemetry following dosing (3, 30, 300mg/kg); sleep stage was determined from visual analysis of EEG level. Sleep parameters were also assessed in NOP receptor knock-out rats (300mg/kg). The side-effect profile for V117957 was determined by functional observation battery, whole-body plethysmography, Morris water maze (MWM) (up to 600mg/kg) and conditioned place preference (CPP) assay (up to 300mg/kg).

Results

V117957 displayed limited distribution into the CNS but achieved a high level of receptor occupancy (75% at 30mg/kg). Administration of V117957 produced dose-dependent and statistically significant increases in non-REM sleep with a minimally efficacious dose of 30mg/kg; a coincident dose-dependent and statistically significant decrease in wakefulness and a non-dose-dependent effect on REM sleep occurred. These changes were not seen in knock-out animals demonstrating effects are via NOP receptors. At doses higher than those that increased non-REM sleep, V117957 had no effects in a functional observational battery, did not affect escape latency in MWM or produce CPP; additionally, V117957 did not affect respiratory parameters.

Conclusion

We conclude that activation of NOP receptors decreases wakefulness and increases non-REM sleep in rats with an improved preclinical profile compared to historical profiles of current treatments and, therefore, may represent a novel and attractive target for the treatment of insomnia.

Support

Funded by Shionogi and Imbrium Therapeutics, a subsidiary of Purdue Pharma L.P.

Differential dynamics and impacts of prophages and plasmids on the pangenome and virulence factor repertoires of Shiga toxin-producing Escherichia coli O145:H28

Phages and plasmids play important roles in bacterial evolution and diversification. Although many draft genomes have been generated, phage and plasmid genomes are usually fragmented, limiting our understanding of their dynamics. Here, we performed a systematic analysis of 239 draft genomes and 7 complete genomes of Shiga toxin (Stx)-producing Escherichia coli O145:H28, the major virulence factors of which are encoded by prophages (PPs) or plasmids. The results indicated that PPs are more stably maintained than plasmids. A set of ancestrally acquired PPs was well conserved, while various PPs, including Stx phages, were acquired by multiple sublineages. In contrast, gains and losses of a wide range of plasmids have frequently occurred across the O145:H28 lineage, and only the virulence plasmid was well conserved. The different dynamics of PPs and plasmids have differentially impacted the pangenome of O145:H28, with high proportions of PP- and plasmid-associated genes in the variably present and rare gene fractions, respectively. The dynamics of PPs and plasmids have also strongly impacted virulence gene repertoires, such as the highly variable distribution of stx genes and the high conservation of a set of type III secretion effectors, which probably represents the core effectors of O145:H28 and the genes on the virulence plasmid in the entire O145:H28 population. These results provide detailed insights into the dynamics of PPs and plasmids, and show the application of genomic analyses using a large set of draft genomes and appropriately selected complete genomes.

P3615The number of coronary risk factors and mortality in patients with acute myocardial infarction from Japanese nation-wide real-world database

Background

Hypertension, diabetes, dyslipidemia and smoking are so-called coronary risk factors for coronary heart disease, which were established by extensive epidemiological research. However, in Japanese patients with acute myocardial infarction (AMI), the impact of number of coronary risk factors on in-hospital morality has not been elucidated.

Methods

The Japan Acute Myocardial Infarction Registry (JAMIR) is a nationwide real-world database integrated form 10 regional registries. We examined the association between number of coronary risk factors and in-hospital mortality from this JAMIR registry.

Results

The data were obtained from total of 20462 AMI patients (mean age, 68.8±13.3 years old; 15281 men, 5181 women). Figure 1 shows the prevalence of each coronary risk factors stratified by sex and decade. The prevalence of hypertension became higher with the advanced age while the prevalence of smoking became lower with the advanced age. Prevalence of diabetes and dyslipidemia were highest in middle age. Majority (76.9%) of the patients with AMI had at least 1 of these coronary risk factors and, 23.1% had none of them. Overall, except women under 50, number of coronary risk factor was relatively less in older age (Figure 2). In-hospital mortality by sex and decades was shown in figure 3. In-hospital mortality rates were 10.7%, 10.5%, 7.2%, 5.0% and 4.5% with 0, 1, 2, 3 and 4 risk factors, respectively (Figure 4A). After adjusting age and sex, there was an inverse association between the number of coronary risk factors and in-hospital mortality (adjusted odds ratio [1.68; 95% CI, 1.20–2.35] among individuals with 0 vs. 4 risk factors, Figure 4B).

Conclusion

In the present study of Japanese patients with AMI, who received modern medical treatment, in-hospital mortality was inversely related to the number of coronary risk factors.

Acknowledgement/Funding

Grant-in-Aid for Scientific Research

Multi-step genomic dissection of a suspected intra-hospital Helicobacter cinaedi outbreak

Helicobacter cinaedi is an emerging pathogen causing bacteraemia and cellulitis. Nosocomial transmission of this microbe has been described, but detailed molecular-epidemiological analyses have not been performed. Here, we describe the results of a multi-step genome-wide phylogenetic analysis of a suspected intra-hospital outbreak of H. cinaedi that occurred in a hospital in Japan. The outbreak was recognized by the infectious control team (ICT) of the hospital as a sudden increase in H. cinaedi bacteraemia. ICT defined this outbreak case based on 16S rRNA sequence data and epidemiological information, but were unable to determine the source and route of the infections. We therefore re-investigated this case using whole-genome sequencing (WGS). We first performed a species-wide analysis using publicly available genome sequences to understand the level of genomic diversity of this under-studied species. The clusters identified were then separately analysed using the genome sequence of a representative strain in each cluster as a reference. These analyses provided a high-level phylogenetic resolution of each cluster, identified a confident set of outbreak isolates, and discriminated them from other closely related but distinct clones, which were locally circulating and invaded the hospital during the same period. By considering the epidemiological data, possible strain transmission chains were inferred, which highlighted the role of asymptomatic carriers or environmental contamination. The emergence of a subclone with increased resistance to fluoroquinolones in the outbreak was also recognized. Our results demonstrate the impact of the use of a closely related genome as a reference to maximize the power of WGS.

Bcl11b controls odorant receptor class choice in mice

Each olfactory sensory neuron (OSN) expresses a single odorant receptor (OR) gene from the class I or class II repertoire in mice. The mechanisms that regulate OR class choice in OSNs remain unknown. Here, we show that the transcription factor Bcl11b determines the OR class to be expressed in OSNs. Both loss- and gain-of-function analyses demonstrate that class I is a default fate of OSNs and that Bcl11b dictates a class II OR choice by suppressing the effect of the J-element, a class I-OR enhancer. We further demonstrate that OSN-specific genetic manipulations of Bcl11b bias the OR class choice, generating mice with "class I-dominant" and "class II-dominant" noses, which display contrasting innate olfactory behaviors to two distinct aversive odorants. Overall, these findings reveal a unique transcriptional mechanism mediating a binary switch for OR class choice that is crucial to both the anatomical and functional organization of the olfactory system.

Curing cytoplasmic male sterility via TALEN-mediated mitochondrial genome editing

Sequence-specific nucleases are commonly used to modify the nuclear genome of plants. However, targeted modification of the mitochondrial genome of land plants has not yet been achieved. In plants, a type of male sterility called cytoplasmic male sterility (CMS) has been attributed to certain mitochondrial genes, but none of these genes has been validated by direct mitochondrial gene-targeted modification. Here, we knocked out CMS-associated genes (orf79 and orf125) of CMS varieties of rice and rapeseed, respectively, using transcription activator-like effector nucleases (TALENs) with mitochondria localization signals (mitoTALENs). We demonstrate that knocking out these genes cures male sterility, strongly suggesting that these genes are causes of CMS. Sequencing revealed that double-strand breaks induced by mitoTALENs were repaired by homologous recombination, and that during this process, the target genes and surrounding sequences were deleted. Our results show that mitoTALENs can be used to stably and heritably modify the mitochondrial genome in plants.

Genomic Characterization of β-Glucuronidase-Positive Escherichia coli O157:H7 Producing Stx2a

Among Shiga toxin (Stx)–producing Escherichia coli (STEC) O157:H7 strains, those producing Stx2a cause more severe diseases. Atypical STEC O157:H7 strains showing a β-glucuronidase–positive phenotype (GP STEC O157:H7) have rarely been isolated from humans, mostly from persons with asymptomatic or mild infections; Stx2a-producing strains have not been reported. We isolated, from a patient with bloody diarrhea, a GP STEC O157:H7 strain (PV15-279) that produces Stx2a in addition to Stx1a and Stx2c. Genomic comparison with other STEC O157 strains revealed that PV15-279 recently emerged from the stx1a/stx2c-positive GP STEC O157:H7 clone circulating in Japan. Major virulence genes are shared between typical (β-glucuronidase–negative) and GP STEC O157:H7 strains, and the Stx2-producing ability of PV15-279 is comparable to that of typical STEC O157:H7 strains; therefore, PV15-279 presents a virulence potential similar to that of typical STEC O157:H7. This study reveals the importance of GP O157:H7 as a source of highly pathogenic STEC clones.

Evaluation of SNP calling methods for closely related bacterial isolates and a novel high-accuracy pipeline: BactSNP

Bacteria are highly diverse, even within a species; thus, there have been many studies which classify a single species into multiple types and analyze the genetic differences between them. Recently, the use of whole-genome sequencing (WGS) has been popular for these analyses, and the identification of single-nucleotide polymorphisms (SNPs) between isolates is the most basic analysis performed following WGS. The performance of SNP-calling methods therefore has a significant effect on the accuracy of downstream analyses, such as phylogenetic tree inference. In particular, when closely related isolates are analyzed, e.g. in outbreak investigations, some SNP callers tend to detect a high number of false-positive SNPs compared with the limited number of true SNPs among isolates. However, the performances of various SNP callers in such a situation have not been validated sufficiently. Here, we show the results of realistic benchmarks of commonly used SNP callers, revealing that some of them exhibit markedly low accuracy when target isolates are closely related. As an alternative, we developed a novel pipeline BactSNP, which utilizes both assembly and mapping information and is capable of highly accurate and sensitive SNP calling in a single step. BactSNP is also able to call SNPs among isolates when the reference genome is a draft one or even when the user does not input the reference genome. BactSNP is available at https://github.com/IEkAdN/BactSNP.