The Chara Genome: Secondary Complexity and Implications for Plant Terrestrialization.

Land plants evolved from charophytic algae, among which Charophyceae possess the most complex body plans. We present the genome of Chara braunii; comparison of the genome to those of land plants identified evolutionary novelties for plant terrestrialization and land plant heritage genes. C. braunii employs unique xylan synthases for cell wall biosynthesis, a phragmoplast (cell separation) mechanism similar to that of land plants, and many phytohormones. C. braunii plastids are controlled via land-plant-like retrograde signaling, and transcriptional regulation is more elaborate than in other algae. The morphological complexity of this organism may result from expanded gene families, with three cases of particular note: genes effecting tolerance to reactive oxygen species (ROS), LysM receptor-like kinases, and transcription factors (TFs). Transcriptomic analysis of sexual reproductive structures reveals intricate control by TFs, activity of the ROS gene network, and the ancestral use of plant-like storage and stress protection proteins in the zygote.

An improved method for the highly specific detection of transcription start sites.

Precise detection of the transcriptional start site (TSS) is a key for characterizing transcriptional regulation of genes and for annotation of newly sequenced genomes. Here, we describe the development of an improved method, designated 'TSS-seq2.' This method is an iterative improvement of TSS-seq, a previously published enzymatic cap-structure conversion method to detect TSSs in base sequences. By modifying the original procedure, including by introducing split ligation at the key cap-selection step, the yield and the accuracy of the reaction has been substantially improved. For example, TSS-seq2 can be conducted using as little as 5 ng of total RNA with an overall accuracy of 96%; this yield a less-biased and more precise detection of TSS. We then applied TSS-seq2 for TSS analysis of four plant species that had not yet been analyzed by any previous TSS method.

MCM10 compensates for Myc-induced DNA replication stress in breast cancer stem-like cells.

Cancer stem-like cells (CSCs) induce drug resistance and recurrence of tumors when they experience DNA replication stress. However, the mechanisms underlying DNA replication stress in CSCs and its compensation remain unclear. Here, we demonstrate that upregulated c-Myc expression induces stronger DNA replication stress in patient-derived breast CSCs than in differentiated cancer cells. Our results suggest critical roles for mini-chromosome maintenance protein 10 (MCM10), a firing (activating) factor of DNA replication origins, to compensate for DNA replication stress in CSCs. MCM10 expression is upregulated in CSCs and is maintained by c-Myc. c-Myc-dependent collisions between RNA transcription and DNA replication machinery may occur in nuclei, thereby causing DNA replication stress. MCM10 may activate dormant replication origins close to these collisions to ensure the progression of replication. Moreover, patient-derived breast CSCs were found to be dependent on MCM10 for their maintenance, even after enrichment for CSCs that were resistant to paclitaxel, the standard chemotherapeutic agent. Further, MCM10 depletion decreased the growth of cancer cells, but not of normal cells. Therefore, MCM10 may robustly compensate for DNA replication stress and facilitate genome duplication in cancer cells in the S-phase, which is more pronounced in CSCs. Overall, we provide a preclinical rationale to target the c-Myc-MCM10 axis for preventing drug resistance and recurrence of tumors.

Monoallelic, antisense and total RNA transcription in an in vitro neural differentiation system based on F1 hybrid mice.

We developed an in vitro system to differentiate embryonic stem cells (ESCs) derived from reciprocally crossed F1 hybrid mice into neurons, and used it to investigate poly(A)+ and total RNA transcription at different stages of cell differentiation. By comparing expression profiles of transcripts assembled from 20 RNA sequencing datasets [2 alleles×(2 cell lines×4 time-points+2 mouse brains)], the relative influence of strain, cell and parent specificities to overall expression could be assessed. Divergent expression profiles of ESCs converged tightly at neural progenitor stage. Patterns of temporal variation of monoallelically expressed transcripts and antisense transcripts were quantified. Comparison of sense and antisense transcript pairs within the poly(A)+ sample, within the total RNA sample, and across poly(A)+ and total RNA samples revealed distinct rates of pairs showing anti-correlated expression variation. Unique patterns of sharing of poly(A)+ and poly(A)- transcription were identified in distinct RNA species. Regulation and functionality of monoallelic expression, antisense transcripts and poly(A)- transcription remain elusive. We demonstrated the effectiveness of our approach to capture these transcriptional activities, and provided new resources to elucidate the mammalian developmental transcriptome.

p53-inducible DPYSL4 associates with mitochondrial supercomplexes and regulates energy metabolism in adipocytes and cancer cells.

The tumor suppressor p53 regulates multiple cellular functions, including energy metabolism. Metabolic deregulation is implicated in the pathogenesis of some cancers and in metabolic disorders and may result from the inactivation of p53 functions. Using RNA sequencing and ChIP sequencing of cancer cells and preadipocytes, we demonstrate that p53 modulates several metabolic processes via the transactivation of energy metabolism genes including dihydropyrimidinase-like 4 (DPYSL4). DPYSL4 is a member of the collapsin response mediator protein family, which is involved in cancer invasion and progression. Intriguingly, DPYSL4 overexpression in cancer cells and preadipocytes up-regulated ATP production and oxygen consumption, while DPYSL4 knockdown using siRNA or CRISPR/Cas9 down-regulated energy production. Furthermore, DPYSL4 was associated with mitochondrial supercomplexes, and deletion of its dihydropyrimidinase-like domain abolished its association and its ability to stimulate ATP production and suppress the cancer cell invasion. Mouse-xenograft and lung-metastasis models indicated that DPYSL4 expression compromised tumor growth and metastasis in vivo. Consistently, database analyses demonstrated that low DPYSL4 expression was significantly associated with poor survival of breast and ovarian cancers in accordance with its reduced expression in certain types of cancer tissues. Moreover, immunohistochemical analysis using the adipose tissue of obese patients revealed that DPYSL4 expression was positively correlated with INFg and body mass index in accordance with p53 activation. Together, these results suggest that DPYSL4 plays a key role in the tumor-suppressor function of p53 by regulating oxidative phosphorylation and the cellular energy supply via its association with mitochondrial supercomplexes, possibly linking to the pathophysiology of both cancer and obesity.

A single female-specific piRNA is the primary determiner of sex in the silkworm.

The silkworm Bombyx mori uses a WZ sex determination system that is analogous to the one found in birds and some reptiles. In this system, males have two Z sex chromosomes, whereas females have Z and W sex chromosomes. The silkworm W chromosome has a dominant role in female determination, suggesting the existence of a dominant feminizing gene in this chromosome. However, the W chromosome is almost fully occupied by transposable element sequences, and no functional protein-coding gene has been identified so far. Female-enriched PIWI-interacting RNAs (piRNAs) are the only known transcripts that are produced from the sex-determining region of the W chromosome, but the function(s) of these piRNAs are unknown. Here we show that a W-chromosome-derived, female-specific piRNA is the feminizing factor of B. mori. This piRNA is produced from a piRNA precursor which we named Fem. Fem sequences were arranged in tandem in the sex-determining region of the W chromosome. Inhibition of Fem-derived piRNA-mediated signalling in female embryos led to the production of the male-specific splice variants of B. mori doublesex (Bmdsx), a gene which acts at the downstream end of the sex differentiation cascade. A target gene of Fem-derived piRNA was identified on the Z chromosome of B. mori. This gene, which we named Masc, encoded a CCCH-type zinc finger protein. We show that the silencing of Masc messenger RNA by Fem piRNA is required for the production of female-specific isoforms of Bmdsx in female embryos, and that Masc protein controls both dosage compensation and masculinization in male embryos. Our study characterizes a single small RNA that is responsible for primary sex determination in the WZ sex determination system.

DBTSS/DBKERO for integrated analysis of transcriptional regulation.

DBTSS (Database of Transcriptional Start Sites)/DBKERO (Database of Kashiwa Encyclopedia for human genome mutations in Regulatory regions and their Omics contexts) is the database originally initiated with the information of transcriptional start sites and their upstream transcriptional regulatory regions. In recent years, we updated the database to assist users to elucidate biological relevance of the human genome variations or somatic mutations in cancers which may affect the transcriptional regulation. In this update, we facilitate interpretations of disease associated genomic variation, using the Japanese population as a model case. We enriched the genomic variation dataset consisting of the 13,368 individuals collected for various genome-wide association studies and the reference epigenome information in the surrounding regions using a total of 455 epigenome datasets (four tissue types from 67 healthy individuals) collected for the International Human Epigenome Consortium (IHEC). The data directly obtained from the clinical samples was associated with that obtained from various model systems, such as the drug perturbation datasets using cultured cancer cells. Furthermore, we incorporated the results obtained using the newly developed analytical methods, Nanopore/10x Genomics long-read sequencing of the human genome and single cell analyses. The database is made publicly accessible at the URL (http://dbtss.hgc.jp/).

Artificial "ping-pong" cascade of PIWI-interacting RNA in silkworm cells.

PIWI-interacting RNAs (piRNAs) play essential roles in the defense system against selfish elements in animal germline cells by cooperating with PIWI proteins. A subset of piRNAs is predicted to be generated via the "ping-pong" cascade, which is mainly controlled by two different PIWI proteins. Here we established a cell-based artificial piRNA production system using a silkworm ovarian cultured cell line that is believed to possess a complete piRNA pathway. In addition, we took advantage of a unique silkworm sex-determining one-to-one ping-pong piRNA pair, which enabled us to precisely monitor the behavior of individual artificial piRNAs. With this novel strategy, we successfully generated artificial piRNAs against endogenous protein-coding genes via the expected back-and-forth traveling mechanism. Furthermore, we detected "primary" piRNAs from the upstream region of the artificial "ping-pong" site in the endogenous gene. This artificial piRNA production system experimentally confirms the existence of the "ping-pong" cascade of piRNAs. Also, this system will enable us to identify the factors involved in both, or each, of the "ping" and "pong" cascades and the sequence features that are required for efficient piRNA production.

Effects of androgen and oestrogen on IGF pathways controlling phallus growth.

The development of the mammalian phallus involves hormone-dependent mesenchymal-epithelial signalling mechanisms that contribute to urethral closure and regulation of phallus elongation and growth. In marsupials, most differentiation and growth of the phallus occurs post-natally, making them amenable to direct hormone treatment. Expression of IGFs, FGFs, EFNB2, MAFB, DLX5 and AP-1 mRNAs in the phallus at day 50 post-partum (pp) were altered after treatment of tammar wallaby young from day 20 to 40 pp with androgen, oestrogen or after castration at day 25 pp. However, the most interesting changes occurred in the IGF pathway genes. Androgen treatment upregulated IGF1 in female phalluses and oestrogen treatment upregulated IGF1 in male phalluses, but it was downregulated by castration. IGFBP3 was higher in female phalluses and downregulated by androgen. IGF1 expression was higher in all untreated male than in female phalluses from day 50 to 150 pp, but IGFBP3 had the reverse pattern. At day 90 pp, when urethral closure in males is progressing and male phallus growth is accelerating. IGF1 and PCNA protein were only detected in the male urorectal septum, suggesting for the first time that closure and elongation may involve IGF1 activation of cell proliferation specifically in male phalluses. These effects of sex steroids on gene expression and on the IGF1 signalling pathway in particular, suggest that the developing phallus may be especially susceptible to perturbation by exogenous hormones.

Restriction glycosylases: involvement of endonuclease activities in the restriction process.

All restriction enzymes examined are phosphodiesterases generating 3΄-OH and 5΄-P ends, but one restriction enzyme (restriction glycosylase) excises unmethylated bases from its recognition sequence. Whether its restriction activity involves endonucleolytic cleavage remains unclear. One report on this enzyme, R.PabI from a hyperthermophile, ascribed the breakage to high temperature while another showed its weak AP lyase activity generates atypical ends. Here, we addressed this issue in mesophiles. We purified R.PabI homologs from Campylobacter coli (R.CcoLI) and Helicobacter pylori (R.HpyAXII) and demonstrated their DNA cleavage, DNA glycosylase and AP lyase activities in vitro at 37°C. The AP lyase activity is more coupled with glycosylase activity in R.CcoLI than in R.PabI. R.CcoLI/R.PabI expression caused restriction of incoming bacteriophage/plasmid DNA and endogenous chromosomal DNA within Escherichia coli at 37°C. The R.PabI-mediated restriction was promoted by AP endonuclease action in vivo or in vitro. These results reveal the role of endonucleolytic DNA cleavage in restriction and yet point to diversity among the endonucleases. The cleaved ends are difficult to repair in vivo, which may indicate their biological significance. These results support generalization of the concept of restriction­modification system to the concept of self-recognizing epigenetic system, which combines any epigenetic labeling and any DNA damaging.

Identification of physical interactions between genomic regions by enChIP-Seq.

Physical interactions between genomic regions play critical roles in the regulation of genome functions, including gene expression. Here, we show the feasibility of using engineered DNA-binding molecule-mediated chromatin immunoprecipitation (enChIP) in combination with next-generation sequencing (NGS) (enChIP-Seq) to detect such interactions. In enChIP-Seq, the target genomic region is captured by an engineered DNA-binding complex, such as a clustered regularly interspaced short palindromic repeats (CRISPR) system consisting of a catalytically inactive form of Cas9 and a single guide RNA. Subsequently, the genomic regions that physically interact with the target genomic region in the captured complex are sequenced by NGS. Using enChIP-Seq, we found that the 5'HS5 locus, which is involved in the regulation of globin genes expression at the ß-globin locus, interacts with multiple genomic regions upon erythroid differentiation in the human erythroleukemia cell line K562. Genes near the genomic regions inducibly associated with the 5'HS5 locus were transcriptionally up-regulated in the differentiated state, suggesting the existence of a coordinated transcription mechanism mediated by physical interactions between these loci. Thus, enChIP-Seq might be a potentially useful tool for detecting physical interactions between genomic regions in a nonbiased manner, which would facilitate elucidation of the molecular mechanisms underlying regulation of genome functions.

Hormone-responsive genes in the SHH and WNT/ß-catenin signaling pathways influence urethral closure and phallus growth.

Environmental endocrine disruptors (EEDs) that affect androgen or estrogen activity may disrupt gene regulation during phallus development to cause hypospadias or a masculinized clitoris. We treated developing male tammar wallabies with estrogen and females with androgen from day 20-40 postpartum (pp) during the androgen imprinting window of sensitivity. Estrogen inhibited phallus elongation but had no effect on urethral closure and did not significantly depress testicular androgen synthesis. Androgen treatment in females did not promote phallus elongation but initiated urethral closure. Phalluses were collected for transcriptome sequencing at day 50 pp when they first become sexually dimorphic to examine changes in two signaling pathways, sonic hedgehog (SHH) and wingless-type MMTV integration site family (WNT)/ß-catenin. SHH mRNA and ß-catenin were predominantly expressed in the urethral epithelium in the tammar phallus, as in eutherian mammals. Estrogen treatment and castration of males induced an upregulation of SHH, while androgen treatment downregulated SHH. These effects appear to be direct since we detected putative estrogen receptor α (ERα) and androgen receptor (AR) binding sites near SHH. WNT5A, like SHH, was downregulated by androgen, while WNT4 was upregulated in female phalluses after androgen treatment. After estrogen treatment, WIF1 and WNT7A were both downregulated in male phalluses. After castration, WNT9A was upregulated. These results suggest that SHH and WNT pathways are regulated by both estrogen and androgen to direct the proliferation and elongation of the phallus during differentiation. Their response to exogenous hormones makes these genes potential targets of EEDs in the etiology of abnormal phallus development including hypospadias.

No novel, high penetrant gene might remain to be found in Japanese patients with unknown MODY.

MODY 5 and 6 have been shown to be low-penetrant MODYs. As the genetic background of unknown MODY is assumed to be similar, a new analytical strategy is applied here to elucidate genetic predispositions to unknown MODY. We examined to find whether there are major MODY gene loci remaining to be identified using SNP linkage analysis in Japanese. Whole-exome sequencing was performed with seven families with typical MODY. Candidates for novel MODY genes were examined combined with in silico network analysis. Some peaks were found only in either parametric or non-parametric analysis; however, none of these peaks showed a LOD score greater than 3.7, which is approved to be the significance threshold of evidence for linkage. Exome sequencing revealed that three mutated genes were common among 3 families and 42 mutated genes were common in two families. Only one of these genes, MYO5A, having rare amino acid mutations p.R849Q and p.V1601G, was involved in the biological network of known MODY genes through the intermediary of the INS. Although only one promising candidate gene, MYO5A, was identified, no novel, high penetrant MODY genes might remain to be found in Japanese MODY.

Local potentiation of stress-responsive genes by upstream noncoding transcription.

It has been postulated that a myriad of long noncoding RNAs (lncRNAs) contribute to gene regulation. In fission yeast, glucose starvation triggers lncRNA transcription across promoter regions of stress-responsive genes including fbp1 (fructose-1,6-bisphosphatase1). At the fbp1 promoter, this transcription promotes chromatin remodeling and fbp1 mRNA expression. Here, we demonstrate that such upstream noncoding transcription facilitates promoter association of the stress-responsive transcriptional activator Atf1 at the sites of transcription, leading to activation of the downstream stress genes. Genome-wide analyses revealed that ∼50 Atf1-binding sites show marked decrease in Atf1 occupancy when cells are treated with a transcription inhibitor. Most of these transcription-enhanced Atf1-binding sites are associated with stress-dependent induction of the adjacent mRNAs or lncRNAs, as observed in fbp1 These Atf1-binding sites exhibit low Atf1 occupancy and high histone density in glucose-rich conditions, and undergo dramatic changes in chromatin status after glucose depletion: enhanced Atf1 binding, histone eviction, and histone H3 acetylation. We also found that upstream transcripts bind to the Groucho-Tup1 type transcriptional corepressors Tup11 and Tup12, and locally antagonize their repressive functions on Atf1 binding. These results reveal a new mechanism in which upstream noncoding transcription locally magnifies the specific activation of stress-inducible genes via counteraction of corepressors.

A Genome Scan for Genes Underlying Microgeographic-Scale Local Adaptation in a Wild Arabidopsis Species.

Adaptive divergence at the microgeographic scale has been generally disregarded because high gene flow is expected to disrupt local adaptation. Yet, growing number of studies reporting adaptive divergence at a small spatial scale highlight the importance of this process in evolutionary biology. To investigate the genetic basis of microgeographic local adaptation, we conducted a genome-wide scan among sets of continuously distributed populations of Arabidopsis halleri subsp. gemmifera that show altitudinal phenotypic divergence despite gene flow. Genomic comparisons were independently conducted in two distinct mountains where similar highland ecotypes are observed, presumably as a result of convergent evolution. Here, we established a de novo reference genome and employed an individual-based resequencing for a total of 56 individuals. Among 527,225 reliable SNP loci, we focused on those showing a unidirectional allele frequency shift across altitudes. Statistical tests on the screened genes showed that our microgeographic population genomic approach successfully retrieve genes with functional annotations that are in line with the known phenotypic and environmental differences between altitudes. Furthermore, comparison between the two distinct mountains enabled us to screen out those genes that are neutral or adaptive only in either mountain, and identify the genes involved in the convergent evolution. Our study demonstrates that the genomic comparison among a set of genetically connected populations, instead of the commonly-performed comparison between two isolated populations, can also offer an effective screening for the genetic basis of local adaptation.

Genomic Landscape of Esophageal Squamous Cell Carcinoma in a Japanese Population.

BACKGROUND & AIMS: Esophageal squamous cell carcinoma (ESCC) is the predominant form of esophageal cancer in Japan. Smoking and drinking alcohol are environmental risk factors for ESCC, whereas single nucleotide polymorphisms in ADH1B and ALDH2, which increase harmful intermediates produced by drinking alcohol, are genetic risk factors. We conducted a large-scale genomic analysis of ESCCs from patients in Japan to determine the mutational landscape of this cancer. METHODS: We performed whole-exome sequence analysis of tumor and nontumor esophageal tissues collected from 144 patients with ESCC who underwent surgery at 5 hospitals in Japan. We also performed single-nucleotide polymorphism array-based copy number profile and germline genotype analyses of polymorphisms in ADH1B and ALDH2. Polymorphisms in CYP2A6, which increase harmful effects of smoking, were analyzed. Functions of TET2 mutants were evaluated in KYSE410 and HEK293FT cells. RESULTS: A high proportion of mutations in the 144 tumor samples were C to T substitution in CpG dinucleotides (called the CpG signature) and C to G/T substitutions with a flanking 5' thymine (called the APOBEC signature). Based on mutational signatures, patients were assigned to 3 groups, which associated with environmental (drinking and smoking) and genetic (polymorphisms in ALDH2 and CYP2A6) factors. Many tumors contained mutations in genes that regulate the cell cycle (TP53, CCND1, CDKN2A, FBXW7); epigenetic processes (MLL2, EP300, CREBBP, TET2); and the NOTCH (NOTCH1, NOTCH3), WNT (FAT1, YAP1, AJUBA) and receptor-tyrosine kinase-phosphoinositide 3-kinase signaling pathways (PIK3CA, EGFR, ERBB2). Mutations in EP300 and TET2 correlated with shorter survival times, and mutations in ZNF750 associated with an increased number of mutations of the APOBEC signature. Expression of mutant forms of TET2 did not increase cellular levels of 5-hydroxymethylcytosine in HEK293FT cells, whereas knockdown of TET2 increased the invasive activity of KYSE410 ESCC cells. Computational analyses associated the mutations in NFE2L2 we identified with transcriptional activation of its target genes. CONCLUSIONS: We associated environmental and genetic factors with base substitution patterns of somatic mutations and provide a registry of genes and pathways that are disrupted in ESCCs. These findings might be used to design specific treatments for patients with esophageal squamous cancers.

Large hypomethylated domains serve as strong repressive machinery for key developmental genes in vertebrates.

DNA methylation is a fundamental epigenetic modification in vertebrate genomes and a small fraction of genomic regions is hypomethylated. Previous studies have implicated hypomethylated regions in gene regulation, but their functions in vertebrate development remain elusive. To address this issue, we generated epigenomic profiles that include base-resolution DNA methylomes and histone modification maps from both pluripotent cells and mature organs of medaka fish and compared the profiles with those of human ES cells. We found that a subset of hypomethylated domains harbor H3K27me3 (K27HMDs) and their size positively correlates with the accumulation of H3K27me3. Large K27HMDs are conserved between medaka and human pluripotent cells and predominantly contain promoters of developmental transcription factor genes. These key genes were found to be under strong transcriptional repression, when compared with other developmental genes with smaller K27HMDs. Furthermore, human-specific K27HMDs show an enrichment of neuronal activity-related genes, which suggests a distinct regulation of these genes in medaka and human. In mature organs, some of the large HMDs become shortened by elevated DNA methylation and associate with sustained gene expression. This study highlights the significance of domain size in epigenetic gene regulation. We propose that large K27HMDs play a crucial role in pluripotent cells by strictly repressing key developmental genes, whereas their shortening consolidates long-term gene expression in adult differentiated cells.

The Endosymbiotic Bacterium Wolbachia Selectively Kills Male Hosts by Targeting the Masculinizing Gene.

Pathogens are known to manipulate the reproduction and development of their hosts for their own benefit. Wolbachia is an endosymbiotic bacterium that infects a wide range of insect species. Wolbachia is known as an example of a parasite that manipulates the sex of its host's progeny. Infection of Ostrinia moths by Wolbachia causes the production of all-female progeny, however, the mechanism of how Wolbachia accomplishes this male-specific killing is unknown. Here we show for the first time that Wolbachia targets the host masculinizing gene of Ostrinia to accomplish male-killing. We found that Wolbachia-infected O. furnacalis embryos do not express the male-specific splice variant of doublesex, a gene which acts at the downstream end of the sex differentiation cascade, throughout embryonic development. Transcriptome analysis revealed that Wolbachia infection markedly reduces the mRNA level of Masc, a gene that encodes a protein required for both masculinization and dosage compensation in the silkworm Bombyx mori. Detailed bioinformatic analysis also elucidated that dosage compensation of Z-linked genes fails in Wolbachia-infected O. furnacalis embryos, a phenomenon that is extremely similar to that observed in Masc mRNA-depleted male embryos of B. mori. Finally, injection of in vitro transcribed Masc cRNA into Wolbachia-infected embryos rescued male progeny. Our results show that Wolbachia-induced male-killing is caused by a failure of dosage compensation via repression of the host masculinizing gene. Our study also shows a novel strategy by which a pathogen hijacks the host sex determination cascade.

Maternal variant in the upstream of FOXP3 gene on the X chromosome is associated with recurrent infertility in Japanese Black cattle.

BACKGROUND: Repeat breeding, which is defined as cattle failure to conceive after three or more inseminations in the absence of clinical abnormalities, is a substantial problem in cattle breeding. To identify maternal genetic variants of repeat breeding in Japanese Black cattle, we selected 29 repeat-breeding heifers that failed to conceive following embryo transfer (ET) and conducted a genome-wide association study (GWAS) using the traits. RESULTS: We found that a single-nucleotide polymorphism (SNP; g.92,377,635A > G) in the upstream region of the FOXP3 gene on the X chromosome was highly associated with repeat breeding and failure to conceive following ET (P = 1.51 × 10-14). FOXP3 is a master gene for differentiation of regulatory T (Treg) cells that function in pregnancy maintenance. Reporter assay results revealed that the activity of the FOXP3 promoter was lower in reporter constructs with the risk-allele than in those with the non-risk-allele by approximately 0.68 fold. These findings suggest that the variant in the upstream region of FOXP3 with the risk-allele decreased FOXP3 transcription, which in turn, could reduce the number of maternal Treg cells and lead to infertility. The frequency of the risk-allele in repeat-breeding heifers is more than that in cows, suggesting that the risk-allele could be associated with infertility in repeat-breeding heifers. CONCLUSIONS: This GWAS identified a maternal variant in the upstream region of FOXP3 that was associated with infertility in repeat-breeding Japanese Black cattle that failed to conceive using ET. The variant affected the level of FOXP3 mRNA expression. Thus, the results suggest that the risk-allele could serve as a useful marker to reduce and eliminate animals with inferior fertility in Japanese Black cattle.

Frequent pathway mutations of splicing machinery in myelodysplasia.

Myelodysplastic syndromes and related disorders (myelodysplasia) are a heterogeneous group of myeloid neoplasms showing deregulated blood cell production with evidence of myeloid dysplasia and a predisposition to acute myeloid leukaemia, whose pathogenesis is only incompletely understood. Here we report whole-exome sequencing of 29 myelodysplasia specimens, which unexpectedly revealed novel pathway mutations involving multiple components of the RNA splicing machinery, including U2AF35, ZRSR2, SRSF2 and SF3B1. In a large series analysis, these splicing pathway mutations were frequent (∼45 to ∼85%) in, and highly specific to, myeloid neoplasms showing features of myelodysplasia. Conspicuously, most of the mutations, which occurred in a mutually exclusive manner, affected genes involved in the 3'-splice site recognition during pre-mRNA processing, inducing abnormal RNA splicing and compromised haematopoiesis. Our results provide the first evidence indicating that genetic alterations of the major splicing components could be involved in human pathogenesis, also implicating a novel therapeutic possibility for myelodysplasia.