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.

Development of sexual dimorphism of skeletal muscles through the adrenal cortex, caused by androgen-induced global gene suppression.

The zona fasciculata (zF) in the adrenal cortex contributes to multiple physiological actions through glucocorticoid synthesis. The size, proliferation, and glucocorticoid synthesis characteristics are all female biased, and sexual dimorphism is established by androgen. In this study, transcriptomes were obtained to unveil the sex differentiation mechanism. Interestingly, both the amount of mRNA and the expressions of nearly all genes were higher in females. The expression of Nr5a1, which is essential for steroidogenic cell differentiation, was also female biased. Whole-genome studies demonstrated that NR5A1 regulates nearly all gene expression directly or indirectly. This suggests that androgen-induced global gene suppression is potentially mediated by NR5A1. Using Nr5a1 heterozygous mice, whose adrenal cortex is smaller than the wild type, we demonstrated that the size of skeletal muscles is possibly regulated by glucocorticoid synthesized by zF. Taken together, considering the ubiquitous presence of glucocorticoid receptors, our findings provide a pathway for sex differentiation through glucocorticoid synthesis.

Variations in gender identity and sexual orientation of university students.

Background: Previous studies have shown that a small percentage of people in the general population have atypical gender identity and/or sexual orientation. Aim: This study aimed to explore variations in gender identity and sexual orientation in university students and determine genetic factors associated with these variations. Methods: Deviations from complete gender congruence and exclusive heterosexual orientation in 736 Japanese university students were quantitatively assessed with self-assessment questionnaires. Next, we conducted genetic tests for 80 participants who showed relatively low gender identity scores and/or atypical sexual orientation. These genetic tests consisted of repeat number analysis of the androgen receptor gene (AR) and a SKAT-O: an optimal unified sequence kernel association test, which is an exome-based rare variant association study. The results of the genetic tests were compared with the Japanese reference data and the results of our 637 control samples. Outcomes: We calculated the gender identity and sexual orientation scores of all participants and analyzed the molecular data of 80 selected participants. Results: The gender identity scores of 736 participants were broadly distributed: only ~15% of natal males and ~5% of natal females had the maximum score that corresponds to complete gender congruence. The sexual orientation scores also varied: ~80% of natal males and ~60% of natal females showed exclusive heterosexual orientation. We found no association between gender characteristics and AR repeat numbers. The SKAT-O showed that rare damaging variants of TDRP and 3 other genes were more common in the 80 participants than in the control group. Clinical Implications: Our data support the view that gender is a phenotypic continuum rather than a binary trait. Strength and Limitations: This study quantitatively assessed the gender characteristics of a large cohort of university students. Moreover, we conducted systematic screening for genetic factors associated with gender variations. The weaknesses of the study were the limited analytic power of the questionnaires, the relatively small sample for molecular analyses, and incomplete clinical information and relatively advanced ages of the control group. Conclusion: This study revealed significant variations in gender identity and sexual orientation in university students, which may be partly associated with variants in TDRP or other genes.

A capture methyl-seq protocol with improved efficiency and cost-effectiveness using pre-pooling and enzymatic conversion.

OBJECTIVE: The opportunities for sequencing-based methylome analysis of clinical samples are increasing. To reduce its cost and the amount of genomic DNA required for library preparation, we aimed to establish a capture methyl-seq protocol, which adopts pre-pooling of multiple libraries before hybridization capture and TET2/APOBEC-mediated conversion of unmethylated cytosine to thymine. RESULTS: We compared a publicly available dataset generated by the standard Agilent protocol of SureSelect XT Human Methyl-Seq Kit and our dataset obtained by our modified protocol, EMCap, that adopted sample pre-pooling and enzymatic conversion. We confirmed that the quality of DNA methylation data was comparable between the two datasets. As our protocol, EMCap, is more cost-effective and reduces the amount of input genomic DNA, it would serve as a better choice for clinical methylome sequencing.

Prevalence of pathogenic variants in cancer-predisposing genes in second cancer after childhood solid cancers.

BACKGROUND: Second malignant neoplasms (SMNs) are one of the most severe late complications after pediatric cancer treatment. However, the effect of genetic variation on SMNs remains unclear. In this study, we revealed germline genetic factors that contribute to the development of SMNs after treatment of pediatric solid tumors. METHODS: We performed whole-exome sequencing in 14 pediatric patients with SMNs, including three brain tumors. RESULTS: Our analysis revealed that five of 14 (35.7%) patients had pathogenic germline variants in cancer-predisposing genes (CPGs), which was significantly higher than in the control cohort (p < 0.01). The identified genes with variants were TP53 (n = 2), DICER1 (n = 1), PMS2 (n = 1), and PTCH1 (n = 1). In terms of the type of subsequent cancer, leukemia and multiple episodes of SMN had an exceptionally high rate of CPG pathogenic variants. None of the patients with germline variants had a family history of SMN development. Mutational signature analysis showed that platinum drugs contributed to the development of SMN in three cases, which suggests the role of platinum agents in SMN development. CONCLUSIONS: We highlight that overlapping effects of genetic background and primary cancer treatment contribute to the development of second cancers after treatment of pediatric solid tumors. A comprehensive analysis of germline and tumor samples may be useful to predict the risk of secondary cancers.

Familial Pseudohypoparathyroidism Type IB Associated with an SVA Retrotransposon Insertion in the GNAS Locus.

Loss of methylation (LOM) at GNAS-A/B:TSS-differentially methylated regions (DMRs) in the GNAS locus is observed in pseudohypoparathyroidism type 1B (PHP1B). Many PHP1B cases are sporadic, but autosomal dominant-PHP1B has a deletion involving NESP55 expressed from the maternal allele or STX16 located upstream of the GNAS locus on the maternal allele. We report the possible first familial PHP1B cases with retrotransposon insertion in the GNAS locus on the maternal allele. To our knowledge, they are the possible first cases with imprinting disorders caused by retrotransposon insertion. The two sibling cases experienced tetany and/or cramps from school age and had hypocalcemia and an increased serum intact parathyroid hormone (PTH) level together with overweight, round face, and normal intellectual levels. Methylation analysis for DMRs in the GNAS locus showed only LOM of the GNAS-A/B:TSS-DMR. Copy number abnormalities at STX16 and the GNAS locus were not detected by array comparative genomic hybridization. Whole-genome sequencing and Sanger sequencing revealed an approximately 1000-bp SVA retrotransposon insertion upstream of the first exon of A/B on the GNAS locus in these siblings. Whole-genome methylome analysis by Enzymatic Methyl-Seq in the siblings showed normal methylation status in the region surrounding the insertion site and mild LOM of the GNAS-A/B:TSS-DMR. We conducted transcriptome analysis using mRNA from skin fibroblasts and induced pluripotent stem cells (iPSCs) derived from the siblings and detected no aberrant NESP55 transcripts. Quantitative reverse-transcriptase PCR (qRT-PCR) analysis in skin fibroblasts showed increased A/B expression in the patients and no NESP55 expression, even in a control. qRT-PCR analysis in iPSCs showed decreased NESP55 expression with normal methylation status of the GNAS-NESP:TSS-DMR in the patients. The retrotransposon insertion in the siblings likely caused decreased NESP55 expression that could lead to increased A/B expression via LOM of the GNAS-A/B:TSS-DMR, subsequent reduced Gsα expression, and finally, PHP1B development. © 2022 American Society for Bone and Mineral Research (ASBMR).

Quantitative assessment of copy number alterations by liquid biopsy for neuroblastoma.

Liquid biopsy, a method of detecting genomic alterations using blood specimens, has recently attracted attention as a noninvasive alternative to surgical tissue biopsy. We attempted quantitative analysis to detect amplification of MYCN (MYCNamp) and loss of heterozygosity at 11q (11qLOH), which are clinical requisites as prognostic factors of neuroblastoma (NB). In this study, cell-free DNA (cfDNA) was extracted from plasma samples from 24 NB patients at diagnosis. Copy numbers of MYCN and NAGK genes were quantitatively analyzed by droplet digital PCR (ddPCR). 11qLOH was also assessed by detecting allelic imbalances of heterozygous single nucleotide polymorphisms in the 11q region. The results obtained were compared to those of specimens from tumor tissues. The correlation coefficient of MYCN copy number of cfDNA and tumor DNA was 0.88 (p < 0.00001). 11qLOH was also accurately detected from cfDNA, except for one case with localized NB. Given the high accuracy of liquid biopsy, to investigate components of cfDNA, the proportion of tumor-derived DNA was estimated by examining the variant allele frequency of tumor-specific mutations in cfDNA. The proportion of tumor-derived DNA in cfDNA was 42.5% (range, 16.9%-55.9%), suggesting sufficient sensitivity of liquid biopsy for NB. In conclusion, MYCN copy number and 11qLOH could be quantitatively analyzed in plasma cfDNA by ddPCR assay. These results suggest that plasma cfDNA can be substituted for tumor DNA and can also be applied for comprehensive genomic profiling analysis.

Genetic features of precursor B-cell phenotype Burkitt leukemia with IGH-MYC rearrangement.

BACKGROUND: An atypical form of Burkitt leukemia/lymphoma (BL), BL with a phenotype of precursor B-cells (preBLL), is listed in the WHO Classification. Recent reports suggested that preBLL and classical BL could be distinguished by the differences in IG-MYC translocation architecture and an additional mutated genes profile. The characteristics of classical BL are IG-MYC by aberrant somatic hypermutation or class switch recombination, and BL-specific gene mutations such as MYC, ID3, and CCND3. Meanwhile, preBLL is characterized by IG-MYC due to aberrant VDJ recombination and mutations in NRAS and KRAS. However, it is not clear whether all preBLL cases can be differentiated. This report investigated the molecular characteristics of an infant preBLL case, with a more advanced stage of maturity than typical preBLL. CASE: The patient showed BL-like morphology with IGH-MYC rearrangement. In the immunophenotyping, CD20 and surface immunoglobulin were negative, whereas other markers were consistent with BL. To evaluate the genetic contribution, we performed whole-exome sequencing. The breakpoint analysis revealed the IG-MYC occurred due to an aberrant VDJ recombination. Meanwhile, additional somatic mutations were detected in FBXO11, one of the mutant genes specific to BL. In the analysis of the specimen in complete remission, mutation in KRAS, frequently mutated in preBLL, was detected with low frequency, suggesting somatic mosaicism. CONCLUSION: The present case showed the characteristics of both typical preBLL and classical BL. Because preBLL includes atypical cases such as the present case, further studies are required to elucidate preBLL features.

NUDT15 variants confer high incidence of second malignancies in children with acute lymphoblastic leukemia.

The effect of genetic variation on second malignant neoplasms (SMNs) remains unclear. First, we identified the pathogenic germline variants in cancer-predisposing genes among 15 children with SMNs after childhood leukemia/lymphoma using whole-exome sequencing. Because the prevalence was low, we focused on the association between SMNs and NUDT15 in primary acute lymphoblastic leukemia (ALL) cases. NUDT15 is one of the 6-mercaptopurine (6-MP) metabolic genes, and its variants are common in East Asian individuals. The prevalence of NUDT15 hypomorphic variants was higher in patients with SMNs (n = 14; 42.9%) than in the general population in the gnomAD database (19.7%; P = .042). In the validation study with a cohort of 438 unselected patients with ALL, the cumulative incidence of SMNs was significantly higher among those with (3.0%; 95% confidence interval [CI], 0.6% to 9.4%) than among those without NUDT15 variants (0.3%; 95% CI, 0.0% to 1.5%; P = .045). The 6-MP dose administered to patients with ALL with a NUDT15 variant was higher than that given to those without SMNs (P = .045). The 6-MP-related mutational signature was observed in SMN specimens after 6-MP exposure. In cells exposed to 6-MP, a higher level of 6-MP induced DNA damage in NUDT15-knockdown induced pluripotent stem cells. Our study indicates that NUDT15 variants may confer a risk of SMNs after treatment with 6-MP in patients with ALL.

Characteristics of genetic alterations of peripheral T-cell lymphoma in childhood including identification of novel fusion genes: the Japan Children's Cancer Group (JCCG).

Peripheral T-cell lymphoma (PTCL) is a group of heterogeneous non-Hodgkin lymphomas showing a mature T-cell or natural killer cell phenotype, but its molecular abnormalities in paediatric patients remain unclear. By employing next-generation sequencing and multiplex ligation-dependent probe amplification of tumour samples from 26 patients, we identified somatic alterations in paediatric PTCL including Epstein-Barr virus (EBV)-negative (EBV- ) and EBV-positive (EBV+ ) patients. As recurrent mutational targets for PTCL, we identified several previously unreported genes, including TNS1, ZFHX3, LRP2, NCOA2 and HOXA1, as well as genes previously reported in adult patients, e.g. TET2, CDKN2A, STAT3 and TP53. However, for other reported mutations, VAV1-related abnormalities were absent and mutations of NRAS, GATA3 and JAK3 showed a low frequency in our cohort. Concerning the association of EBV infection, two novel fusion genes: STAG2-AFF2 and ITPR2-FSTL4, and deletion and alteration of CDKN2A/2B, LMO1 and HOXA1 were identified in EBV- PTCL, but not in EBV+ PTCL. Conversely, alterations of PCDHGA4, ADAR, CUL9 and TP53 were identified only in EBV+ PTCL. Our observations suggest a clear difference in the molecular mechanism of onset between paediatric and adult PTCL and a difference in the characteristics of genetic alterations between EBV- and EBV+ paediatric PTCL.

Three nuclear protein-coding genes corroborate a recent phylogenomic model of the Branchiopoda (Crustacea) and provide estimates of the divergence times of the major branchiopodan taxa.

The class Branchiopoda (Crustacea) shows great diversity in morphology and lifestyle among its constituent higher-level taxa: Anostraca, Notostraca, Laevicaudata, Spinicaudata, Cyclestherida and Cladocera. The phylogenetic relationships among these taxa have long been controversial. We sequenced three orthologous nuclear genes that encode the catalytic subunit of DNA polymerase delta and the largest and second-largest subunits of RNA polymerase II in the expectation that the amino acid sequences encoded by these genes might be effective in clarifying branchiopod phylogeny and estimating the times of divergence of the major branchiopodan taxa. The results of phylogenetic analyses based on these amino acid sequences support the monophyly of Branchiopoda and provide strong molecular evidence in support of the following phylogenetic relationships: (Anostraca, (Notostraca, (Laevicaudata, (Spinicaudata, (Cyclestherida, Cladocera))))). Within Cladocera, comparison of the nucleotide sequences of these same genes shows Ctenopoda to be the sister group of Haplopoda + Anomopoda. Three statistical tests based on the present amino acid sequence data-the approximately unbiased test, Kishino-Hasegawa test and weighted Shimodaira-Hasegawa test-tend to refute most of the previous molecular phylogenetic studies on Branchiopoda, which have placed Notostraca differently than here; however, our results corroborate those of one recent phylogenomic study, thus confirming the effectiveness of these three genes to investigate relationships among branchiopod higher taxa. Divergence time estimates calibrated on the basis of fossil evidence suggest that the first divergence of extant branchiopods occurred about 534 Ma during the early Cambrian period and that diversification within the extant branchiopod lineages started in or after the late Permian.

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.

Molecular phylogenetic analyses support the monophyly of Hexapoda and suggest the paraphyly of Entognatha.

BACKGROUND: Molecular phylogenetic analyses have revealed that Hexapoda and Crustacea form a common clade (the Pancrustacea), which is now widely accepted among zoologists; however, the origin of Hexapoda remains unresolved. The main problems are the unclear relationships among the basal hexapod lineages, Protura (proturans), Collembola (springtails), Diplura (diplurans), and Ectognatha (bristletails, silverfishes, and all winged insects). Mitogenomic analyses have challenged hexapod monophyly and suggested the reciprocal paraphyly of Hexapoda and Crustacea, whereas studies based on nuclear molecular data support the monophyletic origin of hexapods. Additionally, there are significant discrepancies with respect to these issues between the results of morphological and molecular studies. To investigate these problems, we performed phylogenetic analyses of Pancrustacea based on the protein sequences of three orthologous nuclear genes encoding the catalytic subunit of DNA polymerase delta and the largest and second largest subunits of RNA polymerase II from 64 species of arthropods, including representatives of all hexapod orders. RESULTS: Phylogenetic analyses were conducted based on the inferred amino acid (aa) sequences (~3400 aa in total) of the three genes using the maximum likelihood (ML) method and Bayesian inference. Analyses were also performed with additional datasets generated by excluding long-branch taxa or by using different outgroups. These analyses all yielded essentially the same results. All hexapods were clustered into a common clade, with Branchiopoda as its sister lineage, whereas Crustacea was paraphyletic. Within Hexapoda, the lineages Ectognatha, Palaeoptera, Neoptera, Polyneoptera, and Holometabola were each confirmed to be monophyletic with robust support, but monophyly was not supported for Entognatha (Protura + Collembola + Diplura), Ellipura (Protura + Collembola), or Nonoculata (Protura + Diplura). Instead, our results showed that Protura is the sister lineage to all other hexapods and that Diplura or Diplura + Collembola is closely related to Ectognatha. CONCLUSION: This is the first study to include all hexapod orders in a phylogenetic analysis using multiple nuclear protein-coding genes to investigate the phylogeny of Hexapoda, with an emphasis on Entognatha. The results strongly support the monophyletic origin of hexapods but reject the monophyly of Entognatha, Ellipura, and Nonoculata. Our results provided the first molecular evidence in support of Protura as the sister group to other hexapods. These findings are expected to provide additional insights into the origin of hexapods and the processes involved in the adaptation of insects to life on land.

Phylogenetic relationships among insect orders based on three nuclear protein-coding gene sequences.

Many attempts to resolve the phylogenetic relationships of higher groups of insects have been made based on both morphological and molecular evidence; nonetheless, most of the interordinal relationships of insects remain unclear or are controversial. As a new approach, in this study we sequenced three nuclear genes encoding the catalytic subunit of DNA polymerase delta and the two largest subunits of RNA polymerase II from all insect orders. The predicted amino acid sequences (In total, approx. 3500 amino acid sites) of these proteins were subjected to phylogenetic analyses based on the maximum likelihood and Bayesian analysis methods with various models. The resulting trees strongly support the monophyly of Palaeoptera, Neoptera, Polyneoptera, and Holometabola, while within Polyneoptera, the groupings of Isoptera/"Blattaria"/Mantodea (Superorder Dictyoptera), Dictyoptera/Zoraptera, Dermaptera/Plecoptera, Mantophasmatodea/Grylloblattodea, and Embioptera/Phasmatodea are supported. Although Paraneoptera is not supported as a monophyletic group, the grouping of Phthiraptera/Psocoptera is robustly supported. The interordinal relationships within Holometabola are well resolved and strongly supported that the order Hymenoptera is the sister lineage to all other holometabolous insects. The other orders of Holometabola are separated into two large groups, and the interordinal relationships of each group are (((Siphonaptera, Mecoptera), Diptera), (Trichoptera, Lepidoptera)) and ((Coleoptera, Strepsiptera), (Neuroptera, Raphidioptera, Megaloptera)). The sister relationship between Strepsiptera and Diptera are significantly rejected by all the statistical tests (AU, KH and wSH), while the affinity between Hymenoptera and Mecopterida are significantly rejected only by AU and KH tests. Our results show that the use of amino acid sequences of these three nuclear genes is an effective approach for resolving the relationships of higher groups of insects.

A viable Bacillus subtilis strain without functional extracytoplasmic function sigma genes.

We constructed a Bacillus subtilis Marburg strain that harbors deletion mutations in all seven extracytoplasmic function (ECF) sigma genes. The strain shows wild-type growth at 37 degrees C both in a complex and in a synthetic medium and exhibits wild-type sporulation. ECF sigma genes of B. subtilis are dispensable as long as no stress is imposed, although they seem to be required for quick response to stresses.