Brief report: high-throughput sequencing of IL23R reveals a low-frequency, nonsynonymous single-nucleotide polymorphism that is associated with ankylosing spondylitis in a Han Chinese population.

OBJECTIVE: Ankylosing spondylitis (AS) is a highly heritable common inflammatory arthritis that targets the spine and sacroiliac joints of the pelvis, causing pain and stiffness and leading eventually to joint fusion. Although previous studies have shown a strong association of IL23R with AS in white Europeans, similar studies in East Asian populations have shown no association with common variants of IL23R, suggesting either that IL23R variants have no role or that rare genetic variants contribute. The present study was undertaken to screen IL23R to identify rare variants associated with AS in Han Chinese. METHODS: A 170-kb region containing IL23R and its flanking regions was sequenced in 50 patients with AS and 50 ethnically matched healthy control subjects from a Han Chinese population. In addition, the 30-kb region of peak association in white Europeans was sequenced in 650 patients with AS and 1,300 healthy controls. Validation genotyping was undertaken in 846 patients with AS and 1,308 healthy controls. RESULTS: We identified 1,047 variants, of which 729 were not found in the dbSNP genomic build 130. Several potentially functional rare variants in IL23R were identified, including one nonsynonomous single-nucleotide polymorphism (nsSNP), Gly(149) Arg (position 67421184 GA on chromosome 1). Validation genotyping showed that the Gly(149) Arg variant was associated with AS (odds ratio 0.61, P = 0.0054). CONCLUSION: This is the first study to implicate rare IL23R variants in the pathogenesis of AS. The results identified a low-frequency nsSNP with predicted loss-of-function effects that was protectively associated with AS in Han Chinese, suggesting that decreased function of the interleukin-23 (IL-23) receptor protects against AS. These findings further support the notion that IL-23 signaling has an important role in the pathogenesis of AS.

MicroRNAs-140-5p/140-3p modulate Leydig cell numbers in the developing mouse testis.

MicroRNAs (miRNAs) have been shown to play key regulatory roles in a range of biological processes, including cell differentiation and development. To identify miRNAs that participate in gonad differentiation, a fundamental and tightly regulated developmental process, we examined miRNA expression profiles at the time of sex determination and during the early fetal differentiation of mouse testes and ovaries using high-throughput sequencing. We identified several miRNAs that were expressed in a sexually dimorphic pattern, including several members of the let-7 family, miR-378, and miR-140-3p. We focused our analysis on the most highly expressed, sexually dimorphic miRNA, miR-140-3p, and found that both miR-140-3p and its more lowly expressed counterpart, the previously annotated guide strand, miR-140-5p, are testis enriched and expressed in testis cords. Analysis of the miR-140-5p/miR-140-3p-null mouse revealed a significant increase in the number of Leydig cells in the developing XY gonad, strongly suggesting an important role for miR-140-5p/miR-140-3p in testis differentiation in mouse.

Resequencing and fine-mapping of the chromosome 12q13-14 locus associated with multiple sclerosis refines the number of implicated genes.

Multiple sclerosis (MS) is a common chronic inflammatory disease of the central nervous system. Susceptibility to the disease is affected by both environmental and genetic factors. Genetic factors include haplotypes in the histocompatibility complex (MHC) and over 50 non-MHC loci reported by genome-wide association studies. Amongst these, we previously reported polymorphisms in chromosome 12q13-14 with a protective effect in individuals of European descent. This locus spans 288 kb and contains 17 genes, including several candidate genes which have potentially significant pathogenic and therapeutic implications. In this study, we aimed to fine-map this locus. We have implemented a two-phase study: a variant discovery phase where we have used next-generation sequencing and two target-enrichment strategies [long-range polymerase chain reaction (PCR) and Nimblegen's solution phase hybridization capture] in pools of 25 samples; and a genotyping phase where we genotyped 712 variants in 3577 healthy controls and 3269 MS patients. This study confirmed the association (rs2069502, P = 9.9 × 10(-11), OR = 0.787) and narrowed down the locus of association to an 86.5 kb region. Although the study was unable to pinpoint the key-associated variant, we have identified a 42 (genotyped and imputed) single-nucleotide polymorphism haplotype block likely to harbour the causal variant. No evidence of association at previously reported low-frequency variants in CYP27B1 was observed. As part of the study we compared variant discovery performance using two target-enrichment strategies. We concluded that our pools enriched with Nimblegen's solution phase hybridization capture had better sensitivity to detect true variants than the pools enriched with long-range PCR, whilst specificity was better in the long-range PCR-enriched pools compared with solution phase hybridization capture enriched pools; this result has important implications for the design of future fine-mapping studies.

Novel roles for KLF1 in erythropoiesis revealed by mRNA-seq.

KLF1 (formerly known as EKLF) regulates the development of erythroid cells from bi-potent progenitor cells via the transcriptional activation of a diverse set of genes. Mice lacking Klf1 die in utero prior to E15 from severe anemia due to the inadequate expression of genes controlling hemoglobin production, cell membrane and cytoskeletal integrity, and the cell cycle. We have recently described the full repertoire of KLF1 binding sites in vivo by performing KLF1 ChIP-seq in primary erythroid tissue (E14.5 fetal liver). Here we describe the KLF1-dependent erythroid transcriptome by comparing mRNA-seq from Klf1(+/+) and Klf1(-/-) erythroid tissue. This has revealed novel target genes not previously obtainable by traditional microarray technology, and provided novel insights into the function of KLF1 as a transcriptional activator. We define a cis-regulatory module bound by KLF1, GATA1, TAL1, and EP300 that coordinates a core set of erythroid genes. We also describe a novel set of erythroid-specific promoters that drive high-level expression of otherwise ubiquitously expressed genes in erythroid cells. Our study has identified two novel lncRNAs that are dynamically expressed during erythroid differentiation, and discovered a role for KLF1 in directing apoptotic gene expression to drive the terminal stages of erythroid maturation.

Multicentric carpotarsal osteolysis is caused by mutations clustering in the amino-terminal transcriptional activation domain of MAFB.

Multicentric carpotarsal osteolysis (MCTO) is a rare skeletal dysplasia characterized by aggressive osteolysis, particularly affecting the carpal and tarsal bones, and is frequently associated with progressive renal failure. Using exome capture and next-generation sequencing in five unrelated simplex cases of MCTO, we identified previously unreported missense mutations clustering within a 51 base pair region of the single exon of MAFB, validated by Sanger sequencing. A further six unrelated simplex cases with MCTO were also heterozygous for previously unreported mutations within this same region, as were affected members of two families with autosomal-dominant MCTO. MAFB encodes a transcription factor that negatively regulates RANKL-induced osteoclastogenesis and is essential for normal renal development. Identification of this gene paves the way for development of novel therapeutic approaches for this crippling disease and provides insight into normal bone and kidney development.

Whole-exome re-sequencing in a family quartet identifies POP1 mutations as the cause of a novel skeletal dysplasia.

Recent advances in DNA sequencing have enabled mapping of genes for monogenic traits in families with small pedigrees and even in unrelated cases. We report the identification of disease-causing mutations in a rare, severe, skeletal dysplasia, studying a family of two healthy unrelated parents and two affected children using whole-exome sequencing. The two affected daughters have clinical and radiographic features suggestive of anauxetic dysplasia (OMIM 607095), a rare form of dwarfism caused by mutations of RMRP. However, mutations of RMRP were excluded in this family by direct sequencing. Our studies identified two novel compound heterozygous loss-of-function mutations in POP1, which encodes a core component of the RNase mitochondrial RNA processing (RNase MRP) complex that directly interacts with the RMRP RNA domains that are affected in anauxetic dysplasia. We demonstrate that these mutations impair the integrity and activity of this complex and that they impair cell proliferation, providing likely molecular and cellular mechanisms by which POP1 mutations cause this severe skeletal dysplasia.

Integrated genome-wide chromatin occupancy and expression analyses identify key myeloid pro-differentiation transcription factors repressed by Myb.

To gain insight into the mechanisms by which the Myb transcription factor controls normal hematopoiesis and particularly, how it contributes to leukemogenesis, we mapped the genome-wide occupancy of Myb by chromatin immunoprecipitation followed by massively parallel sequencing (ChIP-Seq) in ERMYB myeloid progenitor cells. By integrating the genome occupancy data with whole genome expression profiling data, we identified a Myb-regulated transcriptional program. Gene signatures for leukemia stem cells, normal hematopoietic stem/progenitor cells and myeloid development were overrepresented in 2368 Myb regulated genes. Of these, Myb bound directly near or within 793 genes. Myb directly activates some genes known critical in maintaining hematopoietic stem cells, such as Gfi1 and Cited2. Importantly, we also show that, despite being usually considered as a transactivator, Myb also functions to repress approximately half of its direct targets, including several key regulators of myeloid differentiation, such as Sfpi1 (also known as Pu.1), Runx1, Junb and Cebpb. Furthermore, our results demonstrate that interaction with p300, an established coactivator for Myb, is unexpectedly required for Myb-mediated transcriptional repression. We propose that the repression of the above mentioned key pro-differentiation factors may contribute essentially to Myb's ability to suppress differentiation and promote self-renewal, thus maintaining progenitor cells in an undifferentiated state and promoting leukemic transformation.

Expression of distinct RNAs from 3' untranslated regions.

The 3' untranslated regions (3'UTRs) of eukaryotic genes regulate mRNA stability, localization and translation. Here, we present evidence that large numbers of 3'UTRs in human, mouse and fly are also expressed separately from the associated protein-coding sequences to which they are normally linked, likely by post-transcriptional cleavage. Analysis of CAGE (capped analysis of gene expression), SAGE (serial analysis of gene expression) and cDNA libraries, as well as microarray expression profiles, demonstrate that the independent expression of 3'UTRs is a regulated and conserved genome-wide phenomenon. We characterize the expression of several 3'UTR-derived RNAs (uaRNAs) in detail in mouse embryos, showing by in situ hybridization that these transcripts are expressed in a cell- and subcellular-specific manner. Our results suggest that 3'UTR sequences can function not only in cis to regulate protein expression, but also intrinsically and independently in trans, likely as noncoding RNAs, a conclusion supported by a number of previous genetic studies. Our findings suggest novel functions for 3'UTRs, as well as caution in the use of 3'UTR sequence probes to analyze gene expression.

Induction of pluripotency in human endothelial cells resets epigenetic profile on genome scale.

Reprogramming of a limited number of human cell types has been achieved through ectopic expression of four transcription factors to yield induced pluripotent stem (iPS) cells that closely resemble human embryonic stem cells (ESCs). Here, we determined functional and epigenetic properties of iPS cells generated from human umbilical vein endothelial cells (HUVEC) by conventional method of direct reprogramming. Retroviral overexpression of four transcription factors resets HUVEC to the pluripotency. Human endothelial cell-derived iPS (endo-iPS) cells were similar to human ESCs in morphology, gene expression, in vitro and in vivo differentiation capacity. Endo-iPS cells were efficiently differentiated in vitro into endothelial cells. Using genome-wide methylation profiling we show that promoter elements of endothelial specific genes were methylated following reprogramming whereas pluripotency-related gene promoters were hypomethylated similar to levels observed in ESCs. Genome-wide methylation analysis of CpG sites located in the functional regions of over than 14,000 genes indicated that human endo-iPS cells were highly similar to human ES cells, although differences in methylation levels of 46 genes were found. Overall CpG methylation of promoter regions in the pluripotent cells was higher than in somatic. We also show that during reprogramming female human endo-iPS cells exhibited reactivation of the somatically silenced X chromosome. Our findings demonstrate that iPS cells can be generated from human endothelial cells and reprogramming resets epigenetic status of endothelial cells to pluripotency.

Genome-wide association study of ankylosing spondylitis identifies non-MHC susceptibility loci.

To identify susceptibility loci for ankylosing spondylitis, we undertook a genome-wide association study in 2,053 unrelated ankylosing spondylitis cases among people of European descent and 5,140 ethnically matched controls, with replication in an independent cohort of 898 ankylosing spondylitis cases and 1,518 controls. Cases were genotyped with Illumina HumHap370 genotyping chips. In addition to strong association with the major histocompatibility complex (MHC; P < 10(-800)), we found association with SNPs in two gene deserts at 2p15 (rs10865331; combined P = 1.9 x 10(-19)) and 21q22 (rs2242944; P = 8.3 x 10(-20)), as well as in the genes ANTXR2 (rs4333130; P = 9.3 x 10(-8)) and IL1R2 (rs2310173; P = 4.8 x 10(-7)). We also replicated previously reported associations at IL23R (rs11209026; P = 9.1 x 10(-14)) and ERAP1 (rs27434; P = 5.3 x 10(-12)). This study reports four genetic loci associated with ankylosing spondylitis risk and identifies a major role for the interleukin (IL)-23 and IL-1 cytokine pathways in disease susceptibility.

Characterization of microRNAs encoded by the bovine herpesvirus 1 genome.

Bovine herpesvirus 1 (BoHV-1) is a ubiquitous and important pathogen of cattle worldwide. This study reports the identification of 10 microRNA (miRNA) genes, Bhv1-mir-B1-Bhv1-mir-B10, encoded by the BoHV-1 genome that were processed into 12 detectable mature miRNAs as determined by ultra-high throughput sequencing bioinformatics analyses of small RNA libraries and expression studies. We found that four of the miRNA genes were present as two copies in the BoHV-1 genome, resulting in a total of 14 miRNA encoding loci. Unique features of the BoHV-1 miRNAs include evidence of bidirectional transcription and a close association of two miRNA genes with the origin of replication, including one miRNA that is encoded within the origin of replication. The miRNA gene Bhv1-mir-B5 was encoded on the opposite DNA strand to the latency associated transcript, potentially giving rise to antisense transcripts originating from this locus. The association of herpesvirus miRNAs with latency appears to be a common feature in the alphaherpesviruses. Analyses of the BoHV-5 genome for putative miRNA gene orthologues identified a high degree of evolutionary conservation for nine of the BoHV-1 miRNA genes. The possible roles for BoHV-1 miRNAs in the regulation of known BoHV-1 transcription units and the genetics of the BoHV-1 genotypes are also discussed.

Repertoire of bovine miRNA and miRNA-like small regulatory RNAs expressed upon viral infection.

MicroRNA (miRNA) and other types of small regulatory RNAs play a crucial role in the regulation of gene expression in eukaryotes. Several distinct classes of small regulatory RNAs have been discovered in recent years. To extend the repertoire of small RNAs characterized in mammals and to examine relationship between host miRNA expression and viral infection we used Illumina's ultrahigh throughput sequencing approach. We sequenced three small RNA libraries prepared from cell line derived from the adult bovine kidney under normal conditions and upon infection of the cell line with Bovine herpesvirus 1. We used a bioinformatics approach to distinguish authentic mature miRNA sequences from other classes of small RNAs and short RNA fragments represented in the sequencing data. Using this approach we detected 219 out of 356 known bovine miRNAs and 115 respective miRNA* sequences. In addition we identified five new bovine orthologs of known mammalian miRNAs and discovered 268 new cow miRNAs many of which are not identifiable in other mammalian genomes and thus might be specific to the ruminant lineage. In addition we found seven new bovine mirtron candidates. We also discovered 10 small nucleolar RNA (snoRNA) loci that give rise to small RNA with possible miRNA-like function. Results presented in this study extend our knowledge of the biology and evolution of small regulatory RNAs in mammals and illuminate mechanisms of small RNA biogenesis and function. New miRNA sequences and the original sequencing data have been submitted to miRNA repository (miRBase) and NCBI GEO archive respectively. We envisage that these resources will facilitate functional annotation of the bovine genome and promote further functional and comparative genomics studies of small regulatory RNA in mammals.

Small RNAs derived from snoRNAs.

Small nucleolar RNAs (snoRNAs) guide RNA modification and are localized in nucleoli and Cajal bodies in eukaryotic cells. Components of the RNA silencing pathway associate with these structures, and two recent reports have revealed that a human and a protozoan snoRNA can be processed into miRNA-like RNAs. Here we show that small RNAs with evolutionary conservation of size and position are derived from the vast majority of snoRNA loci in animals (human, mouse, chicken, fruit fly), Arabidopsis, and fission yeast. In animals, sno-derived RNAs (sdRNAs) from H/ACA snoRNAs are predominantly 20-24 nucleotides (nt) in length and originate from the 3' end. Those derived from C/D snoRNAs show a bimodal size distribution at approximately 17-19 nt and >27 nt and predominantly originate from the 5' end. SdRNAs are associated with AGO7 in Arabidopsis and Ago1 in fission yeast with characteristic 5' nucleotide biases and show altered expression patterns in fly loquacious and Dicer-2 and mouse Dicer1 and Dgcr8 mutants. These findings indicate that there is interplay between the RNA silencing and snoRNA-mediated RNA processing systems, and that sdRNAs comprise a novel and ancient class of small RNAs in eukaryotes.

Tiny RNAs associated with transcription start sites in animals.

It has been reported that relatively short RNAs of heterogeneous sizes are derived from sequences near the promoters of eukaryotic genes. In conjunction with the FANTOM4 project, we have identified tiny RNAs with a modal length of 18 nt that map within -60 to +120 nt of transcription start sites (TSSs) in human, chicken and Drosophila. These transcription initiation RNAs (tiRNAs) are derived from sequences on the same strand as the TSS and are preferentially associated with G+C-rich promoters. The 5' ends of tiRNAs show peak density 10-30 nt downstream of TSSs, indicating that they are processed. tiRNAs are generally, although not exclusively, associated with highly expressed transcripts and sites of RNA polymerase II binding. We suggest that tiRNAs may be a general feature of transcription in metazoa and possibly all eukaryotes.

Analysis of the complement and molecular evolution of tRNA genes in cow.

BACKGROUND: Detailed information regarding the number and organization of transfer RNA (tRNA) genes at the genome level is becoming readily available with the increase of DNA sequencing of whole genomes. However the identification of functional tRNA genes is challenging for species that have large numbers of repetitive elements containing tRNA derived sequences, such as Bos taurus. Reliable identification and annotation of entire sets of tRNA genes allows the evolution of tRNA genes to be understood on a genomic scale. RESULTS: In this study, we explored the B. taurus genome using bioinformatics and comparative genomics approaches to catalogue and analyze cow tRNA genes. The initial analysis of the cow genome using tRNAscan-SE identified 31,868 putative tRNA genes and 189,183 pseudogenes, where 28,830 of the 31,868 predicted tRNA genes were classified as repetitive elements by the RepeatMasker program. We then used comparative genomics to further discriminate between functional tRNA genes and tRNA-derived sequences for the remaining set of 3,038 putative tRNA genes. For our analysis, we used the human, chimpanzee, mouse, rat, horse, dog, chicken and fugu genomes to predict that the number of active tRNA genes in cow lies in the vicinity of 439. Of this set, 150 tRNA genes were 100% identical in their sequences across all nine vertebrate genomes studied. Using clustering analyses, we identified a new tRNA-GlyCCC subfamily present in all analyzed mammalian genomes. We suggest that this subfamily originated from an ancestral tRNA-GlyGCC gene via a point mutation prior to the radiation of the mammalian lineages. Lastly, in a separate analysis we created phylogenetic profiles for each putative cow tRNA gene using a representative set of genomes to gain an overview of common evolutionary histories of tRNA genes. CONCLUSION: The use of a combination of bioinformatics and comparative genomics approaches has allowed the confident identification of a set of cow tRNA genes that will facilitate further studies in understanding the molecular evolution of cow tRNA genes.

A microRNA catalog of the developing chicken embryo identified by a deep sequencing approach.

MicroRNA (miRNA) and other types of small regulatory RNAs play a crucial role in the regulation of gene expression in eukaryotes. Several distinct classes of small regulatory RNAs have been discovered in recent years. To extend the repertoire of small regulatory RNAs characterized in chickens we used a deep sequencing approach developed by Solexa (now Illumina Inc.). We sequenced three small RNA libraries prepared from different developmental stages of the chicken embryo (days five, seven, and nine) to produce over 9.5 million short sequence reads. We developed a bioinformatics pipeline to distinguish authentic mature miRNA sequences from other classes of small RNAs and short RNA fragments represented in the sequencing data. Using this approach we detected almost all of the previously known chicken miRNAs and their respective miRNA* sequences. In addition we discovered 449 new chicken miRNAs including 88 miRNA candidates. Of these, 430 miRNAs appear to be specific to the avian lineage. Another six new miRNAs had evidence of evolutionary conservation in at least one vertebrate species outside of the bird lineage. The remaining 13 putative miRNAs appear to represent chicken orthologs of known vertebrate miRNAs. We discovered 39 additional putative miRNA candidates originating from miRNA generating intronic sequences known as mirtrons.

Origin, evolution, and biological role of miRNA cluster in DLK-DIO3 genomic region in placental mammals.

MicroRNAs (miRNAs) are a rapidly growing family of small regulatory RNAs modulating gene expression in plants and animals. In animals, most of the miRNAs discovered in early studies were found to be evolutionarily conserved across the whole kingdom. More recent studies, however, have identified many miRNAs that are specific to a particular group of organisms or even a single species. These present a question about evolution of the individual miRNAs and their role in establishing and maintaining lineage-specific functions and characteristics. In this study, we describe a detailed analysis of the miRNA cluster (hereafter mir-379/mir-656 cluster) located within the imprinted DLK-DIO3 region on human chromosome 14. We show that orthologous miRNA clusters are present in all sequenced genomes of the placental (eutherian) mammals but not in the marsupial (metatherian), monotreme (prototherian), or any other vertebrate genomes. We provide evidence that the locus encompassing this cluster emerged in an early eutherian ancestor prior to the radiation of modern placental mammals by tandem duplication of the ancient precursor sequence. The original amplified cluster may have contained in excess of 250 miRNA precursor sequences, most of which now appear to be inactive. Examination of the eutherian genomes showed that the cluster has been maintained in evolution for approximately 100 Myr. Analysis of genes that contain predicted evolutionarily conserved targets for miRNAs from this cluster revealed significant overrepresentation of the Gene Ontology terms associated with biological processes such as neurogenesis, embryonic development, transcriptional regulation, and RNA metabolism. Consistent with these findings, a survey of the miRNA expression data within the cluster demonstrates a strong bias toward brain and placenta samples from adult organisms and some embryonic tissues. Our results suggest that emergence of the mir-379/mir-656 miRNA cluster was one of the factors that facilitated evolution of the placental mammals. Overrepresentation of genes involved in regulation of neurogenesis among predicted miRNAs targets indicates an important role of the mir-379/mir-656 cluster in this biological process in the placental mammals.

Evidence for control of splicing by alternative RNA secondary structures in Dipteran homothorax pre-mRNA.

In a recent study that identified highly evolutionary conserved sequences in three genomes of Diptera species we described an ultraconserved element found at an internal exon-intron junction of the Drosophila melanogaster homothorax (hth) gene that appeared to be involved in the control of hth pre-mRNA splicing. We also discussed a possible role of RNA secondary structure at this site in the regulation of hth pre-mRNA splicing. In this report we identify a shorter evolutionary conserved intronic element within the hth gene that is located downstream of the first element and has sequence complementarity to it. We demonstrate that intramolecular interactions between these two elements would give rise to alternative RNA secondary structures, which in turn may result in differential control of homothorax pre-mRNA splicing. We also provide additional comparative genomic data from several newly available insect genomes supporting our original conclusion that these conserved elements are important in the post-transcriptional regulation of homothorax gene expression in Diptera.

Ultraconserved elements in insect genomes: a highly conserved intronic sequence implicated in the control of homothorax mRNA splicing.

Recently, we identified a large number of ultraconserved (uc) sequences in noncoding regions of human, mouse, and rat genomes that appear to be essential for vertebrate and amniote ontogeny. Here, we used similar methods to identify ultraconserved genomic regions between the insect species Drosophila melanogaster and Drosophila pseudoobscura, as well as the more distantly related Anopheles gambiae. As with vertebrates, ultraconserved sequences in insects appear to occur primarily in intergenic and intronic sequences, and at intron-exon junctions. The sequences are significantly associated with genes encoding developmental regulators and transcription factors, but are less frequent and are smaller in size than in vertebrates. The longest identical, nongapped orthologous match between the three genomes was found within the homothorax (hth) gene. This sequence spans an internal exon-intron junction, with the majority located within the intron, and is predicted to form a highly stable stem-loop RNA structure. Real-time quantitative PCR analysis of different hth splice isoforms and Northern blotting showed that the conserved element is associated with a high incidence of intron retention in hth pre-mRNA, suggesting that the conserved intronic element is critically important in the post-transcriptional regulation of hth expression in Diptera.