Identification Of Candidate Genes Research Articles

Identification of candidate genes for Fusarium head blight resistance from QTLs using RIL population in wheat.

Fusarium head blight (FHB) stands out as one of the most devastating wheat diseases and leads to significantly grain yield losses and quality reductions in epidemic years. Exploring quantitative trait loci (QTL) for FHB resistance is a critical step for developing new FHB-resistant varieties. We previously constructed a genetic map of unigenes (UG-Map) according to the physical positions using a set of recombinant-inbred lines (RILs) derived from the cross of 'TN18 × LM6' (TL-RILs). Here, the number of diseased spikelets (NDS) and relative disease index (RDI) for FHB resistance were investigated under four environments using TL-RILs, which were distributed across 13 chromosomes. A number of 36 candidate genes for NDS and RDI from of 19 stable QTLs were identified. The average number of candidate genes per QTL was 1.89, with 14 (73.7%), two (10.5%), and three (15.8%) QTLs including one, two, and 3-10 candidate genes, respectively. Among the 24 candidate genes annotated in the reference genome RefSeq v1.1, the homologous genes of seven candidate genes, including TraesCS4B02G227300 for QNds/Rdi-4BL-4553, TraesCS5B02G303200, TraesCS5B02G303300, TraesCS5B02G303700, TraesCS5B02G303800 and TraesCS5B02G304000 for QNds/Rdi-5BL-9509, and TraesCS7A02G568400 for QNds/Rdi-7AL-14499, were previously reported to be related to FHB resistance in wheat, barely or Brachypodium distachyon. These genes should be closely associated with FHB resistance in wheat. In addition, the homologous genes of five genes, including TraesCS1A02G037600LC for QNds-1AS-2225, TraesCS1D02G017800 and TraesCS1D02G017900 for QNds-1DS-527, TraesCS1D02G018000 for QRdi-1DS-575, and TraesCS4B02G227400 for QNds/Rdi-4BL-4553, were involved in plant defense responses against pathogens. These genes should be likely associated with FHB resistance in wheat.

Genome-wide association analysis uncovers rice blast resistance alleles of Ptr and Pia

A critical step to maximize the usefulness of genome-wide association studies (GWAS) in plant breeding is the identification and validation of candidate genes underlying genetic associations. This is of particular importance in disease resistance breeding where allelic variants of resistance genes often confer resistance to distinct populations, or races, of a pathogen. Here, we perform a genome-wide association analysis of rice blast resistance in 500 genetically diverse rice accessions. To facilitate candidate gene identification, we produce de-novo genome assemblies of ten rice accessions with various rice blast resistance associations. These genome assemblies facilitate the identification and functional validation of novel alleles of the rice blast resistance genes Ptr and Pia. We uncover an allelic series for the unusual Ptr rice blast resistance gene, and additional alleles of the Pia resistance genes RGA4 and RGA5. By linking these associations to three thousand rice genomes we provide a useful tool to inform future rice blast breeding efforts. Our work shows that GWAS in combination with whole-genome sequencing is a powerful tool for gene cloning and to facilitate selection of specific resistance alleles for plant breeding.

Identification of soybean mosaic virus strain SC7 resistance loci and candidate genes in soybean [Glycine max (L.) Merr.

Soybean [Glycine max (L.) Merr.] is an important legume crop worldwide, which provides abundant plant protein and oil for human beings. Soybean mosaic virus (SMV) can cause serious damage to the yield and quality of soybean, but it is difficult to control SMV with chemicals, breeding SMV-resistant varieties has become the most effective way to control the disease. Therefore, it is important to identify SMV resistance genes from soybean resources and apply them to soybean breeding. In this study, the disease rates (DRs) of 219 soybean accessions to SMV strain SC7 in two environments were investigated. A high-density NJAU 355K SoySNP array was used for genome-wide association study (GWAS) of DR. A 274kb region on chromosome 15 (1,110,567bp to 1,384,173bp) was repeatedly detected in two environments. Six new significant single nucleotide polymorphisms (SNPs) on chromosome 15 were identified. Four of these six SNPs were located within two candidate genes, Glyma.15G015700 and Glyma.15G015800. The elite haplotype Glyma.15G015700Hap I with low DR exhibited strong resistance to SC7. The expression of Glyma.15G015700 in the SMV-resistant accession increased significantly after inoculation with SC7. Furthermore, most of the proteins predicted to interact with Glyma.15G015700 are heat shock proteins, which have been shown to be related to disease resistance. In summary, new SMV resistance loci and a new candidate gene, Glyma.15G015700, were identified and might be utilized in further soybean disease resistance breeding.

Identification of potential candidate genes for the Huoyan trait in developing Wulong goose embryos by transcriptomic analysis

ABSTRACT 1. The Wulong goose is a Chinese breed and a source of high-quality meat and eggs. A characteristic of the Wulong goose is that a proportion of the birds do not have eyelids, known as the Huoyon trait. 2. Wulong geese exhibiting the Huoyan trait at embryonic stages of 9 days (E9), 12 days (E12) and 14 days (E14) were selected alongside those with normal eyelids for comprehensive transcriptome sequencing. Differentially expressed gene (DEG) and functional enrichment analyses were performed and finally, eight DEG were chosen to verify the accuracy of qPCR sequencing. 3. Overall, 466, 962 and 550 DEG were obtained from the three control groups, D9 vs. N9, D12 vs. N12 and D14 vs. N14, respectively, by differential analysis (p < 0.05). CDKN1C, CRH, CROCC and TYSND1 were significantly expressed in the three groups. Enrichment analysis revealed the enrichment of CROCC and TYSND1 in pathways of cell cycle process, endocytosis, microtubule-based process, microtubule organising centre organisation, protein processing and protein maturation. CDKN1C and CRH were enriched in the cell cycle and cAMP signalling pathway. 4. Some collagen family genes were detected among the DEGs, including COL3A1, COL4A5, COL4A2 and COL4A1. FREM1 and FREM2 genes were detected in both Huoyan and normal eyelids. There was a significant difference (p < 0.01) in FREM1 expression between ED9 and ED14 in female embryos, but this difference was not observed in male embryos.

Genomic regions, candidate genes, and pleiotropic variants associated with physiological and anatomical indicators of heat stress response in lactating sows

BackgroundHeat stress (HS) poses significant threats to the sustainability of livestock production. Genetically improving heat tolerance could enhance animal welfare and minimize production losses during HS events. Measuring phenotypic indicators of HS response and understanding their genetic background are crucial steps to optimize breeding schemes for improved climatic resilience. The identification of genomic regions and candidate genes influencing the traits of interest, including variants with pleiotropic effects, enables the refinement of genotyping panels used to perform genomic prediction of breeding values and contributes to unraveling the biological mechanisms influencing heat stress response. Therefore, the main objectives of this study were to identify genomic regions, candidate genes, and potential pleiotropic variants significantly associated with indicators of HS response in lactating sows using imputed whole-genome sequence (WGS) data. Phenotypic records for 18 traits and genomic information from 1,645 lactating sows were available for the study. The genotypes from the PorcineSNP50K panel containing 50,703 single nucleotide polymorphisms (SNPs) were imputed to WGS and after quality control, 1,622 animals and 7,065,922 SNPs were included in the analyses.ResultsA total of 1,388 unique SNPs located on sixteen chromosomes were found to be associated with 11 traits. Twenty gene ontology terms and 11 biological pathways were shown to be associated with variability in ear skin temperature, shoulder skin temperature, rump skin temperature, tail skin temperature, respiration rate, panting score, vaginal temperature automatically measured every 10 min, vaginal temperature measured at 0800 h, hair density score, body condition score, and ear area. Seven, five, six, two, seven, 15, and 14 genes with potential pleiotropic effects were identified for indicators of skin temperature, vaginal temperature, animal temperature, respiration rate, thermoregulatory traits, anatomical traits, and all traits, respectively.ConclusionsPhysiological and anatomical indicators of HS response in lactating sows are heritable but highly polygenic. The candidate genes found are associated with important gene ontology terms and biological pathways related to heat shock protein activities, immune response, and cellular oxidative stress. Many of the candidate genes with pleiotropic effects are involved in catalytic activities to reduce cell damage from oxidative stress and cellular mechanisms related to immune response.

Study on the Characteristics of Coarse Feeding Tolerance of Ding’an Pigs: Phenotypic and Candidate Genes Identification

Ding’an (DA) pig, a prominent local breed in Hainan Province, exhibits notable advantages in coarse feeding tolerance and high-quality meat. To explore the potential genetic mechanism of coarse feeding tolerance in DA pigs, 60-day-old full sibling pairs of DA and DLY (Duroc-Landrace-Yorkshire) pigs were subjected to fed normal (5%) and high (10%) crude fiber diets for 56 days, respectively. The findings showed that increasing the crude fiber level had no impact on the apparent digestibility of crude fiber, intramuscular fat, and marbling scores in DA pigs, whereas these factors were significantly reduced in DLY pigs (p &lt; 0.05). Through differential expression analysis and Weighted Gene Co-expression Network Analysis (WGCNA) of the colonic mucosal transcriptome data, 65 and 482 candidate genes with coarse feeding tolerance in DA pigs were identified, respectively. Joint analysis screened four key candidate genes, including LDHB, MLC1, LSG1, and ESM1, potentially serving as key regulated genes for coarse feeding tolerance. Functional analysis revealed that the most significant pathway enriched in differential genes associated with coarse feeding tolerance in Ding’an pigs was the signaling receptor binding. The results hold substantial significance for advancing our understanding of the genetic mechanisms governing coarse feeding tolerance in Ding’an pigs.

Transcriptomic responses to drought stress in the Patagonian southern beech Nothofagus alpina

BackgroundDeciphering the genetic architecture of drought tolerance could allow the candidate genes identification responding to water stress. In the Andean Patagonian forest, the genus Nothofagus represents an ecologically relevant species to be included in different genomic studies. These studies are scarce in South American ecosystems however represent an important source of genomic data in order to interpret future climate-change environment scenarios of these emblematic forests. Here, we achieved the assemblage of the transcriptome of N. alpina while searching for key genes of activated or suppressed metabolic pathways in response to drought stress.ResultsDe novo transcriptome assembly resulted in 104,030 transcripts. Following confirmation of drought conditions, based on reduction of leaf water potential and stomatal conductance, a differential gene expression analysis resulted in 2720 significantly expressed genes (1601 up-regulated and 1119 down-regulated). Enrichment analysis (over-representation analysis and gene set enrichment analysis) resulted in more than one hundred stress-responsive term ontologies (i.e. biological processes) and pathways. Terms such as response to abscisic acid and pathways such as plant hormone signal transduction or starch and sucrose metabolism were over-represented. Protein–protein interaction assessment resulted in networks with significantly expressed top common hub gene clusters (e.g. plant-type cell wall biogenesis among down-regulated or ABA-signalling among up-regulated). These networks evidenced important regulators at gene expression such as transcriptional factors.ConclusionsResponses of N. alpina seedlings to drought stress were evidenced by the activation of several genes linked to GO biological processes and KEGG pathways, which were mainly based on over-expression of specific protein kinases, phosphatases, synthases and transcription factors. This suggests an up-regulation of signalling pathways, triggered through plant hormones such as abscisic acid or auxin, which could counteract the osmotic stress created as a probable immediate response to drought. On the other hand, groups of carbon fixation genes related to the galactose metabolism, photosynthesis, secondary wall biogenesis, and fatty acid biosynthesis degradation were down-regulated under drought. Overall, our results provide new genomic data for understanding how non-model long-lived trees of Patagonian forests would acclimate to environmental changes.

Identification of candidate genes for adult plant stripe rust resistance transferred from Aegilops ventricosa 2NvS into wheat via fine mapping and transcriptome analysis.

An adult plant gene for resistance to stripe rust was narrowed down to the proximal one-third of the 2NvS segment translocated from Aegilops ventricosa to wheat chromosome arm 2AS, and based on the gene expression analysis, two candidate genes were identified showing a stronger response at the adult plant stage compared to the seedling stage. The 2NvS translocation from Aegilops ventricosa, known for its resistance to various diseases, has been pivotal in global wheat breeding for more than three decades. Here, we identified an adult plant resistance (APR) gene in the 2NvS segment in wheat line K13-868. Through fine mapping in a segregating near-isogenic line (NIL) derived population of 6389 plants, the candidate region for the APR gene was narrowed down to between 19.36Mb and 33Mb in the Jagger reference genome. Transcriptome analysis in NILs strongly suggested that this APR gene conferred resistance to stripe rust by triggering plant innate immune responses. Based on the gene expression analysis, two disease resistance-associated genes within the candidate region, TraesJAG2A03G00588940 and TraesJAG2A03G00590140, exhibited a stronger response to Puccinia striiformis f. sp. tritici (Pst) infection at the adult plant stage than at the seedling stage, indicating that they could be potential candidates for the resistance gene. Additionally, we developed a co-dominant InDel marker, InDel_31.05, for detecting this APR gene. Applying this marker showed that over one-half of the wheat varieties approved in 2021 and 2022 in Sichuan province, China, carry this gene. Agronomic trait evaluation of NILs indicated that the 2NvS segment effectively mitigated the negative effects of stripe rust on yield without affecting other important agronomic traits. This study provided valuable insights for cloning and breeding through the utilization of the APR gene present in the 2NvS segment.

Deciphering variations, identification of marker-trait associations and candidate genes for seed oil content under terminal heat stress in Indian mustard (Brassica junceaL. Czern & Coss) germplasm stock.

The online version contains supplementary material available at 10.1007/s13205-024-03985-w.

Identification of Candidate Genes for Salt Tolerance at Seedling Stage in Rice Using QTL-Seq and Chromosome Segment Substitution Line-Derived Population

Rice is a staple food for more than half of the world’s population. However, the pervasive problem of salinity is severely undermining rice production, especially in coastal and low-lying areas where soil salinization is widespread. This stress, exacerbated by climate change, necessitates the development of salt-tolerant rice varieties to ensure food security. In this study, an F2:3 population (n = 454) from a cross of KDML105 and its chromosome segment substitution line (CSSL) was used to identify genomic regions associated with salt tolerance at the seedling stage. Using the QTL-seq approach, a QTL significantly associated with salt tolerance was identified on chromosome 1. Annotation of candidate genes in this region revealed the potential regulators of salt tolerance, including MIKC-type MADS domain proteins, calmodulin-binding transcription factors, and NB-ARC domain-containing proteins. These and other identified genes provide insights into the genetic basis of salt tolerance. This study underscores the importance of using advanced genomics tools and CSSL populations in the study of complex traits such as salt tolerance in rice. Several candidate genes identified in this study could be used in further studies on molecular or physiological mechanisms related to the salt response and tolerance mechanism in rice. Additionally, these genes could also be utilized in plant breeding programs for salt tolerance.

Linkage-Mapping and Genome-wide Association Study Identified Two Peanut Late Leaf Spot Resistance Loci, PLLSR-1 and PLLSR-2, Using a Nested Association Mapping.

Identification of candidate genes and molecular markers for late leaf spot (LLS) disease resistance in peanut (Arachis hypogaea) has been a focus of molecular breeding for the U.S. industry funded peanut genome project. Efforts have been hindered by limited mapping resolution due to low levels of genetic recombination and marker density available in traditional biparental mapping populations. To address this, a multi-parental nested association mapping (NAM) population has been genotyped with the peanut 58 K SNP array and phenotyped for LLS severity in the field for three years. Joint linkage-based QTL mapping identified nine QTLs for LLS resistance with significant phenotypic variance explained (PVE) up to 47.7%. A genome-wide association study (GWAS) identified 13 SNPs consistently associated with LLS resistance. Two genomic regions harboring the consistent QTLs and SNPs were identified from 1,336 Kb to 1,520 Kb (184 Kb) on chromosome B02 and from 1,026.9 Kb to 1,793.2 Kb (767 Kb) on chromosome B03, designated as peanut late leaf spot resistance loci, PLLSR-1 and PLLSR-2, respectively. PLLSR-1 contains 10 NBS-LRR disease resistant genes. An NBS-LRR disease resistance gene Arahy.VKVT6A was also identified on homoeologous chromosome A02. PLLSR-2 contains five significant SNPs associated with five different genes encoding callose synthase, pollen defective in guidance protein, pentatricopeptide repeat (PPR), acyl-activating enzyme, and C2 GRAM domains-containing protein. This study highlights the power of multi-parent populations such as NAM for genetic mapping and marker-trait association studies in peanuts. Validation of these two LLS resistance loci will be needed for marker-assisted breeding.

In silico analysis of the wheat BBX gene family and identification of candidate genes for seed dormancy and germination

BackgroundB-box (BBX) proteins are a type of zinc finger proteins containing one or two B-box domains. They play important roles in development and diverse stress responses of plants, yet their roles in wheat remain unclear.ResultsIn this study, 96 BBX genes were identified in the wheat genome and classified into five subfamilies. Subcellular localization prediction results showed that 68 TaBBXs were localized in the nucleus. Protein interaction prediction analysis indicated that interaction was one way that these proteins exerted their functions. Promoter analysis indicated that TaBBXs may play important roles in light signal, hormone, and stress responses. qRT-PCR analysis revealed that 14 TaBBXs were highly expressed in seeds compared with other tissues. These were probably involved in seed dormancy and germination, and their expression patterns were investigated during dormancy acquisition and release in the seeds of wheat varieties Jing 411 and Hongmangchun 21, showing significant differences in seed dormancy and germination phenotypes. Subcellular localization analysis confirmed that the three candidates TaBBX2-2 A, TaBBX4-2 A, and TaBBX11-2D were nuclear proteins. Transcriptional self-activation experiments further demonstrated that TaBBX4-2A was transcriptionally active, but TaBBX2-2A and TaBBX11-2D were not. Protein interaction analysis revealed that TaBBX2-2A, TaBBX4-2A, and TaBBX11-2D had no interaction with each other, while TaBBX2-2A and TaBBX11-2D interacted with each other, indicating that TaBBX4-2A may regulate seed dormancy and germination by transcriptional regulation, and TaBBX2-2A and TaBBX11-2D may regulate seed dormancy and germination by forming a homologous complex.ConclusionsIn this study, the wheat BBX gene family was identified and characterized at the genomic level by bioinformatics analysis. These observations provide a theoretical basis for future studies on the functions of BBXs in wheat and other species.

Weighted gene co-expression network analysis identifies important modules and hub genes involved in the regulation of breast muscle yield in broilers.

Increasing breast meat production is one of the primary goals of the broiler industry. Over the past few decades, tremendous progress has been made in genetic selection and the identification of candidate genes for improving the breast muscle mass. However, the molecular network contributing to muscle production traits in chickens still needs to be further illuminated. A total of 150 1-day-old male 817 broilers were reared in a floor litter system. At the market age of 50 d, eighteen healthy 817 broilers were slaughtered and the left pectoralis major muscle sample from each bird was collected for RNA-seq sequencing. The birds were then plucked and eviscerated and the whole breast muscle was removed and weighed. Breast muscle yield was calculated as the ratio of the breast muscle weight to the eviscerated weight. To identify the co-expression networks and hub genes contributing to breast muscle yield in chickens, we performed weighted gene co-expression network analysis (WGCNA) based on the 18 transcriptome datasets of pectoralis major muscle from eighteen 817 broilers. The WGCNA analysis classified all co-expressed genes in the pectoral muscle of 817 broilers into 44 modules. Among these modules, the turquoise and skyblue3 modules were found to be most significantly positively (r=0.78, p=1e-04) and negatively (r=-0.57, p=0.01) associated with breast meat yield, respectively. Further analysis identified several hub genes (e.g., DLX3, SH3RF2, TPM1, CAV3, MYF6, and CFL2) that involved in muscle structure and muscle development were identified as potential regulators of breast meat production. The present study has advanced our understanding of the molecular regulatory networks contributing to muscle growth and breast muscle production and will contribute to the molecular breeding of chickens in the future.

Identification of candidate genes associated with growth and development of sheep from a crossbred population using genome-wide association studies

The article presents the results of a search for genome-wide associations with phenotypic traits characterizing the growth and development of sheep from a crossbred population obtained from crossing Romanov sheep and F1 hybrid rams (Romanov sheep x Katahdin). The phenotype database included ten body measurements (withers height, sacral height, back height, chest depth, chest width, ischial tuberosity width, body length, oblique body length, chest girth, pastern girth) recorded at the age of 6 days, 3, 6 and 9 months. Genotyping of sheep was carried out using high-density DNA chips containing about 600,000 SNP markers. Genome-wide association studies (GWAS) were performed using regression analysis in the STATISTICA 10 program. The search for candidate genes localized in the SNP region was performed using Ensembl genome browser 110. There was carried out an analysis of the matches of the identified SNPs with known quantitative trait loci (QTLs) described in the Sheep Quantitative Trait Locus Database. There were found SNPs that were significantly associated with the studied phenotypic traits overlapped with the QTLs, among which the most common categories were “Body weight (slaughter)”, “Muscle weight in carcass”, “Body weight (live)” and “Bone weight in carcass”. There has been established that SNPs significantly associated with exterior traits were localized within or in the immediate vicinity of 64 genes. There were found potential candidates regulating the growth of muscle (FOXO3, PRKAG3, MYOZ2, and ANKRD1) and cartilage tissues FGF12) and involved in metabolic processes, which were critical for the growth of lambs (CLDN, ALB, and MRC1). Along with the known in sheep functional candidates (CAST and SCD5) , there were identified genes that were not previously described in sheep, but regulated growth and development processes in other livestock species including genes RAB28, PRKAG3 and FOXO3. The identified SNPs can be recommended for inclusion in marker-guided selection programs in sheep breeding.

Identification and validation of defense related candidate genes in Sesamum under artificial inoculation of Macrophomina phaseolina

Identification and validation of defense related candidate genes in Sesamum under artificial inoculation of Macrophomina phaseolina

Identification of Chromosomal Regions and Candidate Genes for Round leaf Locus in Cucumis melo L.

Leaf morphology plays a crucial role in plant classification and provides a significant model for studying plant diversity while directly impacting photosynthetic efficiency. In the case of melons, leaf shape not only influences production and classification but also represents a key genetic trait that requires further exploration. In this study, we utilized forward genetics to pinpoint a recessive locus, dubbed Cmrl (Round leaf), which is responsible for regulating melon leaf shape. Through bulked segregant analysis sequencing and extensive evaluation of a two-year F2 population, we successfully mapped the Cmrl locus to a 537.07 kb region on chromosome 8 of the melon genome. Subsequent genetic fine-mapping efforts, leveraging a larger F2 population encompassing 1322 plants and incorporating F2:3 phenotypic data, further refined the locus to an 80.27 kb interval housing five candidate genes. Promoter analysis and coding sequence cloning confirmed that one of these candidates, MELO3C019152.2 (Cmppr encoding a pentatricopeptide repeat-containing family protein, Cmppr), stands out as a strong candidate gene for the Cmrl locus. Notably, comparisons of Cmrl expressions across various stages of leaf development and different leaf regions suggest a pivotal role of Cmrl in the morphogenesis of melon leaves.

QTL mapping for the flag leaf-related traits using RILs derived from Trititrigia germplasm line SN304 and wheat cultivar Yannong15 in multiple environments

BackgroundDeveloping and enriching genetic resources plays important role in the crop improvement. The flag leaf affects plant architecture and contributes to the grain yield of wheat (Triticum aestivum L.). The genetic improvement of flag leaf traits faces problems such as a limited genetic basis. Among the various genetic resources of wheat, Thinopyrum intermedium has been utilized as a valuable resource in genetic improvement due to its disease resistance, large spikes, large leaves, and multiple flowers. In this study, a recombinant inbred line (RIL) population was derived from common wheat Yannong15 and wheat-Th. intermedium introgression line SN304 was used to identify the quantitative trait loci (QTL) for flag leaf-related traits.ResultsQTL mapping was performed for flag leaf length (FLL), flag leaf width (FLW) and flag leaf area (FLA). A total of 77 QTLs were detected, and among these, 51 QTLs with positive alleles were contributed by SN304. Fourteen major QTLs for flag leaf traits were detected on chromosomes 2B, 3B, 4B, and 2D. Additionally, 28 QTLs and 8 QTLs for flag leaf-related traits were detected in low-phosphorus and drought environments, respectively. Based on major QTLs of positive alleles from SN304, we identified a pair of double-ended anchor primers mapped on chromosome 2B and amplified a specific band of Th. intermedium in SN304. Moreover, there was a major colocated QTL on chromosome 2B, called QFll/Flw/Fla-2B, which was delimited to a physical interval of approximately 2.9 Mb and contained 20 candidate genes. Through gene sequence and expression analysis, four candidate genes associated with flag leaf formation and growth in the QTL interval were identified.ConclusionThese results promote the fine mapping of QFll/Flw/Fla-2B, which have pleiotropic effects, and will facilitate the identification of candidate genes for flag leaf-related traits. Additionally, this work provides a theoretical basis for the application of Th. intermedium in wheat breeding.

Transcriptome profiling of aerial and subterranean peanut pod development

Peanut (Arachis hypogaea) showcases geocarpic behavior, transitioning from aerial flowering to subterranean seed development. We recently obtained an atavistic variant of this species, capable of producing aerial and subterranean pods on a single plant. Notably, although these pod types share similar vigor levels, they exhibit distinct differences in their physical aspects, such as pod size, color, and shell thickness. We constructed 63 RNA-sequencing datasets, comprising three biological replicates for each of 21 distinct tissues spanning six developmental stages for both pod types, providing a rich tapestry of the pod development process. This comprehensive analysis yielded an impressive 409.36 Gb of clean bases, facilitating the detection of 42,401 expressed genes. By comparing the transcriptomic data of the aerial and subterranean pods, we identified many differentially expressed genes (DEGs), highlighting their distinct developmental pathways. By providing a detailed workflow from the initial sampling to the final DEGs, this study serves as an important resource, paving the way for future research into peanut pod development and aiding transcriptome-based expression profiling and candidate gene identification.

Identification of candidate genes associated with milk production and mastitis based on transcriptome-wide association study.

Genetic research for the assessment of mastitis and milk production traits simultaneously has a long history. The main issue that arises in this context is the known existence of a positive correlation between the risk of mastitis and lactation performance due to selection. The transcriptome-wide association study (TWAS) approach endeavors to combine the expression quantitative trait loci and genome-wide association study summary statistics to decode complex traits or diseases. Accordingly, we used the farmgtex project results as a complete bovine database for mastitis and milk production. The results of colocalization and TWAS approaches were used for the detection of functional associated candidate genes with milk production and mastitis traits on multiple tissue-based transcriptome records. Also, we used the david database for gene ontology to identify significant terms and associated genes. For the identification of interaction networks, the genemania and string databases were used. Also, the available z-scores in TWAS results were used for the calculation of the correlation between tissues. Therefore, the present results confirm that LYNX1, DGAT1, C14H8orf33, and LY6E were identified as significant genes associated with milk production in eight, six, five, and five tissues, respectively. Also, FBXL6 was detected as a significant gene associated with mastitis trait. CLN3 and ZNF34 genes emerged via both the colocalization and TWAS approaches as significant genes for milk production trait. It is expected that TWAS and colocalization can improve our perception of the potential health status control mechanism in high-yielding dairy cows.

Biosynthesis of Hesperetin, Homoeriodictyol, and Homohesperetin in a Transcriptomics-Driven Engineered Strain of Streptomyces albidoflavus.

Flavonoids exhibit various bioactivities including anti-oxidant, anti-tumor, anti-inflammatory, and anti-viral properties. Methylated flavonoids are particularly significant due to their enhanced oral bioavailability, improved intestinal absorption, and greater stability. The heterologous production of plant flavonoids in bacterial factories involves the need for enough biosynthetic precursors to allow for high production levels. These biosynthetic precursors are malonyl-CoA and l-tyrosine. In this work, to enhance flavonoid biosynthesis in Streptomyces albidoflavus, we conducted a transcriptomics study for the identification of candidate genes involved in l-tyrosine catabolism. The hypothesis was that the bacterial metabolic machinery would detect an excess of this amino acid if supplemented with the conventional culture medium and would activate the genes involved in its catabolism towards energy production. Then, by inactivating those overexpressed genes (under an excess of l-tyrosine), it would be possible to increase the intracellular pools of this precursor amino acid and eventually the final flavonoid titers in this bacterial factory. The RNAseq data analysis in the S. albidoflavus wild-type strain highlighted the hppD gene encoding 4-hydroxyphenylpyruvate dioxygenase as a promising target for knock-out, exhibiting a 23.2-fold change (FC) in expression upon l-tyrosine supplementation in comparison to control cultivation conditions. The subsequent knock-out of the hppD gene in S. albidoflavus resulted in a 1.66-fold increase in the naringenin titer, indicating enhanced flavonoid biosynthesis. Leveraging the improved strain of S. albidoflavus, we successfully synthesized the methylated flavanones hesperetin, homoeriodictyol, and homohesperetin, achieving titers of 2.52 mg/L, 1.34 mg/L, and 0.43 mg/L, respectively. In addition, the dimethoxy flavanone homohesperetin was produced as a byproduct of the endogenous metabolism of S. albidoflavus. To our knowledge, this is the first time that hppD deletion was utilized as a strategy to augment the biosynthesis of flavonoids. Furthermore, this is the first report where hesperetin and homoeriodictyol have been synthesized from l-tyrosine as a precursor. Therefore, transcriptomics is, in this case, a successful approach for the identification of catabolism reactions affecting key precursors during flavonoid biosynthesis, allowing the generation of enhanced production strains.