核仁小RNA 在肿瘤发生中的作用

核仁小RNA（small nucleolar RNA,snoRNA）是一类真核细胞核仁中的60～300个核苷酸长度的非编码RNA,主要参与rRNA和其它小RNA转录后的成熟加工过程. 它们与肿瘤的关系曾一度被人们所忽视,然而,近年来有关snoRNA新功能的研究证明,它们与肿瘤的发生、发展密切相关. snoRNA以多种方式参与肿瘤的发生：一些snoRNA（如：U50、SNORD12、SNORD12b、SNORD12c、SNORD44和h5sn2等）具有抑癌活性,而另一些snoRNA（如：SNORD33、SNORD66、SNORD76、SNORD112、SNORD113、SNORD114、SNORA42、U70C和ACA59B等）具有促癌活性. 另外,编码snoRNA基因的异常也被发现与肿瘤的发生有关. 因此,开展snoRNA与肿瘤关系的研究将有可能为肿瘤诊治提供新线索.     

The SNORD115 (H/MBII-52) and SNORD116 (H/MBII-85) gene clusters at the imprinted Prader-Willi locus generate canonical box C/D snoRNAs

The imprinted Snurf-Snrpn chromosomal domain contains two large arrays of tandemly repeated, paternally expressed box C/D small-nucleolar RNA (snoRNA) genes: the SNORD115 (H/MBII-52) and SNORD116 (H/MBII-85) gene clusters believed to play key roles in the fine-tuning of serotonin receptor (5-HT2C) pre-mRNA processing and in the etiology of the Prader-Willi Syndrome (PWS), respectively. SNORD115 and SNORD116 were recently proposed to undergo significant conversion into shorter RNA species, the so-called psnoRNAs. Here, we provide evidence that argues against the existence of abundant psnoRNAs in human or mouse brain. Instead, we characterize a previously unsuspected low-abundance, fibrillarin-associated SNORD115-derived smaller RNA species. Based on these findings, we strongly recommend that PWS-encoded SNORD115 and SNORD116 be considered as bona fide box C/D snoRNAs.     

RNA-Seq of the Nucleolus Reveals Abundant SNORD44-Derived Small RNAs

Small non-coding RNAs represent RNA species that are not translated to proteins, but which have diverse and broad functional activities in physiological and pathophysiological states. The knowledge of these small RNAs is rapidly expanding in part through the use of massive parallel (deep) sequencing efforts. We present here the first deep sequencing of small RNomes in subcellular compartments with particular emphasis on small RNAs (sRNA) associated with the nucleolus. The vast majority of the cellular, cytoplasmic and nuclear sRNAs were identified as miRNAs. In contrast, the nucleolar sRNAs had a unique size distribution consisting of 19–20 and 25 nt RNAs, which were predominantly composed of small snoRNA-derived box C/D RNAs (termed as sdRNA). Sequences from 47 sdRNAs were identified, which mapped to both 5′ and 3′ ends of the snoRNAs, and retained conserved box C or D motifs. SdRNA reads mapping to SNORD44 comprised 74% of all nucleolar sdRNAs, and were confirmed by Northern blotting as comprising both 20 and 25 nt RNAs. A novel 120 nt SNORD44 form was also identified. The expression of the SNORD44 sdRNA and 120 nt form was independent of Dicer/Drosha–mediated processing pathways but was dependent on the box C/D snoRNP proteins/sno-ribonucleoproteins fibrillarin and NOP58. The 120 nt SNORD44-derived RNA bound to fibrillarin suggesting that C/D sno-ribonucleoproteins are involved in regulating the stability or processing of SNORD44. This study reveals sRNA cell-compartment specific expression and the distinctive unique composition of the nucleolar sRNAs.     

The short non-coding transcriptome of the protozoan parasite Trypanosoma cruzi

The pathway for RNA interference is widespread in metazoans and participates in numerous cellular tasks, from gene silencing to chromatin remodeling and protection against retrotransposition. The unicellular eukaryote Trypanosoma cruzi is missing the canonical RNAi pathway and is unable to induce RNAi-related processes. To further understand alternative RNA pathways operating in this organism, we have performed deep sequencing and genome-wide analyses of a size-fractioned cDNA library (16-61 nt) from the epimastigote life stage. Deep sequencing generated 582,243 short sequences of which 91% could be aligned with the genome sequence. About 95-98% of the aligned data (depending on the haplotype) corresponded to small RNAs derived from tRNAs, rRNAs, snRNAs and snoRNAs. The largest class consisted of tRNA-derived small RNAs which primarily originated from the 3' end of tRNAs, followed by small RNAs derived from rRNA. The remaining sequences revealed the presence of 92 novel transcribed loci, of which 79 did not show homology to known RNA classes.     

Global gene methylation profiling of common warts caused by human papillomaviruses infection

Infection with the human papillomaviruses (HPV) often involves the epigenetic modification of the host genome. Despite its prevalence among the population, host genome methylation in HPV-induced warts is not clearly understood. In this study, genome-wide methylation profiling was carried out on paired healthy skin and wart samples in order to investigate the effects that benign HPV infection has on gene methylation status. To overcome this gap in knowledge, paired wart (n = 12) and normal skin (n = 12) samples were obtained from Arab males in order to perform DNA extraction and subsequent genome-wide methylation profiling on the Infinium Methylation EPIC Bead Chip microarray. Analysis of differential methylation revealed a clear pattern of discrimination between the wart and normal skin samples. In warts, the most differentially methylated (DM) genes included long non-coding RNAs (AC005884, AL049646.2, AC126121.2, AP001790.1, and AC107959.3), microRNAs (MIR374B, MIR596, MIR1255B1, MIR26B, and MIR196A2),snoRNAs (SNORD114-22, SNORD70, and SNORD114-31), pseudogenes (AC069366.1, RNU4ATAC11P, AC120057.1, NANOGP3, AC106038.2, TPT1P2, SDC4P, PKMP3, and VN2R3P), and protein-coding genes (AREG, GJB2, C12orf71, AC020909.2, S100A8, ZBED2, FABP7, and CYSLTR1). In addition, pathway analysis revealed that, among the most differentially methylated genes, STAT5A, RARA, MEF2D, MAP3K8, and THRA were the common regulators. It can be observed that HPV-induced warts involve a clear and unique epigenetic alteration to the host genome.     

Direct cloning of double-stranded RNAs from RNase protection analysis reveals processing patterns of C/D box snoRNAs and provides evidence for widespread antisense transcript expression

We describe a new method that allows cloning of double-stranded RNAs (dsRNAs) that are generated in RNase protection experiments. We demonstrate that the mouse C/D box snoRNA MBII-85 (SNORD116) is processed into at least five shorter RNAs using processing sites near known functional elements of C/D box snoRNAs. Surprisingly, the majority of cloned RNAs from RNase protection experiments were derived from endogenous cellular RNA, indicating widespread antisense expression. The cloned dsRNAs could be mapped to genome areas that show RNA expression on both DNA strands and partially overlapped with experimentally determined argonaute-binding sites. The data suggest a conserved processing pattern for some C/D box snoRNAs and abundant expression of longer, non-coding RNAs in the cell that can potentially form dsRNAs.     

The GPSM2/LGN GoLoco motifs are essential for hearing

The planar cell polarity (PCP) pathway is responsible for polarizing and orienting cochlear hair cells during development through movement of a primary cilium, the kinocilium. GPSM2/LGN, a mitotic spindle-orienting protein associated with deafness in humans, is a PCP effector involved in kinocilium migration. Here, we link human and mouse truncating mutations in the GPSM2/LGN gene, both leading to hearing loss. The human variant, p.(Trp326*), was identified by targeted genomic enrichment of genes associated with deafness, followed by massively parallel sequencing. Lgn ^ΔC mice, with a targeted deletion truncating the C-terminal GoLoco motifs, are profoundly deaf and show misorientation of the hair bundle and severe malformations in stereocilia shape that deteriorates over time. Full-length protein levels are greatly reduced in mutant mice, with upregulated mRNA levels. The truncated Lgn ^ΔC allele is translated in vitro, suggesting that mutant mice may have partially functioning Lgn. Gαi and aPKC, known to function in the same pathway as Lgn, are dependent on Lgn for proper localization. The polarization of core PCP proteins is not affected in Lgn mutants; however, Lgn and Gαi are misoriented in a PCP mutant, supporting the role of Lgn as a PCP effector. The kinocilium, previously shown to be dependent on Lgn for robust localization, is essential for proper localization of Lgn, as well as Gαi and aPKC, suggesting that cilium function plays a role in positioning of apical proteins. Taken together, our data provide a mechanism for the loss of hearing found in human patients with GPSM2/LGN variants.     

Evolution of blue-flowered species of genus <Emphasis Type="Italic">Linum</Emphasis> based on high-throughput sequencing of ribosomal RNA genes

Background

The species relationships within the genus Linum have already been studied several times by means of different molecular and phylogenetic approaches. Nevertheless, a number of ambiguities in phylogeny of Linum still remain unresolved. In particular, the species relationships within the sections Stellerolinum and Dasylinum need further clarification. Also, the question of independence of the species of the section Adenolinum still remains unanswered. Moreover, the relationships of L. narbonense and other species of the section Linum require further clarification. Additionally, the origin of tetraploid species of the section Linum (2n?=?30) including the cultivated species L. usitatissimum has not been explored. The present study examines the phylogeny of blue-flowered species of Linum by comparisons of 5S rRNA gene sequences as well as ITS1 and ITS2 sequences of 35S rRNA genes.

Results

High-throughput sequencing has been used for analysis of multicopy rRNA gene families. In addition to the molecular phylogenetic analysis, the number and chromosomal localization of 5S and 35S rDNA sites has been determined by FISH.Our findings confirm that L. stelleroides forms a basal branch from the clade of blue-flowered flaxes which is independent of the branch formed by species of the sect. Dasylinum. The current molecular phylogenetic approaches, the cytogenetic analysis as well as different genomic DNA fingerprinting methods applied previously did not discriminate certain species within the sect. Adenolinum. The allotetraploid cultivated species L. usitatissimum and its wild ancestor L. angustifolium (2n?=?30) could originate either as the result of hybridization of two diploid species (2n?=?16) related to the modern L. gandiflorum and L. decumbens, or hybridization of a diploid species (2n?=?16) and a diploid ancestor of modern L. narbonense (2n?=?14).

Conclusions

High-throughput sequencing of multicopy rRNA gene families allowed us to make several adjustments to the phylogeny of blue-flowered flax species and also reveal intra- and interspecific divergence of the rRNA gene sequences.

     

Associations of rs3740677 within <Emphasis Type="Italic">GAB2</Emphasis> Gene with LOAD in Chinese Han Population

GRB2-associated binding protein 2 (GAB2) has been identified as a crucial factor in Alzheimer’s disease (AD), and ten common variants within GAB2 have been detected to be associated with AD onset risk in genome-wide association studies (GWAS). Here, we first screened a common locus (rs3740677) in 3′ UTR of GAB2 sequence which is targeted by the miRNA-185 and initiatively explored the probable associations of rs3740677 with risk for late-onset AD (LOAD) in a large scale case–control study from Chinese Han populations (992 LOAD patients and 1358 healthy subjects). Eventually, the genotype (P?=?0.024) and allele (P?=?0.008) distribution of rs3740677 showed significant difference between LOAD and control group, and we observed a significant association of T allele in rs3740677 with LOAD risk in multivariate analysis and it decreased the risk for LOAD (dominant: OR?=?0.831, 95 % CI?=?0.702–0.983, P?=?0.031; additive: OR?=?0.855, 95 % CI?=?0.745–0.983, P?=?0.027) adjusted for age, gender, and APOE ε4 status. Our study further confirmed the association of GAB2 and AD. However, the absolute and correct association of rs3740677 with AD still required more investigations in diverse regions and ethnics.     

An expression atlas of miRNAs in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

MicroRNAs (miRNAs) are small non-coding RNAs that regulate a variety of biological processes. MiRNA expression often exhibits spatial and temporal specificity. However, genome-wide miRNA expression patterns in different organs during development of Arabidopsis thaliana have not yet been systemically investigated. In this study, we sequenced small RNA libraries generated from 27 different organ/tissue types, which cover the entire life cycle of Arabidopsis. Analysis of the sequencing data revealed that most miRNAs are ubiquitously expressed, whereas a small set of miRNAs display highly specific expression patterns. In addition, different miRNA members within the same family have distinct spatial and temporal expression patterns. Moreover, we found that some miRNAs are produced from different arms of their hairpin precursors at different developmental stages. This work provides new insights into the regulation of miRNA biogenesis and a rich resource for future investigation of miRNA functions in Arabidopsis.     

Mining small RNA sequencing data: a new approach to identify small nucleolar RNAs in Arabidopsis

Small nucleolar RNAs (snoRNAs) are noncoding RNAs that direct 2′-O-methylation or pseudouridylation on ribosomal RNAs or spliceosomal small nuclear RNAs. These modifications are needed to modulate the activity of ribosomes and spliceosomes. A comprehensive repertoire of snoRNAs is needed to expand the knowledge of these modifications. The sequences corresponding to snoRNAs in 18–26-nt small RNA sequencing data have been rarely explored and remain as a hidden treasure for snoRNA annotation. Here, we showed the enrichment of small RNAs at Arabidopsis snoRNA termini and developed a computational approach to identify snoRNAs on the basis of this characteristic. The approach successfully uncovered the full-length sequences of 144 known Arabidopsis snoRNA genes, including some snoRNAs with improved 5′- or 3′-end annotation. In addition, we identified 27 and 17 candidates for novel box C/D and box H/ACA snoRNAs, respectively. Northern blot analysis and sequencing data from parallel analysis of RNA ends confirmed the expression and the termini of the newly predicted snoRNAs. Our study especially expanded on the current knowledge of box H/ACA snoRNAs and snoRNA species targeting snRNAs. In this study, we demonstrated that the use of small RNA sequencing data can increase the complexity and the accuracy of snoRNA annotation.     

Analysis of antimicrobial and immunomodulatory substances produced by heterofermentative <Emphasis Type="Italic">Lactobacillus reuteri</Emphasis>

Antimicrobial and immunomodulatory potential of various Lactobacillus reuteri strains is closely connected to their metabolite production profile under given cultivation conditions. We determined the in vitro production of antimicrobial substances such as organic acids, ethanol, and reuterin by four strains of L. reuteri (L. reuteri E, L. reuteri KO5, L. reuteri CCM 3625, and L. reuteri ATCC 55730). All studied L. reuteri strains showed the ability to produce lactic acid, acetic acid, and ethanol with concominant consumption of glucose and together with phenyllactic acid—a potent antifungal compound—with concominant consumption of phenylalanine. The reuterin production from glycerol was confirmed for all analyzed lactobacilli strains except L. reuteri CCM 3625. Production of organic acids, ethanol, and reuterin is significantly involved in antimicrobial activity of lactobacilli which was determined using the dual-culture overlay diffusion method against six indicator bacteria and five indicator moulds. In comparison to the referential L. reuteri ATCC 55730, the highest inhibition potential was observed against Escherichia coli CCM 3988 and Pseudomonas aeruginosa CCM 3955. Among analyzed indicators of moulds, the growth of Alternaria alternata CCM F-128 was the most inhibited by all four analyzed L. reuteri strains. Finally, the immunomodulatory potential of analyzed lactobacilli were proven by the determination of the in vitro production of biogenic amines histamine and tyramine. L. reuteri CCM 3625 was able to produce tyramine, and L. reuteri E and L. reuteri KO5 were able to produce histamine under given cultivation conditions.     

Association of <Emphasis Type="Italic">p</Emphasis>53 codon 72 polymorphism and survival of North Indian lung cancer patients treated with platinum-based chemotherapy

p53 helps in maintaining genomic stability by undergoing cellular arrest, DNA repair or cellular apoptosis during DNA damage. So, as to find the association of p53Arg ⁷² Pro towards lung carcinogenesis and overall survival of North Indian lung cancer patients, single nucleotide polymorphic variant (rs1042522) was analyzed. 840 subjects including 420 cases and 420 controls were recruited and genotyped using PCR-RFLP technique for p53Arg ⁷² Pro polymorphic site. Association was analyzed using adjusted odds ratio along with its confidence intervals (95?% CI) and p value predicted from logistic regression whereas overall survival for lung cancer patients was obtained using Kaplan–Meir and Cox regression model for different parameters to obtain hazard ratio and survival time with statistical significance (log-rank p value). None of the variant genotypes for p53Arg ⁷² Pro showed any association towards lung cancer risk or any specific histological subtype. Lung cancer subjects with Pro/Pro genotype had better median survival time as compared to Arg/Pro genotype (10 months; HR?=?0.65; 95?% CI?=?0.45–0.95; p?=?0.03). Furthermore, female lung cancer patients with Arg/Pro (HR?=?0.08; 95?% CI?=?0.02–0.34; p?=?0.0005) and Pro/Pro (HR?=?0.21; 95?% CI?=?0.06–0.67; p?=?0.008) genotypes showed a better overall survival and hence a better prognosis as compared to males. Our data also reveals that lung cancer patients with ECOG scores between 0 and 1 and carrying the Pro/Pro had better chances of survival. p53 codon 72 polymorphism could play a role as a prognostic marker in lung cancer patients.     

Development of user-friendly markers for disease resistance to black root rot of tobacco through genotyping by sequencing

Black root rot (BRR), a disease caused by the hemibiotrophic fungus Thielaviopsis basicola, seriously compromises yield and leaf quality in tobacco (Nicotiana tabacum). Full resistance to black root rot, conferred by the resistance to BRR 1 (RBRR1) locus from Nicotiana debneyi Domin, was transferred to a burley tobacco cultivar through interspecific hybridization. Some undesirable traits potentially caused by linkage drag restrict wider application of RBRR1 in flue-cured tobacco. Therefore, user-friendly molecular markers are needed to assist selection for resistance to black root rot and to break the unfavorable linkage. Genotyping by sequencing (GBS) is a rapid and robust approach for reduced representation sequencing of multiplexed genomic DNA samples that combines genome-wide molecular marker discovery with genotyping. In the present study, we used GBS to identify single-nucleotide polymorphisms (SNPs) linked to the RBRR1 locus, and PCR-based assays for detection of these SNPs were also developed. Sequence analysis of the SNP markers suggested that RBRR1 is located on chromosome 17, providing a basis for map-based cloning of this valuable gene. Co-dominant CAPS markers that co-segregate with the disease-resistant phenotype offer user-friendly tools for tobacco breeding and variety improvement. Furthermore, tested with diverse N. tabacum germplasm, SS192650 displayed 100% selection accuracy for resistance to BRR, suggesting that this marker can be used in diverse tobacco populations.     

In silico identification and computational characterization of endogenous small interfering RNAs from diverse grapevine tissues and stages

Small interfering RNAs (siRNAs) are effectors of regulatory pathways underlying plant development, metabolism, and stress- and nutrient-signaling regulatory networks. The endogenous siRNAs are generally not conserved between plants; consequently, it is necessary and important to identify and characterize siRNAs from various plants. To address the nature and functions of siRNAs, and understand the biological roles of the huge siRNA population in grapevine (Vitis vinifera L.). The high-throughput sequencing technology was used to identify a large set of putative endogenous siRNAs from six grapevine tissues/organs. Subsequently, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed to classify the target genes of siRNA. In total, 520,519 candidate siRNAs were identified and their expression profiles exhibited typical temporal characters during grapevine development. In addition, we identified two grapevine trans-acting siRNA (TAS) gene homologs (VvTAS3 and VvTAS4) and the derived trans-acting siRNAs (tasiRNAs) that could target grapevine auxin response factor (ARF) and myeloblastosis (MYB) genes. Furthermore, the GO and KEGG analysis of target genes showed that most of them covered a broad range of functional categories, especially involving in disease-resistance process. The large-scale and completely genome-wide level identification and characterization of grapevine endogenous siRNAs from the diverse tissues by high throughput technology revealed the nature and functions of siRNAs in grapevine.     

Role of the box C/D motif in localization of small nucleolar RNAs to coiled bodies and nucleoli.

Small nucleolar RNAs (snoRNAs) are a large family of eukaryotic RNAs that function within the nucleolus in the biogenesis of ribosomes. One major class of snoRNAs is the box C/D snoRNAs named for their conserved box C and box D sequence elements. We have investigated the involvement of cis-acting sequences and intranuclear structures in the localization of box C/D snoRNAs to the nucleolus by assaying the intranuclear distribution of fluorescently labeled U3, U8, and U14 snoRNAs injected into Xenopus oocyte nuclei. Analysis of an extensive panel of U3 RNA variants showed that the box C/D motif, comprised of box C', box D, and the 3' terminal stem of U3, is necessary and sufficient for the nucleolar localization of U3 snoRNA. Disruption of the elements of the box C/D motif of U8 and U14 snoRNAs also prevented nucleolar localization, indicating that all box C/D snoRNAs use a common nucleolar-targeting mechanism. Finally, we found that wild-type box C/D snoRNAs transiently associate with coiled bodies before they localize to nucleoli and that variant RNAs that lack an intact box C/D motif are detained within coiled bodies. These results suggest that coiled bodies play a role in the biogenesis and/or intranuclear transport of box C/D snoRNAs.     

Genome-wide association study of body weight in Wenshang Barred chicken based on the SLAF-seq technology

Chicken body weight (BW) is an economically important trait, and many studies have been conducted on genetic selection for BW. However, previous studies have detected functional chromosome mutations or regions using gene chips. The present study used the specific-locus amplified fragment sequencing (SLAF-seq) technology to perform a genome-wide association study (GWAS) on purebred Wengshang Barred chicken. A total of 1,286,715 single-nucleotide polymorphisms (SNPs) were detected, and 175,211 SNPs were selected as candidate SNPs for genome-wide association analysis using TASSEL general linear models. Six SNP markers reached genome-wide significance. Of these, rs732048524, rs735522839, rs738991545, and rs15837818 were significantly associated with body weight at 28 days (BW28), while rs314086457 and rs315694878 were significantly associated with BW120. These SNPs are close to seven genes (PRSS23, ME3, FAM181B, NABP1, SDPR, TSSK6L2, and RBBP8). Moreover, 24 BW-associated SNPs reached “suggestive” genome-wide significance. Of these, 6, 13, 1, and 4 SNPs were associated with BW28, BW56, BW80, and BW120, respectively. These results would enrich the studies on BW and promote the use of Chinese chicken, especially the Wenshang Barred chicken.