Paphia textile is an important, aquaculture bivalve clam species distributed mainly in China, Philippines, and Malaysia. Recent studies of P. textile have focused mainly on artificial breeding and nutrition analysis, and the transcriptome and genome of P. textile have rarely been reported. In this work, the transcriptome of P. textile foot tissue was sequenced on an Illumina HiSeq™ 2000 platform. A total of 20,219,795 reads were generated, resulting in 4.08 Gb of raw data. The raw reads were cleaned and assembled into 54,852 unigenes with an N50 length of 829 bp. Of these unigenes, 38.92% were successfully annotated based on their matches to sequences in seven public databases. Among the annotated unigenes, 14,571 were assigned Gene Ontology terms, 5448 were classified to Clusters of Orthologous Groups categories, and 6738 were mapped to 228 pathways in the Kyoto Encyclopedia of Genes and Genomes database. For functional marker development, 5605 candidate simple sequence repeats were identified in the transcriptome and 80 primer pairs were selected randomly and amplified in a wild population of P. textile. A total of 36 loci that exhibited obvious repeat length polymorphisms were detected. The transcriptomic data and microsatellite markers will provide valuable resources for future functional gene analyses, genetic map construction, and quantitative trait loci mapping in P. textile.
Pax3 and Pax7 are the regulators and markers of muscle progenitors and satellite cells that contribute to the embryonic development and postembryonic growth of skeletal muscle in vertebrates, as well as to its repair and regeneration. However, information regarding them in vertebrate genome model, torafugu Takifugu rubripes, has remained unknown. Therefore, as an initial step, here we characterized Pax3 and Pax7 from torafugu and investigated their expression patterns during different developmental stages by RT-PCR. In silico analysis with the Fugu genome database (ver. 4.0) yielded two distinct genes each for Pax3 (Pax3a and Pax3b) and Pax7 (Pax7a and Pax7b). The 75th amino acid, glutamine (Gln75), from the N-terminus was replaced by proline in the paired box domain (PD) of Pax3a. One single cDNA clone encoding Pax3a had deletion of Gln75 in PD, suggesting the presence of alternatively spliced variants (Q+/Q-). This was further supported by identification of two adjacent alternative 3' splice acceptor sites which produce Pax3b Q+ (aagCAGGGA) and Q- (aagcagGGA) variants. Interestingly, torafugu Pax7a, but not Pax7b, had an insert encoding five amino acid residues (SGEAS) in a C-terminal region of PD in two out of three cDNA clones. Genomic analysis showed two alternate splice donor sites at exon 4 of Pax7a. In synteny analysis, torafugu Pax3a showed syntenic relationship with the corresponding regions in other teleosts only, whereas Pax3b and Pax7b showed high syntenic relationship with the corresponding regions of both mammals and other teleosts. RT-PCR revealed that expression of Pax3a and Pax3b transcripts was restricted to embryonic stages only, whereas those of Pax7a and Pax7b was continued to be expressed in larvae and importantly those of Pax7a were found in adult skeletal muscles. Therefore, Pax3 appears to be most important for primary myogenesis and Pax7 for secondary myogenesis and growth by hyperplasia in fish. In this regard, the transcripts of torafugu Pax3 and Pax7 genes might be used for further investigation as a marker for identification of muscle precursor cells during different phases of growth, and this ambiguity is the next target of our research.
This article explores the genetic history of the various sub-populations currently living in Peninsular Malaysia. This region has received multiple waves of migrants like the Orang Asli in prehistoric times and the Chinese, Indians, Europeans and Arabs during historic times. There are three highly distinct lineages that make up the Orang Asli; Semang, Senoi and Proto-Malays. The Semang, who have 'Negrito' characteristics, represent the first human settlers in Peninsular Malaysia arriving from about 50,000ya. The Senoi later migrated from Indochina and are a mix between an Asian Neolithic population and the Semang. These Asian genomes probably came in before Austroasiatic languages arrived between 5000 and 4000years ago. Semang and Senoi both now speak Austro-Asiatic languages indicative of cultural diffusion from Senoi to Semang. In contrast, the Proto-Malays who came last to the southern part of this region speak Austronesian language and are Austronesians with some Negrito admixture. It is from this group that the contemporary Malays emerged. Here we provide an overview of the best available genetic evidences (single nucleotide polymorphisms, mitochondrial DNA, Y-chromosome, blood groups, human platelet antigen, human leukocyte antigen, human neutrophil antigen and killer-cell immunoglobulin-like receptor) supporting the complex genetic history of Peninsular Malaysia. Large scale sampling and high throughput genetic screening programmes such as those using genome-wide single nucleotide polymorphism analyses have provided insights into various ancestral and admixture genetic fractions in this region. Given the now extensive admixture present in the contemporary descendants of ancient sub-populations in Peninsular Malaysia, improved reconstruction of human migration history in this region will require new evidence from ancient DNA in well-preserved skeletons. All other aspects of the highly diverse and complex genetic makeup in Peninsular Malaysia should be considered carefully for genetic mapping of disease loci and policy formation by health authorities.
Trichohepatoenteric syndrome (THES) is a rare autosomal recessive disorder that is classically associated with intractable diarrhea with an onset within the first few months of life. Herein, we investigated and reported novel mutations in two causal genes in 3 Malaysian cases. Genomic DNA was extracted from peripheral blood obtained from patients in two Malaysian Chinese families. The exons of SKIV2L and TTC37 genes were amplified and sequenced by bi-directional sequencing to identify the point mutations within the coding sequence. Three Chinese boys from two families with characteristic features and clinical course were diagnosed with THES. In family-1, two point mutations were identified in the SKIV2L gene (c.1891G>A and c.3187C>T). In family-2, a single-nucleotide duplication (c.3426dupA) was found in the TTC37 gene. These mutations cause the production of abnormal non-functional gene product leading to the clinical manifestations in the patients. We reported three point mutations, which have not been previously described in other patients with THES in SKIV2L and TTC37 genes, including one nonsense, one frameshift, and one missense mutations.
Alternative pre-mRNA splicing provides a source of vast protein diversity by removing non-coding sequences (introns) and accurately linking different exonic regions in the correct reading frame. The regulation of alternative splicing is essential for various cellular functions in both pathological and physiological conditions. In eukaryotic cells, this process is commonly used to increase proteomic diversity and to control gene expression either co- or post-transcriptionally. Alternative splicing occurs within a megadalton-sized, multi-component machine consisting of RNA and proteins; during the splicing process, this complex undergoes dynamic changes via RNA-RNA, protein-protein and RNA-protein interactions. Co-transcriptional splicing functionally integrates the transcriptional machinery, thereby enabling the two processes to influence one another, whereas post-transcriptional splicing facilitates the coupling of RNA splicing with post-splicing events. This review addresses the structural aspects of spliceosomes and the mechanistic implications of their stepwise assembly on the regulation of pre-mRNA splicing. Moreover, the role of phosphorylation-based, signal-induced changes in the regulation of the splicing process is demonstrated.
DAPK3 belongs to family of DAPK (death-associated protein kinases) and is involved in the regulation of progression of the cell cycle, cell proliferation, apoptosis and autophagy. It is considered as a tumor suppressor kinase, suggesting the loss of its function in case of certain specific mutations. The T112M, D161N and P216S mutations in DAPK3 have been observed in cancer patients. These DAPK3 mutants have been associated with very low kinase activity, which results in the cellular progression towards cancer. However, a clear understanding of the structural and biophysical variations that occur in DAPK3 with these mutations, resulting in the decreased kinase activity has yet not been deciphered. We performed a molecular dynamic simulation study to investigate such structural variations. Our results revealed that mutations caused a significant structural variation in DAPK3, majorly concentrated in the flexible loops that form part of the ATP binding pocket. Interestingly, D161N and P216S mutations collapsed the ATP binding pocket through flexible loops invasion, hindering ATP binding which resulted in very low kinase activity. On the contrary, T112M mutant DAPK3 reduces ATP binding potential through outward distortion of flexible loops. In addition, the mutant lacked characteristic features of the active protein kinase including proper interaction between HR/FD and DFG motifs, well structured hydrophobic spine and Lys42-Glu64 salt bridge interaction. These observations could possibly explain the underlying mechanism associated with the loss of kinase activity with T112M, D161N and P216S mutation in DAPK3.
Tomistoma schlegelii, also referred to as the "false gharial", is one of the most exclusive and least known of the world's fresh water crocodilians, limited to Southeast Asia. Indeed, lack of economic value for its skin has led to neglect the biodiversity of the species. The current study aimed to investigate the mentioned case using 40 simple sequence repeat (SSR) primer pairs and 45 inter-simple sequence repeat (ISSR) primers. DNA analysis of 17 T. schlegelii samples using the SSR and ISSR markers resulted in producing a total of 49 and 108 polymorphic bands, respectively. Furthermore, the SSR- and ISSR-based cluster analyses both generated two main clusters. However, the SSR based results were found to be more in line with the geographical distributions of the crocodile samples collected across the country as compared with the ISSR-based results. The observed heterozygosity (HO) and expected heterozygosity (HE) of the polymorphic SSRs ranged between 0.588-1 and 0.470-0.891, respectively. The present results suggest that the Malaysian T. schlegelii populations had originated from a core population of crocodiles. In cooperation with the SSR markers, the ISSRs showed high potential for studying the genetic variation of T. schlegelii, and these markers are suitable to be employed in conservation genetic programs of this endangered species. Both SSR- and ISSR-based STRUCTURE analyses suggested that all the individuals of T. schlegelii are genetically similar with each other.
MicroRNAs (miRNAs) are a class of small, endogenous non-coding RNAs that negatively regulate gene expression, resulting in the silencing of target mRNA transcripts through mRNA cleavage or translational inhibition. MiRNAs play significant roles in various biological and physiological processes in plants. However, the miRNA-mediated gene regulatory network in pineapple, the model tropical non-climacteric fruit, remains largely unexplored. Here, we report a complete list of pineapple mature miRNAs obtained from high-throughput small RNA sequencing and precursor miRNAs (pre-miRNAs) obtained from ESTs. Two small RNA libraries were constructed from pineapple fruits and leaves, respectively, using Illumina's Solexa technology. Sequence similarity analysis using miRBase revealed 579,179 reads homologous to 153 miRNAs from 41 miRNA families. In addition, a pineapple fruit transcriptome library consisting of approximately 30,000 EST contigs constructed using Solexa sequencing was used for the discovery of pre-miRNAs. In all, four pre-miRNAs were identified (MIR156, MIR399, MIR444 and MIR2673). Furthermore, the same pineapple transcriptome was used to dissect the function of the miRNAs in pineapple by predicting their putative targets in conjunction with their regulatory networks. In total, 23 metabolic pathways were found to be regulated by miRNAs in pineapple. The use of high-throughput sequencing in pineapples to unveil the presence of miRNAs and their regulatory pathways provides insight into the repertoire of miRNA regulation used exclusively in this non-climacteric model plant.
Glycogenes regulate a large number of biological processes such as cancer and development. In this work, we created an interaction network of 923 glycogenes to detect potential hubs from different mouse tissues using RNA-Seq data. DAVID functional cluster analysis revealed enrichment of immune response, glycoprotein and cholesterol metabolic processes. We also explored nsSNPs that may modify the expression and function of identified hubs using computational methods. We observe that the number of nsSNPs predicted by any two methods to affect protein function is 4, 7 and 2 for FLT1, NID2 and TNFRSF1B. Residues in the native and mutant proteins were analyzed for solvent accessibility and secondary structure change. Analysis of hubs can help in determining their degree of conservation and understanding their functions in biological processes. The nsSNPs proposed in this work may be further targeted through experimental methods for understanding structural and functional relationships of hub mutants.
Constitutive androstane receptor (CAR) encoded by the nuclear receptor subfamily 1, group I, member 3 (NR1I3) gene regulates the elimination of bilirubin through activating the components of the bilirubin clearance pathway. Hence, NR1I3 genetic variants may affect bilirubin metabolism and result in neonatal hyperbilirubinemia. Thus far, research which investigates the association between NR1I3 variants and neonatal hyperbilirubinemia has not been undertaken in any population. The present study aimed to evaluate the influence of MPJ6_1I3008 (rs10157822), IVS8+116T>G (rs4073054) and 540A>G (rs2307424) on neonatal hyperbilirubinemia development in the Malay population. Buccal swabs were collected from 232 hyperbilirubinemia and 277 control term newborns with gestational age ≥37weeks and birth weight ≥2500g. The NR1I3 variants were genotyped by using high resolution melting (HRM) assays and verified by DNA sequencing. Gender, mode of delivery and birth weight did not differ between hyperbilirubinemia and control groups. The genotypic and allelic frequencies of MPJ6_1I3008, IVS8+116T>G and 540A>G were not significantly different between the groups. However, stratification by gender revealed a significant inverse association between homozygous variant genotype of MPJ6_1I3008 and risk of neonatal hyperbilirubinemia in the females (OR, 0.44; 95% CI, 0.20-0.95; p=0.034). This study demonstrates that the homozygous variant genotype of MPJ6_1I3008 was associated with a significant reduced risk of neonatal hyperbilirubinemia in the females.
Three new open reading frames were found downstream from cbm71, a toxin gene from Clostridium bifermentans malaysia (Cbm) strain CH18. The first one (91bp downstream) called cbm72, is 1857bp long and encodes a 71727-Da protein (Cbm72) with a sequence similar to that of Bacillus thuringiensis delta-endotoxins. This protein shows no significant toxicity to mosquito larvae. The two others, cbm17.1 (462bp) and cbm17.2 (459bp), are copies of the same gene encoding Cbm P18 and P16 polypeptides and located 426bp and 1022bp downstream from cbm72, respectively. They encode 17189-Da and 17451-Da proteins with sequences 44.6% similar to that of Aspergillus fumigatus hemolysin; however, they were not hemolytic in the conditions tested.
Plants maintain extensive growth flexibility under different environmental conditions, allowing them to continuously and rapidly adapt to alterations in their environment. A large portion of many plant genomes consists of transposable elements (TEs) that create new genetic variations within plant species. Different types of mutations may be created by TEs in plants. Many TEs can avoid the host's defense mechanisms and survive alterations in transposition activity, internal sequence and target site. Thus, plant genomes are expected to utilize a variety of mechanisms to tolerate TEs that are near or within genes. TEs affect the expression of not only nearby genes but also unlinked inserted genes. TEs can create new promoters, leading to novel expression patterns or alternative coding regions to generate alternate transcripts in plant species. TEs can also provide novel cis-acting regulatory elements that act as enhancers or inserts within original enhancers that are required for transcription. Thus, the regulation of plant gene expression is strongly managed by the insertion of TEs into nearby genes. TEs can also lead to chromatin modifications and thereby affect gene expression in plants. TEs are able to generate new genes and modify existing gene structures by duplicating, mobilizing and recombining gene fragments. They can also facilitate cellular functions by sharing their transposase-coding regions. Hence, TE insertions can not only act as simple mutagens but can also alter the elementary functions of the plant genome. Here, we review recent discoveries concerning the contribution of TEs to gene expression in plant genomes and discuss the different mechanisms by which TEs can affect plant gene expression and reduce host defense mechanisms.
Atmospheric CO2 level is one of the most important factors which affect plant growth and crop production. Although many crucial genes and pathways have been identified in response to atmospheric CO2 changes, the integrated and precise mechanisms of plant CO2 response are not well understood. Alternative splicing (AS) is an important gene regulation process that affects many biological processes in plants. However, the AS pattern changes in plants in response to elevated CO2 levels have not yet been investigated. Here, we used RNA-Seq data of Arabidopsis thaliana grown under different CO2 concentration to analyze the global changes in AS. We found that AS increased with the rise in CO2 concentration. Additionally, we identified 345 differentially expressed (DE) genes and 251 differentially alternative splicing (DAS) genes under the elevated CO2 condition. Moreover, the results showed that the expression of most of the DAS genes did not change significantly, indicating that AS can serve as an independent mechanism for gene regulation in response to elevated CO2. Furthermore, our analysis of function categories revealed that the DAS genes were associated mainly with the stimulus response. Overall, this the first study to explore the changes of AS in plants in response to elevated CO2.
Citrus maxima "seedless" is originally from Malaysia, and now is widely cultivated in Hainan province, China. The essential features of this cultivar are thin skin, green epicarp and seedless at the ripening stage. Here, using C. maxima "seedless" as experimental material, we investigated the physical and inclusion indicators, and found the accumulation of storage compounds during 120-210 DAF leading to inconsistent increase between volume and weight. Component analysis of soluble sugar indicated that arabinose and xylose have a high content in early development of pummelo juice sacs (PJS), whereas fructose, glucose and sucrose show a significant increase during PJS maturation. To clarify a global overview of the gene expressing profiles, the PJSs from four periods (60, 120, 180 and 240 DAF) were selected for comparative transcriptome analysis. The resulting 8275 unigenes showed differential expression during PJS development. Also, the stability of 11 housekeeping genes were evaluated by geNorm method, resulting in a set of five genes (UBC, ACT, OR23, DWA2 and CYP21D) used as control for normalization of gene expression. Based on transcriptome data, 5 sucrose synthases (SUSs) and 10 invertases (INVs) were identified to be involved in sucrose degradation. Importantly, SUS4 may be responsible for arabinose and xylose biosynthesis to form the cell wall in early development, while SUS3 and VIN2 may be important in the accumulation of soluble hexose leading to cell expansion through an osmotic-independent pathway in late development. The information provides valuable metabolite and genetic resources in C. maxima "seedless", and is important for achieving high fruit yield and quality.
Short-lived therapeutic gene expression in mammalian cells by DNA methylation is one of the major challenges in gene therapy. In this study, we assessed the implication of DNA methylation on the duration of GFP expression in mouse embryonic stem (ES) and mouse induced pluripotent stem (iPS) cells. The cells were transduced with lentivirus (LV) carrying green fluorescent protein (GFP) driven by either human elongation factor (EF1α) or cytomegalovirus (CMV) promoter. Transduced iPS cells exhibited higher percentage of GFP+ cells with persistent mean fluorescent intensity than transduced ES cells. Analysis on the integrated copy of transgene in the population of the transduced cells demonstrated similar copy number. However, significant increase in GFP intensity following 5-azaC treatment was observed in transduced ES cells only, suggesting the influence of DNA methylation in transgene silencing. Subsequent DNA methylation analysis showed that the promoter and the GFP region of the provirus in iPS cells had negligible methylation profile compared to transduced ES cells. Interestingly, sustained transgene expression was observed upon directed differentiation of transduced iPS cells towards CD34+ CD45+ cells. Hence, this study has shown that favourable transgene activity from lentiviral transduced iPS cells was due to the lack of methylation at the proviral regions.
The true mahseer (Tor spp.) is one of the highest valued fish in the world due to its high nutritional value and great unique taste. Nevertheless, its morphological characterization and single mitochondrial gene phylogeny in the past had yet to resolve the ambiguity in its taxonomical classification. In this study, we sequenced and assembled 11 complete mahseer mitogenomes collected from Java of Indonesia, Pahang and Terengganu of Peninsular Malaysia as well as Sarawak of East Malaysia. The mitogenome evolutionary relationships among closely related Tor spp. samples were investigated based on maximum likelihood phylogenetic tree construction. Compared to the commonly used COX1 gene fragment, the complete COX1, Cytb, ND2, ND4 and ND5 genes appear to be better phylogenetic markers for genetic differentiation at the population level. In addition, a total of six population-specific mitolineage haplotypes were identified among the mahseer samples analyzed, which this offers hints towards its taxonomical landscape.
Rubber materials have greatly contributed to human civilization. However, being a polymeric material does not decompose easily, it has caused huge environmental problems. On the other hand, only few bacteria are known to degrade rubber, with studies pertaining them being intensively focusing on the mechanism involved in microbial rubber degradation. The Streptomyces sp. strain CFMR 7, which was previously confirmed to possess rubber-degrading ability, was subjected to whole genome sequencing using the single molecule sequencing technology of the PacBio® RS II system. The genome was further analyzed and compared with previously reported rubber-degrading bacteria in order to identify the potential genes involved in rubber degradation. This led to the interesting discovery of three homologues of latex-clearing protein (Lcp) on the chromosome of this strain, which are probably responsible for rubber degrading activities. Genes encoding oxidoreductase α-subunit (oxiA) and oxidoreductase β-subunit (oxiB) were also found downstream of two lcp genes which are located adjacent to each other. In silico analysis reveals genes that have been identified to be involved in the microbial degradation of rubber in the Streptomyces sp. strain CFMR 7. This is the first whole genome sequence of a clear-zone-forming natural rubber- degrading Streptomyces sp., which harbours three Lcp homologous genes with the presence of oxiA and oxiB genes compared to the previously reported Gordonia polyisoprenivorans strain VH2 (with two Lcp homologous genes) and Nocardia nova SH22a (with only one Lcp gene).
In the phylum of Proteobacteria, quorum sensing (QS) system is widely driven by synthesis and response of N-acyl homoserine lactone (AHL) signalling molecules. AHL is synthesized by LuxI homologue and sensed by LuxR homologue. Once the AHL concentration achieves a threshold level, it triggers the regulation of target genes. In this study, QS activity of Citrobacter amalonaticus strain YG6 which was isolated from clams was investigated. In order to characterise luxI/R homologues, the genome of C. amalonaticus strain YG6 (4.95 Mbp in size) was sequenced using Illumina MiSeq sequencer. Through in silico analysis, a pair of canonical luxI/R homologues and an orphan luxR homologue were identified and designated as camI, camR, and camR2, respectively. A putative lux box was identified at the upstream of camI. The camI gene was cloned and overexpressed in E. coli BL21 (DE3)pLysS. High-resolution triple quadrupole liquid chromatography mass spectrometry (LC-MS/MS) analysis verified that the CamI is a functional AHL synthase which produced multiple AHL species, namely N‑butyryl‑l‑homoserine lactone (C4-HSL), N‑hexanoyl‑l‑homoserine lactone (C6-HSL), N‑octanoyl‑l‑homoserine lactone (C8-HSL), N‑tetradecanoyl‑l‑homoserine lactone (C14-HSL) and N‑hexadecanoyl‑l‑homoserine lactone (C16-HSL) in C. amalonaticus strain YG6 and camI gene in recombinant E. coli BL21(DE3)pLysS. To our best knowledge, this is the first functional study report of camI as well as the first report describing the production of C14-HSL by C. amalonaticus.
Casuarina equisetifolia L. is an important multi-purpose, fast growing and widely planted tree species native to tropical and subtropical coastlines of Australia, Southeast Asia, Malaysia, Melanesia, Polynesia and New Caledonia. It is a nitrogen-fixing tree mainly used for charcoal making, construction poles, landscaping, timber, pulp, firewood, windbreaks, shelterbelts, soil erosion and sand dune stabilization. Casuarina wood is presently used for paper and pulp production. Raw material with reduced lignin is highly preferred to increase the pulp yield. Hence, understanding the molecular regulation of wood formation in this tree species is vital for selecting industrially suitable phenotypes for breeding programs. The lignin biosynthetic pathway has been extensively studied in tree species like Eucalypts, poplars, pines, Picea, Betula and Acacia sp. However, studies on wood formation at molecular level is presently lacking in casuarinas. Hence, in the present study, the transcriptome of the developing secondary tissues of 15 years old Casuarina equiseitfolia subsp. equisetifolia was sequenced, de novo assembled, annotated and mapped to functional pathways. Transcriptome sequencing generated a total of 26,985 transcripts mapped to 31 pathways. Mining of the annotated data identified nine genes involved in lignin biosynthesis pathway and relative expression of the transcripts in four tissues including scale-like leaves, needle-like brachlets, wood and root were documented. The expression of CeCCR1 and CeF5H were found to be significantly high in wood tissues, while maximum expression of CeHCT was documented in stem. Additionally, CeTUBA and CeH2A were identified as the most stable reference transcript for normalization of qRT-PCR data in C. equisetifolia. The present study is the first wood genomic resource in C. equisetifolia, which will be valuable for functional genomics research in this genus.
MicroRNA-3099 is highly expressed during neuronal differentiation and development of the central nervous system. Here we characterised the role of miR-3099 during neural differentiation and embryonic brain development using a stable and regulatable mouse embryonic stem cell culture system for miR-3099 expression and in utero electroporation of miR-3099 expression construct into E15.5 embryonic mouse brains. In the in vitro system, miR-3099 overexpression upregulated gene related to neuronal markers such as Tuj1, NeuN, Gat1, vGluT1 and vGluT2. In contrast, gene related to astrocyte markers (Gfap, S100β and Slc1a3) were suppressed upon overexpression of miR-3099. Furthermore, miR-3099 overexpression between E15.5 and E18.5 mouse embryonic brains led to disorganised neuronal migration potentially due to significantly decreased Gfap+ cells. Collectively, our results indicated that miR-3099 plays a role in modulating and regulating expression of key markers involved in neuronal differentiation. In silico analysis was also performed to identify miR-3099 homologues in the human genome, and candidates were validated by stem-loop RT-qPCR. Analysis of the miR-3099 seed sequence AGGCUA against human transcriptomes revealed that a potential miRNA, mds21 (Chr21:39186698-39186677) (GenBank accession ID: MK521584), was 100% identical to the miR-3099 seed sequence. Mds21 expression was observed and validated in various human cell lines (293FT, human Wharton's jelly and dental pulp mesenchymal stem cells, and MCF-7, MDA-MB-231, C-Sert, SW780, RT112, 5637, EJ28 and SH-SY5Y cells), with the highest levels detected in human mesenchymal stem cell lines. The analysis validated mds21 as a novel miRNA and a novel homologue of miR-3099 in the human genome.