Single-cell profiling of lncRNAs in human germ cells and molecular analysis reveals transcriptional regulation of LNC1845 on LHX8

  1. Nan Wang
  2. Jing He
  3. Xiaoyu Feng
  4. Shengyou Liao
  5. Yi Zhao
  6. Fuchou Tang
  7. Kehkooi Kee  Is a corresponding author
  1. Center for Stem Cell Biology and Regenerative Medicine, Department of Basic Medical Sciences, School of Medicine, Tsinghua University, China
  2. Key Laboratory of Intelligent Information Processing, Advanced Computer Research Center, Institute of Computing Technology, Chinese Academy of Sciences, China
  3. Biodynamic Optical Imaging Center & Department of Obstetrics and Gynecology, College of Life Sciences, Third Hospital, Peking University, China
  4. Tsinghua-Peking Center for Life Sciences, School of Life Sciences, Tsinghua University, China
8 figures, 1 table and 10 additional files

Figures

Figure 1 with 1 supplement
Figure numerous.

(A) Long non-coding RNAs (lncRNAs) are differentially expressed during human germ cell (hGC) development. (B) Two unannotated lncRNAs and their neighboring protein-coding genes with their genomic …

Figure 1—figure supplement 1
Numerous long non-coding RNAs (lncRNAs) are differentially expressed during human germ cell development.

(A) Numbers of protein-coding genes expressed in individual female (upper) and male (bottom) germ cells or somatic cells from different developmental week. Student’s t-test was used for the …

Figure 2 with 1 supplement
Integrating long non-coding RNA (lncRNA) expression enhances cell-type classifications of human female germ cells.

(A) tSNE plot of female germ cells and somatic cells colored by identified cell types. Left: cluster according to protein-coding genes. Right: cluster according to protein-coding genes and lncRNAs. …

Figure 2—figure supplement 1
Integrating long non-coding RNA (lncRNA) expression enhances cell-type classifications of human female germ cells.

(A) Gene ontology biological processes enrichment of mitotic-4 cells feature genes. (B) Single-cell trajectories of female germ cell states through the pseudotime according to protein-coding genes …

Figure 3 with 1 supplement
Genomic distributions and biotypes of the long non-coding RNAs (lncRNAs) expressed during human gonadal development.

(A) The genomic position distribution of lncRNAs expressed in human germ cells (hGCs) and gonadal somatic cells. (B) The distribution and percentage of the six locus biotypes of each developmental …

Figure 3—figure supplement 1
Genomic distributions and biotypes of the long non-coding RNAs (lncRNAs) expressed during human gonadal development.

(A) The genomic position distribution of annotated and unannotated lncRNAs expressed in human germ cells (hGCs) and gonadal somatic cells. (B) Expression correlations of lncRNA–mRNA pairs of the six …

Figure 4 with 2 supplements
Characterization of a divergent long non-coding RNA (lncRNA), LNC1845, expressed during the meiotic stage and is required for normal expression of LHX8.

(A) Genomic locus and relative position of LNC1845 (in green) and LHX8 (in blue). (B) Expression analysis by RT-qPCR of PA-1 cells, undifferentiated H9 cells, and RNA sample extracted from human …

Figure 4—figure supplement 1
Characterization of a divergent long non-coding RNA (lncRNA), LNC1845, expressed during the meiotic stage and is required for normal expression of LHX8.

(A) LNC1845 was predicted to be non-coding RNA. The RNA sequences of LNC1845 were put into the Coding Potential Calculator (CPC) program and the Coding Potential Assessment Tool. (B) In vitro …

Figure 4—figure supplement 1—source data 1

Characterization of a divergent long non-coding RNA (lncRNA), LNC1845, expressed during the meiotic stage and is required for normal expression of LHX8.

https://cdn.elifesciences.org/articles/78421/elife-78421-fig4-figsupp1-data1-v2.zip
Figure 4—figure supplement 2
Upregulated XH long non-coding RNA (lncRNA)–mRNA pairs in hPGCLCs.

Relative position of upregulated divergent lncRNA–mRNA pairs and their relative expression in primordial germ cell-like cells (PGCLCs) and control cells.

Figure 5 with 1 supplement
LNC1845 RNA transcripts are required for normal expression of LHX8.

(A) Schematic diagram of knockout strategies at the LNC1845 loci. (B) Expression analysis by RT-qPCR of WT primordial germ cell-like cells (PGCLCs) (n = 3) and LNC1845 KO PGCLCs (n = 3). (C) Western …

Figure 5—source data 1

Western blot analysis of LHX8 in PGCLCs from WT, LNC1845 KO, or RNA interfering LNC1845 (KD) hESCs. SD, spontaneous differentiation cells as the negative control.

https://cdn.elifesciences.org/articles/78421/elife-78421-fig5-data1-v2.zip
Figure 5—source data 2

Western blot analysis of GAPDH in PGCLCs from WT, LNC1845 KO, or RNA interfering LNC1845 (KD) hESCs. SD, spontaneous differentiation cells as the negative control.

https://cdn.elifesciences.org/articles/78421/elife-78421-fig5-data2-v2.zip
Figure 5—figure supplement 1
LNC1845 RNA transcripts are required for normal expression of LHX8.

(A) Genotyping and sequencing of LNC1845 KO single colonies to confirm the deletion of LNC1845. (B) Principle component analysis (PCA) of RNA-seq datasets from ESCs, LNC1845 KO primordial germ …

Figure 5—figure supplement 1—source data 1

Genotyping of LNC1845 KO single colonies to confirm the deletion of LNC1845.

https://cdn.elifesciences.org/articles/78421/elife-78421-fig5-figsupp1-data1-v2.zip
Figure 5—figure supplement 1—source data 2

Genotyping of LNC1845-upstream 3×polyA insertion single colonies.

https://cdn.elifesciences.org/articles/78421/elife-78421-fig5-figsupp1-data2-v2.zip
Figure 6 with 2 supplements
LNC1845 regulates LHX8 expression in cis by changing chromatin modifications.

(A) Schematic diagram of LNC1845-upstream EF1α promoter knock-in and the targeted allele after EF1α promoter insertion. (B) Expression analysis by RT-qPCR of LNC1845 and LHX8 in WT primordial germ …

Figure 6—source data 1

Western blot analysis of LHX8 in PGCLCs from WT, LNC1845-upstream 3×polyA knock-in, or LNC1845-upstream EF1α promoter knock-in cells.

https://cdn.elifesciences.org/articles/78421/elife-78421-fig6-data1-v2.zip
Figure 6—source data 2

Western blot analysis of GAPDH in PGCLCs from WT, LNC1845-upstream 3×polyA knock-in, or LNC1845-upstream EF1α promoter knock-in cells.

https://cdn.elifesciences.org/articles/78421/elife-78421-fig6-data2-v2.zip
Figure 6—figure supplement 1
LNC1845 regulates LHX8 expression in cis by changing chromatin modifications.

(A) Genotyping and sequencing of LNC1845-upstream EF1α promoter insertion single colonies. (B) Karyotype analysis of LNC1845 KO human embryonic stem cell (hESC) clone F11, LNC1845 3xpolyA knock in …

Figure 6—figure supplement 1—source data 1

Genotyping of LNC1845-upstream EF1α promoter insertion single colonies.

https://cdn.elifesciences.org/articles/78421/elife-78421-fig6-figsupp1-data1-v2.zip
Figure 6—figure supplement 1—source data 2

Agarose gel analysis of LNC1845 transcripts by in vitro transcription.

https://cdn.elifesciences.org/articles/78421/elife-78421-fig6-figsupp1-data2-v2.zip
Figure 6—figure supplement 2
LNC1845 regulates LHX8 expression in cis by physically interacting with WDR5.

(A) WDR5 expression level in single hPGC or gonadal somatic cells at different developmental stages. (B) Immunofluorescence of WDR5 (green) co-stained with DAZL (red) in 14-week ovary section …

Figure 6—figure supplement 2—source data 1

LNC1845 binding for WDR5 was verified by CHIRP-western blot.

https://cdn.elifesciences.org/articles/78421/elife-78421-fig6-figsupp2-data1-v2.zip
Figure 7 with 1 supplement
FOXP3 upregulates expression of LNC1845 and other long non-coding RNAs (lncRNAs).

(A) RT-qPCR analysis of LNC1845 and LHX8 in primordial germ cell-like cells (PGCLCs) overexpresses FOXP3. (B) Heatmap of differentially expressed ms-lncRNAs in control or OE FOXP3 groups, rep1 and …

Figure 7—figure supplement 1
FOXP3 upregulates expression of LNC1845 and other long non-coding RNAs (lncRNAs).

(A) Schematic diagram of input features for the computational prediction used to predict transcription factors involved in female meiotic stage-specific lncRNAs upregulation. (B) Bar chart …

Author response image 1
Integrative genomic viewer images showing FOXP3 binding sites at the LC1845/LHX8 locus.

Tables

Appendix 1—key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Gene (Homo sapiens)LIM homebox 8(LHX8)Human genome
informatics
Ensembl:ENSG
00000162624
Gene (Homo sapiens)LNC1845(AC099786.3)Human genome
informatics
NONHSAG001845.
2/ENSG00000261213.1
Cell line (Homo sapiens)H9 LNC1845 KO stem cellsThis paperLNC1845 KO cell line generated by CRISPR/Cas9
Cell line (Homo sapiens)H9 LNC1845 polyA
KI stem cells
This paperLNC1845 polyA KI cell line generated by CRISPR/Cas9
Cell line (Homo sapiens)H9 LNC1845 EF1α
KI stem cells
This paperLNC1845 EF1α KI cell line generated by CRISPR/Cas9
Cell line (Homo sapiens)PA-1 cellsThis paperHuman ovary teratocarcinoma Cell Line
Cell line (Homo sapiens)293 FT cellsThermo Fisher ScientificCat# R70007
Cell line (Mus musculus)Mouse embryonic fibroblastsThis paperMouse embryonic fibroblasts isolated from embryonic day 13.5
Recombinant DNA reagentP2k7-EF1α-DAZL2-
P2A-mCherry
This studyMethods
Recombinant DNA reagentP2k7-EF1α-mCherryThis studyMethods
Recombinant DNA reagentP2k7-EF1α-1845This studyMethods
Recombinant DNA reagentP2k7-EF1α-FOXP3This studyMethods
Recombinant DNA reagentP2k7-EF1α-TFAP2AThis studyMethods
Recombinant DNA reagentP2k7-EF1α-FOXP3-
P2A-mCherry
This studyMethods
Recombinant DNA reagentP2k7-EF1α-FOXP3-
3xFLAG-P2A-mCherry
This studyMethods
Recombinant DNA reagentpEASY-T7 promoterThis studyMethods
Recombinant DNA reagentpEASY-T7-1845This studyMethods
Recombinant DNA reagentpX335-U6-Chimeric_BB-CBh-hSpCas9n(D10A)Addgene# 42335
Recombinant DNA reagentPGK-puro-1845KO donorThis studyMethods
Recombinant DNA reagentSV40-puro-3xpolyA 1845 KI donorThis studyMethods
Recombinant DNA reagentSV40-puro-EF1α 1845 KI donorThis studyMethods
Recombinant DNA reagentH1-shLACZ-Ubc-GFPThis studyMethods
Recombinant DNA reagentH1-sh1845①-Ubc-GFPThis studyMethods
Recombinant DNA reagentH1-sh1845②-Ubc-GFPThis studyMethods
Recombinant DNA reagentH1-sh15266①-Ubc-GFPThis studyMethods
Recombinant DNA reagentH1-sh15266②-Ubc-GFPThis studyMethods
Recombinant DNA reagentH1-sh3346①-Ubc-GFPThis studyMethods
Recombinant DNA reagentH1-sh3346②-Ubc-GFPThis studyMethods
Recombinant DNA reagentH1-shFOXP3①-Ubc-GFPThis studyMethods
Recombinant DNA reagentH1-shFOXP3②-Ubc-GFPThis studyMethods
Recombinant DNA reagentH1-shWDR5①-Ubc-GFPThis studyMethods
Recombinant DNA reagentH1-shWDR5②-Ubc-GFPThis studyMethods
Recombinant DNA reagentH1-shWDR5③-Ubc-GFPThis studyMethods
Recombinant DNA reagentLenti-dCAS-VP64-BlastAddgene# 61425
Recombinant DNA reagentLenti-sgRNA(MS2)-zeo backboneAddgene# 61427
Recombinant DNA reagentLenti MS2-P65-HSF1_HygroAddgene# 61426
Recombinant DNA reagentLenti-sgRNA(MS2)–1845 gRNA 1This studyMethods
Recombinant DNA reagentLenti-sgRNA(MS2)–1845 gRNA 2This studyMethods
Recombinant DNA reagentLenti-sgRNA(MS2)–1845 gRNA 3This studyMethods
Biological sample (Homo sapiens)Human 17 w ovary RNABIOPIKE
Antibodyanti-DAZL (Mouse
monoclonal)
Bio-RadCat# MCA2336,
RRID:AB_2292585
(IF 1:50)
Antibodyanti-LHX8 (Rabbit
polyclonal)
abcamCat# ab221882(IF 1:50; WB 1:1000)
Antibodyanti-WDR5 (Rabbit monoclonal)Cell Signaling
Technology
Cat# 13105S,
RRID:AB_2620133
(IF 1:200; WB 1:2000)
Antibodyanti-GAPDH (Mouse monoclonal)CwbioCat# CW0100M,
RRID:AB_2801390
(WB 1:5000)
Antibodyanti-FLAG (Rabbit monoclonal)SigmaCat# F3165,
RRID:AB_259529
(IP 3 μg)
Antibodyanti-H3K27Ac (Rabbit monoclonal)Cell Signaling
Technology
Cat# 8173,
RRID:AB_10949503
(ChIP 1 μg/rxn)
Antibodyanti-H3K4me3 (Rabbit monoclonal)MilliporeCat# 04-745,
RRID:AB_1163444
(ChIP 1 μg/rxn)
AntibodyGoat anti-Rabbit IgG (H + L) Cross-Adsorbed Secondary
Antibody, Alexa Fluor 488
(Goat polyclonal)
Thermo Fisher ScientificCat# A-11008,
RRID:AB_143165
(WB 1:5000)
AntibodyGoat anti-Mouse IgG (H + L)
Highly Cross-Adsorbed Secondary
Antibody, Alexa Fluor Plus
594 (Goat polyclonal)
Thermo Fisher ScientificCat# A32742,
RRID:AB_2762825
(WB 1:5000)
AntibodyHorseradish-labeled goat anti-rabbit IgG (H + L) (Goat polyclonal)ZSGB-BioCat# ZB-2301,
RRID:AB_2747412
(WB 1:5000)
AntibodyPeroxidase-Conjugated Goat anti-Mouse IgG (H + L)
(Goat polyclonal)
ZSGB-BioCat# ZB-2305,
RRID:AB_2747415
(WB 1:5000)
Chemical compound, drugPuromycinSigmaCat# P9620
Chemical compound, drugGeneticinThermo Fisher ScientificCat# 10131035
Chemical compound, drugBlasticidinThermo Fisher ScientificCat# 461120
Chemical compound, drugROCK inhibitorstemRDCat# Y-005
Peptide, recombinant proteinRecombinant human BMP-4R&DCat# 314 BP
Peptide, recombinant proteinRecombinant human BMP-8aR&DCat# 1073-BPC
Peptide, recombinant proteinRecombinant human FGF basicR&DCat# 233-FB-001MG/CF
Commercial assay or kitGateway LR Clonase II KitThermo Fisher ScientificCat# 11791100
Commercial assay or kitPrimeScript RT reagent
Kit with gDNA Eraser
TaKaRaCat# RR047A
Commercial assay or kitFluorescent In Situ
Hybridization Kit
RIBOBIOCat# C10910
Commercial assay or kitRNAmax-T7 kitRIBOBIOCat# R11073
Commercial assay or kitTnT Quick Coupled
Transcription/
Translation Systems
PromegaCat# L1170
Commercial assay or kitFirstChoice RLM-RACE kitAmbionCat# AM1700
Commercial assay or kitAMPure XP beadsBeckman CoulterA63881
Commercial assay or kitNEBNext DNA
Library Prep Kit
NEB#E7645S
Commercial assay or kitNEBNext Multiplex OligosNEB#E7335S
Commercial assay or kitVAHTS mRNA-seq v2
Library Prep Kit
VazymeNR611-02
Software, algorithmGraphPad Prism 6.0GraphPad Prismhttps://www.graphpad.com/
Software, algorithmFlowJoFlowJohttps://www.flowjo.com/
Software, algorithmDAVID Bioinformatics ResourcesHuang et al., 2009https://david.ncifcrf.gov
Software, algorithmCoding Potential CalculatorKong et al., 2007http://cpc.cbi.pku.edu.cn/
Software, algorithmCoding Potential
Assessment Tool
Wang et al., 2013http://rna-cpat.sourceforge.net/
Software, algorithmPROMOMesseguer et al., 2002http://alggen.lsi.upc.es/
Software, algorithmTopHat (v2.1.1)Trapnell et al., 2009http://ccb.jhu.edu/software/tophat/index.shtml
Software, algorithmHTSeq packageAnders et al., 2015https://htseq.readthedocs.io/en/release_0.9.1/
Software, algorithmGO (DAVID)Huang et al., 2009https://david.ncifcrf.gov/home.jsp
Software, algorithmHISAT2 (v2.1.0)Kim et al., 2015http://ccb.jhu.edu/software/hisat2/manual.shtml
Software, algorithmNONCODE (v5)Fang et al., 2017http://www.noncode.org/index.php
Software, algorithmGencode (v29)Uszczynska-Ratajczak et al., 2018https://www.gencodegenes.org/
Software, algorithmSamtools (1.3.1)Li et al., 2009http://samtools.sourceforge.net/
Software, algorithmMACS2Zhang et al., 2008https://pypi.org/project/MACS2/2.1.1.20160309/
Software, algorithmIGVThorvaldsdóttir et al., 2013http://www.igv.org/
Software, algorithmTreeviewSaldanha, 2004http://jtreeview.sourceforge.net/
Software, algorithmR (3.5.1)N/Ahttps://www.r-project.org/
Software, algorithmpheatmapN/AR software
Software, algorithmDESeq2Love et al., 2014R software
Software, algorithmggplot2Wickham, 2016R software
Software, algorithmSeurat (3.2.3)Satija et al., 2015R software
Software, algorithmMonocle (2.10.1)Qiu et al., 2017R software
OtherMatrigelCorningCat# 354277Matrigel Basement Membrane Matrix
OtherDAPIInvitrogenCat# D1306Blue-fluorescent DNA stain
OtherTrypLE ExpressInvitrogenCat# 12605010Trypsin substitution
OtherDynabeads Protein AInvitrogenCat# 10002DDynabeads Protein A for Immunoprecipitation
OtherDynabeads Protein GInvitrogenCat# 10003DDynabeads Protein G for Immunoprecipitation
OtherProLong Diamond
Antifade Mountant
InvitrogenCat# P10144Antifade Solution
OtherTRIzol ReagentInvitrogenCat# 15596026RNA extraction reagents
OtherTransStart Top Green qPCR MixTransGenCat# AQ131-02qPCR reagents
OtherLipofectamine 2000InvitrogenCat# 11668019Transfection reagents
OtherLipofectamine 3000InvitrogenCat# L3000001Transfection reagents

Additional files

Supplementary file 1

Weekly expressed lncRNAs of female and male germ cells and gonadal somatic cells.

https://cdn.elifesciences.org/articles/78421/elife-78421-supp1-v2.xlsx
Supplementary file 2

Human germ cells and gonadal somatic cells specific lncRNAs.

https://cdn.elifesciences.org/articles/78421/elife-78421-supp2-v2.xlsx
Supplementary file 3

Human germ cell each stage-specific lncRNAs.

https://cdn.elifesciences.org/articles/78421/elife-78421-supp3-v2.xlsx
Supplementary file 4

Feature protein-coding genes or lncRNAs of each cluster.

https://cdn.elifesciences.org/articles/78421/elife-78421-supp4-v2.xlsx
Supplementary file 5

Female germ cell (fGC) meiotic and male germ cell (mGC) mitotic arrest lncRNA–mRNA pairs and expression correlation.

https://cdn.elifesciences.org/articles/78421/elife-78421-supp5-v2.xlsx
Supplementary file 6

Upregulated lncRNA genes in primordial germ cell-like cells (PGCLCs).

https://cdn.elifesciences.org/articles/78421/elife-78421-supp6-v2.xlsx
Supplementary file 7

Upregulated protein-coding genes in primordial germ cell-like cells (PGCLCs).

https://cdn.elifesciences.org/articles/78421/elife-78421-supp7-v2.xlsx
Supplementary file 8

Differentially expressed protein-coding genes in LNC1845 KO primordial germ cell-like cells (PGCLCs).

https://cdn.elifesciences.org/articles/78421/elife-78421-supp8-v2.xlsx
Supplementary file 9

Differentially expressed protein-coding genes in LHX8 OE primordial germ cell-like cells (PGCLCs).

https://cdn.elifesciences.org/articles/78421/elife-78421-supp9-v2.xlsx
Supplementary file 10

Oligonucleotides.

https://cdn.elifesciences.org/articles/78421/elife-78421-supp10-v2.xlsx

Download links