Human and mouse gene structure: comparative analysis and application to exon prediction

Xenomake: a pipeline for processing and sorting xenograft reads from spatial transcriptomic experiments

SUMMARY

Xenograft models are attractive models that mimic human tumor biology and permit one to perturb the tumor microenvironment and study its drug response. Spatially resolved transcriptomics (SRT) provides a powerful way to study the organization of xenograft models, but currently there is a lack of specialized pipeline for processing xenograft reads originated from SRT experiments. Xenomake is a standalone pipeline for the automated handling of spatial xenograft reads. Xenomake handles read processing, alignment, xenograft read sorting, and connects well with downstream spatial analysis packages. We additionally show that Xenomake can correctly assign organism-specific reads, reduce sparsity of data by increasing gene counts, while maintaining biological relevance for studies.

AVAILABILITY AND IMPLEMENTATION

Xenomake is an open-source program that is available on Github (https://github.com/qianzhulab/Xenomake). Complete documentation can be found at the link.

5. Collaborative Study on the Genetics of Alcoholism: Functional genomics

Alcohol Use Disorder is a complex genetic disorder, involving genetic, neural, and environmental factors, and their interactions. The Collaborative Study on the Genetics of Alcoholism (COGA) has been investigating these factors and identified putative alcohol use disorder risk genes through genome-wide association studies. In this review, we describe advances made by COGA in elucidating the functional changes induced by alcohol use disorder risk genes using multimodal approaches with human cell lines and brain tissue. These studies involve investigating gene regulation in lymphoblastoid cells from COGA participants and in post-mortem brain tissues. High throughput reporter assays are being used to identify single nucleotide polymorphisms in which alternate alleles differ in driving gene expression. Specific single nucleotide polymorphisms (both coding or noncoding) have been modeled using induced pluripotent stem cells derived from COGA participants to evaluate the effects of genetic variants on transcriptomics, neuronal excitability, synaptic physiology, and the response to ethanol in human neurons from individuals with and without alcohol use disorder. We provide a perspective on future studies, such as using polygenic risk scores and populations of induced pluripotent stem cell-derived neurons to identify signaling pathways related with responses to alcohol. Starting with genes or loci associated with alcohol use disorder, COGA has demonstrated that integration of multimodal data within COGA participants and functional studies can reveal mechanisms linking genomic variants with alcohol use disorder, and potential targets for future treatments.

Heuristic Pairwise Alignment in Database Environments

Biological data have gained wider recognition during the last few years, although managing and processing these data in an efficient way remains a challenge in many areas. Increasingly, more DNA sequence databases can be accessed; however, most algorithms on these sequences are performed outside of the database with different bioinformatics software. In this article, we propose a novel approach for the comparative analysis of sequences, thereby defining heuristic pairwise alignment inside the database environment. This method takes advantage of the benefits provided by the database management system and presents a way to exploit similarities in data sets to quicken the alignment algorithm. We work with the column-oriented MonetDB, and we further discuss the key benefits of this database system in relation to our proposed heuristic approach.

PhIP-Seq Reveals Autoantibodies for Ubiquitously Expressed Antigens in Viral Myocarditis

Simple Summary

Myocarditis is the inflammation of the heart muscle, and viral infections are a common cause of this disease. Myocarditis in some patients can progress to dilated cardiomyopathy (DCM). The mouse model of coxsackievirus B3 (CVB3) is commonly used to understand this disease progression in DCM patients. In this paper, we have attempted to analyze antibodies for heart antigens that could be produced as a result of heart damage in animals infected with CVB3 using a technique called Phage ImmunoPrecipitation Sequencing (PhIP-Seq). The analyses led us to identify antibodies for several proteins that were not previously reported that may have relevance to human disease.

Abstract

Enteroviruses such as group B coxsackieviruses (CVB) are commonly suspected as causes of myocarditis that can lead to dilated cardiomyopathy (DCM), and the mouse model of CVB3 myocarditis is routinely used to understand DCM pathogenesis. Mechanistically, autoimmunity is suspected due to the presence of autoantibodies for select antigens. However, their role continues to be enigmatic, which also raises the question of whether the breadth of autoantibodies is sufficiently characterized. Here, we attempted to comprehensively analyze the autoantibody repertoire using Phage ImmunoPrecipitation Sequencing (PhIP-Seq), a versatile and high-throughput platform, in the mouse model of CVB3 myocarditis. First, PhIP-Seq analysis using the VirScan library revealed antibody reactivity only to CVB3 in the infected group but not in controls, thus validating the technique in this model. Second, using the mouse peptide library, we detected autoantibodies to 32 peptides from 25 proteins in infected animals that are ubiquitously expressed and have not been previously reported. Third, by using ELISA as a secondary assay, we confirmed antibody reactivity in sera from CVB3-infected animals to cytochrome c oxidase assembly factor 4 homolog (COA4) and phosphoinositide-3-kinase adaptor protein 1 (PIK3AP1), indicating the specificity of antibody detection by PhIP-Seq technology. Fourth, we noted similar antibody reactivity patterns in CVB3 and CVB4 infections, suggesting that the COA4- and PIK3AP1-reactive antibodies could be common to multiple CVB infections. The specificity of the autoantibodies was affirmed with influenza-infected animals that showed no reactivity to any of the antigens tested. Taken together, our data suggest that the autoantibodies identified by PhIP-Seq may have relevance to CVB pathogenesis, with a possibility that similar reactivity could be expected in human DCM patients.

Does the Pachytene Checkpoint, a Feature of Meiosis, Filter Out Mistakes in Double-Strand DNA Break Repair and as a side-Effect Strongly Promote Adaptive Speciation?

This essay aims to explain two biological puzzles: why eukaryotic transcription units are composed of short segments of coding DNA interspersed with long stretches of non-coding (intron) DNA, and the near ubiquity of sexual reproduction. As is well known, alternative splicing of its coding sequences enables one transcription unit to produce multiple variants of each encoded protein. Additionally, padding transcription units with non-coding DNA (often many thousands of base pairs long) provides a readily evolvable way to set how soon in a cell cycle the various mRNAs will begin being expressed and the total amount of mRNA that each transcription unit can make during a cell cycle. This regulation complements control via the transcriptional promoter and facilitates the creation of complex eukaryotic cell types, tissues, and organisms. However, it also makes eukaryotes exceedingly vulnerable to double-strand DNA breaks, which end-joining break repair pathways can repair incorrectly. Transcription units cover such a large fraction of the genome that any mis-repair producing a reorganized chromosome has a high probability of destroying a gene. During meiosis, the synaptonemal complex aligns homologous chromosome pairs and the pachytene checkpoint detects, selectively arrests, and in many organisms actively destroys gamete-producing cells with chromosomes that cannot adequately synapse; this creates a filter favoring transmission to the next generation of chromosomes that retain the parental organization, while selectively culling those with interrupted transcription units. This same meiotic checkpoint, reacting to accidental chromosomal reorganizations inflicted by error-prone break repair, can, as a side effect, provide a mechanism for the formation of new species in sympatry. It has been a long-standing puzzle how something as seemingly maladaptive as hybrid sterility between such new species can arise. I suggest that this paradox is resolved by understanding the adaptive importance of the pachytene checkpoint, as outlined above.

The Coming of Age of Preclinical Models of MDS

Myelodysplastic syndromes (MDS) are a heterogeneous group of clonal bone-marrow diseases with ineffective hematopoiesis resulting in cytopenias and morphologic dysplasia of hematopoietic cells. MDS carry a wide spectrum of genetic abnormalities, ranging from chromosomal abnormalities such as deletions/additions, to recurrent mutations affecting the spliceosome, epigenetic modifiers, or transcription factors. As opposed to AML, research in MDS has been hindered by the lack of preclinical models that faithfully replicate the complexity of the disease and capture the heterogeneity. The complex molecular landscape of the disease poses a unique challenge when creating transgenic mouse-models. In addition, primary MDS cells are difficult to manipulate ex vivo limiting in vitro studies and resulting in a paucity of cell lines and patient derived xenograft models. In recent years, progress has been made in the development of both transgenic and xenograft murine models advancing our understanding of individual contributors to MDS pathology as well as the complex primary interplay of genetic and microenvironment aberrations. We here present a comprehensive review of these transgenic and xenograft models for MDS and future directions.

The Cerebellar Involvement in Autism Spectrum Disorders: From the Social Brain to Mouse Models

Autism spectrum disorders (ASD) are pervasive neurodevelopmental disorders that include a variety of forms and clinical phenotypes. This heterogeneity complicates the clinical and experimental approaches to ASD etiology and pathophysiology. To date, a unifying theory of these diseases is still missing. Nevertheless, the intense work of researchers and clinicians in the last decades has identified some ASD hallmarks and the primary brain areas involved. Not surprisingly, the areas that are part of the so-called “social brain”, and those strictly connected to them, were found to be crucial, such as the prefrontal cortex, amygdala, hippocampus, limbic system, and dopaminergic pathways. With the recent acknowledgment of the cerebellar contribution to cognitive functions and the social brain, its involvement in ASD has become unmistakable, though its extent is still to be elucidated. In most cases, significant advances were made possible by recent technological developments in structural/functional assessment of the human brain and by using mouse models of ASD. Mouse models are an invaluable tool to get insights into the molecular and cellular counterparts of the disease, acting on the specific genetic background generating ASD-like phenotype. Given the multifaceted nature of ASD and related studies, it is often difficult to navigate the literature and limit the huge content to specific questions. This review fulfills the need for an organized, clear, and state-of-the-art perspective on cerebellar involvement in ASD, from its connections to the social brain areas (which are the primary sites of ASD impairments) to the use of monogenic mouse models.

Modeling injury and repair in kidney organoids reveals that homologous recombination governs tubular intrinsic repair

Kidneys have the capacity for intrinsic repair, preserving kidney architecture with return to a basal state after tubular injury. When injury is overwhelming or repetitive, however, that capacity is exceeded and incomplete repair results in fibrotic tissue replacing normal kidney parenchyma. Loss of nephrons correlates with reduced kidney function, which defines chronic kidney disease (CKD) and confers substantial morbidity and mortality to the worldwide population. Despite the identification of pathways involved in intrinsic repair, limited treatments for CKD exist, partly because of the limited throughput and predictivity of animal studies. Here, we showed that kidney organoids can model the transition from intrinsic to incomplete repair. Single-nuclear RNA sequencing of kidney organoids after cisplatin exposure identified 159 differentially expressed genes and 29 signal pathways in tubular cells undergoing intrinsic repair. Homology-directed repair (HDR) genes including Fanconi anemia complementation group D2 (FANCD2) and RAD51 recombinase (RAD51) were transiently up-regulated during intrinsic repair but were down-regulated in incomplete repair. Single cellular transcriptomics in mouse models of obstructive and hemodynamic kidney injury and human kidney samples of immune-mediated injury validated HDR gene up-regulation during tubular repair. Kidney biopsy samples with tubular injury and varying degrees of fibrosis confirmed loss of FANCD2 during incomplete repair. Last, we performed targeted drug screening that identified the DNA ligase IV inhibitor, SCR7, as a therapeutic candidate that rescued FANCD2/RAD51-mediated repair to prevent the progression of CKD in the cisplatin-induced organoid injury model. Our findings demonstrate the translational utility of kidney organoids to identify pathologic pathways and potential therapies.

Progress in Vocal Fold Regenerative Biomaterials: An Immunological Perspective

Vocal folds, housed in the upper respiratory tract, are important to daily breathing, speech and swallowing functions. Irreversible changes to the vocal fold mucosae, such as scarring and atrophy, require a regenerative medicine approach to promote a controlled regrowth of the extracellular matrix (ECM)-rich mucosa. Various biomaterial systems have been engineered with an emphasis on stimulating local vocal fold fibroblasts to produce new ECM. At the same time, it is imperative to limit the foreign body reaction and associated immune components that can hinder the integration of the biomaterial into the host tissue. Modern biomaterial designs have become increasingly focused on actively harnessing the immune system to accelerate and optimize the process of tissue regeneration. An array of physical and chemical biomaterial parameters have been reported to effectively modulate local immune cells, such as macrophages, to initiate tissue repair, stimulate ECM production, promote biomaterial-tissue integration, and restore the function of the vocal folds. In this perspective paper, the unique immunological profile of the vocal folds will first be reviewed. Key physical and chemical biomaterial properties relevant to immunomodulation will then be highlighted and discussed. A further examination of the physicochemical properties of recent vocal fold biomaterials will follow to generate deeper insights into corresponding immune-related outcomes. Lastly, a perspective will be offered on the opportunity of integrating material-led immunomodulatory strategies into future vocal fold tissue engineering therapies.

The aryl hydrocarbon receptor activates ceramide biosynthesis in mice contributing to hepatic lipogenesis

Aryl hydrocarbon receptor (AHR) activation via 2,3,7,8-tetrachlorodibenzofuran (TCDF) induces the accumulation of hepatic lipids. Here we report that AHR activation by TCDF (24  μg/kg body weight given orally for five days) induced significant elevation of hepatic lipids including ceramides in mice, was associated with increased expression of key ceramide biosynthetic genes, and increased activity of their respective enzymes. Results from chromatin immunoprecipitation (ChIP), electrophoretic mobility shift assay (EMSA) and cell-based reporter luciferase assays indicated that AHR directly activated the serine palmitoyltransferase long chain base subunit 2 (Sptlc2, encodes serine palmitoyltransferase 2 (SPT2)) gene whose product catalyzes the initial rate-limiting step in de novo sphingolipid biosynthesis. Hepatic ceramide accumulation was further confirmed by mass spectrometry-based lipidomics. Taken together, our results revealed that AHR activation results in the up-regulation of Sptlc2, leading to ceramide accumulation, thus promoting lipogenesis, which can induce hepatic lipid accumulation.

New gene discoveries in skeletal diseases with short stature

In the last decade, the widespread use of massively parallel sequencing has considerably boosted the number of novel gene discoveries in monogenic skeletal diseases with short stature. Defects in genes playing a role in the maintenance and function of the growth plate, the site of longitudinal bone growth, are a well-known cause of skeletal diseases with short stature. However, several genes involved in extracellular matrix composition or maintenance as well as genes partaking in various biological processes have also been characterized. This review aims to describe the latest genetic findings in spondyloepiphyseal dysplasias, spondyloepimetaphyseal dysplasias, and some monogenic forms of isolated short stature. Some examples of novel genetic mechanisms leading to skeletal conditions with short stature will be described. Strategies on how to successfully characterize novel skeletal phenotypes with short stature and genetic approaches to detect and validate novel gene-disease correlations will be discussed in detail. In summary, we review the latest gene discoveries underlying skeletal diseases with short stature and emphasize the importance of characterizing novel molecular mechanisms for genetic counseling, for an optimal management of the disease, and for therapeutic innovations.

Genetic susceptibility to nicotine addiction: Advances and shortcomings in our understanding of the CHRNA5/A3/B4 gene cluster contribution

Over the last decade, robust human genetic findings have been instrumental in elucidating the heritable basis of nicotine addiction (NA). They highlight coding and synonymous polymorphisms in a cluster on chromosome 15, encompassing the CHRNA5, CHRNA3 and CHRNB4 genes, coding for three subunits of the nicotinic acetylcholine receptor (nAChR). They have inspired an important number of preclinical studies, and will hopefully lead to the definition of novel drug targets for treating NA. Here, we review these candidate gene and genome-wide association studies (GWAS) and their direct implication in human brain function and NA-related phenotypes. We continue with a description of preclinical work in transgenic rodents that has led to a mechanistic understanding of several of the genetic hits. We also highlight important issues with regards to CHRNA3 and CHRNB4 where we are still lacking a dissection of their role in NA, including even in preclinical models. We further emphasize the use of human induced pluripotent stem cell-derived models for the analysis of synonymous and intronic variants on a human genomic background. Finally, we indicate potential avenues to further our understanding of the role of this human genetic variation. This article is part of the special issue on 'Contemporary Advances in Nicotine Neuropharmacology'.

Behavioral tests assessing neuropsychiatric phenotypes in adolescent mice reveal strain- and sex-specific effects

In humans, infancy and adolescence are associated with major changes in synaptic functions and ongoing maturation of neural networks, which underlie the major behavioral changes during these periods. Among adult cases with neuropsychiatric disorders including autism spectrum disorder, schizophrenia, attention deficit hyperactivity, and bipolar disorders, 50% have developed behavioral symptoms and received a diagnosis before 15 years of age. However, most of the behavioral studies in mice modeling neuropsychiatric phenotypes are performed in adult animals, missing valuable phenotypic information related to the effect of synaptic maturation during development. Here, we explored which behavioral experiments assessing neuropsychiatric phenotypes can be performed during a specific window of development in adolescent male and female C57BL/6N, DBA/2, and FVB/N mice that are typically used as background strains for generating genetically-modified mouse models. The three wild-type strains were evaluated across anxiety, social behaviors, and cognitive functions in order to cover the main behavioral impairments that occur in neuropsychiatric disorders. During adolescence, the three strains displayed significant differences under certain behavioral paradigms. In addition, C57BL/6N and FVB/N, but not DBA/2 mice revealed some sex-related differences. Our results provide new insights into discrete behaviors during development and emphasize the crucial importance of the genetic background, sex, and experimental settings in the age-dependent regulation of different behaviors.

Inducible aging in Hydra oligactis implicates sexual reproduction, loss of stem cells, and genome maintenance as major pathways

Freshwater polyps of the genus Hydra do not age. However, temperature stress induces aging and a shift from reproduction by asexual budding to sexual gamete production in a cold-sensitive (CS) strain of H. oligactis. We sequenced the transcriptome of a male CS strain before and after this life history shift and compared changes in gene expression relative to those seen in a cold-resistant (CR) strain that does not undergo a life history shift in response to altered temperature. We found that the switch from non-aging asexual reproduction to aging and sexual reproduction involves upregulation of genes not only involved in gametogenesis but also genes involved in cellular senescence, apoptosis, and DNA repair accompanied by a downregulation of genes involved in stem cell maintenance. These results suggest that aging is a byproduct of sexual reproduction-associated cellular reprogramming and underscore the power of these H. oligactis strains to identify intrinsic mechanisms of aging.

Molecular docking and dynamics simulations studies of OmpATb identifies four potential novel natural product-derived anti-Mycobacterium tuberculosis compounds

The outer membrane protein A (OmpATb) of Mycobacterium tuberculosis is a virulence factor that neutralizes the host pH to impede the uptake of hydrophilic antitubercular drugs. Identifying natural compounds with the potential to inhibit OmpATb could allow circumvention of the porin-like activities of OmpATb. Four potential leads comprising ZINC000003958185, ZINC000000157405, ZINC000000001392 and ZINC000034268676 were obtained by virtual screening of 6394 diverse natural products. Characterization of the binding interactions of the potential leads with OmpATb revealed nine critical residues comprising ARG86, LEU110, LEU113, LEU114, ALA115, PHE142, SER145, VAL146, and PHE151. Molecular dynamics simulations also revealed very stable protein-lead complexes. Most residues contributed lower binding energies to the overall molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) binding free energies of the interactions between the molecules and OmpATb protein. Induced Fit Docking (IFD) of the compounds regenerated poses of the molecular docking using AutoDock Vina. These molecules could be starting templates for designing inhibitors to bypass the pore mediating activities of OmpATb. Based on structural similarity, ZINC000034268676 was suggested as a potential scaffold for designing efflux pump inhibitors of the gate mediating activities of OmpATb and may enhance the uptake of hydrophilic drugs to reduce the duration time of tuberculosis treatment. Furthermore, structurally similar compounds available in the DrugBank database with a similarity threshold of 0.7 have been reported to exhibit antitubercular and anti-mycobacterial activities. These biomolecules can be further characterized experimentally to corroborate their antitubercular activity. Also, the skeletons of the molecules can be adopted as sub-structures for the design of future anti-mycobacterial drugs.

Mitigating Effect of 1-Palmitoyl-2-Linoleoyl-3-Acetyl-Rac-Glycerol (PLAG) on a Murine Model of 5-Fluorouracil-Induced Hematological Toxicity

5-Fluorouracil (5-FU) is an antimetabolite chemotherapy widely used for the treatment of various cancers. However, many cancer patients experience hematological side effects following 5-FU treatment. Here, we investigated the protective effects of 1-palmitoyl-2-linoleoyl-3-acetyl-rac-glycerol (PLAG) as a mitigator against 5-FU-induced hematologic toxicity, including neutropenia, monocytopenia, thrombocytopenia, and thrombocytosis, in Balb/c mice injected with 5-FU (100 mg/kg, i.p.). Administration of PLAG significantly and dose-dependently reduced the duration of neutropenia and improved the nadirs of absolute neutrophil counts (ANCs). Moreover, while the ANCs of all mice in the control fell to the severely neutropenic range, none of the mice in the PLAG 200 and 400 mg/kg-treated groups experienced severe neutropenia. Administration of PLAG significantly delayed the mean first day of monocytopenia and reduced the duration of monocytopenia. PLAG also effectively reduced extreme changes in platelet counts induced by 5-FU treatment, thus preventing 5-FU-induced thrombocytopenia and thrombocytosis. PLAG significantly decreased plasma levels of the chemokine (C–X–C motif) ligand 1 (CXCL1), CXCL2, interleukin (IL)-6, and C-reactive protein (CRP), which were elevated consistently with the occurrence time of neutropenia, monocytopenia, and thrombocytopenia. When compared with olive oil and palmitic linoleic hydroxyl glycerol (PLH), only PLAG effectively mitigated 5-FU-induced hematological toxicity, indicating that it has a distinctive mechanism of action. In conclusion, PLAG may have therapeutic potential as a mitigator for 5-FU-induced neutropenia and other hematological disorders.

Myocardial regeneration: role of epicardium and implicated genes

Lower invertebrates' hearts such as those of zebrafish have the capacity for scarless myocardial regeneration which is lost by mammalian hearts as they form a fibrotic scar tissue instead of regenerating the injured area. However, neonatal mammalian hearts have a remarkable capacity for regeneration highlighting conserved evolutionary mechanisms underlying such a process. Studies investigated the underlying mechanism of myocardial regeneration in species capable to do so, to see its applicability on mammals. The epicardium, the mesothelial outer layer of the vertebrate heart, has proven to play an important role in the process of repair and regeneration. It serves as an important source of smooth muscle cells, cardiac fibroblasts, endothelial cells, stem cells, and signaling molecules that are involved in this process. Here we review the role of the epicardium in myocardial regeneration focusing on the different involved; Activation, epithelial to mesenchymal transition, and differentiation. In addition, we will discuss its contributory role to different aspects that support myocardial regeneration. Of these we will discuss angiogenesis and the formation of a regenerate extracellular matrix. Moreover, we will discuss several factors that act on the epicardium to affect regeneration. Finally, we will highlight the utility of the epicardium as a mode of cell therapy in the treatment of myocardial injury.

Genomic data analysis workflows for tumors from patient-derived xenografts (PDXs): challenges and guidelines

BACKGROUND

Patient-derived xenograft (PDX) models are in vivo models of human cancer that have been used for translational cancer research and therapy selection for individual patients. The Jackson Laboratory (JAX) PDX resource comprises 455 models originating from 34 different primary sites (as of 05/08/2019). The models undergo rigorous quality control and are genomically characterized to identify somatic mutations, copy number alterations, and transcriptional profiles. Bioinformatics workflows for analyzing genomic data obtained from human tumors engrafted in a mouse host (i.e., Patient-Derived Xenografts; PDXs) must address challenges such as discriminating between mouse and human sequence reads and accurately identifying somatic mutations and copy number alterations when paired non-tumor DNA from the patient is not available for comparison.

RESULTS

We report here data analysis workflows and guidelines that address these challenges and achieve reliable identification of somatic mutations, copy number alterations, and transcriptomic profiles of tumors from PDX models that lack genomic data from paired non-tumor tissue for comparison. Our workflows incorporate commonly used software and public databases but are tailored to address the specific challenges of PDX genomics data analysis through parameter tuning and customized data filters and result in improved accuracy for the detection of somatic alterations in PDX models. We also report a gene expression-based classifier that can identify EBV-transformed tumors. We validated our analytical approaches using data simulations and demonstrated the overall concordance of the genomic properties of xenograft tumors with data from primary human tumors in The Cancer Genome Atlas (TCGA).

CONCLUSIONS

The analysis workflows that we have developed to accurately predict somatic profiles of tumors from PDX models that lack normal tissue for comparison enable the identification of the key oncogenic genomic and expression signatures to support model selection and/or biomarker development in therapeutic studies. A reference implementation of our analysis recommendations is available at https://github.com/TheJacksonLaboratory/PDX-Analysis-Workflows .

Whole-Genome Alignment and Comparative Annotation

Rapidly improving sequencing technology coupled with computational developments in sequence assembly are making reference-quality genome assembly economical. Hundreds of vertebrate genome assemblies are now publicly available, and projects are being proposed to sequence thousands of additional species in the next few years. Such dense sampling of the tree of life should give an unprecedented new understanding of evolution and allow a detailed determination of the events that led to the wealth of biodiversity around us. To gain this knowledge, these new genomes must be compared through genome alignment (at the sequence level) and comparative annotation (at the gene level). However, different alignment and annotation methods have different characteristics; before starting a comparative genomics analysis, it is important to understand the nature of, and biases and limitations inherent in, the chosen methods. This review is intended to act as a technical but high-level overview of the field that should provide this understanding. We briefly survey the state of the genome alignment and comparative annotation fields and potential future directions for these fields in a new, large-scale era of comparative genomics.

The Challenge to Search for New Nervous System Disease Biomarker Candidates: the Opportunity to Use the Proteogenomics Approach

Alzheimer's disease, Parkinson's disease, prion diseases, schizophrenia, and multiple sclerosis are the most common nervous system diseases, affecting millions of people worldwide. The current scientific literature associates these pathological conditions to abnormal expression levels of certain proteins, which in turn improved the knowledge concerning normal and affected brains. However, there is no available cure or preventive therapy for any of these disorders. Proteogenomics is a recent approach defined as the data integration of both nucleotide high-throughput sequencing and protein mass spectrometry technologies. In the last years, proteogenomics studies in distinct diseases have emerged as a strategy for the identification of uncharacterized proteoforms, which are all the different protein forms derived from a single gene. For many of these diseases, at least one protein used as biomarker presents more than one proteoform, which fosters the analysis of publicly available data focusing proteoforms. Given this context, we describe the most important biomarkers for each neurodegenerative disease and how genomics, transcriptomics, and proteomics separately contributed to unveil them. Finally, we present a selection of proteogenomics studies in which the combination of nucleotide and proteome high-throughput data, from cell lines or brain tissue samples, is used to uncover proteoforms not previously described. We believe that this new approach may improve our knowledge about nervous system diseases and brain function and an opportunity to identify new biomarker candidates.

Human Models Are Needed for Studying Human Neurodevelopmental Disorders

The analysis of animal models of neurological disease has been instrumental in furthering our understanding of neurodevelopment and brain diseases. However, animal models are limited in revealing some of the most fundamental aspects of development, genetics, pathology, and disease mechanisms that are unique to humans. These shortcomings are exaggerated in disorders that affect the brain, where the most significant differences between humans and animal models exist, and could underscore failures in targeted therapeutic interventions in affected individuals. Human pluripotent stem cells have emerged as a much-needed model system for investigating human-specific biology and disease mechanisms. However, questions remain regarding whether these cell-culture-based models are sufficient or even necessary. In this review, we summarize human-specific features of neurodevelopment and the most common neurodevelopmental disorders, present discrepancies between animal models and human diseases, demonstrate how human stem cell models can provide meaningful information, and discuss the challenges that exist in our pursuit to understand distinctively human aspects of neurodevelopment and brain disease. This information argues for a more thoughtful approach to disease modeling through consideration of the valuable features and limitations of each model system, be they human or animal, to mimic disease characteristics.

Fabrication of an anthropomorphic heterogeneous mouse phantom for multimodality medical imaging

This work presents a comprehensive methodology for constructing a tissue equivalent mouse phantom using image modeling and 3D printing technology. The phantom can be used in multimodality imaging and irradiation experiments, quality control, and management. Computed tomography (CT) images of a mouse were acquired and imported into 3D modeling software. A skeleton and skin shell models were segmented in the modeling software and manufactured using 3D printing technology. The bone model was constructed with VERO-WHITE printing material with additional ingredients, including a photosensitive resin, polyurethane epoxy resin, and acrylate. Acrylonitrile butadiene styrene resin material was used to construct the skin shell. The skin shell was attached to the skeleton and filled with a specially formulated gel to act as a soft tissue substitute. The gel consisted of agarose, micro-pearl powder, sodium chloride, and magnevist solution (gadopentetate dimeglumine). A micro-container filled with ¹⁸F-fluorodeoxyglucose (¹⁸F-FDG) radioactive tracer was placed in the abdomen for micro and human positron emission tomography (PET)/CT imaging. The mouse phantom had tissue equivalency in dose attenuation with x-rays and relaxation times with magnetic resonance imaging (MRI). The CT Hounsfield Unit (HU) range for the gel soft tissue material was 31-36 HU. The 3D printed bone mimetic material had equivalent tissue/bone contrast compared with in vivo mouse measurements with a mean value of 130  ±  10 HU. At different magnetic field strengths, the T ₁ relaxation time of the soft tissue was 382.75-506.48 ms, and T ₂ was 51.11-70.76 ms. ¹⁸F-FDG tracer could be clearly observed in PET imaging. The 3D printed mouse phantom was successfully constructed with tissue-equivalent materials. Our model can be used for CT, MRI, and PET as a standard device for small-animal imaging and quality control.

Genetics of Alcohol Use Disorder: A Role for Induced Pluripotent Stem Cells?

Alcohol use disorder (AUD) affects millions of people and costs nearly 250 billion dollars annually. Few effective FDA-approved treatments exist, and more are needed. AUDs have a strong heritability, but only a few genes have been identified with a large effect size on disease phenotype. Genomewide association studies (GWASs) have identified common variants with low effect sizes, most of which are in noncoding regions of the genome. Animal models frequently fail to recapitulate key molecular features of neuropsychiatric disease due to the polygenic nature of the disease, partial conservation of coding regions, and significant disparity in noncoding regions. By contrast, human induced pluripotent stem cells (hiPSCs) derived from patients provide a powerful platform for evaluating genes identified by GWAS and modeling complex interactions in the human genome. hiPSCs can be differentiated into a wide variety of human cells, including neurons, glia, and hepatic cells, which are compatible with numerous functional assays and genome editing techniques. In this review, we focus on current applications and future directions of patient hiPSC-derived central nervous system cells for modeling AUDs in addition to highlighting successful applications of hiPSCs in polygenic neuropsychiatric diseases.

Opportunities and challenges in modeling human brain disorders in transgenic primates

Molecular genetic tools have had a profound impact on neuroscience, but until recently their application has largely been confined to a few model species, most notably mouse, zebrafish, Drosophila melanogaster and Caenorhabditis elegans. With the development of new genome engineering technologies such as CRISPR, it is becoming increasingly feasible to apply these molecular tools in a wider range of species, including nonhuman primates. This will lead to many opportunities for brain research, but it will also pose challenges. Here we identify some of these opportunities and challenges in light of recent and foreseeable technological advances and offer some suggestions. Our main focus is on the creation of new primate disease models for understanding the pathological mechanisms of brain disorders and for developing new approaches to effective treatment. However, we also emphasize that primate genetic models have great potential to address many fundamental questions about brain function, providing an essential foundation for future progress in disease research.

Shortcomings of Animal Models and the Rise of Engineered Human Cardiac Tissue

We provide here an historical context of how studies utilizing engineered human cardiac muscle can complement and in some cases substitute animal and cell models for studies of disease and drug testing. We give an overview of the development of animal models and discuss the ability of novel human tissue models to overcome limited predictive power of cell culture and animal models in studies of drug efficacy and safety. The in vitro generation of cardiac tissue is discussed in the context of state of the art in the field. Finally we describe the assembly of multitissue platforms for more accurate representation of integrated human cardiac physiology and consider the advantages of in silico drug trials to augment our ability to predict drug-drug and organ-organ interactions in humans.

Using hiPSCs to model neuropsychiatric copy number variations (CNVs) has potential to reveal underlying disease mechanisms

Schizophrenia is a neuropsychological disorder with a strong heritable component; genetic risk for schizophrenia is conferred by both common variants of relatively small effect and rare variants with high penetrance. Genetically engineered mouse models can recapitulate rare variants, displaying some behavioral defects associated with schizophrenia; however, these mouse models cannot recapitulate the full genetic architecture underlying the disorder. Patient-derived human induced pluripotent stem cells (hiPSCs) present an alternative approach for studying rare variants, in the context of all other risk alleles. Genome editing technologies, such as CRISPR-Cas9, enable the generation of isogenic hiPSC lines with which to examine the functional contribution of single variants within any genetic background. Studies of these rare variants using hiPSCs have the potential to identify commonly disrupted pathways in schizophrenia and allow for the identification of new therapeutic targets. This article is part of a Special Issue entitled SI:StemsCellsinPsychiatry.

Gene length as a biological timer to establish temporal transcriptional regulation

Transcriptional timing is inherently influenced by gene length, thus providing a mechanism for temporal regulation of gene expression. While gene size has been shown to be important for the expression timing of specific genes during early development, whether it plays a role in the timing of other global gene expression programs has not been extensively explored. Here, we investigate the role of gene length during the early transcriptional response of human fibroblasts to serum stimulation. Using the nascent sequencing techniques Bru-seq and BruUV-seq, we identified immediate genome-wide transcriptional changes following serum stimulation that were linked to rapid activation of enhancer elements. We identified 873 significantly induced and 209 significantly repressed genes. Variations in gene size allowed for a large group of genes to be simultaneously activated but produce full-length RNAs at different times. The median length of the group of serum-induced genes was significantly larger than the median length of all expressed genes, housekeeping genes, and serum-repressed genes. These gene length relationships were also observed in corresponding mouse orthologs, suggesting that relative gene size is evolutionarily conserved. The sizes of transcription factor and microRNA genes immediately induced after serum stimulation varied dramatically, setting up a cascade mechanism for temporal expression arising from a single activation event. The retention and expansion of large intronic sequences during evolution have likely played important roles in fine-tuning the temporal expression of target genes in various cellular response programs.

Evaluation of the mirn23a Cluster through an iTRAQ-based Quantitative Proteomic Approach

MicroRNAs (miRNAs) are post-transcriptional regulators of gene expression that are implicated in a number of disease states. MiRNAs can exist as individual entities or may be clustered and transcribed as a single polycistron. The mirn23a cluster consists of three miRNAs: miR-23a, miR-24-2, and miR-27a. Although these miRNAs are transcribed together, they often exist at varying levels in the cell. Despite the fact that the mirn23a cluster is known to play a role in a number of diseases and developmental processes, few direct targets have been identified. In this study, we examined the effects of miR-23a, miR-24-2, miR-27a, or the mirn23a cluster overexpression on the proteome of 70Z/3 pre-B lymphoblast cells. Quantitative mass spectrometry using isobaric tags for relative and absolute quantification (iTRAQ) allowed for the global profiling of cell lines after miRNA overexpression. We identified a number of targets of each miRNA that contained predicted miRNA seed sequences and are likely direct targets. In addition, we discovered a cohort of shared miRNA targets and cluster targets, demonstrating the importance of studying miRNA clusters in their entirety.

Pseudo-Reference-Based Assembly of Vertebrate Transcriptomes

High-throughput RNA sequencing (RNA-seq) provides a comprehensive picture of the transcriptome, including the identity, structure, quantity, and variability of expressed transcripts in cells, through the assembly of sequenced short RNA-seq reads. Although the reference-based approach guarantees the high quality of the resulting transcriptome, this approach is only applicable when the relevant reference genome is present. Here, we developed a pseudo-reference-based assembly (PRA) that reconstructs a transcriptome based on a linear regression function of the optimized mapping parameters and genetic distances of the closest species. Using the linear model, we reconstructed transcriptomes of four different aves, the white leg horn, turkey, duck, and zebra finch, with the Gallus gallus genome as a pseudo-reference, and of three primates, the chimpanzee, gorilla, and macaque, with the human genome as a pseudo-reference. The resulting transcriptomes show that the PRAs outperformed the de novo approach for species with within about 10% mutation rate among orthologous transcriptomes, enough to cover distantly related species as far as chicken and duck. Taken together, we suggest that the PRA method can be used as a tool for reconstructing transcriptome maps of vertebrates whose genomes have not yet been sequenced.

Evolutionary analysis across mammals reveals distinct classes of long non-coding RNAs

BACKGROUND

Recent advances in transcriptome sequencing have enabled the discovery of thousands of long non-coding RNAs (lncRNAs) across many species. Though several lncRNAs have been shown to play important roles in diverse biological processes, the functions and mechanisms of most lncRNAs remain unknown. Two significant obstacles lie between transcriptome sequencing and functional characterization of lncRNAs: identifying truly non-coding genes from de novo reconstructed transcriptomes, and prioritizing the hundreds of resulting putative lncRNAs for downstream experimental interrogation.

RESULTS

We present slncky, a lncRNA discovery tool that produces a high-quality set of lncRNAs from RNA-sequencing data and further uses evolutionary constraint to prioritize lncRNAs that are likely to be functionally important. Our automated filtering pipeline is comparable to manual curation efforts and more sensitive than previously published computational approaches. Furthermore, we developed a sensitive alignment pipeline for aligning lncRNA loci and propose new evolutionary metrics relevant for analyzing sequence and transcript evolution. Our analysis reveals that evolutionary selection acts in several distinct patterns, and uncovers two notable classes of intergenic lncRNAs: one showing strong purifying selection on RNA sequence and another where constraint is restricted to the regulation but not the sequence of the transcript.

CONCLUSION

Our results highlight that lncRNAs are not a homogenous class of molecules but rather a mixture of multiple functional classes with distinct biological mechanism and/or roles. Our novel comparative methods for lncRNAs reveals 233 constrained lncRNAs out of tens of thousands of currently annotated transcripts, which we make available through the slncky Evolution Browser.

Rotational imaging optical coherence tomography for full-body mouse embryonic imaging

Optical coherence tomography (OCT) has been widely used to study mammalian embryonic development with the advantages of high spatial and temporal resolutions and without the need for any contrast enhancement probes. However, the limited imaging depth of traditional OCT might prohibit visualization of the full embryonic body. To overcome this limitation, we have developed a new methodology to enhance the imaging range of OCT in embryonic day (E) 9.5 and 10.5 mouse embryos using rotational imaging. Rotational imaging OCT (RI-OCT) enables full-body imaging of mouse embryos by performing multiangle imaging. A series of postprocessing procedures was performed on each cross-section image, resulting in the final composited image. The results demonstrate that RI-OCT is able to improve the visualization of internal mouse embryo structures as compared to conventional OCT.

MSuPDA: A Memory Efficient Algorithm for Sequence Alignment

Space complexity is a million dollar question in DNA sequence alignments. In this regard, memory saving under pushdown automata can help to reduce the occupied spaces in computer memory. Our proposed process is that anchor seed (AS) will be selected from given data set of nucleotide base pairs for local sequence alignment. Quick splitting techniques will separate the AS from all the DNA genome segments. Selected AS will be placed to pushdown automata's (PDA) input unit. Whole DNA genome segments will be placed into PDA's stack. AS from input unit will be matched with the DNA genome segments from stack of PDA. Match, mismatch and indel of nucleotides will be popped from the stack under the control unit of pushdown automata. During the POP operation on stack, it will free the memory cell occupied by the nucleotide base pair.

Microfluidic Chip-LC/MS-based Glycomic Analysis Revealed Distinct N-glycan Profile of Rat Serum

The rat is an important alternative for studying human pathology owing to certain similarities to humans. Glycomic studies on rat serum have revealed that variations in the N-glycans of glycoproteins correlated with disease progression, which is consistent with the findings in human serum. Therefore, we comprehensively characterized the rat serum N-glycome using microfluidic chip-LC-ESI-QTOF MS and MS/MS techniques. In total, 282 N-glycans, including isomers, were identified. This study is the first to present comprehensive profiling of N-glycans containing O-acetylated sialic acid, among which 27 N-glycans are novel. In addition, the co-existence of N-acetylneuraminic acid (NeuAc) and N-glycolylneuraminic acid (NeuGc) in a single N-glycan ('mixed' N-glycan) was detected and represents a new type of N-glycan in rat serum. The existence of O-acetylated sialic acid is the characteristic feature of rat serum that distinguishes it from mouse and human sera. Comparisons between the rat, mouse, and human serum glycomes revealed that the rat glycome is more similar to that of human sera than to that of mouse sera. Our findings highlight the similarities between the glycomic profile of rat and human sera and provided important selection criteria for choosing an appropriate animal model for pathological and pharmacological studies.

C-It-Loci: a knowledge database for tissue-enriched loci

MOTIVATION

Increasing evidences suggest that most of the genome is transcribed into RNAs, but many of them are not translated into proteins. All those RNAs that do not become proteins are called 'non-coding RNAs (ncRNAs)', which outnumbers protein-coding genes. Interestingly, these ncRNAs are shown to be more tissue specifically expressed than protein-coding genes. Given that tissue-specific expressions of transcripts suggest their importance in the expressed tissue, researchers are conducting biological experiments to elucidate the function of such ncRNAs. Owing greatly to the advancement of next-generation techniques, especially RNA-seq, the amount of high-throughput data are increasing rapidly. However, due to the complexity of the data as well as its high volume, it is not easy to re-analyze such data to extract tissue-specific expressions of ncRNAs from published datasets.

RESULTS

Here, we introduce a new knowledge database called 'C-It-Loci', which allows a user to screen for tissue-specific transcripts across three organisms: human, mouse and zebrafish. C-It-Loci is intuitive and easy to use to identify not only protein-coding genes but also ncRNAs from various tissues. C-It-Loci defines homology through sequence and positional conservation to allow for the extraction of species-conserved loci. C-It-Loci can be used as a starting point for further biological experiments.

AVAILABILITY AND IMPLEMENTATION

C-It-Loci is freely available online without registration at http://c-it-loci.uni-frankfurt.de.

CONTACT

uchida@med.uni-frankfurt.de

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Mixture models for gene expression experiments with two species

Cross-species research in drug development is novel and challenging. A bivariate mixture model utilizing information across two species was proposed to solve the fundamental problem of identifying differentially expressed genes in microarray experiments in order to potentially improve the understanding of translation between preclinical and clinical studies for drug development. The proposed approach models the joint distribution of treatment effects estimated from independent linear models. The mixture model posits up to nine components, four of which include groups in which genes are differentially expressed in both species. A comprehensive simulation to evaluate the model performance and one application on a real world data set, a mouse and human type II diabetes experiment, suggest that the proposed model, though highly structured, can handle various configurations of differential gene expression and is practically useful on identifying differentially expressed genes, especially when the magnitude of differential expression due to different treatment intervention is weak. In the mouse and human application, the proposed mixture model was able to eliminate unimportant genes and identify a list of genes that were differentially expressed in both species and could be potential gene targets for drug development.

Automatic structural parcellation of mouse brain MRI using multi-atlas label fusion

Multi-atlas segmentation propagation has evolved quickly in recent years, becoming a state-of-the-art methodology for automatic parcellation of structural images. However, few studies have applied these methods to preclinical research. In this study, we present a fully automatic framework for mouse brain MRI structural parcellation using multi-atlas segmentation propagation. The framework adopts the similarity and truth estimation for propagated segmentations (STEPS) algorithm, which utilises a locally normalised cross correlation similarity metric for atlas selection and an extended simultaneous truth and performance level estimation (STAPLE) framework for multi-label fusion. The segmentation accuracy of the multi-atlas framework was evaluated using publicly available mouse brain atlas databases with pre-segmented manually labelled anatomical structures as the gold standard, and optimised parameters were obtained for the STEPS algorithm in the label fusion to achieve the best segmentation accuracy. We showed that our multi-atlas framework resulted in significantly higher segmentation accuracy compared to single-atlas based segmentation, as well as to the original STAPLE framework.

Carotenoids gene markers for sweetpotato (Ipomoea batatas L. Lam): applications in genetic mapping, diversity evaluation and cross-species transference

Carotenoids play essential biological roles in plants, and genes involved in the carotenoid biosynthesis pathway are evolutionarily conserved. Orange sweetpotato is an important source of β-carotene, a precursor of vitamin A. In spite of this, only a few research studies have focussed on the molecular aspects of carotenoid genes regarding their specific sequence and structure. In this study, we used published carotenoid gene sequences from Ipomoea and other species for "exon-primed intron-crossing" approaches. Fifteen pairs of primers representing six carotenoid genes were designed for different introns, eleven of which amplified scorable and reproducible alleles. The sequence of PCR products showed high homology to the original ones. Moreover, the structure and sequence of the introns and exons from five carotenoid structural genes were partially defined. Intron length polymorphism and intron single nucleotide polymorphisms were detected in amplified sequences. Marker dosages and allelic segregations were analysed in a mapping population. The developed markers were evaluated in a set of Ipomoeas batatas accessions so as to analyse genetic diversity and conservation applicability. Using CG strategy combined with EPIC-PCR technique, we developed carotenoid gene markers in sweetpotato. We reported the first set of polymorphic Candidate Gene markers for I. batatas, and demonstrated transferability in seven wild Ipomoea species. We described the sequence and structure of carotenoid genes and introduced new information about genomic constitution and allele dosage.

Syntenic global alignment and its application to the gene prediction problem

Given the increasing number of available genomic sequences, one now faces the task of identifying their protein coding regions. The gene prediction problem can be addressed in several ways, and one of the most promising methods makes use of information derived from the comparison of homologous sequences. In this work, we develop a new comparative-based gene prediction program, called Exon_Finder2. This tool is based on a new type of alignment we propose, called syntenic global alignment, that can deal satisfactorily with sequences that share regions with different rates of conservation. In addition to this new type of alignment itself, we also describe a dynamic programming algorithm that computes a best syntenic global alignment of two sequences, as well as its related score. The applicability of our approach was validated by the promising initial results achieved by Exon_Finder2. On a benchmark including 120 pairs of human and mouse genomic sequences, most of their encoded genes were successfully identified by our program.

Heuristic alignment methods

Computation of multiple sequence alignment (MSA) is usually formulated as a combinatory optimization problem of an objective function. Solving the problem for virtually all sensible objective functions is known to be NP-complete implying that some heuristics must be adopted. Several general strategies have been proven effective to obtain accurate MSAs in reasonable computational costs. This chapter is devoted to a brief summary of most successful heuristic approaches.

Comparison of RefSeq protein-coding regions in human and vertebrate genomes

BACKGROUND

Advances in high-throughput sequencing technology have yielded a large number of publicly available vertebrate genomes, many of which are selected for inclusion in NCBI's RefSeq project and subsequently processed by NCBI's eukaryotic annotation pipeline. Genome annotation results are affected by differences in available support evidence and may be impacted by annotation pipeline software changes over time. The RefSeq project has not previously assessed annotation trends across organisms or over time. To address this deficiency, we have developed a comparative protocol which integrates analysis of annotated protein-coding regions across a data set of vertebrate orthologs in genomic sequence coordinates, protein sequences, and protein features.

RESULTS

We assessed an ortholog dataset that includes 34 annotated vertebrate RefSeq genomes including human. We confirm that RefSeq protein-coding gene annotations in mammals exhibit considerable similarity. Over 50% of the orthologous protein-coding genes in 20 organisms are supported at the level of splicing conservation with at least three selected reference genomes. Approximately 7,500 ortholog sets include at least half of the analyzed organisms, show highly similar sequence and conserved splicing, and may serve as a minimal set of mammalian "core proteins" for initial assessment of new mammalian genomes. Additionally, 80% of the proteins analyzed pass a suite of tests to detect proteins that lack splicing conservation and have unusual sequence or domain annotation. We use these tests to define an annotation quality metric that is based directly on the annotated proteins thus operates independently of other quality metrics such as availability of transcripts or assembly quality measures. Results are available on the RefSeq FTP site [http://ftp.ncbi.nlm.nih.gov/refseq/supplemental/ProtCore/SM1.txt].

CONCLUSIONS

Our multi-factored analysis demonstrates a high level of consistency in RefSeq protein representation among vertebrates. We find that the majority of the RefSeq vertebrate proteins for which we have calculated orthology are good as measured by these metrics. The process flow described provides specific information on the scope and degree of conservation for the analyzed protein sequences and annotations and will be used to enrich the quality of RefSeq records by identifying targets for further improvement in the computational annotation pipeline, and by flagging specific genes for manual curation.

Open access to tree genomes: the path to a better forest

An open-access culture and a well-developed comparative-genomics infrastructure must be developed in forest trees to derive the full potential of genome sequencing in this diverse group of plants that are the dominant species in much of the earth's terrestrial ecosystems.

Mobilomics in Saccharomyces cerevisiae strains

BACKGROUND

Mobile Genetic Elements (MGEs) are selfish DNA integrated in the genomes. Their detection is mainly based on consensus-like searches by scanning the investigated genome against the sequence of an already identified MGE. Mobilomics aims at discovering all the MGEs in a genome and understanding their dynamic behavior: The data for this kind of investigation can be provided by comparative genomics of closely related organisms. The amount of data thus involved requires a strong computational effort, which should be alleviated.

RESULTS

Our approach proposes to exploit the high similarity among homologous chromosomes of different strains of the same species, following a progressive comparative genomics philosophy. We introduce a software tool based on our new fast algorithm, called regender, which is able to identify the conserved regions between chromosomes. Our case study is represented by a unique recently available dataset of 39 different strains of S.cerevisiae, which regender is able to compare in few minutes. By exploring the non-conserved regions, where MGEs are mainly retrotransposons called Tys, and marking the candidate Tys based on their length, we are able to locate a priori and automatically all the already known Tys and map all the putative Tys in all the strains. The remaining putative mobile elements (PMEs) emerging from this intra-specific comparison are sharp markers of inter-specific evolution: indeed, many events of non-conservation among different yeast strains correspond to PMEs. A clustering based on the presence/absence of the candidate Tys in the strains suggests an evolutionary interconnection that is very similar to classic phylogenetic trees based on SNPs analysis, even though it is computed without using phylogenetic information.

CONCLUSIONS

The case study indicates that the proposed methodology brings two major advantages: (a) it does not require any template sequence for the wanted MGEs and (b) it can be applied to infer MGEs also for low coverage genomes with unresolved bases, where traditional approaches are largely ineffective.

Unique profile of ordered arrangements of repetitive elements in the C57BL/6J mouse genome implicating their functional roles

The entirety of all protein coding sequences is reported to represent a small fraction (∼2%) of the mouse and human genomes; the vast majority of the rest of the genome is presumed to be repetitive elements (REs). In this study, the C57BL/6J mouse reference genome was subjected to an unbiased RE mining to establish a whole-genome profile of RE occurrence and arrangement. The C57BL/6J mouse genome was fragmented into an initial set of 5,321 units of 0.5 Mb, and surveyed for REs using unbiased self-alignment and dot-matrix protocols. The survey revealed that individual chromosomes had unique profiles of RE arrangement structures, named RE arrays. The RE populations in certain genomic regions were arranged into various forms of complexly organized structures using combinations of direct and/or inverse repeats. Some of these RE arrays spanned stretches of over 2 Mb, which may contribute to the structural configuration of the respective genomic regions. There were substantial differences in RE density among the 21 chromosomes, with chromosome Y being the most densely populated. In addition, the RE array population in the mouse chromosomes X and Y was substantially different from those of the reference human chromosomes. Conversion of the dot-matrix data pertaining to a tandem 13-repeat structure within the Ch7.032 genome unit into a line map of known REs revealed a repeat unit of ∼11.3 Kb as a mosaic of six different RE types. The data obtained from this study allowed for a comprehensive RE profiling, including the establishment of a library of RE arrays, of the reference mouse genome. Some of these RE arrays may participate in a spectrum of normal and disease biology that are specific for mice.

BOAT: automatic alignment of biomedical ontologies using term informativeness and candidate selection

The biomedical sciences is one of the few domains where ontologies are widely being developed to facilitate information retrieval and knowledge sharing, but there still remains the problem that applications using different ontologies cannot share knowledge without explicit references between overlapping concepts. Ontology alignment is the task of identifying such equivalence relations between concepts across ontologies. Its application to the biomedical domain should address two open issues: (1) determining the equivalence of concept-pairs which have overlapping terms in their names, and (2) the high run-time required to align large ontologies which are typical in the biomedical domain. To address them, we present a novel approach, named the Biomedical Ontologies Alignment Technique (BOAT), which is state-of-the-art in terms of F-measure, precision and speed. A key feature of BOAT is that it considers the informativeness of each component word in the concept labels, which has significant impact on biomedical ontologies, resulting in a 12.2% increase in F-measure. Another important feature of BOAT is that it selects for comparison only concept pairs that show high likelihoods of equivalence, based on the similarity of their annotations. BOAT's F-measure of 0.88 for the alignment of the mouse and human anatomy ontologies is on par with that of another state-of-the-art matcher, AgreementMaker, while taking a shorter time.

A mouse model of lamellar intrastromal femtosecond laser keratotomy: ultra-structural, inflammatory, and wound healing responses

PURPOSE

The availability of knockout mouse species provide a highly versatile platform for critically examining the corneal wound healing response. We aimed to develop and characterize the wound healing response in a mouse model of intrastromal femtosecond laser (FSL) keratotomy.

METHODS

An intrastromal lamellar dissection using a Visumax FSL was performed on 16 wild type mice (C57BL6) . The energy level was optimized at 150nJ. The FSL was programmed to perform a lamellar dissection at 50 µM depth without sidecut. The flap was not lifted. Fellow eyes were used as controls. Slit lamp photography and confocal microscopy were performed immediately before the mice were sacrificed 4 h, 1, 3, and 7 days post surgery. Corneas were harvested for immunocytochemistry, transmission electron microscopy (TEM) and light microscopy (LM).

RESULTS

Confocal microscopy showed an absence of keratocytes in the area immediately surrounding the dissection plane. The dissection plane and individual FSL plasma cavitation bubbles were clearly evident on TEM. There was evidence of Keratocyte cell death along the laser resection plane on TEM. LM revealed the dissection plane at a 20 µM depth, although not all epithelial cell layers were intact. Staining for monocytes using antibodies for CD11b (cluster of differentiation 11b) showed early migration at the peripheries at 4 h that increased at 24 h and became more central in treated corneas (p<0.001). Apoptotic cells were evident on TUNEL (terminal deoxynucleotidyl transferase dUTP nick end labeling) assay in the immediate ablation zone and were significantly raised at 4 and 24 h (p<0.001). Ki67 (Kiel 67 protein) positive proliferating keratocytes are evident at 3 days and increased significantly by 7 days (p<0.001). Minimal fibroblast (cluster of differentiation 90, CD90) transformation was seen at 1 week. No myofibroblasts were detected.

DISCUSSION

We have demonstrated that FSL lamellar cuts can be effectively performed on mice and that this model exhibits typical signs of the corneal wound healing response. This model could provide a ubiquitous platform in which to study corneal wound healing responses in both wild type and knockout mice species. The ability to create such a lamellar pocket may be utilizzd for intrastromal drug delivery.

Approaches to Fungal Genome Annotation

Fungal genome annotation is the starting point for analysis of genome content. This generally involves the application of diverse methods to identify features on a genome assembly such as protein-coding and non-coding genes, repeats and transposable elements, and pseudogenes. Here we describe tools and methods leveraged for eukaryotic genome annotation with a focus on the annotation of fungal nuclear and mitochondrial genomes. We highlight the application of the latest technologies and tools to improve the quality of predicted gene sets. The Broad Institute eukaryotic genome annotation pipeline is described as one example of how such methods and tools are integrated into a sequencing center's production genome annotation environment.