Search

Next-generation sequencing (NGS) is a useful molecular diagnostic tool for genetic diseases. However, due to the presence of highly homologous pseudogenes, it is challenging to use short-read NGS for analyzing mutations of the Shwachman-Bodian-Diamond syndrome (SBDS) gene. The SBDS mutation spectrum was analyzed in the Chinese population, which revealed that SBDS variants were primarily from sequence exchange between SBDS and its pseudogene at the base-pair level, predominantly in the coding region and splice junction of exon two. The c.258+2T>C and c.185_184TA>GT variants were the two most common pathogenic SBDS variants in the Chinese population, resulting in a total carrier frequency of 1.19%. When analyzing pathogenic variants in the SBDS gene from the NGS data, the misalignment was identified as a common issue, and there were different probabilities of misalignment for different pathogenic variants. Here, we present a novel mathematical method for identifying pathogenic variants in the SBDS gene from the NGS data, which utilizes read-depth of the paralogous sequence variant (PSV) loci of SBDS and its pseudogene. Combined with PCR and STR orthogonal experiments, SBDS gene mutation analysis results were improved in 40% of clinical samples, and various types of mutations such as homozygous, compound heterozygous, and uniparental diploid were explored. The findings effectively reduce the impact of misalignment in NGS-based SBDS mutation analysis and are helpful for the clinical diagnosis of SBDS-related diseases, the research into population variation, and the carrier screening. [ABSTRACT FROM AUTHOR]

Multi-marker metabarcoding is increasingly being used to generate biodiversity information across different domains of life from microbes to fungi to animals such as for molecular ecology and biomonitoring applications in different sectors from academic research to regulatory agencies and industry. Current popular bioinformatic pipelines support microbial and fungal marker analysis, while ad hoc methods are often used to process animal metabarcode markers from the same study. MetaWorks provides a harmonized processing environment, pipeline, and taxonomic assignment approach for demultiplexed Illumina reads for all biota using a wide range of metabarcoding markers such as 16S, ITS, and COI. A Conda environment is provided to quickly gather most of the programs and dependencies for the pipeline. Several workflows are provided such as: taxonomically assigning exact sequence variants, provides an option to generate operational taxonomic units, and facilitates single-read processing. Pipelines are automated using Snakemake to minimize user intervention and facilitate scalability. All pipelines use the RDP classifier to provide taxonomic assignments with confidence measures. We extend the functionality of the RDP classifier for taxonomically assigning 16S (bacteria), ITS (fungi), and 28S (fungi), to also support COI (eukaryotes), rbcL (eukaryotes, land plants, diatoms), 12S (fish, vertebrates), 18S (eukaryotes, diatoms) and ITS (fungi, plants). MetaWorks properly handles ITS by trimming flanking conserved rRNA gene regions as well as protein coding genes by providing two options for removing obvious pseudogenes. MetaWorks can be downloaded from https://github.com/terrimporter/MetaWorks and quickstart instructions, pipeline details, and a tutorial for new users can be found at https://terrimporter.github.io/MetaWorksSite. [ABSTRACT FROM AUTHOR]

BackgroundPseudogenes are non-functional sequences in the genome with homologous sequences that are functional (i.e. genes). They are abundant in eukaryotes where they have been extensively investigated, while in prokaryotes they are significantly scarcer and less well studied. Here we conduct a comprehensive analysis of the evolution of orthologs of Mycobacterium leprae pseudogenes in prokaryotes. The leprosy pathogen M. leprae is of particular interest since it contains an unusually large number of pseudogenes, comprising approximately 40% of its entire genome. The analysis is conducted in both broad and narrow phylogenetic ranges.ResultsWe have developed an informatics-based approach to characterize the evolution of pseudogenes. This approach combines tools from phylogenomics, genomics, and transcriptomics. The results we obtain are used to assess the contributions of two mechanisms for pseudogene formation: failed horizontal gene transfer events and disruption of native genes.ConclusionsWe conclude that, although it was reported that in most bacteria the former is most likely responsible for the majority of pseudogenization events, in mycobacteria, and in particular in M. leprae with its exceptionally high pseudogene numbers, the latter predominates. We believe that our study sheds new light on the evolution of pseudogenes in bacteria, by utilizing new methodologies that are applied to the unusually abundant M. leprae pseudogenes and their orthologs.

Transfer RNAs (tRNAs) are intermediate-sized non-coding RNAs found in all organisms that help translate messenger RNA into protein. Recently, the number of sequenced plant genomes has increased dramatically. The availability of this extensive data greatly accelerates the study of tRNAs on a large scale. Here, 8,768,261 scaffolds/chromosomes containing 229,093 giga-base pairs representing whole-genome sequences of 256 plant species were analyzed to identify tRNA genes. As a result, 331,242 nuclear, 3,216 chloroplast, and 1,467 mitochondrial tRNA genes were identified. The nuclear tRNA genes include 275,134 tRNAs decoding 20 standard amino acids, 1,325 suppressor tRNAs, 6,273 tRNAs with unknown isotypes, 48,475 predicted pseudogenes, and 37,873 tRNAs with introns. Efforts also extended to the creation of PltRNAdb (https://bioinformatics.um6p.ma/PltRNAdb/index.php), a data source for tRNA genes from 256 plant species. PltRNAdb website allows researchers to search, browse, visualize, BLAST, and download predicted tRNA genes. PltRNAdb will help improve our understanding of plant tRNAs and open the door to discovering the unknown regulatory roles of tRNAs in plant genomes. [ABSTRACT FROM AUTHOR]

The wolverine (Gulo gulo) in Finland has undergone significant population declines in the past. Since major histocompatibility complex (MHC) genes encode proteins involved in pathogen recognition, the diversity of these genes provides insights into the immunological fitness of regional populations. We sequenced 862 amplicons (242 bp) of MHC class II DRB exon 2 from 32 Finnish wolverines and identified 11 functional alleles and three pseudogenes. A molecular phylogenetic analysis indicated trans-species polymorphism, and PAML and MEME analyses indicated positive selection, suggesting that the Finnish wolverine DRB genes have evolved under balancing and positive selection. In contrast to DRB gene analyses in other species, allele frequencies in the Finnish wolverines clearly indicated the existence of two regional subpopulations, congruent with previous studies based on neutral genetic markers. In the Finnish wolverine, rapid population declines in the past have promoted genetic drift, resulting in a lower genetic diversity of DRB loci, including fewer alleles and positively selected sites, than other mustelid species analyzed previously. Our data suggest that the MHC region in the Finnish wolverine population was likely affected by a recent bottleneck. [ABSTRACT FROM AUTHOR]

Here we present the results of a large-scale bioinformatics annotation of non-coding RNA loci in 48 avian genomes. Our approach uses probabilistic models of hand-curated families from the Rfam database to infer conserved RNA families within each avian genome. We supplement these annotations with predictions from the tRNA annotation tool, tRNAscan-SE and microRNAs from miRBase. We identify 34 lncRNA-associated loci that are conserved between birds and mammals and validate 12 of these in chicken. We report several intriguing cases where a reported mammalian lncRNA, but not its function, is conserved. We also demonstrate extensive conservation of classical ncRNAs (e.g., tRNAs) and more recently discovered ncRNAs (e.g., snoRNAs and miRNAs) in birds. Furthermore, we describe numerous "losses" of several RNA families, and attribute these to either genuine loss, divergence or missing data. In particular, we show that many of these losses are due to the challenges associated with assembling avian microchromosomes. These combined results illustrate the utility of applying homology-based methods for annotating novel vertebrate genomes.

The seminal vesicle secretions of guinea pig and rabbit were analyzed for semen coagulum proteins. Using SDS-PAGE we discovered a previously not fully recognized semen coagulum protein, Svp5, in the guinea pig and a single predominant component, SVP200, in the rabbit. Potential genes of these proteins were identified in genome databases by their homology with human and murine genes. The structure of their fullength transcripts was determined using seminal vesicle cDNA and sequencing primers based on genomic sequences. Homology searching indicated that both Svp5 and SVP200 were synthesized from composite genes that were the result of merger between two genes showing homology with human SEMG2 and PI3. For a deeper understanding of the evolution of the genes, we retrieved and analyzed genome sequences from the REST gene loci, encompassing genes of semen coagulum proteins and related rapidly evolving seminal vesicle-transcribed genes, of 14 rodents and 2 lagomorphs. The analysis showed that rodents of the suborders myomorpha, hystricomorpha, and castorimorpha had unique sets of REST genes, whereas sciuromorpha seemed to be lacking such genes. It also indicated a closer relationship between myomorpha and castorimorpha than to rodents of the two other analyzed suborders. In the lagomorph species, the pika appeared to be devoid of REST genes, whereas the rabbit had a single expressed REST gene, SVP200, and two pseudogenes. The structural similarity of semen coagulum proteins in rabbit and hystricomph species suggests that they are closely related. This was also supported by other similarities at their REST gene loci, e.g. the finding of a PI3-like gene in the rabbit that also had features in common with caltrin2 of hystricomorph rodents. The homologies indicate that hystricomorpha may have separated from myomorpha and castorimorpha before the separation of hystricomorpha from lagomorpha. [ABSTRACT FROM AUTHOR]

ZMAT2 is among the least-studied of mammalian proteins and genes, even though it is the ortholog of Snu23, a protein involved in pre-mRNA splicing in yeast. Here we have used data from genomic and gene expression repositories to examine the Zmat2 gene and locus in 8 terrestrial vertebrates, 10 ray-finned fish, and 1 lobe-finned fish representing > 500 million years of evolutionary diversification. The analyses revealed that vertebrate Zmat2 genes are similar to their mammalian counterparts, as in 16/19 species studied they contain 6 exons, and in 18/19 encode a single conserved protein. However, unlike in mammals, no Zmat2 pseudogenes were identified in these vertebrates, although an expressed Zmat2 paralog was characterized in flycatcher that resembled a DNA copy of a processed and retro-transposed mRNA, and thus could be a proto-pseudogene captured during its evolutionary journey from active to inert. The Zmat2 locus in terrestrial vertebrates, and in spotted gar and coelacanth, also shares additional genes with its mammalian counterparts, including Histidyl-tRNA synthetase (Hars), Hars2, and others, but these are absent from the Zmat2 locus in teleost fish, in which Stem-loop-binding protein 2 (Slbp2) and Lymphocyte cytosolic protein 2a (Lcp2a) are present instead. Taken together, these observations argue that a recognizable Zmat2 was present in the earliest vertebrate ancestors, and postulate that during chromosomal tetraploidization and subsequent re-diploidization during modern teleost evolution, the duplicated Zmat2 gene was retained and the original lost. This study also highlights how information from genomic resources can be leveraged to reveal new biologically significant insights. [ABSTRACT FROM AUTHOR]

For rodents, olfaction is essential for locating food, recognizing mates and competitors, avoiding predators, and navigating their environment. It is thought that rodents may have expanded olfactory receptor repertoires in order to specialize in olfactory behavior. Despite being the largest clade of mammals and depending on olfaction relatively little work has documented olfactory repertoires outside of conventional laboratory species. Here we report the olfactory receptor repertoire of the African giant pouched rat (Cricetomys ansorgei), a Muroid rodent distantly related to mice and rats. The African giant pouched rat is notable for its large cortex and olfactory bulbs relative to its body size compared to other sympatric rodents, which suggests anatomical elaboration of olfactory capabilities. We hypothesized that in addition to anatomical elaboration for olfaction, these pouched rats might also have an expanded olfactory receptor repertoire to enable their olfactory behavior. We examined the composition of the olfactory receptor repertoire to better understand how their sensory capabilities have evolved. We identified 1145 intact olfactory genes, and 260 additional pseudogenes within 301 subfamilies from the African giant pouched rat genome. This repertoire is similar to mice and rats in terms of size, pseudogene percentage and number of subfamilies. Analyses of olfactory receptor gene trees revealed that the pouched rat has 6 expansions in different subfamilies compared to mice, rats and squirrels. We identified 81 orthologous genes conserved among 4 rodent species and an additional 147 conserved genes within the Muroid rodents. The orthologous genes shared within Muroidea suggests that there may be a conserved Muroid-specific olfactory receptor repertoire. We also note that the description of this repertoire can serve as a complement to other studies of rodent olfaction, as the pouched rat is an outgroup within Muroidea. Thus, our data suggest that African giant pouched rats are capable of both natural and trained olfactory behaviors with a typical Muriod olfactory receptor repertoire. [ABSTRACT FROM AUTHOR]

Pathways leading to formation of non-coding RNA and protein genes are varied and complex. We report finding a conserved repeat sequence present in human and chimpanzee genomes that appears to have originated from a common primate ancestor. This sequence is repeatedly copied in human chromosome 22 (chr22) low copy repeats (LCR22) or segmental duplications and forms twenty-one different genes, which include the human long intergenic non-coding RNA (lincRNA) family FAM230, a newly discovered lincRNA gene family termed conserved long intergenic non-coding RNAs (clincRNA), pseudogene families, as well as the gamma-glutamyltransferase (GGT) protein gene family and the RNA pseudogenes that originate from GGT sequences. Of particular interest are the GGT5 and USP18 protein genes that appear to have formed from an homologous repeat sequence that also forms the clincRNA gene family. The data point to ancestral DNA sequences, conserved through evolution and duplicated in humans by chromosomal repeat sequences that may serve as functional genomic elements in the development of diverse genes. [ABSTRACT FROM AUTHOR]

Anaplasma phagocytophilum is a tick-transmitted bacterial pathogen of humans and other animals, and is an obligate intracellular parasite. Throughout the course of infection, hosts acquire temporary resistance to granulocytic anaplasmosis as they develop immunity specific for the major antigen, major surface protein 2 (Msp2). However, the bacterium then utilizes a novel recombination mechanism shuffling functional pseudogenes sequentially into an expression cassette with conserved 5' and 3' ends, bypassing host immunity. Approximately 100 pseudogenes are present in the only fully sequenced human-origin HZ genome, representing the possibility for almost unlimited antigenic diversity. In the present study, we identified a select group of 20% of the A. phagocytophilum HZ msp2 pseudogenes that have matched preferentially to human, canine, and equine expression cassettes. Pseudogenes cluster predominantly in one spatial run limited to a single genomic island in less than 50% of the genome but phylogenetically related pseudogenes are neither necessarily located in close proximity on the genome nor share similar percent identity with expression cassettes. Pseudogenes near the expression cassette (and the origin) are more likely to be expressed than those farther away. Taken together, these findings suggest that there may be natural selection pressure to retain pseudogenes in one cluster near the putative origin of replication, even though global recombination shuffles pseudogenes around the genome, separating pseudogenes that share genetic origins as well as those with similar identities.

Chloroplasts significantly influence species phylogenies because of their maternal inheritance and the moderate evolutionary rate of their genomes. Avocado, which is a member of the family Lauraceae, has received considerable attention from botanists, likely because of its position as a basal angiosperm. However, there is relatively little avocado genomic information currently available. In this study, six complete avocado chloroplast genomes from three ecological races were assembled to examine the sequence diversity among the three avocado ecological races. A comparative genomic analysis revealed that 515 simple sequence repeat loci and 176 repeats belonging to four other types were polymorphic across the six chloroplast genomes. Three highly variable regions (trnC-GCA-petN, petN-psbM, and petA-psbJ) were identified as highly informative markers. A phylogenetic analysis based on 79 common protein-coding genes indicated that the six examined avocado accessions from three ecological races form a monophyletic clade. The other three genera belonging to the Persea group clustered to form a sister clade with a high bootstrap value. These chloroplast genomes provide important genetic information for future attempts at identifying avocado races and for the related biological research. [ABSTRACT FROM AUTHOR]

Campylobacter transducer-like proteins (Tlps), also known as methyl-accepting chemotaxis proteins (MCPs), are associated with virulence as well as niche and host adaptation. While functional attributes of these proteins are being elucidated, little has been published regarding their sequence diversity or chromosomal locations and context, although they appear to define invertible regions within Campylobacter jejuni genomes. Genome assemblies for several species of Campylobacter were obtained from the publicly available NCBI repositories. Genomes from all isolates were obtained from GenBank and assessed for Tlp content, while data from isolates with complete, finished genomes were used to determine the identity of Tlps as well as the gene content of putative invertible elements (IEs) in C. jejuni (Cj) and C. coli (Cc). Tlps from several Campylobacter species were organized into a nomenclature system and novel Tlps were defined and named for Cj and Cc. The content of Tlps appears to be species-specific, though diverse within species. Cj and Cc carried overlapping, related Tlp content, as did the three C. fetus subspecies. Tlp1 was detected in 88% of Cj isolates and approximately 43% of Cc, and was found in a different conserved chromosomal location and genetic context in each species. Tlp1 and Tlp 3 predominated in genomes from Cj whereas other Tlps were detected less frequently. Tlp13 and Tlp20 predominated in genomes from Cc while some Cj/Cc Tlps were not detected at all. Tlps 2–4 and 11–20 were less frequently detected and many showed sequence heterogeneity that could affect substrate binding, signal transduction, or both. Tlps other than Tlp1, 7, and 10 had substantial sequence identity in the C-terminal half of the protein, creating chromosomal repeats potentially capable of mediating the inversion of large chromosomal DNA. Cj and Cc Tlps were both found in association with only 14 different genes, indicating a limited genomic context. In Cj these Tlps defined IEs that were for the most part found at a single chromosomal location and comprised of a conserved set of genes. Cc IEs were situated at very different chromosomal locations, had different structures than Cj IEs, and were occasionally incomplete, therefore not capable of inversion. Tlps may have a role in Campylobacter genome structure and dynamics as well as acting as chemoreceptors mediating chemotactic responses. [ABSTRACT FROM AUTHOR]

Human Multipotent Stromal Cells (MSCs) are a valuable resource for regenerative medicine and are widely studied. They can be isolated from a variety of tissues and differentiate into multiple cell types (multi-potent). Many reports have been published using human MSCs and to be able to compare outcome, or be able to identify differences between MSCs, several cell surface markers have been proposed. Nevertheless, still many differences remain. Gene expression is known to be different between cell stage and origin. Furthermore, cells cultured on a culture dish (2D) show different gene expression profiles as compared to cells grown on scaffolds (3D). Even the RNA extraction method and the selection of genes used for normalisation have a role in gene expression profiling. To be able to compare gene expression data from samples cultured in different dimensions and RNA extracted using a variety of protocols we set out to define a set of reference genes suitable to normalise qPCR data from a very heterogeneous sample set. Hereto, Trizol was used to extract RNA from human MSCs cultured in 3D and 2D to validate newly designed and previously published primer sets. Subsequently, RNA from fresh human MSC samples and samples stored in RLT-buffer, Trizol or RNAlater was extracted using RNeasy and Trizol methods. All samples have been used to rank the candidate reference genes according to their stability after qPCR enabling identification of the most suitable reference gene(s) for normalisation of a heterogeneous sample set. The most stably expressed reference genes indicated superior normalisation of MSC marker gene expression over the least stable reference genes. [ABSTRACT FROM AUTHOR]

Glutathione S-Transferases (GSTs) are phase II detoxification enzymes that may have evolved in response to changes of environmental substrates. GST genes formed a multigene family and in mammals, there are six classes known as Alpha, Mu, Omega, Pi, Theta, and Zeta. Recent studies in phase I detoxification system specifically the cytochrome P450s provided a general explanation on why genes from a common origin such as those in a multigene family have both phylogenetically stable and unstable genes. Genes that participate in core functions of organisms such as development and physiology are stable whereas genes that play a role in detoxification are unstable and evolve in a process known as birth-death evolution, which is characterised by frequent gene gains and losses. The generality of the birth-death model at explaining the evolution of detoxification enzymes beyond the phase I enzyme has not been comprehensively explored. This work utilized 383 Gst genes and 300 pseudogenes across 22 mammalian species to study gene gains and losses. GSTs vary greatly in their phylogenetic stability despite their overall sequence similarity. Stable Gst genes from Omega and Zeta classes do not show fluctuation in gene numbers from human to opossum. These genes play a role in biosynthesis related functions. Unstable genes that include Alpha, Mu, Pi and Theta undergo frequent gene gain and loss in a process known as birth-death evolution. Gene members of these four classes are well known for their roles in detoxification. Our positive selection screen identified five positively selected sites in mouse GSTA3. Previous studies showed two of these sites (108H and 208E) were biochemically tested as important residues that conferred catalytic activity against the toxic aflatoxin B1-8,9-epoxide. The functional significance against aflatoxin of the remaining three positively selected sites warrant further investigation. [ABSTRACT FROM AUTHOR]

As part of large on-going vaccine impact studies in Fiji and Mongolia, we identified 25/2750 (0.9%) of nasopharyngeal swabs by microarray that were positive for Streptococcus pneumoniae contained pneumococci with a divergent 33F capsular polysaccharide locus (designated ‘33F-1’). We investigated the 33F-1 capsular polysaccharide locus to better understand the genetic variation and its potential impact on serotyping results. Whole genome sequencing was conducted on ten 33F-1 pneumococcal isolates. Initially, sequence reads were used for molecular serotyping by PneumoCaT. Phenotypic typing of 33F-1 isolates was then performed using the Quellung reaction and latex agglutination. Genome assemblies were used in phylogenetic analyses of each gene in the capsular locus to investigate genetic divergence. All ten pneumococcal isolates with the 33F-1 cps locus typed as 33F by Quellung and latex agglutination. Unlike the reference 33F capsule locus sequence, DNA microarray and PneumoCaT analyses found that 33F-1 pneumococci lack the wcjE gene, and instead contain wcyO with a frameshift mutation. Phylogenetic analyses found the wzg, wzh, wzd, wze, wchA, wciG and glf genes in the 33F-1 cps locus had higher DNA sequence similarity to homologues from other serotypes than to the 33F reference sequence. We have discovered a novel genetic variant of serotype 33F, which lacks wcjE and contains a wcyO pseudogene. This finding adds to the understanding of molecular epidemiology of pneumococcal serotype diversity, which is poorly understood in low and middle-income countries. [ABSTRACT FROM AUTHOR]

Introduction: Pseudogenes are paralogues of functional genes historically viewed as defunct due to either the lack of regulatory elements or the presence of frameshift mutations. Recent evidence, however, suggests that pseudogenes may regulate gene expression, although the functional role of pseudogenes remains largely unknown. We previously reported that MYLKP1, the pseudogene of MYLK that encodes myosin light chain kinase (MLCK), is highly expressed in lung and colon cancer cell lines and tissues but not in normal lung or colon. The MYLKP1 promoter is minimally active in normal bronchial epithelial cells but highly active in lung adenocarcinoma cells. In this study, we further validate MYLKP1 as an oncogene via elucidation of the functional role of MYLKP1 genetic variants in colon cancer risk. Methods: Proliferation and migration assays were performed in MYLKP1-transfected colon and lung cancer cell lines (H441, A549) and commercially-available normal lung and colon cells. Fourteen MYLKP1 SNPs (MAFs >0.01) residing within the 4 kb MYLKP1 promoter region, the core 1.4 kb of MYLKP1 gene, and a 4 kb enhancer region were selected and genotyped in a colorectal cancer cohort. MYLKP1 SNP influences on activity of MYLKP1 promoter (2kb) was assessed by dual luciferase reporter assay. Results: Cancer cell lines, H441 and A549, exhibited increased MYLKP1 expression, increased MYLKP1 luciferase promoter activity, increased proliferation and migration. Genotyping studies identified two MYLKP1 SNPs (rs12490683; rs12497343) that significantly increase risk of colon cancer in African Americans compared to African American controls. Rs12490683 and rs12497343 further increase MYLKP1 promoter activity compared to the wild type MYLKP1 promoter. Conclusion: MYLKP1 is a cancer-promoting pseudogene whose genetic variants differentially enhance cancer risk in African American populations. [ABSTRACT FROM AUTHOR]

Santalales is a large order, with over 2200 species, most of which are root or aerial (stem) hemiparasites. In this study, we report the newly assembled chloroplast genome of Dendrotrophe varians (140,666 bp) in the family Amphorogynaceae and the cp genomes of Helixanthera parasitica (124,881 bp) and Macrosolen cochinchinensis (122,986 bp), both in the family Loranthaceae. We compared the cp genomes of 11 Santalales including eight currently available cp genomes. Santalales cp genomes are slightly or not reduced in size (119–147 kb), similar to other hemiparasitic species, when compared with typical angiosperm cp genomes (120–170 kb). In a phylogeny examining gene content, the NADH dehydrogenase gene group is the only one among eight functional gene groups that lost complete functionally in all examined Santalales. This supports the idea that the functional loss of ndh genes is the initial stage in the evolution of the plastome of parasitic plants, but the loss has occurred independently multiple times in angiosperms, while they are not found in some parasites. This suggests that the functional loss of ndh genes is not essential for the transition from autotroph to parasite. We additionally examined the correlation between gene content and type of parasitism (obligate/facultative and stem/root parasites) of all hemiparasitic species in which cp genomes have been reported to date. Correlation was not found in any types of parasitism. [ABSTRACT FROM AUTHOR]

Cytochromes P450 are enzymes that participate in a wide range of functions in plants, from hormonal signaling and biosynthesis of structural polymers, to defense or communication with other organisms. They represent one of the largest gene/protein families in the plant kingdom. The manual annotation of cytochrome P450 genes in the genome of Vitis vinifera PN40024 revealed 579 P450 sequences, including 279 complete genes. Most of the P450 sequences in grapevine genome are organized in physical clusters, resulting from tandem or segmental duplications. Although most of these clusters are small (2 to 35, median = 3), some P450 families, such as CYP76 and CYP82, underwent multiple duplications and form large clusters of homologous sequences. Analysis of gene expression revealed highly specific expression patterns, which are often the same within the genes in large physical clusters. Some of these genes are induced upon biotic stress, which points to their role in plant defense, whereas others are specifically activated during grape berry ripening and might be responsible for the production of berry-specific metabolites, such as aroma compounds. Our work provides an exhaustive and robust annotation including clear identification, structural organization, evolutionary dynamics and expression patterns for the grapevine cytochrome P450 families, paving the way to efficient functional characterization of genes involved in grapevine defense pathways and aroma biosynthesis. [ABSTRACT FROM AUTHOR]

NANOG has been demonstrated to play an essential role in the maintenance of embryonic stem cells, and its pseudogene, NANOGP8, is suggested to promote the cancer stem cell phenotype. As the roles of these genes are intimately involved with glioblastoma multiforme progression and exosomes are critical in intercellular communication, we conducted a detailed analysis of the association of the NANOG gene family with exosomes to identify diagnostic markers for cancer. Exosomes were precipitated from conditioned culture media from various cell lines, and NANOG gene fragments were directly amplified without DNA isolation using multiple primer sets. The use of the enzymes AlwNI and SmaI with restriction fragment length polymorphism analysis functioned to distinguish NANOGP8 from other NANOG family members. Collectively, results suggest that the NANOG DNA associated with exosomes is not full length and that mixed populations of the NANOG gene family exist. Furthermore, sequence analysis of exosomal DNA amplified with a NANOGP8 specific primer set frequently showed an insertion of a 22 bp sequence into the 3’ UTR. The occurrence rate of this insertion was significantly higher in exosomal DNA clones from cancer cells as compared to normal cells. We have detected mixed populations of NANOG DNA associated with exosomes and have identified preferential modulations in the sequences from cancer samples. Our findings, coupled with the properties of exosomes, may allow for the detection of traditionally inaccessible cancers (i.e. GBM) through minimally invasive techniques. Further analysis of exosomal DNA sequences of NANOG and other embryonic stemness genes (OCT3/4, SOX2, etc.) may establish a robust collection of exosome based diagnostic markers, and further elucidate the mechanisms of cancer formation, progression, and metastasis. [ABSTRACT FROM AUTHOR]

DNA methylation is mediated by a conserved family of DNA methyltransferases (Dnmts). The human genome encodes three active Dnmts (Dnmt1, Dnmt3a and Dnmt3b), the tRNA methyltransferase Dnmt2, and the regulatory protein Dnmt3L. Despite their high degree of conservation among different species, genes encoding Dnmts have been duplicated and/or lost in multiple lineages throughout evolution, indicating that the DNA methylation machinery has some potential to undergo evolutionary change. However, little is known about the extent to which this machinery, or the methylome, varies among vertebrates. Here, we study the molecular evolution of Dnmt1, the enzyme responsible for maintenance of DNA methylation patterns after replication, in 79 vertebrate species. Our analyses show that all studied species exhibit a single copy of the DNMT1 gene, with the exception of tilapia and marsupials (tammar wallaby, koala, Tasmanian devil and opossum), each of which displays two apparently functional DNMT1 copies. Our phylogenetic analyses indicate that DNMT1 duplicated before the radiation of major marsupial groups (i.e., at least ~75 million years ago), thus giving rise to two DNMT1 copies in marsupials (copy 1 and copy 2). In the opossum lineage, copy 2 was lost, and copy 1 recently duplicated again, generating three DNMT1 copies: two putatively functional genes (copy 1a and 1b) and one pseudogene (copy 1ψ). Both marsupial copies (DNMT1 copies 1 and 2) are under purifying selection, and copy 2 exhibits elevated rates of evolution and signatures of positive selection, suggesting a scenario of neofunctionalization. This gene duplication might have resulted in modifications in marsupial methylomes and their dynamics. [ABSTRACT FROM AUTHOR]

Background: Radiation-therapy (RT) induces mucositis, a clinically challenging condition with limited prophylactic interventions and no predictive tests. In this pilot study, we applied global gene-expression analysis on serial human oral mucosa tissue and blood cells from patients with tonsil squamous cell cancer (TSCC) to identify genes involved in mucositis pathogenesis. Methods and findings: Eight patients with TSCC each provided consecutive buccal biopsies and blood cells before, after 7 days of RT treatment, and 20 days following RT. We monitored clinical mucositis and performed gene-expression analysis on tissue samples. We obtained control tissue from nine healthy individuals. After RT, expression was upregulated in apoptosis inducer and inhibitor genes, EDA2R and MDM2, and in POLH, a DNA-repair polymerase. Expression was downregulated in six members of the histone cluster family, e.g., HIST1H3B. Gene expression related to proliferation and differentiation was altered, including MKI67 (downregulated), which encodes the Ki-67-proliferation marker, and KRT16 (upregulated), which encodes keratin16. These alterations were not associated with the clinical mucositis grade. However, the expression of LY6G6C, which encodes a surface immunoregulatory protein, was upregulated before treatment in three cases of clinical none/mild mucositis, but not in four cases of ulcerative mucositis. Conclusion: RT caused molecular changes related to apoptosis, DNA-damage, DNA-repair, and proliferation without a correlation to the severity of clinical mucositis. LY6G6C may be a potential protective biomarker for ulcerative mucositis. Based on these results, our study model of consecutive human biopsies will be useful in designing a prospective clinical validation trial to characterize molecular mucositis and identify predictive biomarkers. [ABSTRACT FROM AUTHOR]

In the carnivorous plant family Lentibulariaceae, all three genome compartments (nuclear, chloroplast, and mitochondria) have some of the highest rates of nucleotide substitutions across angiosperms. While the genera Genlisea and Utricularia have the smallest known flowering plant nuclear genomes, the chloroplast genomes (cpDNA) are mostly structurally conserved except for deletion and/or pseudogenization of the NAD(P)H-dehydrogenase complex (ndh) genes known to be involved in stress conditions of low light or CO2 concentrations. In order to determine how the cpDNA are changing, and to better understand the evolutionary history within the Genlisea genus, we sequenced, assembled and analyzed complete cpDNA from six species (G. aurea, G. filiformis, G. pygmaea, G. repens, G. tuberosa and G. violacea) together with the publicly available G. margaretae cpDNA. In general, the cpDNA structure among the analyzed Genlisea species is highly similar. However, we found that the plastidial ndh genes underwent a progressive process of degradation similar to the other terrestrial Lentibulariaceae cpDNA analyzed to date, but in contrast to the aquatic species. Contrary to current thinking that the terrestrial environment is a more stressful environment and thus requiring the ndh genes, we provide evidence that in the Lentibulariaceae the terrestrial forms have progressive loss while the aquatic forms have the eleven plastidial ndh genes intact. Therefore, the Lentibulariaceae system provides an important opportunity to understand the evolutionary forces that govern the transition to an aquatic environment and may provide insight into how plants manage water stress at a genome scale. [ABSTRACT FROM AUTHOR]

Background: Multiorgan dysfunction syndrome contributes to adverse outcomes in advanced heart failure (AdHF) patients after mechanical circulatory support (MCS) implantation and is associated with aberrant leukocyte activity. We tested the hypothesis that preoperative peripheral blood mononuclear cell (PBMC) gene expression profiles (GEP) can predict early postoperative improvement or non-improvement in patients undergoing MCS implantation. We believe this information may be useful in developing prognostic biomarkers. Methods & design: We conducted a study with 29 patients undergoing MCS-surgery in a tertiary academic medical center from 2012 to 2014. PBMC samples were collected one day before surgery (day -1). Clinical data was collected on day -1 and day 8 postoperatively. Patients were classified by Sequential Organ Failure Assessment score and Model of End-stage Liver Disease Except INR score (measured eight days after surgery): Group I = improving (both scores improved from day -1 to day 8, n = 17) and Group II = not improving (either one or both scores did not improve from day -1 to day 8, n = 12). RNA-sequencing was performed on purified mRNA and analyzed using Next Generation Sequencing Strand. Differentially expressed genes (DEGs) were identified by Mann-Whitney test with Benjamini-Hochberg correction. Preoperative DEGs were used to construct a support vector machine algorithm to predict Group I vs. Group II membership. Results: Out of 28 MCS-surgery patients alive 8 days postoperatively, one-year survival was 88% in Group I and 27% in Group II. We identified 28 preoperative DEGs between Group I and II, with an average 93% prediction accuracy. Out of 105 DEGs identified preoperatively between year 1 survivors and non-survivors, 12 genes overlapped with the 28 predictive genes. Conclusions: In AdHF patients following MCS implantation, preoperative PBMC-GEP predicts early changes in organ function scores and correlates with long-term outcomes. Therefore, gene expression lends itself to outcome prediction and warrants further studies in larger longitudinal cohorts. [ABSTRACT FROM AUTHOR]

In this paper, we compare ndh genes in the plastid genome of many Cymbidium species and three closely related taxa in Orchidaceae looking for evidence of ndh gene degradation. Among the 11 ndh genes, there were frequently large deletions in directly repeated or AT-rich regions. Variation in these degraded ndh genes occurs between individual plants, apparently at population levels in these Cymbidium species. It is likely that ndh gene transfers from the plastome to mitochondrial genome (chondriome) occurred independently in Orchidaceae and that ndh genes in the chondriome were also relatively recently transferred between distantly related species in Orchidaceae. Four variants of the ycf1-rpl32 region, which normally includes the ndhF genes in the plastome, were identified, and some Cymbidium species contained at least two copies of that region in their organellar genomes. The four ycf1-rpl32 variants seem to have a clear pattern of close relationships. Patterns of ndh degradation between closely related taxa and translocation of ndh genes to the chondriome in Cymbidium suggest that there have been multiple bidirectional intracellular gene transfers between two organellar genomes, which have produced different levels of ndh gene degradation among even closely related species. [ABSTRACT FROM AUTHOR]

Keratins represent a large protein family with essential structural and functional roles in epithelial cells of skin, hair follicles, and other organs. During evolution the genes encoding keratins have undergone multiple rounds of duplication and humans have two clusters with a total of 55 functional keratin genes in their genomes. Due to the high similarity between different keratin paralogs and species-specific differences in gene content, the currently available keratin gene annotation in species with draft genome assemblies such as dog and horse is still imperfect. We compared the National Center for Biotechnology Information (NCBI) (dog annotation release 103, horse annotation release 101) and Ensembl (release 87) gene predictions for the canine and equine keratin gene clusters to RNA-seq data that were generated from adult skin of five dogs and two horses and from adult hair follicle tissue of one dog. Taking into consideration the knowledge on the conserved exon/intron structure of keratin genes, we annotated 61 putatively functional keratin genes in both the dog and horse, respectively. Subsequently, curators in the RefSeq group at NCBI reviewed their annotation of keratin genes in the dog and horse genomes (Annotation Release 104 and Annotation Release 102, respectively) and updated annotation and gene nomenclature of several keratin genes. The updates are now available in the NCBI Gene database (). [ABSTRACT FROM AUTHOR]

We report the complete chloroplast genomes of three Adenophora species, and analyzed these compared them to five published Campanuloid plastomes. The total genome length of Adenophora divaricata, Adenophora erecta, and Adenophora stricta ranged from 159,759 to 176,331 bp. Among the eight Campanuloid species, many inversions were found to be only in the LSC region. IR contraction was also identified in the plastid genome of Adenophora stricta. Phylogenetic analyses based on 76 protein coding genes showed that Campanuloids are monophyletic, and are composed of two major groups: Campanula s. str. and Rapunculus. When we compared each homologous locus among the four Adenophora species, ten regions showed high nucleotide divergence value (>0.03). Among these, nine loci, excepting ycf3-rpoB, are considered to be useful molecular markers for phylogenetic studies and will be helpful to resolve phylogenetic relationships of Adenophora. [ABSTRACT FROM AUTHOR]

By searching a clinical database of over 60,000 individuals referred for array-based CNV analyses and online resources, we identified four males from three families with intellectual disability, developmental delay, hypotonia, joint hypermobility and relative macrocephaly who carried small, overlapping deletions of Xp11.22. The maximum region of overlap between their deletions spanned ~430 kb and included two pseudogenes, CENPVL1 and CENPVL2, whose functions are not known, and two protein coding genes—the G1 to S phase transition 2 gene (GSPT2) and the MAGE family member D1 gene (MAGED1). Deletions of this ~430 kb region have not been previously implicated in human disease. Duplications of GSPT2 have been documented in individuals with intellectual disability, but the phenotypic consequences of a loss of GSPT2 function have not been elucidated in humans or mouse models. Changes in MAGED1 have not been associated with intellectual disability in humans, but loss of MAGED1 function is associated with neurocognitive and neurobehavioral phenotypes in mice. In all cases, the Xp11.22 deletion was inherited from an unaffected mother. Studies performed on DNA from one of these mothers did not show evidence of skewed X-inactivation. These results suggest that deletions of an ~430 kb region on chromosome Xp11.22 that encompass CENPVL1, CENPVL2, GSPT2 and MAGED1 cause a distinct X-linked syndrome characterized by intellectual disability, developmental delay, hypotonia, joint hypermobility and relative macrocephaly. Loss of GSPT2 and/or MAGED1 function may contribute to the intellectual disability and developmental delay seen in males with these deletions. [ABSTRACT FROM AUTHOR]

Traditionally, the tribe Sophoreae sensu lato has been considered a basal but also heterogeneous taxonomic group of the papilionoid legumes. Phylogenetic studies have placed Sophoreae sensu stricto (s.s.) as a member of the core genistoids. The recently suggested new circumscription of this tribe involved the removal of traditional members and the inclusion of Euchresteae and Thermopsideae. Nonetheless, definitions and inter- and intra-taxonomic issues of Sophoreae remain unclear. Within the field of legume systematics, the molecular characteristics of a plastid genome (plastome) have an important role in helping to define taxonomic groups. Here, we examined the plastome of Maackia fauriei, belonging to Sophoreae s.s., to elucidate the molecular characteristics of Sophoreae. Its gene contents are similar to the plastomes of other typical legumes. Putative pseudogene rps16 of Maackia and Lupinus species imply independent functional gene loss from the genistoids. Our overall examination of that loss among legumes suggests that it is common among all major clades of Papilionoideae. The M. fauriei plastome has a novel 24-kb inversion in its large single copy region, as well as previously recognized 50-kb and 36-kb inversions. The 36-kb inversion is shared by the core genistoids. The 24-kb inversion is present in the eight genera belonging to three tribes: Euchresteae, Sophoreae s.s., and Thermopsideae. The phylogenetic distribution of this 24-kb inversion strongly supports the monophyly of members of Sophoreae s.s. with Euchresteae and Thermopsideae. Hence, it can be used as a putative synapomorphic characteristic for the newly circumscribed Sophoreae, including Euchresteae and Thermopsideae. However, plastome conformation suggests a slightly narrower taxonomic group because of heterogeneous results from Bolusanthus and Dicraeopetalum. The phylogenetic analysis, based on plastome sequences from 43 legumes, represents well our understanding of legume systematics while resolving the genistoid clade as a sister group to an Old World clade. It also demonstrates the value that plastomes are powerful marker for systematic studies of basal papilionoid legumes. [ABSTRACT FROM AUTHOR]

A portion of the mitochondrial cytochrome c oxidase I gene was sequenced using both genomic DNA and complement DNA from three planktonic copepod Neocalanus species (N. cristatus, N. plumchrus, and N. flemingeri). Small but critical sequence differences in CO1 were observed between gDNA and cDNA from N. plumchrus. Furthermore, careful observation revealed the presence of recombination between sequences in gDNA from N. plumchrus. Moreover, a chimera of the N. cristatus and N. plumchrus sequences was obtained from N. plumchrus gDNA. The observed phenomena can be best explained by the preferential amplification of the nuclear mitochondrial pseudogenes from gDNA of N. plumchrus. Two conclusions can be drawn from the observations. First, nuclear mitochondrial pseudogenes are pervasive in N. plumchrus. Second, a mating between a female N. cristatus and a male N. plumchrus produced viable offspring, which further backcrossed to a N. plumchrus individual. These observations not only demonstrate intriguing mating behavior in these species, but also emphasize the importance of careful interpretation of species marker sequences amplified from gDNA. [ABSTRACT FROM AUTHOR]

Since the discovery of the high abundance of Alu elements in the human genome, the interest for the functional significance of these retrotransposons has been increasing. Primate Alu and rodent Alu-like elements are retrotransposed by a mechanism driven by the LINE1 (L1) encoded proteins, the same machinery that generates the L1 repeats, the processed pseudogenes (PPs), and other retroelements. Apart from free Alu RNAs, Alus are also transcribed and retrotranscribed as part of cellular gene transcripts, generally embedded inside 3’ untranslated regions (UTRs). Despite different proposed hypotheses, the functional implication of the presence of Alus inside 3’UTRs remains elusive. In this study we hypothesized that Alu elements in 3’UTRs could be involved in the genesis of PPs. By analyzing human genome data we discovered that the existence of 3’UTR-embedded Alu elements is overrepresented in genes source of PPs. In contrast, the presence of other retrotransposable elements in 3’UTRs does not show this PP linked overrepresentation. This research was extended to mouse and rat genomes and the results accordingly reveal overrepresentation of 3’UTR-embedded B1 (Alu-like) elements in PP parent genes. Interestingly, we also demonstrated that the overrepresentation of 3’UTR-embedded Alus is particularly significant in PP parent genes with low germline gene expression level. Finally, we provide data that support the hypothesis that the L1 machinery is also the system that herpesviruses, and possibly other large DNA viruses, use to capture host genes expressed in germline or somatic cells. Altogether our results suggest a novel role for Alu or Alu-like elements inside 3’UTRs as facilitators of the genesis of PPs, particularly in lowly expressed genes. Moreover, we propose that this L1-driven mechanism, aided by the presence of 3’UTR-embedded Alus, may also be exploited by DNA viruses to incorporate host genes to their viral genomes. [ABSTRACT FROM AUTHOR]

Phosphatase and tensin homolog (PTEN) is a tumor-suppressor gene. PTEN pseudogene (PTENp) acts as an endogenous RNA, which regulates its parental gene by competitively binding to the 3’ UTR of PTEN gene in the human. Despite the importance of this pseudogene, little is known about the molecular evolution of PTENp in mammals. In this study, we identified 37 pseudogenes from 65 mammalian genomes. Among them, 32 were from rodents or primates. Phylogenetic analyse showed a complex evolutionary history of this gene family. Some PTENps were shared both in primates and rodents. However, some PTENps were shown to be species-specific, such as the tasmanian devil PTENp1, nine banded armadillo PTENp1 and gibbon PTENp1. Most interestingly, the naked mole rat (NMR), an anticancer model organism, possessed 17 copies of PTENps, which were classified into four clades based on the phylogenetic analyses. Furthermore, we found that all the 3’UTR of PTEN and PTENps shared common microRNA (MicroRNA) binding sites in NMR, based on our prediction of specific MicroRNA binding sites. Our findings suggested that multiple gene duplications have occurred in the formation of PTEN/PTENp gene family during the evolution of mammals. Some PTENps were relatively ancient and were shared by primates and rodents; others were newly originated through species- specific gene duplications. PTENps in NMR may function as competitive endogenous RNAs (ceRNAs) to regulate their counterpart genes by competing for common MicroRNAs, which may be one of the interpretations for the cancer resistance in NMR. [ABSTRACT FROM AUTHOR]

We have developed a new bioinformatics framework for the analysis of rearranged bovine heavy chain immunoglobulin (Ig) variable regions by combining and refining widely used alignment algorithms. This bioinformatics framework allowed us to investigate alignments of heavy chain framework regions (FRHs) and the separate alignments of FRHs and heavy chain complementarity determining regions (CDRHs) to determine their germline origin in the four cattle breeds Aubrac, German Black Pied, German Simmental, and Holstein Friesian. Now it is also possible to specifically analyze Ig heavy chains possessing exceptionally long CDR3Hs. In order to gain more insight into breed specific differences in Ig combinatorial diversity, somatic hypermutations and putative gene conversions of IgG, we compared the dominantly transcribed variable (IGHV), diversity (IGHD), and joining (IGHJ) segments and their recombination in the four cattle breeds. The analysis revealed the use of 15 different IGHV segments, 21 IGHD segments, and two IGHJ segments with significant different transcription levels within the breeds. Furthermore, there are preferred rearrangements within the three groups of CDR3H lengths. In the sequences of group 2 (CDR3H lengths (L) of 11–47 amino acid residues (aa)) a higher number of recombination was observed than in sequences of group 1 (L≤10 aa) and 3 (L≥48 aa). The combinatorial diversity of germline IGHV, IGHD, and IGHJ-segments revealed 162 rearrangements that were significantly different. The few preferably rearranged gene segments within group 3 CDR3H regions may indicate specialized antibodies because this length is unique in cattle. The most important finding of this study, which was enabled by using the bioinformatics framework, is the discovery of strong evidence for gene conversion as a rare event using pseudogenes fulfilling all definitions for this particular diversification mechanism. [ABSTRACT FROM AUTHOR]

Joinvilleaceae is a family of tropical grass-like monocots that comprises only the genus Joinvillea. Previous studies have placed Joinvilleaceae in close phylogenetic proximity to the well-studied grass family. A full plastome sequence was determined and characterized for J. ascendens. The plastome was sequenced with next generation methods, fully assembled de novo and annotated. The assembly revealed two novel inversions specific to the Joinvilleaceae lineage and at least one novel plastid inversion in the Joinvilleaceae-Poaceae lineage. Two previously documented inversions in the Joinvilleaceae-Poaceae lineage and one previously documented inversion in the Poaceae lineage were also verified. Inversion events were identified visually and verified computationally by simulation mutations. Additionally, the loss and subsequent degradation of the accD gene in order Poales was explored extensively in Poaceae and J. ascendens. The two novel inversions along with changes in gene composition between families better delimited lineages in the Poales. The presence of large inversions and subsequent reversals in this small family suggested a high potential for large-scale rearrangements to occur in plastid genomes. [ABSTRACT FROM AUTHOR]

Francisella tularensis is composed of a number of subspecies with varied geographic distribution, host ranges, and virulence. In view of these marked differences, comparative functional genomics may elucidate some of the molecular mechanism(s) behind these differences. In this study a shared probe microarray was designed that could be used to compare the transcriptomes of Francisella tularensis subsp. tularensis Schu S4 (Ftt), Francisella tularensis subsp. holarctica OR960246 (Fth), Francisella tularensis subsp. holarctica LVS (LVS), and Francisella novicida U112 (Fn). To gain insight into expression differences that may be related to the differences in virulence of these subspecies, transcriptomes were measured from each strain grown in vitro under identical conditions, utilizing a shared probe microarray. The human avirulent Fn strain exhibited high levels of transcription of genes involved in general metabolism, which are pseudogenes in the human virulent Ftt and Fth strains, consistent with the process of genome decay in the virulent strains. Genes encoding an efflux system (emrA2 cluster of genes), siderophore (fsl operon), acid phosphatase, LPS synthesis, polyamine synthesis, and citrulline ureidase were all highly expressed in Ftt when compared to Fn, suggesting that some of these may contribute to the relative high virulence of Ftt. Genes expressed at a higher level in Ftt when compared to the relatively less virulent Fth included genes encoding isochorismatases, cholylglycine hydrolase, polyamine synthesis, citrulline ureidase, Type IV pilus subunit, and the Francisella Pathogenicity Island protein PdpD. Fth and LVS had very few expression differences, consistent with the derivation of LVS from Fth. This study demonstrated that a shared probe microarray designed to detect transcripts in multiple species/subspecies of Francisella enabled comparative transcriptional analyses that may highlight critical differences that underlie the relative pathogenesis of these strains for humans. This strategy could be extended to other closely-related bacterial species for inter-strain and inter-species analyses. [ABSTRACT FROM AUTHOR]

Scrophularia dentata is an important Tibetan medicinal plant and traditionally used for the treatment of exanthema and fever in Traditional Tibetan Medicine (TTM). However, there is little sequence and genomic information available for S. dentata. In this paper, we report the complete chloroplast genome sequence of S. dentata and it is the first sequenced member of the Sect. Tomiophyllum within Scrophularia (Scrophulariaceae). The gene order and organization of the chloroplast genome of S. dentata are similar to other Lamiales chloroplast genomes. The plastome is 152,553 bp in length and includes a pair of inverted repeats (IRs) of 25,523 bp that separate a large single copy (LSC) region of 84,058 bp and a small single copy (SSC) region of 17,449 bp. It has 38.0% GC content and includes 114 unique genes, of which 80 are protein-coding, 30 are transfer RNA, and 4 are ribosomal RNA. Also, it contains 21 forward repeats, 19 palindrome repeats and 41 simple sequence repeats (SSRs). The repeats and SSRs within S. dentata were compared with those of S. takesimensis and present certain discrepancies. The chloroplast genome of S. dentata was compared with other five publicly available Lamiales genomes from different families. All the coding regions and non-coding regions (introns and intergenic spacers) within the six chloroplast genomes have been extracted and analysed. Furthermore, the genome divergent hotspot regions were identified. Our studies could provide basic data for the alpine medicinal species conservation and molecular phylogenetic researches of Scrophulariaceae and Lamiales. [ABSTRACT FROM AUTHOR]

Bacillus pumilus strain SAFR-032 is a non-pathogenic spore-forming bacterium exhibiting an anomalously high persistence in bactericidal environments. In its dormant state, it is capable of withstanding doses of ultraviolet (UV) radiation or hydrogen peroxide, which are lethal for the vast majority of microorganisms. This unusual resistance profile has made SAFR-032 a reference strain for studies of bacterial spore resistance. The complete genome sequence of B. pumilus SAFR-032 was published in 2007 early in the genomics era. Since then, the SAFR-032 strain has frequently been used as a source of genetic/genomic information that was regarded as representative of the entire B. pumilus species group. Recently, our ongoing studies of conservation of gene distribution patterns in the complete genomes of various B. pumilus strains revealed indications of misassembly in the B. pumilus SAFR-032 genome. Synteny-driven local genome resequencing confirmed that the original SAFR-032 sequence contained assembly errors associated with long sequence repeats. The genome sequence was corrected according to the new findings. In addition, a significantly improved annotation is now available. Gene orders were compared and portions of the genome arrangement were found to be similar in a wide spectrum of Bacillus strains. [ABSTRACT FROM AUTHOR]

We have previously described a mouse model of human posterior polymorphous corneal dystrophy (PPCD) and localized the causative mutation to a 6.2 Mbp region of chromosome 2, termed Ppcd1. We now show that the gene rearrangement linked to mouse Ppcd1 is a 3.9 Mbp chromosomal inversion flanked by 81 Kbp and 542 bp deletions. This recombination event leads to deletion of Csrp2bp Exons 8 through 11, Dzank1 Exons 20 and 21, and the pseudogene Znf133. In addition, we identified translocation of novel downstream sequences to positions adjacent to Csrp2bp Exon 7 and Dzank1 Exon 20. Twelve novel fusion transcripts involving Csrp2bp or Dzank1 linked to downstream sequences have been identified. Eight are expressed at detectable levels in PPCD1 but not wildtype eyes. Upregulation of two Csrp2bp fusion transcripts, as well as upregulation of the adjacent gene, Ovol2, was observed. Absence of the PPCD1 phenotype in animals haploinsufficient for Csrp2bp or both Csrp2bp and Dzank1 rules out haploinsufficiency of these genes as a cause of mouse PPCD1. Complementation experiments confirm that PPCD1 embryonic lethality is due to disruption of Csrp2bp expression. The ocular expression pattern of Csrp2bp is consistent with a role for this protein in corneal development and pathogenesis of PPCD1. [ABSTRACT FROM AUTHOR]

Background: Prothymosin α (ProTα) (isoform 2: iso2) is a widely distributed, small acidic protein with intracellular and extracellular-associated functions. Recently, we identified two new ProTα variants with potent anti-HIV activity from CD8+ T cells and cervicovaginal lavage. The first is a splice variant of the ProTα gene known as isoB and the second is the product of ProTα pseudogene 7 (p7). Similarly to iso2, the anti-HIV activity of both variants is mediated by type I IFN. Here we tested whether the immunomodulatory activity of isoB and p7 are also TLR4 dependent and determined their kinetic of release in response to HIV-1 infection. Methods: Type I, type III, TNF-α and IL-6 mRNA inducing activity was determined in macrophages from wild type and TLR4 knockout mice treated with recombinant ProTα variants. Supernatants from mock and HIV infected cells were analyzed by mass spectrometry in positive and negative modes for the presence of ProTα variants. In silico structural and functional analysis of ProTα variants were performed. Results: We show that both isoB and p7 upregulate IFN-β, IFN-λ1, IL-6, TNF-α and RANTES mRNAs in primary human macrophages. The potent stimulation of IFN-β by the recombinant ProTα variants in human macrophages is dependent on the TLR4 pathway, whereas the induction of TNF-α and IL-6 may also occur independently of TLR4, suggesting the interaction of ProTα variants with other signaling molecules/receptors. In silico analyses confirmed that the novel isoB and p7 variants are intrinsically disordered proteins, which lack the NLS and mass spectrometry showed release of ProTα variants within minutes post HIV-1 infection. These features are consistent with the function of ProTα variants as damage associate molecular patterns (DAMPs). Conclusions: Our findings indicate that ProTα variants strongly inhibit viral replication mainly, but not exclusively, through TLR4 signaling and that they are released within minutes of viral infection suggesting that they may function as DAMPs. [ABSTRACT FROM AUTHOR]

The phist gene family has members identified across the Plasmodium genus, defined by the presence of a domain of roughly 150 amino acids having conserved aromatic residues and an all alpha-helical structure. The family is highly amplified in P. falciparum, with 65 predicted genes in the genome of the 3D7 isolate. In contrast, in the rodent malaria parasite P. berghei 3 genes are identified, one of which is an apparent pseudogene. Transcripts of the P. berghei phist genes are predominant in schizonts, whereas in P. falciparum transcript profiles span different asexual blood stages and gametocytes. We pursued targeted disruption of P. berghei phist genes in order to characterize a simplistic model for the expanded phist gene repertoire in P. falciparum. Unsuccessful attempts to disrupt P. berghei PBANKA_114540 suggest that this phist gene is essential, while knockout of phist PBANKA_122900 shows an apparent normal progression and non-essential function throughout the life cycle. Epitope-tagging of P. falciparum and P. berghei phist genes confirmed protein export to the erythrocyte cytoplasm and localization with a punctate pattern. Three P. berghei PEXEL/HT-positive exported proteins exhibit at least partial co-localization, in support of a common vesicular compartment in the cytoplasm of erythrocytes infected with rodent malaria parasites. [ABSTRACT FROM AUTHOR]

Ferritin, the major iron storage protein, performs its essential functions in the cytoplasm, nucleus and mitochondria. The variable assembly of 24 subunits of the Heavy (H) and Light (L) type composes the cytoplasmic molecule. In humans, two distinct genes code these subunits, both belonging to complex multigene families. Until now, one H gene has been identified with the coding sequence interrupted by three introns and more than 20 intronless copies widely dispersed on different chromosomes. Two of the intronless genes are actively transcribed in a tissue-specific manner. Herein, we report that FTH1P3, another intronless pseudogene, is transcribed. FTH1P3 transcript was detected in several cell lines and tissues, suggesting that its transcription is ubiquitary, as it happens for the parental ferritin H gene. Moreover, FTH1P3 expression is positively regulated during the cell differentiation process. [ABSTRACT FROM AUTHOR]

Trace amine-associated receptors (TAARs) are a member of the G-protein-coupled receptor superfamily and are known to be expressed in olfactory sensory neurons. A limited number of molecular evolutionary studies have been done for TAARs so far. To elucidate how lineage-specific evolution contributed to their functional divergence, we examined 30 metazoan genomes. In total, 493 TAAR gene candidates (including 84 pseudogenes) were identified from 26 vertebrate genomes. TAARs were not identified from non-vertebrate genomes. An ancestral-type TAAR-like gene appeared to have emerged in lamprey. We found four therian-specific TAAR subfamilies (one eutherian-specific and three metatherian-specific) in addition to previously known nine subfamilies. Many species-specific TAAR gene duplications and losses contributed to a large variation of TAAR gene numbers among mammals, ranging from 0 in dolphin to 26 in flying fox. TAARs are classified into two groups based on binding preferences for primary or tertiary amines as well as their sequence similarities. Primary amine-detecting TAARs (TAAR1-4) have emerged earlier, generally have single-copy orthologs (very few duplication or loss), and have evolved under strong functional constraints. In contrast, tertiary amine-detecting TAARs (TAAR5-9) have emerged more recently and the majority of them experienced higher rates of gene duplications. Protein members that belong to the tertiary amine-detecting TAAR group also showed the patterns of positive selection especially in the area surrounding the ligand-binding pocket, which could have affected ligand-binding activities and specificities. Expansions of the tertiary amine-detecting TAAR gene family may have played important roles in terrestrial adaptations of therian mammals. Molecular evolution of the TAAR gene family appears to be governed by a complex, species-specific, interplay between environmental and evolutionary factors. [ABSTRACT FROM AUTHOR]

Plants from the family Orobanchaceae are widely used as a model to study different aspects of parasitic lifestyle including host–parasite interactions and physiological and genomic adaptations. Among the latter, the most prominent are those that occurred due to the loss of photosynthesis; they include the reduction of the photosynthesis-related gene set in both nuclear and plastid genomes. In Orobanchaceae, the transition to non-photosynthetic lifestyle occurred several times independently, but only one lineage has been in the focus of evolutionary studies. These studies included analysis of plastid genomes and transcriptomes and allowed the inference of patterns and mechanisms of genome reduction that are thought to be general for parasitic plants. Here we report the plastid genome of Lathraea squamaria, a holoparasitic plant from Orobanchaceae, clade Rhinantheae. We found that in this plant the degree of plastome reduction is the least among non-photosynthetic plants. Like other parasites, Lathraea possess a plastome with elevated absolute rate of nucleotide substitution. The only gene lost is petL, all other genes typical for the plastid genome are present, but some of them–those encoding photosystem components (22 genes), cytochrome b6/f complex proteins (4 genes), plastid-encoded RNA polymerase subunits (2 genes), ribosomal proteins (2 genes), ccsA and cemA–are pseudogenized. Genes for cytochrome b6/f complex and photosystems I and II that do not carry nonsense or frameshift mutations have an increased ratio of non-synonymous to synonymous substitution rates, indicating the relaxation of purifying selection. Our divergence time estimates showed that transition to holoparasitism in Lathraea lineage occurred relatively recently, whereas the holoparasitic lineage Orobancheae is about two times older. [ABSTRACT FROM AUTHOR]

Reference genes used in normalizing qRT-PCR data are critical for the accuracy of gene expression analysis. However, many traditional reference genes used in zebrafish early development are not appropriate because of their variable expression levels during embryogenesis. In the present study, we used our previous RNA-Seq dataset to identify novel reference genes suitable for gene expression analysis during zebrafish early developmental stages. We first selected 197 most stably expressed genes from an RNA-Seq dataset (29,291 genes in total), according to the ratio of their maximum to minimum RPKM values. Among the 197 genes, 4 genes with moderate expression levels and the least variation throughout 9 developmental stages were identified as candidate reference genes. Using four independent statistical algorithms (delta-CT, geNorm, BestKeeper and NormFinder), the stability of qRT-PCR expression of these candidates was then evaluated and compared to that of actb1 and actb2, two commonly used zebrafish reference genes. Stability rankings showed that two genes, namely mobk13 (mob4) and lsm12b, were more stable than actb1 and actb2 in most cases. To further test the suitability of mobk13 and lsm12b as novel reference genes, they were used to normalize three well-studied target genes. The results showed that mobk13 and lsm12b were more suitable than actb1 and actb2 with respect to zebrafish early development. We recommend mobk13 and lsm12b as new optimal reference genes for zebrafish qRT-PCR analysis during embryogenesis and early larval stages. [ABSTRACT FROM AUTHOR]

The human RB1 gene is imprinted due to integration of the PPP1R26P1 pseudogene into intron 2. PPP1R26P1 harbors the gametic differentially methylated region of the RB1 gene, CpG85, which is methylated in the female germ line. The paternally unmethylated CpG85 acts as promoter for the alternative transcript 2B of RB1, which interferes with expression of full-length RB1 in cis. In mice, PPP1R26P1 is not present in the Rb1 gene and Rb1 is not imprinted. Assuming that the mechanisms responsible for genomic imprinting are conserved, we investigated if imprinting of mouse Rb1 can be induced by transferring human PPP1R26P1 into mouse Rb1. We generated humanized Rb1_PPP1R26P1 knock-in mice that pass human PPP1R26P1 through the mouse germ line. We found that the function of unmethylated CpG85 as promoter for an alternative Rb1 transcript and as cis-repressor of the main Rb1 transcript is maintained in mouse tissues. However, CpG85 is not recognized as a gametic differentially methylated region in the mouse germ line. DNA methylation at CpG85 is acquired only in tissues of neuroectodermal origin, independent of parental transmission of PPP1R26P1. Absence of CpG85 methylation in oocytes and sperm implies a failure of imprint methylation establishment in the germ line. Our results indicate that site-specific integration of a proven human gametic differentially methylated region is not sufficient for acquisition of DNA methylation in the mouse germ line, even if promoter function of the element is maintained. This suggests a considerable dependency of DNA methylation induction on the surrounding sequence. However, our model is suited to determine the cellular function of the alternative Rb1 transcript. [ABSTRACT FROM AUTHOR]

Background: Vibrio parahaemolyticus is a recognized human enteropathogen. Thermostable direct hemolysin (TDH) and TDH-related hemolysin (TRH) as well as the type III secretion system 2 (T3SS2) are considered as major virulence factors. As tdh positive strains are not detected in coastal waters of Germany, we focused on the characterization of trh positive strains, which were isolated from mussels, seawater and patients in Germany. Results: Ten trh harbouring V. parahaemolyticus strains from Germany were compared to twenty-one trh positive strains from other countries. The complete trh sequences revealed clustering into three different types: trh1 and trh2 genes and a pseudogene Ψtrh. All German isolates possessed alleles of the trh2 gene. MLST analysis indicated a close relationship to Norwegian isolates suggesting that these strains belong to the autochthonous microflora of Northern Europe seawaters. Strains carrying the pseudogene Ψtrh were negative for T3SS2β effector vopC. Transcription of trh and vopC genes was analyzed under different growth conditions. Trh2 gene expression was not altered by bile while trh1 genes were inducible. VopC could be induced by urea in trh2 bearing strains. Most trh1 carrying strains were hemolytic against sheep erythrocytes while all trh2 positive strains did not show any hemolytic activity. TRH variants were synthesized in a prokaryotic cell-free system and their hemolytic activity was analyzed. TRH1 was active against sheep erythrocytes while TRH2 variants were not active at all. Conclusion: Our study reveals a high diversity among trh positive V. parahaemolyticus strains. The function of TRH2 hemolysins and the role of the pseudogene Ψtrh as pathogenicity factors are questionable. To assess the pathogenic potential of V. parahaemolyticus strains a differentiation of trh variants and the detection of T3SS2β components like vopC would improve the V. parahaemolyticus diagnostics and could lead to a refinement of the risk assessment in food analyses and clinical diagnostics. [ABSTRACT FROM AUTHOR]

Serum amyloid A3 (SAA3) possesses characteristics distinct from the other serum amyloid A isoforms, SAA1, SAA2, and SAA4. High density lipoprotein contains the latter three isoforms, but not SAA3. The expression of mouse SAA3 (mSAA3) is known to be up-regulated extrahepatically in inflammatory responses, and acts as an endogenous ligand for the toll-like receptor 4/MD-2 complex. We previously reported that mSAA3 plays an important role in facilitating tumor metastasis by attracting circulating tumor cells and enhancing hyperpermeability in the lungs. On the other hand, human SAA3 (hSAA3) has long been regarded as a pseudogene, which is in contrast to the abundant expression levels of the other isoforms. Although the nucleotide sequence of hSAA3 is very similar to that of the other SAAs, a single oligonucleotide insertion in exon 2 causes a frame-shift to generate a unique amino acid sequence. In the present study, we identified that hSAA3 was transcribed in the hSAA2-SAA3 fusion transcripts of several human cell lines. In the fusion transcript, hSAA2 exon 3 was connected to hSAA3 exon 1 or hSAA3 exon 2, located approximately 130kb downstream from hSAA2 exon 3 in the genome, which suggested that it is produced by alternative splicing. Furthermore, we succeeded in detecting and isolating hSAA3 protein for the first time by an immunoprecipitation-enzyme linked immune assay system using monoclonal and polyclonal antibodies that recognize the hSAA3 unique amino acid sequence. We also demonstrated that hSAA3 bound oxidized low density lipoprotein receptor (oxLDL receptor, LOX-1) and elevated the phosphorylation of ERK, the intracellular MAP-kinase signaling protein. [ABSTRACT FROM AUTHOR]