The TIGR Gene Indices: analysis of gene transcript sequences in highly sampled eukaryotic species

While genome sequencing projects are advancing rapidly, EST sequencing and analysis remains a primary research tool for the identification and categorization of gene sequences in a wide variety of species and an important resource for annotation of genomic sequence. The TIGR Gene Indices (http://www.tigr.org/tdb/tgi.shtml) are a collection of species-specific databases that use a highly refined protocol to analyze EST sequences in an attempt to identify the genes represented by that data and to provide additional information regarding those genes. Gene Indices are constructed by first clustering, then assembling EST and annotated gene sequences from GenBank for the targeted species. This process produces a set of unique, high-fidelity virtual transcripts, or Tentative Consensus (TC) sequences. The TC sequences can be used to provide putative genes with functional annotation, to link the transcripts to mapping and genomic sequence data, to provide links between orthologous and paralogous genes and as a resource for comparative sequence analysis.

Retroposon analysis of major cetacean lineages: The monophyly of toothed whales and the paraphyly of river dolphins

SINE (short interspersed element) insertion analysis elucidates contentious aspects in the phylogeny of toothed whales and dolphins (Odontoceti), especially river dolphins. Here, we characterize 25 informative SINEs inserted into unique genomic loci during evolution of odontocetes to construct a cladogram, and determine a total of 2.8 kb per taxon of the flanking sequences of these SINE loci to estimate divergence times among lineages. We demonstrate that: (i) Odontocetes are monophyletic; (ii) Ganges River dolphins, beaked whales, and ocean dolphins diverged (in this order) after sperm whales; (iii) three other river dolphin taxa, namely the Amazon, La Plata, and Yangtze river dolphins, form a monophyletic group with Yangtze River dolphins being the most basal; and (iv) the rapid radiation of extant cetacean lineages occurred some 28–33 million years B.P., in strong accord with the fossil record. The combination of SINE and flanking sequence analysis suggests a topology and set of divergence times for odontocete relationships, offering alternative explanations for several long-standing problems in cetacean evolution.

Analysis of the genetic differences between Neisseria meningitidis and Neisseria gonorrhoeae: two closely related bacteria expressing two different pathogenicities.

We have investigated genetic differences between the closely related pathogenic Neisseria species, Neisseria meningitidis and Neisseria gonorrhoeae, as a novel approach to the elucidation of the genetic basis for their different pathogenicities. N. meningitidis is a major cause of cerebrospinal meningitis, whereas N. gonorrhoeae is the agent of gonorrhoea. The technique of representational difference analysis was adapted to the search for genes present in the meningococcus but absent from the gonococcus. The libraries achieved are comprehensive and specific in that they contain sequences corresponding to the presently identified meningococcus-specific genes (capsule, frp, rotamase, and opc) but lack genes more or less homologous between the two species, e.g., ppk and pilC1. Of 35 randomly chosen clones specific to N. meningitidis, DNA sequence analysis has confirmed that the large majority have no homology with published neisserial sequences. Mapping of the cloned DNA fragments onto the chromosome of N. meningitidis strain Z2491 has revealed a nonrandom distribution of meningococcus-specific sequences. Most of the genetic differences between the meningococcus and gonococcus appear to be clustered in three distinct regions, one of which (region 1) contains the capsule-related genes. Region 3 was found only in strains of serogroup A, whereas region 2 is present in a variety of meningococci belonging to different serogroups. At a time when bacterial genomes are being sequenced, we believe that this technique is a powerful tool for a rapid and directed analysis of the genetic basis of inter- or intraspecific phenotypic variations.

Restriction fragment length polymorphism-coupled domain-directed differential display: a highly efficient technique for expression analysis of multigene families.

In this paper, a reverse-transcriptase PCR-based protocol suitable for efficient expression analysis of multigene families is presented. The method combines restriction fragment length polymorphism (RFLP) technology with a gene family-specific version of mRNA differential display and hence is called "RFLP-coupled domain-directed differential display. "With this method, expression of all members of a multigene family at many different developmental stages, in diverse tissues and even in different organisms, can be displayed on one gel. Moreover, bands of interest, representing gene family members, are directly accessible to sequence analysis, without the need for subcloning. The method thus enables a detailed, high-resolution expression analysis of known gene family members as well as the identification and characterization of new ones. Here the technique was used to analyze differential expression of MADS-box genes in male and female inflorescences of maize (Zea mays ssp. mays). Six different MADS-box genes could be identified, being either specifically expressed in the female sex or preferentially expressed in male or female inflorescences, respectively. Other possible applications of the method are discussed.

Genetic diversity among community methicillin-resistant Staphylococcus aureus strains causing outpatient infections in Australia

Increasing reports of the appearance of novel nonmultiresistant methicillin-resistant Staphylococcus aureus MRSA (MRSA) strains in the community and of the spread of hospital MRSA strains into the community are cause for public health concern. We conducted two national surveys of unique isolates of S. aureus from clinical specimens collected from nonhospitalized patients commencing in 2000 and 2002, respectively. A total of 11.7% of 2,498 isolates from 2000 and 15.4% of 2,486 isolates from 2002 were MRSA. Approximately 54% of the MRSA isolates were nonmultiresistant (resistant to less than three of nine antibiotics) in both surveys. The majority of multiresistant MRSA isolates in both surveys belonged to two strains (strains AUS-2 and AUS-3), as determined by pulsed-field gel electrophoresis (PFGE) and resistogram typing. The 3 AUS-2 isolates and 10 of the 11 AUS-3 isolates selected for multilocus sequence typing (MLST) and staphylococcal chromosomal cassette mec (SCCmec) analysis were ST239-MRSA-III (where ST is the sequence type) and thus belonged to the same clone as the eastern Australian MRSA strain of the 1980s, which spread internationally. Four predominant clones of novel nonmultiresistant MRSA were identified by PFGE, MLST, and SCCmec analysis: ST22-MRSA-IV (strain EMRSA-15), ST1-MRSA-IV (strain WA-1), ST30-MRSA-IV (strain SWP), and ST93-MRSA-IV (strain Queensland). The last three clones are associated with community acquisition. A total of 14 STs were identified in the surveys, including six unique clones of novel nonmultiresistant MRSA, namely, STs 73, 93, 129, 75, and 80sIv and a new ST. SCCmec types IV and V were present in diverse genetic backgrounds. These findings provide support for the acquisition of SCCmec by multiple lineages of S. aureus. They also confirm that both hospital and community strains of MRSA are now common in nonhospitalized patients throughout Australia.

Predicting the solvent accessibility of transmembrane residues from protein sequence

In this study, we propose a novel method to predict the solvent accessible surface areas of transmembrane residues. For both transmembrane alpha-helix and beta-barrel residues, the correlation coefficients between the predicted and observed accessible surface areas are around 0.65. On the basis of predicted accessible surface areas, residues exposed to the lipid environment or buried inside a protein can be identified by using certain cutoff thresholds. We have extensively examined our approach based on different definitions of accessible surface areas and a variety of sets of control parameters. Given that experimentally determining the structures of membrane proteins is very difficult and membrane proteins are actually abundant in nature, our approach is useful for theoretically modeling membrane protein tertiary structures, particularly for modeling the assembly of transmembrane domains. This approach can be used to annotate the membrane proteins in proteomes to provide extra structural and functional information.

Sequence variation in the Newcastle disease virus genome

Full-length genome sequences of five virulent and five avirulent strains of Newcastle disease virus isolated between 1998 and 2002 in Victoria and New South Wales, Australia were determined. Comparisons between these strains revealed that coding sequence variability in the haemagglutinin-neuraminidase (HN), matrix (M) and phosphoprotein (P) gene sequences appeared to be more variable than in the fusion (F), nucleocapsid (N) and RNA dependent-RNA replicase (L) genes. Sequence analysis of a number of other isolates made during the recent virulent NDV outbreaks, also identified the presence of a number of variants with altered F gene cleavage sites, which resulted in altered biological properties of those viruses. Quasispecies analysis of a number of field isolates indicated the presence of virulent virus in one particular isolate. Gene sequence analysis of the progenitor virus isolated in 1998 showed very little sequence variation when compared to that of a progenitor-like virus isolated in 2001 demonstrating that in the field. viral genome sequence variation appears to be biologically restricted to that of a consensus sequence. (c) 2005 Elsevier B.V. All rights reserved.

Methicillin-resistant Staphylococcus aureus genotyping using a small set of polymorphisms

The aim of this study was to identify a set of genetic polymorphisms that efficiently divides methicillin-resistant Staphylococcus aureus (MRSA) strains into groups consistent with the population structure. The rationale was that such polymorphisms could underpin rapid real-time PCR or low-density array-based methods for monitoring MRSA dissemination in a cost-effective manner. Previously, the authors devised a computerized method for identifying sets of single nucleoticle polymorphisms (SNPs) with high resolving power that are defined by multilocus sequence typing (MLST) databases, and also developed a real-time PCR method for interrogating a seven-member SNP set for genotyping S. aureus. Here, it is shown that these seven SNPs efficiently resolve the major MRSA lineages and define 27 genotypes. The SNP-based genotypes are consistent with the MRSA population structure as defined by eBURST analysis. The capacity of binary markers to improve resolution was tested using 107 diverse MRSA isolates of Australian origin that encompass nine SNP-based genotypes. The addition of the virulence-associated genes cna, pvl and bbplsdrE, and the integrated plasmids pT181, p1258 and pUB110, resolved the nine SNP-based genotypes into 21 combinatorial genotypes. Subtyping of the SCCmec locus revealed new SCCmec types and increased the number of combinatorial genotypes to 24. It was concluded that these polymorphisms provide a facile means of assigning MRSA isolates into well-recognized lineages.

Staphylococcus aureus genotyping using novel real-time PCR formats

One approach to microbial genotyping is to make use of sets of single-nucleotide polymorphisms (SNPs) in combination with binary markers. Here we report the modification and automation of a SNP-plus-binary-marker-based approach to the genotyping of Staphylococcus aureus and its application to 391 S. aureus isolates from southeast Queensland, Australia. The SNPs used were arcC210, tpi243, arcC162, gmk318, pta294, tpi36, tpi241, and pta383. These provide a Simpson's index of diversity (D) of 0.95 with respect to the S. aureus multilocus sequence typing database and define 61 genotypes and the major clonal complexes. The binary markers used were pvl, cna, sdrE, pT181, and pUB110. Two novel real-time PCR formats for interrogating these markers were compared. One of these makes use of light upon extension (LUX) primers and biplexed reactions, while the other is a streamlined modification of kinetic PCR using SYBR green. The latter format proved to be more robust. In addition, automated methods for DNA template preparation, reaction setup, and data analysis were developed. A single SNP-based method for ST-93 (Queensland clone) identification was also devised. The genotyping revealed the numerical importance of the South West Pacific and Queensland community-acquired methicillin-resistant S. aureus (MRSA) clones and the clonal complex 239 Aus-1/Aus-2 hospital-associated MRSA. There was a strong association between the community-acquired clones and pvl.

Methicillin-resistant Staphylococcus aureus in the Australian community: an evolving epidemic

Objective: To describe antimicrobial resistance and molecular epidemiology of methicillin-resistant Staphylococcus aureus (MRSA) isolated in community settings in Australia. Design and setting: Survey of S. aureus isolates collected prospectively Australia-wide between July 2004 and February 2005; results were compared with those of similar surveys conducted in 2000 and 2002. Main outcome measures: Up to 100 consecutive, unique clinical isolates of S. aureus from outpatient settings were collected at each of 22 teaching hospital and five private laboratories from cities in all Australian states and territories. They were characterised by antimicrobial susceptibilities (by agar dilution methods), coagulase gene typing, pulsed-field gel electrophoresis, multilocus sequence typing, SCCmec typing and polymerase chain reaction tests for Panton-Valentine leukocidin (PVL) gene. Results: 2652 S. aureus isolates were collected, of which 395 (14.9%) were MRSA. The number of community-associated MRSA (CA-MRSA) isolates rose from 4.7% (118/2498) of S. aureus isolates in 2000 to 7.3% (194/2652) in 2004 (P=0.001). Of the three major CA-MRSA strains, WA-1 constituted 45/257 (18%) of MRSA in 2000 and 64/395 (16%) in 2004 (P=0.89), while the Queensland (OLD) strain increased from 13/257 (5%) to 58/395 (15%) (P=0.0004), and the south-west Pacific (SWP) strain decreased from 33/257 (13%) to 26/395 (7%) (P=0.01). PVL genes were detected in 90/195 (46%) of CA-MRSA strains, including 5/64 (8%) of WA-1, 56/58 (97%) of OLD, and 25/26 (96%) of SWP strains. Among health care-associated MRSA strains, all AUS-2 and AUS-3 isolates were multidrug-resistant, and UK EMRSA-15 isolates were resistant to ciprofloxacin and erythromycin (50%) or to ciprofloxacin alone (44%). Almost all (98%) of CA-MRSA strains were non-multiresistant. Conclusions: Community-onset MRSA continues to spread throughout Australia. The hypervirulence determinant PVL is often found in two of the most common CA-MRSA strains. The rapid changes in prevalence emphasise the importance of ongoing surveillance.

Identification, mutagenesis, and transcriptional analysis of the methanesulfonate transport operon of methylosulfonomonas methylovora

Recently identified genes located downstream (3') of the msmEF (transport encoding) gene cluster, msmGH, and located 5' of the structural genes for methanesulfonate monooxygenase (MSAMO) are described from Methylosulfonomonas methylovora. Sequence analysis of the derived polypeptide sequences encoded by these genes revealed a high degree of identity to ABC-type transporters. MsmE showed similarity to a putative periplasmic substrate binding protein, MsmF resembled an integral membraneassociated protein, and MsmG was a putative ATP-binding enzyme. MsmH was thought to be the cognate permease component of the sulfonate transport system. The close association of these putative transport genes to the MSAMO structural genes msmABCD suggested a role for these genes in transport of methanesulfonic acid (MSA) into M. methylovora. msmEFGH and msmABCD constituted two operons for the coordinated expression of MSAMO and the MSA transporter systems. Reverse-transcription-PCR analysis of msmABCD and msmEFGH revealed differential expression of these genes during growth on MSA and methanol. The msmEFGH operon was constitutively expressed, whereas MSA induced expression of msmABCD. A mutant defective in msmE had considerably slower growth rates than the wild type, thus supporting the proposed role of MsmE in the transport of MSA into M. methylovora.

Transcriptomic Analysis of the Host Response and Innate Resilience to Enterotoxigenic Escherichia coli Infection in Humans.

BACKGROUND: Enterotoxigenic Escherichia coli (ETEC) is a globally prevalent cause of diarrhea. Though usually self-limited, it can be severe and debilitating. Little is known about the host transcriptional response to infection. We report the first gene expression analysis of the human host response to experimental challenge with ETEC. METHODS: We challenged 30 healthy adults with an unattenuated ETEC strain, and collected serial blood samples shortly after inoculation and daily for 8 days. We performed gene expression analysis on whole peripheral blood RNA samples from subjects in whom severe symptoms developed (n = 6) and a subset of those who remained asymptomatic (n = 6) despite shedding. RESULTS: Compared with baseline, symptomatic subjects demonstrated significantly different expression of 406 genes highlighting increased immune response and decreased protein synthesis. Compared with asymptomatic subjects, symptomatic subjects differentially expressed 254 genes primarily associated with immune response. This comparison also revealed 29 genes differentially expressed between groups at baseline, suggesting innate resilience to infection. Drug repositioning analysis identified several drug classes with potential utility in augmenting immune response or mitigating symptoms. CONCLUSIONS: There are statistically significant and biologically plausible differences in host gene expression induced by ETEC infection. Differential baseline expression of some genes may indicate resilience to infection.

Gene mapping and expression analysis of 16q loss of heterozygosity identifies WWOX and CYLD as being important in determining clinical outcome in multiple myeloma.

We performed fluorescent in situ hybridization (FISH) for 16q23 abnormalities in 861 patients with newly diagnosed multiple myeloma and identified deletion of 16q [del(16q)] in 19.5%. In 467 cases in which demographic and survival data were available, del(16q) was associated with a worse overall survival (OS). It was an independent prognostic marker and conferred additional adverse survival impact in cases with the known poor-risk cytogenetic factors t(4;14) and del(17p). Gene expression profiling and gene mapping using 500K single-nucleotide polymorphism (SNP) mapping arrays revealed loss of heterozygosity (LOH) involving 3 regions: the whole of 16q, a region centered on 16q12 (the location of CYLD), and a region centered on 16q23 (the location of the WW domain-containing oxidoreductase gene WWOX). CYLD is a negative regulator of the NF-kappaB pathway, and cases with low expression of CYLD were used to define a "low-CYLD signature." Cases with 16q LOH or t(14;16) had significantly reduced WWOX expression. WWOX, the site of the translocation breakpoint in t(14;16) cases, is a known tumor suppressor gene involved in apoptosis, and we were able to generate a "low-WWOX signature" defined by WWOX expression. These 2 genes and their corresponding pathways provide an important insight into the potential mechanisms by which 16q LOH confers poor prognosis.