SVAD: A genetic database curates non-ischemic sudden cardiac death-associated variants.

Sudden cardiac death (SCD) is an important cause of mortality worldwide. It accounts for approximately half of all deaths from cardiovascular disease. While coronary artery disease and acute myocardial infarction account for the majority of SCD in the elderly population, inherited cardiac diseases (inherited CDs) comprise a substantial proportion of younger SCD victims with a significant genetic component. Currently, the use of next-generation sequencing enables the rapid analysis to investigate relationships between genetic variants and inherited CDs causing SCD. Genetic contribution to risk has been considered an alternate predictor of SCD. In the past years, large numbers of SCD susceptibility variants were reported, but these results are scattered in numerous publications. Here, we present the SCD-associated Variants Annotation Database (SVAD) to facilitate the interpretation of variants and to meet the needs of data integration. SVAD contains data from a broad screening of scientific literature. It was constructed to provide a comprehensive collection of genetic variants along with integrated information regarding their effects. At present, SVAD has accumulated 2,292 entries within 1,239 variants by manually surveying pertinent literature, and approximately one-third of the collected variants are pathogenic/likely-pathogenic following the ACMG guidelines. To the best of our knowledge, SVAD is the most comprehensive database that can provide integrated information on the associated variants in various types of inherited CDs. SVAD represents a valuable source of variant information based on scientific literature and benefits clinicians and researchers, and it is now available on http://svad.mbc.nctu.edu.tw/.

Force Dependence of Velocity and Run Length of Kinesin-1, Kinesin-2 and Kinesin-5 Family Molecular Motors.

Kinesin-1, kinesin-2 and kinesin-5 are three families of a superfamily of motor proteins; which can walk processively on microtubule filaments by hydrolyzing ATP. It was experimentally shown that while the three kinesin dimers show similar feature on the force dependence of velocity, they show rather different features on the force dependence of run length. However, why the three families of kinesins show these rather different features is unclear. Here, we computationally studied the movement dynamics of the three dimers based on our proposed model. The simulated results reproduce well the available experimental data on the force dependence of velocity and run length. Moreover, the simulated results on the velocity and run length for the three dimers with altered neck linker lengths are also in quantitative agreement with the available experimental data. The studies indicate that the three families of kinesins show much similar movement mechanism and the rather different features on the force dependence of run length arise mainly from the difference in rate constants of the ATPase activity and neck linker docking. Additionally, the asymmetric (limping) movement dynamics of the three families of homodimers with and without altered neck linker lengths are studied, providing predicted results.

Diversity and enterotype in gut bacterial community of adults in Taiwan.

BACKGROUND: Gastrointestinal microbiota, particularly gut microbiota, is associated with human health. The biodiversity of gut microbiota is affected by ethnicities and environmental factors such as dietary habits or medicine intake, and three enterotypes of the human gut microbiome were announced in 2011. These enterotypes are not significantly correlated with gender, age, or body weight but are influenced by long-term dietary habits. However, to date, only two enterotypes (predominantly consisting of Bacteroides and Prevotella) have shown these characteristics in previous research; the third enterotype remains ambiguous. Understanding the enterotypes can improve the knowledge of the relationship between microbiota and human health. RESULTS: We obtained 181 human fecal samples from adults in Taiwan. Microbiota compositions were analyzed using next-generation sequencing (NGS) technology, which is a culture-independent method of constructing microbial community profiles by sequencing 16S ribosomal DNA (rDNA). In these samples, 17,675,898 sequencing reads were sequenced, and on average, 215 operational taxonomic units (OTUs) were identified for each sample. In this study, the major bacteria in the enterotypes identified from the fecal samples were Bacteroides, Prevotella, and Enterobacteriaceae, and their correlation with dietary habits was confirmed. A microbial interaction network in the gut was observed on the basis of the amount of short-chain fatty acids, pH value of the intestine, and composition of the bacterial community (enterotypes). Finally, a decision tree was derived to provide a predictive model for the three enterotypes. The accuracies of this model in training and independent testing sets were 97.2 and 84.0%, respectively. CONCLUSIONS: We used NGS technology to characterize the microbiota and constructed a predictive model. The most significant finding was that Enterobacteriaceae, the predominant subtype, could be a new subtype of enterotypes in the Asian population.

Water-soluble ionic species of coarse and fine particulate matter and gas precursor characteristics at urban and rural sites of central Taiwan.

Coarse and fine particulate matter (PM) were taken by a dichotomous sampler, and gas precursors were determined by a denuder sampler at two stations in central Taiwan. Water-soluble ionic constituents of PM and their precursor gases were analyzed by ionic chromatograph. In summer, the daytime/nighttime PM10 concentrations were 37 ± 10/41 ± 18 µg m(-3) and 36 ± 14/34 ± 18 µg m(-3) for Xitun and Jhushan, respectively. Average PM10 concentration in winter was 1.55 and 1.76 times that of summer for Xitun and Jhushan, respectively. PM mass concentrations were similar for both stations, although one station is located in the downtown area of Taichung, and the other is in a rural area with no heavy pollution sources. Water-soluble ionic species content was 38-53 % of PM2.5 and 43-48 % of PM10 mass concentration. HNO3, HCl, and SO2 were high in the daytime; the daytime-to-nighttime concentration ratio was 3.75-6.88 for HNO3,1.7-7.8 for HCl, and 1.45-2.77 for SO2. High NH3 levels were determined in the area, especially in winter, which could be a precursor of NH4 (+) to form particulate matter. In Xitun, motor vehicles downtown and in the industrial district could be sources of air pollution. In contrast, there are few industrial sources at Jhushan; therefore, the transport of air pollutants from upwind of other regions and the accumulation of pollutants could be important PM sources at Jhushan.

Oxypurinol-Specific T Cells Possess Preferential TCR Clonotypes and Express Granulysin in Allopurinol-Induced Severe Cutaneous Adverse Reactions.

Allopurinol, a first-line drug for treating gout and hyperuricemia, is one of the leading causes of severe cutaneous adverse reactions (SCARs). To investigate the molecular mechanism of allopurinol-induced SCAR, we enrolled 21 patients (13 Stevens-Johnson syndrome (SJS)/toxic epidermal necrolysis (TEN) and 8 drug reaction with eosinophilia and systemic symptoms (DRESS)), 11 tolerant controls, and 23 healthy donors. We performed in vitro T-cell activation assays by culturing peripheral blood mononuclear cells (PBMCs) with allopurinol, oxypurinol, or febuxostat and measuring the expression of granulysin and IFN-γ in the supernatants of cultures. TCR repertoire was investigated by next-generation sequencing. Oxypurinol stimulation resulted in a significant increase in granulysin in the cultures of blood samples from SCAR patients (n=14) but not tolerant controls (n=11) or healthy donors (n=23). Oxypurinol induced T-cell response in a concentration- and time-dependent manner, whereas allopurinol or febuxostat did not. T cells from patients with allopurinol-SCAR showed no crossreactivity with febuxostat. Preferential TCR-V-ß usage and clonal expansion of specific CDR3 (third complementarity-determining region) were found in the blister cells from skin lesions (n=8) and oxypurinol-activated T-cell cultures (n=4) from patients with allopurinol-SCAR. These data suggest that, in addition to HLA-B*58:01, clonotype-specific T cells expressing granulysin upon oxypurinol induction participate in the pathogenesis of allopurinol-induced SCAR.

Mapping intracellular diffusion distribution using single quantum dot tracking: compartmentalized diffusion defined by endoplasmic reticulum.

The crowded intracellular environment influences the diffusion-mediated cellular processes, such as metabolism, signaling, and transport. The hindered diffusion of macromolecules in heterogeneous cytoplasm has been studied over years, but the detailed diffusion distribution and its origin still remain unclear. Here, we introduce a novel method to map rapidly the diffusion distribution in single cells based on single-particle tracking (SPT) of quantum dots (QDs). The diffusion map reveals the heterogeneous intracellular environment and, more importantly, an unreported compartmentalization of QD diffusions in cytoplasm. Simultaneous observations of QD motion and green fluorescent protein-tagged endoplasmic reticulum (ER) dynamics provide direct evidence that the compartmentalization results from micron-scale domains defined by ER tubules, and ER cisternae form perinuclear areas that restrict QDs to enter. The same phenomenon was observed using fluorescein isothiocyanate-dextrans, further confirming the compartmentalized diffusion. These results shed new light on the diffusive movements of macromolecules in the cell, and the mapping of intracellular diffusion distribution may be used to develop strategies for nanoparticle-based drug deliveries and therapeutics.

Systematic analysis of the association between gut flora and obesity through high-throughput sequencing and bioinformatics approaches.

Eighty-one stool samples from Taiwanese were collected for analysis of the association between the gut flora and obesity. The supervised analysis showed that the most, abundant genera of bacteria in normal samples (from people with a body mass index (BMI) ≤ 24) were Bacteroides (27.7%), Prevotella (19.4%), Escherichia (12%), Phascolarctobacterium (3.9%), and Eubacterium (3.5%). The most abundant genera of bacteria in case samples (with a BMI ≥ 27) were Bacteroides (29%), Prevotella (21%), Escherichia (7.4%), Megamonas (5.1%), and Phascolarctobacterium (3.8%). A principal coordinate analysis (PCoA) demonstrated that normal samples were clustered more compactly than case samples. An unsupervised analysis demonstrated that bacterial communities in the gut were clustered into two main groups: N-like and OB-like groups. Remarkably, most normal samples (78%) were clustered in the N-like group, and most case samples (81%) were clustered in the OB-like group (Fisher's P value = 1.61E - 07). The results showed that bacterial communities in the gut were highly associated with obesity. This is the first study in Taiwan to investigate the association between human gut flora and obesity, and the results provide new insights into the correlation of bacteria with the rising trend in obesity.

Oxaliplatin and its enantiomer induce different condensation dynamics of single DNA molecules.

The interactions of DNA with oxaliplatin (Pt(R,R-DACH)) or its enantiomer (Pt(S,S-DACH)) were investigated using magnetic tweezers and atomic force microscope. In the process of DNA condensation induced by Pt-DACH, only diadducts and micro-loops are formed at low Pt-DACH concentrations, while at high Pt-DACH concentrations, besides the diadducts and micro-loops, long-range cross-links are also formed. The diadduct formation rate of Pt(R,R-DACH) is higher than that of Pt(S,S-DACH). However, the proportions of micro-loops and long-range cross-links for Pt(S,S-DACH) are higher than those for Pt(R,R-DACH). We propose a model to explain these differences between the effect of Pt(R,R-DACH) and that of Pt(S,S-DACH) on DNA condensation. The study has strong implications for the understanding of the effect of chirality on the interaction between Pt-DACH and DNA and the kinetics of DNA condensation induced by platinum complexes.

Sterol regulatory element binding protein 2 activation of NLRP3 inflammasome in endothelium mediates hemodynamic-induced atherosclerosis susceptibility.

BACKGROUND: The molecular basis for the focal nature of atherosclerotic lesions is poorly understood. Here, we explored whether disturbed flow patterns activate an innate immune response to form the NLRP3 inflammasome scaffold in vascular endothelial cells via sterol regulatory element binding protein 2 (SREBP2). METHODS AND RESULTS: Oscillatory flow activates SREBP2 and induces NLRP3 inflammasome in endothelial cells. The underlying mechanisms involve SREBP2 transactivating NADPH oxidase 2 and NLRP3. Consistently, SREBP2, NADPH oxidase 2, and NLRP3 levels were elevated in atheroprone areas of mouse aortas, suggesting that the SREBP2-activated NLRP3 inflammasome causes functionally disturbed endothelium with increased inflammation. Mimicking the effect of atheroprone flow, endothelial cell-specific overexpression of the activated form of SREBP2 synergized with hyperlipidemia to increase atherosclerosis in the atheroresistant areas of mouse aortas. CONCLUSIONS: Atheroprone flow induces NLRP3 inflammasome in endothelium through SREBP2 activation. This increased innate immunity in endothelium synergizes with hyperlipidemia to cause topographical distribution of atherosclerotic lesions.

Transplatin enhances effect of cisplatin on both single DNA molecules and live tumor cells.

Cisplatin is the main platinum antitumor drug applied in clinical settings. However, its trans isomer, transplatin, is known to have an ineffective antitumor activity. Despite intensive studies in this field, the structural and biophysical properties of DNA molecules reacting with these two platinum complexes have not been fully elucidated. In the present study, we observed that transplatin made efficient cross-linking of DNA in the vicinity of cisplatin adducts. High-resolution atomic force microscopy studies revealed that the transplatin-induced cross-linkings of nucleotides flanking cisplatin adducts were characterized by kinked-loop structures with rod-like shapes of nanometer scales (∼10-60nm). The results were further confirmed by denaturing gel electrophoresis and single-molecule experiment using magnetic tweezers. In vivo studies revealed that transplatin and cisplatin co-treatment could induce a considerable amount of kinked loops with smaller sizes (∼15nm) in cellular DNA. Furthermore, compared with cisplatin treatment alone, the co-treatment resulted in enhanced cytotoxicity, increased amount of interstrand cross-links, and cell lesions more reluctant to cellular repair system. The results of the present study provide a new clue for understanding the stepwise reactions of DNA with platinum drugs and might serve as a basis for the development of a new antitumor strategy.

Hypoxia-responsive miRNAs target argonaute 1 to promote angiogenesis.

Despite a general repression of translation under hypoxia, cells selectively upregulate a set of hypoxia-inducible genes. Results from deep sequencing revealed that Let-7 and miR-103/107 are hypoxia-responsive microRNAs (HRMs) that are strongly induced in vascular endothelial cells. In silico bioinformatics and in vitro validation showed that these HRMs are induced by HIF1α and target argonaute 1 (AGO1), which anchors the microRNA-induced silencing complex (miRISC). HRM targeting of AGO1 resulted in the translational desuppression of VEGF mRNA. Inhibition of HRM or overexpression of AGO1 without the 3' untranslated region decreased hypoxia-induced angiogenesis. Conversely, AGO1 knockdown increased angiogenesis under normoxia in vivo. In addition, data from tumor xenografts and human cancer specimens indicate that AGO1-mediated translational desuppression of VEGF may be associated with tumor angiogenesis and poor prognosis. These findings provide evidence for an angiogenic pathway involving HRMs that target AGO1 and suggest that this pathway may be a suitable target for anti- or proangiogenesis strategies.

Effects of paclitaxel on EGFR endocytic trafficking revealed using quantum dot tracking in single cells.

Paclitaxel (PTX), a chemotherapeutic drug, affects microtubule dynamics and influences endocytic trafficking. However, the mechanism and the dynamics of altered endocytic trafficking by paclitaxel treatment in single living cells still remain elusive. By labeling quantum dots (QDs) to the epidermal growth factor (EGF), we continuously tracked the endocytosis and post-endocytic trafficking of EGF receptors (EGFRs) in A549 cells for a long time interval. A single-cell analysis method was introduced to quantitatively study the dynamics of endocytic trafficking. Compared with the control cells, the velocity of directed motion was reduced by 30% due to the suppression of high speed movements of EGF-QDs along the microtubules in PTX-treated cells. The endocytic trafficking in PTX-treated cells was mainly via super-diffusive mode of motion, whereas in control cells, it was mostly via sub-diffusive mode of motion. Moreover, PTX shortened endosomal trafficking and prevented EGF-QDs from moving to the perinuclear area via the rapid delivery of EGF-QDs into the peripheral lysosomes. The present study may shed light on the mechanism of the effect of PTX on the treatment of lung cancer.

miR-103/107 promote metastasis of colorectal cancer by targeting the metastasis suppressors DAPK and KLF4.

Metastasis is the major cause of poor prognosis in colorectal cancer (CRC), and increasing evidence supports the contribution of miRNAs to cancer progression. Here, we found that high expression of miR-103 and miR-107 (miR-103/107) was associated with metastasis potential of CRC cell lines and poor prognosis in patients with CRC. We showed that miR-103/107 targeted the known metastasis suppressors death-associated protein kinase (DAPK) and Krüppel-like factor 4 (KLF4) in CRC cells, resulting in increased cell motility and cell-matrix adhesion and decreased cell-cell adhesion and epithelial marker expression. miR-103/107 expression was increased in the presence of hypoxia, thereby potentiating DAPK and KLF4 downregulation and hypoxia-induced motility and invasiveness. In mouse models of CRC, miR-103/107 overexpression potentiated local invasion and liver metastasis effects, which were suppressed by reexpression of DAPK or KLF4. miR-103/107-mediated downregulation of DAPK and KLF4 also enabled the colonization of CRC cells at a metastatic site. Clinically, the signature of a miR-103/107 high, DAPK low, and KLF4 low expression profile correlated with the extent of lymph node and distant metastasis in patients with CRC and served as a prognostic marker for metastasis recurrence and poor survival. Our findings therefore indicate that miR-103/107-mediated repression of DAPK and KLF4 promotes metastasis in CRC, and this regulatory circuit may contribute in part to hypoxia-stimulated tumor metastasis. Strategies that disrupt this regulation might be developed to block CRC metastasis.

Elastic response and length change of single DNA molecules induced by a combination of cisplatin and transplatin.

Magnetic tweezers were employed to investigate the interactions between DNA and cisplatin (transplatin). Cisplatin shortened DNA and reduced its persistence length more significantly than transplatin due to the formation of many more diadducts than those formed by transplatin. Interestingly, the presence of transplatin could enhance the ability of cisplatin in shortening DNA. An optimal concentration ratio of transplatin to cisplatin existed at which cisplatin showed the highest ability in shortening DNA. Moreover, abrupt length changes were also observed when DNA was treated with a mixture of cisplatin and transplatin at high concentrations. A model was proposed to qualitatively explain well these results.

Palladacycles bearing tridentate CNS-type benzamidinate ligands as catalysts for cross-coupling reactions.

Three pendant benzamidines, [Ph-C(=NC(6)H(5))-{NH(E)}] [E = -(CH(2))(2)SMe (1); -(CH(2))(2)S(t)Bu (2); -o-C(6)H(4)SMe (3)], are described. Reactions of 1, 2 or 3 with one molar equivalent of Pd(OAc)(2) in CH(2)Cl(2) give the palladacyclic complexes, [Ph-C{-NH(η(1)-C(6)H(4))}{=N(E)}]Pd(OAc) [E = -(CH(2))(2)SMe (4); -(CH(2))(2)S(t)Bu (5); -o-C(6)H(4)SMe (6)], as mononuclear palladium complexes respectively. A minor product described as 5', {[Ph-C{-N(C(6)H(5))}{-N(CH(2))(2)S(t)Bu}]Pd(OAc)}(2), was isolated as benzamidinate-bridged dinuclear palladium complex upon recrystallizing from Et(2)O/hexane solution. Treatment of 1, 2 or 3 with one molar equivalent of PdCl(2) in the presence of NEt(3) in CH(2)Cl(2) gives the palladacyclic complexes, [Ph-C{-NH(η(1)-C(6)H(4))}{=N(E)}]PdCl [E = -(CH(2))(2)SMe (7); -(CH(2))(2)S(t)Bu (8); -o-C(6)H(4)SMe (9)], as mononuclear palladium complexes respectively. The crystal and molecular structures are reported for compounds 5, 5' and 6-8. The application of these palladacyclic complexes to the Suzuki and Heck coupling reactions was examined.

Flow-Dependent Regulation of Kruppel-Like Factor 2 Is Mediated by MicroRNA-92a.

BACKGROUND: Upregulated by atheroprotective flow, the transcription factor Krüppel-like factor 2 (KLF2) is crucial for maintaining endothelial function. MicroRNAs (miRNAs) are noncoding small RNAs that regulate gene expression at the posttranscriptional level. We examined the role of miRNAs, particularly miR-92a, in the atheroprotective flow-regulated KLF2. METHODS AND RESULTS: Dicer knockdown increased the level of KLF2 mRNA in human umbilical vein endothelial cells, suggesting that KLF2 is regulated by miRNA. In silico analysis predicted that miR-92a could bind to the 3' untranslated region of KLF2 mRNA. Overexpression of miR-92a decreased the expression of KLF2 and the KLF2-regulated endothelial nitric oxide synthase and thrombomodulin at mRNA and protein levels. A complementary finding is that miR-92a inhibitor increased the mRNA and protein expression of KLF2, endothelial nitric oxide synthase, and thrombomodulin. Subsequent studies revealed that atheroprotective laminar flow downregulated the level of miR-92a precursor to induce KLF2, and the level of this flow-induced KLF2 was reduced by miR-92a precursor. Furthermore, miR-92a level was lower in human umbilical vein endothelial cells exposed to the atheroprotective pulsatile shear flow than under atheroprone oscillatory shear flow. Anti-Ago1/2 immunoprecipitation coupled with real-time polymerase chain reaction revealed that pulsatile shear flow decreased the functional targeting of miR-92a precursor/KLF2 mRNA in human umbilical vein endothelial cells. Consistent with these findings, mouse carotid arteries receiving miR-92a precursor exhibited impaired vasodilatory response to flow. CONCLUSIONS: Atheroprotective flow patterns decrease the level of miR-92a, which in turn increases KLF2 expression to maintain endothelial homeostasis.

Clusters of nucleotide substitutions and insertion/deletion mutations are associated with repeat sequences.

The genome-sequencing gold rush has facilitated the use of comparative genomics to uncover patterns of genome evolution, although their causal mechanisms remain elusive. One such trend, ubiquitous to prokarya and eukarya, is the association of insertion/deletion mutations (indels) with increases in the nucleotide substitution rate extending over hundreds of base pairs. The prevailing hypothesis is that indels are themselves mutagenic agents. Here, we employ population genomics data from Escherichia coli, Saccharomyces paradoxus, and Drosophila to provide evidence suggesting that it is not the indels per se but the sequence in which indels occur that causes the accumulation of nucleotide substitutions. We found that about two-thirds of indels are closely associated with repeat sequences and that repeat sequence abundance could be used to identify regions of elevated sequence diversity, independently of indels. Moreover, the mutational signature of indel-proximal nucleotide substitutions matches that of error-prone DNA polymerases. We propose that repeat sequences promote an increased probability of replication fork arrest, causing the persistent recruitment of error-prone DNA polymerases to specific sequence regions over evolutionary time scales. Experimental measures of the mutation rates of engineered DNA sequences and analyses of experimentally obtained collections of spontaneous mutations provide molecular evidence supporting our hypothesis. This study uncovers a new role for repeat sequences in genome evolution and provides an explanation of how fine-scale sequence contextual effects influence mutation rates and thereby evolution.

Arabidopsis Argonaute 2 regulates innate immunity via miRNA393(∗)-mediated silencing of a Golgi-localized SNARE gene, MEMB12.

Argonaute (AGO) proteins are critical components of RNA silencing pathways that bind small RNAs and mediate gene silencing at their target sites. We found that Arabidopsis AGO2 is highly induced by the bacterial pathogen Pseudomonas syringae pv. tomato (Pst). Further genetic analysis demonstrated that AGO2 functions in antibacterial immunity. One abundant species of AGO2-bound small RNA is miR393b(∗), which targets a Golgi-localized SNARE gene, MEMB12. Pst infection downregulates MEMB12 in a miR393b(∗)-dependent manner. Loss of function of MEMB12, but not SYP61, another intracellular SNARE, leads to increased exocytosis of an antimicrobial pathogenesis-related protein, PR1. Overexpression of miR393b(∗) resembles memb12 mutant in resistance responses. Thus, AGO2 functions in antibacterial immunity by binding miR393b(∗) to modulate exocytosis of antimicrobial PR proteins via MEMB12. Since miR393 also contributes to antibacterial responses, miR393(∗)/miR393 represent an example of a miRNA(∗)/miRNA pair that functions in immunity through two distinct AGOs: miR393(∗) through AGO2 and miR393 through AGO1.

Global analyses of small interfering RNAs derived from Bamboo mosaic virus and its associated satellite RNAs in different plants.

BACKGROUND: Satellite RNAs (satRNAs), virus parasites, are exclusively associated with plant virus infection and have attracted much interest over the last 3 decades. Upon virus infection, virus-specific small interfering RNAs (vsiRNAs) are produced by dicer-like (DCL) endoribonucleases for anti-viral defense. The composition of vsiRNAs has been studied extensively; however, studies of satRNA-derived siRNAs (satsiRNAs) or siRNA profiles after satRNA co-infection are limited. Here, we report on the small RNA profiles associated with infection with Bamboo mosaic virus (BaMV) and its two satellite RNAs (satBaMVs) in Nicotiana benthamiana and Arabidopsis thaliana. METHODOLOGY/PRINCIPAL FINDINGS: Leaves of N. benthamiana or A. thaliana inoculated with water, BaMV alone or co-inoculated with interfering or noninterfering satBaMV were collected for RNA extraction, then large-scale Solexa sequencing. Up to about 20% of total siRNAs as BaMV-specific siRNAs were accumulated in highly susceptible N. benthamiana leaves inoculated with BaMV alone or co-inoculated with noninterfering satBaMV; however, only about 0.1% of vsiRNAs were produced in plants co-infected with interfering satBaMV. The abundant region of siRNA distribution along BaMV and satBaMV genomes differed by host but not by co-infection with satBaMV. Most of the BaMV and satBaMV siRNAs were 21 or 22 nt, of both (+) and (-) polarities; however, a higher proportion of 22-nt BaMV and satBaMV siRNAs were generated in N. benthamiana than in A. thaliana. Furthermore, the proportion of non-viral 24-nt siRNAs was greatly increased in N. benthamiana after virus infection. CONCLUSIONS/SIGNIFICANCE: The overall composition of vsiRNAs and satsiRNAs in the infected plants reflect the combined action of virus, satRNA and different DCLs in host plants. Our findings suggest that the structure and/or sequence demands of various DCLs in different hosts may result in differential susceptibility to the same virus. DCL2 producing 24-nt siRNAs under biotic stresses may play a vital role in the antiviral mechanism in N. benthamiana.

miRExpress: analyzing high-throughput sequencing data for profiling microRNA expression.

BACKGROUND: MicroRNAs (miRNAs), small non-coding RNAs of 19 to 25 nt, play important roles in gene regulation in both animals and plants. In the last few years, the oligonucleotide microarray is one high-throughput and robust method for detecting miRNA expression. However, the approach is restricted to detecting the expression of known miRNAs. Second-generation sequencing is an inexpensive and high-throughput sequencing method. This new method is a promising tool with high sensitivity and specificity and can be used to measure the abundance of small-RNA sequences in a sample. Hence, the expression profiling of miRNAs can involve use of sequencing rather than an oligonucleotide array. Additionally, this method can be adopted to discover novel miRNAs. RESULTS: This work presents a systematic approach, miRExpress, for extracting miRNA expression profiles from sequencing reads obtained by second-generation sequencing technology. A stand-alone software package is implemented for generating miRNA expression profiles from high-throughput sequencing of RNA without the need for sequenced genomes. The software is also a database-supported, efficient and flexible tool for investigating miRNA regulation. Moreover, we demonstrate the utility of miRExpress in extracting miRNA expression profiles from two Illumina data sets constructed for the human and a plant species. CONCLUSION: We develop miRExpress, which is a database-supported, efficient and flexible tool for detecting miRNA expression profile. The analysis of two Illumina data sets constructed from human and plant demonstrate the effectiveness of miRExpress to obtain miRNA expression profiles and show the usability in finding novel miRNAs.