LettuceDB: an integrated multi-omics database for cultivated lettuce.

Crop genomics has advanced rapidly during the past decade, which generated a great abundance of omics data from multi-omics studies. How to utilize the accumulating data becomes a critical and urgent demand in crop science. As an attempt to integrate multi-omics data, we developed a database, LettuceDB (https://db.cngb.org/lettuce/), aiming to assemble multidimensional data for cultivated and wild lettuce germplasm. The database includes genome, variome, phenome, microbiome and spatial transcriptome. By integrating user-friendly bioinformatics tools, LettuceDB will serve as a one-stop platform for lettuce research and breeding in the future. Database URL: https://db.cngb.org/lettuce/.

The Molecular Genetics and Genomics of Plant-Pathogen Interactions.

Plants have evolved an intricate immune system to protect themselves from potential pathogens [...].

Metabolomic and genomic insights into TMA degradation by a novel halotolerant strain - Paracoccus sp. PS1.

Trimethylamine N-oxide (TMAO) is a gut metabolite that acts as a biomarker for chronic diseases, and is generated by the oxidation of trimethylamine (TMA) produced by gut microflora. Since, microbial degradation of TMA is predicted to be used to restrict the production of TMAO, we aimed to isolate bacterial strains that could effectively degrade TMA before being oxidized to TMAO. As marine fish is considered to have a rich content of TMAO, we have isolated TMA degrading isolates from fish skin. Out of the fourteen isolates, depending on their rapid TMA utilization capability in mineral salt medium supplemented with TMA as a sole carbon and nitrogen source, isolate PS1 was selected as our desired isolate. Its TMA degrading capacity was further confirmed through spectrophotometric, Electrospray Ionization Time-of-Flight Mass Spectrometry (ESI TOF-MS) and High performance liquid chromatography (HPLC) analysis and in silico analysis of whole genome (WG) gave further insights of protein into its TMA degradation pathways. PS1 was taxonomically identified as Paracoccus sp. based on its 16S rRNA and whole genome sequence analysis. As PS1 possesses the enzymes required for degradation of TMA, clinical use of this isolate has the potential to reduce TMAO generation in the human gut.

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

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

Aspiring toward equitable benefits from genomic advances to individuals of ancestrally diverse backgrounds.

Advancements in genomic technologies have shown remarkable promise for improving health trajectories. The Human Genome Project has catalyzed the integration of genomic tools into clinical practice, such as disease risk assessment, prenatal testing and reproductive genomics, cancer diagnostics and prognostication, and therapeutic decision making. Despite the promise of genomic technologies, their full potential remains untapped without including individuals of diverse ancestries and integrating social determinants of health (SDOHs). The NHGRI launched the 2020 Strategic Vision with ten bold predictions by 2030, including "individuals from ancestrally diverse backgrounds will benefit equitably from advances in human genomics." Meeting this goal requires a holistic approach that brings together genomic advancements with careful consideration to healthcare access as well as SDOHs to ensure that translation of genetics research is inclusive, affordable, and accessible and ultimately narrows rather than widens health disparities. With this prediction in mind, this review delves into the two paramount applications of genetic testing-reproductive genomics and precision oncology. When discussing these applications of genomic advancements, we evaluate current accessibility limitations, highlight challenges in achieving representativeness, and propose paths forward to realize the ultimate goal of their equitable applications.

Comparative genomics of a novel Erwinia species associated with the Highland midge (Culicoides impunctatus).

Erwinia (Enterobacterales: Erwiniaceae) are a group of cosmopolitan bacteria best known as the causative agents of various plant diseases. However, other species in this genus have been found to play important roles as insect endosymbionts supplementing the diet of their hosts. Here, I describe Candidatus Erwinia impunctatus (Erwimp) associated with the Highland midge Culicoides impunctatus (Diptera: Ceratopogonidae), an abundant biting pest in the Scottish Highlands. The genome of this new Erwinia species was assembled using hybrid long and short read techniques, and a comparative analysis was undertaken with other members of the genus to understand its potential ecological niche and impact. Genome composition analysis revealed that Erwimp is similar to other endophytic and ectophytic species in the genus and is unlikely to be restricted to its insect host. Evidence for an additional plant host includes the presence of a carotenoid synthesis operon implicated as a virulence factor in plant-associated members in the sister genus Pantoea. Unique features of Erwimp include several copies of intimin-like proteins which, along with signs of genome pseudogenization and a loss of certain metabolic pathways, suggests an element of host restriction seen elsewhere in the genus. Furthermore, a screening of individuals over two field seasons revealed the absence of the bacteria in Culicoides impunctatus during the second year indicating this microbe-insect interaction is likely to be transient. These data suggest that Culicoides impunctatus may have an important role to play beyond a biting nuisance, as an insect vector transmitting Erwimp alongside any conferred impacts to surrounding biota.

Genomic and transcriptomic studies on flavonoid biosynthesis in Lagerstroemia indica.

BACKGROUND: Lagerstroemia indica is a widely cultivated ornamental woody shrub/tree of the family Lythraceae that is used as a traditional medicinal plant in East Asia and Egypt. However, unlike other ornamental woody plants, its genome is not well-investigated, which hindered the discovery of the key genes that regulate important traits and the synthesis of bioactive compounds. RESULTS: In this study, the genomic sequences of L. indica were determined using several next-generation sequencing technologies. Altogether, 324.01 Mb sequences were assembled and 98.21% (318.21 Mb) of them were placed in 24 pseudo-chromosomes. The heterozygosity, repeated sequences, and GC residues occupied 1.65%, 29.17%, and 38.64% of the genome, respectively. In addition, 28,811 protein-coding gene models, 327 miRNAs, 552 tRNAs, 214 rRNAs, and 607 snRNAs were identified. The intra- and interspecies synteny and Ks analysis revealed that L. indica exhibits a hexaploidy. The co-expression profiles of the genes involved in the phenylpropanoid (PA) and flavonoid/anthocyanin (ABGs) pathways with the R2R3 MYB genes (137 members) showed that ten R2R3 MYB genes positively regulate flavonoid/anthocyanin biosynthesis. The colors of flowers with white, purple (PB), and deep purplish pink (DPB) petals were found to be determined by the levels of delphinidin-based (Dp) derivatives. However, the substrate specificities of LiDFR and LiOMT probably resulted in the different compositions of flavonoid/anthocyanin. In L. indica, two LiTTG1s (LiTTG1-1 and LiTTG1-2) were found to be the homologs of AtTTG1 (WD40). LiTTG1-1 was found to repress anthocyanin biosynthesis using the tobacco transient transfection assay. CONCLUSIONS: This study showed that the ancestor L. indica experienced genome triplication approximately 38.5 million years ago and that LiTTG1-1 represses anthocyanin biosynthesis. Furthermore, several genes such as LiDFR, LiOMTs, and R2R3 LiMYBs are related to anthocyanin biosynthesis. Further studies are required to clarify the mechanisms and alleles responsible for flower color development.

Comprehensive genomic analysis of Burkholderia arboris PN-1 reveals its biocontrol potential against Fusarium solani-induced root rot in Panax notoginseng.

Panax notoginseng (Burkill) F.H. Chen, a valuable traditional Chinese medicine, faces significant yield and quality challenges stemming from root rot primarily caused by Fusarium solani. Burkholderia arboris PN-1, isolated from the rhizosphere soil of P. notoginseng, demonstrated a remarkable ability to inhibit the growth of F. solani. This study integrates phenotypic, phylogenetic, and genomic analyses to enhance our understanding of the biocontrol mechanisms employed by B. arboris PN-1. Phenotype analysis reveals that B. arboris PN-1 effectively suppresses P. notoginseng root rot both in vitro and in vivo. The genome of B. arboris PN-1 comprises three circular chromosomes (contig 1: 3,651,544 bp, contig 2: 1,355,460 bp, and contig 3: 3,471,056 bp), with a 66.81% GC content, housing 7,550 protein-coding genes. Notably, no plasmids were detected. Phylogenetic analysis places PN-1 in close relation to B. arboris AU14372, B. arboris LMG24066, and B. arboris MEC_B345. Average nucleotide identity (ANI) values confirm the PN-1 classification as B. arboris. Comparative analysis with seven other B. arboris strains identified 4,628 core genes in B. arboris PN-1. The pan-genome of B. arboris appears open but may approach closure. Whole-genome sequencing revealed 265 carbohydrate-active enzymes and identified 9 gene clusters encoding secondary metabolites. This comprehensive investigation enhances our understanding of B. arboris genomes, paving the way for their potential as effective biocontrol agents against fungal plant pathogens in the future.

Exploiting integrative metabolomics to study host-parasite interactions in Plasmodium infections.

Despite years of research, malaria remains a significant global health burden, with poor diagnostic tests and increasing antimalarial drug resistance challenging diagnosis and treatment. While 'single-omics'-based approaches have been instrumental in gaining insight into the biology and pathogenicity of the Plasmodium parasite and its interaction with the human host, a more comprehensive understanding of malaria pathogenesis can be achieved through 'multi-omics' approaches. Integrative methods, which combine metabolomics, lipidomics, transcriptomics, and genomics datasets, offer a holistic systems biology approach to studying malaria. This review highlights recent advances, future directions, and challenges involved in using integrative metabolomics approaches to interrogate the interactions between Plasmodium and the human host, paving the way towards targeted antimalaria therapeutics and control intervention methods.

Genomic insights and anti-phytopathogenic potential of siderophore metabolome of endolithic Nocardia mangyaensis NH1.

Actinobacteria are one of the predominant groups that successfully colonize and survive in various aquatic, terrestrial and rhizhospheric ecosystems. Among actinobacteria, Nocardia is one of the most important agricultural and industrial bacteria. Screening and isolation of Nocardia related bacteria from extreme habitats such as endolithic environments are beneficial for practical applications in agricultural and environmental biotechnology. In this work, bioinformatics analysis revealed that a novel strain Nocardia mangyaensis NH1 has the capacity to produce structurally varied bioactive compounds, which encoded by non-ribosomal peptide synthases (NRPS), polyketide synthase (PKS), and post-translationally modified peptides (RiPPs). Among NRPS, five gene clusters have a sequence homology with clusters encoding for siderophore synthesis. We also show that N. mangyaensis NH1 accumulates both catechol- and hydroxamate-type siderophores simultaneously under iron-deficient conditions. Untargeted LC-MS/MS analysis revealed a variety of metabolites, including siderophores, lipopeptides, cyclic peptides, and indole-3-acetic acid (IAA) in the culture medium of N. mangyaensis NH1 grown under iron deficiency. We demonstrate that four CAS (chrome azurol S)-positive fractions display variable affinity to metals, with a high Fe3+ chelating capability. Additionally, three of these fractions exhibit antioxidant activity. A combination of iron scavenging metabolites produced by N. mangyaensis NH1 showed antifungal activity against several plant pathogenic fungi. We have shown that the pure culture of N. mangyaensis NH1 and its metabolites have no adverse impact on Arabidopsis seedlings. The ability of N. mangyaensis NH1 to produce siderophores with antifungal, metal-chelating, and antioxidant properties, when supplemented with phytohormones, has the potential to improve the release of macro- and micronutrients, increase soil fertility, promote plant growth and development, and enable the production of biofertilizers across diverse soil systems.

Assessing small RNA profiles in potato diploid hybrid and its resynthesized allopolyploid reveals conserved abundance with distinct genomic distribution.

KEY MESSAGE: The shock produced by the allopolyploidization process on a potato interspecific diploid hybrid displays a non-random remobilization of the small RNAs profile on a variety of genomic features. Allopolyploidy, a complex process involving interspecific hybridization and whole genome duplication, significantly impacts plant evolution, leading to the emergence of novel phenotypes. Polyploids often present phenotypic nuances that enhance adaptability, enabling them to compete better and occasionally to colonize new habitats. Whole-genome duplication represents a genomic "shock" that can trigger genetic and epigenetic changes that yield novel expression patterns. In this work, we investigate the polyploidization effect on a diploid interspecific hybrid obtained through the cross between the cultivated potato Solanum tuberosum and the wild potato Solanum kurtzianum, by assessing the small RNAs (sRNAs) profile of the parental diploid hybrid and its derived allopolyploid. Small RNAs are key components of the epigenetic mechanisms involved in silencing by RNA-directed DNA Methylation (RdDM). A sRNA sequencing (sRNA-Seq) analysis was performed to individually profile the 21 to 22 nucleotide (21 to 22-nt) and 24-nt sRNA size classes due to their unique mechanism of biogenesis and mode of function. The composition and distribution of different genomic features and differentially accumulated (DA) sRNAs were evaluated throughout the potato genome. We selected a subset of genes associated with DA sRNAs for messenger RNA (mRNA) expression analysis to assess potential impacts on the transcriptome. Interestingly, we noted that 24-nt DA sRNAs that exclusively mapped to exons were correlated with differentially expressed mRNAs between genotypes, while this behavior was not observed when 24-nt DA sRNAs were mapped to intronic regions. These findings collectively emphasize the nonstochastic nature of sRNA remobilization in response to the genomic shock induced by allopolyploidization.

Spatial multi-omics in medicinal plants: from biosynthesis pathways to industrial applications.

With the rapid development of molecular sequencing and imaging technology, the multi-omics of medicinal plants enters the single-cell era. We discuss spatial multi-omics applied in medicinal plants, evaluate the special products' biosynthesis pathways, and highlight the applications, perspectives, and challenges of biomanufacturing natural products (NPs).

Molecular mechanisms of selenite reduction by Lactiplantibacillus plantarum BSe: An integrated genomic and transcriptomic analysis.

The reduction of selenite [Se(Ⅳ)] by microorganisms is a green and efficient detoxification strategy. We found that Se(Ⅳ) inhibited exopolysaccharide and protein secretion by Lactiplantibacillus plantarum BSe and compromised cell integrity. In this study, L. plantarum BSe reduced Se(Ⅳ) by increasing related enzyme activity and electron transfer. Genomic analysis demonstrated that L. plantarum BSe should be able to reduce Se(Ⅳ). Further transcriptome analysis showed that L. plantarum BSe enhanced its tolerance to Se(Ⅳ) by upregulating the expression of surface proteins and transporters, thus reducing the extracellular Se(Ⅳ) concentration through related enzymatic reactions and siderophore-mediated pathways. Lactiplantibacillus plantarum BSe was able to regulate the expression of related genes involved in quorum sensing and a two-component system and then select appropriate strategies for Se(Ⅳ) transformation in response to varying environmental Se(Ⅳ) concentrations. In addition, azo reductase was linked to the reduction of Se(Ⅳ) for the first time. The present study established a multipath model for the reduction of Se(Ⅳ) by L. plantarum, providing new insights into the biological reduction of Se(Ⅳ) and the biogeochemical cycle of selenium.

Comprehensive Genomic Analysis Identifies a Diverse Landscape of Sideroblastic and Nonsideroblastic Iron-Related Anemias with Novel and Pathogenic Variants in an Iron-Deficient Endemic Setting.

Inherited iron metabolism defects are possibly missed or underdiagnosed in iron-deficient endemic settings because of a lack of awareness or a methodical screening approach. Hence, we systematically evaluated anemia cases (2019 to 2021) based on clinical phenotype, normal screening tests (high-performance liquid chromatography, α gene sequencing, erythrocyte sedimentation rate, C-reactive protein, and tissue transglutaminase), and abnormal iron profile by targeted next-generation sequencing (26-gene panel) supplemented with whole-exome sequencing, multiplex ligation probe amplification/mitochondrial DNA sequencing, and chromosomal microarray. Novel variants in ALAS2, STEAP3, and HSPA9 genes were functionally validated. A total of 290 anemia cases were screened, and 41 (14%) enrolled for genomic testing as per inclusion criteria. Comprehensive genomic testing revealed pathogenic variants in 23 of 41 cases (56%). Congenital sideroblastic anemia was the most common diagnosis (14/23; 61%), with pathogenic variations in ALAS2 (n = 6), SLC25A38 (n = 3), HSPA9 (n = 2) and HSCB, SLC19A2, and mitochondrial DNA deletion (n = 1 each). Nonsideroblastic iron defects included STEAP3-related microcytic anemia (2/23; 8.7%) and hypotransferrenemia (1/23; 4.3%). A total of 6 of 22 cases (27%) revealed a non-iron metabolism gene defect on whole-exome sequencing. Eleven novel variants (including variants of uncertain significance) were noted in 13 cases. Genotype-phenotype correlation revealed a significant association of frameshift/nonsense/splice variants with lower presentation age (0.8 months versus 9 years; P < 0.01) compared with missense variants. The systematic evaluation helped uncover an inherited iron defect in 41% (17/41) of cases, suggesting the need for active screening and awareness for these rare diseases in an iron-deficient endemic population.

Genomic insight into the origin, domestication, dispersal, diversification and human selection of Tartary buckwheat.

BACKGROUND: Tartary buckwheat, Fagopyrum tataricum, is a pseudocereal crop with worldwide distribution and high nutritional value. However, the origin and domestication history of this crop remain to be elucidated. RESULTS: Here, by analyzing the population genomics of 567 accessions collected worldwide and reviewing historical documents, we find that Tartary buckwheat originated in the Himalayan region and then spread southwest possibly along with the migration of the Yi people, a minority in Southwestern China that has a long history of planting Tartary buckwheat. Along with the expansion of the Mongol Empire, Tartary buckwheat dispersed to Europe and ultimately to the rest of the world. The different natural growth environments resulted in adaptation, especially significant differences in salt tolerance between northern and southern Chinese Tartary buckwheat populations. By scanning for selective sweeps and using a genome-wide association study, we identify genes responsible for Tartary buckwheat domestication and differentiation, which we then experimentally validate. Comparative genomics and QTL analysis further shed light on the genetic foundation of the easily dehulled trait in a particular variety that was artificially selected by the Wa people, a minority group in Southwestern China known for cultivating Tartary buckwheat specifically for steaming as a staple food to prevent lysine deficiency. CONCLUSIONS: This study provides both comprehensive insights into the origin and domestication of, and a foundation for molecular breeding for, Tartary buckwheat.

Genomic Analysis of Romanian Lycium Genotypes: Exploring BODYGUARD Genes for Stress Resistance Breeding.

Goji berries, long valued in Traditional Chinese Medicine and Asian cuisine for their wide range of medicinal benefits, are now considered a 'superfruit' and functional food worldwide. Because of growing demand, Europe and North America are increasing their goji berry production, using goji berry varieties that are not originally from these regions. European breeding programs are focusing on producing Lycium varieties adapted to local conditions and market demands. By 2023, seven varieties of goji berries were successfully registered in Romania, developed using germplasm that originated from sources outside the country. A broader project focused on goji berry breeding was initiated in 2014 at USAMV Bucharest. In the present research, five cultivated and three wild L. barbarum genotypes were compared to analyse genetic variation at the whole genome level. In addition, a case study presents the differences in the genomic coding sequences of BODYGUARD (BDG) 3 and 4 genes from chromosomes 4, 8, and 9, which are involved in cuticle-related resistance. All three BDG genes show distinctive differences between the cultivated and wild-type genotypes at the SNP level. In the BDG 4 gene located on chromosome 8, 69% of SNPs differentiate the wild from the cultivated genotypes, while in BDG 3 on chromosome 4, 64% of SNPs could tell the difference between the wild and cultivated goji berry. The research also uncovered significant SNP and InDel differences between cultivated and wild genotypes, in the entire genome, providing crucial insights for goji berry breeders to support the development of goji berry cultivation in Romania.

Phenotypic, Technological, Safety, and Genomic Profiles of Gamma-Aminobutyric Acid-Producing Lactococcus lactis and Streptococcus thermophilus Strains Isolated from Cow's Milk.

Gamma-aminobutyric acid (GABA)-producing lactic acid bacteria (LAB) can be used as starters in the development of GABA-enriched functional fermented foods. In this work, four GABA-producing strains each of Lactococcus lactis and Streptococcus thermophilus species were isolated from cow's milk, and their phenotypic, technological, and safety profiles determined. Genome analysis provided genetic support for the majority of the analyzed traits, namely, GABA production, growth in milk, and the absence of genes of concern. The operon harboring the glutamate decarboxylase gene (gadB) was chromosomally encoded in all strains and showed the same gene content and gene order as those reported, respectively, for L. lactis and S. thermophilus. In the latter species, the operon was flanked (as in most strains of this species) by complete or truncated copies of insertion sequences (IS), suggesting recent acquisition through horizontal gene transfer. The genomes of three L. lactis and two S. thermophilus strains showed a gene encoding a caseinolytic proteinase (PrtP in L. lactis and PrtS in S. thermophilus). Of these, all but one grew in milk, forming a coagulum of good appearance and an appealing acidic flavor and taste. They also produced GABA in milk supplemented with monosodium glutamate. Two L. lactis strains were identified as belonging to the biovar. diacetylactis, utilized citrate from milk, and produced significant amounts of acetoin. None of the strains showed any noticeable antibiotic resistance, nor did their genomes harbor transferable antibiotic resistance genes or genes involved in toxicity, virulence, or pathogenicity. Altogether these results suggest that all eight strains may be considered candidates for use as starters or components of mixed LAB cultures for the manufacture of GABA-enriched fermented dairy products.

Perceptions and beliefs of community gatekeepers about genomic risk information in African cleft research.

BACKGROUND: A fundamental ethical issue in African genomics research is how socio-cultural factors impact perspectives, acceptance, and utility of genomic information, especially in stigmatizing conditions like orofacial clefts (OFCs). Previous research has shown that gatekeepers (e.g., religious, political, family or community leaders) wield considerable influence on the decision-making capabilities of their members, including health issues. Thus, their perspectives can inform the design of engagement strategies and increase exposure to the benefits of genomics testing/research. This is especially important for Africans underrepresented in genomic research. Our study aims to investigate the perspectives of gatekeepers concerning genomic risk information (GRI) in the presence of OFCs in a sub-Saharan African cohort. METHODS: Twenty-five focus group discussions (FGDs) consisting of 214 gatekeepers (religious, community, ethnic leaders, and traditional birth attendants) in Lagos, Nigeria, explored the opinions of participants on genomic risk information (GRI), OFC experience, and the possibility of involvement in collaborative decision-making in Lagos, Nigeria. Transcripts generated from audio recordings were coded and analyzed in NVivo using thematic analysis. RESULTS: Three main themes-knowledge, beliefs, and willingness to act-emerged from exploring the perspective of gatekeepers about GRI in this group. We observed mixed opinions regarding the acceptance of GRI. Many participants believed their role is to guide and support members when they receive results; this is based on the level of trust their members have in them. However, participants felt they would need to be trained by medical experts to do this. Also, religious and cultural beliefs were crucial to determining participants' understanding of OFCs and the acceptance and utilization of GRI. CONCLUSIONS: Incorporating cultural sensitivity into public engagement could help develop appropriate strategies to manage conflicting ideologies surrounding genomic information in African communities. This will allow for more widespread access to the advances in genomics research in underrepresented populations. We also recommend a synergistic relationship between community health specialists/scientists, and community leaders, including spiritual providers to better understand and utilize GRI.

Genomic and co-expression network analyses reveal candidate genes for oil accumulation based on an introgression population in Upland cotton (Gossypium hirsutum).

KEY MESSAGE: Integrated QTL mapping and WGCNA condense the potential gene regulatory network involved in oil accumulation. A glycosyl hydrolases gene (GhHSD1) for oil biosynthesis was confirmed in Arabidopsis, which will provide useful knowledge to understand the functional mechanism of oil biosynthesis in cotton. Cotton is an economical source of edible oil for the food industry. The genetic mechanism that regulates oil biosynthesis in cottonseeds is essential for the genetic enhancement of oil content (OC). To explore the functional genomics of OC, this study utilized an interspecific backcross inbred line population to dissect the quantitative trait locus (QTL) interlinked with OC. In total, nine OC QTLs were identified, four of which were novel, and each QTL explained 3.62-34.73% of the phenotypic variation of OC. The comprehensive transcript profiling of developing cottonseeds revealed 3,646 core genes differentially expressed in both inbred parents. Functional enrichment analysis determined 43 genes were annotated with oil biosynthesis processes. Implementation of weighted gene co-expression network analysis showed that 803 differential genes had a significant correlation with the OC phenotype. Further integrated analysis identified seven important genes located in OC QTLs. Of which, the GhHSD1 gene located in stable QTL qOC-Dt3-1 exhibited the highest functional linkages with the other network genes. Phylogenetic analysis showed significant evolutionary differences in the HSD1 sequences between oilseed- and starch- crops. Furthermore, the overexpression of GhHSD1 in Arabidopsis yielded almost 6.78% higher seed oil. This study not only uncovers important genetic loci for oil accumulation in cottonseed, but also provides a set of new candidate genes that potentially influence the oil biosynthesis pathway in cottonseed.