Evolutionary History of Arabidopsis thaliana Aminoacyl-tRNA Synthetase Dual-Targeted Proteins

Aminoacyl-transfer RNA (tRNA) synthetases (aaRS) are key players in translation and act early in protein synthesis by mediating the attachment of amino acids to their cognate tRNA molecules. In plants, protein synthesis may occur in three subcellular compartments (cytosol, mitochondria, and chloroplasts), which requires multiple versions of the protein to be correctly delivered to its proper destination. The organellar aaRS are nuclear encoded and equipped with targeting information at the N-terminal sequence, which enables them to be specifically translocated to their final location. Most of the aaRS families present organellar proteins that are dual targeted to mitochondria and chloroplasts. Here, we examine the dual targeting behavior of aaRS from an evolutionary perspective. Our results show that Arabidopsis thaliana aaRS sequences are a result of a horizontal gene transfer event from bacteria. However, there is no evident bias indicating one single ancestor (Cyanobacteria or Proteobacteria). The dual-targeted aaRS phylogenetic relationship was characterized into two different categories (paralogs and homologs) depending on the state recovered for both dual-targeted and cytosolic proteins. Taken together, our results suggest that the dual-targeted condition is a gain-of-function derived from gene duplication. Selection may have maintained the original function in at least one of the copies as the additional copies diverged.

Strategies for genetic model specification in the screening of genome-wide meta-analysis signals for further replication

Background Meta-analysis is increasingly being employed as a screening procedure in large-scale association studies to select promising variants for follow-up studies. However, standard methods for meta-analysis require the assumption of an underlying genetic model, which is typically unknown a priori. This drawback can introduce model misspecifications, causing power to be suboptimal, or the evaluation of multiple genetic models, which augments the number of false-positive associations, ultimately leading to waste of resources with fruitless replication studies. We used simulated meta-analyses of large genetic association studies to investigate naive strategies of genetic model specification to optimize screenings of genome-wide meta-analysis signals for further replication. Methods Different methods, meta-analytical models and strategies were compared in terms of power and type-I error. Simulations were carried out for a binary trait in a wide range of true genetic models, genome-wide thresholds, minor allele frequencies (MAFs), odds ratios and between-study heterogeneity (tau(2)). Results Among the investigated strategies, a simple Bonferroni-corrected approach that fits both multiplicative and recessive models was found to be optimal in most examined scenarios, reducing the likelihood of false discoveries and enhancing power in scenarios with small MAFs either in the presence or in absence of heterogeneity. Nonetheless, this strategy is sensitive to tau(2) whenever the susceptibility allele is common (MAF epsilon 30%), resulting in an increased number of false-positive associations compared with an analysis that considers only the multiplicative model. Conclusion Invoking a simple Bonferroni adjustment and testing for both multiplicative and recessive models is fast and an optimal strategy in large meta-analysis-based screenings. However, care must be taken when examined variants are common, where specification of a multiplicative model alone may be preferable.

Phylogenetic Analysis of Hepatitis B Virus Genotype F Complete Genome Sequences From Chilean Patients With Chronic Infection

Molecular epidemiological data concerning the hepatitis B virus (HBV) in Chile are not known completely. Since the HBV genotype F is the most prevalent in the country, the goal of this study was to obtain full HBV genome sequences from patients infected chronically in order to determine their subgenotypes and the occurrence of resistance-associated mutations. Twenty-one serum samples from antiviral drug-naive patients with chronic hepatitis B were subjected to full-length PCR amplification, and both strands of the whole genomes were fully sequenced. Phylogenetic analyses were performed along with reference sequences available from GenBank (n = 290). The sequences were aligned using Clustal X and edited in the SE-AL software. Bayesian phylogenetic analyses were conducted by Markov Chain Monte Carlo simulations (MCMC) for 10 million generations in order to obtain the substitution tree using BEAST. The sequences were also analyzed for the presence of primary drug resistance mutations using CodonCode Aligner Software. The phylogenetic analyses indicated that all sequences were found to be the HBV subgenotype F1b, clustered into four different groups, suggesting that diverse lineages of this subgenotype may be circulating within this population of Chilean patients. J. Med. Virol. 83: 1530-1536, 2011. (C) 2011 Wiley-Liss, Inc.

Whole-Genome Characterization of Human and Simian Immunodeficiency Virus Intrahost Diversity by Ultradeep Pyrosequencing

Rapid evolution and high intrahost sequence diversity are hallmarks of human and simian immunodeficiency virus (HIV/SIV) infection. Minor viral variants have important implications for drug resistance, receptor tropism, and immune evasion. Here, we used ultradeep pyrosequencing to sequence complete HIV/SIV genomes, detecting variants present at a frequency as low as 1%. This approach provides a more complete characterization of the viral population than is possible with conventional methods, revealing low-level drug resistance and detecting previously hidden changes in the viral population. While this work applies pyrosequencing to immunodeficiency viruses, this approach could be applied to virtually any viral pathogen.

High frequency of submicroscopic chromosomal imbalances in patients with syndromic craniosynostosis detected by a combined approach of microsatellite segregation analysis, multiplex ligation-dependent probe amplification and array-based comparative genome hybridisation

We present the first comprehensive study, to our knowledge, on genomic chromosomal analysis in syndromic craniosynostosis. In total, 45 patients with craniosynostotic disorders were screened with a variety of methods including conventional karyotype, microsatellite segregation analysis, subtelomeric multiplex ligation-dependent probe amplification) and whole-genome array-based comparative genome hybridisation. Causative abnormalities were present in 42.2% (19/45) of the samples, and 27.8% (10/36) of the patients with normal conventional karyotype carried submicroscopic imbalances. Our results include a wide variety of imbalances and point to novel chromosomal regions associated with craniosynostosis. The high incidence of pure duplications or trisomies suggests that these are important mechanisms in craniosynostosis, particularly in cases involving the metopic suture.

Distinct patterns of somatic alterations in a lymphoblastoid and a tumor genome derived from the same individual

Although patterns of somatic alterations have been reported for tumor genomes, little is known on how they compare with alterations present in non-tumor genomes. A comparison of the two would be crucial to better characterize the genetic alterations driving tumorigenesis. We sequenced the genomes of a lymphoblastoid (HCC1954BL) and a breast tumor (HCC1954) cell line derived from the same patient and compared the somatic alterations present in both. The lymphoblastoid genome presents a comparable number and similar spectrum of nucleotide substitutions to that found in the tumor genome. However, a significant difference in the ratio of non-synonymous to synonymous substitutions was observed between both genomes (P = 0.031). Protein-protein interaction analysis revealed that mutations in the tumor genome preferentially affect hub-genes (P = 0.0017) and are co-selected to present synergistic functions (P < 0.0001). KEGG analysis showed that in the tumor genome most mutated genes were organized into signaling pathways related to tumorigenesis. No such organization or synergy was observed in the lymphoblastoid genome. Our results indicate that endogenous mutagens and replication errors can generate the overall number of mutations required to drive tumorigenesis and that it is the combination rather than the frequency of mutations that is crucial to complete tumorigenic transformation.

The genome of the blood fluke Schistosoma mansoni

Schistosoma mansoni is responsible for the neglected tropical disease schistosomiasis that affects 210 million people in 76 countries. Here we present analysis of the 363 megabase nuclear genome of the blood fluke. It encodes at least 11,809 genes, with an unusual intron size distribution, and new families of micro-exon genes that undergo frequent alternative splicing. As the first sequenced flatworm, and a representative of the Lophotrochozoa, it offers insights into early events in the evolution of the animals, including the development of a body pattern with bilateral symmetry, and the development of tissues into organs. Our analysis has been informed by the need to find new drug targets. The deficits in lipid metabolism that make schistosomes dependent on the host are revealed, and the identification of membrane receptors, ion channels and more than 300 proteases provide new insights into the biology of the life cycle and new targets. Bioinformatics approaches have identified metabolic chokepoints, and a chemogenomic screen has pinpointed schistosome proteins for which existing drugs may be active. The information generated provides an invaluable resource for the research community to develop much needed new control tools for the treatment and eradication of this important and neglected disease.

The Genome Sequence of Taurine Cattle: A Window to Ruminant Biology and Evolution

To understand the biology and evolution of ruminants, the cattle genome was sequenced to about sevenfold coverage. The cattle genome contains a minimum of 22,000 genes, with a core set of 14,345 orthologs shared among seven mammalian species of which 1217 are absent or undetected in noneutherian (marsupial or monotreme) genomes. Cattle-specific evolutionary breakpoint regions in chromosomes have a higher density of segmental duplications, enrichment of repetitive elements, and species-specific variations in genes associated with lactation and immune responsiveness. Genes involved in metabolism are generally highly conserved, although five metabolic genes are deleted or extensively diverged from their human orthologs. The cattle genome sequence thus provides a resource for understanding mammalian evolution and accelerating livestock genetic improvement for milk and meat production.

Complete genome analysis of a rabies virus isolate from Brazilian wild fox

The complete genome sequence of wild-type rabies virus (RABV) isolated from a wild Brazilian hoary fox (Dusicyon sp.), the BR-Pfx1 isolate, was determined and compared with fixed RABV strains. The genome structure and organization of the BR-Pfx1 isolate were composed of 11,924 nt and included the five standard genes of rhabdoviruses. Sequences of mRNA start and stop signals for transcription were highly conserved among all structural protein genes of the BR-Pfx1 isolate. All amino acid residues in the glycoprotein (G) gene associated with pathogenicity were retained in the BR-Pfx1 isolate, while unique amino acid substitutions were found in antigenic region I of the nucleoprotein gene and III of G. These results suggest that although the standard genome structure and organization of the RABV isolate are common between the BR-Pfx1 isolate and fixed RABV strains, the unique amino acid substitutions in functional sites of the BR-Pfx1 isolate may result in different biological characteristics from fixed RABV strains.

Coupling Virus-Induced Gene Silencing to Exogenous Green Fluorescence Protein Expression Provides a Highly Efficient System for Functional Genomics in Arabidopsis and across All Stages of Tomato Fruit Development

Since the advent of the postgenomic era, efforts have focused on the development of rapid strategies for annotating plant genes of unknown function. Given its simplicity and rapidity, virus-induced gene silencing (VIGS) has become one of the preeminent approaches for functional analyses. However, several problems remain intrinsic to the use of such a strategy in the study of both metabolic and developmental processes. The most prominent of these is the commonly observed phenomenon of ""sectoring"" the tissue regions that are not effectively targeted by VIGS. To better discriminate these sectors, an effective marker system displaying minimal secondary effects is a prerequisite. Utilizing a VIGS system based on the tobacco rattle virus vector, we here studied the effect of silencing the endogenous phytoene desaturase gene (pds) and the expression and subsequent silencing of the exogenous green fluorescence protein (gfp) on the metabolism of Arabidopsis (Arabidopsis thaliana) leaves and tomato (Solanum lycopersicum) fruits. In leaves, we observed dramatic effects on primary carbon and pigment metabolism associated with the photobleached phenotype following the silencing of the endogenous pds gene. However, relatively few pleiotropic effects on carbon metabolism were observed in tomato fruits when pds expression was inhibited. VIGS coupled to gfp constitutive expression revealed no significant metabolic alterations after triggering of silencing in Arabidopsis leaves and a mild effect in mature green tomato fruits. By contrast, a wider impact on metabolism was observed in ripe fruits. Silencing experiments with an endogenous target gene of interest clearly demonstrated the feasibility of cosilencing in this system; however, carefully constructed control experiments are a prerequisite to prevent erroneous interpretation.

Genomic Analysis of Wild Tomato Introgressions Determining Metabolism- and Yield-Associated Traits

With the aim of determining the genetic basis of metabolic regulation in tomato fruit, we constructed a detailed physical map of genomic regions spanning previously described metabolic quantitative trait loci of a Solanum pennellii introgression line population. Two genomic libraries from S. pennellii were screened with 104 colocated markers from five selected genomic regions, and a total of 614 bacterial artificial chromosome (BAC)/cosmids were identified as seed clones. Integration of sequence data with the genetic and physical maps of Solanum lycopersicum facilitated the anchoring of 374 of these BAC/cosmid clones. The analysis of this information resulted in a genome-wide map of a nondomesticated plant species and covers 10% of the physical distance of the selected regions corresponding to approximately 1% of the wild tomato genome. Comparative analyses revealed that S. pennellii and domesticated tomato genomes can be considered as largely colinear. A total of 1,238,705 bp from both BAC/cosmid ends and nine large insert clones were sequenced, annotated, and functionally categorized. The sequence data allowed the evaluation of the level of polymorphism between the wild and cultivated tomato species. An exhaustive microsynteny analysis allowed us to estimate the divergence date of S. pennellii and S. lycopersicum at 2.7 million years ago. The combined results serve as a reference for comparative studies both at the macrosyntenic and microsyntenic levels. They also provide a valuable tool for fine-mapping of quantitative trait loci in tomato. Furthermore, they will contribute to a deeper understanding of the regulatory factors underpinning metabolism and hence defining crop chemical composition.

Molecular characterization of a retrotransposon in the Rhynchosciara americana genome and its association with telomere

Non-LTR retrotransposons, also known as long interspersed nuclear elements (LINEs), are transposable elements that encode a reverse transcriptase and insert into genomic locations via RNA intermediates. The sequence analysis of a cDNA library constructed from mRNA of the salivary glands of R. americana showed the presence of putative class I elements. The cDNA clone with homology to a reverse transcriptase was the starting point for the present study. Genomic phage was isolated and sequenced and the molecular structure of the element was characterized as being a non-LTR retrotransposable element. Southern blot analysis indicated that this transposable element is represented by repeat sequences in the genome of R. americana. Chromosome tips were consistently positive when this element was used as probe in in-situ hybridization. Real-time RT-PCR showed that this retrotransposon is transcribed at different periods of larval development. Most interesting, the silencing of this retrotransposon in R. americana by RNA interference resulted in reduced transcript levels and in accelerated larval development.

Mariner-like elements in Rhynchosciara americana (Sciaridae) genome: molecular and cytological aspects

Two mariner-like elements, Ramar1 and Ramar2, are described in the genome of Rhynchosciara americana, whose nucleotide consensus sequences were derived from multiple defective copies containing deletions, frame shifts and stop codons. Ramar1 contains several conserved amino acid blocks which were identified, including a specific D,D(34)D signature motif. Ramar2 is a defective mariner-like element, which contains a deletion overlapping in most of the internal region of the transposase ORF while its extremities remain intact. Predicted transposase sequences demonstrated that Ramar1 and Ramar2 phylogenetically present high identity to mariner-like elements of mauritiana subfamily. Southern blot analysis indicated that Ramar1 is widely represented in the genome of Rhynchosciara americana. In situ hybridizations showed Ramar1 localized in several chromosome regions, mainly in pericentromeric heterochromatin and their boundaries, while Ramar2 appeared as a single band in chromosome A.

Analysis of the genome of Spodoptera frugiperda nucleopolyhedrovirus (SfMNPV-19) and of the high genomic heterogeneity in group II nucleopolyhedroviruses

The genome of the most virulent among 22 Brazilian geographical isolates of Spodoptera frugiperda nucleopolyhedrovirus, isolate 19 (SfMNPV-1 9), was completely sequenced and shown to comprise 132 565 bp and 141 open reading frames (ORFs). A total of 11 ORFs with no homology to genes in the GenBank database were found. Of those, four had typical baculovirus; promoter motifs and polyadenylation sites. Computer-simulated restriction enzyme cleavage patterns of SfMNPV-1 9 were compared with published physical maps of other SfMNPV isolates. Differences were observed in terms of the restriction profiles and genome size. Comparison of SfMNPV-1 9 with the sequence of the SfMNPV isolate 3AP2 indicated that they differed due to a 1427 bp deletion, as well as by a series of smaller deletions and point mutations. The majority of genes of SfMNPV-1 9 were conserved in the closely related Spodoptera exigua NPV (SeMNPV) and Agrotis segetum NPV (AgseMNPV-A), but a few regions experienced major changes and rearrangements. Synthenic maps for the genomes of group 11 NPVs revealed that gene collinearity was observed only within certain clusters. Analysis of the dynamics of gene gain and loss along the phylogenetic tree of the NPVs showed that group 11 had only five defining genes and supported the hypothesis that these viruses form ten highly divergent ancient lineages. Crucially, more than 60% of the gene gain events followed a power-law relation to genetic distance among baculoviruses, indicative of temporal organization in the gene accretion process.

Sugarcane Phosphoribosyl Pyrophosphate Synthetase: Molecular Characterization of a Phosphate-independent PRS

Phosphoribosyl pyrophosphate synthetase (PRS-EC:2.7.6.1) is an important enzyme present in several metabolic pathways, thus forming a complex family of isoenzymes. However, plant PRS enzymes have not been extensively investigated. In this study, a sugarcane prs gene has been characterized from the Sugar Cane Expressed Sequence Tag Genome Project. This gene contains a 984-bp open reading frame encoding a 328-amino acid protein. The predicted amino acid sequence has 77% and 78% amino acid sequence identity to Arabidopsis thaliana and Spinacia oleracea PRS4, respectively. The assignment of sugarcane PRS as a phosphate-independent PRS isoenzyme (Class II PRS) is verified following enzyme assay and phylogenetic reconstruction of PRS homologues. To gain further insight into the structural framework of the phosphate independence of sugarcane PRS, a molecular model is described. This model reveals the formation of two conserved domains elucidating the structural features involved in sugarcane PRS phosphate independence. The recombinant PRS retains secondary structure elements and a quaternary arrangement consistent with known PRS homologues, based on circular dichroism measurements.