De Novo Transcriptome Assembly of the Chinese Swamp Buffalo by RNA Sequencing and SSR Marker Discovery

SSR Marker Acquisition and Application from Transcriptome of Captive Chinese Forest Musk Deer (Moschus berezovskii)

Forest musk deer (Moschus berezovskii) is one of the most endangered medicinally important wild animals in the world. Forest musk deer farming is the main way of production of musk. However, the single provenance and lack of genetic information lead to reduced genetic diversity of forest musk deer. Therefore, more SSR markers need to be developed to identify forest musk deer germplasm. In this study, bone marrow derived mesenchymal cells were used to construct cDNA library for transcriptome sequencing. The datasets were de novo assembled and annotated. 9 polymorphic simple sequence repeat (SSR) markers were finally identified and used to detect population genetic diversity. 6.07 Gb clean data were generated using Illumina sequencing technology, and de novo assembled into 138,591 transcripts and 81,553 unigenes. 5,777 simple sequence repeats (SSRs) were identified, in which there were 578 repeating motif types, with mono-nucleotide and tri-nucleotides comprising 55.88% and 25.60%, respectively. 100 primer pairs were designed to validate amplification and polymorphism using DNA from fecal samples. 9 polymorphic SSRs were developed and used to detect population genetic diversity of 122 forest musk deer in 2 farms. The average number of alleles per locus varied from 4 to 15 (average = 8.3). The observed heterozygosity (HO) per locus ranged from 0.102 to 0.941, while the expected heterozygosity (HE) per locus was from 0.111 to 0.651. All loci deviated significantly from the Hardy-Weinberg equilibrium (p < 0.001). The polymorphism information content (PIC) of these loci varied from 0.108 to 0.619. 9 polymorphic SSR markers were developed in this research. These sites can be used for breeding planning and conservation of germplasm resources.

Comparative transcriptome analysis revealed genes involved in the sexual size dimorphisms and expressed sequence tag-Simple Sequence Repeat loci validation in Odorrana graminea

Sexual size dimorphism (SSD) is widespread among animals and is characterized by differences in body size between sexes. Previous studies suggested SSD might reflect the adaptations of particular sexes to their specific reproductive or ecological roles. The large green cascade frogs (Odorrana graminea) exhibit obvious SSD that females are nearly twice the body size of males. However, the molecular mechanisms underlying SSD of O. graminea are still unknown. In the present study, we first obtained nearly 5 Gb of the transcriptome data through Illumina sequencing, and the de novo transcriptome assembly produced 189,868 unigenes of O. graminea. A total of 774 significantly sex-differentially expressed genes (DEGs) were identified. Of which, 436 DEGs showed significantly higher expression levels in females than those in males, whereas 338 DEGs showed significantly lower expression in females than those in males. We also found 10 sex-differentially expressed genes related to energy metabolism between sexes of O. graminea, and these DEGs were related to the estrogen signaling pathway, oxidative phosphorylation, fatty acid biosynthesis, gastric acid secretion, and nitrogen metabolism. We found that the differences in energy metabolism and steroid hormone synthesis might be the main driving force leading to the sexual growth dimorphism of O. graminea. In addition, a total of 63,269 potential EST-SSR loci and 4,669 EST-SSR loci were detected and validated in different populations of O. graminea and other species within Odorrana. The assembled transcriptome will facilitate functional genomic studies of O. graminea and the developed EST-SSR markers will contribute to the population genetics of the species within Odorrana. The sex-differentially expressed genes involved in energy metabolism might provide insights into the genetic mechanisms underlying the SSD of O. graminea.

Genes Associated with Biological Nitrogen Fixation Efficiency Identified Using RNA Sequencing in Red Clover (Trifolium pratense L.)

Commonly studied in the context of legume-rhizobia symbiosis, biological nitrogen fixation (BNF) is a key component of the nitrogen cycle in nature. Despite its potential in plant breeding and many years of research, information is still lacking as to the regulation of hundreds of genes connected with plant-bacteria interaction, nodulation, and nitrogen fixation. Here, we compared root nodule transcriptomes of red clover (Trifolium pratense L.) genotypes with contrasting nitrogen fixation efficiency, and we found 491 differentially expressed genes (DEGs) between plants with high and low BNF efficiency. The annotation of genes expressed in nodules revealed more than 800 genes not yet experimentally confirmed. Among genes mediating nodule development, four nod-ule-specific cysteine-rich (NCR) peptides were confirmed in the nodule transcriptome. Gene duplication analyses revealed that genes originating from tandem and dispersed duplication are significantly over-represented among DEGs. Weighted correlation network analysis (WGCNA) organized expression profiles of the transcripts into 16 modules linked to the analyzed traits, such as nitrogen fixation efficiency or sample-specific modules. Overall, the results obtained broaden our knowledge about transcriptomic landscapes of red clover's root nodules and shift the phenotypic description of BNF efficiency on the level of gene expression in situ.

De Novo Transcriptome Dataset Generation of the Swamp Buffalo Brain and Non-Brain Tissues

The sequenced data availability opened new horizons related to buffalo genetic control of economic traits and genomic diversity. The visceral organs (brain, liver, etc.) significantly involved in energy metabolism, docility, or social interactions. We performed swamp buffalo transcriptomic profiling of 24 different tissues (brain and non-brain) to identify novel transcripts and analyzed the differentially expressed genes (DEGs) of brain vs. non-brain tissues with their functional annotation. We obtained 178.57 Gb clean transcriptomic data with GC contents 52.77%, reference genome alignment 95.36%, exonic coverage 88.49%. Totally, 26363 mRNAs transcripts including 5574 novel genes were obtained. Further, 7194 transcripts were detected as DEGs by comparing brain vs. non-brain tissues group, of which 3,999 were upregulated and 3,195 downregulated. These DEGs were functionally associated with cellular metabolic activities, signal transduction, cytoprotection, and structural and binding activities. The related functional pathways included cancer pathway, PI3k-Akt signaling, axon guidance, JAK-STAT signaling, basic cellular metabolism, thermogenesis, and oxidative phosphorylation. Our study provides an in-depth understanding of swamp buffalo transcriptomic data including DEGs potentially involved in basic cellular activities and development that helped to maintain their working capacity and social interaction with humans, and also, helpful to disclose the genetic architecture of different phenotypic traits and their gene expression regulation.

Draft genome of the lowland anoa (Bubalus depressicornis) and comparison with buffalo genome assemblies (Bovidae, Bubalina)

Genomic data for wild species of the genus Bubalus (Asian buffaloes) are still lacking while several whole genomes are currently available for domestic water buffaloes. To address this, we sequenced the genome of a wild endangered dwarf buffalo, the lowland anoa (Bubalus depressicornis), produced a draft genome assembly and made comparison to published buffalo genomes. The lowland anoa genome assembly was 2.56 Gbp long and contained 103,135 contigs, the longest contig being 337.39 kbp long. N50 and L50 values were 38.73 and 19.83 kbp, respectively, mean coverage was 44× and GC content was 41.74%. Two strategies were adopted to evaluate genome completeness: (1) determination of genomic features with de novo and homology-based predictions using annotations of chromosome-level genome assembly of the river buffalo and (2) employment of benchmarking against universal single-copy orthologs (BUSCO). Homology-based predictions identified 94.51% complete and 3.65% partial genomic features. De novo gene predictions identified 32,393 genes, representing 97.14% of the reference’s annotated genes, whilst BUSCO search against the mammalian orthologs database identified 71.1% complete, 11.7% fragmented, and 17.2% missing orthologs, indicating a good level of completeness for downstream analyses. Repeat analyses indicated that the lowland anoa genome contains 42.12% of repetitive regions. The genome assembly of the lowland anoa is expected to contribute to comparative genome analyses among bovid species.

Novel Insight into the Potential Role of Acylglycerophosphate Acyltransferases Family Members on Triacylglycerols Synthesis in Buffalo

Acylglycerophosphate acyltransferases (AGPATs) are the rate-limiting enzymes for the de novo pathway of triacylglycerols (TAG) synthesis. Although AGPATs have been extensively explored by evolution, expression and functional studies, little is known on functional characterization of how many members of the AGPAT family are involved in TAG synthesis and their impact on the cell proliferation and apoptosis. Here, 13 AGPAT genes in buffalo were identified, of which 12 AGPAT gene pairs were orthologous between buffalo and cattle. Comparative transcriptomic analysis and real-time quantitative reverse transcription PCR (qRT-PCR) further showed that both AGPAT1 and AGPAT6 were highly expressed in milk samples of buffalo and cattle during lactation. Knockdown of AGPAT1 or AGPAT6 significantly decreased the TAG content of buffalo mammary epithelial cells (BuMECs) and bovine mammary epithelial cells (BoMECs) by regulating lipogenic gene expression (p < 0.05). Knockdown of AGPAT1 or AGPAT6 inhibited proliferation and apoptosis of BuMECs through the expression of marker genes associated with the proliferation and apoptosis (p < 0.05). Our data confirmed that both AGPAT1 and AGPAT6 could regulate TAG synthesis and growth of mammary epithelial cells in buffalo. These findings will have important implications for understanding the role of the AGPAT gene in buffalo milk performance.

Development of 16 novel EST-SSR markers for species identification and cross-genus amplification in sambar, sika, and red deer

Deer genera around the globe are threatened by anthropogenic interference. The translocation of alien species and their subsequent genetic introgression into indigenous deer populations is particularly harmful to the species of greatest conservation concern. Products derived from deer, including venison and antler velvet, are also at risk of fraudulent labeling. The current molecular markers used to genetically identify deer species were developed from genome sequences and have limited applicability for cross-species amplification. The absence of efficacious diagnostic techniques for identifying deer species has hampered conservation and wildlife crime investigation efforts. Expressed sequence tag-simple sequence repeat (EST-SSR) markers are reliable tools for individual and species identification, especially in terms of cross-species genotyping. We conducted transcriptome sequencing of sambar (Rusa unicolor) antler velvet and acquired 11,190 EST-SSRs from 65,074 newly assembled unigenes. We identified a total of 55 unambiguous amplicons in sambar (n = 45), which were selected as markers to evaluate cross-species genotyping in sika deer (Cervus nippon, n = 30) and red deer (Cervus elaphus, n = 46), resulting in cross-species amplification rates of 94.5% and 89.1%, respectively. Based on polymorphic information content (>0.25) and genotyping fidelity, we selected 16 of these EST-SSRs for species identification. This marker set revealed significant genetic differentiation based on the fixation index and genetic distance values. Principal coordinate analysis and STRUCTURE analysis revealed distinct clusters of species and clearly identified red-sika hybrids. These markers showed applicability across different genera and proved suitable for identification and phylogenetic analyses across deer species.

pSATdb: a database of mitochondrial common, polymorphic, and unique microsatellites

The polymorphic microSATellites database (pSATdb) provides information on common, polymorphic, and unique mitochondrial microsatellites.

Microsatellites, also termed as simple sequence repeats, are repetitive tracts in a DNA sequence, typically consisting of one to six nucleotides. These repeats are found in all genomes and play key roles in phylogeny and species identification. Microsatellites are highly polymorphic, and their length may differ from species to species. There are several online resources dedicated to mitochondria; however, comprehensive information is not available about the length variation of mitochondrial microsatellites. Therefore, to explore it between species among a genus, we have developed a database named pSATdb (polymorphic microSATellites database; https://lms.snu.edu.in/pSATdb/). pSATdb contains 28,710 perfect microsatellites identified across 5,976 mitochondrial genome (mt-genome) sequences from 1,576 genera which includes 1,535 (5,846 mt-genome) and 41 (130 mt-genome) genera of Metazoa and Viridiplantae, respectively. pSATdb is the only database which provides genus-wise information about the length variation of mitochondrial microsatellites. Because of the emerging role of microsatellites in genomics studies, the identified common, polymorphic, and unique microsatellites stored in pSATdb will be effectively useful in various studies including genetic diversity, mapping, marker-assisted selection, and comparative population studies.

Genomic Identification, Evolution, and Expression Analysis of Bromodomain Genes Family in Buffalo

Bromodomain (BRD) is an evolutionarily conserved protein-protein interaction module that is critical in gene regulation, cellular homeostasis, and epigenetics. This study aimed to conduct an identification, evolution, and expression analysis of the BRD gene family in the swamp buffalo (Bubalus bubalis). A total of 101 BRD protein sequences deduced from 22 BRD genes were found in the buffalo genome. The BRD proteins were classified into six groups based on phylogenetic relationships, conserved motifs, and conserved domains. The BRD genes were irregularly distributed in 13 chromosomes. Collinearity analysis revealed 20 BRD gene pairs that had remarkable homologous relationships between the buffalo and cattle, although no tandem or segmental duplication event was found in the buffalo BRD genes. Comparative transcriptomics using a 10x sequencing platform analysis showed that 22 BRD genes were identified in the Sertoli cells (SCs) at different developmental stages of buffalo. Further, the mRNA expression levels of bromodomain and the extraterminal (BET) family in SCs at the pubertal stage were higher than that at the prepubertal stage of buffalo. However, the SMARCA2, PHIP, BRD9, and TAF1 genes exhibited the opposite trend. The maturation process of SCs may be regulated by the BRD family members expressed differentially in SCs at different developmental stages of buffalo. In summary, our findings provide an understanding of the evolutionary, structural, and functional properties of the buffalo BRD family members, and further characterize the function of the BRD family in the maturation of SCs. It also provides a theoretical basis for further understanding in the future of the mechanism of SCs regulating spermatogenesis.

De Novo Transcriptome Assembly, Functional Annotation and SSR Marker Discovery of Qinling Takin (Budorcas taxicolor bedfordi)

The takin (Budorcas taxicolor) is an endemic ruminant species belonging to the bovine family. The International Union for Conservation of Nature (IUCN) has listed it as an endangered and vulnerable species. However, little is known about its molecular characterization since it lacks a reference genome. This study used RNA sequencing followed by de novo assembly, annotation and simple sequence repeats (SSRs) prediction to assess the transcriptome of Qinling takin (Budorcas taxicolor bedfordi) muscles. In total, 21,648 unigenes with an N50 and mean length of 1388 bp and 817 bp, respectively, were successfully detected and annotated against the public databases (NR, GO, KEGG, and EggNOG). Furthermore, 6222 SSRs were identified using the MIcroSAtellite (MISA) identification tool software. Taken together, these findings will provide valuable information for genetic, genomic, and evolutionary studies on takin.

Molecular and functional characterization of buffalo nasal epithelial odorant binding proteins and their structural insights by in silico and biochemical approaches

The olfactory system is capable of detecting and distinguishing thousands of environmental odorants that play a key role in reproduction, social behaviours including pheromones influenced classical events. Membrane secretary odorant binding proteins (OBPs) are soluble lipocalins, localized in the nasal membrane of mammals. They bind and carry odorants within the nasal epithelium to putative olfactory transmembrane receptors (ORs). OBP has not yet been exploited to develop a suitable technique to detect oestrus which is being reported as a difficult task in buffalo. In the present study, using molecular biology and protein engineering approaches, we have cloned six novel OBP isoforms from buffalo nasal epithelium odorant-binding proteins (bnOBPs). Furthermore, 3 D models were developed and molecular-docking, dynamics experiments were performed by in silico approaches. In particular, we found four residues (Phe104, Phe134, Phe69 and Asn118) in OBP1a, which contributed to favourable interactions towards two sex pheromones, specifically oleic acid and p-cresol. We expressed this protein in Escherichia coli from female buffalo urine and validated through fluorescence quenching studies to show similar strong binding affinities of OBP1a to oleic acid and p-cresol. By using structural data, the binding specificity was also verified by site-directed mutagenesis of the four residues followed by in vitro binding assays. Our results enable us to better understand the functions of different nasal epithelium OBP isoforms in buffaloes. They also lead to improved understanding of the interaction between olfactory proteins and odorants to develop highly selective biosensing devices for non-invasive detection of oestrus in buffaloes. Communicated by Ramaswamy H. Sarma.

Asian water buffalo: domestication, history and genetics

The domestic Asian water buffalo (Bubalus bubalis) is found on all five continents, with a global population of some 202 million. The livelihoods of more people depend on this species than on any other domestic animal. The two distinct types (river and swamp) descended from different wild Asian water buffalo (Bubalus arnee) populations that diverged some 900 kyr BP and then evolved in separate geographical regions. After domestication in the western region of the Indian subcontinent (ca. 6300 years BP), the river buffalo spread west as far as Egypt, the Balkans and Italy. Conversely, after domestication in the China/Indochina border region ca. 3000-7000 years BP, swamp buffaloes dispersed through south-east Asia and China as far as the Yangtze River valley. Molecular and morphological evidence indicates that swamp buffalo populations have strong geographic genetic differentiation and a lack of gene flow, but strong phenotypic uniformity. In contrast, river buffalo populations show a weaker phylogeographic structure, but higher phenotypic diversity (i.e. many breeds). The recent availability of a high-quality reference genome and of a medium-density marker panel for genotyping has triggered a number of genome-wide investigations on diversity, evolutionary history, production traits and functional elements. The growing molecular knowledge combined with breeding programmes should pave the way to improvements in production, environmental adaptation and disease resistance in water buffalo populations worldwide.

Genomic diversity and selection sweeps identified in Indian swamp buffaloes reveals it's uniqueness with riverine buffaloes

The present investigation was focused to study genomic diversity of Indian swamp buffalo populations through reduced representation approach (ddRAD). The heterozygosity (F_ST) among the swamp buffaloes was 0.11 between Assam and Manipuri; 0.20 between swamp (Manipuri) and riverine buffaloes; 0.30 between swamp (Manipuri) and cattle. The average observed and expected heterozygosity in swamp buffalo populations was 0.254 and 0.221 respectively. The Inbreeding coefficient (F_IS) value was 0.02 among the swamp buffaloes. PCA and structure analysis revealed Manipuri swamp buffalo was genetically distinct and closely related to Nagaland swamp buffalo and least to Assam swamp buffalo. Identification of selective sweeps revealed 1087 regions to have undergone selection related to immune response, adaptation and nervous system. A total of 3451 SSRs were identified in the genome of swamp buffaloes. The study evidenced the genomic diversity in the swamp buffalo populations and its uniqueness in comparison with riverine buffalo and cattle.

A Gene Expression Atlas of the Domestic Water Buffalo (Bubalus bubalis)

The domestic water buffalo (Bubalus bubalis) makes a major contribution to the global agricultural economy in the form of milk, meat, hides, and draught power. The global water buffalo population is predominantly found in Asia, and per head of population more people depend upon the buffalo than on any other livestock species. Despite its agricultural importance, there are comparatively fewer genomic and transcriptomic resources available for buffalo than for other livestock species. We have generated a large-scale gene expression atlas covering multiple tissue and cell types from all major organ systems collected from three breeds of riverine water buffalo (Mediterranean, Pandharpuri and Bhadawari) and used the network analysis tool Graphia Professional to identify clusters of genes with similar expression profiles. Alongside similar data, we and others have generated for ruminants as part of the Functional Annotation of Animal Genomes Consortium; this comprehensive transcriptome supports functional annotation and comparative analysis of the water buffalo genome.

Discovery of novel genic-SSR markers from transcriptome dataset of an important non-human primate, Macaca fascicularis

Macaca fascicularis, also known as the cynomolgus macaque, is an important non-human primate animal model used in biomedical research. It is an Old-World primate widely distributed in Southeast Asia and is one of the most abundant macaque species in Malaysia. However, the genetic structure of wild cynomolgus macaque populations in Malaysia has not been thoroughly elucidated. In this study, we developed genic-simple sequence repeat (genic-SSR) markers from an in-house transcriptome dataset generated from the Malaysian cynomolgus macaque via RNA sequencing, and applied these markers on 26 cynomolgus macaque individuals. A collection of 14,751 genic-SSRs were identified, where 13,709 were perfect SSRs. Dinucleotide repeats were the most common repeat motifs with a frequency of 65.05%, followed by trinucleotide repeats (20.55%). Subsequently, we designed 300 pairs of primers based on perfect di- and trinucleotide SSRs, in which 105 SSRs were associated with functional genes. A subset of 30 SSR markers were randomly selected and validated, yielding 19 polymorphic markers with an average polymorphism information content value of 0.431. The development of genic-SSR markers in this study is indeed timely to provide useful markers for functional and population genetic studies of the cynomolgus macaque and other related non-human primate species.

Buffalo milk transcriptome: A comparative analysis of early, mid and late lactation

The expression of genes and their regulation during lactation in buffaloes remains less understood. To understand the interplay of various genes and pathways, the milk transcriptome from three lactation stages of Murrah buffalo was analyzed by RNA sequencing. The filtered reads were mapped to the Bubalus bubalis as well as Bos taurus reference assemblies. The average mapping rate to water buffalo and Btau 4.6 reference sequence, was 75.5% and 75.7% respectively. Highly expressed genes (RPKM > 3000), throughout lactation included CSN2, CSN1S1, CSN3, LALBA, SPP1 and TPT1. A total of 12833 transcripts were common across all the stages, while 271, 205 and 418 were unique to early, mid and late lactation respectively. Majority of the genes throughout lactation were linked to biological functions like protein metabolism, transport and immune response. A discernible shift from metabolism in early stage to metabolism and immune response in mid stage, and an increase in immune response functions in late lactation was observed. The results provide information of candidate genes and pathways involved in the different stages of lactation in buffalo. The study also identified 14 differentially expressed and highly connected genes across the three lactation stages, which can be used as candidates for future research.

Comparative analysis of the transcriptome of the Amazonian fish species Colossoma macropomum (tambaqui) and hybrid tambacu by next generation sequencing

Background

The C. macropomum is a characiform fish from the Amazon basin that has been hybridized with other pacu species to produce commercial hybrids, such as the tambacu. However, little is known of the functional genomics of the parental species or these hybrid forms. The transcriptome of C. macropomum and tambacu were sequenced using 454 Roche platform (pyrosequencing) techniques to characterize the domains of Gene Ontology (GO) and to evaluate the levels of gene expression in the two organisms.

Results

The 8,188,945 reads were assembled into 400,845 contigs. A total of 58,322 contigs were annotated with a predominance of biological processes for both organisms, as determined by Gene Ontology (GO). Similar numbers of metabolic pathways were identified in both the C. macropomum and the tambacu, with the metabolism category presenting the largest number of transcripts. The BUSCO analysis indicated that our assembly was more than 40% complete. We identified 21,986 genes for the two fishes. The P and Log2FC values indicated significant differences in the levels of gene expression, with a total of 600 up-regulated genes.

Conclusion

In spite of the lack of a reference genome, the functional annotation was successful, and confirmed a considerable difference in the specificity and levels of gene expression between the two organisms. This report provides a comprehensive baseline for the genetic management of these commercially important fishes, in particular for the identification of specific genes that may represent markers involved in the immunity, growth, and fertility of these organisms, with potential practical applications in aquaculture management.

De Novo Transcriptome Assembly and Population Genetic Analyses for an Endangered Chinese Endemic Acer miaotaiense (Aceraceae)

Acer miaotaiense (P. C. Tsoong) is a rare and highly endangered plant in China. Because of the lack of genomic information and the limited number of available molecular markers, there are insufficient tools to determine the genetic diversity of this species. Here, 93,305 unigenes were obtained by multiple assembled contigs with a transcriptome sequencing program. Furthermore, 12,819 expressed sequence tag-derived simple sequence repeat (EST-SSR) markers were generated, 300 were randomly selected and synthesized, 19 primer pairs were identified as highly polymorphic (average number of alleles (N_a) = 8, expected heterozygosity (H_e) = 0.635, polymorphism information content (PIC) = 0.604) and were further used for population genetic analysis. All 261 samples were grouped into two genetic clusters by UPGMA, a principal component analyses and a STRUCTURE analyses. A moderate level of genetic differentiation (genetic differentiation index (F_st) = 0.059–0.116, gene flow = 1.904–3.993) among the populations and the major genetic variance (81.01%) within populations were revealed by the AMOVA. Based on the results, scientific conservation strategies should be established using in situ and ex situ conservation strategies. The study provides useful genetic information for the protection of precious wild resources and for further research on the origin and evolution of this endangered plant and its related species.

De novo transcriptome sequencing of Camellia sasanqua and the analysis of major candidate genes related to floral traits

Camellia sasanqua is one of the most famous horticultural plants in Camellia (Theaceae) due to its aesthetic appeal as landscape plant. Knowledge regarding the genetic basis of flowering time, floral aroma and color in C. sasanqua is limited, but is essential to breed new varieties with desired floral traits. Here, we described the de novo transcriptome of young leaves, flower buds and flowers of C. sasanqua. A total of 60,127 unigenes were functionally annotated based on the sequence similarity. After analysis, we found that two floral integrator genes, SOC1 and AP1, in flowering time pathway showed evidence of gene family expansion. Compared with 1-deoxy-D-xylulose-5-phosphate pathway, some genes in the mevalonate pathway were most highly expressed, suggesting that this might represent the major pathway for terpenoid biosynthesis related to floral aroma in C. sasanqua. In flavonoid biosynthesis pathway, PAL, CHI, DFR and ANS showing significantly higher expression levels in flowers and flower buds might have important role in regulation of floral color. The top five most transcription factors (TFs) families in C. sasanqua transcriptome were MYB, MIKC, C3H, FAR1 and HD-ZIP, many of which have a direct relationship with floral traits. In addition, we also identified 33,540 simple sequence repeats (SSRs) in the C. sasanqua transcriptome. Collectively, the C. sasanqua transcriptome dataset generated from this study along with the SSR markers provide a new resource for the identification of novel regulatory transcripts and will accelerate the genetic improvement of C. sasanqua breeding programs.

De novo assembly, characterization and annotation for the transcriptome of Sphaeroma terebrans and microsatellite marker discovery

Sphaeroma terebrans, an economically and ecologically important marine wood-boring isopod, is mainly distributed in tropical and subtropical mangroves. Nevertheless, available genomic and transcriptomic information for this isopod is extremely deficient. Here, we first performed the assembly of S. terebrans transcriptome by Illumina sequencing. A total of 51,092 high-quality unigenes with an average length of 641 bp were obtained and 19,915 unigenes were successfully annotated in four public databases. Only 9932 out of 19915 unigenes were commonly annotated by all four databases. In addition, 9609 unigenes were categorized into 54 function categories of Gene Ontology (GO), and 14,512 unigenes were successfully grouped into 25 functional categories of the EuKaryotic Orthologous Groups (KOG) database. Moreover, 11,507 unigenes were assigned to 228 pathways of the Kyoto Encyclopedia of Genes and Genomes (KEGG). Out of 51,092 unigenes, a total of 4257 different microsatellites with motifs range from di- to hexa-nucleotide were identified from 3324 unigene sequences. Among 64 primer pairs selected for validation, 35 were successful in PCR amplification and 13 exhibited obvious repeat polymorphisms in the wild population of S. terebrans in Dongzhaigang (Hainan Island). The transcriptome dataset and the identified microsatellite markers (SSRs) will provide abundant information for researches on the discovery of new genes, metabolic mechanism and genetic diversity of S. terebrans.

Development and Evaluation of a Novel Set of EST-SSR Markers Based on Transcriptome Sequences of Black Locust (Robinia pseudoacacia L.)

Black locust (Robinia pseudoacacia L. of the family Fabaceae) is an ecologically and economically important deciduous tree. However, few genomic resources are available for this forest species, and few effective expressed sequence tag-derived simple sequence repeat (EST-SSR) markers have been developed to date. In this study, paired-end sequencing was used to sequence transcriptomes of R. pseudoacacia by the Illumina HiSeq TM2000 platform, and EST-SSR loci were identified by de novo assembly. Furthermore, a total of 1697 primer pairs were successfully designed, from which 286 primers met the selection screening criteria; 94 pairs were randomly selected and tested for validation using polymerase chain reaction amplification. Forty-five primers were verified as polymorphic, with clear bands. The polymorphism information content values were 0.033–0.765, the number of alleles per locus ranged from 2 to 10, and the observed and expected heterozygosities were 0.000–0.931 and 0.035–0.810, respectively, indicating a high level of informativeness. Subsequently, 45 polymorphic EST-SSR loci were tested for amplification efficiency, using the verified primers, in an additional nine species of Leguminosae, 23 loci were amplified in more than three species, of which two loci were amplified successfully in all species. These EST-SSR markers provide a valuable tool for investigating the genetic diversity and population structure of R. pseudoacacia, constructing a DNA fingerprint database, performing quantitative trait locus mapping, and preserving genetic information.

Seasonal differences in the testicular transcriptome profile of free-living European beavers (Castor fiber L.) determined by the RNA-Seq method

The European beaver (Castor fiber L.) is an important free-living rodent that inhabits Eurasian temperate forests. Beavers are often referred to as ecosystem engineers because they create or change existing habitats, enhance biodiversity and prepare the environment for diverse plant and animal species. Beavers are protected in most European Union countries, but their genomic background remains unknown. In this study, gene expression patterns in beaver testes and the variations in genetic expression in breeding and non-breeding seasons were determined by high-throughput transcriptome sequencing. Paired-end sequencing in the Illumina HiSeq 2000 sequencer produced a total of 373.06 million of high-quality reads. De novo assembly of contigs yielded 130,741 unigenes with an average length of 1,369.3 nt, N50 value of 1,734, and average GC content of 46.51%. A comprehensive analysis of the testicular transcriptome revealed more than 26,000 highly expressed unigenes which exhibited the highest homology with Rattus norvegicus and Ictidomys tridecemlineatus genomes. More than 8,000 highly expressed genes were found to be involved in fundamental biological processes, cellular components or molecular pathways. The study also revealed 42 genes whose regulation differed between breeding and non-breeding seasons. During the non-breeding period, the expression of 37 genes was up-regulated, and the expression of 5 genes was down-regulated relative to the breeding season. The identified genes encode molecules which are involved in signaling transduction, DNA repair, stress responses, inflammatory processes, metabolism and steroidogenesis. Our results pave the way for further research into season-dependent variations in beaver testes.

Genome-wide mapping and characterization of microsatellites in the swamp eel genome

We described genome-wide screening and characterization of microsatellites in the swamp eel genome. A total of 99,293 microsatellite loci were identified in the genome with an overall density of 179 microsatellites per megabase of genomic sequences. The dinucleotide microsatellites were the most abundant type representing 71% of the total microsatellite loci and the AC-rich motifs were the most recurrent in all repeat types. Microsatellite frequency decreased as numbers of repeat units increased, which was more obvious in long than short microsatellite motifs. Most of microsatellites were located in non-coding regions, whereas only approximately 1% of the microsatellites were detected in coding regions. Trinucleotide repeats were most abundant microsatellites in the coding regions, which represented amino acid repeats in proteins. There was a chromosome-biased distribution of microsatellites in non-coding regions, with the highest density of 203.95/Mb on chromosome 8 and the least on chromosome 7 (164.06/Mb). The most abundant dinucleotides (AC)n was mainly located on chromosome 8. Notably, genomic mapping showed that there was a chromosome-biased association of genomic distributions between microsatellites and transposon elements. Thus, the novel dataset of microsatellites in swamp eel provides a valuable resource for further studies on QTL-based selection breeding, genetic resource conservation and evolutionary genetics.