Comparisons of host mitochondrial, nuclear and endosymbiont bacterial genes reveal cryptic fig wasp species and the effects of Wolbachia on host mtDNA evolution and diversity

Background Figs and fig-pollinating wasp species usually display a highly specific one-to-one association. However, more and more studies have revealed that the "one-to-one" rule has been broken. Co-pollinators have been reported, but we do not yet know how they evolve. They may evolve from insect speciation induced or facilitated by Wolbachia which can manipulate host reproduction and induce reproductive isolation. In addition, Wolbachia can affect host mitochondrial DNA evolution, because of the linkage between Wolbachia and associated mitochondrial haplotypes, and thus confound host phylogeny based on mtDNA. Previous research has shown that fig wasps have the highest incidence of Wolbachia infection in all insect taxa, and Wolbachia may have great influence on fig wasp biology. Therefore, we look forward to understanding the influence of Wolbachia on mitochondrial DNA evolution and speciation in fig wasps. Results We surveyed 76 pollinator wasp specimens from nine Ficus microcarpa trees each growing at a different location in Hainan and Fujian Provinces, China. We found that all wasps were morphologically identified as Eupristina verticillata, but diverged into three clades with 4.22-5.28% mtDNA divergence and 2.29-20.72% nuclear gene divergence. We also found very strong concordance between E. verticillata clades and Wolbachia infection status, and the predicted effects of Wolbachia on both mtDNA diversity and evolution by decreasing mitochondrial haplotypes. Conclusions Our study reveals that the pollinating wasp E. verticillata on F. microcarpa has diverged into three cryptic species, and Wolbachia may have a role in this divergence. The results also indicate that Wolbachia strains infecting E. verticillata have likely resulted in selective sweeps on host mitochondrial DNA.

Endoparasite communities of five fish species (Labridae : Cheilininae) from Lizard Island: how important is the ecology and phylogeny of the hosts?

The parasite community of animals is generally influenced by host physiology, ecology, and phylogeny. Therefore, sympatric and phylogenetically related hosts with similar ecologies should have similar parasite communities. To test this hypothesis we surveyed the endoparasites of 5 closely related cheilinine fishes (Labridae) from the Great Barrier Reef. They were Cheilinus chlorounts, C. trilobatus, C. fasciatils, Epibulus insidiator and OxYcheilinus diagrainnia. VVe examined the relationship between parasitological variables (richness, abundance and diversity) and host characteristics (bodv weight, diet and phuylogeny). The 5 fishes had 31 parasite species with 9-18 parasite species per fish species. Cestode larvae (mostly Tetraphyllidea) were the most abundant and prevalent parasites followed by nematodes and digeneans. Parasites, body size and diet of hosts differed between fish species. In general, body weight, diet and host phylogeny each explained some of the variation in richness and composition of parasites among the fishes. The 2 most closely related species, Cheilinus chlorourus and C. trilobatus, had broadly similar parasites but the Other fish species differed significantly in all variables. However, there was no all -encompassing pattern. This may, be because different lineages of parasites may react differently to ecological variables. We also argue that adult parasites may respond principally to host diet. In contrast, larval parasite composition may respond both to host diet and predator-prey interactions because this is the path by which many, parasites complete their life-cycles. Finally, variation in parasite phylogeny and parasite life-cycles among hosts likely increase the complexity of the system making it difficult to find all-encompassing patterns between host characteristics and parasites, particularly when all the species in rich parasite communities are considered.

Functional genomics of commensal Lactobacilli

Catabolic flexibility affords a bacterium the ability to utilise different sugar sources as carbon for energy. This is important for commensal lactobacilli like Lactobacillus ruminis which can be exposed to a variety of carbohydrates in vivo. However, little is known about the fermentation capabilities, metabolic pathways, genetic diversity or potential survival mechanisms used by L. ruminis in vivo. A combination of in vitro and in silico techniques was used to identify the catabolic pathways of L. ruminis. I also compared 16 L. ruminis strains using a panel of biochemical and survival assays, genetically, whole genome sequencing and RNA sequencing. Multi locus sequence typing revealed that strains clustered according to their host sources. Transcriptome analysis by RNAseq of two motile strains under three growth conditions, including swarming, identified the up-regulation of carbohydrate-related genes under swarming conditions. This suggests that carbohydrate flexibility may have an uncharacterised role in L. ruminis swarming. Following on from the assessment of L. ruminis catabolic flexibility, the porcine diet was supplemented with galactooligosaccharides or L. ruminis ATCC 25644 plus galactooligosaccharides. Supplementation of the porcine diet with galactooligosaccharide had no effect on microbiota diversity. In contrast, the L. ruminis plus galactooligosaccharide treatment significantly reduced the microbiota diversity. Diet is a major factor that affects the diversity of the gut microbiota. In order to get a more thorough understanding of diet and gut health in animals such as racehorses and domesticated herbivores, I determined the core microbiota of animals consuming different feeds. Interestingly, the gut microbiota diversity correlated with the host phylogeny of the animal. The genome of Lactobacillus equi (2.19 Mb), isolated from a healthy Irish thoroughbred was also sequenced and annotated, and comprised 2,263 predicted genes. The large repertoire of predicted carbohydrate-related genes may offer L. equi an advantage in the complex and harsh hindgut environment. In summary, this thesis uses functional genomics to assess the effect that carbohydrates have on commensal lactobacilli and the microbiota as a whole.

Similarities among ectoparasite fauna of sigmodontine rodents: phylogenetic and geographical influences

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Transposable elements: structure and dynamic of the autonomous retrotransposon R2 in Arthropoda with non-canonical reproduction

Eukaryotic ribosomal DNA constitutes a multi gene family organized in a cluster called nucleolar organizer region (NOR); this region is composed usually by hundreds to thousands of tandemly repeated units. Ribosomal genes, being repeated sequences, evolve following the typical pattern of concerted evolution. The autonomous retroelement R2 inserts in the ribosomal gene 28S, leading to defective 28S rDNA genes. R2 element, being a retrotransposon, performs its activity in the genome multiplying its copy number through a “copy and paste” mechanism called target primed reverse transcription. It consists in the retrotranscription of the element’s mRNA into DNA, then the DNA is integrated in the target site. Since the retrotranscription can be interrupted, but the integration will be carried out anyway, truncated copies of the element will also be present in the genome. The study of these truncated variants is a tool to examine the activity of the element. R2 phylogeny appears, in general, not consistent with that of its hosts, except some cases (e.g. Drosophila spp. and Reticulitermes spp.); moreover R2 is absent in some species (Fugu rubripes, human, mouse, etc.), while other species have more R2 lineages in their genome (the turtle Mauremys reevesii, the Japanese beetle Popilia japonica, etc). R2 elements here presented are isolated in 4 species of notostracan branchiopods and in two species of stick insects, whose reproductive strategies range from strict gonochorism to unisexuality. From sequencing data emerges that in Triops cancriformis (Spanish gonochoric population), in Lepidurus arcticus (two putatively unisexual populations from Iceland) and in Bacillus rossius (gonochoric population from Capalbio) the R2 elements are complete and encode functional proteins, reflecting the general features of this family of transposable elements. On the other hand, R2 from Italian and Austrian populations of T. cancriformis (respectively unisexual and hermaphroditic), Lepidurus lubbocki (two elements within the same Italian population, gonochoric but with unfunctional males) and Bacillus grandii grandii (gonochoric population from Ponte Manghisi) have sequences that encode incomplete or non-functional proteins in which it is possible to recognize only part of the characteristic domains. In Lepidurus couesii (Italian gonochoric populations) different elements were found as in L. lubbocki, and the sequencing is still in progress. Two hypothesis are given to explain the inconsistency of R2/host phylogeny: vertical inheritance of the element followed by extinction/diversification or horizontal transmission. My data support previous study that state the vertical transmission as the most likely explanation; nevertheless horizontal transfer events can’t be excluded. I also studied the element’s activity in Spanish populations of T. cancriformis, in L. lubbocki, in L. arcticus and in gonochoric and parthenogenetic populations of B. rossius. In gonochoric populations of T. cancriformis and B. rossius I found that each individual has its own private set of truncated variants. The situation is the opposite for the remaining hermaphroditic/parthenogenetic species and populations, all individuals sharing – in the so far analyzed samples - the majority of variants. This situation is very interesting, because it isn’t concordant with the Muller’s ratchet theory that hypothesizes the parthenogenetic populations being either devoided of transposable elements or TEs overloaded. My data suggest a possible epigenetic mechanism that can block the retrotransposon activity, and in this way deleterious mutations don’t accumulate.

Recurrent invasion and extinction of a selfish gene

Homing endonuclease genes show super-Mendelian inheritance, which allows them to spread in populations even when they are of no benefit to the host organism. To test the idea that regular horizontal transmission is necessary for the long-term persistence of these genes, we surveyed 20 species of yeasts for the ω-homing endonuclease gene and associated group I intron. The status of ω could be categorized into three states (functional, nonfunctional, or absent), and status was not clustered on the host phylogeny. Moreover, the phylogeny of ω differed significantly from that of the host, strong evidence of horizontal transmission. Further analyses indicate that horizontal transmission is more common than transposition, and that it occurs preferentially between closely related species. Parsimony analysis and coalescent theory suggest that there have been 15 horizontal transmission events in the ancestry of our yeast species, through simulations indicate that this value is probably an underestimate. Overall, the data support a cyclical model of invasion, degeneration, and loss, followed by reinvasion, and each of these transitions is estimated to occur about once every 2 million years. The data are thus consistent with the idea that frequent horizontal transmission is necessary for the long-term persistence of homing endonuclease genes, and further, that this requirement limits these genes to organisms with easily accessible germ lines. The data also show that mitochondrial DNA sequences are transferred intact between yeast species; if other genes do not show such high levels of horizontal transmission, it would be due to lack of selection, rather than lack of opportunity.

Phylogeny of a taxonomically difficult group and evolution of host location mechanism

The parasitic wasps are one of the largest insect groups and their life histories are remarkably variable. Common to all parasitic wasps is that they kill their hosts, which are usually beetles, butterflies and sometimes spiders. Hosts are often at a larval or pupal stage and live in concealed conditions, such as in plant tissue. Parasitic wasps have two main ways of finding their host. 1) They can detect chemical compounds emitted by damaged plant material or released by larvae living in plant tissue, and 2) detect the larvae by sound vibrations. Even though pupae are immobile and silent, and therefore do not cause vibration, parasitoids have, however, adapted to find passive developmental stages by producing vibration themselves by knocking the substrate with their antennae, and then detecting the echoes with their legs. This echolocation allows a parasitoid to locate its potential hosts that are deeply buried in wood. This study focuses on the relationships of the subfamily Cryptinae (Hymenoptera: Ichneumonidae) and related taxa, and the evolution of host location mechanism. There are no earlier studies of the phylogeny of the Cryptinae, and the position of related taxa are unclear. According to the earlier classification, which is entirely intuitional, the Cryptinae is divided into three tribes: Cryptini, Hemigasterini and Phygadeuontini. Further, these tribes are subdiveded into numerous subtribes. This work, based on molecular characters, shows that the cryptine tribes Cryptini, Phygadeuon¬tini and Hemigasterini come out largely as monophyletic groups, thus agreeing with the earlier classification. The earlier subtribal classification had no support. In addition, it is shown that modified antennal structures are associated with host usage of wood-boring coleopteran hosts. The cryptines have a clear modification series on their antennal tips from a simply tip to a hammer-like structure. The species with strongly modified antennae belong mostly to the tribe Cryptini and they utilise wood-boring beetles as hosts. Also, field observations on insect behaviour support this result.

Host range, symptom expression and RNA 3 sequence analyses of six Australian strains of Cucumber mosaic virus

We have characterised six Australian Cucumber mosaic virus (CMV) strains belonging to different subgroups, determined by the sequence of their complete RNA 3 and by their host range and the symptoms they cause on species in the Solanaceae, Cucurbitaceae and on sweet corn. These data allowed classification of strains into the known three CMV subgroups and identification of plant species able to differentiate the Australian strains by symptoms and host range. Western Australian strains 237 and Twa and Queensland strains 207 and 242 are closely related members of CMV subgroup IA, which cause similar severe symptoms on Nicotiana species. Strains 207 and 237 (subgroup IA) were the only strains tested which systemically infected sweet corn. Strain 243 caused the most severe symptoms of all strains on Nicotiana species, tomato and capsicum and appears to be the first confirmed subgroup IB strain reported in Australia. Based on pair-wise distance analysis and phylogeny of RNA 3, as well as mild disease symptoms on Nicotiana species, CMV 241 was assigned to subgroup II, as the previously described Q-CMV and LY-CMV.

Novel Papillomaviruses in Free-Ranging Iberian Bats: No Virus-Host Co-evolution, No Strict Host Specificity, and Hints for Recombination

Papillomaviruses (PVs) are widespread pathogens. However, the extent of PV infections in bats remains largely unknown. This work represents the first comprehensive study of PVs in Iberian bats. We identified four novel PVs in the mucosa of free-ranging Eptesicus serotinus (EserPV1, EserPV2, and EserPV3) and Rhinolophus ferrumequinum (RferPV1) individuals and analyzed their phylogenetic relationships within the viral family. We further assessed their prevalence in different populations of E. serotinus and its close relative E. isabellinus. Although it is frequent to read that PVs co-evolve with their host, that PVs are highly species-specific, and that PVs do not usually recombine, our results suggest otherwise. First, strict virus-host co-evolution is rejected by the existence of five, distantly related bat PV lineages and by the lack of congruence between bats and bat PVs phylogenies. Second, the ability of EserPV2 and EserPV3 to infect two different bat species (E. serotinus and E. isabellinus) argues against strict host specificity. Finally, the description of a second noncoding region in the RferPV1 genome reinforces the view of an increased susceptibility to recombination in the E2-L2 genomic region. These findings prompt the question of whether the prevailing paradigms regarding PVs evolution should be reconsidered.

Novel Papillomaviruses in Free-Ranging Iberian Bats: No Virus-Host Co-evolution, No Strict Host Specificity, and Hints for Recombination

Papillomaviruses (PVs) are widespread pathogens. However, the extent of PV infections in bats remains largely unknown. This work represents the first comprehensive study of PVs in Iberian bats. We identified four novel PVs in the mucosa of free-ranging Eptesicus serotinus (EserPV1, EserPV2, and EserPV3) and Rhinolophus ferrumequinum (RferPV1) individuals and analyzed their phylogenetic relationships within the viral family. We further assessed their prevalence in different populations of E. serotinus and its close relative E. isabellinus. Although it is frequent to read that PVs co-evolve with their host, that PVs are highly species-specific, and that PVs do not usually recombine, our results suggest otherwise. First, strict virus-host co-evolution is rejected by the existence of five, distantly related bat PV lineages and by the lack of congruence between bats and bat PVs phylogenies. Second, the ability of EserPV2 and EserPV3 to infect two different bat species (E. serotinus and E. isabellinus) argues against strict host specificity. Finally, the description of a second noncoding region in the RferPV1 genome reinforces the view of an increased susceptibility to recombination in the E2-L2 genomic region. These findings prompt the question of whether the prevailing paradigms regarding PVs evolution should be reconsidered.

Evidence of host specificity and congruence between phylogenies of bitterling and freshwater mussels

Evidence of host specificity and congruence between phylogenies of bitterling and freshwater mussels. Zoological Studies 45(3): 428-434. Bitterling (Cyprinidae: Acheilognathinae) are freshwater fishes with a unique spawning relationship with freshwater mussels on whose gills they lay their eggs. During the breeding season of bitterling fishes, we collected 843 mussels belonging to 16 species from Lake Qinglan, central China and examined their gill chambers for the presence of bitterling larvae. Three species of bitterling larvae were identified; Acheilognathus tonkinensis, Ach. cf. meridianus, and Ach. barbatulus, in 3 species of mussel: Unio douglasiae, Lamprotula caveata, and L. tortuosa, suggesting host specialization. Using our own and other published data, we compared the respective phylogenies of bitterling and mussels, but failed to show clear congruence. However, broad specializations are evident, with Acheilognathus and Tanakia showing preferences for mussels with a relatively simple gill structure (Ableminae), and Rhodeus spp. showing preferences for mussels of the Anodontinae and Unioninae, which have more-complex gill structures.

Characterisation of bacteriophage-host interactions in Lactococcus lactis

Lactococcus lactis is used extensively world-wide for the production of fermented dairy products. Bacteriophages (phages) infecting L. lactis can result in slow or incomplete fermentations, or may even cause total fermentation failure. Therefore, bacteriophages disrupting L. lactis fermentation are of economic concern. This thesis employed a multifaceted approach to investigate various molecular aspects of phage-host interaction in L. lactis. The genome sequence of an Irish dairy starter strain, the prophage-cured L. lactis subsp. cremoris UC509.9, was studied. The 2,250,427 bp circular chromosome represents the smallest among its sequenced lactococcal equivalents. The genome displays clear genetic adaptation to the dairy niche in the form of extensive reductive evolution. Gene prediction identified 2066 protein-encoding genes, including 104 which showed significant homology to transposase-specifying genes. Over 9 % of the identified genes appear to be inactivated through stop codons or frame shift mutations. Many pseudogenes were found in genes that are assigned to carbohydrate and amino acid transport and metabolism orthologous groups, reflecting L. lactis UC509.9’s adaptation to the lactose and casein-rich dairy environment. Sequence analysis of the eight plasmids of L. lactis revealed extensive adaptation to the dairy environment. Key industrial phenotypes were mapped and novel lactococcal plasmid-associated genes highlighted. In addition to chromosomally-encoded bacteriophage resistance systems, six functional such systems were identified, including two abortive infection systems, AbiB and AbiD1, explaining the observed phage resistance of L. lactis UC509.9 Molecular analysis suggests that the constitutive expression of AbiB is not lethal to cells, suggesting the protein is expressed in an un/inactivated form. Analysis of 936 species phage sk1-escape mutants of AbiB revealed that all such mutants harbour mutations in orf6, which encodes the major capsid protein. Results suggest that the major capsid protein is required for activation of the AbiB system, although this requires furrther investigations. Temporal transcriptomes of L. lactis UC509.9 undergoing lytic infection with either one of two distinct bacteriophages, Tuc2009 and c2, was determined and compared to the transcriptome of uninfected UC509.9 cells. Whole genome microarrays performed at various time-points post-infection demonstrated a rather modest impact on host transcription. Alterations in the UC509.9 transcriptome during lytic infection appear phage-specific, with a relatively small number of differentially transcribed genes shared between infection with either Tuc2009 or c2. Transcriptional profiles of both bacteriophages during lytic infection was shown to generally correlate with previous studies and allowed the confirmation of previously predicted promoter sequences. Bioinformatic analysis of genomic regions encoding the presumed cell wall polysaccharide (CW PS) biosynthesis gene cluster of several strains of L. lactis was performed. Results demonstrate the presence of three dominant genetic types of this gene cluster, termed type A, B and C. These regions were used for the development of a multiplex PCR to identify CW PS genotype of various lactococcal strains. Analysis of 936 species phage receptor binding protein phylogeny (RBP) and CW PS genotype revealed an apparent correlation between RBP phylogeny and CW PS type, thereby providing a partial explanation for the observed narrow host range of 936 phages. Further analysis of the genetic locus encompassing the presumed CW PS biosynthesis operon of eight strains identified as belonging to the CW PS C (geno)type, revealed the presence of a variable region among the examined strains. The obtained comparative analysis allowed for the identification of five subgroups of the C type, named C1 to C5. We purified an acidic polysaccharide from the cell wall of L. lactis 3107 (C2 subtype) and confirmed that it is structurally different from the CW PS of the C1 subtype L. lactis MG1363. Combinations of genes from the variable region of C2 subtype were amplified from L. lactis 3107 and introduced into a mutant of the C1 subtype L. lactis NZ9000 (a direct derivative of MG1363) deficient in CW PS biosynthesis. The resulting recombinant mutant synthesized a CW PS with a composition characteristic for that of the C2 subtype L. lactis 3107 and not the wildtype C1 L. lactis NZ9000. The recombinant mutant exhibited a changed phage resistance/sensitivity profile consistent with that of L. lactis 3107, which unambiguously demonstrated that L. lactis 3107 CW PS is the host cell surface receptor of two bacteriophages belonging to the P335 species as well as phages that are member of the 936 species. The research presented in this thesis has significantly advanced our understanding of L. lactis bacteriophage-host interactions in several ways. Firstly, the examination of plasmidencoded bacteriophage resistance systems has allowed inferences to be made regarding the mode of action of AbiB, thereby providing a platform for further elucidation of the molecular trigger of this system. Secondly, the phage infection transcriptome data presented, in addition to previous work, has made L. lactis a model organism in terms of transcriptomic studies of bacteriophage-host interactions. And finally, the research described in this thesis has for the first time explicitly revealed the nature of a carbohydrate bacteriophage receptor in L. lactis, while also providing a logical explanation for the observed narrow host ranges exhibited by 936 and P335 phages. Future research in discerning the structures of other L. lactis CW PS, combined with the determination of the molecular interplay between receptor binding proteins of these phages and CW PS will allow an in depth understanding of the mechanism by which the most prevalent lactococcal phages identify and adsorb to their specific host.

Low frequency of paleoviral infiltration across the avian phylogeny.

BACKGROUND: Mammalian genomes commonly harbor endogenous viral elements. Due to a lack of comparable genome-scale sequence data, far less is known about endogenous viral elements in avian species, even though their small genomes may enable important insights into the patterns and processes of endogenous viral element evolution. RESULTS: Through a systematic screening of the genomes of 48 species sampled across the avian phylogeny we reveal that birds harbor a limited number of endogenous viral elements compared to mammals, with only five viral families observed: Retroviridae, Hepadnaviridae, Bornaviridae, Circoviridae, and Parvoviridae. All nonretroviral endogenous viral elements are present at low copy numbers and in few species, with only endogenous hepadnaviruses widely distributed, although these have been purged in some cases. We also provide the first evidence for endogenous bornaviruses and circoviruses in avian genomes, although at very low copy numbers. A comparative analysis of vertebrate genomes revealed a simple linear relationship between endogenous viral element abundance and host genome size, such that the occurrence of endogenous viral elements in bird genomes is 6- to 13-fold less frequent than in mammals. CONCLUSIONS: These results reveal that avian genomes harbor relatively small numbers of endogenous viruses, particularly those derived from RNA viruses, and hence are either less susceptible to viral invasions or purge them more effectively.

The phylogeny, structure and function of trematode cysteine proteases, with particular emphasis on the Fasciola hepatica cathepsin L family

Helminth parasites (nematodes, flatworms and cestodes) infect over 1 billion of the world's population causing high morbidity and mortality. The large tissue-dwelling worms express papain-like cysteine peptidases, termed cathepsins that play important roles in virulence including host entry, tissue migration and the suppression of host immune responses. Much of our knowledge of helminth cathepsins comes from studies using flatworms or trematode (fluke) parasites. The developmentally-regulated expression of these proteases correlates with the passage of parasites through host tissues and their encounters with different host macromolecules. Recent phylogenetic, biochemical and structural studies indicate that trematode cathepsins exhibit overlapping but distinct substrate specificities due to divergence within the protease active site. Here we provide an overview of the evolution, biochemistry and structure of these important enzymes and highlight how recent advances in proteomics and gene silencing techniques are allowing researchers to probe their biological functions. We focus mainly on members of the cathepsin L gene family of the animal and human pathogen, Fasciola hepatica, because of our deep understanding of their function, biochemistry and structure.