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.

Evaluation and application of the strand-specific protocol for next-generation sequencing.

Next-generation sequencing (NGS) has become a powerful sequencing tool, applied in a wide range of biological studies. However, the traditional sample preparation protocol for NGS is non-strand-specific (NSS), leading to biased estimates of expression for transcripts overlapped at the antisense strand. Strand-specific (SS) protocols have recently been developed. In this study, we prepared the same RNA sample by using the SS and NSS protocols, followed by sequencing with Illumina HiSeq platform. Using real-time quantitative PCR as a standard, we first proved that the SS protocol more precisely estimates gene expressions compared with the NSS protocol, particularly for those overlapped at the antisense strand. In addition, we also showed that the sequence reads from the SS protocol are comparable with those from conventional NSS protocols in many aspects. Finally, we also mapped a fraction of sequence reads back to the antisense strand of the known genes, originally without annotated genes located. Using sequence assembly and PCR validation, we succeeded in identifying and characterizing the novel antisense genes. Our results show that the SS protocol performs more accurately than the traditional NSS protocol and can be applied in future studies.

miRSeq: a user-friendly standalone toolkit for sequencing quality evaluation and miRNA profiling.

MicroRNAs (miRNAs) present diverse regulatory functions in a wide range of biological activities. Studies on miRNA functions generally depend on determining miRNA expression profiles between libraries by using a next-generation sequencing (NGS) platform. Currently, several online web services are developed to provide small RNA NGS data analysis. However, the submission of large amounts of NGS data, conversion of data format, and limited availability of species bring problems. In this study, we developed miRSeq to provide alternatives. To test the performance, we had small RNA NGS data from four species, including human, rat, fly, and nematode, analyzed with miRSeq. The alignments results indicate that miRSeq can precisely evaluate the sequencing quality of samples regarding percentage of self-ligation read, read length distribution, and read category. miRSeq is a user-friendly standalone toolkit featuring a graphical user interface (GUI). After a simple installation, users can easily operate miRSeq on a PC or laptop by using a mouse. Within minutes, miRSeq yields useful miRNA data, including miRNA expression profiles, 3' end modification patterns, and isomiR forms. Moreover, miRSeq supports the analysis of up to 105 animal species, providing higher flexibility.

Functional evolution of cardiac microRNAs in heart development and functions.

MicroRNAs (miRNAs) are a class of endogenous small noncoding RNAs that regulate gene expression either by degrading target mRNAs or by suppressing protein translation. miRNAs have been found to be involved in many biological processes, such as development, differentiation, and growth. However, the evolution of miRNA regulatory functions and networks has not been well studied. In this study, we conducted a cross-species analysis to study the evolution of cardiac miRNAs and their regulatory functions and networks. We found that conserved cardiac miRNA target genes have maintained highly conserved cardiac functions. Additionally, most of cardiac miRNA target genes in human with annotations of cardiac functions evolved from the corresponding homologous targets, which are also involved in heart development-related functions. On the basis of these results, we investigated the functional evolution of cardiac miRNAs and presented a functional evolutionary map. From this map, we identified the evolutionary time at which the cardiac miRNAs became involved in heart development or function and found that the biological processes of heart development evolved earlier than those of heart functions, for example, heart contraction/relaxation or cardiac hypertrophy. Our study of the evolution of the cardiac miRNA regulatory networks revealed the emergence of new regulatory functional branches during evolution. Furthermore, we discovered that early evolved cardiac miRNA target genes tend to participate in the early stages of heart development. This study sheds light on the evolution of developmental features of genes regulated by cardiac miRNAs.

A simple method using PyrosequencingTM to identify de novo SNPs in pooled DNA samples.

A practical way to reduce the cost of surveying single-nucleotide polymorphism (SNP) in a large number of individuals is to measure the allele frequencies in pooled DNA samples. Pyrosequencing(TM) has been frequently used for this application because signals generated by this approach are proportional to the amount of DNA templates. The Pyrosequencing(TM) pyrogram is determined by the dispensing order of dNTPs, which is usually designed based on the known SNPs to avoid asynchronistic extensions of heterozygous sequences. Therefore, utilizing the pyrogram signals to identify de novo SNPs in DNA pools has never been undertook. Here, in this study we developed an algorithm to address this issue. With the sequence and pyrogram of the wild-type allele known in advance, we could use the pyrogram obtained from the pooled DNA sample to predict the sequence of the unknown mutant allele (de novo SNP) and estimate its allele frequency. Both computational simulation and experimental Pyrosequencing(TM) test results suggested that our method performs well. The web interface of our method is available at http://life.nctu.edu.tw/∼yslin/PSM/.