Candidate Genes Research Articles

Genome-wide association study identifies novel loci and candidate genes for rust resistance in wheat (Triticum aestivum L.)

BackgroundWheat rusts are important biotic stresses, development of rust resistant cultivars through molecular approaches is both economical and sustainable. Extensive phenotyping of large mapping populations under diverse production conditions and high-density genotyping would be the ideal strategy to identify major genomic regions for rust resistance in wheat. The genome-wide association study (GWAS) population of 280 genotypes was genotyped using a 35 K Axiom single nucleotide polymorphism (SNP) array and phenotyped at eight, 10, and, 10 environments, respectively for stem/black rust (SR), stripe/yellow rust (YR), and leaf/brown rust (LR).ResultsForty-one Bonferroni corrected marker-trait associations (MTAs) were identified, including 17 for SR and 24 for YR. Ten stable MTAs and their best combinations were also identified. For YR, AX-94990952 on 1A + AX-95203560 on 4A + AX-94723806 on 3D + AX-95172478 on 1A showed the best combination with an average co-efficient of infection (ACI) score of 1.36. Similarly, for SR, AX-94883961 on 7B + AX-94843704 on 1B and AX-94883961 on 7B + AX-94580041 on 3D + AX-94843704 on 1B showed the best combination with an ACI score of around 9.0. The genotype PBW827 have the best MTA combinations for both YR and SR resistance. In silico study identifies key prospective candidate genes that are located within MTA regions. Further, the expression analysis revealed that 18 transcripts were upregulated to the tune of more than 1.5 folds including 19.36 folds (TraesCS3D02G519600) and 7.23 folds (TraesCS2D02G038900) under stress conditions compared to the control conditions. Furthermore, highly expressed genes in silico under stress conditions were analyzed to find out the potential links to the rust phenotype, and all four genes were found to be associated with the rust phenotype.ConclusionThe identified novel MTAs, particularly stable and highly expressed MTAs are valuable for further validation and subsequent application in wheat rust resistance breeding. The genotypes with favorable MTA combinations can be used as prospective donors to develop elite cultivars with YR and SR resistance.

Insights into the genetic architecture of Phytophthora capsici root rot resistance in chile pepper (Capsicum spp.) from multi-locus genome-wide association study

BackgroundPhytophthora root rot, a major constraint in chile pepper production worldwide, is caused by the soil-borne oomycete, Phytophthora capsici. This study aimed to detect significant regions in the Capsicum genome linked to Phytophthora root rot resistance using a panel consisting of 157 Capsicum spp. genotypes. Multi-locus genome wide association study (GWAS) was conducted using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing (GBS). Individual plants were separately inoculated with P. capsici isolates, ‘PWB-185’, ‘PWB-186’, and ‘6347’, at the 4–8 leaf stage and were scored for disease symptoms up to 14-days post-inoculation. Disease scores were used to calculate disease parameters including disease severity index percentage, percent of resistant plants, area under disease progress curve, and estimated marginal means for each genotype.ResultsMost of the genotypes displayed root rot symptoms, whereas five accessions were completely resistant to all the isolates and displayed no symptoms of infection. A total of 55,117 SNP markers derived from GBS were used to perform multi-locus GWAS which identified 330 significant SNP markers associated with disease resistance. Of these, 56 SNP markers distributed across all the 12 chromosomes were common across the isolates, indicating association with more durable resistance. Candidate genes including nucleotide-binding site leucine-rich repeat (NBS-LRR), systemic acquired resistance (SAR8.2), and receptor-like kinase (RLKs), were identified within 0.5 Mb of the associated markers.ConclusionsResults will be used to improve resistance to Phytophthora root rot in chile pepper by the development of Kompetitive allele-specific markers (KASP®) for marker validation, genomewide selection, and marker-assisted breeding.

Function of MYB8 in larch under PEG simulated drought stress

Larch, a prominent afforestation, and timber species in northeastern China, faces growth limitations due to drought. To further investigate the mechanism of larch’s drought resistance, we conducted full-length sequencing on embryonic callus subjected to PEG-simulated drought stress. The sequencing results revealed that the differentially expressed genes (DEGs) primarily played roles in cellular activities and cell components, with molecular functions such as binding, catalytic activity, and transport activity. Furthermore, the DEGs showed significant enrichment in pathways related to protein processing, starch and sucrose metabolism, benzose-glucuronic acid interconversion, phenylpropyl biology, flavonoid biosynthesis, as well as nitrogen metabolism and alanine, aspartic acid, and glutamic acid metabolism. Consequently, the transcription factor T_transcript_77027, which is involved in multiple pathways, was selected as a candidate gene for subsequent drought stress resistance tests. Under PEG-simulated drought stress, the LoMYB8 gene was induced and showed significantly upregulated expression compared to the control. Physiological indices demonstrated an improved drought resistance in the transgenic plants. After 48 h of PEG stress, the transcriptome sequencing results of the transiently transformed LoMYB8 plants and control plants exhibited that genes were significantly enriched in biological process, cellular component and molecular function. Function analyses indicated for the enrichment of multiple KEGG pathways, including energy synthesis, metabolic pathways, antioxidant pathways, and other relevant processes. The pathways annotated by the differential metabolites mainly encompassed signal transduction, carbohydrate metabolism, amino acid metabolism, and flavonoid metabolism.

Whole-genome resequencing reveals genetic diversity and adaptive evolution in Chinese honeybee (Apis cerana cerana) in Guizhou, China

Introduction: Guizhou Province, characterized by complex and diverse geographic and climatic environments, has rich genetic resources for the Chinese honeybee (Apis cerana cerana) and is one of the main bee-producing areas in China. However, research on the genetic diversity of Chinese honeybee in the Guizhou region is very limited, despite implications for conservation of biodiversity.Methods: In this study, we analyzed the genetic diversity, differentiation, and selection signals based on 116 Chinese honeybees from 12 regions in Guizhou Province using whole-genome sequencing.Results: We identified 1,400,430 high-quality SNPs across all samples. A population structure analysis revealed two independent genetic subgroups of Chinese honeybees in Guizhou, a Yunnan-Guizhou Plateau population in western Guizhou and a hilly-mountainous population in eastern Guizhou. The average nucleotide diversity (Pi) ranged from 0.00138 to 0.00161 and average expected heterozygosity (He) ranged from 0.2592 to 0.2604. The average genetic differentiation index (FST) for Chinese honeybees in pairwise comparisons of 12 regions ranged from 0.0094 to 0.0293. There was clear genetic differentiation between the western plateau and the eastern hilly mountainous areas of Guizhou; however, FST values between the eastern and western populations ranged from 0.0170 to 0.0293, indicating a low degree of differentiation. A genome-wide scan revealed a number of genes under selection in the Yunnan-Guizhou Plateau environment. These genes were related to growth and development, reproduction, and cold resistance, and several candidate genes involved in environmental adaptation were identified, including CTR, MAPK, MAST, HSF, and MKKK.Discussion: The results of the present study provide important theoretical bases for the conservation, evaluation, development, and utilization of genetic resources for Chinese honeybees in the Guizhou region and for further investigations of environmental adaptation and underlying mechanisms in the species.

Fine mapping and candidate gene analysis of CRA8.1.6, which confers clubroot resistance in turnip (Brassica rapa ssp. rapa)

Clubroot disease poses a significant threat to Brassica crops, necessitating ongoing updates on resistance gene sources. In F2 segregants of the clubroot-resistant inbred line BrT18-6-4-3 and susceptible DH line Y510, the genetic analysis identified a single dominant gene responsible for clubroot resistance. Through bulk segregant sequencing analysis and kompetitive allele-specific polymerase chain reaction assays, CRA8.1.6 was mapped within 110 kb (12,255–12,365 Mb) between markers L-CR11 and L-CR12 on chromosome A08. We identified B raA08g015220.3.5C as the candidate gene of CRA8.1.6. Upon comparison with the sequence of disease-resistant material BrT18-6-4-3, we found 249 single-nucleotide polymorphisms, seven insertions, six deletions, and a long terminal repeat (LTR) retrotransposon (5,310 bp) at 909 bp of the first intron. However, the LTR retrotransposon was absent in the coding sequence of the susceptible DH line Y510. Given the presence of a non-functional LTR insertion in other materials, it showed that the LTR insertion might not be associated with susceptibility. Sequence alignment analysis revealed that the fourth exon of the susceptible line harbored two deletions and an insertion, resulting in a frameshift mutation at 8,551 bp, leading to translation termination at the leucine-rich repeat domain’s C-terminal in susceptible material. Sequence alignment of the CDS revealed a 99.4% similarity to Crr1a, which indicate that CRA8.1.6 is likely an allele of the Crr1a gene. Two functional markers, CRA08-InDel and CRA08-KASP1, have been developed for marker-assisted selection in CR turnip cultivars. Our findings could facilitate the development of clubroot-resistance turnip cultivars through marker-assisted selection.

Integrated proteomic and genomic analysis to identify predictive biomarkers for valproate response in bipolar disorder: a 6-month follow-up study

BackgroundSeveral genetic studies have been undertaken to elucidate the intricate interplay between genetics and drug responses in bipolar disorder (BD). However, there has been notably limited research on biomarkers specifically linked to valproate, with only a few studies investigating integrated proteomic and genomic factors in response to valproate treatment. Therefore, this study aimed to identify biological markers for the therapeutic response to valproate treatment in BD. Patients with BD in remission were assessed only at baseline, whereas those experiencing acute mood episodes were evaluated at three points (baseline, 8 ± 2 weeks, and 6 ± 1 months). The response to valproate treatment was measured using the Alda scale, with individuals scoring an Alda A score ≥ 5 categorized into the acute-valproate responder (acute-VPAR) group. We analyzed 158 peptides (92 proteins) from peripheral blood samples using multiple reaction monitoring mass spectrometry, and proteomic result-guided candidate gene association analyses, with 1,627 single nucleotide variants (SNVs), were performed using the Korean chip.ResultsThe markers of 37 peptides (27 protein) showed temporal upregulation, indicating possible association with response to valproate treatment. A total of 58 SNVs in 22 genes and 37 SNVs in 16 genes showed nominally significant associations with the Alda A continuous score and the acute-VPAR group, respectively. No SNVs reached the genome-wide significance threshold; however, three SNVs (rs115788299, rs11563197, and rs117669164) in the secreted phosphoprotein 2 gene reached a gene-based false discovery rate-corrected significance threshold with response to valproate treatment. Significant markers were associated with the pathophysiological processes of bipolar disorders, including the immune response, acute phase reaction, and coagulation cascade. These results suggest that valproate effectively suppresses mechanisms associated with disease progression.ConclusionsThe markers identified in this study could be valuable indicators of the underlying mechanisms associated with response to valproate treatment.

The Features of Shared Genes among Transcriptomes Probed in Atopic Dermatitis, Psoriasis, and Inflammatory Acne: S100A9 Selection as the Target Gene.

Atopic dermatitis (AD), psoriasis (PS), and inflammatory acne (IA) are well-known as inflammatory skin diseases. Studies of the transcriptome with altered expression levels have reported a large number of dysregulated genes and gene clusters, particularly those involved in inflammatory skin diseases. To identify genes commonly shared in AD, PS, and IA that are potential therapeutic targets, we have identified consistently dysregulated genes and disease modules that overlap with AD, PS, and IA. Microarray data from AD, PS, and IA patients were downloaded from Gene Expression Omnibus (GEO), and identification of differentially expressed genes from microarrays of AD, PS, and IA was conducted. Subsequently, gene ontology and gene set enrichment analysis, detection of disease modules with known disease-associated genes, construction of the protein-protein interaction (PPI) network, and PPI sub-mapping analysis of shared genes were performed. Finally, the computational docking simulations between the selected target gene and inhibitors were conducted. We identified 50 shared genes (36 up-regulated and 14 down-regulated) and disease modules for each disease. Among the shared genes, 20 common genes in PPI network were detected such as LCK, DLGAP5, SELL, CEP55, CDC20, RRM2, S100A7, S100A9, MCM10, AURKA, CCNB1, CHEK1, BTC, IL1F7, AGTR1, HABP4, SERPINB13, RPS6KA4, GZMB, and TRIP13. Finally, S100A9 was selected as the target gene for therapeutics. Docking simulations between S100A9 and known inhibitors indicated several key binding residues, and based on this result, we suggested several cannabinoids such as WIN-55212-2, JZL184, GP1a, Nabilone, Ajulemic acid, and JWH-122 could be potential candidates for a clinical study for AD, PS, and IA via inhibition of S100A9-related pathway. Overall, our approach may become an effective strategy for discovering new disease candidate genes for inflammatory skin diseases with a reevaluation of clinical data.

Genome-wide association analysis explores the genetic loci of amino acid content in duck’s breast muscle

BackgroundAmino acids are the basic components of protein and an important index to evaluate meat quality. With the rapid development of genomics, candidate regions and genes affecting amino acid content in livestock and poultry have been gradually revealed. Hence, genome-wide association study (GWAS) can be used to screen candidate loci associated with amino acid content in duck meat.ResultIn the current study, the content of 16 amino acids was detected in 358 duck breast muscles. The proportion of Glu to the total amino acid content was relatively high, and the proportion was 0.14. However, the proportion of Met content was relatively low, at just 0.03. By comparative analysis, significant differences were found between males and females in 3 amino acids, including Ser, Met, and Phe. In addition, 12 SNPs were significantly correlated with Pro content by GWAS analysis, and these SNPs were annotated by 7 protein-coding genes; 8 significant SNPs were associated with Tyr content, and these SNPs were annotated by 6 protein-coding genes. At the same time, linkage disequilibrium (LD) analysis was performed on these regions with significant signals. The results showed that three SNPs in the 55–56 Mbp region of chromosome 3 were highly correlated with the leader SNP (chr3:55526954) that affected Pro content (r2 > 0.6). Similarly, LD analysis showed that there were three SNPs in the 21.2–21.6 Mbp region of chromosome 13, which were highly correlated with leader SNP (chr13:21421661) (r2 > 0.6). Moreover, Through functional enrichment analysis of all candidate genes. The results of GO enrichment analysis showed that several significant GO items were associated with amino acid transport function, including amino acid transmembrane transport and glutamine transport. The results further indicate that these candidate genes are closely associated with amino acid transport. Among them, key candidate genes include SLC38A1. For KEGG enrichment analysis, CACNA2D3 and CACNA1D genes were covered by significant pathways.ConclusionIn this study, GWAS analysis found a total of 28 significant SNPs affecting amino acid content. Through gene annotation, a total of 20 candidate genes were screened. In addition, Through LD analysis and enrichment analysis, we considered that SERAC1, CACNA2D3 and SLC38A1 genes are important candidate genes affecting amino acid content in duck breast muscle.

Mendelian randomization and single-cell expression analyses identify the causal relationship between depression and chronic rhinosinusitis

BackgroundThe causative relationship between chronic rhinosinusitis (CRS) and depression remains unclear. Herein we employed Mendelian randomization (MR) coupled with single-cell analysis to investigate the causality between CRS and depression.MethodsData pertaining to CRS and depression were mined from the genome-wide association study database, and a single-cell dataset was sourced from the literature. To explore causality, we conducted bidirectional MR analysis using MR-Egger, weighted median, inverse variance weighted (IVW), simple mode, and weighted mode, with IVW representing the most important method. Further, sensitivity analysis was performed to evaluate the robustness of MR analysis results. Candidate genes were analyzed via single-cell combined MR analysis.ResultsForward MR analysis indicated depression as a risk factor for CRS when depression was the exposure factor and CRS was the outcome (OR = 1.425, P < 0.001). Reverse MR analysis revealed the same positive relationship between CRS and depression when CRS was the exposure factor and depression was the outcome (OR = 1.012, P = 0.038). Sensitivity analysis validated the robustness of bidirectional MR analysis results. Ten cell types (endothelial, ciliated, basal, myeloid, mast, apical, plasma, glandular, fibroblast, and T cells) were identified in the single-cell dataset. The network of receptor–ligand pairs showed that in normal samples, cell–cell interactions were present among various cell types, such as epithelial, mast, myeloid, and endothelial cells. In contrast, CRS samples featured only one specific receptor–ligand pair, confined to myeloid cells. TCF4 and MEF2C emerged as potentially crucial for CRS-associated depression development.ConclusionsOur findings suggest a bidirectional causal relationship between CRS and depression, offering a new perspective on the association between CRS and depression.

Genome-wide detections for runs of homozygosity and selective signatures reveal novel candidate genes under domestication in chickens

BackgroundIndigenous chickens were developed through a combination of natural and artificial selection; essentially, changes in genomes led to the formation of these modern breeds via admixture events. However, their confusing genetic backgrounds include a genomic footprint regulating complex traits, which is not conducive to modern animal breeding.ResultsTo better evaluate the candidate regions under domestication in indigenous chickens, we considered both runs of homozygosity (ROHs) and selective signatures in 13 indigenous chickens. The genomes of Silkie feather chickens presented the highest heterozygosity, whereas the highest inbreeding status and ROH number were found in Luhua chickens. Short ROH (< 1 Mb), were the principal type in all chickens. A total of 291 ROH islands were detected, and QTLdb mapping results indicated that body weight and carcass traits were the most important traits. An ROH on chromosome 2 covering VSTM2A gene was detected in 12 populations. Combined analysis with the Tajima’s D index revealed that 18 genes (e.g., VSTM2A, BBOX1, and RYR2) were under selection and covered by ROH islands. Transcriptional analysis results showed that RYR2 and BBOX1 were specifically expressed in the heart and muscle tissue, respectively.ConclusionBased on genome-wide scanning for ROH and selective signatures, we evaluated the genomic characteristics and detected significant candidate genes covered by ROH islands and selective signatures. The findings in this study facilitated the understanding of genetic diversity and provided valuable insights for chicken breeding and conservation strategies.

Identification of a novel QTL for lean meat percentage using imputed genotypes.

Lean meat percentage is a critical production trait in pig breeding systems with direct implications for the sustainability of the industry. In this study, we conducted a genome-wide association study for lean meat percentage using a cohort of 850 Duroc × (Landrace × Yorkshire) crossbred pigs and we identified QTL on SSC3 and SSC18. Based on the predicted effect of imputed variants and using the PigGTEx database of molecular QTL, we prioritized candidate genes and SNPs located within the QTL regions, which may be involved in the regulation of porcine leanness. Our results indicate that a nonsense mutation in ZC3HAV1L on SSC18 has a direct effect on lean meat percentage.

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.

Population dynamics of potentially harmful haplotypes: a pedigree analysis

BackgroundThe identification of low-frequency haplotypes, never observed in homozygous state in a population, is considered informative on the presence of potentially harmful alleles (candidate alleles), putatively involved in inbreeding depression. Although identification of candidate alleles is challenging, studies analyzing the dynamics of potentially harmful alleles are lacking. A pedigree of the highly endangered Gochu Asturcelta pig breed, including 471 individuals belonging to 51 different families with at least 5 offspring each, was genotyped using the Axiom PigHDv1 Array (658,692 SNPs). Analyses were carried out on four different cohorts defined according to pedigree depth and at the whole population (WP) level.ResultsThe 4,470 Linkage Blocks (LB) identified in the Base Population (10 individuals), gathered a total of 16,981 alleles in the WP. Up to 5,466 (32%) haplotypes were statistically considered candidate alleles, 3,995 of them (73%) having one copy only. The number of alleles and candidate alleles varied across cohorts according to sample size. Up to 4,610 of the alleles identified in the WP (27% of the total) were present in one cohort only. Parentage analysis identified a total of 67,742 parent-offspring incompatibilities. The number of mismatches varied according to family size. Parent-offspring inconsistencies were identified in 98.2% of the candidate alleles and 100% of the LB in which they were located. Segregation analyses informed that most potential candidate alleles appeared de novo in the pedigree. Only 17 candidate alleles were identified in the boar, sow, and paternal and maternal grandparents and were considered segregants.ConclusionsOur results suggest that neither mutation nor recombination are the major forces causing the apparition of candidate alleles. Their occurrence is more likely caused by Allele-Drop-In events due to SNP calling errors. New alleles appear when wrongly called SNPs are used to construct haplotypes. The presence of candidate alleles in either parents or grandparents of the carrier individuals does not ensure that they are true alleles. Minimum Allele Frequency thresholds may remove informative alleles. Only fully segregant candidate alleles should be considered potentially harmful alleles. A set of 16 candidate genes, potentially involved in inbreeding depression, is described.

Full-length transcriptome and co-expression network analysis reveal molecular mechanisms of seed development in Elymus sibiricus

Abstract Grass seeds play a critical and fundamental role in grass breeding and production. Elymus sibiricus L. is a widespread Poaceae forage grass in northern Eurasia which is used for ecological restoration and forage production. Sucrose is the main source of substrate and energy required for starch synthesis in the seeds, so the hydrolysis of sucrose determines and influences starch synthesis and filling in the seeds, especially Poaceae. However, the process behind carbohydrate metabolism during E. sibiricus seed development remains unclear. This study addresses a significant gap in our understanding of the carbohydrate metabolism during seed development in E. sibiricus by employing full-length transcriptome sequencing across five developmental stages for the first time. Full-length transcriptome sequencing was performed on E. sibiricus seeds at five developmental stages (S5, S9, S15, S20, S25) to get better molecular insights. We identified 13,205 differentially expressed genes, with 7,471 up-regulated and 5,734 down-regulated. Through KEGG enrichment analysis, genes were enriched in ‘starch and sucrose metabolism’, ‘photosynthetic-related’ and ‘hormone signal transduction’ pathways. Gene ontology enrichment analysis showed that genes were enriched in the ‘beta-amylase activity’ term of molecular functions. In addition, top 21 transcription factor families were identified as involved in seed development. The homologous genes of ABSCISIC ACID-INSENSITIVE 3 (ABI3), NUCLEAR FACTOR-YB1 (NF-YB1), STARCH SYNTHASE I (SSI) were identified as candidate genes of seed development in E. sibiricus. Combined with physiological index, transcriptome analyses, weighted gene co-expression network analysis and real-time quantitative PCR, the mechanism of starch and sucrose content of seed development was revealed and ten hub genes were identified. Overall, this study provides the molecular bases to understand seed development and starch and sucrose metabolism at the different seed developmental stages in E. sibiricus.

Genome-wide analysis of NAC transcription factors and exploration of candidate genes regulating selenium metabolism in Broussonetia papyrifera.

Genome-wide identification revealed 79 BpNAC genes belonging to 16 subfamilies, and their gene structures and evolutionary relationships were characterized. Expression analysis highlighted their importance in plant selenium stress responses. Paper mulberry (Broussonetia papyrifera), a deciduous arboreal plant of the Moraceae family, is distinguished by its leaves, which are abundant in proteins, polysaccharides, and flavonoids, positioning it as a novel feedstock. NAC transcription factors, exclusive to plant species, are crucial in regulating growth, development, and response to biotic and abiotic stress. However, extensive characterization of the NAC family within paper mulberry is lacking. In this study, 79 BpNAC genes were identified from the paper mulberry genome, with an uneven distribution across 13 chromosomes. A comprehensive, genome-wide analysis of BpNACs was performed, including investigating gene structures, promoter regions, and chromosomal locations. Phylogenetic tree analysis, alongside comparisons with Arabidopsis thaliana NACs, allowed for categorizing these genes into 16 subfamilies in alignment with gene structure and motif conservation. Collinearity analysis suggested a significant homologous relationship between the NAC genes of paper mulberry and those in Morus notabilis, Ficus hispida, Antiaris toxicaria, and Cannabis sativa. Integrating transcriptome data and Se content revealed that 12 BpNAC genes were associated with selenium biosynthesis. Subsequent RT-qPCR analysis corroborated the correlation between BpNAC59, BpNAC62 with sodium selenate, and BpNAC55 with sodium selenite. Subcellular localization experiments revealed the nuclear functions of BpNAC59 and BpNAC62. This study highlights the potential BpNAC transcription factors involved in selenium metabolism, providing a foundation for strategically breeding selenium-fortified paper mulberry.

Bioinformatics approach combined with experimental verification reveals OAS3 gene implicated in paclitaxel resistance in head and neck cancer.

This study aimed to identify a candidate gene associated with paclitaxel (PTX) resistance and to evaluate functionally its biological role in the PTX-resistant head and neck squamous cell carcinoma (HNSCC) cell lines and clinical specimens. Microarray data series containing samples of different types of cancers resistant to PTX were analyzed and then a candidate gene associated with PTX resistance was identified using various bioinformatics tools. After the suppression of the target gene expression, changes in cell viability and colony-forming ability were evaluated in PTX-resistant FaDu and SCC-9 cell lines. Bioinformatics analyses of upregulated genes in PTX-resistant cancer cells indicated that OAS3 was associated with PTX resistance. The downregulation of OAS3 expression significantly reduced the viability and colony-forming capacity of PTX-resistant SCC-9 cells by inducing apoptosis and cell cycle arrest at G0/G1 phase. The therapeutic targeting of OAS3 may resensitize PTX-resistant HNSCC cells with high OAS3 expression to PTX treatment.

Sequence and epigenetic variations of R2R3-MYB transcription factors determine the diversity of taproot skin and flesh colors in different cultivated types of radish (Raphanus sativus L.).

This study found that three paralogous R2R3-MYB transcription factors exhibit functional divergence among different subspecies and cultivated types in radish. Cultivated radish taproots exhibit a wide range of color variations due to unique anthocyanin accumulation patterns in various tissues. This study investigated the universal principles of taproot color regulation that developed during domestication of different subspecies and cultivated types. The key candidate genes RsMYB1 and RsMYB2, which control anthocyanin accumulation in radish taproots, were identified using bulked segregant analysis in two genetic populations. We introduced the RsMYB1-RsF3'H-RsMYB1Met genetic model to elucidate the complex and unstable genetic regulation of taproot flesh color in Xinlimei radish. Furthermore, we analyzed the expression patterns of three R2R3-MYB transcription factors in lines with different taproot colors and investigated the relationship between RsMYB haplotypes and anthocyanin accumulation in a natural population of 56 germplasms. The results revealed that three paralogous RsMYBs underwent functional divergence during radish domestication, with RsMYB1 regulating the red flesh of Xinlimei radish, and RsMYB2 and RsMYB3 regulating the red skin of East Asian big long radish (R. sativus var. hortensis) and European small radish (R. sativus var. sativus), respectively. Moreover, RsMYB1-H1, RsMYB2-H10, and RsMYB3-H6 were identified as the primary haplotypes exerting regulatory functions on anthocyanin synthesis. These findings provide an understanding of the genetic mechanisms regulating anthocyanin synthesis in radish and offer a potential strategy for early prediction of color variations in breeding programs.

Regulatory Mechanisms of Pollen Development: Transcriptomic and Bioinformatic Insights into the Role of β-1,3 Glucanase Gene (LbGlu1) in Lycium barbarum

Pollen fertility is a critical factor in seed development and crop breeding. Extensive studies have explored the mechanisms of pollen fertility in model plants and economic crops. However, the mechanisms of pollen abortion in medicinal and edible plants, including Lycium barbarum, remain elusive. This study utilized transcriptome analysis to identify key genes and regulatory networks implicated in pollen fertility in L. barbarum. The results demonstrated differential expression of 12,185 genes (DEGs) between the sterile and fertile lines, encompassing 489 genes that exhibited variation across the five stages of pollen development. Additionally, GO and KEGG enrichment analyses indicated that the DEGs were predominantly associated with energy metabolism, carbohydrate metabolism, and notably, hydrolase activity. Co-expression network analysis unveiled two modules intimately associated with fertility, each comprising 908 and 756 hub genes, incorporating β-1,3-glucanase genes (Glu) and co-expressed transcription factors (TFs). Phylogenetic analysis implied that LbGlu1 was a potential candidate gene implicated in regulating pollen abortion in L. barbarum. This work advances a novel understanding of pollen abortion in L. barbarum and offers theoretical support for the utilization of sterility genes to enhance crop improvement.

Understanding genome structure facilitates the use of wild lentil germplasm for breeding: A case study with shattering loci.

Plant breeders are generally reluctant to cross elite crop cultivars with their wild relatives to introgress novel desirable traits due to associated negative traits such as pod shattering. This results in a genetic bottleneck that could be reduced through better understanding of the genomic locations of the gene(s) controlling this trait. We integrated information on parental genomes, pod shattering data from multiple environments, and high-density genetic linkage maps to identify pod shattering quantitative trait loci (QTLs) in three lentil interspecific recombinant inbred line populations. The broad-sense heritability on a multi-environment basis varied from 0.46 (in LR-70, Lens culinaris × Lens odemensis) to 0.77 (in LR-68, Lens orientalis × L. culinaris). Genetic linkage maps of the interspecific populations revealed reciprocal translocations of chromosomal segments that differed among the populations, and which were associated with reduced recombination. LR-68 had a 2-5 translocation, LR-70 had 1-5, 2-6, and 2-7 translocations, and LR-86 had a 2-7 translocation in one parent relative to the other. Segregation distortion was also observed for clusters of single nucleotide polymorphisms on multiple chromosomes per population, further affecting introgression. Two major QTL, on chromosomes 4 and 7, were repeatedly detected in the three populations and contain several candidate genes. These findings will be of significant value for lentil breeders to strategically access novel superior alleles while minimizing the genetic impact of pod shattering from wild parents.

The lack of Tex44 causes severe subfertility with flagellar abnormalities in male mice

By analyzing a mouse Interspecific Recombinant Congenic Strain (IRCS), we previously identified a quantitative trait locus (QTL), called Mafq1 on mouse chromosome 1, that is associated with male hypofertility and ultrastructural sperm abnormalities. Within this locus, we identified a new candidate gene that could be implicated in a reproductive phenotype: Tex44 (Testis-expressed protein 44). We thus performed a CRISPR/Cas9-mediated complete deletion of this gene in mice in order to study its function. Tex44-KO males were severely hypofertile in vivo and in vitro due to a drastic reduction of sperm motility which itself resulted from important morphological sperm abnormalities. Namely, Tex44-KO sperm showed a disorganized junction between the midpiece and the principal piece of the flagellum, leading to a 180° flagellar bending in this region. In addition, the loss of some axonemal microtubule doublets and outer dense fibers in the flagellum’s principal piece has been observed. Our results suggest that, in mice, TEX44 is implicated in the correct set-up of the sperm flagellum during spermiogenesis and its absence leads to flagellar abnormalities and consequently to severe male hypofertility.Graphical