Characterization of the pufferfish (Fugu) genome as a compact model vertebrate genome

Transcriptome Analysis Reveals the Immunosuppression in Tiger Pufferfish (Takifugu rubripes) under Cryptocaryon irritans Infection

The tiger pufferfish (Takifugu rubripes), also known as fugu, has recently suffered from severe C. irritans infections under aquaculture environment, yet the underlying immune mechanisms against the parasite remain poorly understood. In this study, we conducted a comprehensive transcriptome analysis of the gill tissue from infected and uninfected fish using PacBio long-read (one pooled sample each for seriously infected and healthy individuals, respectively) and Illumina short-read (three pools for mildly infected, seriously infected, and healthy individuals, respectively) RNA sequencing technologies. After aligning sequence data to fugu's reference genome, 47,307 and 34,413 known full-length transcripts were identified and profiled in healthy and infected fish, respectively. Similarly, we identified and profiled 1126 and 803 novel genes that were obtained from healthy and infected fish, respectively. Interestingly, we found a decrease in the number of alternative splicing (AS) events and long non-coding RNAs (lncRNAs) after infection with C. irritans, suggesting that they may be involved in the regulation of the immune response in fugu. There were 687 and 1535 differentially expressed genes (DEGs) in moderately and heavily infected fish, respectively, compared to uninfected fish. Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that immune-related DEGs in the two comparison groups were mainly enriched in cytokine-cytokine receptor interactions, ECM-receptor interactions, T-cell receptor signaling pathways, Th1 and Th2 cell differentiation, and Th17 cell differentiation pathways. Further analysis revealed that a large number of immune-related genes were downregulated in infected fish relative to uninfected ones, such as CCR7, IL7R, TNFRSF21, CD4, COL2A1, FOXP3B, and ITGA8. Our study suggests that C. irritans is potentially a highly efficient parasite that may disrupt the defense mechanisms of fugu against it. In addition, in combination of short-read RNA sequencing and previous genome-wide association analyses, we identified five key genes (NDUFB6, PRELID1, SMOX, SLC25A4, and DENND1B) that might be closely associated with C. irritans resistance. This study not only provides valuable resources of novel genic transcripts for further research, but also provides new insights into the immune mechanisms underlying C. irritans infection response in farmed fugu.

Database construction and comparative genomics analysis of genes involved in nutritional metabolic diseases in fish

Nutritional metabolic diseases in fish frequently arise in the setting of intensive aquaculture. The etiology and pathogenesis of these conditions involve energy metabolic disorders influenced by both internal genetic factors and external environmental conditions. The exploration of genes associated with nutritional and metabolic disorder has sparked considerable interest within both the aquaculture scientific community and the industry. High-throughput sequencing technology offers researchers extensive genetic information. Effectively mining, analyzing, and securely storing this data is crucial, especially for advancing disease prevention and treatment strategies. Presently, the exploration and application of gene databases concerning nutritional and metabolic disorders in fish are at a nascent stag. Therefore, this study focused on the model organism zebrafish and five primary economic fish species as the subjects of investigation. Using information from KEGG, OMIM, and existing literature, a novel gene database associated with nutritional metabolic diseases in fish was meticulously constructed. This database encompassed 4583 genes for Danio rerio, 6287 for Cyprinus carpio, 3289 for Takifugu rubripes, 3548 for Larimichthys crocea, 3816 for Oreochromis niloticus, and 5708 for Oncorhynchus mykiss. Through a comparative systems biology approach, we discerned a relatively high conservation of genes linked to nutritional metabolic diseases across these fish species, with over 54.9 % of genes being conserved throughout all six species. Additionally, the analysis pinpointed the existence of 13 species-specific genes within the genomes of large yellow croaker, tilapia, and rainbow trout. These genes exhibit the potential to serve as novel candidate targets for addressing nutritional metabolic diseases.

High Diversity of Long Terminal Repeat Retrotransposons in Compact Vertebrate Genomes: Insights from Genomes of Tetraodontiformes

This study aimed to investigate the evolutionary profile (including diversity, activity, and abundance) of retrotransposons (RTNs) with long terminal repeats (LTRs) in ten species of Tetraodontiformes. These species, Arothron firmamentum, Lagocephalus sceleratus, Pao palembangensis, Takifugu bimaculatus, Takifugu flavidus, Takifugu ocellatus, Takifugu rubripes, Tetraodon nigroviridis, Mola mola, and Thamnaconus septentrionalis, are known for having the smallest genomes among vertebrates. Data mining revealed a high diversity and wide distribution of LTR retrotransposons (LTR-RTNs) in these compact vertebrate genomes, with varying abundances among species. A total of 819 full-length LTR-RTN sequences were identified across these genomes, categorized into nine families belonging to four different superfamilies: ERV (Orthoretrovirinae and Epsilon retrovirus), Copia, BEL-PAO, and Gypsy (Gmr, Mag, V-clade, CsRN1, and Barthez). The Gypsy superfamily exhibited the highest diversity. LTR family distribution varied among species, with Takifugu bimaculatus, Takifugu flavidus, Takifugu ocellatus, and Takifugu rubripes having the highest richness of LTR families and sequences. Additionally, evidence of recent invasions was observed in specific tetraodontiform genomes, suggesting potential transposition activity. This study provides insights into the evolution of LTR retrotransposons in Tetraodontiformes, enhancing our understanding of their impact on the structure and evolution of host genomes.

Parallel gene size and isoform expansion of ancient neuronal genes

How nervous systems evolved is a central question in biology. A diversity of synaptic proteins is thought to play a central role in the formation of specific synapses leading to nervous system complexity. The largest animal genes, often spanning hundreds of thousands of base pairs, are known to be enriched for expression in neurons at synapses and are frequently mutated or misregulated in neurological disorders and diseases. Although many of these genes have been studied independently in the context of nervous system evolution and disease, general principles underlying their parallel evolution remain unknown. To investigate this, we directly compared orthologous gene sizes across eukaryotes. By comparing relative gene sizes within organisms, we identified a distinct class of large genes with origins predating the diversification of animals and, in many cases, the emergence of neurons as dedicated cell types. We traced this class of ancient large genes through evolution and found orthologs of the large synaptic genes potentially driving the immense complexity of metazoan nervous systems, including in humans and cephalopods. Moreover, we found that while these genes are evolving under strong purifying selection, as demonstrated by low dN/dS ratios, they have simultaneously grown larger and gained the most isoforms in animals. This work provides a new lens through which to view this distinctive class of large and multi-isoform genes and demonstrates how intrinsic genomic properties, such as gene length, can provide flexibility in molecular evolution and allow groups of genes and their host organisms to evolve toward complexity.

PiggyBac Transposon Mining in the Small Genomes of Animals

TEs, including DNA transposons, are major contributors of genome expansions, and have played a very significant role in shaping the evolution of animal genomes, due to their capacity to jump from one genomic position to the other. In this study, we investigated the evolution landscapes of PB transposons, including their distribution, diversity, activity and structure organization in 79 species of small (compact) genomes of animals comprising both vertebrate and invertebrates. Overall, 212 PB transposon types were detected from almost half (37) of the total number of the small genome species (79) investigated. The detected PB transposon types, which were unevenly distributed in various genera and phyla, have been classified into seven distinct clades or families with good bootstrap support (>80%). The PB transposon types that were identified have a length ranging from 1.23 kb to 9.51 kb. They encode transposases of approximately ≥500 amino acids in length, and possess terminal inverted repeats (TIRs) ranging from 4 bp to 24 bp. Though some of the transposon types have long TIRs (528 bp), they still maintain the consistent and reliable 4 bp target site duplication (TSD) of TTAA. However, PiggyBac-2_Cvir transposon originating from the Crassostrea virginica species exhibits a unique TSD of TATG. The TIRs of the transposons in all the seven families display high divergence, with a highly conserved 5' end motif. The core transposase domains (DDD) were better conserved among the seven different families compared to the other protein domains, which were less prevalent in the vertebrate genome. The divergent evolution dynamics analysis also indicated that the majority of the PB transposon types identified in this study are either relatively young or old, with some being active. Additionally, numerous invasions of PB transposons were found in the genomes of both vertebrate and invertebrate animals. The data reveals that the PB superfamily is widely distributed in these species. PB transposons exhibit high diversity and activity in the small genomes of animals, and might play a crucial role in shaping the evolution of these small genomes of animals.

Tandem mass tag (TMT) quantitative proteomics and phosphoproteomic of Takifugu rubripes infected with Cryptocaryon irritans

In this study, we identified the differentially expressed proteins in gills stimulated by infected ciliates and analyzed the immune mechanisms of T. rubripes infected with the ciliate Cryptocaryon irritans. Through liquid chromatography analysis, a total of 144 proteins were identified with significant differences, of which 58 were upregulated and 86 were downregulated. Among phosphorylated proteins, we identified a total of 167 significantly different phosphorylated proteins, of which 44 were upregulated, 123 were downregulated, 60 were upregulated, and 208 were downregulated. We analyzed the data of proteomics and Phosphorylated proteome quantification protein omics to finally identify three phosphorylated proteins (RPS27, eNOS and CaM) and two phosphorylated protein kinases(CaMKII and MAPK1) as potential biomarkers for T. rubripes immune responses. We finally identified three phosphorylated proteins (RPS27, eNOS and CaM) and two phosphorylated protein kinases (CaMKII and MAPK1) as potential biomarkers of immune response of T. rubripes. Our research findings provide new insights into the immune mechanism of T. rubripes, which may serve as an effective indicator of C. irritans infection in T. rubripes.

Ancient origins of complex neuronal genes

How nervous systems evolved is a central question in biology. An increasing diversity of synaptic proteins is thought to play a central role in the formation of specific synapses leading to nervous system complexity. The largest animal genes, often spanning millions of base pairs, are known to be enriched for expression in neurons at synapses and are frequently mutated or misregulated in neurological disorders and diseases. While many of these genes have been studied independently in the context of nervous system evolution and disease, general principles underlying their parallel evolution remain unknown. To investigate this, we directly compared orthologous gene sizes across eukaryotes. By comparing relative gene sizes within organisms, we identified a distinct class of large genes with origins predating the diversification of animals and in many cases the emergence of dedicated neuronal cell types. We traced this class of ancient large genes through evolution and found orthologs of the large synaptic genes driving the immense complexity of metazoan nervous systems, including in humans and cephalopods. Moreover, we found that while these genes are evolving under strong purifying selection as demonstrated by low dN/dS scores, they have simultaneously grown larger and gained the most isoforms in animals. This work provides a new lens through which to view this distinctive class of large and multi-isoform genes and demonstrates how intrinsic genomic properties, such as gene length, can provide flexibility in molecular evolution and allow groups of genes and their host organisms to evolve toward complexity.

RNA out of the mist

RNA has long been regarded primarily as the intermediate between genes and proteins. It was a surprise then to discover that eukaryotic genes are mosaics of mRNA sequences interrupted by large tracts of transcribed but untranslated sequences, and that multicellular organisms also express many long 'intergenic' and antisense noncoding RNAs (lncRNAs). The identification of small RNAs that regulate mRNA translation and half-life did not disturb the prevailing view that animals and plant genomes are full of evolutionary debris and that their development is mainly supervised by transcription factors. Gathering evidence to the contrary involved addressing the low conservation, expression, and genetic visibility of lncRNAs, demonstrating their cell-specific roles in cell and developmental biology, and their association with chromatin-modifying complexes and phase-separated domains. The emerging picture is that most lncRNAs are the products of genetic loci termed 'enhancers', which marshal generic effector proteins to their sites of action to control cell fate decisions during development.

Host Hybridization Dominates over Cohabitation in Affecting Gut Microbiota of Intrageneric Hybrid Takifugu Pufferfish

Microbial symbionts are of great importance for macroscopic life, including fish, and both collectively comprise an integrated biological entity known as the holobiont. Yet little is known as to how the normal balance within the fish holobiont is maintained and how it responds to biotic and/or abiotic influences. Here, through amplicon profiling, the genealogical relationship between artificial F1 hybrid pufferfish with growth heterosis, produced from crossing female Takifugu obscurus with male Takifugu rubripes and its maternal halfsibling purebred, was well recapitulated by their gut microbial community similarities, indicating an evident parallelism between host phylogeny (hybridity) and microbiota relationships therein. Interestingly, modest yet significant fish growth promotion and gut microbiota alteration mediated by hybrid-purebred cohabitation were observed, in comparison with their respective monoculture cohorts that share common genetic makeups, implying a certain degree of environmental influences. Moreover, the underlying assemblage patterns of gut microbial communities were found associated with a trade-off between variable selection and dispersal limitation, which are plausibly driven by the augmented social interactions between hybrid and purebred cohabitants differing in behaviors. Results from this study not only can enrich, from a microbial perspective, the sophisticated understanding of complex and dynamic assemblage of the fish holobiont, but will also provide deeper insights into the ecophysiological factors imposed on the diversity-function relationships thereof. Our findings emphasize the intimate associations of gut microbiota in host genetics-environmental interactions and would have deeper practical implications for microbial contributions to optimize performance prediction and to improve the production of farmed fishes. IMPORTANCE Microbial symbionts are of great importance for macroscopic life, including fish, and yet little is known as to how the normal balance within the fish holobiont is maintained and how it responds to the biotic and/or abiotic influences. Through gut microbiota profiling, we show that host intrageneric hybridization and cohabitation can impose a strong disturbance upon pufferfish gut microbiota. Moreover, marked alterations in the composition and function of gut microbiota in both hybrid and purebred pufferfish cohabitants were observed, which are potentially correlated with different metabolic priorities and behaviors between host genealogy. These results can enrich, from a microbial perspective, the sophisticated understanding of the complex and dynamic assemblage of the fish holobiont and would have deeper practical implications for microbial contributions to optimize performance prediction and to improve farmed fish production.

The functional verification and analysis of Fugu promoter of cardiac gene tnni1a in zebrafish

Troponin I type 1b (Tnni1b) is thought to be a novel isoform that is expressed only in the zebrafish heart. Knocking down of tnni1b can lead to cardiac defects in zebrafish. Although both the zebrafish tnni1b and human troponin I1 (TNNI1) genes are thought to be closely associated with fatal cardiac development, the regulatory molecular mechanisms of these genes are poorly understood. Analyzing the functionally conserved sequence, especially in the noncoding regulatory region involved in gene expression, clarified these mechanisms. In this study, we isolated a 3 kb fragment upstream of Fugu tnni1a that can regulate green fluorescence protein (GFP) expression in a heart-specific manner, similar to the pattern of zebrafish homologue expression. Three evolutionarily conserved regions (ECRs) in the 5'-flanking sequence of Fugu tnni1a were identified by sequence alignment. Deletion analysis led to the identification of ECR2 as a core sequence that affects the heart-specific expression function of the Fugu tnni1a promoter. Interestingly, both the Fugu tnni1a promoter and ECR2 sequence were functionally conserved in zebrafish, although they shared no sequence similarity. Together, the findings of our study provided further evidence for the important role of tnni1a homologous in cardiac development and demonstrated that two functionally conserved sequences in the zebrafish and Fugu genomes may be ECRs, despite their lack of similarity.

Transcriptome Analysis Identifies Key Metabolic Changes in the Brain of Takifugu rubripes in Response to Chronic Hypoxia

The brain is considered to be an extremely sensitive tissue to hypoxia, and the brain of fish plays an important role in regulating growth and adapting to environmental changes. As an important aquatic organism in northern China, the economic yield of Takifugu rubripes is deeply influenced by the oxygen content of seawater. In this regard, we performed RNA-seq analysis of T. rubripes brains under hypoxia and normoxia to reveal the expression patterns of genes involved in the hypoxic response and their enrichment of metabolic pathways. Studies have shown that carbohydrate, lipid and amino acid metabolism are significant pathways for the enrichment of differentially expressed genes (DEGs) and that DEGs are significantly upregulated in those pathways. In addition, some biological processes such as the immune system and signal transduction, where enrichment is not significant but important, are also discussed. Interestingly, the DEGs associated with those pathways were significantly downregulated or inhibited. The present study reveals the mechanism of hypoxia tolerance in T. rubripes at the transcriptional level and provides a useful resource for studying the energy metabolism mechanism of hypoxia response in this species.

Life cycle assessment of tiger puffer (Takifugu rubripes) farming: A case study in Dalian, China

In China, energy consumption and carbon emission by the aquaculture industry have become major problems. The tiger puffer (Takifugu rubripes) is an emerging aquaculture species in China, but its environmental impact during the farming process has not yet been evaluated systematically. To the best of our knowledge, this is the first life cycle assessment (LCA) of tiger puffer land-sea relay strategy in Dalian, China. To analyze the environmental impact of the tiger puffer farming process, the following four stages were considered: seed rearing, deep-sea cage farming-1, industrial recirculating aquaculture, and deep-sea cage farming-2. The LCA software GaBi 10.5 academy version and CML-IA-Jan. 2016-world method were used to calculate the environmental impacts. According to the LCA results, marine aquatic ecotoxicity potential was the largest contributor to the environmental impact, and industrial recirculating aquaculture was the largest farming stage in the whole tiger puffer farming process. Energy in the form of electricity, coal, and gasoline was consumed to maintain the power supply in the tiger puffer farming process, and it was a key factor that influenced the environmental performance. Based on the sensitivity and energy analyses, energy consumption for equipment operation at the industrial recirculating aquaculture stage, feed consumption, and gasoline consumption for transportation at the deep-sea cage farming-2 stage need to be carefully considered. The following improvement measures were suggested to improve the environmental performance of tiger puffer farming and the aquaculture industry: establish electricity, wind power, and solar energy integrated management systems; ex-ante LCA for parameter optimization in future technology research and development; and new production strategies such as aquaponics and integrated multi-trophic aquaculture. Moreover, life cycle inventory (LCI) of tiger puffer land-sea relay farming was established to obtain essential information, enrich aquaculture LCI databases, and support aquaculture LCA research.

Comparative Metabolomics and Proteomics Reveal Vibrio parahaemolyticus Targets Hypoxia-Related Signaling Pathways of Takifugu obscurus

Coronavirus disease 2019 (COVID-19) raises the issue of how hypoxia destroys normal physiological function and host immunity against pathogens. However, there are few or no comprehensive omics studies on this effect. From an evolutionary perspective, animals living in complex and changeable marine environments might develop signaling pathways to address bacterial threats under hypoxia. In this study, the ancient genomic model animal Takifugu obscurus and widespread Vibrio parahaemolyticus were utilized to study the effect. T. obscurus was challenged by V. parahaemolyticus or (and) exposed to hypoxia. The effects of hypoxia and infection were identified, and a theoretical model of the host critical signaling pathway in response to hypoxia and infection was defined by methods of comparative metabolomics and proteomics on the entire liver. The changing trends of some differential metabolites and proteins under hypoxia, infection or double stressors were consistent. The model includes transforming growth factor-β1 (TGF-β1), hypoxia-inducible factor-1α (HIF-1α), and epidermal growth factor (EGF) signaling pathways, and the consistent changing trends indicated that the host liver tended toward cell proliferation. Hypoxia and infection caused tissue damage and fibrosis in the portal area of the liver, which may be related to TGF-β1 signal transduction. We propose that LRG (leucine-rich alpha-2-glycoprotein) is widely involved in the transition of the TGF-β1/Smad signaling pathway in response to hypoxia and pathogenic infection in vertebrates as a conserved molecule.

Comparative transcriptomic analysis of the brain in Takifugu rubripes shows its tolerance to acute hypoxia

Hypoxia in water that caused by reduced levels of oxygen occurred frequently, due to the complex aquatic environment. Hypoxia tolerance for fish depends on a complete set of coping mechanisms such as oxygen perception and gene-protein interaction regulation. The present study examined the short-term effects of hypoxia on the brain in Takifugu rubripes. We sequenced the transcriptomes of the brain in T. rubripes to study their response mechanism to acute hypoxia. A total of 167 genes were differentially expressed in the brain of T. rubripes after exposed to acute hypoxia. Gene ontology and KEGG enrichment analysis indicated that hypoxia could cause metabolic and neurological changes, showing the clues of their adaptation to acute hypoxia. As the most complex and important organ, the brain of T. rubripes might be able to create a self-protection mechanism to resist or reduce damage caused by acute hypoxia stress.

Integrated Genomic Analyses From Low-Depth Sequencing Help Resolve Phylogenetic Incongruence in the Bamboos (Poaceae: Bambusoideae)

The bamboos (Bambusoideae, Poaceae) comprise a major grass lineage with a complex evolutionary history involving ancient hybridization and allopolyploidy. About 1700 described species are classified into three tribes, Olyreae (herbaceous bamboos), Bambuseae (tropical woody bamboos), and Arundinarieae (temperate woody bamboos). Nuclear analyses strongly support monophyly of the woody tribes, whereas plastome analyses strongly support paraphyly, with Bambuseae sister to Olyreae. Our objectives were to clarify the origin(s) of the woody bamboo tribes and resolve the nuclear vs. plastid conflict using genomic tools. For the first time, plastid and nuclear genomic information from the same bamboo species were combined in a single study. We sampled 51 species of bamboos representing the three tribes, estimated their genome sizes and generated low-depth sample sequence data, from which plastomes were assembled and nuclear repeats were analyzed. The distribution of repeat families was found to agree with nuclear gene phylogenies, but also provides novel insights into nuclear evolutionary history. We infer two early, independent hybridization events, one between an Olyreae ancestor and a woody ancestor giving rise to the two Bambuseae lineages, and another between two woody ancestors giving rise to the Arundinarieae. Retention of the Olyreae plastome associated with differential dominance of nuclear genomes and subsequent diploidization in some lineages explains the paraphyly observed in plastome phylogenetic estimations. We confirm ancient hybridization and allopolyploidy in the origins of the extant woody bamboo lineages and propose biased fractionation and diploidization as important factors in their evolution.

Response of gut microbiota to feed-borne bacteria depends on fish growth rate: a snapshot survey of farmed juvenile Takifugu obscurus

Environmental bacteria have a great impact on fish gut microbiota, yet little is known as to where fish acquire their gut symbionts, and how gut microbiota response to the disturbance from environmental bacteria. Through the integrative analysis by community profiling and source tracking, we show that feed‐associated bacteria can impose a strong disturbance upon the hindgut microbiota of cultured fugu. Consequently, marked alterations in the composition and function of gut microbiota in slow growth fugu were observed, implying a reduced stability upon bacterial disturbance from feed. Moreover, quantitative ecological analyses indicated that homogeneous selection and dispersal limitation largely contribute to the community stability and partial variations among hosts in the context of lower degree of disturbance. While the disturbance peaked, variable selection leads to an augmented interaction within gut microbiota, entailing community unstability and shift. Our findings emphasized the intricate linkage between feed and gut microbiota and highlighted the importance of resolving the feed source signal before the conclusion of comparative analysis of microbiota can be drawn. Our results provide a deeper insight into aquaculture of fugu and other economically important fishes and have further implications for an improved understanding of host–microbe interactions in the vertebrate gastrointestinal tract.

Conservative route to genome compaction in a miniature annelid

The causes and consequences of genome reduction in animals are unclear because our understanding of this process mostly relies on lineages with often exceptionally high rates of evolution. Here, we decode the compact 73.8-megabase genome of Dimorphilus gyrociliatus, a meiobenthic segmented worm. The D. gyrociliatus genome retains traits classically associated with larger and slower-evolving genomes, such as an ordered, intact Hox cluster, a generally conserved developmental toolkit and traces of ancestral bilaterian linkage. Unlike some other animals with small genomes, the analysis of the D. gyrociliatus epigenome revealed canonical features of genome regulation, excluding the presence of operons and trans-splicing. Instead, the gene-dense D. gyrociliatus genome presents a divergent Myc pathway, a key physiological regulator of growth, proliferation and genome stability in animals. Altogether, our results uncover a conservative route to genome compaction in annelids, reminiscent of that observed in the vertebrate Takifugu rubripes.

This study reports the genome of the miniature segmented annelid Dimorphilus gyrociliatus and reveals no drastic changes in genome architecture and regulation, unlike other cases of genome miniaturization.

Quantitative proteomic analysis in serum of Takifugu rubripes infected with Cryptocaryon irritans

Cryptocaryon irritans can cause cryptocaryonosis (white spot disease) in marine fish but the pathogenesis of the disease is unclear. In this work, we used high-throughput proteomics to identify differentially expressed proteins in the serum of Takifugu rubripes challenged with C. irritans. By using quantitative proteomic assays combined with Tandem Mass Tag-labeled quantitative proteomic analysis, we identified a total of 2088 differentially abundant proteins (1706 proteins were quantified, p < 0.05, fold-change threshold ≥ 2), including 21 up-regulated and 44 down-regulated. Combined with STRING-based functional analysis, we ultimately obtained eight proteins including glucokinase-like, integrin beta-1-like isoform X2, H4, H2A.V, histone H1-like, histone H2AX-like, histone H2B 1/2-like and myosin-9 isoform X1, which could be considered as potential biomarkers for T. rubripes immune responses. Eight proteins that were selected to validate significant differentially expressed genes at the proteomic level were consistent with qPCR at the transcriptomic level. The proteins identified in our work may serve as candidates for elucidating the molecular mechanism of cryptocaryonosis in T. rubripes. Our collective findings could provide new insights into searching for disease-specific targets and biomarkers, which may be effective indicators of C. irritans infection in T. rubripes.

Whole genome-wide chromosome fusion and new gene birth in the Monopterus albus genome

Background

Teleost fishes account for over half of extant vertebrate species. A core question in biology is how genomic changes drive phenotypic diversity that relates to the origin of teleost fishes.

Results

Here, we used comparative genomic analyses with chromosome assemblies of diverse lineages of vertebrates and reconstructed an ancestral vertebrate genome, which revealed phylogenomic trajectories in vertebrates. We found that the whole-genome-wide chromosome fission/fusions took place in the Monopterus albus lineage after the 3-round whole-genome duplication. Four times of genomic fission/fusions events resulted in the whole genome-wide chromosome fusions in the genomic history of the lineage. In addition, abundant recently evolved new genes for reproduction emerged in the Monopterus albus after separated from medaka. Notably, we described evolutionary trajectories of conserved blocks related to sex determination genes in teleosts.

Conclusions

These data pave the way for a better understanding of genomic evolution in extant teleosts.

A chromosome-scale assembly of the smallest Dothideomycete genome reveals a unique genome compaction mechanism in filamentous fungi

Background

The wide variation in the size of fungal genomes is well known, but the reasons for this size variation are less certain. Here, we present a chromosome-scale assembly of ectophytic Peltaster fructicola, a surface-dwelling extremophile, based on long-read DNA sequencing technology, to assess possible mechanisms associated with genome compaction.

Results

At 18.99 million bases (Mb), P. fructicola possesses one of the smallest known genomes sequence among filamentous fungi. The genome is highly compact relative to other fungi, with substantial reductions in repeat content, ribosomal DNA copies, tRNA gene quantity, and intron sizes, as well as intergenic lengths and the size of gene families. Transposons take up just 0.05% of the entire genome, and no full-length transposon was found. We concluded that reduced genome sizes in filamentous fungi such as P. fructicola, Taphrina deformans and Pneumocystis jirovecii occurred through reduction in ribosomal DNA copy number and reduced intron sizes. These dual mechanisms contrast with genome reduction in the yeast fungus Saccharomyces cerevisiae, whose small and compact genome is associated solely with intron loss.

Conclusions

Our results reveal a unique genomic compaction architecture of filamentous fungi inhabiting plant surfaces, and broaden the understanding of the mechanisms associated with compaction of fungal genomes.

Comparative genomic analysis of inbred rat strains reveals the existence of ancestral polymorphisms

In an alignment of closely related genomic sequences, the existence of discordant mutation sites, which do not reflect the phylogenetic relationship of the genomes, is often observed. Although these discordant mutation sites are thought to have emerged by ancestral polymorphism or gene flow, their frequency and distribution in the genome have not yet been analyzed in detail. Using the genome sequences of all protein coding genes of 25 inbred rat strains, we analyzed the frequency and genome-wide distribution of the discordant mutation sites. From the comparison of different substrains, it was found that these loci are not substrain specific, but are common among different groups of substrains, suggesting that the discordant sites might have mainly emerged through ancestral polymorphism. It was also revealed that the discordant sites are not uniformly distributed along chromosomes, but are concentrated at certain genomic loci, such as RT1, major histocompatibility complex of rats, and olfactory receptors, indicating that genes known to be highly polymorphic tend to have more discordant sites. Our results also showed that loci with a high density of discordant sites are also rich in heterozygous variants, even though these are inbred strains.

Construction of a High-Density Genetic Linkage Map and QTL Mapping for Growth-Related Traits in Takifugu bimaculatus

Takifugu bimaculatus is a euryhaline species, distributed ranging from the southern Yellow Sea to the South China Sea. Their tolerance to a wide range of salinity and temperature, coupled with a desirable firm texture, makes T. bimaculatus a strong candidate for Takifugu aquaculture in subtropics areas. Due to the increasing demand in markets and emerging of the Takifugu aquaculture industry, close attention has been paid to improvement on the T. bimaculatus production. In aquaculture, the great effort has been put into marker-assisted selective breeding, and efficient improvement was realized. However, few genetic resources on T. bimaculatus are provided so far. Aiming at understanding the genetic basis underlying important economic growth traits, facilitating genetic improvement and enriching the genetic resource in T. bimaculatus, we constructed the first genetic linkage map for T. bimaculatus via double digestion restriction-site association DNA sequencing and conducted quantitative traits locus (QTL) mapping for growth-related traits. The map comprised 1976 single nucleotide polymorphism markers distributed on 22 linkage groups (LG), with a total genetic distance of 2039.74 cM. Based on the linkage map, a chromosome-level assembly was constructed whereby we carried out comparative genomics analysis, verifying the high accuracy on contigs ordering of the linkage map. On the other hand, 18 QTLs associated with growth traits were detected on LG6, LG7, LG8, LG10, LG20, and LG21 with phenotypical variance ranging from 15.1 to 56.4%. Candidate genes participating in cartilage development, fat accumulation, and other growth-related regulation activities were identified from these QTLs, including col11a1, foxa2, and thrap3. The linkage map provided a solid foundation for chromosomes assembly and refinement. QTLs reported here unraveled the genomic architecture of some growth traits, which will advance the investigation of aquaculture breeding efforts in T. bimaculatus.

Chromosomal-level assembly of Takifugu obscurus (Abe, 1949) genome using third-generation DNA sequencing and Hi-C analysis

The Tetraodontidae family are known to have relatively small and compact genomes compared to other vertebrates. The obscure puffer fish Takifugu obscurus is an anadromous species that migrates to freshwater from the sea for spawning. Thus the euryhaline characteristics of T. obscurus have been investigated to gain understanding of their survival ability, osmoregulation, and other homeostatic mechanisms in both freshwater and seawater. In this study, a high quality chromosome-level reference genome for T. obscurus was constructed using long-read Pacific Biosciences (PacBio) Sequel sequencing and a Hi-C-based chromatin contact map platform. The final genome assembly of T. obscurus is 381 Mb, with a contig N50 length of 3,296 kb and longest length of 10.7 Mb, from a total of 62 Gb of raw reads generated using single-molecule real-time sequencing technology from a PacBio Sequel platform. The PacBio data were further clustered into chromosome-scale scaffolds using a Hi-C approach, resulting in a 373 Mb genome assembly with a contig N50 length of 15.2 Mb and and longest length of 28 Mb. When we directly compared the 22 longest scaffolds of T. obscurus to the 22 chromosomes of the tiger puffer Takifugu rubripes, a clear one-to-one orthologous relationship was observed between the two species, supporting the chromosome-level assembly of T. obscurus. This genome assembly can serve as a valuable genetic resource for exploring fugu-specific compact genome characteristics, and will provide essential genomic information for understanding molecular adaptations to salinity fluctuations and the evolution of osmoregulatory mechanisms.

Transcriptome-wide identification and characterization of the Sox gene family and microsatellites for Corbicula fluminea

The Asian clam, Corbicula fluminea, is a commonly consumed small freshwater bivalve in East Asia. However, available genetic information of this clam is still limited. In this study, the transcriptome of female C. fluminea was sequenced using the Illumina HiSeq 2500 platform. A total of 89,563 unigenes were assembled with an average length of 859 bp, and 36.7% of them were successfully annotated. Six members of Sox gene family namely SoxB1, SoxB2, SoxC, SoxD, SoxE and SoxF were identified. Based on these genes, the divergence time of C. fluminea was estimated to be around  476 million years ago. Furthermore, a total of 3,117 microsatellites were detected with a distribution density of 1:12,960 bp. Fifty of these microsatellites were randomly selected for validation, and 45 of them were successfully amplified with 31 polymorphic ones. The data obtained in this study will provide useful information for future genetic and genomic studies in C. fluminea.

A highly flexible and repeatable genotyping method for aquaculture studies based on target amplicon sequencing using next-generation sequencing technology

Studies using genome-wide single nucleotide polymorphisms (SNPs) have become commonplace in genetics and genomics, due to advances in high-throughput sequencing technologies. Since the numbers of required SNPs and samples vary depending on each research goal, genotyping technologies with high flexibility in the number of SNPs/samples and high repeatability have been intensively investigated. For example, the ultrahigh-multiplexed amplicon sequencing, Ion AmpliSeq, has been used as a high-throughput genotyping method mainly for diagnostic purposes. Here, we designed a custom panel targeting 3,187 genome-wide SNPs of fugu, Takifugu rubripes, and applied it for genotyping farmed fugu to test its feasibility in aquaculture studies. We sequenced two libraries consisting of different pools of individuals (n = 326 each) on the Illumina MiSeq sequencer. Consequently, over 99% target regions (3,178 SNPs) were amplified and 2,655 SNPs were available after filtering steps. Strong correlation was observed in the mean depth of coverage of each SNP between duplicate runs (r = 0.993). Genetic analysis using these genotype data successfully detected the known population structure and the sex determining locus of fugu. These results show the method is superior in repeatability and flexibility, and suits genetic studies including molecular breeding, such as marker assisted and genomic selection.

The Divergent Genomes of Teleosts

Boasting nearly 30,000 species, teleosts account for half of all extant vertebrates and approximately 98% of all ray-finned fish species (Actinopterygii). Teleosts are also the largest and most diverse group of vertebrates, exhibiting an astonishing level of morphological, physiological, and behavioral diversity. Previous studies had indicated that the teleost lineage has experienced an additional whole-genome duplication event. Recent comparative genomic analyses of teleosts and other bony vertebrates using spotted gar (a nonteleost ray-finned fish) and elephant shark (a cartilaginous fish) as outgroups have revealed several divergent features of teleost genomes. These include an accelerated evolutionary rate of protein-coding and nucleotide sequences, a higher rate of intron turnover, loss of many potential cis-regulatory elements and shorter conserved syntenic blocks. A combination of these divergent genomic features might have contributed to the evolution of the amazing phenotypic diversity and morphological innovations of teleosts.

Centromeric Satellite DNA in Flatfish (Order Pleuronectiformes) and Its Relation to Speciation Processes

Two new centromeric satellite DNAs in flatfish (Order Pleuronectiformes) have been characterized. The SacI-family from Hippoglossus hippoglossus, restricted to this species, had a monomeric size of 334 base pair (bp) and was located in most of the centromeres of its karyotype. The PvuII-family, with a monomeric size of 177 bp, was initially isolated from the genome of Solea senegalensis, and fluorescent in situ hybridization (FISH) localized the repeat to centromeres of most of the chromosomes. This family could only be amplified in 2 other species of the genus Solea (Solea solea and Solea lascaris). Molecular features and chromosomal location indicated a possible structural and/or functional role of these sequence repeats. The presence of species-specific satellite-DNA families in the centromeres and their possible role in the speciation processes in this group of fishes is discussed.

Relationships between Gene Structure and Genome Instability in Flowering Plants

Flowering plant (angiosperm) genomes are exceptional in their variability with respect to genome size, ploidy, chromosome number, gene content, and gene arrangement. Gene movement, although observed in some of the earliest plant genome comparisons, has been relatively underinvestigated. We present herein a description of several interesting properties of plant gene and genome structure that are pertinent to the successful movement of a gene to a new location. These considerations lead us to propose a model that can explain the frequent success of plant gene mobility, namely that Small Insulated Genes Move Around (SIGMAR). The SIGMAR model is then compared with known processes for gene mobilization, and predictions of the SIGMAR model are formulated to encourage future experimentation. The overall results indicate that the frequent gene movement in angiosperm genomes is partly an outcome of the unusual properties of angiosperm genes, especially their small size and insulation from epigenetic silencing.

Ultrahigh-Density Linkage Map Construction Using Low-Coverage Whole-Genome Sequencing of a Doubled Haploid Population: Case Study of Torafugu (Takifugu rubripes)

Next-generation sequencing enables genome-wide genotyping of a large population and further facilitates the construction of a genetic linkage map. Low-coverage whole-genome sequencing has been employed for genetic linkage map construction in several species. However, this strategy generally requires available high-quality reference genomes and/or designed inbred pedigree lines, which restrict the scope of application for non-model and unsequenced species. Here, using torafugu (Takifugu rubripes) as a test model, we propose a new strategy for ultrahigh-density genetic linkage map construction using low-coverage whole-genome sequencing of a haploid/doubled haploid (H/DH) population without above requirements. Low-coverage (≈1×) whole-genome sequencing data of 165 DH individuals were used for de novo assembly and further performed single nucleotide polymorphisms (SNPs) calling, resulting in the identification of 1,070,601 SNPs. Based on SNP genotypes and de novo assembly, genotypes were associated with short DNA segments and an ultrahigh-density linkage map was constructed containing information of 802,277 SNPs in 3090 unique positions. Comparative analyses showed near-perfect concordance between the present linkage map and the latest published torafugu genome (FUGU5). This strategy would facilitate ultrahigh-density linkage map construction in various sexually reproducing organisms for which H/DH populations can be generated.

Production of Tiger Puffer Takifugu rubripes Offspring from Triploid Grass Puffer Takifugu niphobles Parents

The tiger puffer Takifugu rubripes is one of the most popular aquacultural fish; however, there are two major obstacles to selective breeding. First, they have a long generation time of 2 or 3 years until maturation. Second, the parental tiger puffer has a body size (2-5 kg) much larger than average market size (0.6-1.0 kg). The grass puffer Takifugu niphobles is closely related to the tiger puffer and matures in half the time. Furthermore, grass puffer can be reared in small areas since their maturation weight is about 1/150 that of mature tiger puffer. Therefore, to overcome the obstacles of maturation size and generation time of tiger puffer, we generated surrogate grass puffer that can produce tiger puffer gametes through germ cell transplantation. Approximately 5000 tiger puffer testicular cells were transplanted into the peritoneal cavity of triploid grass puffer larvae at 1 day post hatching. When the recipient fish matured, both males and females produced donor-derived gametes. Through their insemination, we successfully produced donor-derived tiger puffer offspring presenting the same body surface dot pattern, number of dorsal fin rays, and DNA fingerprint as those of the donor tiger puffer, suggesting that the recipient grass puffer produced functional eggs and sperm derived from the donor tiger puffer. Although fine tunings are still needed to improve efficiencies, surrogate grass puffer are expected to accelerate the breeding process of tiger puffer because of their short generation time and small body size.

Decontaminating eukaryotic genome assemblies with machine learning

BACKGROUND

High-throughput sequencing has made it theoretically possible to obtain high-quality de novo assembled genome sequences but in practice DNA extracts are often contaminated with sequences from other organisms. Currently, there are few existing methods for rigorously decontaminating eukaryotic assemblies. Those that do exist filter sequences based on nucleotide similarity to contaminants and risk eliminating sequences from the target organism.

RESULTS

We introduce a novel application of an established machine learning method, a decision tree, that can rigorously classify sequences. The major strength of the decision tree is that it can take any measured feature as input and does not require a priori identification of significant descriptors. We use the decision tree to classify de novo assembled sequences and compare the method to published protocols.

CONCLUSIONS

A decision tree performs better than existing methods when classifying sequences in eukaryotic de novo assemblies. It is efficient, readily implemented, and accurately identifies target and contaminant sequences. Importantly, a decision tree can be used to classify sequences according to measured descriptors and has potentially many uses in distilling biological datasets.

Embryonic development of the grass pufferfish (Takifugu niphobles): From egg to larvae

Tetraodontidae (pufferfish) family members carry the smallest genomes among vertebrates, and these pocket-sized genomes have directly contributed to our understanding of the structure and evolution of higher animals. The grass pufferfish (Takifugu niphobles) could be considered a potential new model organism for comparative genomics and development due to the potential access to embryos, and availability of sequence data for two similar genomes: that of spotted green pufferfish (Tetraodon nigroviridis) and Fugu (Takifugu rubripes). In this study, we provide the first description of the normal embryonic development of T. niphobles, by drawing comparisons with the closely related species cited above. Embryos were obtained by in vitro fertilization of eggs, and subsequent development was monitored at a constant temperature consistent with natural conditions. T. niphobles development was divided into seven periods of embryogenesis: the zygote, cleavage, blastula, gastrula, segmentation, pharyngula, and hatching periods; and stages subdividing these periods are defined based on morphological characteristics. The developmental stage series described in this study aims to provide the utilization of T. niphobles as an experimental model organism for comparative developmental studies.

New insights on the sister lineage of percomorph fishes with an anchored hybrid enrichment dataset

Percomorph fishes represent over 17,100 species, including several model organisms and species of economic importance. Despite continuous advances in the resolution of the percomorph Tree of Life, resolution of the sister lineage to Percomorpha remains inconsistent but restricted to a small number of candidate lineages. Here we use an anchored hybrid enrichment (AHE) dataset of 132 loci with over 99,000 base pairs to identify the sister lineage of percomorph fishes. Initial analyses of this dataset failed to recover a strongly supported sister clade to Percomorpha, however, scrutiny of the AHE dataset revealed a bias towards high GC content at fast-evolving codon partitions (GC bias). By combining several existing approaches aimed at mitigating the impacts of convergence in GC bias, including RY coding and analyses of amino acids, we consistently recovered a strongly supported clade comprised of Holocentridae (squirrelfishes), Berycidae (Alfonsinos), Melamphaidae (bigscale fishes), Cetomimidae (flabby whalefishes), and Rondeletiidae (redmouth whalefishes) as the sister lineage to Percomorpha. Additionally, implementing phylogenetic informativeness (PI) based metrics as a filtration method yielded this same topology, suggesting PI based approaches will preferentially filter these fast-evolving regions and act in a manner consistent with other phylogenetic approaches aimed at mitigating GC bias. Our results provide a new perspective on a key issue for studies investigating the evolutionary history of more than one quarter of all living species of vertebrates.

De novo assembly, characterization and annotation for the transcriptome of Sarcocheilichthys sinensis

The Chinese lake gudgeon Sarcocheilichthys sinensis is a small cyprinid fish with great aquaculture potential both for its edible and ornamental values. Nevertheless, available genomic and transcriptomic information for this fish is extremely deficient. In this study, a normalized cDNA library was constructed using 13 mixed tissues of an adult male S. sinensis, and was sequenced by the Illumina HiSeq2500 platform. De novo assembly was performed using 38,911,511 obtained clean reads, and a total of 147,282 unigenes with an average length of 900 bp were finally achieved. 96.2% of these unigenes were annotated in 9 public databases, and 16 segments of growth-related genes were identified for future studies. In addition, 28,493 unigenes were assigned to 61 subcategories of Gene Ontology (GO), and 10,483 unigenes were assigned to 25 categories of Cluster of Orthologous Group (COG). Moreover, 14,943 unigenes were classified into 225 pathways of the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. A total of 30,666 microsatellites were detected from 17,627 unigenes with an average distribution density of 1:2405 bp. This transcriptome data set will be valuable for researches on discovery, expression and evolution on genes of interest. Meanwhile, the identified microsatellites would be useful tools for genetic and genomic studies in S. sinensis.

The Tetraodon nigroviridis reference transcriptome: developmental transition, length retention and microsynteny of long non-coding RNAs in a compact vertebrate genome

Pufferfish such as fugu and tetraodon carry the smallest genomes among all vertebrates and are ideal for studying genome evolution. However, comparative genomics using these species is hindered by the poor annotation of their genomes. We performed RNA sequencing during key stages of maternal to zygotic transition of Tetraodon nigroviridis and report its first developmental transcriptome. We assembled 61,033 transcripts (23,837 loci) representing 80% of the annotated gene models and 3816 novel coding transcripts from 2667 loci. We demonstrate the similarities of gene expression profiles between pufferfish and zebrafish during maternal to zygotic transition and annotated 1120 long non-coding RNAs (lncRNAs) many of which differentially expressed during development. The promoters for 60% of the assembled transcripts result validated by CAGE-seq. Despite the extreme compaction of the tetraodon genome and the dramatic loss of transposons, the length of lncRNA exons remain comparable to that of other vertebrates and a small set of lncRNAs appears enriched for transposable elements suggesting a selective pressure acting on lncRNAs length and composition. Finally, a set of lncRNAs are microsyntenic between teleost and vertebrates, which indicates potential regulatory interactions between lncRNAs and their flanking coding genes. Our work provides a fundamental molecular resource for vertebrate comparative genomics and embryogenesis studies.

Identification and expression of a novel carbonic anhydrase isozyme in the pufferfish Takifugu vermicularis

Carbonic anhydrase (CA) is a key element for maintaining acid base balance in fish. In our present experiment, novel CA isozymes were identified from the pear puffer (Takifugu vermicularis). Based on the high homology of two predicted CA sequences of the tiger puffer (Takifugu rubripes), a 1715bp novel cDNA was obtained from T. vermicularis. The open reading frame showed a complete coding sequence of 552bp with a deduced peptide sequence of 183 amino acids that exhibited highest (97%) identity with pufferfish putative CA III and CA IV-like sequences. In addition, this translated protein sequence showed 36-37% identity with zebrafish CA IV-like, CA XVa, CA XVb, and CA XVc proteins. Phylogenetic analysis revealed that the pufferfish novel protein (pCAn) was a membrane-bound CA protein. Alignment of multiple CA sequences illustrated that most of the putative active site residues of the pCAn isozyme were situated at highly conserved regions of the CA sequences. Examination of motif distribution suggested that the pCAn isozyme was very similar to the puffer predicted CA IV-like isozyme. Reverse transcription-polymerase chain reaction (PCR) analysis showed highly differential expression in the brain, gills, kidney, and muscle, whereas CA mRNA expression was almost absent in heart, liver, and intestine. Quantitative PCR expression of CA mRNA abundance suggested several-fold higher expression of pCAn isozymes in the gills compared to other tissues tested. Our results suggest that the pCAn isozyme might be related to CA IV-like isozymes. Further functional studies are needed to investigate the function of the pCAn isozyme in T. vermicularis.

Model Organisms for Studying the Cell Cycle

Regulation of the cell-division cycle is fundamental for the growth, development, and reproduction of all species of life. In the past several decades, a conserved theme of cell cycle regulation has emerged from research in diverse model organisms. A comparison of distinct features of several diverse model organisms commonly used in cell cycle studies highlights their suitability for various experimental approaches, and recaptures their contributions to our current understanding of the eukaryotic cell cycle. A historic perspective presents a recollection of the breakthrough upon unfolding the universal principles of cell cycle control by scientists working with diverse model organisms, thereby appreciating the discovery pathways in this field. A comprehensive understanding is necessary to address current challenging questions about cell cycle control. Advances in genomics, proteomics, quantitative methodologies, and approaches of systems biology are redefining the traditional concept of what constitutes a model organism and have established a new era for development of novel, and refinement of the established model organisms. Researchers working in the field are no longer separated by their favorite model organisms; they have become more integrated into a larger community for gaining greater insights into how a cell divides and cycles. The new technologies provide a broad evolutionary spectrum of the cell-division cycle and allow informative comparisons among different species at a level that has never been possible, exerting unimaginable impact on our comprehensive understanding of cell cycle regulation.

Comparative Analysis of WRKY Genes Potentially Involved in Salt Stress Responses in Triticum turgidum L. ssp. durum

WRKY transcription factors are involved in multiple aspects of plant growth, development and responses to biotic stresses. Although they have been found to play roles in regulating plant responses to environmental stresses, these roles still need to be explored, especially those pertaining to crops. Durum wheat is the second most widely produced cereal in the world. Complex, large and unsequenced genomes, in addition to a lack of genomic resources, hinder the molecular characterization of tolerance mechanisms. This paper describes the isolation and characterization of five TdWRKY genes from durum wheat (Triticum turgidum L. ssp. durum). A PCR-based screening of a T. turgidum BAC genomic library using primers within the conserved region of WRKY genes resulted in the isolation of five BAC clones. Following sequencing fully the five BACs, fine annotation through Triannot pipeline revealed 74.6% of the entire sequences as transposable elements and a 3.2% gene content with genes organized as islands within oceans of TEs. Each BAC clone harbored a TdWRKY gene. The study showed a very extensive conservation of genomic structure between TdWRKYs and their orthologs from Brachypodium, barley, and T. aestivum. The structural features of TdWRKY proteins suggested that they are novel members of the WRKY family in durum wheat. TdWRKY1/2/4, TdWRKY3, and TdWRKY5 belong to the group Ia, IIa, and IIc, respectively. Enrichment of cis-regulatory elements related to stress responses in the promoters of some TdWRKY genes indicated their potential roles in mediating plant responses to a wide variety of environmental stresses. TdWRKY genes displayed different expression patterns in response to salt stress that distinguishes two durum wheat genotypes with contrasting salt stress tolerance phenotypes. TdWRKY genes tended to react earlier with a down-regulation in sensitive genotype leaves and with an up-regulation in tolerant genotype leaves. The TdWRKY transcripts levels in roots increased in tolerant genotype compared to sensitive genotype. The present results indicate that these genes might play some functional role in the salt tolerance in durum wheat.

Micromanaging metabolism-a role for miRNAs in teleost energy metabolism

MicroRNAs (miRNAs) are small, non-protein coding RNA sequences, which are found in most eukaryotes. Since their initial discovery, miRNAs have emerged as important regulators of many biological processes. One of the most important processes profoundly regulated by miRNAs is energy metabolism. Traditionally, metabolic functions of miRNAs have been studied in genome-sequenced mammalian organisms, especially the mouse model. However, partially driven by commercial interest in aquaculture, increasingly feasible large-scale molecular techniques have resulted in the characterization of miRNA repertoires, and importantly, several genome sequences of several (commercially important) teleost species, which also hold important roles as research models in the comparative physiology of energy metabolism. This review aims to introduce the recent advances in miRNA research in teleost fish and to describe the current knowledge of miRNA function in teleost energy metabolism. The most pressing research needs and questions to determine metabolic roles of miRNAs in teleost models are presented, as well as applicable technical approaches and current bottlenecks. Rainbow trout, which possess the advantages of newly available molecular tools and a long history as comparative research model in teleost energy metabolism, are discussed as a promising research model to address these questions.

Acytota - associated kingdom of neglected life

There is a huge variety of RNA- and DNA-containing entities that multiply within and propagate between cells across all kingdoms of life, having no cells of their own. Apart from cellular organisms, these entities (viroids, plasmids, mobile elements and viruses among others) are the only ones with distinct genetic identities but which are not included in any traditional tree of life. We suggest to introduce or, rather, revive the distinct category of acellular organisms, Acytota, as an additional, undeservedly ignored full-fledged kingdom of life. Acytota are indispensable players in cellular life and its evolution. The six traditional kingdoms (Cytota) and Acytota together complete the classification of the biological world (Biota), leaving nothing beyond.

Toward understanding the evolution of vertebrate gene regulatory networks: comparative genomics and epigenomic approaches

Vertebrates, as most animal phyla, originated >500 million years ago during the Cambrian explosion, and progressively radiated into the extant classes. Inferring the evolutionary history of the group requires understanding the architecture of the developmental programs that constrain the vertebrate anatomy. Here, I review recent comparative genomic and epigenomic studies, based on ChIP-seq and chromatin accessibility, which focus on the identification of functionally equivalent cis-regulatory modules among species. This pioneer work, primarily centered in the mammalian lineage, has set the groundwork for further studies in representative vertebrate and chordate species. Mapping of active regulatory regions across lineages will shed new light on the evolutionary forces stabilizing ancestral developmental programs, as well as allowing their variation to sustain morphological adaptations on the inherited vertebrate body plan.

Description of embryonic development of spotted green pufferfish (Tetraodon nigroviridis)

Pufferfish species of the Tetraodontidae family carry the smallest genomes among vertebrates. Their compressed genomes are thought to be enriched for functional DNA compared to larger vertebrate genomes, and they are important models for comparative genomics. The significance of pufferfish as model organisms in comparative genomics is due to the availability of two sequenced genomes, that of spotted green pufferfish (Tetraodon nigroviridis) and fugu (Takifugu rubripes). However, there is only a very limited utilization of pufferfish as an experimental model organism, due to the lack of established husbandry and developmental genetics protocols. In this study, we provide the first description of the normal embryonic development of Tetraodon nigroviridis. Embryos were obtained by in vitro fertilization of eggs, and subsequent development was monitored by brightfield microscopy at constant temperature. Tetraodon development was divided into distinct stages based on diagnostic morphological features, which were adopted from published literature on normal development of other fish species like medaka (Oryzias latipes), zebrafish (Danio rerio), and fugu. Tetraodon embryos show more similar morphologies to medaka than to zebrafish, reflecting its phylogenetic position. The early developmental stage series described in this study forms the foundation for the utilization of tetraodon as an experimental model organism for comparative developmental studies.

Genetic Basis Underlying Behavioral Correlation Between Fugu Takifugu rubripes and a Closely Related Species, Takifugu niphobles

Correlated suits of behaviors (behavioral syndrome) are commonly observed in both inter- and intraspecific studies. In order to understand the genetic basis of such a correlation between species, we compared ten behaviors classified into five categories (acclimation, feeding, normal swimming, reaction to a novel object and activity in a novel environment) between two pufferfish species, Takifugu rubripes and T. niphobles. The two species showed consistent differences in nine behaviors with a significant correlation among behaviors. Quantitative trait locus (QTL) analysis using second generation hybrids revealed that different sets of small effect QTL are associated with the observed interspecific behavioral disparity. This indicates that correlations in temperament traits between them are governed by many genes with small effects, and each behavior has been selected to form particular combination patterns. One of the QTL showing small pleiotropic effect includes the Drd4 gene known for its association with behavioral traits in some animal taxa including mammals.

Characterization of the torafugu (Takifugu rubripes) immunoglobulin heavy chain gene locus

In this study, we investigated the immunoglobulin heavy (IGH) gene locus of torafugu (Takifugu rubripes) from publicly available assembly sequences and presented an annotated locus map, including the IGHV genes, pseudogenes, and IGHC genes. Three new IGHV gene families (IGHV3-IGHV5) were discovered. We observed the interspersion of IGHV1 and IGHV2 family members and that they often intermingled with each other, while other family members were further interspersed. Conservation of the promoter and recombination signal sequences (RSS) was observed in a family-specific manner. In addition to known variable region genes present on chromosome 5 (current torafugu genome assembly), we found 34 additional IGHV genes on scaffold 287 and three novel potentially functional IGHD genes on scaffold 483. In total, the variable region of the torafugu IGH locus consists of at least 48 IGHV genes, seven IGHD genes, and six IGHJ genes. IGHC genes have also been mapped in this study, with three genes encoding immunoglobulin classes: IgT, IgM, and IgD. We confirmed the expression of newly identified IGHV3 family sequences in the spleen and kidney of adult torafugu and found a favorable IGHV segment usage by IgM and IgT. Possible structural variation in the IGHδ locus was observed based on the current torafugu assembly. The complete characterization of the torafugu IGH locus will facilitate detailed studies of large-scale mechanisms associated with the recombination of the variable region genes and will offer insights into the genetic basis of the potential diversity in the antibody response observed in torafugu.