Expression levels and DNA methylation profiles of the growth gene SHOX in cartilage tissues and chondrocytes.

All attempts to identify male-specific growth genes in humans have failed. This study aimed to clarify why men are taller than women. Microarray-based transcriptome analysis of the cartilage tissues of four adults and chondrocytes of 12 children showed that the median expression levels of SHOX, a growth gene in the pseudoautosomal region (PAR), were higher in male samples than in female samples. Male-dominant SHOX expression was confirmed by quantitative RT-PCR for 36 cartilage samples. Reduced representation bisulfite sequencing of four cartilage samples revealed sex-biased DNA methylation in the SHOX-flanking regions, and pyrosequencing of 22 cartilage samples confirmed male-dominant DNA methylation at the CpG sites in the SHOX upstream region and exon 6a. DNA methylation indexes of these regions were positively correlated with SHOX expression levels. These results, together with prior findings that PAR genes often exhibit male-dominant expression, imply that the relatively low SHOX expression in female cartilage tissues reflects the partial spread of X chromosome inactivation into PAR. Altogether, this study provides the first indication that sex differences in height are ascribed, at least in part, to the sex-dependent epigenetic regulation of SHOX. Our findings deserve further validation.

Systematic molecular analyses for 115 karyotypically normal men with isolated non-obstructive azoospermia.

STUDY QUESTION: Do copy-number variations (CNVs) in the azoospermia factor (AZF) regions and monogenic mutations play a major role in the development of isolated (non-syndromic) non-obstructive azoospermia (NOA) in Japanese men with a normal 46, XY karyotype? SUMMARY ANSWER: Deleterious CNVs in the AZF regions and damaging sequence variants in eight genes likely constitute at least 8% and approximately 8% of the genetic causes, respectively, while variants in other genes play only a minor role. WHAT IS KNOWN ALREADY: Sex chromosomal abnormalities, AZF-linked microdeletions, and monogenic mutations have been implicated in isolated NOA. More than 160 genes have been reported as causative/susceptibility/candidate genes for NOA. STUDY DESIGN, SIZE, DURATION: Systematic molecular analyses were conducted for 115 patients with isolated NOA and a normal 46, XY karyotype, who visited our hospital between 2017 and 2021. PARTICIPANTS/MATERIALS, SETTING, METHODS: We studied 115 unrelated Japanese patients. AZF-linked CNVs were examined using sequence-tagged PCR and multiplex ligation-dependent probe amplification, and nucleotide variants were screened using whole exome sequencing (WES). An optimized sequence kernel association test (SKAT-O), a gene-based association study using WES data, was performed to identify novel disease-associated genes in the genome. The results were compared to those of previous studies and our in-house control data. MAIN RESULTS AND THE ROLE OF CHANCE: Thirteen types of AZF-linked CNVs, including the hitherto unreported gr/gr triplication and partial AZFb deletion, were identified in 63 (54.8%) cases. When the gr/gr deletion, a common polymorphism in Japan, was excluded from data analyses, the total frequency of CNVs was 23/75 (30.7%). This frequency is higher than that of the reference data in Japan and China (11.1% and 14.7%, respectively). Known NOA-causative AZF-linked CNVs were found in nine (7.8%) cases. Rare damaging variants in known causative genes (DMRT1, PLK4, SYCP2, TEX11, and USP26) and hemizygous/multiple-heterozygous damaging variants in known spermatogenesis-associated genes (TAF7L, DNAH2, and DNAH17) were identified in nine cases (7.8% in total). Some patients carried rare damaging variants in multiple genes. SKAT-O detected no genes whose rare damaging variants were significantly accumulated in the patient group. LIMITATIONS, REASONS FOR CAUTION: The number of participants was relatively small, and the clinical information of each patient was fragmentary. Moreover, the pathogenicity of identified variants was assessed only by in silico analyses. WIDER IMPLICATIONS OF THE FINDINGS: This study showed that various AZF-linked CNVs are present in more than half of Japanese NOA patients. These results broadened the structural variations of AZF-linked CNVs, which should be considered for the molecular diagnosis of spermatogenic failure. Furthermore, the results of this study highlight the etiological heterogeneity and possible oligogenicity of isolated NOA. STUDY FUNDING/COMPETING INTEREST(S): This study was supported by Grants from the Japan Society for the Promotion of Science (21K19283 and 21H0246), the Japan Agency for Medical Research and Development (22ek0109464h0003), the National Center for Child Health and Development, the Canon Foundation, the Japan Endocrine Society, and the Takeda Science Foundation. The results of this study were based on samples and patient data obtained from the International Center for Reproductive Medicine, Dokkyo Medical University Saitama Medical Center, Koshigaya, Japan. The authors have no conflicts of interest to disclose. TRIAL REGISTRATION NUMBER: N/A.

Nuclear Receptor Gene Variants Underlying Disorders/Differences of Sex Development through Abnormal Testicular Development.

Gonadal development is the first step in human reproduction. Aberrant gonadal development during the fetal period is a major cause of disorders/differences of sex development (DSD). To date, pathogenic variants of three nuclear receptor genes (NR5A1, NR0B1, and NR2F2) have been reported to cause DSD via atypical testicular development. In this review article, we describe the clinical significance of the NR5A1 variants as the cause of DSD and introduce novel findings from recent studies. NR5A1 variants are associated with 46,XY DSD and 46,XX testicular/ovotesticular DSD. Notably, both 46,XX DSD and 46,XY DSD caused by the NR5A1 variants show remarkable phenotypic variability, to which digenic/oligogenic inheritances potentially contribute. Additionally, we discuss the roles of NR0B1 and NR2F2 in the etiology of DSD. NR0B1 acts as an anti-testicular gene. Duplications containing NR0B1 result in 46,XY DSD, whereas deletions encompassing NR0B1 can underlie 46,XX testicular/ovotesticular DSD. NR2F2 has recently been reported as a causative gene for 46,XX testicular/ovotesticular DSD and possibly for 46,XY DSD, although the role of NR2F2 in gonadal development is unclear. The knowledge about these three nuclear receptors provides novel insights into the molecular networks involved in the gonadal development in human fetuses.

Chromosomal microdeletion leading to pituitary gigantism through hormone-gene overexpression.

Pituitary gigantism is a rare endocrinopathy characterized by tall stature due to growth hormone (GH) hypersecretion. This condition is generally linked to a genetic predisposition to tumors that produce GH or GH-releasing hormone (GHRH). Here, we report a Japanese woman who exhibited prominent body growth from infancy to reach an adult height of 197.4 cm (+7.4 standard deviation). Her blood GH levels were markedly elevated. She carried no pathogenic variants in known growth-controlling genes but had a hitherto unreported 752 kb heterozygous deletion at 20q11.23. The microdeletion was located 8.9 kb upstream of GHRH and encompassed exons 2-9 of a ubiquitously expressed gene TTI1 together with 12 other genes, pseudogenes and non-coding RNAs. Transcript analyses of the patient's leukocytes showed that the microdeletion produced chimeric mRNAs consisting of exon 1 of TTI1 and all coding exons of GHRH. In silico analysis detected promoter-associated genomic features around TTI1 exon 1. Genome-edited mice carrying the same microdeletion recapitulated accelerated body growth from a few weeks after birth. The mutant mice developed pituitary hyperplasia and exhibited ectopic Ghrh expression in all tissues examined. Thus, the extreme phenotype of pituitary gigantism in the patient likely reflects GHRH overexpression driven by an acquired promoter. The results of this study indicate that germline submicroscopic deletions have the potential to cause conspicuous developmental abnormalities due to gene overexpression. Furthermore, this study provides evidence that constitutive expression of a hormone-encoding gene can result in congenital disease.

Maternal uniparental disomy of chromosome 7 underlying argininosuccinic aciduria and Silver-Russell syndrome.

We describe a patient presenting with argininosuccinic aciduria and Silver-Russell syndrome (SRS). SRS was caused by maternal uniparental disomy of chromosome 7 (UPD(7)mat). UPD(7)mat also unmasked a maternally inherited splicing variant in ASL on chromosome 7, leading to the onset of argininosuccinic aciduria. The phenotype of the present case was more severe than that of a previous case, demonstrating a phenotypic variation in the combination of argininosuccinic aciduria and SRS.

Quantification of Maternal Microchimeric Cells in the Liver of Children With Biliary Atresia.

ABSTRACT: Biliary atresia (BA) is a rare disorder of unknown etiology. There is a debate as to whether maternal microchimerism plays a significant role in the development of BA or in graft tolerance after liver transplantation. Here, we performed quantitative-PCR-based assays for liver tissues of children with BA and other diseases. Maternal cells were detected in 4/13 and 1/3 of the BA and control groups, respectively. The estimated number of maternal cells ranged between 0 and 34.7 per 106 total cells. The frequency and severity of maternal microchimerism were similar between the BA and control groups, and between patients with and without acute rejection of maternal grafts. These results highlight the high frequency of maternal microchimerism in the liver. This study provides no evidence for roles of microchimerism in the etiology of BA or in graft tolerance. Thus, the biological consequences of maternal microchimerism need to be clarified in future studies.

Optical Genome Mapping for a Patient with a Congenital Disorder and Chromosomal Translocation.

We performed optical genome mapping (OGM), a newly developed cytogenetic technique, for a patient with a disorder of sex development (DSD) and a 46,XX,t(9;11)(p22;p13) karyotype. The results of OGM were validated using other methods. OGM detected a 9;11 reciprocal translocation and successfully mapped its breakpoints to small regions of 0.9-12.3 kb. OGM identified 46 additional small structural variants, only three of which were detected by array-based comparative genomic hybridization. OGM suggested the presence of complex rearrangements on chromosome 10; however, these variants appeared to be artifacts. The 9;11 translocation was unlikely to be associated with DSD, while the pathogenicity of the other structural variants remained unknown. These results indicate that OGM is a powerful tool for detecting and characterizing chromosomal structural variations, although the current methods of OGM data analyses need to be improved.

Practical use of a multicenter clinical research support system connected to electronic medical records.

BACKGROUND: Electronic medical records (EMRs) are widely used, but in many cases, they are used within a network physically separated from the Internet. Multicenter clinical studies use Internet-connected electronic data capture (EDC) systems to collect data, where data entered into the EMR are manually transcribed into the EDC system. In addition, medical images for clinical research are also collected manually. Variations in EMRs and differing data structures among vendors hamper the use of data for clinical research. METHODS: We solved this problem by developing a network infrastructure for clinical research between Osaka University Hospital and affiliated hospitals in the Osaka area and introducing a clinical data collection system (CDCS). In each hospital's EMR network, we implemented a CRF reporter that accumulated data for clinical research using a template and then sent the data to a management server in the Osaka University Hospital Data Center. To organize the patient profile data and clinical laboratory data stored in each EMR for use in clinical research, the data are retrieved from the template by an interface module developed by each vendor, according to our common data output interface specification. The data entered into the CRF reporter template for clinical research are also recorded in the EMR progress notes and sent to the data management server. This network infrastructure can also be used as a medical image collection system that automatically collects images for research from PACS at each hospital. These systems are managed under common subject numbers issued by the CDCS. RESULTS: A network infrastructure was established among 19 hospitals, and a CRF reporter was incorporated into the EMR. A medical image transfer system was introduced in 13 hospitals. Since 2013, 28 clinical studies have been conducted using this system, and data for 9,987 cases have been collected as of December 31, 2020. CONCLUSION: Incorporating a CRF reporter with medical image transfer system into the EMR has proven useful for collecting research data.

Methylation status of genes escaping from X-chromosome inactivation in patients with X-chromosome rearrangements.

BACKGROUND: X-chromosome inactivation (XCI) is a mechanism in which one of two X chromosomes in females is randomly inactivated in order to compensate for imbalance of gene dosage between sexes. However, about 15% of genes on the inactivated X chromosome (Xi) escape from XCI. The methylation level of the promoter region of the escape gene is lower than that of the inactivated genes. Dxz4 and/or Firre have critical roles for forming the three-dimensional (3D) structure of Xi. In mice, disrupting the 3D structure of Xi by deleting both Dxz4 and Firre genes led to changing of the escape genes list. To estimate the impact for escape genes by X-chromosome rearrangements, including DXZ4 and FIRRE, we examined the methylation status of escape gene promoters in patients with various X-chromosome rearrangements. RESULTS: To detect the breakpoints, we first performed array-based comparative genomic hybridization and whole-genome sequencing in four patients with X-chromosome rearrangements. Subsequently, we conducted array-based methylation analysis and reduced representation bisulfite sequencing in the four patients with X-chromosome rearrangements and controls. Of genes reported as escape genes by gene expression analysis using human hybrid cells in a previous study, 32 genes showed hypomethylation of the promoter region in both male controls and female controls. Three patients with X-chromosome rearrangements had no escape genes with abnormal methylation of the promoter region. One of four patients with the most complicated rearrangements exhibited abnormal methylation in three escape genes. Furthermore, in the patient with the deletion of the FIRRE gene and the duplication of DXZ4, most escape genes remained hypomethylated. CONCLUSION: X-chromosome rearrangements are unlikely to affect the methylation status of the promoter regions of escape genes, except for a specific case with highly complex rearrangements, including the deletion of the FIRRE gene and the duplication of DXZ4.

[Prescribing Patterns for Anti-Parkinsonian Drugs in Japan: Prescription-Based Database Study].

The current therapeutic approach for Parkinson's disease (PD) is mainly dopamine replacement with levodopa and other anti-parkinsonian drugs. As PD progresses, the number of these drugs used steadily increases. Using prescription-based database for 10 or more years up to October 2019, we investigated actual prescribing patterns for anti-parkinsonian drugs in Japan. The main analyses included data from patients continuously prescribed levodopa for 1 or more years (n=16,270), and of these, those continuously prescribed adjuvants to levodopa for 1 or more years (n=3,675). The results showed that the number of anti-parkinsonian drugs, their daily dose frequencies, and the number of tablets increased over time. These trends were observed not only for levodopa but also for adjuvants to levodopa; the number of adjuvants, their daily dose frequencies and number of tablets also increased. As the daily number of tablets increased, the proportion of dopamine agonists increased. Moreover, as the daily dosage of levodopa increased, the daily number of tablets increased for both overall anti-parkinsonian drugs and adjuvants to levodopa. This study revealed the process of polypharmacy in PD treatment objectively. Our results are valuable for maintaining and improving therapeutic adherence in PD. (Received 25 August, 2020; Accepted 23 October, 2020; Published 1 March, 2021).

Lupus anticoagulant hypoprothrombinemia syndrome associated with bilateral adrenal haemorrhage in a child: early diagnosis and intervention.

BACKGROUND: Lupus anticoagulant-hypoprothrombinemia syndrome (LAHPS) is characterized by bleeding and thrombosis in patients with autoimmune diseases or infections. Paediatric LAHPS exhibits various degrees of bleeding, ranging from mild to severe; however, adrenal haemorrhage due to LAHPS and its long-term clinical course have not been sufficiently described. CASE PRESENTATION: A 9-year-old boy presented with prolonged abdominal pain and abnormal coagulation screening tests. The laboratory tests showed prolonged activated partial thromboplastin time and subsequently revealed the presence of lupus anticoagulant, anti-nuclear antibodies, and hypoprothrombinemia, leading to diagnosis of LAHPS. An enhanced computed tomogram demonstrated nodular lesions in the adrenal glands bilaterally, suggestive of adrenal haemorrhage. Laboratory and clinical manifestations exhibited life-threatening adrenal insufficiency that required hydrocortisone administration. The patient developed systemic lupus erythematosus, diagnosed 12 months later. CONCLUSIONS: This patient with LAHPS developed rare adrenal failure due to adrenal haemorrhage, a life-threatening event that should be recognized and treated early. In our case, renal dysfunction was also observed when systemic lupus erythematosus was diagnosed 1 year after LAHPS. Our case emphasizes that early recognition of adrenal failure and careful long-term observation is required in patients with autoantibodies.

NDNF variants are rare in patients with congenital hypogonadotropic hypogonadism.

Although NDNF was recently reported as a novel causative gene for congenital hypogonadotropic hypogonadism (CHH), this conclusion has yet to be validated. In this study, we sequenced NDNF in 61 Japanese CHH patients. No variants, except for nine synonymous substitutions that appear to have no effect on splice-site recognition, were identified in NDNF coding exons or flanking intronic sequences. These results indicate the rarity of NDNF variants in CHH patients and highlight the genetic heterogeneity of CHH.

45,X/46,X,psu idic(Y)(q11.2) in a phenotypically normal male with short stature: a case report.

We report a case of 15-yr-old phenotypically normal male with short stature associated with the chromosomal abnormalities of 46,X,psu idic(Y)(q11.2)/45,X. At 3 yr of age, he underwent urethroplasty for scrotal hypospadias. At 15 yr of age, he was referred to our hospital due to short stature (-3.71 SD). The results of blood examination were mostly normal. A radiological examination revealed his bone age was 15.7 yr (based on the TW2-RUS method). Chromosome analysis of peripheral lymphocytes revealed 46,X,psu idic(Y)(q11.2)[16]/45,X[14], and array comparative genomic hybridization (aCGH) showed a large deletion of Yq which was located distal to the Y chromosome growth-control gene (GCY) region. It is likely that these structural abnormalities in the Y chromosome were responsible for the short stature. This case might provide new insights regarding GCY and emphasizes the importance of chromosome analysis in not only females but also males with short stature, especially when associated with genital anomalies.

Rare variant of the epigenetic regulator SMCHD1 in a patient with pituitary hormone deficiency.

Isolated hypogonadotropic hypogonadism (IHH), combined pituitary hormone deficiency (CPHD), and septo-optic dysplasia (SOD) constitute a disease spectrum whose etiology remains largely unknown. This study aimed to clarify whether mutations in SMCHD1, an epigenetic regulator gene, might underlie this disease spectrum. SMCHD1 is a causative gene for Bosma arhinia microphthalmia syndrome characterized by arhinia, microphthalmia and IHH. We performed mutation screening of SMCHD1 in patients with etiology-unknown IHH (n = 31) or CPHD (n = 43, 19 of whom also satisfied the SOD diagnostic criteria). Rare variants were subjected to in silico analyses and classified according to the American College of Medical Genetics and Genomics guidelines. Consequently, a rare likely pathogenic variant, p.Asp398Asn, was identified in one patient. The patient with p.Asp398Asn exhibited CPHD, optic nerve hypoplasia, and a thin retinal nerve fiber layer, and therefore satisfied the criteria of SOD. This patient showed a relatively low DNA methylation level of the 52 SMCHD1-target CpG sites at the D4Z4 locus. Exome sequencing for the patient excluded additional variants in other IHH/CPHD-causative genes. In vitro assays suggested functional impairment of the p.Asp398Asn variant. These results provide the first indication that SMCHD1 mutations represent a rare genetic cause of the HH-related disease spectrum.

An Automatic Image Collection System for Multicenter Clinical Studies.

The acquisition of medical images from multiple medial institutions has become important for high-quality clinical studies. In recent years, electronic data submission has enabled the transmission of image data to independent institutions more quickly and easily than before. However, the selection, anonymization, and transmission of medical images still require human resources in the form of clinical research collaborators. In this study, we developed an image collection system that works with the electronic data capture (EDC) system. In this image collection system, medical images are selected based on EDC input information, patient ID is anonymized to a subject ID issued by the EDC, and the selected anonymized images are transferred to the research institute without human intervention. In the research institute, clinical information registered by the EDC and clinical images collected by the image collection system are managed by the same subject ID and can be used for clinical studies. In October 2019, our image collection system was introduced to 13 medical institutions and has now begun collecting medical images from the in-hospital picture archiving and communication system (PACS) of those institutions.

Established and Novel Mechanisms Leading to de novo Genomic Rearrangements in the Human Germline.

During gametogenesis, the human genome can acquire various de novo rearrangements. Most constitutional genomic rearrangements are created through 1 of the 4 well-known mechanisms, i.e., nonallelic homologous recombination, erroneous repair after double-strand DNA breaks, replication errors, and retrotransposition. However, recent studies have identified 2 types of extremely complex rearrangements that cannot be simply explained by these mechanisms. The first type consists of chaotic structural changes in 1 or a few chromosomes that result from "chromoanagenesis (an umbrella term that covers chromothripsis, chromoanasynthesis, and chromoplexy)." The other type is large independent rearrangements in multiple chromosomes indicative of "transient multifocal genomic crisis." Germline chromoanagenesis (chromothripsis) likely occurs predominantly during spermatogenesis or postzygotic embryogenesis, while multifocal genomic crisis appears to be limited to a specific time window during oogenesis and early embryogenesis or during spermatogenesis. This review article introduces the current understanding of the molecular basis of de novo rearrangements in the germline.

Unbalanced Y;7 Translocation between Two Low-Similarity Sequences Leading to SRY-Positive 45,X Testicular Disorders of Sex Development.

Unbalanced translocations of Y-chromosomal fragments harboring the sex-determining region Y gene (SRY) to the X chromosome or an autosome result in 46,XX and 45,X testicular disorders of sex development (DSD), respectively. Of these, Y;autosome translocation is an extremely rare condition. Here, we identified a 20-year-old man with a 45,X,t(Y;7)(q11.21;q35) karyotype, who exhibited unilateral cryptorchidism, small testis, intellectual disability, and various congenital anomalies. The fusion junction of the translocation was blunt, and the breakpoint-flanking regions shared only 50% similarity. These results indicate that Y;autosome translocations can occur between 2 low-similarity sequences, probably via nonhomologous end joining. Furthermore, translocations of a Ypterq11.21 fragment to 7q35 likely result in normal or only mildly impaired male-type sexual development, along with various clinical features of 7q deletion syndrome, although their effects on adult testicular function remain to be studied.

DNA Methylation Status of SHOX-Flanking CpG Islands in Healthy Individuals and Short Stature Patients with Pseudoautosomal Copy Number Variations.

SHOX resides in the short arm pseudoautosomal region (PAR1) of the sex chromosomes and escapes X inactivation. SHOX haploinsufficiency underlies idiopathic short stature (ISS) and Leri-Weill dyschondrosteosis (LWD). A substantial percentage of cases with SHOX haploinsufficiency arise from pseudoautosomal copy number variations (CNVs) involving putative enhancer regions of SHOX. Our previous study using peripheral blood samples showed that some CpG dinucleotides adjacent to SHOX exon 1 were hypomethylated in a healthy woman and methylated in a woman with gross X chromosomal rearrangements. However, it remains unknown whether submicroscopic pseudoautosomal CNVs cause aberrant DNA methylation of SHOX-flanking CpG islands. In this study, we examined the DNA methylation status of SHOX-flanking CpG islands in 50 healthy individuals and 10 ISS/LWD patients with pseudoautosomal CNVs. In silico analysis detected 3 CpG islands within the 20-kb region from the translation start site of SHOX. Pyrosequencing and bisulfite sequencing of genomic DNA samples revealed that these CpG islands were barely methylated in peripheral blood cells and cultured chondrocytes of healthy individuals, as well as in peripheral blood cells of ISS/LWD patients with pseudoautosomal CNVs. These results, in conjunction with our previous findings, indicate that the DNA methylation status of SHOX-flanking CpG islands can be affected by gross X-chromosomal abnormalities, but not by submicroscopic CNVs in PAR1. Such CNVs likely disturb SHOX expression through DNA methylation-independent mechanisms, which need to be determined in future studies.

Transient multifocal genomic crisis creating chromothriptic and non-chromothriptic rearrangements in prezygotic testicular germ cells.

BACKGROUND: The co-occurrence of multiple de novo copy number variations (CNVs) is a rare phenomenon in the human genome. Recently, an "organismal CNV mutator phenotype" has been reported to result in transient genomic instability introducing multiple de novo CNVs in primary oocytes and early-stage zygotes. These findings opened a new area of human genome research. METHODS: We performed genome-wide copy number analysis for ~ 2100 individuals with various congenital defects. Furthermore, extensive molecular analyses, including synthetic long-read whole-genome sequencing and haplotype-phasing, were carried out for an individual with multiple de novo CNVs. RESULTS: A boy was found to have de novo rearrangements on five chromosomes. The rearrangements comprised simple duplication and inversion as well as chaotic changes, all of which affected paternally derived chromosomes. Postzygotic genomic instability was ruled out. The duplicated regions on 6q and 13q contained both diallelic and triallelic loci, indicating that the genomic rearrangements were initially created during premeiotic mitosis and subsequently modified by physiological cross-over during meiosis I. Breakpoints of the rearrangements were indicative of non-homologous end joining, replication-based errors, and/or chromothripsis. The mutagenic event was independent of specific local DNA motifs or de novo point mutations, but may be driven by spermatogenesis-specific factors. CONCLUSIONS: These results indicate that during spermatogenesis, a transient multifocal genomic crisis can introduce several chromothriptic and non-chromothriptic changes into the genome. These findings broaden the concept of the "organismal CNV mutator phenotype". This study provides insights into mechanisms for altering the global chromosomal architecture of human embryos.