- Department of Human Anatomy, Histology and Embryology, School of Basic Medicine, Qingdao University, Qingdao, China
Purpose: This work aimed to identify differentially expressed circular RNAs (circRNAs) and elucidate their potential function in aquaporin 5 (AQP5) knockout (AQP5–/–) mice with the primary dry eye phenotype.
Methods: A slit lamp examination was performed on AQP5–/– mice to assess corneal epithelial defects using fluorescein sodium staining. Hematoxylin–eosin staining and transmission electron microscopy analysis were performed to identify structural changes in lacrimal gland epithelial cells due to AQP5 deficiency. The expression profiles of circRNA and messenger RNA (mRNA) were determined by a microarray analysis. The selected circRNA was verified by quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed to predict the biological functions and the potential pathways of parental genes involved in lacrimal gland epithelial cell changes. According to the bioinformatics analysis of identified circRNAs, we predicted a circRNA–miRNA–mRNA network of phagosomes.
Results: The AQP5–/– mice spontaneously exhibit dry eye symptoms, wherein the AQP5 deficiency changes the structure of lacrimal gland epithelial cells. The analysis revealed that, compared to AQP5+/+ mice, 30 circRNAs in the lacrimal glands of AQP5–/– mice were differentially expressed (fold change ≥ 2.0, p < 0.05). Nine upregulated circRNAs were identified using qRT-PCR, and nine upregulated validated circRNAs, 40 altered microRNAs (miRNAs), and nine upregulated mRNAs were identified through a network analysis. The KEGG analysis showed that these nine target genes were expressed in phagosomes.
Conclusion: The AQP5–/– mice have primary and stable dry eye phenotypes from birth. We identified differently expressed circRNAs in the lacrimal glands of AQP5–/– and AQP5+/+ mice, predicting a circRNA–miRNA–mRNA network of phagosomes. CircRNA likely plays an important role in lacrimal gland epithelial cell pathogenesis. Therefore, it is reasonable to use circRNA as a potential therapeutic agent for the treatment of dry eyes.
Introduction
Dry eye is currently the most common eye disease except for ametropia, with a worldwide incidence of 5–34% (Messmer, 2015). The main pathophysiological mechanisms are tear film instability, increased tear permeability, ocular surface inflammation and injury, and neurosensory abnormalities (Craig et al., 2017). The current treatment strategies for dry eye include the use of artificial tears, administration of anti-inflammatory drugs, and surgery (Lemp and Foulks, 2007). However, the incidence of dry eye is still rising, emphasizing the urgent need to determine its pathogenesis and to develop effective treatment measures.
Aquaporin 5 (AQP5), a kind of aquaporin, is highly expressed in the corneal epithelium and lacrimal glands and is closely related to eye diseases (Raina et al., 1995; Takata et al., 2004). Compared with AQP5+/+ mice, the lenses of AQP5–/– mice in vitro appear to be slightly turbid at high glucose concentrations, suggesting that AQP5 might maintain crystal transparency by regulating osmotic pressure (Sindhu Kumari and Varadaraj, 2013). AQP5 has also been shown to reduce saliva and airway mucus secretions and increase corneal thickness in experimental animals (Ma et al., 1999; Thiagarajah and Verkman, 2002). In addition, lack of AQP5 affects the migration and the proliferation of cells, leading to slow corneal healing (Kumari et al., 2018). One of the symptoms of dry eye is an epithelial defect in the cornea (Levin and Verkman, 2006). A previous study considered AQP5 level as an index of dry eye (Lin et al., 2019), and the occurrence of dry eye is believed to be accompanied by a decrease in AQP5 level. However, whether AQP5 deficiency can cause dry eye has not been proved yet.
Circular RNA (circRNA), as a specific and specific non-coding RNA, has a closed circular structure and is more stable than long non-coding RNAs (Jiang et al., 2018). A recent study has shown that circRNA molecules are rich in microRNA (miRNA) binding sites, which play the role of an miRNA sponge in cells, thus relieving the inhibition of miRNA on its target genes and increasing the expression levels of those target genes (Kulcheski et al., 2016). CircRNA might be a biomarker for primary Sjögren’s syndrome, correlating with the primary Sjögren’s syndrome etiology such as has-circRNA-001264, has-circRNA-104121, and has-circRNA-045355 (Su et al., 2019). Meanwhile, miR-146a and miR-155 have also been reported to be related to dry eye (Shi et al., 2014). Therefore, the expression and the clinical significance of circRNA in dry eye need to be studied.
In this study, we used aquaporin 5 knockout (AQP5–/–) mice that exhibit dry eye characteristics and performed hematoxylin–eosin staining to determine structural changes in the lacrimal glands of these mice. We found that the lacrimal glands of AQP5–/– mice exhibited abnormal changes compared with those of AQP5+/+ mice. The differential expression of circRNA was also found by high-throughput sequencing, and a circRNA–miRNA–mRNA network related to phagosomes was predicted. These data suggest that the lack of AQP5 may cause a differential expression of circRNA, which may, in turn, lead to the development of primary dry eye.
Materials and Methods
Animals
Using CRISPR/Cas9 technology, AQP5–/– mice were produced by the high-flux electric transfer of fertilized eggs from Cyagen Biosciences Inc. (Guangzhou, China). We used age-matched AQP5+/+ and AQP5–/– mice for our study. All experimental and animal care procedures were followed according to the ARVO Statement for the Use of Animals in Ophthalmic and Vision Research, and the study was approved by the Animal Care and Use Committee of Qingdao University (Qingdao, China). Spontaneous defects of the corneal epithelium were visualized by instilling 0.25% fluorescein sodium and photographing under a slit lamp (66 Vision Tech. Co., Ltd., Suzhou, China).
Tear production was examined using the phenol red thread (Jingming) test as described previously (Stevenson et al., 2014). In brief, the thread was placed on the palpebral conjunctiva of the lower eyelid at one-third of the distance from the lateral canthus for 20 s. The length of the wet portion was measured in millimeters.
Hematoxylin–Eosin Staining and Transmission Electron Microscopy
AQP5+/+ and AQP5–/– mice were sacrificed by cervical dislocation, and the lacrimal glands were removed for further use. The lacrimal glands were fixed in 10% buffered formalin and embedded in paraffin. Paraffin sections (4 μm) were produced using paraffin-embedded tissues. The sections were stained with hematoxylin–eosin and observed under a light microscope (Nikon Eclipse E100, Nikon, Japan). Samples of approximately 1 mm3 were dissected from the lacrimal glands and subsequently fixed with electron microscope fixation fluid (Servicebio, G1102), post-fixed with 1% osmic acid in 0.1 mol/L phosphate buffer, dehydrated in a graded series of ethanol, and embedded in embedding agent (SPI, 90529-77-4). The specimens were cut into 60–80 nm sections by using an ultrathin microtome (Leica UC7, Leica). The sections were double-stained with uranium and lead and observed under a transmission electron microscope (HITACHI, HT7700).
Extraction and Separation of RNA Samples
According to the manufacturer’s instructions (Lu et al., 2018), Trizol (Invitrogen, Carlsbad, CA, United States) was used to extract RNA from the lacrimal glands of both AQP5+/+ and AQP5–/– mice, and the optical density (OD260/280) value was determined by using a NanoDrop ND-2000 instrument (Thermo Fisher Scientific, Waltham, MA, United States) to measure the RNA concentration. RNA integrity was determined by agarose gel electrophoresis.
High-Throughput Sequencing
High-throughput transcriptome sequencing and bioinformatics analysis were performed using Cloud-Seq Biotech (Shanghai, China). Specifically, total RNA was treated with a Ribo-Zero rRNA Removal kit (Illumina, San Diego, CA, United States) to remove ribosomal RNA (rRNA). An RNA library was constructed according to the manufacturer’s instructions using the purified RNA samples and the TruSeq Stranded Total RNA Library Prep kit (Illumina, San Diego, CA, United States). The BioAnalyzer 2100 system was used for ensuring quality control and library quantification. Then, 10-pM libraries were denatured into single-stranded DNA molecules, captured on Illumina flow cells, amplified in situ, and clustered, and finally, 150 cycles of sequencing were performed on the Illumina HiSeq sequencer according to the manufacturer’s instructions (Lu et al., 2018).
CircRNA Sequencing Analysis
The paired terminal readings were obtained from the Illumina HiSeq 4000 sequencer, and quality control was performed through Q30. High-quality reads were primarily screened by Cutadapt software (version 1.9.3) (Martin, 2011), and after performing 3′ adapter trimming, low-quality reads were removed. We used STAR software (version 2.5.1b) (Dobin et al., 2013) to map and align high-quality reads with the reference genome/transcriptome. Next, we selected some nucleotide sequences from the reads as anchor points and input the results into the DCC software (version 0.4.4) (Cheng et al., 2016) that eventually compared connected and unconnected reads to identify possible circRNAs. We used EdgeR software (version 3.16.5) (Robinson et al., 2010) to normalize the data and analyze the differential expression of the identified circRNAs.
Analysis of Differential CircRNAs and mRNAs
The differential expression of circRNAs and mRNAs between the AQP5+/+ and the AQP5–/– mice groups was calculated using standardized readings. CircRNAs and mRNAs with a fold change of ≥2.0 and a p < 0.05 were considered to indicate a differential expression that was statistically significant.
Validation of Differentially Expressed CircRNAs, miRNAs, and mRNAs
The reliability of high-throughput RNA sequencing was verified by quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR). Among all the identified differentially expressed circRNAs and miRNAs, nine upregulated circRNAs and six upregulated mRNAs were selected for validation. Glyceraldehyde-3-phosphate dehydrogenase was used as a reference for standardization. Among all the identified differentially expressed miRNAs, four downregulated miRNAs were selected for validation. U6 was used as a reference for standardization. Total RNA was reverse-transcribed into complementary DNA using the PrimeScript RT kit (Perfect Real Time; Takara, Osaka, Japan), and qRT-PCR was performed using the Applied Biosystems 7500 Fast Real-Time PCR system. Three independent experiments were performed on all samples. The expression was determined using a threshold cycle, and the relative expression level was calculated using the 2–ΔΔCT method. The primers used for all the selected circRNAs are specified in Table 1. The primers for all the selected mRNAs are specified in Table 2. The primers used for miRNA reverse transcription are specified in Table 3. The primers for all the selected miRNAs used in qRT-PCR are specified in Table 4.
GO and KEGG Pathway Analyses
GO and KEGG were used to analyze genes related to differentially expressed circRNAs and mRNAs. GO analysis is divided into three aspects: molecular function, biological process, and cell composition. We ranked the top 10 enriched GO terms according to the p-value. KEGG pathway analysis was used to analyze the possible biological functions of differentially expressed circRNA and mRNA.
Analysis of the CircRNA–miRNA–mRNA Network and Related Prediction
The miRNA binding sites and target mRNAs were predicted using proprietary software based on TargetScan (version 7.0) and MiRanda (version 3.3a). Based on the prediction results, we used the Cycloscape software (version 3.1.0) to construct a network map of circRNA–miRNA–mRNA.
RNA Immunoprecipitation–qPCR
RNA immunoprecipitation (RIP) analysis was performed using the EZ-Magna RIP RNA binding protein immunoprecipitation kit (Millipore, Billerica, MA, United States) and the RIP protocol. The lacrimal glands of AQP5+/+ mice were cleaved and incubated with RIP buffer containing magnetic beads of anti-AGO2 antibody and immunoglobulin G (IgG, Abcam) as a negative control. To analyze the enrichment of circRNAs and mRNAs, qRT-PCR was performed to analyze the extracted RNA. The primers used for all the selected circRNAs and mRNAs are specified in Tables 1, 2.
Data Analysis
The results were expressed as mean ± standard deviation (SD). All the experimental data were analyzed by GraphPad Prism 7.0 (GraphPad Software Inc., La Jolla, CA, United States). A p < 0.05 was considered as statistically significant. The data are representative of at least three experiments.
Results
Dry-Eye-Like Characteristics Found in AQP5–/– Mice
In the course of our routine breeding of AQP5–/– mice, we found that these mice exhibited spontaneous dry eye symptoms. Under a slit lamp, a spontaneous punctate epithelial defect in the cornea of these mice could be observed (Figure 1A). Then, we used a phenol red cotton thread to measure tears in AQP5–/– mice by using AQP5+/+ mice as control. The tear volume of AQP5–/– mice was significantly lesser than that of AQP5+/+ mice (Figure 1B). However, there was no difference between the tear volume produced by male and female AQP5–/– mice (Figure 1C).
Figure 1. Effect of AQP5 knockout on the tear secretion of mice. (A) The results of sodium fluorescein staining in the cornea of AQP5+/+ mice and AQP5–/– mice at different ages. The staining point was a corneal defect. (B,C) The lacrimal secretion of AQP5+/+ mice and AQP5–/– mice was measured using the phenol red cotton thread method. The wet length of phenol red cotton thread in mice of different ages (B) and genders (C). The data are expressed as mean ± SD (n = 14 per group). *p < 0.05.
AQP5 Deficiency Changed the Structure of Epithelial Cells in Lacrimal Glands
Hematoxylin–eosin staining revealed that there were more vacuoles in the epithelial cells of the lacrimal glands of AQP5–/– mice than in those of AQP5+/+ mice, and the arrangement of the acini was disordered in AQP5–/– mice than in AQP5+/+ mice (Figure 2A). The specific observations in AQP5–/– mice were as follows: the number of acini decreased (Figure 2BI), the area of single acinus increased (Figure 2BII), and the number of epithelial cells decreased (Figure 2BIII), and at the same time, the number and the proportion of vacuoles in lacrimal gland epithelial cells per unit area increased significantly (Figures 2BIV,V).
Figure 2. Structural changes of lacrimal glands in AQP5–/– mice compared to AQP5+/+ mice. (A) Hematoxylin–eosin staining of lacrimal glands from AQP5+/+ mice and AQP5–/– mice. Arrows indicate vacuoles inside the cytoplasm of AQP5–/– acinar cells. (B) The number (I) and the area of single acini (II) in the unit area, the number of lacrimal gland epithelial cells in the unit area (III), and the number (IV) and area of vacuoles (V) in the unit area were compared between AQP5–/– and AQP5+/+mice. The data are expressed as mean ± SD (n = 6 per group). *p < 0.05. (C) Transmission electron microscopy of lacrimal glands from AQP5+/+ mice and AQP5–/– mice. N, nuclei; Nu, nucleolus; M, mitochondrion; RER, rough endoplasmic reticulum; SG, secretory granule; TJ, tight junction; DE, desmosome; AP, autophagosome; triangles, intercellular space.
To further study the structural changes in the epithelial cells of lacrimal glands, the lacrimal glands of both groups of mice were subjected to electron microscopy analysis. The results showed that compared with the epithelial cells of AQP5+/+ mice, those of AQP5–/– mice were seriously damaged, the mitochondria were swollen and enlarged, a large number of mitochondrial cristae had disappeared, the membrane had disintegrated, some secretory granules were dissolved, the secretory granules were fused, the levels of autophagy bodies increased, the cell gap was enlarged, and the endoplasmic reticulum around the cell membrane had expanded (Figure 2C).
Overview of CircRNA Expression in Lacrimal Glands
To study the effect of AQP5 deficiency on lacrimal gland performance, we obtained lacrimal glands from AQP5–/– and AQP5+/+ mice and sequenced their RNA with high throughput to determine the expression profile of circRNA in the glands. A total of 2,541 circRNAs were identified in the lacrimal glands of the two groups, 937 of which have not been reported before (Figure 3A).
Figure 3. Expression pattern of circRNA in the lacrimal glands of C57BL/6 mice. (A) The proportion of newly discovered circRNA in all identified circRNA. (B) Distribution of circRNA on chromosomes. (C) The genomic origin of detected circRNA. (D) The length distribution of exonic circRNA.
Most circRNAs were located on chromosomes 1–19 (Figure 3B). There was a variety of catalogs of circRNAs, and most of them were exonic circRNAs (Figure 3C). The size of these circRNAs ranged from 82 nucleotides (nt) to more than 10,000 nt (Figure 3D). The overall average length was 2,824 nt.
Differential Expression of CircRNAs in the Lacrimal Glands of AQP5–/– and AQP5+/+ Mice
Among the 2,541 identified circRNAs, 920 could only be detected in AQP5+/+ mice, 948 could only be detected in AQP5–/– mice, and 673 could be detected in both types of mice (Figure 4A). Thirty circRNAs were differentially expressed in AQP5–/– mice compared to AQP5+/+ mice, wherein nine of them were upregulated, while 21 were downregulated. Hierarchical clustering showed a distinguishable circRNA expression profile between AQP5–/– and AQP5+/+ mice (Figure 4B).
Figure 4. Identification of differential expression circRNA after AQP5 knockout. (A) Venn diagram of differentially expressed circRNA; 673 circRNAs were detected in both mice. (B) Clustering graph. Hierarchical clustering displayed the circRNA expression profile of AQP5+/+ mice (n = 3) vs. AQP5+/+ mice (n = 3). (C) The expression of nine upregulated circRNAs in AQP5+/+ and AQP5+/+ mice was detected by real-time fluorescence quantitative PCR. *p < 0.05.
To verify the high-throughput sequencing results, quantitative real-time polymerase chain reaction was performed for the nine upregulated circRNAs. The relative expression level of the nine circRNAs was consistent with that determined by high-throughput sequencing (Figure 4C).
Function of Differentially Expressed CircRNAs in the Lacrimal Glands of AQP5–/– and AQP5+/+ Mice
To preferably understand the biological function of the differentially expressed circRNAs present in the lacrimal glands of AQP5–/– mice, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. The GO terms mainly covered three areas: biological processes, cellular components, and molecular functions. We found that the most abundant GO terms in the upregulated circRNAs were RNA phosphodiester in the endoplasmic reticulum, cell surface furrow, and ribonuclease activity. According to the KEGG analysis, the pathway related to the upregulation of circRNAs was protein processing in the endoplasmic reticulum (Figure 5A).
Figure 5. GO analysis and KEGG analysis of circRNA. (A) GO analysis and KEGG analysis of upregulated circRNA. GO analysis identified (I) molecular functions, (II) biological processes, and (III) cellular components. (IV) Relevant pathways were identified for upregulated circRNA. (B) GO analysis and KEGG pathway analysis of downregulated circRNA. GO analysis identified (I) molecular functions, (II) biological processes, and (III) cellular components. (IV) Relevant pathways were identified for downregulated circRNA.
In the downregulated circRNAs, the most abundant GO terms were positive regulation of cellular process, nuclear part, and high-voltage-gated calcium channel activity. The most relevant pathway of these circRNAs was the synaptic vesicle cycle (Figure 5B).
Differential Expression of mRNAs in the Lacrimal Glands of AQP5–/– and AQP5+/+ Mice
By high-throughput analysis, 15,243 messenger RNAs (mRNAs) were detected in the lacrimal glands of both AQP5–/– and AQP5+/+ mice. Among these mRNAs, 549 could only be detected in AQP5+/+ mice, 325 could only be detected in AQP5–/– mice, and 14,369 could be detected in both types of mice (Figure 6A). Compared to the mRNAs found in AQP5+/+ mice, 515 mRNAs were differentially expressed in AQP5–/– mice, wherein 229 were upregulated and 286 were downregulated. Hierarchical clustering showed a distinguishable mRNA expression profile between AQP5–/– and AQP5+/+ mice (Figure 6B).
Figure 6. Identification of differential expression of mRNA after AQP5 knockout. (A) Venn diagram of differentially expressed mRNA; 15,243 mRNAs were detected in both mice. (B) Clustering graph. Hierarchical clustering displayed the mRNA expression profile of the AQP5–/– mice (n = 3) vs. the AQP5+/+ mice (n = 3). (C) The expression of six upregulated mRNAs in AQP5+/+ and AQP5+/+ mice was detected by real-time fluorescence quantitative PCR. (D) The expression of four downregulated miRNAs in AQP5+/+ and AQP5+/+ mice was detected by real-time fluorescence quantitative PCR. *p < 0.05.
To verify the high-throughput sequencing results, qRT-PCR was performed for the six upregulated mRNAs and four downregulated miRNAs. The relative expression levels of the six mRNAs (Figure 6C) and the four miRNAs (Figure 6D) were consistent with the sequencing results.
Function of Differentially Expressed mRNAs in the Lacrimal Glands of AQP5–/– and AQP5+/+ Mice
We also performed GO and KEGG pathway analyses on the differentially expressed mRNAs to determine their biological function in the lacrimal glands of AQP5–/– mice. We found that the most abundant GO terms for the upregulated mRNAs were response to external stimulus, extracellular space, and glycosaminoglycan binding. According to the KEGG analysis, the pathway related to circRNA upregulation was phagosome (Figure 7A).
Figure 7. GO analysis and KEGG pathway analysis of mRNA. (A) GO analysis and KEGG pathway analysis of upregulated mRNA. GO analysis identified (I) molecular functions, (II) biological processes, and (III) cellular components. (IV) Relevant pathways were identified for upregulated mRNA. (B) GO analysis and KEGG pathway analysis of downregulated mRNA. GO identified (I) molecular functions, (II) biological processes, and (III) cellular components. (IV) Relevant pathways were identified for downregulated mRNA.
For the downregulated mRNAs, the most abundant GO terms were single-organism process, intracellular part, and metal ion binding. The most relevant pathway pertaining to the downregulated mRNAs was vitamin B6 metabolism (Figure 7B).
CircRNA–miRNA–mRNA Network Analysis
To further determine the mechanism(s) underlying the functions of the identified circRNAs, we constructed a ceRNA regulatory network of circRNA–miRNA–mRNA. We targeted miRNAs and mRNAs using miRNA target prediction software according to TargetScan and miRanda. The top five miRNAs associated with each circRNA were involved in the network. All the targeted mRNAs of each of the top five miRNAs matched with the upregulated mRNAs that were identified by mRNA sequencing.
The network was constructed by using nine upregulated validated circRNAs, 40 altered miRNAs, and nine upregulated mRNAs (Figure 8). In this complex network, one circRNA can associate with multiple miRNAs and one miRNA can regulate multiple mRNAs. According to the comparison with the database, these nine mRNAs and phagosome pathways are also closely related.
Figure 8. circRNA–miRNA–mRNA network. (A) Circles represent circRNA, arrowhead represents miRNA, and squares represent mRNA. (B) RIP-qPCR detected (I) chr5:147450653-147455188+, (II) chr9:3441055-3460131+, (III) 5430435G22Rik, and (IV) Thbs1. IG, negative control immunoglobulin G group; IP, containing magnetic beads of anti-AGO2 antibody group. *p < 0.05.
To further study the function of the previously identified circRNAs and verify the reliability of the predicted network of circRNA–miRNA–mRNA, RIP-qPCR was performed. As shown in Figure 8, compared with the IgG group (control), chr5:147450653-147455188 + (Figure 8BI) and chr9:3441055-3460131 + (Figure 8BII) can combine with the AGO2–miRNA complex, suggesting that chr5:147450653-147455188 + and chr9:3441055-3460131 + can combine with miRNA through a sponge mechanism. In addition, compared with the IgG group (control), 5430435G22Rik (Figure 8BIII) and Thbs1 (Figure 8BIV), identified in the previously predicted network graph, could also be detected, which indicates that our prediction is reliable.
Discussion
Dry eye is a common chronic ocular surface disease caused by many factors such as reduced tear film stability, increase in tear osmolality, and ocular surface inflammation (Nelson et al., 2017). With changes in people’s working environments and living habits, the incidence of dry eye is increasing annually, but there are few effective drugs to treat it. This lack of effective drugs is mainly related to the diverse physiological symptoms associated with dry eye in humans and the lack of corresponding animal models for studying these diverse symptoms in detail.
Previous studies have used animal models of dry eye by administering drugs to inhibit the secretion of tears or by removing the lacrimal or Harderian glands of mice (Gilbard et al., 1987; Viau et al., 2008; Fakih et al., 2019). However, these methods are too complex and cumbersome, the efficiency of the models is not stable, and the damage to the mice is greater. Therefore, selecting a simpler and more stable animal model of dry eye for studying its pathogenesis and for treatment drug screening is particularly important. In this study, we found that AQP5–/– mice naturally develop dry eye symptoms from birth, which makes them a stable, simple, and effective dry eye research animal model that can be used to study the pathogenesis of this condition. The AQP5 level has been considered to decrease with the occurrence of dry eye and is usually regarded as an indicator of dry eye in fundamental research (Wang et al., 2017). However, whether dry eye occurs because of AQP5 knockout has not been reported.
The occurrence of dry eye is usually related to age and gender; the incidence gradually increases with age, and women are more likely to develop dry eye than men (Kaštelan et al., 2013). However, in this study, we did not observe this phenomenon. AQP5–/– mice had dry eyes since birth, and no age- and gender-related trends were observed. This may be because the type of dry eye after AQP5 knockout is primary.
Hematoxylin–eosin staining of the lacrimal gland specimens of AQP5–/– mice revealed that the structure of the acini was disordered, and there were many changes in the vacuolation structure in the acini. The same phenomenon was observed in TSP–/– mice, a mouse model of Sjögren’s syndrome (an autoimmune disease) (Bhattacharya et al., 2018). Electron microscopic observations revealed that a large number of mitochondrial cristae had disappeared, and there was an increase in the number of autophagosomes. These results suggest that the structural changes in the lacrimal gland epithelial cells of AQP5–/– mice may be related to autophagy.
As new non-coding RNAs, circRNAs play a key role in many ophthalmic diseases such as macular degeneration and corneal neovascularization (Liu et al., 2020; Wu et al., 2020). However, the role of circRNAs in dry eye and lacrimal gland development has not been studied before. Thus, to better understand the complex pathogenesis of dry eye, identification of comprehensive circRNA expression profiles in the lacrimal glands of AQP5–/– mice was critical. We identified 2,541 circRNAs, of which 937 were newly identified. We also identified 30 differentially expressed circRNAs in the lacrimal glands of AQP5–/– mice compared with those of control mice. We selected nine upregulated circRNAs for validation. The qRT-PCR results were consistent with those of high-throughput sequencing, which proved the reliability of high-throughput sequencing.
In the circRNA–miRN–mRNA regulatory network or the ceRNA theory, circRNAs play the role of miRNA sponge and influence the expression of target mRNAs by competitive adsorption of miRNAs (Qi et al., 2015). We found the circRNA and the mRNA with the same trend by sequencing circRNA and mRNA and then find the circRNA and mRNA with the same miRNA binding site by bioinformatics analysis. Only by meeting the above two conditions can the circRNA and mRNA be used to construct the circRNA–miRNA–mRNA regulatory network. Based on the nine validated upregulated circRNAs and the top five miRNAs predicted by using circRNAs, we predicted a circRNA–miRNA–mRNA network. In addition, by performing RIP-qPCR experiments, we proved that our predicted network is reliable. As a known member of the Argonaute protein family, Argonaute 2 (AGO2) can combine with miRNA to form the AGO2–miRNA complex (Zhang et al., 2020). This complex can be combined with circRNA via the sponge mechanism, which proves that chr5:147450653-147455188 + and chr9:3441055-3460131 + can play a role through sponge mechanism. In addition, the circRNAs and the mRNAs identified in our previously predicted network can also be detected by combining with the AGO2–miRNA complex, which shows that our network is reliable.
Interestingly, the target genes involved in this network are all related to phagosomes, according to the comparison with the database. These observations highly correlated with the results of our KEGG analysis of upregulated mRNA. Phagosomes, as a highly dynamic organelle, play an important role in innate and adaptive immunity as well as in tissue homeostasis (Dean et al., 2019). Our results indicate that the circRNA–miRNA–mRNA network predicted in our study is a network of phagosomes. What is more surprising is that our prediction is consistent with the structural changes in the lacrimal gland epithelial cells that we mentioned before. Therefore, these circRNAs may be involved in the regulatory networks and the pathways of phagosomes involved in the pathogenesis of dry eye. However, the regulatory mechanism of circRNAs is very complex; hence, it is necessary to further study the function and the mechanism of these circRNAs.
Conclusion
The present study demonstrates that (i) AQP5 knockout can lead to primary dry eye development in mice (AQP5–/– mice have a stable dry eye phenotype since birth) and this phenotype may be produced by altering the structure of lacrimal glands, (ii) circRNA levels are significantly altered in the lacrimal glands of AQP5–/– mice, and (iii) the interactions of the circRNA–miRNA–mRNA network associated with phagosomes may regulate the expression of AQP5 involved in the pathogenesis of dry eye.
Data Availability Statement
The data of sequencing results for this study can be found in the GEO database (GSE149832). Please see https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE149832 for more details.
Ethics Statement
All experimental and animal care procedures followed the ARVO Statement for the Use of Animals in Ophthalmic and Vision Research and were approved by the Animal Care and Use Committee of Qingdao University (Qingdao, China).
Author Contributions
GD and PC conceived and designed the study. YL, SH, XX, and XW performed the experiments. GD and TZ analyzed the data. YL and PC wrote the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (Grant No. 81970782), Shandong Provincial Natural Science Foundation (Grant No. ZR2018MH016), Qingdao Postdoctoral Application Research Project (Grant No. 40518060071), and China Postdoctoral Science Foundation (Grant No. 2017M612211). All the above funds are provided by the corresponding author (PC).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We cordially thank Ms. Meifang Dai for assistance in data analysis and Editage (www.editage.com) for English language editing.
Abbreviations
AGO2, argonauute 2; AQP5, aquaporin 5; circRNA, circular RNA; GAPDH, glyceraldehyde 3-phosphate dehydrogenase; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; miRNA, microRNA; mRNA, messenger RNA; nt, nucleotides; qRT-PCR, quantitative real-time polymerase chain reaction; rRNAs, ribosomal RNA; RIP, RNA immunoprecipitation; SD, standard deviation.
References
Bhattacharya, S., García-Posadas, L., Hodges, R. R., Makarenkova, H. P., Masli, S., and Dartt, D. A. (2018). Alteration in nerves and neurotransmitter stimulation of lacrimal gland secretion in the TSP-1-/- mouse model of aqueous deficiency dry eye. Mucosal Immunol. 11, 1138–1148. doi: 10.1038/s41385-018-0002-y
Cheng, J., Metge, F., and Dieterich, C. (2016). Specific identification and quantification of circular RNAs from sequencing data. Bioinformatics 32, 1094–1096. doi: 10.1093/bioinformatics/btv656
Craig, J. P., Nichols, K. K., Akpek, E. K., Caffery, B., Dua, H. S., Joo, C.-K., et al. (2017). TFOS DEWS II definition and classification report. Ocul. Surf. 15, 276–283.
Dean, P., Heunis, T., Härtlova, A., and Trost, M. (2019). Regulation of phagosome functions by post-translational modifications: a new paradigm. Curr. Opin. Chem. Biol. 48, 73–80. doi: 10.1016/j.cbpa.2018.11.001
Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., et al. (2013). STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21. doi: 10.1093/bioinformatics/bts635
Fakih, D., Zhao, Z., Nicolle, P., Reboussin, E., Joubert, F., Luzu, J., et al. (2019). Chronic dry eye induced corneal hypersensitivity, neuroinflammatory responses, and synaptic plasticity in the mouse trigeminal brainstem. J. Neuroinflamm. 16, 1–20.
Gilbard, J., Rossi, S. R., and Gray, K. L. (1987). A new rabbit model for keratoconjunctivitis sicca. Invest. ophthalmol. Vis. Sci. 28, 225–228.
Jiang, L.-H., Sun, D.-W., Hou, J.-C., and Ji, Z.-L. (2018). CircRNA: a novel type of biomarker for cancer. Breast Cancer 25, 1–7. doi: 10.1007/s12282-017-0793-9
Kaštelan, S., Tomić, M., Salopek-Rabatić, J., and Novak, B. (2013). Diagnostic procedures and management of dry eye. Biomed Res. Int. 2013:309723.
Kulcheski, F. R., Christoff, A. P., and Margis, R. (2016). Circular RNAs are miRNA sponges and can be used as a new class of biomarker. J. Biotechnol. 238, 42–51. doi: 10.1016/j.jbiotec.2016.09.011
Kumari, S. S., Varadaraj, M., Menon, A. G., and Varadaraj, K. (2018). Aquaporin 5 promotes corneal wound healing. Exp. Eye Res. 172, 152–158. doi: 10.1016/j.exer.2018.04.005
Lemp, M. A., and Foulks, G. N. (2007). The definition and classification of dry eye disease. Ocul. Surf. 5, 75–92. doi: 10.1016/s1542-0124(12)70081-2
Levin, M. H., and Verkman, A. (2006). Aquaporin-3-dependent cell migration and proliferation during corneal re-epithelialization. Invest. Ophthalmol. Vis. Sci. 47, 4365–4372.
Lin, H., Liu, Y., and Yiu, S. (2019). Three dimensional culture of potential epithelial progenitor cells in human lacrimal gland. Transl. Vis. Sci. Technol. 8, 32–32.
Liu, X., Zhang, L., Wang, J. H., Zeng, H., Zou, J., Tan, W., et al. (2020). Investigation of circRNA expression profiles and analysis of circRNA-miRNA-mRNA networks in an animal (Mouse) model of age-related macular degeneration. Curr. Eye Res. 6, 1–8. doi: 10.1080/02713683.2020.1722179
Lu, Q., Gong, W., Wang, J., Ji, K., Wang, Y., Xu, C., et al. (2018). Identification of circular RNAs altered in mouse jejuna after radiation. Cell. Physiol. Biochem. 47, 2558–2568. doi: 10.1159/000491652
Ma, T., Song, Y., Gillespie, A., Carlson, E. J., Epstein, C. J., and Verkman, A. S. (1999). Defective secretion of saliva in transgenic mice lacking aquaporin-5 water channels. J. Biol. Chem. 274, 20071–20074. doi: 10.1074/jbc.274.29.20071
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 17, 10–12.
Messmer, E. M. (2015). The pathophysiology, diagnosis, and treatment of dry eye disease. Dtsch. Arztebl. Int. 112, 71–81. doi: 10.3238/arztebl.2015.0071
Nelson, J. D., Craig, J. P., Akpek, E. K., Azar, D. T., Belmonte, C., Bron, A. J., et al. (2017). TFOS DEWS II introduction. Ocul. Surf. 15, 269–275.
Qi, X., Zhang, D. H., Wu, N., Xiao, J. H., Wang, X., and Ma, W. (2015). ceRNA in cancer: possible functions and clinical implications. J. Med. Genet. 52, 710–718. doi: 10.1136/jmedgenet-2015-103334
Raina, S., Preston, G. M., Guggino, W. B., and Agre, P. (1995). Molecular cloning and characterization of an aquaporin cDNA from salivary, lacrimal, and respiratory tissues. J. Biol. Chem. 270, 1908–1912. doi: 10.1074/jbc.270.4.1908
Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616
Shi, H., Zheng, L. Y., Zhang, P., and Yu, C. Q. (2014). miR−146a and miR−155 expression in PBMC s from patients with Sjögren’s syndrome. J. Oral pathol. Med. 43, 792–797. doi: 10.1111/jop.12187
Sindhu Kumari, S., and Varadaraj, K. (2013). Aquaporin 5 knockout mouse lens develops hyperglycemic cataract. Biochem. Biophys. Res. Commun. 441, 333–338. doi: 10.1016/j.bbrc.2013.10.058
Stevenson, W., Chen, Y., Lee, S.-M., Lee, H. S., Hua, J., Dohlman, T., et al. (2014). Extraorbital lacrimal gland excision: a reproducible model of severe aqueous tear-deficient dry eye disease. Cornea 33, 1336–1341. doi: 10.1097/ico.0000000000000264
Su, L.-C., Xu, W.-D., Liu, X.-Y., Fu, L., and Huang, A.-F. (2019). Altered expression of circular RNA in primary Sjögren’s syndrome. Clin. Rheumatol. 38, 3425–3433. doi: 10.1007/s10067-019-04728-6
Takata, K., Matsuzaki, T., and Tajika, Y. (2004). Aquaporins: water channel proteins of the cell membrane. Prog. Histochem. cytochem. 39, 1–83. doi: 10.1016/j.proghi.2004.03.001
Thiagarajah, J. R., and Verkman, A. S. (2002). Aquaporin deletion in mice reduces corneal water permeability and delays restoration of transparency after swelling. J. Biol. Chem. 277, 19139–19144. doi: 10.1074/jbc.M202071200
Viau, S., Maire, M.-A., Pasquis, B., Grégoire, S., Fourgeux, C., Acar, N., et al. (2008). Time course of ocular surface and lacrimal gland changes in a new scopolamine-induced dry eye model. Graefes Arch. Clin. Exp. Ophthalmol. 246, 857–867. doi: 10.1007/s00417-008-0784-9
Wang, S.-Q., Wang, Y.-X., and Hua, H. (2017). Characteristics of labial gland mesenchymal stem cells of healthy individuals and patients with sjögren’s syndrome: a preliminary study. Stem Cell. Dev. 26, 1171–1185. doi: 10.1089/scd.2017.0045
Wu, P., Zhang, D., Geng, Y., Li, R., and Zhang, Y. (2020). Circular RNA-ZNF609 regulates corneal neovascularization by acting as a sponge of miR-184. Exp. Eye Res. 192:107937. doi: 10.1016/j.exer.2020.107937
Keywords: AQP5, lacrimal gland, dry eye, circRNA, miRNA, mRNA
Citation: Liu Y, Di G, Hu S, Zhao T, Xu X, Wang X and Chen P (2020) Expression Profiles of CircRNA and mRNA in Lacrimal Glands of AQP5–/– Mice With Primary Dry Eye. Front. Physiol. 11:1010. doi: 10.3389/fphys.2020.01010
Received: 11 March 2020; Accepted: 24 July 2020;
Published: 03 September 2020.
Edited by:
Maisa Hanna-Maija Seppala, King’s College London, United KingdomCopyright © 2020 Liu, Di, Hu, Zhao, Xu, Wang and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Peng Chen, Y2hlbnBlbmc1OTkyMDVAMTI2LmNvbQ==
†These authors have contributed equally to this work