Molecular mapping of a novel lesion mimic gene (lm4) associated with enhanced resistance to stripe rust in bread wheat

Insights from homozygous signatures of cervus nippon revealed genetic architecture for components of fitness

This study investigates the genomic landscape of Sika deer populations, emphasizing the detection and characterization of runs of homozygosity (ROH) and their contribution towards components of fitness. Using 85,001 high-confidence SNPs, the investigation into ROH distribution unveiled nuanced patterns of autozygosity across individuals especially in 2 out of the 8 farms, exhibiting elevated ROH levels and mean genome coverage under ROH segments. The prevalence of shorter ROH segments (0.5-4 Mb) suggests historical relatedness and potential selective pressures within these populations. Intriguingly, despite observed variations in ROH profiles, the overall genomic inbreeding coefficient (FROH) remained relatively low across all farms, indicating a discernible degree of genetic exchange and effective mitigation of inbreeding within the studied Sika deer populations. Consensus ROH (cROH) were found to harbor genes for important functions viz., EGFLAM gene which is involved in the vision function of the eye, SKP2 gene which regulates cell cycle, CAPSL involved in adipogenesis, SPEF2 which is essential for sperm flagellar assembly, DCLK3 involved in the heat stress. This first ever study on ROH in Sika deer, to shed light on the adaptive role of genes in these homozygous regions. The insights garnered from this study have broader implications in the management of genetic diversity in this vulnerable species.

RiceMetaSys: Drought-miR, a one-stop solution for drought responsive miRNAs-mRNA module in rice

MicroRNAs are key players involved in stress responses in plants and reports are available on the role of miRNAs in drought stress response in rice. This work reports the development of a database, RiceMetaSys: Drought-miR, based on the meta-analysis of publicly available sRNA datasets. From 28 drought stress-specific sRNA datasets, we identified 216 drought-responsive miRNAs (DRMs). The major features of the database include genotype-, tissue- and miRNA ID-specific search options and comparison of genotypes to identify common miRNAs. Co-localization of the DRMs with the known quantitative trait loci (QTLs), i.e., meta-QTL regions governing drought tolerance in rice pertaining to different drought adaptive traits, narrowed down this to 37 promising DRMs. To identify the high confidence target genes of DRMs under drought stress, degradome datasets and web resource on drought-responsive genes (RiceMetaSys: DRG) were used. Out of the 216 unique DRMs, only 193 had targets with high stringent parameters. Out of the 1081 target genes identified by Degradome datasets, 730 showed differential expression under drought stress in at least one accession. To retrieve complete information on the target genes, the database has been linked with RiceMetaSys: DRG. Further, we updated the RiceMetaSys: DRGv1 developed earlier with the addition of DRGs identified from RNA-seq datasets from five rice genotypes. We also identified 759 putative novel miRNAs and their target genes employing stringent criteria. Novel miRNA search has all the search options of known miRNAs and additionally, it gives information on their in silico validation features. Simple sequence repeat markers for both the miRNAs and their target genes have also been designed and made available in the database. Network analysis of the target genes identified 60 hub genes which primarily act through abscisic acid pathway and jasmonic acid pathway. Co-localization of the hub genes with the meta-QTL regions governing drought tolerance narrowed down this to 16 most promising DRGs. Database URL: http://14.139.229.201/RiceMetaSys_miRNA Updated database of RiceMetaSys URL: http://14.139.229.201/RiceMetaSysA/Drought/.

Functional transcriptome analysis revealed upregulation of MAPK-SMAD signalling pathways in chronic heat stress in crossbred cattle

Animal geneticists and breeders have the impending challenge of enhancing the resilience of Indian livestock to heat stress through better selection strategies. Climate change's impact on livestock is more intense in tropical countries like India where dairy cattle crossbreeds are more sensitive to heat stress. The main reason for this study was to find the missing relative changes in transcript levels in thermo-neutral and heat stress conditions in crossbred cattle through whole-transcriptome analysis of RNA-Seq data. Differentially expressed genes (DEGs) identified based on the minimum log twofold change value and false discovery rate 0.05 revealed 468 up-regulated genes and 2273 down-regulated significant genes. Functional annotation and pathway analysis of these significant DEGs were compared based on Gene Ontology (Biological process), Kyoto Encyclopedia of Genes and Genome (KEGG), and Reactome pathways using g: Profiler, ShinyGO v0.76, and iDEP.951 web tools. On finding network visualization, the most over-represented and correlated pathways were neuronal and sensory organ development, calcium signalling pathway, Mitogen-activated protein kinase (MAPK) and Smad signalling pathway, Ras-proximate-1, or Ras-related protein 1 (Rap 1) signalling pathway, apoptosis, and oxidative stress. Similarly, down-regulated genes were most expressed in mRNA processing, immune system, B-cell receptor signalling pathway, Nucleotide oligomerization domain (NOD)-like receptors (NLRs) signalling pathway and nonsense-mediated decay (NMD) pathway. The heat stress-responsive genes identified in this study will facilitate our understanding of the molecular basis for climate resilience and heat tolerance in Indian dairy crossbreeds.

Integration of transcriptomics, metabolomics, and hormone analysis revealed the formation of lesion spots inhibited by GA and CTK was related to cell death and disease resistance in bread wheat (Triticum aestivum L.)

BACKGROUND

Wheat is one of the important grain crops in the world. The formation of lesion spots related to cell death is involved in disease resistance, whereas the regulatory pathway of lesion spot production and resistance mechanism to pathogens in wheat is largely unknown.

RESULTS

In this study, a pair of NILs (NIL-Lm5W and NIL-Lm5M) was constructed from the BC1F4 population by the wheat lesion mimic mutant MC21 and its wild genotype Chuannong 16. The formation of lesion spots in NIL-Lm5M significantly increased its resistance to stripe rust, and NIL-Lm5M showed superiour agronomic traits than NIL-Lm5W under stripe rust infection.Whereafter, the NILs were subjected to transcriptomic (stage N: no spots; stage S, only a few spots; and stage M, numerous spots), metabolomic (stage N and S), and hormone analysis (stage S), with samples taken from normal plants in the field. Transcriptomic analysis showed that the differentially expressed genes were enriched in plant-pathogen interaction, and defense-related genes were significantly upregulated following the formation of lesion spots. Metabolomic analysis showed that the differentially accumulated metabolites were enriched in energy metabolism, including amino acid metabolism, carbohydrate metabolism, and lipid metabolism. Correlation network diagrams of transcriptomic and metabolomic showed that they were both enriched in energy metabolism. Additionally, the contents of gibberellin A7, cis-Zeatin, and abscisic acid were decreased in leaves upon lesion spot formation, whereas the lesion spots in NIL-Lm5M leaves were restrained by spaying GA and cytokinin (CTK, trans-zeatin) in the field.

CONCLUSION

The formation of lesion spots can result in cell death and enhance strip rust resistance by protein degradation pathway and defense-related genes overexpression in wheat. Besides, the formation of lesion spots was significantly affected by GA and CTK. Altogether, these results may contribute to the understanding of lesion spot formation in wheat and laid a foundation for regulating the resistance mechanism to stripe rust.

Recent advancements in the physiological, genetic, and genomic research on Rhododendrons for trait improvement

High species diversity, hybridization potential, broad geographical dispersal range and ornamental characteristics (i.e., attractive size, shape, structure, flowers, and evergreen) have fetched a good international market for Rhododendron. However, most species are restricted to specific geographic areas due to their habitat specificity in acidic soil and cold climates, resulting many species being classified under threat categories of the IUCN. In this review, advances in research on Rhododendron for improvement to floral display quality and stress resistance have been described. The low genetic barrier among species has created opportunities for extensive hybridization and ploidy alteration for introducing quality and adaptive traits during the development of new varieties. Recent technological advances have supported investigations into the mechanism of flower development, as well as cold tolerance and pathogen resistance mechanisms in the Rhododendron. However, most of the species have limited adaptability to drought, line-tolerance, pathogen resistance, and high-temperature conditions and this resistance ability present in few species largely remains unexplored. Additionally, the available genetic diversity and genomic information on species, and possibilities for their application in molecular breeding have been summarized. Overall, genomic resource data are scarce in the majority of the members of this genus. Finally, various research gaps such as genetic mapping of quality traits, understanding the molecular mechanism of quality-related traits and genomic assortment in Rhododendron members have been discussed in the future perspective section.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04006-6.

A review of the role of transcription factors in regulating adipogenesis and lipogenesis in beef cattle

In the past few decades, genomic selection and other refined strategies have been used to increase the growth rate and lean meat production of beef cattle. Nevertheless, the fast growth rates of cattle breeds are often accompanied by a reduction in intramuscular fat (IMF) deposition, impairing meat quality. Transcription factors play vital roles in regulating adipogenesis and lipogenesis in beef cattle. Meanwhile, understanding the role of transcription factors in regulating adipogenesis and lipogenesis in beef cattle has gained significant attention to increase IMF deposition and meat quality. Therefore, the aim of this paper was to provide a comprehensive summary and valuable insight into the complex role of transcription factors in adipogenesis and lipogenesis in beef cattle. This review summarizes the contemporary studies in transcription factors in adipogenesis and lipogenesis, genome-wide analysis of transcription factors, epigenetic regulation of transcription factors, nutritional regulation of transcription factors, metabolic signalling pathways, functional genomics methods, transcriptomic profiling of adipose tissues, transcription factors and meat quality and comparative genomics with other livestock species. In conclusion, transcription factors play a crucial role in promoting adipocyte development and fatty acid biosynthesis in beef cattle. They control adipose tissue formation and metabolism, thereby improving meat quality and maintaining metabolic balance. Understanding the processes by which these transcription factors regulate adipose tissue deposition and lipid metabolism will simplify the development of marbling or IMF composition in beef cattle.

Molecular characterization and evolutionary relationships of DOFs in four cherry species and functional analysis in sweet cherry

The DOF (DNA binding with one finger) has multiple functions in plants. However, it has received little attention in the research field of cherries. In this study, the evolutionary relationship and molecular characterization of DOF in four cherry species were analyzed, revealing its expression pattern in sweet cherry. There are 23 members in Prunus avium cv. 'Tieton', 88 in Prunus cerasus, 53 in Cerasus × yedoensis, and 27 in Cerasus serrulata. Most of these genes are intron-less or non-intron, with a conserved C2-C2 domain. Due to heterozygosity and chromosomal ploidy, whole-genome duplication (WGD) events occur to varying degrees, and DOF genes are contracted during evolution. Furthermore, these genes are affected by purifying selection pressure. Under low-temperature treatment, the expression of PavDOF2 and PavDOF18 were significantly up-regulated, while that of PavDOF16 is significantly down-regulated. The expression of PavDOF9, PavDOF12, PavDOF14, PavDOF16, PavDOF17, PavDOF18, and PavDOF19 exhibits an increasing trend during flower development and varies during sweet cherry fruit development. PavDOF1, PavDOF8, PavDOF9, and PavDOF15 are localized in the nucleus but is not transcriptionally active. The findings systemically demonstrate the molecular characteristics of DOF in different cherry varieties, providing a basis for further research on the functions of these genes.

Whole genome discovery of regulatory genes responsible for the response of chicken to heat stress

Long noncoding RNAs (lncRNAs) are functional bridges connecting the genome with phenotypes by interacting with DNA, mRNA, and proteins. Using publically available acute heat stress (AHS)-related RNA-seq data, we discovered novel lncRNAs and tested their association with AHS along with ~ 8800 known lncRNAs and ~ 28,000 mRNA transcripts. Our pipeline discovered a total of 145 potentially novel-lncRNAs. One of them (Fishcomb_p-value = 0.06) along with another novel transcript (annotated as protein-coding; Fishcomb_p-value = 0.03) were identified as significantly associated with AHS. We found five known-lncRNAs and 134 mRNAs transcripts that were significantly associated with AHS. Four novel lncRNAs interact cis-regulated with 12 mRNA transcripts and are targeted by 11 miRNAs. Also six meta-lncRNAs associate with 134 meta-mRNAs through trans-acting co-expression, each targeted by 15 and 216 miRNAs, respectively. Three of the known-lncRNAs significantly co-expressed with almost 97 of the significant mRNAs (Pearson correlation p-value < 0.05). We report the mentioned three known-lncRNAs (ENSGALT00000099876, ENSGALT00000107573, and ENSGALT00000106323) as the most, significantly regulatory elements of AHS in chicken. It can be concluded that in order to alleviate the adverse effects of AHS on chicken, the manipulation of the three regulatory lncRNAs could lead to a more desirable result than the manipulation of the most significant mRNAs.

Host Response of Arabidopsis thaliana Interaction with Fungal Endophytes Involves microRNAs

Plant and fungus interaction is a complex process involving many molecular factors determining the nature of relationship. The enigmatic methodology by which fungal endophytes are able to colonise a plant harmoniously is still inexplicable. Small RNAs have been identified as major regulatory elements under various biotic interactions. However, their role in endophytic plant-fungal interactions remain to be elucidated. Therefore, transcript expression data available on Gene Expression Omnibus for Arabidopsis thaliana was utilised for miRNAs identification under endophytism. The analysis predicted 15 miRNAs with differential expression of which the ath-miRNA398b modulation was significant. Application of psRNAtarget, C-mii, pmiREN, and TarDB provided a pool of 357 target genes for these miRNAs. Protein-protein interaction analysis identified major hub proteins, including BTB/POZ domain-containing protein, beta-Xylosidase-2 (AtBXL2), and Copper/Zinc Superoxide Dismutase-2 (AtSOD2). The quantitative real-time PCR validated the computational prediction and expression for selected target genes AtSOD2, AtBXL2, and AtRCA along with ath-miRNA398b under endophytism. Overall, results indicate that miRNAs have a significant role in regulating Arabidopsis thaliana-endophytic fungal interaction.

Transcriptome analysis reveals the underlying mechanism for over-accumulation of alkaline protease in Bacillus licheniformis

AIMS

Bacillus licheniformis AQ is an industrial strain with high production of alkaline protease (AprE), which has great industrial application value. However, how to regulate the production of AprE in the process of industrial fermentation is still not completely clear. Therefore, it is important to understand the metabolic process of AprE production in the industrial fermentation medium.

METHODS AND RESULTS

In this study, transcriptome sequencing of the whole fermentation course was performed to explore the synthesis and regulation mechanism of AprE in B. licheniformis AQ. During the fermentation process, the AprE got continuously accumulated, reaching a peak of 42 020 U/mL at the fermentation endpoint (48 h). Meanwhile, the highly expressed genes were observed. Compared with the fermentation endpoint, there were 61 genes in the intersection of differentially expressed genes, functioning as catabolic processes, peptidases and inhibitors, chaperones, and folding catalysts. Furthermore, the protein-protein interactions network of AprE was constructed.

CONCLUSION

This study provides important transcriptome information for B. licheniformis AQ and potential molecular targets for further improving the production of AprE.

Revelation of genetic diversity and genomic footprints of adaptation in Indian pig breeds

In the present study, the genetic diversity measures among four Indian domestic breeds of pig namely Agonda Goan, Ghurrah, Ghungroo, and Nicobari, of different agro-climatic regions of country were explored and compared with European commercial breeds, European wild boar and Chinese domestic breeds. The double digest restriction site-associated DNA sequencing (ddRADseq) data of Indian pigs (102) and Landrace (10 animals) were generated and whole genome sequencing data of exotic pigs (60 animals) from public data repository were used in the study. The principal component analysis (PCA), admixture analysis and phylogenetic analysis revealed that Indian breeds were closer in ancestry to Chinese breeds than European breeds. European breeds exhibited highest genetic diversity measures among all the considered breeds. Among Indian breeds, Agonda Goan and Ghurrah were found to be more genetically diverse than Nicobari and Ghungroo. The selection signature regions in Indian pigs were explored using iHS and XP-EHH, and during iHS analysis, it was observed that genes related to growth, reproduction, health, meat quality, sensory perception and behavior were found to be under selection pressure in Indian pig breeds. Strong selection signatures were recorded in 24.25-25.25 Mb region of SSC18, 123.25-124 Mb region of SSC15 and 118.75-119.5 Mb region of SSC2 in most of the Indian breeds upon pairwise comparison with European commercial breeds using XP-EHH. These regions were harboring some important genes such as EPHA4 for thermotolerance, TAS2R16, FEZF1, CADPS2 and PTPRZ1 for adaptability to scavenging system of rearing, TRIM36 and PGGT1B for disease resistance and CCDC112, PIAS1, FEM1B and ITGA11 for reproduction.

A transcriptomic analysis of skeletal muscle tissues reveals promising candidate genes and pathways accountable for different daily weight gain in Hanwoo cattle

Cattle traits like average daily weight gain (ADG) greatly impact profitability. Selecting based on ADG considering genetic variability can lead to economic and genetic advancements in cattle breeding. This study aimed to unravel genetic influences on ADG variation in Hanwoo cattle at the skeletal muscle transcriptomic level. RNA sequencing was conducted on longissimus dorsi (LD), semimembranosus (SB), and psoas major (PM) muscles of 14 steers assigned to same feed, grouped by low (≤ 0.71 kg) and high (≥ 0.77 kg) ADG. At P ≤ 0.05 and log2fold > 1.5, the distinct pattern of gene expression was identified with 184, 172, and 210 differentially expressed genes in LD, SB, and PM muscles, respectively. Tissue-specific responses to ADG variation were evident, with myogenesis and differentiation associated JAK-STAT signaling pathway and prolactin signaling pathways enriched in LD and SB muscles, while adipogenesis-related PPAR signaling pathways were enriched in PM muscle. Key hub genes (AXIN2, CDKN1A, MYC, PTGS2, FZD5, SPP1) were upregulated and functionally significant in muscle growth and differentiation. Notably, DPP6, CDKN1A, and FZD5 emerged as possible candidate genes linked to ADG variation. These findings enhance our understanding of genetic factors behind ADG variation in Hanwoo cattle, illuminating skeletal muscle mechanisms influencing ADG.

Systematic evaluation of multiple NGS platforms for structural variants detection

Structural variations (SV) are critical genome changes affecting human diseases. Although many hybridization-based methods exist, evaluating SVs through next-generation sequencing (NGS) data is still necessary for broader research exploration. Here, we comprehensively compared the performance of 16 SV callers and multiple NGS platforms using NA12878 whole genome sequencing (WGS) datasets. The results indicated that several SV callers performed well relatively, such as Manta, GRIDSS, LUMPY, TARDIS, FermiKit, and Wham. Meanwhile, all NGS platforms have a similar performance using a single software. Additionally, we found that the source of undetected SVs was mostly from long reads datasets, therefore, the more appropriate strategy for accurate SV detection will be an integration of long and shorter reads in the future. At present, in the period of NGS as a mainstream method in bioinformatics, our study would provide helpful and comprehensive guidelines for specific categories of SV research.

Non-coding RNAs (ncRNAs) in plant: Master regulators for adapting to extreme temperature conditions

Unusual daily temperature fluctuations caused by climate change and climate variability adversely impact agricultural crop production. Since plants are immobile and constantly receive external environmental signals, such as extreme high (heat) and low (cold) temperatures, they have developed complex molecular regulatory mechanisms to cope with stressful situations to sustain their natural growth and development. Among these mechanisms, non-coding RNAs (ncRNAs), particularly microRNAs (miRNAs), small-interfering RNAs (siRNAs), and long-non-coding RNAs (lncRNAs), play a significant role in enhancing heat and cold stress tolerance. This review explores the pivotal findings related to miRNAs, siRNAs, and lncRNAs, elucidating how they functionally regulate plant adaptation to extreme temperatures. In addition, this review addresses the challenges associated with uncovering these non-coding RNAs and understanding their roles in orchestrating heat and cold tolerance in plants.

Maize ZmBES1/BZR1-1 transcription factor negatively regulates drought tolerance

Drought stress is a common abiotic factor and restricts plant growth and development. Exploring maize stress-related genes and their regulatory mechanisms is crucial for ensuring agricultural productivity and food security. The BRI1-EMS1 suppressor (BES1)/brassinazole-resistant 1 (BZR1) transcription factors (TFs) play important roles in plant growth, development, and stress response. However, maize ZmBES1/BZR1s are rarely reported. In the present study, the ZmBES1/BZR1-1 gene was cloned from maize B73 and functionally characterized in transgenic Arabidopsis and rice in drought stress response. The ZmBES1/BZR1-1 protein possessed a conserved bHLH domain characterized by BES1/BZR1 TFs, localized in the nucleus, and showed transcription activation activity. The expression of ZmBES1/BZR1-1 exhibited no tissue specificity but drought-inhibitory expression in maize. Under drought stress, overexpression of ZmBES1/BZR1-1 resulted in the enhancement of drought sensitivity of transgenic Arabidopsis and rice with a lower survival rate, reactive oxygen species (ROS) level and relative water content (RWC) and a higher stomatal aperture and relative electrolyte leakage (REL). The RNA-seq results showed that 56 differentially expressed genes (DEGs) were regulated by ZmBES1/BZR1-1 by binding to E-box elements in their promoters. The GO analysis showed that the DEGs were significantly annotated with response to oxidative stress and oxygen level. The study suggests that the ZmBES1/BZR1-1 gene negatively regulates drought stress, which provides insights into further underlying molecular mechanisms in the drought stress response mediated by BZR1/BES1s.

Role of miRNAs in sucrose stress response, reactive oxygen species, and anthocyanin biosynthesis in Arabidopsis thaliana

In plants, sucrose is the main transported disaccharide that is the primary product of photosynthesis and controls a multitude of aspects of the plant life cycle including structure, growth, development, and stress response. Sucrose is a signaling molecule facilitating various stress adaptations by crosstalk with other hormones, but the molecular mechanisms are not well understood. Accumulation of high sucrose concentrations is a hallmark of many abiotic and biotic stresses, resulting in the accumulation of reactive oxygen species and secondary metabolite anthocyanins that have antioxidant properties. Previous studies have shown that several MYeloBlastosis family/MYB transcription factors are positive and negative regulators of sucrose-induced anthocyanin accumulation and subject to microRNA (miRNA)-mediated post-transcriptional silencing, consistent with the notion that miRNAs may be "nodes" in crosstalk signaling by virtue of their sequence-guided targeting of different homologous family members. In this study, we endeavored to uncover by deep sequencing small RNA and mRNA transcriptomes the effects of exogenous high sucrose stress on miRNA abundances and their validated target transcripts in Arabidopsis. We focused on genotype-by-treatment effects of high sucrose stress in Production of Anthocyanin Pigment 1-Dominant/pap1-D, an activation-tagged dominant allele of MYB75 transcription factor, a positive effector of secondary metabolite anthocyanin pathway. In the process, we discovered links to reactive oxygen species signaling through miR158/161/173-targeted Pentatrico Peptide Repeat genes and two novel non-canonical targets of high sucrose-induced miR408 and miR398b*(star), relevant to carbon metabolic fluxes: Flavonoid 3'-Hydroxlase (F3'H), an important enzyme in determining the B-ring hydroxylation pattern of flavonoids, and ORANGE a post-translational regulator of Phytoene Synthase expression, respectively. Taken together, our results contribute to understanding the molecular mechanisms of carbon flux shifts from primary to secondary metabolites in response to high sugar stress.

Species-specific microRNA discovery and target prediction in the soybean cyst nematode

The soybean cyst nematode (SCN) is a devastating pathogen for economic and food security considerations. Although the SCN genome has recently been sequenced, the presence of any miRNA has not been systematically explored and reported. This paper describes the development of a species-specific SCN miRNA discovery pipeline and its application to the SCN genome. Experiments on well-documented model nematodes (Caenorhabditis elegans and Pristionchus pacificus) are used to tune the pipeline's hyperparameters and confirm its recall and precision. Application to the SCN genome identifies 3342 high-confidence putative SCN miRNA. Prediction specificity within SCN is confirmed by applying the pipeline to RNA hairpins from known exonic regions of the SCN genome (i.e., sequences known to not be miRNA). Prediction recall is confirmed by building a positive control set of SCN miRNA, based on a limited deep sequencing experiment. Interestingly, a number of novel miRNA are predicted to be encoded within the intronic regions of effector genes, known to be involved in SCN parasitism, suggesting that these miRNA may also be involved in the infection process or virulence. Beyond miRNA discovery, gene targets within SCN are predicted for all high-confidence novel miRNA using a miRNA:mRNA target prediction system. Lastly, cross-kingdom miRNA targeting is investigated, where putative soybean mRNA targets are identified for novel SCN miRNA. All predicted miRNA and gene targets are made available in appendix and through a Borealis DataVerse open repository ( https://borealisdata.ca/dataset.xhtml?persistentId=doi:10.5683/SP3/30DEXA ).

Global transcriptome analysis reveals fungal disease responsive core gene regulatory landscape in tea

Fungal infections are the inevitable limiting factor for productivity of tea. Transcriptome reprogramming recruits multiple regulatory pathways during pathogen infection. A comprehensive meta-analysis was performed utilizing previously reported, well-replicated transcriptomic datasets from seven fungal diseases of tea. The study identified a cumulative set of 18,517 differentially expressed genes (DEGs) in tea, implicated in several functional clusters, including the MAPK signaling pathway, transcriptional regulation, and the biosynthesis of phenylpropanoids. Gene set enrichment analyses under each pathogen stress elucidated that DEGs were involved in ethylene metabolism, secondary metabolism, receptor kinase activity, and various reactive oxygen species detoxification enzyme activities. Expressional fold change of combined datasets highlighting 2258 meta-DEGs shared a common transcriptomic response upon fungal stress in tea. Pervasive duplication events caused biotic stress-responsive core DEGs to appear in multiple copies throughout the tea genome. The co-expression network of meta-DEGs in multiple modules demonstrated the coordination of appropriate pathways, most of which involved cell wall organization. The functional coordination was controlled by a number of hub genes and miRNAs, leading to pathogenic resistance or susceptibility. This first-of-its-kind meta-analysis of host-pathogen interaction generated consensus candidate loci as molecular signatures, which can be associated with future resistance breeding programs in tea.

Disruption of the Novel Small Protein RBR7 Leads to Enhanced Plant Resistance to Blast Disease

Plant disease is a threat to global food security. Breeding crops carrying broad-spectrum resistance loci is an effective way to control infectious disease. Disease-resistant mutants are valuable resources for deciphering the underlying mechanisms of plant immunity and could provide genetic loci to generate disease-resistant crops. Here, we identified a rice mutant, rbr7 (rice blast resistance 7), that confers resistance against different strains of Magnaporthe oryzae. Disease-mimicking necrotic lesions started to appear on the leaves of rbr7 four weeks after sowing. Histochemical analysis revealed reactive oxygen species accumulation and cell death accompanied by spontaneous lesion formation in rbr7. Map-based cloning and bulk segregation analysis showed a 2855 bp fragment deletion on chromosome 5, leading to the disruption of the LOC_Os05g28480-coding protein. Transgenic rbr7 complementation plants showed compromised resistance to rice blast, indicating that LOC_Os05g28480, or Rbr7, regulates the rice immune response. Rbr7 encodes a small protein of unknown function with 85 amino acids. Transcriptomic analysis revealed that disruption of RBR7 led to the upregulation of genes responding to salicylic acid, systemic acquired resistance and pathogenesis-related genes. Taken together, our findings reveal insights into a novel small protein involved in regulating plant resistance to rice blast and provide a potential target for crop breeding.

Natural products of pentacyclic triterpenoids: from discovery to heterologous biosynthesis

Covering: up to 2022Pentacyclic triterpenoids are important natural bioactive substances that are widely present in plants and fungi. They have significant medicinal efficacy, play an important role in reducing blood glucose and protecting the liver, and have anti-inflammatory, anti-oxidation, anti-fatigue, anti-viral, and anti-cancer activities. Pentacyclic triterpenoids are derived from the isoprenoid biosynthetic pathway, which generates common precursors of triterpenes and steroids, followed by cyclization with oxidosqualene cyclases (OSCs) and decoration via cytochrome P450 monooxygenases (CYP450s) and glycosyltransferases (GTs). Many biosynthetic pathways of triterpenoid saponins have been elucidated by studying their metabolic regulation network through the use of multiomics and identifying their functional genes. Unfortunately, natural resources of pentacyclic triterpenoids are limited due to their low content in plant tissues and the long growth cycle of plants. Based on the understanding of their biosynthetic pathway and transcriptional regulation, plant bioreactors and microbial cell factories are emerging as alternative means for the synthesis of desired triterpenoid saponins. The rapid development of synthetic biology, metabolic engineering, and fermentation technology has broadened channels for the accumulation of pentacyclic triterpenoid saponins. In this review, we summarize the classification, distribution, structural characteristics, and bioactivity of pentacyclic triterpenoids. We further discuss the biosynthetic pathways of pentacyclic triterpenoids and involved transcriptional regulation. Moreover, the recent progress and characteristics of heterologous biosynthesis in plants and microbial cell factories are discussed comparatively. Finally, we propose potential strategies to improve the accumulation of triterpenoid saponins, thereby providing a guide for their future biomanufacturing.

In ovo probiotic supplementation promotes muscle growth and development in broiler embryos

In chickens, muscle development during embryonic growth is predominantly by myofiber hyperplasia. Following hatch, muscle growth primarily occurs via hypertrophy of the existing myofibers. Since myofiber number is set at hatch, production of more muscle fibers during embryonic growth would provide a greater myofiber number at hatch and potential for posthatch muscle growth by hypertrophy. Therefore, to improve performance in broilers, this study investigated the effect of in ovo spray application of probiotics on overall morphometry and muscle development in broiler embryos. For the study, fertile Ross 308 eggs were sprayed with different probiotics; Lactobacillus paracasei DUP 13076 (LP) and L. rhamnosus NRRL B 442 (LR) prior to and during incubation. The embryos were sacrificed on d 7, 10, 14, and 18 for embryo morphometry and pectoralis major muscle (PMM) sampling. Muscle sections were stained and imaged to quantify muscle fiber density (MFD), myofiber cross-sectional area (CSA), and nuclei density. Additionally, gene expression assays were performed to elucidate the effect of probiotics on myogenic genes. In ovo probiotic supplementation was found to significantly improve embryo weight, breast weight, and leg weight (P < 0.05). Further, histological analysis of PMM revealed a significant increase in MFD and nuclei number in the probiotic-treated embryos when compared to the control (P < 0.05). In 18-day-old broiler embryos, myofibers in the treatment group had a significantly smaller CSA (LP: 95.27 ± 3.28 μm², LR: 178.84 ± 15.1 μm²) when compared to the control (211.41 ± 15.67 μm²). This decrease in CSA was found to be associated with a concomitant increase in MFD (fibers/mm²) in the LP (13,647 ± 482.15) and LR (13,957 ± 463.13) group when compared to the control (7,680 ± 406.78). Additionally, this increase in myofibrillar hyperplasia in the treatment groups was associated with upregulation in the expression of key genes regulating muscle growth including MYF5, MYOD, MYOG, and IGF-1. In summary, in ovo spray application of probiotics promoted overall embryo growth and muscle development in broilers.

The diverging epigenomic landscapes of honeybee queens and workers revealed by multiomic sequencing

The role of the epigenome in phenotypic plasticity is unclear presently. Here we used a multiomics approach to explore the nature of the epigenome in developing honey bee (Apis mellifera) workers and queens. Our data clearly showed distinct queen and worker epigenomic landscapes during the developmental process. Differences in gene expression between workers and queens become more extensive and more layered during the process of development. Genes known to be important for caste differentiation were more likely to be regulated by multiple epigenomic systems than other differentially expressed genes. We confirmed the importance of two candidate genes for caste differentiation by using RNAi to manipulate the expression of two genes that differed in expression between workers and queens were regulated by multiple epigenomic systems. For both genes the RNAi manipulation resulted in a decrease in weight and fewer ovarioles of newly emerged queens compared to controls. Our data show that the distinct epigenomic landscapes of worker and queen bees differentiate during the course of larval development.

Environment-sensitive genic male sterility in rice and other plants

Environment-sensitive genic male sterility is a type of male sterility that is affected by both genetic and environmental factors. Environment-sensitive genic male sterile lines are not only used in two-line hybrid breeding but are also good materials for studying plant-environment interactions. In this study we review the research progress on environment-sensitive genic male sterility in rice from the perspectives of epigenetic, transcriptional, posttranscriptional, posttranslational and metabolic mechanisms as well as signal transduction processes. While significant progress has been made in the genetics, gene cloning and understanding of the molecular mechanisms of environment-sensitive genic male sterility in recent years, the relevant regulatory network is still poorly understood in rice. We therefore also review studies of environment-sensitive genic male sterility in Arabidopsis and other crops, hoping to promote research in this field in rice. Finally, we analyse the challenges posed by environment-sensitive genic male sterility and provide corresponding suggestions. This review will contribute towards an understanding the molecular genetics of environment-sensitive genic male sterility and its application in two-line hybrid breeding in rice and other species.

Integrated single-molecule real-time sequencing and RNA sequencing reveal the molecular mechanisms of salt tolerance in a novel synthesized polyploid genetic bridge between maize and its wild relatives

BACKGROUND

Tripsacum dactyloides (2n = 4x = 72) and Zea perennis (2n = 4x = 40) are tertiary gene pools of Zea mays L. and exhibit many abiotic adaptations absent in modern maize, especially salt tolerance. A previously reported allopolyploid (hereafter referred to as MTP, 2n = 74) synthesized using Zea mays, Tripsacum dactyloides, and Zea perennis has even stronger salt tolerance than Z. perennis and T. dactyloides. This allopolyploid will be a powerful genetic bridge for the genetic improvement of maize. However, the molecular mechanisms underlying its salt tolerance, as well as the key genes involved in regulating its salt tolerance, remain unclear.

RESULTS

Single-molecule real-time sequencing and RNA sequencing were used to identify the genes involved in salt tolerance and reveal the underlying molecular mechanisms. Based on the SMRT-seq results, we obtained 227,375 reference unigenes with an average length of 2300 bp; most of the unigenes were annotated to Z. mays sequences (76.5%) in the NR database. Moreover, a total of 484 and 1053 differentially expressed genes (DEGs) were identified in the leaves and roots, respectively. Functional enrichment analysis of DEGs revealed that multiple pathways responded to salt stress, including "Flavonoid biosynthesis," "Oxidoreductase activity," and "Plant hormone signal transduction" in the leaves and roots, and "Iron ion binding," "Acetyl-CoA carboxylase activity," and "Serine-type carboxypeptidase activity" in the roots. Transcription factors, such as those in the WRKY, B3-ARF, and bHLH families, and cytokinin negatively regulators negatively regulated the salt stress response. According to the results of the short time series-expression miner analysis, proteins involved in "Spliceosome" and "MAPK signal pathway" dynamically responded to salt stress as salinity changed. Protein-protein interaction analysis revealed that heat shock proteins play a role in the large interaction network regulating salt tolerance.

CONCLUSIONS

Our results reveal the molecular mechanism underlying the regulation of MTP in the response to salt stress and abundant salt-tolerance-related unigenes. These findings will aid the retrieval of lost alleles in modern maize and provide a new approach for using T. dactyloides and Z. perennis to improve maize.

Systematic and benchmarking studies of pipelines for mammal WGBS data in the novel NGS platform

BACKGROUND

Whole genome bisulfite sequencing (WGBS), possesses the aptitude to dissect methylation status at the nucleotide-level resolution of 5-methylcytosine (5-mC) on a genome-wide scale. It is a powerful technique for epigenome in various cell types, and tissues. As a recently established next-generation sequencing (NGS) platform, GenoLab M is a promising alternative platform. However, its comprehensive evaluation for WGBS has not been reported. We sequenced two bisulfite-converted mammal DNA in this research using our GenoLab M and NovaSeq 6000, respectively. Then, we systematically compared those data via four widely used WGBS tools (BSMAP, Bismark, BatMeth2, BS-Seeker2) and a new bisulfite-seq tool (BSBolt). We interrogated their computational time, genome depth and coverage, and evaluated their percentage of methylated Cs.

RESULT

Here, benchmarking a combination of pre- and post-processing methods, we found that trimming improved the performance of mapping efficiency in eight datasets. The data from two platforms uncovered ~ 80% of CpG sites genome-wide in the human cell line. Those data sequenced by GenoLab M achieved a far lower proportion of duplicates (~ 5.5%). Among pipelines, BSMAP provided an intriguing representation of 5-mC distribution at CpG sites with 5-mC levels > ~ 78% in datasets from human cell lines, especially in the GenoLab M. BSMAP performed more advantages in running time, uniquely mapped reads percentages, genomic coverage, and quantitative accuracy. Finally, compared with the previous methylation pattern of human cell line and mouse tissue, we confirmed that the data from GenoLab M performed similar consistency and accuracy in methylation levels of CpG sites with that from NovaSeq 6000.

CONCLUSION

Together we confirmed that GenoLab M was a qualified NGS platform for WGBS with high performance. Our results showed that BSMAP was the suitable pipeline that allowed for WGBS studies on the GenoLab M platform.

Global analysis of the association between pig muscle fatty acid composition and gene expression using RNA-Seq

Fatty acids (FAs) play an essential role as mediators of cell signaling and signal transduction, affecting metabolic homeostasis and determining meat quality in pigs. However, FAs are transformed by the action of several genes, such as those encoding desaturases and elongases of FAs in lipogenic tissues. The aim of the current work was to identify candidate genes, biological processes, and pathways involved in the modulation of intramuscular FA profile from longissimus dorsi muscle. FA profile by gas chromatography of methyl esters and gene expression by RNA-Seq were determined in 129 Iberian × Duroc backcrossed pigs. An association analysis between the muscle transcriptome and its FA profile was performed, followed by a concordance and functional analysis. Overall, a list of well-known (e.g., PLIN1, LEP, ELOVL6, SC5D, NCOA2, ACSL1, MDH1, LPL, LGALS12, TFRC, GOT1, and FBP1) and novel (e.g., TRARG1, TANK, ENSSSCG00000011196, and ENSSSCG00000038429) candidate genes was identified, either in association with specific or several FA traits. Likewise, several of these genes belong to biological processes and pathways linked to energy, lipid, and carbohydrate metabolism, which seem determinants in the modulation of FA compositions. This study can contribute to elucidate the complex relationship between gene expression and FA profile in pig muscle.

P-TarPmiR accurately predicts plant-specific miRNA targets

microRNAs (miRNAs) are small non-coding ribonucleic acids that post-transcriptionally regulate gene expression through the targeting of messenger RNA (mRNAs). Most miRNA target predictors have focused on animal species and prediction performance drops substantially when applied to plant species. Several rule-based miRNA target predictors have been developed in plant species, but they often fail to discover new miRNA targets with non-canonical miRNA-mRNA binding. Here, the recently published TarDB database of plant miRNA-mRNA data is leveraged to retrain the TarPmiR miRNA target predictor for application on plant species. Rigorous experiment design across four plant test species demonstrates that animal-trained predictors fail to sustain performance on plant species, and that the use of plant-specific training data improves accuracy depending on the quantity of plant training data used. Surprisingly, our results indicate that the complete exclusion of animal training data leads to the most accurate plant-specific miRNA target predictor indicating that animal-based data may detract from miRNA target prediction in plants. Our final plant-specific miRNA prediction method, dubbed P-TarPmiR, is freely available for use at http://ptarpmir.cu-bic.ca . The final P-TarPmiR method is used to predict targets for all miRNA within the soybean genome. Those ranked predictions, together with GO term enrichment, are shared with the research community.

Genome-wide identification, characterization and gene expression of BES1 transcription factor family in grapevine (Vitis vinifera L.)

BES1, as the most important transcription factor responsible for brassinolide (BR) signaling, has been confirmed to play a significant role in regulating plant growth and the improvement of stress resistance. The transcriptional regulatory mechanism of BES1 has been well elucidated in several plants, such as Arabidopsis thaliana (A. thaliana), Triticum aestivum L. (T. aestivum), and Oryza sativa L. (O. sativa). Nevertheless, the genome-wide analysis of the BES1 family in Vitis vinifera L. (V. vinifera). has not been comprehensively carried out. Thus, we have conducted a detailed analysis and identification of the BES1 transcription factors family in V. vinifera; a total of eight VvBES1 genes was predicted, and the phylogenetic relationships, gene structures, and Cis-acting element in their promoters were also analyzed. BES1 genes have been divided into three groups (I, II and III) based on phylogenetic relationship analysis, and most of VvBES1 genes were in group III. Also, we found that VvBES1 genes was located at seven of the total nineteen chromosomes, whereas VvBES1-2 (Vitvi04g01234) and VvBES1-5 (Vitvi18g00924) had a collinearity relationship, and their three copies are well preserved. In addition, the intron-exon model of VvBES1 genes were mostly conserved, and there existed several Cis-acting elements related to stress resistance responsive and phytohormones responsive in BES1s genes promoter. Moreover, the BES1 expressions were different in different V. vinifera organs, and BES1 expressions were different in different V. vinifera varieties under saline-alkali stress and heat stress, the expression of VvBES1 also changed with the prolongation of saline-alkali stress treatment time. The above findings could not only lay a primary foundation for the further validation of VvBES1 function, but could also provide a reference for molecular breeding in V. vinifera.

Accuracy benchmark of the GeneMind GenoLab M sequencing platform for WGS and WES analysis

Background

GenoLab M is a recently developed next-generation sequencing (NGS) platform from GeneMind Biosciences. To establish the performance of GenoLab M, we present the first report to benchmark and compare the WGS and WES sequencing data of the GenoLab M sequencer to NovaSeq 6000 and NextSeq 550 platform in various types of analysis. For WGS, thirty-fold sequencing from Illumina NovaSeq platform and processed by GATK pipeline is currently considered as the golden standard. Thus this dataset is generated as a benchmark reference in this study.

Results

GenoLab M showed an average of 94.62% of Q20 percentage for base quality, while the NovaSeq was slightly higher at 96.97%. However, GenoLab M outperformed NovaSeq or NextSeq at a duplication rate, suggesting more usable data after deduplication. For WGS short variant calling, GenoLab M showed significant accuracy improvement over the same depth dataset from NovaSeq, and reached similar accuracy to NovaSeq 33X dataset with 22x depth. For 100X WES, the F-score and Precision in GenoLab M were higher than NovaSeq or NextSeq, especially for InDel calling.

Conclusions

GenoLab M is a promising NGS platform for high-performance WGS and WES applications. For WGS, 22X depth in the GenoLab M sequencing platform offers a cost-effective alternative to the current mainstream 33X depth on Illumina.

How do emerging long-read sequencing technologies function in transforming the plant pathology research landscape?

Long-read sequencing technologies are revolutionizing the sequencing and analysis of plant and pathogen genomes and transcriptomes, as well as contributing to emerging areas of interest in plant-pathogen interactions, disease management techniques, and the introduction of new plant varieties or cultivars. Long-read sequencing (LRS) technologies are progressively being implemented to study plants and pathogens of agricultural importance, which have substantial economic effects. The variability and complexity of the genome and transcriptome affect plant growth, development and pathogen responses. Overcoming the limitations of second-generation sequencing, LRS technology has significantly increased the length of a single contiguous read from a few hundred to millions of base pairs. Because of the longer read lengths, new analysis methods and tools have been developed for plant and pathogen genomics and transcriptomics. LRS technologies enable faster, more efficient, and high-throughput ultralong reads, allowing direct sequencing of genomes that would be impossible or difficult to investigate using short-read sequencing approaches. These benefits include genome assembly in repetitive areas, creating more comprehensive and exact genome determinations, assembling full-length transcripts, and detecting DNA and RNA alterations. Furthermore, these technologies allow for the identification of transcriptome diversity, significant structural variation analysis, and direct epigenetic mark detection in plant and pathogen genomic regions. LRS in plant pathology is found efficient for identifying and characterization of effectors in plants as well as known and unknown plant pathogens. In this review, we investigate how these technologies are transforming the landscape of determination and characterization of plant and pathogen genomes and transcriptomes efficiently and accurately. Moreover, we highlight potential areas of interest offered by LRS technologies for future study into plant-pathogen interactions, disease control strategies, and the development of new plant varieties or cultivars.

Genome-wide identification, molecular evolution and expression analysis of the non-specific lipid transfer protein (nsLTP) family in Setaria italica

BACKGROUND

Foxtail millet (Setaria italica L.) is a millet species with high tolerance to stressful environments. Plant non-specific lipid transfer proteins (nsLTPs) are a kind of small, basic proteins involved in many biological processes. So far, the genome of S. italica has been fully sequenced, and a comprehensive understanding of the evolution and expression of the nsLTP family is still lacking in foxtail millet.

RESULTS

Forty-five nsLTP genes were identified in S. italica and clustered into 5 subfamilies except three single genes (SinsLTP38, SinsLTP7, and SinsLTP44). The proportion of SinsLTPs was different in each subfamily, and members within the same subgroup shared conserved exon-intron structures. Besides, 5 SinsLTP duplication events were investigated. Both tandem and segmental duplication contributed to nsLTP expansion in S. italica, and the duplicated SinsLTPs had mainly undergone purifying selection pressure, which suggested that the function of the duplicated SinsLTPs might not diverge much. Moreover, we identified the nsLTP members in 5 other monocots, and 41, 13, 10, 4, and 1 orthologous gene pairs were identified between S. italica and S. viridis, S. bicolor, Z. mays, O. sativa, and B. distachyon, respectively. The functional divergence within the nsLTP orthologous genes might be limited. In addition, the tissue-specific expression patterns of the SinsLTPs were investigated, and the expression profiles of the SinsLTPs in response to abiotic stress were analyzed, all the 10 selected SinsLTPs were responsive to drought, salt, and cold stress. Among the selected SinsLTPs, 2 paired duplicated genes shared almost equivalent expression profiles, suggesting that these duplicated genes might retain some essential functions during subsequent evolution.

CONCLUSIONS

The present study provided the first systematic analysis for the phylogenetic classification, conserved domain and gene structure, expansion pattern, and expression profile of the nsLTP family in S. italica. These findings could pave a way for further comparative genomic and evolution analysis of nsLTP family in foxtail millet and related monocots, and lay the foundation for the functional analysis of the nsLTPs in S. italica.

Multi-Layer Biosensor for Pre-Symptomatic Detection of Puccinia strifformis, the Causal Agent of Yellow Rust

The yellow rust of wheat (caused by Puccinia striiformis f. sp. tritici) is a devastating fungal infection that is responsible for significant wheat yield losses. The main challenge with the detection of this disease is that it can only be visually detected on the leaf surface between 7 and 10 days after infection, and by this point, counter measures such as the use of fungicides are generally less effective. The hypothesis of this study is to develop and use a compact electrochemical-based biosensor for the early detection of P. striiformis, thus enabling fast countermeasures to be taken. The biosensor that was developed consists of three layers. The first layer mimics the wheat leaf surface morphology. The second layer consists of a sucrose/agar mixture that acts as a substrate and contains a wheat-derived terpene volatile organic compound that stimulates the germination and growth of the spores of the yellow rust pathogen P. s. f. sp. tritici. The third layer consists of a nonenzymatic glucose sensor that produces a signal once invertase is produced by P. striiformis, which comes into contact with the second layer, thereby converting sucrose to glucose. The results show the proof that this innovative biosensor can enable the detection of yellow rust spores in 72 h.

A dominant gene Ihrl1 is tightly linked to and inhibits the gene Ndhrl1 mediating nitrogen-dependent hypersensitive reaction-like phenotype in wheat

Identification and mapping of an inhibitor of Ndhrl1 mediating nitrogen-dependent hypersensitive reaction-like phenotype in wheat. Hypersensitive reaction-like (HRL) traits are characteristic of spontaneous lesions including yellowish spots, brown spots or white-stripe that appeared randomly and dispersedly on all the leaves in the absence of plant pathogens. Our previous studies have shown that the wheat line P7001 showed an HRL trait at low nitrogen supply, and this trait was controlled by the gene Ndhrl1 (Nitrogen-dependent hypersensitive reaction-like 1). In order to investigate the robustness of the trait expression mediated by Ndhrl1 under different genetic backgrounds, seven genetic populations, with P7001 being the common female parent, were constructed and analyzed. F₁ plants from six of the seven combinations showed HRL trait and Ndhrl1 segregated in a dominant way of HRL: non-HRL = 3:1 in the six populations (F₂). Exceptionally, the F₁ plants of P7001/Fielder combination showed non-HRL trait and HRL trait in the F₂ population showed a contrasting recessive segregation ratio of HRL: non-HRL = 1:3, suggesting Fielder may have another HRL-related gene. Using 55 K SNP array and PCR-based markers, the HRL-related gene in Fielder was mapped to an interval of 5.63-12.91 Mb on the short arm of chromosome 2B with the flanking markers Yzu660R075552 and Yzu660F075941. A recombinant with genomic region of Fielder at Ndhrl1 locus showing HRL trait demonstrated that Fielder also harbored Ndhrl1 same as P7001. Thus, Fielder carries a single dominant suppressor of Ndhrl1, designated as Ihrl1 (Inhibitor of hypersensitive reaction-like). Interestingly, Ihrl1 is tightly linked to Ndhrl1 and may be also involved in nitrogen metabolic and (or) signaling pathways.

Transcriptome analysis and mining of genes related to shade tolerance in foxtail millet (Setaria italica (L.) P. Beauv.)

A stereo interplanting system with foxtail millet beneath chestnut trees is an effective planting method to raise the utilization of land in chestnut orchards, increase yields and improve quality of chestnut nuts. Consequently, exploration of genes involved in shade tolerance response in foxtail millet and breeding shade-tolerant varieties have become urgent issues. In this study, RNA-seq of leaf samples from two shade-tolerant varieties and three shade-intolerant varieties of foxtail millet at the booting stage was performed. Comparisons between the varieties revealed that 70 genes were commonly differentially expressed. Moreover, the ratio of net photosynthetic rate under shaded environment to that under light environment could be used as an indicator of shade tolerance. Subsequently, weighted gene co-expression network analysis was employed to construct a co-expression network and modules were correlated with this ratio. A total of 375 genes were identified as potentially relevant to shade tolerance, among which nine genes were also present in the 70 differentially expressed genes, which implied that they were good candidates for genes involved in shade tolerance. Our results provide valuable resources for elucidation of the molecular mechanisms underlying shade tolerance and will contribute to breeding of shade-tolerant foxtail millet that are adapted to the shaded environment under chestnut trees.

Regulatory non-coding RNA: The core defense mechanism against plant pathogens

Plant pathogens damage crops and threaten global food security. Plants have evolved complex defense networks against pathogens, using crosstalk among various signaling pathways. Key regulators conferring plant immunity through signaling pathways include protein-coding genes and non-coding RNAs (ncRNAs). The discovery of ncRNAs in plant transcriptomes was first considered "transcriptional noise". Recent reviews have highlighted the importance of non-coding RNAs. However, understanding interactions among different types of noncoding RNAs requires additional research. This review attempts to consider how long-ncRNAs, small-ncRNAs and circular RNAs interact in response to pathogenic diseases within different plant species. Developments within genomics and bioinformatics could lead to the further discovery of plant ncRNAs, knowledge of their biological roles, as well as an understanding of their importance in exploiting the recent molecular-based technologies for crop protection.

Identification and genetic diversity analysis of hybrid offspring of azalea based on EST-SSR markers

Azalea is a world famous flower with high economic and ornamental value. The breeding of new azalea varieties is mainly done by cross breeding. However, there is a risk that cross breeding may cause errors in the hybrid offspring due to contamination by exogenous pollen. Therefore, the identification of hybrid offspring is an important part of azalea breeding. In this study, the parents of three hybrid combinations and their 88 F1 hybrid offspring were selected to screen 15 pairs of EST-SSR primers to identify the authenticity of azalea hybrid offspring. The results showed that the authenticity of 88 azalea F1 hybrid progenies could be determined by at least four primer pairs. Genetic diversity analysis of azalea hybrid progeny revealed that the number of alleles and polymorphic information content of the progeny increased to different degrees, and the more distant the genetic distance between parents, the richer the polymorphic information. It is suggested that EST-SSR molecular marker can be applied for the early identification and genetic diversity analysis of the progeny of azalea hybrids. This method is of positive significance for improving the breeding efficiency of new varieties and exploring the genetic background of azalea.

Genome-wide identification and analysis of ACP gene family in Sorghum bicolor (L.) Moench

BACKGROUND

Acyl carrier proteins (ACP) constitute a very conserved carrier protein family. Previous studies have found that ACP not only takes part in the fatty acid synthesis process of almost all organisms, but also participates in the regulation of plant growth, development, and metabolism, and makes plants adaptable to stresses. However, this gene family has not been systematically studied in sorghum.

RESULTS

Nine ACP family members were identified in the sorghum genome, which were located on chromosomes 1, 2, 5, 7, 8 and 9, respectively. Evolutionary analysis among different species divided the ACP family into four subfamilies, showing that the SbACPs were more closely related to maize. The prediction results of subcellular localization showed that SbACPs were mainly distributed in chloroplasts and mitochondria, while fluorescence localization showed that SbACPs were mainly localized in chloroplasts in tobacco leaf. The analysis of gene structure revealed a relatively simple genetic structure, that there were 1-3 introns in the sorghum ACP family, and the gene structure within the same subfamily had high similarity. The amplification method of SbACPs was mainly large fragment replication, and SbACPs were more closely related to ACPs in maize and rice. In addition, three-dimensional structure analysis showed that all ACP genes in sorghum contained four α helices, and the second helix structure was more conserved, implying a key role in function. Cis-acting element analysis indicated that the SbACPs might be involved in light response, plant growth and development regulation, biotic and abiotic stress response, plant hormone regulation, and other physiological processes. What's more, qRT-PCR analysis uncovered that some of SbACPs might be involved in the adaptive regulation of drought and salt stresses, indicating the close relationship between fatty acids and the resistance to abiotic stresses in sorghum.

CONCLUSIONS

In summary, these results showed a comprehensive overview of the SbACPs and provided a theoretical basis for further studies on the biological functions of SbACPs in sorghum growth, development and abiotic stress responses.

Transcriptome analysis and characterization of genes associated to leaf tannin content in foxtail millet [Setaria italica (L.) P. Beauv.]

Background

Chinese chestnut is an economically important tree species whose yield and quality are seriously affected by red spider attack. Tannins is one of the most important class secondary metabolites in plants, and is closely associated with plant defense mechanisms against insect and herbivory. In our previous studies, it was revealed that several low-tannin foxtail millet varieties growing under the Chinese chestnut trees could attract red spiders to feed on their leaves and protect the chestnut trees from the infestation of red spiders, meanwhile, the growth and yield of foxtail millet plants themselves were not greatly affected.

Results

To identify genes related to leaf tannin content and selection of foxtail millet germplasm resources with low tannin content for interplanting with Chinese chestnut and preventing the red spider attack, the leaves of 4 varieties with different levels of tannin content were harvested for comparative transcriptome analysis. In total, 335 differentially expressed genes (DEGs) were identified. For acquisition of gene functions and biological pathways they involved in, gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analyses were performed, and several DEGs were found to possibly participate in the tannins biosynthesis pathway and transport processes of precursors. In addition, according to the PlantTFDB database, some transcription factors were predicted among the DEGs, suggesting their role in regulation of tannins biosynthesis pathway.

Conclusion

Our results provide valuable gene resources for understanding the biosynthesis and regulation mechanisms of tannins in foxtail millet, and pave the way for speeding up the breeding of low-tannin varieties through marker-assisted selection, which could be utilized for interplanting with Chinese chestnut trees to confer protection against red spider attack.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08746-8.

Underlying mechanism of accelerated cell death and multiple disease resistance in a maize lethal leaf spot 1 allele

Multiple disease resistance (MDR) in maize has attracted increasing attention. However, the interplay between cell death and metabolite changes and their contributions to MDR remains elusive in maize. In this study, we identified a mutant named as lesion mimic 30 (les30) that showed 'suicidal' lesion formation in the absence of disease and had enhanced resistance to the fungal pathogen Curvularia lunata. Using map-based cloning, we identified the causal gene encoding pheophorbide a oxidase (PAO), which is known to be involved in chlorophyll degradation and MDR, and is encoded by LETHAL LEAF SPOT1 (LLS1). LLS1 was found to be induced by both biotic and abiotic stresses. Transcriptomics analysis showed that genes involved in defense responses and secondary metabolite biosynthesis were mildly activated in leaves of the les30 mutant without lesions, whilst they were strongly activated in leaves with lesions. In addition, in les30 leaves with lesions, there was overaccumulation of defense-associated phytohormones including jasmonic acid and salicylic acid, and of phytoalexins including phenylpropanoids, lignin, and flavonoids, suggesting that their biosynthesis was activated in a lesion-dependent manner. Taken together, our study implies the existence of an interactive amplification loop of interrupted chlorophyll degradation, cell death, expression of defense-related genes, and metabolite changes that results in suicidal lesion formation and MDR, and this has the potential to be exploited by genetic manipulation to improve maize disease resistance.

From a large-scale genomic analysis of insertion sequences to insights into their regulatory roles in prokaryotes

Background

Insertion sequences (ISs) are mobile repeat sequences and most of them can copy themselves to new host genome locations, leading to genome plasticity and gene regulation in prokaryotes. In this study, we present functional and evolutionary relationships between IS and neighboring genes in a large-scale comparative genomic analysis.

Results

IS families were located in all prokaryotic phyla, with preferential occurrence of IS3, IS4, IS481, and IS5 families in Alpha-, Beta-, and Gammaproteobacteria, Actinobacteria and Firmicutes as well as in eukaryote host-associated organisms and autotrophic opportunistic pathogens. We defined the concept of the IS-Gene couple (IG), which allowed to highlight the functional and regulatory impacts of an IS on the closest gene. Genes involved in transcriptional regulation and transport activities were found overrepresented in IG. In particular, major facilitator superfamily (MFS) transporters, ATP-binding proteins and transposases raised as favorite neighboring gene functions of IS hotspots. Then, evolutionary conserved IS-Gene sets across taxonomic lineages enabled the classification of IS-gene couples into phylum, class-to-genus, and species syntenic IS-Gene couples. The IS5, IS21, IS4, IS607, IS91, ISL3 and IS200 families displayed two to four times more ISs in the phylum and/or class-to-genus syntenic IGs compared to other IS families. This indicates that those families were probably inserted earlier than others and then subjected to horizontal transfer, transposition and deletion events over time. In phylum syntenic IG category, Betaproteobacteria, Crenarchaeota, Calditrichae, Planctomycetes, Acidithiobacillia and Cyanobacteria phyla act as IS reservoirs for other phyla, and neighboring gene functions are mostly related to transcriptional regulators. Comparison of IS occurrences with predicted regulatory motifs led to ~ 26.5% of motif-containing ISs with 2 motifs per IS in average. These results, concomitantly with short IS-Gene distances, suggest that those ISs would interfere with the expression of neighboring genes and thus form strong candidates for an adaptive pairing.

Conclusions

All together, our large-scale study provide new insights into the IS genetic context and strongly suggest their regulatory roles.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08678-3.

Characterization and fine mapping of a lesion mimic mutant (Lm5) with enhanced stripe rust and powdery mildew resistance in bread wheat (Triticum aestivum L.)

A novel light intensity-dependent lesion mimic mutant with enhanced disease resistance was physiologically, biochemically, and genetically characterized, and the causative gene was fine mapped to a 1.28 Mbp interval containing 17 high-confidence genes. Lesion mimic mutants are ideal for studying disease resistance and programmed cell death photosynthesis in plants to improve crop yield. In this study, a novel light intensity-dependent lesion mimic mutant (MC21) was obtained from the wheat variety Chuannong16 (CN16) by ethyl methane sulfonate treatment. The mutant initially developed tiny lesion spots on the basal part of the leaves, which then gradually proceeded down to leaf sheaths, stems, shells, and awns at the flowering stage. The major agronomic traits were significantly altered in the mutant compared to that in the wild-type CN16. Furthermore, the mutant exhibited a lesion phenotype with degenerated chloroplast structure, decreased chlorophyll content, increased level of reactive oxygen species, and increased resistance to stripe rust and powdery mildew. Genetic analysis indicated that the lesion phenotype was controlled by a novel single semi-dominant nuclear gene. The target gene was mapped on chromosome arm 2AL located between Kompetitive Allele Specific PCR (KASP) markers, KASP-4211 and KASP-5353, and tentatively termed as lesion mimic 5 (Lm5). The fine mapping suggested that Lm5 was located in a 1.28 Mbp interval between markers KASP-5825 and KASP-9366; 17 high-confidence candidate genes were included in this genomic region. This study provides an important foundational step for further cloning of Lm5 using a map-based approach.

Mapping and CRISPR homology-directed repair of a recessive white eye mutation in Plodia moths

The pantry moth Plodia interpunctella is a worldwide pest of stored food products and a promising laboratory model system for lepidopteran functional genomics. Here we describe efficient methods for precise genome editing in this insect. A spontaneous recessive white-eyed phenotype maps to a frameshift deletion (c.737delC) in the white gene. CRISPR NHEJ mutagenesis of white replicates this phenotype with high rates of somatic biallelic knockout. G₀ individuals with mutant clones on both eyes produced 100% mutant progeny, making white an ideal marker for co-conversion when targeting other genes. CRISPR HDR experiments corrected c.737delC and reverted white eyes to a pigmented state in 37% of G₀ mosaic adults. These repaired alleles showed practical rates of germline transmission in backcrosses, demonstrating the potential of the technique for precise genome editing. Plodia offers a promising avenue for research in this taxon because of its lab-ready features, egg injectability, and editability.

•

Plodia pantry moths are an emerging model organism for functional genomics in Lepidoptera

•

Spontaneous and CRISPR-induced white mutations yield recessive-white eye phenotypes

•

CRISPR HDR repair with ssODN donor result in practical rates of base editing

•

We provide optimized protocols for Plodia handling and genome editing

A dominant spotted leaf gene TaSpl1 activates endocytosis and defense-related genes causing cell death in the absence of dominant inhibitors

The spotted leaf lesion mimic trait simulates cell death in a plant responding to pathogen infection. Some spotted leaf genes are recessive, while others are dominant. A small number of plants with a lesion mimic phenotype appeared in a segregating population obtained by crossing two normal green wheat strains, XN509 and N07216. Here, we clarified the genetic model and its breeding value. Phenotyping of the consecutive progeny populations over six cropping seasons showed that the spotted leaf lesion mimic phenotype was controlled by a dominant gene designated TaSpl1, which was inhibited by two other dominant genes, designated TaSpl1-I₁ and TaSpl1-I₂. Using bulked segregant analysis RNA-seq (BSR-Seq) and newly developed KASP-PCR markers, the TaSpl1 and TaSpl1-I₁ loci in N07216 were mapped to the end of chromosomes 3DS and 3BS, respectively. Plants with the spotted phenotype showed lower levels of stripe rust and powdery mildew than those with the normal green phenotype. Compared with normal leaves, the differentially expressed genes in spotted leaves were significantly enriched in plant-pathogen interaction and endocytosis pathways. There were no differences in the yield parameters of the F₅ and F₆ sister lines, N13039S with TaSpl1 and N13039 N without TaSpl1. These results provide a greater understanding of spotted leaf phenotyping and the breeding value of the lesion mimic allele in developing disease-resistance varieties.