Mosaic evolution of ruminant stomach lysozyme genes

Adaptive Evolution of C-Type Lysozyme in Vampire Bats

In mammals, chicken-type (c-type) lysozymes are part of the innate immune system, killing bacteria by degrading peptidoglycan in their cell walls. Many of the studies on the evolution of c-type lysozymes have focused on its new digestive function, including the duplicated stomach lysozymes in ruminants. Similarly, in bats, gene duplications and subsequent adaptive evolution of c-type lysozyme have been reported in a clade of insectivorous species, which might have been driven by the need to digest chitin. However, no studies on the evolution of c-type lysozyme have been carried out in the second largest and dietary diverse bat family Phyllostomidae, which includes insectivorous, frugivorous, nectarivorous and sanguivorous species. Here, we sequenced and analyzed c-type lysozyme genes from four phyllostomid bats, the common vampire bat, the white-winged vampire bat, the lesser long-nosed bat and the big fruit-eating bat. Only a single lysozyme gene was identified in each of these species. Evidence for positive selection on mature lysozyme was found on lineages leading to vampire bats, but not other bats with single copy lysozyme genes. Moreover, several amino acid substitutions found in mature lysozymes from the sanguivorous clade are predicted to have functional impacts, adding further evidence for the adaptive evolution of lysozyme in vampire bats. Functional adaptation of vampire bat lysozymes could be associated with anti-microbial defense, possibly driven by the specialized sanguivory-related habits of vampire bats.

Cultivation and sequencing of rumen microbiome members from the Hungate1000 Collection

Rumen microbiome biology gets a boost with the release of 410 high-quality reference genomes from the Hungate1000 project.

Productivity of ruminant livestock depends on the rumen microbiota, which ferment indigestible plant polysaccharides into nutrients used for growth. Understanding the functions carried out by the rumen microbiota is important for reducing greenhouse gas production by ruminants and for developing biofuels from lignocellulose. We present 410 cultured bacteria and archaea, together with their reference genomes, representing every cultivated rumen-associated archaeal and bacterial family. We evaluate polysaccharide degradation, short-chain fatty acid production and methanogenesis pathways, and assign specific taxa to functions. A total of 336 organisms were present in available rumen metagenomic data sets, and 134 were present in human gut microbiome data sets. Comparison with the human microbiome revealed rumen-specific enrichment for genes encoding de novo synthesis of vitamin B₁₂, ongoing evolution by gene loss and potential vertical inheritance of the rumen microbiome based on underrepresentation of markers of environmental stress. We estimate that our Hungate genome resource represents ∼75% of the genus-level bacterial and archaeal taxa present in the rumen.

Genomic organization and evolution of ruminant lysozyme c genes

Ruminant stomach lysozyme is a long established model of adaptive gene evolution. Evolution of stomach lysozyme function required changes in the site of expression of the lysozyme c gene and changes in the enzymatic properties of the enzyme. In ruminant mammals, these changes were associated with a change in the size of the lysozyme c gene family. The recent release of near complete genome sequences from several ruminant species allows a more complete examination of the evolution and diversification of the lysozyme c gene family. Here we characterize the size of the lysozyme c gene family in extant ruminants and demonstrate that their pecoran ruminant ancestor had a family of at least 10 lysozyme c genes, which included at least two pseudogenes. Evolutionary analysis of the ruminant lysozyme c gene sequences demonstrate that each of the four exons of the lysozyme c gene has a unique evolutionary history, indicating that they participated independently in concerted evolution. These analyses also show that episodic changes in the evolutionary constraints on the protein sequences occurred, with lysozyme c genes expressed in the abomasum of the stomach of extant ruminant species showing the greatest levels of selective constraints.

Comparative digestive physiology

In vertebrates and invertebrates, morphological and functional features of gastrointestinal (GI) tracts generally reflect food chemistry, such as content of carbohydrates, proteins, fats, and material(s) refractory to rapid digestion (e.g., cellulose). The expression of digestive enzymes and nutrient transporters approximately matches the dietary load of their respective substrates, with relatively modest excess capacity. Mechanisms explaining differences in hydrolase activity between populations and species include gene copy number variations and single-nucleotide polymorphisms. Transcriptional and posttranscriptional adjustments mediate phenotypic changes in the expression of hydrolases and transporters in response to dietary signals. Many species respond to higher food intake by flexibly increasing digestive compartment size. Fermentative processes by symbiotic microorganisms are important for cellulose degradation but are relatively slow, so animals that rely on those processes typically possess special enlarged compartment(s) to maintain a microbiota and other GI structures that slow digesta flow. The taxon richness of the gut microbiota, usually identified by 16S rRNA gene sequencing, is typically an order of magnitude greater in vertebrates than invertebrates, and the interspecific variation in microbial composition is strongly influenced by diet. Many of the nutrient transporters are orthologous across different animal phyla, though functional details may vary (e.g., glucose and amino acid transport with K+ rather than Na+ as a counter ion). Paracellular absorption is important in many birds. Natural toxins are ubiquitous in foods and may influence key features such as digesta transit, enzymatic breakdown, microbial fermentation, and absorption.

Novel gene expression responses in the ovine abomasal mucosa to infection with the gastric nematode Teladorsagia circumcincta

Infection of sheep with the gastric nematode Teladorsagia circumcincta results in distinct Th2-type changes in the mucosa, including mucous neck cell and mast cell hyperplasia, eosinophilia, recruitment of IgA/IgE producing cells and neutrophils, altered T-cell subsets and mucosal hypertrophy. To address the protective mechanisms generated in animals on previous exposure to this parasite, gene expression profiling was carried out using samples of abomasal mucosa collected pre- and post- challenge from animals of differing immune status, using an experimental model of T. circumcincta infection. Recently developed ovine cDNA arrays were used to compare the abomasal responses of sheep immunised by trickle infection with worm-naïve sheep, following a single oral challenge of 50 000 T. circumcincta L3. Key changes were validated using qRT-PCR techniques. Immune animals demonstrated highly significant increases in levels of transcripts normally associated with cytotoxicity such as granulysin and granzymes A, B and H, as well as mucous-cell derived transcripts, predominantly calcium-activated chloride channel 1 (CLCA1). Challenge infection also induced up-regulation of transcripts potentially involved in initiating or modulating the immune response, such as heat shock proteins, complement factors and the chemokine CCL2. In contrast, there was marked infection-associated down-regulation of gene expression of members of the gastric lysozyme family. The changes in gene expression levels described here may reflect roles in direct anti-parasitic effects, immuno-modulation or tissue repair. (Funding; DEFRA/SHEFC (VT0102) and the BBSRC (BB/E01867X/1)).

Evolution of the mammalian lysozyme gene family

BACKGROUND

Lysozyme c (chicken-type lysozyme) has an important role in host defense, and has been extensively studied as a model in molecular biology, enzymology, protein chemistry, and crystallography. Traditionally, lysozyme c has been considered to be part of a small family that includes genes for two other proteins, lactalbumin, which is found only in mammals, and calcium-binding lysozyme, which is found in only a few species of birds and mammals. More recently, additional testes-expressed members of this family have been identified in human and mouse, suggesting that the mammalian lysozyme gene family is larger than previously known.

RESULTS

Here we characterize the extent and diversity of the lysozyme gene family in the genomes of phylogenetically diverse mammals, and show that this family contains at least eight different genes that likely duplicated prior to the diversification of extant mammals. These duplicated genes have largely been maintained, both in intron-exon structure and in genomic context, throughout mammalian evolution.

CONCLUSIONS

The mammalian lysozyme gene family is much larger than previously appreciated and consists of at least eight distinct genes scattered around the genome. Since the lysozyme c and lactalbumin proteins have acquired very different functions during evolution, it is likely that many of the other members of the lysozyme-like family will also have diverse and unexpected biological properties.

Duplication, concerted evolution and purifying selection drive the evolution of mosquito vitellogenin genes

Background

Mosquito vitellogenin (Vtg) genes belong to a small multiple gene family that encodes the major yolk protein precursors required for egg production. Multiple Vtg genes have been cloned and characterized from several mosquito species, but their origin and molecular evolution are poorly understood.

Results

Here we used in silico and molecular cloning techniques to identify and characterize the evolution of the Vtg gene family from the genera Culex, Aedes/Ochlerotatus, and Anopheles. We identified the probable ancestral Vtg gene among different mosquito species by its conserved association with a novel gene approximately one kilobase upstream of the start codon. Phylogenetic analysis indicated that the Vtg gene family arose by duplication events, but that the pattern of duplication was different in each mosquito genera. Signatures of purifying selection were detected in Culex, Aedes and Anopheles. Gene conversion is a major driver of concerted evolution in Culex, while unequal crossover is likely the major driver of concerted evolution in Anopheles. In Aedes, smaller fragments have undergone gene conversion events.

Conclusions

The study shows concerted evolution and purifying selection shaped the evolution of mosquito Vtg genes following gene duplication. Additionally, similar evolutionary patterns were observed in the Vtg genes from other invertebrate and vertebrate organisms, suggesting that duplication, concerted evolution and purifying selection may be the major evolutionary forces driving Vtg gene evolution across highly divergent taxa.

Identification and expression analysis of the g-type and c-type lysozymes in grass carp Ctenopharyngodon idellus

Lysozyme is an important molecule of innate immune system for the defense against bacterial infections. We identified two genes encoding g-type and c-type lysozymes from grass carp Ctenopharyngodon idellus by the RACE method. The deduced amino acids of both lysozymes possessed typical structural residues and conserved catalytic sites similar to their counterparts across the species. In contrast, there was only 8.6% similarity of amino acid sequence between these two lysozymes. Phylogenetic analyses revealed that these two genes evolved at different rate. C-type lysozyme of grass carp was diverged early in the evolutionary history. Moreover, the expression patterns of these two genes differed. The mRNA levels of both genes were increased after bacterial infection, but the up-regulation of g-type lysozyme was much stronger than that of c-type lysozyme. We also showed that the c-type and g-type recombinant lysozymes possessed different lytic activities against fish bacterial pathogens. These results confirmed that both lysozymes play important roles in the defense of grass carp against bacterial infections. The g-type lysozyme may be induced for the defense against bacterial infections, while c-type lysozyme might be the main molecule for the house-keeping defense under normal conditions. These two types of lysozymes likely use different mechanisms to regulate their expressions.

Elucidating the Evolutionary Relationships among Bos taurus Digestive Organs Using Unigene Expression Data

Although the nature of ruminant evolution is still disputed, current theory based on physiology and genetic analysis suggests that the abomasum is the evolutionarily oldest stomach compartment, the rumen evolved some time after the abomasum, and the omasum is the evolutionarily youngest stomach compartment. In addition, there is some evidence of relaxed selective constraint in the stomach-like organ and the foregut shortly after the foregut formation event. Along with the assumption of a mean, stochastic rate of evolution, analysis of differences in genetic profiles among digestive body organs can give clues to the relationships among these organs. The presence of large numbers of uniquely expressed entries in the abomasum and rumen indicates either a period of relaxed selective constraint or greater evolutionary age. Additionally, differences in expression profiles indicate that the abomasum, rumen, and intestine are more closely related to each other, while the reticulum and omasum are more closely related to the rumen. Functional analysis using Gene Ontology (GO) categories also supports the proposed evolutionary relationships by identifying shared functions, such as muscle activity and development, lipid transport, and urea metabolism, between all sections of the digestive tract investigated.

Identifying concerted evolution and gene conversion in mammalian gene pairs lasting over 100 million years

BACKGROUND

Concerted evolution occurs in multigene families and is characterized by stretches of homogeneity and higher sequence similarity between paralogues than between orthologues. Here we identify human gene pairs that have undergone concerted evolution, caused by ongoing gene conversion, since at least the human-mouse divergence. Our strategy involved the identification of duplicated genes with greater similarity within a species than between species. These genes were required to be present in multiple mammalian genomes, suggesting duplication early in mammalian divergence. To eliminate genes that have been conserved due to strong purifying selection, our analysis also required at least one intron to have retained high sequence similarity between paralogues.

RESULTS

We identified three human gene pairs undergoing concerted evolution (BMP8A/B, DDX19A/B, and TUBG1/2). Phylogenetic investigations reveal that in each case the duplication appears to have occurred prior to eutherian mammalian radiation, with exactly two paralogues present in all examined species. This indicates that all three gene duplication events were established over 100 million years ago.

CONCLUSION

The extended duration of concerted evolution in multiple distant lineages suggests that there has been prolonged homogenization of specific segments within these gene pairs. Although we speculate that selection for homogenization could have been utilized in order to maintain crucial homo- or hetero- binding domains, it remains unclear why gene conversion has persisted for such extended periods of time. Through these analyses, our results demonstrate additional examples of a process that plays a definite, although unspecified, role in molecular evolution.

Stomach lysozymes of the three-toed sloth (Bradypus variegatus), an arboreal folivore from the Neotropics

Lysozymes are antimicrobial defences that act as digestive enzymes when expressed in the stomach of herbivores with pre-gastric fermentation. We studied this enzyme in the complex stomach of the three-toed sloth (Bradypus variegatus), a folivore with pre-gastric fermentation. Lysozymes were identified by SDS-PAGE and immunoblotting in all portions: diverticulum, pouch, glandular and muscular prepyloric area with 14.3 kDa of molecular mass. Purified lysozymes from all areas but the diverticulum were characterized by MALDI-TOF, optimal pH, optimal ionic strength, and specific activity. The differences observed suggested at least three isoforms. The optimal pHs were similar to the pH of the stomach portion where the enzymes were isolated. The lysozyme from the pouch (fermentation chamber) exhibited higher specific activity and concentration than the others. The specific activity of the enzyme from the acid muscular prepyloric portion was comparable to that reported in the cow abomasums; however, its concentration was lower than that observed in cow. This distinctive pattern of secretion/specific activity and overall low concentration suggests different roles for the lysozymes in this herbivore compared to Artiodactyla. We postulate that sloth stomach lysozymes may still be antimicrobial defences by protecting the microbial flora of the fermentation chamber against foreign bacteria.

Comparative genomics and concerted evolution of beta-tubulin paralogs in Leishmania spp

Background

Tubulin isotypes and expression patterns are highly regulated in diverse organisms. The genome sequence of the protozoan parasite Leishmania major contains three distinct β-tubulin loci. To investigate the diversity of β-tubulin genes, we have compared the published genome sequence to draft genome sequences of two further species L. infantum and L. braziliensis. Untranscribed regions and coding sequences for each isoform were compared within and between species in relation to the known diversity of β-tubulin transcripts in Leishmania spp.

Results

All three β-tubulin loci were present in L. infantum and L. braziliensis, showing conserved synteny with the L. major sequence, hence confirming that these loci are paralogous. Flanking regions suggested that the chromosome 21 locus is an amastigote-specific isoform and more closely related (either structurally or functionally) to the chromosome 33 'array' locus than the chromosome 8 locus. A phylogenetic network of all isoforms indicated that paralogs from L. braziliensis and L. mexicana were monophyletic, rather than clustering by locus.

Conclusion

L. braziliensis and L. mexicana sequences appeared more similar to each other than each did to its closest relative in another species; this indicates that these sequences have evolved convergently in each species, perhaps through ectopic gene conversion; a process not yet evident among the more recently derived L. major and L. infantum isoforms. The distinctive non-coding regions of each β-tubulin locus showed that it is the regulatory regions of these loci that have evolved most during the diversification of these genes in Leishmania, while the coding regions have been conserved and concerted. The various loci in Leishmania satisfy a need for innovative expression of β-tubulin, rather than elaboration of its structural role.

Duplication and divergence of 2 distinct pancreatic ribonuclease genes in leaf-eating African and Asian colobine monkeys

Unique among primates, the colobine monkeys have adapted to a predominantly leaf-eating diet by evolving a foregut that utilizes bacterial fermentation to breakdown and absorb nutrients from such a food source. It has been hypothesized that pancreatic ribonuclease (pRNase) has been recruited to perform a role as a digestive enzyme in foregut fermenters, such as artiodactyl ruminants and the colobines. We present molecular analyses of 23 pRNase gene sequences generated from 8 primate taxa, including 2 African and 2 Asian colobine species. The pRNase gene is single copy in all noncolobine primate species assayed but has duplicated more than once in both the African and Asian colobine monkeys. Phylogenetic reconstructions show that the pRNase-coding and noncoding regions are under different evolutionary constraints, with high levels of concerted evolution among gene duplicates occurring predominantly in the noncoding regions. Our data suggest that 2 functionally distinct pRNases have been selected for in the colobine monkeys, with one group adapting to the role of a digestive enzyme by evolving at an increased rate with loss of positive charge, namely arginine residues. Conclusions relating our data to general hypotheses of evolution following gene duplication are discussed.

A tandem gene duplication followed by recruitment of a retrotransposon created the paralogous bucentaur gene (bcntp97) in the ancestral ruminant

Retrotransposable element-1 (RTE-1) is a class of long interspersed nucleotide elements that contain in its open reading frame an apurinic/apyrimidinic endonuclease domain (AP-END) and a reverse transcriptase domain. Ruminants have a clade-specific RTE-1 (BovB/RTE). The bovine bcnt gene (bucentaur or craniofacial developmental protein 1) has a duplicated paralog (bcntp97) in tandem that recruited an AP-END of BovB/RTE as a coding exon (RTE exon). We obtained sequence of the bcnt region from several animals and showed that other ruminants also have the bcntp97 with a conserved RTE exon while camels and pigs do not. Genomic Southern analysis showed that camels and pigs have multiple bcnt-related sequences but not BovB/RTE which bovines and lesser mouse deer have abundantly. These results indicate that the bcnt gene duplication followed by the creation of bcntp97 including recruitment of the RTE exon occurred in the ancestral ruminant about 55 MYA. The indication of time frame is supported by a phylogenetic analysis. Taken together with a result of differential tissue expression of the two bcnt paralogs, we conclude that bcntp97 was created concurrently with the early radiation of BovB/RTE in an ancestral ruminant and then acquired a novel function.

Expression patterns of duplicate genes in the developing root in Arabidopsis thaliana

Data on gene expression in the development of the root in Arabidopsis thaliana were used to test for expression profile differences among multi-gene families and to examine the extent to which expression differences accompanied coding sequences divergence within families. Significant differences among families were observed on two principal axes, accounting for over 80% of the variance in the expression data. The number of synonymous nucleotide substitutions per synonymous site (d(S)) and the number of nonsynonymous nucleotide substitutions per nonsynonymous site (d(N)) were estimated between the members of two-member families (N = 428) and between phylogenetically independent sister pairs (N = 190) of sequences within larger families. Ribosomal proteins and a few other proteins were exceptional in showing highly divergent expression patterns in spite of very low levels of amino acid sequence divergence, as indicated by the low d(N) relative to d(S). However, the majority of gene duplicates showed relatively high levels of amino acid sequence divergence without appreciable change in expression pattern in the cell types analyzed.

Evolution of cow nonstomach lysozyme genes

Expansion of the lysozyme gene family is associated with the evolution of the ruminant lifestyle in ruminant artiodactyls such as the cow. Gene duplications allowed recombination between stomach lysozyme genes that may have assisted in the evolution of an enzyme adapted to survive and function in the stomach environment. Despite amplification of lysozyme genes, cow tears, milk, and blood are considered to be lysozyme deficient. Here we have identified 2 new cow lysozyme cDNA sequences and show that at least 4 different lysozymes are expressed in cows in nonstomach tissues and probably function as antibacterial defence enzymes. These 4 lysozyme genes are in addition to the 4 digestive lysozyme genes expressed in the stomach, yielding a number of expressed lysozyme genes in the cow larger than that found in most nonlysozyme-deficient mammals. In contrast to expectations, evidence for recombination between stomach and nonstomach lysozyme genes was found. Recombination, through concerted evolution, may have allowed some lysozymes to acquire the ability to survive in occasional acidic environments.Key words: molecular evolution, adaptive evolution, lysozyme, ruminants, gene duplication, gene evolution.

A novel Giraffidae-specific interspersed repeat with a microsatellite, originally found in an intron of a ruminant paralogous p97bcnt gene

The ruminant-specific p97bcnt gene (bcntp97) is a paralogous gene that includes a region derived from a retrotransposable element 1 (RTE-1). The region comprises an exon (RTE-1 exon) encoding 325 amino acids in the middle of the p97bcnt protein. To understand how the bcntp97 paralog evolved, we examined its organization in several ruminants. We found a 700-base pair (bp) insert in the 5' intron of the RTE-1 exon in giraffe bcntp97. This insert is missing in the corresponding regions of bovine and sika deer. Furthermore, the sequence of the insert is interspersed in the genome of giraffe but not bovine and also contains a (GA)n microsatellite. A highly homologous insert harboring significantly different (GA)n microsatellite was detected in the corresponding region of okapi bcntp97. Therefore, the interspersed fragments with (GA)n microsatellite might serve as a marker for tracking how duplicated genes evolve in a family-specific manner.

Human rDNA: Evolutionary Patterns within the Genes and Tandem Arrays Derived from Multiple Chromosomes

Human rDNA forms arrays on five chromosome pairs and is homogenized by concerted evolution through recombination and gene conversion (loci RNR1, RNR2, RNR3, RNR4, RNR5, OMIM: 180450). Homogenization is not perfect, however, so that it becomes possible to study its efficiency within genes, within arrays, and between arrays by measuring and comparing DNA sequence variation. Previous studies with randomly cloned genomic DNA fragments showed that different parts of the gene evolve at different rates but did not allow comparison of rDNA sequences derived from specific chromosomes. We have now cloned and sequenced rDNA fragments from specific acrocentric chromosomes to (1) study homogenization along the rDNA and (2) compare homogenization within chromosomes and between homologous and nonhomologous chromosomes. Our results show high homogeneity among regulatory and coding regions of rDNA on all chromosomes, a surprising homogeneity among adjacent distal non-rDNA sequences, and the existence of one to three very divergent intergenic spacer classes within each array.

Evolution of gene families

In eukaryote genomes, there are many kinds of gene families. Gene duplication and conversion are sources of the evolution of gene families, including those with uniform members and those with diverse functions. Population genetics theory on identity coefficients among gene members of a gene family shows that the balance between diversification by mutation, and homogenization by unequal crossing over and gene conversion, is important. Also, evolution of new functions is due to gene duplication followed by differentiation. Positive selection is necessary for the evolution of novel functions. However, many examples of current gene families suggest that both drift and selection are at work on their evolution.