Single-Cell DNA Methylation Analysis of Chicken Lampbrush Chromosomes

DNA methylation is an essential epigenetic regulation mechanism implicated in transcription and replication control, developmental reprogramming, retroelements silencing and other genomic processes. During mammalian development, a specific DNA methylation pattern should be established in germ cells to allow embryonic development. Less is known about germ cell DNA methylation in other species. To close this gap, we performed a single-cell methylome analysis of chicken diplotene oocytes. We comprehensively characterized methylation patterns in these cells, obtained methylation-based chicken genome segmentation and identified oocyte-specific methylated gene promoters. Our data show that despite the formation of specific transcriptionally hyperactive genome architecture in chicken diplotene oocytes, methylation patterns in these cells closely resemble genomic distribution observed in somatic tissues.


Introduction
DNA methylation in vertebrates is implicated in endogenous and exogenous gene silencing, heterochromatin formation, establishment, and maintenance of chromatin architecture, DNA repair, and other genomic processes. In contrast to other epigenetic modifications, DNA methylation can be inherited through cell divisions, allowing the transfer of epigenetic information across the mitotic cycle in somatic cells. To reset this epigenetic memory, there are specific mechanisms governing epigenetic reprogramming and the establishment of germ-line-specific methylation patterns in germ cells and early embryos.
During mammalian development, sperm progenitors undergo two rounds of epigenetic reprogramming, and at the end of gametogenesis, the mature sperm genome exhibits 80-90% CpG methylation. The genomic distribution of DNA methylation in sperm is roughly similar to that of somatic cells, although with a slightly higher DNA methylation content [1]. In oocytes, methylation gradually increases during development, starting from practically unmethylated non-growing oocytes present in the primordial follicle before folliculogenesis, and reaching the level of~40-50% in mature gametes. The establishment of methylation patterns in oocytes depends on transcription, thus methylated regions are mainly restricted to the transcribed gene bodies, leaving intergenic regions and nontranscribing genes hypomethylated. Overall, mammalian gametes show specific patterns and levels of DNA methylation, which is required to maintain the developmental program of these species [2].
There is substantially less information about DNA methylation in gametes of nonmammalian vertebrate species. Here, we focus on epigenetic profiling of germ cells in domestic chicken, Gallus gallus. Chicken is one of the most characterized non-mammalian vertebrate species, and chicken models play a pivotal role in animal research as an alternative and outbreed experimental species. In addition, domestic chicken is the globally most important source of commercially produced meat. Chicken models were extensively studied using modern genomic approaches, providing rich data about genome architecture [3], transcription dynamics [4,5], chromatin accessibility [6], and methylation patterns [7] in somatic cells. Moreover, DNA methylation in chicken sperm cells was recently profiled, showing that sperm DNA is hypomethylated, which correlates with the absence of the DNMT3L cofactor in the chicken genome [7]. No genomic data are available to study DNA methylation in growing chicken oocytes. During the diplotene stage of prophase I of meiosis, the chicken oocyte genome acquires specific organization, characterized by elevated transcription and formation of lampbrush chromosome structures. Although cytological evidence suggests that DNA methylation in growing chicken oocytes correlates with chromatin organization and transcription, there are no molecular data available [8,9].
Here, we profile DNA methylation in two stages of chicken oocyte diplotene. We show that oocyte methylation at both stages shows patterns previously observed in somatic cells, including hypomethylation of CpG islands and hypermethylation of transposons and other repetitive DNA elements. In contrast to sperm, average levels of DNA methylation in oocytes are on par with somatic cells. We developed a new HMM-based segmentation tool that identifies hypomethylated regions in the chicken genome and showed that there is a limited number of oocyte-specific hypomethylated gene promoters. Overall, our data suggest that DNA methylation does not significantly contribute to the developmental reprogramming of the oocyte and lampbrush chromosome formation in birds.

Chicken Oocytes Isolation and Bisulfite Sequencing
To profile oocyte DNA methylation in adult chickens and study its dynamics, we chose SWF (small white follicles) and LWF (large white follicles) stages, which can be easily distinguished using microscopy analyses. Both SWF and LWF represent diplotene stages, roughly similar to 10-15 days post-partum stages of mammalian oocyte development [2]. We defined SWF and LWF stages based on the follicle size (1-2 mm follicle for the SWF stage, 4-7 mm follicle for the LWF stage). To ensure that follicles with different sizes contain oocytes at different developmental stages, we isolated oocyte nuclei and accessed their size and chromosome morphology ( Figure 1). During the SWF stage, a typical avian diplotene nucleus can be observed-with lampbrush bivalents evenly distributed across the nucleus ( Figure 1A), which corresponds to previously published data [10]. Nuclei at later stages slightly increase in size, but the main difference between SWF and LWF stages is chromosome morphology. At the LWF stage, chromosomes tend to gather together, and their morphology is altered to a far more condensed state ( Figure 1C) [11]. Thus, we confirmed that the stratification of oocytes by follicle size allows for unambiguous discrimination between lampbrush and post-lampbrush stages.
In order to avoid contamination by follicular cells surrounding the oocyte, we isolated the oocyte nuclei as described previously [12]. We performed bisulfite sequencing of 16 oocyte nuclei from SWF and LWF stages (7 SWF and 9 LWF oocytes collected from 3 animals) generating from 3 to 7 million read pairs for each oocyte which results in coverage of 2-4 million CpG dinucleotides per oocyte (Table S1). Merged data for all oocytes allowed us to access methylation in 69.73% of CpG dinucleotides. Our optimized mapping setup for bismark allowed us to reach average mapping efficiency of 69.7% which is on par with existing single-cell bisulfite sequencing data on mammalian oocytes [13].
Given the fact that during the diplotene stage there are four chromatids of each chromosome and each converted DNA strand is sequenced independently, we can expect no more than eight unique reads covering each CpG dinucleotide. Thus, cases, where dinucleotides are covered more than eight times, can be considered artifacts or contamination from somatic cells. Our data show that there is a very small amount (~0.1-0.2%) of artifacts across all libraries ( Figure 2A). Thus, our approach for bisulfite sequencing of individual chicken oocyte nuclei performs on par with existing single-cell approaches by showing high mapping efficiency and reproducibility. In order to avoid contamination by follicular cells surrounding the oocyte, we isolated the oocyte nuclei as described previously [12]. We performed bisulfite sequencing of 16 oocyte nuclei from SWF and LWF stages (7 SWF and 9 LWF oocytes collected from 3 animals) generating from 3 to 7 million read pairs for each oocyte which results in coverage of 2-4 million CpG dinucleotides per oocyte (Table S1). Merged data for all oocytes allowed us to access methylation in 69.73% of CpG dinucleotides. Our optimized mapping setup for bismark allowed us to reach average mapping efficiency of 69.7% which is on par with existing single-cell bisulfite sequencing data on mammalian oocytes [13].
Given the fact that during the diplotene stage there are four chromatids of each chromosome and each converted DNA strand is sequenced independently, we can expect no more than eight unique reads covering each CpG dinucleotide. Thus, cases, where dinucleotides are covered more than eight times, can be considered artifacts or contamination from somatic cells. Our data show that there is a very small amount (~0.1-0.2%) of artifacts across all libraries ( Figure 2A). Thus, our approach for bisulfite sequencing of individual chicken oocyte nuclei performs on par with existing single-cell approaches by showing high mapping efficiency and reproducibility.

Methylation Profile of Chicken Oocytes
To compare levels and patterns of DNA methylation in oocytes with other chicken cell types we uniformly processed publicly available data for fibroblasts [14], spleen, jejunum, ileum [7], muscles [15], and sperm cells [16]. We found that the average level of CpG methylation in oocytes (53.2%) was comparable to those observed in somatic cells (52.45-

Methylation Profile of Chicken Oocytes
To compare levels and patterns of DNA methylation in oocytes with other chicken cell types we uniformly processed publicly available data for fibroblasts [14], spleen, jejunum, ileum [7], muscles [15], and sperm cells [16]. We found that the average level of CpG methylation in oocytes (53.2%) was comparable to those observed in somatic cells (52.45-62.8%). We also confirmed that the sperm genome is hypomethylated (40.5%), as was previously reported [7].
We next compared methylation levels on micro-and macro-chromosomes, and methylation patterns of different genomic features. In general, there are no major differences in DNA methylation patterns between chicken lampbrush chromosomes and somatic cells ( Figure 2B). Promoters and CpG islands showed bimodal distribution: whereas the majority of these loci were hypomethylated, a small fraction displayed high methylation levels. Exons also showed a bimodal distribution of methylation, where the majority of the loci were highly methylated and few percentages of exons showed no methylation. Similarly, 3 -untranslated regions were highly methylated, whereas promoter-proximal 5 -untranslated regions were almost completely demethylated. Intergenic sequences, repetitive elements and introns showed increased methylation levels ( Figure 2C).
We did not detect any substantial difference between micro-and macrochromosomes, except slightly higher average methylation levels on macrochromosomes, which could be explained by a higher portion of introns, repetitive elements and intergenic features. Sperm cells show decreased methylation levels for all genomic features analyzed, confirming again the previously described hypomethylation of the sperm genome ( Figure 2D). We next measured the methylation of all 5-kb genomic bins and compared the methylation of individual loci between chicken cell types. This analysis showed a moderate correlation between methylation profiles in oocytes and somatic cells (Spearman's R~0.6-0.7, Figure 3A). We noticed a difference between hypomethylated (<40%) and hypermethylated (>60%) loci: the latter always show more dissimilarities between cell types (except muscle cells that differ in methylation pattern from all somatic cells; Figure S1). A pairwise comparison of cell types showed that methylation level within hypermethylated loci is less consistent between cell types than in hypomethylated regions ( Figure S1). Nevertheless, methylation of the CpG islands and gene promoters, which often fall within hypomethylated regions, was very similar. Moreover, this methylation pattern does not reflect known features of oocyte-specific transcription. For example, the histone genes cluster was demethylated in oocytes similarly to somatic cell types ( Figure S2), despite the absence of histone locus bodies and inactivation of the histone gene cluster [10]. HoxD gene cluster contains an extremely long (~80-kb, chr7:16,335,000-16,407,000) demethylated region, which was observed in all analyzed cell types ( Figure S3).
Since lampbrush chromosomes display unique chromatin organization compared to somatic cell interphase nuclei, we examined whether the methylation of CTCF binding sites differs between diplotene oocytes and fibroblasts. CTCF is a major vertebrate insulator that plays an important role in genome folding by mediating the formation of chromatin loops [17,18]. Reportedly, domain borders in lampbrush chromosomes do not correspond with somatic TAD boundaries [19]. However, we did not find any significant differences in levels of DNA methylation between oocytes and fibroblasts in CTCF sites suggesting that there are different mechanisms involved in chromatin domain insulation in lampbrush chromosomes ( Figure S4). It is known that some repetitive sequences, such as CNM, PO41, and Z-microsatellite are transcriptionally active during the lampbrush stage [20], whereas others, such as chromosome-specific centromeric repeats Cen1, Cen2, Cen3, Cen4, Cen7, Cen8, and Cen11, are transcriptionally silent [21]. Due to their repetitive nature, methylation analysis of these sequences cannot be performed using standard genomic alignment. We manually designed sequences representing methylated and demethylated states of these repetitive elements after bisulfite conversion. Next, we aligned reads using these references and counted a portion of alignments corresponding to the methylated and demethylated states. The results of this analysis ( Figure S5) showed that all repeats except PO41 display relatively high methylation levels (50-60%). However, the distinctive results of PO41 methylation (10-20% in most of the cell types) should be interpreted with caution, because this repetitive element was covered by a minimal amount of reads (~10 reads in each cell type). We did not observe any correlation between the transcription status of the repeat and its methylation: actively transcribed CNM and Z-microsatellite repeats show relatively the same methylation level as non-transcribed Cen-repeats. Interestingly, methylation levels of all repeat sequences in oocytes and sperm cells were lower than in somatic cells.
Finally, we compared methylation of SWF and LWF stages of oocyte development. We found no difference in methylation levels between the stages (p-value = 0.11); moreover, the difference between individual animals was substantially higher than the variation between different oocyte stages ( Figure S6). We did not observe any difference in the methylation patterns for the analyzed genomic features, such as promoters, enhancers, and CpG islands, specific for SWF or LWF. Thus, we concluded that, in contrast to mammals, there are no changes in methylation during SWF to LWF oocyte transition. Overall, It is known that some repetitive sequences, such as CNM, PO41, and Z-microsatellite are transcriptionally active during the lampbrush stage [20], whereas others, such as chromosome-specific centromeric repeats Cen1, Cen2, Cen3, Cen4, Cen7, Cen8, and Cen11, are transcriptionally silent [21]. Due to their repetitive nature, methylation analysis of these sequences cannot be performed using standard genomic alignment. We manually designed sequences representing methylated and demethylated states of these repetitive elements after bisulfite conversion. Next, we aligned reads using these references and counted a portion of alignments corresponding to the methylated and demethylated states. The results of this analysis ( Figure S5) showed that all repeats except PO41 display relatively high methylation levels (50-60%). However, the distinctive results of PO41 methylation (10-20% in most of the cell types) should be interpreted with caution, because this repetitive element was covered by a minimal amount of reads (~10 reads in each cell type). We did not observe any correlation between the transcription status of the repeat and its methylation: actively transcribed CNM and Z-microsatellite repeats show relatively the same methylation level as non-transcribed Cen-repeats. Interestingly, methylation levels of all repeat sequences in oocytes and sperm cells were lower than in somatic cells.
Finally, we compared methylation of SWF and LWF stages of oocyte development. We found no difference in methylation levels between the stages (p-value = 0.11); moreover, the difference between individual animals was substantially higher than the variation between different oocyte stages ( Figure S6). We did not observe any difference in the methylation patterns for the analyzed genomic features, such as promoters, enhancers, and CpG islands, specific for SWF or LWF. Thus, we concluded that, in contrast to mammals, there are no changes in methylation during SWF to LWF oocyte transition. Overall, our data showed that chicken oocytes at both SWF and LWF stages display a methylation pattern very similar to those previously described in chicken and mammalian somatic cells.

Methylation-Informed Segmentation of the Chicken Genome Identifies Oocyte-Specific Hypomethylated Regions
In contrast to the limited difference in methylation patterns between cell types, we observed well-pronounced differences in methylation levels between different loci ( Figure 3B). As was mentioned above, low methylation was often attributed to promoters and CpG islands; however, the overlap between demethylated regions and these genomic features was not complete (the quantification of the overlap size is provided below).
To identify differentially methylated genomic regions, we developed a segmentation method based on the Hidden Markov model (HMM). We modified standard multinomial HMM to account for the unequal distance between neighboring CpG sites and unequal coverage of each CpG site (see a more detailed description of modifications in Methods). Iterating over model parameters showed that the observed methylation levels should be explained by two or three different methylation regimes, therefore we used three methylation regimes for all subsequent studies. These regimes, according to the average methylation levels of the corresponding loci, were assigned to hypo-, intermediate-and hypermethylated regimes ( Figure 3C).
Genomic segmentation obtained using HMM corresponds well to the visual observations ( Figure 3B). All three methylation regimes can be well-discriminated based on the average methylation level of the locus; however, the difference between hypomethylated and other methylation regimes was the most pronounced. The vast majority of hypomethylated loci show methylation level that does not exceed 20% (Table S2). We identified more than 40,000 hypomethylated loci in oocytes, with a mean size of~650 b.p. (the exception was chromosome Z, where the mean hypomethylated loci length exceeds 1 kb). There was slight, although significant enrichment of gene promoters across hypomethylated loci (656 vs. 455 ± 21 expected at random), and more than two-times enrichment of CpG islands (634 vs. 286 ± 19 expected at random), suggesting that hypomethylated loci represent regulatory genomic elements. However, many hypomethylated loci do not match any known promoter or CpG island. Consistent with our previous results, loci identified as hypomethylated in oocytes often show low average methylation levels in other cell types. Using a genome segmentation obtained on oocyte data and a threshold of 20% CpG methylation as a discriminative feature of the hypomethylated loci, we compared the distribution of hypomethylated loci in oocytes and other cell types. We defined oocyte-specific loci as regions that are either hypomethylated in oocytes but not in other cell types, or hypomethylated in all cell types but not oocytes. Note that in this analysis we excluded sperm cells from the comparison because the methylation of these cells strongly differs from other cell types ( Figure 2B). Using this definition, we found~197 loci specific for oocytes (Table S3).
Only 656 of the identified loci overlapped promoters, and after visual inspection, we filtered out many of them based on the signal quality (Table S4). This results in 194 genes that show oocyte-specific patterns of methylation in promoters. Interestingly, several genes implicated in early development showed distinct methylation profiles in oocytes. For example, the promoter area of the NKX2-6 gene (homeobox-containing early development transcription regulator) contained an oocyte-specific hypomethylated region ( Figure 3D and Figure S7), whereas for KLF4 (pluripotency-associated transcription factor) we observed oocyte-specific gain of promoter methylation. Moreover, we identified an oocyte-specific decrease in promoter methylation of the TET2 (Tet Methylcytosine Dioxygenase 2) gene. Increased expression of this demethylation enzyme may be responsible for oocyte-specific hypomethylation. However, gene ontology analysis on the complete set of oocyte-specific hypomethylated promoters did not show any enrichment of developmental terms. This, again, suggests that methylation patterns observed in oocytes are largely consistent between cell types and there is only a limited number of loci showing oocyte-specific methylation.

Discussion
Here, we present a comprehensive characterization of chicken lampbrush-stage oocyte methylome. Computational analysis and concordance between our results and previous cytological analyses suggest a high quality of the obtained data. For example, a high level of methylation of transposons and other repetitive DNA correlates with previous cytological data where the distribution of 5-methylcytosine along the axes of chicken lampbrush chromosomes was characterized by immunostaining [9]. Enrichment of 5mC was demonstrated in compact chromomeres associated with the clusters of certain tandem repeats, leading to the general repression of these genomic elements during oogenesis [8].
To our surprise, we did not find a substantial genome-wide difference between oocyte and somatic cells methylome, despite the specific genome architecture of diplotene oocytes. All differences that we observed were limited to a small number of loci, and the number of differentially methylated loci was on par with a variability of somatic cell types. Similarly, transcriptional silencing that occurs during SWF to LWF transition does not change the methylation pattern that appeared to be very similar at SWF and LWF stages. Thus, we concluded that epigenetic factors other than DNA methylation are responsible for lampbrush chromosome formation. Among these factors-hyperacetylation of histone H4 was observed not only on the transcription loops but also at places where they are anchored into chromomeres [9,22]. Increased chromatin accessibility followed by initiation of hypertranscription may lead to the appearance of visible extended lateral loops. In addition, we note that these observations are completely opposite to what can be seen for mammals, where during the oocyte growth phase, DNA methylation increases and imprinted loci acquire a specific methylation profile.
It was suggested that CpG methylation along lampbrush chromosomes can carry certain epigenetic information and regulate gene expression after fertilization [8]. However, it was not known whether the DNA methylation profile is preserved after bivalent condensation. Although limited in number, we identified several gene promoters that showed oocyte-specific methylation profiles. Presumably, demethylation of transcription start sites can act as labeling for activation of early development genes during zygotic activation. Further epigenetic and transcriptomic studies of SWF, LWF and later oocyte stages can be performed to test this hypothesis.
Finally, we highlight the new HMM-based approach that we developed to segment the genome onto differentially methylated states. Although HMM is a well-established genome segmentation technique [23], here we implemented some specific modifications that are important for methylation data analysis. This allowed us to produce methylation-based segmentation of the chicken genome. Although the obtained segmentation was based on the oocytes' data, due to the relatively small variability of methylation levels between cell types, we assume that this segmentation can be used as a reference in future research focused on chicken genomics.

Nuclei Isolation
Nuclei were extracted from oocytes of domestic chicken with 1-2 mm (SWF stage) and 4-7 mm (LWF stage) diameter by using a modified protocol [11]. Oocytes were isolated from the chicken ovary and transferred to a drop of "5:1+" buffer (83.0 мМ KCl, 17.0 мМ NaCl, 6.5 мМ Na 2 HPO 4 , 3.5 мМ KH 2 PO 4 , 1.0 мМ MgCl 2 , 1.0 мМ DTT), tungsten needle was used to make holes in oocyte envelope, then exposed nuclei were carefully transferred to a clean drop with a micropipette to wash away the remaining yolk. All animal experiments were approved by the bioethics committee of the Institute of Cytology and Genetics SB RAS (Protocol №66, 9 October 2020). International guidelines were followed during experimental procedures ("Guide for the Care and Use of Laboratory Animals" [24]).
Confocal laser scanning microscopy was performed by using LSM 780 (Zeiss, Jena, Germany) inverted confocal system supplied with a 488 mm argon laser and Plan-Apochromat 20x/0.8 M27 objective. A Series of optical slices was captured with 1024 × 1024 format followed by 3D reconstruction and image maximum projection views in ZEN Black software (Zeiss, Jena, Germany). Additional processing of digital images was performed by using ImageJ 1.53 (National Institutes of Health, Bethesda, MD, USA) and Photoshop CC 2021 (Adobe, San Jose, CA, USA) software.

Library Preparation and Sequencing
For the preparation of single-cell bisulfite sequencing libraries, we used a modified version of a default protocol for the Pico Methyl-Seq Library Prep Kit (Zymo Research, Irvine, CA, USA). Each isolated nuclei in a volume of 3 µL were transferred from a drop of "5:1" buffer to individual tubes with 16 µL of nuclease-free water and 0.8 units of Proteinase K (New England Biolabs, Ipswich, MA, USA) and incubated for 2 h at 65 • C. After incubation, 130 µL of Lightning Conversion Reagent was added. Further steps were performed according to standard protocol with adjustments: incubation time with L-Desulphonation Buffer was reduced to 17 min, PrepAmp primers final concentration was reduced to 2 µM and the amount of amplification cycles with LibraryAmp primers was set to 10.
Libraries were quantified using a Qubit fluorometer (Invitrogen, Waltham, MA, USA). Sequencing was performed on the Illumina platform in 150 nt paired-end modes in the ICG SB RAS.

HMM
We modified the python HMM implementation proposed by [25]. We modeled the HMM state emission as a Binomial distribution function Bin(nt,βk), where nt-t's CpG site read-depth and βk-average methylation level of HMM state k. In addition, the Binomial distribution function has been raised to the power of (nt + 2) during HMM training so that cytosines with greater read depth contribute more to HMM learning. These modifications were based on [26]. Additionally, we modified transition probability to model the decrease in a correlation of methylation levels with increasing genomic distance [27]): where S t -state of CpG site t, S t−1 -state of CpG site t − 1, k-state number, d t−1 -distance from CpG site t − 1 and t, d 0 -distance parameter.
In contrast with standard multinomial HMM that works with directive data series (such as time-dependent events), the methylation state of t's CpG site should be connected with both t − 1 and t + 1 CpG site states. To take this into account, we add the reversedirected data in the training set and decoding set. After the whole data were decoded, we intersect the forward and reverse decodings and named all CpG sites with a difference in forward and reverse states as "un" (unknown) and excluded them from further analysis.
All data sets were prepared by the set_maker.py script. For parameter optimization, the train and test subsamples of the oocyte data were used. To define train and test subsets, we selected 54 non-overlapping loci, 300 kb length each, with randomly chosen start positions so that all three types of chicken chromosomes (macro-chromosomes, interchromosomes and micro-chromosomes) had the same representation in train and test sets. The AIC and BIC metrics from HMM implementation [25,28] were used to find the optimal parameters for MethylationHMM (n_states = 3, distance_param = 1250). These optimal parameters were used to make further decoding of oocyte methylation states.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/ijms232012601/s1. Author Contributions: V.F., M.G. and A.K. conceived the study. M.G. and A.N. isolated oocytes, A.K. supervised oocytes isolation. A.N. and M.G. prepared and sequenced Hi-C libraries. T.L. performed NGS data analysis with help from A.N., T.L. developed and implemented HMM tool. V.F. supervised data analysis. All authors contributed to manuscript preparation. All authors have read and agreed to the published version of the manuscript.
Funding: This work, including cell isolation and library preparation, was supported by Russian Science Foundation grant #20-64-46021. NGS sequencing was supported by the strategic academic leadership program "Priority 2030" at Novosibirsk State University. Computation analysis was performed using equipment of the Novosibirsk State University, supported by the Ministry of Education and Science of the Russian Federation, grant #2019-0546 (FSUS-2020-0040). The microscopy was conducted at the Joint Access Center for Microscopy of Biological Objects with the Siberian Branch of the Russian Academy of Sciences, supported by project No. 121031800061-7.
Institutional Review Board Statement: The animal study protocol was approved by the bioethics committee of the Institute of Cytology and Genetics SB RAS (Protocol №66, 9 October 2020).

Informed Consent Statement: Not applicable.
Data Availability Statement: The data presented in this study are available in Methods section.

Conflicts of Interest:
The authors declare no conflict of interest.