Monocytes/macrophages have begun to emerge as key cellular modulators of brain homeostasis and central nervous system (CNS) disease. In the healthy brain, resident microglia are the predominant macrophage cell population; however, under conditions of blood-brain barrier leakage, peripheral monocytes/macrophages can infiltrate the brain and participate in CNS disease pathogenesis. Distinguishing these two populations is often challenging, owing to a paucity of universally accepted and reliable markers. To identify discriminatory marker sets for microglia and peripheral monocytes/macrophages, we employed a large meta-analytic approach using five published murine transcriptional datasets. Following hierarchical clustering, we filtered the top differentially expressed genes (DEGs) through a brain cell type-specific sequencing database, which led to the identification of eight microglia and eight peripheral monocyte/macrophage markers. We then validated their differential expression, leveraging a published single cell RNA sequencing dataset and quantitative RT-PCR using freshly isolated microglia and peripheral monocytes/macrophages from two different mouse strains. We further verified the translation of these DEGs at the protein level. As top microglia DEGs, we identified P2ry12, Tmem119, Slc2a5 and Fcrls, whereas Emilin2, Gda, Hp and Sell emerged as the best DEGs for identifying peripheral monocytes/macrophages. Lastly, we evaluated their utility in discriminating monocyte/macrophage populations in the setting of brain pathology (glioma), and found that these DEG sets distinguished glioma-associated microglia from macrophages in both RCAS and GL261 mouse models of glioblastoma. Taken together, this unbiased bioinformatic approach facilitated the discovery of a robust set of microglia and peripheral monocyte/macrophage expression markers to discriminate these monocyte populations in both health and disease.

We took advantage of the unusual genomic organization of the ciliate Oxytricha trifallax to screen for eukaryotic non-coding RNA (ncRNA) genes. Ciliates have two types of nuclei: a germ line micronucleus that is usually transcriptionally inactive, and a somatic macronucleus that contains a reduced, fragmented and rearranged genome that expresses all genes required for growth and asexual reproduction. In some ciliates including Oxytricha, the macronuclear genome is particularly extreme, consisting of thousands of tiny 'nanochromosomes', each of which usually contains only a single gene. Because the organism itself identifies and isolates most of its genes on single-gene nanochromosomes, nanochromosome structure could facilitate the discovery of unusual genes or gene classes, such as ncRNA genes. Using a draft Oxytricha genome assembly and a custom-written protein-coding genefinding program, we identified a subset of nanochromosomes that lack any detectable protein-coding gene, thereby strongly enriching for nanochromosomes that carry ncRNA genes. We found only a small proportion of non-coding nanochromosomes, suggesting that Oxytricha has few independent ncRNA genes besides homologs of already known RNAs. Other than new members of known ncRNA classes including C/D and H/ACA snoRNAs, our screen identified one new family of small RNA genes, named the Arisong RNAs, which share some of the features of small nuclear RNAs.

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor "silent" germline micronuclear genome by a process of "unscrambling" and fragmentation. The tiny macronuclear "nanochromosomes" typically encode single, protein-coding genes (a small portion, 10%, encode 2-8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.

Exome sequencing (ES) has become an important tool in pediatric genomic medicine, improving identification of disease-associated variation due to assay breadth. Depth is also afforded by ES, enabling detection of lower-frequency mosaic variation compared to Sanger sequencing in the studied tissue, thus enhancing diagnostic yield. Within a pediatric tertiary-care hospital, we report two years of clinical ES data from probands evaluated for genetic disease to assess diagnostic yield, characteristics of causal variants, and prevalence of mosaicism among disease-causing variants. Exome-derived, phenotype-driven variant data from 357 probands was analyzed concurrent with parental ES data, when available. Blood was the source of nucleic acid. Sequence read alignments were manually reviewed for all assessed variants. Sanger sequencing was used for suspected de novo or mosaic variation. Clinical provider notes were reviewed to determine concordance between laboratory-reported data and the ordering provider's interpretation of variant-associated disease causality. Laboratory-derived diagnostic yield and provider-substantiated diagnoses had 91.4% concordance. The cohort returned 117 provider-substantiated diagnoses among 115 probands for a diagnostic yield of 32.2%. De novo variants represented 64.9% of disease-associated variation within trio analyses. Among the 115 probands, five harbored disease-associated somatic mosaic variation. Two additional probands were observed to inherit a disease-associated variant from an unaffected mosaic parent. Among inheritance patterns, de novo variation was the most frequent disease etiology. Somatic mosaicism is increasingly recognized as a significant contributor to genetic disease, particularly with increased sequence depth attainable from ES. This report highlights the potential and importance of detecting mosaicism in ES.

Bacterial binding to platelets is a key step in the development of infective endocarditis (IE). Sialic acid, a common terminal carbohydrate on host glycans, is the major receptor for streptococci on platelets. So far, all defined interactions between streptococci and sialic acid on platelets are mediated by serine-rich repeat proteins (SRRPs). However, we identified Streptococcus oralis subsp. oralis IE-isolates that bind sialic acid but lack SRRPs. In addition to binding sialic acid, some SRRP- isolates also bind the cryptic receptor β-1,4-linked galactose through a yet unknown mechanism. Using comparative genomics, we identified a novel sialic acid-binding adhesin, here named AsaA (associated with sialic acid adhesion A), present in IE-isolates lacking SRRPs. We demonstrated that S. oralis subsp. oralis AsaA is required for binding to platelets in a sialic acid-dependent manner. AsaA comprises a non-repeat region (NRR), consisting of a FIVAR/CBM and two Siglec-like and Unique domains, followed by 31 DUF1542 domains. When recombinantly expressed, Siglec-like and Unique domains competitively inhibited binding of S. oralis subsp. oralis and directly interacted with sialic acid on platelets. We further demonstrated that AsaA impacts the pathogenesis of S. oralis subsp. oralis in a rabbit model of IE. Additionally, we found AsaA orthologues in other IE-causing species and demonstrated that the NRR of AsaA from Gemella haemolysans blocked binding of S. oralis subsp. oralis, suggesting that AsaA contributes to the pathogenesis of multiple IE-causing species. Finally, our findings provide evidence that sialic acid is a key factor for bacterial-platelets interactions in a broader range of species than previously appreciated, highlighting its potential as a therapeutic target.

The advantages of Pacific Biosciences (PacBio) single-molecule real-time (SMRT) technology include long reads, low systematic bias, and high consensus read accuracy. Here we use these attributes to improve on the genome annotation of the parasitic hookworm Ancylostoma ceylanicum using PacBio RNA-Seq.

'Orang-utan' is derived from a Malay term meaning 'man of the forest' and aptly describes the southeast Asian great apes native to Sumatra and Borneo. The orang-utan species, Pongo abelii (Sumatran) and Pongo pygmaeus (Bornean), are the most phylogenetically distant great apes from humans, thereby providing an informative perspective on hominid evolution. Here we present a Sumatran orang-utan draft genome assembly and short read sequence data from five Sumatran and five Bornean orang-utan genomes. Our analyses reveal that, compared to other primates, the orang-utan genome has many unique features. Structural evolution of the orang-utan genome has proceeded much more slowly than other great apes, evidenced by fewer rearrangements, less segmental duplication, a lower rate of gene family turnover and surprisingly quiescent Alu repeats, which have played a major role in restructuring other primate genomes. We also describe a primate polymorphic neocentromere, found in both Pongo species, emphasizing the gradual evolution of orang-utan genome structure. Orang-utans have extremely low energy usage for a eutherian mammal, far lower than their hominid relatives. Adding their genome to the repertoire of sequenced primates illuminates new signals of positive selection in several pathways including glycolipid metabolism. From the population perspective, both Pongo species are deeply diverse; however, Sumatran individuals possess greater diversity than their Bornean counterparts, and more species-specific variation. Our estimate of Bornean/Sumatran speciation time, 400,000 years ago, is more recent than most previous studies and underscores the complexity of the orang-utan speciation process. Despite a smaller modern census population size, the Sumatran effective population size (N(e)) expanded exponentially relative to the ancestral N(e) after the split, while Bornean N(e) declined over the same period. Overall, the resources and analyses presented here offer new opportunities in evolutionary genomics, insights into hominid biology, and an extensive database of variation for conservation efforts.

Efforts to characterize the human transcriptome have largely been limited to blood, urine, and tissue analyses (i.e., normally sterile materials). We report here an extraction protocol using commercially available reagents to obtain high-yield, reverse-transcribable RNA from human stool. Quantitative reverse transcriptase polymerase chain reactions demonstrated minimal intra-specimen but considerable intra-subject variability over time of transcripts for interleukin-6 (IL-6), IL-8, epidermal growth factor (EGF), calprotectin, and glyceraldehyde-3-phosphate dehydrogenase (GAPDH). This technique now expands opportunities to use the human fecal transcriptome to characterize gastrointestinal pathophysiology.

Many genetic studies of intractable epilepsy in pediatric patients primarily focus on inherited, constitutional genetic deficiencies identified in patient blood. Recently, studies have revealed somatic mosaicism associated with epilepsy in which genetic variants are present only in a subset of brain cells. We hypothesize that tissue-specific, somatic mosaicism represents an important genetic etiology in epilepsy and aim to discover somatic alterations in epilepsy-affected brain tissue.

Expression profiling of distinct central nervous system (CNS) cell populations has been employed to facilitate disease classification and to provide insights into the molecular basis of brain pathology. One important cell type implicated in a wide variety of CNS disease states is the resident brain macrophage (microglia). In these studies, microglia are often isolated from dissociated brain tissue by flow sorting procedures [fluorescence-activated cell sorting (FACS)] or from postnatal glial cultures by mechanic isolation. Given the highly dynamic and state-dependent functions of these cells, the use of FACS or short-term culture methods may not accurately capture the biology of brain microglia. In the current study, we performed RNA-sequencing using Cx3cr1(+/GFP) labeled microglia isolated from the brainstem of 6-week-old mice to compare the transcriptomes of FACS-sorted versus laser capture microdissection (LCM). While both isolation techniques resulted in a large number of shared (common) transcripts, we identified transcripts unique to FACS-isolated and LCM-captured microglia. In particular, ∼50% of these LCM-isolated microglial transcripts represented genes typically associated with neurons and glia. While these transcripts clearly localized to microglia using complementary methods, they were not translated into protein. Following the induction of murine experimental autoimmune encephalomyelitis, increased oligodendrocyte and neuronal transcripts were detected in microglia, while only the myelin basic protein oligodendrocyte transcript was increased in microglia after traumatic brain injury. Collectively, these findings have implications for the design and interpretation of microglia transcriptome-based investigations.

Next-generation sequencing has been used to infer the clonality of heterogeneous tumor samples. These analyses yield specific predictions-the population frequency of individual clones, their genetic composition, and their evolutionary relationships-which we set out to test by sequencing individual cells from three subjects diagnosed with secondary acute myeloid leukemia, each of whom had been previously characterized by whole genome sequencing of unfractionated tumor samples. Single-cell mutation profiling strongly supported the clonal architecture implied by the analysis of bulk material. In addition, it resolved the clonal assignment of single nucleotide variants that had been initially ambiguous and identified areas of previously unappreciated complexity. Accordingly, we find that many of the key assumptions underlying the analysis of tumor clonality by deep sequencing of unfractionated material are valid. Furthermore, we illustrate a single-cell sequencing strategy for interrogating the clonal relationships among known variants that is cost-effective, scalable, and adaptable to the analysis of both hematopoietic and solid tumors, or any heterogeneous population of cells.

Low-grade brain tumors (pilocytic astrocytomas) arising in the neurofibromatosis type 1 (NF1) inherited cancer predisposition syndrome are hypothesized to result from a combination of germline and acquired somatic NF1 tumor suppressor gene mutations. However, genetically engineered mice (GEM) in which mono-allelic germline Nf1 gene loss is coupled with bi-allelic somatic (glial progenitor cell) Nf1 gene inactivation develop brain tumors that do not fully recapitulate the neuropathological features of the human condition. These observations raise the intriguing possibility that, while loss of neurofibromin function is necessary for NF1-associated low-grade astrocytoma development, additional genetic changes may be required for full penetrance of the human brain tumor phenotype. To identify these potential cooperating genetic mutations, we performed whole-genome sequencing (WGS) analysis of three NF1-associated pilocytic astrocytoma (PA) tumors. We found that the mechanism of somatic NF1 loss was different in each tumor (frameshift mutation, loss of heterozygosity, and methylation). In addition, tumor purity analysis revealed that these tumors had a high proportion of stromal cells, such that only 50%-60% of cells in the tumor mass exhibited somatic NF1 loss. Importantly, we identified no additional recurrent pathogenic somatic mutations, supporting a model in which neuroglial progenitor cell NF1 loss is likely sufficient for PA formation in cooperation with a proper stromal environment.

Brain tumors (gliomas) contain large populations of infiltrating macrophages and recruited microglia, which in experimental murine glioma models promote tumor formation and progression. Among the barriers to understanding the contributions of these stromal elements to high-grade glioma (glioblastoma; GBM) biology is the relative paucity of tools to characterize infiltrating macrophages and resident microglia. In this study, we leveraged multiple RNA analysis platforms to identify new monocyte markers relevant to GBM patient outcome.

Gangliogliomas (WHO grade I) are rare tumors affecting the central nervous system and are most frequently observed in children. Next-generation sequencing of tumors is being utilized at an increasing rate in both research and clinical settings to characterize the genetic factors that drive tumorigenesis. Here, we report a rare BRAF somatic mutation (NM_004333.4:c.1794_1796dupTAC; p.Thr599dup) in the tumor genome from a pediatric patient in her late teens, who was initially diagnosed with low-grade ganglioglioma at age 13. This duplication of 3 nt introduces a second threonine residue at amino acid 599 of the BRAF protein. Based on previous studies, this variant is likely to increase kinase activity, similar to the well-characterized BRAF p.Val600Glu (V600E) pathogenic variant. In addition, although the p.T599dup somatic mutation has been documented rarely in human cancers, the variant has not been previously reported in ganglioglioma. The identification of this variant presents an opportunity to consider targeted therapy (e.g., BRAF inhibitor) for this patient.

Because of the constant threat posed by emerging infectious diseases and the limitations of existing approaches used to identify new pathogens, there is a great demand for new technological methods for viral discovery. We describe herein a DNA microarray-based platform for novel virus identification and characterization. Central to this approach was a DNA microarray designed to detect a wide range of known viruses as well as novel members of existing viral families; this microarray contained the most highly conserved 70mer sequences from every fully sequenced reference viral genome in GenBank. During an outbreak of severe acute respiratory syndrome (SARS) in March 2003, hybridization to this microarray revealed the presence of a previously uncharacterized coronavirus in a viral isolate cultivated from a SARS patient. To further characterize this new virus, approximately 1 kb of the unknown virus genome was cloned by physically recovering viral sequences hybridized to individual array elements. Sequencing of these fragments confirmed that the virus was indeed a new member of the coronavirus family. This combination of array hybridization followed by direct viral sequence recovery should prove to be a general strategy for the rapid identification and characterization of novel viruses and emerging infectious disease.

Genomic studies in non-domestic avian models, such as the California condor and white-throated sparrow, can lead to more comprehensive conservation plans and provide clues for understanding mechanisms affecting genetic variation, adaptation and evolution.Developing genomic tools and resources including genomic libraries and a genetic map of the California condor is a prerequisite for identification of candidate loci for a heritable embryonic lethal condition. The white-throated sparrow exhibits a stable genetic polymorphism (i.e. chromosomal rearrangements) associated with variation in morphology, physiology, and behavior (e.g., aggression, social behavior, sexual behavior, parental care).In this paper we outline the utility of these species as well as report on recent advances in the study of their genomes.

The human reference genome assembly plays a central role in nearly all aspects of today's basic and clinical research. GRCh38 is the first coordinate-changing assembly update since 2009; it reflects the resolution of roughly 1000 issues and encompasses modifications ranging from thousands of single base changes to megabase-scale path reorganizations, gap closures, and localization of previously orphaned sequences. We developed a new approach to sequence generation for targeted base updates and used data from new genome mapping technologies and single haplotype resources to identify and resolve larger assembly issues. For the first time, the reference assembly contains sequence-based representations for the centromeres. We also expanded the number of alternate loci to create a reference that provides a more robust representation of human population variation. We demonstrate that the updates render the reference an improved annotation substrate, alter read alignments in unchanged regions, and impact variant interpretation at clinically relevant loci. We additionally evaluated a collection of new de novo long-read haploid assemblies and conclude that although the new assemblies compare favorably to the reference with respect to continuity, error rate, and gene completeness, the reference still provides the best representation for complex genomic regions and coding sequences. We assert that the collected updates in GRCh38 make the newer assembly a more robust substrate for comprehensive analyses that will promote our understanding of human biology and advance our efforts to improve health.

Chromatin diminution is the programmed elimination of specific DNA sequences during development. It occurs in diverse species, but the function(s) of diminution and the specificity of sequence loss remain largely unknown. Diminution in the nematode Ascaris suum occurs during early embryonic cleavages and leads to the loss of germline genome sequences and the formation of a distinct genome in somatic cells. We found that ∼43 Mb (∼13%) of genome sequence is eliminated in A. suum somatic cells, including ∼12.7 Mb of unique sequence. The eliminated sequences and location of the DNA breaks are the same in all somatic lineages from a single individual and between different individuals. At least 685 genes are eliminated. These genes are preferentially expressed in the germline and during early embryogenesis. We propose that diminution is a mechanism of germline gene regulation that specifically removes a large number of genes involved in gametogenesis and early embryogenesis.

Genome-wide association studies (GWAS) have identified >500 common variants associated with quantitative metabolic traits, but in aggregate such variants explain at most 20-30% of the heritable component of population variation in these traits. To further investigate the impact of genotypic variation on metabolic traits, we conducted re-sequencing studies in >6,000 members of a Finnish population cohort (The Northern Finland Birth Cohort of 1966 [NFBC]) and a type 2 diabetes case-control sample (The Finland-United States Investigation of NIDDM Genetics [FUSION] study). By sequencing the coding sequence and 5' and 3' untranslated regions of 78 genes at 17 GWAS loci associated with one or more of six metabolic traits (serum levels of fasting HDL-C, LDL-C, total cholesterol, triglycerides, plasma glucose, and insulin), and conducting both single-variant and gene-level association tests, we obtained a more complete understanding of phenotype-genotype associations at eight of these loci. At all eight of these loci, the identification of new associations provides significant evidence for multiple genetic signals to one or more phenotypes, and at two loci, in the genes ABCA1 and CETP, we found significant gene-level evidence of association to non-synonymous variants with MAF<1%. Additionally, two potentially deleterious variants that demonstrated significant associations (rs138726309, a missense variant in G6PC2, and rs28933094, a missense variant in LIPC) were considerably more common in these Finnish samples than in European reference populations, supporting our prior hypothesis that deleterious variants could attain high frequencies in this isolated population, likely due to the effects of population bottlenecks. Our results highlight the value of large, well-phenotyped samples for rare-variant association analysis, and the challenge of evaluating the phenotypic impact of such variants.

Our efforts to prevent and treat breast cancer are significantly impeded by a lack of knowledge of the biology and developmental genetics of the normal mammary gland. In order to provide the specimens that will facilitate such an understanding, The Susan G. Komen for the Cure Tissue Bank at the IU Simon Cancer Center (KTB) was established. The KTB is, to our knowledge, the only biorepository in the world prospectively established to collect normal, healthy breast tissue from volunteer donors. As a first initiative toward a molecular understanding of the biology and developmental genetics of the normal mammary gland, the effect of the menstrual cycle and hormonal contraceptives on DNA expression in the normal breast epithelium was examined.