KOBAS-i: intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis

Study on the Multitarget Mechanism and Active Compounds of Essential Oil from Artemisia argyi Treating Pressure Injuries Based on Network Pharmacology

In order to comprehensively explore multitarget mechanism and key active compounds of Artemisia argyi essential oil (AAEO) in the treatment of pressure injuries (PIs), we analyzed the biological functions and pathways involved in the intersection targets of AAEO and PIs based on network pharmacology, and the affinity of AAEO active compounds and core targets was verified by molecular docking finally. In our study, we first screened 54 effective components according to the relative content and biological activity. In total, 103 targets related to active compounds of AAEO and 2760 targets associated with PIs were obtained, respectively, and 50 key targets were overlapped by Venny 2.1.0. The construction of key targets-compounds network was achieved by the STRING database and Cytoscape 3.7.2 software. GO analysis from Matespace shows that GO results are mainly enriched in biological processes, including adrenergic receptor activity, neurotransmitter clearance, and neurotransmitter metabolic process. KEGG analysis by the David and Kobas website shows that the key targets can achieve the treatment on PIs through a pathway in cancer, PI3K-Akt signaling pathway, human immunodeficiency virus 1 infection, MAPK signaling pathway, Wnt signaling pathway, etc. In addition, molecular docking results from the CB-Dock server indicated that active compounds of AAEO had good activity docking with the first 10 key targets. In conclusion, the potential targets and regulatory molecular mechanisms of AAEO in the treatment of PIs were analyzed by network pharmacology and molecular docking. AAEO can cure PIs through the synergistic effect of multicomponent, multitarget, and multipathway, providing a theoretical basis and new direction for further study.

Construction and Comprehensive Analysis of a circRNA-miRNA-mRNA Regulatory Network to Reveal the Pathogenesis of Hepatocellular Carcinoma

Background: Circular RNAs (circRNAs) have been demonstrated to be closely related to the carcinogenesis of human cancer in recent years. However, the molecular mechanism of circRNAs in the pathogenesis of hepatocellular carcinoma (HCC) has not been fully elucidated. We aimed to identify critical circRNAs and explore their potential regulatory network in HCC.

Methods: The robust rank aggregation (RRA) algorithm and weighted gene co-expression network analysis (WGCNA) were conducted to unearth the differentially expressed circRNAs (DEcircRNAs) in HCC. The expression levels of DEcircRNAs were validated by quantitative real-time polymerase chain reaction (qRT-PCR). A circRNA-miRNA-mRNA regulatory network was constructed by computational biology, and protein-protein interaction (PPI) network, functional enrichment analysis, survival analysis, and infiltrating immune cells analysis were performed to uncover the potential regulatory mechanisms of the network.

Results: A total of 22 DEcircRNAs were screened out from four microarray datasets (GSE94508, GSE97332, GSE155949, and GSE164803) utilizing the RRA algorithm. Meanwhile, an HCC-related module containing 404 circRNAs was identified by WGCNA analysis. After intersection, only four circRNAs were recognized in both algorithms. Following qRT-PCR validation, three circRNAs (hsa_circRNA_091581, hsa_circRNA_066568, and hsa_circRNA_105031) were chosen for further analysis. As a result, a circRNA-miRNA-mRNA network containing three circRNAs, 17 miRNAs, and 222 mRNAs was established. Seven core genes (ESR1, BUB1, PRC1, LOX, CCT5, YWHAZ, and DDX39B) were determined from the PPI network of 222 mRNAs, and a circRNA-miRNA-hubgene network was also constructed. Functional enrichment analysis suggested that these seven hub genes were closely correlated with several cancer related pathways. Survival analysis revealed that the expression levels of the seven core genes were significantly associated with the prognosis of HCC patients. In addition, we also found that these seven hub genes were remarkably related to the infiltrating levels of immune cells.

Conclusion: Our research identified three pivotal HCC-related circRNAs and provided novel insights into the underlying mechanisms of the circRNA-miRNA-mRNA regulatory network in HCC.

The proteomic landscape of ovarian cancer cells in response to melatonin

Ovarian cancer (OC) is the most lethal gynecological malignancy with a highly negative prognosis. Melatonin is an indoleamine secreted by the pineal gland during darkness and has shown antitumor activity in both in vitro and in vivo experiments. Herein, we investigated the influence of melatonin on the proteome of human ovarian carcinoma cells (SKOV-3 cell line) using the Ultimate 3000 LC Liquid NanoChromatography equipment coupled to a Q-Exactive mass spectrometry. After 48 h of treatment, melatonin induced a significant cytotoxicity especially with the highest melatonin concentration. The proteomic profile revealed 639 proteins in the control group, and 98, 110, and 128 proteins were altered by melatonin at the doses of 0.8, 1.6, and 2.4 mM, respectively. Proteins associated with the immune system and tricarboxylic acid cycle were increased in the three melatonin-exposed groups of cells. Specifically, the dose of 2.4 mM led to a reduction in molecules associated with protein synthesis, especially those of the ribosomal protein family. We also identified 28 potential genes shared between normal ovarian tissue and OC in all experimental groups, and melatonin was predicted to alter genes encoding ribosomal proteins. Notably, the set of proteins changed by melatonin was linked to a better prognosis for OC patients. We conclude that melatonin significantly alters the proteome of SKOV-3 cells by changing proteins involved with the immune response and mitochondrial metabolism. The concentration of 2.4 mM of melatonin promoted the largest number of protein changes. The evidence suggests that melatonin may be an effective therapeutic strategy against OC.

Identification and Validation of Quantitative Trait Loci Mapping for Spike-Layer Uniformity in Wheat

Spike-layer uniformity (SLU), the consistency of the spike distribution in the vertical space, is an important trait. It directly affects the yield potential and appearance. Revealing the genetic basis of SLU will provide new insights into wheat improvement. To map the SLU-related quantitative trait loci (QTL), 300 recombinant inbred lines (RILs) that were derived from a cross between H461 and Chinese Spring were used in this study. The RILs and parents were tested in fields from two continuous years from two different pilots. Phenotypic analysis showed that H461 was more consistent in the vertical spatial distribution of the spike layer than in Chinese Spring. Based on inclusive composite interval mapping, four QTL were identified for SLU. There were two major QTL on chromosomes 2BL and 2DL and two minor QTL on chromosomes 1BS and 2BL that were identified. The additive effects of QSlu.sicau-1B, Qslu.sicau-2B-2, and QSlu.sicau-2D were all from the parent, H461. The major QTL, QSlu.sicau-2B-2 and QSlu.sicau-2D, were detected in each of the conducted trials. Based on the best linear unbiased prediction values, the two loci explained 23.97% and 15.98% of the phenotypic variation, respectively. Compared with previous studies, the two major loci were potentially novel and the two minor loci were overlapped. Based on the kompetitive allele-specific PCR (KASP) marker, the genetic effects for QSlu.sicau-2B-2 were validated in an additional RIL population. The genetic effects ranged from 26.65% to 32.56%, with an average value of 30.40%. In addition, QSlu.sicau-2B-2 showed a significant (p < 0.01) and positive influence on the spike length, spikelet number, and thousand kernel weight. The identified QTL and the developed KASP marker will be helpful for fine-mapping these loci, finally contributing to wheat breeding programs in a marker-assisted selection way.

Integrated Analysis of Basic Helix Loop Helix Transcription Factor Family and Targeted Terpenoids Reveals Candidate AarbHLH Genes Involved in Terpenoid Biosynthesis in Artemisia argyi

Artemisia argyi is a valuable traditional medicinal plant in Asia. The essential oil from its leaves is rich in terpenoids and has been used to enhance health and well-being. In China, the market scale of industries related to A. argyi has attained tens of billions of Chinese Yuan. The basic helix-loop-helix (bHLH) family is one of the largest transcription factors families in plants that plays crucial roles in diverse biological processes and is an essential regulatory component of terpenoid biosynthesis. However, the bHLH TFs and their regulatory roles in A. argyi remain unknown. Here, 53 AarbHLH genes were identified from the transcriptome of A. argyi and were classified into 15 subfamilies based on the classification of bHLH proteins in Arabidopsis thaliana. The MEME analysis showed that the conserved motif 1 and motif 2 constituted the most conserved bHLH domain and distributed in most AarbHLH proteins. Additionally, integrated analysis of the expression profiles of AarbHLH genes and the contents of targeted terpenoids in different tissues group and JA-treated group were performed. Eleven up-regulated AarbHLHs and one down-regulated AarbHLH were screened as candidate genes that may participate in the regulation of terpenoid biosynthesis (TPS-AarbHLHs). Correlation analysis between gene expression and terpenoid contents indicated that the gene expression of these 12 TPS-AarbHLHs was significantly correlated with the content changes of 1,8-cineole or β-caryophyllene. Protein-protein interaction networks further illustrated that these TPS-AarbHLHs might be involved in terpenoid biosynthesis in A. argyi. This finding provides a basis to further investigate the regulation mechanism of AarbHLH genes in terpenoid biosynthesis, and will be helpful to improve the quality of A. argyi.

Genome analysis of Klebsiella oxytoca complex for antimicrobial resistance and virulence genes

Klebsiella oxytoca complex comprises nine closely-related species causing human infections. We curated genomes labeled Klebsiella (n=14,256) in GenBank and identified 588 belonging to the complex, which were examined for precise species, sequence types, K- and O-antigen types, virulence and antimicrobial resistance genes. The complex and Klebsiella pneumoniae share many K- and O-antigen types. Of the complex, K. oxytoca and Klebsiella michiganensis appear to carry more virulence genes and be more commonly associated with human infections.

Frequency and Spectrum of Mutations Induced by Gamma Rays Revealed by Phenotype Screening and Whole-Genome Re-Sequencing in Arabidopsis thaliana

Genetic variations are an important source of germplasm diversity, as it provides an allele resource that contributes to the development of new traits for plant breeding. Gamma rays have been widely used as a physical agent for mutation creation in plants, and their mutagenic effect has attracted extensive attention. However, few studies are available on the comprehensive mutation profile at both the large-scale phenotype mutation screening and whole-genome mutation scanning. In this study, biological effects on M₁ generation, large-scale phenotype screening in M₂ generation, as well as whole-genome re-sequencing of seven M₃ phenotype-visible lines were carried out to comprehensively evaluate the mutagenic effects of gamma rays on Arabidopsis thaliana. A total of 417 plants with visible mutated phenotypes were isolated from 20,502 M₂ plants, and the phenotypic mutation frequency of gamma rays was 2.03% in Arabidopsis thaliana. On average, there were 21.57 single-base substitutions (SBSs) and 11.57 small insertions and deletions (InDels) in each line. Single-base InDels accounts for 66.7% of the small InDels. The genomic mutation frequency was 2.78 × 10⁻¹⁰/bp/Gy. The ratio of transition/transversion was 1.60, and 64.28% of the C > T events exhibited the pyrimidine dinucleotide sequence; 69.14% of the small InDels were located in the sequence with 1 to 4 bp terminal microhomology that was used for DNA end rejoining, while SBSs were less dependent on terminal microhomology. Nine genes, on average, were predicted to suffer from functional alteration in each re-sequenced line. This indicated that a suitable mutation gene density was an advantage of gamma rays when trying to improve elite materials for one certain or a few traits. These results will aid the full understanding of the mutagenic effects and mechanisms of gamma rays and provide a basis for suitable mutagen selection and parameter design, which can further facilitate the development of more controlled mutagenesis methods for plant mutation breeding.

TcMYC2 regulates Pyrethrin biosynthesis in Tanacetum cinerariifolium

Pyrethrins constitute a class of terpene derivatives with high insecticidal activity and are mainly synthesized in the capitula of the horticulturally important plant, Tanacetum cinerariifolium. Treatment of T. cinerariifolium with methyl jasmonate (MeJA) in the field induces pyrethrin biosynthesis, but the mechanism linking MeJA with pyrethrin biosynthesis remains unclear. In this study, we explored the transcription factors involved in regulating MeJA-induced pyrethrin biosynthesis. A single spray application of MeJA to T. cinerariifolium leaves rapidly upregulated the expression of most known pyrethrin biosynthesis genes and subsequently increased the total pyrethrin content in the leaf. A continuous 2-week MeJA treatment resulted in enhanced pyrethrin content and increased trichome density. TcMYC2, a key gene in jasmonate signaling, was screened at the transcriptome after MeJA treatment. TcMYC2 positively regulated expression of the pyrethrin biosynthesis genes TcCHS, TcAOC, and TcGLIP by directly binding to E-box/G-box motifs in the promoters. The stable overexpression of TcMYC2 in T. cinerariifolium hairy roots significantly increased the expression of TcAOC and TcGLIP. Further transient overexpression and viral-induced gene-silencing experiments demonstrated that TcMYC2 positively promoted pyrethrin biosynthesis. Collectively, the results reveal a novel molecular mechanism for MeJA-induced pyrethrin biosynthesis in T. cinerariifolium involving TcMYC2.

CG and CHG Methylation Contribute to the Transcriptional Control of OsPRR37-Output Genes in Rice

Plant circadian clock coordinates endogenous transcriptional rhythms with diurnal changes of environmental cues. OsPRR37, a negative component in the rice circadian clock, reportedly regulates transcriptome rhythms, and agronomically important traits. However, the underlying regulatory mechanisms of OsPRR37-output genes remain largely unknown. In this study, whole genome bisulfite sequencing and high-throughput RNA sequencing were applied to verify the role of DNA methylation in the transcriptional control of OsPRR37-output genes. We found that the overexpression of OsPRR37 suppressed rice growth and altered cytosine methylations in CG and CHG sequence contexts in but not the CHH context (H represents A, T, or C). In total, 35 overlapping genes were identified, and 25 of them showed negative correlation between the methylation level and gene expression. The promoter of the hexokinase gene OsHXK1 was hypomethylated at both CG and CHG sites, and the expression of OsHXK1 was significantly increased. Meanwhile, the leaf starch content was consistently lower in OsPRR37 overexpression lines than in the recipient parent Guangluai 4. Further analysis with published data of time-course transcriptomes revealed that most overlapping genes showed peak expression phases from dusk to dawn. The genes involved in DNA methylation, methylation maintenance, and DNA demethylation were found to be actively expressed around dusk. A DNA glycosylase, namely ROS1A/DNG702, was probably the upstream candidate that demethylated the promoter of OsHXK1. Taken together, our results revealed that CG and CHG methylation contribute to the transcriptional regulation of OsPRR37-output genes, and hypomethylation of OsHXK1 leads to decreased starch content and reduced plant growth in rice.

Research on Pathogenic Hippocampal Voxel Detection in Alzheimer's Disease Using Clustering Genetic Random Forest

Alzheimer's disease (AD) is an age-related neurological disease, which is closely associated with hippocampus, and subdividing the hippocampus into voxels can capture subtle signals that are easily missed by region of interest (ROI) methods. Therefore, studying interpretable associations between voxels can better understand the effect of voxel set on the hippocampus and AD. In this study, by analyzing the hippocampal voxel data, we propose a novel method based on clustering genetic random forest to identify the important voxels. Specifically, we divide the left and right hippocampus into voxels to constitute the initial feature set. Moreover, the random forest is constructed using the randomly selected samples and features. The genetic evolution is used to amplify the difference in decision trees and the clustering evolution is applied to generate offspring in genetic evolution. The important voxels are the features that reach the peak classification. The results demonstrate that our method has good classification and stability. Particularly, through biological analysis of the obtained voxel set, we find that they play an important role in AD by affecting the function of the hippocampus. These discoveries demonstrate the contribution of the voxel set to AD.

Morroniside ameliorates inflammatory skeletal muscle atrophy via inhibiting canonical and non-canonical NF-κB and regulating protein synthesis/degradation

No drug options exist for skeletal muscle atrophy in clinical, which poses a huge socio-economic burden, making development on drug interventions a general wellbeing need. Patients with a variety of pathologic conditions associated with skeletal muscle atrophy have systemically elevated inflammatory factors. Morroniside, derived from medicinal herb Cornus officinalis, possesses anti-inflammatory effect. However, whether and how morroniside combat muscle atrophy remain unknown. Here, we identified crucial genetic associations between TNFα/NF-κB pathway and grip strength based on population using 377,807 European participants from the United Kingdom Biobank dataset. Denervation increased TNFα in atrophying skeletal muscles, which inhibited myotube formation in vitro. Notably, morroniside treatment rescued TNFα-induced myotube atrophy in vitro and impeded skeletal muscle atrophy in vivo, resulting in increased body/muscles weights, No. of satellite cells, size of type IIA, IIX and IIB myofibers, and percentage of type IIA myofibers in denervated mice. Mechanistically, in vitro and/or in vivo studies demonstrated that morroniside could not only inhibit canonical and non-canonical NF-κB, inflammatory mediators (IL6, IL-1b, CRP, NIRP3, PTGS2, TNFα), but also down-regulate protein degradation signals (Follistatin, Myostatin, ALK4/5/7, Smad7/3), ubiquitin-proteasome molecules (FoxO3, Atrogin-1, MuRF1), autophagy-lysosomal molecules (Bnip3, LC3A, and LC3B), while promoting protein synthesis signals (IGF-1/IGF-1R/IRS-1/PI3K/Akt, and BMP14/BMPR2/ALK2/3/Smad5/9). Moreover, morroniside had no obvious liver and kidney toxicity. This human genetic, cells and mice pathological evidence indicates that morroniside is an efficacious and safe inflammatory muscle atrophy treatment and suggests its translational potential on muscle wasting.

Genetic Basis of Sexual Maturation Heterosis: Insights From Ovary lncRNA and mRNA Repertoire in Chicken

Sexual maturation is fundamental to the reproduction and production performance, heterosis of which has been widely used in animal crossbreeding. However, the underlying mechanism have long remained elusive, despite its profound biological and agricultural significance. In the current study, the reciprocal crossing between White Leghorns and Beijing You chickens were performed to measure the sexual maturation heterosis, and the ovary lncRNAs and mRNAs of purebreds and crossbreeds were profiled to illustrate molecular mechanism of heterosis. Heterosis larger than 20% was found for pubic space and oviduct length, whereas age at first egg showed negative heterosis in both crossbreeds. We identified 1170 known lncRNAs and 1994 putative lncRNAs in chicken ovary using a stringent pipeline. Gene expression pattern showed that nonadditivity was predominant, and the proportion of nonadditive lncRNAs and genes was similar between two crossbreeds, ranging from 44.24% to 49.15%. A total of 200 lncRNAs and 682 genes were shared by two crossbreeds, respectively. GO and KEGG analysis showed that the common genes were significantly enriched in the cell cycle, animal organ development, gonad development, ECM-receptor interaction, calcium signaling pathway and GnRH signaling pathway. Weighted gene co-expression network analysis (WGCNA) identified that 7 out of 20 co-expressed lncRNA-mRNA modules significantly correlated with oviduct length and pubic space. Interestingly, genes harbored in seven modules were also enriched in the similar biological process and pathways, in which nonadditive lncRNAs, such as MSTRG.17017.1 and MSTRG.6475.20, were strongly associated with nonadditive genes, such as CACNA1C and TGFB1 to affect gonad development and GnRH signaling pathway, respectively. Moreover, the results of real-time quantitative PCR (RT-qPCR) correlated well with the transcriptome data. Integrated with positive heterosis of serum GnRH and melatonin content detected in crossbreeds, we speculated that nonadditive genes involved in the GnRH signaling pathway elevated the gonad development, leading to the sexual maturation heterosis. We characterized a systematic landscape of ovary lncRNAs and mRNAs related to sexual maturation heterosis in chicken. The quantitative exploration of hybrid transcriptome changes lays foundation for genetic improvement of sexual maturation traits and provides insights into endocrine control of sexual maturation.

Genome-wide association mapping of Fusarium crown rot resistance in Aegilops tauschii

Fusarium crown rot (FCR), caused by various Fusarium species, is a primary fungal disease in most wheat-growing regions worldwide. A. tauschii, the diploid wild progenitor of the D-genome of common wheat, is a reservoir of genetic diversity for improving bread wheat biotic and abiotic resistance/tolerance. A worldwide collection of 286 A. tauschii accessions was used to evaluate FCR resistance. Population structure analysis revealed that 115 belonged to the A. tauschii ssp. strangulata subspecies, and 171 belonged to the A. tauschii ssp. tauschii subspecies. Five accessions with disease index values lower than 20 showed moderate resistance to FCR. These five originated from Afghanistan, China, Iran, Uzbekistan, and Turkey, all belonging to the tauschii subspecies. Genome-wide association mapping using 6,739 single nucleotide polymorphisms (SNPs) revealed that two SNPs on chromosome 2D and four SNPs on chromosome 7D were significantly associated with FCR resistance. Almost all FCR resistance alleles were presented in accessions from the tauschii subspecies, and only 4, 11, and 19 resistance alleles were presented in accessions from the strangulata subspecies. Combining phenotypic correlation analysis and genome-wide association mapping confirmed that FCR resistance loci were independent of flowering time, heading date, and plant height in this association panel. Six genes encoding disease resistance-related proteins were selected as candidates for further validation. The identified resistant A. tauschii accessions will provide robust resistance gene sources for breeding FCR-resistant cultivars. The associated loci/genes will accelerate and improve FCR in breeding programs by deploying marker-assisted selection.

Transcriptome-wide N 6-Methyladenosine Methylome Profiling Reveals m6A Regulation of Skeletal Myoblast Differentiation in Cattle (Bos taurus)

N ⁶ -methyladenosine (m⁶A) is the most prevalent methylation modification of eukaryotic mRNA, and it plays an important role in regulating gene expression. Previous studies have found that m6A methylation plays a role in mammalian skeletal muscle development. However, the effect of m⁶A on bovine skeletal myogenesis are still unclear. Here, we selected proliferating myoblasts (GM) and differentiated myotubes (on the 4th day of differentiation, DM) for m⁶A-seq and RNA-seq to explore the m⁶A methylation modification pattern during bovine skeletal myogenesis. m⁶A-seq analysis revealed that m⁶A methylation was an abundant modification of the mRNA in bovine myoblasts and myotubes. We scanned 5,691-8,094 m⁶A-modified transcripts, including 1,437 differentially methylated genes (DMGs). GO and KEGG analyses revealed that DMGs were primarily involved in transcriptional regulation and RNA metabolism, as well as insulin resistance and metabolic pathways related to muscle development. The combined analysis further identified 268 genes that had significant changes at both m⁶A and mRNA levels, suggesting that m⁶A modification may regulate myoblast differentiation by mediating the expression of these genes. Furthermore, we experimentally confirmed four genes related to myogenesis, including MYOZ2, TWIST1, KLF5 and MYOD1, with differential changes in both m⁶A and mRNA levels during bovine myoblast differentiation, indicating that they can be potential candidate targets for m⁶A regulation of skeletal myogenesis. Our results may provide new insight into molecular genetics and breeding of beef cattle, and provide a reference for investigating the mechanism of m⁶A regulating skeletal muscle development.

Signatures of selection reveal candidate genes involved in production traits in Chinese crossbred buffaloes

Identification of selection signature is important for a better understanding of genetic mechanisms that affect phenotypic differentiation in livestock. However, the genome-wide selection responses have not been investigated for the production traits of Chinese crossbred buffaloes. In this study, an SNP data set of 133 buffaloes (Chinese crossbred buffalo, n = 45; Chinese local swamp buffalo, n = 88) was collected from the Dryad Digital Repository database (https://datadryad.org/stash/). Population genetics analysis showed that these buffaloes were divided into the following 2 groups: crossbred buffalo and swamp buffalo. The crossbred group had higher genetic diversity than the swamp group. Using 3 complementary statistical methods (integrated haplotype score, cross population extended haplotype homozygosity, and composite likelihood ratio), a total of 31 candidate selection regions were identified in the Chinese crossbred population. Here, within these candidate regions, 25 genes were under the putative selection. Among them, several candidate genes were reported to be associated with production traits. In addition, we identified 13 selection regions that overlapped with bovine QTLs that were mainly involved in milk production and composition traits. These results can provide useful insights regarding the selection response for production traits of Chinese crossbred buffalo, as identified candidate genes influence production performance.

Different Regulatory Modes of Synechocystis sp. PCC 6803 in Response to Photosynthesis Inhibitory Conditions

Cyanobacteria are promising industrial platforms owing to their ability to produce diverse natural secondary metabolites and nonnative value-added biochemicals from CO₂ and light. To fully utilize their industrial potency, it is critical to understand their photosynthetic efficiency under various environmental conditions. In this study, we elucidated the inhibitory mechanisms of photosynthesis under high-light and low-temperature stress conditions in the model cyanobacterium Synechocystis sp. PCC 6803. Under each stress condition, the transcript abundance and translation efficiency were measured using transcriptome sequencing (RNA-seq) and ribosome profiling, and the genome-wide transcription unit architecture was constructed by data integration of transcription start sites and transcript 3′-end positions obtained from differential RNA-seq and sequencing of 3′-ends (Term-seq), respectively. Our results suggested that the mode of photosynthesis inhibition differed between the two stress conditions; high light stress induced photodamage responses, while low temperature stress impaired the translation efficiency of photosynthesis-associated genes. In particular, poor translation of photosystem I resulted from ribosome stalling at the untranslated regions, affecting the overall photosynthetic yield under low temperature stress. Our comprehensive multiomics analysis with transcription unit architecture provides foundational information on photosynthesis for future industrial strain development.

IMPORTANCE Cyanobacteria are a compelling biochemical production platform for their ability to propagate using light and atmospheric CO₂ via photosynthesis. However, the engineering of strains is hampered by limited understanding of photosynthesis under diverse environmental conditions such as high-light and low-temperature stresses. Herein, we decipher the transcriptomic and translatomic responses of the photosynthetic efficiency to stress conditions using the integrative analysis of multiomic data generated by RNA-seq and ribosome profiling, respectively. Through the generated massive data, along with the guide of the genome-wide transcription unit architecture constructed by transcription start sites and transcript 3′-end positions, we identified the factors affecting photosynthesis at transcription, posttranscription, and translation levels. Importantly, the high-light stress induces photodamage responses, and the low-temperature stress cripples the translation efficiency of photosynthesis-associated genes. The resulting insights provide pivotal information for future cyanobacterial cell factories powered by the engineering toward robust photosynthesis ability.

Transcriptome-Wide Analyses Identify Dominant as the Predominantly Non-Conservative Alternative Splicing Inheritance Patterns in F1 Chickens

Transcriptome analysis has been used to investigate many economically traits in chickens; however, alternative splicing still lacks a systematic method of study that is able to promote proteome diversity, and fine-tune expression dynamics. Hybridization has been widely utilized in chicken breeding due to the resulting heterosis, but the dynamic changes in alternative splicing during this process are significant yet unclear. In this study, we performed a reciprocal crossing experiment involving the White Leghorn and Cornish Game chicken breeds which exhibit major differences in body size and reproductive traits, and conducted RNA sequencing of the brain, muscle, and liver tissues to identify the inheritance patterns. A total of 40 515 and 42 612 events were respectively detected in the brain and muscle tissues, with 39 843 observed in the liver; 2807, 4242, and 4538 events significantly different between two breeds were identified in the brain, muscle, and liver tissues, respectively. The hierarchical cluster of tissues from different tissues from all crosses, based on the alternative splicing profiles, suggests high tissue and strain specificity. Furthermore, a comparison between parental strains and hybrid crosses indicated that over one third of alternative splicing genes showed conserved patterns in all three tissues, while the second prevalent pattern was non-additive, which included both dominant and transgressive patterns; this meant that the dominant pattern plays a more important role than suppression. Our study provides an overview of the inheritance patterns of alternative splicing in layer and broiler chickens, to better understand post-transcriptional regulation during hybridization.

Two novel RNA-binding proteins identification through computational prediction and experimental validation

Since RBPs play important roles in the cell, it's particularly important to find new RBPs. We performed iRIP-seq and CLIP-seq to verify two proteins, CLIP1 and DMD, predicted by RBPPred whether are RBPs or not. The experimental results confirm that these two proteins have RNA-binding activity. We identified significantly enriched binding motifs UGGGGAGG, CUUCCG and CCCGU for CLIP1 (iRIP-seq), DMD (iRIP-seq) and DMD (CLIP-seq), respectively. The computational KEGG and GO analysis show that the CLIP1 and DMD share some biological processes and functions. Besides, we found that the SNPs between DMD and its RNA partners may be associated with Becker muscular dystrophy, Duchenne muscular dystrophy, Dilated cardiomyopathy 3B and Cardiovascular phenotype. Among the thirteen cancers data, CLIP1 and another 300 oncogenes always co-occur, and 123 of these 300 genes interact with CLIP1. These cancers may be associated with the mutations occurred in both CLIP1 and the genes it interacts with.

Chromosome-Level Genome Assembly Provides New Insights into Genome Evolution and Tuberous Root Formation of Potentilla anserina

Potentilla anserina is a perennial stoloniferous plant with edible tuberous roots in Rosaceae, served as important food and medicine sources for Tibetans in the Qinghai-Tibetan Plateau (QTP), China, over thousands of years. However, a lack of genome information hindered the genetic study. Here, we presented a chromosome-level genome assembly using single-molecule long-read sequencing, and the Hi-C technique. The assembled genome was 454.28 Mb, containing 14 chromosomes, with contig N50 of 2.14 Mb. A total of 46,495 protein-coding genes, 169.74 Mb repeat regions, and 31.76 Kb non-coding RNA were predicted. P. anserina diverged from Potentilla micrantha ∼28.52 million years ago (Mya). Furthermore, P. anserina underwent a recent tetraploidization ∼6.4 Mya. The species-specific genes were enriched in Starch and sucrose metabolism and Galactose metabolism pathways. We identified the sub-genome structures of P. anserina, with A sub-genome was larger than B sub-genome and closer to P. micrantha phylogenetically. Despite lacking significant genome-wide expression dominance, the A sub-genome had higher homoeologous gene expression in shoot apical meristem, flower and tuberous root. The resistance genes was contracted in P. anserina genome. Key genes involved in starch biosynthesis were expanded and highly expressed in tuberous roots, which probably drives the tuber formation. The genomics and transcriptomics data generated in this study advance our understanding of the genomic landscape of P. anserina, and will accelerate genetic studies and breeding programs.

Analysis of the resistance of small peptides from Periplaneta americana to hydrogen peroxide-induced apoptosis in human ovarian granular cells based on RNA-seq

Apoptosis of ovarian granular cells is closely related with weakening fertility of women. Hence, resisting apoptosis of human ovarian granular cells is of important significance. According to previous studies, DAPI fluorescence staining experiment and Western Blot test of Caspase-3 demonstrate that small peptides from Periplaneta americana (SPPA) can improve hydrogen peroxide (H₂O₂) -induced apoptosis of human ovarian granular cells (KGN cells). However, the molecular mechanism of SPPA resistance against apoptosis of granular cells still remains unknown. In this study, key genes and signaling pathways for SPPA to resist H₂O₂-induced apoptosis of KGN cells were determined through transcriptome sequencing (RNA-seq). Experiments were divided into three groups, namely, the control group, H₂O₂ group and H₂O₂ + SPPA group. A total of 1196 differentially expressed genes (DEGs) were screened by comparing the control group and the H₂O₂ group, and 2805 DEGs were screened by comparing the H₂O₂ group and H₂O₂ + SPPA group. It is important to note that 87 overlapping genes were identified upregulating in H₂O₂ exposure, but downregulating in SPPA repair. Another 151 overlapping genes were identified downregulating in H₂O₂ exposure, but upregulating in SPPA repair. These 238 overlapping genes have significant enrichment in multiple KEGG pathways. Among them, 13 genes play significant roles in SPPA resistance process of cell apoptosis: EIF3D, RAN, UPF1 and EIF2B4 participate in RNA transport; ACTG1, SIPA1 and CTNND1 participate in Leukocyte transendothelial migration; S100A7, S100A9, RELA and IL17RE participate in IL-17 signaling pathway; BCL2L13, EIF2AK3 and RELA participate in Mitophapy-animal. Ten genes were selected for florescence quantitative PCR (qPCR) verification and the expression level was consistent with sequencing results. Finally, a control network of SPPA resistance against the H₂O₂-induced KGN cell apoptosis was built based on the target genes screened by the RNA-seq technology. This study provides a direction and some references to further understand the molecular mechanism of SPPA resistance against the H₂O₂-induced KGN cell apoptosis.

Hippocampal Transcriptome-Wide Association Study Reveals Correlations Between Impaired Glutamatergic Synapse Pathway and Age-Related Hearing Loss in BXD-Recombinant Inbred Mice

Age-related hearing loss (ARHL) is associated with cognitive dysfunction; however, the detailed underlying mechanisms remain unclear. The aim of this study is to investigate the potential underlying mechanism with a system genetics approach. A transcriptome-wide association study was performed on aged (12–32 months old) BXD mice strains. The hippocampus gene expression was obtained from 56 BXD strains, and the hearing acuity was assessed from 54 BXD strains. Further correlation analysis identified a total of 1,435 hearing-related genes in the hippocampus (p < 0.05). Pathway analysis of these genes indicated that the impaired glutamatergic synapse pathway is involved in ARHL (p = 0.0038). Further gene co-expression analysis showed that the expression level of glutamine synthetase (Gls), which is significantly correlated with ARHL (n = 26, r = −0.46, p = 0.0193), is a crucial regulator in glutamatergic synapse pathway and associated with learning and memory behavior. In this study, we present the first systematic evaluation of hippocampus gene expression pattern associated with ARHL, learning, and memory behavior. Our results provide novel potential molecular mechanisms involved in ARHL and cognitive dysfunction association.

Small RNA sequencing reveals various microRNAs involved in piperine biosynthesis in black pepper (Piper nigrum L.)

Background

Black pepper (Piper nigrum L.), an important and long-cultivated spice crop, is native to South India and grown in the tropics. Piperine is the main pungent and bioactive alkaloid in the berries of black pepper, but the molecular mechanism for piperine biosynthesis has not been determined. MicroRNAs (miRNAs), which are classical endogenous noncoding small RNAs, play important roles in regulating secondary metabolism in many species, but less is known regarding black pepper or piperine biosynthesis.

Results

To dissect the functions of miRNAs in secondary metabolism especially in piperine biosynthesis, 110 known miRNAs, 18 novel miRNAs and 1007 individual targets were identified from different tissues of black pepper by small RNA sequencing. qRT-PCR and 5′-RLM-RACE experiments were conducted to validate the reliability of the sequencing data and predicted targets. We found 3 miRNAs along with their targets including miR166-4CL, miR396-PER and miR397-CCR modules that are involved in piperine biosynthesis.

Conclusion

MiRNA regulation of secondary metabolism is a common phenomenon in plants. Our study revealed new miRNAs that regulate piperine biosynthesis, which are special alkaloids in the piper genus, and they might be useful for future piperine genetic improvement of black pepper.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08154-4.

Transcriptome and metabolome analysis to reveal major genes of saikosaponin biosynthesis in Bupleurum chinense

BACKGROUND

Bupleurum chinense DC. is a widely used traditional Chinese medicinal plant. Saikosaponins are the major bioactive constituents of B. chinense, but relatively little is known about saikosaponin biosynthesis. In the present study, we performed an integrated analysis of metabolic composition and the expressed genes involved in saikosaponin biosynthetic pathways among four organs (the root, flower, stem, and leaf) of B. chinense to discover the genes related to the saikosaponin biosynthetic pathway.

RESULTS

Transcript and metabolite profiles were generated through high-throughput RNA-sequencing (RNA-seq) data analysis and liquid chromatography tandem mass spectrometry, respectively. Evaluation of saikosaponin contents and transcriptional changes showed 152 strong correlations (P < 0.05) over 3 compounds and 77 unigenes. These unigenes belonged to eight gene families: the acetoacetyl CoA transferase (AACT) (6), HMG-CoA synthase (HMGS) (2), HMG-CoA reductase (HMGR) (2), mevalonate diphosphate decarboxylase (MVD) (1), 1-deoxy-D-xylulose-5-phosphate synthase (DXS) (3), farnesyl diphosphate synthase (FPPS) (11), β-amyrin synthase (β-AS) (13) and cytochrome P450 enzymes (P450s) (39) families.

CONCLUSIONS

Our results investigated the diversity of the saikosaponin triterpene biosynthetic pathway in the roots, stems, leaves and flowers of B. chinese by integrated transcriptomic and metabolomic analysis, implying that manipulation of P450s genes such as Bc95697 and Bc35434 might improve saikosaponin biosynthesis. This is a good candidate for the genetic improvement of this important medicinal plant.

Butyrate and Forskolin Augment Host Defense, Barrier Function, and Disease Resistance Without Eliciting Inflammation

Host defense peptides (HDPs) are an integral part of the innate immune system with both antimicrobial and immunomodulatory activities. Induction of endogenous HDP synthesis is being actively explored as an antibiotic-alternative approach to disease control and prevention. Butyrate, a short-chain fatty acid, and forskolin, a phytochemical, have been shown separately to induce HDP gene expression in human cells. Here, we investigated the ability of butyrate and forskolin to induce the expressions of chicken HDP genes and the genes involved in barrier function such as mucin 2 and claudin 1 both in vitro and in vivo. We further evaluated their efficacy in protecting chickens from Clostridium perfringens-induced necrotic enteritis. Additionally, we profiled the transcriptome and global phosphorylation of chicken HD11 macrophage cells in response to butyrate and forskolin using RNA sequencing and a kinome peptide array, respectively. Our results showed a strong synergy between butyrate and forskolin in inducing the expressions of several, but not all, HDP genes. Importantly, dietary supplementation of butyrate and a forskolin-containing plant extract resulted in significant alleviation of intestinal lesions and the C. perfringens colonization in a synergistic manner in a chicken model of necrotic enteritis. RNA sequencing revealed a preferential increase in HDP and barrier function genes with no induction of proinflammatory cytokines in response to butyrate and forskolin. The antiinflammatory and barrier protective properties of butyrate and forskolin were further confirmed by the kinome peptide array. Moreover, we demonstrated an involvement of inducible cAMP early repressor (ICER)-mediated negative feedback in HDP induction by butyrate and forskolin. Overall, these results highlight a potential for developing butyrate and forskolin, two natural products, as novel antibiotic alternatives to enhance intestinal health and disease resistance in poultry and other animals.

Association of exposure to deoxynivalenol with DNA methylation in white blood cells in children in China

Deoxynivalenol (DON) is a mycotoxin that commonly contaminates cereals worldwide. Dietary exposure to DON is a subject of great public health concern, but studies on the health effects of chronic exposure to DON are not available. In this study, we investigated the connection between DNA methylation levels and DON exposure in children. The DNA methylation status of white blood cells from 32 children aged 2~15 years old in Henan, China, was profiled. A total of 378 differentially methylated CpGs were identified between the high and low DON exposure groups, and 8 KEGG pathways were found to be significantly enriched among the differentially methylated genes. In addition, the quantitative methylation of EIF2AK4, EMID2 and GNASAS was analysed using the Sequenom MassARRAY platform. The results showed that the methylation level of EIF2AK4 was significantly different between the two groups, and the methylation levels were associated with exposure to DON. Conclusively, our study found that chronic exposure to DON during childhood could affect DNA methylation levels.

Impact of Bicarbonate-β-Lactam Exposures on Methicillin-Resistant Staphylococcus aureus (MRSA) Gene Expression in Bicarbonate-β-Lactam-Responsive vs. Non-Responsive Strains

Methicillin-resistant Staphylococcus aureus (MRSA) infections represent a difficult clinical treatment issue. Recently, a novel phenotype was discovered amongst selected MRSA which exhibited enhanced β-lactam susceptibility in vitro in the presence of NaHCO₃ (termed ‘NaHCO₃-responsiveness’). This increased β-lactam susceptibility phenotype has been verified in both ex vivo and in vivo models. Mechanistic studies to-date have implicated NaHCO₃-mediated repression of genes involved in the production, as well as maturation, of the alternative penicillin-binding protein (PBP) 2a, a necessary component of MRSA β-lactam resistance. Herein, we utilized RNA-sequencing (RNA-seq) to identify genes that were differentially expressed in NaHCO₃-responsive (MRSA 11/11) vs. non-responsive (COL) strains, in the presence vs. absence of NaHCO₃-β-lactam co-exposures. These investigations revealed that NaHCO₃ selectively repressed the expression of a cadre of genes in strain 11/11 known to be a part of the sigB-sarA-agr regulon, as well as a number of genes involved in the anchoring of cell wall proteins in MRSA. Moreover, several genes related to autolysis, cell division, and cell wall biosynthesis/remodeling, were also selectively impacted by NaHCO₃-OXA exposure in the NaHCO₃-responsive strain MRSA 11/11. These outcomes provide an important framework for further studies to mechanistically verify the functional relevance of these genetic perturbations to the NaHCO₃-responsiveness phenotype in MRSA.

Genomic Insight Into the Population Structure and Admixture History of Tai-Kadai-Speaking Sui People in Southwest China

Sui people, which belong to the Tai-Kadai-speaking family, remain poorly characterized due to a lack of genome-wide data. To infer the fine-scale population genetic structure and putative genetic sources of the Sui people, we genotyped 498,655 genome-wide single-nucleotide polymorphisms (SNPs) using SNP arrays in 68 Sui individuals from seven indigenous populations in Guizhou province and Guangxi Zhuang Autonomous Region in Southwest China and co-analyzed with available East Asians via a series of population genetic methods including principal component analysis (PCA), ADMIXTURE, pairwise Fst genetic distance, f-statistics, qpWave, and qpAdm. Our results revealed that Guangxi and Guizhou Sui people showed a strong genetic affinity with populations from southern China and Southeast Asia, especially Tai-Kadai- and Hmong-Mien-speaking populations as well as ancient Iron Age Taiwan Hanben, Gongguan individuals supporting the hypothesis that Sui people came from southern China originally. The indigenous Tai-Kadai-related ancestry (represented by Li), Northern East Asian-related ancestry, and Hmong-Mien-related lineage contributed to the formation processes of the Sui people. We identified the genetic substructure within Sui groups: Guizhou Sui people were relatively homogeneous and possessed similar genetic profiles with neighboring Tai-Kadai-related populations, such as Maonan. While Sui people in Yizhou and Huanjiang of Guangxi might receive unique, additional gene flow from Hmong-Mien-speaking populations and Northern East Asians, respectively, after the divergence within other Sui populations. Sui people could be modeled as the admixture of ancient Yellow River Basin farmer-related ancestry (36.2-54.7%) and ancient coastal Southeast Asian-related ancestry (45.3-63.8%). We also identified the potential positive selection signals related to the disease susceptibility in Sui people via integrated haplotype score (iHS) and number of segregating sites by length (nSL) scores. These genomic findings provided new insights into the demographic history of Tai-Kadai-speaking Sui people and their interaction with neighboring populations in Southern China.

MapZ deficiency leads to defects in the envelope structure and changes stress tolerance of Streptococcus mutans

Cell division is a central process in bacteria and a prerequisite for pathogenicity. Several proteins are involved in this process to ensure the accurate localization and proper function of the division machinery. In Streptococcus mutans, MapZ marks the division sites and position of the Z-ring to regulate cell division; however, whether MapZ deficiency can impair the cariogenic virulence of S. mutans remains unclear. Here, using a phenotypic assay and RNA-seq, we investigated the role of MapZ in cell envelope maintenance, biofilm formation, and stress tolerance in S. mutans. The results show that MapZ is important for normal cell shape and envelope structure, and its deletion causes abnormal septum structure and a thin cell wall. Subsequently, we found that the absence of MapZ leads to a greater level of cell death within 12 h biofilms, but it does not seem to affect biofilm architecture and accumulation. mapZ deletion also results in a decreased acid and osmotic stress tolerance. Furthermore, RNA-seq data reveal that MapZ deficiency causes changes in the expression levels of genes involved in transport systems, sugar metabolism, nature competence, and bacteriocin synthesis. Interestingly, we found that mapZ mutation renders S. mutans more sensitive to chlorhexidine. Taken together, our study suggests that MapZ plays a role in maintaining cell envelope structure and stress tolerance in S. mutans, showing a potential application as a drug target for caries prevention.

Leber's hereditary optic neuropathy: Current approaches and future perspectives on Mesenchymal stem cell-mediated rescue

Leber's Hereditary Optic Neuropathy (LHON) is an inherited optic nerve disorder. It is a mitochondrially inherited disease due to point mutation in the MT-ND1, MT-ND4, and MT-ND6 genes of mitochondrial DNA (mtDNA) coding for complex I subunit proteins. These mutations affect the assembly of the mitochondrial complex I and hence the electron transport chain leading to mitochondrial dysfunction and oxidative damage. Optic nerve cells like retinal ganglion cells (RGCs) are more sensitive to mitochondrial loss and oxidative damage which results in the progressive degeneration of RGCs at the axonal region of the optic nerve leading to bilateral vision loss. Currently, gene therapy using Adeno-associated viral vector (AAV) is widely studied for the therapeutics application in LHON. Our review highlights the application of cell-based therapy for LHON. Mesenchymal stem cells (MSCs) are known to rescue cells from the pre-apoptotic stage by transferring healthy mitochondria through tunneling nanotubes (TNT) for cellular oxidative function. Empowering the transfer of healthy mitochondria using MSCs may replace the mitochondria with pathogenic mutation and possibly benefit the cells from progressive damage. This review discusses the ongoing research in LHON and mitochondrial transfer mechanisms to explore its scope in inherited optic neuropathy.

Identification of Circular RNA Expression Profiles in White Adipocytes and Their Roles in Adipogenesis

Obesity and its related metabolic diseases have become great public health threats worldwide. Although accumulated evidence suggests that circRNA is a new type of non-coding RNAs regulating various physiological and pathological processes, little attention has been paid to the expression profiles and functions of circRNAs in white adipose tissue. In this study, 3,771 circRNAs were detected in three stages of white adipogenesis (preadipocyte, differentiating preadipocyte, and mature adipocyte) by RNA-seq. Experimental validation suggested that the RNA-seq results are highly reliable. We found that nearly 10% of genes which expressed linear RNAs in adipocytes could also generate circRNAs. In addition, 40% of them produced multiple circRNA isoforms. We performed correlation analysis and found that a great deal of circRNAs (nearly 50%) and their parental genes were highly correlated in expression levels. A total of 41 differential expression circRNAs (DECs) were detected during adipogenesis and an extremely high ratio of them (80%) were correlated with their parental genes, indicating these circRNAs may potentially play roles in regulating the expression of their parental genes. KEGG enrichment and GO annotation of the parental genes suggesting that the DECs may participate in several adipogenesis-related pathways. Following rigorous selection, we found that many up-regulated circRNAs contain multiple miRNAs binding sites, such as miR17, miR-30c, and miR-130, indicating they may potentially facilitate their regulatory functions by acting as miRNA sponges. These results suggest that plenty of circRNAs are expressed in white adipogenesis and the DECs may serve as new candidates for future adipogenesis regulation.

Workflows for Rapid Functional Annotation of Diverse Arthropod Genomes

Genome sequencing of a diverse array of arthropod genomes is already underway, and these genomes will be used to study human health, agriculture, biodiversity, and ecology. These new genomes are intended to serve as community resources and provide the foundational information required to apply 'omics technologies to a more diverse set of species. However, biologists require genome annotation to use these genomes and derive a better understanding of complex biological systems. Genome annotation incorporates two related, but distinct, processes: Demarcating genes and other elements present in genome sequences (structural annotation); and associating a function with genetic elements (functional annotation). While there are well-established and freely available workflows for structural annotation of gene identification in newly assembled genomes, workflows for providing the functional annotation required to support functional genomics studies are less well understood. Genome-scale functional annotation is required for functional modeling (enrichment, networks, etc.). A first-pass genome-wide functional annotation effort can rapidly identify under-represented gene sets for focused community annotation efforts. We present an open-source, open access, and containerized pipeline for genome-scale functional annotation of insect proteomes and apply it to various arthropod species. We show that the performance of the predictions is consistent across a set of arthropod genomes with varying assembly and annotation quality.

LncGSEA: a versatile tool to infer lncRNA associated pathways from large-scale cancer transcriptome sequencing data

Background

Long non-coding RNAs (lncRNAs) are a growing focus in cancer research. Deciphering pathways influenced by lncRNAs is important to understand their role in cancer. Although knock-down or overexpression of lncRNAs followed by gene expression profiling in cancer cell lines are established approaches to address this problem, these experimental data are not available for a majority of the annotated lncRNAs.

Results

As a surrogate, we present lncGSEA, a convenient tool to predict the lncRNA associated pathways through Gene Set Enrichment Analysis of gene expression profiles from large-scale cancer patient samples. We demonstrate that lncGSEA is able to recapitulate lncRNA associated pathways supported by literature and experimental validations in multiple cancer types.

Conclusions

LncGSEA allows researchers to infer lncRNA regulatory pathways directly from clinical samples in oncology. LncGSEA is written in R, and is freely accessible at https://github.com/ylab-hi/lncGSEA.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07900-y.

Comparative Transcriptome Analysis of Wheat Lines in the Field Reveals Multiple Essential Biochemical Pathways Suppressed by Obligate Pathogens

Mildew and rust are the most devastating cereal pathogens, and in wheat they can cause up to 50% yield loss every year. Wheat lines containing resistance genes are used to effectively control fungal diseases, but the molecular mechanisms underlying the interaction between wheat and its fungal pathogens are poorly understood. Here, we used RNA sequencing (RNA-Seq) to compare the transcriptomic landscape of susceptible and resistant wheat lines to identify genes and pathways that are targeted by obligate biotrophic fungal pathogens. The five lines differed in the expression of thousands of genes under infection as well as control conditions. Generally, mixed infection with powdery mildew and leaf rust resulted in downregulation of numerous genes in susceptible lines. Interestingly, transcriptomic comparison between the nearly isogenic lines Thatcher and Thatcher-Lr34 identified 753 genes that are uniquely downregulated in the susceptible line upon infection. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis, revealed the suppression of six major biochemical pathways, namely nuclear transport, alternative splicing, DNA damage response, ubiquitin-mediated proteolysis, phosphoinositol signaling, and photosynthesis. We conclude that powdery mildew and leaf rust evade the wheat defense system by suppression of programmed cell death (PCD) and responses to cellular damage. Considering the broad range of the induced changes, we propose that the pathogen targets "master regulators" at critical steps in the respective pathways. Identification of these wheat genes targeted by the pathogen could inspire new directions for future wheat breeding.