Mitochondrial Genomics: structure, inheritance and phylogenetic utility of the Mitochondrial genome in Bivalvia and Insecta

This PhD Thesis is the result of my research activity in the last three years. My main research interest was centered on the evolution of mitochondrial genome (mtDNA), and on its usefulness as a phylogeographic and phylogenetic marker at different taxonomic levels in different taxa of Metazoa. From a methodological standpoint, my main effort was dedicated to the sequencing of complete mitochondrial genomes, and the approach to whole-genome sequencing was based on the application of Long-PCR and shotgun sequences. Moreover, this research project is a part of a bigger sequencing project of mtDNAs in many different Metazoans’ taxa, and I mostly dedicated myself to sequence and analyze mtDNAs in selected taxa of bivalves and hexapods (Insecta). Sequences of bivalve mtDNAs are particularly limited, and my study contributed to extend the sampling. Moreover, I used the bivalve Musculista senhousia as model taxon to investigate the molecular mechanisms and the evolutionary significance of their aberrant mode of mitochondrial inheritance (Doubly Uniparental Inheritance, see below). In Insects, I focused my attention on the Genus Bacillus (Insecta Phasmida). A detailed phylogenetic analysis was performed in order to assess phylogenetic relationships within the genus, and to investigate the placement of Phasmida in the phylogenetic tree of Insecta. The main goal of this part of my study was to add to the taxonomic coverage of sequenced mtDNAs in basal insects, which were only partially analyzed.

Fine-scale population structure analysis of seven local Swiss sheep breeds using genome-wide SNP data

As part of the global sheep Hapmap project, 24 individuals from each of seven indigenous Swiss sheep breeds (Bundner Oberländer sheep (BOS), Engadine Red sheep (ERS), Swiss Black-Brown Mountain sheep (SBS), Swiss Mirror sheep (SMS), Swiss White Alpine (SWA) sheep, Valais Blacknose sheep (VBS) and Valais Red sheep (VRS)), were genotyped using Illumina’s Ovine SNP50 BeadChip. In total, 167 animals were subjected to a detailed analysis for genetic diversity using 45 193 informative single nucleotide polymorphisms. The results of the phylogenetic analyses supported the known proximity between populations such as VBS and VRS or SMS and SWA. Average genomic relatedness within a breed was found to be 12 percent (BOS), 5 percent (ERS), 9 percent (SBS), 10 percent (SMS), 9 percent (SWA), 12 percent (VBS) and 20 percent (VRS). Furthermore, genomic relationships between breeds were found for single individuals from SWA and SMS, VRS and VBS as well as VRS and BOS. In addition, seven out of 40 indicated parent–offspring pairs could not be confirmed. These results were further supported by results from the genome-wide population cluster analysis. This study provides a better understanding of fine-scale population structures within and between Swiss sheep breeds. This relevant information will help to increase the conservation activities of the local Swiss sheep breeds.

Large-scale protein structure modeling of the Saccharomyces cerevisiae genome

The function of a protein generally is determined by its three-dimensional (3D) structure. Thus, it would be useful to know the 3D structure of the thousands of protein sequences that are emerging from the many genome projects. To this end, fold assignment, comparative protein structure modeling, and model evaluation were automated completely. As an illustration, the method was applied to the proteins in the Saccharomyces cerevisiae (baker’s yeast) genome. It resulted in all-atom 3D models for substantial segments of 1,071 (17%) of the yeast proteins, only 40 of which have had their 3D structure determined experimentally. Of the 1,071 modeled yeast proteins, 236 were related clearly to a protein of known structure for the first time; 41 of these previously have not been characterized at all.

Effects of Genome Position and the DNA Damage Checkpoint on the Structure and Frequency of sod2 Gene Amplification in Fission Yeast

The Schizosaccharomyces pombe sod2 gene, located near the telomere on the long arm of chromosome I, encodes a Na+ (or Li+)/H+ antiporter. Amplification of sod2 has previously been shown to confer resistance to LiCl. We analyzed 20 independent LiCl-resistant strains and found that the only observed mechanism of resistance is amplification of sod2. The amplicons are linear, extrachromosomal elements either 225 or 180 kb long, containing both sod2 and telomere sequences. To determine whether proximity to a telomere is necessary for sod2 amplification, a strain was constructed in which the gene was moved to the middle of the same chromosomal arm. Selection of LiCl-resistant strains in this genetic background also yielded amplifications of sod2, but in this case the amplified DNA was exclusively chromosomal. Thus, proximity to a telomere is not a prerequisite for gene amplification in S. pombe but does affect the mechanism. Relative to wild-type cells, mutants with defects in the DNA damage aspect of the rad checkpoint control pathway had an increased frequency of sod2 amplification, whereas mutants defective in the S-phase completion checkpoint did not. Two models for generating the amplified DNA are presented.

The four hexamerin genes in the honey bee: structure, molecular evolution and function deduced from expression patterns in queens, workers and drones

Background: Hexamerins are hemocyanin-derived proteins that have lost the ability to bind copper ions and transport oxygen; instead, they became storage proteins. The current study aimed to broaden our knowledge on the hexamerin genes found in the honey bee genome by exploring their structural characteristics, expression profiles, evolution, and functions in the life cycle of workers, drones and queens. Results: The hexamerin genes of the honey bee (hex 70a, hex 70b, hex 70c and hex 110) diverge considerably in structure, so that the overall amino acid identity shared among their deduced protein subunits varies from 30 to 42%. Bioinformatics search for motifs in the respective upstream control regions (UCRs) revealed six overrepresented motifs including a potential binding site for Ultraspiracle (Usp), a target of juvenile hormone (JH). The expression of these genes was induced by topical application of JH on worker larvae. The four genes are highly transcribed by the larval fat body, although with significant differences in transcript levels, but only hex 110 and hex 70a are re-induced in the adult fat body in a caste-and sex-specific fashion, workers showing the highest expression. Transcripts for hex 110, hex 70a and hex70b were detected in developing ovaries and testes, and hex 110 was highly transcribed in the ovaries of egg-laying queens. A phylogenetic analysis revealed that HEX 110 is located at the most basal position among the holometabola hexamerins, and like HEX 70a and HEX 70c, it shares potential orthology relationship with hexamerins from other hymenopteran species. Conclusions: Striking differences were found in the structure and developmental expression of the four hexamerin genes in the honey bee. The presence of a potential binding site for Usp in the respective 5' UCRs, and the results of experiments on JH level manipulation in vivo support the hypothesis of regulation by JH. Transcript levels and patterns in the fat body and gonads suggest that, in addition to their primary role in supplying amino acids for metamorphosis, hexamerins serve as storage proteins for gonad development, egg production, and to support foraging activity. A phylogenetic analysis including the four deduced hexamerins and related proteins revealed a complex pattern of evolution, with independent radiation in insect orders.

Changes in integument structure during the imaginal molt of the honey bee

The changing pattern of developing cuticle and associated epidermis is described during the imaginal molt in the honey bee. Observations began immediately after the pupal molt, and included histological analyses of the integument during apolysis and the subsequent deposition and differentiation of the adult cuticle. Apolysis coincides with a marked increase in the thickness and reorganization of the epidermal layer, reflecting changes in cell structure. The epidermis remains thickened during the period of cuticle deposition, suggesting intense biosynthetic activity, but turns into a very thin layer during cuticle differentiation, clearly indicating that secretory activity for cuticle formation is terminating. The thoracic cuticle differentiates earlier and becomes thicker than the abdominal. The observed changes in integument structure provide insights that permit an improved physiological characterization for staging pupal and pharate adult development.

An empirical evaluation of imputation accuracy for association statistics reveals increased type-I error rates in genome-wide associations

Background: Genome wide association studies (GWAS) are becoming the approach of choice to identify genetic determinants of complex phenotypes and common diseases. The astonishing amount of generated data and the use of distinct genotyping platforms with variable genomic coverage are still analytical challenges. Imputation algorithms combine directly genotyped markers information with haplotypic structure for the population of interest for the inference of a badly genotyped or missing marker and are considered a near zero cost approach to allow the comparison and combination of data generated in different studies. Several reports stated that imputed markers have an overall acceptable accuracy but no published report has performed a pair wise comparison of imputed and empiric association statistics of a complete set of GWAS markers. Results: In this report we identified a total of 73 imputed markers that yielded a nominally statistically significant association at P < 10(-5) for type 2 Diabetes Mellitus and compared them with results obtained based on empirical allelic frequencies. Interestingly, despite their overall high correlation, association statistics based on imputed frequencies were discordant in 35 of the 73 (47%) associated markers, considerably inflating the type I error rate of imputed markers. We comprehensively tested several quality thresholds, the haplotypic structure underlying imputed markers and the use of flanking markers as predictors of inaccurate association statistics derived from imputed markers. Conclusions: Our results suggest that association statistics from imputed markers showing specific MAF (Minor Allele Frequencies) range, located in weak linkage disequilibrium blocks or strongly deviating from local patterns of association are prone to have inflated false positive association signals. The present study highlights the potential of imputation procedures and proposes simple procedures for selecting the best imputed markers for follow-up genotyping studies.

Near full-length genome analysis of low prevalent human immunodeficiency virus type 1 subclade F1 in Sao Paulo, Brazil

Background: The genetic diversity of the human immunodeficiency virus type 1 (HIV-1) is critical to lay the groundwork for the design of successful drugs or vaccine. In this study we aimed to characterize and define the molecular prevalence of HIV-1 subclade F1 currently circulating in Sao Paulo, Brazil. Methods: A total of 36 samples were selected from 888 adult patients residing in Sao Paulo who had previously been diagnosed in two independent studies in our laboratory as being infected with subclade F1 based on pol subgenomic fragment sequencing. Proviral DNA was amplified from the purified genomic DNA of all 36 blood samples by 5 fragments overlapping PCR followed by direct sequencing. Sequence data were obtained from the 5 fragments of pure subclade F1 and phylogenetic trees were constructed and compared with previously published sequences. Subclades F1 that exhibited mosaic structure with other subtypes were omitted from any further analysis Results: Our methods of fragment amplification and sequencing confirmed that only 5 sequences inferred from pol region as subclade F1 also holds true for the genome as a whole and, thus, estimated the true prevalence at 0.56%. The results also showed a single phylogenetic cluster of the Brazilian subclade F1 along with non-Brazilian South American isolates in both subgenomic and the full-length genomes analysis with an overall intrasubtype nucleotide divergence of 6.9%. The nucleotide differences within the South American and Central African F1 strains, in the C2-C3 env, were 8.5% and 12.3%, respectively. Conclusion: All together, our findings showed a surprisingly low prevalence rate of subclade F1 in Brazil and suggest that these isolates originated in Central Africa and subsequently introduced to South America.

Genome Size, Karyotype Polymorphism and Chromosomal Evolution in Trypanosoma cruzi

Background: The Trypanosoma cruzi genome was sequenced from a hybrid strain (CL Brener). However, high allelic variation and the repetitive nature of the genome have prevented the complete linear sequence of chromosomes being determined. Determining the full complement of chromosomes and establishing syntenic groups will be important in defining the structure of T. cruzi chromosomes. A large amount of information is now available for T. cruzi and Trypanosoma brucei, providing the opportunity to compare and describe the overall patterns of chromosomal evolution in these parasites. Methodology/Principal Findings: The genome sizes, repetitive DNA contents, and the numbers and sizes of chromosomes of nine strains of T. cruzi from four lineages (TcI, TcII, TcV and TcVI) were determined. The genome of the TcI group was statistically smaller than other lineages, with the exception of the TcI isolate Tc1161 (Jose-IMT). Satellite DNA content was correlated with genome size for all isolates, but this was not accompanied by simultaneous amplification of retrotransposons. Regardless of chromosomal polymorphism, large syntenic groups are conserved among T. cruzi lineages. Duplicated chromosome-sized regions were identified and could be retained as paralogous loci, increasing the dosage of several genes. By comparing T. cruzi and T. brucei chromosomes, homologous chromosomal regions in T. brucei were identified. Chromosomes Tb9 and Tb11 of T. brucei share regions of syntenic homology with three and six T. cruzi chromosomal bands, respectively. Conclusions: Despite genome size variation and karyotype polymorphism, T. cruzi lineages exhibit conservation of chromosome structure. Several syntenic groups are conserved among all isolates analyzed in this study. The syntenic regions are larger than expected if rearrangements occur randomly, suggesting that they are conserved owing to positive selection. Mapping of the syntenic regions on T. cruzi chromosomal bands provides evidence for the occurrence of fusion and split events involving T. brucei and T. cruzi chromosomes.

Genome-wide SNP genotyping highlights the role of natural selection in Plasmodium falciparum population divergence

Background: The malaria parasite Plasmodium falciparum exhibits abundant genetic diversity, and this diversity is key to its success as a pathogen. Previous efforts to study genetic diversity in P. falciparum have begun to elucidate the demographic history of the species, as well as patterns of population structure and patterns of linkage disequilibrium within its genome. Such studies will be greatly enhanced by new genomic tools and recent large-scale efforts to map genomic variation. To that end, we have developed a high throughput single nucleotide polymorphism (SNP) genotyping platform for P. falciparum. Results: Using an Affymetrix 3,000 SNP assay array, we found roughly half the assays (1,638) yielded high quality, 100% accurate genotyping calls for both major and minor SNP alleles. Genotype data from 76 global isolates confirm significant genetic differentiation among continental populations and varying levels of SNP diversity and linkage disequilibrium according to geographic location and local epidemiological factors. We further discovered that nonsynonymous and silent (synonymous or noncoding) SNPs differ with respect to within-population diversity, interpopulation differentiation, and the degree to which allele frequencies are correlated between populations. Conclusions: The distinct population profile of nonsynonymous variants indicates that natural selection has a significant influence on genomic diversity in P. falciparum, and that many of these changes may reflect functional variants deserving of follow-up study. Our analysis demonstrates the potential for new high-throughput genotyping technologies to enhance studies of population structure, natural selection, and ultimately enable genome-wide association studies in P. falciparum to find genes underlying key phenotypic traits.

The mitochondrial genome of the phytopathogenic basidiomycete Moniliophthora perniciosa is 109 kb in size and contains a stable integrated plasmid

We present here the sequence of the mitochondrial genome of the basidiomycete phytopathogenic hemibiotrophic fungus Moniliophthora perniciosa, causal agent of the Witches` Broom Disease in Theobroma cacao. The DNA is a circular molecule of 109103 base pairs, with 31.9 % GC, and is the largest sequenced so far. This size is due essentially to the presence of numerous non-conserved hypothetical ORFs. It contains the 14 genes coding for proteins involved in the oxidative phosphorylation, the two rRNA genes, one ORF coding for a ribosomal protein (rps3), and a set of 26 tRNA genes that recognize codons for all amino acids. Seven homing endonucleases are located inside introns. Except atp8, all conserved known genes are in the same orientation. Phylogenetic analysis based on the cox genes agrees with the commonly accepted fungal taxonomy. An uncommon feature of this mitochondrial genome is the presence of a region that contains a set of four, relatively small, nested, inverted repeats enclosing two genes coding for polymerases with an invertron-type structure and three conserved hypothetical genes interpreted as the stable integration of a mitochondrial linear plasmid. The integration of this plasmid seems to be a recent evolutionary event that could have implications in fungal biology. This sequence is available under GenBank accession number AY376688. (c) 2008 The British Mycological Society. Published by Elsevier Ltd. All rights reserved.

Gulliver, a long terminal repeat retrotransposon from the genome of the oriental blood fluke Schistosoma japonicum

We characterized the consensus sequence and structure of a long terminal repeat (LTR) retrotransposon from the genome of the human blood fluke, Schistosoma japonicum, and have earned this element, Gulliver. The full length, consensus Gulliver LTR retrotransposon was 4788 bp, and it was flanked at its 5'- and 3'-ends by LTRs of 259 bp. Each LTR included RNA polymerase II promoter sequences, a CAAT signal and a TATA box, Gulliver exhibited features characteristic of a functional LTR retrotransposon including two read through (termination) ORFs encoding retroviral gag and pol proteins of 312 and 1071 amino acid residues, respectively. The gag ORF encoded motifs conserved in nucleic acid binding proteins, while the pol ORF encoded conserved domains of aspartic protease, reverse transcriptase (RT), RNaseH and integrase, in that order, a pol pattern conserved in the gypsy lineage of LTR retrotransposons. Whereas the sequence and structure of Gulliver was similar to that of gypsy, phylogenetic analysis revealed that Gulliver did not group particularly closely with the gypsy family. Rather, its closest relatives were a LTR retrotransposon from Caenorhabditis elegans, mag from Bombyx mori and, to a lesser extent, easel from the salmon Oncorhynchus keta. Dot blot hybridizations indicated that Gulliver was present at between 100 and several thousand copies in the S. japonicum genome, and Southern hybridization analysis suggested its probable presence in the genome of Schistosoma mansoni. Transcripts encoding the RT domain of Gulliver were detected by RT-PCR in larval and adult stages of S. japonicum, indicating that (at least) the RT domain of Gulliver is transcribed. This is the first report of the sequence and structure of an LTR retrotransposon from any schistosome or indeed from any species belonging to the phylum Platyhelminthes. (C) 2001 Elsevier Science B.V. All rights reserved.

Complete DNA sequence and gene organization of the mitochondrial genome of the liverfluke, Fasciola hepatica L. (Platyhelminthes; Trematoda)

The complete nucleotide sequence of the mitochondrial (mt) DNA molecule of the liverfluke, Fasciola hepatica (phylum Platyhelminthes, class Trematoda, family Fasciolidae), was determined, It comprises 14462 bp, contains 12 protein-encoding, 2 ribosomal and 22 transfer RNA genes, and is the second complete flatworm (and the first trematode) mitochondrial sequence to be described in detail. All of the genes are transcribed from the same strand. Of the genes typically found in mitochondrial genomes of eumetazoans, only atp8 is absent. The nad4L and nad4 genes overlap by 40 nt. Most intergenic sequences are very short. Two larger non-coding regions are present. The longer one (817 nt) is located between trnG and cox3 and consists of 8 identical tandem repeats of 85 nt, rich in G and C, followed by 1 imperfect repeat. The shorter non-coding region (187 nt) exhibits no special features and is separated from the longer region by trnG. The gene arrangement resembles that of some other trematodes including the eastern Asian Schistosoma species (and cyclophyllidean cestode species) but it is strikingly different from that of the African schistosomes, represented by Schistosoma mansoni. The genetic code is as inferred previously for flatworms. Transfer RNA genes range in length from 58 to 70 nt, their products producing characteristic 'clover leaf' structures, except for tRNA(S-VON) and tRNA(S-AGN) lacking the DHU arm.

pido, a non-long terminal repeat retrotransposon of the chicken repeat 1 family from the genome of the Oriental blood fluke, Schistosoma japonicum

A newly described non-long terminal repeat (non-LTR) retrotransposon element was isolated from the genome of the Oriental schistosome, Schistosoma japonicum. At least 1000 partial copies of the element, which was named pido, were dispersed throughout the genome of S. japonicum. As is usual with non-LTR retrotransposons, it is expected that many pido elements will be 5'-truncated. A consensus sequence of 3564 bp of the truncated pido element was assembled from several genomic fragments that contained pido-hybridizing sequences. The sequence encoded part of the first open reading frame (ORF), the entire second ORF and, at its 3'-terminus, a tandemly repetitive, A-rich (TA(6)TA(5)TA(8)) tail, The ORF1 of pido encoded a nucleic acid binding protein and ORF2 encoded a retroviral-like polyprotein that included apurinic/apyrimidinic endonuclease (EN) and reverse transcriptase (RT) domains, in that order. Based on its sequence and structure, and phylogenetic analyses of both the RT and EN domains, pido belongs to the chicken repeat 1 (CR1)-like lineage of elements known from the chicken, turtle, puffer fish, mosquitoes and other taxa. pido shared equal similarity with CRI from chicken, an uncharacterized retrotransposon from Caenorhabditis elegans and SR1 (a non-LTR retrotransposon) from the related blood fluke Schistosoma mansoni; the level of similarity between pido and SR1 indicated that these two schistosome retrotransposons were related but not orthologous. The findings indicate that schistosomes have been colonized by at least two discrete CRI-like elements. Whereas pido did not appear to have a tight target site specificity, at least one copy of pido has inserted into the 3'-untranslated region of a protein-encoding gene (GeriBank AW736757) of as yet unknown identity. mRNA encoding the RT of pido was detected by reverse transcription-polymerase chain reaction in the egg, miracidium. and adult developmental stages of S. japonicum, indicating that the RT domain was transcribed and suggesting that pido was replicating actively and mobile within the S. japonicum genome. (C) 2002 Elsevier Science B.V. All rights reserved.

Significance in replication of the terminal nucleotides of the Flavivirus genome

Point mutations that resulted in a substitution of the conserved 3'-penultimate cytidine in genomic RNA or the RNA negative strand of the self-amplifying replicon of the Flavivirus Kunjin virus completely blocked in vivo replication. Similarly, substitutions of the conserved 3'-terminal uridine in the RNA negative or positive strand completely blocked replication or caused much-reduced replication, respectively. The same preference for cytidine in the 3'-terminal dinucleotide was noted in reports of the in vitro activity of the RNA-dependent RNA polymerase (RdRp) for the other genera of Flaviviridae that also employ a double-stranded RNA (dsRNA) template to initiate asymmetric semiconservative RNA positive-strand synthesis. The Kunjin virus replicon results were interpreted in the context of a proposed model for initiation of RNA synthesis based on the solved crystal structure of the RdRp of phi6 bacteriophage, which also replicates efficiently using a dsRNA template with conserved 3'-penultimate cytidines and a 3'-terminal pyrimidine. A previously untested substitution of the conserved pentanucleotide at the top of the 3'-terminal stem-loop of all Flavivirus species also blocked detectable in vivo replication of the Kunjin virus replicon RNA.