library_construction_protocol: Seminiferous tubules were isolated from the testes of mice after removal of tunica  albuginea. Single cell suspension were prepared by passing the cells through the 40 μm  strainer. Cells were collected in collected in PBS and washed one time. Then pellet was  resuspended in FACS incubation buffer (HBSS supplemented with 20 mM HEPES (pH  7.2), 1.2 mM MgSO4 , 1.3 mM CaCl 2 , 6.6 mM sodium pyruvate, 0.05% lactate, glutamine,  and 1% fetal calf serum) at a density of 2 million cells/ml. Then 5 μg/ml of Hoechst 33342  was added and cells were incubated for 1hr at 32°C in a water bath. Before the sorting 2  μg/ml of propidium iodide (PI) is added to exclude the dead cells. Meiocytes were sorted  with 85 μm in an ARIA II (BD Biosciences). PI and Hoechst were excited using lasers 488  nm and 355 nm respectively. 685 LP and 710/40 BP filter were used to detect PI signal.  Hoechst blue signal was detected using 600 LP and 620/10 BP filters in front of first  detector and Hoechst red was detected with 440/40 BP filter in front of the second  detector. Cells were collected in PBS and frozen at −80°C for subsequent RNA  sequencing, qPCR and protein extraction. For ATAC sequencing 5000 cells per genotype  were collected. FlowJo (Tree Star Inc.) and FACSDiva software were used for cell sorting  and data analysis. Library preparation was prepared according to the 10x genomics guidelines. 5000 spematocytes were sorted into 30 µl PBS and transposition is carried in 10 µl of reaction buffer.  Following transposition DNA was purified either with bead based methods or with Qiagen MinElute Kit. Index PCR was performed in 25 µl of reaction buffer mixture to prepare the libraries. Libraries were purified twice with 1.3x volume of GE beads. Size distribution were analyzed with fragment analyzer. Library size distribution were analyzed by the Fragment Analyzer (company) and samples diluted in 1:3 ratio.