Comparative genomic analysis of eutherian Mas-related G protein-coupled receptor genes

The present study made attempts to update comprehensive eutherian Mas-related G protein-coupled receptor gene data sets, using public eutherian genomic sequence data sets and new genomics and molecular evolution tests. Among 254 potential coding sequences, the most comprehensive gene data set of eutherian Mas-related G protein-coupled receptor genes included 119 complete coding sequences that described eight major gene clusters. The present analysis integrated gene annotations, phylogenetic analysis and protein molecular evolution analysis and first explained differential gene expansion patterns of eutherian Mas-related G protein-coupled receptor genes. The updated classification and nomenclature of eutherian Mas-related G protein-coupled receptor genes were proposed as new framework of future experiments.     

Comparative genomic analysis of eutherian interferon-γ-inducible GTPases

The interferon-γ-inducible GTPases, IFGGs, are intracellular proteins involved in immune response against pathogens. A comprehensive comparative genomic review and analysis of eutherian IFGGs was carried out using public genomic sequences. The 64 eutherian IFGG genes were examined in detail and annotated. The eutherian IFGG promoter types were first catalogued followed by a phylogenetic analysis of eutherian IFGGs, which described five major IFGG clusters. The patterns of differential gene expansions and protein regions that may regulate IFGG catalytic features suggested a new classification of eutherian IFGGs. This mini-review has also provided new tests of reliability of public genomic sequences as well as tests of protein molecular evolution.     

GST-PRIME: a genome-wide primer design software for the generation of gene sequence tags

The availability of sequenced genomes has generated a need for experimental approaches that allow the simultaneous analysis of large, or even complete, sets of genes. To facilitate such analyses, we have developed GST-PRIME, a software package for retrieving and assembling gene sequences, even from complex genomes, using the NCBI public database, and then designing sets of primer pairs for use in gene amplification. Primers were designed by the program for the direct amplification of gene sequence tags (GSTs) from either genomic DNA or cDNA. Test runs of GST-PRIME on 2000 randomly selected Arabidopsis and Drosophila genes demonstrate that 93 and 88% of resulting GSTs, respectively, fulfilled imposed length criteria. GST-PRIME primer pairs were tested on a set of 1900 Arabidopsis genes coding for chloroplast-targeted proteins: 95% of the primer pairs used in PCRs with genomic DNA generated the correct amplicons. GST-PRIME can thus be reliably used for large-scale or specific amplification of intron-containing genes of multicellular eukaryotes.     

Molecular evolution of the HBII-52 snoRNA cluster

HBII-52 is a human brain-specific C/D box snoRNA that potentially regulates the editing and/or alternative splicing of the serotonin receptor. Forty-two nearly identical copies of the HBII-52 gene are located immediately downstream of the SNRPN protein-coding gene in an imprinted locus associated with Prader-Willi syndrome. Other eutherian mammals, with genomic assemblies covering the corresponding locus, also have multiple orthologous copies of HBII-52. The SNRPB gene, which is known to have given rise to SNRPN through gene duplication, expresses a C/D box snoRNA, SNORD119, from its fifth intron. Here we show that, despite the fact that they lie in different positions relative to the orthologous SNRPB/SNRPN coding sequences, there are significant sequence similarities between SNORD119 and HBII-52, including the antisense element and the stem-forming regions. By analysing these snoRNAs in marsupial and eutherian mammal genomes, we reconstruct the likely evolutionary history of the HBII-52 cluster and SNORD119 and suggest that they have evolved from a common ancestor.     

Genomic structure and sequence of the pufferfish (Fugu rubripes) growth hormone-encoding gene: a comparative analysis of teleost growth hormone genes

A nested polymerase chain reaction (PCR) technique for amplifying a fragment of the gene (GH) encoding teleost growth hormone has been developed. Using this technique, a fragment of the pufferfish, Fugu rubripes and Arothron maculatus; dwarf gourami, Colisa lalia; guppy, Poecilia reticulata; and goldfish, Carassius auratus GH genes were cloned. The Fugu rubripes (Fugu) gene fragment was used to isolate the GH gene from a Fugu genomic library. The complete nucleotide sequence of a 8.5-kb SacI genomic fragment containing the Fugu GH gene has been determined. The GH gene spans 2.5 kb from the first codon to polyadenylation signal, and contains six exons and five introns similar to the GH genes of salmonids, tilapia, barramundi, flounder and yellowtail. The GH introns contain microsatellite and satellite sequences. The microsatellites found in the fifth intron of the GH gene are also present in the corresponding introns of tilapia, barramundi and flounder GH genes. Southern analysis revealed that the GH gene is a single-copy gene in the Fugu. The promoter region of the Fugu GH gene contains conserved sequences that are likely to be involved in the pituitary-specific expression of the gene. A phylogenetic tree of nucleotide (nt) sequences of all known teleost GH genes has been inferred using the distance matrix method. The topology of this tree reflects the major phylogenetic groupings of teleosts. The intron patterns and repetitive sequences of GH genes can serve as useful natural markers for the classification and phylogenetic studies of teleosts.     

Cloning and characterization of ribonuclease T2 gene (RNHe30) from the basidiomycete,Hericium erinaceum

A gene encoding a ribonuclease T2 (RNase T2) family enzyme, RNHe30, was cloned from Hericium erinaceum by PCR. The deduced amino acid sequence from the complimentary DNA (cDNA) (1074 bp) encodes a 302-aa protein (RNase He30) that has the consensus amino acid sequences of RNase T2 family enzymes including the putative signal peptide. The presence of five introns in the genomic DNA was confirmed by comparison of the cDNA and genomic DNA sequences. The promoter region contains a putative CAAT box and a consensus TATA box. Genes coding homologous enzymes were also identified in various other basidiomycetes. A phylogenetic tree of RNase T2s from these fungi was constructed from a multiple alignment of the deduced amino acid sequences. The tree showed that the enzymes were divided into two main groups.     

Structural organization and classification of cytochrome P450 genes in flax (Linum usitatissimum L.)

Flax CYPome analysis resulted in the identification of 334 putative cytochrome P450 (CYP450) genes in the cultivated flax genome. Classification of flax CYP450 genes based on the sequence similarity with Arabidopsis orthologs and CYP450 nomenclature, revealed 10 clans representing 44 families and 98 subfamilies. CYP80, CYP83, CYP92, CYP702, CYP705, CYP708, CYP728, CYP729, CYP733 and CYP736 families are absent in the flax genome. The subfamily members exhibited conserved sequences, length of exons and phasing of introns. Similarity search of the genomic resources of wild flax species Linum bienne with CYP450 coding sequences of the cultivated flax, revealed the presence of 127 CYP450 gene orthologs, indicating amplification of novel CYP450 genes in the cultivated flax. Seven families CYP73, 74, 75, 76, 77, 84 and 709, coding for enzymes associated with phenylpropanoid/fatty acid metabolism, showed extensive gene amplification in the flax. About 59% of the flax CYP450 genes were present in the EST libraries.     

Development of intron length polymorphism markers in cowpea [Vigna unguiculata (L.) Walp.] and their transferability to other Vigna species

Expressed sequence tag (EST) sequences available in the public databases provide a cost-effective and valuable genomic resource for the development of molecular markers. Introns which are non-coding DNA sequences of the gene could be used as potential molecular markers as they are highly variable compared to the coding sequences. This study reports the development of intron length polymorphism markers in cowpea [Vigna unguiculata (L.) Walp.]. The ESTs of cowpea were aligned with genomic sequences of Arabidopsis and soybean to predict the position and number of introns in cowpea. Of the 110 PCR primer pairs designed to amplify the intronic regions, 98 primer pairs resulted in successful amplification and were identified as cowpea intron length polymorphism (CILP) markers. Out of the 45 randomly selected CILP markers, 36?% markers produced length variation in the ten cowpea genotypes, collectively yielding 33 alleles with an average of 2.0 alleles/locus. The polymorphism information content of the CILP markers ranged from 0.18 to 0.64 with an average of 0.34. Of the 98 CILP markers, 93 markers (95?%) showed transferability to other Vigna species. Dendrograms based on CILP markers clearly distinguished the cowpea genotypes as well as other Vigna species, demonstrating the utility of CILP markers in genetic diversity and phylogenetic studies. These CILP markers will be very useful in the genome analysis and marker-assisted breeding of cowpea and other Vigna species.     

Sequence analysis of duplicated actin genes in Lagenidium giganteum and Pythium irregulare (Oomycota)

Southern analysis of genomic DNA identified multiple-copy actin gene families in Lagenidium giganteum and Pythium irregulare (Oomycota). Polymerase chain reaction (PCR) protocols were used to amplify members of these actin gene families. Sequence analysis of genomic coding regions demonstrated five unique actin sequences in L. giganteum (Lg-Ac 1, 2, 3, 4, 5) and four unique actin sequences in P. irregulare (Pi-Acl, 2, 3, 4); none were interrupted by introns. Maximum parsimony analysis of the coding regions demonstrated a close phylogenetic relationship between oomycetes and the chromophyte alga Costaria costata. Three types of actin coding regions were identified in the chromophyte/oomycete lineage. The type 1 actin is the single-copy coding region found in C. costata. The type 2 and type 3 actins are found in the oomycetes and are the result of a gene duplication which occurred soon after the divergence of the oomycetes from the chromophyte algae. The type 2 coding regions are the single-copy sequence of Phytophthora megasperma, the Phytophthora infestans actB gene, Lg-Ac5 and Pi-Ac2. The type 3 coding regions are the single-copy sequence of Achlya bisexualis, the P. infestans actA gene, Lg-Ac1, 2, 3, 4 and Pi-Acl, 3, 4. Correspondence to: D. Bhattacharya     

Phylogenetic position of the Tenrecs (Mammalia: Tenrecidae) of Madagascar based on analysis of the complete mitochondrial genome sequence of Echinops telfairi

The phylogenetic position of the lesser hedgehog tenrec, Echinops telfairi, was studied on the basis of analysis of the concatenated sequences of 12 mitochondrial protein‐coding genes. In addition to the tenrec, the analysis included two other representatives of the insectivore order Lipotyphla, the hedgehog and the mole. The eutherian tree was rooted with three non‐eutherian mammalian taxa. The analysis joined the tenrec, the African elephant (order Proboscidea) and the aardvark (order Tubulidentata) on a common branch. The three lipotyphlan taxa, the tenrec (Tenrecidae), the mole (Talpidae) and the hedgehog (Erinaceidae) were dispersed in the eutherian tree, demonstrating lipotyphlan polyphyly.     

Phylogeny of Cucurbitaceae species in Korea based on 5S rDNA non-transcribed spacer

The 5S rDNA coding region and its spacer have been successfully utilized in phylogenetic studies of plants. However, it has not been utilized in the phylogenetic analysis of Cucurbitaceae. Here, we obtained the 5S rDNA sequences of 12 Cucurbitaceae species by direct PCR or cloning. The 5S rDNA sequences ranged from 275 to 359 bp, and the coding regions of all species were 120 bp long, except for that of Cucurbita, which was 119 bp. Some genus-specific SNPs were observed in the coding regions of Cucurbita, Lagenaria, Melothria, and Tricosanthes. The GC content of the coding regions was generally higher than that of the NTS regions, and the difference in GC content between the coding and NTS regions varied among species, with Gynostemma pentaphyllum having the greatest difference of 20.3. The phylogenetic trees generated using maximum parsimony and maximum likelihood were congruent and well supported by the recently published classification of Cucurbitaceae. These results demonstrated the utility of the 5S rDNA sequence in inferring phylogenetic relationships among 12 Cucurbitaceae species, and its utility could be extended by using a greater number of species in future studies.     

The cephalochordate Branchiostoma genome contains 26 intermediate filament (IF) genes: Implications for evolution of chordate IF proteins

We analyzed the draft genome of the cephalochordate Branchiostoma floridae (B. floridae) for genes encoding intermediate filament (IF) proteins. From 26 identified IF genes 13 were not reported before. Four of the new IF genes belong to the previously established Branchiostoma IF group A, four to the Branchiostoma IF group B, one is homologous to the type II keratin E2 while the remaining four new IF sequences N1 to N4 could not be readily classified in any of the previously established Branchiostoma IF groups. All eleven identified A and B2-type IF genes are located on the same genomic scaffold and arose due to multiple cephalochordate-specific duplications. Another IF gene cluster, identified in the B. floridae genome, contains three keratins (E1, Y1, D1), two keratin-like IF genes (C2, X1), one new IF gene (N1) and one IF unrelated gene, but does not show any similarities to the well defined vertebrate type I or type II keratin gene clusters. In addition, some type III sequence features were documented in the new IF protein N2, which, however, seems to share a common ancestry with the Branchiostoma keratins D1 and two keratin-related genes C. Thus, a few type I and type II keratin genes existed in a common ancestor of cephalochordates and vertebrates, which after separation of these two lineages gave rise to the known complexities of the vertebrate cytoplasmic type I–IV IF proteins, as well as to the multiple keratin and related IF genes in cephalochordates, due to multiple gene duplications, deletions and sequence divergences.     

RNA-Seq quantification of the human small airway epithelium transcriptome

Background

The feline genome is valuable to the veterinary and model organism genomics communities because the cat is an obligate carnivore and a model for endangered felids. The initial public release of the Felis catus genome assembly provided a framework for investigating the genomic basis of feline biology. However, the entire set of protein coding genes has not been elucidated.

Results

We identified and characterized 1227 protein coding feline sequences, of which 913 map to public sequences and 314 are novel. These sequences have been deposited into NCBI's genbank database and complement public genomic resources by providing additional protein coding sequences that fill in some of the gaps in the feline genome assembly. Through functional and comparative genomic analyses, we gained an understanding of the role of these sequences in feline development, nutrition and health. Specifically, we identified 104 orthologs of human genes associated with Mendelian disorders. We detected negative selection within sequences with gene ontology annotations associated with intracellular trafficking, cytoskeleton and muscle functions. We detected relatively less negative selection on protein sequences encoding extracellular networks, apoptotic pathways and mitochondrial gene ontology annotations. Additionally, we characterized feline cDNA sequences that have mouse orthologs associated with clinical, nutritional and developmental phenotypes. Together, this analysis provides an overview of the value of our cDNA sequences and enhances our understanding of how the feline genome is similar to, and different from other mammalian genomes.

Conclusions

The cDNA sequences reported here expand existing feline genomic resources by providing high-quality sequences annotated with comparative genomic information providing functional, clinical, nutritional and orthologous gene information.     

Genome Sequence-Based Discriminator for Vancomycin-Intermediate Staphylococcus aureus

Vancomycin is the mainstay of treatment for patients with Staphylococcus aureus infections, and reduced susceptibility to vancomycin is becoming increasingly common. Accordingly, the development of rapid and accurate assays for the diagnosis of vancomycin-intermediate S. aureus (VISA) will be critical. We developed and applied a genome-based machine-learning approach for discrimination between VISA and vancomycin-susceptible S. aureus (VSSA) using 25 whole-genome sequences. The resulting machine-learning model, based on 14 gene parameters, including 3 molecular typing markers and 11 genes implicated in reduced vancomycin susceptibility, is able to unambiguously distinguish between the VISA and VSSA isolates analyzed here despite the fact that they do not form evolutionarily distinct groups. As such, the model is able to discriminate based on specific genomic markers of antibiotic susceptibility rather than overall sequence relatedness. Subsequent evaluation of the model using leave-one-out validation yielded a classification accuracy of 84%. The machine-learning approach described here provides a generalized framework for the application of genome sequence analysis to the classification of bacteria that differ with respect to clinically relevant phenotypes and should be particularly useful in defining the genomic features that underlie antibiotic resistance.     

Molecular delineation of the Y-borne Sry gene in the Formosan pangolin (Manis pentadactyla pentadactyla) and its phylogenetic implications for Pholidota in extant mammals

The systematic status of Pholidota has been a matter of debate, particularly regarding the apparent inconsistency between morphological and molecular studies. The Sry gene, a master regulator of male sex determination in eutherian mammals, has not yet been used for phylogenetic analyses of extant mammals. The objective of the present study was to clone and characterize the complete gene (1300 base pairs; bp) and amino acid sequences (229 residues) of Sry from the Formosan pangolin (Manis pentadactyla pentadactyla), a member of Pholidota. The Sry amino acid identity between pangolin and other reported species ranged from 42.5% (mouse, Mus musculus) to 84.1% (European hare, Lepus europaeus). Sequence conservation was primarily in the high motility group (HMG) box (234 bp), whereas homology outside the HMG box was low. The cloned Sry was mapped to the pangolin Y chromosome by fluorescence in situ hybridization (FISH); this was confirmed to be the first Y-borne molecular marker identified in Pholidota. Based on Bayesian phylogenetic analysis for Sry HMG sequences from 36 representative taxa, including the Formosan pangolin, Pholidota was more closely related to Carnivora than to Xenarthra, consistent with the emerging molecular tree inferred from markers not located on the Y chromosome. In conclusion, this study characterized the gene structure of Sry of the Formosan pangolin and provided insights into the phylogenetic position of Pholidota.     

Rapid Evolution by Positive Selection and Gene Gain and Loss: PLA<Subscript>2</Subscript> Venom Genes in Closely Related <Emphasis Type="Italic">Sistrurus</Emphasis> Rattlesnakes with Divergent Diets

Rapid evolution of snake venom genes by positive selection has been reported previously but key features of this process such as the targets of selection, rates of gene turnover, and functional diversity of toxins generated remain unclear. This is especially true for closely related species with divergent diets. We describe the evolution of PLA₂ gene sequences isolated from genomic DNA from four taxa of Sistrurus rattlesnakes which feed on different prey. We identified four to seven distinct PLA₂ sequences in each taxon and phylogenetic analyses suggest that these sequences represent a rapidly evolving gene family consisting of both paralogous and homologous loci with high rates of gene gain and loss. Strong positive selection was implicated as a driving force in the evolution of these protein coding sequences. Exons coding for amino acids that make up mature proteins have levels of variation two to three times greater than those of the surrounding noncoding intronic sequences. Maximum likelihood models of coding sequence evolution reveal that a high proportion (∼30%) of all codons in the mature protein fall into a class of codons with an estimated d _N /d _S (ω) ratio of at least 2.8. An analysis of selection on individual codons identified nine residues as being under strong (p < 0.01) positive selection, with a disproportionately high proportion of these residues found in two functional regions of the PLA₂ protein (surface residues and putative anticoagulant region). This is direct evidence that diversifying selection has led to high levels of functional diversity due to structural differences in proteins among these snakes. Overall, our results demonstrate that both gene gain and loss and protein sequence evolution via positive selection are important evolutionary forces driving adaptive divergence in venom proteins in closely related species of venomous snakes.     

A unique DNA repair and recombination gene (recN) sequence for identification and intraspecific molecular typing of bacterial wilt pathogen Ralstonia solanacearum and its comparative analysis with ribosomal DNA sequences

Ribosomal gene sequences are a popular choice for identification of bacterial species and, often, for making phylogenetic interpretations. Although very popular, the sequences of 16S rDNA and 16-23S intergenic sequences often fail to differentiate closely related species of bacteria. The availability of complete genome sequences of bacteria, in the recent years, has accelerated the search for new genome targets for phylogenetic interpretations. The recently published full genome data of nine strains of R. solanacearum, which causes bacterial wilt of crop plants, has provided enormous genomic choices for phylogenetic analysis in this globally important plant pathogen. We have compared a gene candidate recN, which codes for DNA repair and recombination function, with 16S rDNA/16-23S intergenic ribosomal gene sequences for identification and intraspecific phylogenetic interpretations in R. solanacearum. recN gene sequence analysis of R. solanacearum revealed subgroups within phylotypes (or newly proposed species within plant pathogenic genus, Ralstonia), indicating its usefulness for intraspecific genotyping. The taxonomic discriminatory power of recN gene sequence was found to be superior to ribosomal DNA sequences. In all, the recN-sequence-based phylogenetic tree generated with the Bayesian model depicted 21 haplotypes against 15 and 13 haplotypes obtained with 16S rDNA and 16-23S rDNA intergenic sequences, respectively. Besides this, we have observed high percentage of polymorphic sites (S 23.04%), high rate of mutations (Eta 276) and high codon bias index (CBI 0.60), which makes the recN an ideal gene candidate for intraspecific molecular typing of this important plant pathogen.