SEARCH · PubMed Health

Results for “Small noncoding RNA”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21Linked to original sources

PubMed2002

The chicken beta-globin insulator element conveys chromatin boundary activity but not imprinting at the mouse Igf2/H19 domain.

Imprinting of the mouse insulin-like growth factor 2 (Igf2) and H19 genes is regulated by an imprinting control region (ICR). The hypomethylated maternal copy functions as a chromatin insulator through the binding of CTCF and prevents Igf2 activation in cis, while hypermethylation of the paternal copy inactivates insulator function and leads to inactivation of H19 in cis. The specificity of the ICR sequence for mediating imprinting and chromatin insulation was investigated by substituting it for two copies of the chicken beta-globin insulator element, (Ch beta GI)(2), in mice. This introduced sequence resembles the ICR in size, and in containing CTCF-binding sites and CpGs, but otherwise lacks homology. On maternal inheritance, the (Ch beta GI)(2) was hypomethylated and displayed full chromatin insulator activity. Monoallelic expression of Igf2 and H19 was retained and mice were of normal size. These results suggest that the ICR sequence, aside from CTCF-binding sites, is not uniquely specialized for chromatin insulation at the Igf2/H19 region. On paternal inheritance, the (Ch beta GI)(2) was also hypomethylated and displayed strong insulator activity--fetuses possessed very low levels of Igf2 RNA and were greatly reduced in size, being as small as Igf2-null mutants. Furthermore, the paternal H19 allele was active. These results suggest that differential ICR methylation in the female and male germ lines is not acquired through differential binding of CTCF. Rather, it is likely to be acquired through a separate or downstream process.

Animals↗

PubMed1997

A 1.1-Mb transcript map of the hereditary hemochromatosis locus.

In the process of positionally cloning a candidate gene responsible for hereditary hemochromatosis (HH), we constructed a 1.1-Mb transcript map of the region of human chromosome 6p that lies 4.5 Mb telomeric to HLA-A. A combination of three gene-finding techniques, direct cDNA selection, exon trapping, and sample sequencing, were used initially for a saturation screening of the 1.1-Mb region for expressed sequence fragments. As genetic analysis further narrowed the HH candidate locus, we sequenced completely 0.25 Mb of genomic DNA as a final measure to identify all genes. Besides the novel MHC class 1-like HH candidate gene HLA-H, we identified a family of five butyrophilin-related sequences, two genes with structural similarity to a type 1 sodium phosphate transporter, 12 novel histone genes, and a gene we named RoRet based on its strong similarity to the 52-kD Ro/SSA lupus and Sjogren's syndrome auto-antigen and the RET finger protein. Several members of the butyrophilin family and the RoRet gene share an exon of common evolutionary origin called B30-2. The B30-2 exon was originally isolated from the HLA class 1 region, yet has apparently "shuffled" into several genes along the chromosome telomeric to the MHC. The conservation of the B30-2 exon in several novel genes and the previously described amino acid homology of HLA-H to MHC class 1 molecules provide further support that this gene-rich region of 6p21.3 is related to the MHC. Finally, we performed an analysis of the four approaches for gene finding and conclude that direct selection provides the most effective probes for cDNA screening, and that as much as 30% of ESTs in this 1.1-Mb region may be derived from noncoding genomic DNA.

Amino Acid Sequence↗

PubMed1998

GBV-C/HGV infection in hepatitis C virus-infected deferred Swedish blood donors.

Sera from 62 hepatitis C virus (HCV)-infected Swedish blood donors were tested by a nested polymerase chain reaction using primers targeting the 5'-noncoding region of the GB virus-C/hepatitis G (GBV-C/HGV) genome and an enzyme-linked immunosorbent assay that detects antibodies to the envelope protein E2 of GBV-C/HGV (anti-E2). Fourteen (22%) and 21 (34%) of the 62 blood donors were found to be GBV-C/HGV RNA and anti-E2 positive, respectively. None of the blood donors was positive for both GBV-C/HGV RNA and anti-E2. Thus, 35 of 62 (56%) HCV-infected donors had been exposed to GBV-C/HGV infection. At sequencing of the 14 GBV-C/HGV isolates, 12 were identified as subtype 2a and 2 as subtype 2b. One of 7 (14%) donors with mild liver disease such as steatosis and nonspecific reactive hepatitis had been exposed to GBV-C/HGV vs. 34 of 55 (62%) with chronic hepatitis with or without cirrhosis (P = 0.04). All other differences in histology were small between HCV and dual HCV GBV-C/HGV-infected donors. In conclusion, more than half of HCV-infected Swedish blood donors in this study were positive for either GBV-C/HGV RNA or anti-E2. GBV-C/HGV viremia and seropositivity were mutually exclusive.

Adult↗

PubMed1996

Human phenol sulfotransferase gene contains two alternative promoters: Structure and expression of the gene.

Phenol sulfotransferases catalyze the transfer of a sulfonate moiety from 3'-phosphoadenosyl 5'-phosphosulfate to a phenolic group of lipophylic substrates to generate soluble sulfate esters. Using a phenol sulfotransferase cDNA as probe to screen a human leukocyte genomic DNA library constructed in lambda EMBL3, we obtained a clone containing a complete gene sequence. Comparison of the gene sequence with that of the corresponding cDNAs, namely phenol-sulfating phenol sulfotransferase (P-PST) or thermostable sulfotransferase (TS-PST), and human aryl sulfotransferase 1 and 2 (HAST1 and HAST2) indicates that the gene possesses eight short exons separated by seven introns included in approximately 5 kb. HAST2 has a different 5' untranslated sequence, and thus is encoded by a different mRNA species. While the nucleotide sequence corresponding to the 5' noncoding region of P-PST (TS-PST and HAST1) is included in the exon I, the 5' untranslated sequence of HAST2 is located in the beginning of exon IIa. The remaining sequence in exon II that is identical to both P-PST and HAST2 was termed exon IIb. Exons III to VIII, which cover the coding region and the 3' untranslated region, are almost identical in all types of PST or AST cDNAs. These results suggest that the phenol sulfotransferase gene possesses two alternate promoters that drive the expression of the two different mRNA species in a tissue-specific manner. Transfection of chloramphenicol acetyl transferase (CAT) reporter gene vectors containing the 5'-flanking sequence upstream from exon I and exon II, respectively, in transformed human embryonal kidney (293) cells indicate that both sequences possess promoter activity with higher activity for promoter 1. RNA blot analysis indicates that human phenol sulfotransferase gene is expressed in kidney, liver, lung, leukocyte, colon, small intestine, and spleen.

Amino Acid Sequence↗

PubMed1988

Sigma region located between C mu and C delta genes of human immunoglobulin heavy chain: possible involvement of tRNA-like structure in RNA splicing.

Noncoding regions within the cluster of immunoglobulin heavy chain constant genes in the human genome contained a number of repeats. In the mu-delta intron, two repeating units were contained. One 442-base-long fragment located JH-mu intron (defined as "sigma mu(sigma mu)") occupied the position in the mu-delta intron. The other 1166-base-long fragment located somewhere in front of S (class switch) region of C gamma gene was also found in the mu-delta intron. We defined the repeats in the mu-delta intron as "SIGMA (sigma)". The polarities of the longer repeats in the genome were opposite between the mu-delta intron and the upstreams of C gamma genes. These inverted copies (defined as sigma gamma 3 and sigma gamma 4), located 6 kb upstream of their respective C gamma's, were apparently transcribed in vitro, via RNA polymerase III and transcripts should have contained tRNA-like structures. Small DNA fragments capable of encoding tRNA-like structures were also found in corresponding regions of mouse Ig C gamma cluster.

Animals↗

PubMed1989

B-G cDNA clones have multiple small repeats and hybridize to both chicken MHC regions.

We used rabbit antisera to the chicken MHC erythrocyte molecule B-G and to the class I alpha chain (B-F) to screen lambda gt11 cDNA expression libraries made with RNA selected by oligo-dT from bone marrow cells of anemic B19 homozygous chickens. Eight clones were found to encode B-G molecules which hybridize with sequences in the chicken MHC as defined by congenic strains; the fusion proteins react with multiple immune but not preimmune sera, they select antibodies from the antisera to B-G, which then react with distinct erythrocyte B-G protein patterns, and they elicit antibodies from mice which in turn react with authentic B-G proteins. None of the clones represent a complete message, some--if not all--bear introns, and none of them match with any sequences presently stored in the data banks. The following new information did, however, emerge. At least two homologous transcripts are present in this homozygous chicken, thereby formally proving the existence of an expressed multigene family. The 3' ends (3'UT) are simple sequences with 80% nucleotide identity between clones, while the 5' ends (either coding or noncoding) are composed of multiple short repeats which are far less similar. These repeats could explain the bewildering variation in size of B-G proteins within and between haplotypes. Southern blots of genomic chicken DNA gave complex patterns for most probes, with many bands in common using different probes, but few bands in common between haplotypes. The sequences detected are all present in the MHC, based on the congenic lines CB and CC. Most of these sequences map into the B-G region, but some map into the B-F/B-L region as defined by the haplotypes B15, B21, and their apparently reciprocal recombinants B21r3 and B15r1.

Animals↗

PubMed1995

The small cysteine-rich protein P14 of beet necrotic yellow vein virus regulates accumulation of RNA 2 in cis and coat protein in trans.

The effect of null mutations of the small cysteine-rich protein P14 encoded by RNA 2 of beet necrotic yellow vein virus has been investigated using in vitro transcripts of viral RNA to infect Chenopodium quinoa protoplasts. The P14 mutations down-regulated RNA 2 accumulation by approximately 10- to 50-fold. Accumulation of minus-strand RNA 2 was also diminished but RNA 1 accumulation was much less affected. The inhibition of RNA 2 accumulation could not be complemented in trans by providing P14 from another source (either a second molecule of RNA 2 or an RNA 3-based replicon) containing and expressing the P14 gene. The P14 null mutations dramatically inhibited accumulation of viral coat protein, which is encoded by the 5'-proximal gene on RNA 2, but this effect could be complemented in trans, indicating that it occurs by a mechanism distinct from that affecting RNA 2 accumulation. Transient expression experiments were also carried out in which a plasmid expressing P14 and plasmids expressing a reporter gene placed downstream of potential translational control sequences (the 5'-noncoding sequences of RNAs 2, 3, or 4) were introduced into C. quinoa or Nicotiana tabacum leaves by microprojectile bombardment. Coexpression of P14 produced a 3- to 4-fold stimulation of reporter gene expression levels for all the constructs. The lack of sequence specificity suggests that this phenomenon is not directly related to the RNA 2-specific stimulation of coat protein accumulation observed in a viral infection.

Base Sequence↗

PubMed1997

bic, a novel gene activated by proviral insertions in avian leukosis virus-induced lymphomas, is likely to function through its noncoding RNA.

The bic locus is a common retroviral integration site in avian leukosis virus (ALV)-induced B-cell lymphomas originally identified by infection of chickens with ALVs of two different subgroups (Clurman and Hayward, Mol. Cell. Biol. 9:2657-2664, 1989). Based on its frequent association with c-myc activation and its preferential activation in metastatic tumors, the bic locus is thought to harbor a gene that can collaborate with c-myc in lymphomagenesis and presumably plays a role in late stages of tumor progression. In the present study, we have cloned and characterized two novel genes, bdw and bic, at the bic locus. bdw encoded a putative novel protein of 345 amino acids. However, its expression did not appear to be altered in tumor tissues, suggesting that it is not involved in oncogenesis. The bic gene consisted of two exons and was expressed as two spliced and alternatively polyadenylated transcripts at low levels in lymphoid/hematopoietic tissues. In tumors harboring bic integrations, proviruses drove bic gene expression by promoter insertion, resulting in high levels of expression of a chimeric RNA containing bic exon 2. Interestingly, bic lacked an extensive open reading frame, implying that it may function through its RNA. Computer analysis of RNA from small exon 2 of bic predicted extensive double-stranded structures, including a highly ordered RNA duplex between nucleotides 316 and 461. The possible role of bic in cell growth and differentiation is discussed in view of the emerging evidence that untranslated RNAs play a role in growth control.

Amino Acid Sequence↗

PubMed1996

Antibody reactivity profiles following immunization with diverse peptides of the PERB11 (MIC) family.

PERB11 (MIC) is a gene family possessing multiple copies located within the MHC. Structurally, PERB11 is related to the MHC class I, neonatal IgG Fc receptor (FcRn) and Zn-alpha 2-glycoprotein molecules. The MHC class I family is complex in terms of its genomic arrangement, expression and function, and available evidence suggests that the PERB11 family may be similarly complex. We have adopted an approach to study the expression of such complex gene families by immunizing with multiple peptides and by screening the resulting antibodies against a large range of tissues. The amino acid sequences of PERB11.1 and PERB11.2 as well as those of other related molecules were analysed and compared. Peptides were chosen for immunization based upon (i) loop formation within the equivalent known structure of the MHC class I molecules; (ii) immunogenicity by computer analysis; and (iii) evolutionary relationships. Antibodies in serum from immunized rabbits bound to three out of six peptides used for immunization. ELISA and immunoprecipitation demonstrated binding both to the peptides and to the PERB11.2 recombinant protein. By immunofluorescent staining of various tissues of several species, the three antisera generated overlapping profiles of activity. These included reactions with kidney, small and large intestine, oesophagus, testis, ovary and human neutrophils. This is the first description of antibodies induced by the PERB11 peptides. The extreme complexity of these profiles requires further investigation, but may be explained in terms of antibodies against diverse products of the PERB11 gene family and/or related molecules.

Amino Acid Sequence↗

PubMed1991

The primary structure of human glutaminyl-tRNA synthetase. A highly conserved core, amino acid repeat regions, and homologies with translation elongation factors.

We describe the nucleotide sequences of several overlapping cDNA clones specific for human glutaminyl-tRNA synthetase. The identified open reading frame indicates that the enzyme is composed of 1440 amino acids. A stretch of about 360 amino acids of the human enzyme is highly conserved in bacterial and yeast glutaminyl-tRNA synthetases. However, the human enzyme is three times larger than the bacterial and twice as large as the yeast enzyme suggesting that a considerable part of human glutaminyl-tRNA synthetase has evolved to perform functions other than the charging of tRNA. The sequence outside of the conserved core region includes three 57-amino acid repeats followed by a consecutive stretch of 11 charged amino acids. A computer assisted search of two protein data banks reveals that the human glutaminyl-tRNA synthetase shares small blocks of amino acid similarities with several other synthetases of different amino acid specificities. Interestingly, the enzyme also possesses some regions of similarities with eukaryotic translation elongation factor EF-1 but not with any other sequence stored in the protein data banks. The coding regions of human and mouse glutaminyl-tRNA synthetase cDNAs are identical at 94% of the codons. However, the 3'-noncoding regions of mouse and human mRNAs are more divergent (approximately 68%) but both possess the potential to form stable secondary structures of similar general architecture.

Amino Acid Sequence↗

PubMed2000

Complex structure and regulation of expression of the rat gene for inward rectifier potassium channel Kir7.1.

Genomic organization of the rat inward rectifier K(+) channel Kir7.1 was determined in an attempt to clarify how multiple species of its mRNA are generated in a tissue-specific manner and how its expression is regulated. The rat Kir7.1 gene spans >40 kilobases (kb) and consists of eight exons; the first four exons encode the 5'-untranslated region that is unusually long ( approximately 3 kb). The coding region is located in exons 5 and 6. In the testis, exon 4 is processed as four exons (4a-4d), whereas it is recognized as a single exon in the small intestine. The three major species of rat Kir7.1 mRNA (1.4, 2.2, and 3.2 kb) were found to arise from alternative usage of the two promoters and polyadenylation signals and by alternative splicing of the 5'-noncoding exons. The splicing pattern of the 5'-noncoding exons is quite complex and highly tissue-specific, suggesting that complex mechanisms may operate to regulate the Kir7.1 expression. Deletion and mutational analysis of the promoter activity indicated that the rat Kir7.1 gene is regulated by cAMP through a CCAAT element. The cAMP induction was also demonstrated using the rat follicular cell line FRTL-5 endogenously expressing Kir7.1.

5' Untranslated Regions↗

PubMed2006

A novel class of small RNAs in mouse spermatogenic cells.

Small noncoding RNAs, including small interfering RNAs (siRNAs) and micro RNAs (miRNAs) of approximately 21 nucleotides (nt) in length, have emerged as potent regulators of gene expression at both transcriptional and post-transcriptional levels in diverse organisms. Here we report the identification of a novel class of small RNAs in the mouse male germline termed piwi-interacting RNAs (piRNAs). piRNAs are approximately 30 nt in length. They are expressed during spermatogenesis, mostly in spermatids. piRNAs are associated with MIWI, a spermatogenesis-specific PIWI subfamily member of the Argonaute protein family, and depend on MIWI for their biogenesis and/or stability. Furthermore, a subpopulation of piRNAs are associated with polysomes, suggesting their potential role in translational regulation.

Animals↗

PubMedDate not supplied

Turner syndrome female with a small ring X chromosome lacking the XIST, an unexpectedly mild phenotype and an atypical association with alopecia universalis.

Rearranged X chromosome in Turner syndrome (TS) are generally well tolerated but in cases of ring X chromosomes and of X/autosome translocations the incidence of mental retardation and other congenital abnormalities can be significantly higher. These abnormal phenotypes can be ascribed to failed or partial X inactivation. Here, we report a 10-year-old female who was referred for a cytogenetic analysis because she developed an alopecia universalis. The patient, of normal intelligence, had been found to have traits of TS, especially short stature. A first cytogenetic analysis showed a no mosaic 45,X karyotype. Since, the risk of developing gonadoblastoma in TS patients with mosaicism for a Y derivative chromosome and because association of alopecia universalis and TS is uncommon, fluorescence in situ hybridization (FISH) was performed to search for a second cell population. Our patient was found to have a mosaic 45,X/46,X,+r. FISH analysis using sex chromosome probes permitted us to identify the very small marker as a ring X chromosome, detected in 90% of cells. The ring appeared to be formed almost totally of alphoid sequences with breakpoints in the juxtacentromeric region. The r(X) does not include the XIST locus and may, therefore, not be subject to X-inactivation. Unexpectedly mild phenotype in our patient and its association with alopecia universalis will be discussed.

Abnormalities, Multiple↗

PubMed1987

The 3' untranslated region of the human interferon-beta mRNA has an inhibitory effect on translation.

In vitro-transcribed human interferon-beta (IFN-beta) mRNA, which contains all the sequence of the natural molecule, is poorly translated in a reticulocyte lysate or when injected in Xenopus oocytes. This low level of translation is due to an inhibition by the 5' and 3' untranslated regions (UTRs). Indeed, the replacement of these sequences by those of Xenopus beta-globin mRNA dramatically increases the translational efficiency of the mRNA, especially in oocytes. This phenomenon is not due to a difference in mRNA stability since both native and chimeric mRNAs remain undegraded, at least during the translation period considered. Construction of different chimeric molecules having various combinations of 5' and 3' UTRs from IFN-beta or Xenopus beta-globin mRNA or a small sequence of SP6 polylinker as 5' UTR has revealed that the 3' UTR of IFN-beta in itself has a pronounced inhibitory effect on translation in the two translation systems from animal cells. Indeed, the addition of this 3' UTR at the 3' end of the coding region of a chicken lysozyme mRNA also causes a large decrease of its translational capacity in both systems. However, the nature of the 5' noncoding sequence influences the degree of translation inhibition exerted by the 3' UTR. Remarkably, we observed no difference in translation level when the different mRNAs were tested in a wheat germ extract.

Animals↗

PubMed2006

X chromosomal abnormalities in basal-like human breast cancer.

Sporadic basal-like cancers (BLC) are a distinct class of human breast cancers that are phenotypically similar to BRCA1-associated cancers. Like BRCA1-deficient tumors, most BLC lack markers of a normal inactive X chromosome (Xi). Duplication of the active X chromosome and loss of Xi characterized almost half of BLC cases tested. Others contained biparental but nonheterochromatinized X chromosomes or gains of X chromosomal DNA. These abnormalities did not lead to a global increase in X chromosome transcription but were associated with overexpression of a small subset of X chromosomal genes. Other, equally aneuploid, but non-BLC rarely displayed these X chromosome abnormalities. These results suggest that X chromosome abnormalities contribute to the pathogenesis of BLC, both inherited and sporadic.

Alleles↗

PubMed2022

Expansive and Diverse Phenotypic Landscape of Field Aedes aegypti (Diptera: Culicidae) Larvae with Differential Susceptibility to Temephos: Beyond Metabolic Detoxification.

Arboviruses including dengue, Zika, and chikungunya are amongst the most significant public health concerns worldwide. Arbovirus control relies on the use of insecticides to control the vector mosquito Aedes aegypti (Linnaeus), the success of which is threatened by widespread insecticide resistance. The work presented here profiled the gene expression of Ae. aegypti larvae from field populations of Ae. aegypti with differential susceptibility to temephos originating from two Colombian urban locations, Bello and Cúcuta, previously reported to have distinctive disease incidence, socioeconomics, and climate. We demonstrated that an exclusive field-to-lab (Ae. aegypti strain New Orleans) comparison generates an over estimation of differential gene expression (DGE) and that the inclusion of a geographically relevant field control yields a more discrete, and likely, more specific set of genes. The composition of the obtained DGE profiles is varied, with commonly reported resistance associated genes including detoxifying enzymes having only a small representation. We identify cuticle biosynthesis, ion exchange homeostasis, an extensive number of long noncoding RNAs, and chromatin modelling among the differentially expressed genes in field resistant Ae. aegypti larvae. It was also shown that temephos resistant larvae undertake further gene expression responses when temporarily exposed to temephos. The results from the sampling triangulation approach here contribute a discrete DGE profiling with reduced noise that permitted the observation of a greater gene diversity, increasing the number of potential targets for the control of insecticide resistant mosquitoes and widening our knowledge base on the complex phenotypic network of the Ae. aegypti response to insecticides.

Aedes↗

PubMed2005

Assembly of severe acute respiratory syndrome coronavirus RNA packaging signal into virus-like particles is nucleocapsid dependent.

The severe acute respiratory syndrome coronavirus (SARS-CoV) was recently identified as the etiology of SARS. The virus particle consists of four structural proteins: spike (S), small envelope (E), membrane (M), and nucleocapsid (N). Recognition of a specific sequence, termed the packaging signal (PS), by a virus N protein is often the first step in the assembly of viral RNA, but the molecular mechanisms involved in the assembly of SARS-CoV RNA are not clear. In this study, Vero E6 cells were cotransfected with plasmids encoding the four structural proteins of SARS-CoV. This generated virus-like particles (VLPs) of SARS-CoV that can be partially purified on a discontinuous sucrose gradient from the culture medium. The VLPs bearing all four of the structural proteins have a density of about 1.132 g/cm(3). Western blot analysis of the culture medium from transfection experiments revealed that both E and M expressed alone could be released in sedimentable particles and that E and M proteins are likely to form VLPs when they are coexpressed. To examine the assembly of the viral genomic RNA, a plasmid representing the GFP-PS580 cDNA fragment encompassing the viral genomic RNA from nucleotides 19715 to 20294 inserted into the 3' noncoding region of the green fluorescent protein (GFP) gene was constructed and applied to the cotransfection experiments with the four structural proteins. The SARS-CoV VLPs thus produced were designated VLP(GFP-PS580). Expression of GFP was detected in Vero E6 cells infected with the VLP(GFP-PS580), indicating that GFP-PS580 RNA can be assembled into the VLPs. Nevertheless, when Vero E6 cells were infected with VLPs produced in the absence of the viral N protein, no green fluorescence was visualized. These results indicate that N protein has an essential role in the packaging of SARS-CoV RNA. A filter binding assay and competition analysis further demonstrated that the N-terminal and C-terminal regions of the SARS-CoV N protein each contain a binding activity specific to the viral RNA. Deletions that presumably disrupt the structure of the N-terminal domain diminished its RNA-binding activity. The GFP-PS-containing SARS-CoV VLPs are powerful tools for investigating the tissue tropism and pathogenesis of SARS-CoV.

Animals↗

PubMed2025

Expression regulation network in papillae of sea cucumbers: Whole-transcriptome and DNA methylation datasets.

To elucidate the expression regulation network of papilla size of sea cucumbers (Apostichopus japonicus), the whole-transcriptome and DNA methylome datasets of different sizes of papillae in sea cucumbers were generated. Average clean bases of whole-transcriptome (16.35 G) and DNA methylome (28.92 G) were obtained using RNA sequencing and whole-genome bisulfite sequencing techniques. A total of 3,188 ceRNA networks were also identified including 3,081 long non-coding RNAs (lncRNA)/microRNAs (miRNA)/mRNA networks and 107 circular RNA (circRNA)/miRNA/mRNA networks. Methylome data indicate that there were 3,307 and 3,776 differentially methylated regions (DMRs) with high-level methylation as well as 3,125 and 3,016 DMRs with low-level methylation in big papillae compared to small papillae. The identified DMRs were mainly distributed in introns, promotors, or exons. The whole-transcriptome and DNA methylome datasets generated from this study not only established a robust theoretical foundation (especially from the epigenetic aspect) for elucidating expression regulation network determining papilla size in sea cucumbers but also can be a valuable resource of biomarker mining for papilla appearance-based selective breeding in sea cucumbers.

DNA Methylation↗