Jeremy D. DeBarry

MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity

MCScan is an algorithm able to scan multiple genomes or subgenomes in order to identify putative homologous chromosomal regions, and align these regions using genes as anchors. The MCScanX toolkit implements an adjusted MCScan algorithm for detection of synteny and collinearity that extends the original software by incorporating 14 utility programs for visualization of results and additional downstream analyses. Applications of MCScanX to several sequenced plant genomes and gene families are shown as examples. MCScanX can be used to effectively analyze chromosome structural changes, and reveal the history of gene family expansions that might contribute to the adaptation of lineages and taxa. An integrated view of various modes of gene duplication can supplement the traditional gene tree analysis in specific families. The source code and documentation of MCScanX are freely available at http://chibba.pgml.uga.edu/mcscan2/.

De novo genome sequencing and comparative genomics of date palm (Phoenix dactylifera)

Date palm is one of the most economically important woody crops cultivated in the Middle East and North Africa and is a good candidate for improving agricultural yields in arid environments. Nonetheless, long generation times (5–8 years) and dioecy (separate male and female trees) have complicated its cultivation and genetic analysis. To address these issues, we assembled a draft genome for a Khalas variety female date palm, the first publicly available resource of its type for a member of the order Arecales. The ∼380 Mb sequence, spanning mainly gene-rich regions, includes >25,000 gene models and is predicted to cover ∼90% of genes and ∼60% of the genome. Sequencing of eight other cultivars, including females of the Deglet Noor and Medjool varieties and their backcrossed males, identified >3.5 million polymorphic sites, including >10,000 genic copy number variations. A small subset of these polymorphisms can distinguish multiple varieties. We identified a region of the genome linked to gender and found evidence that date palm employs an XY system of gender inheritance.

Sequence Analysis of Bacterial Artificial Chromosome Clones from the Apospory-Specific Genomic Region of Pennisetum and Cenchrus

Apomixis, asexual reproduction through seed, is widespread among angiosperm families. Gametophytic apomixis in Pennisetum squamulatum and Cenchrus ciliaris is controlled by the apospory-specific genomic region (ASGR), which is highly conserved and macrosyntenic between these species. Thirty-two ASGR bacterial artificial chromosomes (BACs) isolated from both species and one ASGR-recombining BAC from P. squamulatum, which together cover approximately 2.7 Mb of DNA, were used to investigate the genomic structure of this region. Phrap assembly of 4,521 high-quality reads generated 1,341 contiguous sequences (contigs; 730 from the ASGR and 30 from the ASGR-recombining BAC in P. squamulatum, plus 580 from the C. ciliaris ASGR). Contigs containing putative protein-coding regions unrelated to transposable elements were identified based on protein similarity after Basic Local Alignment Search Tool X analysis. These putative coding regions were further analyzed in silico with reference to the rice (Oryza sativa) and sorghum (Sorghum bicolor) genomes using the resources at Gramene (www.gramene.org) and Phytozome (www.phytozome.net) and by hybridization against sorghum BAC filters. The ASGR sequences reveal that the ASGR (1) contains both gene-rich and gene-poor segments, (2) contains several genes that may play a role in apomictic development, (3) has many classes of transposable elements, and (4) does not exhibit large-scale synteny with either rice or sorghum genomes but does contain multiple regions of microsynteny with these species.

The dynamics of LTR retrotransposon accumulation across 25 million years of panicoid grass evolution

Sample sequence analysis was employed to investigate the repetitive DNAs that were most responsible for the evolved variation in genome content across seven panicoid grasses with >5-fold variation in genome size and different histories of polyploidy. In all cases, the most abundant repeats were LTR retrotransposons, but the particular families that had become dominant were found to be different in the Pennisetum, Saccharum, Sorghum and Zea lineages. One element family, Huck, has been very active in all of the studied species over the last few million years. This suggests the transmittal of an active or quiescent autonomous set of Huck elements to this lineage at the founding of the panicoids. Similarly, independent recent activity of Ji and Opie elements in Zea and of Leviathan elements in Sorghum and Saccharum species suggests that members of these families with exceptional activation potential were present in the genome(s) of the founders of these lineages. In a detailed analysis of the Zea lineage, the combined action of several families of LTR retrotransposons were observed to have approximately doubled the genome size of Zea luxurians relative to Zea mays and Zea diploperennis in just the last few million years. One of the LTR retrotransposon amplification bursts in Zea may have been initiated by polyploidy, but the great majority of transposable element activations are not. Instead, the results suggest random activation of a few or many LTR retrotransposons families in particular lineages over evolutionary time, with some families especially prone to future activation and hyper-amplification.Heredity (2013) 110, 194–204; doi:10.1038/hdy.2012.99

Theileria parasites secrete a prolyl isomerase to maintain host leukocyte transformation

Infectious agents develop intricate mechanisms to interact with host cell pathways and hijack their genetic and epigenetic machinery to change host cell phenotypic states. Among the Apicomplexa phylum of obligate intracellular parasites, which cause veterinary and human diseases, Theileria is the only genus that transforms its mammalian host cells. Theileria infection of bovine leukocytes induces proliferative and invasive phenotypes associated with activated signalling pathways, notably JNK and AP-1 (ref. 2). The transformed phenotypes are reversed by treatment with the theilericidal drug buparvaquone. We used comparative genomics to identify a homologue of the peptidyl-prolyl isomerase PIN1 in T. annulata (TaPIN1) that is secreted into the host cell and modulates oncogenic signalling pathways. Here we show that TaPIN1 is a bona fide prolyl isomerase and that it interacts with the host ubiquitin ligase FBW7, leading to its degradation and subsequent stabilization of c-JUN, which promotes transformation. We performed in vitro and in silico analysis and in vivo zebrafish xenograft experiments to demonstrate that TaPIN1 is directly inhibited by the anti-parasite drug buparvaquone (and other known PIN1 inhibitors) and is mutated in a drug-resistant strain. Prolyl isomerization is thus a conserved mechanism that is important in cancer and is used by Theileria parasites to manipulate host oncogenic signalling.

Jumbled Genomes: Missing Apicomplexan Synteny

Whole-genome comparisons provide insight into genome evolution by informing on gene repertoires, gene gains/losses, and genome organization. Most of our knowledge about eukaryotic genome evolution is derived from studies of multicellular model organisms. The eukaryotic phylum Apicomplexa contains obligate intracellular protist parasites responsible for a wide range of human and veterinary diseases (e.g., malaria, toxoplasmosis, and theileriosis). We have developed an in silico protein-encoding gene based pipeline to investigate synteny across 12 apicomplexan species from six genera. Genome rearrangement between lineages is extensive. Syntenic regions (conserved gene content and order) are rare between lineages and appear to be totally absent across the phylum, with no group of three genes found on the same chromosome and in the same order within 25 kb up- and downstream of any orthologous genes. Conserved synteny between major lineages is limited to small regions in Plasmodium and Theileria/Babesia species, and within these conserved regions, there are a number of proteins putatively targeted to organelles. The observed overall lack of synteny is surprising considering the divergence times and the apparent absence of transposable elements (TEs) within any of the species examined. TEs are ubiquitous in all other groups of eukaryotes studied to date and have been shown to be involved in genomic rearrangements. It appears that there are different criteria governing genome evolution within the Apicomplexa relative to other well-studied unicellular and multicellular eukaryotes.

Real-Time Loop-Mediated Isothermal Amplification (RealAmp) for the Species-Specific Identification of Plasmodium vivax

Plasmodium vivax infections remain a major source of malaria-related morbidity and mortality. Early and accurate diagnosis is an integral component of effective malaria control programs. Conventional molecular diagnostic methods provide accurate results but are often resource-intensive, expensive, have a long turnaround time and are beyond the capacity of most malaria-endemic countries. Our laboratory has recently developed a new platform called RealAmp, which combines loop-mediated isothermal amplification (LAMP) with a portable tube scanner real-time isothermal instrument for the rapid detection of malaria parasites. Here we describe new primers for the detection of P. vivax using the RealAmp method. Three pairs of amplification primers required for this method were derived from a conserved DNA sequence unique to the P. vivax genome. The amplification was carried out at 64°C using SYBR Green or SYTO-9 intercalating dyes for 90 minutes with the tube scanner set to collect fluorescence signals at 1-minute intervals. Clinical samples of P. vivax and other human-infecting malaria parasite species were used to determine the sensitivity and specificity of the primers by comparing with an 18S ribosomal RNA-based nested PCR as the gold standard. The new set of primers consistently detected laboratory-maintained isolates of P. vivax from different parts of the world. The primers detected P. vivax in the clinical samples with 94.59% sensitivity (95% CI: 87.48–98.26%) and 100% specificity (95% CI: 90.40–100%) compared to the gold standard nested-PCR method. The new primers also proved to be more sensitive than the published species-specific primers specifically developed for the LAMP method in detecting P. vivax.

Applied genomics: Data mining reveals species-specific malaria diagnostic targets more sensitive than 18S rRNA

ABSTRACT Accurate and rapid diagnosis of malaria infections is crucial for implementing species-appropriate treatment and saving lives. Molecular diagnostic tools are the most accurate and sensitive method of detecting Plasmodium, differentiating between Plasmodium species, and detecting subclinical infections. Despite available whole-genome sequence data for Plasmodium falciparum and P. vivax, the majority of PCR-based methods still rely on the 18S rRNA gene targets. Historically, this gene has served as the best target for diagnostic assays. However, it is limited in its ability to detect mixed infections in multiplex assay platforms without the use of nested PCR. New diagnostic targets are needed. Ideal targets will be species specific, highly sensitive, and amenable to both single-step and multiplex PCRs. We have mined the genomes of P. falciparum and P. vivax to identify species-specific, repetitive sequences that serve as new PCR targets for the detection of malaria. We show that these targets (Pvr47 and Pfr364) exist in 14 to 41 copies and are more sensitive than 18S rRNA when utilized in a single-step PCR. Parasites are routinely detected at levels of 1 to 10 parasites/μl. The reaction can be multiplexed to detect both species in a single reaction. We have examined 7 P. falciparum strains and 91 P. falciparum clinical isolates from Tanzania and 10 P. vivax strains and 96 P. vivax clinical isolates from Venezuela, and we have verified a sensitivity and specificity of ∼100% for both targets compared with a nested 18S rRNA approach. We show that bioinformatics approaches can be successfully applied to identify novel diagnostic targets and improve molecular methods for pathogen detection. These novel targets provide a powerful alternative molecular diagnostic method for the detection of P. falciparum and P. vivax in conventional or multiplex PCR platforms.

A New Single-Step PCR Assay for the Detection of the Zoonotic Malaria Parasite Plasmodium knowlesi

Background Recent studies in Southeast Asia have demonstrated substantial zoonotic transmission of Plasmodium knowlesi to humans. Microscopically, P. knowlesi exhibits several stage-dependent morphological similarities to P. malariae and P. falciparum. These similarities often lead to misdiagnosis of P. knowlesi as either P. malariae or P. falciparum and PCR-based molecular diagnostic tests are required to accurately detect P. knowlesi in humans. The most commonly used PCR test has been found to give false positive results, especially with a proportion of P. vivax isolates. To address the need for more sensitive and specific diagnostic tests for the accurate diagnosis of P. knowlesi, we report development of a new single-step PCR assay that uses novel genomic targets to accurately detect this infection. Methodology and Significant Findings We have developed a bioinformatics approach to search the available malaria parasite genome database for the identification of suitable DNA sequences relevant for molecular diagnostic tests. Using this approach, we have identified multi-copy DNA sequences distributed in the P. knowlesi genome. We designed and tested several novel primers specific to new target sequences in a single-tube, non-nested PCR assay and identified one set of primers that accurately detects P. knowlesi. We show that this primer set has 100% specificity for the detection of P. knowlesi using three different strains (Nuri, H, and Hackeri), and one human case of malaria caused by P. knowlesi. This test did not show cross reactivity with any of the four human malaria parasite species including 11 different strains of P. vivax as well as 5 additional species of simian malaria parasites. Conclusions The new PCR assay based on novel P. knowlesi genomic sequence targets was able to accurately detect P. knowlesi. Additional laboratory and field-based testing of this assay will be necessary to further validate its utility for clinical diagnosis of P. knowlesi.

Discovery and assembly of repeat family pseudomolecules from sparse genomic sequence data using the Assisted Automated Assembler of Repeat Families (AAARF) algorithm

BackgroundHigher eukaryotic genomes are typically large, complex and filled with both genes and multiple classes of repetitive DNA. The repetitive DNAs, primarily transposable elements, are a rapidly evolving genome component that can provide the raw material for novel selected functions and also indicate the mechanisms and history of genome evolution in any ancestral lineage. Despite their abundance, universality and significance, studies of genomic repeat content have been largely limited to analyses of the repeats in fully sequenced genomes.ResultsIn order to facilitate a broader range of repeat analyses, the Assisted Automated Assembler of Repeat Families algorithm has been developed. This program, written in PERL and with numerous adjustable parameters, identifies sequence overlaps in small shotgun sequence datasets and walks them out to create long pseudomolecules representing the most abundant repeats in any genome. Testing of this program in maize indicated that it found and assembled all of the major repeats in one or more pseudomolecules, including coverage of the major Long Terminal Repeat retrotransposon families. Both Sanger sequence and 454 datasets were appropriate.ConclusionThese results now indicate that hundreds of higher eukaryotic genomes can be efficiently characterized for the nature, abundance and evolution of their major repetitive DNA components.