Search

Background: Secretory carrier membrane proteins (SCAMPs) form a family of integral membrane proteins and play a crucial role in mediating exocytosis in both animals and plants. While SCAMP genes have been studied in several plant species, their functions in cotton, particularly in response to abiotic stress, have not yet been reported., Results: In this study, a total of 53 SCAMP genes were identified in G. arboreum, G. raimondii, G. hirsutum, and G. barbadense. These genes were classified into five groups based on a phylogenetic analysis with SCAMPs from Arabidopsis thaliana. The main factor driving the expansion of the SCAMP gene family in G. hirsutum is tandem and segmental duplication events. Using MEME, in addition to the conserved SCAMP domain, we identified 3-13 other domains in each GhSCAMP. The cis-element analysis suggested that GhSCAMPs were widely involved in cotton growth and development, and responses to abiotic stresses. RNA sequencing (RNA-Seq) and reverse transcription-quantitative polymerase chain reaction (RT-qPCR) results showed that most GhSCAMPs were expressed highly in many tissues and had differential expression responses to drought, cold, and heat stresses. Knock-down of GhSCAMP2 and GhSCAMP4 by virus-induced gene silencing (VIGS) lead to a salt-sensitive phenotype and had a lower content of CAT, POD, and SOD., Conclusions: This study identified SCAMP genes in four cotton species, enhancing our understanding of the potential biological functions of SCAMPs. Additionally, we demonstrated that GhSCAMP2 and GhSCAMP4 positively regulate cotton tolerance to salt stress., (© 2024. The Author(s).)

Background: Phenylalanine ammonia-lyase (PAL) serves as a key gateway enzyme, bridging primary metabolism and the phenylpropanoid pathway, and thus playing an indispensable role in flavonoid, anthocyanin and lignin biosynthesis. PAL gene families have been extensively studied across species using public genomes. However, a comprehensive exploration of PAL genes in Epimedium species, especially those involved in prenylated flavonol glycoside, anthocyanin, or lignin biosynthesis, is still lacking. Moreover, an in-depth investigation into PAL gene family evolution is warranted., Results: Seven PAL genes (EpPAL1-EpPAL7) were identified. EpPAL2 and EpPAL3 exhibit low sequence identity to other EpPALs (ranging from 61.09 to 64.38%) and contain two unique introns, indicating distinct evolutionary origins. They evolve at a rate ~ 10 to ~ 54 times slower compared to EpPAL1 and EpPAL4-7, suggesting strong purifying selection. EpPAL1 evolved independently and is another ancestral gene. EpPAL1 formed EpPAL4 through segmental duplication, which lead to EpPAL5 and EpPAL6 through tandem duplications, and EpPAL7 through transposed duplication, shaping modern EpPALs. Correlation analysis suggests EpPAL1, EpPAL2 and EpPAL3 play important roles in prenylated flavonol glycosides biosynthesis, with EpPAL2 and EpPAL3 strongly correlated with both Epimedin C and total prenylated flavonol glycosides. EpPAL1, EpPAL2 and EpPAL3 may play a role in anthocyanin biosynthesis in leaves. EpPAL2, EpPAL3, EpPAL6, and EpPAL7 might be engaged in anthocyanin production in petals, and EpPAL2 and EpPAL3 might also contribute to anthocyanin synthesis in sepals. Further experiments are needed to confirm these hypotheses. Novel insights into the evolution of PAL gene family suggest that it might have evolved from a monophyletic group in bryophytes to large-scale sequence differentiation in gymnosperms, basal angiosperms, and Magnoliidae. Ancestral gene duplications and vertical inheritance from gymnosperms to angiosperms likely occurred during PAL evolution. Most early-diverging eudicotyledons and monocotyledons have distinct histories, while modern angiosperm PAL gene families share similar patterns and lack distant gene types., Conclusions: EpPAL2 and EpPAL3 may play crucial roles in biosynthesis of prenylated flavonol glycosides and anthocyanins in leaves and flowers. This study provides novel insights into PAL gene family evolution. The findings on PAL genes in E. pubescens will aid in synthetic biology research on prenylated flavonol glycosides production., (© 2024. The Author(s).)

Objective: The eukaryotic tree of life has been subject of numerous studies ever since the nineteenth century, with more supergroups and their sister relations being decoded in the last years. In this study, we reconstructed the phylogeny of eukaryotes using complete 18S rDNA sequences and their individual secondary structures simultaneously. After the sequence-structure data was encoded, it was automatically aligned and analyzed using sequence-only as well as sequence-structure approaches. We present overall neighbor-joining trees of 211 eukaryotes as well as the respective profile neighbor-joining trees, which helped to resolve the basal branching pattern. A manually chosen subset was further inspected using neighbor-joining, maximum parsimony, and maximum likelihood analyses. Additionally, the 75 and 100 percent consensus structures of the subset were predicted. Results: All sequence-structure approaches show improvements compared to the respective sequence-only approaches: the average bootstrap support per node of the sequence-structure profile neighbor-joining analyses with 90.3, was higher than the average bootstrap support of the sequence-only profile neighbor-joining analysis with 73.9. Also, the subset analyses using sequence-structure data were better supported. Furthermore, more subgroups of the supergroups were recovered as monophyletic and sister group relations were much more comparable to results as obtained by multi-marker analyses. [ABSTRACT FROM AUTHOR]

The composition and diversity of animal-associated microbial communities are shaped by multiple ecological and evolutionary processes acting at different spatial and temporal scales. Skin microbiomes are thought to be strongly influenced by the environment due to the direct interaction of the host's skin with the external media. As expected, the diversity of amphibian skin microbiomes is shaped by climate and host sampling habitats, whereas phylogenetic effects appear to be weak. However, the relative strength of phylogenetic and environmental effects on salamander skin microbiomes remains poorly understood. Here, we analysed sequence data from 1164 adult salamanders of 44 species to characterise and compare the diversity and composition of skin bacteria. We assessed the relative contribution of climate, host sampling habitat, and host phylogeny to the observed patterns of bacterial diversity. We found that bacterial alpha diversity was mainly associated with host sampling habitat and climate, but that bacterial beta diversity was more strongly associated with host taxonomy and phylogeny. This phylogenetic effect predominantly occurred at intermediate levels of host divergence (0–50 Mya). Our results support the importance of environmental factors shaping the diversity of salamander skin microbiota, but also support host phylogenetic history as a major factor shaping these bacterial communities. [ABSTRACT FROM AUTHOR]

Background: The receptor-like kinase (RLK) gene families in plants contains a large number of members. They are membrane proteins with an extracellular receptor domain and participate in biotic and abiotic stress responses. Results: In this study, we identified RLKs in 15 representative plant genomes, including wheat, and classified them into 64 subfamilies by using four types of phylogenetic trees and HMM models. Conserved exon‒intron structures with conserved exon phases in the kinase domain were found in many RLK subfamilies from Physcomitrella patens to Triticum aestivum. Domain distributions of RLKs were also diagrammed. Collinearity events and tandem gene clusters suggested that polyploidization and tandem duplication events contributed to the member expansions of T. aestivum RLKs. Global expression pattern analysis was performed by using public transcriptome data. These analyses were involved in T. aestivum, Aegilops tauschii and Brachypodium distachyon RLKs under biotic and abiotic stresses. We also selected 9 RLKs to validate the transcriptome prediction by using qRT‒PCR under drought treatment and with Fusarium graminearum infection. The expression trends of these 9 wheat RLKs from public transcriptome data were consistent with the results of qRT‒PCR, indicating that they might be stress response genes under drought or F. graminearum treatments. Conclusion: In this study, we identified, classified, evolved, and expressed RLKs in wheat and related plants. Thus, our results will provide insights into the evolutionary history and molecular mechanisms of wheat RLKs. [ABSTRACT FROM AUTHOR]

Background: Sargassum polycystum C. Agardh and Sargassum plagiophyllum C. Agardh are inhabitants of tropical coastal areas, their populations are negatively influenced by global warming and marine environment changes. The mitochondrial and chloroplast genomes of these species have not been sequenced. Results: The mitochondrial genomes of S. polycystum and S. plagiophyllum were 34,825 bp and 34,862 bp, respectively, and their corresponding chloroplast genomes were 124,493 bp and 124,536 bp, respectively. The mitochondrial and chloroplast genomes of these species share conserved synteny, sequence regions and gene number when compared with the organellar genomes of other Sargassum species. Based on sequence analysis of 35 protein-coding genes, we deduced that S. polycystum and S. plagiophyllum were closely related with S. ilicifolium; these species diverged approximately 0.3 million years ago (Ma; 0.1–0.53 Ma) during the Pleistocene period (0.01–2.59 Ma). Rates of synonymous and non-synonymous substitutions in the mitochondrial genome of the Sargassum genus were 3 times higher than those in the chloroplast genome. In the mitochondrial genome, rpl5, rpl31 and rps11 had the highest synonymous substitution rates. In the chloroplast genome, psaE, rpl14 and rpl27 had the highest synonymous substitution rates. Conclusions: Phylogenetic analysis confirms the close relationship between the two sequenced species and S. ilicifolium. Both synonymous and non-synonymous substitution rates show significant divergence between the group of mitochondrial genomes versus the group of chloroplast genomes. The deciphering of complete mitochondrial and chloroplast genomes is significant as it advances our understanding of the evolutionary and phylogenetic relationships between species of brown seaweeds. [ABSTRACT FROM AUTHOR]

Background: Walnuts (Juglans regia L.) are known for their nutrient-rich nuts and are one of the important economic tree species in the world. However, due to global warming and soil salinization, walnuts suffer from various abiotic stresses. TIFY (TIF[F/Y]XG) proteins play an essential role in the growth and development of plants, signal transduction, and stress response in plants. At present, although the TIFY gene family of a number of plants has been identified and studied, how TIFY takes part in stress tolerance remains obscure and many functions of TIFY require further investigation. Result: In this study, twenty-one TIFY transcription factors were identified in the walnut genome database, and they were divided into four subfamilies (TIFY, JAZ, ZML, and PPD) by bioinformatics analysis. Chromosome location revealed tandem duplication of some genes. Phylogenetic tree analysis showed JrTIFYs were closely related to the TIFY gene family of Arabidopsis thaliana (A. thaliana). qRT-PCR (quantitative real-time PCR) analysis revealed the TIFY genes have different expression patterns in 'Qingxiang' and 'Xiangling' walnut varieties under drought, heat, and salt stress. JAZ subfamily was more expressed in different abiotic stress than other subfamilies. The expressions of JrTIFY14 under heat and salt stress were significantly higher than those under drought stress. However, the expression of JrTIFYs was not significant in 'Xiangling'. Conclusion: This study reveals the TIFY gene family plays an important role in walnuts facing abiotic stresses and provides a theoretical basis for walnut breeding. [ABSTRACT FROM AUTHOR]

Background: Heat shock proteins (HSPs) perform a fundamental role in protecting plants against abiotic stresses. Individual family members have been analyzed in previous studies, but there has not yet been a comprehensive analysis of the HSP70 gene family in Pyropia yezoensis. Results: We investigated 15 putative HSP70 genes in Py. yezoensis. These genes were classified into two sub-families, denoted as DnaK and Hsp110. In each sub-family, there was relative conservation of the gene structure and motif. Synteny-based analysis indicated that seven and three PyyHSP70 genes were orthologous to HSP70 genes in Pyropia haitanensis and Porphyra umbilicalis, respectively. Most PyyHSP70s showed up-regulated expression under different degrees of dehydration stress. PyyHSP70-1 and PyyHSP70-3 were expressed in higher degrees compared with other PyyHSP70s in dehydration treatments, and then expression degrees somewhat decreased in rehydration treatment. Subcellular localization showed PyyHSP70-1-GFP and PyyHSP70-3-GFP were in the cytoplasm and nucleus/cytoplasm, respectively. Similar expression patterns of paired orthologs in Py. yezoensis and Py. haitanensis suggest important roles for HSP70s in intertidal environmental adaptation during evolution. Conclusions: These findings provide insight into the evolution and modification of the PyyHSP70 gene family and will help to determine the functions of the HSP70 genes in Py. yezoensis growth and development. [ABSTRACT FROM AUTHOR]

Background: The mitochondrial (mt) genomes of Sarcoptiformes mites typically contain 37 genes. Although the loss of genes is rare in Sarcoptiformes mite mitogenomes, two of the six previously reported oribatid mites (Acariforms: Sarcoptiformes) are reported to have lost parts of their tRNA genes. To confirm whether the tRNA genes were indeed lost and whether the loss is universal, we re-annotated the available oribatid mite sequences and sequenced the mitogenome of Oribatula sakamorii. Methods: The mitogenome of O. sakamorii was sequenced using an Illumina HiSeq sequencer. The mt tRNA gene was annotated using multi-software combined with a manual annotation approach. Phylogenetic analyses were performed using the maximum likelihood and Bayesian inference methods with concatenated nucleotide and amino acid sequences. Results: The mitogenomes of O. sakamorii contained 37 genes, including 22 tRNA genes. We identified all mt tRNA genes that were reported as "lost" in Steganacarus magnus and Paraleius leontonychus and revealed certain atypical tRNA annotation errors in oribatid mite sequences. Oribatid mite mitogenomes are characterized by low rates of genetic rearrangement, with six or seven gene blocks conserved between the mitogenome of all species and that of ancestral arthropods. Considering the relative order of the major genes (protein-coding genes and rRNAs), only one or two genes were rearranged with respect to their positions in the ancestral genome. We explored the phylogenetic relationships among the available oribatid mites, and the results confirmed the systematic position of Hermannia in the Crotonioidea superfamily. This was also supported by the synapomorphic gene-derived boundaries. Conclusions: The tRNA "lost" phenomenon is not universal in oribatid mites. Rather, highly atypical secondary structure of the inferred mt tRNA genes made them unidentifiable using a single type of tRNA search program. The use of multi-software combined with a manual annotation approach can improve the accuracy of tRNA gene annotation. In addition, we identified the precise systematic position of Hermannia and validated that Astigmata is nested in Oribatida. [ABSTRACT FROM AUTHOR]

Background: The genus Brassica mainly comprises three diploid and three recently derived allotetraploid species, most of which are highly important vegetable, oil or ornamental crops cultivated worldwide. Despite being extensively studied, the origination of B. napus and certain detailed interspecific relationships within Brassica genus remains undetermined and somewhere confused. In the current high-throughput sequencing era, a systemic comparative genomic study based on a large population is necessary and would be crucial to resolve these questions. Results: The chloroplast DNA and mitochondrial DNA were synchronously resequenced in a selected set of Brassica materials, which contain 72 accessions and maximally integrated the known Brassica species. The Brassica genomewide cpDNA and mtDNA variations have been identified. Detailed phylogenetic relationships inside and around Brassica genus have been delineated by the cpDNA- and mtDNA- variation derived phylogenies. Different from B. juncea and B. carinata, the natural B. napus contains three major cytoplasmic haplotypes: the cam-type which directly inherited from B. rapa, polima-type which is close to cam-type as a sister, and the mysterious but predominant nap-type. Certain sparse C-genome wild species might have primarily contributed the nap-type cytoplasm and the corresponding C subgenome to B. napus, implied by their con-clustering in both phylogenies. The strictly concurrent inheritance of mtDNA and cpDNA were dramatically disturbed in the B. napus cytoplasmic male sterile lines (e.g., mori and nsa). The genera Raphanus, Sinapis, Eruca, Moricandia show a strong parallel evolutional relationships with Brassica. Conclusions: The overall variation data and elaborated phylogenetic relationships provide further insights into genetic understanding of Brassica, which can substantially facilitate the development of novel Brassica germplasms. [ABSTRACT FROM AUTHOR]

Background: Herpes simplex viruses form a genus within the alphaherpesvirus subfamily, with three identified viral species isolated from Old World monkeys (OWM); Macacine alphaherpesvirus 1 (McHV-1; herpes B), Cercopithecine alphaherpesvirus 2 (SA8), and Papiine alphaherpesvirus 2 (PaHV-2; herpes papio). Herpes B is endemic to macaques, while PaHV-2 and SA8 appear endemic to baboons. All three viruses are genetically and antigenically similar, with SA8 and PaHV-2 thought to be avirulent in humans, while herpes B is a biosafety level 4 pathogen. Recently, next-generation sequencing (NGS) has resulted in an increased number of published OWM herpes simplex genomes, allowing an encompassing phylogenetic analysis. Results: In this study, phylogenetic networks, in conjunction with a genome-based genetic distance cutoff method were used to examine 27 OWM monkey herpes simplex isolates. Genome-based genetic distances were calculated, resulting in distances between lion and pig-tailed simplex viruses themselves, and versus herpes B core strains that were higher than those between PaHV-2 and SA8 (approximately 14 and 10% respectively). The species distance cutoff was determined to be 8.94%, with the method recovering separate species status for PaHV-2 and SA8 and showed that lion and pig-tailed simplex viruses (vs core herpes B strains) were well over the distance species cutoff. Conclusions: We propose designating lion and pig-tailed simplex viruses as separate, individual viral species, and that this may be the first identification of viral cryptic species. [ABSTRACT FROM AUTHOR]

Backgrounds: Cyclic nucleotide gated channels (CNGCs) play multifaceted roles in plant physiological processes, especially with respect to signalling processes, plant development, and responses to environmental stresses. However, little information is known about the CNGC family in the large cosmopolitan family Rhamnaceae, which has strong tolerance to biotic and abiotic stresses. Results: In the current study, a total of 15 ZjCNGCs which located on 7 chromosomes were firstly identified in Chinese jujube (Ziziphus jujuba Mill.), the most important species of Rhamnaceae in terms of economic and ecological values. Phylogenetic analysis showed that these ZjCNGCs could be classified into four groups, ZjCNGC12 belonged to group IVA, and ZjCNGC13, 14, 15 belonged to group IVB. In addition, the paralogous and orthologous homology duplication of ZjCNGC15 occurred during the evolutionary process. The characteristics of ZjCNGCs regarding to exon-intron numbers and post-translational modifications showed diversified structures and functions. Motif composition and protein sequence analysis revealed that the phosphate-binding cassette and hinge regions were conserved among ZjCNGCs. Prediction of the cis-acting regulatory elements and expression profiles by real-time quantitative PCR analysis showed that some of the ZjCNGCs responded to environmental changes, especially ZjCNGC2, which was significantly downregulated in response to cold stress, and ZjCNGC4 was highly induced in response to cold, salt and alkaline stresses. ZjCNGC13 and 14 were highly induced in the phytoplasma-resistant cultivar and downregulated in the susceptible cultivar. Furthermore, ZjCNGC2 could be regulated by cAMP treatment, microtubule changes and interact with ZjMAPKK4, which suggested that cAMP and microtubule might play important roles in ZjCNGC2 mediated ZjMAPKK4 signalling transduction involved in cold stress. Conclusions: The identification and classification analysis of ZjCNGCs were firstly reported, and some key individual ZjCNGCs might play essential roles in the response to biotic and abiotic stresses, especially ZjCNGC2 mediated ZjMAPKK4 signalling transduction involved in cold stress. This systematic analysis could provide important information for further functional characterization of ZjCNGCs with the aim of breeding stress-resistant cultivars. [ABSTRACT FROM AUTHOR]

Background: Continuous vector pathogen surveillance is essential for preventing outbreaks of mosquito-borne diseases. Several mosquito species acting as vectors of Japanese encephalitis virus (JEV), dengue virus, Zika virus, malaria parasites and other pathogens are primary mosquito species in Shanghai, China. However, few surveys of human pathogenic arboviruses in mosquitoes in Shanghai have been reported in the last ten years. Therefore, in this study, we evaluated mosquito activity in Shanghai, China during 2016 and tested for the presence of alphaviruses, flaviviruses, orthobunyaviruses and several parasitic pathogens. Results: Five pooled samples were JEV-positive [4/255 pools of Culex tritaeniorhynchus and 1/256 pools of Cx. pipiens (s.l.)] based on analysis of the NS5 gene. Alphaviruses, orthobunyaviruses, Plasmodium and filariasis were not found in this study. Phylogenetic and molecular analyses revealed that the JEV strains belonged to genotype I. Moreover, newly detected Shanghai JEV strains were genetically close to previously isolated Shandong strains responsible for transmission during the 2013 Japanese encephalitis (JE) outbreak in Shandong Province, China but were more distantly related to other Shanghai strains detected in the early 2000s. The E proteins of the newly detected Shanghai JEV strains differed from that in the live attenuated vaccine SA14-14-2-derived strain at six amino residues: E130 (Ile→Val), E222 (Ala→Ser), E327 (Ser→Thr), E366 (Arg→Ser/Pro), E393 (Asn→Ser) and E433 (Val→Ile). However, no differences were observed in key amino acid sites related to antigenicity. Minimum JEV infection rates were 1.01 and 0.65 per 1000 Cx. tritaeniorhynchus and Cx. pipiens (s.l.), respectively. Conclusions: Five new Shanghai JEV genotype I strains, detected after a ten-year hiatus in local mosquito surveillance, were genetically close to strains involved in the 2013 Shandong JE outbreak. Because JEV is still circulating, vaccination in children should be extensively and continuously promoted. Moreover, JEV mosquito surveillance programmes should document the genotype variation, intensity and distribution of circulating viruses for use in the development and implementation of disease prevention and control strategies. [ABSTRACT FROM AUTHOR]

Background: In sharks, chickens, rats, frogs, medaka and zebrafish there is haplotypic variation in MHC class I and closely linked genes involved in antigen processing, peptide translocation and peptide loading. At least in chicken, such MHCIa haplotypes of MHCIa, TAP2 and Tapasin are shown to influence the repertoire of pathogen epitopes being presented to CD8+ T-cells with subsequent effect on cell-mediated immune responses. Results: Examining MHCI haplotype variation in Atlantic salmon using transcriptome and genome resources we found little evidence for polymorphism in antigen processing genes closely linked to the classical MHCIa genes. Looking at other genes involved in MHCI assembly and antigen processing we found retention of functional gene duplicates originating from the second vertebrate genome duplication event providing cyprinids, salmonids, and neoteleosts with the potential of several different peptide-loading complexes. One of these gene duplications has also been retained in the tetrapod lineage with orthologs in frogs, birds and opossum. Conclusion: We postulate that the unique salmonid whole genome duplication (SGD) is responsible for eliminating haplotypic content in the paralog MHCIa regions possibly due to frequent recombination and reorganization events at early stages after the SGD. In return, multiple rounds of whole genome duplications has provided Atlantic salmon, other teleosts and even lower vertebrates with alternative peptide loading complexes. How this affects antigen presentation remains to be established. [ABSTRACT FROM AUTHOR]

Background: The species relationships within the genus Linum have already been studied several times by means of different molecular and phylogenetic approaches. Nevertheless, a number of ambiguities in phylogeny of Linum still remain unresolved. In particular, the species relationships within the sections Stellerolinum and Dasylinum need further clarification. Also, the question of independence of the species of the section Adenolinum still remains unanswered. Moreover, the relationships of L. narbonense and other species of the section Linum require further clarification. Additionally, the origin of tetraploid species of the section Linum (2n = 30) including the cultivated species L. usitatissimum has not been explored. The present study examines the phylogeny of blue-flowered species of Linum by comparisons of 5S rRNA gene sequences as well as ITS1 and ITS2 sequences of 35S rRNA genes. Results: High-throughput sequencing has been used for analysis of multicopy rRNA gene families. In addition to the molecular phylogenetic analysis, the number and chromosomal localization of 5S and 35S rDNA sites has been determined by FISH. Our findings confirm that L. stelleroides forms a basal branch from the clade of blue-flowered flaxes which is independent of the branch formed by species of the sect. Dasylinum. The current molecular phylogenetic approaches, the cytogenetic analysis as well as different genomic DNA fingerprinting methods applied previously did not discriminate certain species within the sect. Adenolinum. The allotetraploid cultivated species L. usitatissimum and its wild ancestor L. angustifolium (2n = 30) could originate either as the result of hybridization of two diploid species (2n = 16) related to the modern L. gandiflorum and L. decumbens, or hybridization of a diploid species (2n = 16) and a diploid ancestor of modern L. narbonense (2n = 14). Conclusions: High-throughput sequencing of multicopy rRNA gene families allowed us to make several adjustments to the phylogeny of blue-flowered flax species and also reveal intra- and interspecific divergence of the rRNA gene sequence. [ABSTRACT FROM AUTHOR]

Background: Babesia spp. are hemoparasites which infect the red blood cells of a large variety of mammals. In bats, the only known species of the genus is Babesia vesperuginis. However, except a few old reports, the host range and geographical distribution of this bat parasite have been poorly studied. This study aimed to investigate the presence of piroplasms in tissues of bats collected in four different countries from eastern and central Europe: Austria, Czech Republic, Hungary and Romania. Methods: A total of 461 bat carcasses (24 species) were collected between 2001 and 2016 from caves, mines and buildings. PCR was performed using specific primers targeting a portion of the 18S rDNA nuclear gene and cytochrome c oxidase subunit 1 mitochondrial gene, followed by sequencing. Results: The results of this study show for the first time the presence of B. vesperuginis in bats in central and eastern Europe. The phylogenetic analysis of the 18S rDNA nuclear gene revealed no variability between the sequences and the phylogenetic analysis of the cox1 mitochondrial gene proved that B. vesperuginis could be divided into two subclades. Conclusion: Our study showed a broad geographical distribution of B. vesperuginis in European bats, reporting its presence in five new host species (M. cf. alcathoe, M. bechsteinii, M. myotis, Pi. nathusii and V. murinus) and three new countries. [ABSTRACT FROM AUTHOR]

Background: The noble crayfish (Astacus astacus) displays a complex historical and contemporary genetic status in Europe. The species divergence has been shaped by geological events (i.e. Pleistocene glaciations) and humanly induced impacts (i.e. translocations, pollution, etc.) on its populations due to species commercial value and its niche degradation. Until now, limited genetic information has been procured for the Balkan area and especially for the southernmost distribution of this species (i.e. Greece). It is well known that the rich habitat diversity of the Balkan Peninsula offers suitable conditions for genetically diversified populations. Thus, the present manuscript revisits the phylogenetic relationships of the noble crayfish in Europe and identifies the genetic make-up and the biogeographical patterns of the species in its southern range limit. Results: Mitochondrial markers (i.e. COI and 16S) were used in order to elucidate the genetic structure and diversity of the noble crayfish in Europe. Two of the six European haplotypic lineages, were found exclusively in Greece. These two lineages exhibited greater haplotypic richness when compared with the rest four (of "Central European" origin) while they showed high genetic diversity. Divergence time analysis identified that the majority of this divergence was captured through Pleistocene, suggesting a southern glacial refugium (Greece, southern Balkans). Furthermore, six microsatellite markers were used in order to define the factors affecting the genetic structure and demographic history of the species in Greece. The population structure analysis revealed six to nine genetic clusters and eight putative genetic barriers. Evidence of bottleneck effects in the last ~5000 years (due to climatic and geological events and human activities) is also afforded. Findings from several other research fields (e.g. life sciences, geology or even archaeology) have been utilized to perceive the genetic make-up of the noble crayfish. Conclusions: The southernmost part of Balkans has played a major role as a glacial refugium for A. astacus. Such refugia have served as centres of expansion to northern regions. Recent history of the noble crayfish in southern Balkans reveals the influence of environmental (climate, geology and/or topology) and anthropogenic factors. [ABSTRACT FROM AUTHOR]

Background: Dispersal ability, population genetic structure and species divergence in marine nematodes are still poorly understood, especially in remote areas such as the Southern Ocean. We investigated genetic differentiation of species and populations of the free-living endobenthic nematode genera Sabatieria and Desmodora using nuclear 18S rDNA, internal transcribed spacer (ITS) rDNA, and mitochondrial cytochrome oxidase I (COI) gene sequences. Specimens were collected at continental shelf depths (200-500 m) near the Antarctic Peninsula, Scotia Arc and eastern side of the Weddell Sea. The two nematode genera co-occurred at all sampled locations, but with different vertical distribution in the sediment. A combination of phylogenetic (GMYC, Bayesian Inference, Maximum Likelihood) and population genetic (AMOVA) analyses were used for species delimitation and assessment of gene flow between sampling locations. Results: Sequence analyses resulted in the delimitation of four divergent species lineages in Sabatieria, two of which could not be discriminated morphologically and most likely constitute cryptic species. Two species were recognised in Desmodora, one of which showed large intraspecific morphological variation. Both genera comprised species that were restricted to one side of the Weddell Sea and species that were widely spread across it. Population genetic structuring was highly significant and more pronounced in the deeper sediment-dwelling Sabatieria species, which are generally less prone to resuspension and passive dispersal in the water column than surface Desmodora species. Conclusions: Our results indicate that gene flow is restricted at large geographic distance in the Southern Ocean, which casts doubt on the efficiency of the Weddell gyre and Antarctic Circumpolar Current in facilitating circum-Antarctic nematode species distributions. We also show that genetic structuring and cryptic speciation can be very different in nematode species isolated from the same geographic area, but with different habitat preferences (surface versus deeper sediment layers). [ABSTRACT FROM AUTHOR]

Background: Glass sponges (Class Hexactinellida) are important components of deep-sea ecosystems and are of interest from geological and materials science perspectives. The reconstruction of their phylogeny with molecular data has only recently begun and shows a better agreement with morphology-based systematics than is typical for other sponge groups, likely because of a greater number of informative morphological characters. However, inconsistencies remain that have far-reaching implications for hypotheses about the evolution of their major skeletal construction types (body plans). Furthermore, less than half of all described extant genera have been sampled for molecular systematics, and several taxa important for understanding skeletal evolution are still missing. Increased taxon sampling for molecular phylogenetics of this group is therefore urgently needed. However, due to their remote habitat and often poorly preserved museum material, sequencing all 126 currently recognized extant genera will be difficult to achieve. Utilizing morphological data to incorporate unsequenced taxa into an integrative systematics framework therefore holds great promise, but it is unclear which methodological approach best suits this task. Results: Here, we increase the taxon sampling of four previously established molecular markers (18S, 28S, and 16S ribosomal DNA, as well as cytochrome oxidase subunit I) by 12 genera, for the first time including representatives of the order Aulocalycoida and the type genus of Dactylocalycidae, taxa that are key to understanding hexactinellid body plan evolution. Phylogenetic analyses suggest that Aulocalycoida is diphyletic and provide further support for the paraphyly of order Hexactinosida; hence these orders are abolished from the Linnean classification. We further assembled morphological character matrices to integrate so far unsequenced genera into phylogenetic analyses in maximum parsimony (MP), maximum likelihood (ML), Bayesian, and morphology-based binning frameworks. We find that of these four approaches, total-evidence analysis using MP gave the most plausible results concerning congruence with existing phylogenetic and taxonomic hypotheses, whereas the other methods, especially ML and binning, performed more poorly. We use our total-evidence phylogeny of all extant glass sponge genera for ancestral state reconstruction of morphological characters in MP and ML frameworks, gaining new insights into the evolution of major hexactinellid body plans and other characters such as different spicule types. Conclusions: Our study demonstrates how a comprehensive, albeit in some parts provisional, phylogeny of a larger taxon can be achieved with an integrative approach utilizing molecular and morphological data, and how this can be used as a basis for understanding phenotypic evolution. The datasets and associated trees presented here are intended as a resource and starting point for future work on glass sponge evolution. [ABSTRACT FROM AUTHOR]

Background: Heterophyidiasis is now a major public health threat in many tropical countries. Species in the trematode family Heterophyidae infecting humans include Centrocestus formosanus, Haplorchis pumilio, H. taichui, H. yokogawai, Procerovum varium and Stellantchasmus falcatus. For molecular phylogenetic and systematic studies on trematodes, we need more prospective markers for taxonomic identification and classification. This study provides near-complete ribosomal transcription units (rTU) from Haplorchis pumilio and H. taichui and demonstrates the use of 28S rDNA sequences for identification and phylogenetic analysis. Results: The near-complete ribosomal transcription units (rTU), consisting of 18S, ITS1, 5.8S, ITS2 and 28S rRNA genes and spacers, from H. pumilio and H. taichui from human hosts in Vietnam, were determined and annotated. Sequence analysis revealed tandem repetitive elements in ITS1 in H. pumilio and in ITS2 in H. taichui. A phylogenetic tree inferred from 28S rDNA sequences of 40 trematode strains/species, including 14 Vietnamese heterophyid individuals, clearly confirmed the status of each of the Vietnamese species: Centrocestus formosanus, Haplorchis pumilio, H. taichui, H. yokogawai, Procerovum varium and Stellantchasmus falcatus. However, the family Heterophyidae was clearly not monophyletic, with some genera apparently allied with other families within the superfamily Opisthorchioidea (i.e. Cryptogonimidae and Opisthorchiidae). These families and their constituent genera require substantial re-evaluation using a combination of morphological and molecular data. Our new molecular data will assist in such studies. Conclusions: The 28S rDNA sequences are conserved among individuals within a species but varied between genera. Based on analysis of 40 28S rDNA sequences representing 19 species in the superfamily Opisthorchioidea and an outgroup taxon (Alaria alata, family Diplostomidae), six common human pathogenic heterophyids were identified and clearly resolved. The phylogenetic tree inferred from these sequences again confirmed anomalies in molecular placement of some members of the family Heterophyidae and demonstrates the need for reappraisal of the entire superfamily Opisthorchioidea. The new sequences provided here supplement those already available in public databases and add to the array of molecular tools that can be used for the diagnosis of heterophyid species in human and animal infections. [ABSTRACT FROM AUTHOR]

Background: Phage genome analysis is a rapidly growing field. Recurrent obstacles include software access and usability, as well as genome sequences that vary in sequence orientation and/or start position. Here we describe modifications to the phage comparative genomics software program, Phamerator, provide public access to the code, and include instructions for creating custom Phamerator databases. We further report genomic analysis techniques to determine phage packaging strategies and identification of the physical ends of phage genomes. Results: The original Phamerator code can be successfully modified and custom databases can be generated using the instructions we provide. Results of genome map comparisons within a custom database reveal obstacles in performing the comparisons if a published genome has an incorrect complementarity or an incorrect location of the first base of the genome, which are common issues in GenBank-downloaded sequence files. To address these issues, we review phage packaging strategies and provide results that demonstrate identification of the genome start location and orientation using raw sequencing data and software programs such as PAUSE and Consed to establish the location of the physical ends of the genome. These results include determination of exact direct terminal repeats (DTRs) or cohesive ends, or whether phages may use a headful packaging strategy. Phylogenetic analysis using ClustalO and phamily circles in Phamerator demonstrate that the large terminase gene can be used to identify the phage packaging strategy and thereby aide in identifying the physical ends of the genome. Conclusions: Using available online code, the Phamerator program can be customized and utilized to generate databases with individually selected genomes. These databases can then provide fruitful information in the comparative analysis of phages. Researchers can identify packaging strategies and physical ends of phage genomes using raw data from high-throughput sequencing in conjunction with phylogenetic analyses of large terminase proteins and the use of custom Phamerator databases. We promote publication of phage genomes in an orientation consistent with the physical structure of the phage chromosome and provide guidance for determining this structure. [ABSTRACT FROM AUTHOR]

Background: Nuclear factor Y (NF-Y) transcription factor is composed of three distinct subunits: NF-YA, NF-YB and NF-YC. Many members of NF-Y family have been reported to be key regulators in plant development, phytohormone signaling and drought tolerance. However, the function of the NF-Y family is less known in grape (Vitis vinifera L.). Results: A total of 34 grape NF-Y genes that distributed unevenly on grape (V. vinifera) chromosomes were identified in this study. Phylogenetic analysis was performed to predict functional similarities between Arabidopsis thaliana and grape NF-Y genes. Comparison of the structures of grape NF-Y genes (VvNF-Ys) revealed their functional conservation and alteration. Furthermore, we investigated the expression profiles of VvNF-Ys in response to various stresses, phytohormone treatments, and in leaves and grape berries with various sugar contents at different developmental stages. The relationship between VvNF-Y transcript levels and sugar content was examined to select candidates for exogenous sugar treatments. Quantitative real-time PCR (qPCR) indicated that many VvNF-Ys responded to different sugar stimuli with variations in transcript abundance. qPCR and publicly available microarray data suggest that VvNF-Ys exhibit distinct expression patterns in different grape organs and developmental stages, and a number of VvNF-Ys may participate in responses to multiple abiotic and biotic stresses, phytohormone treatments and sugar accumulation or metabolism. Conclusions: In this study, we characterized 34 VvNF-Ys based on their distributions on chromosomes, gene structures, phylogenetic relationship with Arabidopsis NF-Y genes, and their expression patterns. The potential roles of VvNF-Ys in sugar accumulation or metabolism were also investigated. Altogether, the data provide significant insights on VvNF-Ys, and lay foundations for further functional studies of NF-Y genes in grape. [ABSTRACT FROM AUTHOR]

Background: Sperm size and quality are key factors for fertilization success. There is increasing empirical evidence demonstrating that sperm form and function are influenced by selective pressures. Theoretical models predict that sperm competition could favour the evolution of longer sperm. In hermaphrodites, self-fertilizing species are expected to have shorter sperm than cross-fertilizing species, which use sperm stored from several mating partners for the fertilization of their eggs and thus are exposed to intense sperm competition. We tested this hypothesis by comparing original data on sperm length in 57 species of simultaneously hermaphroditic stylommatophoran gastropods from Europe and South America with respect to the species' breeding system. We used 28S rRNA nuclear and COI mitochondrial sequence data to construct a molecular phylogeny. Phylogenetic generalized linear models were applied to examine the potential influence of morphological and life-history characters. Results: The best-fit model revealed that the breeding system and age at sexual maturity influence sperm length in gastropods. In general, species with predominant cross-fertilization had longer sperm than species with predominant self-fertilization or a mixed breeding system. Across species with shells (snails), sperm length also increased with shell size. Conclusions: Our study provides evidence that sperm length in stylommatophoran gastropods is influenced by the risk of sperm competition, as well as by age at sexual maturity and shell size. This finding extends present knowledge of sperm evolution to a group of so far poorly studied simultaneous hermaphrodites. [ABSTRACT FROM AUTHOR]

Background: The type-2C protein phosphatases (PP2Cs), negatively regulating ABA responses and MAPK cascade pathways, play important roles in stress signal transduction in plants. Brachypodium distachyon is a new model plant for exploring the functional genomics of temperate grasses, cereals and biofuel crops. To date, genome-wide identification and analysis of the PP2C gene family in B. distachyon have not been investigated. Results: In this study, 86 PP2C genes in B. distachyon were identified. Domain-based analyses of PP2C proteins showed that they all contained the phosphatase domains featured as 11 conserved signature motifs. Although not all phosphatase domains of BdPP2C members included all 11 motifs, tertiary structure analysis showed that four residues contributing to magnesium/manganese ions (Mg2+/Mn2+) coordination were conserved, except for two noncanonical members. The analysis of their chromosomal localizations showed that most of the BdPP2C genes were located within the low CpG density region. Phylogenetic tree and synteny blocks analyses among B. distachyon, Arabidopsis thaliana and Oryza sativa revealed that all PP2C members from the three species can be phylogenetically categorized into 13 subgroups (A-M) and BdPP2Cs were evolutionarily more closely related to OsPP2Cs than to AtPP2Cs. Segmental duplications contributed particularly to the expansion of the BdPP2C gene family and all duplicated BdPP2Cs evolved mainly from purifying selection. Real-time quantitative reverse transcription PCR (qRT-PCR) analysis showed that BdPP2Cs were broadly expressed in disparate tissues. We also found that almost all members displayed up-regulation in response to abiotic stresses such as cold, heat, PEG and NaCl treatments, but down-regulation to biotic stresses such as Ph14, Guy11 and F0968 infection. Conclusions: In the present study, a comprehensive analysis of genome-wide identification and characterization of protein domains, phylogenetic relationship, gene and protein structure, chromosome location and expression pattern of the PP2C gene family was carried out for the first time in a new model monocot, i.e., B. distachyon. Our results provide a reference for genome-wide identification of the PP2C gene family of other species and also provide a foundation for future functional research on PP2C genes in B. distachyon. [ABSTRACT FROM AUTHOR]

Background: The accurate annotation of genes in newly sequenced genomes remains a challenge. Although sophisticated comparative pipelines are available, computationally derived gene models are often less than perfect. This is particularly true when multiple similar paralogs are present. The issue is aggravated further when genomes are assembled only at a preliminary draft level to contigs or short scaffolds. However, these genomes deliver valuable information for studying gene families. High accuracy models of protein coding genes are needed in particular for phylogenetics and for the analysis of gene family histories. Results: We present a pipeline, ExonMatchSolver, that is designed to help the user to produce and curate high quality models of the protein-coding part of genes. The tool in particular tackles the problem of identifying those coding exon groups that belong to the same paralogous genes in a fragmented genome assembly. This paralog-tocontig assignment problem is shown to be NP-complete. It is phrased and solved as an Integer Linear Programming problem. Conclusions: The ExonMatchSolver-pipeline can be employed to build highly accurate models of protein coding genes even when spanning several genomic fragments. This sets the stage for a better understanding of the evolutionary history within particular gene families which possess a large number of paralogs and in which frequent gene duplication events occurred. [ABSTRACT FROM AUTHOR]

Background: Maternal care (egg-nymph guarding behavior) has been recorded in some genera of Acanthosomatidae. However, the origin of the maternal care in the family has remained unclear due to the lack of phylogenetic hypotheses. Another reproductive mode is found in non-caring species whose females smear their eggs before leaving them. They possess pairs of complex organs on the abdominal venter called Pendergrast's organ (PO) and spread the secretion of this organ onto each egg with their hind legs, which is supposed to provide a protective function against enemies. Some authors claim that the absence of PO may be associated with the presence of maternal care. No study, however, has tested this hypothesis of a correlated evolution between the two traits. Results: We reconstructed the molecular phylogeny of the subfamily Acanthosomatinae using five genetic markers sequenced from 44 species and one subspecies with and without maternal care. Eight additional species from the other two acanthosomatid subfamilies were included as outgroups. Our results indicated that maternal care has evolved independently at least three times within Acanthosomatinae and once in the outgroup species. Statistical tests for correlated evolution showed that the presence of maternal care is significantly correlated with the secondary loss or reduction of PO. Ancestral state reconstruction for the node of Acanthosoma denticaudum (a non-caring species in which egg smearing with developed POs occurs) and A. firmatum (a caring species with reduced POs) suggested egg smearing was still present in their most recent common ancestor and that maternal care in A. firmatum has evolved relatively recently. Conclusions: We showed that maternal care is an apomorphic trait that has arisen multiple times from the presence of PO within the subfamily Acanthosomatinae. The acquisition of maternal care is correlated with the reduction or loss of PO, which suggests an evolutionary trade-off between the two traits resulting from physiological costs. This prediction also implies that presence of maternal care can be highly expected for those groups lacking behavioral data, which invariably also lack the organ. No secondary loss of maternal care was detected in the present tree. We suggest that the loss of maternal care may be suppressed due to the vulnerability of the PO-free condition, which thus maintains maternal care. [ABSTRACT FROM AUTHOR]

Background: Daldinia eschscholtzii is a wood-inhabiting fungus that causes wood decay under certain conditions. It has a broad host range and produces a large repertoire of potentially bioactive compounds. However, there is no extensive genome analysis on this fungal species. Results: Two fungal isolates (UM 1400 and UM 1020) from human specimens were identified as Daldinia eschscholtzii by morphological features and ITS-based phylogenetic analysis. Both genomes were similar in size with 10,822 predicted genes in UM 1400 (35.8 Mb) and 11,120 predicted genes in UM 1020 (35.5 Mb). A total of 751 gene families were shared among both UM isolates, including gene families associated with fungus-host interactions. In the CAZyme comparative analysis, both genomes were found to contain arrays of CAZyme related to plant cell wall degradation. Genes encoding secreted peptidases were found in the genomes, which encode for the peptidases involved in the degradation of structural proteins in plant cell wall. In addition, arrays of secondary metabolite backbone genes were identified in both genomes, indicating of their potential to produce bioactive secondary metabolites. Both genomes also contained an abundance of gene encoding signaling components, with three proposed MAPK cascades involved in cell wall integrity, osmoregulation, and mating/filamentation. Besides genomic evidence for degrading capability, both isolates also harbored an array of genes encoding stress response proteins that are potentially significant for adaptation to living in the hostile environments. Conclusions: Our genomic studies provide further information for the biological understanding of the D. eschscholtzii and suggest that these wood-decaying fungi are also equipped for adaptation to adverse environments in the human host. [ABSTRACT FROM AUTHOR]

Background: -glutamyltranspeptidase (GGT) is a bi-substrate enzyme conserved in all three domains of life. It catalyzes the cleavage and transfer of -glutamyl moiety of glutathione to either water (hydrolysis) or substrates like peptides (transpeptidation). GGTs exhibit great variability in their enzyme kinetics although the mechanism of catalysis is conserved. Recently, GGT has been shown to be a virulence factor in microbes like Helicobacter pylori and Bacillus anthracis. In mammalian cells also, GGT inhibition prior to chemotherapy has been shown to sensitize tumors to the therapy. Therefore, lately both bacterial and eukaryotic GGTs have emerged as potential drug targets, but the efforts directed towards finding suitable inhibitors have not yielded any significant results yet. We propose that delineating the residues responsible for the functional diversity associated with these proteins could help in design of species/clade specific inhibitors. Results: In the present study, we have carried out phylogenetic analysis on a set of 47 GGT-like proteins to address the functional diversity. These proteins segregate into various subfamilies, forming separate clades on the tree. Sequence conservation and motif prediction studies show that even though most of the highly conserved residues have been characterized biochemically in previous studies, a significant number of novel putative sites and motifs are discovered that vary in a clade specific manner. Many of the putative sites predicted during the functional divergence type I and type II analysis, lie close to the known catalytic residues and line the walls of the substrate binding cavity, reinforcing their role in modulating the substrate specificity, catalytic rates and stability of this protein. Conclusion: The study offers interesting insights into the evolution of GGT-like proteins in pathogenic vs. nonpathogenic bacteria, archaea and eukaryotes. Our analysis delineates residues that are highly specific to each GGT subfamily. We propose that these sites not only explain the differences in stability and catalytic variability of various GGTs but can also aid in design of specific inhibitors against particular GGTs. Thus, apart from the commonly used in-silico inhibitor screening approaches, evolutionary analysis identifying the functional divergence hotspots in GGT proteins could augment the structure based drug design approaches. [ABSTRACT FROM AUTHOR]

Background: Species of Bryopsidales form ecologically important components of seaweed communities worldwide. These siphonous macroalgae are composed of a single giant tubular cell containing millions of nuclei and chloroplasts, and harbor diverse bacterial communities. Little is known about the diversity of chloroplast genomes (cpDNAs) in this group, and about the possible consequences of intracellular bacteria on genome composition of the host. We present the complete cpDNAs of Bryopsis plumosa and Tydemania expeditiones, as well as a re-annotated cpDNA of B. hypnoides, which was shown to contain a higher number of genes than originally published. Chloroplast genomic data were also used to evaluate phylogenetic hypotheses in the Chlorophyta, such as monophyly of the Ulvophyceae (the class in which the order Bryopsidales is currently classified). Results: Both DNAs are circular and lack a large inverted repeat. The cpDNA of B. plumosa is 106,859 bp long and contains 115 unique genes. A 13 kb region was identified with several freestanding open reading frames (ORFs) of putative bacterial origin, including a large ORF (>8 kb) closely related to bacterial rhs-family genes. The cpDNA of T. expeditiones is 105,200 bp long and contains 125 unique genes. As in B. plumosa, several regions were identified with ORFs of possible bacterial origin, including genes involved in mobile functions (transposases, integrases, phage/plasmid DNA primases), and ORFs showing close similarity with bacterial DNA methyltransferases. The cpDNA of B. hypnoides differs from that of B. plumosa mainly in the presence of long intergenic spacers, and a large tRNA region. Chloroplast phylogenomic analyses were largely inconclusive with respect to monophyly of the Ulvophyceae, and the relationship of the Bryopsidales within the Chlorophyta. Conclusions: The cpDNAs of B. plumosa and T. expeditiones are amongst the smallest and most gene dense chloroplast genomes in the core Chlorophyta. The presence of bacterial genes, including genes typically found in mobile elements, suggest that these have been acquired through horizontal gene transfer, which may have been facilitated by the occurrence of obligate intracellular bacteria in these siphonous algae. [ABSTRACT FROM AUTHOR]

Background: A central problem of computational metagenomics is determining the correct placement into an existing phylogenetic tree of individual reads (nucleotide sequences of varying lengths, ranging from hundreds to thousands of bases) obtained using next-generation sequencing of DNA samples from a mixture of known and unknown species. Correct placement allows us to easily identify or classify the sequences in the sample as to taxonomic position or function. Results: Here we propose a novel method (PhyClass), based on the Minimum Evolution (ME) phylogenetic inference criterion, for determining the appropriate phylogenetic position of each read. Without using heuristics, the new approach efficiently finds the optimal placement of the unknown read in a reference phylogenetic tree given a sequence alignment for the taxa in the tree. In short, the total resulting branch length for the tree is computed for every possible placement of the unknown read and the placement that gives the smallest value for this total is the best (optimal) choice. By taking advantage of computational efficiencies and mathematical formulations, we are able to find the true optimal ME placement for each read in the phylogenetic tree. Using computer simulations, we assessed the accuracy of the new approach for different read lengths over a variety of data sets and phylogenetic trees. We found the accuracy of the new method to be good and comparable to existing Maximum Likelihood (ML) approaches. Conclusions: In particular, we found that the consensus assignments based on ME and ML approaches are more correct than either method individually. This is true even when the statistical support for read assignments was low, which is inevitable given that individual reads are often short and come from only one gene. [ABSTRACT FROM AUTHOR]

Background: The correct taxonomic assignment of bacterial genomes is a primary and challenging task. With the availability of whole genome sequences, the gene content based approaches appear promising in inferring the bacterial taxonomy. The complete genome sequencing of a bacterial genome often reveals a substantial number of unique genes present only in that genome which can be used for its taxonomic classification. Results: In this study, we have proposed a comprehensive method which uses the taxon-specific genes for the correct taxonomic assignment of existing and new bacterial genomes. The taxon-specific genes identified at each taxonomic rank have been successfully used for the taxonomic classification of 2,342 genomes present in the NCBI genomes, 36 newly sequenced genomes, and 17 genomes for which the complete taxonomy is not yet known. This approach has been implemented for the development of a tool 'Microtaxi' which can be used for the taxonomic assignment of complete bacterial genomes. Conclusion: The taxon-specific gene based approach provides an alternate valuable methodology to carry out the taxonomic classification of newly sequenced or existing bacterial genomes. [ABSTRACT FROM AUTHOR]

Background: The human (Homo sapiens) chemokine-like protein macrophage migration inhibitory factor (HsMIF) is a pivotal mediator of inflammatory, infectious and immune diseases including septic shock, colitis, malaria, rheumatoid arthritis, and atherosclerosis, as well as tumorigenesis. HsMIF has been found to exhibit several sequential and three-dimensional sequence motifs that in addition to its receptor binding sites include catalytic sites for oxidoreductase and tautomerase activity, which provide this 12.5 kDa protein with a remarkable functional complexity. A human MIF paralog, D-dopachrome tautomerase (HsDDT), has been identified, but its physiological relevance is incompletely understood. MIF/ DDT-like proteins have been described in animals, protists and bacteria. Although based on sequence data banks the presence of MIF/DDT-like proteins has also been recognized in the model plant species Arabidopsis thaliana, details on these plant proteins have not been reported. Results: To broaden the understanding of the biological role of these proteins across kingdoms we performed a comprehensive in silico analysis of plant MIF/DDT-like (MDL) genes/proteins. We found that the A. thaliana genome harbors three MDL genes, of which two are chiefly constitutively expressed in aerial plant organs, while the third gene shows stress-inducible transcript accumulation. The product of the latter gene likely localizes to peroxisomes. Structure prediction suggests that all three Arabidopsis proteins resemble the secondary and tertiary structure of human MIF. MIF-like proteins are found in all species across the plant kingdom, with an increasing family complexity towards evolutionarily advanced plant taxa. Plant MDL proteins are predicted to lack oxidoreductase activity, but possibly share tautomerase activity with human MIF/DDT. Conclusions: Peroxisome localization seems to be a specific feature of a subset of MIF/DDT orthologs found in dicotyledonous plant species, which together with its stress-inducible gene expression might point to convergent evolution in higher plants and vertebrates towards neofunctionalization of MIF/MDL proteins in stress response pathways including innate immunity. [ABSTRACT FROM AUTHOR]

Background: MHC class I (MHCI) molecules are the key presenters of peptides generated through the intracellular pathway to CD8-positive T-cells. In fish, MHCI genes were first identified in the early 1990's, but we still know little about their functional relevance. The expansion and presumed sub-functionalization of cod MHCI and access to many published fish genome sequences provide us with the incentive to undertake a comprehensive study of deduced teleost fish MHCI molecules. Results: We expand the known MHCI lineages in teleosts to five with identification of a new lineage defined as P. The two lineages U and Z, which both include presumed peptide binding classical/typical molecules besides more derived molecules, are present in all teleosts analyzed. The U lineage displays two modes of evolution, most pronouncedly observed in classical-type alpha 1 domains; cod and stickleback have expanded on one of at least eight ancient alpha 1 domain lineages as opposed to many other teleosts that preserved a number of these ancient lineages. The Z lineage comes in a typical format present in all analyzed ray-finned fish species as well as lungfish. The typical Z format displays an unprecedented conservation of almost all 37 residues predicted to make up the peptide binding groove. However, also co-existing atypical Z sub-lineage molecules, which lost the presumed peptide binding motif, are found in some fish like carps and cavefish. The remaining three lineages, L, S and P, are not predicted to bind peptides and are lost in some species. Conclusions: Much like tetrapods, teleosts have polymorphic classical peptide binding MHCI molecules, a number of classical-similar non-classical MHCI molecules, and some members of more diverged MHCI lineages. Different from tetrapods, however, is that in some teleosts the classical MHCI polymorphism incorporates multiple ancient MHCI domain lineages. Also different from tetrapods is that teleosts have typical Z molecules, in which the residues that presumably form the peptide binding groove have been almost completely conserved for over 400 million years. The reasons for the uniquely teleost evolution modes of peptide binding MHCI molecules remain an enigma. [ABSTRACT FROM AUTHOR]

Background The Central and Southern European mountain ranges represent important biodiversity hotspots and show high levels of endemism. In the land snail genus Orcula Held, 1837 nine species are distributed in the Alps and a few taxa inhabit the Carpathians, the Dinarids and the Western Black Sea region. In order to elucidate the general patterns of temporal and geographic diversification, mitochondrial and nuclear markers were analyzed in all 13 Orcula species. We particularly aimed to clarify whether the Alpine taxa represent a monophyletic group and if the local species diversity is rather the result of isolation in geographically separated Pleistocene glacial refuges or earlier Tertiary and Quaternary palaeogeographic events. In order to test if patterns of molecular genetic and morphological differentiation were congruent and/or if hybridization had occurred, shell morphometric investigations were performed on the Orcula species endemic to the Alps. Results The phylogenetic trees resulting from the analyses of both the mitochondrial (COI, 12S and 16S) and the nuclear (H4/H3) data sets revealed three main groups, which correspond to the three subgenera Orcula, Illyriobanatica and Hausdorfia. The reconstruction of the historic geographic ranges suggested that the genus originated in the Dinarides during the Middle Miocene and first colonized the Alps during the Late Miocene, giving rise to the most diverse subgenus Orcula. Within the latter subgenus (including all Alpine endemics) almost all species were differentiated by both molecular genetic markers and by shell morphometrics, except O. gularis and O. pseudodolium. Conclusions The present study confirms the importance of the Alps as biodiversity hotspot and origin center of land snail diversity. The species diversity in the subgenus Orcula was likely promoted by Miocene to Pliocene palaeogeographic events and the insular distribution of preferred limestone areas. In some cases, speciation events could be linked to the divergence of populations in glacial refuges during the Pleistocene. Sporadic contact between geographically separated and reproductively not yet isolated populations led to intermixture of haplogroups within species and even hybridization and mitochondrial capture between species. [ABSTRACT FROM AUTHOR]

Background: The animal mitochondrial genome is generally considered to be under selection for both compactness and gene order conservation. As more mitochondrial genomes are sequenced, mitochondrial duplications and gene rearrangements have been frequently identified among diverse animal groups. Although several mechanisms of gene rearrangement have been proposed thus far, more observational evidence from major taxa is needed to validate specific mechanisms. In the current study, the complete mitochondrial DNA of sixteen bird species from the family Ardeidae was sequenced and the evolution of mitochondrial gene rearrangements was investigated. The mitochondrial genomes were then used to review the phylogenies of these ardeid birds. Results: The complete mitochondrial genome sequences of the sixteen ardeid birds exhibited four distinct mitochondrial gene orders in which two of them, named as “duplicate tRNAGlu–CR” and “duplicate tRNAThr–tRNAPro and CR”, were newly discovered. These gene rearrangements arose from an evolutionary process consistent with the tandem duplication - random loss model (TDRL). Additionally, duplications in these gene orders were near identical in nucleotide sequences within each individual, suggesting that they evolved in concert. Phylogenetic analyses of the sixteen ardeid species supported the idea that Ardea ibis, Ardea modesta and Ardea intermedia should be classified as genus Ardea, and Ixobrychus flavicollis as genus Ixobrychus, and indicated that within the subfamily Ardeinae, Nycticorax nycticorax is closely related to genus Egretta and that Ardeola bacchus and Butorides striatus are closely related to the genus Ardea. Conclusions: The duplicate tRNAThr–CR gene order is found in most ardeid lineages, suggesting this gene order is the ancestral pattern within these birds and persisted in most lineages via concerted evolution. In two independent lineages, when the concerted evolution stopped in some subsections due to the accumulation of numerous substitutions and deletions, the duplicate tRNAThr–CR gene order was transformed into three other gene orders. The phylogenetic trees produced from concatenated rRNA and protein coding genes have high support values in most nodes, indicating that the mitochondrial genome sequences are promising markers for resolving the phylogenetic issues of ardeid birds when more taxa are added. [ABSTRACT FROM AUTHOR]

Background NASP is an essential protein in mammals that functions in histone transport pathways and maintenance of a soluble reservoir of histones H3/H4. NASP has been studied exclusively in Opisthokonta lineages where some functional diversity has been reported. In humans, growing evidence implicates NASP miss-regulation in the development of a variety of cancers. Although a comprehensive phylogenetic analysis is lacking, NASP-family proteins that possess four TPR motifs are thought to be widely distributed across eukaryotes. Results We characterize the molecular evolution of NASP by systematically identifying putative NASP orthologs across diverse eukaryotic lineages ranging from excavata to those of the crown group. We detect extensive silent divergence at the nucleotide level suggesting the presence of strong purifying selection acting at the protein level. We also observe a selection bias for high frequencies of acidic residues which we hypothesize is a consequence of their critical function(s), further indicating the role of functional constraints operating on NASP evolution. Our data indicate that TPR1 and TPR4 constitute the most rapidly evolving functional units of NASP and may account for the functional diversity observed among well characterized family members. We also show that NASP paralogs in ray-finned fish have different genomic environments with clear differences in their GC content and have undergone significant changes at the protein level suggesting functional diversification. Conclusion We draw four main conclusions from this study. First, wide distribution of NASP throughout eukaryotes suggests that it was likely present in the last eukaryotic common ancestor (LECA) possibly as an important innovation in the transport of H3/H4. Second, strong purifying selection operating at the protein level has influenced the nucleotide composition of NASP genes. Further, we show that selection has acted to maintain a high frequency of functionally relevant acidic amino acids in the region that interrupts TPR2. Third, functional diversity reported among several well characterized NASP family members can be explained in terms of quickly evolving TPR1 and TPR4 motifs. Fourth, NASP fish specific paralogs have significantly diverged at the protein level with NASP2 acquiring a NNR domain. [ABSTRACT FROM AUTHOR]

Background Nucleotide and amino acid substitution tendencies are characteristic of each species, organelle, and protein family. Hence, various empirical amino acid substitution rate matrices have needed to be estimated for phylogenetic analysis: JTT, WAG, and LG for nuclear proteins, mtREV for mitochondrial proteins, cpREV10 and cpREV64 for chloroplast-encoded proteins, and FLU for influenza proteins. On the other hand, in a mechanistic codon substitution model, in which each codon substitution rate is proportional to the product of a codon mutation rate and the ratio of fixation depending on the type of amino acid replacement, mutation rates and the strength of selective constraint on amino acids can be tailored to each protein family with additional 11 parameters. As a result, in the evolutionary analysis of codon sequences it outperforms codon substitution models equivalent to empirical amino acid substitution matrices. Is it superior even for amino acid sequences, among which synonymous substitutions cannot be identified? Results Nucleotide mutations are assumed to occur independently of codon positions but multiple nucleotide changes in infinitesimal time are allowed. Selective constraints on the respective types of amino acid replacements are tailored to each gene with a linear function of a given estimate of selective constraints, which were estimated by maximizing the likelihood of an empirical amino acid or codon substitution frequency matrix, each of JTT, WAG, LG, and KHG. It is shown that the mechanistic codon substitution model with the assumption of equal codon usage yields better values of Akaike and Bayesian information criteria for all three phylogenetic trees of mitochondrial, chloroplast, and influenza-A hemagglutinin proteins than the empirical amino acid substitution models with mtREV, cpREV64, and FLU, which were designed specifically for those protein families, respectively. The variation of selective constraint across sites fits the datasets significantly better than variable codon mutation rates, confirming that substitution rate variations across sites detected by amino acid substitution models are caused primarily by the variation of selective constraint against amino acid substitutions rather than the variation of codon mutation rate. Conclusions The mechanistic codon substitution model is superior to amino acid substitution models even in the evolutionary analysis of protein sequences. [ABSTRACT FROM AUTHOR]

Background Lyme disease is caused by spirochete bacteria from the Borrelia burgdorferi sensu lato (B. burgdorferi s.l.) species complex. To reconstruct the evolution of B. burgdorferi s.l. and identify the genomic basis of its human virulence, we compared the genomes of 23 B. burgdorferi s.l. isolates from Europe and the United States, including B. burgdorferi sensu stricto (B. burgdorferi s.s., 14 isolates), B. afzelii (2), B. garinii (2), B. "bavariensis" (1), B. spielmanii (1), B. valaisiana (1), B. bissettii (1), and B. "finlandensis" (1). Results Robust B. burgdorferi s.s. and B. burgdorferi s.l. phylogenies were obtained using genomewide single-nucleotide polymorphisms, despite recombination. Phylogeny-based pan-genome analysis showed that the rate of gene acquisition was higher between species than within species, suggesting adaptive speciation. Strong positive natural selection drives the sequence evolution of lipoproteins, including chromosomally-encoded genes 0102 and 0404, cp26- encoded ospC and b08, and lp54-encoded dbpA, a07, a22, a33, a53, a65. Computer simulations predicted rapid adaptive radiation of genomic groups as population size increases. Conclusions Intra- and inter-specific pan-genome sizes of B. burgdorferi s.l. expand linearly with phylogenetic diversity. Yet gene-acquisition rates in B. burgdorferi s.l. are among the lowest in bacterial pathogens, resulting in high genome stability and few lineage-specific genes. Genome adaptation of B. burgdorferi s.l. is driven predominantly by copy-number and sequence variations of lipoprotein genes. New genomic groups are likely to emerge if the current trend of B. burgdorferi s.l. population expansion continues. [ABSTRACT FROM AUTHOR]

Background: The species Pectobacterium carotovorum includes a diverse subspecies of bacteria that cause disease on a wide variety of plants. In Morocco, approximately 95% of the P. carotovorum isolates from potato plants with tuber soft rot are P. carotovorum subsp. carotovorum. However, identification of this pathogen is not always related to visual disease symptoms. This is especially true when different pathogen cause similar diseases on potato, citing as an example, P. carotovorum, P. atrosepticum and P. wasabiae. Numerous conventional methods were used to characterize Pectobacterium spp., including biochemical assays, specific PCR-based tests, and construction of phylogenetic trees by using gene sequences. In this study, an alternative method is presented using a gene linked to pathogenicity, in order to allow accuracy at subspecies level. The pmrA gene (response regulator) has been used for identification and analysis of the relationships among twenty nine Pectobacterium carotovorum subsp. carotovorum and other Pectobacterium subspecies. Results: Phylogenetic analyses of pmrA sequences compared to ERIC-PCR and 16S rDNA sequencing, demonstrated that there is considerable genetic diversity in P. carotovorum subsp. carotovorum strains, which can be divided into two distinct groups within the same clade. Conclusions: pmrA sequence analysis is likely to be a reliable tool to identify the subspecies Pectobacterium carotovorum subsp. carotovorum and estimate their genetic diversity. [ABSTRACT FROM AUTHOR]

Background: Phylogenetic comparison of bacteriophages requires whole genome approaches such as dotplot analysis, genome pairwise maps, and gene content analysis. Currently mycobacteriophages, a highly studied phage group, are categorized into related clusters based on the comparative analysis of whole genome sequences. With the recent explosion of phage isolation, a simple method for phage cluster prediction would facilitate analysis of crude or complex samples without whole genome isolation and sequencing. The hypothesis of this study was that mycobacteriophage-cluster prediction is possible using comparison of a single, ubiquitous, semi-conserved gene. Tape Measure Protein (TMP) was selected to test the hypothesis because it is typically the longest gene in mycobacteriophage genomes and because regions within the TMP gene are conserved. Results: A single gene, TMP, identified the known Mycobacteriophage clusters and subclusters using a Gepard dotplot comparison or a phylogenetic tree constructed from global alignment and maximum likelihood comparisons. Gepard analysis of 247 mycobacteriophage TMP sequences appropriately recovered 98.8% of the subcluster assignments that were made by whole-genome comparison. Subcluster-specific primers within TMP allow for PCR determination of the mycobacteriophage subcluster from DNA samples. Using the single-gene comparison approach for siphovirus coliphages, phage groupings by TMP comparison reflected relationships observed in a whole genome dotplot comparison and confirm the potential utility of this approach to another widely studied group of phages. Conclusions: TMP sequence comparison and PCR results support the hypothesis that a single gene can be used for distinguishing phage cluster and subcluster assignments. TMP single-gene analysis can quickly and accurately aid in mycobacteriophage classification. [ABSTRACT FROM AUTHOR]

Background: RNAalifold, a popular computational method for RNA consensus structure prediction, incorporates covarying mutations into a thermodynamic model to fold the aligned RNA sequences. When quantifying covariance, it evaluates conserved signals of two aligned columns with base-pairing rules. This scoring scheme performs better than some other approaches, such as mutual information. However it ignores the phylogenetic history of the aligned sequences, which is an important criterion to evaluate the level of sequence covariance. Results: In this article, in order to improve the accuracy of consensus structure folding, we propose a novel approach named PhyloRNAalifold. It incorporates the number of covarying mutations on the phylogenetic tree of the aligned sequences into the covariance scoring of RNAalifold. The benchmarking results show that the new scoring scheme of PhyloRNAalifold can improve the consensus structure detection of RNAalifold. Conclusion: Incorporating additional phylogenetic information of aligned sequences into the covariance scoring of RNAalifold can improve its performance of consensus structures folding. This improvement is correlated with alignment characteristics, such as pair-wise identity and the number of sequences in the alignment. [ABSTRACT FROM AUTHOR]

Background: The extant squamates (>9400 known species of lizards and snakes) are one of the most diverse and conspicuous radiations of terrestrial vertebrates, but no studies have attempted to reconstruct a phylogeny for the group with large-scale taxon sampling. Such an estimate is invaluable for comparative evolutionary studies, and to address their classification. Here, we present the first large-scale phylogenetic estimate for Squamata. Results: The estimated phylogeny contains 4161 species, representing all currently recognized families and subfamilies. The analysis is based on up to 12896 base pairs of sequence data per species (average = 2497 bp) from 12 genes, including seven nuclear loci (BDNF, c-mos, NT3, PDC, R35, RAG-1, and RAG-2), and five mitochondrial genes (12S, 16S, cytochrome b, ND2, and ND4). The tree provides important confirmation for recent estimates of higher-level squamate phylogeny based on molecular data (but with more limited taxon sampling), estimates that are very different from previous morphology-based hypotheses. The tree also includes many relationships that differ from previous molecular estimates and many that differ from traditional taxonomy. Conclusions: We present a new large-scale phylogeny of squamate reptiles that should be a valuable resource for future comparative studies. We also present a revised classification of squamates at the family and subfamily level to bring the taxonomy more in line with the new phylogenetic hypothesis. This classification includes new, resurrected, and modified subfamilies within gymnophthalmid and scincid lizards, and boid, colubrid, and lamprophiid snakes [ABSTRACT FROM AUTHOR]

Background: Three species of the genus Cercopithifilaria have been morphologically and molecularly characterized in dog populations in southern Europe: Cercopithifilaria grassii (Noè, 1907), Cercopithifilaria sp. sensu Otranto et al., 2011 (reported as Cercopithifilaria sp. I), and Cercopithifilaria sp. II sensu Otranto et al., 2012. The adults of Cercopithifilaria sp. I have remained unknown until the present study. Methods: The material originated from a dog from Sardinia (Italy) diagnosed with dermal microfilariae of Cercopithifilaria sp. I. The holotype and three paratypes of Cercopithifilaria bainae Almeida & Vicente, 1984, described from dogs in Brazil, were studied as comparative material. A cox1 (∼689 bp) and 12S (∼330 bp) gene fragments were amplified and phylogenetic analysis carried out. Results: The highest numbers of adult nematodes (82%) were collected in the sediment of the subcutaneous tissues of the trunk (n = 37) and forelimbs (n = 36). The morphology of the adult nematodes and microfilariae collected from the dog in Sardinia corresponded to those of C. bainae. All cox1 and 12S gene sequences showed a high homology (99-100%) with sequences from microfilariae of Cercopithifilaria sp. I. Conclusions: The morphological and molecular identity of the microfilariae of C. bainae overlap those described previously as Cercopithifilaria sp. sensu Otranto et al., 2011 (=Cercopithifilaria sp. I). Therefore, the present study reports the occurrence of C. bainae in Europe, for the first time after its description and the single record in Brazil. C. bainae appears to be highly diffused in dog populations in southern Europe. The phylogenetic analyses based on cox1 and 12S do not reveal the three species of Cercopithifilaria parasitizing dogs as a monophyletic group, which suggests that they have derived independently by host switching. [ABSTRACT FROM AUTHOR]

Background: Strongylocentrotid sea urchins have a long tradition as model organisms for studying many fundamental processes in biology including fertilization, embryology, development and genome regulation but the phylogenetic relationships of the group remain largely unresolved. Although the differing isolating mechanisms of vicariance and rapidly evolving gamete recognition proteins have been proposed, a stable and robust phylogeny is unavailable. Results: We used a phylogenomic approach with mitochondrial and nuclear genes taking advantage of the whole-genome sequencing of nine species in the group to establish a stable (i.e. concordance in tree topology among multiple lies of evidence) and robust (i.e. high nodal support) phylogenetic hypothesis for the family Strongylocentrotidae. We generated eight draft mitochondrial genome assemblies and obtained 13 complete mitochondrial genes for each species. Consistent with previous studies, mitochondrial sequences failed to provide a reliable phylogeny. In contrast, we obtained a very well-supported phylogeny from 2301 nuclear genes without evidence of positive Darwinian selection both from the majority of most-likely gene trees and the concatenated fourfold degenerate sites: ((P. depressus, (M. nudus, M. franciscanus), (H. pulcherrimus, (S. purpuratus, (S. fragilis, (S. pallidus, (S. droebachiensis, S. intermedius)). This phylogeny was consistent with a single invasion of deep-water environments followed by a holarctic expansion by Strongylocentrotus. Divergence times for each species estimated with reference to the divergence times between the two major clades of the group suggest a correspondence in the timing with the opening of the Bering Strait and the invasion of the holarctic regions. Conclusions: Nuclear genome data contains phylogenetic signal informative for understanding the evolutionary history of this group. However, mitochondrial genome data does not. Vicariance can explain major patterns observed in the phylogeny. Other isolating mechanisms are appropriate to explore in this system to help explain divergence patterns not well supported by vicariance, such as the effects of rapidly evolving gamete recognition proteins on isolating populations. Our findings of a stable and robust phylogeny, with the increase in mitochondrial and nuclear comparative genomic data, provide a system in which we can enhance our understanding of molecular evolution and adaptation in this group of sea urchins. [ABSTRACT FROM AUTHOR]

Background: After the occurrence of 2009 pandemic H1N1, close attention has been paid to the H1N1 subtype swine influenza viruses (H1N1 SIV) by scientific communities in many countries. A large-scale sequence analysis of the NCBI Influenza Virus Resource Database on H1N1 SIVs submitted primarily by scientists in China during 1992 to 2011 was performed. The aims of this study were to elucidate the genetic and evolutionary characteristics of H1N1 SIVs, to identify and unify the lineages and genetic characteristics of the H1N1 SIVs isolated in mainland China. Results: Most of the strains were isolated during the period of 2008 to 2010 from Guangdong and Shandong provinces, China. Based on the phylogenetic and genotypic analyses, all of the H1N1 SIV strains can be classified into 8 lineages and 10 genotypes. All strains were of the characteristics of low pathogenic influenza viruses. The viruses of different lineage are characterized with different amino acid residues at the receptor-binding sites. Viruses containing PB2 genes of the classical swine, early seasonal human and recent seasonal human lineage might be more infectious to human. Some genotypes were directly related with human influenza viruses, which include strains that harbored genes derived from human influenza viruses. Conclusions: Phylogenetic diversity and complexity existed in H1N1 SIVs isolated in mainland China. These H1N1 SIV strains were closely related to other subtype influenza viruses, especially to human influenza viruses. Moreover, it was shown that, novel lineages and genotypes of H1N1 SIVs emerged recently in mainland China. These findings provided new and essential information for further understanding of the genetic and evolutionary characteristics and monitoring the H1N1 SIVs in mainland China. [ABSTRACT FROM AUTHOR]

Background: Anguillicolidae Yamaguti, 1935 is a family of parasitic nematode infecting fresh-water eels of the genus Anguilla, comprising five species in the genera Anguillicola and Anguillicoloides. Anguillicoloides crassus is of particular importance, as it has recently spread from its endemic range in the Eastern Pacific to Europe and North America, where it poses a significant threat to new, naïve hosts such as the economic important eel species Anguilla anguilla and Anguilla rostrata. The Anguillicolidae are therefore all potentially invasive taxa, but the relationships of the described species remain unclear. Anguillicolidae is part of Spirurina, a diverse clade made up of only animal parasites, but placement of the family within Spirurina is based on limited data.Results: We generated an extensive DNA sequence dataset from three loci (the 5' one-third of the nuclear small subunit ribosomal RNA, the D2-D3 region of the nuclear large subunit ribosomal RNA and the 5' half of the mitochondrial cytochrome c oxidase I gene) for the five species of Anguillicolidae and used this to investigate specific and generic boundaries within the family, and the relationship of Anguillicolidae to other spirurine nematodes. Neither nuclear nor mitochondrial sequences supported monophyly of Anguillicoloides. Genetic diversity within the African species Anguillicoloides papernai was suggestive of cryptic taxa, as was the finding of distinct lineages of Anguillicoloides novaezelandiae in New Zealand and Tasmania. Phylogenetic analysis of the Spirurina grouped the Anguillicolidae together with members of the Gnathostomatidae and Seuratidae.Conclusions: The Anguillicolidae is part of a complex radiation of parasitic nematodes of vertebrates with wide host diversity (chondrichthyes, teleosts, squamates and mammals), most closely related to other marine vertebrate parasites that also have complex life cycles. Molecular analyses do not support the recent division of Anguillicolidae into two genera. The described species may hide cryptic taxa, identified here by DNA taxonomy, and this DNA barcoding approach may assist in tracking species invasions. The propensity for host switching, and thus the potential for invasive behaviour, is found in A. crassus, A. novaezelandiae and A. papernai, and thus may be common to the group. [ABSTRACT FROM AUTHOR]

Background: The evolution of mycoplasmas from a common ancestor with Firmicutes has been characterized not only by genome down-sizing but also by horizontal gene transfer between mycoplasma species sharing a common host. The mechanisms of these gene transfers remain unclear because our knowledge of the mycoplasma mobile genetic elements is limited. In particular, only a few plasmids have been described within the Mycoplasma genus. Results: We have shown that several species of ruminant mycoplasmas carry plasmids that are members of a large family of elements and replicate via a rolling-circle mechanism. All plasmids were isolated from species that either belonged or were closely related to the Mycoplasma mycoides cluster; none was from the Mycoplasma bovis-Mycoplasma agalactiae group. Twenty one plasmids were completely sequenced, named and compared with each other and with the five mycoplasma plasmids previously reported. All plasmids share similar size and genetic organization, and present a mosaic structure. A peculiar case is that of the plasmid pMyBK1 from M. yeatsii; it is larger in size and is predicted to be mobilizable. Its origin of replication and replication protein were identified. In addition, pMyBK1 derivatives were shown to replicate in various species of the M. mycoides cluster, and therefore hold considerable promise for developing gene vectors. The phylogenetic analysis of these plasmids confirms the uniqueness of pMyBK1 and indicates that the other mycoplasma plasmids cluster together, apart from the related replicons found in phytoplasmas and in species of the clade Firmicutes. Conclusions: Our results unraveled a totally new picture of mycoplasma plasmids. Although they probably play a limited role in the gene exchanges that participate in mycoplasma evolution, they are abundant in some species. Evidence for the occurrence of frequent genetic recombination strongly suggests they are transmitted between species sharing a common host or niche. [ABSTRACT FROM AUTHOR]

Roseophage RDJLϕ1 is a siphovirus isolated from South China Sea on Roseobacter denitrificans OCh114. Its virion encapsulates 62.7 kb genome that encodes 87 gene products. RDJLϕ1 shares similar genome organization and gene content with the marine bacteriophage FJL001 and Pseudomonas phages YuA and M6, which are different from those of typical l- or Mu-like phages. Four hallmark genes (ORFs 81 to 84) of RDJLϕ1 were highly homologous to RcGTA-like genes 12 to 15. The largest gene (ORF 84) was predicted to encode a tail fibre protein that could be involved in host recognition. Extended phylogenetic and comparative genomic analyses based on 77 RcGTA-like element-containing bacterial genomes revealed that RcGTA-like genes 12 to 15 together appear to be a conserved modular element that could also be found in some phage or prophage genomes. Our study suggests that RcGTA-like genes-containing phages and prophages and complete RcGTAs possibly descended from a same prophage ancestor that had diverged and then evolved vertically. The complete genome of RDJLϕ1 provides evidence into the hypothesis that extant RcGTA may be a prophage remnant. [ABSTRACT FROM AUTHOR]

Background: We present a novel method to encode ambiguously aligned regions in fixed multiple sequence alignments by 'Pairwise Identity and Cost Scores Ordination' (PICS-Ord). The method works via ordination of sequence identity or cost scores matrices by means of Principal Coordinates Analysis (PCoA). After identification of ambiguous regions, the method computes pairwise distances as sequence identities or cost scores, ordinates the resulting distance matrix by means of PCoA, and encodes the principal coordinates as ordered integers. Three biological and 100 simulated datasets were used to assess the performance of the new method. Results: Including ambiguous regions coded by means of PICS-Ord increased topological accuracy, resolution, and bootstrap support in real biological and simulated datasets compared to the alternative of excluding such regions from the analysis a priori. In terms of accuracy, PICS-Ord performs equal to or better than previously available methods of ambiguous region coding (e.g., INAASE), with the advantage of a practically unlimited alignment size and increased analytical speed and the possibility of PICS-Ord scores to be analyzed together with DNA data in a partitioned maximum likelihood model. Conclusions: Advantages of PICS-Ord over step matrix-based ambiguous region coding with INAASE include a practically unlimited number of OTUs and seamless integration of PICS-Ord codes into phylogenetic datasets, as well as the increased speed of phylogenetic analysis. Contrary to word- and frequency-based methods, PICS-Ord maintains the advantage of pairwise sequence alignment to derive distances, and the method is flexible with respect to the calculation of distance scores. In addition to distance and maximum parsimony, PICS-Ord codes can be analyzed in a Bayesian or maximum likelihood framework. RAxML (version 7.2.6 or higher that was developed for this study) allows up to 32-state ordered or unordered characters. A GTR, MK, or ORDERED model can be applied to analyse the PICS-Ord codes partition, with GTR performing slightly better than MK and ORDERED. Availability: An implementation of the PICS-Ord algorithm is available from http://scit.us/projects/ngila/wiki/PICSOrd. It requires both the statistical software, R http://www.r-project.org and the alignment software Ngila http://scit. us/projects/ngila. [ABSTRACT FROM AUTHOR]

Background: Wnt5 genes belong to the large Wnt family, encoding proteins implicated into several tumorigenic and developmental processes. Phylogenetic analyses showed that Wnt5 gene has been duplicated at the divergence time of gnathostomata from agnatha. Interestingly, experimental data for some species indicated that only one of the two Wnt5 paralogs participates in the development of the endocrine pancreas. The purpose of this paper is to reexamine the phylogenetic history of the Wnt5 developmental regulators and investigate the functional shift between paralogs through comparative genomics. Results: In this study, the phylogeny of Wnt5 genes was investigated in species belonging to protostomia and deuterostomia. Furthermore, an in silico regulatory region analysis of Wnt5 paralogs was conducted, limited to those species with insulin producing cells and pancreas, covering the evolutionary distance from agnatha to gnathostomata. Our results confirmed the Wnt5 gene duplication and additionally revealed that this duplication event included also the upstream region. Moreover, within this latter region, a conserved module was detected to which a complex of transcription factors, known to be implicated in embryonic pancreas formation, bind. Conclusions: Results and observations presented in this study, allow us to conclude that during evolution, the Wnt5 gene has been duplicated in early vertebrates, and that some paralogs conserved a module within their regulatory region, functionally related to embryonic development of pancreas. Interestingly, our results allowed advancing a possible explanation on why the Wnt5 orthologs do not share the same function during pancreas development. As a final remark, we suggest that an in silico comparative analysis of regulatory regions, especially when associated to published experimental data, represents a powerful approach for explaining shift of roles among paralogs. Reviewers: This article was reviewed by Sarath Janga (nominated by Sarah Teichmann), Ran Kafri (nominated by Yitzhak Pilpel), and Andrey Mironov (nominated by Mikhail Gelfand). [ABSTRACT FROM AUTHOR]