Housekeeping protein-coding genes interrogated with tissue and individual variations.

Housekeeping protein-coding genes are stably expressed genes in cells and tissues that are thought to be engaged in fundamental cellular biological functions. They are often utilized as normalization references in molecular biology research and are especially important in integrated bioinformatic investigations. Prior studies have examined human housekeeping protein-coding genes by analyzing various gene expression datasets. The inclusion of different tissue types significantly impacted the discovery of housekeeping genes. In this report, we investigated particularly individual human subject expression differences in protein-coding genes across different tissue types. We used GTEx V8 gene expression datasets obtained from more than 16,000 human normal tissue samples. Furthermore, the Gini index is utilized to investigate the expression variations of protein-coding genes between tissue and individual donor subjects. Housekeeping protein-coding genes found using Gini index profiles may vary depending on the tissue subtypes investigated, particularly given the diverse sample size collections across the GTEx tissue subtypes. We subsequently selected major tissues and identified subsets of housekeeping genes with stable expression levels among human donors within those tissues. In this work, we provide alternative sets of housekeeping protein-coding genes that show more consistent expression patterns in human subjects across major solid organs. Weblink: https://hpsv.ibms.sinica.edu.tw .

Distinctive physical properties of DNA shared by RNA polymerase II gene promoters and 5'-flanking regions of tRNA genes.

Numerous noncoding (nc)RNAs have been identified. Similar to the transcription of protein-coding (mRNA) genes, long noncoding (lnc)RNA genes and most of micro (mi)RNA genes are transcribed by RNA polymerase II (Pol II). In the transcription of mRNA genes, core promoters play an indispensable role; they support the assembly of the preinitiation complex (PIC). However, the structural and/or physical properties of the core promoters of lncRNA and miRNA genes remain largely unexplored, in contrast with those of mRNA genes. Using the core promoters of human genes, we analyzed the repertoire and population ratios of residing core promoter elements (CPEs) and calculated the following five DNA physical properties (DPPs): duplex DNA free energy, base stacking energy, protein-induced deformability, rigidity and stabilizing energy of Z-DNA. Here, we show that their CPE and DPP profiles are similar to those of mRNA gene promoters. Importantly, the core promoters of these three classes of genes have two highly distinctive sites in their DPP profiles around the TSS and position -27. Similar characteristics in DPPs are also found in the 5'-flanking regions of tRNA genes, indicating their common essential roles in transcription initiation over the kingdom of RNA polymerases.

The evolutionary characterization of Gamasida based on mitochondrial genes codon usage pattern.

Mites belonging to the suborder Gamasida are species-rich and habitat-diverse, with a worldwide distribution. To adapt to the environment and obtain better living conditions, all species of the suborder Gamasida have been undergoing constant evolution. The complete mitochondrial genome (mitogenome) is an invaluable molecular marker for studying the origin of species, genetic differentiation between closely related species, and between intraspecific groups. In some species of the suborder Gamasida, mitochondrial tRNA genes are truncated and carried unstable genetic information. This study presents a comparative analysis of codon usage pattern and preference of 13 protein-coding genes of 24 species in 17 genera and 10 families of the suborder Gamasida. Results showed that have an obvious AT preference (0.664-0.829) for codon usage in the suborder Gamasida. Most of the optimal and high-frequency codons also end in A/T. The degree of natural selection varies between the same protein-coding genes of different gamasid mites or among different protein-coding genes within the same gamasid mites. Base and codon usage pattern and preference are very similar between the same species and genus, namely the closer species, the more similar their bases and codons usage patterns and preference are. T bases and C bases were the preference bases for codon usage of 24 species in the suborder Gamasida. Evolution of the suborder Gamasida was dominated by natural selection (64.1%). This study provides the first comprehensive analysis of codon usage in the suborder Gamasida, which will greatly improve our understanding of codon usage patterns and preference, genetics, and evolution of the suborder Gamasida. It will help to evaluate the degree of molecular adaptation in the suborder Gamasida and to further explore evolutionary features of the suborder Gamasida.

Galba: genome annotation with miniprot and AUGUSTUS.

BACKGROUND: The Earth Biogenome Project has rapidly increased the number of available eukaryotic genomes, but most released genomes continue to lack annotation of protein-coding genes. In addition, no transcriptome data is available for some genomes. RESULTS: Various gene annotation tools have been developed but each has its limitations. Here, we introduce GALBA, a fully automated pipeline that utilizes miniprot, a rapid protein-to-genome aligner, in combination with AUGUSTUS to predict genes with high accuracy. Accuracy results indicate that GALBA is particularly strong in the annotation of large vertebrate genomes. We also present use cases in insects, vertebrates, and a land plant. GALBA is fully open source and available as a docker image for easy execution with Singularity in high-performance computing environments. CONCLUSIONS: Our pipeline addresses the critical need for accurate gene annotation in newly sequenced genomes, and we believe that GALBA will greatly facilitate genome annotation for diverse organisms.

Extant thrips diverged in the early tertiary period.

Thysanoptera, commonly known as thrips, are diverse tiny insects whose earliest fossil record dates back to the Triassic period. Since there are few studies on the divergence time of taxa under Thysanoptera, this study used 13 mitochondrial coding protein genes to reconstruct the phylogenetic tree with divergence time of 26 species of this order and show a reliable phylogenetic relationship of thrips species. The time tree of this study shows that most extant thrips diverged in the early Tertiary period, while the fossil records also support that most extinct thrips appeared in this period. This study expands our understanding of the evolution of thrips and provides a feasible way of using multiple mitochondrial genes to establish robust phylogenetic relationships and explore divergence time between species.

Transcriptome Screening of Long Noncoding RNAs and Their Target Protein-Coding Genes Unmasks a Dynamic Portrait of Seed Coat Coloration Associated with Anthocyanins in Tibetan Hulless Barley.

Many plants have the capability to accumulate anthocyanins for coloration, and anthocyanins are advantageous to human health. In the case of hulless barley (Hordeum vulgare L. var. nudum), investigation into the mechanism of anthocyanin formation is limited to the level of protein-coding genes (PCGs). Here, we conducted a comprehensive bioinformatics analysis to identify a total of 9414 long noncoding RNAs (lncRNAs) in the seed coats of purple and white hulless barley along a developmental gradient. Transcriptome-wide profiles of lncRNAs documented several properties, including GC content fluctuation, uneven length, a diverse range of exon numbers, and a wide variety of transcript classifications. We found that certain lncRNAs in hulless barley possess detectable sequence conservation with Hordeum vulgare and other monocots. Furthermore, both differentially expressed lncRNAs (DElncRNAs) and PCGs (DEPCGs) were concentrated in the later seed development stages. On the one hand, DElncRNAs could potentially cis-regulate DEPCGs associated with multiple metabolic pathways, including flavonoid and anthocyanin biosynthesis in the late milk and soft dough stages. On the other hand, there was an opportunity for trans-regulated lncRNAs in the color-forming module to affect seed coat color by upregulating PCGs in the anthocyanin pathway. In addition, the interweaving of hulless barley lncRNAs and diverse TFs may function in seed coat coloration. Notably, we depicted a dynamic portrait of the anthocyanin synthesis pathway containing hulless barley lncRNAs. Therefore, this work provides valuable gene resources and more insights into the molecular mechanisms underlying anthocyanin accumulation in hulless barley from the perspective of lncRNAs, which facilitate the development of molecular design breeding in crops.

Comprehensive analysis of candidate signatures of long non-coding RNA LINC01116 and related protein-coding genes in patients with hepatocellular carcinoma.

BACKGROUND: Hepatocellular carcinoma (HCC) is a long-term malignancy that causes high morbidities and mortalities worldwide. Notably, long non-coding RNAs (LncRNAs) have been identified as candidate targets for malignancy treatments. METHODS: LncRNA LINC01116 and its Pearson-correlated genes (PCGs) were identified and analyzed in HCC patients. The diagnostic and prognostic value of the lncRNA was evaluated using data from The Cancer Genome Atlas (TCGA). Further, we explored the target drugs of LINC01116 for clinical application. Relationships between immune infiltration and PCGs, methylation and PCGs were explored. The diagnostic potentials were then validated by Oncomine cohorts. RESULTS: LINC01116 and the PCG OLFML2B are differentially and highly expressed in tumor tissues (both P ≤ 0.050). We found that LINC01116, TMSB15A, PLAU, OLFML2B, and MRC2 have diagnostic potentials (all AUC ≥ 0.700, all P ≤ 0.050) while LINC01116 and TMSB15A have prognostic significance (both adjusted P ≤ 0.050). LINC01116 was enriched in the vascular endothelial growth factor (VEGF) receptor signaling pathway, mesenchyme morphogenesis, etc. After that, candidate target drugs with potential clinical significance were identified: Thiamine, Cromolyn, Rilmenidine, Chlorhexidine, Sulindac_sulfone, Chloropyrazine, and Meprylcaine. Analysis of immune infiltration revealed that MRC2, OLFML2B, PLAU, and TMSB15A are negatively associated with the purity but positively associated with the specific cell types (all P < 0.050). Analysis of promoter methylation demonstrated that MRC2, OLFML2B, and PLAU have differential and high methylation levels in primary tumors (all P < 0.050). Validation results of the differential expressions and diagnostic potential of OLFML2B (Oncomine) were consistent with those obtained in the TCGA cohort (P < 0.050, AUC > 0.700). CONCLUSIONS: Differentially expressed LINC01116 could be a candidate diagnostic and an independent prognostic signature in HCC. Besides, its target drugs may work for HCC therapy via the VEGF receptor signaling pathway. Differentially expressed OLFML2B could be a diagnostic signature involved in HCC via immune infiltrates.

Tissue heterogeneity of mitochondrial activity, biogenesis and mitochondrial protein gene expression in buffalo.

BACKGROUND: Cellular metabolism is most invariant process, occurring in all living organisms, which involves mitochondrial proteins from both nuclear and mitochondrial genomes. The mitochondrial DNA (mtDNA) copy number, protein-coding genes (mtPCGs) expression, and activity vary between various tissues to fulfill specific energy demands across the tissues. METHODS AND RESULTS: In present study, we investigated the OXPHOS complexes and citrate synthase activity in isolated mitochondria from various tissues of freshly slaughtered buffaloes (n = 3). Further, the evaluation of tissue-specific diversity based on the quantification of mtDNA copy numbers was performed and also comprised an expression study of 13 mtPCGs. We found that the functional activity of individual OXPHOS complex I was significantly higher in the liver compared to muscle and brain. Additionally, OXPHOS complex III and V activities was observed significantly higher levels in liver compared to heart, ovary, and brain. Similarly, CS-specific activity differs between tissues, with the ovary, kidney, and liver having significantly greater. Furthermore, we revealed the mtDNA copy number was strictly tissue-specific, with muscle and brain tissues exhibiting the highest levels. Among 13 PCGs expression analyses, mRNA abundances in all genes were differentially expressed among the different tissue. CONCLUSIONS: Overall, our results indicate the existence of a tissue-specific variation in mitochondrial activity, bioenergetics, and mtPCGs expression among various types of buffalo tissues. This study serves as a critical first stage in gathering vital comparable data about the physiological function of mitochondria in energy metabolism in distinct tissues, laying the groundwork for future mitochondrial based diagnosis and research.

Evolutionary and phylogenetic analyses of 11 Cerasus species based on the complete chloroplast genome.

The subgenus Cerasus, one of the most important groups in the genus Prunus sensu lato, comprises over 100 species; however, the taxonomic classification and phylogenetic relationships of Cerasus remain controversial. Therefore, it is necessary to reconstruct the phylogenetic tree for known Cerasus species. Here, we report the chloroplast (cp) genome sequences of 11 Cerasus species to provide insight into evolution of the plastome. The cp genomes of the 11 Cerasus species (157,571-158,830 bp) displayed a typical quadripartite circular structure. The plastomes contain 115 unique genes, including 80 protein-coding genes, four ribosomal RNAs, and 31 transfer RNAs. Twenty genes were found to be duplicated in inverted repeats as well as at the boundary. The conserved non-coding sequences showed significant divergence compared with the coding regions. We found 12 genes and 14 intergenic regions with higher nucleotide diversity and more polymorphic sites, including matK, rps16, rbcL, rps16-trnQ, petN-psbM, and trnL-trnF. During cp plastome evolution, the codon profile has been strongly biased toward the use of A/T at the third base, and leucine and isoleucine codons appear the most frequently. We identified strong purifying selection on the rpoA, cemA, atpA, and petB genes; whereas ccsA, rps19, matK, rpoC2, ycf2 and ndhI showed a signature of possible positive selection during the course of Cerasus evolution. In addition, we further analyzed the phylogenetic relationships of these species with 57 other congenic related species.Through reconstructing the Cerasus phylogeny tree, we found that true cherry is similar to the flora of China forming a distinct group, from which P. mahaleb was separated as an independent subclade. Microcerasus was genetically closer to Amygdalus, Armeniaca, and Prunus (sensu stricto) than to members of true cherry, whereas P. japonica and P. tomentosa were most closely related to P. triloba and P. pedunculata. However, P. tianshanica formed a clade with P. cerasus, P. fruticosa, P. cerasus × P. canescens 'Gisela 6', and P. avium as a true cherry group. These results provide new insights into the plastome evolution of Cerasus, along with potential molecular markers and candidate DNA barcodes for further phylogenetic and phylogeographic analyses of Cerasus species.

Identification of Long Non-Coding RNAs and Their Target Genes from Mycelium and Primordium in Model Mushroom Schizophyllum commune.

Schizophyllum commune has emerged as the most promising model mushroom to study developmental stages (mycelium, primordium), which are two primary processes of fruit body development. Long non-coding RNA (lncRNA) has been proved to participate in fruit development and sex differentiation in fungi. However, potential lncRNAs have not been identified in S. commune from mycelium to primordium developmental stages. In this study, lncRNA-seq was performed in S. commune and 61.56 Gb clean data were generated from mycelium and primordium developmental stages. Furthermore, 191 lncRNAs had been obtained and a total of 49 lncRNAs were classified as differently expressed lncRNAs. Additionally, 26 up-regulated differently expressed lncRNAs and 23 down-regulated between mycelium and primordia libraries were detected. Further, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that differentially expressed lncRNAs target genes from the MAPK pathway, phosphatidylinositol signal, ubiquitin-mediated proteolysis, autophagy, and cell cycle. This study provides a new resource for further research on the relationship between lncRNA and two developmental stages (mycelium, primordium) in S. commune.

Bioinformatics analysis identified RGS4 as a potential tumor promoter in glioma.

Gliomas is the most common type of intracranial primary malignant tumor and it accounts for ∼80% of primary malignant tumors of the central nervous system. At present, surgical resection with adjuvant radiotherapy and temozolomide adjuvant chemotherapy combined with radiotherapy are the only standard treatments for glioma. However, but overall survival of patients is only 15 months. Glioma is resistant to radiotherapy and chemotherapy, and this malignant behavior leads to a high recurrence rate. Therefore, the use of therapeutics is usually ineffective. As a result, patients with glioma do not significantly benefit from standard treatment. There is therefore an urgent need to develop novel diagnostic approaches and, in particular, more effective treatment strategies. The application of gene expression microarrays provides a feasible and effective way to study gliomas. The present study therefore aimed to identify the key protein-coding genes of glioma using bioinformatics methods and thereby search, for novel biomarkers and therapeutic targets for the treatment of glioma. First, mRNA microarray datasets were selected and obtained from the Gene Expression Omnibus database to identify differentially expressed genes (DEGs) between gliomas and normal tissues. The DEGs were clarified using Gene Ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG), Protein-Protein Interaction (PPI) network and statistical analysis. Subsequently, reverse transcription-quantitative PCR (RT-qPCR)and western blot were used to verify the results of the bioinformatics analysis. A total of 400 DEGs were identified in glioma and they were enriched in several cancer-related GO and KEGG pathways. In the PPI network, it was observed that G-protein signal regulatory protein 4 (RGS4), thymidine phosphorylase, collagen type VI alpha-1, Src homology 2 domain-containing transforming protein1(SHC1) and ring finger protein 135 exhibited a strong protein-protein interaction. Furthermore, . Subsequently, brain damaged tissues and glioma cell lines were selected for RT-qPCR and western blotting analysis. The results demonstrated that RGS4 was highly expressed in glioma cell lines. In conclusion, RGS4 may be a key protein-coding gene in glioma. RGS4 should therefore be studied further to verify its feasibility and effectiveness as a potential glioma biomarker and therapeutic target.

DNA Methylome and LncRNAome Analysis Provide Insights Into Mechanisms of Genome-Dosage Effects in Autotetraploid Cassava.

Whole genome duplication (WGD) increases the dosage of all coding and non-coding genes, yet the molecular implications of genome-dosage effects remain elusive. In this study, we generated integrated maps of the methylomes and lncRNAomes for diploid and artificially generated autotetraploid cassava (Manihot esculenta Crantz). We found that transposable elements (TEs) suppressed adjacent protein coding gene (PCG)-expression levels, while TEs activated the expression of nearby long non-coding RNAs (lncRNAs) in the cassava genome. The hypermethylation of DNA transposons in mCG and mCHH sites may be an effective way to suppress the expression of nearby PCGs in autotetraploid cassava, resulting in similar expression levels for most of PCGs between autotetraploid and diploid cassava. In the autotetraploid, decreased methylation levels of retrotransposons at mCHG and mCHH sites contributed to reduced methylation of Gypsy-neighboring long intergenic non-coding RNAs, potentially preserving diploid-like expression patterns in the major of lncRNAs. Collectively, our study highlighted that WGD-induced DNA methylation variation in DNA transposons and retrotransposons may be as direct adaptive responses to dosage of all coding-genes and lncRNAs, respectively.

Three Complete Mitochondrial Genomes of Orestes guangxiensis, Peruphasma schultei, and Phryganistria guangxiensis (Insecta: Phasmatodea) and Their Phylogeny.

Insects of the order Phasmatodea are mainly distributed in the tropics and subtropics and are best known for their remarkable camouflage as plants. In this study, we sequenced three complete mitochondrial genomes from three different families: Orestes guangxiensis, Peruphasma schultei, and Phryganistria guangxiensis. The lengths of the three mitochondrial genomes were 15,896 bp, 16,869 bp, and 17,005 bp, respectively, and the gene composition and structure of the three stick insects were identical to those of the most recent common ancestor of insects. The phylogenetic relationships among stick insects have been chaotic for a long time. In order to discuss the intra- and inter-ordinal relationship of Phasmatodea, we used the 13 protein-coding genes (PCGs) of 85 species for maximum likelihood (ML) and Bayesian inference (BI) analyses. Results showed that the internal topological structure of Phasmatodea had a few differences in both ML and BI trees and long-branch attraction (LBA) appeared between Embioptera and Zoraptera, which led to a non-monophyletic Phasmatodea. Consequently, after removal of the Embioptera and Zoraptera species, we re-performed ML and BI analyses with the remaining 81 species, which showed identical topology except for the position of Tectarchus ovobessus (Phasmatodea). We recovered the monophyly of Phasmatodea and the sister-group relationship between Phasmatodea and Mantophasmatodea. Our analyses also recovered the monophyly of Heteropterygidae and the paraphyly of Diapheromeridae, Phasmatidae, Lonchodidae, Lonchodinae, and Clitumninae. In this study, Peruphasma schultei (Pseudophasmatidae), Phraortes sp. YW-2014 (Lonchodidae), and species of Diapheromeridae clustered into the clade of Phasmatidae. Within Heteropterygidae, O. guangxiensis was the sister clade to O. mouhotii belonging to Dataminae, and the relationship of (Heteropteryginae + (Dataminae + Obriminae)) was recovered.

Comprehensive evaluation of protein-coding sORFs prediction based on a random sequence strategy.

Background: Small open reading frames (sORFs) with protein-coding ability present unprecedented challenge for genome annotation because of their short sequence and low expression level. In the past decade, only several prediction methods have been proposed for discovery of protein-coding sORFs and lack of objective and uniform negative datasets has become an important obstacle to sORFs prediction. The prediction efficiency of current sORFs prediction methods needs to be further evaluated to provide better research strategies for protein-coding sORFs discovery. Methods: In this work, nine mainstream existing methods for predicting protein-coding potential of ORFs are comprehensively evaluated based on a random sequence strategy. Results: The results show that the current methods perform poorly on different sORFs datasets. For comparison, a sequence based prediction algorithm trained on prokaryotic sORFs is proposed and its better prediction performance indicates that the random sequence strategy can provide feasible ideas for protein-coding sORFs predictions. Conclusions: As a kind of important functional genomic element, discovery of protein-coding sORFs has shed light on the dark proteomes. This evaluation work indicates that there is an urgent need for developing specialized prediction tools for protein-coding sORFs in both eukaryotes and prokaryotes. It is expected that the present work may provide novel ideas for future sORFs researches.

A Methodology to Study Pseudogenized lincRNAs.

Long intergenic noncoding RNAs (lincRNAs) are known to be tissue specifically expressed and able to regulate functional protein-coding genes: some can even act as competing endogenous RNAs (ceRNAs), because microRNAs can bind to them instead of the corresponding mRNA binding sites. Some lincRNAs contain remnants of protein-coding sequences and it has been hypothesized that they might arise after a pseudogenization processes. However, a major limitation in the study of such phenomenon is the lack of proper computational tools designed to align/analyze protein-coding sequences and noncoding sequences. To overcome this limitation, we published a method that finds the remnants of protein-coding sequences within the sequence of lincRNAs, as well as the corresponding sequences in parental proteins. This method, together with the visualization platform for tracing frameshifts and single point mutations within this type of sequences, are described here.

Characterization of the Complete Mitochondrial Genome of Basidiomycete Yeast Hannaella oryzae: Intron Evolution, Gene Rearrangement, and Its Phylogeny.

In this study, the mitogenome of Hannaella oryzae was sequenced by next-generation sequencing (NGS) and successfully assembled. The H. oryzae mitogenome comprised circular DNA molecules with a total size of 26,444 bp. We found that the mitogenome of H. oryzae partially deleted the tRNA gene transferring cysteine. Comparative mitogenomic analyses showed that intronic regions were the main factors contributing to the size variations of mitogenomes in Tremellales. Introns of the cox1 gene in Tremellales species were found to have undergone intron loss/gain events, and introns of the H. oryzae cox1 gene may have different origins. Gene arrangement analysis revealed that H. oryzae contained a unique gene order different from other Tremellales species. Phylogenetic analysis based on a combined mitochondrial gene set resulted in identical and well-supported topologies, wherein H. oryzae was closely related to Tremella fuciformis. This study represents the first report of mitogenome for the Hannaella genus, which will allow further study of the population genetics, taxonomy, and evolutionary biology of this important phylloplane yeast and other related species.

Genes on the circular code alphabet.

The X motifs, motifs from the circular code X, are enriched in the (protein coding) genes of bacteria, archaea, eukaryotes, plasmids and viruses, moreover, in the minimal gene set belonging to the three domains of life, as well as in tRNA and rRNA sequences. They allow to retrieve, maintain and synchronize the reading frame in genes, and contribute to the regulation of gene expression. These results lead here to a theoretical study of genes based on the circular code alphabet. A new occurrence relation of the circular code X under the hypothesis of an equiprobable (balanced) strand pairing is given. Surprisingly, a statistical analysis of a large set of bacterial genes retrieves this relation on the circular code alphabet, but not on the DNA alphabet. Furthermore, the circular code X has the strongest balanced circular code pairing among 216 maximal C3 self-complementary trinucleotide circular codes, a new property of this circular code X. As an application of this theory, different tRNAs studied on the circular code alphabet reveal an unexpected stem structure. Thus, the circular code X would have constructed a coding stem in tRNAs as an outline of the future gene structure and the future DNA double helix.

Genome-wide transcriptome study using deep RNA sequencing for myocardial infarction and coronary artery calcification.

BACKGROUND: Coronary artery calcification (CAC) is a noninvasive measure of coronary atherosclerosis, the proximal pathophysiology underlying most cases of myocardial infarction (MI). We sought to identify expression signatures of early MI and subclinical atherosclerosis in the Framingham Heart Study (FHS). In this study, we conducted paired-end RNA sequencing on whole blood collected from 198 FHS participants (55 with a history of early MI, 72 with high CAC without prior MI, and 71 controls free of elevated CAC levels or history of MI). We applied DESeq2 to identify coding-genes and long intergenic noncoding RNAs (lincRNAs) differentially expressed in MI and high CAC, respectively, compared with the control. RESULTS: On average, 150 million paired-end reads were obtained for each sample. At the false discovery rate (FDR) < 0.1, we found 68 coding genes and 2 lincRNAs that were differentially expressed in early MI versus controls. Among them, 60 coding genes were detectable and thus tested in an independent RNA-Seq data of 807 individuals from the Rotterdam Study, and 8 genes were supported by p value and direction of the effect. Immune response, lipid metabolic process, and interferon regulatory factor were enriched in these 68 genes. By contrast, only 3 coding genes and 1 lincRNA were differentially expressed in high CAC versus controls. APOD, encoding a component of high-density lipoprotein, was significantly downregulated in both early MI (FDR = 0.007) and high CAC (FDR = 0.01) compared with controls. CONCLUSIONS: We identified transcriptomic signatures of early MI that include differentially expressed protein-coding genes and lincRNAs, suggesting important roles for protein-coding genes and lincRNAs in the pathogenesis of MI.

The first eleven mitochondrial genomes from the ectomycorrhizal fungal genus (Boletus) reveal intron loss and gene rearrangement.

In the present study, eleven novel complete mitogenomes of Boletus were assembled and compared. The eleven complete mitogenomes were all composed of circular DNA molecules, with sizes ranging from 32,883 bp to 48,298 bp. The mitochondrial gene arrangement of Boletus varied greatly from other Boletales mitogenomes, and gene position reversal were observed frequently in the evolution of Boletus. Across the 15 core protein-coding genes (PCGs) tested, atp9 had the least and rps3 had the largest genetic distances among the eleven Boletus species, indicating varied evolution rates of core PCGs. In addition, the Ka/Ks value for nad3 gene was >1, suggesting that this gene was subject to possible positive selection pressure. Comparative mitogenomic analysis indicated that the intronic region was significantly correlated with the size of mitogenomes in Boletales. Two large-scale intron loss events were detected in the evolution of Boletus. Phylogenetic analyses based on a combined mitochondrial gene dataset yielded a well-supported (BPP ≥ 0.99; BS =100) phylogenetic tree for 72 Agaricomycetes, and the Boletus species had a close relationship with Paxillus. This study served as the first report on complete mitogenomes in Boletus, which will further promote investigations of the genetics, evolution and phylogeny of the Boletus genus.