Anopheles stephensi is the key vector of malaria throughout the Indian subcontinent and Middle East and an emerging model for molecular and genetic studies of mosquito-parasite interactions. The type form of the species is responsible for the majority of urban malaria transmission across its range.

Many animal phyla have no representatives within the catalog of whole metazoan genome sequences. This dataset fills in one gap in the genome knowledge of animal phyla with a draft genome of Bugula neritina (phylum Bryozoa). Interest in this species spans ecology and biomedical sciences because B. neritina is the natural source of bioactive compounds called bryostatins. Here we present a draft assembly of the B. neritina genome obtained from PacBio and Illumina HiSeq data, as well as genes and proteins predicted de novo and verified using transcriptome data, along with the functional annotation. These sequences will permit a better understanding of host-symbiont interactions at the genomic level, and also contribute additional phylogenomic markers to evaluate Lophophorate or Lophotrochozoa phylogenetic relationships. The effort also fits well with plans to ultimately sequence all orders of the Metazoa.

Asian seabass (or commonly known as barramundi), Lates calcarifer, is a bony euryhaline teleost from the Family Latidae, inhabiting nearshore, estuarine, and marine connected freshwaters throughout the tropical Indo-West Pacific region. The species is catadromous, whereby adults spawn in salinities between 28 and 34 ppt at the mouth of estuaries, with resultant juveniles usually moving into brackish and freshwater systems to mature, before returning to the sea to spawn again as adults. The species lives in both marine and freshwater habitats and can move quickly between the two; thus, the species' ability to tolerate changes in salinity makes it a good candidate for studying the salinity acclimation response in teleosts. In this study, the transcriptome of two major osmoregulatory organs (gills and kidneys) of young juvenile Asian seabass reared in freshwater and seawater were compared. The euryhaline nature of Asian seabass was found to be highly pliable and the moldability of the trait was further confirmed by histological analyses of gills and kidneys. Differences in major expression pathways were observed, with differentially expressed genes including those related to osmoregulation, tissue/organ morphogenesis, and cell volume regulation as central to the osmo-adaptive response. Additionally, genes coding for mucins were upregulated specifically under saline conditions, whereas several genes important for growth and development, as well as circadian entrainment were specifically enriched in fish reared in freshwater. Routing of the circadian rhythm mediated by salinity changes could be the initial step in salinity acclimation and possibly migration in euryhaline fish species such as the Asian seabass.

Receptor-like cytoplasmic kinases (RLCKs) in plants belong to the super family of receptor-like kinases (RLKs). These proteins show homology to RLKs in kinase domain but lack the transmembrane domain. Some of the functionally characterized RLCKs from plants have been shown to play roles in development and stress responses. Previously, 149 and 187 RLCK encoding genes were identified from Arabidopsis and rice, respectively. By using HMM-based domain structure and phylogenetic relationships, we have identified 379 OsRLCKs from rice. OsRLCKs are distributed on all 12 chromosomes of rice and some members are located on duplicated chromosomal segments. Several OsRLCKs probably also undergo alternative splicing, some having evidence only in the form of gene models. To understand their possible functions, expression patterns during landmark stages of vegetative and reproductive development as well as abiotic and biotic stress using microarray and MPSS-based data were analyzed. Real-time PCR-based expression profiling for a selected few genes confirmed the outcome of microarray analysis. Differential expression patterns observed for majority of OsRLCKs during development and stress suggest their involvement in diverse functions in rice. Majority of the stress-responsive OsRLCKs were also found to be localized within mapped regions of abiotic stress QTLs. Outcome of this study would help in selecting organ/development stage specific OsRLCK genes/targets for functional validation studies.

The Asian seabass is an important marine food fish that has been cultured for several decades in Asia Pacific. However, the lack of a high quality reference genome has hampered efforts to improve its selective breeding. A 3D BAC pool set generated in this study was screened using 22 SSR markers located on linkage group 2 which contains a growth-related QTL region. Seventy-two clones corresponding to 22 FPC contigs were sequenced by Illumina MiSeq technology. We co-assembled the MiSeq-derived scaffolds from each FPC contig with error-corrected PacBio reads, resulting in 187 sequences covering 9.7 Mb. Eleven genes annotated within this region were found to be potentially associated with growth and their tissue-specific expression was investigated. Correlation analysis demonstrated that SNPs in ctsb, skp1 and ppp2ca can be potentially used as markers for selecting fast-growing fingerlings. Conserved syntenies between seabass LG2 and five other teleosts were identified. This study i) provided a 10 Mb targeted genome assembly; ii) demonstrated NGS of BAC pools as a potential approach for mining candidates underlying QTLs of this species; iii) detected eleven genes potentially responsible for growth in the QTL region; and iv) identified useful SNP markers for selective breeding programs of Asian seabass.

The first complete mitochondrial genome sequence of parthenogenetic Caucasian rock lizard Darevskia unisexualis (Lacertidae family) is determined by hybrid assembly with Illumina HiSeq and PacBio RS II platforms. The circular 21.4 kbp mitogenome contains 13 protein-coding genes, 12S and 16S rRNA genes, 20 tRNAs, two pseudogenized tRNAs, and one long tandem repeats with 4.1 kbp length formed by 59 bp monomer repeated x70.6 times located before control region. This finding represents a new example of mitogenome variation in lizards of hybrid origin, providing fundamental data for following study of a unique hybridization system formed by parthenogenetic and bisexual species in the mountain steppe of central Armenia.

The Asian seabass (Lates calcarifer) is a bony fish from the Latidae family, which is widely distributed in the tropical Indo-West Pacific region. The karyotype of the Asian seabass contains 24 pairs of A chromosomes and a variable number of AT- and GC-rich B chromosomes (Bchrs or Bs). Dot-like shaped and nucleolus-associated AT-rich Bs were microdissected and sequenced earlier. Here we analyzed DNA fragments from Bs to determine their repeat and gene contents using the Asian seabass genome as a reference. Fragments of 75 genes, including an 18S rRNA gene, were found in the Bs; repeats represented 2% of the Bchr assembly. The 18S rDNA of the standard genome and Bs were similar and enriched with fragments of transposable elements. A higher nuclei DNA content in the male gonad and somatic tissue, compared to the female gonad, was demonstrated by flow cytometry. This variation in DNA content could be associated with the intra-individual variation in the number of Bs. A comparison between the copy number variation among the B-related fragments from whole genome resequencing data of Asian seabass individuals identified similar profiles between those from the South-East Asian/Philippines and Indian region but not the Australian ones. Our results suggest that Bs might cause variations in the genome among the individuals and populations of Asian seabass. A personalized copy number approach for segmental duplication detection offers a suitable tool for population-level analysis across specimens with low coverage genome sequencing.

Domestic cats enjoy an extensive veterinary medical surveillance which has described nearly 250 genetic diseases analogous to human disorders. Feline infectious agents offer powerful natural models of deadly human diseases, which include feline immunodeficiency virus, feline sarcoma virus and feline leukemia virus. A rich veterinary literature of feline disease pathogenesis and the demonstration of a highly conserved ancestral mammal genome organization make the cat genome annotation a highly informative resource that facilitates multifaceted research endeavors.

As part of our Asian seabass genome project, we are generating an inventory of repeat elements in the genome and transcriptome. The karyotype showed a diploid number of 2n = 24 chromosomes with a variable number of B-chromosomes. The transcriptome and genome of Asian seabass were searched for repetitive elements with experimental and bioinformatics tools. Six different types of repeats constituting 8-14% of the genome were characterized. Repetitive elements were clustered in the pericentromeric heterochromatin of all chromosomes, but some of them were preferentially accumulated in pretelomeric and pericentromeric regions of several chromosomes pairs and have chromosomes specific arrangement. From the dispersed class of fish-specific non-LTR retrotransposon elements Rex1 and MAUI-like repeats were analyzed. They were wide-spread both in the genome and transcriptome, accumulated on the pericentromeric and peritelomeric areas of all chromosomes. Every analyzed repeat was represented in the Asian seabass transcriptome, some showed differential expression between the gonads. The other group of repeats analyzed belongs to the rRNA multigene family. FISH signal for 5S rDNA was located on a single pair of chromosomes, whereas that for 18S rDNA was found on two pairs. A BAC-derived contig containing rDNA was sequenced and assembled into a scaffold containing incomplete fragments of 18S rDNA. Their assembly and chromosomal position revealed that this part of Asian seabass genome is extremely rich in repeats containing evolutionarily conserved and novel sequences. In summary, transcriptome assemblies and cDNA data are suitable for the identification of repetitive DNA from unknown genomes and for comparative investigation of conserved elements between teleosts and other vertebrates.

We report here the ~670 Mb genome assembly of the Asian seabass (Lates calcarifer), a tropical marine teleost. We used long-read sequencing augmented by transcriptomics, optical and genetic mapping along with shared synteny from closely related fish species to derive a chromosome-level assembly with a contig N50 size over 1 Mb and scaffold N50 size over 25 Mb that span ~90% of the genome. The population structure of L. calcarifer species complex was analyzed by re-sequencing 61 individuals representing various regions across the species' native range. SNP analyses identified high levels of genetic diversity and confirmed earlier indications of a population stratification comprising three clades with signs of admixture apparent in the South-East Asian population. The quality of the Asian seabass genome assembly far exceeds that of any other fish species, and will serve as a new standard for fish genomics.

Patterns of genetic and genomic variance are informative in inferring population history for human, model species and endangered populations.

As nutrition-related expenses constitute the majority of the costs for aquaculture farms, it is essential for them to use feeds that provide an ideal combination of nutrients for the species of choice. In this study, the relative effect of consuming three different pelleted feeds (B, C and D) in comparison to frozen baitfish (A; control) were compared on juvenile Asian seabass (77.3 ± 22.4g) that were selected for increased growth rate over two generations. Our objectives were: 1) to evaluate the effects of different pelleted feeds based on overall physiological changes and nutritional quality of fillets; 2) improve our understanding of the underlying mechanisms with transcriptomic analysis; 3) if possible, identify the feed type that supports the growth of these fishes without substantially reducing the nutritional quality of fillet. The growth performance, fatty acid composition of fillet, hepatic histology and transcriptome of the fishes (Groups A-D) were analyzed. The majority of fatty acids of the fillets, except γ-linolenic acid (GLA, C18:3n6), correlated significantly with the respective diets. Asian seabass fed Feed C showed highest specific growth rate (SGR) and feed conversion efficiency (FCE) with closest histology and transcriptomic profile to control, but their fillet contained the highest n6/n3 ratio. When the liver-based transcriptomes were analyzed, a complex set of differentially expressed genes were detected between groups fed pelleted feeds and controls as well as among the pellet-fed groups themselves. Significant enrichment of genes with growth-related function tallied with the morphological data measured. When compared with control (Group A), 'Biosynthesis of unsaturated fatty acids' and 'Steroid biosynthesis' pathways were significantly enriched in pellet-fed groups. Reduced goblet cell numbers were observed in the gut of pellet-fed fish compared to controls and fads6 was found to be a suitable candidate gene to separate wild-caught Asian seabass, from pellet-fed ones. These results provide insights for researchers on the various effects of feeds on the biochemistry and global gene expression of the fish and potentially for seabass farms to make more informed feed choices.

Solenodons are insectivores that live in Hispaniola and Cuba. They form an isolated branch in the tree of placental mammals that are highly divergent from other eulipothyplan insectivores The history, unique biology, and adaptations of these enigmatic venomous species could be illuminated by the availability of genome data. However, a whole genome assembly for solenodons has not been previously performed, partially due to the difficulty in obtaining samples from the field. Island isolation and reduced numbers have likely resulted in high homozygosity within the Hispaniolan solenodon (Solenodon paradoxus). Thus, we tested the performance of several assembly strategies on the genome of this genetically impoverished species. The string graph-based assembly strategy seemed a better choice compared to the conventional de Bruijn graph approach due to the high levels of homozygosity, which is often a hallmark of endemic or endangered species. A consensus reference genome was assembled from sequences of 5 individuals from the southern subspecies (S. p. woodi). In addition, we obtained an additional sequence from 1 sample of the northern subspecies (S. p. paradoxus). The resulting genome assemblies were compared to each other and annotated for genes, with an emphasis on venom genes, repeats, variable microsatellite loci, and other genomic variants. Phylogenetic positioning and selection signatures were inferred based on 4,416 single-copy orthologs from 10 other mammals. We estimated that solenodons diverged from other extant mammals 73.6 million years ago. Patterns of single-nucleotide polymorphism variation allowed us to infer population demography, which supported a subspecies split within the Hispaniolan solenodon at least 300 thousand years ago.

Genome-wide assessment of genetic diversity has the potential to increase the ability to understand admixture, inbreeding, kinship and erosion of genetic diversity affecting both captive (ex situ) and wild (in situ) populations of threatened species. The sable antelope (Hippotragus niger), native to the savannah woodlands of sub-Saharan Africa, is a species that is being managed ex situ in both public (zoo) and private (ranch) collections in the United States. Our objective was to develop whole genome sequence resources that will serve as a foundation for characterizing the genetic status of ex situ populations of sable antelope relative to populations in the wild. Here we report the draft genome assembly of a male sable antelope, a member of the subfamily Hippotraginae (Bovidae, Cetartiodactyla, Mammalia). The 2.596 Gb draft genome consists of 136,528 contigs with an N50 of 45.5 Kbp and 16,927 scaffolds with an N50 of 4.59 Mbp. De novo annotation identified 18,828 protein-coding genes and repetitive sequences encompassing 46.97% of the genome. The discovery of single nucleotide variants (SNVs) was assisted by the re-sequencing of seven additional captive and wild individuals, representing two different subspecies, leading to the identification of 1,987,710 bi-allelic SNVs. Assembly of the mitochondrial genomes revealed that each individual was defined by a unique haplotype and these data were used to infer the mitochondrial gene tree relative to other hippotragine species. The sable antelope genome constitutes a valuable resource for assessing genome-wide diversity and evolutionary potential, thereby facilitating long-term conservation of this charismatic species.

Vaccination is an important strategy in the protection of aquaculture species from major diseases. However, we still do not have a good understanding of the mechanisms underlying vaccine-induced disease resistance. This is further complicated by the presence of several lymphoid organs that play different roles when mounting an immune response. In this study, we attempt to elucidate some of these mechanisms using a microarray-based approach. Asian seabass (Lates calcarifer) were vaccinated against Streptococcus iniae and the transcriptomic changes within the spleen and head kidney at one and seven days post-vaccination were profiled. We subsequently challenged the seabass at three weeks post-vaccination with live S. iniae and similarly profiled the transcriptomes of the two organs after the challenge. We found that vaccination induced an early, but transient transcriptomic change in the spleens and a delayed response in the head kidneys, which became more similar to one another compared to un-vaccinated ones. When challenged with the pathogen, the spleen, but not the head kidneys, responded transcriptomically at 25-29 hours post-challenge. A unique set of genes, in particular those involved in the activation of NF-κB signaling, was up-regulated in the vaccinated spleens upon pathogen challenge but not in the un-vaccinated spleens. A semi-quantitative PCR detection of S. iniae using metagenomic DNA extracted from the water containing the seabass also revealed that vaccination resulted in reduction of pathogen shedding. This result indicated that vaccination not only led to a successful immune defense against the infection, but also reduced the chances for horizontal transmission of the pathogen. In conclusion, we have provided a transcriptomic analysis of how the teleost spleen and head kidneys responded to vaccination and subsequent infection. The different responses from the two organs are suggestive of their unique roles in establishing a vaccine-induced disease resistance.

As complete absence of germ cells leads to sterile males in zebrafish, we explored the relationship between primordial germ cell (PGC) number and sexual development. Our results revealed dimorphic proliferation of PGCs in the early zebrafish larvae, marking the beginning of sexual differentiation. We applied morpholino-based gene knockdown and cell transplantation strategies to demonstrate that a threshold number of PGCs is required for the stability of ovarian fate. Using histology and transcriptomic analyses, we determined that zebrafish gonads are in a meiotic ovarian stage at 14 days postfertilization and identified signaling pathways supporting meiotic oocyte differentiation and eventual female fate. The development of PGC-depleted gonads appears to be restrained and delayed, suggesting that PGC number may directly regulate the variability and length of gonadal transformation and testicular differentiation in zebrafish. We propose that gonadal transformation may function as a developmental buffering mechanism to ensure the reproductive outcome.

It is generally believed that the last eukaryotic common ancestor (LECA) was a unicellular organism with motile cilia. In the vertebrates, the winged-helix transcription factor FoxJ1 functions as the master regulator of motile cilia biogenesis. Despite the antiquity of cilia, their highly conserved structure, and their mechanism of motility, the evolution of the transcriptional program controlling ciliogenesis has remained incompletely understood. In particular, it is presently not known how the generation of motile cilia is programmed outside of the vertebrates, and whether and to what extent the FoxJ1-dependent regulation is conserved. We have performed a survey of numerous eukaryotic genomes and discovered that genes homologous to foxJ1 are restricted only to organisms belonging to the unikont lineage. Using a mis-expression assay, we then obtained evidence of a conserved ability of FoxJ1 proteins from a number of diverse phyletic groups to activate the expression of a host of motile ciliary genes in zebrafish embryos. Conversely, we found that inactivation of a foxJ1 gene in Schmidtea mediterranea, a platyhelminth (flatworm) that utilizes motile cilia for locomotion, led to a profound disruption in the differentiation of motile cilia. Together, all of these findings provide the first evolutionary perspective into the transcriptional control of motile ciliogenesis and allow us to propose a conserved FoxJ1-regulated mechanism for motile cilia biogenesis back to the origin of the metazoans.

A reduction in daily caloric or nutrient intake has been observed to promote health benefits in mammals and other vertebrates. Feed Restriction (FR), whereby the overall food intake of the organism is reduced, has been explored as a method to improve metabolic and immune health, as well as to optimize productivity in farming. However, less is known regarding the molecular and physiological consequences of FR. Using the model organism, Danio rerio, we investigated the impact of a short-term (month-long) FR on growth, gut morphology and gene expression. Our data suggest that FR has minimal effects on the average growth rates, but it may affect weight and size heterogeneity in a sex-dependent manner. In the gut, we observed a significant reduction in gut circumference and generally lower mucosal heights, whereas other parameters remained unchanged. Gene Ontology (GO), EuKaryotic Orthologous Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified numerous metabolic, reproductive, and immune response pathways that were affected by FR. These results broaden our understanding of FR and contribute towards growing knowledge of its effects on vertebrate health.

The Puma lineage within the family Felidae consists of 3 species that last shared a common ancestor around 4.9 million years ago. Whole-genome sequences of 2 species from the lineage were previously reported: the cheetah (Acinonyx jubatus) and the mountain lion (Puma concolor). The present report describes a whole-genome assembly of the remaining species, the jaguarundi (Puma yagouaroundi). We sequenced the genome of a male jaguarundi with 10X Genomics linked reads and assembled the whole-genome sequence. The assembled genome contains a series of scaffolds that reach the length of chromosome arms and is similar in scaffold contiguity to the genome assemblies of cheetah and puma, with a contig N50 = 100.2 kbp and a scaffold N50 = 49.27 Mbp. We assessed the assembled sequence of the jaguarundi genome using BUSCO, aligned reads of the sequenced individual and another published female jaguarundi to the assembled genome, annotated protein-coding genes, repeats, genomic variants and their effects with respect to the protein-coding genes, and analyzed differences of the 2 jaguarundis from the reference mitochondrial genome. The jaguarundi genome assembly and its annotation were compared in quality, variants, and features to the previously reported genome assemblies of puma and cheetah. Computational analyzes used in the study were implemented in transparent and reproducible way to allow their further reuse and modification.

The mitochondrial genome of the long-spined black sea urchin, Diadema antillarum, was sequenced using Illumina next-generation sequencing technology. The complete mitogenome is 15,708 bp in length, containing two rRNA, 22 tRNA and 13 protein-coding genes, plus a noncoding control region of 133 bp. The nucleotide composition is 18.37% G, 23.79% C, 26.84% A and 30.99% T. The A + T bias is 57.84%. Phylogenetic analysis based on 12 complete mitochondrial genomes of sea urchins, including four species of the family Diadematidae, supported familial monophyly; however, the two Diadema species, D. antillarum and D. setosum were not recovered as sister taxa.