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.

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.

Basophil-derived interleukin-4 controls the function of natural helper cells, a member of ILC2s, in lung inflammation.

Allergic asthma is an inflammatory disease characterized by lung eosinophilia controlled by type 2 cytokines. Cysteine proteases are potent triggers of allergic inflammation by causing barrier disruption in lung epithelial cells inducing the elevation of interleukin-5 (IL-5) and IL-13 from natural helper (NH) cells, a member of ILC2s, which leads to lung eosinophilia. In this study, we found that basophils play a crucial role in NH cell-mediated eosinophilic inflammation induced by protease allergens. Conditional deletion of basophils caused a resolution of the papain-induced eosinophilia and mucus production. Resolution of eosinophilia was also observed in mice lacking IL-4 specifically in basophils, indicating that basophil-derived IL-4 enhanced expression of the chemokine CCL11, as well as IL-5, IL-9, and IL-13 in NH cells, thus attracting eosinophils. These results demonstrate that IL-4 from basophils has an important role in the NH-derived cytokine and chemokine expression, subsequently leading to protease allergen-induced airway inflammation.

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.

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.

Genome evolution in the allotetraploid frog Xenopus laevis.

To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We characterize the allotetraploid origin of X. laevis by partitioning its genome into two homoeologous subgenomes, marked by distinct families of 'fossil' transposable elements. On the basis of the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged around 34 million years ago (Ma) and combined to form an allotetraploid around 17-18 Ma. More than 56% of all genes were retained in two homoeologous copies. Protein function, gene expression, and the amount of conserved flanking sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.

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.

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.

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.

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.

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.

Identification of CD34+ and CD34- leukemia-initiating cells in MLL-rearranged human acute lymphoblastic leukemia.

Translocation of the mixed-lineage leukemia (MLL) gene with AF4, AF9, or ENL results in acute leukemia with both lymphoid and myeloid involvement. We characterized leukemia-initiating cells (LICs) in primary infant MLL-rearranged leukemia using a xenotransplantation model. In MLL-AF4 patients, CD34(+)CD38(+)CD19(+) and CD34(-)CD19(+) cells initiated leukemia, and in MLL-AF9 patients, CD34(-)CD19(+) cells were LICs. In MLL-ENL patients, either CD34(+) or CD34(-) cells were LICs, depending on the pattern of CD34 expression. In contrast, in patients with these MLL translocations, CD34(+)CD38(-)CD19(-)CD33(-) cells were enriched for normal hematopoietic stem cells (HSCs) with in vivo long-term multilineage hematopoietic repopulation capacity. Although LICs developed leukemic cells with clonal immunoglobulin heavy-chain (IGH) rearrangement in vivo, CD34(+)CD38(-)CD19(-)CD33(-) cells repopulated recipient bone marrow and spleen with B cells, showing broad polyclonal IGH rearrangement and recipient thymus with CD4(+) single positive (SP), CD8(+) SP, and CD4(+)CD8(+) double-positive (DP) T cells. Global gene expression profiling revealed that CD9, CD32, and CD24 were over-represented in MLL-AF4, MLL-AF9, and MLL-ENL LICs compared with normal HSCs. In patient samples, these molecules were expressed in CD34(+)CD38(+) and CD34(-) LICs but not in CD34(+)CD38(-)CD19(-)CD33(-) HSCs. Identification of LICs and LIC-specific molecules in primary human MLL-rearranged acute lymphoblastic leukemia may lead to improved therapeutic strategies for MLL-rearranged leukemia.

Ancestral Y-linked genes were maintained by translocation to the X and Y chromosomes fused to an autosomal pair in the Okinawa spiny rat Tokudaia muenninki.

Two species of the genus Tokudaia lack the Y chromosome and SRY, but several Y-linked genes have been rescued by translocation or transposition to other chromosomes. Tokudaia muenninki is the only species in the genus that maintains the Y owing to sex chromosome-autosome fusions. According to previous studies, many SRY pseudocopies and other Y-linked genes have evolved by excess duplication in this species. Using RNA-seq and RT-PCR, we found that ZFY, EIF2S3Y, TSPY, UTY, DDX3Y, USP9Y, and RBMY, but not UBA1Y, had high deduced amino acid sequence similarity and similar expression patterns with other rodents, suggesting that these genes were functional. Based on FISH and quantitative real-time PCR, all of the genes except for UTY and DDX3Y were amplified on the X and Y chromosomes with approximately 10-66 copies in the male genome. In a comparative analysis of the 372.4-kb BAC sequence and Y-linked gene transcripts from T. muenninki with the mouse Y genomic sequence, we observed that multiple-copy genes in the ancestral Y genome were nonfunctional, indicating that the gene functions were assumed by amplified copies. We also found a LTR sequence at the distal end of a SRY duplication unit, suggesting that unequal sister chromatid exchange mediated by retrotransposable elements could have been involved in SRY amplification. Our results revealed that the Y-linked genes were rescued from degeneration via translocations to other sex chromosomal regions and amplification events in T. muenninki.

Initiation of recombination suppression and PAR formation during the early stages of neo-sex chromosome differentiation in the Okinawa spiny rat, Tokudaia muenninki.

BACKGROUND: Sex chromosomes of extant eutherian species are too ancient to reveal the process that initiated sex-chromosome differentiation. By contrast, the neo-sex chromosomes generated by sex-autosome fusions of recent origin in Tokudaia muenninki are expected to be evolutionarily 'young', and therefore provide a good model in which to elucidate the early phases of eutherian sex chromosome evolution. Here we describe the genomic evolution of T. muenninki in neo-sex chromosome differentiation. RESULTS: FISH mapping of a T. muenninki male, using 50 BAC clones as probes, revealed no chromosomal rearrangements between the neo-sex chromosomes. Substitution-direction analysis disclosed that sequence evolution toward GC-richness, which positively correlates with recombination activity, occurred in the peritelomeric regions, but not middle regions of the neo-sex chromosomes. In contrast, the sequence evolution toward AT-richness was observed in those pericentromeric regions. Furthermore, we showed genetic differentiation between the pericentromeric regions as well as an accelerated rate of evolution in the neo-Y region through the detection of male-specific substitutions by gene sequencing in multiple males and females, and each neo-sex-derived BAC sequencing. CONCLUSIONS: Our results suggest that recombination has been suppressed in the pericentromeric region of neo-sex chromosomes without chromosome rearrangement, whereas high levels of recombination activity is limited in the peritelomeric region of almost undifferentiated neo-sex chromosomes. We conclude that PAR might have been formed on the peritelomeric region of sex chromosomes as an independent event from spread of recombination suppression during the early stages of sex chromosome differentiation.

Comparative genome and transcriptome analyses of the social amoeba Acytostelium subglobosum that accomplishes multicellular development without germ-soma differentiation.

BACKGROUND: Social amoebae are lower eukaryotes that inhabit the soil. They are characterized by the construction of a starvation-induced multicellular fruiting body with a spore ball and supportive stalk. In most species, the stalk is filled with motile stalk cells, as represented by the model organism Dictyostelium discoideum, whose developmental mechanisms have been well characterized. However, in the genus Acytostelium, the stalk is acellular and all aggregated cells become spores. Phylogenetic analyses have shown that it is not an ancestral genus but has lost the ability to undergo cell differentiation. RESULTS: We performed genome and transcriptome analyses of Acytostelium subglobosum and compared our findings to other available dictyostelid genome data. Although A. subglobosum adopts a qualitatively different developmental program from other dictyostelids, its gene repertoire was largely conserved. Yet, families of polyketide synthase and extracellular matrix proteins have not expanded and a serine protease and ABC transporter B family gene, tagA, and a few other developmental genes are missing in the A. subglobosum lineage. Temporal gene expression patterns are astonishingly dissimilar from those of D. discoideum, and only a limited fraction of the ortholog pairs shared the same expression patterns, so that some signaling cascades for development seem to be disabled in A. subglobosum. CONCLUSIONS: The absence of the ability to undergo cell differentiation in Acytostelium is accompanied by a small change in coding potential and extensive alterations in gene expression patterns.

Evolutionary history of novel genes on the tammar wallaby Y chromosome: Implications for sex chromosome evolution.

We report here the isolation and sequencing of 10 Y-specific tammar wallaby (Macropus eugenii) BAC clones, revealing five hitherto undescribed tammar wallaby Y genes (in addition to the five genes already described) and several pseudogenes. Some genes on the wallaby Y display testis-specific expression, but most have low widespread expression. All have partners on the tammar X, along with homologs on the human X. Nonsynonymous and synonymous substitution ratios for nine of the tammar XY gene pairs indicate that they are each under purifying selection. All 10 were also identified as being on the Y in Tasmanian devil (Sarcophilus harrisii; a distantly related Australian marsupial); however, seven have been lost from the human Y. Maximum likelihood phylogenetic analyses of the wallaby YX genes, with respective homologs from other vertebrate representatives, revealed that three marsupial Y genes (HCFC1X/Y, MECP2X/Y, and HUWE1X/Y) were members of the ancestral therian pseudoautosomal region (PAR) at the time of the marsupial/eutherian split; three XY pairs (SOX3/SRY, RBMX/Y, and ATRX/Y) were isolated from each other before the marsupial/eutherian split, and the remaining three (RPL10X/Y, PHF6X/Y, and UBA1/UBE1Y) have a more complex evolutionary history. Thus, the small marsupial Y chromosome is surprisingly rich in ancient genes that are retained in at least Australian marsupials and evolved from testis-brain expressed genes on the X.

Comparative analysis of chimpanzee and human Y chromosomes unveils complex evolutionary pathway.

The mammalian Y chromosome has unique characteristics compared with the autosomes or X chromosomes. Here we report the finished sequence of the chimpanzee Y chromosome (PTRY), including 271 kb of the Y-specific pseudoautosomal region 1 and 12.7 Mb of the male-specific region of the Y chromosome. Greater sequence divergence between the human Y chromosome (HSAY) and PTRY (1.78%) than between their respective whole genomes (1.23%) confirmed the accelerated evolutionary rate of the Y chromosome. Each of the 19 PTRY protein-coding genes analyzed had at least one nonsynonymous substitution, and 11 genes had higher nonsynonymous substitution rates than synonymous ones, suggesting relaxation of selective constraint, positive selection or both. We also identified lineage-specific changes, including deletion of a 200-kb fragment from the pericentromeric region of HSAY, expansion of young Alu families in HSAY and accumulation of young L1 elements and long terminal repeat retrotransposons in PTRY. Reconstruction of the common ancestral Y chromosome reflects the dynamic changes in our genomes in the 5-6 million years since speciation.

The amphioxus genome and the evolution of the chordate karyotype.

Lancelets ('amphioxus') are the modern survivors of an ancient chordate lineage, with a fossil record dating back to the Cambrian period. Here we describe the structure and gene content of the highly polymorphic approximately 520-megabase genome of the Florida lancelet Branchiostoma floridae, and analyse it in the context of chordate evolution. Whole-genome comparisons illuminate the murky relationships among the three chordate groups (tunicates, lancelets and vertebrates), and allow not only reconstruction of the gene complement of the last common chordate ancestor but also partial reconstruction of its genomic organization, as well as a description of two genome-wide duplications and subsequent reorganizations in the vertebrate lineage. These genome-scale events shaped the vertebrate genome and provided additional genetic variation for exploitation during vertebrate evolution.