Widespread divergence between incipient Anopheles gambiae species revealed by whole genome sequences.

The Afrotropical mosquito Anopheles gambiae sensu stricto, a major vector of malaria, is currently undergoing speciation into the M and S molecular forms. These forms have diverged in larval ecology and reproductive behavior through unknown genetic mechanisms, despite considerable levels of hybridization. Previous genome-wide scans using gene-based microarrays uncovered divergence between M and S that was largely confined to gene-poor pericentromeric regions, prompting a speciation-with-ongoing-gene-flow model that implicated only about 3% of the genome near centromeres in the speciation process. Here, based on the complete M and S genome sequences, we report widespread and heterogeneous genomic divergence inconsistent with appreciable levels of interform gene flow, suggesting a more advanced speciation process and greater challenges to identify genes critical to initiating that process.

A survey of SL1-spliced transcripts from the root-lesion nematode Pratylenchus penetrans.

Plant-parasitic nematodes are important and cosmopolitan pathogens of crops. Here, we describe the generation and analysis of 1928 expressed sequence tags (ESTs) of a splice-leader 1 (SL1) library from mixed life stages of the root-lesion nematode Pratylenchus penetrans. The ESTs were grouped into 420 clusters and classified by function using the Gene Ontology (GO) hierarchy and the Kyoto KEGG database. Approximately 80% of all translated clusters show homology to Caenorhabditis elegans proteins, and 37% of the C. elegans gene homologs had confirmed phenotypes as assessed by RNA interference tests. Use of an SL1-PCR approach, while ensuring the cloning of the 5' ends of mRNAs, has demonstrated bias toward short transcripts. Putative nematode-specific and Pratylenchus -specific genes were identified, and their implications for nematode control strategies are discussed.

An expression profile of human pancreatic islet mRNAs by Serial Analysis of Gene Expression (SAGE).

AIMS/HYPOTHESIS: The Human Genome Project seeks to identify all genes with the ultimate goal of evaluation of relative expression levels in physiology and in disease states. The purpose of the current study was the identification of the most abundant transcripts in human pancreatic islets and their relative expression levels using Serial Analysis of Gene Expression. METHODS: By cutting cDNAs into small uniform fragments (tags) and concatemerizing them into larger clones, the identity and relative abundance of genes can be estimated for a cDNA library. Approximately 49,000 SAGE tags were obtained from three human libraries: (i) ficoll gradient-purified islets (ii) islets further individually isolated by hand-picking, and (iii) pancreatic exocrine tissue. RESULTS: The relative abundance of each of the genes identified was approximated by the frequency of the tags. Gene ontology functions showed that all three libraries contained transcripts mostly encoding secreted factors. Comparison of the two islet libraries showed various degrees of contamination from the surrounding exocrine tissue (11 vs 25%). After removal of exocrine transcripts, the relative abundance of 2180 islet transcripts was determined. In addition to the most common genes (e.g. insulin, transthyretin, glucagon), a number of other abundant genes with ill-defined functions such as proSAAS or secretagogin, were also observed. CONCLUSION/INTERPRETATION: This information could serve as a resource for gene discovery, for comparison of transcript abundance between tissues, and for monitoring gene expression in the study of beta-cell dysfunction of diabetes. Since the chromosomal location of the identified genes is known, this SAGE expression data can be used in setting priorities for candidate genes that map to linkage peaks in families affected with diabetes.

Complete genome sequence of Salmonella enterica serovar Typhimurium LT2.

Salmonella enterica subspecies I, serovar Typhimurium (S. typhimurium), is a leading cause of human gastroenteritis, and is used as a mouse model of human typhoid fever. The incidence of non-typhoid salmonellosis is increasing worldwide, causing millions of infections and many deaths in the human population each year. Here we sequenced the 4,857-kilobase (kb) chromosome and 94-kb virulence plasmid of S. typhimurium strain LT2. The distribution of close homologues of S. typhimurium LT2 genes in eight related enterobacteria was determined using previously completed genomes of three related bacteria, sample sequencing of both S. enterica serovar Paratyphi A (S. paratyphi A) and Klebsiella pneumoniae, and hybridization of three unsequenced genomes to a microarray of S. typhimurium LT2 genes. Lateral transfer of genes is frequent, with 11% of the S. typhimurium LT2 genes missing from S. enterica serovar Typhi (S. typhi), and 29% missing from Escherichia coli K12. The 352 gene homologues of S. typhimurium LT2 confined to subspecies I of S. enterica-containing most mammalian and bird pathogens-are useful for studies of epidemiology, host specificity and pathogenesis. Most of these homologues were previously unknown, and 50 may be exported to the periplasm or outer membrane, rendering them accessible as therapeutic or vaccine targets.

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.

Complete genome sequence of an M1 strain of Streptococcus pyogenes.

The 1,852,442-bp sequence of an M1 strain of Streptococcus pyogenes, a Gram-positive pathogen, has been determined and contains 1,752 predicted protein-encoding genes. Approximately one-third of these genes have no identifiable function, with the remainder falling into previously characterized categories of known microbial function. Consistent with the observation that S. pyogenes is responsible for a wider variety of human disease than any other bacterial species, more than 40 putative virulence-associated genes have been identified. Additional genes have been identified that encode proteins likely associated with microbial "molecular mimicry" of host characteristics and involved in rheumatic fever or acute glomerulonephritis. The complete or partial sequence of four different bacteriophage genomes is also present, with each containing genes for one or more previously undiscovered superantigen-like proteins. These prophage-associated genes encode at least six potential virulence factors, emphasizing the importance of bacteriophages in horizontal gene transfer and a possible mechanism for generating new strains with increased pathogenic potential.

Comparison of the Escherichia coli K-12 genome with sampled genomes of a Klebsiella pneumoniae and three salmonella enterica serovars, Typhimurium, Typhi and Paratyphi.

The Escherichia coli K-12 genome (ECO) was compared with the sampled genomes of the sibling species Salmonella enterica serovars Typhimurium, Typhi and Paratyphi A (collectively referred to as SAL) and the genome of the close outgroup Klebsiella pneumoniae (KPN). There are at least 160 locations where sequences of >400 bp are absent from ECO but present in the genomes of all three SAL and 394 locations where sequences are present in ECO but close homologs are absent in all SAL genomes. The 394 sequences in ECO that do not occur in SAL contain 1350 (30.6%) of the 4405 ECO genes. Of these, 1165 are missing from both SAL and KPN. Most of the 1165 genes are concentrated within 28 regions of 10-40 kb, which consist almost exclusively of such genes. Among these regions were six that included previously identified cryptic phage. A hypothetical ancestral state of genomic regions that differ between ECO and SAL can be inferred in some cases by reference to the genome structure in KPN and the more distant relative Yersinia pestis. However, many changes between ECO and SAL are concentrated in regions where all four genera have a different structure. The rate of gene insertion and deletion is sufficiently high in these regions that the ancestral state of the ECO/SAL lineage cannot be inferred from the present data. The sequencing of other closely related genomes, such as S.bongori or Citrobacter, may help in this regard.

Novel sequence organization and insertion specificity of IS605 and IS606: chimaeric transposable elements of Helicobacter pylori.

IS605, an insertion sequence (IS) that is unusual in containing homologs of genes for the single putative transposases of two other unrelated IS elements (IS200 and IS1341), was found in nearly one-third of a set of 238 independent isolates of the gastric pathogen Helicobacter pylori. Hybridization and PCR tests indicated that any strain carrying one of these ORFs also carried the other, which implies that both ORFs are in the same unit of transposition. The IS605 ends and target sites for insertion were identified by sequencing eight preexisting insertions in strain NCTC11638, corresponding empty sites in other strains, and new transpositions in E. coli of an IS605 derivative marked with a selectable chloramphenicol-resistance gene. These tests showed that IS605 is also unusual in: (1) having unique, not inverted repeat, ends; (2) not duplicating (or deleting) target sequences during transposition; and (3) inserting with its left (IS200-homolog) end next to 5'-TTTAA or 5'-TTTAAC. IS605 was implicated in at least two genome rearrangements in strain NCTC11638. A second member of the IS605 family, called IS606 (25% amino acid identity to IS605 in inferred proteins) was found in one-third of 38 H. pylori strains tested, many of which did not carry IS605. The features of these two chimaeric IS elements are discussed in terms of possible transposition mechanisms, IS element evolution, and effects of IS elements on genome organization and evolution in the microbes that they inhabit.

Cloning of crustacean ecdysteroid receptor and retinoid-X receptor gene homologs and elevation of retinoid-X receptor mRNA by retinoic acid.

We report the cloning and analysis of ecdysteroid receptor (bpEcR) and retinoid-X receptor (UpRXR) cDNA homologs from the fiddler crab Uca pugilator. The deduced amino acid sequence of this crustacean EcR most closely resembles the insect EcRs within the DNA binding and ligand binding domains (LBDs). For UpRXR, the DNA binding domain (DBD) shares greatest identity to the insect USPs. The ligand binding domain, however, is closer to vertebrate RXRs but may have a nonfunctional AF-2 domain. Probes derived from these clones were used to examine transcript levels in blastemas during early limb regeneration. Both UpEcR and UpRXR transcripts were detected in low amounts 1 day after limb loss, but increased during the next 4 days. Immersion of crabs in sea water containing all-trans retinoic acid increased the steady state concentrations of UpRXR transcript and altered the pattern of circulating ecdysteroids. These effects correlate with the disruptive effects of retinoic acid on blastemal differentiation observed in earlier studies.

Analyses of the cag pathogenicity island of Helicobacter pylori.

Most strains of Helicobacter pylori from patients with peptic ulcer disease or intestinal-type gastric cancer carry cagA, a gene that encodes an immunodominant protein of unknown function, whereas many of the strains from asymptomatically infected persons lack this gene. Recent studies showed that the cagA gene lies near the right end of a approximately 37kb DNA segment (a pathogenicity island, or PAI) that is unique to cagA+ strains and that the cag PAI was split in half by a transposable element insertion in the reference strain NCTC11638. In complementary experiments reported here, we also found the same cag PAI, and sequenced a 39 kb cosmid clone containing the left 'cagII' half of this PAI. Encoded in cagII were four proteins each with homology to four components of multiprotein complexes of Bordetella pertussis ('Ptl'), Agrobacterium tumefaciens ('Vir'), and conjugative plasmids ('Tra') that help deliver pertussis toxin and T (tumour inducing) and plasmid DNA, respectively, to target eukaryotic or prokaryotic cells, and also homologues of eukaryotic proteins that are involved in cytoskeletal structure. To the left of cagII in this cosmid were genes for homologues of HsIU (heat-shock protein) and Era (essential GTPase); to the right of cagII were homologues of genes for a type I restriction endonuclease and ion transport functions. Deletion of the cag PAI had no effect on synthesis of the vacuolating cytotoxin, but this deletion and several cag insertion mutations blocked induction of synthesis of proinflammatory cytokine IL-8 in gastric epithelial cells. Comparisons among H. pylori strains indicated that cag PAI gene content and arrangement are rather well conserved. We also identified two genome rearrangements with end-points in the cag PAI. One, in reference strain NCTC11638, involved IS605, a recently described transposable element (as also found by others). Another rearrangement, in 3 of 10 strains tested (including type strain NCTC11637), separated the normally adjacent cagA and picA genes and did not involve IS605. Our results are discussed in terms of how cag-encoded proteins might help trigger the damaging inflammatory responses in the gastric epithelium and possible contributions of DNA rearrangements to genome evolution.

Characterization of the mouse beta-prime adaptin gene; cDNA sequence, genomic structure, and chromosomal localization.

Adaptins are important subunits of heterotetrameric complexes called adaptors, which participate in the clathrin-coated, vesicle-mediated endocytosis and intracellular receptor transport. The gene family of adaptins is divided into three classes, alpha, beta, and gamma, with further subdivision into beta- and beta-prime components. Two beta-prime adaptins, the rat AP105a and the human BAM22, have previously been characterized. The BAM22 gene is located on human Chromosome (Chr) 22q12 and can be considered a candidate meningioma tumor suppressor gene. We report here the characterization of the mouse ortholog of the BAM22 gene, and we suggest the name adtb1 for the mouse gene. Like the BAM22 gene, the adtb1 transcript is highly and ubiquitously expressed. We provide 3885-bp cDNA sequence, which entirely covers the open reading frame of the adtb1, capable of encoding a protein of 943 amino acids. The adtb1 protein is highly conserved (>96% identity) when compared with AP105a and BAM22 proteins. We also report the genomic organization of adtb1, which is similar to the BAM22 gene. The adtb1 gene has been assigned to mouse Chr 11, band 11A2, which confirms the synteny between human Chr 22q12 and mouse Chr 11.

Multicellular ascomycetous fungal genomes contain more than 8000 genes.

Fungi comprise a large monophyletic group of uni- and multicellular eukaryotic organisms in which many species are of economic or medical importance. Fungal genomes are variable in size (13-42 Mb), and multicellular species support true spatial and temporal cell-type-specific regulation of gene expression. In a 38.8-kb Aspergillus nidulans contiguous genomic DNA region, a transposable element and 12 potential genes were identified, 7 similar to genes in other organisms. This observation is consistent with the prediction that multicellular ascomycetous fungi harbor 8000-9000 genes in a 36-Mb average genome. Thus, the genomic DNA sequence of filamentous fungi will provide substantial amounts of genetic and functional information that is not available in yeast, for the human and other metazoan minimal gene complement.

Structure of the promoter and genomic organization of the human beta'-adaptin gene (BAM22) from chromosome 22q12.

Adaptins are major structural components of heterotetrameric protein complexes called adaptors, which are essential in intracellular receptor transport via clathrin-coated vesicles. beta-adaptins constitute one of three known classes (alpha, beta, gamma) of adaptins, including beta and beta' subtypes. We previously cloned the human beta'-adaptin gene (BAM22) (GDB symbol, ADTB1) from chromosome 22q12 and proposed its involvement in the development of meningiomas. Here we describe the genomic organization of this gene, which consists of 22 exons spanning over approximately 100 kb. We also report results from point mutation screening of 7 randomly chosen exons analyzed in 110 sporadic meningiomas. As part of the genomic characterization of the BAM22 locus, we sequenced 40 kb covering exons 1-4 and 12 kb upstream from the start of gene transcription. Analysis of the sequence suggests that the BAM22 gene has a CG-rich promoter.

Sequence and analysis of the human ABL gene, the BCR gene, and regions involved in the Philadelphia chromosomal translocation.

The complete human BCR gene (152-141 nt) on chromosome 22 and greater than 80% of the human ABL gene (179-512 nt) on chromosome 9 have been sequenced from mapped cosmid and plasmid clones via a shotgun strategy. Because these two chromosomes are translocated with breakpoints within the BCR and ABL genes in Philadelphia chromosome-positive leukemias, knowledge of these sequences also might provide insight into the validity of various theories of chromosomal rearrangements. Comparison of these genes with their cDNA sequences reveal the positions of 23 BCR exons and putative alternative BCR first and second exons, as well as the common ABL exons 2-11, respectively. Additionally, these regions include the alternative ABL first exons 1b and 1a, a new gene 5' to the first ABL exon, and an open reading frame with homology to an EST within the BCR fourth intron. Further analysis reveals an Alu homology of 38.83 and 39.35% for the BCR and ABL genes, respectively, with other repeat elements present to a lesser extent. Four new Philadelphia chromosome translocation breakpoints from chronic myelogenous leukemia patients also were sequenced, and the positions of these and several other previously sequenced breakpoints now have been mapped precisely, although no consistent breakpoint features immediately were apparent. Comparative analysis of genomic sequences encompassing the murine homologues to the human ABL exons 1b and 1a, as well as regions encompassing the ABL exons 2 and 3, reveals that although there is a high degree of homology in their corresponding exons and promoter regions, these two vertebrate species show a striking lack of homology outside these regions.

The esg locus of Myxococcus xanthus encodes the E1 alpha and E1 beta subunits of a branched-chain keto acid dehydrogenase.

The esg locus of Myxococcus xanthus appears to control the production of a signal that must be transmitted between cells for the completion of multicellular development. DNA sequence analysis suggested that the esg locus encodes the E1 decarboxylase (composed of E1 alpha and E1 beta subunits) of a branched-chain keto acid dehydrogenase (BCKAD) that is involved in branched-chain amino acid (BCAA) metabolism. The properties of an esg::Tn5 insertion mutant supported this conclusion. These properties include: (i) the growth yield of the mutant was reduced with increasing concentrations of the BCAAs in the medium while the growth yield of wild-type cells increased, (ii) mutant extracts were deficient in BCKAD activity, and (iii) growth of the mutant in media with short branched-chain fatty acids related to the expected products of the BCKAD helped to correct the mutant defects in growth, pigmentation and development. The esg BCKAD appears to be involved in the synthesis of long branched-chain fatty acids since the mutant contained reduced levels of this class of compounds. Our results are consistent with a model in which the esg-encoded enzyme is involved in the synthesis of branched-chain fatty acids during vegetative growth, and these compounds are used later in cell-cell signalling during development.

Sequence of the rec-2 locus of Haemophilus influenzae: homologies to comE-ORF3 of Bacillus subtilis and msbA of Escherichia coli.

The nucleotide sequence of a 4243-bp PstI fragment containing the rec-2 gene of Haemophilus influenzae was determined. The amino acid (aa) sequences of four putative proteins were deduced from the corresponding open reading frames (ORFs). The 2400-bp ORF2 accounted for rec-2, based on the sequences of DNA fragments that contained rec-2::mini-Tn10 mutations. The rec-2 gene encoded a putative 800-aa protein with a M(r) of 90,561. Sequence analysis suggested that the rec-2 product contained nine highly probable integral membrane-spanning segments. Database searches showed that rec-2 was homologous to the comE-ORF3 gene of Bacillus subtilis. This hypothesis is consistent with the known involvement of both of these genes in the passage of transforming DNA through the competent-cell envelope. Although the sequences of the other three ORFs were incomplete, sufficient data were available to allow inferences about their homologies to other genes. ORF4, which overlapped ORF1, was homologous to the Escherichia coli dnaK suppressor gene, dksA, and therefore was named dsh-1 (dnaK suppressor homolog). Mutations in dsh-1 and its putative promoter region caused a mild sensitivity to UV light, but did not affect DNA recombination. ORF3, located downstream from rec-2, was homologous to msbA, an essential gene of E. coli with extensive similarity to the ATP-dependent translocators. ORF3 was named msh-1 (msbA homolog). Mutations in msh-1 had no effects on genetic transformation. The close juxtaposition of rec-2 and msh-1 implied that the expression of msh-1 could be linked to the translation of the rec-2 ORF.