BMC Genomics - Latest Articles

Comparative transcriptome analysis of obligately asexual and cyclically sexual rotifers reveals genes with putative functions in sexual reproduction, dormancy, and asexual egg production

Sara Hanson — 2013-06-19T00:00:00Z

Background: Sexual reproduction is a widely studied biological process because it is critically important to the genetics, evolution, and ecology of eukaryotes. Despite decades of study on this topic, no comprehensive explanation has been accepted that explains the evolutionary forces underlying its prevalence and persistence in nature. Monogonont rotifers offer a useful system for experimental studies relating to the evolution of sexual reproduction due to their rapid reproductive rate and close relationship to the putatively ancient asexual bdelloid rotifers. However, little is known about the molecular underpinnings of sex in any rotifer species. Results: We generated mRNA-seq libraries for obligate parthenogenetic (OP) and cyclical parthenogenetic (CP) strains of the monogonont rotifer, Brachionus calyciflorus, to identify genes specific to both modes of reproduction. Our differential expression analysis identified receptors with putative roles in signaling pathways responsible for the transition from asexual to sexual reproduction. Differential expression of a specific copy of the duplicated cell cycle regulatory gene CDC20 and specific copies of histone H2A suggest that such duplications may underlie the phenotypic plasticity required for reproductive mode switch in monogononts. We further identified differential expression of genes involved in the formation of resting eggs, a process linked exclusively to sex in this species. Finally, we identified transcripts from the bdelloid rotifer Adineta ricciae that have significant sequence similarity to genes with higher expression in CP strains of B. calyciflorus. Conclusions: Our analysis of global gene expression differences between facultatively sexual and exclusively asexual populations of B. calyciflorus provides insights into the molecular nature of sexual reproduction in rotifers. Furthermore, our results offer insight into the evolution of obligate asexuality in bdelloid rotifers and provide indicators important for the use of monogononts as a model system for investigating the evolution of sexual reproduction.

Comparative transcriptome analysis of microsclerotia development in Nomuraea rileyi

Zhangyong Song — 2013-06-19T00:00:00Z

Background: Nomuraea rileyi is used as an environmental-friendly biopesticide. However, mass production and commercialization of this organism are limited due to its fastidious growth and sporulation requirements. When cultured in amended medium, we found that N. rileyi could produce microsclerotia bodies, replacing conidiophores as the infectious agent. However, little is known about the genes involved in microsclerotia development. In the present study, the transcriptomes were analyzed using next-generation sequencing technology to find the genes involved in microsclerotia development. Results: A total of 4.69 Gb of clean nucleotides comprising 32,061 sequences was obtained, and 20,919 sequences were annotated (about 65%). Among the annotated sequences, only 5928 were annotated with 34 gene ontology (GO) functional categories, and 12,778 sequences were mapped to 165 pathways by searching against the Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) database. Furthermore, we assessed the transcriptomic differences between cultures grown in minimal and amended medium. In total, 4808 sequences were found to be differentially expressed; 719 differentially expressed unigenes were assigned to 25 GO classes and 1888 differentially expressed unigenes were assigned to 161 KEGG pathways, including 25 enrichment pathways. Subsequently, we examined the up-regulation or uniquely expressed genes following amended medium treatment, which were also expressed on the enrichment pathway, and found that most of them participated in mediating oxidative stress homeostasis. To elucidate the role of oxidative stress in microsclerotia development, we analyzed the diversification of unigenes using quantitative reverse transcription-PCR (RT-qPCR). Conclusion: Our findings suggest that oxidative stress occurs during microsclerotia development, along with a broad metabolic activity change. Our data provide the most comprehensive sequence resource available for the study of N. rileyi. We believe that the transcriptome datasets will serve as an important public information platform to accelerate studies on N. rileyi microsclerotia.

Phage cluster relationships identified through single gene analysis

Kyle Smith — 2013-06-19T00:00:00Z

Background: Phylogenetic comparison of bacteriophages requires whole genome approaches such as dotplot analysis, genome pairwise maps, and gene content analysis. Currently mycobacteriophages, a highly studied phage group, are categorized into related clusters based on the comparative analysis of whole genome sequences. With the recent explosion of phage isolation, a simple method for phage cluster prediction would facilitate analysis of crude or complex samples without whole genome isolation and sequencing. The hypothesis of this study was that mycobacteriophage-cluster prediction is possible using comparison of a single, ubiquitous, semi-conserved gene. Tape Measure Protein (TMP) was selected to test the hypothesis because it is typically the longest gene in mycobacteriophage genomes and because regions within the TMP gene are conserved. Results: A single gene, TMP, identified the known Mycobacteriophage clusters and subclusters using a Gepard dotplot comparison or a phylogenetic tree constructed from global alignment and maximum likelihood comparisons. Gepard analysis of 247 mycobacteriophage TMP sequences appropriately recovered 98.8% of the subcluster assignments that were made by whole-genome comparison. Subcluster-specific primers within TMP allow for PCR determination of the mycobacteriophage subcluster from DNA samples. Using the single-gene comparison approach for siphovirus coliphages, phage groupings by TMP comparison reflected relationships observed in a whole genome dotplot comparison and confirm the potential utility of this approach to another widely studied group of phages. Conclusions: TMP sequence comparison and PCR results support the hypothesis that a single gene can be used for distinguishing phage cluster and subcluster assignments. TMP single-gene analysis can quickly and accurately aid in mycobacteriophage classification.

The complete mitochondrial genome of Solemya velum (Mollusca: Bivalvia) and its relationships with Conchifera

Federico Plazzi — 2013-06-18T00:00:00Z

Background: Bivalve mitochondrial genomes exhibit a wide array of uncommon features, like extensive gene rearrangements, large sizes, and unusual ways of inheritance. Species pertaining to the order Solemyida (subclass Opponobranchia) show many peculiar evolutionary adaptations, f.i. extensive symbiosis with chemoautotrophic bacteria. Despite Opponobranchia are central in bivalve phylogeny, being considered the sister group of all Autobranchia, a complete mitochondrial genome has not been sequenced yet. Results: In this paper, we characterized the complete mitochondrial genome of the Atlantic awning clam Solemya velum: A-T content, gene arrangement and other features are more similar to putative ancestral mollusks than to other bivalves. Two supranumerary open reading frames are present in a large, otherwise unassigned, region, while the origin of replication could be located in a region upstream to the cox3 gene. Conclusions: We show that S. velum mitogenome retains most of the ancestral conchiferan features, which is unusual among bivalve mollusks, and we discuss main peculiarities of this first example of an organellar genome coming from the subclass Opponobranchia. Mitochondrial genomes of Solemya (for bivalves) and Haliotis (for gastropods) seem to retain the original condition of mollusks, as most probably exemplified by Katharina.

Salmon lice (Lepeophtheirus salmonis) showing varying emamectin benzoate susceptibilities differ in neuronal acetylcholine receptor and GABA-gated chloride channel mRNA expression

Stephen Carmichael — 2013-06-18T00:00:00Z

Background: Caligid copepods, also called sea lice, are fish ectoparasites, some species of which cause significant problems in the mariculture of salmon, where the annual cost of infection is in excess of [euro sign]300 million globally. At present, caligid control on farms is mainly achieved using medicinal treatments. However, the continued use of a restricted number of medicine actives potentially favours the development of drug resistance. Here, we report transcriptional changes in a laboratory strain of the caligid Lepeophtheirus salmonis (Kr[latin small letter o with stroke]yer, 1837) that is moderately (~7-fold) resistant to the avermectin compound emamectin benzoate (EMB), a component of the anti-salmon louse agent SLICE(R) (Merck Animal Health). Results: Suppression subtractive hybridisation (SSH) was used to enrich transcripts differentially expressed between EMB-resistant (PT) and drug-susceptible (S) laboratory strains of L. salmonis. SSH libraries were subjected to 454 sequencing. Further L. salmonis transcript sequences were available as expressed sequence tags (EST) from GenBank. Contiguous sequences were generated from both SSH and EST sequences and annotated. Transcriptional responses in PT and S salmon lice were investigated using custom 15 K oligonucleotide microarrays designed using the above sequence resources. In the absence of EMB exposure, 359 targets differed in transcript abundance between the two strains, these genes being enriched for functions such as calcium ion binding, chitin metabolism and muscle structure. gamma-aminobutyric acid (GABA)-gated chloride channel (GABA-Cl) and neuronal acetylcholine receptor (nAChR) subunits showed significantly lower transcript levels in PT lice compared to S lice. Using RT-qPCR, the decrease in mRNA levels was estimated at ~1.4-fold for GABA-Cl and ~2.8-fold for nAChR. Salmon lice from the PT strain showed few transcriptional responses following acute exposure (1 or 3 h) to 200 mug L-1 of EMB, a drug concentration tolerated by PT lice, but toxic for S lice. Conclusions: Avermectins are believed to exert their toxicity to invertebrates through interaction with glutamate-gated and GABA-gated chloride channels. Further potential drug targets include other Cys-loop ion channels such as nAChR. The present study demonstrates decreased transcript abundances of GABA-Cl and nAChR subunits in EMB-resistant salmon lice, suggesting their involvement in avermectin toxicity in caligids.

Transcriptome sequencing and de novo annotation of the critically endangered Adriatic sturgeon

Michele Vidotto — 2013-06-18T00:00:00Z

Background: Sturgeons are a group of Condrostean fish with very high evolutionary, economical and conservation interest. The eggs of these living fossils represent one of the most high prized foods of animal origin. The intense fishing pressure on wild stocks to harvest caviar has caused in the last decades a dramatic decline of their distribution and abundance leading the International Union for Conservation of Nature to list them as the more endangered group of species. As a direct consequence, world-wide efforts have been made to develop sturgeon aquaculture programmes for caviar production. In this context, the characterization of the genes involved in sex determination could provide relevant information for the selective farming of the more profitable females. Results: The 454 sequencing of two cDNA libraries from the gonads and brain of one male and one female full-sib A. naccarii, yielded 182,066 and 167,776 reads respectively, which, after strict quality control, were iterative assembled into more than 55,000 high quality ESTs. The average per-base coverage reached by assembling the two libraries was 4X. The multi-step annotation process resulted in 16% successfully annotated sequences with GO terms. We screened the transcriptome for 32 sex-related genes and highlighted 7 genes that are potentially specifically expressed, 5 in male and 2 in females, at the first life stage at which sex is histologically identifiable. In addition we identified 21,791 putative EST-linked SNPs and 5,295 SSRs. Conclusions: This study represents the first large massive release of sturgeon transcriptome information that we organized into the public database AnaccariiBase, which is freely available at http://compgen.bio.unipd.it/anaccariibase/. This transcriptomic data represents an important source of information for further studies on sturgeon species. The hundreds of putative EST-linked molecular makers discovered in this study will be invaluable for sturgeon reintroduction and breeding programs.

Combined methylation mapping of 5mC and 5hmC during early embryonic stages in bovine

Béatrice de Montera — 2013-06-18T00:00:00Z

Background: It was recently established that changes in methylation during development are dynamic and involve both methylation and demethylation processes. Yet, which genomic sites are changing and what are the contributions of methylation (5mC) and hydroxymethylation (5hmC) to this epigenetic remodeling is still unknown. When studying early development, options for methylation profiling are limited by the unavailability of sufficient DNA material from these scarce samples and limitations are aggravated in non-model species due to the lack of technological platforms. We therefore sought to obtain a representation of differentially 5mC or 5hmC loci during bovine early embryo stages through the use of three complementary methods, based on selective methyl-sensitive restriction and enrichment by ligation-mediated PCR or on subtractive hybridization. Using these strategies, libraries of putative methylation and hydroxymethylated sites were generated from Day-7 and Day-12 bovine embryos. Results: Over 1.2 million sequencing reads were analyzed, resulting in 151,501 contigs, of which 69,136 were uniquely positioned on the genome. A total of 101,461 putative methylated sites were identified. The output of the three methods differed in genomic coverage as well as in the nature of the identified sites. The classical MspI/HpaII combination of restriction enzymes targeted CpG islands whereas the other methods covered 5mC and 5hmC sites outside of these regions. Data analysis suggests a transition of these methylation marks between Day-7 and Day-12 embryos in specific classes of repeat-containing elements. Conclusions: Our combined strategy offers a genomic map of the distribution of cytosine methylation/hydroxymethylation during early bovine embryo development. These results support the hypothesis of a regulatory phase of hypomethylation in repeat sequences during early embryogenesis.

Cubilin expression is monoallelic and epigenetically augmented via PPARs

Obaidullah Aseem — 2013-06-18T00:00:00Z

Background: Cubilin is an endocytic receptor that is necessary for renal and intestinal absorption of a range of ligands. Endocytosis mediated by cubilin and its co-receptor megalin is the principal mechanism for proximal tubule reabsorption of proteins from the glomerular filtrate. Cubilin is also required for intestinal endocytosis of intrinsic factor-vitamin B12 complex. Despite its importance, little is known about the regulation of cubilin expression. Results: Here we show that cubilin expression is under epigenetic regulation by at least two processes. The first process involves inactivation of expression of one of the cubilin alleles. This monoallelic expression state could not be transformed to biallelic by inhibiting DNA methylation or histone deacetylation. The second process involves transcriptional regulation of cubilin by peroxisome proliferator-activated receptor (PPAR) transcription factors that are themselves regulated by DNA methylation and histone deacetylation. This is supported by findings that inhibitors of DNA methylation and histone deacetylation, 5Aza and TSA, increase cubilin mRNA and protein in renal and intestinal cell lines. Not only was the expression of PPARalpha and gamma inducible by 5Aza and TSA, but the positive effects of TSA and 5Aza on cubilin expression were also dependent on both increased PPAR transcription and activation. Additionally, 5Aza and TSA had similar effects on the expression of the cubilin co-receptor, megalin. Conclusions: Together, these findings reveal that cubilin and megalin mRNA expression is under epigenetic control and thus point to new avenues for overcoming pathological suppression of these genes through targeting of epigenetic regulatory processes.

Genetic heterogeneity revealed by sequence analysis of Mycobacterium tuberculosis isolates from extra-pulmonary tuberculosis patients

Sarbashis Das — 2013-06-17T00:00:00Z

Background: Tuberculosis remains a major public health problem. Clinical tuberculosis manifests often as pulmonary and occasionally as extra-pulmonary tuberculosis. The emergence of drug resistant tubercle bacilli and its association with HIV is a formidable challenge to curb the spread of tuberculosis. There have been concerted efforts by whole genome sequencing and bioinformatics analysis to identify genomic patterns and to establish a relationship between the genotype of the organism and clinical manifestation of tuberculosis. Extra-pulmonary TB constitutes 15--20 percent of the total clinical cases of tuberculosis reported among immunocompetent patients, whereas among HIV patients the incidence is more than 50 percent. Genomic analysis of M. tuberculosis isolates from extra pulmonary patients has not been explored. Results: The genomic DNA of 5 extra-pulmonary clinical isolates of M. tuberculosis derived from cerebrospinal fluid, lymph node fine needle aspirates (FNAC) / biopsies, were sequenced. Next generation sequencing approach (NGS) was employed to identify Single Nucleotide Variations (SNVs) and computational methods used to predict their consequence on functional genes. Analysis of distribution of SNVs led to the finding that there are mixed genotypes in patient isolates and that many SNVs are likely to influence either gene function or their expression. Phylogenetic relationship between the isolates correlated with the origin of the isolates. In addition, insertion sites of IS elements were identified and their distribution revealed a variation in number and position of the element in the 5 extra-pulmonary isolates compared to the reference M. tuberculosis H37Rv strain. Conclusions: The results suggest that NGS sequencing is able to identify small variations in genomes of M. tuberculosis isolates including changes in IS element insertion sites. Moreover, variations in isolates of M. tuberculosis from non-pulmonary sites were documented. The analysis of our results indicates genomic heterogeneity in the clinical isolates.

A unique nucleosome arrangement, maintained actively by chromatin remodelers facilitates transcription of yeast tRNA genes

Yatendra Kumar — 2013-06-17T00:00:00Z

Background: RNA polymerase (pol) III transcribes a unique class of genes with intra-genic promoters and high transcriptional activity. The major contributors to the pol III transcriptome, tRNAs genes are found scattered on all chromosomes of yeast. A prototype tDNA of <150 bp length, is generally considered nucleosome-free while some pol III-transcribed genes have been shown to have nucleosome-positioning properties. Results: Using high resolution ChIP-chip and ChIP-seq methods, we found several unique features associated with nucleosome profiles on all tRNA genes of budding yeast, not seen on nucleosome-dense counterparts in fission yeast and resting human CD4+ T cells. The nucleosome-free region (NFR) on all but three yeast tDNAs is found bordered by an upstream (US) nucleosome strongly positioned at -140 bp position and a downstream (DS) nucleosome at variable positions with respect to the gene terminator. Perturbation in this nucleosomal arrangement interferes with the tRNA production. Three different chromatin remodelers generate and maintain the NFR by targeting different gene regions. Isw1 localizes to the gene body and makes it nucleosome-depleted, Isw2 maintains periodicity in the upstream nucleosomal array, while RSC targets the downstream nucleosome. Direct communication of pol III with RSC serves as a stress-sensory mechanism for these genes. In its absence, the downstream nucleosome moves towards the gene terminator. Levels of tRNAs from different families are found to vary considerably as different pol III levels are seen even on isogenes within a family. Pol III levels show negative correlation with the nucleosome occupancies on different genes. Conclusions: Budding yeast tRNA genes maintain an open chromatin structure, which is not due to sequence-directed nucleosome positioning or high transcription activity of genes. Unlike 5[prime] NFR on pol II-transcribed genes, the tDNA NFR, which facilitates tDNA transcription, results from action of chromatin remodeler Isw1, aided by Isw2 and RSC. The RSC-regulated nucleosome dynamics at the 3[prime] gene-end serves as a novel regulatory mechanism for pol III transcription in vivo, probably by controlling terminator-dependent facilitated recycling of pol III. Salient features of yeast tDNA chromatin structure reported in this study can explain the basis of the novel non-transcriptional roles ascribed to tDNAs.