- 1Center for Reproductive Medicine, Renji Hospital, School of Medicine, Shanghai Jiao Tong University, Shanghai, China
- 2Shanghai Key Laboratory for Assisted Reproduction and Reproductive Genetics, Shanghai, China
- 3Animal Laboratory, Renji Hospital, School of Medicine, Shanghai Jiao Tong University, Shanghai, China
- 4International Peace Maternity and Child Health Hospital, School of Medicine, Shanghai Jiao Tong University, Shanghai, China
- 5Shanghai Key Laboratory of Embryo Original Diseases, Shanghai, China
- 6Hospital of Obstetrics and Gynecology, Fudan University, Shanghai, China
High maternal estradiol is reported to induce metabolic disorders by modulating hypothalamic gene expression in offspring. Since neurogenesis plays a crucial role during hypothalamus development, we explored whether prenatal high estradiol exposure (HE) affects proliferation and differentiation of fetal hypothalamic neural stem/progenitor cells (NSC/NPCs) in mice and performed RNA sequencing to identify the critical genes involved. NSC/NPCs in HE mice presented attenuated cell proliferation but increased neuronal differentiation in vitro compared with control (NC) cells. Gene set enrichment analysis of mRNA profiles indicated that genes downregulated in HE NSC/NPCs were enriched in neurogenesis-related Gene Ontology (GO) terms, while genes upregulated in HE NSC/NPCs were enriched in response to estradiol. Protein-protein interaction analysis of genes with core enrichment in GO terms of neurogenesis and response to estradiol identified 10 Hub mRNAs, among which three were potentially correlated with six differentially expressed (DE) lncRNAs based on lncRNA profiling and co-expression analysis. These findings offer important insights into developmental modifications in hypothalamic NSC/NPCs and may provide new clues for further investigation on maternal environment programmed neural development disorders.
Introduction
The theory that the intrauterine environment can influence prenatal development and the future health of offspring (Bateson et al., 2004) has resulted in increased interest in the developmental origin of chronic disease. Ovulation induction clinically used in assisted reproductive technology generates a supraphysiologic level of blood estradiol, which may predispose offspring to an abnormal intrauterine environment after fresh embryo transfer (Hu et al., 2014). Our previous study reported a programming effect by high maternal estradiol during early pregnancy on the hypothalamic glucoregulatory system of male mouse offspring, which induces adult metabolic disorders (Wang et al., 2018). This discovery indicates that high prenatal estradiol levels alter hypothalamus development, though the underlying mechanism is not yet well understood.
Neural stem/progenitor cells (NSC/NPCs) generate neurons during a process called neurogenesis (Bond et al., 2015). In the central nervous system (CNS), neurogenesis begins in the embryonic stage and continues throughout life. NSC/NPCs are preprogrammed to form specific types of functional neurons even before birth (Fuentealba et al., 2015); therefore, the study of prenatal neurogenesis may help researchers understand the mechanisms underlying adult neural disorders.
The hypothalamus regulates the metabolic homeostasis of the whole body and is sensitive to adverse prenatal environment (Ralevski and Horvath, 2015), so the development of hypothalamic neurons could affect metabolism in later life. Hypothalamic NSC/NPCs are first generated between embryonic day 10.5 (E10.5) and E14.5 in rodents (Padilla et al., 2010), and can proliferate to form neurospheres and differentiate into neurons in vitro (Desai et al., 2011). Substantial evidence has shown that dysfunctions of hypothalamic NSC/NPCs are associated with metabolic disorders, such as obesity and glucose intolerance (Li et al., 2012, 2014; Livesey, 2012), but the corresponding whole genomic features of these NSC/NPCs are rarely reported.
RNA sequencing (RNA-Seq) is an effective approach to revealing genome alterations which provides expression information for all transcripts, including mRNAs and non-coding RNAs. LncRNAs may play a considerable regulatory role by interacting with mRNAs, so exploring the link between them may provide more clues for elucidating molecular pathways. Gene set enrichment analysis (GSEA) is a robust and tractable analytical method for interpreting RNA-Seq data, as it can detect differential enrichment of biological functions across an entire network of genes (Subramanian et al., 2005), without the limitations associated with the single-gene method.
Because a prenatally programmed adult hypothalamic disorder resulting from high maternal estradiol has been identified (Wang et al., 2018), here we explore whether neurogenesis potential is affected in fetal hypothalamic NSC/NPCs and attempt to identify the key lncRNA-mRNA network through integrated bioinformatic analysis. These findings may help us to understand molecular modifications of fetal hypothalamic neurogenesis resulting from an adverse intrauterine environment.
Materials and Methods
Animal Model and 5-Bromodeoxyuridine Labeling in Fetal Brain
A mouse model of prenatal high estradiol exposure was created based on our previously published method (Wang et al., 2018). Briefly, 8-week-old pregnant C57BL/6 mice received 100 μg/kg/d estradiol valerate [Sigma; high estradiol (HE) group] or an equal amount of blank solvent [corn oil; control (NC) group] via gavage from E5.5 to E11.5. For bromodeoxyuridine (BrdU) labeling, pregnant mice at E14.5 received a single intraperitoneal injection of 100 mg/kg BrdU (Sigma) 2 h before euthanized; then, the fetuses were extracted and decapitated. Male fetuses were identified by visual identification of testes next to the bladder using a dissection microscope. The brains of male fetuses were removed and fixed in 4% paraformaldehyde (PFA) for 24 h and then infiltrated with 20–30% sucrose. Brain sections of 20 μm were made using a freezing microtome (Leica) for immunofluorescence staining.
Tissue Immunofluorescence
Brain sections were blocked with 5% bovine serum albumin/0.3% Triton X-100 for 1 h at room temperature and incubated with primary antibodies mouse anti-Nestin (1:200, Millipore, catalog no. MAB353) and rabbit anti-BrdU (1:100, Abcam, catalog no. ab152095) overnight at 4°C, followed by reaction with secondary antibodies anti-rabbit Alexa Fluor 488 (1:200, Invitrogen, catalog no. A-11008) and anti-mouse Alexa Fluor 594 (1:200, Invitrogen, catalog no. A-11005) for 2 h at room temperature before counterstaining with 4',6-diamidino-2-phenylindole. The BrdU+Nestin+ cells were counted in five serial sections across the hypothalamus in each mouse.
NSC/NPC Isolation and Neurosphere Assay
The brains of E14.5 male fetuses were dissected quickly on ice to remove the hypothalami, which were then fragmented in Neurobasal-A (Gibco), digested with TrypLE (Gibco) in 37°C for 15 min, and gently triturated into single cells with tips. The cells were then washed twice in Neurobasal-A and suspended in a proliferation medium containing Neurobasal-A, 2% B27 (Gibco), 10 ng/ml EGF (PeproTech), 10 ng/ml bFGF (PeproTech), and 1% GlutaMAX (Gibco), seeded in ultralow adhesion 6-well plates at a density of 105/ml (Li et al., 2012), and incubated in 5% CO2 at 37°C. The neurospheres were photographed under a microscope for 4 days (Marshall et al., 2007), and the number was counted and diameter measured using the software ImageJ on the fourth day after isolation.
NSC/NPC Proliferation and Differentiation Assay
To assess the proliferation ability of NSC/NPCs, the primary neurospheres were digested to count single cells and passaged at a density of 105/ml in ultralow adhesion 6-well plates. The total cell number in each of the first four passages was calculated based on the assumption that all of the cells from the previous passage were replated. For the BrdU incorporation assay, primary NCS/NPCs were plated on Matrigel (BD)-coated coverslips at a density of 105/ml in proliferation medium and cultured for 24 h. Then, the cells were treated with 10 μm BrdU for 2 h before immunofluorescence staining.
For induced differentiation, primary NSC/NPCs were seeded as single cells at a density of 3 × 105/ml in Matrigel-coated coverslips placed in 24-well plates in differentiation medium containing Neurobasal-A, 2% B27, 1% fetal bovine serum (Gibco), and 1 μm retinoid acid (Sigma). The medium was changed every second day for 10 days, and then, the coverslips were removed to receive immunofluorescence detection of target neuron markers.
NSC/NPC Immunofluorescence
For immunofluorescence staining of NSC/NPCs, neurospheres were moved using tips to seed on Matrigel-coated coverslips for 20 min before detection, and cells on coverslips were fixed with 4% PFA for 15 min and blocked with 5% bovine serum albumin/0.3% Triton X-100 for 1 h at room temperature. Cells were then incubated with primary antibodies rabbit anti-Sox2 (1:400, Cell Signaling Technology, catalog no. 23064), mouse anti-Nestin (1:200, Millipore, catalog no. MAB353), rabbit anti-BrdU (1:100, Abcam, catalog no. ab152095), or mouse anti-Tuj1 (1:200, Cell Signaling Technology, catalog no. 4466) overnight at 4°C and with secondary antibodies anti-rabbit Alexa Fluor 488 (1:200, Invitrogen, catalog no. A-11008), anti-mouse Alexa Fluor 594 (1:200, Invitrogen, catalog no. A-11005), or anti-mouse Alexa Fluor 488 (1:200, Invitrogen, catalog no. A-11001) for 2 h at room temperature before counterstaining with 4',6-diamidino-2-phenylindole. The BrdU/Nestin and Tuj1 positive cells were counted in each group.
RNA-Seq Analysis
Three samples of first-passage NC and HE NSC/NPCs were harvested separately, each containing cells obtained from two mice. RNA was extracted with TRIzol (Invitrogen), its quality valued by spectrophotometer, and its integrity checked by Agilent 2,100 bioanalyzer. Total RNA was enriched by oligo beads, fragmented into small pieces, and reverse transcribed into cDNA. Second-strand cDNA was synthesized by DNA polymerase I with dUTP to construct a strand-specific library. The cDNA was then purified, end-repaired, poly A-added, and ligated to Illumina adaptor. The libraries were size-selected by agarose gel electrophoresis, PCR-amplified, and sequenced by Illumina NextSeq 500 by Personal Bio Co. (Shanghai, China).
The raw RNA-Seq data were filtered by removing low-quality and adaptor-related reads. The clean reads were then aligned to the mouse reference genome (10 mm) using Tophat2. Coding and non-coding transcripts were distinguished by Coding Potential Calculator, Coding-Non-Coding Index, and Pfam-scan. Non-coding RNAs with length >200 nt and exon number ≥2 were considered to be lncRNAs. Expression values were expressed as reads per kilobase per million reads. Differential expression analysis was conducted using DESeq2 (Love et al., 2014). LncRNAs and mRNAs with a log2 (fold change) ≥1 or ≤−1 and FDR < 0.05 were considered differentially expressed.
Bioinformatics Analysis
We conducted enrichment analyses using GSEA with the standard procedure obtained from the GSEA Web site.1 The number of permutations was set to 1,000, and FDR < 0.25 with p < 0.05 was considered statistically significant. We downloaded gene sets needed for Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analysis from the GSEA Web site. The enrichment bubble diagrams were made with R software (version 4.0.3).
A protein-protein interaction (PPI) network of selected mRNAs was constructed using STRING 11.0,2 and the Hub mRNAs [the top 10 nodes ranked by Maximal Clique Centrality (MCC); Chin et al., 2014] were identified using the plugin cytoHubba in Cytoscape software.3
To explore the lncRNA-mRNA regulatory network, Pearson’s correlation coefficient (PCC) between DE lncRNAs and Hub mRNAs was calculated and plotted using R software, and gene pairs with PCC ≥0.990 or ≤−0.990 and p < 0.05 were considered to be potentially correlated. The interaction network was visualized using Cytoscape software.
Quantitative Real-Time Polymerase Chain Reaction
Total RNA of NSC/NPCs was extracted using TRIzol (Invitrogen) and reverse-transcribed into cDNA using Primer Script RT Reagent Kit (Takara) and amplified with QuantiNova SYBR Green PCR Kit (QIAGEN) according to the manufacturer’s instructions. The thermocycling conditions were 95°C for 2 min, followed by 40 cycles of 95°C for 5 s and 60°C for 10 s. The primers are listed in Table 1. GAPDH was used as an endogenous control, and the relative expression level was analyzed using the 2−ΔΔCT method.
Statistical Analysis
Data aside from RNA-Seq and bioinformatics analysis were analyzed using the Statistical Package for Sciences Software, version 21.0 (IBM), and are presented as the mean ± standard error of the mean. Unpaired Student’s t-tests were used for comparisons between two groups, and p < 0.05 was considered statistically significant.
Results
Prenatal High Estradiol Affects Neurogenesis Potential of Fetal Hypothalamic NSC/NPCs
Serum estradiol after gavage in the HE pregnant mice reaches a peak value of four times that of the control group (Wang et al., 2018), forming a high maternal estradiol environment. Fetal brain sections were made on E14.5, and hypothalamic NSC/NPCs were isolated at the same time (Figure 1A). Authentic biomarkers Sox2 and Nestin were used to label NSC/NPCs (Suh et al., 2007; Gilyarov, 2008). Tissue immunofluorescence presented decreased number of BrdU+Nestin+ cells in HE fetal hypothalami after BrdU injection (Figures 1B,C), indicating a reduction of proliferating NSC/NPCs in vivo. Immunofluorescence staining of Sox2 and Nestin in neurospheres was performed for NSC/NPC identification after cell isolation (Figure 1D). The neurosphere assay showed a decreased amount of neurospheres with shorter average diameters in HE NSC/NPCs compared with NC on the fourth day (D4) of the first passage (Figures 1E,F), and the proliferation curve presented the accumulated NC NSC/NPC number significantly exceeded that in the HE group from the second to fourth passage (P2 to P4; Figure 1G), indicating attenuated proliferation ability in HE NSC/NPCs. We also performed BrdU incorporation assay in primary NSC/NPCs and found the decreased proportion of proliferating NSC/NPCs in vitro in HE group (Figures 1H,I).
Figure 1. Isolation and evaluation of fetal hypothalamic neural stem/progenitor cell (NSC/NPC) neurogenesis potential in a mouse model. (A) Schematic of the method used to generate fetal brain sections and hypothalamic NSC/NPCs in a mouse model. (B) Representative images of bromodeoxyuridine (BrdU) and Nestin immunofluorescence in fetal hypothalamic tissue; scale bar: 200 μm. Arrows show BrdU+Nestin+ cells. (C) Quantification of BrdU+Nestin+ cells in five serial fetal hypothalamic tissue sections (n = 3 mice per group). (D) Representative images of Sox2 and Nestin immunofluorescence in neurospheres; scale bar: 100 μm. (E) Representative images of neurosphere formation across 4 days in the two experimental groups; scale bar: 400 μm. (F) Top: quantification of neurospheres on D4 of passage in two groups (n = 3 mice per group). Bottom: measurement of average diameters of neurospheres on D4 of passage in two groups (n = 3 mice per group). (G) Accumulated cell number of NSC/NPCs from P1 to P4 in the two experimental groups (n = 3 mice per group). (H) Representative images of BrdU and Nestin immunofluorescence in NSC/NPCs; scale bar: 100 μm. Arrows show BrdU+Nestin+ cells. (I) Quantification of BrdU+Nestin+ cells (%) in NSC/NPCs (n = 3 mice per group). (J) Representative images of Tuj1 immunofluorescence in neurons differentiated in vitro in the two experimental groups; scale bar: 100 μm. (K) Quantification of Tuj1 positive cells (%) in the two experimental groups (n = 3 mice per group). Error bars represent the standard error of the mean. Significance was determined by Student’s t-test. *p < 0.05; **p < 0.01; and ****p < 0.0001.
NSC/NPCs from two groups were induced to differentiate into neurons, and the neuronic marker Tuj1 (Lee et al., 1990) was stained for quantification. The results showed neurons formed in both groups after a 10-day induction (Figure 1J); however, in contrast with proliferation assay, the number of neurons significantly increased in HE NSC/NPCs (Figure 1K).
Transcriptional Analysis of NSC/NPCs Reveals Hub mRNAs Involved in Neurogenesis
To elucidate the transcriptional changes related to altered neurogenesis, we compared the transcriptional profile of HE and NC NSC/NPCs by RNA-Seq. Heatmap of mRNAs showed distinctly different clustering between NSC/NPCs from two groups (Figure 2A). mRNAs with a log2 (fold change) ≥1 or ≤−1 and FDR < 0.05 were considered DE mRNAs, the volcano plot presented a total of 117 DE mRNAs, including 45 upregulated and 72 downregulated in HE NSC/NPCs compared with NC (Figure 2B). We conducted GSEA afterward, aiming to find neurogenesis-related gene sets (Figures 2C–F). The results revealed that genes downregulated in HE NSC/NPCs were enriched in neurogenesis-related GO biological processes (BP), such as neuroblast division, neuroblast proliferation, stem cell division, and neuron fate commitment, while genes upregulated in HE NSC/NPCs were enriched in response to estradiol (Figure 2D). The enrichment plots of these gene sets are presented in Figure 2G.
Figure 2. mRNA profiling and Hub mRNA identification. (A) Heatmap of 117 DE mRNAs in HE NSC/NPCs compared to NC. Red indicates upregulation, and blue indicates downregulation; row scale is from 0 to 1. (B) Volcano plot of DE mRNAs in HE NSC/NPCs compared with NC, red dots represent 45 upregulated mRNAs, and blue dots represent 72 downregulated mRNAs. (C) Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment results from gene set enrichment analysis (GSEA). The x-axis represents normalized enrichment score (NES), and the y-axis represents KEGG terms. The size of the dot indicates gene count, and the color indicates normalized value of p. Positive and negative NES indicate upregulation and downregulation in HE. (D-F) GO BP, CC (cellular component), and MF (molecular function) enrichment results from GSEA. The x-axis represents NES, and the y-axis represents GO terms. The size of the dot indicates gene count, and the color indicates normalized value of p. Positive and negative NESs indicate upregulation and downregulation in HE, respectively. (G) The enrichment plot of gene sets involved in neurogenesis and estradiol response. (H) Protein-protein interaction analysis of genes with core enrichment and Hub mRNA identification. Colored nodes indicate Hub mRNAs, with their shade positively correlated with the Maximal Clique Centrality score.
To further explore key mRNAs in the gene sets above and their interactions, mRNAs with core enrichment in each gene set (found in GSEA details) were picked for PPI analysis in STRING followed by Hub gene identification using Cytoscape. The top 10 genes ranked by MCC score were identified as Hub mRNAs, including Tbr1, Six3, Foxg1, Pou3f2, Dlx2, Fezf2, Dlx1, Nkx2-1, Sox1, and Notch1 (Figure 2H).
To validate the mRNA profiling and Hub gene identification results, the expression of the Hub mRNAs above was screened by quantitative real-time polymerase chain reaction with NSC/NPC samples used in RNA-Seq (six NC mice and six HE mice). The relative gene expression indicated all 10 Hub mRNAs decreased in the HE NSC/NPCs compared to the NC group and were identical to expression trends in RNA-Seq (Figure 3).
Figure 3. Verification of selected mRNAs by qPCR. Fold change of 10 Hub mRNAs in HE NSC/NPCs compared with NC (n = 6 mice per group). Significance was determined by Student’s t-test. *p < 0.05; **p < 0.01; and ***p < 0.001.
Identification of DE lncRNA-Hub mRNA Interaction Network
Epigenetic modification is recognized to regulate early life neurodevelopment (LaSalle et al., 2013; Yao et al., 2016), and dysregulation of lncRNAs leads to impaired development or neural dysfunction (Ng et al., 2013); therefore, we investigated whether lncRNA profiles were affected in NSC/NPCs with prenatal high estradiol stimulation. The resultant heatmap showed separated clustering of lncRNA transcripts in two groups (Figure 4A). Transcripts with a log2 (fold change) ≥1 or ≤−1 and FDR < 0.05 were considered DE lncRNA transcripts, and the volcano plot showed a total of 85 DE lncRNA transcripts, including 58 upregulated and 27 downregulated transcripts in the HE NSC/NPCs compared with the NC group (Figure 4B). The correlation between these DE lncRNA transcripts and the Hub mRNAs identified above were evaluated by constructing an expression matrix and calculating the PCC of each gene pair (Figure 4C). Gene pairs with a PCC ≥0.990 or ≤−0.990 and p < 0.05 were considered potentially correlated (Figure 4D). The co-expression network of these correlated genes was constructed (Figure 4E), including 6 lncRNA transcripts (ENSMUST00000037953, ENSMUST00000136217, ENSMUST00000138077, ENSMUST00000145804, ENSMUST 00000170557, and ENSMUST00000189763) and 3 Hub mRNAs (Sox1, Fezf2, Foxg1). Among them, ENSMUST00000189763, ENSMUST00000037953, and ENSMUST00000145804 were upregulated and the rest were downregulated. Therefore, we estimated these lncRNA-mRNA interactions may play a part in mediating the less proliferative and more neurogenic potential of fetal hypothalamic NSC/NPCs resulting from high maternal estradiol exposure.
Figure 4. LncRNA profiling and DE lncRNA-Hub mRNA interaction network construction. (A) Heatmap of 85 DE lncRNA transcripts in HE NSC/NPCs compared to NC. Red indicates upregulation, and blue indicates downregulation; row scale is from 0 to 1. (B) Volcano plot of DE lncRNA transcripts in HE NSC/NPCs compared to NC. Red dots represent 58 upregulated transcripts, and blue dots represent 27 downregulated transcripts. (C) Heatmap of DE lncRNA-Hub mRNA Pearson’s correlation coefficient (PCC) score. Orange indicates positive PCC, and blue indicates negative PCC. PCC ≥0.990 or ≤−0.990 and p < 0.05 are labeled with ⋆. (D) List of correlated DE lncRNA transcript-Hub mRNA pairs, with their PCC and value of p. (E) The interaction network of six DE lncRNA transcripts and three Hub mRNAs. Red indicates upregulation, and blue indicates downregulation in HE NSC/NPCs compared to NC. Triangles represent lncRNAs, and rectangles represent mRNAs.
Discussion
Mounting evidence suggests that an adverse intrauterine hormonal environment could impair the health of offspring. High maternal estradiol is usually induced by ovarian stimulation in assisted reproductive technology, and it can disrupt neurodevelopment, resulting in metabolic disorders and diminished verbal abilities in offspring (Wang et al., 2018; Zhou et al., 2020). In our previous study, high maternal estradiol led to insulin resistance and disordered eating in mouse offspring due to decreased insulin receptor and elevated neuropeptide Y expression in the hypothalamus (Wang et al., 2018). That these alterations were detected postnatally prompted us to search for corresponding events in earlier developmental stages.
Because fetal NSC/NPCs proliferate and differentiate actively, they are vulnerable to exogenous stimulators like high maternal estradiol, which can pass through the placental barrier (Gude et al., 2004) and may directly affect their biological properties. We isolated primary neural stem cells from the fetal hypothalamus and cultured them as neurospheres. Actually, neurospheres may derive from mixed cells with dynamic markers, and even the purified cells remain heterogeneous (Pastrana et al., 2011), resulting in limitations in neurosphere assay if applied alone. Therefore, we examined proliferation in Sox2 and Nestin-positive cells both in vivo and in vitro, which are recognized to be NSCs and NPCs (Bani-Yaghoub et al., 2006; Ernst and Christie, 2006; Hernandez et al., 2007; Campbell et al., 2015). Since the vast majority of the primary cells were Sox2 and/or Nestin positive, as shown in Figures 1D,H, we consider the researches were performed mainly in the same population of cells in two groups. Although a small proportion of other type of cells might exist, their effect on gene expression appeared insignificant when compared with NCS/NPCs. In consideration of the limitations of our study, we believe single-cell RNA-Seq should make a more precise method in future research. In fact, the cell composition in neurospheres would possibly change during different passages, and we speculate that the diminished proliferation might be a continuation of properties in primary NSC/NPCs, and further studies are expected to reveal whether increased cell apoptosis and senescence could occur.
The proliferation medium contained EGF and bFGF, which are essential factors for NSC/NPC growth. To explore whether cell proliferation changed due to different levels of EGF and bFGF receptors, we examined the mRNA expressions of them in primary NSC/NPCs and found no significant difference between the two groups (Supplementary Figure S1). This result indicated that the decreased proliferation in HE NSC/NPCs was not attributed to different levels of growth factor receptors in the culture medium, but more likely to the intrauterine programming effect.
Neurogenesis is a complex process, and how NSC/NPCs are maintained, divided, and differentiated remain controversial (Lazutkin et al., 2019). Our results indicated decreased proliferative activity and increased neuronal production in HE NSC/NPCs; however, whether or not this resulted from premature exhaustion of the stem cell pool requires further study. A comprehensive evaluation of hypothalamic neurogenesis from the prenatal period to adulthood may show us a more precise effect of high maternal estradiol on neurodevelopment. It should be noted that the NSC/NPC differentiation assay was carried out in vitro, which generates early and immature neurons rather than functional neurons, and an in vivo labeling of NSC/NPCs at the embryonic stage followed by detections like lineage tracing may provide a more accurate indication of their differentiation directions.
Several published studies show that estrogen stimulates both proliferation and differentiation of NSC/NPCs (Okada et al., 2010; Li et al., 2020) and attenuates damage to neurogenesis in the developing brain caused by chemical drug exposure (Li et al., 2019); however, another study shows that 10 nM estradiol increases NSC/NPCs proliferation and stimulates differentiation into neurons in vitro, but 50 nM estradiol markedly decreases NSC/NPCs proliferation (Zhang et al., 2019a). Thus, the effects of maternal estradiol on neurogenesis of fetal hypothalamic NSC/NPCs may be dose dependent, and the dose that caused metabolic disorder in our mouse model exerted a different effect on proliferation and differentiation. The stimulative effect on neuron formation may explain our previous finding that hypothalamic neuropeptide Y increases in HE offspring (Wang et al., 2018); that is, prenatal high estradiol probably promotes orexigenic neuron generation, leading to disordered eating.
We used log2 (fold change) ≥1 or ≤−1 and FDR < 0.05 as the cutoffs to define DE genes in our study, and looser criteria of log2 (fold change) ≥1 or ≤−1 with p < 0.05 were also tried to identify DE genes and predict DE lncRNA-Hub mRNA regulatory network. This method presented 567 DE mRNAs (383 upregulated and 184 downregulated) and 148 DE lncRNA transcripts (89 upregulated and 59 downregulated) in HE NSC/NPCs compared with NC group (Supplementary Figures S2A–D). In spite of this, the potentially correlated DE lncRNA-Hub mRNA pairs remained unchanged compared with those marked in Figure 4C (Supplementary Figure S2E).
Gene set enrichment analysis (GSEA) revealed both enrichment and expression of mRNA profiles in two groups; positive NES indicated upregulation in HE compared with NC NSC/NPCs, and negative NES indicated downregulation. Although KEGG enrichment did not reveal items directly involved in cell proliferation and differentiation, it showed differential enrichment of cell cycles, DNA replication, nucleotide excision repair, and RNA polymerase. Moreover, we found upregulated mRNAs that are enriched in Huntington’s and Parkinson’s disease in the HE NSC/NPCs, suggesting an increased risk of neurodegenerative disease in offspring exposed to high maternal estradiol, which requires further validation. The GO CC and MF enrichment revealed downregulated genes enriched in autophagosome and dopamine receptor binding, and upregulated genes enriched in neurotransmitter binding and glutamate receptor binding in HE NSC/NPCs. These discoveries may help illuminate the mechanisms of maternal estradiol-induced neurodevelopmental disorders.
The GO BP enrichment result revealed downregulated genes in HE NSC/NPCs enriched in stem cell division, neuroblast proliferation, and neuroblast division, which consisted with their less proliferative potential observed. However, genes enriched in neuron fate commitment also decreased in HE NSC/NPCs, which seemingly contradicted the more neurogenic activity in vitro. To find a rational explanation for this, we focused on Hub mRNAs of core enrichment in these gene sets.
The Hub mRNAs were identified based on a PPI network consisting of genes enriched in neurogenesis and response to estradiol; however, the top 10 Hub mRNAs ranked by MCC in Cytoscape were all genes involved in neurogenesis, and they were all downregulated in HE NSC/NPCs. This result may be explained by the fact that Hub genes are highly connected genes in a co-expression network, and genes enriched in response to estradiol failed to present such close connections with those enriched in altered neurogenesis, implying a probable indirect effect of estradiol stimulation on neurogenesis in our study. Most of these Hub mRNAs are transcription factors except Notch1. Pou3f2 influences multiple stages of neurogenesis by promoting neural transcription factor Tbr1 (Dominguez et al., 2013), which regulates cell differentiation and migration and involves glutamatergic neurogenesis (Mihalas and Hevner, 2017). Six3, Foxg1, and Sox1 maintain the balance between proliferation and neuronal differentiation of NSC/NPCs. Specifically, upregulation of Six3 plays a role in keeping NSC/NPCs in an undifferentiated state (Appolloni et al., 2008); Foxg1 deficiency leads to premature differentiation of neurons, and its overexpression increases the NSC/NPC pool (Hanashima et al., 2004; Brancaccio et al., 2010); and Sox1 loss induces depletion of proliferating NSC/NPCs with increased cell cycle exit (Bylund et al., 2003). Dlx1 and Dlx2 drive GABAergic neuron generation (Lindtner et al., 2019; Barretto et al., 2020), and Fezf2 is involved in the dopaminergic neuron generation (Eckler and Chen, 2014); moreover, knockdown of Fezf2 leads to decreased Foxg1 and Six3 in mouse embryonic stem cells (Wang et al., 2011). Notch1 signaling is reported to promote NSC/NPC proliferation but decrease neuronal differentiation during meningitis and spinal cord injury (Peng et al., 2019; Zhang et al., 2019b). Nkx2-1 is a critical factor maintaining the anorectic gene Pomc expression from early development to adulthood (Orquera et al., 2019). To sum up, the published information above supports our findings that decreased Six3, Foxg1, Sox1, Fezf2, and Notch1 in HE NSC/NPCs directly correlated with decreased proliferation and enhanced neuronal generation, which probably reflected a premature differentiation. As the upstream regulator of Pomc, decreased Nkx2-1 in HE offspring could contribute to the verified orexigenic phenotype in later life (Wang et al., 2018). Since Pou3f2, Tbr1, Dlx1, and Dlx2 are involved in the generation of glutamatergic or dopaminergic or GABAergic neurons, figuring out the neuron types that these NSC/NPCs tended to form would help validate the effect on neurogenesis of these genes in our experiment.
The published studies fail to specify the effect of estradiol on Hub mRNAs above during neurogenesis. One research shows estradiol stimulation does not affect Notch1 expression during hippocampus development but reduces the level of its transcriptionally active domain (Bender et al., 2010). Since we previously found DNA methylation programs hypothalamic gene expression in HE offspring (Wang et al., 2018), it may support the hypothesis that expression changes of these Hub mRNAs could be attributed to epigenetic regulators, such as DNA methylation and non-coding RNAs.
LncRNAs are important components of regulatory networks in CNS development whose dysregulation leads to neurological disorders (Ng et al., 2013). Since lncRNAs exert functions mainly through regulating mRNA expression, we constructed the interaction network between DE lncRNAs and Hub mRNAs, aiming to discover lncRNAs possibly connected with the altered neurogenesis. Our study predicted six DE lncRNA transcripts correlated with three Hub mRNAs based on their expression level. These lncRNA transcripts (ENSMUST 00000037953, ENSMUST00000136217, ENSMUST 00000138077, ENSMUST00000145804, ENSMUST00000170557, and ENSMUST00000189763) are coded by genes 2810032 G03Rik, Prdm16os, Gm13110, Ppp1r18os, Gm17035, and D130058E05Rik, respectively, according to Ensemble genome browser.4 We also checked the expression correlations between DE lncRNAs and all mRNAs in Figure 2H (mRNAs with core enrichment in neurogenesis and response to estradiol, core mRNAs for short), which predicted a potential regulatory network of 14 DE lncRNA transcripts and 11 core mRNAs (Supplementary Figures S3A–C). In consideration of our limited sample size, we repeated the co-expression analysis of these lncRNAs and mRNAs using the pubic dataset GSE65487 in Gene Expression Omnibus, which assessed the RNA profiles of proliferating progenitors, differentiating progenitors and neurons from E14.5 mouse cortex. However, five of DE lncRNAs were not detected in GSE65487 (Supplementary Figure S3D), and four of the rest presented potential correlations with three mRNAs each, when the cutoffs were set to be PCC ≥0.900 or ≤−0.900 and p < 0.05 (PCC ≥0.990 or ≤−0.990 identified no significant correlation; Supplementary Figures S3E,F). Although we failed to discover the same lncRNA-mRNA pair as demonstrated in our study by using this public dataset, it still reflected possible connections of these lncRNAs with neurogenesis. Since there are not yet any published literatures about the functions of these lncRNAs, future work focused on their specific roles is expected to help answer questions regarding development-originated neuroendocrine disorders.
In short, our research presents the cytologic changes in early neural development under a high maternal estradiol environment and reveals the corresponding whole genomic features with a prediction of the underlying molecular modifications. This study demonstrates comprehensive information about fetal hypothalamic NSC/NPCs with prenatal high estradiol exposure and contributes to our understanding of the fetal-programmed adult diseases.
Data Availability Statement
The datasets presented in this study can be found in an online repository. The name of the repository and accession number can be found at https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE168075.
Ethics Statement
This animal study was reviewed and approved by the Institutional Animal Care and Use Committee of Shanghai Jiao Tong University.
Author Contributions
HW and CZ designed the experiments, collected the data, analyzed the data, and drafted the manuscript. MH, HH, and YS revised the final manuscript. All authors have reviewed the manuscript before submitting it and approved the final version.
Funding
This work was supported by the National Natural Science Foundation of China (Grant Nos. 81901552 and 82001650).
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.
Supplementary Material
The Supplementary Material for this article can be found online at https://www.frontiersin.org/articles/10.3389/fgene.2021.677935/full#supplementary-material
Supplementary Figure S1 | mRNA expression of EGF and bFGF receptors examined by qPCR. Fold change of EGF and bFGF receptor mRNAs in HE NSC/NPCs compared with NC (n = 6 mice per group). Significance was determined by Student’s t-test; ns, not significant.
Supplementary Figure S2 | DE gene identification and DE lncRNA-Hub mRNA interaction network analysis using log2 (fold change) ≥1 or ≤−1 with p < 0.05 as the cutoff. (A) Heatmap of 567 DE mRNAs in HE NSC/NPCs compared to NC. Red indicates upregulation, and blue indicates downregulation; row scale is from 0 to 1. (B) Volcano plot of DE mRNAs in HE NSC/NPCs compared to NC, red dots represent 383 upregulated mRNAs, and blue dots represent 184 downregulated mRNAs. (C) Heatmap of 148 DE lncRNA transcripts in HE NSC/NPCs compared to NC. Red indicates upregulation, and blue indicates downregulation; row scale is from 0 to 1. (D) Volcano plot of DE lncRNA transcripts in HE NSC/NPCs compared with NC. Red dots represent 89 upregulated transcripts, and blue dots represent 59 downregulated transcripts. (E) Heatmap of DE lncRNA-Hub mRNA PCC score. Orange indicates positive PCC, and blue indicates negative PCC. PCC ≥0.990 or ≤−0.990 and p < 0.05 are labeled with ⋆.
Supplementary Figure S3 | DE lncRNA-core mRNA interaction network construction and verification with public dataset. (A) Heatmap of DE lncRNA-core mRNA PCC score. Orange indicates positive PCC, and blue indicates negative PCC. PCC ≥0.990 or ≤−0.990 and p < 0.05 are labeled with ⋆ (DE lncRNA-Hub mRNA) or ● (DE lncRNA-non Hub mRNA). (B) List of correlated DE lncRNA transcript-core mRNA pairs, with their PCC and value of p. (C) The interaction network of 14 DE lncRNA transcripts and 11 core mRNAs. Red indicates upregulation, and blue indicates downregulation in HE NSC/NPCs compared to NC. Triangles represent lncRNAs, and rectangles represent mRNAs. (D) List of lncRNAs in panel B with the corresponding lncRNA genes found in GSE65487. (E) Heatmap of lncRNA-mRNA PCC score graphed using data from GSE65487. Orange indicates positive PCC, and blue indicates negative PCC. PCC ≥0.900 or ≤−0.900 and p < 0.05 are labeled with ⋆. (F) List of correlated lncRNA-mRNA pairs presented in panel E, with their PCC and value of p.
Footnotes
References
Appolloni, I., Calzolari, F., Corte, G., Perris, R., and Malatesta, P. (2008). Six3 controls the neural progenitor status in the murine CNS. Cereb. Cortex 18, 553–562. doi: 10.1093/cercor/bhm092
Bani-Yaghoub, M., Tremblay, R. G., Lei, J. X., Zhang, D., Zurakowski, B., Sandhu, J. K., et al. (2006). Role of Sox2 in the development of the mouse neocortex. Dev. Biol. 295, 52–66. doi: 10.1016/j.ydbio.2006.03.007
Barretto, N., Zhang, H., Powell, S. K., Fernando, M. B., Zhang, S., Flaherty, E. K., et al. (2020). ASCL1- and DLX2-induced GABAergic neurons from hiPSC-derived NPCs. J. Neurosci. Methods 334:108548. doi: 10.1016/j.jneumeth.2019.108548
Bateson, P., Barker, D., Clutton-Brock, T., Deb, D., D'Udine, B., Foley, R. A., et al. (2004). Developmental plasticity and human health. Nature 430, 419–421. doi: 10.1038/nature02725
Bender, R. A., Zhou, L., Wilkars, W., Fester, L., Lanowski, J. S., Paysen, D., et al. (2010). Roles of 17ss-estradiol involve regulation of reelin expression and synaptogenesis in the dentate gyrus. Cereb. Cortex 20, 2985–2995. doi: 10.1093/cercor/bhq047
Bond, A. M., Ming, G. L., and Song, H. (2015). Adult mammalian neural stem cells and neurogenesis: five decades later. Cell Stem Cell 17, 385–395. doi: 10.1016/j.stem.2015.09.003
Brancaccio, M., Pivetta, C., Granzotto, M., Filippis, C., and Mallamaci, A. (2010). Emx2 and Foxg1 inhibit gliogenesis and promote neuronogenesis. Stem Cells 28, 1206–1218. doi: 10.1002/stem.443
Bylund, M., Andersson, E., Novitch, B. G., and Muhr, J. (2003). Vertebrate neurogenesis is counteracted by Sox1-3 activity. Nat. Neurosci. 6, 1162–1168. doi: 10.1038/nn1131
Campbell, J. G., Miller, D. C., Cundiff, D. D., Feng, Q., and Litofsky, N. S. (2015). Neural stem/progenitor cells react to non-glial CNS neoplasms. Springerplus 4:53. doi: 10.1186/s40064-015-0807-z
Chin, C. H., Chen, S. H., Wu, H. H., Ho, C. W., Ko, M. T., and Lin, C. Y. (2014). cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst. Biol. 8(Suppl. 4):S11. doi: 10.1186/1752-0509-8-S4-S11
Desai, M., Li, T., and Ross, M. G. (2011). Fetal hypothalamic neuroprogenitor cell culture: preferential differentiation paths induced by leptin and insulin. Endocrinology 152, 3192–3201. doi: 10.1210/en.2010-1217
Dominguez, M. H., Ayoub, A. E., and Rakic, P. (2013). POU-III transcription factors (Brn1, Brn2, and Oct6) influence neurogenesis, molecular identity, and migratory destination of upper-layer cells of the cerebral cortex. Cereb. Cortex 23, 2632–2643. doi: 10.1093/cercor/bhs252
Eckler, M. J., and Chen, B. (2014). Fez family transcription factors: controlling neurogenesis and cell fate in the developing mammalian nervous system. BioEssays 36, 788–797. doi: 10.1002/bies.201400039
Ernst, C., and Christie, B. R. (2006). The putative neural stem cell marker, nestin, is expressed in heterogeneous cell types in the adult rat neocortex. Neuroscience 138, 183–188. doi: 10.1016/j.neuroscience.2005.10.065
Fuentealba, L. C., Rompani, S. B., Parraguez, J. I., Obernier, K., Romero, R., Cepko, C. L., et al. (2015). Embryonic origin of postnatal neural stem cells. Cell 161, 1644–1655. doi: 10.1016/j.cell.2015.05.041
Gilyarov, A. V. (2008). Nestin in central nervous system cells. Neurosci. Behav. Physiol. 38, 165–169. doi: 10.1007/s11055-008-0025-z
Gude, N. M., Roberts, C. T., Kalionis, B., and King, R. G. (2004). Growth and function of the normal human placenta. Thromb. Res. 114, 397–407. doi: 10.1016/j.thromres.2004.06.038
Hanashima, C., Li, S. C., Shen, L., Lai, E., and Fishell, G. (2004). Foxg1 suppresses early cortical cell fate. Science 303, 56–59. doi: 10.1126/science.1090674
Hernandez, P. P., Olivari, F. A., Sarrazin, A. F., Sandoval, P. C., and Allende, M. L. (2007). Regeneration in zebrafish lateral line neuromasts: expression of the neural progenitor cell marker sox2 and proliferation-dependent and-independent mechanisms of hair cell renewal. Dev. Neurobiol. 67, 637–654. doi: 10.1002/dneu.20386
Hu, X. L., Feng, C., Lin, X. H., Zhong, Z. X., Zhu, Y. M., Lv, P. P., et al. (2014). High maternal serum estradiol environment in the first trimester is associated with the increased risk of small-for-gestational-age birth. J. Clin. Endocrinol. Metab. 99, 2217–2224. doi: 10.1210/jc.2013-3362
LaSalle, J. M., Powell, W. T., and Yasui, D. H. (2013). Epigenetic layers and players underlying neurodevelopment. Trends Neurosci. 36, 460–470. doi: 10.1016/j.tins.2013.05.001
Lazutkin, A., Podgorny, O., and Enikolopov, G. (2019). Modes of division and differentiation of neural stem cells. Behav. Brain Res. 374:112118. doi: 10.1016/j.bbr.2019.112118
Lee, M. K., Tuttle, J. B., Rebhun, L. I., Cleveland, D. W., and Frankfurter, A. (1990). The expression and posttranslational modification of a neuron-specific beta-tubulin isotype during chick embryogenesis. Cell Motil. Cytoskeleton 17, 118–132. doi: 10.1002/cm.970170207
Li, J., Tang, Y., and Cai, D. (2012). IKKbeta/NF-kappaB disrupts adult hypothalamic neural stem cells to mediate a neurodegenerative mechanism of dietary obesity and pre-diabetes. Nat. Cell Biol. 14, 999–1012. doi: 10.1038/ncb2562
Li, J., Tang, Y., Purkayastha, S., Yan, J., and Cai, D. (2014). Control of obesity and glucose intolerance via building neural stem cells in the hypothalamus. Mol. Metab. 3, 313–324. doi: 10.1016/j.molmet.2014.01.012
Li, W., Li, H., Wei, H., Lu, Y., Lei, S., Zheng, J., et al. (2019). 17beta-estradiol treatment attenuates neurogenesis damage and improves behavior performance after ketamine exposure in neonatal rats. Front. Cell. Neurosci. 13:251. doi: 10.3389/fncel.2019.00251
Li, W., Lu, P., Lu, Y., Wei, H., Niu, X., Xu, J., et al. (2020). 17beta-estradiol protects neural stem/progenitor cells against ketamine-induced injury through estrogen receptor beta pathway. Front. Neurosci. 14:576813. doi: 10.3389/fnins.2020.576813
Lindtner, S., Catta-Preta, R., Tian, H., Su-Feher, L., Price, J. D., Dickel, D. E., et al. (2019). Genomic resolution of DLX-orchestrated transcriptional circuits driving development of forebrain GABAergic neurons. Cell Rep. 28, 2048.e8–2063.e8. doi: 10.1016/j.celrep.2019.07.022
Livesey, F. J. (2012). A potential link between obesity and neural stem cell dysfunction. Nat. Cell Biol. 14, 987–989. doi: 10.1038/ncb2599
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8
Marshall, G. P. 2nd, Reynolds, B. A., and Laywell, E. D. (2007). Using the neurosphere assay to quantify neural stem cells in vivo. Curr. Pharm. Biotechnol. 8, 141–145. doi: 10.2174/138920107780906559
Mihalas, A. B., and Hevner, R. F. (2017). Control of neuronal development by T-box genes in the brain. Curr. Top. Dev. Biol. 122, 279–312. doi: 10.1016/bs.ctdb.2016.08.001
Ng, S. Y., Lin, L., Soh, B. S., and Stanton, L. W. (2013). Long noncoding RNAs in development and disease of the central nervous system. Trends Genet. 29, 461–468. doi: 10.1016/j.tig.2013.03.002
Okada, M., Makino, A., Nakajima, M., Okuyama, S., Furukawa, S., and Furukawa, Y. (2010). Estrogen stimulates proliferation and differentiation of neural stem/progenitor cells through different signal transduction pathways. Int. J. Mol. Sci. 11, 4114–4123. doi: 10.3390/ijms11104114
Orquera, D. P., Tavella, M. B., de Souza, F. S. J., Nasif, S., Low, M. J., and Rubinstein, M. (2019). The homeodomain transcription factor NKX2.1 is essential for the early specification of melanocortin neuron identity and activates Pomc expression in the developing hypothalamus. J. Neurosci. 39, 4023–4035. doi: 10.1523/JNEUROSCI.2924-18.2019
Padilla, S. L., Carmody, J. S., and Zeltser, L. M. (2010). Pomc-expressing progenitors give rise to antagonistic neuronal populations in hypothalamic feeding circuits. Nat. Med. 16, 403–405. doi: 10.1038/nm.2126
Pastrana, E., Silva-Vargas, V., and Doetsch, F. (2011). Eyes wide open: a critical review of sphere-formation as an assay for stem cells. Cell Stem Cell 8, 486–498. doi: 10.1016/j.stem.2011.04.007
Peng, Z., Li, X., Fu, M., Zhu, K., Long, L., Zhao, X., et al. (2019). Inhibition of Notch1 signaling promotes neuronal differentiation and improves functional recovery in spinal cord injury through suppressing the activation of Ras homolog family member A. J. Neurochem. 150, 709–722. doi: 10.1111/jnc.14833
Ralevski, A., and Horvath, T. L. (2015). Developmental programming of hypothalamic neuroendocrine systems. Front. Neuroendocrinol. 39, 52–58. doi: 10.1016/j.yfrne.2015.09.002
Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U. S. A. 102, 15545–15550. doi: 10.1073/pnas.0506580102
Suh, H., Consiglio, A., Ray, J., Sawai, T., D'Amour, K. A., and Gage, F. H. (2007). In vivo fate analysis reveals the multipotent and self-renewal capacities of Sox2+ neural stem cells in the adult hippocampus. Cell Stem Cell 1, 515–528. doi: 10.1016/j.stem.2007.09.002
Wang, H. H., Zhou, C. L., Lv, M., Yang, Q., Li, J. X., Hou, M., et al. (2018). Prenatal high estradiol exposure induces sex-specific and dietarily reversible insulin resistance through decreased hypothalamic INSR. Endocrinology 159, 465–476. doi: 10.1210/en.2017-03017
Wang, Z. B., Boisvert, E., Zhang, X., Guo, M., Fashoyin, A., Du, Z. W., et al. (2011). Fezf2 regulates telencephalic precursor differentiation from mouse embryonic stem cells. Cereb. Cortex 21, 2177–2186. doi: 10.1093/cercor/bhr006
Yao, B., Christian, K. M., He, C., Jin, P., Ming, G. L., and Song, H. (2016). Epigenetic mechanisms in neurogenesis. Nat. Rev. Neurosci. 17, 537–549. doi: 10.1038/nrn.2016.70
Zhang, L., Ma, Y., Liu, M., Ma, Y., and Guo, H. (2019a). The effects of various estrogen doses on the proliferation and differentiation of cultured neural stem cells. Gen. Physiol. Biophys. 38, 417–425. doi: 10.4149/gpb_2019022
Zhang, Z., Xu, D., Zhao, S., Lian, D., Wu, J., He, D., et al. (2019b). Notch1 signaling pathway promotes proliferation and mediates differentiation direction in hippocampus of streptococcus pneumonia meningitis rats. J. Infect. Dis. 220, 1977–1988. doi: 10.1093/infdis/jiz414
Keywords: prenatal exposure, estradiol, neural stem/progenitor cells, neurogenesis, RNA sequencing, gene set enrichment analysis, protein-protein interaction, interaction network
Citation: Wang H, Zhou C, Hou M, Huang H and Sun Y (2021) Neurogenesis Potential Evaluation and Transcriptome Analysis of Fetal Hypothalamic Neural Stem/Progenitor Cells With Prenatal High Estradiol Exposure. Front. Genet. 12:677935. doi: 10.3389/fgene.2021.677935
Edited by:
Carmem Gottfried, Federal University of Rio Grande do Sul, BrazilReviewed by:
Eduardo Bouth Sequerra, Federal University of Rio Grande do Norte, BrazilJubao Duan, University of Chicago, United States
Copyright © 2021 Wang, Zhou, Hou, Huang and Sun. 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: Hefeng Huang, aHVhbmdoZWZnQGZ1ZGFuLmVkdS5jbg==; Yun Sun, c3l1bjE2M0AxNjMuY29t
†These authors have contributed equally to this work and share first authorship