Preselection of shotgun clones by oligonucleotide fingerprinting: an efficient and high throughput strategy to reduce redundancy in large-scale sequencing projects

Microbiota in cancer: molecular mechanisms and therapeutic interventions

The diverse bacterial populations within the symbiotic microbiota play a pivotal role in both health and disease. Microbiota modulates critical aspects of tumor biology including cell proliferation, invasion, and metastasis. This regulation occurs through mechanisms like enhancing genomic damage, hindering gene repair, activating aberrant cell signaling pathways, influencing tumor cell metabolism, promoting revascularization, and remodeling the tumor immune microenvironment. These microbiota-mediated effects significantly impact overall survival and the recurrence of tumors after surgery by affecting the efficacy of chemoradiotherapy. Moreover, leveraging the microbiota for the development of biovectors, probiotics, prebiotics, and synbiotics, in addition to utilizing antibiotics, dietary adjustments, defensins, oncolytic virotherapy, and fecal microbiota transplantation, offers promising alternatives for cancer treatment. Nonetheless, due to the extensive and diverse nature of the microbiota, along with tumor heterogeneity, the molecular mechanisms underlying the role of microbiota in cancer remain a subject of intense debate. In this context, we refocus on various cancers, delving into the molecular signaling pathways associated with the microbiota and its derivatives, the reshaping of the tumor microenvironmental matrix, and the impact on tolerance to tumor treatments such as chemotherapy and radiotherapy. This exploration aims to shed light on novel perspectives and potential applications in the field.

Degenerate primer design: theoretical analysis and the HYDEN program

A polymerase chain reaction (PCR) primer sequence is called degenerate if some of its positions have several possible bases. The degeneracy of the primer is the number of unique sequence combinations it contains. We study the problem of designing a pair of primers with prescribed degeneracy that match a maximum number of given input sequences. Such problems occur, for example, when studying a family of genes that is known only in part or is known in a related species. We discuss the complexity of several versions of the problem and give approximation algorithms for one simplified variant. On the basis of these algorithms, we developed a program called HYDEN for designing highly degenerate primers for a set of genomic sequences. We describe HYDEN, and report on its success in several applications for identifying olfactory receptor genes in mammals.

Analyses of the vrl gene cluster in Desulfococcus multivorans: homologous to the virulence-associated locus of the ovine footrot pathogen Dichelobacter nodosus strain A198

Major parts of the virulence-associated vrl locus known from the gammaproteobacterium Dichelobacter nodosus, the causative agent of ovine footrot, were analyzed in the genome of the sulfate-reducing deltaproteobacterium Desulfococcus multivorans. In the genome of D. multivorans 13 of the 19 vrl genes described for D. nodosus are present and highly conserved with respect to gene sequence and order. The vrl locus and its flanking regions suggest a bacteriophage-mediated transfer into the genome of D. multivorans. Comparative analysis of the deduced Vrl proteins reveals a wide distribution of parts of the virulence-associated vrl locus in distantly related bacteria. Horizontal transfer is suggested as driving mechanism for the circulation of the vrl genes in bacteria. Except for the vrlBMN genes D. multivorans and Desulfovibrio desulfuricans G20 together contain all vrl genes displaying a high degree of similarity. For D. multivorans it could be shown that guanine plus cytosine (GC) content, GC skew, di-, tri- or tetranucleotide distribution did not differ between the vrl locus and its flanking sequences. This could be a hint that the vrl locus originated from a related organism or at least a genome with similar characteristics. The conspicuous high degree of conservation of the analyzed vrl genes may result from a recent transfer event or reflect a function of the vrl genes, which is still unknown and not necessarily disease associated. The latter is supported by the evidence for expression of the vrl genes in D. multivorans, which has not been described as pathogen or to be associated to any disease pattern before.

The degenerate primer design problem: theory and applications

A PCR primer sequence is called degenerate if some of its positions have several possible bases. The degeneracy of the primer is the number of unique sequence combinations it contains. We study the problem of designing a pair of primers with prescribed degeneracy that match a maximum number of given input sequences. Such problems occur when studying a family of genes that is known only in part, or is known in a related species. We prove that various simplified versions of the problem are hard, show the polynomiality of some restricted cases, and develop approximation algorithms for one variant. Based on these algorithms, we implemented a program called HYDEN for designing highly degenerate primers for a set of genomic sequences. We report on the success of the program in several applications, one of which is an experimental scheme for identifying all human olfactory receptor (OR) genes. In that project, HYDEN was used to design primers with degeneracies up to 10(10) that amplified with high specificity many novel genes of that family, tripling the number of OR genes known at the time.

A 6374 unigene set corresponding to low abundance transcripts expressed following fertilization in Solanum chacoense Bitt, and characterization of 30 receptor-like kinases

In order to characterize regulatory genes that are expressed in ovule tissues after fertilization we have undertaken an EST sequencing project in Solanum chacoense, a self-incompatible wild potato species. Two cDNA libraries made from ovule tissues covering embryo development from zygote to late torpedo-stage were constructed and plated at high density on nylon membranes. To decrease EST redundancy and enrich for transcripts corresponding to weakly expressed genes a self-probe subtraction method was used to select the colonies harboring the genes to be sequenced. 7741 good sequences were obtained and, from these, 6374 unigenes were isolated. Thus, the self-probe subtraction resulted in a strong enrichment in singletons, a decrease in the number of clones per contigs, and concomitantly, an enrichment in the total number of unigenes obtained (82%). To gain insights into signal transduction events occurring during embryo development all the receptor-like kinases (or protein receptor kinases) were analyzed by quantitative real-time RT-PCR. Interestingly, 28 out of the 30 RLK isolated were predominantly expressed in ovary tissues or young developing fruits, and 23 were transcriptionaly induced following fertilization. Thus, the self-probe subtraction did not select for genes weakly expressed in the target tissue while being highly expressed elsewhere in the plant. Of the receptor-like kinases (RLK) genes isolated, the leucine-rich repeat (LRR) family of RLK was by far the most represented with 25 members covering 11 LRR classes.

Establishment of T cell-specific and natural killer cell-specific unigene sets: towards high-throughput genomics of leukaemia

We report the establishment of highly non-redundant unigene sets consisting of cDNA clones derived from T lymphocytes and natural killer cells. Each set consists of 10 506 and 13 409 clones, respectively, arrayed on nylon membranes in duplicate. The sets provide an excellent tool for genome-wide gene expression analysis studies in immunology research.

The mouse homeobox gene Not is required for caudal notochord development and affected by the truncate mutation

The floating head (flh) gene in zebrafish encodes a homeodomain protein, which is essential for notochord formation along the entire body axis. flh orthologs, termed Not genes, have been isolated from chick and Xenopus, but no mammalian ortholog has yet been identified. Truncate (tc) is an autosomal recessive mutation in mouse that specifically disrupts the development of the caudal notochord. Here, we demonstrate that truncate arose by a mutation in the mouse Not gene. The truncate allele (Nottc) contains a point mutation in the homeobox of Not that changes a conserved Phenylalanine residue in helix 1 to a Cysteine (F20C), and significantly destabilizes the homeodomain. Reversion of F20C in one allele of homozygous tc embryonic stem (ES) cells is sufficient to restore normal notochord formation in completely ES cell-derived embryos. We have generated a targeted mutation of Not by replacing most of the Not coding sequence, including the homeobox with the eGFP gene. The phenotype of NoteGFP/eGFP, NoteGFP/tc, and Nottc/tc embryos is very similar but slightly more severe in NoteGFP/eGFP than in Nottc/tc embryos. This confirms allelism of truncate and Not, and indicates that tc is not a complete null allele. Not expression is abolished in Foxa2 and T mutant embryos, suggesting that Not acts downstream of both genes during notochord development. This is in contrast to zebrafish embryos, in which flh interacts with ntl (zebrafish T) in a regulatory loop and is essential for development of the entire notochord, and suggests that different genetic control circuits act in different vertebrate species during notochord formation.

Characterization of a spontaneous nonmagnetic mutant of Magnetospirillum gryphiswaldense reveals a large deletion comprising a putative magnetosome island

Frequent spontaneous loss of the magnetic phenotype was observed in stationary-phase cultures of the magnetotactic bacterium Magnetospirillum gryphiswaldense MSR-1. A nonmagnetic mutant, designated strain MSR-1B, was isolated and characterized. The mutant lacked any structures resembling magnetosome crystals as well as internal membrane vesicles. The growth of strain MSR-1B was impaired under all growth conditions tested, and the uptake and accumulation of iron were drastically reduced under iron-replete conditions. A large chromosomal deletion of approximately 80 kb was identified in strain MSR-1B, which comprised both the entire mamAB and mamDC clusters as well as further putative operons encoding a number of magnetosome-associated proteins. A bacterial artificial chromosome clone partially covering the deleted region was isolated from the genomic library of wild-type M. gryphiswaldense. Sequence analysis of this fragment revealed that all previously identified mam genes were closely linked with genes encoding other magnetosome-associated proteins within less than 35 kb. In addition, this region was remarkably rich in insertion elements and harbored a considerable number of unknown gene families which appeared to be specific for magnetotactic bacteria. Overall, these findings suggest the existence of a putative large magnetosome island in M. gryphiswaldense and other magnetotactic bacteria.

Honeybee (Apis mellifera L.) mrjp gene family: computational analysis of putative promoters and genomic structure of mrjp1, the gene coding for the most abundant protein of larval food

Mrjp1 gene belongs to the honeybee mrjp gene family encoding the major royal jelly proteins (MRJPs), secreted by nurse bees into the royal jelly. In this study, we have isolated the genomic clone containing the entire mrjp1 gene and determined its sequence. The mrjp1 gene sequence spans over 3038 bp and contains six exons separated by five introns. Seven mismatches between the mrjp1 gene sequence and two previously independently published cDNA sequences were found, but these differences do not lead to any change in the deduced amino acid sequence of MRJP1. With the aid of inverse polymerase chain reaction we obtained sequences flanking the 5' ends of other mrjp genes (mrjp2, mrjp3, mrjp4 and mrjp5). Putative promoters were predicted upstream of all mrjp genes (including mrjp1). The predicted promoters contain the TATA motif (TATATATT), highly conserved both in sequence and position. Ultraspiracle (USP) transcription factor (TF) binding sites in putative promoter regions and clusters of dead ringer TF binding sites upstream of these promoters were predicted computationally. We propose that USP, as a juvenile hormone (JH) binding TF, might possibly act as a mediator of mrjp expression in response to JH. Mrjp1's genomic locus is predicted to encode an antisense transcript, partially overlapping with five mrjp1 exons and entirely overlapping with the putative promoter and predicted transcriptional start point of mrjp1. This finding may shed light on the mechanisms of regulation of mrjps expression. Southern blot analysis of genomic DNA revealed that all so far known members of mrjp gene family (mrjp1, mrjp2, mrjp3, mrjp4 and mrjp5) are present as single-copy genes per haploid honeybee genome. Although MRJPs and the yellow protein of Drosophila melanogaster share a certain degree of similarity in aa sequence and although it has been shown that they share a common evolutionary origin, neither structural similarities in the gene organization, nor significant similarities between intron sequences of mrjp1 gene and fourteen yellow-like genes of D. melanogaster were found.

Construction of a 'unigene' cDNA clone set by oligonucleotide fingerprinting allows access to 25 000 potential sugar beet genes

Access to the complete gene inventory of an organism is crucial to understanding physiological processes like development, differentiation, pathogenesis, or adaptation to the environment. Transcripts from many active genes are present at low copy numbers. Therefore, procedures that rely on random EST sequencing or on normalisation and subtraction methods have to produce massively redundant data to get access to low-abundance genes. Here, we present an improved oligonucleotide fingerprinting (ofp) approach to the genome of sugar beet (Beta vulgaris), a plant for which practically no molecular information has been available. To identify distinct genes and to provide a representative 'unigene' cDNA set for sugar beet, 159 936 cDNA clones were processed utilizing large-scale, high-throughput data generation and analysis methods. Data analysis yielded 30 444 ofp clusters reflecting the number of different genes in the original cDNA sample. A sample of 10 961 cDNA clones, each representing a different cluster, were selected for sequencing. Standard sequence analysis confirmed that 89% of these EST sequences did represent different genes. These results indicate that the full set of 30 444 ofp clusters represent up to 25 000 genes. We conclude that the ofp analysis pipeline is an accurate and effective way to construct large representative 'unigene' sets for any plant of interest with no requirement for prior molecular sequence data.

DEFOG: a practical scheme for deciphering families of genes

We developed a novel efficient scheme, DEFOG (for "deciphering families of genes"), for determining sequences of numerous genes from a family of interest. The scheme provides a powerful means to obtain a gene family composition in species for which high-throughput genomic sequencing data are not available. DEFOG uses two key procedures. The first is a novel algorithm for designing highly degenerate primers based on a set of known genes from the family of interest. These primers are used in PCR reactions to amplify the members of the gene family. The second combines oligofingerprinting of the cloned PCR products with clustering of the clones based on their fingerprints. By selecting members from each cluster, a low-redundancy clone subset is chosen for sequencing. We applied the scheme to the human olfactory receptor (OR) genes. OR genes constitute the largest gene superfamily in the human genome, as well as in the genomes of other vertebrate species. DEFOG almost tripled the size of the initial repertoire of human ORs in a single experiment, and only 7% of the PCR clones had to be sequenced. Extremely high degeneracies, reaching over a billion combinations of distinct PCR primer pairs, proved to be very effective and yielded only 0.4% nonspecific products.

An oligonucleotide fingerprint normalized and expressed sequence tag characterized zebrafish cDNA library

The zebrafish is a powerful system for understanding the vertebrate genome, allowing the combination of genetic, molecular, and embryological analysis. Expressed sequence tags (ESTs) provide a rapid means of identifying an organism's genes for further analysis, but any EST project is limited by the availability of suitable libraries. Such cDNA libraries must be of high quality and provide a high rate of gene discovery. However, commonly used normalization and subtraction procedures tend to select for shorter, truncated, and internally primed inserts, seriously affecting library quality. An alternative procedure is to use oligonucleotide fingerprinting (OFP) to precluster clones before EST sequencing, thereby reducing the re-sequencing of common transcripts. Here, we describe the use of OFP to normalize and subtract 75,000 clones from two cDNA libraries, to a minimal set of 25,102 clones. We generated 25,788 ESTs (11,380 3' and 14,408 5') from over 16,000 of these clones. Clustering of 10,654 high-quality 3' ESTs from this set identified 7232 clusters (likely genes), corresponding to a 68% gene diversity rate, comparable to what has been reported for the best normalized human cDNA libraries, and indicating that the complete set of 25,102 clones contains as many as 17,000 genes. Yet, the library quality remains high. The complete set of 25,102 clones is available for researchers as glycerol stocks, filters sets, and as individual EST clones. These resources have been used for radiation hybrid, genetic, and physical mapping of the zebrafish genome, as well as positional cloning and candidate gene identification, molecular marker, and microarray development.

RAI1 is a novel polyglutamine encoding gene that is deleted in Smith-Magenis syndrome patients

The human chromosomal band 17p11.2 is a genetically unstable interval. It has been shown to be deleted in patients suffering from Smith-Magenis syndrome. Previous efforts of physical and transcriptional mapping in 17p11.2 and subsequent genomic sequencing of the candidate interval allowed the identification of new genes that might be responsible for the Smith-Magenis syndrome. In this report, one of these genes named RAI1, the human homologue of the mouse Rai1 gene, has been investigated for its contribution to the syndrome. Expression analysis on different human adult and fetal tissues has shown the existence of at least three splice variants. Moreover, the most interesting feature of the gene is the presence of a polymorphic CAG repeat coding for a polyglutamine stretch in the amino terminal domain of the protein.

Tissue gene expression analysis using arrayed normalized cDNA libraries

We have used oligonucleotide-fingerprinting data on 60,000 cDNA clones from two different mouse embryonic stages to establish a normalized cDNA clone set. The normalized set of 5,376 clones represents different clusters and therefore, in almost all cases, different genes. The inserts of the cDNA clones were amplified by PCR and spotted on glass slides. The resulting arrays were hybridized with mRNA probes prepared from six different adult mouse tissues. Expression profiles were analyzed by hierarchical clustering techniques. We have chosen radioactive detection because it combines robustness with sensitivity and allows the comparison of multiple normalized experiments. Sensitive detection combined with highly effective clustering algorithms allowed the identification of tissue-specific expression profiles and the detection of genes specifically expressed in the tissues investigated. The obtained results are publicly available (http://www.rzpd.de) and can be used by other researchers as a digital expression reference.

Large-scale clustering of cDNA-fingerprinting data

Clustering is one of the main mathematical challenges in large-scale gene expression analysis. We describe a clustering procedure based on a sequential k-means algorithm with additional refinements that is able to handle high-throughput data in the order of hundreds of thousands of data items measured on hundreds of variables. The practical motivation for our algorithm is oligonucleotide fingerprinting-a method for simultaneous determination of expression level for every active gene of a specific tissue-although the algorithm can be applied as well to other large-scale projects like EST clustering and qualitative clustering of DNA-chip data. As a pairwise similarity measure between two p-dimensional data points, x and y, we introduce mutual information that can be interpreted as the amount of information about x in y, and vice versa. We show that for our purposes this measure is superior to commonly used metric distances, for example, Euclidean distance. We also introduce a modified version of mutual information as a novel method for validating clustering results when the true clustering is known. The performance of our algorithm with respect to experimental noise is shown by extensive simulation studies. The algorithm is tested on a subset of 2029 cDNA clones coming from 15 different genes from a cDNA library derived from human dendritic cells. Furthermore, the clustering of these 2029 cDNA clones is demonstrated when the entire set of 76,032 cDNA clones is processed.

Construction of physical maps from oligonucleotide fingerprints data

A new algorithm for the construction of physical maps from hybridization fingerprints of short oligonucleotide probes has been developed. Extensive simulations in high-noise scenarios show that the algorithm produces an essentially completely correct map in over 95% of trials. Tests for the influence of specific experimental parameters demonstrate that the algorithm is robust to both false positive and false negative experimental errors. The algorithm was also tested in simulations using real DNA sequences of C. elegans, E. coli, S. cerevisiae, and H. sapiens. To overcome the non-randomness of probe frequencies in these sequences, probes were preselected based on sequence statistics and a screening process of the hybridization data was developed. With these modifications, the algorithm produced very encouraging results.

Toward the gene catalogue of sea urchin development: the construction and analysis of an unfertilized egg cDNA library highly normalized by oligonucleotide fingerprinting

We describe the use of oligonucleotide fingerprinting for the generation of a normalized cDNA library from unfertilized sea urchin eggs and report the preliminary analysis of this library, which resulted in the establishment of a partial gene catalogue of the sea urchin egg. In an analysis of 21,925 cDNA clones by hybridization with 217 oligonucleotide probes, we were able to identify 6291 clusters corresponding to different transcripts, ranging in size from 1 to 265 clones. This corresponds to an average 3.5-fold normalization of the starting library. The normalized library represents about one-third of all genes expressed in the sea urchin egg. To generate sequence information for the transcripts represented by the clusters, representative clones selected from 711 clusters were sequenced. The construction and preliminary analysis of the normalized library are the first steps in the assembly of an increasingly complete collection of maternal genes expressed in the sea urchin egg, which will provide a number of insights into the early development of this well-characterized model organism.