Members of the teleost family Syngnathidae (seahorses, pipefishes and seadragons) (Extended Data Fig. 1), comprising approximately 300 species, display a complex array of morphological innovations and reproductive behaviours. This includes specialized morphological phenotypes such as an elongated snout with a small terminal mouth, fused jaws, absent pelvic and caudal fins, and an extended body covered with an armour of bony plates instead of scales 1 (Fig. 1a). Syngnathids are also unique among vertebrates due to their 'male pregnancy' , whereby males nourish developing embryos in a brood pouch until hatching and parturition occurs 2,3 . In addition, members of the subfamily Hippocampinae (seahorses) exhibit other derived features such as the lack of a caudal fin, a characteristic prehensile tail, and a vertical body axis 4 (Fig. 1a). To understand the genetic basis of the specialized morphology and reproductive system of seahorses, we sequenced the genome of the tiger tail seahorse, H. comes, and carried out comparative genomic analyses with the genome sequences of other ray-finned fishes (Actinopterygii). Genome assembly and annotationThe genome of a male H. comes individual was sequenced using the Illumina HiSeq 2000 platform. After filtering low-quality and duplicate reads, 132.13 Gb (approximately 190-fold coverage of the estimated 695 Mb genome) of reads from libraries with insert sizes ranging from 170 bp to 20 kb were retained for assembly. The filtered reads were assembled using SOAPdenovo (version 2.04) to yield a 501.6 Mb assembly with an N50 contig size and N50 scaffold size of 34.7 kb and 1.8 Mb, respectively. Total RNA from combined soft tissues of H. comes was sequenced using RNA-sequencing (RNA-seq) and assembled de novo. The H. comes genome assembly is of high quality, as > 99% of the de novo assembled transcripts (76,757 out of 77,040) could be mapped to the assembly; and 243 out of 248 core eukaryotic genes mapping approach (CEGMA) genes are complete in the assembly.We predicted 23,458 genes in the genome of H. comes based on homology and by mapping the RNA-seq data of H. comes and a closely related species, the lined seahorse, Hippocampus erectus, to the genome assembly (see Methods and Supplementary Information). More than 97% of the predicted genes (22,941 genes) either have homologues in public databases (Swissprot, Trembl and the Kyoto Encyclopedia of Genes and Genomes (KEGG)) or are supported by assembled RNAseq transcripts. Analysis of gene family evolution using a maximum likelihood framework identified an expansion of 25 gene families (261 genes; 1.11%) and contraction of 54 families (96 genes; 0.41%) in the H. comes lineage (Extended Data Fig. 2 and Supplementary Tables 4.1, 4.2). Transposable elements comprise around 24.8% (124.5 Mb) of the H. comes genome, with class II DNA transposons being the most abundant class (9%; 45 Mb). Only one wave of transposable element expansion was identified, with no evidence for a recent transposable element burst (Kimura divergence ≤ 5) (Extended D...

BackgroundBlunt snout bream (Megalobrama amblycephala) is an herbivorous freshwater fish species native to China and has been recognized as a main aquaculture species in the Chinese freshwater polyculture system with high economic value. Right now, only limited EST resources were available for M. amblycephala. Recent advances in large-scale RNA sequencing provide a fast, cost-effective, and reliable approach to generate large expression datasets for functional genomic analysis, which is especially suitable for non-model species with un-sequenced genomes.Methodology and Principal FindingsUsing 454 pyrosequencing, a total of 1,409,706 high quality reads (total length 577 Mbp) were generated from the normalized cDNA of pooled M. amblycephala individuals. These sequences were assembled into 26,802 contigs and 73,675 singletons. After BLAST searches against the NCBI non-redundant (NR) and UniProt databases with an arbitrary expectation value of E−10, over 40,000 unigenes were functionally annotated and classified using the FunCat functional annotation scheme. A comparative genomics approach revealed a substantial proportion of genes expressed in M. amblycephala tanscriptome to be shared across the genomes of zebrafish, medaka, tetraodon, fugu, stickleback, human, mouse, and chicken, and identified a substantial number of potentially novel M. amblycephala genes. A total number of 4,952 SSRs were found and 116 polymorphic loci have been characterized. A significant number of SNPs (25,697) and indels (23,287) were identified based on specific filter criteria in the M. amblycephala.ConclusionsThis study is the first comprehensive transcriptome analysis for a fish species belonging to the genus Megalobrama. These large EST resources are expected to be valuable for the development of molecular markers, construction of gene-based linkage map, and large-scale expression analysis of M. amblycephala, as well as comparative genome analysis for the genus Megalobrama fish species. The identified SSR and SNP markers will greatly benefit its breeding program and whole genome association studies.

A high-density genetic map constructed with a wheat 55 K SNP array was highly consistent with the physical map of this species and it facilitated the identification of a novel major QTL for productive tiller number. Productive tiller number (PTN) plays a key role in wheat grain yield. In this study, a recombinant inbred line population with 199 lines derived from a cross between '20828' and 'Chuannong16' was used to construct a high-density genetic map using wheat 55 K single nucleotide polymorphism (SNP) array. The constructed genetic map contains 12,109 SNP markers spanning 3021.04 cM across the 21 wheat chromosomes. The orders of the genetic and physical positions of these markers are generally in agreement, and they also match well with those based on the 660 K SNP array from which the one used in this study was derived. The ratios of SNPs located in each of the wheat deletion bins were similar among the wheat 9 K, 55 K, 90 K, 660 K and 820 K SNP arrays. Based on the constructed maps, a novel major quantitative trait locus QPtn.sau-4B for PTN was detected across multi-environments in a 0.55 cM interval on 4B and it explained 17.23-45.46% of the phenotypic variance. Twenty common genes in the physical interval between the flanking markers were identified on chromosome 4B of 'Chinese Spring' and wild emmer. These results indicate that wheat 55 K SNP array could be an ideal tool in primary mapping of target genes and the identification of QPtn.sau-4B laid a foundation for the following fine mapping and cloning work.

BackgroundBlunt snout bream (Megalobrama amblycephala) is an economically important fish species in the Chinese freshwater polyculture system for its delicacy and high economic value. MicroRNAs (miRNAs) play important roles in regulation of almost all biological processes in eukaryotes. Although previous studies have identified thousands of miRNAs from many species, little information is known for miRNAs of M. amblycephala. To investigate functions of miRNAs associated with growth of M. amblycephala, we adopted the Solexa sequencing technology to sequence two small RNA libraries prepared from four growth related tissues (brain, pituitary, liver and muscle) of M. amblycephala using individuals with relatively high and low growth rates.ResultsIn this study, we have identified 347 conserved miRNAs (belonging to 123 families) and 22 novel miRNAs in M. amblycephala. Moreover, we observed sequence variants and seed edits of the miRNAs. Of the 5,166 single nucleotide substitutions observed in two libraries, the most abundant were G-to-U (15.9%), followed by U-to-C (12.1%), G-to-A (11.2%), and A to G (11.2%). Subsequently, we compared the expression patterns of miRNAs in the two libraries (big-size group with high growth rate versus small-size group with low growth rate). Results indicated that 27 miRNAs displayed significant differential expressions between the two libraries (p < 0.05). Of these, 16 were significantly up-regulated and 11 were significantly down-regulated in the big-size group compared to the small-size group. Furthermore, stem-loop RT-PCR was applied to validate and profile the expression of the differentially expressed miRNAs in ten tissues, and the result revealed that the conserved miRNAs expressed at higher levels than the novel miRNAs, especially in brain, liver and muscle. Also, targets prediction of differentially expressed miRNAs and KEGG pathway analysis suggested that differentially expressed miRNAs are involved in growth and metabolism, signal transduction, cell cycle, neural development and functions.ConclusionsThe present study provides the first large-scale characterization of miRNAs in M. amblycephala and miRNA profile related to different growth performances. The discovery of miRNA resource from this study is expected to contribute to a better understanding of the miRNAs roles playing in regulating the growth biological processes and the study of miRNA function and phenotype-associated miRNA identification in fish.

The blunt snout bream Megalobrama amblycephala is the economically most important cyprinid fish species. As an herbivore, it can be grown by eco-friendly and resource-conserving aquaculture. However, the large number of intermuscular bones in the trunk musculature is adverse to fish meat processing and consumption. As a first towards optimizing this aquatic livestock, we present a 1.116-Gb draft genome of M. amblycephala, with 779.54 Mb anchored on 24 linkage groups. Integrating spatiotemporal transcriptome analyses, we show that intermuscular bone is formed in the more basal teleosts by intramembranous ossification and may be involved in muscle contractibility and coordinating cellular events. Comparative analysis revealed that olfactory receptor genes, especially of the beta type, underwent an extensive expansion in herbivorous cyprinids, whereas the gene for the umami receptor T1R1 was specifically lost in M. amblycephala. The composition of gut microflora, which contributes to the herbivorous adaptation of M. amblycephala, was found to be similar to that of other herbivores. As a valuable resource for the improvement of M. amblycephala livestock, the draft genome sequence offers new insights into the development of intermuscular bone and herbivorous adaptation.

No abstract

The MADS-box genes encode transcription factors with key roles in plant growth and development. A comprehensive analysis of the MADS-box gene family in bread wheat (Triticum aestivum) has not yet been conducted, and our understanding of their roles in stress is rather limited. Here, we report the identification and characterization of the MADS-box gene family in wheat. A total of 180 MADS-box genes classified as 32 Mα, 5 Mγ, 5 Mδ, and 138 MIKC types were identified. Evolutionary analysis of the orthologs among T. urartu, Aegilops tauschii and wheat as well as homeologous sequences analysis among the three sub-genomes in wheat revealed that gene loss and chromosomal rearrangements occurred during and/or after the origin of bread wheat. Forty wheat MADS-box genes that were expressed throughout the investigated tissues and development stages were identified. The genes that were regulated in response to both abiotic stresses (i.e., phosphorus deficiency, drought, heat, and combined drought and heat) and biotic stresses (i.e., Fusarium graminearum, Septoria tritici, stripe rust and powdery mildew) were detected as well. A few notable MADS-box genes were specifically expressed in a single tissue and those showed relatively higher expression differences between the stress and control treatment. The expression patterns of considerable MADS-box genes differed from those of their orthologs in Brachypodium, rice, and Arabidopsis. Collectively, the present study provides new insights into the possible roles of MADS-box genes in response to stresses and will be valuable for further functional studies of important candidate MADS-box genes.

Background: The lined seahorse, Hippocampus erectus, is an Atlantic species and mainly inhabits shallow sea beds or coral reefs. It has become very popular in China for its wide use in traditional Chinese medicine. In order to improve the aquaculture yield of this valuable fish species, we are trying to develop genomic resources for assistant selection in genetic breeding. Here, we provide whole genome sequencing, assembly, and gene annotation of the lined seahorse, which can enrich genome resource and further application for its molecular breeding. Findings: A total of 174.6 Gb (Gigabase) raw DNA sequences were generated by the Illumina Hiseq2500 platform. The final assembly of the lined seahorse genome is around 458 Mb, representing 94% of the estimated genome size (489 Mb by k-mer analysis). The contig N50 and scaffold N50 reached 14.57 kb and 1.97 Mb, respectively. Quality of the assembled genome was assessed by BUSCO with prediction of 85% of the known vertebrate genes and evaluated using the de novo assembled RNA-seq transcripts to prove a high mapping ratio (more than 99% transcripts could be mapped to the assembly). Using homology-based, de novo and transcriptome-based prediction methods, we predicted 20 788 protein-coding genes in the generated assembly, which is less than our previously reported gene number (23 458) of the tiger tail seahorse (H. comes). Conclusion: We report a draft genome of the lined seahorse. These generated genomic data are going to enrich genome resource of this economically important fish, and also provide insights into the genetic mechanisms of its iconic morphology and male pregnancy behavior.