A panel of polymorphic EST-derived microsatellite loci for the small yellow croaker (Larimichthys polyactis)

Small yellow croaker (Larimichthys polyactis) is an important economic species of marine fishery. We developed and evaluated simple sequence repeat (SSR) markers from expressed sequence tags (ESTs) of Pseudosciaena crocea, Paralichthys olivaceus and Psetta maxima. Characteristics of nine EST–SSR loci were investigated using 46 L. polyactis individuals. The number of alleles per locus ranged from 2 to 6. The observed heterozygosity ranged from 0.0652 to 0.7391, while the expected heterozygosity ranged from 0.0638 to 0.7754. Seven loci departed from Hardy–Weinberg equilibrium (P < 0.01) significantly. These loci and markers will be useful for population genetics and systemic evolution among species of small yellow croaker.     

Complete mitochondrial genome sequence of bighead croaker Collichthys niveatus (Perciformes, Sciaenidae): a mitogenomic perspective on the phylogenetic relationships of Pseudosciaeniae

The monophyly and phylogenetic relationships of Pseudosciaeniae have long been controversial. Here we describe the mitochondrial genome (mitogenome) sequence of Collichthys niveatus. It is a circular double-stranded DNA molecule of 16,450 base pairs (bp) in length with a standard set of 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), 13 protein-coding genes as well as a non-coding control region. The mitogenome of C. niveatus shared common features with those of other bony fishes in terms of gene arrangement, base composition, and tRNA structures. The C. niveatus mitogenome exhibited pronounced strand-specific asymmetry in nucleotide composition, which was also reflected in the codon usage of genes oriented in opposite directions. Contrary to the typical structure of the control region, the central conserved blocks (CSB-D, -E, and -F) could not be detected in C. niveatus mitogenome. Phylogenetic analysis based on whole mitogenome sequences provided strong support for the monophyly of Pseudosciaeniae, and sister-group relationships of C. niveatus + Collichthys lucidus and Larimichthys crocea + Larimichthys polyactis, which was consistent with the traditional taxonomy. Unexpected divergence was found in two C. niveatus mitogenomes and several hypotheses were proposed to explain this observation including misidentification and introgressive hybridization between C. niveatus and L. polyactis, and polyphyletic origin of C. niveatus. We considered species misidentification to be the main hypothesis. However, additional data is essential to test these proposed hypotheses.     

Identification of a large SNP dataset in Larimichthys crocea using specific‐locus amplified fragment sequencing

The large yellow croaker, Larimichthys crocea, is a commercially important drum fish (Family: Sciaenidae) native to the East and South China Sea. Habitat deterioration and overfishing have led to significant population decline and the collapse of its fishery over the past decades. Today, the market supply of L. crocea depends solely on stocks produced in hatcheries and farms. Common issues that occur in the culture of L. crocea include germplasm degradation, precocious puberty, elevated disease susceptibility and growth retardation. In this study, we employed SLAF‐seq (specific‐locus amplified fragment sequencing) technology to identify single nucleotide polymorphism (SNP) loci across the L. crocea genome. Sixty samples were selected for SLAF analysis out of 1000 progeny in the same cohort of a cultured stock. Our analysis obtained a total of 151 253 SLAFs, of which 65.88% (99 652) were identified to be polymorphic, scoring a total of 710 567 putative SNPs. Further filtration resulted in a final panel of 1782 SNP loci. The data derived from this work could be beneficial for understanding the genetics of complex phenotypic traits as well as for developing marker‐selection‐assisted breeding programs in L. crocea.     

Phylogenetic estimation of Sciaenidae in the East China Sea inferred from nuclear EPIC DNA sequence variation

Here we describe a phylogenetic analysis of sciaenids of the East China Sea based on nuclear exon-primed intron-crossing genes (EPIC markers) and a mitochondrial gene (CO1). Separate analyses of the two data partitions resulted in mostly congruent trees. Although there were some differences in the classification of these species, the main difference between trees obtained by the mitochondrial gene (CO1) and nuclear DNA sequences was the position of Miichthys miiuy and Johnius belangerii. In the mitochondrial phylogeny, Johnius belangerii was placed at the most basal position forming an individual clade, while other species formed another large cluster. Miichthys miiuy formed an independent basal sub-clade grouped with Larimichthys and Collichthys. Collichthys lucidus was grouped with Larimichthys crocea and Larimichthys polyactis. Trees based on the nuclear genes differed somewhat from those based on the CO1 mitochondrial gene. In this analysis, two groups resulted, the Larimichthys and Collichthys clade, and another clade including a total of five species: Johnius belangerii, Nibea albiflora, Pennahia argentata, Sciaenops ocellatus, and Argyrosomus japonicus; Johnius belangerii clustered with Nibea albiflora. Miichthys miiuy was placed at the basal position of the other cluster because it was an independent basal sub-clade grouped with Johnius belangerii, Nibea albiflora, Pennahia argentata, Sciaenops ocellatus, and Argyrosomus japonicus. Many aspects of the phylogeny of the Sciaenidae remain unresolved, and further analysis based on more molecular information and extensive taxon sampling is necessary to elucidate the phylogenetic relationships among the major lineages within Sciaenidae.     

A high‐quality genome of taro (Colocasia esculenta (L.) Schott), one of the world's oldest crops

Taro (Colocasia esculenta (L.), Schott), from the Araceae family, is one of the oldest crops with important edible, medicinal, nutritional and economic value. Taro is a highly polymorphic species including diverse genotypes adapted to a broad range of environments, but the taro genome has rarely been investigated. Here, a high‐quality chromosome‐level genome of C. esculenta was assembled using data sequenced by Illumina, PacBio and Nanopore platforms. The assembled genome size was 2,405 Mb with a contig N50 of 400.0 kb and a scaffold N50 of 159.4 Mb. In total, 2,311 Mb (96.09%) of the contig sequences was anchored onto 14 chromosomes to form pseudomolecules, and 2,126 Mb (88.43%) was annotated as repetitive sequences. Of the 28,695 predicted protein‐coding genes, 26,215 genes (91.4%) could be functionally annotated. On the basis of phylogenetic analysis using 769 genes, C. esculenta and Spirodela polyrhiza were placed on one branch of the tree that diverged approximately 73.23 million years ago. The synteny analyses showed that there have been two whole‐genome duplication events in C. esculenta separated by a relatively short gap. According to comparative genome analysis, a larger number (1,189) of distinct gene families and long terminal repeats were enriched in C. esculenta. Our high‐quality taro genome will provide valuable resources for further genetic, ecological and evolutionary analyses of taro or other species in the Araceae.     

Identification of a large dataset of SNPs in Larimichthys polyactis using high-throughput 2b-RAD sequencing

The small yellow croaker (Larimichthys polyactis) is one of the most commercially exploited marine fishes along the coast of the Yellow-Bohai Sea and the East China Sea. In this study, we used next-generation high-throughput 2b-RAD-seq technology to identify novel SNPs in L. polyactis. We scored a total of 1 374 008 putative SNPs genome-wide. Further filtration yielded a final dataset of 6457 high-quality SNPs. These SNP markers presented sufficient power in detecting genetic distinction between two wild samples from the Yellow-Bohai Sea and one captive sample from the East China Sea, and inbreeding reflected by strong heterozygote deficiency within the samples; some of the genetic polymorphisms are diagnostic among the samples and related to biological functions. The panel of SNPs could be used as powerful tools in further population genetic and stock assessment research in L. polyactis as it has relatively scarce genomic resources. The findings from this study will advance our understanding of population and functional genomics for facilitating fishery resource management and developing desirable characteristics for the benefit of culture and farming of L. polyactis.     

Genome size and phylogenetic analysis of the A and L races of <Emphasis Type="Italic">Botryococcus braunii</Emphasis>

Botryococcus braunii (Chlorophyta, Botryococcaceae) is a colony-forming green microalga that produces large amounts of liquid hydrocarbons, which can be converted into transportation fuels. There are three different races of B. braunii, A, B, and L, that are distinguished based on the type of hydrocarbon each produces. Each race also has many strains that are distinguished by the location from which they were collected. While B. braunii has been well studied for the chemistry of the hydrocarbon production, very little is known about the molecular biology of B. braunii. To begin to address this problem, we determined the genome size of the A race, Yamanaka strain, and the L race, Songkla Nakarin strain, of B. braunii. Flow cytometry analysis indicates that the A race of B. braunii has a genome size of 166.0 ± 0.4 Mb, while the L race has a substantially larger genome size at 211.3 ± 1.7 Mb. We also used phylogenetic analysis with the nuclear small subunit (18S) rRNA gene to classify strains of the A and B races that have not yet been compared evolutionarily to previously published B. braunii phylogenetics. The analysis suggests that the evolutionary relationship between B. braunii races is correlated with the type of liquid hydrocarbon they produce.     

Impacts of ship noise on the growth and immunophysiological response in the juveniles of two Sciaenidae species,Larimichthys crocea and Nibea albiflora

The wild stocks of large yellow croaker Larimichthys crocea and yellow drum Nibea albiflora have been depleted since the 1980. Although the Chinese government has put a lot of effort (e.g., implementing fishing ban and releasing artificially bred fish into naturally distributed sea areas) into restoring their natural resources, the restoration still fails to expectation. We speculate that the slow restoration may be related to the increasing ship noise near the spawning grounds of these two species of fish. To test the harmfulness of ship noise on the L. crocea and N. albiflora, in our previous study, we studied the impacts of ship noise on their behaviors, parallelly, in the present study, we focused on the impacts of noise on their immunophysiological responses. We designed two experiments. In the first experiment, the juveniles were exposed to 120 dB ship noise for only once, then the plasma physiological indices were monitored every 1 hr within 4 hr. In the second experiment, the juveniles were exposed to 120 dB noise twice a day for 30 days, afterwards, the growth, plasma immune indices and intestinal microbiota were analyzed. The results showed that stimulated by ship noise, the physiological indices of L. crocea and N. albiflora both increased sharply within 3 hr. After a month noise stimulation, the growth and immune indices decreased significantly, and the proportion of intestinal microbiota was seriously imbalanced, showing Vibrio and Pseudomonas were dominant and other genera's abundance was quite low, especially some common intestinal probiotics. In addition, we also found that N. albiflora may be more sensitive to ship noise than L. crocea in terms of the required duration of physiological indices recovering to unstimulated level, and growth. This study highlights that ship noise negatively affects L. crocea and N. albiflora, which is helpful for some departments taking measures to protect the natural resource of these two species of fish.     

Genetic diversity in the mtDNA control region and population structure in the small yellow croaker <Emphasis Type="Italic">Larimichthys polyactis</Emphasis>

The genetic diversity and population genetic structure of the small yellow croaker (Larimichthys polyactis) were investigated. One hundred and fourteen individuals were sampled from 8 localities of the Yellow Sea and the northern East China Sea. Genetic variation in DNA sequences were examined from the first hypervariable region (HVR-1) of the mitochondrial DNA control region. High levels of haplotype diversity (h = 0.98 ± 0.87%) in the HVR-1 region were detected, indicating a high level of genetic diverstiy. A total of 84 polymorphic sites were found, and 87 haplotypes were defined. The pairwise nucleotide differences between samples ranged from 3.83 ± 2.19 to 6.56 ± 3.25. The demographic history of L. polyactis was examined by using neutrality tests and mismatch distribution analysis, which indicated a Pleistocene population expansion at about 49,300–197,000 years. The star burst structure of the minimum spanning tree also suggestted a very recent origin for most haplotypes. Hierarchical molecular variance analysis (AMOVA) and conventional population Fst comparisons revealed no significant genetic structure throughout the examined range, which is inconsistent with previous findings based on the morphological and ecological studies. Long-term dispersal and high gene flow likely have contributed to the genetically homogeneous population structure of the species. The knowledge on genetic diversity and genetic structure will be crucial to establish appropriate fishery management stocks for the species.     

Chromosome-level genome assembly of an agricultural pest,the rice leaffolder Cnaphalocrocis exigua (Crambidae,Lepidoptera)

The rice leaffolder Cnaphalocrocis exigua (Crambidae, Lepidoptera) is an important agricultural pest that damages rice crops and other members of related grass families. C. exigua exhibits a very similar morphological phenotype and feeding behaviour to C. medinalis, another species of rice leaffolder whose genome was recently reported. However, genomic information for C. exigua remains extremely limited. Here, we used a hybrid strategy combining different sequencing technologies, including Illumina, PacBio, 10× Genomics, and Hi – C scaffolding, to generate a high-quality chromosome-level genome assembly of C. exigua. We initially obtained a 798.8 Mb assembly with a contig N50 size of 2.9 Mb, and the N50 size was subsequently increased to 25.7 Mb using Hi – C technology to anchor 1413 scaffolds to 32 chromosomes. We detected a total of 97.7% Benchmarking Universal Single-Copy Orthologues (BUSCO) in the genome assembly, which was comprised of ~52% repetitive sequence and annotated 14,922 protein-coding genes. Of note, the Z and W sex chromosomes were assembled and identified. A comparative genomic analysis demonstrated that despite the high synteny observed between the two rice leaffolders, the species have distinct genomic features associated with expansion and contraction of gene families and selection pressure. In summary, our chromosome-level genome assembly and comparative genomic analysis of C. exigua provide novel insights into the evolution and ecology of this rice insect pests and offer useful information for pest control.     

A chromosome‐level genome assembly of the parasitoid wasp Pteromalus puparum

Parasitoid wasps represent a large proportion of hymenopteran species. They have complex evolutionary histories and are important biocontrol agents. To advance parasitoid research, a combination of Illumina short‐read, PacBio long‐read and Hi‐C scaffolding technologies was used to develop a high‐quality chromosome‐level genome assembly for Pteromalus puparum, which is an important pupal endoparasitoid of caterpillar pests. The chromosome‐level assembly has aided in studies of venom and detoxification genes. The assembled genome size is 338 Mb with a contig N50 of 38.7 kb and a scaffold N50 of 1.16 Mb. Hi‐C analysis assembled scaffolds onto five chromosomes and raised the scaffold N50 to 65.8 Mb, with more than 96% of assembled bases located on chromosomes. Gene annotation was assisted by RNA sequencing for the two sexes and four different life stages. Analysis detected 98% of the BUSCO (Benchmarking Universal Single‐Copy Orthologs) gene set, supporting a high‐quality assembly and annotation. In total, 40.1% (135.6 Mb) of the assembly is composed of repetitive sequences, and 14,946 protein‐coding genes were identified. Although venom genes play important roles in parasitoid biology, their spatial distribution on chromosomes was poorly understood. Mapping has revealed venom gene tandem arrays for serine proteases, pancreatic lipase‐related proteins and kynurenine–oxoglutarate transaminases, which have amplified in the P. puparum lineage after divergence from its common ancestor with Nasonia vitripennis. In addition, there is a large expansion of P450 genes in P. puparum. These examples illustrate how chromosome‐level genome assembly can provide a valuable resource for molecular, evolutionary and biocontrol studies of parasitoid wasps.     

Chromosomal‐level genomes of three rice planthoppers provide new insights into sex chromosome evolution

The brown planthopper Nilaparvata lugens, white‐backed planthopper Sogatella furcifera, and small brown planthopper Laodelphax striatellus are three major insect pests of rice. They are genetically close; however, they differ in several ecological traits such as host range, migration capacity, and in their sex chromosomes. Though the draft genome of these three planthoppers have been previously released, the quality of genome assemblies need to be improved. The absence of chromosome‐level genome resources has hindered in‐depth research of these three species. Here, we performed a de novo genome assembly for N. lugens to increase its genome assembly quality with PacBio and Illumina platforms, increasing the contig N50 to 589.46 Kb. Then, with the new N. lugens genome and previously reported S. furcifera and L. striatellus genome assemblies, we generated chromosome‐level scaffold assemblies of these three planthopper species using HiC scaffolding technique. The scaffold N50s significantly increased to 77.63 Mb, 43.36 Mb and 29.24 Mb for N. lugens, S. furcifera and L. striatellus, respectively. To identify sex chromosomes of these three planthopper species, we carried out genome re‐sequencing of males and females and successfully determined the X and Y chromosomes for N. lugens, and X chromosome for S. furcifera and L. striatellus. The gene content of the sex chromosomes showed high diversity among these three planthoppers suggesting the rapid evolution of sex‐linked genes, and all chromosomes showed high synteny. The chromosome‐level genome assemblies of three planthoppers would provide a valuable resource for a broad range of future research in molecular ecology, and subsequently benefits development of modern pest control strategies.     

Factors affecting morphological development of the sagittal otolith in juvenile and adult small yellow croaker (Larimichthys polyactis Bleeker, 1877)

The morphology and morphometrics of the sagittal otoliths of small yellow croaker (Larimichthys polyactis) from the southern Yellow Sea were investigated. Study objectives were to evaluate the shape variability and morphometric variables of sagittae of juveniles and adults as related to developmental changes and habit shift. A total of 152 fish were sampled from April to June of 2012 and 2013, along the coastal waters of the Lüsi spawning ground in the southern Yellow Sea. Changes in otolith shapes from the juvenile to the adult were presented with the rim development through the entire‐lobed‐entire transition and with the curvature of the cauda toward the ventral margin. The otolith elongation in the juvenile stage occurred at 10–20 mm standard length (SL) and was likely associated with the formation of otolith accessory growth centers from larvae to juvenile. The L. polyactis sagittal otoliths acquired their definitive shape at 130 mm SL maturity. Ontogeny on otolith shape might be related, for example, to the factors of metamorphic development, feeding habitat, and ambient water salinity, which varied throughout the growth of L. polyactis.     

The complete mitochondrial genome of bighead croaker, <Emphasis Type="Italic">Collichthys niveatus</Emphasis> (Perciformes,Sciaenidae): structure of control region and phylogenetic considerations

Sciaenidae is a diverse, commercially important family. To understand the phylogenetic position of Collichthys niveatus in this family, we present its complete mitochondrial genome sequence. The genome is 16469 bp in length and contains 37 mitochondrial genes (13 protein-coding genes, 2 ribosomal RNA genes and 22 transfer RNA genes) and a control region (CR) as in other bony fishes. Further sequencing for the complete control region was performed on Collichthys lucida. Although the conserved sequence domains such as extend termination associated sequence (ETAS) and conserved sequence block domains (CSB-1, CSB-2 and CSB-3) are recognized in the control region of the two congeneric species, the typical central conserved blocks (CSB-F, CSB-E and CSB-D) could not be detected, while they are found in Miichthys miiuy and Cynoscion acoupa of Sciaenidae and other Percoidei fishes. Phylogenetic analyses do not support the monophyly of Pseudosciaeniae, which is against with the morphological results. C. niveatus is most closely related to Larimichthys polyactis, and Collichthys and Larimichthys may be merged into one genus, based on the current datasets.     

The pomegranate (Punica granatum L.) draft genome dissects genetic divergence between soft‐ and hard‐seeded cultivars

Complete and highly accurate reference genomes and gene annotations are indispensable for basic biological research and trait improvement of woody tree species. In this study, we integrated single‐molecule sequencing and high‐throughput chromosome conformation capture techniques to produce a high‐quality and long‐range contiguity chromosome‐scale genome assembly of the soft‐seeded pomegranate cultivar ‘Tunisia’. The genome covers 320.31 Mb (scaffold N50 = 39.96 Mb; contig N50 = 4.49 Mb) and includes 33 594 protein‐coding genes. We also resequenced 26 pomegranate varieties that varied regarding seed hardness. Comparative genomic analyses revealed many genetic differences between soft‐ and hard‐seeded pomegranate varieties. A set of selective loci containing SUC8‐like, SUC6, FoxO and MAPK were identified by the selective sweep analysis between hard‐ and soft‐seeded populations. An exceptionally large selective region (26.2 Mb) was identified on chromosome 1. Our assembled pomegranate genome is more complete than other currently available genome assemblies. Our results indicate that genomic variations and selective genes may have contributed to the genetic divergence between soft‐ and hard‐seeded pomegranate varieties.