Search

The honey bee trypanosomatid parasite, Lotmaria passim, contains two genes that encode the flagellar calcium binding protein (FCaBP) through tandem duplication in its genome. FCaBPs localize in the flagellum and entire body membrane of L. passim through specific N-terminal sorting sequences. This finding suggests that this is an example of protein subcellular relocalization resulting from gene duplication, altering the intracellular localization of FCaBP. However, this phenomenon may not have occurred in Leishmania, as one or both of the duplicated genes have become pseudogenes. Multiple copies of the FCaBP gene are present in several Trypanosoma species and Leptomonas pyrrhocoris, indicating rapid evolution of this gene in trypanosomatid parasites. The N-terminal flagellar sorting sequence of L. passim FCaBP1 is in close proximity to the BBSome complex, while that of Trypanosoma brucei FCaBP does not direct GFP to the flagellum in L. passim. Deletion of the two FCaBP genes in L. passim affected growth and impaired flagellar morphogenesis and motility, but it did not impact host infection. Therefore, FCaBP represents a duplicated gene with a rapid evolutionary history that is essential for flagellar structure and function in a trypanosomatid parasite. Author summary: Protein subcellular relocalization (PSR) was proposed as a mechanism for the functional divergence and retention of duplicate genes. We explored this hypothesis using flagellar calcium binding protein (FCaBP) genes duplicated in many trypanosomatid parasites. FCaBPs localize in the flagellum and entire body membrane of honey bee trypanosomatid parasite, Lotmaria passim through specific N-terminal sorting sequences. However, this is not common with all trypanosomatid parasites since one or both of the duplicated genes have become pseudogenes in Leishmania. N-terminal flagellar sorting sequence of FCaBP1 interacts with BBSome complex, and thus they must have co-evolved in each trypanosomatid species. FCaBPs are essential for the normal growth, flagellar morphogenesis, and motility but not host infection of L. passim. Our findings demonstrate that duplicated FCaBP genes have undergone PSR by changing amino acids in the N-terminal end with some but not all trypanosomatid parasites. In the parasites with multiple FCaBP genes, they play critical roles for the growth control as well as flagellar structure and function. [ABSTRACT FROM AUTHOR]

Small nuclear RNAs (snRNAs) are structural and functional cores of the spliceosome. In metazoan genomes, each snRNA has multiple copies/variants, up to hundreds in mammals. However, the expressions and functions of each copy/variant in one organism have not been systematically studied. Focus on U1 snRNA genes, we investigated all five copies in Drosophila melanogaster using two series of constructed strains. Analyses of transgenic flies that each have a U1 promoter-driven gfp revealed that U1:21D is the major and ubiquitously expressed copy, and the other four copies have specificities in developmental stages and tissues. Mutant strains that each have a precisely deleted copy of U1-gene exhibited various extents of defects in fly morphology or mobility, especially deletion of U1:82Eb. Interestingly, splicing was changed at limited levels in the deletion strains, while large amounts of differentially-expressed genes and alternative polyadenylation events were identified, showing preferences in the down-regulation of genes with 1–2 introns and selection of proximal sites for 3'-end polyadenylation. In vitro assays suggested that Drosophila U1 variants pulled down fewer SmD2 proteins compared to the canonical U1. This study demonstrates that all five U1-genes in Drosophila have physiological functions in development and play regulatory roles in transcription and 3'-end formation. Author summary: Small nuclear RNAs (snRNAs) are the core of the spliceosome, which catalyzes the removal of introns. In metazoan genomes, each snRNA has multiple copies with variant sequences, up to hundreds in mammals, in which many singular copies/variants have shown function and are related to human diseases. However, expressions and functions of all snRNA copies in one organism have remained unknown. U1 snRNA is critical for early intron recognition and spliceosome assembly. The Drosophila genome has five U1 genes that generate three U1 variants. Here, we report that the U1:21D is ubiquitously high in all tissues and developmental stages, 95Ca is high in the salivary gland and male accessory gland, and 82Eb, 95Cb, and 95Cc are expressed in the ovary and embryo. Deletions of each singular U1-gene cause various developmental defects, demonstrating that none of them are pseudogenes. The absence of a U1-gene results in more serious defects in transcription and 3'-end formation than in pre-mRNA splicing, showing large amounts of expression-changed genes and alternative polyadenylation changes but much fewer changes in alternative splicing. Thus, our findings advance the understanding of in vivo functions of all the U1 snRNA genes in one organism. [ABSTRACT FROM AUTHOR]

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]

CYP2D6 is a very important pharmacogene as it is responsible for the metabolization or bioactivation of 20 to 30% of the clinically used drugs. However, despite its relatively small length of only 4.4 kb, it is one of the most challenging pharmacogenes to genotype due to the high similarity with its neighboring pseudogenes and the frequent occurrence of CYP2D6-CYP2D7 hybrids. Unfortunately, most current genotyping methods are therefore not able to correctly determine the complete CYP2D6-CYP2D7 sequence. Therefore, we developed a genotyping assay to generate complete allele-specific consensus sequences of complex regions by optimizing the PCR-free nanopore Cas9-targeted sequencing (nCATS) method combined with adaptive sequencing, and developing a new comprehensive long read genotyping (CoLoRGen) pipeline. The CoLoRGen pipeline first generates consensus sequences of both alleles and subsequently determines both large structural and small variants to ultimately assign the correct star-alleles. In reference samples, our genotyping assay confirms the presence of CYP2D6-CYP2D7 large structural variants, single nucleotide variants (SNVs), and small insertions and deletions (INDELs) that go undetected by most current assays. Moreover, our results provide direct evidence that the CYP2D6 genotype of the NA12878 DNA should be updated to include the CYP2D6-CYP2D7 *68 hybrid and several additional single nucleotide variants compared to existing references. Ultimately, the nCATS-CoLoRGen genotyping assay additionally allows for more accurate gene function predictions by enabling the possibility to detect and phase de novo mutations in addition to known large structural and small variants. Author summary: During the last decades, the usefulness of personalized medicine has become increasingly apparent. Directly linked to that is the need for accurate genotyping assays to determine the pharmacogenetic profile of patients. Continuing research has led to the development of genotyping assays that perform quite robustly. However, complex genes remain an issue when it comes to determining the complete sequence correctly. An example of such a complex but very important pharmacogene is CYP2D6. Therefore, we developed a genotyping assay in an attempt to generate complete allele-specific consensus sequences of CYP2D6, by optimizing a targeted amplification-free long-read sequencing method and developing a new analysis pipeline. In reference samples, we showed that our genotyping assay performed accurately and confirmed the presence of variants that go undetected by most current assays. However, the implementation of this assay in practice is still hampered as the selected enrichment strategies inherently lead to a low percentage of on-target reads, resulting in low on-target sequencing depths. Further optimization and validation of the assay is thus needed, but definitely worth considering for follow-up research as we already demonstrated the added value for generating more complete genotypes, which on its turn will result in more accurate gene function predictions. [ABSTRACT FROM AUTHOR]

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]

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]

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]

Background: Pseudogenes 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. Results: We 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. Conclusions: We 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. [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]

Emerging pathogens are a major threat to public health, however understanding how pathogens adapt to new niches remains a challenge. New methods are urgently required to provide functional insights into pathogens from the massive genomic data sets now being generated from routine pathogen surveillance for epidemiological purposes. Here, we measure the burden of atypical mutations in protein coding genes across independently evolved Salmonella enterica lineages, and use these as input to train a random forest classifier to identify strains associated with extraintestinal disease. Members of the species fall along a continuum, from pathovars which cause gastrointestinal infection and low mortality, associated with a broad host-range, to those that cause invasive infection and high mortality, associated with a narrowed host range. Our random forest classifier learned to perfectly discriminate long-established gastrointestinal and invasive serovars of Salmonella. Additionally, it was able to discriminate recently emerged Salmonella Enteritidis and Typhimurium lineages associated with invasive disease in immunocompromised populations in sub-Saharan Africa, and within-host adaptation to invasive infection. We dissect the architecture of the model to identify the genes that were most informative of phenotype, revealing a common theme of degradation of metabolic pathways in extraintestinal lineages. This approach accurately identifies patterns of gene degradation and diversifying selection specific to invasive serovars that have been captured by more labour-intensive investigations, but can be readily scaled to larger analyses. [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]

The repeated rDNA array gives rise to the nucleolus, an organelle that is central to cellular processes as varied as stress response, cell cycle regulation, RNA modification, cell metabolism, and genome stability. The rDNA array is also responsible for the production of more than 70% of all cellular RNAs (the ribosomal RNAs). The rRNAs are produced from two sets of loci: the 5S rDNA array resides exclusively on human chromosome 1 while the 45S rDNA arrays reside on the short arm of five human acrocentric chromosomes. These critical genome elements have remained unassembled and have been excluded from all Hi-C analyses to date. Here we built the first high resolution map of 5S and 45S rDNA array contacts with the rest of the genome combining over 15 billion Hi-C reads from several experiments. The data enabled sufficiently high coverage to map rDNA-genome interactions with 1MB resolution and identify rDNA-gene contacts. The map showed that the 5S and 45S arrays display preferential contact at common sites along the genome but are not themselves sufficiently close to yield 5S-45S Hi-C contacts. Ribosomal DNA contacts are enriched in segments of closed, repressed, and late replicating chromatin, as well as CTCF binding sites. Finally, we identified functional categories whose dispersed genes coalesced in proximity to the rDNA arrays or instead avoided proximity with the rDNA arrays. The observations further our understanding of the spatial localization of rDNA arrays and their contribution to the architecture of the cell nucleus. [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]

Mycobacteriaceae comprises pathogenic species such as Mycobacterium tuberculosis, M. leprae and M. abscessus, as well as non-pathogenic species, for example, M. smegmatis and M. thermoresistibile. Genome comparison and annotation studies provide insights into genome evolutionary relatedness, identify unique and pathogenicity-related genes in each species, and explore new targets that could be used for developing new diagnostics and therapeutics. Here, we present a comparative analysis of ten-mycobacterial genomes with the objective of identifying similarities and differences between pathogenic and non-pathogenic species. We identified 1080 core orthologous clusters that were enriched in proteins involved in amino acid and purine/pyrimidine biosynthetic pathways, DNA-related processes (replication, transcription, recombination and repair), RNA-methylation and modification, and cell-wall polysaccharide biosynthetic pathways. For their pathogenicity and survival in the host cell, pathogenic species have gained specific sets of genes involved in repair and protection of their genomic DNA. M. leprae is of special interest owing to its having the smallest genome (1600 genes and ~1300 psuedogenes), yet poor genome annotation. More than 75% of the pseudogenes were found to have a functional ortholog in the other mycobacterial genomes and belong to protein families such as transferases, oxidoreductases and hydrolases. [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]

Retroduplications come from reverse transcription of mRNAs and their insertion back into the genome. Here, we performed comprehensive discovery and analysis of retroduplications in a large cohort of 2,535 individuals from 26 human populations, as part of 1000 Genomes Phase 3. We developed an integrated approach to discover novel retroduplications combining high-coverage exome and low-coverage whole-genome sequencing data, utilizing information from both exon-exon junctions and discordant paired-end reads. We found 503 parent genes having novel retroduplications absent from the reference genome. Based solely on retroduplication variation, we built phylogenetic trees of human populations; these represent superpopulation structure well and indicate that variable retroduplications are effective population markers. We further identified 43 retroduplication parent genes differentiating superpopulations. This group contains several interesting insertion events, including a SLMO2 retroduplication and insertion into CAV3, which has a potential disease association. We also found retroduplications to be associated with a variety of genomic features: (1) Insertion sites were correlated with regular nucleosome positioning. (2) They, predictably, tend to avoid conserved functional regions, such as exons, but, somewhat surprisingly, also avoid introns. (3) Retroduplications tend to be co-inserted with young L1 elements, indicating recent retrotranspositional activity, and (4) they have a weak tendency to originate from highly expressed parent genes. Our investigation provides insight into the functional impact and association with genomic elements of retroduplications. We anticipate our approach and analytical methodology to have application in a more clinical context, where exome sequencing data is abundant and the discovery of retroduplications can potentially improve the accuracy of SNP calling. [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]

Accurate mapping of next-generation sequencing (NGS) reads to reference genomes is crucial for almost all NGS applications and downstream analyses. Various repetitive elements in human and other higher eukaryotic genomes contribute in large part to ambiguously (non-uniquely) mapped reads. Most available NGS aligners attempt to address this by either removing all non-uniquely mapping reads, or reporting one random or "best" hit based on simple heuristics. Accurate estimation of the mapping quality of NGS reads is therefore critical albeit completely lacking at present. Here we developed a generalized software toolkit "AlignerBoost", which utilizes a Bayesian-based framework to accurately estimate mapping quality of ambiguously mapped NGS reads. We tested AlignerBoost with both simulated and real DNA-seq and RNA-seq datasets at various thresholds. In most cases, but especially for reads falling within repetitive regions, AlignerBoost dramatically increases the mapping precision of modern NGS aligners without significantly compromising the sensitivity even without mapping quality filters. When using higher mapping quality cutoffs, AlignerBoost achieves a much lower false mapping rate while exhibiting comparable or higher sensitivity compared to the aligner default modes, therefore significantly boosting the detection power of NGS aligners even using extreme thresholds. AlignerBoost is also SNP-aware, and higher quality alignments can be achieved if provided with known SNPs. AlignerBoost’s algorithm is computationally efficient, and can process one million alignments within 30 seconds on a typical desktop computer. AlignerBoost is implemented as a uniform Java application and is freely available at . [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]

The study of genomic regions that contain gene copies and structural variation is a major challenge in modern genomics. Unlike variation involving single nucleotide changes, data on the variation of copy number is difficult to collect and few tools exist for analyzing the variation between individuals. The immunoglobulin heavy variable (IGHV) locus, which plays an integral role in the adaptive immune response, is an example of a complex genomic region that varies in gene copy number. Lack of standard methods to genotype this region prevents it from being included in association studies and is holding back the growing field of antibody repertoire analysis. Here we develop a method that takes short reads from high-throughput sequencing and outputs a genetic profile of the IGHV locus with the read coverage depth and a putative nucleotide sequence for each operationally defined gene cluster. Our operationally defined gene clusters aim to address a major challenge in studying the IGHV locus: the high sequence similarity between gene segments in different genomic locations. Tests on simulated data demonstrate that our approach can accurately determine the presence or absence of a gene cluster from reads as short as 70 bp. More detailed resolution on the copy number of gene clusters can be obtained from read coverage depth using longer reads (e.g., ≥ 100 bp). Detail at the nucleotide resolution of single copy genes (genes present in one copy per haplotype) can be determined with 250 bp reads. For IGHV genes with more than one copy, accurate nucleotide-resolution reconstruction is currently beyond the means of our approach. When applied to a family of European ancestry, our pipeline outputs genotypes that are consistent with the family pedigree, confirms existing multigene variants and suggests new copy number variants. This study paves the way for analyzing population-level patterns of variation in IGHV gene clusters in larger diverse datasets and for quantitatively handling regions of copy number variation in other structurally varying and complex loci. [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]

Most humans harbor both CD177neg and CD177pos neutrophils but 1–10% of people are CD177null, placing them at risk for formation of anti-neutrophil antibodies that can cause transfusion-related acute lung injury and neonatal alloimmune neutropenia. By deep sequencing the CD177 locus, we catalogued CD177 single nucleotide variants and identified a novel stop codon in CD177null individuals arising from a single base substitution in exon 7. This is not a mutation in CD177 itself, rather the CD177null phenotype arises when exon 7 of CD177 is supplied entirely by the CD177 pseudogene (CD177P1), which appears to have resulted from allelic gene conversion. In CD177 expressing individuals the CD177 locus contains both CD177P1 and CD177 sequences. The proportion of CD177hi neutrophils in the blood is a heritable trait. Abundance of CD177hi neutrophils correlates with homozygosity for CD177 reference allele, while heterozygosity for ectopic CD177P1 gene conversion correlates with increased CD177neg neutrophils, in which both CD177P1 partially incorporated allele and paired intact CD177 allele are transcribed. Human neutrophil heterogeneity for CD177 expression arises by ectopic allelic conversion. Resolution of the genetic basis of CD177null phenotype identifies a method for screening for individuals at risk of CD177 isoimmunisation. [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]