The evolution and structure of snake venom phosphodiesterase (svPDE) highlight its importance in venom actions.

For decades, studies of snake venoms focused on the venom-ome-specific toxins (VSTs). VSTs are dominant soluble proteins believed to contribute to the main venomous effects and emerged into gene clusters for fast adaptation and diversification of snake venoms. However, the conserved minor venom components, such as snake venom phosphodiesterase (svPDE), remain largely unexplored. Here, we focus on svPDE by genomic and transcriptomic analysis across snake clades and demonstrate that soluble svPDE is co-opted from the ancestral membrane-attached ENPP3 (ectonucleotide pyrophosphatase/phosphodiesterase 3) gene by replacing the original 5' exon with the exon encoding a signal peptide. Notably, the exons, promoters, and transcription/translation starts have been replaced multiple times during snake evolution, suggesting the evolutionary necessity of svPDE. The structural and biochemical analyses also show that svPDE shares the similar functions with ENPP family, suggesting its perturbation to the purinergic signaling and insulin transduction in venomous effects.

A genome-wide association study (GWAS) of the personality constructs in CPAI-2 in Taiwanese Hakka populations.

Here in this study we adopted genome-wide association studies (GWAS) to investigate the genetic components of the personality constructs in the Chinese Personality Assessment Inventory 2 (CPAI-2) in Taiwanese Hakka populations, who are likely the descendants of a recent admixture between a group of Chinese immigrants with high emigration intention and a group of the Taiwanese aboriginal population generally without it. A total of 279 qualified participants were examined and genotyped by an Illumina array with 547,644 SNPs to perform the GWAS. Although our sample size is small and that unavoidably limits our statistical power (Type 2 error but not Type 1 error), we still found three genomic regions showing strong association with Enterprise, Diversity, and Logical vs. Affective Orientation, respectively. Multiple genes around the identified regions were reported to be nervous system related, which suggests that genetic variants underlying the certain personalities should indeed exist in the nearby areas. It is likely that the recent immigration and admixture history of the Taiwanese Hakka people created strong linkage disequilibrium between the emigration intention-related genetic variants and their neighboring genetic markers, so that we could identify them despite with only limited statistical power.

MicroRNA retrocopies generated via L1-mediated retrotransposition in placental mammals help to reveal how their parental genes were transcribed.

In mammalian genomes, most retrocopies emerged via the L1 retrotransposition machinery. The hallmarks of an L1-mediated retrocopy, i.e., the intronlessness, the presence of a 3' poly-A tail, and the TSDs at both ends, were frequently used to identify retrotransposition events. However, most previous studies only focused on protein-coding genes as their possible parental sources and thus only a few retrocopies derived from non-coding genes were reported. Remarkably, none of them was from microRNAs. Here in this study, we found several retrocopies generated from the mir-302-367 cluster gene (MIR302CHG), and identified a novel alternatively spliced exon encoding mir-302a. The other recognized microRNA retrotransposition events are primate-specific with mir-373 and mir-498 as their parental genes. The 3' poly-A tracts of these two retrocopy groups were directly attached to the end of the microRNA precursor homologous regions, which suggests that their parental transcripts might alternatively terminate at the end of mir-373 and mir-498. All the three parental microRNAs are highly expressed in specific tissues with elevated retrotransposon activity, such as the embryonic stem cells and the placenta. This might be the reason that our first microRNA retrocopy findings were derived from these three microRNA genes.

Global Investigation of Immune Repertoire Suggests Kawasaki Disease Has Infectious Cause.

BACKGROUND: Kawasaki disease (KD) severely threatens young children's health worldwide. The pathogenic mechanism of KD has not yet been solved, so there is still debate over whether KD is an infectious disease or an autoimmune disease.Methods and Results:To solve this problem, an immune repertoire analysis of KD was conducted. We collected blood cell RNA samples and prepared them into amplicons with iRepertoire kits. The amplicons were sequenced and analyzed with the iRepertoire pipeline. We first identified KD-specific VJ and VDJ forms that had the potential to serve as biomarkers of KD. In addition, the KD-specific VDJ forms were contributed mostly by immunoglobulin G. The D50 value analysis showed that B-cell diversity in KD is decreased, suggesting unique immunoglobulins are produced in KD. Moreover, V, D and J segment usage in IgA, IgG and IgM was consistent with previous KD studies. Further comparison showed no difference in CDR3 peptide length between KD and fever controls (subjects with fever but not diagnosed as KD), indicting KD had B-cell selection phenomenon that has a non-autoimmune pattern. The comparison of amino acid usage of the CDR3 region demonstrated a preference for hydrophilic amino acids in KD. CONCLUSIONS: The results of D50 value, VDJ usage and CDR3 peptide length analyses suggested the characteristics of infectious disease for KD.

Multiomics analyses identified epigenetic modulation of the S100A gene family in Kawasaki disease and their significant involvement in neutrophil transendothelial migration.

BACKGROUND: Kawasaki disease (KD) is a prevalent pediatric disease worldwide and can cause coronary artery aneurysm as a severe complication. Typically, DNA methylation is thought to repress the expression of nearby genes. However, the cases in which DNA methylation promotes gene expression have been reported. In addition, globally, to what extent DNA methylation affects gene expression and how it contributes to the pathogenesis of KD are not yet well understood. METHODS: To address these important biological questions, we enrolled subjects, collected DNA and RNA samples from the subjects' total white blood cells, and performed DNA methylation (M450K) and gene expression (HTA 2.0) microarray assays. RESULTS: By analyzing the variation ratios of CpG beta values (methylation percentage) and gene expression intensities, we first concluded that the CpG markers close (- 1500 bp to + 500 bp) to the transcription start sites had higher variation ratios, reflecting significant regulation capacities. Next, we observed that, globally speaking, gene expression was modestly negatively correlated (correlation rho ≈ - 0.2) with the DNA methylation status of both upstream and downstream CpG markers in the promoter region. Third, we found that specific CpG markers were hypo-methylated in disease samples compared with healthy samples and hyper-methylated in convalescent samples compared with disease samples, promoting and repressing S100A genes' expressions, respectively. Finally, using an in vitro cell model, we demonstrated that S100A family proteins enhanced leukocyte transendothelial migration in KD. CONCLUSIONS: This is the first study to integrate genome-wide DNA methylation with gene expression assays in KD and showed that the S100A family plays important roles in the pathogenesis of KD.

Evolutionary gain of highly divergent tRNA specificities by two isoforms of human histidyl-tRNA synthetase.

The discriminator base N73 is a key identity element of tRNAHis. In eukaryotes, N73 is an "A" in cytoplasmic tRNAHis and a "C" in mitochondrial tRNAHis. We present evidence herein that yeast histidyl-tRNA synthetase (HisRS) recognizes both A73 and C73, but somewhat prefers A73 even within the context of mitochondrial tRNAHis. In contrast, humans possess two distinct yet closely related HisRS homologues, with one encoding the cytoplasmic form (with an extra N-terminal WHEP domain) and the other encoding its mitochondrial counterpart (with an extra N-terminal mitochondrial targeting signal). Despite these two isoforms sharing high sequence similarities (81% identity), they strongly preferred different discriminator bases (A73 or C73). Moreover, only the mitochondrial form recognized the anticodon as a strong identity element. Most intriguingly, swapping the discriminator base between the cytoplasmic and mitochondrial tRNAHis isoacceptors conveniently switched their enzyme preferences. Similarly, swapping seven residues in the active site between the two isoforms readily switched their N73 preferences. This study suggests that the human HisRS genes, while descending from a common ancestor with dual function for both types of tRNAHis, have acquired highly specialized tRNA recognition properties through evolution.

The complete mitochondrial genome of the crocodile shark, Pseudocarcharias kamoharai (Chondrichthyes, Lamnidae).

The complete mitochondrial genome of the crocodile shark consists of 16,688 bp and includes 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, 1 replication origin region, and 1 control region. The mitochondrial gene arrangement of the crocodile shark is the same as that of most vertebrates. Base composition of the genome is A (32.0%), T (31.0%), C (23.7%) and G (13.3%).

Complete mitochondrial genome of the longfin mako shark, Isurus paucus (Chondrichthyes, Lamnidae).

Here we describe the complete mitochondrial genome sequence of the longfin mako, Isurus paucus, which is a pelagic shark found in temperate and tropical waters. The circle genome (16,704 bp) consists of 13 protein coding, 22 tRNA, 2 rRNA genes and 1 control region. It has the typical vertebrate mitochondrial gene arrangement.

The complete mitochondrial genome of the salmon shark, Lamna ditropis (Chondrichthyes, Lamnidae).

The complete mitochondrial genome of the salmon shark consists of 16,699 bp and includes 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, 1 replication origin region and 1 control region. The mitochondrial gene arrangement of the salmon shark is the same as that of most vertebrates. Base composition of the genome is A (29.6%), T (28.6%), C (27.1%), and G (14.8%).

Triple-layer dissection of the lung adenocarcinoma transcriptome: regulation at the gene, transcript, and exon levels.

Lung adenocarcinoma is one of the most deadly human diseases. However, the molecular mechanisms underlying this disease, particularly RNA splicing, have remained underexplored. Here, we report a triple-level (gene-, transcript-, and exon-level) analysis of lung adenocarcinoma transcriptomes from 77 paired tumor and normal tissues, as well as an analysis pipeline to overcome genetic variability for accurate differentiation between tumor and normal tissues. We report three major results. First, more than 5,000 differentially expressed transcripts/exonic regions occur repeatedly in lung adenocarcinoma patients. These transcripts/exonic regions are enriched in nicotine metabolism and ribosomal functions in addition to the pathways enriched for differentially expressed genes (cell cycle, extracellular matrix receptor interaction, and axon guidance). Second, classification models based on rationally selected transcripts or exonic regions can reach accuracies of 0.93 to 1.00 in differentiating tumor from normal tissues. Of the 28 selected exonic regions, 26 regions correspond to alternative exons located in such regulators as tumor suppressor (GDF10), signal receptor (LYVE1), vascular-specific regulator (RASIP1), ubiquitination mediator (RNF5), and transcriptional repressor (TRIM27). Third, classification systems based on 13 to 14 differentially expressed genes yield accuracies near 100%. Genes selected by both detection methods include C16orf59, DAP3, ETV4, GABARAPL1, PPAR, RADIL, RSPO1, SERTM1, SRPK1, ST6GALNAC6, and TNXB. Our findings imply a multilayered lung adenocarcinoma regulome in which transcript-/exon-level regulation may be dissociated from gene-level regulation. Our described method may be used to identify potentially important genes/transcripts/exonic regions for the tumorigenesis of lung adenocarcinoma and to construct accurate tumor vs. normal classification systems for this disease.

Divergent Alanyl-tRNA Synthetase Genes of Vanderwaltozyma polyspora Descended from a Common Ancestor through Whole-Genome Duplication Followed by Asymmetric Evolution.

Cytoplasmic and mitochondrial forms of a eukaryotic aminoacyl-tRNA synthetase (aaRS) are generally encoded by two distinct nuclear genes, one of eukaryotic origin and the other of mitochondrial origin. However, in most known yeasts, only the mitochondrial-origin alanyl-tRNA synthetase (AlaRS) gene is retained and plays a dual-functional role. Here, we present a novel scenario of AlaRS evolution in the yeast Vanderwaltozyma polyspora. V. polyspora possesses two significantly diverged AlaRS gene homologues, one encoding the cytoplasmic form and the other its mitochondrial counterpart. Clever selection of transcription and translation initiation sites enables the two isoforms to be localized and thus functional in their respective cellular compartments. However, the two isoforms can also be stably expressed and function in the reciprocal compartments by insertion or removal of a mitochondrial targeting signal. Synteny and phylogeny analyses revealed that the AlaRS homologues of V. polyspora arose from a dual-functional common ancestor through whole-genome duplication (WGD). Moreover, the mitochondrial form had higher synonymous (1.6-fold) and nonsynonymous (2.8-fold) substitution rates than did its cytoplasmic counterpart, presumably due to a lesser constraint imposed on components of the mitochondrial translational apparatus. Our study suggests that asymmetric evolution confers the divergence between the AlaRS paralogues of V. polyspora.

The complete mitochondrial genome of the sand tiger shark, Carcharias taurus (Chondrichthyes, Odontaspididae).

The complete mitochondrial genome of the sand tiger shark consists of 16,773 bp and including 13 protein-coding genes, 2 ribosomal RNA, 22 transfer RNA genes, 1 replication origin region and 1 control region. The mitochondrial gene arrangement of the sand tiger shark is the same as the one observed in most vertebrates. Base composition of the genome is A (31.8%), T (28.7%), C (26.3%) and G (13.2%).

The complete mitochondrial genome of the Rhodeus shitaiensis (Teleostei, Cypriniformes, Acheilognathidae).

The complete mitochondrial genome of the Rhodeus shitaiensis was determined by using a PCR-based method. The total length of mitochondrial DNA of this bitterling is 16,774 bp and includes 13 protein-coding genes, 2 ribosomal RNA, 22 transfer RNA genes, 1 replication origin region and 1 control region. The mitochondrial gene arrangement of the R. shitaiensis is also matching the one observed in the most vertebrate creatures. Base composition of the genome is A (28.7%), T (26.5%), C (27.4%) and G (17.4%) with an A + T rich hallmark as that of other vertebrate mitochondrial genomes.

The complete mitochondrial genome of the shortfin mako, Isurus oxyrinchus (Chondrichthyes, Lamnidae).

The complete mitochondrial genome of the shortfin mako (Isurus oxyrinchus) was determined by using a PCR-based method. The total length of mitochondrial DNA is 16,701 bp and includes 13 protein-coding genes, 2 ribosomal RNA, 22 transfer RNA genes, 1 replication origin region, and 1 control region. The mitochondrial gene arrangement of the tiger tail seahorse is also matching the one observed in the most vertebrate creatures. Base composition of the genome is A (28.8%), T (28.0%), C (28.0%), and G (15.2%) with an A + T rich hallmark as that of other vertebrate mitochondrial genomes.

Differentiation of visual spectra and nuptial colorations of two Paratanakia himantegus subspecies (Cyprinoidea: Acheilognathidae) in response to the distinct photic conditions of their habitats.

BACKGROUND: Vision, an important sensory modality of many animals, exhibits plasticity in that it adapts to environmental conditions to maintain its sensory efficiency. Nuptial coloration is used to attract mates and hence should be tightly coupled to vision. In Taiwan, two closely related bitterlings (Paratanakia himantegus himantegus and Paratanakia himantegus chii) with different male nuptial colorations reside in different habitats. We compared the visual spectral sensitivities of these subspecies with the ambient light spectra of their habitats to determine whether their visual abilities correspond with photic parameters and correlate with nuptial colorations. RESULTS: Theelectroretinogram (ERG) results revealed that the relative spectral sensitivity of P.h. himantegus was higher at 670 nm, but lower at 370 nm, than the sensitivity of P. h. chii. Both bitterlings could perceive and reflect UV light, but the UV reflection patterns differed between genders. Furthermore, the relative irradiance intensity of the light spectra in the habitat of P. h. himantegus was higher at long wavelengths (480-700 nm), but lower at short wavelengths (350-450 nm), than the light spectra in the habitats of P. h.chii. CONCLUSIONS: Two phylogenetically closely related bitterlings, P. h. himantegus and P. h. chii, dwell in different waters and exhibit different nuptial colorations and spectral sensitivities, which may be the results of speciation by sensory drive. Sensory ability and signal diversity accommodating photic environment may promote diversity of bitterling fishes. UV light was demonstrated to be a possible component of bitterling visual communication. The UV cue may assist bitterlings in genderidentification.

Phylogenetic relationships of Acheilognathidae (Cypriniformes: Cyprinoidea) as revealed from evidence of both nuclear and mitochondrial gene sequence variation: evidence for necessary taxonomic revision in the family and the identification of cryptic species.

Bitterlings are relatively small cypriniform species and extremely interesting evolutionarily due to their unusual reproductive behaviors and their coevolutionary relationships with freshwater mussels. As a group, they have attracted a great deal of attention in biological studies. Understanding the origin and evolution of their mating system demands a well-corroborated hypothesis of their evolutionary relationships. In this study, we provide the most comprehensive phylogenetic reconstruction of species relationships of the group based on partitioned maximum likelihood and Bayesian methods using DNA sequence variation of nuclear and mitochondrial genes on 41 species, several subspecies and three undescribed species. Our findings support the monophyly of the Acheilognathidae. Two of the three currently recognized genera are not monophyletic and the family can be subdivided into six clades. These clades are further regarded as genera based on both their phylogenetic relationships and a reappraisal of morphological characters. We present a revised classification for the Acheilognathidae with five genera/lineages: Rhodeus, Acheilognathus (new constitution), Tanakia (new constitution), Paratanakia gen. nov., and Pseudorhodeus gen. nov. and an unnamed clade containing five species currently referred to as "Acheilognathus". Gene trees of several bitterling species indicate that the taxa are not monophyletic. This result highlights a potentially dramatic underestimation of species diversity in this family. Using our new phylogenetic framework, we discuss the evolution of the Acheilognathidae relative to classification, taxonomy and biogeography.

The complete mitochondrial genome of the great white shark, Carcharodon carcharias (Chondrichthyes, Lamnidae).

The complete mitochondrial genome of the great white shark having 16,744 bp and including 13 protein-coding genes, 2 ribosomal RNA, 22 transfer RNA genes, 1 replication origin region and 1 control region. The mitochondrial gene arrangement of the great white shark is the same as the one observed in the most vertebrates. Base composition of the genome is A (30.6%), T (28.7%), C (26.9%) and G (13.9%).

The complete mitochondrial genome of the big-eye thresher shark, Alopias superciliosus (Chondrichthyes, Alopiidae).

The complete mitochondrial genome of the big-eye thresher shark was sequenced using a polymerase chain reaction (PCR)-based method. The total length of mitochondrial DNA is 16,719 bp and includes 13 protein-coding genes, 2 ribosomal RNA, 22 transfer RNA genes, 1 replication origin region and 1 control region. The mitochondrial gene arrangement of the big-eye thresher shark is the same as the one observed in the most vertebrates. Base composition of the genome is A (31.8%), T (28.9%), C (25.8%) and G (13.5%).

Complete mitochondrial genome of the megamouth shark Megachasma pelagios (Chondrichthyes, Megachasmidae).

Here we describe the complete mitochondrial genome sequence of the megamouth shark, Megachasma pelagios, which is an extremely rare species of deepwater shark. The circle genome (16,694 bp) consists of 13 protein coding, 22 tRNA, 2 rRNA genes and 1 control region. It has the typical vertebrate mitochondrial gene arrangement.

The complete mitochondrial genome of the three-spot seahorse, Hippocampus trimaculatus (Teleostei, Syngnathidae).

The complete mitochondrial genome of the three-spot seahorse was sequenced using a polymerase chain reaction-based method. The total length of mitochondrial DNA is 16,535 bp and includes 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region. The mitochondrial gene order of the three-spot seahorse also conforms to the distinctive vertebrate mitochondrial gene order. The base composition of the genome is A (32.7%), T (29.3%), C (23.4%), and G (14.6%) with an A + T-rich hallmark as that of other vertebrate mitochondrial genomes.