Micro-ribonucleic acid expression profiling and bioinformatic target gene analyses in laryngeal carcinoma.

Abnormal expression of micro-ribonucleic acid (miRNA) might be clinically valuable as a biomarker or treatment target in the early diagnosis, treatment, and prognosis of tumors. However, little is known concerning abnormal miRNA expression of laryngeal carcinoma, one of the most commonly encountered head and neck tumors. Microarray analysis was used to obtain miRNA-expression profiles of ten pairs of freshly frozen laryngeal carcinoma tissue and surrounding normal tissue specimens. Characteristic miRNAs that were significantly related to laryngeal carcinoma were identified. Verification was performed using an additional 32 pairs of samples. The expression of two miRNAs (miR-21-3p and miR-106b-3p) was upregulated in both microarray and quantitative real-time polymerase chain-reaction analyses, whereas the expression of six miRNAs (let-7f-5p, miR-10a-5p, miR-125a-5p, miR-144-3p, miR-195-5p, and miR-203) was downregulated. The decreased expression of let-7f-5p and miR-195-5p is a novel finding in head and neck cancer. The target genes of these miRNAs were also predicted through multiple software programs. The differential expression of miRNAs might be related to the early onset and development of laryngeal carcinoma, and may be exploited as new biomarkers and therapeutic targets in the treatment of laryngeal carcinoma.

Introduction

Head and neck tumors constitute one of the main cancers endangering human health, ranking eighth in incidence among all cancer types. Laryngeal carcinoma is one of the most common types of head and neck tumors, and is the second-most prevalent of the respiratory tract tumors.1 Most laryngeal carcinoma cases are squamous cell carcinoma. The incidence of laryngeal carcinoma has been increasing in recent years. In the People’s Republic of China, the incidence of laryngeal carcinoma has increased by about 25% annually in an era of rapid industrialization.2,3 Although new surgical procedures, chemotherapeutic drugs, and more advanced radiotherapy have been applied to the treatment of laryngeal carcinoma in the past 30 years, the overall survival of patients with laryngeal carcinoma and other head and neck tumors alike has not improved, and in some cases has reduced.4 The survival rate of late-stage laryngeal carcinoma is only 30%–40%.5 Possible epidemiologically related factors for laryngeal carcinoma include smoking, drinking, air pollution, and occupational factors.6,7 The pathogenesis of laryngeal carcinoma is related to the activation of oncogenes, including BCL2 and c-MYC,8,9 and the inactivation of tumor-suppressor genes, including p53 and RB.10,11 However, the pathogenic mechanism of laryngeal carcinoma at the molecular level is still unclear, which has hindered therapeutic interventions and efforts to increase the survival rate.

Micro-ribonucleic acids (miRNAs) are small (19–25 nucleotides) noncoding RNAs.12 At the post-transcriptional level, miRNAs regulate the expression level of target genes, which influences the regulation of biological events, including individual development, apoptosis, proliferation, and differentiation.13,14 Recent studies have implicated miRNAs in the regulation of oncogenes or tumor-suppressor genes that are germane in the occurrence, development, invasion, and metastasis of a variety of tumors.15,16 Furthermore, miRNAs might be related to the pathogenesis of head and neck tumors, such as laryngeal carcinoma.17 The differential expression of miRNAs might have potential clinical value as biomarkers or treatment targets in the early diagnosis, treatment, and prognosis of laryngeal carcinoma.

Although miRNA-expression profiling of many head and neck tumors, including oral and thyroid carcinoma, has been reported,18,19 the situation is far less clear for laryngeal carcinoma, one of the most common head and neck tumors. The present study was undertaken to address this shortcoming. The aim of this study was to identify new diagnostic biomarkers for the early diagnosis of laryngeal carcinoma or potential treatment targets for laryngeal carcinoma.

Materials and methods

Tissue sampling and clinical data collection

A total of 42 pairs of surgically removed laryngeal carcinoma tissue and paracancerous tissue samples were collected from patients who received surgical resection treatment at Guangdong General Hospital from October 2007 to March 2009. Informed consent for tissue use was provided beforehand by all patients, and the study was approved by the ethics committee of Guangdong General Hospital. Laryngeal carcinoma was confirmed by pathological study. Normal laryngeal mucosa specimens were retrieved 10 mm outside the negative margin. Samples were placed in RNAlater TissueProtect Tubes (Qiagen, Hilden, Germany) and stored at −80°C. Postsurgical pathology confirmed the diagnosis of squamous cell carcinoma in all patients. The patients did not receive radio- or chemotherapy before surgery. Samples were randomly allocated into two groups for analysis using the miRNA microarray (ten pairs) or quantitative real-time polymerase chain reaction (qRT-PCR; 32 pairs).

RNA extraction

Total RNA was isolated using TRIzol® (Life Technologies, Carlsbad, CA, USA) and miRNeasy mini kit (Qiagen) according to the manufacturers’ instructions, which efficiently recovered all RNA species, including miRNAs. RNA quality and quantity was measured using an ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA), and RNA integrity was determined by gel electrophoresis.

Microarray analysis of miRNA

Total RNA of ten pairs of laryngeal carcinoma tissue samples was randomly selected, and microarray analysis was completed with the sixth generation of miRCURY LNA™ (version 16.0) by Exiqon (Vedbaek, Denmark).

RNA labeling and array hybridization

After RNA isolation, the miRCURY Hy3™/Hy5™ power labeling kit (Exiqon) was used according to the manufacturer’s guidelines for miRNA labeling. After stopping the labeling procedure, the Hy3-labeled samples were hybridized on the miRCURY LNA array according to the instructions in the array manual. Following hybridization, the slides were washed several times using a wash buffer kit (Exiqon) and dried by centrifugation for 5 minutes at 400 rpm. The slides were scanned using a GenePix 4000B microarray scanner (Molecular Devices, Sunnyvale, CA, USA).

Data analyses

Scanned images were imported into the GenePix Pro 6.0 software (Molecular Devices) for grid alignment and data extraction. Hierarchical clustering was performed using MEV (multiexperiment viewer) version 4.6 (TIGR, Rockville, MD, USA).

Validation using qRT-PCR

Total RNA was extracted from 32 pairs of laryngeal carcinoma tissue and paracancerous tissue samples. The SYBR green assays (Life Technologies) with U6 as internal control was used, and qRT-PCR was performed to validate the results.

Reverse transcription

After total RNA extraction, complementary deoxyribonucleic acid was synthesized using 2 ng to 2 μg of total RNA by Moloney Murine Leukemia Virus Reverse Transcriptase (M-MLV RT) according to the manufacturer’s instructions (Life Technologies). Stem-loop reverse-transcriptase primer was added to total RNA (synthesized by RiboBio, Guangzhou, People’s Republic of China). The volume was expanded to 24 μL with ribonuclease (RNase)-free water and mixed well. After incubation at 70°C for 10 minutes, the mixture was placed in an ice bath for 2 minutes. The following reagents were added to the mixture in sequence: 2 μL of RNase inhibitor, 8 μL of 5× M-MLV RT buffer, 2 μL of M-MLV RT reverse transcriptase (RNase H-), and 2 μL of dNTP mixture (Life Technologies) and RNase-free water (TaKaRa Biotechnology, Dalian, People’s Republic of China) to a final volume of 40 μL. The mixture was placed at 30°C for 10 minutes, followed by 42°C incubation for 1 hour, and then the mixture was placed at 70°C for 15 minutes. The product of reverse transcription was used as the template of PCR.

qPCR

Amplification was performed on a fluorescent qPCR instrument (Bio-Rad Laboratories, Hercules, CA, USA). For each reaction, three duplicates were prepared. The total reaction volume was 20 μL, which included 10 μL of 2× Mix SYBR Green I (Life Technologies) fluorescent solution, 0.25 μL each of forward and reverse primers (10 μM), 1 μL of sample template, and RNase-free water that brought the volume to 20 μL. The reaction cycle was predenaturation at 95°C for 20 seconds, followed by 40 cycles of denaturation at 95°C for 10 seconds, annealing at 60°C for 20 seconds and extension at 70°C for 10 seconds. An amplification curve was plotted. A fluorescence signal was collected during the heating process from 60°C to 95°C after 40 cycles, and a melting curve was plotted. After each reaction, the threshold cycle (Ct) was calculated. Relative quantitation of gene expression was analyzed with the 2−ΔΔCt method using the obtained Ct values.

Statistical analyses

The results of the miRNA microarray were analyzed with significance analysis of microarrays (SAM) software. The q-value representing the false discovery rate (FDR) was set at 0.05, and miRNAs showing differential expression in laryngeal carcinoma tissue and paracancerous tissue were screened and identified. The results of qRT-PCR were processed with the 2−ΔΔCt method and the difference between two groups was analyzed with the independent-samples t-test using SPSS 13.0 (SPSS, Chicago, IL, USA).

Results

Demographic and clinical data of the 42 laryngeal cancer patients

The results are summarized in Table 1.

Table 1

Demographic and clinical features of the laryngeal cancer patients

Characteristic	n (%)
Age (years), mean ± standard deviation	59.8±11.7
Sex
Male	39 (92.9)
Female	3 (7.1)
Type
Supraglottic	12 (28.6)
Glottic	28 (66.7)
Subglottic	2 (4.8)
Clinical stage
I	19 (45.2)
II	8 (19.0)
III	2 (4.8)
IV	13 (31.0)
T stage
T1	19 (45.2)
T2	10 (23.8)
T3	5 (11.9)
T4	8 (19.1)
N stage
N0	33 (78.6)
N1–3	9 (21.4)
Pathological differentiation
Poorly differentiated	3 (7.1)
Moderately differentiated	23 (54.8)
Highly differentiated	16 (38.1)

Abbreviations: n, number of patients; T, tumor size; N, lymph node status.

Results of miRNA microarray

The miRNA microarray analysis screened out 780 miRNAs. After SAM analysis with FDR at 0.05, expressions of eleven miRNAs were upregulated in tumor tissue and 114 were downregulated, compared to expression in normal mucosa tissue (Tables 2, 3 and Figure 1).

Table 2

Significantly upregulated microRNAs (SAM, FDR ≤5%)

Name	Fold change	q-value (%)
hsa-miR-21-3p	1.823432255	0
hsa-miR-205-3p	1.942721294	0
hsa-miR-106b-3p	1.213735489	0
hsa-let-7f-1-3p	1.075200118	0
hsa-miR-1246	185.9579442	1.153846154
hsa-miR-491-3p	4,835.799572	1.153846154
hsa-miR-130b-3p	1.172966978	1.978021978
hsa-miR-3182	1.816258072	2.492307692
hsa-miR-1290	157.9054397	2.492307692
hsa-miR-2115-3p	1.252350132	2.492307692
hsa-miR-720	44,171.33165	4.984615385

Abbreviations: RNA, ribonucleic acid; SAM, significance analysis of microarrays; FDR, false discovery rate.

Table 3

Significantly downregulated microRNAs (SAM, FDR ≤5%)

Name	Fold change	q-value (%)
hsa-miR-99a-5p	0.229416618	0
hsa-miR-1228-5p	0.969412922	0
hsv1-miR-H8	0.822813522	0
hsa-miR-139-5p	0.745290916	0
hsa-miR-99a-3p	0.943420998	0
hsa-miR-125b-5p	0.000504141	0
hsa-miR-208b	0.856998652	0
hsa-let-7c	0.163369241	0
hsa-miR-26a-5p	0.41283036	0
hsa-miR-101-3p	0.011547838	0
hsa-miR-26b-5p	0.044496706	0
hsa-let-7a-5p	0.115018014	0
hsa-miR-30a-5p	0.555963217	0
hsa-miR-10a-5p	0.730429293	0
hsa-miR-100-5p	0.395246953	0
hsa-miR-1297	0.675961512	0
hsa-miR-378a-5p	0.959315913	0
hsa-miR-3656	0.201641898	0
hsa-miR-337-5p	0.949658024	0
hsa-miR-375	0.536732979	0
hsa-miR-338-3p	0.651622173	0
hsa-miR-451a	3.05E-05	0
hsa-miR-3185	0.985728913	0
hsa-miR-126-3p	0.019184657	0
hsa-miR-497-5p	0.83274609	0
hsa-miR-101-5p	0.970538127	0
hsa-miR-29c-3p	0.141944203	0
hsa-miR-126-5p	0.332925548	0
hsa-miR-625-5p	0.839062834	0
hsa-miR-10b-5p	0.450653681	0
hsa-miR-30c-5p	0.226168321	0
hsa-miR-135a-5p	0.93013101	0
hsa-miR-125a-5p	0.09744404	0
hsa-miRPlus-C1070	0.732527419	0
hsa-miR-199b-5p	0.286949641	0
hsa-miR-29a-5p	0.911032396	0
hsa-miR-323b-5p	0.772414981	0
hsa-miR-29c-5p	0.862849161	1.416083916
hsa-miR-608	0.982843666	1.416083916
hsa-miR-4328	0.313675459	1.416083916
hsa-miR-129-1-3p	0.937204791	1.416083916
hsa-miR-3133	0.940025146	1.416083916
hsa-miR-23c	0.940121803	1.416083916
hsa-miR-30e-5p	0.359399404	1.416083916
hsa-let-7f-5p	0.653879295	1.416083916
hsa-miR-30e-3p	0.734655027	1.416083916
hsa-miR-144-3p	0.039732846	1.416083916
hsa-miR-98	0.728012549	1.416083916
hsa-miR-203	0.001955606	1.416083916
hsa-miR-486-5p	0.714682911	1.416083916
hsa-miR-3607-5p	0.839520455	1.416083916
hsa-miR-3195	0.367848181	1.978021978
hsa-miR-29a-3p	0.005497941	1.978021978
hsa-miR-378a-3p/hsa-miR-378c/hsa-miR-378d	0.680647548	1.978021978
hsa-miR-374a-5p	0.629791599	1.978021978
hsa-miR-30b-5p	0.195820835	1.978021978
hsa-miRPlus-C1100	0.969513061	1.978021978
hsa-miR-4321	0.973550709	1.978021978
hsa-miR-335-5p	0.832271437	1.978021978
hsa-miR-125b-1-3p	0.95000675	1.978021978
hsa-miR-378c	0.787146273	1.978021978
hsa-miR-1243	0.860868473	2.492307692
hsa-miR-191-5p	0.624357656	2.492307692
hsa-miR-20b-5p	0.810769473	2.492307692
hsa-miR-299-3p	0.910445284	2.492307692
hsa-miR-10b-3p	0.990773272	2.760467381
hsa-miR-145-5p	0.421421711	2.760467381
hsa-let-7g-5p	0.113667898	2.760467381
hsa-miR-23b-5p	0.977354147	2.760467381
mcv-miR-M1-3p	0.895527792	2.760467381
hsa-miR-29b-2-5p	0.988923772	2.760467381
hsa-miR-3176	0.934463957	2.760467381
hsa-miR-299-5p	0.554176248	2.760467381
hsa-miR-3934	0.965320293	2.760467381
hsa-miR-4317	0.570519513	2.760467381
hsa-miR-195-5p	0.434556198	2.760467381
hsa-miR-629-5p	0.921854174	2.760467381
hsa-miR-4796-3p	0.987307839	2.760467381
hsa-miR-363-3p	0.953261477	2.760467381
hsa-miR-208a	0.952193195	2.760467381
hsv1-miR-H6-5p	0.955520715	2.760467381
hsa-miR-3065-5p	0.984253829	2.760467381
hsa-miR-1287	0.98230182	2.760467381
hsa-miR-143-3p	0.023686522	3.279352227
hsa-miR-3650	0.953404765	3.279352227
hsa-miRPlus-I107*	0.928817452	3.279352227
hcmv-miR-UL148D	0.884013565	3.279352227
hsa-miR-20a-5p	0.519873511	3.279352227
hsa-miR-135b-5p	0.882995911	3.279352227
hsa-miR-200a-5p	0.982126485	3.279352227
hsa-miR-200a-3p	0.125234644	3.279352227
hsa-miR-660-5p	0.919752055	3.279352227
hsa-miR-3126-3p	0.979864468	3.279352227
hsa-miR-648	0.968124533	3.279352227
hsa-miR-3620	0.967574384	4.984615385
hsa-miR-130a-3p	0.473237407	4.984615385
hsa-miR-604	0.975573233	4.984615385
hsa-miR-138-2-3p	0.826911224	4.984615385
hsa-miRPlus-I137*	0.909617132	4.984615385
hsa-miR-2053	0.984695021	4.984615385
hsa-miR-374c-5p	0.901568699	4.984615385
hsa-miR-148a-5p	0.967409666	4.984615385
hsa-miR-659-3p	0.92450892	4.984615385
hsa-miR-27b-3p	0.176138633	4.984615385
hsa-miR-652-3p	0.944030002	4.984615385
hsa-miR-557	0.944551902	4.984615385
hsa-miR-34a-3p	0.978989503	4.984615385
hsa-miR-7-1-3p	0.992453904	4.984615385
hsa-miR-376a-5p	0.97514717	4.984615385
hsa-miR-148a-3p	0.568972471	4.984615385
hsa-miR-28-5p	0.932157662	4.984615385
hsa-miR-23b-3p	0.002949212	4.984615385
hsa-miR-551b-3p	0.976024864	4.984615385
hsa-miR-559	0.917545663	4.984615385

Abbreviations: RNA, ribonucleic acid; SAM, significance analysis of microarrays; FDR, false discovery rate.

Figure 1

Cluster-analysis diagram of laryngeal carcinoma micro-ribonucleic acid (miRNA) genes. The heat map of hierarchical clustering shows an miRNA gene differential expression profiling between paired laryngeal tumor and surrounding normal tissue samples. Red represents upregulated miRNAs, and green represents downregulated mirnas.

Abbreviations: lt, laryngeal tumor; lnt, laryngeal normal tissue.

qRT-PCR validation

Eight miRNAs (miR-21-3p, miR-106b-3p, let-7f-5p, miR-10a-5p, miR-125a-5p, miR-144-3p, miR-195-5p, and miR-203) were validated with qRT-PCR. The expression difference of these eight miRNAs was statistically significant between laryngeal carcinoma and normal mucosa tissue in the 32 pairs of specimens (P<0.05). Further comparison between the microarray and qRT-PCR results (Figure 2) showed that the expression pattern of these eight miRNAs in microarray was consistent with that in qRT-PCR.

Figure 2

Comparison of quantitative real-time polymerase chain-reaction (qRT-PCR) results and microarray results of eight micro-ribonucleic acids. The open bars represent qRT-PCR analysis, while the solid bars represent the microarray data.

Abbreviation: RNA, ribonucleic acid.

Prediction of target genes

Potential target genes of the eight miRNAs were preliminarily predicted using bioinformatics software, including PicTar, miRanda, and TargetScan (Table 4 and Figure 3).

Table 4

Putative target genes of differentially expressed micro-ribonucleic acids

MicroRNA	Chromosome location	Putative target
miR-21-3p	17q23.1	LUM, CHD7, NFAT5, WWP1, TGFBI
miR-106b-3p	7q22.1	TNFRSF21, ACPL2, EPHA5, RAB22A, SLITRK3
let-7f-5p	9q22.32	ABCB9, CCNJ, RDH10, GATM, HIC2
miR-10a-5p	17q21.32	TNRC6B, BCL6, WNK3, RAP2A, GATA6
miR-125a-5p	19q13.41	SMARCD2, BMF, ERBB4, MTF1, ZFYVE1
miR-144-3p	17q11.2	ZNF800, ACBD3, TRIO, ETS1, FBXL3
miR-195-5p	17p13.1	TLE4, ST8SIA3, LRP6, WNT3A, KIF21A
miR-203	14q32.33	MAP3K13, ZNF281, FAT3, KIAA1429, ADAMTS6

Abbreviation: RNA, ribonucleic acid.

Figure 3

Target-gene prediction of selected micro-ribonucleic acids (miRNAs) through bioinformatic techniques. Red-box nodes represent miRNA, and green-cycle nodes represent putative target gene. Edges describe the inhibitory effects of miRNA on potential target gene.

Discussion

In this study, the latest version of microarray mercury LNA (version 16.0) was used to analyze ten pairs of freshly frozen laryngeal carcinoma tissue and paracancerous tissue samples, and the miRNA-expression profile of laryngeal carcinoma was obtained. The results were verified by qRT-PCR using another 32 pairs of samples. The results between the two analyses were comparable. miR-21-3p and miR-106b-3p were upregulated, and let-7f-5p, miR-10a-5p, miR-125a-5p, miR-144-3p, miR-195-5p, and miR-203 were downregulated in the microarray and qRT-PCR analyses. The target genes of these miRNAs were also predicted through multiple software programs in this study.

Studies of miRNA-expression profiling of laryngeal carcinoma have suffered from some shortcomings. Most studies were conducted on head and neck squamous cell carcinoma,17,20,21 and few studies focused simply on laryngeal carcinoma.22 Although the disease locations of all kinds of tumors in the head and neck area are in close proximity, the larynx is part of the respiratory system, and so is different from such components of digestive system as the oral and pharyngeal areas.23,24 Laryngeal carcinoma has significantly distinctive pathological characteristics and prognosis from oral cancer, oropharyngeal cancer, and hypopharyngeal carcinoma. For example, the incidence of laryngeal carcinoma in men is significantly higher than that of other head and neck tumors.24 In addition, a comparative genomic hybridization study showed that the chromosomal patterns of laryngeal carcinoma and other head and neck squamous cell carcinomas were significantly different.25 Thus, without distinguishing laryngeal carcinoma from other head and neck squamous cell carcinoma, the results will be unavoidably biased.f

In this study, the differentially expressed miRNAs detected in laryngeal carcinoma were different from reports of other head and neck tumors or laryngeal carcinoma. The numbers of differentially expressed miRNAs were larger, and were predominantly accounted for by down-regulated miRNAs. The sample size of laryngeal carcinoma samples in this study was the largest among studies on head and neck tumors reported to date. Due to the heterogeneity of tumor samples, a large sample size might be beneficial for the accurate detection of abnormally expressed miRNA.26 In addition, the present study utilized the latest version of microarray chip with locked nucleic acid (LNA) technology,27 which guaranteed the maximal degree of reliability of the results. This is the first report on data from southern People’s Republic of China, and the different geographic locales of patients might also have contributed to the difference of the results. Furthermore, some studies only performed microarray detection without further verification by qRT-PCR, which has been recognized as a gold standard. This might also affect study results.

Among the eight differentially expressed miRNAs in laryngeal carcinoma verified by qRT-PCR, miR-21-3p,22 miR-106b-3p,20 miR-10a-5p,20 miR-125a-5p,17 miR-203,28 and miR-144-3p29 have been reported in other head and neck cancer studies. miR-125a-5p, miR-203, and miR-144-3p, which were only reported in other studies utilizing microarrays, were further verified by qRT-PCR in our study.

In addition, miRNAs not previously reported or inconsistent with previous head and neck squamous cell carcinoma or laryngeal carcinoma studies were also identified, and included let-7f-5p and miR-195-5p. Lethal-7 (let-7) was initially discovered in Caenorhabditis elegans. Being one of the first tumor-suppressive miRNAs identified, let-7 is downregulated in lung, ovarian, and breast cancer. Studies showed that let-7i is highly expressed17 in head and neck squamous cell carcinoma, and the expressions of let-7a, c, and e are low.20 In our study, let-7f-5p showed low expression, which was opposite to a report from Tran et al.30 The samples used in the prior study were from the cell lines of different head and neck squamous cell carcinomas, including tongue, tonsil, hypopharynx, and larynx, and no further verification of let-7f was conducted.30 Our study used freshly frozen laryngeal carcinoma tissue as samples, and qRT-PCR verification was performed. This further confirmed the advantage of singling out such diseases as laryngeal carcinoma, instead of studying a mixture of different types of head and neck tumors. miR-195 belongs to the miR-15/miR-16 family. miR-195 expression was downregulated in some tumor types, including hepatocellular carcinoma and breast cancer,31 suggesting that MIR195 might be an important tumor-suppressor gene. Also, it has been reported to be involved in the cell cycle, apoptosis, and proliferation through regulating such target genes as WEE1, CDK6, and BCL2.31 Our study is the first to demonstrate the possible presence of abnormally decreased expression of miR-195 in laryngeal carcinoma.

The molecular mechanism of miRNA in the regulation of laryngeal carcinoma-related genes is still unclear. Through such bioinformatic tools as PicTar, miRanda, and TargetScan,32,33 target genes of these miRNAs were preliminarily screened. The results showed that the target genes are those closely related to proliferation, cell-cycle regulation, and signal transduction, which are involved in aspects of tumorigenesis, tumor development, and metastasis. Bioinformatic computational methods for predicting target genes are necessary to guide experimental validation.

One limitation of this study was the lack of cell-function experiments and target-gene function validation. Whether the abnormally expressed miRNAs identified by microarray and verified by qRT-PCR are involved in the pathogenesis of laryngeal carcinoma requires further verification by function experiments.

In summary, miRNA-expression profiling of laryngeal carcinoma was performed using microarray analysis, and qRT-PCR verification was conducted. The decreased expression of let-7f-5p and miR-195-5p is a novel finding in head and neck cancer. In addition, target genes of identified miRNA were predicted using multiple databases. The identified miRNAs may have potential value as novel biomarkers and treatment targets of laryngeal carcinoma.