Mitochondrial genomes from major lizard families suggest their phylogenetic relationships and ancient radiations

In placental mammals and birds, molecular data generally support a view that they diverged into their ordinal groups in good response to mid-Cretaceous continental fragmentations. However, such divergence patterns have rarely been studied for reptiles for which phylogenetic relationships among their major groups have not yet been established molecularly. Here, I determined complete or nearly complete mitochondrial DNA sequences from seven lizard families and reconstructed phylogenetic relationships between major lizard families. When snakes were included, maximum likelihood analysis did not support a morphological view of the snakes-varanoids affinity, although several other competing hypotheses on the position of snakes still cannot be discriminated presumably due to extremely long branches of the snake lineages. I also conducted clock-free Bayesian analyses to show that divergence times between major lizard families were centered in Triassic-Jurassic times. Thus, lizards include much deeper divergences than the mammals and birds and they appear to have already radiated into various families prior to the mid-Cretaceous major continental fragmentation.     

In silico analysis of SSRs in mitochondrial genomes of plants

Simple sequence repeats (SSRs) or microsatellites constitute a countable portion of genomes. However, the significance of SSRs in organelle genomes has not been completely understood. The availability of organelle genome sequences allows us to understand the organization of SSRs in their genic and intergenic regions. In the current study we surveyed the patterns of SSRs in mitochondrial genomes of different taxa of plants. A total of 16 mitochondrial genomes, from algae to angiosperms, have been considered to analyze the pattern of simple sequence repeats present in them. Based on study, the mononucleotide repeats of A/T were found to be more prevalent in mitochondrial genomes over other repeat types. The dinucleotides repeats, TA/AT, were the second most numerous, whereas tri-, tetra-, and pentanucleotide repeats were in less number and present in intronic or intergenic portions only. Mononucleotide repeats prevailed in protein-coding exonic portions of all organisms. These results indicates that microsatellite pattern in mitochondrial genomes is different from nuclear genomes and also focuses on organization and diversity at SSR locuses in mitochondrial genomes. This is the novel report of microsatellite polymorphism in plant mitochondrion on whole genome level.     

Immunological relationships among the major seed proteins of cereals

The structural relationships among the major seed proteins of cereals was evaluated by Western blot analyses using antibodies raised against the wheat gliadin, rice glutelin acidic and basic subunits, and rice prolamine polypeptide. Consitent with the conservation of the primary sequences of these proteins, antibodies to the acidic and basic glutelin subunits cross-reacted with homologous polypeptides from oat as well as pea. The rice glutelin antibodies did not react with the major seed proteins from barley, rye, maize and sorghum. Antibodies raised against the acidic glutelin subunit reacted with the wheat glutenins but antibodies to the basic glutelin subunit did not. A comparison of the published primary sequences of a high molecular weight glutenin and rice glutelin showed little similarity except for a conserved peptide with the motif arg-gln-leu-gln-cys. The possible significance of this conserved element shared by these widely different proteins is discussed. Similar studies with the wheat gliadin antibody showed immunologically related components in plants of the subfamily Festucoideae except for rice. Antibodies raised against the rice prolamine recognized only the rice prolamine, indicating that this polypeptide was structurally distinct from other cereal prolamines. Overall, these results support and help clarify the evolutionary relationship of the cereals.     

In silico whole-genome EST analysis reveals 2322 novel microsatellites for the silver-lipped pearl oyster, Pinctada maxima

Molecular stock improvement techniques such as marker assisted selection have great potential in accelerating selective breeding programmes for animal production industries. However, the discovery and application of trait/marker associations usually requires a large number of genome-wide polymorphic loci. Here, we present 2322 unique microsatellites for the silver-lipped pearl oyster, Pinctada maxima, a species of aquaculture importance throughout the Indo-Australian Archipelago for production of the highly valued South Sea pearl. More than 1.2 million Roche 454 expressed sequence tag (EST) reads were screened for microsatellite repeat motifs. A total of 12,604 sequences contained either a di, tri, tetra, penta or hexa microsatellite repeat motif (n ≥ 6), with 6435 of these sequences having sufficient flanking regions for primer development. All identified microsatellites with designed primers were condensed into 2322 unique clusters (i.e., unique loci) of which 360 were shown to be polymorphic based on multiple sequence reads with different repeat motifs. Genotyping of five microsatellite loci demonstrated that in silico evaluation of polymorphism levels was a very useful method for identification of polymorphic loci, with the variation uncovered being a lower bound. Gene Ontology annotations of sequences containing microsatellites suggest that most are derived from a diverse array of unique genes. This EST derived microsatellite database will be a valuable resource for future studies in genetic map construction, diversity analysis, quantitative trait loci analysis, association mapping and marker assisted selection, not only for P. maxima, but also closely related species within the genus Pinctada.     

In silico analysis of complete bacterial genomes: PCR, AFLP-PCR and endonuclease restriction

We have developed a website, www.in-silico.com, which runs a software program that performs three basic tasks in completely sequenced bacterial genomes by in silico analysis: PCR amplification, amplified fragment length polymorphism (AFLP-PCR) and endonuclease restriction. For PCR, after selection of the genome and introduction of primers, fragment size, DNA sequence and corresponding open reading frame (ORF) identity of the resulting PCR product is computed. Plasmids of sequenced species may be included in the analysis. Theoretical AFLP-PCR analyzes similar parameters, and includes a suggestion tool providing a list of commercial restriction enzyme pairs yielding up to 50 amplicons in the selected genome. Endonuclease restriction analysis of complete genomes and plasmids calculates the number of restriction sites for endonucleases in a given genome. If the number of fragments is 50 or fewer, pulsed field gel electrophoresis image and restriction maps are illustrated. Other tools that have been included in this site are ORF search by name and DNA to protein translation as well as restriction digestion of user-defined DNA sequences. AVAILABILITY: This is a new molecular biology resource freely available over the Internet at http://www.in-silico.com     

The analysis of microsatellites and compound microsatellites in 56 complete genomes of Herpesvirales

Simple sequence repeats (SSRs), or microsatellites, are special DNA/RNA sequences with repeated unit of 1–6 bp. The genomes of Herpesvirales have many repeating structures, which is an excellent system to study the evolution and roles of microsatellites and compound microsatellites in viruses. Therefore, 56 genomes of Herpesvirales were selected and the occurrence, composition and complexity of different repeats were investigated in the genomes. A total of 63,939 microsatellites and 5825 compound microsatellites were extracted from 56 genomes. It found that GC content has a significant strong correlation with both the counts of microsatellites (CM) and the counts of compound microsatellites (CCM). However, genome size has a moderate correlation only with CM and almost no correlation with CCM. The compound microsatellites occurring in genic regions are obviously more than that in intergenic regions. In general, the number of compound microsatellite decreases with the increase of complexity (C) (the count of individual microsatellites being part of a compound microsatellite) and the complexity hardly exceeds C = 4. The vast majority of compound microsatellites exist in intergenic regions, when C ≥ 10. The distributions of SSRs tend to be organism-specific rather than host-specific in herpesvirus genomes. The diversity of microsatellites and compound microsatellites may be helpful for a better understanding of the viral genetic diversity, genotyping, and evolutionary biology in herpesviruses genomes.     

In silico mining in expressed sequences of Neurospora crassa for identification and abundance of microsatellites

In the present study, 3217 UniGene sequences of Neurospora crassa downloaded from the National Center for Biotechnology Information (NCBI) were mined for the identification of microsatellites or simple sequence repeats (SSRs). A total of 287 SSRs detected gives density of 1SSR/14.6 kb of 4187.86 kb sequences mined suggests that only 250 (7.8%) of sequences contained SSRs. Depending on the repeat units, the length of SSRs ranged from 14 to 17 bp for mono-, 14 to 48 bp for di-, 18 to 90 bp for tri-, 24 to 48 bp for tetra-, 30 for penta- and 42 to 48 bp for hexa-nucleotide repeats. Tri-nucleotide repeats were the most frequent repeat type (88.8%) followed by di-nucleotide repeats (5.9%). An attempt was also made with the help of bioinformatics approach to find out primer pairs for identified SSRs and primers were found only for 239 sequences. But, this part needs experimental validation. Annotation of SSRs containing sequences was also carried out.     

Occurrence and analysis of imperfect microsatellites in diverse potyvirus genomes

Simple sequence repeats (SSRs) or microsatellites are known to exhibit ubiquitous across all kingdoms of life including viruses. However, imperfections in simple sequence repeats have been analyzed in genomes of human, Escherichia coli and Human Immunodeficiency virus. The assessment of compound microsatellites in plant viral genomes is yet to be studied. Potyviruses severely affect crop plant growth and reduce economic yield in diverse cropping systems worldwide. Hence, we analyze the nature and distribution of compound microsatellites present in complete genome of 45 potyvirus species. The results indicate that compound microsatellites accounted for about 0% to 15.15% of all microsatellites and have low complexity as compared to that of prokaryotic genomes. Overall, 14% of compound microsatellites were of similar motifs and such motif duplications were observed for CA, TA and AG repeats. Among all 45 potyvirus genomes analyzed, SSR couple (AG)-x-(AC) was found to be the most abundant one. Hence it is apparent that in contrast to eukaryotes, majority of compound microsatellites in potyviruses were composed of variant motifs. We also highlight the relative frequency of different classes of compound microsatellites as well as their patterns of distribution and correlate with biology of potyviruses. Further characterization of such variation is important for elucidating the origin, mutational processes, and structure of these widely used, but incompletely understood sequences.     

Phylogenetic positioning of the Antarctic alga Prasiola crispa (Trebouxiophyceae) using organellar genomes and their structural analysis

Antarctica is one of the most difficult habitats for sustaining life on earth; organisms that live there have developed different strategies for survival. Among these organisms is the green alga Prasiola crispa, belonging to the class Trebouxiophyceae. The literature on P. crispa taxonomy is scarce, and many gaps in the evolutionary relationship with its closest relatives remain. The goal of this study was to analyze the evolutionary relationships between P. crispa and other green algae using plastid and mitochondrial genomes. In addition, we analyzed the synteny conservation of these genomes of P. crispa with those of closely related species. Based on the plastid genome, P. crispa grouped with Prasiolopsis sp. SAG 84.81, another Trebouxiophyceaen species from the Prasiola clade. Based on the mitochondrial genome analysis, P. crispa grouped with other Trebouxiophyceaen species but had a basal position. The structure of the P. crispa chloroplast genome had low synteny with Prasiolopsis sp. SAG 84.81, despite some conserved gene blocks. The same was observed in the mitochondrial genome compared with Coccomyxa subellipsoidea C‐169. We were able to establish the phylogenetic position of P. crispa with other species of Trebouxiophyceae using its genomes. In addition, we described the plasticity of these genomes using a structural analysis. The plastid and mitochondrial genomes of P. crispa will be useful for further genetic studies, phylogenetic analysis and resource protection of P. crispa as well as for further phylogenetic analysis of Trebouxiophyceaen green algae.     

In silico identification and Bayesian phylogenetic analysis of multiple new mammalian kallikrein gene families

Kallikrein gene families have been identified previously in genomes of the human, the mouse, and the rat, and individual kallikrein-like genes have been found in many more species. This study presents the in silico identification of kallikrein gene families in the recently sequenced genomes of four additional mammalian species, the chimpanzee, the dog, the pig, and the opossum. Phylogenies were constructed with gene sequences from all seven mammalian families, using Bayesian analysis, which clarified the evolutionary relationships between these genes. Individual gene sequences, as well as concatenated constructs of multiple sequences, were used. Fifteen kallikrein genes were located in the chimpanzee (Pan troglodytes) genome, while only 14 were identified in the canine (Canis familiaris) genome as no orthologue to human KLK3 was found. Thirteen genes were identified from the pig (Sus scrofa) genome, which lacked homologues to KLK2 and KLK3, and 11 genes, orthologous to human KLK5 through KLK15, were found in the opossum (Monodelphis domestica) genome. No kallikrein genes were identified from the available genome sequences of the chicken (Gallus gallus) or African clawed frog (Xenopus tropicalis). Within the family of kallikreins several subfamilies were suggested by phylogenetic analysis. One consisted of KLK4, KLK5, and KLK14; another of KLK9, KLK11, and KLK15; a third of KLK10 and KLK12; a fourth of KLK6 and KLK13; and finally one of KLK8 and the classical kallikreins (KLK1, KLK2, and KLK3).     

In vitro and in silico cloning of Xenopus laevis SOD2 cDNA and its phylogenetic analysis

By using the methodology of both wet and dry biology (i.e., RT-PCR and cycle sequencing, and biocomputational technology, respectively) and the data obtained through the Genome Projects, we have cloned Xenopus laevis SOD2 (MnSOD) cDNA and determined its nucleotide sequence. These data and the deduced protein primary structure were compared with all the other SOD2 nucleotide and amino acid sequences from eukaryotes and prokaryotes, published in public databases. The analysis was performed by using both Clustal W, a well known and widely used program for sequence analysis, and AntiClustAl, a new algorithm recently created and implemented by our group. Our results demonstrate a very high conservation of the enzyme amino acid sequence during evolution, which proves a close structure-function relationship. This is to be expected for very ancient molecules endowed with critical biological functions, performed through a specific structural organization. The nucleotide sequence conservation is less pronounced: this too was foreseeable, due to neutral mutations and to the species-specific codon usage. The data obtained by using AntiClustAl are comparable with those produced with Clustal W, which validates this algorithm as an important new tool for biocomputational analysis. Finally, it is noteworthy that evolutionary trees, drawn by using all the available data on SOD2 nucleotide sequences and amino acid and either Clustal W or AntiClustAl, are comparable to those obtained through phylogenetic analysis based on fossil records.     

Housekeeping genes for phylogenetic analysis of eutherian relationships

The molecular relationship of placental mammals has attracted great interest in recent years. However, 2 crucial and conflicting hypotheses remain, one with respect to the position of the root of the eutherian tree and the other the relationship between the orders Rodentia, Lagomorpha (rabbits, hares), and Primates. Although most mitochondrial (mt) analyses have suggested that rodents have a basal position in the eutherian tree, some nuclear data in combination with mt-rRNA genes have placed the root on the so-called African clade or on a branch that includes this clade and the Xenarthra (e.g., anteater and armadillo). In order to generate a new and independent set of molecular data for phylogenetic analysis, we have established cDNA sequences from different tissues of various mammalian species. With this in mind, we have identified and sequenced 8 housekeeping genes with moderately fast rate of evolution from 22 placental mammals, representing 11 orders. In order to determine the root of the eutherian tree, the same genes were also sequenced for 3 marsupial species, which were used as outgroup. Inconsistent with the analyses of nuclear + mt-rRNA gene data, the current data set did not favor a basal position of the African clade or Xenarthra in the eutherian tree. Similarly, by joining rodents and lagomorphs on the same basal branch (Glires hypothesis), the data set is also inconsistent with the tree commonly favored in mtDNA analyses. The analyses of the currently established sequences have helped examination of problematic parts in the eutherian tree at the same time as they caution against suggestions that have claimed that basal eutherian relationships have been conclusively settled.     

In silico prediction of proteins related to xyloglucan fucosyltransferases in Solanaceae genomes

Two independent studies have shown that the cell wall of pollen tubes from tobacco and tomato species contained fucosylated xyloglucan (XyG). These findings are intriguing as many reports have shown that XyG of somatic cells of these species is not fucosylated but instead is arabinosylated. In order to produce fucosylated XyG, plants must express a functional galactoside α-2-fucosyltransferase. Here, using a bioinformatics approach, we show that several candidate genes coding for XyG fucosyltransferases are present in the genome of coffee and several Solanaceae species including tomato, tobacco, potato, eggplant and pepper. BLAST and protein alignments with the 2 well-characterized XyG fucosyltransferases from Arabidopsis thaliana and Pisum sativum revealed that at least 6 proteins from different Solanaceae species and from coffee displayed the 3 conserved motifs required for XyG fucosyltransferase activity.     

Mitochondrial genomes of five Hyphessobrycon tetras and their phylogenetic implications

To date, the taxonomic status and phylogenetic affinities within Hyphessobrycon, even among other genera in Characidae, remain unclear. Here, we determined five new mitochondrial genomes (mitogenomes) of Hyphessobrycon species (H. elachys, H. flammeus, H. pulchripinnis, H. roseus, and H. sweglesi). The mitogenomes were all classical circular structures, with lengths ranging from 16,008 to 17,224 bp. The type of constitutive genes and direction of the coding strand that appeared in the mitogenomes were identical to those of other species in Characidae. The highest value of the Ka/Ks ratio within 13 protein‐coding genes (PCGs) was found in ND2 with 0.83, suggesting that they were subject to purifying selection in the Hyphessobrycon genus. Comparison of the control region sequences among seven Hyphessobrycon fish revealed that repeat units differ in length and copy number across different species, which led to sharp differences in mitogenome sizes. Phylogenetic trees based on the 13 PCGs did not support taxonomic relationships, as the Hyphessobrycon fish mixed with those from other genera. These data were combined to explore higher level relationships within Characidae and could aid in the understanding of the evolution of this group.     

The isochore patterns of mammalian genomes and their phylogenetic implications

The compositional distributions of high molecular weight DNA fragments from 20 species belonging to 9 out of the 17 eutherian orders were investigated by analytical CsCl density gradient centrifugation and by preparative fractionation in Cs₂SO₄/BAMD density gradients followed by analysis of the fractions in CsCl. These compositional distributions reflect those of the isochores making up the corresponding genomes. A “general distribution” was found in species belonging to eight mammalian orders. A “myomorph distribution” was found in Myomorpha, but not in the other rodent infraorders Sciuromorpha and Histricomorpha, which share the general distribution. Two other distributions were found in a megachiropteran (but not in microchiropteran, which, again, shares the general distribution) and in pangolin (a species from the only genus of the order Pholidota), respectively. The main difference between the general distribution and all other distributions is that the former contains sizable amounts (6–10%) of GC-rich isochores (detected as DNA fragments equal to, or higher than, 1.710 g/cm³ in modal buoyant density), which are scarce, or absent, in the other distributions. This difference is remarkable because gene concentrations in mammalian genomes are paralleled by GC levels, the highest gene concentrations being present in the GC-richest isochores. The compositional distributions of mammalian genomes reported here shed light on mammalian phylogeny. Indeed, all orders investigated, with the exception of Pholidota, seem to share a common ancestor. The compositional patterns of the megachiropteran and of Myomorpha may be derived from the general pattern or have independent origins.     

Mitochondrial genomes in the Hymenoptera and their utility as phylogenetic markers

Abstract.  In this study, we assessed the ability of mitochondrial genome sequences to recover a test phylogeny of five hymenopteran taxa from which phylogenetic relationships are well accepted. Our analyses indicated that the test phylogeny was well recovered in all nucleotide Bayesian analyses when all the available holometabolan (i.e. outgroup) taxa were included, but only in Bayesian analyses excluding third codon positions when only the hymenopteran representatives and a single outgroup were included. This result suggests that taxon sampling of the outgroup might be as important as taxon sampling of the ingroup when recovering hymenopteran phylogenetic relationships using whole mitochondrial genomes. Parsimony analyses were more sensitive to both taxon sampling and the analytical model than Bayesian analyses, and analyses using the protein dataset did not recover the test phylogeny. In general, mitochondrial genomes did not resolve the position of the Hymenoptera within the Holometabola with confidence, suggesting that an increased taxon sampling, both within the Holometabola and among outgroups, is necessary.     

Genome-wide analysis of SINA family in plants and their phylogenetic relationships(dagger).

SINA genes in plants are part of a multigene family with 5 members in Arabidopsis thaliana, 10 members in Populus trichocarpa, 6 members in Oryza sativa, at least 6 members in Zea mays and at least 1 member in Physcomitrella patens. Six members in maize were confirmed by RT-PCR. All SINAs have one RING domain and one SINA domain. These two domains are highly conserved in plants. According to the motif organization and phylogenetic tree, SINA family members were divided into 2 groups. In addition, through semi-quantitative RT-PCR analysis of maize members and Digital Northern analysis of Arabidopsis and rice members, we found that the tissue expression patterns are more diverse in monocot than in Arabidopsis.     

GMATo: A novel tool for the identification and analysis of microsatellites in large genomes

Simple Sequence Repeats (SSR), also called microsatellite, is very useful for genetic marker development and genome application. The increasing whole sequences of more and more large genomes provide sources for SSR mining in silico. However currently existing SSR mining tools can’t process large genomes efficiently and generate no or poor statistics. Genome-wide Microsatellite Analyzing Tool (GMATo) is a novel tool for SSR mining and statistics at genome aspects. It is faster and more accurate than existed tools SSR Locator and MISA. If a DNA sequence was too long, it was chunked to short segments at several Mb followed by motifs generation and searching using Perl powerful pattern match function. Matched loci data from each chunk were then merged to produce final SSR loci information. Only one input file is required which contains raw fasta DNA sequences and output files in tabular format list all SSR loci information and statistical distribution at four classifications. GMATo was programmed in Java and Perl with both graphic and command line interface, either executable alone in platform independent manner with full parameters control. Software GMATo is a powerful tool for complete SSR characterization in genomes at any size.

Availability

The soft GMATo is freely available at http://sourceforge.net/projects/gmato/files/?source=navbar or on contact     

Reconstruction of certain phylogenetic networks from the genomes at their leaves

A network N is a rooted acyclic digraph. A base-set X for N is a subset of vertices including the root (or outgroup), all leaves, and all vertices of outdegree 1. A simple model of evolution is considered in which all characters are binary and in which back-mutations occur only at hybrid vertices. It is assumed that the genome is known for each member of the base-set X. If the network is known and is assumed to be “normal,” then it is proved that the genome of every vertex is uniquely determined and can be explicitly reconstructed. Under additional hypotheses involving time-consistency and separation of the hybrid vertices, the network itself can also be reconstructed from the genomes of all members of X. An explicit polynomial-time procedure is described for performing the reconstruction.