Use Of Comparative Genomics Research Articles

Use of Comparative Genomics To Resolve an Unusual Case of Aminoglycoside Susceptibility in the Melioidosis Pathogen Burkholderia pseudomallei in Bangladesh.

Melioidosis is an emerging tropical infectious disease with a rising global burden caused by the environmental bacterium Burkholderia pseudomallei. It is endemic in Southeast and South Asia, including Bangladesh. A rare aminoglycoside-susceptible B. pseudomallei isolate (Y2019) has recently been reported from a melioidosis patient in Dhaka, Bangladesh. To understand the geographical origins of Y2019, we subjected it and 10 other isolates from Bangladesh to whole-genome sequencing. In a phylogenetic tree with a global set of B. pseudomallei genomes, most Bangladeshi genomes clustered tightly within the Asian clade. In contrast, Y2019 was closely related to ST881 isolates from Sarawak, Malaysian Borneo, a gentamicin-sensitive sequence type, suggesting infection in Borneo. Y2019 also contained the same gentamicin sensitivity conferring nonsynonymous mutation in the drug efflux pump encoding the amrB gene. In the absence of a full travel history, whole-genome sequencing and bioinformatics tools have revealed the likely origin of this rare isolate.

A comparative transcriptomics and eQTL approach identifies SlWD40 as a tomato fruit ripening regulator.

Although multiple vital genes with strong effects on the tomato (Solanum lycopersicum) ripening process have been identified via the positional cloning of ripening mutants and cloning of ripening-related transcription factors (TFs), recent studies suggest that it is unlikely that we have fully characterized the gene regulatory networks underpinning this process. Here, combining comparative transcriptomics and expression QTLs, we identified 16 candidate genes involved in tomato fruit ripening and validated them through virus-induced gene silencing analysis. To further confirm the accuracy of the approach, one potential ripening regulator, SlWD40 (WD-40 repeats), was chosen for in-depth analysis. Co-expression network analysis indicated that master regulators such as RIN (ripening inhibitor) and NOR (nonripening) as well as vital TFs including FUL1 (FRUITFUL1), SlNAC4 (NAM, ATAF1,2, and CUC2 4), and AP2a (Activating enhancer binding Protein 2 alpha) strongly co-expressed with SlWD40. Furthermore, SlWD40 overexpression and RNAi lines exhibited substantially accelerated and delayed ripening phenotypes compared with the wild type, respectively. Moreover, transcriptome analysis of these transgenics revealed that expression patterns of ethylene biosynthesis genes, phytoene synthase, pectate lyase, and branched chain amino transferase 2, in SlWD40-RNAi lines were similar to those of rin and nor fruits, which further demonstrated that SlWD40 may act as an important ripening regulator in conjunction with RIN and NOR. These results are discussed in the context of current models of ripening and in terms of the use of comparative genomics and transcriptomics as an effective route for isolating causal genes underlying differences in genotypes.

Functional pangenome analysis reveals high virulence plasticity of Aliarcobacter butzleri and affinity to human mucus

Aliarcobacter butzleri is an emerging pathogen that may cause enteritis in humans, however, the incidence of disease caused by this member of the Campylobacteriaceae family is still underestimated. Furthermore, little is known about the precise virulence mechanism and behavior during infection. Therefore, in the present study, through complementary use of comparative genomics and physiological tests on human gut models, we sought to elucidate the genetic background of a set of 32 A. butzleri strains of diverse origin and to explore the correlation with the ability to colonize and invade human intestinal cells in vitro.The simulated infection of human intestinal models showed a higher colonization rate in presence of mucus-producing cells. For some strains, human mucus significantly improved the resistance to physical removal from the in vitro mucosa, while short time-frame growth was even observed. Pangenome analysis highlighted a hypervariable accessory genome, not strictly correlated to the isolation source. Likewise, the strain phylogeny was unrelated to their shared origin, despite a certain degree of segregation was observed among strains isolated from different segments of the intestinal tract of pigs. The putative virulence genes detected in all strains were mostly encompassed in the accessory fraction of the pangenome. The LPS biosynthesis and in particular the chain glycosylation of the O-antigen is harbored in a region of high plasticity of the pangenome, which would indicate frequent horizontal gene transfer phenomena, as well as the involvement of this hypervariable structure in the adaptive behavior and sympatric evolution of A. butzleri.Results of the present study deepen the current knowledge on A. butzleri pangenome by extending the pool of genes regarded as virulence markers and provide bases to develop new diagnostic approaches for the detection of those strains with a higher virulence potential.

Comparative genomics of the gut commensal Bifidobacterium bifidum reveals adaptation to carbohydrate utilization

Bifidobacterium bifidum is one of the most abundant members of the gut microbiota at the early stage of life. The established association of the bacterium with the human gut confers health benefits. Such successful persistence of B. bifidum necessitates metabolic adaptation to the host–derived carbohydrates, a process which is poorly understood. The current study focuses on revealing the genomic-based phylogeny (phylogenomics) of B. bifidum and utilizing comparative genomics to decipher the glycolytic abilities of bifidobacterial strains isolated from different human body niches (feces, human gut, vagina, and breast milk). When the phylogenomic analysis was performed on 95 B. bifidum strains, currently available on the RefSeq database, the bacterium was clearly distinguished from other members of the Bifidobacterium genus. Furthermore, a pairwise genomic comparison indicated that a large proportion of orthologous gene families were shared among the B. bifidum strains. These findings highlight the notion that the B. bifidum species is genetically similar and may perform similar functions in their host. When 15 B. bifidum genomes representing strains from different human body niches were annotated, the resulting functional profile showed the presence of enriched proteins involved in carbohydrate utilization. Moreover, mining the 15 B. bifidum genomes for the presence of Carbohydrate–Active Enzyme (CAZY) systems, the analysis found the existence of diverse protein families which include glycosyl hydrolases, glycosyl transferases, carbohydrate-binding modules, and carbohydrate esterases. Collectively, these CAZY systems enables B. bifidum to utilize host–derived glycans (e.g., mucin) and diet–derived carbohydrates (e.g., starch). In contrast, a correlation analysis revealed that B. bifidum strains isolated from the different body niches were indistinguishable in the context of presence–absence of CAZY systems. These findings emphasize the valuable use of comparative genomics in deciphering the glycolytic abilities of B. bifidum and consequently its adaptation to carbohydrate utilization in the human gut environment.

Chromosome-Scale Assembly and Annotation of the Macadamia Genome (Macadamia integrifolia HAES 741).

Macadamia integrifolia is a representative of the large basal eudicot family Proteaceae and the main progenitor species of the Australian native nut crop macadamia. Since its commercialisation in Hawaii fewer than 100 years ago, global production has expanded rapidly. However, genomic resources are limited in comparison to other horticultural crops. The first draft assembly of M. integrifolia had good coverage of the functional gene space but its high fragmentation has restricted its use in comparative genomics and association studies. Here we have generated an improved assembly of cultivar HAES 741 (4,094 scaffolds, 745 Mb, N50 413 kb) using a combination of Illumina paired and PacBio long read sequences. Scaffolds were anchored to 14 pseudo-chromosomes using seven genetic linkage maps. This assembly has improved contiguity and coverage, with >120 Gb of additional sequence. Following annotation, 34,274 protein-coding genes were predicted, representing 90% of the expected gene content. Our results indicate that the macadamia genome is repetitive and heterozygous. The total repeat content was 55% and genome-wide heterozygosity, estimated by read mapping, was 0.98% or an average of one SNP per 102 bp. This is the first chromosome-scale genome assembly for macadamia and the Proteaceae. It is expected to be a valuable resource for breeding, gene discovery, conservation and evolutionary genomics.

PPanGGOLiN: Depicting microbial diversity via a partitioned pangenome graph.

The use of comparative genomics for functional, evolutionary, and epidemiological studies requires methods to classify gene families in terms of occurrence in a given species. These methods usually lack multivariate statistical models to infer the partitions and the optimal number of classes and don't account for genome organization. We introduce a graph structure to model pangenomes in which nodes represent gene families and edges represent genomic neighborhood. Our method, named PPanGGOLiN, partitions nodes using an Expectation-Maximization algorithm based on multivariate Bernoulli Mixture Model coupled with a Markov Random Field. This approach takes into account the topology of the graph and the presence/absence of genes in pangenomes to classify gene families into persistent, cloud, and one or several shell partitions. By analyzing the partitioned pangenome graphs of isolate genomes from 439 species and metagenome-assembled genomes from 78 species, we demonstrate that our method is effective in estimating the persistent genome. Interestingly, it shows that the shell genome is a key element to understand genome dynamics, presumably because it reflects how genes present at intermediate frequencies drive adaptation of species, and its proportion in genomes is independent of genome size. The graph-based approach proposed by PPanGGOLiN is useful to depict the overall genomic diversity of thousands of strains in a compact structure and provides an effective basis for very large scale comparative genomics. The software is freely available at https://github.com/labgem/PPanGGOLiN.

Comparative Genomics to Develop a Specific Multiplex PCR Assay for Detection of Clavibacter michiganensis.

Clavibacter michiganensis is a Gram-positive bacterial pathogen that proliferates in the xylem vessels of tomato, causing bacterial wilt and canker symptoms. Accurate detection is a crucial step in confirming outbreaks of bacterial canker and developing management strategies. A major problem with existing detection methods are false-positive and -negative results. Here, we report the use of comparative genomics of 37 diverse Clavibacter strains, including 21 strains sequenced in this study, to identify specific sequences that are C. michiganensis detection targets. Genome-wide phylogenic analyses revealed additional diversity within the genus Clavibacter. Pathogenic C. michiganensis strains varied in plasmid composition, highlighting the need for detection methods based on chromosomal targets. We utilized sequences of C. michiganensis-specific loci to develop a multiplex PCR-based diagnostic platform using two C. michiganensis chromosomal genes (rhuM and tomA) and an internal control amplifying both bacterial and plant DNA (16s ribosomal RNA). The multiplex PCR assay specifically detected C. michiganensis strains from a panel of 110 additional bacteria, including other Clavibacter spp. and bacterial pathogens of tomato. The assay was adapted to detect the presence of C. michiganensis in seed and tomato plant materials with high sensitivity and specificity. In conclusion, the described method represents a robust, specific tool for detection of C. michiganensis in tomato seed and infected plants.

A chromosome-scale genome assembly of cucumber (Cucumis sativus L.).

BackgroundAccurate and complete reference genome assemblies are fundamental for biological research. Cucumber is an important vegetable crop and model system for sex determination and vascular biology. Low-coverage Sanger sequences and high-coverage short Illumina sequences have been used to assemble draft cucumber genomes, but the incompleteness and low quality of these genomes limit their use in comparative genomics and genetic research. A high-quality and complete cucumber genome assembly is therefore essential.FindingsWe assembled single-molecule real-time (SMRT) long reads to generate an improved cucumber reference genome. This version contains 174 contigs with a total length of 226.2 Mb and an N50 of 8.9 Mb, and provides 29.0 Mb more sequence data than previous versions. Using 10X Genomics and high-throughput chromosome conformation capture (Hi-C) data, 89 contigs (∼211.0 Mb) were directly linked into 7 pseudo-chromosome sequences. The newly assembled regions show much higher guanine-cytosine or adenine-thymine content than found previously, which is likely to have been inaccessible to Illumina sequencing. The new assembly contains 1,374 full-length long terminal retrotransposons and 1,078 novel genes including 239 tandemly duplicated genes. For example, we found 4 tandemly duplicated tyrosylprotein sulfotransferases, in contrast to the single copy of the gene found previously and in most other plants.ConclusionThis high-quality genome presents novel features of the cucumber genome and will serve as a valuable resource for genetic research in cucumber and plant comparative genomics.

Detecting Adaptation with Genome-Scale Molecule Evolutionary Analysis: An Educational Primer for Use with "RNA Interference Pathways Display High Rates of Adaptive Protein Evolution in Multiple Invertebrates".

Hosts and pathogens impose coevolutionary pressure on each other as pathogens strive to establish themselves and hosts seek to suppress infection. RNA interference (RNAi) is a mechanism by which cells repress viruses and transposable elements, thereby serving as a form of immune defense. Previous studies have shown that antiviral RNAi genes evolve extraordinarily quickly in the fruit fly Drosophila melanogaster, suggesting that they may adaptively coevolve with viruses and transposable elements. An article by Palmer and colleagues extends this observation to nematodes and multiple insects. Their article can be combined with this Primer to demonstrate the use of comparative genomics and molecular evolutionary analyses in the measurement of natural selection.Related article in GENETICS: Palmer, W. H., J. D. Hadfield, and D. J. Obbard, 2018 RNA-Interference pathways display high rates of adaptive protein evolution in multiple invertebrates. Genetics 208: 1585–1599.

Development and application of genomic resources for comparative and translational genomics in legumes through leveraging genomic sequence of Medicago truncatula

The expressed sequence tags (ESTs) of common bean were BLAST aligned with barred medic genome sequence and developed 1196 conserved intron spanning primers (CISPs) to facilitate genetic studies in legumes. Randomly selected 288 CISPs, representing loci on barrel medic genome, were tested on 10 selected members of legume family. On the source taxa, the highest single copy amplification success rates of 61.8% (barrel medic) and 56.2% (common bean) was obtained. The success rate of markers was 54.5% in cowpea followed by 53.5% in pigeonpea and chickpea, signifying cross taxon amplification and their potential use in comparative genomics. However, relatively low percentages of primer set amplified (40-43%) in soybean, urdbean and peanut. Further, these primers were tested on different varieties of chickpea, pigeonpea and cowpea. The PCR products were sequenced and aligned which resulted in detection of 26 SNPs and eight INDeLs in cowpea, seven SNPs and two INDeLs in chickpea and 27 SNPs and 14 INDeLs in pigeonpea. These SNPs were successfully converted in to size variation for gel-based genotyping. The CISP markers developed in this study are expected to aid in map saturation of legumes and in marker-assisted selection for accelerated crop improvement.

Indoleacetate decarboxylase is a glycyl radical enzyme catalysing the formation of malodorant skatole

Skatole is a malodorous compound that contributes to the characteristic smell of animal faeces. Although skatole has long been known to originate from bacterial tryptophan fermentation, the enzyme catalysing its formation has so far remained elusive. Here we report the use of comparative genomics for the discovery of indoleacetate decarboxylase, an O2-sensitive glycyl radical enzyme catalysing the decarboxylation of indoleacetate to form skatole as the terminal step of tryptophan fermentation in certain anaerobic bacteria. We describe its biochemical characterization and compare it to other glycyl radical decarboxylases. Indoleacetate decarboxylase may serve as a genetic marker for the identification of skatole-producing environmental and human-associated bacteria, with impacts on human health and the livestock industry.

Varied Genomic Responses to Maladaptive Gene Flow and Their Evidence.

Adaptation to a local environment often occurs in the face of maladaptive gene flow. In this perspective, I discuss several ideas on how a genome may respond to maladaptive gene flow during adaptation. On the one hand, selection can build clusters of locally adaptive alleles at fortuitously co-localized loci within a genome, thereby facilitating local adaptation with gene flow (‘allele-only clustering’). On the other hand, the selective pressure to link adaptive alleles may drive co-localization of the actual loci relevant for local adaptation within a genome through structural genome changes or an evolving intra-genomic crossover rate (‘locus clustering’). While the expected outcome is, in both cases, a higher frequency of locally adaptive alleles in some genome regions than others, the molecular units evolving in response to gene flow differ (i.e., alleles versus loci). I argue that, although making this distinction is important, we commonly lack the critical empirical evidence to do so. This is mainly because many current approaches are biased towards detecting local adaptation in genome regions with low crossover rates. The importance of low-crossover genome regions for adaptation with gene flow, such as in co-localizing relevant loci within a genome, thus remains unclear. Future empirical investigations should address these questions by making use of comparative genomics, where multiple de novo genome assemblies from species evolved under different degrees of genetic exchange are compared. This research promises to advance our understanding of how a genome adapts to maladaptive gene flow, thereby promoting adaptive divergence and reproductive isolation.

New Insights into the Diversity of the Genus Faecalibacterium.

Faecalibacterium prausnitzii is a commensal bacterium, ubiquitous in the gastrointestinal tracts of animals and humans. This species is a functionally important member of the microbiota and studies suggest it has an impact on the physiology and health of the host. F. prausnitzii is the only identified species in the genus Faecalibacterium, but a recent study clustered strains of this species in two different phylogroups. Here, we propose the existence of distinct species in this genus through the use of comparative genomics. Briefly, we performed analyses of 16S rRNA gene phylogeny, phylogenomics, whole genome Multi-Locus Sequence Typing (wgMLST), Average Nucleotide Identity (ANI), gene synteny, and pangenome to better elucidate the phylogenetic relationships among strains of Faecalibacterium. For this, we used 12 newly sequenced, assembled, and curated genomes of F. prausnitzii, which were isolated from feces of healthy volunteers from France and Australia, and combined these with published data from 5 strains downloaded from public databases. The phylogenetic analysis of the 16S rRNA sequences, together with the wgMLST profiles and a phylogenomic tree based on comparisons of genome similarity, all supported the clustering of Faecalibacterium strains in different genospecies. Additionally, the global analysis of gene synteny among all strains showed a highly fragmented profile, whereas the intra-cluster analyses revealed larger and more conserved collinear blocks. Finally, ANI analysis substantiated the presence of three distinct clusters—A, B, and C—composed of five, four, and four strains, respectively. The pangenome analysis of each cluster corroborated the classification of these clusters into three distinct species, each containing less variability than that found within the global pangenome of all strains. Here, we propose that comparison of pangenome subsets and their associated α values may be used as an alternative approach, together with ANI, in the in silico classification of new species. Altogether, our results provide evidence not only for the reconsideration of the phylogenetic and genomic relatedness among strains currently assigned to F. prausnitzii, but also the need for lineage (strain-based) differentiation of this taxon to better define how specific members might be associated with positive or negative host interactions.

Phylogenomic, Pan-genomic, Pathogenomic and Evolutionary Genomic Insights into the Agronomically Relevant Enterobacteria Pantoea ananatis and Pantoea stewartii.

Pantoea ananatis is ubiquitously found in the environment and causes disease on a wide range of plant hosts. By contrast, its sister species, Pantoea stewartii subsp. stewartii is the host-specific causative agent of the devastating maize disease Stewart’s wilt. This pathogen has a restricted lifecycle, overwintering in an insect vector before being introduced into susceptible maize cultivars, causing disease and returning to overwinter in its vector. The other subspecies of P. stewartii subsp. indologenes, has been isolated from different plant hosts and is predicted to proliferate in different environmental niches. Here we have, by the use of comparative genomics and a comprehensive suite of bioinformatic tools, analyzed the genomes of ten P. stewartii and nineteen P. ananatis strains. Our phylogenomic analyses have revealed that there are two distinct clades within P. ananatis while far less phylogenetic diversity was observed among the P. stewartii subspecies. Pan-genome analyses revealed a large core genome comprising of 3,571 protein coding sequences is shared among the twenty-nine compared strains. Furthermore, we showed that an extensive accessory genome made up largely by a mobilome of plasmids, integrated prophages, integrative and conjugative elements and insertion elements has resulted in extensive diversification of P. stewartii and P. ananatis. While these organisms share many pathogenicity determinants, our comparative genomic analyses show that they differ in terms of the secretion systems they encode. The genomic differences identified in this study have allowed us to postulate on the divergent evolutionary histories of the analyzed P. ananatis and P. stewartii strains and on the molecular basis underlying their ecological success and host range.

Use of Comparative Genomics To Characterize the Diversity of Acinetobacter baumannii Surveillance Isolates in a Health Care Institution

Despite the increasing prevalence of the nosocomial pathogen Acinetobacter baumannii, little is known about which genomic components contribute to clinical presentation of this important pathogen. Most whole-genome comparisons of A. baumannii have focused on specific genomic regions associated with phenotypes in a limited number of genomes. In this work, we describe the results of a whole-genome comparative analysis of 254 surveillance isolates of Acinetobacter species, 203 of which were A. baumannii, isolated from perianal swabs and sputum samples collected as part of an infection control active surveillance program at the University of Maryland Medical Center. The collection of surveillance isolates includes both carbapenem-susceptible and -resistant isolates. Based on the whole-genome phylogeny, the A. baumannii isolates collected belong to two major phylogenomic lineages. Results from multilocus sequence typing indicated that one of the major phylogenetic groups of A. baumannii was comprised solely of strains from the international clonal lineage 2. The genomic content of the A. baumannii isolates was examined using large-scale BLAST score ratio analysis to identify genes that are associated with carbapenem-susceptible and -resistant isolates, as well as genes potentially associated with the source of isolation. This analysis revealed a number of genes that were exclusive or at greater frequency in each of these classifications. This study is the most comprehensive genomic comparison of Acinetobacter isolates from a surveillance study to date and provides important information that will contribute to our understanding of the success of A. baumannii as a human pathogen.

Paradoxical Hypersusceptibility of Drug-resistant Mycobacteriumtuberculosis to β-lactam Antibiotics

Mycobacterium tuberculosis (M. tuberculosis) is considered innately resistant to β-lactam antibiotics. However, there is evidence that susceptibility to β-lactam antibiotics in combination with β–lactamase inhibitors is variable among clinical isolates, and these may present therapeutic options for drug-resistant cases. Here we report our investigation of susceptibility to β-lactam/β–lactamase inhibitor combinations among clinical isolates of M. tuberculosis, and the use of comparative genomics to understand the observed heterogeneity in susceptibility. Eighty-nine South African clinical isolates of varying first and second-line drug susceptibility patterns and two reference strains of M. tuberculosis underwent minimum inhibitory concentration (MIC) determination to two β-lactams: amoxicillin and meropenem, both alone and in combination with clavulanate, a β–lactamase inhibitor. 41/91 (45%) of tested isolates were found to be hypersusceptible to amoxicillin/clavulanate relative to reference strains, including 14/24 (58%) of multiple drug-resistant (MDR) and 22/38 (58%) of extensively drug-resistant (XDR) isolates. Genome-wide polymorphisms identified using whole-genome sequencing were used in a phylogenetically-aware linear mixed model to identify polymorphisms associated with amoxicillin/clavulanate susceptibility. Susceptibility to amoxicillin/clavulanate was over-represented among isolates within a specific clade (LAM4), in particular among XDR strains. Twelve sets of polymorphisms were identified as putative markers of amoxicillin/clavulanate susceptibility, five of which were confined solely to LAM4. Within the LAM4 clade, ‘paradoxical hypersusceptibility’ to amoxicillin/clavulanate has evolved in parallel to first and second-line drug resistance. Given the high prevalence of LAM4 among XDR TB in South Africa, our data support an expanded role for β-lactam/β-lactamase inhibitor combinations for treatment of drug-resistant M. tuberculosis.

Linking Biomarker and Comparative Omics to Pathogens in Legumes

It is envisioned that a more precise study of the association between the traits and biomarkers will dramatically decrease the time and costs required to bring new improved disease resistance lines to market. The field of omics has an enormous potential to assess diseases more precise, including the identification and understanding of pathogenic mechanisms in legume crops, and have been exemplified by a relatively large number of studies. Recently, molecular genetic studies have accumulated a huge amount of genotypic data, through a more affordable next generation sequencing (NGS) technology, causing the omics approaches to fall behind. In this paper I provide an overview of genomics and proteomics and their use in legume crops, including the use of comparative genomics to identify homologous markers within legume crops.

Five decades of invasion genetics.

Five decades of invasion genetics.

A Novel Candidate Vaccine for Cytauxzoonosis Inferred from Comparative Apicomplexan Genomics

Cytauxzoonosis is an emerging infectious disease of domestic cats (Felis catus) caused by the apicomplexan protozoan parasite Cytauxzoon felis. The growing epidemic, with its high morbidity and mortality points to the need for a protective vaccine against cytauxzoonosis. Unfortunately, the causative agent has yet to be cultured continuously in vitro, rendering traditional vaccine development approaches beyond reach. Here we report the use of comparative genomics to computationally and experimentally interpret the C. felis genome to identify a novel candidate vaccine antigen for cytauxzoonosis. As a starting point we sequenced, assembled, and annotated the C. felis genome and the proteins it encodes. Whole genome alignment revealed considerable conserved synteny with other apicomplexans. In particular, alignments with the bovine parasite Theileria parva revealed that a C. felis gene, cf76, is syntenic to p67 (the leading vaccine candidate for bovine theileriosis), despite a lack of significant sequence similarity. Recombinant subdomains of cf76 were challenged with survivor-cat antiserum and found to be highly seroreactive. Comparison of eleven geographically diverse samples from the south-central and southeastern USA demonstrated 91–100% amino acid sequence identity across cf76, including a high level of conservation in an immunogenic 226 amino acid (24 kDa) carboxyl terminal domain. Using in situ hybridization, transcription of cf76 was documented in the schizogenous stage of parasite replication, the life stage that is believed to be the most important for development of a protective immune response. Collectively, these data point to identification of the first potential vaccine candidate antigen for cytauxzoonosis. Further, our bioinformatic approach emphasizes the use of comparative genomics as an accelerated path to developing vaccines against experimentally intractable pathogens.

Understanding the Dynamics of Gene Regulatory Systems; Characterisation and Clinical Relevance of cis-Regulatory Polymorphisms.

Modern genetic analysis has shown that most polymorphisms associated with human disease are non-coding. Much of the functional information contained in the non-coding genome consists of cis-regulatory sequences (CRSs) that are required to respond to signal transduction cues that direct cell specific gene expression. It has been hypothesised that many diseases may be due to polymorphisms within CRSs that alter their responses to signal transduction cues. However, identification of CRSs, and the effects of allelic variation on their ability to respond to signal transduction cues, is still at an early stage. In the current review we describe the use of comparative genomics and experimental techniques that allow for the identification of CRSs building on recent advances by the ENCODE consortium. In addition we describe techniques that allow for the analysis of the effects of allelic variation and epigenetic modification on CRS responses to signal transduction cues. Using specific examples we show that the interactions driving these elements are highly complex and the effects of disease associated polymorphisms often subtle. It is clear that gaining an understanding of the functions of CRSs, and how they are affected by SNPs and epigenetic modification, is essential to understanding the genetic basis of human disease and stratification whilst providing novel directions for the development of personalised medicine.