Search

Staphylococcus aureus is a pathogenic microorganism of humans and animals, able to cause foodborne intoxication due to the production of staphylococcal enterotoxins (SEs) and to resist antibiotic treatment as in the case of methicillin-resistant S. aureus (MRSA). In this study, we performed a genomic characterisation of 12 genetically diverse S. aureus strains isolated from ready-to-eat foods in Algiers (Algeria). Moreover, their ability to produce some classical and new staphylococcal enterotoxins (SEs) was investigated. The 12 S. aureus strains resulted to belong to nine known sequence types (STs) and to the novel ST7199 and ST7200. Furthermore, S. aureus SA46 was assigned to the European clone MRSA-ST80-SCCmec-IV. The 12 strains showed a wide endowment of se and sel (staphylococcal enterotoxin-like toxin) genes (sea, seb, sed, seg, seh, sei, selj, sek, sem, sen, seo, seq, ser, selu2, selw, selx, sey, sel30; ψent1-ψent2), including variants and pseudogenes, and harboured the enterotoxin gene cluster (egc) types 1 and 5. Additionally, they produced various amounts of SEA (64.54–345.02 ng/mL), SEB (2871.28–14739.17 ng/mL), SED (322.70–398.94 ng/mL), SEH (not detectable–239.48 ng/mL), and SER (36,720.10–63,176.06 ng/mL) depending on their genotypes. The genetic determinants related to their phenotypic resistance to β-lactams (blaZ, mecA), ofloxacin (gyrA-S84L), erythromycin (ermB), lincomycin (lmrS), kanamycin (aph(3′)-III, ant(6)-I), and tetracyclin (tet(L), tet(38)) were also detected. A plethora of virulence-related genes, including major virulence genes such as the tst gene, determinant for the toxic shock syndrome toxin-1, and the lukF-PV and lukS-PV genes, encoding the panton-valentine leukocidin (PVL), were present in the S. aureus strains, highlighting their pathogenic potential. Furthermore, a phylogenomic reconstruction including worldwide foodborne S. aureus showed a clear clustering based on ST and geographical origin rather than the source of isolation. [ABSTRACT FROM AUTHOR]

Burkholderia glumae is an important rice pathogen, thus the genomic and evolutionary history may be helpful to control this notorious pathogen. Here, we present two complete genomes of the B. glumae strains HN1 and HN2, which were isolated from diseased rice seed in China. Average nucleotide identity (ANI) analysis shows greater than 99% similarity of the strains HN1 and HN2 with other published B. glumae genomes. Genomic annotation revealed that the genome of strain HN1 consists of five replicons (6,680,415 bp) with an overall G + C content of 68.06%, whereas the genome of strain HN2 comprises of three replicons (6,560,085 bp) with an overall G + C content of 68.34%. The genome of HN1 contains 5434 protein-coding genes, 351 pseudogenes, and 1 CRISPR, whereas the genome of HN2 encodes 5278 protein-coding genes, 357 pseudogenes, and 2 CRISPR. Both strains encode many pathogenic-associated genes (143 genes in HN1 vs. 141 genes in HN2). Moreover, comparative genomic analysis shows the extreme plasticity of B. glumae, which may contribute to its pathogenicity. In total, 259 single-copy genes were affected by positive selection. These genes may contribute to the adaption to different environments. Notably, six genes were characterized as virulence factors which may be an additional way to assist the pathogenicity of B. glumae. [ABSTRACT FROM AUTHOR]

Background: Processed pseudogenes (PΨgs) are disabled gene copies that are transcribed and may affect expression of paralogous genes. Moreover, their insertion in the genome can disrupt the structure or the regulatory region of a gene, affecting its expression level. These events have been identified as occurring mutations during cancer development, thus being able to identify PΨgs and their location will improve their impact on diagnostic testing, not only in cancer but also in inherited disorders. Results: We have implemented PΨFinder (P-psy-finder), a tool that identifies PΨgs, annotates known ones and predicts their insertion site(s) in the genome. The tool screens alignment files and provides user-friendly summary reports and visualizations. To demonstrate its applicability, we scanned 218 DNA samples from patients screened for hereditary colorectal cancer. We detected 423 PΨgs distributed in 96% of the samples, comprising 7 different parent genes. Among these, we confirmed the well-known insertion site of the SMAD4-PΨg within the last intron of the SCAI gene in one sample. While for the ubiquitous CBX3-PΨg, present in 82.6% of the samples, we found it reversed inserted in the second intron of the C15ORF57 gene. Conclusions: PΨFinder is a tool that can automatically identify novel PΨgs from DNA sequencing data and determine their location in the genome with high sensitivity (95.92%). It generates high quality figures and tables that facilitate the interpretation of the results and can guide the experimental validation. PΨFinder is a complementary analysis to any mutational screening in the identification of disease-causing mutations within cancer and other diseases. [ABSTRACT FROM AUTHOR]

Colorectal cancer (CRC) is one of the most common and malignant carcinomas. Many long noncoding RNAs (lncRNAs) have been reported to play important roles in the tumorigenesis of CRC by influencing the expression of some mRNAs via competing endogenous RNA (ceRNA) networks and interacting with miRNAs. Pseudogene is one kind of lncRNA and can act as RNA sponges for miRNAs and regulate gene expression via ceRNA networks. However, there are few studies about pseudogenes in CRC. In this study, 31 differentially expressed (DE) pseudogenes, 17 DE miRNAs and 152 DE mRNAs were identified by analyzing the expression profiles of colon adenocarcinoma obtained from The Cancer Genome Atlas. A ceRNA network was constructed based on these RNAs. Kaplan–Meier analysis showed that 7 pseudogenes, 4 miRNAs and 30 mRNAs were significantly associated with overall survival. Then multivariate Cox regression analysis of the ceRNA-related DE pseudogenes was performed and a 5-pseudogene signature with the greatest prognostic value for CRC was identified. Moreover, the results were validated by the Gene Expression Omnibus database, and quantitative real-time PCR in 113 pairs of CRC tissues and colon cancer cell lines. This study provides a pseudogene-associated ceRNA network, 7 prognostic pseudogene biomarkers, and a 5-pseudogene prognostic risk signature that may be useful for predicting the survival of CRC patients. [ABSTRACT FROM AUTHOR]

Background: Mycobacteroides abscessus complex (MABC), an emerging pathogen, causes human infections resistant to multiple antibiotics. In this study, the genome data of 1,581 MABC strains were downloaded from NCBI database for phylogenetic relatedness inference, resistance profile identification and the estimation of evolutionary pressure on resistance genes in silico. Results: From genes associated with resistance to 28 antibiotic classes, 395 putative proteins (ARPs) were identified, based on the information in two antibiotic resistance databases (CARD and ARG-ANNOT). The ARPs most frequently identified in MABC were those associated with resistance to multiple antibiotic classes, beta-lactams and aminoglycosides. After excluding ARPs that had undergone recombination, two ARPs were predicted to be under diversifying selection and 202 under purifying selection. This wide occurrence of purifying selection suggested that the diversity of commonly shared ARPs in MABC have been reduced to achieve stability. The unequal distribution of ARPs in members of the MABC could be due to horizontal gene transfer or ARPs pseudogenization events. Most (81.5%) of the ARPs were observed in the accessory genome and 72.2% ARPs were highly homologous to proteins associated with mobile genetic elements such as plasmids, prophages and viruses. On the other hand, with TBLASTN search, only 18 of the ARPs were identified as pseudogenes. Conclusion: Altogether, our results suggested an important role of horizontal gene transfer in shaping the resistome of MABC. [ABSTRACT FROM AUTHOR]

Recurrent epidemics of cholera denote robust adaptive mechanisms of Vibrio cholerae for ecological shifting and persistence despite variable stress conditions. Tracking the evolution of pathobiological traits requires comparative genomic studies of isolates from endemic areas. Here, we investigated the genetic differentiation among V. cholerae clinical and environmental isolates by highlighting the genomic divergence associated with gene decay, genome plasticity, and the acquisition of virulence and adaptive traits. The clinical isolates showed high phylogenetic relatedness due to a higher frequency of shared orthologs and fewer gene variants in contrast to the evolutionarily divergent environmental strains. Divergence of the environmental isolates is linked to extensive genomic rearrangements in regions containing mobile genetic elements resulting in numerous breakpoints, relocations, and insertions coupled with the loss of virulence determinants acf, zot, tcp, and ctx in the genomic islands. Also, four isolates possessed the CRISPR-Cas systems with spacers specific for Vibrio phages and plasmids. Genome synteny and homology analysis of the CRISPRCas systems suggest horizontal acquisition. The marked differences in the distribution of other phage and plasmid defense systems such as Zorya, DdmABC, DdmDE, and type-I Restriction Modification systems among the isolates indicated a higher propensity for plasmid or phage disseminated traits in the environmental isolates. Our results reveal that V. cholerae strains undergo extensive genomic rearrangements coupled with gene acquisition, reflecting their adaptation during ecological shifts and pathogenicity. [ABSTRACT FROM AUTHOR]

To investigate the role of Pseudogenes (PG) in juvenile nasopharyngeal angiofibroma (JNA) that were once considered to be junk-DNA or 'genomic fossils'. Five to 10 fresh JNA samples were analyzed for molecular expressions of 5 PG/respective parent genes (VEGFR1P1/VEGFR; FGFR3P1/FGFR3; PDGFAP1/PDGFA; IL6RP1/IL6R; POU5F1B/POU5F1) and clinical details correlated. IL6R, PDGFA, VEGFR2, FGF3 and their respective PG (VEGFR1P1, PDGFAP1, IL6RP1, FGFR3P1) were highly expressed but POU5F1/POU5F1B were not. The difference in expression between IL6R & IL6RP1 was substantially larger compared with other 4 pairs. While VEGFR1P1, PDGFAP1 & POU5F1B were expressed more than their respective parent genes; IL6RP1 & FGFR3P1 showed reduced expression. No clinical significance was apparent in any parameter nor was any difference noted between recurrent and upfront cases. A definite implication of pseudogene in JNA is evident in this ever first global study but future studies are needed to validate the current findings as well as further characterize its role/profile in larger sample. This may explain extreme variability of JNA, its heterogenous etiopathogenesis, evolving patterns and molecular characterization for possible targeted therapy. [ABSTRACT FROM AUTHOR]

Trypanosomatids belong to a remarkable group of unicellular, parasitic organisms of the order Kinetoplastida, an early diverging branch of the phylogenetic tree of eukaryotes, exhibiting intriguing biological characteristics affecting gene expression (intronless polycistronic transcription, trans-splicing, and RNA editing), metabolism, surface molecules, and organelles (compartmentalization of glycolysis, variation of the surface molecules, and unique mitochondrial DNA), cell biology and life cycle (phagocytic vacuoles evasion and intricate patterns of cell morphogenesis). With numerous genomic-scale data of several trypanosomatids becoming available since 2005 (genomes, transcriptomes, and proteomes), the scientific community can further investigate the mechanisms underlying these unusual features and address other unexplored phenomena possibly revealing biological aspects of the early evolution of eukaryotes. One fundamental aspect comprises the processes and mechanisms involved in the acquisition and loss of genes throughout the evolutionary history of these primitive microorganisms. Here, we present a comprehensive in silico analysis of pseudogenes in three major representatives of this group: Leishmania major , Trypanosoma brucei , and Trypanosoma cruzi. Pseudogenes, DNA segments originating from altered genes that lost their original function, are genomic relics that can offer an essential record of the evolutionary history of functional genes, as well as clues about the dynamics and evolution of hosting genomes. Scanning these genomes with functional proteins as proxies to reveal intergenic regions with protein-coding features, relying on a customized threshold to distinguish statistically and biologically significant sequence similarities, and reassembling remnant sequences from their debris, we found thousands of pseudogenes and hundreds of open reading frames, with particular characteristics in each trypanosomatid: mutation profile, number, content, density, codon bias, average size, single- or multi-copy gene origin, number and type of mutations, putative primitive function, and transcriptional activity. These features suggest a common process of pseudogene formation, different patterns of pseudogene evolution and extant biological functions, and/or distinct genome organization undertaken by those parasites during evolution, as well as different evolutionary and/or selective pressures acting on distinct lineages. [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]

The internal transcribed spacer (ITS) is one of the most extensively sequenced molecular markers in plant systematics due to its generally concerted evolution. While non-concerted evolution has been found in some plant taxa, such information is missing in Lycium. Molecular studies of six species and two variants of the genus Lycium revealed high levels of intra- and inter-individual polymorphism in the ITS, indicating non-concerted evolution. All genomic DNA ITS paralogues were identified as putative pseudogenes or functional paralogues through a series of comparisons of sequence features, including length and substitution variation, GC content, secondary structure stability, and the presence of conserved motifs in the 5.8S gene, and the rate of evolution. Approximately, 60% of ITS pseudogenes could be easily detected. Based on phylogenetic analysis, all pseudogenes were highly distinct fromtheir corresponding functional copies, tended to evolve neutrally, and clustered randomly together in the evolutionary tree. The results probably suggest that this ITS non-concerted evolution is related to the recent divergence between tandemrepeats within the Lyciumgenome and hybridization between species. Our study complements those of pseudogenes in plant taxa and provides a theoretical basis for the phylogeny and genetic origin of the genus Lycium while having important implications for the use of ITS molecular markers for phylogenetic reconstruction. [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]

Hevea pauciflora belongs to the Euphorbiaceae family, an important wild relative of the rubber tree. This study sequenced, assembled, and annotated the complete chloroplast genome of H. pauciflora. The complete chloroplast genome is 161,123 bp with a canonical quadripartite structure containing a large single-copy (LSC) region (89,109 bp), a small single-copy (SSC) region (18,376 bp), and two inverted repeat regions (IRa and IRb) (26,819 bp, each). A total of 134 genes were annotated, including 86 protein-coding genes, four pseudogenes, 36 tRNA genes, and eight rRNA genes. The 134 genes include four major groups: 'self-replication', 'photosynthesis', 'unknown function', and 'others'. A phylogenetic analysis clustered H. pauciflora, H. brasiliensis, H. camargoana, and H. benthamiana into one clade, consistent with traditional taxonomy. This study provides useful data for further studies of Hevea genus and the phylogenetic relationships of Euphorbiaceae species. [ABSTRACT FROM AUTHOR]

We report the finding of two copy number variants (CNVs) in a 12-year-old boy presenting both with autism spectrum disorder (ASD) and attention deficit/hyperactivity disorder (ADHD). Clinical features included aggressive behavior, mood instability, suicidal statements, repetitive and restrictive behavior, sensitivity to noise, learning problems and dyslexia, though no intellectual disability was present. Using array-based comparative genomic hybridization (array-CGH), we identified two CNVs, both triplex duplications of 324 kb on 3p26.3, and 284 kb on 4q13.1, respectively. One of the CNVs is located on chromosome 4q13.1 in the region of the gene encoding for adhesion G protein-coupled receptor L3 (ADGRL3, former name: latrophilin-3, LPHN3), the other on chromosome 3p26.3 in the region of the two pseudogenes AC090043.1 and RPL23AP39. The patient described in the present study showed increased symptoms under methylphenidate treatment but responded positively to 3 mg per day of the atypical neuroleptic drug aripiprazole. To our knowledge, this is the first report of a CNV in the ADGRL3 gene and its first association with ASD in humans. [ABSTRACT FROM AUTHOR]

Abstract Background Naked mole-rats (Heterocephalus glaber, NMRs) and blind mole-rats (Spalax galili, BMRs) are representative subterranean rodents that have evolved many extraordinary traits, including hypoxia tolerance, longevity, and cancer resistance. Although multiple candidate loci responsible for these traits have been uncovered by genomic studies, many of them are limited to functional changes to amino acid sequence and little is known about the contributions of other genetic events. To address this issue, we focused on gene losses (unitary pseudogenes) and systematically analyzed gene losses in NMRs and BMRs, aiming to elucidate the potential roles of pseudogenes in their adaptation to subterranean lifestyle. Results We obtained the pseudogene repertoires in NMRs and BMRs, as well as their respective aboveground relatives, guinea pigs and rats, on a genome-wide scale. As a result, 167, 139, 341, and 112 pseudogenes were identified in NMRs, BMRs, guinea pigs, and rats, respectively. Functional enrichment analysis identified 4 shared and 2 species-specific enriched functional groups (EFGs) in subterranean lineages. Notably, the pseudogenes in these EFGs might be associated with either regressive (e.g., visual system) or adaptive (e.g., altered DNA damage response) traits. In addition, several pseudogenes including TNNI3K and PDE5A might be associated with specific cardiac features observed in subterranean lineages. Interestingly, we observed 20 convergent gene losses in NMRs and BMRs. Given that the functional investigations of these genes are generally scarce, we provided functional evidence that independent loss of TRIM17 in NMRs and BMRs might be beneficial for neuronal survival under hypoxia, supporting the positive role of eliminating TRIM17 function in hypoxia adaptation. Our results also suggested that pseudogenes, together with positively selected genes, reinforced subterranean adaptations cooperatively. Conclusions Our study provides new insights into the molecular underpinnings of subterranean adaptations and highlights the importance of gene losses in mammalian evolution.

Avian sex chromosomes evolved after the divergence of birds and crocodilians from their common ancestor, so they are younger than the better-studied chromosomes of mammals. It has long been recognized that there may have been several stages to the evolution of avian sex chromosomes. For example, the CHD1 undergoes recombination in paleognaths but not neognaths. Genome assemblies have suggested that there may be variation in the timing of barriers to recombination among Neognathae, but there remains little understanding of the extent of this variability. Here, we look at partial sequences of ATP5F1A, which is on the avian Z and W chromosomes. It is known that recombination of this gene has independently ceased in Galliformes, Anseriformes, and at least five neoavian orders, but whether there are other independent cessations of recombination among Neoaves is not understood. We analyzed a combination of data extracted from published chromosomal-level genomes with data collected using PCR and cloning to identify Z and W copies in 22 orders. Our results suggest that there may be at least 19 independent cessations of recombination within Neognathae, and 3 clades that may still be undergoing recombination (or have only recently ceased recombination). Analyses of ATP5F1A protein sequences revealed an increased amino acid substitution rate for W chromosome gametologs, suggesting relaxed purifying selection on the W chromosome. Supporting this hypothesis, we found that the increased substitution rate was particularly pronounced for buried residues, which are expected to be more strongly constrained by purifying selection. This highlights the dynamic nature of avian sex chromosomes, and that this level of variation among clades means they should be a good system to understand sex chromosome evolution.

Statins are the most commonly prescribed cardiovascular disease drug, but their inter-individual efficacy varies considerably. Genetic factors uncovered to date have only explained a small proportion of variation in low-density lipoprotein cholesterol (LDLC) lowering. To identify novel markers and determinants of statin response, we used whole transcriptome sequence data collected from simvastatin and control incubated lymphoblastoid cell lines (LCLs) established from participants of the Cholesterol and Pharmacogenetics (CAP) simvastatin clinical trial. We looked for genes whose statin-induced expression changes were most different between LCLs derived from individuals with high versus low plasma LDLC statin response during the CAP trial. We created a classification model of 82 "signature" gene expression changes that distinguished high versus low LDLC statin response. One of the most differentially changing genes was zinc finger protein 542 pseudogene (ZNF542P), the signature gene with changes most correlated with statin-induced change in cellular cholesterol ester, an in vitro marker of statin response. ZNF542P knock-down in a human hepatoma cell line increased intracellular cholesterol ester levels upon simvastatin treatment. Together, these findings imply a role for ZNF542P in LDLC response to simvastatin and, importantly, highlight the potential significance of noncoding RNAs as a contributing factor to variation in drug response.

Loblolly pine (LP; Pinus taeda L.) is an economically and ecologically important tree in the southeastern U.S. To advance understanding of the loblolly pine (LP; Pinus taeda L.) genome, we sequenced and analyzed 100 BAC clones and performed a Cot analysis. The Cot analysis indicates that the genome is composed of 57, 24, and 10% highly-repetitive, moderately-repetitive, and single/low-copy sequences, respectively (the remaining 9% of the genome is a combination of fold back and damaged DNA). Although single/low-copy DNA only accounts for 10% of the LP genome, the amount of single/low-copy DNA in LP is still 14 times the size of the Arabidopsis genome. Since gene numbers in LP are similar to those in Arabidopsis, much of the single/low-copy DNA of LP would appear to be composed of DNA that is both gene- and repeat-poor. Macroarrays prepared from a LP bacterial artificial chromosome (BAC) library were hybridized with probes designed from cell wall synthesis/wood development cDNAs, and 50 of the "targeted" clones were selected for further analysis. An additional 25 clones were selected because they contained few repeats, while 25 more clones were selected at random. The 100 BAC clones were Sanger sequenced and assembled. Of the targeted BACs, 80% contained all or part of the cDNA used to target them. One targeted BAC was found to contain fungal DNA and was eliminated from further analysis. Combinations of similarity-based and ab initio gene prediction approaches were utilized to identify and characterize potential coding regions in the 99 BACs containing LP DNA. From this analysis, we identified 154 gene models (GMs) representing both putative protein-coding genes and likely pseudogenes. Ten of the GMs (all of which were specifically targeted) had enough support to be classified as intact genes. Interestingly, the 154 GMs had statistically indistinguishable (α = 0.05) distributions in the targeted and random BAC clones (15.18 and 12.61 GM/Mb, respectively), whereas the low-repeat BACs contained significantly fewer GMs (7.08 GM/Mb). However, when GM length was considered, the targeted BACs had a significantly greater percentage of their length in GMs (3.26%) when compared to random (1.63%) and low-repeat (0.62%) BACs. The results of our study provide insight into LP evolution and inform ongoing efforts to produce a reference genome sequence for LP, while characterization of genes involved in cell wall production highlights carbon metabolism pathways that can be leveraged for increasing wood production.

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.

Recurrent epidemics of cholera denote robust adaptive mechanisms of Vibrio cholerae for ecological shifting and persistence despite variable stress conditions. Tracking the evolution of pathobiological traits requires comparative genomic studies of isolates from endemic areas. Here, we investigated the genetic differentiation among V. cholerae clinical and environmental isolates by highlighting the genomic divergence associated with gene decay, genome plasticity, and the acquisition of virulence and adaptive traits. The clinical isolates showed high phylogenetic relatedness due to a higher frequency of shared orthologs and fewer gene variants in contrast to the evolutionarily divergent environmental strains. Divergence of the environmental isolates is linked to extensive genomic rearrangements in regions containing mobile genetic elements resulting in numerous breakpoints, relocations, and insertions coupled with the loss of virulence determinants acf, zot, tcp, and ctx in the genomic islands. Also, four isolates possessed the CRISPR-Cas systems with spacers specific for Vibrio phages and plasmids. Genome synteny and homology analysis of the CRISPR-Cas systems suggest horizontal acquisition. The marked differences in the distribution of other phage and plasmid defense systems such as Zorya, DdmABC, DdmDE, and type-I Restriction Modification systems among the isolates indicated a higher propensity for plasmid or phage disseminated traits in the environmental isolates. Our results reveal that V. cholerae strains undergo extensive genomic rearrangements coupled with gene acquisition, reflecting their adaptation during ecological shifts and pathogenicity.

ABSTRACT In Colombia Aedes aegypti is present in 80% of the country up to 2,300 m; however, little is known of its genetic relations within a country context and, hence, within a broader context, for example, America. The aforementioned, herein, analyzed the gene flow within a context of the Americas, its directionality and genetic diversity of the mitochondrial lineages in the A. aegypti populations for Colombia. This called for the use of the sequences for A. aegypti available of the mitochondrial ND4 gene in the GenBank for Colombia and the American continent. No presence was found of nuclear mitochondrial pseudogenes (NUMTs) for Colombia. It is estimated that in Colombia the gene flow of the A. aegypti populations is occurring from the southeast and northeast toward the center of the country. In comparison with the mitochondrial sequences for America, the vector’s haplotypes in Colombia suggest connections between the populations of mosquitoes from the south with those from the north of the continent. The gene flow model at continental scale suggests bidirectional connections between the populations from the north of the continent with those from the south, while at South American scale it proposes the gene flow in all the directions with respect to the Colombian. The Colombian A. aegypti vector monitoring and control strategies must pay special attention to the vector’s points of entry into Colombia related with Peru, Venezuela, Brazil, Mexico, and North America to avoid the entry of populations with characteristics like resistance to insecticides or vector competition.

BackgroundAnalyzing close species with diverse developmental modes is instrumental for investigating the evolutionary significance of physiological, anatomical and behavioral features at a molecular level. Many examples of trait loss are known in metazoan populations living in dark environments. Tunicates are the closest living relatives of vertebrates and typically present a lifecycle with distinct motile larval and sessile adult stages. The nervous system of the motile larva contains melanized cells associated with geotactic and light-sensing organs. It has been suggested that these are homologous to vertebrate neural crest-derived melanocytes. Probably due to ecological adaptation to distinct habitats, several species of tunicates in the Molgulidae family have tailless (anural) larvae that fail to develop sensory organ-associated melanocytes. Here we studied the evolution of Tyrosinase family genes, indispensible for melanogenesis, in the anural, unpigmented Molgula occulta and in the tailed, pigmented Molgula oculata by using phylogenetic, developmental and molecular approaches.ResultsWe performed an evolutionary reconstruction of the tunicate Tyrosinase gene family: in particular, we found that M. oculata possesses genes predicted to encode one Tyrosinase (Tyr) and three Tyrosinase-related proteins (Tyrps) while M. occulta has only Tyr and Tyrp.a pseudogenes that are not likely to encode functional proteins. Analysis of Tyr sequences from various M. occulta individuals indicates that different alleles independently acquired frameshifting short indels and/or larger mobile genetic element insertions, resulting in pseudogenization of the Tyr locus. In M. oculata, Tyr is expressed in presumptive pigment cell precursors as in the model tunicate Ciona robusta. Furthermore, a M. oculata Tyr reporter gene construct was active in the pigment cell precursors of C. robusta embryos, hinting at conservation of the regulatory network underlying Tyr expression in tunicates. In contrast, we did not observe any expression of the Tyr pseudogene in M. occulta embryos. Similarly, M. occulta Tyr allele expression was not rescued in pigmented interspecific M. occulta × M. oculata hybrid embryos, suggesting deleterious mutations also to its cis-regulatory sequences. However, in situ hybridization for transcripts from the M. occulta Tyrp.a pseudogene revealed its expression in vestigial pigment cell precursors in this species.ConclusionsWe reveal a complex evolutionary history of the melanogenesis pathway in tunicates, characterized by distinct gene duplication and loss events. Our expression and molecular data support a tight correlation between pseudogenization of Tyrosinase family members and the absence of pigmentation in the immotile larvae of M. occulta. These results suggest that relaxation of purifying selection has resulted in the loss of sensory organ-associated melanocytes and core genes in the melanogenesis biosynthetic pathway in M. occulta.

BackgroundGenome duplication has played a pivotal role in the evolution of many eukaryotic lineages, including the vertebrates. A relatively recent vertebrate genome duplication is that in Xenopus laevis, which resulted from the hybridization of two closely related species about 17 million years ago. However, little is known about the consequences of this duplication at the level of the genome, the epigenome, and gene expression.ResultsThe X. laevis genome consists of two subgenomes, referred to as L (long chromosomes) and S (short chromosomes), that originated from distinct diploid progenitors. Of the parental subgenomes, S chromosomes have degraded faster than L chromosomes from the point of genome duplication until the present day. Deletions appear to have the largest effect on pseudogene formation and loss of regulatory regions. Deleted regions are enriched for long DNA repeats and the flanking regions have high alignment scores, suggesting that non-allelic homologous recombination has played a significant role in the loss of DNA. To assess innovations in the X. laevis subgenomes we examined p300-bound enhancer peaks that are unique to one subgenome and absent from X. tropicalis. A large majority of new enhancers comprise transposable elements. Finally, to dissect early and late events following interspecific hybridization, we examined the epigenome and the enhancer landscape in X. tropicalis × X. laevis hybrid embryos. Strikingly, young X. tropicalis DNA transposons are derepressed and recruit p300 in hybrid embryos.ConclusionsThe results show that erosion of X. laevis genes and functional regulatory elements is associated with repeats and non-allelic homologous recombination and furthermore that young repeats have also contributed to the p300-bound regulatory landscape following hybridization and whole-genome duplication.

A group of aggressive pathogens have evolved to colonize the plant xylem. In this vascular tissue, where water and nutrients are transported from the roots to the rest of the plant, pathogens must be able to thrive under acropetal xylem sap flow and scarcity of nutrients while having direct contact only with predominantly dead cells. Nevertheless, a few bacteria have adapted to exclusively live in the xylem, and various pathogens may colonize other plant niches without causing symptoms unless they reach the xylem. Once established, the pathogens modulate its physicochemical conditions to enhance their growth and virulence. Adaptation to the restrictive lifestyle of the xylem leads to genome reduction in xylem-restricted bacteria, as they have a higher proportion of pseudogenes in their genome. The basis of xylem adaptation is not completely understood; therefore, a need still exists for model systems to advance the knowledge on this topic. [ABSTRACT FROM AUTHOR]

Nuclear mitochondrial pseudogenes (numts) are short, nonfunctional DNA sequences resulting from mitochondrial gene transfer to the nuclear genome; these sequences seriously affect the accuracy of mitochondrial gene for population genetic and phylogenetic analyses. The nuclear genome of Dendroctonus valens (Coleoptera: Curculionidae: Scolytidae), a major forest quarantine pest that seriously threatens the ecological security of pine forests, contains a large number of numts, and these interfere with analyses of genetic variation based on the mitochondrial cox1 gene. To eliminate interference by numts, we sequenced, assembled and analysed the complete mitochondrial genome of D. valens. The mitogenome is a circular molecule of 16,547 bp and shows typical gene content and arrangement. We identified its codon preference, tRNA secondary structure and nucleotide evolution rate. Furthermore, we designed primers for long fragment amplification to eliminate numts in the mitochondrial cox1 gene. Finally, amplified products were sequenced with universal primers, and cox1 gene sequences (464 bp) without numts were successfully obtained. The Heilihe Nature Reserve (harbouring Pinus tabulaeformis) in Inner Mongolia, is in the early stage of D. valens northward invasion, with a high population genetic diversity. We hypothesized that the high genetic diversity was related to human activities around the reserve and the reproductive characteristics in the early stage of D. valens invasion. The newly developed method resolves the issue of numt interference and lays a foundation for the analyses of genetic diversity in D. valens populations in the region as well as its invasion routes. [ABSTRACT FROM AUTHOR]

Type 2 diabetes (T2D) is a complex metabolic derangement that has a strong genetic basis. There is substantial population-specificity in the association of genetic variants with T2D. The Indian urban Sindhi population is at a high risk of T2D. The genetic basis of T2D in this population is unknown. We interrogated 28 pooled whole blood genomes of 1402 participants from the Diabetes In Sindhi Families In Nagpur (DISFIN) study using Illumina's Global Screening Array. From a total of 608,550 biallelic variants, 140 were significantly associated with T2D after adjusting for comorbidities, batch effects, pooling error, kinship status and pooling variation in a random effects multivariable logistic regression framework. Of the 102 well-characterized genes that these variants mapped onto, 70 genes have been previously reported to be associated with T2D to varying degrees with known functional relevance. Excluding open reading frames, intergenic non-coding elements and pseudogenes, our study identified 22 novel candidate genes in the Sindhi population studied. Our study thus points to the potential, interesting candidate genes associated with T2D in an ethnically endogamous population. These candidate genes need to be fully investigated in future studies. [ABSTRACT FROM AUTHOR]

tRNA, as well as their derived products such as short interspersed nuclear elements (SINEs), pseudogenes, and transfer RNA (tRNA)-derived fragments (tRFs), have now been shown to be vital for cellular life, functioning, and adaptation during different stress conditions in all diverse life forms. In this study, we have developed PtRNAdb (www.nipgr.ac.in/PtRNAdb), a plant-exclusive tRNA database containing 113,849 tRNA gene sequences from phylogenetically diverse plant species. We have analyzed a total of 106 nuclear, 89 plastidial, and 38 mitochondrial genomes of plants by the tRNAscan-SE software package, and after careful curation of the output data, we integrated the data and developed this database. The information about the tRNA gene sequences obtained was further enriched with a consensus sequence-based study of tRNA genes based on their isoacceptors and isodecoders. We have also built covariance models based on the isoacceptors and isodecoders of all the tRNA sequences using the infernal tool. The user can also perform BLAST not only against PtRNAdb entries but also against all the tRNA sequences stored in the PlantRNA database and annotated tRNA genes across the plant kingdom available at NCBI. This resource is believed to be of high utility for plant researchers as well as molecular biologists to carry out further exploration of the plant tRNAome on a wider spectrum, as well as for performing comparative and evolutionary studies related to tRNAs, and their derivatives across all domains of life. [ABSTRACT FROM AUTHOR]

Pseudogenes have been considered as non-functional genes. However, peptides and long non-coding RNAs produced by pseudogenes are expressed in different tumors. Moreover, the dysregulation of pseudogenes is associated with cancer, and their expressions are higher in tumors compared to normal tissues. Recent studies show that pseudogenes can influence the liquid phase condensates formation. Liquid phase separation involves regulating different epigenetic stages, including transcription, chromatin organization, 3D DNA structure, splicing, and post-transcription modifications like m6 A. Several membrane-less organelles, formed through the liquid phase separate, are also involved in the epigenetic regulation, and their defects are associated with cancer development. However, the association between pseudogenes and liquid phase separation remains unrevealed. The current study sought to investigate the relationship between pseudogenes and liquid phase separation in cancer development, as well as their therapeutic implications. [ABSTRACT FROM AUTHOR]

This is a response to a preprint version of “A re-analysis of the data in Sharkey et al.’s (2021) minimalist revision reveals that BINs do not deserve names, but BOLD Systems needs a stronger commitment to open science”, https://www.biorxiv.org/content/10.1101/2021.04.28.441626v2. Meier et al. strongly criticized Sharkey et al.’s publication in which 403 new species were deliberately minimally described, based primarily on COI barcode sequence data. Here we respond to these criticisms. The following points are made: 1) Sharkey et al. did not equate BINs with species, as demonstrated in several examples in which multiple species were found to be in single BINs. 2) We reiterate that BINs were used as a preliminary sorting tool, just as preliminary morphological identification commonly sorts specimens based on color and size into unit trays; despite BINs and species concepts matching well over 90% of species, this matching does not equate to equality. 3) Consensus barcodes were used only to provide a diagnosis to conform to the rules of the International Code of Zoological Nomenclature just as consensus morphological diagnoses are. The barcode of a holotype is definitive and simply part of its cellular morphology. 4) Minimalist revisions will facilitate and accelerate future taxonomic research, not hinder it. 5) We refute the claim that the BOLD sequences of Plesiocoelus vanachterbergi are pseudogenes and demonstrate that they simply represent a frameshift mutation. 6) We reassert our observation that morphological evidence alone is insufficient to recognize species within species-rich higher taxa and that its usefulness lies in character states that are congruent with molecular data. 7) We show that in the cases in which COI barcodes code for the same amino acids in different putative species, data from morphology, host specificity, and other ecological traits reaffirm their utility as indicators of genetically distinct lineages. [ABSTRACT FROM AUTHOR]

DNA barcoding is based on the premise that the barcode sequences can distinguish individuals (strains) of different species because their sequence variation between species exceeds that within species. The primary barcodes used in fungal and yeast taxonomy are the ITS segments and the LSU (large subunit) D1/D2 domain of the homogenized multicopy rDNA repeats. The secondary barcodes are conserved segments of protein‐encoding genes, which usually have single copies in haploid genomes. This study shows that the analysis of barcode sequences fails to reconstruct accurate species trees and differentiate species when the organisms have chimeric genomes composed of admixed mosaics of different origins. It is shown that the type strains of 10 species of the pulcherrima clade of the ascomycetous yeast genus Metschnikowia cannot be differentiated with standard barcodes because their intragenomic diversity is comparable to or even higher than the interstrain diversity. The analysis of a large group of genes of the sequenced genomes of the clade and the viability and segregation of the hybrids of ex‐type strains indicate that the high intragenomic barcode differences can be attributed to admixed genome structures. Because of the mosaic structures of the genomes, the rDNA repeats do not form continuous arrays and thus cannot be homogenized. Since the highly diverse ITS and D1/D2 sequences of the type strains form a continuous pool including pseudogenes, the evolution of their rDNA appears to involve reticulation. The secondary barcode sequences and the nonbarcode genes included in the analysis show incongruent phylogenetic relationships among the type strains, which can also be attributed to differences in the phylogenetic histories of the genes. [ABSTRACT FROM AUTHOR]

Apart from I GBA i and I PRKRA i , a potential role of pseudogenes has not (yet) caught our attention in the field of movement disorders, although several well-established disease-linked genes also have pseudogenes. The definition of a pseudogene has now broadened to include any DNA sequence that is similar to a known gene and has lost some of its original functionality.25 Therefore, pseudogenes can be translated into proteins, as recently proven for 140 human pseudogenes, and also act at the protein level.26 For instance, pseudogene-derived proteins may have the same activity as the parental proteins, but function in different tissues, cellular compartments, or pathophysiological conditions. Pseudogenes can be divided into two major groups, that is, unprocessed pseudogenes and processed pseudogenes, with the latter representing the majority of pseudogenes in humans (~70%, https://www.gencodegenes.org/human/stats.html). Therefore, the sequence identity for these pseudogenes should be too less to seriously trouble sequence analyses.27 Overall, pseudogenes have emerged as a hot topic because they can challenge genetic testing and trigger mutational events even in biologically functional parental genes. [Extracted from the article]

Prokaryotic genomes are usually densely packed with intact and functional genes. However, in certain contexts, such as after recent ecological shifts or extreme population bottlenecks, broken and nonfunctional gene fragments can quickly accumulate and form a substantial fraction of the genome. Identification of these broken genes, called pseudogenes, is a critical step for understanding the evolutionary forces acting upon, and the functional potential encoded within, prokaryotic genomes. Here, we present Pseudofinder, an open-source software dedicated to pseudogene identification and analysis in bacterial and archaeal genomes. We demonstrate that Pseudofinder's multi-pronged, reference-based approach can detect a wide variety of pseudogenes, including those that are highly degraded and typically missed by gene-calling pipelines, as well newly formed pseudogenes containing only one or a few inactivating mutations. Additionally, Pseudofinder can detect genes that lack inactivating substitutions but experiencing relaxed selection. Implementation of Pseudofinder in annotation pipelines will allow more precise estimations of the functional potential of sequenced microbes, while also generating new hypotheses related to the evolutionary dynamics of bacterial and archaeal genomes. [ABSTRACT FROM AUTHOR]

Objective. Alterations in the methylation state of pseudogenes may serve as clinically useful biomarkers of glioblastomas (GBMs) that do not have glioma-CpG island methylator phenotype (G-CIMP). Methods. Non-G-CIMP GBM datasets were included for evaluation, and a RISK-score signature was determined from the methylation state of pseudogene loci. Both bioinformatic and experimental analyses were performed for biological validation. Results. By integrating clinical information with DNA methylation microarray data, we screened a panel of eight CpGs from discovery cohorts of non-G-CIMP GBMs. Each CpG could accurately and independently predict the prognosis of patients under a treatment regime that combined radiotherapy (RT) and temozolomide (TMZ). The 8-CpG signature appeared to show opposite prognostic correlations between patients treated with RT/TMZ and those treated with RT monotherapy. The analyses further indicated that this signature had predictive value for TMZ efficacy because different survival benefits between RT/TMZ and RT therapies were observed in each risk subgroup. The incorporation of other risk factors, such as age and O-6-methylguanine-DNA methyltransferase (MGMT) promoter methylation status, with our pseudogene methylation signature could provide precise risk classification. In vitro experimental data revealed that two locus-specific pseudogenes (ZNF767P and CLEC4GP1) may modulate TMZ resistance via distinct mechanisms in GBM cells. Conclusion. The biologically and clinically relevant RISK-score signature, based on pseudogene methylation loci, may offer information for predicting TMZ responses of non-G-CIMP GBMs, that is independent from, but complementary to, MGMT-based approaches. [ABSTRACT FROM AUTHOR]

CRC is the third most common cancer occurring worldwide and the second leading cause of cancer deaths. In the year 2020, 1,931,590 new cases of CRC and 935,173 deaths were reported. The last two decades have witnessed an intensive study of noncoding RNAs and their implications in various pathological conditions including cancer. Noncoding RNAs such as miRNAs, tsRNAs, piRNAs, lncRNAs, pseudogenes, and circRNAs have emerged as promising prognostic and diagnostic biomarkers in preclinical studies of cancer. Some of these noncoding RNAs have also been shown as promising therapeutic targets for cancer treatment. In this review, we have discussed the emerging roles of various types of noncoding RNAs in CRC and their future implications in colorectal cancer management and research. [ABSTRACT FROM AUTHOR]

Numerous animal lineages have expanded and diversified the opsin-based photoreceptors in their eyes underlying color vision behavior. However, the selective pressures giving rise to new photoreceptors and their spectral tuning remain mostly obscure. Previously, we identified a violet receptor (UV2) that is the result of a UV opsin gene duplication specific to Heliconius butterflies. At the same time the violet receptor evolved, Heliconius evolved UV-yellow coloration on their wings, due to the pigment 3-hydroxykynurenine (3-OHK) and the nanostructure architecture of the scale cells. In order to better understand the selective pressures giving rise to the violet receptor, we characterized opsin expression patterns using immunostaining (14 species) and RNA-Seq (18 species), and reconstructed evolutionary histories of visual traits in five major lineages within Heliconius and one species from the genus Eueides. Opsin expression patterns are hyperdiverse within Heliconius. We identified six unique retinal mosaics and three distinct forms of sexual dimorphism based on ommatidial types within the genus Heliconius. Additionally, phylogenetic analysis revealed independent losses of opsin expression, pseudogenization events, and relaxation of selection on UVRh2 in one lineage. Despite this diversity, the newly evolved violet receptor is retained across most species and sexes surveyed. Discriminability modeling of behaviorally preferred 3-OHK yellow wing coloration suggests that the violet receptor may facilitate Heliconius color vision in the context of conspecific recognition. Our observations give insights into the selective pressures underlying the origins of new visual receptors.

Background: China is the hotspot of global freshwater crab diversity, but their wild populations are facing severe pressures associated with anthropogenic factors, necessitating the need to map their taxonomic and genetic diversity and design conservation policies., Results: Herein, we sequenced the mitochondrial genome of a Chinese freshwater crab species Bottapotamon fukienense, and found that it is fragmented into two chromosomes. We confirmed that fragmentation was not limited to a single specimen or population. Chromosome 1 comprised 15,111 base pairs (bp) and there were 26 genes and one pseudogene (pseudo-nad1) encoded on it. Chromosome 2 comprised 8,173 bp and there were 12 genes and two pseudogenes (pseudo-trnL2 and pseudo-rrnL) encoded on it. Combined, they comprise the largest mitogenome (23,284 bp) among the Potamidae. Bottapotamon was the only genus in the Potamidae dataset exhibiting rearrangements of protein-coding genes. Bottapotamon fukienense exhibited average rates of sequence evolution in the dataset and did not differ in selection pressures from the remaining Potamidae., Conclusions: This is the first experimentally confirmed fragmentation of a mitogenome in crustaceans. While the mitogenome of B. fukienense exhibited multiple signs of elevated mitogenomic architecture evolution rates, including the exceptionally large size, duplicated genes, pseudogenisation, rearrangements of protein-coding genes, and fragmentation, there is no evidence that this is matched by elevated sequence evolutionary rates or changes in selection pressures., (© 2024. The Author(s).)

Recently, substantial evidence has demonstrated that pseudogene-derived long noncoding RNAs (lncRNAs) as regulatory RNAs have been implicated in basic physiological processes and disease development through multiple modes of functional interaction with DNA, RNA, and proteins. Here, we report an important role for GBP1P1, the pseudogene of guanylate-binding protein 1, in regulating influenza A virus (IAV) replication in A549 cells. GBP1P1 was dramatically upregulated after IAV infection, which is controlled by JAK/STAT signaling. Functionally, ectopic expression of GBP1P1 in A549 cells resulted in significant suppression of IAV replication. Conversely, silencing GBP1P1 facilitated IAV replication and virus production, suggesting that GBP1P1 is one of the interferon-inducible antiviral effectors. Mechanistically, GBP1P1 is localized in the cytoplasm and functions as a sponge to trap DHX9 (DExH-box helicase 9), which subsequently restricts IAV replication. Together, these studies demonstrate that GBP1P1 plays an important role in antagonizing IAV replication.IMPORTANCELong noncoding RNAs (lncRNAs) are extensively expressed in mammalian cells and play a crucial role as regulators in various biological processes. A growing body of evidence suggests that host-encoded lncRNAs are important regulators involved in host-virus interactions. Here, we define a novel function of GBP1P1 as a decoy to compete with viral mRNAs for DHX9 binding. We demonstrate that GBP1P1 induction by IAV is mediated by JAK/STAT activation. In addition, GBP1P1 has the ability to inhibit IAV replication. Importantly, we reveal that GBP1P1 acts as a decoy to bind and titrate DHX9 away from viral mRNAs, thereby attenuating virus production. This study provides new insight into the role of a previously uncharacterized GBP1P1, a pseudogene-derived lncRNA, in the host antiviral process and a further understanding of the complex GBP network., Competing Interests: The authors declare no conflict of interest.

Pseudogenes have been considered as non-functional genes. However, peptides and long non-coding RNAs produced by pseudogenes are expressed in different tumors. Moreover, the dysregulation of pseudogenes is associated with cancer, and their expressions are higher in tumors compared to normal tissues. Recent studies show that pseudogenes can influence the liquid phase condensates formation. Liquid phase separation involves regulating different epigenetic stages, including transcription, chromatin organization, 3D DNA structure, splicing, and post-transcription modifications like m6A. Several membrane-less organelles, formed through the liquid phase separate, are also involved in the epigenetic regulation, and their defects are associated with cancer development. However, the association between pseudogenes and liquid phase separation remains unrevealed. The current study sought to investigate the relationship between pseudogenes and liquid phase separation in cancer development, as well as their therapeutic implications.

Mitochondrial DNA (mtDNA) translocated into the nuclear genome (numt), when co-analysed with genuine mtDNA, could plague phylogeographic studies. To evaluate numt-related biases in population genetics parameters in birds, which are prone to accumulating numts, we targeted the mitochondrial mt-cytb gene. We looked at 13 populations of Audubon's shearwater (Puffinus lherminieri), including five mitochondrial lineages. mt-cytb homologue and paralogue (numt) sequences were determined by Sanger sequencing with and without prior exonuclease digestion of nuclear DNA. Numts formed monophyletic clades corresponding to three of the five mitochondrial lineages tested (the remaining two forming a paraphyletic group). Nineteen percent of numt alleles fell outside of their expected mitochondrial clade, a pattern consistent with multiple translocation events, incomplete lineage sorting (ILS), and/or introgression. When co-analysing mt-cytb paralogues and homologues, excluding individuals with ambiguities underestimates genetic diversity (4%) and differentiation (11%) among least-sampled populations. Removing ambiguous sites drops the proportion of inter-lineage genetic variance by 63%. While co-analysing numts with mitochondrial sequences can lead to severe bias and information loss in bird phylogeographic studies, the separate analysis of genuine mitochondrial loci and their nuclear paralogues can shed light on numt molecular evolution, as well as evolutionary processes such as ILS and introgression.

Background: The role of copy number variants (CNVs) in susceptibility to asthma is not well understood. This is, in part, due to the difficulty of accurately measuring CNVs in large enough sample sizes to detect associations. The recent availability of whole-exome sequencing (WES) in large biobank studies provides an unprecedented opportunity to study the role of CNVs in asthma. Methods: We called common CNVs in 49,953 individuals in the first release of UK Biobank WES using ClinCNV software. CNVs were tested for association with asthma in a stage 1 analysis comprising 7098 asthma cases and 36,578 controls from the first release of sequencing data. Nominally-associated CNVs were then meta-analysed in stage 2 with an additional 17,280 asthma cases and 115,562 controls from the second release of UK Biobank exome sequencing, followed by validation and fine-mapping. Results: Five of 189 CNVs were associated with asthma in stage 2, including a deletion overlapping the HLA-DQA1 and HLA-DQB1 genes, a duplication of CHROMR/PRKRA, deletions within MUC22 and TAP2, and a duplication in FBRSL1. The HLA-DQA1, HLA-DQB1, MUC22 and TAP2 genes all reside within the human leukocyte antigen (HLA) region on chromosome 6. In silico analyses demonstrated that the deletion overlapping HLA-DQA1 and HLA-DQB1 is likely to be an artefact arising from under-mapping of reads from non-reference HLA haplotypes, and that the CHROMR/PRKRA and FBRSL1 duplications represent presence/absence of pseudogenes within the HLA region. Bayesian fine-mapping of the HLA region suggested that there are two independent asthma association signals. The variants with the largest posterior inclusion probability in the two credible sets were an amino acid change in HLA-DQB1 (glutamine to histidine at residue 253) and a multi-allelic amino acid change in HLA-DRB1 (presence/absence of serine, glycine or leucine at residue 11). Conclusions: At least two independent loci characterised by amino acid changes in the HLA-DQA1, HLA-DQB1 and HLA-DRB1 genes are likely to account for association of SNPs and CNVs in this region with asthma. The high divergence of haplotypes in the HLA can give rise to spurious CNVs, providing an important, cautionary tale for future large-scale analyses of sequencing data. [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 major histocompatibility complex (MHC) enables vertebrates to cope with pathogens and maintain healthy populations, thus making it a unique set of loci for addressing ecology and evolutionary biology questions. The aim of our study was to examine the variability of Heermann's Gull MHC class II (MHCIIB) and compare these loci with other Charadriiformes. Fifty-nine MHCIIB haplotypes were recovered from sixty-eight Heermann's Gulls by cloning, of them, twelve were identified as putative true alleles, forty-five as unique alleles, and two as pseudogenes. Intra and interspecific relationships indicated at least two loci in Heermann's Gull MHCIIB and trans-species polymorphism among Charadriiformes (coinciding with the documented evidence of two ancient avian MHCIIB lineages, except in the Charadriidae family). Additionally, sites under diversifying selection revealed a better match with peptide-binding sites inferred in birds than those described in humans. Despite the negative anthropogenic activity reported on Isla Rasa, Heermann's Gull showed MHCIIB variability consistent with population expansion, possibly due to a sudden growth following conservation efforts. Duplication must play an essential role in shaping Charadriiformes MHCIIB variability, buffering selective pressures through balancing selection. These findings suggest that MHC copy number and protected islands can contribute to seabird conservation. [ABSTRACT FROM AUTHOR]

The hypothesis that the evolution of humans involves hybridization between diverged species has been actively debated in recent years. We present the following novel evidence in support of this hypothesis: the analysis of nuclear pseudogenes of mtDNA ("NUMTs"). NUMTs are considered "mtDNA fossils" as they preserve sequences of ancient mtDNA and thus carry unique information about ancestral populations. Our comparison of a NUMT sequence shared by humans, chimpanzees, and gorillas with their mtDNAs implies that, around the time of divergence between humans and chimpanzees, our evolutionary history involved the interbreeding of individuals whose mtDNA had diverged as much as ~4.5 Myr prior. This large divergence suggests a distant interspecies hybridization. Additionally, analysis of two other NUMTs suggests that such events occur repeatedly. Our findings suggest a complex pattern of speciation in primate/human ancestors and provide one potential explanation for the mosaic nature of fossil morphology found at the emergence of the hominin lineage. A preliminary version of this manuscript was uploaded to the preprint server BioRxiv in 2017 (10.1101/134502). [ABSTRACT FROM AUTHOR]

HMGN2 have functions in inflammatory response. However, the role of HMGN2 in severe acute pancreatitis (SAP) remains unclear. Here, our study was to discuss the role and regulatory mechanism ofHMGN2 in SAP. In this study, the SAP cell model of AR42J was used to study the function and mechanism of HMGN2 in SAP. The protein expression in cells and serums were examined by western blot and ELISA assay. qPCR was used to test the transcriptional RNA level. Cell viability were examined by MTT assay. Luciferase assay was used to evaluate the interaction between gene and gene. Our results showed that HMGN2 was significantly upregulated in SAP patients. The database predicted and luciferase assay data indicated the HMGN2 was directly binding with miR-590-3p. ELISA, MTT and western blot experiments showed that the HMGN2 were promoted the cell proliferation, reduced the inflammation, and repressed the cell autophagy. Mechanism studies showed that the pseudogene HMGN2P46 level was positively correlated with HMGN2 and upregulated HMGN2 expression by competing for miR-590-3p in SAP. Taken together, all over these results showed upregulation of HMGN2 alleviates SAP, this process was regulated by HMGN2P46 competitively binding with miR-590-3p, which may provide a new insight for the treatment and intervention in SAP. Pseudogene HMGN2P46 was a miRNA sponge to regulate HMGN2 level by competing for miR-590-3p to alleviate the process of SAP. It provided a novel strategy for the diagnosis and treatment of severe acute pancreatitis. [ABSTRACT FROM AUTHOR]

Processed pseudogenes, also known as retrocopy genes, are copies of messenger RNAs that have been reverse transcribed into DNA and inserted into the genome. In this study, we identified 62 processed pseudogene candidates as intron‐less genes from whole‐genome sequencing (WGS) data of Thoroughbred horses using delly structural variation software. The 62 processed pseudogene candidates were confirmed by PCR amplification of intron‐less products. A total of 11 processed pseudogenes were confirmed in the genome of all 23 analysed horses, whereas three processed pseudogenes with structures of ATP11B, DPH3 and RPL17 were detected in only one of 115 horses by PCR amplification of intron‐less products. Currently, most of the gene doping tests proposed in human and horse sports are adapted PCR‐based methods using hydrolysis probes to detect exon/exon junctions in transgenes because the operation is simple and economical. However, when the pseudogene is present in the host genome, the PCR‐based methods may have a potential risk of detecting false positives. In this study, because processed pseudogenes that exist less frequently in the horse genome may affect PCR‐based transgene detection in gene‐doping tests, we propose and demonstrate that PCR amplification and sequencing using primers designed on transgene and promotors and/or polyadenylation signal for gene expression are useful for gene‐doping detection as an additional confirmatory test to prevent false positives. [ABSTRACT FROM AUTHOR]

Peach trees showing witches'‐broom disease symptoms in the northwest of Iran were sampled for phytoplasma detection. PCR assays and Sanger sequence analyses indicated that ʻCandidatus Phytoplasma phoeniciumʼ was associated with peach witchesʼ‐broom disease. Virtual RFLP analyses of the 16S rRNA gene indicated that ʻCa. Phytoplasma phoeniciumʼ strain, which was prevalent in the northwest of Iran belonged to 16SrIX‐C subgroup. For the genomic characterization of Iranian ʻCa. Phytoplasma phoeniciumʼ strain, total DNA extracted from the infected peach trees was subjected to Illumina next‐generation sequencing. The NGS sequencing resulted in 41157647 read pairs of raw data. The raw reads length was 150 bp and the insert size was 350 bp. Assembly and genes clustering finally resulted in 750 contigs with a total size of 228345 bp. The GC content of sequence was 33.2%, N50 value was 259 and L50 value was 256. According to NCBI‐ PGAP annotation the genome of Prunus persica phytoplasma PP2 contained 410 genes including 377 CDSs with protein, 10 rRNAs, nine tRNAs and 14 pseudogenes. [ABSTRACT FROM AUTHOR]

The major histocompatibility complex (MHC) class I region of cattle is both highly polymorphic and, unlike many species, highly variable in gene content between haplotypes. Cattle MHC class I alleles were historically grouped by sequence similarity in the more conserved 3′ end of the coding sequence to form phylogenetic allele groups. This has formed the basis of current cattle MHC class I nomenclature. We presently describe and compare five fully assembled MHC class I haplotypes using the latest cattle and yak genome assemblies. Of the five previously described "pseudogenes" in the cattle MHC class I region, Pseudogene 3 is putatively functional in all haplotypes and Pseudogene 6 and Pseudogene 7 are putatively functional in some haplotypes. This was reinforced by evidence of transcription. Based on full gene sequences as well as 3′ coding sequence, we identified distinct subgroups of BoLA‐3 and BoLA‐6 that represent distinct genetic loci. We further examined allele‐specific expression using transcriptomic data revealing that certain alleles are consistently weakly expressed compared to others. These observations will help to inform further studies into how MHC class I region variability influences T cell and natural killer cell functions in cattle. [ABSTRACT FROM AUTHOR]

Baculoviruses, like their holometabolous host insects, comprise a monophyletic clade and, since the evolutionary origin of this virus family, have coevolved with their hosts. Despite this intimate coevolutionary relationship, baculoviruses have occasionally diverged independently of the host phylogeny by acquiring new infection mechanisms. Given that current baculovirus genomes contain various anti‐apoptotic genes and pseudogenes, the possibility arises that ancient baculoviruses had the genetic means to counter antiviral host apoptosis mechanisms. In this review, we propose an evolutionary scenario for baculovirus diversification in which the development of novel antiviral apoptosis mechanisms in host insects provided the necessary pressure to produce new baculovirus lineages possessing novel anti‐apoptotic genes. This apoptotic arms race and subsequent viral host adaptation likely play a crucial role in the establishment of baculovirus diversity. [ABSTRACT FROM AUTHOR]