Functional Metabolomics Describes the Yeast Biosynthetic Regulome.

Genome-metabolism interactions enable cell growth. To probe the extent of these interactions and delineate their functional contributions, we quantified the Saccharomyces amino acid metabolome and its response to systematic gene deletion. Over one-third of coding genes, in particular those important for chromatin dynamics, translation, and transport, contribute to biosynthetic metabolism. Specific amino acid signatures characterize genes of similar function. This enabled us to exploit functional metabolomics to connect metabolic regulators to their effectors, as exemplified by TORC1, whose inhibition in exponentially growing cells is shown to match an interruption in endomembrane transport. Providing orthogonal information compared to physical and genetic interaction networks, metabolomic signatures cluster more than half of the so far uncharacterized yeast genes and provide functional annotation for them. A major part of coding genes is therefore participating in gene-metabolism interactions that expose the metabolism regulatory network and enable access to an underexplored space in gene function.

An RNA-informed dosage sensitivity map reflects the intrinsic functional nature of genes.

Understanding dosage sensitivity or why Mendelian diseases have dominant vs. recessive modes of inheritance is crucial for uncovering the etiology of human disease. Previous knowledge of dosage sensitivity is mainly based on observations of rare loss-of-function mutations or copy number changes, which are underpowered due to ultra rareness of such variants. Thus, the functional underpinnings of dosage constraint remain elusive. In this study, we aim to systematically quantify dosage perturbations from cis-regulatory variants in the general population to yield a tissue-specific dosage constraint map of genes and further explore their underlying functional logic. We reveal an inherent divergence of dosage constraints in genes by functional categories with signaling genes (transcription factors, protein kinases, ion channels, and cellular machinery) being dosage sensitive, while effector genes (transporters, metabolic enzymes, cytokines, and receptors) are generally dosage resilient. Instead of being a metric of functional dispensability, we show that dosage constraint reflects underlying homeostatic constraints arising from negative feedback. Finally, we employ machine learning to integrate DNA and RNA metrics to generate a comprehensive, tissue-specific map of dosage sensitivity (MoDs) for autosomal genes.

A genome-scale deep learning model to predict gene expression changes of genetic perturbations from multiplex biological networks.

Systematic characterization of biological effects to genetic perturbation is essential to the application of molecular biology and biomedicine. However, the experimental exhaustion of genetic perturbations on the genome-wide scale is challenging. Here, we show TranscriptionNet, a deep learning model that integrates multiple biological networks to systematically predict transcriptional profiles to three types of genetic perturbations based on transcriptional profiles induced by genetic perturbations in the L1000 project: RNA interference, clustered regularly interspaced short palindromic repeat, and overexpression. TranscriptionNet performs better than existing approaches in predicting inducible gene expression changes for all three types of genetic perturbations. TranscriptionNet can predict transcriptional profiles for all genes in existing biological networks and increases perturbational gene expression changes for each type of genetic perturbation from a few thousand to 26 945 genes. TranscriptionNet demonstrates strong generalization ability when comparing predicted and true gene expression changes on different external tasks. Overall, TranscriptionNet can systemically predict transcriptional consequences induced by perturbing genes on a genome-wide scale and thus holds promise to systemically detect gene function and enhance drug development and target discovery.

Efficient gene editing of a model fern species through gametophyte-based transformation.

The clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated nuclease (Cas) system allows precise and easy editing of genes in many plant species. However, this system has not yet been applied to any fern species through gametophytes due to the complex characteristics of fern genomes, genetics, and physiology. Here, we established a protocol for gametophyte-based screening of single-guide RNAs (sgRNAs) with high efficiency for CRISPR/Cas9-mediated gene knockout in a model fern species, Ceratopteris richardii. We utilized the C. richardii ACTIN promoter to drive sgRNA expression and the enhanced CaMV 35S promoter to drive the expression of Streptococcus pyogenes Cas9 in this CRISPR-mediated editing system, which was employed to successfully edit a few genes, such as Nucleotidase/phosphatase 1 (CrSAL1) and Phytoene Desaturase (CrPDS), which resulted in an albino phenotype in C. richardii. Knockout of CrSAL1 resulted in significantly (P<0.05) reduced stomatal conductance (gs), leaf transpiration rate (E), guard cell length, and abscisic acid (ABA)-induced reactive oxygen species (ROS) accumulation in guard cells. Moreover, CrSAL1 overexpressing plants showed significantly increased net photosynthetic rate (A), gs, and E as well as most of the stomatal traits and ABA-induced ROS production in guard cells compared to the wild-type (WT) plants. Taken together, our optimized CRISPR/Cas9 system provides a useful tool for functional genomics in a model fern species, allowing the exploration of fern gene functions for evolutionary biology, herbal medicine discovery, and agricultural applications.

Variant graph craft (VGC): a comprehensive tool for analyzing genetic variation and identifying disease-causing variants.

BACKGROUND: The variant call format (VCF) file is a structured and comprehensive text file crucial for researchers and clinicians in interpreting and understanding genomic variation data. It contains essential information about variant positions in the genome, along with alleles, genotype calls, and quality scores. Analyzing and visualizing these files, however, poses significant challenges due to the need for diverse resources and robust features for in-depth exploration. RESULTS: To address these challenges, we introduce variant graph craft (VGC), a VCF file visualization and analysis tool. VGC offers a wide range of features for exploring genetic variations, including extraction of variant data, intuitive visualization, and graphical representation of samples with genotype information. VGC is designed primarily for the analysis of patient cohorts, but it can also be adapted for use with individual probands or families. It integrates seamlessly with external resources, providing insights into gene function and variant frequencies in sample data. VGC includes gene function and pathway information from Molecular Signatures Database (MSigDB) for GO terms, KEGG, Biocarta, Pathway Interaction Database, and Reactome. Additionally, it dynamically links to gnomAD for variant information and incorporates ClinVar data for pathogenic variant information. VGC supports the Human Genome Assembly Hg37 and Hg38, ensuring compatibility with a wide range of data sets, and accommodates various approaches to exploring genetic variation data. It can be tailored to specific user needs with optional phenotype input data. CONCLUSIONS: In summary, VGC provides a comprehensive set of features tailored to researchers working with genomic variation data. Its intuitive interface, rapid filtering capabilities, and the flexibility to perform queries using custom groups make it an effective tool in identifying variants potentially associated with diseases. VGC operates locally, ensuring data security and privacy by eliminating the need for cloud-based VCF uploads, making it a secure and user-friendly tool. It is freely available at https://github.com/alperuzun/VGC .

The tree peony nuclear factor Y transcription factor PrNF-YC2 promotes seed oil accumulation.

Seed oil not only provides energy for seed postgermination development but also provides essential nutrients and raw materials for human products. However, the transcriptional regulatory mechanism controlling seed oil accumulation remains largely unknown. Tree peony (Paeonia rockii) is an emerging woody oilseed crop in China that is known for its high-quality seed oil. Here, we revealed that a tree peony nuclear factor Y transcription factor, PrNF-YC2, is expressed predominantly in developing seeds and functions as an essential positive regulator of seed oil accumulation. PrNF-YC2 promoted oil accumulation in both transient ectopic overexpression Nicotiana benthamiana leaves and stable transgenic Arabidopsis thaliana seeds, globally upregulating the expression of genes involved in oil accumulation. In contrast, PrNF-YC2-silenced tree peony leaves using a virus-induced gene silencing system showed reduced oil content and expression of oil synthesis-related genes, including four master positive regulators contributing to oil accumulation, namely, LEAFY COTYLEDON1 (LEC1), ABSCISIC ACID INSENSITIVE3 (ABI3), FUSCA3 (FUS3), and WRINKLED1 (WRI1). We demonstrated that PrNF-YC2 directly activates PrLEC1 and PrABI3 alone and indirectly activates PrFUS3 and PrWRI1 by interacting with PrLEC1. Moreover, interaction with PrLEC1 also enhances the activation capacity of PrNF-YC2. The activation of these four master positive regulators by PrNF-YC2 triggered the upregulation of numerous oil synthesis-related genes, thus promoting oil accumulation. These findings provide new insights into the regulatory mechanism of seed oil accumulation and manipulation of PrNF-YC2 may be beneficial for enhancing oil yield in tree peony and other oilseed crops.

Genes divided according to the relative position of the longest intron show increased representation in different KEGG pathways.

Despite the fact that introns mean an energy and time burden for eukaryotic cells, they play an irreplaceable role in the diversification and regulation of protein production. As a common feature of eukaryotic genomes, it has been reported that in protein-coding genes, the longest intron is usually one of the first introns. The goal of our work was to find a possible difference in the biological function of genes that fulfill this common feature compared to genes that do not. Data on the lengths of all introns in genes were extracted from the genomes of six vertebrates (human, mouse, koala, chicken, zebrafish and fugu) and two other model organisms (nematode worm and arabidopsis). We showed that more than 40% of protein-coding genes have the relative position of the longest intron located in the second or third tertile of all introns. Genes divided according to the relative position of the longest intron were found to be significantly increased in different KEGG pathways. Genes with the longest intron in the first tertile predominate in a range of pathways for amino acid and lipid metabolism, various signaling, cell junctions or ABC transporters. Genes with the longest intron in the second or third tertile show increased representation in pathways associated with the formation and function of the spliceosome and ribosomes. In the two groups of genes defined in this way, we further demonstrated the difference in the length of the longest introns and the distribution of their absolute positions. We also pointed out other characteristics, namely the positive correlation between the length of the longest intron and the sum of the lengths of all other introns in the gene and the preservation of the exact same absolute and relative position of the longest intron between orthologous genes.

DhuFAP: a platform for gene functional analysis in Dendrobium huoshanense.

BACKGROUND: Dendrobium huoshanense, a traditional medicinal and food plant, has a rich history of use. Recently, its genome was decoded, offering valuable insights into gene function. However, there is no comprehensive gene functional analysis platform for D. huoshanense. RESULT: To address this, we created a platform for gene function analysis and comparison in D. huoshanense (DhuFAP). Using 69 RNA-seq samples, we constructed a gene co-expression network and annotated D. huoshanense genes by aligning sequences with public protein databases. Our platform contained tools like Blast, gene set enrichment analysis, heatmap analysis, sequence extraction, and JBrowse. Analysis revealed co-expression of transcription factors (C2H2, GRAS, NAC) with genes encoding key enzymes in alkaloid biosynthesis. We also showcased the reliability and applicability of our platform using Chalcone synthases (CHS). CONCLUSION: DhuFAP ( www.gzybioinformatics.cn/DhuFAP ) and its suite of tools represent an accessible and invaluable resource for researchers, enabling the exploration of functional information pertaining to D. huoshanense genes. This platform stands poised to facilitate significant biological discoveries in this domain.

Identification and functional characterization of bidirectional gene pairs and their intergenic regions in cotton.

BACKGROUND: In research to improve the quality of transgenic crops, it is often necessary to introduce multiple functionally related genes into recipient plants simultaneously to improve crop genetic traits effectively. Compared with unidirectional promoters, bidirectional promoters simultaneously regulate the expression of multiple genes and improve the efficiency of biotechnology. Therefore, in this study, bidirectional gene pairs were systematically analyzed in Gossypium hirsutum TM-1, and the structure, function and evolutionary relationships of the bidirectional genes were analyzed. The endogenous bidirectional promoters of cotton were mined, and their specific regulatory elements and biological functions were explored to provide useful promoter resources and a theoretical basis for cultivating new cotton germplasms with excellent fiber quality. RESULTS: Using an improved search model, a total of 1,383 bidirectional transcript pairs were identified in the Gossypium hirsutum TM-1 genome, and their gene structure and functional annotations were systematically analyzed. Thirty bidirectional intergenic sequences were randomly screened for promoter activity analysis via a transient expression system, and 25 intergenic sequences were found to have bidirectional promoter activity. Comparative analysis of the bidirectional gene profiles of the four cotton subspecies revealed that these subspecies presented abundant bidirectional gene pairs with high homology and that the bidirectional genes in the cotton subspecies were more similar in terms of their molecular functions, cellular components and biological processes. In addition, parallel analysis of bidirectional genes in dicotyledons and monocotyledons revealed that abundant bidirectional gene pairs exist in different species. Although the total number of orthologous bidirectional genes was similar, there was a significant difference in the number of orthologous bidirectional gene pairs between dicotyledons and monocotyledons. This evolutionary analysis of the function and structure of homologous bidirectional gene pairs in different varieties and different subspecies of the same species revealed potential pathways by which these gene pairs originated, which may be necessary for the evolution of a new species. CONCLUSION: In this study, many bidirectional gene pairs in Gossypium hirsutum TM-1 were identified using computer programming, and systematic analysis was conducted to explore their functions and evolutionary relationships. In addition, the promoter activity of the bidirectional intergenic sequences was verified. The combination of computer programming screening, experimental validation and other methods is expected to provide preferred bidirectional promoters for transgenic breeding work via multigene cotransformation methods, and this information is valuable for genetic engineering research and applications.

MGEGFP: a multi-view graph embedding method for gene function prediction based on adaptive estimation with GCN.

In recent years, a number of computational approaches have been proposed to effectively integrate multiple heterogeneous biological networks, and have shown impressive performance for inferring gene function. However, the previous methods do not fully represent the critical neighborhood relationship between genes during the feature learning process. Furthermore, it is difficult to accurately estimate the contributions of different views for multi-view integration. In this paper, we propose MGEGFP, a multi-view graph embedding method based on adaptive estimation with Graph Convolutional Network (GCN), to learn high-quality gene representations among multiple interaction networks for function prediction. First, we design a dual-channel GCN encoder to disentangle the view-specific information and the consensus pattern across diverse networks. By the aid of disentangled representations, we develop a multi-gate module to adaptively estimate the contributions of different views during each reconstruction process and make full use of the multiplexity advantages, where a diversity preservation constraint is designed to prevent the over-fitting problem. To validate the effectiveness of our model, we conduct experiments on networks from the STRING database for both yeast and human datasets, and compare the performance with seven state-of-the-art methods in five evaluation metrics. Moreover, the ablation study manifests the important contribution of the designed dual-channel encoder, multi-gate module and the diversity preservation constraint in MGEGFP. The experimental results confirm the superiority of our proposed method and suggest that MGEGFP can be a useful tool for gene function prediction.

Root hairs: an underexplored target for sustainable cereal crop production.

To meet the demands of a rising human population, plant breeders will need to develop improved crop varieties that maximize yield in the face of increasing pressure on crop production. Historically, the optimization of crop root architecture has represented a challenging breeding target due to the inaccessibility of the root systems. Root hairs, single cell projections from the root epidermis, are perhaps the most overlooked component of root architecture traits. Root hairs play a central role in facilitating water, nutrient uptake, and soil cohesion. Current root hair architectures may be suboptimal under future agricultural production regimes, coupled with an increasingly variable climate. Here, we review the genetic control of root hair development in the world's three most important crops-rice, maize, and wheat-and highlight conservation of gene function between monocots and the model dicot species Arabidopsis. Advances in genomic techniques including gene editing combined with traditional plant breeding methods have the potential to overcome many inherent issues associated with the design of improved root hair architectures. Ultimately, this will enable detailed characterization of the effects of contrasting root hair morphology strategies on crop yield and resilience, and the development of new varieties better adapted to deliver future food security.

Genes involved in the adhesion and invasion of Arcobacter butzleri.

Arcobacter butzleri is a foodborne pathogen that mainly causes enteritis in humans, but the number of cases of bacteraemia has increased in recent years. However, there is still limited knowledge on the pathogenic mechanisms of this bacterium. To investigate how A. butzleri causes disease, single knockout mutants were constructed in the cadF, ABU_RS00335, ciaB, and flaAB genes, which might be involved in adhesion and invasion properties. These mutants and the isogenic wild-type (WT) were then tested for their ability to adhere and invade human Caco-2 and HT29-MTX cells. The adhesion and invasion of A. butzleri RM4018 strain was also visualized by a Leica CTR 6500 confocal microscope. The adhesion and invasion abilities of mutants lacking the invasion antigen CiaB or a functional flagellum were lower than those of the WTs. However, the extent of the decrease varied depending on the strain and/or cell line. Mutants lacking the fibronectin (FN)-binding protein CadF consistently exhibited reduced abilities, while the inactivation of the other studied FN-binding protein, ABU_RS00335, led to a reduction in only one of the two strains tested. Therefore, the ciaB and flaAB genes appear to be important for A. butzleri adhesion and invasion properties, while cadF appears to be indispensable.

Identification of squalene epoxidase in triterpenes biosynthesis in Poria cocos by molecular docking and CRISPR-Cas9 gene editing.

BACKGROUND: Squalene epoxidase is one of the rate-limiting enzymes in the biosynthetic pathway of membrane sterols and triterpenoids. The enzyme catalyzes the formation of oxidized squalene, which is a common precursor of sterols and triterpenoids. RESULT: In this study, the squalene epoxidase gene (PcSE) was evaluated in Poria cocos. Molecular docking between PcSE and squalene was performed and the active amino acids were identified. The sgRNA were designed based on the active site residues. The effect on triterpene synthesis in P. cocos was consistent with the results from ultra-high-performance liquid chromatography-quadruplex time-of-flight-double mass spectrometry (UHPLC-QTOF-MS/MS) analysis. The results showed that deletion of PcSE inhibited triterpene synthesis. In vivo verification of PcSE function was performed using a PEG-mediated protoplast transformation approach. CONCLUSION: The findings from this study provide a foundation for further studies on heterologous biosynthesis of P. cocos secondary metabolites.

Transcriptome-wide identification of ARF gene family in medicinal plant Polygonatum kingianum and expression analysis of PkARF members in different tissues.

BACKGROUND: Polygonatum kingianum holds significant importance in Traditional Chinese Medicine due to its medicinal properties, characterized by its diverse chemical constituents including polysaccharides, terpenoids, flavonoids, phenols, and phenylpropanoids. The Auxin Response Factor (ARF) is a pivotal transcription factor known for its regulatory role in both primary and secondary metabolite synthesis. However, our understanding of the ARF gene family in P. kingianum remains limited. METHODS AND RESULTS: We employed RNA-Seq to sequence three distinct tissues (leaf, root, and stem) of P. kingianum. The analysis revealed a total of 31,558 differentially expressed genes (DEGs), with 43 species of transcription factors annotated among them. Analyses via gene ontology and the Kyoto Encyclopedia of Genes and Genomes demonstrated that these DEGs were predominantly enriched in metabolic pathways and secondary metabolite biosynthesis. The proposed temporal expression analysis categorized the DEGs into nine clusters, suggesting the same expression trends that may be coordinated in multiple biological processes across the three tissues. Additionally, we conducted screening and expression pattern analysis of the ARF gene family, identifying 12 significantly expressed PkARF genes in P. kingianum roots. This discovery lays the groundwork for investigations into the role of PkARF genes in root growth, development, and secondary metabolism regulation. CONCLUSION: The obtained data and insights serve as a focal point for further research studies, centred on genetic manipulation of growth and secondary metabolism in P. kingianum. Furthermore, these findings contribute to the understanding of functional genomics in P. kingianum, offering valuable genetic resources.

Identification of genomic binding sites and direct target genes for the transcription factor DDIT3/CHOP.

DDIT3 is a tightly regulated basic leucine zipper (bZIP) transcription factor and key regulator in cellular stress responses. It is involved in a variety of pathological conditions and may cause cell cycle block and apoptosis. It is also implicated in differentiation of some specialized cell types and as an oncogene in several types of cancer. DDIT3 was originally believed to act as a dominant-negative inhibitor by forming heterodimers with other bZIP transcription factors, preventing their DNA binding and transactivating functions. DDIT3 has, however, been reported to bind DNA and regulate target genes. Here, we employed ChIP sequencing combined with microarray-based expression analysis to identify direct binding motifs and target genes of DDIT3. The results reveal DDIT3 binding to motifs similar to other bZIP transcription factors, known to form heterodimers with DDIT3. Binding to a class III satellite DNA repeat sequence was also detected. DDIT3 acted as a DNA-binding transcription factor and bound mainly to the promotor region of regulated genes. ChIP sequencing analysis of histone H3K27 methylation and acetylation showed a strong overlap between H3K27-acetylated marks and DDIT3 binding. These results support a role for DDIT3 as a transcriptional regulator of H3K27ac-marked genes in transcriptionally active chromatin.

SNP-based heritability and selection analyses: Improved models and new results.

Complex-trait genetics has advanced dramatically through methods to estimate the heritability tagged by SNPs, both genome-wide and in genomic regions of interest such as those defined by functional annotations. The models underlying many of these analyses are inadequate, and consequently many SNP-heritability results published to date are inaccurate. Here, we review the modelling issues, both for analyses based on individual genotype data and association test statistics, highlighting the role of a low-dimensional model for the heritability of each SNP. We use state-of-art models to present updated results about how heritability is distributed with respect to functional annotations in the human genome, and how it varies with allele frequency, which can reflect purifying selection. Our results give finer detail to the picture that has emerged in recent years of complex trait heritability widely dispersed across the genome. Confounding due to population structure remains a problem that summary statistic analyses cannot reliably overcome. Also see the video abstract here: https://youtu.be/WC2u03V65MQ.

Functional analysis of AgJHAMT gene related to developmental period in Aphis gossypii Glover.

Aphis gossypii is one of the most economically important agricultural pests that cause serious crop losses worldwide, and the indiscriminate chemical application causes resistance development in A. gossypii, a major obstacle to successful control. In this study, we selected the up-regulated expression gene AgJHAMT, which was enriched into juvenile hormone pathway though transcriptome sequencing analysis of the cotton aphids that fed on transgenic cotton lines expressing dsAgCYP6CY3 (the TG cotton). The AgJHAMT gene was overexpressed in cotton aphids which fed on the TG cotton, and its expression profile during the nymphs was clarified. Then, silencing AgJHAMT could advance the developmental period of cotton aphids by 0.5 days compared with control groups. The T and t values of cotton aphids in the dsJHAMT treatment group (6.88 ± 0.15, 1.65 ± 0.06) were significantly shorter than that of the sprayed H2O control group (7.6 ± 0.14, 1.97 ± 0.09) (P < 0.05), respectively. The fast growth caused by AgJHAMT silencing was rescued by applying the JH analogue, methoprene. Overall, these findings clarified the function of AgJHAMT in the developmental period of A. gossypii. This study contributes to further clarify the molecular mechanisms of delaying the growth and development of cotton aphids by the transgenic cotton lines expressing dsAgCYP6CY3.

OsWNK9 regulates cadmium concentration in brown rice by restraining cadmium transport from straw to brown rice.

Selecting and breeding rice cultivars that enable strong cadmium (Cd) accumulation in rice straw but low accumulation in brown rice is a promising way to achieve Cd phytoremediation as well as to ensure the food safety of rice. Herein, we isolated a gene OsWNK9 from the quantitative trait locus associated with reducing Cd translocation from rice straw to brown rice and decreasing the Cd concentration in brown rice (BRCdC). Continuous strong expression of OsWNK9 was observed in nodes and internode and was induced after Cd supply. OsWNK9 was localized in the rice cell nucleus and participated in the regulation of Cd transport in yeast. Two independent oswnk9 rice mutants were generated via CRISPR/Cas9 gene-editing and showed significantly higher BRCdC than that of the wild type (WT). The BRCdC of knockout oswnk9 mutants was 0.227 mg kg-1and 0.238 mg kg-1, increased by 14 % and 19 % compared with that of the WT due to the lower Cd allocation in the basal stem, internode, and node III, which was unrelated to Cd uptake. Interestingly, OsWNK9 could promote iron (Fe) accumulation in rice under Cd-contaminated conditions, suggesting that OsWNK9 is an ideal gene for Cd phytoremediation and Fe biofortification in rice to support safe food production.

Functional studies of McSTE24, McCYP305a1, and McJHEH, three essential genes act in cantharidin biosynthesis in the blister beetle (Coleoptera: Meloidae).

Cantharidin is a toxic defensive substance secreted by most blister beetles when attacked. It has been used to treat many complex diseases since ancient times and has recently regained popularity as an anticancer agent. However, the detailed mechanism of the cantharidin biosynthesis has not been completely addressed. In this study, we cloned McSTE24 (encoding STE24 endopeptidase) from terpenoid backbone pathway, McCYP305a1 (encoding cytochrome P450, family 305) and McJHEH [encoding subfamily A, polypeptide 1 and juvenile hormone (JH) epoxide hydrolase] associated to JH synthesis/degradation in the blister beetle Mylabris cichorii (Linnaeus, 1758, Coleoptera: Meloidae). Expression pattern analyses across developmental stages in adult males revealed that the expressions of 3 transcripts were closely linked to cantharidin titer exclusively during the peak period of cantharidin synthesis (20-25 days old). In contrast, at other stages, these genes may primarily regulate different biological processes. When RNA interference with double-stranded RNA suppressed the expressions of the 3 genes individually, significant reductions in cantharidin production were observed in males and also in females following McJHEH knockdown, indicating that these 3 genes might primarily contribute to cantharidin biosynthesis in males, but not in females, while females could self-synthesis a small amount of cantharidin. These findings support the previously hypothesized sexual dimorphism in cantharidin biosynthesis during the adult phase. McCYP305a1 collaborates with its upstream gene McSTE24 in cantharidin biosynthesis, while McJHEH independently regulates cantharidin biosynthesis in males.

Integrating Multiple Database Resources to Elucidate the Gene Flow in Southeast Asian Pig Populations.

The domestic pig (Sus scrofa) and its subfamilies have experienced long-term and extensive gene flow, particularly in Southeast Asia. Here, we analyzed 236 pigs, focusing on Yunnan indigenous, European commercial, East Asian, and Southeast Asian breeds, using the Pig Genomics Reference Panel (PGRP v1) of Pig Genotype-Tissue Expression (PigGTEx) to investigate gene flow and associated complex traits by integrating multiple database resources. In this study, we discovered evidence of admixtures from European pigs into the genome of Yunnan indigenous pigs. Additionally, we hypothesized that a potential conceptual gene flow route that may have contributed to the genetic composition of the Diannan small-ear pig is a gene exchange from the Vietnamese pig. Based on the most stringent gene introgression scan using the fd statistic, we identified three specific loci on chromosome 8, ranging from 51.65 to 52.45 Mb, which exhibited strong signatures of selection and harbored the NAF1, NPY1R, and NPY5R genes. These genes are associated with complex traits, such as fat mass, immunity, and litter weight, in pigs, as supported by multiple bio-functionalization databases. We utilized multiple databases to explore the potential dynamics of genetic exchange in Southeast Asian pig populations and elucidated specific gene functionalities.