Incognito rRNA and rDNA in databases and libraries

Non-canonical ribosomal DNA segments in the human genome, and nucleoli functioning

Ribosomal DNA (rDNA) in the human genome is represented by tandem repeats of 43 kb nucleotide sequences that form nucleoli organizers (NORs) on each of five pairs of acrocentric chromosomes. RDNA-similar segments of different lengths are also present on (NOR)(-) chromosomes. Many of these segments contain nucleotide substitutions, supplementary microsatellite clusters, and extended deletions. Recently, it was shown that, in addition to ribosome biogenesis, nucleoli exhibit additional functions, such as cell-cycle regulation and response to stresses. In particular, several stress-inducible loci located in the ribosomal intergenic spacer (rIGS) produce stimuli-specific noncoding nucleolus RNAs. By mapping the 5'/3' ends of the rIGS segments scattered throughout (NOR)(-) chromosomes, we discovered that the bonds in the rIGS that were most often susceptible to disruption in the rIGS were adjacent to, or overlapped with stimuli-specific inducible loci. This suggests the interconnection of the two phenomena - nucleoli functioning and the scattering of rDNA-like sequences on (NOR)(-) chromosomes.

Genomic characterization of the mouse ribosomal DNA locus

The transcription of rRNA is critical to all living cells and is tightly controlled at the level of chromatin structure. Although the widespread adoption of genomic technologies including chromatin immunoprecipitation with massively parallel short-read sequencing (ChIP-seq) has allowed for the interrogation of chromatin structure on a genome-wide scale, until recently rDNA has not been analyzed by this technique. We extended genomic analysis of rDNA to mouse (Mus musculus), in which rDNA is similar in structure but highly divergent in sequence compared with human rDNA. Comparison of rDNA histone marks between mouse embryonic stem cells (mESCs) and more differentiated mouse cell types revealed differences between pluripotent and differentiated states. We also observed substantial divergence in rDNA histone modification patterns between mESCs and human embryonic stem cells (hESCs). Surprisingly, we found that the pluripotency factor OCT4 was bound to rDNA in similar patterns in mESCs and hESCs. Extending this analysis, we found that an additional 17 pluripotency-associated factors were bound to rDNA in mESCs, suggesting novel modes of rDNA regulation in pluripotent cells. Taken together, our results provide a detailed view of rDNA chromatin structure in an important model system and enable high-resolution comparison of rDNA regulation between mouse and human.

Expression of stress-related genes in zebrawood (Astronium fraxinifolium, Anacardiaceae) seedlings following germination in microgravity

Seeds of a tropical tree species from Brazil, Astronium fraxinifolium, or zebrawood, were germinated, for the first time in microgravity, aboard the International Space Station for nine days. Following three days of subsequent growth under normal terrestrial gravitational conditions, greater root length and numbers of secondary roots was observed in the microgravity-treated seedlings compared to terrestrially germinated controls. Suppression subtractive hybridization of cDNA and EST analysis were used to detect differential gene expression in the microgravity-treated seedlings in comparison to those initially grown in normal gravity (forward subtraction). Despite their return to, and growth in normal gravity, the subtracted library derived from microgravity-treated seedlings was enriched in known microgravity stress-related ESTs, corresponding to large and small heat shock proteins, 14-3-3-like protein, polyubiquitin, and proteins involved in glutathione metabolism. In contrast, the reverse-subtracted library contained a comparatively greater variety of general metabolism-related ESTs, but was also enriched for peroxidase, possibly indicating the suppression of this protein in the microgravity-treated seedlings. Following continued growth for 30 days, higher concentrations of total chlorophyll were detected in the microgravity-exposed seedlings.

Gut pathology and responses to the microsporidium Nosema ceranae in the honey bee Apis mellifera

The microsporidium Nosema ceranae is a newly prevalent parasite of the European honey bee (Apis mellifera). Although this parasite is presently spreading across the world into its novel host, the mechanisms by it which affects the bees and how bees respond are not well understood. We therefore performed an extensive characterization of the parasite effects at the molecular level by using genetic and biochemical tools. The transcriptome modifications at the midgut level were characterized seven days post-infection with tiling microarrays. Then we tested the bee midgut response to infection by measuring activity of antioxidant and detoxification enzymes (superoxide dismutases, glutathione peroxidases, glutathione reductase, and glutathione-S-transferase). At the gene-expression level, the bee midgut responded to N. ceranae infection by an increase in oxidative stress concurrent with the generation of antioxidant enzymes, defense and protective response specifically observed in the gut of mammals and insects. However, at the enzymatic level, the protective response was not confirmed, with only glutathione-S-transferase exhibiting a higher activity in infected bees. The oxidative stress was associated with a higher transcription of sugar transporter in the gut. Finally, a dramatic effect of the microsporidia infection was the inhibition of genes involved in the homeostasis and renewal of intestinal tissues (Wnt signaling pathway), a phenomenon that was confirmed at the histological level. This tissue degeneration and prevention of gut epithelium renewal may explain early bee death. In conclusion, our integrated approach not only gives new insights into the pathological effects of N. ceranae and the bee gut response, but also demonstrate that the honey bee gut is an interesting model system for studying host defense responses.

Identification of a tubulin-α gene specifically expressed in testis and adductor muscle during stable reference gene selection in the hermaphrodite gonad of the lion's paw scallop Nodipecten subnodosus

For non-model species, as many used for aquaculture, with minimal or no genomic information, relative quantification of gene expression studies requires preliminary research including the isolation of potential reference genes and the identification of those stably expressed under the biological conditions of interest. Here we report on the isolation of five partial gene sequences from gonad tissue cDNA in the functional hermaphrodite scallop Nodipecten subnodosus to be evaluated as reference genes: 18S-rRNA, riboprotein l8 (rp-l8), actin-β (act-β), elongation factor 1α (ef-1α) and alpha-tubulin-α (tub-α). We found that 18S-rRNA was stably expressed independently of the priming method used to reverse transcribe RNA to cDNA, oligo-dT or random hexamer. Stability analysis for the five putative reference genes with geNorm and NormFinder indicated that 18S together with rp-l8 were the most stable genes for normalization of gene expression during gonad development in both, male and female sexual regions of the hermaphrodite N. subnodosus. The least stable gene was tub-α, showing a biased expression profile between sexual regions of the gonad, therefore this gene was analyzed thereafter as a target gene together with vitellogenin (vit) and a DEAD-box RNA helicase (dbx) gene. Relative expression, estimated by normalization with the combination of 18S and rp-l8 as reference genes, indicated that as gonad development advanced two of the target genes were up-regulated, tub-α in the male region and vit in the female region. Whereas an increased expression was expected during development for vit for its known role in vitellogenesis, the increased expression of tub-α in the male sexual region was unexpected, and pointed toward this gene being a testis-specific α-tubulin isotype. Further analyses of gene expression among tissues indicated that tub-α is specifically and highly expressed in the male gonad, although expression in adductor muscle was also observed at significantly lower levels. The existence of testis specific α- and β-tubulins has been previously reported in other taxa, relating their function to sperm axoneme formation. Tissue-specific tubulin genes, particularly their promoters, have recently found an application as native promoters for transgene tissue-specific expression in research and reproductive control of insect plagues. The third target gene, a putative member of the DEAD-box RNA helicase family (dbx), showed no changes in expression during gonad development or between sexual regions, therefore it was chosen to discuss the different statistical inferences resulting from the arbitrary use of 'randomly chosen' reference genes when normalizing gene expression.

Use of methylation filtration and C(0)t fractionation for analysis of genome composition and comparative genomics in bread wheat

We investigated the compositional and structural differences in sequences derived from different fractions of wheat genomic DNA obtained using methylation filtration and C(0)t fractionation. Comparative analysis of these sequences revealed large compositional and structural variations in terms of GC content, different structural elements including repeat sequences (e.g., transposable elements and simple sequence repeats), protein coding genes, and non-coding RNA genes. A correlation between methylation status [determined on the basis of selective inclusion/exclusion in methylation-filtered (MF) library] of different repeat elements and expression level was observed. The expression levels were determined by comparing MF sequences with expressed sequence tags (ESTs) available in the public domain. Only a limited overlap among MF, high C(0)t (HC), and ESTs was observed, suggesting that these sequences may largely either represent the low-copy non-transcribed sequences or include genes with low expression levels. Thus, these results indicated a need to study MF and HC sequences along with ESTs to fully appreciate complexity of wheat gene space.

Changes in transcript abundance relating to colony collapse disorder in honey bees (Apis mellifera)

Colony collapse disorder (CCD) is a mysterious disappearance of honey bees that has beset beekeepers in the United States since late 2006. Pathogens and other environmental stresses, including pesticides, have been linked to CCD, but a causal relationship has not yet been demonstrated. Because the gut acts as a primary interface between the honey bee and its environment as a site of entry for pathogens and toxins, we used whole-genome microarrays to compare gene expression between guts of bees from CCD colonies originating on both the east and west coasts of the United States and guts of bees from healthy colonies sampled before the emergence of CCD. Considerable variation in gene expression was associated with the geographical origin of bees, but a consensus list of 65 transcripts was identified as potential markers for CCD status. Overall, elevated expression of pesticide response genes was not observed. Genes involved in immune response showed no clear trend in expression pattern despite the increased prevalence of viruses and other pathogens in CCD colonies. Microarray analysis revealed unusual ribosomal RNA fragments that were conspicuously more abundant in the guts of CCD bees. The presence of these fragments may be a possible consequence of picorna-like viral infection, including deformed wing virus and Israeli acute paralysis virus, and may be related to arrested translation. Ribosomal fragment abundance and presence of multiple viruses may prove to be useful diagnostic markers for colonies afflicted with CCD.

A pig multi-tissue normalised cDNA library: large-scale sequencing, cluster analysis and 9K micro-array resource generation

Background

Domestic animal breeding and product quality improvement require the control of reproduction, nutrition, health and welfare in these animals. It is thus necessary to improve our knowledge of the major physiological functions and their interactions. This would be greatly enhanced by the availability of expressed gene sequences in the databases and by cDNA arrays allowing the transcriptome analysis of any function.

The objective within the AGENAE French program was to initiate a high-throughput cDNA sequencing program of a 38-tissue normalised library and generate a diverse microarray for transcriptome analysis in pig species.

Results

We constructed a multi-tissue cDNA library, which was normalised and subtracted to reduce the redundancy of the clones. Expressed Sequence Tags were produced and 24449 high-quality sequences were released in EMBL database. The assembly of all the public ESTs (available through SIGENAE website) resulted in 40786 contigs and 54653 singletons. At least one Agenae sequence is present in 11969 contigs (12.5%) and in 9291 of the deeper-than-one-contigs (22.8%). Sequence analysis showed that both normalisation and subtraction processes were successful and that the initial tissue complexity was maintained in the final libraries. A 9K nylon cDNA microarray was produced and is available through CRB-GADIE. It will allow high sensitivity transcriptome analyses in pigs.

Conclusion

In the present work, a pig multi-tissue cDNA library was constructed and a 9K cDNA microarray designed. It contributes to the Expressed Sequence Tags pig data, and offers a valuable tool for transcriptome analysis.

Evaluating and improving cDNA sequence quality with cQC

SUMMARY

Errors are prevalent in cDNA sequences but the extent to which sequence collections differ in frequencies and types of errors has not been investigated systematically. cDNA quality control, or cQC, was developed to evaluate the quality of cDNA sequence collections and to revise those sequences that differ from a higher quality genomic sequence. After removing rRNA, vector, bacterial insertion sequence and chimeric cDNA contaminants, small-scale nucleotide discrepancies were found in 51% of cDNA sequences from one Arabidopsis cDNA collection, 89% from a second Arabidopsis collection and 75% from a rice collection. These errors created premature termination codons in 4 and 42% of cDNA sequences in the respective Arabidopsis collections and in 7% of the rice cDNA sequences.

Sequencing, transcript identification, and quantitative gene expression profiling in the breast cancer loss of heterozygosity region 16q24.3 reveal three potential tumor-suppressor genes

Loss of heterozygosity (LOH) of chromosome 16q24.3 is a common genetic alteration observed in invasive ductal and lobular breast carcinomas. We constructed a physical map and generated genomic DNA sequence data spanning 2.4 Mb in this region. Detailed in silico and in vitro analyses of the genomic sequence data enabled the identification of 104 genes. It was hypothesized that tumor-suppressor genes would exhibit marked mRNA expression variability in a panel of breast cancer cell lines as a result of downregulation due to mutation or hypermethylation. We examined the mRNA expression profiles of the genes identified at 16q24.3 in normal breast, a normal breast epithelial cell line, and several breast cancer cell lines exhibiting 16q24.3 LOH. Three of the genes, CYBA, Hs.7970, and CBFA2T3, exhibited variability ten times higher than the baseline. The possible role of these genes as tumor suppressors is discussed.

Gene2EST: a BLAST2 server for searching expressed sequence tag (EST) databases with eukaryotic gene-sized queries

Expressed sequence tags (ESTs) are randomly sequenced cDNA clones. Currently, nearly 3 million human and 2 million mouse ESTs provide valuable resources that enable researchers to investigate the products of gene expression. The EST databases have proven to be useful tools for detecting homologous genes, for exon mapping, revealing differential splicing, etc. With the increasing availability of large amounts of poorly characterised eukaryotic (notably human) genomic sequence, ESTs have now become a vital tool for gene identification, sometimes yielding the only unambiguous evidence for the existence of a gene expression product. However, BLAST-based Web servers available to the general user have not kept pace with these developments and do not provide appropriate tools for querying EST databases with large highly spliced genes, often spanning 50 000-100 000 bases or more. Here we describe Gene2EST (http://woody.embl-heidelberg.de/gene2est/), a server that brings together a set of tools enabling efficient retrieval of ESTs matching large DNA queries and their subsequent analysis. RepeatMasker is used to mask dispersed repetitive sequences (such as Alu elements) in the query, BLAST2 for searching EST databases and Artemis for graphical display of the findings. Gene2EST combines these components into a Web resource targeted at the researcher who wishes to study one or a few genes to a high level of detail.

Identification of differentially expressed genes in Con A-activated carp (Cyprinus carpio L.) leucocytes

A cDNA library was constructed from the message RNA (mRNA) obtained from Con A-induced head kidney (HK) leucocytes of carp (Cyprinus carpio L.). Differential screening of the cDNA was carried out by hybridization against the total cDNA probes from normal, Con A-uninduced HK leucocytes or Con A-induced HK leucocytes of carp. The differential expression patterns of certain cDNA clones were confirmed by Southern-blot and Northern-blot analysis. Single-pass of the sequencing analysis and homology search in Genbank (EMBL) revealed those differentially expressed cDNA clones encode for cytochrome c oxidase sub-unit II and III (COII and COIII), elongation factor-1 beta (EF-1 beta), bleomycin hydrolase (BH), heat shock cognate protein 70 (HSC70) and 16S ribosomal RNA (16S rRNA).

Beyond ribosomal DNA: on towards the telomere

We have sequenced and analyzed 8.3 kb of sequence adjacent and distal to the human ribosomal DNA (rDNA); this distal sequence connects to the rDNA cluster just 4 kb upstream of the first promoter and is shared among the acrocentric chromosomes and, at least in part, it is also present in other primates. The sequence differs in character from that of the rDNA intergenic spacer (IGS) in that it does not contain long stretches of either polypyrimidine or polypurine. However, just like the IGS, it contains numerous repetitive elements, including retroposed fragments of 28S rRNA and large pieces of the IGS. In addition, we show that the rDNA clusters are not interrupted by other sequences and do not recombine with this distal segment.