Evolutionary Genomics and Conservation of the Endangered Przewalski's Horse.

Przewalski's horses (PHs, Equus ferus ssp. przewalskii) were discovered in the Asian steppes in the 1870s and represent the last remaining true wild horses. PHs became extinct in the wild in the 1960s but survived in captivity, thanks to major conservation efforts. The current population is still endangered, with just 2,109 individuals, one-quarter of which are in Chinese and Mongolian reintroduction reserves [1]. These horses descend from a founding population of 12 wild-caught PHs and possibly up to four domesticated individuals [2-4]. With a stocky build, an erect mane, and stripped and short legs, they are phenotypically and behaviorally distinct from domesticated horses (DHs, Equus caballus). Here, we sequenced the complete genomes of 11 PHs, representing all founding lineages, and five historical specimens dated to 1878-1929 CE, including the Holotype. These were compared to the hitherto-most-extensive genome dataset characterized for horses, comprising 21 new genomes. We found that loci showing the most genetic differentiation with DHs were enriched in genes involved in metabolism, cardiac disorders, muscle contraction, reproduction, behavior, and signaling pathways. We also show that DH and PH populations split ∼45,000 years ago and have remained connected by gene-flow thereafter. Finally, we monitor the genomic impact of ∼110 years of captivity, revealing reduced heterozygosity, increased inbreeding, and variable introgression of domestic alleles, ranging from non-detectable to as much as 31.1%. This, together with the identification of ancestry informative markers and corrections to the International Studbook, establishes a framework for evaluating the persistence of genetic variation in future reintroduced populations.

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.

A short critical history of the application of genomics to animal breeding

Two scientific schools have been in coexistence from the beginning of genetics, one of them searching for factors of inheritance and the other one applying biometrical models to study the relationships between relatives. With the development of molecular genetics, the possibilities of detecting genes having a noticeable effect in traits augmented. Some genes with large or medium effects were localized in animals, although the most common result was to detect markers linked to these genes, allowing the possibility of assisting selection programs with markers. When a large amount of simple and inexpensive markers were available, the SNPs, new possibilities were opened since they did not need the presence of genes of large or medium effect controlling a trait, because the whole genome was scanned. Using a large amount of SNPs permits having a prediction of the breeding value at birth accurate enough to be used in some cases, like dairy cattle, to halve its generation interval. In other animal breeding programs, the implementation of genomic selection is less clear and the way in which it can be useful should be carefully studied. The need for large populations for associating phenotypic data and markers, plus the need for repeating the process continuously, complicates its application in some cases. The implementation of the information provided by the SNPs in current genetic programs has led to the development of complex statistical tools, joining the efforts of the two schools, factorial and biometrical, that nowadays work closely related.

Structuring research methods and data with the research object model: genomics workflows as a case study

Background: One of the main challenges for biomedical research lies in the computer-assisted integrative study of large and increasingly complex combinations of data in order to understand molecular mechanisms. The preservation of the materials and methods of such computational experiments with clear annotations is essential for understanding an experiment, and this is increasingly recognized in the bioinformatics community. Our assumption is that offering means of digital, structured aggregation and annotation of the objects of an experiment will provide necessary meta-data for a scientist to understand and recreate the results of an experiment. To support this we explored a model for the semantic description of a workflow-centric Research Object (RO), where an RO is defined as a resource that aggregates other resources, e.g., datasets, software, spreadsheets, text, etc. We applied this model to a case study where we analysed human metabolite variation by workflows. Results: We present the application of the workflow-centric RO model for our bioinformatics case study. Three workflows were produced following recently defined Best Practices for workflow design. By modelling the experiment as an RO, we were able to automatically query the experiment and answer questions such as “which particular data was input to a particular workflow to test a particular hypothesis?”, and “which particular conclusions were drawn from a particular workflow?”. Conclusions: Applying a workflow-centric RO model to aggregate and annotate the resources used in a bioinformatics experiment, allowed us to retrieve the conclusions of the experiment in the context of the driving hypothesis, the executed workflows and their input data. The RO model is an extendable reference model that can be used by other systems as well.

Functional genomics tools to decipher the pathogenicity mechanisms of the necrotrophic fungus Plectosphaerella cucumerina in Arabidopsis thaliana

The analysis of the interaction between Arabidopsis thaliana and adapted (PcBMM) and nonadapted (Pc2127) isolates of the necrotrophic fungus Plectosphaerella cucumerina has contributed to the identification of molecular mechanisms controlling plant resistance to necrotrophs.To characterize the pathogenicity bases of the virulence of necrotrophic fungi in Arabidopsis, we developed P. cucumerina functional genomics tools using Agrobacterium tumefaciens-mediated transformation.We generated PcBMM-GFP and Pc2127-GFP transformants constitutively expressing the green fluorescence protein (GFP), and a collection of random T-DNA insertional PcBMM transformants. Confocal microscopy analyses of the initial stages of PcBMM-GFP infection revealed that this pathogen, like other necrotrophic fungi, does not form an appressorium or penetrate into plant cells, but causes successive degradation of leaf cell layers

Populations genomics analysis of legume host preference for specific rhizobial genotypes in the Rhizobium leguminosarum bv. viciae symbioses

Rhizobium leguminosarum bv. viciae establishes root nodule symbioses with several legume genera. Although most isolates are equally effective in establishing symbioses with all host genera, previous evidence suggests that hosts select specific rhizobial genotypes among those present in the soil. We have used population genomics to further investigate this observation. P. sativum, L. culinaris, V. sativa, and V. faba plants were used to trap rhizobia from a well-characterized soil, and pooled genomic DNAs from one-hundred isolates from each plant were sequenced. Sequence reads were aligned to the R. leguminosarum bv. viciae 3841 reference genome. High overall conservation of sequences was observed in all subpopulations, although several multigenic regions were absent from the soil population. A large fraction (16-22%) of sequence reads could not be recruited to the reference genome, suggesting that they represent sequences specific to that particular soil population. Although highly conserved, the 16S-23S rRNA gene region presented single nucleotide polymorphisms (SNPs) regarding the reference genome, but no striking differences could be found among plant-selected subpopulations. Plant-specific SNP patterns were, however, clearly observed within the nod gene cluster, supporting the existence of a plant preference for specific rhizobial genotypes. This was also shown after genome-wide analysis of SNP patterns.

Genomics of host specificity in the Rhizobium-legume simbiosis

Most Rhizobium leguminosarum bv. viciae isolates are able to specifically nodulate plants of any of four different legume genera: Pisum, Lens, Vicia, and Lathyrus. However, previous evidence suggests that some genotypes are more adapted to a given plant host than others, and that the plant host can select specific genotypes among those present in a given soil population. We have used a population genomics approach to confirm that this is indeed the case, and to analyze the specific genotypic characteristics that each plant host selects

Population genomics of host specificity in Rhizobium leguminosarum bv. viciae

Legumes establish a root-nodule symbiosis with soil bacteria collectively known as rhizobia. This symbiosis allows legumes to benefit from the nitrogen fixation capabilities of rhizobia and thus to grow in the absence of any fixed nitrogen source. This is especially relevant for Agriculture, where intensive plant growth depletes soils of useable, fixed nitrogen sources. One of the main features of the root nodule symbiosis is its specificity. Different rhizobia are able to nodulate different legumes. Rhizobium leguminosarum bv. viciae is able to establish an effective symbiosis with four different plant genera (Pisum, Lens, Vicia, Lathyrus), and any given isolate will nodulate any of the four plant genera. A population genomics study with rhizobia isolated from P. sativum, L. culinaris, V. sativa or V. faba, all originating in the same soil, showed that plants select specific genotypes from those available in that soil. This was demonstrated at the genome-wide level, but also for specific genes. Accelerated mesocosm studies with successive plant cultures provided additional evidence on this plant selection and on the nature of the genotypes selected. Finally, representatives from the major rhizobial genotypes isolated from these plants allowed characterization of the size and nature of the respective pangenome and specific genome compartments. These were compared to the different genotypes ?symbiotic and non-symbiotic?present in rhizobial populations isolated directly from the soil without plant intervention.

Hox cluster genomics in the horn shark, Heterodontus francisci

Reconstructing the evolutionary history of Hox cluster origins will lead to insights into the developmental and evolutionary significance of Hox gene clusters in vertebrate phylogeny and to their role in the origins of various vertebrate body plans. We have isolated two Hox clusters from the horn shark, Heterodontus francisci. These have been sequenced and compared with one another and with other chordate Hox clusters. The results show that one of the horn shark clusters (HoxM) is orthologous to the mammalian HoxA cluster and shows a structural similarity to the amphioxus cluster, whereas the other shark cluster (HoxN) is orthologous to the mammalian HoxD cluster based on cluster organization and a comparison with noncoding and Hox gene-coding sequences. The persistence of an identifiable HoxA cluster over an 800-million-year divergence time demonstrates that the Hox gene clusters are highly integrated and structured genetic entities. The data presented herein identify many noncoding sequence motifs conserved over 800 million years that may function as genetic control motifs essential to the developmental process.

Toward functional genomics in bacteria: Analysis of gene expression in Escherichia coli from a bacterial artificial chromosome library of Bacillus cereus

As the study of microbes moves into the era of functional genomics, there is an increasing need for molecular tools for analysis of a wide diversity of microorganisms. Currently, biological study of many prokaryotes of agricultural, medical, and fundamental scientific interest is limited by the lack of adequate genetic tools. We report the application of the bacterial artificial chromosome (BAC) vector to prokaryotic biology as a powerful approach to address this need. We constructed a BAC library in Escherichia coli from genomic DNA of the Gram-positive bacterium Bacillus cereus. This library provides 5.75-fold coverage of the B. cereus genome, with an average insert size of 98 kb. To determine the extent of heterologous expression of B. cereus genes in the library, we screened it for expression of several B. cereus activities in the E. coli host. Clones expressing 6 of 10 activities tested were identified in the library, namely, ampicillin resistance, zwittermicin A resistance, esculin hydrolysis, hemolysis, orange pigment production, and lecithinase activity. We analyzed selected BAC clones genetically to identify rapidly specific B. cereus loci. These results suggest that BAC libraries will provide a powerful approach for studying gene expression from diverse prokaryotes.

Identification of eIF2Bγ and eIF2γ as cofactors of hepatitis C virus internal ribosome entry site-mediated translation using a functional genomics approach

The 5′-untranslated region of hepatitis C virus (HCV) is highly conserved, folds into a complex secondary structure, and functions as an internal ribosome entry site (IRES) to initiate translation of HCV proteins. We have developed a selection system based on a randomized hairpin ribozyme gene library to identify cellular factors involved in HCV IRES function. A retroviral vector ribozyme library with randomized target recognition sequences was introduced into HeLa cells, stably expressing a bicistronic construct encoding the hygromycin B phosphotransferase gene and the herpes simplex virus thymidine kinase gene (HSV-tk). Translation of the HSV-tk gene was mediated by the HCV IRES. Cells expressing ribozymes that inhibit HCV IRES-mediated translation of HSV-tk were selected via their resistance to both ganciclovir and hygromycin B. Two ribozymes reproducibly conferred the ganciclovir-resistant phenotype and were shown to inhibit IRES-mediated translation of HCV core protein but did not inhibit cap-dependent protein translation or cell growth. The functional targets of these ribozymes were identified as the gamma subunits of human eukaryotic initiation factors 2B (eIF2Bγ) and 2 (eIF2γ), respectively. The involvement of eIF2Bγ and eIF2γ in HCV IRES-mediated translation was further validated by ribozymes directed against additional sites within the mRNAs of these genes. In addition to leading to the identification of cellular IRES cofactors, ribozymes obtained from this cellular selection system could be directly used to specifically inhibit HCV viral translation, thereby facilitating the development of new antiviral strategies for HCV infection.

A chemical genomics approach toward understanding the global functions of the target of rapamycin protein (TOR)

The target of rapamycin protein (TOR) is a highly conserved ataxia telangiectasia-related protein kinase essential for cell growth. Emerging evidence indicates that TOR signaling is highly complex and is involved in a variety of cellular processes. To understand its general functions, we took a chemical genomics approach to explore the genetic interaction between TOR and other yeast genes on a genomic scale. In this study, the rapamycin sensitivity of individual deletion mutants generated by the Saccharomyces Genome Deletion Project was systematically measured. Our results provide a global view of the rapamycin-sensitive functions of TOR. In contrast to conventional genetic analysis, this approach offers a simple and thorough analysis of genetic interaction on a genomic scale and measures genetic interaction at different possible levels. It can be used to study the functions of other drug targets and to identify novel protein components of a conserved core biological process such as DNA damage checkpoint/repair that is interfered with by a cell-permeable chemical compound.

Structure-based assignment of the biochemical function of a hypothetical protein: A test case of structural genomics

Many small bacterial, archaebacterial, and eukaryotic genomes have been sequenced, and the larger eukaryotic genomes are predicted to be completely sequenced within the next decade. In all genomes sequenced to date, a large portion of these organisms’ predicted protein coding regions encode polypeptides of unknown biochemical, biophysical, and/or cellular functions. Three-dimensional structures of these proteins may suggest biochemical or biophysical functions. Here we report the crystal structure of one such protein, MJ0577, from a hyperthermophile, Methanococcus jannaschii, at 1.7-Å resolution. The structure contains a bound ATP, suggesting MJ0577 is an ATPase or an ATP-mediated molecular switch, which we confirm by biochemical experiments. Furthermore, the structure reveals different ATP binding motifs that are shared among many homologous hypothetical proteins in this family. This result indicates that structure-based assignment of molecular function is a viable approach for the large-scale biochemical assignment of proteins and for discovering new motifs, a basic premise of structural genomics.

A rapid classification protocol for the CATH Domain Database to support structural genomics

In order to support the structural genomic initiatives, both by rapidly classifying newly determined structures and by suggesting suitable targets for structure determination, we have recently developed several new protocols for classifying structures in the CATH domain database (http://www.biochem.ucl.ac.uk/bsm/cath). These aim to increase the speed of classification of new structures using fast algorithms for structure comparison (GRATH) and to improve the sensitivity in recognising distant structural relatives by incorporating sequence information from relatives in the genomes (DomainFinder). In order to ensure the integrity of the database given the expected increase in data, the CATH Protein Family Database (CATH-PFDB), which currently includes 25 320 structural domains and a further 160 000 sequence relatives has now been installed in a relational ORACLE database. This was essential for developing more rigorous validation procedures and for allowing efficient querying of the database, particularly for genome analysis. The associated Dictionary of Homologous Superfamilies [Bray,J.E., Todd,A.E., Pearl,F.M.G., Thornton,J.M. and Orengo,C.A. (2000) Protein Eng., 13, 153–165], which provides multiple structural alignments and functional information to assist in assigning new relatives, has also been expanded recently and now includes information for 903 homo­logous superfamilies. In order to improve coverage of known structures, preliminary classification levels are now provided for new structures at interim stages in the classification protocol. Since a large proportion of new structures can be rapidly classified using profile-based sequence analysis [e.g. PSI-BLAST: Altschul,S.F., Madden,T.L., Schaffer,A.A., Zhang,J., Zhang,Z., Miller,W. and Lipman,D.J. (1997) Nucleic Acids Res., 25, 3389–3402], this provides preliminary classification for easily recognisable homologues, which in the latest release of CATH (version 1.7) represented nearly three-quarters of the non-identical structures.

PlantsP: a functional genomics database for plant phosphorylation

The PlantsP database is a curated database that combines information derived from sequences with experimental functional genomics information. PlantsP focuses on plant protein kinases and protein phosphatases. The database will specifically provide a resource for information on a collection of T-DNA insertion mutants (knockouts) in each protein kinase and phosphatase in Arabidopsis thaliana. PlantsP also provides a curated view of each protein that includes a comprehensive annotation of functionally related sequence motifs, sequence family definitions, alignments and phylogenetic trees, and descriptive information drawn directly from the literature. PlantsP is available at http://PlantsP.sdsc.edu.