Search

Bunostomum trigonocephalum and Bunostomum phlebotomum are blood-feeding hookworms of sheep and cattle, causing considerable economic losses to the live stock industries. Studying genetic variability within and among hookworm populations is critical to addressing epidemiological and ecological questions. Mitochondrial (mt) DNA is known to provide useful markers for investigations of population genetics of hookworms, but mt genome sequence data are scant. In the present study, the complete mitochondrial DNA (mtDNA) sequences of the sheep and goat hookworm B. trigonocephalum were determined for the first time, and the mt genome of B. phlebotomum from yak in China was also sequenced for comparative analyses of their gene contents and genome organizations. The lengths of mt DNA sequences of B. trigonocephalum sheep isolate, B. trigonocephalum goat isolate and B. phlebotomum China yak isolate were 13,764bp, 13,771bp and 13,803bp in size, respectively. The identity of the mt genomes was 99.7% between B. trigonocephalum sheep isolate and B. trigonocephalum goat isolate. The identity of B. phlebotomum China yak isolate mt genomes was 85.3% with B. trigonocephalum sheep isolate, and 85.2% with B. trigonocephalum goat isolate. All the mt genes of the two hookworms were transcribed in the same direction and gene arrangements were consistent with those of the GA3 type, including 12 protein-coding genes, 2 rRNA genes and 22 tRNA genes, but lacking ATP synthetase subunit 8 gene. The mt genomes of B. trigonocephalum and B. phlebotomum were similar to prefer bases A and T, the contents of A+T are 76.5% (sheep isolate), 76.4% (goat isolate) and 76.9% (China yak isolate), respectively. Phylogenetic relationships reconstructed using concatenated amino acid sequences of 12 protein-coding genes with three methods (maximum likelihood, Bayesian inference and neighbor joining) revealed that the B. trigonocephalum and B. phlebotomum represent distinct but closely-related species. These data provide novel and useful genetic markers for studying the systematics, and population genetics of the two ruminant hookworms.

Freshwater fishes in India are poorly known and plagued by many unresolved cryptic species complexes that masks some latent and endemic species. Limitations in traditional taxonomy have resulted in this crypticism. Hence, molecular approaches like DNA barcoding, are needed to diagnose these latent species. We have analyzed 1383 barcode sequences of 175 Indian freshwater fish species available in the databases, of which 172 sequences of 70 species were generated. The congeneric and conspecific genetic divergences were calculated using Kimura's 2 parameter distance model followed by the construction of a Neighbor Joining tree using the MEGA 5.1. DNA barcoding principle at its first hand approach, led to the straightforward identification of 82% of the studied species with 2.9% (S.E=0.2) divergence between the nearest congeners. However, after validating some cases of synonymy and mislabeled sequences, 5% more species were found to be valid. Sequences submitted to the database under different names were found to represent single species. On the other hand, some sequences of the species like Barilius barna, Barilius bendelisis and Labeo bata were submitted to the database under a single name but were found to represent either some unexplored species or latent species. Overall, 87% of the available Indian freshwater fish barcodes were diagnosed as true species in parity with the existing checklist and can act as reference barcode for the particular taxa. For the remaining 13% (21 species) the correct species name was difficult to assign as they depicted some erroneous identification and cryptic species complex. Thus, these barcodes will need further assay and inclusion of barcodes of more specimens from same and sister species.

Many microRNAs (miRNAs) are clustered on chromosomes and co-transcribed as polycistronic transcripts. Here, an integrated evolutionary analysis of human miRNA gene clusters and families was performed. Generally, miRNA gene clusters include 2-8 members, but some larger clusters have been found to have more members (over 40 miRNAs). 62.22% of them have been shown to be involved in homologous miRNA genes, including multicopy pre-miRNAs and sense/antisense homologous miRNAs. Multicopy pre-miRNAs can enrich the distribution and relationship between miRNA clusters and families. An miRNA family may be located in one or more clusters, and a cluster may be involved in one or more families. Members of different families have been shown to be prone to appear in clusters, and vice versa. Reconstructed phylogenetic trees and networks may indicate potential evolutionary relationships, which also indicate duplication history in specific related gene clusters and families. Related miRNA families are always found to share common target mRNAs and biological pathways. Some clusters containing non-homologous miRNAs also tend to be clustered together as well as homologous miRNAs. In the present work, it is shown that homologous miRNAs are prone to appear in clusters based on functional and evolutionary pressures. The phenomenon of miRNA clusters containing homologous or genetic relationships is quite common. The integrative evolutionary analysis will provide more potential evolutionary and functional relationships between homologous and clustered miRNAs.

An effective DNA marker for authenticating the genus Salvia was screened using seven DNA regions (rbcL, matK, trnL-F, and psbA-trnH from the chloroplast genome, and ITS, ITS1, and ITS2 from the nuclear genome) and three combinations (rbcL+matK, psbA-trnH+ITS1, and trnL-F+ITS1). The present study collected 232 sequences from 27 Salvia species through DNA sequencing and 77 sequences within the same taxa from the GenBank. The discriminatory capabilities of these regions were evaluated in terms of PCR amplification success, intraspecific and interspecific divergence, DNA barcoding gaps, and identification efficiency via a tree-based method. ITS1 was superior to the other marker for discriminating between species, with an accuracy of 81.48%. The three combinations did not increase species discrimination. Finally, we found that ITS1 is a powerful barcode for identifying Salvia species, especially Salvia miltiorrhiza.

At an average altitude of 4000m above sea level, the Tibetan plateau is one of the highest plains on the planet. It is surrounded on three sides by massive mountain ranges: the Kunlun, the Karakoram and the Himalayas. These natural barriers have kept Tibet relatively isolated. In the present study, 15 autosomal STR loci were genotyped in 338 unrelated individuals from three traditional provinces of historical Tibet: Amdo (86), Kham (101) and U-Tsang (151). All the studied loci were in Hardy-Weinberg equilibrium except for the D19S433 locus in the Kham province. FGA, D21S11 and D2S1338 show the highest observed heterozygosity values in Amdo (0.8954), Kham (0.9208) and U-Tsang (0.8940), respectively, whereas TPOX is the least variable marker displaying the lowest value for the same parameter. U-Tsang exhibits the highest total numbers of alleles (139) followed by Kham (130) and Amdo (128) groups. The allele frequency data from this study were compared to relevant global reference populations. Our results indicate that although these three Tibetan populations group together in both the Correspondence Analysis (CA) plot and the Neighbor Joining (NJ) tree, they exhibit some degree of genetic differentiation among themselves congruent with their unique dialects, cultures and traditions. The 15 autosomal STR loci studied were found to be informative and discriminating, thereby providing a useful set of markers for population genetic studies.

We analyzed 22 clinical isolates of Plasmodium vivax from Thailand and 17 from Brazil to investigate the extent of sequence variation in the thrombospondin-related adhesive protein of Plasmodium vivax (PvTRAP), a homologue of P. falciparum TRAP (PfTRAP) which has been considered to be a promising vaccine candidate. In total 54 haplotypes were identified from 73 distinct gene clones. Coexistence of different PvTRAP in circulation occurred in 10 and 13 isolates from Thailand and Brazil, respectively. Forty out of 48 substituted nucleotides are non-synonymous changes. Most of the substituted residues reside in the von Willebrand factor type A-domain (region II), a sulfated glycosaminoglycan-binding domain (region III) and a proline-rich region (region IV). All nucleotide substitutions are dimorphic. Two haplotypes from Thailand contain an inserted sequence encoding aspartic acid-serine-proline in the proline-rich region. Sequence analysis has revealed that nucleotide diversity in PvTRAP is low although Brazilian isolates display a higher degree of variation than those from Thailand. Phylogenetic construction using the neighbor joining method has shown that most of the Thai and the Brazilian isolates appear to be mainly clustered into distinct groups. Significantly greater than expected values of the mean number of non-synonymous (d(n)) than synonymous (d(s)) nucleotide substitutions per site were observed in regions II and III of PvTRAP. Analysis of the published PfTRAP sequences has shown a similar finding in regions II and IV suggesting that positive selection operates on the regions. Hence, different regions in PvTRAP and PfTRAP could be under different pressures in terms of immune selection, structural and/or functional constraints.

The halotolerant alkane-assimilating yeast Debaryomyces hansenii was examined for P450 alkane hydroxylase genes known to be required for alkane assimilation in Candida. Four distinct P450alk gene segments and an allelic segment were isolated using PCR based on degenerate primers derived from the CYP52 family of alkane-inducible P450 genes. A screen of a genomic library (15–20 kb inserts) constructed for this study, using a probe based on the PCR-isolated segments, yielded seven clones. This has led to the isolation and sequence of two full-length genes DH-ALK1 and DH-ALK2. These genes, each with an ORF of 1557 bp (519 aa), contained no apparent introns and showed 64% nucleotide sequence homology (61% based on the deduced amino acid sequences). The deduced proteins had predicted molecular weights of 59,254 Da (DH-ALK1) and 59,614 Da (DH-ALK2) and have been designated CYP52A12 and CYP52A13 by the P450 Nomenclature Committee. Phylogenetic analysis based on Neighbor Joining Tree showed that DH-ALK1 and DH-ALK2 constitute new genes located on two distinct branches and are most related to the gene CYP52A3 (60% deduced aa homology) and are least related to the gene CYP52C2 (41% deduced aa homology), both of C. maltosa. The isolated genes will provide tools to better understand the diversity of the P450alk family in eukaryotic microorganisms adapted to varied environmental conditions.

MicroRNAs (miRNAs) are crucial negative regulators of gene expression at the post-transcriptional level. Next-generation sequencing technologies have identified a series of miRNA variants (named isomiRs). In this study, paralogous isomiR assemblies (from the miRNA locus) were systematically analyzed based on data acquired from deep sequencing data sets. Evolutionary analysis of paralogous (members in miRNA gene family in a specific species) and orthologues (across different animal species) miRNAs was also performed. The sequence diversity of paralogous isomiRs was found to be similar to the diversity of paralogous and orthologues miRNAs. Additionally, both isomiRs and paralogous/orthologues miRNAs were implicated in 5' and 3' ends (especially 3' ends), nucleotide substitutions, and insertions and deletions. Generally, multiple isomiRs can be produced from a single miRNA locus, but most of them had lower enrichment levels, and only several dominant isomiR sequences were detected. These dominant isomiR groups were always stable, and one of them would be selected as the most abundant miRNA sequence in specific animal species. Some isomiRs might be consistent to miRNA sequences in some species but not the other. Homologous miRNAs were often detected in similar isomiR repertoires, and showed similar expression patterns, while dominant isomiRs showed complex evolutionary patterns from miRNA sequences across the animal kingdom. These results indicate that the phenomenon of multiple isomiRs is not a random event, but rather the result of evolutionary pressures. The existence of multiple isomiRs enables different species to express advantageous sequences in different environments. Thus, dominant sequences emerge in response to functional and evolutionary pressures, allowing an organism to adapt to complex intra- and extra-cellular events.

TRAF3 interacting protein 2 (TRAF3IP2) is important for immune responses to pathogens, inflammatory signals and autoimmunity in mammals. In the present study, we collected 19 mammalian TRAF3IP2 sequences and investigated the various types of selection pressure acting on them. Maximum likelihood estimations of nonsynonymous (dN) to synonymous (dS) substitution (dN/dS) ratios for the aligned coding sequences indicated that, as a whole, TRAF3IP2 has been subject to purifying selection. However, the N-terminus of the protein has been subject to higher selection pressure than the C-terminal domain. While eight amino acid residues within the N-terminus appear to have evolved under positive selection, no evidence for such selection was found in the C-terminus. The positively selected residues, which fall outside the currently known functional sites within TRAF3IP2, may have novel functions. The different selection pressures acting on the N- and C-terminal regions are consistent with their protein structures: the C-terminal structure is an ordered structure, whereas the N-terminus is disordered. Taken together with the results of previous studies, it is plausible that positive selection on the N-terminus of TRAF3IP2 may have occurred by competitive coevolution between mammalian hosts and viruses.

Disintegrins are small, non-enzymatic proteins produced in snake venom. PCR and DNA sequencing analysis of genomic DNA for all subspecies of the copperhead snake (Agkistrodon contortrix) were analyzed for the presence of a disintegrin gene. Four samples each of the subspecies: A. c. contortrix, A. c. laticinctus, A. c. mokasen, A. c. phaeogaster, and A. c. pictigaster were collected from different locations across their geographic range and analyzed. A single PCR fragment from each sample was obtained, containing exon and intron sequences. The disintegrins identified in this study shared the highest amino acid identity to contortrostatin and acostatin b chain. Neighbor joining analysis of the disintegrin haplotypes and bootstrap tests of significance grouped the A. contortrix subspecies into two clades. The A. c. mokasen samples collected in Kentucky were grouped in one clade, while the A. c. contortrix, A. c. laticinctus, A. c. phaeogaster, and A. c. pictigaster samples collected in Texas, Louisiana, and Missouri were grouped in a different clade. Analysis of molecular variance (AMOVA) and PhiST pairwise comparisons showed significant genetic variation between subspecies. Nucleotide substitution analysis suggests the rapid evolution of disintegrin genes in A. contortrix subspecies.