• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Long non-coding RNA Databases in Cardiovascular Research

    2016-11-17 08:41:14FrankRuhleMonikaStoll
    Genomics,Proteomics & Bioinformatics 2016年4期

    Frank Ru¨hleMonika Stoll*b

    1Institute of Human Genetics,Genetic Epidemiology,University of Muenster,48149 Muenster,Germany

    2Cardiovascular Research Institute Maastricht(CARIM),Genetic Epidemiology and Statistical Genetics,Maastricht Center for Systems Biology(MaCSBio),Maastricht University,6211 LK Maastricht,The Netherlands

    RESOURCE REVIEW

    Long non-coding RNA Databases in Cardiovascular Research

    Frank Ru¨hle1,a,Monika Stoll1,2,*,b

    1Institute of Human Genetics,Genetic Epidemiology,University of Muenster,48149 Muenster,Germany

    2Cardiovascular Research Institute Maastricht(CARIM),Genetic Epidemiology and Statistical Genetics,Maastricht Center for Systems Biology(MaCSBio),Maastricht University,6211 LK Maastricht,The Netherlands

    Received 18 January 2016;revised 16 March 2016;accepted 17 March 2016 Available online 2 April 2016

    Handled by Andreas Keller

    Database;

    Non-coding;

    lncRNA;

    Gene regulation;

    Cardiovascular disease;ANRIL

    With the rising interest in the regulatory functions of long non-coding RNAs(lncRNAs)in complex human diseases such as cardiovascular diseases,there is an increasing need in public databases offering comprehensive and integrative data for all aspects of these versatile molecules. Recently,a variety of public data repositories that specialized in lncRNAs have been developed,which make use of huge high-throughput data particularly from next-generation sequencing(NGS)approaches.Here,we provide an overview of current lncRNA databases covering basic and functional annotation,lncRNA expression and regulation,interactions with other biomolecules,and genomic variants influencing the structure and function of lncRNAs.The prominent lncRNA antisense noncoding RNA in the INK4 locus(ANRIL),which has been unequivocally associated with coronary artery disease through genome-wide association studies(GWAS),serves as an example to demonstrate the features of each individual database.

    Introduction

    Although substantial genetic heritability is estimated for complex cardiovascular diseases,e.g.,40%in coronary artery disease(CAD)[1],and extraordinary efforts have been made in genome-wide association studies(GWAS)and meta-analyses to identify genetic variants leading to CAD,only a small fraction of genetic variance of CAD of~10%can be explained by genetic variants in protein-coding genes[2].Additionally,the high proportion of GWAS associations in non-coding genome regions contradicts the simple view of potentiallydeleterious protein mutations and indicates a complex regulatory network driven by non-coding RNAs(ncRNAs)[3,4]. Since only 1%of the mammalian genome is translated into proteins,but approximately 85%of the genome is transcribed into RNA,ncRNAs potentially represent an additional layer of epigenetic regulation.Especially long ncRNAs(lncRNAs,RNA>200 nucleotides in length)provide a wide range of regulatory functions including interactions with DNA,RNAs,and proteins[5,6].

    For instance,the lncRNA X-inactive-specific transcript(XIST)directly binds to the polycomb repressive complex 2(PRC2)and thereby downregulates the entire chromosome during X-chromosome inactivation[7].Other lncRNAs influence gene activity by RNA-directed chromatin remodeling[8],RNA-directed DNA methylation[9]or as activator or repressor molecules for transcription factors(TFs)[10,11]. By recruiting splicing factors or by masking splice junctions of mRNAs,lncRNAs can influence alternative splicing of coding genes[12].Various lncRNA interactions with microRNAs(miRNAs)impact mRNA stability by masking miRNA-responsive elements or by competing for miRNA binding in competing endogenous RNA(ceRNA)networks[13,14]. Additionally,discrimination between coding and non-coding genes is sometimes ambiguous,because functional lncRNA transcripts containing open reading frames may also be translated to(small)proteins[15].

    Dysregulated expression or function of lncRNAs has been recognized to contribute to heart development and complex cardiovascular diseases[6].For instance,transcript levels of the antisense noncoding RNA in the INK4 locus(ANRIL,alias CDKN2B-AS1)lncRNA,which is encoded on chromosome 9p21 at the strongest genetic susceptibility locus for CAD,are directly correlated with the severity of atherosclerosis[16].The locus at chromosome 5q31 carrying the noncoding steroid receptor RNA activator(SRA1)as well as human leukocyte antigen(HLA)complex group 22(HCG22)at chromosome 6p21 have been significantly associated with dilated cardiomyopathy(DCM)[17,18].The myocardial infarction(MI)-associated transcript(MIAT)encoded on chromosome 22q12 is implicated to play a role in MI[19].

    To discover potentially harmful lncRNA functions,it is important to understand the complex interaction networks of these molecules.In general,expression of lncRNAs is more specific for cell type and developmental stage than that of protein-coding genes[20].Functional prediction of lncRNAs is more difficult than,e.g.,that for smaller miRNAs,because function of a lncRNA is not solely determined by its nucleotide sequence,but by the resulting secondary structure enabling it to interact with other biomolecules[21].This is supported by the fact that lncRNA sequences are less conserved than miRNAs or protein-coding genes except for their promoter regions[22].Genomic variants in lncRNA sequences may induce abnormal expression and function of their harboring lncRNAs,e.g.,by gaining or losing binding sites for interaction partners or by altering the secondary structure even at distant positions of the RNA molecule,possibly explaining part of genetic susceptibility to certain diseases[23].Many of the aforementioned disease-associated lncRNAs like ANRIL,MIAT,and HCG22 have gene variants,whose structural impact is not yet understood.Generally,there is a large gap between the number of identified lncRNAs and their known functional impact.

    Therefore there is a need for comprehensive lncRNA databases to utilize the huge experimental datasets from current high-throughput technologies joined by massively parallel sequencing such as RNA-Seq,chromatin immunoprecipitation(ChIP)-Seq,RNA immunoprecipitation(RIP)-Seq,crosslinking immunoprecipitation(CLIP)-Seq or chromatin isolation by RNA purification(ChIRP)-Seq[24].In addition to the main genomics data portals from NCBI,EMBL and UCSC,which also provide data on non-coding genes,several specialized databases have been developed that collect and integrate data in the context of lncRNAs[25,26].Alldatabases discussed here are accessible via a web-based interface and have been published in peer-reviewed journals(Table 1 and Table 2).Apart from these,there are further data repositories with down loadable data files like the Human lincRNA Catalog from Broad Institute[20].While some databases have been performing well for several years,many specialized databases have been developed in very recent time,highlighting the strong momentum of this research field.

    In the following,we will give an overview of selected databases for different kinds of lncRNA-related information(Figure 1).The suggested analysis outline is exemplified by the CAD-related human lncRNA ANRIL.All specifications of database contents and query results refer to status of March,2016(Table 1 and Table 2).To start a query for a lncRNA of interest,basic information about lncRNA type,chromosomal location,nucleotide sequence,expression profiles,and functional annotation may be retrieved at NONCODE[27-31]and lncRNAdb[32].Roughly,the classification of lncRNA types is based on their genomic context concerning sense,antisense,bidirectional,intergenic(lincRNAs),or intronic lncRNAs[33].Known biological functions such as gene ontology(GO)annotation and disease associations are documented in LncRNAD is ease[34].In addition to providing lncRNA expression profiles,lncRNAtor[35]offers coexpression data for protein-coding genes to identify potential functional connections between coding and non-noncoding transcripts.To gain insights in regulation of lncRNA expression,ChIPBase[36]contains information on TFs that regulate the expression of non-coding genes.In the next step,the interactions of lncRNAs with other biomolecules may be examined by using star Base[37]and NPInter[38].These databases provide experimentally-validated data on interactions with proteins,DNA,and other RNA types,especially miRNAs. Finally,genomic variations within the lncRNA gene sequence can be analyzed to explore their potentialfunctional influence on the lncRNA transcript using lncRNASNP[39].

    NONCODE 2016

    NONCODE was first published in 2005 as an integrated knowledge database of ncRNAs[27]and has repeatedly been updated since then[28-31].Its latest version NONCODE 2016 offers data for 16 species including 167,150 human lncRNAs[31].In addition to lncRNA class,chromosomal location,sequence,Coding-Non-Coding Index(CNCI)for protein coding prediction and expression profiles,the database included conservation annotation and disease association as new features in its latest version.The collected data are curated from published literature and include input from other public databases such as Ensembl[40],RefSeq[41],lncRNAdb[32],and GENCODE[42].The database established a lncRNA nomenclature consisting of‘NON”,a three character code that specifies the species,‘T”or‘G”for transcript or gene,respectively,followed by six sequential numbers and a version number where applicable.For ANRIL(NONHSAG051899),we find 22 transcript isoforms of type‘a(chǎn)ntisense’on chromosome 9,which are mostly expressed in lung,lymph nodes,prostate,skeletal muscle,and white blood cells.However,although this molecule has been linked in the literature tocardiovascular diseases and several other pathologies,we don’t find any disease association for ANRIL in NONCODE.

    Table1 Overview of current lncRNA databases

    lncRNAdb v2.0

    Unlike NONCODE,lncRNAdb[32,43]contains only functionally-annotated entries manually curated from referenced literature.lncRNAs that have been associated with diseases but have not been further characterized by knockdown or over expression experiments,are not included in the repository.To date,the database contains 295 function ally annotated lncRNAs covering 71 species,including 183 lncRNAs annotated in human.The database gives information on lncRNA type,sequence,chromosomal and subcellular localization,characteristics and functional annotation complemented by literature references,evolutionary conservation,interactions with other biomolecules,as well as expression profiles based on the Illumina body map[44].In lncRNAdb,ANRIL is implicated in a range of complex diseases including

    cancer,T2D,and coronary heart disease.Its expression in tissues and cell types that are affected by atherosclerosis,such as peripheral blood mononuclear cells,whole blood,and athero sclerotic plaque tissue,is directly correlated with the severity of atherosclerosis.Functional interaction of ANRIL is described for chromobox 7(CBX7),a component of the PRC1[8].According to the deposited expression profiles,ANRIL is mainly expressed in colon tissue.

    Table2 Web links and data content of the presented lnc RNA databases

    Figure1 Types of information curated in lncRNA databases

    LncRNADisease

    LncRNADisease[34]collects experimentally-validated disease associations of human lncRNAs extracted from the literature. By now,the database contains more than 1000 lncRNA-disease annotations including 321 lncRNAs and 221 diseases from about 500 publications.LncRNADisease also curates 475 entries of validated lncRNA interactions with other biomolecules including protein,RNA,and DNA.In addition,a computational method has been developed to predict new potential disease associations for a given lncRNA based on its genomic context.The database can be queried for either lncRNAs or diseases.Since ANRIL belongs to the wellannotated lncRNAs contributing to disease,we find a total of 134 lncRNA-disease associations described in 65 publications and affecting 37 disease phenotypes including CAD,MI,T2D,and several cancer types.Additionally,25 interaction entries of ANRIL with 9 different biomolecules are annotated,including co-expression and regulatory interactions with its protein-coding counterparts CDKN2A and CDKN2B which encode cyclin-dependent kinase inhibitor 2A and B,and protein binding interaction with PRC1 and PRC2.

    lncRNAtor

    Expression data from 243 RNA-Seq experiments comprising 5237 samples of various tissues and development alstages have been collected from the public databases,including Gene Expression Omnibus(NCBI GEO)[45],ENCODE[46],mod-ENCODE[47],and The Cancer Genome Atlas(TCGA)[48],and are updated on an annual basis[35].The lncRNA compendium was taken from Ensembl[40],HUGO Gene Nomenclature Committee(HGNC)[49],Mouse Genome Database(MGD),[50]and lncRNAdb[32],and comprises a total of 21,575 lncRNA genes from human,mouse,zebrafish,fruit fly,worm,and yeast.In addition to visualizing tissue-specific expression profiles of lncRNAs,expression data can be searched for co-expression of mRNAs to identify putative lncRNA-mRNA pairs.Functional investigation of lncRNAs is complemented by CLIP-Seq and RIP-Seq data included from public data repositories to identify potential proteinlncRNA interactions.As most of the included human datasets are cancer-related,we find ANRIL to be upregulated in several cancer types compared to normaltissue,namely in kidney-and liver-related carcinoma.Co-expression analysis in a dataset of kidney renal clear cell carcinoma reveals ARF5(encoding ADP ribosylation factor 5)as highly-correlated coding gene,which possibly interacts with ANRIL in trans.Protein interactions are displayed for 12 proteins,including argonaute 2(AGO2,RNA-induced silencing complex catalytic component)and DiGeorge syndrome critical region 8(DGCR8),indicating potential involvement in miRNA regulation.

    ChIPBase

    ChIPBase[36]aims at analysis of the trans criptional regulation of lncRNAs and miRNAs.It contains TF-lncRNA and TF-miRNA regulatory relationships identified by data coming from 543 ChIP-Seq experiments for 252 different TFs retrieved from the respective research articles and the NCBI GEO[45],ENCODE[46],and modENCODE[47]databases.The collection comprises diverse tissues and cell lines from human,mouse,dog,chicken,fruit fly,and nematodes(TF-lncRNA relationships are not available for dog,chicken,and nematode).Additionally the database is complemented by human expression profiles from 22 tissues.ChIPBase can be queried by lncRNAs,miRNAs or TFs,and the results are visualized by an integrative genome browser.For ANRIL,ChIPBase displays experimentally-supported binding sites of 18 different TFs in human,including androgen receptor(AR),vets avian erythroblastosis virus E26 oncogene homolog(ERG),and signal transducer and activator of transcription 1(STAT1).

    NPInter v3.0

    NPInter v3.0[38,51]provides experimentally-verified functional interactions between ncRNAs and other biomolecules such as proteins,RNAs and genomic DNA.Interaction data for 23 different species(mainly human and mouse)are collected from literature datasets and related databases such as lncRNADisease[34].ncRNAs are screened against NONCODE[31],which serves as ncRNA reference database.The majority of included data stems from systematic identification of protein-binding sites by CLIP-Seq experiments,while other interactions such as ncRNA-RNA and TF-ncRNA are obtained mainly from interaction studies on individual ncRNAs.NPInter classifies allinteractions as‘binding’,‘regulatory’,or‘co-expression’.Every interaction entry includes a description of the kind of interaction and the interacting partner,complemented by the source of experimental data and the corresponding PubMed ID.Additionally,computational tools have been added to its latest version to predict further RNARNA and RNA-protein interactions.For ANRIL(query for NONCODE ID NONHSAG051899),NPInter displays 73 interactions,including RNA-DNA and RNA-protein binding interactions to its protein-coding counterparts CDKN2A and CDKN2B,regulatory interaction with miRNA hsa-miR-106a,and binding to the TF STAT1.

    starBase 2.0

    starBase 2.0[37,52]collected 111 CLIP-Seq data sets from various tissues and cell lines generated by 40 independent studies from the NCBIGEO[45]to explore protein-RNA and various RNA-RNA interactions as well as ceRNA regulatory networks involving miRNAs,lncRNAs and mRNAs.miRNA data and gene annotations were retrieved from miRBase[53],GENCODE[54],Ensembl[40],and RefSeq[41],respectively. miRNA target sites on lncRNAs are predicted by miRanda[55]and subsequently filtered for CLIP-supported interactions. For ANRIL(query for CDKN2B-AS1 because the gene symbol ANRIL is not found by starBase),21 human miRNA-lncRNA interactions are annotated in the database. Interestingly,these do not include hsa-miR-106a identified by NPInter mentioned above.In addition,expression profiles are given for miRNAs and lncRNAs if available.When searching for ceRNA networks involving ANRIL and a minimum of 5 common miRNAs,we find TMEM41A coding for trans membrane protein 41A to be part of the network.

    lncRNASNP

    lncRNASNP aims at the influence of genetic variants on the expression and function of the encoded lncRNAs.This influence may arise from gain or loss of binding sites for miRNAs or induction of conformational changes within the secondary structure of a lncRNA.Therefore,lncRNASNP collected SNP data and lncRNAs from dbSNP[56],LNCipedia[57],and NONCODE[31],respectively.Changes in secondary structure are predicted by RNAfold[58]based on the minimal free energy of the alternative transcript sequence.miRNA sequences were downloaded from mirBase[53]and are used to predict target sites on lncRNAs using the TargetScan[59]and miRanda[55]algorithms.Furthermore,experiment ally supported lncRNA-miRNA interactions from star Base[37]and disease associations from the National Human Genome Research Institute(NHGRI)GWAS Catalog[60]are embedded in the database.lncRNASNP is divided into human and mouse sub-databases and can be queried for SNPs,lncRNAs,miRNAs,or genomic regions.It returns 17 transcripts for ANRIL(query for CDKN2B-AS1 because the gene symbol ANRIL is not found by lncRNASNP).Transcript CDKN2BAS1-001 contains 20 SNPs and 90 predicted miRNA-binding sites.Two binding sites are gained due to alternative SNP alleles while 9 binding sites are lost.The secondary structures of wild type and variant sequence can be visualized for each SNP of interest.

    Other resources

    In addition to the data repositories presented above,several other public data resources for lncRNA research are listed in Table 1 and Table 2.For instance,LNCipedia[57,61]summarizes 111,685 human lncRNA transcripts from Ensembl[40],lncRNAdb[32],NONCODE[31],RefSeq[41],the Human lincRNA Catalog[20],and two further datasets published by Hangauer et al and Nielsen et al[62,63].It offers transcript and structure information as well as computational scores for protein-coding potential and miRNA-binding sites. lncRNome[64]also provides a range of general annotations for sequence,structure,function,variation,and epigenetic modifications for more than 17,000 human lncRNAs derived from public databases.For protein-lncRNA interactions,the database included published photo activatable-ribonucleosideenhanced CLIP(PAR-CLIP)experiments and computational prediction methods.More specialized databases exist for evolutionary conservation(PhyloNONCODE[65])and functional annotation based on ceRNA interaction networks(Linc2GO[66]).Expression profiles of non-coding and coding genes arefurther available from RNA-Seq experiments(lncRNAMap[67])and microarray platforms(NRED[68]).Co-LncRNA identifies co-expressed coding genes from RNA-Seq data,which are then functionally annotated.Potential influence of lncRNAs on target gene expression may be identified with LncRNA2Target[69],which contains manually-curated differential expression data from 217 lncRNA knockdown or over expression experiments for human and mouse.Further interaction data between lncRNAs and other biomolecules can be found in LncReg[70],DIANA-LncBase[71],or lncRNAMap[67].Many of these databases provide information for human and murines only,but there are also databases specialized in other model organisms such as zflncRNApedia[72]for zebrafish(Danio rerio)or PLncDB[73]for Arabidopsis thaliana.

    Concluding remarks

    The growing number of interconnected lncRNA databases reflects the immense research interest in lncRNAs,which is increasingly gaining momentum in the quest to understanding the(dys)function of biomolecular networks potentially contributing to complex human diseases[74].Current high-throughput technologies joined with massive parallel sequencing generate data for non-coding transcripts at an unprecedented scale.To date,there is still a strong disconnection between the large number of identified transcripts and the small amount of lncRNA functional data,which is illustrated best by two of the most cited lncRNA databases,NONCODE and lncRNAdb.While NONCODE contains as much as 167,150 known human lncRNA transcripts,lncRNAdb is dedicated to functionally-characterized lncRNAs,restricting its content to 183 human lncRNAs.However,even for well characterized lncRNAs,such as ANRIL,further investigation is warranted.Despite the wealth of information from public databases,the exact mechanisms of ANRIL functionality remain enigmatic.Another drawback are occasional discrepancies across databases for similar queries,which force researchers to use and compare several databases[25].When choosing a database,researchers should also assure that the database of interest is curated and regularly updated as novel information becomes available.For instance,the Functional lncRNADatabase[75]was not considered for this review since it has last been updated in March 2012.Nevertheless,current databases offer valuable resources for integration and interpretation of various kinds of experimental lncRNA data.This is essentialfor understanding the function and relevance of these versatile molecules and may pave the way to new translational applications in cardiovascular research.

    Competing interests

    The authors declare that there are no conflicts of interests.

    Acknowledgments

    We gratefully thank our bioinformatics colleagues at the Institute of Human Genetics,University of Muenster,Germany and Professor Dr.Leon de Windt,Department of Cardiology,Maastricht University,the Netherlands for helpful discussions.

    [1]Marenberg ME,Risch N,Berkman LF,F(xiàn)loderus B,de Faire U. Genetic susceptibility to death from coronary heart disease in a study of twins.N Engl J Med 1994;330:1041-6.

    [2]Schunkert H,Konig IR,Kathiresan S,Reilly MP,Assimes TL,Holm H,et al.Large-scale association analysis identifies 13 new susceptibility loci for coronary artery disease.Nat Genet 2011;43:333-8.

    [3]Dunham I,Kundaje A,Aldred SF,Collins PJ,Davis CA,Doyle F,et al.An integrated encyclopedia of DNA elements in the human genome.Nature 2012;489:57-74.

    [4]Manolio TA,Collins FS,Cox NJ,Goldstein DB,Hindorff LA,Hunter DJ,et al.Finding the missing heritability of complex diseases.Nature 2009;461:747-53.

    [5]Lee JT.Epigenetic regulation by long noncoding RNAs.Science 2012;338:1435-9.

    [6]Schonrock N,Harvey RP,Mattick JS.Long noncoding RNAs in cardiac development and pathophysiology.Circ Res 2012;111:1349-62.

    [7]Wutz A,Rasmussen TP,Jaenisch R.Chromosomal silencing and localization are mediated by different domains of Xist RNA.Nat Genet 2002;30:167-74.

    [8]Yap KL,Li S,Munoz-Cabello AM,Raguz S,Zeng L,Mujtaba S,et al.Molecular interplay of the noncoding RNA ANRIL and methylated histone H3 lysine 27 by polycomb CBX7 in transcriptional silencing of INK4a.Mol Cell 2010;38:662-74.

    [9]Imamura T,Yamamoto S,Ohgane J,HattoriN,Tanaka S,Shiota K.Non-coding RNA directed DNAdemethylation of Sphk1 CpG island.Biochem Biophys Res Commun 2004;322:593-600.

    [10]Feng J,Bi C,Clark BS,Mady R,Shah P,Kohtz JD.The Evf-2 noncoding RNA is transcribed from the Dlx-5/6 ultraconserved region and functions as a Dlx-2 transcriptional coactivator.Genes Dev 2006;20:1470-84.

    [11]Martianov I,Ramadass A,Serra Barros A,Chow N,Akoulitchev A.Repression of the human dihydrofolate reductase gene by a non-coding interfering transcript.Nature 2007;445:666-70.

    [12]Anko ML,Neugebauer KM.Long noncoding RNAs add another layer to pre-mRNA splicing regulation.Mol Cell 2010;39:833-4.

    [13]Mercer TR,Dinger ME,Mattick JS.Long non-coding RNAs: insights into functions.Nat Rev Genet 2009;10:155-9.

    [14]Tay Y,Kats L,Salmena L,Weiss D,Tan SM,Ala U,et al. Coding-independent regulation of the tumor suppressor PTENby competing endogenous mRNAs.Cell 2011;147:344-57.

    [15]Nelson BR,Makarewich CA,Anderson DM,Winders BR,Troupes CD,Wu F,et al.A peptide encoded by a transcript annotated as long noncoding RNA enhances SERCA activity in muscle.Science 2016;351:271-5.

    [16]Holdt LM,Beutner F,Scholz M,Gielen S,Gabel G,Bergert H,et al.ANRIL expression is associated with atherosclerosis risk at chromosome 9p21.Arterioscler Thromb Vasc Biol2010;30:620-7.

    [17]Friedrichs F,Zugck C,Rauch GJ,Ivandic B,Weichenhan D,Muller-Bardorff M,et al.HBEGF,SRA1,and IK:Three cosegregating genes as determinants of cardiomyopathy.Genome Res 2009;19:395-403.

    [18]Meder B,Ruhle F,Weis T,Homuth G,Keller A,F(xiàn)ranke J,et al. A genome-wide association study identifies 6p21 as novel risk locus for dilated cardiomyopathy.Eur Heart J 2014;35:1069-77.

    [19]Ishii N,Ozaki K,Sato H,Mizuno H,Saito S,Takahashi A,et al. Identification of a novel non-coding RNA,MIAT,that confers risk of myocardial infarction.J Hum Genet 2006;51:1087-99.

    [20]Cabili MN,Trapnell C,Goff L,Koziol M,Tazon-Vega B,Regev A,et al.Integrative annotation of human large intergenic noncoding RNAs reveals globalproperties and specific subclasses. Genes Dev 2011;25:1915-27.

    [21]Johnsson P,Lipovich L,Grander D,Morris KV.Evolutionary conservation of long non-coding RNAs;sequence,structure,function.Biochim Biophys Acta 2014;1840:1063-71.

    [22]Necsulea A,Kaessmann H.Evolutionary dynamics of coding and non-coding transcriptomes.Nat Rev Genet 2014;15: 734-48.

    [23]Solem AC,Halvorsen M,Ramos SB,Laederach A.The potential of the riboSNitch in personalized medicine.Wiley Interdiscip Rev RNA 2015;6:517-32.

    [24]Chu C,Spitale RC,Chang HY.Technologies to probe functions and mechanisms of long noncoding RNAs.Nat Struct Mol Biol 2015;22:29-35.

    [25]Fritah S,Niclou SP,Azuaje F.Databases for lncRNAs: a comparative evaluation of emerging tools.RNA 2014;20: 1655-65.

    [26]Jalali S,Kapoor S,Sivadas A,Bhartiya D,Scaria V.Computationalapproaches towards understanding human long non-coding RNA biology.Bioinformatics 2015;31:2241-51.

    [27]Liu C,Bai B,Skogerbo G,Cai L,Deng W,Zhang Y,et al. NONCODE:an integrated knowledge database of non-coding RNAs.Nucleic Acids Res 2005;33:D112-5.

    [28]He S,Liu C,Skogerbo G,Zhao H,Wang J,Liu T,et al. NONCODE v2.0:decoding the non-coding.Nucleic Acids Res 2008;36:D170-2.

    [29]Bu D,Yu K,Sun S,Xie C,Skogerbo G,Miao R,et al. NONCODE v3.0:integrative annotation of long noncoding RNAs.Nucleic Acids Res 2012;40:D210-5.

    [30]Xie C,Yuan J,Li H,Li M,Zhao G,Bu D,et al.NONCODEv4: exploring the world of long non-coding RNA genes.Nucleic Acids Res 2014;42:D98-103.

    [31]Zhao Y,Li H,F(xiàn)ang S,Kang Y,Wu W,Hao Y,etal.NONCODE 2016:an informative and valuable data source of long non-coding RNAs.Nucleic Acids Res 2016;44:D203-8.

    [32]Quek XC,Thomson DW,Maag JL,Bartonicek N,Signal B,Clark MB,et al.lncRNAdb v2.0:expanding the reference database for functional long noncoding RNAs.Nucleic Acids Res 2015;43:D168-73.

    [33]Ma L,Bajic VB,Zhang Z.On the classification of long noncoding RNAs.RNA Biol 2013;10:925-33.

    [34]Chen G,Wang Z,Wang D,Qiu C,Liu M,Chen X,et al. LncRNADisease:a database for long-non-coding RNA-associated diseases.Nucleic Acids Res 2013;41:D983-6.

    [35]Park C,Yu N,Choi I,Kim W,Lee S.lncRNAtor:a comprehensive resource for functional investigation of long non-coding RNAs.Bioinformatics 2014;30:2480-5.

    [36]Yang JH,Li JH,Jiang S,Zhou H,Qu LH.ChIPBase:a database for decoding the transcriptional regulation of long non-coding RNA and microRNA genes from ChIP-Seq data.Nucleic Acids Res 2013;41:D177-87.

    [37]Li JH,Liu S,Zhou H,Qu LH,Yang JH.starBase v2.0:decoding miRNA-ceRNA,miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data.Nucleic Acids Res 2014;42:D92-7.

    [38]Wu T,Wang J,Liu C,Zhang Y,Shi B,Zhu X,et al.NPInter:the noncoding RNAs and protein related biomacromolecules interaction database.Nucleic Acids Res 2006;34:D150-2.

    [39]Gong J,Liu W,Zhang J,Miao X,Guo AY.lncRNASNP:a database of SNPs in lncRNAs and their potential functions in human and mouse.Nucleic Acids Res 2015;43:D181-6.

    [40]Cunningham F,Amode MR,Barrell D,Beal K,Billis K,Brent S,et al.Ensembl 2015.Nucleic Acids Res 2015;43:D662-9.

    [41]Pruitt KD,Brown GR,Hiatt SM,Thibaud-Nissen F,Astashyn A,Ermolaeva O,et al.RefSeq:an update on mammalian reference sequences.Nucleic Acids Res 2014;42:D756-63.

    [42]Derrien T,Johnson R,Bussotti G,Tanzer A,Djebali S,Tilgner H,et al.The GENCODE v7 catalog of human long noncoding RNAs:analysis of their gene structure,evolution,and expression. Genome Res 2012;22:1775-89.

    [43]Amaral PP,Clark MB,Gascoigne DK,Dinger ME,Mattick JS. lncRNAdb:a reference database for long noncoding RNAs. Nucleic Acids Res 2011;39:D146-51.

    [44]Petryszak R,Burdett T,F(xiàn)iorelliB,F(xiàn)onseca NA,Gonzalez-Porta M,Hastings E,et al.Expression Atlas update-a database of gene and transcript expression from microarray-and sequencing-based functional genomics experiments.Nucleic Acids Res 2014;42:D926-32.

    [45]Barrett T,Wilhite SE,Ledoux P,Evangelista C,Kim IF,Tomashevsky M,et al.NCBI GEO:archive for functional genomics data sets-update.Nucleic Acids Res 2013;41:D991-5.

    [46]Consortium EP.A user’s guide to the encyclopedia of DNA elements(ENCODE).PLoS Biol 2011;9:e1001046.

    [47]Muers M.Functional genomics:the modENCODE guide to the genome.Nat Rev Genet 2011;12:80.

    [48]Cancer Genome Atlas Research N,Weinstein JN,Collisson EA,Mills GB,Shaw KR,et al.The Cancer Genome Atlas Pan-Cancer analysis project.Nat Genet 2013;45:1113-20.

    [49]Gray KA,Yates B,Seal RL,Wright MW,Bruford EA. Genenames.org:the HGNC resources in 2015.Nucleic Acids Res 2015;43:D1079-85.

    [50]Eppig JT,Blake JA,Bult CJ,Kadin JA,Richardson JE,Mouse Genome Database G,et al.The Mouse Genome Database(MGD):facilitating mouse as a model for human biology and disease.Nucleic Acids Res 2015;43:D726-36.

    [51]Yuan J,Wu W,Xie C,Zhao G,Zhao Y,Chen R.NPInter v2.0: an updated database of ncRNA interactions.Nucleic Acids Res 2014;42:D104-8.

    [52]Yang JH,Li JH,Shao P,Zhou H,Chen YQ,Qu LH.StarBase:a database for exploring microRNA-mRNA interaction maps from Argonaute CLIP-Seq and Degradome-Seq data.Nucleic Acids Res 2011;39:D202-9.

    [53]Kozomara A,Griffiths-Jones S.MiRBase:integrating microRNA annotation and deep-sequencing data.Nucleic Acids Res 2011;39: D152-7.

    [54]Harrow J,F(xiàn)rankish A,Gonzalez JM,Tapanari E,Diekhans M,Kokocinski F,et al.GENCODE:the reference human genome annotation for The ENCODE Project.Genome Res 2012;22:1760-74.

    [55]Betel D,Wilson M,Gabow A,Marks DS,Sander C.The microRNA.org resource:targets and expression.Nucleic Acids Res 2008;36:D149-53.

    [56]Sherry ST,Ward MH,Kholodov M,Baker J,Phan L,Smigielski EM,et al.DbSNP:the NCBI database of genetic variation. Nucleic Acids Res 2001;29:308-11.

    [57]Volders PJ,Verheggen K,Menschaert G,Vandepoele K,Martens L,Vandesompele J,et al.An update on LNCipedia:a database for annotated human lncRNA sequences.Nucleic Acids Res 2015;43:D174-80.

    [58]Hofacker IL.Vienna RNA secondary structure server.Nucleic Acids Res 2003;31:3429-31.

    [59]Friedman RC,F(xiàn)arh KK,Burge CB,BartelDP.Most mammalian mRNAs are conserved targets of microRNAs.Genome Res 2009;19:92-105.

    [60]Welter D,MacArthur J,Morales J,Burdett T,Hall P,Junkins H,et al.The NHGRI GWAS Catalog,a curated resource of SNP-trait associations.Nucleic Acids Res 2014;42:D1001-6.

    [61]Volders PJ,Helsens K,Wang X,Menten B,Martens L,GevaertK,et al.LNCipedia:a database for annotated human lncRNA transcript sequences and structures.Nucleic Acids Res 2013;41: D246-51.

    [62]Hangauer MJ,Vaughn IW,McManus MT.Pervasive transcription of the human genome produces thousands of previously unidentified long intergenic noncoding RNAs.PLoS Genet 2013;9:e1003569.

    [63]Nielsen MM,Tehler D,Vang S,Sudzina F,Hedegaard J,Nordentoft I,et al.Identification of expressed and conserved human noncoding RNAs.RNA 2014;20:236-51.

    [64]Bhartiya D,Pal K,Ghosh S,Kapoor S,Jalali S,Panwar B,et al. lncRNome:a comprehensive knowledgebase of human long noncoding RNAs.Database(Oxford)2013;2013,bat034.

    [65]Bu D,Luo H,Jiao F,F(xiàn)ang S,Tan C,Liu Z,et al.Evolutionary annotation of conserved long non-coding RNAs in major mammalian species.Sci China Life Sci 2015;58:787-98.

    [66]Liu K,Yan Z,Li Y,Sun Z.Linc2GO:a human LincRNA function annotation resource based on ceRNA hypothesis. Bioinformatics 2013;29:2221-2.

    [67]Chan WL,Huang HD,Chang JG.lncRNAMap:a map of putative regulatory functions in the long non-coding transcriptome.Comput Biol Chem 2014;50:41-9.

    [68]Dinger ME,Pang KC,Mercer TR,Crowe ML,Grimmond SM,Mattick JS.NRED:a database of long noncoding RNA expression.Nucleic Acids Res 2009;37:D122-6.

    [69]Jiang Q,Wang J,Wu X,Ma R,Zhang T,Jin S,et al. LncRNA2Target:a database for differentially expressed genes after lncRNA knockdown or over expression.Nucleic Acids Res 2015;43:D193-6.

    [70]Zhou Z,Shen Y,Khan MR,Li A.LncReg:a reference resource for lncRNA-associated regulatory networks.Database(Oxford)2015;2015,bav083.

    [71]Paraskevopoulou MD,Georgakilas G,Kostoulas N,Reczko M,Maragkakis M,Dalamagas TM,et al.DIANA-LncBase:experimentally verified and computationally predicted microRNA targets on long non-coding RNAs.Nucleic Acids Res 2013;41: D239-45.

    [72]Dhiman H,Kapoor S,Sivadas A,Sivasubbu S,Scaria V. ZflncRNApedia:a comprehensive online resource for Zebrafish long non-coding RNAs.PLoS One 2015;10:e0129997.

    [73]Jin J,Liu J,Wang H,Wong L,Chua NH.PLncDB:plant long non-coding RNA database.Bioinformatics 2013;29:1068-71.

    [74]Qi P,Du X.The long non-coding RNAs,a new cancer diagnostic and therapeutic gold mine.Mod Pathol 2013;26:155-65.

    [75]Niazi F,Valadkhan S.Computational analysis of functional long noncoding RNAs reveals lack of peptide-coding capacity and parallels with 3′UTRs.RNA 2012;18:825-43.

    [76]Weirick T,John D,Dimmeler S,Uchida S.C-It-Loci:a knowledge database for tissue-enriched loci.Bioinformatics 2015;31:3537-43.

    [77]Zhao Z,Bai J,Wu A,Wang Y,Zhang J,Wang Z,et al.Co-LncRNA:investigating the lncRNA combinatorial effects in GO annotations and KEGG pathways based on human RNA-Seq data.Database(Oxford)2015;2015,bav082.

    [78]Luk AC,Gao H,Xiao S,Liao J,Wang D,Tu J,et al. GermlncRNA:a unique catalogue of long non-coding RNAs and associated regulations in male germ cell development.Database(Oxford)2015;2015,bav044.

    [79]Ning S,Zhao Z,Ye J,Wang P,Zhi H,Li R,et al.LincSNP:a database of linking disease-associated SNPs to human large intergenic non-coding RNAs.BMC Bioinformatics 2014;15:152.

    [80]Das S,Ghosal S,Sen R,Chakrabarti J.lnCeDB:database of human long noncoding RNA acting as competing endogenous RNA.PLoS One 2014;9:e98965.

    [81]Jiang Q,Ma R,Wang J,Wu X,Jin S,Peng J,et al. LncRNA2Function:a comprehensive resource for functional investigation of human lncRNAs based on RNA-seq data. BMC Genomics 2015;16:S2.

    [82]Ma L,Li A,Zou D,Xu X,Xia L,Yu J,et al.LncRNAWiki: harnessing community knowledge in collaborative curation of human long non-coding RNAs.Nucleic Acids Res 2015;43: D187-92.

    [83]Ning S,Zhao Z,Ye J,Wang P,Zhi H,Li R,et al. SNP@lincTFBS:an integrated database of polymorphisms in human LincRNA transcription factor binding sites.PLoS One 2014;9:e103851.

    [84]Jiang Q,Wang J,Wang Y,Ma R,Wu X,Li Y.TF2LncRNA: identifying common transcription factors for a list of lncRNA genes from ChIP-Seq data.Biomed Res Int 2014;2014:317642.

    *Corresponding author.

    E-mail:mstoll@uni-muenster.de(Stoll M).

    aORCID:0000-0001-7688-2394.

    bORCID:0000-0002-2711-4281.

    Peer review under responsibility of Beijing Institute of Genomics,Chinese Academy of Sciences and Genetics Society of China.

    http://dx.doi.org/10.1016/j.gpb.2016.03.001

    1672-0229?2016 The Authors.Production and hosting by Elsevier B.V.on behalf of Beijing Institute of Genomics,Chinese Academy of Sciences and Genetics Society of China.

    This is an open access article under the CC BY license(http://creativecommons.org/licenses/by/4.0/).

    svipshipincom国产片| 国产乱人伦免费视频| 日本在线视频免费播放| 午夜福利高清视频| 女性被躁到高潮视频| 女人爽到高潮嗷嗷叫在线视频| 精品久久久久久久久久免费视频| 欧美日韩一级在线毛片| 日日干狠狠操夜夜爽| 国产一卡二卡三卡精品| 成人国产一区最新在线观看| 欧美日韩黄片免| 午夜福利在线在线| ponron亚洲| 99国产精品一区二区三区| 美女扒开内裤让男人捅视频| 俺也久久电影网| 一边摸一边抽搐一进一小说| 欧美在线一区亚洲| 在线永久观看黄色视频| 午夜福利在线在线| 亚洲 欧美 日韩 在线 免费| 黄网站色视频无遮挡免费观看| 中文字幕久久专区| av片东京热男人的天堂| 母亲3免费完整高清在线观看| 在线天堂中文资源库| 久久久久久国产a免费观看| 日韩三级视频一区二区三区| 啦啦啦 在线观看视频| 国产黄a三级三级三级人| 日韩三级视频一区二区三区| 一二三四社区在线视频社区8| 侵犯人妻中文字幕一二三四区| 超碰成人久久| netflix在线观看网站| 国产精品久久电影中文字幕| 夜夜夜夜夜久久久久| 午夜精品在线福利| 真人一进一出gif抽搐免费| 97超级碰碰碰精品色视频在线观看| 国产欧美日韩一区二区精品| 日本 欧美在线| 久久九九热精品免费| 一二三四社区在线视频社区8| 村上凉子中文字幕在线| 亚洲 欧美 日韩 在线 免费| 国产亚洲av高清不卡| 中文字幕人妻熟女乱码| 亚洲精品中文字幕一二三四区| 天堂√8在线中文| 久久精品aⅴ一区二区三区四区| 日韩精品青青久久久久久| av免费在线观看网站| 男男h啪啪无遮挡| 亚洲欧洲精品一区二区精品久久久| 在线av久久热| 久久热在线av| 9191精品国产免费久久| 亚洲中文字幕日韩| 国产精品影院久久| 日韩精品青青久久久久久| 日本成人三级电影网站| 成人欧美大片| 国内精品久久久久久久电影| 国产高清有码在线观看视频 | 色播在线永久视频| 正在播放国产对白刺激| 在线看三级毛片| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品 欧美亚洲| 欧美日韩福利视频一区二区| 国产国语露脸激情在线看| 99热6这里只有精品| 国产单亲对白刺激| 制服丝袜大香蕉在线| 国产黄色小视频在线观看| 18禁美女被吸乳视频| 亚洲av美国av| 日韩精品免费视频一区二区三区| 午夜福利18| 在线观看66精品国产| 黄色毛片三级朝国网站| 免费一级毛片在线播放高清视频| aaaaa片日本免费| 18禁黄网站禁片午夜丰满| 亚洲九九香蕉| 夜夜爽天天搞| 亚洲国产欧美一区二区综合| 免费高清视频大片| 人人妻,人人澡人人爽秒播| 久久久久久九九精品二区国产 | 日韩欧美国产一区二区入口| 91国产中文字幕| 成人亚洲精品一区在线观看| av在线天堂中文字幕| 日本免费a在线| 一个人观看的视频www高清免费观看 | a在线观看视频网站| 日韩三级视频一区二区三区| 国产爱豆传媒在线观看 | 亚洲自偷自拍图片 自拍| 一卡2卡三卡四卡精品乱码亚洲| 非洲黑人性xxxx精品又粗又长| 久久久国产成人精品二区| 免费在线观看日本一区| 亚洲午夜精品一区,二区,三区| 精品电影一区二区在线| 变态另类丝袜制服| 久久精品成人免费网站| 亚洲成a人片在线一区二区| 久久 成人 亚洲| 一区二区三区高清视频在线| 天堂√8在线中文| 桃红色精品国产亚洲av| 欧美成人免费av一区二区三区| 国产精品亚洲美女久久久| av天堂在线播放| 国产一区二区三区在线臀色熟女| 嫁个100分男人电影在线观看| 国产久久久一区二区三区| 99久久99久久久精品蜜桃| 亚洲自偷自拍图片 自拍| av中文乱码字幕在线| 90打野战视频偷拍视频| 一区二区三区精品91| 亚洲av电影不卡..在线观看| 校园春色视频在线观看| 熟女少妇亚洲综合色aaa.| 欧美一级a爱片免费观看看 | 欧美国产精品va在线观看不卡| 精品一区二区三区视频在线观看免费| 妹子高潮喷水视频| 日韩欧美 国产精品| 亚洲国产精品999在线| 女生性感内裤真人,穿戴方法视频| 精品乱码久久久久久99久播| 国产午夜精品久久久久久| 精品国产乱码久久久久久男人| 亚洲狠狠婷婷综合久久图片| 叶爱在线成人免费视频播放| 成人午夜高清在线视频 | aaaaa片日本免费| 亚洲国产欧美一区二区综合| 国产一区二区三区在线臀色熟女| 日韩欧美一区二区三区在线观看| 波多野结衣高清无吗| 亚洲中文字幕一区二区三区有码在线看 | 亚洲人成网站在线播放欧美日韩| www.精华液| 好男人在线观看高清免费视频 | 真人做人爱边吃奶动态| 精品午夜福利视频在线观看一区| 一区二区三区精品91| 国产熟女午夜一区二区三区| 欧美丝袜亚洲另类 | 国产区一区二久久| 亚洲av电影不卡..在线观看| 麻豆国产av国片精品| 黄色视频,在线免费观看| 久久久水蜜桃国产精品网| 在线观看一区二区三区| 日韩欧美 国产精品| 女同久久另类99精品国产91| 热re99久久国产66热| 亚洲男人的天堂狠狠| 亚洲免费av在线视频| 又大又爽又粗| 国产亚洲欧美在线一区二区| 日韩 欧美 亚洲 中文字幕| 99久久久亚洲精品蜜臀av| 国产精品一区二区精品视频观看| 丝袜在线中文字幕| 在线天堂中文资源库| 一边摸一边抽搐一进一小说| 99久久精品国产亚洲精品| 中文字幕久久专区| 日本黄色视频三级网站网址| 99国产极品粉嫩在线观看| 亚洲美女黄片视频| 哪里可以看免费的av片| 一二三四社区在线视频社区8| 一进一出抽搐gif免费好疼| 波多野结衣高清作品| 午夜免费激情av| 制服人妻中文乱码| 国产一区二区激情短视频| videosex国产| 久久这里只有精品19| 亚洲国产欧美网| 午夜福利在线在线| 欧美性猛交╳xxx乱大交人| 99久久精品国产亚洲精品| 好看av亚洲va欧美ⅴa在| 99re在线观看精品视频| 国产又爽黄色视频| 久久婷婷人人爽人人干人人爱| 色播亚洲综合网| 成熟少妇高潮喷水视频| 久久精品亚洲精品国产色婷小说| 国产精品久久久久久亚洲av鲁大| 女性生殖器流出的白浆| 一进一出抽搐动态| 成人欧美大片| 亚洲真实伦在线观看| 亚洲九九香蕉| 国产单亲对白刺激| 不卡av一区二区三区| 搡老岳熟女国产| 欧美成人性av电影在线观看| 精品欧美一区二区三区在线| 中文字幕精品免费在线观看视频| 日韩成人在线观看一区二区三区| 免费看美女性在线毛片视频| 成人三级做爰电影| 亚洲精品在线美女| а√天堂www在线а√下载| 亚洲一区高清亚洲精品| 国产精品1区2区在线观看.| 99在线人妻在线中文字幕| 欧美乱妇无乱码| 国产成人精品久久二区二区免费| 亚洲,欧美精品.| 欧美黄色淫秽网站| 成人国产一区最新在线观看| tocl精华| 亚洲国产精品合色在线| 久久中文字幕人妻熟女| 国产激情偷乱视频一区二区| 黄片小视频在线播放| 久久久国产成人精品二区| 亚洲精品中文字幕在线视频| 免费看十八禁软件| 亚洲午夜精品一区,二区,三区| 国产一区二区激情短视频| svipshipincom国产片| 人人妻,人人澡人人爽秒播| 99riav亚洲国产免费| 无限看片的www在线观看| 国产欧美日韩精品亚洲av| 精品久久久久久成人av| 亚洲国产精品999在线| 99在线视频只有这里精品首页| 曰老女人黄片| 精品免费久久久久久久清纯| 在线观看舔阴道视频| 国产成年人精品一区二区| 色av中文字幕| 在线视频色国产色| 国产日本99.免费观看| 久久久久久免费高清国产稀缺| 国产成人精品无人区| 免费一级毛片在线播放高清视频| 18禁国产床啪视频网站| svipshipincom国产片| 国产亚洲精品第一综合不卡| 国产极品粉嫩免费观看在线| 国产成人欧美| 在线观看免费日韩欧美大片| 国产精品98久久久久久宅男小说| 亚洲av熟女| 亚洲中文日韩欧美视频| 日韩三级视频一区二区三区| 91老司机精品| 亚洲激情在线av| 成年女人毛片免费观看观看9| 每晚都被弄得嗷嗷叫到高潮| 麻豆av在线久日| 超碰成人久久| 日韩欧美免费精品| 国产亚洲精品av在线| 天堂动漫精品| 久久午夜综合久久蜜桃| 久久精品aⅴ一区二区三区四区| 99久久精品国产亚洲精品| 丝袜在线中文字幕| 制服诱惑二区| 国产99白浆流出| 国产精品 国内视频| 成人永久免费在线观看视频| 国产精品精品国产色婷婷| 精品国内亚洲2022精品成人| 国产私拍福利视频在线观看| 日本在线视频免费播放| 久久精品aⅴ一区二区三区四区| 欧美丝袜亚洲另类 | 女人爽到高潮嗷嗷叫在线视频| 久久久久久久久久黄片| 国产片内射在线| 男女午夜视频在线观看| 国产真实乱freesex| 中文亚洲av片在线观看爽| 欧美在线一区亚洲| 高清在线国产一区| 999精品在线视频| 亚洲精品美女久久久久99蜜臀| 99riav亚洲国产免费| 99国产精品99久久久久| 国产成人一区二区三区免费视频网站| 后天国语完整版免费观看| 久久草成人影院| 中文亚洲av片在线观看爽| 一级毛片高清免费大全| 在线观看一区二区三区| 国产亚洲av嫩草精品影院| 国产免费男女视频| 亚洲精品美女久久久久99蜜臀| www国产在线视频色| 999精品在线视频| av视频在线观看入口| svipshipincom国产片| 亚洲电影在线观看av| 中文字幕人妻丝袜一区二区| 欧美中文日本在线观看视频| 一边摸一边做爽爽视频免费| 国产一区二区三区在线臀色熟女| 成人亚洲精品av一区二区| 国产精品免费视频内射| 国内揄拍国产精品人妻在线 | 久久中文字幕一级| 99精品久久久久人妻精品| 国产精品永久免费网站| 国产亚洲精品av在线| 欧美三级亚洲精品| 在线观看免费视频日本深夜| av中文乱码字幕在线| 此物有八面人人有两片| 少妇被粗大的猛进出69影院| 俄罗斯特黄特色一大片| 别揉我奶头~嗯~啊~动态视频| 亚洲第一电影网av| 日本熟妇午夜| 日本五十路高清| 国产精品九九99| 成年人黄色毛片网站| 国内久久婷婷六月综合欲色啪| 亚洲自偷自拍图片 自拍| 男人舔奶头视频| 99在线人妻在线中文字幕| 中文字幕人妻熟女乱码| 色播在线永久视频| 久久精品夜夜夜夜夜久久蜜豆 | 最近最新免费中文字幕在线| 国产熟女xx| 在线天堂中文资源库| 午夜久久久久精精品| 嫩草影视91久久| 亚洲av熟女| 免费一级毛片在线播放高清视频| 一级a爱片免费观看的视频| 男女床上黄色一级片免费看| 国产亚洲精品久久久久5区| 久久久久九九精品影院| 色综合欧美亚洲国产小说| 中文字幕av电影在线播放| 免费无遮挡裸体视频| 午夜福利高清视频| 久久天堂一区二区三区四区| 国产精品av久久久久免费| 久久久久九九精品影院| 久久精品国产99精品国产亚洲性色| 国产精品一区二区免费欧美| 亚洲最大成人中文| www.自偷自拍.com| 18禁观看日本| 欧美午夜高清在线| 制服丝袜大香蕉在线| 欧美人与性动交α欧美精品济南到| 成人三级做爰电影| 亚洲av片天天在线观看| 自线自在国产av| 国产精品免费一区二区三区在线| 叶爱在线成人免费视频播放| 精品欧美国产一区二区三| 亚洲av电影不卡..在线观看| 精品久久久久久久人妻蜜臀av| 变态另类成人亚洲欧美熟女| 免费看a级黄色片| 19禁男女啪啪无遮挡网站| 麻豆av在线久日| 国产v大片淫在线免费观看| 精品久久久久久久人妻蜜臀av| 亚洲成人免费电影在线观看| 搡老岳熟女国产| 国产一区二区三区在线臀色熟女| 嫁个100分男人电影在线观看| 亚洲成国产人片在线观看| 精品国内亚洲2022精品成人| 女同久久另类99精品国产91| 巨乳人妻的诱惑在线观看| 久久久久久大精品| 欧美乱妇无乱码| 国产人伦9x9x在线观看| 女性被躁到高潮视频| 2021天堂中文幕一二区在线观 | 色尼玛亚洲综合影院| 免费在线观看视频国产中文字幕亚洲| ponron亚洲| 长腿黑丝高跟| 成人手机av| 91大片在线观看| 国产黄片美女视频| 丁香欧美五月| 欧美性猛交╳xxx乱大交人| 99精品久久久久人妻精品| 真人做人爱边吃奶动态| 老汉色∧v一级毛片| 香蕉av资源在线| 亚洲五月天丁香| 一本综合久久免费| 亚洲精品美女久久久久99蜜臀| a在线观看视频网站| 午夜成年电影在线免费观看| 超碰成人久久| 亚洲成人精品中文字幕电影| 免费女性裸体啪啪无遮挡网站| x7x7x7水蜜桃| 午夜日韩欧美国产| 天天躁夜夜躁狠狠躁躁| 久久久水蜜桃国产精品网| 青草久久国产| 久久国产精品人妻蜜桃| 亚洲人成网站在线播放欧美日韩| 国产主播在线观看一区二区| 午夜亚洲福利在线播放| 久久这里只有精品19| 99国产精品99久久久久| 自线自在国产av| 国产三级黄色录像| 国产精品美女特级片免费视频播放器 | 男女床上黄色一级片免费看| 国产亚洲精品久久久久5区| 国产日本99.免费观看| 中文字幕高清在线视频| 日本免费一区二区三区高清不卡| 亚洲成av片中文字幕在线观看| 国产成人影院久久av| www.熟女人妻精品国产| 婷婷精品国产亚洲av| 岛国视频午夜一区免费看| 一区二区三区精品91| 91大片在线观看| 亚洲国产高清在线一区二区三 | 在线观看免费日韩欧美大片| 精品久久久久久久久久免费视频| 伊人久久大香线蕉亚洲五| 搡老妇女老女人老熟妇| 高清在线国产一区| 欧美亚洲日本最大视频资源| or卡值多少钱| 欧美成人免费av一区二区三区| 久久香蕉精品热| 国产精品日韩av在线免费观看| 国产精品亚洲av一区麻豆| 精品福利观看| 99在线人妻在线中文字幕| 国产97色在线日韩免费| 国产精品99久久99久久久不卡| www.999成人在线观看| 两人在一起打扑克的视频| 国产精品1区2区在线观看.| 一二三四社区在线视频社区8| 97超级碰碰碰精品色视频在线观看| 免费无遮挡裸体视频| 免费高清视频大片| 成人国语在线视频| 日本黄色视频三级网站网址| 又大又爽又粗| 午夜两性在线视频| 大香蕉久久成人网| 国产精品1区2区在线观看.| 亚洲激情在线av| 亚洲人成伊人成综合网2020| 亚洲午夜精品一区,二区,三区| 亚洲成av片中文字幕在线观看| 亚洲成a人片在线一区二区| 变态另类成人亚洲欧美熟女| 美女高潮到喷水免费观看| 成人亚洲精品av一区二区| 在线观看www视频免费| 免费看美女性在线毛片视频| 91大片在线观看| 国产三级在线视频| 嫁个100分男人电影在线观看| 亚洲av成人一区二区三| 国产蜜桃级精品一区二区三区| 国产国语露脸激情在线看| 老司机在亚洲福利影院| 亚洲熟妇中文字幕五十中出| 久久草成人影院| 悠悠久久av| 欧美在线一区亚洲| 看片在线看免费视频| 日本 av在线| 女警被强在线播放| 亚洲欧美精品综合一区二区三区| videosex国产| 亚洲欧洲精品一区二区精品久久久| 神马国产精品三级电影在线观看 | 日韩欧美免费精品| 淫妇啪啪啪对白视频| 美女午夜性视频免费| 免费在线观看黄色视频的| 国产精品98久久久久久宅男小说| 久久久国产精品麻豆| 级片在线观看| 久久久国产欧美日韩av| tocl精华| 午夜激情av网站| 黄色视频,在线免费观看| 黄片播放在线免费| 一a级毛片在线观看| 99久久国产精品久久久| 精品不卡国产一区二区三区| 国内毛片毛片毛片毛片毛片| 我的亚洲天堂| 日韩欧美一区视频在线观看| 国产日本99.免费观看| 国产一级毛片七仙女欲春2 | 久久久久免费精品人妻一区二区 | 色精品久久人妻99蜜桃| 国产熟女xx| 日韩三级视频一区二区三区| 午夜成年电影在线免费观看| 精品人妻1区二区| 99riav亚洲国产免费| 在线十欧美十亚洲十日本专区| 国产亚洲欧美在线一区二区| 久久午夜综合久久蜜桃| 国产一区在线观看成人免费| 非洲黑人性xxxx精品又粗又长| 亚洲 欧美 日韩 在线 免费| 亚洲av第一区精品v没综合| 99国产精品99久久久久| 亚洲av美国av| 国产久久久一区二区三区| 91成人精品电影| 18禁黄网站禁片午夜丰满| 动漫黄色视频在线观看| 一个人观看的视频www高清免费观看 | 一级毛片女人18水好多| 一进一出好大好爽视频| 亚洲欧美精品综合久久99| 身体一侧抽搐| 桃红色精品国产亚洲av| 欧美不卡视频在线免费观看 | 久久久久久大精品| 麻豆一二三区av精品| 老汉色av国产亚洲站长工具| 亚洲欧洲精品一区二区精品久久久| 久久国产精品男人的天堂亚洲| 日韩一卡2卡3卡4卡2021年| 日韩欧美 国产精品| 桃红色精品国产亚洲av| 国产精品98久久久久久宅男小说| 精品电影一区二区在线| 色精品久久人妻99蜜桃| 国产精品国产高清国产av| 免费观看人在逋| 9191精品国产免费久久| 黄色视频不卡| 中文字幕久久专区| 88av欧美| 观看免费一级毛片| 国产免费av片在线观看野外av| 婷婷丁香在线五月| 亚洲欧美日韩无卡精品| 午夜激情福利司机影院| 国产爱豆传媒在线观看 | 又大又爽又粗| av免费在线观看网站| 69av精品久久久久久| 久久精品亚洲精品国产色婷小说| 亚洲av电影在线进入| 国产1区2区3区精品| 变态另类丝袜制服| 免费观看精品视频网站| 88av欧美| 精品久久久久久久久久免费视频| 宅男免费午夜| 制服人妻中文乱码| 国产私拍福利视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 男人舔女人下体高潮全视频| 99久久综合精品五月天人人| 国产亚洲精品第一综合不卡| 日韩一卡2卡3卡4卡2021年| 日韩精品青青久久久久久| 人人妻人人澡人人看| 老熟妇仑乱视频hdxx| 国产真实乱freesex| 男人舔女人的私密视频| 非洲黑人性xxxx精品又粗又长| 韩国精品一区二区三区| 久久久精品国产亚洲av高清涩受| 久久精品人妻少妇| 中文字幕精品亚洲无线码一区 | 长腿黑丝高跟| 啦啦啦韩国在线观看视频| 国产午夜精品久久久久久| 免费一级毛片在线播放高清视频| 久久久久久久久久黄片| 淫妇啪啪啪对白视频| 长腿黑丝高跟| 国产1区2区3区精品| 午夜福利免费观看在线| 9191精品国产免费久久| 欧美性猛交╳xxx乱大交人| av电影中文网址| 性色av乱码一区二区三区2| 少妇熟女aⅴ在线视频| 51午夜福利影视在线观看| 欧美激情久久久久久爽电影| 夜夜爽天天搞| 成人特级黄色片久久久久久久|