Coexpression network analysis of the genes regulated by two types of resistance responses to powdery mildew in wheat

StCoExpNet: a global co-expression network analysis facilitates identifying genes underlying agronomic traits in potatoes

KEY MESSAGE

We constructed a gene expression atlas and co-expression network for potatoes and identified several novel genes associated with various agronomic traits. This resource will accelerate potato genetics and genomics research. Potato (Solanum tuberosum L.) is the world's most crucial non-cereal food crop and ranks third in food production after wheat and rice. Despite the availability of several potato transcriptome datasets at public databases like NCBI SRA, an effort has yet to be put into developing a global transcriptome atlas and a co-expression network for potatoes. The objectives of our study were to construct a global expression atlas for potatoes using publicly available transcriptome datasets, identify housekeeping and tissue-specific genes, construct a global co-expression network and identify co-expression clusters, investigate the transcriptional complexity of genes involved in various essential biological processes related to agronomic traits, and provide a web server (StCoExpNet) to easily access the newly constructed expression atlas and co-expression network to investigate the expression and co-expression of genes of interest. In this study, we used data from 2299 publicly available potato transcriptome samples obtained from 15 different tissues to construct a global transcriptome atlas. We found that roughly 87% of the annotated genes exhibited detectable expression in at least one sample. Among these, we identified 281 genes with consistent and stable expression levels, indicating their role as housekeeping genes. Conversely, 308 genes exhibited marked tissue-specific expression patterns. We exemplarily linked some co-expression clusters to important agronomic traits of potatoes, such as self-incompatibility, anthocyanin biosynthesis, tuberization, and defense responses against multiple pathogens. The dataset compiled here constitutes a new resource (StCoExpNet), which can be accessed at https://stcoexpnet.julius-kuehn.de . This transcriptome atlas and the co-expression network will accelerate potato genetics and genomics research.

PlantNexus: A Gene Co-expression Network Database and Visualization Tool for Barley and Sorghum

Global gene co-expression networks (GCNs) are powerful tools for functional genomics whereby putative functions and regulatory mechanisms can be inferred by gene co-expression. Cereal crops, such as Hordeum vulgare (barley) and Sorghum bicolor (sorghum), are among the most important plants to civilization. However, co-expression network tools for these plants are lacking. Here, we have constructed global GCNs for barley and sorghum using existing RNA-seq data sets. Meta-information was manually curated and categorized by tissue type to also build tissue-specific GCNs. To enable GCN searching and visualization, we implemented a website and database named PlantNexus. PlantNexus is freely available at https://plantnexus.ohio.edu/.

The decreased expression of GW2 homologous genes contributed to the increased grain width and thousand‑grain weight in wheat-Dasypyrum villosum 6VS·6DL translocation lines

This study demonstrated that the aberrant transcription of DvGW2 contributed to the increased grain width and thousand-grain weight in wheat-Dasypyrum villosum T6VS·6DL translocation lines. Due to the high immunity to powdery mildew, Dasypyrum villosum 6VS has been one of the most successful applications of the wild relatives in modern wheat breeding. Along with the desired traits, side-effects could be brought when large alien chromosome fragments are introduced into wheat, but little is known about effects of 6VS on agronomic traits. Here, we found that T6VS·6DL translocation had significantly positive effects on grain weight, plant heightand spike length, and small negative effects on total spikelet number and spikelet compactness using recipient and wheat-D. villosum T6VS·6DL allohexaploid wheats, Wan7107 and Pm97033. Further analysis showed that the 6VS segment might exert direct genetic effect on grain width, then driving the increase of thousand-grain weight. Furthermore, comparative transcriptome analysis identified 2549 and 1282 differentially expressed genes (DEGs) and 2220 and 1496 specifically expressed genes (SEGs) at 6 days after pollination (DAP) grains and 15 DAP endosperms, respectively. Enrichment analysis indicated that the process of cell proliferation category was over-represented in the DEGs. Notably, two homologous genes, TaGW2-D1 and DvGW2, were identified as putative candidate genes associated with grain weight and yield. The expression analysis showed that DvGW2 had an aberrant expression in Pm97033, resulting in significantly lower total expression level of GW2 than Wan7107, which drives the increase of grain weight and width in Pm97033. Collectively, our data indicated that the compromised expression of DvGW2 is critical for increased grain width and weight in T6VS·6DL translocation lines.

Network biology to uncover functional and structural properties of the plant immune system

In the last two decades, advances in network science have facilitated the discovery of important systems' entities in diverse biological networks. This graph-based technique has revealed numerous emergent properties of a system that enable us to understand several complex biological processes including plant immune systems. With the accumulation of multiomics data sets, the comprehensive understanding of plant-pathogen interactions can be achieved through the analyses and efficacious integration of multidimensional qualitative and quantitative relationships among the components of hosts and their microbes. This review highlights comparative network topology analyses in plant-pathogen co-expression networks and interactomes, outlines dynamic network modeling for cell-specific immune regulatory networks, and discusses the new frontiers of single-cell sequencing as well as multiomics data integration that are necessary for unraveling the intricacies of plant immune systems.

Efficient expression and function of a receptor-like kinase in wheat powdery mildew defence require an intron-located MYB binding site

The LRK10-like receptor kinases (LRK10L-RLKs) are ubiquitously present in higher plants, but knowledge of their expression and function is still limited. Here, we report expression and functional analysis of TtdLRK10L-1, a typical LRK10L-RLK in durum wheat (Triticum turgidum L. ssp. durum). The introns of TtdLRK10L-1 contained multiple kinds of predicted cis-elements. To investigate the potential effect of these cis-elements on TtdLRK10L-1 expression and function, two types of transgenic wheat lines were prepared, which expressed a GFP-tagged TtdLRK10L-1 protein (TtdLRK10L-1:GFP) from the cDNA or genomic DNA (gDNA) sequence of TtdLRK10L-1 under the native promoter. TtdLRK10L-1:GFP expression was up-regulated by the powdery mildew pathogen Blumeria graminis f. sp. tritici (Bgt) in both types of transgenic plants, with the scale of the elevation being much stronger in the gDNA lines. Both types of transgenic plants exhibited enhanced resistance to Bgt infection relative to wild type control. Notably, the Bgt defence activated in the gDNA lines was significantly stronger than that in the cDNA lines. Further analysis revealed that a putative MYB transcription factor binding site (MYB-BS, CAGTTA) located in TtdLRK10L-1 intron I was critical for the efficient expression and function of TtdLRK10L-1 in Bgt defence. This MYB-BS could also increase the activity of a superpromoter widely used in ectopic gene expression studies in plants. Together, our results deepen the understanding of the expression and functional characteristics of LRK10L-RLKs. TtdLRK10L-1 is likely useful for further dissecting the molecular processes underlying wheat defence against Bgt and for developing Bgt resistant wheat crops.

Structure and predictive metabolic contribution of intestinal microbiota of Longfin yellowtail (Seriola rivoliana) juveniles in aquaculture systems

Seriola rivoliana intestinal microbiota (IM) was characterised under aquaculture conditions through 16S rRNA amplicon sequencing. Specimens of 30 days after hatching (DAH) were maintained in three tanks and fed under the same environmental conditions for characterisation 15 days prior to sampling. Three fish were randomly taken from each tank; total DNA extraction of the gut microbiota was performed to characterise microbial composition and its metabolic prediction. The V3 hypervariable region of the 16S rRNA was amplified and sequenced with Illumina pair-end technology. The prokaryotic components in the S. rivoliana intestine were dominated mainly by the phyla Proteobacteria, Firmicutes, Bacteroidetes, Cyanobacteria and Actinobacteria. No significant differences in beta diversity were detected in the three samples (tanks). However in alpha diversity, they were detected in juveniles of the same cohort within the same group, as exemplified by enrichment of certain bacterial groups, mainly of the Clostridia class, which were specific in each fish within the same tank. The metabolic prediction analyses suggested that S. rivoliana IM contribute to the metabolism of amino acids, carbohydrates, lipids, and immune system. This study provides the first IM characterisation under rearing conditions of S. rivoliana-a species with broad economic potential-and contributes to novel information for potential use of probiotics in future trials.

Epigenetic Activation of Enoyl-CoA Reductase By An Acetyltransferase Complex Triggers Wheat Wax Biosynthesis

The epidermal surface of bread wheat (Triticum aestivum) is coated with a hydrophobic cuticular wax layer that protects plant tissues against environmental stresses. However, the regulatory mechanism of cuticular wax biosynthesis remains to be uncovered in bread wheat. Here, we identified wheat Enoyl-CoA Reductase (TaECR) as a core component responsible for biosynthesis of wheat cuticular wax. Silencing of TaECR in bread wheat resulted in a reduced cuticular wax load and attenuated conidia germination of the adapted fungal pathogen powdery mildew (Blumeria graminis f.sp. tritici). Furthermore, we established that TaECR genes are direct targets of TaECR promoter-binding MYB transcription factor1 (TaEPBM1), which could interact with the adapter protein Alteration/Deficiency in Activation2 (TaADA2) and recruit the histone acetyltransferase General Control Nonderepressible5 (TaGCN5) to TaECR promoters. Most importantly, we demonstrated that the TaEPBM1-TaADA2-TaGCN5 ternary protein complex activates TaECR transcription by potentiating histone acetylation and enhancing RNA polymerase II enrichment at TaECR genes, thereby contributing to the wheat cuticular wax biosynthesis. Finally, we identified very-long-chain aldehydes as the wax signals provided by the TaECR-TaEPBM1-TaADA2-TaGCN5 circuit for triggering B graminis f.sp. tritici conidia germination. These results demonstrate that specific transcription factors recruit the TaADA2-TaGCN5 histone acetyltransferase complex to epigenetically regulate biosynthesis of wheat cuticular wax, which is required for triggering germination of the adapted powdery mildew pathogen.

A spontaneous wheat-Aegilops longissima translocation carrying Pm66 confers resistance to powdery mildew

A spontaneous Robertsonian T4S^lS·4BL translocation chromosome carrying Pm66 for powdery mildew resistance was discovered and confirmed by RNA-seq, molecular marker, and in situ hybridization analyses. Powdery mildew caused by Blumeria graminis f. sp. tritici (Bgt) is a severe disease of bread wheat worldwide. Discovery and utilization of resistance genes to powdery mildew from wild relatives of wheat have played important roles in wheat improvement. Aegilops longissima, one of the S-genome diploid wild relatives of wheat, is a valuable source of disease and pest resistance for wheat. Chromosome 4S^l from Ae. longissima confers moderate resistance to powdery mildew. In this study, we conducted RNA-seq on a putative Chinese Spring (CS)-Ae. longissima 4S^l(4B) disomic substitution line (TA3465) to develop 4S^l-specific markers to assist the transfer of a Bgt resistance gene from 4S^l by induced homoeologous recombination. A pairwise comparison of genes between CS and TA3465 demonstrated that a number of genes on chromosome 4BS in CS were not expressed in TA3465. Analysis of 4B- and 4S^l-specific molecular markers showed that 4BS and 4S^lL were both missing in TA3465, whereas 4BL and 4S^lS were present. Further characterization by genomic and fluorescent in situ hybridization confirmed that TA3465 carried a spontaneous Robertsonian T4S^lS·4BL translocation. Powdery mildew tests showed that TA3465 was resistant to 10 of 16 Bgt isolates collected from different regions of China, whereas CS was susceptible to all those Bgt isolates. The powdery mildew resistance gene(s) in TA3465 was further mapped to the short arm of 4S^l and designated as Pm66.

The TuMYB46L-TuACO3 module regulates ethylene biosynthesis in einkorn wheat defense to powdery mildew

Powdery mildew disease, elicited by the obligate fungal pathogen Blumeria graminis f.sp. tritici (Bgt), causes widespread yield losses in global wheat crop. However, the molecular mechanisms governing wheat defense to Bgt are still not well understood. Here we found that TuACO3, encoding the 1-aminocyclopropane-1-carboxylic acid (ACC) oxidase functioning in ethylene (ET) biosynthesis, was induced by Bgt infection of the einkorn wheat Triticum urartu, which was accompanied by increased ET content. Silencing TuACO3 decreased ET production and compromised wheat defense to Bgt, whereas both processes were enhanced in the transgenic wheat overexpressing TuACO3. TuMYB46L, phylogenetically related to Arabidopsis MYB transcription factor AtMYB46, was found to bind to the TuACO3 promoter region in yeast-one-hybrid and EMSA experiments. TuMYB46L expression decreased rapidly following Bgt infection. Silencing TuMYB46L promoted ET content and Bgt defense, but the reverse was observed when TuMYB46L was overexpressed. Hence, decreased expression of TuMYB46L permits elevated function of TuACO3 in ET biosynthesis in Bgt-infected wheat. The TuMYB46L-TuACO3 module regulates ET biosynthesis to promote einkorn wheat defense against Bgt. Furthermore, we found four chitinase genes acting downstream of the TuMYB46L-TuACO3 module. Collectively, our data shed a new light on the molecular mechanisms underlying wheat defense to Bgt.

Physical Mapping of Pm57, a Powdery Mildew Resistance Gene Derived from Aegilops searsii

Powdery mildew caused by Blumeria graminis f. sp. tritici (Bgt) is one of many severe diseases that threaten bread wheat (Triticum aestivum L.) yield and quality worldwide. The discovery and deployment of powdery mildew resistance genes (Pm) can prevent this disease epidemic in wheat. In a previous study, we transferred the powdery mildew resistance gene Pm57 from Aegilops searsii into common wheat and cytogenetically mapped the gene in a chromosome region with the fraction length (FL) 0.75-0.87, which represents 12% segment of the long arm of chromosome 2Ss#1. In this study, we performed RNA-seq using RNA extracted from leaf samples of three infected and mock-infected wheat-Ae. searsii 2Ss#1 introgression lines at 0, 12, 24, and 48 h after inoculation with Bgt isolates. Then we designed 79 molecular markers based on transcriptome sequences and physically mapped them to Ae. searsii chromosome 2Ss#1- in seven intervals. We used these markers to identify 46 wheat-Ae. searsii 2Ss#1 recombinants induced by ph1b, a deletion mutant of pairing homologous (Ph) genes. After analyzing the 46 ph1b-induced 2Ss#1L recombinants in the region where Pm57 is located with different Bgt-responses, we physically mapped Pm57 gene on the long arm of 2Ss#1 in a 5.13 Mb genomic region, which was flanked by markers X67593 (773.72 Mb) and X62492 (778.85 Mb). By comparative synteny analysis of the corresponding region on chromosome 2B in Chinese Spring (T. aestivum L.) with other model species, we identified ten genes that are putative plant defense-related (R) genes which includes six coiled-coil nucleotide-binding site-leucine-rich repeat (CNL), three nucleotide-binding site-leucine-rich repeat (NL) and a leucine-rich receptor-like repeat (RLP) encoding proteins. This study will lay a foundation for cloning of Pm57, and benefit the understanding of interactions between resistance genes of wheat and powdery mildew pathogens.

Identification of key genes involved in the biosynthesis of triterpenic acids in the mint family

Triterpenic acids (TAs), a large group of natural compounds with diverse biological activity, are produced by several plant taxa. Betulinic, oleanolic, and ursolic acids are the most medicinally important TAs and are mainly found in plants of the mint family. Metabolic engineering is strongly dependent on identifying the key genes in biosynthetic pathways toward the products of interest. In this study, gene expression tracking was performed by transcriptome mining, co-expression network analysis, and tissue-specific metabolite-expression analysis in order to identify possible key genes involved in TAs biosynthetic pathways. To this end, taxa-specific degenerate primers of six important genes were designed using an effective method based on the MEME algorithm in a phylogenetically related group of sequences and successfully applied in three members of the Lamiaceae (Rosmarinus officinalis, Salvia officinalis, and Thymus persicus). Based on the results of in-depth data analysis, genes encoding squalene epoxidase and oxido squalene cyclases are proposed as targets for boosting triterpene production. The results emphasize the importance of identifying key genes in triterpene biosynthesis, which may facilitate genetic manipulation or overexpression of target genes.

Genome-Wide Association and Gene Co-expression Network Analyses Reveal Complex Genetics of Resistance to Goss's Wilt of Maize

Goss’s bacterial wilt and leaf blight is a disease of maize caused by the gram positive bacterium Clavibacter michiganensis subsp. nebraskensis (Cmn). First discovered in Nebraska, Goss’s wilt has now spread to major maize growing states in the United States and three provinces in Canada. Previous studies conducted using elite maize inbred lines and their hybrids have shown that resistance to Goss’s wilt is a quantitative trait. The objective of this study was to further our understanding of the genetic basis of resistance to Goss’s wilt by using a combined approach of genome-wide association mapping and gene co-expression network analysis. Genome-wide association analysis was accomplished using a diversity panel consisting of 555 maize inbred lines and a set of 450 recombinant inbred lines (RILs) from three bi-parental mapping populations, providing the most comprehensive screening of Goss’s wilt resistance to date. Three SNPs in the diversity panel and 10 SNPs in the combined dataset, including the diversity panel and RILs, were found to be significantly associated with Goss’s wilt resistance. Each significant SNP explained 1–5% of the phenotypic variation for Goss’s wilt (total of 8–11%). To augment the results of genome-wide association mapping and help identify candidate genes, a time course RNA sequencing experiment was conducted using resistant (N551) and susceptible (B14A) maize inbred lines. Gene co-expression network analysis of this time course experiment identified one module of 141 correlated genes that showed differential regulation in response to Cmn inoculations in both resistant and susceptible lines. SNPs inside and flanking these genes explained 13.3% of the phenotypic variation. Among 1,000 random samples of genes, only 8% of samples explained more phenotypic variance for Goss’s wilt resistance than those implicated by the co-expression network analysis. While a statistically significant enrichment was not observed (P < 0.05), these results suggest a possible role for these genes in quantitative resistance at the field level and warrant more research on combining gene co-expression network analysis with quantitative genetic analyses to dissect complex disease resistance traits. The results of the GWAS and co-expression analysis both support the complex nature of resistance to this important disease of maize.

Development of PCR markers specific to Dasypyrum villosum genome based on transcriptome data and their application in breeding Triticum aestivum-D. villosum#4 alien chromosome lines

Background

Dasypyrum villosum is an important wild species of wheat (Triticum aestivum L.) and harbors many desirable genes that can be used to improve various traits of wheat. Compared with other D. villosum accessions, D. villosum#4 still remains less studied. In particular, chromosomes of D. villosum#4 except 6V#4 have not been introduced into wheat by addition or substitution and translocation, which is an essential step to identify and apply the alien desired genes. RNA-seq technology can generate large amounts of transcriptome sequences and accelerate the development of chromosome-specific molecular markers and assisted selection of alien chromosome line.

Results

We obtained the transcriptome of D. villosum#4 via a high-throughput sequencing technique, and then developed 76 markers specific to each chromosome arm of D. villosum#4 based on the bioinformatic analysis of the transcriptome data. The D. villosum#4 sequences containing the specific DNA markers were expected to be involved in different genes, among which most had functions in metabolic processes. Consequently, we mapped these newly developed molecular markers to the homologous chromosome of barley and obtained the chromosome localization of these markers on barley genome. Then we analyzed the collinearity of these markers among D. villosum, wheat, and barley. In succession, we identified six types of T. aestivum-D. villosum#4 alien chromosome lines which had one or more than one D. villosum#4 chromosome in the cross and backcross BC₃F₅ populations between T. durum–D. villosum#4 amphidiploid TH3 and wheat cv. Wan7107 by employing the selected specific markers, some of which were further confirmed to be translocation or addition lines by genomic in situ hybridization (GISH).

Conclusion

Seventy-six PCR markers specific to chromosomes of D. villosum#4 based on transcriptome data were developed in the current study and their collinearity among D. villosum, wheat, and barley were carried out. Six types of Triticum aestivum-D. villosum#4 alien chromosome lines were identified by using 12 developed markers and some of which were further confirmed by GISH. These novel T. aestivum-D. villosum#4 chromosome lines have great potential to be used for the introduction of desirable genes from D. villosum#4 into wheat by chromosomal translocation to breed new wheat varieties.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5630-4) contains supplementary material, which is available to authorized users.

Wheat WD40-repeat protein TaHOS15 functions in a histone deacetylase complex to fine-tune defense responses to Blumeria graminis f.sp. tritici

Powdery mildew caused by Blumeria graminis f.sp. tritici (Bgt) seriously threatens the production of common wheat (Triticum aestivum). In eukaryotes, WD40-repeat (WDR) proteins usually participate in assembling protein complexes involved in a wide range of cellular processes, including defense responses. However, the potential function of WDR proteins in regulating crop resistance to biotrophic fungal pathogens, such as Bgt, remains unclear. In this study, we isolated TaHOS15, encoding a WDR protein, from the Bgt-susceptible wheat cultivar Jing411 and demonstrated that knockdown of TaHOS15 expression using virus- or transient-induced gene-silencing attenuated wheat susceptibility to Bgt. Biochemical and molecular-biological assays revealed that TaHOS15 interacts with TaHDA6, a wheat homolog of Arabidopsis histone deacetylase AtHDA6, to constitute a transcriptional repressor complex. We determined the role of TaHOS15, which might act as an adaptor protein recruiting TaHDA6 to the chromatin of wheat defense-related genes including TaPR1, TaPR2, TaPR5, and TaWRKY45, where they repress histone acetylation. Reduced TaHOS15 or TaHDA6 transcript levels led to decreased susceptibility to Bgt together with enhanced defense-related transcription under Bgt infection. Collectively, these results demonstrate that TaHOS15 functions in a histone deacetylase complex with TaHDA6 to fine-tune the defense response to Bgt in common wheat.

High-throughput sequencing analysis of lncRNAs in hippocampus tissues with hypoxic-ischemic brain damage

LncRNAs abundantly expressed in the brain have vital and wide-ranging functions in different biological processes. However, little is currently known regarding the influence of lncRNAs in developing brains after hypoxic-ischemic brain damage (HIBD). In this study, to investigate the lncRNAs expression signatures and the co-expression network of lncRNAs and mRNAs in the brain after HIBD, we established a neonatal rat HIBD model and detected the expression profiles of lncRNAs in the HIBD brain and a sham control using high-throughput sequencing. Further, highly differentially expressed lncRNAs were selected and validated by qRT-PCR. Finally, the biological functions of the selected lncRNAs were investigated by over-expressing or silencing the target genes through lentivirus transfection in hippocampal neuron cells. Our results revealed that the expression profile of lncRNAs was dramatically different between the HIBD brains and the sham control, showing as the aberrant expression of 617 lncRNA transcripts and 441 mRNA transcripts at 24 hours after HIBD. GO and KEGG analyses indicated that the differentially expressed mRNAs were mostly involved in the apoptosis signaling pathway. After validating the expression of 8 randomly selected lncRNA transcripts by qRT-PCR, we found that the TNFRSF17 gene (ID: ENSRNOG00000021987) was down-regulated in HI brains. After stable over-expression and silencing of TNFRSF17, the apoptosis rate of hippocampal neuron cells exhibited obvious changes under hypoxia or normaxia. The over-expression of TNFRSF17 could significantly up-regulate Bcl-2 but down-regulate Bax, caspase-3, and caspase-9 at the mRNA and protein levels, while the silencing of TNFRSF17 led to just the opposite phenomenon. Notably, the regulation effects of TNFRSF17 on apoptotic related genes and proteins under hypoxia were more obvious than those under normaxia. Moreover, the over-expression of TNFRSF17 reduced the apoptotic rate, but the loss of TNFRSF17 led to a high rate of apoptosis under hypoxia. Taken together, the silencing of TNFRSF17 exacerbated, while over-expression attenuated, neuron apoptosis induced by HI injury, suggesting that TNFRSF17 may be a target for the prognosis, diagnosis, and treatment of HIBD.

Comparison of gene co-networks analysis provide a systems view of rice (Oryza sativa L.) response to Tilletia horrida infection

The biotrophic soil-borne fungus Tilletia horrida causes rice kernel smut, an important disease affecting the production of rice male sterile lines in most hybrid rice growing regions of the world. There are no successful ways of controlling this disease and there has been little study of mechanisms of resistance to T. horrida. Based on transcriptional data of different infection time points, we found 23, 782 and 23, 718 differentially expressed genes (fragments per kilobase of transcript sequence per million, FPKM >1) in Jiangcheng 3A (resistant to T. horrida) and 9311A (susceptible to T. horrida), respectively. In order to illuminate the differential responses of the two rice male sterile lines to T. horrida, we identified gene co-expression modules using the method of weighted gene co-expression network analysis (WGCNA) and compared the different biological functions of gene co-expression networks in key modules at different infection time points. The results indicated that gene co-expression networks in the two rice genotypes were different and that genes contained in some modules of the two groups may play important roles in resistance to T. horrida, such as DTH8 and OsHop/Sti1a. Furthermore, these results provide a global view of the responses of two different phenotypes to T. horrida, and assist our understanding of the regulation of expression changes after T. horrida infection.

Gene coexpression network analysis combined with metabonomics reveals the resistance responses to powdery mildew in Tibetan hulless barley

Powdery mildew is a fungal disease that represents a ubiquitous threat to crop plants. Transcriptomic and metabolomic analyses were used to identify molecular and physiological changes in Tibetan hulless barley in response to powdery mildew. There were 3418 genes and 405 metabolites differentially expressed between the complete resistance cultivar G7 and the sensitive cultivar Z13. Weighted gene coexpression network analysis was carried out, and the differentially expressed genes were enriched in five and four major network modules in G7 and Z13, respectively. Further analyses showed that phytohormones, photosynthesis, phenylpropanoid biosynthesis, and flavonoid biosynthesis pathways were altered during Qingke-Blumeria graminis (DC.) f.sp. hordei (Bgh) interaction. Comparative analyses showed a correspondence between gene expression and metabolite profiles, and the activated defenses resulted in changes of metabolites involved in plant defense response, such as phytohormones, lipids, flavone and flavonoids, phenolamides, and phenylpropanoids. This study enabled the identification of Bgh responsive genes and provided new insights into the dynamic physiological changes that occur in Qingke during response to powdery mildew. These findings greatly improve our understanding of the mechanisms of induced defense response in Qingke and will provide new clues for the development of resistant Tibetan hulless barley varieties.

Comparison of gene co-networks reveals the molecular mechanisms of the rice (Oryza sativa L.) response to Rhizoctonia solani AG1 IA infection

Rhizoctonia solani causes rice sheath blight, an important disease affecting the growth of rice (Oryza sativa L.). Attempts to control the disease have met with little success. Based on transcriptional profiling, we previously identified more than 11,947 common differentially expressed genes (TPM > 10) between the rice genotypes TeQing and Lemont. In the current study, we extended these findings by focusing on an analysis of gene co-expression in response to R. solani AG1 IA and identified gene modules within the networks through weighted gene co-expression network analysis (WGCNA). We compared the different genes assigned to each module and the biological interpretations of gene co-expression networks at early and later modules in the two rice genotypes to reveal differential responses to AG1 IA. Our results show that different changes occurred in the two rice genotypes and that the modules in the two groups contain a number of candidate genes possibly involved in pathogenesis, such as the VQ protein. Furthermore, these gene co-expression networks provide comprehensive transcriptional information regarding gene expression in rice in response to AG1 IA. The co-expression networks derived from our data offer ideas for follow-up experimentation that will help advance our understanding of the translational regulation of rice gene expression changes in response to AG1 IA.

Genome sequence of the progenitor of wheat A subgenome Triticum urartu

Triticum urartu (diploid, AA) is the progenitor of the A subgenome of tetraploid (Triticum turgidum, AABB) and hexaploid (Triticum aestivum, AABBDD) wheat^1,2. Genomic studies of T. urartu have been useful for investigating the structure, function and evolution of polyploid wheat genomes. Here we report the generation of a high-quality genome sequence of T. urartu by combining bacterial artificial chromosome (BAC)-by-BAC sequencing, single molecule real-time whole-genome shotgun sequencing ³ , linked reads and optical mapping^4,5. We assembled seven chromosome-scale pseudomolecules and identified protein-coding genes, and we suggest a model for the evolution of T. urartu chromosomes. Comparative analyses with genomes of other grasses showed gene loss and amplification in the numbers of transposable elements in the T. urartu genome. Population genomics analysis of 147 T. urartu accessions from across the Fertile Crescent showed clustering of three groups, with differences in altitude and biostress, such as powdery mildew disease. The T. urartu genome assembly provides a valuable resource for studying genetic variation in wheat and related grasses, and promises to facilitate the discovery of genes that could be useful for wheat improvement.

Genome sequence of the progenitor of wheat A subgenome Triticum urartu

Triticum urartu (diploid, AA) is the progenitor of the A subgenome of tetraploid (Triticum turgidum, AABB) and hexaploid (Triticum aestivum, AABBDD) wheat^1,2. Genomic studies of T. urartu have been useful for investigating the structure, function and evolution of polyploid wheat genomes. Here we report the generation of a high-quality genome sequence of T. urartu by combining bacterial artificial chromosome (BAC)-by-BAC sequencing, single molecule real-time whole-genome shotgun sequencing³, linked reads and optical mapping^4,5. We assembled seven chromosome-scale pseudomolecules and identified protein-coding genes, and we suggest a model for the evolution of T. urartu chromosomes. Comparative analyses with genomes of other grasses showed gene loss and amplification in the numbers of transposable elements in the T. urartu genome. Population genomics analysis of 147 T. urartu accessions from across the Fertile Crescent showed clustering of three groups, with differences in altitude and biostress, such as powdery mildew disease. The T. urartu genome assembly provides a valuable resource for studying genetic variation in wheat and related grasses, and promises to facilitate the discovery of genes that could be useful for wheat improvement.

The genome sequence of Triticum urartu, the progenitor of the A subgenome of hexaploid wheat, provides insight into genome duplication during grass evolution.

The NB-LRR gene Pm60 confers powdery mildew resistance in wheat

Powdery mildew is one of the most devastating diseases of wheat. To date, few powdery mildew resistance genes have been cloned from wheat due to the size and complexity of the wheat genome. Triticum urartu is the progenitor of the A genome of wheat and is an important source for powdery mildew resistance genes. Using molecular markers designed from scaffolds of the sequenced T. urartu accession and standard map-based cloning, a powdery mildew resistance locus was mapped to a 356-kb region, which contains two nucleotide-binding and leucine-rich repeat domain (NB-LRR) protein-encoding genes. Virus-induced gene silencing, single-cell transient expression, and stable transformation assays demonstrated that one of these two genes, designated Pm60, confers resistance to powdery mildew. Overexpression of full-length Pm60 and two allelic variants in Nicotiana benthamiana leaves induced hypersensitive cell death response, but expression of the coiled-coil domain alone was insufficient to induce hypersensitive response. Yeast two-hybrid, bimolecular fluorescence complementation and luciferase complementation imaging assays showed that Pm60 protein interacts with its neighboring NB-containing protein, suggesting that they might be functionally related. The identification and cloning of this novel wheat powdery mildew resistance gene will facilitate breeding for disease resistance in wheat.

Development and comparative genomic mapping of Dasypyrum villosum 6V#4S-specific PCR markers using transcriptome data

Twenty-five Dasypyrum villosum 6V#4S-specific PCR markers were developed using transcriptome data and further assigned to comparative genomic maps of wheat chromosome 6A, 6B, and 6D and barley chromosome 6H contrasting their homologous genes in these genomes. Two Dasypyrum villosum accessions, D.v#2 and No. 1026 from England and Russia, respectively, contain Pm21 on chromosome 6V#2S and PmV on chromosome 6V#4S. Both genes confer high resistance to powdery mildew (PM) in wheat. Even though several molecular markers have been developed to detect Pm21 and PmV, only the MBH1 marker can simultaneously detect both Pm21 and PmV. In this study, we first used a high-throughput sequencing technique to obtain the transcriptome sequences of a wheat-D. villosum translocation line, Pm97033-which contains chromosome 6V#4S carrying the PmV locus, under wheat PM pathogen induction. Twenty-five 6V#4S chromosome-specific markers were developed. Three of them were able to clearly distinguish chromosomes 6V#4S and 6V#2S by product size, four amplified the product specific for chromosome 6V#4S only, and the remaining 18 markers identified chromosome 6VS in wheat backgrounds. Two different D. villosum accessions, their derived translocation lines and wheat varieties carrying different chromosome 6VS were identified using these specific markers. The 25 newly developed markers together with the known PM resistance gene Stpk-V were used to construct comparative genomic maps with the homoeologous chromosome arms of wheat and barley. The colinearity of the identified gene sequences amplified by the 25 markers among wheat chromosomes 6A, 6B, and 6D and barley chromosome 6H was not very conserved and interrupted frequently by inversion and insertion. Our markers have potential in marker assisted selection for PM resistance breeding, and for locating other potential important genes and cloning the PmV gene on chromosome 6V#4S.

Connectivity in gene coexpression networks negatively correlates with rates of molecular evolution in flowering plants

Gene coexpression networks are a useful tool for summarizing transcriptomic data and providing insight into patterns of gene regulation in a variety of species. Though there has been considerable interest in studying the evolution of network topology across species, less attention has been paid to the relationship between network position and patterns of molecular evolution. Here, we generated coexpression networks from publicly available expression data for seven flowering plant taxa (Arabidopsis thaliana, Glycine max, Oryza sativa, Populus spp., Solanum lycopersicum, Vitis spp., and Zea mays) to investigate the relationship between network position and rates of molecular evolution. We found a significant negative correlation between network connectivity and rates of molecular evolution, with more highly connected (i.e., “hub”) genes having significantly lower nonsynonymous substitution rates and dN/dS ratios compared to less highly connected (i.e., “peripheral”) genes across the taxa surveyed. These findings suggest that more centrally located hub genes are, on average, subject to higher levels of evolutionary constraint than are genes located on the periphery of gene coexpression networks. The consistency of this result across disparate taxa suggests that it holds for flowering plants in general, as opposed to being a species-specific phenomenon.

High-Throughput Sequencing and Co-Expression Network Analysis of lncRNAs and mRNAs in Early Brain Injury Following Experimental Subarachnoid Haemorrhage

Subarachnoid haemorrhage (SAH) is a fatal neurovascular disease following cerebral aneurysm rupture with high morbidity and mortality rates. Long non-coding RNAs (lncRNAs) are a type of mammalian genome transcript, are abundantly expressed in the brain and are involved in many nervous system diseases. However, little is currently known regarding the influence of lncRNAs in early brain injury (EBI) after SAH. This study analysed the expression profiles of lncRNAs and mRNAs in SAH brain tissues of mice using high-throughput sequencing. The results showed a remarkable difference in lncRNA and mRNA transcripts between SAH and control brains. Approximately 617 lncRNA transcripts and 441 mRNA transcripts were aberrantly expressed at 24 hours after SAH. Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that the differentially expressed mRNAs were mostly involved in inflammation. Based on the lncRNA/mRNA co-expression network, knockdown of fantom3_F730004F19 reduced the mRNA and protein levels of CD14 and toll-like receptor 4 (TLR4) and attenuated inflammation in BV-2 microglia cells. These results indicate that lncRNA fantom3_F730004F19 may be associated with microglia induced inflammation via the TLR signaling pathway in EBI following SAH. LncRNA represent a potential therapeutic target for the prognosis, diagnosis, and treatment of SAH.

Constitutive heterologous overexpression of a TIR-NB-ARC-LRR gene encoding a putative disease resistance protein from wild Chinese Vitis pseudoreticulata in Arabidopsis and tobacco enhances resistance to phytopathogenic fungi and bacteria

Plants use resistance (R) proteins to detect pathogen effector proteins and activate their innate immune response against the pathogen. The majority of these proteins contain an NB-ARC (nucleotide-binding adaptor shared by APAF-1, R proteins, and CED-4) domain along with a leucine-rich repeat (LRR), and some also bear a toll interleukin 1 receptor (TIR) domain. In this study, we characterized a gene encoding a TIR-NB-ARC-LRR R protein (VpTNL1) (GenBank accession number KX649890) from wild Chinese grapevine Vitis pseudoreticulata accession "Baihe-35-1", which was identified previously from a transcriptomic analysis of leaves inoculated with powdery mildew (PM; Erysiphe necator (Schw.)). The VpTNL1 transcript was found to be highly induced in V. pseudoreticulata following inoculation with E. necator, as well as treatment with salicylic acid (SA). Sequence analysis demonstrated that the deduced amino acid sequence contained a TIR domain at the N-terminus, along with an NB-ARC and four LRRs domains within the C-terminus. Constitutive expression of VpTNL1 in Arabidopsis thaliana resulted in either a wild-type or dwarf phenotype. Intriguingly, the phenotypically normal transgenic lines displayed enhanced resistance to Arabidopsis PM, Golovinomyces cichoracearum, as well as to the virulent bacterial pathogen Pseudomonas syringae pv. tomato DC3000. Similarly, constitutive expression of VpTNL1 in Nicotiana tabacum was found to confer enhanced resistance to tobacco PM, Erysiphe cichoacearum DC. Subsequent isolation of the VpTNL1 promoter and deletion analysis indicated that TC-rich repeats and TCA elements likely play an important role in its response to E. necator and SA treatment, respectively. Taken together, these results indicate that VpTNL1 contributes to PM resistance in grapevine and provide an interesting gene target for the future amelioration of grape via breeding and/or biotechnology.