Comparison of automated candidate gene prediction systems using genes implicated in type 2 diabetes by genome-wide association studies

Background
Automated candidate gene prediction systems allow geneticists to hone in on disease genes more rapidly by identifying the most probable candidate genes linked to the disease phenotypes under investigation. Here we assessed the ability of eight different candidate gene prediction systems to predict disease genes in intervals previously associated with type 2 diabetes by benchmarking their performance against genes implicated by recent genome-wide association studies.

Results
Using a search space of 9556 genes, all but one of the systems pruned the genome in favour of genes associated with moderate to highly significant SNPs. Of the 11 genes associated with highly significant SNPs identified by the genome-wide association studies, eight were flagged as likely candidates by at least one of the prediction systems. A list of candidates produced by a previous consensus approach did not match any of the genes implicated by 706 moderate to highly significant SNPs flagged by the genome-wide association studies. We prioritized genes associated with medium significance SNPs.

Conclusion
The study appraises the relative success of several candidate gene prediction systems against independent genetic data. Even when confronted with challengingly large intervals, the candidate gene prediction systems can successfully select likely disease genes. Furthermore, they can be used to filter statistically less-well-supported genetic data to select more likely candidates. We suggest consensus approaches fail because they penalize novel predictions made from independent underlying databases. To realize their full potential further work needs to be done on prioritization and annotation of genes.

Identification and disruption of the gene encoding the third member of the low-molecular-mass rhoptry complex in Plasmodium falciparum

The low-molecular-mass rhoptry complex of Plasmodium falciparum consists of three proteins, rhoptry-associated protein 1 (RAP1), RAP2, and RAP3. The genes encoding RAP1 and RAP2 are known; however, the RAP3 gene has not been identified. In this study we identify the RAP3 gene from the P. falciparum genome database and show that this protein is part of the low-molecular-mass rhoptry complex. Disruption of RAP3 demonstrated that it is not essential for merozoite invasion, probably because RAP2 can complement the loss of RAP3. RAP3 has homology with RAP2, and the genes are encoded on chromosome 5 in a head-to-tail fashion. Analysis of the genome databases has identified homologous genes in all Plasmodium spp., suggesting that this protein plays a role in merozoite invasion. The region surrounding the RAP3 homologue in the Plasmodium yoelii genome is syntenic with the same region in P. falciparum; however, there is a single gene. Phylogenetic comparison of the RAP2/3 protein family from Plasmodium spp. suggests that the RAP2/3 duplication occurred after divergence of these parasite species.

Analysis of genome-wide association study data using the protein knowledge base

Article 98

The conformation of the mature dimeric human immunodeficiency virus type 1 RNA genome requires packaging of pol protein

The packaging of a mature dimeric RNA genome is an essential step in human immunodeficiency virus type 1 (HIV-1) replication. We have previously shown that overexpression of a protease (PR)-inactive HIV-1 Gag-Pro-Pol precursor protein generates noninfectious virions that contain mainly monomeric RNA (M. Shehu-Xhilaga, S. M. Crowe, and J. Mak, J. Virol. 75:1834-1841, 2001). To further define the contribution of HIV-1 Gag and Gag-Pro-Pol to RNA maturation, we analyzed virion RNA dimers derived from Gag particles in the absence of Gag-Pro-Pol. Compared to wild-type (WT) dimeric RNAs, these RNA dimers have altered mobility and low stability under electrophoresis conditions, suggesting that the HIV-1 Gag precursor protein alone is not sufficient to stabilize the dimeric virion RNA structure. The inclusion of an active viral PR, without reverse transcriptase (RT) and integrase (IN), rescued the stability of the virion RNA dimers in the Gag particles but did not restore the mobility of the RNAs, suggesting that RT and IN are also required for virion RNA dimer maturation. Thin-section electron microscopy showed that viral particles deficient in RT and IN contain empty cone-shaped cores. The abnormal core structure indicates a requirement for Gag-Pro-Pol packaging during core maturation. Supplementing viral particles with either RT or IN via Vpr-RT or Vpr-IN alone did not correct the conformation of the dimer RNAs, whereas expression of both RT and IN in trans as a Vpr-RT-IN fusion restored RNA dimer conformation to that of the WT virus and also restored the electron-dense, cone-shaped virion core characteristic of WT virus. Our data suggest a role for RT-IN in RNA dimer conformation and the formation of the electron-dense viral core.

Complete genome sequence of the Medicago microsymbiont Ensifer (Sinorhizobium) medicae strain WSM419

Ensifer (Sinorhizobium) medicae is an effective nitrogen fixing microsymbiont of a diverse range of annual Medicago (medic) species. Strain WSM419 is an aerobic, motile, non-spore forming, Gram-negative rod isolated from a M. murex root nodule collected in Sardinia, Italy in 1981. WSM419 was manufactured commercially in Australia as an inoculant for annual medics during 1985 to 1993 due to its nitrogen fixation, saprophytic competence and acid tolerance properties. Here we describe the basic features of this organism, together with the complete genome sequence, and annotation. This is the first report of a complete genome sequence for a microsymbiont of the group of annual medic species adapted to acid soils. We reveal that its genome size is 6,817,576 bp encoding 6,518 protein-coding genes and 81 RNA only encoding genes. The genome contains a chromosome of size 3,781,904 bp and 3 plasmids of size 1,570,951, 1,245,408 and 219,313 bp. The smallest plasmid is a feature unique to this medic microsymbiont.

Complete genome sequence of Rhizobium leguminosarum bv trifolii strain WSM2304, an effective microsymbiont of the South American clover Trifolium polymorphum

Rhizobium leguminosarum bv. trifolii is the effective nitrogen fixing microsymbiont of a diverse range of annual and perennial Trifolium (clover) species. Strain WSM2304 is an aerobic, motile, non-spore forming, Gram-negative rod isolated from Trifolium polymorphum in Uruguay in 1998. This microsymbiont predominated in the perennial grasslands of Glencoe Research Station, in Uruguay, to competitively nodulate its host, and fix atmospheric nitrogen. Here we describe the basic features of WSM2304, together with the complete genome sequence, and annotation. This is the first completed genome sequence for a nitrogen fixing microsymbiont of a clover species from the American centre of origin. We reveal that its genome size is 6,872,702 bp encoding 6,643 protein-coding genes and 62 RNA only encoding genes. This multipartite genome was found to contain 5 distinct replicons; a chromosome of size 4,537,948 bp and four circular plasmids of size 4,537,948, 1,266,105, 501,946, 308,747 and 257,956 bp.

Complete genome sequence of rhizobium leguminosarum bv. trifolii strain WSM1325, an effective microsymbiont of annual Mediterranean clovers

Rhizobium leguminosarum bv trifolii is a soil-inhabiting bacterium that that has the capacity to be an effective nitrogen fixing microsymbiont of a diverse range of annual Trifolium (clover) species. Strain WSM1325 is an aerobic, motile, non-spore forming, Gram-negative rod isolated from root nodules collected in 1993 from the Greek Island of Serifos. WSM1325 is manufactured commercially in Australia as an inoculant for a broad range of annual clovers of Mediterranean origin due to its superior attributes of saprophytic competence, nitrogen fixation and acid-tolerance. Here we describe the basic features of this organism, together with the complete genome sequence, and annotation. This is the first completed genome sequence for a microsymbiont of annual clovers. We reveal that its genome size is 7,418,122 bp encoding 7,232 protein-coding genes and 61 RNA-only encoding genes. This multipartite genome contains 6 distinct replicons; a chromosome of size 4,767,043 bp and 5 plasmids of size 828,924, 660,973, 516,088, 350,312 and 294,782 bp.

Consideration of multiple template switches is critical to the determination of the recombination rate on the HIV genome

Retroviral recombination drives viral diversity and facilitates the emergence of immune escape and drug resistant mutants that contribute to disease progression. Current estimates of retroviral recombination rates are based on indirect measurements that do not take into account the effects of multiple recombination events. In the presence of multiple template switches, any even number of template switches result in no observed recombination and any odd number is detected as a single recombination event. We demonstrate that ignoring multiple recombination events consistently underestimates the true recombination rate, especially over large genetic distances and high rates of recombination. Here, we present a novel approach to measure rates of recombination across different gene segments regardless of the effects of genetic distance and the overall rate of recombination. We apply these tools to a novel HIV-1 marker system, which mimics the recombination process between closely related genomes, analogous to those found within the quasispecies of an infected individual. We directly measure the recombination rate in gag, correcting for the effects of multiple template switches and background recombination. Furthermore, our analysis indicates that recombination rates are likely to vary across the viral genome. This system is applicable to other studies to accurately measure the recombination rate that is critical for the diversification of retroviruses.

Extracting DNA from museum bird eggs, and whole genome amplification of archive DNA

We present a comprehensive protocol for extracting DNA from egg membranes and other internal debris recovered from the interior of blown museum bird eggs. A variety of commercially available DNA extraction methods were found to be applicable. DNA sequencing of polymerase chain reaction (PCR) products for a 176-bp fragment of mitochondrial DNA was successful for most egg samples (> 78%) even though the amount of DNA extracted (mean = 14.71 ± 4.55 ng/µL) was significantly less than that obtained for bird skin samples (mean = 67.88 ± 4.77 ng/µL). For PCR and sequencing of snipe (Gallinago) DNA, we provide eight new primers for the ‘DNA barcode’ region of COI mtDNA. In various combinations, the primers target a range of PCR products sized from 72 bp to the full ‘barcode’ of 751 bp. Not all possible combinations were tested with archive snipe DNA, but we found a significantly better success rate of PCR amplification for a shorter 176-bp target compared with a larger 288-bp fragment (67% vs. 39%). Finally, we explored the feasibility of whole genome amplification (WGA) for extending the use of archive DNA in PCR and sequencing applications. Of two WGA approaches, a PCR-based method was found to be able to amplify whole genomic DNA from archive skins and eggs from museum bird collections. After WGA, significantly more archive egg samples produced visible PCR products on agarose (56.9% before WGA vs. 79.0% after WGA). However, overall sequencing success did not improve significantly (78.8% compared with 83.0%).

Candidate disease gene prediction using Gentrepid: application to a genome-wide association study on coronary artery disease

Current single-locus-based analyses and candidate disease gene prediction methodologies used in genome-wide association studies (GWAS) do not capitalize on the wealth of the underlying genetic data, nor functional data available from molecular biology. Here, we analyzed GWAS data from the Wellcome Trust Case Control Consortium (WTCCC) on coronary artery disease (CAD). Gentrepid uses a multiple-locus-based approach, drawing on protein pathway- or domain-based data to make predictions. Known disease genes may be used as additional information (seeded method) or predictions can be based entirely on GWAS single nucleotide polymorphisms (SNPs) (ab initio method). We looked in detail at specific predictions made by Gentrepid for CAD and compared these with known genetic data and the scientific literature. Gentrepid was able to extract known disease genes from the candidate search space and predict plausible novel disease genes from both known and novel WTCCC-implicated loci. The disease gene candidates are consistent with known biological information. The results demonstrate that this computational approach is feasible and a valuable discovery tool for geneticists.