NR0B1 Supporting Information Biorxiv
Total Page:16
File Type:pdf, Size:1020Kb
Supporting Information: Supporting Figures - Figure S1. Comparison of Xp21.2 CNVs. Schematic representation comparing NR0B1 locus copy number variations in patients presented in this study with those previously described in the literature and associated with 46, XY sex reversal. The green bars represent the duplications and the triplication described in the patients of this study (P1 -P3) and the blue bars represent duplications previously reported by others researchers (Barbaro, Oscarson et al. 2007, Barbaro, Cicognani et al. 2008, Ledig, Hiort et al. 2010, White, Ohnesorg et al. 2011, Barbaro, Cook et al. 2012, Dong, Yi et al. 2016, Garcia-Acero, Molina et al. 2019) . All previously reported duplications include the NR0B1 (red underlined), CXorf21 and GK gene. The red bar marks the deletion reported by Smyk et al. (Smyk, Berg et al. 2007). There is a mutual 35kb region present in all Xp21.2 copy number variations associated with 46,XY GD, highlighted by the green box. The representation was generated using the UCSC Genome Browser (GRCh37/hg19)(Kent, Sugnet et al. 2002). 1 Figure S2. Sox9 activation during testicular development. Mechanism of Sox9 regulation, as suggested by Ludbrook et al. (2012). Sox9 gene with upstream enhancer Enh13 as grey bar. Grey spheres symbolize SF1 and SRY/SOX9 binding sites within the enhancer. Besides Enh13 (Gonen, Futtner et al. 2018) also an enhancer termed TESCO (Sekido and Lovell-Badge 2008) has previously been proposed as especially relevant for testis specific Sox9 expression in mice. Sox9 enhancers appear to be activated in three phases (i-iii) (Sekido and Lovell-Badge 2008). The individual enhancers may have varying phases of predominant activation and relevance for sex determination. (i) Sox9 expression is initiated by SF1, independently of SRY. (ii) Following initiation, Sox9 is upregulated synergistically by SRY and SF1. Sry is only transiently expressed in mice (Sekido, Bar et al. 2004), due to repression by SOX9 (Chaboissier, Kobayashi et al. 2004). (iii) Sox9 enhancer activation is maintained in SRY absence through autoregulation by SOX9. SOX9 substitutes SRY as they share the same DNA-Binding motif and it also synergizes with SF1. Thus, SRY is only briefly needed to increase Sox9 expression above a threshold allowing self-sustained expression (Sekido and Lovell-Badge 2009). If this threshold is not reached in time, this leads to the commitment of the supporting cell line precursors to the ovary determining pathway through accumulating ovary-determining/anti- testis factors (Gonen, Futtner et al. 2018). It is of note however, that in humans SOX9 enhancer duplications were found in several 46,XX SRY-negative males. This manifestation is attributed to enhancers termed eSR-A (enhancer sex reversal A) and eSR-B, which seem to be predominantely activated through SOX9 and SF1 (Croft, Ohnesorg et al. 2018). SOX9 is essential for Sertoli cell differentiation and activates downstream targets including Fgf9/Fgfr2 and Sox8 (Sekido and Lovell- Badge 2008, Ludbrook, Bernard et al. 2012). NR0B1 (DAX1) interferes with SF1 activation of target genes possibly through protein-protein interactions (Nachtigal, Hirokawa et al. 1998) or displacement from DNA binding sites in a dose-dependent manner (Ludbrook, Bernard et al. 2012) SRY is insufficient to activate Sox9 expression alone, thus the necessary threshold to induce Sox9 autoregulation and subsequent testis formation is not reached (Sekido and Lovell-Badge 2008). 2 Figure S3: Sex development genes in patients 1 & 2 show intermediate expression between ovary and testis tissue. The combined beeswarm and box-whisker plots depict the log2-transformed transcript per million (TPM) expression for selected sex development genes of the GTEx ovary (left) and testis (middle) data (Consortium 2013, Carithers, Ardlie et al. 2015). Expression in the two patients is depicted on the right. The boxes and thick black bar denote the first to third quartile of the data and the whiskers extend the box by 1.5 of the interquartile range. A moderated t-test shows significant differential regulation of genes between all groups. 3 Figure S4. MLPA and CMA Analysis of P1. (A) MLPA of P1 revealing the duplication of CXorf21 probes at Xp21.2 [30,595,621-30,615,321] (GRCh37/hg19); signals for NR0B1, SRY, SOX9 and NR5A1 are normal. (B) Normal probe signals for DMRT1, CYP17A1, SRD5A2 and HSD17B3 genes in MLPA. (C) Chromosome micro array (CMA) analysis using software Chromosome Analysis Suite (ChAS - Affymetrix®) revealing a ≈277kb duplication of Xp21.2 (arr[GRCh37] Xp21.2(30580693_30857187)x2) covered by 165 markers. The minimal duplicated region includes the fully duplicated GK gene and parts of TAB3 and CXorf21. 4 Figure S5. aCGH and GS comparison in P2 at NR0B1 locus (A) aCGH results for P2. Green dots indicate probes used for hybridization. Results show two separate duplications as illustrated by increased hybridization of probes and emphasized by green bars. To the left a 389kb duplication (chrX:29,924,420-30,313,761; GRCh37/hg19) was detected by increased hybridization of 17 probes. The second duplication was detected as 447kb of size (chrX: 30,401,819-30,848,988) indicated by 23 probes. The duplications are separated by four probes labelled 1 – 4. Probes 1, 3 and 4 hybridized only to a single copy. Probe 2 showed no hybridization. (B) Probes 1 to 4 have been compared to genome sequencing data. For visualisation of WGS data of P2 (region ChrX:30,292,562-30,423,695) the Integrative Genome Viewer (IGV) was used (Robinson, Thorvaldsdottir et al. 2011). Increased number of reads (grey bars) at both sides of the viewing window signify duplication detected by WGS (green bars). Deletions are evident through the lack of reads in the two regions marked with red bars. Single hybridized CGH probes 1, 3 and 4 map to areas not duplicated according to WGS data. Probe 2 maps to a deletion and correspondingly did not hybridize. 5 Figure S6. aCGH in P3 at NR0B1 locus. Green dots indicate probes used for hybridization. Results show a 1,2 Mb triplication (ChrX:29,851,537-31,069,736) emphasized by green bar. The genes MAGEB1-4, NR0B1, CXorf21 and TAB3 are fully triplicated whereas IL1RAPL1 is only partially triplicated. 6 Figure S7. Structural Variation of Patient 3 including proposed TAD structure at NR0B1 locus. Genomic distances are not to scale. The chromosome segment is drawn with the distal chromosome arm to the left and centromere to the right. WGS showed the triplication originally mapped to chrX: 29,849,782 – 31,088,713 and is arranged in a tandem manner. * marks two identical 49 bp inserts at two of the breakpoints. The sequence originates from an intron of the DMD gene, located to the centromere of the SV event. Grey arrow heads represent binding sites of the insulator CCCTC-binding factor (CCTF). The arrangement is such, that between every NR0B1 gene and closest putative enhancer there is an insulator (CTCF) binding site. However, these sites only perform insulation when bound to a partnering CTCF. In lack of sufficient number of partnering CTCFs we propose three different possible TAD conformations (dashed grey lines). Each conformation enables NR0B1 enhancer adoption in different sets of gene/enhancer pairs (arrows). 7 Figure S8. NR0B1 TAD conformation of different human cell lines according to PennState 3D GenomeBrowser data. Green bars represent duplications of patients described in this study (P1 & P2). The red bar shows the deletion reported by Smyk et al. (2007). The TAD encompassing NR0B1 (NR0B1 TAD) and 8 the adjacent centromeric and telomeric TADs identified by Hi-C sequence data and made available through the PennState 3D Genome Browser (Wang et al., 2018) have been extracted for 36 cell lines. TAD position and size is visualized as blue bars for corresponding cell lines using the UCSC Genome Browser (GRCh37/hg19) (Kent, Sugnet et al. 2002). Most cell lines show a conserved NR0B1 TAD boundary at CXorf21, which is in a region shared by the duplications of P1 and P2 as well as the deletion reported by Smyk et al. (2007). 9 Figure S9. Binding sites in mutual region between P1 - P2 duplications, P3 triplication and the deletion reported by Smyk et al. (2007). Minimal region of overlap (35kb) between duplication of P1 reported in this study (green bars) and the deletion reported by Smyk et al. (2007) (red bar). This consensus region also covered by P2 and P3 encompasses the whole CXorf21 gene. ENCODE Chip-seq Data (Consortium 2012, Davis, Hitz et al. 2018) available in the UCSC genome browser (Kent, Sugnet et al. 2002) shows a CTCF binding site in this region together with binding sites for the cohesin complex components RAD21 and SMC3. We propose this to be the centromeric NR0B1 TAD border. 10 A B S PERMATOGENES IS 0.6 PHOX2B WNT BETA CATENIN SIGNALING GMNN LMX1A EPITHE LIAL MES ENCHYMAL TRANSITION 0.4 NOTCH S IGNALING SIX3 0.5 E E TGF BETA S IGNALING MEIS2 n n r i c 0.2 r UV RESPONSE DN i TCF15 c h h m ADIPOGENESIS m E2F6 e e n n P53 PATHWAY 0 ZNF143 0 t S t S MYOGE NESIS NKX6.1 c c o o r RFX2 e PI3K AKT MTOR S IGNALING r -0.2 e IL2 STAT5 S IGNALING CREB3 TNFA SIGNALING VIA NFKB PHOX2A -0.5 INTERFERON ALPHA RES PONSE -0.4 FOSL2.JUN IL6 J AK STAT3 S IGNALING HOXB9 INTERFERON GAMMA RES PONSE TBX3 -0.6 ANGIOGENES IS DLX2 COAGULATION MAF APICAL J UNCTION TEAD4 HYPOXIA NR1I3 APOPTOSIS AHR PROTEIN S ECRETION EPAS1 OXIDATIVE PHOSPHORYLATION PRDM1 MYC TARGETS V1 SMAD1 UNFOLDED PROTE IN RESPONSE TBX2 MYC TARGETS V2 DBP UV RESPONSE UP DLX5 G G G G G G G P P G T T T T T T e e