Long-read next-generation sequencing for molecular diagnosis of pediatric endocrine disorders.

Recent advances in long-read next-generation sequencing (NGS) have enabled researchers to identify several pathogenic variants overlooked by short-read NGS, array-based comparative genomic hybridization, and other conventional methods. Long-read NGS is particularly useful in the detection of structural variants and repeat expansions. Furthermore, it can be used for mutation screening in difficultto- sequence regions, as well as for DNA-methylation analyses and haplotype phasing. This mini-review introduces the usefulness of long-read NGS in the molecular diagnosis of pediatric endocrine disorders.

Isodicentric Y chromosome with multiple breakpoints in the pseudoautosomal region 1.

INTRODUCTION: Isodicentric Y chromosomes are relatively common structural variants of the human genome. The underlying mechanism of isodicentric Y chromosomes with short arm breakpoints [idic(Yq)] remains to be clarified. CASE PRESENTATION: We encountered a Japanese man with azoospermia and mild short stature. G-banding and array-based comparative genomic hybridization indicated that his karyotype was 45,X/46,X,idic(Y)(qter→p11.32::p11.32→qter) with a ~1.8 Mb terminal deletion. Whole genome sequencing suggested that the Y chromosome had four breakpoints in a ~7 kb region of the pseudoautosomal region 1 (PAR1). CONCLUSION: This case was assumed to have an idic(Yq) resulting from multiple DNA double-strand breaks in PAR1. This rearrangement may have been facilitated by the PAR1-specific chromatin architecture. The clinical features of the patient can be ascribed to SHOX haploinsufficiency and the presence of a 45,X cell line, although copy-number gains of some Yq genes and the size reduction of PAR1 may also contribute to his spermatogenic failure.

Exome-wide benchmark of difficult-to-sequence regions using short-read next-generation DNA sequencing.

Next-generation DNA sequencing (NGS) in short-read mode has recently been used for genetic testing in various clinical settings. NGS data accuracy is crucial in clinical settings, and several reports regarding quality control of NGS data, primarily focusing on establishing NGS sequence read accuracy, have been published thus far. Variant calling is another critical source of NGS errors that remains unexplored at the single-nucleotide level despite its established significance. In this study, we used a machine-learning-based method to establish an exome-wide benchmark of difficult-to-sequence regions at the nucleotide-residue resolution using 10 genome sequence features based on real-world NGS data accumulated in The Genome Aggregation Database (gnomAD) of the human reference genome sequence (GRCh38/hg38). The newly acquired metric, designated the 'UNMET score,' along with additional lines of structural information from the human genome, allowed us to assess the sequencing challenges within the exonic region of interest using conventional short-read NGS. Thus, the UNMET score could provide a basis for addressing potential sequential errors in protein-coding exons of the human reference genome sequence GRCh38/hg38 in clinical sequencing.

Y Chromosome Genomic Variations and Biological Significance in Human Diseases and Health.

The Y chromosome is a haploid genome unique to males with no genes essential for life. It is easily transmitted to the next generation without being repaired by recombination, even if a major genomic structural alteration occurs. On the other hand, the Y chromosome genome is basically a region transmitted only from father to son, reflecting a male-specific inheritance between generations. The Y chromosome exhibits genomic structural differences among different ethnic groups and individuals. The Y chromosome was previously thought to affect only male-specific phenotypes, but recent studies have revealed associations between the Y chromosomes and phenotypes common to both males and females, such as certain types of cancer and neuropsychiatric disorders. This evidence was discovered with the finding of the mosaic loss of the Y chromosome in somatic cells. This phenomenon is also affected by environmental factors, such as smoking and aging. In the past, functional analysis of the Y chromosome has been elucidated by assessing the function of Y chromosome-specific genes and the association between Y chromosome haplogroups and human phenotypes. These studies are currently being conducted intensively. Additionally, the recent advance of large-scale genome cohort studies has increased the amount of Y chromosome genomic information available for analysis, making it possible to conduct more precise studies of the relationship between genome structures and phenotypes. In this review, we will introduce recent analyses using large-scale genome cohort data and previously reported association studies between Y chromosome haplogroups and human phenotypes, such as male infertility, cancer, cardiovascular system traits, and neuropsychiatric disorders. The function and biological role of the Y chromosome in human phenotypes will also be discussed.

Genetic variants of G-protein coupled receptors associated with pubertal disorders.

Background: The human hypothalamic-pituitary-gonadal (HPG) axis is the regulatory center for pubertal development. This axis involves six G-protein coupled receptors (GPCRs) encoded by KISS1R, TACR3, PROKR2, GNRHR, LHCGR, and FSHR. Methods: Previous studies have identified several rare variants of the six GPCR genes in patients with pubertal disorders. In vitro assays and animal studies have provided information on the function of wild-type and variant GPCRs. Main Findings: Of the six GPCRs, those encoded by KISS1R and TACR3 are likely to reside at the top of the HPG axis. Several loss-of-function variants in the six genes were shown to cause late/absent puberty. In particular, variants in KISS1R, TACR3, PROKR2, and GNRHR lead to hypogonadotropic hypogonadism in autosomal dominant, recessive, and oligogenic manners. Furthermore, a few gain-of-function variants of KISS1R, PROKR2, and LHCGR have been implicated in precocious puberty. The human HPG axis may contain additional GPCRs. Conclusion: The six GPCRs in the HPG axis govern pubertal development through fine-tuning of hormone secretion. Rare sequence variants in these genes jointly account for a certain percentage of genetic causes of pubertal disorders. Still, much remains to be clarified about the molecular network involving the six GPCRs.

Identification of a New Enhancer That Promotes Sox9 Expression by a Comparative Analysis of Mouse and Sry-Deficient Amami Spiny Rat.

INTRODUCTION: Testis differentiation is initiated by the SRY gene on the Y chromosome in mammalian species. However, the Amami spiny rat, Tokudaia osimensis, lacks both the Y chromosome and the Sry gene and acquired a unique Sox9 regulatory mechanism via a male-specific duplication upstream of Sox9, without Sry. In general mammalian species, the SRY protein binds to a testis-specific enhancer to promote SOX9 gene expression. Several enhancers located upstream of Sox9/SOX9 have been reported in mice and humans. In particular, the binding of SRY to the highly conserved enhancer Enh13 is thought to be a common mechanism underlying testis differentiation and sex determination in mammals. METHODS: Sequences of T. osimensis homologues of three Sox9 enhancers that were previously reported in mice, Enh8, Enh14, and Enh13, were determined. We performed in vitro assays to confirm enhancer activity involved in Sox9 regulation in T. osimensis. RESULTS: T. osimensis Enh13 showed enhancer activity when co-transfected with NR5A1 and SOX9. Mouse Enh13 was activated by NR5A1 and SRY; however, T. osimensis Enh13 did not respond to SRY, even though the binding sites of SRY and NR5A1 were conserved. To identify the key sequence that is present in mouse but absent from T. osimensis, we performed reporter gene assays using vectors in which partial sequences of T. osimensis Enh13 were replaced with mouse sequences. For T. osimensis Enh13 in which the second half (approximately 430 bp) was replaced with the corresponding mouse sequence, activity in response to NR5A1 and SRY was recovered. Further, reporter assays revealed that multiple regions in the second half of the mouse Enh13 sequence are required for the response to NR5A1 and SRY. The latter 49 bp was particularly important and contained four binding sites for three transcription factors, POU2F1, HOXA3, and GATA1. CONCLUSION: We showed that there are unknown sequences responsible for the interaction between NR5A1 and SRY and mEnh13 based on comparative analyses of Sry-dependent and Sry-independent species. Our comparative analyses revealed new molecular mechanisms underlying mammalian sex determination.

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.

The Neo-X Does Not Form a Barr Body but Shows a Slightly Condensed Structure in the Okinawa Spiny Rat (Tokudaia muenninki).

X chromosome inactivation (XCI) is an essential mechanism for gene dosage compensation between male and female cells in mammals. The Okinawa spiny rat (Tokudaia muenninki) is a native rodent in Japan with XX/XY sex chromosomes, like most mammals; however, the X chromosome has acquired a neo-X region (Xp) by fusion with an autosome. We previously reported that dosage compensation has not yet evolved in the neo-X region; however, X-inactive-specific transcript (Xist) RNA (long non-coding RNA required for the initiation of XCI) is partially localized in the region. Here, we show that the neo-X region represents an early chromosomal state in the acquisition of XCI by analyses of heterochromatin and Barr body formation. We found no evidence for heterochromatin formation in the neo-X region by R-banding by acridine orange (RBA) assays and immunostaining of H3K27me3. Double-immunostaining of H3K27me3 and HP1, a component of the Barr body, revealed that the entire ancestral X chromosome region (Xq) showed a bipartite folded structure. By contrast, HP1 was not localized to the neo-X region. However, BAC-FISH revealed that the signals of genes on the neo-X region of the inactive X chromosome were concentrated in a narrow region. These findings indicated that although the neo-X region of the inactive X chromosome does not form a complete Barr body structure (e.g., it lacks HP1), it forms a slightly condensed structure. These findings combined with the previously reported partial binding of Xist RNA suggest that the neo-X region exhibits incomplete inactivation. This may represent an early chromosomal state in the acquisition of the XCI mechanism.

Biodistribution and radiation dosimetry of the positron emission tomography probe for AMPA receptor, [11C]K-2, in healthy human subjects.

[11C]K-2, a radiotracer exhibiting high affinity and selectivity for α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptors (AMPARs), is suitable for the quantification of AMPARs in living human brains and potentially useful in the identification of epileptogenic foci in patients. This study aimed to estimate the radiation doses of [11C]K-2 in various organs and calculate the effective dose after injection of [11C]K-2 in healthy human subjects. Twelve healthy male subjects were registered and divided into two groups (370 or 555 MBq of [11C]K-2), followed by 2 h whole-body scans. We estimated the radiation dose of each organ and then calculated the effective dose for each subject. The highest uptake of [11C]K-2 was observed in the liver, while the brain also showed relatively high uptake. The urinary bladder exhibited the highest radiation dose. The kidneys and liver also showed high radiation doses after [11C]K-2 injections. The effective dose of [11C]K-2 ranged from 5.0 to 5.2 µSv/MBq. Our findings suggest that [11C]K-2 is safe in terms of the radiation dose and adverse effects. The injection of 370-555 MBq (10 to 15 mCi) for PET studies using this radiotracer is applicable in healthy human subjects and enables serial PET scans in a single subject.

Suppression of Non-Random Fertilization by MHC Class I Antigens.

Hermaphroditic invertebrates and plants have a self-recognition system on the cell surface of sperm and eggs, which prevents their self-fusion and enhances non-self-fusion, thereby contributing to genetic variation. However, the system of sperm-egg recognition in mammals is under debate. To address this issue, we explored the role of major histocompatibility complex class I (MHC class I, also known as histocompatibility 2-Kb or H2-Kb and H2-Db in mice) antigens by analyzing H2-Kb-/-H2-Db-/-ß2-microglobulin (ß2M)-/- triple-knockout (T-KO) male mice with full fertility. T-KO sperm exhibited an increased sperm number in the perivitelline space of wild-type (WT) eggs in vitro. Moreover, T-KO sperm showed multiple fusion with zona pellucida (ZP)-free WT eggs, implying that the ability of polyspermy block for sperm from T-KO males was weakened in WT eggs. When T-KO male mice were intercrossed with WT female mice, the percentage of females in progeny increased. We speculate that WT eggs prefer fusion with T-KO sperm, more specifically X-chromosome-bearing sperm (X sperm), suggesting the presence of preferential (non-random) fertilization in mammals, including humans.

Correction: Extra-mitochondrial citrate synthase initiates calcium oscillation and suppresses age-dependent sperm dysfunction.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Visualization of AMPA receptors in living human brain with positron emission tomography.

Although aberrations in the number and function of glutamate AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptors are thought to underlie neuropsychiatric disorders, no methods are currently available for visualizing AMPA receptors in the living human brain. Here we developed a positron emission tomography (PET) tracer for AMPA receptors. A derivative of 4-[2-(phenylsulfonylamino)ethylthio]-2,6-difluoro-phenoxyacetamide radiolabeled with 11C ([11C]K-2) showed specific binding to AMPA receptors. Our clinical trial with healthy human participants confirmed reversible binding of [11C]K-2 in the brain according to Logan graphical analysis (UMIN000020975; study design: non-randomized, single arm; primary outcome: dynamics and distribution volumes of [11C]K-2 in the brain; secondary outcome: adverse events of [11C]K-2 during the 4-10 d following dosing; this trial met prespecified endpoints). In an exploratory clinical study including patients with epilepsy, we detected increased [11C]K-2 uptake in the epileptogenic focus of patients with mesial temporal lobe epilepsy, which was closely correlated with the local AMPA receptor protein distribution in surgical specimens from the same individuals (UMIN000025090; study design: non-randomized, single arm; primary outcome: correlation between [11C]K-2 uptake measured with PET before surgery and AMPA receptor protein density examined by biochemical study after surgery; secondary outcome: adverse events during the 7 d following PET scan; this trial met prespecified endpoints). Thus, [11C]K-2 is a potent PET tracer for AMPA receptors, potentially providing a tool to examine the involvement of AMPA receptors in neuropsychiatric disorders.

Extra-mitochondrial citrate synthase initiates calcium oscillation and suppresses age-dependent sperm dysfunction.

Men and women become infertile with age, but the mechanism of declining male fertility, more specifically, the decrease in in sperm quality, is not well known. Citrate synthase (CS) is a core enzyme of the mitochondrial tricarboxylic acid (TCA) cycle, which directly controls cellular function. Extra-mitochondrial CS (eCS) is produced and abundant in the sperm head; however, its role in male fertility is unknown. We investigated the role of eCS in male fertility by producing eCs-deficient (eCs-KO) mice. The initiation of the first spike of Ca2+ oscillation was substantially delayed in egg fused with eCs-KO sperm, despite normal expression of sperm factor phospholipase C zeta 1. The eCs-KO male mice were initially fertile, but the fertility dropped with age. Metabolomic analysis of aged sperm revealed that the loss of eCS enhances TCA cycle in the mitochondria with age, presumably leading to depletion of extra-mitochondrial citrate. The data suggest that eCS suppresses age-dependent male infertility, providing insights into the decline of male fertility with age.

Construction of JRG (Japanese reference genome) with single-molecule real-time sequencing.

In recent genome analyses, population-specific reference panels have indicated important. However, reference panels based on short-read sequencing data do not sufficiently cover long insertions. Therefore, the nature of long insertions has not been well documented. Here, we assembled a Japanese genome using single-molecule real-time sequencing data and characterized insertions found in the assembled genome. We identified 3691 insertions ranging from 100 bps to ~10,000 bps in the assembled genome relative to the international reference sequence (GRCh38). To validate and characterize these insertions, we mapped short-reads from 1070 Japanese individuals and 728 individuals from eight other populations to insertions integrated into GRCh38. With this result, we constructed JRGv1 (Japanese Reference Genome version 1) by integrating the 903 verified insertions, totaling 1,086,173 bases, shared by at least two Japanese individuals into GRCh38. We also constructed decoyJRGv1 by concatenating 3559 verified insertions, totaling 2,536,870 bases, shared by at least two Japanese individuals or by six other assemblies. This assembly improved the alignment ratio by 0.4% on average. These results demonstrate the importance of refining the reference assembly and creating a population-specific reference genome. JRGv1 and decoyJRGv1 are available at the JRG website.

Construction of full-length Japanese reference panel of class I HLA genes with single-molecule, real-time sequencing.

Human leukocyte antigen (HLA) is a gene complex known for its exceptional diversity across populations, importance in organ and blood stem cell transplantation, and associations of specific alleles with various diseases. We constructed a Japanese reference panel of class I HLA genes (ToMMo HLA panel), comprising a distinct set of HLA-A, HLA-B, HLA-C, and HLA-H alleles, by single-molecule, real-time (SMRT) sequencing of 208 individuals included in the 1070 whole-genome Japanese reference panel (1KJPN). For high-quality allele reconstruction, we developed a novel pipeline, Primer-Separation Assembly and Refinement Pipeline (PSARP), in which the SMRT sequencing and additional short-read data were used. The panel consisted of 139 alleles, which were all extended from known IPD-IMGT/HLA sequences, contained 40 with novel variants, and captured more than 96.5% of allelic diversity in 1KJPN. These newly available sequences would be important resources for research and clinical applications including high-resolution HLA typing, genetic association studies, and analyzes of cis-regulatory elements.

Genome analyses for the Tohoku Medical Megabank Project towards establishment of personalized healthcare.

Personalized healthcare (PHC) based on an individual's genetic make-up is one of the most advanced, yet feasible, forms of medical care. The Tohoku Medical Megabank (TMM) Project aims to combine population genomics, medical genetics and prospective cohort studies to develop a critical infrastructure for the establishment of PHC. To date, a TMM CommCohort (adult general population) and a TMM BirThree Cohort (birth+three-generation families) have conducted recruitments and baseline surveys. Genome analyses as part of the TMM Project will aid in the development of a high-fidelity whole-genome Japanese reference panel, in designing custom single-nucleotide polymorphism (SNP) arrays specific to Japanese, and in estimation of the biological significance of genetic variations through linked investigations of the cohorts. Whole-genome sequencing from >3,500 unrelated Japanese and establishment of a Japanese reference genome sequence from long-read data have been done. We next aim to obtain genotype data for all TMM cohort participants (>150,000) using our custom SNP arrays. These data will help identify disease-associated genomic signatures in the Japanese population, while genomic data from TMM BirThree Cohort participants will be used to improve the reference genome panel. Follow-up of the cohort participants will allow us to test the genetic markers and, consequently, contribute to the realization of PHC.

Androgen and Oestrogen Affect the Expression of Long Non-Coding RNAs During Phallus Development in a Marsupial.

There is increasing evidence that long non-coding RNAs (lncRNAs) are important for normal reproductive development, yet very few lncRNAs have been identified in phalluses so far. Unlike eutherians, phallus development in the marsupial tammar wallaby occurs post-natally, enabling manipulation not possible in eutherians in which differentiation occurs in utero. We treated with sex steroids to determine the effects of androgen and oestrogen on lncRNA expression during phallus development. Hormonal manipulations altered the coding and non-coding gene expression profile of phalluses. We identified several predicted co-regulatory lncRNAs that appear to be co-expressed with the hormone-responsive candidate genes regulating urethral closure and phallus growth, namely IGF1, AR and ESR1. Interestingly, more than 50% of AR-associated coding genes and lncRNAs were also associated with ESR1. In addition, we identified and validated three novel co-regulatory and hormone-responsive lncRNAs: lnc-BMP5, lnc-ZBTB16 and lncRSPO4. Lnc-BMP5 was detected in the urethral epithelium of male phalluses and was downregulated by oestrogen in males. Lnc-ZBTB16 was downregulated by oestrogen treatment in male phalluses at day 50 post-partum (pp). LncRSPO4 was downregulated by adiol treatment in female phalluses but increased in male phalluses after castration. Thus, the expression pattern and hormone responsiveness of these lncRNAs suggests a physiological role in the development of the phallus.

Biallelic mutations in SZT2 cause a discernible clinical entity with epilepsy, developmental delay, macrocephaly and a dysmorphic corpus callosum.

Mutations in SZT2 were first reported in 2013 as a cause of early-onset epileptic encephalopathy. Because only five reports have been published to date, the clinical features associated with SZT2 remain unclear. We herein report an additional patient with biallelic mutations in SZT2. The proband, a four-year-old girl, showed developmental delay and seizures from two years of age. Her seizures were not intractable and readily controlled by valproate. She showed mildly dysmorphic facies with macrocephaly, high forehead, and hypertelorism, and also had pectus carinatum. An EEG showed epileptic discharges which rarely occurred. A brain MRI revealed a short and thick corpus callosum. Whole-exome sequencing detected compound heterozygous biallelic mutations (c.8596dup (p.Tyr2866Leufs∗42) and c.2930-17_2930-3delinsCTCGTG) in SZT2, both of which were novel and predicted to be truncating. This case suggested a broad phenotypic spectrum arises from SZT2 mutations, forming a continuum from epileptic encephalopathy and severe developmental delay to mild intellectual disability without epilepsy. The characteristic thick and short corpus callosum observed in 7/8 cases with epilepsy, including the proband, but not in three non-syndromic cases, appears to be specific, and thus useful for indicating the possibility of SZT2 mutations. This feature has the potential to make loss of SZT2 a clinically discernible disorder despite a wide clinical spectrum.

Identification of somatic mutations in postmortem human brains by whole genome sequencing and their implications for psychiatric disorders.

AIM: Somatic mutations in the human brain are hypothesized to contribute to the functional diversity of brain cells as well as the pathophysiology of neuropsychiatric diseases. However, there are still few reports on somatic mutations in non-neoplastic human brain tissues. This study attempted to unveil the landscape of somatic mutations in the human brain. METHODS: We explored the landscape of somatic mutations in human brain tissues derived from three individuals with no neuropsychiatric diseases by whole-genome deep sequencing at a depth of around 100. The candidate mutations underwent multi-layered filtering, and were validated by ultra-deep target amplicon sequencing at a depth of around 200 000. RESULTS: Thirty-one somatic mutations were identified in the human brain, demonstrating the utility of whole-genome sequencing of bulk brain tissue. The mutations were enriched in neuron-expressed genes, and two-thirds of the identified somatic single nucleotide variants in the brain tissues were cytosine-to-thymine transitions, half of which were in CpG dinucleotides. CONCLUSION: Our developed filtering and validation approaches will be useful to identify somatic mutations in the human brain. The vulnerability of neuron-expressed genes to mutational events suggests their potential relevance to neuropsychiatric diseases.