The transcriptome of the bowhead whale Balaena mysticetus reveals adaptations of the longest‐lived mammal

Mammals vary dramatically in lifespan, by at least two-orders of magnitude, but the molecular basis for this difference remains largely unknown. The bowhead whale Balaena mysticetus is the longest-lived mammal known, with an estimated maximal lifespan in excess of two hundred years. It is also one of the two largest animals and the most cold-adapted baleen whale species. Here, we report the first genome-wide gene expression analyses of the bowhead whale, based on the de novo assembly of its transcriptome. Bowhead whale or cetacean-specific changes in gene expression were identified in the liver, kidney and heart, and complemented with analyses of positively selected genes. Changes associated with altered insulin signaling and other gene expression patterns could help explain the remarkable longevity of bowhead whales as well as their adaptation to a lipid-rich diet. The data also reveal parallels in candidate longevity adaptations of the bowhead whale, naked mole rat and Brandt's bat. The bowhead whale transcriptome is a valuable resource for the study of this remarkable animal, including the evolution of longevity and its important correlates such as resistance to cancer and other diseases.

Total transcriptome, proteome, and allergome of Johnson grass pollen, which is important for allergic rhinitis in subtropical regions

Background Genomic data are lacking for many allergen sources. To circumvent this limitation, we implemented a strategy to reveal the repertoire of pollen allergens of a grass with clinical importance in subtropical regions, where an increasing proportion of the world's population resides. Objective We sought to identify and immunologically characterize the allergenic components of the Panicoideae Johnson grass pollen (JGP; Sorghum halepense). Methods The total pollen transcriptome, proteome, and allergome of JGP were documented. Serum IgE reactivities with pollen and purified allergens were assessed in 64 patients with grass pollen allergy from a subtropical region. Results Purified Sor h 1 and Sor h 13 were identified as clinically important allergen components of JGP with serum IgE reactivity in 49 (76%) and 28 (43.8%), respectively, of patients with grass pollen allergy. Within whole JGP, multiple cDNA transcripts and peptide spectra belonging to grass pollen allergen families 1, 2, 4, 7, 11, 12, 13, and 25 were identified. Pollen allergens restricted to subtropical grasses (groups 22-24) were also present within the JGP transcriptome and proteome. Mass spectrometry confirmed the IgE-reactive components of JGP included isoforms of Sor h 1, Sor h 2, Sor h 13, and Sor h 23. Conclusion Our integrated molecular approach revealed qualitative differences between the allergenic components of JGP and temperate grass pollens. Knowledge of these newly identified allergens has the potential to improve specific diagnosis and allergen immunotherapy treatment for patients with grass pollen allergy in subtropical regions and reduce the burden of allergic respiratory disease globally.

Combining transcriptome assemblies from multiple de novo assemblers in the allo-tetraploid plant Nicotiana benthamiana

Background Nicotiana benthamiana is an allo-tetraploid plant, which can be challenging for de novo transcriptome assemblies due to homeologous and duplicated gene copies. Transcripts generated from such genes can be distinct yet highly similar in sequence, with markedly differing expression levels. This can lead to unassembled, partially assembled or mis-assembled contigs. Due to the different properties of de novo assemblers, no one assembler with any one given parameter space can re-assemble all possible transcripts from a transcriptome. Results In an effort to maximise the diversity and completeness of de novo assembled transcripts, we utilised four de novo transcriptome assemblers, TransAbyss, Trinity, SOAPdenovo-Trans, and Oases, using a range of k-mer sizes and different input RNA-seq read counts. We complemented the parameter space biologically by using RNA from 10 plant tissues. We then combined the output of all assemblies into a large super-set of sequences. Using a method from the EvidentialGene pipeline, the combined assembly was reduced from 9.9 million de novo assembled transcripts to about 235,000 of which about 50,000 were classified as primary. Metrics such as average bit-scores, feature response curves and the ability to distinguish paralogous or homeologous transcripts, indicated that the EvidentialGene processed assembly was of high quality. Of 35 RNA silencing gene transcripts, 34 were identified as assembled to full length, whereas in a previous assembly using only one assembler, 9 of these were partially assembled. Conclusions To achieve a high quality transcriptome, it is advantageous to implement and combine the output from as many different de novo assemblers as possible. We have in essence taking the ‘best’ output from each assembler while minimising sequence redundancy. We have also shown that simultaneous assessment of a variety of metrics, not just focused on contig length, is necessary to gauge the quality of assemblies.

The multiple sclerosis whole blood mRNA transcriptome and genetic associations indicate dysregulation of specific T cell pathways in pathogenesis

Multiple sclerosis (MS) is an autoimmune disease with a genetic component, caused at least in part by aberrant lymphocyte activity. The whole blood mRNA transcriptome was measured for 99 untreated MS patients: 43 primary progressive MS, 20 secondary progressive MS, 36 relapsing remitting MS and 45 age-matched healthy controls. The ANZgene Multiple Sclerosis Genetics Consortium genotyped more than 300 000 SNPs for 115 of these samples. Transcription from genes on translational regulation, oxidative phosphorylation, immune synapse and antigen presentation pathways was markedly increased in all forms of MS. Expression of genes tagging T cells was also upregulated (P < 10-12) in MS. A T cell gene signature predicts disease state with a concordance index of 0.79 with age and gender as co-variables, but the signature is not associated with clinical course or disability. The ANZgene genome wide association screen identified two novel regions with genome wide significance: one encoding the T cell co-stimulatory molecule, CD40; the other a region on chromosome 12q13-14. The CD40 haplotype associated with increased MS susceptibility has decreased gene expression in MS (P < 0.0007). The second MS susceptibility region includes 17 genes on 12q13-14 in tight linkage disequilibrium. Of these, only 13 are expressed in leukocytes, and of these the expression of one, FAM119B, is much lower in the susceptibility haplotype (P tdthomlt; 10-14). Overall, these data indicate dysregulation of T cells can be detected in the whole blood of untreated MS patients, and supports targeting of activated T cells in therapy for all forms of MS.

Transcriptome analysis of Neotyphodium and Epichloë grass endophytes

Large-scale gene discovery has been performed for the grass fungal endophytes Neotyphodium coenophialum, Neotyphodium lolii, and Epichloë festucae. The resulting sequences have been annotated by comparison with public DNA and protein sequence databases and using intermediate gene ontology annotation tools. Endophyte sequences have also been analysed for the presence of simple sequence repeat and single nucleotide polymorphism molecular genetic markers. Sequences and annotation are maintained within a MySQL database that may be queried using a custom web interface. Two cDNA-based microarrays have been generated from this genome resource. They permit the interrogation of 3806 Neotyphodium genes (NchipTM microarray), and 4195 Neotyphodium and 920 Epichloë genes (EndoChipTM microarray), respectively. These microarrays provide tools for high-throughput transcriptome analysis, including genome-specific gene expression studies, profiling of novel endophyte genes, and investigation of the host grass–symbiont interaction. Comparative transcriptome analysis in Neotyphodium and Epichloë was performed

Genome-wide mapping of cystitis due to Streptococcus agalactiae and Escherichia coli in mice identifies a unique bladder transcriptome that signifies pathogen-specific antimicrobial defense against urinary tract infection

The most common causes of urinary tract infections (UTIs) are Gram-negative pathogens such as Escherichia coli; however, Gram-positive organisms including Streptococcus agalactiae, or group B streptococcus (GBS), also cause UTI. In GBS infection, UTI progresses to cystitis once the bacteria colonize bladder, but the host responses triggered in the bladder immediately following infection are largely unknown. Here, we used genome-wide expression profiling to map the bladder transcriptome of GBS UTI in mice infected transurethrally with uropathogenic GBS that was cultured from a 35 year-old women with cystitis. RNA from bladders was applied to Affymetrix Gene-1.0ST microarrays; qRT-PCR was used to analyze selected gene responses identified in array datasets. A surprisingly small significant gene list of 172 genes was identified at 24h; this compared to 2507 genes identified in a side-by-side comparison with uropathogenic E. coli (UPEC). No genes exhibited significantly altered expression at 2h in GBS-infected mice according to arrays despite high bladder bacterial loads at this early time point. The absence of a marked early host response to GBS juxtaposed with broad-based bladder responses activated by UPEC at 2h. Bioinformatics analyses including integrative systems-level network mapping revealed multiple activated biological pathways in the GBS cystitis transcriptome that regulate leukocyte activation, inflammation, apoptosis, and cytokine-chemokine biosynthesis. These findings define a novel, minimalistic type of bladder host response triggered by GBS UTI, which comprises collective antimicrobial pathways that differ dramatically from those activated by UPEC. Overall, this study emphasizes the unique nature of bladder immune activation mechanisms triggered by distinct uropathogens.

Sexual selection in true fruit flies (Diptera: Tephritidae): transcriptome and experimental evidences for phytochemicals increasing male competitive ability

In male tephritid fruit flies of the genus Bactrocera, feeding on secondary plant compounds (sensu lato male lures = methyl eugenol, raspberry ketone and zingerone) increases male mating success. Ingested male lures alter the male pheromonal blend, normally making it more attractive to females and this is considered the primary mechanism for the enhanced mating success. However, the male lures raspberry ketone and zingerone are known, across a diverse range of other organisms, to be involved in increasing energy metabolism. If this also occurs in Bactrocera, then this may represent an additional benefit to males as courtship is metabolically expensive and lure feeding may increase a fly's short-term energy. We tested this hypothesis by performing comparative RNA-seq analysis between zingerone-fed and unfed males of Bactrocera tryoni. We also carried out behavioural assays with zingerone- and cuelure-fed males to test whether they became more active. RNA-seq analysis revealed, in zingerone-fed flies, up-regulation of 3183 genes with homologues transcripts to those known to regulate intermale aggression, pheromone synthesis, mating and accessory gland proteins, along with significant enrichment of several energy metabolic pathways and gene ontology terms. Behavioural assays show significant increases in locomotor activity, weight reduction and successful mating after mounting; all direct/indirect measures of increased activity. These results suggest that feeding on lures leads to complex physiological changes, which result in more competitive males. These results do not negate the pheromone effect, but do strongly suggest that the phytochemical-induced sexual selection is governed by both female preference and male competitive mechanisms.

Transcriptome and glycome profiling of human hematopoietic CD133+ and CD34+ cells

New blood cells are continuously provided by self-renewing multipotent hematopoietic stem cells (HSC). The capacity of HSCs to regenerate the hematopoietic system is utilized in the treatment of patients with hematological malignancies. HSCs can be enriched using an antibody-based recognition of CD34 or CD133 glycoproteins on the cell surface. The CD133+ and CD34+ cells may have partly different roles in hematopoiesis. Furthermore, each cell has a glycome typical for that cell type. Knowledge of HSC glycobiology can be used to design therapeutic cells with improved cell proliferation or homing properties. The present studies characterize the global gene expression profile of human cord blood-derived CD133+ and CD34+ cells, and demonstrate the differences between CD133+ and CD34+ cell populations that may have an impact in transplantation when CD133+ and CD34+ selected cells are used. In addition, these studies unravel the glycome profile of primitive hematopoietic cells and reveal the transcriptional regulation of N-glycan biosynthesis in CD133+ and CD34+ cells. The gene expression profile of CD133+ cells represents 690 differentially expressed transcripts between CD133+ cells and CD133- cells. CD34+ cells have 620 transcripts differentially expressed when compared to CD34- cells. The integrated CD133+/CD34+ cell gene expression profiles proffer novel transcripts to specify HSCs. Furthermore, the differences between the gene expression profiles of CD133+ and CD34+ cells indicate differences in the transcriptional regulation of CD133+ and CD34+ cells. CD133+ cells express a lower number of hematopoietic lineage differentiation marker genes than CD34+ cells. The expression profiles suggest a more primitive nature of CD133+ cells. Moreover, CD133+ cells have characteristic glycome that differ from the glycome of CD133- cells. High mannose-type and biantennary complex-type N-glycans are enriched in CD133+ cells. N-glycosylation-related gene expression pattern of CD133+ cells identify the key genes regulating the CD133+ cell-specific glycosylation including the overexpression of MGAT2 and underexpression of MGAT4. The putative role of MAN1C1 in the increase of unprocessed high mannose-type N-glycans in CD133+ cells is also discussed. These studies provide new information on the characteristics of HSCs. Improved understanding of HSC biology can be used to design therapeutic cells with improved cell proliferation and homing properties. As a result, HSC engineering could further their clinical use.

Transcriptome and proteome analysis of xylose-metabolising Saccharomyces cerevisiae

Increasing concern about global climate warming has accelerated research into renewable energy sources that could replace fossil petroleum-based fuels and materials. Bioethanol production from cellulosic biomass by fermentation with baker s yeast Saccharomyces cerevisiae is one of the most studied areas in this field. The focus has been on metabolic engineering of S. cerevisiae for utilisation of the pentose sugars, in particular D-xylose that is abundant in the hemicellulose fraction of biomass. Introduction of a heterologous xylose-utilisation pathway into S. cerevisiae enables xylose fermentation, but ethanol yield and productivity do not reach the theoretical level. In the present study, transcription, proteome and metabolic flux analyses of recombinant xylose-utilising S. cerevisiae expressing the genes encoding xylose reductase (XR) and xylitol dehydrogenase (XDH) from Pichia stipitis and the endogenous xylulokinase were carried out to characterise the global cellular responses to metabolism of xylose. The aim of these studies was to find novel ways to engineer cells for improved xylose fermentation. The analyses were carried out from cells grown on xylose and glucose both in batch and chemostat cultures. A particularly interesting observation was that several proteins had post-translationally modified forms with different abundance in cells grown on xylose and glucose. Hexokinase 2, glucokinase and both enolase isoenzymes 1 and 2 were phosphorylated differently on the two different carbon sources studied. This suggests that phosphorylation of glycolytic enzymes may be a yet poorly understood means to modulate their activity or function. The results also showed that metabolism of xylose affected the gene expression and abundance of proteins in pathways leading to acetyl-CoA synthesis and altered the metabolic fluxes in these pathways. Additionally, the analyses showed increased expression and abundance of several other genes and proteins involved in cellular redox reactions (e.g. aldo-ketoreductase Gcy1p and 6-phosphogluconate dehydrogenase) in cells grown on xylose. Metabolic flux analysis indicated increased NADPH-generating flux through the oxidative part of the pentose phosphate pathway in cells grown on xylose. The most importantly, results indicated that xylose was not able to repress to the same extent as glucose the genes of the tricarboxylic acid and glyoxylate cycles, gluconeogenesis and some other genes involved in the metabolism of respiratory carbon sources. This suggests that xylose is not recognised as a fully fermentative carbon source by the recombinant S. cerevisiae that may be one of the major reasons for the suboptimal fermentation of xylose. The regulatory network for carbon source recognition and catabolite repression is complex and its functions are only partly known. Consequently, multiple genetic modifications and also random approaches would probably be required if these pathways were to be modified for further improvement of xylose fermentation by recombinant S. cerevisiae strains.

Transcriptome analysis and characterisation of gill-expressed carbonic anhydrase and other key osmoregulatory genes in freshwater crayfish Cherax quadricarinatus

The pH and salinity balance mechanisms of crayfish are controlled by a set of transport-related genes. We identified a set of the genes from the gill transcriptome from a freshwater crayfish Cherax quadricarinatus using the Illumina NGS-sequencing technology. We identified and characterized carbonic anhydrase (CA) genes and some other key genes involved in systematic acid-base balance and osmotic/ionic regulation. We also examined expression patterns of some of these genes across different sublethal pH levels [1]. A total of 72,382,710 paired-end Illumina reads were assembled into 36,128 contigs with an average length of 800 bp. About 37% of the contigs received significant BLAST hits and 22% were assigned gene ontology terms. These data will assist in further physiological-genomic studies in crayfish.

Studies of the Human Transcriptome

Gene expression is one of the most critical factors influencing the phenotype of a cell. As a result of several technological advances, measuring gene expression levels has become one of the most common molecular biological measurements to study the behaviour of cells. The scientific community has produced enormous and constantly increasing collection of gene expression data from various human cells both from healthy and pathological conditions. However, while each of these studies is informative and enlighting in its own context and research setup, diverging methods and terminologies make it very challenging to integrate existing gene expression data to a more comprehensive view of human transcriptome function. On the other hand, bioinformatic science advances only through data integration and synthesis. The aim of this study was to develop biological and mathematical methods to overcome these challenges and to construct an integrated database of human transcriptome as well as to demonstrate its usage. Methods developed in this study can be divided in two distinct parts. First, the biological and medical annotation of the existing gene expression measurements needed to be encoded by systematic vocabularies. There was no single existing biomedical ontology or vocabulary suitable for this purpose. Thus, new annotation terminology was developed as a part of this work. Second part was to develop mathematical methods correcting the noise and systematic differences/errors in the data caused by various array generations. Additionally, there was a need to develop suitable computational methods for sample collection and archiving, unique sample identification, database structures, data retrieval and visualization. Bioinformatic methods were developed to analyze gene expression levels and putative functional associations of human genes by using the integrated gene expression data. Also a method to interpret individual gene expression profiles across all the healthy and pathological tissues of the reference database was developed. As a result of this work 9783 human gene expression samples measured by Affymetrix microarrays were integrated to form a unique human transcriptome resource GeneSapiens. This makes it possible to analyse expression levels of 17330 genes across 175 types of healthy and pathological human tissues. Application of this resource to interpret individual gene expression measurements allowed identification of tissue of origin with 92.0% accuracy among 44 healthy tissue types. Systematic analysis of transcriptional activity levels of 459 kinase genes was performed across 44 healthy and 55 pathological tissue types and a genome wide analysis of kinase gene co-expression networks was done. This analysis revealed biologically and medically interesting data on putative kinase gene functions in health and disease. Finally, we developed a method for alignment of gene expression profiles (AGEP) to perform analysis for individual patient samples to pinpoint gene- and pathway-specific changes in the test sample in relation to the reference transcriptome database. We also showed how large-scale gene expression data resources can be used to quantitatively characterize changes in the transcriptomic program of differentiating stem cells. Taken together, these studies indicate the power of systematic bioinformatic analyses to infer biological and medical insights from existing published datasets as well as to facilitate the interpretation of new molecular profiling data from individual patients.

Studies on Causes of Multiple Sclerosis : From Genes to Transcriptome

Multiple sclerosis (MS) is the most common cause of neurological disability in young adults, affecting more than two million people worldwide. It manifests as a chronic inflammation in the central nervous system (CNS) and causes demyelination and neurodegeneration. Depending on the location of the demyelinated plaques and axonal loss, a variety of symptoms can be observed including deficits in vision, coordination, balance and movement. With a typical age of onset at 20-40 years, the social and economic impacts of MS on lives of the patients and their families are considerable. Unfortunately the current treatments are relatively inefficient and the development of more effective treatments has been impeded by our limited understanding of the causes and pathogenesis of MS. Risk of MS is higher in biological relatives of MS patients than in the general population. Twin and adoption studies have shown that familial clustering of MS is explained by shared genetic factors rather than by shared familial environment. While the involvement of the human leukocyte antigen (HLA) genes was first discovered four decades ago, additional genetic risk factors have only recently been identified through genome-wide association studies (GWAS). Current evidence suggests that MS is a highly polygenic disease with perhaps hundreds of common variants with relatively modest effects contributing to susceptibility. Despite extensive research, the majority of these risk factors still remain to be identified. In this thesis the aim was to identify novel genes and pathways involved in MS. Using genome-wide microarray technology, gene expression levels in peripheral blood mononuclear cells (PBMC) from 12 MS patients and 15 controls were profiled and more than 600 genes with altered expression in MS were identified. Three of five selected findings, DEFA1A3, LILRA4 and TNFRSF25, were successfully replicated in an independent sample. Increased expression of DEFA1A3 in MS is a particularly interesting observation, because its elevated levels have previously been reported also in several other autoimmune diseases. A systematic review of seven microarray studies was then performed leading to identification of 229 genes, in which either decreased or increased expression in MS had been reported in at least two studies. In general there was relatively little overlap across the experiments: 11 of the 229 genes had been reported in three studies and only HSPA1A in four studies. Nevertheless, these 229 genes were associated with several immunological pathways including interleukin pathways related to type 2 and type 17 helper T cells and regulatory T cells. However, whether these pathways are involved in causing MS or related to secondary processes activated after disease onset remains to be investigated. The 229 genes were also compared with loci identified in published MS GWASs. Single nucleotide polymorphisms (SNP) in 17 of the 229 loci had been reported to be associated with MS with P-value less than 0.0001 including variants in CXCR4 and SAPS2, which were the only loci where evidence for correlation between the associated variant and gene expression was found. The CXCR4 variant was further tested for association with MS in a large case-control sample and the previously reported suggestive association was replicated (P-value is 0.0004). Finally, common genetic variants in candidate genes, which had been selected on the basis of showing association with other autoimmune diseases (MYO9B) or showing differential expression in MS in our study (DEFA1A3, LILRA4 and TNFRSF25), were tested for association with MS, but no evidence of association was found. In conclusion, through a systematic review of genome-wide expression studies in MS we have identified several promising candidate genes and pathways for future studies. In addition, we have replicated a previously suggested association of a SNP variant upstream of CXCR4 with MS. Keywords: autoimmune disease, common variant, CXCR4, DEFA1A3, HSPA1A,gene expression, genetic association, GWAS, MS, multiple sclerosis, systematic review

Analysis of the Genome and Transcriptome of Cryptococcus neoformans var. grubii Reveals Complex RNA Expression and Microevolution Leading to Virulence Attenuation

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Profiling of Luteal Transcriptome during Prostaglandin F2-Alpha Treatment in Buffalo Cows: Analysis of Signaling Pathways Associated with Luteolysis

In several species including the buffalo cow, prostaglandin (PG) F-2 alpha is the key molecule responsible for regression of corpus luteum (CL). Experiments were carried out to characterize gene expression changes in the CL tissue at various time points after administration of luteolytic dose of PGF(2 alpha) in buffalo cows. Circulating progesterone levels decreased within 1 h of PGF(2 alpha) treatment and evidence of apoptosis was demonstrable at 18 h post treatment. Microarray analysis indicated expression changes in several of immediate early genes and transcription factors within 3 h of treatment. Also, changes in expression of genes associated with cell to cell signaling, cytokine signaling, steroidogenesis, PG synthesis and apoptosis were observed. Analysis of various components of LH/CGR signaling in CL tissues indicated decreased LH/CGR protein expression, pCREB levels and PKA activity post PGF(2 alpha) treatment. The novel finding of this study is the down regulation of CYP19A1 gene expression accompanied by decrease in expression of E-2 receptors and circulating and intra luteal E-2 post PGF(2 alpha) treatment. Mining of microarray data revealed several differentially expressed E-2 responsive genes. Since CYP19A1 gene expression is low in the bovine CL, mining of microarray data of PGF(2 alpha)-treated macaques, the species with high luteal CYP19A1 expression, showed good correlation between differentially expressed E-2 responsive genes between both the species. Taken together, the results of this study suggest that PGF(2 alpha) interferes with luteotrophic signaling, impairs intraluteal E-2 levels and regulates various signaling pathways before the effects on structural luteolysis are manifest.

Comparative sequence analyses of genome and transcriptome reveal novel transcripts and variants in the Asian elephant Elephas maximus

The Asian elephant Elephas maximus and the African elephant Loxodonta africana that diverged 5-7 million years ago exhibit differences in their physiology, behaviour and morphology. A comparative genomics approach would be useful and necessary for evolutionary and functional genetic studies of elephants. We performed sequencing of E. maximus and map to L. africana at similar to 15X coverage. Through comparative sequence analyses, we have identified Asian elephant specific homozygous, non-synonymous single nucleotide variants (SNVs) that map to 1514 protein coding genes, many of which are involved in olfaction. We also present the first report of a high-coverage transcriptome sequence in E. maximus from peripheral blood lymphocytes. We have identified 103 novel protein coding transcripts and 66-long non-coding (lnc)RNAs. We also report the presence of 181 protein domains unique to elephants when compared to other Afrotheria species. Each of these findings can be further investigated to gain a better understanding of functional differences unique to elephant species, as well as those unique to elephantids in comparison with other mammals. This work therefore provides a valuable resource to explore the immense research potential of comparative analyses of transcriptome and genome sequences in the Asian elephant.