Search

Abstract Background Zygotic transcription in fish embryos initiates around the time of gastrulation, and all prior development is initiated and controlled by maternally derived messenger RNAs. Atlantic cod egg and embryo viability is variable, and it is hypothesized that the early development depends upon the feature of these maternal RNAs. Both the length and the presence of specific motifs in the 3’UTR of maternal RNAs are believed to regulate expression and stability of the maternal transcripts. Therefore, the aim of this study was to characterize the overall composition and 3’UTR structure of the most common maternal RNAs found in cod eggs and pre-zygotic embryos. Results 22229 Sanger-sequences were obtained from 3’-end sequenced cDNA libraries prepared from oocyte, 1-2 cell, blastula and gastrula stages. Quantitative PCR revealed that EST copy number below 9 did not reflect the gene expression profile. Consequently genes represented by less than 9 ESTs were excluded from downstream analyses, in addition to sequences with low-quality gene hits. This provided 12764 EST sequences, encoding 257 unique genes, for further analysis. Mitochondrial transcripts accounted for 45.9-50.6% of the transcripts isolated from the maternal stages, but only 12.2% of those present at the onset of zygotic transcription. 3’UTR length was predicted in nuclear sequences with poly-A tail, which identified 191 3’UTRs. Their characteristics indicated a more complex regulation of transcripts that are abundant prior to the onset of zygotic transcription. Maternal and stable transcripts had longer 3’UTR (mean 187.1 and 208.8 bp) and more 3’UTR isoforms (45.7 and 34.6%) compared to zygotic transcripts, where 15.4% had 3’UTR isoforms and the mean 3’UTR length was 76 bp. Also, diversity and the amount of putative polyadenylation motifs were higher in both maternal and stable transcripts. Conclusions We report on the most pronounced processes in the maternally transferred cod transcriptome. Maternal stages are characterized by a rich abundance of mitochondrial transcripts. Maternal and stable transcripts display longer 3'UTRs with more variation of both polyadenylation motifs and 3'UTR isoforms. These data suggest that cod eggs possess a complex array of maternal RNAs which likely act to tightly regulate early developmental processes in the newly fertilized egg.

Abstract Background The complex genome of rapeseed (Brassica napus) is not well understood despite the economic importance of the species. Good knowledge of sequence variation is needed for genetics approaches and breeding purposes. We used a diversity set of B. napus representing eight different germplasm types to sequence genome-wide distributed restriction-site associated DNA (RAD) fragments for polymorphism detection and genotyping. Results More than 113,000 RAD clusters with more than 20,000 single nucleotide polymorphisms (SNPs) and 125 insertions/deletions were detected and characterized. About one third of the RAD clusters and polymorphisms mapped to the Brassica rapa reference sequence. An even distribution of RAD clusters and polymorphisms was observed across the B. rapa chromosomes, which suggests that there might be an equal distribution over the Brassica oleracea chromosomes, too. The representation of Gene Ontology (GO) terms for unigenes with RAD clusters and polymorphisms revealed no signature of selection with respect to the distribution of polymorphisms within genes belonging to a specific GO category. Conclusions Considering the decreasing costs for next-generation sequencing, the results of our study suggest that RAD sequencing is not only a simple and cost-effective method for high-density polymorphism detection but also an alternative to SNP genotyping from transcriptome sequencing or SNP arrays, even for species with complex genomes such as B. napus.

Abstract Background The Manila clam, Ruditapes philippinarum, is one of the major aquaculture species in the world and a potential sentinel organism for monitoring the status of marine ecosystems. However, genomic resources for R. philippinarum are still extremely limited. Global analysis of gene expression profiles is increasingly used to evaluate the biological effects of various environmental stressors on aquatic animals under either artificial conditions or in the wild. Here, we report on the development of a transcriptomic platform for global gene expression profiling in the Manila clam. Results A normalized cDNA library representing a mixture of adult tissues was sequenced using a ultra high-throughput sequencing technology (Roche 454). A database consisting of 32,606 unique transcripts was constructed, 9,747 (30%) of which could be annotated by similarity. An oligo-DNA microarray platform was designed and applied to profile gene expression of digestive gland and gills. Functional annotation of differentially expressed genes between different tissues was performed by enrichment analysis. Expression of Natural Antisense Transcripts (NAT) analysis was also performed and bi-directional transcription appears a common phenomenon in the R. philippinarum transcriptome. A preliminary study on clam samples collected in a highly polluted area of the Venice Lagoon demonstrated the applicability of genomic tools to environmental monitoring. Conclusions The transcriptomic platform developed for the Manila clam confirmed the high level of reproducibility of current microarray technology. Next-generation sequencing provided a good representation of the clam transcriptome. Despite the known limitations in transcript annotation and sequence coverage for non model species, sufficient information was obtained to identify a large set of genes potentially involved in cellular response to environmental stress.

Abstract Background The bumblebee Bombus terrestris is an ecologically and economically important pollinator and has become an important biological model system. To study fundamental evolutionary questions at the genomic level, a high resolution genetic linkage map is an essential tool for analyses ranging from quantitative trait loci (QTL) mapping to genome assembly and comparative genomics. We here present a saturated linkage map and match it with the Apis mellifera genome using homologous markers. This genome-wide comparison allows insights into structural conservations and rearrangements and thus the evolution on a chromosomal level. Results The high density linkage map covers ~ 93% of the B. terrestris genome on 18 linkage groups (LGs) and has a length of 2'047 cM with an average marker distance of 4.02 cM. Based on a genome size of ~ 430 Mb, the recombination rate estimate is 4.76 cM/Mb. Sequence homologies of 242 homologous markers allowed to match 15 B. terrestris with A. mellifera LGs, five of them as composites. Comparing marker orders between both genomes we detect over 14% of the genome to be organized in synteny and 21% in rearranged blocks on the same homologous LG. Conclusions This study demonstrates that, despite the very high recombination rates of both A. mellifera and B. terrestris and a long divergence time of about 100 million years, the genomes' genetic architecture is highly conserved. This reflects a slow genome evolution in these bees. We show that data on genome organization and conserved molecular markers can be used as a powerful tool for comparative genomics and evolutionary studies, opening up new avenues of research in the Apidae.

Abstract Background Calcium ion is tightly regulated in body fluids and for euryhaline fish, which are exposed to rapid changes in environmental [Ca2+], homeostasis is especially challenging. The gill is the main organ of active calcium uptake and therefore plays a crucial role in the maintenance of calcium ion homeostasis. To study the molecular basis of the short-term responses to changing calcium availability, the whole gill transcriptome obtained by Super Serial Analysis of Gene Expression (SuperSAGE) of the euryhaline teleost green spotted puffer fish, Tetraodon nigroviridis, exposed to water with altered [Ca2+] was analysed. Results Transfer of T. nigroviridis from 10 ppt water salinity containing 2.9 mM Ca2+ to high (10 mM Ca2+ ) and low (0.01 mM Ca2+) calcium water of similar salinity for 2-12 h resulted in 1,339 differentially expressed SuperSAGE tags (26-bp transcript identifiers) in gills. Of these 869 tags (65%) were mapped to T. nigroviridis cDNAs or genomic DNA and 497 (57%) were assigned to known proteins. Thirteen percent of the genes matched multiple tags indicating alternative RNA transcripts. The main enriched gene ontology groups belong to Ca2+ signaling/homeostasis but also muscle contraction, cytoskeleton, energy production/homeostasis and tissue remodeling. K-means clustering identified co-expressed transcripts with distinct patterns in response to water [Ca2+] and exposure time. Conclusions The generated transcript expression patterns provide a framework of novel water calcium-responsive genes in the gill during the initial response after transfer to different [Ca2+]. This molecular response entails initial perception of alterations, activation of signaling networks and effectors and suggests active remodeling of cytoskeletal proteins during the initial acclimation process. Genes related to energy production and energy homeostasis are also up-regulated, probably reflecting the increased energetic needs of the acclimation response. This study is the first genome-wide transcriptome analysis of fish gills and is an important resource for future research on the short-term mechanisms involved in the gill acclimation responses to environmental Ca2+ changes and osmoregulation.

Abstract Background The European sea bass (Dicentrarchus labrax) is a marine fish of great importance for fisheries and aquaculture. Functional genomics offers the possibility to discover the molecular mechanisms underlying productive traits in farmed fish, and a step towards the application of marker assisted selection methods in this species. To this end, we report here on the development of an oligo DNA microarray for D. labrax. Results A database consisting of 19,048 unique transcripts was constructed, of which 12,008 (63%) could be annotated by similarity and 4,692 received a GO functional annotation. Two non-overlapping 60mer probes were designed for each unique transcript and in-situ synthesized on glass slides using Agilent SurePrint™ technology. Probe design was positively completed for 19,035 target clusters; the oligo microarray was then applied to profile gene expression in mandibles and whole-heads of fish affected by prognathism, a skeletal malformation that strongly affects sea bass production. Statistical analysis identified 242 transcripts that are significantly down-regulated in deformed individuals compared to normal fish, with a significant enrichment in genes related to nervous system development and functioning. A set of genes spanning a wide dynamic range in gene expression level were selected for quantitative RT-PCR validation. Fold change correlation between microarray and qPCR data was always significant. Conclusions The microarray platform developed for the European sea bass has a high level of flexibility, reliability, and reproducibility. Despite the well known limitations in achieving a proper functional annotation in non-model species, sufficient information was obtained to identify biological processes that are significantly enriched among differentially expressed genes. New insights were obtained on putative mechanisms involved on mandibular prognathism, suggesting that bone/nervous system development might play a role in this phenomenon.

Abstract Background The genus Erwinia includes plant-associated pathogenic and non-pathogenic Enterobacteria. Important pathogens such as Erwinia amylovora, the causative agent of fire blight and E. pyrifoliae causing bacterial shoot blight of pear in Asia belong to this genus. The species E. tasmaniensis and E. billingiae are epiphytic bacteria and may represent antagonists for biocontrol of fire blight. The presence of genes that are putatively involved in virulence in E. amylovora and E. pyrifoliae is of special interest for these species in consequence. Results Here we provide the complete genome sequences of the pathogenic E. pyrifoliae strain Ep1/96 with a size of 4.1 Mb and of the non-pathogenic species E. billingiae strain Eb661 with a size of 5.4 Mb, de novo determined by conventional Sanger sequencing and next generation sequencing techniques. Genome comparison reveals large inversions resulting from homologous recombination events. Furthermore, comparison of deduced proteins highlights a relation of E. billingiae strain Eb661 to E. tasmaniensis strain Et1/99 and a distance to E. pyrifoliae for the overall gene content as well as for the presence of encoded proteins representing virulence factors for the pathogenic species. Pathogenicity of E. pyrifoliae is supposed to have evolved by accumulation of potential virulence factors. E. pyrifoliae carries factors for type III secretion and cell invasion. Other genes described as virulence factors for E. amylovora are involved in the production of exopolysaccharides, the utilization of plant metabolites such as sorbitol and sucrose. Some virulence-associated genes of the pathogenic species are present in E. tasmaniensis but mostly absent in E. billingiae. Conclusion The data of the genome analyses correspond to the pathogenic lifestyle of E. pyrifoliae and underlines the epiphytic localization of E. tasmaniensis and E. billingiae as a saprophyte.

Abstract Background The bumblebee, Bombus terrestris (Order Hymenoptera), is of widespread importance. This species is extensively used for commercial pollination in Europe, and along with other Bombus spp. is a key member of natural pollinator assemblages. Furthermore, the species is studied in a wide variety of biological fields. The objective of this project was to create a B. terrestris EST resource that will prove to be valuable in obtaining a deeper understanding of this significant social insect. Results A normalised cDNA library was constructed from the thorax and abdomen of B. terrestris workers in order to enhance the discovery of rare genes. A total of 29'428 ESTs were sequenced. Subsequent clustering resulted in 13'333 unique sequences. Of these, 58.8 percent had significant similarities to known proteins, with 54.5 percent having a "best-hit" to existing Hymenoptera sequences. Comparisons with the honeybee and other insects allowed the identification of potential candidates for gene loss, pseudogene evolution, and possible incomplete annotation in the honeybee genome. Further, given the focus of much basic research and the perceived threat of disease to natural and commercial populations, the immune system of bumblebees is a particularly relevant component. Although the library is derived from unchallenged bees, we still uncover transcription of a number of immune genes spanning the principally described insect immune pathways. Additionally, the EST library provides a resource for the discovery of genetic markers that can be used in population level studies. Indeed, initial screens identified 589 simple sequence repeats and 854 potential single nucleotide polymorphisms. Conclusion The resource that these B. terrestris ESTs represent is valuable for ongoing work. The ESTs provide direct evidence of transcriptionally active regions, but they will also facilitate further functional genomics, gene discovery and future genome annotation. These are important aspects in obtaining a greater understanding of this key pollinator species.

Abstract Background Food supply from the ocean is constrained by the shortage of domesticated and selected fish. Development of genomic models of economically important fishes should assist with the removal of this bottleneck. European sea bass Dicentrarchus labrax L. (Moronidae, Perciformes, Teleostei) is one of the most important fishes in European marine aquaculture; growing genomic resources put it on its way to serve as an economic model. Results End sequencing of a sea bass genomic BAC-library enabled the comparative mapping of the sea bass genome using the three-spined stickleback Gasterosteus aculeatus genome as a reference. BAC-end sequences (102,690) were aligned to the stickleback genome. The number of mappable BACs was improved using a two-fold coverage WGS dataset of sea bass resulting in a comparative BAC-map covering 87% of stickleback chromosomes with 588 BAC-contigs. The minimum size of 83 contigs covering 50% of the reference was 1.2 Mbp; the largest BAC-contig comprised 8.86 Mbp. More than 22,000 BAC-clones aligned with both ends to the reference genome. Intra-chromosomal rearrangements between sea bass and stickleback were identified. Size distributions of mapped BACs were used to calculate that the genome of sea bass may be only 1.3 fold larger than the 460 Mbp stickleback genome. Conclusions The BAC map is used for sequencing single BACs or BAC-pools covering defined genomic entities by second generation sequencing technologies. Together with the WGS dataset it initiates a sea bass genome sequencing project. This will allow the quantification of polymorphisms through resequencing, which is important for selecting highly performing domesticated fish.

Abstract Background Although bivalves are among the most-studied marine organisms because of their ecological role and economic importance, very little information is available on the genome sequences of oyster species. This report documents three large-scale cDNA sequencing projects for the Pacific oyster Crassostrea gigas initiated to provide a large number of expressed sequence tags that were subsequently compiled in a publicly accessible database. This resource allowed for the identification of a large number of transcripts and provides valuable information for ongoing investigations of tissue-specific and stimulus-dependant gene expression patterns. These data are crucial for constructing comprehensive DNA microarrays, identifying single nucleotide polymorphisms and microsatellites in coding regions, and for identifying genes when the entire genome sequence of C. gigas becomes available. Description In the present paper, we report the production of 40,845 high-quality ESTs that identify 29,745 unique transcribed sequences consisting of 7,940 contigs and 21,805 singletons. All of these new sequences, together with existing public sequence data, have been compiled into a publicly-available Website http://public-contigbrowser.sigenae.org:9090/Crassostrea_gigas/index.html. Approximately 43% of the unique ESTs had significant matches against the SwissProt database and 27% were annotated using Gene Ontology terms. In addition, we identified a total of 208 in silico microsatellites from the ESTs, with 173 having sufficient flanking sequence for primer design. We also identified a total of 7,530 putative in silico, single-nucleotide polymorphisms using existing and newly-generated EST resources for the Pacific oyster. Conclusion A publicly-available database has been populated with 29,745 unique sequences for the Pacific oyster Crassostrea gigas. The database provides many tools to search cleaned and assembled ESTs. The user may input and submit several filters, such as protein or nucleotide hits, to select and download relevant elements. This database constitutes one of the most developed genomic resources accessible among Lophotrochozoans, an orphan clade of bilateral animals. These data will accelerate the development of both genomics and genetics in a commercially-important species with the highest annual, commercial production of any aquatic organism.

Abstract Background The European seabass (Dicentrarchus labrax), one of the most extensively cultured species in European aquaculture productions, is, along with the gilthead sea bream (Sparus aurata), a prospective model species for the Perciformes which includes several other commercially important species. Massive mortalities may be caused by bacterial or viral infections in intensive aquaculture production. Revealing transcripts involved in immune response and studying their relative expression enhances the understanding of the immune response mechanism and consequently also the creation of vaccines. The analysis of expressed sequence tags (EST) is an efficient and easy approach for gene discovery, comparative genomics and for examining gene expression in specific tissues in a qualitative and quantitative way. Results Here we describe the construction, analysis and comparison of a total of ten cDNA libraries, six from different tissues infected with V. anguillarum (liver, spleen, head kidney, gill, peritoneal exudates and intestine) and four cDNA libraries from different tissues infected with Nodavirus (liver, spleen, head kidney and brain). In total 9605 sequences representing 3075 (32%) unique sequences (set of sequences obtained after clustering) were obtained and analysed. Among the sequences several immune-related proteins were identified for the first time in the order of Perciformes as well as in Teleostei. Conclusion The present study provides new information to the Gene Index of seabass. It gives a unigene set that will make a significant contribution to functional genomic studies and to studies of differential gene expression in relation to the immune system. In addition some of the potentially interesting genes identified by in silico analysis and confirmed by real-time PCR are putative biomarkers for bacterial and viral infections in fish.

Abstract Background Microscopic monogonont rotifers, including the euryhaline species Brachionus plicatilis, are typically found in water bodies where environmental factors restrict population growth to short periods lasting days or months. The survival of the population is ensured via the production of resting eggs that show a remarkable tolerance to unfavorable conditions and remain viable for decades. The aim of this study was to generate Expressed Sequence Tags (ESTs) for molecular characterisation of processes associated with the formation of resting eggs, their survival during dormancy and hatching. Results Four normalized and four subtractive libraries were constructed to provide a resource for rotifer transcriptomics associated with resting-egg formation, storage and hatching. A total of 47,926 sequences were assembled into 18,000 putative transcripts and analyzed using both Blast and GO annotation. About 28–55% (depending on the library) of the clones produced significant matches against the Swissprot and Trembl databases. Genes known to be associated with desiccation tolerance during dormancy in other organisms were identified in the EST libraries. These included genes associated with antioxidant activity, low molecular weight heat shock proteins and Late Embryonic Abundant (LEA) proteins. Real-time PCR confirmed that LEA transcripts, small heat-shock proteins and some antioxidant genes were upregulated in resting eggs, therefore suggesting that desiccation tolerance is a characteristic feature of resting eggs even though they do not necessarily fully desiccate during dormancy. The role of trehalose in resting-egg formation and survival remains unclear since there was no significant difference between resting-egg producing females and amictic females in the expression of the tps-1 gene. In view of the absence of vitellogenin transcripts, matches to lipoprotein lipase proteins suggest that, similar to the situation in dipterans, these proteins may serve as the yolk proteins in rotifers. Conclusion The 47,926 ESTs expand significantly the current sequence resource of B. plicatilis. It describes, for the first time, genes putatively associated with resting eggs and will serve as a database for future global expression experiments, particularly for the further identification of dormancy related genes.

Abstract Background Aquaculture represents the most sustainable alternative of seafood supply to substitute for the declining marine fisheries, but severe production bottlenecks remain to be solved. The application of genomic technologies offers much promise to rapidly increase our knowledge on biological processes in farmed species and overcome such bottlenecks. Here we present an integrated platform for mRNA expression profiling in the gilthead sea bream (Sparus aurata), a marine teleost of great importance for aquaculture. Results A public data base was constructed, consisting of 19,734 unique clusters (3,563 contigs and 16,171 singletons). Functional annotation was obtained for 8,021 clusters. Over 4,000 sequences were also associated with a GO entry. Two 60mer probes were designed for each gene and in-situ synthesized on glass slides using Agilent SurePrint™ technology. Platform reproducibility and accuracy were assessed on two early stages of sea bream development (one-day and four days old larvae). Correlation between technical replicates was always > 0.99, with strong positive correlation between paired probes. A two class SAM test identified 1,050 differentially expressed genes between the two developmental stages. Functional analysis suggested that down-regulated transcripts (407) in older larvae are mostly essential/housekeeping genes, whereas tissue-specific genes are up-regulated in parallel with the formation of key organs (eye, digestive system). Cross-validation of microarray data was carried out using quantitative qRT-PCR on 11 target genes, selected to reflect the whole range of fold-change and both up-regulated and down-regulated genes. A statistically significant positive correlation was obtained comparing expression levels for each target gene across all biological replicates. Good concordance between qRT-PCR and microarray data was observed between 2- and 7-fold change, while fold-change compression in the microarray was present for differences greater than 10-fold in the qRT-PCR. Conclusion A highly reliable oligo-microarray platform was developed and validated for the gilthead sea bream despite the presently limited knowledge of the species transcriptome. Because of the flexible design this array will be able to accommodate additional probes as soon as novel unique transcripts are available.

Abstract Background The goldfish (Carassius auratus) uses steroids and prostaglandins as pheromone cues at different stages of the reproductive cycle to facilitate spawning synchronization. Steroid progestin pheromone binding has been detected in goldfish olfactory membranes but the receptors responsible for this specific binding remain unknown. In order to shed some light on the olfactory epithelium transcriptome and search for possible receptor candidates a large set of EST from this tissue were analysed and compared to and combined with a similar zebrafish (Danio rerio) resource. Results We generated 4,797 high quality sequences from a normalized cDNA library of the goldfish olfactory epithelium, which were clustered in 3,879 unique sequences, grouped in 668 contigs and 3,211 singletons. BLASTX searches produced 3,243 significant (E-value < e-10) hits and Gene Ontology (GO) analysis annotated a further 1,223 of these genes (37.7%). Comparative analysis with zebrafish olfactory epithelium ESTs revealed 1,088 identical unigenes. The transcriptome size of both species was estimated at about 16,400 unigenes, based on the proportion of genes identified involved in Glucose Metabolic Process. Of 124 G-protein coupled receptors identified in the olfactory epithelium of both species, 56 were olfactory receptors. Beta and gamma membrane progestin receptors were also isolated by subcloning of RT-PCR products from both species and an olfactory epithelium specific splice form identified. Conclusion The high similarity between the goldfish and zebrafish olfactory systems allowed the creation of a 'cyprinid' olfactory epithelium library estimated to represent circa 70% of the transcriptome. These results are an important resource for the identification of components of signalling pathways involved in olfaction as well as putative targets for pharmacological and histochemical studies. The possible function of the receptors identified in the olfactory system is described. Moreover, the role of olfactory epithelium specific isoforms of classical membrane progestin receptor genes as candidates for preovulatory pheromone sensing is discussed.

Abstract Background Phytoplasmas are insect-transmitted, uncultivable bacterial plant pathogens that cause diseases in hundreds of economically important plants. They represent a monophyletic group within the class Mollicutes (trivial name mycoplasmas) and are characterized by a small genome with a low GC content, and the lack of a firm cell wall. All mycoplasmas, including strains of 'Candidatus (Ca.) Phytoplasma asteris' and 'Ca. P. australiense', examined so far have circular chromosomes, as is the case for almost all walled bacteria. Results Our work has shown that 'Ca. Phytoplasma mali', the causative agent of apple proliferation disease, has a linear chromosome. Linear chromosomes were also identified in the closely related provisional species 'Ca. P. pyri' and 'Ca. P. prunorum'. The chromosome of 'Ca. P. mali' strain AT is 601,943 bp in size and has a GC content of 21.4%. The chromosome is further characterized by large terminal inverted repeats and covalently closed hairpin ends. Analysis of the protein-coding genes revealed that glycolysis, the major energy-yielding pathway supposed for 'Ca. P. asteris', is incomplete in 'Ca. P. mali'. Due to the apparent lack of other metabolic pathways present in mycoplasmas, it is proposed that maltose and malate are utilized as carbon and energy sources. However, complete ATP-yielding pathways were not identified. 'Ca. P. mali' also differs from 'Ca. P. asteris' by a smaller genome, a lower GC content, a lower number of paralogous genes, fewer insertions of potential mobile DNA elements, and a strongly reduced number of ABC transporters for amino acids. In contrast, 'Ca. P. mali' has an extended set of genes for homologous recombination, excision repair and SOS response than 'Ca. P. asteris'. Conclusion The small linear chromosome with large terminal inverted repeats and covalently closed hairpin ends, the extremely low GC content and the limited metabolic capabilities reflect unique features of 'Ca. P. mali', not only within phytoplasmas, but all mycoplasmas. It is expected that the genome information obtained here will contribute to a better understanding of the reduced metabolism of phytoplasmas, their fastidious nutrition requirements that prevented axenic cultivation, and the mechanisms involved in pathogenicity.

Abstract Background Triosephosphate isomerase (TPI) is a central and conserved glycolytic enzyme. In humans, TPI is encoded by a single gene on 12p13, and associated with a rare genetic disorder, TPI deficiency. Reduced TPI activity can increase specific oxidant resistances of model organisms and TPI null-alleles have been hypothesized to promote a heterozygote advantage in man. However, comprehensive genetic information about the TPI1 locus is still lacking. Results Here, we sequenced the TPI1 locus in a sample of 357 German long-lived individuals (LLI) aged 95 to 110 years. We identified 17 different polymorphisms, of which 15 were rare and previously unknown. The two remaining SNPs occurred at much higher frequency and were tested for association with the longevity phenotype in larger samples of LLI (n = 1422) and younger controls (n = 967). Neither of the two markers showed a statistically significant difference in allele or genotype frequency between LLI and control subjects. Conclusion This study marks the TPI1 locus as extraordinarily conserved, even when analyzing intronic and non-coding regions of the gene. None of the identified sequence variations affected the amino acid composition of the TPI protein and hence, are unlikely to impact the catalytic activity of the enzyme. Thus, TPI variants occur less frequent than expected and inactive alleles are not enriched in German centenarians.

Abstract Background Ice, snow and temperatures of -14°C are conditions which most animals would find difficult, if not impossible, to survive in. However this exactly describes the Arctic winter, and the Arctic springtail Onychiurus arcticus regularly survives these extreme conditions and re-emerges in the spring. It is able to do this by reducing the amount of water in its body to almost zero: a process that is called "protective dehydration". The aim of this project was to generate clones and sequence data in the form of ESTs to provide a platform for the future molecular characterisation of the processes involved in protective dehydration. Results Five normalised libraries were produced from both desiccating and rehydrating populations of O. arcticus from stages that had previously been defined as potentially informative for molecular analyses. A total of 16,379 EST clones were generated and analysed using Blast and GO annotation. 40% of the clones produced significant matches against the Swissprot and trembl databases and these were further analysed using GO annotation. Extraction and analysis of GO annotations proved an extremely effective method for identifying generic processes associated with biochemical pathways, proving more efficient than solely analysing Blast data output. A number of genes were identified, which have previously been shown to be involved in water transport and desiccation such as members of the aquaporin family. Identification of these clones in specific libraries associated with desiccation validates the computational analysis by library rather than producing a global overview of all libraries combined. Conclusion This paper describes for the first time EST data from the arctic springtail (O. arcticus). This significantly enhances the number of Collembolan ESTs in the public databases, providing useful comparative data within this phylum. The use of GO annotation for analysis has facilitated the identification of a wide variety of ESTs associated with a number of different biochemical pathways involved in the dehydration and recovery process in O. arcticus.

Abstract Background The sheep is an important model animal for testing novel fracture treatments and other medical applications. Despite these medical uses and the well known economic and cultural importance of the sheep, relatively little research has been performed into sheep genetics, and DNA sequences are available for only a small number of sheep genes. Results In this work we have sequenced over 47 thousand expressed sequence tags (ESTs) from libraries developed from healing bone in a sheep model of fracture healing. These ESTs were clustered with the previously available 10 thousand sheep ESTs to a total of 19087 contigs with an average length of 603 nucleotides. We used the newly identified sequences to develop RT-PCR assays for 78 sheep genes and measured differential expression during the course of fracture healing between days 7 and 42 postfracture. All genes showed significant shifts at one or more time points. 23 of the genes were differentially expressed between postfracture days 7 and 10, which could reflect an important role for these genes for the initiation of osteogenesis. Conclusion The sequences we have identified in this work are a valuable resource for future studies on musculoskeletal healing and regeneration using sheep and represent an important head-start for genomic sequencing projects for Ovis aries, with partial or complete sequences being made available for over 5,800 previously unsequenced sheep genes.

In addition to the BAC-based reference sequence of the accession Columbia-0 from the year 2000, several short read assemblies of THE plant model organism Arabidopsis thaliana were published during the last years. Also, a SMRT-based assembly of Landsberg erecta has been generated that identified translocation and inversion polymorphisms between two genotypes of the species. Here we provide a chromosome-arm level assembly of the A. thaliana accession Niederzenz-1 (AthNd-1_v2c) based on SMRT sequencing data. The best assembly comprises 69 nucleome sequences and displays a contig length of up to 16 Mbp. Compared to an earlier Illumina short read-based NGS assembly (AthNd-1_v1), a 75 fold increase in contiguity was observed for AthNd-1_v2c. To assign contig locations independent from the Col-0 gold standard reference sequence, we used genetic anchoring to generate a de novo assembly. In addition, we assembled the chondrome and plastome sequences. Detailed analyses of AthNd-1_v2c allowed reliable identification of large genomic rearrangements between A. thaliana accessions contributing to differences in the gene sets that distinguish the genotypes. One of the differences detected identified a gene that is lacking from the Col-0 gold standard sequence. This de novo assembly extends the known proportion of the A. thaliana pan-genome. [ABSTRACT FROM AUTHOR]

Community-acquired (CA) Staphylococcus aureus cause various diseases even in healthy individuals. Enhanced virulence of CA-strains is partly attributed to increased production of toxins such as phenol-soluble modulins (PSM). The pathogen is internalized efficiently by mammalian host cells and intracellular S. aureus has recently been shown to contribute to disease. Upon internalization, cytotoxic S. aureus strains can disrupt phagosomal membranes and kill host cells in a PSM-dependent manner. However, PSM are not sufficient for these processes. Here we screened for factors required for intracellular S. aureus virulence. We infected escape reporter host cells with strains from an established transposon mutant library and detected phagosomal escape rates using automated microscopy. We thereby, among other factors, identified a non-ribosomal peptide synthetase (NRPS) to be required for efficient phagosomal escape and intracellular survival of S. aureus as well as induction of host cell death. By genetic complementation as well as supplementation with the synthetic NRPS product, the cyclic dipeptide phevalin, wild-type phenotypes were restored. We further demonstrate that the NRPS is contributing to virulence in a mouse pneumonia model. Together, our data illustrate a hitherto unrecognized function of the S. aureus NRPS and its dipeptide product during S. aureus infection. [ABSTRACT FROM AUTHOR]

The utility of genome assemblies does not only rely on the quality of the assembled genome sequence, but also on the quality of the gene annotations. The Pacific Biosciences Iso-Seq technology is a powerful support for accurate eukaryotic gene model annotation as it allows for direct readout of full-length cDNA sequences without the need for noisy short read-based transcript assembly. We propose the implementation of the TeloPrime Full Length cDNA Amplification kit to the Pacific Biosciences Iso-Seq technology in order to enrich for genuine full-length transcripts in the cDNA libraries. We provide evidence that TeloPrime outperforms the commonly used SMARTer PCR cDNA Synthesis Kit in identifying transcription start and end sites in Arabidopsis thaliana. Furthermore, we show that TeloPrime-based Pacific Biosciences Iso-Seq can be successfully applied to the polyploid genome of bread wheat (Triticum aestivum) not only to efficiently annotate gene models, but also to identify novel transcription sites, gene homeologs, splicing isoforms and previously unidentified gene loci. [ABSTRACT FROM AUTHOR]

The identification of body fluids is an essential tool for clarifying the course of events at a criminal site. The analytical problem is the fact that the biological material has been very often exposed to detrimental exogenous influences. Thereby, the molecular substrates used for the identification of the traces may become degraded. So far, most protocols utilize cell specific proteins or RNAs. Instead of measuring these more sensitive compounds this paper describes the application of the differential DNA-methylation. As a result of two genome wide screenings with the Illumina HumanMethylation BeadChips 27 and 450k we identified 150 candidate loci revealing differential methylation with regard to the body fluids venous blood, menstrual blood, vaginal fluid, saliva and sperm. Among them we selected 9 loci as the most promising markers. For the final determination of the methylation degree we applied the SNuPE-method. Because the degree of methylation might be modified by various endogenous and exogenous factors, we tested each marker with approximately 100 samples of each target fluid in a validation study. The stability of the detection procedure is proved in various simulated forensic surroundings according to standardized conditions. We studied the potential influence of 12 relatively common tumors on the methylation of the 9 markers. For this purpose the target fluids of 34 patients have been analysed. Only the cervix carcinoma might have an remarkable effect because impairing the signal of both vaginal markers. Using the Illumina MiSeq device we tested the potential influence of cis acting sequence variants on the methylation degree of the 9 markers in the specific body fluid DNA of 50 individuals. For 4 marker loci we observed such an influence either by sole SNPs or haplotypes. The identification of each target fluid is possible in arbitrary mixtures with the remaining four body fluids. The sensitivity of the individual body fluid tests is in the same range as for the forensic STR-analysis. It is the first forensic body fluid protocol which considers the exogenic and endogenic parameters potentially interfering with the true results. [ABSTRACT FROM AUTHOR]

The Epithelial Cell Adhesion Molecule (EpCAM) is overexpressed in many cancers including ovarian cancer and EpCAM overexpression correlates with decreased survival of patients. It was the aim of this study to achieve a targeted methylation of the EpCAM promoter and silence EpCAM gene expression using an engineered zinc finger protein that specifically binds the EpCAM promoter fused to the catalytic domain of the Dnmt3a DNA methyltransferase. We show that transient transfection of this construct increased the methylation of the EpCAM promoter in SKOV3 cells from 4–8% in untreated cells to 30%. Up to 48% methylation was observed in stable cell lines which express the chimeric methyltransferase. Control experiments confirmed that the methylation was dependent on the fusion of the Zinc finger and the methyltransferase domains and specific for the target region. The stable cell lines with methylated EpCAM promoter showed a 60–80% reduction of EpCAM expression as determined at mRNA and protein level and exhibited a significantly reduced cell proliferation. Our data indicate that targeted methylation of the EpCAM promoter could be an approach in the therapy of EpCAM overexpressing cancers. [ABSTRACT FROM AUTHOR]

Campylobacter jejuni is currently the leading cause of bacterial gastroenteritis in humans. Comparison of multiple Campylobacter strains revealed a high genetic and phenotypic diversity. However, little is known about differences in transcriptome organization, gene expression, and small RNA (sRNA) repertoires. Here we present the first comparative primary transcriptome analysis based on the differential RNA–seq (dRNA–seq) of four C. jejuni isolates. Our approach includes a novel, generic method for the automated annotation of transcriptional start sites (TSS), which allowed us to provide genome-wide promoter maps in the analyzed strains. These global TSS maps are refined through the integration of a SuperGenome approach that allows for a comparative TSS annotation by mapping RNA–seq data of multiple strains into a common coordinate system derived from a whole-genome alignment. Considering the steadily increasing amount of RNA–seq studies, our automated TSS annotation will not only facilitate transcriptome annotation for a wider range of pro- and eukaryotes but can also be adapted for the analysis among different growth or stress conditions. Our comparative dRNA–seq analysis revealed conservation of most TSS, but also single-nucleotide-polymorphisms (SNP) in promoter regions, which lead to strain-specific transcriptional output. Furthermore, we identified strain-specific sRNA repertoires that could contribute to differential gene regulation among strains. In addition, we identified a novel minimal CRISPR-system in Campylobacter of the type-II CRISPR subtype, which relies on the host factor RNase III and a trans-encoded sRNA for maturation of crRNAs. This minimal system of Campylobacter, which seems active in only some strains, employs a unique maturation pathway, since the crRNAs are transcribed from individual promoters in the upstream repeats and thereby minimize the requirements for the maturation machinery. Overall, our study provides new insights into strain-specific transcriptome organization and sRNAs, and reveals genes that could modulate phenotypic variation among strains despite high conservation at the DNA level. [ABSTRACT FROM AUTHOR]

X-linked intellectual disability (XLID), also known as X-linked mental retardation, is a highly genetically heterogeneous condition for which mutations in >90 different genes have been identified. In this study, we used a custom-made sequencing array based on the Affymetrix 50k platform for mutation screening in 17 known XLID genes in patients from 135 families and found eight single-nucleotide changes that were absent in controls. For four mutations affecting ATRX (p.1761M>T), PQBP1 (p.155R>X) and SLC6A8 (p.390P>L and p.477S>L), we provide evidence for a functional involvement of these changes in the aetiology of intellectual disability. [ABSTRACT FROM AUTHOR]

Filamentous cyanobacteria of the genus Lyngbya are important contributors to coral reef ecosystems, occasionally forming dominant cover and impacting the health of many other co-occurring organisms. Moreover, they are extraordinarily rich sources of bioactive secondary metabolites, with 35% of all reported cyanobacterial natural products deriving from this single pantropical genus. However, the true natural product potential and life strategies of Lyngbya strains are poorly understood because of phylogenetic ambiguity, lack of genomic information, and their close associations with heterotrophic bacteria and other cyanobacteria. To gauge the natural product potential of Lyngbya and gain insights into potential microbial interactions, we sequenced the genome of Lyngbya majuscula 3L, a Caribbean strain that produces the tubulin polymerization inhibitor curacin A and the molluscicide barbamide, using a combination of Sanger and 454 sequencing approaches. Whereas ∼293,000 nucleotides of the draft genome are putatively dedicated to secondary metabolism, this is far too few to encode a large suite of Lyngbya metabolites, suggesting Lyngbya metabolites are strain specific and may be useful in species delineation. Our analysis revealed a complex gene regulatory network, including a large number of sigma factors and other regulatory proteins, indicating an enhanced ability for environmental adaptation or microbial associations. Although Lyngbya species are reported to fix nitrogen, nitrogenase genes were not found in the genome or by PCR of genomic DNA. Subsequent growth experiments confirmed that L. majuscula 3L is unable to fix atmospheric nitrogen. These unanticipated life history characteristics challenge current views of the genus Lyngbya. [ABSTRACT FROM AUTHOR]

Primary ciliary dyskinesia (PCD) is a genetically heterogeneous autosomal recessive disorder characterized by recurrent infections of the respiratory tract associated with the abnormal function of motile cilia. Approximately half of individuals with PCD also have alterations in the left-right organization of their internal organ positioning, including situs inversus and situs ambiguous (Kartagener's syndrome). Here, we identify an uncharacterized coiled-coil domain containing a protein, CCDC40, essential for correct left-right patterning in mouse, zebrafish and human. In mouse and zebrafish, Ccdc40 is expressed in tissues that contain motile cilia, and mutations in Ccdc40 result in cilia with reduced ranges of motility. We further show that CCDC40 mutations in humans result in a variant of PCD characterized by misplacement of the central pair of microtubules and defective assembly of inner dynein arms and dynein regulatory complexes. CCDC40 localizes to motile cilia and the apical cytoplasm and is required for axonemal recruitment of CCDC39, disruption of which underlies a similar variant of PCD. [ABSTRACT FROM AUTHOR]

Background: Food supply from the ocean is constrained by the shortage of domesticated and selected fish. Development of genomic models of economically important fishes should assist with the removal of this bottleneck. European sea bass Dicentrarchus labrax L. (Moronidae, Perciformes, Teleostei) is one of the most important fishes in European marine aquaculture; growing genomic resources put it on its way to serve as an economic model. Results: End sequencing of a sea bass genomic BAC-library enabled the comparative mapping of the sea bass genome using the three-spined stickleback Gasterosteus aculeatus genome as a reference. BAC-end sequences (102,690) were aligned to the stickleback genome. The number of mappable BACs was improved using a two-fold coverage WGS dataset of sea bass resulting in a comparative BAC-map covering 87% of stickleback chromosomes with 588 BAC-contigs. The minimum size of 83 contigs covering 50% of the reference was 1.2 Mbp; the largest BAC-contig comprised 8.86 Mbp. More than 22,000 BAC-clones aligned with both ends to the reference genome. Intra-chromosomal rearrangements between sea bass and stickleback were identified. Size distributions of mapped BACs were used to calculate that the genome of sea bass may be only 1.3 fold larger than the 460 Mbp stickleback genome. Conclusions: The BAC map is used for sequencing single BACs or BAC-pools covering defined genomic entities by second generation sequencing technologies. Together with the WGS dataset it initiates a sea bass genome sequencing project. This will allow the quantification of polymorphisms through resequencing, which is important for selecting highly performing domesticated fish. [ABSTRACT FROM AUTHOR]

In the last decade, a new gene family encoding non-rearranging receptors, called novel immune-type receptors (NITRs), has been discovered in teleost fish. NITRs belong to the immunoglobulin superfamily and represent an extraordinarily divergent and rapidly evolving gene complex. Genomic analysis of a region spanning 270 kb led to the discovery of a NITR gene cluster in the European sea bass ( Dicentrarchus labrax). In total, 27 NITR genes and three putative pseudogenes, organised in a tandemly arrayed cluster, were identified. Sea bass NITR genes maintain the three major genomic organisations that appear to be essentially conserved among fish species along with new features presumably involving processes of intron loss, exon deletion and acquisition of new exons. Comparative and evolutionary analyses suggest that these receptors have evolved following a “birth-and-death” model of gene evolution in which duplication events together with lineage-specific gain and loss of individual members contributed to the rapid diversification of individual gene families. In this study, we demonstrate that species-specific gene expansions provide the raw material for diversifying, positive Darwinian selection favouring the evolution of a highly diverse array of molecules. [ABSTRACT FROM AUTHOR]

trans splicing of a spliced-leader RNA (SL RNA) to the 5′ ends of mRNAs has been shown to have a limited and sporadic distribution among eukaryotes. Within metazoans, only nematodes are known to process polycistronic pre-mRNAs, produced from operon units of transcription, into mature monocistronic mRNAs via an SL RNA trans-splicing mechanism. Here we demonstrate that a chordate with a highly compact genome, Oikopleura dioica, now joins Caenorhabditis elegans in coupling trans splicing with processing of polycistronic transcipts. We identified a single SL RNA which associates with Sm proteins and has a trimethyl guanosine cap structure reminiscent of spliceosomal snRNPs. The same SL RNA, estimated to be trans-spliced to at least 25% of O. dioica mRNAs, is used for the processing of both isolated or first cistrons and downstream cistrons in a polycistronic precursor. Remarkably, intercistronic regions in O. dioica are far more reduced than those in either nematodes or kinetoplastids, implying minimal cis-regulatory elements for coupling of 3′-end formation and trans splicing. [ABSTRACT FROM AUTHOR]

Nephronophthisis (NPHP), a group of autosomal recessive cystic kidney disorders, is the most common genetic cause of progressive renal failure in children and young adults. NPHP may be associated with Leber congenital amaurosis, tapeto-retinal degeneration, cerebellar ataxia, cone-shaped epiphyses, congenital oculomotor apraxia and hepatic fibrosis. Loci associated with an infantile type of NPHP on 9q22-q31 (NPHP2), juvenile types of NPHP on chromosomes 2q12-q13 (NPHP1) and 1p36 (NPHP4) and an adolescent type of NPHP on 3q21-q22 (NPHP3) have been mapped. NPHP1 and NPHP4 have been identified, and interaction of the respective encoded proteins nephrocystin and nephrocystin-4 has been shown. Here we report the identification of NPHP3, encoding a novel 1,330-amino acid protein that interacts with nephrocystin. We describe mutations in NPHP3 in families with isolated NPHP and in families with NPHP with associated hepatic fibrosis or tapeto-retinal degeneration. We show that the mouse ortholog Nphp3 is expressed in the node, kidney tubules, retina, respiratory epithelium, liver, biliary tract and neural tissues. In addition, we show that a homozygous missense mutation in Nphp3 is probably responsible for the polycystic kidney disease (pcy) mouse phenotype. Interventional studies in the pcy mouse have shown beneficial effects by modification of protein intake and administration of methylprednisolone, suggesting therapeutic strategies for treating individuals with NPHP3. [ABSTRACT FROM AUTHOR]

The sequence-based map of a part of the rhesus macaque major histocompatibility complex (MHC) extended class II region is presented. The sequenced region encompasses 67,401 bp and contains the SACM2L, RING1, FABGL and KE4 genes, as well as the HTATSF1-like and ZNF-like pseudogenes. Similar to human, but different from rat and mouse, no class I genes are found in the SACM2L-RING1 interval. The rhesus macaque extended MHC class II region shows a high degree of conservation of exonic as well as intronic and intergenic sequences compared with the respective human region. It is concluded that this particular genomic organization of the extended class II region—i.e., the absence of class I genes and the presence of the HTATSF1-like and ZNF-like pseudogenes—can be traced back to a common ancestor of humans and rhesus macaques about 23 million years ago. [ABSTRACT FROM AUTHOR]

Ras-related C3 toxin substrate 1 (Rac1) is a small Rho-GTPase with important functions in fundamental cellular processes such as cytoskeleton rearrangements, signal transduction, cell cycle progression and malignant transformation. Using Rac1 primer, we identified a 5.5-kb DNA sequence on chromosome 4 (Chr. 4) in the human genome, containing the intronless protein coding sequence of Rac1. Sequence analysis revealed features of a processed pseudogene, which we named ψ1Rac1, that could be detected by Southern blot and polymerase chain reaction (PCR) on genomic DNA. A ψ1Rac1 pseudogene transcript was not detected by reverse transcription-polymerase chain reaction (RT-PCR), nor had the ψ1Rac1 promoter any transcriptional activity. In addition, three other intronless pseudogenes of Rac1 on chromosomes 4, 13 and X were identified (ψ1Rac1–ψ4Rac1) sharing an 86–96% sequence similarity with Rac1. Neither RT-PCR with pseudogene specific restriction enzymes, nor the sequencing of 130 cDNA clones from benign and malignant breast tissue and cell lines, detected the transcription of any of the Rac1 pseudogenes (ψ2Rac1–ψ4Rac1). Existence of Rac1 pseudogenes should be taken into consideration when analyzing genomic alterations of the human Rac1 gene. [Copyright &y& Elsevier]

We have determined the sequence of a 4-Mb interval on rat chromosome 20p12 that encompasses the rat major histocompatibility complex (MHC). This is the first report of a finished sequence for a segment of the rat genome and constitutes one of the largest contiguous sequences thus far for rodent genomes in general. The rat MHC is, next to the human MHC, the second mammalian MHC sequenced to completion. Our analysis has resulted in the identification of at least 220 genes located within the sequenced interval. Although gene content and order are well conserved in the class II and class III gene intervals as well as the framework gene regions, profound rat-specific features were encountered within the class I gene regions, in comparison to human and mouse. Class I region-associated differences were found both at the structural level, the number, and organization of class I genes and gene families, and, in a more global context, in the way that evolution worked to shape the present-day rat MHC. [ABSTRACT FROM AUTHOR]

Homeodomain transcription factors are involved in many developmental processes and have been intensely studied in a few model organisms, such as mouse, Drosophila and Caenorhabditis elegans. Homeobox genes fall into 10 classes (ANTP, PRD, POU, LIM, TALE, SIX, Cut, ZFH, HNF1, Prox) and 89 different families/groups, all of which are present in vertebrates. Additional groups may be uncovered by further genome annotation, particularly of complex vertebrate genomes. Eight of these groups have been found only in vertebrates, but not in the genome of the tunicate Ciona intestinalis. The other 81 groups of homeobox gene that have been detected in vertebrates so far probably appeared during the early evolution of bilaterians or earlier, as they are also present outside the chordates. How the homeobox genes evolved during and after the main radiation of the bilaterians remains poorly understood, as only a few animal genomes have been sequenced completely. However, drastic changes have occurred at least in the lineage of C. elegans , such as loss of several Hox genes and Hox cluster fragmentation . Here we report considerable alterations of the homeobox gene complement in the tunicate lineage. [Copyright &y& Elsevier]