Polymorphism, natural selection, and structural modeling of class Ia MHC in the African clawed frog (Xenopus laevis)

In the African clawed frog (Xenopus laevis), two deeply divergent allelic lineages of multiple genes of the class I MHC region have been discovered. For the MHC class I UAA locus, functional differences and the molecular basis for lineages maintenance are unknown. Alleles of linked class I region genes also exhibit strong disequilibrium with specific MHC alleles, but the underlying cause is not clear. We use MHC class Ia sequence data to estimate substitution rates and investigate structural differences between allelic lineages from protein models. Results indicate the operation of natural selection, and differences in the steric properties in the F pocket of the peptide-binding region among lineages. Variability in this pocket likely enables allelic lineages to bind very different sets of peptides and to interact differently with MHC chaperones in the endoplasmic reticulum. These results constitute evidence of the molecular evolutionary basis for 1) the maintenance of allelic lineages, 2) functional differences among lineages, and 3) strong linkage disequilibrium of allelic variants of class I region genes in X. laevis.Electronic Supplementary Material Supplementary material is available for this article at     

Evidence for convergent evolution in the antimicrobial peptide system in anuran amphibians

Amphibians are characterized in part by their highly specialized and glandular skin that enables key physiological functions such as cutaneous respiration and defense against a variety of micro- and macroscopic predators via toxic components (e.g., alkaloids and bufodienolids), biogenic amines, neuropeptides and antimicrobial peptides (AMPs). To date, DNA sequence information regarding AMP genes in anurans is restricted to only a few anuran families and largely to “higher frogs” (Neobatrachia). Here, we analyze the DNA information for the signal sequences of the AMP precursors in anuran amphibians available to the end of 2009 in an explicit phylogenetic framework to characterize the evolution of this large, diverse gene family. Comparison of cDNA sequences suggests that there are at least three different motifs within the signal peptide sequence of the AMP-precursor corresponding to the evolutionary lineages Neobatrachia, Bombinatoridae (Bombina spp.) and Pipidae (Xenopus laevis). The signal sequences are strongly conserved within each lineage (as previously noted for Neobatrachia), but highly divergent between them. Together with the lack of a linear relationship between the degree of sequence divergence and evolutionary time, we hypothesize that the anuran AMP system has evolved convergently on at least three occasions. However, additional sampling, especially among the largely poorly sampled non-neobatrachian lineages, is required to confirm this hypothesis and could reveal the existence of additional signal sequence motifs.     

Comparative gene expression analysis and fate mapping studies suggest an early segregation of cardiogenic lineages in Xenopus laevis

Retrospective clonal analysis in mice suggested that the vertebrate heart develops from two sources of cells called first and second lineages, respectively. Cells of the first lineage enter the linear heart tube and initiate terminal differentiation earlier than cells of the second lineage. It is thought that both heart lineages arise from a common progenitor cell population prior to the cardiac crescent stage (E7.5 of mouse development). The timing of segregation of different lineages as well as the molecular mechanisms underlying this process is not yet known. Furthermore, gene expression data for those lineages are very limited. Here we provide the first comparative study of cardiac marker gene expression during Xenopus laevis embryogenesis complemented by single cell RT-PCR analysis. In addition we provide fate mapping data of cardiac progenitor cells at different stages of development. Our analysis indicates an early segregation of cardiac lineages and a fairly complex heterogeneity of gene expression in the cardiac progenitor cells. Furthermore, this study sets a reference for all further studies analyzing cardiac development in X. laevis.     

Geographic independence and phylogenetic diversity of red shiner introductions

Identifying areas at risk of invasion can be difficult when the distribution of a non-native species encompasses geographically disjunct regions. Understanding genealogical relationships among native and non-native populations can clarify the origins of fragmented distributions, which in turn can clarify how fast and far a non-native species may spread. We evaluated genetic variation across the native and invasive ranges of red shiner (Cyprinella lutrensis), a minnow known to displace and hybridize with native species, to reconstruct invasion pathways across the United States (USA). Examination of mitochondrial cytochrome-b variation found that native range populations of red shiner fall into four highly divergent lineages that likely warrant species recognition. Introduced red shiner populations in the eastern and western USA are derived from only two of these lineages. Western USA populations originate from the mid-western and western genetic lineages, whereas eastern introductions derive only from the mid-western lineage. Western USA invasive populations exhibit fewer, but more diverse haplotypes compared to eastern USA invasive populations. We also recovered an undescribed, divergent lineage of Cyprinella that has been cryptically introduced into the western USA, which raises the possibility that hybridization has proceeded following secondary contact between previously allopatric lineages. Approximate Bayesian Computation modeling suggests that the disjunct distribution of red shiner across North America is an agglomeration of independent regional invasions with distinct origins, rather than stepwise advance of an invasion front or secondary introductions across regions. Thus localized control may be effective in managing non-native red shiner, including further spread to areas of conservation concern.     

Evolutionary History of the Morocco lizard-Fingered Geckos of the <Emphasis Type="Italic">Saurodactylus brosseti</Emphasis> Complex

Studies of biodiversity in the Maghreb have revealed high genetic diversity and divergent genetic lineages among many taxa including squamates. Geographic barriers such as the Atlas Mountains are one of the key factors promoting genetic differentiation and the high levels of endemism. The lizard-fingered gecko Saurodactylus brosseti is endemic to Morocco. Its range includes both sides of the Atlas Mountains, and although high diversity was previously identified within the species, much of the range was unsampled. To understand the evolutionary and biogeographical history of this species, we used mitochondrial and nuclear DNA sequence data from 64 populations sampled across the entire species range. We employed phylogenetic methods based on gene trees and species trees as well as a time calibrated Bayesian genealogy and coalescent species delimitation approaches. We uncovered four highly divergent and allopatric mitochondrial lineages that did not share haplotypes at variable nuclear loci, suggesting the four groups have been evolving independently since the Miocene, according to our molecular dating estimates. These results coupled with the geographic pattern of genetic diversity suggest a possible role of the Atlas Mountains for the divergence observed between the four lineages of S. brosseti, while each lineage probably later underwent several episodes of fragmentation followed by (re-) expansion during Pleistocene climatic oscillations. Bayesian species delimitation analysis indicates that the four lineages may well be distinct species but we suggest that detailed morphological analyses are needed prior to taxonomic changes. The four lineages represent ancient independent evolutionary units, and deserve conservation management as distinct entities.     

Duplication of the MHC-linked Xenopus complement factor B gene

We have previously reported the molecular cloning of the mammalian major histocompatibility complex (MHC) class III gene, complement factor B (Bf) from Xenopus laevis, and linkage of the gene to the frog MHC. Here, we estimated the copy number of the Xenopus Bf gene by genomic Southern blotting analysis and demonstrated that Xenopus laevis has two copies of the Bf gene. Both genes co-segregated with the MHC-linked HSP70 genes among 19 offspring of an f/r × f/r cross, indicating a close linkage of the two Bf genes to the frog MHC. Both genes are transcribed and contain open reading frames. When compared with the previously determined cDNA sequence (Xenopus Bf A), the predicted amino acid sequence of the second cDNA species (Xenopus Bf B) shows 82% overall identity. Polymerase chain reaction analysis indicated that all of the partially inbred frogs with the f, r, g, and j MHC haplotypes, as well as 12 outbred frogs tested have both Bf genes, suggesting that the duplicated Bf genes are stable genetic traits in Xenopus laevis.     

Prediction of Functionally Important Phospho-Regulatory Events in Xenopus laevis Oocytes

The African clawed frog Xenopus laevis is an important model organism for studies in developmental and cell biology, including cell-signaling. However, our knowledge of X. laevis protein post-translational modifications remains scarce. Here, we used a mass spectrometry-based approach to survey the phosphoproteome of this species, compiling a list of 2636 phosphosites. We used structural information and phosphoproteomic data for 13 other species in order to predict functionally important phospho-regulatory events. We found that the degree of conservation of phosphosites across species is predictive of sites with known molecular function. In addition, we predicted kinase-protein interactions for a set of cell-cycle kinases across all species. The degree of conservation of kinase-protein interactions was found to be predictive of functionally relevant regulatory interactions. Finally, using comparative protein structure models, we find that phosphosites within structured domains tend to be located at positions with high conformational flexibility. Our analysis suggests that a small class of phosphosites occurs in positions that have the potential to regulate protein conformation.     

Comprehensive analysis of medaka major histocompatibility complex (MHC) class II genes: Implications for evolution in teleosts

The major histocompatibility complex (MHC) class II molecules play central roles in adaptive immunity by regulating immune response via the activation of CD4 T cells. The full complement of the MHC class II genes has been elucidated only in mammalian species to date. To understand the evolution of these genes, we performed their first comprehensive analysis in nonmammalian species using a teleost, medaka (Oryzias latipes). Based on a database search, cDNA cloning, and genomic PCR, medaka was shown to possess five pairs of expressed class II genes, comprising one IIA and one IIB gene. Each pair was located on a different chromosome and was not linked to the class I genes. Only one pair showed a high degree of polymorphism and was considered to be classical class II genes, whereas the other four pairs were nonclassical. Phylogenetic analysis of all medaka class II genes and most reported teleost class II genes revealed that the IIA and IIB genes formed separate clades, each containing three well-corresponding lineages. One lineage contained three medaka genes and all known classical class II genes of Ostariophysi and Euteleostei and was presumed to be an original lineage of the teleost MHC class II genes. The other two lineages contained one nonclassical medaka gene each and some Euteleostei genes. These results indicate that multiple lineages of the teleost MHC class II genes have been conserved for hundreds of millions of years and that the tightly linked IIA and IIB genes have undergone concerted evolution.     

Revisiting phylogenetic diversity and cryptic species of <Emphasis Type="Italic">Cenococcum geophilum</Emphasis> sensu lato

The fungus Cenococcum geophilum Fr. (Dothideomycetes, Ascomycota) is one of the most common ectomycorrhizal fungi in boreal to temperate regions. A series of molecular studies has demonstrated that C. geophilum is monophyletic but a heterogeneous species or a species complex. Here, we revisit the phylogenetic diversity of C. geophilum sensu lato from a regional to intercontinental scale by using new data from Florida (USA) along with existing data in GenBank from Japan, Europe, and North America. The combination of internal transcribed spacer (ITS) ribosomal DNA and the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene resolved six well-supported lineages (87–100 % bootstrap values) that are closely related to each other and a seventh lineage that is phylogenetically distinct. A multi-locus analysis (small subunit (SSU), large subunit (LSU), translational elongation factor (TEF), and the largest and second-largest subunits of RNA polymerase II (RPB1 and RPB2)) revealed that the divergent lineage is the sister group to all other known Cenococcum isolates. Isolates of the divergent lineage grow fast on nutrient media and do not form ectomycorrhizas on seedlings of several pine and oak species. Our results indicate that C. geophilum sensu lato includes more phylogenetically distinct cryptic species than have previously been reported. Furthermore, the divergent lineage appears to be a non-mycorrhizal sister group. We discuss the phylogenetic diversity of C. geophilum sensu lato and argue in favor of species recognition based on phylogenetic and ecological information in addition to morphological characteristics. A new genus and species (Pseudocenococcum floridanum gen. et sp. nov.) is proposed to accommodate a divergent and putatively non-mycorrhizal lineage.     

Cryptic genetic diversity and spatial patterns of admixture within Belizean marine reserves

Despite the potential for long-distance gene flow in the sea, there is growing evidence of cryptic genetic diversity in many marine taxa. Understanding the geographic distribution of cryptic lineages, as well as the spatial patterns of admixture among them, can have important implications for conservation planning. Here, we explore patterns of divergence in a coral reef fish, the neon goby Elacatinus lori, across the species’ range. First, we use targeted amplicon sequencing to describe the spatial pattern of genetic divergence using two marker types (57 anonymous ddRAD-derived loci and mtDNA cytb). Second, we quantify the degree of admixture and hybridization between two previously-unidentified divergent lineages within Belize. Third, we assess whether the existing group of marine protected areas (MPAs) in Belize protects this cryptic genetic diversity. The results provide strong evidence for two divergent genetic lineages of E. lori within Belize, separated geographically by only 30 km of low-suitability habitat. There is a sharp genetic cline across these 30 km, and evidence of admixture and introgression at the boundary regions of the habitat break. We also show that the broadly-distributed arrangement of MPAs within Belize protects both major lineages as well as subtle structure within-lineages, and therefore may confer protection to co-distributed species that exhibit similar spatial patterns of divergence.     

Widespread Presence of Human BOULE Homologs among Animals and Conservation of Their Ancient Reproductive Function

Sex-specific traits that lead to the production of dimorphic gametes, sperm in males and eggs in females, are fundamental for sexual reproduction and accordingly widespread among animals. Yet the sex-biased genes that underlie these sex-specific traits are under strong selective pressure, and as a result of adaptive evolution they often become divergent. Indeed out of hundreds of male or female fertility genes identified in diverse organisms, only a very small number of them are implicated specifically in reproduction in more than one lineage. Few genes have exhibited a sex-biased, reproductive-specific requirement beyond a given phylum, raising the question of whether any sex-specific gametogenesis factors could be conserved and whether gametogenesis might have evolved multiple times. Here we describe a metazoan origin of a conserved human reproductive protein, BOULE, and its prevalence from primitive basal metazoans to chordates. We found that BOULE homologs are present in the genomes of representative species of each of the major lineages of metazoans and exhibit reproductive-specific expression in all species examined, with a preponderance of male-biased expression. Examination of Boule evolution within insect and mammalian lineages revealed little evidence for accelerated evolution, unlike most reproductive genes. Instead, purifying selection was the major force behind Boule evolution. Furthermore, loss of function of mammalian Boule resulted in male-specific infertility and a global arrest of sperm development remarkably similar to the phenotype in an insect boule mutation. This work demonstrates the conservation of a reproductive protein throughout eumetazoa, its predominant testis-biased expression in diverse bilaterian species, and conservation of a male gametogenic requirement in mice. This shows an ancient gametogenesis requirement for Boule among Bilateria and supports a model of a common origin of spermatogenesis.     

Comparative morphological analysis of two parallel mycoheterotrophic transitions reveals divergent and convergent traits in the genus <Emphasis Type="Italic">Pyrola</Emphasis> (Pyroleae,Ericaceae)

The genus Pyrola includes species with different degree of mycoheterotrophy; some species possess individuals that rely on all carbon through their associations with fungi (full mycoheterotrophy, FM), whereas some species obtain carbon through both fungi and photosynthesis by itself (partial mycoheterotrophy, PM). To investigate how plant functional traits of photosynthesis and reproduction are related to the degree of mycoheterotrophy in the initial stage of the transition from PM to FM, we determined morphological traits in FM (or nearly FM) and PM species in two independent lineages, P. picta and P. japonica complexes. We used herbarium specimens and examined leaf number, leaf area, flower number, and scape length in FM or nearly FM species (P. aphylla and P. subaphylla) and PM species (P. picta s.l. and P. japonica). We found a leaf area reduction in FM (or nearly FM) species in both lineages, suggesting that this is a convergent trait. The number of flowers was not significantly different between FM (or nearly FM) and PM species in both lineages. On the other hand, differences in the variation between FM (or nearly FM) and PM species were found in some traits between the two lineages. The FM (or nearly FM) species in one lineage only possessed rudimentary leaves, whereas that in the other linage possessed a few small, ordinary leaves in addition to those with only rudimentary leaves. The scape length of the FM (or nearly FM) species was significantly longer than that of PM species in one lineage, whereas it was shorter in the other lineage. The different and common variations are divergent and convergent traits, respectively, that could be associated with the transition to FM in Pylora. In addition, shoots of both PM species occasionally lacked ordinary leaves, possibly indicating possession of these shoots is preadaptation for the transition to FM in Pyrola.     

The hidden diversity of Coleodactylus amazonicus (Sphaerodactylinae,Gekkota) revealed by molecular data

Coleodactylus amazonicus, a small leaf-litter diurnal gecko widely distributed in Amazon Basin, has been considered a single species with no significant morphological differences between populations along its range. A recent molecular study, however, detected large genetic differences between populations of central Amazonia and those in the easternmost part of the Amazon Basin, suggesting the presence of taxonomically unrecognised diversity. In this study, DNA sequences of three mitochondrial (16S, cytb, and ND4) and two nuclear genes (RAG-1, c-mos) were used to investigate whether the species currently identified as C.amazonicus contains morphologically cryptic species lineages. The present phylogenetic analysis reveals further genetic subdivision including at least five potential species lineages, restricted to northeastern (lineage A), southeastern (lineage B), central-northern (lineage E) and central-southern (lineages C and D) parts of Amazon Basin. All clades are characterized by exclusive groups of alleles for both nuclear genes and highly divergent mitochondrial haplotype clades, with corrected pairwise net sequence divergence between sister lineages ranging from 9.1% to 20.7% for the entire mtDNA dataset. Results of this study suggest that the real diversity of “C. amazonicus” has been underestimated due to its apparent cryptic diversification.     

Regulation of TGF-β receptor activity

Background

Thyroid hormone (T3) is important for adult organ function and vertebrate development. Amphibian metamorphosis is totally dependent on T3 and offers a unique opportunity to study how T3 controls postembryonic development in vertebrates. Earlier studies have demonstrated that TR mediates the metamorphic effects of T3 in Xenopus laevis. Liganded TR recruits histone modifying coactivator complexes to target genes during metamorphosis. This leads to nucleosomal removal and histone modifications, including methylation of histone H3 lysine (K) 79, in the promoter regions, and the activation of T3-inducible genes.

Results

We show that Dot1L, the only histone methyltransferase capable of methylating H3K79, is directly regulated by TR via binding to a T3 response element in the promoter region during metamorphosis in Xenopus tropicalis, a highly related species of Xenopus laevis. We further show that Dot1L expression in both the intestine and tail correlates with the transformation of the organs.

Conclusions

Our findings suggest that TR activates Dot1L, which in turn participates in metamorphosis through a positive feedback to enhance H3K79 methylation and gene activation by liganded TR.     

Homoeologous chromosomes of Xenopus laevis are highly conserved after whole-genome duplication

It has been suggested that whole-genome duplication (WGD) occurred twice during the evolutionary process of vertebrates around 450 and 500 million years ago, which contributed to an increase in the genomic and phenotypic complexities of vertebrates. However, little is still known about the evolutionary process of homoeologous chromosomes after WGD because many duplicate genes have been lost. Therefore, Xenopus laevis (2n=36) and Xenopus (Silurana) tropicalis (2n=20) are good animal models for studying the process of genomic and chromosomal reorganization after WGD because X. laevis is an allotetraploid species that resulted from WGD after the interspecific hybridization of diploid species closely related to X. tropicalis. We constructed a comparative cytogenetic map of X. laevis using 60 complimentary DNA clones that covered the entire chromosomal regions of 10 pairs of X. tropicalis chromosomes. We consequently identified all nine homoeologous chromosome groups of X. laevis. Hybridization signals on two pairs of X. laevis homoeologous chromosomes were detected for 50 of 60 (83%) genes, and the genetic linkage is highly conserved between X. tropicalis and X. laevis chromosomes except for one fusion and one inversion and also between X. laevis homoeologous chromosomes except for two inversions. These results indicate that the loss of duplicated genes and inter- and/or intrachromosomal rearrangements occurred much less frequently in this lineage, suggesting that these events were not essential for diploidization of the allotetraploid genome in X. laevis after WGD.     

Multiple divergent haplotypes express completely distinct sets of class I MHC genes in zebrafish

The zebrafish is an important animal model for stem cell biology, cancer, and immunology research. Histocompatibility represents a key intersection of these disciplines; however, histocompatibility in zebrafish remains poorly understood. We examined a set of diverse zebrafish class I major histocompatibility complex (MHC) genes that segregate with specific haplotypes at chromosome 19, and for which donor-recipient matching has been shown to improve engraftment after hematopoietic transplantation. Using flanking gene polymorphisms, we identified six distinct chromosome 19 haplotypes. We describe several novel class I U lineage genes and characterize their sequence properties, expression, and haplotype distribution. Altogether, ten full-length zebrafish class I genes were analyzed, mhc1uba through mhc1uka. Expression data and sequence properties indicate that most are candidate classical genes. Several substitutions in putative peptide anchor residues, often shared with deduced MHC molecules from additional teleost species, suggest flexibility in antigen binding. All ten zebrafish class I genes were uniquely assigned among the six haplotypes, with dominant or codominant expression of one to three genes per haplotype. Interestingly, while the divergent MHC haplotypes display variable gene copy number and content, the different genes appear to have ancient origin, with extremely high levels of sequence diversity. Furthermore, haplotype variability extends beyond the MHC genes to include divergent forms of psmb8. The many disparate haplotypes at this locus therefore represent a remarkable form of genomic region configuration polymorphism. Defining the functional MHC genes within these divergent class I haplotypes in zebrafish will provide an important foundation for future studies in immunology and transplantation.     

Evolution of Chemical Diversity in Echinocandin Lipopeptide Antifungal Metabolites

The echinocandins are a class of antifungal drugs that includes caspofungin, micafungin, and anidulafungin. Gene clusters encoding most of the structural complexity of the echinocandins provided a framework for hypotheses about the evolutionary history and chemical logic of echinocandin biosynthesis. Gene orthologs among echinocandin-producing fungi were identified. Pathway genes, including the nonribosomal peptide synthetases (NRPSs), were analyzed phylogenetically to address the hypothesis that these pathways represent descent from a common ancestor. The clusters share cooperative gene contents and linkages among the different strains. Individual pathway genes analyzed in the context of similar genes formed unique echinocandin-exclusive phylogenetic lineages. The echinocandin NRPSs, along with the NRPS from the inp gene cluster in Aspergillus nidulans and its orthologs, comprise a novel lineage among fungal NRPSs. NRPS adenylation domains from different species exhibited a one-to-one correspondence between modules and amino acid specificity that is consistent with models of tandem duplication and subfunctionalization. Pathway gene trees and Ascomycota phylogenies are congruent and consistent with the hypothesis that the echinocandin gene clusters have a common origin. The disjunct Eurotiomycete-Leotiomycete distribution appears to be consistent with a scenario of vertical descent accompanied by incomplete lineage sorting and loss of the clusters from most lineages of the Ascomycota. We present evidence for a single evolutionary origin of the echinocandin family of gene clusters and a progression of structural diversification in two fungal classes that diverged approximately 290 to 390 million years ago. Lineage-specific gene cluster evolution driven by selection of new chemotypes contributed to diversification of the molecular functionalities.