Comparative genomics of metabolic networks of free-living and parasitic eukaryotes

Background Obligate endoparasites often lack particular metabolic pathways as compared to free-living organisms. This phenomenon comprises anabolic as well as catabolic reactions. Presumably, the corresponding enzymes were lost in adaptation to parasitism. Here we compare the predicted core metabolic graphs of obligate endoparasites and non-parasites (free living organisms and facultative parasites) in order to analyze how the parasites' metabolic networks shrunk in the course of evolution. Results Core metabolic graphs comprising biochemical reactions present in the presumed ancestor of parasites and non-parasites were reconstructed from the Kyoto Encyclopedia of Genes and Genomes. While the parasites' networks had fewer nodes (metabolites) and edges (reactions), other parameters such as average connectivity, network diameter and number of isolated edges were similar in parasites and non-parasites. The parasites' networks contained a higher percentage of ATP-consuming reactions and a lower percentage of NAD-requiring reactions. Control networks, shrunk to the size of the parasites' by random deletion of edges, were scale-free but exhibited smaller diameters and more isolated edges. Conclusions The parasites' networks were smaller than those of the non-parasites regarding number of nodes or edges, but not regarding network diameters. Network integrity but not scale-freeness has acted as a selective principle during the evolutionary reduction of parasite metabolism. ATP-requiring reactions in particular have been retained in the parasites' core metabolism while NADH- or NADPH-requiring reactions were lost preferentially.

Mitochondrial preprotein translocase of trypanosomatids has a bacterial origin

Mitochondria are found in all eukaryotic cells and derive from a bacterial endosymbiont [1, 2]. The evolution of a protein import system was a prerequisite for the conversion of the endosymbiont into a true organelle. Tom40, the essential component of the protein translocase of the outer membrane, is conserved in mitochondria of almost all eukaryotes but lacks bacterial orthologs [3-6]. It serves as the gateway through which all mitochondrial proteins are imported. The parasitic protozoa Trypanosoma brucei and its relatives do not have a Tom40-like protein, which raises the question of how proteins are imported by their mitochondria [7, 8]. Using a combination of bioinformatics and in vivo and in vitro studies, we have discovered that T. brucei likely employs a different import channel, termed ATOM (archaic translocase of the outer mitochondria! membrane). ATOM mediates the import of nuclear-encoded proteins into mitochondria and is essential for viability of trypanosomes. It is not related to Tom40 but is instead an ortholog of a subgroup of the 0mp85 protein superfamily that is involved in membrane translocation and insertion of bacterial outer membrane proteins [9]. This suggests that the protein import channel in trypanosomes is a relic of an archaic protein transport system that was operational in the ancestor of all eukaryotes.

Analysis of gene order data supports vertical inheritance of the leukotoxin operon and genome rearrangements in the 5' flanking region in genus Mannheimia

BACKGROUND: The Mannheimia subclades belong to the same bacterial genus, but have taken divergent paths toward their distinct lifestyles. For example, M. haemolytica + M. glucosida are potential pathogens of the respiratory tract in the mammalian suborder Ruminantia, whereas M. ruminalis, the supposed sister group, lives as a commensal in the ovine rumen. We have tested the hypothesis that vertical inheritance of the leukotoxin (lktCABD) operon has occurred from the last common ancestor of genus Mannheimia to any ancestor of the diverging subclades by exploring gene order data. RESULTS: We examined the gene order in the 5' flanking region of the leukotoxin operon and found that the 5' flanking gene strings, hslVU-lapB-artJ-lktC and xylAB-lktC, are peculiar to M. haemolytica + M. glucosida and M. granulomatis, respectively, whereas the gene string hslVU-lapB-lktC is present in M. ruminalis, the supposed sister group of M. haemolytica + M. glucosida, and in the most ancient subclade M. varigena. In M. granulomatis, we found remnants of the gene string hslVU-lapB-lktC in the xylB-lktC intergenic region. CONCLUSION: These observations indicate that the gene string hslVU-lapB-lktC is more ancient than the hslVU-lapB-artJ-lktC and xylAB-lktC gene strings. The presence of (remnants of) the ancient gene string hslVU-lapB-lktC among any subclades within genus Mannheimia supports that it has been vertically inherited from the last common ancestor of genus Mannheimia to any ancestor of the diverging subclades, thus reaffirming the hypothesis of vertical inheritance of the leukotoxin operon. The presence of individual 5' flanking regions in M. haemolytica + M. glucosida and M. granulomatis reflects later genome rearrangements within each subclade. The evolution of the novel 5' flanking region in M. haemolytica + M. glucosida resulted in transcriptional coupling between the divergently arranged artJ and lkt promoters. We propose that the chimeric promoter have led to high level expression of the leukotoxin operon which could explain the increased potential of certain M. haemolytica + M. glucosida strains to cause a particular type of infection.

Characterization of the type I secretion system of the RTX toxin ApxII in "Actinobacillus porcitonsillarum"

Strains of Actinobacillus porcitonsillarum are regularly isolated from the tonsils of healthy pigs. A. porcitonsillarum is non pathogenic but phenotypically it strongly resembles the pathogenic species Actinobacillus pleuropneumoniae, thereby interfering with the diagnosis of the latter. A. porcitonsillarum is hemolytic but unlike A. pleuropneumoniae, it contains only apxII genes and not apxI or apxIII genes. In contrast to the truncated apxII operon of A. pleuropneumoniae, which lacks the type I secretion genes BD, characterization of the apxII operon in A. porcitonsillarum revealed that it contains an intact and complete apxII operon. This shows a typical RTX operon structure with the gene arrangement apxIICABD. The region upstream of the apxII operon is also different from that in A. pleuropneumoniae and contains an additional gene, aspC, encoding a putative aspartate aminotransferase. Trans-complementation experiments in Escherichia coli and A. pleuropneumoniae indicated that the entire apxII operon of A. porcitonsillarum is sufficient to express and secrete the ApxIIA toxin and that the ApxIIA toxin of A. pleuropneumoniae can be secreted by the type I secretion system encoded by apxIIBD. These findings suggest that the complete apxII operon found in A. porcitonsillarum might be an ancestor of the truncated homologue found in A. pleuropneumoniae. The genetic context of the apxII locus in A. porcitonsillarum and A. pleuropneumoniae suggests that in the latter, the contemporary truncated operon is the result of a recombination event within the species, rather than a horizontal transfer of an incomplete operon.

A non-coding genomic duplication at the HMX1 locus is associated with crop ears in highland cattle.

Highland cattle with congenital crop ears have notches of variable size on the tips of both ears. In some cases, cartilage deformation can be seen and occasionally the external ears are shortened. We collected 40 cases and 80 controls across Switzerland. Pedigree data analysis confirmed a monogenic autosomal dominant mode of inheritance with variable expressivity. All affected animals could be traced back to a single common ancestor. A genome-wide association study was performed and the causative mutation was mapped to a 4 Mb interval on bovine chromosome 6. The H6 family homeobox 1 (HMX1) gene was selected as a positional and functional candidate gene. By whole genome re-sequencing of an affected Highland cattle, we detected 6 non-synonymous coding sequence variants and two variants in an ultra-conserved element at the HMX1 locus with respect to the reference genome. Of these 8 variants, only a non-coding 76 bp genomic duplication (g.106720058_106720133dup) located in the conserved region was perfectly associated with crop ears. The identified copy number variation probably results in HMX1 misregulation and possible gain-of-function. Our findings confirm the role of HMX1 during the development of the external ear. As it is sometimes difficult to phenotypically diagnose Highland cattle with slight ear notches, genetic testing can now be used to improve selection against this undesired trait.

Nuclear markers reveal unexpected genetic variation and a Congolese-Nilotic origin of the Lake Victoria cichlid species flock

Phylogenetic analyses based on mitochondrial (mt) DNA have indicated that the cichlid species flock of the Lake Victoria region is derived from a single ancestral species found in East African rivers, closely related to the ancestor of the Lake Malawi cichlid species flock. The Lake Victoria flock contains ten times less mtDNA variation than the Lake Malawi radiation, consistent with current estimates of the ages of the lakes. We present results of a phylogenetic investigation using nuclear (amplified fragment length polymorphism) markers and a wider coverage of riverine haplochromines. We demonstrate that the Lake Victoria–Edward flock is derived from the morphologically and ecologically diverse cichlid genus Thoracochromis from the Congo and Nile, rather than from the phenotypically conservative East African Astatotilapia. This implies that the ability to express much of the morphological diversity found in the species flock may by far pre–date the origin of the flock. Our data indicate that the nuclear diversity of the Lake Victoria–Edward species flock is similar to that of the Lake Malawi flock, indicating that the genetic diversity is considerably older than the 15 000 years that have passed since the lake began to refill. Most of this variation is manifested in trans–species polymorphisms, indicating very recent cladogenesis from a genetically very diverse founder stock. Our data do not confirm strict monophyly of either of the species flocks, but raise the possibility that these flocks have arisen from hybrid swarms.

Epidermolysis bullosa in Danish Hereford calves is caused by a deletion in LAMC2 gene.

BACKGROUND Heritable forms of epidermolysis bullosa (EB) constitute a heterogeneous group of skin disorders of genetic aetiology that are characterised by skin and mucous membrane blistering and ulceration in response to even minor trauma. Here we report the occurrence of EB in three Danish Hereford cattle from one herd. RESULTS Two of the animals were necropsied and showed oral mucosal blistering, skin ulcerations and partly loss of horn on the claws. Lesions were histologically characterized by subepidermal blisters and ulcers. Analysis of the family tree indicated that inbreeding and the transmission of a single recessive mutation from a common ancestor could be causative. We performed whole genome sequencing of one affected calf and searched all coding DNA variants. Thereby, we detected a homozygous 2.4 kb deletion encompassing the first exon of the LAMC2 gene, encoding for laminin gamma 2 protein. This loss of function mutation completely removes the start codon of this gene and is therefore predicted to be completely disruptive. The deletion co-segregates with the EB phenotype in the family and absent in normal cattle of various breeds. Verifying the homozygous private variants present in candidate genes allowed us to quickly identify the causative mutation and contribute to the final diagnosis of junctional EB in Hereford cattle. CONCLUSIONS Our investigation confirms the known role of laminin gamma 2 in EB aetiology and shows the importance of whole genome sequencing in the analysis of rare diseases in livestock.

Maize Domestication and Anti-Herbivore Defences: Leaf-Specific Dynamics during Early Ontogeny of Maize and Its Wild Ancestors

As a consequence of artificial selection for specific traits, crop plants underwent considerable genotypic and phenotypic changes during the process of domestication. These changes may have led to reduced resistance in the cultivated plant due to shifts in resource allocation from defensive traits to increased growth rates and yield. Modern maize (Zea mays ssp. mays) was domesticated from its ancestor Balsas teosinte (Z. mays ssp. parviglumis) approximately 9000 years ago. Although maize displays a high genetic overlap with its direct ancestor and other annual teosintes, several studies show that maize and its ancestors differ in their resistance phenotypes with teosintes being less susceptible to herbivore damage. However, the underlying mechanisms are poorly understood. Here we addressed the question to what extent maize domestication has affected two crucial chemical and one physical defence traits and whether differences in their expression may explain the differences in herbivore resistance levels. The ontogenetic trajectories of 1,4-benzoxazin-3-ones, maysin and leaf toughness were monitored for different leaf types across several maize cultivars and teosinte accessions during early vegetative growth stages. We found significant quantitative and qualitative differences in 1,4-benzoxazin-3-one accumulation in an initial pairwise comparison, but we did not find consistent differences between wild and cultivated genotypes during a more thorough examination employing several cultivars/accessions. Yet, 1,4-benzoxazin-3-one levels tended to decline more rapidly with plant age in the modern maize cultivars. Foliar maysin levels and leaf toughness increased with plant age in a leaf-specific manner, but were also unaffected by domestication. Based on our findings we suggest that defence traits other than the ones that were investigated are responsible for the observed differences in herbivore resistance between teosinte and maize. Furthermore, our results indicate that single pairwise comparisons may lead to false conclusions regarding the effects of domestication on defensive and possibly other traits.

Bovine Staphylococcus aureus: Subtyping, evolution, and zoonotic transfer.

Staphylococcus aureus is globally one of the most important pathogens causing contagious mastitis in cattle. Previous studies using ribosomal spacer (RS)-PCR, however, demonstrated in Swiss cows that Staph. aureus isolated from bovine intramammary infections are genetically heterogeneous, with Staph. aureus genotype B (GTB) and GTC being the most prominent genotypes. Furthermore, Staph. aureus GTB was found to be contagious, whereas Staph. aureus GTC and all the remaining genotypes were involved in individual cow disease. In addition to RS-PCR, other methods for subtyping Staph. aureus are known, including spa typing and multilocus sequence typing (MLST). They are based on sequencing the spa and various housekeeping genes, respectively. The aim of the present study was to compare the 3 analytic methods using 456 strains of Staph. aureus isolated from milk of bovine intramammary infections and bulk tanks obtained from 12 European countries. Furthermore, the phylogeny of animal Staph. aureus was inferred and the zoonotic transfer of Staph. aureus between cattle and humans was studied. The analyzed strains could be grouped into 6 genotypic clusters, with CLB, CLC, and CLR being the most prominent ones. Comparing the 3 subtyping methods, RS-PCR showed the highest resolution, followed by spa typing and MLST. We found associations among the methods but in many cases they were unsatisfactory except for CLB and CLC. Cluster CLB was positive for clonal complex (CC)8 in 99% of the cases and typically positive for t2953; it is the cattle-adapted form of CC8. Cluster CLC was always positive for t529 and typically positive for CC705. For CLR and the remaining subtypes, links among the 3 methods were generally poor. Bovine Staph. aureus is highly clonal and a few clones predominate. Animal Staph. aureus always evolve from human strains, such that every human strain may be the ancestor of a novel animal-adapted strain. The zoonotic transfer of IMI- and milk-associated strains of Staph. aureus between cattle and humans seems to be very limited and different hosts are not considered as a source for mutual, spontaneous infections. Spillover events, however, may happen.

Lethal chondrodysplasia in a family of Holstein cattle is associated with a de novo splice site variant of COL2A1

Background Lethal chondrodysplasia (bulldog syndrome) is a well-known congenital syndrome in cattle and occurs sporadically in many breeds. In 2015, it was noticed that about 12 % of the offspring of the phenotypically normal Danish Holstein sire VH Cadiz Captivo showed chondrodysplasia resembling previously reported bulldog calves. Pedigree analysis of affected calves did not display obvious inbreeding to a common ancestor, suggesting the causative allele was not a rare recessive. The normal phenotype of the sire suggested a dominant inheritance with incomplete penetrance or a mosaic mutation. Results Three malformed calves were examined by necropsy, histopathology, radiology, and computed tomography scanning. These calves were morphologically similar and displayed severe disproportionate dwarfism and reduced body weight. The syndrome was characterized by shortening and compression of the body due to reduced length of the spine and the long bones of the limbs. The vicerocranium had severe dysplasia and palatoschisis. The bones had small irregular diaphyses and enlarged epiphyses consisting only of chondroid tissue. The sire and a total of four affected half-sib offspring and their dams were genotyped with the BovineHD SNP array to map the defect in the genome. Significant genetic linkage was obtained for several regions of the bovine genome including chromosome 5 where whole genome sequencing of an affected calf revealed a COL2A1 point mutation (g.32473300 G > A). This private sequence variant was predicted to affect splicing as it altered the conserved splice donor sequence GT at the 5’-end of COL2A1 intron 36, which was changed to AT. All five available cases carried the mutant allele in heterozygous state and all five dams were homozygous wild type. The sire VH Cadiz Captivo was shown to be a gonadal and somatic mosaic as assessed by the presence of the mutant allele at levels of about 5 % in peripheral blood and 15 % in semen. Conclusions The phenotypic and genetic findings are comparable to a previously reported COL2A1 missense mutation underlying lethal chondrodysplasia in the offspring of a mosaic French Holstein sire (Igale Masc). The identified independent spontaneous splice site variant in COL2A1 most likely caused chondrodysplasia and must have occurred during the early foetal development of the sire. This study provides a first example of a dominant COL2A1 splice site variant as candidate causal mutation of a severe lethal chondrodysplasia phenotype. Germline mosaicism is a relatively frequent mechanism in the origin of genetic disorders and explains the prevalence of a certain fraction of affected offspring. Paternal dominant de novo mutations are a risk in cattle breeding, especially because the ratio of defective offspring may be very high and be associated with significant animal welfare problems.