Search

Next-generation DNA sequencing (NGS) in short-read mode has been recently 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, focusing mostly on establishing NGS sequence read accuracy, have been published thus far. Variant calling is another critical source of NGS errors that remains mostly unexplored 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 using 10 genome sequence features on the basis of real-world NGS data accumulated in The Genome Aggregation Database (gnomAD) of the human reference genome sequence (GRCh38/hg38). We used the obtained metrics, designated “UNMET score,” along with other lines of structural information of the human genome to identify difficult-to-sequence genomic regions using conventional NGS. Thus, the UNMET score could provide appropriate caveats to address potential sequential errors in protein-coding exons of the human reference genome sequence GRCh38/hg38 in clinical sequencing.

[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.

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., Citrate synthase is a core enzyme of tricarboxylic acid cycle, but the role of extra-mitochondrial CS (eCS) is less clear. eCS, abundant in the sperm head, triggers a first spike of egg Ca2+ oscillation in a phospholipase C zeta 1 -independent manner. Moreover, the fertility of eCs-knockout male mice dropped with age, indicating that eCS suppresses age-dependent sperm dysfunction by accelerating Ca2+ oscillation.

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., Reference genome: Reading longer sequences improves Japanese reference Researchers in Japan have assembled a Japanese reference genome, which includes sequences missing from the international reference genome, as well as others specific to East Asian populations. A team led by Masao Nagasaki and Masayuki Yamamoto sequenced a Japanese individual using a method, which produces longer sequences than previous technologies. Using this approach, they identified thousands of sequences spanning 2.5 million bases, which were absent in the international reference genome. Many of these were sequences able to move within the genome. They showed that the majority of these sequences are also present in early humans and chimpanzees, demonstrating that their absence from the current reference is due to deletions or limitations of earlier sequencing methodologies. In addition to providing a population-specific reference, these findings demonstrate the importance of continually improving the international reference genome.

[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.

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.

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. A newly developed PET tracer allows visualization of AMPA receptors in the living human brain, providing a new tool to study their potential involvement in neurological or psychiatric disorders.

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&ndash, 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-/-&beta, 2-microglobulin (&beta, 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.

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.

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.

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.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.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.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.

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.

Reproductive success depends on a robust and appropriately timed preovulatory luteinizing hormone (LH) surge, which is induced by the activation of gonadotropin-releasing hormone (GnRH) neurons in response to positive feedback from increasing estrogen levels. Here we document an increase in postsynaptic GluR2-lacking Ca2+-permeable AMPA-type glutamate receptors (CP-AMPARs) at synapses on GnRH neurons on the day of proestrus in rats, coincident with the increase in estrogen levels. Functional blockade of CP-AMPARs depressed the synaptic responses only on the day of proestrus and concomitantly attenuated the LH surge. Thus, the phasic synaptic incorporation of postsynaptic CP-AMPARs on GnRH neurons is involved in the generation of the LH surge.

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.

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.

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.

In humans and mice, there are 11 genes derived from sushi-ichi related retrotransposons, some of which are known to play essential roles in placental development. Interestingly, this family of retrotransposons was thought to exist only in eutherian mammals, indicating their significant contributions to the eutherian evolution, but at least one, PEG10, is conserved between marsupials and eutherians. Here we report a novel sushi-ichi retrotransposon-derived gene, SIRH12, in the tammar wallaby, an Australian marsupial species of the kangaroo family. SIRH12 encodes a protein highly homologous to the sushi-ichi retrotransposon Gag protein in the tammar wallaby, while SIRH12 in the South American short-tailed grey opossum is a pseudogene degenerated by accumulation of multiple nonsense mutations. This suggests that SIRH12 retrotransposition occurred only in the marsupial lineage but acquired and retained some as yet unidentified novel function, at least in the lineage of the tammar wallaby.

Quantitative phosphoproteome and transcriptome analysis of ligand-stimulated MCF-7 human breast cancer cells was performed to understand the mechanisms of tamoxifen resistance at a system level. Phosphoproteome data revealed that WT cells were more enriched with phospho-proteins than tamoxifen-resistant cells after stimulation with ligands. Surprisingly, decreased phosphorylation after ligand perturbation was more common than increased phosphorylation. In particular, 17β-estradiol induced down-regulation in WT cells at a very high rate. 17β-Estradiol and the ErbB ligand heregulin induced almost equal numbers of up-regulated phospho-proteins in WT cells. Pathway and motif activity analyses using transcriptome data additionally suggested that deregulated activation of GSK3β (glycogen-synthase kinase 3β) and MAPK1/3 signaling might be associated with altered activation of cAMP-responsive element-binding protein and AP-1 transcription factors in tamoxifen-resistant cells, and this hypothesis was validated by reporter assays. An examination of clinical samples revealed that inhibitory phosphorylation of GSK3β at serine 9 was significantly lower in tamoxifen-treated breast cancer patients that eventually had relapses, implying that activation of GSK3β may be associated with the tamoxifen-resistant phenotype. Thus, the combined phosphoproteome and transcriptome data set analyses revealed distinct signal transcription programs in tumor cells and provided a novel molecular target to understand tamoxifen resistance.

We constructed and characterized arrayed bacterial artificial chromosome (BAC) libraries of five Drosophila species (D. melanogaster, D. simulans, D. sechellia, D. auraria, and D. ananassae), which are genetically well characterized in the studies of meiosis, evolution, population genetics, and developmental biology. The BAC libraries comprise 8,000 to 12,500 clones for each species, estimated to cover the most of the genomes. We sequenced both ends of most of these BAC clones with a success rate of 91%. Of these, 53,701 clones consisting of non-repetitive BAC end sequences (BESs) were mapped with reference of the public D. melanogaster genome sequences. The BES mapping estimated that the BAC libraries of D. auraria and D. ananassae covered 47% and 57% of the D. melanogaster genome, respectively, and those of D. melanogaster, D. sechellia, and D. simulans covered 94-97%. The low coverage by BESs of D. auraria and D. ananassae may be due to the high sequence divergence with D. melanogaster. From the comparative BES mapping, 111 possible breakpoints of chromosomal rearrangements were identified in these four species. The breakpoints of the major chromosome rearrangement between D. simulans and D. melanogaster on the third chromosome were determined within 20 kb in 84E and 30 kb in 93E/F. Corresponding breakpoints were also identified in D. sechellia. The BAC clones described here will be an important addition to the Drosophila genomic resources.

The laboratory rat is one of the most extensively studied model organisms. Inbred laboratory rat strains originated from limited Rattus norvegicus founder populations, and the inherited genetic variation provides an excellent resource for the correlation of genotype to phenotype. Here, we report a survey of genetic variation based on almost 3 million newly identified SNPs. We obtained accurate and complete genotypes for a subset of 20,238 SNPs across 167 distinct inbred rat strains, two rat recombinant inbred panels and an F2 intercross. Using 81% of these SNPs, we constructed high-density genetic maps, creating a large dataset of fully characterized SNPs for disease gene mapping. Our data characterize the population structure and illustrate the degree of linkage disequilibrium. We provide a detailed SNP map and demonstrate its utility for mapping of quantitative trait loci. This community resource is openly available and augments the genetic tools for this workhorse of physiological studies. This work was supported by European Union grants LSGH-2004-005235 and LSHG-CT-2005-019015. D.G. is supported by a Wellcome Trust Senior Fellowship in Basic Biomedical Science (057733/Z/99/A). M.-T.B. and D.G. acknowledge support from the Wellcome Cardiovascular Functional Genomics Initiative (066780/Z/01/Z)

Medaka (Oryzias latipes) is a small egg-laying freshwater teleost native to East Asia that has become an excellent model system for developmental genetics and evolutionary biology. The draft medaka genome sequence (700 Mb) was reported in June 2007, and its substantial genomic resources have been opened to the public through the University of Tokyo Genome Browser Medaka (UTGB/medaka) database. This database provides basic genomic information, such as predicted genes, expressed sequence tags (ESTs), guanine/cytosine (GC) content, repeats and comparative genomics, as well as unique data resources including (i) 2473 genetic markers and experimentally confirmed PCR primers that amplify these markers, (ii) 142,414 bacterial artificial chromosome (BAC) and 217,344 fosmid end sequences that amount to 15.0- and 11.1-fold clone coverage of the entire genome, respectively, and were used for draft genome assembly, (iii) 16,519,460 single nucleotide polymorphisms (SNPs), and 2 859 905 insertions/deletions detected between two medaka inbred strain genomes and (iv) 841 235 5'-end serial analyses of gene-expression (SAGE) tags that identified 344 266 transcription start sites on the genome. UTGB/medaka is available at: http://medaka.utgenome.org/.

The first draft of the common marmoset (Callithrix jacchus) genome was published by the Marmoset Genome Sequencing and Analysis Consortium. The draft was based on whole-genome shotgun sequencing and the current assembly version is Callithrix_jacches-3.2.1, but there still exist 187,214 undetermined gap regions and supercontigs and relatively short contigs that are unmapped to chromosomes in the draft genome. We performed resequencing and assembly of the genome of common marmoset by deep sequencing with high-throughput sequencing technology. Several different sequence runs using Illumina sequencing platforms were executed and 181 Gbp of high-quality bases including mate-pairs with long insert lengths of 3, 8, 20 and 40 Kbp were obtained, that is, approximately 60× coverage. The resequencing significantly improved the MGSAC draft genome sequence. The N50 of the contigs, which is a statistical measure used to evaluate assembly quality, doubled. As a result, 51% of the contigs (total length: 299 Mbp) that were unmapped to chromosomes in the MGSAC draft were merged with chromosomal contigs and the improved genome sequence helped to detect 5,288 new genes that are homologous to human cDNAs and the gaps in 5,187 transcripts of the Ensembl gene annotations were completely filled.

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.

The Tohoku Medical Megabank Organization reports the whole-genome sequences of 1,070 healthy Japanese individuals and construction of a Japanese population reference panel (1KJPN). Here we identify through this high-coverage sequencing (32.4 × on average), 21.2 million, including 12 million novel, single-nucleotide variants (SNVs) at an estimated false discovery rate of, The Tohoku Medical Megabank Organization establishes a biobank with detailed patient health care and genome information. Here the authors analyse whole-genome sequences of 1,070 Japanese individuals, allowing them to catalogue 21 million single-nucleotide variants including 12 million novel ones.

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.

Chromosome 18 appears to have the lowest gene density of any human chromosome and is one of only three chromosomes for which trisomic individuals survive to term. There are also a number of genetic disorders stemming from chromosome 18 trisomy and aneuploidy. Here we report the finished sequence and gene annotation of human chromosome 18, which will allow a better understanding of the normal and disease biology of this chromosome. Despite the low density of protein-coding genes on chromosome 18, we find that the proportion of non-protein-coding sequences evolutionarily conserved among mammals is close to the genome-wide average. Extending this analysis to the entire human genome, we find that the density of conserved non-protein-coding sequences is largely uncorrelated with gene density. This has important implications for the nature and roles of non-protein-coding sequence elements.

SRY (sex-determining region Y) is widely conserved in eutherian mammals as a sex-determining gene located on the Y chromosome. SRY proteins bind to the testis-specific enhancer of SOX9 (TES) with SF1 to upregulate SOX9 expression in undifferentiated gonads of XY embryos of humans and mice. The core region within TES, named TESCO, is an important enhancer for mammalian sex determination. We show that TESCO of the genus Tokudaia lost enhancer activity caused by mutations in its SRY and SF1 binding sites. Two species of Tokudaia do not have the Y chromosome or SRY, and one species has multiple SRYs located on the neo-Y chromosome consisting of the Y fused with an autosome. The sequence of Tokudaia TESCO exhibited more than 83% identity with mouse TESCO, however, nucleotide substitution(s) were found in two out of three SRY binding sites and in five out of six SF1 binding sites. TESCO of all species showed low enhancer activity in cells co-transfected with SRY and SF1, and SOX9 and SF1 in reporter gene assays. Mutated TESCO, in which nucleotide substitutions found in SRY and SF1 binding sites were replaced with mouse sequence, recovered the activity. Furthermore, SRYs of the SRY-positive species could not activate the mutated TESCO or mouse TESCO, suggesting that SRYs lost function as a sex-determining gene any more. Our results indicate that the SRY dependent sex-determining mechanism was lost in a common ancestor of the genus Tokudaia caused by nucleotide substitutions in SRY and SF1 binding sites after emergence of a new sex-determining gene. We present the first evidence for an intermediate stage of the switchover from SRY to a new sex-determining gene in the evolution of mammalian sex-determining mechanism.

Sex chromosomes harbour a primary sex-determining signal that triggers sexual development of the organism. However, diverse sex chromosome systems have been evolved in vertebrates. Here we use positional cloning to identify the sex-determining locus of a medaka-related fish, Oryzias dancena, and find that the locus on the Y chromosome contains a cis-regulatory element that upregulates neighbouring Sox3 expression in developing gonad. Sex-reversed phenotypes in Sox3Y transgenic fish, and Sox3Y loss-of-function mutants all point to its critical role in sex determination. Furthermore, we demonstrate that Sox3 initiates testicular differentiation by upregulating expression of downstream Gsdf, which is highly conserved in fish sex differentiation pathways. Our results not only provide strong evidence for the independent recruitment of Sox3 to male determination in distantly related vertebrates, but also provide direct evidence that a novel sex determination pathway has evolved through co-option of a transcriptional regulator potentially interacted with a conserved downstream component. Sex chromosomes harbour specific sequences that determine the sexual development of the organism; yet these sequences remain unknown for many species. Here, Takehana et al. show that, similarly to mammals, Sox3 on the Y chromosome is the male-determining factor in the medaka-related fish Oryzias dancena.

Four cases of late onset clear cell renal cell carcinoma (RCC), a case of gastric cancer, and a case of exocrine pancreatic cancer were identified in a Japanese family. In order to elucidate the underlying mechanism for tumorigenesis in this family, extensive genetic studies were performed including routine and spectral karyotyping (SKY), fluorescence in situ hybridisation (FISH), comparative genomic hybridisation (CGH), loss of heterozygosity studies (LOH), and VHL mutation analysis. A germline translocation t(1;3)(q32-q41;q13-q21) was identified by karyotyping in five members of the family including all three RCC cases tested. The translocation was refined to t(1;3)(q32;q13.3) by FISH analysis using locus specific genomic clones, and the two breakpoints were mapped to a 5 cM region in 3q13.3 and a 3.6 cM region in 1q32. Both CGH and allelotyping using microsatellite markers showed loss of the derivative chromosome 3 carrying a 1q segment in the three familial RCCs analysed. Additional chromosomal imbalances were identified by CGH, including amplifications of chromosomes 5 and 7 and loss of 8p and 9. No germline VHL mutation was found but two different somatic mutations, a splice (IVS1-2A>C) and a frameshift (726delG), were identified in two RCCs from the same patient confirming their distinct origin.Taken together, these results firmly support a three step model for tumorigenesis in this family. A constitutional translocation t(1q;3q) increased the susceptibility to loss of the derivative chromosome 3 which is then followed by somatic mutations of the RCC related tumour suppressor gene VHL located in the remaining copy of chromosome 3. Keywords: familial renal cell carcinoma; translocation; von Hippel-Lindau disease; loss of heterozygosity

Two newly developed microsatellite markers on Yp11 (DXYS265) and Yq11.21 (DXYS266) and our previously reported marker, on Yp11 (DXYS241), were typed by triplex polymerase chain reaction (PCR) in 102 Japanese, 18 white American, and 17 black American males. The DXYS265 locus revealed three alleles, the DXYS266 locus showed two alleles, while the DXYS241 locus showed five alleles. Nine different compound haplotypes were observed among the males. Of these, two haplotypes were common to all groups, while four were limited to Japanese. Pedigree analysis of 61 Japanese families revealed no mutations of these loci. The triplex PCR developed in this study, as well as the new loci, are useful for tracing paternal lineages in human migration studies and population analysis, in addition to Y chromosome evolutionary studies.

A polymorphism in the coding sequence of the SRY gene was found by single-strand conformation polymorphism (SSCP) and direct sequencing analysis. The new allele of the SRY gene, which is raised by a C-to-T transition in the 155th codon, was found in 24% of Honshu, 35% of Okinawan, and 51% of Korean males respectively, whereas it was not observed among 16 Caucasian and 18 Negroid males. A haplotype analysis of the Y chromosome was carried out in Japanese, Korean, Caucasian and Negroid populations, using a combination of the polymorphisms in SRY, DXYS5Y, DYS287, and DXYS241Y loci. The results indicated that the Y chromosomes can be classified into seven heplotypes (Ia, Ib, Ic, IIa, IIb, III, IV). However, of these seven, only four (Ia, IIa, III, IV) were observed in the Japanese population. Furthermore, the presumed haplotype C, Y1, YAP, (CA)14, from which haplotype III was probably derived, was not found in any populations in this study. The regional distribution of each haplotype revealed that type III is more frequently observed in Okinawa (16%) and in Korea (21%) than in Honshu (4.4%). The haplotype analysis of the Y chromosome may contribute to the exploration of the origin of Japanese and the relationship between east Asian populations.

We describe a 28 year old male with a pseudodicentric Y chromosome who suffered from azoospermia attributed to maturation arrest of the primary spermatocyte, as diagnosed by testicular biopsy. Chromosome analysis, using G, Q and C banding techniques, revealed an abnormal karyotype of 45,X[7]/46,X,psu dic (Y)(pter-->q11.2::q11.2-->pter)[33]. Polymerase chain reaction (PCR) DNA analysis did not detect the absence of DAZ and RBM1 which are candidates for azoospermic factor (AZF) genes. Therefore, it is suggested that the maturation arrest of the primary spermatocyte in this patient was caused either by a pairing dysfunction between the X and Y chromosomes during meiosis or by deletions in the autosomal or the Y chromosomal spermatogenesis controlling genes, excluding DAZ and RBM1.

Embryogenesis in placental mammals is sustained by exquisite interplay between the embryo proper and placenta. UTF1 is a developmentally regulated gene expressed in both cell lineages. Here, we analyzed the consequence of loss of the UTF1 gene during mouse development. We found that homozygous UTF1 mutant newborn mice were significantly smaller than wild-type or heterozygous mutant mice, suggesting that placental insufficiency caused by the loss of UTF1 expression in extra-embryonic ectodermal cells at least in part contributed to this phenotype. We also found that the effects of loss of UTF1 expression in embryonic stem cells on their pluripotency were very subtle. Genome structure and sequence comparisons revealed that the UTF1 gene exists only in placental mammals. Our analyses of a family of genes with homology to UTF1 revealed a possible mechanism by which placental mammals have evolved the UTF1 genes.

The entire human protective protein gene has been cloned, and structural analysis revealed that the gene spans 7.5kb and comprises 15 exons. Furthermore, it partially overlaps on the opposite strand with the gene encoding phospholipid transfer protein. This region of DNA on chromosome 20 appears to encode two distinct mRNAs expressing defined functional products, and the mRNAs overlap by 58 nucleotides at their 3′-untranslated ends.

Three sex-determining (SD) genes, SRY (mammals), Dmy (medaka), and DM-W (Xenopus laevis), have been identified to date in vertebrates. However, how and why a new sex-determining gene appears remains unknown, as do the switching mechanisms of the master sex-determining gene. Here, we used positional cloning to search for the sex-determining gene in Oryzias luzonensis and found that GsdfY (gonadal soma derived growth factor on the Y chromosome) has replaced Dmy as the master sex-determining gene in this species. We found that GsdfY showed high expression specifically in males during sex differentiation. Furthermore, the presence of a genomic fragment that included GsdfY converts XX individuals into fertile XX males. Luciferase assays demonstrated that the upstream sequence of GsdfY contributes to the male-specific high expression. Gsdf is downstream of Dmy in the sex-determining cascade of O. latipes, suggesting that emergence of the Dmy-independent Gsdf allele led to the appearance of this novel sex-determining gene in O. luzonensis.

BACKGROUND: We present the genome sequence of the tammar wallaby, Macropus eugenii, which is a member of the kangaroo family and the first representative of the iconic hopping mammals that symbolize Australia to be sequenced. The tammar has many unusual biological characteristics, including the longest period of embryonic diapause of any mammal, extremely synchronized seasonal breeding and prolonged and sophisticated lactation within a well-defined pouch. Like other marsupials, it gives birth to highly altricial young, and has a small number of very large chromosomes, making it a valuable model for genomics, reproduction and development. RESULTS: The genome has been sequenced to 2 × coverage using Sanger sequencing, enhanced with additional next generation sequencing and the integration of extensive physical and linkage maps to build the genome assembly. We also sequenced the tammar transcriptome across many tissues and developmental time points. Our analyses of these data shed light on mammalian reproduction, development and genome evolution: there is innovation in reproductive and lactational genes, rapid evolution of germ cell genes, and incomplete, locus-specific X inactivation. We also observe novel retrotransposons and a highly rearranged major histocompatibility complex, with many class I genes located outside the complex. Novel microRNAs in the tammar HOX clusters uncover new potential mammalian HOX regulatory elements. CONCLUSIONS: Analyses of these resources enhance our understanding of marsupial gene evolution, identify marsupial-specific conserved non-coding elements and critical genes across a range of biological systems, including reproduction, development and immunity, and provide new insight into marsupial and mammalian biology and genome evolution.

The carotenoid-binding protein (CBP) of the domesticated silkworm, Bombyx mori, a major determinant of cocoon color, is likely to have been substantially influenced by domestication of this species. We analyzed the structure of the CBP gene in multiple strains of B. mori, in multiple individuals of the wild silkworm, B. mandarina (the putative wild ancestor of B. mori), and in a number of other lepidopterans. We found the CBP gene copy number in genomic DNA to vary widely among B. mori strains, ranging from 1 to 20. The copies of CBP are of several types, based on the presence of a retrotransposon or partial deletion of the coding sequence. In contrast to B. mori, B. mandarina was found to possess a single copy of CBP without the retrotransposon insertion, regardless of habitat. Several other lepidopterans were found to contain sequences homologous to CBP, revealing that this gene is evolutionarily conserved in the lepidopteran lineage. Thus, domestication can generate significant diversity of gene copy number and structure over a relatively short evolutionary time.

Background The glutamate receptors (GluRs) play a vital role in the mediation of excitatory synaptic transmission in the central nervous system. To clarify the evolutionary dynamics and mechanisms of the GluR genes in the lineage leading to humans, we determined the complete sequences of the coding regions and splice sites of 26 chimpanzee GluR genes. Results We found that all of the reading frames and splice sites of these genes reported in humans were completely conserved in chimpanzees, suggesting that there were no gross structural changes in humans after their divergence from the human-chimpanzee common ancestor. We observed low K A /K S ratios in both humans and chimpanzees, and we found no evidence of accelerated evolution. We identified 30 human-specific "fixed" amino acid substitutions in the GluR genes by analyzing 80 human samples of seven different populations worldwide. Grantham's distance analysis showed that GRIN2C and GRIN3A are the most and the second most diverged GluR genes between humans and chimpanzees. However, most of the substitutions are non-radical and are not clustered in any particular region. Protein motif analysis assigned 11 out of these 30 substitutions to functional regions. Two out of these 11 substitutions, D71G in GRIN3A and R727H in GRIN3B, caused differences in the functional assignments of these genes between humans and other apes. Conclusion We conclude that the GluR genes did not undergo drastic changes such as accelerated evolution in the human lineage after the divergence of chimpanzees. However, there remains a possibility that two human-specific "fixed" amino acid substitutions, D71G in GRIN3A and R727H in GRIN3B, are related to human-specific brain function.